CN104950821B - 一种基于分数阶微积分的数控***速度规划方法 - Google Patents
一种基于分数阶微积分的数控***速度规划方法 Download PDFInfo
- Publication number
- CN104950821B CN104950821B CN201510358118.6A CN201510358118A CN104950821B CN 104950821 B CN104950821 B CN 104950821B CN 201510358118 A CN201510358118 A CN 201510358118A CN 104950821 B CN104950821 B CN 104950821B
- Authority
- CN
- China
- Prior art keywords
- speed
- alpha
- controlled object
- fractional
- control system
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 47
- 230000001133 acceleration Effects 0.000 claims abstract description 56
- 238000006073 displacement reaction Methods 0.000 claims abstract description 15
- 238000009795 derivation Methods 0.000 claims description 10
- 230000010354 integration Effects 0.000 claims description 9
- 239000000203 mixture Substances 0.000 claims description 3
- 230000006698 induction Effects 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 abstract description 2
- 230000035939 shock Effects 0.000 abstract 1
- 230000006870 function Effects 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000003754 machining Methods 0.000 description 1
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
- G05B19/00—Programme-control systems
- G05B19/02—Programme-control systems electric
- G05B19/18—Numerical 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/416—Numerical 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 control of velocity, acceleration or deceleration
-
- 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
- G05B2219/00—Program-control systems
- G05B2219/30—Nc systems
- G05B2219/43—Speed, acceleration, deceleration control ADC
- G05B2219/43065—Limitation of jerk
Landscapes
- Engineering & Computer Science (AREA)
- Human Computer Interaction (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
- Numerical Control (AREA)
Abstract
一种基于分数阶微积分的数控***速度规划方法,本发明涉及数控***速度规划方法。本发明是要解决现有的速度规划方法导致机床产生剧烈振动、计算量巨大,编程复杂以及出现加速度突变,产生刚性冲击的问题。本发明是通过1、利用连续函数y=f(t),推导分数阶导数;2、求得t1、t2和t3对应的名义加速度分别为a1、a2和a3;3、得到所需的加速段速度规划曲线;4、计算总位移增量;5、根据公式(3)求得位移为D(m);6、被控对象未达到最大速度Vmax则从减速点第i步开始减速直到速度为0;7、从减速点第j步开始减速到速度为0;8、被控对象按照类S型曲线进行减速到速度为0等步骤实现的。本发明应用于数控***速度规划领域。
Description
技术领域
本发明涉及数控***速度规划方法,特别涉及一种基于分数阶微积分的数控***速度规划方法。
背景技术
速度规划是数控***的核心技术。常用的速度规划方法有梯形速度规划方法,S型曲线速度规划方法,指数型速度规划方法等。梯形速度规划方法虽然计算量小,编程简单,但是在加减速阶段存在加速度突变的现象,导致机床产生剧烈振动,不适合于高速高精度加工。常用S型曲线速度规划方法通过限制加加速度来控制加速度的突变现象,然而传统的S形曲线速度规划方法利用多项式表示法将整个速度规划分为7个阶段,然后在每个阶段内进行讨论,在实现过程中不仅要对相邻阶段之间要进行边界判断,还要判断在实际运动中有几个阶段存在,比如运动距离很小(几毫米),那么就有匀速等一个或几个阶段不存在,所以计算量巨大,编程复杂。同样的指数型速度规划方法,在加减速切换的瞬间会出现加速度突变,产生刚性冲击,不适合于高速高精度场合。
发明内容
本发明的目的是为了解决现有的速度规划方法导致机床产生剧烈振动、计算量巨大,编程复杂以及出现加速度突变,产生刚性冲击的问题而提出的一种基于分数阶微积分的数控***速度规划方法。
上述的发明目的是通过以下技术方案实现的:
步骤一、根据被控对象运动时间t,利用连续函数y=f(t),推导分数阶导数:
其中,h为时间步长,j为连加运算∑的索引,分数阶积分的次数为α;c为微积分起始时间,通常c=0;当α>0时为分数阶微分运算,α<0时为分数阶积分运算;
步骤二、将被控对象加速段运动时间T分为三段分别为t1、t2和t3,并且t1、t2和t3分别对应的名义加速度分别为a1、a2和a3;其中,a1、a2和a3的约束是a2>a1>0,a3<0;
步骤三、令α<0,利用公式(2)对名义加速度a1,a2,a3进行分数阶积分,即得到所需的加速段速度规划曲线,此曲线为类S型曲线;
步骤四、被控对象第i个运动周期的速度为v(i),运动周期即时间步长为h;位移增量为d(i)则到第N周期结束,即总位移增量;
步骤五、t1时刻达到最大速度Vmax,第N周期走过的步数为m=t1/h,利用m=N根据公式(3)求得位移为D(m);
步骤六、被控对象按照类S型曲线进行运动,若行程s<D(m)时,被控对象未达到最大速度Vmax则从减速点第i步开始减速直到速度为0,满足s/2=D(i);
步骤七、被控对象按照类S型曲线进行运动;若行程s>D(m)时,被控对象达到了最大速度Vmax后,经过匀速Vmax过程后再从减速点第j步开始减速到速度为0,减速点满足s-D(j)=D(m);
步骤八、若行程s=D(m)时,被控对象达到最大速度Vmax后,被控对象按照类S型曲线进行减速到速度为0;即完成了一种基于分数阶微积分的数控***速度规划方法。
发明效果
针对现有技术中存在的不足之处,本发明提供一种数控***速度曲线规划的新方法,实现被控对象的平稳加减速。所采用的技术为分数阶微积分技术。本发明采用的分数阶微积分是对整数阶微积分的广义推广。它的研究几乎同时于传统的牛顿微积分体系,历史上许多大数学家都对此进行过相关研究。人们从各自领域出发曾经给出了多种定义,其中最适合我们本发明研究的是Grünwald-Letnikov(G-L)定义。
本发明采用的分数阶微积分具有记忆性和非局部性,能够很好地对数据进行平滑,做到精细化。采用分数阶微积分技术对加速度进行积分,可以得到平滑的加减速速度曲线,实现平稳速度规划的要求。如图3c所示在长行程下的加减速曲线,拥有完整的加速阶段,匀速阶段和减速阶段。
附图说明
图1为具体实施方式一提出的名义加速度示意图;
图2为具体实施方式一提出的速度规划曲线示意图;
图3a为具体实施方式一提出的不同程下的加减速曲线示意图;
图3b为具体实施方式一提出的不同程下的加减速曲线示意图;
图3c为具体实施方式一提出的不同程下的加减速曲线示意图。
具体实施方式
具体实施方式一:本实施方式的一种基于分数阶微积分的数控***速度规划方法,具体是按照以下步骤制备的:
一种基于分数阶微积分的数控***速度规划方法利用分数阶微积分的记忆特性,对预定义的名义加速度进行一次分数阶积分,即可得到平滑的加减速速度规划曲线
步骤一、根据被控对象运动时间t,利用连续函数y=f(t),推导分数阶导数:
其中,h为时间步长,j为连加运算∑的索引,分数阶积分的次数为α;c为微积分起始时间,通常c=0;当α>0时为分数阶微分运算,α<0时为分数阶积分运算;本发明即是利用公式(2)对被积函数进行分数阶积分得到期望的加减速曲线的;
步骤二、将被控对象加速段运动时间T分为三段分别为t1、t2和t3,并且t1、t2和t3分别对应的名义加速度分别为a1、a2和a3;其中,a1、a2和a3的约束是a2>a1>0,a3<0;
步骤三、令α<0,利用公式(2)对名义加速度a1,a2,a3进行分数阶积分,即得到所需的加速段速度规划曲线,此曲线为类S型曲线,最终整个过程的加速段速度规划曲线见图2;
步骤四、速度规划共分为三个阶段:加速阶段、匀速阶段和减速阶段;减速阶段与加速阶段是镜像的逆过程,而匀速阶段为以最大速度Vmax匀速运行,直到到达减速点;如此,整个速度规划即是对加速阶段的规划,下面介绍加速阶段速度规划方法;加速阶段速度规划有两个性能指标:一是加速时间T,二是最大速度Vmax,这两个参数作为用户的输入,是已知参量;然后是定义加速阶段的名义加速度,如图1所示;实际运动过程中,需要根据行程的大小,被控对象第i个运动周期的速度为v(i),运动周期即时间步长为h;位移增量为d(i)(位移增量就是每一步的走的位移,所以总位移增量为位移增量之和),则到第N周期结束,即总位移增量;
步骤五、图2所示,t1时刻达到最大速度Vmax,第N周期走过的步数为m=t1/h,利用m=N根据公式(3)求得位移为D(m);
步骤六、被控对象按照类S型曲线进行运动,若行程s<D(m)时,被控对象未达到最大速度Vmax则从减速点第i步开始减速直到速度为0,满足s/2=D(i);如图3a所示为短行程下的加减速曲线,此时没有达到最大速度就开始减速;
步骤七、被控对象按照类S型曲线进行运动;若行程s>D(m)时,被控对象达到了最大速度Vmax后,经过匀速Vmax过程后再从减速点第j步开始减速到速度为0,减速点满足s-D(j)=D(m)(s-D(j)=D(m)根据速度规划曲线加速段与减速段是对称镜像的);如图3c为长行程下的加减速曲线,拥有完整的加速阶段,匀速阶段和减速阶段;
步骤八、若行程s=D(m)时,被控对象达到最大速度Vmax后,被控对象按照类S型曲线进行减速到速度为0如图3b为刚好达到最大速度,但是没有匀速阶段,直接进入减速阶段;即完成了一种基于分数阶微积分的数控***速度规划方法。
本实施方式效果:
针对现有技术中存在的不足之处,本实施方式提供一种数控***速度曲线规划的新方法,实现被控对象的平稳加减速。所采用的技术为分数阶微积分技术。本实施方式采用的分数阶微积分是对整数阶微积分的广义推广。它的研究几乎同时于传统的牛顿微积分体系,历史上许多大数学家都对此进行过相关研究。人们从各自领域出发曾经给出了多种定义,其中最适合我们本实施方式研究的是Grünwald-Letnikov(G-L)定义。
本实施方式采用的分数阶微积分具有记忆性和非局部性,能够很好地对数据进行平滑,做到精细化。采用分数阶微积分技术对加速度进行积分,可以得到平滑的加减速速度曲线,实现平稳速度规划的要求。如图3c所示在长行程下的加减速曲线,拥有完整的加速阶段,匀速阶段和减速阶段。
具体实施方式二:本实施方式与具体实施方式一不同的是:步骤一中根据被控对象运动时间t,利用连续函数y=f(t),推导分数阶导数具体过程为:
(1)对于连续函数y=f(t),根据经典的求导公式,函数y=f(t)的一阶求导公式为:
其物理意义可以理解为上述函数的变化率;
(2)根据步骤(1)得到的f'(t)求导为二阶导数具体为:
其物理意义可以理解为曲线的曲率或者某个物理量的加速度;
(3)根据步骤(2)得到的f″(t)求导为三阶导数具体为:
(4)利用数学归纳法,根据步骤(1)~(3)的推导过程推导出f(t)的n阶导数具有以下形式:
其中,
由此,推论到分数阶导数如下:
其它步骤及参数与具体实施方式一相同。
具体实施方式三:本实施方式与具体实施方式一或二不同的是:步骤三中令α<0,步骤二中定义的名义加速度的数学表达式为:
利用公式(2)对名义加速度a(t)进行α阶分数阶积分,即可得到加速段的速度规划曲线。其它步骤及参数与具体实施方式一或二相同。
具体实施方式四:本实施方式与具体实施方式一至三之一不同的是:步骤三中α设定为-1~-2之间调整,以得到满足需要的加速段速度规划曲线。其它步骤及参数与具体实施方式一至三之一相同。
Claims (4)
1.一种基于分数阶微积分的数控***速度规划方法,其特征在于一种基于分数阶微积分的数控***速度规划方法,具体是按照以下步骤进行的:
步骤一、根据被控对象运动时间t,利用连续函数y=f(t),推导分数阶导数:
其中,h为时间步长,j为连加运算∑的索引,分数阶积分的次数为α;c为微积分起始时间,通常c=0;当α>0时为分数阶微分运算,α<0时为分数阶积分运算;
步骤二、将被控对象加速段运动时间T分为三段分别为t1、t2和t3,并且t1、t2和t3分别对应的名义加速度分别为a1、a2和a3;其中,a1、a2和a3的约束是a2>a1>0,a3<0;
步骤三、令α<0,利用公式(2)对名义加速度a1,a2,a3进行分数阶积分,即得到所需的加速段速度规划曲线,此曲线为类S型曲线;
步骤四、被控对象第i个运动周期的速度为v(i),运动周期即时间步长为h;位移增量为d(i)则到第N周期结束,即总位移增量;
步骤五、t1时刻达到最大速度Vmax,第N周期走过的步数为m=t1/h,利用m=N根据公式(3)求得位移为D(m);
步骤六、被控对象按照类S型曲线进行运动,若行程s<D(m)时,被控对象未达到最大速度Vmax则从减速点第i步开始减速直到速度为0,满足s/2=D(i);
步骤七、被控对象按照类S型曲线进行运动;若行程s>D(m)时,被控对象达到了最大速度Vmax后,经过匀速Vmax过程后再从减速点第j步开始减速到速度为0,减速点满足s-D(j)=D(m);
步骤八、若行程s=D(m)时,被控对象达到最大速度Vmax后,被控对象按照类S型曲线进行减速到速度为0;即完成了一种基于分数阶微积分的数控***速度规划方法。
2.根据权利要求1所述一种基于分数阶微积分的数控***速度规划方法,其特征在于:步骤一中根据被控对象运动时间t,利用连续函数y=f(t),推导分数阶导数具体过程为:
(1)对于连续函数y=f(t),根据经典的求导公式,函数y=f(t)的一阶求导公式为:
(2)根据步骤(1)得到的f'(t)求导为二阶导数具体为:
(3)根据步骤(2)得到的f″(t)求导为三阶导数具体为:
(4)利用数学归纳法,根据步骤(1)~(3)的推导过程推导出f(t)的n阶导数具有以下形式:
其中,
由此,推论到分数阶导数如下:
3.根据权利要求1所述一种基于分数阶微积分的数控***速度规划方法,其特征在于:步骤三中令α<0,利用公式(2)对名义加速度a1,a2,a3的数学表达式为:
利用公式(2)对名义加速度a(t)进行α阶分数阶积分,即可得到加速段的速度规划曲线。
4.根据权利要求1所述一种基于分数阶微积分的数控***速度规划方法,其特征在于:步骤三中α设定为-1~-2之间。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510358118.6A CN104950821B (zh) | 2015-06-25 | 2015-06-25 | 一种基于分数阶微积分的数控***速度规划方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510358118.6A CN104950821B (zh) | 2015-06-25 | 2015-06-25 | 一种基于分数阶微积分的数控***速度规划方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104950821A CN104950821A (zh) | 2015-09-30 |
CN104950821B true CN104950821B (zh) | 2017-07-28 |
Family
ID=54165545
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510358118.6A Expired - Fee Related CN104950821B (zh) | 2015-06-25 | 2015-06-25 | 一种基于分数阶微积分的数控***速度规划方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104950821B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2017113069A1 (zh) * | 2015-12-28 | 2017-07-06 | 深圳配天智能技术研究院有限公司 | S型曲线规划方法、装置及数控机床 |
CN106898034B (zh) * | 2017-02-28 | 2020-04-03 | 北京润科通用技术有限公司 | 一种目标运动轨迹模拟方法和装置 |
CN110248082B (zh) * | 2018-11-20 | 2020-12-08 | 浙江大华技术股份有限公司 | 一种规划摄像机运动控制曲线的方法及装置 |
CN109765851B (zh) * | 2019-02-15 | 2021-02-05 | 大族激光科技产业集团股份有限公司 | 轨迹运动的速度计算方法、装置、存储介质及计算机设备 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101853013A (zh) * | 2009-04-01 | 2010-10-06 | 中国科学院沈阳计算技术研究所有限公司 | 一种用于数控机床高速加工的加减速控制方法 |
CN104181860A (zh) * | 2014-08-25 | 2014-12-03 | 浙江理工大学 | 数控机床s型加减速控制方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011123777A (ja) * | 2009-12-14 | 2011-06-23 | Panasonic Corp | 数値制御データの作成方法 |
-
2015
- 2015-06-25 CN CN201510358118.6A patent/CN104950821B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101853013A (zh) * | 2009-04-01 | 2010-10-06 | 中国科学院沈阳计算技术研究所有限公司 | 一种用于数控机床高速加工的加减速控制方法 |
CN104181860A (zh) * | 2014-08-25 | 2014-12-03 | 浙江理工大学 | 数控机床s型加减速控制方法 |
Non-Patent Citations (5)
Title |
---|
Combining differential evolution and particle swarm optimization to tune and realize fractional-order controllers;Guido Maione等;《Mathematical & Computer Modelling of Dynamical Systems》;20121119(第19期);277-299 * |
分数阶混沌***的控制及同步研究;孙光辉;《中国博士学位论文全文数据库 基础学科辑》;20110815(第08期);29-32页 * |
基于PC的PCB钻孔机数控***的研究与开发;郭钊;《中国优秀硕士学位论文全文数据库 工程科技I辑》;20130415(第04期);52-58页 * |
数控***S曲线加减速规划研究;石川等;《中国机械工程》;20070625;第18卷(第12期);1421-1424 * |
数控***的插补算法及加减速控制方法研究;赵巍;《中国博士学位论文全文数据库 工程科技I辑》;20061115(第11期);72-81页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104950821A (zh) | 2015-09-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104950821B (zh) | 一种基于分数阶微积分的数控***速度规划方法 | |
CN102540978B (zh) | 面向高速加工的表面质量优先的样条曲线实时插补方法 | |
CN102419570B (zh) | 数控机床高速加工的加减速前瞻控制方法 | |
CN102799146B (zh) | 用于数控设备控制***的s形加减速控制的速度规划方法 | |
CN104181860B (zh) | 数控机床s型加减速控制方法 | |
CN104678899A (zh) | 曲线速度规划方法、装置及其数控加工路径数据处理方法 | |
CN103941647B (zh) | 用于嵌入式数控设备的柔性加减速控制插补方法 | |
CN106020122B (zh) | 基于牛顿迭代的数控轨迹控制方法 | |
CN110196590A (zh) | 一种用于机器人路径跟踪的时间最优轨迹规划***及方法 | |
CN111341104B (zh) | 一种交通流跟驰模型的速度时滞反馈控制方法 | |
CN101833306A (zh) | 一种基于前瞻-滤波技术的多程序段连续加减速控制方法 | |
CN103076759A (zh) | 一种获取数控***的s型曲线的方法及其装置 | |
CN103279069A (zh) | 用于数控机床的三角函数二阶连续可导加减速算法 | |
CN102566496B (zh) | 适用于数控装置的进给速度动态实时前瞻控制方法 | |
Wei et al. | Smooth transonic steady states of hydrodynamic model for semiconductors | |
CN103246236B (zh) | 一种考虑小直线段和连接点速度的数控***速度控制方法 | |
CN103312248A (zh) | 一种基于dsp的直线加减速拐点误差补偿方法 | |
CN105278462A (zh) | 机器人控制***轨迹加减速插补算法 | |
CN107368041B (zh) | 一种基于移动平均算法的s曲线加减速控制方法 | |
CN111381565B (zh) | 三次多项式s曲线加减速的加加速度递归计算限制方法 | |
CN109506606B (zh) | 一种预测边坡临滑时间的线性速度倒数法 | |
CN103870655A (zh) | 一种适用于伺服转塔冲床的静音加工控制方法 | |
Lee et al. | S-curve profile switching method using fuzzy system for position control of DC motor under uncertain load | |
CN103246237B (zh) | 基于局部结束点的数控***速度控制方法 | |
CN105955029A (zh) | 一种保鲁棒性的pid控制参数优化方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170728 Termination date: 20210625 |
|
CF01 | Termination of patent right due to non-payment of annual fee |