ANSYS三维瞬态仿真电抗器，铁芯加入BH曲线，发现结果不收敛，求大神指导
各位大神求指导啊，我的ansys仿真问题大概是这样的：
我采用场路耦合+瞬态进行的电抗器的电磁仿真。仿真方法采用的是节点矢量法。单元采用的是solid97，铁芯加入了BH曲线。铁芯的keyopt(1)=0,绕组的keyopt(1)=3用于模拟载压绞线圈，线圈激励是通过场路耦合的外部电路给入。然后仿真时出现了不收敛，提示节点的EMF过大。改变仿真的时间步长感觉依然出现不收敛，不知道为什么。但如果铁芯的磁导率变为常数，则收敛性很好。看帮助文件说solid97单元，如果keyopt(1)=3,则必须磁导率为常数。但不知道是不是说含有载压绞线圈时，铁芯不能加BH曲线。而实际上solid97是可以含有BH曲线的。所以想知道应该怎么调节参数才能使得仿真能够收敛？跪求各位大神指导啊
以下为ansys仿真程序：
finish/clear/TITLE,reductor/FILNAM,reductor_3D,1
/prep7
t1=0.03ii5=60jiedian=(t1/0.02)*10
ET,1,solid97,0 !!!!!空气ET,2,solid97,0 !!!!!铁芯ET,3,solid97,3 !!!!!绕组1ET,4,solid97,3 !!!!!绕组2ET,5,solid97,3 !!!!!绕组3ET,6,solid97,3 !!!!!绕组4
MP,MURX,1,1 !空气磁导率
!铁芯tb,bh,2,,25TBPT,,100,0.46512 TBPT,,200,0.72993TBPT,,300,0.90090TBPT,,400,1.0204TBPT,,500,1.1086TBPT,,600,1.1765TBPT,,700,1.2302TBPT,,800,1.2739TBPT,,900,1.3100TBPT,,1000,1.3405TBPT,,1400,1.4257TBPT,,1800,1.4778TBPT,,2200,1.5131TBPT,,2600,1.5385TBPT,,3000,1.5576TBPT,,3400,1.5726TBPT,,3800,1.5847TBPT,,4200,1.5945TBPT,,4600,1.6028TBPT,,5000,1.6098TBPT,,7000,1.6332TBPT,,9000,1.6465TBPT,,11000,1.6551TBPT,,13000,1.6611TBPT,,15000,1.6656
MP,MURX,3,1 !!!绕组1材料属性MP,RSVX,3,1.79e-8MP,MURX,4,1 !!!绕组2材料属性MP,RSVX,4,1.79e-8MP,MURX,5,1 !!!绕组3材料属性MP,RSVX,5,1.79e-8MP,MURX,6,1 !!!绕组4材料属性MP,RSVX,6,1.79e-8
a_rz=0.02*0.1 !!!!!!!!!!!!绕组横截面积
v_rz=0.628E-3 !!!!!!!!!!!!绕组体积
!参数n=80r1=0.035 !铁芯半径r2=0.04 !绕组内半径r3=0.06 !绕组外半径d1=0.15 !铁芯间距d2=0.04 !旁轭间距h1=0.3 !铁芯高度h2=0.1 !铁轭高度
r,3,a_rz,n,v_rz,0,1,0,1,0.95r,4,a_rz,n,v_rz,0,1,0,1,0.95r,5,a_rz,n,v_rz,0,1,0,1,0.95r,6,a_rz,n,v_rz,0,1,0,1,0.95
!!!!!!!!!!!!!!铁芯cylind,r1,0,0,h1 !主铁芯1
csys,1wpoffs,d1cylind,r1,0,0,h1 !主铁芯2
wpcsys,-1 csys,1block,-0.15,0.3,-0.04,0.04,0,-0.1block,-0.15,0.3,-0.04,0.04,0.3,0.4block,-0.15,-0.1,-0.04,0.04,0,0.3block,0.25,0.3,-0.04,0.04,0,0.3 !铁轭和旁轭allsvglue,allallsvadd,all !!!!!!将铁芯的体相加numcmp,volu!!cm,iron,voluallsvatt,2,,2
!!!!!!!!!!!!绕组cylind,r2,r3,0.03,0.13 !绕组1allsvsel,u,mat,,2cm,rz1,voluvatt,3,3,3
cylind,r2,r3,0.17,0.27 !绕组2allsvsel,u,mat,,2,3cm,rz2,voluvatt,4,4,4
csys,1wpoffs,d1
cylind,r2,r3,0.03,0.13 !绕组3allsvsel,u,mat,,2,4cm,rz3,voluvatt,5,5,5
cylind,r2,r3,0.17,0.27 !绕组4allsvsel,u,mat,,2,5cm,rz4,voluvatt,6,6,6
wpcsys,-1 csys,1!!!!!!!!!!!!!!绕组附近空气1cylind,0.035,0.06,0.03,0.27allsvsel,u,mat,,2vovlap,allallsvsel,u,mat,,2,6vatt,1,,1
!!!!!!!!!!!!!!绕组附近空气1csys,1wpoffs,d1cylind,0.035,0.06,0.03,0.27allsvsel,s,,,4,6 !!!!此处坐标为原坐标vovlap,allallsvsel,u,mat,,1,6vatt,1,,1
wpoffs,-0.075 !找到中性线csys,1
allsvplot!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!外层空气allscylind,0,0.4,-0.2,0.5cylind,0,1,-1,2allsvovlap,allallsvsel,s,,,12vatt,2,,2allsvsel,u,mat,,1,6cm,air,voluvatt,1,,1
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!meshallsvsel,s,mat,,3,6esize,0.01,0mshape,1,3dmshkey,0vmesh,all
allsvsel,s,,,10,11esize,0.01,0mshape,1,3dmshkey,0vmesh,all
allsvsel,s,mat,,2esize,0.02,0mshape,1,3dmshkey,0vmesh,all
allsvsel,s,,,14esize,0.1,0,mshape,1,3dmshkey,0vmesh,all
allsvsel,s,,,13esize,0.2,0mshape,1,3dmshkey,0vmesh,all
save,'mesh',db!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!寻找路耦合节点的编号*do,i,1,4,1allsvsel,s,mat,,i+2nslv,s,1*get,nn_rz%i%,node,,num,min !!!!!!!!!!!!!!!!!!寻找绕组耦合节点编号*ENDDO
et,111,124,4,1 ! ind vlgt src sinet,112,124,4 ! dc piecewise linearet,7,124,5 ! scal sourceet,8,124,0 ! resistoret,11,124,1 ! inductor
r,111,,311,50,,90 !220*1.414=311r,112,ii5 !改变此处可以加上直流r,7,1 !scal=1 r,8,1 !Rs=1r,9,1 !Rload=1r,10,2 !Rd=2r,11,3.4 !Ls=0.0136r,12,37.72 !Lload=0.1233r,13,0.1 !Ld=0.02
allscsys,0
!生成节点，用于连接电路*get,nm,node,,num,maxnm=nm+1
n,nm+1,-10,0,1n,nm+2,-10,0,0n,nm+3,-10,0,-1n,nm+4,-9,0,1n,nm+5,-9,0,-1n,nm+6,-6,0,1n,nm+7,-8,0,0n,nm+8,-4,0,0n,nm+9,-6,0,-1n,nm+10,-7,0,0n,nm+11,-6,0,0n,nm+12,-5,0,0
real,111 !交流电源type,111e,nm+1,nm+3,nm+2
real,112 !直流电源type,112e,nm+11,nm+8,nm+12
real,8 !resistortype,8e,nm+1,nm+4
real,9type,8e,nm+3,nm+5
real,10type,8e,nm+7,nm+10
real,11 !inductortype,8e,nm+4,nm+6
real,12 !inductortype,8e,nm+5,nm+9
real,13type,8e,nm+10,nm+11
real,7type,7e,nm+6,nm+7,nn_rz2 !绕组耦合连接e,nm+6,nm+8,nn_rz4e,nm+7,nm+9,nn_rz3e,nm+8,nm+9,nn_rz1
!!!!!!!!!!!!!!!!!!!!!!!!!!!!耦合体单元节点*do,i,3,6,1vsel,s,mat,,inslv,s,1cp,next,curr,allcp,next,emf,all*enddo
!!!!!!!!!!!!外边界allsvsel,s,mat,,1nslv,s,1nsel,r,extd,all,ax,0d,all,ay,0d,all,az,0
allsd,nm+3,volt,0
alls!!!!!!!!!!!!求解/SOLUantype,transeqslv,sparseneqit,100outres,all,alltime,t1deltim,0.0005,0.0005,0.001cnvtol,volt,,,,1e-12allssolve
