基于组合多模型的电液伺服系统重构控制设计

刘舒畅 杨忠 张钊 陈爽 张小恺

刘舒畅, 杨忠, 张钊, 等. 基于组合多模型的电液伺服系统重构控制设计 [J]. 应用科技, 2023, 50(1): 79-86. doi: 10.11991/yykj.202110002
引用本文: 刘舒畅, 杨忠, 张钊, 等. 基于组合多模型的电液伺服系统重构控制设计 [J]. 应用科技, 2023, 50(1): 79-86. doi: 10.11991/yykj.202110002
LIU Shuchang, YANG Zhong, ZHANG Zhao, et al. Reconfiguration control design for electro-hydraulic servo system based on combined multiple models [J]. Applied Science and Technology, 2023, 50(1): 79-86. doi: 10.11991/yykj.202110002
Citation: LIU Shuchang, YANG Zhong, ZHANG Zhao, et al. Reconfiguration control design for electro-hydraulic servo system based on combined multiple models [J]. Applied Science and Technology, 2023, 50(1): 79-86. doi: 10.11991/yykj.202110002

基于组合多模型的电液伺服系统重构控制设计

doi: 10.11991/yykj.202110002
基金项目: 航空科学基金项目(重点实验室类)(201928052006, 20162852031).
详细信息
    作者简介:

    刘舒畅,女,硕士研究生;

    杨忠, 男,教授,博士生导师.

    通讯作者:

    杨忠,E-mail: YangZhong@nuaa.edu.cn.

  • 中图分类号: TP277

