CN113794411B - 航空柱塞泵用内嵌式永磁同步电机的多重抗扰控制方法 - Google Patents

航空柱塞泵用内嵌式永磁同步电机的多重抗扰控制方法 Download PDF

Info

Publication number
CN113794411B
CN113794411B CN202111013500.5A CN202111013500A CN113794411B CN 113794411 B CN113794411 B CN 113794411B CN 202111013500 A CN202111013500 A CN 202111013500A CN 113794411 B CN113794411 B CN 113794411B
Authority
CN
China
Prior art keywords
load
plunger pump
aviation
aviation plunger
permanent magnet
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
Application number
CN202111013500.5A
Other languages
English (en)
Other versions
CN113794411A (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.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical University
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 Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN202111013500.5A priority Critical patent/CN113794411B/zh
Publication of CN113794411A publication Critical patent/CN113794411A/zh
Application granted granted Critical
Publication of CN113794411B publication Critical patent/CN113794411B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/05—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation specially adapted for damping motor oscillations, e.g. for reducing hunting
    • 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/20—Estimation of torque
    • 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
    • H02P25/00—Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details
    • H02P25/02—Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details characterised by the kind of motor
    • H02P25/022—Synchronous motors
    • H02P25/024—Synchronous motors controlled by supply frequency
    • H02P25/026—Synchronous motors controlled by supply frequency thereby detecting the rotor position
    • 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
    • H02P27/00—Arrangements or methods for the control of AC motors characterised by the kind of supply voltage
    • H02P27/04—Arrangements or methods for the control of AC motors characterised by the kind of supply voltage using variable-frequency supply voltage, e.g. inverter or converter supply voltage
    • H02P27/06—Arrangements or methods for the control of AC motors characterised by the kind of supply voltage using variable-frequency supply voltage, e.g. inverter or converter supply voltage using dc to ac converters or inverters
    • H02P27/08—Arrangements or methods for the control of AC motors characterised by the kind of supply voltage using variable-frequency supply voltage, e.g. inverter or converter supply voltage using dc to ac converters or inverters with pulse width modulation
    • 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
    • H02P2207/00—Indexing scheme relating to controlling arrangements characterised by the type of motor
    • H02P2207/05—Synchronous machines, e.g. with permanent magnets or DC excitation

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Control Of Positive-Displacement Pumps (AREA)

Abstract

本发明涉及一种航空柱塞泵用内嵌式永磁同步电机的多重抗扰控制方法，针对航空柱塞泵负载特性建立了离散域模型，并提出了泵用内嵌式永磁同步电机的高性能转速控制方法。本发明设计的泵电机转速控制器由自适应准比例积分谐振转速控制器和Luenberger动态带宽负载转矩观测器组成，能够分别解决航空柱塞泵周期脉动载荷造成的驱动电机稳态转速高频波动问题和航空柱塞泵非周期冲击载荷造成的驱动电机转速跟踪误差突增问题，具有系统高可靠性、高稳定性、强抗扰性等优点。

Description

航空柱塞泵用内嵌式永磁同步电机的多重抗扰控制方法
技术领域
本发明属于电机驱动控制技术领域，涉及一种航空柱塞泵用内嵌式永磁同步电机的多重抗扰控制方法。
背景技术
航空柱塞泵电机系统背景：为了使航空发动机运转顺滑，及时传导出摩擦产生的高温，需要航空滑油系统快速、不间断地将航空润滑油输送到发动机转动部件，保证发动机长时间无故障工作。由于柱塞泵具有排量易调节、高功率、高可靠性等优点，滑油供油系统中通常采用柱塞泵作为主供压泵源。在多电飞机“功率电传”概念的影响下，使用高功率密度、高可靠性的永磁同步电机作为柱塞泵的动力源是目前的发展趋势。这种由永磁同步电机、柱塞泵、数字控制器、高性能控制策略构成的分布式电动泵系统，相比于使用集中式传统的机械泵系统，避免了复杂、冗长的联轴器和运动转换机构，有节省建造成本和使用空间的优点。然而，考虑到柱塞泵工作在高压高速条件下，任何扰动都有可能破坏系统的稳定性，从而造成严重的航空事故，因此，为了维持航空柱塞泵用内嵌式永磁同步电机的高性能、高可靠控制特性，需针对航空柱塞泵负载特性设计高性能控制策略。
柱塞泵电机负载特性与控制技术现状：柱塞泵的柱塞在缸体柱塞孔中往复运动实现吸油和排油。油液的瞬时流量与柱塞的数量和缸体的转速相关。缸体转速与时间的乘积表示缸体转角，当角度变化时柱塞的轴向位移距离和移动速度将相应变化，因此泵在一个工作周期中的瞬时输出流量将随转角变化产生周期性脉动现象。周期性的流量脉动反映在电机转轴上为高频的负载力矩。同时，考虑到柱塞泵斜盘倾角造成的排量突变以及外部干扰的影响，将柱塞泵的负载特性归结为周期脉动载荷和非周期冲击载荷。目前，柱塞泵周期性流量脉动和压力脉动的主要解决方法是通过改变柱塞泵机械结构或增加机械补偿装置实现，这种方法的实现成本高，且将降低柱塞泵功重比，不利于航空滑油系统级设计。为了实时补偿泵电机的非周期冲击载荷，提高控制系统抗扰性能，通常采用滑模控制器、模糊控制器等非线性控制器。其中，滑模控制器具有快速响应性能，而当控制参数到达滑模面后将产生围绕滑模面的抖动，不利于稳定控制。此外，模糊控制器需要依靠经验设计模糊规则，且稳定性无法保证。
航空柱塞泵电机系统控制策略存在以下问题：1)航空柱塞泵负载特性建模问题，2)航空柱塞泵周期脉动载荷扰动补偿问题，3)航空柱塞泵非周期冲击载荷扰动补偿问题。
发明内容
要解决的技术问题
为了避免现有技术的不足之处，本发明提出一种航空柱塞泵用内嵌式永磁同步电机的多重抗扰控制方法，要解决：1)航空柱塞泵负载特性建模问题，2)航空柱塞泵周期脉动载荷扰动补偿问题，3)航空柱塞泵非周期冲击载荷扰动补偿问题。
技术方案
一种航空柱塞泵用内嵌式永磁同步电机的多重抗扰控制方法，其特征在于步骤如下：
步骤1：将永磁同步电机的转子角速度ω和转速角度θ，采集电机的定子三相电流ia、ib、ic，进行坐标变换获得旋转坐标系下的电流分量id、iq，并计算电磁转矩值Te；
步骤2：由于航空柱塞泵的瞬时流量具有周期脉动特性，同时考虑到瞬时冲击性载荷的影响，对航空柱塞泵的负载特性进行建模为：
Figure BDA0003239804310000021
其中，Ttotal(k)为所述航空柱塞泵负载特性；
Figure BDA0003239804310000022
为航空柱塞泵非周期性冲击载荷；
Figure BDA0003239804310000023
为航空柱塞泵周期脉动载荷；
步骤3、设计自适应准比例积分谐振控制器：依据航空柱塞泵的机械参数计算得到周期脉动载荷的实时谐振频率，再使用所述实时谐振频率设计自适应准比例积分谐振控制器，针对所述航空柱塞泵周期脉动载荷
Figure BDA0003239804310000031
设计自适应准比例积分谐振控制器：
(3.1)依据航空柱塞泵的机械参数计算获得航空柱塞泵周期脉动载荷实时谐振频率：
Figure BDA0003239804310000032
其中，z为航空柱塞泵柱塞数量，ω为航空柱塞泵旋转角速度，f1、f2分别为偶数柱塞和奇数柱塞的周期脉动载荷实时谐振频率；
(3.2)将步骤1获得的转子角速度ω和角速度给定值ω*做差，得到角速度反馈误差Δω，在s域内设计自适应准比例积分谐振控制器，其中，控制器输入为角速度反馈误差Δω，输出为定子电流给定值
Figure BDA0003239804310000033
表达式如下：
Figure BDA0003239804310000034
式中，谐振角频率ω0可根据所述实时谐振频率计算得到，ω0＝2πfi,i＝1or 2，Kp为比例参数，Ki为积分参数，Kr为谐振参数，ωb为谐振带宽，
Figure BDA0003239804310000035
为定子电流给定值；
通过控制定子电流给定值
Figure BDA0003239804310000036
实现对所述航空柱塞泵周期脉动载荷
Figure BDA0003239804310000037
的控制，在稳态工况下，
Figure BDA0003239804310000038
与
Figure BDA0003239804310000039
的关系式为：
Figure BDA00032398043100000310
式中
Figure BDA00032398043100000311
其中p为电机极对数，
Figure BDA00032398043100000312
为永磁体磁链；
步骤4：依据步骤1获得的电磁转矩Te和永磁同步电机的运动方程设计观测器模型，使用转速观测误差驱动增益变化，从而获得冲击载荷的观测值和电机角速度的观测值；
针对所述航空柱塞泵非周期冲击载荷
Figure BDA0003239804310000041
设计Luenberger动态带宽负载转矩观测器对其进行观测，观测器表达式为：
Figure BDA0003239804310000042
Figure BDA0003239804310000043
式中，
Figure BDA0003239804310000044
分别为角速度观测值和负载转矩观测值，J、Bm分别为永磁同步电机的惯性系数和摩擦系数，Te为电磁转矩，k1、k2为观测器增益系数；
非周期冲击载荷的观测值
Figure BDA0003239804310000045
即为所述航空柱塞泵非周期性冲击载荷
Figure BDA0003239804310000046
Figure BDA0003239804310000047
步骤5：依据步骤3和步骤4获得的
Figure BDA0003239804310000048
和
Figure BDA0003239804310000049
经过计算获得修正的定子电流给定值
Figure BDA00032398043100000410
Figure BDA00032398043100000411
再将
Figure BDA00032398043100000412
通过MTPA运算获得dq轴定子电流给定值
Figure BDA00032398043100000413
和
Figure BDA00032398043100000414
其中，MTPA运算公式如下：
Figure BDA00032398043100000415
式中，
Figure BDA00032398043100000416
为永磁体磁链，Lq、Ld分别表示q轴电感和d轴电感，sgn表示符号函数；
步骤6：将步骤5获得的dq轴定子电流给定值
Figure BDA00032398043100000417
和
Figure BDA00032398043100000418
经过电流控制器输出dq轴电压给定值
Figure BDA00032398043100000419
和
Figure BDA00032398043100000420
再经过坐标变换得到静止旋转坐标系下的电压给定值
Figure BDA00032398043100000421
和
Figure BDA00032398043100000422
并将
Figure BDA00032398043100000423
和
Figure BDA00032398043100000424
输入SVPWM进行脉宽调制，输出PWM波触发逆变器，实现航空柱塞泵用内嵌式永磁同步电机的多重抗扰控制。
所述步骤2的航空柱塞泵周期性脉动载荷由瞬时脉动流量对时间求导数获得，离散域的表达式为：
Figure BDA0003239804310000051
式中，K为比例系数，
Figure BDA0003239804310000052
为单周期脉动流量变化值，其表达式为
Figure BDA0003239804310000053
Figure BDA0003239804310000054
和
Figure BDA0003239804310000055
分别表示k时刻和k-1时刻的脉动流量，Δt为周期间隔时间。
所述步骤2的非周期性冲击载荷的离散域模型为：航空柱塞泵非周期性冲击载荷由外部冲击主载荷和内部附加载荷构成，其中内部附加载荷包括滑动摩擦载荷、粘性摩擦载荷和滚动摩擦载荷，表达式如下：
Figure BDA0003239804310000056
式中，d(k)为所述航空柱塞泵外部冲击主载荷，ΔTt(k)为所述内部附加载荷，表达式为ΔTt(k)＝Tf(k)+Tv(k)+Ts(k)，其中Tf(k)为滑动摩擦载荷，Tv(k)为粘性摩擦载荷，Ts(k)为滚动摩擦载荷。
所述步骤4的Luenberger动态带宽负载转矩观测器的增益系数k1、k2获取方式如下：
1)建立观测器特征方程
Figure BDA0003239804310000057
依据特征方程获得初始增益系数k′1和k′2的表达式：
k′1＝Jω‘
Figure BDA0003239804310000058
式中，ω’为设定的观测器带宽，Bm为电机摩擦系数，J为电机惯性系数；
2)依据所述初始增益系数k′1和k′2，建立基于双曲正切函数的增益系数k1、k2表达式为：
Figure BDA0003239804310000059
Figure BDA00032398043100000510
式中，β1、β2、c1、c2均为比例系数，
Figure BDA00032398043100000511
为观测器的转速估计误差。
有益效果
本发明提出的一种航空柱塞泵用内嵌式永磁同步电机的多重抗扰控制方法，针对航空柱塞泵负载特性建立了离散域模型，并提出了泵用内嵌式永磁同步电机的高性能转速控制方法。本发明设计的泵电机转速控制器由自适应准比例积分谐振转速控制器和Luenberger动态带宽负载转矩观测器组成，能够分别解决航空柱塞泵周期脉动载荷造成的驱动电机稳态转速高频波动问题和航空柱塞泵非周期冲击载荷造成的驱动电机转速跟踪误差突增问题，具有系统高可靠性、高稳定性、强抗扰性等优点。
本发明主要是根据航空柱塞泵结构和运动特征，将负载特性划分为周期性脉动载荷和非周期性冲击载荷，从而针对性设计高性能控制器进行补偿，实现航空柱塞泵电机系统的强稳定性、强抗扰、高可靠性控制。
与现有技术相比，该方法的有益效果在于：
(1)结合航空柱塞泵的特征，将泵用驱动电机的负载转矩归结为周期脉动载荷和非周期冲击载荷，便于针对不同负载特性设计高性能控制器。
(2)针对周期脉动载荷，使用自适应方法实时获取脉动载荷的谐振频率，采用准比例积分谐振控制器降低了周期性脉动载荷对泵电机造成的转速波动，提高了电动泵系统的运行可靠性。
(3)针对非周期冲击载荷，设计了基于Luenberger动态带宽负载转矩观测器，实时补偿泵电机的非周期冲击载荷，提高了泵电机对冲击载荷的转速跟踪性能和抗扰性能。
附图说明
图1：航空柱塞泵用内嵌式永磁同步电机控制架构
图2：航空柱塞泵电机周期性脉动载荷波形
图3：准比例谐振控制器Bode图
图4：驱动电机稳态转速对比
图5：基于双曲正切函数的动态带宽负载转矩观测器的增益变化
图6：基于Luenberger动态带宽负载转矩观测器的估计负载转矩波形
图7：驱动电机起动和加载工况转速波形对比
图8：驱动电机突加干扰时转速波形对比
具体实施方式
现结合实施例、附图对本发明作进一步描述：
航空柱塞泵电机系统控制策略存在以下问题：1)航空柱塞泵负载特性建模问题，2)航空柱塞泵周期脉动载荷扰动补偿问题，3)航空柱塞泵非周期冲击载荷扰动补偿问题。
航空柱塞泵用电机的多重抗扰控制方法，其技术特征包括以下几方面：
(1)航空柱塞泵负载特性建模。基本原理如下：
若轴向柱塞泵的柱塞数为z，柱塞直径为d，柱塞孔的分布圆直径为D，斜盘轴线与缸体轴线间的夹角为β，当缸体转动一周时，泵的排量可表示为
Figure BDA0003239804310000071
由式可知调节泵的排量只能通过改变斜盘倾角实现。在定排量柱塞泵中斜盘倾角不发生变化，因此在以下推导过程中皆假设β为恒值。
当缸体转角
Figure BDA0003239804310000072
时，柱塞的轴向位移S为
Figure BDA0003239804310000073
将上式对时间变量t求导数，得柱塞的瞬时移动速度V
Figure BDA0003239804310000074
故单个柱塞的瞬时流量为
Figure BDA0003239804310000075
由(4)可知，单个柱塞的瞬时流量依照正弦规律变化。当液压泵有z个柱塞，柱塞间的夹角即为
Figure BDA0003239804310000076
则所有压油区的各柱塞的瞬时理论流量为
Figure BDA0003239804310000081
式中，q1、q2、qm为压油区各柱塞的瞬时理论流量；m为处于压油过程的柱塞数。
柱塞泵的总瞬时流量是压油区的所有柱塞瞬时流量的总和。所以总瞬时流量为
Figure BDA0003239804310000082
假设正弦脉动函数
Figure BDA0003239804310000083
经过化简可知，当柱塞数为偶数时，柱塞泵的瞬时脉动流量为
Figure BDA0003239804310000084
当柱塞数为奇数时，柱塞泵的瞬时脉动流量为
Figure BDA0003239804310000085
由式(6)可知，瞬时流量与角速度成正比。通过对式(7)(8)求导，可得
Figure BDA0003239804310000086
由式(9)可知，瞬时脉动流量的导数与柱塞泵缸体角加速度a′成正比。根据转动定律可知，在柱塞泵缸体角加速度a′的影响下将产生周期性脉动载荷，且角加速度与脉动载荷成比例关系，因此可建立电机周期脉动载荷模型
Figure BDA0003239804310000087
其中K为比例系数，
Figure BDA0003239804310000088
为相邻时间间隔的流量差值，Δt为间隔时间。
柱塞泵非周期性冲击载荷具有时间和幅值不确定的特点，同时考虑到柱塞泵内部机械和液体间的摩擦力，将航空柱塞泵非周期冲击载荷归结为外部冲击主载荷和内部附加载荷。假设冲击载荷作用时间为第k采样周期，则基于离散域的泵电机非周期冲击载荷表达式为
Figure BDA0003239804310000091
式中，d(k)为所述航空柱塞泵外部冲击主载荷，ΔTt(k)为所述外部冲击影响下的内部附加载荷，表达式为ΔTt(k)＝Tf(k)+Tv(k)+Ts(k)，其中Tf(k)为滑动摩擦载荷，Tv(k)为粘性摩擦载荷，Ts(k)为滚动摩擦载荷。
通过使用欧拉方程将周期性脉动载荷离散化，因此离散域内的航空柱塞泵电机总的负载转矩表达式为
Figure BDA0003239804310000092
式中单周期脉动流量变化值
Figure BDA0003239804310000093
Δt为间隔时间。
(2)由航空柱塞泵负载特性建模可知，周期性脉动载荷的变化频率与电机机械角速度成正比，即电机角速度越快，脉动载荷变化频率越快。因此，电机转子侧被施加周期性脉动载荷后将产生周期性转速波动，易造成驱动电机的控制系统不稳定。为了解决以上问题，在此针对柱塞泵周期性脉动载荷设计转速环的自适应准比例积分谐振转速控制器。
基于自适应准比例积分谐振控制器的航空柱塞泵用电机周期性脉动载荷补偿方法，基本原理如下：
当柱塞泵的柱塞数为z且角速度为ω时，柱塞泵周期性脉动载荷谐振频率为
Figure BDA0003239804310000094
考虑到理想比例谐振控制器仅对单一谐振角频率起作用，而在实际应用中，信号采集频率将受到噪声、电路延迟等因素发生偏移，因此围绕主谐振频率f，设计自适应准比例谐振控制器，其传递函数为
Figure BDA0003239804310000095
式中ωb为谐振带宽，能够对谐振频率ω0周围ωb频率范围内的信号起作用；Kp为比例系数；Kr为谐振系数。其中，参考柱塞泵周期性脉动载荷特性，谐振控制器的谐振频率取ω0＝2πfi，i＝1 or 2。
为了实现转速无静差控制，将设计的自适应准比例谐振控制器与比例积分控制器组合为自适应准比例积分谐振控制器，其传递函数为
Figure BDA0003239804310000101
其中Kp和Ki分别为控制器的比例和积分系数。
由式(15)建立自适应准比例积分谐振控制器的表达式为
Figure BDA0003239804310000102
其中Δω为给定转速和反馈转速的误差，
Figure BDA0003239804310000103
为定子电流给定值。自适应准比例积分谐振控制器通过控制定子电流给定值削弱周期性脉动载荷
Figure BDA0003239804310000104
对转速的影响，在稳态工况下，电流给定值与周期性脉动载荷的关系为
Figure BDA0003239804310000105
Kt为参数项。
(3)基于Luenberger动态带宽负载转矩观测器的航空柱塞泵用电机非周期冲击载荷扰动补偿方法，其原理如下：
由式(11)可知，对于驱动电机而言，柱塞泵的非周期冲击载荷可看作阶跃信号。为了实现小冲击载荷强跟踪性能，大冲击载荷强抗扰性能，设计了Luenberger动态带宽负载转矩观测器可实时估计出负载转矩并进行前馈补偿，提高电动泵系统控制性能。
电机运动方程如下：
Figure BDA0003239804310000106
式中，J为转动惯量；ω为电机机械角速度；Bm为摩擦系数；TL为负载转矩；Te为电磁转矩。电机角速度ω与机械位置角度θ的关系为
Figure BDA0003239804310000107
根据运动方程，将电机负载转矩、电机机械角速度作为状态变量，速度作为输出变量得到以下状态方程
Figure BDA0003239804310000111
基于式(17)，建立负载转矩观测器模型如下
Figure BDA0003239804310000112
其中负载转矩观测值与非周期冲击载荷的关系为
Figure BDA0003239804310000113
由观测值
Figure BDA0003239804310000114
即可通过公式推导出定子电流给定值的补偿值。
通过将式(18)(19)做差可得观测器误差状态方程为
Figure BDA0003239804310000115
可得观测器的特征方程为
Figure BDA0003239804310000116
根据式(20)设定初始增益为k′1＝Jω‘，
Figure BDA0003239804310000117
保证观测器稳定性。由于冲击载荷的幅值和施加时间未知，基于初始增益系数，设计的动态带宽观测器的增益分别为
Figure BDA0003239804310000118
b,c为系数，k′1和k′2分别为初始增益。这样设计的好处是能够保证大冲击载荷的条件下快速减小转速误差，实现大负载转矩强抗扰；小冲击载荷的条件下避免转速超调，保持转速跟踪，实现转速强追踪控制性能。
现结合附图和实施例对本发明作进一步描述：
由于周期性脉动载荷只与柱塞数和流量脉动频率相关，因此本发明实施例以柱塞数为10的航空柱塞泵为研究对象，驱动用内嵌式永磁同步电机具体参数如表1所示，控制系统框架如图1所示。
表1驱动用内嵌式永磁同步电机参数
Figure BDA0003239804310000121
实施例包含的具体步骤如下：
1.根据式(10)，建立航空柱塞泵离散化的周期脉动载荷模型如下式所示
Figure BDA0003239804310000122
式中
Figure BDA0003239804310000123
且
Figure BDA0003239804310000124
Δt＝t(k)-t(k-1)。
若驱动电机的给定转速为分段线性函数，如下式所示，则周期性脉动载荷变化波形如图2所示。
Figure BDA0003239804310000125
2.由式(13)可知柱塞泵在第k采样周期内的谐振频率如下
Figure BDA0003239804310000126
基于谐振频率设计的准比例谐振控制器如下式
Figure BDA0003239804310000127
式中ω0＝2πf1。假设准比例谐振控制器参数为Kp＝5，Kr＝30，ωb＝50。该控制器的Bode图如图3所示。由图可知设计的准比例谐振控制器在谐振频率处的幅值较大，对谐波抑制效果较好，且系统相位裕度为-45deg，满足工业要求的系统稳定相位裕度。
将设计的准比例谐振控制器和比例积分控制器相加，获得自适应准比例积分谐振控制器。当给定转速为3000rpm时，稳态工况下的比例积分控制器和自适应准比例积分谐振控制器的转速控制波形如图4所示。由图可知，比例积分控制器在稳态时受到周期脉动载荷的影响转速波动为1rpm，而自适应准比例积分谐振控制器的转速波动只有0.4rpm，能够有效减小转速周期性高频波动。
3.根据式(18)设计基于Luenberger动态带宽负载转矩观测器，模型如下
Figure BDA0003239804310000131
电机的电磁转矩和机械角速度为输入值，负载转矩估计值为输出值，对内环控制器输入进行前馈补偿，以提高速度环的抗扰性能。其中，基于双曲正切函数构造动态带宽负载转矩观测器的增益系数，变化趋势如图5所示，可知，增益随转速误差变大而增大，有助于快速补偿高幅值非周期冲击载荷。
假设在9s处突加16Nm负载扰动，负载转矩观测器的估计负载转矩值如图6所示。可知该负载转矩观测器能够快速估计出突加非周期冲击载荷，准确度高。起动和加载工况下的转速波形如图7所示。由图可知，起动时自适应准比例积分谐振控制器的起动超调更小，更有利于柱塞泵的高精度转速跟踪控制。图8为非周期冲击载荷影响下的电机转速波形。由图可知，在突加16Nm负载时，自适应准比例积分谐振控制器和Luenberger动态带宽负载转矩观测器的组合控制具有更小的转速跌落幅值，仅为5rpm，且转速回调速度更快；只有自适应准比例积分谐振控制器的作用时转速跌落为7rpm，而转速回调时间较长；比例积分控制器的转速跌落幅度更大，为14rpm。由此可知，驱动用永磁同步电机的自适应准比例积分谐振控制器和Luenberger动态带宽负载转矩观测器的组合控制对航空柱塞泵的周期脉动载荷和非周期冲击载荷的补偿效果最好，具有强稳定性和强抗扰控制性能。

