finish
/clear
/NOPR
/title,the transient anslysis of the reactor
/filname,reactor_tran,1
keyw,magedg,1
/prep7
et,1,53
ET,2,117
KEYOPT,2,1,1
Mp,rsvx,2,1.673e-8
Mp,murx,2,1
Mp,rsvx,3,2.83e-8
Mp,murx,3,1
Mp,murx,4,1 !空气
turn=5
lastpoint1=turn*4
lastpoint=turn*4-2
rin=3
rout=6
m=3
n=3
*do,i,1,lastpoint,4
k,i,rin+m*(i-1)
k,i+1,,rin+m*i
k,i+2,-rin-(m*(i+1))
k,i+3,,-rin-(m*(i+2))
*ENDDO
csys,1
*do,i,1,lastpoint+1
l,i,i+1
*enddo
lsel,s,line,,1,lastpoint+3
csys,0
firpoint2=turn*4+1
lastpoint2=firpoint2+lastpoint
*do,i,firpoint2,lastpoint2,4
k,i,rout+m*(i-firpoint2)
k,i+1,,rout+m*(i-firpoint2+1)
k,i+2,-rout-(m*(i-firpoint2+2))
k,i+3,,-rout-(m*(i-firpoint2+3))
*ENDDO
csys,1
*do,i,firpoint2,lastpoint2
l,i,i+1
*enddo
*get,kmax,kp,,num,max
*get,lmax,kp,kmax,loc,y
csys,0
lstr,1,firpoint2
lstr,lastpoint1,lastpoint2+1
type,1
lsel,all
al,all
type,2
mat,2
vext,all,,,,,-60
!-------------建切割面,用面切割体
*get,lleft,kp,37,loc,x
*get,lright,kp,39,loc,x
*get,lup,kp,38,loc,y
*get,ldown,kp,40,loc,y
wprot,0,90
rectng,lleft,lright,0,-60
*get,areanum1,area,,num,max
wprot,0,0,90
rectng,lup,ldown,0,-60
*get,areanum2,area,,num,max
vsel,s,,,1
vsba,1,areanum1
vsel,s,,,2,11
cm,volucoil,volum
vsba,volucoil,areanum2
!----------空气建模
wpcsys,-1
allsel,all
type,2
mat,4
block,lmax*3,-lmax*3,lmax*3,-lmax*3,-60*4,60*3
!----------空气建模完毕
vovlap,all !overlap线圈体和空气体
!-------------------对线圈的网格划分
vsel,s,loc,z,0,-60
esize,1,0
mshkey,1
asel,s,loc,z,0
amesh,all
lsel,s,loc,z,-30
lsel,r,length,,60
lesize,all,10
vsweep,all
!-------------------对线圈的网格划分完毕
这是我的命令,做到线圈的网格划分完毕这里,下面再进行空气体的网格划分就不知道该怎么做了,尝试了各种方法都不行,有没有牛人帮我看看,这个空气体该怎么划分? |