电机磁链的问题

azedd

根据网上找到的计算永磁电机磁链的例子，计算结果为
****** COMPUTED INDUCTANCE, FLUX ARRAY ******
  
        I        FLUX
  
        1.    -0.25324E+00
        2.     0.31564E+00
        3.    -0.23154E+00
  
  
        I       J         INDUC
  
        1.      1.     0.18526E+01
        1.      2.     0.18383E+01
        1.      3.     0.18387E+01
        2.      2.     0.18518E+01
        2.      3.     0.18383E+01
        3.      3.     0.18525E+01
  
****** END OF INDUCTANCE ARRAY ******

为什么三相磁链之和不为零呢？模型中是否存在什么问题呢

以下是这个例子的命令流：  
!********************************
! Time-stamp: "2008-04-21 19:16:18"     
!********************************
! Summary: This is a model without loop structure. It works well.
! change nolinear material(B-H) to linear material(murx)

! Problems:
! 1 Contact command is for deformation, not for rotating parts
! 2 When Ansys generates new areas, it will assign that area a new number.
! I can't control that number. Without exact number, I can't manipulate the new areas


finish
/clear


/PREP7
! Enters the general input data preprocessor (PREP7).
/TITLE,Brushless permanent magnet DC motor

! presetting

emunit,mks        ! standard international units μ0=4 Pi e-7 henries/meter
csys,0            ! 0 (default)  —  Cartesian
*afun,deg         ! use degree as unit for angle

!******************************************************************************
!define materials properties
!******************************************************************************

ET,1,PLANE53
MP,MURX,1,1
! material 1 ,air

ET,2,PLANE53,
mp,murx,2,0.431
mp,mgxx,2,-330e3

ET,3,PLANE53,
mp,murx,3,0.431
mp,mgxx,3,+330e3
! permanent magnet,take linear for example,value is intended to be bigger
! material 2 and 3 are both the same magenet. We define it as positive and
! negative mgxx, so that we can ! achieve different polarity. For example N--S--N--S

ET,4,PLANE53,
MP,MURX,4,1  
! copper

ET,5,PLANE53,


TB,BH,5,,30
! TB,material property,number of material,working temperature,30 degree,material number 4
TBPT,, 23.200000 , 0.1
TBPT,, 46.400000 , 0.200000000
TBPT,, 69.600000 , 0.300000000
TBPT,, 92.800000 , 0.400000000
TBPT,, 116.000000 , 0.500000000
TBPT,, 128.000000 , 0.600000000
TBPT,, 141.000000 , 0.700000000
TBPT,, 156.000000 , 0.80000000
TBPT,, 174.000000 , 0.90000000
TBPT,, 197.000000 , 1.00000000
TBPT,, 230.000000 , 1.10000000
TBPT,, 278.000000 , 1.20000000
TBPT,, 362.000000 , 1.30000000
TBPT,, 540.000000 , 1.40000000
TBPT,, 1034.000000 , 1.50000000
TBPT,, 2439.000000 , 1.60000000
TBPT,, 5370.000000 , 1.70000000
TBPT,, 9948.000000 , 1.80000000
TBPLOT,BH,5,,,
/IMAGE,SAVE,BH2,JPEG       !将材料2的B-H曲线存储成bh2.jpg
! steel
! /IMAGE, Label, Fname, Ext, --
! Allows graphics data to be captured and saved.



!*****************************************************************************
! parameters input
!****************************************************************************

! MULTIPRO cannot be called from UIDL. You cannot use MULTIPRO within a DO loop.
! MULTIPRO,'start',Prompt_Num (the inputs that you need, here is 6, 6 inputs)
! *CSET,Strt_Loc (1+3*n,n from 0 to...),End_Loc (3+3*n,n from 0 to... 3*n),Param_Name,
! 'Prompt_String',Def_Value
! MULTIPRO,'end'


