本发明涉及一种基于三角函数的改进s曲线加减速控制方法，属于运动控制领域。
背景技术：
：：加减速控制是数控系统中的关键技术，良好的加减速控制方法可以有效的避免数控设备运动过程中出现的冲击、失步或振动等现象，提高加工精度和表明质量，同时实现数控机床的快速反应，短时间内达到指定速度，缩短加减速时间，提高生产效率。目前数控系统中常用的加减速控制方法有：直线型加减速、指数型加减速、三角函数加减速和s形加减速。直线型加减速算法简单、耗时短，但是加速度曲线不连续，存在速度突变，易出现各种振动和噪声，对数控设运动造成柔性冲击，适用于对运动精度要求不高的低速、低成本的数控系统。指数型加减速算法平滑性比直线型加减速要好，运动精度高，但其加减速起点和终点加速度突变，同样存在柔性冲击。三角函数加减速，主要是利用0-π之间的正弦曲线构建速度曲线，进而实现加速度与加加速度的连续可导，从而达到速度平稳光滑过度的目的，但是由于当前的三角函数加减速为保证速度线型，无法充分发挥数控设备加加速度、加速度最大允许值的性能，造成速度无法在较短时间、距离内达到期望值。也有学者使用三角函数构建加速度和加加速度曲线，通过积分获得速度与位移曲线，同样也存在为保证三角函数线型而无法充分发挥数控设备的性能的问题。s形曲线加减速具有加速度曲线连续、速度曲线光滑等优点，适合高速高精加工场的应用，但目前的s曲线仍存在加加速度的阶跃变化，在阶跃变化的上升沿和下降沿，对数控设备进给系统造成振动和冲击，影响运行平稳性。中国专利文件(申请号201310095677.3)公开了一种用于数控机床的三角函数二阶连续可导加减速算法，结合通用s型加减速曲线算法的优势，提出了改进型的三角函数算法，但该方案在对加加速度处理时，加加速度直接从最大加加速度以三角函数形式降低至0，从开始便有较大冲击，并且加加速度到达最大值时并未维持，加速效率较低，耗时长。从以上分析可知，s曲线加减速与三角函数加减速是采用不同的曲线实现来完成速度的平稳光滑过度。在曲线的连续性上三角函数具有较优的能力，但无法充分发挥数控设备最大加速度、加加速度的性能，执行效率低。而s曲线具有良好的速度光滑性，但加加速度不连续，对数控设备进给系统有振动和冲击。技术实现要素：：针对现有数控系统使用的加减速控制方法存在的问题，本发明的目的是提供一种基于三角函数的改进s曲线加减速控制方法，实现加加速度的连续控制，保证运动过程平稳、无冲击；根据加加速度、加速度、速度和位移之间的积分关系，推导出计算公式，并给出该加减速控制方法的速度规划方法。为实现以上目的，本发明的技术方案如下：定义：t为时间，j(t)、a(t)、v(t)、s(t)分别为加加速度、加速度、速度和位移随时间的变化关系，jcom、acom、vcom分别为数控设备允许最大加加速度、最大加速度和最大速度，vmax为实际能够达到的最大速度，vs为起点速度，ve为终点速度，s为待插补的位移。一种基于三角函数的改进s曲线加减速控制方法，将曲线按照不同的运动特点分阶段，规划过程中通过参数限制和实际情况讨论各阶段是否存在，可以获得各阶段的运行时间，进而可以得到规划后的位移曲线。一种基于三角函数的改进s曲线加减速控制方法，具有与传统s曲线相同的7段式结构，其曲线按时间顺序包括七个阶段，分别是：加加速段l1、匀加速段l2、减加速段l3、匀速段l4、加减速段l5、匀减速段l6和减减速段l7；对应这7个运动阶段，总运动时间t也被分成了7部分，每个阶段对应的时间为：加加速段时间δt1，δt1＝t1；匀加速段时间δt2，δt2＝t2-t1；减加速段时间δt3，δt3＝t3-t2；匀速段时间δt4，δt4＝t4-t3；加减速段时间δt5，δt5＝t5-t4；匀减速段时间δt6，δt6＝t6-t5；减减速段时间δt7，δt7＝t7-t6；其中ti为各阶段过渡点的时刻，i＝1,2...,7，δti为各阶段时间，i＝1,2...,7；(1.1)将l1、l3、l5、l7段的时间分别分割为三部分，tij为分割后过渡点的时刻，i＝1,3,5,7；j＝1,2,3；通过改进传统s曲线在l1、l3、l5、l7段的加加速度变化规律，引入三角函数形式的曲线形式，能够保证加加速度在整个运动过程中的连续性，同时保留匀加加速度段，能够充分发挥数控设备的加减速性能。但引入三角函数之后，l1、l3、l5、l7段的时间将会继续分割，以l1段为例，在进行修改之后，原加加速段时间t1被分割成加加加速段t11、匀加加速段t12-t11和减加加速段t13-t12三部分(t13＝t1，t13所对应的时刻即为分割前t1所指的时刻)。因此致使总时间t被分成15段，大大增加了计算的难度。为简化计算，优选的，步骤(1.1)中，将l1、l3、l5、l7段的时间分别进行三等分分割，设t1为l1、l3段内时间三等分分割之后的分段时间长度，设t2为l5、l7段内时间三等分分割之后的分段时间长度，由于加速度的对称关系，存在如下关系t11＝t12-t11＝t13-t12＝t31-t2＝t32-t31＝t33-t32＝t1，t51＝t52-t51＝t53-t52＝t71-t6＝t72-t71＝t73-t72＝t2(1)。以l1段为例，t11＝t12-t11＝t13-t12＝t1，既能降低计算难度，又能保证存在匀加加速段，充分发挥数控设备的性能。同时有如下关系t13＝t1，t33＝t3，t53＝t5，t73＝t7(2)，t13所对应的时刻即为分割前t1所指的时刻，t33所对应的时刻即为分割前t3所指的时刻，t53所对应的时刻即为分割前t5所指的时刻，t73所对应的时刻即为分割前t7所指的时刻；分割时间段后，其加加速度表达式如下：τij(i＝1,2,..,7；j＝1,2,3)为以各个区间段的起点作为时间零点的时间表示。存在如下关系：τij＝t-tij(i＝1,3,5,7；j＝1,2,3)或τi＝t-ti(i＝2,4,6)(4)上述加加速度与加速度、速度、位移的积分关系如下：根据加加速度与加速度之间的积分关系，可得加速度表达式如下：其中a12＝a11+jcomt1，a32＝a31-jcomt，a52＝a51-jcomt2，a72＝a71+jcomt2。a2为加速段所能达到的最大加速度，后续规划计算中，在不同假设情况下使用aacci(i＝1,2,3,4)表示a2，即aacci＝a2。同理-a6为减速段能达到的最大加速度，后续规划计算中使用adeci(i＝1,2,3,4)表示-a6，即adeci＝(-a2)。根据速度与加速度之间的积分关系，可得速度表达式为：其中v2＝v13+a2(t2-t13)；v6＝v53+a6(t6-t53)；根据速度与位移之间的积分关系，可得位移表达式为：其中s4＝s33+v4(t4-t33)；具体规划包括步骤如下：(2.1)假设能够达到的最大速度为vmax1，且vmax1＝vcom，(2.1.1)加速段时间计算，所述加速段包括加加速段l1、匀加速段l2、减加速段l3；首先分析匀加速段l2的存在性，假设加速段实际能够达到的最大加速度aacc1为数控设备允许的最大加速度acom，即aacc1＝acom(11)根据加加速度与加速度的积分关系，以及加加速段、减加速段的时间等分关系，可获得下式：3t1＝δt1(13)；联立(12)-(13)，可获得加加速段时间且存在关系δt1＝δt3(15)；假设匀加速段l2不存在时，仅通过加加速段和减加速段，根据速度关系，速度增量为δv＝aacc1δt1(16)通过比较速度增量δv与(vmax1-vs)的关系，可知匀加速段l2的存在性，判断如下(2.1.1.a)δv＜(vmax-vs)此时说明存在匀加速段l2，匀加速段的时间为：(2.1.1.b)δv＝(vmax-vs)此时说明恰好不需要匀加速段，无匀加速段l2；(2.1.1.c)δv＞(vmax-vs)此时说明不存在匀加速段l2，且需要重新计算加加速段l1和减加速段l3的时间，由速度关系可得自此，已求出三种不同情况a)、b)、c)的δt1，δt2，δt3；此时加速段实际能够达到的最大加速度为加速段运动位移为(2.1.2)减速段时间计算，所述减速段包括加减速段l5、匀减速段l6和减减速段l7；减速段的处理方式与加速段相同，即减速段可以看成反向的加速段，通过上述算法也可获得减速段各段时间δt5、δt6、δt7，计算过程不再赘述，先分析匀减速段l6的存在性，假设减速段实际能够达到的最大加速度adec1为数控设备允许的最大加速度acom，根据加加速度与加速度的积分关系，以及加减速段、减减速段的时间等分关系，可求得δt5，δt7＝δt5；假设匀减速段l6不存在，仅通过加减速段和减减速段，根据速度关系，速度增量为δv＝aacc1δt1，通过比较速度增量δv与vmax-ve的关系，可知匀减速段l6的存在性，进而求出三种不同情况的δt5、δt6、δt7。减速段位移为(2.1.3)匀速段时间计算根据速度与位移的积分关系，可推导出加速段与减速段的运动总位移为：s1＝sacc1+sdec1(22)通过比较s1与给定的待插补位移s的关系可知匀速段的存在性，如下：(2.1.3.a)s1＜s此时存在匀速段，匀速段时间为至此，所有的7段的时间规划已经完成，实际能够达到的最大速度为vmax1＝vcom，结束计算；(2.1.3.b)s1＝s此时存在匀速段不存在，匀速段时间δt4＝0(24)至此，所有的7段的时间规划已经完成，实际能够达到的最大速度为vmax1＝vcom，结束计算；(2.1.3.c)s1＞s此时说明假设的最大速度vmax1大于实际所能达到的最大速度vmax，需要继续减小vmax1，进入(2.2)；(2.2)假设能够达到的最大速度为vmax2，且vmax2为vs与ve中的较大者进行不包含匀加速段的加速过程所能达到的速度；此时不存在匀速段，即δt4＝0；设定vs≥ve；(下面的算法有几个需要选择的部分，设定一个较大值只是为了表述具体直观)(2.2.1)加速段时间计算对于加速段，根据(2.2)假设条件，不存在匀加速段，即δt2＝0；能够到达的加速度aacc2为数控设备允许的最大加速度acom，aacc2＝acom(25)；根据加速度关系，加加速段与减加速段时间相等，即δt1＝δt3，有求得加加速段和减加速段时间为根据速度关系，可得最大速度为vmax2＝aacc2×δt1+vs(28)比较vmax2与vcom的关系，如下：(2.2.1.a)vmax2＞vcom此时说明假设的最大速度大于实际所能达到的最大速度，若vs＞ve，进入(2.3)；若vs＝ve，进入(2.4)；(2.2.1.b)vmax2≤vcom此时说明假设最大速度可能是满足条件的，继续进行计算；由位移条件，加速段位移为进入(2.2.2)；(2.2.2)减速段时间计算ve≤vs，因此减速段实际能够达到的最大加速度adec2为数控设备允许的最大加速度acom，adec2＝acom(30)减减速段与加减速段时间为如果ve＜vs，则存在匀减速段，有adec(δt5+δt6)＝vmax2-ve(32)则匀减速段时间为如果ve＝vs，则不存在匀减速段，匀减速段时间δt6＝0；由位移关系，减速段位移为判断sacc2+sdec2与给定的待插补距离s的关系，如下：(2.2.2.a)sacc2+sdec2＝s此时说明假设的最大速度vmax2恰好为实际能够达到的最大速度vmax，所有7段时间都已经计算完成，结束计算；(2.2.2.b)sacc2+sdec2＞s此时说明假设的最大速度vmax2仍然大于实际能够达到的最大速度vmax，进入(2.3)；(2.2.2.c)sacc2+sdec2＜s此时说明假设的最大速度vmax2小于实际能达到的最大速度vmax，其取值区间为vmax∈(vmax2,vcom)，同时存在匀加速段和匀减速段，但不包含匀速段；具体计算如下：加速段与减速段实际能够达到的最大加速度为aacc2＝adec2＝acom(35)加加速段、减加速段、减减速段、加减速段的时间为由速度关系，有vmax-vs＝acom(δt1+δt2)(37)可得匀加速段时间为同理可得匀减速段时间vmax-ve＝acom(δt5+δt6)(39)；由位移关系，有联立(35)-(41)，由总插补位移关系构造方程其中b＝t1，f(vmax)是vmax二次函数，由于b＞0，所以根据求根公式得通过将vmax带入前述各段时间的计算式，即可计算出各段时间，(各段时间指的是δt1，δt2，δt3、δt4、δt5、δt6、δt7)，完成计算；(2.3)假设能够达到的最大速度为vmax3，且vmax3为vs与ve中的较小者进行不包含匀加速段(此处是设定一个最大速度，两端速度都看成向最大速度加速的过程)的加速过程所能达到的速度；此时不存在匀速段，δt4＝0；设定vs＜ve，(2.3.1)减速段时间计算对于减速段，根据(2.3)假设条件，不存在匀减速段，即δt6＝0；能够到达的加速度adec3为数控设备允许的最大加速度acom，即adec3＝acom(44)根据加速度关系，减减速段与加减速段时间相等，即δt5＝δt7有可求得减减速段与加减速段时间为根据速度关系，可得最大速度为vmax3＝adec2×δt5+ve(47)比较vmax3与vcom的关系，如下：(2.3.1.a)vmax3＞vcom此时说明假设的最大速度大于实际所能达到的最大速度，进入(2.4)；(2.3.1.b)vmax3≤vcom此时说明假设最大速度可能是满足条件的，继续进行计算，由位移条件，减速段位移为(2.3.2)加速段时间计算由于ve＜vs，因此没有匀加速段，即δt2＝0，设加速段所能达到的最大加速度设为aacc3，由加速度关系，加加速段与减加速段时间为由位移关系，加速速段位移为判断sacc3+sdec3与给定的待插补距离s的关系，如下：(2.3.2.a)sacc3+sdec3＝s此时说明假设的最大速度vmax3恰好为实际能够达到的最大速度vmax，所有7段时间都已经计算完成，结束计算；(2.3.2.b)sacc3+sdec3＞s此时说明假设的最大速度vmax3仍然大于实际能够达到的最大速度vmax，进入(2.4)；(2.3.2.c)sacc3+sdec3＜s此时说明假设的最大速度vmax3小于实际能达到的最大速度vmax，其取值区间为vmax∈(vmax3,vmax2)，包含匀减速段，但不包含匀加速段和匀速段，具体计算如下：由速度关系，加加速段与减加速段时间为加速段所能达到的最大加速度为加速段位移为对于减速段，所能达到的最大加速度为acom，减减速段和减加速段时间为由速度关系，有vmax-ve＝acom(δt5+δt6)(55)可得匀减速段时间为减速段位移为由位移关系，可得s3(vmax)＝sacc3+sdec3(58)联立(51)-(58)，由总位移关系，构造方程因为s3(vmax)是vmax的单调递增函数，所以方程f3(vmax)在(vmax3,vmax2)内有且有唯一解，可使用二分法或牛顿迭代法进行求解，根据求解出的vmax可计算出δt6，至此，已求解出各段时间，完成计算；(2.4)最大速度vmax的取值范围为[max(vs,ve),vmax3)此时匀速段、匀加速段和匀减速段都不存在，设定vs≥ve，进行推导计算：首先假设最大速度为起点速度，即vmax4＝vs，此时加加速段与减加速段都不存在，加速段位移sacc4为0，仅存在减减速段和加减速段，由加速度关系，减速段实际能达到的最大加速度为由速度关系，减减速段与加减速段时间相等，即δt5＝δt7，有vmax4-ve＝adec4δt5(61)联立(60)-(61)可得减减速段与加减速段时间为由位移关系，减速段位移为比较sdec4与待插补位移s的关系，如下：(2.4.a)sdec4＞s此时说明假设的最大速度vmax4要小于起点速度vs，因此在给定条件下不存在最大速度，无解，结束计算；(2.4.b)sdec4＝s此时说明假设的最大速度vmax4恰好等于实际能够达到的最大速度vmax，即vmax＝vmax4，结束计算；(2.4.c)sdec4＜s此时说明实际能够达到的最大速度vmax应大于vmax4，其取值范围为vmax∈(vmax4,vmax3)，由加速度关系，加加速段与减加速段时间相等，即δt1＝δt3，加速段实际能够达到的加速度为由速度关系，有vmax-vs＝aacc4δt1(65)由位移关系，可得加速段位移为同理，有减速段相关关系vmax-ve＝adec4δt5(68)；由位移关系得s4(vmax)＝sacc4+sdec4(70)由总位移关系，联立(63)-(70)，构建方程s(vmax)为vmax的单调递增函数，因此在(vmax4,vmax3)内有且有唯一解，可使用二分法或牛顿迭代法求解，获得vmax后，根据求解出的vmax可计算出δt6，至此，已求解出各段时间，完成计算。基于三角函数的改进s曲线规划原理实际使用基于三角函数的改进s曲线加减速方法进行任务规划时，可能并不完全包含上述给出的全部7段，会根据给定的最大加加速度jcom、最大加速度acom、最大速度vcom、起点速度vs、终点速度ve、待规划位移s的不同，存在以下表1所示几种情况：表1情况列表包含不包含l1、l2、l3、l4、l5、l6、l7l1、l2、l3、l5、l6、l7l4l1、l3、l4、l5、l6、l7l2l1、l2、l3、l4、l5、l7l6l1、l3、l4、l5、l7l2、l6l1、l3、l5、l6、l7l2、l4l1、l2、l3、l5、l7l6、l4l1、l3、l5、l7l2、l4、l6l5、l7l1、l2、l3、l4、l6l1、l3l2、l4、l5、l6、l7但利用本申请的规划方法可适应上述所有的情况。本发明的优点在于：1)本发明在传统s曲线加减速的基础上，引入基于三角函数的加加速度控制方法，能够有效的保证加加速度的连续性，最大限度的减小对数控设备造成的冲击。2)基于三角函数的加加速度的控制中，保留了匀加加速度段和匀减减速度段，能够充分发挥数控设备最大允许加加速度和加速度的性能，保证较高的生产效率。3)基于三角函数的加加速度的控制中，提出对加加加速度段、匀加加速度段和减加加速度段三部分运动时间进行等分控制的方法，能够有效简化计算，保证算法的实时性能。4)本申请是加加速度是从0以三角函数形式增大至最大加加速度，因此无冲击，且本方案可将加加速度维持在最大值，直到合适时再下降，速度变化效率高。附图说明图1为根据本发明基于三角函数的改进s曲线加减速控制方法改进后的标准曲线示意图图2为不包含匀速段的基于三角函数的改进s曲线加减速加速度曲线图3为不包含匀加速段的基于三角函数的改进s曲线加减速加速度曲线图4为不包含匀减速段的基于三角函数的改进s曲线加减速加速度曲线图5为不包含匀加速段和匀减速段的基于三角函数的改进s曲线加减速加速度曲线图6为不包含匀加速段和匀速段的基于三角函数的改进s曲线加减速加速度曲线图7为不包含匀减速段和匀速段的基于三角函数的改进s曲线加减速加速度曲线图8为不包含匀加速段、匀速段和匀减速段的基于三角函数的改进s曲线加减速加速度曲线图9为仅包含减减速段和加减速段的基于三角函数的改进s曲线加减速加速度曲线图10为仅包含加加速段和减加速段的基于三角函数的改进s曲线加减速加速度曲线具体实施方式下面通过实施例并结合附图对本发明做进一步的说明，但不仅于此。实施例1：如图1所示，基于三角函数的改进s曲线加减速控制方法具有7段式结构，分别是：加加速段l1、匀加速度段l2、减加速段l3、匀速段l4、加减速段l5、匀减速段l6和减减速度段l7。对应这7个运动阶段，总运动时间t也被分成了7部分，分别是加加速段时间t1(设为δt1)、匀加速段时间t2-t1(设为δt2)、减加速段时间t3-t2(设为δt3)、匀速段时间t4-t3(设为δt4)、加减速段时间t5-t4(设为δt5)、匀减速段时间t6-t5(设为δt6)和减减速段时间t7-t6(设为δt7)，其中ti(i＝1,2...,7)为各阶段过渡点的时刻，δti(i＝1,2...,7)为各阶段时间。为简化计算，对l1、l3、l5、l7段时间进行三等分，由于加速度的对称关系，存在如下关系t11＝t12-t11＝t13-t12＝t31-t2＝t32-t31＝t33-t32＝t1，t51＝t52-t51＝t53-t52＝t71-t6＝t72-t71＝t73-t72＝t2(1)设t1为l1、l3段内时间三等分分割之后的分段时间长度，设t2为l5、l7段内时间三等分分割之后的分段时间长度，其中tij(i＝1,3,5,7；j＝1,2,3)为通过三等分方法新获得的过度时刻，同时有如下关系t13＝t1，t33＝t3，t53＝t5，t73＝t7(2)分割时间段后，其加加速度表达式如前述公式(3)：τij(i＝1,2,..,7；j＝1,2,3)为以各个区间段的起点作为时间零点的时间表示。存在如下关系：τij＝t-tij(i＝1,3,5,7；j＝1,2,3)或τi＝t-ti(i＝2,4,6)(4)上述加加速度与加速度、速度、位移的积分关系如前述公式(5)、(6)、(7)。根据加加速度与加速度之间的积分关系，可得加速度表达式如前述公式(8)。a2为加速段所能达到的最大加速度，后续规划计算中，在不同假设情况下使用aacci(i＝1,2,3,4)表示a2，即aacci＝a2。同理-a6为减速段能达到的最大加速度，后续规划计算中使用adeci(i＝1,2,3,4)表示-a6，即adeci＝(-a2)。根据速度与加速度之间的积分关系，可得速度表达式如前述公式(9)。根据速度与位移之间的积分关系，可得位移表达式如前述公式(10)。如图1-图10所示，实际使用基于三角函数的改进s曲线加减速方法进行任务规划时，可能并不完全所有的7段，会根据给定的最大加加速度jcom、最大加速度acom、最大速度vcom、起点速度vs、终点速度ve、待规划位移s的不同而存在多种情况。规划过程中，通过求解能够真实到达的最大加速度vmax，就可以获得各段的运行时间及每个时刻的速度特性。具体规划如下：(2.1)假设能够达到的最大速度为vmax1，且vmax1＝vcom。(2.1.1)加速段时间计算，所述加速段包括加加速段l1、匀加速段l2、减加速段l3；首先分析匀加速段的存在性，假设加速段实际能够达到的最大加速度aacc1为设备允许的最大加速度acom，即aacc1＝acom(11)根据加加速度与加速度的积分关系，以及加加速段、减加速段的时间等分关系，可获得下式：3t1＝δt1(13)联立(12)-(13)，可获得加加速段时间且存在关系δt1＝δt3(15)；假设匀加速段不存在时，仅通过加加速段和减加速段，根据速度关系，速度增量为δv＝aacc1δt1(16)通过比较速度增量δv与(vmax1-vs)的关系，可知匀加速段的存在性，如下(2.1.1.a)δv＜(vmax-vs)此时说明存在匀加速段，匀加速时间为：(2.1.1.b)v＝(vmax-vs)此时说明恰好不需要匀加速段。(2.1.1.c)v＞(vmax-vs)此时说明不存在匀速段，且需要重新计算加加速段和减加速段时间，由速度关系可得此时加速段实际能够达到的最大加速度为加速段运动位移为(2.1.2)减速段时间计算，所述减速段包括加减速段l5、匀减速段l6和减减速段l7；减速段的处理方式与加速段相同，可获得减速段各段时间δt5、δt6、δt7，计算过程不再赘述。减速段位移为(2.1.3)匀速段时间计算根据速度与位移的积分关系，可推导出加速段与减速段的运动总位移为：s1＝sacc1+sdec1(22)通过比较s1与给定的待插补位移s的关系可知匀速段的存在性，如下：(2.1.3.a)s1＜s此时存在匀速段，匀速段时间为至此，所有的7段的时间规划已经完成，实际能够达到的最大速度为vmax1＝vcom，结束计算。此时规划得到的各段加速度会出现如图1、图3、图4或图5所示的形式。(2.1.3.b)s1＝s此时存在匀速段不存在，匀速段时间δt4＝0(24)至此，所有的7段的时间规划已经完成，实际能够达到的最大速度为vmax1＝vcom，结束计算。此时规划得到的各段加速度会出现如图2、图6、图7或图8所示的形式。(2.1.3.c)s1＞s此时说明假设的最大速度vmax1大于实际所能达到的最大速度vmax，需要继续减小vmax1，进入(2.2)。(2.2)假设能够达到的最大速度为vmax2，且vmax2为vs与ve中的较大者进行不包含匀加速段的加速过程所能达到的速度。此时不存在匀速段，即δt4＝0。下面以vs≥ve为例进行推导计算。(2.2.1)加速段时间计算对于加速段，根据(2.2)假设条件，不存在匀加速段，即δt2＝0。能够到达的加速度aacc2为数控设备允许的最大加速度acom，即aacc2＝acom(25)根据加速度关系，加加速段与减加速段时间相等，即δt1＝δt3。有可求得加加速段和减加速段时间为由此速度关节，可得最大速度为vmax2＝aacc2×δt1+vs(28)比较vmax2与vcom的关系，如下：(2.2.1.a)vmax2＞vcom此时说明假设的最大速度大于实际所能达到的最大速度，若vs＞ve，进入(2.3)；若vs＝ve，进入(2.4)。(2.2.1.b)vmax2≤vcom此时说明假设最大速度可能是满足条件的，继续进行计算；由位移条件，加速段位移为进入(2.2.2)。(2.2.2)减速段时间计算由于ve≤vs，因此减速段实际能够达到的最大加速度adec2为设备允许的最大加速度acom，即adec2＝acom(30)减减速段与加减速段时间为如果ve＜vs，则存在匀减速段，有adec(δt5+δt6)＝vmax2-ve(32)则匀减速段时间为如果ve＝vs，则不存在匀减速段，匀减速段时间δt6＝0。由位移关系，减速段位移为判断sacc2+sdec2与给定的待插补距离s的关系，如下：(2.2.2.a)sacc2+sdec2＝s此时说明假设的最大速度vmax2恰好为实际能够达到的最大速度vmax，所有7段时间都已经计算完成，结束计算。此时规划得到的各段加速度会出现如图6或图8所示的形式。如果vs≤ve，各段加速度可能会出现如图7或图8所示的形式。(2.2.2.b)sacc2+sdec2＞s此时说明假设的最大速度vmax2仍然大于实际能够达到的最大速度vmax，进入(2.3)。(2.2.2.c)sacc2+sdec2＜s此时说明假设的最大速度vmax2小于于实际能达到的最大速度vmax，其取值区间为vmax∈(vmax2,vcom)，同时存在匀加速段和匀减速段，但不包含匀速段。具体计算如下：加速段与减速段实际能够达到的最大加速度为aacc2＝adec2＝acom(35)加加速段、减加速段、减减速段、加减速段的时间为由速度关系，有vmax-vs＝acom(δt1+δt2)(37)可得匀加速段时间为同理可得匀减速段时间vmax-ve＝acom(δt5+δt6)(39)；由位移关系，有联立(35)-(41)，由总插补位移关系构造方程其中b＝t1，f(vmax)是vmax二次函数，由于b＞0，所以根据求根公式得通过将vmax带入前述各段时间的计算式，即可计算出各段时间，(各段时间指的是δt1，δt2，δt3、δt4、δt5、δt6、δt7)，完成计算。此时规划得到的各段加速度可能出现如图2所示的形式。(2.3)假设能够达到的最大速度为vmax3，且vmax3为vs与ve中的较小者进行不包含匀加速段的加速过程所能达到的速度。此时不存在匀速段，即δt4＝0。下面以vs＜ve为例进行推导计算。(2.3.1)减速段时间计算对于减速段，根据(2.3)假设条件，不存在匀减速段，即δt6＝0。能够到达的加速度adec3为数控设备允许的最大加速度acom，即adec3＝acom(44)根据加速度关系，减减速段与加减速段时间相等，即δt5＝δt7有可求得减减速段与加减速段时间为由此速度关节，可得最大速度为vmax3＝adec2×δt5+ve(47)比较vmax3与vcom的关系，如下：(2.3.1.a)vmax3＞vcom此时说明假设的最大速度大于实际所能达到的最大速度，进入(2.4)。(2.3.1.b)vmax3≤vcom此时说明假设最大速度可能是满足条件的，继续进行计算。由位移条件，减速段位移为(2.3.2)加速段时间计算由于ve＜vs，因此没有匀加速段，即δt2＝0。设加速段所能达到的最大加速度为aacc3。由加速度关系，加加速段与减加速段时间为由位移关系，加速速段位移为判断sacc3+sdec3与给定的待插补距离s的关系，如下：(2.3.2.a)sacc3+sdec3＝s此时说明假设的最大速度vmax3恰好为实际能够达到的最大速度vmax，所有7段时间都已经计算完成，结束计算。此时规划得到的各段加速度可能出现如图8所示的形式。(2.3.2.b)sacc3+sdec3＞s此时说明假设的最大速度vmax3仍然大于实际能够达到的最大速度vmax，进入(2.4)。(2.3.2.c)sacc3+sdec3＜s此时说明假设的最大速度vmax3小于实际能达到的最大速度vmax，其取值区间为vmax∈(vmax3,vmax2)，包含匀减速段，但不包含匀加速段和匀速段。具体计算如下：由速度关系，加加速段与减加速段时间为加速段所能达到的最大加速度为加速段位移为对于减速段，所能达到的最大加速度为acom，减减速段和减加速段时间为由速度关系，有vmax-ve＝acom(δt5+δt6)(55)可得匀减速段时间为减速段位移为由位移关系，可得s3(vmax)＝sacc3+sdec3(58)联立(51)-(58)，由总位移关系，构造方程因为s3(vmax)是vmax的单调递增函数，所以方程f3(vmax)在(vmax3,vmax2)内有且有唯一解，可使用二分法或牛顿迭代法进行求解。根据求解出的vmax可计算出δt6，至此，已求解出各段时间，完成计算。此时规划得到的各段加速度可能出现如图6所示的形式。如果vs＜ve，各段加速度可能出现如图7所示的形式。(2.4)最大速度vmax的取值范围为[max(vs,ve),vmax3)此时匀速段、匀加速段和匀减速段都不存在，下面以vs≥ve为例进行推导计算。首先假设最大速度为起点速度，即vmax4＝vs，此时加加速段与减加速段都不存在，加速段位移sacc4为0，仅存在减减速段和加减速段。由加速度关系，减速段实际能达到的最大加速度为由速度关系，减减速段与加减速段时间相等，即δt5＝δt7，有vmax4-ve＝adec4δt5(61)联立(60)-(61)可得减减速段与加减速段时间为由位移关系，减速段位移为比较sdec4与待插补位移s的关系，如下：(2.4.a)sdec4＞s此时说明假设的最大速度vmax4要小于起点速度vs，因此在给定条件下不存在最大速度无解，结束计算。(2.4.b)sdec4＝s此时说明假设的最大速度vmax4恰好等于实际能够达到的最大速度vmax，即vmax＝vmax4，结束计算。此时规划得到的各段加速度可能出现如图9所示的形式。如果vs≤ve，各段加速度可能出现如图10所示的形式。(2.4.c)sdec4＜s此时说明实际能够达到的最大速度vmax应大于vmax4，其取值范围为vmax∈(vmax4,vmax3)。由加速度关系，加加速段与减加速段时间相等，即δt1＝δt3。加速段实际能够达到的加速度为由速度关系，有vmax-vs＝aacc4δt1(65)由位移关系，可得加速段位移为同理，有减速段相关关系vmax-ve＝adec4δt5(68)；由位移关系得s4(vmax)＝sacc4+sdec4(70)由总位移关系，联立(63)-(70)，构建方程s(vmax)为vmax的单调递增函数，因此在(vmax4，vmax3)内有且有唯一解，可使用二分法或牛顿迭代法求解。获得vmax后，即可计算出各段时间，完成计算。此时规划得到的各段加速度可能出现如图8所示的形式。实验例1利用实施例1所述的控制方法，用c++编程，放入matlab优化算法求得对应图像。本实验例为不包含匀速段的情况，设定表2数据如下。表2初始参数jmaxamaxvcomvsves参数值5×104mm/s31500mm/s2100mm/s10mm/s20mm/s10mm时间(ms)t1＝39.59t2＝18.62t3＝39.59t4＝0t5＝39.59t6＝11.95t7＝39.59位移(mm)0.741.725.025.028.459.1410在matlab中运行后，得出对应图2的曲线图，由表中可知，运行完成后，位移完成设定，由图2可知，整个过程速度平滑。实验例2利用实施例1所述的控制方法，用c++编程，放入matlab优化算法求得对应图像。本实验例为不包含匀加速段的情况，设定表3数据如下。表3初始参数jmaxamaxvcomvsves参数值5×104mm/s31500mm/s2100mm/s60mm/s20mm/s15mm时间(ms)t1＝32.49t2＝0t3＝32.49t4＝42.26t5＝39.59t6＝13.74t7＝39.59位移(mm)2.112.115.039.2212.7613.5715在matlab中运行后，得出对应图3的曲线图，由表中可知，运行完成后，位移完成设定，由图3可知，整个过程速度平滑。实验例3利用实施例1所述的控制方法，用c++编程，放入matlab优化算法求得对应图像。本实验例为不包含匀减速段的情况，设定表4数据如下。表4初始参数jmaxamaxvcomvsves参数值5×104mm/s31500mm/s2100mm/s10mm/s50mm/s20mm时间(ms)t1＝39.59t2＝20.41t3＝39.59t4＝90.73t5＝36.33t6＝0t7＝36.33位移(mm)0.741.845.2414.217.4917.4920在matlab中运行后，得出对应图4的曲线图，由表中可知，运行完成后，位移完成设定，由图4可知，整个过程速度平滑。实验例4利用实施例1所述的控制方法，用c++编程，放入matlab优化算法求得对应图像。本实验例为不包含匀加速段和匀减速段的情况，设定表5数据如下。表5初始参数jmaxamaxvcomvsves参数值5×104mm/s31500mm/s2100mm/s50mm/s55mm/s15mm时间(ms)t1＝36.33t2＝0t3＝36.33t4＝42.09t5＝34.46t6＝0t7＝34.46位移(mm)2.072.075.279.4312.5912.5915在matlab中运行后，得出对应图5的曲线图，由表中可知，运行完成后，位移完成设定，由图5可知，整个过程速度平滑。实验例5利用实施例1所述的控制方法，用c++编程，放入matlab优化算法求得对应图像。本实验例为不包含匀加速段和匀速段的情况，设定表6数据如下。表6初始参数jmaxamaxvcomvsves参数值5×104mm/s31500mm/s2100mm/s70mm/s25mm/s10mm时间(ms)t1＝30.13t2＝0t3＝30.13t4＝0t5＝39.59t6＝16.67t7＝39.59位移(mm)2.262.265.175.178.919.9510在matlab中运行后，得出对应图6的曲线图，由表中可知，运行完成后，位移完成设定，由图6可知，整个过程速度平滑。实验例6利用实施例1所述的控制方法，用c++编程，放入matlab优化算法求得对应图像。本实验例为不包含匀减速段和匀速段的情况，设定表7数据如下。表7初始参数jmaxamaxvcomvsves参数值5×104mm/s31500mm/s2100mm/s20mm/s60mm/s10mm时间(ms)t1＝39.59t2＝13.33t3＝39.59t4＝0t5＝32.24t6＝0t7＝32.24位移(mm)1.121.895.225.228.128.1210在matlab中运行后，得出对应图7的曲线图，由表中可知，运行完成后，位移完成设定，由图7可知，整个过程速度平滑。实验例7利用实施例1所述的控制方法，用c++编程，放入matlab优化算法求得对应图像。本实验例为不包含匀加速段、匀速段和匀减速段的情况，设定表8数据如下。表8初始参数jmaxamaxvcomvsves参数值5×104mm/s31500mm/s2100mm/s20mm/s30mm/s7mm时间(ms)t1＝37.53t2＝0t3＝37.53t4＝0t5＝33.83t6＝0t7＝33.83位移(mm)1.041.043.333.335.555.557在matlab中运行后，得出对应图8的曲线图，由表中可知，运行完成后，位移完成设定，由图8可知，整个过程速度平滑。实验例8利用实施例1所述的控制方法，用c++编程，放入matlab优化算法求得对应图像。本实验例为仅包含减减速段和加减速段的情况，设定表9数据如下。表9初始参数jmaxamaxvcomvsves参数值5×104mm/s31500mm/s250mm/s60mm/s20mm/s8mm时间(ms)t1＝0t2＝0t3＝0t4＝90.01t5＝32.49t6＝0t7＝32.49位移(mm)0005.47.147.148在matlab中运行后，得出对应图9的曲线图，由表中可知，运行完成后，位移完成设定，由图9可知，整个过程速度平滑。实验例9利用实施例1所述的控制方法，用c++编程，放入matlab优化算法求得对应图像。本实验例为仅包含加加速段和减加速段的情况，设定表10数据如下。表10初始参数jmaxamaxvcomvsves参数值5×104mm/s31500mm/s260mm/s30mm/s60mm/s7mm时间(ms)t1＝28.13t2＝0t3＝28.13t4＝74.46t5＝0t6＝0t7＝0位移(mm)0.960.962.457777在matlab中运行后，得出对应图10的曲线图，由表中可知，运行完成后，位移完成设定，由图10可知，整个过程速度平滑。当前第1页12