求助：三维六面体试件热应力分析出错，如何解决

xukyle

一三维六面体试件，尺寸30mm×10mm×10mm，顶底面位移固定，底面温度200度，顶面温度20度，求在试件内部产生的热应力。下面是我的apdl流程文件，出现错误信息：Large negative pivot value(-0.458137512) in Eqn.system.
May be because of a bad temperature-dependent material property used in the model.   请教各位高手，该如何解决？

/BATCH  
finish
/clear
/PREP7  
!参数定义*  
ET,1,SOLID70
MP,KXX,1,0.0234
BLOCK,0,30,0,10,0,10,   

/RGB,INDEX,100,100,100, 0   
/RGB,INDEX, 80, 80, 80,13   
/RGB,INDEX, 60, 60, 60,14   
/RGB,INDEX, 0, 0, 0,15   

!网格划分
Lsel,s,line,,1,3,2  
LESIZE,all, , ,10   
Lsel,s,line,,6,10,2
LESIZE,all, , ,10   
Lsel,s,line,,11,12,1
LESIZE,all, , ,10   
Lsel,s,line,,2,4,2  
LESIZE,all, , ,30   
Lsel,s,line,,5,7,2  
LESIZE,all, , ,30   
vmesh,1   
FINISH  

/solu
antype,static

!施加温度载荷
nsel,s,loc,x,0
d,all,temp,100
nsel,s,loc,x,30
d,all,temp,10
nsel,all

Outpr,basic,all
solve
finish

/PREP7
!单元类型转换
ETCHG,TTS   
!单元力学参数
MP,EX,1,5e10
MP,PRXY,1,0.25   
MP,ALPX,1,5e-6

nsel,s,loc,x,30
cp,7,ux,all
nsel,s,loc,x,0
cp,8,ux,all
nsel,all
finish
  
/SOL
!*  
ANTYPE,0

!施加边界条件
nsel,s,loc,x,0
d,all,ux,0
nsel,s,loc,x,30
d,all,ux,0

tref,20   !定义参考温度
LDREAD,TEMP,,, ,,,rth
allsel
solve
finish

yatto

帮你分析了下有点约束过度，产生奇异解了。
nsel,s,loc,x,30
cp,7,ux,all
！nsel,s,loc,x,0
！cp,8,ux,all
nsel,all
finish
  
/SOL
!*  
ANTYPE,0

!施加边界条件
nsel,s,loc,x,0
d,all,ux,0
！nsel,s,loc,x,30
！d,all,ux,0

屏蔽掉以这几条就行了。
不知道你具体的位置约束条件是怎样的，建议在x=0处加全约束。

我现在计算也碰上产生奇异解的情况了，不过问题还没解决，郁闷ing。

[ 本帖最后由 yatto 于 2006-11-20 11:31 编辑 ]