基于Abaqus的DLOAD子程序实现移动载荷

有限元模型中，当载荷比较复杂时，难以通过ABAQUS/CAE界面直接进行设置，这时候就需要使用DLOAD子程序。在工程实际应用中，经常会遇到移动载荷的例子，如车辙实验，汽车过桥等。本案例介绍基于ABAQUS子程序的DLOAD实现移动载荷。

DLOAD子程序接口：

         SUBROUTINE DLOAD(F,KSTEP,KINC,TIME,NOEL,NPT,LAYER,KSPT,

     1 COORDS,JLTYP,SNAME)

C

      INCLUDE 'ABA_PARAM.INC'

C

      DIMENSION TIME(2), COORDS (3)

      CHARACTER*80 SNAME

      user coding to define F 

      RETURN

      END

子程序中的各个参数的意义可以参考ABAQUS帮助文档

F——分布负载的大小

KSTEP——分析步

TIME(1)——步进时间的当前值或负载比例因子的当前值

TIME(2)——当前总时间值

NOEL——单元编号

NPT——根据负载类型，在元素或元素表面加载积分点数。

LAYER——层数

KSPT——当前图层中的Section point

COORDS——积分点坐标（数组）

JLTYP——DLOAD调用的载荷类型的标识

 SNAME——表面加载定义的表面名称（JLTYP = 0）。 对于体力或基于单元的表面载荷，表面名称以空白形式传递

 下面以一个简单例子为例进行说明，最终的效果如下图所示。

有限元模型主要设置：

1）建立“工字形”梁几何模型（3D可变形体），单位为mm；

2）建立静态-通用分析步，为使得载荷等间距移动，采用固定增量步；

3）“工字形”梁一端建立固定约束，上表面建立Pressure载荷，Distribution采用User-defined；

4）网格全局尺寸为0.5mm。