d1=127.4  ! 转子外直径
d2=120.4  ! 转子内直径
temp1=5.3   !磁铁厚度
poles=8
d4=51       ! 定子yoke外直径
d5=34       ! 定子yoke内直径
temp2=3  
gap=0.9  ! air gap


! *set only allowed 8 charaters maximum
*DIM,name,char,1,14  ! Index begins with 1, not zero!

! *SET,name(1,1),'Angle(degree)'
! *SET,name(1,2),'Iu [A]'
! *SET,name(1,3),'Iv [A]'
! *SET,name(1,4),'I  [A]'
! *SET,name(1,5),'Luu [H/m]'
! *SET,name(1,6),'Lvv [H/m]'
! *SET,name(1,7),'Lww [H/m]'
! *SET,name(1,8),'Lvu [H/m]'
! *SET,name(1,9),'Lwu [H/m]'
! *SET,name(1,10),'Lwv [H/m]'
! *SET,name(1,11),'Torque [Nm/m]'
! *SET,name(1,12),'energy [J/m]'
! *SET,name(1,13),'coenergy [J/m]'
! *SET,name(1,14),'--------'


*SET,name(1,1),'Angle'
*SET,name(1,2),'Iu'
*SET,name(1,3),'Iv '
*SET,name(1,4),'I  '
*SET,name(1,5),'Luu '
*SET,name(1,6),'Lvv '
*SET,name(1,7),'Lww '
*SET,name(1,8),'Lvu '
*SET,name(1,9),'Lwu '
*SET,name(1,10),'Lwv '
*SET,name(1,11),'Torque '
*SET,name(1,12),'energy '
*SET,name(1,13),'coenergy '
*SET,name(1,14),'--------'

! *DIM,energy,array,1,2



! multipro,'start',8
! *cset, 1, 3,d1, 'Outer Diameter of rotor(mm)',133.1 !转子外直径
! *cset, 4, 6,d2, 'Inner Diameter of rotor(mm)',127.1 !转子内直径
! *cset, 7, 9,temp1, 'Thickness of magnet(mm)',9 !磁铁厚度
! *cset,10,12,poles,'Pole number of magnet(on Rotor)',10 !磁铁的极数
! *cset,13,15,d4,'Outer Diameter of stator yoke(mm)', 41 !定子yoke外直径
! *cset,16,18,d5, 'Inner Diameter of stator yoke(mm)', 34 !定子yoke内直径
! *cset,19,21,temp2, 'Thickeness of stator shoe(mm)', 3 !定子 shoe thickness
! *cset,22,24,gap, 'Air gap(mm)', 0.8 ! Air gap
! multipro,'end'

r1=d1/2/1000 ! 转子外半径,转化为国际单位制 yoke of rotor
r2=d2/2/1000 ! 转子内半径,转化为国际单位制 yoke of rotor
r3=(d2-2*temp1)/2/1000 ! 转子上的磁铁内半径，转化为国际单位制
r4=d4/2/1000 ! 定子yoke外半径,转化为国际单位制
r5=d5/2/1000 ! 定子yoke内半径,转化为国际单位制
r6=(d2-2*temp1-gap*2)/2/1000 ! 定子shoe外半径,转化为国际单位制，air gap 0.8 mm
r7=(d2-2*temp1-gap*2-temp2*2)/2/1000 ! 定子shoe内半径,转化为国际单位制

! basic loop structure
! current in phase from -9A to 9A, step size 3A
! mechanical rotating degree from 0 tp 72 degree, step size 3 degrees

m_angle=72  ! mechanical rotating angle
a_step=3    ! step size for rotating angle, 3 degrees/step
i= m_angle/a_step  ! increment of angle loop

current_range=9+9  ! current in every phase is from -9 A to 9 A
c_step=3           ! step sie for current, 3 A/step
j=current_range/c_step  ! increment of current loop

