直流电机例子

/PREP7
ET,1,PLANE53 !设置单元类型 plane53
/UNITS,SI !采用国际单位制
EMUNIT,MKS !采用公斤米制
LOCAL,11,0 !定义局部卡式坐标
MP,MURX,1,1 !定义空气材料属性
/com D21-0.5 !注释
TB,BH,2,,170 !定义D21-0.5硅钢片磁化曲线 15个点
TBPT,,4000,1.4
TBPT,,4100,1.43
TBPT,,4200,1.46
TBPT,,4300,1.49
TBPT,,4400,1.52
TBPT,,4500,1.55
TBPT,,4600,1.58
TBPT,,4700,1.61
TBPT,,4800,1.64
TBPT,,4900,1.67

TBPT,,5000,1.7
TBPT,,5100,1.73
TBPT,,5200,1.76
TBPT,,5300,1.79
TBPT,,5400,1.82
TBPT,,5500,1.85
TBPT,,5600,1.88
TBPT,,5700,1.91
TBPT,,5800,1.94
TBPT,,5900,1.97

TBPT,,6000,2
TBPT,,6100,2.03
TBPT,,6200,2.06
TBPT,,6300,2.09
TBPT,,6400,2.12
TBPT,,6500,2.15
TBPT,,6600,2.18
TBPT,,6700,2.21
TBPT,,6800,2.24
TBPT,,6900,2.27

TBPT,,7000,2.3
TBPT,,7100,2.33
TBPT,,7200,2.36
TBPT,,7300,2.39
TBPT,,7400,2.42
TBPT,,7500,2.45
TBPT,,7600,2.48
TBPT,,7700,2.51
TBPT,,7800,2.54
TBPT,,7900,2.57

TBPT,,8000,2.6
TBPT,,8100,2.64
TBPT,,8200,2.68
TBPT,,8300,2.72
TBPT,,8400,2.76
TBPT,,8500,2.8
TBPT,,8600,2.84
TBPT,,8700,2.88
TBPT,,8800,2.92
TBPT,,8900,2.96

TBPT,,9000,3
TBPT,,9100,3.04
TBPT,,9200,3.08
TBPT,,9300,3.12
TBPT,,9400,3.16
TBPT,,9500,3.2
TBPT,,9600,3.24
TBPT,,9700,3.28
TBPT,,9800,3.32
TBPT,,9900,3.36

TBPT,,10000,3.4
TBPT,,10100,3.44
TBPT,,10200,3.48
TBPT,,10300,3.52
TBPT,,10400,3.56
TBPT,,10500,3.6
TBPT,,10600,3.64
TBPT,,10700,3.68
TBPT,,10800,3.72
TBPT,,10900,3.76

TBPT,,11000,3.8
TBPT,,11100,3.84
TBPT,,11200,3.88
TBPT,,11300,3.92
TBPT,,11400,3.96
TBPT,,11500,4
TBPT,,11600,4.05
TBPT,,11700,4.1
TBPT,,11800,4.15
TBPT,,11900,4.2

TBPT,,12000,4.25
TBPT,,12100,4.3
TBPT,,12200,4.35
TBPT,,12300,4.4
TBPT,,12400,4.45
TBPT,,12500,4.5
TBPT,,12600,4.6
TBPT,,12700,4.7
TBPT,,12800,4.8
TBPT,,12900,4.9

TBPT,,13000,5
TBPT,,13100,5.2
TBPT,,13200,5.4
TBPT,,13300,5.6
TBPT,,13400,5.8
TBPT,,13500,6
TBPT,,13600,6.6
TBPT,,13700,7.2
TBPT,,13800,7.8
TBPT,,13900,8.4

TBPT,,14000,9
TBPT,,14100,10
TBPT,,14200,11
TBPT,,14300,12
TBPT,,14400,13
TBPT,,14500,14
TBPT,,14600,15
TBPT,,14700,16
TBPT,,14800,17
TBPT,,14500,18.5

TBPT,,15000,20
TBPT,,15100,21.5
TBPT,,15200,23
TBPT,,15300,24.5
TBPT,,15400,26
TBPT,,15500,28
TBPT,,15600,30
TBPT,,15700,32
TBPT,,15800,34
TBPT,,15000,36

TBPT,,16000,39
TBPT,,16100,42
TBPT,,16200,45
TBPT,,16300,48
TBPT,,16400,51
TBPT,,16500,54
TBPT,,16600,57.5
TBPT,,16700,61
TBPT,,16800,64.5
TBPT,,16900,68

TBPT,,17000,71.5
TBPT,,17100,75
TBPT,,17200,79
TBPT,,17300,83
TBPT,,17400,87
TBPT,,17500,91
TBPT,,17600,95
TBPT,,17700,99
TBPT,,17800,103
TBPT,,17900,108

