CN107825424A - 一种高速机械手减少残余振动的非对称s型轨迹规划方法 - Google Patents

一种高速机械手减少残余振动的非对称s型轨迹规划方法 Download PDF

Info

Publication number
CN107825424A
CN107825424A CN201711053996.2A CN201711053996A CN107825424A CN 107825424 A CN107825424 A CN 107825424A CN 201711053996 A CN201711053996 A CN 201711053996A CN 107825424 A CN107825424 A CN 107825424A
Authority
CN
China
Prior art keywords
acceleration
displacement
asymmetric
impact
stage
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
Application number
CN201711053996.2A
Other languages
English (en)
Other versions
CN107825424B (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.)
Hebei University of Technology
Original Assignee
Hebei University of 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 Hebei University of Technology filed Critical Hebei University of Technology
Priority to CN201711053996.2A priority Critical patent/CN107825424B/zh
Publication of CN107825424A publication Critical patent/CN107825424A/zh
Application granted granted Critical
Publication of CN107825424B publication Critical patent/CN107825424B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/1656Programme controls characterised by programming, planning systems for manipulators
    • B25J9/1664Programme controls characterised by programming, planning systems for manipulators characterised by motion, path, trajectory planning

Abstract

本发明公开了一种高速机械手减少残余振动的非对称S型轨迹规划方法,现有的S型轨迹规划方法只考虑了运行轨迹的平滑性,没有考虑到机械手自身特性对操作精度的影响及不同工况对机械手定位精度的要求,导致实际操作精度无法满足要求,不能兼顾不同工况对平稳性和生产效率最大化的问题,并且需要单独增加抑制残余振动的阶段,增加了方法的复杂性。本发明充分考虑到机械手自身特性对运行精度的影响和在不同工况下对定位精度的要求,兼顾运行效率和操作精度,灵活性高。

Description

一种高速机械手减少残余振动的非对称S型轨迹规划方法
技术领域
本发明涉及机械工程控制技术领域,具体而言,尤其涉及一种高速机械手减少残余振动的非对称S型轨迹规划方法。
背景技术
为了提高工业生产效率,大量高速搬运机械手应用于制造业加工。为了降低启停阶段由于速度的突变引起的振动,高速机械手快速启停主要采用S曲线进行速度规划,通过选取不同的冲击值(加加速度),获得不同形状的S型轨迹曲线。如专利“一种高速机构减少残余振动的S型运动曲线规划方法”所述,常规S曲线为了降低机械手的残余振动需要单独增加抑制残余振动的阶段。目前的S型曲线规划方法加速行程和减速行程各参数相同属于对称S曲线,当前的S型曲线侧重于保证运动的平稳性,很少考虑驱动负载、运行速度、工况等因素的影响。由于运行速度的提高会增大机械手的柔性残余振动对机械手的运动精度的影响。对于高速高精度的机械手,在进行速度规划时不仅要减少运行时间提高运行效率,还要保证所需的定位精度要求,通常运行时间越短冲击值越大,残余振动要经过更长的时间才能衰减到指定的位置。
当前的S型轨迹规划方法,即主要考虑了曲线运行平稳性的S型曲线规划方法,未充分考虑不同机械手自身特性对运行精度的影响和在不同工况下对定位精度的要求。导致规划灵活性不够,不能够兼顾运行效率和操作精度。
发明内容
为解决现有技术的不足,本发明提供一种高速机械手减少残余振动的非对称S型轨迹规划方法。现有的S型轨迹规划方法只考虑了运行轨迹的平滑性,没有考虑到机械手自身特性对操作精度的影响及不同工况对机械手定位精度的要求,导致实际操作精度无法满足要求,不能兼顾不同工况对平稳性和生产效率最大化的问题,并且需要单独增加抑制残余振动的阶段,增加了方法的复杂性。
本发明解决所述技术问题的技术方案是:设计一种高速机械手减少残余振动的非对称S型轨迹规划方法,该方法包括以下步骤:
步骤一、首先确定机械手运行的最大速度、最大加速度以及加速形成的冲击值、减速形成的冲击值的限制条件。
步骤二、求出非对称S曲线的各阶段时间节点的速度、加速度值。
步骤三、由各阶段之间的连续性以及边界条件确定出各阶段的运动学参数。
步骤四、通过已知条件确定出短位移、中位移、长位移下的数学模型。
步骤五、通过改变非对称系数R获得各种非对称度系数下的非对称S曲线。
步骤六、通过合理的选取采样点进行三次样条插值拟合出各个关节的输入七阶非对称S曲线分为冲击J1进行的加加速运行阶段(T1);匀加速阶段(T2);以冲击J1进行的减加速运行阶段(T3),匀加速阶段(T4);以冲击J2进行的减加速运行阶段(T5);匀减速阶段(T6);以冲击J2进行的减减速阶段(T7)。
进一步,根据非对称S曲线的运动特点,加速行程和减速行程的冲击值为常数,匀加速阶段和匀减速阶段的冲击值为0,匀速运行阶段的加速度值为0,减减速阶段结束时速度、加速度为零;因此,有以下的等式约束:
T7=RT1
T6=RT2
T5=RT3
J2=R2J1
因此通过调节参数R获得不同非对称度下的非对称S曲线。评价指标和残余振动的振幅Ap、残余振动衰减时间Tp的关系如下式所示:
在目标位移大小发生变化时,匀速阶段和匀加速阶段的时间可能等于零。通过调节位移获得在短位移、中位移下的数学模型,依据工况和机械手的动力学特性对相关参数进行设置,然后通过实验检测装置不同位移下的最优的非对称系数R,最后确定机械手的轨迹规划最优参数。
与现有技术相比,本发明的有益效果在于:现有的S型轨迹规划方法只考虑了运行轨迹的平滑性,没有考虑到机械手自身特性对操作精度的影响及不同工况对机械手定位精度的要求,导致实际操作精度无法满足要求,不能兼顾不同工况对平稳性和生产效率最大化的问题,并且需要单独增加抑制残余振动的阶段,增加了方法的复杂性。本发明充分考虑到机械手自身特性对运行精度的影响和在不同工况下对定位精度的要求,兼顾运行效率和操作精度,灵活性高。
附图说明
图1是本发明一种实施例提供的减少残余振动的非对称S型规划曲线图;
图2是中位移和短位移下非对称S型曲线与对称S型曲线的对比图,其中图2(a)是中位移下非对称S型曲线与对称S型曲线的对比图,图2(b)是短位移下非对称S型曲线与对称S型曲线的对比图;
图3是本发明一种实施例的非对称S曲线轨迹规划方法的流程图;
图4是本发明一种实施例的减少残余振动的非对称S曲线规划方法的装置图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实例,对本发明进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
图2演示了本发明一种实施例提供的减少残余振动的非对称S曲线的曲线图,图3为曲线规划方法的流程图,为了便于说明仅仅展示了与本发明相关的部分。
本发明提供的一种高速机械手减少残余振动的非对称S型轨迹规划方法,该方法包括以下步骤:
步骤一、首先确定最大速度Ax,最大加速度值Vx,加速形成的冲击值(加加速度)为J1,减速形成的冲击值为J2,由关系值得到J1=J2/R2
步骤二、求出各阶段时间节点的加速度值、速度值。
加加速阶段T1
A1=J1t
匀加速阶段T2
A2=Ax
减加速阶段T3
A3=-J1t+Ax
匀速阶段T4
A4=0
V4=Vx
加减速阶段T5
匀减速阶段T6
减减速阶段T7
A7=0
步骤三、由各阶段之间的连续性及边界条件确定出各阶段的参数
T7=RT1
T6=RT2
T5=RT3
J2=R2J1
步骤四、通过已知参数确定出短位移、中位移、长位移下的数学模型
在Vx和Ax限制条件不变的情况下,目标位移P发生变化时T1和T4可能等于零。为T1,T4恰好同时等于零时的目标位移,为恰好T4等于零时的目标位移。
时为短位移,T2=T4=0
时为中位移,
时为长位移,
步骤五、通过改变非对称系数R便获得了不同目标位移下的非对称曲线。
步骤六、通过合理的选取非对称S曲线的采样点进行三次样条插值拟合出机械手各个关节的输入。
长位移的非对称S曲线的行程为:以冲击(加加速度)J1进行的加加速阶段(T1),匀加速阶段(T2),以冲击J1进行的减加速阶段(T3),匀速阶段(T4),以冲击J2进行的减加速阶段(T5),匀减速阶段(T6),以冲击J2进行的减减速阶段(T7)。
中位移的非对称S曲线的行程为:以冲击(加加速度)J1进行的加加速阶段(T1),匀加速阶段(T2),以冲击J1进行的减加速阶段(T3),以冲击J2进行的减加速阶段(T5),匀减速阶段(T6),以冲击J2进行的减减速阶段(T7)。
短位移的非对称S曲线的行程为:以冲击(加加速度)J1进行的加加速阶段(T1),以冲击J1进行的减加速阶段(T3),以冲击J2进行的减加速阶段(T5),以冲击J2进行的减减速阶段(T7)。
本发明的轨迹规划方案,依据非对称S曲线的运行特点得知加速行程和减速行程的冲击值为常数,匀速运行阶段的加速度为零,有如下等式的约束条件:
T7=RT1
T6=RT2
T5=RT3
J2=R2J1
T1=T3,T5=T7
因此减速行程的参数可以由加速行程和非对称参数R确定。
为了详细叙述发明内容,下面结合具体实例来对规划方法进行说明:
本文中涉及到的所有变量的名称含义见表1(各阶段的图像请参照配图1)
表1变量总表
首先根据机械手的机械结构特性和搬运流程确定出相关的运行参数极值如最大加速度、最大加速度、加速行程的冲击等。得到的参数列表如表2所示:
表2参数变量取值表
为不同冲击下的加加速阶段和匀加速阶段的时间
得到
得到
表示匀加速阶段刚好为0时的位移,即ta=tv=0时的位移。
表示匀速阶段刚好等于0时的位移,即时的位移。
得到
得到
为了计算方便,上述结果均做了取整处理。
时为短位移,当时为中位移,当时为长位移。
短位移下ta=tv=0,
中位移下
长位移下
为了计算方便短位移取值为100000,中位移取值为300000,长位移取值为700000。将上面计算得到的各种位移情况下的P、J1、tj、ta、tv分别代入各阶段的表达式,便得到了各种位移情况下的非对称S曲线的数学模型。
非对称S曲线的各阶段的具体表达式:
<1>加加速阶段0~t1(tj=t1-t0):
A1=J1t
<2>匀加速阶段t1~t2(ta=t2-t1):
A2=Ax
<3>减加速阶段t2~t3(tj=t3-t2):
A3=-J1T1+A,
其中T1=t-ta-tj
<4>匀速阶段t3~t4(tv=t4-t3):
A4=0
V4=Vx
<5>加减速阶段t4~t5(Rtj=t5-t4):
其中T2=t-ta-2tj-tv
<6>匀减速阶段t5~t6(Rta=t6-t5):
其中T3=t-ta-(2+R)tj-tv
<7>减减速阶段t6~t7(Rtj=t7-t6):
其中T4=t-(1+R)ta-(2+R)tj-tv
通过上述计算求得短位移、中位移、长位移下加速度、速度、位移关于非对称系数R的数学模型。通过调节R便获得了不同非对称度下的轨迹规划模型。实施二:
由于不同机械手的固有运动学参数不同,因此对应的最优非对称系数R的取值也有所不同,所以有必要通过一套装置来检测系统最优的非对称系数。
系统的评价指标和残余振动的振幅Ap、残余振动衰减时间Tp的关系如下式所示:
由于不同的系统对精度的要求不同,所以通过调整权重系数α不同的评价指标,使该方法有更广泛的应用。
在目标位移发生变化时,匀速阶段T4,匀加速阶段T2可能等于零,当目标位移为长位移时T2≠0,T4≠0,当目标位移为中位移时T2≠0,T4=0,当目标位移为短位移时T2=0,T4=0。通过轨迹规划方法获得不同目标位移下的非对称S曲线的数学模型。通过调整参数获得不同的R对应的关节输入。由于检测过程中会存在干扰因素,所以实验装置需要对输入信号中的数据进行筛选,剔除掉不合理的数据。
由非对称S曲线的对称性可知T1=T3,T5=T7,高速机械手减少残余振动的非对称S曲线的规划方法的模型为:
Find(T1,T2,T4,R)
Objective:
Subject to:
下面结合附图及具体应用对本发明的具体原理进一步进行描述。
本发明的高速机械手的非对称S型速度曲线的模型如图1所示,在不同的位移下非对称S曲线的速度曲线的如图2所示。
长位移的非对称S曲线的行程为:以冲击(加加速度)J1进行的加加速阶段(T1),匀加速阶段(T2),以冲击J1进行的减加速阶段(T3),匀速阶段(T4),以冲击J2进行的减加速阶段(T5),匀减速阶段(T6),以冲击J2进行的减减速阶段(T7)。为了降低残余振动的影响,本发明考虑了不同位移下非对称系数R对残余振动的影响,并提出了定位精度的评价指标,通过实验测定出机械手不同工况下的最优指标。
依据非对称S曲线的运行特点得知加速行程和减速行程的冲击值为常数,匀速运行阶段的加速度为零。
因此减速行程的参数可以由加速行程和非对称参数R确定。评价指标和残余振动的振幅Ap、残余振动衰减时间Tp的关系如下式所示:
在目标位移发生变化时,匀速阶段T4,匀加速阶段T2可能等于零,当目标位移为长位移时T2≠0,T4≠0,当目标位移为中位移时T2≠0,T4=0,当目标位移为短位移时T2=0,T4=0。
由非对称S曲线的对称性可知T1=T3,T5=T7,高速机械手减少残余振动的非对称S曲线的规划方法的模型为:
Find(T1,T2,T4,R)
Objective:
Subject to:
本发明未述及之处适用于现有技术。