! We need to calculate the current density from the phase current
! 4 coil pairs per phase, they are parallel connected. 27 turns per coil pair
! square of coil=??
! current density in every coil=phase current/(4*square of coil)
! one good example is in script of small motor, page 6-16
! we must define one variable for the square of coil
! pi*(50.75^2-20.5^2)*(28/360)-30*8= 286.6 mm^2= 286.6/1e6   transfer unit into m^2

! area_coil=4*286.6/1e6

! from datasheet
! area_coil=4*31.53/1e6
! area_coil=4*31.53/1e8

! current desnity in phase = phase current / area_coil

! jdens1=1300
! jdens2=900
! angle1=15
angle2=3.62  ! the angle between 2 near magnets


! multipro,'start',4
! *cset, 1, 3,jdens1, 'Current density(winding U)', 1300 !定子线圈上的电流密度
! *cset, 4, 6,jdens2, 'Current density(winding V)', 900 !
! *cset, 7, 9,angle1, 'rotating angle of rotor(degree)', 5 ! 转子的转动角度
! *cset, 10, 12,angle2, 'magnet displacement(degree)',1 ! 转子磁铁之间的角度displacement
! multipro,'end'

! winding W, summ of all 3 phase current should be zero
! 12 windings total, one phase has 4 windings, they are all parallel connected
! current density is positive in 2 windings, negative in other 2 windings

! jdens3=-(jdens1+jdens2)

! current in winding U: jdens1
! current in winding V: jdens2
! current in winding W: jdens3

!*************************************************************************************
! building model, stator part
!*************************************************************************************

/PNUM,AREA,1   
! display area number


! PCIRC,r6,r4,0,27.5,

PCIRC,r6,r7,0,27.5,  ! shoe part
wpoff,37.15/1000,9.09/1000  
wprot,13.75   ! change the WP to center, for plotting rect
BLC5,0/1000,0/1000,30/1000,8.2/1000
wprot,-13.75
wpoff,-37.15/1000,-9.09/1000  


PCIRC,r7,r4,4.17,23.33,
aovlap,all

ASEL,S,AREA,,7  
ASUM, fine
!  TOTAL SURFACE AREA OF ALL SELECTED AREAS =  0.58295E-04
! unit in m^2, it is 58.295 mm^2. This meets the datasheet.

allsel
aplot,all

ASEL,S,AREA,,4
ASEL,A,AREA,,9
ASEL,A,AREA,,8
ASEL,A,AREA,,5
AADD,ALL
! build the tooth and shoe

allsel
aplot,all

NUMCMP,AREA
! compress the area number

aplot


csys,1
! 1  —  Cylindrical with Z as the axis of rotation, for agen!

ASEL,S,AREA,,1
ASEL,A,AREA,,2
ASEL,A,AREA,,3




agen,12,all,,,,360/12
! generate all 12 teeth and shoes on stator

allsel
aplot,all
! have a full view of what we have drawn


PCIRC,r4,r5,0,360,
! draw the stator yoke,area 38

! allsel,
! aplot
! only for easy selection of area 2

ASEL,S,AREA,,37
ASEL,A,AREA,,1
ASEL,A,AREA,,4
ASEL,A,AREA,,7
ASEL,A,AREA,,10
ASEL,A,AREA,,13
ASEL,A,AREA,,16
ASEL,A,AREA,,19
ASEL,A,AREA,,22
ASEL,A,AREA,,25
ASEL,A,AREA,,28
ASEL,A,AREA,,31
ASEL,A,AREA,,34
AADD,ALL
! adding yoke and teeth together

allsel
aplot,all
! have a full view of what we have drawn

!*************************************************************************************
! building model, air gap part, 2 layers; stator part(fixed) & rotor part(totated)
!*************************************************************************************
! We split air gap into 2 parts. We are going to build the air gap in stator part.
! It is easier for the later meshing & unmeshing


! part 1,stator part, fixed, radius from(small, big) r6,r6+(r3-r6)/2, through overlap r4,r6+(r3-r6)/2
PCIRC,r5,r6+(r3-r6)/2,0,360,
aovlap,all
! air gap A4



