CN105245041B - 一种求解大型电机定转子通风沟表面散热系数的方法 - Google Patents
一种求解大型电机定转子通风沟表面散热系数的方法 Download PDFInfo
- Publication number
- CN105245041B CN105245041B CN201510666861.8A CN201510666861A CN105245041B CN 105245041 B CN105245041 B CN 105245041B CN 201510666861 A CN201510666861 A CN 201510666861A CN 105245041 B CN105245041 B CN 105245041B
- Authority
- CN
- China
- Prior art keywords
- rotor
- heat transfer
- coefficient
- stator
- size machine
- 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
Abstract
本发明涉及一种求解大型电机定转子通风沟表面散热系数的方法，包括如下步骤：(1)建立大型电机的二维瞬态电磁场求解模型，计算出大型电机定转子各部件的损耗；(2)建立该大型电机定转子直线段部分的三维稳态温度场求解模型；(3)计算大型电机定转子径向通风沟表面散热系数，并将其作为通风沟初始表面散热系数；(4)在步骤(1)、步骤(3)的基础上，对步骤(2)中建立的大型电机定转子直线段部分的三维稳态温度场求解模型进行迭代计算，得到最终的定转子径向通风沟表面散热系数。本发明所述的方法，实现精确确定定转子通风沟表面散热系数，从而准确计算铁心温升。
Description
技术领域
本发明涉及电机技术领域，具体说是一种求解大型电机定转子通风沟表面散热系数的方法。所述定转子通风沟指定转子径向通风沟。
背景技术
在具有轴径向通风系统的大中型电机中，定转子径向通风沟风速和表面散热系数对电机的定转子的温升影响较大，而定转子径向通风沟(简称为通风沟)的气流和传热是在非常复杂的情况下进行的。
由于通风沟表面粗糙、入口较陡、弯曲截面非常短且有突变，再加上旋转的转子通风沟所产生的离心力和科里奥力，故从转子流入定子的空气气流是没有规律的，在通风沟内构成复杂的热交换边界条件。
正是由于通风沟内空气流动状态的复杂性，目前，大型电机的定转子通风沟表面散热系数的计算仍然由经验公式给出，如何准确确定定转子通风沟表面散热系数成为准确计算铁心温升的一大关键。
发明内容
针对现有技术中存在的缺陷，本发明的目的在于提供一种求解大型电机定转子通风沟表面散热系数的方法，实现精确确定定转子通风沟表面散热系数，从而准确计算铁心温升。
为达到以上目的，本发明采取的技术方案是：
一种求解大型电机定转子通风沟表面散热系数的方法，其特征在于，包括如下步骤：
(1)根据大型电机的实际结构、尺寸、额定数据等已知参数，建立该大型电机的二维瞬态电磁场求解模型，对整个二维瞬态电磁场求解模型进行网格划分，计算出大型电机定转子各部件的基本损耗；
(2)根据《电工专业指导性技术文件(汽轮发电机电磁计算公式)》的解析方法与瞬态电磁场相结合的方法计算出大型电机定转子各部件的谐波损耗；
(3)根据大型电机定转子的直线段部分的实际尺寸，建立该大型电机定转子直线段部分的三维稳态温度场求解模型；
(4)根据流体相似理论及牛顿放热定律，采用解析方法计算出大型电机定转子径向通风沟表面散热系数，并将该表面散热系数作为大型电机定转子直线段部分的三维稳态温度场求解模型的初始值，即通风沟初始表面散热系数；
(5)在步骤(1)、(2)计算出大型电机定转子各部件的基本损耗和谐波损耗、步骤(4)计算出通风沟初始表面散热系数的基础上，对步骤(3)中建立的大型电机定转子直线段部分的三维稳态温度场求解模型进行迭代计算，得到最终的定转子径向通风沟表面散热系数。
在上述技术方案的基础上，步骤(1)中，二维瞬态电磁场求解模型包括：定子铁心1，定子线圈2，定子槽楔3，气隙4，转子槽楔5，励磁绕组6，副槽7，转子铁心8，阻尼绕组9以及包含计算模型的空气域。
在上述技术方案的基础上，步骤(1)中，二维瞬态电磁场求解模型如下：
式中：σ为电导率(S/m)；μ为磁导率(H/m)；t为时间(s)；Jz为电密的轴向分量；Az为矢量磁位(Wb/m)，只有z轴分量。
在上述技术方案的基础上，步骤(1)中，计算出大型电机定转子各部件的基本损耗具体包括至少以下三项：
三相定子绕组中基本铜损，
转子励磁线圈损耗，
电机内定转子铁心的基本铁损。
在上述技术方案的基础上，步骤(2)中，计算出大型电机定转子各部件的谐波损耗具体包括至少以下五项：
转子磁场高次谐波在定子表面产生的损耗，
转子齿谐波在定子表面产生的损耗，
转子齿谐波在定子齿中产生的脉动损耗，
定子磁场高次谐波在转子表面产生的损耗，
定子齿谐波在转子表面产生的损耗。
在上述技术方案的基础上，步骤(3)中，大型电机定转子直线段部分的三维稳态温度场求解模型分为定子部分求解模型和转子部分求解模型，包括：定子铁心1，定子线圈2，定子槽楔3，气隙4，转子槽楔5，励磁绕组6，副槽7，转子铁心8，阻尼绕组9，定子径向通风沟10，转子径向通风沟11。
在上述技术方案的基础上，步骤(3)中，大型电机定转子直线段部分的三维稳态温度场求解模型如下：
式中，T为物体的温度；λx、λy、λz分别为x、y、z方向的导热系数；q为热源的发热密度；α为S2表面散热系数；S1、S2分别为发热体的边界面，Tf为环境温度，T1为S1面上给定温度。
在上述技术方案的基础上，步骤(4)的具体步骤为：
步骤4.1，根据流体相似理论，得到流体呈紊流状态下径向通风沟内流体满足的相似准则方程：
式中，Re为流体的雷诺系数；υ为径向通风沟内流体的流速(m/s)；v为流体的粘度系数(s/m2)；d为径向通风沟等效直径(m)；Nu为努谢尔数；λ为流体的导热系数(W/(m·K))，α为表面散热系数(W/(m2·K))；
步骤4.2，当流体在径向通风沟内呈紊流运动时，努谢尔数相似准则方程与其他相似准则方程之间具有如下形式：
式中，Pr(CP)为通风沟内流体平均温度所对应的普朗特数；Pr(CT)为通风沟内壁面温度所对应的普朗特数；ε为修正系数；
步骤4.3，将式(3)、(4)联立，得到径向通风沟内流体在紊流状态下的初始表面散热系数：
在上述技术方案的基础上，步骤(5)的具体步骤为：
步骤5.1，根据步骤(4)得到的通风沟初始表面散热系数，获取定转子铁心温度；
步骤5.2，确定大型电机定转子径向通风沟表面散热系数及定转子铁心温度的收敛条件；
确定大型电机定转子径向通风沟表面散热系数的具体步骤如下：
由通风沟初始表面散热系数与定转子铁心温度之间的非线性关系，根据牛顿放热定律，当物体表面与周围存在温度差时，单位时间从单位面积散失的热量与温度差成正比，得到计算定转子铁心温度计算公式：
由于散热系数与电机铁心温度之间的非线性关系，因此对通风沟表面的散热系数，可得：
式中，T1为定转子靠近边界面外侧表面温度值；T2为定转子内侧的温度值；Tf为环境温度值，单位为(K)；X1为定子靠近边界面外侧表面节点横坐标值；X2为定子内表面节点横坐标值；α′i为第i个节点的散热系数近似值；α′i+1为第i-1个节点的散热系数近似值；α′i+1为第i+1个节点的散热系数近似值；αi为第i个节点的散热系数修正值；β为松弛因子；
确定定转子铁心温度的收敛条件的具体步骤如下：
所述通风沟表面散热系数及定转子铁心温度的收敛条件为：
式中：Ts为铁心的实测温度值；Tw为对应此处铁心的计算温度值；αi为第i个节点的散热系数；αi+1为第i+1个节点的散热系数；ε1和ε2分别为满足迭代要求的残差；
当迭代结果同时满足收敛条件(8)和(9)时，停止迭代，此时，根据式(7)计算的到的通风沟表面散热系数即为大型电机额定运行时通风沟表面散热系数。
本发明所述的求解大型电机定转子通风沟表面散热系数的方法，与传统方法相比具有以下优点：
1、给出了新方法，通过迭代的方式，将有限元温度场计算与解析公式相结合；
2、该新方法对通风沟散热系数进行计算，可准确地计算大型电机定转子通风沟表面散热系数；
3、进一步提高了对大型电机铁心温升计算的精度；
4、该新方法具有计算准确且计算时间短的优点。
本发明所述的求解大型电机定转子通风沟表面散热系数的方法，以大型汽轮发电机定转子径向通风沟散热系数计算为例，但该方法不仅仅局限于大型汽轮发电机定转子径向通风沟散热系数的计算，同时该方法同样适合于任意种类带径向通风结构的大型电机。但是在计算损耗的过程中，不同电机需要考虑的损耗应该以该电机直线段部分所需计算风沟所在部位的损耗为准，例如计算转子附槽通风结构散热系数时需要计算转子部分的附加损耗等，计算永磁电机转子通风沟散热系数时需要计算永磁体上的涡流损耗等。
附图说明
本发明有如下附图：
图1本发明的流程图。
图2大型电机的二维瞬态电磁场求解模型。
图3大型电机的三维稳态温度场定子部分求解模型。
图4大型电机的三维稳态温度场转子部分求解模型。
具体实施方式
以下结合附图对本发明作进一步详细说明。
如图1所示，本发明所述的求解(计算)大型电机定转子通风沟表面散热系数的方法，包括如下步骤：
(1)根据大型电机的实际结构、尺寸、额定数据等已知参数，建立该大型电机的二维瞬态电磁场求解模型，对整个二维瞬态电磁场求解模型进行网格划分，计算出大型电机定转子各部件的基本损耗；
求解模型亦可称为数值计算模型；
(2)根据《电工专业指导性技术文件(汽轮发电机电磁计算公式)》的解析方法与瞬态电磁场相结合的方法计算出大型电机定转子各部件的谐波损耗；
(3)根据大型电机定转子的直线段部分的实际尺寸，建立该大型电机定转子直线段部分的三维稳态温度场求解模型；
(4)根据流体相似理论及牛顿放热定律，采用解析方法计算出大型电机定转子径向通风沟表面散热系数，并将该表面散热系数作为大型电机定转子直线段部分的三维稳态温度场求解模型的初始值，即通风沟初始表面散热系数；
(5)在步骤(1)、(2)计算出大型电机定转子各部件的基本损耗和谐波损耗、步骤(4)计算出通风沟初始表面散热系数的基础上，对步骤(3)中建立的大型电机定转子直线段部分的三维稳态温度场求解模型进行迭代计算，得到最终的定转子径向通风沟表面散热系数。
在上述技术方案的基础上，如图2所示，步骤(1)中，二维瞬态电磁场求解模型包括：定子铁心1，定子线圈2，定子槽楔3，气隙4，转子槽楔5，励磁绕组6，副槽7，转子铁心8，阻尼绕组9以及包含计算模型的空气域。
在上述技术方案的基础上，步骤(1)中，二维瞬态电磁场求解模型如下：
式中：σ为电导率(S/m)；μ为磁导率(H/m)；t为时间(s)；Jz为电密的轴向分量；Az为矢量磁位(Wb/m)，只有z轴分量。
在上述技术方案的基础上，步骤(1)中，计算出大型电机定转子各部件的基本损耗具体包括至少以下三项：
三相定子绕组中基本铜损，
转子励磁线圈损耗，
电机内定转子铁心的基本铁损。
在上述技术方案的基础上，步骤(2)中，计算出大型电机定转子各部件的谐波损耗具体包括至少以下五项：
转子磁场高次谐波在定子表面产生的损耗，
转子齿谐波在定子表面产生的损耗，
转子齿谐波在定子齿中产生的脉动损耗，
定子磁场高次谐波在转子表面产生的损耗，
定子齿谐波在转子表面产生的损耗。
在上述技术方案的基础上，如图3、4所示，步骤(3)中，大型电机定转子直线段部分的三维稳态温度场求解模型分为定子部分求解模型和转子部分求解模型，包括：定子铁心1，定子线圈2，定子槽楔3，气隙4，转子槽楔5，励磁绕组6，副槽7，转子铁心8，阻尼绕组9，定子径向通风沟10，转子径向通风沟11。
在上述技术方案的基础上，步骤(3)中，大型电机定转子直线段部分的三维稳态温度场求解模型如下：
式中，T为物体的温度；λx、λy、λz分别为x、y、z方向的导热系数；q为热源的发热密度；α为S2表面散热系数；S1、S2分别为发热体的边界面，Tf为环境温度，T1为S1面上给定温度。
在上述技术方案的基础上，步骤(4)的具体步骤为：
步骤4.1，根据流体相似理论，得到流体呈紊流状态下径向通风沟内流体满足的相似准则方程：
式中，Re为流体的雷诺系数；υ为径向通风沟内流体的流速(m/s)；v为流体的粘度系数(s/m2)；d为径向通风沟等效直径(m)；Nu为努谢尔数；λ为流体的导热系数(W/(m·K))，α为表面散热系数(W/(m2·K))；
步骤4.2，当流体在径向通风沟内呈紊流运动时，努谢尔数相似准则方程与其他相似准则方程之间具有如下形式：
式中，Pr(CP)为通风沟内流体平均温度所对应的普朗特数；Pr(CT)为通风沟内壁面温度所对应的普朗特数；ε为修正系数；
步骤4.3，将式(3)、(4)联立，得到径向通风沟内流体在紊流状态下的初始表面散热系数：
在上述技术方案的基础上，步骤(5)的具体步骤为：
步骤5.1，根据步骤(4)得到的通风沟初始表面散热系数，获取定转子铁心温度；
步骤5.2，确定大型电机定转子径向通风沟表面散热系数及定转子铁心温度的收敛条件；
确定大型电机定转子径向通风沟表面散热系数的具体步骤如下：
由通风沟初始表面散热系数与定转子铁心温度之间的非线性关系，根据牛顿放热定律，当物体表面与周围存在温度差时，单位时间从单位面积散失的热量与温度差成正比，得到计算定转子铁心温度计算公式：
由于散热系数与电机铁心温度之间的非线性关系，因此对通风沟表面的散热系数，可得：
式中，T1为定转子靠近边界面外侧表面温度值；T2为定转子内侧的温度值；Tf为环境温度值，单位为(K)；X1为定子靠近边界面外侧表面节点横坐标值；X2为定子内表面节点横坐标值；α′i为第i个节点的散热系数近似值；α′i+1为第i-1个节点的散热系数近似值；α′i+1为第i+1个节点的散热系数近似值；αi为第i个节点的散热系数修正值；β为松弛因子；
确定定转子铁心温度的收敛条件的具体步骤如下：
所述通风沟表面散热系数及定转子铁心温度的收敛条件为：
式中：Ts为铁心的实测温度值；Tw为对应此处铁心的计算温度值；αi为第i个节点的散热系数；αi+1为第i+1个节点的散热系数；ε1和ε2分别为满足迭代要求的残差；
当迭代结果同时满足收敛条件(8)和(9)时，停止迭代，此时，根据式(7)计算的到的通风沟表面散热系数即为大型电机额定运行时通风沟表面散热系数。
本发明所述的求解方法是为了解决现有大型电机在进行温度场计算时，需要准确确定大型电机定转子表面散热系数大小的问题。本发明基于流体相似理论和牛顿放热定律，根据解析方法与二维瞬态电磁场与三维稳态温度场相结合的数值计算方法，在计算出大型电机定转子通风沟初始散热系数情况下，通过一系列数学迭代方法，求解出大型电机温度场分布情况，并与实测温度值进行比较，最终得到满足工程要求精度的散热系数。该方法有收敛速度快，计算量小的显著优势。为确定定转子通风沟散热系数分布的求解方法，提供了新的科学思路。
为了更清楚地说明本发明，以大型汽轮发电机定子径向通风沟表面散热系数计算为例，对本发明作进一步详细说明。
步骤1，根据大型电机的实际结构和尺寸，建立了大型电机的二维瞬态电磁场求解模型(数值计算模型)，二维瞬态电磁场数值计算模型包括：定子铁心1，定子线圈2，定子槽楔3，气隙4，转子槽楔5，励磁绕组6，副槽7，转子铁心8，阻尼绕组9以及包含计算模型的空气域。
根据二维瞬态电磁场中各个部件的实际尺寸及对计算量和计算时间的考虑，确定了各个部件关键剖分区域以及最小划分网格尺寸后，对整个二维瞬态电磁场计算模型进行了网格划分。
通过对瞬态电磁场数学方程(1-1)-(1-9)的求解，确定了定子绕组铜损、铁心部分的涡流损耗、大型电机存在的定子表面损耗及各部分的附加损耗等各部分的损耗值；
式中：σ为电导率(S/m)；μ为磁导率(H/m)；t为时间(s)；Jz为轴向上的电密分量；Az为矢量磁位(Wb/m)，只有z轴分量。
1三相定子绕组中基本铜损采用如下公式计算：
其中，为定子相电流，R(75)为75°时定子每相电阻。
2额定负载情况下转子励磁线圈损耗采用如下公式计算：
当采用有刷励磁时：
当采用无刷励磁时：
其中，IfN为负载励磁电流，R2(75)为75°时每对极绕组电阻，ηf为励磁机的效率。
3转子磁场高次谐波在定子表面产生的损耗采用如下公式计算：
式中，为转子磁场高次谐波在定子表面产生损耗的计算系数，AWo为空载额定电压时每极的总磁动势，Kδ为总的气隙系数，δ是气隙实际长度，Di为定子铁心内径，p为电机极对数，ρ为电阻率，Z1为定子槽数。
4转子齿谐波在定子表面产生的损耗采用如下公式计算：
其中，为转子齿谐波在定子表面产生损耗的计算系数，GFez为定子齿部质量，Z′2为转子槽分度数。
5转子齿谐波在定子齿中产生的脉动损耗采用如下公式计算：
其中，为转子齿谐波在定子齿中产生的脉动损耗的计算系数。
6定子磁场高次谐波在转子表面产生的损耗采用如下公式计算：
其中，l2为转子本体长，S1为铁心外圆面积，
7定子齿谐波在转子表面产生的损耗采用如下公式计算：
其中，为定子齿谐波在转子表面产生的损耗的计算系数，t1为定子齿距，
8电机内部产生的基本铁损采用如下公式计算：
pFe＝ph+pc＝khfB2+kc(fB)2 (1-9)
pFe为单位铁心损耗；ph为单位磁滞损耗；pc为单位涡流损耗；kh为磁滞损耗系数；kc为涡流损耗系数；f为电机运行频率；B为正弦磁密幅值。
步骤2，对大型电机定转子直线段部分建立三维稳态温度场数值计算模型，该计算模型以定子计算模型部分为例，定子计算模型部分包含了定子铁心1，定子线圈2，定子槽楔3，气隙4，定子径向通风沟10。根据定子上轴径向通风沟的实际尺寸、位置及数量，确定了各个通风沟所有边上的网格节点数后对三维电机直线段定子表面散热系数的计算模型进行网格剖分，；
所述步骤2中计算电机定转子直线段部分三维稳态温度场数学模型如下：
式中，T为求解域内温度；λx、λy、λz分别为x、y、z方向的导热系数；q为热源的发热密度；α为S2表面散热系数；S1、S2分别为发热体的边界面，Tf为环境温度，T1为S1面上给定温度。
步骤3，根据流体相似理论及牛顿放热定律计算得到通风沟初始表面散热系数。
根据流体相似理论得到流体呈紊流状态下定子径向通风沟内流体满足的相似准则方程：
式中，Re为流体的雷诺系数；υ为径向通风沟内流体的流速(m/s)；v为流体的粘度系数(s/m2)；d为径向通风沟等效直径(m)；Nu为努谢尔数；λ为流体的导热系数(W/(m·K))，α为表面散热系数(W/(m2·K))。
当流体在径向通风沟内呈紊流运动时，努谢尔数相似准则方程与其他相似准则方程之间具有如下形式：
式中Pr(CP)为通风沟内流体平均温度所对应的普朗特数；Pr(CT)为通风沟内壁面温度所对应的普朗特数；ε为修正系数。
考虑到在径向通风沟内流体的物性参数受温度变化的影响很小，将式(3)、(4)联立，得到通风沟内流体在紊流状态下的初始表面散热系数：
步骤4，根据步骤1中确定的大型发电机定转子各个部分的损耗分、步骤3中确定的大型发电机定转子径向通风沟初始表面散热系数的基础上，对步骤2中建立的大型电机定转子直线段部分三维稳态温度场模型进行迭代计算，得到最终的定转子径向通风沟表面散热系数。
确定定子铁心温度的收敛条件，获取定子铁心温度，得到计算定子铁心温度计算公式
根据散热系数与电机铁心温度之间的非线性关系，可得：
式中，T1为定转子靠近边界面外侧表面温度值；T2为定转子内侧的温度值；Tf为环境温度值，单位为(K)；X1为定子靠近边界面外侧表面节点横坐标值；X2为定子内表面节点横坐标值；α′i为第i个节点的散热系数近似值；α′i-1为第i-1个节点的散热系数近似值；α′i+1为第i+1个节点的散热系数近似值；αi为第i个节点的散热系数修正值；β为松弛因子。
为了能够更快的收敛和求得更为准确的定子通风沟表面散热系数，通风沟表面网格节点的环境温度按照整个叠片段区域流体流体温度的近似线型分布进行赋值，然后求解非线性方程组。常用的求解方法有牛顿拉佛逊迭代法和欠松弛迭代法。采用欠松弛方法迭代，β取为0.8求解非线性方程组。当收敛条件同时满足式(8)-(9)中两个条件时停止迭代。
式中：Ts为铁心的实测温度值；Tw为对应此处铁心的计算温度值；αi为第i个节点的散热系数；αi+1为第i+1个节点的散热系数；ε1和ε2分别为满足迭代要求的残差；
当迭代结果同时满足收敛条件(8)和(9)时，停止迭代，此时，根据式(7)计算的到的通风沟表面散热系数即为大型电机额定运行时通风沟表面散热系数。
迭代过程中定转子的铁心温度可按照如下公式进行计算：
式中，T为物体的温度；λx、λy、λz分别为x、y、z方向的导热系数；q为热源的发热密度。
若发热体的边界面S由S1和S2两部分组成，则S1和S2上的边界条件分别为：
式中T1为S1面上的给定温度；Tf为环境温度；α为S2面上的散热系数。
显然，本发明的上述实施例仅仅是为了清除说明本发明所做的一种可能性举例，而并非是对本发明的实施方式的限定，对于所属领域的普通技术人员来说，在上述说明的基础上，本发明还可以用于电机其他带有通风结构的部位的通风沟散热系数的计算以及不同类型带有通风结构的大中型电机的通风沟散热系数的计算。这里无法对所有的实施方式予以穷举，凡是属于本发明的技术方案所引申出的显而易见的变化或变动仍处于本发明的保护范围之内。
本说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。
Claims (7)
1.一种求解大型电机定转子通风沟表面散热系数的方法，其特征在于，包括如下步骤：
(1)根据大型电机的实际结构、尺寸、额定数据及其他已知参数，建立该大型电机的二维瞬态电磁场求解模型，对整个二维瞬态电磁场求解模型进行网格划分，计算出大型电机定转子各部件的基本损耗；
(2)根据《电工专业指导性技术文件(汽轮发电机电磁计算公式)》的解析方法与瞬态电磁场相结合的方法计算出大型电机定转子各部件的谐波损耗；
(3)根据大型电机定转子的直线段部分的实际尺寸，建立该大型电机定转子直线段部分的三维稳态温度场求解模型；
(4)根据流体相似理论及牛顿放热定律，采用解析方法计算出大型电机定转子径向通风沟表面散热系数，并将该表面散热系数作为大型电机定转子直线段部分的三维稳态温度场求解模型的初始值，即通风沟初始表面散热系数；
(5)在步骤(1)、(2)计算出大型电机定转子各部件的基本损耗和谐波损耗、步骤(4)计算出通风沟初始表面散热系数的基础上，对步骤(3)中建立的大型电机定转子直线段部分的三维稳态温度场求解模型进行迭代计算，得到最终的定转子径向通风沟表面散热系数；
步骤(4)的具体步骤为：
步骤4.1，根据流体相似理论，得到流体呈紊流状态下径向通风沟内流体满足的相似准则方程：
式中，Re为流体的雷诺系数；υ为径向通风沟内流体的流速(m/s)；v为流体的粘度系数(s/m2)；d为径向通风沟等效直径(m)；Nu为努谢尔数；λ为流体的导热系数(W/(m·K))，α为表面散热系数(W/(m2·K))；
步骤4.2，当流体在径向通风沟内呈紊流运动时，努谢尔数相似准则方程与其他相似准则方程之间具有如下形式：
式中，Pr(CP)为通风沟内流体平均温度所对应的普朗特数；Pr(CT)为通风沟内壁面温度所对应的普朗特数；ε为修正系数；
步骤4.3，将式(3)、(4)联立，得到径向通风沟内流体在紊流状态下的初始表面散热系数：
步骤(5)的具体步骤为：
步骤5.1，根据步骤(4)得到的通风沟初始表面散热系数，获取定转子铁心温度；
步骤5.2，确定大型电机定转子径向通风沟表面散热系数及定转子铁心温度的收敛条件；
确定大型电机定转子径向通风沟表面散热系数的具体步骤如下：
由通风沟初始表面散热系数与定转子铁心温度之间的非线性关系，根据牛顿放热定律，当物体表面与周围存在温度差时，单位时间从单位面积散失的热量与温度差成正比，得到计算定转子铁心温度计算公式：
由于散热系数与电机铁心温度之间的非线性关系，因此对通风沟表面的散热系数，可得：
式中，T1为定转子靠近边界面外侧表面温度值；T2为定转子内侧的温度值；Tf为环境温度值，单位为(K)；X1为定子靠近边界面外侧表面节点横坐标值；X2为定子内表面节点横坐标值；α′i为第i个节点的散热系数近似值；α′i-1为第i-1个节点的散热系数近似值；α′i+1为第i+1个节点的散热系数近似值；αi为第i个节点的散热系数修正值；β为松弛因子；
确定定转子铁心温度的收敛条件的具体步骤如下：
所述通风沟表面散热系数及定转子铁心温度的收敛条件为：
式中：Ts为铁心的实测温度值；Tw为对应此处铁心的计算温度值；αi为第i个节点的散热系数；αi+1为第i+1个节点的散热系数；ε1和ε2分别为满足迭代要求的残差；
当迭代结果同时满足收敛条件(8)和(9)时，停止迭代，此时，根据式(7)计算的到的通风沟表面散热系数即为大型电机额定运行时通风沟表面散热系数。
2.如权利要求1所述的求解大型电机定转子通风沟表面散热系数的方法，其特征在于：步骤(1)中，二维瞬态电磁场求解模型包括：定子铁心(1)，定子线圈(2)，定子槽楔(3)，气隙(4)，转子槽楔(5)，励磁绕组(6)，副槽(7)，转子铁心(8)，阻尼绕组(9)以及包含计算模型的空气域。
3.如权利要求1所述的求解大型电机定转子通风沟表面散热系数的方法，其特征在于：步骤(1)中，二维瞬态电磁场求解模型如下：
式中：σ为电导率(S/m)；μ为磁导率(H/m)；t为时间(s)；Jz为电密的轴向分量；Az为矢量磁位(Wb/m)，只有z轴分量。
4.如权利要求1所述的求解大型电机定转子通风沟表面散热系数的方法，其特征在于：步骤(1)中，计算出大型电机定转子各部件的基本损耗具体包括至少以下三项：
三相定子绕组中基本铜损，
转子励磁线圈损耗，
电机内定转子铁心的基本铁损。
5.如权利要求1所述的求解大型电机定转子通风沟表面散热系数的方法，其特征在于：步骤(2)中，计算出大型电机定转子各部件的谐波损耗具体包括至少以下五项：
转子磁场高次谐波在定子表面产生的损耗，
转子齿谐波在定子表面产生的损耗，
转子齿谐波在定子齿中产生的脉动损耗，
定子磁场高次谐波在转子表面产生的损耗，
定子齿谐波在转子表面产生的损耗。
6.如权利要求1所述的求解大型电机定转子通风沟表面散热系数的方法，其特征在于：步骤(3)中，大型电机定转子直线段部分的三维稳态温度场求解模型分为定子部分求解模型和转子部分求解模型，包括：定子铁心(1)，定子线圈(2)，定子槽楔(3)，气隙(4)，转子槽楔(5)，励磁绕组(6)，副槽(7)，转子铁心(8)，阻尼绕组(9)，定子径向通风沟(10)，转子径向通风沟(11)。
7.如权利要求1所述的求解大型电机定转子通风沟表面散热系数的方法，其特征在于：步骤(3)中，大型电机定转子直线段部分的三维稳态温度场求解模型如下：
式中，T为物体的温度；λx、λy、λz分别为x、y、z方向的导热系数；q为热源的发热密度；α为S2表面散热系数；S1、S2分别为发热体的边界面，Tf为环境温度，T1为S1面上给定温度。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510666861.8A CN105245041B (zh) | 2015-10-15 | 2015-10-15 | 一种求解大型电机定转子通风沟表面散热系数的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510666861.8A CN105245041B (zh) | 2015-10-15 | 2015-10-15 | 一种求解大型电机定转子通风沟表面散热系数的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105245041A CN105245041A (zh) | 2016-01-13 |
| CN105245041B true CN105245041B (zh) | 2018-11-06 |
Family
ID=55042572
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510666861.8A Active CN105245041B (zh) | 2015-10-15 | 2015-10-15 | 一种求解大型电机定转子通风沟表面散热系数的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105245041B (zh) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106777459A (zh) * | 2016-11-10 | 2017-05-31 | 北京交通大学 | 发电机轴径向转子温度场的计算方法 |
| CN107301302B (zh) * | 2017-07-12 | 2019-10-25 | 北京交通大学 | 一种汽轮发电机定子绕组附加损耗的计算方法 |
| CN109101693B (zh) * | 2018-07-16 | 2020-10-27 | 北京交通大学 | 基于风路-流体-传热协同耦合计算转子多物理场的方法 |
| CN109359339A (zh) * | 2018-09-17 | 2019-02-19 | 北京交通大学 | 带有通风沟的核电汽轮发电机转子温度场分布计算方法 |
| CN109738709B (zh) * | 2018-11-26 | 2021-01-05 | 上海电气电站设备有限公司 | 一种大型汽轮发电机端部电磁场、温度场计算方法 |
| CN110031506A (zh) * | 2019-04-25 | 2019-07-19 | 北京交通大学 | 永磁电机气隙导热系数的计算方法 |
| CN112560367B (zh) * | 2020-11-27 | 2022-05-27 | 上海电气电站设备有限公司 | 一种发电机轴向双孔铜排转子线圈热点温度计算方法 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104008228A (zh) * | 2014-04-28 | 2014-08-27 | 北京交通大学 | 一种求解发电机定子端部绕组表面散热系数的新方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103246764B (zh) * | 2013-04-19 | 2016-02-03 | 哈尔滨工业大学 | 基于Ansys Workbench的潜油电机温度场模拟方法 |
| CN104462624A (zh) * | 2013-09-16 | 2015-03-25 | 同济大学 | 一种基于多场耦合的电机温度场数据处理方法 |
-
2015
- 2015-10-15 CN CN201510666861.8A patent/CN105245041B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104008228A (zh) * | 2014-04-28 | 2014-08-27 | 北京交通大学 | 一种求解发电机定子端部绕组表面散热系数的新方法 |
Non-Patent Citations (1)
| Title |
|---|
| 异步牵引电机温度场分析;夏正泽;《中国优秀硕士学位论文工程科技Ⅱ辑》;20080831;正文第4-8、21、35-41、53-63页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105245041A (zh) | 2016-01-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105245041B (zh) | 一种求解大型电机定转子通风沟表面散热系数的方法 | |
| Dong et al. | Thermal optimization of a high-speed permanent magnet motor | |
| Li et al. | Thermal analysis of the permanent-magnet spherical motor | |
| Wang et al. | Ultrafast steady-state multiphysics model for PM and synchronous reluctance machines | |
| Fatemi et al. | Computationally efficient strand eddy current loss calculation in electric machines | |
| Acquaviva et al. | Computationally efficient modeling of electrical machines with cooling jacket | |
| Li et al. | Calculating the electromagnetic field and losses in the end region of a large synchronous generator under different operating conditions with 3-D transient finite-element analysis | |
| Nerg et al. | Lumped-parameter-based thermal analysis of a doubly radial forced-air-cooled direct-driven permanent magnet wind generator | |
| CN109101693A (zh) | 基于风路-流体-传热协同耦合计算转子多物理场的方法 | |
| Lu et al. | Simulation and analysis of thermal fields of rotor multislots for nonsalient-pole motor | |
| CN109359339A (zh) | 带有通风沟的核电汽轮发电机转子温度场分布计算方法 | |
| Xia et al. | Analyzing temperature rise and fluid flow of high-power-density and high-voltage induction motor in the starting process | |
| Liang et al. | Influence of critical parameters in lumped-parameter thermal models for electrical machines | |
| Bersch et al. | Thermal and electromagnetic stator vent design optimisation for synchronous generators | |
| Zhang et al. | Electromagnetic–Thermal Coupling Analysis of Permanent Magnet Synchronous Machines for Electric Vehicle Applications Based on Improved (${\mu+ 1} $) Evolution Strategy | |
| Wang et al. | Thermal analysis of a water-cooled interior permanent magnet traction machine | |
| Kowal et al. | The effect of the electrical steel properties on the temperature distribution in direct-drive PM synchronous generators for 5 MW wind turbines | |
| Zhiyong et al. | Thermal analysis of PMSM based on lumped parameter thermal network method | |
| Darabi et al. | Design of the forced water cooling system for a claw pole transverse flux permanent magnet synchronous motor | |
| Wang et al. | Rotor loss calculation and thermal analysis of a dual-stator axial-flux permanent magnet machine with combined rectangle-shaped magnets | |
| Kang et al. | Thermal analysis of the cooling system with the circulation between rotor holes of enclosed PMSMs based on modified models | |
| Min et al. | Coupling calculation of 3D whole domain steady flow and temperature field for underwater oil-filled brushless DC motors | |
| Naskar et al. | Numerical analysis of three dimensional steady state heat conduction in the rotor of an induction motor by finite element method | |
| Luu et al. | Effect of axial-layered permanent-magnet on operating temperature in outer rotor machine | |
| Akawung et al. | Thermal analysis of air cooling system for electric machines using lumped parameter and flow resistance network |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |