CN109669395B - Variable-radius circular interpolation method for axisymmetric aspheric surface - Google Patents

Variable-radius circular interpolation method for axisymmetric aspheric surface Download PDF

Info

Publication number
CN109669395B
CN109669395B CN201811471921.0A CN201811471921A CN109669395B CN 109669395 B CN109669395 B CN 109669395B CN 201811471921 A CN201811471921 A CN 201811471921A CN 109669395 B CN109669395 B CN 109669395B
Authority
CN
China
Prior art keywords
interpolation
radius
aspheric surface
interval
axisymmetric aspheric
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
CN201811471921.0A
Other languages
Chinese (zh)
Other versions
CN109669395A (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.)
Tianjin Jinhang Institute of Technical Physics
Original Assignee
Tianjin Jinhang Institute of Technical Physics
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 Tianjin Jinhang Institute of Technical Physics filed Critical Tianjin Jinhang Institute of Technical Physics
Priority to CN201811471921.0A priority Critical patent/CN109669395B/en
Publication of CN109669395A publication Critical patent/CN109669395A/en
Application granted granted Critical
Publication of CN109669395B publication Critical patent/CN109669395B/en
Active 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
    • 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)
  • Numerical Control (AREA)

Abstract

The invention discloses a variable radius circular interpolation method of an axisymmetric aspheric surface, which comprises the following steps: a. determining an interpolated diameter D; b. determining an interpolation interval dx; c. calculating an interpolation radius ri of any interpolation interval; d. and d, calculating the interpolation radius value of each interpolation interval according to the established method in the steps a, b and c, and finishing the interpolation. On the basis of the existing interpolation mode, the invention improves the mode of adopting variable-radius circular interpolation, and effectively reduces the influence of interpolation errors on the surface shape precision of the axisymmetric aspheric surface; the method has simple calculation, can directly carry out parameterized optimization after programming through programming software, does not need to supplement hardware, can improve the processing precision, has high efficiency and low cost, and is a quick and effective axisymmetric aspheric surface programming interpolation method.

Description

Variable-radius circular interpolation method for axisymmetric aspheric surface
Technical Field
The invention belongs to the technical field of ultra-precision machining, and relates to a variable-radius circular interpolation method for an axisymmetric aspheric surface.
Background
With the continuous improvement of the requirements of optical systems, the application of the axisymmetric aspheric surface is more and more extensive, wherein ultra-precision milling and single-point diamond turning are one of two main processing methods for processing the axisymmetric aspheric surface. The common point of the two is that the cutter and the processed surface are in point-to-point contact mode in the processing process, and the processing precision is influenced by various factors such as cutter profile precision, machine tool precision, tool setting precision, clamping stress, numerical control interpolation precision and the like in the contact mode. Among the above influencing factors, numerical control interpolation accuracy is the only soft influencing factor, i.e., one of the ways of improving the machining accuracy by optimizing the interpolation program.
At present, the interpolation mode commonly used in the process of the axisymmetric aspheric surface interpolation is mainly a linear interpolation mode, and in commercial software, two interpolation point selection modes are provided, one is an equal linear quantity interval interpolation point selection method, and the other is an equal arc length interval selection method. However, in any selection method, in the linear interpolation method, a straight line connection method is adopted between two adjacent points, and a rise difference, that is, a surface shape accuracy error, is generated between the connection straight line and the aspheric surface curve between the two points in the optical axis direction. Therefore, it is one of the directions of research to change the interpolation method between two points after the interpolation point is determined to reduce the surface shape accuracy error.
Disclosure of Invention
Objects of the invention
The purpose of the invention is: the variable-radius circular interpolation method for the axisymmetric aspheric surface adopts a variable-radius circular interpolation mode, so that the adjustable-radius circular interpolation is adopted between two adjacent interpolation points, the vector height difference is reduced, and the fitting precision of the axisymmetric aspheric surface is improved.
(II) technical scheme
In order to solve the above technical problem, the present invention provides a variable radius circular interpolation method for an axisymmetric aspheric surface, which includes the following steps:
a. determining an interpolated diameter D;
b. determining an interpolation interval dx;
c. calculating an interpolation radius ri of any interpolation interval;
d. and d, calculating the interpolation radius value of each interpolation interval according to the established method in the steps a, b and c, and finishing the interpolation.
(III) advantageous effects
The variable-radius circular interpolation method for the axisymmetric aspheric surface provided by the technical scheme improves the mode of adopting variable-radius circular interpolation on the basis of the conventional interpolation mode, and effectively reduces the influence of interpolation errors on the surface shape precision of the axisymmetric aspheric surface; the method has simple calculation, can directly carry out parameterized optimization after programming through programming software, does not need to supplement hardware, can improve the processing precision, has high efficiency and low cost, and is a quick and effective axisymmetric aspheric surface programming interpolation method.
Drawings
Fig. 1 is a schematic diagram of an axisymmetric aspheric surface and a circular interpolation.
FIG. 2 is a schematic diagram of the process of the present invention.
Fig. 3 is a view of the arc interpolation according to the present invention.
Fig. 4 is a schematic view of an aspheric surface.
Detailed Description
In order to make the objects, contents, and advantages of the present invention clearer, the following detailed description of the embodiments of the present invention will be made in conjunction with the accompanying drawings and examples.
The invention relates to a variable radius circular interpolation method of an axisymmetric aspheric surface, which comprises the following steps:
a. determining an interpolated diameter D
And determining the interpolation diameter D according to the axial symmetry aspheric surface equation z ═ f (x), drawing requirement indexes, the processed material, the blank and other factors.
The interpolation diameter D should be larger than the design diameter D0 of the axisymmetric aspheric surface that needs to be processed actually.
Typically, D is 0.5mm to 2mm greater than D0.
b. The interpolation interval dx is determined.
And determining an interpolation interval dx according to parameters such as an aspheric equation z ═ f (x), drawing requirement indexes, an interpolation diameter D and the like.
Further, the interpolation interval dx is determined to be equal in the x-axis direction, i.e., dx is constant.
Further, D/(2dx) should be an integer, i.e., the patch point must pass through the aspheric vertex.
Further, dx is selected as: when interpolation is carried out according to the straight line, the maximum value of the difference value of the straight line and the axial symmetry aspheric surface along the Z direction does not exceed the PV value required in the technical index, and the PV value is the difference between the maximum value and the minimum value in a comparison curve of the actual surface and the theoretical surface.
c. An interpolation radius ri of an arbitrary interpolation interval is calculated.
And (3) regarding the aspheric surface between any intervals as an off-axis aspheric surface, and determining the best fitting spherical surface according to the principle of minimum height difference along the z-axis. The sphere at this time is a circle passing through the two interpolation points and having a center on the z-axis, and the radius of the circle is the interpolation radius riThe specific calculation method is as follows:
suppose that two interpolation points at an arbitrary interpolation interval are Pi(xi，zi) And Pi+1(xi+1，zi+1) The midpoint of the line between the two points is Pm(xm，zm) Then, then
xi+1-xi＝dx (2)
According to the attached figure 1, let
Figure BDA0001891169490000031
Figure BDA0001891169490000032
Then the process of the first step is carried out,
zm＝k(xm-xi)+zi (5)
suppose OiHas a coordinate value of (0, Z)oi) Then, then
Figure BDA0001891169490000033
From the trigonometric relationship, r can be obtainediIs expressed as
Figure BDA0001891169490000041
Further, the joint type (3) to (7) can obtain the value of the arc interpolation curvature radius required for the interval.
d. And d, calculating the interpolation radius value of each interpolation interval according to the established method in the steps a, b and c, and finishing the interpolation.
Examples
With reference to the schematic interpolation diagrams of aspheric surfaces shown in fig. 1 and fig. 2 and the flowchart of the present invention, the following steps are described for determining the specific implementation of the method of the present invention, taking the aspheric surface shown in fig. 4 as an example, with reference to other drawings:
a. an interpolated radius D is determined. In the aspheric drawing shown in fig. 4, D0 requires 11mm, parameters of the aspheric surface are given, and the interpolation radius D is 12mm according to the requirements of D0.
b. The interpolation interval dx is determined. According to the determination principle that passes through the origin of the aspheric surface and that D/(2dx) should be an integer, in this example, dx is determined as:
dx＝1mm
c. an interpolation radius of an arbitrary interpolation interval is calculated. Because of the axisymmetric aspheric surface, the aspheric surface is completely consistent in the areas of-6-0 and 0-6 in the x-axis direction. Therefore, only the directions of 0 to 6 are considered as xi2 and xi+1Two points corresponding to 3 are calculated as an interval.
X is to bei2 and xi+1Substituting the aspheric equations with the parameters shown in fig. 4 into 3, the coordinate between the two adjacent points can be obtained as Pi(2, 0.2338) and Pi+1(3，0.537)
The coordinates of the two points obtained are substituted for equations (3) to (7), and the radius of the interpolated arc at that time can be determined to be ri＝8.6306mm。
d. According to the above calculation method, when dx is 1, the corresponding ri value is:
r1 r2 r3 r4 r5 r6
8.6859mm 8.6689mm 8.6306mm 8.5596mm 8.4436mm 8.3mm
by simulation, it can be seen that the interpolation curve is shown in fig. 3.
Since the present invention is directed to an axisymmetric aspherical surface, only the x-axis positive portion may be considered in the calculation.
The above description is only a preferred embodiment of the present invention, and it should be noted that, for those skilled in the art, several modifications and variations can be made without departing from the technical principle of the present invention, and these modifications and variations should also be regarded as the protection scope of the present invention.

