CN104932259A - 一种电液位置伺服系统的增益自调节的超螺旋滑模控制方法 - Google Patents

一种电液位置伺服系统的增益自调节的超螺旋滑模控制方法 Download PDF

Info

Publication number
CN104932259A
CN104932259A CN201510261197.9A CN201510261197A CN104932259A CN 104932259 A CN104932259 A CN 104932259A CN 201510261197 A CN201510261197 A CN 201510261197A CN 104932259 A CN104932259 A CN 104932259A
Authority
CN
China
Prior art keywords
centerdot
alpha
beta
formula
kappa
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.)
Pending
Application number
CN201510261197.9A
Other languages
English (en)
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.)
Nanjing University of Science and Technology
Original Assignee
Nanjing University of Science and Technology
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 Nanjing University of Science and Technology filed Critical Nanjing University of Science and Technology
Priority to CN201510261197.9A priority Critical patent/CN104932259A/zh
Publication of CN104932259A publication Critical patent/CN104932259A/zh
Pending legal-status Critical Current

Links

Landscapes

  • Feedback Control In General (AREA)

Abstract

本发明公开了一种电液位置伺服系统的增益自调节的超螺旋滑模控制方法。利用已知的系统模型信息,在传统超螺旋滑模控制算法中引入基于模型的前馈控制律,提升系统伺服精度。采用自适应律实时更新控制器增益,无需先验地获知系统建模不确定性的确切界,避免了传统算法中由人为设定与该界相关的控制增益造成的保守性。基于Lyapunov稳定性理论证明了闭环系统全局稳定,系统跟踪误差可在有限时间内渐近收敛到零附近任意小的范围内,且收敛的速度和稳态误差的界可通过参数进行调节。

Description

