求助高手帮我看看为什么仿真不出磁感线

/title, 3-D Transformer Harmonic Response

/vie,1,2,-1,3

/pnu,mat,1

/num,1

pi=acos(-1)

! *** Model parameters

a_core=0.040 ! core cross-section width

w_core=0.100 ! overall core width

h_core=0.100 ! overall core height

r1_coil=0.025 ! inner radius, both coils

r2_coil=0.029 ! outer radius, both coils

p_coil=0.020 ! height, primary coils

s_coil=0.030 ! height, second coils

d_dmn=0.020 ! depth of surrounding domain

esz1=a_core/8 ! element size, components

esz2=esz1*2 ! element size, surrounding domain

mu_core=100000 ! core permeability

! *** Primary coil

Np=100 ! # of turns

Rp=10 ! DC resistance (ohms)

care_left=(r2_coil-r1_coil)*p_coil ! cross-sectional area

volu_left=pi*(r2_coil**2-r1_coil**2)*p_coil ! volume

Vp=220 ! voltage (V)

! *** Secondary coil

Ns=200 ! # of turns

Rs=Rp*(Ns/Np)**2 ! DC resistance (ohms)

care_right=(r2_coil-r1_coil)*s_coil ! cross-sectional area

volu_right=pi*(r2_coil**2-r1_coil**2)*s_coil ! volume

R=1e3 ! resistance attached to the secondary coil

frqncy=50 ! operating frequency (Hz)

symm=2 ! symmetry factor

/nopr

! *** Geometry

/PREP7

vsel,none ! core

bloc,-w_core/2,w_core/2,-h_core/2,h_core/2,-a_core/2,0

cm,scrap1_v,volu

vsel,none

bloc,a_core/2-w_core/2,a_core-w_core/2,a_core/2-h_core/2,h_core/2-a_core/2,-a_core/2,0

bloc,w_core/2-a_core,w_core/2-a_core/2,a_core/2-h_core/2,h_core/2-a_core/2,-a_core/2,0

cm,scrap2_v,volu

cmse,s,scrap1_v

cmse,a,scrap2_v

vsbv,scrap1_v,scrap2_v

cm,core_v,volu

vatt,2,2,2

wpcs,-1,0 ! left coil ESYS

wpof,,h_core/2-a_core/2

wpro,,-90

cswp,11,1

wpcs,-1,0 ! right coil ESYS

wpof,,a_core/2-h_core/2

wpro,,-90

cswp,12,1

csys

vsel,none ! left coil

wpcs,-1,11

cyli,r1_coil,r2_coil,0,-p_coil,0,180

vatt,3,3,3,11

vsel,none ! right coil

wpcs,-1,12

cyli,r1_coil,r2_coil,s_coil,0,0,180

vatt,4,4,4,12

alls

cm,keep_v,volu

wpcs,-1,0 ! surrounding domain

vsel,none

bloc,-w_core/2-d_dmn,w_core/2+d_dmn,-h_core/2-d_dmn,h_core/2+d_dmn,-a_core/2-d_dmn,0

cm,scrap_v,volu

cmse,all

vsbv,scrap_v,keep_v,,dele,keep

cmse,u,keep_v

cm,air_v,volu

vatt,1,1,1

alls

vplo

! *** FE model

et,1,236 ! air

mp,murx,1,1

et,2,236 ! core (laminated, non-conducting)

mp,murx,2,mu_core

et,3,236,2 ! left primary coil

mp,murx,3,1

r,3,care_left,Np,volu_left,0,1,0 ! left coil data

rmore,Rp,symm

et,4,236,2 ! right secondary coil

mp,murx,4,1

r,4,care_right,Ns,volu_right,0,-1,0 ! right coil data

rmore,Rs,symm

! *** Mesh

numm,kp,1e-8,1e-8

esiz,esz1

vsel,s,mat,,3,4

vmes,all

vsel,s,mat,,2

vswe,all

vsel,s,mat,,1

msha,1

esiz,esz2

vmes,all

alls

! *** Boundary conditions and loads

asel,s,ext ! flux parallel exterior

csys

asel,u,loc,y

da,all,az

vsel,s,mat,,3 ! left primary coil

alls,belo,volu

cp,1,emf,all

cp,2,volt,all

nd_p=ndnext(0)

alls

vsel,s,mat,,4 ! right secondary coil

alls,belo,volu

cp,3,emf,all

cp,4,volt,all

nd_s=ndnext(0)

alls

d,nd_p,volt,Vp

! *** Circuit

et,5,124,0 ! resistor connected to the secondary coil

r,5,R

*get,nmax,node,,num,max

n,nmax+1,0,h_coil/2

type,5

real,5

e,nd_s,nmax+1

d,nmax+1,volt,0 ! ground

csys

eplo

fini

! *** Solution

/solu

antype,harmonic

harf,frqncy

solve

fini