『求助』热电耦合问题

这是我正在做的一个热电耦合的例子，可以就是算不出结果，请各位大侠指教!
finish
/clear
/prep7
/title,polyfuse
/COM, Structural
/COM, Thermal
/COM, Electric
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
ET,1,SOLID69
!WU
MPTEMP,1,0
MPDATA,DENS,1,,19300

MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,C,1,,134

MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,KXX,1,,163.3
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,RSVX,1,,5.65e-8

!!!!!!!!!!!!!!SIO2
MPTEMP,1,0
MPDATA,DENS,2,,2203

MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,C,2,,743.157

MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,KXX,2,,1.3808005
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,RSVX,2,,1.00E+12

!!!!!!!!!!!!!!!PTSIO2
MPTEMP,1,0
MPDATA,DENS,3,,12400

MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,C,3,,769
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,KXX,3,,158.55
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,RSVX,3,,5.10E-07

!POLYsilicon
MPTEMP,1,0
MPDATA,DENS,4,,2330

MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,C,4,,753
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,KXX,4,,400
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,RSVX,4,,1.49E-05

!silicon
MPTEMP,1,0
MPDATA,DENS,5,,2330

MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,C,5,,707.5692

MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,KXX,5,,148.12224
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,RSVX,5,,4.4

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!WU
lau=0.6 ! Length of the WU
wau=0.6 !Width of the WU
hau=1.8 !Height of WU

pitch=1.2 ! Pitch of WU in Rows of Length

L1=2.3
L2=4.5
L3=2
L4=0.7 !Length of the ploy silicon in the ptsi2

W1=9.6 !Width of the whoel
W2=0.5 !Width of the ploy silicon in the ptsi2

h1=4 !Height of sio2
h2=0.065 !Height of ptsi2
h3=0.27 !Height of ptsi2 and ploy silicon
h4=20 !Height of silicon
!moding
!!!!!!!!!!!!!!!!!!!!!!!!!!mode au
wpoff,l2+l3/2+(l1-pitch)/2,pitch/2,h4+h3
BLOCK,-lau/2,lau/2,-wau/2,wau/2,0,hau
VGEN,4,-1, , , ,pitch, , ,0
VGEN,2,-4, , ,pitch ,-pitch/2, , ,0
wpcsys,-1,0
VSYMM,X,all, , , ,0,0
wprota,,90
vsbw,all
csys,4
vsel,s,loc,z,0,pitch
VDELE,all, , ,1
allsel
cm,vAu,volu
wpcsys,-1,0

vsel,none
allsel,below,volu

k,1001,
k,1002,l4/2
k,1003,l4/2,w2/2
k,1004,,w2/2
k,1005,,w1/2
k,1006,l4/2,w1/2
k,1007,l3/2+l2,w1/2
k,1008,(l2+l3)/2,(w1+w2)/4
*AFUN,DEG
k,1009,l4/2,(w1+w2)/4+tan(15)*(l3-l4+l2)/2
k,1010,l3/2,w2/2
k,1011,l3/2
*AFUN,DEG
k,1012,(l2+l3)/2+tan(15)*(w1+w2)/4
k,1013,l3/2+l2+(l1-pitch-lau)/2,pitch+(3*wau)/2
k,1014,l3/2+l2+(l1-pitch-lau)/2,w1/2
k,1015,l3/2+l2+(l1-pitch-lau)/2,
k,1016,l1+l2+l3/2,w1/2,
k,1017,l1+l2+l3/2,0,

l,1001,1002
l,1003,1004
l,1001,1004
l,1002,1003
l,1004,1005
l,1005,1006
l,1006,1009
l,1009,1003
l,1003,1010
l,1010,1011
l,1002,1011
l,1006,1007
l,1007,1008
l,1008,1010
l,1008,1009
l,1011,1012
l,1008,1012
l,1012,1015
l,1015,1013
l,1013,1008
l,1014,1013
l,1007,1014
l,1014,1016
l,1017,1016
l,1015,1017

numcmp,all

al,49,51,52,53
al,53,145,147,148
al,51,58,57,56,54
al,57,149,150,153
al,153,152,145,58
al,152,193,154,147
al,193,194,196,197
al,197,150,205,206
al,196,217,209,208,205

ARSYM,X,all, , , ,0,0
aadd,99,108
aadd,97,106
CM,plant,AREA
VEXT,plant, , ,0,0,h3+h1+h4,,,,

allsel
VOVLAP,all
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
wpoffs,,,h4
vsbw,all
wpoffs,,,h3-h2
vsbw,all
wpoffs,,,h2
vsbw,all
wpoffs,,,hau
vsbw,all
!!!!!!!!!!!!!!!!!!!纵向分割

WPCSYS,-1,0
wpro,,,90.000000

wpoffs,,,l2+l3/2+(l1-pitch-lau)/2+lau
vsbw,all
wpoffs,,,pitch/2
vsbw,all
wpoffs,,,lau
vsbw,all