!*************************************************************************************
! building model, rotor part; yoke
!*************************************************************************************

PCIRC,r1,r2,0,360,
! yoke of rotor,without magnet

!*************************************************************************************
! building model, rotor part; magnets
!*************************************************************************************

angle1=0
! rotating angle begins from 3,value 3,6,9,12,........72

PCIRC,r2,r3,0+angle1,(360/poles)+angle1-angle2,
! the first magnet

ASEL,S, , ,      10
! select the newly generated magnet area,A10

agen,poles,all,,,,360/poles

allsel,
aplot,
! select all areas, so that we can have a clear view of whole structure.

!*************************************************************************************
! building model, air gap part, 2 layers; stator part(fixed) & rotor part(rotated)
!*************************************************************************************
! We split air gap into 2 parts. We are going to build the air gap in stator part.
! It is easier for the later meshing & unmeshing


! part 2, rotor part,rotating,radius from(small, big) r6+(r3-r6)/2,r2 through overlap
PCIRC,r6+(r3-r6)/2,r2,0,360,
aovlap,all

! air gap, rotor part, A46
! air gap, stator part, A4
! stator yoke+teeth+shoes, A7

!*************************************************************************************
! compress the area numbers
!*************************************************************************************
NUMCMP,AREA

! air gap, rotor part, A36
! air gap, stator part, A3
! stator yoke+teeth+shoes, A6
! rotor yoke, A35


!*************************************************************************************
!glue all areas, so that flux lines can pass through
!*************************************************************************************

aglu,all
! glue all areas together, so that flux line can pass through
! area number
! air gap, rotor part, A36
! air gap, stator part, A3
! stator yoke+teeth+shoes, A6
! rotor yoke, A35


!*************************************************************************************
! create 8 local coordination systems, for 8 poles
! polarities magnets. Our goal is to achieve switched polarity N--s--N--s......
!*************************************************************************************

! LOCAL, KCN, KCS, XC, YC, ZC, THXY, THYZ, THZX, PAR1, PAR2   KCS=1, cylindarical
! Defines a local coordinate system by a location and orientation.
! the first pole center is at (36-angle2)/2+angle1+n*36, n from 0 to 9

local,11,1,,,,((360/poles)-angle2)/2+angle1

local,12,1,,,,((360/poles)-angle2)/2+angle1+(360/poles)

local,13,1,,,,((360/poles)-angle2)/2+angle1+(360/poles)*2

local,14,1,,,,((360/poles)-angle2)/2+angle1+(360/poles)*3

local,15,1,,,,((360/poles)-angle2)/2+angle1+(360/poles)*4

local,16,1,,,,((360/poles)-angle2)/2+angle1+(360/poles)*5

local,17,1,,,,((360/poles)-angle2)/2+angle1+(360/poles)*6

local,18,1,,,,((360/poles)-angle2)/2+angle1+(360/poles)*7



!*******************************************************************************
! applly material properties
!*****************************************************************************

ASEL, , , ,       6
AATT,       5, ,   5,       0,  
! select stator yoke+teeth+shoes, apply material property No.5,steel

ASEL, , , ,      35
AATT,       5, ,   5,       0,  
! select rotor yoke, apply material property No.5,steel

ASEL, , , ,      36
AATT,       1, ,   1,       0,  
! select rotor air gap, apply material property No.1, air

ASEL, , , ,      3
AATT,       1, ,   1,       0,  
! select stator air gap, apply material property No.1, air

