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 PDFInfo
- 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
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/404—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 control arrangements for compensation, e.g. for backlash, overshoot, tool offset, tool wear, temperature, machine construction errors, load, inertia
-
- 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/35—Nc in input of data, input till input file format
- G05B2219/35408—Calculate 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
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.
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)
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)
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 |
-
2016
- 2016-12-30 CN CN201611258996.1A patent/CN106774154B/en not_active Expired - Fee Related
Patent Citations (4)
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)
Title |
---|
杨振玲,等: "等弓高误差变步长的插补点求取方法优化", 《工具技术》 * |
Cited By (7)
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 |