CN106483927A - 一种点位运动s曲线生成方法 - Google Patents

一种点位运动s曲线生成方法 Download PDF

Info

Publication number
CN106483927A
CN106483927A CN201611170134.3A CN201611170134A CN106483927A CN 106483927 A CN106483927 A CN 106483927A CN 201611170134 A CN201611170134 A CN 201611170134A CN 106483927 A CN106483927 A CN 106483927A
Authority
CN
China
Prior art keywords
time
function
max
curve
acceleration
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
CN201611170134.3A
Other languages
English (en)
Other versions
CN106483927B (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.)
Victor (dongguan) Intelligent Control Co Ltd
Original Assignee
Victor (dongguan) Intelligent Control Co Ltd
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 Victor (dongguan) Intelligent Control Co Ltd filed Critical Victor (dongguan) Intelligent Control Co Ltd
Priority to CN201611170134.3A priority Critical patent/CN106483927B/zh
Publication of CN106483927A publication Critical patent/CN106483927A/zh
Application granted granted Critical
Publication of CN106483927B publication Critical patent/CN106483927B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B19/00Programme-control systems
    • G05B19/02Programme-control systems electric
    • G05B19/18Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form
    • G05B19/19Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form characterised by positioning or contouring control systems, e.g. to control position from one programmed point to another or to control movement along a programmed continuous path
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00Program-control systems
    • G05B2219/30Nc systems
    • G05B2219/36Nc in input of data, input key till input tape
    • G05B2219/36063During machining, compare simulated with detected profile, correct, modify program

Abstract

本发明属于运动控制领域,具体是一种点位运动S曲线生成方法,包括以下步骤:S1:将S曲线位置‑时间关系通过无量纲化拟合成5次方多项式的函数,并将位置‑时间函数进行微分,得到速度‑时间函数;对速度‑时间函数进行微分,得到了加速度‑时间函数;S2:要求满足条件一以及满足条件二,并通过条件一和条件二求得6个函数系数;S3:通过无量纲值和真实值之间的转换关系,求出真实值位置相对于时间的函数,以及求出真实值速度相对于时间的函数。本发明使S曲线生成方法简单,速度函数为4次方函数,加速度为3次方函数,表明曲线光滑、柔和;可以进行设定值由加减速距离转换为加减速时间,对加减速距离要求严格的场合将得到精确的加减速距离控制。

Description

