一种基于滑模控制的高铁低频振荡抑制方法与流程

本发明涉及高铁供电技术领域没截图为一种基于滑模控制的高铁低频振荡抑制方法。

背景技术：

随着中国高铁的迅猛发展，高速、重载、大密度已经成为当下电气化铁路的重要特性，因此，对高铁的供电可靠性提出了越来越高的要求。越来越多的新型交直交传动机车投入运营，大大增加了牵引供电系统的复杂性，我国牵引供电系统的低频振荡现象就是在这样的背景下产生的，这种现象根本上是车网参数不匹配造成的不稳定现象。许多研究结果均表明，网侧整流器的控制策略是影响车网系统稳定性的一大重要因素。交流传动机车、高速动车组中所用的牵引变流器主要由四象限脉冲整流器和牵引逆变器组成。四象限脉冲整流即网侧整流器，属于pwm整流器。动车组网侧整流器通常采用电压外环和电流内环相结合的双闭环控制方式，目前主要的控制方式有瞬态电流控制策略、可变相位角控制策略、预测直接电流控制策略、dq解耦控制策略等。传统的线性控制方法的控制效果已经很难得到提升，因此将非线性控制方法，例如预测控制、无源控制以及滑模控制等引入变流器的控制中是有必要的。

技术实现要素：

基于上述原因，本发明的目的在于提供一种能够提高整流器的动态响应和鲁棒性能，并能抑制牵引网-动车组电气量低频振荡的发生的基于滑模控制的高铁低频振荡抑制方法，技术方案如下：

一种基于滑模控制的高铁低频振荡抑制方法，包括以下步骤：

步骤a：定义网侧电流正交量，构建d-q旋转坐标系下动车组网侧整流器的状态空间模型；

步骤b：结合控制目标，选取控制系统的输出，建立基于外环电压控制的两个滑模面，得到网侧电流无功分量的参考值，完成滑模面表达式的求取；

步骤c：对上述两个滑模面选取滑模控制趋近率，将步骤b建立的两个滑模面代入滑模控制趋近率，将得到的结果代入d-q旋转坐标系下动车组网侧整流器的状态空间模型，得到开关函数，将该开关函数与直流侧电压相乘，得到控制电压；

步骤d：将步骤c中得到的控制电压经过坐标变换得到α-β坐标系分量，再通过正弦脉宽调制输出控制脉冲。

进一步的，所述步骤a的具体过程如下：

针对单相两电平拓扑结构，将交流侧电流in表示为iα，交流侧电动势en表示为uα，构建与iα、uα正交的虚拟量iβ、eβ，根据坐标变换遵循的等功率变换原则，得到

uαiα+uβiβ＝udcidc

式中，udc为直流侧稳定电压，idc为直流侧稳定负载电流；

定义uα＝sαudc，uα是uab在α轴上的分量，uβ＝sβudc，uβ是uab在β轴上的分量，sα、sβ是α-β坐标系下的理想开关函数；得到idc的表达式：

结合交流侧的kvl方程，将交流量、虚拟正交量表示在同一方程组中，得到

式中，l和r分别为车载变压器等效到次边的漏电感和漏电阻，c为直流侧支撑电容；进行两相静止坐标系到两相旋转坐标系的坐标变换，得到d-q旋转坐标系下动车组网侧整流器的数学模型：

式中，ed、eq分别为网侧电动势与虚拟正交分量分别在旋转坐标系中d轴与q轴的分量；id、iq分别为网侧电动势与虚拟正交分量分别在旋转坐标系中d、q分量；根据定义，d轴分量表示有功分量，q轴分量表示无功分量，sdudc、squdc则是uab相应的分量，表示成状态空间模型的矩阵形式：

得到动车组网侧整流器在两相旋转坐标系下矩阵形式的数学模型。

更进一步的，所述步骤b的具体过程如下：

选择udc、iq为控制系统的输出，设则根据步骤a中得到的d-q旋转坐标系下动车组网侧整流器状态空间模型的矩阵形式，推出其能控标准型

udc、iq、的误差值eφ分别表示为：

建立两个滑模面s1、s2分别与udc、iq对应：

其中，α、α1、α2为放大增益，并用β代替作为滑模控制的反馈系数；

将误差值eφ的具体表达式代入上式，得到

将d-q旋转坐标系下动车组网侧整流器的数学模型以及直流侧电压参考值代入上式得到

