Pseudospectral Method Based Time Optimal Anti-swing Trajectory Planning for Double Pendulum Crane Systems
-
摘要: 在工业生产过程中,桥式吊车系统经常会体现出双摆系统的特性,导致更多欠驱动状态量的出现,增大控制难度.基于此,论文提出了一种针对双摆桥式吊车系统的时间最优轨迹规划方法,可以得到全局时间最优且具有消摆能力的轨迹.具体而言,为方便地构造以时间为代价函数的优化问题,首先对系统运动学模型进行相应的变换;在此基础上,考虑包括两级摆角及台车速度和加速度上限值在内的多种约束,构造出相应的优化问题;然后,利用高斯伪谱法(Gauss-pseudospectral method, GPM)将该带约束的优化问题转化为更易于求解的非线性规划问题,且在转化过程中,可以非常方便地考虑轨迹约束.求解该非线性规划问题,即可得到时间最优的台车轨迹.不同于已有的大多数方法,该方法可获得全局时间最优的结果.最后,通过仿真与实验结果验证了这种时间最优轨迹规划方法具有满意的控制性能.Abstract: In practice, an overhead crane system may behave like a double pendulum, which has more unactuated states and is more difficult to be controlled properly. Motivated by this observation, we present a time optimal trajectory planning scheme for double pendulum crane systems, which can yield a global time-optimal swing-free trajectory. Specifically, to facilitate the optimization problem creation process, we first implement some basic transformations on the system kinematics. Then, various constraints, including upper and lower bounds of the two pendulum angles and upper bounds of the trolley velocity and acceleration, are taken into consideration to set up the optimization problem. After that, the Gauss-pseudospectral method(GPM) is utilized to convert the constrained optimization problem into a nonlinear programming problem, which can be solved more conveniently, while the trajectory constraints are also considered during the transformation. By solving the constructed nonlinear programming problem, a global time-optimal result is obtained, which is different from most existing methods. Finally, numerical simulation and experimental results are given to illustrate the satisfactory performance of the proposed method.
-
在工业生产过程中,为运送负载到所期望的位置,包括桥式吊车、悬臂式吊车、塔式吊车、船用吊车在内的各类吊车系统,有着非常广泛的应用. 为了简化吊车系统的机械结构,往往不直接控制负载,而是通过台车的运动,间接地拖动负载至目标位置.这种结构带来的结果是,吊车系统的控制输入维数小于待控自由度维数.具有该特性的系统即所谓的欠驱动系统[1]. 相比全驱动系统,由于欠驱动特性的存在,欠驱动系统往往更难控制. 而在实际生产中,吊车系统往往是由经验丰富的工人操作. 如果发生误操作,可能引起负载剧烈摆动导致碰撞,甚至发生安全事故. 因此,对吊车系统自动控制方法的研究具有现实意义与广泛应用价值,得到了广大学者的关注.
对于桥式吊车系统,主要的控制目标包括两个方面,即快速且精确的台车定位与负载摆动的抑制与消除. 然而,这两方面通常是互相矛盾的,即过快的台车运动往往导致较大的负载摆动.因此,同时实现这两方面的控制目标具有较高的难度.为了获得较好的控制效果,目前国内外学者已经提出了很多吊车系统自动控制方法.Tuan等提出了基于部分反馈线性化的控制方法[2-3],可以简化桥式吊车系统的控制算法设计. 在文献[4-5]中,Singhose等和Blackburn等利用输入整形的思想对吊车系统进行控制,可有效地抑制负载残余摆动. 为处理不确定的外界干扰,研究人员利用滑模算法控制吊车系统[6-7],可以获得很好的鲁棒性.胡洲等提出了一种非线性信息融合控制方法[8],可以处理控制器输入饱和的问题,实现对吊车系统的高性能控制.文献[9-10]提出了基于能量与无源性的控制策略,可获得较好的效果.除此之外,近年来,包括遗传算法[11]、模糊控制[12]等一些智能控制方法同样在吊车控制领域有着一定的应用.
众所周知,吊车系统的负载摆动由台车的加减速运动引起,在台车运动与负载摆动之间存在着较强的耦合. 基于此,可以为台车规划一条合适的轨迹,当台车按照该轨迹运动时,即可实现对其快速精确定位的目标. 同时,考虑到摆角抑制与消除的要求,在轨迹规划过程中,通过深入分析与合理利用台车运动与负载摆动之间的耦合关系,可以规划出一条具有消摆能力的台车轨迹.这样即可完成台车快速精确定位与负载消摆的双重目标. 基于该思想,研究人员已提出了很多吊车轨迹规划方法[13-17]. 例如,在文献[13]中,Uchiyama 等提出了一种针对悬臂吊车的开环控制策略,为了消除残余摆动,他们为悬臂的水平运动规划了一条S型轨迹.文献[14]则提出了一种基于相平面分析的轨迹规划方法,它可以较好地抑制负载摆动,并消除了残余摆动.
对于吊车系统,虽然上述各种方法在理想情况下可以得到较好的控制效果,但是这些方法均假设吊钩的质量可以忽略且负载可看作质点,并将吊车系统的负载摆动视为单摆系统. 若吊钩的质量较大,无法忽略,或者负载形状很大,不能简单看作质点,在这种情况下,吊车系统的摆动将会呈现出双摆现象,即吊钩绕台车进行一级摆动,同时负载绕吊钩发生二级摆动. 此时,上述的单摆吊车控制方法无法取得令人满意的控制性能,目前仅有少量文献报道了双摆效应的吊车系统控制策略.文献[18-20]通过分析双摆吊车系统的固有频率,将输入整形的方法成功地扩展到双摆吊车系统.孙宁等在考虑系统摆角约束、台车速度约束等一系列约束的前提下,提出了一种基于微分平坦理论的双摆吊车最优轨迹规划方法[21].郭卫平等通过分析吊车系统的能量,提出了一种基于无源性的双摆吊车控制策略[22].尽管上述控制策略可以实现对双摆吊车系统的控制,但它们难以保证吊车系统运行效率的最大化.
为提高双摆吊车系统的运送效率,本文提出了一种全局时间最优轨迹规划方法,可以同时实现台车精确定位以及负载快速消摆的控制目标.与现有方法不同的是,本文所提方法可得到全局时间最优的台车轨迹.具体而言,首先对吊车系统的运动学模型进行转化,得到一种加速度驱动的系统模型. 接着,利用该模型,考虑吊车系统工作过程中的各种物理及安全约束,构造出以运送时间为优化目标的待优化函数. 为了便于求解该优化问题,利用高斯伪谱法(Gauss-pseudospectral method),将所得的优化问题以及相应约束,在勒让德-高斯(Legendre-Gauss,LG)点处进行离散化与近似化. 这样,原优化问题即转化为一种具有代数约束的非线性规划问题.通过求解该非线性规划问题,便可得到所需的时间最优台车轨迹.在整个规划过程中,可方便地考虑包括两级摆角约束以及台车速度和加速度约束在内的各种实际物理约束,因此得到的轨迹光滑合理,便于实际跟踪. 最后,通过数值仿真以及实验结果验证了所提方法的有效性.
本文剩余部分组织如下: 第1节给出了双摆吊车运动学模型,以及所需考虑的各种轨迹约束; 第2节具体描述了所设计轨迹规划的过程,包括模型变换与基于高斯伪谱法的优化问题转化与求解;第3节给出了所提方法的仿真与实验结果,并与已有方法进行分析比较;第4节对本文进行了总结与展望.
1. 问题描述
在本文中,考虑如图 1所示的双摆桥式吊车系统,其运动学模型如下[21]:
\begin{align}&(m_1+m_2) l_1 \cos{\theta_1} \ddot{x}+ (m_1+m_2) l_1^2\ddot{\theta}_1 +\notag \\&\quad m_2 l_1 l_2 \cos{(\theta_1-\theta_2) } \ddot{\theta}_2+m_2 l_1 l_2 \sin{(\theta_1-\theta_2) }\dot{\theta}_2^2+ \notag \\&\quad (m_1+m_2) g l_1 \sin{\theta_1}=0 \end{align}
(1) \begin{align}& m_2 l_2 \cos\theta_2 \ddot{x} + m_2 l_1 l_2 \cos(\theta_1-\theta_2) \ddot{\theta}_1 + m_2 l_2^2 \ddot{\theta}_2- \notag \\&\quad m_2 l_1 l_2 \dot{\theta}_1^2 \sin(\theta_1-\theta_2) + m_2g l_2 \sin\theta_2=0\end{align}
(2) 其中, $M,~m_1,~m_2$ 分别代表台车质量、吊钩质量以及负载质量; $l_1$ 表示吊绳的长度, $l_2$ 为等效绳长,即负载质心与吊钩质心之间的距离; $x(t)$ 为台车位移, $\theta_1(t),~\theta_2(t)$ 表示一级与二级摆角(吊钩摆角与负载绕吊钩的摆角); $g$ 代表重力加速度常数.
对该双摆桥式吊车系统,拟规划一条时间最优轨迹,以实现快速与精确的台车定位,同时完成两级摆动的抑制与消除的控制目标. 考虑到该系统的欠驱动特性,无法对两级摆动进行直接控制,只能通过对台车的运动进行合理的规划,以此来间接地控制负载摆动. 为此,对式(1) 和(2) 分别除以 $(m_1+m_2) l_1$ 和 $m_2 l_2$ ,化简可得:
\begin{align}&\cos\theta_1 \ddot{x}+l_1 \ddot{\theta}_1 + \frac{m_2 l_2}{m_1+m_2} \cos(\theta_1-\theta_2) \ddot{\theta}_2+ \notag \\&\quad \frac{m_2 l_2}{m_1+m_2}\sin(\theta_1-\theta_2) \dot{\theta}_2^2+g\sin\theta_1=0 \end{align}
(3) \begin{align}&\cos\theta_2\ddot{x}+l_1 \cos(\theta_1-\theta_2) \ddot{\theta}_1+l_2 \ddot{\theta}_2- \notag \\&\quad l_1 \dot{\theta}_1^2\sin(\theta_1-\theta_2) +g\sin\theta_2=0\end{align}
(4) 式(3) 和(4) 描述了台车位移 $x(t)$ 与系统两级摆角 $\theta_1(t),~\theta_2(t)$ 之间的耦合关系,即台车的运动对负载摆动的影响. 通过深入分析该耦合关系,规划一条具有消摆能力的台车轨迹,是本方法的基础.
为完成时间最优轨迹规划,考虑到吊车系统在实际工作时的目标、物理约束及安全性,本文将系统地考虑如下几个方面的轨迹约束[21]:
1) 为实现快速而精确的台车定位,台车从初始位置 $x_0$ 开始运动, $T$ 时刻达到目标位置 $x_f$ ,且开始时刻与结束时刻的台车速度、加速度均为0,即
\begin{align}&x(0) =\dot{x}(0) =\ddot{x}(0) =0 \notag \\&x(T)=x_f,~\dot{x}(T)=\ddot{x}(T)=0 \notag\end{align}
其中, $T$ 表示运送过程所需要时间; 对于初始位置,不失一般性,本文选取 $x_0=0$ .
2) 考虑到系统输入饱和问题,不同于文献[23-24]中的方法,这里等价于考虑在运送过程中,台车的速度、加速度均应保持在适当的范围内,即
\begin{align}|\dot{x}(t)| \leq v_{\max},\quad |\ddot{x}(t)| \leq a_{\max}\notag\end{align}
其中, $v_{\max},~a_{\max}$ 分别表示允许的台车最大速度和加速度.
3) 为保证在运送结束时可直接对负载进行下一步的处理,台车达到目标位置后应无残余摆动且角速度也均为0,即
\begin{align}&\theta_1(0) =\dot{\theta}_1(0) =0,\quad \theta_1(T)=\dot{\theta}_1(T)=0 \notag \\&\theta_2(0) =\dot{\theta}_2(0) =0,\quad\theta_2(T)=\dot{\theta}_2(T)=0\notag\end{align}
4) 为避免由于负载的剧烈摆动引起的碰撞,在运送过程中,两级摆动的摆角及角速度都应保持在允许的范围内,即
\begin{align}&|\theta_1(t)| \leq \theta_{1\max},\quad |\theta_2(t)| \leq \theta_{2\max} \end{align}
(5) \begin{align}&|\dot{\theta}_1(t)| \leq \omega_{1\max},\quad |\dot{\theta}_2(t)|\leq \omega_{2\max}\end{align}
(6) 其中, $\theta_{1\max},~\theta_{2\max}$ 分别表示一级与二级摆角允许的最大角度; $\omega_{1\max},~\omega_{2\max}$ 代表允许的一级与二级最大角速度.
综上,可以构造如下的优化问题:
\begin{align}{\rm{min}}\qquad &T \\{\rm{s.t.}}\qquad &x(0) =\dot{x}(0) =\ddot{x}(0) =0 \notag\\&x(T)=x_f,\quad \dot{x}(T)=\ddot{x}(T)=0 \notag\\&|\dot{x}(t)| \leq v_{\max},\quad |\ddot{x}(t)| \leq a_{\max} \notag \\&\theta_1(0) =\dot{\theta}_1(0) =0,\quad \theta_1(T)=\dot{\theta}_1(T)=0 \notag \\&\theta_2(0) =\dot{\theta}_2(0) =0,\quad \theta_2(T)=\dot{\theta}_2(T)=0 \notag \\&|\theta_1(t)| \leq \theta_{1\max},\quad |\theta_2(t)| \leq \theta_{2\max} \notag \\&|\dot{\theta}_1(t)| \leq \omega_{1\max},\quad|\dot{\theta}_2(t)| \leq \omega_{2\max} \notag\end{align}
(7) 接下来,将通过伪谱法求解该优化问题,并为台车规划出一条时间最优轨迹.
2. 时间最优轨迹规划方法
本节将提出一种基于伪谱法的时间最优轨迹规划策略,并通过求解式(7) 得到台车轨迹. 具体而言,首先将吊车系统运动学模型转化成一种加速度驱动模型,其中台车加速度可看作系统输入量; 随后,基于该加速度驱动模型,可将原优化问题改写成一种新的形式; 接着利用高斯伪谱法,将带约束优化问题转化为一系列非线性规划问题;最后通过求解该非线性优化问题,即可得到时间最优的台车轨迹.
2.1 系统模型及优化问题的转化
为方便后续的轨迹规划,这里首先对双摆吊车系统模型以及优化问题(7) 进行转化.为此定义系统全状态向量 ${\pmb \zeta}(t)$ 如下:
\begin{align}{\pmb \zeta}=[\begin{array}{cccccc} x & \dot{x} & \theta_1 &\dot{\theta}_1 & \theta_2 & \dot{\theta}_2\end{array}]^{\rm{T}}\end{align}
(8) 根据系统运动学模型(3) 和(4) ,可将台车的加速度作为系统的输入. 此时,运动学模型可转化为如下形式:
\begin{align}\dot{{\pmb \zeta}}={\pmb f}({\pmb \zeta})+{\pmb h}({\pmb \zeta})u\end{align}
(9) 其中, $u(t)$ 即为台车加速度 $\ddot{x}(t)$ ; ${\pmb f}({\pmb \zeta}),{\pmb h}({\pmb \zeta})$ 表示关于 ${\pmb \zeta}(t)$ 的辅助函数,具体形式如下:
\begin{align}&{\pmb f}({\pmb \zeta})=\left[\begin{array}{cccccc} \dot{x} & 0 &\dot{\theta}_1 & A & \dot{\theta}_2 & B \end{array} \right]^{\rm{T}} \end{align}
(10) \begin{align}&{\pmb h}({\pmb \zeta})=\left[\begin{array}{cccccc} 0 & 1 & 0 & C &0 & D\end{array} \right]^{\rm{T}}\end{align}
(11) 其中,为方便描述,定义了如下的辅助变量 $A,~B,~C,~D$ :
\begin{align}A=~&-\frac{m_2 C_{1-2}}{l_1 (m_1+m_2) -m_2 l_1 C_{1-2}^2} \bigg[l_1 S_{1-2} \dot{\theta}_1^2 +\notag \\&\frac{m_2 l_2}{m_1+m_2}S_{1-2}C_{1-2}\dot{\theta}_2^2-g(S_2-S_1 C_{1-2})\bigg]- \notag \\&\frac{1}{l_1}gS_1-\frac{m_2 l_2}{l_1(m_1+m_2) }S_{1-2}\dot{\theta}_2^2 \notag \\B=~& \frac{m_1+m_2}{l_2 (m_1+m_2) -m_2 l_2 C_{1-2}^2}\bigg[\frac{m_2 l_2}{m_1+m_2}S_{1-2}× \notag \\&C_{1-2}\dot{\theta}_2^2+l_1 S_{1-2}\dot{\theta}_1^2-g(S_2-S_1C_{1-2})\bigg] \notag\\C=~& \frac{m_2 C_{1-2}}{l_1 (m_1+m_2) -m_2 l_1 C_{1-2}^2} (C_2-C_1 C_{1-2})- \frac{1}{l_1} C_1 \notag \\D=~&-\frac{m_1+m_2}{l_2 (m_1+m_2) -m_2 l_2 C_{1-2}^2}(C_2-C_1C_{1-2}) \notag\end{align}
式中,使用了如下的简化形式:
\begin{align}&S_1=\sin\theta_1,\quad S_2=\sin\theta_2,\quad C_1=\cos\theta_1,\quad C_2=\cos\theta_2 \notag \\&S_{1-2}=\sin(\theta_1-\theta_2) ,\quad C_{1-2}=\cos(\theta_1-\theta_2) \notag\end{align}
利用所得的加速度驱动系统模型(9) ,原优化问题(7) 可以转化为
$\text{min}\qquad \text{ }T$
(12a) \begin{align}{\rm{s.t.}} \qquad &\dot{\pmb \zeta}={\pmb f}({\pmb \zeta})+{\pmb h}({\pmb \zeta})u \end{align}
(12b) \begin{align}&{\pmb \zeta}(0) =[\begin{array}{cccccc} 0 & 0 & 0 & 0 & 0 &0\end{array}]^{\rm{T}} \end{align}
(12c) \begin{align}&{\pmb \zeta}(T)=[\begin{array}{cccccc} x_f & 0 & 0 & 0 & 0 & 0\end{array}]^{\rm{T}} \end{align}
(12d) \begin{align}&|\dot{x}(t)| \leq v_{\max},\quad |u(t)| \leq a_{\max} \end{align}
(12e) \begin{align} &|\theta_1(t)| \leq \theta_{1\max},\quad |\theta_2(t)| \leq \theta_{2\max}\end{align}
(12f) \begin{align}&|\dot{\theta}_1(t)| \leq \omega_{1\max},\quad|\dot{\theta}_2(t)| \leq \omega_{2\max} \end{align}
(12g) 求解该优化问题,即可得到完成控制目标所需的最优时间 $T^*$ ,同时得到对应的最优台车轨迹.
2.2 基于高斯伪谱法的轨迹规划
为得到时间最优的台车轨迹,关键是如何求解所构造的带约束优化问题.本文通过高斯伪谱法对优化问题进行处理,可方便地得到时间最优解以及最优轨迹. 与现有大部分方法不同的是,这种方法可更为直接地对原系统进行分析和处理,并可以获得全局时间最优解.
本文方法的主要思想如图 2中的流程图所示,具体描述如下:首先利用拉格朗日插值方法,选择LG点处的离散系统状态轨迹以及输入轨迹,通过离散轨迹与拉格朗日插值多项式,表示相应的近似轨迹模型. 接着,通过对近似后的轨迹模型进行求导,可将系统状态的导数用拉格朗日多项式导数表示. 随后,利用离散的轨迹模型及其导数,可将原系统模型(微分方程约束)转化为一系列多项式方程. 利用高斯积分,优化问题中的边界条件同样可以表示成多项式方程的形式. 最后,时间最优轨迹规划问题即转化为一种具有代数约束的非线性规划问题,通过对其求解即可得到最优时间及最优轨迹.
为适应高斯伪谱法的要求,首先需要利用坐标变换,将轨迹对应的时间区间由 $t\in[0,T]$ 转化到区间 $\tau \in[-1,1]$ 上,即
\begin{align}\tau=\frac{2t}{T}-1\end{align}
(13) 随后选取 $K$ 个节点 $\{\tau_1,\tau_2,\cdots,\tau_K\}\in (-1,1) $ 组成点列,即LG点.LG点的选取可以通过求解 $K$ 阶的勒让德多项式的零点获得. 同时,把 $\tau_0=-1$ 加到点列的首位.待规划的系统状态量及输入量可离散表示成如下的形式:
\begin{align}&{\pmb \zeta}(\tau_0) ,{\pmb \zeta}(\tau_1) ,{\pmb\zeta}(\tau_2) ,\cdots,{\pmb \zeta}(\tau_K) \notag \\&u(\tau_0) ,u(\tau_1) ,u(\tau_2) ,\dots,u(\tau_K)\notag\end{align}
利用该 $K+1$ 个节点,可以构造出 $K+1$ 个拉格朗日插值多项式,具体形式如下:
\begin{align}\mathcal{L}_i(\tau)=\prod^K_{j=0,j\neq i}{\frac{\tau-\tau_j}{\tau_i-\tau_j}}\end{align}
(14) 其中, $\tau \in [-1,1]$ ; $\mathcal{L}_i(\tau)$ 代表第 $i$ 个拉格朗日插值多项式函数; $i\in {0,1,\cdots,K}$ . 利用式(14) 以及LG点处的系统状态量和输入量的值,系统的状态量轨迹以及输入量轨迹可通过下面的方式近似表出:
\begin{align}{\pmb \zeta}(\tau) \approx \sum_{i=0}^{K}{{\pmb\zeta}(\tau_i)\mathcal{L}_i(\tau)},~~~u(\tau) \approx\sum_{i=0}^{K}{u(\tau_i)\mathcal{L}_i(\tau)}\end{align}
(15) 其中, ${\pmb\zeta}(\tau_i),~u(\tau_i)$ 分别表示 $\tau=\tau_i$ 处的系统状态量以及输入量.对式(15) 求导,并利用式(14) 中插值函数的具体形式,计算并化简,可得状态量轨迹的导数如下:
\begin{align}\dot{\pmb \zeta}(\tau_k)=\sum^K_{i=0}{\mathcal{D}_{ki}(\tau_k){\pmb\zeta}(\tau_i)}\end{align}
(16) 其中, $\dot{\pmb\zeta}(\tau_k)$ 表示 $\tau=\tau_k$ 处的状态轨迹导数值; $\mathcal{D}_{ki}(\tau_k)$ 代表 $\tau=\tau_k$ 处, $\mathcal{L}_i$ 的导数值,具体形式如下:
\begin{align}\mathcal{D}_{ki}(\tau_k)=\dot{\mathcal{L}}_i(\tau_k)=\sum^K_{l=0,l\neq i}{\frac{\prod\limits_{j=0,j\neq i,l}^K(\tau_k-\tau_j)}{\prod\limits_{j=0,j\neq i}^K(\tau_i-\tau_j)}}\end{align}
(17) 利用式(16) 和(17) 以及LG点处的轨迹值,可对优化问题(12) 中的微分方程约束(12b)进行离散化与近似化处理,具体结果如下:
\begin{align}\sum^K_{i=0}{\mathcal{D}_{ki}(\tau_k){\pmb\zeta}(\tau_i)}-&\frac{2}{T}[{\pmb f}({\pmb \zeta}(\tau_k))+{\pmb h}({\pmb \zeta}(\tau_k))u(\tau_k)]=0\end{align}
(18) 其中, $k\in\{0,1,2,\cdots,K\}$ . 式(18) 可以看作代数约束.接下来,优化问题中的边界条件约束也需要转化成代数约束的形式,其中式(12c)可直接改写如下:
\begin{align}{\pmb \zeta}(\tau_0) =[\begin{array}{cccccc} 0 & 0 & 0 & 0 & 0 & 0\end{array}]^{\rm{T}}\end{align}
(19) 为表示运送过程结束时刻的边界约束,定义 $\tau_{K+1}=1$ ,由式(13) 可知, $\tau_{K+1}=1$ 即对应运送结束时刻 $t=T$ . 利用高斯积分,式(12d)可表示为
\begin{align}{\pmb \zeta}(\tau_{K+1})=&{\pmb\zeta}(\tau_0) +\frac{2}{T}\sum^K_{k=1}w_k[{\pmb f}({\pmb\zeta}(\tau_k))+{\pmb h}({\pmb \zeta}(\tau_k))u(\tau_k)]\end{align}
(20) 其中, ${\pmb \zeta}(\tau_0) $ 即上述的系统初始状态向量; $w_k$ 表示第 $k$ 个勒让德权值(Legendre weight),具体值可在求解LG点时一并求得.
综上,优化问题中所有的约束都可以通过代数约束的形式表出,基于此,原优化问题可以转化成一种具有代数约束的非线性规划问题,具体如下所示:
\begin{align}&{\rm{min}}~~~~~T \end{align}
(21) \begin{align}&{\rm{s.t.}} \notag\\&\sum^K_{i=0}{\mathcal{D}_{ki}(\tau_k){\pmb \zeta}(\tau_i)}-\frac{2}{T}[{\pmb f}({\pmb \zeta}(\tau_k))+ \notag \\&{\pmb h}({\pmb \zeta}(\tau_k))u(\tau_k)]=0,~~~~~~k\in\{0,1,2,\cdots,K\} \notag \\&{\pmb \zeta}(\tau_0) =[\begin{array}{cccccc} 0 & 0 & 0 & 0 & 0 & 0\end{array}]^{\rm{T}} \notag\\&{\pmb \zeta}(\tau_0) +\frac{2}{T}\sum^K_{k=1}{w_k[{\pmb f}({\pmb \zeta}(\tau_k))+{\pmb h}({\pmb \zeta}(\tau_k))u(\tau_k)]}=\notag \\&[\begin{array}{cccccc} x_f & 0 & 0 & 0 & 0 & 0\end{array}]^{\rm{T}} \notag\\&{\pmb \zeta}(\tau)- {\pmb \gamma} \leq 0,~ -{\pmb \zeta}(\tau)-{\pmb \gamma} \leq 0 \notag\\& u(\tau)-a_{\max} \leq 0,~ -u(\tau)-a_{\max} \leq 0 \notag\end{align}
其中,向量 ${\pmb \gamma}$ 的具体形式如下:
\begin{align}{\pmb \gamma}=[\begin{array}{cccccc} \infty & v_{\max} &\theta_{1\max} & \omega_{1\max} & \theta_{2\max} & \omega_{2\max}\end{array}]^{\rm{T}} \notag\end{align}
对于上述带约束非线性规划问题,本文选择连续二次型规划方法(Sequential quadratic programming,SQP)求解,可得到如下的时间最优状态向量序列:
\begin{align}{\pmb \zeta}(\tau_0) ,{\pmb \zeta}(\tau_1) ,{\pmb\zeta}(\tau_2) ,\cdots,{\pmb \zeta}(\tau_K),{\pmb\zeta}(\tau_{K+1})\end{align}
(22) 上式即时间离散的最优状态向量序列.取每个向量的前两项(台车位移与台车速度),并进行插值,即可得到相应的时间最优的台车位移与速度轨迹.
3. 仿真与实验
在本节中,为验证所设计的轨迹规划方法的有效性,给出了相应的仿真与实验结果. 具体而言,首先,利用Matlab及其相关工具箱,按照所设计的轨迹规划方法,离线计算出时间最优的台车轨迹; 随后,利用Matlab/Simulink环境对其进行仿真测试; 最后,将构造的最优轨迹在双摆吊车实验平台上进行实验验证,并与现有的方法进行对比.所得结果均验证了本文所设计的轨迹规划方法具有良好的性能.
3.1 仿真结果
为实现本文所提出的基于伪谱法的轨迹规划策略,这里使用GPOPS软件工具箱[25]以及SNOPT工具箱[26]离线求解优化问题并得到相应的时间最优轨迹.相应的系统参数按照实验平台中的实际参数进行选取,具体值如下:
\begin{align}&M=6. 5~{\rm{kg}},\quad m_1=2. 003~{\rm{kg}},\quad m_2=0. 559~{\rm{kg}} \notag\\&g=9. 8~{\rm{m/s^2}},\quad l_1=0. 53~{\rm{m}},\quad l_2=0.4~{\rm{m}} \notag\end{align}
台车的目标位置选取为 $x_f=0. 6~{\rm{m}}$ ,同时轨迹约束选择如下:
\begin{align}&\theta_{1\max}=\theta_{2\max}=2~{\rm{deg}},\quad v_{\max}=0. 3~{\rm{m/s}} \notag\\&\omega_{1\max}=\omega_{2\max}=5~{\rm{deg/s}},\quad a_{\max}=015~{\rm{m/s^2}} \notag\end{align}
仿真结果如图 3和4所示. 从图 3中可以看出,当台车按照所规划轨迹进行运动时,完成给定运送任务仅需时间 $T^*=4.2776~{\rm{s}}$ ,且台车可以精确达到目标位置,整个过程中,台车速度未超过给定限制 $v_{\max}=0. 3~{\rm{m/s}}$ . 同时,由图 3可知,一级摆动和二级摆动的最大摆角均未超过给定值 $2~\rm{deg}$ ,且运送结束时不存在残余摆动. 这也就保证了负载运送过程中的安全性.同样地,两级摆动对应的角速度也均保持在给定范围内. 综上所述,仿真结果验证了所提最优轨迹规划方法的高效性与安全性.
图 3 本文方法仿真结果(台车位置以及速度) (实线:仿真结果; 虚线: 目标位置 $x_f=0. 6~{\rm{m}}$ ; 点画线:台车速度约束 $v_{\max}=0. 3~{\rm{m/s}})$Fig. 3 Simulation results (trolley position and velocity)(Solid line: simulation results; Dashed line: target position $x_f=0. 6~{\rm{m}}$ ; Dotted-dashed line: trolley velocity constraint $v_{\max}=0. 3~{\rm{m/s}}$ )图 4 本文方法仿真结果(两级摆动对应的摆角及角速度) (实线: 仿真结果;虚线: 摆角约束 $\theta_{1\max}=\theta_{2\max}=2~{\rm{deg}}$ ;点画线: 角速度约束 $\omega_{1\max}=\omega_{2\max}=5~{\rm{deg/s}})$Fig. 4 Simulation results (first and second order swing angles and angular velocities) (Solid line: simulation results;Dashed line: swing angle constraint $\theta_{1\max}=\theta_{2\max}=2~{\rm{deg}}$ ; Dotted-dashed line:angular velocity constraint $\omega_{1\max}=\omega_{2\max}=5~{\rm{deg/s})}$3.2 实验结果
为进一步验证所提方法的性能,这里将所规划轨迹在双摆实验平台上进行验证.本方法轨迹规划与跟踪控制的整体流程图如图 5所示. 实验过程中,为使台车沿所规划的轨迹运动,本文选择工业中有广泛应用的比例-微分(Proportional-derivative,PD)控制器作为跟踪控制器,以实现台车对目标轨迹的跟踪.\hfill实验平台的参数与仿真中相同,
因此通过离线计算所得的最优轨迹即仿真中给出的轨迹. 同时,作为对比,这里给出了文献[21]的最优轨迹规划方法,以及线性二次型调节器(Linear quadratic regulator,LQR)方法的实验结果. 其中,文献[21]中轨迹规划方法的约束选取与本文所提方法的约束一致;而对于LQR方法,其控制器表达式如下:
\begin{align}F=-k_1(x-x_f)-k_2\dot{x}-k_3\theta_1-k_4\dot{\theta}_1-k_5\theta_2-k_6\dot{\theta}_2\notag\end{align}
代价函数选取如下:
\begin{align}J=\int^{\infty}_{0} ({\pmb X}^{\rm{T}}Q{\pmb X}+RF^2) {\rm d}t\notag\end{align}
其中, ${\pmb X}$ 的定义如下:
\begin{align}{\pmb X}=[\begin{array}{cccccc}e(t) & \dot{x}(t) & \theta_1(t) & \dot{\theta}_1(t) & \theta_2(t)& \dot{\theta}_2(t)\end{array}]^{\rm{T}} \notag\end{align}
$e(t)$ 代表台车定位误差, $e(t)=x(t)-x_f$ ; 矩阵 $Q,~R$ 的选择如下:{align}Q={\rm diag}\{200,1,200,1,200,1\},~{ R}=0. 05 \notag{align}计算得到控制器增益如下:
\begin{align}&k_1=63. 2456,~k_2=50. 7765,~k_3=-129. 3086 \notag\\&k_4=-6. 9634,~k_5=19. 9137,~k_6=-6. 7856 \notag\end{align}
本文方法、LQR方法、文献[21]中方法的实验结果如图 6~8所示.从图 6中可以看出,当台车跟踪所规划最优轨迹运行时,完成运送过程仅需 $4. 095~\rm{s}$ ,且整个过程中两级摆动均保持在给定约束 $2~\rm{deg}$ 以内,运送完成时基本无残余摆动. 而对于文献[21]所提方法,完成运送过程需要 $5. 445~\rm{s}$ ; 对于LQR方法,完成运送过程需要 $7. 425~\rm{s}$ . 同时,LQR方法导致的一级摆动最大摆角达到 $6. 5~\rm{deg}$ ,二级摆动最大摆角达到 $11. 5~\rm{deg}$ ,远大于本文所提方法的摆角.综上可知,本文所设计的方法可以实现台车的精确定位以及系统两级摆动的快速消除,获得良好的控制性能. 从图 3、图 4和图 6可以看出,由于摩擦等多种干扰的影响,实验与仿真结果存在着一定程度上的区别.同时,对比图 4与图 6中的摆角结果,可知摆角曲线的形状以及趋势大致相同,且对摆角的限幅及消摆的目标也基本完成. 因此我们可得如下结论,即本文中的仿真与实验结果基本一致,均验证了所提方法可以实现较好的控制效果. 需要指出的是,由于外界干扰的影响,一般来说,实验结果与仿真结果均存在一定程度上的区别[13-14].
图 6 本文方法实验结果(台车位置、 一级摆角、二级摆角) (实线: 实验结果; 虚线: 待跟踪最优轨迹; 点画线:摆角约束 $\theta_{1\max}=\theta_{2\max}=2~{\rm{deg})}$Fig. 6 Experimental results of proposed method (trolley position,first and second order swing angles) (Solid line:experimental results; Dashed line: planned trajectory; Dotted-dashed line: swing angle constraint $\theta_{1\max}=\theta_{2\max}=2~{\rm{deg}}$ )注1. 图 6中,本文所提方法的实验结果,在台车停止运动后存在着一定的残摆,但其幅度很小,仅为 $0.2~\rm{deg}$ .这与图 4中本方法仿真结果有一定的差异. 需要指出的是,残摆的出现并非由于本问方法存在缺陷,而是因为在实验过程中,存在着摩擦等多种外界干扰的影响. 同时,考虑到残摆的幅度很小,在实际吊车工作过程中,这种幅度的残摆基本可以忽略. 需要注意的是,虽然本文方法以及对比方法是在相同的实验平台上进行测试,但每次测试时,轨道摩擦等外界干扰均会略有差别. 只会在几次实验中,得到基本无残摆的结果. 除此之外,一般来说,实验结果存在小幅度的残摆属于正常的情况[13, 18].
图 8 文献[21]方法实验结果(台车位置、一级摆角、 二级摆角) (实线: 实验结果; 虚线: 待跟踪轨迹; 点画线:摆角约束 $\theta_{1\max}=\theta_{2\max}=2~{\rm{deg}}$ )Fig. 8 Experimental results of the method in [21] (trolley position,first and second order swing angles) (Solid line:experimental results; Dashed line: planned trajectory; Dotted-dashed line: swing angle constraint$\theta_{1\max}=\theta_{2\max}=2~{\rm{deg})4. 结论
本文针对具有双摆效应的桥式吊车,提出了一种基于伪谱法的时间最优轨迹规划策略. 具体而言,首先,将吊车系统的运动学模型转化为一种加速度驱动模型,并基于此模型,考虑各种约束,构造出带约束的优化问题; 随后,利用高斯伪谱法对所得优化问题进行处理,将其转化为更方便求解的非线性规划问题. 在此基础上,即可得到时间最优台车轨迹.本文提出的这种轨迹规划方法除了考虑消摆的目标之外,还可以非常方便地处理摆角约束、角速度约束、台车速度约束、加速度约束等实际物理约束.与现有方法不同的是,本文提出的方法可获得全局时间最优解,极大地提高了吊车系统的工作效率. 最后,通过仿真与实验,验证了本文所提方法的有效性. 在未来的工作中,将考虑在轨迹规划过程中加入对角加速度的约束,以获得更好的控制效果,并考虑为具有双摆效应的吊车系统设计闭环控制算法,以提高吊车系统在运行过程中对外界干扰的鲁棒性.
-
图 3 本文方法仿真结果(台车位置以及速度) (实线:仿真结果; 虚线: 目标位置 $x_f=0. 6~{\rm{m}}$ ; 点画线:台车速度约束 $v_{\max}=0. 3~{\rm{m/s}})$
Fig. 3 Simulation results (trolley position and velocity)(Solid line: simulation results; Dashed line: target position $x_f=0. 6~{\rm{m}}$ ; Dotted-dashed line: trolley velocity constraint $v_{\max}=0. 3~{\rm{m/s}}$ )
图 4 本文方法仿真结果(两级摆动对应的摆角及角速度) (实线: 仿真结果;虚线: 摆角约束 $\theta_{1\max}=\theta_{2\max}=2~{\rm{deg}}$ ;点画线: 角速度约束 $\omega_{1\max}=\omega_{2\max}=5~{\rm{deg/s}})$
Fig. 4 Simulation results (first and second order swing angles and angular velocities) (Solid line: simulation results;Dashed line: swing angle constraint $\theta_{1\max}=\theta_{2\max}=2~{\rm{deg}}$ ; Dotted-dashed line:angular velocity constraint $\omega_{1\max}=\omega_{2\max}=5~{\rm{deg/s})}$
图 6 本文方法实验结果(台车位置、 一级摆角、二级摆角) (实线: 实验结果; 虚线: 待跟踪最优轨迹; 点画线:摆角约束 $\theta_{1\max}=\theta_{2\max}=2~{\rm{deg})}$
Fig. 6 Experimental results of proposed method (trolley position,first and second order swing angles) (Solid line:experimental results; Dashed line: planned trajectory; Dotted-dashed line: swing angle constraint $\theta_{1\max}=\theta_{2\max}=2~{\rm{deg}}$ )
图 8 文献[21]方法实验结果(台车位置、一级摆角、 二级摆角) (实线: 实验结果; 虚线: 待跟踪轨迹; 点画线:摆角约束 $\theta_{1\max}=\theta_{2\max}=2~{\rm{deg}}$ )
Fig. 8 Experimental results of the method in [21] (trolley position,first and second order swing angles) (Solid line:experimental results; Dashed line: planned trajectory; Dotted-dashed line: swing angle constraint$\theta_{1\max}=\theta_{2\max}=2~{\rm{deg})
-
[1] Liu Y, Yu H N. A survey of underactuated mechanical systems. IET Control Theory & Applications, 2013, 7(7):921-935 https://www.researchgate.net/profile/Hongnian_Yu/publication/260586292_A_survey_of_underactuated_mechanical_systems/links/02e7e535974fe51e60000000.pdf [2] Tuan L A, Lee S G, Dang V H, Moon S, Kim B. Partial feedback linearization control of a three-dimensional overhead crane. International Journal of Control, Automation and Systems, 2013, 11(4):718-727 doi: 10.1007/s12555-012-9305-z [3] Tuan L A, Kim G H, Kim M Y, Lee S G. Partial feedback linearization control of overhead cranes with varying cable lengths. International Journal of Precision Engineering and Manufacturing, 2012, 13(4):501-507 doi: 10.1007/s12541-012-0065-8 [4] Singhose W, Kim D, Kenison M. Input shaping control of double-pendulum bridge crane oscillations. ASME Journal of Dynamic Systems, Measurement, and Control, 2008, 130(3):034504 doi: 10.1115/1.2907363 [5] Blackburn D, Singhose W, Kitchen J, Patrangenaru V, Lawrence J, Kamoi T, Taura A. Command shaping for nonlinear crane dynamics. Journal of Vibration and Control, 2010, 16(4):477-501 doi: 10.1177/1077546309106142 [6] 王伟, 易建强, 赵冬斌, 刘殿通. 桥式吊车系统的分级滑模控制方法. 自动化学报, 2004, 30(5):784-788 http://www.aas.net.cn/CN/abstract/abstract16252.shtmlWang Wei, Yi Jian-Qiang, Zhao Dong-Bin, Liu Dian-Tong. Hierarchical sliding-mode control method for overhead cranes. Acta Automatica Sinica, 2004, 30(5):784-788 http://www.aas.net.cn/CN/abstract/abstract16252.shtml [7] Xi Z, Hesketh T. Discrete time integral sliding mode control for overhead crane with uncertainties. IET Control Theory & Applications, 2010, 4(10):2071-2081 http://ieeexplore.ieee.org/xpl/articleDetails.jsp?reload=true&arnumber=5611727&contentType=Journals+%26+Magazines [8] 胡洲, 王志胜, 甄子洋. 带输入饱和的欠驱动吊车非线性信息融合控制. 自动化学报, 2014, 40(7):1522-1527 http://www.aas.net.cn/CN/abstract/abstract18422.shtmlHu Zhou, Wang Zhi-Sheng, Zhen Zi-Yang. Nonlinear information fusion control for underactuated cranes with input saturation. Acta Automatica Sinica, 2014, 40(7):1522-1527 http://www.aas.net.cn/CN/abstract/abstract18422.shtml [9] Sun N, Fang Y C, Zhang X B. Energy coupling output feedback control of 4-DOF underactuated cranes with saturated inputs. Automatica, 2013, 49(5):1318-1325 doi: 10.1016/j.automatica.2013.01.039 [10] Sun N, Fang Y C. New energy analytical results for the regulation of underactuated overhead cranes:an end-effector motion-based approach. IEEE Transactions on Industrial Electronics, 2012, 59(12):4723-4734 doi: 10.1109/TIE.2012.2183837 [11] Nakazono K, Ohnishi K, Kinjo H, Yamamoto T. Load swing suppression for rotary crane system using direct gradient descent controller optimized by genetic algorithm. Transactions of the Institute of Systems, Control and Information Engineers, 2011, 22(8):303-310 [12] Zhao Y, Gao H J. Fuzzy-model-based control of an overhead crane with input delay and actuator saturation. IEEE Transactions on Fuzzy Systems, 2012, 20(1):181-186 doi: 10.1109/TFUZZ.2011.2164083 [13] Uchiyama N, Ouyang H M, Sano S. Simple rotary crane dynamics modeling and open-loop control for residual load sway suppression by only horizontal boom motion. Mechatronics, 2013, 23(8):1223-1236 doi: 10.1016/j.mechatronics.2013.09.001 [14] Sun N, Fang Y, Zhang X, Yuan Y. Transportation task-oriented trajectory planning for underactuated overhead cranes using geometric analysis. IET Control Theory & Applications, 2012, 6(10):1410-1423 http://ieeexplore.ieee.org/iel5/4079545/6257070/06257077.pdf [15] 孙宁, 方勇纯, 王鹏程, 张雪波. 欠驱动三维桥式吊车系统自适应跟踪控制器设计. 自动化学报, 2010, 36(9):1287-1294 doi: 10.3724/SP.J.1004.2010.01287Sun Ning, Fang Yong-Chun, Wang Peng-Cheng, Zhang Xue-Bo. Adaptive trajectory tracking control of underactuated 3-dimensional overhead crane systems. Acta Automatica Sinica, 2010, 36(9):1287-1294 doi: 10.3724/SP.J.1004.2010.01287 [16] Sun N, Fang Y C, Zhang Y D, Ma B J. A novel kinematic coupling-based trajectory planning method for overhead cranes. IEEE/ASME Transactions on Mechatronics, 2012, 17(1):166-173 doi: 10.1109/TMECH.2010.2103085 [17] Wang P C, Fang Y C, Jiang Z Y. A direct swing constraint-based trajectory planning method for underactuated overhead cranes. Acta Automatica Sinica, 2014, 40(11):2414-2419 doi: 10.1016/S1874-1029(14)60397-9 [18] Vaughan J, Kim D, Singhose W. Control of tower cranes with double-pendulum payload dynamics. IEEE Transactions on Control Systems Technology, 2010, 18(6):1345-1358 [19] Singhose W, Kim D. Manipulation with tower cranes exhibiting double-pendulum oscillations. In:Proceedings of the 2007 IEEE International Conference on Robotics and Automation. Roma, Italy:IEEE, 2007. 4550-4555 [20] Masoud Z, Alhazza K, Abu-Nada E, Majeed M. A hybrid command-shaper for double-pendulum overhead cranes. Journal of Vibration and Control, 2014, 20(1):24-37 doi: 10.1177/1077546312461371 [21] 孙宁, 方勇纯, 钱彧哲. 带有状态约束的双摆效应吊车轨迹规划. 控制理论与应用, 2014, 31(7):974-980 http://www.cnki.com.cn/Article/CJFDTOTAL-KZLY201407019.htmSun Ning, Fang Yong-Chun, Qian Yu-Zhe. Motion planning for cranes with double pendulum effects subject to state constraints. Control Theory and Applications, 2014, 31(7):974-980 http://www.cnki.com.cn/Article/CJFDTOTAL-KZLY201407019.htm [22] 郭卫平, 刘殿通. 二级摆型吊车系统动态及基于无源的控制. 系统仿真学报, 2008, 20(18):4945-4948 http://www.cnki.com.cn/Article/CJFDTOTAL-XTFZ200818045.htmGuo Wei-Ping, Liu Dian-Tong. Double-pendulum-type crane dynamics and passivity based control. Journal of System Simulation, 2008, 20(18):4945-4948 http://www.cnki.com.cn/Article/CJFDTOTAL-XTFZ200818045.htm [23] Xu B, Huang X Y, Wang D W, Sun F C. Dynamic surface control of constrained hypersonic flight models with parameter estimation and actuator compensation. Asian Journal of Control, 2014, 16(1):162-174 doi: 10.1002/asjc.2014.16.issue-1 [24] Xu B. Robust adaptive neural control of flexible hypersonic flight vehicle with dead-zone input nonlinearity. Nonlinear Dynamics, 2015, 80(3):1509-1520 doi: 10.1007/s11071-015-1958-8 [25] Garg D, Patterson M, Hager W W, Rao A V, Benson D A, Huntington G T. A unified framework for the numerical solution of optimal control problems using pseudospectral methods. Automatica, 2010, 46(11):1843-1851 doi: 10.1016/j.automatica.2010.06.048 [26] Gill P E, Murray W, Saunders M A. SNOPT:an SQP algorithm for large-scale constrained optimization. SIAM Review, 2005, 47:99-131 doi: 10.1137/S0036144504446096 期刊类型引用(20)
1. 梁越永,晏志维,杨耕,周道付. 基于多运行状态参数切换的冶金无人天车防摇速度控制. 冶金自动化. 2025(01): 1-10 . 百度学术
2. 李刚,马昕,李轾,李贻斌. 带有双球面摆和变绳长效应的桥式起重机轨迹规划. 自动化学报. 2024(05): 911-923 . 本站查看
3. 邵雪卷,张井岗,刘丽琴,陈志梅. 基于T-S模糊模型的桥吊系统积分状态反馈控制. 实验技术与管理. 2023(03): 113-118+124 . 百度学术
4. 魏炳胜,许安定. 基于Kalman滤波和LQG控制的起重机控制研究. 重庆工商大学学报(自然科学版). 2023(03): 20-26 . 百度学术
5. 邱泽昊,孙宁,刘卓清,杨桐,吴庆祥,方勇纯. 双桅杆式起重机的时变输入整形控制方法. 控制理论与应用. 2023(08): 1509-1518 . 百度学术
6. 吴易鸣,孙宁,杨钦朝,尹健宇. 双摆桥式起重机抗扰防摆跟踪控制. 振动与冲击. 2023(21): 36-42 . 百度学术
7. 文学,李华强,朱梅玉,朱科军,刘志辉,陈志刚. 多功能桥式行车实验台控制系统设计. 装备制造技术. 2023(10): 16-20 . 百度学术
8. 王辉堤,肖友刚,李蔚,韩锟. 基于光滑整形自抗扰的双摆起重机消摆控制. 铁道科学与工程学报. 2022(03): 831-840 . 百度学术
9. 石怀涛,黄建起,王延臣. 基于能量强耦合的双摆型塔机消摆控制. 控制工程. 2022(08): 1395-1403 . 百度学术
10. 肖友刚,朱铖臻,李蔚,谢劲松. 基于负载摆动状态估计的双摆型吊车滑模控制. 中南大学学报(自然科学版). 2021(04): 1129-1137 . 百度学术
11. 刘惠康,孙博文,柴琳,鄢梦伟. 分组自适应模糊神经网络在双摆吊车上的应用. 科学技术与工程. 2021(15): 6285-6290 . 百度学术
12. 王岳,孙宁,吴易鸣,梁潇,陈鹤,方勇纯. 深海起重机系统的实时轨迹规划方法. 自动化学报. 2021(12): 2761-2770 . 本站查看
13. 彭海军,施博洋,王昕炜,谢小辉,孙立宁. 考虑区间不确定性的双摆吊车运动轨迹规划. 机械工程学报. 2019(02): 204-213 . 百度学术
14. 欧阳慧珉,王健,张广明,梅磊,邓歆. 基于新型滑模算法的双摆旋转起重机消摆跟踪控制. 自动化学报. 2019(07): 1344-1353 . 本站查看
15. 王杰,强宝民,何祯鑫,杜文正. 基于改进非线性趋近律的桥式起重机滑模控制. 兵器装备工程学报. 2019(07): 148-152+175 . 百度学术
16. 张炜炜,范波,廖志明. 基于平滑过渡轨迹的桥式起重机定位防摆. 计算机仿真. 2019(07): 405-410 . 百度学术
17. 欧阳慧珉,王健,张广明,梅磊,邓歆. 双摆旋转起重机轨迹生成. 控制理论与应用. 2019(08): 1265-1274 . 百度学术
18. 张天成,邵雪卷,张井岗. 基于摆角约束的双摆吊车轨迹规划. 自动化仪表. 2019(09): 40-45 . 百度学术
19. 王杰,强宝民,何祯鑫,董思岐,关标. 欠驱动桥式起重机系统控制技术研究. 兵器装备工程学报. 2019(11): 116-121+159 . 百度学术
20. 罗淑贞,孙青林,檀盼龙,陶金,贺应平,罗浩文. 基于高斯伪谱法的翼伞系统复杂多约束轨迹规划. 航空学报. 2017(03): 220-230 . 百度学术
其他类型引用(23)
-