CN107885166B - 一种空间圆锥曲线通用插补方法 - Google Patents

一种空间圆锥曲线通用插补方法 Download PDF

Info

Publication number
CN107885166B
CN107885166B CN201711045309.2A CN201711045309A CN107885166B CN 107885166 B CN107885166 B CN 107885166B CN 201711045309 A CN201711045309 A CN 201711045309A CN 107885166 B CN107885166 B CN 107885166B
Authority
CN
China
Prior art keywords
interpolation
curve
point
conic section
end point
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
Application number
CN201711045309.2A
Other languages
English (en)
Other versions
CN107885166A (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.)
Beijing University of Aeronautics and Astronautics
Original Assignee
Beijing University of Aeronautics and Astronautics
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 Beijing University of Aeronautics and Astronautics filed Critical Beijing University of Aeronautics and Astronautics
Priority to CN201711045309.2A priority Critical patent/CN107885166B/zh
Publication of CN107885166A publication Critical patent/CN107885166A/zh
Application granted granted Critical
Publication of CN107885166B publication Critical patent/CN107885166B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/41—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 interpolation, e.g. the computation of intermediate points between programmed end points to define the path to be followed and the rate of travel along that path
    • G05B19/4103—Digital interpolation
    • 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/34—Director, elements to supervisory
    • G05B2219/34083—Interpolation general

