【技术实现步骤摘要】
本专利技术隶属于电机拓扑优化设计领域，具体涉及一种考虑空载情况下的基于变密度法的电机拓扑优化方法。
技术介绍
1、拓扑优化作为最高层次的优化方法，突破了传统参数优化的设计理念，通过对有限设计域内有限材料的充分利用及合理布局，对电机结构进行创新构型设计，在保障结构实现高效率、高性能等要求的同时满足结构轻量化的要求。针对电机优化设计领域存在的复杂工况、多约束、高维度的问题，电机拓扑优化方法主要包括离散方法、变密度方法(simp)、水平集方法等。因变密度方法可有效地解决超高维度设计变量的数学规划问题被广泛的应用，该方法的核心思想是将离散设计变量松弛为[0,1]区间的连续设计变量，通过构建虚拟材料密度的非线性插值模型来描述材料属性，使得其转化为连续性优化问题。
2、在多物理场的复杂工况条件下，电机优化设计的轻量化、高性能、高效率等方面的性能指标对基于simp的电机拓扑优化提出了更高的要求。正因电机结构的多种高性能指标，基于simp的电机拓扑优化通常是一个多目标问题，需要对优化目标与约束条件进行有效处理与组合。一般地，基于simp的电机拓扑优化主要是计算电机带载性能，但因输入的正弦电流将会增加其复杂性，导致收敛速度缓慢。并且，优化后的电机拓扑构型具有锯齿状边缘，导致实际制造难度剧增。
技术实现思路
1、为解决综合考虑带载及空载条件下永磁电机拓扑优化的有限元效率问题，本专利技术提供一种考虑空载情况下的基于变密度法的电机拓扑优化方法。首先，确定转子结构设计域和有限元模型，定义设计变量
2、本专利技术提供了一种考虑空载情况下的基于变密度法的电机拓扑优化方法，本专利技术包括如下步骤：
3、步骤1、构建ipmsm几何模型，以1/8的ipmsm转子作为设计域，并利用有限元进行网格剖分；
4、步骤2、定义单元密度ρi为设计变量，即对于各向同性材料，引入单元密度作为设计变量，其变化区间为(0,1]；
5、步骤3、对单元密度ρi进行helmholtz过滤和heaviside投影过滤后得到单元物理密度。helmholtz过滤方程可表示为：
6、
7、其中，为中间密度，ni,j为单元i在过滤半径内的单元集合，ρj为前述集合中的第j个单元的密度，wi,j为权函数可表示为：
8、wi,j＝max(0,r-ri,j)q (2)
9、其中，r为过滤半径，ri,j为单元i质心到单元j质心的距离，q为过滤指数。经过helmholtz过滤后，对中间密度进行heaviside投影过滤后得到单元物理密度，heaviside投影过滤方程可表示为：
10、
11、其中，为单元物理密度，β为投影参数，η为过滤密度的阈值；
12、步骤4、基于simp密度插值模型，构造电磁与力学耦合的材料特性插值模型，其中电磁特性包括磁导率，力学特性包括杨氏模量和物理质量。材料密度与磁导率的插值关系可表示为：
13、
14、其中，ν为磁导率，νfe和νair分别为铁磁材料和空气的磁导率，p为惩罚因子。材料密度与杨氏模量和物理质量的插值关系可表示为：
15、
16、其中，e为杨氏模量，efe为铁磁材料的杨氏模量，m为物理质量，mfe为铁磁材料物理质量向量，⊙为哈达玛积；
17、步骤5、基于电磁与力学耦合的材料特性插值模型，进行有限元分析。电磁场的非线性有限元分析可表示为：
18、sa＝q (7)
19、其中，s为非线性材料的刚度矩阵，a为未知节点矢量磁位，q为激励源向量。静力学场的有限元分析可表示为：
20、ku＝f (8)
21、其中，k为静力学刚度矩阵，u为未知位移向量，f为载荷向量；
22、步骤6、考虑空载情况下，构造关于齿槽转矩的方差为目标函数。利用麦克斯韦张量法计算齿槽转矩，计算公式可表示为：
23、
24、其中，ti为第i个转子位置的齿槽转矩，θelec,i为第i个位置的转子电气角位移，r为气隙圆周半径，l为电机轴向长度，μ0为真空磁导率，br和bt分别为径向和切向磁密，sag为气隙弧度。
25、以齿槽转矩的方差的目标函数可表示为：
26、
27、其中，tvr为齿槽转矩方差，n为转子旋转位置的数量，tav为平均转矩；
28、步骤7、利用离散傅里叶变换技术，构造空载情况下关于磁链的基波和谐波的约束条件。基于有限元分析的结果，每一项绕组磁链的计算公式可表示为：
29、
30、其中，ψk为第k项绕组磁链，θelec为任意一个转子电气角位移，nz为线圈边数，ns为第s个线圈边的导体数，azs为第s个线圈边所在区域的矢量磁位，ds为±1，表示第s个线圈边电流的方向。通过离散傅里叶变换技术，将第k项绕组磁链用不同频次h的磁链谐波之和来表示，转换公式可表示为：
31、
32、其中，h为磁链谐波频次，表示为h＝(0,1,...,m-1)。如果转子旋转位置的数量n为偶数或n＝2m，磁链谐波振幅和可表示为：
33、
34、磁链谐波的幅值可表示为：
35、
36、经过离散傅里叶变换技术后，根据谐波的基础定义，即频次h＝1的波形则为磁链基波，频次h≥2的波形则为高次磁链谐波，构造以磁链基波幅值和总谐波失真率的约束条件，可分别表示为：
37、
38、步骤8、利用冻结磁导率技术，构造空载情况下关于凸极比的约束条件。通过冻结磁导率技术获得电机dq轴电感表征电机凸极比，计算公式可表示为：
39、s＝lq/ld (19)
40、其中，s为凸极比，lq为q轴电感，ld为d轴电感；
41、步骤9、结合以上构造的目标函数和约束条件，构建空载情况下的拓扑优化模型。以齿槽转矩方差为目标函数，应力、质量、柔度、磁链基波幅值、总谐波失真率和凸极比为约束函数的一般拓扑优化模型可表示为：
42、
43、其中，ne为设计域内单元数量，应力、质量、柔度、磁链基波幅值、总谐波失真率和凸极比的约束函数可表示为：
44、
45、其中，为第j个单元的冯·米塞斯应力，σlim为材料的许用应力，i为单位向量，ar为质量百分比，c为转子结构柔度，c0为柔度参考值，ψf0为磁链基波幅值参考值，thd0为总谐波失真率参考值，s0为凸极比参考值。
46、基于一般拓扑优化模型，利用归一化拉格朗日方法构建增广拉格朗日优化框架，可表示为：<本文档来自技高网...
【技术保护点】
1.一种考虑空载情况下的基于变密度法的电机拓扑优化方法，其特征在于，所述方法步骤如下：
2.根据权利要求1所述的一种考虑空载情况下的基于变密度法的电机拓扑优化方法，其特征在于：步骤5中的电磁场的非线性有限元方程SA＝q，采用牛顿拉夫森迭代法求解；静力学场的有限元方程KU＝f中的单元载荷的求解方程可表示为：
3.根据权利要求2所述的一种考虑空载情况下的基于变密度法的电机拓扑优化方法，其特征在于：步骤6中的平均转矩Tav的求解方程可表示为：
4.根据权利要求3所述的一种考虑空载情况下的基于变密度法的电机拓扑优化方法，其特征在于：步骤8中的凸极比的求解过程中运用了基于电磁关系的绕组磁链求解模型和冻结磁导率技术求解dq轴电感，具体步骤包括：
5.根据权利要求4所述的一种考虑空载情况下的基于变密度法的电机拓扑优化方法，其特征在于：步骤9中的一般优化拓扑优化模型中的单元冯·米塞斯应力和转子结构柔度C的计算公式可表示为：
6.根据权利要求5所述的一种考虑空载情况下的基于变密度法的电机拓扑优化方法，其特征在于：步骤10中的空载情况下的拓
7.根据权利要求6所述的一种考虑空载情况下的基于变密度法的电机拓扑优化方法，其特征在于：步骤13中的水平集的分割法则，具体步骤包括：
...【技术特征摘要】
1.一种考虑空载情况下的基于变密度法的电机拓扑优化方法，其特征在于，所述方法步骤如下：
2.根据权利要求1所述的一种考虑空载情况下的基于变密度法的电机拓扑优化方法，其特征在于：步骤5中的电磁场的非线性有限元方程sa＝q，采用牛顿拉夫森迭代法求解；静力学场的有限元方程ku＝f中的单元载荷的求解方程可表示为：
3.根据权利要求2所述的一种考虑空载情况下的基于变密度法的电机拓扑优化方法，其特征在于：步骤6中的平均转矩tav的求解方程可表示为：
4.根据权利要求3所述的一种考虑空载情况下的基于变密度法的电机拓扑优化方法，其特征在于：步骤8中的凸极比的求解过...
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。