法律状态公告日
法律状态信息
法律状态
2022-05-27
未缴年费专利权终止 IPC(主分类):G06F17/50 专利号:ZL201610439958X 申请日:20160615 授权公告日:20190322
专利权的终止
2019-03-22
授权
授权
2016-12-07
实质审查的生效 IPC(主分类):G06F17/50 申请日:20160615
实质审查的生效
2016-11-09
公开
公开
技术领域
本发明属于电机温度场计算特性分析领域,尤其涉及一种基于将冷却介质与电机自身结构的热网络结合为一体模型的电机温度场计算方法。
背景技术
永磁同步电机与传统电励磁电机相比,具有结构简单、可靠性强、功率密度高等诸多优点,因而应用广泛。然而温度过高会导致永磁材料退磁,影响电机的效率、使用寿命和可靠性,所以温度场计算是永磁同步电机设计过程中必不可少的环节。在计算电机温度场时,常用有限元法和等效热网络法。有限元法计算结果准确,但是计算量大,耗时长,较大的计算负担不利于在电机优化设计大规模使用;等效热网络法具有计算量小,耗时短的特点,虽然计算精度没有有限元法高,但是计算精度在可以接受的范围之内,在电机温度场计算当中仍然有重要的作用。
最初的热网络法是将电机温度场离散化为网格,将分布参数转换为集中参数,构成电机等效热网络,可根据电机的具体结构以及期望得到温度的部位灵活设置节点和热阻,其网络参数设置的合理性直接影响电机温度的计算精度,需要有丰富可靠的经验设置热路结构参数才能得到符合工程要求精度的分析结果。后来T型等效热网络模型被学者广泛应用,因为该模型以固定的热路形式表示电机的各个部件,能同时实现轴向和径向的热流分析和温度计算,能够简化网络的设计过程。目前,电机等效热网络模型存在的问题为:建立网络时,通常假设电机轴向呈对称结构,忽略了冷却介质的轴向温差对电机温度场分布的影响,认为电机温度以中间横截面呈对称分布,这在冷却介质轴向温差较大的电机中是不合理的。
发明内容
本发明的目的是为了克服上述热网络模型计算电机温度场存在的不足,改进热网络模型,提供一种轴向通风永磁同步电机温度场计算方法。本发明通过建立考虑到冷却介质温差影响的永磁同步电机三维等效热网络模型,提供更为精确的温度场计算方法。为了达到上述目的,本发明所采用的技术方案为:
一种轴向通风永磁同步电机温度场计算方法,包括下列步骤:
1)以传热学为基础分别建立电机各部件的热网络,电机部件包括外壳、定子轭、定子齿、定子绕组、绕组端部、定转子间气隙、永磁体、转子轭和转轴。
2)将各部件的热网络相互连接构成整个电机自身结构的热网络模型。
3)对电机中的冷却气流进行建模,将冷却气流模拟为由冷却气流量控制的热压源,设气流温度在端部和转子通风孔内呈线性变化,把各自中点位置的气流温度视为气流的平均温升:
①冷却气流从电机一侧端部流入,吸收掉绕组端部由定子轭,绕组端部,转子轭和永磁体通过热对流散发的热功率,使气流温度升高,得到第一个温升值;
②冷却气流吸收掉转子轭的部分损耗,使温度升高,得到第二个温升值;
③冷却气流在高温一侧端部吸收的由定子轭,绕组端部,转子轭和永磁体通过热对流散发的热功率,使其温度升高,得到第三个温升值。
用以上三个部位的温升值代表冷却气流在三个部位的平均温升,以热压源的形式在热网络模型中表示。
4)将表示冷却气流的热压源与电机自身结构热网络结合:
①电机两侧端部代表冷却气流的热压源分别通过一个对流热阻与外壳、定子轭、定子齿、永磁体、转子轭、转轴的侧面以及绕组端部相连接;
②电机转子通风孔内代表冷却气流的热压源通过一个对流热阻与转子轭的平均温度点相连接,最终形成完整的轴向通风永磁同步电机热网络模型。
5)列写节点方程求解各个节点温度。
本发明有益效果如下:
(1)所建立的三维热网络模型使电机自身结构形成的无源热网络与通风系统冷却气流形成一个整体,能够在同一方程组中求解,可以得到轴向非对称的温度场分布,弥补了传统的利用T型热网络模型建立电机热网络方法的缺点。
(2)本发明主要针对电机内部主要部件的温度进行计算,省略了外壳、端盖与轴承之间复杂的热连接,建模和计算较简单。
(3)冷却气流的建模方法是普适的,可以灵活地应用于更多结构类型的电机温度场计算当中。建模的方式也多种多样,不一定需要三个节点代表冷却气流温度,还可以设置更多节点代表冷却气流以得到更加精确的温度分布。
附图说明
图1为(a)永磁同步电机电机内的圆环柱体结构示意图及其对应的(b)T型等效热网络模型。
图2为定子齿热网络模型示意图,(a)为热网络,(b)为部件结构。
图3为外壳热网络模型示意图,(a)为热网络,(b)为部件结构。
图4为转轴热网络模型示意图,(a)为热网络,(b)为部件结构。
图5为定子绕组热网络模型示意图,(a)为热网络,(b)为部件结构。
图6为绕组端部热网络模型示意图,(a)为热网络,(b)为部件结构。
图7为定转子间气隙热网络模型示意图,(a)为热网络,(b)为部件结构。
图8为永磁同步电机中冷却气流温升示意图。
图9为永磁同步电机三维热网络模型。
图10为通风孔流速矢量图。
图11为通风孔冷却气流温度分布图。
图12为(a)为有冷却气流电机转子温度分布与(b)为无冷却气流电机转子温度分布的比较。
图13为电机定子的温度分布图。
图14为热网络模型求解结果与有限元仿真结果的比较。每组柱状图里,中间的柱状图为网络计算值,两边的柱状图分别为仿真最小值和最大值。
具体实施方式
本发明首先利用MATLAB软件编写程序建立永磁同步电机三维热网络模型,对网络进行求解,得到轴向非对称的温度场分布。然后利用有限元软件进行温度场和流体场耦合仿真,求解得到冷却气流的温度分布以及电机各部件的温度分布,比较两者得到的结果,证明了冷却通风系统模型以及整个三维热网络模型的正确性和可行性。
下面结合实施例对本发明进一步说明。本发明的具体实施步骤如下:
1)把永磁同步电机三维热网络模型一分为二,分别对冷却介质和电机自身结构进行建模。
2)按照电机各部件的结构特性将电机分为外壳、定子轭、定子齿、定子绕组、绕组端部、永磁体、转子轭、转轴定转子间气隙9个部分分成三组,其中,定子轭、转子轭为第一组;永磁体,定子齿为第二组;外壳、转轴、定子绕组、绕组端部和定转子间气隙为第三组。
3)第一组中的定子轭、转子轭均视为标准的圆环主体结构,可以用如图1所示的T型等效网络来表示,各个热阻的计算公式如下:
式中:kr和ka分别是圆环柱体径向和轴向的热导率(W/(m·℃);图1中:l为轴向长度;r1,r2分别是外圆环半径和内圆环半径;Taxial,left,Taxial,right分别为轴向两个截面的温度;Tradial,in,Tradial,out分别为内环和外环表面的温度;Tm表示整个柱体的平均温度;P是该柱体的生热功率。
4)第二组中的定子齿和永磁体并非圆环柱体,用等效截面积面积的方法在以上公式的基础上乘以一个系数,公式形式如下:
式中:Rcylindrical为把永磁体和定子齿按照圆环柱体的计算公式得到的热阻;Scylindrical为把永磁体和定子齿当做圆环柱体,用它们的内外环半径算得的横截面面积;Sreal为永磁体和定子齿各自实际的横截面面积。该公式实际上与电阻的定义类似,导体的横截面面积越大,电阻越小。
拿定子齿作为例子,如图2所示的定子齿尺寸,其T型等效网络六个热阻的计算公式如下:
式中,φp为齿距角,φe为单个齿的角度,如图所示。
5)第三组中的外壳、转轴、定子绕组和绕组端部,定转子间气隙结构特殊,下面分别给出这五个部分热网络建立过程:
(1)外壳:如图3所示,代表外壳的两个热阻R1,R2,其计算公式如下:
式中:α为外壳与外界空气的对流换热系数;S为外壳与外界的接触面积;hc为接触热阻系数。
(2)转轴:转轴内部无热量产生,将其热阻网络设计为图4所示,图中l为转轴与转子接触的轴向长度,la为转轴与端部气隙接触轴向长度的一半。R3为转轴与转子轭的接触热阻,R4,R5包含了径向边界到转轴平均温度点之间的热阻、平均温度点到轴向边界(此处为转轴与端部气隙接触面轴向的中点)的热阻、轴向边界与端部气体的对流热阻。
三个热阻计算公式如下:
式中:α为转轴与端部空气的对流换热系数;k为转轴的热传导系数;hc为接触热阻系数。
(3)定子绕组:定子绕组实际上由三部分组成:铜导体,导体绝缘漆,定子齿表面的绝缘层,假设只有铜导体进行轴向传热,另外两者只进行径向传热,导体与定子齿、定子轭、气隙的接触面积之比根据实际情况定义为2:1:1,导体径向传热系数根据经验定义为一个导热系数F与导体绝缘漆导热系数的乘积。其热网络模型如图5所示,四个热阻计算公式如下:
式中:di为导体绝缘层的厚度;kv为导体绝缘层的导热系数;ki为齿槽内绝缘层的导热系数;rw为导体的等效半径;l为导体长度;F为导体径向导热系数;n为槽数;Ac为导体的横截面积;kc为铜导体的热导率。
(4)绕组端部:绕组端部视为如图6所示,由环形结构和36个圆柱体两部分组成,环形结构代表绕组端部的导体及导体的绝缘,36个圆柱体代表定子绕组悬空部分。假设该环形结构和36个圆柱体的平均温度节点为同一个节点,且两者内部产生的铜耗均匀分布。图中的三个热阻计算公式如下:
式中:半径R视为定子槽的平均半径;re为绕组端部截面的半径;rw为悬空导体的等效半径;Ac为导体的截横面积;di为导体绝缘层的厚度;n为槽数;F为导体径向导热系数;kv为导体绝缘层的导热系数;kc为铜导体的热导率。
(5)定转子间气隙:定转子间气隙连接定子齿、定子绕组和永磁体,三个热阻均由气隙处的对流换热系数决定。其网络模型如图7所示,三个热阻计算公式如下:
式中:φp为齿距角;φe为单个齿的角度;r1为定子齿内半径;r2为永磁体外半径;hr为气隙对流换热系数;l为气隙的轴向长度。
至此,所有部件的热网络模型建立完成,连接各个部件,形成电机自身固体结构的热网络,尽量将串联的热阻相加合并为一个,简化热网络的节点数和支路数。
6)冷却气流可以模拟为由冷却气流量控制的热压源,与电路中的流控电压源形成类比。图8展示了永磁同步电机中冷却气流的温升,假设气流温度在端部和转子通风孔内呈线性变化,因此把各自中点位置的温度视为气流的平均温度。冷却气流从电机一侧端部流入,吸收掉绕组端部由定子轭,绕组端部,转子轭和永磁体通过热对流散发的热功率Pec1,使气流温度升高,温升为:
θ1end=Pec1/(ρcpq)=2RqPec1>
其中Rq具有热阻的量纲,其表达式为:
Rq=1/(2ρcpq)>
式中:ρ为气体密度(kg/m3);cp为气体热容量(J/kg·K);q为气流速度(m3/s).
各部位温度线性上升,低温一侧气体平均温升:
θ1=θ1end/2=RqPec1>
冷却气流吸收掉转子轭的部分损耗Pyr,使温度升高:
θ2end-θ1end=Pyr/(ρcpq)>
将(1)、(2)两式代入(4)式中,得:
θ2end=2RqPec1+2RqPyr>
与θ1的算法相似,分别得到转子通风孔内和高温一侧端部冷却气流的平均温升θ2与θ3:
θ2=(θ1end+θ2end)/2=Rq(2Pec1+Pyr)>
θ3=Rq(2Pec1+2Pyr+Pec2)>
式中:Pec2为冷却气流在高温一侧端部吸收的由定子轭,绕组端部,转子轭和永磁体通过热对流散发的热功率。冷却通风系统建立如上所述,用三个节点来表示冷却气流。类似于电路中的流控电压源,这里的冷却气流表示为气体流量q控制的位于三个节点的热压源。
7)得到了三个热压源的表达式,可以建立完整的永磁同步电机三维热网络模型,如图9所示。求解网络时要注意,加入冷却气流的热网络模型中存在热压源,方程形式如下:
式中:θ为节点温升列向量;G为节点热导纳矩阵;Pcf=[Pec1>yr>ec2]T为流入三个热压源正极的热流列向量;P=[P4>5 …>22]T为节点注入热流列向量,其中灰色节点为该部位注入的热功率数值,黑色节点没有功率注入,其值均为0。
热导纳矩阵的表达式如下:
式中:n为热网络的节点数,此处为22;1/Ri,j为节点i和节点j之间的热导。式(8)中热流列向量Pcf中有三个未知量,需要额外补充三个方程才能求得定解,将描述冷却气流的式(3)、(6)、(7)三式作为附加的三个方程,其矩阵形式如下:
简写为:
θcf=Rq·Pcf>
将(11)式代入(8)式中,得:
移项并合并同类项得到热网络各节点的温升:
至此,轴向通风永磁同步电机的三维热网络模型建立完成,建立了描述导纳矩阵、网络节点温度、冷却气流温度之间关系的节点方程,得到了轴向非对称的温度分布。
8)利用有限元仿真软件进行有限元仿真。电机模型为3000r/min,2对极,表贴式永磁同步电机,转子设有轴向通风孔,存在轴向冷却通风系统。仿真过程中做如下假设:冷却气流仅从转子通风孔内流入流出,不经过定转子间气隙,电机中各部件产生的热量均匀地分布于固体之中且忽略气隙风摩擦损耗,下表1为各部位施加的热载荷。
表1 各部件施加热载荷
9)建立永磁同步电机八分之一模型流体场和温度场的耦合场。得到仿真结果如图10至图13所示。
10)将三维热网络模型的求解结果与有限元耦合场仿真得到的结果相比较,如图14所示,两者结果非常接近。很好地证明了冷却通风系统模型以及整个三维热网络模型的可行性。可以在电机的优化设计过程中提供一些帮助。
机译: 包括轴向通风叶片的旋转电机的转子法兰的组装,该轴向通风叶片促进转子内部的轴向空气流动,并且涉及一种电机
机译: 一种用于风力发电缓慢旋转的永磁同步电机,在磁轮上具有作为固定件的配对件,该磁电机由两个盘形单元和轴向垂直延伸的对齐杆组成
机译: 一种用于在封闭的mantelgekuehlte机器中将敞开的druchzugsbeluefteten的外壳划分为轴向的电机的通风转换装置