-
摘要: 针对不确定机械系统中普遍存在的摩擦力,由于其非线性和不确定性,传统基于摩擦模型的补偿控制方法难以达到满意的系统性能要求.本文提出基于自适应区间二型(Type-2)模糊逻辑系统对系统摩擦进行补偿建模,并在该摩擦补偿方法的基础上设计出鲁棒自适应控制器,保证系统输出精度,且对摩擦环境的变化具有较强自适应性.区间二型模糊逻辑系统相对于传统一型模糊逻辑系统具有较强的处理不确定性问题的能力,在本文中使用自适应区间二型模糊逻辑系统不断逼近摩擦力,根据李雅普诺夫稳定性理论求出自适应律并证明系统跟踪误差的有界性.在不同摩擦环境下的仿真结果验证了本文所提摩擦建模方法与控制策略的有效性与实用性.Abstract: Due to the nonlinearity and uncertainty of friction in uncertain mechanical systems, it is difficult to establish an accurate friction model. The traditional compensation control method based on friction model is difficult to meet the requirements of system performance. The adaptive interval Type-2 fuzzy logic system is proposed to model system friction, then the robust adaptive controller is designed on the basis of the friction compensation method, which ensures the output precision of the system and has strong adaptability to variational environment. Interval Type-2 fuzzy logic system has stronger ability to deal with the uncertainty problem than the traditional Type-1 fuzzy logic system, so the adaptive interval Type-2 fuzzy logic system is used to approximate friction. The adaptive law is derived and the boundedness of the tracking error is proved based on the Lyapunov stability theory. The effectiveness and practicability of the friction modeling method and control strategy are verified by simulations under different friction environments.
-
Key words:
- Nonlinear friction /
- interval Type-2 fuzzy /
- Lyapunov stability /
- adaptability
-
摩擦现象在机械动力系统中普遍存在, 随着工业化加快, 动力系统的位移精度、速度精度等有更高的要求, 例如高精度数控机床、机械手臂、电子产品等行业.由于摩擦力的非线性和不确定性, 导致系统产生输出误差, 降低系统控制精度和工作性能, 已经成为控制领域的难题[1-6].
近年来, 研究人员主要从两方面对摩擦补偿控制策略进行研究, 一方面是基于已知摩擦数学模型, 设置控制系统中控制量包含抵消摩擦项; 另一方面是在摩擦模型未知的情况下进行补偿建模控制.摩擦数学模型是根据摩擦机理的研究和实验建立具体表达方式, 但是这种数学模型只是真实摩擦力的近似表达式, 并不完全准确; 另外系统实际运行的工作环境比较复杂, 摩擦模型自适应特性较差、也会导致模型的不精确.因此不依赖摩擦模型的补偿方法受到广泛关注.文献[7]分别从摩擦模型、摩擦系统分析工具、摩擦补偿控制方法等方面进行详细阐述.从文献[7]可以看出摩擦研究在各个领域取得重大成果.但是, 近20年来, 工业控制的精度要求越来越高, 不确定性摩擦在机械系统中普遍存在, 摩擦补偿方法仍是人们所关心的问题, 研究人员继续对摩擦进行深入研究.如文献[8]针对汽车电动助力制动系统存在的摩擦现象, 建立了LuGre摩擦模型来表征系统的摩擦特性, 并采用遗传算法进行摩擦模型的参数辨识.文献[9]使用时延技术估计双框架磁悬浮控制力矩陀螺框架伺服系统存在的摩擦力矩.文献[10]利用动态增益的自适应扩展状态观测器估计磁性无杆气缸的静摩擦力.
近年来, 研究人员对基于自适应模糊逻辑系统的摩擦建模与控制方法进行了一系列研究.文献[11]采用自适应模糊逻辑系统对动力系统摩擦力进行建模, 并在此基础上构建PD (Proportion differentiation)控制器, 保证闭环系统跟踪误差的有界性.文献[12]提出用模糊状态估计器估计摩擦模型中的不可测变量, 并基于李雅普诺夫稳定性理论证明系统跟踪误差的有界性.文献[13]针对多轴伺服系统中不确定摩擦力, 采用模糊逻辑系统逼近摩擦力作为控制补偿, 并引入了全局滑模控制提高系统对参数变化和外界干扰的鲁棒性.上述文献均采用传统一型模糊逻辑系统对摩擦进行建模估计, 然而用二型模糊逻辑系统对摩擦补偿建模方面很少报道.
1975年, Zadeh教授首次提出Type-2模糊集合概念[14].在近20年来, Type-2模糊理论得到快速的发展, 国内外研究人员不断进行深入研究丰富其理论基础, 逐渐形成一个新兴的研究领域. Type-2模糊集合本质是对一型模糊集合的扩展, 对一型模糊的隶属度再次进行模糊化, 扩大隶属度空间, 可以处理语言规则不确定性问题, 能够综合同一事物的多个专家经验.文献[15-20]对Type-2模糊逻辑系统做了更深入详细的研究, 并形成完整的Type-2模糊逻辑系统理论框架.区间Type-2模糊逻辑系统具有Type-2模糊逻辑系统处理不确定问题能力, 并且避免Type-2模糊逻辑系统繁琐的计算量, 所以区间Type-2模糊逻辑系统的实际应用面比较广泛[21-22].
本文的主要贡献有以下几个方面: 1)提出基于自适应区间Type-2模糊逻辑系统对摩擦进行补偿建模; 2)在该摩擦补偿方法的基础上设计出鲁棒自适应控制器, 保证系统输出精度满足要求, 并在不同摩擦环境下仿真验证该方法的自适应性; 3)根据Lyapunov稳定性导出自适应律, 并证明跟踪误差的一致最终有界.
本文结构安排如下:第1节对传统摩擦模型进行了简要介绍.第2节对Type-2模糊集合与Type-2模糊逻辑系统进行详细介绍, 在此理论基础上使用自适应区间Type-2模糊逻辑系统对摩擦进行补偿建模, 并设计出鲁棒自适应控制器, 根据Lyapunov稳定性导出自适应律, 保证系统的输出精度.第3节在不同的摩擦环境下进行仿真, 验证本文所提方法的有效性与自适应性, 并合理设计衰减系数, 保证系统跟踪误差在适当范围之内.第4节是结论.
1. 考虑摩擦特性的不确定机械系统
参考不确定机械系统的微分方程如下:
$ \begin{equation} m\frac{\mathrm{d^{2}}x}{\mathrm{d}t^{2}} = \tau-F_{f} \end{equation} $
(1) 式中, $ m $为运动物体质量, $ x $为物体运动位移, $ \tau $为输入控制力, $ F_{f} $为不确定摩擦力.
为保证系统的控制精度, 有效方法是降低摩擦对系统影响.由于摩擦的不确定性和非线性, 摩擦建模至今仍然是一个研究难题.研究人员目前为止提出库伦摩擦模型、Stribeck摩擦模型、LuGre摩擦模型、Karnopp摩擦模型等.这里对常用摩擦模型做简要介绍[23-25]:
1) 库伦摩擦模型
该模型是一种理想干摩擦模型, 该模型表达方式简洁, 但只局限于非零速度下的摩擦, 不能描述相对速度为零时的摩擦力情况.摩擦力与作用在摩擦面的法向力成正比, 摩擦力与相对滑动速度$ v $关系如图 1(a)所示, 其数学公式为:
$ \begin{equation} F_{f} = F_{c}\cdot \mathrm{sgn} {(v)} \end{equation} $
(2) 式中, $ F_{f} $是摩擦力, $ F_{c} = \mu|F_{N}| $, $ \mu $是摩擦系数, $ F_{N} $为法向力.
2) 库伦摩擦+静摩擦模型
物体从静止加速到一定速度, 所受到的最大静摩擦力大于库伦摩擦力, 该模型相对于单一的库伦摩擦具有较高准确度, 实际上并不能真正仿真静摩擦现象.其摩擦力与相对滑动速度$ v $关系如图 1(b)所示, 其数学公式为:
$ \begin{equation} F_{f} = \left\{ \begin{array}{lll} F_{e}, \ \ \ \ \ \ \ \ \ \ \ \ \ v = 0, |F_{e}|<F_{s}\\ F_{s}\cdot \mathrm{sgn}(F_{e}), \ v = 0, |F_{e}|\geq{F_{s}}\\ F_{c}\cdot \mathrm{sgn}(v), \ \ \ v\neq0 \end{array} \right. \end{equation} $
(3) 式中, $ F_{e} $是外力, $ F_{s} $为最大静摩擦, 单位为$ \mathrm{N} $.
3) 库伦摩擦+粘滞摩擦模型
随着流体动力学的发展, 研究人员推出粘滞摩擦模型, 该模型经常与库伦摩擦配合使用, 该模型比较简单, 但该摩擦模型具有不连续性, 且无法描述摩擦的动态特性, 其控制效果受到一定程度限制.其摩擦力与相对滑动速度$ v $关系如图 1(c)所示.
4) 静摩擦+库伦摩擦+粘滞摩擦
随着粘滞摩擦理论的引入, 该模型能够改善系统在零速附近的响应, 但该模型要求准确的速度信号, 且只能描述摩擦静态特性.其摩擦力与相对滑动速度$ v $关系如图 1(d)所示.
5) Stribeck摩擦模型
该模型在描述摩擦现象时, 能够将静摩擦与粘滞摩擦之间用过渡曲线连接, 更能准确描述系统低速时摩擦特性.该摩擦模型虽然能够实现摩擦的连续性, 但不能描述摩擦力的动态特性.其摩擦力与相对滑动速度$ v $关系如图 1(e)所示, 其数学公式为:
$ \begin{align} F_{f}(v) = \, &\left(F_{c}+(F_{s}-F_{c})\mathrm{exp} \left(-\left|\frac{v}{v_{s}}\right|^{\delta_{s}}\right)\right) \cdot\\&\mathrm{sgn}(v)+F_{v}v \end{align} $
(4) 式中, $ F_{v} $为摩擦力的粘性摩擦系数, 单位为$ \mathrm{Ns/m} $; $ v_{s} $为Stribeck速度, 单位为$ \mathrm{m/s} $; $ \delta_{s} $为经验参数.
6) Karnopp摩擦模型
该模型能够比较准确描述相对速度$ v $为零时系统的摩擦特性, 且该模型对速度检测精度要求不高, 实用性比较大.但该模型只能描述摩擦模型静态特性, 不能描述动态特性.其摩擦力与相对滑动速度$ v $关系如图 1(f)所示, 其数学公式为:
$ \begin{equation} F_{f} = \left\{ \begin{array}{lll} F_{f}{(v)}, \ \ \ \ |v|\geq {\rm d}v\\ F_{e}, \ \ \ \ \ \ \ \ |v|<{\rm d}v\; \text{且}\; |F_{e}|<{F_{s}}\\ F_{s}, \ \ \ \ \ \ \ \ \text{其他} \end{array} \right. \end{equation} $
(5) 7) LuGre摩擦模型
该模型认为刚体接触面存在弹性的鬃毛, 下表面刚体的材料刚度大于上表面.该模型能够准确体现出摩擦的动态特性和静态特性, 可以更加准确描述摩擦现象, 但其参数不易确定.其数学公式为:
$ \begin{align} \left\{\begin{aligned} &F_{f} = \sigma_{0}z+\sigma_{1}\dot{z}+\sigma_{2}v\\ &\dot{z} = v-\frac{|v|z}{g(v)}\\ &\sigma_{0}g(v) = F_{c}+(F_{s}-F_{c}){\rm e}^{-|\frac{v}{v_{s}}|^{2}} \end{aligned}\right. \end{align} $
(6) 式中, $ z $为鬃毛的平均变形, $ v $为速度, $ \sigma_{0} $, $ \sigma_{1} $为动摩擦参数, $ F_{c} $, $ F_{s} $, $ \sigma_{2} $, $ v_{s} $为静摩擦参数.
综上, 摩擦模型具有一定研究成果, 但仍然存在一些问题有待解决:现提出的摩擦模型不能精确地描述摩擦的动态特性和静态特性; 模型参数不易辨识, 难以准确描述摩擦力; 由于工作环境的复杂特性, 摩擦模型未考虑其他因素的影响, 如温度、湿度、润滑条件的变化影响.
2. Ⅱ型模糊逻辑系统与控制器设计
Type-2模糊集合$ \tilde{A} $是Type-1模糊集合的扩展, 在处理复杂不确定性时具有更多的隶属度空间, 提高模糊逻辑系统对未知函数逼近的准确性, 如图 2所示. Type-2模糊集合$ \tilde{A} $表示为:
$ \begin{equation} \tilde{A} = \{((x, u), \mu_{\tilde{A}}(x, u))|x\in {X}, u\in J_{x}\subseteq[0, 1]\} \end{equation} $
(7) 式中, $ x $为二型模糊集合$ \tilde{A} $的第一变量, $ u $为二型模糊集合$ \tilde{A} $的第二变量, $ J_{x} $为$ x $的第一隶属度, $ \mu_{\tilde{A}}(x, u)) $是$ x $的第二隶属度.
二型模糊集合$ \tilde{A} $的不确定覆盖域(Footprint of uncertainty, FOU)为论域上每一点与在该点的主隶属度的笛卡尔积之并[18], 如图 2阴影部分所示, 其表达式为:
$ \begin{equation} \mathrm{FOU}(\tilde{A}) = \bigcup\limits_{x\in {X}} x\times J_{x} \end{equation} $
(8) FOU的边界由上隶属度函数(Upper membership function, UMF)和下隶属度函数(Lower membership function, LMF)构成.
区间Type-2模糊集合是Type-2模糊集合的一种特殊情况, 其第二隶属度$ \mu_{\tilde{A}}(x, u) $都等于1.区间Type-2模糊逻辑系统相对于一般的Type-2模糊逻辑系统来说计算比较简便, 更加实用, 所以在本文中使用区间Type-2模糊逻辑系统对摩擦进行建模.
以$ n $输入单输出的Mamdani区间Type-2模糊逻辑系统为例, 输入为$ x_{i}\in X_{i} $, 输出为$ y\in Y $, $ i = 1, \cdots, n $.模糊逻辑系统由下面“IF-THEN”语句构成:
$ \begin{eqnarray} R^{(j)}: {\rm{IF}}\ x_{1}\ \ {\rm is} \ \tilde{A}_{1}^{j}, \ \cdots, x_{n}\ \ {\rm is}\ \ \tilde{A}_{n}^{j}, \ \\ \ {\rm{THEN}}\ \ y \ {\rm is} \ B^{j} \end{eqnarray} $
(9) 式中, $ j = 1, \cdots, M $是规则数, $ \tilde{A}_{i}^{j} $为前件区间二型模糊集合, 与$ B^{j} $为后件一型模糊集合.在本文中, $ \pmb{x} = (x_{1}, x_{2})^{\rm{T}} = (x, \dot{x})^{\rm{T}} $, $ x $为位移信号, $ \dot{x} $为速度信号, $ y = \hat{F}_{f} $为估计摩擦力.
对于上面的模糊逻辑系统, 基于KM (Kuhn-Munkes)算法的降型在线计算过程如下:
步骤1. 计算每条规则输入向量$ \pmb{x} $的激活区间:
$ \begin{equation} \tilde{A}^{j}:\left\{ \begin{array}{lll} \tilde{A}^{j}(\pmb{x})\equiv[\underline{a}^{j}(\pmb{x}), \overline{a}^{j}(\pmb{x})]\\ \underline{a}^{j}(\pmb{x})\equiv T_{i = 1}^{n}\underline{\mu}_{\tilde{A}_{i}^{j}}(x_{i}) \\ \overline{a}^{j}(\pmb{x})\equiv T_{i = 1}^{n}\overline{\mu}_{\tilde{A}_{i}^{j}}(x_{i}) \end{array} \right. \end{equation} $
(10) 步骤2. 每条规则输出变量$ y $的最大隶属度所对应的值为:
$ \begin{equation} C_{B^{j}} = \bar{y}_{j} = \Theta_{j} \end{equation} $
(11) 步骤3. 基于KM算法降型:
$ \begin{equation} Y_{ACOS}(\pmb{x}) = \frac{1}{[y_{l}(\pmb{x}), y_{r}(\pmb{x})]} \end{equation} $
(12) 式中, $ [y_{l}(\pmb{x}), y_{r}(\pmb{x})] $为用KM算法得到的输出区间,
$ \begin{align} \left\{\begin{aligned} y_{l}(\pmb{x}) = \min\limits_{\forall \upsilon^{l}_{j}\in [\underline{a}^{j}(\pmb{x}), \overline{a}^{j}(\pmb{x})]} \frac{\sum\limits_{j = 1}^{M}\Theta_{j}\upsilon^{l}_{j}} {\sum\limits_{j = 1}^{M}\upsilon^{l}_{j}} \\ y_{r}(\pmb{x}) = \max\limits_{\forall \upsilon^{r}_{j}\in [\underline{a}^{j}(\pmb{x}), \overline{a}^{j}(\pmb{x})]} \frac{\sum\limits_{j = 1}^{M}\Theta_{j}\upsilon^{r}_{j}} {\sum\limits_{j = 1}^{ M}\upsilon^{r}_{j}} \end{aligned}\right. \end{align} $
(13) 步骤4. 计算区间Type-2模糊逻辑系统输出:
$ \begin{equation} y = \frac{1}{2}(y_{l}(\pmb{x})+y_{r}(\pmb{x})) = \frac{1}{2}\pmb{\Theta}^{\rm{T}}(\pmb{\xi}_{l}(\pmb{x})+\pmb{\xi}_{r}(\pmb{x})) \end{equation} $
(14) 式中, $ \pmb{\Theta} = [{\Theta}_{1}, \cdots, {\Theta}_{\mathrm{M}}]^{\rm{T}} $为区间Type-2模糊逻辑系统的自适应参数向量, $ \pmb{\xi}_{l}(\pmb{x}) $为区间Type-2模糊逻辑系统的左基函数, $ \pmb{\xi}_{r}(\pmb{x}) $为区间Type-2模糊逻辑系统的右基函数, 表达式如下:
$ \begin{align} \left\{\begin{aligned} \pmb{\xi}_{l}(\pmb{x}) = \left[\frac{\upsilon^{l}_{1}}{\sum\limits_{j = 1}^{ M}\upsilon^{l}_{j}}, \cdots, \frac{\upsilon^{l}_{ M}}{\sum\limits_{j = 1}^{ M}\upsilon^{l}_{j}}\right]^{\rm{T}}\\ \pmb{\xi}_{r}(\pmb{x}) = \left[\frac{\upsilon^{r}_{1}}{\sum\limits_{j = 1}^{ M}\upsilon^{r}_{j}}, \cdots, \frac{\upsilon^{r}_{ M}}{\sum\limits_{j = 1}^{ M}\upsilon^{r}_{j}}\right]^{\rm{T}} \end{aligned}\right. \end{align} $
(15) 系统控制目标为系统的输出位移追踪目标位移$ x_{r} $, 假设真实摩擦力$ F_{f} $已知, 设置控制器为:
$ \begin{equation} \tau = m(x_{r}^{(2)}+\pmb{k}^{\rm{T}}\pmb{e})+{F}_{f} \end{equation} $
(16) 式中, $ e = x_{r}-x $, $ \pmb{e} = [e, \dot{e}]^{\rm{T}} $, $ \pmb k = [k_{2}, k_{1}]^{\rm{T}} $.
将式(16)代入式(1)可得跟踪误差动态方程为
$ \begin{equation} \ddot{e}+k_{1}\dot{e}+k_{2}e = 0 \end{equation} $
(17) 通过选择$ k_{1} $和$ k_{2} $值, 使多项式$ h(s) = s^{2}+k_{1}s+k_{2} $所有根都在左半开平面上, 满足Hurwitz条件, 跟踪误差$ \pmb e $趋近于0.
由于真实摩擦力$ F_{f} $是未知的, 本文采用区间Type-2模糊逻辑系统$ \hat{F}_{f} $估计摩擦力值, 可表示为:
$ \begin{equation} \hat{F}_{f}(\pmb{x}|\pmb{\Theta}) = \frac{1}{2}\pmb{\Theta}^{\rm{T}}(\pmb{\xi}_{l}(\pmb{x})+\pmb{\xi}_{r}(\pmb{x})) \end{equation} $
(18) 并在控制器(16)中引入补偿项$ \tau_{a} $来抵消模糊逻辑系统建模误差, 构建鲁棒自适应控制器:
$ \begin{equation} \tau = m(x_{r}^{(2)}+\pmb{k}^{\rm{T}}\pmb{e})+\hat{F}_{f}-\tau_{a} \end{equation} $
(19) 将式(19)代入式(1)中, 可得跟踪误差动态方程:
$ \begin{equation} \dot{\pmb{e}} = A\pmb{e}+B\tau_{a}+B(F_{f}-\hat{F}_{f}) \end{equation} $
(20) 式中, $ A = \left[ {array}{cc} 0 & 1\\ -k_{2} & -k_{1} \\ {array} \right], B = \left[ {array}{cc} 0\\ \dfrac{1}{m} \\ {array} \right] $
设置最优参数$ \pmb{\Theta}^{\ast} $:
$ \begin{equation} \pmb{\Theta}^{\ast} = \mathrm{arg}\min\limits_{\pmb{\Theta}\in\Omega_{\pmb\Theta}}(\mathrm{sup}_{\pmb{x}\in \Omega_{\pmb x}}\parallel{F_{f}}(\pmb{x})-\hat{F}_{f}(\pmb{x}|\pmb{\Theta})\parallel) \end{equation} $
(21) 式中, $ \Omega_{\pmb\Theta}: = \{{\pmb{\Theta}|\parallel\pmb{\Theta}\parallel\leq}M_{\pmb\Theta}\}, \Omega_{\pmb x}: = \{\pmb{x}|\parallel\pmb{x}\parallel\leq\ M_{\pmb x} $$ \} $, $ M_{\pmb\Theta} $, $ M_{\pmb x} $是正数.
模糊逻辑系统最小逼近误差可表示为:
$ \begin{equation} \omega = F_{f}(\pmb{x})-\hat{F}_{f}(\pmb{x}|\pmb{\Theta}^{\ast}) \end{equation} $
(22) 跟踪误差动态方程(20)可以写为:
$ \begin{equation} \dot{\pmb{e}} = A\pmb{e}+B\tau_{a}+B(\hat{F}_{f}(\pmb{x}|\pmb{\Theta}^{\ast})-\hat{F}_{f}(\pmb{x}|\pmb{\Theta}))+B\omega \end{equation} $
(23) 经过进一步计算上式可以简化为:
$ \begin{equation} \dot{\pmb{e}} = A\pmb{e}+B\tau_{a}+\frac{1}{2}B\tilde{\pmb {\Theta}}^{\rm{T}}(\pmb{\xi}_{l}(\pmb{x})+\pmb{\xi}_{r}(\pmb{x}))+B\omega \end{equation} $
(24) 式中, $ \tilde{\pmb {\Theta}} = \pmb{\Theta^{\ast}}-\pmb{\Theta} $.
定理1. 对于不确定机械系统(1), 采用鲁棒自适应控制器为式(19), 其中补偿项$ \tau_{a} $为:
$ \begin{equation} \tau_{a} = -\frac{1}{\lambda}B^{\rm{T}}P\pmb{e} \end{equation} $
(25) 二型模糊逻辑系统的自适应律为
$ \begin{align} \dot{\pmb{\Theta}} = \left\{\begin{aligned} & \frac{\gamma}{2}\pmb{e}^{\rm{T}}PB(\pmb{\xi}_{l}(\pmb{x})+ \pmb{\xi}_{r}(\pmb{x})), \\ & \ \ \ \ \ \ \ \ \ {\rm{if}} \; ( \| \pmb{\Theta}\|<M_{\pmb\Theta}) \ {\rm{or}}\; \{\|\pmb{\Theta}\| = M_{\pmb\Theta}\\ & \ \ \ \ \ \ \ \ \ {\rm{and}}\; \ B^{\rm{T}}P\pmb{e}(\pmb{\xi}_{l}(\pmb{x})+\pmb{\xi}_{r}(\pmb{x}))^{\rm{T}}\pmb{\Theta}\leq0\}\\ &P[\cdot], \ \ {\rm{if}} \ (\| \pmb{\Theta}\| = M_{\pmb\Theta}) \ {\rm{and}}\; B^{\rm{T}}P\pmb{e}(\pmb{\xi}_{l}(\pmb{x}) + \\ & \ \ \ \ \ \ \ \ \ \ \pmb{\xi}_{r}(\pmb{x}))^{\rm{T}} \pmb{\Theta}>0 \end{aligned}\right. \end{align} $
(26) 式中, $ \gamma $是给定正常数, $ P[\cdot] $投影算子可表示为:
$ \begin{align} P[\cdot] = \, & \frac{\gamma}{2}\pmb{e}^{\rm{T}}PB (\pmb{\xi}_{l}(\pmb{x})+\pmb{\xi}_{r}(\pmb{x}))- \\ & \frac{\gamma}{2}\frac{B^{\rm{T}}P\pmb e(\pmb{\xi}_{l}(\pmb{x})+\pmb{\xi}_{r}(\pmb{x}))^{\rm{T}}\pmb \Theta}{\|\pmb \Theta\|^{2}}\pmb \Theta \end{align} $
(27) 矩阵$ P = {P}^{\rm{T}}>0 $是Riccati方程的解, Riccati方程表示为:
$ \begin{align} PA+A^{\rm{T}}P+ Q-\frac{2}{\lambda}PBB^{\rm{T}}P +\frac{1}{\alpha^{2}}PBB^{\rm{T}}P = 0 \end{align} $
(28) 式中, $ Q $为正定对称矩阵, $ \lambda>0 $, $ \alpha $为衰减系数.则可以保证系统的稳定性和跟踪误差有界性.
证明. 考虑Lyapunov函数:
$ \begin{align} V = \frac{1}{2}\pmb{e}^{\rm{T}} P\pmb{e}+\frac{1}{2\gamma}\tilde{\pmb{\Theta}}^{\rm{T}}\tilde{\pmb{\Theta}} \end{align} $
(29) 则Lyapunov函数的时间导数$ {\dot{V}} $可以表示为:
$ \begin{align} \dot{V} = \frac{1}{2}\dot{\pmb{e}}^{\rm{T}} P\pmb{e}+\frac{1}{2}\pmb{e}^{\rm{T}} P\dot{\pmb{e}}+\frac{1}{\gamma}\dot{\tilde{\pmb{\Theta}}}^{\rm{T}}\tilde{\pmb{\Theta}} \end{align} $
(30) 将式(24)和式(25)代入上式得到:
$ \begin{align} \dot{V} = \, & \frac{1}{2}\Bigg[{\pmb{e}}^{\rm{T}} A^{\rm{T}}{P}\pmb e-\frac{1}{\lambda}{\pmb{e}}^{\rm{T}}{P}B B^{\rm{T}}P{\pmb{e}}+ \frac{1}{2}(\pmb{\xi}_{l}(\pmb{x})+ \\ & \pmb{\xi}_{r}(\pmb{x}))^{\rm{T}}\tilde{\pmb{\Theta }} B^{\rm{T}}P{\pmb{e}}+\omega^{\rm{T}} B^{\rm{T}}P{\pmb{e}}+{\pmb{e}}^{ \mathrm{T}}{ P}{ A}{\pmb{e}} - \\ &\frac{1}{\lambda}{\pmb{e}}^{\rm{T}}{P}B B^{\rm{T}}P{\pmb{e}}+\pmb{e}^{\rm{T}}{P}B\omega+ \\ & \frac{1}{2}{\pmb{e}}^{\rm{T}}{P} B\tilde{\pmb{\Theta }}^{\rm{T}}(\pmb{\xi}_{l}(\pmb{x})+\pmb{\xi}_{r}(\pmb{x}))\Bigg]+\frac{1}{\gamma}\dot{\tilde{\pmb{\Theta}}}^{\rm{T}}\tilde{\pmb{\Theta}} \end{align} $
(31) 由式(28)代入式(31), 并经过一系列计算之后得到:
$ \begin{align} \dot{V} = \, & -\frac{1}{2}\pmb{e}^{\rm{T}}Q\pmb{e} -\frac{1}{2\alpha^{2}}\pmb{e}^{\rm{T}}PBB^{\rm{T}}P\pmb{e} +\\ &\frac{1}{2}(\pmb e^{\rm{T}}P B\omega+\omega^{\rm{T}}B^{\rm{T}} P\pmb{e})+\\ &\left[\frac{\gamma}{2}(\pmb{\xi}_{l}(\pmb{x}) +\pmb{\xi}_{r}(\pmb{x}))^{\rm{T}}B^{\rm{T}} P\pmb{e}+\dot{\tilde{\pmb{\Theta }}}^{\rm{T}}\right]\frac{\tilde{\pmb{\Theta}}}{\gamma} \end{align} $
(32) 由于$ \dot{\tilde{\pmb{\Theta}}} = -\dot{{\pmb{\Theta}}} $, 根据自适应律(26), 上式可以写为:
$ \begin{align} \dot{V} = \, & -\frac{1}{2}\pmb{e}^{\rm{T}}Q\pmb{e} -\frac{1}{2}\left(\frac{1}{\alpha}B^{\rm{T}}P\pmb{e}- \alpha\omega\right)^{\rm{T}}\times\\ & \left(\frac{1}{\alpha}B^{\rm{T}}P\pmb{e}-\alpha\omega\right)+\frac{1}{2}\alpha^{2}\omega^{\rm{T}}\omega \end{align} $
(33) 上式可推导为:
$ \begin{align} \dot{V}\leq-\frac{1}{2}\pmb{e}^{\rm{T}}Q\pmb{e}+\frac{1}{2}\alpha^{2}\omega^{\rm{T}}\omega. \end{align} $
(34) 由于矩阵$ Q $为正定矩阵, 特征值均大于零, 则上式可推导为:
$ \begin{align} \dot{V}\leq-\frac{1}{2}\delta_{\rm min}(Q)\|\pmb{e}\|^2+ \frac{1}{2}\alpha^{2}\omega_{0}^{2} \end{align} $
(35) 式中, $ \delta_{\rm min}(Q) $为矩阵$ Q $的最小特征值, $ \omega_{0} $为逼近误差$ \omega $的界, 即$ \|\omega\|\leq\omega_{0} $.
可以看出式(35)为变号函数, 当在约束集$ \Gamma $之外时, 则$ \dot{{ V}}<0 $
$ \begin{align} \Gamma = \left\{\pmb{e} | 0\leq||\pmb{e}||\leq\frac{\alpha\omega_{0}}{\sqrt{\delta_{\rm min}( Q)}}\right\} \end{align} $
(36) 由标准Lyapunov理论可知, 跟踪误差$ \pmb e $一致最终有界(Uniformly ultimately bounded, UUB).
3. 仿真与结果分析
仿真采用机械动力系统模型为:
$ \begin{equation} m\frac{\mathrm{d^{2}}x}{\mathrm{d}t^{2}} = \tau-F_{f} \end{equation} $
(37) 步骤1. 建立二型模糊逻辑系统
摩擦力$ F_{f} $由区间Type-2模糊逻辑系统估计, 输入向量$ \pmb{x} = (x_{1}, x_{2})^{\rm{T}} = (x, \dot{x})^{\rm{T}} $, 模糊规则为:
$ \begin{align} R^{(j)}: {\rm{IF}}\ x_{1}\ \ {\rm is} \ \tilde{A}_{1}^{j}, x_{2}\ \ {\rm is}\ \ \tilde{A}_{2}^{j}, \\ \ \ {\rm{THEN}}\ \ \hat{F}_{f} \ {\rm is}\ B^{j} \end{align} $
(38) 式中, $ j = 1, \cdots, 6 $, $ \tilde{A}_{1}^{j}, \tilde{A}_{2}^{j} $分别是$ x_{1}, x_{2} $的区间二型模糊集合, 如图 3所示, $ B^{j} $是$ \hat{F}_{f} $的传统一型模糊集合.
设置模糊规则前件隶属度函数为:
$ \begin{align} \left\{\begin{aligned} \overline{\mu}_{\tilde{A}_{i}^{1}}(x_{i}) = \, &\frac{1}{1+\mathrm{exp}(5\times(x_{i}+2))}\\ \underline{\mu}_{\tilde{A}_{i}^{1}}(x_{i}) = \, &\frac{1}{1+\mathrm{exp}(5\times(x_{i}+2.2))} \end{aligned}\right. \end{align} $
(39) $ \begin{align} \left\{\begin{aligned} \overline{\mu}_{\tilde{A}_{i}^{2}}(x_{i}) = \, &\mathrm{exp}(-(x_{i}+1.5)^{2})\\ \underline{\mu}_{\tilde{A}_{i}^{2}}(x_{i}) = \, &\mathrm{exp}(-2\times(x_{i}+1.5)^2) \end{aligned}\right. \end{align} $
(40) $ \begin{align} \left\{\begin{aligned} \overline{\mu}_{\tilde{A}_{i}^{3}}(x_{i}) = \, &\mathrm{exp}(-(x_{i}+0.5)^{2})\\ \underline{\mu}_{\tilde{A}_{i}^{3}}(x_{i}) = \, &\mathrm{exp}(-2\times(x_{i}+0.5)^2) \end{aligned}\right. \end{align} $
(41) $ \begin{align} \left\{\begin{aligned} \overline{\mu}_{\tilde{A}_{i}^{4}}(x_{i}) = \, &\mathrm{exp}(-(x_{i}-0.5)^{2})\\ \underline{\mu}_{\tilde{A}_{i}^{4}}(x_{i}) = \, &\mathrm{exp}(-2\times(x_{i}-0.5)^2) \end{aligned}\right. \end{align} $
(42) $ \begin{align} \left\{\begin{aligned} \overline{\mu}_{\tilde{A}_{i}^{5}}(x_{i}) = \, &\mathrm{exp}(-(x_{i}-1.5)^{2})\\ \underline{\mu}_{\tilde{A}_{i}^{5}}(x_{i}) = \, &\mathrm{exp}(-2\times(x_{i}-1.5)^2) \end{aligned}\right. \end{align} $
(43) $ \begin{align} \left\{\begin{aligned} \overline{\mu}_{\tilde{A}_{i}^{6}}(x_{i}) = \, &\frac{1}{1+\mathrm{exp}(-5\times(x_{i}-2))}\\ \underline{\mu}_{\tilde{A}_{i}^{6}}(x_{i}) = \, &\frac{1}{1+\mathrm{exp}(-5\times(x_{i}-2.2))} \end{aligned}\right. \end{align} $
(44) 根据第2节中计算步骤可建立区间二型模糊逻辑系统$ \hat{F}_{f} $.
步骤2. 参数选择
在仿真中, 选择仿真环境参数: $ m = 1 \, \mathrm{kg} $, $ F_{N} = 5\, \mathrm{N} $, $ F_{v} = 0.04F_{N} $, $ v_{s} = 1 $, $ \delta_{s} = 1 $.系统状态初始值$ x_{1}(0) = 0.7 $, $ x_{2}(0) = 0.5 $.选择控制器参数: $ \gamma = 0.1 $, $ k_{1} = 2 $, $ k_{2} = 1 $, 自适应参数初始值$ \pmb \Theta(0) = [-8, -4, -1, 2, 3, 4]^{\rm{T}} $.目标位移$ x_{r} = \mathrm{sin}(t) $, 目标速度$ \dot{x}_{r} = \mathrm{cos}(t) $.取$ Q = \left[ {array}{cc} 10 & 0\\ 0 & 10 \\ {array} \right] $, 考虑两种衰减系数情况$ \alpha = 0.2, 0.4, $设置$ \lambda = 0.08, 0.32 $, 通过解Riccati方程(28)得$ P = \left[ {array}{cc} 15 & 5\\ 5 & 5 \\ {array} \right] $.
求得鲁棒自适应控制器$ \tau $为:
$ \begin{equation} \tau = (-\mathrm{sin}(t)+2\dot{e}+e)+\hat{F}_{f}-\tau_{a} \end{equation} $
(45) 求得自适应律为:
$ \begin{align} \dot{\pmb{\Theta}} = \frac{1}{2}\gamma(5e+5\dot{e})(\pmb{\xi}_{l}(\pmb{x})+\pmb{\xi}_{r}(\pmb{x})) \end{align} $
(46) 求得补偿$ \tau_{a} $为:
$ \begin{align} \tau_{a} = -\frac{1}{\lambda}(5e+5\dot{e}) \end{align} $
(47) 步骤3. 仿真环境设置
为验证本文所提建模方法的有效性与自适应性, 设置仿真环境如表 1所示.环境1为Stribeck摩擦模型, 不同时间段参数发生变化, 但摩擦结构不变; 环境2中摩擦力在不同时间段不仅摩擦力参数发生变化, 摩擦结构也发生变化.
表 1 摩擦模型参数与结构Table 1 Friction model parameters and structures仿真环境 摩擦模型Ff(t/s) 0<t≤20 20<t≤40 40<t≤60 环境1 Stribeck Stribeck Stribeck Fc=0.40FN Fc=0.30FN Fc=0.20FN Fs=0.60FN Fs=0.50FN Fs=0.40FN 环境2 Stribeck Stribeck Coulomb Fc=0.40FN Fc=0.30FN Fc=0.20FN Fs=0.60FN Fs=0.50FN 步骤4. 仿真结果与分析
为验证基于区间Type-2模糊逻辑系统对摩擦进行建模补偿的控制效果, 在表 1所示的两种仿真环境下, 对不确定机械系统进行仿真, 系统位移与速度追踪情况如图 4、图 5所示.
图 4和图 5是分别是在两种摩擦环境下取不同的衰减系数时系统输出位移与速度追踪曲线, 可以看出本文提出摩擦建模方法和鲁棒自适应控制器能够使不确定机械系统输出的实际位移$ x_{1}(t) $与速度$ x_{2}(t) $准确追踪目标位移$ \mathrm{sin}(t) $与目标速度$ \mathrm{cos}(t) $, 追踪误差在一定范围内, 且具有较强的自适应性.从数学公式(36)与仿真结果两方面可以看出系统跟踪误差受衰减系数影响, 且成正比关系, 所以在设置控制器时应根据实际情况选取合适衰减系数.
注1. 在实际工程应用中, 衰减系数过小, 会导致控制变量剧烈变化, 对系统执行器不利.
4. 结论
针对传统摩擦模型建模不准确、参数难以辨识以及不具有自适应特性等问题, 本文运用区间Type-2模糊逻辑系统对不确定机械系统的摩擦力进行补偿建模, 自适应律根据李雅普诺夫稳定性理论求出.设计鲁棒自适应控制器使系统输出具有较好追踪性能, 根据李雅普诺夫理论证明该系统跟踪误差一致最终有界.经过数值仿真证明:本文所提方法能够适用于不同的摩擦环境, 自适应能力较强, 具有准确的追踪性能, 且通过设置合理的衰减系数, 保证系统跟踪误差满足精度要求.
-
表 1 摩擦模型参数与结构
Table 1 Friction model parameters and structures
仿真环境 摩擦模型Ff(t/s) 0<t≤20 20<t≤40 40<t≤60 环境1 Stribeck Stribeck Stribeck Fc=0.40FN Fc=0.30FN Fc=0.20FN Fs=0.60FN Fs=0.50FN Fs=0.40FN 环境2 Stribeck Stribeck Coulomb Fc=0.40FN Fc=0.30FN Fc=0.20FN Fs=0.60FN Fs=0.50FN -
[1] Bisoffi A, Lio M D, Teel A R, Zaccarian L. Global asymptotic stability of a PID control system with coulomb friction. IEEE Transactions on Automatic Control, 2018, 63(8):2654-2661 doi: 10.1109/TAC.2017.2774443 [2] Kim S. Moment of inertia and friction torque coefficient identification in a servo drive system. IEEE Transactions on Industrial Electronics, 2019, 66(1):60-70 doi: 10.1109/TIE.41 [3] 刘强, 尔联洁, 刘金琨.参数不确定机械伺服系统的鲁棒非线性摩擦补偿控制.自动化学报, 2003, 29(4):628-632 http://www.aas.net.cn/CN/abstract/abstract13916.shtmlLiu Qiang, Er Lian-Jie, Liu Jin-Kun. Robust nonlinear friction compensation of mechanical servo system with time variable parameters. Acta Automatica Sinica, 2003, 29(4):628-632 http://www.aas.net.cn/CN/abstract/abstract13916.shtml [4] 丛爽, De Carli Alessandro.两种补偿动态摩擦力的先进控制策略.自动化学报, 1998, 24(2):236-240 http://www.aas.net.cn/CN/abstract/abstract16891.shtmlCong Shuang, De Carli-Alessandro. Two advanced control strategies for dynamic friction compensation. Acta Automatica Sinica, 1998, 24(2):236-240 http://www.aas.net.cn/CN/abstract/abstract16891.shtml [5] Verbert K A J, Toth R, Babuska R. Adaptive friction compensation:a globally stable approach. IEEE/ASME Transactions on Mechatronics, 2016, 21(1):351-363 http://d.old.wanfangdata.com.cn/NSTLQK/NSTL_QKJJ0234015635/ [6] Roy R, Wang L, Simaan N. Modeling and estimation of friction, extension, and coupling effects in multisegment continuum robots. IEEE/ASME Transactions on Mechatronics, 2017, 22(2):909-920 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=7a4445d072724a96787132c1742ba9a7 [7] Armstrong-Helouvry B, Dupont P, Wit C C D. A survey of models, analysis tools and compensation methods for the control of machines with friction. Automatica, 1994, 30(7):1083-1138 doi: 10.1016/0005-1098(94)90209-7 [8] 何睿, 吴坚, 高吉.汽车电动助力制动系统摩擦建模与补偿控制.汽车工程, 2017, 39(6):683-688 http://d.old.wanfangdata.com.cn/Periodical/qcgc201706012He Rui, Wu Jian, Gao Ji. Modeling and compensation control for friction in vehicle power assisted braking system, Automotive Engineering, 2017, 39(6):683-688 http://d.old.wanfangdata.com.cn/Periodical/qcgc201706012 [9] Cui P L, Zhang D C, Yang S, Li H T. Friction compensation based on time delay control and internal model control for gimbal system in MSCMG. IEEE Transactions on Industrial Electronics, 2017, 64(5):3798-3807 doi: 10.1109/TIE.2016.2644620 [10] Yang H J, Sun J H, Xia Y Q, Zhao L. Position control for magnetic rodless cylinders with strong static friction. IEEE Transactions on Industrial Electronics, 2018, 65(7):5806-5815 doi: 10.1109/TIE.2017.2782198 [11] 王永富, 王殿辉, 柴天佑.基于数据挖掘与系统理论建立摩擦模糊模型与控制补偿.自动化学报, 2010, 36(3):412-420 http://www.aas.net.cn/CN/abstract/abstract13682.shtmlWang Yong-Fu, Wang Dian-Hui, Chai Tian-You. Data mining and systems theory based fuzzy modeling and control compensation for friction. Acta Automatica Sinica, 2010, 36(3):412-420 http://www.aas.net.cn/CN/abstract/abstract13682.shtml [12] 王永富, 王殿辉, 柴天佑.基于状态估计的摩擦模糊建模与鲁棒自适应控制.自动化学报, 2011, 37(2):245-252 http://www.aas.net.cn/CN/abstract/abstract17430.shtmlWang Yong-Fu, Wang Dian-Hui, Chai Tian-You. State estimate-based friction fuzzy modeling and robust adaptive control. Acta Automatica Sinica, 2011, 37(2):245-252 http://www.aas.net.cn/CN/abstract/abstract17430.shtml [13] Zhong G L, Shao Z Z, Deng H, Ren J L. Precise position synchronous control for multi-axis servo systems. IEEE Transactions on Industrial Electronics, 2017, 64(5):3707-3717 doi: 10.1109/TIE.2017.2652343 [14] Zadeh L A. The Concept of a Linguistic Variable and its Application to Approximate Reasoning-I. Information Sciences, 1975, 8(3):199-249 doi: 10.1016/0020-0255(75)90036-5 [15] 王飞跃, 莫红.关于二型模糊集合的一些基本问题.自动化学报, 2017, 43(16):1114-1141 http://www.aas.net.cn/CN/abstract/abstract19087.shtmlWang Fei-Yue, Mo Hong. Some fundamental issues on type-2 fuzzy sets. Acta Automatica Sinica, 2017, 43(16):1114-1141 http://www.aas.net.cn/CN/abstract/abstract19087.shtml [16] Wang Jia-Jun. A new type of fuzzy membership function designed for interval type-2 fuzzy neural network. Acta Automatica Sinica, 2017, 43(8):1425-1433 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=zdhxb201708014 [17] 莫红, 王飞跃, 肖志权, 陈茜.基于区间二型模糊集合的语言动力系统稳定性.自动化学报, 2011, 37(8):1018-1024 http://www.aas.net.cn/CN/abstract/abstract17522.shtmlMo Hong, Wang Fei-Yue, Xiao Zhi-Quan, Chen Qian.Stabilities of linguistic dynamic systems based on interval type-2 fuzzy sets. Acta Automatica Sinica, 2011, 37(8):1018-1024 http://www.aas.net.cn/CN/abstract/abstract17522.shtml [18] Mo Hong, Wang Fei-Yue, Zhou Min, Li R M, Xiao Z Q. Footprint of uncertainty for type-2 fuzzy sets. Information Sciences, 2014, 272:96-110 doi: 10.1016/j.ins.2014.02.092 [19] Mendel J M. On KM Algorithms for solving type-2 fuzzy set problems. IEEE Transactions on Fuzzy Systems, 2013, 21(3):426-446 doi: 10.1109/TFUZZ.2012.2227488 [20] Mendel J M. General type-2 fuzzy logic systems made simple:a tutorial. IEEE Transactions on Fuzzy Systems, 2014, 22(5):1162-1182 doi: 10.1109/TFUZZ.2013.2286414 [21] 李润梅, 梁秋鸿.基于区间二型模糊集合的人工交通系统可信度评估.自动化学报, DOI: 10.16383/j.aas.c180105Li Run-Mei, Wang Qiu-Hong. Artiflcial traffic system credibility evaluation based on interval type-2 fuzzy sets. Acta Automatica Sinica, DOI: 10.16383/j.aas.c180105 [22] 唐晓铭, 邓梨, 虞继敏, 屈洪春.基于区间二型T-S模糊模型的网络控制系统的输出反馈预测控制.自动化学报, DOI: 10.16383/j.aas.c170554Tang Xiao-Ming, Deng Li, Yu Ji-Min, Qu Hong-Chun. Output feedback model predictive control for interval type-2 T-S fuzzy networked control systems with bounded sisturbance. Acta Automatica Sinica, DOI: 10.16383/j.aas.c170554 [23] 丁千, 翟红梅.机械系统摩擦动力学研究进展.力学进展, 2013, 43(1):112-131 http://d.old.wanfangdata.com.cn/Conference/8280928Ding Qian, Zhai Hong-Mei. The adavance in resarches of friction dynamics in mechanics system, Advances in Mechanics, 2013, 43(1):112-131 http://d.old.wanfangdata.com.cn/Conference/8280928 [24] 张新刚, 基于扩展Stribeck效应的摩擦实验建模及系统动力学研究[博士学位论文], 上海交通大学, 中国, 2009 http://cdmd.cnki.com.cn/Article/CDMD-10248-2010033193.htmZhang X G. Researches of Experimental Modeling and System Dynamics on Frictions Concerning Extended Stribeck Effect[Ph.D. thesis], Shanghai Jiao Tong University, China, 2009 http://cdmd.cnki.com.cn/Article/CDMD-10248-2010033193.htm [25] 向红标, 开放式伺服系统的摩擦建模与补偿研究[博士学位论文], 天津大学, 中国, 2010 http://cdmd.cnki.com.cn/Article/CDMD-10056-1011266411.htmXiang H B. Research on Model and Compensation for Open Servo System[Ph.D. thesis], Tianjin University, China, 2010 http://cdmd.cnki.com.cn/Article/CDMD-10056-1011266411.htm 期刊类型引用(7)
1. 王骏骋,吕林峰,王法慧. 基于区间二型模糊逻辑的制动防抱死控制. 交通运输工程学报. 2024(03): 238-250 . 百度学术
2. 付佳倩,李珊瑚,胡俊宇. 直线伺服系统位移跟踪误差的分析与迭代补偿. 动力学与控制学报. 2024(12): 87-95 . 百度学术
3. 汤健,夏恒,余文,乔俊飞. 城市固废焚烧过程智能优化控制研究现状与展望. 自动化学报. 2023(10): 2019-2059 . 本站查看
4. 黄金,陈纪旸,霍江航,马晓静. 无人飞艇纵向L1自适应控制. 计算机仿真. 2023(09): 1-6+255 . 百度学术
5. 张燕,杨安杰,孙善乐,李璇,李小觅,孙慧. 基于IT2-PCM的上肢外骨骼系统T-S模糊建模与预测控制. 兰州大学学报(自然科学版). 2023(05): 576-584 . 百度学术
6. 程国扬,刘阳阳. 带摩擦补偿的运动控制探索型实验设计. 实验室研究与探索. 2022(03): 8-13 . 百度学术
7. 石暕,程国扬. 基于自适应摩擦补偿的复合非线性轨迹跟踪控制. 电机与控制应用. 2021(09): 14-20 . 百度学术
其他类型引用(14)
-