变压器场路耦合不收敛问题
三相变压器采用三维建模，副边通过场路耦合有电阻，模拟负载，在原边通过bfe命令加入电密，采用瞬态分析，但随着时间步，副边线圈中的电密和铁芯中的磁密会越来越大，超出ansys的最大值限制而出错。
下面是命令流和报错，请大家指导一下。
et,2,solid97,0
et,3,solid97,0
et,4,solid97,0
mp,murx,1,1 !定义交流线圈相对磁导率为1
mp,rsvx,1,1.75e-2 !定义电阻率
!MP,MURX,2,10000
mp,murx,3,1!定义空气相对磁导率为1
et,1,solid97,3
r,1,Ht*(Rt2-Rt1),N,PI*(Rt2*Rt2-Rt1*Rt1)*Ht,0,1,0,1,1 !副边线圈实常数
r,8,Ht*(Rt2-Rt1),2*N,PI*(Rt2*Rt2-Rt1*Rt1)*Ht,0,-1,0,1,1!原边线圈实常数
...................建模省略
tb,bh,2,,29 !定义铁心的磁化曲线 fx
tbpt,,7.5,0.42
tbpt,,8.25,0.53
tbpt,,10,0.8
tbpt,,12,0.98
tbpt,,15,1.17
tbpt,,20,1.31
tbpt,,32,1.46
tbpt,,38,1.51
tbpt,,44,1.55
tbpt,,60,1.62
tbpt,,70,1.64
tbpt,,80,1.66
tbpt,,90,1.67
tbpt,,100,1.68
tbpt,,150,1.71
tbpt,,200,1.73
tbpt,,290,1.76
tbpt,,380,1.78
tbpt,,500,1.8
tbpt,,750,1.83
tbpt,,1000,1.85
tbpt,,1500,1.88
tbpt,,2000,1.9
tbpt,,3500,1.935
tbpt,,5000,1.96
tbpt,,10000,2
ALLSEL,ALL
vovlap,all
VSEL,S,volu , ,1,5,2
vatt,1,1,1!赋线圈材料属性
VSEL,S,VOLU,,2,6,2
VATT,1,8,2
VSEL,S,VOLU,,37,72,1
VATT,2,,3 !赋铁芯材料属性
VSEL,S,VOLU,,8
VATT,3,,4 !赋边界空气属性
VSEL,S,MAT,,1 !剖分线圈
smrt,6
mshkey,0
mshape,1,3d
vmesh,all
VSEL,S,MAT,,2
smrt,6
mshkey,0
mshape,1,3d
vmesh,all !剖分铁芯
VSEL,S,MAT,,3
vmesh,all !剖分空气
wpcsys,-1,0 !定义线圈柱坐标系
wproat,,-90
wpoff,Wi1/2,-Lbg-Lcd-Lef/2,Hy+Ht+0.011
CSWPLA,12,1,1,1,
wpcsys,-1
vsel,s,,,1,2
ESLV,S
emodif,all,esys,12,
wpcsys,-1,0
wproat,,-90
wpoff,Wi1*1.5+Ww1,-Lbg-Lcd-Lef/2,Hy+Ht+0.011
CSWPLA,13,1,1,1,
wpcsys,-1
vsel,s,,,3,4
ESLV,S
emodif,all,esys,13,
wpcsys,-1,0
wproat,,-90
wpoff,Wi1*2.5+Ww1+Ww2+Ww3+Wi2,-Lbg-Lcd-Lef/2,Hy+Ht+0.011
CSWPLA,14,1,1,1,
wpcsys,-1
vsel,s,,,5,6
ESLV,S
emodif,all,esys,14,
VSEL,S,VOLU , ,1
NSLV,S,1
CP,1,EMF,all
CP,2,CURR,all
CM,W1,ELEM
R,2, !场路耦合
*GET,_zz2,active,,csys
*GET,_zz3,active,,dsys
CSYS,0
DSYS,0
N,20000,-0.2,0.7,0
N,20001,0.2,0.7,0
CSYS,_zz2
DSYS,_zz3
RMOD,2,15,0,5
ET,5,CIRCU124,5,0
TYPE,5
REAL,2
MAT,1
E,20000,20001,1397
R,5,10000, !副边电阻
*GET,_zz2,active,,csys
*GET,_zz3,active,,dsys
CSYS,0
DSYS,0
CSYS,_zz2
DSYS,_zz3
RMOD,5,15,-0.1,8
ET,8,CIRCU124,0,0
TYPE,8
REAL,5
MAT,1
E,20000,20001
VSEL,S,VOLU , ,3
NSLV,S,1
CP,3,EMF,all
CP,4,CURR,all
CM,W3,ELEM
R,3,
*GET,_zz2,active,,csys
*GET,_zz3,active,,dsys
CSYS,0
DSYS,0
N,20001,0.2,0.7,0
N,20002,0.6,0.7,0
CSYS,_zz2
DSYS,_zz3
RMOD,3,15,0,6
ET,6,CIRCU124,5,0
TYPE,6
REAL,3
MAT,1
E,20001,20002,2853
!E,20001,20002,5003 !自由剖分
R,6,10000,
*GET,_zz2,active,,csys
*GET,_zz3,active,,dsys
CSYS,0
DSYS,0
CSYS,_zz2
DSYS,_zz3
RMOD,6,15,-0.1,9
ET,9,CIRCU124,0,0
TYPE,9
REAL,6
MAT,1
E,20001,20002
VSEL,S,VOLU , ,5
NSLV,S,1
CP,5,EMF,all
CP,6,CURR,all
CM,W3,ELEM
R,4
*GET,_zz2,active,,csys
*GET,_zz3,active,,dsys
CSYS,0
DSYS,0
N,20001,0.2,0.7,0
N,20003,0.2,0.3,0
CSYS,_zz2
DSYS,_zz3
RMOD,4,15,0,7
ET,7,CIRCU124,5,0
TYPE,7
REAL,4
MAT,1
E,20001,20003,4240
!E,20001,20003,7952
R,7,10000,
*GET,_zz2,active,,csys
*GET,_zz3,active,,dsys
CSYS,0
DSYS,0
CSYS,_zz2
DSYS,_zz3
RMOD,7,15,-0.1,10
ET,10,CIRCU124,0,0
TYPE,10
REAL,7
MAT,1
E,20001,20003
WPCSYS,-1,0
finish
/solu
NSEL,S,NODE,,20001
D,ALL,VOLT,0
ASEL,S,,,48,53,1
DA,all,az,0,0 !边界条件
allsel,all
antype,4
trnopt,full
neqit,30,
cnvtol,csg,1000,0.1,2, ,
outres,all,all,
*DO,I,1,50
IA=10*COS(800*PI*t*i)
IB=10*COS(800*PI*t*i+PI*2/3)
IC=10*COS(800*PI*t*i+PI*4/3)
ID_A=N*IA/Ht/(Rt2-Rt1)!三相电流密度
ID_B=N*IB/Ht/(Rt2-Rt1)
ID_C=N*IC/Ht/(Rt2-Rt1)
vsel,S,volu,,2 !加载电流密度
ESLV,S
BFE,all,JS,1,,ID_A
vsel,S,volu,,4
ESLV,S
BFE,all,JS,1,,ID_B
vsel,S,volu,,6
ESLV,S
BFE,all,JS,1,,ID_C
time,t*i
autots,0
deltim,t, , ,1
kbc,1
tsres,erase
allsel,all
solve
不好意思，报错的图片没发成功，我接着发一下！
*** ERROR *** CP = 159.406 TIME= 12:09:14
The value of EMF at node 6798 is 1.754987335E+26.It is greater than
the current limit of 1.E+25.This generally indicates rigid body
motion as a result of an unconstrained model.Verify that your model
is properly constrained. 这个错误一般都是因为约束加的不对，你再查查约束 没解决
页:
[1]