磁场-热耦合问题

zhs1984

在ANSYS10.0工程电磁分析技术与实例详解中看到一个磁场热耦合的例子，按例操作却报错，无法获得结果文件，现附上命令流，请高手指点！

zhs1984

/prep7
ET,1,13,,,1             !定义PLANE13单元,轴对称,AZ自由度
ET,2,13,,,1
ET,3,151,,,1,1,1        !定义SURF151单元,热辐射
R,3,0                   !定义SURF151单元实常数
row=0.015               !部件的外径
ric=0.0175              !线圈内径
roc=0.02                !线圈外径
ro=0.05                 !模型外径
t=0.001                 !模型厚度
freq=150000
pi=4*atan(1)
cond=3.92e6             !最大电导率
muzero=4e-7*pi          !真空磁导率
mur=200                 !最大相对磁导率
skind=sqrt(1/(pi*freq*cond*muzero*mur))  ! 趋肤深度

ftime=3                 !结束时间
tinc=0.05               !谐波分析时间增量
time=0                  !起始时间
delt=0.01               !最大时间步
emunit,mks              !设置磁场单位
mp,murx,1,1             !空气相对导磁率
mp,murx,3,1             !线圈相对导磁率
mptemp,1,25.5,160,291.5,477.6,635,698     !钢的相对磁导率与温度的关系
mptemp,7,709,720.3,742,761,1000
mpdata,murx,2,1,200,190,182,161,135,104   !钢的相对磁导率与温度的关系
mpdata,murx,2,7,84,35,17,1,1
mptemp
mptemp,1,0,125,250,375,500,625,           !钢的电阻率与温度的关系
mptemp,7,750,875,1000                     !钢的电阻率与温度的关系
mpdata,rsvx,2,1,1.84e-7,2.72e-7,3.84e-7,5.12e-7,6.56e-7,8.24e-7  
mpdata,rsvx,2,7,1.032e-6,1.152e-6,1.2e-6
rectng,0,row,0,t                           !钢条
rectng,row,ric,0,t                         !空气隙
rectng,ric,roc,0,t                         !线圈
rectng,roc,ro,0,t                          !外空气区
aglue,all
numcmp,area
asel,s,,,1                                  !设置钢条属性
aatt,2,,1
asel,s,,,3                                 !设置线圈属性
aatt,3,,2
asel,s,,,2,4,2                             !设置空气区属性
aatt,1,,2
allsel
mshape,0,2d
mshk,1
amesh,1
lsel,s,loc,y,o
lsel,a,loc,y,t
lsel,u,loc,x,row
lesize,all,0.001
lsel,all
amesh,all
allsel
n                                        !为SURF151单元建立空间节点
*get,nmax,node,,num,max
lsel,s,loc,row
type,3
real,3
mat,2
lmesh,all
*get,emax,elem,,num,max
emodif,emax,3,nmax                       !把单元修改成用于辐射的空间节点
et,3,0                                   !把3号单元转换成0号单元
nsel,s,loc,x
d,all,az,0                               !施加磁力线平等边界条件
nsel,all
esel,s,mat,,3                            !选择线圈-按材料3号
bfe,all,js,,,,15e6                       !对线圈施加电流密度
esel,all
antyp,harm
harfrq,150000
physics,write,emag                       !写入emag物理环境
!单元转换,转入热分析
lsclear,all              !清除所有的边界条件和选项
et,1,55,,,1              !PLANE55热单元,轴对称
et,2,0                   !线圈和空气区为0号单元
et,3,151,,,1,1,1         !用于辐射的SURF151单元
keyopt,3,9,1            
r,3,1,5.67e-8            !定义factor,Stefan-Boltzmann常数
mptemp
mptemp,1,0,730,930,1000                     !定义温度与传导率的关系
mpdata,kxx,2,1,60.64,29.5,28,28
mptemp                                      !定义温度与焓的关系
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,0.68                              !辐射率
!设置热求解选项
antype,trans                 !定义瞬态分析
toffst,273
tunif,100                    !初始归一化温度
d,nmax,temp,25               !环境温度
cnvtol,heat,1                !收敛公差
kbc,1
trnopt,full
autos,on                     !自动时间-步
deltim,1e-5,1e-6,delt,on     !时间步控制
outres,basic,all             !保存所有载荷步信息
!写入物理环境thermal
physics,write,thermal
!用命令流实现循环求解
*do,i,1,ftime/tinc
time=time+tinc
physics,read,emag            !读入物理文件
/solu
*if,i,eq,1,then
tunif,100                   !初始化温度
*else
ldread,temp,last,,,,,rth    !读入热分析的温度
*endif
solve                        !求谐波分析
finish
physics,read,thermal         !读入物理文件

！从这里开始出错

/assign,esav,therm,esav      !redirect files for use in thermal restart
/assign,emat,therm,emat
/solu
*if,i,gt,i,then
antype,trans,rest           !热重启动
*endif
time,time                    !热分析结束时间
ldread,hgen,,,,2,,rmg        !从emag分析中读入焦耳热载荷
esel,all
solve
finish

/assign,esav                  !
/assgin,emat
*enddo                        !循环结束
finish
save
/post26
/show
nsol,2,1,temp,,tempcl         !储存钢条中心线的温度
nsol,3,2,temp,,tempsurf       !储存钢条外径的温度
plvar,2,3                     !画出温度随时间上升的曲线
prvar,2,3
finish

yinflying

楼主出错信息是什么啊？
我的也出过错 不过保存一下文件 就不出错了

zhs1984

是保存什么文件呀？ 3# yinflying

zoe122208

我做了两遍，如楼主所说，晕，不会。

xulang1013

呵呵 尤其是多场耦合法求解时在指定把单元分成 磁场1  热场2 时 会出错 说啥AZ平行啥的  但保存下 就不会弹错  虽然不弹 但是我想这其实已经出错了  只是我们找不到错的原因啊 请问有知道到底是怎么回事吗

weizumin

n                                        !为SURF151单元建立空间节点
这句有点问题，好象把他建在铁区右侧中心处就能解决