本帖最后由 twt716 于 2010-5-26 13:11 编辑
我做感应加热的分析,用物理环境耦合法,需要将电磁场分析的数据导入热分析中,但电磁场求解完后没有生成.rmg文件,没法读取啊,请问高手这是什么回事?
我大部分用的GUI操作,所以log文件比较乱,下面是一部分,希望大家帮我解答一下这个问题,急!谢谢!
/prep7
/title,bearing induction heating simulatin
/com,
et,1,13,,,1 !设置物理场分析单元
et,2,53,,,1
!*
!*
ET,3,PLANE53
!*
ET,4,INFIN110
!*
KEYOPT,3,1,2
KEYOPT,3,2,0
KEYOPT,3,3,1
KEYOPT,3,4,0
KEYOPT,3,5,0
KEYOPT,3,7,0
!*
KEYOPT,4,1,0
KEYOPT,4,2,0
KEYOPT,4,3,1
!*
EMUNIT,MKS
!*
MPTEMP,,,,,,,, !设置物理材料属性
MPTEMP,1,0
MPDATA,MURX,1,,1
*CSET,1,2, 1,2,
MPCOPY, ,1,2
TBCOPY,ALL,1,2
!*
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,RSVX,2,,0.2e-6
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,MURX,3,,190
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,RSVX,3,,2e-7
!*
R,1,4e-3,800, ,1,1,
!*
rectng,0,.1,0,.15 !实体建模
rectng,.1,.115,0,.05
rectng,0.165,0.245,0,0.05
rectng,0.345,0.6,0,0.15
rectng,0.1,0.345,0,0.15
FLST,2,5,5,ORDE,2
FITEM,2,1
FITEM,2,-5
AOVLAP,P51X
/DIST,1,1.08222638492,1
/REP,FAST
/DIST,1,0.924021086472,1
/REP,FAST
/PNUM,KP,0
/PNUM,LINE,0
/PNUM,AREA,1
/PNUM,VOLU,0
/PNUM,NODE,0
/PNUM,TABN,0
/PNUM,SVAL,0
/NUMBER,0
!*
/PNUM,ELEM,0
/REPLOT
!*
CM,_Y,AREA
ASEL, , , , 7
CM,_Y1,AREA
CMSEL,S,_Y
!*
CMSEL,S,_Y1
AATT, 3, 1, 1, 0,
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!*
CM,_Y,AREA
ASEL, , , , 6
CM,_Y1,AREA
CMSEL,S,_Y
!*
CMSEL,S,_Y1
AATT, 3, 1, 1, 0,
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!*
CM,_Y,AREA
ASEL, , , , 8
CM,_Y1,AREA
CMSEL,S,_Y
!*
CMSEL,S,_Y1
AATT, 1, 1, 2, 0,
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!*
CM,_Y,AREA
ASEL, , , , 3
CM,_Y1,AREA
CMSEL,S,_Y
!*
CMSEL,S,_Y1
AATT, 2, 1, 3, 0,
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!*
CM,_Y,AREA
ASEL, , , , 4
CM,_Y1,AREA
CMSEL,S,_Y
!*
CMSEL,S,_Y1
AATT, 1, 1, 4, 0,
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!*
!!!!!!网格划分部分,省!!!!!!!!!!!!!!!!
*get,nmax,node,,num,max
*SET,n1,node(0.165,0,0) !在(0.165,0,0)处设置一节点,以耦合电流自由度
ESEL,S,MAT,,2 !耦合电流自由度
NSLE,S
FLST,4,225,1,ORDE,2
FITEM,4,1
FITEM,4,-225
CP,1,CURR,P51X
ALLSEL,ALL
FINISH
/SOL
NSEL,S,LOC,X,0.6
FLST,2,9,1,ORDE,2
FITEM,2,227
FITEM,2,-235
SF,P51X,INF !远场单元加远场标志
NSEL,S,LOC,X,0
FLST,2,16,1,ORDE,4
FITEM,2,244
FITEM,2,247
FITEM,2,269
FITEM,2,-282
D,P51X,AZ,0 !设x=0处磁势平行边界条件
ESEL,S,MAT,,2
FLST,2,64,2,ORDE,2
FITEM,2,1
FITEM,2,-64
BFE,P51X,VLTG,1,380, , !加电压载荷
esel,all
finish
/solu
antyp,harm !设置谐波分析
harfrq,50
ALLSEL,ALL
physics,write,emag !写物理环境
finish
/prep7
lsclear,all !设置热分析环境
et,1,55,,,1
et,2,0
et,3,0
et,4,0
!*
MPTEMP,,,,,,,, !设置材料属性,为简单此处均设成常数
MPTEMP,1,0
MPDATA,KXX,1,,48
MPDE,KXX,1
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,KXX,3,,48
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,3,,7000
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,ENTH,3,,1e9
/REPLOT,RESIZE
/REPLOT,RESIZE
finish
/solu
antype,trans !设置瞬态分析
toffst,273 !温度摄氏单位
tunif,100
tunif,25 !节点初始温度
d,nmax,temp,25
cnvtol,heat,1
kbc,1 !阶跃方式
trnopt,full !full方法
autots,on
deltim,1e-5,1e-6,.01,on
outres,basic,all
physics,write,thermal
finish
/REPLOT,RESIZE
physics,read,emag
/solu
solve !?????????此处未生成.rmg文件
finish
physics,read,thermal
/solu
antype,trans,rest
esel,s,mat,,3
ldread,hgen,,,,3,,rmg ! 上面没有生成.rmg文件,导致这里无法读取数据 |