求助，，简单热应力问题出现了奇怪的结果，请高手指点

zhang0925

小弟在进行热应力分析时，计算得到了很奇怪的结果，想请教这方面的达人给指点一下
问题是这样的：
一个长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

Chelsea

”small equation solver pivot terms“
表明系数矩阵中出现小主元，最大的可能是约束不足，另外还有可能是网格质量太差等