请教大牛 热电耦合瞬态

karengao

问题描述：模型分为三部分，中间是一根银导体，该导体上开了两个狭颈，一个方的一个圆的，外面两部分是沙子，加载一个正弦变化的斜坡载荷，计算方形狭劲的发热，要求算到银的熔点1223K 。
      存在的问题：程序算出的温度分布形状不对称，在方形狭劲的中心本应温度分布最高，但是出先了灰色区域，增加了迭代次数后，仍然存在这个不对称的温度分布，请教这是什么原因造成的？谢谢！~

流文件：

  /prep7  
/title,fuse_thermal  
et,1,solid69
*SET,d,1.62e-3
*SET,h,0.35e-3
*SET,r,0.6e-3  
*SET,dh,0.81e-3
*SET,lh,1.18e-3   
*SET,df,0.81e-3
*SET,lf,2e-3
*SET,a,9e-3
*SET,c,4e-3
*SET,b,3e-3
mp,kxx,1,418.68
mp,dens,1,10490
mp,c,1,232
mp,rsvx,1,15.86e-9,0.060268e-9
mp,kxx,2,1.83
mp,dens,2,2230
mp,c,2,710
mp,rsvx,2,2.05e4
block,0,b+a+(lh-lf)/2,0,d,0,h
block,b+a+(lh-lf)/2,b+a+(lh+lf)/2,df/2,1.5*df,0,h   
block,b+a+(lh+lf)/2,b+c+a+(lh+lf)/2,0,d,0,h
wpave,b+lh/2,(d-dh-2*r)/2,0
cylind,0,r,0,h,0,360
wpave,b+lh/2,d-(d-dh-2*r)/2,0   
cylind,0,r,0,h,0,360
vsbv,1,5
vsbv,6,4
wpave,0,0.81e-3,0.175e-3,
WPROTA,0,0,-90
cylind,0,1.62e-3,-(b+c+a+(lh+lf)/2),0
cylind,1.62e-3,3.24e-3,-(b+c+a+(lh+lf)/2),0
WPROTA,0,0,90
allsel,all
vovlap,all
aglue,all
allsel,all
numcmp,all

vsel,s,,,4
vsel,a,,,1
vsel,a,,,3
vatt,1,,1
vsel,s,,,5
vsel,a,,,2
vatt,2,,1
vsel,s,mat,,,1
esize,2e-4
vsweep,all
vsel,s,volu,,5
mshkey,0
mshape,1,3D
esize,7e-4
vmesh,5
vsel,s,volu,,2
mshkey,0
mshape,1,3D
smrtsize,6
vmesh,2

/solu
asel,s,,,13
nsla,s,1
cp,1,volt,all
allsel,all
asel,s,,,4
nsla,s,1
d,all,volt,0
allsel,all
ic,all,temp,293
t=1e-4
c=0
n=10
pi=2*asin(1)
con=6e-3/n
nsubst,2
neqit,20
allsel,all
*creat,load
time,t
KBC,0
autots,on
i=1080*1.414*sin(100*pi*t)
nsel,s,,,1651
f,1651,amps,i
t=t+con
c=c+1
allsel,all
antype,4
eqslv,sparse
ncnv,0
trnopt,full
solve
*end
*do,i,1,n
*use,load
*enddo
finish
save
/post26
NSOL,2,6661,TEMP, ,max   
PLVAR,2, , , , , , , , , ,