CN113742980B - 双斜槽感应电机多物理场模型的优化设计方法 - Google Patents
双斜槽感应电机多物理场模型的优化设计方法 Download PDFInfo
- Publication number
- CN113742980B CN113742980B CN202111097471.5A CN202111097471A CN113742980B CN 113742980 B CN113742980 B CN 113742980B CN 202111097471 A CN202111097471 A CN 202111097471A CN 113742980 B CN113742980 B CN 113742980B
- Authority
- CN
- China
- Prior art keywords
- motor
- double
- model
- rotor
- chute
- 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
- 230000006698 induction Effects 0.000 title claims abstract description 31
- 238000000034 method Methods 0.000 title claims abstract description 19
- 238000004088 simulation Methods 0.000 claims abstract description 20
- 230000017525 heat dissipation Effects 0.000 claims abstract description 17
- 230000008878 coupling Effects 0.000 claims abstract description 7
- 238000010168 coupling process Methods 0.000 claims abstract description 7
- 238000005859 coupling reaction Methods 0.000 claims abstract description 7
- 230000002457 bidirectional effect Effects 0.000 claims abstract description 3
- 230000035882 stress Effects 0.000 claims description 16
- 230000008646 thermal stress Effects 0.000 claims description 14
- 241000555745 Sciuridae Species 0.000 claims description 5
- 230000004323 axial length Effects 0.000 claims 1
- 230000006641 stabilisation Effects 0.000 claims 1
- 238000011105 stabilization Methods 0.000 claims 1
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 4
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 4
- 229910052802 copper Inorganic materials 0.000 description 4
- 239000010949 copper Substances 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000005299 abrasion Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 229910052742 iron Inorganic materials 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000004804 winding Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 239000004020 conductor Substances 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 230000020169 heat generation Effects 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000000087 stabilizing effect Effects 0.000 description 1
- 238000009423 ventilation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Manufacture Of Motors, Generators (AREA)
Abstract
本发明公开了双斜槽感应电机多物理场模型的优化设计方法,包括根据已经确定的电机内外圆及长度,构建几何模型;初步构建中环参数化模型,根据电机模型的接触面散热系数和电机模型的电机内部损耗对电机几何模型进行温度场‑应力场耦合仿真,应力场数据对电磁性能进行双向耦合,选取最优的双斜槽中环厚度。本发明提出的双斜槽感应电机多物理场模型的优化设计方法,确定最优中环厚度,实现双斜槽电机的温度场‑应力场模拟并结合模型选取电磁方案最优点。
Description
技术领域
本发明涉及电机技术领域,特别是涉及双斜槽感应电机多物理场模型的优化设计方法。
背景技术
双斜槽转子的感应电机相比于单斜槽转子的感应电机具有诸多优点,其具体的转子结构可以分为上下层端环、上下层导条,中间环,通过设置合适的上下层导条的倾斜角度与相互交错的角度可以有效的削弱感应电机的转子上的轴向力以及特定阶次的气隙谐波从而有效的降低感应电机的振动与噪声。中间环的引入也可以有效降低由于感应电机转子斜槽所带来的横向电流引起的附加损耗,但同时,中间环的设计也对电磁方案带来了新的影响:
1.因为中环会导致气隙等效长度变长,励磁电抗降低,所以为了削弱中环对电机性能的影响,中环应尽可能的薄。但是过薄的中环将会因电密大而导致温度高,进而因热应力过大而存在于导条连接处产生断条风险。
2.若中环太厚虽避免了热应力的问题,但将使电磁性能下降。中环已导致气隙磁导轴向分布不均,过厚的中环会激发出更多次的谐波。
发明内容
本发明专利旨在提双斜槽感应电机多物理场模型的优化设计方法,该方法比传统电磁方案多考虑了温度场与应力场,将温度与热应力也作为电机性能的一部分对电磁方案进行考核,进而选取优化后的电磁方案。
本发明的目的通过以下技术方案实现:
双斜槽感应电机多物理场模型的优化设计方法,包括如下步骤:
步骤10:根据双斜槽感应电机机座号确定电机外壳内径D、长度L等参数,构建几何模型,构建电机几何模型,该电机几何模型包括机壳、定子、双斜槽转子、轴、转子端部和机尾风扇,并且本模型机壳考虑了接线盒。其中,双斜槽转子包括中上转子部分、中间环和下转子部分之间,中间环固定在上转子部分和下转子之间;
步骤20:仿真得到电机额定工况下的电磁性能,根据电磁性能仿真电机的电机内部损耗;
步骤30:根据步骤10中的几何模型计算温度场中接触面散热系数;
步骤40:初步构建中环参数化模型,并根据步骤30计算的接触面散热系数和步骤20得出的电机内部损耗对电机模型进行温度场-应力场耦合仿真,计及应力场数据对电磁性能进行双向耦合,选取最优的双斜槽中环厚度。
作为本发明进一步优化的方案,电机几何模型包括接线盒和机壳,进而使本方法确定的增大中间环参数确定的准确性。
作为本发明进一步优化的方案,其中所述步骤30中接触面散热系数包括机壳表面换热系数、转子端环表面换热系数、转子铁心端面换热系数、转轴表面换热系数。
作为本发明进一步优化的方案,其中所述步骤30中各接触面散热系数计算方法为:
(1)机壳表面换热系数采用经验公式:
α1=9.73+14v0.62
其中,v为机壳表面散热翅风速,单位为m/s;
(2)转子端环表面换热系数采用如下公式计算:
其中:λα为空气导热系数,单位为W/m·K;D2为转子直径,单位为m;n为转子转速,单位为r/min;hrf为风叶高度,单位为m;γ为空气运动黏度,单位为m2/s;
(3)转子铁心端面换热系数,其换热系数计算公式如下:
其中:Nur为转子铁心端面努塞尔特常数;Rer为转子铁心端面气流雷诺数;
(4)转轴表面换热系数,其换热系数计算公式如下:
α4=7.68(0.75ω)0.78
此处,ω为转子旋转角速度。
作为本发明进一步优化的方案,根据已经设置好的温度场数据,得到一组参数化应力场分布并采集鼠笼热应力最大点数据,取应力变化趋于稳定点为最优中间环厚度。
本发明公开的双斜槽感应电机多物理场模型的优化设计方法,因为以往的双斜槽电机电磁方案仅就电磁进行考虑,电磁方案最优但热与热应力方案不一定最优,根据感应电机产热原理,采用仿真进行电机稳态温度场和热应力场模拟,综合考虑最优中环厚度参数,实现双斜槽感应电机的温度场与热应力场模拟并优化散热。本方法适应性强,可优化在双斜槽感应电机不同工作点下的温度场分析,保证了电机关键零部件温度均匀,降低电机最高温度,优化电机工作环境。
附图说明
图1是双斜槽感应电机多物理场模型的优化设计方法流程图;
图2是双斜槽感应电机多物理场模型仿真模型图;
图3是双斜槽感应电机多物理场模型的优化设计方法程序框图;
图4是双斜槽感应电机温度场模型各接触面换热系数标识图;
图5是未优化双斜槽感应电机应力场分布云图;
图6是未优化双斜槽感应电机最大应力分布局部放大云图;
图7是中环厚度参数化优化图。
图中标号:步骤10:根据已经确定的内外圆及长度,构建几何模型进行电磁仿真;步骤20:仿真得到电磁方案数据,为将损耗作为载荷与温度场耦合仿真;步骤30:计算电机内部各接触面散热系数;步骤40:构建中环参数化模型,并进行电机发热稳态温度场仿真,根据电机温度场分布进行热应力场分布,根据温度场和应力场分布确定最优中环厚度;
201定子;202双斜槽转子;203转子端部;204转子转轴;205机尾风扇;206定子绕组;207接线盒;208机壳;
S1机壳表面接触面;S2转子端环表面接触面;S3转子铁心端面接触面;S4转轴表面接触面;
具体实施方式
下面详细描述本发明的实施例,所述实施例的示例在附图中表示,其中自始至终相同或类似的符号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,仅用于解释本发明,而不能理解对本发明的限制。
为使本发明的目的、技术方案和优点更加清楚,下面将结合实例及附图作进一步详细的描述。
如图1所示,双斜槽感应电机多物理场模型的优化设计方法,包括:
步骤10:根据双斜槽感应电机机座号确定电机外壳内径D、长度L,构建几何模型;
如图2所示为双斜槽感应电机的模型,模型由定子201与双斜槽转子202构成,其中双斜槽转子选择了最优的方案即斜一个定子槽距且上下笼错开半个齿距,为考虑更为精确的温度分布,考虑了转子端部203与轴204,机尾风扇205,定子绕组206则不考虑端部排列,该模型也考虑了接线盒207与机壳208,双斜槽转子202包括中上转子部分、中间环和下转子部分之间,中间环固定在上转子部分和下转子之间。
表1为一台25kW的四极电机主要参数
额定电压(V) | 230 |
额定频率(Hz) | 118 |
定子外径(mm) | 260 |
定子内径(mm) | 170 |
转子内径(mm) | 60 |
每槽导体数 | 16 |
中环厚度(mm) | 2 |
表1
步骤20:有限元仿真得到电机额定工况下的电磁性能,根据电磁性能有限元仿真电机的电机内部损耗电机内部损耗包括铁耗、定子铜耗、转子铜耗、杂散损耗和风磨耗。
如图3所示,上述步骤20中各损耗由基础模型有限元仿真得到;
表2为该电机满载运行时通过有限元仿真得到的电机内部损耗;
铁耗(W) | 890.32 |
定子铜耗(W) | 442.26 |
转子铜耗(W) | 324.69 |
杂散损耗(W) | 400 |
风磨耗(W) | 304 |
表2
步骤30:计算电机内部各接触面散热系数,为温度场与热应力场仿真做铺垫;
如图4所示,上述步骤30电机内部散热系数计算方法为:
1.机壳表面S1换热系数
电机采用机壳表面冷却通风系统,机壳表面的换热系数主要与机壳表面流体介质流速有关,采用经验公式:
α1=9.73+14v0.62
其中,v为机壳表面散热翅风速,m/s,本实施例中的机壳表面换热系数为78.371W/(m2·K)
2.转子端环表面S2换热系数
转子端环通过端环端面和风叶将导条产生的热量散向端腔内的空气,其换热系数常用如下公式计算:
其中,λα为空气导热系数,单位为W/m·K;D2为转子直径,单位为m;n为转子转速,单位为r/min;hrf为风叶高度,单位为m;γ为空气运动黏度,单位为m2/s。
本实施例中的转子端环表面换热系数为73.262W/(m2·K)。
3.转子铁心端面S3换热系数
转子铁心端面包括端环和转轴之间的环形面积,一部分热量通过此面向机内空气散出,其换热系数计算公式如下:
其中,Nur为转子铁心端面努塞尔特常数,Rer为转子铁心端面气流雷诺数。本实施例中的转子铁心端面换热系数为69.495W/(m2·K)。
4.转轴表面S4换热系数
转轴处几乎没有轴向气流影响,此处散热系数可按气流掠过大平面壁的相对运动来计算,又因旋转体表面气流存在着离心力和哥里奥里加速力影响,利用公式计算时,应按圆周速的75%来计算:其计算公式为:
α4=7.68(0.75ω)0.78,
本实施例中的转轴表面换热系数为:64.117W/(m2·K)。
步骤40:构建中环厚度参数化模型,根据现有的对中间环的外径和内径进行初步设计,对于此模型中环外径选择169.1mm、内径选择107.1mm,将上述步骤20与步骤30的计算结果导入温度场对步骤10的电机模型进行仿真,然后将温度场数据导入应力场并继续对电机模型进行仿真,拟固定中环外径与内径不变,根据电机应力场分布确定最优中环厚度。
将温度场数据导入应力场仿真中,将电机外壳与转轴两端设为固定支撑,即可得到热应力分布图,如图5所示,双斜槽感应电机的热应力分布图。
图6可以看出图5中A处由温升图得到的热应力图左可以看出鼠笼最大应力点出现在中环与鼠笼条连接处。
参数化中环厚度,根据已经设置好的温度场数据,将温度场数据导入应力场仿真,将得到一组参数化应力场分布并采集鼠笼热应力最大点数据,如图6所示。从图中可以看出,在电磁性能基本不受影响的前提下,对该模型来说,1.9mm为最佳中环厚度,中环再厚热应力对中环厚度就不再敏感。
虽然本发明所揭露的实施方法如上,但所述的内容只是为了便于理解本发明的实施方式,并非用于限定本发明。任何本发明所属技术领域内的技术人员,在不脱离本发明所揭露的精神和范围前提下,可以在实施的形式上及细节上作任何修改与变化,但本发明的专利保护范围,仍需以所附的权利要求书所界定的范围为准。
Claims (3)
1.双斜槽感应电机多物理场模型的优化设计方法,其特征在于,包括如下步骤:
步骤10:根据确定的内外圆及轴向长度,构建电机几何模型;
步骤20:仿真得到电机额定工况下的电磁性能,根据电磁性能仿真电机的电机内部损耗;
步骤30:根据步骤10中的几何模型计算温度场中接触面散热系数,接触面散热系数包括机壳表面换热系数、转子端环表面换热系数、转子铁心端面换热系数、转轴表面换热系数;
各接触面散热系数计算方法为:
(1)机壳表面换热系数采用经验公式:
α1=9.73+14V0.62
其中,v为机壳表面散热翅风速;
(2)转子端环表面换热系数采用用如下公式计算:
其中:λα为空气导热系数;D2为转子直径;n为转子转速;hrf为风叶高度;γ为空气运动黏度;
(3)转子铁心端面换热系数,其换热系数计算公式如下:
其中:Nurf为转子铁心端面努塞尔特常数;Rer为转子铁心端面气流雷诺数;
(4)转轴表面换热系数,其换热系数计算公式如下:
α4=7.68(0.75ω)0.78
其中:ω为转子旋转角速度;
步骤40:初步构建中环参数化模型,根据步骤30得出的接触面散热系数和步骤20得出的电机内部损耗对电机几何模型进行温度场-应力场耦合仿真,应力场数据对电磁性能进行双向耦合,选取最优的双斜槽中环厚度。
2.根据权利要求1所述的双斜槽感应电机多物理场模型的优化设计方法,其特征在于,电机几何模型包括接线盒和机壳。
3.如权利要求1中所述的双斜槽感应电机多物理场模型的优化设计方法,其特征在于,所述步骤40中参数化中环最优厚度的确定为:根据已经设置好的温度场数据,得到一组参数化应力场分布并采集鼠笼热应力最大点数据,取应力变化趋于稳定点为最优中间环厚度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111097471.5A CN113742980B (zh) | 2021-09-18 | 2021-09-18 | 双斜槽感应电机多物理场模型的优化设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111097471.5A CN113742980B (zh) | 2021-09-18 | 2021-09-18 | 双斜槽感应电机多物理场模型的优化设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113742980A CN113742980A (zh) | 2021-12-03 |
CN113742980B true CN113742980B (zh) | 2024-03-29 |
Family
ID=78739900
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111097471.5A Active CN113742980B (zh) | 2021-09-18 | 2021-09-18 | 双斜槽感应电机多物理场模型的优化设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113742980B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010105471A1 (zh) * | 2009-03-20 | 2010-09-23 | 西安电子科技大学 | 基于机电热三场耦合的电子设备机箱结构优化设计方法 |
WO2020056405A1 (en) * | 2018-09-14 | 2020-03-19 | Northwestern University | Data-driven representation and clustering discretization method and system for design optimization and/or performance prediction of material systems and applications of same |
CN112134410A (zh) * | 2020-09-28 | 2020-12-25 | 揭阳市汇宝昌电器有限公司 | 一种永磁电机的稳态温度场模型的优化散热设计方法 |
CN112528540A (zh) * | 2020-12-08 | 2021-03-19 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种基于力热耦合计算的结构优化方法 |
CN112800685A (zh) * | 2021-03-26 | 2021-05-14 | 湖南大学 | 一种高精度的电机多物理场耦合仿真计算方法 |
-
2021
- 2021-09-18 CN CN202111097471.5A patent/CN113742980B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010105471A1 (zh) * | 2009-03-20 | 2010-09-23 | 西安电子科技大学 | 基于机电热三场耦合的电子设备机箱结构优化设计方法 |
WO2020056405A1 (en) * | 2018-09-14 | 2020-03-19 | Northwestern University | Data-driven representation and clustering discretization method and system for design optimization and/or performance prediction of material systems and applications of same |
CN112134410A (zh) * | 2020-09-28 | 2020-12-25 | 揭阳市汇宝昌电器有限公司 | 一种永磁电机的稳态温度场模型的优化散热设计方法 |
CN112528540A (zh) * | 2020-12-08 | 2021-03-19 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种基于力热耦合计算的结构优化方法 |
CN112800685A (zh) * | 2021-03-26 | 2021-05-14 | 湖南大学 | 一种高精度的电机多物理场耦合仿真计算方法 |
Non-Patent Citations (2)
Title |
---|
基于磁热耦合法的电机模型参数计算;吴尧辉;方鑫;付豪;;电机与控制应用(06);全文 * |
防爆型三相异步电机中高精度多物理场耦合模型的研究和分析;程寿国;赵勇;陈虎威;王荷芬;;计量与测试技术(01);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113742980A (zh) | 2021-12-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ulbrich et al. | Cooling fin optimization on a TEFC electrical machine housing using a 2-D conjugate heat transfer model | |
JP5358667B2 (ja) | 永久磁石式発電機 | |
US10072672B2 (en) | Fan | |
JP3279258B2 (ja) | 車両用交流発電機 | |
Zhang et al. | Electrothermal combined optimization on notch in air-cooled high-speed permanent-magnet generator | |
CN108111079B (zh) | 基于永磁电机转子分段护套涡流损耗的热传递比计算方法 | |
Valenzuela et al. | Heat transfer and thermal design of finned frames for TEFC variable-speed motors | |
CN207320974U (zh) | 一种含非对称结构转子铁芯的自启动同步磁阻电机 | |
CN112134410B (zh) | 一种永磁电机的稳态温度场模型的优化散热设计方法 | |
CN104659964A (zh) | 交流发电机 | |
Moon et al. | Numerical investigation on thermal-flow characteristics of a totally enclosed fan cooled induction motor | |
Choi et al. | Enhanced efficiency of the brushless direct current motor by introducing air flow for cooling | |
Fang et al. | Effects of eddy current in electrical connection surface of laminated cores on high-speed PM motor supported by active magnetic bearings | |
US11218056B2 (en) | Main electric-motor for vehicle | |
CN113742980B (zh) | 双斜槽感应电机多物理场模型的优化设计方法 | |
Xu et al. | Global fluid flow and heat transfer characteristics analysis of an open air-cooled drive motor for drilling application | |
Kepsu et al. | Modeling of a 30 000 Rpm bearingless SPM drive with loss and thermal analyses for a 0.5 MW high-temperature heat pump | |
Priyanka et al. | Thermal Analysis of Multiphase Induction Motor for Electric Vehicle Applications | |
Vu et al. | New cooling system design of bldc motor for electric vehicle using computation fluid dynamics modeling | |
Kimotho | Optimizing heat sink geometry for electric vehicle BLDC motor using CFD | |
CN113890234B (zh) | 一种具有转子孔间循环对流的封闭式电机冷却装置 | |
CN215934547U (zh) | 电机和电器设备 | |
JP4640681B2 (ja) | 回転電機 | |
JP2017050913A (ja) | 回転電機 | |
JPH0698511A (ja) | 車両用交流発電機 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |