CN104613888B - Smog deformation of body measuring method is passed through under a kind of flame smoky environment - Google Patents

Smog deformation of body measuring method is passed through under a kind of flame smoky environment Download PDF

Info

Publication number
CN104613888B
CN104613888B CN201510056306.3A CN201510056306A CN104613888B CN 104613888 B CN104613888 B CN 104613888B CN 201510056306 A CN201510056306 A CN 201510056306A CN 104613888 B CN104613888 B CN 104613888B
Authority
CN
China
Prior art keywords
speckle
deformation
flame
smog
strain
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
CN201510056306.3A
Other languages
Chinese (zh)
Other versions
CN104613888A (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.)
Tsinghua University
Original Assignee
Tsinghua University
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 Tsinghua University filed Critical Tsinghua University
Priority to CN201510056306.3A priority Critical patent/CN104613888B/en
Publication of CN104613888A publication Critical patent/CN104613888A/en
Application granted granted Critical
Publication of CN104613888B publication Critical patent/CN104613888B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Length Measuring Devices By Optical Means (AREA)

Abstract

Through the deformation of body measuring method of smog under a kind of flame smoky environment, belong to malformation measurement and optical measurement mechanics field of measuring technique.The present invention irradiates surface and prepares the testee surface for having high temperature reflective mirror speckle by LASER Light Source or narrow-band light source, with CCD camera of the camera lens with optical filter object under test surface is shot through flame smog, obtain body surface minute surface speckle blurred picture, the clear original image that testee surface speckle is not obscured by flame smog is recovered by Phase Retrieve Algorithm, body surface is calculated by the body surface speckle shape profile variations during testee temperature distortion to deform, realize real-time measurement of full field testee areal deformation strain and displacement.The present invention is solved through flame smog to the real-time measurement of full field of deformation of body, be can be used under flame smoky environment object in stress or is heated etc. and to cause to carry out real-time full field deformation measure in deformation even destructive process.

Description