将id分离，得到

式中，c、是已知的常数，udc、iq、il可以通过对应的测量工具测得，β为设计者设计的参数，sd与sq是开关函数，将作为id的参考值；

根据d-q旋转坐标系下遵循功率不变原则

udid+uqiq＝udcidc

根据控制目标，无功电流分量iq＝0，则

udid＝udcidc

定义ud＝sdudc，由于ud＝ed-rid，上式表示为

(ed-rid)id＝udcsdid

则sd的表达式为

由d-q旋转坐标系下动车组网侧整流器的数学模型

由于容量较大的变压器的电阻r极小，故忽略，又因为在理想情况下，eq＝0，则

得到

将上式代入id的参考值的表达式，得到

则

将滑模面改写为：

更进一步的，所述步骤c的具体过程如下：

对上所述两个滑模面s1、s2选取指数趋近率，得到

ε1为常数，表示系统状态点趋近于滑模面s1(x)＝0的速率，速率与ε1的大小成正比；-k1s为指数趋近项，代表系统状态趋近滑模面的将滑模面的过程；为滑模面s1的求导；

将滑模面的表达式代入

由于在理想状态下为常数，故则上式表示为

代入d-q旋转坐标系下动车组网侧整流器的状态空间模型，得到

将开关函数sd、sq表示出来

得到脉冲宽度调制所需控制电压

本发明的有益效果是：

1)本发明提出的适用于动车组网侧整流器的滑模控制策略因其本身是变结构控制的特性，符合pwm整流器在本质上也是变结构系统的特征，提高了该系统的鲁棒性能，能够阻尼牵引网-动车组电气量振荡的问题。

2)本发明所需要设置的参数较少，控制系统且对参数变化不敏感，在有效地范围内该控制具有较好的鲁棒性和动态性能。

3)本发明控制策略结构简单，能够有效提高系统的响应速度。

附图说明

图1为本发明的控制结构图。

图2为本发明的整流器等效电路图。

图3位本发明的整流器包含开关函数的简化电路图。

图4为本发明的仿真模型图，图中from与goto标号相同为一组，表示两个信号实际上相连；un为交流侧等效电源，rn为交流侧等效电阻，ln为交流侧等效电感；u_n为交流侧电压表，ln_m为交流侧电流表，ln_m1为负载电流表，u_dc为负载电压表，cd为稳压电容，rd为负载等效电阻。

图5为本发明的控制器仿真模型图。

图6a为本发明直流侧电压udc波形图。

图6b交流侧电流in、交流侧电流un波形图。

图6c为本发明车网级联仿真在0s时刻首先接入5台动车，接着在3s、6s以及9s时刻分别增加1台动车，一共接入8台动车的网侧电压u、网侧电流i波形图。

具体实施方式

下面结合附图和具体实施例对本发明做进一步详细说明，本实施例以crh5型动车组为例。

步骤a：定义网侧电流正交量，构建d-q旋转坐标系下动车组网侧整流器的状态空间模型。

动车组网侧整流器采用的是四象限脉冲整流器，本专利针对的是单相两电平拓扑结构。受电弓从接触网取流，经车载变压器降压后作为整流器的输入，整流器则将输入的单相交流电压变换成稳定的直流电压。

将交流侧电流in表示为iα，交流侧电动势en表示为uα，构建与iα、uα正交的虚拟量iβ、eβ，根据坐标变换遵循的等功率变换原则，可得到

uαiα+uβiβ＝udcidc

式中，udc为直流侧稳定电压，idc为直流侧稳定负载电流。

定义uα＝sαudc，uα是uab在α轴上的分量，uβ＝sβudc，uβ是uab在β轴上的分量，sα、sβ是α-β坐标系下的理想开关函数。得到idc的表达式：

如图2所示通过对交流侧列写kvl方程，并将交流量、虚拟正交量表示在同一方程组中，可获得d-q旋转坐标系下动车组网侧整流器的数学模型

式中，l和r分别为车载变压器等效到次边的漏电感和漏电阻，c为直流侧支撑电容。动车组网侧整流器中，相应的参数取值分别为：牵引网网压有效值un＝27500v,网侧电阻r＝0.145ω,网侧电感l＝0.0054h，直流侧支撑电容c＝0.009f，直流侧负载等效电阻r＝25ω。给定直流侧电压udc＝3600v。