一种电液位置伺服系统的增益自调节的超螺旋滑模控制方法
技术领域
本发明涉及电液伺服控制技术领域,主要涉及一种电液位置伺服系统的增益自调节的超螺旋滑模控制方法。
背景技术
电液伺服系统具有功率密度大、输出力/力矩大及抗负载刚性强等突出优点,在众多重要领域有着广泛的应用。随着现代工业技术的快速发展,对于电液伺服系统的伺服精度要求也越来越高。电液伺服系统是一个复杂的非线性控制对象,存在由于工作环境变化(如温度变化)和组件磨损等引起的系统参数不确定性以及非线性摩擦和外负载干扰等不确定性非线性,这些建模不确定性会严重恶化所设计的控制器的性能,导致系统出现跟踪误差、极限环振荡甚至失稳。为提升伺服系统的性能,需要设计高性能的控制器以抑制或消除建模不确定性对系统性能的影响。针对电液伺服系统高性能控制问题,许多方法被广泛研究。自适应控制通过设计参数自适应律实时更新系统的参数值以提高模型补偿精度,从而使系统获得良好的跟踪性能。然而自适应控制器设计都是基于系统只存在参数不确定性的前提,然而实际系统不可避免地存在难以建模的动态及外干扰等不确定性非线性,将会使所设计的自适应控制器性能降阶。
滑模控制是另一类有效的控制方法。传统的滑模控制可以处理所有有界的建模不确定性并获得渐近跟踪的性能,但其最大的缺点是存在的不连续的符号函数会造成控制器抖振,这对于实际系统是不允许的。为解决传统滑模控制器不连续的问题,可采用连续的饱和函数代替不连续的符号函数可有效地避免控制输入抖振,然而却只能保证跟踪误差有界,丧失了渐近跟踪的性能。另外,高阶滑模控制器在保证控制器的连续性的同时还可获得渐近跟踪,但是控制器的设计需要滑模变量的导数的信息,这在实际中往往是认为不可获知的,因此不易工程实现。超螺旋滑模控制是一种特殊的二阶滑模控制方法,其设计只需要滑模变量本身的信息。
发明内容
基于以上分析,针对存在各种建模不确定性的电液位置伺服系统,本发明提出一种能有效地抑制建模不确定性对系统性能影响的电液位置伺服系统的增益自调节的超螺旋滑模控制方法。
实现本发明目的的技术解决方案为:一种电液位置伺服系统的增益自调节的超螺旋滑模控制方法,包括以下步骤:
步骤1,建立电液位置伺服系统的数学模型;
步骤2,设计具有自适应增益的超螺旋滑模控制器;
步骤3,具有自适应增益的超螺旋滑模控制器性能及稳定性测试。
本发明利用已知的系统模型信息,在传统超螺旋滑模控制算法中引入基于模型的前馈控制律,提升系统伺服精度。该方法无需知道建模不确定性的确切界,而是设计自适应律不断调整与该界相关的控制器增益。与现有技术相比,其显著优点是:控制器的设计无需先验地获知系统建模不确定性的确切界,而是采用自适应律实时更新控制器增益,避免了由人为设定与该界相关的控制器增益造成的保守性,同时所设计的控制器可保证跟踪误差在有限时间内渐近收敛到零附近任意小的范围。仿真结果验证了其有效性。
附图说明
图1是本发明电液位置伺服系统的原理图;
图2是具有自适应增益的电液位置伺服系统超螺旋滑模控制(ASTSC)方法原理示意图;
图3是ASTSC控制器作用下系统位置输出对期望输入指令的跟踪过程;
图4是三种控制器作用下系统的跟踪误差对比曲线;
图5是增益α随时间变化的曲线;
图6是传统滑模控制方法的控制输入信号;
图7是ASTSC和SSMC控制方法的控制输入信号。
具体实施方式
下面结合附图及具体实施例对本发明作进一步详细说明。
结合图1~2本发明电液位置伺服系统的增益自调节的超螺旋滑模控制方法,包括以下步骤:
步骤1,建立液电液位置伺服系统的数学模型;
(1.1)本发明所考虑的液电液位置伺服系统如图1所示,是通过伺服阀控制的液压马达驱动惯性负载。
图1左侧是电液位置伺服系统结构,右侧是液压马达结构示意图。
因此,根据牛顿第二定律,惯性负载的运动方程为:
J θ · · m = P L D m - B θ · m - f ( t ) - - - ( 1 )
式(1)中J和θm分别为负载转动惯量和转角;PL=P1-P2为液压马达的负载压力,P1和P2分别为液压马达进油腔和回油腔压力;Dm为液压马达的体积排量;B为粘性摩擦系数;f(t)为系统建模不确定性,包含未建模的非线性摩擦、外负载干扰等。
忽略液压马达的外泄漏,则液压马达负载压力动态方程为:
V t 4 β e P · L = - D m θ · m - C t P L + Q L - - - ( 2 )
式(2)中Vt为马达两腔的总控制容积;βe为有效油液弹性模量;Ct为马达腔室的内泄漏系数;QL为负载流量。
将伺服阀的动态近似为比例环节,即伺服阀阀芯位移与控制输入成比例,因此负载流量可表示如下:
Q L = k t u P s - sign ( u ) P L - - - ( 3 )
式(3)中kt为总的流量增益;Ps为系统供油压力且sign(u)定义如下:
sign ( u ) = 1 , u &GreaterEqual; 0 - 1 , u < 0 - - - ( 4 )
(1.2)假设系统建模不确定性f(t)是连续可微的,则基于式(1)、(2)和(3),定义状态变量那么系统模型可以写成如下的状态空间的形式:
x &CenterDot; 1 = x 2
x &CenterDot; 2 = x 3 - - - ( 5 )
x &CenterDot; 3 = f 1 ( P L , u ) u - f 2 ( x 2 ) - f 3 ( x 3 ) + d ( t )
式(5)中
f 1 ( P L , u ) = 4 D m &beta; e k t J V t P s - sign ( u ) P L ,
f 2 ( x 2 ) = 4 D m &beta; e J V t ( D m + C t B D m ) x 2 , f 3 ( x 3 ) = 4 &beta; e J V t ( C t J + V t B 4 &beta; e ) x 3 , - - - ( 6 )
d ( t ) = - C t D m k t f ( t ) - V t 4 D m &beta; e k t f &CenterDot; ( t )
对于系统物理参数J、Dm、B、βe、Ct、kt和Ps,尽管无法获取其精确值,但是可假设其名义值已知并用于控制器的设计,而将参数名义值与其真值之间的偏差归并到建模不确定性d(t)中。
系统控制器的设计目标为:给定系统参考信号yd(t)=x1d(t),设计一个连续且有界的控制输入u使系统输出y=x1尽可能地跟踪系统的参考信号。
步骤2,设计具有自适应增益的超螺旋滑模控制器,步骤如下:
(2.1)定义系统误差变量如下:
z1=x1-x1d
z 2 = x 2 - x &CenterDot; 1 d - - - ( 7 )
z 3 = x 3 - x &CenterDot; &CenterDot; 1 d
式(7)中z1、z2和z3分别为惯性负载位置、速度和加速度跟踪误差。
定义如下的滑模变量:
s=k1z1+k2z2+z3    (8)
为便于控制器设计,假设如下:
假设1:系统参考指令信号x1d(t)是三阶连续的,且系统期望位置指令、速度指令、加速度指令及加加速度指令都是有界的。液压马达位置伺服系统在一般工况下工作,即液压马达两腔压力P1,P2均小于供油压力Ps,且|PL|也小于Ps以保证式(6)中的f1>0。
假设2:系统建模不确定性d(t)满足以下条件
|d(t)|≤δ|s|1/2    (9)
式中:δ为未知的正数。
(2.2)对滑模变量求导可得:
s &CenterDot; = k 1 z 2 + k 2 z 3 + z &CenterDot; 3 = k 1 z 2 + k 2 z 3 + f 1 ( P L , u ) u - f 2 ( x 2 ) - f 3 ( x 3 ) + d ( t ) - x &CenterDot; &CenterDot; &CenterDot; 1 d - - - ( 10 )
基于式(10)中的滑模动态,设计超螺旋滑模控制器如下:
u = 1 f 1 ( P L , u ) ( u a + u s ) , u s = u s 1 + u s 2
u a = f 2 ( x 2 ) + f 3 ( x 3 ) + x &CenterDot; &CenterDot; &CenterDot; 1 d u s 1 = - k 1 z 2 - k 3 z 3 - &alpha; | s | 1 / 2 sign ( s ) - - - ( 11 )
u s 2 = - &Integral; 0 t &beta; 2 sign ( s ) d&tau;
式中:ua为用于改善模型补偿的前馈控制律;us为用于抑制建模不确定性d(t)对伺服系统性能影响的鲁棒控制律;为时变的控制器增益。
(2.3)设计如下的自适应律实时更新增益α和β:
&alpha; &CenterDot; = &gamma; 1 &kappa; 1 2 sign ( | s | - v ) - - - ( 12 )
β=2εα
式(12)中γ1、к1和ν都是任意的正数;ε为任意的实数。
由式(11)中控制器的结构可知,超螺旋滑模控制器的设计只依赖于滑模变量本身,而不需要其导数的信息,这是超螺旋算法与其他高阶滑模控制算法的本质区别。而且,由于控制算法中含有的符号函数经由积分运算以及与滑模变量绝对值函数的乘积运算,使得控制输入连续化,更利于工程中的实际执行。
步骤3,具有自适应增益的超螺旋滑模控制器性能及稳定性测试,具体如下:
考虑由式(5)描述的电液位置伺服系统,满足假设1和2,则利用式(11)中的超螺旋滑模控制器,且其控制增益α和β由式(12)中的自适应律进行更新,则对于任意的系统初始条件,存在有限时间tF>0使得滑模变量s渐近趋近于零附近任意小的范围内。
稳定性测速:
将式(11)代入(10)中可得
s &CenterDot; = - &alpha; | s | 1 / 2 sign ( s ) + &eta; + d ( t ) &eta; &CenterDot; = - &beta; 2 sign ( s ) - - - ( 13 )
定义新的状态变量
&xi; = [ &xi; 1 , &xi; 2 ] T = [ | s | 1 / 2 sign ( s ) , ] &eta; T - - - ( 14 )
结合式(13)可得
&xi; &CenterDot; = 1 2 | &xi; 1 | ( - &alpha; &xi; 1 + &xi; 2 + d ( t ) ) &xi; &CenterDot; 2 = - &beta; 2 | &xi; 1 | &xi; 1 - - - ( 15 )
由假设2可知,满足条件(9)的建模不确定性d(t)必然具有如下形式[9]
d(t)=ρ(x,t)|s|1/2sign(s)=ρ(x,t)ξ1    (16)
式中:ρ(x,t)是正的函数且满足
0<ρ(x,t)<δ    (17)
联立式(15)和(16)可得
&xi; &CenterDot; = A&xi;
A = 1 2 | &xi; 1 | &rho; ( x , t ) - &alpha; 1 - &beta; 0 - - - ( 18 )
定义Lyapunov函数
V = V 0 + 1 2 &kappa; 1 ( &alpha; - &alpha; 0 ) 2 + 1 2 &kappa; 2 ( &beta; - &beta; 0 ) 2 - - - ( 19 )
式中:α0和β0为正数,V0的定义如下
V 0 = &xi; T P&xi; , P = &lambda; + 4 &epsiv; 2 - 2 &epsiv; - 2 &epsiv; 1 - - - ( 20 )
式(20)中λ为任意正数;矩阵P是正定对称的,从而保证了函数V0为正。
先求函数V0对时间t的导数,并联立式(18)和(20)可得
V &CenterDot; 0 = &xi; &CenterDot; T P&xi; + &xi; T P &xi; &CenterDot; = &xi; T ( A T P + PA ) &xi; = - 1 2 | &xi; 1 | &xi; T &Lambda;&xi; - - - ( 21 )
式(21)中对称矩阵Λ定义如下
&Lambda; = &Lambda; 11 &Lambda; 12 &Lambda; 21 &Lambda; 22 - - - ( 22 )
式(22)中Λ11=-2(λ+4ε2)(ρ-α)-4εβ,Λ22=4ε,Λ12=Λ21=2ε(ρ-α)+β-λ-4ε2
取β=2εα,当α满足如下条件时,
&alpha; > &delta; ( &lambda; + 4 &epsiv; 2 ) - &epsiv; &lambda; ( 1 - &kappa; 1 ) + ( 2 &epsiv;&delta; - &lambda; - 4 &epsiv; 2 ) 2 12 &epsiv;&lambda; ( 1 - &kappa; 1 ) - - - ( 23 )
矩阵Λ的最小特征值λmin(Λ)≥2ε,即保证了其正定性。
因此,当α满足式(23)中的条件时,对于式(21)有
V &CenterDot; 0 &le; - &epsiv; | &xi; 1 | &xi; T &xi; = - &epsiv; | &xi; 1 | | | &xi; | | 2 - - - ( 24 )
由于有以下不等式成立,
λmin(P)||ξ||2≤ξTPξ≤λmax(P)||ξ||2    (25)
| &xi; 1 | = | s | 1 / 2 &le; | | &xi; | | 2 &le; V 0 1 / 2 &lambda; min 1 / 2 ( P ) - - - ( 26 )
式(26)中λmax(P)和λmin(P)分别为矩阵P的最大和最小特征值。因此对于式(24)有
V &CenterDot; 0 &le; - &zeta; V 0 1 / 2 , &zeta; = &epsiv; &lambda; min 1 / 2 ( P ) &lambda; max ( P ) - - - ( 27 )
求式(19)中函数V对时间t的导数,并结合式(27)可得
V &CenterDot; = V &CenterDot; 0 + 1 &kappa; 1 ( &alpha; - &alpha; 0 ) &alpha; &CenterDot; + 1 &kappa; 2 ( &beta; - &beta; 0 ) &beta; &CenterDot; &le; - &zeta; V 0 1 / 2 - &gamma; 1 2 &kappa; 1 | &alpha; - &alpha; 0 | - &gamma; 2 2 &kappa; 2 | &beta; - &beta; 0 | + 1 &kappa; 1 ( &alpha; - &alpha; 0 ) &alpha; &CenterDot; + 1 &kappa; 2 ( &beta; - &beta; 0 ) &beta; &CenterDot; + &gamma; 1 2 &kappa; 1 | &alpha; - &alpha; 0 | + &gamma; 2 2 &kappa; 2 | &beta; - &beta; 0 | - - - ( 28 )
利用不等式性质(a2+b2+c2)1/2≤|a|+|b|+|c|,可得
V &CenterDot; &le; - &sigma; V 1 / 2 + 1 &kappa; 1 ( &alpha; - &alpha; 0 ) &alpha; &CenterDot; + 1 &kappa; 2 ( &beta; - &beta; 0 ) &beta; &CenterDot; + &gamma; 1 2 &kappa; 1 | &alpha; - &alpha; 0 | + &gamma; 2 2 &kappa; 2 | &beta; - &beta; 0 | - - - ( 29 )
式(29)中σ=min{ζ,γ12}。
采用式(12)中的自适应律对增益α和β进行更新所得到的α和β都是有界的,这是由于当|s|>ν,0≤t≤tF时,
&alpha; = &alpha; ( 0 ) + &lambda; 1 &kappa; 1 2 &CenterDot; t - - - ( 30 )
tF是有限时间,因此α和β是有界的;而当|s|≤ν时,α和β都在递减直到|s|>ν重新满足。
基于对α和β的有界性分析可知,必定存在正数α0和β0使得α-α0<0和β-β0<0成立,因此式(29)可化为
V &CenterDot; &le; - &sigma; V - | &alpha; - &alpha; 0 | ( 1 &kappa; 1 &alpha; &CenterDot; - &gamma; 1 2 &kappa; 1 ) - | &beta; - &beta; 0 | ( 1 &kappa; 2 &beta; &CenterDot; - &gamma; 2 2 &kappa; 2 ) - - - ( 31 )
当|s|>ν时,由式(12)知
&alpha; &CenterDot; = &gamma; 1 &kappa; 1 2 - - - ( 32 )
γ2、к2为任意正数,则
&beta; &CenterDot; = 2 &epsiv; &alpha; = &CenterDot; &gamma; 2 &kappa; 2 2 - - - ( 33 )
将式(32)和(33)代入(31)中可得
V &CenterDot; &le; - &sigma; V 1 / 2 - - - ( 34 )
上式表明当滑模变量s远离零点即|s|>ν时,控制增益α和β分别以自适应律(32)和(33)递增直至条件(23)满足,使得矩阵Λ正定,从而保证滑模变量在有限时间tF内渐近收敛到零点附近,即|s|≤ν。
而当|s|≤ν时,控制增益α和β分别以如下的自适应律递减,
&alpha; &CenterDot; = - &gamma; 1 &kappa; 1 2 , &beta; &CenterDot; = - &gamma; 2 &kappa; 2 2 - - - ( 35 )
直至|s|>ν重新满足,如此循环往复。
因此,采用超螺旋滑模控制器(11)和增益自适应律(12)可使得滑模变量s在有限时间tF内渐近收敛到零附近任意小的范围内。由滑模变量s的定义可知系统输出的跟踪误差z1在有限时间tF内渐近收敛到零附近任意小的范围内,且该范围可通过参数进行调节。下面给出滑模变量s和系统输出的跟踪误差z1趋近于零点附近的有限时间tF的范围。
考虑|s|>ν时,增益α和β分别以自适应律(32)和(33)递增,条件(23)必在有限时间内满足。若ν=0,
对式(34)两边积分可知,在有限时间ts后,s→0,且
t s &le; 2 V 1 / 2 ( t 0 ) &sigma; - - - ( 36 )
则对于ν>0的情况,其有限时间tF一定满足
t F &le; t s &le; 2 V 1 / 2 ( t 0 ) &sigma; - - - ( 37 )
电液位置伺服系统的增益自调节的超螺旋滑模控制方法原理示意图如图2所示。
为考核所设计的控制器性能,采用式(5)描述的电液位置伺服系统模型进行仿真。电液位置伺服系统参数的真值如表1所示。
用于控制器设计的系统参数名义值则取为:转动惯量J=40kg·m2,粘性阻尼系数B=85N·m/(rad·s),马达排量Dm=3×10-4m3/rad,供油压力Ps=7×106Pa,回油压力Pr=0Pa,马达腔总控制容积Vt=2.2×10-3m3,泄漏系数Ct=1×10-11m5/(N·s),总流量系数kt=3.8×10-8m4/(s·V·N1/2),油液有效弹性模量βe=2.2×108Pa。系统参数的名义值与其真值之间的偏差归并到系统建模不确定性中。给定满足假设2的系统建模不确定性f(t)的取值为:f(t)=10sint·|s|1/2,期望跟踪的系统输入指令信号:x1d(t)=0.2sinπt·[1-exp(-0.01t3)],仿真采样时间为Ts=0.2ms。
表1电液位置伺服系统参数真值
为了便于分析本发明所提出控制器的有效性,分别使用本发明提出的具有自适应增益的超螺旋滑模控制(ASTSC)方法、传统的滑模控制(SMC)方法以及连续化近似的滑模控制(SSMC)方法进行仿真对比。三种控制器的参数选取如下:
1)ASTSC:控制器增益k1=5×104,k2=2×103;增益α和β的自适应律参数:γ1=150,к1=2,ε=1,ν=0.1;增益α和β的初值取为:α(0)=30,β(0)=0。
2)SMC:针对式(5)中的电液位置伺服系统,设计传统的滑模控制器如下:
u = 1 f 1 ( P L , u ) [ f 2 ( x 2 ) + f 3 ( x 3 ) + x &CenterDot; &CenterDot; &CenterDot; 1 d - k 1 z 2 - k 2 z 3 - k 3 sign ( s ) ] - - - ( 38 )
控制器参数取为:k1=5×104,k2=2×103,k3=200。
3)SSMC:运用饱和函数代替不连续的符号函数以使控制输入连续化,针对电液位置伺服系统(5)设计连续化近似的滑模控制器如下:
u = 1 f 1 ( P L , u ) [ f 2 ( x 2 ) + f 3 ( x 3 ) + x &CenterDot; &CenterDot; &CenterDot; 1 d - k 1 z 2 - k 2 z 3 - k 3 sat ( k 3 4 &chi; s ) ] - - - ( 39 )
式中:sat(·)为饱和函数。控制器参数取为:k1=5×104,k2=2×103,k3=200,χ=1。
ASTSC控制器作用下系统的实际位置输出跟踪输入指令信号的结果以及两种控制器跟踪误差对比曲线分别如图3和图4所示。由图4可以看出,在建模不确定性存在的情况下,ASTSC和SMC两种控制器都可以有效地抑制其对跟踪性能的影响,获得了很高的跟踪精度,而连续化的近似使得SSMC控制器的跟踪性能最差。由于选取的增益α和β的初值比较小使得ASTSC控制器的暂态性能不如SMC,但是在自适应律的作用下,增益α和β逐渐递增使ASTSC控制器的稳态性能明显优于SMC控制器。由于增益α和β呈比例关系,因此只给出增益α的自适应变化过程如图5所示。
尽管传统滑模控制器也获得了良好的跟踪精度,但是其付出的代价是控制输入发生如图6中所示的严重的颤振现象。由图7可见,在本发明提出的ASTSC控制器作用下,系统的控制输入信号是连续的,避免了SMC中的颤振问题。而且,值得注意的是,ASTSC的控制输入信号幅值比SMC小,却获得了更好的跟踪性能,这也说明了所提出的ASTSC控制方法更加高效。
综上所述,本发明公开的电液位置伺服系统的增益自调节的超螺旋滑模控制方法利用已知的系统模型信息,在传统超螺旋滑模控制算法中引入基于模型的前馈控制律,提升系统伺服精度。采用自适应律实时更新控制器增益,无需先验地获知系统建模不确定性的确切界,避免了传统算法中由人为设定与该界相关的控制增益造成的保守性。基于Lyapunov稳定性理论证明了闭环系统全局稳定,系统跟踪误差可在有限时间内渐近收敛到零附近任意小的范围内,且收敛的速度和稳态误差的界可通过参数进行调节。
虽然本发明已以较佳实施例揭露如上,然其并非用以限定本发明。本发明所属技术领域中具有通常知识者,在不脱离本发明的精神和范围内,当可作各种的更动与润饰。因此,本发明的保护范围当视权利要求书所界定者为准。

