求助，，简单热应力问题出现了奇怪的结果，请高手指点
小弟在进行热应力分析时，计算得到了很奇怪的结果，想请教这方面的达人给指点一下问题是这样的：
一个长20cm，宽3mm，厚1.8mm的材料，上表面受激光辐照，加热流密度120W/cm2.两端面固定。
采用载荷传递的方法计算时，结构分析的时间步较大时出现small equation solver pivot terms的信息，在改小时间步长后，不再出现这样的警告，只是计算结果却是变形不是向着受辐照的方向，温度低的一侧突出去了，与常理有些不合。
在将板的宽度变为20cm时，板的变形方向又是向温度高的一侧突起。
下面是我用的方法，请各位达人给看看，是不是哪儿设置错了。
/prep7
ET,1,SOLID70
*SET,a,0.1! 板长度200mm
*SET,b,0.0015! 板宽度3mm
*SET,h,0.0018! 厚度1.8mm
! 材料参数
*SET,Toffset,273
MP,DENS,1,2800
MP,KXX,1,121
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,C,1,,900
!* 划分网格
BLOCK,-a,a,-b,b,0,h,
lsel,s,,,2
lsel,a,,,4,5,1
lsel,a,,,7
LESIZE,all, , ,80, , , , ,1
lsel,s,,,9,12,1
LESIZE,all, , ,4, , , , ,1!厚度
lsel,s,,,1,3,2
lsel,a,,,6,8,2
LESIZE,all, , ,4, , , , ,1!宽度
vsel,s,,,1
VMESH,all
NUMCMP,ALL
FINISH
/SOL
!* 设置求解选项并求解
ANTYPE,TRANS !*Perform a transient analysis
SOLCONTROL,0
KBC,1
DELTIM,0.0025
TIME,0.2 ! TIME AT END OF TRANSIENT
OUTRES,,ALL
AUTOTS,Off
*SET,heatflux,120*10000
nsel,s,loc,z,h
sf,all,hflux,heatflux
allsel
TUNIF,293
allsel
SOLVE
FINISH
SAVE
/PREP7
allsel
etchg,tts
!*
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,68e9
MPDATA,PRXY,1,,0.3
MPTEMP,,,,,,,,
MPTEMP,1,0
UIMP,1,REFT,,,
MPDATA,ALPX,1,,21.4*1e-6
asel,s,,,5
DA,all,UX,0
DA,all,UY,0
DA,all,UZ,0
asel,s,,,6
DA,all,UX,0
DA,all,UY,0
DA,all,UZ,0
allsel
FINISH
/SOL
antype,static
nlgeom,1 ! 大变形
LNSRCH,1
tref,293 ! 参考温度
DELTIM,0.0025,0,0
EQSLV,-1
allsel
*do,it,0.0025,0.2,0.0025
time,it
ldread,temp,,,it,,,rth
solve
*enddo
FINISH
”small equation solver pivot terms“
表明系数矩阵中出现小主元，最大的可能是约束不足，另外还有可能是网格质量太差等
页:
[1]