CN105423948A - Splicing-interference-detection aspheric surface shape apparatus using distorting lens and method thereof - Google Patents

Splicing-interference-detection aspheric surface shape apparatus using distorting lens and method thereof Download PDF

Info

Publication number
CN105423948A
CN105423948A CN201510922163.XA CN201510922163A CN105423948A CN 105423948 A CN105423948 A CN 105423948A CN 201510922163 A CN201510922163 A CN 201510922163A CN 105423948 A CN105423948 A CN 105423948A
Authority
CN
China
Prior art keywords
sub
aperture
distorting lens
aspheric surface
interferometer
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201510922163.XA
Other languages
Chinese (zh)
Other versions
CN105423948B (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.)
Changchun Institute of Optics Fine Mechanics and Physics of CAS
Original Assignee
Changchun Institute of Optics Fine Mechanics and Physics of CAS
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 Changchun Institute of Optics Fine Mechanics and Physics of CAS filed Critical Changchun Institute of Optics Fine Mechanics and Physics of CAS
Priority to CN201510922163.XA priority Critical patent/CN105423948B/en
Publication of CN105423948A publication Critical patent/CN105423948A/en
Priority to PCT/CN2016/100862 priority patent/WO2017101557A1/en
Application granted granted Critical
Publication of CN105423948B publication Critical patent/CN105423948B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/24Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/24Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
    • G01B11/2441Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures using interferometry

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Instruments For Measurement Of Length By Optical Means (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

A splicing-interference-detection aspheric surface shape apparatus using a distorting lens and a method thereof belong to the optical detection technology field. In the invention, an interference detection system, a distorting lens compensation system, a mechanical adjusting mechanism and a computer data processing module are included. The method comprises the following steps of step1, constructing a splicing interference detection apparatus using the distorting lens; step2, dividing a subaperture; step3, removing to a position of the subaperture to be detected and adjusting a distorting lens surface shape or a refractive index according to modeling and interference fringes; step4, carrying out subaperture interferogram acquisition processing; step5, carrying out subaperture backhaul error correction; and step6, carrying out full aperture surface shape splicing. In the invention, the distorting lens and the subaperture splicing interference detection method are combined; the subaperture number needed by full aperture coverage can be effectively reduced and effective areas of each subaperture and an overlapping region are increased; a problem that the overlapping region is small so that splicing precision is influenced is solved.

Description

The stitching interferometer of distorting lens is adopted to detect the device and method of aspheric surface
Technical field
The invention belongs to technical field of optical detection, be specifically related to a kind of device and method adopting the stitching interferometer of distorting lens to detect aspheric surface.
Background technology
Interfere detection technique to be conventional surface testing technology, it is tested surface and the relative detection of reference surface.Reference surface is all generally plane or sphere, and there will not be more interference fringe when the plane after detecting polishing with them and sphere, and detect aspheric surface with spherical reference face, when aspheric slope is larger, interference fringe density is also more.When interference fringe density exceedes the nyquist frequency of interferometer detector, interference fringe cannot be detected device and differentiate, and causes aspheric surface cannot be measured.For the aspheric surface that slope is larger, adopt sub-aperture stitching interferometer detection method, aspheric surface larger for slope is divided into multiple small-bore sub-aperture, makes the interference fringe in each sub-aperture to be detected resolution.
But when aspheric surface bore slope is excessive, the area that each sub-aperture from center to edge can be detected is more and more less, now need more sub-aperture could cover whole aspheric surface, reduce detection efficiency, add machine error, detection time and data processing amount etc.
Summary of the invention
The object of the invention is to propose a kind of device and method adopting the stitching interferometer of distorting lens to detect aspheric surface, solve the problem that detection efficiency is low, machine error is large, detection time is long and data processing amount is large that prior art exists.In conjunction with distorting lens technology and annular sub-aperture stitching interferometer detection method, effectively can reduce and cover unified required sub-aperture number, increase each sub-aperture useful area, solve and the difficult problem of impact splicing precision little due to overlay region.
For achieving the above object, the device of the stitching interferometer detection aspheric surface of employing distorting lens of the present invention comprises distorting lens bucking-out system, interferes detection system, mechanical adjusting mechanism and computer digital animation module;
Described computer digital animation module comprises distorting lens control module, multidimensional adjusting table control module, interferogram sampling processing unit, hysterisis error correcting unit and sub-aperture stitching algorithm unit;
Aspheric surface to be detected is fixed on described mechanical adjusting mechanism, and described interference detection system compensates the optical path difference with aspheric surface reflected light to be detected by distorting lens bucking-out system, and described interferogram sampling processing unit is connected with the detector of interfering in detection system;
Described interferogram sampling processing unit is connected with described distorting lens control module and hysterisis error correcting unit, distorting lens control module carries out close-loop feedback according to interferogram to the adjustment of distorting lens, control multidimensional adjusting table control module and close-loop feedback control is carried out to the motion of multidimensional adjusting table, each sub-aperture interferogram that detector collects is through described interferogram sampling processing unit processes, obtain each sub-aperture and return Wave-front phase, each sub-aperture returns Wave-front phase and obtains each sub-aperture diametric plane shape information through hysterisis error correcting unit, each sub-aperture diametric plane shape information obtains unified shape information through the process of sub-aperture stitching algorithm unit.
Described mechanical adjusting mechanism comprises clamping device and multidimensional adjusting table, and described aspheric surface to be detected is fixed on described multidimensional adjusting table by clamping device.
Described interference detection system is Tai Man-Green's type interferometer or Feisuo type interferometer.
The method adopting the stitching interferometer of distorting lens to detect aspheric surface comprises the following steps:
Step one: build the device adopting the stitching interferometer of distorting lens to detect aspheric surface;
Step 2: system modelling, divides sub-aperture, and makes the reference paper of each sub-aperture, be specially:
1) sub-aperture planning is carried out according to the bore of the radius-of-curvature of aspheric best-fit ball to be detected, best-fit ball, the radius-of-curvature with reference to sphere, the bore with reference to sphere and sub-aperture overlapping region;
2) interfere the parameter of the concrete device of detection system according to aspheric surface in pick-up unit, in optical design software, interfere detection system to carry out modeling to aspheric surface, obtain system model;
3) according to step 1) planning carried out sub-aperture is in step 2) in carry out sub-aperture division in the system model that obtains;
4) detect the known high precision sphere of face shape error ε by interferometer, obtain the measured value W of high precision spherical surface shape sphere-test; In Zemax optical design software, faced by simulate ideal reference sphere, high precision sphere detects, and high precision sphere surface figure accuracy testing result is W sphere-ideal, obtain interferometer system error E by formula (1) irigging error E is manufactured with collimating mirror f:
E I+E F=W sphere-teat-W sphere-ideal-2ε(1)
5) by any one sub-aperture diametric plane shape of system model analog detection desired aspheric W asphere-ideal; By any one sub-aperture diametric plane shape W asphere-idealdepart from the error caused with reference to sphere as aspheric surface, obtain the fixed phase W of sub-aperture according to formula (2) reference:
W reference=W asphere-ideal+W sphere-test-W sphere-ideal-2ε(2);
Step 3: divide according to sub-aperture and show that the amount of exercise detecting each axle of multidimensional adjusting table needed for each sub-aperture is accurately located aspheric surface, and adjust distoring mirror shape or refractive index, until interference fringe picture meets testing requirement according to system model and interference fringe;
Step 4: gather each sub-aperture interferogram by interferogram sampling processing unit, utilizes unwrapping algorithm to carry out interferogram phase demodulating, obtains each sub-aperture Wave-front phase W that in testing, detector receives asphere-test;
Step 5: calculating each sub-aperture diametric plane shape ε according to formula (3) asphere:
ϵ a s p h e r e = 1 2 [ W a s p h e r e - t e s t - W r e f e r e n c e ] (3)
Wherein: W asphere-testfor sub-aperture Wave-front phase
W referencefor the fixed phase of sub-aperture;
Step 6: according to each sub-aperture diametric plane shape ε obtained in step 5 asphere, reduce the impact of machine error on splicing result by stitching algorithm, the most each sub-aperture diametric plane graphic data splicing, draws unified shape.
Described interference detection system is Tai Man-Green's type interferometer or Feisuo type interferometer.
Described in step 3 according to system model and interference fringe adjustment distoring mirror shape or refractive index, until the concrete steps that interference fringe picture meets testing requirement are:
1) detecting the Zelnick aberration at the detector place obtained according to emulation in optical design software, when interfering detection system to select Twyman Green Interferometer, the face shape of distorting lens being set to 1/2 of Zelnick aberration; When interfering detection system to select Feisuo interferometer, adopt the method for " inversion " to calculate the index distribution compensated needed for aspheric surface, the phase place needed for calculation compensation aspheric surface, P is any point in aspheric surface, cross P and do aspheric normal, hand over distorting lens rear surface to be M; Light PM, after distorting lens diffraction, hands over distorting lens front surface to be N, then to reflect last convergent point be F; Wherein the refractive index of distorting lens is ng, and the refractive index of surrounding air is na; Aspheric surface is d1 apart from the distance of distorting lens rear surface, and distorting lens thickness is the distance of d, F distance distorting lens front surface is d2;
According to Fermat principle, any point in aspheric surface be made all to converge to a F along the light of normal emergence, then have:
n a|PM|+n g|MN|+n a|NF|=n ad 1+n gd+n ad 2
Distorting lens refractive index ng is calculated by formula (4):
n g = n a d 1 + n g d + n a d 2 - n a | N F | - n a | P M | | M N | - - - ( 4 )
2) by interfering detection system to detect sub-aperture, interference fringe picture is obtained;
3) according to step 2) in the interference fringe picture that obtains distorting lens carry out step 1) in feedback regulation, until interference fringe picture meets testing requirement, distorting lens has regulated.
Stitching algorithm described in step 6 is specially:
1) according to sub-aperture diametric plane shape after formula (5) calculating introducing compensation rate
E asphere i = ϵ asphere i + Σ k = 1 L F i k f k ( x , y ) , k = 1 , 2 ... L - - - ( 5 )
Wherein: f kthe compensation factor that (x, y) is sub-aperture;
L is the number compensating factor;
for the face shape of i-th sub-aperture that step 5 obtains;
F ikfor sub-aperture compensates the penalty coefficient of factor, obtained the penalty coefficient F of each sub-aperture by the least square method shown in formula (6) ik;
min = Σ i = 1... N - 1 Σ j = 0... N - 1 j ∩ i ( ( ϵ asphere i ( x , y ) + Σ k = 1 L F i k f k ( x , y ) ) - ( ϵ asphere j ( x , y ) + Σ k = 1 L F j k f k ( x , y ) ) ) 2 - - - ( 6 )
Wherein, N is the number of sub-aperture;
Sub-aperture counts from zero, i.e. sub-aperture centered by the 0th sub-aperture;
representing only has the jth of an overlapping region sub-aperture to calculate to i-th sub-aperture;
2) according to step 1) each sub-aperture diametric plane shape of obtaining splicing obtains unified shape, is about to all be added, the face shape of overlapping region gets the weighted mean value of each sub-aperture data, and weight is relevant to the accuracy of measurement data.
Described optical design software is Zemax software.
Beneficial effect of the present invention is: the present invention is in conjunction with distorting lens and sub-aperture stitching interferometer detection method, distorting lens is utilized to produce the different sub-aperture region of non-spherical wavefront coupling tested surface, make to detect unified required sub-aperture number to greatly reduce, add sub-aperture area, thus add the overlay region of adjacent sub-aperture, reduce the accumulation of detection time, machine error, improve the efficiency and precision that detect.Efficiently solve the problem that the splicing precision that causes because sub-aperture area is too low is low, add splicing efficiency simultaneously.
Accompanying drawing explanation
Fig. 1 is that the stitching interferometer of employing distorting lens of the present invention detects in the device of aspheric surface structural representation when adopting Tai Man-Green's type interferometer;
Fig. 2 is that the stitching interferometer of employing distorting lens of the present invention detects in the device of aspheric surface structural representation when adopting Feisuo type interferometer;
Fig. 3 is the method overhaul flow chart of the stitching interferometer detection aspheric surface of employing distorting lens of the present invention;
Fig. 4 is the method sub-aperture Structure Planning schematic diagram of the stitching interferometer detection aspheric surface of employing distorting lens of the present invention;
Fig. 5 is that the stitching interferometer of employing distorting lens of the present invention detects distoring mirror shape or refractive index adjust structure schematic diagram in the method for aspheric surface;
Wherein: 1, frequency stabilized laser, 2, collimating and beam expanding system, 3, Amici prism, 4, distorting lens, 5, imaging lens, 6, detector, 7, converging lenses group, 8, aspheric surface, 9, clamping device, 10, computing machine, 11, multidimensional adjusting table, 12, distorting lens control module, 13, multidimensional adjusting table control module, 14, interferogram sampling processing unit, 15, catoptron, 16, standard mirror, 17, reference surface.
Embodiment
Below in conjunction with accompanying drawing, embodiments of the present invention are described further.
The device of the stitching interferometer detection aspheric surface of employing distorting lens of the present invention comprises distorting lens bucking-out system, interferes detection system, mechanical adjusting mechanism and computer digital animation module;
Described computer digital animation module comprises distorting lens control module 12, multidimensional adjusting table control module 13, interferogram sampling processing unit 14, hysterisis error correcting unit and sub-aperture stitching algorithm unit;
Aspheric surface 8 to be detected is fixed on described mechanical adjusting mechanism, described interference detection system compensates the optical path difference with aspheric surface 8 reflected light to be detected by distorting lens bucking-out system, and described interferogram sampling processing unit 14 is connected with the detector 6 of interfering in detection system;
Described interferogram sampling processing unit 14 is connected with described distorting lens control module 12 and hysterisis error correcting unit, distorting lens control module 12 carries out close-loop feedback according to interferogram to the adjustment of distorting lens 4, close-loop feedback control is carried out in the motion controlling multidimensional adjusting table control module 13 pairs of multidimensional adjusting tables 11, each sub-aperture interferogram that detector 6 collects processes through described interferogram sampling processing unit 14, obtain each sub-aperture and return Wave-front phase, each sub-aperture returns Wave-front phase and obtains each sub-aperture diametric plane shape information through hysterisis error correcting unit, each sub-aperture diametric plane shape information obtains unified shape information through the process of sub-aperture stitching algorithm unit.
Described mechanical adjusting mechanism comprises clamping device 9 and multidimensional adjusting table 11, and described aspheric surface 8 to be detected is fixed on described multidimensional adjusting table 11 by clamping device 9.
See accompanying drawing 1, described interference detection system is Tai Man-Green's type interferometer, and the light pencil of frequency stabilized laser 1 outgoing is expanded as angle pencil of ray directional light through collimating and beam expanding system 2, and angle pencil of ray directional light forward direction is divided into two-way to Amici prism 3 place; Wherein a road returns as reference ripple through Amici prism 3 reflections propagate to distorting lens 4 Hou Yuan road; Another road is first assembled and is dispersed afterwards after Amici prism 3 transmission forward direction to converging lenses group 7, and diverging light again after converging lenses group 7, is formed and detects ripple after aspheric surface 8 reflects.Detector 6 place is imaged in through imaging lens 5 after reference wave and detection wave interference.
See accompanying drawing 2, described interference detection system is Feisuo type interferometer, the light pencil of the frequency stabilized laser 1 outgoing light forward direction after collimating and beam expanding system 2 is to Amici prism 3 place, Amici prism 3 reflections propagate is to catoptron 15, through standard mirror 16, a part is as reference ripple after reference surface 17 reflects, and another part transmission forward direction is to distorting lens 4, after aspheric surface 8 reflects again after distorting lens 4, formed and detect ripple.Detector 6 place is imaged in through imaging lens 5 after reference wave and detection wave interference.
The need of work distorting lens drive system of distorting lens bucking-out system, interference detection system, computer control system complete jointly.By the interference fringe picture of interferometer measurement, close-loop feedback control is carried out to distorting lens 4.Distorting lens bucking-out system mainly comprises the distorting lens control module 12 in distorting lens 4 and computer digital animation module, draw required distorting lens 4 deflection according to computing machine 10 simulation results or interference detection system measurement result, by computing machine 10 controlling distortion mirror driver, distorting lens 4 is out of shape.
Mechanical adjusting mechanism comprises the multidimensional adjusting table control module 13 in the locating devices such as multidimensional adjusting table 11, grating scale and computer digital animation module.Aspheric surface 8 is fixed on clamping device 9, and clamping device 9 is fixed on multidimensional adjusting table 11.The mirror pose to be measured providing and detect needed for each sub-aperture is planned by sub-aperture, computing machine 10 can drive multidimensional adjusting table 11 to make clamping device 9 along the translation of x, y, z axle, x, y-axis tilt, and rotate around z-axis, the position of multi-dimensional movement can be carried out close-loop feedback by locating devices such as all directions grating scales and be controlled.
The image that detector 6 collects by interferogram sampling processing unit 14 is after machine 10 processes as calculated, and Output rusults, to hysterisis error correcting unit, then obtains unified shape information through sub-aperture stitching unit.
See accompanying drawing 3, the method adopting the stitching interferometer of distorting lens to detect aspheric surface comprises the following steps:
Step one: build the device adopting the stitching interferometer of distorting lens to detect aspheric surface;
Step 2: system modelling, divides sub-aperture, and makes the reference paper of each sub-aperture, be specially:
1) sub-aperture planning is carried out according to the bore of the radius-of-curvature of the best-fit ball of aspheric surface 8 to be detected, best-fit ball, the radius-of-curvature with reference to sphere, the bore with reference to sphere and sub-aperture overlapping region;
See accompanying drawing 4, such as aspheric surface 8 best-fit ball, reference sphere, overlapping region are as following table:
Sub-aperture distribution is as follows:
2) interfere the parameter of the concrete device of detection system according to aspheric surface in pick-up unit 8, in optical design software, interfere detection system to carry out modeling to aspheric surface 8, obtain system model;
Detection system is interfered for Tai Man-Green's type, the parameter of described concrete device comprises: the optical maser wavelength of frequency stabilized laser 1 outgoing, the clear aperture of collimating and beam expanding system 2 and enlargement factor, the parameter of Amici prism 3, the bore of distorting lens 4, the bore of imaging lens 5 and surface curvature radius and thickness, the nominal face shape equation of tested aspheric surface 8 and bore;
3) according to step 1) planning carried out sub-aperture is in step 2) in carry out sub-aperture division in the system model that obtains;
4) detect the known high precision sphere of face shape error ε by interferometer, obtain the measured value W of high precision spherical surface shape sphere-test; In Zemax optical design software, faced by simulate ideal reference sphere, high precision sphere detects, and high precision sphere surface figure accuracy testing result is W sphere-ideal, obtain interferometer system error E by formula (1) irigging error E is manufactured with collimating mirror f:
E I+E F=W sphere-test-W sphere-ideal-2ε(1)
5) by any one sub-aperture diametric plane shape of system model analog detection desired aspheric 8 W asphere-ideal; By any one sub-aperture diametric plane shape W asphere-idealdepart from the error caused with reference to sphere as aspheric surface 8, obtain the fixed phase W of sub-aperture according to formula (2) reference:
W reference=W asphere-ideal+W sphere-test-W sphere-ideal-2ε(2);
Step 3: divide according to sub-aperture and show that the amount of exercise detecting each axle of multidimensional adjusting table 11 needed for each sub-aperture is accurately located aspheric surface 8, and adjust distorting lens 4 shapes or refractive index, until interference fringe picture meets testing requirement according to system model and interference fringe;
Step 4: gather each sub-aperture interferogram by interferogram sampling processing unit 14, utilizes unwrapping algorithm to carry out interferogram phase demodulating, obtains each sub-aperture Wave-front phase W that in testing, detector 6 receives asphere-test;
Step 5: calculating each sub-aperture diametric plane shape ε according to formula (3) asphere:
ϵ a s p h e r e = 1 2 [ W a s p h e r e - t e s t - W r e f e r e n c e ] - - - ( 3 )
Wherein: W asphere-testfor sub-aperture Wave-front phase
W referencefor the fixed phase of sub-aperture;
Step 6: according to each sub-aperture diametric plane shape ε obtained in step 5 asphere, reduce the impact of machine error on splicing result by stitching algorithm, the most each sub-aperture diametric plane graphic data splicing, draws unified shape.
Described interference detection system is Tai Man-Green's type interferometer or Feisuo type interferometer.
Described in step 3 according to system model and interference fringe adjustment distorting lens 4 shapes or refractive index, until the concrete steps that interference fringe picture meets testing requirement are:
1) detecting the Zelnick aberration at detector 6 place obtained according to emulation in optical design software, when interfering detection system to select Twyman Green Interferometer, the face shape of distorting lens 4 being set to 1/2 of Zelnick aberration; When interfering detection system to select Feisuo interferometer, the method of " inversion " is adopted to calculate the index distribution compensated needed for aspheric surface 8, phase place needed for calculation compensation aspheric surface 8, see accompanying drawing 5, P is any point in aspheric surface 8, cross the normal that P does aspheric surface 8, hand over distorting lens 4 rear surface to be M; Light PM, after distorting lens 4 diffraction, hands over distorting lens 4 front surface to be N, then to reflect last convergent point be F; Wherein the refractive index of distorting lens 4 is n g, the refractive index of surrounding air is n a; Aspheric surface 8 is d apart from the distance of distorting lens 4 rear surface 1, distorting lens 4 thickness is the distance of d, F distance distorting lens 4 front surface is d 2;
According to Fermat principle, any point in aspheric surface 8 be made all to converge to a F along the light of normal emergence, then have:
n a|PM|+n g|MN|+n a|NF|=n ad 1+n gd+n ad 2
Distorting lens 4 refractive index n is calculated by formula (4) g:
n g = n a d 1 + n g d + n a d 2 - n a | N F | - n a | P M | | M N | - - - ( 4 )
2) by interfering detection system to detect sub-aperture, interference fringe picture is obtained;
3) according to step 2) in the interference fringe picture that obtains distorting lens 4 carry out step 1) in feedback regulation, until interference fringe picture meets testing requirement, distorting lens 4 has regulated.
Stitching algorithm described in step 6 is specially:
1) according to sub-aperture diametric plane shape after formula (5) calculating introducing compensation rate
E asphere i = ϵ asphere i + Σ k = 1 L F i k f k ( x , y ) , k = 1 , 2 ... L - - - ( 5 )
Wherein: f kthe compensation factor that (x, y) is sub-aperture;
L is the number compensating factor;
for the face shape of i-th sub-aperture that step 5 obtains;
F ikfor sub-aperture compensates the penalty coefficient of factor, obtain the penalty coefficient F of each sub-aperture by the such as least square method shown in formula (6) ik;
min = Σ i = 1... N - 1 Σ j = 0... N - 1 j ∩ i ( ( ϵ asphere i ( x , y ) + Σ k = 1 L F i k f k ( x , y ) ) - ( ϵ asphere j ( x , y ) + Σ k = 1 L F j k f k ( x , y ) ) ) 2 - - - ( 6 )
Wherein, N is the number of sub-aperture;
Sub-aperture counts from zero, i.e. sub-aperture centered by the 0th sub-aperture;
representing only has the jth of an overlapping region sub-aperture to calculate to i-th sub-aperture;
2) according to step 1) in each sub-aperture diametric plane shape of obtaining splicing obtains unified shape, by used be added, the face shape of overlapping region gets the weighted mean value of each sub-aperture data, and weight is relevant to the accuracy of measurement data.
Described optical design software is Zemax software, is the optical design software developed by RadiantZemax company of the U.S..

