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 |