为什么算出的磁场不对称？
各位大侠，为什么我算出的永磁磁场波形不对称呀（如图），我的ANSYS程序如下：finish
/clear,start
/PREP7
/units,si
*afun,deg ! Degree units for trig. functions
ET,1,PLANE53
stator_outer_radius=128/2/1000 !定子外半径
stator_pole_arc=22 !定子极弧角度
stator_pole_height=13/1000 !定子齿高（或极高）
stator_inner_radius=75/2/1000 !定子内半径
air_gap=0.45/1000 !气隙宽度
rotor_pole_arc=26 !转子极弧角度
rotor_outer_radius=stator_inner_radius-air_gap !转子外径半
rotor_pole_height=10/1000 !转子齿高（或极高）
rotor_inner_radius=22/2/1000 !转子内半径
rotor_yoke_radius=rotor_outer_radius-rotor_pole_height !计算转子轭半径
stack_length=75/1000 !叠片总厚度
winding_turns=220 !每相绕组匝数
pm_wd=37/1000 !永磁体宽度
pm_height=6/1000 !永磁体高度
magnet_remanence=1.08 !永磁体剩磁（T）
r1=rotor_outer_radius !转子内径，以下将气隙分割为4部分
r2=r1+air_gap/4
r3=r2+air_gap/4
r4=r3+air_gap/4
csys,0
LOCAL,11,0,0,0,0,90, , ,1,1,
csys,0
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,MURX,1,,1!定义第一种材料为气隙
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,MURX,2,,1!定义第二种材料为电枢绕组
!定义第三种材料为硅钢片，非线性B-H曲线，共45个点
TB,BH,3,1,45 !set B-H curve for material 3:stator and rotor iron
TBPT,,138,0.4 !Hc: A/m; Br:T!d23数据用了cm单位
TBPT,,148,0.45
TBPT,,158,0.5
TBPT,,169,0.55
TBPT,,181,0.6
TBPT,,194,0.65
TBPT,,210,0.7
TBPT,,228,0.75
TBPT,,250,0.8
TBPT,,276,0.85
TBPT,,306,0.9
TBPT,,341,0.95
TBPT,,383,1
TBPT,,433,1.05
TBPT,,493,1.1
TBPT,,568,1.15
TBPT,,652,1.2
TBPT,,762,1.25
TBPT,,890,1.3
TBPT,,1050,1.35
TBPT,,1260,1.4
TBPT,,1550,1.45
TBPT,,2010,1.5
TBPT,,2670,1.55
TBPT,,3780,1.6
TBPT,,5340,1.65
TBPT,,7200,1.7
TBPT,,9500,1.75
TBPT,,12200,1.8
TBPT,,15200,1.85
TBPT,,18000,1.89
TBPT,,22000,1.915
TBPT,,27300,1.95
TBPT,,43000,2
TBPT,,65500,2.1
TBPT,,112000,2.15
TBPT,,144000,2.2
TBPT,,184000,2.25
TBPT,,224000,2.3
TBPT,,264000,2.35
TBPT,,304000,2.4
TBPT,,344000,2.45
TBPT,,384000,2.5
TBPT,,424000,2.55 !定义第3种材料为硅钢片，非线性B-H曲线，共45个点
TBPT,,20000000,25.5
MPTEMP,,,,,,,, !定义第4种材料为永磁体，相对磁导率为1.02
MPTEMP,1,0
MPDATA,MURX,4,,1.242
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDE,MGXX,4
MPDE,MGYY,4
MPDE,MGZZ,4
MPDATA,MGxx,4,,756000!永磁矫顽力
!设计转子结构
csys,0
k,1,0,0,0
csys,1
cyl4,0,0,rotor_inner_radius,,rotor_yoke_radius,,,
lsel,s,,,1,8,1
lesize,all,,1
k,10,rotor_outer_radius,rotor_pole_arc/2,0 !rotor_outer_radius=0.03705
csys,0
k,11,sqrt(rotor_yoke_radius**2-ky(10)**2),ky(10),0
ksymm,y,10,11,1,0,0,0
lstr,10,11
lstr,13,12
lsel,s,,,9,10,1
lesize,all,,,30
csys,1
larc,11,13,1,rotor_yoke_radius
larc,10,12,1,rotor_outer_radius
lsel,s,,,11,12,1
lesize,all,,1
lsel,s,,,9,12,1
al,all
asel,s,,,2
agen,6,all,,,,60,,,0
allsel,all
aadd,1,2,3,4,5,6,7
allsel,all
NUMCMP,AREA
NUMCMP,KP
NUMCMP,LINE
!设计气隙
cyl4,0,0,rotor_inner_radius,,r2,,,
cyl4,0,0,rotor_inner_radius,,r3,,,
cyl4,0,0,rotor_inner_radius,,r4,,,
lsel,s,,,31,54,1
lesize,all,,1
AOVLAP,all
allsel,all
NUMCMP,AREA
NUMCMP,KP
NUMCMP,LINE
lsel,s,,,23,30,1
lesize,all,,1
!设计定子结构
!定子齿
k,44,stator_inner_radius,stator_pole_arc/2,0
csys,0
k,45,sqrt(stator_inner_radius**2-ky(44)**2)+stator_pole_height,ky(44),0
ksymm,y,44,45,1,0,0,0
lstr,44,45
lstr,46,47
lsel,s,,,43,44,1
lesize,all,,,35
csys,1
larc,44,46,1,stator_inner_radius
larc,45,47,1,stator_inner_radius+stator_pole_height
lsel,s,,,45,46,1
lesize,all,,1
lsel,s,,,43,46,1
al,all
!asel,s,,,5
agen,2,5,,,,67.5,,,1
adele,5,,,1
agen,2,6,,,,-45,,,1
asel,s,,,5,6,1
agen,4,all,,,,90,,,0
!定子轭
csys,0
k,76,sqrt((stator_inner_radius+stator_pole_height)**2-(pm_height/2)**2),pm_height/2,0
k,77,sqrt((stator_inner_radius+stator_pole_height)**2-(pm_height/2)**2)+pm_wd,pm_height/2,0
ksymm,y,76,77,1,0,0,0
k,80,sqrt((stator_inner_radius+stator_pole_height)**2-(pm_height/2)**2),0,0
lstr,76,77
lstr,78,79
lsel,s,,,75,76,1
lesize,all,,,50
lstr,76,80
lstr,80,78
lsel,s,,,77,78,1
lesize,all,,,5
lstr,77,79
lesize,79,,,10
lsel,s,,,75,79,1
al,all
csys,1
agen,2,13,,,,180,,,0
csys,0
theta=acos(stator_outer_radius/sqrt((sqrt((stator_inner_radius+stator_pole_height)**2-(pm_height/2)**2)+pm_wd)**2+(pm_height/2)**2))
k,86,stator_outer_radius*cos(theta),stator_outer_radius*sin(theta),0
ksymm,y,86,,,0,0,0
ksymm,x,86,87,1,0,0,0
lstr,77,86
lstr,83,88
lstr,82,89
lstr,79,87
lsel,s,,,85,88,1
lesize,all,,,70
csys,1
larc,86,88,1,stator_outer_radius
larc,89,87,1,stator_outer_radius
lsel,s,,,89,90,1
lesize,all,,1,
k,90,stator_inner_radius,0,0
k,91,stator_inner_radius,45,0
k,92,stator_inner_radius,90,0
k,93,stator_inner_radius,135,0
k,94,stator_inner_radius,180,0
k,95,stator_inner_radius,225,0
k,96,stator_inner_radius,270,0
k,97,stator_inner_radius,315,0
k,98,stator_inner_radius+stator_pole_height,45,0
k,99,stator_inner_radius+stator_pole_height,90,0
k,100,stator_inner_radius+stator_pole_height,135,0
k,101,stator_inner_radius+stator_pole_height,225,0
k,102,stator_inner_radius+stator_pole_height,270,0
k,103,stator_inner_radius+stator_pole_height,315,0
lstr,80,90
lstr,91,98
lstr,92,99
lstr,93,100
lstr,94,85
lstr,95,101
lstr,96,102
lstr,97,103
lsel,s,,,91,98,1
lesize,all,,,30
larc,46,76,1,stator_inner_radius+stator_pole_height
larc,45,98,1,stator_inner_radius+stator_pole_height
larc,98,50,1,stator_inner_radius+stator_pole_height
larc,49,99,1,stator_inner_radius+stator_pole_height
larc,99,54,1,stator_inner_radius+stator_pole_height
larc,53,100,1,stator_inner_radius+stator_pole_height
larc,100,58,1,stator_inner_radius+stator_pole_height
larc,57,84,1,stator_inner_radius+stator_pole_height
larc,81,62,1,stator_inner_radius+stator_pole_height
larc,61,101,1,stator_inner_radius+stator_pole_height
larc,101,66,1,stator_inner_radius+stator_pole_height
larc,65,102,1,stator_inner_radius+stator_pole_height
larc,102,70,1,stator_inner_radius+stator_pole_height
larc,69,103,1,stator_inner_radius+stator_pole_height
larc,103,74,1,stator_inner_radius+stator_pole_height
larc,73,78,1,stator_inner_radius+stator_pole_height
lsel,s,,,99,114,1
lesize,all,,1,
lsel,s,,,99,106,1
lsel,a,,,44,56,4
lsel,a,,,85,86,1
lsel,a,,,75,89,7
al,all
lsel,s,,,107,114,1
lsel,a,,,60,76,4
lsel,a,,,87,88,1
lsel,a,,,80,90,10
al,all
larc,90,47,1,stator_inner_radius
larc,44,91,1,stator_inner_radius
larc,91,51,1,stator_inner_radius
larc,48,92,1,stator_inner_radius
larc,92,55,1,stator_inner_radius
larc,52,93,1,stator_inner_radius
larc,93,59,1,stator_inner_radius
larc,56,94,1,stator_inner_radius
larc,94,63,1,stator_inner_radius
larc,60,95,1,stator_inner_radius
larc,95,67,1,stator_inner_radius
larc,64,96,1,stator_inner_radius
larc,96,71,1,stator_inner_radius
larc,68,97,1,stator_inner_radius
larc,97,75,1,stator_inner_radius
larc,72,90,1,stator_inner_radius
lsel,s,,,115,130,1
lesize,all,,1,
lsel,s,,,115,130,1
lsel,a,,,46,74,4
al,all
allsel,all
AOVLAP,all
adele,18,
NUMCMP,AREA
NUMCMP,KP
NUMCMP,LINE
al,45,77,91,99,115
al,43,92,100,116
al,49,92,101,117
al,47,93,102,118
al,53,93,103,119
al,51,94,104,120
al,57,94,105,121
al,55,83,95,106,122
al,61,84,95,107,123
al,59,96,108,124
al,65,96,109,125
al,63,97,110,126
al,69,97,111,127
al,67,98,112,128
al,73,98,113,129
al,71,78,91,114,130
cyl4,0,0,rotor_inner_radius,,sqrt((sqrt((stator_inner_radius+stator_pole_height)**2-(pm_height/2)**2)+pm_wd)**2+(pm_height/2)**2)+0.01,,,
allsel,all
AOVLAP,all
lsel,s,,,131,134,1
lesize,all,,1,
aadd,5,6,7,8,15
aadd,9,10,11,12,16
allsel,all
NUMCMP,AREA
NUMCMP,KP
NUMCMP,LINE
!定义区域材料属性
asel,s,,,2,4,1 !select air gap
asel,a,,,8,26,18!select air space
aatt,1,1 !assign material 1 to air areas
asel,s,,,9,24,1 !select slots
aatt,2,2 !assign material 2 to the areas of coils
asel,s,,,5,25,20!select stator yoke,select stator tooth
asel,a,,,1 !select rotor yoke，select rotor tooth
aatt,3,3 !assign material 3 to the iron of stator anf rotor
asel,s,,,6,7,1 !select 2 PMs and assign the material 4
aatt,4,11,,11
!施加边界条件
lsel,s,,,123,126,1
DL,all, ,AZ,0,
lsel,s,,,1,4,1
DL,all, ,AZ,0, !apply az=0 to the selected lines
allsel,all
!划分网格
SMRT,6 !intelligent meshing
MSHAPE,1,2D
MSHKEY,0
amesh,all
/SOL
FINISH
/POST1
magsolv
/gst,on
FINISH
/POST1
senergy
senergy,1
FINISH
/POST1
plf2d
大侠看看这些程序错在哪里呀？ 直接上波形 直接上波形？什么意思？ 回复 1# sunqiang960424
老兄，你的图呢？M
PTEMP,,,,,,,,
MPTEMP,1,0
MPDE,MGXX,4
MPDE,MGYY,4
MPDE,MGZZ,4 ？？是什么意思？
磁极充磁方向好像没弄对，是不？用局部坐标系1 我也有过这个问题，最后发现两种永磁体，磁化方向我只改了一种的方向，另一种默认方向x了，所以不对称。
页:
[1]