WPCSYS,-1,0
wpro,,,90.000000
wpoffs,,,-l2-l3/2-(l1-pitch-lau)/2-lau
vsbw,all
wpoffs,,,-pitch/2
vsbw,all
wpoffs,,,-lau
vsbw,all

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!横向分割
WPCSYS,-1,0
vsel,u,loc,x,-l2-l3/2,l2+l3/2
wpro,,-90.000000,
wpoffs,,,wau/2
vsbw,all
wpoffs,,,wau
vsbw,all
wpoffs,,,wau
vsbw,all
wpoffs,,,wau
vsbw,all
wpoffs,,,wau
vsbw,all
wpoffs,,,wau
vsbw,all
wpoffs,,,wau
vsbw,all
wpoffs,,,wau
vsbw,all
WPCSYS,-1,0

allsel

!!!!!!!!!!!!!!!!!!!!!!!!! 画au 的网格
/PNUM,MAT,1
/REPLOT
CMSEL,S,VAU
VATT, 1, , 1, 0
esize,0.3
MSHAPE,0,3D
MSHKEY,1
VMESH,VAU
!!!!!!!!!!!!!!!!!!!!!!!!!

allsel ! 画sio2 left&right
vsel,s,loc,z,h4+h3,h4+h3+h1
vsel,u,loc,x,-l2-l3/2,l2+l3/2
CMSEL,u,VAU
CM,Vsi1,VOLU
VATT, 2, , 1, 0
VSWEEP,Vsi1

allsel !
vsel,s,loc,z,h4+h3+hau,h4+h3+h1
vsel,r,loc,x,-l2-l3/2,l2+l3/2
CM,Vsi2,VOLU
VATT, 2, , 1, 0

MSHAPE,0,3D
MSHKEY,1
VMESH,79
VMESH,20
VMESH,25
VSWEEP,71
VSWEEP,80
VMESH,18

MSHAPE,0,3D
MSHKEY,1
VMESH,74
VMESH,75
VMESH,76
VSWEEP,38
VSWEEP,77
VMESH,73
LESIZE,313, , ,4, , , , ,0
LESIZE,314, , ,4, , , , ,0
VSWEEP,107
VSWEEP,101

!!!!!!!!!!
allsel !画sio2 和 含au层
vsel,s,loc,z,h4+h3,h4+h3+hau
vsel,r,loc,x,-l2-l3/2,l2+l3/2
CM,Vsi3,VOLU
VATT, 2, , 1, 0
VSWEEP,Vsi3

!!!!!!!!!!!!!!!!!!!!!!!!sio2inpolysilicon
allsel
vsel,s,loc,z,h4,h4+h3
vsel,r,loc,x,-l3/2-l2/3,l3/2+l2/3
vsel,u,loc,y,0,1.5*w2
cm,Vsi4,VOLU
VATT, 2, , 1, 0
VSWEEP,Vsi4

allsel
CMsel,s,Vsi1
CMsel,a,Vsi2
CMsel,a,Vsi3
cmsel,a,Vsi4
CM,Vsi,VOLU

!!!!!!!!!!!!!!!!!!!!ptsi2
allsel
vsel,s,loc,z,h4+h3-h2,h4+h3
CMSEL,U,Vsi4
vsel,u,loc,x,-l4/2,l4/2
CM,Vpt,VOLU
VATT, 3, , 1, 0
VSWEEP,Vpt

!!!!!!!!!!!!!POLYSILICON
allsel
vsel,s,loc,z,h4,h4+h3
CMSEL,U,Vsi4
CMSEL,U,Vpt
CM,VPSI ,VOLU
VATT, 4, , 1, 0
VSWEEP,VPSI

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!SILICON
allsel
vsel,s,loc,z,0,h4
VATT, 5, , 1, 0
CM,VSILICON,VOLU
VSWEEP,VSILICON
save
allsel

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

CMSEL,S,VAU
ALLSEL,BELOW,VOLU
ASEL,R,LOC,Z,h4+h3
asel,r,loc,x,0,-l3/2-l2-l1
CM,aload,AREA

CMSEL,S,VAU
ALLSEL,BELOW,VOLU
ASEL,R,LOC,Z,h4+h3
asel,r,loc,x,0,l3/2+l2+l1
CM,load0,AREA

allsel
ASEL,s,LOC,x,l1+l2+l3/2
ASEL,a,LOC,x,-l1-l2-l3/2
ASEL,a,LOC,y,w1/2
ASEL,a,LOC,z,h1+h4+h3
ASEL,a,LOC,z,0
CM,aconv,AREA
save
allsel
VLSCAL,all, , ,1/1000000,1/1000000,1/1000000, ,0,1
save

CMSEL,s,load0
NSLA,S,1
D,all,VOLT,0
CMSEL,s,aload
NSLA,S,1
CP,1,VOLT,all
*get,nmax,node,,num,max
nsel,s,,,nmax
f,all,amps,60

CMSEL,s,aconv
SFA,aconv,1,CONV,20,273+50

allsel
save

/SOLU
solve