leejack 发表于 2016-11-14 14:13

ABAQUS中的DLOAD子程序应用

  (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。由此可证明我们程序计算的结果是正确的。


页: [1]
查看完整版本: ABAQUS中的DLOAD子程序应用