Ansys 中用Lmatrix，计算inductance 的问题,特别带有一个brushless motor 例子

acer007

要做一个brushless motor的simulation，平台Ansys. 其中，coil 是27 turns，所有的数据都按照datasheet输入，但是在计算self inductance 的时候，结果不能和datasheet符合，计算结果45 mH, 而在datasheet上，为45 uH,结果差了1000倍啊，究竟哪里的错误？
各位高人可否指点一下？谢谢了！

coil 是27 turns，所以在计算 current density 使用了

27*current(single turnm)/size of coil(58 mm^2)    coil的面积为58mm^2



! Summary: This is a model without loop structure. It works well.

! 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








FINISH

perch88

注释太多了，不方便看
材料换线性的看差多少，现在是非线性，如果你的电流由于换算成1极、串并联的时候换算错了，电流加的跟实际相差太大。导致工作点跟实际的不一样，也可能是个原因。

koyan725

acer007你好，谢谢你分享这个例子，我最近对BLDC motor 也很感兴趣。能问问你上面程序是怎么计算出来self inductance=45mH吗？
上面的例子是2D模型，程序计算出来的self inductance=1.85H，你是怎么转换到3D中的self inductance=45mH的？直接乘以线圈的Z向长度？

zhendfliuzp

我也有楼上遇到的问题，过了这么久了，还有没有人能解决这个问题？

sunxikai007

没有乘以轴向长度吧。。。

zengxiaodong

是要乘以轴向长度的，我这里有个BLDC的微分电感，随着转子转动360度电角度的变化波形图：

自感：


互感：