一种多电机总量协同有限时间抗饱和控制方法与流程

[0001]
本发明涉及多电机牵引系统领域，更具体地，涉及一种多电机总量协同有限时间抗饱和控制方法。


背景技术：

[0002]
多永磁同步电机共同提供动力是电力机车常见的驱动方式。机车复杂多变的运行环境常常引起一个或多个电机牵引性能的损失。经典的控制是运用同步控制或一致性算法，实现系统中个体状态(如速度、位置等)的同步。而近几年为了保证机车整体牵引性能，总量一致的理论在多电机牵引系统中被提出，即实现各电机输出转矩之和与总牵引特性曲线的一致。但以上所提出的控制率大多只能达到渐近稳定的结果，而有限时间控制则收敛更快，精度更高和鲁棒性更强。结合加幂积分终端滑模技术，保证了控制器的连续与非奇异，又可严格求取有限时间上界。但加幂积分技术由于其在证明过程中大量使用了不等式放缩，导致其参数的强约束。因此，针对加幂积分参数简化方面进行研究，并保证系统全局有限时间稳定，就可有效增强控制策略的工程实用性。
[0003]
同时，考虑工程中各电机的参数摄动和负载转矩扰动，构建了基于滑模变结构的牵引转矩控制，以提高各电机的动态响应性能。但其较大的滑模切换增益和加幂积分中参数的约束，极易导致控制输入饱和问题。饱和问题是一个非线性问题，会影响系统的动态性能，致使系统不稳定，甚至造成器件的损坏。目前针对输入饱和问题的研究非常广泛：利用数学上相关函数处理饱和；设计静态或动态抗饱和补偿器削弱饱和影响；将饱和约束转化为带有线性矩阵不等式约束的优化问题等等。针对复杂的多电机牵引系统，滑模切换增益及加幂积分参数的强约束均易引起控制输入饱和问题，尤其在基于总量一致的理论中，多电机的输入饱和问题更加突出，更加复杂。因此，结合先进的抗饱和技术，解决复杂多电机系统中的饱和问题就极具有工程意义。


技术实现要素：

