ABAQUS非线性用户单元的开发例析.ppt

ABAQUS非线性用户单元的开发;ABAQUS的用户子程序;用户子程序;ABAQUS中的用户子程序： CREEP－定义时间相关的粘塑性行为（蠕变和膨胀）的用户子程序 DFLOW－在压实分析中定义非均匀孔隙 流体速度的用户子程序 DFLUX－在热传导或质量扩散分析中定义非均匀分布流的用户子程序 DISP－定义边界条件的用户子程序;DLOAD－定义非均匀分布载荷的用户子程序 FILM－在热传导分析中定义非均匀膜系数和相关的下沉温度的用户子程序 FLOW－在压实分析中定义非均匀渗流系数和相关的下沉孔隙压力的用户子程序 FRIC－定义接触面摩擦行为的用户子程序 GAPCON－在温度－位移耦合或者纯热传导分析中定义接触面或节点之间的热传导系数的用户子程序;GAPELECTR－在热－电耦合分析中定义表面间导电系数的用户子程序 HARDINI－定义初始等效塑性应变和初始背应力张量的用户子程序 HETVAL－在热传导分析???定义内部热产生的用户子程序 MPC－定义多点约束的用户子程序 ORIENT－为运动耦合约束定义局部方向或者定义局部材料取向的用户子程序;RSURFU－定义刚性表面的用户子程序 SDVINI －定义依赖于解的状态变量场初始值的用户子程序 SIGINI－定义初应力场的用户子程序 UCORR－定义随机响应载荷的交叉 －相关属性的用户子程序 UEL－定义单元的用户子程序;非线性用户单元的开发;概要;动机;在分析程序如ABAQUS中开发用户单元而不是编写一个完全的分析程序的优点是很明显的： ABAQUS提供大量的结构单元、分析程序和建模工具。 ABAQUS提供前处理和后处理。 许多第三方软件提供带有ABAQUS接口的前后处理器。 维护和移植子程序比维护和移植一个完全的有限元程序容易的多。 ;定义用户单元;其它重要的单元性质 另外，必须定义下面的单元性质： 需要在UEL外定义的单元性质数目 每一个单元需要存储的依赖于解的状态变量数目（SDVs） 可能会作用于单元上的（分布）载荷类型数目 这些数目不必立即定义：它们可以很容易地在完成基本的UEL子程序后加入。 ;定义用户单元行为 在一般的分析步中单元对模型的主要贡献是提供作用于节点上的“力”FN－依赖于节点的自由度uN的值。 如果自由度是位移，那么相关的“力”为节点力。类似地，旋转对应于转矩，温度对应于热流。 在非线性用户单元中“力”常常依赖于自由度的增量?uN和内部状态变量H?。 状态变量必须在用户子程序中进行更新。;在一般分析步中求解（非线性）方程组需要你定义单元的刚度矩阵： 刚度矩阵应包含 FN对uN所有的直接和间接的依赖性，例如 在一般分析步中刚度矩阵定义的越精确，收敛性就越好。 刚度矩阵确定线性扰动分析步的解，所以其定义必须精确。 刚度矩阵可以为对称或非对称。;用户单元的复杂程度可以变化很大。 在由常规单元构成的分析中可以采用简单的用户单元来达到“控制”和“反馈”的功能。 复杂的非线性结构单元的编制常常需要很大的工作量。 如果单元由非线性材料构成，你最好用单独的子程序描述材料行为。 如果材料模型是在子程序UMAT中实现的，那么可以在UEL中包含对UMAT的调用。 前面对UMAT中积分的讨论同样适用于UEL中的材料模型。;UEL接口;UEL接口 用户子程序UEL的接口如下： SUBROUTINE UEL(RHS,AMATRX,SVARS,ENERGY,NDOFEL,NRHS,NSVARS, 1 PROPS,NPROPS,COORDS,MCRD,NNODE,U,DU,V,A,JTYPE,TIME,DTIME, 2 KSTEP,KINC,JELEM,PARAMS,NDLOAD,JDLTYP,ADLMAG,PREDEF, 3 NPREDF,LFLAGS,MLVARX,DDLMAG,MDLOAD,PNEWDT,JPROPS,NJPROP, 4 PERIOD) C INCLUDE ABA_PARAM.INC C DIMENSION RHS(MLVARX,*),AMATRX(NDOFEL,NDOFEL), 1 SVARS(NSVARS),ENERGY(8),PROPS(*),COORDS(MCRD,NNODE), 2 U(NDOFEL),DU(MLVARX,*),V(NDOFEL),A(NDOFEL),TIME(2), 3 PARAMS(3),JDLTYP(MDLOAD,*),ADLMAG(MDLOAD,*), 4 DDLMAG(MDLOAD,*),PREDEF(2,NPREDF,NNODE),LFLAGS