[go: up one dir, main page]

CN114425732B - Automatic optimization method, system and medium for sub-caliber processing technology - Google Patents

Automatic optimization method, system and medium for sub-caliber processing technology Download PDF

Info

Publication number
CN114425732B
CN114425732B CN202210355185.2A CN202210355185A CN114425732B CN 114425732 B CN114425732 B CN 114425732B CN 202210355185 A CN202210355185 A CN 202210355185A CN 114425732 B CN114425732 B CN 114425732B
Authority
CN
China
Prior art keywords
spectral density
density function
volume
processing technology
function
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
CN202210355185.2A
Other languages
Chinese (zh)
Other versions
CN114425732A (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.)
Laser Fusion Research Center China Academy of Engineering Physics
Original Assignee
Laser Fusion Research Center China Academy of Engineering 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 Laser Fusion Research Center China Academy of Engineering Physics filed Critical Laser Fusion Research Center China Academy of Engineering Physics
Priority to CN202210355185.2A priority Critical patent/CN114425732B/en
Publication of CN114425732A publication Critical patent/CN114425732A/en
Application granted granted Critical
Publication of CN114425732B publication Critical patent/CN114425732B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B24GRINDING; POLISHING
    • B24BMACHINES, DEVICES, OR PROCESSES FOR GRINDING OR POLISHING; DRESSING OR CONDITIONING OF ABRADING SURFACES; FEEDING OF GRINDING, POLISHING, OR LAPPING AGENTS
    • B24B13/00Machines or devices designed for grinding or polishing optical surfaces on lenses or surfaces of similar shape on other work; Accessories therefor
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Operations Research (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Mechanical Engineering (AREA)
  • Length Measuring Devices By Optical Means (AREA)
  • Constituent Portions Of Griding Lathes, Driving, Sensing And Control (AREA)

Abstract

The invention is suitable for the technical field of optical processing, and provides an automatic optimization method, a system and a medium for a sub-caliber processing technology, wherein the automatic optimization method for the sub-caliber processing technology comprises the following steps: obtaining effective removal rate spectrums of different removal functions, wherein the effective removal rate spectrums are convergence rates of correction of error volumes of all spatial frequencies by the removal functions; acquiring a volume spectral density function of the optical element, wherein the volume spectral density function is the density of the volume of residual error material contained in the surface shape error of the optical element at each frequency; and obtaining the optimal processing technology through the effective removal rate spectrum and the volume spectral density function. The automatic optimization method, system and medium for the sub-caliber processing technology provided by the invention have the advantages of high processing efficiency and low production cost.

Description

Automatic optimization method, system and medium for sub-caliber processing technology
Technical Field
The invention relates to the technical field of optical processing, in particular to an automatic optimization method, system and medium for a sub-caliber processing technology.
Background
The sub-aperture polishing technology is proposed by U.S. R.A. Jones in 60 s, and the basic principle is that a tool far smaller than the aperture of a processing element is adopted, the residence time is controlled by a numerical control technology to remove local materials, the error convergence is realized by combining high-precision surface shape detection iterative processing, and the surface shape error RMS of the large-aperture optical element, which is better than 0.01 μm, can be obtained. Due to the advantages of high precision and high certainty, the sub-aperture polishing technology is widely applied to the final precise shaping stage of ultra-precision machining, and is the most feasible means for realizing high-precision machining of large-aperture optical elements at present. Advanced optical manufacturing technologies such as small tool numerical control, air bags, atmospheric plasma, magnetorheological, jet flow, ion beams and the like all belong to sub-aperture polishing technologies, the processing precision and efficiency of optical elements are greatly improved, but the ultra-precision processing cost of the existing large-aperture optical elements is still very high, and tens of thousands and millions of times are used.
In the sub-aperture polishing digital automatic control model, the surface shape error represents the distribution of the spatial frequency band error of the function to be removed of the element to be processed, the removal function represents the capability of a processing tool for correcting the spatial frequency band error, and the removal function can represent the removal rate distribution of the tool in a fixed-point processing influence range by using a two-dimensional matrix.
Although the sub-aperture polishing technologies such as magnetorheological, ion beam and small tool numerical control are adopted, the processing of the large-aperture optical element with the RMS value less than 0.1 μm can be realized, and as the surface shape error distribution has the characteristics of randomness and diversity of removal functions, in order to obtain the final precision, the low removal rate and small-size removal function which are conservative as much as possible are generally adopted, so that the processing working hours are greatly increased.
The method for representing the correction capability of the removal function mainly comprises two methods, namely, the convergence capability of the removal function on the error distribution of a processed element is indirectly reflected according to a processing experiment result, a large number of experiment results are required to be used as a premise, and an experimenter can have a certain degree of real operability only by having a certain experience accumulation; and secondly, the Fourier transform of the removal function is directly carried out, and the correction capability of removing the spatial frequency band error can be qualitatively represented. The two methods can only evaluate the correction capability of the removal function through experience or qualitative characterization among a plurality of removal functions, and the targeted process optimization cannot be made according to different conditions and different surface shape error distributions of the optical element.
In summary, none of the automatic optimization methods for the sub-aperture processing techniques in the prior art solves the problem that the most suitable process parameters (removal functions) for the current element cannot be obtained by combining the surface shape error distribution of the current element in actual processing, and the comprehensive efficiency of sub-aperture polishing is seriously reduced.
Disclosure of Invention
The invention aims to provide an automatic optimization method, system and medium for a sub-caliber machining process, which have high machining efficiency and low production cost.
The invention provides an automatic optimization method of a sub-caliber processing technology in a first aspect, which comprises the following steps:
step S10: obtaining effective removal rate spectrums of different removal functions, wherein the effective removal rate spectrums are convergence rates of correction of error volumes of all spatial frequencies by the removal functions;
step S20: acquiring a volume spectral density function of the optical element, wherein the volume spectral density function is the density of the volume of residual error material contained in the surface shape error of the optical element at each frequency;
step S30: and obtaining the optimal processing technology through the effective removal rate spectrum and the volume spectral density function.
Further, step S10 further includes: and obtaining an effective removal rate spectrum according to the ratio of the volume variation of the surface shape error to the processing time.
Further, a preferred processing technique is obtained by the effective removal rate spectrum and the volume spectral density function.
Further, step S20 further includes calculating a time spectral density function according to the effective removal rate spectrum and the volume spectral density function, and integrating the time spectral density function to obtain the processing time of the removal function; the removal function is a preferred machining process when the machining time of the removal function is minimal.
Further, in step S20, the temporal spectral density function is a volume spectral density function divided by the effective removal rate spectrum.
Further, step S20 further includes: and integrating the volume spectral density function at a preset space frequency band to obtain the surface shape error volume of the optical element under the preset space frequency band, calculating the surface shape error volume outside the cut-off frequency of the removal function, and selecting the removal function with the minimum volume difference value as the optimal processing technology.
Further, the volume spectral density function is a two-dimensional volume spectral density function or a one-dimensional volume spectral density function.
The invention provides a preferred system of a sub-caliber processing technology, which comprises an effective removal rate spectrum acquisition module, a volume spectrum density function acquisition module and a preferred processing technology acquisition module; the effective removal rate spectrum acquisition module is used for acquiring effective removal rate spectrums of different removal functions; the volume spectral density function acquisition module is used for acquiring a volume spectral density function of the optical element; the optimal processing technology obtaining module is used for obtaining an optimal processing technology through the effective removal rate spectrum and the volume spectrum density function.
A third aspect of the invention provides a readable storage medium for storing a program which, when executed, implements the automated preferred method of sub-aperture machining process.
In summary, the present invention has at least the following technical effects:
1. according to the invention, the quantitative characterization and the frequency-band division characterization of the surface shape error volume of the optical element on a frequency domain are realized through the volume spectrum density function of the optical element, the surface shape error volume quantity of the optical element on different frequency bands can be obtained, and accurate data reference is provided for the processing of the optical element;
2. according to the invention, the optimized processing technology is obtained by effectively removing the rate spectrum and the volume spectrum density function, so that the volume data of the removed function and the surface shape error volume data of the optical element can be matched with each other, thereby effectively solving the problem that the removed function in the prior art cannot be combined with the surface shape error distribution of the optical element to obtain the most suitable technological parameters of the optical element, effectively improving the comprehensive efficiency of sub-aperture polishing, reducing the production cost and bringing huge economic benefits;
3. according to the method, the effective removal rate spectrum is obtained through the ratio of the surface shape error volume variation to the processing time, the concept of the effective removal rate spectrum directly related to the form, the size and the like of the removal function is provided, the capability of the removal function for correcting the spatial frequency error can be quantitatively represented by utilizing the effective removal rate spectrum, and the quantitative representation of the removal function is realized;
4. according to the method, the optimized processing technology is obtained through the effective removal rate spectrum and the volume spectrum density function, the effective removal rate of the removal function is combined with the surface shape error distribution to be removed of the optical element, and the combination of the removal functions can be used for improving the processing rate or the processing precision of the optical element, so that experimenters can perform optimization according to actual conditions;
5. according to the invention, the processing time of the removal function is obtained by integrating the time spectrum density function, so that the removal function with the minimum processing time is selected as an optimal processing technology, the processing speed of the optical element is effectively increased, the comprehensive efficiency of sub-aperture polishing is improved, and the method has great economic benefits and practical production application value;
6. according to the invention, the volume spectral density function is integrated in a preset spatial frequency band to obtain the surface shape error volume, the volume difference value between the residual error material volume and the surface shape error volume is calculated, and then the removal function with the minimum volume difference value is selected as the optimal processing technology.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings needed to be used in the embodiments of the present invention or in the description of the prior art will be briefly described below, and it is obvious that the drawings described below are only some embodiments of the present invention, and it is obvious for those skilled in the art that other drawings can be obtained according to the drawings without creative efforts.
FIG. 1 is a schematic diagram of an automatically preferred method and system for the sub-caliber manufacturing process of the present invention;
FIG. 2 is a plot of the areal error 1 volume spectral density function of the present invention;
FIG. 3 is a graph of the effective removal rate spectrum for removal function 1 in accordance with the present invention;
FIG. 4 is a first schematic diagram of an automatically preferred method of the sub-caliber manufacturing process of the present invention;
FIG. 5 is a graph of time spectral density in the present invention;
FIG. 6 is a graph of the error distribution and values of profile 1 in the present invention;
FIG. 7 is a graph of the error distribution and values for profile 2 of the present invention;
FIG. 8 is a graph of the morphological distribution and parameters of the removal function 1 of the present invention;
FIG. 9 is a graph of the morphological distribution and parameters of the removal function 2 of the present invention;
FIG. 10 is a second schematic diagram of an automatically preferred method of the sub-caliber manufacturing process of the present invention;
fig. 11 is a schematic diagram of the rand rotation transformation function of the present invention.
Detailed Description
The following description provides many different embodiments, or examples, for implementing different features of the invention. The particular examples set forth below are intended as a brief description of the invention and are not intended as limiting the scope of the invention.
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. The components of embodiments of the present invention generally described and illustrated in the figures herein may be arranged and designed in a wide variety of different configurations. Thus, the following detailed description of the embodiments of the present invention, presented in the figures, is not intended to limit the scope of the invention, as claimed, but is merely representative of selected embodiments of the invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments of the present invention without making any creative effort, shall fall within the protection scope of the present invention.
It should be noted that: like reference numbers and letters refer to like items in the following figures, and thus, once an item is defined in one figure, it need not be further defined and explained in subsequent figures. Meanwhile, in the description of the present invention, the terms "first", "second", and the like are used only for distinguishing the description, and are not construed as indicating or implying relative importance.
The first embodiment is as follows:
as shown in fig. 1, a first embodiment of the present invention provides an automatic optimization method for a sub-caliber processing technology, including the following steps:
step S10: obtaining effective removal rate spectrums of different removal functions, wherein the effective removal rate spectrums are convergence rates of correction of error volumes of all spatial frequencies by the removal functions;
step S20: acquiring a volume spectral density function of the optical element, wherein the volume spectral density function is the density of the volume of residual error material contained in the surface shape error of the optical element at each frequency;
step S30: and obtaining the optimal processing technology through the effective removal rate spectrum and the volume spectral density function.
The embodiment of the automatic optimization method for the neutron caliber processing technology enables quantitative representation of matching of the removal function and the surface shape error, solves the problem that in the prior art, in the machining process of the neutron caliber, the requirement of the surface shape error of an element to be processed needs to be met, reasonable technological parameters can be selected for machining, automatic optimization of the technological parameters can be achieved, the removal function which is most matched with the spatial distribution of the surface shape error of the current element to be processed is obtained, high-efficiency machining is achieved, and production cost is reduced.
Volume spectral density functionVSD(ω)The quantitative characterization mathematical model of the error volume content of the surface shape error in each spatial frequency band can be characterized, and the frequency domain quantitative characterization of the error volume content of the surface shape is obtained, as shown in fig. 2.
The quantitative characterization and the frequency-division characterization of the surface shape error volume of the optical element on the frequency domain are realized through the volume spectral density function of the optical element, the quantity of the surface shape error volumes of the optical element on different frequency bands can be obtained, and accurate data reference is provided for the processing of the optical element.
The optimal processing technology is obtained by effectively removing the rate spectrum and the volume spectrum density function, so that the volume data of the removed function and the surface shape error volume data of the optical element can be matched with each other, the matching degree of the removed function and the surface shape error of the element to be processed is realized, and an effective and reliable technical approach is provided for automatic intelligent optimal selection of the ultra-precision processing technology parameters of the optical element, thereby effectively solving the problem that the removed function in the prior art cannot be combined with the surface shape error distribution of the optical element to obtain the most suitable technological parameters of the optical element, effectively improving the comprehensive efficiency of sub-aperture polishing, reducing the production cost and bringing huge economic benefits;
further, step S10 further includes: and obtaining an effective removal rate spectrum through the ratio of the volume change of the surface shape error to the processing time.
The effective removal rate spectrum is calculated based on airspace simulation to obtain an effective removal rate spectrum of the error of each space frequency band trimmed by the removal function, and quantitative representation of the error correction capability of each space frequency band by the removal function is obtained.
The surface shape error volume variation is the volume variation of the surface shape error before and after processing under the preset spatial frequency, and the effective removal rate spectrum is obtained according to the ratio of the surface shape error volume variation to the processing timeRE(ω)As shown in fig. 3, the quantitative characterization of the removal function of all sub-aperture polishing such as small tool numerical control polishing, air bag polishing, jet polishing, atmospheric plasma polishing, magnetorheological polishing, ion beam polishing, etc. can be realized, the concept of the effective removal rate spectrum directly related to the form, size, etc. of the removal function is provided, the capability of the removal function to correct the spatial frequency error can be quantitatively characterized by using the effective removal rate spectrum, and the quantitative characterization of the removal function is realized.
Further, a preferred processing technique is obtained by the effective removal rate spectrum and the volume spectral density function.
The optimal processing technology is obtained through the effective removal rate spectrum and the volume spectrum density function, the effective removal rate of the removal function is combined with the surface shape error distribution to be removed of the optical element, the combination of the removal functions can be used for improving the processing rate or the processing precision of the optical element, and experimenters can perform optimization according to actual conditions conveniently.
Further, as shown in fig. 4, step S20 further includes calculating a time spectral density function according to the effective removal rate spectrum and the volume spectral density function, and integrating the time spectral density function to obtain a processing time of the removal function; the removal function is a preferred machining process when the machining time of the removal function is minimal
By efficient removal rate spectroscopyRE(ω)And characterizing the correction capability of the removal function, and quantitatively characterizing the effective removal rate spectrum of the removal function for correcting the error of each frequency band. By volume spectral density functionVSD(ω)Quantitatively characterizing the volume content of each frequency band.
After the correction capability of the removal function and the representation of the surface shape error on the frequency domain are respectively obtained, the information is combined to quantitatively represent the matching degree of the removal function and the surface shape error, namely a time spectral density functionT(ω). And the time spectrum density function is used for representing the removal function and the sub-band matching representation of the surface shape error.
As shown in fig. 5, the function of the volume spectral density can be knownVSD(ω)And effective removal Rate SpectrumRE(ω)Curve obtaining time density spectrum functionT(ω)The time spectrum density curve is integrated, so that the theoretical total time for removing the function correction surface shape error can be obtained, and the theoretical time for removing the function correction surface shape error of each spatial frequency band of the element to be processed, namely the processing time, is obtained.
For the error distribution of the element to be processed, the error distribution of all the removing functions to be selected can be calculated quicklyT i (ω)Correcting the processing time corresponding to the target frequency band error by comparisont i iIn order to remove the serial number of the function, a simple sorting is performed, that is, the removal function with the highest processing efficiency can be selected preferably, so as to improve the processing efficiency, as shown in fig. 4.
In order to verify whether the obtained time spectrum density is in accordance with the real machining process, verification is performed by comparing whether the time of the time spectrum density is in accordance with the actual simulation machining. Orthogonal tests are performed through two groups of surface shape errors (surface shape 1 shown in figure 6 and surface shape 2 shown in figure 7) and two groups of removal functions (removal function 1 shown in figure 8 and removal function 2 shown in figure 9), four groups of time spectrum densities are obtained through calculation and simulation processing, detailed information and feature description of specific parameters are shown in table 1, and removal functions with different sizes and surface shape error data with obvious difference in error frequency domain distribution are selected as cross validation parameters.
Table 1 simulation process input parameter conditions
Figure 346600DEST_PATH_IMAGE001
The cross validation results are shown in table 2, which indicates that the total processing time obtained by calculation by using the time spectrum density curve is basically consistent with the time of simulation processing, the deviation is less than 1%, and the correctness of the model is fully proved.
TABLE 2 simulation verification results
Figure 607948DEST_PATH_IMAGE002
In the first embodiment, the processing time of the removal function is obtained by integrating the time spectrum density function, so that the removal function with the minimum processing time is selected as the optimal processing technology, the processing speed of the optical element is effectively increased, the comprehensive efficiency of sub-aperture polishing is improved, and the method has great economic benefits and practical production and application values.
Further, in step S20, the temporal spectral density function is a volume spectral density function divided by the effective removal rate spectrum.
Volume spectral density functionVSD(ω)And effective removal rate spectrumRE(ω)The division can obtain the time spectrum density functionT(ω)The physical meaning of which is to be modified by a removal functionThe surface shape error corresponds to the time consumed by the unit frequency volume on the frequency band, and the obtained full-band curve is the time spectral density curve.
Further, as shown in fig. 10, step S20 further includes: and integrating the volume spectral density function at a preset space frequency band to obtain the surface shape error volume of the optical element under the preset space frequency band, calculating the surface shape error volume outside the cut-off frequency of the removal function, and selecting the removal function with the minimum volume difference value as the optimal processing technology.
And integrating the volume spectral density function at a preset spatial frequency band to obtain a surface shape error volume, calculating the volume difference value between the volume of the residual error material and the surface shape error volume, and selecting the removal function with the minimum volume difference value as the optimal processing technology. The volume spectral density function is integrated in a preset spatial frequency band to obtain a surface shape error volume, the volume difference value between the volume of the residual error material and the surface shape error volume is calculated, and then the removal function with the minimum volume difference value is selected as an optimal processing technology.
It should be noted that, when a person skilled in the art selects the removal function, the general rule is: the smaller the size of the machining tool, the better the machining accuracy. However, this rule also has a certain limitation, and in practical application, the machining precision of the large-size machining tool is better than that of the small-size machining tool, because the fundamental reason for determining the machining precision is the distribution form of the removal function, not the size of the machining tool. Therefore, the embodiment can effectively avoid errors caused by the fact that a person skilled in the art only selects the size of the machining tool, and improves the polishing efficiency.
Further, the volume spectral density function is a two-dimensional volume spectral density function or a one-dimensional volume spectral density function.
The two-dimensional volume spectral density function
Figure 715581DEST_PATH_IMAGE003
Is calculated by the formula
Figure 314053DEST_PATH_IMAGE004
Wherein, in the step (A),Sthe area of the element is defined as the area of the element,ΔLin order to be the sampling interval of the sample,xin a first direction, the first direction is,yin the second direction, the first direction is the first direction,Nis the number of sampling points in the first direction,Mis the number of sample points in the second direction,
Figure 971168DEST_PATH_IMAGE005
is a fourier transform of the arithmetic square root of the surface error,
Figure 121527DEST_PATH_IMAGE006
is the frequency domain spatial coordinate of the first direction,
Figure 857401DEST_PATH_IMAGE007
is the frequency domain spatial coordinate of the second direction.
Fourier transform of arithmetic square root of surface shape error
Figure 400509DEST_PATH_IMAGE005
The calculation formula of (c) is:
Figure 272650DEST_PATH_IMAGE008
wherein, in the step (A),
Figure 859490DEST_PATH_IMAGE009
is a two-dimensional matrix of surface-shaped error amplitude values, i is an imaginary number unit,jis the serial number of the sampling point in the second direction,kis the number of the first direction sampling points,j、kis a positive integer.
The surface shape error amplitude two-dimensional matrixz(j,k) The optical element can be detected by a detection instrument such as a three-coordinate measuring machine and an interferometer.
The one-dimensional volume spectral density function
Figure 193912DEST_PATH_IMAGE010
Is calculated by the formula
Figure 399766DEST_PATH_IMAGE011
Wherein, in the process,
Figure 251047DEST_PATH_IMAGE012
for the rotational transformation of a two-dimensional volumetric spectral density function,θin order to be the angle of rotation,
Figure 149733DEST_PATH_IMAGE013
is the sample length of the first direction after the rotation transformation.
The rotational transformation two-dimensional volume spectral density function
Figure 735566DEST_PATH_IMAGE014
Is calculated by the formula
Figure 338586DEST_PATH_IMAGE015
Wherein, in the step (A),randis a rotational transformation function.
randThe rotation transformation is as shown in FIG. 11 whenθ=0, obtainXThe calculation formula of the one-dimensional volume spectral density of the direction is as follows:
Figure 450898DEST_PATH_IMAGE016
wherein, in the step (A),L y a second direction sample length in cartesian coordinates. When in useθ=(ii) 90 °, obtainYThe calculation formula of the one-dimensional volume spectral density of the direction is as follows:
Figure 628808DEST_PATH_IMAGE017
wherein, in the step (A),L x the first direction sample length in cartesian coordinates.
The one-dimensional volume spectral density function
Figure 560992DEST_PATH_IMAGE018
A certain frequency band of
Figure 233281DEST_PATH_IMAGE019
Is integrated as the error volume in that frequency band
Figure 200100DEST_PATH_IMAGE020
I.e. by
Figure 50376DEST_PATH_IMAGE021
Wherein, in the step (A),
Figure 594490DEST_PATH_IMAGE022
the lower form error amplitude of the spatial frequency band is preset.
The two-dimensional volume spectral density function
Figure 945837DEST_PATH_IMAGE023
Or one-dimensional volume spectral density function
Figure 878414DEST_PATH_IMAGE024
The predetermined spatial frequency band integral of (a) is equal to the surface shape error volume of the optical element, i.e.
Figure 289804DEST_PATH_IMAGE025
Wherein, in the step (A),
Figure 321214DEST_PATH_IMAGE026
is the average height of the surface-shaped errors,Vis the surface error volume.
Example two:
as shown in fig. 1, a second embodiment of the present invention provides a preferred system for a sub-aperture processing technology, including an effective removal rate spectrum obtaining module, a volume spectrum density function obtaining module, and a preferred processing technology obtaining module; the effective removal rate spectrum acquisition module is used for acquiring effective removal rate spectrums of different removal functions; the volume spectral density function acquisition module is used for acquiring a volume spectral density function of the optical element; the optimal processing technology obtaining module is used for obtaining an optimal processing technology through the effective removal rate spectrum and the volume spectrum density function.
Example three:
the third embodiment of the invention provides a readable storage medium for storing a program, and the program is used for realizing the automatic optimization method of the sub-aperture machining process when being executed.
Example four:
an embodiment of the present invention provides an electronic device, including: one or more processors; a memory having one or more programs stored thereon; the one or more programs, when executed by the one or more processors, cause the one or more processors to implement the automated preferred method of sub-aperture machining.
The automatic optimization method, the system and the medium for the sub-caliber processing technology can realize the quantitative characterization of the surface shape error and the removal function matching degree in the sub-caliber polishing processing process, can be suitable for processing any surface shape error distribution by all sub-caliber polishing technologies, and have the advantages of simple and rapid realization process and good stability. Meanwhile, by quantitative representation of the matching degree of the removal function and the surface shape error, the removal function with the optimal matching degree of the surface shape error of the element to be processed can be automatically optimized, the uncertain influence of human factors on the production process is avoided, the processing efficiency of the large-caliber optical element is greatly improved, the production cost is reduced, and great economic benefit can be brought.
The above description is intended to be illustrative of the preferred embodiment of the present invention and should not be taken as limiting the invention, but rather, the intention is to cover all modifications, equivalents, and alternatives falling within the spirit and scope of the invention.

Claims (8)

1.一种子口径加工工艺的自动优选方法,其特征在于,包括如下步骤:1. an automatic optimal method of sub-caliber processing technology, is characterized in that, comprises the steps: 步骤S10:获取不同去除函数的有效去除速率谱,所述有效去除速率谱为去除函数修正各空间频率误差体积的收敛速率;Step S10: obtaining effective removal rate spectra of different removal functions, where the effective removal rate spectrum is the convergence rate at which the removal function corrects each spatial frequency error volume; 步骤S20:获取光学元件的体积谱密度函数,所述体积谱密度函数为光学元件的面形误差在各频率下所含残余误差材料体积的密度;Step S20: obtaining a volume spectral density function of the optical element, where the volume spectral density function is the density of the residual error material volume contained in the surface shape error of the optical element at each frequency; 步骤S30:通过所述有效去除速率谱和体积谱密度函数得到优选加工工艺;Step S30: obtaining a preferred processing technology through the effective removal rate spectrum and the volume spectral density function; 步骤S10还包括:通过面形误差体积变化量与加工时间的比值,得到有效去除速率谱;Step S10 further includes: obtaining an effective removal rate spectrum by the ratio of the surface shape error volume change and the processing time; 所述体积谱密度函数为二维体积谱密度函数或一维体积谱密度函数;The volume spectral density function is a two-dimensional volume spectral density function or a one-dimensional volume spectral density function; 所述二维体积谱密度函数
Figure 257700DEST_PATH_IMAGE001
的计算公式为
Figure 670226DEST_PATH_IMAGE002
,其中,S为元件面积,ΔL为采样间隔,x为第一方向,y为第二方向,N为第一方向的采样点数,M为第二方向的采样点数,
Figure 963804DEST_PATH_IMAGE003
为面形误差算术平方根的傅里叶变换,
Figure 410966DEST_PATH_IMAGE004
为第一方向的频域空间坐标,
Figure 866218DEST_PATH_IMAGE005
为第二方向的频域空间坐标;
The two-dimensional volume spectral density function
Figure 257700DEST_PATH_IMAGE001
The calculation formula is
Figure 670226DEST_PATH_IMAGE002
, where S is the element area, ΔL is the sampling interval, x is the first direction, y is the second direction, N is the number of sampling points in the first direction, M is the number of sampling points in the second direction,
Figure 963804DEST_PATH_IMAGE003
is the Fourier transform of the arithmetic square root of the surface error,
Figure 410966DEST_PATH_IMAGE004
is the frequency domain spatial coordinate of the first direction,
Figure 866218DEST_PATH_IMAGE005
is the frequency domain spatial coordinate in the second direction;
所述一维体积谱密度函数
Figure 766041DEST_PATH_IMAGE006
的计算公式为
Figure 597731DEST_PATH_IMAGE007
,其中,
Figure 915711DEST_PATH_IMAGE008
为旋转变换二维体积谱密度函数,θ为旋转角度,
Figure 541865DEST_PATH_IMAGE009
为旋转变换后第一方向的采样长度。
The one-dimensional volume spectral density function
Figure 766041DEST_PATH_IMAGE006
The calculation formula is
Figure 597731DEST_PATH_IMAGE007
,in,
Figure 915711DEST_PATH_IMAGE008
is the rotation transform two-dimensional volume spectral density function, θ is the rotation angle,
Figure 541865DEST_PATH_IMAGE009
is the sampling length of the first direction after the rotation transformation.
2.如权利要求1所述的一种子口径加工工艺的自动优选方法,所述面形误差算术平方根的傅里叶变换
Figure 928984DEST_PATH_IMAGE003
的计算公式为:
Figure 564364DEST_PATH_IMAGE010
,其中,
Figure 720539DEST_PATH_IMAGE011
为面形误差幅值二维矩阵,i为虚数单位,j为第二方向采样点的序号,k为第一方向采样点的序号,j、k为正整数。
2. the automatic optimization method of a kind of sub-caliber processing technology as claimed in claim 1, the Fourier transform of the arithmetic square root of described surface shape error
Figure 928984DEST_PATH_IMAGE003
The calculation formula is:
Figure 564364DEST_PATH_IMAGE010
,in,
Figure 720539DEST_PATH_IMAGE011
is a two-dimensional matrix of surface error amplitudes, i is an imaginary unit, j is the sequence number of the sampling point in the second direction, k is the sequence number of the sampling point in the first direction, and j and k are positive integers.
3.如权利要求1所述的一种子口径加工工艺的自动优选方法,所述旋转变换二维体积谱密度函数
Figure 783173DEST_PATH_IMAGE012
的计算公式为
Figure 392009DEST_PATH_IMAGE013
,其中,rand为旋转变换函数。
3. the automatic optimization method of a kind of sub-caliber processing technology as claimed in claim 1, described rotation transforms two-dimensional volume spectral density function
Figure 783173DEST_PATH_IMAGE012
The calculation formula is
Figure 392009DEST_PATH_IMAGE013
, where rand is the rotation transformation function.
4.如权利要求1所述的一种子口径加工工艺的自动优选方法,其特征在于,步骤S20还包括,通过所述有效去除速率谱和所述体积谱密度函数计算得到时间谱密度函数,对所述时间谱密度函数进行积分得到所述去除函数的加工时间;当所述去除函数的加工时间最小时,所述去除函数为优选加工工艺。4. the automatic optimization method of a kind of sub-caliber processing technology as claimed in claim 1, is characterized in that, step S20 also comprises, obtains time spectral density function by described effective removal rate spectrum and described volume spectral density function calculation, to The time spectral density function is integrated to obtain the processing time of the removal function; when the processing time of the removal function is the smallest, the removal function is a preferred processing technology. 5.如权利要求4所述的一种子口径加工工艺的自动优选方法,其特征在于,步骤S20中,时间谱密度函数为体积谱密度函数与有效去除速率谱相除。5 . The automatic optimization method for a sub-caliber processing technology according to claim 4 , wherein, in step S20 , the time spectral density function is the division of the volume spectral density function and the effective removal rate spectrum. 6 . 6.如权利要求1所述的一种子口径加工工艺的自动优选方法,其特征在于,步骤S20还包括:在预设空间频段对所述体积谱密度函数进行积分得到光学元件在该预设空间频段下的面形误差体积,计算所述去除函数截止频率外的面形误差体积,选择体积差值最小的去除函数作为优选加工工艺。6. The automatic optimization method of a sub-aperture processing technology according to claim 1, wherein step S20 further comprises: integrating the volume spectral density function in a preset space frequency band to obtain the optical element in the preset space The surface shape error volume under the frequency band is calculated, and the surface shape error volume outside the cut-off frequency of the removal function is calculated, and the removal function with the smallest volume difference is selected as the preferred processing technology. 7.一种子口径加工工艺的优选系统,其特征在于,包括有效去除速率谱获取模块、体积谱密度函数获取模块和优选加工工艺获取模块;所述有效去除速率谱获取模块用于获取不同去除函数的有效去除速率谱,通过面形误差体积变化量与加工时间的比值,得到有效去除速率谱;所述体积谱密度函数获取模块用于获取光学元件的体积谱密度函数,所述体积谱密度函数为二维体积谱密度函数或一维体积谱密度函数;所述优选加工工艺获取模块用于通过所述有效去除速率谱和所述体积谱密度函数得到优选加工工艺;7. a preferred system of sub-caliber processing technology, is characterized in that, comprises effective removal rate spectrum acquisition module, volume spectral density function acquisition module and preferred processing technology acquisition module; Described effective removal rate spectrum acquisition module is used to acquire different removal functions The effective removal rate spectrum is obtained through the ratio of the volume change of the surface error to the processing time to obtain the effective removal rate spectrum; the volume spectral density function acquisition module is used to obtain the volume spectral density function of the optical element, and the volume spectral density function is a two-dimensional volume spectral density function or a one-dimensional volume spectral density function; the preferred processing technology acquisition module is used to obtain the preferred processing technology through the effective removal rate spectrum and the volume spectral density function; 所述二维体积谱密度函数
Figure 831080DEST_PATH_IMAGE001
的计算公式为
Figure 107341DEST_PATH_IMAGE002
,其中,S为元件面积,ΔL为采样间隔,x为第一方向,y为第二方向,N为第一方向的采样点数,M为第二方向的采样点数,
Figure 340876DEST_PATH_IMAGE003
为面形误差算术平方根的傅里叶变换,
Figure 437008DEST_PATH_IMAGE004
为第一方向的频域空间坐标,
Figure 675178DEST_PATH_IMAGE005
为第二方向的频域空间坐标;
The two-dimensional volume spectral density function
Figure 831080DEST_PATH_IMAGE001
The calculation formula is
Figure 107341DEST_PATH_IMAGE002
, where S is the element area, ΔL is the sampling interval, x is the first direction, y is the second direction, N is the number of sampling points in the first direction, M is the number of sampling points in the second direction,
Figure 340876DEST_PATH_IMAGE003
is the Fourier transform of the arithmetic square root of the surface error,
Figure 437008DEST_PATH_IMAGE004
is the frequency domain spatial coordinate of the first direction,
Figure 675178DEST_PATH_IMAGE005
is the frequency domain spatial coordinate in the second direction;
所述一维体积谱密度函数
Figure 805945DEST_PATH_IMAGE006
的计算公式为
Figure 210381DEST_PATH_IMAGE007
,其中,
Figure 793809DEST_PATH_IMAGE008
为旋转变换二维体积谱密度函数,θ为旋转角度,
Figure 309104DEST_PATH_IMAGE009
为旋转变换后第一方向的采样长度。
The one-dimensional volume spectral density function
Figure 805945DEST_PATH_IMAGE006
The calculation formula is
Figure 210381DEST_PATH_IMAGE007
,in,
Figure 793809DEST_PATH_IMAGE008
is the rotation transform two-dimensional volume spectral density function, θ is the rotation angle,
Figure 309104DEST_PATH_IMAGE009
is the sampling length of the first direction after the rotation transformation.
8.一种可读存储介质,其特征在于,用于存储程序,所述程序被执行时,用于实现如权利要求1-6之一所述的子口径加工工艺的自动优选方法。8 . A readable storage medium, characterized in that it is used for storing a program, and when the program is executed, it is used to realize the automatic optimization method of the sub-caliber machining process according to any one of claims 1 to 6 .
CN202210355185.2A 2022-04-06 2022-04-06 Automatic optimization method, system and medium for sub-caliber processing technology Active CN114425732B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210355185.2A CN114425732B (en) 2022-04-06 2022-04-06 Automatic optimization method, system and medium for sub-caliber processing technology

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210355185.2A CN114425732B (en) 2022-04-06 2022-04-06 Automatic optimization method, system and medium for sub-caliber processing technology

Publications (2)

Publication Number Publication Date
CN114425732A CN114425732A (en) 2022-05-03
CN114425732B true CN114425732B (en) 2022-06-03

Family

ID=81314403

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210355185.2A Active CN114425732B (en) 2022-04-06 2022-04-06 Automatic optimization method, system and medium for sub-caliber processing technology

Country Status (1)

Country Link
CN (1) CN114425732B (en)

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101240998A (en) * 2008-03-14 2008-08-13 中国人民解放军国防科学技术大学 Analysis method of frequency band error distribution characteristics under the condition of deterministic optical processing
CN102848287A (en) * 2012-09-14 2013-01-02 中国人民解放军国防科学技术大学 Combination machining method for removing high-frequency errors in optical elements
CN103144004A (en) * 2013-03-22 2013-06-12 哈尔滨工业大学 Edge precision control method of large aperture optical element being processed through air bag polishing
CN103342476A (en) * 2013-07-03 2013-10-09 中国科学院光电技术研究所 Ion beam sacrificial layer processing method for suppressing high frequency error in optical surface
CN103395000A (en) * 2013-07-25 2013-11-20 中国科学院光电技术研究所 Evaluation method for inhibiting error capability of different frequency bands by CCOS polishing process
CN103517784A (en) * 2011-05-13 2014-01-15 依视路国际集团(光学总公司) Method for determining positional parameters of a fabrication surface relative to a reference surface
CN105328535A (en) * 2015-09-29 2016-02-17 中国人民解放军国防科学技术大学 Nanometer-precision optical curved-face ion beam processing method based on non-linear modeling
CN109227226A (en) * 2018-11-12 2019-01-18 中国科学院光电技术研究所 Uniform-sliding method for residence time in optical element processing process
CN109909815A (en) * 2019-03-28 2019-06-21 中国人民解放军国防科技大学 Magnetorheological polishing compensation processing method, system and medium for optical complex curved surface element
CN110059407A (en) * 2019-04-18 2019-07-26 中国科学院光电技术研究所 A kind of method of improved evaluation CCOS polishing removal function error rejection ability
CN110245317A (en) * 2019-05-16 2019-09-17 中国工程物理研究院激光聚变研究中心 A kind of extracting method and device of Magnetorheological Polishing removal function
CN112428026A (en) * 2020-11-13 2021-03-02 中国人民解放军国防科技大学 Pulse control beam diameter adjustable ion beam processing method based on surface shape error frequency band

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6709314B2 (en) * 2001-11-07 2004-03-23 Applied Materials Inc. Chemical mechanical polishing endpoinat detection
US6866793B2 (en) * 2002-09-26 2005-03-15 University Of Florida Research Foundation, Inc. High selectivity and high planarity dielectric polishing
US8992286B2 (en) * 2013-02-26 2015-03-31 Applied Materials, Inc. Weighted regression of thickness maps from spectral data
US11451735B2 (en) * 2018-11-20 2022-09-20 Teledyne Flir, Llc High dynamic range micromirror imaging array systems and methods

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101240998A (en) * 2008-03-14 2008-08-13 中国人民解放军国防科学技术大学 Analysis method of frequency band error distribution characteristics under the condition of deterministic optical processing
CN103517784A (en) * 2011-05-13 2014-01-15 依视路国际集团(光学总公司) Method for determining positional parameters of a fabrication surface relative to a reference surface
CN102848287A (en) * 2012-09-14 2013-01-02 中国人民解放军国防科学技术大学 Combination machining method for removing high-frequency errors in optical elements
CN103144004A (en) * 2013-03-22 2013-06-12 哈尔滨工业大学 Edge precision control method of large aperture optical element being processed through air bag polishing
CN103342476A (en) * 2013-07-03 2013-10-09 中国科学院光电技术研究所 Ion beam sacrificial layer processing method for suppressing high frequency error in optical surface
CN103395000A (en) * 2013-07-25 2013-11-20 中国科学院光电技术研究所 Evaluation method for inhibiting error capability of different frequency bands by CCOS polishing process
CN105328535A (en) * 2015-09-29 2016-02-17 中国人民解放军国防科学技术大学 Nanometer-precision optical curved-face ion beam processing method based on non-linear modeling
CN109227226A (en) * 2018-11-12 2019-01-18 中国科学院光电技术研究所 Uniform-sliding method for residence time in optical element processing process
CN109909815A (en) * 2019-03-28 2019-06-21 中国人民解放军国防科技大学 Magnetorheological polishing compensation processing method, system and medium for optical complex curved surface element
CN110059407A (en) * 2019-04-18 2019-07-26 中国科学院光电技术研究所 A kind of method of improved evaluation CCOS polishing removal function error rejection ability
CN110245317A (en) * 2019-05-16 2019-09-17 中国工程物理研究院激光聚变研究中心 A kind of extracting method and device of Magnetorheological Polishing removal function
CN112428026A (en) * 2020-11-13 2021-03-02 中国人民解放军国防科技大学 Pulse control beam diameter adjustable ion beam processing method based on surface shape error frequency band

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
强光光学零件加工误差频谱分析与控制方法研究;李富仁等;《精密加工》;20130830;第49卷(第4期);第1-4页 *

Also Published As

Publication number Publication date
CN114425732A (en) 2022-05-03

Similar Documents

Publication Publication Date Title
CN101456680B (en) Processing method for correcting low steepness optical mirror surface error
CN104143210B (en) Multi-scale normal feature point cloud registering method
CN109887015A (en) A Point Cloud Automatic Registration Method Based on Local Surface Feature Histogram
CN109982361B (en) Signal interference analysis method, device, equipment and medium
CN117094232B (en) Method and system for updating deep oil gas accurate navigation three-dimensional lithology model in real time
CN107369140A (en) High-accuracy target ball center extraction method under unstructured moving grids
CN114663373A (en) Point cloud registration method and device for detecting surface quality of part
CN119740058B (en) Background signal extraction method and device
CN114425732B (en) Automatic optimization method, system and medium for sub-caliber processing technology
CN120351954A (en) Calibration method and system of high-precision inclination angle measurement sensor
CN108228534A (en) Curve fitting test method
CN116542944B (en) Wafer warpage detection method, device, equipment and medium
CN110458106A (en) A kind of intelligent analysis method and intellectual analysis device of tomato growth state
CN106873315B (en) A kind of via layer OPC modeling methods
CN119579581B (en) Optical element surface shape full-frequency-band error extraction method based on neural network
CN111913945A (en) Data management method and device and storage medium
CN110705132A (en) Guidance control system performance fusion evaluation method based on multi-source heterogeneous data
EP4105672A1 (en) Systems and methods for provisioning training data to enable neural networks to analyze signals in nmr measurements
CN118747725B (en) A LiDAR Point Cloud Noise Reduction Method Based on Improved DBSCAN
CN113343492A (en) Theoretical spectral data optimization method and system and optical measurement method
CN114662047B (en) Characterization method, system and medium for correction capability of removing function error
CN114722643B (en) Virtual-real consistency verification method based on complex system simulation model
CN112200252A (en) Joint dimension reduction method based on probability box global sensitivity analysis and active subspace
US12057336B2 (en) Estimating heights of defects in a wafer by scaling a 3D model using an artificial neural network
CN105488521B (en) A kind of dilatation screening sample method based on kernel function

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