Claims (1)

1. A variable radius circular interpolation method of an axisymmetric aspheric surface is characterized by comprising the following steps:
a. determining an interpolated diameter D;
b. determining an interpolation interval dx;
c. calculating an interpolation radius r of an arbitrary interpolation intervali；
d. According to the method established in the steps a, b and c, calculating the interpolation radius value of each interpolation interval to finish the interpolation;
in the step a, the interpolation diameter D is larger than the design diameter D0 of the axisymmetric aspheric surface which needs to be processed actually;
in the step a, D is 0.5mm-2mm larger than D0;
in the step b, the interpolation interval dx is determined as equal interval along the x axis, that is, dx is a constant;
in the step b, D/(2dx) is an integer;
in the step c, the aspheric surface between any intervals is regarded as an off-axis aspheric surface, and the best fitting spherical surface is determined according to the principle of minimum height difference along the z-axis; the sphere at this time is a circle passing through the two interpolation points and having a center on the z axis, and the radius of the circle is the interpolation radius ri;
the interpolated radius riThe calculation method comprises the following steps:
let P be the two interpolation points at any interpolation intervali(xi，zi) And Pi+1(xi+1，zi+1) The midpoint of the line between the two points is Pm(xm，zm) Then, then
xi+1-xi＝dx (2)
Order:
Figure FDA0003064469450000011
Figure FDA0003064469450000012
then:
zm＝k(xm-xi)+zi (5)
suppose OiHas a coordinate value of (0, Z)oi) And then:
Figure FDA0003064469450000021
from the trigonometric relationship, r is obtainediIs expressed as
Figure FDA0003064469450000022
Further, the joint type (3) to (7) can obtain the value of the arc interpolation curvature radius required for the interval.
CN201811471921.0A 2018-12-04 2018-12-04 Variable-radius circular interpolation method for axisymmetric aspheric surface Active CN109669395B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811471921.0A CN109669395B (en) 2018-12-04 2018-12-04 Variable-radius circular interpolation method for axisymmetric aspheric surface

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811471921.0A CN109669395B (en) 2018-12-04 2018-12-04 Variable-radius circular interpolation method for axisymmetric aspheric surface

