CN106774154A - A kind of space curve interpolating method theoretical based on osculating plane - Google Patents

A kind of space curve interpolating method theoretical based on osculating plane Download PDF

Info

Publication number
CN106774154A
CN106774154A CN201611258996.1A CN201611258996A CN106774154A CN 106774154 A CN106774154 A CN 106774154A CN 201611258996 A CN201611258996 A CN 201611258996A CN 106774154 A CN106774154 A CN 106774154A
Authority
CN
China
Prior art keywords
point
space curve
bow
error
equation
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.)
Granted
Application number
CN201611258996.1A
Other languages
Chinese (zh)
Other versions
CN106774154B (en
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.)
Taiyuan University of Science and Technology
Original Assignee
Taiyuan University of Science and 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 Taiyuan University of Science and Technology filed Critical Taiyuan University of Science and Technology
Priority to CN201611258996.1A priority Critical patent/CN106774154B/en
Publication of CN106774154A publication Critical patent/CN106774154A/en
Application granted granted Critical
Publication of CN106774154B publication Critical patent/CN106774154B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B19/00Programme-control systems
    • G05B19/02Programme-control systems electric
    • G05B19/18Numerical 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/404Numerical 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 control arrangements for compensation, e.g. for backlash, overshoot, tool offset, tool wear, temperature, machine construction errors, load, inertia
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00Program-control systems
    • G05B2219/30Nc systems
    • G05B2219/35Nc in input of data, input till input file format
    • G05B2219/35408Calculate new position data from actual data to compensate for contour error

Landscapes

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

Abstract

The invention belongs to Computerized Numerical Control processing technology field, it is proposed that a kind of space curve interpolating method theoretical based on osculating plane, it is mainly used in tool-path planning operation in Multi-axis Machining.According to the local geometric shape at interpolated point on space curve, using the interpolated point of the bow error high such as equidistant surface construction satisfaction on the principal normal direction from tangent plane of each interpolated point.But the bow error high that this side up bow error high not necessarily truly, bow error function high is set up followed by analytic geometry midpoint to straight line apart from implicationThe most longbow error high calculated between two adjacent discrete points is verified.Using new method energy self-adaptative adjustment cutting step, bow error high truly is realized, be evenly distributed in the bow error high respectively approached in straightway, machining accuracy is high, and surface quality is good.

Description