Claims (2)

1.一种高速机械手减少残余振动的非对称S型轨迹规划方法,其特征在于,该方法包括以下步骤:
步骤一、首先确定机械手运行的最大速度、最大加速度以及加速形成的冲击值、减速形成的冲击值的限制条件;
步骤二、求出非对称S曲线的各阶段时间节点的速度、加速度值;
步骤三、由各阶段之间的连续性以及边界条件确定出各阶段的运动学参数;
步骤四、通过已知条件确定出短位移、中位移、长位移下的数学模型;
步骤五、通过改变非对称系数R获得各种非对称度系数下的非对称S曲线;
步骤六、通过合理的选取采样点进行三次样条插值拟合出各个关节的输入七阶非对称S曲线分为冲击J1进行的加加速运行阶段(T1);匀加速阶段(T2);以冲击J1进行的减加速运行阶段(T3),匀加速阶段(T4);以冲击J2进行的减加速运行阶段(T5);匀减速阶段(T6);以冲击J2进行的减减速阶段(T7);
根据非对称S曲线的运动特点,加速行程和减速行程的冲击值为常数,匀加速阶段和匀减速阶段的冲击值为0,匀速运行阶段的加速度值为0,减减速阶段结束时速度、加速度为零;因此,有以下的等式约束:
T7=RT1
T6=RT2
T5=RT3
J2=R2J1
因此通过调节参数R获得不同非对称度下的非对称S曲线;评价指标和残余振动的振幅Ap、残余振动衰减时间Tp的关系如下式所示:
在目标位移大小发生变化时,匀速阶段和匀加速阶段的时间可能等于零;通过调节位移获得在短位移、中位移下的数学模型,依据工况和机械手的动力学特性对相关参数进行设置,然后通过实验检测装置不同位移下的最优的非对称系数R,最后确定机械手的轨迹规划最优参数。
2.根据权利要求1所述的一种高速机械手减少残余振动的非对称S型轨迹规划方法,其特征在于,该方法包括以下步骤:
步骤一、首先确定最大速度Ax,最大加速度值Vx,加速形成的冲击值为J1,减速形成的冲击值为J2,由关系值得到J1=J2/R2
步骤二、求出各阶段时间节点的加速度值、速度值
加加速阶段T1
A1=J1t
<mrow> <msub> <mi>V</mi> <mn>1</mn> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msub> <mi>J</mi> <mn>1</mn> </msub> <msup> <mi>t</mi> <mn>2</mn> </msup> </mrow>
匀加速阶段T2
A2=Ax
<mrow> <msub> <mi>V</mi> <mn>2</mn> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msub> <mi>J</mi> <mn>1</mn> </msub> <msubsup> <mi>t</mi> <mi>j</mi> <mn>2</mn> </msubsup> <mo>+</mo> <msub> <mi>A</mi> <mi>x</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <msub> <mi>t</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> </mrow>
减加速阶段T3
A3=-J1t+Ax
<mrow> <msub> <mi>V</mi> <mn>3</mn> </msub> <mo>=</mo> <mo>-</mo> <msub> <mi>J</mi> <mn>1</mn> </msub> <mfrac> <msup> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <msub> <mi>T</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>T</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mn>2</mn> </mfrac> <mo>+</mo> <msub> <mi>A</mi> <mi>x</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <msub> <mi>T</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>T</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>A</mi> <mi>x</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>t</mi> <mi>a</mi> </msub> <mo>+</mo> <mfrac> <msub> <mi>t</mi> <mi>j</mi> </msub> <mn>2</mn> </mfrac> <mo>)</mo> </mrow> </mrow>
匀速阶段T4
A4=0
V4=Vx
加减速阶段T5
<mrow> <msub> <mi>A</mi> <mn>5</mn> </msub> <mo>=</mo> <mo>-</mo> <mfrac> <msub> <mi>J</mi> <mn>1</mn> </msub> <msup> <mi>R</mi> <mn>2</mn> </msup> </mfrac> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mn>2</mn> <msub> <mi>T</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>T</mi> <mn>3</mn> </msub> <mo>-</mo> <msub> <mi>T</mi> <mn>4</mn> </msub> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mi>V</mi> <mn>5</mn> </msub> <mo>=</mo> <msub> <mi>V</mi> <mi>x</mi> </msub> <mo>-</mo> <mfrac> <mi>J</mi> <mrow> <mn>2</mn> <msup> <mi>R</mi> <mn>2</mn> </msup> </mrow> </mfrac> <msup> <mrow> <mo>(</mo> <mi>T</mi> <mo>-</mo> <mn>2</mn> <msub> <mi>T</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>T</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>T</mi> <mn>4</mn> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow>
匀减速阶段T6
<mrow> <msub> <mi>A</mi> <mn>6</mn> </msub> <mo>=</mo> <mo>-</mo> <mfrac> <msub> <mi>A</mi> <mi>x</mi> </msub> <mi>R</mi> </mfrac> </mrow>
<mrow> <msub> <mi>V</mi> <mn>6</mn> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msub> <mi>A</mi> <mi>x</mi> </msub> <msub> <mi>T</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>A</mi> <mi>x</mi> </msub> <msub> <mi>T</mi> <mn>2</mn> </msub> <mo>-</mo> <mfrac> <msub> <mi>A</mi> <mi>x</mi> </msub> <mi>R</mi> </mfrac> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <msub> <mi>T</mi> <mn>2</mn> </msub> <mo>-</mo> <mo>(</mo> <mrow> <mn>2</mn> <mo>+</mo> <mi>R</mi> </mrow> <mo>)</mo> <msub> <mi>T</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>T</mi> <mn>4</mn> </msub> <mo>)</mo> </mrow> </mrow>
减减速阶段T7
A7=0
<mrow> <msub> <mi>V</mi> <mn>7</mn> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msub> <mi>A</mi> <mi>x</mi> </msub> <msub> <mi>T</mi> <mn>1</mn> </msub> <mo>+</mo> <mfrac> <msub> <mi>J</mi> <mn>1</mn> </msub> <mrow> <mn>2</mn> <msup> <mi>R</mi> <mn>2</mn> </msup> </mrow> </mfrac> <msup> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mo>(</mo> <mrow> <mn>1</mn> <mo>+</mo> <mi>R</mi> </mrow> <mo>)</mo> <msub> <mi>T</mi> <mn>2</mn> </msub> <mo>-</mo> <mo>(</mo> <mrow> <mn>2</mn> <mo>+</mo> <mi>R</mi> </mrow> <mo>)</mo> <msub> <mi>T</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>T</mi> <mn>4</mn> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>-</mo> <mfrac> <msub> <mi>A</mi> <mi>x</mi> </msub> <mi>R</mi> </mfrac> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mo>(</mo> <mrow> <mn>1</mn> <mo>+</mo> <mi>R</mi> </mrow> <mo>)</mo> <msub> <mi>T</mi> <mn>2</mn> </msub> <mo>-</mo> <mo>(</mo> <mrow> <mn>2</mn> <mo>+</mo> <mi>R</mi> </mrow> <mo>)</mo> <msub> <mi>T</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>T</mi> <mn>4</mn> </msub> <mo>)</mo> </mrow> </mrow>
步骤三、由各阶段之间的连续性及边界条件确定出各阶段的参数
T7=RT1
T6=RT2
T5=RT3
J2=R2J1
步骤四、通过已知参数确定出短位移、中位移、长位移下的数学模型
在Vx和Ax限制条件不变的情况下,目标位移P发生变化时T1和T4可能等于零;为T1,T4恰好同时等于零时的目标位移,为恰好T4等于零时的目标位移;
<mrow> <msubsup> <mi>T</mi> <mn>1</mn> <mo>*</mo> </msubsup> <mo>=</mo> <mfrac> <msub> <mi>A</mi> <mi>x</mi> </msub> <msub> <mi>J</mi> <mn>1</mn> </msub> </mfrac> <mo>=</mo> <mfrac> <mrow> <msub> <mi>&amp;beta;V</mi> <mi>x</mi> </msub> </mrow> <msub> <mi>A</mi> <mi>x</mi> </msub> </mfrac> <mo>,</mo> <msubsup> <mi>T</mi> <mn>2</mn> <mo>*</mo> </msubsup> <mo>=</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mi>&amp;beta;</mi> <mo>)</mo> </mrow> <mfrac> <msub> <mi>V</mi> <mi>x</mi> </msub> <msub> <mi>A</mi> <mi>x</mi> </msub> </mfrac> </mrow>
<mrow> <msubsup> <mi>P</mi> <mi>s</mi> <mo>*</mo> </msubsup> <mo>=</mo> <mi>P</mi> <msub> <mo>|</mo> <mrow> <msub> <mi>T</mi> <mn>1</mn> </msub> <mo>=</mo> <msubsup> <mi>T</mi> <mn>1</mn> <mo>*</mo> </msubsup> </mrow> </msub> <mo>=</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mi>R</mi> <mo>)</mo> </mrow> <msup> <mi>&amp;beta;</mi> <mn>2</mn> </msup> <mfrac> <msubsup> <mi>V</mi> <mi>x</mi> <mn>2</mn> </msubsup> <msub> <mi>A</mi> <mi>x</mi> </msub> </mfrac> <mo>,</mo> <msubsup> <mi>P</mi> <mi>l</mi> <mo>*</mo> </msubsup> <mo>=</mo> <mi>P</mi> <msub> <mo>|</mo> <mrow> <msub> <mi>T</mi> <mn>2</mn> </msub> <mo>=</mo> <msubsup> <mi>T</mi> <mn>2</mn> <mo>*</mo> </msubsup> </mrow> </msub> <mo>=</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mi>R</mi> <mo>)</mo> </mrow> <mfrac> <mrow> <mn>1</mn> <mo>+</mo> <mi>&amp;beta;</mi> </mrow> <mn>2</mn> </mfrac> <mfrac> <msubsup> <mi>V</mi> <mi>x</mi> <mn>2</mn> </msubsup> <msub> <mi>A</mi> <mi>x</mi> </msub> </mfrac> </mrow>
时为短位移,T2=T4=0
<mrow> <msub> <mi>T</mi> <mn>1</mn> </msub> <mo>=</mo> <mroot> <mrow> <mfrac> <msubsup> <mi>T</mi> <mn>1</mn> <mo>*</mo> </msubsup> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mi>R</mi> <mo>)</mo> <msub> <mi>A</mi> <mi>x</mi> </msub> </mrow> </mfrac> <mi>P</mi> </mrow> <mn>3</mn> </mroot> </mrow>
时为中位移,
<mrow> <msub> <mi>T</mi> <mn>2</mn> </msub> <mo>=</mo> <mfrac> <mrow> <msqrt> <mrow> <mo>&amp;lsqb;</mo> <msub> <mi>A</mi> <mi>x</mi> </msub> <mrow> <mo>(</mo> <mi>R</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> <mo>&amp;lsqb;</mo> <msubsup> <mi>A</mi> <mi>x</mi> <mn>3</mn> </msubsup> <mrow> <mo>(</mo> <mi>R</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>+</mo> <mn>8</mn> <msub> <mi>PJ</mi> <mn>1</mn> </msub> <mo>&amp;rsqb;</mo> </mrow> </msqrt> <mo>-</mo> <mn>3</mn> <mrow> <mo>(</mo> <mi>R</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> <msubsup> <mi>A</mi> <mi>x</mi> <mn>2</mn> </msubsup> </mrow> <mrow> <mn>2</mn> <mrow> <mo>(</mo> <mi>R</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> <msub> <mi>J</mi> <mn>1</mn> </msub> <msub> <mi>A</mi> <mi>x</mi> </msub> </mrow> </mfrac> </mrow>
时为长位移,
<mrow> <msub> <mi>T</mi> <mn>4</mn> </msub> <mo>=</mo> <mfrac> <mrow> <mi>P</mi> <mo>-</mo> <msubsup> <mi>P</mi> <mi>l</mi> <mo>*</mo> </msubsup> </mrow> <msub> <mi>V</mi> <mi>x</mi> </msub> </mfrac> </mrow>
步骤五、通过改变非对称系数R便获得了不同目标位移下的非对称曲线;
步骤六、通过合理的选取非对称S曲线的采样点进行三次样条插值拟合出机械手各个关节的输入;
长位移的非对称S曲线的行程为:以冲击(加加速度)J1进行的加加速阶段(T1),匀加速阶段(T2),以冲击J1进行的减加速阶段(T3),匀速阶段(T4),以冲击J2进行的减加速阶段(T5),匀减速阶段(T6),以冲击J2进行的减减速阶段(T7);
中位移的非对称S曲线的行程为:以冲击(加加速度)J1进行的加加速阶段(T1),匀加速阶段(T2),以冲击J1进行的减加速阶段(T3),以冲击J2进行的减加速阶段(T5),匀减速阶段(T6),以冲击J2进行的减减速阶段(T7);
短位移的非对称S曲线的行程为:以冲击(加加速度)J1进行的加加速阶段(T1),以冲击J1进行的减加速阶段(T3),以冲击J2进行的减加速阶段(T5),以冲击J2进行的减减速阶段(T7);
依据非对称S曲线的运行特点得知加速行程和减速行程的冲击值为常数,匀速运行阶段的加速度为零,有如下等式的约束条件:
T7=RT1
T6=RT2
T5=RT3
J2=R2J1
T1=T3,T5=T7
因此减速行程的参数可以由加速行程和非对称参数R确定;评价指标和残余振动的振幅Ap、残余振动衰减时间Tp的关系如下式所示:
在目标位移发生变化时,匀速阶段T4,匀加速阶段T2,当目标位移为长位移时T2≠0,T4≠0,当目标位移为中位移时T2≠0,T4=0,当目标位移为短位移时T2=0,T4=0;
由非对称S曲线的对称性可知T1=T3,T5=T7,高速机械手减少残余振动的非对称S曲线的规划方法的模型为:
Find(T1,T2,T4,R)
Objective:
Subject to:
CN201711053996.2A 2017-10-31 2017-10-31 一种高速机械手减少残余振动的非对称s型轨迹规划方法 Active CN107825424B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711053996.2A CN107825424B (zh) 2017-10-31 2017-10-31 一种高速机械手减少残余振动的非对称s型轨迹规划方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711053996.2A CN107825424B (zh) 2017-10-31 2017-10-31 一种高速机械手减少残余振动的非对称s型轨迹规划方法

Publications (2)

Publication Number Publication Date
CN107825424A true CN107825424A (zh) 2018-03-23
CN107825424B CN107825424B (zh) 2020-05-26

Family

ID=61651360

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711053996.2A Active CN107825424B (zh) 2017-10-31 2017-10-31 一种高速机械手减少残余振动的非对称s型轨迹规划方法

Country Status (1)

Country Link
CN (1) CN107825424B (zh)

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109542134A (zh) * 2018-10-24 2019-03-29 广州霞光技研有限公司 一种解决机构间隙引起抖动的速度控制算法
CN110000794A (zh) * 2019-05-06 2019-07-12 江苏集萃智能制造技术研究所有限公司 一种基于协作机器人的截断式非对称速度规划方法
CN110398993A (zh) * 2018-04-24 2019-11-01 北京京东尚科信息技术有限公司 速度控制方法、装置及计算机可读存储介质
CN110455290A (zh) * 2019-07-17 2019-11-15 爱克斯维智能科技(苏州)有限公司 一种智能液压挖掘机的最优轨迹规划方法
CN110451408A (zh) * 2019-07-24 2019-11-15 江苏天宙检测有限公司 一种基于能耗优化的塔机运动规划方法
CN111015668A (zh) * 2019-12-27 2020-04-17 深圳市越疆科技有限公司 一种加加速度连续的速度规划方法、装置、控制器及机器人
CN111015674A (zh) * 2020-01-13 2020-04-17 上海景吾智能科技有限公司 一种机器人运动关节控制方法、系统、存储介质及机器人
CN112327954A (zh) * 2020-11-14 2021-02-05 大连理工大学 非对称s型速度曲线控制的直线电机高精定位方法
CN113110568A (zh) * 2021-03-13 2021-07-13 浙江御穹电子科技有限公司 一种无人机运动轨迹规划系统及方法
CN113246175A (zh) * 2021-05-07 2021-08-13 配天机器人技术有限公司 机器人末端震动分类方法、机器人及存储介质
NL2029652B1 (en) * 2021-11-05 2023-06-05 Flex Ind B V Handling device
CN116700151A (zh) * 2023-07-13 2023-09-05 哈尔滨工业大学 一种精密运动平台运动轨迹规划系统及其参数整定方法
CN116700150A (zh) * 2023-07-13 2023-09-05 哈尔滨工业大学 精密运动平台点位运动鲁棒轨迹规划系统及其规划方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103513575A (zh) * 2013-10-08 2014-01-15 广东工业大学 一种高速机构减小残余振动的s型运动曲线规划方法
CN106444635A (zh) * 2015-08-04 2017-02-22 深圳市雷赛智能控制股份有限公司 一种运动控制系统的非对称s曲线加减速控制方法和装置
CN108287527A (zh) * 2017-07-18 2018-07-17 山东大学 一种基于三角函数的改进s曲线加减速控制方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103513575A (zh) * 2013-10-08 2014-01-15 广东工业大学 一种高速机构减小残余振动的s型运动曲线规划方法
CN106444635A (zh) * 2015-08-04 2017-02-22 深圳市雷赛智能控制股份有限公司 一种运动控制系统的非对称s曲线加减速控制方法和装置
CN108287527A (zh) * 2017-07-18 2018-07-17 山东大学 一种基于三角函数的改进s曲线加减速控制方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
FENGSHAN ZOU等: "《Asymmetric S-curve Trajectory Planning for Robot Point-to-point Motion》", 《PROCEEDINGS OF THE 2009 IEEE》 *
KEUN-HO REW等: "《Using Asymmetric S-Curve Profile for Fast and Vibrationless Motion》", 《INTERNATIONAL CONFERENCE ON CONTROL, AUTOMATION AND SYSTEMS 2007》 *
吴军等: "《平面并联机构的残余振动控制》", 《机械工程学报》 *
吴明月等: "《基于最优S型曲线轨迹的硅片传输机器人末端振动抑制方法》", 《机器人》 *
潘海鸿等: "《一种适用于前瞻的高精度7段式非对称S曲线加减速控制》", 《机械科学与技术》 *

Cited By (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110398993A (zh) * 2018-04-24 2019-11-01 北京京东尚科信息技术有限公司 速度控制方法、装置及计算机可读存储介质
CN110398993B (zh) * 2018-04-24 2022-09-06 北京京东尚科信息技术有限公司 速度控制方法、装置及计算机可读存储介质
CN109542134A (zh) * 2018-10-24 2019-03-29 广州霞光技研有限公司 一种解决机构间隙引起抖动的速度控制算法
CN110000794A (zh) * 2019-05-06 2019-07-12 江苏集萃智能制造技术研究所有限公司 一种基于协作机器人的截断式非对称速度规划方法
CN110455290B (zh) * 2019-07-17 2020-12-08 爱克斯维智能科技(苏州)有限公司 一种智能液压挖掘机的最优轨迹规划方法
CN110455290A (zh) * 2019-07-17 2019-11-15 爱克斯维智能科技(苏州)有限公司 一种智能液压挖掘机的最优轨迹规划方法
CN110451408A (zh) * 2019-07-24 2019-11-15 江苏天宙检测有限公司 一种基于能耗优化的塔机运动规划方法
CN111015668A (zh) * 2019-12-27 2020-04-17 深圳市越疆科技有限公司 一种加加速度连续的速度规划方法、装置、控制器及机器人
CN111015674A (zh) * 2020-01-13 2020-04-17 上海景吾智能科技有限公司 一种机器人运动关节控制方法、系统、存储介质及机器人
CN111015674B (zh) * 2020-01-13 2021-05-11 上海景吾智能科技有限公司 一种机器人运动关节控制方法、系统、存储介质及机器人
CN112327954A (zh) * 2020-11-14 2021-02-05 大连理工大学 非对称s型速度曲线控制的直线电机高精定位方法
CN113110568A (zh) * 2021-03-13 2021-07-13 浙江御穹电子科技有限公司 一种无人机运动轨迹规划系统及方法
CN113110568B (zh) * 2021-03-13 2022-09-20 浙江御穹电子科技有限公司 一种无人机运动轨迹规划方法
CN113246175A (zh) * 2021-05-07 2021-08-13 配天机器人技术有限公司 机器人末端震动分类方法、机器人及存储介质
CN113246175B (zh) * 2021-05-07 2024-05-03 配天机器人技术有限公司 机器人末端震动分类方法、机器人及存储介质
NL2029652B1 (en) * 2021-11-05 2023-06-05 Flex Ind B V Handling device
CN116700151A (zh) * 2023-07-13 2023-09-05 哈尔滨工业大学 一种精密运动平台运动轨迹规划系统及其参数整定方法
CN116700150A (zh) * 2023-07-13 2023-09-05 哈尔滨工业大学 精密运动平台点位运动鲁棒轨迹规划系统及其规划方法
CN116700150B (zh) * 2023-07-13 2024-01-30 哈尔滨工业大学 精密运动平台点位运动鲁棒轨迹规划系统及其规划方法

Also Published As

Publication number Publication date
CN107825424B (zh) 2020-05-26

Similar Documents

Publication Publication Date Title
CN107825424B (zh) 一种高速机械手减少残余振动的非对称s型轨迹规划方法
CN111788040B (zh) 机器人的动力学参数辨识方法、机器人和存储装置
JP6499720B2 (ja) 機械学習装置、サーボ制御装置、サーボ制御システム、及び機械学習方法
CN103563245B (zh) 驱动机械的负荷特性推测装置
CN109471361B (zh) 一种异构多机器人系统的多目标跟踪控制方法
CN111176306B (zh) 一种水下航行器变增益自抗扰纵向控制方法
WO2022160894A1 (zh) 基于迭代学习的数控机床快速误差补偿控制系统及方法
CN112255913B (zh) 一种用于飞行拍摄效率最优的柔性速度规划方法
JP2018073097A (ja) 工具経路生成装置、工具経路生成方法及び工具経路生成プログラム
CN110703684B (zh) 一种端点速度不限的轨迹规划方法及装置
CN116774648A (zh) 速度规划方法、装置、机床控制系统和存储介质
CN109164699A (zh) 一种贴片机机头运行位置精度控制方法
CN107544247B (zh) 一种抑制机械运动结构振动的方法及系统
CN115556103A (zh) 一种基于动力学前馈的分数阶阻抗控制器及其设计方法
CN109324503A (zh) 基于鲁棒积分的多层神经网络电机系统控制方法
Śleszycki et al. Parameter Identification of the Two-Mass System with the help of Multi-layer Estimator
CN107918356A (zh) 用于数控机床轴联动加工过程中的过象限摩擦补偿的方法
Ebrahimi Regulated model-based and non-model-based sliding mode control of a MEMS vibratory gyroscope
CN109129470A (zh) 机器人直线路径的规划方法及系统
CN111046510B (zh) 一种基于轨迹分段优化的柔性机械臂的振动抑制方法
CN109309468B (zh) 一种直线电机调速系统及方法
CN116107263B (zh) 一种消除末端器件残余振动的方法、装置、工控机及介质
CN113467234B (zh) 一种确定飞机极限环振荡幅值与频率的方法及装置
CN114237152B (zh) 一种用于激光切割的柔性速度规划及位移补偿方法
TWI777395B (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