基于局部代理模型的电机多目标鲁棒性优化方法

1.本发明公开基于局部代理模型的电机多目标鲁棒性优化方法，涉及电机本体优化设计技术，属于计算、推算或计数的技术领域。


背景技术：

2.在大批量电机生产中，产品的合格率至关重要。然而，在电机实际加工生产中，由于材料的分散性以及加工公差的存在，产品和设计之间的出入无可避免。当某些电机性能对于不确定性敏感时，就有可能造成设计违背约束条件，进而导致次品的发生。
3.电机鲁棒性优化的初衷，就是在设计阶段考虑不确定性对电机性能的影响，从而使得参数变化时，电机性能仍然能够满足要求，进而保证批量产品的合格率。gang lei在《robust design optimization of electrical machines:a comparative study and space reduction strategy》一文中提出了一种基于性能函数均值和方差的电机鲁棒性优化方法，即六西格玛设计(six sigma design，ssd)。ssd方法以性能函数的均值和方差关系式代替原始的优化目标和约束条件，但是这种方法默认将约束条件的可靠度设为100％，设计者无法调整，同时该方法过分看重了鲁棒性而非电机性能。
4.bo ma在《robust design optimization of electrical machines considering hybrid random and interval uncertainties》一文中提出了一种基于基函数的混合不确定性分析方法。通过将电机性能函数在不确定性处分解成多项式级数形式，以均值和方差的上界代替原始优化模型。但是，这种方法天然存在着区间缩放导致的固有偏差，并且实施方式较为复杂。
5.gang lei在《robust multiobjective and multidisciplinary design optimization of electrical drive systems》一文中提出了一种基于kriging模型的鲁棒性优化设计方法。通过代理模型计算电机性能，以减小整体优化的计算负担。但是，这种通过空间采样方法建立的代理模型，在不确定性所在的局部小区间内的精度难以达到要求，将会影响最终的优化结果。
6.本发明旨在提出基于局部代理模型的电机多目标鲁棒性优化方法，在不需要增加额外采样数据点的情况下逐步提高代理模型精度。


技术实现要素：

