西莫电机圈

 找回密码
 立即注册

QQ登录

只需一步,快速开始

手机号码,快捷登录

手机号码,快捷登录

查看: 1689|回复: 4

为什么算出的磁场不对称?

[复制链接]

该用户从未签到

发表于 2011-1-26 11:45 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?立即注册

x
各位大侠,为什么我算出的永磁磁场波形不对称呀(如图),我的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
大侠看看这些程序错在哪里呀?
西莫电机论坛微信公众平台正式上线!★详情请点击★ 西莫电机论坛会员交流专用群欢迎您西莫电机论坛加群请注明论坛用户名及所从事专业,否则不予通过

该用户从未签到

发表于 2011-1-26 13:35 | 显示全部楼层
直接上波形
西莫电机论坛微信公众平台正式上线!★详情请点击★ 西莫电机论坛会员交流专用群欢迎您西莫电机论坛加群请注明论坛用户名及所从事专业,否则不予通过
回复

使用道具 举报

该用户从未签到

 楼主| 发表于 2011-2-1 14:59 | 显示全部楼层
直接上波形?什么意思?
西莫电机论坛微信公众平台正式上线!★详情请点击★ 西莫电机论坛会员交流专用群欢迎您西莫电机论坛加群请注明论坛用户名及所从事专业,否则不予通过
回复

使用道具 举报

该用户从未签到

发表于 2011-3-4 20:49 | 显示全部楼层
回复 1# sunqiang960424
老兄,你的图呢?M
PTEMP,,,,,,,,  
MPTEMP,1,0  
MPDE,MGXX,4
MPDE,MGYY,4                  
MPDE,MGZZ,4 ??是什么意思?
磁极充磁方向好像没弄对,是不?用局部坐标系1
西莫电机论坛微信公众平台正式上线!★详情请点击★ 西莫电机论坛会员交流专用群欢迎您西莫电机论坛加群请注明论坛用户名及所从事专业,否则不予通过
回复

使用道具 举报

该用户从未签到

发表于 2017-10-12 17:03 | 显示全部楼层
我也有过这个问题,最后发现两种永磁体,磁化方向我只改了一种的方向,另一种默认方向x了,所以不对称。
西莫电机论坛微信公众平台正式上线!★详情请点击★ 西莫电机论坛会员交流专用群欢迎您西莫电机论坛加群请注明论坛用户名及所从事专业,否则不予通过
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

QQ|小黑屋|手机版|Archiver|西莫电机圈 ( 浙ICP备10025899号-3 浙公网安备:33028202000436号

GMT+8, 2024-11-22 16:23 , Processed in 0.080438 second(s), 21 queries .

Powered by Discuz! X3.4

© 2001-2023 Discuz! Team.

快速回复 返回顶部 返回列表