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