CN105843174B - 一种样条曲线插补参数计算方法 - Google Patents

一种样条曲线插补参数计算方法 Download PDF

Info

Publication number
CN105843174B
CN105843174B CN201610165047.2A CN201610165047A CN105843174B CN 105843174 B CN105843174 B CN 105843174B CN 201610165047 A CN201610165047 A CN 201610165047A CN 105843174 B CN105843174 B CN 105843174B
Authority
CN
China
Prior art keywords
interpolation
parameter
value
curve
period
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.)
Active
Application number
CN201610165047.2A
Other languages
English (en)
Other versions
CN105843174A (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.)
Dalian University of Technology
Original Assignee
Dalian University of Technology
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 Dalian University of Technology filed Critical Dalian University of Technology
Priority to CN201610165047.2A priority Critical patent/CN105843174B/zh
Publication of CN105843174A publication Critical patent/CN105843174A/zh
Application granted granted Critical
Publication of CN105843174B publication Critical patent/CN105843174B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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

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)
  • Numerical Control (AREA)

Abstract

本发明一种样条曲线插补参数计算方法属于精密高效数控加工技术领域，特别涉及一种样条曲线插补过程中减小速度波动的插补点曲线参数计算方法。该方法首先将各插补周期的指令参数值写成关于时间的函数，并利用二阶Runge‑Kutta法计算下一插补周期的指令参数初始值；其次，通过比较利用参数初始值计算得到的进给速度与理想进给速度，计算参数补偿值；最后，根据参数初始值和参数补偿值，确定下一插补点曲线参数值，实现速度波动小的样条曲线实时插补；本发明可实现在不进行迭代计算的前提下有效抑制曲线插补过程中的速度波动，对提高样条曲线插补加工质量具有实用价值。

Description

