CN105068426B - 基于干扰补偿的电液位置伺服系统连续滑模控制方法 - Google Patents

基于干扰补偿的电液位置伺服系统连续滑模控制方法 Download PDF

Info

Publication number
CN105068426B
CN105068426B CN201510522508.2A CN201510522508A CN105068426B CN 105068426 B CN105068426 B CN 105068426B CN 201510522508 A CN201510522508 A CN 201510522508A CN 105068426 B CN105068426 B CN 105068426B
Authority
CN
China
Prior art keywords
msub
mrow
mover
centerdot
mtd
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
CN201510522508.2A
Other languages
English (en)
Other versions
CN105068426A (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.)
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 CN201510522508.2A priority Critical patent/CN105068426B/zh
Publication of CN105068426A publication Critical patent/CN105068426A/zh
Application granted granted Critical
Publication of CN105068426B publication Critical patent/CN105068426B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Steroid Compounds (AREA)
  • Feedback Control In General (AREA)

Abstract

本发明公开了一种基于干扰补偿的电液位置伺服系统连续滑模控制方法,方法包括:建立电液位置伺服系统的数学模型;设计电液位置伺服系统两个通道的干扰观测器;设计基于干扰观测器的连续滑模控制器。本发明消除滑模控制方法中的不连续项的同时保证了该方法的强鲁棒性,使得滑模控制器的输出连续化,彻底消除了滑模控制的抖振问题;在控制器中不使用系统的加速度信号,削弱了测量噪声对跟踪性能的恶化,利于在工程实际中运用;在系统同时存在匹配和不匹配干扰和强非线性的情况下使用滑模控制方法依旧获得了渐近跟踪的稳态性能,保证了电液位置伺服系统良好的控制性能。

Description

