变压器场路耦合不收敛问题

luomu

三相变压器采用三维建模，副边通过场路耦合有电阻，模拟负载，在原边通过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                                                                                                                                                                              
*enddo
FINISHfile:///d:/My

luomu

不好意思，报错的图片没发成功，我接着发一下！

*** 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.

lilulilu

这个错误一般都是因为约束加的不对，你再查查约束

chens3109160027

请问楼主这个问题最终怎么解决？

wtsyansheng

我也遇到相同的问题不知道楼主怎么解决的啊