一种点位运动S曲线生成方法
技术领域
本发明属于运动控制领域,具体是一种点位运动S曲线生成方法。
背景技术
运动控制中的点位运动就是电机从一个位置走到另一个位置的运动。速度曲线的生成是点位运动的关键。相较于梯形曲线,S型曲线具有加速度连续,速度曲线光滑、均匀,运动平稳、无冲击等优势,能够获得较理想的运动控制效果。
传统PLC是靠发送脉冲控制伺服或步进系统实现点位运动,其速度曲线分为梯形曲线和多段速曲线。
运动控制器靠模拟量或发送脉冲控制伺服或步进系统实现点位运动,其速度曲线分为梯形曲线和七段速S型曲线“。如专利“运动控制系统S曲线加减速的实现方法”所述,七段速S型曲线加减速的完整过程分为7个运动段,它们是加加速段、匀加速段、减加速段、匀速段、加减速段、匀减速段与减减速段。为生成该曲线需设置加速时间、加加速时间、加减速时间、减加速时间、减减速时间、减速时间、最高运行速度、点位运行距离。上述S型曲线的计算方法较为复杂繁琐,在特定场合才需要用到。
发明内容
为了克服现有技术中存在的缺点和不足,本发明的目的在于提供一种点位运动S曲线生成方法,将使得曲线计算简单化;此算法经过转变,也可以通过设置加速距离和减速距离来生成S曲线,来精确的控制加减速距离。
本发明的技术方案如下:一种点位运动S曲线生成方法,包括以下步骤:
S1:将S曲线的位置-时间关系通过无量纲化拟合成5次方多项式的函数:Y=FY(T)=D0+D1T+D2T2+D3T3+D4T4+D5T5,并将位置-时间函数进行微分,得到了速度-时间函数:V=FV(T)=D1+2D2T+3D3T2+4D4T3+5D5T4;对速度-时间函数进行微分,得到了加速度-时间函数:A=FA(T)=2D2+6D3T+12D4T2+20D5T3
S2:要求满足条件一:T=0时,Y=0,V=0,A=0,以及满足条件二:T=1时,Y=1,V=0,A=0,并通过条件一和条件二求得D0、D1、D2、D3、D4、D56个函数系数;
S3:通过无量纲值和真实值之间的转换关系,求出真实值位置相对于时间的函数,以及求出真实值速度相对于时间的函数。
进一步的,还包括以下步骤:
S1:设定预值加速时间t1、减速时间t3、运行距离s和最高运行速度vset
S2:计算出加速段距离S1:S1=16t1vmax/30以及减速段距离S3:S3=16t3vmax/30,判断s与S1+S3的大小关系:
若s>S1+S3,则为第一种曲线,第一种曲线包括加速段、匀速段和减速段,并且vmax=vset;计算出匀速距离S2:S2=s-S1-S3和匀速过程运行时间t2:t2=S2/vmax
若s<=S1+S3,则为第二种曲线,第二种曲线包括加速段和减速段,重新计算最高运行速度vmax:vmax=30S11/16t1=30S13/16t3=30s/16(t1+t3)、加速距离S11:S11=S1s/(S1+S3)、减速距离S13:S13=S3s/(S1+S3);
S3:通过坐标偏移得到真实值位置相对于时间为变量的分段函数以及真实值速度相对于时间为变量的分段函数。
进一步的,还包括以下步骤:
S1:设定预值加速距离S1、减速距离S3、运行距离s和最高运行速度vset
S2:直接判断s与S1+S3的大小关系:
若s>S1+S3,则为第一种曲线,即vmax=vset,通过S1=16t1vmax/30计算出加速时间t1和通过S3=16t3vmax/30计算出减速时间t3,并且计算出匀速距离S2:S2=s-S1-S3,和计算出匀速运行时间t2:t2=S2/vmax
若s<=S1+S3,则为第二种曲线,则通过:S1=16t1vmax/30、S3=16t3vmax/30、S11=S1s/(S1+S3)、S13=S3s/(S1+S3)以及vmax=30S11/16t1=30S13/16t3=30s/16(t1+t3)重新计算最高运行速度vmax、加速时间t1和减速时间t3
S3:通过坐标偏移得到真实值位置相对于时间为变量的分段函数以及真实值速度相对于时间为变量的分段函数。
本发明的有益效果:本发明提出的方法为一种点位运动S型速度曲线的生成方法,使得S曲线生成方法简单,速度函数为4次方函数,加速度为3次方函数,表明曲线光滑、柔和;并且可以进行设定值由加减速距离转换为加减速时间,对加减速距离要求严格的场合将得到精确的加减速距离控制。
附图说明
图1是本发明的实现流程图。
图2是本发明的无量纲化的位置-时间、速度-时间、加速度-时间曲线。
图3是本发明的带有匀速段的速度-时间曲线。
图4是本发明的没有匀速段的速度-时间曲线。
具体实施方式
下面结合附图1与具体实施方式对本技术方案进一步说明如下:
第一步:
设定预值为加速时间t1、减速时间t3、运行距离s和最高运行速度vset
将S曲线的位置-时间关系通过无量纲化拟合成5次方多项式的函数:
Y=FY(T)=D0+D1T+D2T2+D3T3+D4T4+D5T5 (1)
对其微分得到速度-时间函数、对速度-时间函数进行微分得到加速度-时间函数:
V=FV(T)=D1+2D2T+3D3T2+4D4T3+5D5T4 (2)
A=FA(T)=2D2+6D3T+12D4T2+20D5T3 (3)
要求满足条件一:T=0时,Y=0,V=0,A=0,以及满足条件二T=1时,Y=1,V=0,A=0。
带入可得D0、D1、D2、D3、D4、D5,将求得的系数带入式(1)、式(2)、式(3)可得到3个系数已知的函数方程,如图2所示。
对上述3个函数进行分析,T∈[0,0.5],A>0,为加速段,T∈[0.5,1],A<0,为减速段,T=0.5时,速度到达T∈[0,1]的最高点。
第二步:
设无量纲值和真实值之间的转换关系为Y=y/y0,T=t/t0
求出真实值位置相对于时间t的函数、真实值速度相对于时间t的函数表达如下:
y=y0FY(t/t0) (4)
v=v0FV(t/t0) (5)
根据第一步的分析,公式(5)要满足当t=t0/2时v=vmax,并且公式(5)是公式(4)的微分那么满足v0=y0/t0
根据以上两个条件,并且由式(2)可得T=0.5时,V到达最高点30/16,可得以下两个结论:
v0=16vmax/30, (6)
y0=16t0vmax/30 (7)
因为由公式(2)、(3)可以分析出无量纲化公式的加速段取T∈[0,0.5],减速段取T∈[0.5,1]。由此可知,在无量纲值和真实值转换后在加速段取:
t0=2t1 (8)
在减速段取:
t0=2t3。 (9)
第三步:
由式(1)、(4)、(8),并将T∈[0,0.5]作为加速段可知:
加速段距离S1=16t1vmax/30 (10)
由式(1)、(4)、(9),并将T∈[0.5,1]做为减速段可知:
减速距离S3=16t3vmax/30 (11)
判断s与S1+S3的关系。
如果s>S1+S3,如图3所示,则为第一种曲线,即完整曲线,第一种曲线包括加速段、匀速段和减速段并且:
vmax=vset (12)
匀速距离S2=s-S1-S3 (13)
匀速过程运行时间t2=S2/vmax (14)
如果s<=S1+S3,如图4所示,则为第二种曲线,第二种曲线包括加速段和减速段,第二种曲线没有匀速段,重新计算最高运行速度vmax、加速距离t1、减速距离t3
加速距离S11=S1s/(S1+S3) (15)
减速距离S13=S3s/(S1+S3) (16)
最高运行速度vmax=30S11/16t1=30S13/16t3=30s/16(t1+t3) (17)
第四歩:
根据以上推导对公式(4)、公式(5)进行坐标偏移得到真实值位置相对于时间为变量的分段函数以及真实值速度相对于时间为变量的分段函数。其中,式(18)、(19)为第一种曲线,即有匀速段时候的真实值位置相对于时间为变量的分段函数以及真实值速度相对于时间为变量的分段函数,式(20)、(21)为第二种曲线,即无匀速段时候的真实值位置相对于时间为变量的分段函数以及真实值速度相对于时间为变量的分段函数。
实施例二:
第一步:
设定预值为加速度距离s1、减速度距离s3运行距离s和最高运行速度vset
将S曲线的位置-时间关系通过无量纲化拟合成5次方多项式的函数:
Y=FY(T)=D0+D1T+D2T2+D3T3+D4T4+D5T5 (1)
对其微分得到速度-时间函数、对速度-时间函数进行微分得到加速度-时间函数:
V=FV(T)=D1+2D2T+3D3T2+4D4T3+5D5T4 (2)
A=FA(T)=2D2+6D3T+12D4T2+20D5T3 (3)
要求满足条件一:T=0时,Y=0,V=0,A=0,以及满足条件二T=1时,Y=1,V=0,A=0。
带入可得D0、D1、D2、D3、D4、D5,将求得的系数带入式(1)、式(2)、式(3)可得到3个系数已知的函数方程,如图2所示。
对上述3个函数进行分析,T∈[0,0.5],A>0,为加速段,T∈[0.5,1],A<0,为减速段,T=0.5时,速度到达T∈[0,1]的最高点。
第二步:
设无量纲值和真实值之间的转换关系为Y=y/y0,T=t/t0
求出真实值位置相对于时间t的函数、真实值速度相对于时间t的函数表达如下:
y=y0FY(t/t0) (4)
v=v0FV(t/t0) (5)
根据第一步的分析,公式(5)要满足当t=t0/2时v=vmax,并且公式(5)是公式(4)的微分那么满足v0=y0/t0
根据以上两个条件,并且由式(2)可得T=0.5时,V到达最高点30/16,可得以下两个结论:
v0=16vmax/30, (6)
y0=16t0vmax/30 (7)
因为由式(2)、(3)可以分析出无量纲化公式的加速段取T∈[0,0.5],减速段取T∈[0.5,1]。由此可知,在无量纲值和真实值转换后在加速段取:
t0=2t1 (8)
在减速段取:
t0=2t3。 (9)
第三步:
由式(1)、(4)、(8),并将T∈[0,0.5]作为加速段可知:
加速段距离S1=16t1vmax/30 (10)
由式(1)、(4)、(9),并将T∈[0.5,1]做为减速段可知:
减速距离S3=16t3vmax/30 (11)
直接判断s与s1+s3的关系。
若s>S1+S3,则为第一种曲线,根据式(10)、计算出加速时间t1,根据式(11)算出减速时间t3
第一种曲线包括加速段、匀速段和减速段并且:
vmax=vset (12)
匀速距离S2=s-S1-S3 (13)
匀速过程运行时间t2=S2/vmax (14)
若s<=S1+S3,则为第二种曲线,第二种曲线包括加速段和减速段,第二种曲线没有匀速段,重新计算最高运行速度vmax、加速距离t1、减速距离t3
加速距离S11=S1s/(S1+S3) (15)
减速距离S13=S3s/(S1+S3) (16)
最高运行速度vmax=30S11/16t1=30S13/16t3=30s/16(t1+t3) (17)
根据公式(10)、(11)、(15)、(16)以及(17)重新计算最高运行速度vmax、加速时间t1和减速时间t3
第四歩:
根据式(4)、(5)进行坐标偏移得到真实值位置相对于时间为变量的分段函数以及真实值速度相对于时间为变量的分段函数。其中,式(18)、(19)为第一种曲线,即有匀速段时候的真实值位置相对于时间为变量的分段函数以及真实值速度相对于时间为变量的分段函数,式(20)、(21)为第二种曲线,即无匀速段时候的真实值位置相对于时间为变量的分段函数以及真实值速度相对于时间为变量的分段函数。
综上所述,本发明使得S曲线生成方法简单,速度函数为4次方函数,加速度为3次方函数,表明曲线光滑、柔和,并且根据实施例二可以进行设定值由加减速距离转换为加减速时间,对加减速距离要求严格的场合将得到精确的加减速距离控制。
最后应当说明的是,以上实施例仅用以说明本发明的技术方案,而非对本发明保护范围的限制,尽管参照较佳实施例对本发明作了详细地说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的实质和范围。