一种样条曲线插补参数计算方法
技术领域
本发明属于精密高效数控加工技术领域，特别涉及一种样条曲线插补过程中减小速度波动的插补点曲线参数计算方法。
背景技术
具有复杂面形结构的复杂曲面零件在航空航天、能源动力等领域应用越来越广。随着高端装备领域迅速发展，制造业对该类复杂曲面零件的加工质量、加工效率要求不断提高。鉴于复杂曲面零件加工刀轨多为具有复杂几何特征的曲线，采用传统的直线、圆弧插补数控加工时，存在诸多缺陷，如利用微小直线、圆弧段代替曲线时产生的逼近误差，刀位文件庞大，存在一阶不连续点导致的频繁加减速等，不利于该类复杂曲面零件的精密高效加工。因此，随着计算机辅助制造技术的发展，参数样条曲线插补技术可克服直线、圆弧插补的上述缺陷，故得到了广泛关注。在样条曲线插补过程中，由于曲线参数和曲线弧长具有非线性对应关系，故插补周期内理想运行距离对应的曲线参数精确增量值获取困难，从而导致实际计算参数对应的插补点与理想插补点不一致，诱发实际进给速度与理想进给速度产生偏差，即速度波动。因此，研究样条曲线插补过程中，速度波动小的插补参数计算方法对实现精密、可靠的曲线插补数控加工具有重要意义。
文献“Parametric interpolator versus linear interpolator for precisionCNC machining”，Yang等，Computer-Aided Design，1994，26(3)：225-234，该文献提出一种二阶泰勒技术展开参数计算方法，然而，该方法需要计算曲线的二阶导失，且存在较大舍入误差；文献“A real-time predictor-corrector interpolator for CNC machining”，Tsai等，Journal of Manufacturing Science and Engineering，2003，125(3)：449-460，该文献提出一种基于参数预估-校正的闭环曲线插补参数计算方法，然而，该方法计算曲线参数时，需要进行迭代循环，直到满足精度条件时才能输出参数值，计算时间无法预测，不利于实时插补。
发明内容
本发明旨在克服现有技术缺陷，发明一种样条曲线插补参数计算方法，该方法基于二阶Runge-Kutta法和参数补偿的样条曲线插补参数计算方法，在避免迭代计算、避免大幅增大计算量的情况下，有效抑制曲线插补过程中的进给速度波动，提高曲线插补质量。
本发明的技术方案是一种样条曲线插补参数计算方法，其特性在于，该方法，首先将各插补周期的指令参数值写成关于时间的函数，并利用二阶Runge-Kutta法计算下一插补周期的指令参数初始值；其次，通过比较利用参数初始值计算得到的进给速度与理想进给速度，计算参数补偿值；最后，根据参数初始值和参数补偿值，确定下一插补点曲线参数值，实现速度波动小的样条曲线实时插补；方法具体步骤如下：
第一步计算插补参数初始值
将插补参数u写成时间t的函数：
u＝u(t) (1)
则参数u对时间t的导数为：
其中，s为曲线弧长，ti为第i个插补周期时的总插补时间，ui为第i个插补周期的指令参数值；令第i个插补周期的理想进给速度为vi，则：
设样条曲线参数方程为C＝C(u)，则：
其中，C′(ui)为曲线在参数u＝ui位置处的一阶导失；将公式(3)、(4)代入公式(2)得：
因此，利用二阶Runge-Kutta法展开公式(1)得到下一个，即第(i+1)个插补周期的指令插补参数初始值us,i+1为：
其中，Ts为插补周期，
第二步计算插补参数补偿值
为保证实际插补进给速度和理想进给速度相等，第(i+1)个插补周期的指令参数值ui+1应满足：
该插补点的插补参数补偿值Δui+1应满足：
Δui+1＝ui+1-us,i+1 (8)
根据一阶泰勒级数展开，可得到：
C(ui+1)＝C(us,i+1)+C′(us,i+1)Δui+1 (9)
将公式(9)代入公式(7)得：
整理得：
其中，A＝||C′(us,i+1)||2，B＝2(C′(us,i+1),(C(us,i+1)-C(ui))),
解方程(11)得插补参数的两个根Δui+1,1和Δui+1,2分别为：
由于因此Δui+1,1≈0，Δui+1,2≈-B/A；因此，为保证插补过程的稳定性，当方程(11)有实数根，即判别式B2-4AD≥0时，取绝对值较小的Δui+1,1作为插补参数补偿值Δui+1；当方程(11)无实数根，即判别式B2-4AD<0时，令参数补偿值Δui+1为零；因此，参数补偿值Δui+1的计算方法为：
第三步根据插补参数初始值us,i+1和插补参数补偿值Δui+1确定第(i+1)个插补周期的插补参数值ui+1：
ui+1＝us,i+1+Δui+1 (15)
利用该参数值作为数控系统插补过程中第i个插补周期的输入指令进行实际插补加工；判断是否到达曲线终点，若到达终点，则结束插补，否则，令i＝i+1，返回第一步。
本发明的有益效果是：建立了基于二阶Runge-Kutta法和参数补偿的样条曲线插补参数计算方法，可有效抑制曲线插补过程中的进给速度波动；参数计算过程中无需迭代循环，对高速精密实时参数曲线插补具有重要意义。
附图说明
图1——计算方法整体流程图；
图2——“倒尖八”形非均匀有理B样条曲线几何模型图；
图3——采用本计算方法的速度波动图；其中，A轴表示样条曲线参数，B轴表示进给速度波动值，单位为％；
图4——采用二阶泰勒级数展开法插补参数计算方法的速度波动图；其中，A轴表示样条曲线参数，B轴表示进给速度波动值，单位为％；
具体实施方式
结合技术方案与附图详细说明本发明的具体实施方式。
参数样条曲线插补过程中，由于曲线参数与弧长具有非线性对应关系，若不采用耗时无法确定的迭代方法计算插补参数，极易引起实际进给速度与理想速度产生偏差，即速度波动，影响零件加工表面质量。据此，发明一种样条曲线插补参数计算方法。
以“倒尖八”形非均匀有理B样条曲线插补为例，详细说明本发明实施过程，整体流程参见附图1。“倒尖八”形非均匀有理B样条曲线参数为：阶数：2；控制点：{(0,0)；(-50,-50)；(-50,50)；(0,0)；(50,-50)；(50,50)；(0,0)}；权因子：{5,5,10,1,10,5,5}；节点向量：{0,0,0,0.25,0.5,0.5,0.75,1,1,1}；其几何模型如附图2所示；
附图1是方法整体流程图，如附图1所示，实施例的具体步骤如下：
第一步计算插补参数初始值：首先，令u1＝0，i＝1，i为当前插补点序号；采用公式(6)计算下一插补点即第(i+1)个插补点的样条曲线插补参数初始值us,i+1；本实例中，插补周期Ts＝0.002s，理想进给速度vi恒为100mm/s；
第二步计算插补参数补偿值：将第一步中计算得到的插补参数初始值us,i+1代入公式(14)中，计算该插补点处的指令插补参数补偿值Δui+1；
第三步根据插补参数初始值us,i+1和插补参数补偿值Δui+1确定第(i+1)个插补周期的插补参数值ui+1：将第一步中计算得到的插补参数初始值us,i+1和第二步中计算得到的插补参数补偿值Δui+1代入公式(15)，从而确定下一插补点曲线参数值ui+1，并利用该参数值作为数控系统输入指令进行实际插补加工；判断是否到达曲线终点，若到达终点，则结束插补，否则，令i＝i+1，返回第一步。
利用上述方法进行曲线插补加工，计算得到曲线各位置处实际插补进给速度与理想进给速度的偏差，即速度波动值如附图3所示，可见，最大进给速度波动值为0.00000388％，几乎将速度波动完全消除；为说明本发明实施效果，将本发明方法与传统二阶泰勒级数展开法的速度波动值进行对比，采用传统二阶泰勒级数展开法进行曲线插补加工，得到曲线上各位置处速度波动值如附图4所示，可见，最大速度波动为0.1746％，远大于本方法的速度波动值。
通过将本发明基于二阶Runge-Kutta法和参数补偿的插补参数计算方法与传统二阶泰勒级数展开法对比可见，本发明插补参数计算方法的加工进给速度波动值远小于传统方法，可以在不进行迭代计算的情况下有效降低速度波动，提高曲线插补加工质量。
本发明针对样条曲线插补过程中，由于曲线参数与弧长具有非线性对应关系所导致的进给速度波动问题，发明了基于二阶Runge-Kutta法和参数补偿的一种样条曲线插补参数计算方法，用来在不进行迭代计算的前提下有效抑制曲线插补过程中的速度波动，对提高样条曲线插补加工质量具有重大实用意义。

