CN108818541B - 一种柔性关节机器人的自适应神经网络跟踪控制方法 - Google Patents

一种柔性关节机器人的自适应神经网络跟踪控制方法 Download PDF

Info

Publication number
CN108818541B
CN108818541B CN201811057430.1A CN201811057430A CN108818541B CN 108818541 B CN108818541 B CN 108818541B CN 201811057430 A CN201811057430 A CN 201811057430A CN 108818541 B CN108818541 B CN 108818541B
Authority
CN
China
Prior art keywords
equation
error
control
formula
matrix
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
CN201811057430.1A
Other languages
English (en)
Other versions
CN108818541A (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.)
Minjiang University
Original Assignee
Minjiang University
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 Minjiang University filed Critical Minjiang University
Priority to CN201811057430.1A priority Critical patent/CN108818541B/zh
Publication of CN108818541A publication Critical patent/CN108818541A/zh
Application granted granted Critical
Publication of CN108818541B publication Critical patent/CN108818541B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B25HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
    • B25JMANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
    • B25J9/00Programme-controlled manipulators
    • B25J9/16Programme controls
    • B25J9/1628Programme controls characterised by the control loop
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B25HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
    • B25JMANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
    • B25J13/00Controls for manipulators

Landscapes

  • Engineering & Computer Science (AREA)
  • Robotics (AREA)
  • Mechanical Engineering (AREA)
  • Feedback Control In General (AREA)

Abstract

本发明涉及一种柔性关节机器人的自适应神经网络跟踪控制方法,包括超声波电机,所述超声波电机一侧输出轴与光电编码器相连接,另一侧输出轴与飞轮惯性负载相连接,所述飞轮惯性负载的输出轴经联轴器与力矩传感器相连接,所述光电编码器的信号输出端、所述力矩传感器的信号输出端分别接至控制***;所述控制***建立在自适应神经网络基础上,使用补偿器使得跟踪误差最小,从而能获得更好的控制效能。

Description

