五轴加工走刀轨迹奇异点避免方法技术

本发明专利技术公开了一种五轴加工走刀轨迹奇异点避免方法，用于解决现有五轴加工走刀轨迹奇异点避免方法加工精度低的技术问题。技术方案是使用四元数表示刀轴方向与参考方向之间的变换关系，再将四元数投影到P平面，“奇异锥”在P平面内对应于一个圆心在原点的圆“奇异圆”；在二维平面中使用B样条对投影点进行插值；再次建立约束优化问题，目标函数定义为控制顶点增量的平方和，约束定义为B样条上的点到P平面原点的距离大于“奇异圆”半径，控制顶点的初值定义为插值得到的B样条控制顶点；将约束优化问题的约束用B样条控制顶点的增量线性表示，使用二次规划理论求解优化问题的最优解。本发明专利技术将一般约束优化问题转换为二次规划问题，保证了加工精度。

【技术实现步骤摘要】
五轴加工走刀轨迹奇异点避免方法
本专利技术属于数控加工领域，特别涉及一种五轴加工走刀轨迹奇异点避免方法。
技术介绍
文献1“A.Affouard,E.Duc,C.Lartigue,J.M.Langeron,P.Bourdet,Avoiding5-axissingularitiesusingtoolpathdeformation,InternationalJournalofMachineToolsandManufacture44(2004)415-425.”公开了一种适用于对偶B样条路径的五轴数控机床奇异点避免方法，提出使用“奇异锥”的概念来代替奇异点的概念，同时通过修改B样条的控制顶点避免刀轴方向进入“奇异锥”。但这种刀具路径变形无法保证刀具路径变形引起的加工误差最小。文献2“J.X.Yang,Y.Altintas,Generalizedkinematicsoffive-axisserialmachineswithnonsingulartoolpathgeneration,InternationalJournalofMachineToolsandManufacture75(2013)119-132.”公开了一种使用四元数避免五轴数控机床奇异点的方法，该方法接受了“奇异锥”的概念，通过将刀轴方向样条进入“奇异锥”部分的控制顶点旋转一个“奇异锥”的锥角来避免刀轴方向进入“奇异锥”。但是该方法对刀轴方向样条的修改过大，会造成不必要的加工精度损失。以上文献的典型特点是：均对奇异点附近的刀轴方向进行了修改且均不能保证由刀轴方向变形引起的加工误差最小，会造成加工精度不必要的损失。
技术实现思路
为了克服现有五轴加工走刀轨迹奇异点避免方法加工精度低的不足，本专利技术提供一种五轴加工走刀轨迹奇异点避免方法。该方法首先使用四元数表示刀轴方向与参考方向之间的变换关系，再将得到的四元数投影到一个二维平面(P平面)，“奇异锥”在P平面内对应于一个圆心在原点的圆“奇异圆”；其次，在二维平面中使用B样条对投影点进行插值；再次建立约束优化问题，目标函数定义为控制顶点增量的平方和，约束定义为B样条上的点到P平面原点的距离大于“奇异圆”半径，控制顶点的初值定义为插值得到的B样条控制顶点；最后，将约束优化问题的约束用B样条控制顶点的增量线性表示，使用二次规划理论求解优化问题的最优解。本专利技术将一般约束优化问题转换为正定的二次规划问题，保证了优化问题有且仅有一个全局最优解。故本专利技术能够保证避免奇异点后得到的路径的变形最小，路径变形导致的加工误差最小，解决了现有的奇异点避免方法带来的较大的加工精度损失的问题。本专利技术解决其技术问题所采用的技术方案：一种五轴加工走刀轨迹奇异点避免方法，其特点是包括以下步骤：(1)选择参考方向or为机床运动链中第一个转轴坐标轴的方向，刀轴方向oi相对于参考方向的关系用四元数qu,i表示为：i＝0,1,…,N，N表示刀位点的个数，*表示四元数的共轭。四元数qu,i计算如下：(2)根据步骤(1)得到的四元数qu,i的特点，将其投影到P平面，四元数qu,i对应于P平面上的点pi，奇异锥对应于P平面上一个半径的奇异圆，α表示奇异锥的半锥角。(3)使用B样条对点集{pi}，i＝0,1,…,N，进行插值，节点矢量U＝{u0,u1,...,uN+k,uN+k+1}计算如下：B样条基函数计算如下：k表示B样条的阶数。控制顶点{Pi}，i＝0,1,…,N通过求解如下(N+1)×(N+1)的线性方程组得到：(4)考虑几何约束奇异圆对步骤(3)得到的B样条进行优化。目标函数定义为B样条控制顶点增量的平方和，约束条件定义为B样条上点到P平面原点的距离大于奇异圆半径R的c倍,c为一个大于1的常数。ΔPi为控制顶点Pi的增量，(xj,yj)为B样条上一点,i＝0,1,…N。m为B样条上检查点的个数，所有检查点处的约束共同组成了优化问题的约束优化后的B样条为(5)忽略ΔPi的高阶小项，将步骤(4)得到的约束条件用控制顶点的增量线性表示为，AΔP<B其中，和Nj＝[N0,k(uj),N1,k(uj),…,NN,k(uj)]N+1，j＝1,…,m表示步骤(3)插值得到的控制顶点向量和B样条基函数向量。(6)对步骤(5)的优化过程进行迭代，优化函数及约束条件整理为，AiterΔPiter<Biter其中，Piter-1表示上一次迭代得到的B样条控制顶点，ΔPiter表示当前迭代过程的控制顶点增量。迭代结束后得到最优的避免奇异点后的B样条控制顶点。(7)根据步骤(6)得到的P平面上的B样条，四元数q(u)计算如下，p＝(x,y)为P平面上一点，q为p对应的四元数，inf表示P平面上的无穷远点。根据公式o＝qorq*得到避免奇异点后的刀轴方向o。本专利技术的有益效果是：该方法首先使用四元数表示刀轴方向与参考方向之间的变换关系，再将得到的四元数投影到一个二维平面(P平面)，“奇异锥”在P平面内对应于一个圆心在原点的圆“奇异圆”；其次，在二维平面中使用B样条对投影点进行插值；再次建立约束优化问题，目标函数定义为控制顶点增量的平方和，约束定义为B样条上的点到P平面原点的距离大于“奇异圆”半径，控制顶点的初值定义为插值得到的B样条控制顶点；最后，将约束优化问题的约束用B样条控制顶点的增量线性表示，使用二次规划理论求解优化问题的最优解。本专利技术将一般约束优化问题转换为正定的二次规划问题，保证了优化问题有且仅有一个全局最优解。因此，本专利技术能够保证避免奇异点后得到的路径的变形最小，路径变形导致的加工误差最小，解决了现有的奇异点避免方法带来的较大的加工精度损失的问题。下面结合附图和具体实施方式对本专利技术作详细说明。附图说明图1是本专利技术五轴加工走刀轨迹奇异点避免方法与
技术介绍
方法在避免奇异点后的结果对比图。具体实施方式参照图1。本专利技术五轴加工走刀轨迹奇异点避免方法具体步骤如下：(1)选择数控机床结构为AC摇篮式五轴数控机床，下表给出的刀轴方向包含奇异点，用于验证本专利技术得到解的最优性及有效性。oiojok-0.212570.5174080.828918-0.104970.5519950.8272140.0393480.5524760.83260.1745320.499170.8487450.2296760.4508210.862560.2410740.4110430.8791620.2360420.3681780.8992940.2203080.3320030.9171910.1921370.2982750.9349410.1139860.231680.9660910.0207980.1281320.991539-0.004340.0851680.996357-0.013570.0479880.998756-0.012670.0221780.9996740010.02313-0.010380.9996790.061559-0.009760.9980560.1319440.0175070.9911030.303070.0717380.9502640.3986340.0826440.9133790.4854680.0590770.8722560.549128-0.0255本文档来自技高网...