Claims (4)

1.一种航空柱塞泵用内嵌式永磁同步电机的多重抗扰控制方法，其特征在于步骤如下：
步骤1：将永磁同步电机的转子角速度ω和转速角度θ，采集电机的定子三相电流ia、ib、ic，进行坐标变换获得旋转坐标系下的电流分量id、iq，并计算电磁转矩值Te；
步骤2：对航空柱塞泵的负载特性进行建模为：
Figure FDA00032398043000000110
其中，Ttotal(k)为所述航空柱塞泵负载特性；
Figure FDA0003239804300000019
为航空柱塞泵非周期性冲击载荷；
Figure FDA00032398043000000111
为航空柱塞泵周期脉动载荷；
步骤3、设计自适应准比例积分谐振控制器：针对所述航空柱塞泵周期脉动载荷
Figure FDA00032398043000000112
设计自适应准比例积分谐振控制器：
(3.1)依据航空柱塞泵的机械参数计算获得航空柱塞泵周期脉动载荷实时谐振频率：
Figure FDA0003239804300000011
其中，z为航空柱塞泵柱塞数量，ω为航空柱塞泵旋转角速度，f1、f2分别为偶数柱塞和奇数柱塞的周期脉动载荷实时谐振频率；
(3.2)将步骤1获得的转子角速度ω和角速度给定值ω*做差，得到角速度反馈误差Δω，在s域内设计自适应准比例积分谐振控制器，其中，控制器输入为角速度反馈误差Δω，输出为定子电流给定值
Figure FDA0003239804300000012
表达式如下：
Figure FDA0003239804300000013
式中，谐振角频率ω0可根据所述实时谐振频率计算得到，ω0＝2πfi,i＝1 or 2，Kp为比例参数，Ki为积分参数，Kr为谐振参数，ωb为谐振带宽，
Figure FDA0003239804300000014
为定子电流给定值；
通过控制定子电流给定值
Figure FDA0003239804300000015
实现对所述航空柱塞泵周期脉动载荷
Figure FDA0003239804300000016
的控制，在稳态工况下，
Figure FDA0003239804300000017
与
Figure FDA0003239804300000018
的关系式为：
Figure FDA0003239804300000021
式中
Figure FDA0003239804300000022
其中p为电机极对数，
Figure FDA0003239804300000023
为永磁体磁链；
步骤4：针对所述航空柱塞泵非周期冲击载荷
Figure FDA0003239804300000024
设计Luenberger动态带宽负载转矩观测器对其进行观测，观测器表达式为：
Figure FDA0003239804300000025
Figure FDA0003239804300000026
式中，
Figure FDA0003239804300000027
分别为角速度观测值和负载转矩观测值，J、Bm分别为永磁同步电机的惯性系数和摩擦系数，Te为电磁转矩，k1、k2为观测器增益系数；
非周期冲击载荷的观测值
Figure FDA0003239804300000028
即为所述航空柱塞泵非周期性冲击载荷
Figure FDA0003239804300000029
Figure FDA00032398043000000210
步骤5：依据步骤3和步骤4获得的
Figure FDA00032398043000000211
和
Figure FDA00032398043000000212
经过计算获得修正的定子电流给定值
Figure FDA00032398043000000213
Figure FDA00032398043000000214
再将
Figure FDA00032398043000000215
通过MTPA运算获得dq轴定子电流给定值
Figure FDA00032398043000000216
和
Figure FDA00032398043000000217
其中，MTPA运算公式如下：
Figure FDA00032398043000000218
式中，
Figure FDA00032398043000000219
为永磁体磁链，Lq、Ld分别表示q轴电感和d轴电感，sgn表示符号函数；
步骤6：将步骤5获得的dq轴定子电流给定值
Figure FDA00032398043000000220
和
Figure FDA00032398043000000221
经过电流控制器输出dq轴电压给定值
Figure FDA00032398043000000222
和
Figure FDA00032398043000000223
再经过坐标变换得到静止旋转坐标系下的电压给定值
Figure FDA00032398043000000224
和
Figure FDA00032398043000000225
并将
Figure FDA00032398043000000226
和
Figure FDA00032398043000000227
输入SVPWM进行脉宽调制，输出PWM波触发逆变器，实现航空柱塞泵用内嵌式永磁同步电机的多重抗扰控制。
2.根据权利要求1所述航空柱塞泵用内嵌式永磁同步电机的多重抗扰控制方法，其特征在于：所述步骤2的航空柱塞泵周期性脉动载荷由瞬时脉动流量对时间求导数获得，离散域的表达式为：
Figure FDA0003239804300000031
式中，K为比例系数，
Figure FDA0003239804300000032
为单周期脉动流量变化值，其表达式为
Figure FDA0003239804300000033
Figure FDA0003239804300000034
和
Figure FDA0003239804300000035
分别表示k时刻和k-1时刻的脉动流量，Δt为周期间隔时间。
3.根据权利要求1所述航空柱塞泵用内嵌式永磁同步电机的多重抗扰控制方法，其特征在于：所述步骤2的非周期性冲击载荷的离散域模型为：航空柱塞泵非周期性冲击载荷由外部冲击主载荷和内部附加载荷构成，其中内部附加载荷包括滑动摩擦载荷、粘性摩擦载荷和滚动摩擦载荷，表达式如下：
Figure FDA0003239804300000036
式中，d(k)为所述航空柱塞泵外部冲击主载荷，ΔTt(k)为所述内部附加载荷，表达式为ΔTt(k)＝Tf(k)+Tv(k)+Ts(k)，其中Tf(k)为滑动摩擦载荷，Tv(k)为粘性摩擦载荷，Ts(k)为滚动摩擦载荷。
4.根据权利要求1所述航空柱塞泵用内嵌式永磁同步电机的多重抗扰控制方法，其特征在于：所述步骤4的Luenberger动态带宽负载转矩观测器的增益系数k1、k2获取方式如下：
1)建立观测器特征方程
Figure FDA0003239804300000037
依据特征方程获得初始增益系数k′1和k′2的表达式：
k′1＝Jω‘
Figure FDA0003239804300000038
式中，ω’为设定的观测器带宽，Bm为电机摩擦系数，J为电机惯性系数；
2)依据所述初始增益系数k′1和k′2，建立基于双曲正切函数的增益系数k1、k2表达式为：
Figure FDA0003239804300000041
Figure FDA0003239804300000042
式中，β1、β2、c1、c2均为比例系数，
Figure FDA0003239804300000043
为观测器的转速估计误差。
CN202111013500.5A 2021-08-31 2021-08-31 航空柱塞泵用内嵌式永磁同步电机的多重抗扰控制方法 Active CN113794411B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111013500.5A CN113794411B (zh) 2021-08-31 2021-08-31 航空柱塞泵用内嵌式永磁同步电机的多重抗扰控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111013500.5A CN113794411B (zh) 2021-08-31 2021-08-31 航空柱塞泵用内嵌式永磁同步电机的多重抗扰控制方法