Smog deformation of body measuring method is passed through under a kind of flame smoky environment
Technical field
The present invention relates to pass through smog deformation of body measuring method under a kind of flame smoky environment, belong to malformation measurement With optical measurement mechanics field of measuring technique.
Background technology
It is the problem that a people pay close attention to research for a long time through scattering medium blur ed image restoration, by some image procossings Method, be attempted to by object by scattering medium into blur ed image restoration through scattering medium the clear preimage of object, This has been made that many research work in astronomical observation and medically through bio-tissue imaging research field.From the upper world seven Since the ten's, recovering blurred picture by Phase Retrieve Algorithm turns into a study hotspot, seeks the laser sent from object Or the method that narrow band light goes out the clear artwork of object through the blur ed image restoration that nonisotropic scattering medium is obtained.Subsequent more than 30 years In, people propose many methods and respective algorithms for different situations, for example relax average alternating reflex method (Luke, And phase controlling HIO methods (Harder et al. .2010) etc. 2005).
Optical measurement mechanics is the important method of deformation measurement, has the advantages that noncontact whole-field measuring.Modern optical measurement mechanics is borrowed The development of digital image processing techniques and computer technology is helped, has developed a set of Digital Image Correlation Method measurement test specimen or structure Deform under load.Digital Image Correlation Method makes random speckle image by test specimen surface, by digital correlation Algorithm calculates the displacement of each speckle subregion on the front and rear test specimen of deformation, so as to obtain the overall strain of test specimen.By Simple in party's subtraction unit, simple operation obtains extensive use in engineering structure part deformation measurement.Rim detection is directed to number The computer digital image algorithm at word image detection testee edge, object structures part deformation in can to body surface some Mark detects its edge, by the deformation for marking the deformation of profile to embody testee.However, Digital Image Correlation Method In it is very sensitive to the change of the brightness fluctuation of light source lighting speckle, it is and due to environmental disturbances in some test environments, such as high Warm heating environment Flame smog wraps up test specimen, and camera shoots the surface of test piece speckle image for obtaining and is distorted even It is smudgy, cannot directly be deformed with digital image methods using blurred picture and calculated.
The content of the invention
The purpose of the present invention is to propose to pass through smog deformation of body measuring method under a kind of flame smoky environment, with phase Recovery algorithms go out the clear artwork of object from object through the blur ed image restoration that nonisotropic scattering medium is obtained, with reference to digital picture Edge detection method, deformation of body is calculated using testee surface speckle profile variations.
Technical scheme is as follows:Flame smog deformation of body measuring method is passed through under a kind of hot environment, including Following steps:
1) reflective high temperature resistant minute surface speckle is made on testee surface;
2) testee of the surface with speckle is put into the chamber filled with flame smog, is placed in the hot loading bench of power On, object under test is wrapped up by flame and smog, with light source lighting testee surface, clapped with the CCD camera with optical filter The body surface minute surface speckle wrapped up by flame smog is taken the photograph, the body surface blurred picture of minute surface speckle reflection is obtained, will be shot The incoming image computing module of blurred picture for obtaining;
3) image computing module uses Phase Retrieve Algorithm, and it is not fiery that body surface blur ed image restoration is gone out into testee The fuzzy clear initial surface minute surface speckle image of flame smog, Phase Retrieve Algorithm uses iterative method, iteration until final acquisition The clear initial surface minute surface speckle image that testee is not obscured by flame smog;
4) the hot loading bench of power is loaded to testee, and body surface is obtained with edge detection algorithm respectively before and after loading Minute surface speckle edge contour, is shaped as reference pattern, by minute surface speckle wheel before taking loading for each minute surface speckle edge contour Wide change calculations go out body surface minute surface speckle minute surface speckle deformation gradient F preceding relative to deformation after deformation, further obtain Green Strain tensor E:
F in formulaTIt is the transposition of deformation gradient F, I is unit tensor, and object is obtained by the strain interpolation of each minute surface speckle The strain of the whole audience, in the case of small deformation, small strain strain tensor is that Cauchy strain tensor ε is equal to above-mentioned Green strain tensor E, Tri- component ε=(ε of ε ≈ E, whole audience strain field εxyxy), wherein x directions normal strain εx, y directions normal strain εyShould with shearing Become γxyDefinition be respectively:
Wherein u, v be followed successively by body surface in the x-direction, the displacement of y directions;Thing is obtained by body surface whole audience strain field integration Body surface facial disfigurement displacement field (u, v):
Wherein C is path of integration, uO、vOPath of integration starting point x directions, the displacement constant in y directions are followed successively by,C1It is integral constant,
Reflective high temperature resistant minute surface speckle of the present invention is sprayed at testee surface, minute surface using aluminum oxide or silver Being shaped as speckle is square, rhombus or triangle.The transmission peak wavelength of CCD camera optical filter of the present invention is in light source used Cardiac wave is long.
The present invention compared with prior art, with advantages below and high-lighting effect:By digital image processing method application Measured in deformation of body under flame smoky environment, overcome Fuzzy Influence of the flame smog to image, solve extremely complicated ring The body surface deformation measurement wrapped up by flame smog under border.
Brief description of the drawings
Fig. 1 is shown through the structure of smog deformation of body measuring method embodiment under the flame smoky environment that the present invention is provided It is intended to.
Fig. 2 is the body surface speckle shape example used by the present invention.
Fig. 3 is that testee surface speckle profile variations of the present invention survey body surface deformation principle schematic diagram.
In accompanying drawing:1- chambers;2- watch windows;3- flame fog guns;The hot loading bench of 4- power;5- light sources;6-CCD phases Machine;7- image computing modules;8- testees.
Specific embodiment
Specific embodiment of the invention is further illustrated below in conjunction with the accompanying drawings, but protection of the invention should not be limited with this Scope.
Through the structural representation of smog deformation of body measuring method embodiment under the flame smoky environment that Fig. 1 present invention is provided Figure.Testee surface aluminum oxide or silvery make can reflection light minute surface speckle, speckle shape is square, rhombus or triangle Shape, speckle surface it is smooth present minute surface can reflected illumination to its surface light, the body surface used by the present embodiment as shown in Figure 2 Special shape speckle example.
Testee of the surface with minute surface speckle is put into the chamber filled with flame smog, power heat loading is positioned over On platform, object under test is wrapped up by flame and smog, and testee surface is illuminated with LASER Light Source, and testee is received on loading bench Power or temperature distortion, the body surface wrapped up by flame smog is shot with the CCD camera with optical filter, obtains minute surface speckle The body surface blurred picture of reflection, the incoming image computing module of blurred picture that shooting is obtained, image computing module uses phase Bit recovery algorithm recovers the clear initial surface speckle image that testee is not obscured by flame smog, uses edge detection algorithm Body surface speckle edge contour is obtained, body surface speckle profile is obtained respectively before and after being loaded to testee, for every Individual speckle is shaped as reference pattern before taking deformation, and body surface whole audience strain field and position after deforming are calculated by speckle profile variations Move field.
Operation principle of the invention is as follows:
Shot through non-uniform dielectric and obtain the picture rich in detail that blur ed image restoration is not blurred, calculated using phase recovery Method, its algorithm is:
Wherein operatorFourier transformation operator and inverse Fourier transform operator are followed successively by, f (x, y) is survey The clear original image gray scale of body surface speckle is taken the photograph, (x, y) is the pixel point coordinates of gray level image, and F (u, v) is f's (x, y) Fourier transformation valueI is imaginary unitK is iterative steps, φk(u, v) is gk(x, Y) phase of Fourier transformation value, θ 'k(u, v) is G'kThe phase of (u, v) inverse Fourier transform value, γ is decision space on (x, y), For example (x, y) in γ (x, y) ∈ γ g'k(x, y) value is negative, and (x, y) be not in γG'k(x, y) value It is nonnegative number, iterative initial value g1(x, y) is random distribution gray-scale map, and the variance of calculation error is:
As calculation error variance EkGo to zero, iterative calculation terminates, it is believed that the g for calculatingk+1(x, y) is and is not blurred The clear original image of object.It is f'(x that camera directly shoots the original image that testee surface obtains, y), can by following formula To obtain the mould of Fourier transformation value F (u, v) of f (x, y) | F (u, v) |:
W (x, y) is calculation window function in above formula, and typically using rectangular window, window size recovers speckle with body surface Area size is roughly the same.I.e. according to recovery region filter function is calculated, the value outside window is taken as zero for its effect.
On the basis of above-mentioned Phase Retrieve Algorithm formula (I), some iteration qualifications can be increased and accelerate convergence rate, Different calculation windows for example is used in different iterative calculation steps, is recovered image size and is fallen in calculation window, iterated to calculate Only within calculation window, iterated conditional can also be exported using mixed iteration, i.e., the g in γ domainsk+1(x, y)=gk(x,y)-β g'k(x, y), wherein β are the real number between 0 to 1.According to actual conditions, can be changed using different in different iteration steps For condition, increase convergence rate.
By the body surface image with speckle to recovering, go out surface speckle profile using rim detection, by thing Deform caused surface speckle profile variations in body thermal histories to be calculated the displacement field of body surface deformation, by position Move field further calculate testee temperature distortion in real time, the strain field of the whole audience.
The speckle picture rich in detail that testee surface is not obscured by smog is restored with phase recovery method, in object stand under load In deformation process, speckle profile also changes.Speckle profile is obtained by rim detection, is calculated by speckle profile variations and obtained Speckle deformation gradient F, as shown in figure 3, can directly calculate strain field by deformation gradient, it is to avoid by displacement field once differentiation The error for calculating strain field is amplified.Real deformation process is original minute surface speckle shape D0It is changed into shape D through deformation gradient F, After trying to achieve deformation gradient F by minute surface speckle profile variations, for Finite-Deformation Problems, typically directly Green strain is obtained using F Tensor E:
F in formulaTIt is the transposition of deformation gradient F, I is unit tensor.Object is obtained by the strain interpolation of each minute surface speckle The strain of the surface whole audience, in the case of small deformation, conventional small strain strain tensor is that Cauchy strain tensor ε is equal to above-mentioned Green Tri- component ε=(ε of strain tensor E, ε ≈ E, whole audience strain field εxyxy), wherein x directions normal strain εx, y directions normal strain εyWith shear strain γxyDefinition be respectively:
Wherein u, v be followed successively by body surface in the x-direction, the displacement of y directions.Thing is obtained by body surface whole audience strain field integration Body surface facial disfigurement displacement field (u, v):
Wherein C is path of integration, uO、vOPath of integration starting point x directions, the displacement constant in y directions are followed successively by,C1It is integral constant,

