电磁场与热场耦合时候，读取涡流生成焦耳热，两个单元类型必须要对应吗？
电磁场与热场耦合时候，读取涡流生成焦耳热，两个单元类型必须要对应吗？ 看看下边的命令流：/batch,list
/filenam, induc
/prep7
shpp,off
/title, induction heating of a solid cylinder billet
/com,
et,1,13,,,1 ! PLANE13, axisymmetric, AZ dof
et,2,13,,,1
et,3,151,,,1,1,1 ! SURF151, thermal, radiation
r,3,0 ! Real constant set for SURF151
row=.015 ! outer radius of workpiece
ric=.0175 ! inner radius of coil
roc=.0200 ! outer radius of coil
ro=.05 ! outer radius of model
t=.001 ! model thickness
freq=150000 ! frequency (Hz.)
pi=4*atan(1) ! pi
cond=.392e7 ! maximum conductivity
muzero=4e-7*pi ! free-space permeability
mur=200 ! maximum relative permeability
skind=sqrt(1/(pi*freq*cond*muzero*mur)) ! skin depth
ftime=3 ! final time
tinc=.05 ! time increment for harmonic analysis
time=0 ! initialize time
delt=.01 ! maximum delta time step
emunit,mks ! set magnetic units
mp,murx,1,1 ! air relative permeability
mp,murx,3,1 ! coil relative permeability
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
rectng,0,row,0,t ! billet
rectng,row,ric,0,t ! air-gap
rectng,ric,roc,0,t ! coil
rectng,roc,ro,0,t ! outer air
aglue,all
numcmp,area
ksel,s,loc,x,row ! select keypoints at outer radius of workpiece
kesize,all,skind/2 ! set meshing size to 1/2 skin depth
ksel,s,loc,x,0 ! select keypoints at center
kesize,all,40*skind ! set meshing size
lsel,s,loc,y,t/2 ! select vertical lines
lesize,all,,,1 ! set 1 division through thickness
lsel,all
asel,s,area,,1
aatt,2,1,1 ! set attributes for billet region
asel,s,area,,3
aatt,3,1,2 ! set attributes for coil region
asel,s,area,,2,4,2
aatt,1,1,2 ! set attributes for air region
asel,all
mshape,0,2d
mshk,1
amesh,1 ! mesh billet area
lsel,s,loc,y,0
lsel,a,loc,y,t
lsel,u,loc,x,row/2
lesize,all,.001
lsel,all
amesh,all ! mesh remaining areas
n ! create space node for SURF151
*get,nmax,node,,num,max
lsel,s,loc,x,row
type,3
real,3
mat,2
lmesh,all ! mesh billet outer radius with SURF151
*get,emax,elem,,num,max
emodif,emax,3,nmax ! modify element to add space node for radiation
et,3,0 ! reset type 3 to null element
nsel,s,loc,x
d,all,az,0 ! apply flux-normal b.c.
nsel,all
esel,s,mat,,3
bfe,all,js,,,,15e6 ! apply current density to coil
esel,all
finish
/solu
antyp,harm
harfrq,150000
physics,write,emag ! write emag physics file
finish
/prep7
lsclear,all ! clear all b.c.'s and options
et,1,55,,,1 ! PLANE55 thermal element, axisymmetric
et,2,0 ! null element type for coil and air region
et,3,151,,,1,1,1 ! SURF151 element for radiation
keyopt,3,9,1
r,3,1,5.67e-8 ! form factor, Stefan-Boltzmann constant
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
finish
/solu
antype,trans
toffst,273
tunif,100 ! initial uniform temperature
d,nmax,temp,25 ! ambient temperature
cnvtol,heat,1 ! convergence tolerance
kbc,1 ! step loads
trnopt,full
autos,on ! auto time-stepping
deltim,1e-5,1e-6,delt,on ! time step control
outres,basic,all ! save all load step information
physics,write,thermal ! write thermal physics file
finish
*do,i,1,ftime/tinc ! solution *do loop
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 restart
/assign,emat,therm,emat
/solu
*if,i,gt,1,then
antype,trans,rest ! thermal restart
*endif
time,time ! time at end of thermal run
esel,s,mat,,2 ! select billet region
ldread,hgen,,,,2,,rmg ! apply coupled joule heating load from emag
esel,all
solve
finish
/assign,esav ! reassign files to default
/assign,emat
*enddo ! end of solution looping
finish
看看上边的红线部分 ，我在查询ETCHG 命令中看到：
ETCHG, Cnv
Changes element types to their corresponding types
MTT -- Magnetic to Thermal53 > 7762 > 7096 > 7097 > 70117 > 90119 > 87120 > 90
在这个当中没有说到电磁场中的 13单元 对应 热场中的55单元，这样的话我就有问题了:
1) 在划分网格的时候是先画电磁场的，那么在电磁场的划分的基础上，分配热单元，两者是否总是共容的呢？
2）电磁场分析后的涡流产生的焦耳热是在每个电磁场单元内的，在热场分析的时候作为载荷读入，那么它是否就会无条件的分配到每个热场单元中去呢？？而这种无条件是以共同的单元划分为前提（即在对几何相同网格，可以分配不同单元，在耦合时候是按照几何单元进行读取和分配的）
回复 #1 csuzt 的帖子
磁热耦合的时候，要求磁环境和热环境的单元划分是一样的，即在磁环境中定义好单元数，单元坐标系等，在接下来的热分析中不能变，即模型几何参数不能变，单元格也不能变，否则 两个物理场之间不能正确的传递数据。 查电磁和热单元的类型，如果耦合没问题，再查看node数，比如20-node就要匹配20-node的，保证单元和节点完全相同回复 #2 linlin2001450 的帖子
也就是说在不同场之间传递数据的条件是他们的几何参数是一致的，一一对应的。是吧 对，帮助里面都有提到耦合场那本书有个比较的表格，就是2d和3d情况，不同节点数的单元，电磁热结构单元都应该用哪个 after define magnetic physic enviroment, you should clear the load
PHYSICS,CLEAR after define magnetic physic enviroment, you should clear the load
PHYSICS,CLEAR 楼主的头像太有才了;P，路过学习下 perch88 发表于 2007-7-17 21:54 static/image/common/back.gif
对，帮助里面都有提到
耦合场那本书有个比较的表格，就是2d和3d情况，不同节点数的单元，电磁热结构单元都 ...
哥，那本书啊，找不到电磁线圈磁热耦合的资料啊，
页:
[1]