ASEL,S,AREA,,1
ASEL,A,AREA,,2
ASEL,A,AREA,,25
ASEL,A,AREA,,26
ASEL,A,AREA,,23
ASEL,A,AREA,,24
ASEL,A,AREA,,21
ASEL,A,AREA,,22
ASEL,A,AREA,,19
ASEL,A,AREA,,20
ASEL,A,AREA,,17
ASEL,A,AREA,,18
ASEL,A,AREA,,15
ASEL,A,AREA,,16
ASEL,A,AREA,,13
ASEL,A,AREA,,14
ASEL,A,AREA,,11
ASEL,A,AREA,,12
ASEL,A,AREA,,9
ASEL,A,AREA,,10
ASEL,A,AREA,,7
ASEL,A,AREA,,8
ASEL,A,AREA,,4
ASEL,A,AREA,,5
AATT,       4, ,   4,       0,
! selct 24 coils, apply material property No.4, copper


! magnets
ASEL, , , ,      27
AATT,       2, ,   2,      11,  

ASEL, , , ,      28
AATT,       3, ,   3,      12,  

ASEL, , , ,      29
AATT,       2, ,   2,      13,  

ASEL, , , ,      30
AATT,       3, ,   3,      14,  

ASEL, , , ,      31
AATT,       2, ,   2,      15,  

ASEL, , , ,      32
AATT,       3, ,   3,      16,  
  
ASEL, , , ,      33
AATT,       2, ,   2,      17,  

ASEL, , , ,      34
AATT,       3, ,   3,      18,  
  

! select all magnets

allsel  
aplot,all   


!****************************************************************************
!meshing
!*******************************************************************************
/PREP7

SMRT,5  
MSHAPE,1,2D
! 1  —  Mesh with triangle-shaped  
MSHKEY,0
! 0  —  Use free meshing (the default)
ASEL,all
AMESH,all

!************************************************************************************
! define 3 components for 3 phases, U, V, W
! CM_1(phase U),CM_2(phase V),CM_3(phase W)
!************************************************************************************

! Phase U

! winding U, positive current density, 4 areas:
ASEL,S,AREA,,1
ASEL,A,AREA,,9
ASEL,A,AREA,,15
ASEL,A,AREA,,21

! winding U, negative current density, 4 areas:

ASEL,A,AREA,,2
ASEL,A,AREA,,10
ASEL,A,AREA,,16
ASEL,A,AREA,,22


ESLA,S  
! ESLA, Type     ! select all elements associated with the areas
! Selects those elements associated with the selected areas.

CM,CM_1,ELEM
! name selected elements as CM_1

CMSEL,S,CM_1
CMPLOT  
! plot the selected windings(phase U) to see the position



! winding V, positive current density, 4 areas:
ASEL,S,AREA,,4
ASEL,A,AREA,,11
ASEL,A,AREA,,17
ASEL,A,AREA,,23

! winding V, negative current density, 4 areas:
ASEL,A,AREA,,5
ASEL,A,AREA,,12
ASEL,A,AREA,,18
ASEL,A,AREA,,24

ESLA,S  
! ESLA, Type     ! select all elements associated with the areas
! Selects those elements associated with the selected areas.

CM,CM_2,ELEM
! name selected elements as CM_2

CMSEL,S,CM_2
CMPLOT  
! plot the selected windings(phase V) to see the position

! winding W, positive current density, 4 areas:
ASEL,S,AREA,,7
ASEL,A,AREA,,13
ASEL,A,AREA,,19
ASEL,A,AREA,,25

! winding W, negative current density, 4 areas:
ASEL,A,AREA,,8
ASEL,A,AREA,,14
ASEL,A,AREA,,20
ASEL,A,AREA,,26

ESLA,S  
! ESLA, Type     ! select all elements associated with the areas
! Selects those elements associated with the selected areas.

CM,CM_3,ELEM
! name selected elements as CM_3

CMSEL,S,CM_3
CMPLOT  
! plot the selected windings(phase W) to see the position

Allsel
aplot,all

!************************************************************************************
!Define the boundaries, MVP=0
!************************************************************************************

/solu
ANTYPE,0
! static analysis

csys,1
nsel,s,loc,x,r1
D,ALL,AZ,0
! define the boundary of outside diameteer(rotor)

ALLSEL,ALL

