halfsunny 发表于 2007-10-15 22:17:53

大侠赐教:关于变压器电磁场2D瞬态分析

出现一些错误,可是找不出来,急,谢谢啦,bow

*create,SetData
/prep7
*SET,NL,19                        ! NUMBER OF TURNS IN COIL1
*SET,NH,38                     ! NUMBER OF TURNS IN COIL2
*SET,SYMFAC,1                ! SYMMETRIC FACTOR FOR INDUCTANCE COMPUTATION
*SET,NC,2                  ! NUMBER OF COILS
*SET,NN,10                        ! MESHING PARAMETER
*SET,pi,3.1415926   
*SET,H,0.83                     !
*SET,W,0.093                  !
*SET,Rcore,0.105                  !
*SET,He,0.10                  !
*SET,Wpan,0.10               !
*SET,HcoilL,0.758                !
*SET,HcoilH,0.774                !

*SET,Ht,0.185               !
*SET,Wt,0.079               !
*SET,Td,0.01                !油箱壁厚度

Htin=H+He+He+Ht+Ht          !油箱内侧高度
Wtin=W+Rcore+Wpan+Wt      !油箱内侧宽度
Htout=Htin+Td+Td            !油箱外侧高度
Wtout=Wtin+Td               !油箱外侧宽度

RcoilLin=0.12
RcoilLout=0.1415            !
Hc1=0.036                   !
RcoilHin=0.1585             !
RcoilHout=0.177             !
Hc2=0.028

/PREP7
EMUNIT,MKS
MP,MURX,1,1             ! 空气
MP,MURX,2,1             ! 低压COIL1
MP,RSVX,2,2.135E-8
!MP,RSVX,2,2.135E-20      !低压线圈电导率取接近0
MP,MURX,3,1             ! 高压COIL2
MP,RSVX,3,1.72E-8
!MP,RSVX,3,1.72E-20       !高压线圈电导率取接近0
!MP,MURX,4,3572.85          ! 铁心
TB,BH,4,,
TBPT,,6.40000000    , 0.300000000   
TBPT,,8.80000000    , 0.600000000   
TBPT,,11.0000000    , 0.900000000   
TBPT,,14.0000000    , 1.100000000   
TBPT,,25.0000000    , 1.400000000   
TBPT,,70.0000000    , 1.650000000   
TBPT,,200.000000    , 1.740000000   
TBPT,,400.000000    , 1.795000000   
TBPT,,600.000000    , 1.820000000   
TBPT,,1000.00000    , 1.850000000   
TBPT,,2000.00000    , 1.900000000   
TBPT,,5000.00000    , 1.960000000   
TBPT,,6000.00000    , 1.975000000   
TBPT,,10000.0000    , 2.010000000   
tbplot,bh,4
!MP,MURX,4,1                !空心变压器
MP,MURX,5,1000             !定义油箱壁铁芯材料
MP,RSVX,5,1.7E-7
/PNUM,AREA,1
/PNUM,LINE,1
RECTNG,0,Wtout,-Ht-Td,Htout-Ht-Td
RECTNG,0,Wtin,-Ht,Htin-Ht
RECTNG,0,W+Wpan+Rcore,0,H+He*2
RECTNG,Rcore,Rcore+W,He,He+H
RECTNG,RcoilLin,RcoilLout,Hc1+He,He+Hc1+HcoilL
RECTNG,RcoilHin,RcoilHout,Hc2+He,He+Hc2+HcoilH
AOVLAP,ALL
NUMCMP,area
NUMCMP,line
NUMCMP,kp

ASEL,S,AREA,,1   !低压线圈COIL1
AATT,2,1,2
ASEL,S,AREA,,2   !高压线圈COIL2
AATT,3,2,2

ASEL,S,AREA,,5
AATT,4,,1
ASEL,S,AREA,,6
AATT,1,,1
ASEL,S,AREA,,3
AATT,5,,1
ASEL,S,AREA,,4
AATT,1,,1
*end

*create,meshin
/prep7

ET,1,53,,,1         ! AIR
ET,2,53,3,,1      ! CIRCUIT-COUPLED STRANDED COIL
*use,SetData
ALLSEL,ALL
SMRTSIZE,2
MSHAPE,1,2d
MSHKEY,0
AMESH,ALL
R,1,A,NL,,1,1
R,2,B,NH,,1,1

*end


/prep7
*use,meshin

ESEL,S,MAT,,3            ! GET COIL ELEMENTS
CM,COIL_2,ELEM
NSLE,S
CP,2,CURR,ALL            ! COUPLE CURR DOF IN H COIL
!CP,3,EMF,ALL             ! COUPLE EMF DOF IN H COIL



ESEL,S,MAT,,2            ! GET COIL ELEMENTS
CM,COIL_1,ELEM
NSLE,S
!CP,1,CURR,ALL            ! COUPLE CURR DOF IN L COIL
CP,1,EMF,ALL             ! COUPLE EMF DOF IN L COIL

FINISH

/sol
antype,trans
alls
esel,,ename,,53
nsle
nsel,r,ext
d,all,,az,0                              !为节点施加dof约束
alls

*DO,t_increase,0.001,0.001,0.001

*SET,V,10000*sqrt(2)*sin(100*pi*t_increase)
CMSEL,S,COIL_2
BFE,ALL,VLTG,1,V
TIME,0.001+t_increase/10000

alls
solve
OUTRES,ALL,ALL
*ENDDO
save
fini

/post1
SET,1,LAST,1,0, , ,
PLF2D,27,0,10,1             !画磁力线分部图

FINISH

=========
问题是什么?

[ 本帖最后由 北极熊甲 于 2007-10-15 22:30 编辑 ]

sagna 发表于 2012-2-28 16:40:55

你好解决了没
页: [1]
查看完整版本: 大侠赐教:关于变压器电磁场2D瞬态分析