《UEL——ABAQUS非线性用户单元的开发》.pptVIP

* * ABAQUS非线性用户单元的开发 郭永进 庄 茁 清华大学工程力学系 北京 99年11月8日 ABAQUS 6.1版本发布会暨99中国地区用户会议 ABAQUS的用户子程序 ABAQUS 6.1版本发布会暨99中国地区用户会议 用户子程序 用户子程序为分析提供了一个非常强大、灵活的工具。 用户子程序： 可以提高几个ABAQUS选项的功能，仅用数据行的方式使用这些选项可能太过严格； 用FORTRAN语言写成； 在某些情况下可以利用ABAQUS的功能程序。 ABAQUS 6.1版本发布会暨99中国地区用户会议 ABAQUS中的用户子程序： CREEP－定义时间相关的粘塑性行为（蠕变和膨胀）的用户子程序 DFLOW－在压实分析中定义非均匀孔隙 流体速度的用户子程序 DFLUX－在热传导或质量扩散分析中定义非均匀分布流的用户子程序 DISP－定义边界条件的用户子程序 ABAQUS 6.1版本发布会暨99中国地区用户会议 DLOAD－定义非均匀分布载荷的用户子程序 FILM－在热传导分析中定义非均匀膜系数和相关的下沉温度的用户子程序 FLOW－在压实分析中定义非均匀渗流系数和相关的下沉孔隙压力的用户子程序 FRIC－定义接触面摩擦行为的用户子程序 GAPCON－在温度－位移耦合或者纯热传导分析中定义接触面或节点之间的热传导系数的用户子程序 ABAQUS 6.1版本发布会暨99中国地区用户会议 GAPELECTR－在热－电耦合分析中定义表面间导电系数的用户子程序 HARDINI－定义初始等效塑性应变和初始背应力张量的用户子程序 HETVAL－在热传导分析中定义内部热产生的用户子程序 MPC－定义多点约束的用户子程序 ORIENT－为运动耦合约束定义局部方向或者定义局部材料取向的用户子程序 ABAQUS 6.1版本发布会暨99中国地区用户会议 RSURFU－定义刚性表面的用户子程序 SDVINI －定义依赖于解的状态变量场初始值的用户子程序 SIGINI－定义初应力场的用户子程序 UCORR－定义随机响应载荷的交叉 －相关属性的用户子程序 UEL－定义单元的用户子程序 ABAQUS 6.1版本发布会暨99中国地区用户会议 非线性用户单元的开发 动机 定义用户单元 UEL接口 例：具有非线性截面行为的平面梁单元 例：新开发的两种用户单元 ABAQUS 6.1版本发布会暨99中国地区用户会议 概要 ABAQUS/Standard提供接口，允许用户实现线性和非线性有限单元。 非线性有限单元在用户子程序UEL中实现。 接口允许定义任意复杂的单元。 如果编写适当，用户单元可以在ABAQUS/Standard的大部分程序中应用。 可以在一个UEL子程序中实现多个用户单元，这些单元可以同时使用。 这里只讨论非线性用户单元的实现，并采用实例说明。 ABAQUS 6.1版本发布会暨99中国地区用户会议 动机 ABAQUS/Standard是一个功能强大的分析工具，拥有一个大单元库，有500多种单元－可以进行绝大部分复杂结构问题的分析。 然而，在有些情况下需要开发用户单元来补充ABAQUS的单元库： 所需的单元几何信息与ABAQUS已有的单元不同 模拟与结构行为耦合的非结构物理过程 施加依赖于解的载荷 模拟主动控制机理 ABAQUS 6.1版本发布会暨99中国地区用户会议 在分析程序如ABAQUS中开发用户单元而不是编写一个完全的分析程序的优点是很明显的： ABAQUS提供大量的结构单元、分析程序和建模工具。 ABAQUS提供前处理和后处理。 许多第三方软件提供带有ABAQUS接口的前后处理器。 维护和移植子程序比维护和移植一个完全的有限元程序容易的多。 ABAQUS 6.1版本发布会暨99中国地区用户会议 定义用户单元 用户单元的主要特征 在编写UEL子程序前，必须先定义单元的以下主要特征： 单元的节点数 节点的坐标数 每一个节点处的自由度数 ABAQUS 6.1版本发布会暨99中国地区用户会议 其它重要的单元性质 另外，必须定义下面的单元性质： 需要在UEL外定义的单元性质数目 每一个单元需要存储的依赖于解的状态变量数目（SDVs） 可能会作用于单元上的（分布）载荷类型数目 这些数目不必立即定义：它们可以很容易地在完成基本的UEL子程序后加入。 ABAQUS 6.1版本发布会暨99中国地区用户会议 定义用户单元行为 在一般的分析步中单元对模型的主要贡献是提供作用于节点上的“力”FN－依赖于节点的自由度uN的值。 如果自由度是位移，那么相关的“力”为节点力。类似地，旋转对应于转矩，温度对应于热流。 在非线性用户单元中“力”常常依赖于自由度的增量?uN和内部状态变量H?。 状态变量必须在用户子程序中进行更新。