Claims (3)

1.一种点位运动S曲线生成方法,其特征在于:包括以下步骤:
S1:将S曲线的位置-时间关系通过无量纲化拟合成5次方多项式的函数:Y=FY(T)=D0+D1T+D2T2+D3T3+D4T4+D5T5,并将位置-时间函数进行微分,得到了速度-时间函数:V=FV(T)=D1+2D2T+3D3T2+4D4T3+5D5T4;对速度-时间函数进行微分,得到了加速度-时间函数:A=FA(T)=2D2+6D3T+12D4T2+20D5T3
S2:要求满足条件一:T=0时,Y=0,V=0,A=0,以及满足条件二:T=1时,Y=1,V=0,A=0,并通过条件一和条件二求得D0、D1、D2、D3、D4、D56个函数系数;
S3:通过无量纲值和真实值之间的转换关系,求出真实值位置相对于时间的函数,以及求出真实值速度相对于时间的函数。
2.根据权利要求1所述的一种点位运动S曲线生成方法,其特征在于:还包括以下步骤:
S1:设定预值加速时间t1、减速时间t3、运行距离s和最高运行速度vset
S2:计算出加速段距离S1:S1=16t1vmax/30以及减速段距离S3:S3=16t3vmax/30,判断s与S1+S3的大小关系:
若s>S1+S3,则为第一种曲线,第一种曲线包括加速段、匀速段和减速段,并且vmax=vset;计算出匀速距离S2:S2=s-S1-S3和匀速过程运行时间t2:t2=S2/vmax
若s<=S1+S3,则为第二种曲线,第二种曲线包括加速段和减速段,重新计算最高运行速度vmax:vmax=30S11/16t1=30S13/16t3=30s/16(t1+t3)、加速距离S11:S11=S1s/(S1+S3)、减速距离S13:S13=S3s/(S1+S3);
S3:通过坐标偏移得到真实值位置相对于时间为变量的分段函数以及真实值速度相对于时间为变量的分段函数。
3.根据权利要求1所述的一种点位运动S曲线生成方法,其特征在于:还包括以下步骤:
S1:设定预值加速距离S1、减速距离S3、运行距离s和最高运行速度vset
S2:直接判断s与S1+S3的大小关系:
若s>S1+S3,则为第一种曲线,即vmax=vset,通过S1=16t1vmax/30计算出加速时间t1和通过S3=16t3vmax/30计算出减速时间t3,并且计算出匀速距离S2:S2=s-S1-S3,和计算出匀速运行时间t2:t2=S2/vmax
若s<=S1+S3,则为第二种曲线,则通过:S1=16t1vmax/30、S3=16t3vmax/30、S11=S1s/(S1+S3)、S13=S3s/(S1+S3)以及vmax=30S11/16t1=30S13/16t3=30s/16(t1+t3)重新计算最高运行速度vmax、加速时间t1和减速时间t3
S3:通过坐标偏移得到真实值位置相对于时间为变量的分段函数以及真实值速度相对于时间为变量的分段函数。
CN201611170134.3A 2016-12-16 2016-12-16 一种点位运动s曲线生成方法 Active CN106483927B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611170134.3A CN106483927B (zh) 2016-12-16 2016-12-16 一种点位运动s曲线生成方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611170134.3A CN106483927B (zh) 2016-12-16 2016-12-16 一种点位运动s曲线生成方法

