CN115309044A - 一种基于模型预测控制的机械臂角速度控制方法 - Google Patents
一种基于模型预测控制的机械臂角速度控制方法 Download PDFInfo
- Publication number
- CN115309044A CN115309044A CN202210887459.2A CN202210887459A CN115309044A CN 115309044 A CN115309044 A CN 115309044A CN 202210887459 A CN202210887459 A CN 202210887459A CN 115309044 A CN115309044 A CN 115309044A
- Authority
- CN
- China
- Prior art keywords
- control
- signal
- matrix
- mode
- state
- 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
- 238000000034 method Methods 0.000 title claims abstract description 23
- 239000011159 matrix material Substances 0.000 claims description 41
- 230000004044 response Effects 0.000 claims description 11
- 230000010355 oscillation Effects 0.000 claims description 7
- 230000008859 change Effects 0.000 claims description 6
- 230000001133 acceleration Effects 0.000 claims description 3
- 238000009795 derivation Methods 0.000 claims description 3
- 239000007788 liquid Substances 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 239000003921 oil Substances 0.000 description 5
- 238000013016 damping Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000013178 mathematical model Methods 0.000 description 3
- 238000006073 displacement reaction Methods 0.000 description 2
- 239000010720 hydraulic oil Substances 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 230000001052 transient effect Effects 0.000 description 2
- 101100346189 Caenorhabditis elegans mpc-1 gene Proteins 0.000 description 1
- 101710165595 Mitochondrial pyruvate carrier 2 Proteins 0.000 description 1
- 102100025031 Mitochondrial pyruvate carrier 2 Human genes 0.000 description 1
- 101710101698 Probable mitochondrial pyruvate carrier 2 Proteins 0.000 description 1
- 230000002457 bidirectional effect Effects 0.000 description 1
- 230000009123 feedback regulation Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
Abstract
本发明公开了一种基于模型预测控制的机械臂角速度控制方法,包括以下步骤:S1、获取***的指令信号ucmd(t);S2、对***进行分段线性化近似划分,构造不同线性化阶段的模型预测控制器;S3、依据预测的下一阶段的指令信号ucmd(k+1)‑uout(k+1)对模型预测控制器进行逻辑切换;S4、获得压力反馈信号,所述压力反馈信号经模型预测控制器优化得到输出信号uout(t),将输出信号uout(t)作用于指令信号ucmd(t),得到优化后的输入信号u(t);S5、基于优化后的输入信号u(t)对直驱液压机械臂进行控制。
Description
技术领域
本发明涉及机械臂的动态性能改善技术领域,更具体的说是涉及一种基于模型预测控制的机械臂角速度控制方法。
背景技术
目前液压挖掘机、起重机、高空作业等机械臂大都采用阀控***,在运行过程中存在较大的泄漏、压力超调、压力波动且能量损失较大,极大影响机械臂的动态性能;其他控制方法都是非线性***直接进行反馈调节,存在较大且不可避免的非线性扰动。
因此,针对机械臂液压执行器在运行过程中存在的振荡问题,通过对其***非线性模型在不同工作点处线性化,进而近似将***由多个分段线性模型表示,最后针对每个线性化模型分别设计不同的增益规划模型预测控制器(MPC),并定义控制器之间切换算法,在不同运行速度和负载条件下降低压力超调量,缩短到达稳定的时间,从而极大改善机械臂的动态性能是本领域技术人员亟需解决的问题。
发明内容
有鉴于此,本发明提供了一种基于模型预测控制的机械臂角速度控制方法;通过当前***状态预测机械臂的未来状态,将增益规划模型预测控制器(MPC)与直驱泵控***相结合,应用于不同液压机械臂,根据不同速度、负载工况下,调节压力波动,进而在机械臂工作中充当振动的阻尼,执行器压力和动臂角度位置通过在适当位置安装传感器进行监控,并用作反馈信号;将非线性模型在不同操作点处线性化,获得一个分段线性模型,以表示整个操作空间内的***;通过分段线性化建立不同增益效果的增益规划模型预测控制器(MPC),设计切换逻辑,进而在不同运行速度和负载条件下的线性化阶段选择合适的增益规划模型预测控制器(MPC)来降低压力超调量,缩短到达稳定的时间,从而极大改善机械臂的动态性能。
为了实现上述目的,本发明采用如下技术方案:
一种基于模型预测控制的机械臂角速度控制方法,包括以下步骤:
S1、获取***的指令信号ucmd(t);
S2、对***进行分段线性化近似划分,构造不同线性化阶段的模型预测控制器;
S3、依据预测的下一阶段的指令信号ucmd(k+1)-uout(k+1)对模型预测控制器进行逻辑切换;
S4、获得压力反馈信号,所述压力反馈信号经模型预测控制器优化得到输出信号uout(t),将输出信号uout(t)作用于指令信号ucmd(t),得到优化后的输入信号u(t)=ucmd(t)-uout(t);
S5、基于优化后的输入信号u(t)对机械臂进行控制。
优选的,所述步骤S2具体包括:
S21、构建状态模型,分段线性化***方程的状态空间形式为:
Xi(k+1)=AiXi(k)+Biu(k);i=(1,2,3…n)
y(k)=CXi(k)
其中,Xi(k+1)为预测***k+1阶段状态向量,Xi(k)是***k阶段状态向量,Ai是***矩阵,Bi是控制输入矩阵,u(k)是输入变量矩阵,y(k)为***输出,C为二阶单位矩阵;
***阶段状态向量表示为:
其中,θ为机械臂运行角度;为机械臂运行角度变化率;pA,pB分别为液压缸有杆腔和无杆腔反馈压力值;p'A,p'B分别为液压缸有杆腔和无杆腔忽略单向阀、管路等压降的两腔压力,理想情况下p'A=pA,p'B=pB;xV为液压缸运行速度;为液压缸运行变化率即加速度;βe为油液有效体积弹性模量;
在不同的线性化阶段考虑不同状态变量的影响,因此根据具体的控制信号大小,对***的状态向量进行划分,根据上述分段线性化***方程,得到增量式方程:
Δx(k+1)=Xi(k+1)-Xi(k)=AiΔX(k)+BiΔu(k)
其中,Δx(k+1)=x(k+1)-x(k),Δu(k+1)=u(k+1)-u(k)
同理,输出的增量方程为:
y(k)=CΔx(k)+y(k-1)
其中,y(k)为***当前阶段输出量,C为二阶单位矩阵,Δx(k)为***状态向量变化量;y(k-1)为上一阶段***输出量;
以最新的测量值为初始条件,预测步长为l,控制步长为m,并假设外部扰动对***影响很小,可得:
Δu(k+i)=0,i=m,m+1,…,l-1;
S22、构建状态预测模型:
其中,Y0为模型初始输入量;Su为***状态预测空间方程,C为二阶单位矩阵;A为***输入矩阵;Bi为控制输入矩阵;Δu(k)为输入控制增量;l为预测步长,m为控制步长
S23、在目标函数建立中,设置针对状态量代价Q矩阵和控制量的代价矩阵R权重,以对控制量、状态量的控制,使***的功能靠近理想目标,通过对目标函数进行求导,得出对于控制量输入的最优解:
o(k+i)=αiyk+(1-αi)yr
其中,o(k+i)为控制量输入最优解,即期望轨迹,αi为期望轨迹的权重因子,其值越小,响应速度越迅速,0<α<1,yr为目标值,yk为k时刻下***输出量;yr为***目标值;
模型预测控制器建立的目标函数为:
其中,rk+i是参考输出,yk+i是预测输出;Δu是控制信号与其标称值的偏差;Q为状态量代价矩阵;R为控制量代价矩阵;
则成本函数为:
其中Q权重大小与控制量有关,R权重大小与使用多少能量达到目标值有关,Δzk+i***输出跟踪与输入之间误差;Δuk+i-1为控制变化量;
控制目标目的在于使压力振荡最小化,同时输出信号与控制输入信号偏离达到最小,若目标函数最小,则其必要条件为:
通过求导,可得到最优解Δu(k),将最优解作为控制增量,得到下一时刻的信号输入;
则补偿函数为:
其中,H为补偿误差权重值,e为误差值;
Ycor=[Ycor(k+1)Ycor(k+2)…Ycor(k+p)]T;H=[h1h2…hp]T
将求得的补偿函数Ycor代入状态预测模型,即可得到新的加入补偿过后的预测函数:
优选的,所述步骤S3具体包括:
定义切换函数s(k),s(k)在闭合区间内取整数值,其表示在第k个阶段处于活动状态的模型预测控制器的控制模式将成本函数j(k)作为预测模型控制器的候选李雅普诺夫函数,只有当候选李雅普诺夫函数在当前模式域外的其他模式下具有较低的值时,切换函数s(k)才改变其值,则模式切换判断方式为:
(1)ucmd(k+1)-uout(k+1)在当前模式的限制范围内,则控制器保持在当前模式,即s(k+1)=s(k);
(2)ucmd(k+1)-uout(k+1)超过当前模式的上限,则切换块在当前模式以上的模式中搜索j(k+1)最小的模式,并选择j(k+1)最小的模式;
(3)ucmd(k+1)-uout(k+1)低于当前模式的下限,则切换块在当前模式以下的模式中搜索j(k+1)最小的模式,并选择j(k+1)最小的模式。
优选的,所述限制范围为根据非线性模型在操作空间得到的响应面,划分不同阶段有着对应的控制命令区间范围;所以最小模式为向上递增搜索切换到满足控制要求最靠近的模式。
经由上述的技术方案可知,与现有技术相比,本发明公开提供了一种基于模型预测控制的机械臂角速度控制方法;通过当前***状态预测机械臂的未来状态,将增益规划模型预测控制器(MPC)与直驱泵控***相结合,应用于不同液压机械臂,根据不同速度、负载工况下,调节压力波动,进而在机械臂工作中充当振动的阻尼,执行器压力和动臂角度位置通过在适当位置安装传感器进行监控,并用作反馈信号;将非线性模型在不同操作点处线性化,获得一个分段线性模型,以表示整个操作空间内的***;通过分段线性化建立不同增益规划模型预测控制器(MPC),设计切换逻辑,进而在不同运行速度和负载条件下的线性化阶段选择合适的增益规划模型预测控制器(MPC)来降低压力超调量,缩短到达稳定的时间,从而极大改善机械臂的动态性能。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图获得其他的附图。
图1附图为本发明提供的控制方法流程结构示意图。
图2附图为本发明提供的***原理框架结构示意图。
图3附图为本发明提供的非线性模型在操作空间上的响应面示意图。
图4附图为本发明提供的非线性模型的近似分段线性划分示意图。
图5附图为本发明提供的控制过程逻辑示意图。
图6附图为本发明提供的模型预测控制方法示意图。
其中,1为伺服电机,2、3为定量泵,4为蓄能器,5为液压缸,6、9为溢流阀,7、8为液控单向阀。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明实施例公开了一种基于模型预测控制的机械臂角速度控制方法,包括以下步骤:
S1、获取***的指令信号ucmd(t);
S2、对***进行分段线性化近似划分,构造不同线性化阶段的模型预测控制器;
S3、依据预测的下一阶段的指令信号ucmd(k+1)-uout(k+1)对模型预测控制器进行逻辑切换;
S4、获得压力反馈信号,压力反馈信号经模型预测控制器优化得到输出信号uout(t),将输出信号uout(t)作用于指令信号ucmd(t),得到优化后的输入信号u(t)=ucmd(t)-uout(t);
S5、基于优化后的输入信号u(t)对机械臂进行控制。
为进一步优化上述技术方案,步骤S2具体包括:
S21、构建状态模型,分段线性化***方程的状态空间形式为:
Xi(k+1)=AiXi(k)+Biu(k);i=(1,2,3…n)
y(k)=CXi(k)
其中,Xi(k+1)为预测***k+1阶段状态向量,Xi(k)是***k阶段状态向量,Ai是***矩阵,Bi是控制输入矩阵,u(k)是输入变量矩阵,y(k)为***输出,C为二阶单位矩阵;
***阶段状态向量表示为:
其中,θ为机械臂运行角度;为机械臂运行角度变化率;pA,pB分别为液压缸有杆腔和无杆腔反馈压力值;p′A,p′B分别为液压缸有杆腔和无杆腔忽略单向阀、管路等压降的两腔压力,理想情况下p'A=pA,p'B=pB;xV为液压缸运行速度;为液压缸运行变化率即加速度;βe为油液有效体积弹性模量;
在不同的线性化阶段考虑不同状态变量的影响,因此根据具体的控制信号大小,对***的状态向量进行划分,根据上述分段线性化***方程,得到增量式方程:
Δx(k+1)=Xi(k+1)-Xi(k)=AiΔX(k)+BiΔu(k)
其中,Δx(k+1)=x(k+1)-x(k),Δu(k+1)=u(k+1)-u(k)
同理,输出的增量方程为:
y(k)=CΔx(k)+y(k-1)
其中,y(k)为***当前阶段输出量,C为二阶单位矩阵,Δx(k)为***状态向量变化量;y(k-1)为上一阶段***输出量;
以最新的测量值为初始条件,预测步长为l,控制步长为m,并假设外部扰动对***影响很小,可得:
Δu(k+i)=0,i=m,m+1,…,k-1;
S22、构建状态预测模型:
其中,Y0为模型初始输入量;Su为***状态预测空间方程,C为二阶单位矩阵;A为***输入矩阵;Bi为控制输入矩阵;Δu(k)为输入控制增量;l为预测步长,m为控制步长
S23、在目标函数建立中,设置针对状态量代价Q矩阵和控制量的代价矩阵R权重,以对控制量、状态量的控制,使***的功能靠近理想目标,通过对目标函数进行求导,得出对于控制量输入的最优解:
o(k+i)=αiyk+(1-αi)yr
其中,o(k+i)为控制量输入最优解,即期望轨迹,αi为期望轨迹的权重因子,其值越小,响应速度越迅速,0<α<1,yr为目标值,yk为k时刻下***输出量;yr为***目标值;
模型预测控制器建立的目标函数为:
其中,rk+i是参考输出,yk+i是预测输出;Δu是控制信号与其标称值的偏差;Q为状态量代价矩阵;R为控制量代价矩阵;
则成本函数为:
其中Q权重大小与控制量有关,R权重大小与使用多少能量达到目标值有关,Δzk+i***输出跟踪与输入之间误差;Δuk+i-1为控制变化量;
控制目标目的在于使压力振荡最小化,同时输出信号与控制输入信号偏离达到最小,若目标函数最小,则其必要条件为:
通过求导,可得到最优解Δu(k),将最优解作为控制增量,得到下一时刻的信号输入;
则补偿函数为:
其中,H为补偿误差权重值,e为误差值;
Ycor=[Ycor(k+1)Ycor(k+2)…Ycor(k+p)]T;H=[h1h2…hp]T
将求得的补偿函数Ycor代入状态预测模型,即可得到新的加入补偿过后的预测函数:
为进一步优化上述技术方案,步骤S3具体包括:
定义切换函数s(k),s(k)在闭合区间内取整数值,其表示在第k个阶段处于活动状态的模型预测控制器的控制模式将成本函数j(k)作为预测模型控制器的候选李雅普诺夫函数,只有当候选李雅普诺夫函数在当前模式域外的其他模式下具有较低的值时,切换函数s(k)才改变其值,则模式切换判断方式为:
(1)ucmd(k+1)-uout(k+1)在当前模式的限制范围内,则控制器保持在当前模式,即s(k+1)=s(k);
(2)ucmd(k+1)-uout(k+1)超过当前模式的上限,则切换块在当前模式以上的模式中搜索j(k+1)最小的模式,并选择j(k+1)最小的模式;
(3)ucmd(k+1)-uout(k+1)低于当前模式的下限,则切换块在当前模式以下的模式中搜索j(k+1)最小的模式,并选择j(k+1)最小的模式。
优选的,限制范围为根据非线性模型在操作空间得到的响应面,划分不同阶段有着对应的控制命令区间范围;所以最小模式为向上递增搜索切换到满足控制要求最靠近的模式。
具体的,如2所示,直驱泵控液压***主要包括差动液压缸、一对为低压侧补油或排油的液控单向阀、双向定量泵/马达、用于代替低压油箱的液压蓄能器以及伺服电动机/发电机。
直驱泵控液压***是将两台与液压缸两侧面积成比例的定量泵A、B(2,3)和伺服电机(1)同轴连接,两泵出口分别与液压缸两侧相连。泵的进出口方向相反,当伺服电机给定正向旋转时,定量泵A(2)将液压油输送到液压缸A腔,B腔液压油从泵B(3)排出,通过此种方式解决了缸的流量不匹配问题。同时为了保证***安全运行,增加单向阀(7,8)与安全阀(6,9),根据***设定最大安全压力值。蓄能器(4)充当低压油箱。
通过压力传感器(10,11)分别监测液压缸A、B两腔的动态压力;通过角度位置传感器(12)测量动臂运行角度(θ)。
根据各部件数学公式原理建立***非线性数学模型。***的输入命令是命令信号(ucmd),***的输出是动臂角速度(ω)。控制命令信号(ucmd)在[-1,1]中间隔取值,动臂角度(θ)从26°变化到90°。在已经建立好的非线性数学模型运行模拟,并记录稳态角速度。针对整个操作空间生成非线性响应曲面。
根据得到的非线性响应曲面对***围绕多个工作点进行线性化。基于响应面,建立***的分段线性模型表示非线性模型,如图4所示。
分段线性化***方程可以用状态空间形式表示为:
根据划分的7个阶段,设计不同的模型预测控制器(MPC1,MPC2,…,MPC7),根据算法对未来的一组控制输入,基于二次成本函数最小化选择最优输入。
如图6所示,模型预测控制器的预测步长(l)是预测***未来状态的时间步长。控制器为其生成未来控制运动序列的时间步长称为控制步长(m)。预测步长(l)和控制步长(m)(l≥m)的值是根据开环装置响应的动力学和过程的时间常数来选择的。较长的预测步长(l)将增加计算成本,而较短的预测步长(l)将使控制达不到理想状况,同时存在不稳定现象。预测步长(l)的选择是一个在计算能力和控制器精度之间找到平衡的问题。
通过状态模型建立、预测模型建立、滚动优化、误差补偿四个步骤,根据调整预测步长(l)、控制步长(m)、状态量代价矩阵(Q)及控制量代价矩阵(R)设计不同的增益效果模型预测控制器。
将控制算法在电脑中编译完成,将不同增益效果模型预测控制器的采样时间(s)、预测范围(s)、控制范围(s)设置完成。并与建立好的液压***数学模型进行匹配。
在使用动态压力反馈的主动振荡中,来自执行机构的压力信号被反馈到控制器以产生控制信号。控制信号产生的流量将补偿执行机构中的压力波动,但必须首先过滤压力信号,以便在反馈至控制器之前,将瞬态部分与稳态值隔离。高通滤波器用于将压力信号的瞬态部分与稳态值隔离。
如图5所示,当液压缸压力(p)经过高通滤波器将压力信号(pHPF)反馈到合适的增益规划模型预测控制器,经过增益规划模型预测控制器优化得到输出命令uout(t),进而得到优化后输入命令u(t)=ucmd(t)-uout(t)。同步伺服电机控制器(i)对伺服电机转速(n)进行控制,最终对液压泵流量(Q)进行调节,对液压缸在不同负载、速度下的四象限工况引起的压力振荡进行阻尼抑制调节,极大减少机械臂振荡,提升***的动态性能。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的装置而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
Claims (4)
1.一种基于模型预测控制的机械臂角速度控制方法,其特征在于,包括以下步骤:
S1、获取***的指令信号ucmd(t);
S2、对***进行分段线性化近似划分,构造不同线性化阶段的模型预测控制器;
S3、依据预测的下一阶段的指令信号ucmd(k+1)-uout(k+1)对模型预测控制器进行逻辑切换;
S4、获得压力反馈信号,所述压力反馈信号经模型预测控制器优化得到输出信号uout(t),将输出信号uout(t)作用于指令信号ucmd(t),得到优化后的输入信号u(t)=ucmd(t)-uout(t);
S5、基于优化后的输入信号u(t)对机械臂进行控制。
2.根据权利要求1所述的一种基于模型预测控制的机械臂角速度控制方法,其特征在于,所述步骤S2具体包括:
S21、构建状态模型,分段线性化***方程的状态空间形式为:
Xi(k+1)=AiXi(k)+Biu(k);i=(1,2,3...n)
y(k)=CXi(k)
其中,Xi(k+1)为预测***k+1阶段状态向量,Xi(k)是***k阶段状态向量,Ai是***矩阵,Bi是控制输入矩阵,u(k)是输入变量矩阵,y(k)为***输出,C为二阶单位矩阵;
***阶段状态向量表示为:
其中,θ为机械臂运行角度;为机械臂运行角度变化率;pA,pB分别为液压缸有杆腔和无杆腔反馈压力值;分别为液压缸有杆腔和无杆腔忽略单向阀、管路等压降的两腔压力,理想情况下p′A=pA,p′B=pB;xV为液压缸运行速度;为液压缸运行变化率即加速度;βe为油液有效体积弹性模量;
在不同的线性化阶段考虑不同状态变量的影响,因此根据具体的控制信号大小,对***的状态向量进行划分,根据上述分段线性化***方程,得到增量式方程:
Δx(k+1)=Xi(k+1)-Xi(k)=AiΔX(k)+BiΔu(k)
其中,Δx(k+1)=x(k+1)-x(k),Δu(k+1)=u(k+1)-u(k)
同理,输出的增量方程为:
y(k)=CΔx(k)+y(k-1)
其中,y(k)为***当前阶段输出量,C为二阶单位矩阵,Δx(k)为***状态向量变化量;y(k-1)为上一阶段***输出量;
以最新的测量值为初始条件,预测步长为l,控制步长为m,并假设外部扰动对***影响很小,可得:
Δu(k+i)=0,i=m,m+1,...,l—1;
S22、构建状态预测模型:
其中,Y0为模型初始输入量;Su为***状态预测空间方程,C为二阶单位矩阵;A为***输入矩阵;Bi为控制输入矩阵;Δu(k)为输入控制增量;l为预测步长,m为控制步长
S23、在目标函数建立中,设置针对状态量代价Q矩阵和控制量的代价矩阵R权重,以对控制量、状态量的控制,使***的功能靠近理想目标,通过对目标函数进行求导,得出对于控制量输入的最优解:
o(k+i)=αiyk+(1-αi)yr
其中,o(k+i)为控制量输入最优解,即期望轨迹,αi为期望轨迹的权重因子,其值越小,响应速度越迅速,0<α<1,yr为目标值,yk为k时刻下***输出量;yr为***目标值;
模型预测控制器建立的目标函数为:
其中,rk+i是参考输出,yk+i是预测输出;Δu是控制信号与其标称值的偏差;Q为状态量代价矩阵;R为控制量代价矩阵;
则成本函数为:
其中Q权重大小与控制量有关,R权重大小与使用多少能量达到目标值有关,Δzk+i***输出跟踪与输入之间误差;Δuk+i-1为控制变化量;
控制目标目的在于使压力振荡最小化,同时输出信号与控制输入信号偏离达到最小,若目标函数最小,则其必要条件为:
通过求导,可得到最优解Δu(k),将最优解作为控制增量,得到下一时刻的信号输入;
则补偿函数为:
其中,H为补偿误差权重值,e为误差值;
Ycor=[Ycor(k+1)Ycor(k+2)...Ycor(k+p)]T;H=[h1h2...hp]T
将求得的补偿函数Ycor代入状态预测模型,即可得到新的加入补偿过后的预测函数:
3.根据权利要求1所述的一种基于模型预测控制的机械臂角速度控制方法,其特征在于,所述步骤S3具体包括:
定义切换函数s(k),s(k)在闭合区间内取整数值,其表示在第k个阶段处于活动状态的模型预测控制器的控制模式将成本函数j(k)作为预测模型控制器的候选李雅普诺夫函数,只有当候选李雅普诺夫函数在当前模式域外的其他模式下具有较低的值时,切换函数s(k)才改变其值,则模式切换判断方式为:
(1)ucmd(k+1)-uout(k+1)在当前模式的限制范围内,则控制器保持在当前模式,即s(k+1)=s(k);
(2)ucmd(k+1)-uout(k+1)超过当前模式的上限,则切换块在当前模式以上的模式中搜索j(k+1)最小的模式,并选择j(k+1)最小的模式;
(3)ucmd(k+1)-uout(k+1)低于当前模式的下限,则切换块在当前模式以下的模式中搜索j(k+1)最小的模式,并选择j(k+1)最小的模式。
4.根据权利要求3所述的一种基于模型预测控制的机械臂角速度控制方法,其特征在于,所述限制范围为根据非线性模型在操作空间得到的响应面,划分不同阶段有着对应的控制命令区间范围;所以最小模式为向上递增搜索切换到满足控制要求最靠近的模式。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210887459.2A CN115309044B (zh) | 2022-07-26 | 2022-07-26 | 一种基于模型预测控制的机械臂角速度控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210887459.2A CN115309044B (zh) | 2022-07-26 | 2022-07-26 | 一种基于模型预测控制的机械臂角速度控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115309044A true CN115309044A (zh) | 2022-11-08 |
CN115309044B CN115309044B (zh) | 2024-05-17 |
Family
ID=83859636
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210887459.2A Active CN115309044B (zh) | 2022-07-26 | 2022-07-26 | 一种基于模型预测控制的机械臂角速度控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115309044B (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080071397A1 (en) * | 2006-08-01 | 2008-03-20 | Rawlings James B | Partial Enumeration Model Predictive Controller |
US20130041482A1 (en) * | 2010-04-19 | 2013-02-14 | Abb Research Ltd | Method and system for updating a model in a model predictive controller |
CN106094524A (zh) * | 2016-07-07 | 2016-11-09 | 西北工业大学 | 基于输入趋势补偿的快速模型预测控制方法 |
CN108614435A (zh) * | 2018-08-08 | 2018-10-02 | 西安思优特自动控制技术有限公司 | 非线性无模型预测控制方法 |
JP2019113926A (ja) * | 2017-12-21 | 2019-07-11 | 株式会社Ihi | モデル予測制御装置 |
CN111555684A (zh) * | 2020-04-03 | 2020-08-18 | 浙江工业大学 | 一种可变开关点的无权重因子永磁同步电机有限集多步模型预测转矩控制方法 |
CN112213946A (zh) * | 2020-10-13 | 2021-01-12 | 辽宁石油化工大学 | 针对时变轨迹注塑成型异步切换过程的鲁棒预测控制方法 |
CN112428278A (zh) * | 2020-10-26 | 2021-03-02 | 北京理工大学 | 机械臂的控制方法、装置及人机协同模型的训练方法 |
CN113681543A (zh) * | 2021-08-25 | 2021-11-23 | 深圳亿嘉和科技研发有限公司 | 一种基于模型预测的机械臂零力控制方法 |
-
2022
- 2022-07-26 CN CN202210887459.2A patent/CN115309044B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080071397A1 (en) * | 2006-08-01 | 2008-03-20 | Rawlings James B | Partial Enumeration Model Predictive Controller |
US20130041482A1 (en) * | 2010-04-19 | 2013-02-14 | Abb Research Ltd | Method and system for updating a model in a model predictive controller |
CN106094524A (zh) * | 2016-07-07 | 2016-11-09 | 西北工业大学 | 基于输入趋势补偿的快速模型预测控制方法 |
JP2019113926A (ja) * | 2017-12-21 | 2019-07-11 | 株式会社Ihi | モデル予測制御装置 |
CN108614435A (zh) * | 2018-08-08 | 2018-10-02 | 西安思优特自动控制技术有限公司 | 非线性无模型预测控制方法 |
CN111555684A (zh) * | 2020-04-03 | 2020-08-18 | 浙江工业大学 | 一种可变开关点的无权重因子永磁同步电机有限集多步模型预测转矩控制方法 |
CN112213946A (zh) * | 2020-10-13 | 2021-01-12 | 辽宁石油化工大学 | 针对时变轨迹注塑成型异步切换过程的鲁棒预测控制方法 |
CN112428278A (zh) * | 2020-10-26 | 2021-03-02 | 北京理工大学 | 机械臂的控制方法、装置及人机协同模型的训练方法 |
CN113681543A (zh) * | 2021-08-25 | 2021-11-23 | 深圳亿嘉和科技研发有限公司 | 一种基于模型预测的机械臂零力控制方法 |
Non-Patent Citations (4)
Title |
---|
KILLPACK, MD等: "Fast Reaching in Clutter While Regulating Forces Using Model Predictive Control", 《13TH IEEE-RAS INTERNATIONAL CONFERENCE ON HUMANOID ROBOTS》, 31 December 2013 (2013-12-31) * |
刘国晴: "基于提高车辆稳定性的分布式电动汽车转矩协调控制研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》, no. 01, 15 January 2020 (2020-01-15) * |
孙浩然: "面向智能工厂的轮式移动机械臂运动控制研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》, no. 01, 15 January 2022 (2022-01-15) * |
熊清军: "基于模型预测控制的机械臂轨迹跟踪算法", 《工业技术创新》, vol. 8, no. 3, 31 December 2021 (2021-12-31) * |
Also Published As
Publication number | Publication date |
---|---|
CN115309044B (zh) | 2024-05-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Sarkar et al. | GA-optimized feedforward-PID tracking control for a rugged electrohydraulic system design | |
Grabbel et al. | An investigation of swash plate control concepts for displacement controlled actuators | |
JP5164047B2 (ja) | アクチュエータ、アクチュエータの制御方法及びアクチュエータの制御プログラム | |
Jing et al. | A novel architecture of electro-hydrostatic actuator with digital distribution | |
Williamson et al. | Pump mode prediction for four-quadrant velocity control of valueless hydraulic actuators | |
Ren et al. | Finite-time command filtered backstepping algorithm-based pitch angle tracking control for wind turbine hydraulic pitch systems | |
Wang et al. | Experimental investigation of an electro-hydrostatic actuator based on the novel active compensation method | |
CN111308889A (zh) | 一种喷杆***的自适应积分鲁棒控制方法 | |
CN115309044B (zh) | 一种基于模型预测控制的机械臂角速度控制方法 | |
CN108333919B (zh) | 一种考虑路面波动干扰的非平衡身管俯仰位置控制方法 | |
Zhou et al. | High precise fuzzy control for piezoelectric direct drive electro-hydraulic servo valve | |
CN110725817A (zh) | 一种基于比例溢流阀的进出口独立调节电液位置伺服***高能效控制方法 | |
CN115808876A (zh) | 发动机尾喷管执行机构自适应控制方法及装置 | |
CN113238476A (zh) | 一种非平衡炮管高精度抗扰控制方法 | |
CN114439805A (zh) | 调平***、调平方法及工程机械 | |
CN114439804A (zh) | 调平***、调平方法及工程机械 | |
Alaydi | Mathematical modeling for pump controlled system of hydraulic drive unit of single bucket excavator digging mechanism | |
CN112068447A (zh) | 大型装备电液***高精度鲁棒位置控制方法 | |
Yang et al. | Adaptive backstepping control of uncertain electro-hydrostatic actuator with unknown dead-zone nonlinearity | |
Pencelli et al. | Accurate position control for hydraulic servomechanisms | |
CN112631133B (zh) | 一种基于双蓄能器的液压位置伺服***控制方法 | |
CN112780637B (zh) | 一种举升液压伺服***节能及位置跟踪多目标控制方法 | |
Liu et al. | Flow match and precision motion control of asymmetric electro-hydrostatic actuators with complex external force in four-quadrants | |
CN117863491B (zh) | 一种高压树脂传递模塑工艺液压机滑块悬停控制方法 | |
CN114625006B (zh) | 一种高速大惯量电液位置伺服***的输出反馈控制方法 |
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 |