abaqus自定义载荷子程序------Dload使用
abaqus子程序Dload的主要作用：
（1）可用于定义作为位置、时间、单元编号、被加载积分点数量等的函数分布载荷大小的变化。
（2）在应力分析期间，将在每个基于单元或基于表面的非均匀分布载荷定义的载荷积分点处调用；
（3）将在每个积分点调用，以计算承受不均匀荷载类型PENU和PINU的管道元件的有效轴向力ESF1；
（4）不能在基于模态的程序中用于描述负载的时间变化；并且忽略可能与相关联的阶跃定义或非均匀分布负载定义一起出现的任何幅度参考。
子程序接口界面
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 F RETURN END
待定义变量
F：分布载荷的大小。表面载荷的单位为FL−2，体力的单位为FL−3。F将作为基于单元或基于表面的分布式载荷定义的一部分指定的载荷大小传递到例程中。如果未定义大小，F将作为零传入。对于使用修正Riks法（静态应力分析）的静态分析，F必须定义为荷载比例系数λ的函数。分布式负载大小不可用于输出目的。
用于传递信息的变量
KSTEP：Step 编号
KINC：增量数
TIME(1)：当前分析步对应的当前时间
TIME(2)：所有分析步对应的当前时间
NOEL：单元编号
NPT：根据荷载类型，构件内或构件表面上的荷载积分点编号。
LAYER：层编号
KSPT：当前层截面点编号
COORDS：包含被加载积分点坐标的数组。如果在步骤中考虑几何非线性，则这些是当前坐标（参见定义分析）；否则，数组包含点的原始坐标。对于允许非轴对称变形的轴对称单元，COORDS（3）是积分点的角位置，单位为度
JLTYP：标识要调用DLOAD的加载类型。载荷类型可以是体力、基于表面的载荷或基于单元的表面载荷。对于基于单元的表面加载，此变量标识要调用DLOAD的元素面。当多个不同的非均匀分布载荷同时施加在一个元件上时，该信息非常有用。有关图元面标识，请参见关于图元库。关键如下：
SNAME:基于曲面的载荷定义的曲面名称（JLTYP=0）。对于体力或基于元素的表面载荷，表面名称将作为空白传递。
案例演示
三维六面体中心圆形面载荷加载模拟
几何模型20*20*1mm
材料采用纯弹性的钢参数，杨氏模量210Gpa，泊松比0.33，使用线性静力分析，平面四周被完全固定，中心半径为2.0mm区域施加100Mpa的表面力，并应用Dload子程序进行实现
结果与后处理
应力分布情况：
位移分布情况：
类似的可以生成其他形式的载荷，如随着时间移动的载荷，如铁轨上的火车的移动等
