- 积分
- 0
- 注册时间
- 2010-9-14
- 仿真币
-
- 最后登录
- 1970-1-1
|
很显然结果不对,应该是在热分析的时候加载以及约束的数值或条件不对,不知道空气四周的温度约束是否必要?
/prep7
cyl4,0.020,,0.001
agen,10,1,,,0.003
rectng,0,0.050,0.013,0.014
rectng,0.049,0.050,0.014,0.030
aadd,11,12
rectng,0,0.049,0.014,0.025
rectng,0,0.080,-0.030,0.070
aovlap,all
numcmp,area
! Electromagnetic model
*SET,freq,25000 ! frequency (Hz.)
*SET,pi,4*atan(1) ! pi
*SET,cond,.392e7 ! maximum conductivity
*SET,muzero,4e-7*pi ! free-space permeability
*SET,mur,200 ! maximum relative permeability
*SET,skind,sqrt(1/(pi*freq*cond*muzero*mur)) ! skin depth
*SET,ftime,300 ! final time
*SET,tinc,.05 ! time increment for harmonic analysis
*SET,time,0 ! initialize time
*SET,delt,.01 ! maximum delta time step
emunit,mks
! Electromagnetic model
et,1,53,,,1 ! PLANE53, axisymmetric, AZ dof
et,2,53,,,1
emunit,mks ! set magnetic units
mp,murx,1,1 ! air relative permeability
mp,murx,3,1 ! coil relative permeability
mp,murx,4,1
mptemp
mptemp,1,25.5,160,291.5,477.6,635,698 ! temps for relative permeability
mptemp,7,709,720.3,742,761,1000
mpdata,murx,2,1,200,190,182,161,135,104 ! steel relative permeability
mpdata,murx,2,7,84,35,17,1,1
mptemp
mptemp,1,0,125,250,375,500,625 ! temps for resistivity
mptemp,7,750,875,1000
mpdata,rsvx,2,1,.184e-6,.272e-6,.384e-6,.512e-6,.656e-6,.824e-6
mpdata,rsvx,2,7,1.032e-6,1.152e-6,1.2e-6 ! steel resistivity
mptemp
mptemp,1,0,730,930,1000 ! temps for conductivity
mpdata,kxx,2,1,60.64,29.5,28,28
mptemp ! temps for enthalpy
mptemp,1,0,27,127,327,527,727
mptemp,7,765,765.001,927
mpdata,enth,2,1,0,91609056,453285756,1.2748e9,2.2519e9,3.3396e9
mpdata,enth,2,7,3.548547e9,3.548556e9,4.3520e9
mp,emis,2,.68 ! emissivity
asel,s,,,13
aatt,1,,2
asel,s,,,11
aatt,4,,2
asel,s,,,12
aatt,2,,1
asel,s,,,1,10,1
aatt,3,1,2
asel,s,,,all
asel,s,,,1,10,1
smrt,1
mshape,0,2d
mshkey,0
amesh,all
asel,s,,,12
lsla,r
lesize,all,0.0002
mshape,0,2d
mshkey,0
amesh,12
asel,s,,,11
lsla,r
lesize,all,0.0002
mshape,0,2d
mshkey,0
amesh,11
asel,s,,,13
lsla,r
lesize,all,0.00002
mshape,0,2d
mshkey,0
amesh,13
nsel,s,loc,x,0
nsel,a,loc,x,0.08
nsel,a,loc,y,-0.03
nsel,a,loc,y,0.07
d,all,az,0
allsel,all
allsel,all
esel,s,mat,,3
bfe,all,js,,,,3.18e6
allsel,all
finish
/solu
antyp,harm ! Emag analysis options
harfrq,freq
solve
save
finish
!theraml
/prep7
etchg,mtt
mptemp
mptemp,1,0,730,930,1000 ! temps for conductivity
mpdata,kxx,2,1,60.64,29.5,28,28
mptemp ! temps for enthalpy
mptemp,1,0,27,127,327,527,727
mptemp,7,765,765.001,927
mpdata,enth,2,1,0,91609056,453285756,1.2748e9,2.2519e9,3.3396e9
mpdata,enth,2,7,3.548547e9,3.548556e9,4.3520e9
mp,emis,2,.68 ! emissivity
mp,kxx,1,0.023 ! 空气水及铜的密度、导热系数和比热
mp,c,1,1.4
mp,dens,1,1.29
mp,kxx,3,270
mp,c,3,390
mp,dens,3,8900
mp,kxx,4,0.61
mp,c,4,4200
mp,dens,4,1000
finish
/solu
LDREAD,HGEN,,, ,0,'file','rmg',' '
allsel
asel,s,,,12
nsel,r,loc,y,0.013
sf,all,rdsf,0.68,1
allsel,all
asel,s,,,12
nsel,r,loc,x,0.050
sf,all,rdsf,0.68,1
allsel,all
stef,5.67e-8 ! Stefan-Boltzman constant
radopt,,0.01
spctemp,1,25 ! Space temperature
v2dopt,1
esel,s,mat,,2 ! select billet material
bfe,all,fvin,,1 ! define volumetric interface
allsel,all
allsel
asel,s,,,12 !铁区表面对流系数为70
lsel,r,,,all
lsel,u,loc,x,0
nsll,r,,,all
sf,all,conv,70,25
allsel
nsel,s,loc,x,0 !对称轴约束对流为0
sf,all,conv,0,25
lsel,s,loc,y,0.025 !水上表面对流系数为0.54
nsll,r,,,all
sf,all,conv,0.54,25
allsel
asel,s,,,1,10,1 !线圈四周对流系数为100
lsla,r,,,all
nsll,r,,,all
sf,all,conv,100,25
allsel !空气场外围对流系数为0.024
lsel,s,loc,x,0.08
lsel,a,loc,y,-0.03
lsel,a,loc,y,0.07
nsel,r,,,all
sf,all,conv,0.024,25
allsel !空气场外围约束温度为25度
nsel,s,loc,x,0.08
nsel,a,loc,x,0
nsel,a,loc,y,-0.03
nsel,a,loc,y,0.07
d,all,temp,25
allsel
antype,trans ! Thermal analysis options
trnopt,full
toffst,273
tunif,25 ! initial uniform temperature
time,ftime
autots,on ! auto time-stepping
deltim,1,.5,1,on ! time step control
kbc,1
outres,all,all
solve
save |
|