Claims (1)

1.一种样条曲线插补参数计算方法，其特性在于，该方法首先将各插补周期的指令参数值写成关于时间的函数，并利用二阶Runge-Kutta法计算下一插补周期的指令参数初始值；其次，通过比较利用参数初始值计算得到的进给速度与理想进给速度，计算参数补偿值；最后，根据参数初始值和参数补偿值，确定下一插补点曲线参数值，实现速度波动小的样条曲线实时插补；方法具体步骤如下：
第一步 计算插补参数初始值
将插补参数u写成时间t的函数：
u＝u(t) (1)
则参数u对时间t的导数为：
<mrow> <mfrac> <mrow> <mi>d</mi> <mi>u</mi> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <msub> <mo>|</mo> <mrow> <mi>t</mi> <mo>=</mo> <msub> <mi>t</mi> <mi>i</mi> </msub> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <mi>d</mi> <mi>u</mi> </mrow> <mrow> <mi>d</mi> <mi>s</mi> </mrow> </mfrac> <msub> <mo>|</mo> <mrow> <mi>u</mi> <mo>=</mo> <msub> <mi>u</mi> <mi>i</mi> </msub> </mrow> </msub> <mo>&amp;CenterDot;</mo> <mfrac> <mrow> <mi>d</mi> <mi>s</mi> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <msub> <mo>|</mo> <mrow> <mi>t</mi> <mo>=</mo> <msub> <mi>t</mi> <mi>i</mi> </msub> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <mfrac> <mrow> <mi>d</mi> <mi>s</mi> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <msub> <mo>|</mo> <mrow> <mi>t</mi> <mo>=</mo> <msub> <mi>t</mi> <mi>i</mi> </msub> </mrow> </msub> </mrow> <mrow> <mfrac> <mrow> <mi>d</mi> <mi>s</mi> </mrow> <mrow> <mi>d</mi> <mi>u</mi> </mrow> </mfrac> <msub> <mo>|</mo> <mrow> <mi>u</mi> <mo>=</mo> <msub> <mi>u</mi> <mi>i</mi> </msub> </mrow> </msub> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
其中，s为曲线弧长，ti为第i个插补周期时的总插补时间，ui为第i个插补周期的指令参数值；令第i个插补周期的理想进给速度为vi，则：
<mrow> <mfrac> <mrow> <mi>d</mi> <mi>s</mi> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <msub> <mo>|</mo> <mrow> <mi>t</mi> <mo>=</mo> <msub> <mi>t</mi> <mi>i</mi> </msub> </mrow> </msub> <mo>=</mo> <msub> <mi>v</mi> <mi>i</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
设样条曲线参数方程为C＝C(u)，则：
<mrow> <mfrac> <mrow> <mi>d</mi> <mi>s</mi> </mrow> <mrow> <mi>d</mi> <mi>u</mi> </mrow> </mfrac> <msub> <mo>|</mo> <mrow> <mi>u</mi> <mo>=</mo> <msub> <mi>u</mi> <mi>i</mi> </msub> </mrow> </msub> <mo>=</mo> <mo>|</mo> <mo>|</mo> <msup> <mi>C</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>|</mo> <mo>|</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>
其中，C′(ui)为曲线在参数u＝ui位置处的一阶导失；将公式(3)、(4)代入公式(2)得：
<mrow> <mfrac> <mrow> <mi>d</mi> <mi>u</mi> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <msub> <mo>|</mo> <mrow> <mi>t</mi> <mo>=</mo> <msub> <mi>t</mi> <mi>i</mi> </msub> </mrow> </msub> <mo>=</mo> <mfrac> <msub> <mi>v</mi> <mi>i</mi> </msub> <mrow> <mo>|</mo> <mo>|</mo> <msup> <mi>C</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>|</mo> <mo>|</mo> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>
因此，利用二阶Runge-Kutta法展开公式(1)得到下一个，即第(i+1)个插补周期的指令插补参数初始值us,i+1为：
<mrow> <msub> <mi>u</mi> <mrow> <mi>s</mi> <mo>,</mo> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>=</mo> <msub> <mi>u</mi> <mi>i</mi> </msub> <mo>+</mo> <mfrac> <mrow> <mi>T</mi> <mi>s</mi> </mrow> <mn>2</mn> </mfrac> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>k</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>
其中，Ts为插补周期，
第二步 计算插补参数补偿值
为保证实际插补进给速度和理想进给速度相等，第(i+1)个插补周期的指令参数值ui+1应满足：
<mrow> <mfrac> <mrow> <mo>|</mo> <mo>|</mo> <mi>C</mi> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mi>C</mi> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>|</mo> <mo>|</mo> </mrow> <mrow> <mi>T</mi> <mi>s</mi> </mrow> </mfrac> <mo>=</mo> <msub> <mi>v</mi> <mi>i</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>
该插补点的插补参数补偿值Δui+1应满足：
Δui+1＝ui+1-us,i+1 (8)
根据一阶泰勒级数展开，可得到：
C(ui+1)＝C(us,i+1)+C′(us,i+1)Δui+1 (9)
将公式(9)代入公式(7)得：
<mrow> <mo>|</mo> <mo>|</mo> <mi>C</mi> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mrow> <mi>s</mi> <mo>,</mo> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msup> <mi>C</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mrow> <mi>s</mi> <mo>,</mo> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> <msub> <mi>&amp;Delta;u</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <mi>C</mi> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>|</mo> <msup> <mo>|</mo> <mn>2</mn> </msup> <mo>=</mo> <msubsup> <mi>v</mi> <mi>i</mi> <mn>2</mn> </msubsup> <msup> <mi>Ts</mi> <mn>2</mn> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow>
整理得：
<mrow> <msubsup> <mi>A&amp;Delta;u</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> <mn>2</mn> </msubsup> <mo>+</mo> <msub> <mi>B&amp;Delta;u</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>+</mo> <mi>D</mi> <mo>=</mo> <mn>0</mn> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo> </mrow> </mrow>
其中，A＝||C′(us,i+1)||2，B＝2(C′(us,i+1)×(C(us,i+1)-C(ui)))，
<mrow> <mi>D</mi> <mo>=</mo> <mo>|</mo> <mo>|</mo> <mi>C</mi> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mrow> <mi>s</mi> <mo>,</mo> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mi>C</mi> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>|</mo> <msup> <mo>|</mo> <mn>2</mn> </msup> <mo>-</mo> <msubsup> <mi>v</mi> <mi>i</mi> <mn>2</mn> </msubsup> <msup> <mi>Ts</mi> <mn>2</mn> </msup> <mo>;</mo> </mrow>
解方程(11)得插补参数的两个根Δui+1,1和Δui+1,2分别为：
<mrow> <msub> <mi>&amp;Delta;u</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> <mo>,</mo> <mn>1</mn> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <mo>-</mo> <mi>B</mi> <mo>+</mo> <msqrt> <mrow> <msup> <mi>B</mi> <mn>2</mn> </msup> <mo>-</mo> <mn>4</mn> <mi>A</mi> <mi>D</mi> </mrow> </msqrt> </mrow> <mrow> <mn>2</mn> <mi>A</mi> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>12</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mi>&amp;Delta;u</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <mo>-</mo> <mi>B</mi> <mo>-</mo> <msqrt> <mrow> <msup> <mi>B</mi> <mn>2</mn> </msup> <mo>-</mo> <mn>4</mn> <mi>A</mi> <mi>D</mi> </mrow> </msqrt> </mrow> <mrow> <mn>2</mn> <mi>A</mi> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>13</mn> <mo>)</mo> </mrow> </mrow>
由于因此Δui+1,1≈0，Δui+1,2≈-B/A；因此，为保证插补过程的稳定性，当方程(11)有实数根，即判别式B2-4AD≥0时，取绝对值较小的Δui+1,1作为插补参数补偿值Δui+1；当方程(11)无实数根，即判别式B2-4AD<0时，令参数补偿值Δui+1为零；因此，参数补偿值Δui+1的计算方法为：
<mrow> <msub> <mi>&amp;Delta;u</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mfrac> <mrow> <mo>-</mo> <mi>B</mi> <mo>+</mo> <msqrt> <mrow> <msup> <mi>B</mi> <mn>2</mn> </msup> <mo>-</mo> <mn>4</mn> <mi>A</mi> <mi>D</mi> </mrow> </msqrt> </mrow> <mrow> <mn>2</mn> <mi>A</mi> </mrow> </mfrac> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <msup> <mi>B</mi> <mn>2</mn> </msup> <mo>-</mo> <mn>4</mn> <mi>A</mi> <mi>D</mi> <mo>&amp;GreaterEqual;</mo> <mn>0</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mn>0</mn> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <msup> <mi>B</mi> <mn>2</mn> </msup> <mo>-</mo> <mn>4</mn> <mi>A</mi> <mi>D</mi> <mo>&lt;</mo> <mn>0</mn> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>14</mn> <mo>)</mo> </mrow> </mrow>
第三步根据插补参数初始值us,i+1和插补参数补偿值Δui+1确定第(i+1)个插补周期的插补参数值ui+1：
ui+1＝us,i+1+Δui+1 (15)
利用该参数值作为数控系统插补过程中第i个插补周期的输入指令进行实际插补加工；判断是否到达曲线终点，若到达终点，则结束插补，否则，令i＝i+1，返回第一步。
CN201610165047.2A 2016-03-22 2016-03-22 一种样条曲线插补参数计算方法 Active CN105843174B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610165047.2A CN105843174B (zh) 2016-03-22 2016-03-22 一种样条曲线插补参数计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610165047.2A CN105843174B (zh) 2016-03-22 2016-03-22 一种样条曲线插补参数计算方法

