CN109753737B - 用于交流牵引电机温度场分析的定子绕组气隙建模方法 - Google Patents

用于交流牵引电机温度场分析的定子绕组气隙建模方法 Download PDF

Info

Publication number
CN109753737B
CN109753737B CN201910023474.0A CN201910023474A CN109753737B CN 109753737 B CN109753737 B CN 109753737B CN 201910023474 A CN201910023474 A CN 201910023474A CN 109753737 B CN109753737 B CN 109753737B
Authority
CN
China
Prior art keywords
air gap
stator winding
traction motor
current traction
stator
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
Application number
CN201910023474.0A
Other languages
English (en)
Other versions
CN109753737A (zh
Inventor
张小平
刘苹
匡斯建
张铸
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hunan University of Science and Technology
Original Assignee
Hunan University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hunan University of Science and Technology filed Critical Hunan University of Science and Technology
Priority to CN201910023474.0A priority Critical patent/CN109753737B/zh
Publication of CN109753737A publication Critical patent/CN109753737A/zh
Priority to PCT/CN2019/122279 priority patent/WO2020143361A1/zh
Application granted granted Critical
Publication of CN109753737B publication Critical patent/CN109753737B/zh
Priority to US17/037,873 priority patent/US11196374B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02K—DYNAMO-ELECTRIC MACHINES
    • H02K3/00—Details of windings
    • H02K3/04—Windings characterised by the conductor shape, form or construction, e.g. with bar conductors
    • H02K3/12—Windings characterised by the conductor shape, form or construction, e.g. with bar conductors arranged in slots
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • 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]
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02K—DYNAMO-ELECTRIC MACHINES
    • H02K3/00—Details of windings
    • H02K3/32—Windings characterised by the shape, form or construction of the insulation
    • H02K3/34—Windings characterised by the shape, form or construction of the insulation between conductors or between conductor and core, e.g. slot insulation
    • H02K3/345—Windings characterised by the shape, form or construction of the insulation between conductors or between conductor and core, e.g. slot insulation between conductor and core, e.g. slot insulation
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02K—DYNAMO-ELECTRIC MACHINES
    • H02K3/00—Details of windings
    • H02K3/46—Fastening of windings on the stator or rotor structure
    • H02K3/47—Air-gap windings, i.e. iron-free windings
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P29/00—Arrangements for regulating or controlling electric motors, appropriate for both AC and DC motors
    • H02P29/60—Controlling or determining the temperature of the motor or of the drive
    • H02P29/64—Controlling or determining the temperature of the winding
    • 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
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P23/00—Arrangements or methods for the control of AC motors characterised by a control method other than vector control
    • H02P23/14—Estimation or adaptation of motor parameters, e.g. rotor time constant, flux, speed, current or voltage

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Power Engineering (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Computer Hardware Design (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Manufacture Of Motors, Generators (AREA)

Abstract

本发明公开了一种用于交流牵引电机温度场分析的定子绕组气隙建模方法，是通过按间距设定值改变定子绕组等效模型的气隙宽度，建立带有定子绕组气隙的交流牵引电机的三维有限元模型，根据不同气隙宽度的三维有限元模型经温度场分析得交流牵引电机的温度场分布图，根据温度场分布图中数据进行数值拟合，得定子绕组等效模型气隙宽度与温度间函数关系，通过实测电机定子绕组与其等效模型表面某处对应的实际温度，计算当前交流牵引电机定子绕组建模所对应的最佳气隙宽度。与现有技术相比，采用本发明可建立交流牵引电机定子绕组气隙的准确模型，从而有效提高交流牵引电机有限元模型的准确性。

Description

用于交流牵引电机温度场分析的定子绕组气隙建模方法
技术领域
本发明涉及交流牵引电机的温度场分析领域，具体涉及一种用于交流牵引电机温度场分析的定子绕组气隙建模方法。
背景技术
交流牵引电机因具有结构简单、运行可靠、坚固耐用、功率大、转速高等系列优点而在众多领域得到了日益广泛的应用。交流牵引电机运行时因各种损耗而使其内部温度上升，如果温升过高则会严重影响电机的使用寿命。因此，针对电机内部的温度场展开分析，以指导电机结构的优化设计，对于降低电机温升具有重要意义。
现有技术中，由于采用有限元法对交流牵引电机进行温度场分析时，能准确反映出电机内部各单元的温度分布情况，分析结果准确，因而获得了广泛应用。但采用有限元法对交流牵引电机进行温度场分析时，对模型的准确性要求高；如果完全按照电机的实际结构来建立其有限元模型，则存在模型复杂、建模时间长、计算工作量大、对计算机性能要求高等诸多问题，因而难以推广应用。因此，在实际采用有限元法对交流牵引电机进行温度场分析时，往往对电机模型进行适当简化，尤其是对于定子绕组中存在的气隙，由于绕组中各导线间及相关导线与定子槽壁间气隙的不规则，因而使得绕组气隙的建模变得十分困难，故在实际建模时往往将其忽略；如此处理虽然简化了模型，但影响了温度场分析的准确性。
发明内容
针对现有技术存在的上述问题，本发明提供一种用于交流牵引电机温度场分析中原理简单、算法精度高、占用系统资源少的定子绕组气隙建模方法。
本发明提出的技术方案为：
本发明所提供的用于交流牵引电机温度场分析的定子绕组气隙建模方法，是通过按间距设定值改变定子绕组等效模型的气隙宽度，建立带有定子绕组气隙的交流牵引电机的三维有限元模型，根据不同气隙宽度的三维有限元模型经温度场分析得交流牵引电机的温度场分布图，根据温度场分布图中数据进行数值拟合，得定子绕组等效模型气隙宽度与温度间函数关系，通过实测电机定子绕组与其等效模型表面某处对应的实际温度，计算当前交流牵引电机定子绕组建模所对应的最佳气隙宽度。
是通过按间距设定值改变定子绕组等效模型的气隙宽度，建立带有定子绕组气隙的交流牵引电机的三维有限元模型，经温度场分析得交流牵引电机的温度场分布图，求解不同气隙宽度下交流牵引电机的温度场分布，根据上述定子绕组等效模型的气隙宽度与其相应的温度数据进行数值拟合，得定子绕组等效模型气隙宽度与温度间函数关系，通过实测电机定子绕组与其等效模型表面某处对应的实际温度，计算得与该实测温度对应的气隙宽度，该气隙宽度即作为当前交流牵引电机定子绕组建模所对应的最佳气隙宽度。
本发明所提供的用于交流牵引电机温度场分析的定子绕组气隙建模方法，包括以下步骤：
步骤1)根据交流牵引电机定子绕组的结构建立定子绕组等效模型；
步骤2)根据步骤1所得定子绕组等效模型确定该模型与定子铁芯槽壁及槽楔间隔离层的结构与宽度；
步骤3)根据步骤2所得隔离层的结构与宽度确定隔离层中对应绝缘层与气隙的结构与初始宽度，再根据电机其余部件的结构参数，建立交流牵引电机的三维有限元模型；
步骤4)根据步骤3所得交流牵引电机的三维有限元模型进行温度场分析，得到交流牵引电机的温度场分布图，并选取定子绕组等效模型表面某处的温度作为待测值；
步骤5)按间距设定值改变定子绕组等效模型的气隙宽度，再分别建立交流牵引电机的三维有限元模型，求解不同气隙宽度下交流牵引电机的温度场分布，并分别测取定子绕组等效模型表面与步骤4中某处对应的温度；
步骤6)针对步骤5所得N组定子绕组等效模型的气隙宽度与其相应的温度数据进行数值拟合，得到定子绕组等效模型气隙宽度与温度间的计算公式(8)：
h1(T)＝aebT+cedT (8)
式中：h1(T)为绕组气隙宽度函数；T为绕组等效模型表面某处的温度；a、b、c、d分别为气隙宽度函数的系数；e为自然对数函数的底数；具体地，a、b、c、d根据最小二乘法进行确定；
步骤7)实测电机定子绕组与其等效模型表面某处对应的实际温度，并将该实测温度代入上述计算公式(8)，得到与该实测温度对应的气隙宽度，该气隙宽度即作为当前交流牵引电机定子绕组建模所对应的最佳气隙宽度。
优选地，步骤1在建模过程中将定子绕组中多匝铜线等效为单匝铜线，单匝铜线的截面积与多匝铜线的截面积相等，且单匝铜线的截面形状与定子槽型相同，单匝铜线的截面积通过公式(1)求得：
s1＝nπr1 2 (1)
式中：s1为单匝铜线的截面积，r1为定子绕组多匝铜线中每匝铜线的半径，n为定子绕组多匝铜线的匝数。
优选地，步骤2中定子绕组等效模型与定子铁芯槽壁及槽楔间隔离层的宽度处处相等，所述隔离层的宽度通过公式(2)计算而得：
Figure BDA0001941629540000031
式中：h为隔离层的宽度，r为定子槽底半径，s2为定子槽截面面积。
优选地，步骤3中隔离层中对应气隙的初始宽度设定为定子槽内单层铜线的等效气隙宽度，具体为：
Figure BDA0001941629540000032
式中：h1为隔离层中气隙的初始宽度。
优选地，步骤3中隔离层中对应绝缘层的初始宽度为：
h2＝h-h1 (4)
式中：h2为隔离层中绝缘层的初始宽度。
优选地，步骤4的具体操作是：
步骤4-1)针对交流牵引电机三维有限元模型进行网格划分；
步骤4-2)针对网格划分后的交流牵引电机三维有限元模型施加热源；
步骤4-3)设置交流牵引电机三维有限元模型的边界条件和对流换热系数；
步骤4-4)针对三维有限元模型进行温度场有限元计算，从而获得交流牵引电机温度场分布图，并选取绕组等效模型表面某处的温度作为待测值。
更优选地，步骤4-3的具体操作为：
(a)机壳、散热筋与空气间的对流换热系数设置为处处相等；
(b)根据公式(5)和公式(6)计算定子和转子之间气隙的雷诺数Re及其临界雷诺数Rel；
Figure BDA0001941629540000033
Figure BDA0001941629540000034
式中：d1为转子半径，d2为定子半径，δ为定子和转子之间气隙的长度，ωg为转子转速，v为空气的运动粘度；
(c)根据步骤b所得雷诺数Re及其临界雷诺数Rel的大小确定相应的对流换热系数，将定子和转子之间气隙中对流换热系数的设置为计算所得数值，具体为：
当Re<Rel时，表明此时气隙内空气流动为层流，取对流换热系数α为空气的导热系数；
当Re>Rel时，表明此时气隙内空气流动为湍流，则对流换热系数α通过公式(7)计算而得：
Figure BDA0001941629540000041
式中：λ为空气的导热系数。
优选地，步骤5中的间距设定值，是指以初始气隙宽度为初始值，以隔离层宽度为最终值，根据数值拟合的需要确定间距的大小，从初始值开始按间距逐步增大气隙宽度，直至最终值为止。
与现有技术相比，本发明提供的用于交流牵引电机温度场分析的定子绕组气隙建模方法具有如下优点：
采用本发明提供的用于交流牵引电机温度场分析的定子绕组气隙建模方法，可建立交流牵引电机定子绕组气隙的准确模型，从而有效提高交流牵引电机有限元模型的准确性，进而可有效提高交流牵引电机温度场分析的准确性，为实现交流牵引电机结构的优化设计提供有益的指导。
附图说明
图1为本发明实施例提供的交流牵引电机定子槽内建模示意图；
图2为本发明实施例提供的根据定子槽模型建立的交流牵引电机三维有限元模型示意图；
图3为本发明实施例提供的用于交流牵引电机温度场分析的定子绕组气隙建模方法流程图；
图4为本发明实施例提供的用于交流牵引电机温度场分析的定子绕组气隙建模方法详细流程图；
图5为本发明实施例提供的交流牵引电机温度场分布示意图；
表1为本发明实施例提供的定子绕组等效模型1的气隙宽度与其相应的温度数据；
其中，1，定子绕组等效模型；2，绝缘层2；3，定子绕组等效模型表面绝缘层与定子铁芯槽壁及槽楔间的气隙；4，槽壁；5，槽楔；6，机壳；7，定子铁芯；8，定子楔部；9，转子导条；10，转子铁芯；11，定子和转子之间的气隙；12，转子端环；13，轴承；14，待测点。
具体实施方式
以下结合附图对本发明的具体实施方式进行详细说明。应当理解的是，此处所描述的具体实施方式仅用于说明和解释本发明，并不用于限制本发明。
图1为本实施例提供的交流牵引电机定子槽内建模示意图。参见图1，该模型具体为：定子铁芯7内端面设置一绝缘槽，绝缘槽内设有定子线圈组(作为定子绕组等效模型1)，定子绕组等效模型1设有一隔离层，隔离层包括定子绕组等效模型1表面的绝缘层2、气隙3，绝缘层2设于定子绕组等效模型1表面，绝缘槽的外端面有两槽壁4、绝缘槽槽口设置槽楔5，槽壁4和槽楔5与绝缘层2之间为气隙3。其中，定子绕组等效模型1表面的绝缘层2是指将定子绕组多匝铜线表面包裹的聚酰亚胺绝缘层、浸渍漆、漆膜等多种绝缘材料等效为定子绕组等效模型1表面的单一绝缘层2；气隙3是指将定子绕组中多匝铜线间及铜线与定子铁芯槽壁4、槽楔5间的气隙3等效为单匝铜线与定子铁芯槽壁4、槽楔5间的气隙3，并设气隙3的宽度处处相等。
图2为本实施例提供的交流牵引电机三维有限元模型示意图。图2为该交流牵引电机有限元模型的一部分，至少包括定子绕组等效模型1，定子绕组等效模型表面的绝缘层2，定子绕组等效模型表面绝缘层与定子铁芯槽壁及槽楔间的气隙3、定子铁芯7、定子楔部8、转子铁芯10、转子导条9、转子端环12、机壳6、轴承13以及定子和转子之间的气隙11；具体的，定子铁芯7固定装配在机壳6上，定子绕组等效模型1嵌于定子铁芯7槽内，在定子绕组等效模型1表面覆盖有绝缘层2，在绝缘层2与定子铁芯7槽壁及槽楔间存在气隙3，定子楔部8固定装配在定子铁芯7槽的端口上；转子铁芯10固定装配在轴承13上，转子导条9固定装配在转子铁芯10槽内，转子端环12固定装配在转子铁芯10轴向端面上。
基于上述定子槽内建模部位以及交流牵引电机三维有限元模型，下面给出一种用于交流牵引电机温度场分析的定子绕组气隙建模方法的可能实现方式。具体的，图3和图4为本发明实施例提供的用于交流牵引电机温度场分析的定子绕组气隙建模方法流程图及进一步详解图，该方法包括如下步骤(图中S1表示步骤1，S2表示步骤2，S3表示步骤3，S4表示步骤4，S4-1表示步骤4-1，S4-2表示步骤4-2，S4-3表示步骤4-3，S4-4表示步骤4-4，S5表示步骤5，S6表示步骤6，S7表示步骤7)：
步骤1)根据交流牵引电机定子绕组的结构建立定子绕组等效模型1；
在建模过程中将定子绕组中多匝铜线等效为单匝铜线，单匝铜线的截面积与多匝铜线的截面积相等，且单匝铜线的截面形状与定子槽型相同，单匝铜线的截面积通过公式(1)求得：
s1＝nπr1 2 (1)
式中：s1为单匝铜线的截面积，r1为定子绕组多匝铜线中每匝铜线的半径，n为定子绕组多匝铜线的匝数。
步骤2)根据步骤1所得定子绕组等效模型1确定该模型与定子铁芯槽壁4及槽楔5间隔离层的结构与宽度，其中所述隔离层为定子绕组等效模型1与定子铁芯槽壁4及槽楔5间的间隔，且隔离层的宽度处处相等，所述隔离层的宽度通过公式(2)计算而得：
Figure BDA0001941629540000061
式中：h为隔离层的宽度，r为定子槽底半径，s2为定子槽截面面积。
步骤3)根据步骤2所得隔离层的结构与宽度分别确定隔离层中绝缘层2与气隙3的结构与初始宽度，再根据电机其余部件的结构参数(即机壳6，定子铁芯7，定子楔部8，转子导条9，转子铁芯10，定子和转子之间的气隙11，转子端环12，轴承13)，建立交流牵引电机的三维有限元模型。
步骤3中隔离层中对应气隙3处在绝缘层2与定子铁芯槽壁4及槽楔5之间，其宽度处处相等，其初始宽度设定为定子槽内单层铜线的等效气隙宽度，具体为：
Figure BDA0001941629540000062
式中：h1为隔离层中气隙3的初始宽度。
又步骤3中隔离层中对应绝缘层2紧贴于定子绕组等效模型1的外表面，其宽度也是处处相等，其初始宽度为：
h2＝h-h1 (4)
式中：h2为隔离层中绝缘层2的初始宽度。
步骤4)根据步骤3所得交流牵引电机的三维有限元模型进行温度场分析，得到交流牵引电机的温度场分布图，并选取定子绕组等效模型表面某处的温度作为待测值，具体为：
步骤4-1)针对交流牵引电机三维有限元模型进行网格划分。
步骤4-2)针对网格划分后的交流牵引电机三维有限元模型施加热源；具体地，热源指交流牵引电机运转中产生的损耗，交流牵引电机的损耗包括定子绕组铜耗、转子导条铝耗、铁芯中的铁耗、机械损耗。
步骤4-3)设置交流牵引电机三维有限元模型的边界条件和对流换热系数，具体为：
(a)机壳、散热筋与空气间的对流换热系数设置为处处相等；
(b)根据公式(5)和公式(6)计算定子和转子之间气隙的雷诺数Re及其临界雷诺数Rel；
Figure BDA0001941629540000071
Figure BDA0001941629540000072
式中：d1为转子半径，d2为定子半径，δ为定子和转子之间气隙的长度，ωg为转子转速，v为空气的运动粘度；
(c)根据步骤b所得雷诺数Re及其临界雷诺数Rel的大小确定相应的对流换热系数，将定子和转子之间气隙中的对流换热系数设置为计算所得数值，具体为：
当Re<Rel时，表明此时气隙内空气流动为层流，取对流换热系数α为空气的导热系数，且α＝0.0267W/m·K；
当Re>Rel时，表明此时气隙内空气流动为湍流，则对流换热系数α通过公式(7)计算而得：
Figure BDA0001941629540000073
式中，λ为空气的导热系数。
其中，边界条件设置具体包括：交流牵引电机机壳的边界温度设置，交流牵引电机三维有限元模型两轴向对称侧面散热方式设置；交流牵引电机机壳的边界温度可设置为环境温度，交流牵引电机三维有限元模型两轴向对称侧面散热方式设置为绝热。
步骤4-4)针对三维有限元模型进行温度场有限元计算，从而获得交流牵引电机温度场分布图，具体的，图5为本发明实施例提供的交流牵引电机温度场分布示意图，并选取绕组等效模型表面待测点14的温度作为待测值。
步骤5)按一定间距改变定子绕组等效模型1的气隙3宽度，再分别建立交流牵引电机的三维有限元模型，求解不同气隙3宽度下交流牵引电机的温度场分布，并分别测取定子绕组等效模型表面1与步骤4中待测点14对应的温度；其中，按一定间距改变绕组等效模型1的气隙3宽度，是指以初始气隙3宽度为初始值，以隔离层宽度为最终值，根据数值拟合的需要确定间距的大小，从初始值开始按间距逐步增大气隙3宽度，直至最终值为止。
具体的，表1为本发明实施例提供的定子绕组等效模型1的8组气隙宽度与其相应的温度数据。如初始气隙3宽度为0.1mm，隔离层宽度为0.45mm，故绕组等效模型1的气隙3宽度范围为0.1～0.45mm，在该绕组等效模型1的气隙3宽度范围内按0.05mm间距取8个绕组等效模型1的气隙3宽度数据，并分别建立交流牵引电机的三维有限元模型，求解不同气隙3宽度下交流牵引电机的温度场分布，并分别在定子绕组等效模型1表面的待测点14测取相应的温度数据。
步骤6)针对步骤5所得N组定子绕组等效模型1的气隙宽度与其相应的温度数据进行数值拟合，数值拟合方法优选采用最小二乘法，得到定子绕组等效模型1气隙3宽度与温度间的计算公式(8)：
h1(T)＝aebT+cedT (8)
式中：h1(T)为绕组气隙宽度函数；T为绕组等效模型表面某处的温度；a、b、c、d分别为气隙宽度函数的系数，a、b、c、d根据最小二乘法进行确定，e为自然对数函数的底数。
具体地，根据表1中的数据，并根据最小二乘法利用MATLAB分析软件即可计算出系数a、b、c、d分别为：a＝0.1737；b＝0.007736；c＝-217.7；d＝-0.07486。
步骤7)实测电机定子绕组与其等效模型表面待测点14对应处的实际温度，并将该实测温度代入上述计算公式(8)，得到与该实测温度对应的气隙宽度，该气隙宽度即作为当前交流牵引电机定子绕组建模所对应的最佳气隙宽度。
表1
Figure BDA0001941629540000081
上述只是本发明的较佳实施例，并非对本发明作任何形式上的限制。虽然本发明已以较佳实施例揭露如上，然而并非用以限定本发明。因此，凡是未脱离本发明技术方案的内容，依据本发明技术实质对以上实施例所做的任何简单修改、等同变化及修饰，均应落在本发明技术方案保护的范围内。

