18.1.1 User subroutines: overview。
PNEWDT
PNEWDT指的是建议新的时间增量与所使用时间增量间的比例(DTIME,如下)。此变量能够允许你在Abaqus/Standard中对自动时间增量提供输入(假如选择了自动时间增量的话)。 在每次调用USDFLD之前,将PNEWDT设定为一个大数值
如果PNEWDT被重新定义为小于1.0的数,Abaqus/Standard必须舍弃时间增量,用一个更小时间增量再次尝试。提供给自动时间积分算法的建议新的时间增量是PNEWDT × DTIME,其中PNEWDT是允许迭代中重新定义PNEWDT的所有调用用户子程序中最小值。
如果PNEWDT是调用用户子程序中一个给定大于1.0的值,Abaqus/Standard可能会增加时间增量。提供给自动时间积分算法的建议新的时间增量是PNEWDT × DTIME,其中PNEWDT是允许迭代中重新定义PNEWDT的所有调用用户子程序中最小值。
如果在分析程序中没有选择自动时间增量,将忽略大于1.0的PNEWDT,小于1.0的PNEWDT会造成任务停止。
传递信息的变量 DIRECT (3,3)
包含物质点关于全局方向的方向余弦的数组。DIRECT(1,1) , DIRECT(2,1) , DIRECT(3,1) 给出了第一材料方向上的(1, 2, 3)分量;DIRECT(1,2) , DIRECT(2,2) , DIRECT(3,2) 给出了第二材料方向,等等。对于壳单元与膜单元来说,前两个方向是在单元面内,第三个方向是法向。对于梁单元来说,这些都是不可用的。 T (3,3)
包含相对于单元基本方向的材料方向分量的方向余弦的数组。它定义了材料方向(DIRECT)关于单元基本方向的方向。对于连续单元来说,T与DIRECT是完全一致的。对于壳单元与膜单元,T 1,1 =????????,T 1,2 =?sinθ,T 2,1 =????????,T 2,2 =????????,T 3,3 =1.0,其余分量为0,其中??为法向矢量逆时针旋转角度。如果没有使用方形,那么T是单位矩阵。方向并不适用于梁单元。
CELENT
特征单元长度。对于一阶单元来说,是沿单元一条线的长度;对于二阶单元来说,是同样典型长度的一半。对于梁与桁架,是沿单元轴的特征长度。对于膜与壳单元,它是参考面的特征长度。对于轴对称单元,它只是(r,z)平面内的特征长度。
TIME (1):在当前增量步开始时的分析步时间值 TIME (2):在当前增量步开始时总的时间值 DTIME:时间增量
CMNAME:用户指定的材料名称,左对齐 ORNAME:用户指定的局部方向名称,左对齐
NFIELD:定义在材料点的场变量的数目
NSTATV:用户定义的解依赖的状态变量(详见Abaqus用户分析手册Section 18.1.1 “Allocating space‖ in ―User subroutines: overview) NOEL:单元数目 NPT:积分点数目
LAYER:层数(对于复合壳体与分层固体) KSPT:当前层内截面点的数目 KSTEP:步数 KINC:增量数
NDI:此点法向应力分量的数目 NSHR:此点切向应力分量的数目 COORD:在材料点的坐标
JMAC:必须将变量传到GETVRM实用程序来访问输出变量 JMATYP:必须将变量传到GETVRM实用程序来访问输出变量 MATLAYO:必须将变量传到GETVRM实用程序来访问输出变量 LACCFLA:必须将变量传到GETVRM实用程序来访问输出变量
实例:损伤弹性模型
下面介绍一个用户子程序USDFLD的例子。在例子中,桁架单元承受拉应力。介绍了损伤弹性模型:模量随在加载历史中发生的最大拉应变的函数而减少。最大拉应变存储为解依赖的状态变量,详见Abaqus Analysis User‘s Manual S