一个线圈通过感应涡流加热钢板的模型,这个问题已经困惑我很久了。。也做了很多次的模拟,但是却一直弄不出速度效应来,命令流文件如下。神曲了划分网格的部分。该程序完全可以运算,但算完之后却没有速度效应出来,和静止的居然一摸一样。。
shpp,off
vovlap,all
VGLUE,ALL
numcmp,all
EMUNIT,MKS !
*AFUN,DEG !
ET,1,SOLID117 ! 空气
ET,2,SOLID117,1,1 ! 钢板属性
r,2
rmore, , ,10
ET,3,SOLID117,1 ! 线圈属性
mp,murx,1,1 ! air relative permeability
mp,murx,3,1 ! coil relative permeability
mp,rsvx,1,1e-9
mptemp,1,25.5,160,291.5,477.6,635,698 ! temps for relative permea
mptemp,7,709,720.3,742,761,1000
mpdata,murx,2,1,200,190,182,161,135,104 ! steel relative permeabili
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 !
*SET,FTIME,0.2 !!!
*SET,TINC,0.05 !
*SET,TIME,0 !
*SET,DELT,0.01
VSEL,S, , , 4,5,1
vatt,2,2,2 !钢的属性
VSEL,S, , , 1
vatt,1,1,3 !线圈属性
VSEL,S,, ,2,3,1
VSEL,A, , ,6,7,1
vatt,3,1,1 !空气属性
网格划分部分
allsel
eplot
NSEL,EXT
D,ALL,AZ,0
FINISH
/SOL
ANTYP,HARM !
HARFR,50000
ASEL,S, , , 5
NSLA,S,1
d,all,volt,0
ASEL,S, , , 9
nsla,s,1
cp,1,volt,all
*get,n1,node, ,num,max
F,n1,AMPS,14000
allsel
PHYSICS,WRITE,EMAG ! 写入emag
FINISH
/PREP7
LSCLEAR,ALL !
ET,1,0 ! 将1号材料转换成0
ET,2,SOLID90 ! 将2号材料转换为SOLID90
ET,3,0 ! 将3号单元转换成0
MP,KXX,2,30 !!!!!传导率各向同性Thermal conductivities (also KYY, KZZ
!MP,KYY,2,30
!MP,KZZ,2,30
MP,C,2,500 !
MP,DENS,2,7900 !
FINISH
/SOLU
ANTYPE,TRANS !
TOFFST,273
TUNIF,100 !
CNVTOL,HEAT,1 !
KBC,1 !
TRNOPT,FULL
AUTOS,ON ! 自动时间-
DELTIM,5E-3,5E-4,DELT,ON !
OUTRES,BASIC,ALL !
PHYSICS,WRITE,THERMAL ! 写入THERMAL
FINISH
*do,i,1,ftime/tinc ! solution *do loop
*SET,time,time+tinc ! increment time
physics,read,emag ! read emag physics file
/solu
*if,i,eq,1,then
tunif,100 ! initial temperature
*else
ldread,temp,last,,,,,rth ! read thermal analysis temperatures
*endif
solve ! solve harmonic analysis
finish
physics,read,thermal ! read thermal physics file
/assign,esav,therm,esav ! redirect files for use in thermal res
/assign,emat,therm,emat
/solu
parsav,scalar,parameter,sav !save parameters before multiframe rest
*if,i,gt,1,then
antype,trans,rest ! thermal restart
*endif
parres,new,parameter,sav !restore parameters after multiframe re
time,time ! time at end of thermal run
esel,s,mat,,2 ! select billet region
ldread,hgen,,,,2,,rst ! apply coupled joule heating load from
esel,all
solve
finish
/assign,esav ! reassign files to default
/assign,emat
*enddo ! end of solution looping
|