Claims (9)

1.用于交流牵引电机温度场分析的定子绕组气隙建模方法，其特征在于，包括如下步骤：
步骤1)根据交流牵引电机定子绕组的结构建立定子绕组等效模型；
步骤2)根据步骤1所得定子绕组等效模型确定该模型与定子铁芯槽壁及槽楔间隔离层的结构与宽度；
步骤3)根据步骤2所得隔离层的结构与宽度确定隔离层中对应绝缘层与气隙的结构与初始宽度，再根据电机其余部件的结构参数，建立交流牵引电机的三维有限元模型；
步骤4)根据步骤3所得交流牵引电机的三维有限元模型进行温度场分析，得到交流牵引电机的温度场分布图，并选取定子绕组等效模型表面某处的温度作为待测值；
步骤5)按间距设定值改变定子绕组等效模型的气隙宽度，再分别建立交流牵引电机的三维有限元模型，求解不同气隙宽度下交流牵引电机的温度场分布，并分别测取定子绕组等效模型表面与步骤4中某处对应的温度；
步骤6)针对步骤5所得N组定子绕组等效模型的气隙宽度与其相应的温度数据进行数值拟合，得到定子绕组等效模型气隙宽度与温度间的计算公式(8)：
h1(T)＝aebT+cedT (8)
式中：h1(T)为绕组气隙宽度函数；T为绕组等效模型表面某处的温度；a、b、c、d分别为气隙宽度函数的系数；e为自然对数函数的底数；
步骤7)实测电机定子绕组与其等效模型表面某处对应的实际温度，并将实测的实际温度代入上述计算公式(8)，得到与实测的实际温度对应的气隙宽度，该气隙宽度即作为当前交流牵引电机定子绕组建模所对应的最佳气隙宽度。
2.根据权利要求1所述的用于交流牵引电机温度场分析的定子绕组气隙建模方法，其特征在于，步骤1在建模过程中将定子绕组中多匝铜线等效为单匝铜线，单匝铜线的截面积与多匝铜线的截面积相等，且单匝铜线的截面形状与定子槽型相同，单匝铜线的截面积通过公式(1)求得：
s1＝nπr1 2 (1)
式中：s1为单匝铜线的截面积，r1为定子绕组多匝铜线中每匝铜线的半径，n为定子绕组多匝铜线的匝数。
3.根据权利要求1所述的用于交流牵引电机温度场分析的定子绕组气隙建模方法，其特征在于，步骤2中定子绕组等效模型与定子铁芯槽壁及槽楔间隔离层的宽度处处相等，所述隔离层的宽度通过公式(2)计算而得：
Figure FDA0002495799750000021
式中：h为隔离层的宽度，r为定子槽底半径，s2为定子槽截面面积。
4.根据权利要求1所述的用于交流牵引电机温度场分析的定子绕组气隙建模方法，其特征在于，步骤3中隔离层中对应气隙的初始宽度设定为定子槽内单层铜线的等效气隙宽度，具体为：
Figure FDA0002495799750000022
式中：h1为隔离层中气隙的初始宽度。
5.根据权利要求1所述的用于交流牵引电机温度场分析的定子绕组气隙建模方法，其特征在于，步骤3中隔离层中对应绝缘层的初始宽度为：
h2＝h-h1 (4)
式中：h2为隔离层中绝缘层的初始宽度。
6.根据权利要求1所述的用于交流牵引电机温度场分析的定子绕组气隙建模方法，其特征在于，步骤4的具体操作是：
步骤4-1)针对交流牵引电机三维有限元模型进行网格划分；
步骤4-2)针对网格划分后的交流牵引电机三维有限元模型施加热源；
步骤4-3)设置交流牵引电机三维有限元模型的边界条件和对流换热系数；
步骤4-4)针对三维有限元模型进行温度场有限元计算，从而获得交流牵引电机温度场分布图，并选取绕组等效模型表面某处的温度作为待测值。
7.根据权利要求6所述的用于交流牵引电机温度场分析的定子绕组气隙建模方法，其特征在于，步骤4-3的具体操作为：
(a)机壳、散热筋与空气间的对流换热系数设置为处处相等；
(b)根据公式(5)和公式(6)计算定子和转子之间气隙的雷诺数Re及其临界雷诺数Rel；
Figure FDA0002495799750000023
Figure FDA0002495799750000024
式中：d1为转子半径，d2为定子半径，δ为定子和转子之间气隙的长度，ωg为转子转速，v为空气的运动粘度；
(c)根据步骤b所得雷诺数Re及其临界雷诺数Rel的大小确定相应的对流换热系数，将定子和转子之间气隙中的对流换热系数设置为计算所得数值。
8.根据权利要求7所述的用于交流牵引电机温度场分析的定子绕组气隙建模方法，其特征在于，步骤c中确定相应的对流换热系数的条件是：
当Re<Rel时，表明此时气隙内空气流动为层流，取对流换热系数α为空气的导热系数；
当Re>Rel时，表明此时气隙内空气流动为湍流，则对流换热系数α通过公式(7)计算而得：
Figure FDA0002495799750000031
式中：λ为空气的导热系数。
9.根据权利要求1所述的用于交流牵引电机温度场分析的定子绕组气隙建模方法，其特征在于，步骤5中的间距设定值，是指以初始气隙宽度为初始值，以隔离层宽度为最终值，根据数值拟合的需要确定间距的大小，从初始值开始按间距逐步增大气隙宽度，直至最终值为止。
CN201910023474.0A 2019-01-10 2019-01-10 用于交流牵引电机温度场分析的定子绕组气隙建模方法 Active CN109753737B (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN201910023474.0A CN109753737B (zh) 2019-01-10 2019-01-10 用于交流牵引电机温度场分析的定子绕组气隙建模方法
PCT/CN2019/122279 WO2020143361A1 (zh) 2019-01-10 2019-11-30 用于交流牵引电机温度场分析的定子绕组气隙建模方法
US17/037,873 US11196374B2 (en) 2019-01-10 2020-09-30 Modeling method of stator winding air gap for temperature field analysis of AC traction motor

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910023474.0A CN109753737B (zh) 2019-01-10 2019-01-10 用于交流牵引电机温度场分析的定子绕组气隙建模方法

Publications (2)

Publication Number Publication Date
CN109753737A CN109753737A (zh) 2019-05-14
CN109753737B true CN109753737B (zh) 2020-07-10

Family

ID=66405424

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910023474.0A Active CN109753737B (zh) 2019-01-10 2019-01-10 用于交流牵引电机温度场分析的定子绕组气隙建模方法

Country Status (3)

Country Link
US (1) US11196374B2 (zh)
CN (1) CN109753737B (zh)
WO (1) WO2020143361A1 (zh)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109753737B (zh) * 2019-01-10 2020-07-10 湖南科技大学 用于交流牵引电机温度场分析的定子绕组气隙建模方法
CN110896265B (zh) * 2019-08-20 2021-04-27 湖南科技大学 开关磁阻电机温度场分析中定子双绕组等效气隙建模方法
CN112564358B (zh) * 2020-11-10 2022-06-28 湖南科技大学 交流牵引电机转子齿部轴向通风槽的设计方法
CN112446150B (zh) * 2020-11-25 2024-06-04 国网江苏省电力有限公司检修分公司 一种双水内冷同步调相机温度场三维可视化方法和系统
CN112966388B (zh) * 2021-03-15 2022-04-15 武汉船用电力推进装置研究所（中国船舶重工集团公司第七一二研究所） 一种扁线电机温度场分析中定子绕组等效方法
CN112803630B (zh) * 2021-03-22 2021-12-31 湖南科技大学 交流牵引电机降噪优化设计方法
CN112886587B (zh) * 2021-03-29 2022-09-13 北京世纪百合科技有限公司 一种输电线路杆塔塔头空气间隙的校核和表示的方法
CN113569454B (zh) * 2021-08-04 2024-04-12 大连交通大学 一种柴油机配气机构的模拟分析方法
CN114547909B (zh) * 2022-03-03 2024-07-05 北京信息科技大学 一种湿式摩擦副温度场确定方法及系统
CN115422780A (zh) * 2022-10-12 2022-12-02 中车长春轨道客车股份有限公司 牵引系统的温度仿真方法和装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8543365B1 (en) * 2004-10-25 2013-09-24 Novatorque, Inc. Computer-readable medium, a method and an apparatus for designing and simulating electrodynamic machines implementing conical and cylindrical magnets
CN106096157A (zh) * 2016-06-15 2016-11-09 天津大学 一种轴向通风永磁同步电机温度场计算方法
CN107844647A (zh) * 2017-11-03 2018-03-27 中国矿业大学 一种单相圆筒型开关磁阻直线电机二维磁路建模方法
CN107885955A (zh) * 2017-11-29 2018-04-06 中昇创举（天津）科技有限公司 电机温度场的仿真方法和装置
CN108733887A (zh) * 2018-04-12 2018-11-02 湖南科技大学 一种基于变密度对称网格划分的开关磁阻电机热分析方法

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1357275A1 (en) * 2002-04-26 2003-10-29 Visteon Global Technologies, Inc. Modelling of the thermal behaviour of a switched reluctance motor driving a supercharger of an internal combustion engine
TW200635179A (en) * 2005-03-29 2006-10-01 Delta Electronics Inc Starting method for brushless direct current motor
US7570074B2 (en) * 2005-05-09 2009-08-04 Square D Company Electronic overload relay for mains-fed induction motors
WO2012006605A2 (en) * 2010-07-08 2012-01-12 Tesla Motors, Inc. Rotor temperature estimation
FR2975502B1 (fr) * 2011-05-16 2014-02-28 Electricite De France Procede et systeme de controle de l'instant de mise sous tension d'un dispositif comprenant un circuit magnetique
GB201116643D0 (en) * 2011-09-27 2011-11-09 Potter Clarkson Llp Rotating machines
DE102014011911A1 (de) * 2014-08-12 2016-02-18 Hans-Jürgen Remus Berechnungsverfahren zur Auslegung von Reluktanzsystemen und ein Computerprogramm
CN105868485A (zh) * 2016-04-08 2016-08-17 中国矿业大学 一种开关磁阻直线电机磁路建模方法
CN105808887B (zh) * 2016-04-08 2018-10-23 中国矿业大学 一种气隙不对称开关磁阻直线电机磁路建模方法
CN106446364B (zh) * 2016-09-08 2019-06-25 东南大学 一种温度场-热路直接耦合的电机热分析方法
CN107391884B (zh) * 2017-08-28 2020-06-26 天津大学 基于等效热网络模型的双余度永磁同步电机温升计算方法
CN109753737B (zh) * 2019-01-10 2020-07-10 湖南科技大学 用于交流牵引电机温度场分析的定子绕组气隙建模方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8543365B1 (en) * 2004-10-25 2013-09-24 Novatorque, Inc. Computer-readable medium, a method and an apparatus for designing and simulating electrodynamic machines implementing conical and cylindrical magnets
CN106096157A (zh) * 2016-06-15 2016-11-09 天津大学 一种轴向通风永磁同步电机温度场计算方法
CN107844647A (zh) * 2017-11-03 2018-03-27 中国矿业大学 一种单相圆筒型开关磁阻直线电机二维磁路建模方法
CN107885955A (zh) * 2017-11-29 2018-04-06 中昇创举（天津）科技有限公司 电机温度场的仿真方法和装置
CN108733887A (zh) * 2018-04-12 2018-11-02 湖南科技大学 一种基于变密度对称网格划分的开关磁阻电机热分析方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
"基于Ansoft的永磁起动电机的建模与性能优化";崔天翔;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20160115;论文第1-28页 *

Also Published As

Publication number Publication date
CN109753737A (zh) 2019-05-14
US11196374B2 (en) 2021-12-07
WO2020143361A1 (zh) 2020-07-16
US20210013820A1 (en) 2021-01-14

Similar Documents

Publication Publication Date Title
CN109753737B (zh) 用于交流牵引电机温度场分析的定子绕组气隙建模方法
CN106126861B (zh) 空冷汽轮发电机齿部和轭背部的散热系数的计算方法
Boglietti et al. Analysis of the endwinding cooling effects in TEFC induction motors
Sixel et al. Ceramic 3D printed direct winding heat exchangers for improving electric machine thermal management
Tong et al. Totally enclosed self-circulation axial ventilation system design and thermal analysis of a 1.65-MW direct-drive PMSM
CN109086514B (zh) 一种基于多物理场耦合的电抗器温度场计算方法
CN108388700A (zh) 一种动态确定高压电缆暂态热路中绝缘最佳分层数的方法
CN108896209B (zh) 一种油浸式变压器热点温度监测方法
Perez et al. A stationary thermal model for smooth air-gap rotating electric machines
CN108920774A (zh) 一种油浸式变压器内部温度监测方法
CN112100809B (zh) 一种基于多物理场耦合仿真三相同轴超导电缆设计的方法
Sixel et al. Ceramic 3-D printed direct winding heat exchangers for thermal management of concentrated winding electric machines
Torriano et al. Heat transfer coefficient distribution on the pole face of a hydrogenerator scale model
Langheck et al. Evaluation of an integral injection molded housing for high power density synchronous machines with concentrated single-tooth winding
CN103353926A (zh) 一种电机温度分布实时监测方法
CN109359339A (zh) 带有通风沟的核电汽轮发电机转子温度场分布计算方法
Wang et al. Convective heat transfer characteristics on end-winding of stator immersed oil-cooled electrical machines for aerospace applications
Wang et al. Investigation of lumped-parameter thermal model and thermal parameters test for IPMSM
CN108733887B (zh) 一种基于变密度对称网格划分的开关磁阻电机热分析方法
Liu et al. Model calibration of oil jet and oil spray cooling in electrical machines with hairpin windings
CN110896265B (zh) 开关磁阻电机温度场分析中定子双绕组等效气隙建模方法
CN109060176B (zh) 一种油浸式变压器顶层油温监测方法
Ren et al. Multi-Physical Field Analysis and Optimization Design of the High-Speed Motor of an Air Compressor for Hydrogen Oxygen Fuel Cells
Tindall et al. Loss-at-source thermal modelling in salient pole alternators using 3-dimensional finite difference techniques
CN112199869A (zh) 一种考虑电缆沟道通风特性的电缆载流量计算方法

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