基于双参数奇异摄动的液压刚柔机械臂控制方法

摘要

基于双参数奇异摄动的液压刚柔机械臂控制方法涉及智能机械与机器人控制技术领域，该方法通过双参数奇异摄动原理，将液压刚柔机械臂系统降阶为表征液压伺服驱动的快变子系统、表征弹性振动的次快变子系统及表征大范围刚性运动的慢变子系统；其中，快变子系统控制模块采用自适应滑模控制器，次快变子系统控制模块采用最优控制器抑制弹性振动，慢变子系统控制模块采用二阶滑模控制器实现关节轨迹跟踪控制。本发明的控制方法基于双参数奇异摄动将复杂的液压刚柔机械臂系统模型进行降阶，设计多时标子系统控制器，使得控制方法大大简化，且更具实用性。

主权项

1.基于双参数奇异摄动的液压刚柔机械臂控制方法，其特征在于，该方法包括如下步骤：第一步、模型分解：三自由度液压刚柔机械臂动力学方程如式(13)所示：$\left\{\begin{array}{c}M(\theta ,q)\left[\begin{array}{c}\stackrel{··}{\theta }\\ \stackrel{··}{q}\end{array}\right]+K\left[\begin{array}{c}\theta \\ q\end{array}\right]+G(\theta ,\stackrel{·}{\theta },q,\stackrel{·}{q})=\left[\begin{array}{c}\tau \\ 0\end{array}\right]\\ \stackrel{·}{\tau }+\mathrm{A\tau }+B\stackrel{·}{\theta }=\mathrm{CI}\end{array}\right.---\left(13\right)$其中，M(θ,q)∈R^5×5为对称、正定的惯性矩阵；θ=[θ₁ θ₂ θ₃]^T为关节角；q=[q₁ q₂]^T为模态坐标；为包含哥氏力、离心力、重力的非线性项；柔性臂的刚度矩阵为K=diag(0,0,0,k₁,k₂)；作用在关节上的广义力矩为τ=[τ₁ τ₂ τ₃]^T；I=[i₁ i₂ i₃]^T为伺服阀控制电流；A=diag(a₁,a₂,a₃)，a₁=4β_e(C_tm1+K_c1)V_t1，B=diag(b₁,b₂,b₃)，b₁=4β_eD²V_t1，C=diag(c₁,c₂,c₃)，c₁=4β_eDK_q1K_i1V_t1,c_i=4β_eA_piDK_qiK_iJ_iV_ti,i=2,3；采用双参数奇异摄动技术将式(13)分解，取第一个小参数并且满足0<ε₁<<1，并在边界层上引入第一级快变时标得到液压刚柔机械臂的第一级快变子系统动力学方程如式(20)$\frac{{\mathrm{d\tau }}_{f1}}{{\mathrm{d\sigma }}_1}=-\tilde{A}{\tau }_{f1}+\tilde{C}{I}_{f1}---\left(20\right)$其中，I_f1为第一级快变时标σ₁下的控制电流，下标f1表示系统处在第一级快变时标下；下面取第二个小参数其中k=min(k₁,k₂)，且两个小参数满足0<ε₁<<ε₂<<1，得到系统的第二级慢变流形表达式为$M_{1,s2}(\theta ,0)\stackrel{··}{\theta }+G_{1,s2}(\theta ,\stackrel{·}{\theta },\mathrm{0,0})=-{\tilde{A}}^{-1}\tilde{B}\stackrel{·}{\theta }+{\tilde{A}}^{-1}\tilde{C}{I}_{s2}---\left(28\right)$其中，I_s2为第二级慢变时标t下的控制电流，下标s2表示系统处在第二级慢变时标下；在边界层上引入第二级快变时标得液压刚柔机械臂系统第二级快变子系统的动力学方程为$\frac{d^2{y}_{f2}}{{\mathrm{d\sigma }}_2^2}=-D_{4,s2}(\theta ,{ϵ}_{2}y)\tilde{K}{y}_{f2}+D_{3,s2}(\theta ,{ϵ}_{2}y){\tilde{A}}^{-1}\tilde{C}{I}_{f2}---\left(31\right)$其中，I_f2为第二级快变时标σ₂下的控制电流，下标f2表示系统处在第二级快变时标下；第二步、根据第一步得到的第二级慢变子系统、第二级快变子系统以及第一级快变子系统动力学方程，设计相应的子系统控制模块：针对第二级慢变子系统，取一阶滑动模态面及二阶滑动模态面设计二阶滑模控制律如式(37)$I_{s2}={\tilde{C}}^{-1}\tilde{A}[M_{1,s}({\stackrel{··}{\theta }}_d+\mathrm{\alpha e}+\beta \stackrel{·}{e}+u)+{\tilde{A}}^{-1}\tilde{B}\stackrel{·}{\theta }+G_{1,s}]---\left(37\right)$针对第二级快变子系统，取二次型性能指标函数及Ricatti方程$A_k^{T}P+{\mathrm{PA}}_k-{\mathrm{PB}}_k{R}^{-1}{B}_k^{T}P+Q=0,$设计最优控制律如式(41)$I_{f2}=-K_f{X}_k={-R}^{-1}{B}_k^T{\mathrm{PX}}_k---\left(41\right)$针对第一级快变子系统，取滑模面s_f=e_f及自适应更新律设计自适应滑模控制律如式(50)$I_{f1}={\tilde{C}}^{-1}[{\stackrel{·}{\tau }}_{f1d}+{\tilde{A}}_0{\tau }_{f1}+\hat{F}+D_{0}E+\mathrm{\xi sgn}\left(s_f\right)]---\left(50\right)$第三步、根据多重时间尺度理论，将第二步得到的各子系统控制器组合得到液压刚柔机械臂组合控制器如式(52)I=I_f1+I_s1=I_f1+I_f2+I_s2 (52)进而完成基于双参数奇异摄动的液压刚柔机械臂的控制方法。