CN114386238A - 一种考虑温度变化的异步电机半实物建模方法 - Google Patents

一种考虑温度变化的异步电机半实物建模方法 Download PDF

Info

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
Application number
CN202111503865.6A
Other languages
English (en)
Other versions
CN114386238B (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.)
CRRC Yongji Electric Co Ltd
Original Assignee
CRRC Yongji Electric Co Ltd
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 CRRC Yongji Electric Co Ltd filed Critical CRRC Yongji Electric Co Ltd
Priority to CN202111503865.6A priority Critical patent/CN114386238B/zh
Publication of CN114386238A publication Critical patent/CN114386238A/zh
Application granted granted Critical
Publication of CN114386238B publication Critical patent/CN114386238B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08Thermal 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

一种考虑温度变化的异步电机半实物建模方法
技术领域
本发明属于异步电机的仿真模型构建技术领域,涉及半实物仿真,具体为一种考虑温度变化的异步电机半实物仿真建模方法,为半实物仿真应用提供了一种更贴近实际的电机模型。
背景技术
半实物仿真是一种硬件在环仿真技术,通过快速建立被控对象或控制器仿真模型,对整个控制系统进行实时在线测试验证,可覆盖控制器产品全生命周期内的测试需求,尤其是电机缺相、短路、接地等许多真实试验环境无法验证的难题,在降低试验风险、节约成本的同时,大大缩短了产品开发周期,已经成为电机控制系统研发流程中不可或缺的重要环节。半实物仿真的精度主要取决于被控对象仿真模型的精度,仿真模型依据现有的等效数学方程搭建,而对于复杂的强耦合、非线性系统进行数学方程描述时往往会忽略一些本体结构、环境等因素的影响,与实际物理特性存在偏差,这就导致仿真模型搭建的不完善,也是影响半实物仿真精度的重要因素。
异步电机等效数学方程是研究异步电机动静态特性及其控制技术的理论基础,异步电机半实物仿真模型根据其数学方程建立。由于电机自身的复杂性,包括齿槽效应、磁路饱和、温度场效应等非线性因素,使得异步电机成为一个多变量、强耦合、非线性的高阶系统。为了简化电机模型及实现解耦控制,异步电机数学方程组通常转化到两相静止坐标系下,其电压矩阵方程为:
Figure BDA0003402666840000011
其磁链方程为:
Figure BDA0003402666840000012
此为电机统一理论,其核心在于线性变换,将复杂的非线性系统转换为具有常系数微分方程,而不考虑外部因素和特殊工况下的运行状态。使用此理论方程的前提假定条件是:电机结构对称;忽略磁饱和;忽略时间和空间谐波;忽略铁损、铜损等;忽略温度、湿度等环境因素的影响。
异步电机半实物模型使用电机统一理论数学方程搭建,表现为输入输出的线性变换关系,不能体现电机在运行过程中的参数变化特性,尤其是温度变化、磁饱和、转子时间常数等对电机参数的影响。所以,目前的异步电机半实物仿真模型不完善,不能解决电机运行过程中参数变化而引起的仿真结果与真实电机运行结果存在偏差的问题。
发明内容
本发明旨在解决目前的异步电机半实物仿真模型不完善,不能解决电机运行过程中参数变化而引起的仿真结果与真实电机运行结果存在偏差的技术问题,提供了一种考虑温度变化的异步电机半实物建模方法。
本发明解决其技术问题采用的技术手段是:一种考虑温度变化的异步电机半实物建模方法,包括以下步骤:
步骤一、建立定子绕组温度热网络模型:建立定子铁芯到电机机壳、电机机壳到环境温度两层热阻抗的电机定子温度热阻抗串联模型,设Rx_c、Cx_c分别是电机定子铁芯到电机机壳的热阻与热容,Rc_s、Cc_s分别是电机机壳到环境的热阻与热容,Ps为电机定子上产生的损耗功率,其等效为一个电流源;Tsx为定子绕组温度,Ts为环境温度,定子温度的时域模型表达式为:
Figure BDA0003402666840000021
式(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)可以看出,电机定子上产生的热损耗是与流过的电流值和电阻值相关的,在计算中电流作为实时值是可以获取的;电机定子铜绕组的电阻值在运行过程中会随着温度的变化而变化,规律为:
Figure BDA0003402666840000022
根据式(5)进一步将式(4)写为:
Figure BDA0003402666840000031
式(6)中,Tsb为冷态测量定子电阻值Rsb时的温度,Tsx(t)为定子电阻实时温度;
式(3)与式(6)组成了定子温度的时域模型,定子温度的时域模型的特点在于温度模型与损耗模型之间存在耦合,计算定子温度需要损耗值,计算损耗又需要定子温度值,实时定子电阻值采用式(5)即可计算;步骤一主要是通过计算电机定子温度,进而计算电子定子阻值Rs
步骤二、建立转子绕组温度热网络模型,转子绕组的温度建模采用与定子绕组温度相同的建模方法:设Trx为转子绕组温度,Tsx为定子绕组温度,Rr_s为转子绕组到定子绕组的热阻,Cr_s转子绕组到定子绕组的热容,对应步骤一,转子温度模型为:
Figure BDA0003402666840000032
式(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为时间;电机转子铜绕组的电阻值在运行过程中会随着温度的变化而变化,规律为:
Figure BDA0003402666840000033
根据式(9)进一步将式(8)写为:
Figure BDA0003402666840000034
式(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,从而使电机定、转子绕组温度热网络模型和异步电机电气模型实现数据同步交互。
本发明针对温度引起的异步电机定、转子阻值的变化,提出了阻值随温度变化的热阻热容并联模型,用迭代运算解决了定转子损耗与温度之间的耦合问题,同时提出了针对该类包含慢过程与快过程刚性系统的同步交互方法,保证了仿真运算的精度。
优选的,步骤四中,电机定、转子电流平均有效值的计算公式分别为:
Figure BDA0003402666840000041
Figure BDA0003402666840000042
式(11)和(12)中,k对应于当前运算时刻,即以h2为时间间隔的第k步运算。这是为了电机定子电流Is和电机转子电流Ir更加准确,使电机定、转子绕组温度热网络模型的各参数计算更加精确,进而更贴近实际中的电机运行状况。
优选的,电机定子绕组温度热网络模型相当于两个相串联的一阶惯性模型,其S域传递函数表达为:
Figure BDA0003402666840000043
式(13)中,ΔTsx_s为定子绕组相对于环境的温升,ΔTsx_c为定子绕组相对于电机机壳的温升,ΔTsc_s为电机机壳相对于环境的温升;
电机转子绕组温度热网络模型相当于一个一阶惯性模型,其S域传递函数表达为:
Figure BDA0003402666840000051
式(14)中,ΔTrx_sx为转子绕组到定子绕组的温升;
以h1为离散采样步长,将式(13)和式(14)采用Tustin法离散化为:
Figure BDA0003402666840000052
Figure BDA0003402666840000053
Figure BDA0003402666840000054
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为环境温度,定子温度的时域模型表达式为:
Figure BDA0003402666840000061
式(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)可以看出,电机定子上产生的热损耗是与流过的电流值和电阻值相关的,在计算中电流作为实时值是可以获取的;电机定子铜绕组的电阻值在运行过程中会随着温度的变化而变化,规律为:
Figure BDA0003402666840000071
根据式(5)进一步将式(4)写为:
Figure BDA0003402666840000072
式(6)中,Tsb为冷态测量定子电阻值Rsb时的温度,Tsx(t)为定子电阻实时温度;
式(3)与式(6)组成了定子温度的时域模型,定子温度的时域模型的特点在于温度模型与损耗模型之间存在耦合,计算定子温度需要损耗值,计算损耗又需要定子温度值,实时定子温度值采用式(5)即可计算;步骤一主要是通过计算电机定子温度,进而计算电子定子阻值Rs;步骤一提供了两个相互串联的一阶惯性模型,可以动态模拟异步电机运行过程中定子温度的变化;
步骤二、建立转子绕组温度热网络模型,如图3所示,转子绕组的温度建模采用与定子绕组温度相同的建模方法:设Trx为转子绕组温度,Tsx为定子绕组温度,Rr_s为转子绕组到定子绕组的热阻,Cr_s转子绕组到定子绕组的热容,对应步骤一,转子温度模型为:
Figure BDA0003402666840000073
式(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为时间;电机转子铜绕组的电阻值在运行过程中会随着温度的变化而变化,规律为:
Figure BDA0003402666840000081
根据式(9)进一步将式(8)写为:
Figure BDA0003402666840000082
式(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,从而使电机定、转子绕组温度热网络模型和异步电机电气模型实现数据同步交互;
其中,电机定、转子电流平均有效值的计算公式分别为:
Figure BDA0003402666840000083
Figure BDA0003402666840000091
式(11)和(12)中,k对应于当前运算时刻,即以h2为时间间隔的第k步运算;这是为了电机定子电流Is和电机转子电流Ir更加准确,使电机定、转子绕组温度热网络模型的各参数计算更加精确,进而更贴近实际中的电机运行状况;
电机定子绕组温度热网络模型相当于两个相串联的一阶惯性模型,其S域传递函数表达为:
Figure BDA0003402666840000092
式(13)中,ΔTsx_s为定子绕组相对于环境的温升,ΔTsx_c为定子绕组相对于电机机壳的温升,ΔTsc_s为电机机壳相对于环境的温升;
电机转子绕组温度热网络模型相当于一个一阶惯性模型,其S域传递函数表达为:
Figure BDA0003402666840000093
式(14)中,ΔTrx_sx为转子绕组相对于定子绕组的温升;
以h1离散采样步长为基础,将电机定、转子绕组温度热网络模型放入CPU运行,将异步电机电气模型放入FPGA中运行,采用Tustin法对式(13)和式(14)进行离散化处理来实现模型编程,将式(13)和式(14)离散化为:
Figure BDA0003402666840000094
Figure BDA0003402666840000095
Figure BDA0003402666840000096
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为环境温度,定子温度的时域模型表达式为:
Figure FDA0003402666830000011
式(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为时间;电机定子铜绕组的电阻值在运行过程中会随着温度的变化而变化,规律为:
Figure FDA0003402666830000012
根据式(5)进一步将式(4)写为:
Figure FDA0003402666830000013
式(6)中,Tsb为冷态测量定子电阻值Rsb时的温度,Tsx(t)为定子电阻实时温度;
式(3)与式(6)组成了定子温度的时域模型,定子温度的时域模型的特点在于温度模型与损耗模型之间存在耦合,计算定子温度需要损耗值,计算损耗又需要定子温度值,实时定子电阻值采用式(5)即可计算;
步骤二、建立转子绕组温度热网络模型:设Trx为转子绕组温度,Tsx为定子绕组温度,Rr_s为转子绕组到定子绕组的热阻,Cr_s转子绕组到定子绕组的热容,对应步骤一,转子温度模型为:
Figure FDA0003402666830000014
式(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为时间;电机转子铜绕组的电阻值在运行过程中会随着温度的变化而变化,规律为:
Figure FDA0003402666830000021
根据式(9)进一步将式(8)写为:
Figure FDA0003402666830000022
式(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,从而使电机定、转子绕组温度热网络模型和异步电机电气模型实现数据同步交互。
2.根据权利要求1所述的一种考虑温度变化的异步电机半实物建模方法,其特征在于,步骤四中,电机定、转子电流平均有效值的计算公式分别为:
Figure FDA0003402666830000031
Figure FDA0003402666830000032
式(11)和(12)中,k对应于当前运算时刻,即以h2为时间间隔的第k步运算。
3.根据权利要求2所述的一种考虑温度变化的异步电机半实物建模方法,其特征在于,电机定子绕组温度热网络模型相当于两个相串联的一阶惯性模型,其S域传递函数表达为:
Figure FDA0003402666830000033
式(13)中,ΔTsx_s为定子绕组相对于环境的温升,ΔTsx_c为定子绕组相对于电机机壳的温升,ΔTsc_s为电机机壳相对于环境的温升;
电机转子绕组温度热网络模型相当于一个一阶惯性模型,其S域传递函数表达为:
Figure FDA0003402666830000034
式(14)中,ΔTrx_sx为转子绕组相对于定子绕组的温升;
以h1为离散采样步长,将式(13)和式(14)采用Tustin法离散化为:
Figure FDA0003402666830000035
Figure FDA0003402666830000036
Figure FDA0003402666830000037
Tsx(k)=ΔTsx_c(k)+ΔTsc_s(k)+Ts(k) (18),
Trx(k)=ΔTrx_sx(k)+Tsx(k) (19)。
4.根据权利要求3所述的一种考虑温度变化的异步电机半实物建模方法,其特征在于,步骤四中,将电机定、转子绕组温度热网络模型放入CPU运行,将异步电机电气模型放入FPGA中运行。
CN202111503865.6A 2021-12-09 2021-12-09 一种考虑温度变化的异步电机半实物建模方法 Active CN114386238B (zh)

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)

* Cited by examiner, † Cited by third party
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 盖耀辉 一种永磁电机转子永磁体温度的实时监测方法及模型

Patent Citations (4)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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