【求教】移动热源分析中材料属性变化为什么会使问题无法求解？

hunterma

　　我仿真的问题是在硅片上用激光进行加热，首先假设材料属性不随温度变化时能有好的求解结果，但当把材料属性改变时就出现了问题，不知道如何解决。我的命令流如下：
　　tm=0.5
tinc=0.02
esize=0.1e-3
v=0.005
R=0.1e-3
pi=3.1415926
length=0.0005
width=0.01
thickness=0.5e-3
!max_time=time/tinc+1=26
!max_x=x/esize+1=0.5e-3/0.1e-3+1=6
!max_y=y/esize+1=1e-2/0.1e-3+1=101
max_time=26
max_x=6
max_y=101
power=40
flux_max=40/(pi*R*R)


/prep7
/units,si
et,1,solid70
!properties of silicon
mptemp,1,300,400,500,600,700,800,900,1000,1100,1200,1300,1400,1500,1600
mpdata,kxx,1,1,135,97.8,76.2,62.1,52.3,45,39.4,35.1,31.5,28.6,26.1,24,22.3,20.7
mp,dens,1,2.4e3
mp,c,1,1,715,740,765,789,814,839,864,888,913,938,962,987,1012,1037

blc4,0,0,length,width,thickness !定义si实体模型
blc4,length,0,0.0045,width,thickness

vglue,all
numcmp,all

lsel,s,,,6,8,1
lesize,all,0.1e-3
lsel,s,,,9
lesize,all,,,4,2
lsel,s,,,10
lesize,all,,,4,0.5
lsel,s,,,17,20,1
lesize,all,,,15,2
allsel
lsel,s,,,15,16
lesize,all,,,4,0.5
lsel,s,,,13,14
lesize,all,width/30
allsel
mat,1
mshkey,1
vmesh,1
allsel

mat,1
mshkey,1
vmesh,2

allsel
fini
/psf,hflu, ,2,0  !Shows surface load symbols on model displays.
/solu
antype,trans,new
trnopt,full
!toffst,273
!thopt,quasi    !using the quasi solver option
!eqslv,jcg
tunif,300
tref,300
kbc,1
*SET,DELT,1  
timint,on
autots,off

*DIM,laser,TABLE,MAX_X,MAX_Y,MAX_TIME,X,Y,TIME
*DO,K,1,MAX_TIME,1
  *DO,I,1,MAX_X,1
   laser(I,0,K)=(I-1)*ESIZE
   *ENDDO
   *DO,J,1,MAX_Y,1
     laser(0,J,K)=(J-1)*ESIZE
   *ENDDO
  *DO,I,1,MAX_X,1
  *DO,J,1,MAX_Y,1
     XCENTER=0
     YCENTER=V*(K-1)*TINC+3e-3
     DISTANCE=SQRT(((I-1)*ESIZE-XCENTER)**2+((J-1)*ESIZE-YCENTER)**2)
       *IF,DISTANCE,LE,R,THEN
         laser(I,J,K)=flux_max*EXP(-DISTANCE**2/(R**2))
       *ELSE
         laser(I,J,K)=0
       *ENDIF
  *ENDDO
  *ENDDO
     laser(0,0,K)=(K-1)*TINC
   
*ENDDO


*do,t,1e-6,tm,tinc
!sfdele,all,hflux
time,t
autots,off  !
deltim,,,,1
kbc,0
tsres,erase
!nsubst,10,20,5
nsel,s,loc,x,0,length
nsel,r,loc,z,thickness,thickness

sf,all,hflux,%laser%
allsel
solve
outres,all,all
*enddo
finish
/post1
set,last
plnsol,temp

hunterma

当属性不变化时的温度分布如下

hunterma

我尝试用ansys thermal analysis guide 推荐在材料非线性使用的thopt,quasi  求解器时也是出现问题，真的不知道怎么解决这个问题了。希望大家知道的能帮帮忙，在此谢过了！

hunterma

修改了一下求解设置/solu
antype,trans,new
trnopt,full
!toffst,273
tunif,300
tref,300
!timint,on
autots,-1
cnvtol,HEAT,2e-6
SOLCONTROL,ON
现在出现这种问题：The calculated reference HT FLOW CONVERGENCE VALUE= 7.377837119E-08 is   less than the default minimum of 1E-6 (MINREF field of CNVTOL           
command).  Check results carefully.      
有知道原因的吗？我已经把cnvtol,heat,2e-6了怎么还是出现这个问题呢

hunterma

问题已解决！

yuxiaoyuh1138

敢问 ，错误原因可否相告？谢谢！

glaciery

能告知错误原因了，我也出现了类似的错误，谢谢了

sigh08

5# hunterma
不知道楼主能否告知解决办法？