下面是我的命令流,有几个问题:1、本来铁芯与样板都是q235钢,但因为是谐态的,材料不能为非线性,在网上找的将其材料改为磁导率4000 2、样板的单元选项是AZ+volt,在对成面进行volt耦合,不知道这样对不 3、网格的问题,sweep划分时与free差别很大,为嘛、
最重要的是想知道这样命令流有没什么问题,能力有限,感觉没问题,但是路径出来的结果不对
/title, 3-D Harmonic Response, 1/2 symmetry
/vie,1,2,-1,3
/pnu,mat,1
/num,1
! *** Model parameters
a1_core=0.02 ! core cross-section width
a2_core=0.04 ! core cross-section depth
w1_core=0.08 ! medium core width
w2_core=a1_core*2+w1_core ! overall core width
h_core=0.08 ! overall core height
d_core=0
thk_coil1=0.01 !left right thick
thk_coil2=0.015 !UP thick
h_coil1=0.04 ! height,left right coils
h_coil2=0.058 ! height,UP coils
d_plate=0.001 !plate
l_plate=0.5
h_plate=0.01
w_plate=0.2
l_flaw=0.0005 !length
h_flaw=0.002 !depth
w_flaw=0.005
d_dmn=0.025 ! depth of surrounding domain
esz1=a1_core/2 ! element size, components
esz2=2*esz1 ! element size, surrounding domain
esz3=thk_coil2/3
esz4=h_flaw/3
esz5=h_plate/2
rs_q235=0.12E-06
rs_coil=1.724e-8
! *** Primary coil
I=10 ! current (A)
N1=100 ! # of turns
R1=0.0668 ! DC resistance (ohms)
care_left=thk_coil1*h_coil1 ! cross-sectional area
volu_left=((thk_coil1*2+a1_core)*(thk_coil1*2+a2_core)-a1_core*a2_core)*h_coil1 ! volume left
js1=n1*I/care_left
N2=200 ! # of turns
R2=0.1336 ! DC resistance (ohms)
care_up=thk_coil2*h_coil2 ! cross-sectional area
volu_up=((thk_coil2*2+a1_core)*(thk_coil2*2+a2_core)-a1_core*a2_core)*h_coil2 ! volume up
js2=n2*I/care_up
frqncy=50 ! operating frequency (Hz)
symm=2 ! symmetry factor
/nopr
! *** Geometry
/PREP7
wpcs,-1,0
bloc,-(a1_core*2+w1_core)/2,(a1_core*2+w1_core)/2,h_core-a1_core,h_core,-a2_core/2,0 !core
bloc,-(a1_core*2+w1_core)/2,-(a1_core*2+w1_core)/2+a1_core,0,h_core-a1_core,-a2_core/2,0
bloc,(a1_core*2+w1_core)/2-a1_core,(a1_core*2+w1_core)/2,0,h_core-a1_core,-a2_core/2,0
!up coil
wpcs,-1,0 ! up coil ESYS1
wpof,,h_core-a1_core/2
wpro,,,90
cswp,13,1
csys !up coil
wpcs,-1,13
bloc,0,a2_core/2,-(a1_core/2+thk_coil2),-a1_core/2,-h_coil2/2,h_coil2/2,
bloc,a2_core/2,a2_core/2+thk_coil2,-a1_core/2,a1_core/2,-h_coil2/2,h_coil2/2,
bloc,0,a2_core/2,a1_core/2,(a1_core/2+thk_coil2),-h_coil2/2,h_coil2/2,
wpcs,-1,13 ! up coil ESYS2
wpof,a2_core/2,-a1_core/2
cswp,20,1
cylind,0,thk_coil2,-h_coil2/2,h_coil2/2,-90,0,
wpcs,-1,13 ! up coil ESYS3
wpof,a2_core/2,a1_core/2
cswp,21,1
cylind,0,thk_coil2,-h_coil2/2,h_coil2/2,0,90,
!left coil
wpcs,-1,0 ! left coil ESYS1
wpof,-(w1_core+a1_core)/2
wpro,,-90,
cswp,11,1
csys !left coil
wpcs,-1,11
bloc,-a1_core/2,-(a1_core/2+thk_coil1),0,a2_core/2,0,h_coil1,
bloc,-a1_core/2,a1_core/2,a2_core/2,a2_core/2+thk_coil1,0,h_coil1,
bloc,a1_core/2,(a1_core/2+thk_coil1),0,a2_core/2,0,h_coil1,
wpcs,-1,11 ! left coil ESYS2
wpof,-a1_core/2,a2_core/2
cswp,22,1
cylind,0,thk_coil1,0,h_coil1,90,180,
wpcs,-1,11 ! left coil ESYS3
wpof,a1_core/2,a2_core/2
cswp,23,1
cylind,0,thk_coil1,0,h_coil1,0,90,
!right coil
wpcs,-1,0 ! right coil ESYS1
wpof,(w1_core+a1_core)/2
wpro,,90,
cswp,12,1
csys !right coil
wpcs,-1,12
bloc,-a1_core/2,-(a1_core/2+thk_coil1),-a2_core/2,0,-h_coil1,0,
bloc,-a1_core/2,a1_core/2,-(a2_core/2+thk_coil1),-a2_core/2,-h_coil1,0,
bloc,a1_core/2,(a1_core/2+thk_coil1),-a2_core/2,0,-h_coil1,0,
wpcs,-1,12 ! right coil ESYS2
wpof,-a1_core/2,-a2_core/2
cswp,24,1
cylind,0,thk_coil1,-h_coil1,0,-90,-180,
wpcs,-1,12 ! right coil ESYS3
wpof,a1_core/2,-a2_core/2
cswp,25,1
cylind,0,thk_coil1,-h_coil1,0,-90,0,
vglue,all
vsel,none ! plate
wpcs,-1,0
bloc,-l_plate/2,l_plate/2,-h_plate-d_plate,-d_plate,-w_plate/2,0
cm,plate_1,volu
vsel,none !flaw
bloc,-l_flaw/2,l_flaw/2,-h_flaw-d_plate,-d_plate,-w_flaw/2,0
cm,flaw_v,volu
vatt,1,1,1
cmse,s,plate_1
cmse,a,flaw_v
vsbv,plate_1,flaw_v,,dele,keep
cmse,u,flaw_v
cm,plate_v,volu
vatt,4,4,4
alls
csys
cm,keep_v,volu
*get,xmin,kp,,mnloc,x
*get,xmax,kp,,mxloc,x
*get,ymax,kp,,mxloc,y
*get,ymin,kp,,mnloc,y
*get,zmin,kp,,mnloc,z
wpcs,-1,0 ! surrounding domain
bloc,xmin-d_dmn,xmax+d_dmn,ymin-d_dmn,ymax+d_dmn,zmin-d_dmn,0
allsel,all
vovlap,all
numcmp,all
vsel,s,volu,,1
vatt,4,1,4
vsel,s,volu,,2
vsel,a,volu,,11,18
vatt,3,1,3
vsel,s,volu,,5
vatt,3,1,3,20
vsel,s,volu,,6
vatt,3,1,3,21
vsel,s,volu,,7
vatt,3,1,3,22
vsel,s,volu,,8
vatt,3,1,3,23
vsel,s,volu,,9
vatt,3,1,3,24
vsel,s,volu,,10
vatt,3,1,3,25
vsel,s,volu,,19,21
vatt,2,1,2
vsel,s,volu,,3,4
vatt,1,1,1
alls
! *** FE model
et,1,236 ! air
mp,murx,1,1
et,2,236 ! core (laminated, non-conducting)
mp,murx,2,4000
et,3,236 ! left primary coil
mp,murx,3,1
mp,rsvx,3,rs_coil
et,4,236,1 ! plate
mp,murx,4,4000
mp,rsvx,4,rs_q235
! *** Mesh
numm,kp,1e-8,1e-8
msha,0
esiz,esz3
vsel,s,mat,,3
vswe,all
esiz,esz3
vsel,s,mat,,2
vswe,all
esiz,esz4
vsel,s,volu,,3
vswe,all
msha,1
esiz,esz5
vsel,s,mat,,4
vmes,all
vsel,s,mat,,1
msha,1
esiz,esz2
vmes,all
alls
! *** Boundary conditions and loads
asel,s,ext ! flux parallel exterior
csys
da,all,az
vsel,s,mat,,4
aslv,s
asel,r,loc,z,0
nsla,s
d,all,volt,0
nsel,all
vsel,s,volu,,2
alls,belo,volu
bfe,all,js,,0,0,-js2
vsel,s,volu,,5,6
alls,belo,volu
bfe,all,js,,0,js2,0
vsel,s,volu,,11
alls,belo,volu
bfe,all,js,,0,js2,0
vsel,s,volu,,12
alls,belo,volu
bfe,all,js,,0,0,js2
vsel,s,volu,,13
alls,belo,volu
bfe,all,js,,0,0,js1
vsel,s,volu,,7,8
alls,belo,volu
bfe,all,js,,0,js1,0
vsel,s,volu,,14
alls,belo,volu
bfe,all,js,,-js1,0,0
vsel,s,volu,,15
alls,belo,volu
bfe,all,js,,0,0,-js1
vsel,s,volu,,16
alls,belo,volu
bfe,all,js,,0,0,-js1
vsel,s,volu,,9,10
alls,belo,volu
bfe,all,js,,0,js1,0
vsel,s,volu,,17
alls,belo,volu
bfe,all,js,,js1,0,0
vsel,s,volu,,18
alls,belo,volu
bfe,all,js,,0,0,js1
csys
alls
eplo
fini
! *** Solution
/solu
antype,harmonic
harf,frqncy
solve
fini
! *** Post-processing
/post1
set,,,,0 ! Real solution set
/com,
vsel,s,mat,,2,4
alls,belo,volu
plve,jt,,,,vect,,on
plve,b,,,,vect,,on
plnsol,emf
alls
alls
path,path1,2,10,40
ppath,1,,-0.01,-0.001,0
ppath,2,,0.01,-0.001,0
pdef,By,B,y,avg
pdef,Bx,B,x,avg
plpath,BY,BX,
prpath,BY,BX,
fini
希望路过的能帮帮忙,谢谢啦 |