请教一个变压器仿真的问题

哪位高手能帮我看看为什么我的变压器仿真结果二次侧电压图形是这样的，我铁芯是用的非线性的BH曲线

指令如下：

/title, TEsT_1

/vie,1,2,-1,3

/pnu,mat,1

/num,1

pi=acos(-1)

! *** Model parameters

a_core=0.072 ! core cross-section width

w_core=0.090 ! overall core width

h_core=0.075 ! overall core height

b_core=0.030

d_core=0.004 !漏磁冲片总厚度

n_l=4 !漏磁冲片片数

Air_core=0.0005 !硅钢片间距

d_l=(d_core-(n_l-1)*air_core)/n_l !漏磁冲片每片厚度

N_core=20 !硅钢片数

S_core=(a_core-(n_core-1)*air_core)/n_core !硅钢片厚度

w_coil=0.013 ! 线圈宽度

d_coil=0.003 ! 线圈与铁芯距离

h_p_coil=0.015 ! height, primary coils

h_s_coil=0.022 ! height, second coils

d_air=0.001 ! 气隙间距

h_coil=h_p_coil+h_s_coil+d_air*4+d_core

d_dmn=0.020 ! 空气范围值

esz1=a_core/12 ! element size, components

esz2=a_core/12 ! element size, surrounding domain

mu_core=500000 ! core permeability

! *** primary coil

Np=126 ! # of turns

Rp=6 ! dC resistance (ohms)

care_p=w_coil*h_p_coil

volu_p_c=pi*((w_coil+d_air)**2-d_air**2)*h_p_coil

volu_p_b=(2*a_core+4*d_coil)*care_p

volu_p_b1=(2*b_core+4*d_air)*care_p

Vp=220 ! voltage (V)

! *** secondary coil

Ns=1260 ! # of turns

Rs=200 ! dC resistance (ohms)

care_s=w_coil*h_s_coil

volu_s_c=pi*((w_coil+d_air)**2-d_air**2)*h_s_coil

volu_s_b=(2*a_core+4*d_coil)*care_s

volu_s_b1=(2*b_core+4*d_air)*care_s

R=1e3 ! resistance attached to the secondary coil

symm=2 ! symmetry factor

! --- Excitation parameters

Vp=220*sqrt(2) ! voltage (V)

ncyc=1.0 ! number of cycles

nsubT=40 ! number of time steps per cycle

frqncy=50 ! frequency, Hz

! --- Derived parameters

T=1/frqncy
! time period

tend=ncyc*T ! solution time

nsubs=nsubT*ncyc ! total number of time steps

*dim,U,table,nsubs

*do,i,0,nsubs

U(i,1)=Vp*sin(2*pi*frqncy*(i/nsubs)*T)

U(i,0)=(i/nsubT)*T

*enddo

/axl,x,Time (s)

/axl,y,Volt (V)

*vplo,U(1,0),U(1,1)

/nopr

! *** Geometry

/pREp7

vsel,none ! core

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

cm,scrap1_v,volu

vsel,none

bloc,-w_coil-b_core/2-2*d_air,-b_core/2,-h_coil/2,h_coil/2,-S_core,0

bloc,b_core/2,b_core/2+w_coil+2*d_air,-h_coil/2,h_coil/2,-S_core,0

cm,scrap2_v,volu

cmse,s,scrap1_v

cmse,a,scrap2_v

vsbv,scrap1_v,scrap2_v

Vgen,N_core/2,all,,,,,-Air_core-S_core

vsel,none

bloc,-w_coil-b_core/2-d_air,-b_core/2-d_air,2*d_air+h_p_coil-h_coil/2,2*d_air+h_p_coil+d_l-h_coil/2,-a_core/2,0

bloc,b_core/2+d_air,b_core/2+w_coil+d_air,2*d_air+h_p_coil-h_coil/2,2*d_air+h_p_coil+d_l-h_coil/2,-a_core/2,0

Vgen,n_l,all,,,,Air_core+d_l,

vsel,all

cm,core_v,volu

vatt,2,2,2

wpcs,-1,0 ! primary coil left

wpof,-b_core/2,d_air-h_coil/2,-a_core/2-d_coil

wpro,,-90

cswp,11,1

wpcs,-1,0 ! primary coil right

wpof,b_core/2,d_air-h_coil/2,-a_core/2-d_coil

wpro,,-90

cswp,12,1

wpcs,-1,0 ! second coil left

wpof,-b_core/2,h_coil/2-d_air,-a_core/2-d_coil

wpro,,-90

cswp,13,1

wpcs,-1,0 ! second coil right

wpof,b_core/2,h_coil/2-d_air,-a_core/2-d_coil

wpro,,-90

cswp,14,1

csys

wpcs,-1,0
! second coil

vsel,none ! second coil left block

bloc,-w_coil-b_core/2-d_air,-b_core/2-d_air,h_coil/2-d_air-h_s_coil,h_coil/2-d_air,0,-a_core/2-d_coil

vatt,4,6,4,0

vsel,none
! second coil right block

bloc,w_coil+b_core/2+d_air,b_core/2+d_air,h_coil/2-d_air-h_s_coil,h_coil/2-d_air,0,-a_core/2-d_coil

vatt,4,8,4,0

vsel,none
! second coil behind block

bloc,-w_coil-d_air*2,w_coil+d_air*2,h_coil/2-d_air-h_s_coil,h_coil/2-d_air,-a_core/2-d_coil-d_air-w_coil,-a_core/2-d_coil-d_air

vatt,4,7,4,0

