这是应力部分的一些程序,请问哪里出错了?
tunif,25 !设置温度升高25度,产生应力
tref,25 !参考温度25度
!antype,trans !瞬态分析类型
!*******************************
!位移边界条件
!*******************************
tref,25
nsel,s,loc,x,-0.09,0
nsel,r,loc,y,0
d,all,,,,,,uy,,,,,
allsel,all
nsel,s,loc,y,0,0.36
nsel,r,loc,x,-0.09
d,all,,,,,,ux,,,,,
DA,6,SYMM
allsel,all
nsel,s,loc,x,-0.09,0
nsel,r,loc,Z,0
d,all,,,,,,uz,,,,,
!*******************************
!非线性求解
!*******************************
/sol
allsel
antype,STATIC !稳态分析
!trnopt,full !完全分析
lumpm,on !集中质量近似形式
nlgeom,on !大应变开关
autots,on !自动时间步长
lnsrch,on
neqit,150
TREF,25 !参考温度
AUTOTS,ON !打开自动时间步长
ESEL,S,TYPE,,2 !选择类型为2的单元节点
NSLE,S !选择所有单元
!IC,ALL,TEMP,25 !指定初始温度
timint,on
*do,i,0,2,0.2
OUTRES,ALL,LAST
time,i+0.0001
kbc,0
NSUBST,20,200,10 !子步数,默认150,若不收敛,子步数范围为10-10000
LDREAD,TEMP,,,i+0.0001,0,'LONG-1001','rth'
solve
*enddo
![薄板tig焊的残余应力,做出来是图中的样子,非常小,请问是程序哪里出错了吗?]()
LONG-1005000.jpg |