Publications (2)

Publication Number Publication Date
CN105843174A CN105843174A (zh) 2016-08-10
CN105843174B true CN105843174B (zh) 2018-05-08

Family

ID=56588427

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610165047.2A Active CN105843174B (zh) 2016-03-22 2016-03-22 一种样条曲线插补参数计算方法

Country Status (1)

Country Link
CN (1) CN105843174B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108594765A (zh) * 2018-04-08 2018-09-28 哈尔滨理工大学 一种逐点比较法平面直线插补计算方法
CN109164759B (zh) * 2018-08-01 2021-12-21 深圳市汇川技术股份有限公司 曲线插补方法、设备及计算机可读存储介质
CN109283892B (zh) * 2018-11-08 2020-11-17 北京航空航天大学 一种基于参数化曲线几何特征和弓高误差限制的进给速率自适应插补算法
CN109597962B (zh) * 2018-12-07 2021-04-30 山东大学 一种基于预测和迭代补偿的参数插值方法
CN111983975B (zh) * 2020-07-23 2023-04-18 江苏威泽智能科技股份有限公司 一种参数增量补偿的样条曲线插补方法
CN114035512B (zh) * 2021-11-22 2022-09-23 广东工业大学 一种基于两次泰勒展开的最小速度波动插补方法
CN114221825B (zh) * 2022-02-23 2022-05-24 武汉海牛泰克信息技术有限公司 一种基于数字信息传输的加密系统
CN114740803B (zh) * 2022-04-06 2022-12-27 南华大学 一种nurbs曲线插补速度控制方法、装置及设备

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010104916A1 (en) * 2009-03-10 2010-09-16 Complex Data Technologies Hardware and method of numerical integration
CN103101053A (zh) * 2011-11-07 2013-05-15 精工爱普生株式会社 机器人控制系统、机器人系统以及传感器信息处理装置
CN103809521A (zh) * 2012-11-14 2014-05-21 中国科学院沈阳计算技术研究所有限公司 基于弦截法的样条曲线插补方法
CN103926881A (zh) * 2014-04-21 2014-07-16 北京航空航天大学 一种基于割线法的无速度波动参数曲线直接插补方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010104916A1 (en) * 2009-03-10 2010-09-16 Complex Data Technologies Hardware and method of numerical integration
CN103101053A (zh) * 2011-11-07 2013-05-15 精工爱普生株式会社 机器人控制系统、机器人系统以及传感器信息处理装置
CN103809521A (zh) * 2012-11-14 2014-05-21 中国科学院沈阳计算技术研究所有限公司 基于弦截法的样条曲线插补方法
CN103926881A (zh) * 2014-04-21 2014-07-16 北京航空航天大学 一种基于割线法的无速度波动参数曲线直接插补方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
CNC加工中直线插补刀路的误差评估修正研究;张永超等;《组合机床与自动化加工技术》;20110630(第6期);第61-64页 *
开放数控参数化样条曲线插补指令开发;高雄;《中国优秀硕士学位论文全文数据库 工程科技I辑》;20110531(第5期);全文 *
高性能复杂曲面零件测量-再设计-数字加工一体_贾振元化加工方法;贾振元等;《机械工程学报》;20131031;第49卷(第19期);第126-132页 *

