abaqus2用户单元子程序.docx

20 ABAQUS用户单元子程序（UEL） 在这一章屮将列举两个在这些年里发展过的ABAQUS/Standard用户单元子程序 （UEL）。第一个例子是一个非线性的索单元，我们的目的是通过这个比较简单的例子 让读者了解用户单元子程序的基本开发过程；第二个例子是一个用于计算应变梯度理论 的单元，应变梯度是当今比较热点的一个科研前沿问题，有各种理论，我们为了验证新 的理论，需要数值结果与实验对照来进行评价，整个例子的目的是通过它说明用户子单 元可以求解的问题范围很广，但是由于内容比较艰深，程序也很长，所以这个例子我们 并没有给出最后的全部程序。 另外，到目前为止，ABAQUS还只有隐式求解器ABAQUS/Standard支持用户自定 义单元，而显式求解器ABAQUS/Explicit中还不支持这一功能。 20.1非线性索单元 20.1.1背景 钢索斜拉桥和斜拉索结构广泛应用于土木工程建筑上。索力的计算分析是设计和施 工的关键环节。清华大学工程力学系在采用ABAQUS进行荆沙长江斜拉桥的计算机仿 真分析（这个项目我们已在第15章“ABAQUS在土木工程屮的应用（一）一一荆州长 江大桥南汉斜拉桥结构三维仿真分析”屮讨论过）时，也曾进行了自行建立索单元的尝 试。本节介绍的就是这方面的工作。 香港理工大学土木与结构工程系采用ABAQUS有限元软件进行计算，完成了香港 Ting Kau斜拉桥和Tsing Ma悬索桥的结构计算和分析。对于钢索计算，他们采用梁单 元进行模拟。由于梁单元含有弯曲刚度，计算的高阶频率值偏高，周期较低。 一般假设索是单向受拉力的构件。随着应变的非线性增加，索力呈非线性增加。尽 管ABAQUS单元库屮有500个以上的单元类型，但是，还没有索单元。本文发展了三 维非线性索单元模型，形成ABAQUS的用户单元子程序，可以利用ABAQUS输入文件 调入到具体的分析屮。通过静态和动态例题的计算比较，索单元工作良好。 20.1.2基本公式 在三维索单元计算屮，如图20?1所示，坐标x和位移”的变量表达式为: Xji=Xj_X, (xfytz) (u,vfw) Uy = Uj 应变的公式为： 1 1 / 2 2 2 ￡ = 7 XjMji + y心 + ZjiWji + - W + + wjf 公式(20?2)中，厶为索的长度，索的张力为： N = eAE + N() (20-1) (20-2) (20-3) 在总体坐标系下，单元刚度矩阵为: [K]加 K -K -K K 单元刚度矩阵屮的子阵K分别由线性和非线性矩阵项组成: (20-4) xji x xji x 訂 ” xjizji y； y^ji 2 5 二 T L zm j 0_ 0 1 (20-6) (20-7) (20-8) 求解运动方程, 在公式（20-5）中的位和K也均是3X3的对称矩阵，分别为: 索单元的节点质量为: m = — pAL 在公式（20-6）屮，。为密度。索单元的质塑矩阵为: [耐 结构的运动方程为: 公式 ⑵?8）中严为作用在结构上的外力。在不断变化的索的变形中, 得到节点的位移值。 20.1.3应用举例 图19-2由五个单元组成的两端饺接的索杆结构 由5个单兀组成的两端钱接的索朴结构，咼5m长10m, 6个节点号码依次为 101?106,如图19-2所示。计算自由振动的频率和周期。 输入文件中的用户单元界面 ABAQUS输入文件(.inp)中的用户单元界面如下： ★HEADING Two dimensional overhead hoist frame using 2 nodes self-developed truss element. Initial force N is defined in property(5) and referenced by user element SI Units 1-axis horizontal, 2-axis vertical *USER ELEMENT, NODES=2, TYPE=U1, PROPERTIES=5, COORDINATES=3, VARIABLES=12 1,2,3 *UEL PROPERTY, ELSET=UTRUSS 1.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阶模态开始，随着梁单元弯曲变形的增加，梁的 弯曲刚度逐渐发挥作用并和轴向刚度耦合，与同