[0004]
本发明针对现有技术中复杂的多电机牵引系统，特别是基于总量一致下的输入饱和问题，提供一种多电机总量协同有限时间抗饱和控制方法。
[0005]
为实现上述目的，本发明采用以下技术方案：
[0006]
一种多电机总量协同有限时间抗饱和控制方法，包括以下步骤：
[0007]
s1、依据总量一致的理论以及永磁同步电机旋转坐标系下的方程，考虑实际中电阻、电感变化引起的参数摄动及负载转矩扰动，推导出电机的状态方程；
[0008]
s2、根据电机状态方程的参数，设计有限时间收敛的辅助抗饱和系统；
[0009]
s3、根据s1中电机状态方程和总量一致的理论，构造了误差动力学方程；
[0010]
s4、根据s2中有限时间收敛的辅助抗饱和系统和s3中所述误差动力学方程，基于非奇异终端滑模设计总量协同有限时间控制器，并简化了加幂积分参数；
[0011]
总量协同有限时间控制器为：
[0012][0013]
其中，抗饱和系数c
2j
＞0，滑模参数k
j
≥α
11/q-2
d
j
，加幂积分参数简化为α1＞0，α2＞0。
[0014]
s5、根据加幂积分技术和有限时间李雅普诺夫稳定定理完成总量协同有限时间控制器的稳定性证明并且求得有限时间上界。
[0015]
进一步地，在步骤s1中电机的状态方程为：
[0016][0017]
其中，其中，，
[0018]
d
2j
＝-δa
1j
x
1j-δa
2j
x
2j
+δb
j
u
qj
+f
j
(t)。
[0019]
进一步地，在步骤s2中有限时间收敛的辅助抗饱和系统为：
[0020][0021]
其中，x
aj
为辅助状态，a
aj
为待设计的正常数，满足τ是一个小的正常数，δu
qj
＝u
qj-v
qj
，常数k4＞0。
[0022]
进一步地，在步骤s3在中误差动力学方程为：
[0023][0024]
进一步地，在步骤s5中求得的时间上界为：
[0025]
t
s
≤2q/[k8(q-1)]
·
v(x0)
(q-1)/(2q)
[0026]
进一步地，在步骤s1中永磁同步电机数学模型为隐极式永磁同步电机。
[0027]
进一步的，隐极式永磁同步电机初始转矩值设置为0.25，控制输入饱和设置为+220—-220。
[0028]
进一步地，在步骤s5中采用有限时间李雅普诺夫稳定定理完成总量协同有限时间抗饱和控制器的稳定性证明。
[0029]
本发明的有益效果为：
[0030]
1、将多电机协同控制从个体状态一致扩展到总牵引转矩一致，保证了机车在实际运行环境下的整体牵引性能；
[0031]
2、简化了加幂积分参数，降低了控制输入饱和的可能；
[0032]
3、设计了有限时间辅助抗饱和系统，提高了抗饱和的快速性；
[0033]
4、设计了多电机总量协同有限时间抗饱和控制器，实现了全局有限时间稳定，削
弱了输入饱和对整体牵引性能的影响。
附图说明
[0034]
图1为系统框架图；
[0035]
图2为加幂积分参数强约束的控制输入曲线图；
[0036]
图3为加幂积分参数弱约束的控制输入曲线图；
[0037]
图4为加幂积分参数简化后的控制输入曲线图；
[0038]
图5为无抗饱和的跟踪效果曲线图
[0039]
图6为无抗饱和的跟踪误差曲线图；
[0040]
图7为渐近稳定抗饱和的跟踪效果曲线图；
[0041]
图8为渐近稳定抗饱和的跟踪误差曲线图；
[0042]
图9为有限时间抗饱和的跟踪效果曲线图；
[0043]
图10为有限时间抗饱和的跟踪误差曲线图；
[0044]
图11为有限时间抗饱和后的控制输入曲线图；
[0045]
图12为rt-lab半实物实验平台连接图；
[0046]
图13为实验的跟踪效果曲线图；
[0047]
图14为实验的跟踪误差曲线图；
[0048]
图15为实验的控制输入曲线图
具体实施方式
[0049]
下面结合具体实施方式对本发明作进一步的说明。
[0050]
一种多电机总量协同有限时间抗饱和控制方法：
[0051]
s1、考虑参数摄动和负载转矩扰动，依据永磁同步电机旋转坐标系下的方程，推导出电机的状态方程；旋转坐标系下的方程为：
[0052][0053]
在多电机系统中对单个电机采用i
dj
＝0的控制方式，则可得如下方程：
[0054]
[0055]
选取状态变量为输出角速度ω
j
＝x
1j
，角加速度电磁转矩t
ej
＝x
3j
，化简可得：
[0056][0057]
其中，为时变扰动，
[0058]
考虑电机在实际运行过程中电阻、电感参数随温度而变化。所以，将参数摄动与负载转矩扰动统一归结为未知复合干扰，则状态方程为：
[0059][0060]
其中，其中，其中，d
2j
＝-δa
1j
x
1j-δa
2j
x
2j
+δb
j
u
qj
+f
j
(t)。
[0061]
s2、设计有限时间收敛的辅助抗饱和系统；
[0062][0063]
其中，x
aj
为辅助状态，a
aj
为待设计的正常数，满足τ是一个小的正常数，δu
qj
＝u
qj-v
qj
，常数k4＞0。
[0064]
s3、基于总量一致理论定义n个电机的输出转矩总和与给定的牵引特性曲线的偏差为：
[0065][0066]
令则可得如下误差动力学方程为：
[0067][0068]
s4、根据有限时间收敛的辅助抗饱和系统，基于非奇异终端滑模设计总量协同有限时间控制器，并简化了加幂积分参数；选取非奇异终端滑模面为：
[0069][0070]
因此，结合加幂积分技术，设计基于非奇异终端滑模的总量协同抗饱和有限时间控制器为：
[0071][0072]
其中，系数c
2j
＞0，参数加幂积分参数简化为α1＞0，α2＞0即可。
[0073]
s5、根据加幂积分技术和有限时间李雅普诺夫稳定定理完成总量协同有限时间控制器的稳定性证明并且求得有限时间上界；
[0074]
构造正定的李雅普诺夫函数如下：
[0075]
v＝v1+v2+v3[0076]
其中，并且，虚拟控制律设计为然后，对v2化简可得到
[0077]
步骤1：对v1求导，
[0078][0079]
化简可得：
[0080][0081]
其中，λ1＝2
1/q-1
γ1(1+q)/q，γ1＞0，
[0082]
步骤2：对v2求导，
[0083][0084]
先对上式右边第一项代入并化简可得：
[0085][0086]
其中，
[0087]
然后对右边第二项代入与控制器，并化简可得：
[0088][0089]
其中，
[0090]
合并公式可得：
[0091][0092]
步骤3：对v3求导，
[0093]
通过不等式并代入辅助抗饱和系统可得：
[0094][0095]
综上三步：对v求导，
[0096][0097]
其中，k＝k
2-k3，
[0098]
当参数满足γ1+γ2+γ3＜1，k1＞0，k2＞k3，k4＞0则：
[0099][0100]
对v2化简分析得考虑|σ1|＝2
1/2
v
11/2
和可得：
[0101][0102]
其中，k5＝2
(q+1)/(2q)
k1，k7＝2
(q+1)/(2q)
k4，k8＝min{k5,k6,k7}。
[0103]
由上式满足有限时间李雅普诺夫稳定定理，则可得有限时间上界为：
[0104]
t
s
≤2q/[k8(q-1)]
·
v(x0)
(q-1)/(2q)
[0105]
因此，σ1，s，x
aj
及σ2均能在有限时间t
s
内收敛到0。
[0106]
同时，仿真和实验验证了本方法的有效性和可行性。
[0107]
本发明仿真采用4台参数不同的永磁同步电机组成的多电机牵引系统为仿真对象，各电机参数如下表所示。其中，为便于与本团队前期成果进行对比，各电机控制输入饱和仍设置为
±
220，各电机初始转矩值设置为0.25。
[0108]
参数电机1电机2电机3电机4电阻r
s
/ω2.8732.8562.8672.861电感l/mh8.78.68.98.8极对数n
p
2222转动惯量j/kg
·
m24.5*10-4
4.47*10-4
4.51*10-4
4.49*10-4
摩擦系数r
ω
/n
·
m
·
s4.831*10-5
4.846*10-5
4.827*10-5
4.838*10-5
永磁体磁链ψ
f
/wb0.1750.1730.1780.177
[0109]
考虑实际中未知复合干扰的不确定性，采用4种不同的干扰信号(突变干扰、缓变干扰、高频噪声和均匀噪声)分别作用于4台电机中。同时，给定牵引特性曲线由分段函数组成，如下面式子所示：在0
→
5s模拟电机启动阶段；在5
→
10s模拟电机匀速运转阶段；在10
→
15s模拟电机减速停机阶段。
[0110][0111]
通过不同的加幂积分参数α1，α2，分析各电机控制输入波形的变化，从而验证本文简化加幂积分参数对降低控制输入饱和方面的有效性。
[0112]
通过图2，图3和图4可以看出加幂积分参数约束会严重影响控制输入曲线的变化，
为方便分析，将不同加幂积分参数对应的各电机控制输入峰值的绝对值整理成下表所示。
[0113][0114]
上表中粗体字表示控制输出超出了饱和限幅220。对上表进行分析，可明显看出传统的加幂积分方法对于参数α1,α2均有较强的约束，如第一栏和第二栏，在多电机系统中极易引起饱和问题。而本文放宽了加幂积分参数的约束，即当取α1＝1时，由参数设置依据k2＞k3可得α2＞5.35时就可保证有限时间稳定，如上表的第三栏。该方法可以明显降低多电机中出现的饱和问题。但由于参数设置的依据，仍可能造成小幅的超限，所以进一步结合辅助抗饱和系统，就可有效解决输入饱和的问题。
[0115]
通过无抗饱和、渐近稳定的抗饱和，以及本申请有限时间抗饱和的总量协同跟踪控制三者的对比仿真，展现了本申请所提控制方法的优良性。同时，总量协同有限时间抗饱和的控制中加幂积分参数按α1＝1和α2＝6来设置，进一步体现了本申请结合两个角度解决输入饱和问题的优越性。图5-图10分别为无抗饱和、渐近稳定的抗饱和，以及本申请有限时间抗饱和的总量协同跟踪控制效果图。
[0116]
分析图5-图10可得，当总量协同控制中出现电机输入饱和时，图5、图6的无抗饱和策略会造成大约5倍左右的冲击，严重影响跟踪性能；图7、图8的渐近稳定抗饱和策略会造成大约0.2倍左右的冲击；而由图9、图10可得，在本申请的控制策略下，输入饱和引起的冲击在千分之五左右，几乎可以忽略不计。所以，本申请提出的总量协同有限时间抗饱和控制器可以有效削弱输入饱和对整体牵引性能的影响。同时，图9、图10抗饱和控制参数a
aj
＝50，c
2j
＝100相比于图7、图8中抗饱和控制参数a
aj
＝1000，c
2j
＝1900更具有工程实用性。最后，由图11有限时间抗饱和后的控制输入曲线可知，各个电机均能有效控制在输入限幅220以内。仔细分析，在11.22秒左右电机1出现饱和后，在辅助抗饱和系统和总量协同控制的作用下，于11.69秒左右电机2又出现了饱和，所以这体现了在多电机总量协同控制中，输入饱和问题的复杂性，而本申请所提的控制策略仍可有效削弱多电机输入饱和对整体牵引性能的影响。
[0117]
为了验证本文所提控制策略的实用性，保证系统仿真尽可能贴近实际工程环境，本申请进行了半实物实验。rt-lab半实物实验平台如图12所示，它包括：tms320f2812数字信号生成器、rt-lab op5600仿真电机、个人计算机、示波器及相关连接线等。实验参数和条件的设置与上述仿真参数相同。实验结果如图13-15所示。
[0118]
在有限时间抗饱和的总量协同跟踪控制策略下，图13为实验的跟踪效果图，图14为实验的跟踪误差图。可以很明显看出，输入饱和引起的千分之五左右的冲击对于整体牵引性能几乎毫无影响，这与simulink仿真结果一致。图15为实验的各个电机控制输入曲线图，表明在本文的控制策略下，各电机控制输入均能有效保持在饱和受限幅值以内。总而言之，图13-15表明，本文所提方法的实验结果与simulink仿真结果基本一致，为实际工程应
用提供了一定的理论依据。
[0119]
以上所述实施例仅是为充分说明本发明而所举的较佳的实施例，本发明的保护范围不限于此。本技术领域的技术人员在本发明基础上所作的等同替代或变换，均在本发明的保护范围之内。