Claims (4)

1.一种电液位置伺服系统的增益自调节的超螺旋滑模控制方法,其特征在于,包括以下步骤:
步骤1,建立电液位置伺服系统的数学模型;
步骤2,设计具有自适应增益的超螺旋滑模控制器;
步骤3,具有自适应增益的超螺旋滑模控制器性能及稳定性测试。
2.根据权利要求1所述的电液位置伺服系统的增益自调节的超螺旋滑模控制方法,其特征在于,步骤1所述建立电液位置伺服系统的数学模型,具体如下:
(1.1)根据牛顿第二定律,电液位置伺服系统的惯性负载的运动方程为:
J &theta; &CenterDot; &CenterDot; m = P L D m - B &theta; &CenterDot; m - f ( t ) - - - ( 1 )
式(1)中J和θm分别为负载转动惯量和转角;PL=P1-P2为液压马达的负载压力,P1和P2分别为液压马达进油腔和回油腔压力;Dm为液压马达的体积排量;B为粘性摩擦系数;f(t)为系统建模不确定性,包含未建模的非线性摩擦、外负载干扰;
忽略液压马达的外泄漏,则液压马达负载压力动态方程为:
V t 4 &beta; e P &CenterDot; L = - D m &theta; &CenterDot; m - C t P L + Q L - - - ( 2 )
式(2)中Vt为马达两腔的总控制容积;βe为有效油液弹性模量;Ct为马达腔室的内泄漏系数;QL为负载流量;
将伺服阀的动态近似为比例环节,即伺服阀阀芯位移与控制输入成比例,因此负载流量可表示如下:
Q L = k t u P s - sign ( u ) P L - - - ( 3 )
式(3)中kt为总的流量增益;Ps为系统供油压力且sign(u)定义如下:
sign ( u ) = 1 , u &GreaterEqual; 0 - 1 , u < 0 - - - ( 4 )
(2.2)假设系统建模不确定性f(t)是连续可微的,则基于式(1)、(2)和(3),定义状态变量那么系统模型可以写成如下的状态空间的形式:
x &CenterDot; 1 = x 2
x &CenterDot; 2 = x 3 - - - ( 5 )
x &CenterDot; 3 = f 1 ( P L , u ) u - f 2 ( x 2 ) - f 3 ( x 3 ) + d ( t )
式(5)中
f 1 ( P L , u ) = 4 D m &beta; e k t J V t P s - sign ( u ) P L ,
f 2 ( x 2 ) = 4 D m &beta; e JV t ( D m + C t B D m ) x 2 ,               (6)
f 3 ( x 3 ) = 4 &beta; e JV t ( C t J + V t B 4 &beta; e ) x 3 ,
d ( t ) = - C t D m k t f ( t ) - V t 4 D m &beta; e k t f &CenterDot; ( t )
对于系统物理参数J、Dm、B、βe、Ct、kt和Ps,尽管无法获取其精确值,但是可假设其名义值已知并用于控制器的设计,而将参数名义值与其真值之间的偏差归并到建模不确定性d(t)中;
系统控制器的设计目标为:给定系统参考信号yd(t)=x1d(t),设计一个连续且有界的控制输入u使系统输出y=x1尽可能地跟踪系统的参考信号。
3.根据权利要求2所述的电液位置伺服系统的增益自调节的超螺旋滑模控制方法,其特征在于,步骤2所述设计具有自适应增益的超螺旋滑模控制器,步骤如下:
(2.1)定义系统误差变量如下:
z1=x1-x1d
z 2 = x 2 - x &CenterDot; 1 d - - - ( 7 )
z 3 = x 3 - x &CenterDot; &CenterDot; 1 d
式(7)中z1、z2和z3分别为惯性负载位置、速度和加速度跟踪误差;
定义如下的滑模变量:
s=k1z1+k2z2+z3             (8)
假设如下:
假设1:系统参考指令信号x1d(t)是三阶连续的,且系统期望位置指令、速度指令、加速度指令及加加速度指令都是有界的。液压马达位置伺服系统在一般工况下工作,即液压马达两腔压力P1,P2均小于供油压力Ps,且|PL|也小于Ps以保证式(6)中的f1>0;
假设2:系统建模不确定性d(t)满足以下条件
|d(t)|≤δ|s|1/2          (9)
式中:δ为未知的正数;
(3.2)对滑模变量求导可得:
s &CenterDot; = k 1 z 2 + k 2 z 3 + z &CenterDot; 3 = k 1 z 2 + k 2 z 3 + f 1 ( P L , u ) u - f 2 ( x 2 ) - f 3 ( x 3 ) + d ( t ) - x &CenterDot; &CenterDot; &CenterDot; 1 d - - - ( 10 )
基于式(10)中的滑模动态,设计超螺旋滑模控制器如下:
u = 1 f 1 ( P L , u ) ( u a + u s ) , u s = u s 1 + u s 2
u a = f 2 ( x 2 ) + f 3 ( x 3 ) + x &CenterDot; &CenterDot; &CenterDot; 1 d u s 1 = - k 1 z 2 - k 3 z 3 - &alpha; | s | 1 / 2 sign ( s ) - - - ( 11 )
u s 2 = - &Integral; 0 t &beta; 2 sign ( s ) d&tau;
式中:ua为用于改善模型补偿的前馈控制律;us为用于抑制建模不确定性d(t)对伺服系统性能影响的鲁棒控制律;为时变的控制器增益;
(3.3)设计如下的自适应律实时更新增益α和β:
&alpha; &CenterDot; = &gamma; 1 &kappa; 1 2 sign ( | s | - v ) - - - ( 12 )
β=2εα
式中:γ1、к1和ν都是任意的正数;ε为任意的实数。
4.根据权利要求3所述的电液位置伺服系统的增益自调节的超螺旋滑模控制方法,其特征在于,步骤3所述具有自适应增益的超螺旋滑模控制器性能及稳定性测试,具体如下:
将式(11)代入(10)中可得
s &CenterDot; = - &alpha; | s | 1 / 2 sign ( s ) + &eta; + d ( t )
&eta; &CenterDot; = - &beta; 2 sign ( s ) - - - ( 13 )
定义新的状态变量
ξ=[ξ12]T=[|s|1/2sign(s),η]T           (14)
结合式(13)可得
&xi; &CenterDot; 1 = 1 2 | &xi; 1 | ( - &alpha; &xi; 1 + &xi; 2 + d ( t ) ) &xi; &CenterDot; 2 = - &beta; 2 | &xi; 1 | &xi; 1 - - - ( 15 )
由假设2可知,满足条件(9)的建模不确定性d(t)必然具有如下形式[9]
d(t)=ρ(x,t)|s|1/2sign(s)=ρ(x,t)ξ1         (16)
式中:ρ(x,t)是正的函数且满足
0<ρ(x,t)<δ              (17)
联立式(15)和(16)可得
&xi; &CenterDot; = A&xi;
A = 1 2 | &xi; 1 | &rho; ( x , t ) - &alpha; 1 - &beta; 0 - - - ( 18 )
定义Lyapunov函数
V = V 0 + 1 2 &kappa; 1 ( &alpha; - &alpha; 0 ) 2 + 1 2 &kappa; 2 ( &beta; - &beta; 0 ) 2 - - - ( 19 )
式中:α0和β0为正数,V0的定义如下
V 0 = &xi; T P&xi; , P = &lambda; + 4 &epsiv; 2 - 2 &epsiv; - 2 &epsiv; 1 - - - ( 20 )
式(20)中λ为任意正数;矩阵P是正定对称的,从而保证了函数V0为正;
先求函数V0对时间t的导数,并联立式(18)和(20)可得
V &CenterDot; 0 = &xi; &CenterDot; T P&xi; + &xi; T P &xi; &CenterDot; = &xi; T ( A T P + PA ) &xi; = - 1 2 | &xi; 1 | &xi; T &Lambda;&xi; - - - ( 21 )
式(21)中对称矩阵Λ定义如下
&Lambda; = &Lambda; 11 &Lambda; 12 &Lambda; 21 &Lambda; 22 - - - ( 22 )
式(22)中Λ11=-2(λ+4ε2)(ρ-α)-4εβ,Λ22=4ε,Λ12=Λ21=2ε(ρ-α)+β-λ-4ε2
取β=2εα,当α满足如下条件时,
&alpha; > &delta; ( &lambda; + 4 &epsiv; 2 ) - &epsiv; &lambda; ( 1 - &kappa; 1 ) + ( 2 &epsiv;&delta; - &lambda; - 4 &epsiv; 2 ) 2 12 &epsiv;&lambda; ( 1 - &kappa; 1 ) - - - ( 23 )
矩阵Λ的最小特征值λmin(Λ)≥2ε,即保证了其正定性;
因此,当α满足式(23)中的条件时,对于式(21)有
V &CenterDot; 0 &le; - &epsiv; | &xi; 1 | &xi; T &xi; = - &epsiv; | &xi; 1 | | | &xi; | | 2 - - - ( 24 )
由于有以下不等式成立,
λmin(P)||ξ||2≤ξTPξ≤λmax(P)||ξ||2    (25)
| &xi; 1 | = | s | 1 / 2 &le; | | &xi; | | 2 &le; V 0 1 / 2 &lambda; min 1 / 2 ( P ) - - - ( 26 )
式(26)中λmax(P)和λmin(P)分别为矩阵P的最大和最小特征值。因此对于式(24)有
V &CenterDot; 0 &le; - &zeta; V 0 1 / 2 , &zeta; = &epsiv; &lambda; min 1 / 2 ( P ) &lambda; max ( P ) - - - ( 27 )
求式(19)中函数V对时间t的导数,并结合式(27)可得
V &CenterDot; = V &CenterDot; 0 + 1 &kappa; 1 ( &alpha; - &alpha; 0 ) &alpha; &CenterDot; + 1 &kappa; 2 ( &beta; - &beta; 0 ) &beta; &CenterDot; &le; - &zeta; V 0 1 / 2 - &gamma; 1 2 &kappa; 1 | &alpha; - &alpha; 0 | - &gamma; 2 2 &kappa; 2 | &beta; - &beta; 0 | + 1 &kappa; 1 ( &alpha; - &alpha; 0 ) &alpha; &CenterDot; + 1 &kappa; 2 ( &beta; - &beta; 0 ) &beta; &CenterDot; + &gamma; 1 2 &kappa; 1 | &alpha; - &alpha; 0 | + &gamma; 2 2 &kappa; 2 | &beta; - &beta; 0 | - - - ( 28 )
利用不等式性质(a2+b2+c2)1/2≤|a|+|b|+|c|,可得
V &CenterDot; &le; - &sigma; V 1 / 2 + 1 &kappa; 1 ( &alpha; - &alpha; 0 ) &alpha; &CenterDot; + 1 &kappa; 2 ( &beta; - &beta; 0 ) &beta; &CenterDot; + &gamma; 1 2 &kappa; 1 | &alpha; - &alpha; 0 | + &gamma; 2 2 &kappa; 2 | &beta; - &beta; 0 | - - - ( 29 )
式(29)中σ=min{ζ,γ12};
采用式(12)中的自适应律对增益α和β进行更新所得到的α和β都是有界的,这是由于当|s|>ν,0≤t≤tF时,
&alpha; = &alpha; ( 0 ) + &lambda; 1 &kappa; 1 2 &CenterDot; t - - - ( 30 )
tF是有限时间,因此α和β是有界的;而当|s|≤ν时,α和β都在递减直到|s|>ν重新满足;
基于对α和β的有界性分析可知,必定存在正数α0和β0使得α-α0<0和β-β0<0成立,因此式(29)可化为
V &CenterDot; &le; - &sigma; V - | &alpha; - &alpha; 0 | ( 1 &kappa; 1 &alpha; &CenterDot; - &gamma; 1 2 &kappa; 1 ) - | &beta; - &beta; 0 | ( 1 &kappa; 2 &beta; &CenterDot; - &gamma; 2 2 &kappa; 2 ) - - - ( 31 )
当|s|>ν时,由式(12)知
&alpha; &CenterDot; = &gamma; 1 &kappa; 1 2 - - - ( 32 )
γ2、к2为任意正数,则
&beta; &CenterDot; = 2 &epsiv; &alpha; &CenterDot; = &gamma; 2 &kappa; 2 2 - - - ( 33 )
将式(32)和(33)代入(31)中可得
V &CenterDot; &le; - &sigma; V 1 / 2 - - - ( 34 )
上式表明当滑模变量s远离零点即|s|>ν时,控制增益α和β分别以自适应律(32)和(33)递增直至条件(23)满足,使得矩阵Λ正定,从而保证滑模变量在有限时间tF内渐近收敛到零点附近,即|s|≤ν;
而当|s|≤ν时,控制增益α和β分别以如下的自适应律递减,
&alpha; &CenterDot; = - &gamma; 1 &kappa; 1 2 , &beta; &CenterDot; = - &gamma; 2 &kappa; 2 2 - - - ( 35 )
直至|s|>ν重新满足,如此循环往复;
因此,采用超螺旋滑模控制器(11)和增益自适应律(12)可使得滑模变量s在有限时间tF内渐近收敛到零附近任意小的范围内。由滑模变量s的定义可知系统输出的跟踪误差z1在有限时间tF内渐近收敛到零附近任意小的范围内,且该范围可通过参数进行调节。下面给出滑模变量s和系统输出的跟踪误差z1趋近于零点附近的有限时间tF的范围。
考虑|s|>ν时,增益α和β分别以自适应律(32)和(33)递增,条件(23)必在有限时间内满足;若ν=0,对式(34)两边积分可知,在有限时间ts后,s→0,且
t s &le; 2 V 1 / 2 ( t 0 ) &sigma; - - - ( 36 )
则对于ν>0的情况,其有限时间tF一定满足
t F &le; t s &le; 2 V 1 / 2 ( t 0 ) &sigma; - - - ( 37 )
CN201510261197.9A 2015-05-20 2015-05-20 一种电液位置伺服系统的增益自调节的超螺旋滑模控制方法 Pending CN104932259A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510261197.9A CN104932259A (zh) 2015-05-20 2015-05-20 一种电液位置伺服系统的增益自调节的超螺旋滑模控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510261197.9A CN104932259A (zh) 2015-05-20 2015-05-20 一种电液位置伺服系统的增益自调节的超螺旋滑模控制方法