一种柔性关节机器人的自适应神经网络跟踪控制方法
技术领域
本发明涉及电机控制领域,特别是一种柔性关节机器人的自适应神经网络跟踪控制方法。
背景技术
许多机器人,关节的灵活性可以通过谐波传动轴,造成饱和,轴承变形,在液压机械的液压流体的可压缩性,等等。在机器人模型介绍及灵活性相当复杂的运动方程。特别是相关动力学的规则变为刚性机器人的两倍,自由度的数目大于控制输入的数量。因此很难控制这样的***。
本发明讨论了一类由直流电机驱动的柔性关节机器人的鲁棒跟踪控制问题。采用自适应神经网络***来逼近不确定的机电动力学行为。构造了一个降阶观测器来估计速度信号。反馈信号只需要测量链路位置和电枢电流。利用反推技术,建立了一种基于观测器的鲁棒无速度反馈动态跟踪控制器,使闭环***的所有状态和信号都有界,并可使轨迹跟踪误差尽可能小。
发明内容
有鉴于此,本发明的目的是提出一种柔性关节机器人的自适应神经网络跟踪控制方法,整个控制器的***建立在自适应神经网络基础上,使用补偿器使得跟踪误差最小,从而能获得更好的控制效能。
本发明采用以下方案实现:一种柔性关节机器人的自适应神经网络跟踪控制方法,包括超声波电机,所述超声波电机一侧输出轴与光电编码器相连接,另一侧输出轴与飞轮惯性负载相连接,所述飞轮惯性负载的输出轴经联轴器与力矩传感器相连接,所述光电编码器的信号输出端、所述力矩传感器的信号输出端分别接至控制***;所述控制***建立在自适应神经网络基础上,使用补偿器使得跟踪误差最小,从而能获得更好的控制效能,具体设计为:
所述超声波电机采用由直流伺服电机驱动的柔性关节机械臂运动方程:
Figure BDA0001796163990000021
Figure BDA0001796163990000022
Figure BDA0001796163990000023
式中,其中q1∈Rn、q2∈Rn分别表示连杆角度和电机角度,M(q1)∈Rn×n是对称正定惯性矩阵,
Figure BDA0001796163990000024
表示向心力和科里奥利力的矢量,G(q1)∈Rn表示引力矢量,K=diag{Ki}∈Rn×n是关节刚度系数矩阵,J∈Rn×n是电机惯量矩阵,B∈Rn×n是执行器阻尼矩阵,KT=diag{Kti}∈Rn×n是电机的转矩常数矩阵,I∈Rn是电枢电流矢量,L=diag{Li}∈Rn×n是电感矩阵,R=diag{Ri}∈Rn×n是电阻矩阵,KB=diag{Kbi}∈Rn×n表示电动机反电动势常数矩阵,v∈Rn是输入电压的向量,d1、d2、d3均为外在的干扰;
给出期望的参考轨迹qr1(t),并设计包括下式的基于观测器的动态反馈跟踪控制器:
Figure BDA0001796163990000025
Figure BDA0001796163990000026
Figure BDA0001796163990000027
Figure BDA0001796163990000028
Figure BDA0001796163990000029
Figure BDA00017961639900000210
Figure BDA00017961639900000211
对所有有界的初始条件,都确定一个选择增益:k01、k02、k03、k04、k05,使得闭环***的所有状态和信号都是有界的,跟踪误差是一致且最终有界的。
进一步地,所述M(q1)及
Figure BDA0001796163990000031
满足如下性质与假设:
性质1:矩阵M(q1)是对称正定的;
性质2:矩阵
Figure BDA0001796163990000032
是斜对称的;式(1)至式(3)是由未知时变不确定性引起的扰动;
假设1:参数矩阵M(q1)、
Figure BDA0001796163990000033
和G(q1)在式(1)设定为完全未知的;
假设2:式(2)中的参数矩阵表示为:
J(·)=J0+ΔJ(t);
B(·)=B0+ΔB(t);
K(·)=K0+ΔK(t);
KT(·)=KT0+ΔKT(t);
式中,(·)表示一个已知的标称矩阵,Δ(·)表示一个小扰动;此外,还有一个常数0≤εK<1,使得
Figure BDA0001796163990000034
假设3:式(3)中的参数矩阵表示为L(·)=L0+ΔL(t);R和KB假定是完全未知且有界。此外,还有一个常数0≤εL<1,使得||ΔLL-1||2≤εL
进一步地,所述给出期望的参考轨迹qr1(t),目的是找到一个在有不确定性、未建模扰动和外部扰动的情况下柔性关节机器人***的自适应神经网络控制器,引入状态变量:
x1=q1
Figure BDA0001796163990000035
x3=q2
Figure BDA0001796163990000036
x5=I;
将式(1)、式(2)以及式(3)改写成:
Figure BDA0001796163990000037
Figure BDA0001796163990000038
Figure BDA0001796163990000039
Figure BDA0001796163990000041
Figure BDA0001796163990000042
转速表的速度测量容易受噪声影响,不可直接用于控制器的设计,因此使用一个2n维状态估计器来重构缺失的速度信号x2、x4;采用反推设计方法得到控制器并分析其稳定性,从反推技术的观点看,状态变量x2,x3,x4,x5分别视为(4)、(5)、(6)和(7)中的虚拟控制输入。
进一步地,所述基于观测器的动态反馈跟踪控制器的设计具体为:
定义位置跟踪误差
Figure BDA0001796163990000043
和滤波环节的跟踪误差
Figure BDA0001796163990000044
其中
Figure BDA0001796163990000045
表示位置跟踪误差对时间的导数,对于增益k1>0,误差动力学方程
Figure BDA0001796163990000046
改写成:
Figure BDA0001796163990000047
Figure BDA0001796163990000048
其中:
Figure BDA0001796163990000049
Figure BDA00017961639900000410
定义:
Figure BDA00017961639900000411
式中,
Figure BDA00017961639900000412
则F(xe)表示为:
Figure BDA00017961639900000413
则式(10)中的误差动态方程修改为:
Figure BDA00017961639900000414
式中,Fd(qre)=Fd0(qre)+K0qr1
假设上式的不确定项Fd(qre)完全未知,自适应逼近***
Figure BDA0001796163990000051
近似于不确定项Fd(qre);在此
Figure BDA0001796163990000052
包含可调近似参数;根据线性参数化的神经网络模型,将
Figure BDA0001796163990000053
写成:
Figure BDA0001796163990000054
式中,Yd(qre)表示神经网络模型里近似参数的修正系数;
假设4:存在一个最佳逼近参数
Figure BDA0001796163990000055
Figure BDA0001796163990000056
能够无限近似于Fd(qre),将
Figure BDA0001796163990000057
当做最佳逼近误差;存在一个常数εd>0使得||Δωd(·)||2≤εd;由于速度信号x2未测量,滤波环节的跟踪误差
Figure BDA0001796163990000058
不能直接用于控制设计,取η1表示观察状态,
Figure BDA0001796163990000059
表示观测器输出:
Figure BDA00017961639900000510
Figure BDA00017961639900000511
对于控制器增益k01>0,定义观测误差值
Figure BDA00017961639900000512
及一个辅助功能
Figure BDA00017961639900000513
在此γ1>0,经过推导可得:
Figure BDA00017961639900000514
误差动力学方程表示为:
Figure BDA00017961639900000515
使用反推技术,将x3作为式(6)的虚拟控制输入,并选择所需的轨迹,得到:
Figure BDA00017961639900000516
对于控制增益k2>0,令
Figure BDA00017961639900000517
通过将假设4代入式(14)、式(18),推导得到:
Figure BDA00017961639900000518
Figure BDA0001796163990000061
使用表示逼近误差。电机的位置误差动态方程为:
Figure BDA0001796163990000062
由于速度信号x4不能直接用于控制设计,需要开发一个观测器来重建速度信号,取观测器的状态为η2,观测器的输出为
Figure BDA0001796163990000063
其目的是估计x4的值;
定义观察器误差:
Figure BDA0001796163990000064
因此式(22)变为:
Figure BDA0001796163990000065
Figure BDA0001796163990000066
作为式(23)中的虚拟控件输入,用
Figure BDA0001796163990000067
表示期望的轨迹,定义
Figure BDA0001796163990000068
则有:
Figure BDA0001796163990000069
对于控制增益k3>0,分别将式(24)和
Figure BDA00017961639900000610
代入式(23)得:
Figure BDA00017961639900000611
取观察器的输出为:
Figure BDA00017961639900000612
对于观测器增益,若k02>0则误差动力学方程
Figure BDA00017961639900000613
表示为:
Figure BDA00017961639900000614
Figure BDA00017961639900000615
对于控制增益k4>0,修改误差动态方程(27)为:
Figure BDA00017961639900000616
此外,还计算了误差动力学方程
Figure BDA0001796163990000071
Figure BDA0001796163990000072
对于假设2,将J-1B、J-1K、J-1KT表示为:
Figure BDA0001796163990000073
对于ΔJB,ΔJK,ΔJKT,综合(26)和(28),将误差动力学方程(30)写成:
Figure BDA0001796163990000074
从反推技术,将I(t)作为一个虚拟控制输入,I*(t)表示所需电流;设
Figure BDA0001796163990000075
则:
Figure BDA0001796163990000076
将式(32)代入式(31)得:
Figure BDA0001796163990000077
因此,
Figure BDA0001796163990000078
最后,考虑误差动力学方程
Figure BDA0001796163990000079
Figure BDA00017961639900000710
在此
Figure BDA00017961639900000711
考虑到等式
Figure BDA00017961639900000712
Figure BDA00017961639900000713
将误差动力学方程(34)表示为:
Figure BDA0001796163990000081
假设3中
Figure BDA0001796163990000082
又有
Figure BDA0001796163990000083
所以有:
Figure BDA0001796163990000084
因此,
Figure BDA0001796163990000085
Figure BDA0001796163990000086
自适应逼近***
Figure BDA0001796163990000087
近似于
Figure BDA0001796163990000088
所以将
Figure BDA0001796163990000089
写成:
Figure BDA00017961639900000810
假设5:存在一个最佳逼近参数
Figure BDA00017961639900000811
使
Figure BDA00017961639900000812
使用
Figure BDA00017961639900000813
表示最佳逼近误差加上时变不确定项
Figure BDA00017961639900000814
存在常数ε>0使得||ΔωI(·)||2≤εI;选择输入电压:
Figure BDA00017961639900000815
对于控制增益k5>0,则跟踪误差的误差动力学方程
Figure BDA00017961639900000816
表示为:
Figure BDA00017961639900000817
使用
Figure BDA00017961639900000818
表示逼近误差,得到:
Figure BDA00017961639900000819
考虑设备不确定性和外界干扰下柔性关节机器人***,式(1)-式(3)给出期望的参考轨迹,使用qr1(t)设计包括式(38)-式(44)一种基于观测器的动态反馈跟踪控制器,其中,
Figure BDA0001796163990000091
进一步地,在上述过程中,有不确定项ΔF5(·)和
Figure BDA0001796163990000092
假定其是有界的,然后通过适当设计控制增益来补偿这些不确定项的影响;作为一种替代方法,设计滑模控制算法来消除这些不确定项的影响,具体为:
定义在式(38)中的输入电压v修改为:
Figure BDA0001796163990000093
这里
Figure BDA0001796163990000094
由于δ>0,σ>0,且:
Figure BDA0001796163990000095
误差动态方程(37)计算为:
Figure BDA0001796163990000096
进一步地,所述控制***包括超声波电机驱动控制电路,所述超声波电机驱动控制电路包括控制芯片电路和驱动芯片电路,所述光电编码器的信号输出端与所述控制芯片电路的相应输入端相连接,所述控制芯片电路的输出端与所述驱动芯片电路的相应输入端相连接,以驱动所述驱动芯片电路,所述驱动芯片电路的驱动频率调节信号输出端和驱动半桥电路调节信号输出端分别与所述超声波电机的相应输入端相连接。
进一步地,所述联轴器为弹性联轴器。
进一步地,所述超声波电机、光电编码器、力矩传感器分别经超声波电机固定支架、光电编码器固定支架、力矩传感器固定支架固定于一基座上。
考虑不确定性和外界干扰下柔性关节机器人***式(1)-式(3),给出期望qr1(t)的参考轨迹。并设计了一种基于观测器的动态反馈跟踪控制器(39)-(46)。然后,对于所有有界的初始条件,都存在一个选择增益:k01,k02,k03,k04,k05,这样闭环***的所有状态和信号都是有界的,跟踪误差是一致且最终有界的。
与现有技术相比,本发明有以下有益效果:本发明采用自适应神经网络***来逼近不确定的机电动力学行为。构造了一个降阶观测器来估计速度信号。反馈信号只需要测量链路位置和电枢电流。同时本发明利用反推技术,建立了一种基于观测器的鲁棒无速度反馈动态跟踪控制器,使闭环***的所有状态和信号都有界,并可使轨迹跟踪误差尽可能小。
附图说明
图1为本发明实施例方法配套的装置示意图。
图2为本发明实施例的控制***电路示意图。
图中,1为光电编码器,2为光电编码器固定支架,3为超声波电机输出轴,4为超声波电机,5为超声波电机固定支架,6为超声波电机输出轴,7为飞轮惯性负载,8为飞轮惯性负载输出轴,9为弹性联轴器,10为力矩传感器,11为力矩传感器固定支架,12为基座,13为控制芯片电路,14为驱动芯片电路,15、16、17分别为光电编码器输出的A、B、Z相信号,18、19、20、21分别为驱动芯片电路产生的驱动频率调节信号,22为驱动芯片电路产生的驱动半桥电路调节信号,23、24、25、26、27、28分别为控制芯片电路产生的驱动芯片电路的信号,29为超声波电机驱动控制电路。
具体实施方式
下面结合附图及实施例对本发明做进一步说明。
应该指出,以下详细说明都是示例性的,旨在对本申请提供进一步的说明。除非另有指明,本文使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本申请的示例性实施方式。如在这里所使用的,除非上下文另外明确指出,否则单数形式也意图包括复数形式,此外,还应当理解的是,当在本说明书中使用术语“包含”和/或“包括”时,其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
如图1所示,本实施例基于的超声波电机控制***如图1所示,包括基座12和设于基座12上的超声波电机4,所述超声波电机4一侧输出轴3与光电编码器1相连接,另一侧输出轴6与飞轮惯性负载7相连接,所述飞轮惯性负载7的输出轴8经弹性联轴器9与力矩传感器10相连接,所述光电编码器1的信号输出端、所述力矩传感器10的信号输出端分别接至控制***。
上述超声波电机4、光电编码器1、力矩传感器10分别经超声波电机固定支架5、光电编码器固定支架2、力矩传感器固定支架11固定于所述基座12上。
如图2所示,上述控制***包括超声波电机驱动控制电路29,所述超声波电机驱动控制电路29包括控制芯片电路13和驱动芯片电路14,所述光电编码器1的信号输出端与所述控制芯片电路13的相应输入端相连接,所述控制芯片电路13的输出端与所述驱动芯片电路14的相应输入端相连接,以驱动所述驱动芯片电路14,所述驱动芯片电路14的驱动频率调节信号输出端和驱动半桥电路调节信号输出端分别与所述超声波电机4的相应输入端相连接。所述驱动芯片电路14产生驱动频率调节信号和驱动半桥电路调节信号,对超声波电机输出A、B两相PWM的频率、相位及通断进行控制。通过开通及关断PWM波的输出来控制超声波电机的启动和停止运行;通过调节输出的PWM波的频率及两相的相位差来调节电机的最佳运行状态。
本实施例的一种柔性关节机器人的自适应神经网络跟踪控制,整个控制器的***建立在自适应神经网络的基础上,引入滑模控制器抑制消除动态误差,鲁棒自适应控制器确保所有闭环信号都是有界的,最后状态向量的误差和期望的轨迹收敛到给定曲线,从而能获得更好的控制效能。如上所述,在本实施例中,所述控制***的硬件电路包括超声波电机驱动控制电路,所述超声波电机驱动控制电路包括控制芯片电路和驱动芯片电路,所述基于神经网络的超声波电机伺服控制***设于所述控制芯片电路中。
本实施例提供了一种柔性关节机器人的自适应神经网络跟踪控制方法,包括超声波电机,所述超声波电机一侧输出轴与光电编码器相连接,另一侧输出轴与飞轮惯性负载相连接,所述飞轮惯性负载的输出轴经联轴器与力矩传感器相连接,所述光电编码器的信号输出端、所述力矩传感器的信号输出端分别接至控制***;所述控制***建立在自适应神经网络基础上,使用补偿器使得跟踪误差最小,从而能获得更好的控制效能,具体设计为:
所述超声波电机采用由直流伺服电机驱动的柔性关节机械臂运动方程:
Figure BDA0001796163990000121
Figure BDA0001796163990000122
Figure BDA0001796163990000123
式中,其中q1∈Rn、q2∈Rn分别表示连杆角度和电机角度,M(q1)∈Rn×n是对称正定惯性矩阵,
Figure BDA0001796163990000124
表示向心力和科里奥利力的矢量,G(q1)∈Rn表示引力矢量,K=diag{Ki}∈Rn×n是关节刚度系数矩阵,J∈Rn×n是电机惯量矩阵,B∈Rn×n是执行器阻尼矩阵,KT=diag{Kti}∈Rn×n是电机的转矩常数矩阵,I∈Rn是电枢电流矢量,L=diag{Li}∈Rn×n是电感矩阵,R=diag{Ri}∈Rn×n是电阻矩阵,KB=diag{Kbi}∈Rn×n表示电动机反电动势常数矩阵,v∈Rn是输入电压的向量,d1、d2、d3均为外在的干扰;
给出期望的参考轨迹qr1(t),并设计包括下式的基于观测器的动态反馈跟踪控制器:
Figure BDA0001796163990000131
Figure BDA0001796163990000132
Figure BDA0001796163990000133
Figure BDA0001796163990000134
Figure BDA0001796163990000135
Figure BDA0001796163990000136
Figure BDA0001796163990000137
对所有有界的初始条件,都确定一个选择增益:k01、k02、k03、k04、k05,使得闭环***的所有状态和信号都是有界的,跟踪误差是一致且最终有界的。
在本实施例中,所述M(q1)及
Figure BDA0001796163990000138
满足如下性质与假设:
性质1:矩阵M(q1)是对称正定的;
性质2:矩阵
Figure BDA0001796163990000139
是斜对称的;式(1)至式(3)是由未知时变不确定性引起的扰动;
假设1:参数矩阵M(q1)、
Figure BDA00017961639900001310
和G(q1)在式(1)设定为完全未知的;
假设2:式(2)中的参数矩阵表示为:
J(·)=J0+ΔJ(t);
B(·)=B0+ΔB(t);
K(·)=K0+ΔK(t);
KT(·)=KT0+ΔKT(t);
式中,(·)表示一个已知的标称矩阵,Δ(·)表示一个小扰动;此外,还有一个常数0≤εK<1,使得
Figure BDA0001796163990000141
假设3:式(3)中的参数矩阵表示为L(·)=L0+ΔL(t);R和KB假定是完全未知且有界。此外,还有一个常数0≤εL<1,使得||ΔLL-1||2≤εL
在本实施例中,所述给出期望的参考轨迹qr1(t),目的是找到一个在有不确定性、未建模扰动和外部扰动的情况下柔性关节机器人***的自适应神经网络控制器,引入状态变量:
x1=q1
Figure BDA0001796163990000142
x3=q2
Figure BDA0001796163990000143
x5=I;
将式(1)、式(2)以及式(3)改写成:
Figure BDA0001796163990000144
Figure BDA0001796163990000145
Figure BDA0001796163990000146
Figure BDA0001796163990000147
Figure BDA0001796163990000148
转速表的速度测量容易受噪声影响,不可直接用于控制器的设计,因此使用一个2n维状态估计器来重构缺失的速度信号x2、x4;采用反推设计方法得到控制器并分析其稳定性,从反推技术的观点看,状态变量x2,x3,x4,x5分别视为(4)、(5)、(6)和(7)中的虚拟控制输入。
在本实施例中,所述基于观测器的动态反馈跟踪控制器的设计具体为:
定义位置跟踪误差
Figure BDA0001796163990000149
和滤波环节的跟踪误差
Figure BDA0001796163990000151
其中
Figure BDA0001796163990000152
表示位置跟踪误差对时间的导数,对于增益k1>0,误差动力学方程
Figure BDA0001796163990000153
改写成:
Figure BDA0001796163990000154
Figure BDA0001796163990000155
其中:
Figure BDA0001796163990000156
Figure BDA0001796163990000157
定义:
Figure BDA0001796163990000158
式中,
Figure BDA0001796163990000159
则F(xe)表示为:
Figure BDA00017961639900001510
则式(10)中的误差动态方程修改为:
Figure BDA00017961639900001511
式中,Fd(qre)=Fd0(qre)+K0qr1
假设上式的不确定项Fd(qre)完全未知,自适应逼近***
Figure BDA00017961639900001512
近似于不确定项Fd(qre);在此
Figure BDA00017961639900001513
包含可调近似参数;根据线性参数化的神经网络模型,将
Figure BDA00017961639900001514
写成:
Figure BDA00017961639900001515
式中,Yd(qre)表示神经网络模型里近似参数的修正系数;
假设4:存在一个最佳逼近参数
Figure BDA00017961639900001516
Figure BDA00017961639900001517
能够无限近似于Fd(qre),将
Figure BDA00017961639900001518
当做最佳逼近误差;存在一个常数εd>0使得||Δωd(·)||2≤εd;由于速度信号x2未测量,滤波环节的跟踪误差
Figure BDA00017961639900001519
不能直接用于控制设计,取η1表示观察状态,
Figure BDA00017961639900001520
表示观测器输出:
Figure BDA0001796163990000161
Figure BDA0001796163990000162
对于控制器增益k01>0,定义观测误差值
Figure BDA0001796163990000163
及一个辅助功能
Figure BDA0001796163990000164
在此γ1>0,经过推导可得:
Figure BDA0001796163990000165
误差动力学方程表示为:
Figure BDA0001796163990000166
使用反推技术,将x3作为式(6)的虚拟控制输入,并选择所需的轨迹,得到:
Figure BDA0001796163990000167
对于控制增益k2>0,令
Figure BDA0001796163990000168
通过将假设4代入式(14)、式(18),推导得到:
Figure BDA0001796163990000169
Figure BDA00017961639900001610
使用表示逼近误差。电机的位置误差动态方程为:
Figure BDA00017961639900001611
由于速度信号x4不能直接用于控制设计,需要开发一个观测器来重建速度信号,取观测器的状态为η2,观测器的输出为
Figure BDA00017961639900001612
其目的是估计x4的值;
定义观察器误差:
Figure BDA0001796163990000171
因此式(22)变为:
Figure BDA0001796163990000172
Figure BDA0001796163990000173
作为式(23)中的虚拟控件输入,用
Figure BDA0001796163990000174
表示期望的轨迹,定义
Figure BDA0001796163990000175
则有:
Figure BDA0001796163990000176
对于控制增益k3>0,分别将式(24)和
Figure BDA0001796163990000177
代入式(23)得:
Figure BDA0001796163990000178
取观察器的输出为:
Figure BDA0001796163990000179
对于观测器增益,若k02>0则误差动力学方程
Figure BDA00017961639900001710
表示为:
Figure BDA00017961639900001711
Figure BDA00017961639900001712
对于控制增益k4>0,修改误差动态方程(27)为:
Figure BDA00017961639900001713
此外,还计算了误差动力学方程
Figure BDA00017961639900001714
Figure BDA00017961639900001715
对于假设2,将J-1B、J-1K、J-1KT表示为:
Figure BDA00017961639900001716
对于ΔJB,ΔJK,ΔJKT,综合(26)和(28),将误差动力学方程(30)写成:
Figure BDA00017961639900001717
Figure BDA0001796163990000181
从反推技术,将I(t)作为一个虚拟控制输入,I*(t)表示所需电流;设
Figure BDA0001796163990000182
则:
Figure BDA0001796163990000183
将式(32)代入式(31)得:
Figure BDA0001796163990000184
因此,
Figure BDA0001796163990000185
最后,考虑误差动力学方程
Figure BDA00017961639900001819
Figure BDA0001796163990000186
在此
Figure BDA0001796163990000187
考虑到等式
Figure BDA0001796163990000188
Figure BDA0001796163990000189
将误差动力学方程(34)表示为:
Figure BDA00017961639900001810
假设3中
Figure BDA00017961639900001811
又有
Figure BDA00017961639900001812
所以有:
Figure BDA00017961639900001813
因此,
Figure BDA00017961639900001814
Figure BDA00017961639900001815
自适应逼近***
Figure BDA00017961639900001816
近似于
Figure BDA00017961639900001817
所以将
Figure BDA00017961639900001818
写成:
Figure BDA0001796163990000191
假设5:存在一个最佳逼近参数
Figure BDA0001796163990000192
使
Figure BDA0001796163990000193
使用
Figure BDA0001796163990000194
表示最佳逼近误差加上时变不确定项
Figure BDA0001796163990000195
存在常数ε>0使得||ΔωI(·)||2≤εI;选择输入电压:
Figure BDA0001796163990000196
对于控制增益k5>0,则跟踪误差的误差动力学方程
Figure BDA0001796163990000197
表示为:
Figure BDA0001796163990000198
使用
Figure BDA0001796163990000199
表示逼近误差,得到:
Figure BDA00017961639900001910
考虑设备不确定性和外界干扰下柔性关节机器人***,式(1)-式(3)给出期望的参考轨迹,使用qr1(t)设计包括式(38)-式(44)一种基于观测器的动态反馈跟踪控制器,其中,
Figure BDA00017961639900001911
在本实施例中,在上述过程中,有不确定项ΔF5(·)和
Figure BDA00017961639900001912
假定其是有界的,然后通过适当设计控制增益来补偿这些不确定项的影响;作为一种替代方法,设计滑模控制算法来消除这些不确定项的影响,具体为:
定义在式(38)中的输入电压v修改为:
Figure BDA00017961639900001913
这里
Figure BDA0001796163990000201
由于δ>0,σ>0,且:
Figure BDA0001796163990000202
误差动态方程(37)计算为:
Figure BDA0001796163990000203
本实施例考虑不确定性和外界干扰下柔性关节机器人***式(1)-式(3),给出期望qr1(t)的参考轨迹。并设计了一种基于观测器的动态反馈跟踪控制器(39)-(46)。然后,对于所有有界的初始条件,都存在一个选择增益:k01,k02,k03,k04,k05,这样闭环***的所有状态和信号都是有界的,跟踪误差是一致且最终有界的。
以上所述仅为本发明的较佳实施例,凡依本发明申请专利范围所做的均等变化与修饰,皆应属本发明的涵盖范围。

Claims (4)

1.一种柔性关节机器人的自适应神经网络跟踪控制方法,其特征在于:包括超声波电机,所述超声波电机一侧输出轴与光电编码器相连接,另一侧输出轴与飞轮惯性负载相连接,所述飞轮惯性负载的输出轴经联轴器与力矩传感器相连接,所述光电编码器的信号输出端、所述力矩传感器的信号输出端分别接至控制***;所述控制***建立在自适应神经网络基础上,使用补偿器使得跟踪误差最小,从而能获得更好的控制效能,具体设计为:
所述超声波电机采用由直流伺服电机驱动的柔性关节机械臂运动方程:
Figure FDA0002430281540000011
Figure FDA0002430281540000012
Figure FDA0002430281540000013
式中,其中q1∈Rn、q2∈Rn分别表示连杆角度和电机角度,M(q1)∈Rn×n是对称正定惯性矩阵,
Figure FDA0002430281540000014
表示向心力和科里奥利力的矢量,G(q1)∈Rn表示引力矢量,K=diag{Ki}∈Rn×n是关节刚度系数矩阵,J∈Rn×n是电机惯量矩阵,B∈Rn×n是执行器阻尼矩阵,KT=diag{Kti}∈Rn×n是电机的转矩常数矩阵,I∈Rn是电枢电流矢量,L=diag{Li}∈Rn×n是电感矩阵,R=diag{Ri}∈Rn×n是电阻矩阵,KB=diag{Kbi}∈Rn×n表示电动机反电动势常数矩阵,v∈Rn是输入电压的向量,d1、d2、d3均为外在的干扰;
给出期望的参考轨迹qr1(t),并设计包括下式的基于观测器的动态反馈跟踪控制器:
Figure FDA0002430281540000015
Figure FDA0002430281540000016
Figure FDA0002430281540000017
Figure FDA0002430281540000018
Figure FDA0002430281540000019
Figure FDA00024302815400000110
Figure FDA0002430281540000021
对所有有界的初始条件,都确定一个选择增益:k01、k02、k03、k04、k05,使得闭环***的所有状态和信号都是有界的,跟踪误差是一致且最终有界的;
其中,所述M(q1)及
Figure FDA0002430281540000022
满足如下性质与假设:
性质1:矩阵M(q1)是对称正定的;
性质2:矩阵
Figure FDA0002430281540000023
是斜对称的;式(1)至式(3)是由未知时变不确定性引起的扰动;
假设1:参数矩阵M(q1)、
Figure FDA0002430281540000024
和G(q1)在式(1)设定为完全未知的;
假设2:式(2)中的参数矩阵表示为:
J(·)=J0+ΔJ(t);
B(·)=B0+ΔB(t);
K(·)=K0+ΔK(t);
KT(·)=KT0+ΔKT(t);
式中,(·)表示一个已知的标称矩阵,Δ(·)表示一个小扰动;此外,还有一个常数0≤εK<1,使得
Figure FDA0002430281540000025
假设3:式(3)中的参数矩阵表示为L(·)=L0+ΔL(t);R和KB假定是完全未知且有界;此外,还有一个常数0≤εL<1,使得||ΔLL-1||2≤εL
其中,所述给出期望的参考轨迹qr1(t),目的是找到一个在有不确定性、未建模扰动和外部扰动的情况下柔性关节机器人***的自适应神经网络控制器,引入状态变量:
Figure FDA0002430281540000026
将式(1)、式(2)以及式(3)改写成:
Figure FDA0002430281540000027
Figure FDA0002430281540000028
Figure FDA0002430281540000031
Figure FDA0002430281540000032
Figure FDA0002430281540000033
转速表的速度测量容易受噪声影响,不可直接用于控制器的设计,因此使用一个2n维状态估计器来重构缺失的速度信号x2、x4;采用反推设计方法得到控制器并分析其稳定性,从反推技术的观点看,状态变量x2,x3,x4,x5分别视为(4)、(5)、(6)和(7)中的虚拟控制输入;
其中,所述基于观测器的动态反馈跟踪控制器的设计具体为:
定义位置跟踪误差
Figure FDA0002430281540000034
和滤波环节的跟踪误差
Figure FDA0002430281540000035
其中
Figure FDA0002430281540000036
表示位置跟踪误差对时间的导数,对于增益k1>0,误差动力学方程
Figure FDA0002430281540000037
改写成:
Figure FDA0002430281540000038
Figure FDA0002430281540000039
其中:
Figure FDA00024302815400000310
Figure FDA00024302815400000311
定义:
Figure FDA00024302815400000312
式中,
Figure FDA00024302815400000313
则F(xe)表示为:
Figure FDA00024302815400000314
则式(10)中的误差动态方程修改为:
Figure FDA00024302815400000315
式中,Fd(qre)=Fd0(qre)+K0qr1
假设上式的不确定项Fd(qre)完全未知,自适应逼近***
Figure FDA00024302815400000316
近似于不确定项Fd(qre);在此
Figure FDA0002430281540000041
包含可调近似参数;根据线性参数化的神经网络模型,将
Figure FDA0002430281540000042
写成:
Figure FDA0002430281540000043
式中,Yd(qre)表示神经网络模型里近似参数的修正系数;
假设4:存在一个最佳逼近参数
Figure FDA0002430281540000044
能够无限近似于Fd(qre),将
Figure FDA0002430281540000045
当做最佳逼近误差;存在一个常数εd>0使得||Δωd(·)||2≤εd;由于速度信号x2未测量,滤波环节的跟踪误差
Figure FDA0002430281540000046
不能直接用于控制设计,取η1表示观察状态,
Figure FDA0002430281540000047
表示观测器输出:
Figure FDA0002430281540000048
Figure FDA0002430281540000049
对于控制器增益k01>0,定义观测误差值
Figure FDA00024302815400000410
及一个辅助功能
Figure FDA00024302815400000411
在此γ1>0,经过推导可得:
Figure FDA00024302815400000412
误差动力学方程表示为:
Figure FDA00024302815400000413
使用反推技术,将x3作为式(6)的虚拟控制输入,并选择所需的轨迹,得到:
Figure FDA00024302815400000414
对于控制增益k2>0,令
Figure FDA00024302815400000415
通过将假设4代入式(14)、式(18),推导得到:
Figure FDA00024302815400000416
Figure FDA00024302815400000417
电机的位置误差动态方程为:
Figure 1
由于速度信号x4不能直接用于控制设计,需要开发一个观测器来重建速度信号,取观测器的状态为η2,观测器的输出为
Figure FDA0002430281540000051
其目的是估计x4的值;
定义观察器误差:
Figure FDA0002430281540000052
因此式(22)变为:
Figure FDA0002430281540000053
Figure FDA0002430281540000054
作为式(23)中的虚拟控件输入,用
Figure FDA0002430281540000055
表示期望的轨迹,定义
Figure FDA0002430281540000056
则有:
Figure FDA0002430281540000057
对于控制增益k3>0,分别将式(24)和
Figure FDA0002430281540000058
代入式(23)得:
Figure FDA0002430281540000059
取观察器的输出为:
Figure FDA00024302815400000510
对于观测器增益,若k02>0则误差动力学方程
Figure FDA00024302815400000511
表示为:
Figure FDA00024302815400000512
Figure FDA00024302815400000513
对于控制增益k4>0,修改误差动态方程(27)为:
Figure FDA00024302815400000514
此外,还计算了误差动力学方程
Figure FDA00024302815400000515
Figure FDA00024302815400000516
对于假设2,将J-1B、J-1K、J-1KT表示为:
Figure FDA00024302815400000517
对于ΔJB,ΔJK,ΔJKT,综合(26)和(28),将误差动力学方程(30)写成:
Figure FDA00024302815400000518
从反推技术,将I(t)作为一个虚拟控制输入,I*(t)表示所需电流;设
Figure FDA0002430281540000061
则:
Figure FDA0002430281540000062
将式(32)代入式(31)得:
Figure FDA0002430281540000063
因此,
Figure FDA0002430281540000064
最后,考虑误差动力学方程
Figure FDA0002430281540000065
Figure FDA0002430281540000066
在此
Figure FDA0002430281540000067
考虑到等式
Figure FDA0002430281540000068
Figure FDA0002430281540000069
将误差动力学方程(34)表示为:
Figure FDA00024302815400000610
假设3中
Figure FDA00024302815400000611
又有
Figure FDA00024302815400000612
所以有:
Figure FDA00024302815400000613
因此,
Figure FDA00024302815400000614
Figure FDA00024302815400000615
自适应逼近***
Figure FDA00024302815400000616
近似于
Figure FDA00024302815400000617
所以将
Figure FDA00024302815400000618
写成:
Figure FDA00024302815400000619
假设5:存在一个最佳逼近参数
Figure FDA00024302815400000620
使
Figure FDA00024302815400000621
使用
Figure FDA00024302815400000622
表示最佳逼近误差加上时变不确定项
Figure FDA00024302815400000623
存在常数ε>0使得||ΔωI(·)||2≤εI;选择输入电压:
Figure FDA0002430281540000071
对于控制增益k5>0,则跟踪误差的误差动力学方程
Figure FDA0002430281540000072
表示为:
Figure FDA0002430281540000073
使用
Figure FDA0002430281540000074
表示逼近误差,得到:
Figure FDA0002430281540000075
考虑设备不确定性和外界干扰下柔性关节机器人***,式(1)-式(3)给出期望的参考轨迹,使用qr1(t)设计包括式(38)-式(44)一种基于观测器的动态反馈跟踪控制器,其中,
Figure FDA0002430281540000076
γ1>0,γ2>0;
其中,在上述过程中,有不确定项ΔF5(·)和
Figure FDA0002430281540000077
假定其是有界的,然后通过适当设计控制增益来补偿这些不确定项的影响;作为一种替代方法,设计滑模控制算法来消除这些不确定项的影响,具体为:
定义在式(38)中的输入电压v修改为:
Figure FDA0002430281540000078
这里
Figure FDA0002430281540000079
由于δ>0,σ>0,且:
Figure FDA00024302815400000710
误差动态方程(37)计算为:
Figure 3
2.根据权利要求1所述的一种柔性关节机器人的自适应神经网络跟踪控制方法,其特征在于:所述控制***包括超声波电机驱动控制电路,所述超声波电机驱动控制电路包括控制芯片电路和驱动芯片电路,所述光电编码器的信号输出端与所述控制芯片电路的相应输入端相连接,所述控制芯片电路的输出端与所述驱动芯片电路的相应输入端相连接,以驱动所述驱动芯片电路,所述驱动芯片电路的驱动频率调节信号输出端和驱动半桥电路调节信号输出端分别与所述超声波电机的相应输入端相连接。
3.根据权利要求1所述的一种柔性关节机器人的自适应神经网络跟踪控制方法,其特征在于:所述联轴器为弹性联轴器。
4.根据权利要求1所述的一种柔性关节机器人的自适应神经网络跟踪控制方法,其特征在于:所述超声波电机、光电编码器、力矩传感器分别经超声波电机固定支架、光电编码器固定支架、力矩传感器固定支架固定于一基座上。
CN201811057430.1A 2018-09-11 2018-09-11 一种柔性关节机器人的自适应神经网络跟踪控制方法 Active CN108818541B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811057430.1A CN108818541B (zh) 2018-09-11 2018-09-11 一种柔性关节机器人的自适应神经网络跟踪控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811057430.1A CN108818541B (zh) 2018-09-11 2018-09-11 一种柔性关节机器人的自适应神经网络跟踪控制方法

