[求助]收敛问题

bijiaoshuai

为什么结果不收敛？请达人指点！

width=1e-3
dx=0.1*width
sel_elem=4806
hf_cent1=600000000
hf_cent2=360000000   
step=0
/PREP7
/UNIT,SI
ET,1,SOLID70   
MP,DENS,1,6202
MPTEMP,1,25,100,500,1000,1500
MPDATA,C,1,1,423.6,453.2,488.4,537.2,587.2
MPTEMP,1,25,500,1000,1500,3000
MPDATA,KXX,1,1,2,5,50,700,800
MPDATA,ENTH,1,1,0,9.9E5,1.8E6,3.25E7,4.3E7
MPTEMP,1,25,100,500,1000,1500
MP,DENS,2,7824
MPDATA,KXX,2,1,47.68,34.16,24.02,225,235
MPDATA,C,2,1,472,480,615,700,500
CSYS,0
BLOCK,-0.005,0.005,-0.0024,0.0024,0,0.001
BLOCK,-0.004,0.004,-0.00135,0.00135,0.001,0.0013
VGLUE,1,2
NUMCMP,ALL
LESIZE,13,dx
LESIZE,14,dx
LESIZE,15,dx
LESIZE,16,dx
LESIZE,17,dx
LESIZE,18,dx
LESIZE,19,dx
LESIZE,20,dx
LESIZE,21,dx
LESIZE,22,dx
LESIZE,23,dx
LESIZE,24,dx   
mat,1
VMESH,1
MAT,2
ESIZE,0.001
MSHAPE,1,3d   
VMESH,2
IC,ALL,TEMP,25
SFA,2,,CONV,12.5,25
SFA,3,,CONV,12.5,25
SFA,4,,CONV,12.5,25
SFA,5,,CONV,12.5,25
SFA,6,,CONV,12.5,25
SFA,8,,CONV,12.5,25
SFA,9,,CONV,12.5,25
SFA,10,,CONV,12.5,25
SFA,11,,CONV,12.5,25
SFA,1,,CONV,12.5,25
FINISH
/SOLU
ANTYPE,4
TRNOPT,FULL
*do,i,0,66
sel_elem=sel_elem+1
sfe,sel_elem,6,hflux,,hf_cent1
sfe,sel_elem+1,6,hflux,,hf_cent1
sfe,sel_elem,2,hflux,,hf_cent1
sfe,sel_elem+1,2,hflux,,hf_cent1
sfe,sel_elem-1,6,hflux,,hf_cent2
sfe,sel_elem-1,2,hflux,,hf_cent2
sfe,sel_elem+2,6,hflux,,hf_cent2
sfe,sel_elem+2,2,hflux,,hf_cent2
sfe,sel_elem+80,6,hflux,,hf_cent2
sfe,sel_elem+81,6,hflux,,hf_cent2
sfe,sel_elem-2160,2,hflux,,hf_cent2
sfe,sel_elem-2159,2,hflux,,hf_cent2
sfe,sel_elem-80,6,hflux,,hf_cent1
sfe,sel_elem+1-80,6,hflux,,hf_cent1
sfe,sel_elem-80,4,hflux,,hf_cent1
sfe,sel_elem+1-80,4,hflux,,hf_cent1
sfe,sel_elem-1-80,6,hflux,,hf_cent2
sfe,sel_elem-1-80,4,hflux,,hf_cent2
sfe,sel_elem+2-80,6,hflux,,hf_cent2
sfe,sel_elem+2-80,4,hflux,,hf_cent2
sfe,sel_elem-80-80,6,hflux,,hf_cent2
sfe,sel_elem-79-80,6,hflux,,hf_cent2
sfe,sel_elem-80-2160,2,hflux,,hf_cent2
sfe,sel_elem-80-2159,2,hflux,,hf_cent2
TIME,0.0025+step*0.0025
nsubst,2
KBC,1
AUTOTS,ON
TIMINT,ON
OUTRES,ALL,ALL
step=step+1   
LSWRITE,step
*ENDDO
LSSOLVE,1,step
save
FINISH

[ 本帖最后由 bijiaoshuai 于 2008-11-21 15:27 编辑 ]

zhoulon

能不能具体讲一下模型，出错提示什么的。

bijiaoshuai

在矩形板上烧结，为什么不收敛？？急？？
BLOCK,-0.005,0.005,-0.0024,0.0024,0,0.001
BLOCK,-0.004,0.004,-0.00135,0.00135,0.001,0.0013
错误信息如下：

*** WARNING ***                         CP =     325.531   TIME= 21:43:32
Both solid model and finite element model boundary conditions have been
applied to this model.  As solid loads are transferred to the nodes or  
elements, they can overwrite directly applied loads.                    
*** WARNING ***                         CP =     353.281   TIME= 21:44:01
INITIAL CONDITIONS SPECIFIED - TUNIF IGNORED.                           
*** WARNING ***                         CP =     354.344   TIME= 21:44:02
Previous testing revealed that 3 of the 14886 selected elements violate
shape warning limits.  To review warning messages, please see the      
output or error file, or issue the CHECK command.                       
*** WARNING ***                         CP =     356.172   TIME= 21:44:04
Material property C of material 1 of element 4724 is evaluated at a     
temperature of 1537.24344, which is above the supplied temperature      
range.  Temperature range checking terminates.