Claims (8)

1. adopt the stitching interferometer of distorting lens to detect the device of aspheric surface, it is characterized in that, comprise distorting lens bucking-out system, interfere detection system, mechanical adjusting mechanism and computer digital animation module;
Described computer digital animation module comprises distorting lens control module (12), multidimensional adjusting table control module (13), interferogram sampling processing unit (14), hysterisis error correcting unit and sub-aperture stitching algorithm unit;
Aspheric surface to be detected (8) is fixed on described mechanical adjusting mechanism, described interference detection system compensates the optical path difference with aspheric surface to be detected (8) reflected light by distorting lens bucking-out system, and described interferogram sampling processing unit (14) is connected with the detector (6) of interfering in detection system;
Described interferogram sampling processing unit (14) is connected with described distorting lens control module (12) and hysterisis error correcting unit, distorting lens control module (12) carries out close-loop feedback according to the adjustment of interferogram to distorting lens (4), control multidimensional adjusting table control module (13) motion to multidimensional adjusting table (11) and carry out close-loop feedback control, each sub-aperture interferogram that detector (6) collects processes through described interferogram sampling processing unit (14), obtain each sub-aperture and return Wave-front phase, each sub-aperture returns Wave-front phase and obtains each sub-aperture diametric plane shape information through hysterisis error correcting unit, each sub-aperture diametric plane shape information obtains unified shape information through the process of sub-aperture stitching algorithm unit.
2. the stitching interferometer of employing distorting lens according to claim 1 detects the device of aspheric surface, it is characterized in that, described mechanical adjusting mechanism comprises clamping device (9) and multidimensional adjusting table (11), and described aspheric surface to be detected (8) is fixed on described multidimensional adjusting table (11) by clamping device (9).
3. the stitching interferometer of employing distorting lens according to claim 1 detects the device of aspheric surface, and it is characterized in that, described interference detection system is Tai Man-Green's type interferometer or Feisuo type interferometer.
4. the stitching interferometer of employing distorting lens according to claim 1 detects the detection method of the device of aspheric surface, it is characterized in that, comprises the following steps:
Step one: build the device adopting the stitching interferometer of distorting lens to detect aspheric surface;
Step 2: system modelling, divides sub-aperture, and makes the reference paper of each sub-aperture, be specially:
1) sub-aperture planning is carried out according to the bore of the radius-of-curvature of the best-fit ball of aspheric surface to be detected (8), best-fit ball, the radius-of-curvature with reference to sphere, the bore with reference to sphere and sub-aperture overlapping region;
2) interfere the parameter of the concrete device of detection system according to aspheric surface in pick-up unit, in optical design software, interfere detection system to carry out modeling to aspheric surface, obtain system model;
3) according to step 1) planning carried out sub-aperture is in step 2) in carry out sub-aperture division in the system model that obtains;
4) detect the known high precision sphere of face shape error ε by interferometer, obtain the measured value W of high precision spherical surface shape sphere-test; In Zemax optical design software, faced by simulate ideal reference sphere, high precision sphere detects, and high precision sphere surface figure accuracy testing result is W sphere-ideal, obtain interferometer system error E by formula (1) irigging error E is manufactured with collimating mirror f:
E I+E F=W sphere-test-W sphere-ideal-2ε(1)
5) by any one sub-aperture diametric plane shape of system model analog detection desired aspheric (8) W asphere-ideal; By any one sub-aperture diametric plane shape W asphere-idealdepart from the error caused with reference to sphere as aspheric surface (8), obtain the fixed phase W of sub-aperture according to formula (2) reference:
W reference=W asphere-ideal+W sphere-test-W sphere-ideal-2ε(2);
Step 3: divide according to sub-aperture and show that the amount of exercise detecting multidimensional adjusting table (11) each axle needed for each sub-aperture is accurately located aspheric surface (8), and adjust distorting lens (4) face shape or refractive index, until interference fringe picture meets testing requirement according to system model and interference fringe;
Step 4: gather each sub-aperture interferogram by interferogram sampling processing unit (14), utilizes unwrapping algorithm to carry out interferogram phase demodulating, obtains each sub-aperture Wave-front phase W that in testing, detector (6) receives asphere-test;
Step 5: calculating each sub-aperture diametric plane shape ε according to formula (3) asphere:
ϵ a s p h e r e = 1 2 [ W a s p h e r e - t e s t - W r e f e r e n c e ] - - - ( 3 )
Wherein: W asphere-testfor sub-aperture Wave-front phase
W referencefor the fixed phase of sub-aperture;
Step 6: according to each sub-aperture diametric plane shape ε obtained in step 5 asphere, reduce the impact of machine error on splicing result by stitching algorithm, the most each sub-aperture diametric plane graphic data splicing, draws unified shape.
5. the stitching interferometer of employing distorting lens according to claim 4 detects the device of aspheric surface, and it is characterized in that, described interference detection system is Tai Man-Green's type interferometer or Feisuo type interferometer.
6. the stitching interferometer of employing distorting lens according to claim 5 detects the device of aspheric surface, it is characterized in that, described in step 3 according to system model and interference fringe adjustment distorting lens (4) face shape or refractive index, until the concrete steps that interference fringe picture meets testing requirement are:
1) the Zelnick aberration at detector (6) place obtained is detected according to emulation in optical design software, when interfering detection system to select Twyman Green Interferometer, the face shape of distorting lens (4) is set to 1/2 of Zelnick aberration; When interfering detection system to select Feisuo interferometer, the method of " inversion " is adopted to calculate the index distribution compensated needed for aspheric surface (8), phase place needed for calculation compensation aspheric surface (8), P is any point in aspheric surface (8), cross the normal that P does aspheric surface (8), hand over distorting lens (4) rear surface to be M; Light PM, after distorting lens (4) diffraction, hand over distorting lens (4) front surface to be N, then to reflect last convergent point is F; Wherein the refractive index of distorting lens (4) is n g, the refractive index of surrounding air is n a; Aspheric surface (8) is d apart from the distance of distorting lens (4) rear surface 1, distorting lens (4) thickness is the distance of d, F distance distorting lens (4) front surface is d 2;
According to Fermat principle, the upper any point of aspheric surface (8) be made all to converge to a F along the light of normal emergence, then have:
n a|PM|+n g|MN|+n a|NF|=n ad 1+n gd+n ad 2
Distorting lens (4) refractive index ng is calculated by formula (4):
n g = n a d 1 + n g d + n a d 2 - n a | N F | - n a | P M | | M N | - - - ( 4 )
2) by interfering detection system to detect sub-aperture, interference fringe picture is obtained;
3) according to step 2) in the interference fringe picture that obtains distorting lens (4) carry out step 1) in feedback regulation, until interference fringe picture meets testing requirement, distorting lens (4) has regulated.
7. the stitching interferometer of employing distorting lens according to claim 4 detects the device of aspheric surface, and it is characterized in that, the stitching algorithm described in step 6 is specially:
1) according to sub-aperture diametric plane shape after formula (5) calculating introducing compensation rate
E a s p h e r e i = ϵ a s p h e r e i + Σ k = 1 L F i k f k ( x , y ) , k = 1 , 2 ... L - - - ( 5 )
Wherein: f kthe compensation factor that (x, y) is sub-aperture;
L is the number compensating factor;
for the face shape of i-th sub-aperture that step 5 obtains;
F ikfor sub-aperture compensates the penalty coefficient of factor, obtain the penalty coefficient F of each sub-aperture by the such as least square method shown in (6) ik;
min = Σ i = 1 ... N - 1 Σ j = 0 ... N - 1 j ∩ i ( ( ϵ asphere i ( x , y ) + Σ k = 1 L F i k f k ( x , y ) ) - ( ϵ asphere j ( x , y ) + Σ k = 1 L F j k f k ( x , y ) ) ) 2 - - - ( 6 )
Wherein: N is the number of sub-aperture;
Sub-aperture counts from zero, i.e. sub-aperture centered by the 0th sub-aperture;
representing only has the jth of an overlapping region sub-aperture to calculate to i-th sub-aperture;
2) according to step 1) in each sub-aperture diametric plane shape of obtaining splicing obtains unified shape, by used be added, the face shape of overlapping region gets the weighted mean value of each sub-aperture data.
8. the stitching interferometer of employing distorting lens according to claim 4 detects the device of aspheric surface, and it is characterized in that, described optical design software is Zemax software.
CN201510922163.XA 2015-12-14 2015-12-14 The device of aspheric surface is detected using the stitching interferometer of distorting lens Expired - Fee Related CN105423948B (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201510922163.XA CN105423948B (en) 2015-12-14 2015-12-14 The device of aspheric surface is detected using the stitching interferometer of distorting lens
PCT/CN2016/100862 WO2017101557A1 (en) 2015-12-14 2016-09-29 Surface shape detection device and method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510922163.XA CN105423948B (en) 2015-12-14 2015-12-14 The device of aspheric surface is detected using the stitching interferometer of distorting lens

Publications (2)

Publication Number Publication Date
CN105423948A true CN105423948A (en) 2016-03-23
CN105423948B CN105423948B (en) 2018-10-16

Family

ID=55502328

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510922163.XA Expired - Fee Related CN105423948B (en) 2015-12-14 2015-12-14 The device of aspheric surface is detected using the stitching interferometer of distorting lens

Country Status (2)

Country Link
CN (1) CN105423948B (en)
WO (1) WO2017101557A1 (en)

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105783780A (en) * 2016-04-29 2016-07-20 浙江大学 Free-form surface unconventional subaperture splicing interference detection device and method
CN106091978A (en) * 2016-06-01 2016-11-09 西安工程大学 The joining method of interference fringe image in inclined in type measurements by laser interferometry
WO2017101557A1 (en) * 2015-12-14 2017-06-22 中国科学院长春光学精密机械与物理研究所 Surface shape detection device and method
CN107782254A (en) * 2017-09-29 2018-03-09 中国科学院长春光学精密机械与物理研究所 A kind of mixed compensating mode sub-aperture stitching surface testing method
CN107796334A (en) * 2017-11-29 2018-03-13 许昌学院 A kind of surface testing system of aspherical optical element
WO2018068199A1 (en) * 2016-10-11 2018-04-19 中国科学院长春光学精密机械与物理研究所 Measurement apparatus and method for transmission wavefront equation of aspherical surface compensator
WO2018068225A1 (en) * 2016-10-12 2018-04-19 中国科学院长春光学精密机械与物理研究所 Measurement apparatus and measurement method for surface-shape error of rotating-axis symmetric curved surface
CN108267094A (en) * 2018-01-12 2018-07-10 暨南大学 A kind of non-cylindrical interference splicing measuring systems and method based on rotation CGH
CN109506589A (en) * 2018-12-25 2019-03-22 东南大学苏州医疗器械研究院 A kind of measuring three-dimensional profile method based on light field imaging
CN110109230A (en) * 2019-05-24 2019-08-09 西安交通大学 A kind of intelligent splicing assembly method of aspheric noodles complex curved surface part
CN110332883A (en) * 2019-07-22 2019-10-15 中国科学院上海光学精密机械研究所 The removing method of fizeau interferometer hysterisis error
CN110487212A (en) * 2019-08-02 2019-11-22 中北大学 One kind being based on vortex light helical phase phase shift interference detectable substance dignity type device
WO2020172897A1 (en) * 2019-02-26 2020-09-03 中国科学院光电技术研究所 Lens surface shape detection device in stitching interferometer
CN111623957A (en) * 2020-05-11 2020-09-04 中国科学院光电技术研究所 Point cloud registration and splicing method for X-ray focusing lens splicing interference detection
CN111811429A (en) * 2020-07-14 2020-10-23 北京理工大学 Sub-aperture splicing interference measurement method and device
CN111895934A (en) * 2020-07-31 2020-11-06 北京理工大学 Optical element surface local gradient surface shape error interferometry method and device
CN113029022A (en) * 2021-02-25 2021-06-25 中国人民解放军国防科技大学 Shape and position error interference measurement device and method for transparent hemispherical shell part
CN113091644A (en) * 2021-06-09 2021-07-09 中国工程物理研究院激光聚变研究中心 Large-aperture optical element surface shape detection method based on stacked coherent diffraction imaging
CN113091638A (en) * 2021-03-25 2021-07-09 中国科学院光电技术研究所 Device and method for calibrating return error in surface shape measurement by interferometry
CN116577931A (en) * 2023-07-14 2023-08-11 中国科学院长春光学精密机械与物理研究所 Optical element splicing detection method based on instrument transfer function
CN117249912A (en) * 2023-11-20 2023-12-19 苏州致将智能光电有限公司 Method and system for detecting large-caliber optical element

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110243306B (en) * 2019-07-22 2024-06-11 中国工程物理研究院激光聚变研究中心 Plane surface shape sub-aperture splicing interferometry device and method based on robot
CN110779443B (en) * 2019-11-04 2021-04-02 中国科学院国家天文台南京天文光学技术研究所 Edge sensor for splicing mirror surface based on interference principle and working method thereof
CN113465539B (en) * 2021-07-06 2024-03-19 上海大学 Automatic cylindricity measuring device and method based on sub-aperture interference splicing
CN115790442B (en) * 2022-11-15 2024-05-14 南京理工大学 Interferometry method based on large-caliber micro-displacement adjusting frame
CN118293787B (en) * 2024-06-05 2024-08-16 杭州晶耐科光电技术有限公司 Spherical dynamic interferometer based on four-wave common-path interference and detection method thereof

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10132531A (en) * 1996-11-05 1998-05-22 Canon Inc Measuring method for aspherical form and its device
JP2002267426A (en) * 2001-03-12 2002-09-18 Canon Inc Shape-measuring instrument and method
CN1746648A (en) * 2005-10-31 2006-03-15 中国科学院光电技术研究所 Heavy-calibre deep type aspherical mirror detecting system
CN101709955A (en) * 2009-11-24 2010-05-19 中国科学院长春光学精密机械与物理研究所 Device for detecting surface shape of optical aspheric surface by sub-aperture stitching interferometer
CN102506750A (en) * 2011-10-28 2012-06-20 中国科学院长春光学精密机械与物理研究所 Partial-compensation aspherical reflector surface shape detection method
CN103776389A (en) * 2014-01-10 2014-05-07 浙江大学 High-precision aspheric combined interference detection device and high-precision aspheric combined interference detection method

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4917088B2 (en) * 2005-04-05 2012-04-18 キューイーディー・テクノロジーズ・インターナショナル・インコーポレーテッド High resolution precision measurement method for non-spherical surfaces
CN100567932C (en) * 2008-05-28 2009-12-09 中国科学院光电技术研究所 Fan-shaped off-axis aspherical mirror splicing measurement system
JP6000577B2 (en) * 2012-03-09 2016-09-28 キヤノン株式会社 Aspherical surface measuring method, aspherical surface measuring device, optical element processing apparatus, and optical element manufacturing method
CN102661719B (en) * 2012-04-16 2014-03-26 中国人民解放军国防科学技术大学 Near-null compensator, surface shape measuring instrument and measuring method for matching measurement of sub-apertures of aspheric surfaces
CN105423948B (en) * 2015-12-14 2018-10-16 中国科学院长春光学精密机械与物理研究所 The device of aspheric surface is detected using the stitching interferometer of distorting lens

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10132531A (en) * 1996-11-05 1998-05-22 Canon Inc Measuring method for aspherical form and its device
JP2002267426A (en) * 2001-03-12 2002-09-18 Canon Inc Shape-measuring instrument and method
CN1746648A (en) * 2005-10-31 2006-03-15 中国科学院光电技术研究所 Heavy-calibre deep type aspherical mirror detecting system
CN101709955A (en) * 2009-11-24 2010-05-19 中国科学院长春光学精密机械与物理研究所 Device for detecting surface shape of optical aspheric surface by sub-aperture stitching interferometer
CN102506750A (en) * 2011-10-28 2012-06-20 中国科学院长春光学精密机械与物理研究所 Partial-compensation aspherical reflector surface shape detection method
CN103776389A (en) * 2014-01-10 2014-05-07 浙江大学 High-precision aspheric combined interference detection device and high-precision aspheric combined interference detection method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
《中国博士学位论文全文数据库 基础科学辑》 *
田超: "《非球面非零位环形子孔径拼接干涉检测技术与系统研究》", 《中国博士学位论文全文数据库》 *

Cited By (31)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017101557A1 (en) * 2015-12-14 2017-06-22 中国科学院长春光学精密机械与物理研究所 Surface shape detection device and method
CN105783780A (en) * 2016-04-29 2016-07-20 浙江大学 Free-form surface unconventional subaperture splicing interference detection device and method
CN106091978A (en) * 2016-06-01 2016-11-09 西安工程大学 The joining method of interference fringe image in inclined in type measurements by laser interferometry
WO2018068199A1 (en) * 2016-10-11 2018-04-19 中国科学院长春光学精密机械与物理研究所 Measurement apparatus and method for transmission wavefront equation of aspherical surface compensator
WO2018068225A1 (en) * 2016-10-12 2018-04-19 中国科学院长春光学精密机械与物理研究所 Measurement apparatus and measurement method for surface-shape error of rotating-axis symmetric curved surface
CN107782254A (en) * 2017-09-29 2018-03-09 中国科学院长春光学精密机械与物理研究所 A kind of mixed compensating mode sub-aperture stitching surface testing method
CN107782254B (en) * 2017-09-29 2019-05-21 中国科学院长春光学精密机械与物理研究所 A kind of mixed compensating mode sub-aperture stitching surface testing method
CN107796334A (en) * 2017-11-29 2018-03-13 许昌学院 A kind of surface testing system of aspherical optical element
CN108267094A (en) * 2018-01-12 2018-07-10 暨南大学 A kind of non-cylindrical interference splicing measuring systems and method based on rotation CGH
CN109506589B (en) * 2018-12-25 2020-07-28 东南大学苏州医疗器械研究院 Three-dimensional profile measuring method based on structural light field imaging
CN109506589A (en) * 2018-12-25 2019-03-22 东南大学苏州医疗器械研究院 A kind of measuring three-dimensional profile method based on light field imaging
WO2020172897A1 (en) * 2019-02-26 2020-09-03 中国科学院光电技术研究所 Lens surface shape detection device in stitching interferometer
US11268808B2 (en) 2019-02-26 2022-03-08 The Institute Of Optics And Electronics, The Chinese Academy Of Sciences Detection device for detecting lens surface in stitching interferometer
CN110109230A (en) * 2019-05-24 2019-08-09 西安交通大学 A kind of intelligent splicing assembly method of aspheric noodles complex curved surface part
CN110332883A (en) * 2019-07-22 2019-10-15 中国科学院上海光学精密机械研究所 The removing method of fizeau interferometer hysterisis error
CN110487212B (en) * 2019-08-02 2021-04-16 中北大学 Device for detecting object surface shape based on vortex optical spiral phase shift interference
CN110487212A (en) * 2019-08-02 2019-11-22 中北大学 One kind being based on vortex light helical phase phase shift interference detectable substance dignity type device
CN111623957A (en) * 2020-05-11 2020-09-04 中国科学院光电技术研究所 Point cloud registration and splicing method for X-ray focusing lens splicing interference detection
CN111811429A (en) * 2020-07-14 2020-10-23 北京理工大学 Sub-aperture splicing interference measurement method and device
CN111811429B (en) * 2020-07-14 2021-04-20 北京理工大学 Sub-aperture splicing interference measurement method and device
CN111895934B (en) * 2020-07-31 2021-07-16 北京理工大学 Optical element surface local gradient surface shape error interferometry method and device
CN111895934A (en) * 2020-07-31 2020-11-06 北京理工大学 Optical element surface local gradient surface shape error interferometry method and device
CN113029022A (en) * 2021-02-25 2021-06-25 中国人民解放军国防科技大学 Shape and position error interference measurement device and method for transparent hemispherical shell part
CN113091638A (en) * 2021-03-25 2021-07-09 中国科学院光电技术研究所 Device and method for calibrating return error in surface shape measurement by interferometry
CN113091638B (en) * 2021-03-25 2022-06-28 中国科学院光电技术研究所 Device and method for calibrating return error in surface shape measurement by interferometry
CN113091644B (en) * 2021-06-09 2021-09-07 中国工程物理研究院激光聚变研究中心 Large-aperture optical element surface shape detection method based on stacked coherent diffraction imaging
CN113091644A (en) * 2021-06-09 2021-07-09 中国工程物理研究院激光聚变研究中心 Large-aperture optical element surface shape detection method based on stacked coherent diffraction imaging
CN116577931A (en) * 2023-07-14 2023-08-11 中国科学院长春光学精密机械与物理研究所 Optical element splicing detection method based on instrument transfer function
CN116577931B (en) * 2023-07-14 2023-09-22 中国科学院长春光学精密机械与物理研究所 Optical element splicing detection method based on instrument transfer function
CN117249912A (en) * 2023-11-20 2023-12-19 苏州致将智能光电有限公司 Method and system for detecting large-caliber optical element
CN117249912B (en) * 2023-11-20 2024-02-13 苏州致将智能光电有限公司 Method and system for detecting large-caliber optical element

Also Published As

Publication number Publication date
WO2017101557A1 (en) 2017-06-22
CN105423948B (en) 2018-10-16

Similar Documents

Publication Publication Date Title
CN105423948A (en) Splicing-interference-detection aspheric surface shape apparatus using distorting lens and method thereof
CN102997863B (en) Direct detection system for surface-shape errors in full-aperture optical aspheric surfaces
CN110188321B (en) Primary and secondary mirror calibration method based on neural network algorithm
CN103591888B (en) The measuring method of large-caliber off-axis non-spherical optical element geometric parameter
CN103134660B (en) Method for obtaining alignment error of primary and secondary telescope mirrors based on astigmatic decomposition
CN103335610B (en) Detection system for large-caliber high-order convex aspheric surface
CN101539397B (en) Method for measuring three-dimensional attitude of object on precision-optical basis
CN101290218B (en) Method for correcting principle error of aspherical non-zero digit detection
CN107421436B (en) Aspherical interferometer measuration system and method based on the spatial light modulator plane of reference
JP5896792B2 (en) Aspherical surface measuring method, aspherical surface measuring device, and optical element processing device
CN102620683B (en) Sub-aperture stitching detects the compensation method of aspheric surface alignment error
CN102288132B (en) Method for measuring vertex curvature radius deviation of aspheric surface by using laser tracking instrument
CN104374334B (en) Free form surface pattern method for three-dimensional measurement and device
CN103353388A (en) Method and device for calibrating binocular integrated microscopy imaging system with camera shooting function
CN106574834A (en) Method for calibrating a measuring device
CN105318847A (en) Aspheric non-zero digit circular subaperture stitching method based on system modeling
CN103234480A (en) Rapid surface shape detection method for circular convex aspheric surfaces
CN106404354A (en) Device and method for measurement of aspheric compensator transmission wavefront equation
CN102937421A (en) Real-time detection method of symmetrical optical non-spherical face of rotary shaft
CN106705888B (en) CCD coordinate system and mirror surface coordinate system non-linear relation scaling method in interference detection
CN105352451B (en) A kind of accurate omnipotent compensating glass and design method based on deformable mirror
CN109002567B (en) Computer-aided debugging method based on actual optical system debugging detection optical path
CN110108271B (en) Star deflection compensation method caused by pneumatic optical effect
CN102927930B (en) Method for detecting ultra-large-diameter reflector surface errors in splicing mode by adopting collimator
CN108267094B (en) Non-cylindrical surface interference splicing measurement system and method based on rotary CGH

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
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: 20181016

Termination date: 20201214