磁热耦合分析的实例及问题
运用MFS代码分析磁热耦合问题，谐波电磁分析+瞬态热分析，出现如下提示：*** WARNING *** SUPPRESSED MESSAGE CP = 131.844 TIME= 21:33:14
Element 92297 has element ( HGEN ) input representing two different
values ( INFINITY, INFINITY) at the same node number 126343.
请问这是什么问题啊？？？
代码附上：
/title,MFS
/filname,MFS
/vie,1,2,-1,3
/pnu,mat,1
/num,1
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!电磁模型
!***1-电磁参数
rho_copper=1.68e-8 !铜的电阻率,欧米
rho_steel=1.18e-7 !铁在高温下的电阻率，欧米
freq=250000 !电流频率HZ
I=500 !电流大小A
/prep7
emunit,mks
!***2-定义单元类型
et,1,solid117,0 !线圈,抑制涡流
et,2,solid117,1 !铁样品球
et,3,solid117,0 !空气
!***3-定义材料属性
!线圈
mp,murx,1,1
mp,rsvx,1,rho_copper
!空气
mp,murx,3,1
!!!!!!铁的电磁及热属性
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
!***4-几何建模
! *** 线圈
nturn=7 !线圈匝数
r1_coil=0.010 ! inner radius, coils
r2_coil=0.016 ! outer radius, coils
h_coil=0.035 ! height, coils
sc_area=(r2_coil-r1_coil)*h_coil !线圈截面积
Jdens=nturn*I/sc_area !电流密度
cyl4,,,r1_coil,0,r2_coil,360,h_coil
local,11,1,,,,
! *** 铁球
r_ball=0.005
h_ball=h_coil+0.01
wpoffs,,,h_ball
sphere,,r_ball,,360
wpcsys,-1,0
csys,0
! *** 空气
sp=0.04
!*空气几何建模
allsel
*get,xmin,kp,,mnloc,x
*get,xmax,kp,,mxloc,x
*get,ymin,kp,,mnloc,y
*get,ymax,kp,,mxloc,y
*get,zmin,kp,,mnloc,z
*get,zmax,kp,,mxloc,z
rmax=sqrt(xmax**2+ymax**2)
wpoffs,,,zmin-sp
cyl4,,,rmax+sp,,,360,zmax+2*sp
vsbv,3,1,,dele,keep
numcmp,all
vsbv,3,2,,dele,keep
numcmp,all
wpcsys,-1,0
csys,0
!***5-划分网格
esz1=0.002 ! element size, components
esz2=5*esz1 ! element size, surrounding domain
!lsel,s,loc,z,0
!lsel,a,loc,z,h_coil
!lesize,all,,,8
vsel,s,,,1
vatt,1,1,1,11 !柱坐标
ESIZE,esz1
MSHAPE,0,3d
vsweep,1
ESIZE,esz1/2
vsel,s,,,2
vatt,2,2,2
MSHAPE,1,3d
MSHKEY,0
vmesh,2
vsel,s,,,3
vatt,3,3,3
ESIZE,esz2
MSHAPE,1,3d
MSHKEY,0
vmesh,3
allsel
!***7-边界条件及激励
allsel
nsel,s,ext
d,all,az,0 !外边界磁力线平行
nsel,none
vsel,s,,,1
ALLSEL,BELOW,VOLU
bfe,all,js,1,,Jdens !线圈电流密度,峰值
ALLSEL
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!热模型
et,4,solid90
VSEL,S,,,2
vgen,2,2,,,,,,,1 !复制体，但不复制其单元及节点
ESIZE,esz1/2
vsel,s,,,4
vatt,2,,4 !铁样品
MSHAPE,1,3d
MSHKEY,0
vmesh,4
vsel,s,,,4
NSLV,S,1
NSEL,R,EXT
SF,ALL,CONV,35,25
esel,s,mat,,2 ! 选择样品材料
bfe,all,fvin,,1 ! 定义体积载荷界面
finish
/solu
mfan,on ! Activate ANSYS Multi-field solver analysis
mfel,1,1,2,3 ! Field #1 ET;s, Emag
mfel,2,4 ! Field #2 ET's, Thermal
mfor,1,2 ! Field solution order
mfti,3 ! Final time
mfdt,0.05 ! Stagger time increment
mfco,all,1e-3 ! Convergence criteria
antyp,harm ! 电磁场分析选项
harfrq,freq
outres,all,all
tunif,25
mfcm,1, ! Write Emag analysis options
mfclear,solu ! Clear analysis options
antype,trans ! 热分析选项
toffst,273
tunif,25 ! initial uniform temperature
kbc,1 ! step loads
trnopt,full
nropt,full,,on
autos,on
deltim,.01,.005,.01,on
mfcm,2, ! Write Thermal analysis options
mfvo,1,1,hgen,2 ! 载荷传递
mfvo,1,2,temp,1 ! 载荷传递
solve
save
finish 类似的警告有3000多个，计算一段时间就挂掉了。
页:
[1]