一种考虑电磁与力学性能的电机转子结构拓扑优化方法技术

本发明专利技术公开了一种考虑电磁与力学性能的电机转子结构拓扑优化方法，该考虑电磁与力学性能的电机转子结构拓扑优化方法，对转子结构上每个单元材料作为设计变量，来决定材料属性，如果是两种材料用离散的方式表示，即空气和铁磁材料可以用0和1表示。通过SIMP方法将离散问题转化为连续问题，定义[0,1]范围内变化的连续变量，建立考虑电机转子力学性能约束的目标方程，操作简单，首先只需设定转子大小即可，例如转子的内径和外径。然后通过拓扑优化算法得到最优拓扑结构，进行不断迭代优化的有限元计算，前期结构设计无需人为的计算，在电机设计结构上，大大节省时间和试验成本。大大节省时间和试验成本。大大节省时间和试验成本。

【技术实现步骤摘要】
一种考虑电磁与力学性能的电机转子结构拓扑优化方法


[0001]本专利技术涉及一种电机转子结构设计
，具体是一种考虑电磁与力学性能的电机转子结构拓扑优化方法。

技术介绍

[0002]早期历史上的电机结构设计同其它传统设计一样是经验和半经验的产品，传统的结构设计，一般是在经验数据的基础上按照几何作图法来不断地调整结构参数，设计周期比较长，而且结果比较粗糙。随着科技进步和社会发展，以及资源短缺问题的突出，致使单纯依靠工程设计人员的经验和模拟实验的传统优化设计方式己难以胜任。一般来说：电机结构优化设计由以下步骤：初始计算粗糙的几何结构模型(经验法)建立形状目标优化目标函数方程有限元电磁计算判断是否满足最优解，未满足则返回第二步，直至迭代满足最优解为止。一般来说：电机结构优化设计由以下步骤：初始计算粗糙的几何结构模型(经验法)
→
建立形状目标优化目标函数方程
→
有限元电磁计算
→
判断是否满足最优解，未满足则返回第二步，直至迭代满足最优解为止。
[0003]在形状优化设计中，电机的初始结构由设计者给定，通过电机结构的尺寸参数为设计变量，以达到目标函数最优的目的。但是，对经验丰富的工程师所设计的结构，仅仅通过尺寸结构的改善是很难实现对原设计的较大优化，无法获得具有创新性的结构形式，因此提出一种考虑电磁与力学性能的电机转子结构拓扑优化方法，以克服经验法无法获得具有创新性的结构形式。

技术实现思路