Claims (5)

1. smog deformation of body measuring method is passed through under a kind of flame smoky environment, it is characterised in that the method includes following step Suddenly:
1) reflective high temperature resistant minute surface speckle is made on testee (8) surface;
2) testee (8) of the surface with speckle is put into the chamber filled with flame smog (1), is placed in power heat loading On platform (4), object under test is wrapped up by flame and smog, testee surface is illuminated with light source (5), with optical filter CCD camera (6) shoots the body surface minute surface speckle wrapped up by flame smog, and the body surface for obtaining minute surface speckle reflection is obscured Image, will shoot the incoming image computing module (7) of blurred picture for obtaining;
3) image computing module (7) uses Phase Retrieve Algorithm, and it is not fiery that body surface blur ed image restoration is gone out into testee The fuzzy clear initial surface minute surface speckle image of flame smog, Phase Retrieve Algorithm uses iterative method, iteration until final acquisition The clear initial surface minute surface speckle image that testee is not obscured by flame smog;
4) testee (8) is loaded using the hot loading bench of power (4), thing is obtained with edge detection algorithm respectively before and after loading Body surface face minute surface speckle edge contour, is shaped as reference pattern, by minute surface before taking loading for each minute surface speckle edge contour Speckle profile variations calculate body surface minute surface speckle minute surface speckle deformation gradient F preceding relative to deformation after deformation, further Obtain Green strain tensor E:
E = 1 2 ( F T F - I )
F in formulaTIt is the transposition of deformation gradient F, I is unit tensor, and the object whole audience is obtained by the strain interpolation of each minute surface speckle Strain, in the case of small deformation, small strain strain tensor is that Cauchy strain tensor ε is equal to above-mentioned Green strain tensor E, ε=E, Body surface whole audience strain field Cauchy strain tensor ε has three component ε=(εxyxy), wherein x directions normal strain εx, y side To normal strain εyWith shear strain γxyDefinition be respectively:
ϵ x = ∂ u ∂ x
ϵ y = ∂ v ∂ y
γ x y = ∂ u ∂ y + ∂ v ∂ x
Wherein u, v be followed successively by body surface in the x-direction, the displacement of y directions;Object table is obtained by body surface whole audience strain field integration Facial disfigurement displacement field (u, v):
u = ∫ C ( ϵ x d x + ∂ u ∂ y d y ) + u O
v = ∫ C ( ϵ y d y + ∂ v ∂ x d x ) + v O
Wherein C is path of integration, uO、vOPath of integration starting point x directions, the displacement constant in y directions are followed successively by,C1It is integral constant,
2., according to smog deformation of body measuring method is passed through under a kind of flame smoky environment described in claim 1, its feature exists Testee surface is sprayed at using aluminum oxide or silver in the reflective high temperature resistant minute surface speckle.
3., according to smog deformation of body measuring method is passed through under a kind of flame smoky environment described in claim 1, its feature exists Square, rhombus or triangle are shaped as in the reflective high temperature resistant minute surface speckle.
4., according to smog deformation of body measuring method is passed through under a kind of flame smoky environment described in claim 1, its feature exists LASER Light Source or narrow-band light source are used in the light source.
5., according to smog deformation of body measuring method is passed through under a kind of flame smoky environment described in claim 1, its feature exists In:In step 2) described in the transmission peak wavelength of CCD camera optical filter be light source center wavelength used.
CN201510056306.3A 2015-02-03 2015-02-03 Smog deformation of body measuring method is passed through under a kind of flame smoky environment Active CN104613888B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510056306.3A CN104613888B (en) 2015-02-03 2015-02-03 Smog deformation of body measuring method is passed through under a kind of flame smoky environment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510056306.3A CN104613888B (en) 2015-02-03 2015-02-03 Smog deformation of body measuring method is passed through under a kind of flame smoky environment