Publications (2)

Publication Number Publication Date
CN108818541A CN108818541A (zh) 2018-11-16
CN108818541B true CN108818541B (zh) 2020-09-04

Family

ID=64149317

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811057430.1A Active CN108818541B (zh) 2018-09-11 2018-09-11 一种柔性关节机器人的自适应神经网络跟踪控制方法

Country Status (1)

Country Link
CN (1) CN108818541B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110107231B (zh) * 2019-06-06 2024-03-29 吉林大学 一种自适应神经网络自动猫道及其控制方法
CN111965976B (zh) * 2020-08-06 2021-04-23 北京科技大学 基于神经网络观测器的机器人关节滑模控制方法及***
CN112180729B (zh) * 2020-10-09 2023-03-03 电子科技大学中山学院 一种机械臂自适应轨迹跟踪控制方法
US11772264B2 (en) * 2020-11-18 2023-10-03 Dibi (Chongqing) Intelligent Technology Research Institute Co., Ltd. Neural network adaptive tracking control method for joint robots
CN113459097B (zh) * 2021-06-28 2022-11-18 北京机械设备研究所 无传感器柔性拉线阻抗控制方法、装置、介质及电子设备
CN114280924B (zh) * 2021-12-21 2023-10-31 南京理工大学 基于模糊观测器的柔性关节机器人***事件触发控制方法
CN114895564B (zh) * 2022-05-20 2023-05-02 南昌大学 一种电驱动柔性关节机械臂自适应神经网络控制器设计方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0949869A (ja) * 1995-08-09 1997-02-18 Fujitsu Ltd 目標追尾制御装置
CN103336523B (zh) * 2013-06-19 2016-03-16 闽江学院 超声波电机伺服控制***的控制策略
CN103513181B (zh) * 2013-07-10 2016-03-02 闽江学院 超声波电机瞬态特性测试装置及其控制***
CN105099263B (zh) * 2015-09-28 2017-03-08 闽江学院 基于广义回归神经网络辨识器的超声波电机位置与速度控制***
CN106078741B (zh) * 2016-06-21 2018-04-13 华南理工大学 基于确定学习理论的性能受限柔性机械臂控制方法