TBPT,,18000,113
TBPT,,18100,118
TBPT,,18200,123
TBPT,,18300,128
TBPT,,18400,133
TBPT,,18500,138
TBPT,,18600,143
TBPT,,18700,148
TBPT,,18800,154
TBPT,,18900,160

TBPT,,19000,166
TBPT,,19100,172
TBPT,,19200,179
TBPT,,19300,186
TBPT,,19400,193
TBPT,,19500,200
TBPT,,19600,207
TBPT,,19700,214
TBPT,,19800,221
TBPT,,19900,228

TBPT,,20000,236
TBPT,,20100,244
TBPT,,20200,254
TBPT,,20300,265
TBPT,,20400,276
TBPT,,20500,287
TBPT,,20600,298
TBPT,,20700,311
TBPT,,20800,324
TBPT,,20900,345
!TBPLOT,BH,2

!输入参数
THETA1=0
THETA2=180
THETA3=11.25 !调整转子位置的角度
HC=292470 !待定，钕铁硼的HC随温度升高而升高
S=0.001 !单位换算因子（mm->M)
R1=7.75*S !轴直径
R2=26.5*S
R3=27.5*S
R4=30.75*S !磁瓦内径，要偏心////
R5=36.6*S !磁瓦外径
D1=13.305*S
*AFUN,DEG !转换成角度
MP,MURX,3,1.15 !定义磁的相对磁导率
MP,MGXX,3,-HC*COS(THETA1) !磁钢线密
MP,MGYY,3,-HC*SIN(THETA1)

*AFUN,DEG
MP,MURX,4,1.15
MP,MGXX,4,HC*COS(THETA2)
MP,MGYY,4,HC*SIN(THETA2)

/PNUM,AREA,1 !创建模型并标注面编号及将面着色
LOCAL,11,1 !定义局部圆柱坐标
K,1,0,0+THETA3
!定义关键点
K,2,R1,11.25+THETA3
K,3,D1,11.25+THETA3
K,4,14.52*S,4.9403+THETA3
K,5,25.21*S,2.8419+THETA3
K,6,R2,5.5571+THETA3
K,7,R2,5.9372+THETA3
K,8,26.96*S,6.9954+THETA3
K,9,R3,7.0794+THETA3

K,10,R3,-7.0794+THETA3
K,11,26.96*S,-6.9954+THETA3
K,12,R2,-5.9372+THETA3
K,13,R2,-5.5571+THETA3
K,14,25.21*S,-2.8419+THETA3
K,15,14.52*S,-4.9403+THETA3
K,16,D1,-11.25+THETA3
K,17,R1,-11.25+THETA3

!磁瓦关键点
K,18,29.51*S,-62.2974+90
K,19,31.2*S,-66.788+90
K,20,35.31*S,-70+90
K,21,36.5*S,-65+90
K,22,36.5*S,65+90
K,23,35.31*S,70+90
K,24,31.2*S,66.788+90
K,25,29.51*S,62.2974+90

!圆心
K,26,15*S,11.25+THETA3 !3，4圆心
K,27,25.3*S,5.5571+THETA3 !5，6圆心
K,28,27*S,5.9372+THETA3 !7，8圆心
K,29,27*S,-5.9372+THETA3 !11，12圆心
K,30,25.3*S,-5.5571+THETA3 !13，14圆心
K,31,15*S,-11.25+THETA3 !15，16圆心
K,32,5*S,-90 !18，25圆心
K,33,32.5*S,-61.9167+90 !18，19圆心
K,34,32.5*S,61.9167+90 !24，25圆心

LSTR,1,2 !连接1，2关键点生成直线，第一条直线
LSTR,2,3 !L2
LARC,3,4,26,1.695*S !L3，三点和一半径生成弧线
LSTR,4,15 !L4
LSTR,4,5 !L5
LARC,5,6,27,1.2*S !L6
LARC,6,7,1,R2 !L7
LARC,7,8,28,0.5*S !L8
LSTR,8,9 !L9
LARC,9,10,1,R3 !L10
LSTR,10,11 !L11
LARC,11,12,29,0.5*S !L12
LARC,12,13,1,R2 !L13
LARC,13,14,30,1.2*S !L14
LSTR,14,15 !L15
LARC,15,16,31,1.695*S !L16
LSTR,16,17 !L17
LARC,2,17,1,R1 !L18
LSTR,1,17 !L19

!18，19，20，21，22，23，24，25条线构成磁瓦线条
LARC,18,19,33,3*S !L20
LSTR,19,20 !L21
LSTR,20,21 !L22
LARC,21,22,1,R5 !L23
LSTR,22,23 !L24
LSTR,23,24 !L25
LARC,24,25,34,3*S !L26
LARC,18,25,32,R4 !L27