Publications (1)

Publication Number Publication Date
CN104932259A true CN104932259A (zh) 2015-09-23

Family

ID=54119474

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510261197.9A Pending CN104932259A (zh) 2015-05-20 2015-05-20 一种电液位置伺服系统的增益自调节的超螺旋滑模控制方法

Country Status (1)

Country Link
CN (1) CN104932259A (zh)

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105388758A (zh) * 2015-11-10 2016-03-09 南京航空航天大学 一种基于Delta算子的液位控制系统的自适应滑模控制方法
CN105573125A (zh) * 2016-03-22 2016-05-11 西安理工大学 一种针对控制方向未知的气动系统位置自适应控制方法
CN106059413A (zh) * 2016-05-30 2016-10-26 中国人民解放军国防科学技术大学 一种直流电机驱动的飞轮系统转速控制方法
CN107165892A (zh) * 2017-06-14 2017-09-15 武汉科技大学 一种电液伺服系统的滑模控制方法
CN107215329A (zh) * 2017-05-10 2017-09-29 江苏大学 一种基于atsm的分布式驱动电动汽车横向稳定性控制方法
CN107861382A (zh) * 2017-09-28 2018-03-30 西北工业大学深圳研究院 一种多执行水下机器人鲁棒自适应运动控制装置及其方法
CN109372832A (zh) * 2018-12-21 2019-02-22 合肥工业大学 一种工况变化下的双变量液压系统能耗优化方法
CN110212834A (zh) * 2019-06-17 2019-09-06 北京润科通用技术有限公司 一种永磁同步电机矢量控制方法、装置及系统
CN111371357A (zh) * 2020-02-20 2020-07-03 江苏大学 一种基于自适应超螺旋算法的永磁同步电机调速控制方法
CN112720456A (zh) * 2020-12-02 2021-04-30 重庆大学 基于双时间尺度的柔性关节柔性杆机械臂控制方法及系统
CN113054883A (zh) * 2021-03-29 2021-06-29 上海海事大学 双凸极永磁电机控制方法
CN113110069A (zh) * 2021-05-24 2021-07-13 武汉大学 一种基于磁悬浮平面电机迭代神经网络鲁棒控制方法
CN113377029A (zh) * 2021-06-25 2021-09-10 中国民航大学 一种飞机舵机电动伺服系统多余力矩抑制方法
CN114035449A (zh) * 2021-11-11 2022-02-11 中国船舶重工集团公司第七0三研究所 一种机带海水冷却系统的多维联合仿真方法

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105388758B (zh) * 2015-11-10 2018-04-03 南京航空航天大学 一种基于Delta算子的液位控制系统的自适应滑模控制方法
CN105388758A (zh) * 2015-11-10 2016-03-09 南京航空航天大学 一种基于Delta算子的液位控制系统的自适应滑模控制方法
CN105573125A (zh) * 2016-03-22 2016-05-11 西安理工大学 一种针对控制方向未知的气动系统位置自适应控制方法
CN105573125B (zh) * 2016-03-22 2018-04-10 西安理工大学 一种针对控制方向未知的气动系统位置自适应控制方法
CN106059413A (zh) * 2016-05-30 2016-10-26 中国人民解放军国防科学技术大学 一种直流电机驱动的飞轮系统转速控制方法
CN107215329A (zh) * 2017-05-10 2017-09-29 江苏大学 一种基于atsm的分布式驱动电动汽车横向稳定性控制方法
CN107165892A (zh) * 2017-06-14 2017-09-15 武汉科技大学 一种电液伺服系统的滑模控制方法
CN107165892B (zh) * 2017-06-14 2018-05-29 武汉科技大学 一种电液伺服系统的滑模控制方法
CN107861382B (zh) * 2017-09-28 2019-12-13 西北工业大学深圳研究院 一种多执行水下机器人鲁棒自适应运动控制装置及其方法
CN107861382A (zh) * 2017-09-28 2018-03-30 西北工业大学深圳研究院 一种多执行水下机器人鲁棒自适应运动控制装置及其方法
CN109372832A (zh) * 2018-12-21 2019-02-22 合肥工业大学 一种工况变化下的双变量液压系统能耗优化方法
CN110212834A (zh) * 2019-06-17 2019-09-06 北京润科通用技术有限公司 一种永磁同步电机矢量控制方法、装置及系统
CN110212834B (zh) * 2019-06-17 2021-08-24 北京润科通用技术有限公司 一种永磁同步电机矢量控制方法、装置及系统
CN111371357A (zh) * 2020-02-20 2020-07-03 江苏大学 一种基于自适应超螺旋算法的永磁同步电机调速控制方法
CN112720456A (zh) * 2020-12-02 2021-04-30 重庆大学 基于双时间尺度的柔性关节柔性杆机械臂控制方法及系统
CN113054883A (zh) * 2021-03-29 2021-06-29 上海海事大学 双凸极永磁电机控制方法
CN113110069A (zh) * 2021-05-24 2021-07-13 武汉大学 一种基于磁悬浮平面电机迭代神经网络鲁棒控制方法
CN113377029A (zh) * 2021-06-25 2021-09-10 中国民航大学 一种飞机舵机电动伺服系统多余力矩抑制方法
CN113377029B (zh) * 2021-06-25 2023-02-24 中国民航大学 一种飞机舵机电动伺服系统多余力矩抑制方法
CN114035449A (zh) * 2021-11-11 2022-02-11 中国船舶重工集团公司第七0三研究所 一种机带海水冷却系统的多维联合仿真方法
CN114035449B (zh) * 2021-11-11 2024-02-20 中国船舶重工集团公司第七0三研究所 一种机带海水冷却系统的多维联合仿真方法