Reconfiguration control design for electro-hydraulic servo system based on combined multiple models

  • 摘要: 针对电液伺服系统潜在的执行器、传感器等组件故障问题,提出了自适应控制与组合多模型相结合的重构控制策略。首先,根据电液伺服系统的物理组成和工作原理建立数学模型,并在其基础上对典型组件故障进行建模;其次,设计了基于模型跟随自适应控制的重构控制器,并给出其稳定性证明;然后,根据故障类型搭建其固定模型以及自适应模型,组成多模型集,模型集中的每个模型都有与之相对应的控制器,且每个控制器间相互独立;接着,提出合适的切换机制对故障系统进行模型匹配,使其自动切换到性能指标最小的控制器。仿真结果表明,该重构策略能够有效地补偿组件故障产生的不利影响,具有良好的鲁棒性和重构效果;同时,与基于单一模型的自适应控制进行对比,显示出其在故障出现瞬间具有更好的瞬态性能,抗扰性相对更强,电液伺服系统的安全性与可靠性得到提高。

     

    Abstract: In order to solve the problem of actuator and sensor faults in the electro-hydraulic servo system, we propose a reconfiguration control strategy based on combined multiple models and adaptive control. Firstly, the mathematical model of the electro-hydraulic servo system is established according to its physical composition and working principle, on which the fault modeling is carried out aiming at typical component faults; Then, a reconfiguration controller based on model following adaptive control is designed and its stability is proved; And further, according to the fault types, a combined multi-model set is introduced, which includes several fixed models and an adaptive model, with each model having its own independent controller. Then, an appropriate switching scheme is proposed to match the model of the faulty system and make it automatically switch to the controller with the minimum performance index. The simulation results show that the reconfiguration strategy can effectively compensate for the adverse effects of component faults, validating good robustness and reconfiguration effect. Compared with single model based adaptive control, better transient performance and better disturbance rejection of the proposed method are shown at the instantaneous moment of failure occurrence. The safety and reliability of the electro-hydraulic servo system are improved.

     

  • 作为典型的航空机电系统,电液伺服系统是集机、电、液于一身的复杂控制系统,具有控制精度高、响应速度快、信号处理灵活以及各种参数的反馈系统等优点[1]。电液伺服系统在航空航天工程中应用越来越广泛,特别是在大型飞机舵面控制以及飞行模拟器控制等[2]。然而电液伺服系统结构复杂,不可避免地存在组件故障隐患。一旦其某个部件发生故障,很容易引发连锁反应,甚至造成更严重的后果,系统的可靠性和安全性会因此大大降低。故障发生后如何保证系统稳定性以及令人满意的跟踪性能,使得电液伺服系统控制器设计成为一个有趣且具有挑战性的课题。

    在设计电液伺服系统时,有几个需要着重考虑的问题。电液伺服系统容易受到参数不确定性和外部干扰影响,例如外负载力、供油压力变化、液压油中混入空气以及内泄漏等。国内外针对电液伺服系统控制器的研究,主要是提出利用各种先进控制算法来克服其非线性和不确定性问题,例如模糊PI控制[3]、动态滑模控制[4]、自适应最优补偿控制[5]、基于反步法的非线性鲁棒控制[6]以及混合模型预测控制[7]等。尽管控制效果和抗干扰能力有所提高,但传感器故障和执行器退化等部件故障并未得到重视,电液伺服系统的可靠性和安全性没有被真正加强。因此,有必要为其设计一个能够自动容忍潜在故障的控制系统,即容错控制系统。文献[8]为了解决比例阀失效的问题,采用了由标准直驱比例阀组成的双阀系统。若1个阀门发生故障,将会被检测到,内部的控制结构将会转移到另一个阀门上。文献[9]针对电液伺服系统内泄露故障,在定量反馈理论的基础上设计了1个固定增益的线性定常控制律,在系统发生的故障的情况下,对故障进行补偿同时拥有良好的跟踪性能。文献[10]设计了基于免疫原理的电液伺服系统故障自愈策略,针对液压内泄漏故障和传感器故障实现了故障自修复。显然,上述研究所针对故障类型单一且方法不具有普适性。

    重构控制技术是确保一些系统安全运行的关键技术之一,在保证控制系统经济性的基础上,提高系统可靠性,多应用于航空航天领域[11-12]。自适应控制具有对系统参数和环境不确定性的适应性,使得其成为主流的重构控制方法[13-14]。然而,单一模型的控制器瞬态性能差,在面对比较严重的故障时,响应时间长,系统参数与无故障时相差甚远,从而导致较大的瞬态误差产生[15]。基于单一模型的自适应控制器需要更长的时间才能使故障系统接近新的运行状态,电液伺服系统的性能可能严重受损,甚至不可控。这使得相关的重构问题变得十分棘手。

    基于上述分析,本文采用基于多模型的自适应控制对电液伺服系统进行重构控制设计,建立多模型集来逼近系统故障模型。研究者们已经将多模型自适应控制广泛应用于重构问题上并取得不错的效果。针对飞控系统的执行器故障问题,文献[16]提出利用基于多模型自适应控制对其进行重构;文献[17]将多模型自适应控制应用于近空间飞行器以补偿执行器故障;文献[18]基于自适应控制和多模型组合,为具有电机动力部分失效故障的四旋翼直升机设计了最优重构控制方案。多模型自适应控制在其他领域的应用研究也取得了良好的成果[19-21]

    针对电液伺服系统潜在的执行器、传感器等组件故障问题,本文提出了自适应控制与组合多模型相结合的重构控制策略。首先建立电液伺服系统数学模型,然后设计了基于模型跟随自适应控制的重构控制器并给出其稳定性证明;紧接着设计1个由多个固定模型和1个自适应模型构成的组合多模型集,用以改善系统故障重构后的瞬态性能,并根据合适的切换机制完成相应控制器切换;最后通过仿真验证了所提出方法的优秀控制重构能力以及良好的鲁棒性和抗扰性。

    本文选择负载位移作为电液伺服系统输出。电液伺服系统的结构框图如图1所示,其中包括控制器、比例放大器、电液伺服阀、液压缸、位移传感器以及负载[22]图1中两电位器组成桥式电路,用以测量输入信号和输出信号之间的位置偏差。位置偏差信号经过控制器得到控制电压,经过放大器放大得到输出电流,用以控制电液伺服阀的阀芯运动。伺服阀的输出压力油推动液压缸活塞杆移动,从而带动工作台移动。

    图  1  电液伺服系统原理
    下载: 全尺寸图片
    1.2.1   液压缸数学模型

    液压缸的物理特性由以下方程描述。

    1)伺服阀的压力流量方程为

    $$ q_{L}=C_{d}A(x_v) \sqrt{\frac{p_s-p_L{\rm{sign}}(x_v)}{\rho}} $$ (1)

    式中: $ q_{L} $ 为负载流量, $ C_d $ 为流量系数,A为电液伺服阀的阀口面积函数, $ x_{v} $ 为伺服阀阀芯位移, $ p_{s} $ 为系统供油压力, $ p_{L} $ 为电液伺服阀的负载压力, $ \rho $ 为液压油密度。

    式(1)可以线性化为[23]

    $$ q_{L}=K_{q} x_{v}-K_{c} p_{L} $$

    式中: $ K_{q} $ 为流量增益, $K_{c}$ 为流量–压力系数。

    2)液压缸的连续性方程可表示为

    $$ q_{L}=A_{t} \frac{{\rm{d}} x}{{\rm{d}} t}+C_{tp} p_{L}+\frac{V_{t}}{4 \beta_{t}} \frac{{\rm{d}} p_{t}}{{\rm{d }}t} $$

    式中: $ x $ 为活塞杆的位移, $A_t$ 为液压缸的有效作用面积, $C_{tp}$ 为液压缸的等效泄露系数, $ \beta_{e} $ 为油液的等效体积弹性模量, $ V_{t} $ 为液压缸的总压缩容积。

    3)液压缸与负载的动力学方程可表示为

    $$ A_{t} p_{L}=m \frac{{\rm{d}}^{2} x}{{\rm{d}} t^{2}}+B_{v} \frac{{\rm{d}} x}{{\rm{d}} t}+k x+F_{L} $$

    式中: $ m $ 为负载作用在活塞杆上的等效质量, $ B_{v} $ 为负载和活塞的粘性阻尼系数, $ k $ 为负载作用在活塞上的弹簧刚度, $ F_{L} $ 为外负载力。

    实际应用中,由于系统的弹性负载和粘性阻尼系数都很小,可以忽略不计[24],即 $ k=0 $ $ B_{v}=0 $ 。令干扰作用 $ F_{L}=0 $ ,可得液压缸状态空间方程:

    $$ \left[ \begin{gathered} {\dddot x} \\ {\ddot x} \\ {\dot x} \\ \end{gathered} \right] = \left[ {\begin{array}{*{20}{c}} { - 2{\zeta _h}{\omega _h}}&{ - \omega _h^2}&0 \\ 1&0&0 \\ 0&1&0 \end{array}} \right]\left[ \begin{gathered} {\ddot x} \\ {\dot x} \\ x \\ \end{gathered} \right] + \left[ {\begin{array}{*{20}{c}} {\dfrac{{{K_q}\omega _h^2}}{{{A_h}}}} \\ 0 \\ 0 \end{array}} \right]{x_v} $$

    式中: $ \omega_{h} $ 为液压固有频率, $\omega_{h}=\sqrt{\dfrac{4 \beta_{e} A_{e}^{2}}{m V_{t}}}$ $ \zeta_{ {h }} $ 为液压阻尼比, $\zeta_{h}=\dfrac{K_{ce}}{A_{e}} \sqrt{\dfrac{\beta_{e} m}{V_{t}}}$

    此外,作为液压缸的一种典型故障,内泄漏是指液压缸的2个压力之间的液压油的相互渗透,主要是由活塞的磨损或液压缸内壁的磨损导致密封性减弱造成的。在仿真中,可以通过改变内泄漏系数 $C_{tp}$ 的值来进行故障的仿真分析[1]

    1.2.2   电液伺服阀的数学模型

    电液伺服阀可视为工作频率接近 $ \omega_{k} $ 的二阶震荡环节,其状态空间方程可表示为

    $$ \left[\begin{array}{l} \ddot{x}_{v} \\ \dot{x}_{v} \end{array}\right]=\left[\begin{array}{cc} -2 \zeta_{sv} \omega_{sv} & -\omega_{sv}^{2} \\ 1 & 0 \end{array}\right]\left[\begin{array}{l} \dot{x}_{v} \\ x_{v} \end{array}\right]+\left[\begin{array}{c} K_{sv} \omega_{sv}^{2} \\ 0 \end{array}\right] I_{r} $$

    式中: $I_r$ 为输入电流, $K_{ {sv }}$ 为伺服阀增益, $\omega_{sv}$ 为电液伺服阀的固有频率, $\zeta_{sv}$ 为电液伺服阀的阻尼比。

    电液伺服阀作为电液伺服系统的执行器,不可避免地受到潜在故障影响,液压油混入空气、内泄漏以及振动等问题严重影响电液伺服阀的工作效率。因此,本文引入了效率损失(loss of efficiency, LOE)来代表电液伺服系统典型的执行器故障。LOE故障的特点是执行器的增益从其额定值下降。在执行器LOE故障的情况下,伺服阀阀芯位移偏离了控制器所期望的指令输出,即

    $$ {x'_v} = {k_{{\text{LOE}}}}{x_v} $$

    式中: $ x_{v}^{\prime} $ 为伺服阀阀芯的实际位移; ${k_{{\text{LOE}}}}$ 为LOE故障增益, ${k_{{\text{LOE}}}} \in \left( {0,1} \right]$ $ n \% $ LOE等同于LOE故障增益 ${k_{{\text{LOE}}}} = 1 - \dfrac{n}{{100}}$ ${k_{{\text{LOE}}}} = 1$ 表示执行器无故障。

    因此,液压缸的数学模型可以改写为

    $$ \left[\begin{array}{c} \dddot{x} \\ \ddot{x} \\ \dot{x} \end{array}\right]=\left[\begin{array}{ccc} -2 \zeta_{h} \omega_{h} & -\omega_{h}^{2} & 0 \\ 1 & 0 & 0 \\ 0 & 1 & 0 \end{array}\right]\left[\begin{array}{c} \ddot{x} \\ \dot{x} \\ x \end{array}\right]+\left[\begin{array}{c} \dfrac{K_{q} \omega_{h}^{2}}{A_{h}} \\ 0 \\ 0 \end{array}\right] x_{v}^{\prime} $$

    需要注意的是,如果组件并不总是具有与无故障相同的特性,就有必要建立故障模型。它不仅清楚地说明了如何将组件故障注入系统,而且为接下来设计多模型集提供了一些准确的故障模型。

    1.2.3   位移传感器数学模型

    位移传感器视为比例环节:

    $$ x_{s}=K_{f} x $$

    式中: $ x_{s} $ 为传感器测量位移, $ K_{f} $ 为传感器反馈增益。

    与执行器故障类似,传感器增益漂移、失效等典型故障也会严重影响系统性能,甚至导致系统不可控。传感器故障可以表示为

    $$ {x'_s} = {k_{{\text{DFT}}}}{x_s} $$

    式中: $ x_{s}^{\prime} $ 为传感器的实际输出; ${k_{{\text{DFT}}}}$ 为传感器漂移系数, ${k_{{\text{DFT}}}} \in \left[ {0,2} \right]$ ${k_{{\text{DFT}}}} = 1$ 表示传感器无故障, ${k_{{\text{DFT}}}} = 0$ 表示传感器失效。

    1.2.4   控制器数学模型

    本文采用PI控制器来保证电液伺服系统的基本跟踪性能,其传递函数如下:

    $$ {u_{s e}} = {K_p}{x_e} + {K_i}\int {{x_e}} {\rm{d}}t $$

    式中: ${u_{s e}}$ 为控制器输出电压, $ K $ $ K_{i} $ 分别为比例增益和积分增益, $ x_{e}=x_{d}-x_{s} $ 是期望位移与测量位移之差。

    1.2.5   伺服放大器数学模型

    伺服放大器方程为

    $$ I_{r}=K_{a} u_{s e} $$

    式中 $ K_{a} $ 是伺服放大器增益。

    无故障的电液伺服系统状态空间表示形式为

    $$ \left\{ \begin{gathered} {\boldsymbol{\dot x}}\left({{t}} \right) = {\boldsymbol{Ax}}\left({{t}} \right){\boldsymbol{ + Bu}}\left({{t}} \right) \\ {\boldsymbol{y}}\left({{t}} \right) = {\boldsymbol{Cx}}\left({{t}} \right) \\ \end{gathered} \right. $$

    式中: ${\boldsymbol{x}}\left({{t}} \right) \in {{\boldsymbol{R}}^n}$ 为状态向量, ${\boldsymbol{u}}\left({{t}} \right) \in {{\boldsymbol{R}}^m}$ 为控制输入向量, ${\boldsymbol{y}}\left({{t}} \right) \in {{\boldsymbol{R}}^q}$ 为输出向量, ${\boldsymbol{A}} \in {{\boldsymbol{R}}^{n \times n}}$ ${\boldsymbol{B}} \in {{\boldsymbol{R}}^{n \times m}}$ ${\boldsymbol{C}} \in {{\boldsymbol{R}}^{q \times n}}$

    本文主要针对位移传感器、执行器(电液伺服阀)以及液压缸等组件故障。故障系统需要进行统一表示来设计重构控制器,可以将其描述为

    $$ \left\{ \begin{gathered} {{{\boldsymbol{\dot x}}}_{{p}}}\left({{t}} \right) = {{\boldsymbol{A}}_{{p}}}{{\boldsymbol{x}}_{{p}}}\left({{t}} \right) + {{\boldsymbol{B}}_{{p}}}{{\boldsymbol{u}}_{{p}}}\left({{t}} \right) \\ {{\boldsymbol{y}}_{{p}}}\left({{t}} \right) = {{\boldsymbol{C}}_{{p}}}{{\boldsymbol{x}}_{{p}}}\left({{t}} \right) \\ \end{gathered} \right. $$ (2)
    $$ \left\{ \begin{gathered} {{\boldsymbol{A}}_{{p}}} = {\boldsymbol{A}} + \sigma {\boldsymbol{A}} \\ {{\boldsymbol{B}}_{{p}}} = {\boldsymbol{B}} + \sigma {\boldsymbol{B}} \\ {{\boldsymbol{C}}_{{p}}} = {\boldsymbol{C}} + \sigma {\boldsymbol{C}} \\ \end{gathered} \right. $$

    式中: ${{\boldsymbol{x}}_{{p}}}\left({{t}} \right) \in {{\boldsymbol{R}}^{{n_p}}}$ ${{\boldsymbol{u}}_{{p}}}\left({{t}} \right) \in {{\boldsymbol{R}}^{{m_p}}}$ ${{\boldsymbol{y}}_{{p}}}\left({{t}} \right) \in {{\boldsymbol{R}}^{{q_p}}}$ ${{\boldsymbol{A}}_{{p}}} \in {{\boldsymbol{R}}^{{n_p} \times {n_p}}}$ ${{\boldsymbol{B}}_{{p}}} \in {{\boldsymbol{R}}^{{n_p} \times {m_p}}}$ ${{\boldsymbol{C}}_{{p}}} \in {{\boldsymbol{R}}^{{q_p} \times {n_p}}}$ $\left\{ {\sigma {\boldsymbol{A}},\sigma {\boldsymbol{B}},\sigma {\boldsymbol{C}}} \right\}$ 是由故障引起的具有适当维数的有界参数摄动矩阵。

    模型跟随自适应控制针对可以用参数未知确定性线性模型描述的被控对象,所采用的控制器是基于参考模型输出与被控对象输出之间的广义误差和稳定性理论进行设计的[25]

    结合模型跟随自适应控制思想,为故障系统设计重构控制器。取参考模型状态方程为

    $$ \left\{ \begin{gathered} {{{\boldsymbol{\dot x}}}_{{m}}}\left({{t}} \right) = {{\boldsymbol{A}}_{{m}}}{{\boldsymbol{x}}_{{m}}}\left({{t}} \right) + {{\boldsymbol{B}}_{{m}}}{\boldsymbol{r}}\left({{t}} \right) \\ {{\boldsymbol{y}}_{{m}}}\left({{t}} \right) = {{\boldsymbol{C}}_{{m}}}{{\boldsymbol{x}}_{{m}}}\left({{t}} \right) \\ \end{gathered} \right. $$

    式中: ${{\boldsymbol{x}}_{{m}}}\left({{t}} \right) \in {{\boldsymbol{R}}^{{n_m}}}$ ${\boldsymbol{r}}\left({{t}} \right) \in {{\boldsymbol{R}}^{{m_m}}}$ ${{\boldsymbol{y}}_{{m}}}\left({{t}} \right) \in {{\boldsymbol{R}}^{{q_m}}}$ ${{\boldsymbol{A}}_{{p}}} \in {{\boldsymbol{R}}^{{n_m} \times {n_m}}}$ ${{\boldsymbol{B}}_{{m}}} \in {{\boldsymbol{R}}^{{n_m} \times {m_{{m}}}}}$ ${{\boldsymbol{C}}_{{m}}} \in {{\boldsymbol{R}}^{{q_m} \times {n_{{m}}}}}$

    最终的控制目标是保持闭环系统里的所有信号均有界的同时,使电液伺服系统的输出 ${{\boldsymbol{y}}_{{p}}}\left({{t}} \right)$ 渐进跟踪参考模型输出 ${{\boldsymbol{y}}_{{m}}}\left({{t}} \right)$

    因此,设计的重构控制器如下:

    $$ {{\boldsymbol{u}}_{{{ad}}}}\left({{t}} \right) = {\boldsymbol{Kr}} + {\boldsymbol{F}}{{\boldsymbol{x}}_{{p}}} $$ (3)

    式中 ${\boldsymbol{K}} \in {{\boldsymbol{R}}^{{m_p} \times {m_{{p}}}}}$ ${\boldsymbol{F}} \in {{\boldsymbol{R}}^{{m_p} \times {m_{{p}}}}}$ 为自适应控制增益矩阵。

    将式(3)代入式(2)中,得:

    $$ {{\boldsymbol{\dot x}}_{{p}}} = {{\boldsymbol{A}}_{{p}}}{{\boldsymbol{x}}_{{p}}} + {{\boldsymbol{B}}_{{p}}}\left( {{\boldsymbol{Kr}} + {\boldsymbol{F}}{{\boldsymbol{x}}_{{p}}}} \right) $$

    定义状态误差为 ${\boldsymbol{e}} = {{\boldsymbol{x}}_{{p}}} - {{\boldsymbol{x}}_{{m}}}$ ,则:

    $$ \begin{split} &\qquad\qquad {\boldsymbol{\dot e}} = {{{\boldsymbol{\dot x}}}_{{p}}} - {{{\boldsymbol{\dot x}}}_{{m}}} = \\&\;\;\; {{\boldsymbol{A}}_{{p}}}{{\boldsymbol{x}}_{{p}}} + {{\boldsymbol{B}}_{{p}}}\left( {{\boldsymbol{Kr}} + {\boldsymbol{F}}{{\boldsymbol{x}}_{{p}}}} \right) - {{\boldsymbol{A}}_{{m}}}{{\boldsymbol{x}}_{{m}}} - {{\boldsymbol{B}}_{{m}}}{\boldsymbol{r}} = \\&\;\;\; \left( {{{\boldsymbol{A}}_{{p}}} + {{\boldsymbol{B}}_{{p}}}{\boldsymbol{F}}} \right){{\boldsymbol{x}}_{{p}}} + \left( {{{\boldsymbol{B}}_{{p}}}{\boldsymbol{K}} - {{\boldsymbol{B}}_{{m}}}} \right){\boldsymbol{r}} - {{\boldsymbol{A}}_{{m}}}{{\boldsymbol{x}}_{{m}}} - \\& \left( {{{\boldsymbol{A}}_{{p}}} + {{\boldsymbol{B}}_{{p}}}{\boldsymbol{F}} - {{\boldsymbol{A}}_{{m}}}} \right){{\boldsymbol{x}}_{{p}}} + \left( {{{\boldsymbol{B}}_{{p}}}{\boldsymbol{K}} - {{\boldsymbol{B}}_{{m}}}} \right){\boldsymbol{r}} + {{\boldsymbol{A}}_{{m}}}{\boldsymbol{e}} \end{split} $$ (4)

    为了保证 ${\boldsymbol{e}} = {{\boldsymbol{x}}_{{p}}} - {{\boldsymbol{x}}_{{m}}}$ 渐进收敛,需满足下列匹配条件:

    $$ \left\{ \begin{gathered} {{\boldsymbol{A}}_{{p}}} + {{\boldsymbol{B}}_{{p}}}{\boldsymbol{F}} = {{\boldsymbol{A}}_{{m}}} \\ {{\boldsymbol{B}}_{{p}}}{\boldsymbol{K}} = {{\boldsymbol{B}}_{{m}}} \\ \end{gathered} \right. $$ (5)

    将式(5)代入式(4)得到:

    $$ {\boldsymbol{\dot e}} = {{\boldsymbol{A}}_{{m}}}{\boldsymbol{e}} + {{\boldsymbol{B}}_{{m}}}{{\boldsymbol{K}}^{{{ - 1}}}}{{\tilde {\boldsymbol{F}}}}{{\boldsymbol{x}}_{{p}}} + {{\boldsymbol{B}}_{{m}}}{{\boldsymbol{K}}^{{{ - 1}}}}{{\tilde {\boldsymbol{Kr}}}} $$ (6)

    式中: ${{\tilde {\boldsymbol{F}}}} = {{\boldsymbol{F}}_{{0}}} - {\boldsymbol{F}}$ ${{\tilde {\boldsymbol{K}}}} = {{\boldsymbol{K}}_{{0}}} - {\boldsymbol{K}}$

    为保证 $\mathop {\lim }\limits_{x \to \infty } {\boldsymbol{e}}=0$ ,设计如下自适应律:

    $$ \left\{ \begin{split} {{\dot {\boldsymbol{F}}}} = - {{\boldsymbol{\gamma }}_{{F}}}{\boldsymbol{B}}_m^{\rm{T}}{\boldsymbol{Pex}}_p^{\rm{T}} \\ {{\dot {\boldsymbol{K}}}} = - {{\boldsymbol{\gamma }}_{{K}}}{\boldsymbol{B}}_m^{\rm{T}}{\boldsymbol{Pe}}{{\boldsymbol{r}}^{\rm{T}}} \\ \end{split} \right. $$ (7)

    式中 ${{\boldsymbol{\gamma }}_{{F}}} \in {{\boldsymbol{R}}^{{m_p} \times {m_p}}}$ ${{\boldsymbol{\gamma }}_{{K}}} \in {{\boldsymbol{R}}^{{m_p} \times {m_p}}}$ 均为对称的正定矩阵。

    定理  将所设计的自适应控制器(3)和自适应律(7)应用于故障后的电液伺服系统(2),能够保证闭环系统的所有信号均有界,同时满足: $\mathop {\lim }\limits_{t \to \infty } {\boldsymbol{e}}\left({{t}} \right) = 0$ $\mathop {\lim }\limits_{t \to \infty } {{\boldsymbol{e}}_{{y}}}\left({{t}} \right) = \mathop {\lim }\limits_{t \to \infty } {{\boldsymbol{C}}_{{p}}}{{\boldsymbol{e}}_{{x}}}\left({{t}} \right) = 0$

    证明  构造如下Lyapunov函数:

    $$ {\boldsymbol{V}} = {{\boldsymbol{e}}^{\text{T}}}{\boldsymbol{Pe}} + {\text{tr}}\left( {{{{{\tilde {\boldsymbol{F}}}}}^{\text{T}}}{\boldsymbol{P}}_{{F}}^{{{ - }}{{1}}}{{\tilde {\boldsymbol{F}}}}} \right) + {\text{tr}}\left( {{{{{\tilde {\boldsymbol{K}}}}}^{\text{T}}}{\boldsymbol{P}}_{{K}}^{{{ - 1}}}{{\tilde {\boldsymbol{K}}}}} \right) $$ (8)

    式中: ${\boldsymbol{P = }}{{\boldsymbol{P}}^{\text{T}}} \in {{\boldsymbol{R}}^{{n_p} \times {n_p}}} > 0$ 满足对于任意的常数矩阵 ${\boldsymbol{Q = }}{{\boldsymbol{Q}}^{\text{T}}} \in {{\boldsymbol{R}}^{{n_p} \times {n_p}}} > 0$ ,都有:

    $$ {\boldsymbol{A}}_{{m}}^{\text{T}}{\boldsymbol{P}} + {\boldsymbol{P}}{{\boldsymbol{A}}_{{m}}} = - {\boldsymbol{Q}} $$ (9)

    结合式(8)和式(9),对式(8)求导可得

    $$\begin{split}&\;\; \dot {\boldsymbol{V}}=\dot {\boldsymbol{e}}^{\rm{T}}{\boldsymbol{P}}{\boldsymbol{e}}+{\boldsymbol{e}}^{\rm{T}}{\boldsymbol{P}}{\dot {\boldsymbol{e}}}+{\rm{tr}}({\dot {\tilde{\boldsymbol{F}^{\rm{T}}}}}{\boldsymbol{P}}_{F}^{-1} {\tilde{{\boldsymbol{F}}}}+{ {\tilde{\boldsymbol{F}^{\rm{T}}}}}{\boldsymbol{P}}_{F}^{-1} \dot{\tilde{{\boldsymbol{F}}}})+ \\ &{\boldsymbol{e}}^{\rm{T}}{\rm{tr}}({\dot {\tilde{\boldsymbol{K}^{\rm{T}}}}}{\boldsymbol{P}}_{K}^{-1} {\tilde{{\boldsymbol{K}}}}+{ {\tilde{\boldsymbol{K}^{\rm{T}}}}}{\boldsymbol{P}}_{K}^{-1} \dot{\tilde{{\boldsymbol{K}}}})={\boldsymbol{e}}^{\rm{T}}({\boldsymbol{A}}_m^{\rm{T}}{\boldsymbol{P}}+{\boldsymbol{PA}}_m){\boldsymbol{e}}+\\ &\qquad{\rm{2tr}}({\dot {\tilde{\boldsymbol{F}^{\rm{T}}}}}{\boldsymbol{P}}_{F}^{-1} {\tilde{{\boldsymbol{F}}}}+{\boldsymbol{x}}_p{\boldsymbol{e}}^{\rm{T}}{\boldsymbol{PB}}_m{\boldsymbol{K}}^{-1} {\tilde{{\boldsymbol{F}}}})+\\ &\qquad{\rm{2tr}}({\dot {\tilde{\boldsymbol{K}^{\rm{T}}}}}{\boldsymbol{P}}_{K}^{-1} {\tilde{{\boldsymbol{K}}}}+{\boldsymbol{re}}^{\rm{T}}{\boldsymbol{PB}}_m{\boldsymbol{K}}^{-1} {\tilde{{\boldsymbol{K}}}})=\\ &\;\;-{\boldsymbol{e}}^{\rm{T}}{\boldsymbol{Qe}}+{\rm{2tr}}[({\dot {\tilde{\boldsymbol{F}^{\rm{T}}}}}{\boldsymbol{P}}_{F}^{-1} {\tilde{{\boldsymbol{F}}}}+{\boldsymbol{x}}_p{\boldsymbol{e}}^{\rm{T}}{\boldsymbol{PB}}_m{\boldsymbol{K}}^{-1} {\tilde{{\boldsymbol{F}}}})+\\ &\;\;{\rm{2tr}}({\dot {\tilde{\boldsymbol{K}^{\rm{T}}}}}{\boldsymbol{P}}_{K}^{-1} {\tilde{{\boldsymbol{K}}}}+{\boldsymbol{re}}^{\rm{T}}{\boldsymbol{PB}}_m{\boldsymbol{K}}^{-1} {\tilde{{\boldsymbol{K}}}})]=-{\boldsymbol{e}}^{\rm{T}}{\boldsymbol{Qe}}\leqslant0 \end{split}$$ (10)

    由此可得: ${\text{0}} \leqslant {\boldsymbol{V}}\left({{t}} \right) \leqslant {\boldsymbol{V}}\left( {\boldsymbol{0}} \right)$ ,则 ${\boldsymbol{V}}\left({{t}} \right) \in {{\boldsymbol{L}}_\infty }$

    对式(10)进行积分,得到:

    $$ \int\limits_{\text{0}}^\infty { - {{\boldsymbol{e}}^{\text{T}}}{\boldsymbol{Qe}}dt = - \int\limits_{\text{0}}^\infty {{\boldsymbol{\dot V}}dt = {\boldsymbol{V}}\left( 0 \right) - {\boldsymbol{V}}\left( \infty \right)} } < \infty $$

    因此 ${\boldsymbol{e}} \in {{\boldsymbol{L}}_{\text{2}}},{{\boldsymbol{L}}_\infty }$ ${\boldsymbol{\dot e}} \in {{\boldsymbol{L}}_\infty }$ ,证毕。

    组件故障会引起较大范围的参数跳变,基于单模型的自适应控制器可能需要较长时间来适应新的工作状况。在这段时间间隔内,被控系统可能会产生在实际情况中不能接受的大的瞬态误差。因此,使用单模型自适应控制器可能无法达到预期的控制目标。在这种情况下,设计了如图2所示的基于组合多模型的自适应控制器,来更好地补偿故障产生的不利影响。该方法主要是通过建立多模型集来逼近被控系统当前状态,模型集中的每个模型都有与之相对应的控制器,且每个控制器间相互独立,依据一定的切换性能指标完成对系统的稳定控制。多模型的建立、控制器的设计和切换策略的选择是设计的核心内容。

    图  2  基于组合多模型的自适应重构系统原理
    下载: 全尺寸图片

    定义参数向量 ${{\boldsymbol{\theta }}_{{i}}} = \left[ {{{\boldsymbol{A}}_{{{pi}}}},{{\boldsymbol{B}}_{{{pi}}}}} \right]$ ${\boldsymbol{W}}\left({{t}} \right){\boldsymbol{ = }}\Big[ {\boldsymbol{x}}_{{p}}^{\text{T}}\left({{t}} \right), {\boldsymbol{u}}_{{p}}^{\text{T}} \left({{t}} \right) \Big]^{\text{T}}$ ,则第 $ f $ 个固定模型可表示为

    $$ {{\boldsymbol{\dot x}}_{{{pi}}}}\left({{t}} \right) = {{\boldsymbol{\theta }}_{{i}}}{\boldsymbol{W}}\left({{t}} \right) $$

    式中 $i = 1,2, \cdots ,n$

    要使被控对象的输出跟踪参考模型的输出,即 ${{\boldsymbol{y}}_{{p}}}\left({{t}} \right) = {{\boldsymbol{y}}_{{m}}}\left({{t}} \right)$ ,则有:

    $$ \begin{split} {{\boldsymbol{u}}_{{{pi}}}}\left({{t}} \right) =& {\left( {{{\boldsymbol{C}}_{{{pi}}}}{{\boldsymbol{B}}_{{{pi}}}}} \right)^{{{ - 1}}}}\left[ {{{\boldsymbol{C}}_{{m}}}{{\boldsymbol{A}}_{{m}}}{{\boldsymbol{x}}_{{m}}}\left({{t}} \right) + {{\boldsymbol{C}}_{{m}}}{{\boldsymbol{B}}_{{m}}}{\boldsymbol{r}}\left({{t}} \right)} -\right. \\ &\left.\qquad\qquad { {{\boldsymbol{C}}_{{{pi}}}}{{\boldsymbol{A}}_{{{pi}}}}{{\boldsymbol{x}}_{{p}}}\left({{t}} \right)} \right] \\ \end{split} $$ (11)

    由式(11)即可求得各种固定模型所对应的控制器。

    针对故障系统,可以构建如下的自适应辨识模型集:

    $$ \begin{split}\begin{cases} \dot{\hat {\boldsymbol{x}}}_{{pk}}={\boldsymbol{A}}_p\hat {\boldsymbol{x}}_{pk}+{\boldsymbol{B}}_p\hat {\boldsymbol{u}}_k+{\boldsymbol{L}}({\boldsymbol{y}}_p-\hat {\boldsymbol{y}}_{pk})\\ \hat {\boldsymbol{y}}_{pk}={\boldsymbol{C}}_p\hat {\boldsymbol{x}}_{pk}(k=1,2,\cdots,n) \end{cases} \end{split} $$ (12)

    定义状态误差 ${{\boldsymbol{\hat e}}_{{k}}} = {{\boldsymbol{\hat x}}_{{{pk}}}} - {{\boldsymbol{x}}_{{p}}}$ ,求导可得

    $$ \dot{\hat {\boldsymbol{e}}}_{k}=({\boldsymbol{A}}_p-{\boldsymbol{LC}}_{p})\hat {\boldsymbol{e}}_k+{\boldsymbol{B}}_p{\boldsymbol{\psi}}_{u_k} $$

    式中 ${{\boldsymbol{\psi }}_{{{{u}}_{{k}}}}} = {{\boldsymbol{\hat u}}_{{k}}} - {{\boldsymbol{u}}_{{p}}}$ ,且有:

    $$ \dot{\hat {\boldsymbol{u}}}_{k}={\rm{Proj}}({\hat {\boldsymbol{u}}}_{k},-\gamma_{u_k}{\hat {\boldsymbol{e}}}_k^{\rm{T}}{\boldsymbol{P}}_0{\boldsymbol{B}}_P)$$ (13)

    式中: $-\gamma_{u_k}>0$ 是权值系数; $ \text { Proj } $ 为投影算子,保证参数有界性; ${{\boldsymbol{P}}_{{0}}}$ 是如下Lyapunov方程的1个对称正定解。

    $$ \begin{split}\begin{cases} {\left( {{{\boldsymbol{A}}_{{p}}} - {\boldsymbol{L}}{{\boldsymbol{C}}_{{p}}}} \right)^{\text{T}}}{{\boldsymbol{P}}_{{0}}} + {{\boldsymbol{P}}_{{0}}}\left( {{{\boldsymbol{A}}_{{p}}} - {\boldsymbol{L}}{{\boldsymbol{C}}_{{p}}}} \right) = - {{\boldsymbol{Q}}_{{0}}} \\ {{\boldsymbol{P}}_{{0}}}{{\boldsymbol{B}}_{{p}}} = {\boldsymbol{C}}_{{p}}^{\text{T}}\end{cases} \end{split} $$

    ${{\boldsymbol{Q}}_{{0}}} = {\boldsymbol{Q}}_{{0}}^{\text{T}} > 0$

    故障引起系统参数产生大跳变,虽然固定模型具有不错的快速性,但需要固定模型的数量足够多,才能覆盖故障系统的全参数域;自适应模型可以在被控对象发生故障后动态调整自身参数,用以提高控制精度,但由于计算量较大导致时间滞后的缺点也显而易见。由此可见,多模型集中元素种类和数量的选择直接影响着重构控制的性能。

    为在多模型中切换以获取“最佳”模型,选择如下的切换性能指标[26]

    $$ J_{k}(t)=c_{1}\left\|\varepsilon_{k}\right\|^{2}+c_{2} \int_{0}^{t} \exp (-\lambda t)\left\|\varepsilon_{k}\right\|^{2} {\rm{d}} t $$ (14)

    式中: $\varepsilon_{k}=\hat{y}_{pk}-y_{p}$ 为残差, $ c_{1}, c_{2}>0 $ 分别为瞬时误差和累积误差的加权系数, $ \lambda \geq 0 $ 为遗忘因子。在每个时刻对 $ J_{k}(t) $ 进行比较以搜索最小值,以搜索与当前系统最匹配的模型,这时 $\mathop {\lim }\limits_{t \to \infty } \varepsilon_{k}(t)=0$ $\mathop {\lim }\limits_{t \to \infty } J_{k}(t)=0$

    定理  切换逻辑式(14)保证系统式(12)和调节律式(13)是稳定的,并且有 $\mathop {\lim }\limits_{t \to \infty } J_{k}(t)=0$

    证明  为了证明,给出下面的Lyapunov函数[27]:

    $$ {{\boldsymbol{V}}_{{k}}}\left( {{{{\boldsymbol{\hat e}}}_{{k}}}{\boldsymbol{,}}{{\boldsymbol{\psi }}_{{{{u}}_{{k}}}}}} \right) = - \frac{1}{2}\left( {{\boldsymbol{\hat e}}_{{k}}^{{{\rm{T}}}}{{\boldsymbol{P}}_{{0}}}{{\boldsymbol{e}}_{{k}}} + \frac{{{\boldsymbol{\psi }}_{{{{u}}_{{k}}}}^{\boldsymbol{2}}}}{{{\gamma _k}}}} \right) $$ (15)

    对式(15)进行微分可得:

    $$ {{\boldsymbol{\dot V}}_{{k}}}\left( {{{{\boldsymbol{\hat e}}}_{{k}}}{\boldsymbol{,}}{{\boldsymbol{\psi }}_{{{{u}}_{{k}}}}}} \right) \leqslant - \frac{1}{2}\left( {{{\boldsymbol{\lambda }}_{{m}}}\left\| {{\boldsymbol{\hat e}}_{{k}}^{\boldsymbol{2}}} \right\|} \right) \leqslant 0 $$

    考虑到自适应映射算法的特性 ${{\boldsymbol{\psi }}_{{{{u}}_{{k}}}}}{{\boldsymbol{\dot \psi }}_{{{{u}}_{{k}}}}} \leqslant - {\gamma _k}{{\boldsymbol{\psi }}_{{{{u}}_{{k}}}}} \left( {{{{\boldsymbol{\hat e}}}_{{k}}}{{\boldsymbol{P}}_{{0}}}{{\boldsymbol{B}}_{{p}}} + {{{\boldsymbol{\dot u}}}_{{k}}}} \right)$ ,则上述结果意味着 ${{\boldsymbol{\hat e}}_{{k}}} \in {{\boldsymbol{L}}^\infty } \cap {{\boldsymbol{L}}^2}$ , $ e_kL \in {{\boldsymbol{L}}^\infty } $ $\mathop {\lim }\limits_{t \to \infty } {{\boldsymbol{J}}_{{k}}}\left({{t}} \right) = 0$ 。当系统发生故障后,切换函数将切换到相应模型,因此 ${{\boldsymbol{\hat e}}_{{k}}}$ 是有界的(在映射算法下 ${{\boldsymbol{\psi }}_{{{{u}}_{{k}}}}}$ 是有界的),进而 ${{\boldsymbol{\hat x}}_{{{pk}}}}$ ${{\boldsymbol{x}}_{{{pk}}}}$ 是有界的。所以 ${{\boldsymbol{u}}_{{k}}}$ 是有界的,这意味着所有信号是有界的,所以 $\mathop {\lim }\limits_{t \to \infty } {{\boldsymbol{\hat e}}_{{k}}}\left({{t}} \right) = 0$ ,进而 $\mathop {\lim }\limits_{t \to \infty } {\boldsymbol{\hat e}}\left({{t}} \right) = 0$

    电液伺服系统主要参数值见表1

    表  1  模型参数
    符号 物理意义 数值
    ${K_{sv}}$ 伺服阀增益 16.5
    $ {\zeta _{sv}} $ 伺服阀阻尼比 0.7
    $ {\omega _{sv}}/{\text{Hz}} $ 伺服阀固有频率 386.94
    $ {K_q}/{{\text{m}}^2} \cdot {{\text{s}}^{ - 1}} $ 流量增益 7.25×10−4
    $ {K_{ce}}{\text{/}}{{\text{m}}^5} \cdot {\left( {N \cdot {\text{s}}} \right)^{ - 1}} $ 伺服阀流量–压力系数 0.008
    $ {\zeta _h} $ 液压缸阻尼比 0.23
    $ {\omega _h}/{\text{Hz}} $ 液压缸固有频率 891.72
    ${A_t}/{{\text{m}}^2}$ 液压缸有效作用面积 $2.206 \times {10^{ - 4}}$
    ${K_f}$ 位移传感器放大系数 12.71
    ${K_p}$ 比例增益 10
    ${K_i}$ 积分增益 5
    $ {K_a} $ 伺服放大器增益 0.0344

    根据本文第2节建模部分,系统可以表示为

    $$ \left\{ \begin{gathered} {\boldsymbol{\dot x}}\left( {{t}} \right) = {\boldsymbol{Ax}}\left( {{t}} \right) + {\boldsymbol{Bu}}\left( {{t}} \right) \\ {\boldsymbol{y}}\left( {{t}} \right) = {\boldsymbol{Cx}}\left( {{t}} \right) \\ \end{gathered} \right. $$

    式中: ${\boldsymbol{x}}\left( {{t}} \right) = \left[ {\ddot x,\dot x,x,{{\dot x}_v},{x_v},{u_{se}}} \right]$ ${\boldsymbol{y}}\left( {{t}} \right) = \left[ x \right]$ ${\boldsymbol{u}}\left( {{t}} \right) = \left[ {{x_d}} \right]$

    $$ {\boldsymbol{A}} = \left[ {\begin{array}{*{20}{c}} { - 2\;576}&{ - 3.136 \times {{10}^7}}&0 \\ 1&0&0 \\ 0&1&0 \\ 0&0&{ - 4.263 \times {{10}^8}} \\ 0&0&0 \\ 0&0&{ - 63.555} \end{array}} \right. \\ \left. {{\text{ }}\begin{array}{*{20}{c}} 0&{1.031 \times {{10}^8}}&0 \\ 0&0&0 \\ 0&0&0 \\ { - 3\;402}&{ - 5.905 \times {{10}^6}}&{6.708 \times {{10}^6}} \\ 1&0&0 \\ 0&0&0 \end{array}} \right] \\ $$
    $$\begin{split}{\boldsymbol{B}} = \left[ {\begin{array}{*{20}{c}} 0 \\ 0 \\ 0 \\ {4.263 \times {{10}^8}} \\ 0 \\ {63.55} \end{array}} \right]\end{split} $$
    $$ {\boldsymbol{C}}{\text{ = }}\left[ {\begin{array}{*{20}{c}} 0&0&1&0&0&0 \end{array}} \right] $$

    为了使跟踪性能更符合实际情况,本文将无故障电液伺服系统模型作为参考模型,即 ${{\boldsymbol{A}}_{{m}}} = {\boldsymbol{A}}$ ${{\boldsymbol{B}}_{{m}}} = {\boldsymbol{B}}$ ${{\boldsymbol{C}}_{{m}}} = {\boldsymbol{C}}$ 。相关控制器参数取值为: $ K=1 $ $ F=3 $ $ y=1 \times 10^{3} $ $ c_{1}=0.9 $ $ c_{2}=0.1 $ $ \lambda=0 $ $ \eta=1 $

    本文提出的重构方案主要是针对电液伺服系统的传感器故障、执行器故障和内泄漏故障。执行器故障为执行器效率损失,传感器故障为传感器恒增益漂移。所建立的模型集包括4个固定模型、1个自适应模型。

    首先,仅考虑液压缸内泄漏故障,系统参考输入为幅值为1的阶跃函数,系统响应如图3所示。

    图  3  内泄漏仿真曲线
    下载: 全尺寸图片

    图3可以看出,在存在内泄漏的情况下,电液伺服系统输出动态性能差,超调严重并伴有剧烈振荡。自适应控制器仍然能够稳定地保持系统输出渐近跟踪参考输入,瞬态误差较小,控制效果与无干扰时基本一致,有效地消除了不良因素的影响,具有较强的鲁棒性和抗扰性。

    其次,讨论所设计的自适应控制器在单一组件故障情况下的重构效果。系统参考输入设置为 $\sin (2 \text{π} t)$ ,分以下2种情况仿真验证:

    故障1情况  0 s≤t<0.75 s,无故障;0.75 s≤t<1.25 s,传感器恒增益漂移 $ {k_{{\text{DFT}}}}{\text{ = 0}}{\text{.8}} $ ;1.25 s≤t<2 s,传感器恒增益漂移 $ {k_{{\text{DFT}}}}{\text{ = 0}}{\text{.3}} $ 。仿真结果如图4所示。

    图  4  故障1系统仿真曲线
    下载: 全尺寸图片

    图4可以看出,传感器故障使系统的输出严重偏离参考输入,可能会造成更严重的后果,仅仅依靠系统基本控制器的鲁棒性无法实现渐进跟踪。加入单一模型自适应控制器后,系统性能得到了改善,输出在一段时间后可以渐进跟踪参考输入。但是,在故障出现瞬间,系统输出产生了明显的瞬态误差,这在一些安全需求较高的应用中是不被接受的。而多模型自适应控制器使得该现象得到了明显的改善,收敛速度更快,保证系统在传感器故障情况下仍能保持良好的跟踪性能与稳定性。

    故障2情况  0 s≤t<0.75 s,无故障;0.75 s≤t<1.25 s,执行器LOE ${k_{{\text{LOE}}}} = 0.7$ ;1.25 s≤t<2 s,执行器LOE ${k_{{\text{LOE}}}} = 0.{\text{3}}$ 。仿真结果如图5所示。

    图  5  故障2系统仿真曲线
    下载: 全尺寸图片

    图5可以看出,和传感器故障类似,执行器故障也会使系统的输出受到严重影响。相比于单一模型自适应控制,多模自适应控制能够更好地补偿故障带来的不利影响,从而提高系统的鲁棒性。

    最后,考虑伴有外部扰动的混合故障的情况。系统参考输入设置为 $ \sin (2{\text{π}} t) $

    混合故障如下:0 s≤t<0.25 s,无故障;0.25 s≤t<1.25 s,传感器恒增益漂移 $ {k_{{\text{DFT}}}}{\text{ = 0}}{\text{.7}} $ ,同时伴有外部扰动 $ \omega=0.1 \sin (10\; t)+0.15 $ ;1.25 s≤t<2 s,在传感器恒增益漂移的基础上,执行器LOE ${k_{{\text{LOE}}}} = 0.7$ ,同时扰动增为 $ \omega=0.3 \sin (10\; t)+0.2 $

    仿真结果如图6所示。由图6可知,在混合故障情况下,多模型自适应控制器仍然能够保持系统输出渐进跟踪参考输入,瞬态误差较小,有效地消除了不良因素的影响,具有较强的鲁棒性和抗扰性,重构控制效果令人满意。

    图  6  混合故障系统仿真曲线
    下载: 全尺寸图片

    综上,本方法满足了系统在故障条件下增加的性能需求以及安全性和可靠性的要求,能够保持系统的整体稳定性和可接受性能。

    针对电液伺服系统潜在的组件故障问题,本文提出了一种自适应控制与组合多模型相结合的重构控制方案,根据电液伺服系统的组成结构和工作原理,建立其数学模型,并对典型组件故障进行建模,设计了基于模型跟随自适应控制的重构控制器,并证明其稳定性,提出了一个覆盖故障范围的组合多模型集,通过合适的切换机制,保证输出渐进跟踪并且改善了系统的瞬态性能。对多种故障情况进行了仿真,其结果验证了该方法较强的鲁棒性、抗扰性以及良好的重构能力。未来工作将重点放在具有非线性项的电液伺服系统重构控制研究上。

  • 图  1   电液伺服系统原理

    下载: 全尺寸图片

    图  2   基于组合多模型的自适应重构系统原理

    下载: 全尺寸图片

    图  3   内泄漏仿真曲线

    下载: 全尺寸图片

    图  4   故障1系统仿真曲线

    下载: 全尺寸图片

    图  5   故障2系统仿真曲线

    下载: 全尺寸图片

    图  6   混合故障系统仿真曲线

    下载: 全尺寸图片

    表  1   模型参数

    符号 物理意义 数值
    ${K_{sv}}$ 伺服阀增益 16.5
    $ {\zeta _{sv}} $ 伺服阀阻尼比 0.7
    $ {\omega _{sv}}/{\text{Hz}} $ 伺服阀固有频率 386.94
    $ {K_q}/{{\text{m}}^2} \cdot {{\text{s}}^{ - 1}} $ 流量增益 7.25×10−4
    $ {K_{ce}}{\text{/}}{{\text{m}}^5} \cdot {\left( {N \cdot {\text{s}}} \right)^{ - 1}} $ 伺服阀流量–压力系数 0.008
    $ {\zeta _h} $ 液压缸阻尼比 0.23
    $ {\omega _h}/{\text{Hz}} $ 液压缸固有频率 891.72
    ${A_t}/{{\text{m}}^2}$ 液压缸有效作用面积 $2.206 \times {10^{ - 4}}$
    ${K_f}$ 位移传感器放大系数 12.71
    ${K_p}$ 比例增益 10
    ${K_i}$ 积分增益 5
    $ {K_a} $ 伺服放大器增益 0.0344
  • [1] 沈晨晖, 汪伟. 基于Simulink的电液伺服系统研究和故障分析[J]. 精密制造与自动化, 2018, 1(4): 11−15.
    [2] 胡良谋, 曹克强, 徐浩军. 基于LS-SVM的电液位置伺服系统多故障诊断[J]. 系统仿真学报, 2007, 19(10): 2252−2255. doi: 10.3969/j.issn.1004-731X.2007.10.027
    [3] 神英淇, 曹太强, 司国雷, 等. 电液伺服系统的模糊自适应复合控制研究[J]. 机床与液压, 2021, 49(19): 38−42.
    [4] TANG Rui, ZHANG Qi. Dynamic sliding mode control scheme for electro-hydraulic position servo system[J]. Procedia engineering, 2011, 24: 28−32. doi: 10.1016/j.proeng.2011.11.2596
    [5] YAO Jianyong, JIAO Zongxia, SHANG Yaoxing, et al. Adaptive nonlinear optimal compensation control for electrohydraulic load simulator[J]. Chinese journal of aeronautics, 2010, 23(6): 101−114.
    [6] KIM H M, PARK S H, SONG J H, et al. Robust control of electro-hydraulic actuator systems using the adaptive back-stepping control scheme[M]. Croatia: InTech, 2010.
    [7] YUAN Haibo, HONG Cheol. System identification and robust position control for electro-hydraulic servo system using hybrid model predictive control[J]. Journal of vibration and control, 2018, 24(18): 4145−4159. doi: 10.1177/1077546317721417
    [8] BECK M, SCHWUNG A, MUENCHHOF M, et al. Active fault tolerant control of an electro-hydraulic servo axis with a duplex-valve-system[J]. IFAC proceedings volumes, 2010, 43(18): 660−668. doi: 10.3182/20100913-3-US-2015.00042
    [9] KARPENKO M, SEPEHRI N. Fault-tolerant control of a servo-hydraulic positioning system with crossport leakage[J]. IEEE transactions on control systems technology, 2005, 13(1): 155−161. doi: 10.1109/TCST.2004.838570
    [10] SUN Wenli, JIAN Dandan, YUAN Yuhao, et al. Fault simulation of electro-hydraulic servo system for fault self-healing based on immune principle[C]// International Symposium on Computational Intelligence & Design. Hangzhou: IEEE, 2017.
    [11] 王大轶, 屠园园, 符方舟, 等. 航天器控制系统的自主诊断重构技术[J]. 控制理论与应用, 2019, 36(12): 1966−1978.
    [12] SHIN D H, KIM Y, Reconfigurable flight control system design using adaptive neural networks[J]. IEEE transactions on control systems technology, 2004, 12(1): 87-100.
    [13] 李俊领, 杨光红. 自适应容错控制的发展与展望[J]. 控制与决策, 2014, 29(11): 1921−1926. doi: 10.13195/j.kzyjc.2013.0905
    [14] ZHANG Zhao, YANG Zhong, XIONG Si, et al. Simple adaptive control-based reconfiguration design of cabin pressure control system[J]. Complexity, 2021(3): 1−16.
    [15] GUO Yuying, JIANG Bing. Multiple model-based adaptive reconfiguration control for actuator fault[J]. Acta automatica sinica, 2009, 35(11): 1452−1458.
    [16] 郭玉英, 姜斌, 张友民, 等. 基于多模型的飞控系统执行器故障调节[J]. 宇航学报, 2009, 30(2): 795−800. doi: 10.3873/j.issn.1000-1328.2009.02.068
    [17] TAN Chang, YAO Xuelian, TAO Gang, et al. A multiple-model based adaptive actuator failure compensation scheme for control of near-space vehicles[J]. IFAC proceedings volumes, 2012, 45(20): 594−599. doi: 10.3182/20120829-3-MX-2028.00241
    [18] CHEN Fuyang, WU Qingbo, TAO Gang, et al. A reconfiguration control scheme for a quadrotor helicopter via combined multiple models[J]. International journal of advanced robotic systems, 2014, 11(8): 122−132. doi: 10.5772/58833
    [19] CILIZ M K, TUNCAY M M. Comparative experiments with a multiple model based adaptive controller for a SCARA type direct drive manipulator[J]. Robotica, 2005, 23(6): 721−729. doi: 10.1017/S026357470500158X
    [20] SOFIANOS N A, BOUTAILS Y S, et al. Robust adaptive multiple models based fuzzy control of nonlinear systems[J]. Neurocomputing, 2016, 173: 1733−1742. doi: 10.1016/j.neucom.2015.09.047
    [21] TAN Chang, TAO Gang, QI Ruiyun. An adaptive control scheme using multiple reference models[J]. International journal of adaptive control and signal processing, 2014, 28(11): 1290−1298. doi: 10.1002/acs.2435
    [22] 司国雷, 神英淇, 王嘉磊, 等. 电液位置伺服系统自抗扰控制研究[J]. 液压与气动, 2020, 352(12): 19−26. doi: 10.11832/j.issn.1000-4858.2020.12.003
    [23] 王海杰. 电液位置伺服系统自适应滑模控制策略研究[D]. 上海: 上海交通大学, 2017.
    [24] 盛夕正. 电液位置伺服系统的设计与控制性能研究[D]. 上海: 上海交通大学, 2019.
    [25] 柴天佑, 岳恒. 自适应控制[M]. 北京: 清华大学出版社, 2016.
    [26] NARENDRA K S, BALAKRISHNAN J. Adaptive control using multiple models[J]. IEEE transactions on automatic control, 1997, 42(2): 171−187. doi: 10.1109/9.554398
    [27] 郭玉英. 基于多模型的飞机舵面故障诊断与主动容错控制[D]. 南京: 南京航空航天大学, 2009.
WeChat 点击查看大图
图(6)  /  表(1)
出版历程
  • 收稿日期:  2021-10-03
  • 网络出版日期:  2022-11-23

目录

    /

    返回文章
    返回