wpcs,-1,13
! second coil left cylinder

vsel,none

cyli,d_air,w_coil+d_air,-h_s_coil,0,90,180

vatt,4,9,4,13

wpcs,-1,14
! second coil right cylinder

vsel,none

cyli,d_air,w_coil+d_air,-h_s_coil,0,0,90

vatt,4,10,4,14

wpcs,-1,0 ! primary coil

vsel,none ! primary coil left block

bloc,-w_coil-b_core/2-d_air,-b_core/2-d_air,-h_coil/2+d_air+h_p_coil,-h_coil/2+d_air,0,-a_core/2-d_coil

vatt,3,1,3,0

vsel,none ! primary coil right block

bloc,w_coil+b_core/2+d_air,b_core/2+d_air,-h_coil/2+d_air+h_p_coil,-h_coil/2+d_air,0,-a_core/2-d_coil

vatt,3,3,3,0

vsel,none ! primary coil behind block

bloc,-w_coil-d_air*2,w_coil+d_air*2,-h_coil/2+d_air+h_p_coil,-h_coil/2+d_air,-a_core/2-d_coil-d_air-w_coil,-a_core/2-d_coil-d_air

vatt,3,2,3,0

wpcs,-1,11

vsel,none ! primary coil left cylinder

Cyli,d_air,w_coil+d_air,h_p_coil,0,90,180

vatt,3,4,3,11

wpcs,-1,12

vsel,none ! primary coil right cylinder

Cyli,d_air,w_coil+d_air,h_p_coil,0,0,90

vatt,3,5,3,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*2,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

alls

vplo

! *** FE model

et,1,236 ! air

mp,murx,1,1

et,2,236 ! core

Tb,bh,2,,40

Tbpt,,100,0.46512

,,200,0.72993

,,300,0.9009

,,400,1.0204

,,500,1.1086

,,600,1.1765

,,700,1.2302

,,800,1.2739

,,900,1.3

,,1000,1.3405

,,1400,1.4257

,,1800,1.4778

,,2200,1.5131

,,2600,1.5385

,,3000,1.5576

,,3400,1.5726

,,3800,1.5847

,,4200,1.5945

,,4600,1.6028

,,5000,1.6098

,,7000,1.6332

,,9000,1.6465

,,11000,1.6551

,,13000,1.6611

,,15000,1.6656

et,3,236,2 ! primary coil

mp,murx,3,1

r,1,care_p,Np,volu_p_b/2,0,0,1

rmore,Rp*volu_p_b/(volu_p_c+volu_p_b1+volu_p_b)/2,symm

r,2,care_p,Np,volu_p_b1,-1,0,0

rmore,Rp*volu_p_b1/(volu_p_c+volu_p_b1+volu_p_b),symm

r,3,care_p,Np,volu_p_b/2,0,0,-1

rmore,Rp*volu_p_b/(volu_p_c+volu_p_b1+volu_p_b)/2,symm

r,4,care_p,Np,volu_p_c/2,0,1,0

rmore,Rp*volu_p_c/(volu_p_c+volu_p_b1+volu_p_b)/2,symm

r,5,care_p,Np,volu_p_c/2,0,1,0

rmore,Rp*volu_p_c/(volu_p_c+volu_p_b1+volu_p_b)/2,symm

et,4,236,2 ! secondary coil

mp,murx,4,1

r,6,care_s,Ns,volu_s_b/2,0,0,-1

rmore,Rs*volu_s_b/(volu_s_c+volu_s_b1+volu_s_b)/2,symm

r,7,care_s,Ns,volu_s_b1,1,0,0

rmore,Rs*volu_s_b1/(volu_s_c+volu_s_b1+volu_s_b),symm

r,8,care_s,Ns,volu_s_b/2,0,0,1

rmore,Rs*volu_s_b/(volu_s_c+volu_s_b1+volu_s_b)/2,symm

r,9,care_s,Ns,volu_s_c/2,0,-1,0

rmore,Rs*volu_s_c/(volu_s_c+volu_s_b1+volu_s_b)/2,symm

r,10,care_s,Ns,volu_s_c/2,0,-1,0

rmore,Rs*volu_s_c/(volu_s_c+volu_s_b1+volu_s_b)/2,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

Smrtsize,6

vMESH,ALL

alls

! *** boundary conditions and loads

Csys ! flux parallel exterior

asel,s,loc,y,-h_core/2-d_dmn

asel,a,loc,y,h_core/2+d_dmn

asel,a,loc,x,-w_core/2-d_dmn

asel,a,loc,x,w_core/2+d_dmn

asel,a,loc,z,-a_core/2-d_dmn*2

asel,a,loc,z,0

da,all,az

vsel,s,mat,,3 ! primary coil

alls,belo,volu

cp,1,emf,all

cp,2,volt,all

nd_p=ndnext(0)

alls

vsel,s,mat,,4 ! secondary coil

alls,belo,volu

cp,3,emf,all

cp,4,volt,all

nd_s=ndnext(0)

alls

d,nd_p,volt,%u%

! *** Circuit

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

r,11,R

*get,nmax,node,,num,max

n,nmax+1,0,h_coil/2

type,5

real,11

e,nd_s,nmax+1

d,nmax+1,volt,0 ! ground

csys

eplo

fini

! *** solution

/solu

antype,trans

NROPT,FULL

time,tend

nsub,nsubs

outres,all,all

kbc,1 ! step applied load

solve

fini

! *** post-solution

/post26

Nsol,2,nd_p,volt,Vp

Nsol,3,nd_s,volt,Vs