基于干扰补偿的电液位置伺服系统连续滑模控制方法
技术领域
本发明属于电液伺服控制技术领域,特别是一种基于干扰补偿的电液位置伺服系统连续滑模控制方法。
背景技术
液压伺服系统具有功重比大、响应快及抗负载刚性强等突出优点,在众多重要领域内得到广泛运用。电液伺服系统是一个典型的非线性系统,包含许多非线性特性和建模不确定性干扰。随着电液伺服系统向高精度、高频响发展时,系统呈现的非线性特性对系统性能的影响越显著,而且干扰的存在会使以系统名义模型设计的控制器不稳定或降阶,因此电液伺服系统非线性特性和干扰是限制系统性能提升的重要因素。随着工业及国防领域技术水平的不断进步,以往基于传统线性理论设计的控制器已逐渐不能满足系统的高性能需求,因此必须针对电液伺服系统中的非线性特性和存在的干扰研究更加先进的非线性控制策略。
针对电液伺服系统的干扰和非线性控制的问题,许多方法相继被提出。在液压位置伺服系统控制器的设计中,针对电液伺服系统存在的多种干扰,反馈线性化控制的基本思想是通过在控制器中对非线性函数进行精确补偿以使误差动态线性化。虽然理论上可以获得完美的渐近跟踪性能,但是实际系统的模型是不可能精确已知的,总会存在建模误差,因此会恶化理论分析获得的跟踪性能。而由于系统中存在不匹配干扰,传统的滑模控制方法的基本思路是通过增大控制器的鲁棒性来克服不匹配和匹配干扰从而到达滑模面,但是,即使到达滑模面后,系统地跟踪误差在不匹配干扰的影响下仍然无法为零,只能得到一个和不匹配干扰上确界相关的一个一致有界的稳态跟踪误差。并且,通过增大不连续项增益的方法来增加控制器的鲁棒性,在实际运用中很可能激发系统高频动态,使系统失稳。同时,由于滑模控制方法在进行控制器设计中使用了符号函数,从而使得控制器的输出不连续,虽然有的滑模控制策略补偿了系统存在的干扰使得不连续项增益大为减少,削弱滑模控制的抖振,但因为控制器表达式中始终存在不连续项,从而使得滑模控制方法始终存在抖振,无法根除抖振,因而传统的滑模控制方法具有很大的工程局限性。
发明内容
本发明的目的在于提供一种基于干扰补偿的电液位置伺服系统连续滑模控制方法。
实现本发明目的的技术解决方案为:一种基于干扰补偿的电液位置伺服系统连续滑模控制方法,包括以下步骤:
步骤1、建立电液位置伺服系统的数学模型;
步骤2、设计电液位置伺服系统两通道的干扰观测器;
步骤3、设计基于干扰观测器的连续滑模控制器。
本发明与现有技术相比,其显著优点是:
(1)本发明消除滑模控制策略中的不连续项的同时保证了该方法的强鲁棒性,使得滑模控制器的输出连续化,彻底消除了滑模控制的抖振问题;
(2)在系统同时存在匹配和不匹配干扰和强非线性的情况下使用滑模控制方法依旧获得了渐近跟踪的稳态性能,保证了电液位置伺服系统良好的控制性能;
(3)本发明设计干扰观测器观测了液压位置伺服系统第二通道和第三通道的干扰,并在控制器中将其完全补偿,抑制了干扰对控制性能的影响;
(4)在控制器中不使用系统的加速度信号,削弱了测量噪声对跟踪性能的恶化,利于在工程实际中运用。
附图说明
图1为本发明的基于干扰补偿的电液位置伺服系统连续滑模控制方法流程图。
图2为本发明的液压位置伺服系统的原理图。
图3为本发明的基于干扰补偿的电液位置伺服系统连续滑模控制方法原理示意图。
图4为本发明实施例中系统输出对期望指令的跟踪曲线图。
图5为本发明实施例中系统的位置跟踪误差随时间变化的曲线图。
图6为本发明实施例中基于干扰补偿的电液位置伺服系统连续滑模控制器(C‐SMC)作用下和不确定性补偿的滑模控制器(UC‐SMC)作用下及反馈线性化控制器(FLC)作用下的位置跟踪误差随时间变化曲线图。
图7为本发明实施例中基于干扰补偿的电液位置伺服系统连续滑模控制器(C‐SMC)作用下系统的控制输入随时间变化的曲线图。
图8为本发明实施例中不确定性补偿的滑模控制器(UC‐SMC)作用下系统的控制输入随时间变化的曲线图。
图9为本发明实施例中基于干扰补偿的电液位置伺服系统连续滑模控制器作用下系统的第二通道干扰观测曲线图。
图10为本发明实施例中基于干扰补偿的电液位置伺服系统连续滑模控制器作用下系统的第二通道干扰观测误差随时间变化曲线图。
图11为本发明实施例中基于干扰补偿的电液位置伺服系统连续滑模控制器作用下系统的第三通道干扰观测曲线图。
图12为本发明实施例中基于干扰补偿的电液位置伺服系统连续滑模控制器作用下系统的第三通道干扰观测误差随时间变化曲线图。
具体实施方式
下面结合附图及具体实施例对本发明作进一步详细说明。
结合图1,本发明的液压位置伺服系统的不确定性补偿的滑模控制方法,包括以下步骤:
步骤1、建立液压位置伺服系统的数学模型;
步骤1‐1、如图2右半部分所示,液压位置伺服系统为通过伺服阀控制的液压马达驱动惯性负载的系统构成,图2左半部分为伺服阀控制的液压马达的原理示意图;根据牛顿第二定律,惯性负载的运动方程为:
式(1)中m为惯性负载参数;PL为液压马达两腔压差;A为液压马达的排量;B为粘性摩擦系数;为建模误差,包括m、PL、B的名义值与真实值之间的偏差以及外负载干扰;y为惯性负载的位移;为惯性负载的速度,为惯性负载的加速度;t为时间变量;
液压马达左右两腔的压力动态方程为:
式(2)中P1和P2分别为液压马达两腔的压力,分别为P1和P2的导数;V1=V01+Ay,V2=V02-Ay,V1和V2分别表示液压马达两腔的控制容积;V01和V02分别为液压马达两腔的初始容积;βe为有效油液弹性模量;Ct为内泄漏系数;q1(t)和q2(t)分别为P1和P2动态方程的建模误差;Q1和Q2分别为液压马达的进油腔流量和回油腔流量;Q1和Q2与伺服阀位移xv的关系为:
式(3)中s(xv)的定义为:
其中,kq为流量增益,Cd流量系数;ω为阀芯面积梯度;ρ为油液密度;Ps为供油压力,Pr为回油压力;液压马达两腔压力满足0<Pr<P1<Ps,0<Pr<P2<Ps,|PL|<<PS
由于考虑伺服阀动态需要安装额外的位移传感器来获取伺服阀阀芯的位移,而且对于跟踪性能只有微小的提升;因此大量相关的研究都忽略伺服阀的动态,假设采用的是高响应的伺服阀,阀芯位移与控制输入近似为比例环节即xv=kiu,故式(3)可以写成
式(5)中kt=kqki代表总的流量增益,ki为伺服阀增益,u为液压位置伺服系统输入,
步骤1‐2、定义状态变量:则系统的状态方程为:
式(6)中m、B、βe、A、kt、V01、V02和Ct在观测器和控制器的设计中为名义值,其与真实值之间的偏差集中归类到系统干扰中处理,在第二通中是d(x,t),在第三通道中是q(t);其中:
由于液压系统参数m,B,βe,kt和Ct受各种因素(如温度、组件磨损程度等)影响变化很大,因此为了简化系统状态方程,定义:
因为|PL|<<PS,从而g(x)≠0;第二通道干扰d1(x,t)和第三通道干扰d2(x,t)都是有界的,即:|d1(x,t)|≤D1,|d2(x,t)|≤D2,其中D1、D2分别为|d1(x,t)|和|d1(x,t)|上界,都是已知正数,并且d1(x,t)一阶导数存在;则液压位置伺服系统模型为
步骤2,设计两个通道的干扰观测器,步骤如下:
步骤2‐1、设计第二通道干扰观测器:
定义第二通道干扰观测器滑模面s1为:
s1=z1-x2 (10)
其中,z1为第二通道干扰观测器内动态;
式(11)中,k1、β1、ε1、p1和q1均为第二通道干扰观测器参数;p1<q1,均为正奇数,k1、β1、ε1均为正数,β1≥D1
由式(10)、(11)有:
故第二通道干扰d1(x,t)的观测值设计为其表达式如下:
定义第二通道干扰观测器的Lyapunov方程:
由于β1≥D1,故:
又因,若存在一正定函数V0(t)满足下面不等式:
则有,V0(t)在时间ts内收敛到平衡点,其中
式(18)中,α>0,λ>0,0<γ<1;
由式(17)、(18)式有,V1(t)将在有限时间内收敛到平衡点,即存在一个时刻t2,在经过有限时间t2后,s1为零,此时也将收敛到零,而d1(x,t)的观测误差为:
则干扰的观测误差也将在有限时间t2内为0;即经过时间t2
综上,得到第二通道干扰观测器表达式为:
步骤2‐2、设计第三通道干扰观测器:
定义第三通道干扰观测器滑模面s2为:
s2=z2-x3 (20)
其中,z2为第三通道干扰观测器内部动态;
式(21)中,k2、β2、ε2、p2和q2均为第三通道干扰观测器参数;其中p2<q2,均为正奇数,k2、β2、ε2均为正数,β2≥D2
由式(20)、(21)有:
故第三通道干扰d2(x,t)的观测值设计为其表达式如下:
定义第三通道干扰观测器的Lyapunov方程:
又因β2>D2,则有,
此时,由式(17)、(18)式有,V2(t)将在有限时间内为零,即存在一个时刻t3,在经过t3后,s2为零,此时也将收敛到零,又因d2(x,t)估计误差
则干扰的估计误差也将在有限时间t3内为零,即经过时间t3后有:
综上,得到第三通道干扰观测器表达式为:
步骤3、设计基于干扰观测器的连续滑模控制器,具体如下:
定义液压伺服系统的位置跟踪误差e0(t)、速度的跟踪误差e1(t)、加速度的跟踪误差e2(t)、加速度导数的跟踪误差e3(t):
e0(t)=x1-xd(t) (27)
其中,xd(t)为系统参考信号,xd(t)是三阶连续可导的,且系统参考位置信号xd(t)、系统参考速度信号系统参考加速度信号及系统参考加速度导数的信号都是有界的;
定义连续滑模控制器滑模面s为:
其中c1、c2、c3为滑模控制器参数,且均大于零,并且使得表达式是Hurwitz的,则有:
设计连续滑模控制器u的表达式为:
其中k3、k4、ζ为控制器参数,且k3>0、k4>0、0<ζ<1。
步骤4、连续滑模控制器的稳定性测试;具体如下:
步骤4‐1、将式(33)代入式(32)有:
定义滑模控制器Lyapunov方程:
则有:
故,由式(17)、(18)可知,V(t)将在有限时间内为零,即存在一个时刻t1,在经过有限时间t1后,s为零,即:
又因s1、s2也是有限时间内为0,t1为s为零的时刻,t2为s1为零的时刻,t3为s2为零的时刻,则存在t4=max{t1,t2,t3},经过t4时刻后有:
又因:
则:
e2(t)+c1e1(t)+c2e0(t)+c3∫e0(t)dt=d1(t) (40)即:
步骤4‐2、当时:
则有,当t→∞有:
故e0(t)=x1-xd(t)在时间趋于无穷的条件下趋于零;
δ为一已知正数时:
则有,当t→∞有:
故e0(t)=x1-xd(t)在时间趋于无穷的条件下收敛于一个一致稳定界内。
综上可知,针对液压位置伺服系统设计的连续滑模控制器可以使系统得到全局渐近稳定的结果。调节观测器系数k1、β1、ε1、p1、q1、k2、β2、ε2、p2、q2可以使观测器的跟踪误差在有限时间内趋于零,调节增益c1、c2、c3、k3、k4、ζ可以使系统的跟踪误差在时间趋于无穷的条件下趋于零。液压位置伺服系统的基于干扰补偿的电液位置伺服系统连续滑模控制方法原理示意图如图3所示。
下面结合具体实施例对本发明做进一步说明。
实施例1
结合图1‐图3,为考核所设计的控制器性能,在MATLAB/Simulink中搭建系统仿真模型,在仿真中取如下参数,仿真步长设置为0.0002s:
负载转动惯量m=40kg·m2,马达排量A=2×10-4m3/rad,粘性摩擦系数B=80N·m·s/rad,供油压力Ps=7MPa,回油压力Pr=0,油液弹性模量βe=2×108Pa,马达两腔初始容积V01=V02=1×10-3m3,泄漏系数Ct=9×10-12m3/s/Pa,总流量增益压力动态建模误差q1(t)=q2(t)=6×10-6m3·rad/s,外负载干扰f(t)=200[1-exp(-0.1t3)](N·m),添加未建模摩擦项Fc(x2)=20/π[arctan(100x2)]N·m。
给定系统的期望指令为xd(t)=arctan(sin(t))[1‐exp(‐0.01t3)]rad。
取如下的控制器以作对比:
连续滑模控制器(Continuous‐Sliding Mode Controller,C‐SMC):取第二通道干扰观测器参数k1=10000,β1=250,ε1=0.06,p1=3,q1=5,第三通道干扰观测器参数k2=10000,β2=60,ε2=0.05,p2=5,q2=7,控制器参数c1=24,c2=192,c3=512,k3=3、k4=5、ζ=0.5。
匹配和不匹配不确定性补偿的滑模控制器(Uncertainties Compensation‐Sliding Mode Controller,UC‐SMC):不匹配不确定性观测器参数k1=10000,β1=250,ε1=0.06,p1=3,q1=5,第三通道干扰观测器参数k2=10000,β2=60,ε2=0.05,p2=5,q2=7,控制器参数c1=24,c2=192,c3=512,k=2。
反馈线性化控制器(Feedback Linearization Controller,FLC):l1=20、l2=20、l3=20。
连续滑模控制器作用下电液位置伺服系统的位置输出对期望指令的跟踪曲线如图4所示,跟踪误差如图5所示;由图4可知,系统的位置输出曲线基本与期望指令曲线重合,说明连续滑模控制器具有较好的控制精度,在该控制器的作用下,系统表现出良好的跟踪性能;由图5可知,在连续滑模控制器作用下,液压伺服系统的位置输跟踪精度很高,稳态跟踪误差收敛到幅值为2×10-3(rad)的界内。
图6为连续滑模控制器(Continuous‐Sliding Mode Controller,C‐SMC)、匹配和不匹配不确定性补偿的滑模控制器(Uncertainties Compensation‐Sliding ModeController,UC‐SMC)、反馈线性化控制器(Feedback Linearization Controller,FLC)分别作用下系统的跟踪误差曲线;由图6可知,在系统存在未建模干扰时,反馈线性化控制器的控制性能比连续滑模控制器和不确定性补偿的滑模控制器差,因为电液位置伺服系统的模型不可能精确已知,总会存在未建模的干扰,而存在的未建模干扰将严重恶化了反馈线性化控制器的控制性能。而连续滑模控制器和不确定性补偿的滑模控制器因为滑模面参数选取是一样的,故两者具有相同的误差动态曲线,同时说明连续滑模控制器和不确定性补偿的滑模控制器具有一样的强鲁棒性。
图7为连续滑模控制器的控制输入曲线,图8为不确定性补偿的滑模控制器的输入曲线;由图7,图8及其对应的局部放大图对比可知,连续滑模控制器的控制器输入曲线为一低频连续曲线,便于在工程实际中运用,而不确定性补偿的滑模控制器虽然控制输入曲线的抖振不大,但是依然存在高频抖振,这在实际的工程运用中,易激发系统潜在的未建模高频动态,导致系统发散失稳。同时结合图6可知,虽然连续滑模控制器的控制输入曲线不存在高频抖振,但是仍然和不确定性补偿的滑模控制器具有一样的强鲁棒性。
图9为第二通道干扰和第二通道干扰观测曲线,图10为第二通道干扰和第二通道干扰观测误差随时间变化曲线,由图9可知,实际存在的第二通道干扰曲线和其观测曲线基本重合,所设计的干扰观测器对液压位置伺服系统的第二通道干扰估计非常准确,同时结合图10可知,第二通道干扰的估计误差在经过很短的时间后迅速收敛到零。
图11为第三通道干扰和第三通道干扰观测曲线,图12为第三通道干扰和第三通道干扰观测误差随时间变化曲线;在图11中,实际存在的第三通道干扰曲线和其估计曲线基本重合,而由图12可知,其估计误差在有限的时间内迅速收敛到零,说明所设计的干扰观测器对系统存在的匹配干扰具有较好的估计能力。