[0004]本专利技术的目的在于提供一种考虑电磁与力学性能的电机转子结构拓扑优化方法，以解决上述
技术介绍
中提出的问题。
[0005]为实现上述目的，本专利技术提供如下技术方案：
[0006]一种考虑电磁与力学性能的电机转子结构拓扑优化方法，所述方法步骤如下：
[0007]步骤一：建立电机几何模型，对转子设计域内进行剖分；
[0008]步骤二：转子设计域内进行初始的材料布局，考虑电磁性能方面，在电磁计算中，计算铁磁材料和空气的区别，定义设计域中的材料分布，建立应用材料插值函数，插值函数为：ν＝ν
fe
x
p
+ν
air
(1
‑
x
p
)，考虑力学性能方面，在应力计算中，计算铁磁材料和空气的区别，定义设计域中的材料分布，建立应用材料插值函数，插值函数为：E＝E
fe
x
p
，w＝w
fe
x
p
；
[0009]步骤三：进行电机电磁场的控制方程计算，计算公式为：SA
‑
F＝0；
[0010]步骤四：进行线性弹性力学的控制方程计算，计算公式为：KU
‑
f＝0；
[0011]步骤五：构建电机转矩为优化对象，离心应力为约束的控制方程组，方程组为：
[0012][0013][0014]g(x)＝Γ(Γ2+1)≤0
[0015][0016][0017]步骤六：增广拉格朗日方程框架的建立，计算公式为：
[0018][0019][0020][0021]步骤七：对于基于MMA梯度的优化算法的应用，增广拉格朗日方程对设计变量的敏度的计算，计算公式为：
[0022][0023]步骤八：通过MMA梯度算法求解增广拉格朗日方程最小值，计算公式为：
[0024]μ
(k+1)
＝min[αμ
(k)
,μ
max
][0025][0026]以达到最优目的；
[0027]步骤九：迭代最优为止，获取拓扑结构形状。
[0028]作为本专利技术进一步的方案：所述步骤二中电磁计算中铁磁材料和空气的区别在于磁导率或磁阻率。
[0029]作为本专利技术再进一步的方案：所述步骤二中应力计算中铁磁材料和空气的区别在于杨氏模量。
[0030]与现有技术相比，本专利技术的有益效果是：
[0031]本专利技术操作简单，首先只需设定转子大小即可，例如转子的内径和外径。然后通过拓扑优化算法得到最优拓扑结构，前期结构设计无需人为的计算，在电机设计结构上，大大节省时间和试验成本。
附图说明
[0032]图1为一种考虑电磁与力学性能的电机转子结构拓扑优化方法的流程图。
[0033]图2为一种考虑电磁与力学性能的电机转子结构拓扑优化方法的效果图。
具体实施方式
[0034]下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。
[0035]请参阅图1～2，本专利技术实施例中，一种考虑电磁与力学性能的电机转子结构拓扑
优化方法，对转子结构上每个单元材料作为设计变量，来决定材料属性，如果是两种材料用离散的方式表示，即空气和铁磁材料可以用0和1表示。通过SIMP方法将离散问题转化为连续问题，定义[0,1]范围内变化的连续变量，建立考虑电机转子力学性能约束的目标方程，进行不断迭代优化的有限元计算，方法步骤如下：
[0036]步骤一：建立电机几何模型，对转子设计域内进行剖分；
[0037]步骤二：转子设计域内进行初始的材料布局，考虑电磁性能方面，在电磁计算中，计算铁磁材料和空气的区别，定义设计域中的材料分布，建立应用材料插值函数，插值函数为：ν＝ν
fe
x
p
+ν
air
(1
‑
x
p
)，考虑力学性能方面，在应力计算中，计算铁磁材料和空气的区别，定义设计域中的材料分布，建立应用材料插值函数，插值函数为：E＝E
fe
x
p
，w＝w
fe
x
p
；
[0038]步骤三：进行电机电磁场的控制方程计算，计算公式为：SA
‑
F＝0；
[0039]步骤四：进行线性弹性力学的控制方程计算，计算公式为：KU
‑
f＝0；
[0040]步骤五：构建电机转矩为优化对象，离心应力为约束的控制方程组，方程组为：
[0041][0042][0043]g(x)＝Γ(Γ2+1)≤0
[0044][0045][0046]步骤六：增广拉格朗日方程框架的建立，计算公式为：
[0047][0048][0049][0050]步骤七：对于基于MMA梯度的优化算法的应用，增广拉格朗日方程对设计变量的敏度的计算，计算公式为：
[0051][0052]步骤八：通过MMA梯度算法求解增广拉格朗日方程最小值，计算公式为：
[0053]μ
(k+1)
＝min[αμ
(k)
,μ
max
][0054][00本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种考虑电磁与力学性能的电机转子结构拓扑优化方法，其特征在于：所述方法步骤如下：步骤一：建立电机几何模型；步骤二：转子设计域内进行初始的材料布局，考虑电磁性能方面，计算铁磁材料和空气的区别，定义设计域中的材料分布，建立应用材料插值函数，插值函数为：ν＝ν
fe
x
p
+ν
air
(1
‑
x
p
)，考虑力学性能方面，计算铁磁材料和空气的区别，定义设计域中的材料分布，建立应用材料插值函数，插值函数为：E＝E
fe
x
p
，w＝w
fe
x
p
；步骤三：进行电机电磁场的控制方程计算，计算公式为：SA
‑
F＝0；步骤四：进行线性弹性力学的控制方程计算，计算公式为：KU
‑
...

【专利技术属性】
技术研发人员：马伯，李宇，郑静，杨江涛，黄守道，吴公平，
申请(专利权)人：湖南大学，
类型：发明
国别省市：