ABAQUS DLOAD子程序
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 体力(例如重力,向心力,加速度)。
)