Claims (2)

1.一种基于干扰补偿的电液位置伺服系统连续滑模控制方法,其特征在于,包括以下步骤:
步骤1、建立电液位置伺服系统的数学模型;具体如下:
步骤1-1、电液位置伺服系统为通过伺服阀控制的液压马达驱动惯性负载的系统;根据牛顿第二定律,惯性负载的运动方程为:
<mrow> <mi>m</mi> <mover> <mi>y</mi> <mo>&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mo>=</mo> <msub> <mi>P</mi> <mi>L</mi> </msub> <mi>A</mi> <mo>-</mo> <mi>B</mi> <mover> <mi>y</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>+</mo> <mi>f</mi> <mrow> <mo>(</mo> <mi>y</mi> <mo>,</mo> <mover> <mi>y</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
式(1)中m为惯性负载参数;PL为液压马达两腔压差;A为液压马达的排量;B为粘性摩擦系数;为建模误差,包括m、PL、B的名义值与真实值之间的偏差以及外负载干扰;y为惯性负载的位移;为惯性负载的速度,为惯性负载的加速度;t为时间变量;
液压马达左右两腔的压力动态方程为:
<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mover> <mi>P</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>1</mn> </msub> <mo>=</mo> <mfrac> <msub> <mi>&amp;beta;</mi> <mi>e</mi> </msub> <msub> <mi>V</mi> <mn>1</mn> </msub> </mfrac> <mo>&amp;lsqb;</mo> <mo>-</mo> <mi>A</mi> <mover> <mi>y</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>-</mo> <msub> <mi>C</mi> <mi>t</mi> </msub> <msub> <mi>P</mi> <mi>L</mi> </msub> <mo>+</mo> <msub> <mi>q</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>Q</mi> <mn>1</mn> </msub> <mo>&amp;rsqb;</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mover> <mi>P</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>2</mn> </msub> <mo>=</mo> <mfrac> <msub> <mi>&amp;beta;</mi> <mi>e</mi> </msub> <msub> <mi>V</mi> <mn>2</mn> </msub> </mfrac> <mo>&amp;lsqb;</mo> <mi>A</mi> <mover> <mi>y</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>+</mo> <msub> <mi>C</mi> <mi>t</mi> </msub> <msub> <mi>P</mi> <mi>L</mi> </msub> <mo>-</mo> <msub> <mi>q</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>Q</mi> <mn>2</mn> </msub> <mo>&amp;rsqb;</mo> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
式(2)中P1和P2分别为液压马达两腔的压力,分别为P1和P2的导数;V1=V01+Ay,V2=V02-Ay,V1和V2分别表示液压马达两腔的控制容积;V01和V02分别为液压马达两腔的初始容积;βe为有效油液弹性模量;Ct为内泄漏系数;q1(t)和q2(t)分别为P1和P2动态方程的建模误差;Q1和Q2分别为液压马达的进油腔流量和回油腔流量;Q1和Q2与伺服阀位移xv的关系为:
<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mi>Q</mi> <mn>1</mn> </msub> <mo>=</mo> <msub> <mi>k</mi> <mi>q</mi> </msub> <msub> <mi>x</mi> <mi>v</mi> </msub> <mo>&amp;lsqb;</mo> <mi>s</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>v</mi> </msub> <mo>)</mo> </mrow> <msqrt> <mrow> <msub> <mi>P</mi> <mi>s</mi> </msub> <mo>-</mo> <msub> <mi>P</mi> <mn>1</mn> </msub> </mrow> </msqrt> <mo>+</mo> <mi>s</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>x</mi> <mi>v</mi> </msub> <mo>)</mo> </mrow> <msqrt> <mrow> <msub> <mi>P</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>P</mi> <mi>r</mi> </msub> </mrow> </msqrt> <mo>&amp;rsqb;</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>Q</mi> <mn>2</mn> </msub> <mo>=</mo> <msub> <mi>k</mi> <mi>q</mi> </msub> <msub> <mi>x</mi> <mi>v</mi> </msub> <mo>&amp;lsqb;</mo> <mi>s</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>v</mi> </msub> <mo>)</mo> </mrow> <msqrt> <mrow> <msub> <mi>P</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>P</mi> <mi>r</mi> </msub> </mrow> </msqrt> <mo>+</mo> <mi>s</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>x</mi> <mi>v</mi> </msub> <mo>)</mo> </mrow> <msqrt> <mrow> <msub> <mi>P</mi> <mi>s</mi> </msub> <mo>-</mo> <msub> <mi>P</mi> <mn>1</mn> </msub> </mrow> </msqrt> <mo>&amp;rsqb;</mo> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
式(3)中s(xv)的定义为:
<mrow> <mi>s</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>v</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mn>1</mn> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <mi>i</mi> <mi>f</mi> <mi> </mi> <msub> <mi>x</mi> <mi>v</mi> </msub> <mo>&amp;GreaterEqual;</mo> <mn>0</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mn>0</mn> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <mi>i</mi> <mi>f</mi> <mi> </mi> <msub> <mi>x</mi> <mi>v</mi> </msub> <mo>&lt;</mo> <mn>0</mn> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>
其中,kq为流量增益,Cd流量系数;ω为阀芯面积梯度;ρ为油液密度;Ps为供油压力,Pr为回油压力;液压马达两腔压力满足0<Pr<P1<Ps,0<Pr<P2<Ps,|PL|<<PS
阀芯位移与控制输入近似为比例环节即xv=kiu,故式(3)写成:
<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mi>Q</mi> <mn>1</mn> </msub> <mo>=</mo> <msub> <mi>k</mi> <mi>t</mi> </msub> <mi>u</mi> <mo>&amp;lsqb;</mo> <mi>s</mi> <mrow> <mo>(</mo> <mi>u</mi> <mo>)</mo> </mrow> <msqrt> <mrow> <msub> <mi>P</mi> <mi>s</mi> </msub> <mo>-</mo> <msub> <mi>P</mi> <mn>1</mn> </msub> </mrow> </msqrt> <mo>+</mo> <mi>s</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>u</mi> <mo>)</mo> </mrow> <msqrt> <mrow> <msub> <mi>P</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>P</mi> <mi>r</mi> </msub> </mrow> </msqrt> <mo>&amp;rsqb;</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>Q</mi> <mn>2</mn> </msub> <mo>=</mo> <msub> <mi>k</mi> <mi>t</mi> </msub> <mi>u</mi> <mo>&amp;lsqb;</mo> <mi>s</mi> <mrow> <mo>(</mo> <mi>u</mi> <mo>)</mo> </mrow> <msqrt> <mrow> <msub> <mi>P</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>P</mi> <mi>r</mi> </msub> </mrow> </msqrt> <mo>+</mo> <mi>s</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>u</mi> <mo>)</mo> </mrow> <msqrt> <mrow> <msub> <mi>P</mi> <mi>s</mi> </msub> <mo>-</mo> <msub> <mi>P</mi> <mn>1</mn> </msub> </mrow> </msqrt> <mo>&amp;rsqb;</mo> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>
式(5)中kt=kqki代表总的流量增益,ki为伺服阀增益,u为液压位置伺服系统输入,
步骤1-2、定义状态变量:则系统的状态方程为:
<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mover> <mi>x</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>1</mn> </msub> <mo>=</mo> <msub> <mi>x</mi> <mn>2</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mover> <mi>x</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>2</mn> </msub> <mo>=</mo> <msub> <mi>x</mi> <mn>3</mn> </msub> <mo>-</mo> <msub> <mi>bx</mi> <mn>2</mn> </msub> <mo>+</mo> <mi>d</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mover> <mi>x</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>3</mn> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mi>A&amp;beta;</mi> <mi>e</mi> </msub> <msub> <mi>k</mi> <mi>t</mi> </msub> </mrow> <mi>m</mi> </mfrac> <mrow> <mo>(</mo> <mfrac> <msub> <mi>R</mi> <mn>1</mn> </msub> <msub> <mi>V</mi> <mn>1</mn> </msub> </mfrac> <mo>+</mo> <mfrac> <msub> <mi>R</mi> <mn>2</mn> </msub> <msub> <mi>V</mi> <mn>2</mn> </msub> </mfrac> <mo>)</mo> </mrow> <mi>u</mi> <mo>-</mo> <mfrac> <mrow> <msup> <mi>A</mi> <mn>2</mn> </msup> <msub> <mi>&amp;beta;</mi> <mi>e</mi> </msub> </mrow> <mi>m</mi> </mfrac> <mrow> <mo>(</mo> <mfrac> <mn>1</mn> <msub> <mi>V</mi> <mn>1</mn> </msub> </mfrac> <mo>+</mo> <mfrac> <mn>1</mn> <msub> <mi>V</mi> <mn>2</mn> </msub> </mfrac> <mo>)</mo> </mrow> <msub> <mi>x</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>&amp;beta;</mi> <mi>e</mi> </msub> <msub> <mi>C</mi> <mi>t</mi> </msub> <mrow> <mo>(</mo> <mfrac> <mn>1</mn> <msub> <mi>V</mi> <mn>1</mn> </msub> </mfrac> <mo>+</mo> <mfrac> <mn>1</mn> <msub> <mi>V</mi> <mn>2</mn> </msub> </mfrac> <mo>)</mo> </mrow> <msub> <mi>x</mi> <mn>3</mn> </msub> <mo>+</mo> <mi>q</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>
式(6)中m、B、A、βe、kt、V01、V02和Ct在观测器和控制器的设计中为名义值,其与真实值之间的偏差集中归类到系统干扰中处理,在第二通道中是d(x,t),在第三通道中是q(t);其中:
<mrow> <mtable> <mtr> <mtd> <mrow> <mi>b</mi> <mo>=</mo> <mfrac> <mi>B</mi> <mi>m</mi> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>d</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>f</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>/</mo> <mi>m</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>q</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msub> <mi>A&amp;beta;</mi> <mi>e</mi> </msub> </mrow> <mi>m</mi> </mfrac> <mrow> <mo>(</mo> <mfrac> <mrow> <msub> <mi>q</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <msub> <mi>V</mi> <mn>1</mn> </msub> </mfrac> <mo>+</mo> <mfrac> <mrow> <msub> <mi>q</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <msub> <mi>V</mi> <mn>2</mn> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>R</mi> <mn>1</mn> </msub> <mo>=</mo> <mi>s</mi> <mrow> <mo>(</mo> <mi>u</mi> <mo>)</mo> </mrow> <msqrt> <mrow> <msub> <mi>P</mi> <mi>s</mi> </msub> <mo>-</mo> <msub> <mi>P</mi> <mn>1</mn> </msub> </mrow> </msqrt> <mo>+</mo> <mi>s</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>u</mi> <mo>)</mo> </mrow> <msqrt> <mrow> <msub> <mi>P</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>P</mi> <mi>r</mi> </msub> </mrow> </msqrt> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>R</mi> <mn>2</mn> </msub> <mo>=</mo> <mi>s</mi> <mrow> <mo>(</mo> <mi>u</mi> <mo>)</mo> </mrow> <msqrt> <mrow> <msub> <mi>P</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>P</mi> <mi>r</mi> </msub> </mrow> </msqrt> <mo>+</mo> <mi>s</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>u</mi> <mo>)</mo> </mrow> <msqrt> <mrow> <msub> <mi>P</mi> <mi>s</mi> </msub> <mo>-</mo> <msub> <mi>P</mi> <mn>2</mn> </msub> </mrow> </msqrt> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>
为简化系统状态方程,定义:
因为|PL|<<PS,则g(x)≠0;第二通道干扰d1(x,t)和第三通道干扰d2(x,t)都是有界的,即:|d1(x,t)|≤D1,|d2(x,t)|≤D2,其中D1、D2分别为|d1(x,t)|和|d1(x,t)|上界,都是已知正数,并且d1(x,t)一阶导数存在;则液压位置伺服系统模型为:
步骤2、设计电液位置伺服系统两通道的干扰观测器;步骤如下:
步骤2-1、设计第二通道干扰观测器:
定义第二通道干扰观测器滑模面s1为:
s1=z1-x2 (10)
其中,z1为第二通道干扰观测器内部动态;
<mrow> <msub> <mover> <mi>z</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>1</mn> </msub> <mo>=</mo> <mo>-</mo> <msub> <mi>k</mi> <mn>1</mn> </msub> <msub> <mi>s</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>&amp;beta;</mi> <mn>1</mn> </msub> <mi>s</mi> <mi>i</mi> <mi>g</mi> <mi>n</mi> <mrow> <mo>(</mo> <msub> <mi>s</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>&amp;epsiv;</mi> <mn>1</mn> </msub> <msubsup> <mi>s</mi> <mn>1</mn> <mrow> <msub> <mi>p</mi> <mn>1</mn> </msub> <mo>/</mo> <msub> <mi>q</mi> <mn>1</mn> </msub> </mrow> </msubsup> <mo>-</mo> <mo>|</mo> <mrow> <msub> <mi>bx</mi> <mn>2</mn> </msub> </mrow> <mo>|</mo> <mi>s</mi> <mi>i</mi> <mi>g</mi> <mi>n</mi> <mrow> <mo>(</mo> <msub> <mi>s</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>x</mi> <mn>3</mn> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo> </mrow> </mrow>
式(11)中,k1、β1、ε1、p1和q1均为第二通道干扰观测器参数;p1<q1,均为正奇数,k1、β1、ε1均为正数,β1≥D1
由式(10)、(11)有:
<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mover> <mi>s</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>1</mn> </msub> <mo>=</mo> <msub> <mover> <mi>z</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>1</mn> </msub> <mo>-</mo> <msub> <mover> <mi>x</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>2</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <mo>-</mo> <msub> <mi>k</mi> <mn>1</mn> </msub> <msub> <mi>s</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>&amp;beta;</mi> <mn>1</mn> </msub> <mi>s</mi> <mi>i</mi> <mi>g</mi> <mi>n</mi> <mrow> <mo>(</mo> <msub> <mi>s</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>&amp;epsiv;</mi> <mn>1</mn> </msub> <msubsup> <mi>s</mi> <mn>1</mn> <mrow> <msub> <mi>p</mi> <mn>1</mn> </msub> <mo>/</mo> <msub> <mi>q</mi> <mn>1</mn> </msub> </mrow> </msubsup> <mo>-</mo> <mo>|</mo> <mrow> <msub> <mi>bx</mi> <mn>2</mn> </msub> </mrow> <mo>|</mo> <mi>s</mi> <mi>i</mi> <mi>g</mi> <mi>n</mi> <mrow> <mo>(</mo> <msub> <mi>s</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>bx</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>d</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>13</mn> <mo>)</mo> </mrow> </mrow>
故第二通道干扰d1(x,t)的观测值设计为其表达式如下:
<mrow> <msub> <mover> <mi>d</mi> <mo>^</mo> </mover> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mo>-</mo> <msub> <mi>k</mi> <mn>1</mn> </msub> <msub> <mi>s</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>&amp;beta;</mi> <mn>1</mn> </msub> <mi>s</mi> <mi>i</mi> <mi>g</mi> <mi>n</mi> <mrow> <mo>(</mo> <msub> <mi>s</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>&amp;epsiv;</mi> <mn>1</mn> </msub> <msubsup> <mi>s</mi> <mn>1</mn> <mrow> <msub> <mi>p</mi> <mn>1</mn> </msub> <mo>/</mo> <msub> <mi>q</mi> <mn>1</mn> </msub> </mrow> </msubsup> <mo>-</mo> <mo>|</mo> <mrow> <msub> <mi>bx</mi> <mn>2</mn> </msub> </mrow> <mo>|</mo> <mi>s</mi> <mi>i</mi> <mi>g</mi> <mi>n</mi> <mrow> <mo>(</mo> <msub> <mi>s</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>bx</mi> <mn>2</mn> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>14</mn> <mo>)</mo> </mrow> </mrow>
定义第二通道干扰观测器的Lyapunov方程:
<mrow> <msub> <mi>V</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msubsup> <mi>s</mi> <mn>1</mn> <mn>2</mn> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>15</mn> <mo>)</mo> </mrow> </mrow>
由于β1≥D1,故:
<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>s</mi> <mn>1</mn> </msub> <msub> <mover> <mi>s</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>1</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <msub> <mi>s</mi> <mn>1</mn> </msub> <mo>&amp;lsqb;</mo> <mo>-</mo> <msub> <mi>k</mi> <mn>1</mn> </msub> <msub> <mi>s</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>&amp;beta;</mi> <mn>1</mn> </msub> <mi>s</mi> <mi>i</mi> <mi>g</mi> <mi>n</mi> <mrow> <mo>(</mo> <msub> <mi>s</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>&amp;epsiv;</mi> <mn>1</mn> </msub> <msubsup> <mi>s</mi> <mn>1</mn> <mrow> <msub> <mi>p</mi> <mn>1</mn> </msub> <mo>/</mo> <msub> <mi>q</mi> <mn>1</mn> </msub> </mrow> </msubsup> <mo>-</mo> <mo>|</mo> <mrow> <msub> <mi>bx</mi> <mn>2</mn> </msub> </mrow> <mo>|</mo> <mi>s</mi> <mi>i</mi> <mi>g</mi> <mi>n</mi> <mrow> <mo>(</mo> <msub> <mi>s</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>bx</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>d</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <mo>-</mo> <msub> <mi>k</mi> <mn>1</mn> </msub> <msubsup> <mi>s</mi> <mn>1</mn> <mn>2</mn> </msubsup> <mo>-</mo> <msub> <mi>&amp;beta;</mi> <mn>1</mn> </msub> <msub> <mi>s</mi> <mn>1</mn> </msub> <mi>s</mi> <mi>i</mi> <mi>g</mi> <mi>n</mi> <mrow> <mo>(</mo> <msub> <mi>s</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>&amp;epsiv;</mi> <mn>1</mn> </msub> <msubsup> <mi>s</mi> <mn>1</mn> <mrow> <mo>(</mo> <msub> <mi>p</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>q</mi> <mn>1</mn> </msub> <mo>)</mo> <mo>/</mo> <msub> <mi>q</mi> <mn>1</mn> </msub> </mrow> </msubsup> <mo>-</mo> <mo>|</mo> <mrow> <msub> <mi>bx</mi> <mn>2</mn> </msub> </mrow> <mo>|</mo> <mo>|</mo> <msub> <mi>s</mi> <mn>1</mn> </msub> <mo>|</mo> <mo>+</mo> <msub> <mi>bx</mi> <mn>2</mn> </msub> <msub> <mi>s</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>d</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> <msub> <mi>s</mi> <mn>1</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>&amp;le;</mo> <mo>-</mo> <msub> <mi>k</mi> <mn>1</mn> </msub> <msubsup> <mi>s</mi> <mn>1</mn> <mn>2</mn> </msubsup> <mo>-</mo> <msub> <mi>&amp;beta;</mi> <mn>1</mn> </msub> <mo>|</mo> <msub> <mi>s</mi> <mn>1</mn> </msub> <mo>|</mo> <mo>-</mo> <msub> <mi>&amp;epsiv;</mi> <mn>1</mn> </msub> <msubsup> <mi>s</mi> <mn>1</mn> <mrow> <mo>(</mo> <msub> <mi>p</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>q</mi> <mn>1</mn> </msub> <mo>)</mo> <mo>/</mo> <msub> <mi>q</mi> <mn>1</mn> </msub> </mrow> </msubsup> <mo>+</mo> <msub> <mi>d</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> <msub> <mi>s</mi> <mn>1</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>&amp;le;</mo> <mo>-</mo> <msub> <mi>k</mi> <mn>1</mn> </msub> <msubsup> <mi>s</mi> <mn>1</mn> <mn>2</mn> </msubsup> <mo>-</mo> <msub> <mi>&amp;epsiv;</mi> <mn>1</mn> </msub> <msubsup> <mi>s</mi> <mn>1</mn> <mrow> <mo>(</mo> <msub> <mi>p</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>q</mi> <mn>1</mn> </msub> <mo>)</mo> <mo>/</mo> <msub> <mi>q</mi> <mn>1</mn> </msub> </mrow> </msubsup> <mo>=</mo> <mo>-</mo> <mn>2</mn> <msub> <mi>k</mi> <mn>1</mn> </msub> <msub> <mi>V</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>-</mo> <msup> <mn>2</mn> <mrow> <mo>(</mo> <msub> <mi>p</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>q</mi> <mn>1</mn> </msub> <mo>)</mo> <mo>/</mo> <mn>2</mn> <msub> <mi>q</mi> <mn>1</mn> </msub> </mrow> </msup> <msub> <mi>&amp;epsiv;</mi> <mn>1</mn> </msub> <msubsup> <mi>V</mi> <mn>1</mn> <mrow> <mo>(</mo> <msub> <mi>p</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>q</mi> <mn>1</mn> </msub> <mo>)</mo> <mo>/</mo> <mn>2</mn> <msub> <mi>q</mi> <mn>1</mn> </msub> </mrow> </msubsup> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>16</mn> <mo>)</mo> </mrow> </mrow>
又因,若存在一正定函数V0(t)满足下面不等式:
<mrow> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>&amp;alpha;V</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <msubsup> <mi>&amp;lambda;V</mi> <mn>0</mn> <mi>&amp;gamma;</mi> </msubsup> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>&amp;le;</mo> <mn>0</mn> <mo>,</mo> <mo>&amp;ForAll;</mo> <mi>t</mi> <mo>&gt;</mo> <msub> <mi>t</mi> <mn>0</mn> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>17</mn> <mo>)</mo> </mrow> </mrow>
则有,V0(t)在时间ts内收敛到平衡点,其中,
<mrow> <msub> <mi>t</mi> <mi>s</mi> </msub> <mo>&amp;le;</mo> <msub> <mi>t</mi> <mn>0</mn> </msub> <mo>+</mo> <mfrac> <mn>1</mn> <mrow> <mi>&amp;alpha;</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mi>&amp;gamma;</mi> <mo>)</mo> </mrow> </mrow> </mfrac> <mi>l</mi> <mi>n</mi> <mfrac> <mrow> <msubsup> <mi>&amp;alpha;V</mi> <mn>0</mn> <mrow> <mn>1</mn> <mo>-</mo> <mi>&amp;gamma;</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msub> <mi>t</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mi>&amp;lambda;</mi> </mrow> <mi>&amp;lambda;</mi> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>18</mn> <mo>)</mo> </mrow> </mrow>
式(18)中,α>0,λ>0,0<γ<1;
由式(17)、(18)式有,V1(t)将在有限时间内收敛到平衡点,即存在一个时刻t2,在经过有限时间t2后,s1为零,此时也将收敛到零,而d1(x,t)的观测误差为:
<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mover> <mi>d</mi> <mo>~</mo> </mover> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mover> <mi>d</mi> <mo>^</mo> </mover> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>d</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <mo>-</mo> <msub> <mi>k</mi> <mn>1</mn> </msub> <msub> <mi>s</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>&amp;beta;</mi> <mn>1</mn> </msub> <mi>s</mi> <mi>i</mi> <mi>g</mi> <mi>n</mi> <mrow> <mo>(</mo> <msub> <mi>s</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>&amp;epsiv;</mi> <mn>1</mn> </msub> <msubsup> <mi>s</mi> <mn>1</mn> <mrow> <msub> <mi>p</mi> <mn>1</mn> </msub> <mo>/</mo> <msub> <mi>q</mi> <mn>1</mn> </msub> </mrow> </msubsup> <mo>-</mo> <mo>|</mo> <mrow> <msub> <mi>bx</mi> <mn>2</mn> </msub> </mrow> <mo>|</mo> <mi>s</mi> <mi>i</mi> <mi>g</mi> <mi>n</mi> <mrow> <mo>(</mo> <msub> <mi>s</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>bx</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>d</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <mo>-</mo> <msub> <mi>k</mi> <mn>1</mn> </msub> <msub> <mi>s</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>&amp;beta;</mi> <mn>1</mn> </msub> <mi>s</mi> <mi>i</mi> <mi>g</mi> <mi>n</mi> <mrow> <mo>(</mo> <msub> <mi>s</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>&amp;epsiv;</mi> <mn>1</mn> </msub> <msubsup> <mi>s</mi> <mn>1</mn> <mrow> <msub> <mi>p</mi> <mn>1</mn> </msub> <mo>/</mo> <msub> <mi>q</mi> <mn>1</mn> </msub> </mrow> </msubsup> <mo>-</mo> <mo>|</mo> <mrow> <msub> <mi>bx</mi> <mn>2</mn> </msub> </mrow> <mo>|</mo> <mi>s</mi> <mi>i</mi> <mi>g</mi> <mi>n</mi> <mrow> <mo>(</mo> <msub> <mi>s</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>bx</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mover> <mi>x</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>2</mn> </msub> <mo>+</mo> <msub> <mi>x</mi> <mn>3</mn> </msub> <mo>-</mo> <msub> <mi>bx</mi> <mn>2</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <mo>-</mo> <msub> <mi>k</mi> <mn>1</mn> </msub> <msub> <mi>s</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>&amp;beta;</mi> <mn>1</mn> </msub> <mi>s</mi> <mi>i</mi> <mi>g</mi> <mi>n</mi> <mrow> <mo>(</mo> <msub> <mi>s</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>&amp;epsiv;</mi> <mn>1</mn> </msub> <msubsup> <mi>s</mi> <mn>1</mn> <mrow> <msub> <mi>p</mi> <mn>1</mn> </msub> <mo>/</mo> <msub> <mi>q</mi> <mn>1</mn> </msub> </mrow> </msubsup> <mo>-</mo> <mo>|</mo> <mrow> <msub> <mi>bx</mi> <mn>2</mn> </msub> </mrow> <mo>|</mo> <mi>s</mi> <mi>i</mi> <mi>g</mi> <mi>n</mi> <mrow> <mo>(</mo> <msub> <mi>s</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>x</mi> <mn>3</mn> </msub> <mo>-</mo> <msub> <mover> <mi>x</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>2</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <msub> <mover> <mi>z</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>1</mn> </msub> <mo>-</mo> <msub> <mover> <mi>x</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>2</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <msub> <mover> <mi>s</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>1</mn> </msub> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>19</mn> <mo>)</mo> </mrow> </mrow>
则干扰的观测误差也将在有限时间t2内为0;即经过时间t2
综上,得到第二通道干扰观测器表达式为:
<mrow> <msub> <mover> <mi>d</mi> <mo>^</mo> </mover> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mo>-</mo> <msub> <mi>k</mi> <mn>1</mn> </msub> <msub> <mi>s</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>&amp;beta;</mi> <mn>1</mn> </msub> <mi>s</mi> <mi>i</mi> <mi>g</mi> <mi>n</mi> <mrow> <mo>(</mo> <msub> <mi>s</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>&amp;epsiv;</mi> <mn>1</mn> </msub> <msubsup> <mi>s</mi> <mn>1</mn> <mrow> <msub> <mi>p</mi> <mn>1</mn> </msub> <mo>/</mo> <msub> <mi>q</mi> <mn>1</mn> </msub> </mrow> </msubsup> <mo>-</mo> <mo>|</mo> <mrow> <msub> <mi>bx</mi> <mn>2</mn> </msub> </mrow> <mo>|</mo> <mi>s</mi> <mi>i</mi> <mi>g</mi> <mi>n</mi> <mrow> <mo>(</mo> <msub> <mi>s</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>bx</mi> <mn>2</mn> </msub> </mrow>
步骤2-2、设计第三通道干扰观测器:
定义第三通道干扰观测器滑模面s2为:
s2=z2-x3 (20)
其中,z2为第三通道干扰观测器内部动态;
式(21)中,k2、β2、ε2、p2和q2均为第三通道干扰观测器参数;其中p2<q2,均为正奇数,k2、β2、ε2均为正数,β2≥D2
由式(20)、(21)有:
故第三通道干扰d2(x,t)的观测值设计为其表达式如下:
定义第三通道干扰观测器的Lyapunov方程:
<mrow> <msub> <mi>V</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msubsup> <mi>s</mi> <mn>2</mn> <mn>2</mn> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>24</mn> <mo>)</mo> </mrow> </mrow>
又因β2>D2,则有,
由式(17)、(18)式有,V2(t)将在有限时间内为零,即存在一个时刻t3,在经过t3后,s2为零,此时也将收敛到零,又因d2(x,t)估计误差
则干扰的估计误差也将在有限时间t3内为零,即经过时间t3后有:
综上,得到第三通道干扰观测器表达式为:
步骤3、设计基于干扰观测器的连续滑模控制器,具体如下:
定义液压伺服系统的位置跟踪误差e0(t)、速度的跟踪误差e1(t)、加速度的跟踪误差e2(t)、加速度导数的跟踪误差e3(t):
e0(t)=x1-xd(t) (27)
<mrow> <msub> <mi>e</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mover> <mi>e</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mover> <mi>x</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>1</mn> </msub> <mo>-</mo> <msub> <mover> <mi>x</mi> <mo>&amp;CenterDot;</mo> </mover> <mi>d</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>x</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mover> <mi>x</mi> <mo>&amp;CenterDot;</mo> </mover> <mi>d</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>28</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mi>e</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mover> <mi>e</mi> <mo>&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mover> <mi>x</mi> <mo>&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mn>1</mn> </msub> <mo>-</mo> <msub> <mover> <mi>x</mi> <mo>&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mi>d</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mover> <mi>x</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>2</mn> </msub> <mo>-</mo> <msub> <mover> <mi>x</mi> <mo>&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mi>d</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>29</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mi>e</mi> <mn>3</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mover> <mi>e</mi> <mo>&amp;CenterDot;&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mover> <mi>x</mi> <mo>&amp;CenterDot;&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mn>1</mn> </msub> <mo>-</mo> <msub> <mover> <mi>x</mi> <mo>&amp;CenterDot;&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mi>d</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mover> <mi>x</mi> <mo>&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mn>2</mn> </msub> <mo>-</mo> <msub> <mover> <mi>x</mi> <mo>&amp;CenterDot;&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mi>d</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>30</mn> <mo>)</mo> </mrow> </mrow>
其中,xd(t)为系统参考信号,xd(t)是三阶连续可导的,且系统参考位置信号xd(t)、系统参考速度信号系统参考加速度信号及系统参考加速度导数的信号都是有界的。
定义连续滑模控制器滑模面s为:
<mrow> <mi>s</mi> <mo>=</mo> <msub> <mi>x</mi> <mn>3</mn> </msub> <mo>-</mo> <msub> <mi>bx</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mover> <mi>x</mi> <mo>&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mi>d</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>c</mi> <mn>1</mn> </msub> <msub> <mi>e</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>c</mi> <mn>2</mn> </msub> <msub> <mi>e</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>c</mi> <mn>3</mn> </msub> <mo>&amp;Integral;</mo> <msub> <mi>e</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mi>d</mi> <mi>t</mi> <mo>-</mo> <mrow> <mo>(</mo> <mi>b</mi> <mo>-</mo> <msub> <mi>c</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <msub> <mi>s</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>s</mi> <mn>2</mn> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>31</mn> <mo>)</mo> </mrow> </mrow>
其中c1、c2、c3为滑模控制器参数,且均大于零,并且使得表达式是Hurwitz的,则有:
设计连续滑模控制器u的表达式为:
其中k3、k4、ζ为控制器参数,且k3>0、k4>0、0<ζ<1。
2.根据权利要求1所述的基于干扰补偿的电液位置伺服系统连续滑模控制方法,其特征在于,对步骤3设计的连续滑模控制器进行稳定性测试,具体如下:
将式(33)代入式(32)有:
<mrow> <mtable> <mtr> <mtd> <mrow> <mover> <mi>s</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>=</mo> <msub> <mi>d</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mover> <mi>d</mi> <mo>^</mo> </mover> <mn>2</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>-</mo> <mrow> <mo>(</mo> <mi>b</mi> <mo>-</mo> <msub> <mi>c</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>1</mn> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>-</mo> <msub> <mover> <mi>d</mi> <mo>^</mo> </mover> <mn>1</mn> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>)</mo> </mrow> <mo>-</mo> <mrow> <mo>(</mo> <mi>b</mi> <mo>-</mo> <msub> <mi>c</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <msub> <mover> <mi>s</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>1</mn> </msub> <mo>+</mo> <msub> <mover> <mi>s</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>k</mi> <mn>3</mn> </msub> <mi>s</mi> <mo>-</mo> <msub> <mi>k</mi> <mn>4</mn> </msub> <mi>s</mi> <mi>i</mi> <mi>g</mi> <mi>n</mi> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> <msup> <mrow> <mo>|</mo> <mi>s</mi> <mo>|</mo> </mrow> <mi>&amp;zeta;</mi> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <mo>-</mo> <msub> <mi>k</mi> <mn>3</mn> </msub> <mi>s</mi> <mo>-</mo> <msub> <mi>k</mi> <mn>4</mn> </msub> <mi>s</mi> <mi>i</mi> <mi>g</mi> <mi>n</mi> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> <msup> <mrow> <mo>|</mo> <mi>s</mi> <mo>|</mo> </mrow> <mi>&amp;zeta;</mi> </msup> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>34</mn> <mo>)</mo> </mrow> </mrow>
定义滑模控制器Lyapunov方程:
<mrow> <mi>V</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msup> <mi>s</mi> <mn>2</mn> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>35</mn> <mo>)</mo> </mrow> </mrow>
则有:
<mrow> <mtable> <mtr> <mtd> <mrow> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>s</mi> <mover> <mi>s</mi> <mo>&amp;CenterDot;</mo> </mover> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <mi>s</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>k</mi> <mn>4</mn> </msub> <mi>s</mi> <mo>-</mo> <msub> <mi>k</mi> <mn>5</mn> </msub> <mi>s</mi> <mi>i</mi> <mi>g</mi> <mi>n</mi> <mo>(</mo> <mi>s</mi> <mo>)</mo> <msup> <mrow> <mo>|</mo> <mi>s</mi> <mo>|</mo> </mrow> <mi>&amp;zeta;</mi> </msup> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <mo>-</mo> <mn>2</mn> <msub> <mi>k</mi> <mn>4</mn> </msub> <mi>V</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>-</mo> <msup> <mn>2</mn> <mrow> <mo>(</mo> <mi>&amp;zeta;</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> <mo>/</mo> <mn>2</mn> </mrow> </msup> <msub> <mi>k</mi> <mn>5</mn> </msub> <msup> <mi>V</mi> <mrow> <mo>(</mo> <mi>&amp;zeta;</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> <mo>/</mo> <mn>2</mn> </mrow> </msup> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>36</mn> <mo>)</mo> </mrow> </mrow>
故,由式(17)、(18)可知,V(t)将在有限时间内为零,即存在一个时刻t1,在经过有限时间t1后,s为零,即:
<mrow> <mi>s</mi> <mo>=</mo> <msub> <mi>x</mi> <mn>3</mn> </msub> <mo>-</mo> <msub> <mi>bx</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mover> <mi>x</mi> <mo>&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mi>d</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>c</mi> <mn>1</mn> </msub> <msub> <mi>e</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>c</mi> <mn>2</mn> </msub> <msub> <mi>e</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>c</mi> <mn>3</mn> </msub> <mo>&amp;Integral;</mo> <msub> <mi>e</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mi>d</mi> <mi>t</mi> <mo>-</mo> <mrow> <mo>(</mo> <mi>b</mi> <mo>-</mo> <msub> <mi>c</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <msub> <mi>s</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>s</mi> <mn>2</mn> </msub> <mo>=</mo> <mn>0</mn> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>37</mn> <mo>)</mo> </mrow> </mrow>
又因s1、s2也是有限时间内为0,t1为s为零的时刻,t2为s1为零的时刻,t3为s2为零的时刻,则存在t4=max{t1,t2,t3},经过t4时刻后有:
<mrow> <mi>s</mi> <mo>=</mo> <msub> <mi>x</mi> <mn>3</mn> </msub> <mo>-</mo> <msub> <mi>bx</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mover> <mi>x</mi> <mo>&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mi>d</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>c</mi> <mn>1</mn> </msub> <msub> <mi>e</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>c</mi> <mn>2</mn> </msub> <msub> <mi>e</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>c</mi> <mn>3</mn> </msub> <mo>&amp;Integral;</mo> <msub> <mi>e</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mi>d</mi> <mi>t</mi> <mo>=</mo> <mn>0</mn> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>38</mn> <mo>)</mo> </mrow> </mrow>
又因:
<mrow> <msub> <mi>x</mi> <mn>3</mn> </msub> <mo>-</mo> <msub> <mi>bx</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mover> <mi>x</mi> <mo>&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mi>d</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>e</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>d</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>39</mn> <mo>)</mo> </mrow> </mrow>
则:
e2(t)+c1e1(t)+c2e0(t)+c3∫e0(t)dt=d1(t) (40)
即:
<mrow> <msub> <mover> <mi>e</mi> <mo>&amp;CenterDot;&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>c</mi> <mn>1</mn> </msub> <msub> <mover> <mi>e</mi> <mo>&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>c</mi> <mn>2</mn> </msub> <msub> <mover> <mi>e</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>c</mi> <mn>3</mn> </msub> <msub> <mi>e</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mover> <mi>d</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>41</mn> <mo>)</mo> </mrow> </mrow>
时:
则有,当t→∞有:
<mrow> <msub> <mover> <mi>e</mi> <mo>&amp;CenterDot;&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>c</mi> <mn>1</mn> </msub> <msub> <mover> <mi>e</mi> <mo>&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>c</mi> <mn>2</mn> </msub> <msub> <mover> <mi>e</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>c</mi> <mn>3</mn> </msub> <msub> <mi>e</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mn>0</mn> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>42</mn> <mo>)</mo> </mrow> </mrow>
故e0(t)=x1-xd(t)在时间趋于无穷的条件下趋于零;
δ为一已知正数时:
则有,当t→∞有:
<mrow> <msub> <mover> <mi>e</mi> <mo>&amp;CenterDot;&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>c</mi> <mn>1</mn> </msub> <msub> <mover> <mi>e</mi> <mo>&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>c</mi> <mn>2</mn> </msub> <msub> <mover> <mi>e</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>c</mi> <mn>3</mn> </msub> <msub> <mi>e</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>&amp;le;</mo> <mi>&amp;delta;</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>43</mn> <mo>)</mo> </mrow> </mrow>
故e0(t)=x1-xd(t)在时间趋于无穷的条件下收敛于一个一致稳定的界内。
CN201510522508.2A 2015-08-24 2015-08-24 基于干扰补偿的电液位置伺服系统连续滑模控制方法 Active CN105068426B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510522508.2A CN105068426B (zh) 2015-08-24 2015-08-24 基于干扰补偿的电液位置伺服系统连续滑模控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510522508.2A CN105068426B (zh) 2015-08-24 2015-08-24 基于干扰补偿的电液位置伺服系统连续滑模控制方法

Publications (2)

Publication Number Publication Date
CN105068426A CN105068426A (zh) 2015-11-18
CN105068426B true CN105068426B (zh) 2018-03-09

Family

ID=54497817

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510522508.2A Active CN105068426B (zh) 2015-08-24 2015-08-24 基于干扰补偿的电液位置伺服系统连续滑模控制方法

Country Status (1)

Country Link
CN (1) CN105068426B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107165892B (zh) * 2017-06-14 2018-05-29 武汉科技大学 一种电液伺服系统的滑模控制方法
CN109541945B (zh) * 2019-01-10 2021-11-02 中国科学院光电技术研究所 一种基于复合型扰动观测器的扰动抑制方法
CN112068447A (zh) * 2020-09-30 2020-12-11 长沙航空职业技术学院 大型装备电液系统高精度鲁棒位置控制方法
CN112486021B (zh) * 2020-12-07 2021-10-08 燕山大学 一种针对非对称伺服液压位置跟踪系统的低复杂控制方法
CN115291507B (zh) * 2022-01-06 2023-05-26 兰州理工大学 一种矿山充填料浆浓度滑模控制方法及系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104238572A (zh) * 2014-07-23 2014-12-24 南京理工大学 基于扰动补偿的电机伺服系统无抖动滑模位置控制方法
CN104360635A (zh) * 2014-10-21 2015-02-18 南京理工大学 一种电机位置伺服系统的抗干扰控制方法
CN104614984A (zh) * 2014-11-20 2015-05-13 南京理工大学 一种电机位置伺服系统的高精度控制方法
CN104698844A (zh) * 2015-02-09 2015-06-10 南京理工大学 液压位置伺服系统的不确定性补偿的滑模控制方法
CN104808653A (zh) * 2015-04-24 2015-07-29 南京理工大学 基于滑模的电机伺服系统加性故障检测和容错控制方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104238572A (zh) * 2014-07-23 2014-12-24 南京理工大学 基于扰动补偿的电机伺服系统无抖动滑模位置控制方法
CN104360635A (zh) * 2014-10-21 2015-02-18 南京理工大学 一种电机位置伺服系统的抗干扰控制方法
CN104614984A (zh) * 2014-11-20 2015-05-13 南京理工大学 一种电机位置伺服系统的高精度控制方法
CN104698844A (zh) * 2015-02-09 2015-06-10 南京理工大学 液压位置伺服系统的不确定性补偿的滑模控制方法
CN104808653A (zh) * 2015-04-24 2015-07-29 南京理工大学 基于滑模的电机伺服系统加性故障检测和容错控制方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于干扰观测器的非线性不确定系统自适应滑模控制;于靖,等;《控制理论与应用》;20140831;第31卷(第8期);第993-995页 *
基于滑模状态观测器的电液位置伺服系统控制;乔继红;《计算机仿真》;20141130;第31卷(第11期);第367-371页 *

Also Published As

Publication number Publication date
CN105068426A (zh) 2015-11-18

Similar Documents

Publication Publication Date Title
CN105068426B (zh) 基于干扰补偿的电液位置伺服系统连续滑模控制方法
CN104698844B (zh) 液压位置伺服系统的不确定性补偿的滑模控制方法
CN104345639B (zh) 一种电液位置伺服系统鲁棒自适应控制方法
CN104635490B (zh) 一种单出杆液压缸位置伺服系统的输出反馈控制方法
CN104111607B (zh) 一种考虑输入时滞的电机位置伺服系统的控制方法
US11454177B2 (en) Method of aero-engine on-line optimization and multivariable control based on model prediction
CN104345638B (zh) 一种液压马达位置伺服系统的自抗扰自适应控制方法
CN105093936B (zh) 电液位置伺服系统的有限时间连续滑模控制方法
CN105700347B (zh) 一种含磁滞补偿的液压马达预设性能跟踪控制方法
CN106483844B (zh) 基于非线性鲁棒的电液伺服系统自适应位置控制器的实现方法
CN110308651B (zh) 基于扩张状态观测器的电液伺服系统全状态约束控制方法
CN104317198A (zh) 具有时变输出约束的电液伺服系统非线性鲁棒位置控制方法
CN104360635A (zh) 一种电机位置伺服系统的抗干扰控制方法
CN104898428A (zh) 一种基于干扰估计的电液伺服系统自适应鲁棒控制方法
CN107165892A (zh) 一种电液伺服系统的滑模控制方法
CN108107728A (zh) 一种基于干扰补偿的电液位置伺服系统控制方法
CN107092189A (zh) 基于模型预测控制的多变量输入eha系统的控制方法
CN107340714A (zh) 一种纳米伺服系统的自适应逆控制方法
CN114296346A (zh) 基于模糊自抗扰控制的电液位置伺服系统重构控制方法及装置
Zhao et al. A novel method of parameter identification of nonlinear electrohydraulic servo systems
Liang-kuan et al. Compound control strategy for MDF continuous hot pressing electrohydraulic servo system with uncertainties and input saturation
Man et al. Hydraulic impulse-testing system with pressure difference feedback from a second-stage valve
CN105739311A (zh) 基于预设回声状态网络的机电伺服系统受限控制方法
CN108345268B (zh) 考虑输入时滞约束的电液伺服系统位置跟踪控制方法
Zhu et al. Prescribed performance control for MDF continuous hot pressing hydraulic system

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
EE01 Entry into force of recordation of patent licensing contract
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20151118

Assignee: Liheng Technology (Jiangsu) Co.,Ltd.

Assignor: NANJING University OF SCIENCE AND TECHNOLOGY

Contract record no.: X2022980020673

Denomination of invention: A Continuous Sliding Mode Control Method for Electro hydraulic Position Servo System Based on Interference Compensation

Granted publication date: 20180309

License type: Exclusive License

Record date: 20221103