3D非线性场路耦合不收敛的问题

doublefracture

最近做一些关于可调电抗的仿真计算，发现模型的emf conv总是不收敛。决定简化模型找问题，因此建立了一个简单的线圈-铁芯-空气的场路耦合模型，模型目的是让一个线圈给铁芯交流励磁。但是仍旧出现计算不收敛的情况，circu124电流源和电压源是直流时，求解没有问题，但一旦换成交流，计算则不收敛（具体表现为电流源就电流收敛，电压源就电压收敛）。改了时间步长，收敛准则，均无效。请各位版主大大们帮我看看问题出在哪儿，无限感激~！
命令流如下：（此例网格很粗，但是我复杂电抗器模型是映射划分的，仍旧不收敛）
/prep7
T=1
pi=2*asin(1)
!**********铁芯铁轭**********
rx1=0.5
!中间铁芯的截面半径
rx2=rx1+0.45
!边沿铁芯内侧与中心轴的距离
rx3=rx2+rx1
!边沿铁芯外侧与中心轴的距离
he1=rx1+0.76
!上铁轭上沿高度
he2=he1-rx1
!铁轭半窗高
be=rx1
!铁轭厚度
!**********直流线圈**********
rzi=1.28/2
!直流线圈内半径
rzo=1.34/2
!直流线圈外半径
dze=he2-0.44
!直流线圈上表面与上铁轭下沿距离
hz=0.44
!直流线圈半高度
vz=pi*(rzo*rzo-rzi*rzi)*hz*2
nz=470
!**********交流线圈**********
rji=0.4
!交流线圈内半径
rjo=0. 5
!交流线圈外半径
hj=1.07
!交流线圈高度
xj=rx2+rx1/2
!交流线圈中心x坐标
yj=0
!交流线圈中心y坐标
zj=0
!交流线圈中心z坐标
vj=pi*(rjo*rjo-rji*rji)*hj*2
nj=40
!**********建模**************
block,-rx3,rx3,be/2,-be/2,he1,-he1
block,rx1,rx2,be/2,-be/2,he2,-he2
vsbv,1,2,,,delete
numcmp,all
block,-rx1,-rx2,be/2,-be/2,he2,-he2
vsbv,1,2,,,delete
numcmp,all
cylind,rzi,rzo,hz,-hz,0,360,
block,-3,3,-1,1,-2,2
allsel,all
vovlap,all
numcmp,all
et,1,97
!mp,murx,1,2000
tb,bh,1,,25
tbpt,,100,0.3321
tbpt,,200,0.7134
tbpt,,300,0.9
tbpt,,400,0.9983
tbpt,,500,1.0936
tbpt,,600,1.1598
tbpt,,700,1.2287
tbpt,,800,1.2725
tbpt,,900,1.3008
tbpt,,1000,1.3385
tbpt,,1400,1.4220
tbpt,,1800,1.4699
tbpt,,2200,1.5131
tbpt,,2600,1.5345
tbpt,,3000,1.5528
tbpt,,3400,1.5678
tbpt,,3800,1.5811
tbpt,,4200,1.5890
tbpt,,4600,1.6
tbpt,,5000,1.6045
tbpt,,7000,1.6300
tbpt,,9000,1.6415
tbpt,,11000,1.6550
tbpt,,13000,1.6610
tbpt,,15000,1.6656


vsel,s,,,1
vatt,1,,1
allsel,all
et,2,97,3
mp,murx,2,1
mp,rsvx,2,1e-3
r,2,(rzo-rzi)*2*hz,nz,vz,0,1,0,1,1
vsel,s,,,2
csys,0
local,11,1
vatt,2,2,2,11
allsel,all
et,3,97
mp,murx,3,1
vsel,s,,,3
vatt,3,,3,0
allsel,all

!***********分网**************
Csys,0
MSHKEY,0
!设置自由网格划分
MSHAPE,1,3D
!设置网格和单元为四面体单元




Smrt,6
VMESH,1
Vmesh,2
Vmesh,3
Allsel,all
!划分网格


!****************************电路建模
csys,0
et,6,circu124,3,1
r,6,,300,50
n,,-2,-1.5,0
n,,-2,-2.5,0
type,6
real,6
n,,-2,-2,0
*get,nmax,node,0,num,max
e,nmax-2,nmax-1,nmax
et,7,circu124,0
r,7,1
n,,0,-1.5,0
type,7
real,7
e,nmax-2,nmax+1
et,8,circu124,5
r,8,1
type,8
real,8
e,nmax+1,nmax-1,1108
/sol
allsel,all
ESEL,S,TYPE,,2


NSLE,S

cp,1,curr,all
cp,2,emf,all
allsel,all
/sol
!ESEL,S,TYPE,,1

!NSLE,S

!d,all,volt,0
!allsel,all
!d,nmax-1,volt,0


/sol

CSYS,0
alls
asel,s,loc,x,3
asel,a,loc,z,-3
asel,a,loc,y,-1
asel,a,loc,y,1
asel,a,loc,z,-2
asel,a,loc,z,2
da,all,asym
allsel,all
/solue
antype,4
kbc,0
neqit,100
autots,on
cnvtol,volt,,0.05,,1e-6
cnvtol,curr,,0.05,,1e-6
cnvtol,emf,,0.05,,1e-6
eqslv,sparse
time,t
solve
finish


*do,c,1,10
/sol
Con=0.00001
antype,,rest
t=t+con
kbc,0
c=c+1
deltim,0.000001
outres,all,last
outpr,all,all
alls
time,t
!cnvtol,volt,,0.05,,1e-6
!cnvtol,curr,,0.05,,1e-6
!cnvtol,emf,,0.05,,1e-6
solve
finish
*enddo