低转矩波动连续极永磁同步电机磁极结构稳健性设计方法与流程

1.本发明属于电机优化设计领域的一种永磁同步电机结构设计和布置方法，具体涉及一种能够减小转矩波动的连续极永磁同步电机(cppmsm)磁极结构稳健性设计方法。


背景技术：

2.稀土永磁同步电机(pmsm)因具有结构简单、运行可靠、体积小、重量轻、损耗少、效率高、结构灵活多样等显著特点，被广泛用于航空航天、国防及工农业生产和日常生活的各个领域。但是稀土材料价格昂贵，导致电机成本过高。为了降低电机成本，cppmsm被提出。cppmsm与表贴式永磁同步电机(spmsm)相比，用铁极代替了极性相同的一半永磁体，因此节约了大量永磁体，降低了电机成本。但是cppmsm的特殊结构导致电机磁场分布不对称，偶次谐波较高，使得电机转矩波动较大，不利于电机的平稳运行。
3.为了使电机各项性能达到综合最优，现有技术需要对电机进行优化设计。


技术实现要素：

4.为了解决背景技术中存在的问题，本发明提出了一种用于cppmsm磁极结构的稳健性设计方法，以在保证较高输出转矩的同时，减小电机的转矩波动，并保证最终结构方案有较好的稳健性。
5.本发明优化了各优化变量变化对转矩平均值和转矩波动的影响，对改进的电磁结构进行优化，使用优化后的电磁结构能够在保持较高转矩平均值的同时降低电机转矩波动，有效提升了电机运行的平稳性。
6.本发明所采用的技术方案是：
7.1)在连续极永磁同步电机(cppmsm)的磁极结构基础上进行改进，提出电机的磁极改进方案，采用在铁极处添加一块永磁体且增加极间气隙而形成混合连续极永磁同步电机(hcppmsm)磁极结构；
8.2)选择优化变量，建立优化目标，对改进后的混合连续极永磁同步电机磁极结构进行优化设置，实现磁极结构稳健性设计。
9.所述的连续极永磁同步电机的磁极结构包括环形的铁芯和布置在铁芯外周围的铁极和磁极，数量相同的多个铁极和多个磁极沿铁芯圆周交替紧密布置，铁极和磁极之间不具有气隙；
10.在此基础上，磁极保持不变，每个铁极和磁极之间设置气隙形成极间气隙，每个铁极在沿同一周向的一半侧改设为辅助永磁体，且辅助永磁体和铁极的另一半侧固定连接形成新的铁极，新的铁极中的辅助永磁体和铁极部分的两侧均和磁极之间设置间隙。
11.所述磁极采用永磁体。铁极和磁极沿径向的厚度相同。
12.所述的铁极的磁极方向和磁极的充磁方向均沿径向方向，铁极中新增设的辅助永磁体的充磁方向也沿径向方向。
13.所述的连续极永磁同步电机的磁极结构中，铁极和磁极的磁方向交替相反布置，
磁方向包括n极和s极，以n极弧长圆心角θ1和s极弧长圆心角θ2的比值为优化变量a，以n极和s极之间的气隙长度圆心角θ3为优化变量b，以辅助永磁体弧长圆心角θ4与s极弧长圆心角θ2的比值为优化变量c，以磁极沿径向的厚度为优化变量d；以减小转矩波动和提高转矩平均值作为优化目标。
14.所述步骤(2)具体如下：
15.2.1)结合电机的结构参数确定各个优化变量的取值，设置水平建立控制因素水平表，控制因素水平表是由每个优化变量在每个水平下的取值构成；
16.同时根据选取的优化变量的个数以及各优化变量对应的取值个数获得试验的不同排列组合情况，进而建立控制因素正交表，控制因素正交表是由每个优化变量在每种控制因素试验下的不同水平组合搭配构成，所有优化变量的一种水平组合作为一种控制因素试验，不同控制因素试验下各个优化变量所对应的水平组合不同；
17.所述的水平是指优化变量的不同取值。
18.2.2)以可能存在加工误差的变量作为噪声因素确定误差范围，根据控制因素水平表再建立和控制因素水平表相吻合的噪声因素水平表，噪声因素水平表是由每个优化变量在每个带有噪声因素的水平下的取值构成；
19.并根据噪声因素的个数及噪声因素水平表建立噪声因素正交表，噪声因素正交表是由每个优化变量在每种噪声因素试验下的不同带有噪声因素的水平组合搭配构成，所有优化变量的一种带有噪声因素的水平组合作为一种噪声因素试验，不同噪声因素试验下各个优化变量所对应的带有噪声因素的水平组合不同；
20.具体实施中，将不同水平设置不同的误差范围，并施加不同的系数到不同水平上。
21.2.3)以控制因素正交表为外表，噪声因素正交表为内表，将内表和外表直积进行多次正交的总试验，即在控制因素正交表的每次控制因素试验下进行噪声因素正交表的各次噪声因素试验，，一种控制因素试验下进行一种噪声因素试验即作为一种总试验，每一次总试验对电机运行在额定点时进行有限元分析，计算获得每次总试验下电机的转矩波动tr和转矩平均值t；
22.2.4)根据2.3)获得所有总试验的转矩波动tr，依据信噪比的望小特性计算每种控制因素试验下转矩波动tr的信噪比s/n，同时获得每种控制因素试验下转矩波动tr的平均值；
23.根据2.3)获得所有总试验的转矩平均值t，依据信噪比的望大特性计算每种控制因素试验下转矩平均值t的信噪比s/n，同时获得每种控制因素试验下转矩平均值t的平均值；
24.2.5)对各种控制因素试验的转矩波动tr和转矩平均值t的s/n值和平均值进行平均值分析，得到转矩波动、转矩平均值随各优化变量取值的变化情况，得到每个控制因素在每个水平下的转矩平均值及其信噪比的平均值以及每个控制因素在每个水平下的转矩波动及其信噪比的平均值，得到转矩波动的平均值最小时各优化变量的水平取值和转矩平均值的平均值最大时各优化变量的水平取值；
25.2.6)在2.5)平均值分析的基础上对信噪比s/n进行贡献率分析：计算每个优化变量的波动平方和，根据每个优化变量的波动平方和得到各个优化变量分别对转矩波动、转矩平均值影响的贡献率；
26.2.7)根据步骤2.6)得到的各优化变量对转矩波动、转矩平均值影响的贡献率，综合考虑后得到各优化变量的水平取值，对步骤2.5)得到的转矩波动最小、转矩平均值最大时各优化变量的水平取值进行综合处理，确定磁极结构中优化变量的最终水平取值。
27.所述2.7)中，对于每个优化变量，在转矩波动最小时的水平取值和在转矩平均值最大的的水平取值的两种不同优化目标下获得的结果中，选取贡献率较大的所对应的水平取值作为该优化变量的最终水平取值；
28.将所有优化变量的最终水平取值施加到混合连续极永磁同步电机(hcppmsm)磁极结构上进行工作，从而实现混合连续极永磁同步电机(hcppmsm)磁极结构的稳健工作。
29.所述的步骤2.6)具体为：
30.计算转矩平均值的信噪比s/n、转矩波动的信噪比s/n分别在各优化变量下的波动平方和，以转矩平均值的信噪比s/n在优化变量a下波动平方和计算为例说明各个优化变量的波动平方和计算，计算方法如下公式：
[0031][0032]
其中，s
ηta
表示优化变量a的波动平方和，j表示水平的序数，l
ηta(j)
为优化变量a在水平j下的转矩平均值的信噪比的平均值，l
ηt
为优化变量a在所有水平下的转矩平均值的信噪比的平均值，q为优化变量a涉及的水平总数。
[0033]
求取转矩平均值的信噪比s/n在所有优化变量下的波动平方和之和作为转矩平均值的信噪比s/n的总波动平方和，求取转矩波动的信噪比s/n在所有优化变量下的波动平方和之和作为转矩波动的信噪比s/n的总波动平方和，计算转矩平均值的信噪比s/n在某优化变量下的波动平方和除以转矩平均值的信噪比s/n的总波动平方和作为该优化变量对转矩平均值影响的贡献率，计算转矩波动的信噪比s/n在某优化变量下的波动平方和除以转矩波动的信噪比s/n的总波动平方和作为该优化变量对转矩波动影响的贡献率，从而获得了各个优化变量分别对转矩平均值、转矩波动影响的贡献率。
[0034]
本发明的核心是提出一种能有效降低cppmsm转矩波动的磁极结构，并选择控制因素正交表及噪声因素正交表方法处理对提出的磁极结构进行优化。
[0035]
电机为多变量耦合系统，对其进行优化设计时需要同时考虑多个变量，随着变量个数的增加，试验次数将呈指数增加，耗费较多的计算资源与时间。在进行多目标优化设计时，采用正交试验表，可极大的减少试验次数，提升电机优化设计效率，且能达到较优的优化设计效果。
[0036]
本发明改进了cppmsm磁极结构，同时对改进后的hcppmsm磁极结构进行优化，通过对hcppmsm进行稳健性设计使电机在保持较高输出转矩的同时，减小了电机的转矩波动。具有以下有益效果：
[0037]
1、本发明改进了cppmsm磁极结构提出了hcppmsm磁极结构，改进后的结构有效降低了气隙磁场中的谐波分量，使电机转矩波动得到了显著降低，提高了电机运行平稳性。
[0038]
2、对hcppmsm磁极结构进行了优化，分析了转矩平均值和转矩波动随各优化变量
变化的规律，和各优化变量对转矩平均值和转矩波动影响的贡献率，进而得到了hcppmsm磁极结构的最终优化方案，使得电机在保持较高输出转矩的同时，减小了电机的转矩波动。
附图说明
[0039]
图1为使用磁极改进结构前的cppmsm磁极结构图；
[0040]
图2为使用磁极改进结构后的hcppmsm磁极结构图；
[0041]
图3为磁极改进结构的优化变量示意图。
具体实施方式
[0042]
下面结合附图和具体实施对本发明作进一步说明。
[0043]
本发明的实施例及其实施情况如下：
[0044]
下面以一台10极12槽电机为例，对本发明的实施方式进行详细说明，电机的参数如表1所示。
[0045]
表1分数槽spmsm参数
[0046]
参数符号数值单位额定转速n
n
3000r/min额定转矩t
n
2.4068nm极对数p5
‑‑
槽数q12
‑‑
转子气隙处半径r
ra
15.75mm气隙长度δ0.5mm铁心长度l82mm额定电流i4.7a
[0047]
cppmsm的初始磁极结构由磁极和铁极组成，如图1所示；
[0048]
提出如图2所示的hcppmsm电机改进磁极结构，即在cppmsm的铁极一侧添加一块永磁体，且增加了极间气隙；以减小电机的转矩波动和增大电机转矩平均值为优化目标，其中优化变量的选取如图3所示：
[0049]
1)结合电机的结构参数，确定优化变量a、b、c和d的取值范围，在每个优化变量取值范围内均匀选取3个水平的值，建立控制因素水平表，如表2所示。根据优化变量个数及各优化变量水平数建立控制因素正交表l9(34)，如表3所示；
[0050]
表2控制因素水平表
[0051][0052]
其中，1、2、3分别表示水平的序数。
[0053]
表3控制因素l9(34)正交表
[0054][0055][0056]
2)受电机加工工艺的影响，上述四个优化变量均可能存在误差，以上四个优化变量在作为控制因素的同时作为噪声因素考虑，选择
±
1％作为加工误差，确定误差范围，建立噪声因素水平表和噪声因素正交表，如表4和表5所示；
[0057]
表4噪声因素水平表
[0058][0059]
表5噪声因素正交表
[0060]
试验次数abcd111112122231333421235223162312731328321393321
[0061]
3)以控制因素正交表为外表，噪声因素正交表为内表，将内外表直积进行试验，对电机运行在额定点时的每一次试验进行有限元分析，计算电机的转矩平均值t如表6所示和转矩波动t
r
如表7所示；
[0062]
表6各组试验方案转矩平均值
[0063] t1t2t3t4t5t6t7t8t912.34332.35332.35632.33272.32482.33242.36832.38122.366322.36942.37212.37762.39072.38022.38422.39892.40272.389932.32592.32692.35922.34002.32882.35832.32682.32972.320642.37362.38852.38882.38872.38552.40062.36442.39382.393352.37292.37572.37852.37022.36012.36390.36742.37802.369262.32532.32682.36012.34602.33042.36182.33432.33612.3214
72.33722.34932.37412.36192.34962.35102.31862.34812.327682.39352.39852.40022.40742.39652.39922.41272.41752.406292.30142.30942.31772.33462.35052.32072.34842.35702.3484
[0064]
表7各组试验方案转矩波动
[0065][0066][0067]
4)以依据信噪比的望大特性和望小特性，分别计算各组试验对应的转矩平均值和转矩波动结果的信噪比s/n值和平均值，如表8所示；
[0068]
表8各组试验方案结果平均值及sn值
[0069][0070]
其中，η
t
表示转矩平均值的信噪比s/n，η
tr
表示转矩波动的信噪比s/n。
[0071]
针对转矩平均值，根据望大特性求取信噪比；针对转矩波动，根据望小特性求取信噪比。
[0072]
5)对各组试验得到的转矩平均值的信噪比s/n和平均值及转矩波动的信噪比s/n和平均值进行平均值分析处理，结果如表9所示。
[0073]
表9转矩平均值、转矩波动及相应sn值在各优化变量取值下的平均值
[0074][0075]
其中，l
t
表示转矩平均值的平均值，l
ηt
表示转矩平均值的信噪比的平均值，l
tr
表示转矩波动的平均值，l
ηtr
表示转矩波动的信噪比的平均值。例如2.3571表示优化变量a在水平1下的转矩平均值的平均值，0.0243表示优化变量b在水平3下的转矩波动的平均值。
[0076]
表9中例如，2.3571是由表8中的2.3512、2.3851、2.3351求取平均获得，2.3864是由表8中的2.3851、2.3707、2.4035求取平均获得。
[0077]
由表9可以看出各优化变量变化对转矩平均值和转矩波动的影响规律，随着优化变量a的增大，转矩平均值先增大后降低，转矩波动逐渐增大；随优化变量b的增大，转矩平均值先增大后降低，转矩波动先降低后增大；随优化变量c的增大，转矩平均值先增大后降低，转矩波动逐渐降低；随优化变量d的增大，转矩平均值逐渐增大，转矩波动先增大后降低。同时可以看出，转矩平均值增大时，其对应的s/n值也相应增大，转矩波动降低时，其对应的信噪比s/n值也相应增大，因此可通过选择s/n值大的组合来确定最优方案。
[0078]
由表9可以看出，使转矩平均值的信噪比s/n值最大的各优化变量的水平组合为a(2)b(2)c(2)d(3)；使转矩波动的信噪比sn值最大的各优化变量的水平组合为a(1)b(2)c(3)d(3)。
[0079]
4)在平均值分析的基础上对转矩平均值和转矩波动的信噪比s/n值进行贡献率分析，得到各个优化变量对转矩平均值和转矩波动结果的贡献率，贡献率越高，说明结果受此优化变量变化的影响越大。
[0080]
5)计算各优化变量的波动平方和，然后根据各个优化变量的波动平方和和与总波动平方和的比值计算出各优化变量对转矩平均值和转矩波动影响的贡献率。
[0081]
优化变量的波动平方和是分为转矩平均值的信噪比sn、转矩波动的信噪比sn在优化变量下波动平方和。
[0082]
以转矩平均值的sn值在因素a下波动平方和计算为例来说明各个优化变量的波动平方和计算，计算方法如公式(1)所示。同理可计算出转矩平均值和转矩波动的sn值在所有因素下的波动平方和。总波动平方和为各因素波动平方和之和。计算结果如表10所示。
[0083][0084]
其中，l
ηta(j)
为表9中l
ηt
列因素a水平j的值，l
ηt
为表9中l
ηt
列因素的平均值，q为因素a的水平数，此处为3。
[0085]
表10转矩平均值和转矩波动sn值贡献率分析
[0086][0087][0088]
其中，s
t
表示优化变量对转矩平均值的贡献率，s
tr
表示优化变量对转矩波动的贡献率，p表示贡献率的百分比。
[0089]
由表10可得，4个优化变量对转矩平均值影响的贡献率从大到小依为bdca而对转矩波动影响的贡献率从大到小依次为abcd。
[0090]
6)根据步骤5)得到的转矩平均值最大、转矩波动最小时各优化变量的取值，以及步骤6)得到的各优化变量对转矩平均值和转矩波动结果的贡献率，综合考虑确定各优化变量的取值：
[0091]
优化变量b和优化变量d选取的水平值相同，分别为b(2)和d(3)，优化变量a、c对应的水平值不同，优化变量a影响最大的是转矩波动，因此优化变量a选取使转矩波动最小的水平值，即a(1)，优化变量c对两个优化目标的影响均不显著，本文选取使转矩平均值最大的水平，即c(2)。
[0092]
最终结构选取的四个水平值分别为a(1)b(2)c(2)d(3)。
[0093]
对使用优化后磁极结构的hcppm进行有限元分析，计算其转矩平均值及转矩波动，结果如表11所示。
[0094]
表11磁极结构优化前后转矩平均值、转矩波动对比
[0095] t(nm)t
r
(％)优化前2.40680.0427优化后2.39680.0198变化率0.42％53.63％
[0096]
由表中可以看出，优化后的电机在保持较高转矩平均值的基础上，转矩波动得到了显著降低，提高了电机运行的平稳性。