磁热耦合

wenfang

我现在正在模拟磁热顺序耦合过程，先做了一个书上的例题，可能是书上写得简单，我在模拟过程中遇见了问题。
命令流如下：

/prep7

et,1,13,,,1
et,2,13,,,1
et,3,151,,,1,1,1                   ！选择磁场单元               
r,3,0                              ！SURF151单元实常数

*SET,row,0.015                     ！钢盘直径
*SET,ric,0.0175                    ！线圈靠钢盘的内径
*SET,roc,0.0200                    ！线圈外径
*SET,ro,0.05                        空气的最远距离
*SET,t,0.001                        模型厚度
*SET,freq,150000                    频率
*SET,pi,4*atan(1)   
*SET,cond,0.392e7   
*SET,muzero,4e-7*pi                  真空磁导率
*SET,mur,200            
*SET,skind,sqrt(1/(pi*freq*cond*muzero*mur))   趋肤深度
*SET,ftime,3                                     结束时间
*SET,tinc,0.05                                  时间增量
*SET,time,0                                    开始时间
*SET,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,0.184e-6,0.272e-6,0.384e-6,0.512e-6,0.656e-6,0.824e

-6   
mpdata,rsvx,2,7,1.032e-6,1.152e-6,1.2e-6


rectng,0,row,0,t
rectng,row,ric,o,t  
rectng,ric,roc,0,t  
rectng,roc,ro,o,t  
aglue,all                           ！创建2D模型
numcmp,area                         ！压缩面编号

asel,s,area,,1                     
aatt,2,1,1  
asel,s,area,,3  
aatt,3,1,2  
asel,s,area,,2,4,2  
aatt,1,1,2  
asel,all                         ！分别对钢、线圈和空气赋特性
  
mshape,0,2d
mshk,1  
amesh,1
lsel,s,loc,y,0  
lsel,a,loc,y,t  
lsel,u,loc,x,row/2  
lesize,all,0.001
lsel,all
amesh,all                          ！网格划分

n   ！为SURF151单元建立空间节点
*get,nmax,node,,num,max
lsel,s,loc,x,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   
bfe,all,js,,,,15e6  
esel,all            施加的边界条件和载荷


antyp,harm  
harfrq,150000   ！定义谐波分析
physics,write,emag  ！写入emag物理文件


！下面转入热分析：
lsclear,all ！清除所有边界条件和选项

et,1,55,,,1
et,2,0  
et,3,151,,,1,1,1

keyopt,3,9,1
r,3,1,5.67e-8   
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
autots,on   
deltim,1e-5,1e-6,delt,on
outres,basic,all
physics,write,thermal   ！写入thermal物理文件


做do循环
*do,i,1,ftime/tinc  
*SET,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
/assign,emat,therm,emat

/solu
*if,i,gt,1,then
antyp,trans,rest
*endif  
time,time   
esel,s,mat,,2   
ldread,hgen,,,,2,,rmg          ！读入emag分析中的焦耳热载荷
esel,all
solve   
finish
/assign,esav
/assign,emat

*enddo
finish
save
/post26
/show
nsol,2,1,tem,,tempcl
nsol,3,2,tem,,tempself
plvar,2,3
prvar,2,3
finish



但是在加载过程中
*** WARNING ***                       
Specified DOF constraint AZ at unused node 193.  
*** ERROR ***  
Shape testing revealed that 1 of the 3 new or modified elements violate
shape error limits.  To review test results, please see the output      
file or issue the CHECK command.                                       
*** ERROR ***                           CP =      32.837   TIME= 16:35:02
Element 142 has an undefined node number 0.                             
*** ERROR ***                           CP =      33.017   TIME= 16:35:09
Element 143 has an undefined node number 0.                             
*** ERROR ***                           CP =      33.198   TIME= 16:35:11
The model data was checked and there were errors found.                 
  Please check output or errors file ( C:\Documents and                  
Settings\wenfang\file.err ) for messages.                              
*** ERROR ***                           CP =      33.388   TIME= 16:35:13
No results file ( file.rth ) available.                                 
  The LDREAD command is ignored.                                         
*** WARNING ***                         CP =      33.648   TIME= 16:36:02
Specified DOF constraint AZ at unused node 193.                        
*** WARNING ***                         CP =      33.698   TIME= 16:36:02
Material property MURX of material 2 of element 1 is evaluated at a     
temperature of 0, which is below the supplied temperature range.        
Temperature range checking terminates.                                 
*** ERROR ***                           CP =      33.949   TIME= 16:36:02
Shape testing revealed that 1 of the 3 new or modified elements violate
shape error limits.  To review test results, please see the output      
file or issue the CHECK command.                                       
*** ERROR ***                           CP =      34.530   TIME= 16:36:09
File therm.esav not found.  Previous solution (required for restarting)
cannot be recovered.                                                   
*** ERROR ***                           CP =      34.750   TIME= 16:36:13
Error in restart, RDB file does not exist.                              
  The ANTYPE command is ignored.                                         
*** WARNING ***                         CP =      34.970   TIME= 16:36:38
File file.rth ( needed for appending ) does not exist.                  
*** WARNING ***                         CP =      34.970   TIME= 16:36:38
File file.rth ( needed for appending ) does not exist.                  
*** ERROR ***                           CP =      34.970   TIME= 16:36:38
File therm.esav not found.  Previous solution (required for restarting)
cannot be recovered.                                                   
*** ERROR ***                           CP =      35.181   TIME= 16:36:50
No results file ( file.rth ) available.                                 
  The LDREAD command is ignored.                                         
*** WARNING ***                         CP =      35.361   TIME= 16:36:54
Specified DOF constraint AZ at unused node 193.                        
*** ERROR ***                           CP =      35.631   TIME= 16:36:55
Shape testing revealed that 1 of the 3 new or modified elements violate
shape error limits.  To review test results, please see the output      
file or issue the CHECK command.                                       
*** ERROR ***                           CP =      36.202   TIME= 16:37:00
File therm.esav not found.  Previous solution (required for restarting)
cannot be recovered.                                                   
*** ERROR ***                           CP =      36.462   TIME= 16:37:04
Error in restart, RDB file does not exist.                              
  The ANTYPE command is ignored.                                         
*** WARNING ***                         CP =      36.693   TIME= 16:37:39
File file.rth ( needed for appending ) does not exist.                  
*** WARNING ***                         CP =      36.693   TIME= 16:37:39
File file.rth ( needed for appending ) does not exist.                  
*** ERROR ***                           CP =      36.693   TIME= 16:37:39
File therm.esav not found.  Previous solution (required for restarting)
cannot be recovered.                                                   
*** ERROR ***                           CP =      36.953   TIME= 16:37:55
No results file ( file.rth ) available.                                 
  The LDREAD command is ignored.                                         
*** WARNING ***                         CP =      37.163   TIME= 16:37:59
Specified DOF constraint AZ at unused node 193.                        
*** ERROR ***                           CP =      37.474   TIME= 16:37:59
Shape testing revealed that 1 of the 3 new or modified elements violate
shape error limits.  To review test results, please see the output      
file or issue the CHECK command.                                       
*** ERROR ***                           CP =      38.055   TIME= 16:38:01
File therm.esav not found.  Previous solution (required for restarting)
cannot be recovered.                                                   
*** ERROR ***                           CP =      38.325   TIME= 16:38:05
Error in restart, RDB file does not exist.                              
  The ANTYPE command is ignored.                                         
*** WARNING ***                         CP =      38.505   TIME= 16:38:06
File file.rth ( needed for appending ) does not exist.                  
*** WARNING ***                         CP =      38.505   TIME= 16:38:06
File file.rth ( needed for appending ) does not exist.                  
*** ERROR ***                           CP =      38.505   TIME= 16:38:06
File therm.esav not found.  Previous solution (required for restarting)
cannot be recovered.     



请您帮我看看哪儿出了问题，谢谢！

twt716

1# wenfang

LZ,你这个问题解决没有？

simwesy

楼主你好强大啊，呵呵。