ABAQUS中的DLOAD子程序应用

leejack

　　(1) ABAQUS中的子程序介绍
　　ABAQUS中随然提供了很多的单元，很多的材料本构模型，但是有时候我们的分析这些单元无法达到我们的要求。比如我们要分析一个编织型的复合材料的时候、我们要考虑层合板复合材料的层间力的时候、我们要计算一个特殊材料橡胶的大变形的时候，我们需要建立自己的本构。在我们加载的时候，当这个力的幅值与坐标有关的时候、当这个力的表达式是一个非线性的函数的时候，我们就需要使用DLOAD的子程序。
　　(2) DLOAD子程序参数介绍

　　其中kstep为分析步参数、kinc为增量步参数、time为主程序传入进来的分析时间参数noel为单元编号、npt为单元积分点编号。coords为积分点的坐标。Jltyp为载荷识别号。f为所加载的载荷。其中jltyp这里我们说明一下。ABAQUS中的单元面的编号一般为逆时针。如下图所示：

　　(3) DLOAD应用实例

　　我们要编写的程序用在下面模型上：

　　在分析中我们使用4节点双线性轴对称单元，通过命令行方式提交分析，如下所示：

　　call abq6141 job=ss user=ss int ss为INP文件名与写的子程序文件名。子程序的主要内容如下图所示：

　　在abaqus中观察分析结果，如下图所示：


　　由上两图可知，S11径向最大应力28620pa。S22周向最大应力7885pa。理论计算得S11精确值30000pa。S22精确值7890pa。由此可证明我们程序计算的结果是正确的。