CN109911773A - 一种单参数调整的欠驱动吊车作业全过程自抗扰控制方法 - Google Patents
一种单参数调整的欠驱动吊车作业全过程自抗扰控制方法 Download PDFInfo
- Publication number
- CN109911773A CN109911773A CN201910306471.8A CN201910306471A CN109911773A CN 109911773 A CN109911773 A CN 109911773A CN 201910306471 A CN201910306471 A CN 201910306471A CN 109911773 A CN109911773 A CN 109911773A
- Authority
- CN
- China
- Prior art keywords
- trolley
- load
- hunting
- control
- crane
- 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.)
- Granted
Links
Landscapes
- Control And Safety Of Cranes (AREA)
- Feedback Control In General (AREA)
Abstract
本发明公布了一种单参数调整的欠驱动吊车作业全过程自抗扰控制方法,其特征在于:根据台车起始及目标位置,考虑台车速度、加速度及作业环境等安全性约束条件,生成效率优先的台车理想作业轨迹,并和台车实际位置进行对比后,构造台车运动的误差反馈控制律;根据吊车负载摆动状态方程设计扩张状态观测器,并构造抑制负载摆动的误差反馈控制律,进而形成不依赖于系统模型参数、且能有效抑制扰动的台车作业轨迹和负载摆动自抗扰控制器;使用Hurwitz稳定矩阵特征值和控制系统增益产生关联,使繁琐的控制系统参数调整转化为很容易实施的单参数调整,实现了在模型参数不确定及外界干扰下吊车能全过程按照设定的理想轨迹运行,并使摆角尽可能小。
Description
技术领域
本发明涉及一种单参数调整的欠驱动吊车作业全过程自抗扰控制方法,属于工程机械自 动控制领域。
背景技术
在港口、电站、地铁、路桥建设、各行业车间等,吊车被广泛应用于搬运重物,其主要 控制要求是将货物安全、快速、准确地搬到目标位置。然而,在货物搬运过程中,易发生摆 动,若到达目标点后残余摆动过大,将导致摆动衰减过程耗时过长,工作效率低,并给司机 的落吊操作带来困难,甚至还可能与周围人员或物体发生碰撞而引发安全事故。现场有经验 的司机可通过控制台车运动来抑制货物摆动,但效率低、技术要求高、还易发生误操作。为 减少货物的残余摆动,现场常采用变频调速方式使货物蚁速慢就位,但低下的效率常遭用户 抱怨。如何保证吊车就位准确、快速,并使货物摆动受到最大抑制并快速衰减,是吊车作业 中尚待解决的重要问题。
吊车系统在工作时,应尽可能快速准确地到达目标位置,并使货物摆幅尽可能小。因此, 对吊车系统进行控制时,需要同时控制吊车位移与货物(负载)摆动,吊车位移可以通过施加 在在台车上的驱动力来控制,而负载的欠驱动特性,使其摆动不能靠驱动力来抑制,只能通 过台车运动来控制,但外界干扰及不当的台车运动又会激发大的负载摆动。因此需要采用科 学的控制方式来控制台车运动,使负载摆幅最小。开环控制和闭环控制是吊车常用的控制方 式,开环控制包括输入整形、离线轨迹规划等,该类方法通过台车运动与负载摆动之间的耦 合关系来规划台车运动,从而实现吊车的防摆、定位,但不能应对外界干扰,鲁棒性差。闭 环控制包括PID控制、最优控制、线性/非线性反馈控制、部分状态反馈控制、跟踪控制、自 适应滑模控制及模糊控制等。这类方法将台车的实时位置或负载的实时摆角等作为反馈信号, 根据理想值和实际值的差异来规划台车的运动,以提高系统的鲁棒性及抗干扰能力。然而, 这些方法大都存在如下的一种或多种缺陷:①台车的工作起点默认为零,目标点的值为初始 误差,当目标点远离初始点时,初始控制量将很大,导致负载初始摆幅过大;②难以考虑吊 绳绳长、负载质量、摩擦系数、初始条件等不确定性因素及外界干扰对台车定位及负载摆动 控制的影响;③需要全状态反馈,但负载摆动角速度难以直接测量,且安装速度传感器会增 加成本,使系统结构复杂;④针对不同塔机需事先设定不同的参考轨迹,并要求参考轨迹满 足一些苛刻的约束条件,难以推广应用;⑤参数整定困难,需要专业人员经过大量试调后 才能确定可行的控制参数;⑥须通过在控制律中添加复杂的非线性耦合项以改善暂态控制性 能,使控制器设计复杂。
发明内容
针对上述问题和不足,为了在模型参数、绳长变化、初始条件等存在不确定性时,使台 车沿满足特定指标的理想轨迹运行,并取得良好的防摆控制效果,本发明公布了一种单参数 调整的欠驱动吊车作业全过程自抗扰控制方法,该方法按照以下步骤实施:
步骤A,根据欠驱动吊车的运动特点,建立欠驱动吊车负载摆动的状态方程,其具体流 程为:
根据由台车、负载、吊绳、轨道组成的三维欠驱动吊车的工作特点,分析吊车系统动态 下的受力状况,并考虑系统未知扰动w(包括摩擦力、模型不确定部分和外扰)对欠驱动吊车 负载摆动的影响,建立吊车运动的拉格朗日方程,在平衡点附近进行线性化处理后,得吊车 的线性化运动微分方程组:
其中,M为台车质量,m为负载质量,θx为吊绳在XZ平面上的投影与Z轴负方向的夹角,θy为吊绳与XZ平面的夹角,l为吊绳长度,fx为X方向的拉力,fy为Y方向的拉力,Dx为X方向的摩擦系数,Dy为Y方向的摩擦系数,Cp为风力系数,Sx为X方向的迎风面积,Sy为Y方向的迎风面积。
为描述方便,将X方向负载摆角θx记为θ1,负载摆动角速度动记为θ2,对X方向负载摆 动的动力学方程(1)进行扩张,得:
其中,θ3为X方向系统总和扰动,u为X方向的小车加速度作为控制量,b为控制量增益。
步骤B、根据台车起始及目标位置,考虑台车的速度、加速度及作业环境等安全性约束 条件,生成效率优先的台车理想作业轨迹,和台车实际位置进行对比后,构造台车作业轨迹 的误差反馈控制律,其具体流程为:
根据环境风速、周围障碍物情况,选择台车的作业方式。在不同作业方式下,为提高作 业效率,设定台车均以该作业方式下的额定速度和加速度运行,具体作业指标为:台车额定 加速度为ae,额定速度为ve,台车作业起点和终点处的速度均为0,台车作业距离为s。根据 台车的作业指标确定台车的理想作业轨迹方程,分如下2种情况:
当s≤ve 2/ae时,台车的理想作业轨迹方程为:
当s>ve 2/ae时,台车的理想作业轨迹方程为:
根据台车的理想作业轨迹方程,可确定任意时刻台车的理想位移信号v1和理想速度信号 v2,而台车的实际位移信号x1和实际速度信号x2很容易测得,则任意时刻台车位置信号与理 想位置信号的误差为ε1=x1-v1,台车速度信号与台车理想速度信号的误差为ε2=x2-v2,以台车 加速度u1为控制量,设计台车作业轨迹的误差控制律:
u1=k1·ε1+k2·ε2=k1·(x1-v1)+k2·(x2-v2) (5)
其中,k1和k2为相应的控制量增益系数。
步骤C、在设计扩张状态观测器观测负载摆动状态的基础上,构造抑制负载摆动的误差 反馈控制律,其具体流程为:
台车运动导致了负载摆动,要抑制负载摆动,就需将摆动的实时值和目标值的差反馈给 台车,并据此施加考虑干扰影响的控制量,使其能通过台车和负载间的耦合运动,在外界干 扰及系统参数不确定情况下抑制负载摆动。为达此目标,需要明确负载摆动的实时状态,但 测量负载摆动角度和角速度不仅会使控制装置复杂,成本增加,而且难以实施,因此,设计 吊车摆动状态的三阶线性扩张状态观测器LESO,实时估计负载摆动角度、角速度及总和扰 动量,具体设计如下:
其中,z1为负载摆角θ1的状态估计,z2为负载摆动角速度θ2的状态估计,z3为负载摆动 状态总和扰动θ3的状态估计,β01,β02,β03为三阶线性扩张状态观测器的参数,b0为扩张状态观测器增益,采用实际控制量增益,即b0=-1/l,并根据实际吊绳长度实时更新。
线性扩张状态观测器的参数β01、β02、β03可根据数值计算步长h确定,即:
吊车作业过程中,要求任意时刻负载摆动的角度和角速度尽可能小,因此其目标值均为 0,据此将抑制负载摆动的误差控制律设计为:
步骤D、将台车运动控制律和负载摆动控制律进行合成,进而设计不依赖于系统模型参 数、结构简单且能有效抑制扰动影响的台车运行轨迹和负载摆动自抗扰控制器,其具体流程 为:
负载摆动状态是通过台车运动来控制的,将台车运动控制量和负载摆动控制量进行合成, 作为台车的综合控制量,可使台车按预定的理想作业轨迹运动,且负载摆动尽可能小,即:
u=u1+u2=k1·(v1-x1)+k2·(v2-x2)+k3·z1+k4·z2-z3/b0 (9)
根据步骤B所述的台车位置控制流程和步骤C所述的负载摆动控制流程,设计吊车作业 全过程自抗扰控制器,该控制器包括台车位置控制环及负载摆动控制环,位置控制环使台车 精准地停在目标位置上,摆动控制环使吊车在运行过程中的摆动趋于0。位置控制环由目标 轨迹自动生成器TTAG和线性反馈控制律组成,TTAG的作用是根据台车起点和终点,考虑 台车安全及运行效率,生成台车运行的理想作业轨迹,台车位置误差反馈控制律根据台车实 时状态和理想作业轨迹的差值来设计,其作用是控制台车精确地按理想作业轨迹运行。负载 摆动控制环由线性扩张状态观测器LESO和负载摆动误差反馈控制律组成,线性扩张状态观 测器的作用是实时估计负载摆动角度、角速度及摆动环的总和扰动,负载摆动误差反馈控制 律根据LESO的观测结果及无负载摆动目标来设计,其作用是在外界干扰及系统不确定情况 下有效抑制负载摆动。最后将台车作业轨迹控制环和负载摆动控制环通过台车运动综合控制 量进行合成,形成欠驱动吊车作业全过程自抗扰控制器。
步骤E、通过Hurwitz稳定矩阵的特征值与控制系统增益产生关联,使繁琐复杂的吊车作 业全过程控制参数调整转化为单参数调整,其具体流程为:
将台车的综合控制量代入到负载摆动的状态方程中,得:
欠驱动吊车的实际目标状态为(x1,x2,θ1,θ2),理想目标状态为(s,0,0,0),其误差可定义为:
ψ1=x1-s,ψ2=x2,ψ3=θ1,ψ4=θ2 (11)
扩张状态观测器的观测误差可定义为:
e1=θ1-z1,e2=θ2-z2,e3=θ3-z3 (12)
令ψ(t)=[ψ1,ψ2,ψ3,ψ4]T,e(t)=[e1,e2,e3]T,ε(t)=[ε1,ε2]T,得控制系统的误差方程:
其中,ψ(t)=[ψ1,ψ2,ψ3,ψ4]T,e(t)=[e1,e2,e3]T,ε(t)=[ε1,ε2]T
其中,Aψ是Hurwitz稳定矩阵,参数k1,k2,k3,k4的选取需符合这一要求。
为使Aψ是Hurwitz稳定矩阵,将Aψ的特征值都配置在点(-ωc,0),ωc>0上,即
|λE-Aψ|=(λ+ωc)4 (14)
求解式(14),可得
由式(15)可知,台车作业轨迹误差反馈控制参数k1,k2和负载摆动误差反馈控制参数k3,k4都由ωc决定,并且当ωc>0时,即可保证矩阵Aψ是Hurwitz稳定矩阵,这样吊车作业全过程 控制的4个待整定参数就转化为ωc一个参数,实现了吊车作业全过程控制的单参数调整,可 大大减少参数整定工作量和难度。
本发明有益效果是:欠驱动吊车作业全过程自抗扰控制通过纠正任一时刻实际位置与理 想位置的偏差,使台车按照设定的轨迹运行,能最大限度地提高吊车作业效率,准确控制台 车的作业时间,为吊车全过程自动作业提供技术支撑;解决了模型参数不确定、外界干扰、 大初始摆角及台车长距离作业情况下欠驱动吊车的防摆定位问题,使吊车在系统模型参数不 确定性及外界环境的干扰下,实际作业曲线都能全过程按照设定的理想轨迹运行,实现吊车 长距离作业时任意位置处的摆动状态控制;实现了吊车作业全过程控制的单参数调整,显著 减少了参数调整工作量和难度。
附图说明
图1为三维欠驱动吊车系统示意图;
图2为二段式台车运行速度-时间曲线;
图3为三段式台车运行速度-时间曲线;
图4为欠驱动吊车作业全过程自抗扰控制框图;
图5为欠驱动吊车长距离作业时本发明的实验结果;
图6为欠驱动吊车在不同初始摆角下本发明的实验结果;
图7为欠驱动吊车吊绳长度变化时本发明的实验结果;
图8为欠驱动吊车在不同脉冲信号干扰下本发明的实验结果。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明了,下面结合附图与具体实施方式,对 本发明进行详细描述。
本发明的欠驱动吊车作业全过程自抗扰控制方法,其基本思路是:首先根据欠驱动吊车 运动特点,建立欠驱动吊车的状态方程,然后,根据台车起始及目标位置,考虑台车的速度、 加速度及作业环境等安全性约束条件,生成效率优先的台车理想作业轨迹,并和台车实际位 置进行对比后,构造台车运动的误差反馈控制律。考虑绳长、负载、空气阻力、系统参数等 不确定性因素影响负载摆动,且负载摆动角度和角速度难以测量的问题,设计扩张状态观测 器,在此基础上构造负载摆动误差反馈控制律,进而形成不依赖于系统模型参数、结构简单 且考虑扰动影响的台车运行轨迹和负载摆动自抗扰控制器;并使用Hurwitz稳定矩阵特征值 和控制系统增益产生关联,将繁多的控制系统参数调整转化为很容易实施的单参数调整,实 现吊车作业全过程自抗扰控制。
本发明的欠驱动吊车作业全过程自抗扰控制方法,按照以下步骤实施:
步骤A,根据欠驱动吊车运动特点,建立欠驱动吊车的状态方程
根据由台车、吊绳、负载、轨道组成的三维欠驱动吊车的工作特点,画出吊车系统模型, 如图1所示。图中,M为台车质量,m为负载质量,θx为吊绳在XZ平面上的投影与Z轴负方向所成的夹角,θy为吊绳与XZ平面的夹角,l为吊绳长度,(xm,ym,zm)表示负载质心坐标。根据图1,(xm,ym,zm)可表示为(x+lsinθxcosθy,y+lsinθy,zm-lcosθxcosθy)。
吊车系统的动能T等于吊车动能与负载动能之和,即:
系统势能V等于吊车势能与负载势能之和,设吊车所在平面为零势能面,得:
V=-mglcosθxcosθy (2)
假设台车摩擦力与台车速度成正比,负载所受空气阻力与负载速度成正比,Dx,Dy分别为X,Y方向的摩擦系数,Cp为空气阻力系数,Sx,Sy分别为X,Y方向负载的迎风面积,则 系统广义坐标X方向的广义力为Y方向的广义力为
令qi=(x,θx,y,θy),Qi表示广义力(i=x,y),L=T-V,可得吊车运动的拉格朗日方程:
由于在实际控制中吊车运动的加速度远小于重力加速度g,且吊车在运动过程中吊绳的 摆角较小(一般不超过10°),对吊车运动的拉格朗日方程展开后,在平衡点附近进行线性化处 理,并考虑系统未知扰动w(包括摩擦力、模型不确定部分和外扰)对欠驱动吊车负载摆动的 影响,得吊车的线性化微分方程组:
从式(4)可以看出,台车X方向与Y方向的动力学方程形式完全相同,只需设计X方向的 控制器即可。为描述方便,将X方向负载摆角θx记为θ1,负载摆动角速度动记为θ2,对X方 向负载摆动的动力学方程(6)进行扩张,得:
其中,θ3为X方向系统总和扰动,u为X方向的小车加速度作为控制量,b为控制量增益。
步骤B、根据台车起始及目标位置,考虑台车的速度、加速度及作业环境等安全性约束 条件,设计台车理想作业轨迹自动生成器TTAG,生成效率优先的台车理想作业轨迹,并和台 车实际位置进行对比后,设计台车作业轨迹的误差反馈控制律;
吊车作业环境复杂,在不同场合,环境风速、周围障碍物差异很大,因此应根据操作规 程,在保障安全的前提下,设定台车的额定运行速度及加速度:风力在3级至5级之间或周 围障碍密集的复杂作业环境,选择低速作业方式;风力在3级至1级之间或周围障碍较多的 普通作业环境,选择中速作业方式;风力在1级以下且无障碍的良好作业环境,选择高速作 业方式。在台车不同作业方式下,为提高作业效率,设定台车均以该作业方式下的额定速度 和加速度运行,具体作业指标为:台车额定加速度为ae,台车额定速度为ve,台车作业起点 和终点处的速度均为0,台车作业距离为s。
在确定了台车运动的额定加速度ae、额定速度ve后,根据起始点和目标点间的距离s, 计算台车以额定加速度运行0.5s所耗的时间t11,并比较台车以ae加速到ve时所经距离D与 0.5s的大小,若D≥0.5s,则说明起始点和目标点的距离s不足以使台车加速到ve,因而台车 以额定加速度运行t11后,就要以ae减速运行至停止,形成如图2所示的二段式台车运行速度 -时间曲线。若D<0.5s,则需计算台车以ae加速到ve所耗的时间t21,及这段时间台车运行的 距离s21,台车加速到ve后,就以ve匀速运行s-2s21距离,再以ae减速运行至停止,形成如 图3所示的三段式台车运行速度-时间曲线。
在第一种情况下,台车在加速段和减速段就走完了全程,这就意味着台车在加速和减速 阶段运行的时间相同、距离均为0.5s、且加速到额定速度ve时所运行的距离不少于0.5s,即:
解得:
因此,当s≤ve 2/ae时,台车的理想轨迹方程只包含加速段和减速段,其数学表达式为:
在第二种情况下,台车以ae加速到ve时,台车的位移还不到0.5s,因此需要以额定速度 ve匀速运动一段时间,在快要到达目标位置时,以额定加速度ae做减速运动,即:
解得:
因此,当s>ve 2/ae时,台车的理想轨迹方程包含加速段、匀速段和减速段,其数学表达 式为:
根据式(8)、式(11)、及台车的起始位置、目标位置、额定加速度、额定速度,设计台车 全过程理想作业轨迹生成器(TTAG)。设t时刻TTAG生成的台车理想位置信号为v1,速度信 号为v2,用相关传感器测定的台车X方向的实际位移信号为x1,速度信号为x2,则台车实际 位置信号与理想位置信号的误差为:
ε1=x1-v1 (12)
台车实际速度信号与理想速度信号的误差为:
ε2=x2-v2 (13)
为使台车精确地按照设定的理想作业轨迹达到目标位置,应使台车在任意时刻的实际值 和理想作业轨迹的误差尽可能小。以台车加速度u1为控制量,设计吊车作业轨迹的误差反馈 控制律:
u1=k1·ε1+k2·ε2=k1·(x1-v1)+k2·(x2-v2) (14)
其中,k1和k2相应的控制量增益系数。
从台车理想作业轨迹的生成过程看出,在任意时刻t,台车都有一个理想的位置信号和速 度信号,通过误差反馈,实时纠正任意时刻实际值与理想值的偏差,控制台车沿着理想轨迹 平稳、准时地达到目标位置,实现对台车运动的全过程控制。相比仅仅根据台车实际位置和 目标位置的轨迹跟踪控制方法,其优势在于:台车作业距离长,运动平稳、时间确定,消除 了目标位置距离初始位置过长而导致的初始控制量过大问题,能为吊车全过程自动作业提供 技术支撑。
步骤C、考虑吊绳绳长、负载、空气阻力、系统参数等不确定性因素影响负载摆动,且 负载摆动角度和角速度难以测量的问题,根据负载摆动状态方程设计扩张状态观测器,并设 计抑制负载摆动的误差反馈控制律;
台车运动导致了负载摆动,要抑制负载摆动,就需将摆动状态的实时值和目标值的差反 馈给台车,并据此施加考虑干扰影响的控制量,使其能通过台车和负载间的耦合运动,在外 界干扰及系统不确定情况下抑制负载摆动。为达此目标,需要明确负载摆动的实时状态,但 测量负载摆动角度和角速度不仅会使控制装置复杂,成本增加,而且难以实施,因此,设计 负载摆动线性扩张状态观测器LESO,实时估计负载摆动角度、角速度,及系统总和扰动量。
具体设计过程如下:
其中,z1是负载摆角θ1的状态估计,z2是对负载摆动角速度θ2的状态估计,z3是对负载 摆动状态总和扰动θ3的状态估计,包括系统未建模部分、不确定参数及外部干扰,β01,β02, β03是三阶线性扩张状态观测器的参数,b0是扩张状态观测器增益,是控制量增益b的估计值。
从式(5)可知,负载摆动的控制量增益为:
因此,负载摆动的控制量增益可以用实际吊绳长度表示。在吊车的控制过程中,吊绳长 度很容易实测得到,因此扩张状态观测器增益b0可根据实际摆线长度实时更新,即:
在计算机控制中,线性扩张状态观测器的参数β01,β02,β03可根据数值计算采样步长h 确定,即:
吊车作业过程中,要求任意时刻负载摆动的角度和角速度尽可能小,因此其目标值均为 0,据此将抑制负载摆动的误差控制律设计为:
其中,k3和k4相应的控制量增益系数。
步骤D、将台车作业轨迹的误差反馈控制律和抑制负载摆动的误差反馈控制律进行合成, 形成不依赖于系统模型参数、结构简单且能有效抑制扰动影响的台车运行轨迹和负载摆动自 抗扰控制器;
负载摆动状态是通过台车运动来控制的,要使台车按预定的理想作业轨迹运行时,负载 摆动还尽可能小,就应该将台车运动控制量和负载摆动控制量进行合成,作为台车的综合控 制量,即:
u=u1+u2=k1·(v1-x1)+k2·(v2-x2)+k3·z1+k4·z2-z3/b0 (20)
根据步骤B所述的台车位置控制流程和步骤C所述的负载摆动控制流程,设计如图4所 示的吊车作业全过程自抗扰控制器,该控制器包括台车位置控制环及负载摆动控制环,位置 控制环使吊车精准地停在目标位置上,摆动控制环使吊车在运行过程中的摆动趋于0。作业 轨迹控制环由目标轨迹自动生成器TTAG和线性反馈控制律组成,TTAG的作用是根据台车 起点和终点,考虑台车安全及运行效率,生成台车的理想作业轨迹,台车位置误差反馈控制 律根据台车实时状态和理想作业轨迹的差值来设计,其作用是控制台车精确地按理想作业轨 迹运行。负载摆动控制环由线性扩张状态观测器LESO和负载摆动误差反馈控制律组成,线 性扩张状态观测器的作用是实时估计负载摆动角度、角速度及总和扰动,负载摆动误差反馈 控制律根据LESO的观测结果及无负载摆动目标来设计,其作用是在外界干扰及系统不确定 情况下有效抑制负载摆动。最后将台车作业轨迹控制环和负载摆动控制环通过台车运动综合 控制量进行合成,形成欠驱动吊车作业全过程自抗扰控制器。
步骤E、使Hurwitz稳定矩阵的特征值和控制系统增益产生关联,将繁琐的控制系统参数 调整转化为很容易实施的单参数调整;
吊车作业全过程自抗扰控制器包括两个误差反馈控制器,待整定的参数包括k1、k2、k3、 k4,且位置控制环和摆动控制环相互耦合,靠试凑的方法显然不合理,且会耗费大量的时间 和精力。由于自抗扰控制参数的取值必须满足系统误差收敛这一基本条件,因此从这一基本 条件出发,来分析如何进一步精简并选取合适的控制参数。
将式(20)代入式(5)中可得:
欠驱动吊车的理想目标状态为(x1,x2,θ1,θ2)=(s,0,0,0),其误差可定义为:
ψ1=x1-s,ψ2=x2,ψ3=θ1,ψ4=θ2 (22)
扩张状态观测器的观测误差可定义为:
e1=θ1-z1,e2=θ2-z2,e3=θ3-z3 (23)
台车理想轨迹与目标状态的误差可定义为:
ε1=v1-s,ε2=v2 (24)
由式(8)及式(11)可知,ε=[ε1,ε2]T严格收敛于0。
令ψ=[ψ1,ψ2,ψ3]T,e=[e1,e2,e3]T,将式(22)(23)(24)代入式(21)得,整个控制系统的误差方 程:
其中,
要使系统误差收敛,Aψ必须是Hurwitz稳定矩阵,因此将Aψ的特征值都配置在点(-ωc,0), ωc>0上,即
|λE-Aψ|=(λ+ωc)4 (26)
求解式(26)可得
由式(27)可知,台车作业轨迹误差反馈控制参数k1,k2和负载摆动误差反馈控制参数k3,k4都由Hurwitz稳定矩阵的特征值-ωc决定,并且当ωc>0时,即可保证矩阵Aψ是Hurwitz稳定 矩阵,保障了系统误差的收敛性,且吊车作业全过程控制的4个待整定参数就转化为ωc一个 参数,实现了吊车作业全过程控制的单参数调整,可大大减少参数整定工作量和难度。
步骤F、使用Lyapunov函数证明吊车作业全过程自抗扰控制器的稳定性,具体过程如下:
假设1变量x1,x2,θ1,θ2,θ3在作业时间内是连续可微的。
假设2系统总和扰动及其导数是有界的,即:(δ1,δ2均为正实数)
将扩张状态观测器的误差e对时间t求导,再将式(15)代入,可得线性扩张状态观测器误 差的微分方程为:
式中,
矩阵E的特征值分别为当h>0时,其特征值的实部均小于0,则矩阵E为Hurwitz稳定矩阵。对于任意给定的对称 正定矩阵Q,存在对称正定矩阵P满足Lyapunov方程:
ETP+PE+Q=0 (36)
将扩张状态观测器的Lyapunov函数定义为:
V0=eTPe (37)
则
其中,λmin(Q)是矩阵Q的最小特征值。
当δ2=0时,因此线性扩张状态观测器LESO在零点大范围渐进稳定。当δ2≠0时, 系统的观测值会出现一定误差。由于系统达到稳态时有根据式(27)得到观 测值的误差范围:
将式(25)进一步表示为:
其中,g(e,ε,w)=Aee+Aεε+Flw。
因为矩阵Aψ是Hurwitz稳定矩阵,令式(40)的Lyapunov函数为:
V3=ψTAψψ (41)
则:
当扰动w=0时,因此,系统在稳定点附近大范围稳定。
当w≠0时,
其中,ζ1=-k3e1-k4e2-le3+w1,
当系统稳定时,根据式(28),可得到系统的误差范 围:
根据欠驱动吊车的实际情况,当系统稳定时,负载摆角会在自身重力作用下回到竖直向 下的位置,因此欠驱动吊车的摆角不可能存在稳态误差,即有ψ3=0,因此式(46)所示的欠驱 动吊车控制系统的稳态误差可更新为:
从式(47)可以看出,当扩张状态观测器的观测值与实际值没有误差,或观测误差足够小 时,欠驱动吊车位置环的稳态误差也可以足够小。
步骤G、通过实验说明吊车定位消摆全过程自抗扰控制器有效
为了验证欠驱动吊车定位消摆全过程自抗扰控制器的有效性,设计一实验平台,其参数 为:M=6.157kg,m=1kg,g=9.81m/s2,CP=0.5N·s·m-3,S=0.0036m2,ae=0.05m/s2,ve=1m/s。 实验中,将台车目标位置设定为100m,计算机步长取0.05,经过充分整定后,式(31)中自抗 扰控制参数ωc=4。
为验证本发明的控制效果,设计绳长l保持0.56m不变的仿真实验,得到的实验结果如 图5所示,图中的曲线自上至下分别表示台车位移x(t)、负载摆角θ(t)。从台车的位移曲线可 以看出,台车够精确按照设定的理想曲线运动,且在规定的时间内到达目标点。从负载摆角 曲线可以看出,台车加速和减速过程中负载都会摆动,但最大摆角只有0.4°,且在加速或减 速结束后负载摆角迅速衰减为0,在台车匀速平移过程中,负载几乎无摆动。说明在负载长 距离运送过程中,本发明能使台车按设定的理想曲线将负载平稳运送到目的地,基本无残余 摆动,且运送时间也和设定时间保持一致,为120s。
吊车搬运大件物体时,物体重心很难控制在吊钩正下方,导致物体起吊时会产生初始摆 角,为分析负载初始摆角变化时,本发明对吊车运动的控制效果,设计台车目标位置为100m, 负载初始摆角为2°、4°、6°三组仿真实验,所得实验结果分别用图6中的实线(2°)、虚线(4°)、 点划线(6°)表示。从图中曲线可以看出,在不同初始摆角下,台车都能够精确按照设定的理想 曲线运动,且到达目标点的时间均为120s,受不同负载初始摆角的影响,0~3s内负载摆角 有较大区别,但随后的负载摆角趋于一致,且被抑制在很小范围内,当台车到达目标位置后 负载几乎无残余摆动,说明本发明对不确定的初始条件有很强的鲁棒性。
为提高系统工作效率,在有些情况下会同步进行负载升降与水平运送动作,此时,吊绳 长度从常值转变为变量,为检验本发明对绳长变化的敏感性,让负载在5~10s内从0.5m上升 到2.5m,在50~70s内从2.5m上升到10.5m,在95~120s从10.5m下降到0.5m进行仿真实验, 得到如图7所示的实验结果。从图中可以看出,虽然在负载上升和下降过程中,吊绳长度经 历了缩短、变长和不变三种状态,但台车还是能够精确按照设定的理想曲线运动,且到达目 标点的时间仍为120s,在台车加、减速过程中,负载的升、降运动虽然使负载摆角出现了一 定波动,但波动幅度极小,在台车匀速平移过程中,负载的升、降运动对负载摆角几乎没有 影响,说明本发明对绳长变化有极强的鲁棒性。
吊车工作环境复杂,周围存在各种各样的干扰,为检验本发明的抗干扰能力,在0s时对 摆角环施加强度为1的脉冲信号(pulse),在50s时对摆角环施加强度为5的脉冲信号(pulse2), 在105s时对台车位置环施加强度为1的脉冲信号(pulse3),在130s时对台车位置环施加强度 为5的脉冲信号(pulse4),进行仿真实验,得到的实验结果如图8所示。从图中可以看出,在 脉冲信号干扰下,台车仍能够按照设定的理想曲线运动,且到达目标点的时间仍为120s,在 脉冲信号作用处,负载摆角出现了一定幅度波动,但波动幅度仍然很小,说明本发明具有良 好的抗干扰能力。
如上所述,结合附图和说明所给出的方案内容,可以衍生出类似的技术方案。但凡是依 据本发明的技术实质所作的任何简单修改、等同变化与修饰,均仍属于本发明技术方案的范 围内。
Claims (5)
1.一种单参数调整的欠驱动吊车作业全过程自抗扰控制方法,其特征在于,该方法包括:
步骤A、根据欠驱动吊车运动特点,建立欠驱动吊车的动力学方程;
步骤B、根据台车起始及目标位置,考虑台车的速度、加速度及作业环境等安全性约束条件,设计台车理想轨迹自动生成器TTAG,生成效率优先的台车理想作业轨迹,并和台车实际位置进行对比后,设计台车作业轨迹的误差反馈控制律;
步骤C、考虑绳长、负载、空气阻力、系统参数等不确定性因素影响负载摆动,且负载摆动角度和角速度难以测量的问题,根据负载摆动的状态方程设计实时估计负载摆动状态的扩张状态观测器,并设计抑制负载摆动的误差反馈控制律;
步骤D、将台车作业轨迹的误差反馈控制律和抑制负载摆动的误差反馈控制律进行合成,得出台车运动的综合控制律,并形成不依赖于系统模型参数、结构简单且能有效抑制扰动影响的台车作业轨迹和吊车摆动自抗扰控制器;
步骤E、使用Hurwitz稳定矩阵的特征值和控制系统增益产生关联,将繁琐的控制系统参数调整转化为很容易实施的单参数调整。
2.根据权利要求1所述的一种单参数调整的欠驱动吊车作业全过程自抗扰控制方法,其特征在于,在所述步骤B中,台车理想作业轨迹及台车作业轨迹误差反馈控制律的具体设计过程为:
根据环境风速、周围障碍物情况,选择台车的作业方式,在不同作业方式下,为提高作业效率,设定台车均以该作业方式下的额定速度和额定加速度运行,具体作业指标为:台车额定加速度为ae,额定速度为ve,台车作业起点和终点位置处的速度均为0,台车作业距离为s,根据台车的作业指标确定台车的理想作业轨迹方程,分如下二种情况:
当s≤ve 2/ae时,台车的理想作业轨迹方程为:
当s>ve 2/ae时,台车的理想作业轨迹方程为:
根据台车的理想作业轨迹方程,确定任意时刻台车的理想位移信号v1和理想速度信号v2,而台车的实际位移信号x1和实际速度信号x2很容易测得,则任意时刻台车实际位置信号与理想位置信号的误差为ε1=x1-v1,实际速度信号与理想速度信号的误差为ε2=x2-v2,以台车加速度u1为控制量,设计台车作业轨迹的误差反馈控制律:
u1=k1·ε1+k2·ε2=k1·(x1-v1)+k2·(x2-v2) (3)
式中,k1和k2为相应的控制量增益系数。
3.根据权利要求1所述的一种单参数调整的欠驱动吊车作业全过程自抗扰控制方法,其特征在于,在所述步骤C中,负载摆动状态的扩张状态观测器及误差反馈控制律的具体设计过程为:
根据欠驱动吊车的动力学方程得出负载摆动的状态方程:
其中,θ1为吊绳X方向摆角,θ2为吊绳X方向摆动角速度,θ3为系统总和扰动,w1为系统扰动,u为X方向的小车加速度作为控制量,b为控制量增益;
根据负载摆动状态方程设计线性扩张状态观测器LESO:
其中,z1是负载摆角θ1的状态估计,z2是负载摆动角速度θ2的状态估计,z3是负载摆动状态总和扰动θ3的状态估计,β01,β02,β03是三阶线性扩张状态观测器参数,b0是扩张状态观测器增益,采用实际控制量增益,即b0=-1/l,并根据实际绳长实时更新;
线性扩张状态观测器的参数β01、β02、β03根据数值计算步长h确定,即:
吊车作业过程中,要求任意时刻负载摆动的角度和角速度尽可能小,因此其目标值均为0,据此将抑制负载摆动的误差控制律设计为:
4.根据权利要求1所述的一种单参数调整的欠驱动吊车作业全过程自抗扰控制方法,其特征在于,在所述步骤D中,所述的台车运动综合控制律及自抗扰控制器的设计过程为:
将台车作业轨迹的误差反馈控制律和抑制负载摆动的误差反馈控制律进行合成,作为台车运动综合控制量:
u=u1+u2=k1·(v1-x1)+k2·(v2-x2)+k3·z1+k4·z2-z3/b0 (8)
吊车作业全过程自抗扰控制器设计包括台车作业轨迹控制环设计及负载摆动控制环设计,作业轨迹控制环由目标轨迹自动生成器TTAG和线性反馈控制律组成,TTAG的作用是根据台车起点和终点,考虑台车安全及运行效率,生成台车的理想作业轨迹,台车位置误差反馈控制律根据台车实时状态和理想轨迹的差值来设计,其作用是控制台车精确地按理想作业轨迹运行;负载摆动控制环由线性扩张状态观测器LESO和负载摆动误差反馈控制律组成,线性扩张状态观测器的作用是实时估计负载摆动角度、角速度及摆动环的总和扰动,负载摆动误差反馈控制律根据LESO的估计结果及无负载摆动目标来设计,其作用是在外界干扰及系统参数不确定情况下有效抑制负载摆动;将台车作业轨迹控制环和负载摆动控制环通过台车运动综合控制量进行合成,形成欠驱动吊车作业全过程自抗扰控制器。
5.根据权利要求1所述的一种单参数调整的欠驱动吊车作业全过程自抗扰控制方法,其特征在于,在所述步骤E中,所述的单参数调整是将将欠驱动吊车待整定的参数k1、k2、k3、k4精简成一个,具体过程为:
将台车的综合控制量代入到负载摆动的状态方程中,可得:
欠驱动吊车的理想目标状态可表示为(x1,x2,θ1,θ2)=(s,0,0,0),其误差可定义为:
ψ1=x1-s,ψ2=x2,ψ3=θ1,ψ4=θ2 (10)
扩张状态观测器的观测误差可定义为:
e1=θ1-z1,e2=θ2-z2,e3=θ3-z3 (11)
台车理想轨迹与目标状态的误差可定义为:
ε1=v1-s ε2=v2 (12)
令ψ(t)=[ψ1,ψ2,ψ3,ψ4]T,e(t)=[e1,e2,e3]T,ε(t)=[ε1,ε2]T,得整个控制系统的误差方程:
其中,
式中,Aψ是Hurwitz稳定矩阵,参数k1,k2,k3,k4的选取需符合这一要求;
为使Aψ是Hurwitz稳定矩阵,将Aψ的特征值都配置在点(-ωc,0),ωc>0上,即
|λE-Aψ|=(λ+ωc)4 (14)
求解式(14)可得:
因此,台车作业轨迹环误差反馈控制参数k1,k2和负载摆动误差反馈控制参数k3,k4都由ωc决定,并且当ωc>0时,即可保证矩阵Aψ是Hurwitz稳定矩阵,这样吊车作业全过程控制的四个待整定参数就转化为ωc一个参数,实现了吊车作业全过程控制的单参数调整,大大减少了参数整定工作量和难度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910306471.8A CN109911773B (zh) | 2019-04-17 | 2019-04-17 | 一种单参数调整的欠驱动吊车作业全过程自抗扰控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910306471.8A CN109911773B (zh) | 2019-04-17 | 2019-04-17 | 一种单参数调整的欠驱动吊车作业全过程自抗扰控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109911773A true CN109911773A (zh) | 2019-06-21 |
CN109911773B CN109911773B (zh) | 2020-09-01 |
Family
ID=66977568
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910306471.8A Active CN109911773B (zh) | 2019-04-17 | 2019-04-17 | 一种单参数调整的欠驱动吊车作业全过程自抗扰控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109911773B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111176318A (zh) * | 2020-02-10 | 2020-05-19 | 南开大学 | 面向消除稳态误差的飞行吊运系统定位消摆方法及系统 |
CN111762688A (zh) * | 2020-07-02 | 2020-10-13 | 北京化工大学 | 桥式吊车控制器生成方法、控制方法以及控制器生成系统 |
CN113093541A (zh) * | 2021-03-31 | 2021-07-09 | 中南大学 | 一种欠驱动吊车微分平坦跟踪控制方法 |
CN113126502A (zh) * | 2021-04-28 | 2021-07-16 | 天津理工大学 | 一种欠驱动吊车系统的控制方法及控制系统 |
WO2021196937A1 (zh) * | 2020-04-02 | 2021-10-07 | 上海驭矩信息科技有限公司 | 一种基于lqr的吊装系统的防摇控制方法及系统 |
CN113608440A (zh) * | 2021-08-04 | 2021-11-05 | 哈尔滨理工大学 | 一种考虑绳长变化的船用吊艇系统减摆控制方法 |
CN113896111A (zh) * | 2021-09-28 | 2022-01-07 | 中国人民解放军火箭军工程大学 | 一种欠驱动起重吊装设备控制系统及其方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101414156A (zh) * | 2008-11-18 | 2009-04-22 | 哈尔滨工业大学 | 欠驱动机械装置acrobot的动态伺服控制方法 |
CN102795544A (zh) * | 2012-08-16 | 2012-11-28 | 南开大学 | 基于轨迹在线规划的桥式吊车高效消摆控制方法 |
CN105152020A (zh) * | 2015-09-30 | 2015-12-16 | 山东大学 | 带有跟踪误差约束的桥式吊车自适应跟踪控制器及方法 |
CN108358062A (zh) * | 2018-03-02 | 2018-08-03 | 上海应用技术大学 | 欠驱动吊车全局稳定控制方法 |
CN108875253A (zh) * | 2018-07-03 | 2018-11-23 | 曲阜师范大学 | 基于干扰观测器的欠驱动吊车系统的终端滑模消摆控制方法及系统 |
-
2019
- 2019-04-17 CN CN201910306471.8A patent/CN109911773B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101414156A (zh) * | 2008-11-18 | 2009-04-22 | 哈尔滨工业大学 | 欠驱动机械装置acrobot的动态伺服控制方法 |
CN102795544A (zh) * | 2012-08-16 | 2012-11-28 | 南开大学 | 基于轨迹在线规划的桥式吊车高效消摆控制方法 |
CN105152020A (zh) * | 2015-09-30 | 2015-12-16 | 山东大学 | 带有跟踪误差约束的桥式吊车自适应跟踪控制器及方法 |
CN108358062A (zh) * | 2018-03-02 | 2018-08-03 | 上海应用技术大学 | 欠驱动吊车全局稳定控制方法 |
CN108875253A (zh) * | 2018-07-03 | 2018-11-23 | 曲阜师范大学 | 基于干扰观测器的欠驱动吊车系统的终端滑模消摆控制方法及系统 |
Non-Patent Citations (4)
Title |
---|
李树江等: "基于逆系统方法的内模控制在吊车消摆运动中的应用", 《信息与控制》 * |
王鹏程等: "一种直接基于摆角约束的欠驱动桥式吊车轨迹规划方法(英文)", 《自动化学报》 * |
谭莹莹等: "基于动态滑模结构的桥式吊车防摇定位控制器设计", 《控制工程》 * |
马博军等: "欠驱动桥式吊车系统自适应控制", 《控制理论与应用》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111176318A (zh) * | 2020-02-10 | 2020-05-19 | 南开大学 | 面向消除稳态误差的飞行吊运系统定位消摆方法及系统 |
CN111176318B (zh) * | 2020-02-10 | 2022-12-06 | 南开大学 | 面向消除稳态误差的飞行吊运系统定位消摆方法及系统 |
WO2021196937A1 (zh) * | 2020-04-02 | 2021-10-07 | 上海驭矩信息科技有限公司 | 一种基于lqr的吊装系统的防摇控制方法及系统 |
US11708248B2 (en) | 2020-04-02 | 2023-07-25 | Shanghai Master Matrix Information Technology Co., Ltd. | LQR-based anti-sway control method and system for lifting system |
CN111762688A (zh) * | 2020-07-02 | 2020-10-13 | 北京化工大学 | 桥式吊车控制器生成方法、控制方法以及控制器生成系统 |
CN111762688B (zh) * | 2020-07-02 | 2021-06-22 | 北京化工大学 | 桥式吊车控制器生成方法、控制方法以及控制器生成系统 |
CN113093541A (zh) * | 2021-03-31 | 2021-07-09 | 中南大学 | 一种欠驱动吊车微分平坦跟踪控制方法 |
CN113126502A (zh) * | 2021-04-28 | 2021-07-16 | 天津理工大学 | 一种欠驱动吊车系统的控制方法及控制系统 |
CN113608440A (zh) * | 2021-08-04 | 2021-11-05 | 哈尔滨理工大学 | 一种考虑绳长变化的船用吊艇系统减摆控制方法 |
CN113608440B (zh) * | 2021-08-04 | 2022-02-22 | 哈尔滨理工大学 | 一种考虑绳长变化的船用吊艇系统减摆控制方法 |
CN113896111A (zh) * | 2021-09-28 | 2022-01-07 | 中国人民解放军火箭军工程大学 | 一种欠驱动起重吊装设备控制系统及其方法 |
CN113896111B (zh) * | 2021-09-28 | 2022-07-29 | 中国人民解放军火箭军工程大学 | 一种欠驱动起重吊装设备控制系统及其方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109911773B (zh) | 2020-09-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109911773A (zh) | 一种单参数调整的欠驱动吊车作业全过程自抗扰控制方法 | |
CN110436347B (zh) | 一种桥式起重机双摆系统优化防摇控制方法 | |
CN107215777B (zh) | 一种起重机智能防摇摆控制系统及其精确定位方法 | |
CN111153328B (zh) | 一种基于lqr的吊装系统的防摇控制方法及系统 | |
CN106959610B (zh) | 桥式吊车系统apd-smc控制器、桥式吊车系统及控制方法 | |
CN105934401B (zh) | 用于控制从电动可滑动件悬吊载荷的摆动设备、工序和起重装置的控制装置 | |
CN101384503B (zh) | 吊装载荷摆动防止装置 | |
CN107487717A (zh) | 基于模糊pid控制的起重机智能防摇控制系统及方法 | |
CN110407094B (zh) | 一种基于动态平滑轨迹的桥式起重机定位防摆控制方法 | |
CN110316657A (zh) | 一种起重机重物的防摇摆控制系统及其控制方法 | |
CN108584700B (zh) | 起重机自适应pid闭环防摇控制方法 | |
CN107826978B (zh) | 一种双摆桥式起重机的速度轨迹规划消摆方法 | |
CN110407095A (zh) | 一种基于在线轨迹规划的桥式起重机定位消摆控制方法 | |
CN106865416A (zh) | 一种桥式起重机行走过程中的自动防摆控制方法 | |
CN112429640B (zh) | 一种控制天车防摇的方法 | |
CN206172785U (zh) | 风力发电机组吊具 | |
CN113321122B (zh) | 起重机双摆系统分层非奇异终端滑模控制方法 | |
CN114594683A (zh) | 基于Hurwitz稳定的动基座桥式起重机防摆滑模控制方法 | |
CN102530729B (zh) | 控制吊物振摆的方法和系统 | |
CN210620028U (zh) | 一种起重机重物的防摇摆控制系统 | |
CN103274298B (zh) | 控制塔式起重机小车驱动力的方法、装置和系统 | |
RU2744647C1 (ru) | Способ адаптивного управления мостовым краном | |
Wang et al. | A new approach of anti-swing control system based on run-to-run control and fuzzy control for overhead crane | |
Li et al. | Experimental study on rapid oscillation suppression of a pendulum using phase delay motion of pivot | |
Prabaningtyas et al. | Linear Quadratic Gaussian (LQG) for stability control of single payload overhead crane system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |