abaqus2用户单元子程序

《abaqus2用户单元子程序》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、20ABAQUS用户单元子程序（UEL）在这一章屮将列举两个在这些年里发展过的ABAQUS/Standard用户单元子程序（UEL）。第一个例子是一个非线性的索单元，我们的目的是通过这个比较简单的例子让读者了解用户单元子程序的基本开发过程；第二个例子是一个用于计算应变梯度理论的单元，应变梯度是当今比较热点的一个科研前沿问题，有各种理论，我们为了验证新的理论，需要数值结果与实验对照来进行评价，整个例子的目的是通过它说明用户子单元可以求解的问题范围很广，但是由于内容比较艰深，程序也很长，所以这个例子我们并没有给出最后的全部程序。另外，到目前为止，ABAQUS还只有隐式求解器ABAQUS/Stan

2、dard支持用户自定义单元，而显式求解器ABAQUS/Explicit中还不支持这一功能。20.1非线性索单元20.1.1背景钢索斜拉桥和斜拉索结构广泛应用于土木工程建筑上。索力的计算分析是设计和施工的关键环节。清华大学工程力学系在采用ABAQUS进行荆沙长江斜拉桥的计算机仿真分析（这个项目我们已在第15章“ABAQUS在土木工程屮的应用（一）一一荆州长江大桥南汉斜拉桥结构三维仿真分析”屮讨论过）时，也曾进行了自行建立索单元的尝试。本节介绍的就是这方面的工作。香港理工大学土木与结构工程系采用ABAQUS有限元软件进行计算，完成了香港TingKau斜拉桥和TsingMa悬索桥的结构计算和分析。

3、对于钢索计算，他们采用梁单元进行模拟。由于梁单元含有弯曲刚度，计算的高阶频率值偏高，周期较低。一般假设索是单向受拉力的构件。随着应变的非线性增加，索力呈非线性增加。尽管ABAQUS单元库屮有500个以上的单元类型，但是，还没有索单元。本文发展了三维非线性索单元模型，形成ABAQUS的用户单元子程序，可以利用ABAQUS输入文件调入到具体的分析屮。通过静态和动态例题的计算比较，索单元工作良好。20.1.2基本公式在三维索单元计算屮，如图20・1所示，坐标x和位移”的变量表达式为:Xji=Xj_X,(xfytz)(u,vfw)Uy=Uj应变的公式为：11/222£=7XjMji+y心+ZjiWj

4、i+-W"++wjf公式(20・2)中，厶为索的长度，索的张力为：N=eAE+N()(20-1)(20-2)(20-3)在总体坐标系下，单元刚度矩阵为:[K]加K-K-KK单元刚度矩阵屮的子阵K分别由线性和非线性矩阵项组成:(20-4)xjix訂”xjizjiy；y^ji25二TLzmj0_01(20-6)(20-7)(20-8)求解运动方程,在公式（20-5）中的位和K也均是3X3的对称矩阵，分别为:索单元的节点质量为:m=—pAL在公式（20-6）屮，。为密度。索单元的质塑矩阵为:[耐结构的运动方程为:公式⑵・8）中严'为作用在结构上的外力。在不断变化的索的变形中,得到节点的位移值。20

5、.1.3应用举例图19-2由五个单元组成的两端饺接的索杆结构由5个单兀组成的两端钱接的索朴结构，咼5m长10m,6个节点号码依次为101〜106,如图19-2所示。计算自由振动的频率和周期。输入文件中的用户单元界面ABAQUS输入文件(.inp)中的用户单元界面如下：★HEADINGTwodimensionaloverheadhoistframeusing2nodesself-developedtrusselement.InitialforceNisdefinedinproperty(5)andreferencedbyuserelementSIUnits1-axishorizontal,2-

6、axisvertical*USERELEMENT,NODES=2,TYPE=U1,PROPERTIES=5,COORDINATES=3,VARIABLES=121,2,3*UELPROPERTY,ELSET=UTRUSS1.963E-5,2.0E11,0.3,7800,10.0E5★ELEMENT,TYPE=U1.ELSET=UTRUSS计算结果和比较表20-1列出了由用户索单元计算的图20-2所示结构的固有周期，并与应用ABAQUS梁单元B31的计算结杲进行了比较。索单元与梁单元前4阶模态的周期基本一致；索单元的第6〜9阶模态与梁单元第7〜10阶模态的周期基本一致。从第11阶模态开始，随着

7、梁单元弯曲变形的增加，梁的弯曲刚度逐渐发挥作用并和轴向刚度耦合，与同阶模态的索单元相比，梁单元的振动周期显著降低，而频率高于索单元。表20-1ABAQUS用户索单元和梁单元B31计算的频率比较振动模态索单元固有周期(Cycles/time)梁单元固有周期(Cycles/time)1112.42112.482112.42112.483213.83213.954213.83213.955249.51222.556