lsel,s,loc,x,r5
dl,all,,AZ,0
! define inner diameter(stator yoke, innner side) as boundary

ALLSEL,ALL
aplot,all



!************************************************************************************
!  calculate the current density with phase current
!************************************************************************************

*DIM,CURRENT,ARRAY,3,1,1, , ,   
temp_J=2


current_U=c_step*temp_J-9  ! -3

current_V=-c_step*temp_J+9 ! 3

current_W=-(current_U+current_V)

*if,current_U,EQ,0,then
current_U=1e-3     
! Zero current values must be approximated by a negligibly small applied current.
! 1e-6 is still too small for Lmatrix,how about 1e-3
*endif

*if,current_V,EQ,0,then
current_V=1e-3     
! Zero current values must be approximated by a negligibly small applied current.
*endif

*if,current_W,EQ,0,then
current_W=1e-3     
! Zero current values must be approximated by a negligibly small applied current.
*endif

area_coil=58e-6

jdens1=27*current_U/area_coil

jdens2=27*current_V/area_coil

jdens3=27*current_W/area_coil
! 1.85 H/m, value is too big?

! jdens1=current_U/area_coil

! jdens2=current_V/area_coil

! jdens3=current_W/area_coil
! 0.25556E-02 H/m, but torque is just too small. Doesn't match


*SET,CURRENT(1,1,1) , current_U
*SET,CURRENT(2,1,1) , current_V
*SET,CURRENT(3,1,1) , current_W


!************************************************************************************
!apply force
!************************************************************************************
! force can only be applied on element and node or area line keypoints, but never on
! components

! winding U, positive current density, 4 areas:
ASEL,S,AREA,,1
ASEL,A,AREA,,9
ASEL,A,AREA,,15
ASEL,A,AREA,,21

BFA,all,JS, , ,jdens1,0  

! winding U, negative current density, 4 areas:

ASEL,A,AREA,,2
ASEL,A,AREA,,10
ASEL,A,AREA,,16
ASEL,A,AREA,,22

BFA,all,JS, , ,-jdens1,0

! winding V, positive current density, 4 areas:
ASEL,S,AREA,,4
ASEL,A,AREA,,11
ASEL,A,AREA,,17
ASEL,A,AREA,,23

BFA,all,JS, , ,jdens2,0

! winding V, negative current density, 4 areas:
ASEL,A,AREA,,5
ASEL,A,AREA,,12
ASEL,A,AREA,,18
ASEL,A,AREA,,24

BFA,all,JS, , ,-jdens2,0

! winding W, positive current density, 4 areas:
ASEL,S,AREA,,7
ASEL,A,AREA,,13
ASEL,A,AREA,,19
ASEL,A,AREA,,25

BFA,all,JS, , ,jdens3,0

! winding W, negative current density, 4 areas:
ASEL,A,AREA,,8
ASEL,A,AREA,,14
ASEL,A,AREA,,20
ASEL,A,AREA,,26

BFA,all,JS, , ,-jdens3,0

allsel,
aplot,all

!************************************************************************************
! choose solver and solve for MVP
!************************************************************************************
/solu

NROPT,AUTO, ,     
! NROPT, Option, --, Adptky
! Specifies the Newton-Raphson options in a static or full transient analysis
! auto, let the program to optimize nonlinear calculation

EQSLV,FRONT, ,0,
! EQSLV, Lab, TOLER, MULT
! Specifies the type of equation solver.

PRECISION,0
! PRECISION, LABEL   ! 0 or double    —  Specifies double precision (default).
! Specifies machine precision for solvers (currently valid only for PCG solvers).

MSAVE,0
! MSAVE, Key    ! 0 or OFF  —  Use global assembly for the stiffness matrix of
! the entire model.
! Sets the solver memory saving option. This option only applies to the PCG and
! DPCG solvers.

PIVCHECK,1  
! PIVCHECK, KEY
! Prevents a batch mode, linear static analysis from stopping when a negative or
! zero equation solver pivot value is encountered.