Publications (2)

Publication Number Publication Date
CN113794411A CN113794411A (zh) 2021-12-14
CN113794411B true CN113794411B (zh) 2023-02-07

Family

ID=78876772

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111013500.5A Active CN113794411B (zh) 2021-08-31 2021-08-31 航空柱塞泵用内嵌式永磁同步电机的多重抗扰控制方法

Country Status (1)

Country Link
CN (1) CN113794411B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09121580A (ja) * 1995-10-26 1997-05-06 Meidensha Corp 低慣性化制御による２慣性共振系の振動抑制装置
JPH10155292A (ja) * 1996-11-22 1998-06-09 Meidensha Corp ２慣性系制御回路
JPH11341852A (ja) * 1998-02-25 1999-12-10 Yoshihiro Matsui 電動機の速度制御装置
CN108832863A (zh) * 2018-06-20 2018-11-16 无锡信捷电气股份有限公司 一种双观测器的伺服系统谐振抑制方法
CN111342720A (zh) * 2020-03-06 2020-06-26 南京理工大学 基于负载转矩观测的永磁同步电机自适应连续滑模控制方法
CN112054729A (zh) * 2020-09-11 2020-12-08 洛阳中重自动化工程有限责任公司 一种适用于低速直驱提升机的永磁电机控制方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09121580A (ja) * 1995-10-26 1997-05-06 Meidensha Corp 低慣性化制御による２慣性共振系の振動抑制装置
JPH10155292A (ja) * 1996-11-22 1998-06-09 Meidensha Corp ２慣性系制御回路
JPH11341852A (ja) * 1998-02-25 1999-12-10 Yoshihiro Matsui 電動機の速度制御装置
CN108832863A (zh) * 2018-06-20 2018-11-16 无锡信捷电气股份有限公司 一种双观测器的伺服系统谐振抑制方法
CN111342720A (zh) * 2020-03-06 2020-06-26 南京理工大学 基于负载转矩观测的永磁同步电机自适应连续滑模控制方法
CN112054729A (zh) * 2020-09-11 2020-12-08 洛阳中重自动化工程有限责任公司 一种适用于低速直驱提升机的永磁电机控制方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Equivalence of resonance ratio and I-PD controllers in zero stiffness torque control for soft robot;Yusuke Kawai et al.;《IECON 2015 - 41st Annual Conference of the IEEE Industrial Electronics Society》;20160128;第1424-1429页 *
基于PR控制的负载模拟系统;王哲 等;《中国电机工程学报》;20150120;第35卷(第2期);第486-493页 *

