DLOAD子程序说明
————————————————————————————————————
Use this subroutine to define nonuniform, distributed
mechanical loads (pressures and body forces).
使用此子程序定义非均匀的分布式机械负载（压力和体力）。
（
当负载是时间和/或位置的复杂函数时，通常使用用户子程序DLOAD。
通常可以用* AMPLITUDE选项建模简单的时间函数。
子程序也可以用来定义一个随元素号和/或积分点数而变化的负载。
）
————————————————————————————————————
一般用法，在INP文件中定义载荷施加的位置（关键字），求解时让其与子DLOAD子程序关联。
（
当* DLOAD或* DSLOAD选项包含不均匀的加载类型标签时，子程序被调用。
）
*DLOAD 1, P4NU
（Apply nonuniform DLOAD to face 4 of element 1，）
（
“ABAQUS /Standard”手册中给出了可用于任何特定单元的非均匀分布载荷类型列表。
）
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）。 对于身体力或基于单元的表面载荷，表面名称以空白形式传递
供参考的子程序
————————————————————————————————————
SUBROUTINE DLOAD(F, KSTEP, KINC, TIME ,NOEL, NPT, 1 LAYER, KSPT, COORDS, JLTYP, SNAME) C C EXPONENTIAL PRESSURE LOAD C INCLUDE ‘ABA_PARAM.INC’ C DIMENSION COORDS(3),TIME(2) CHARACTER*80 SNAME DATA TEN,ONE,CONST /10.,1.,-23.03/ F=TEN*(ONE-(EXP(CONST*TIME(1)))) IF(NPT.EQ.1) WRITE(6,*) ‘ LOAD APPLIED’,F,’AT TIME=’,TIME(1) RETURN END
————————————————————————————————————
（
用户只能定义变量F 这是分布式负载的大小，单位取决于分布式负载的类型：
FL–1 对于沿一维（梁）单元的线载荷，
FL–2 表面负荷（例如压力）和
FL–3 体力（例如重力，向心力，加速度）。
）