Landscapes

  • Engineering & Computer Science (AREA)
  • Computing Systems (AREA)
  • Theoretical Computer Science (AREA)
  • Human Computer Interaction (AREA)
  • Manufacturing & Machinery (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Complex Calculations (AREA)
  • Image Generation (AREA)

Abstract

本发明公开了一种空间圆锥曲线通用插补方法，包括步骤一：设定空间圆锥曲线插补代码格式；步骤二：译码获取空间圆锥曲线数据信息；步骤三：获取圆锥曲线的常系数微分方程；步骤四：判断插补方向；步骤五：根据插补步长计算插补系数；步骤六：计算预插补点；步骤七：校正插补点；步骤八：终点判断。本发明对于各种圆锥曲线具有通用的NC代码格式，并将各种圆锥曲线转化为统一的常系数微分方程形式，采用相同的插补流程进行插补，完成任意空间平面的任意圆锥曲线弧的插补，本发明具有理论上的零径向插补误差，且通过步骤7参考点校正方式，大大减小了计算累积误差，能达到10‑4％以下的速度误差和极高的插补精度。

Description

一种空间圆锥曲线通用插补方法
技术领域
本发明涉及一种高精度空间圆锥曲线通用插补方法，属于数控加工技术领域。
背景技术
现有的数控系统一般只具备G01空间直线插补与G02/03平面圆弧插补，且插补圆弧时必须同时指定圆弧插补平面，即G17代码指定的XY平面、G18代码指定的XZ平面或G19代码指定的YZ平面。因此，空间自由曲线在加工之前一般都被离散成大量微小的空间直线段，形成大量连续的G01代码，数控系统根据G01代码使用空间直线插补方法完成空间自由曲线的加工。实际上，使用空间圆锥曲线(圆弧、椭圆、抛物线和双曲线)逼近空间自由曲线不仅能大量降低代码的数量、提高曲线逼近精度，同时能使得机床运动更加平稳，避免了直线段连接处的拐角突变。另外，空间圆锥曲线在光学元器件和椭圆齿轮等领域有重要应用，其直接插补算法能够很大程度上提高其加工质量。因此，需要开发一种插补精度高、计算速度快的空间圆锥曲线通用插补方法来实现各种形式的圆锥曲线的直接插补加工方式。
发明内容
本发明的目的是为了解决上述问题，提出了一种高精度空间圆锥曲线通用插补方法。该插补方法能根据给定的圆锥曲线的平面方向轴矢量和标准几何参数完成任意空间平面内的圆锥曲线插补。
一种空间圆锥曲线通用插补方法，包括以下几个步骤：
步骤一：设定空间圆锥曲线插补代码格式；
步骤二：译码获取空间圆锥曲线数据信息；
步骤三：获取圆锥曲线的常系数微分方程；
步骤四：判断插补方向；
步骤五：根据插补步长Li计算插补系数hi；
步骤六：计算预插补点；
步骤七：校正插补点；
步骤八：终点判断。
本发明的优点在于：
(1)通用性。对于各种圆锥曲线具有通用的NC代码格式，并将各种圆锥曲线转化为统一的常系数微分方程形式，采用相同的插补流程进行插补，完成任意空间平面的任意圆锥曲线弧的插补；
(2)高精度。该方法具有理论上的零径向插补误差，且通过步骤7参考点校正方式，大大减小了计算累积误差，能达到10-4％以下的速度误差和极高的插补精度；
(3)高效率。本发明在计算插补增量和插补点时除了5次开方计算外均为简单四则运算，避免了三角函数、双曲函数和坐标变换的计算，插补计算效率高；(4)简明性。本发明规定的G代码格式是标准的圆锥曲线几何表达方式，具有很强的几何意义，便于编程设计。
附图说明
图1是本发明中涉及到的圆锥曲线的定义；
图2是本发明插补方法的流程图；
图3是实施示例的插补结果图；
图4是实施示例插补过程中的径向相对误差图；
图5是实施示例插补过程中的速度相对误差曲线图。
缺少图2、图5的解释说明，请在具体实施方式中进行添加说明。另外，图2引用的时候，需要把其中的字母含义添加到相应地方。
具体实施方式
下面将结合附图和实施例对本发明作进一步的详细说明。
如图1所示，为表示一段圆锥曲线弧，需要其中心、起点、终点和方向矢量的信息：
其中Xs为圆锥曲线插补起点，Xe为圆锥曲线插补终点，Xc为圆锥曲线中心；
a,b分别为圆锥曲线长轴长与短轴长，U,V分别为圆锥曲线长轴方向向量与短轴方向向量；
Xs′为圆锥曲线插补起点的切矢，Xe′为圆锥曲线插补终点的切矢。
Xi,Xj为圆锥曲线上任意两点，Xi′,Xj′为对应的切矢。
下面对此种圆锥曲线弧的插补算法进行说明。
本发明是一种空间圆锥曲线通用插补方法流程如图2所示，包括以下几个步骤：
步骤一：设定空间圆锥曲线插补代码格式；
设定圆锥曲线代码格式：G02.*X(U)_Y(V)_Z(W)_I_J_K_UX_UY_UZ_AL_VX_VY_VZ_BL_F_
其中，G02.1为空间椭圆曲线(圆弧)，G02.2为空间双曲线，G02.3为空间抛物线曲线；X,Y,Z为圆锥曲线终点绝对坐标，U,V,W为圆锥曲线终点相对坐标(相对于起点)；I,J,K为圆锥曲线中心的相对坐标(相对于起点)；UX,UY,UZ为圆锥曲线平面的第一矢量，VX,VY,VZ为圆锥曲线对称轴矢量；AL为U方向长短轴的长度，BL为V方向长短轴的长度；F为插补进给速度，单位mm/min。
步骤二：译码获取空间圆锥曲线数据信息；
令圆锥曲线的起点为Xs＝(xs,ys,zs)，即前一段插补指令的终点，为已知；
若读取的圆锥曲线代码如下：
G02.*Xxe Yye Zze Iic Jjc Kkc UXux UYuy UZuz ALa VXvx VYvy VZvz BLb Ff，其中xe,ye,ze,ic,jc,kc,ux,uy为步骤一中定义的插补代码格式里的数值。
则对于椭圆有，终点Xe＝(xe,ye,ze)，中心Xc＝(xc,yc,zc)＝(xs+ic,ys+jc,zs+kc)，长轴方向单位矢量
Figure BDA0001452102500000031
短轴方向单位矢量
Figure BDA0001452102500000032
长轴长度为a，短轴长度为b，最大进给速度为
对于双曲线，终点Xe＝(xe,ye,ze)，中心Xc＝(xc,yc,zc)＝(xs+ic,ys+jc,zs+kc)，长轴方向单位矢量
Figure BDA0001452102500000034
短轴方向单位矢量
Figure BDA0001452102500000035
长轴长度为a，短轴长度为b，最大进给速度为
对于抛物线，终点Xe＝(xe,ye,ze)，中心Xc＝(xc,yc,zc)＝(xs+ic,ys+jc,zs+kc)，长轴方向单位矢量
Figure BDA0001452102500000037
短轴方向单位矢量通径2p＝b/a，最大进给速度为
Figure BDA0001452102500000039
步骤三：获取圆锥曲线的常系数微分方程；
用如下的参数化表示方法分别表示椭圆、双曲线和抛物线：
X＝(acosθ)U+(bsinθ)V+Xc
X＝(achθ)U+(bshθ)V+Xc
其中，X为圆锥曲线上的点，θ为圆锥曲线的参数，W为圆锥曲线平面的单位法线并与U、V组成右手坐标系。
由上述参数化形式可以得到圆锥曲线统一的常系数微分方程表达形式：
Figure BDA0001452102500000041
其中：A为3×3矩阵，微分方程的常系数矩阵。
步骤四：判断插补方向；
设X′s＝A(Xs-Xc),X′e＝A(Xe-Xc)，通过判断(Xe-Xs)与(X′e+X′s)的方向来确定插补方向：
当(Xe-Xs)·(X′e+X′s)≥0时插补参数增量为正，记flag＝1；
当(Xe-Xs)·(X′e+X′s)＜0时插补参数增量为负，记flag＝-1；
其中，flag为插补方向的标记，flag＝1时沿θ的正向插补，flag＝-1时沿θ的负向插补。
步骤五：根据插补步长计算插补系数；
对于第i个插补周期，根据规划的插补速度计算得到插补步长Li。
利用二阶泰勒展开估计圆锥曲线的参数θ的增量
Figure BDA0001452102500000042
X′i+1≈X′i+X″iθi+X″′i(Δθi)2＝X′i+AX′iΔθi+A2X′i(Δθi)2
其中，Δθi为第i个周期里参数θ的增量估计，Xi为第i个插补周期的起点，Xi+1为第i个周期的插补点，Xi+1′为曲线在Xi+1处的切矢，Xi′,Xi″,Xi″′为曲线在Xi处的一阶、二阶和三阶切矢。
可以证明，在步骤二给定的参数化下，三种圆锥曲线均有Xi+1-Xi//X′i+1+X′i，令第i个周期的插补系数
Figure BDA0001452102500000043
则有Xi+1-Xi＝hi(X′i+1+X′i)。
步骤六：计算预插补点；
对第i个周期，根据步骤五中当前插补点和下一插补点的递增关系和步骤中的微分表达形式，得到
(E-hiA)(Xi+1-Xc)＝(E+hiA)(Xi-Xc)，其中E为单位矩阵
当(E-hiA)可逆时有
Figure BDA0001452102500000051
其中
Figure BDA0001452102500000052
为预插补点。
步骤七：校正插补点；
预设插补参考点XR＝Xs；
根据步骤六求得的预插补点
Figure BDA0001452102500000053
和参考点XR计算参考插补系数
Figure BDA0001452102500000054
利用方程Xi+1＝(E-hRA)-1(E+hRA)(XR-Xc)+Xc，计算最终插补点Xi+1；
当|hR|＞0.9时，更新XR＝Xi+1，以保证(E-hiA)可逆和方程求解的数值稳定。
步骤八：终点判断；
求剩余插补投影距离
Figure BDA0001452102500000055
若则认为插补结束；否则，继续向前插补。
根据以上步骤，本发明的一种空间圆锥曲线插补方法适用于任意起终点、任意空间平面的连续圆锥曲线段的插补。插补方向由起终点位置得出，对于椭圆弧和圆弧为劣弧方向，双曲线和抛物线为起点至终点方向。
实施例：
以一个半椭圆弧为插补对象，实施流程图如图1所示，有以下具体步骤：
以一个长轴长为200，短轴长为100，长轴方向方向量为(3,4,0)，短轴方向向量为(-4,3,0)，椭圆中心在原点(0,0,0)的椭圆弧为插补对象。定义椭圆弧插补起点为(120,160,0)，即长轴的一个端点，终点为(-143.9230,-108.5641,0)，椭圆的标准表达式为
Figure BDA0001452102500000056
进给速度为6000mm/min。
步骤一：定义空间圆锥曲线插补代码格式：
定义圆锥曲线(椭圆弧)代码格式：G02.1 X-143.9230 Y-108.5641 Z0 I-120 J-160 K0 UX3 UY4 UZ0 AL200 VX-4 VY3 VZ0 BL100 F6000
步骤二：译码获取空间圆锥曲线数据信息：
起点为Xs＝(120,160,0)；
椭圆插补时，终点为Xe＝(-143.9230,-108.5641,0)，中心为Xc＝(0,0,0)，长轴方向单位矢量U＝(0.6,0.8,0)，短轴方向单位矢量V＝(-0.8,0.6,)0，长轴长度为a＝200，短轴长度为b＝100，最大进给速度为100mm/min；
步骤三：圆锥曲线的常系数微分方程表示：
圆锥曲线统一的常系数微分方程表达形式：
Figure BDA0001452102500000061
其中
Figure BDA0001452102500000062
步骤四：判断插补方向:
X′s＝A(Xs-Xc)＝(-80,60,0),X′e＝A(Xe-Xc)＝(9.2820,-131.9615,0)
(Xe-Xs)·(X′e+X′s)＝(-263.9230,-268.5641,0)·(-70.7180,-71.9615,0)＝37990＞0，故插补参数增量为正，记flag＝1；
步骤五：根据插补步长计算插补系数：
设插补周期Ts＝0.001s，则插补步长Li＝0.1mm
利用二阶泰勒展开估计第1个插补参数增量，当前点为起点，则
X1＝Xs＝(120,160,0)
X′1＝A(X1-Xc)＝(-80,60,0)
X″1＝AX′1＝(-120,-160,0)
X″′1＝AX″1＝(80,-60,0)
Figure BDA0001452102500000063
X′2≈X′1+X″1Δθ1+X″′1(Δθ1)2＝(-80.1199,59.8399,0)
步骤六：计算预插补点：
对第1个周期
Figure BDA0001452102500000065
步骤七：校正终插补点：
预设插补参考点XR＝Xs；
根据步骤6求得的
Figure BDA0001452102500000066
和参考点XR计算
Figure BDA0001452102500000067
利用方程X2＝(E-hRA)-1(E+hRA)(XR-Xc)+Xc计算得到X2＝(119.9199,160.0599,0)作为新的插补点；
因|hR|＜0.9时，故不更新XR。
步骤八：终点判断：
求当前插补点距终点的距离
Figure BDA0001452102500000071
继续向前插补。
反复执行步骤5至步骤8，直至插补到终点，得到的插补点如下：
插补周期i X Y Z
0 120 160 0
1 119.919940040045 160.059919970060 0
2 119.839760400957 160.119679701282 0
3 119.759461405061 160.179278956798 0
4 119.679043376215 160.238717501838 0
5 119.598506639798 160.297995103740 0
6 119.517851522699 160.357111531961 0
7 119.437078353299 160.416066558082 0
8 119.356187461461 160.474859955821 0
9 119.275179178518 160.533491501036 0
10 119.194053837253 160.591960971742 0
11 119.112811771893 160.650268148109 0
12 119.031453318090 160.708412812478 0
13 118.949978812910 160.766394749368 0
14 118.868388594817 160.824213745479 0
15 118.786683003659 160.881869589705 0
16 118.704862380657 160.939362073138 0
17 118.622927068386 160.996690989077 0
…… …… …… ……
上表中(x,y,z)表示第i个插补周期生成的插补点。
图3所示为插补结果，可以看出，本实例为一个在XY平面内的椭圆；图4为插补点的径向相对误差，均在10-13％量级；图5所示为插补过程中产生的相对速度误差曲线，可以看出，本发明方法产生的速度误差在0.075％以下。由插补结果可以得出，本发明提出的插补算法符合插补要求。

Claims (6)

1.一种空间圆锥曲线通用插补方法，包括以下几个步骤：
步骤一：设定空间圆锥曲线插补代码格式；
步骤二：译码获取空间圆锥曲线数据信息；
令圆锥曲线的起点为Xs＝(xs,ys,zs)，即前一段插补指令的终点，为已知；
若读取的圆锥曲线代码如下：
G02.*Xxe Yye Zze Iic Jjc Kkc UXux UYuy UZuz ALa VXvx VYvy VZvz BLb Ff；
则对于椭圆有，终点Xe＝(xe,ye,ze)，中心Xc＝(xc,yc,zc)＝(xs+ic,ys+jc,zs+kc)，长轴方向单位矢量
Figure FDA0002285250490000011
短轴方向单位矢量
Figure FDA0002285250490000012
长轴长度为a，短轴长度为b，最大进给速度为
Figure FDA0002285250490000013
对于双曲线，终点Xe＝(xe,ye,ze)，中心Xc＝(xc,yc,zc)＝(xs+ic,ys+jc,zs+kc)，长轴方向单位矢量
Figure FDA0002285250490000014
短轴方向单位矢量
Figure FDA0002285250490000015
长轴长度为a，短轴长度为b，最大进给速度为
Figure FDA0002285250490000016
对于抛物线，终点Xe＝(xe,ye,ze)，中心Xc＝(xc,yc,zc)＝(xs+ic,ys+jc,zs+kc)，长轴方向单位矢量
Figure FDA0002285250490000017
短轴方向单位矢量通径2p＝b/a，最大进给速度为
Figure FDA0002285250490000019
步骤三：获取圆锥曲线的常系数微分方程；
用如下的参数化表示方法分别表示椭圆、双曲线和抛物线：
X＝(acosθ)U+(bsinθ)V+Xc
X＝(achθ)U+(bshθ)V+Xc
Figure FDA00022852504900000110
其中，X为圆锥曲线上的点，θ为圆锥曲线的参数，W为圆锥曲线平面的单位法线并与U、V组成右手坐标系；
由上述参数化形式得到圆锥曲线统一的常系数微分方程表达形式：
Figure FDA0002285250490000021
其中
其中：A为3×3矩阵，微分方程的常系数矩阵；
步骤四：判断插补方向；
设X′s＝A(Xs-Xc),X′e＝A(Xe-Xc)，通过判断(Xe-Xs)与(X′e+X′s)的方向来确定插补方向：
当(Xe-Xs)·(X′e+X′s)≥0时插补参数增量为正，记flag＝1；
当(Xe-Xs)·(X′e+X′s)＜0时插补参数增量为负，记flag＝-1；
其中：flag为插补方向的标记，flag＝1时沿θ的正向插补，flag＝-1时沿θ的负向插补；Xs为圆锥曲线插补起点，Xe为圆锥曲线插补终点，Xc为圆锥曲线中心；Xs′为圆锥曲线插补起点的切矢，Xe′为圆锥曲线插补终点的切矢；
步骤五：根据插补步长计算插补系数；
步骤六：计算预插补点；
步骤七：校正插补点；
步骤八：终点判断。
2.根据权利要求1所述的一种空间圆锥曲线通用插补方法，所述的步骤一中：
设定圆锥曲线代码格式：G02.*X(U)_Y(V)_Z(W)_I_J_K_UX_UY_UZ_AL_VX_VY_VZ_BL_F_；*取值为1,2,3；
其中，G02.1为空间椭圆曲线，G02.2为空间双曲线，G02.3为空间抛物线曲线；X,Y,Z为圆锥曲线终点绝对坐标，X(U),Y(V),Z(W)为圆锥曲线终点相对坐标；I,J,K为圆锥曲线中心的相对坐标；UX,UY,UZ为圆锥曲线平面的第一矢量，VX,VY,VZ为圆锥曲线对称轴矢量；AL为U方向长短轴的长度，BL为V方向长短轴的长度；F为插补进给速度，单位mm/min。
3.根据权利要求1所述的一种空间圆锥曲线通用插补方法，所述的步骤五中：
对于第i个插补周期，根据规划的插补速度计算得到插补步长Li；
利用二阶泰勒展开估计插补参数增量：
Figure FDA0002285250490000031
X′i+1≈X′i+X″iΔθi+X″′i(Δθi)2＝X′i+AX′iΔθi+A2X′i(Δθi)2
其中，Δθi为第i个周期里参数θ的增量估计，Xi为第i个插补周期的起点，Xi+1为第i个周期的插补点，Xi+1′为曲线在Xi+1处的切矢，Xi′,Xi″,Xi″′为曲线在Xi处的一阶、二阶和三阶切矢；
令第i个周期的插补系数
4.根据权利要求1所述的一种空间圆锥曲线通用插补方法，所述的步骤六中：
对第i个周期，当前插补点和下一插补点有如下关系：
(E-hiA)(Xi+1-Xc)＝(E+hiA)(Xi-Xc)，其中E为单位矩阵；
当(E-hiA)可逆时有
Figure FDA0002285250490000033
其中
Figure FDA0002285250490000034
为预插补点。
5.根据权利要求1所述的一种空间圆锥曲线通用插补方法，所述的步骤七中：
预设插补参考点XR＝Xs；
根据步骤六求得的预插补点
Figure FDA0002285250490000035
和参考点XR计算参考插补系数
Figure FDA0002285250490000036
利用方程Xi+1＝(E-hRA)-1(E+hRA)(XR-Xc)+Xc，计算最终插补点Xi+1；
当|hR|＞0.9时，更新XR＝Xi+1。
6.根据权利要求1所述的一种空间圆锥曲线通用插补方法，所述的步骤八中：
求剩余插补投影距离
Figure FDA0002285250490000037
若
Figure FDA0002285250490000038
则认为插补结束；否则，继续向前插补；
其中，Li为插补步长，Xe′为圆锥曲线插补终点的切矢。
CN201711045309.2A 2017-10-31 2017-10-31 一种空间圆锥曲线通用插补方法 Expired - Fee Related CN107885166B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711045309.2A CN107885166B (zh) 2017-10-31 2017-10-31 一种空间圆锥曲线通用插补方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711045309.2A CN107885166B (zh) 2017-10-31 2017-10-31 一种空间圆锥曲线通用插补方法

Publications (2)

Publication Number Publication Date
CN107885166A CN107885166A (zh) 2018-04-06
CN107885166B true CN107885166B (zh) 2020-01-21

Family

ID=61783080

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711045309.2A Expired - Fee Related CN107885166B (zh) 2017-10-31 2017-10-31 一种空间圆锥曲线通用插补方法

Country Status (1)

Country Link
CN (1) CN107885166B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113467368A (zh) * 2021-07-15 2021-10-01 苏州谋迅智能科技有限公司 S形速度曲线的调整方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2011108902A (ru) * 2011-03-09 2012-09-20 Государственное образовательное учреждение высшего профессионального образования Томский государственный университет систем управл Способ обнаружения и селекции радиолокационных сигналов по поляризационному признаку и устройство для его осуществления
CN103149879A (zh) * 2011-12-07 2013-06-12 沈阳高精数控技术有限公司 一种基于弧长的数控系统椭圆插补方法
CN103728980A (zh) * 2014-01-08 2014-04-16 哈尔滨工业大学 航天器相对轨道的控制方法
CN104133423A (zh) * 2014-07-16 2014-11-05 北京航空航天大学 一种空间椭圆弧插补方法
KR101576426B1 (ko) * 2015-06-24 2015-12-11 중앙대학교 산학협력단 어안 렌즈를 이용한 감시 장치 및 방법
CN105676787A (zh) * 2015-12-28 2016-06-15 龙兵 一种椭圆弧插补算法
CN106950920A (zh) * 2017-04-18 2017-07-14 大连奥托股份有限公司 基于数控类设备的空间圆弧插补方法
CN107102617A (zh) * 2017-06-26 2017-08-29 北京艾利特科技有限公司 一种高精度空间椭圆曲线实时插补方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2011108902A (ru) * 2011-03-09 2012-09-20 Государственное образовательное учреждение высшего профессионального образования Томский государственный университет систем управл Способ обнаружения и селекции радиолокационных сигналов по поляризационному признаку и устройство для его осуществления
CN103149879A (zh) * 2011-12-07 2013-06-12 沈阳高精数控技术有限公司 一种基于弧长的数控系统椭圆插补方法
CN103728980A (zh) * 2014-01-08 2014-04-16 哈尔滨工业大学 航天器相对轨道的控制方法
CN104133423A (zh) * 2014-07-16 2014-11-05 北京航空航天大学 一种空间椭圆弧插补方法
KR101576426B1 (ko) * 2015-06-24 2015-12-11 중앙대학교 산학협력단 어안 렌즈를 이용한 감시 장치 및 방법
CN105676787A (zh) * 2015-12-28 2016-06-15 龙兵 一种椭圆弧插补算法
CN106950920A (zh) * 2017-04-18 2017-07-14 大连奥托股份有限公司 基于数控类设备的空间圆弧插补方法
CN107102617A (zh) * 2017-06-26 2017-08-29 北京艾利特科技有限公司 一种高精度空间椭圆曲线实时插补方法

Also Published As

Publication number Publication date
CN107885166A (zh) 2018-04-06

Similar Documents

Publication Publication Date Title
Feng et al. Constant scallop-height tool path generation for three-axis sculptured surface machining
US4648024A (en) Curvilinear interpolation system and method
Bi et al. A general, fast and robust B-spline fitting scheme for micro-line tool path under chord error constraint
CN108170101B (zh) 面向多项式样条曲线的插补方法及系统
CN108073138B (zh) 适用于高速高精加工的椭圆弧平滑压缩插补算法
CN107102617B (zh) 一种高精度空间椭圆曲线实时插补方法
CN104155916B (zh) 一种高精度快速空间圆弧插补方法
CN109597359B (zh) 数值控制装置
CN112109084A (zh) 基于机器人关节角度补偿的末端位置补偿方法及其应用
US20020049514A1 (en) Method for compensating for static position errors and orientation errors
CN107885166B (zh) 一种空间圆锥曲线通用插补方法
Li et al. The design of a NURBS pre-interpolator for five-axis machining
JP2015121966A (ja) テーブル形式データによる運転でのスムージング機能を備えた数値制御装置
Liu et al. Weld seam fitting and welding torch trajectory planning based on NURBS in intersecting curve welding
Yang et al. A novel contouring error estimation method for contouring control
EP3776111B1 (en) System and method for generating optimal lattice tool paths
US4581698A (en) Method and system for generating interpolation pulses
CN113111404B (zh) 一种加工轨迹的空间连续小线段圆弧与直线拟合方法
CN114035512B (zh) 一种基于两次泰勒展开的最小速度波动插补方法
CN106292531B (zh) 一种计算加工zn1蜗杆盘状成形刀具廓形边界的算法
CN116551457B (zh) 一种通过样条曲线连接加工直线段与圆弧段的方法与设备
CN111444250A (zh) 一种弯管机cad图纸导入方法
Zhao et al. A novel trajectory interpolation algorithm for wedm-unit generalized arc length increment method
CN113515091B (zh) 一种拐角平滑样条的轨迹实时插补方法和装置
JPS5822410A (ja) 数値制御方式

Legal Events

Date Code Title Description
PB01 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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20200121

Termination date: 20201031