Also Published As

Publication number Publication date
CN105843174A (zh) 2016-08-10

Similar Documents

Publication Publication Date Title
CN105843174B (zh) 一种样条曲线插补参数计算方法
Jia et al. A NURBS interpolator with constant speed at feedrate-sensitive regions under drive and contour-error constraints
Du et al. An accurate adaptive NURBS curve interpolator with real-time flexible acceleration/deceleration control
Lee et al. The feedrate scheduling of NURBS interpolator for CNC machine tools
Tsai et al. A real-time predictor-corrector interpolator for CNC machining
CN102591257B (zh) 面向参数曲线刀具轨迹的数控系统轮廓误差控制方法
Chen et al. Augmented Taylor's expansion method for B-spline curve interpolation for CNC machine tools
Hou et al. Geometric modelling of thin-walled blade based on compensation method of machining error and design intent
CN104597845A (zh) 一种用于高质量加工的样条曲线插补算法
CN104898564A (zh) 一种降低三轴联动轮廓误差的方法
CN107479497A (zh) 一种五轴加工轨迹轮廓误差双闭环补偿方法
Dong et al. An accurate NURBS curve interpolation algorithm with short spline interpolation capacity
Song et al. Estimation and compensation for continuous-path running trajectory error in high-feed-speed machining
CN104317251A (zh) 基于Obrechkoff算法的三次NURBS曲线实时插补方法
Du et al. Sliding mode control with third-order contour error estimation for free-form contour following
Wang et al. Research on adaptive non-uniform rational B-spline real-time interpolation technology based on acceleration constraints
Jia et al. Adaptive estimation and nonlinear variable gain compensation of the contouring error for precise parametric curve following
CN114740803A (zh) 一种nurbs曲线插补速度控制方法、装置及设备
Zhao et al. Contour Error Coupled-Control Strategy based on Line Interpolation and Curve Interpolation.
Simba et al. Iterative contouring controller design for biaxial feed drive systems
Song et al. A precise direct parametric interpolation method for NURBS-toolpath generation
Jiang et al. A real-time feedrate planning method and efficient interpolator with minimal feedrate fluctuation for parametric toolpath
Liu et al. A cross-coupled iterative learning control design for biaxial systems based on natural local approximation of contour error
Bhattacharjee et al. Development of a CNC interpolation scheme for CNC controller based on Runge-Kutta method
Sun et al. Smoothing interpolation of five-axis tool path with less feedrate fluctuation and higher computation efficiency

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