CN114386238A - 一种考虑温度变化的异步电机半实物建模方法 - Google Patents
一种考虑温度变化的异步电机半实物建模方法 Download PDFInfo
- Publication number
- CN114386238A CN114386238A CN202111503865.6A CN202111503865A CN114386238A CN 114386238 A CN114386238 A CN 114386238A CN 202111503865 A CN202111503865 A CN 202111503865A CN 114386238 A CN114386238 A CN 114386238A
- Authority
- CN
- China
- Prior art keywords
- motor
- stator
- temperature
- model
- rotor
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 54
- 230000008859 change Effects 0.000 title claims abstract description 39
- 238000004804 winding Methods 0.000 claims abstract description 97
- 230000008569 process Effects 0.000 claims abstract description 23
- 230000003993 interaction Effects 0.000 claims abstract description 14
- 230000001360 synchronised effect Effects 0.000 claims abstract description 14
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 claims description 32
- 229910052802 copper Inorganic materials 0.000 claims description 32
- 239000010949 copper Substances 0.000 claims description 32
- 238000004364 calculation method Methods 0.000 claims description 18
- 239000004576 sand Substances 0.000 claims description 16
- 238000005070 sampling Methods 0.000 claims description 12
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical group [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 claims description 10
- 238000012546 transfer Methods 0.000 claims description 8
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 6
- 238000004088 simulation Methods 0.000 abstract description 22
- 238000010168 coupling process Methods 0.000 abstract description 7
- 238000005859 coupling reaction Methods 0.000 abstract description 7
- 230000008878 coupling Effects 0.000 abstract description 6
- 238000005516 engineering process Methods 0.000 abstract description 6
- 238000010276 construction Methods 0.000 abstract description 2
- 238000012360 testing method Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 230000003068 static effect Effects 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 3
- 230000004907 flux Effects 0.000 description 2
- 230000017525 heat dissipation Effects 0.000 description 2
- 238000010438 heat treatment Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000012356 Product development Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000005669 field effect Effects 0.000 description 1
- 229910052742 iron Inorganic materials 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Control Of Electric Motors In General (AREA)
Abstract
本发明属于异步电机控制系统的仿真模型构建技术领域,具体为一种考虑温度变化的异步电机半实物建模方法,解决了背景技术中的技术问题,其包括分别建立定子绕组温度热网络模型、转子绕组温度热网络模型和异步电机电气特性模型,对各模型进行离散化处理,选定热模型离散化运算步长为h1,选定电气特性模型离散化步长为h2,将定子、转子热网络模型放入CPU运行,将电气特性模型放入FPGA运行,实现CPU与FPGA之间的数据交互及同步。本发明提出了阻值随温度变化的热阻热容并联模型,用迭代运算解决了定转子损耗与温度之间的耦合问题,同时提出了针对该类包含慢过程与快过程刚性系统的同步交互方法,保证了仿真运算的精度。
Description
技术领域
本发明属于异步电机的仿真模型构建技术领域,涉及半实物仿真,具体为一种考虑温度变化的异步电机半实物仿真建模方法,为半实物仿真应用提供了一种更贴近实际的电机模型。
背景技术
半实物仿真是一种硬件在环仿真技术,通过快速建立被控对象或控制器仿真模型,对整个控制系统进行实时在线测试验证,可覆盖控制器产品全生命周期内的测试需求,尤其是电机缺相、短路、接地等许多真实试验环境无法验证的难题,在降低试验风险、节约成本的同时,大大缩短了产品开发周期,已经成为电机控制系统研发流程中不可或缺的重要环节。半实物仿真的精度主要取决于被控对象仿真模型的精度,仿真模型依据现有的等效数学方程搭建,而对于复杂的强耦合、非线性系统进行数学方程描述时往往会忽略一些本体结构、环境等因素的影响,与实际物理特性存在偏差,这就导致仿真模型搭建的不完善,也是影响半实物仿真精度的重要因素。
异步电机等效数学方程是研究异步电机动静态特性及其控制技术的理论基础,异步电机半实物仿真模型根据其数学方程建立。由于电机自身的复杂性,包括齿槽效应、磁路饱和、温度场效应等非线性因素,使得异步电机成为一个多变量、强耦合、非线性的高阶系统。为了简化电机模型及实现解耦控制,异步电机数学方程组通常转化到两相静止坐标系下,其电压矩阵方程为:
其磁链方程为:
此为电机统一理论,其核心在于线性变换,将复杂的非线性系统转换为具有常系数微分方程,而不考虑外部因素和特殊工况下的运行状态。使用此理论方程的前提假定条件是:电机结构对称;忽略磁饱和;忽略时间和空间谐波;忽略铁损、铜损等;忽略温度、湿度等环境因素的影响。
异步电机半实物模型使用电机统一理论数学方程搭建,表现为输入输出的线性变换关系,不能体现电机在运行过程中的参数变化特性,尤其是温度变化、磁饱和、转子时间常数等对电机参数的影响。所以,目前的异步电机半实物仿真模型不完善,不能解决电机运行过程中参数变化而引起的仿真结果与真实电机运行结果存在偏差的问题。
发明内容
本发明旨在解决目前的异步电机半实物仿真模型不完善,不能解决电机运行过程中参数变化而引起的仿真结果与真实电机运行结果存在偏差的技术问题,提供了一种考虑温度变化的异步电机半实物建模方法。
本发明解决其技术问题采用的技术手段是:一种考虑温度变化的异步电机半实物建模方法,包括以下步骤:
步骤一、建立定子绕组温度热网络模型:建立定子铁芯到电机机壳、电机机壳到环境温度两层热阻抗的电机定子温度热阻抗串联模型,设Rx_c、Cx_c分别是电机定子铁芯到电机机壳的热阻与热容,Rc_s、Cc_s分别是电机机壳到环境的热阻与热容,Ps为电机定子上产生的损耗功率,其等效为一个电流源;Tsx为定子绕组温度,Ts为环境温度,定子温度的时域模型表达式为:
式(3)中,τc_s和τx_c均为惯性时间,τc_s=Rc_s×Cc_s,τx_c=Rx_c×Cx_c,Ps(t)为电机定子铜耗功率,电机定子温度建模的方法主要是根据电机的发热机理,将热量在电机上的散热过程用热力学理论表达出来,在热力学理论中,选用热阻抗理论来描述定转子温度,步骤一中提出了考虑定子铁芯到电机机壳、电机机壳到环境温度两层热阻抗的电机定子温度热阻抗串联模型;在计算电机定子上产生的损耗之前,首先分析电机上的能量损耗,电机在带载运行过程中通常会产生铜耗、铁耗、杂散损耗、机械损耗等几种损耗,其中电机定子铁芯上产生的损耗主要为铜耗,故电机定子铜耗功率等同于电机定子铁芯上产生的损耗功率,电机定子铜耗功率Ps(t)的计算公式为:
Ps(t)=3Is(t)2×Rs(t) (4),
式(4)中,Is(t)为电机定子电流有效值,Rs(t)为电机定子电阻值,t为时间,从式(4)可以看出,电机定子上产生的热损耗是与流过的电流值和电阻值相关的,在计算中电流作为实时值是可以获取的;电机定子铜绕组的电阻值在运行过程中会随着温度的变化而变化,规律为:
根据式(5)进一步将式(4)写为:
式(6)中,Tsb为冷态测量定子电阻值Rsb时的温度,Tsx(t)为定子电阻实时温度;
式(3)与式(6)组成了定子温度的时域模型,定子温度的时域模型的特点在于温度模型与损耗模型之间存在耦合,计算定子温度需要损耗值,计算损耗又需要定子温度值,实时定子电阻值采用式(5)即可计算;步骤一主要是通过计算电机定子温度,进而计算电子定子阻值Rs;
步骤二、建立转子绕组温度热网络模型,转子绕组的温度建模采用与定子绕组温度相同的建模方法:设Trx为转子绕组温度,Tsx为定子绕组温度,Rr_s为转子绕组到定子绕组的热阻,Cr_s转子绕组到定子绕组的热容,对应步骤一,转子温度模型为:
式(7)中,Pr(t)为电机转子铜耗功率,惯性时间τr_s=Rr_s×Cr_s;电机转子铜耗功率Pr(t)的计算公式为:
Pr(t)=3Ir(t)2×Rr(t) (8),
式(8)中,Ir(t)为电机转子电流有效值,Rr(t)为电机转子电阻值,t为时间;电机转子铜绕组的电阻值在运行过程中会随着温度的变化而变化,规律为:
根据式(9)进一步将式(8)写为:
式(10)中,Ir为转子绕组上流过的电流有效值,Trb为冷态测量转子绕组电阻值Rrb时的温度;步骤二主要是通过计算电机转子温度,进而计算电子转子阻值Rr;
式(7)与式(10)组成了转子温度的时域模型,转子温度的时域模型的特点在于温度模型与损耗模型之间存在耦合,计算转子温度需要损耗值,计算损耗又需要转子温度值,实时转子电阻值采用式(9)即可计算;
步骤三、构建异步电机电气模型,该电气模型属于现有模型,如背景技术中提到的内容:异步电机电气模型表达式为式(1)和式(2)所示;步骤三主要是异步电机在两相静止坐标下电压和磁链的运算;
步骤四、将步骤一、二中的电机定、转子绕组温度热网络模型和步骤三中的异步电机电气模型进行离散运算并实现数据同步交互:
电机定、转子绕组温度热网络模型是热模型,惯性较大,特性变化缓慢,而异步电机电气模型惯性较小,变化迅速,两种模型组成的系统属于刚性系统,为了准确进行离散运算,减少运算误差,节省运算资源,两部分模型分别以不同的运算步长计算,故需要进行同步数据交互;电机定、转子绕组温度热网络模型采用h1离散采样步长进行运算,异步电机电气模型采用h2离散采样步长进行运算,且h1是h2的整数倍,即h1=n*h2,运算步长具体数值根据各部分模型的运算误差(包含截断误差及舍入误差)最小及运算量尽可能少来综合确定,h1和h2之间的约束主要是为了之后便于同步算法的实现;电机定、转子绕组温度热网络模型以h1为时间间隔向电气模型传递电机定子阻值Rs和电机转子阻值Rr,异步电机电气模型所引用的电阻值每间隔n个h2运算离散采样步长后更新一次,并以n个h2为时间间隔向电机定、转子绕组温度热网络模型传递电机定子电流Is和电机转子电流Ir,其中异步电机电气模型每一运算步都计算过去h1时间段内的电机定、转子电流平均有效值,电机定子绕组温度热网络模型将电机定子电流Is带入式(6)即可计算出电机定子铜耗功率Ps(t),最终通过式(3)和式(5)计算出电机定子阻值Rs,电机转子绕组温度热网络模型将电机转子电流Ir带入式(10)即可计算出电机定子铜耗功率Pr(t),最终通过式(7)和式(9)计算出电机定子阻值Rr,从而使电机定、转子绕组温度热网络模型和异步电机电气模型实现数据同步交互。
本发明针对温度引起的异步电机定、转子阻值的变化,提出了阻值随温度变化的热阻热容并联模型,用迭代运算解决了定转子损耗与温度之间的耦合问题,同时提出了针对该类包含慢过程与快过程刚性系统的同步交互方法,保证了仿真运算的精度。
优选的,步骤四中,电机定、转子电流平均有效值的计算公式分别为:
式(11)和(12)中,k对应于当前运算时刻,即以h2为时间间隔的第k步运算。这是为了电机定子电流Is和电机转子电流Ir更加准确,使电机定、转子绕组温度热网络模型的各参数计算更加精确,进而更贴近实际中的电机运行状况。
优选的,电机定子绕组温度热网络模型相当于两个相串联的一阶惯性模型,其S域传递函数表达为:
式(13)中,ΔTsx_s为定子绕组相对于环境的温升,ΔTsx_c为定子绕组相对于电机机壳的温升,ΔTsc_s为电机机壳相对于环境的温升;
电机转子绕组温度热网络模型相当于一个一阶惯性模型,其S域传递函数表达为:
式(14)中,ΔTrx_sx为转子绕组到定子绕组的温升;
以h1为离散采样步长,将式(13)和式(14)采用Tustin法离散化为:
Tsx(k)=ΔTsx_c(k)+ΔTsc_s(k)+Ts(k) (18),
Trx(k)=ΔTrx_sx(k)+Tsx(k) (19)。
优选的,步骤四中,将电机定、转子绕组温度热网络模型放入CPU运行,将异步电机电气模型放入FPGA中运行。通过上述离散化模型公式处理来实现CPU和FPGA编程。异步电机电气模型计算的平均电流有效值存储于FPGA的寄存器内,每h2的步长更新一次,并以n个h2为时间间隔向电机定、转子绕组温度热网络模型传递电机定子电流Is和电机转子电流Ir。热模型每隔h1步长时间从FPGA的寄存器读取一次平均电流有效值,该方法能够保证CPU和FPGA之间数据交互及同步。
本发明的有益效果是:本发明针对温度引起的异步电机定、转子阻值的变化,提出了阻值随温度变化的热阻热容并联模型,用迭代运算解决了定转子损耗与温度之间的耦合问题,同时提出了针对该类包含慢过程与快过程刚性系统的同步交互方法,保证了仿真运算的精度;采用该模型描述影响异步电机控制输出的转子电阻变化,同时做转子温度的状态监测;采用该模型能研究更多的运行工况,得到全工况范围内的定转子温度模型,采用算法替代定子温度传感器,降低硬件成本;在半实物仿真平台模拟电机温升过程,为半实物仿真平台研究电机温度变化对控制系统的影响提供支撑,本专利提出的建模方法,不仅适用于异步电机,同样适用于永磁同步电机。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案,下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明所述一种考虑温度变化的异步电机半实物建模方法中建立的异步电机模型示意图。
图2为本发明所述定子绕组温度热网络模型的结构示意图。
图3为本发明所述转子绕组温度热网络模型的结构示意图。
图4为本发明所述一种考虑温度变化的异步电机半实物建模方法中建立的异步电机模型在温升过程中的定子电阻值变化示意图。
图5为本发明所述一种考虑温度变化的异步电机半实物建模方法中建立的异步电机模型在温升过程中的定子损耗功率变化示意图。
图6为本发明所述一种考虑温度变化的异步电机半实物建模方法中计算的定子温度与试验所采实际温度的对比示意图。
具体实施方式
下面将结合附图1至图6对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
一种考虑温度变化的异步电机半实物建模方法,如图1所示,包括以下步骤:
步骤一、建立定子绕组温度热网络模型,如图2所示:建立定子铁芯到电机机壳、电机机壳到环境温度两层热阻抗的电机定子温度热阻抗串联模型,设Rx_c、Cx_c分别是电机定子铁芯到电机机壳的热阻与热容,Rc_s、Cc_s分别是电机机壳到环境的热阻与热容,Ps为电机定子上产生的损耗功率,其等效为一个电流源;Tsx为定子绕组温度,Ts为环境温度,定子温度的时域模型表达式为:
式(3)中,τc_s和τx_c均为惯性时间,τc_s=Rc_s×Cc_s,τx_c=Rx_c×Cx_c,Ps(t)为电机定子铜耗功率,电机定子温度建模的方法主要是根据电机的发热机理,将热量在电机上的散热过程用热力学理论表达出来,在热力学理论中,选用热阻抗理论来描述定转子温度,步骤一中提出了考虑定子铁芯到电机机壳、电机机壳到环境温度两层热阻抗的电机定子温度热阻抗串联模型;在计算电机定子上产生的损耗之前,首先分析电机上的能量损耗,电机在带载运行过程中通常会产生铜耗、铁耗、杂散损耗、机械损耗等几种损耗,其中电机定子铁芯上产生的损耗主要为铜耗,故电机定子铜耗功率等同于电机定子铁芯上产生的损耗功率,电机定子铜耗功率Ps(t)的计算公式为:
Ps(t)=3Is(t)2×Rs(t) (4),
式(4)中,Is(t)为电机定子电流有效值,Rs(t)为电机定子电阻值,t为时间,从式(4)可以看出,电机定子上产生的热损耗是与流过的电流值和电阻值相关的,在计算中电流作为实时值是可以获取的;电机定子铜绕组的电阻值在运行过程中会随着温度的变化而变化,规律为:
根据式(5)进一步将式(4)写为:
式(6)中,Tsb为冷态测量定子电阻值Rsb时的温度,Tsx(t)为定子电阻实时温度;
式(3)与式(6)组成了定子温度的时域模型,定子温度的时域模型的特点在于温度模型与损耗模型之间存在耦合,计算定子温度需要损耗值,计算损耗又需要定子温度值,实时定子温度值采用式(5)即可计算;步骤一主要是通过计算电机定子温度,进而计算电子定子阻值Rs;步骤一提供了两个相互串联的一阶惯性模型,可以动态模拟异步电机运行过程中定子温度的变化;
步骤二、建立转子绕组温度热网络模型,如图3所示,转子绕组的温度建模采用与定子绕组温度相同的建模方法:设Trx为转子绕组温度,Tsx为定子绕组温度,Rr_s为转子绕组到定子绕组的热阻,Cr_s转子绕组到定子绕组的热容,对应步骤一,转子温度模型为:
式(7)中,Pr(t)为电机转子铜耗功率,惯性时间τr_s=Rr_s×Cr_s;电机转子铜耗功率Pr(t)的计算公式为:
Pr(t)=3Ir(t)2×Rr(t) (8),
式(8)中,Ir(t)为电机转子电流有效值,Rr(t)为电机转子电阻值,t为时间;电机转子铜绕组的电阻值在运行过程中会随着温度的变化而变化,规律为:
根据式(9)进一步将式(8)写为:
式(10)中,Ir为转子绕组上流过的电流有效值,Trb为冷态测量转子绕组电阻值Rrb时的温度;步骤二主要是通过计算电机转子温度,进而计算电子转子阻值Rr;
式(7)与式(10)组成了转子温度的时域模型,转子温度的时域模型的特点在于温度模型与损耗模型之间存在耦合,计算转子温度需要损耗值,计算损耗又需要转子温度值,实时转子电阻值采用式(9)即可计算;步骤二提供了一阶惯性模型,可以动态模拟异步电机运行过程中转子温度的变化;
步骤三、构建异步电机电气模型,该电气模型属于现有模型,如背景技术中提到的内容:异步电机电气模型表达式为式(1)和式(2)所示;步骤三主要是异步电机在两相静止坐标下电压和磁链的运算;
步骤四、将步骤一、二中的电机定、转子绕组温度热网络模型和步骤三中的异步电机电气模型进行离散运算并实现数据同步交互:
电机定、转子绕组温度热网络模型是热模型,惯性较大,特性变化缓慢,而异步电机电气模型惯性较小,变化迅速,两种模型组成的系统属于刚性系统,为了准确进行离散运算,减少运算误差,节省运算资源,故需要进行同步交互;电机定、转子绕组温度热网络模型采用h1离散采样步长进行运算,异步电机电气模型采用h2离散采样步长进行运算,且h1是h2的整数倍,即h1=n*h2,运算步长具体数值根据各部分模型的运算误差(包含截断误差及舍入误差)最小及运算量尽可能少来综合确定,h1和h2之间的约束主要是为了之后便于同步算法的实现;电机定、转子绕组温度热网络模型以h1为时间间隔向电气模型传递电机定子阻值Rs和电机转子阻值Rr,异步电机电气模型所引用的电阻值每间隔n个h2运算离散采样步长后更新一次,并以n个h2为时间间隔向电机定、转子绕组温度热网络模型传递电机定子电流Is和电机转子电流Ir,其中异步电机电气模型每一运算步都计算过去h1时间段内的电机定、转子电流平均有效值,电机定子绕组温度热网络模型将电机定子电流Is带入式(6)即可计算出电机定子铜耗功率Ps(t),最终通过式(3)和式(5)计算出电机定子阻值Rs,电机转子绕组温度热网络模型将电机转子电流Ir带入式(10)即可计算出电机定子铜耗功率Pr(t),最终通过式(7)和式(9)计算出电机定子阻值Rr,从而使电机定、转子绕组温度热网络模型和异步电机电气模型实现数据同步交互;
其中,电机定、转子电流平均有效值的计算公式分别为:
式(11)和(12)中,k对应于当前运算时刻,即以h2为时间间隔的第k步运算;这是为了电机定子电流Is和电机转子电流Ir更加准确,使电机定、转子绕组温度热网络模型的各参数计算更加精确,进而更贴近实际中的电机运行状况;
电机定子绕组温度热网络模型相当于两个相串联的一阶惯性模型,其S域传递函数表达为:
式(13)中,ΔTsx_s为定子绕组相对于环境的温升,ΔTsx_c为定子绕组相对于电机机壳的温升,ΔTsc_s为电机机壳相对于环境的温升;
电机转子绕组温度热网络模型相当于一个一阶惯性模型,其S域传递函数表达为:
式(14)中,ΔTrx_sx为转子绕组相对于定子绕组的温升;
以h1离散采样步长为基础,将电机定、转子绕组温度热网络模型放入CPU运行,将异步电机电气模型放入FPGA中运行,采用Tustin法对式(13)和式(14)进行离散化处理来实现模型编程,将式(13)和式(14)离散化为:
Tsx(k)=ΔTsx_c(k)+ΔTsc_s(k)+Ts(k) (18),
Trx(k)=ΔTrx_sx(k)+Tsx(k) (19),
离散化处理是为了便于计算机编程。异步电机电气模型计算的平均电流有效值存储于FPGA的寄存器内,每h2的步长更新一次,并以n个h2为时间间隔向电机定、转子绕组温度热网络模型传递电机定子电流Is和电机转子电流Ir。步骤四提出了异步电机热模型与电气模型耦合的数学模型;提出了使用CPU处理器和FPGA联合建模方式,将电机热模型和电气动态模型分别置于CPU处理器和FPGA的解算模式,解决了刚性系统解算精度问题;提出了CPU模型与FPGA模型之间数据交互同步的方法;提出了热模型的离散迭代计算方法;提出了两部分模型运算步长的选择依据,以总运算误差(阶段误差和舍入误差)最小为准则。
采用上述方法对某型号电机进行建模,并以实际温升试验数据来验证所述模型的有效性。通过图4、图5和图6可以看出,本专利所述模型计算出的定子温度与试验所采数据非常吻合,验证了所述异步电机热特性与电气特性耦合模型的有效性。
最后应说明的是:以上各实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述各实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
Claims (4)
1.一种考虑温度变化的异步电机半实物建模方法,其特征在于,包括以下步骤:
步骤一、建立定子绕组温度热网络模型:建立定子铁芯到电机机壳、电机机壳到环境温度两层热阻抗的电机定子温度热阻抗串联模型,设Rx_c、Cx_c分别是电机定子铁芯到电机机壳的热阻与热容,Rc_s、Cc_s分别是电机机壳到环境的热阻与热容,Ps为电机定子上产生的损耗功率,其等效为一个电流源;Tsx为定子绕组温度,Ts为环境温度,定子温度的时域模型表达式为:
式(3)中,τc_s和τx_c均为惯性时间,τc_s=Rc_s×Cc_s,τx_c=Rx_c×Cx_c,Ps(t)为电机定子铜耗功率,电机定子铁芯上产生的损耗主要为铜耗,故电机定子铜耗功率等同于电机定子铁芯上产生的损耗功率,电机定子铜耗功率Ps(t)的计算公式为:
Ps(t)=3Is(t)2×Rs(t) (4),
式(4)中,Is(t)为电机定子电流有效值,Rs(t)为电机定子电阻值,t为时间;电机定子铜绕组的电阻值在运行过程中会随着温度的变化而变化,规律为:
根据式(5)进一步将式(4)写为:
式(6)中,Tsb为冷态测量定子电阻值Rsb时的温度,Tsx(t)为定子电阻实时温度;
式(3)与式(6)组成了定子温度的时域模型,定子温度的时域模型的特点在于温度模型与损耗模型之间存在耦合,计算定子温度需要损耗值,计算损耗又需要定子温度值,实时定子电阻值采用式(5)即可计算;
步骤二、建立转子绕组温度热网络模型:设Trx为转子绕组温度,Tsx为定子绕组温度,Rr_s为转子绕组到定子绕组的热阻,Cr_s转子绕组到定子绕组的热容,对应步骤一,转子温度模型为:
式(7)中,Pr(t)为电机转子铜耗功率,惯性时间τr_s=Rr_s×Cr_s;电机转子铜耗功率Pr(t)的计算公式为:
Pr(t)=3Ir(t)2×Rr(t) (8),
式(8)中,Ir(t)为电机转子电流有效值,Rr(t)为电机转子电阻值,t为时间;电机转子铜绕组的电阻值在运行过程中会随着温度的变化而变化,规律为:
根据式(9)进一步将式(8)写为:
式(10)中,Ir为转子绕组上流过的电流有效值,Trb为冷态测量转子绕组电阻值Rrb时的温度;
式(7)与式(10)组成了转子温度的时域模型,转子温度的时域模型的特点在于温度模型与损耗模型之间存在耦合,计算转子温度需要损耗值,计算损耗又需要转子温度值,实时转子电阻值采用式(9)即可计算;
步骤三、构建异步电机电气模型:异步电机电气模型表达式为式(1)和式(2)所示;
步骤四、将步骤一、二中的电机定、转子绕组温度热网络模型和步骤三中的异步电机电气模型进行离散运算并实现数据同步交互:
电机定、转子绕组温度热网络模型采用h1离散采样步长进行运算,异步电机电气模型采用h2离散采样步长进行运算,且h1是h2的整数倍,即h1=n*h2,电机定、转子绕组温度热网络模型以h1为时间间隔向电气模型传递电机定子阻值Rs和电机转子阻值Rr,异步电机电气模型所引用的电阻值每间隔n个h2运算离散采样步长后更新一次,并以n个h2为时间间隔向电机定、转子绕组温度热网络模型传递电机定子电流Is和电机转子电流Ir,其中异步电机电气模型每一运算步都计算过去h1时间段内的电机定、转子电流平均有效值,电机定子绕组温度热网络模型将电机定子电流Is带入式(6)即可计算出电机定子铜耗功率Ps(t),最终通过式(3)和式(5)计算出电机定子阻值Rs,电机转子绕组温度热网络模型将电机转子电流Ir带入式(10)即可计算出电机定子铜耗功率Pr(t),最终通过式(7)和式(9)计算出电机定子阻值Rr,从而使电机定、转子绕组温度热网络模型和异步电机电气模型实现数据同步交互。
3.根据权利要求2所述的一种考虑温度变化的异步电机半实物建模方法,其特征在于,电机定子绕组温度热网络模型相当于两个相串联的一阶惯性模型,其S域传递函数表达为:
式(13)中,ΔTsx_s为定子绕组相对于环境的温升,ΔTsx_c为定子绕组相对于电机机壳的温升,ΔTsc_s为电机机壳相对于环境的温升;
电机转子绕组温度热网络模型相当于一个一阶惯性模型,其S域传递函数表达为:
式(14)中,ΔTrx_sx为转子绕组相对于定子绕组的温升;
以h1为离散采样步长,将式(13)和式(14)采用Tustin法离散化为:
Tsx(k)=ΔTsx_c(k)+ΔTsc_s(k)+Ts(k) (18),
Trx(k)=ΔTrx_sx(k)+Tsx(k) (19)。
4.根据权利要求3所述的一种考虑温度变化的异步电机半实物建模方法,其特征在于,步骤四中,将电机定、转子绕组温度热网络模型放入CPU运行,将异步电机电气模型放入FPGA中运行。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111503865.6A CN114386238B (zh) | 2021-12-09 | 2021-12-09 | 一种考虑温度变化的异步电机半实物建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111503865.6A CN114386238B (zh) | 2021-12-09 | 2021-12-09 | 一种考虑温度变化的异步电机半实物建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114386238A true CN114386238A (zh) | 2022-04-22 |
CN114386238B CN114386238B (zh) | 2024-04-19 |
Family
ID=81196181
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111503865.6A Active CN114386238B (zh) | 2021-12-09 | 2021-12-09 | 一种考虑温度变化的异步电机半实物建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114386238B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110050141A1 (en) * | 2009-08-31 | 2011-03-03 | Gm Global Technology Operations, Inc. | Electric motor stator winding temperature estimation |
CN107017666A (zh) * | 2017-04-06 | 2017-08-04 | 中国人民解放军后勤工程学院 | 小型单相同步发电机并列运行系统有限元模型建立方法 |
CN108334687A (zh) * | 2018-01-29 | 2018-07-27 | 扬州大学 | 一种大中型电机运行温升可靠度的预测方法 |
CN112234911A (zh) * | 2020-09-21 | 2021-01-15 | 盖耀辉 | 一种永磁电机转子永磁体温度的实时监测方法及模型 |
-
2021
- 2021-12-09 CN CN202111503865.6A patent/CN114386238B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110050141A1 (en) * | 2009-08-31 | 2011-03-03 | Gm Global Technology Operations, Inc. | Electric motor stator winding temperature estimation |
CN107017666A (zh) * | 2017-04-06 | 2017-08-04 | 中国人民解放军后勤工程学院 | 小型单相同步发电机并列运行系统有限元模型建立方法 |
CN108334687A (zh) * | 2018-01-29 | 2018-07-27 | 扬州大学 | 一种大中型电机运行温升可靠度的预测方法 |
CN112234911A (zh) * | 2020-09-21 | 2021-01-15 | 盖耀辉 | 一种永磁电机转子永磁体温度的实时监测方法及模型 |
Non-Patent Citations (1)
Title |
---|
陈起旭;梁得亮;徐俊;曹秉刚;: "盘式交流永磁同步电机温升影响因素研究", 电机与控制学报, no. 03, 12 January 2018 (2018-01-12) * |
Also Published As
Publication number | Publication date |
---|---|
CN114386238B (zh) | 2024-04-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107086606B (zh) | 一种配电网综合等效异步电动机负荷模型参数辨识方法 | |
Chen et al. | Digital hardware emulation of universal machine and universal line models for real-time electromagnetic transient simulation | |
Demiray | Simulation of power system dynamics using dynamic phasor models | |
CN111709107B (zh) | 一种含大规模电力电子设备的电力网络仿真方法及系统 | |
CN104052363B (zh) | 自动化电机控制 | |
Yang et al. | Real-time FPGA-RTDS co-simulator for power systems | |
CN103678900B (zh) | 一种用于区域配电网实时仿真的网络解耦计算方法 | |
CN102609575A (zh) | 一种基于隐式数值积分的电力系统暂态稳定仿真方法 | |
Stifter et al. | Steady-state co-simulation with PowerFactory | |
Rehtanz et al. | Real-time and co-simulations for the development of power system monitoring, control and protection | |
Gao et al. | A model reduction perspective on thermal models for induction machine overload relays | |
CN111881611A (zh) | 电机物理场的仿真处理方法及装置 | |
CN112234911A (zh) | 一种永磁电机转子永磁体温度的实时监测方法及模型 | |
WO2006054784A1 (ja) | モータモデルの演算方法、モータシミュレーション方法、モータシミュレーション装置、モータモデル演算プログラム、シミュレーション方法及びシミュレーションプログラム | |
Haugdal et al. | An open source power system simulator in python for efficient prototyping of wampac applications | |
CN102664397A (zh) | 一种基于隐式精细数值积分的电力系统暂态稳定仿真方法 | |
WO2018102720A1 (en) | System and method for a fast power network simulator | |
CN114254557A (zh) | 一种基于数字孪生的变压器在线监测方法 | |
Najafabadi et al. | Real time estimation of sensitive parameters of composite power system load model | |
CN104950261B (zh) | 电池的硬件在环仿真测试方法及系统 | |
CN114386238A (zh) | 一种考虑温度变化的异步电机半实物建模方法 | |
CN113221298A (zh) | 一种机电暂态过程的仿真方法及系统 | |
CN116595395A (zh) | 一种基于深度学习的逆变器输出电流预测方法及系统 | |
Mongrain et al. | Multi‐platform real‐time microgrid simulation testbed with hierarchical control of distributed energy resources featuring energy storage balancing | |
CN115422808A (zh) | 一种基于Krylov子空间的变压器温度场模型降阶方法 |
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 |