Similar Documents

Publication Publication Date Title
CN104932259A (zh) 一种电液位置伺服系统的增益自调节的超螺旋滑模控制方法
CN104345639B (zh) 一种电液位置伺服系统鲁棒自适应控制方法
CN105700347B (zh) 一种含磁滞补偿的液压马达预设性能跟踪控制方法
CN104111607B (zh) 一种考虑输入时滞的电机位置伺服系统的控制方法
Baghestan et al. An energy-saving nonlinear position control strategy for electro-hydraulic servo systems
CN105700352B (zh) 一种电液负载模拟器误差符号积分鲁棒控制方法
CN104317198A (zh) 具有时变输出约束的电液伺服系统非线性鲁棒位置控制方法
CN106483844B (zh) 基于非线性鲁棒的电液伺服系统自适应位置控制器的实现方法
CN104252134A (zh) 基于扩张状态观测器的电机伺服系统自适应鲁棒位置控制方法
CN105867133B (zh) 一种电液力矩伺服系统自调节误差符号积分鲁棒控制方法
He et al. Adaptive robust dead-zone compensation control of electro-hydraulic servo systems with load disturbance rejection
CN104635490A (zh) 一种单出杆液压缸位置伺服系统的输出反馈控制方法
CN104345638A (zh) 一种液压马达位置伺服系统的自抗扰自适应控制方法
CN105843043B (zh) 一种电液负载模拟装置自适应鲁棒力控制方法
CN104698844A (zh) 液压位置伺服系统的不确定性补偿的滑模控制方法
CN104238361A (zh) 电机伺服系统自适应鲁棒位置控制方法与系统
CN108181818A (zh) 含未建模摩擦动态的电液位置伺服系统鲁棒自适应控制方法
CN106154833A (zh) 一种电液负载模拟装置输出反馈控制方法
CN105003495A (zh) 一种非对称液压缸实时动态补偿加载装置
Hasan Disturbance attenuation of n+ 1 coupled hyperbolic PDEs
CN104360596A (zh) 机电伺服系统有限时间摩擦参数辨识和自适应滑模控制方法
CN110308651A (zh) 基于扩张状态观测器的电液伺服系统全状态约束控制方法
CN106125553A (zh) 一种考虑状态约束的液压系统自适应控制方法
Zad et al. Robust Model Predictive position Control of direct drive electro-hydraulic servo system
CN105068426A (zh) 基于干扰补偿的电液位置伺服系统连续滑模控制方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20150923