Also Published As

Publication number Publication date
CN108818541A (zh) 2018-11-16

Similar Documents

Publication Publication Date Title
CN108818541B (zh) 一种柔性关节机器人的自适应神经网络跟踪控制方法
Wang et al. Robust sliding mode control for robots driven by compliant actuators
Xie Sliding-mode-observer-based adaptive control for servo actuator with friction
US7211979B2 (en) Torque-position transformer for task control of position controlled robots
JP6097174B2 (ja) ロボット制御装置
US9166512B2 (en) Motor control device provided with feedforward control
CN110467111B (zh) 桥式起重机的控制
Zain et al. Hybrid learning control schemes with input shaping of a flexible manipulator system
JP2005301508A (ja) 制御装置
JP3981773B2 (ja) ロボット制御装置
Hung et al. Adaptive control for nonlinearly parameterized uncertainties in robot manipulators
CN108942928A (zh) 一种基于约束力鲁棒伺服控制的欠驱动柔性机械臂***
Szczepanski et al. PMSM drive with adaptive state feedback speed controller
KR20000048977A (ko) 수치 제어되는 탄성 다중 매스 시스템의 모우먼트를 제어하기 위한 방법 및 제어 장치
Katsura et al. Absolute stabilization of multimass resonant system by phase-lead compensator based on disturbance observer
JP3936242B2 (ja) サーボモータ制御用コントローラにおけるゲイン設定法、コントローラの有効性検証法およびロボット制御法
Izadbakhsh et al. Endpoint perfect tracking control of robots—a robust non inversion-based approach
Aksman et al. Force estimation based compliance control of harmonically driven manipulators
Rsetam et al. Sliding mode control based on high-order extended state observer for flexible joint robot under time-varying disturbance
JP4008207B2 (ja) ロボット制御方法および制御装置
Lukichev et al. Speed control in PMSM drive with non-stiff load and unknown parameters using PI-and fuzzy adaptive PID controllers
Yakub et al. Vibration control of two-mass rotary system using improved NCTF controller for positioning systems
Yakub et al. Practical control for two-mass positioning systems in presence of saturation
CN112147894A (zh) 基于运动学和动力学模型的轮式移动机器人主动控制方法
Abdullah et al. Vision based automatic steering control using a PID controller

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