CN113676108B - 一种感应电机转子磁链的观测方法 - Google Patents
一种感应电机转子磁链的观测方法 Download PDFInfo
- Publication number
- CN113676108B CN113676108B CN202110981676.3A CN202110981676A CN113676108B CN 113676108 B CN113676108 B CN 113676108B CN 202110981676 A CN202110981676 A CN 202110981676A CN 113676108 B CN113676108 B CN 113676108B
- Authority
- CN
- China
- Prior art keywords
- flux linkage
- rotor
- rotor flux
- axis component
- observer
- 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
- 230000004907 flux Effects 0.000 title claims abstract description 127
- 238000000034 method Methods 0.000 title claims abstract description 101
- 230000006698 induction Effects 0.000 title claims abstract description 61
- 239000011159 matrix material Substances 0.000 claims description 9
- 238000011438 discrete method Methods 0.000 claims description 4
- 230000007704 transition Effects 0.000 claims description 4
- 239000004576 sand Substances 0.000 claims description 3
- 230000008569 process Effects 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 3
- 239000006185 dispersion Substances 0.000 description 2
- 238000001514 detection method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- 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
-
- 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
- H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
- H02P21/13—Observer control, e.g. using Luenberger observers or Kalman filters
-
- 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
- H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
- H02P21/14—Estimation or adaptation of machine parameters, e.g. flux, current or voltage
- H02P21/141—Flux estimation
-
- 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/12—Observer control, e.g. using Luenberger observers or Kalman filters
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Ac Motors In General (AREA)
Abstract
一种感应电机转子磁链的观测方法,涉及电机控制技术领域。本发明是为了解决现有对磁链观测器进行离散化时,离散化误差大,导致磁链观测精度低的问题。本发明所述的一种感应电机转子磁链的观测方法,采集感应电机定子电流,建立转子磁链观测器,利用改进欧拉法对转子磁链观测器进行离散化并获得观测值,所述观测值用于对感应电机进行控制,所述改进欧拉法为前向欧拉法与梯形离散法的结合离散化方法。本发明适用于感应电机控制中的转子磁链观测部分,以提高磁链观测精度。
Description
技术领域
本发明属于电机控制技术领域。
背景技术
感应电机因其结构简单,强耐用性,高可靠性及低成本等优点被广泛应用于变频调速系统。磁链是导电线圈或电流回路所链环的磁通量。由于感应电机高精度的控制过程离不开对磁链的控制,所以获取高精确度的磁链信息在感应电机控制中是至关重要的一环。
在实际的感应电机数字化驱动系统中,感应电机αβ坐标系下转子磁链观测电流模型(连续域)如下:
由于感应电机控制信号是离散的,因而要对磁链观测器进行离散化来近似拟合其连续域的观测电流模型。
具体的:将式(1)所示的转子磁链观测电流模型写成如下形式:
定义Ts为离散步长,并令t=(k+1)Ts,根据非齐次矩阵方程定理可将式(2)离散化为:
对式(3)进行求解可得:
同样可以得到:
若假设输入信号在迭代周期内不变,可将式(6)等效为:
为了在实际的数字控制系统中对磁链观测器进行离散化,现有技术中常采用基于前向欧拉法的离散转子磁链电流模型、基于后向欧拉法的离散转子磁链电流模型和基于双线性法的离散转子磁链电流模型。
1、基于前向欧拉法的离散转子磁链电流模型(最常采用)。
前向欧拉法结构简单,易于编程实现,是应用最为广泛的离散化方法。首先,假设有一阶常微分方程为:
dx/dt=f(x,y) (1-1)
式中,f(x,y)是关于x和y的函数,a和b是实数且a≤x≤b。根据前向欧拉法,可以将式(1-1)离散化为:
x(k+1)≈x(k)+Tsf[x(k),y(k)] (1-2)
其中,k=0,1,2...。
根据式(3)与式(1-2),可以得到基于前向欧拉法离散化的磁链电流模型观测器为:
从式(1-3)可以看出,前向欧拉法计算量小易于数字实现,但它存在离散误差随转速升高而增大的问题。
2、基于后向欧拉法的离散转子磁链电流模型。
根据后向欧拉法,可以将式(1-1)离散化为:
x(k+1)≈x(k)+Tsf[x(k+1),y(k+1)] (2-1)
根据式(2)与式(2-1),可以得到基于后向欧拉法离散化的磁链电流模型观测器为:
式(2-2)的完整形式可以表示为:
可以看出,基于后向欧拉法的离散过程中存在矩阵的求逆运算,所以该离散方法在每个迭代周期内都伴随着大量的乘除法运算,这显然加重了处理器的负担,各项系数的物理意义并不明确,实际调试困难,不利于实际应用。
3、基于双线性法的离散转子磁链电流模型。
根据双线性法,可以将式(1-1)离散化为:
根据式(1)与式(3-1),可以得到基于双线性法离散化的磁链电流模型观测器为:
式(3-2)的完整形式可以表示为:
由上式可见,基于双线性法的离散过程也存在矩阵的求逆运算,而且比后向欧拉法更加复杂,计算量更大,因此也更加不实用。
综上所述,现有对磁链观测器进行离散化时,后向欧拉法和双线性法计算量巨大,不利于实际应用。而前向欧拉法虽然计算量小、易于数字实现,但其离散化误差大,导致磁链观测精度低。
发明内容
本发明是为了解决现有对磁链观测器进行离散化时,离散化误差大,导致磁链观测精度低的问题,现提供一种感应电机转子磁链的观测方法。
一种感应电机转子磁链的观测方法,采集感应电机定子电流,建立转子磁链观测器,利用改进欧拉法对转子磁链观测器进行离散化,并获得离散化后的转子磁链观测器的观测值,上述观测值用于对感应电机进行控制,上述改进欧拉法为前向欧拉法与梯形离散法的结合离散化方法。
进一步的,上述转子磁链观测器为αβ坐标系下感应电机的转子磁链电流模型。
进一步的,上述转子磁链观测器的微分表达式如下:
进一步的,上述转子磁链观测器在αβ坐标系下的表达式为:
进一步的,上述利用改进欧拉法对转子磁链观测器进行离散化包括:
利用前向欧拉法对转子磁链观测器进行离散化,获得一次离散结果:
以及,基于一次离散结果,获得离散后的感应电机转子磁链观测器。
进一步的,上述基于一次离散结果,获得离散后的感应电机转子磁链观测器包括:
利用梯形离散法对一次离散结果再次进行离散化,获得离散后的感应电机转子磁链观测器:
进一步的,上述离散后的感应电机转子磁链观测器在αβ坐标系下的表达式为:
其中,和分别为kTs时刻转子磁链观测值的α轴分量和β轴分量,和分别为(k+1)Ts时刻转子磁链观测值的α轴分量和β轴分量,Kα(k)和Kβ(k)分别为kTs时刻转子磁链观测值α轴分量和β轴分量的微分,Kα(k+1)和Kβ(k+1)分别为(k+1)Ts时刻转子磁链观测值α轴分量和β轴分量的微分。
进一步的,上述Kα(k)和Kβ(k)的具体表达式如下:
其中,isα(k)和isβ(k)分别为kTs时刻定子电流α轴分量和β轴分量,ωr(k)为kTs时刻转子角速度。
进一步的,上述Kα(k+1)和Kβ(k+1)的具体表达式如下:
其中,isα(k+1)和isβ(k+1)分别为(k+1)Ts时刻定子电流α轴分量和β轴分量,ωr(k+1)为(k+1)Ts时刻转子角速度。
进一步的,上述离散后的感应电机转子磁链观测器的离散误差ε表达式如下:
其中,φ为离散后转子磁链观测器的状态转移矩阵。
本发明的有益效果为:
1、离散化复杂程度低,易于数字化实现;
2、迭代过程中运算次数较少,对处理器要求不高;
3、离散化精度很高,在感应电机的全速范围内的磁链观测几乎不存在误差;
4、离散化方程中各项物理意义明确,便于实际过程的调试。
综上所述,本发明所述的一种感应电机转子磁链的观测方法,采用改进欧拉法对转子磁链观测器进行离散化,相较于目前已存在的离散化方法具有数字化实现简单,离散化精度高的优点,能够取得离散复杂度与离散误差之间的平衡,使得本发明能够保证离散化精度,进一步的还能够提高感应电机的控制精度。本发明适用于感应电机控制中的转子磁链观测部分,以提高磁链观测精度。
附图说明
图1为四种离散化方法离散误差随电机运行频率变化的曲线图;
图2为基于前向欧拉法的感应电机转子磁链观测结果曲线图;
图3为基于改进欧拉法的感应电机转子磁链观测结果曲线图。
具体实施方式
为了兼顾离散误差与离散复杂程度,进而提高对感应电机的控制精度,本实施方式改进了转子磁链观测器的离散化方法,使得感应电机转子磁链的观测更加精确,具体如以下实施方式。
具体实施方式一:本实施方式所述的一种感应电机转子磁链的观测方法,
首先,采集感应电机定子电流,建立转子磁链观测器;然后,利用改进欧拉法对转子磁链观测器进行离散化,并根据离散后的转子磁链观测器获得观测值。改进欧拉法为前向欧拉法与梯形离散法的结合离散化方法。观测值用于对感应电机进行控制。
在电机从0速阶梯式加速至额定转速50Hz的条件下,基于前向欧拉法与改进欧拉法的感应电机转子磁链观测结果对比,如图2和3所示。由两图可以看到:改进欧拉离散化方法在全速范围内都能够实现对转子磁链的精确观测,几乎没有误差。而传统的前向欧拉离散化方法所观测的转子磁链误差随转速升高而增大,最大观测误差已超过100%,电机控制系统无法正常运行。
本实施方式在实际应用时,能够使用在动车故障检测领域里AGV小车的驱动控制中。具体的,采集AGV小车中感应电机的定子电流,建立转子磁链观测器;利用改进欧拉法对转子磁链观测器进行离散化,并根据离散化后的转子磁链观测器获得精确的磁链观测值;之后利用该磁链观测值获得磁链角,根据变换电机电压和电流的坐标,用变换坐标后的电压和电流对AGV小车的电机进行精确控制。
具体实施方式二:本实施方式是对具体实施方式一所述的一种感应电机转子磁链的观测方法作进一步说明,本实施方式中,转子磁链观测器为αβ坐标系下感应电机的转子磁链电流模型。
具体实施方式三:本实施方式是对具体实施方式一所述的一种感应电机转子磁链的观测方法作进一步说明,本实施方式中,转子磁链观测器的微分表达式如下:
具体实施方式四:本实施方式是对具体实施方式三所述的一种感应电机转子磁链的观测方法作进一步说明,本实施方式中,转子磁链观测器在αβ坐标系下的表达式为:
具体实施方式五:本实施方式是对具体实施方式三所述的一种感应电机转子磁链的观测方法作进一步说明,本实施方式中,利用改进欧拉法对转子磁链观测器进行离散化包括:
利用前向欧拉法对转子磁链观测器进行离散化,获得一次离散结果:
以及,基于一次离散结果,获得离散后的感应电机转子磁链观测器。
具体实施方式六:本实施方式是对具体实施方式五所述的一种感应电机转子磁链的观测方法作进一步说明,本实施方式中,基于一次离散结果,获得离散后的感应电机转子磁链观测器包括:
利用梯形离散法对一次离散结果再次进行离散化,获得离散后的感应电机转子磁链观测器:
具体实施方式七:本实施方式是对具体实施方式六所述的一种感应电机转子磁链的观测方法作进一步说明,本实施方式中,离散后的感应电机转子磁链观测器在αβ坐标系下的表达式为:
其中,和分别为kTs时刻转子磁链观测值的α轴分量和β轴分量,和分别为(k+1)Ts时刻转子磁链观测值的α轴分量和β轴分量,Kα(k)和Kβ(k)分别为kTs时刻转子磁链观测值α轴分量和β轴分量的微分,Kα(k+1)和Kβ(k+1)分别为(k+1)Ts时刻转子磁链观测值α轴分量和β轴分量的微分。
具体实施方式八:本实施方式是对具体实施方式七所述的一种感应电机转子磁链的观测方法作进一步说明,本实施方式中,Kα(k)和Kβ(k)的具体表达式如下:
其中,isα(k)和isβ(k)分别为kTs时刻定子电流α轴分量和β轴分量,ωr(k)为kTs时刻转子角速度。
具体实施方式九:本实施方式是对具体实施方式七或八所述的一种感应电机转子磁链的观测方法作进一步说明,本实施方式中,Kα(k+1)和Kβ(k+1)的具体表达式如下:
其中,isα(k+1)和isβ(k+1)分别为(k+1)Ts时刻定子电流α轴分量和β轴分量,ωr(k+1)为(k+1)Ts时刻转子角速度。
具体实施方式十:本实施方式是对具体实施方式六所述的一种感应电机转子磁链的观测方法作进一步说明,本实施方式中,离散后的感应电机转子磁链观测器的离散误差ε表达式如下:
其中,φ为离散后转子磁链观测器的状态转移矩阵。
取离散步长分别为1/3000秒与1/6000秒,电机运行频率为0-200Hz。目前已存在的三种离散化方法与本实施方式中的离散化方法(改进欧拉法)的离散化误差随电机转速升高而变化的曲线图,如附图1所示能够看到:
四种方法的离散化误差均随电机运行频率的增加而增大,且随离散化步长的增大而增大。并且当离散步长相同时,双线性法和改进欧拉法的离散化误差要明显小于前向欧拉法和后向欧拉法。当离散化步长为1/3000秒、电机运行频率在200Hz时,前向欧拉法与后向欧拉法的离散化误差达到了7%以上;而改进欧拉法和双线性法的离散化误差始终保持在2%以内。因此,本实施方式中改进的欧拉法在不过多增加离散运算的同时有效地减小了离散误差,在离散误差与离散复杂度之间取得平衡,适用于感应电机数字控制系统。
Claims (7)
1.一种感应电机转子磁链的观测方法,其特征在于,
采集感应电机定子电流,建立转子磁链观测器,
利用改进欧拉法对转子磁链观测器进行离散化,并获得离散化后的转子磁链观测器的观测值,
所述观测值用于对感应电机进行控制,
所述改进欧拉法为前向欧拉法与梯形离散法的结合离散化方法;
所述转子磁链观测器的微分表达式如下:
所述利用改进欧拉法对转子磁链观测器进行离散化包括:
利用前向欧拉法对转子磁链观测器进行离散化,获得一次离散结果:
以及,基于一次离散结果,获得离散后的感应电机转子磁链观测器;
所述基于一次离散结果,获得离散后的感应电机转子磁链观测器包括:
利用梯形离散法对一次离散结果再次进行离散化,获得离散后的感应电机转子磁链观测器:
2.根据权利要求1所述的一种感应电机转子磁链的观测方法,其特征在于,所述转子磁链观测器为αβ坐标系下感应电机的转子磁链电流模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110981676.3A CN113676108B (zh) | 2021-08-25 | 2021-08-25 | 一种感应电机转子磁链的观测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110981676.3A CN113676108B (zh) | 2021-08-25 | 2021-08-25 | 一种感应电机转子磁链的观测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113676108A CN113676108A (zh) | 2021-11-19 |
CN113676108B true CN113676108B (zh) | 2022-03-08 |
Family
ID=78546114
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110981676.3A Active CN113676108B (zh) | 2021-08-25 | 2021-08-25 | 一种感应电机转子磁链的观测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113676108B (zh) |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9729093B2 (en) * | 2015-12-15 | 2017-08-08 | Whirlpool Corporation | Observer based sensorless control for U-shape single phase synchronous permanent magnet motors |
CN109274303B (zh) * | 2018-11-29 | 2021-06-08 | 合肥工业大学 | 一种基于扩展控制集的异步电机模型预测控制方法 |
CN109818541B (zh) * | 2019-03-13 | 2020-10-02 | 东南大学 | 一种用于磁链观测的记忆电机绕组复用控制方法及系统 |
CN110445444A (zh) * | 2019-07-12 | 2019-11-12 | 西安理工大学 | 一种异步电机改进的模型预测控制方法 |
-
2021
- 2021-08-25 CN CN202110981676.3A patent/CN113676108B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN113676108A (zh) | 2021-11-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108319148B (zh) | 一种控制力矩陀螺框架伺服系统低转速高精度控制方法 | |
CN104639003B (zh) | 一种交流伺服系统的转动惯量辨识方法 | |
CN103592946B (zh) | 一种基于视加速度测量的主动段程序自适应纵向制导方法 | |
CN106385211B (zh) | 一种步进电机负载转矩估计方法 | |
CN110401391B (zh) | 异步电动机随机系统模糊自适应动态面控制方法 | |
CN104201967B (zh) | 一种采用自抗扰控制技术的网络化永磁同步电机时延补偿和控制方法 | |
CN104993765A (zh) | 一种无刷直流电机的转速估计方法 | |
CN105353789A (zh) | 连续振动信号时间历程复现控制方法 | |
CN112422002B (zh) | 一种鲁棒性的永磁同步电机单电流传感器预测控制方法 | |
CN109586645A (zh) | 一种永磁同步电机惯量识别方法及设备 | |
CN113676108B (zh) | 一种感应电机转子磁链的观测方法 | |
CN114035588B (zh) | 一种移动机器人轨迹跟踪事件触发控制方法 | |
CN107681937A (zh) | 一种基于神经网络的超高速永磁同步电机速度观测方法 | |
CN113283090B (zh) | 基于蜂群算法的飞机舵机电动伺服系统摩擦补偿方法 | |
CN112468034B (zh) | 永磁同步电机弱磁区效率最优控制电流轨迹搜索方法及在线控制方法 | |
CN107707169B (zh) | 一种直线感应电机无速度传感器控制系统及方法 | |
CN112083349A (zh) | 一种永磁同步电机定子绕组匝间短路故障诊断方法 | |
CN113890442B (zh) | 一种粒子群滑模观测器风力发电机转子电流故障检测方法 | |
CN108073087B (zh) | 一种多学科综合仿真平台环境下的机电伺服系统控制方法 | |
CN112468032B (zh) | 永磁同步电机的全速域效率map图生成方法 | |
CN104038132B (zh) | 具有时变测量延迟输出和噪声的伺服电机的状态观测方法 | |
CN104022707A (zh) | 基于新型转子磁链观测器的异步电机速度控制装置及系统 | |
CN114826071A (zh) | 基于霍尔传感器的永磁同步电机转子位置和速度估计方法 | |
Sun et al. | A linearization control scheme for permanent magnet linear synchronous motors | |
CN114172429B (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 |