Publications (2)

Publication Number Publication Date
CN106483927A true CN106483927A (zh) 2017-03-08
CN106483927B CN106483927B (zh) 2018-08-31

Family

ID=58285416

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611170134.3A Active CN106483927B (zh) 2016-12-16 2016-12-16 一种点位运动s曲线生成方法

Country Status (1)

Country Link
CN (1) CN106483927B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109240215A (zh) * 2018-11-22 2019-01-18 山东易码智能科技股份有限公司 一种基于改进型s曲线加减速的pvt控制方法
CN110825025A (zh) * 2019-10-24 2020-02-21 威科达(东莞)智能控制有限公司 一种用于瓦楞纸前缘送纸的免编程电子凸轮曲线生成方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101770220A (zh) * 2008-12-30 2010-07-07 台达电子工业股份有限公司 数值控制系统的命令回算计算方法
CN103279069A (zh) * 2013-03-25 2013-09-04 昆山天大精益数控科技发展有限公司 用于数控机床的三角函数二阶连续可导加减速算法
CN105652794A (zh) * 2016-02-29 2016-06-08 东北大学 基于focas的数控机床数据采集和管理系统及其方法
CN105700466A (zh) * 2015-11-23 2016-06-22 上海交通大学 高速数控加工轨迹的曲率光顺方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101770220A (zh) * 2008-12-30 2010-07-07 台达电子工业股份有限公司 数值控制系统的命令回算计算方法
CN103279069A (zh) * 2013-03-25 2013-09-04 昆山天大精益数控科技发展有限公司 用于数控机床的三角函数二阶连续可导加减速算法
CN105700466A (zh) * 2015-11-23 2016-06-22 上海交通大学 高速数控加工轨迹的曲率光顺方法
CN105652794A (zh) * 2016-02-29 2016-06-08 东北大学 基于focas的数控机床数据采集和管理系统及其方法

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109240215A (zh) * 2018-11-22 2019-01-18 山东易码智能科技股份有限公司 一种基于改进型s曲线加减速的pvt控制方法
CN110825025A (zh) * 2019-10-24 2020-02-21 威科达(东莞)智能控制有限公司 一种用于瓦楞纸前缘送纸的免编程电子凸轮曲线生成方法
CN110825025B (zh) * 2019-10-24 2022-06-17 威科达(东莞)智能控制有限公司 一种用于瓦楞纸前缘送纸的免编程电子凸轮曲线生成方法

