CN113485105B - 基于非对称输出的欧拉-伯努利梁自适应迭代控制方法 - Google Patents

基于非对称输出的欧拉-伯努利梁自适应迭代控制方法 Download PDF

Info

Publication number
CN113485105B
CN113485105B CN202110752838.6A CN202110752838A CN113485105B CN 113485105 B CN113485105 B CN 113485105B CN 202110752838 A CN202110752838 A CN 202110752838A CN 113485105 B CN113485105 B CN 113485105B
Authority
CN
China
Prior art keywords
euler
bernoulli beam
beam system
bernoulli
control
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.)
Active
Application number
CN202110752838.6A
Other languages
English (en)
Other versions
CN113485105A (zh
Inventor
刘屿
邬晓奇
万伟伟
横井浩史
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
South China University of Technology SCUT
Guangzhou Institute of Modern Industrial Technology
Original Assignee
South China University of Technology SCUT
Guangzhou Institute of Modern Industrial Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by South China University of Technology SCUT, Guangzhou Institute of Modern Industrial Technology filed Critical South China University of Technology SCUT
Priority to CN202110752838.6A priority Critical patent/CN113485105B/zh
Publication of CN113485105A publication Critical patent/CN113485105A/zh
Application granted granted Critical
Publication of CN113485105B publication Critical patent/CN113485105B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive 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/042Adaptive 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
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

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

本发明公开了一种基于非对称输出限制的欧拉‑伯努利梁边界自适应迭代控制方法,该方法过程如下:根据欧拉‑伯努利梁的动力学特征构建欧拉‑伯努利梁系统;根据欧拉‑伯努利梁系统,基于反步技术设计虚拟控制;制定辅助系统用于补偿欧拉‑伯努利梁输入饱和特性带来的影响;设计Lyapunov函数,设计迭代项;构建合适的Lyapunov函数,根据Lyapunov稳定性原理,设计基于非对称输出限制的自适应迭代控制方法。本发明能够有效抑制欧拉‑伯努利梁的振动偏移,并且在设计过程中考虑到了系统输入饱和特性和非对称输出限制。

Description

基于非对称输出的欧拉-伯努利梁自适应迭代控制方法
技术领域
本发明涉及振动控制技术领域,具体涉及一种基于非对称输出限制的欧拉-伯努利梁边界自适应迭代控制方法。
背景技术
凭借重量轻,效率高和能耗低等优良特点,柔性结构设备广泛被应用于机械臂、海洋立管和航天器等工程领域。在这些领域的研究中,欧拉-伯努利梁是这些工程领域设备的基础模型。在外部扰动的作用下,欧拉-伯努利梁会不断产生弹性形变,进而导致系统的长时间振动,这会影响系统的正常工作,降低系统的工作效率甚至会导致系统瘫痪。因此,如何有效地降低欧拉-伯努利梁的弹性形变和振动,是一个亟待解决的问题。
在现有的研究下,边界控制方法是一种有效能抑制欧拉-伯努利梁振动的控制方法;但在设计的过程中,很少考虑到欧拉-伯努利梁系统的输入饱和特性和非对称输出特性,这些特性在实际中无处不在,忽视这些特性,欧拉-伯努利梁很容易出现不稳定现象。
发明内容
本发明的目的是为了解决现有技术中的上述缺陷,提供一种基于非对称输出限制的欧拉-伯努利梁边界自适应迭代控制方法。
本发明的目的可以通过采取如下技术方案达到:
一种基于非对称输出限制的欧拉-伯努利梁边界自适应迭代控制方法,所述控制方法包括以下步骤:
根据欧拉-伯努利梁的动力学特征,构建欧拉-伯努利梁系统的动力学模型;
将所述欧拉-伯努利梁系统的动力学模型转为下三角形式,基于反步技术设计虚拟控制;
基于欧拉-伯努利梁系统输入饱和特性的影响,提出辅助系统用于补偿输入饱和该欧拉-伯努利梁带来的影响;
基于欧拉-伯努利梁系统受到外部周期性扰动d(t)的影响,构建欧拉- 伯努利梁系统的边界控制的迭代项:构建Lyapunov函数,并对该Lyapunov函数求导,得到迭代项的控制方法;
基于反步技术选取Lyapunov函数,得到边界自适应迭代控制方法,其中,所述边界控制方法包括参数自适应控制律,用于补偿参数不确定性。
进一步地,所述动力学特征包括欧拉-伯努利梁系统的动能、势能以及非保守力对欧拉-伯努利梁系统所做的虚功,将动能、势能、虚功代入哈密顿原理,得到欧拉-伯努利梁系统的动力学模型为:
Figure BDA0003145628590000021
其中s(x,t)为欧拉-伯努利梁系统在长度为x、时间t时产生的偏移量,
Figure BDA0003145628590000022
表示s(x,t)对时间t的二次导数,s″(x,t)和s″″(x,t)分别表示s(x,t)对 x的二阶导数和四阶导数,
其边界条件为:
Figure BDA0003145628590000023
Figure BDA0003145628590000024
其中,l为欧拉-伯努利梁系统的长度,ρ为欧拉-伯努利梁系统的密度,EI为欧拉-伯努利梁系统的弯曲刚度,T为欧拉-伯努利梁系统的张力,Ms为欧拉-伯努利梁系统末端负载的质量,s(0,t)为欧拉- 伯努利梁在长度为0时间为t时的偏移量,s′(0,t)为s(0,t)关于t的变化速率,s′(l,t)为s(x,t)对x一阶导数s′(x,t)在x为l的值,s″(l,t)为s(x,t)对x 二阶导数s″(x,t)在x为l的值,s″′(l,t)为s(x,t)对x三阶导数s″′(x,t)在x为l 的值,
Figure BDA0003145628590000031
表示欧拉-伯努利梁系统在长度为l处位移的加速度,u(t)为控制输入,d(t)为欧拉-伯努利梁系统边界处的周期性扰动。
进一步地,将所述欧拉-伯努利梁系统的动力学模型转为下三角形式,基于反步技术设计虚拟控制的过程如下:
Figure BDA0003145628590000032
其中x1(t)为欧拉-伯努利梁系统在l处的振动偏移量,
Figure BDA0003145628590000033
和x2(t)都为欧拉-伯努利梁系统在l处位移量的变化速率,
Figure BDA0003145628590000034
描述欧拉-伯努利梁系统在l处速率变化的快慢;
定义v(t)为x2(t)的虚拟控制,具体为v(t)=-λ1s′(l,t)+λ2s″′(l,t),其中λ12为控制参数,且λ12>0;
定义e(t)为x2(t)与虚拟控制v(t)之间的误差,具体为:e(t)=x2(t)-v(t)。
进一步地,提出辅助系统用于补偿输入饱和特性给该欧拉-伯努利梁系统带来的影响,具体如下:
将欧拉-伯努利梁系统的饱和输入特性表示为:
Figure BDA0003145628590000035
其中u0(t)表示设计的欧拉-伯努利梁系统输入,um表示饱和输入的最大绝对值;
考虑到欧拉-伯努利梁系统非对称输出,提出如下Lyapunov函数:
Figure BDA0003145628590000036
c1≠c2,c1>0,c2<0且c2<s(l,t)<c1,c1为该欧拉-伯努利梁系统输出的上限值,c2为该欧拉-伯努利梁系统输出的下限值,
Figure BDA0003145628590000041
设计辅助系统用来补偿输入饱和特性,具体为
Figure BDA0003145628590000042
其中λ3>0,为辅助系统的控制参数;Δu(t)=u(t)-u0(t),是控制输入u(t) 与设计的控制输入u0(t)间的差值;μ(t)为辅助系统的状态,
Figure BDA0003145628590000043
表示辅助系统状态的变化速率,μ0>0,为辅助系统的状态临界值,限定函数
Figure BDA0003145628590000044
Figure BDA0003145628590000045
起到限定偏移量位于上限值c1和下限值c2区间之内的作用,如果x1(t)>0时,J(x1(t))的值为1,此时起作用的是限定函数
Figure BDA0003145628590000046
的第一项,使得欧拉-伯努利梁系统的偏移量的值不高于上限值c1;如果x1(t)<0时,J(x1(t))的值为0,此时起作用的是限定函数
Figure BDA0003145628590000047
的第二项,使得欧拉-伯努利梁系统的偏移量的值不低于下限值c2
进一步地,基于欧拉-伯努利梁系统受到外部周期性扰动d(t)的影响,构建欧拉-伯努利梁系统的边界控制的迭代项,该迭代项用于消除外部周期性扰动d(t)的影响,过程如下:
选取如下的Lyapunov函数,表达式为:
Figure BDA0003145628590000048
其中,
Figure BDA0003145628590000049
是一个大于0的估计值且满足
Figure BDA00031456285900000410
To是d(t)的周期,θ是控制参数且0<θ<1,φ(t)为迭代项;
对Vh(t)进行求导,根据Lyapunov稳定性原理,设计出迭代项,迭代项具体为
Figure BDA00031456285900000411
其中
Figure BDA0003145628590000051
是d1(t)关于时间的导数,
Figure BDA0003145628590000052
为迭代项φ(t)关于时间的导数,φ(t-To)为t-To时刻的迭代项的值。
进一步地,基于反步技术选取新的Lyapunov函数,得到边界控制方法,过程如下:
Figure BDA0003145628590000053
其中,V(t)为新构造的Lyapunov函数,从能量函数角度来看,V1(t)为能量项,V2(t)交叉项,Ve(t)为误差控制项,用于处理非对称输出限制,
Figure BDA0003145628590000054
为辅助系统附加项,具体表达式为:
Figure BDA0003145628590000055
Figure BDA0003145628590000056
Figure BDA0003145628590000057
其中,α为能量项V1(t)的加权常数,η为辅助项的加权常数。
进一步地,所述边界自适应迭代控制方法,基于选取的Lyapunov函数,对V(t)求导,根据Lyapunov稳定性原理,保证其负定性,具体如下:
Figure BDA0003145628590000058
其中,k1和k均为大于0的控制参数,φ(t)为迭代项,
Figure BDA0003145628590000059
Figure BDA00031456285900000510
Figure BDA00031456285900000511
是s′(l,t)关于时间t的导数,
Figure BDA00031456285900000512
是s″′(l,t)关于时间t的导数;
Figure BDA00031456285900000513
Figure BDA00031456285900000514
分别为EI、T、Ms的估计值,
Figure BDA00031456285900000515
分别为对应参数的误差值,定义P=[EI、T、Ms],相应的参数误差矩阵为P-P(t)=[EI T Ms]T
Figure BDA00031456285900000516
所述参数自适应控制律为
Figure BDA00031456285900000517
Figure BDA0003145628590000061
关于时间t的导数,RT(t)是R(t)矩阵的转置矩阵;Γ是一个属于实数集空间上的三行三列正定对角矩阵,且ξ>0,该参数自适应控制律用于补偿欧拉-伯努利梁系统的参数的不确定性。
本发明相对于现有技术具有如下的优点及效果:
本发明提出的基于非对称输出限制的欧拉-伯努利梁边界自适应迭代控制方法与传统的控制方法相比,该控制方法易于实现,控制精度高,适应性强。本设计方法中含有迭代项,可以利用先前相关信息来产生期望输出,改善控制质量。并且,可以同时解决欧拉-伯努利梁系统输入饱和现象和处理非对称输出特性。
附图说明
图1是本发明公开的一种基于非对称输出限制的欧拉-伯努利梁边界自适应迭代控制方法的流程示意图;
图2是本发明公开的欧拉-伯努利梁系统的结构示意图;
图3是本发明实施例中未施加控制的欧拉-伯努利梁的弹性形变s(x, t)仿真结果示意图;
图4是本发明实施例中施加控制后,迭代次数k=9的欧拉-伯努利梁的弹性形变s(x,t)的仿真结果示意图;
图5是本发明实施例中施加控制后,迭代次数k=9的欧拉-伯努利梁的末端弹性形变s(l,t)的仿真结果示意图;
图6是本发明实施例中的施加控制后,欧拉-伯努利梁的边界弹性形变 s(l,t)的最大数值与迭代次数k的关系示意图;
图7是本发明提供的控制设计命令的输入波形的分析示意图;
图8是本发明提供的具有饱和非线性特征的控制输入分析示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例一
参阅图1,本实施例提供了一种基于非对称输出限制的欧拉-伯努利梁边界自适应迭代控制方法,针对该方法,至少包含以下步骤:
S1、根据欧拉-伯努利梁的动力学特征,构建欧拉-伯努利梁系统的动力学模型。
如图2所示,一种典型的欧拉-伯努利梁系统,梁的左侧边界固定于坐标原点,边界控制器u(t)和外部周期扰动d(t)作用于柔性梁的右侧。柔性梁长度为l,其振动偏移量为s(l,t)。
欧拉-伯努利梁式结构的动力学方程为:
Figure BDA0003145628590000071
其中,s(x,t)为欧拉-伯努利梁系统在长度为x,时间t时产生的偏移量;
Figure BDA0003145628590000072
表示为s(x,t)对时间t的二次导数,s″′(x,t)和s″″(x,t)分别表示s(x,t)对x的三阶导数和四阶导数。
其边界条件为:
Figure BDA0003145628590000073
Figure BDA0003145628590000074
其中,l为欧拉-伯努利梁的长度,ρ为该系统的密度, EI为所述系统的弯曲刚度,T为欧拉-伯努利梁的张力,Ms为系统末端负载的质量。
S2、基于反步法,将欧拉-伯努利梁的动力学方程转为下三角形式,并构建虚拟控制。
在本实施例中,下三角形式为:
Figure BDA0003145628590000081
定义v(t)为x2(t)的虚拟控制,具体为v(t)=-λ1s′(l,t)+λ2s″′(l,t),其中λ12为控制参数,且λ12>0。
x2(t)与虚拟控制v(t)之间的误差为:e(t)=x2(t)-v(t)。
S3、设计辅助系统补偿欧拉-伯努利梁输入饱和带来的影响。大多数边界控制方法没考虑非线性,因此本发明通过设置辅助系统来补偿欧拉-伯努利梁输入饱和特性。
考虑到欧拉-伯努利梁系统非对称输出,提出如下Lyapunov函数:
Figure BDA0003145628590000082
c1≠c2,c1>0,c2<0且c2<s(l,t)<c1,c1为该欧拉-伯努利梁系统输出的上限值,c2为该欧拉-伯努利梁系统输出的下限值,
Figure BDA0003145628590000083
系统的饱和输入特性表示为:
Figure BDA0003145628590000084
其中u0(t)表示设计的系统输入,um表示饱和输入的最大绝对值。
设计辅助系统用来补偿输入饱和特性,具体为
Figure BDA0003145628590000085
其中λ3>0,Δu(t)=u(t)-u0(t),μ(t)为辅助系统的状态函数,μ0为设计参数,
Figure BDA0003145628590000091
S4、构建Lyapunov函数,并对该函数求导,得到迭代部分的控制方法来处理外部的周期性扰动d(t)。
选取如下的Lyapunov函数,具体为:
Figure BDA0003145628590000092
其中,
Figure BDA0003145628590000093
是一个估计值且满足
Figure BDA0003145628590000094
To是d(t)的周期,θ是控制参数且0<θ<1;φ(t)为迭代项。
对Vh(t)进行求导,根据Lyapunov稳定性原理,即保证
Figure BDA0003145628590000095
负定性,设计出迭代项,迭代项具体为
Figure BDA0003145628590000096
S5、基于反步法,选取Lyapunov函数,对V(t)求导,根据Lyapunov 稳定性原理,保证其负定性,设计边界控制方法,具体为:
Figure BDA0003145628590000097
其中,V1(t)为能量项,V2(t)交叉项,Ve(t)为误差控制项,用来处理非对称输出特性,
Figure BDA0003145628590000098
为辅助系统附加项;具体表达式为:
Figure BDA0003145628590000099
Figure BDA00031456285900000910
Figure BDA00031456285900000911
对选取的Lyapunov函数,对V(t)求导,根据Lyapunov稳定性原理,保证其负定性,设计边界控制方法,具体如下:
Figure BDA00031456285900000912
其中,φ(t)为迭代项,用来处理外部周期性扰动,
Figure BDA0003145628590000101
Figure BDA0003145628590000102
P(t)=[EI T Ms]T,EI,T,Ms分别为EI,T,Ms的估计值,相应的参数误差矩阵为P-P(t)=[EI T Ms]T
Figure BDA0003145628590000103
参数自适应控制律为
Figure BDA0003145628590000104
Γ是一个属于实数集空间上的三行三列正定对角矩阵,ξ>0;将参数自适应控制律用于补偿该系统的参数的不确定性。
边界控制器中的所有信号均可由传感器直接测得或由有限差分法间接获得。
实施例二
本实施例基于Lyapunov稳定性原理,验证在自适应迭代控制方法的作用下,欧拉-伯努利梁系统的稳定性。
Lyapunov函数更新为
Figure BDA0003145628590000105
其中Vp(t)为
Figure BDA0003145628590000106
在本实施例中,验证V(t)的正定性,方法如下:
根据不等式的放缩原理,对于V2(t)有
Figure BDA0003145628590000107
其中
Figure BDA0003145628590000108
-a1V1(t)≤V2(t)≤a1V1(t),为了保证0≤a1≤1,需满足
Figure BDA0003145628590000109
定义
Figure BDA0003145628590000111
记Va(t)=V1(t)+V2(t),
Figure BDA0003145628590000112
则有Va(t)≤(1+a1)V1(t)≤a3V1(t),则V(t)可写成为:
a2V1(t)+Ve(t)+Vb(t)≤V(t)≤a3V1(t)+Ve(t)+Vb(t),其中V1(t)≥0,Ve(t)≥0, Vb(t)≥0,所以可知V(t)≥0,V(t)的正定性得到证明。
验证V(t)的一阶导数
Figure BDA0003145628590000113
的负定性方法如下:
V1(t),V2(t)对时间求导可得:
Figure BDA0003145628590000114
Ve(t)对时间求导可得:
Figure BDA0003145628590000115
Vb(t)对时间求导,
Figure BDA0003145628590000116
其中,λb是矩阵Γ的最小特征值,将迭代项
Figure BDA0003145628590000117
代入(3)式可得
Figure BDA0003145628590000121
Figure BDA0003145628590000122
将式(1)、式(2)、式(3)、式(4)代入式(5),并结合设计的边界控制方法u0(t),化简可得:
Figure 1
其中δ1~δ4均为根据实际情况选取的控制参数,是大于0的常数,且满足
Figure BDA0003145628590000124
Figure BDA0003145628590000125
其中,
Figure BDA0003145628590000126
λ3=min(1-θ,θ)。
以上证明了
Figure BDA0003145628590000127
的负定性,系统在设计的控制方法下是渐近稳定的。
分析系统状态的有界性:
将式(6)乘以e-ζt并且积分可得:
Figure BDA0003145628590000131
结合V1(t)和上式可得:
Figure BDA0003145628590000132
进而有
Figure BDA0003145628590000133
当判断欧拉-伯努利梁系统满足预设的稳定性要求时,利用Matlab仿真软件对该欧拉-伯努利梁系统进行数字仿真,得到仿真结果来验证自适应迭代控制方法的有效性。
在本实施例中,请参阅图3,图3为本发明实施例中未施加控制的欧拉-伯努利梁系统的弹性形变s(x,t)仿真结果示意图。在未施加控制时,欧拉-伯努利梁各处均存在振动(横向位移)。图4为本发明实施例中施加控制作用后,迭代次数为9后的欧拉-伯努利梁系统的末端弹性形变仿真结果示意图,当施加控制后,其末端位置的横向位移的最大值不超过0.012,欧拉-伯努利梁的弹性形变趋于相对平稳。图5是本发明实施例中施加控制作用后,迭代次数为9后的欧拉-伯努利梁系统的弹性形变数值示意图,当施加控制后,欧拉-伯努利梁系统的弹性形变的数值始终维持在-0.003到0.011 之间,欧拉-伯努利梁系统的非对称输出限制符合要求。图6是本发明实施例中施加控制作用后,欧拉-伯努利梁系统在每次迭代过程中的末端振动量最大值随次数的变化情况,可以看出,随着迭代次数的增加,欧拉-伯努利梁系统的振动幅度逐渐减小,迭代项产生了控制作用;从图7和图8可以看出控制设计命令的输入波形在饱和非线性特征作用下依然能将输入信号控制在有效输入范围内,非线性输入饱和特性在辅助系统作用下得到妥善处理。
上述实施例为本发明较佳的实施方式,但本发明的实施方式并不受上述实施例的限制,其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化,均应为等效的置换方式,都包含在本发明的保护范围之内。

Claims (4)

1.一种基于非对称输出限制的欧拉-伯努利梁边界自适应迭代控制方法,其特征在于,所述控制方法包括以下步骤:
根据欧拉-伯努利梁的动力学特征,构建欧拉-伯努利梁系统的动力学模型;所述动力学特征包括欧拉-伯努利梁系统的动能、势能以及非保守力对欧拉-伯努利梁系统所做的虚功,将动能、势能、虚功代入哈密顿原理,得到欧拉-伯努利梁系统的动力学模型为:
Figure FDA0003955939490000011
其中s(x,t)为欧拉-伯努利梁系统在长度为x、时间t时产生的偏移量,
Figure FDA0003955939490000012
表示s(x,t)对时间t的二次导数,s″(x,t)和s″″(x,t)分别表示s(x,t)对x的二阶导数和四阶导数,
其边界条件为:
Figure FDA0003955939490000013
Figure FDA0003955939490000014
其中,l为欧拉-伯努利梁系统的长度,ρ为欧拉-伯努利梁系统的密度,EI为欧拉-伯努利梁系统的弯曲刚度,T为欧拉-伯努利梁系统的张力,Ms为欧拉-伯努利梁系统末端负载的质量,s(0,t)为欧拉-伯努利梁在长度为0时间为t时的偏移量,s′(0,t)为s(0,t)关于t的变化速率,s′(l,t)为s(x,t)对x一阶导数s′(x,t)在x为l的值,s″(l,t)为s(x,t)对x二阶导数s″(x,t)在x为l的值,s″′(l,t)为s(x,t)对x三阶导数s″′(x,t)在x为l的值,
Figure FDA0003955939490000015
表示欧拉-伯努利梁系统在长度为l处位移的加速度,u(t)为控制输入,d(t)为欧拉-伯努利梁系统边界处的周期性扰动;
将所述欧拉-伯努利梁系统的动力学模型转为下三角形式,基于反步技术设计虚拟控制,过程如下:
Figure FDA0003955939490000021
其中x1(t)为欧拉-伯努利梁系统在l处的振动偏移量,
Figure FDA0003955939490000022
和x2(t)都为欧拉-伯努利梁系统在l处位移量的变化速率,
Figure FDA0003955939490000023
描述欧拉-伯努利梁系统在l处速率变化的快慢;
定义v(t)为x2(t)的虚拟控制,具体为v(t)=-λ1s′(l,t)+λ2s″′(l,t),其中λ12为控制参数,且λ12>0;
定义e(t)为x2(t)与虚拟控制v(t)之间的误差,具体为:e(t)=x2(t)-v(t);
基于欧拉-伯努利梁系统输入饱和特性的影响,提出辅助系统用于补偿输入饱和该欧拉-伯努利梁带来的影响,具体如下:
将欧拉-伯努利梁系统的饱和输入特性表示为:
Figure FDA0003955939490000024
其中u0(t)表示设计的欧拉-伯努利梁系统输入,um表示饱和输入的最大绝对值;
考虑到欧拉-伯努利梁系统非对称输出,提出如下Lyapunov函数:
Figure FDA0003955939490000025
c1≠c2,c1>0,c2<0且c2<s(l,t)<c1,c1为该欧拉-伯努利梁系统输出的上限值,c2为该欧拉-伯努利梁系统输出的下限值,
Figure FDA0003955939490000026
设计辅助系统用来补偿输入饱和特性,具体为
Figure FDA0003955939490000027
其中λ3>0,为辅助系统的控制参数;Δu(t)=u(t)-u0(t),是控制输入u(t)与设计的控制输入u0(t)间的差值;μ(t)为辅助系统的状态,
Figure FDA0003955939490000031
表示辅助系统状态的变化速率,μ0>0,为辅助系统的状态临界值,限定函数
Figure FDA0003955939490000032
起到限定偏移量位于上限值c1和下限值c2区间之内的作用,如果x1(t)>0时,J(x1(t))的值为1,此时起作用的是限定函数
Figure FDA0003955939490000033
的第一项,使得欧拉-伯努利梁系统的偏移量的值不高于上限值c1;如果x1(t)<0时,J(x1(t))的值为0,此时起作用的是限定函数
Figure FDA0003955939490000034
的第二项,使得欧拉-伯努利梁系统的偏移量的值不低于下限值c2
基于欧拉-伯努利梁系统受到外部周期性扰动d(t)的影响,构建欧拉-伯努利梁系统的边界控制的迭代项:构建Lyapunov函数,并对该Lyapunov函数求导,得到迭代项的控制方法;
基于反步技术选取Lyapunov函数,得到边界自适应迭代控制方法,其中,所述边界控制方法包括参数自适应控制律,用于补偿参数不确定性。
2.根据权利要求1所述的基于非对称输出限制的欧拉-伯努利梁边界自适应迭代控制方法,其特征在于,基于欧拉-伯努利梁系统受到外部周期性扰动d(t)的影响,构建欧拉-伯努利梁系统的边界控制的迭代项,该迭代项用于消除外部周期性扰动d(t)的影响,过程如下:
选取如下的Lyapunov函数,表达式为:
Figure FDA0003955939490000035
其中,
Figure FDA0003955939490000036
Figure FDA0003955939490000037
是一个大于0的估计值且满足
Figure FDA0003955939490000038
To是d(t)的周期,θ是控制参数且0<θ<1,φ(t)为迭代项;
对Vh(t)进行求导,根据Lyapunov稳定性原理,设计出迭代项,迭代项具体为
Figure FDA0003955939490000039
其中
Figure FDA0003955939490000041
是d1(t)关于时间的导数,
Figure FDA0003955939490000042
为迭代项φ(t)关于时间的导数,φ(t-To)为t-To时刻的迭代项的值。
3.根据权利要求2所述的基于非对称输出限制的欧拉-伯努利梁边界自适应迭代控制方法,其特征在于,基于反步技术选取新的Lyapunov函数,得到边界控制方法,过程如下:
Figure FDA0003955939490000043
其中,V(t)为新构造的Lyapunov函数,从能量函数角度来看,V1(t)为能量项,V2(t)交叉项,Ve(t)为误差控制项,用于处理非对称输出限制,
Figure FDA0003955939490000044
为辅助系统附加项,具体表达式为:
Figure FDA0003955939490000045
Figure FDA0003955939490000046
Figure FDA0003955939490000047
其中,α为能量项V1(t)的加权常数,η为辅助项的加权常数。
4.根据权利要求3所述的基于非对称输出限制的欧拉-伯努利梁边界自适应迭代控制方法,其特征在于,所述边界自适应迭代控制方法,基于选取的Lyapunov函数,对V(t)求导,根据Lyapunov稳定性原理,保证其负定性,具体如下:
Figure FDA0003955939490000048
其中,k1和k均为大于0的控制参数,φ(t)为迭代项,
Figure FDA0003955939490000049
Figure FDA00039559394900000410
Figure FDA00039559394900000411
是s′(l,t)关于时间t的导数,
Figure FDA00039559394900000412
是s″′(l,t)关于时间t的导数;
Figure FDA00039559394900000413
Figure FDA00039559394900000414
Figure FDA00039559394900000415
分别为EI、T、Ms的估计值,
Figure FDA00039559394900000416
分别为对应参数的误差值,定义P=[EI、T、Ms],相应的参数误差矩阵为P-P(t)=[EI T Ms]T
Figure FDA0003955939490000051
所述参数自适应控制律为
Figure FDA0003955939490000052
Figure FDA0003955939490000053
Figure FDA0003955939490000054
关于时间t的导数,RT(t)是R(t)矩阵的转置矩阵;Γ是一个属于实数集空间上的三行三列正定对角矩阵,且ξ>0,该参数自适应控制律用于补偿欧拉-伯努利梁系统的参数的不确定性。
CN202110752838.6A 2021-07-02 2021-07-02 基于非对称输出的欧拉-伯努利梁自适应迭代控制方法 Active CN113485105B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110752838.6A CN113485105B (zh) 2021-07-02 2021-07-02 基于非对称输出的欧拉-伯努利梁自适应迭代控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110752838.6A CN113485105B (zh) 2021-07-02 2021-07-02 基于非对称输出的欧拉-伯努利梁自适应迭代控制方法

Publications (2)

Publication Number Publication Date
CN113485105A CN113485105A (zh) 2021-10-08
CN113485105B true CN113485105B (zh) 2023-02-28

Family

ID=77940438

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110752838.6A Active CN113485105B (zh) 2021-07-02 2021-07-02 基于非对称输出的欧拉-伯努利梁自适应迭代控制方法

Country Status (1)

Country Link
CN (1) CN113485105B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114442483B (zh) * 2022-01-18 2023-10-20 曲阜师范大学 一种自适应非周期间歇边界控制方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU2018100710A4 (en) * 2018-03-28 2018-06-28 Guizhou University Adaptive Neural Network Control Method For Arch-Shaped Microelectromechanical System
CN110609471A (zh) * 2019-07-23 2019-12-24 华南理工大学 基于反步技术的海洋柔性立管系统的边界迭代控制方法
CN110673469A (zh) * 2019-07-23 2020-01-10 华南理工大学 一种基于反步迭代学习的欧拉-伯努利梁的振动控制方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU2018100710A4 (en) * 2018-03-28 2018-06-28 Guizhou University Adaptive Neural Network Control Method For Arch-Shaped Microelectromechanical System
CN110609471A (zh) * 2019-07-23 2019-12-24 华南理工大学 基于反步技术的海洋柔性立管系统的边界迭代控制方法
CN110673469A (zh) * 2019-07-23 2020-01-10 华南理工大学 一种基于反步迭代学习的欧拉-伯努利梁的振动控制方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
基于Lyapunov直接法的柔性梁振动控制;刘玙等;《华南理工大学学报(自然科学版)》;20130215;第41卷(第02期);第24-29页 *
基于内流动力学的海洋输油柔性立管鲁棒边界控制;高红霞等;《控制理论与应用》;20120615;第29卷(第06期);第785-791页 *
基于时滞的柔性机械臂轨迹追踪;任豪;《现代机械》;20200228(第01期);全文 *

Also Published As

Publication number Publication date
CN113485105A (zh) 2021-10-08

Similar Documents

Publication Publication Date Title
Meng et al. Iterative learning control of a robotic arm experiment platform with input constraint
CN110673469B (zh) 一种基于反步迭代学习的欧拉-伯努利梁的振动控制方法
Sheng et al. Observer-based adaptive backstepping control for fractional order systems with input saturation
US20220234198A1 (en) Method for controlling vibration of flexible mechanical arm based on cooperative tracking
CN110647104B (zh) 一种基于边界扰动观测器的柔性立管反步边界控制方法
CN113485105B (zh) 基于非对称输出的欧拉-伯努利梁自适应迭代控制方法
Zhao et al. Robust output feedback stabilization for a flexible marine riser system
CN108388136B (zh) 一种具有非线性输入的柔性弦线系统的振动控制方法
CN111399379B (zh) 一种基于输出限制的柔性卫星系统的自适应振动控制方法
CN112068444B (zh) 一种采用非线性自适应滑模的飞行器攻角控制方法
Zhai et al. Analysis and design of H∞ feedback control systems with two quantized signals
CN109995278B (zh) 一种考虑输入受限的电机伺服系统自调节控制方法
CN112859948B (zh) 基于分布扰动的柔性杆振动抑制方法及系统
Niknam et al. On the passive control of friction-induced instability due to mode coupling
CN115840362A (zh) 一种非线性奇异时滞系统的有限时间鲁棒镇定控制方法
Khan et al. Three control strategies for unified chaotic system
CN112859949B (zh) 基于Euler-Bernoulli梁的边界抗饱和振动抑制方法及系统
Oya et al. Robust non-fragile controllers for uncertain linear continuous-time systems
CN115576193A (zh) 基于delta算子描述的pid控制方法
CN113211448B (zh) 基于扰动观测器的柔性臂自适应边界振动控制方法及系统
Yang et al. Robust Reference Tracking of Linear Uncertain Systems via Uncertainty Estimation and Composite Control
Tian et al. Adaptive stabilization for fractional-order system with nonsymmetrical dead-zone input via backstepping-based sliding mode control
Zhao et al. Vibration control of an experimental flexible manipulator against input saturation
Okajima et al. Analysis of dynamic quantizer in 2-DOF internal model control system with dead-time
CN114660954B (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