Publications (2)

Publication Number Publication Date
CN104613888A CN104613888A (en) 2015-05-13
CN104613888B true CN104613888B (en) 2017-06-13

Family

ID=53148449

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510056306.3A Active CN104613888B (en) 2015-02-03 2015-02-03 Smog deformation of body measuring method is passed through under a kind of flame smoky environment

Country Status (1)

Country Link
CN (1) CN104613888B (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108414529A (en) * 2018-03-09 2018-08-17 天津商业大学 A kind of detection method and its system of rubber product air entrapment
CN109030546B (en) * 2018-07-23 2019-09-20 清华大学 High temperature deformation and temperature measurement system and method
CN109945776B (en) * 2019-02-01 2021-05-25 天津大学 Displacement/strain measurement method based on fluorescent particle marking and active optical measurement in electrochemistry
CN113808029A (en) * 2021-05-25 2021-12-17 南京航空航天大学 Strain smoothing method in digital image correlation
CN114018177B (en) * 2021-11-18 2023-12-26 南昌航空大学 Three-dimensional measurement method for mirror surface object based on speckle pattern
CN115616027B (en) * 2022-10-13 2023-10-13 江苏瑞亿扬材料科技有限公司 PVC calendering membrane temperature resistance detection device and detecting system

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE4338321C1 (en) * 1993-11-10 1995-01-12 Fraunhofer Ges Forschung Method for detecting a surface deformation by means of interferometry
CN101603819A (en) * 2009-07-23 2009-12-16 中国林业科学研究院木材工业研究所 A kind of real-time measurement method for wood deformation microstructure characteristics
CN102374849A (en) * 2011-09-22 2012-03-14 清华大学 High-temperature aerobic loading optical measurement system
CN102788558A (en) * 2012-07-30 2012-11-21 山东师范大学 Three-dimensional deformation measuring system and three-dimensional deformation measuring method combining speckle correlation and speckle interference
CN103363919A (en) * 2013-07-15 2013-10-23 国家电网公司 Device and method for measuring strain of high-temperature pipeline of power plant based on white light speckle technology
CN103558243A (en) * 2013-11-19 2014-02-05 北京航空航天大学 Optical method based high-speed aircraft hot surface full-field deformation measuring device
CN103559686A (en) * 2013-10-16 2014-02-05 中国科学院深圳先进技术研究院 Coaxial phase contrast imaging phase retrieval method and system based on multi-aspect image information
CN203534533U (en) * 2013-10-28 2014-04-09 卓力特光电仪器(苏州)有限公司 One-dimensional in-plane phase shift electronic speckle interferometer
KR101445239B1 (en) * 2013-03-25 2014-09-29 국방과학연구소 Apparatus and method for the noncontact deformation measurement at high temperature

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5637738B2 (en) * 2010-06-04 2014-12-10 キヤノン株式会社 Deformation measuring apparatus and deformation measuring method
JP2013170830A (en) * 2012-02-17 2013-09-02 Mitsubishi Heavy Ind Ltd Strain measuring device

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE4338321C1 (en) * 1993-11-10 1995-01-12 Fraunhofer Ges Forschung Method for detecting a surface deformation by means of interferometry
CN101603819A (en) * 2009-07-23 2009-12-16 中国林业科学研究院木材工业研究所 A kind of real-time measurement method for wood deformation microstructure characteristics
CN102374849A (en) * 2011-09-22 2012-03-14 清华大学 High-temperature aerobic loading optical measurement system
CN102788558A (en) * 2012-07-30 2012-11-21 山东师范大学 Three-dimensional deformation measuring system and three-dimensional deformation measuring method combining speckle correlation and speckle interference
KR101445239B1 (en) * 2013-03-25 2014-09-29 국방과학연구소 Apparatus and method for the noncontact deformation measurement at high temperature
CN103363919A (en) * 2013-07-15 2013-10-23 国家电网公司 Device and method for measuring strain of high-temperature pipeline of power plant based on white light speckle technology
CN103559686A (en) * 2013-10-16 2014-02-05 中国科学院深圳先进技术研究院 Coaxial phase contrast imaging phase retrieval method and system based on multi-aspect image information
CN203534533U (en) * 2013-10-28 2014-04-09 卓力特光电仪器(苏州)有限公司 One-dimensional in-plane phase shift electronic speckle interferometer
CN103558243A (en) * 2013-11-19 2014-02-05 北京航空航天大学 Optical method based high-speed aircraft hot surface full-field deformation measuring device

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
1200℃高温热环境下全场变形的非接触光学测量方法研究;潘兵_等;《强度与环境 》;20110215;第38卷(第1期);第52-59页 *
SUSAN角点检测和匹配算法在高温变形测量中的应用;于合龙_等;《航空学报》;20130525(第5期);第1064-1071页 *
Two-dimensional digital image correlation for in-plane displacement and strain measurement:A Review;Bing Pan_et al;《Measurement science and technology》;20090427;第20卷(第6期);第1-17页 *
基于数字图像相关方法的非接触高温热变形测量***;潘兵_等;《航空学报 》;20101025;第31卷(第10期);第1960-1967页 *
构造变形力学中有限应变张量几何解析;韩玉英;《中国科学(B辑 化学 生命科学 地学) 》;19910531(第5期);第519-532页 *

Also Published As

Publication number Publication date
CN104613888A (en) 2015-05-13

Similar Documents

Publication Publication Date Title
CN104613888B (en) Smog deformation of body measuring method is passed through under a kind of flame smoky environment
CN109839072B (en) DIC-based temperature field and deformation field synchronous measurement method and device
Yang et al. Shadow removal using bilateral filtering
JP5136965B2 (en) Image processing apparatus, image processing method, and image processing program
CN106949836B (en) Device and method for calibrating same-side target position of stereoscopic camera
CN104457603B (en) Object deformation measurement method under high-temperature environment
CN108572181A (en) A kind of mobile phone bend glass defect inspection method based on streak reflex
CN104535412A (en) Ultraviolet illumination DIC (digital image correction)-based mechanical property loading measuring system and method for high-temperature material
CN110910456B (en) Three-dimensional camera dynamic calibration method based on Harris angular point mutual information matching
CN113160339A (en) Projector calibration method based on Samm's law
Cofaru et al. A three-frame digital image correlation (DIC) method for the measurement of small displacements and strains
Nam et al. Multispectral photometric stereo for acquiring high-fidelity surface normals
CN108550160A (en) Non-homogeneous striation characteristic area extracting method based on light intensity template
Ciortan et al. A practical reflectance transformation imaging pipeline for surface characterization in cultural heritage
Wang et al. Improvement method of high-temperature digital image correlation measurement accuracy based on image processing
CN108256166A (en) A kind of data processing method for thermo-mapping technique
Pei et al. Profile measurement of non-Lambertian surfaces by integrating fringe projection profilometry with near-field photometric stereo
Suresh et al. PMENet: phase map enhancement for Fourier transform profilometry using deep learning
CN110910457B (en) Multispectral three-dimensional camera external parameter calculation method based on angular point characteristics
Crawford et al. Quantitative boundary-layer transition measurements using IR thermography
CN115655151B (en) Mobile phone rear cover plate detection device and method based on color phase measurement deflection technology
Yuan et al. Temperature measurement of metal surface at normal temperatures by visible images and machine learning
Liu et al. Neutralizing the impact of heat haze on digital image correlation measurements via deep learning
Stolz et al. Real time polarization imaging of weld pool surface
Liu et al. In-situ 3D shape measurement system for a high-temperature object in a plasma wind tunnel via DMD-projection laser structured light and adaptive speckle filtering

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