Publications (2)

Publication Number Publication Date
CN109669395A CN109669395A (en) 2019-04-23
CN109669395B true CN109669395B (en) 2021-08-06

Family

ID=66144693

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811471921.0A Active CN109669395B (en) 2018-12-04 2018-12-04 Variable-radius circular interpolation method for axisymmetric aspheric surface

Country Status (1)

Country Link
CN (1) CN109669395B (en)

Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0348762A (en) * 1989-07-17 1991-03-01 Olympus Optical Co Ltd Manufacture of acoustic lens
CN1738697A (en) * 2003-10-29 2006-02-22 精工爱普生株式会社 Aspherical surface processing method, aspherical surface forming method and aspherical surface processing apparatus
CN1928755A (en) * 2006-09-27 2007-03-14 沈阳飞机工业（集团）有限公司 Spherical and cylindrical surface numeral control processing method
CN101125411A (en) * 2007-07-30 2008-02-20 厦门大学 Parallel grinding and cutting method for non-axial-symmetry and non-ball-surface optical element
CN101428399A (en) * 2008-12-20 2009-05-13 厦门大学 Abrasive machining method for cuneiform aspheric surface
CN101508025A (en) * 2009-03-13 2009-08-19 厦门大学 Processing control method of axial symmetry free-form surface of aspheric surface optical elements
CN102087517A (en) * 2010-07-19 2011-06-08 长春理工大学 Method for reducing speed interpolation error and hardware system
CN102608956A (en) * 2012-03-05 2012-07-25 浙江工业大学 NURBS (Non-Uniform Rational B-Spline) curve adaptive interpolation control method based on de Boor algorithm
CN202922332U (en) * 2012-07-24 2013-05-08 浙江陀曼精密机械有限公司 Numerical control composite shaft non-circular cross section turning processing control system
CN103862380A (en) * 2014-03-26 2014-06-18 上海交通大学 Analyzing and compensating system for optical mirror face few-axle grinding inclining cambered face grinding wheel errors
CN104668643A (en) * 2014-12-26 2015-06-03 潍柴动力股份有限公司 Method for machining oil guide groove of diesel engine
CN105867311A (en) * 2016-04-18 2016-08-17 哈尔滨工业大学 Method for precisely turning high-gradient aspheric surface through arc-edge diamond cutter
CN107272597A (en) * 2017-07-14 2017-10-20 福建工程学院 A kind of nurbs curve interpolation based on advance and retreat method quickly pre-reads processing method

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3578634B2 (en) * 1998-07-17 2004-10-20 東芝機械株式会社 Creating an arc command
JP4537134B2 (en) * 2004-07-12 2010-09-01 セイコーオプティカルプロダクツ株式会社 Aspheric eyeglass lens design method
CN105965349B (en) * 2016-06-29 2018-12-28 成都精密光学工程研究中心 Large-caliber off-axis non-spherical mirror accurate grinding system

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0348762A (en) * 1989-07-17 1991-03-01 Olympus Optical Co Ltd Manufacture of acoustic lens
CN1738697A (en) * 2003-10-29 2006-02-22 精工爱普生株式会社 Aspherical surface processing method, aspherical surface forming method and aspherical surface processing apparatus
CN1928755A (en) * 2006-09-27 2007-03-14 沈阳飞机工业（集团）有限公司 Spherical and cylindrical surface numeral control processing method
CN101125411A (en) * 2007-07-30 2008-02-20 厦门大学 Parallel grinding and cutting method for non-axial-symmetry and non-ball-surface optical element
CN101428399A (en) * 2008-12-20 2009-05-13 厦门大学 Abrasive machining method for cuneiform aspheric surface
CN101508025A (en) * 2009-03-13 2009-08-19 厦门大学 Processing control method of axial symmetry free-form surface of aspheric surface optical elements
CN102087517A (en) * 2010-07-19 2011-06-08 长春理工大学 Method for reducing speed interpolation error and hardware system
CN102608956A (en) * 2012-03-05 2012-07-25 浙江工业大学 NURBS (Non-Uniform Rational B-Spline) curve adaptive interpolation control method based on de Boor algorithm
CN202922332U (en) * 2012-07-24 2013-05-08 浙江陀曼精密机械有限公司 Numerical control composite shaft non-circular cross section turning processing control system
CN103862380A (en) * 2014-03-26 2014-06-18 上海交通大学 Analyzing and compensating system for optical mirror face few-axle grinding inclining cambered face grinding wheel errors
CN104668643A (en) * 2014-12-26 2015-06-03 潍柴动力股份有限公司 Method for machining oil guide groove of diesel engine
CN105867311A (en) * 2016-04-18 2016-08-17 哈尔滨工业大学 Method for precisely turning high-gradient aspheric surface through arc-edge diamond cutter
CN107272597A (en) * 2017-07-14 2017-10-20 福建工程学院 A kind of nurbs curve interpolation based on advance and retreat method quickly pre-reads processing method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Numerical Control Interpretation Algorithm of Aspheric Surface Based on the Genetic Algorithms and Neural Network;F. H. Zhang;《Advances in Optics Design and Precision Manufacturering Technologies pt.1》;20070113;第25-29页 *
基于轴对称非球面子午线的步长不变式双圆弧插补算法;张学忱，等;《机械工程学报》;20130531;第49卷(第9期);第144-150页 *