7.本发明的发明目的是针对上述背景技术的不足，提供基于局部代理模型的电机多目标鲁棒性优化方法，通过确定性设计的数据集训练全局代理模型后进行初次鲁棒性优化，在初次鲁棒性优化最优解集临近域内构建局部数据集，根据局部数据集训练局部代理模型后再次进行鲁棒性优化，通过选取局部数据集、训练局部代理模型的迭代鲁棒性优化过程，获取鲁棒性优化的最终解集，解决了全局代理模型在局部小区间内精度不足的技术问题。
8.本发明为实现上述发明目的采用如下技术方案：一种提高电机多目标鲁棒性优化
精度的局部代理模型策略。首先，暂时不考虑电机不确定性，由确定性优化数据构建训练全局代理模型k0的数据库，同时确定电机参数多目标优化模型的优化目标及约束条件，再基于全局代理模型预测电机性能，完成初次多目标鲁棒性优化，进而获取较为准确的初次鲁棒性优化最优解集p0；其次，选择初次鲁棒性优化最优解集p0临近域内的样本构建训练局部代理模型k1，获取更为准确的最优解集p1；然后，比较p0和p1的超体积参数h0和h1，在h0和h1的相对差别满足阈值ε时，终止迭代过程，在h0和h1的相对差别不满足阈值时，则选择最优解集p1临近域内的样本训练局部代理模型k2，获取更为准确的最优解集p2，重复上述过程，直至相邻两次最优解集的超体积参数的差别满足阈值ε；最后，相邻两次最优解集的超体积参数的差别满足阈值ε时，即可得到高精度的鲁棒性最优结果。
9.进一步地，在局部代理模型策略中，确定性优化数据通过求解确定性电机参数多目标优化模型获得。
10.进一步地，在局部代理模型策略中，训练全局代理模型k0的方法为：采用遗传算法对全局代理模型k0进行超参数调优得到具备良好预测能力的全局代理模型k0，再应用nsga
ꢀⅱ
算法进行电机鲁棒性优化产生候选设计，进而结合mca法对候选设计进行鲁棒性评估。
11.进一步地，在局部代理模型策略中，结合mca法对候选设计进行鲁棒性评估的具体方法为：对于nsga
ꢀⅱ
算法得到的每一代候选设计，计算各候选设计满足电机多目标优化模型约束条件的概率，将满足各个约束条件的概率乘积作为评价候选设计鲁棒性的总体数值指标，对各候选设计满足电机多目标优化模型约束条件的概率与满足鲁棒性设计的概率设定值进行比较，根据比较结果给nsgaⅱ算法中的约束函数矩阵赋值，根据赋值后的约束函数矩阵进行下一代候选设计的生成，nsga
ꢀⅱ
算法迭代终止时获取的候选设计构成初次鲁棒性优化最优解集，给出了结合mca进行电机鲁棒性优化的实现方式。
12.进一步地，在局部代理模型策略中，最优解集的临近域是通过对解集中数据进行线性回归以及计算回归函数与解集中数据的最远距离确定的。
13.本发明采用上述技术方案，具有以下有益效果：
14.(1)本发明提供一种提高电机多目标鲁棒性优化精度的局部代理模型策略，首先依据确定性优化数据训练表征电机参数与电机性能之间近似关系的全局代理模型，由训练后的全局代理模型预测电机性能进而实现电机参数的鲁棒性优化，在已确定最优解集的临近域内选取训练局部代理模型的数据，采用迭代训练局部代理模型后进行鲁棒性优化的过程替代有限元计算，极大地削减了优化耗时，且最终获得的鲁棒性优化解集仅位于设计空间的一个相对较小的区域内，能够有效提高代理模型对不确定性影响下电机性能在小区间内的预测精度，不增加额外计算负担，过程简单，易于实现，适用于任意种类任意工作情况下的电机多目标鲁棒性优化设计问题。
15.(2)本发明以确定性设计的部分样本为训练代理模型的初始数据集，在后续迭代鲁棒性优化过程中，在前一次鲁棒性优化最优解集临近域内选取下一次鲁棒优化的数据集，使得下一次鲁棒性优化的数据尽可能地聚集在电机参数不确定性设计空间范围内，局部代理模型在迭代优化的过程中逐渐逼近最优解集附近的电机性能，相较于通过拉丁超立方采样方式获取数据集的鲁棒性优化方法，本发明在不需要增加额外采样数据点的情况下逐步提高代理模型精度，并且可以通过mca法准确求解设计的可靠度数值，最终显著提高电机多目标鲁棒性优化的精度。
附图说明
16.图1为本发明局部代理模型策略的实施流程图。
17.图2为基于确定性设计前25代数据构成的数据库在设计空间的分布图。
18.图3为确定最优解集相邻域的示意图。
19.图4为计算超体积参数的示意图。
20.图5为超体积相对差别与迭代次数相关的折线图。
21.图6为基于全局/局部代理模型的鲁棒性优化结果的示意图。
22.图7为基于全局代理模型的最优解集误差的分布图。
23.图8为基于局部代理模型的最优解集误差的分布图。
具体实施方式
24.下面结合附图对发明的技术方案进行详细说明。
25.电机参数不确定性在实际加工生产过程中不可避免，且会导致部分电机性能偏离要求区间，造成大批量生产中产生大量次品。电机鲁棒性优化受限于庞大的计算负担，通常基于代理模型进行，而采用全局代理模型预测不确定性小区间内电机性能存在着精度不足的缺陷。由于最优解集附近的样本点比分布在整个设计空间内的样本点更能表征最优设计的特点，因此可以在每次优化结束后，用相邻域内的数据点训练局部代理模型，逐次更新最优解集，即可得到高精度的鲁棒性优化结果。
26.现以一个永磁同步电机多目标鲁棒性优化设计为例，详细介绍本发明实施过程。
27.本发明整体实施过程如图1所示，首先，通过确定性设计的数据集建立全局代理模型k0的数据库，应用遗传算法对全局代理模型k0进行超参数调优，得到具备较好预测能力的全局代理模型k0，全局代理模型k0用于预测电机性能，结合nsga
ꢀⅱ
(non-dominated sorting genetic algorithmⅱ，非支配排序遗传算法ⅱ)和mca(monte carlo analysis，蒙特卡洛法分析)指导电机多目标优化并计算电机鲁棒性指标，获取初次鲁棒性优化最优解集p0；其次，对p0进行线性回归，并且筛选出位于最优解集临近域内的局部数据集；然后，利用局部数据集，训练新的局部代理模型k1，并且完成新一次多目标鲁棒性优化，更新最优解集p1；最后，比较相邻两次最优解集的超参数指标，若相对差别大于预先设置的值ε，则继续更新局部代理模型，直到相对差别满足要求，否则优化中止，即可得到鲁棒性优化的最终解集。
28.在发明所提局部代理模型策略中，首先，需要以部分确定性优化的数据构建数据库，并且通过部分确定性优化数据组成的数据库训练全局代理模型。由确定性设计前25代数据构成的数据库在设计空间的分布如图2所示。确定性优化设计的多目标优化模型的表达式为：
[0029][0030]
其中，torque表示电机均值转矩，cost表示电机材料成本，ripple表示转矩波动，
efficiency表示电机效率，t
winding
表示绕组平均温度，t
pm
表示永磁体平均温度。
[0031]
以部分确定性优化数据组成的数据库分别训练转矩、成本、转矩波动、效率、绕组平均温度和永磁体平均温度的全局kriging代理模型k0，全局kriging代理模型k0用于多目标鲁棒性优化，此时多目标优化模型的表达式为：
[0032][0033]
其中，p(gi)为第i个约束条件成立的概率，c是设计者设置的满足鲁棒性设计的概率，在本次示例中c设置为100％。
[0034]
通过全局代理模型k0得到的最优解集为p0，通过线性回归算法拟合最优解集p0得到拟合函数l0。接着计算p0内所有数据点到l0的距离，以直线l0上下距离最远的数据点分别确定边界直线l1和l2，从而得出相邻域范围(l1与l2之间)，如图3所示。
[0035]
数据点(xi，yi)到回归直线ax by c＝0的距离计算表达式为：
[0036][0037]
接着以最优解集p0相邻域内的数据点分别训练转矩、成本、转矩波动、效率、绕组平均温度和永磁体平均温度的局部代理模型k1，基于k1完成新一轮电机多目标鲁棒性优化，相应最优解集为p1，比较p0和p1的超体积参数h0和h1。其中超体积参数h在二维平面的几何含义如图4所示，计算方法为式(4)。
[0038][0039]
其中，hn为第n代局部代理模型kn的超体积指标，m为最优解集pn中的设计数目，xi和yi分别第i个最优设计的横纵坐标。
[0040]
相邻两次优化的超体积参数差别将作为优化迭代是否终止的条件，该指标的计算表达式为：
[0041][0042]
在该示例优化问题中，经过了1次全局代理模型优化和6次局部代理模型优化，收敛指标超体积差别ε低于1％，达到了收敛要求，其收敛过程见图5。
[0043]
为了验证局部代理模型(ki，i≥1)在小区间内的预测精度高于全局代理模型k0，抽取了500组不确定性数据，比较了代理模型对转矩、转矩波动、损耗、绕组平均温度和永磁体平均温度的预测能力，比较结果见表1。其中，采用取值范围为0～1的决定系数r2作为精度指标，r2越接近1表明预测能力越好。可以看出局部代理模型精度较全局代理模型有了很大提升。
[0044]
表1代理模型精度指标比较表
[0045]
r2转矩转矩波动损耗绕组温度永磁体温度k00.99370.92100.99740.93700.9834k10.99960.96150.99900.99500.9926k20.99940.98730.99940.99050.9891k30.99970.98580.99950.98700.9825
k40.99970.98060.99980.98600.9915k50.99970.98220.99940.98790.9871k60.99950.94180.99920.97910.9825
[0046]
全局代理模型k0和局部代理k6模型的最优解集如图6所示。通过有限元法验证最优解集p0和p6的误差分布分别如图7和图8所示。可见，采用了本发明所提出的策略之后，鲁棒性优化的最大误差由36％下降到1.1％，优化精度得到了显著提升。
[0047]
以上实施方式只是对本发明的示例性说明，并不限定它的保护范围，本领域技术人员还可以对其局部进行改变，例如，可以依据实际情况选择待优化的电机类型和电机性能，可以根据实际情况设计电机参数多目标优化模型，可以采用其它寻优算法优化代理模型，可以采用其它统计模拟法评价候选设计的可靠性，符合发明宗旨的任意形式的等同替换都落入本发明的保护范围。