对上式进行两相静止坐标系到两相旋转坐标系的坐标变换，可得到d-q旋转坐标系下动车组网侧整流器的数学模型

式中，ed、eq分别为网侧电动势与虚拟正交分量分别在旋转坐标系中d轴与q轴的分量，同理，id、iq分别为网侧电动势与虚拟正交分量分别在旋转坐标系中d、q分量。根据定义，d轴分量表示有功分量，q轴分量表示无功分量。sdudc、squdc则是uab相应的分量，表示成状态空间模型的矩阵形式

得到动车组网侧整流器在两相旋转坐标系下矩阵形式的数学模型。

步骤b：结合控制目标，选取控制系统的输出，建立基于外环电压控制的两个滑模面，得到网侧电流无功分量的参考值，完成滑模面表达式的求取。具体如下：

选择udc、iq为控制系统的输出。设则根据步骤a中得到的d-q旋转坐标系下动车组网侧整流器状态空间模型的矩阵形式，推出其能控标准型

根据动车组网侧整流器的两个控制目标：其一，保证网压网流保持功率因数接近于1；其二，保证直流侧电压稳定在允许的波动范围内，选择udc、iq为控制系统的输出。udc、iq、的误差值eφ分别表示为

建立两个滑模面s1、s2分别与两个控制目标即udc、iq对应

其中，α、α1、α2为放大增益，并使用β代替为滑模控制的反馈系数。

将d-q旋转坐标系下动车组网侧整流器的数学模型以及人为设置的直流侧电压参考值代入上式可得，此处由于crh5型车直流侧电压的额定值3600v，则控制目标要直流侧电压在3600v附近一定范围内波动。

将id分离，得到

式中c、是已知的常数，udc、iq、il可以通过对应的测量工具测得，β是需要设计者设计的参数，sd与sq是开关函数，将作为id的参考值。

根据d-q旋转坐标系下遵循功率不变原则

udid+uqiq＝udcidc

根据控制目标，无功电流分量iq＝0，则

udid＝udcidc

定义ud＝sdudc，根据电路知识可得ud＝ed-rid,得到sd的表达式

由d-q旋转坐标系下动车组网侧整流器的数学模型。且对于容量较大的变压器，其电阻r极小，可以忽略，又因为在理想情况下，eq＝0，则

将上式代入id的参考值的表达式，得到

则

可将滑模面改写成以下简单的形式

步骤c：对上述两个滑模面选取滑模控制趋近率，将步骤b建立的两个滑模面代入滑模控制趋近率，将得到的结果代入d-q旋转坐标系下动车组网侧整流器的状态空间模型，得到开关函数，将该开关函数与直流侧电压相乘，得到控制电压。

对上文所取得两个滑模面s1、s2选取指数趋近率，得到

ε1为常数，表示系统状态点趋近于滑模面s1(x)＝0的速率，速率与ε1的大小成正比；-k1s是指数趋近项，代表系统状态趋近滑模面的将滑模面的过程；为滑模面s1的求导值；ε2、s2、-k2s2、含义规律同上。

将滑模面的表达式代入

由于都是常数，所以上式可表示为

代入d-q旋转坐标系下动车组网侧整流器的状态空间模型，得到

将开关函数sd、sq表示出来

得到脉冲宽度调制所需控制电压

d.将步骤c中得到的控制电压经过坐标变换得到α-β坐标系分量，再通过正弦脉宽调制输出控制脉冲。

最后在matlab/simulink中搭建仿真模型如图4所示，所得电压、电流波形如图6a、6b所示，所加负载为额定负载，等效电阻为r＝25ω。波形显示直流侧电压超调量为0，没有峰值时间，调节时间为0.05s，电压波动60v。相比常用的传统比例积分控制而言性能指标得到较好改善，交流电流从启动到稳定仅需要3个周波,且thd有所减小。

将该滑模控制算法应用于牵引网-动车组级联仿真模型中，依次增加接入牵引网的动车组数量，在传统瞬态直接电流控制下，接入动车组达到8台时动车组和牵引网电压、电流发生明显的波动，及产生车网低频振荡现象。在基于滑模控制高铁低频振荡抑制方法的控制下，接入动车组达到甚至超过8台时，电气量基本稳定，如图6c所示，未发生低频振荡问题。