Also Published As

Publication number Publication date
CN109669395A (en) 2019-04-23

Similar Documents

Publication Publication Date Title
Takeuchi et al. Generation of 5-axis control collision-free tool path and postprocessing for NC data
CN111347294B (en) High-gradient optical mirror surface error polishing correction processing method
CN102078968B (en) Ultraprecise turning method of optical lens with incline
CN109918807B (en) Local tool path smoothing method for optimizing curvature
CN106774154B (en) A kind of space curve interpolating method based on osculating plane theory
CN103056731A (en) Five-axis precision ultrasonic milling machining method of large-aperture off-axis aspheric mirror
CN103862380B (en) The few axle grinding of optical mirror plane tilts cambered surface grinding wheel error and analyzes and make-up system
CN102319921B (en) Hierarchical machining method of tiltable main shaft numerical control milling machine
CN108747602B (en) Method for grinding high-order curve section cylindrical surface array structure by using circular arc array grinding wheel
CN109213083A (en) A kind of ruled surface processing path generation method, device and equipment
CN106774153A (en) A kind of profile tolerance compensation method of S-shaped rose cutter
CN113319656B (en) Grinding method for special-shaped workpiece
CN101839149A (en) Computing method of revolution section profile of turbomachinery blade
CN109669395B (en) Variable-radius circular interpolation method for axisymmetric aspheric surface
Weinert et al. Milling simulation for process optimization in the field of die and mould manufacturing
US20170343982A1 (en) Method for machining a workpiece by means of a chip-removing tool on a numerically-controlled machine tool
CN112526925B (en) Profile finish machining method based on three-dimensional cam profile materialized model deviation compensation
CN117666475A (en) Continuous short line segment corner processing path smoothing method
CN113111404A (en) Fitting method for circular arc and straight line of small continuous spatial line segment
CN102350522B (en) Processing method of numerical control milling machine with tiltable main shaft
CN109491321B (en) Contour error estimation method based on H-shaped precision motion platform
CN104385084A (en) Five-axis grinding method for variably-formed base circle planar envelope convex surface workpiece
CN111026035B (en) Method for solving cyclone milling blade tool location point based on curvature change
CN109669396B (en) Variable interval linear interpolation method for axisymmetric aspheric surface
CN114818071A (en) Railway line plane curve calculating and drawing method

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