Also Published As

Publication number Publication date
CN106483927B (zh) 2018-08-31

Similar Documents

Publication Publication Date Title
CN101833306B (zh) 一种基于前瞻-滤波技术的多程序段连续加减速控制方法
CN103941647B (zh) 用于嵌入式数控设备的柔性加减速控制插补方法
CN108806252A (zh) 一种高速公路混合交通流协同优化控制方法
CN104570735B (zh) 一种码垛机器人弧形运动轨迹算法
CN104635494B (zh) 一种基于遗传算法优化的车辆追尾碰撞模糊控制方法
CN106483927A (zh) 一种点位运动s曲线生成方法
CN106406219A (zh) 一种用于横切的免编程电子凸轮曲线生成方法
CN108163018A (zh) 一种基于固定距离生成停车策略的列车精确停车方法
CN106950923A (zh) 一种定插补周期的速度规划方法
CN103466452B (zh) 一种起重机运行控制方法及起重机运行控制系统
CN103676787B (zh) 一种用于运动控制系统的圆心模式空间圆弧插补方法
CN106843207A (zh) 穿梭车速度规划控制方法、装置及穿梭车
CN107844058A (zh) 一种运动曲线离散动态规划方法
CN102756954B (zh) 一种电梯速度补偿方法
CN105882664B (zh) 索道车厢的速度控制方法
CN103312248B (zh) 一种基于dsp的直线加减速拐点误差补偿方法
CN110398993B (zh) 速度控制方法、装置及计算机可读存储介质
CN103246236B (zh) 一种考虑小直线段和连接点速度的数控系统速度控制方法
CN111596543B (zh) 一种快速规划6阶s型运动曲线的方法
CN105270948A (zh) 电梯的轿厢内气压控制装置
CN107826577A (zh) 一种自动码放喷涂件的智能货架
CN102880118B (zh) 基于插补精度和加速度限制的变插补周期曲线插补方法
CN205166193U (zh) 一种激光切割机双工作台交换机构
CN203409194U (zh) 制芯机移动小车的驱动控制装置
CN107368041A (zh) 一种基于移动平均算法的s曲线加减速控制方法

Legal Events

Date Code Title Description
C06 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