一种离散加工路径的拐点平滑方法【专利摘要】一种离散加工路径的拐点平滑方法，属于高速高性能数控加工领域。该方法采用外切式的且控制点数为9的5次B样条曲线对离散加工路径的拐点进行平滑处理，解决了常规离散路径拐点加工时速度低、易波动的问题。首先，借助重顶点方法和凸包性质构造出由9个控制点组成的B样条曲线的特征多边形及其节点矢量的具体形式，并在逼近误差允值约束下反算出控制点；进而，基于过渡段长度约束，采用比例调节算法对控制点进行一次性修正；最后，根据求得的特征多边形和节点矢量得到最终用于拐点平滑的五次B样条曲线。本方法可使平滑后的路径达到G2连续，较之内切式方法，在相同逼近误差条件下平滑曲线在拐点处具有更大的曲率半径。【专利说明】-种离散加工路径的拐点平滑方法
技术领域：
[0001]本发明设及一种离散加工路径的拐点平滑方法，属于高速高性能数控加工领域。【
背景技术：
】[0002]在数控加工中，离散加工路径是由众多GO连续直线段构成。在拐点处的路径用GO直线段描述会导致切矢和曲率的不连续，使数控系统在加工过程中频繁的启停，引起加工速度波动，加工表面质量下降W及加工效率低等问题。因此，传统的离散加工路径并不适用于高速高性能数控加工领域，需对刀具路径在拐点处进行平滑处理，其要求有=:1)保证加工直线段与拐点过渡曲线在公共连接点处G2连续;2)保证平滑之后的路径在逼近误差允值内；3)在平滑连续拐点时，过渡段不会出现干设现象。对现有的技术文献检索，找到了几种解决方法。目前，比较常用的一种方法是通过在拐点处构造圆弧，用圆弧作为平滑路径(JouanehMK,WangZ,DornfeldDA(1990)Trajectoryplanningforcoordinatedmotionofarobotandapositioningtable.PartI.Pathspecification.IEEETransRobotAutom6:735-745)。该方法只能保证加工直线段与拐点过渡曲线在公共连接点处Gl连续，即速度连续。Chester等人发明的专利是用四阶多项式曲线来平滑路径(YutkowitzSJ,ChesterW(2005)ApparatusandMethodforSmoothCorneringinaMotionControlSystem.UnitedStates,SiemensEnergyMutomation,Inc,Alpharetta,GA,化SPatent6922606))。运种方法实现了加工直线段与拐点过渡曲线在公共连接点处G2连续，但该方法运算量较大，增加了数据处理的负担。上述方法构造的平滑路径都是内切于原拐点路径，并保证了加工直线段与拐点过渡曲线在公共连接点处一阶或二阶几何连续。当在拐点处两相邻加工直线段的转角比较尖锐时，即转角较小时，在逼近误差允值下拐点平滑曲线的曲率仍然较大，平滑效果相对并不明显，还需要牺牲速度来保证加工质量，不利于高速加工。到目前为止，相关文献还未能很好地解决尖锐拐点路径的平滑问题。【
发明内容】[0003]针对现有技术中的缺陷和不足，本发明提供了一种离散加工路径的拐点平滑方法。[0004]本发明采用的技术方案是:一种离散加工路径的拐点平滑方法，首先，借助重顶点方法和凸包性质构造出由9个控制顶点组成的B样条曲线的特征多边形及其节点矢量的具体形式，并在逼近误差允值约束下反算出控制顶点；进而，基于过渡段长度约束，采用比例调节算法对控制顶点进行一次性修正;最后，根据求得的特征多边形和节点矢量得到最终用于拐点平滑的五次B样条曲线;采用的具体步骤为：[0005](a)构造出由9个控制点组成的外切式B样条曲线的特征多边形及其节点矢量，并通过逼近误差允值约束来反算出控制点；其中，已知拐点化处两相邻线性加工直线段顶点分别为(Qo,化)和(Qi，化），拐点处的最大逼近误差允值为e;拟构造的拐点过渡段的外切式五次B样条曲线表示为：[0006](I)[0007]式I中，P功特征多边形的控制点，Bi日为基函数，U为五次B样条曲线的节点矢量，确定祁0，Pl，...，Ps}为待求的特征多边形的9个控制点，且P4=化、?1=?2、？6=?7而、？1和?2共线并位于直线段9()9止爪、？7和?8共线并位于直线段化92上，从保证加工直线段与拐点过渡曲线在公共连接点处G2连续;控制点P3、PS位于离散点Q〇、Ql、化所围立角形的外部，且Po和P8、Pl和P6、P3和Ps分别关于角ZQo化化的平分线对称，点P3到线段Qo化和点Ps到线段9狐的距离均为最大逼近误差e;[000引由构造的特征多边形对称性和B样条凸包的性质知，最大逼近误差的计算公式为：[0009]E=I|P4-C(0.5)I(2)[0010]进而，由式1得C(0.5)的表达式：[0011]C(0.5)=0.0185Pi+0.化93P3+0.4444P4+0.2593P5+0.0185P6(3)[0012]整理式3后推出：[0013]P广C(0.5)=0.0185P4P1+0.2593P4P3+0.0185P4P6+0.化93P4P5(4)[0014]将式4带入到式2中，则最大逼近误差进一步表达为：[001引倒[0016]式5中，ZQo化化=目，（目声180°):P3P4=L由于点P3到线段PiP4的距离为e，且P3在PiP4的垂直平分线上，则推出下列关系：[0017](6)(7)[001引1234567联立式5、6和7,求解出变量的单位向量为化咕此，计算出控制点P3和Pi,即：2P3=P4+LNi(8)3饼4f式导出：5(10)6由式8、9和10求得修正后的控制点？0、？1爪，同理求得控制点？5瓜爪；7(b)考虑过渡段长度约束，采用比例调节算法对控制顶点进行一次性修正；当数控加工连续拐点路径时，为保证路径的连续可调并且不发生干设现象，即需要对原特征多边形进行约束，比例调节算法如下：[0026](11)[0027](I巧[002引k=max(ki,k2)(13)[0029]由式11、12、13求出比例系数k;若k>I，则需对控制多边形进行等比例缩小，即，古此由式14、15、16求得修正后的控制点Po'、Pi'、P3'：[0030]P3，=P4+b^(14)[0031](I巧[0032](1旬[0033]同理，求得控制点?5'瓜'爪'；修正之后，拐点处逼近误差61也相应变为：[0034](17)[0035]式17表明修正之后的逼近误差Ei与最大逼近误差e相比缩小了k倍化〉1);[0036](C)求出用于拐点平滑的五次B样条曲线;若P3P4=L，则直接将步骤(a)中求得的9个控制点带入式1中，得到五次B样条曲结i若P3'P4=Li，则需将步骤(b)中求得的9个控制点带入式1中，得到五次B样条曲续即得到拐点处的平滑路径。[0037]本发明的有益效果是:采用外切式的且控制点数为9的5次B样条曲线对离散加工路径的拐点进行平滑处理，在逼近误差允值约束下生成曲率连续的加工路径，解决了常规离散路径拐点加工时速度低、易波动的问题。与现有技术相比，在满足同等精度的条件下，本方法在处理尖锐拐点路径时所构造的过渡曲线的曲率半径更大、更平滑，能够有效地改善数控机床在高速加工拐点路径时的运动学性能。【附图说明】[0038]图1是一种应用离散加工路径的拐点平滑方法流程图。[0039]图2是离散加工路径的拐点平滑方法模型示意图。[0040]图3是离散加工路径的拐点平滑方法应用在连续拐点示意图。[0041]图4是内切式的拐点平滑路径示意图。[0042]图5是采用本发明方法的拐点平滑路径示意图。[0043]图6是内切式的拐点平滑路径曲率示意图。[0044]图7是采用本发明方法的拐点平滑路径曲率示意图。【具体实施方式】[0045]-种离散加工路径的拐点平滑方法的流程图如图I所示。下面结合具体的实施例对本发明进行详细说明。W下实施例将有助于本领域的技术人员进一步理解本发明，但不W任何形式限制本发明。应当指出的是，对本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可W做出若干变形和改进。运些都属于本发明的保护范围。[0046](a)构造出由9个控制点组成的外切式B样条曲线的特征多边形及其节点矢量，并通过逼近误差允值约束来反算出控制点。如图2所示(虚线表示加工直线段，实线表示拐点处的平滑路径）：已知拐点化处两相邻线性加工直线段顶点分别为(Q〇，Qi)和(Qi，Q2)，拐点处的最大逼近误差允值为e;拟构造的拐点过渡段的外切式五次B样条曲线可表示为：[0047][004引上述B样条曲线方程中Pi为特征多边形的控制点，BiS为基函数，U为五次B样条曲线的节点矢量，可确定^{P0，Pi，...，Ps}为待求的特征多边形的9个控制点，且P4=化、口1=口2、口6=口7瓜、口1和口2共线并位于直线段舶化上术、口7和口8共线并位于直线段化化上，W保证加工直线段与拐点过渡曲线在公共连接点处G2连续;控制点P3、P5位于离散点舶、化、Q2所围立角形的外部，且Po和P8、Pl和P6、P3和Ps分别关于角ZQoQl化的平分线对称，线段PiP冲点为To，线段PsP冲点为Tl，点P3到线段QoQi和点Ps到线段QiQ2的距离均为巧大Js近"i;??差e,即PsTo=E,P日Ti=E。[0049]由构造的特征多边形对称性和B样条凸包的性质可知，最大逼近误差的计算公式为：[0050]E=I|P4-C(0.5)I[0051]进而，可得C(0.5)的表达式：[0052]C(0.5)=0.0185P1+0.化93P3+0.4444P4+0.2593P5+0.0185P6[0化3]整理上式后可推出：[0054]P广C(0.5)=0.0185P4P1+0.2593P4P3+0.0185P4P6+0.化93P4P已[0055]最大逼近误差可进一步表达为：[0化6][0化7]上式中，ZQo化化=目，（目声180°):P3P4=L。由于PsTo=E,且P3在PlP4的垂直平分线上，则可推出下列关系：出控制点P3和Pl,即：[0061]P3=P4+LNi[0化引[0化9][0060]联立上述3式，可求解出变量S和L。化Qo的单位向量为化，由此，可计算[0062][0063]又令P〇Pl=P3P4,则控制点Po可由下式导出：[0064][0065]由上述立式求得修正后的控制点?〇、&爪，控制点?5瓜爪同理可求得。[0066](b)考虑过渡段长度约束，采用比例调节算法对控制顶点进行一次性修正。在数控加工连续拐点路径中，如图3所示(虚线表示加工直线段，实线表示拐点处的平滑路径）：拐点化处两相邻线性加工直线段顶点分别为(Q〇，Qi)和(Qi，Q2)，拐点化处两相邻线性加工直线段顶点分别为(Qi，Q2)和(Q2，Q3)，讯^...：P8。巧ZQ0化化上的特征多边形的控制点，{妃.:.:.，转}为Z化化化上的特征多边形的控制点，为保证路径的连续可调并且不发生干设现象，即需要对原特征多边形进行约束，比例调节算法如下：[0067][006引[0069]k=max化i,k2)[0070]比例调节算法公式可W求出比例系数k，若k〉l，则需要对控制多边形等比例缩小，良P占此可由下列公式求得修正后的控制点Po'、Pi'、P3'：[0071]-[0072][0073][0074]同理，可求得控制点巧'、巧、巧'。修正之后，拐点处逼近误差Cl也相应变为：[0075][0076]上式表明修正之后的逼近误差Ei与最大逼近误差e相比缩小了k倍化〉1)。[0077](C)求出用于拐点平滑的五次B样条曲线。若，则直接将步骤(a)中求得的9个控制点带入五次B样条曲线公式中，得到五次B样条曲I则需要将步骤（b)中求得的9个控制点带入五次B样条曲线公式中，得到五次B样条曲线，即得到拐点处的平滑路径。[0078]图4、图5分别为采用内切式和本发明方法的拐点平滑路径示意图。内切式为采用3次Bgzier多项式曲线的方法，其中虚线表示原始的拐点路径，实线表示拐点处的平滑路径，拐点角度9=14.25°，最大逼近误差e=〇.15mm。图6、图7分别为采用内切式和本发明方法的拐点平滑路径曲率示意图。通过对比可W看出，在在满足同等精度的条件下，内切式的方法曲率最大值约为-284,采用本发明方法曲率最大值约为-15,采用本发明方法后，平滑路径的曲率缩小了近19倍，有更大曲率半径，所W采用本发明方法更适用于高速高性能加工领域。【主权项】1.一种离散加工路径的拐点平滑方法，其特征在于:首先，借助重顶点方法和凸包性质构造出由9个控制顶点组成的B样条曲线的特征多边形及其节点矢量的具体形式，并在逼近误差允值约束下反算出控制顶点;进而，基于过渡段长度约束，采用比例调节算法对控制顶点进行一次性修正;最后，根据求得的特征多边形和节点矢量得到最终用于拐点平滑的五次B样条曲线;采用的具体步骤为：(a)构造出由9个控制点组成的外切式B样条曲线的特征多边形及其节点矢量，并通过逼近误差允值约束来反算出控制点；其中，已知拐点Q:处两相邻线性加工直线段顶点分别为(Q^Qi)和(QiQd，拐点处的最大逼近误差允值为ε;拟构造的拐点过渡段的外切式五次B样条曲线表示为：式1中，Pi为特征多边形的控制点，Bi5为基函数，u为五次Β样条曲线的节点矢量，确定{Ρο,Ρι...，P8}为待求的特征多边形的9个控制点，且P4=Qi、Pl·=卩2、？6=?7，？()、？1和?2共线并位于直线段0(^1上，？6、？7和?8共线并位于直线段〇1〇2上，以保证加工直线段与拐点过渡曲线在公共连接点处G2连续;控制点p3、p5位于离散点所围三角形的外部，且P〇和Ps、Pi和P6、P3和P5分别关于角ZQoQiQ〗的平分线对称，点P3到线段QoQi和点P5到线段Q!Q2的距离均为最大逼近误差ε;由构造的特征多边形对称性和Β样条凸包的性质知，最大逼近误差的计算公式为：ε=||P4-C(0.5)|(2)进而，由式1得c(0.5)的表达式：C(0.5)=0.0185Pi+0.2593P3+0.4444P4+0.2593P5+0.0185Pe(3)整理式3后推出：P4_C(0·5)=0·0185P4Pi+0·2593P4P3+0·0185P4P6+0·2593P4P5(4)将式4带入到式2中，则最大逼近误差进一步表达为：(5)式5中，ZQoQiQ〗=θ，（Θ乒180°)，ZP;Ρ4Ρ5=?，P3P4=L;由于点P3到线段PiP4的距离为ε，且Ρ3在PiP4的垂直平分线上，则推出下列关系：联立式5、6和7,求解出变量0和UGhQo的单位向量为他，由此，计算出控制点P3和Ρι，即：P3=P4+LNi(8)(10)(9)又令PoPi=P3P4，则控制点Ρο由下式导出：由式8、9和10求得修正后的控制点Po、Pi、p3，同理求得控制点p5、p6、p8;(b)考虑过渡段长度约束，采用比例调节算法对控制顶点进行一次性修正；当数控加工连续拐点路径时，为保证路径的连续可调并且不发生干涉现象，即需要对原特征多边形进行约束，比例调节算法如下：k=max(ki，k2)(13)由式11、12、13求出比例系数k;若k>l，则需对控制多边形进行等比例缩小，即丨此由式14、15、16求得修正后的控制点Y3:P/3=P4+LiNi(14)同理，求得控制点6、^8;修正之后，拐点处逼近误差^也相应变为：a4(17)Κ式17表明修正之后的逼近误差￡1与最大逼近误差ε相比缩小了k倍(k>l);(c)求出用于拐点平滑的五次B样条曲线;若P3P4=L，则直接将步骤(a)中求得的9个控8制点带入式1中，得到五次B样条曲线〔'(11)=^?^诺则需将步骤(b)中求得/-0的9个控制点带入式1中，得到五次B样条曲线_，即得到拐点处的平滑路径。?=：0·【文档编号】G05B19/19GK105955194SQ201610308616【公开日】2016年9月21日【申请日】2016年5月10日【发明人】孙玉文,徐金亭,徐富阳,郭东明【申请人】大连理工大学