A kind of space curve interpolating method theoretical based on osculating plane
Technical field:
The invention belongs to Computerized Numerical Control processing technology field, more particularly to a kind of space curve interpolation side theoretical based on osculating plane Method, is mainly used in tool-path planning operation in Machining of Curved Surface.
Background technology:
In the process industries such as aviation, aerospace industry, automobile and mould, existing largely has complex-curved part, These parts not only complex geometry, and the accuracy of manufacture and surface quality requirements are very high.With computer technology and certainly The development of dynamic control technology is more and more extensive with application of the complex-curved part in industrial products.At present, it is complex-curved The processing of part is mostly realized by Digit Control Machine Tool.Tool-path planning is the important composition of design for NC Part, its designing quality directly affects the machining accuracy and processing efficiency of part.The cutter path of Machining of Curved Surface is by one Bar cutter-contact point trace curve is constituted, and from geometrically saying, each cutter-contact point trace curve is essentially all a space curve, because This tool-path planning is converted to the interpolated point computational problem of a rule space curve.Different interpolation point calculating methods are with regard to shape Into different interpolating methods.
At present, range of application wider interpolation point calculating method such as mainly has at parameter step length method, equidistant step length, the step-length Screening method etc..
It is that parameter and curve equidistant are split Deng parameter step length method and equidistant step length, although both approaches are simple, But required in order that obtaining every section of approximate error of straight line and meeting, step-length or parameter increase can only be selected using adverse conditions the most Take, so as to cause remaining processing sequences huge, actual bow error skewness high.
Step-length screening method uses less equidistance or parameter discrete cutter rail before this using based on equal steps Mark, then checks every section of actual approximate error of straight line, finally needs to delete unnecessary discrete point, so that remaining trajectory error It is more uniform.Although step-length screening method the shortcoming of parameter step length method and equidistant step length such as overcomes, but check every section it is straight All it is approximate bow error high during bow error high in line, so as to influence the surface quality of curved surface, reduces the processing essence of workpiece Degree.
The content of the invention:
In order to solve the problems, such as conventional method in terms of Multi-axis Machining tool-path planning step size computation, the present invention is carried A kind of space curve interpolating method of the satisfaction bow high franchise theoretical based on osculating plane is gone out.It can be according to P on space curve0Point The local geometric shape at place, is adaptively adjusted cutting step, and actual bow in every section of straight line is calculated using error function high is bent Error high so that be evenly distributed in the bow error high respectively approached in straightway, machining accuracy is high, and surface quality is good.
According to the relevant knowledge of Differential Geometry, at any on space curve, three orthogonal tangent lines, principal normals With the fundamental three-vector form that binormal constitutes curve, the osculating plane of this point is generated by them, normal plane and from tangent plane. In these three planes, osculating plane is the plane of " most pressing close to curve ", is tangent line and perpendicular to osculating plane from tangent plane.It is based on This conclusion, obtains discrete interpolated point, due at each point in the bow such as satisfaction error requirements high from the principal normal direction of tangent plane Bent curvature of a curve is different with torsion, error now bow error high not necessarily truly, so also needing to utilize The most longbow error high bent between error function calculating two-phase neighbor interpolation point high is verified.
According to being described above, the space curve interpolating method based on the theoretical satisfaction bow franchise high of osculating plane includes following step Suddenly:
Step (1) assumes parametric equation r=r (t) of space curve Γ={ x (t), y (t), z (t) } t ∈ [a, b], When interpolation starts, the 0th interpolated point P0Parameter be t0=a, bends franchise ε high, the bow franchise size influence free curve high Surface quality and machining accuracy, it is better to bend the smaller surface quality of franchise high, and the bow franchise high is by manufacturer according to being processed Surface set, be known quantity;
Step (2) is in P0At point, r (t are calculated0),r′(t0),r″(t0), try to achieve equation the r " (t from tangent plane π0)·[ρ-r (t0)]=0, ρ be any point from tangent plane vector, to bend franchise ε high be distance construction from tangent plane π on principal normal direction Equidistant surface π1
Step (3) simultaneous equidistant surface π1The equation of equation and space curve Γ, solves equidistant surface π1With the friendship of space curve Γ The parameter t of point1
Step (4) calculates t0And t1Two corresponding interpolated point P0And P1Between most longbow error max δ high.Solve most longbow Error high has mid-point method, dichotomy and Fibonacci method, but is not most longbow error high truly, error degree ratio Higher, the present invention establishes P by using the knowledge of the distance at analytic geometry midpoint to straight line0And P1Any point between 2 points Bow error function highThe method being most worth is asked to try to achieve max δ using function;
Step (5) is if max δ<ε, then P1As next interpolated point, satisfaction bow Gao Yun is otherwise found using line search method The point that difference is required is used as next interpolated point;
Step (6) return to step (2);Until an end of job for space curve.
Brief description of the drawings
Fig. 1 is equidistant surface schematic diagram.
Fig. 2 is flow chart of the method for the present invention.
Fig. 3 is the interpolated point schematic diagram of space curve Γ.
Fig. 4 is the most longbow error curve diagram high between space curve Γ interpolated points.
Specific embodiment
In order to make the purpose , technical scheme and advantage of the present invention be clearer, below in conjunction with accompanying drawing 1-4 and embodiment Son, is further described in detail to the present invention.Specific implementation example described herein is only used for the explanation present invention, not For limiting the present invention.
1. assume that a parametric equation of space curve Γ is r=r (t)={ t, t2,t3T ∈ [0,100], the 0th interpolation Point P0Parameter be t0=0, bend franchise ε=0.001 high.
2. in P0At point, r (t are calculated0),r′(t0),r″(t0),
Then r (t0)={ t0,t0 2,t0 3, r, ' (t0)={ 1,2t0,3t0 2, r " (t0)={ 0,2,6t0,
Obtain and be from tangent plane π equations:r″(t0)·[ρ-r(t0)]=0, i.e. y+3t0z-t0 2-3t0 4=0, ρ are to take up an official post from tangent plane The vector of a bit.
To bend franchise ε high it is being distance construction from the equidistant surface π of tangent plane π from the principal normal direction of tangent plane π1
3. simultaneous equidistant surface π1The equation of equation and space curve Γ, obtains intersection point P1Parameter t1It is equationRoot, using matlab softwares solve.
4. P is calculated0And P1Most longbow error max δ high between two-phase neighbor interpolation.Using distance between beeline and dot in analytic geometry Implication, establish P0And P1The bow error function high of any point between 2 pointsProfit The method being most worth is asked to try to achieve max δ with function;
If 5. max δ<ε, then P1As next interpolated point, interval is otherwise progressively reduced using line search method and is found completely The point of arch of foot tolerance requirements high is used as next interpolated point;
6. return to step (2);Until all interpolated points on a space curve have determined.
Space curve Γ interpolated points schematic diagram and most longbow error curve diagram high are as shown in figs. 34.

Claims (2)

1. a kind of space curve interpolating method theoretical based on osculating plane, it is characterised in that comprise the following steps that:
Step (1) assumes that a parametric equation of space curve Γ is r=r (t) t ∈ [a, b], when interpolation starts, the 0th interpolation Point P0Parameter be t0=a, bends franchise ε high;
Step (2) is in P0At point, r (t are calculated0),r′(t0),r″(t0), try to achieve equation the r " (t from tangent plane π0)·[ρ-r(t0)] =0, ρ are the vector of any point from tangent plane, to bend franchise ε high are distance construction from the equidistant of tangent plane π on principal normal direction Face π1
Step (3) simultaneous equidistant surface π1The equation of equation and space curve Γ, solves equidistant surface π1With the intersection point of space curve Γ Parameter t1
Step (4) calculates t0And t1Two corresponding interpolated point P0And P1Between most longbow error max δ high;
Step (5) is if max δ<ε, then P1As next interpolated point, interval is otherwise progressively reduced using line search method and is found completely The point of arch of foot tolerance requirements high is used as next interpolated point;
Step (6) return to step (2);Until an end of job for space curve.
2. the space curve interpolating method theoretical based on osculating plane according to claim 1, it is characterised in that
Step (4) sets up P0And P1The bow error function high of any point between 2 points The method being most worth is asked to try to achieve max δ using function.
CN201611258996.1A 2016-12-30 2016-12-30 A kind of space curve interpolating method based on osculating plane theory Expired - Fee Related CN106774154B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611258996.1A CN106774154B (en) 2016-12-30 2016-12-30 A kind of space curve interpolating method based on osculating plane theory

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611258996.1A CN106774154B (en) 2016-12-30 2016-12-30 A kind of space curve interpolating method based on osculating plane theory

Publications (2)

Publication Number Publication Date
CN106774154A true CN106774154A (en) 2017-05-31
CN106774154B CN106774154B (en) 2018-11-13

Family

ID=58953769

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611258996.1A Expired - Fee Related CN106774154B (en) 2016-12-30 2016-12-30 A kind of space curve interpolating method based on osculating plane theory

Country Status (1)

Country Link
CN (1) CN106774154B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107145126A (en) * 2017-06-29 2017-09-08 南京航空航天大学 Consider the numerical control machining knife rail subregion mapping method of error band distribution
CN110116353A (en) * 2019-06-03 2019-08-13 武汉理工大学 A kind of blade front and rear edge robot abrasive band grinding and polishing step-length optimization method
CN111983975A (en) * 2020-07-23 2020-11-24 江苏大学 Spline curve interpolation method for parameter increment compensation
CN112987647A (en) * 2021-01-04 2021-06-18 苏州科技大学 Equal-error tool path generation method for five-axis machining of parametric curved surface flat-bottom tool
CN114563979A (en) * 2022-02-25 2022-05-31 苏州谷夫道自动化科技有限公司 Compensation method and system for single-channel semi-closed-loop multi-spindle multi-station machining

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6225771B1 (en) * 1999-12-01 2001-05-01 General Electric Company Probe chord error compensation
CN105045211A (en) * 2015-07-13 2015-11-11 太原科技大学 Equal-chord-error variable-step tangent interpolation method
CN105676787A (en) * 2015-12-28 2016-06-15 龙兵 Elliptic arc interpolation algorithm
CN105759725A (en) * 2016-03-22 2016-07-13 大连理工大学 Speed-sensitive section constant-speed curve interpolation speed planning method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6225771B1 (en) * 1999-12-01 2001-05-01 General Electric Company Probe chord error compensation
CN105045211A (en) * 2015-07-13 2015-11-11 太原科技大学 Equal-chord-error variable-step tangent interpolation method
CN105676787A (en) * 2015-12-28 2016-06-15 龙兵 Elliptic arc interpolation algorithm
CN105759725A (en) * 2016-03-22 2016-07-13 大连理工大学 Speed-sensitive section constant-speed curve interpolation speed planning method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
杨振玲,等: "等弓高误差变步长的插补点求取方法优化", 《工具技术》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107145126A (en) * 2017-06-29 2017-09-08 南京航空航天大学 Consider the numerical control machining knife rail subregion mapping method of error band distribution
CN110116353A (en) * 2019-06-03 2019-08-13 武汉理工大学 A kind of blade front and rear edge robot abrasive band grinding and polishing step-length optimization method
CN111983975A (en) * 2020-07-23 2020-11-24 江苏大学 Spline curve interpolation method for parameter increment compensation
CN112987647A (en) * 2021-01-04 2021-06-18 苏州科技大学 Equal-error tool path generation method for five-axis machining of parametric curved surface flat-bottom tool
CN112987647B (en) * 2021-01-04 2022-03-29 苏州科技大学 Equal-error tool path generation method for five-axis machining of parametric curved surface flat-bottom tool
CN114563979A (en) * 2022-02-25 2022-05-31 苏州谷夫道自动化科技有限公司 Compensation method and system for single-channel semi-closed-loop multi-spindle multi-station machining
CN114563979B (en) * 2022-02-25 2023-10-13 苏州谷夫道自动化科技有限公司 Compensation method and system for single-channel semi-closed-loop multi-spindle multi-station processing

Also Published As

Publication number Publication date
CN106774154B (en) 2018-11-13

Similar Documents

Publication Publication Date Title
CN106774154B (en) A kind of space curve interpolating method based on osculating plane theory
CN102637216B (en) Method for generating numerical-control side milling machining tool path for complicated curved surfaces
Feng et al. Constant scallop-height tool path generation for three-axis sculptured surface machining
CN102289534A (en) Method for modeling involute helical gear accurately
CN108062073B (en) Circular arc smooth compression interpolation method for high-quality machining
CN102581705B (en) Division compensation method for machining error during milling of workpiece
CN103412978A (en) Stretch-bending process model correction method based on workpiece three-dimension scanning
CN108073138B (en) Elliptical arc smooth compression interpolation algorithm suitable for high-speed high-precision machining
CN110069041A (en) A kind of Work piece processing method and system based on on-machine measurement
CN109597359B (en) Numerical controller
CN106774153A (en) A kind of profile tolerance compensation method of S-shaped rose cutter
CN110688709A (en) Workpiece point cloud model-based skin process model correction method
CN105045211A (en) Equal-chord-error variable-step tangent interpolation method
CN103092134A (en) Spiral interpolation numerical control (NC) machining tool path generation method of three dimensional trimming deflashing cutting edges
CN104317251A (en) Three-order NURBS curve real-time interpolation method based on Obrechkoff algorithm
CN107037779A (en) Free form surface NC process tool track optimizing methods under non-homogeneous tolerance
CN107322008A (en) The three-dimensional modeling method of machined surface in a kind of turning technique
CN104035379B (en) Compensation method for perpendicularity error between feed shafts
CN106292531B (en) Algorithm for calculating profile boundary of ZN1 worm disc-shaped forming cutter
CN112051803B (en) Small line segment fairing method based on space plane normal vector
Lo et al. An improved method for scheduling the tool paths for three-axis surface machining
CN105116840A (en) Compensation method for slowly varying error of primary interpolation of space curve
Du et al. The avoidance of cutter gouging in five-axis machining with a fillet-end milling cutter
Li et al. Analysis and optimization of the properties of the grinding contact line for form-grinding modified helical gears
CN103543691A (en) Rotary shaft position optimizing method applicable to high-speed high-precision processing of free-form curved surface

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: 20181113

Termination date: 20201230