temperature expansion analysis

dioder

I simulated a  simple temperature expansion application, but I found the results are not the same between Ansys and Femlab, I thought maybe I was wrong in Ansys script. Do you have any idea about my script? I will appreciate your answers. Thank you very much



/PREP7   
ET,1,PLANE13
!*  
KEYOPT,1,1,4
KEYOPT,1,2,0
KEYOPT,1,3,2
KEYOPT,1,4,0
KEYOPT,1,5,0
!*  
MP,DENS,1,3.21e-18  !SiC density g/cm^3
MP,KXX,1,300e6
MP,KXX,2,200e6
!MPTEMP, ,20,100,200,300,400,500,600
!MPDATA, KXX,1, ,

MP,EX,1,4.7E8
MP,PRXY,1,0.14
MP,ALPX,1,3.7E-6
MP,EX,2,3.2E8
MP,PRXY,2,0.31
MP,ALPX,2,4.1E-6

rhole=600e-6
lprotect=100e-6
rdie1=700e-6   !800
rdie2=1000e-6
thmembrane=40e-6
thdie=525e-6
rcavity=800e-6
lcavity=200e-6
rpack=2000e-6
lpack=1800e-6

K,,0,0
CYL4,,,0,0,rdie2
CYL4,,,rdie2,0,rpack

NUMMRG,ALL, , , ,LOW
NUMCMP,ALL
LSEL, , , ,5,8,1,0   
LESIZE,ALL, , ,24, , , , ,1   
LSEL, , , ,1,4,1,0   
LESIZE,ALL, , ,24, , , , ,1   
   

ASEL, , , ,1
aatt,1,1,1
MSHKEY,1
AMESH,ALL  
ASEL, , , ,2  
aatt,2,1,1
MSHKEY,0
AMESH,ALL   
ALLSEL,ALL  
  
/SOLU
DK,1,UX,0
DK,1,UY,0
ANTYPE,STATIC
NROPT,FULL
!IC,ALL,TEMP, 20
TREF,20


ALLSEL,ALL  
BFUNIF,TEMP,500
SOLVE   
FINISH

dioder

I did it, and the problem was meshing

inkey

/PREP7   
ET,1,PLANE13
!*  
KEYOPT,1,1,4
KEYOPT,1,2,0
KEYOPT,1,3,2
KEYOPT,1,4,0
KEYOPT,1,5,0
!*  
MP,DENS,1,3.21e-8  !SiC density g/cm^3
MP,KXX,1,300e6
MP,KXX,2,200e6
MPTEMP, ,20,100,200,300,400,500,600
MPDATA, KXX,1, ,

MP,EX,1,4.7E8
MP,PRXY,1,0.14
MP,ALPX,1,3.7E-6
MP,EX,2,3.2E8
MP,PRXY,2,0.31
MP,ALPX,2,4.1E-6

rhole=600e-6
lprotect=100e-6
rdie1=700e-6   !800
rdie2=1000e-6
thmembrane=40e-6
thdie=525e-6
rcavity=800e-6
lcavity=200e-6
rpack=2000e-6
lpack=1800e-6

K,,0,0
CYL4,,,0,0,rdie2
CYL4,,,rdie2,0,rpack

NUMMRG,ALL, , , ,LOW
NUMCMP,ALL
LSEL, , , ,5,8,1,0   
LESIZE,ALL, , ,24, , , , ,1   
LSEL, , , ,1,4,1,0   
LESIZE,ALL, , ,24, , , , ,1   
   

ASEL, , , ,1
aatt,1,1,1
MSHKEY,1
AMESH,ALL  
ASEL, , , ,2  
aatt,2,1,1
MSHKEY,0
AMESH,ALL   
ALLSEL,ALL  
  
/SOLU
D,1,ALL
ANTYPE,STATIC
NROPT,FULL
IC,ALL,TEMP, 20
TREF,20
ALLSEL,ALL  
BFUNIF,TEMP,500
SOLVE

[ 本帖最后由 inkey 于 2006-12-12 22:21 编辑 ]