模拟过程中出现了如下错误:
The value of TEMP at node 434 is 1.594205862E+25. It is greater than
the current limit of 1.E+25. This generally indicates rigid body
motion as a result of an unconstrained model. Verify that your model
is properly constrained.
查看节点温度,时该节点的温度超高,且该区域的温度都很高。以下是我的命令流,求高手帮忙看看,哪里出了问题,谢了~
/CLEAR !清除数据库,新建文件
/FILNAME, SY_FANG_02 !改变工作名
!开始热分析
!建立有限元模型
/NERR,,10000000
/PREP7 !进入前处理器
ET,1,PLANE77 !平面热单元
ET,2,SOLID5 !空间热单元
V=0.007 !焊接线速度
LENGTH=0.035 !钢板沿焊缝方向长度
L_NUMBER=20 !沿焊缝方向单元段数
TIME=0.01
TIME_INC=LENGTH/(L_NUMBER*V) !第一层焊缝时间增量
!
!
!
!
MPTEMP,1,0,300,600,1000,1300,2000 !定义材料模型,温度表
MPDATA,DENS,1,1,7859,7770,7659,7600,6750,6750 !密度
MPDATA,DENS,2,1,7859,7770,7659,7600,6750,6750
MPDATA,C,1,1,450,514,967,1370,2800,2000 !比热
MPDATA,C,2,1,450,514,967,1370,2800,2000
MPDATA,KXX,1,1,67,53,39,31,20,20 !传导系数
MPDATA,KXX,2,1,67,53,39,31,20,20
MPTEMP,1,0,300,600,1000,1300,1400,2000 !温度表
MPDATA,EX,1,1,2E11,1.86E11,1.35E11,2E10,2E7, 2E7, 2E7 !弹性模量
MPDATA,EX,2,1,2E11,1.86E11,1.35E11,2E10,2E7, 2E7, 2E7
MPDATA,PRXY,1,1,0.33,0.33,0.35,0.36,0.4,0.4 !泊松比
MPDATA,PRXY,2,1,0.33,0.33,0.35,0.36,0.4,0.4
MPDATA,ALPX,1,1,12E-6,12.6E-6,13.0E-6,13.4E-6,13.7E-6,14E-6,13.7E-6 !线膨胀系数
MPDATA,ALPX,1,1,12E-6,12.6E-6,13.0E-6,13.4E-6,13.7E-6,14E-6,13.7E-6
TB,BKIN,1,5
TBTEMP,15
TBDATA,1,280E6,2E9 !屈服极限,切向模量
TBTEMP,800
TBDATA,1,10E6,1E7
TBTEMP,2500
TBDATA,1,10E6,1E6
MPTEMP,,,,,,,, !初始化随温度变化的材料参数表格
MPTEMP,1,20 !定义表格的温度数值
MPTEMP,2,250
MPTEMP,3,400
MPTEMP,4,500
MPTEMP,5,750
MPTEMP,6,800
MPTEMP,7,1000
MPTEMP,8,1200
MPTEMP,9,1500
MPTEMP,10,1700
MPTEMP,11,2500
MPDATA,HF,1,,100 !定义对应温度下的对流换热系数数值
MPDATA,HF,1,,350
MPDATA,HF,1,,480
MPDATA,HF,1,,520
MPDATA,HF,1,,1000
MPDATA,HF,1,,1050
MPDATA,HF,1,,1500
MPDATA,HF,1,,2000
MPDATA,HF,1,,3000
MPDATA,HF,1,,3100
MPDATA,HF,1,,3500
MPTEMP,,,,,,,, !初始化随温度变化的材料参数表格
MPTEMP,1,20 !定义表格的温度数值
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,400
MPTEMP,5,500
MPTEMP,6,1000
MPTEMP,7,1100
MPTEMP,8,1500
MPTEMP,9,1700
MPTEMP,10,2500
MPDATA,HF,2,,120 !定义对应温度下的对流换热系数值
MPDATA,HF,2,,210
MPDATA,HF,2,,300
MPDATA,HF,2,,460
MPDATA,HF,2,,540
MPDATA,HF,2,,1200
MPDATA,HF,2,,1200
MPDATA,HF,2,,2800
MPDATA,HF,2,,3000
MPDATA,HF,2,,3100
!
TBCOPY,BKIN,1,2 !复制材料BKIN数据
MP,REFT,1,20 !参考温度
MP,REFT,2,20
!
!
!
K,1,0,0
K,2,0.015,0
K,3,0.015,0.004
K,4,0,0.004
K,5,0.002,0.004
K,6,0.002,0
A,1,6,5,4
A,6,2,3,5
AGLUE, ALL
LESIZE,1,,,3
LESIZE,3,,,3
LESIZE,4,,,4
LESIZE,2,,,4
LESIZE,6,,,4
LESIZE,5,,,10,4
LESIZE,7,,,10,0.25
VEXT,1,2,1,0,0,0.035
VGLUE,ALL
LESIZE,14,,,30
TYPE,2
MAT,1 !2号材料为母材,1号材料为焊缝
REAL,
VSWEEP,1
TYPE,2
MAT,2
REAL,
VSWEEP,2
!
FINISH
/SOLU !进入求解器
DA,7,SYMM !对称面A12施加对称载荷!已改
DL,5,,UZ
DL,1,,UZ
DL,4,,UX
DL,12,,UX
TUNIF,20,
TOFFST,273,
!
!
!!!!!!!!求解设置!!!!!!!!!!!!!!!!!!!!!!
ANTYPE,4 !瞬态分析
TRNOPT,FULL
LUMPM,0
NLGEOM,1
PSTRES,1 !打开应力刚化
NSUBST,5,20,1
OUTRES,ERASE
OUTRES,ALL,ALL !输出控制
AUTOTS,1 !打开自动时间步长
PRED,ON,,ON
LNSRCH,1 !打开线性搜索
!
!开始激活焊缝单元,进行瞬态热分析
TIMINT,0,MAG
TIMINT,0,ELECT
TIMINT,ON !打开时间积分
!
!
!!!!!!!!!!!!!!定义数组!!!!!!!!!!!!!!!!!!!!!!!
*DIM,QQ,ARRAY,8 !定义数组,存储单元节点热流密度
*DIM,E_NODE,ARRAY,20 !定义数组,存储单元所属节点编号
*GET,ELEM_MAX,ELEM,0,NUM,MAX !获取最大单元号
*GET,ELEM_MIN,ELEM,0,NUM,MIN !获取最小单元号
*GET,NODE_MAX,NODE,0, NUM,MAX !获取最大节点号
*GET,NODE_MIN,NODE,0,NUM,MIN !获取最小节点号
*DIM,NODE_AVE,ARRAY,NODE_MAX !定义数组,存储节点相邻单元数量
!
!
!
!
ESEL,S,MAT,,1 !选择焊缝单元!!!!!!!!!!!!!!!!!
EKILL,ALL !杀死单元!!!!!!!!!!!!!!!
!
!
alls!!!!!!!!!!!!!!
!
!
R=0.007 !电弧有效加热半径
Q=140*20*0.8 !输入热量
QMAX=3*Q/(3.1415926*R*R) !变量
XC= 0 !第一层焊缝质心X坐标
YC= 0.007 !第一层焊缝质心Y坐标
!
!
!
!!!!!!!!!!!!从头至尾激活焊缝,开始
*DO,J,1,L_NUMBER
alls
SFEDELE,ALL,ALL,HFLUX !删除以前施加的热流密度
VSEL,S,,,1 !选择第KK层焊缝所对应的体!!!!!!!!!!!!!!!!!!!!!!
ESLV,S,1 !选择体上单元即第KK层焊缝单元!!!!!!!!!!!!!!!!!!!!
NSEL,R,LOC,Z,(J-1)*LENGTH/L_NUMBER,(J+1)*LENGTH/L_NUMBER!!!!!!!!!!!!!!!!!!!!!
!选择第J行单元的中间节点!!!!!!!!!!!!!!!!!!!!!
ESLN,R !选择第KK层第J行单元!!!!!!!!!!!!!!!!
EALIVE,ALL !激活单元!!!!!!!!!!!!!!!!!
ALLS !!!!!!!!!!!!!!!
ZC=J*LENGTH/L_NUMBER !第KK层第J行单元质心Z坐标
*DO,I,ELEM_MIN,ELEM_MAX !计算施加单元热流密度,循环开始
*IF,ESEL(I),EQ,1,THEN
ESEL,S,,,I !选择第I个单元
*DO,K,1,8 !计算第I单元前8个节点热流密度
*GET,E_NODE(K),ELEM,I,NODE,K !获取第I个单元第K个节点的编号
SS=E_NODE(K)
DISP=SQRT((NX(SS)-XC)*(NX(SS)-XC)+(NY(SS)-YC)*(NY(SS)-YC)+(NZ(SS)-ZC)*(NZ(SS)-ZC))
EEE=3*DISP*DISP/(R*R) !变量
*IF,DISP,LT,R,THEN
QQ(K)=QMAX*EXP(-EEE)/8 !第K个节点的热流密度
*ELSE
QQ(K)=0
*ENDIF
*ENDDO
*IF,QQ(1)+QQ(2)+QQ(3)+QQ(4),NE,0,THEN
SFE,I,1,HFLUX,,QQ(1), QQ(2), QQ(3),QQ(4) !在第I个单元第一个面上加热流密度
*ENDIF
*IF,QQ(1)+QQ(2)+QQ(5)+QQ(6),NE,0,THEN
SFE,I,2,HFLUX,,QQ(1), QQ(2), QQ(6),QQ(5) !在第I个单元第二个面上加热流密度
*ENDIF
*IF,QQ(6)+QQ(2)+QQ(3)+QQ(7),NE,0,THEN
SFE,I,3,HFLUX,,QQ(2), QQ(3), QQ(7),QQ(6) !在第I个单元第三个面上加热流密度
*ENDIF
*IF,QQ(1)+QQ(4)+QQ(5)+QQ(7),NE,0,THEN
SFE,I,5,HFLUX,,QQ(1), QQ(4), QQ(7),QQ(5) !在第I个单元第五个面上加热流密度
*ENDIF
*IF,QQ(5)+QQ(6)+QQ(7)+QQ(8),NE,0,THEN
SFE,I,6,HFLUX,,QQ(5), QQ(6), QQ(7),QQ(8) !在第I个单元第六个面上加热流密度
*ENDIF
ALLS!!!!!!!!!!!!!!!
*ENDIF
*ENDDO !计算施加单元热流密度,循环结束
ALLS
!
!
!
TIME=TIME+TIME_INC !时间变量
!
TIME,TIME !设定分析时间
SOLVE
*ENDDO
SAVE !保存
FINI |