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

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

没解决