- 积分
- 0
- 注册时间
- 2006-10-22
- 仿真币
-
- 最后登录
- 1970-1-1
|
要做一个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 |
评分
-
1
查看全部评分
-
|