Also Published As

Publication number Publication date
CN113794411A (zh) 2021-12-14

Similar Documents

Publication Publication Date Title
CN102969968B (zh) 一种永磁同步电机控制方法
EP3014123B1 (en) A pump system
CN110323974B (zh) 一种基于比例谐振控制器优化的自抗扰控制方法
CN103701386B (zh) 基于观测磁链误差的异步电机无速度传感器的全阶磁链观测器的获取方法
CN108459497B (zh) 一种基于adrc和nlpid的稳瞄伺服系统的控制方法
CN108336935B (zh) 一种反步控制协同eso的直线电机控制方法
JP2019527532A (ja) モータのベクトル制御方法、装置及び航空機
EP3014122B1 (en) Anti-ripple injection method and pump system
CN109347391A (zh) 一种考虑系统噪声的朗道自适应转动惯量辨识方法
CN114687899B (zh) 一种基于滑模自抗扰的航空电动燃油泵流量控制方法
CN112671288A (zh) 一种记忆电机调磁转矩脉动抑制方法
CN111726048B (zh) 基于滑模观测器的永磁同步电机转子位置和速度估算方法
CN113794411B (zh) 航空柱塞泵用内嵌式永磁同步电机的多重抗扰控制方法
CN111082720B (zh) 一种直驱式航空电动燃油泵鲁棒控制器
CN113630052A (zh) 一种直线振荡电机无位置传感器控制方法及装置
CN112072973A (zh) 一种基于预测自适应律的永磁同步电机超扭滑模控制方法
CN115694284A (zh) 一种无刷双馈电机转矩脉动抑制方法
CN112727745B (zh) 压缩机的控制方法、装置、制冷设备和可读存储介质
JP2004064949A (ja) 超音波モータの位置制御方式
CN114465543A (zh) 一种永磁同步电机无位置传感器控制方法
CN109586622B (zh) 基于滑模控制的水下推进器控制方法
Han et al. Research on PMSM sensor-less system based on ADRC-PBC strategy
Wang et al. Design on sliding mode controller with adaptive fuzzy compensation for switched reluctance motor drive systems
Zhu et al. Simulation Research on Ship-borne PMSM Speed Regulation Control System Based on Fuzzy PID Control
CN114721293B (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