【技术保护点】
一种五轴加工走刀轨迹奇异点避免方法，其特征在于包括以下步骤：(1)选择参考方向o

【技术特征摘要】
1.一种五轴加工走刀轨迹奇异点避免方法，其特征在于包括以下步骤：(1)选择参考方向or为机床运动链中第一个转轴坐标轴的方向，刀轴方向oi相对于参考方向的关系用四元数qu,i表示为：i＝0,1,…,N，N表示刀位点的个数，*表示四元数的共轭；四元数qu,i计算如下：(2)根据步骤(1)得到的四元数qu,i的特点，将其投影到P平面，四元数qu,i对应于P平面上的点pi，奇异锥对应于P平面上一个半径的奇异圆，α表示奇异锥的半锥角；(3)使用B样条对点集{pi}，i＝0,1,…,N，进行插值，节点矢量U＝{u0,u1,...,uN+k,uN+k+1}计算如下：B样条基函数计算如下：k表示B样条的阶数；控制顶点{Pi}，i＝0,1,…,N通过求解如下(N+1)×(N+1)的线性方程组得到：(4)考虑几何约束...

【专利技术属性】
技术研发人员：万敏，刘洋，张卫红，邢婉静，
申请(专利权)人：西北工业大学，
类型：发明
国别省市：陕西,61