MAGSOLV
! solve the MVP

!************************************************************************************
! calculate inductance
!************************************************************************************

LMATRIX,1, 'CM_' , 'CURRENT' , 'Inductance'
! LMATRIX, SYMFAC, Coilname, Curname, Indname
! Calculates an inductance matrix and the total flux linkage for an N-winding coil system.
! inductance is an array, differential inductance and flux linkage
! careful, the values in inductance array are arranged like
!             1         2          3     flux linage
! phase1    self
! phase2    Mutual     self
! phase3                          self
! You must set each vector array entry equal to the nominal current per turn in the
! corresponding coil at the operating point. Zero current values must be approximated
! by a negligibly small applied current.

!************************************************************************************
! calculate torque, select arears in air gap
!************************************************************************************
/POST1  

SET,LAST
TORQC2D,r6+(r3-r6)/2,180,0, ,
! TORQC2D, RAD, NUMN, LCSYS
! Calculates torque on a body in a magnetic field based on a circular path.
! 36.8498017 N*m/m.
! 180 is a good value for number of node. My computer can't support 360 or higher value
! result is stored in parameter torque

!************************************************************************************
! calculate energy and co-energy
!************************************************************************************


! SENERGY, OPT, ANTYPE
! Determines the stored magnetic energy or co-energy.
! OPT
! 0  —  Stored magnetic energy.
! 1  —  Stored magnetic co-energy.
! ANTYPE
! 0  —  Static or transient.
! 1  —  Harmonic.

SENERGY,1,0
! get co-energy
! 0.12973E-05  in parameter  C_ENG   

! /prep7
CO_energy=C_ENG

! Attention please! When SENERGY macro is recalled, it will clear out table
! value. It is better to save the result into a file
! We also need to have a clear idea, which part of energy or co-energy we need?
! air gap? copper winding? steel? magnets?

SENERGY,0,0
! get energy
!  0.24512E-06   in parameter   S_ENG

Energy=S_ENG

!************************************************************************************
! File operation, write down the results to a file
!************************************************************************************


*cfopen,'test','txt','..\..\Project\Stud\Ansys_model\modified\',APPEND
! *CFOPEN, Fname, Ext, --, Loc
! Opens a "command" file.
!          angle     Iu        Iv         Iw       Luu       Lvv      Lww        Lvu       Lwu      Lwv         torque   energy    coenergy
!*VWRITE,name(1,1),name(1,2),name(1,3),name(1,4),name(1,5),name(1,6),name(1,7),name(1,8),name(1,9),name(1,10),name(1,11),name(1,12),name(1,13)
! %C,    %C,    %C,    %C,     %C,    %C,    %C,     %C,     %C,     %C,     %C     ,%C     ,%C


!*VWRITE,ANGLE1,current_U,current_V,current_W,INDUCTANCE(1,1),INDUCTANCE(2,2),INDUCTANCE(3,3),INDUCTANCE(1,2),INDUCTANCE(1,3),INDUCTANCE(2,3),  !TORQUE,S_ENG,CO_ENERGY
! %G,    %G,    %G,   %G,     %G,    %G,     %G,    %G,    %G,    %G,    %G,    %G,   %G   

!*CFCLOS
! close file


!************************************************************************************
! Remove the force
!************************************************************************************
! BFADELE, AREA, Lab(JS)
! Deletes body force loads on an area.
! BFA,all,JS, , ,-jdens3,0

! use component name, CM_1(phase U),CM_2(phase V),CM_3(phase W)
/PREP7
BFADELE, CM_1, JS

BFADELE, CM_2, JS

BFADELE, CM_3, JS

该命令流是网上一位老兄贴出来的..

infrom

三相磁链之和为什么会为零？

azedd

附件是对三相磁链的一个推导,根据结果三相磁链之和应为0.

zengxiaodong

三相磁链之和应该为零！假设正弦的情况下