CN105373014B - 计及边缘扩散效应的bsrm建模方法 - Google Patents
计及边缘扩散效应的bsrm建模方法 Download PDFInfo
- Publication number
- CN105373014B CN105373014B CN201510916257.6A CN201510916257A CN105373014B CN 105373014 B CN105373014 B CN 105373014B CN 201510916257 A CN201510916257 A CN 201510916257A CN 105373014 B CN105373014 B CN 105373014B
- Authority
- CN
- China
- Prior art keywords
- air gap
- path
- rotor
- theta
- phase
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B17/00—Systems involving the use of models or simulators of said systems
- G05B17/02—Systems involving the use of models or simulators of said systems electric
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Iron Core Of Rotating Electric Machines (AREA)
Abstract
本发明公开了一种计及边缘扩散效应的BSRM建模方法,包括:根据BSRM气隙的定子和转子活动区域,分别选取转子齿极偏离定子齿极角度θ=0°和θ≠0°时的,计及边缘扩散效应的麦克斯韦应力法积分路径;根据所述麦克斯韦应力法积分路径,计算θ=0°和θ≠0°时的气隙磁密;根据所述气隙磁密计算得到转子齿极所受α、β方向悬浮力和电磁转矩;建立BSRM的α、β方向悬浮力和电磁转矩解析模型。本发明实现了BSRM解析建模,计及边缘扩散效应,建模精度高。
Description
技术领域
本发明涉及无轴承开关磁阻电机及其控制领域,特别涉及一种计及边缘扩散效应的BSRM建模方法。
背景技术
无轴承开关磁阻电机(bearingless switched reluctance motor,BSRM)是发展迅速的磁悬浮技术与开关磁阻电机(switched reluctance motor,SRM)的结合,兼有结构简单坚固、成本低、调速范围宽、运行可靠性高和允许转速高、摩擦功耗小、无需润滑和寿命长等优点,在高速、超高速运行场合具有突出优势,是高速电机研究领域的热点之一。
实际电机中电磁场的磁路分布极为复杂,磁通并非完全限制在气隙磁路中,气隙内磁通分布和磁力线分布也并非完全均匀,在气隙边缘部分有向外扩散的趋势,使得气隙中磁力线产生一定的弯曲变形,将此现象称之为“气隙边缘效应”。
BSRM将主绕组和悬浮绕组叠绕在定子极上,利用主绕组电流和悬浮绕组电流相互作用形成径向力实现转子悬浮。与普通SRM相比,BSRM不仅存在因磁路饱和引起的非线性,而且电磁转矩与径向悬浮力之间、α和β方向的径向悬浮力之间均存在非线性的耦合关系,进而导致BSRM精确建模困难。
发明内容
本发明所要解决的技术问题在于,提供一种计及边缘扩散效应的BSRM建模方法,实现计及边缘扩散效应的BSRM解析建模,建模精度高。
为了解决上述技术问题,本发明提供了一种计及边缘扩散效应的BSRM建模方法,包括:
根据BSRM气隙的定子和转子活动区域,分别选取转子齿极偏离定子齿极角度θ=0°和θ≠0°时的,计及边缘扩散效应的麦克斯韦应力法积分路径;
根据所述麦克斯韦应力法积分路径,计算θ=0°和θ≠0°时的气隙磁密;
根据所述气隙磁密计算得到转子齿极所受α、β方向悬浮力和电磁转矩;
建立BSRM的α、β方向悬浮力和电磁转矩解析模型。
实施本发明,具有如下有益效果:本发明实现了BSRM精确建模,建模精度高。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是三相12/8结构BSRM的结构示意图;
图2是本发明提供的计及边缘扩散效应的BSRM建模方法的一个实施例的流程示意图;
图3是θ=0°时气隙a1处定、转子相对位置示意图;
图4是θ<0°时气隙a1处定、转子相对位置示意图;
图5是θ>0°时气隙a1处定、转子相对位置示意图;
图6所示为α方向悬浮力Fα和电磁转矩与虚位移法、传统麦克斯韦应力法和有限元仿真结果对比图;
图7所示为β方向悬浮力Fβ和电磁转矩与虚位移法、传统麦克斯韦应力法和有限元仿真结果对比图;
图8所示为电磁转矩T和电磁转矩与虚位移法、传统麦克斯韦应力法和有限元仿真结果对比图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示,本发明实施例研究对象为三相12/8结构BSRM,图中仅画出了A相的主绕组(Nma)、α方向悬浮绕组(Nsa1)和β方向悬浮绕组(Nsa2)。Nma由相隔90°的4个极绕组正串而成,而Nsa1、Nsa2分别由α、β方向径向相对的2个悬浮绕组反串而成。Nma产生偏置磁通,Nsa1分别对此偏置磁通产生增强(气隙a1处)、削弱(气隙a3处)的作用,由此产生不平衡磁拉力。同理,通过Nsa2对Nma偏置磁场在气隙a2处的增强及在气隙a4处的削弱作用,也产生不平衡磁拉力。对磁拉力进行分解,可得到α、β方向悬浮力Fα、Fβ。
B相和C相具有和A相相同的绕组结构、连接方式和悬浮机理。利用三相绕组每隔15°轮流导通和转子位移的负反馈控制,可实现转轴稳定悬浮。
如图2所示,本发明实施例具体包括步骤:
S101、根据BSRM气隙的定子和转子活动区域,分别选取转子齿极偏离定子齿极角度θ=0°和θ≠0°时的,计及边缘扩散效应的麦克斯韦应力法积分路径。
具体的,步骤S101包括步骤:
S1011、θ=0°时,从定子和转子的交叠活动区域的边界线上选取点1,2,3,4,路径1→4为麦克斯韦应力法积分路径。具体路径如图3所示。图3所示为θ=0°气隙a1处定、转子相对位置示意图。
S1012、θ≠0°时,即θ<0°或θ>0°时,从BSRM气隙的定子和转子的非交叠活动区域的转子边界线上选取点1、2、3、4、10、11,其中,点4是边缘气隙与主气隙在转子边界线上的分界点,点2及点11分别是两个边缘气隙边界线与转子边界线的交点,点3及点10分别是两个边缘气隙平均磁密线与转子边界线的交点,边缘气隙平均磁密线将两个边缘气隙区域分别划分为f1区域和f2区域;从定子和转子的交叠活动区域的主气隙中心线上选取点5、6,其中,点5是主气隙中心线与转子边界线上的交点,其中,点6是主气隙中心线与定子边界线上的交点;从定子和转子的非交叠活动区域的定子边界线上选取点7、8,其中,点7是边缘气隙与主气隙在定子边界线上的分界点,点8是边缘气隙平均磁密线与定子边界线的交点;边缘气隙边界线近似为主气隙长度的直线加1/4圆弧轨迹线,从边缘气隙平均磁密线上选取两段轨迹线的连接点9;选取路径1→11为麦克斯韦应力法积分路径。具体路径如图4、图5所示,其中图4所示为θ<0°气隙a1处定、转子相对位置示意图;图5所示为θ>0°气隙a1处定、转子相对位置示意图。
S102、根据所述麦克斯韦应力法积分路径,计算θ=0°和θ≠0°时的气隙磁密。
具体的,步骤S102包括步骤:
S1021、θ=0°时,根据路径1→4,计算得到A相气隙的各气隙磁密。
其中,根据路径1→4计算得到的A相气隙aj在路径2→3处的主气隙平均磁密为
式中,η为气隙边缘系数,μ0为真空磁导率,l0为定、转子间气隙长度,气隙a1~a4处的磁路磁动势为Nm为主绕组匝数,Ns为悬浮绕组匝数,ima为A相主绕组电流,isa1为A相α方向悬浮绕组电流,isa2为A相β方向悬浮绕组电流。
其中,边缘磁力线弯曲变形区域的磁路长度数据,在悬浮力建模中,考虑引入气隙边缘系数η,用于对磁路平均长度进行归算:即采用有限元分析软件Maxwell2D,在θ=0°时,改变电机转子半径r、气隙长度l0、定子极弧比和转子极弧比,仿真计算电机转子所受α方向悬浮力Fα,得到气隙边缘系数η值随气隙长度l0和定、转子齿极弧度比τs/τr变化有一定的变化,但变化范围较小。因此,本发明专利选取η值为一恒定值1.02。
S1022、θ≠0°时,根据路径1→11,计算得到A相气隙的各气隙磁密。
其中,计算得到的A相气隙aj在路径5→6处的主气隙磁密为Bmaj=μ0Uaj/l0,j=1,2,3,4;
计算得到的A相气隙aj在路径4→5和路径6→7处的主气隙平均磁密为
计算得到的A相气隙aj在路径8→9和路径9→10处的边缘气隙平均磁密为式中,r为转子半径;
计算得到的A相气隙aj在路径2→3和路径10→11处的边缘气隙平均磁密为式中,为路径2→3和路径10→11所在的边缘气隙f1区域的气隙磁路平均长度,
计算得到的A相气隙aj在路径3→4和路径7→8处的边缘气隙平均磁密为式中,为路径3→4和路径7→8所在的边缘气隙f2区域的气隙磁路平均长度。
S103、根据所述气隙磁密计算得到转子齿极所受α、β方向悬浮力和电磁转矩。
具体的,步骤S103包括步骤:
S1031、根据计算到的θ=0°时的气隙磁密,计算得到A相气隙a1~a4处转子齿极所受α方向悬浮力Fα1~Fα4、所受β方向悬浮力Fβ1~Fβ4、电磁转矩T1~T4分别为:
式中,h为转子叠片长度,τr为转子齿极弧度。
其中,A相气隙a1~a4处转子齿极所受α方向悬浮力Fα1~Fα4、所受β方向悬浮力Fβ1~Fβ4、电磁转矩T1~T4的计算过程为:
Tj=0,j=1,2,3,4
整理后得到Fα1~Fα4、Fβ1~Fβ4、T1~T4。
S1032、根据计算到的θ<0°时的气隙磁密,计算得到A相气隙a1~a4处转子齿极所受α方向悬浮力Fα1~Fα4、所受β方向悬浮力Fβ1~Fβ4、电磁转矩T1~T4分别为:
式中,γ为积分路径法向量与转子齿极中心线方向的夹角,
S1033、整理后得到
式中,
S1034、根据θ<0°时的A相气隙a1~a4处转子齿极所受α方向悬浮力、所受β方向悬浮力和电磁转矩,同理计算得到θ>0°时的A相气隙a1~a4处转子齿极所受α方向悬浮力Fα1~Fα4、所受β方向悬浮力Fβ1~Fβ4、电磁转矩T1~T4分别为
S104、建立BSRM的α、β方向悬浮力和电磁转矩解析模型。
具体的,步骤S104包括步骤:
S1041、当BSRM在线性工作状态下且仅A相通电时,α方向悬浮力Fα、β方向悬浮力Fβ、电磁转矩T分别为:
S1042、根据θ=0°、θ<0°、θ>0°时的A相气隙a1~a4处转子齿极所受α方向悬浮力、所受β方向悬浮力和电磁转矩,对上式转换整理得到BSRM的α、β方向悬浮力和电磁转矩解析模型为
式中,im为通电相主绕组电流;is1为通电相α方向悬浮绕组电流;is2为通电相β方向悬浮绕组电流;悬浮力比例系数Kα、Kβ和电磁转矩系数Jt分别为:
下面对本发明实施例进行仿真分析,本实施例中所采用BSRM的相关参数如表1所示。
表1仿真分析所用BSRM参数
将本实施例所建BSRM线性解析模型中的α、β方向上的悬浮力和电磁转矩与虚位移法、传统麦克斯韦应力法、有限元仿真结果对比,分别如图5~图7中所示。
图6所示为α方向悬浮力Fα和电磁转矩与虚位移法、传统麦克斯韦应力法和有限元仿真结果对比图;
图7所示为β方向悬浮力Fβ和电磁转矩与虚位移法、传统麦克斯韦应力法和有限元仿真结果对比图;
图8所示为电磁转矩T和电磁转矩与虚位移法、传统麦克斯韦应力法和有限元仿真结果对比图。
可以看出,本发明的建模精度较高。
实施本发明,具有如下有益效果:本发明实现了计及边缘扩散效应的BSRM解析建模,建模精度高。
需要说明的是,在本文中,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者装置不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者装置所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括该要素的过程、方法、物品或者装置中还存在另外的相同要素。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
Claims (6)
1.一种计及边缘扩散效应的BSRM建模方法,其特征在于,包括:
根据BSRM气隙的定子和转子活动区域,分别选取转子齿极偏离定子齿极角度θ=0°和θ≠0°时的,计及边缘扩散效应的麦克斯韦应力法积分路径;
根据所述麦克斯韦应力法积分路径,计算θ=0°和θ≠0°时的气隙磁密;
根据所述气隙磁密计算得到转子齿极所受α、β方向悬浮力和电磁转矩;
建立BSRM的α、β方向悬浮力和电磁转矩解析模型;
所述根据BSRM气隙的定子和转子活动区域,分别选取转子齿极偏离定子齿极角度θ=0°和θ≠0°时的,计及边缘扩散效应的麦克斯韦应力法积分路径,具体包括:
θ=0°时,从定子和转子的交叠活动区域的边界线上选取点1,2,3,4,路径1→4为麦克斯韦应力法积分路径;
θ≠0°时,即θ<0°或θ>0°时,从BSRM气隙的定子和转子的非交叠活动区域的转子边界线上选取点1、2、3、4、10、11,其中,点4是边缘气隙与主气隙在转子边界线上的分界点,点2及点11分别是两个边缘气隙边界线与转子边界线的交点,点3及点10分别是两个边缘气隙平均磁密线与转子边界线的交点,边缘气隙平均磁密线将两个边缘气隙区域分别划分为f1区域和f2区域;从定子和转子的交叠活动区域的主气隙中心线上选取点5、6,其中,点5是主气隙中心线与转子边界线上的交点,其中,点6是主气隙中心线与定子边界线上的交点;从定子和转子的非交叠活动区域的定子边界线上选取点7、8,其中,点7是边缘气隙与主气隙在定子边界线上的分界点,点8是边缘气隙平均磁密线与定子边界线的交点;边缘气隙边界线近似为主气隙长度的直线加1/4圆弧轨迹线,从边缘气隙平均磁密线上选取两段轨迹线的连接点9;选取路径1→11为麦克斯韦应力法积分路径。
2.如权利要求1所述的计及边缘扩散效应的BSRM建模方法,其特征在于,所述根据所述麦克斯韦应力法积分路径,计算θ=0°和θ≠0°时的气隙磁密,具体包括:
θ=0°时,根据路径1→4,计算得到A相气隙的各气隙磁密;
θ≠0°时,根据路径1→11,计算得到A相气隙的各气隙磁密。
3.如权利要求2所述的计及边缘扩散效应的BSRM建模方法,其特征在于,所述θ=0°时,根据路径1→4,计算得到A相气隙的各气隙磁密,具体包括:
θ=0°时,根据路径1→4,计算得到A相气隙的各气隙磁密,其中,计算得到的A相气隙aj在路径2→3处的主气隙平均磁密为
<mrow>
<msub>
<mover>
<mi>B</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>j</mi>
<mn>0</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>&mu;</mi>
<mn>0</mn>
</msub>
<msub>
<mi>U</mi>
<mrow>
<mi>a</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>/</mo>
<mover>
<msub>
<mi>l</mi>
<mn>0</mn>
</msub>
<mo>&OverBar;</mo>
</mover>
<mo>,</mo>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>3</mn>
<mo>,</mo>
<mn>4</mn>
<mo>;</mo>
</mrow>
式中, 为主气隙磁路平均长度,η为气隙边缘系数,μ0为真空磁导率,l0为定、转子间气隙长度,气隙a1~a4处的磁路磁动势为
Nm为主绕组匝数,Ns为悬浮绕组匝数,ima为A相主绕组电流,isa1为A相α方向悬浮绕组电流,isa2为A相β方向悬浮绕组电流。
4.如权利要求3所述的计及边缘扩散效应的BSRM建模方法,其特征在于,所述θ≠0°时,根据路径1→11,计算得到A相气隙的各气隙磁密,具体包括:
θ≠0°时,根据路径1→11,计算得到A相气隙的各气隙磁密,其中,
计算得到的A相气隙aj在路径5→6处的主气隙磁密为Bmaj=μ0Uaj/l0,j=1,2,3,4;
计算得到的A相气隙aj在路径4→5和路径6→7处的主气隙平均磁密为
计算得到的A相气隙aj在路径8→9和路径9→10处的边缘气隙平均磁密为式中, 为路径8→9和路径9→10处的边缘气隙磁路平均长度,r为转子半径;
计算得到的A相气隙aj在路径2→3和路径10→11处的边缘气隙平均磁密为式中,为路径2→3和路径10→11所在的边缘气隙f1区域的气隙磁路平均长度,
计算得到的A相气隙aj在路径3→4和路径7→8处的边缘气隙平均磁密为式中,为路径3→4和路径7→8所在的边缘气隙f2区域的气隙磁路平均长度,为路径3→4和路径7→8所在的边缘气隙f2区域的极弧角。
5.如权利要求4所述的计及边缘扩散效应的BSRM建模方法,其特征在于,所述根据所述气隙磁密计算得到转子齿极所受α、β方向悬浮力和电磁转矩,具体包括:
根据计算到的θ=0°时的气隙磁密,计算得到A相气隙a1~a4处转子齿极所受α方向悬浮力Fα1~Fα4、所受β方向悬浮力Fβ1~Fβ4、电磁转矩T1~T4分别为:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>&alpha;</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>k</mi>
<mn>0</mn>
</msub>
<msubsup>
<mi>U</mi>
<mrow>
<mi>a</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>&alpha;</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>&alpha;</mi>
<mn>3</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>k</mi>
<mn>0</mn>
</msub>
<msubsup>
<mi>U</mi>
<mrow>
<mi>a</mi>
<mn>3</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>&alpha;</mi>
<mn>4</mn>
</mrow>
</msub>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>&beta;</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>&beta;</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>k</mi>
<mn>0</mn>
</msub>
<msubsup>
<mi>U</mi>
<mrow>
<mi>a</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>&beta;</mi>
<mn>3</mn>
</mrow>
</msub>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>&beta;</mi>
<mn>4</mn>
</mrow>
</msub>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>k</mi>
<mn>0</mn>
</msub>
<msubsup>
<mi>U</mi>
<mrow>
<mi>a</mi>
<mn>4</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<msub>
<mi>T</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<msub>
<mi>T</mi>
<mn>3</mn>
</msub>
<mo>=</mo>
<msub>
<mi>T</mi>
<mn>4</mn>
</msub>
<mo>=</mo>
<mn>0</mn>
</mrow>
式中,k0为θ=0°时悬浮力比例系数,h为转子叠片长度,τr为转子齿极弧度;
根据计算到的θ<0°时的气隙磁密,计算得到A相气隙a1~a4处转子齿极所受α方向悬浮力Fα1~Fα4、所受β方向悬浮力Fβ1~Fβ4、电磁转矩T1~T4分别为:
<mrow>
<msub>
<mi>T</mi>
<mi>j</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>h</mi>
<mi>r</mi>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>&mu;</mi>
<mn>0</mn>
</msub>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msubsup>
<mo>&Integral;</mo>
<mn>5</mn>
<mn>6</mn>
</msubsup>
<msubsup>
<mi>B</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>j</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mi>d</mi>
<mi>l</mi>
<mo>-</mo>
<msubsup>
<mo>&Integral;</mo>
<mn>8</mn>
<mn>10</mn>
</msubsup>
<msubsup>
<mover>
<mi>B</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>f</mi>
<mi>a</mi>
<mi>j</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mi>d</mi>
<mi>l</mi>
<mo>+</mo>
<msubsup>
<mo>&Integral;</mo>
<mn>10</mn>
<mn>11</mn>
</msubsup>
<msubsup>
<mover>
<mi>B</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>f</mi>
<mi>l</mi>
<mi>a</mi>
<mi>j</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mi>d</mi>
<mi>l</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>3</mn>
<mo>,</mo>
<mn>4</mn>
</mrow>
式中,γ为积分路径法向量与转子齿极中心线方向的夹角, 为路径8→9和定子齿极中线的极弧角;
整理后得到
<mrow>
<mo>{</mo>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>&alpha;</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>k</mi>
<mi>&alpha;</mi>
</msub>
<msubsup>
<mi>U</mi>
<mrow>
<mi>a</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>&alpha;</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>k</mi>
<mi>&beta;</mi>
</msub>
<msubsup>
<mi>U</mi>
<mrow>
<mi>a</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>&alpha;</mi>
<mn>3</mn>
</mrow>
</msub>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>k</mi>
<mi>&alpha;</mi>
</msub>
<msubsup>
<mi>U</mi>
<mrow>
<mi>a</mi>
<mn>3</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>&alpha;</mi>
<mn>4</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>k</mi>
<mi>&beta;</mi>
</msub>
<msubsup>
<mi>U</mi>
<mrow>
<mi>a</mi>
<mn>4</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>,</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>&beta;</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>k</mi>
<mi>&beta;</mi>
</msub>
<msubsup>
<mi>U</mi>
<mrow>
<mi>a</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>&beta;</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>k</mi>
<mi>&alpha;</mi>
</msub>
<msubsup>
<mi>U</mi>
<mrow>
<mi>a</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>&beta;</mi>
<mn>3</mn>
</mrow>
</msub>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>k</mi>
<mi>&beta;</mi>
</msub>
<msubsup>
<mi>U</mi>
<mrow>
<mi>a</mi>
<mn>3</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>&beta;</mi>
<mn>4</mn>
</mrow>
</msub>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>k</mi>
<mi>&alpha;</mi>
</msub>
<msubsup>
<mi>U</mi>
<mrow>
<mi>a</mi>
<mn>4</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<msubsup>
<mi>&epsiv;U</mi>
<mrow>
<mi>a</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>T</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<msubsup>
<mi>&epsiv;U</mi>
<mrow>
<mi>a</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>T</mi>
<mn>3</mn>
</msub>
<mo>=</mo>
<msubsup>
<mi>&epsiv;U</mi>
<mrow>
<mi>a</mi>
<mn>3</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>T</mi>
<mn>4</mn>
</msub>
<mo>=</mo>
<msubsup>
<mi>&epsiv;U</mi>
<mrow>
<mi>a</mi>
<mn>4</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mrow>
式中,
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>k</mi>
<mi>&alpha;</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&mu;</mi>
<mn>0</mn>
</msub>
<mi>h</mi>
</mrow>
<mn>2</mn>
</mfrac>
<mo>{</mo>
<mfrac>
<mrow>
<msub>
<mi>l</mi>
<mn>0</mn>
</msub>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<mo>|</mo>
<mi>&theta;</mi>
<mo>|</mo>
<mo>/</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
<msubsup>
<mi>l</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mi>r</mi>
<mo>&lsqb;</mo>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&tau;</mi>
<mi>r</mi>
</msub>
<mo>/</mo>
<mn>2</mn>
<mo>-</mo>
<mo>|</mo>
<mi>&theta;</mi>
<mo>|</mo>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&tau;</mi>
<mi>r</mi>
</msub>
<mo>/</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<msup>
<msub>
<mover>
<mi>l</mi>
<mo>&OverBar;</mo>
</mover>
<mn>0</mn>
</msub>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mi>l</mi>
<mi>sin</mi>
<msub>
<mover>
<mi>&theta;</mi>
<mo>&OverBar;</mo>
</mover>
<mi>f</mi>
</msub>
<mo>+</mo>
<mi>r</mi>
<mover>
<mi>&theta;</mi>
<mo>~</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>cos</mi>
<msub>
<mover>
<mi>&theta;</mi>
<mo>&OverBar;</mo>
</mover>
<mi>f</mi>
</msub>
<mo>-</mo>
<mi>sin</mi>
<mo>|</mo>
<mi>&theta;</mi>
<mo>|</mo>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msup>
<msub>
<mover>
<mi>l</mi>
<mo>&OverBar;</mo>
</mover>
<mi>f</mi>
</msub>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mi>r</mi>
<mo>&lsqb;</mo>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&tau;</mi>
<mi>r</mi>
</msub>
<mo>/</mo>
<mn>2</mn>
<mo>+</mo>
<mo>|</mo>
<mi>&theta;</mi>
<mo>|</mo>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&tau;</mi>
<mi>r</mi>
</msub>
<mo>/</mo>
<mn>2</mn>
<mo>+</mo>
<mover>
<mi>&theta;</mi>
<mo>~</mo>
</mover>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mo>|</mo>
<mi>&theta;</mi>
<mo>|</mo>
<mo>-</mo>
<mover>
<mi>&theta;</mi>
<mo>~</mo>
</mover>
<mo>)</mo>
</mrow>
<mi>sin</mi>
<mo>|</mo>
<mi>&theta;</mi>
<mo>|</mo>
<mo>&rsqb;</mo>
</mrow>
<msubsup>
<mover>
<mi>l</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>f</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mi>r</mi>
<mo>&lsqb;</mo>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&tau;</mi>
<mi>r</mi>
</msub>
<mo>/</mo>
<mn>2</mn>
<mo>+</mo>
<mover>
<mi>&theta;</mi>
<mo>~</mo>
</mover>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&tau;</mi>
<mi>r</mi>
</msub>
<mo>/</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>sin</mi>
<msub>
<mover>
<mi>&theta;</mi>
<mo>&OverBar;</mo>
</mover>
<mi>f</mi>
</msub>
<mo>-</mo>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&tau;</mi>
<mi>r</mi>
</msub>
<mo>/</mo>
<mn>2</mn>
<mo>-</mo>
<mo>|</mo>
<mi>&theta;</mi>
<mo>|</mo>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<msubsup>
<mover>
<mi>l</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>f</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>}</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>k</mi>
<mi>&beta;</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&mu;</mi>
<mn>0</mn>
</msub>
<mi>h</mi>
</mrow>
<mn>2</mn>
</mfrac>
<mo>{</mo>
<mfrac>
<mrow>
<msub>
<mi>l</mi>
<mn>0</mn>
</msub>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<mo>|</mo>
<mi>&theta;</mi>
<mo>|</mo>
<mo>/</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
<msubsup>
<mi>l</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mi>r</mi>
<mo>&lsqb;</mo>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&tau;</mi>
<mi>r</mi>
</msub>
<mo>/</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&tau;</mi>
<mi>r</mi>
</msub>
<mo>/</mo>
<mn>2</mn>
<mo>-</mo>
<mo>|</mo>
<mi>&theta;</mi>
<mo>|</mo>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<msup>
<msub>
<mover>
<mi>l</mi>
<mo>&OverBar;</mo>
</mover>
<mn>0</mn>
</msub>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mi>l</mi>
<mi>cos</mi>
<msub>
<mover>
<mi>&theta;</mi>
<mo>&OverBar;</mo>
</mover>
<mi>f</mi>
</msub>
<mo>+</mo>
<mi>r</mi>
<mover>
<mi>&theta;</mi>
<mo>~</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>cos</mi>
<mo>|</mo>
<mi>&theta;</mi>
<mo>|</mo>
<mo>-</mo>
<mi>sin</mi>
<msub>
<mover>
<mi>&theta;</mi>
<mo>&OverBar;</mo>
</mover>
<mi>f</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msup>
<msub>
<mover>
<mi>l</mi>
<mo>&OverBar;</mo>
</mover>
<mi>f</mi>
</msub>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mi>r</mi>
<mo>&lsqb;</mo>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&tau;</mi>
<mi>r</mi>
</msub>
<mo>/</mo>
<mn>2</mn>
<mo>+</mo>
<mo>|</mo>
<mi>&theta;</mi>
<mo>|</mo>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&tau;</mi>
<mi>r</mi>
</msub>
<mo>/</mo>
<mn>2</mn>
<mo>+</mo>
<mover>
<mi>&theta;</mi>
<mo>~</mo>
</mover>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mo>|</mo>
<mi>&theta;</mi>
<mo>|</mo>
<mo>-</mo>
<mover>
<mi>&theta;</mi>
<mo>~</mo>
</mover>
<mo>)</mo>
</mrow>
<mi>cos</mi>
<mo>|</mo>
<mi>&theta;</mi>
<mo>|</mo>
<mo>&rsqb;</mo>
</mrow>
<msubsup>
<mover>
<mi>l</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>f</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mi>r</mi>
<mo>&lsqb;</mo>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&tau;</mi>
<mi>r</mi>
</msub>
<mo>/</mo>
<mn>2</mn>
<mo>+</mo>
<mover>
<mi>&theta;</mi>
<mo>~</mo>
</mover>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&tau;</mi>
<mi>r</mi>
</msub>
<mo>/</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&tau;</mi>
<mi>r</mi>
</msub>
<mo>/</mo>
<mn>2</mn>
<mo>-</mo>
<mo>|</mo>
<mi>&theta;</mi>
<mo>|</mo>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>cos</mi>
<msub>
<mover>
<mi>&theta;</mi>
<mo>&OverBar;</mo>
</mover>
<mi>f</mi>
</msub>
<mo>&rsqb;</mo>
</mrow>
<msubsup>
<mover>
<mi>l</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>f</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>}</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<mi>&epsiv;</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&mu;</mi>
<mn>0</mn>
</msub>
<mi>h</mi>
</mrow>
<mn>2</mn>
</mfrac>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>l</mi>
<mn>0</mn>
</msub>
</mfrac>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<msub>
<mover>
<mi>l</mi>
<mo>&OverBar;</mo>
</mover>
<mi>f</mi>
</msub>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mo>|</mo>
<mi>&theta;</mi>
<mo>|</mo>
<mo>-</mo>
<mover>
<mi>&theta;</mi>
<mo>~</mo>
</mover>
<mo>)</mo>
</mrow>
</mrow>
<msubsup>
<mover>
<mi>l</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>f</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
根据θ<0°时的A相气隙a1~a4处转子齿极所受α方向悬浮力、所受β方向悬浮力和电磁转矩,同理计算得到θ>0°时的A相气隙a1~a4处转子齿极所受α方向悬浮力Fα1~Fα4、所受β方向悬浮力Fβ1~Fβ4、电磁转矩T1~T4分别为
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>&alpha;</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>k</mi>
<mi>&alpha;</mi>
</msub>
<msubsup>
<mi>U</mi>
<mrow>
<mi>a</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>&alpha;</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>k</mi>
<mi>&beta;</mi>
</msub>
<msubsup>
<mi>U</mi>
<mrow>
<mi>a</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>&alpha;</mi>
<mn>3</mn>
</mrow>
</msub>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>k</mi>
<mi>&alpha;</mi>
</msub>
<msubsup>
<mi>U</mi>
<mrow>
<mi>a</mi>
<mn>3</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>&alpha;</mi>
<mn>4</mn>
</mrow>
</msub>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>k</mi>
<mi>&beta;</mi>
</msub>
<msubsup>
<mi>U</mi>
<mrow>
<mi>a</mi>
<mn>4</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>&beta;</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>k</mi>
<mi>&beta;</mi>
</msub>
<msubsup>
<mi>U</mi>
<mrow>
<mi>a</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>&beta;</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>k</mi>
<mi>&alpha;</mi>
</msub>
<msubsup>
<mi>U</mi>
<mrow>
<mi>a</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>&beta;</mi>
<mn>3</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>k</mi>
<mi>&beta;</mi>
</msub>
<msubsup>
<mi>U</mi>
<mrow>
<mi>a</mi>
<mn>3</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>&beta;</mi>
<mn>4</mn>
</mrow>
</msub>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>k</mi>
<mi>&alpha;</mi>
</msub>
<msubsup>
<mi>U</mi>
<mrow>
<mi>a</mi>
<mn>4</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<msubsup>
<mi>&epsiv;U</mi>
<mrow>
<mi>a</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>T</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<msubsup>
<mi>&epsiv;U</mi>
<mrow>
<mi>a</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>T</mi>
<mn>3</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<msubsup>
<mi>&epsiv;U</mi>
<mrow>
<mi>a</mi>
<mn>3</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>T</mi>
<mn>4</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<msubsup>
<mi>&epsiv;U</mi>
<mrow>
<mi>a</mi>
<mn>4</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
kα和kβ分别为θ≠0°时相互垂直方向上悬浮力比例系数,ε为θ≠0°时电磁转矩系数。
6.如权利要求5所述的计及边缘扩散效应的BSRM建模方法,其特征在于,所述建立BSRM的α、β方向悬浮力和电磁转矩解析模型,具体包括:
当BSRM在线性工作状态下且仅A相通电时,α方向悬浮力Fα、β方向悬浮力Fβ、电磁转矩T分别为:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>F</mi>
<mi>&alpha;</mi>
</msub>
<mo>=</mo>
<msub>
<mi>F</mi>
<mrow>
<mi>&alpha;</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>F</mi>
<mrow>
<mi>&alpha;</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>F</mi>
<mrow>
<mi>&alpha;</mi>
<mn>3</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>F</mi>
<mrow>
<mi>&alpha;</mi>
<mn>4</mn>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>F</mi>
<mi>&beta;</mi>
</msub>
<mo>=</mo>
<msub>
<mi>F</mi>
<mrow>
<mi>&beta;</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>F</mi>
<mrow>
<mi>&beta;</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>F</mi>
<mrow>
<mi>&beta;</mi>
<mn>3</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>F</mi>
<mrow>
<mi>&beta;</mi>
<mn>4</mn>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>T</mi>
<mo>=</mo>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>T</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>T</mi>
<mn>3</mn>
</msub>
<mo>+</mo>
<msub>
<mi>T</mi>
<mn>4</mn>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
根据θ=0°、θ<0°、θ>0°时的A相气隙a1~a4处转子齿极所受α方向悬浮力、所受β方向悬浮力和电磁转矩,对上式转换整理得到BSRM的α、β方向悬浮力和电磁转矩解析模型为
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>F</mi>
<mi>&alpha;</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>F</mi>
<mi>&beta;</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<msub>
<mi>i</mi>
<mi>m</mi>
</msub>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>K</mi>
<mi>&alpha;</mi>
</msub>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>K</mi>
<mi>&beta;</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>K</mi>
<mi>&beta;</mi>
</msub>
</mtd>
<mtd>
<msub>
<mi>K</mi>
<mi>&alpha;</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>i</mi>
<mrow>
<mi>s</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>i</mi>
<mrow>
<mi>s</mi>
<mn>2</mn>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>T</mi>
<mo>=</mo>
<msub>
<mi>J</mi>
<mi>t</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>2</mn>
<msubsup>
<mi>N</mi>
<mi>m</mi>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>i</mi>
<mi>m</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>N</mi>
<mi>s</mi>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>i</mi>
<mrow>
<mi>s</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>N</mi>
<mi>s</mi>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>i</mi>
<mrow>
<mi>s</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
式中,im为通电相主绕组电流;is1为通电相α方向悬浮绕组电流;is2为通电相β方向悬浮绕组电流;悬浮力比例系数Kα、Kβ和电磁转矩系数Jt分别为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510916257.6A CN105373014B (zh) | 2015-12-11 | 2015-12-11 | 计及边缘扩散效应的bsrm建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510916257.6A CN105373014B (zh) | 2015-12-11 | 2015-12-11 | 计及边缘扩散效应的bsrm建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105373014A CN105373014A (zh) | 2016-03-02 |
CN105373014B true CN105373014B (zh) | 2018-01-16 |
Family
ID=55375308
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510916257.6A Active CN105373014B (zh) | 2015-12-11 | 2015-12-11 | 计及边缘扩散效应的bsrm建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105373014B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107769623B (zh) * | 2017-12-04 | 2019-04-26 | 河海大学 | 一种绕组关断角同时调节的bsrm期望电流计算方法 |
CN108955984B (zh) * | 2018-05-31 | 2020-07-14 | 沈阳工业大学 | 一种可直接指示方向的主应力方向判定方法 |
CN108875168B (zh) * | 2018-06-04 | 2021-03-26 | 同济大学 | 一种考虑饱和的开关磁阻电机磁场解析计算方法 |
CN110298098B (zh) * | 2019-06-21 | 2023-03-31 | 河海大学 | 定转子非重叠区域的bsrm麦克斯韦应力解析建模方法 |
CN110287586B (zh) * | 2019-06-21 | 2022-11-08 | 河海大学 | 适用于麦克斯韦应力法bsrm精确建模的气隙细分方法 |
CN110661390B (zh) * | 2019-09-24 | 2021-05-25 | 江苏大学 | 一种12/14极磁悬浮开关磁阻电机悬浮力精确建模方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6504337B1 (en) * | 2001-06-08 | 2003-01-07 | General Electric Company | Motor modeling using harmonic ampere-turn saturation method |
CN101355281B (zh) * | 2008-03-10 | 2010-12-08 | 江苏大学 | 无轴承电机中磁场等效气隙虚拟绕组电流的分析方法 |
-
2015
- 2015-12-11 CN CN201510916257.6A patent/CN105373014B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN105373014A (zh) | 2016-03-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105373014B (zh) | 计及边缘扩散效应的bsrm建模方法 | |
Rahideh et al. | Optimal brushless DC motor design using genetic algorithms | |
Lubin et al. | Analytic calculation of eddy currents in the slots of electrical machines: Application to cage rotor induction motors | |
CN105245041B (zh) | 一种求解大型电机定转子通风沟表面散热系数的方法 | |
Seo et al. | Analysis of overhang effect for a surface-mounted permanent magnet machine using a lumped magnetic circuit model | |
Farooq et al. | Analytical modeling approach to detect magnet defects in permanent-magnet brushless motors | |
Putek et al. | Topology optimization of rotor poles in a permanent-magnet machine using level set method and continuum design sensitivity analysis | |
CN109768685B (zh) | 基于单槽模型的表贴式永磁电机电枢磁场重构方法 | |
Kappatou et al. | Design optimization of axial flux permanent magnet (AFPM) synchronous machine using 3D FEM analysis | |
Ghasemi | Cogging torque reduction and optimization in surface-mounted permanent magnet motor using magnet segmentation method | |
Pina et al. | Modeling of synchronous reluctance motors aided by permanent magnets with asymmetric rotor poles | |
CN104361227B (zh) | 一种定子铁心背部支持筋感应电压计算方法 | |
Martin et al. | Homogenization technique for axially laminated rotors of synchronous reluctance machines | |
Barré et al. | Concentrated windings in compact permanent magnet synchronous generators: managing efficiency | |
Petrichenko et al. | Design and simulation of turbo-alternators using a coupled permeance network model | |
Bischof et al. | Synthesized magnetic flux density in three-phase cage induction machines with MATLAB & Simulink | |
Varaticeanu et al. | Performance evaluation of permanent magnet assisted synchronous reluctance motor for micro electric vehicle | |
CN109033613A (zh) | 一种无轴承磁通切换电机转子悬浮力快速建模方法 | |
Zhilichev | Analysis of permanent magnet machines using crossing macro-elements | |
Yang et al. | Design and characteristic analysis of a novel bearingless SRM considering decoupling between torque and suspension force | |
Zhu et al. | Decoupling Control for Dual‐Winding Bearingless Switched Reluctance Motor Based on Improved Inverse System Method | |
Imamori et al. | Relation between magnetic properties of stator core and cogging torque in 8-pole 12-slot SPM synchronous motors | |
Li et al. | Analysis of magnetic field and levitation force characteristics for 3-DOF deflection type PM motors | |
Nakatsugawa et al. | Proposal of two‐dimensional magnetic field analysis model for transverse magnetic flux type linear synchronous motor | |
JP5905206B2 (ja) | コギングトルクを低減するモータおよびリニアモータ |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |