求救:一热应力计算的问题（请大家进来看看）

fantry2006

我在进行热应力耦合计算的时候,算到中途的时候
弹出了许多警告信息，越来越多:
/COM,ANSYS RELEASE  9.0    UP20041104       17:27:10    07/18/2006

*** WARNING ***                         CP =       4.391   TIME= 17:27:21
The ETCHG command has changed some element types.  Please examine not   
only the KEYOPTs, but also the real constants, material properties,     
boundary conditions, and loadings, and then modify as needed.           
/COM,ANSYS RELEASE  9.0    UP20041104       22:03:37    07/19/2006

*** WARNING ***                         CP =    2582.359   TIME= 22:50:26
There are 1 small equation solver pivot terms.                          

*** WARNING ***                         CP =    8778.984   TIME= 00:37:13
There are 1 small equation solver pivot terms.                          

*** WARNING ***                         CP =    8819.703   TIME= 00:37:54
There are 2 small equation solver pivot terms.                          

*** WARNING ***                         CP =    9532.203   TIME= 00:50:10
There are 1 small equation solver pivot terms.                          

*** WARNING ***                         CP =    9579.875   TIME= 00:50:59
There are 2 small equation solver pivot terms.                          

*** WARNING ***                         CP =   13850.125   TIME= 02:04:34
There are 1 small equation solver pivot terms.                          

*** WARNING ***                         CP =   16687.906   TIME= 02:53:24
There are 1 small equation solver pivot terms.                          

*** WARNING ***                         CP =   19428.266   TIME= 03:40:31
There are 1 small equation solver pivot terms.                          

*** WARNING ***                         CP =   19470.312   TIME= 03:41:13
There are 2 small equation solver pivot terms.                          

*** WARNING ***                         CP =   31940.234   TIME= 07:15:51
There are 1 small equation solver pivot terms.                          

*** WARNING ***                         CP =   39735.094   TIME= 09:30:26
There are 1 small equation solver pivot terms.                          

*** WARNING ***                         CP =   39789.297   TIME= 09:31:22
There are 4 small equation solver pivot terms.                          

*** WARNING ***                         CP =   39836.734   TIME= 09:32:10
There are 7 small equation solver pivot terms.                          

*** WARNING ***                         CP =   39884.047   TIME= 09:32:59
There are 5 small equation solver pivot terms.                          

*** WARNING ***                         CP =   39932.578   TIME= 09:33:48
There are 9 small equation solver pivot terms.                          

*** ERROR ***                           CP =   39932.828   TIME= 09:33:48
One or more elements have become highly distorted.  Excessive           
distortion of elements is usually a symptom indicating the need for     
corrective action elsewhere.  Try ramping the load up instead of step   
applying the load (KBC,1).  You may need to improve your mesh to        
obtain elements with better aspect ratios.  Also consider the behavior  
of materials, contact pairs, and/or constraint equations.  If this      
message appears in the first iteration of first substep, be sure to     
run shape checking of elements.   
继而就不收敛了
　我的非线形选项的设置如下：
ANTYPE,,REST
TRNOPT,FULL
NLGEOM,ON                !包含预应力效应
NSUBST,3,0,0      
KBC,1                  
CUTCONTROL,PLSLIMIT,10
CUTCONTROL,CRPLIMITexp,0.2,0
CUTCONTROL,CRPLIMITimp,0.25,1  
CUTCONTROL,NPOINT,15
LNSRCH,ON               
TIME,81*0.25
NLGEOM,ON      
SSTIF,1                  
!NROPT,AUTO,,OFF         
NROPT,FULL,,ON           
EQSLV,,,0,               
PRECISION,0              
TOFFST,Wtoffst
OUTRES,ALL,LAST         
SOLVE
　　请各位高手能否指点一下
我都做了好多遍了总是这个结果，郁闷的很，等着大家救火呢，哪位好心人能否给予无私的帮助啊
　　谢谢拉

jun722

KBC=0 代替KBC=1

zhkr2001

原因很多，1加载步太大，也就是载荷施加增量太大；2时间子步太大；3单元太大，畸变太多。个人感觉最大原因可能是载荷步施加的太快，如果是瞬态分析，应减小载荷施加增量，如是静态（稳态），则可减小单元尺寸，优化单元形状！