AL,1,18,19 !第1，16，17条线生成面积1
AL,2,3,4,16,17,18 !面积2
AL,4,5,6,7,8,9,10,11,12,13,14,15 !面积3
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!1
/PREP7
ET,1,PLANE53 ! plane53
/UNITS,SI !
EMUNIT,MKS
*SET,THETA1,0
*SET,THETA2,180
*SET,THETA3,11.25 !
*SET,HC,292470 ! HC
*SET,S,0.001 ! mm->M)
*SET,R1,7.75*S !
*SET,R2,26.5*S
*SET,R3,27.5*S
*SET,R4,30.75*S ! ////
*SET,R5,36.6*S !
*SET,D1,13.305*S
*AFUN,DEG !
MP,MURX,3,1.15 !
MP,MGXX,3,-HC*COS(THETA1) !
MP,MGYY,3,-HC*SIN(THETA1)
*AFUN,DEG
MP,MURX,4,1.15
MP,MGXX,4,HC*COS(THETA2)
MP,MGYY,4,HC*SIN(THETA2)
/PNUM,AREA,1 !
LOCAL,11,1 !
K,1,0,0+THETA3
!
K,2,R1,11.25+THETA3
K,3,D1,11.25+THETA3
K,4,14.52*S,4.9403+THETA3
K,5,25.21*S,2.8419+THETA3
K,6,R2,5.5571+THETA3
K,7,R2,5.9372+THETA3
K,8,26.96*S,6.9954+THETA3
K,9,R3,7.0794+THETA3
K,10,R3,-7.0794+THETA3
K,11,26.96*S,-6.9954+THETA3
K,12,R2,-5.9372+THETA3
K,13,R2,-5.5571+THETA3
K,14,25.21*S,-2.8419+THETA3
K,15,14.52*S,-4.9403+THETA3
K,16,D1,-11.25+THETA3
K,17,R1,-11.25+THETA3
!
K,18,29.51*S,-62.2974+90
K,19,31.2*S,-66.788+90
K,20,35.31*S,-70+90
K,21,36.5*S,-65+90
K,22,36.5*S,65+90
K,23,35.31*S,70+90
K,24,31.2*S,66.788+90
K,25,29.51*S,62.2974+90
!
K,26,15*S,11.25+THETA3 !3 4
K,27,25.3*S,5.5571+THETA3 !5 6
K,28,27*S,5.9372+THETA3 !7 8
K,29,27*S,-5.9372+THETA3 !11 12
K,30,25.3*S,-5.5571+THETA3 !13 14
K,31,15*S,-11.25+THETA3 !15 16
K,32,5*S,-90 !18 25
K,33,32.5*S,-61.9167+90 !18 19
K,34,32.5*S,61.9167+90 !24 25
LSTR,1,2 ! 1 2
LSTR,2,3 !L2
LARC,3,4,26,1.695*S !L3
LSTR,4,15 !L4
LSTR,4,5 !L5
LARC,5,6,27,1.2*S !L6
LARC,6,7,1,R2 !L7
LARC,7,8,28,0.5*S !L8
LSTR,8,9 !L9
LARC,9,10,1,R3 !L10
LSTR,10,11 !L11
LARC,11,12,29,0.5*S !L12
LARC,12,13,1,R2 !L13
LARC,13,14,30,1.2*S !L14
LSTR,14,15 !L15
LARC,15,16,31,1.695*S !L16
LSTR,16,17 !L17
LARC,2,17,1,R1 !L18
LSTR,1,17 !L19
!18 19 20 21 22 23 24 25
LARC,18,19,33,3*S !L20
LSTR,19,20 !L21
LSTR,20,21 !L22
LARC,21,22,1,R5 !L23
LSTR,22,23 !L24
LSTR,23,24 !L25
LARC,24,25,34,3*S !L26
LARC,18,25,32,R4 !L27

!*

!*
ALLSEL,ALL
LPLOT

KL, 4, 0.000000000000000

KL, 4, 1.000000000000000
LSTR, 35, 36
FLST,2,2,4,ORDE,2
FITEM,2,5
FITEM,2,15
LDELE,P51X, , ,1
LSTR, 35, 5
LSTR, 36, 14
AL,1,18,19 ! 1 16 17 1
AL,2,3,4,16,17,18 ! 2

/REP,FAST
FLST,2,12,4
FITEM,2,28
FITEM,2,5
FITEM,2,6
FITEM,2,7
FITEM,2,8
FITEM,2,9
FITEM,2,10
FITEM,2,11
FITEM,2,12
FITEM,2,13
FITEM,2,14
FITEM,2,15
AL,P51X
APLOT