CN107024535A - A kind of multiple index depth detection method of the vertical defect based on surface wave - Google Patents

A kind of multiple index depth detection method of the vertical defect based on surface wave Download PDF

Info

Publication number
CN107024535A
CN107024535A CN201610195694.8A CN201610195694A CN107024535A CN 107024535 A CN107024535 A CN 107024535A CN 201610195694 A CN201610195694 A CN 201610195694A CN 107024535 A CN107024535 A CN 107024535A
Authority
CN
China
Prior art keywords
defect
depth
wave
coefficient
echo
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
CN201610195694.8A
Other languages
Chinese (zh)
Other versions
CN107024535B (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.)
Beijing University of Technology
Original Assignee
Beijing University of Technology
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 Beijing University of Technology filed Critical Beijing University of Technology
Priority to CN201610195694.8A priority Critical patent/CN107024535B/en
Publication of CN107024535A publication Critical patent/CN107024535A/en
Application granted granted Critical
Publication of CN107024535B publication Critical patent/CN107024535B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/04Analysing solids
    • G01N29/041Analysing solids on the surface of the material, e.g. using Lamb, Rayleigh or shear waves
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00Indexing codes associated with group G01N29/00
    • G01N2291/02Indexing codes associated with the analysed material
    • G01N2291/023Solids
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00Indexing codes associated with group G01N29/00
    • G01N2291/04Wave modes and trajectories
    • G01N2291/042Wave modes
    • G01N2291/0423Surface waves, e.g. Rayleigh waves, Love waves
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00Indexing codes associated with group G01N29/00
    • G01N2291/04Wave modes and trajectories
    • G01N2291/042Wave modes
    • G01N2291/0426Bulk waves, e.g. quartz crystal microbalance, torsional waves
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00Indexing codes associated with group G01N29/00
    • G01N2291/04Wave modes and trajectories
    • G01N2291/042Wave modes
    • G01N2291/0428Mode conversion

Landscapes

  • Physics & Mathematics (AREA)
  • Acoustics & Sound (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

A kind of multiple index depth detection method of the vertical defect based on surface wave, belongs to supersonic guide-wave Non-Destructive Testing and evaluation areas.Surface wave is in communication process, the bulk wave of scattering is had when being interacted with defect below defect, the lower surface that bodily wave propagation runs into workpiece can be reflected, surface wave such as (namely mode conversion echo) can be converted to once again and be propagated along defect two ends by running into defect, and it is with the total distance that the time difference of flaw echo is that shear wave is propagated in thickness of workpiece direction.The mode conversion echo at defect two ends with flaw echo and transmitted wave synthesis can as depth of defect characterization parameter, the depth detection method of the vertical defect of multiple index is solved brings inaccuracy and unworthiness for the detection more than 0.45 times of wavelength depth of defect at present.In health monitoring and Nondestructive Evaluation field, with great application value and potentiality.

Description

A kind of multiple index depth detection method of the vertical defect based on surface wave
Technical field
The present invention relates to a kind of multiple index depth detection method of the vertical defect based on surface wave, belong to supersonic guide-wave without Damage detection and evaluation areas.
Background technology
Using surface wave to carry out, Surface Flaw Detection is especially vertical to surface or the vertical defect of class is carried out at present When depth detection, according to parameter informations such as the size of workpiece, materials, pass through pop one's head in transmitting specific frequency or specific wavelength Surface wave carries out defects detection.Excitating surface ripple and flaw echo are calculated according to resulting signal waveform (self excitation and self receiving) Time difference Δ t, is multiplied by time difference Δ t with the known surface wave velocity of wave v in this workpiece and draws the total distance that surface wave is propagated, The position of defect is determined with the half of total distance.For the detection of depth of defect, often a form received is swashed using one.Enter first The foundation of row emulation coefficient curve, the signal waveform received according to receiving transducer records and calculates the amplitude of direct wave, defect The amplitude of echo and the amplitude of transmitted wave, on the basis of through wave amplitude, with the amplitude and the width of transmitted wave of flaw echo The coefficient for obtaining characterizing depth of defect information is normalized in value divided by through wave amplitude.Work is determined further according to related experiment The characterization parameter of part depth of defect, the detection that coefficient of correlation curve carries out depth of defect is determined.But, depth of defect characterizes coefficient Simply show preferable monotonicity in the range of depth of defect is 0.45 times of wavelength, depth of defect be 0.45 times of wavelength with The opposite shape condition of monotonicity is showed in scope afterwards, that is to say, that whole piece coefficient curve is in general 0.45 times of wavelength depth of defect There is flex point at place.This just brings inaccuracy and unworthiness for the detection more than 0.45 times of wavelength depth of defect.
The content of the invention
For above-mentioned existing problem, the multiple index that the present invention proposes a kind of vertical defect based on surface wave is deep Detection method is spent, solves and brings inaccuracy and inapplicable for the detection more than 0.45 times of wavelength depth of defect at present Property.Surface wave has the bulk wave of scattering in communication process below defect when being interacted with defect, bodily wave propagation runs into work The lower surface of part can be reflected, and surface wave such as accompanying drawing 1-2 (namely mode conversion echo) can be converted to once again by running into defect And propagated along defect two ends, it is with the total distance that the time difference of flaw echo is that shear wave is propagated in thickness of workpiece direction.Defect The mode conversion echo at two ends with flaw echo and transmitted wave synthesis can as depth of defect characterization parameter, multiple index The depth detection method of vertical defect solve brought at present for the detection more than 0.45 times of wavelength depth of defect it is inaccurate Property and unworthiness.
Step 1:Establish emulation polyphyly number curve
The parameter informations such as size and material for target workpiece, it is determined that suitable simulation software, sets up suitable emulation Model.According to required, depth of defect size range is determined, serial emulation is carried out.According to simulation model, extract at observation station or face Shape information, with resulting flaw echo value, transmission wave number and defect two ends at mode conversion echo (such as accompanying drawing 1) As ordinate after value divided by the through wave number normalized of benchmark, then with depth of defect or depth of defect divided by wavelength normalization Processing sets up multiple index depth of defect and characterizes curve (such as accompanying drawing 2) as abscissa.
Step 2:Carry out experiment detection
The stimulating frequency set up according to the size of target workpiece and material parameter information and emulation polyphyly number curve and swash Originating party formula, the suitable probe of selection is tested.As shown in figure 3, experimental system includes test piece, a sharp receipts probe, shown Ripple device, excitation set;One swashs a receipts probe is connected with test piece, and one sharp one, which receives probe, passes through connecting line and excitation set company Connect, excitation set is connected with oscillograph.The mode that experiment measures at the defect clawback, transmitted wave and defect two ends of test piece turns Ripple (such as accompanying drawing 1) is gained, record obtains the waveforms amplitude of four sign depth of defects of the above.
Step 3:Determine depth of defect
Four defect informations of mode conversion echo at defect clawback, transmitted wave and defect two ends measured according to experiment, Divided by a reference value direct wave obtains one group of four coefficient value.One group of four coefficient depth of defect obtained by contrast simulation characterizes bent Line, determines the depth of defect.
Multiple index includes flaw echo depth of defect and characterizes coefficient, transmitted wave depth of defect sign coefficient, defect left end mould State conversion echo depth of defect characterizes coefficient, defect right-hand member mode conversion echo depth of defect and characterizes coefficient;Flaw echo defect Depth characterizes coefficient to run into the flaw echo that defect reflection is returned during surface wave propagation, and it is according to for different depth The flaw echo depth of defect sign coefficient that the flaw echo amplitude that defect reflection is returned is different and sets up;Transmitted wave depth of defect Coefficient is characterized to run into the transmitted wave that defect continues to propagate through defect during surface wave propagation, it is according to for different depths The transmitted wave depth of defect spent the transmission wave amplitude difference of defect transmissive and set up characterizes coefficient;Defect left end mode conversion Echo depth of defect characterize coefficient be surface wave in communication process, scattering is had below defect when being interacted with defect Bulk wave, the lower surface that bodily wave propagation runs into workpiece can be reflected, and running into defect can be converted to along the propagation of defect left end once again Surface wave, it is according to the defect left end mode turn that the defect defect left end mode conversion echo amplitude of different depth is different and sets up Gain ripple depth of defect and characterize coefficient;Defect right-hand member mode conversion echo depth of defect sign coefficient is surface wave in communication process In, the bulk wave of scattering is had when being interacted with defect below defect, the lower surface that bodily wave propagation runs into workpiece can be reflected back Come, the surface wave propagated along defect right-hand member can be converted to once again by running into defect, and it is according to the defect defect right-hand member mould of different depth The defect right-hand member mode conversion echo depth of defect that state converts back wave amplitude difference and set up characterizes coefficient.
Brief description of the drawings
Mode conversion echo simulation figure of Fig. 1 surface waves at defect two ends;
Figure 25 00kHz multiple indexs depth of defect characterizes curve;
Fig. 3 tests schematic diagram;
Fig. 4 COMSOL simulation models;
Propagation of Figure 55 00kHz lower surface ripples in aluminium sheet;
Defect clawback and mode conversion echo information under Figure 65 00kHz;
Fig. 7 Electromagnetic Acoustic Transducer schematic diagrames;
Fig. 8 experiments measure mode conversion echo on the left of flaw echo and defect;
Fig. 9 experiments measure mode conversion echo on the right side of transmitted wave and defect;
Embodiment
Present disclosure is described in further detail below in conjunction with 500kHz examples:
Step 1:Establish emulation polyphyly number curve
The parameter informations such as size and material for target workpiece, determine COMSOL simulation softwares, set up suitable emulation Model (such as accompanying drawing 4).According in target workpiece surface wave velocity of wave be 3000m/s, calculatings wavelength be 6mm, determine depth of defect Size range is 0-7.2mm, carries out serial emulation (such as accompanying drawing 5).According to simulation model, the shape information at observation station is extracted, The shape information at observation station is read with MATlab and to carry out Hilbert transform (such as accompanying drawing 6) more accurate to ensure to obtain Amplitude information.Extract the maximum at ripple bag:Turned with the mode at obtained flaw echo value, transmission wave number and defect two ends Handled after gaining wave number divided by the through wave number normalized of benchmark as ordinate, then with depth of defect divided by wavelength normalization As abscissa, set up multiple index depth of defect and characterize curve (such as accompanying drawing 2).
Step 2:Carry out experiment detection
The stimulating frequency that the parameter informations such as size and material according to target workpiece and emulation polyphyly number curve are set up 500kHz and point mode of excitation, selection Electromagnetic Acoustic Transducer (such as accompanying drawing 7) carry out experiment detection.As shown in Figure 3, experimental system Steel plate including 400mm*650mm*25mm, one swash one and receive receiving transducer, oscillograph, excitation set.Experiment measures workpiece, defect Mode conversion echo (such as accompanying drawing 8-9) at clawback, transmitted wave and defect two ends, record obtains four sign depth of defects of the above Waveforms amplitude.
Step 3:Determine depth of defect
Four defect informations of mode conversion echo at defect clawback, transmitted wave and defect two ends measured according to experiment, It is read out and is handled with MATlab, is most worth with amplitude divided by a reference value direct wave obtains one group of four coefficient value.Contrast simulation One group of resulting four coefficient depth of defect characterizes curve, determines the depth 1mm of defect.
A kind of depth of the vertical defect of plate based on surface wave provided above with 500kHz exemplary operations the present invention is examined Survey method is described.The explanation of above example, which is served mainly to facilitate, understands that surface wave is produced in being interacted with defect Mode conversion echo it is very important, returned for the mode conversion added in traditional defect clawback, transmitted wave at defect two ends Ripple can preferably characterize the depth of defect, to break through the limitation of flex point at 0.45 times of wavelength depth of defect;Simultaneously for ability The those skilled in the art in domain, the method according to the present invention is had in variation, this specification in embodiment and scope Appearance should not be construed as limiting the invention.

Claims (1)

1. a kind of multiple index depth detection method of the vertical defect based on surface wave, it is characterised in that:This method includes flow It is as follows:
Step 1:Establish emulation polyphyly number curve
The parameter informations such as size and material for target workpiece, it is determined that suitable simulation software, sets up suitable simulation model; According to required, depth of defect size range is determined, serial emulation is carried out;According to simulation model, the ripple at observation station or face is extracted Shape information is straight with the mode conversion Echo Rating divided by benchmark at resulting flaw echo value, transmission wave number and defect two ends As ordinate after up to wave number normalized, then using depth of defect or depth of defect divided by wavelength normalization processing be used as it is horizontal sit Mark, sets up multiple index depth of defect and characterizes curve;
Step 2:Carry out experiment detection
The stimulating frequency set up according to the size of target workpiece and material parameter information and emulation polyphyly number curve and the side of exciting Formula, the suitable probe of selection is tested;Experimental system includes test piece, a sharp receipts probe, oscillograph, excitation set; One swash one receive probe is connected with test piece, one swash one receive probe is connected by connecting line with excitation set, excitation set and Oscillograph is connected;Experiment measures the mode conversion echo at the defect clawback, transmitted wave and defect two ends of test piece, records The waveforms amplitude of four sign depth of defects more than;
Step 3:Determine depth of defect
Four defect informations of mode conversion echo at defect clawback, transmitted wave and defect two ends measured according to experiment, divided by A reference value direct wave obtains one group of four coefficient value;One group of four coefficient depth of defect obtained by contrast simulation characterizes curve, Determine the depth of defect;
Multiple index includes flaw echo depth of defect and characterizes coefficient, transmitted wave depth of defect sign coefficient, defect left end mode turn Gain ripple depth of defect and characterize coefficient, defect right-hand member mode conversion echo depth of defect sign coefficient;Flaw echo depth of defect Coefficient is characterized to run into the flaw echo that defect reflection is returned during surface wave propagation, it is according to for different deep defects The flaw echo depth of defect sign coefficient that the flaw echo amplitude reflected is different and sets up;Transmitted wave depth of defect is characterized Coefficient is runs into the transmitted wave that defect continues to propagate through defect during surface wave propagation, it is according to scarce for different depth The transmitted wave depth of defect for falling into the transmission wave amplitude difference of transmissive and setting up characterizes coefficient;Defect left end mode conversion echo Depth of defect characterize coefficient be surface wave in communication process, the body of scattering is had when being interacted with defect below defect Ripple, the lower surface that bodily wave propagation runs into workpiece can be reflected, and the table propagated along defect left end can be converted to once again by running into defect Face ripple, it is according to the defect left end mode conversion that the defect defect left end mode conversion echo amplitude of different depth is different and sets up Echo depth of defect characterizes coefficient;Defect right-hand member mode conversion echo depth of defect sign coefficient is surface wave in communication process In, the bulk wave of scattering is had when being interacted with defect below defect, the lower surface that bodily wave propagation runs into workpiece can be reflected back Come, the surface wave propagated along defect right-hand member can be converted to once again by running into defect, and it is according to the defect defect right-hand member mould of different depth The defect right-hand member mode conversion echo depth of defect that state converts back wave amplitude difference and set up characterizes coefficient.
CN201610195694.8A 2016-03-30 2016-03-30 A kind of polyphyly number depth detection method of the vertical defect based on surface wave Expired - Fee Related CN107024535B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610195694.8A CN107024535B (en) 2016-03-30 2016-03-30 A kind of polyphyly number depth detection method of the vertical defect based on surface wave

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610195694.8A CN107024535B (en) 2016-03-30 2016-03-30 A kind of polyphyly number depth detection method of the vertical defect based on surface wave

Publications (2)

Publication Number Publication Date
CN107024535A true CN107024535A (en) 2017-08-08
CN107024535B CN107024535B (en) 2019-08-09

Family

ID=59525087

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610195694.8A Expired - Fee Related CN107024535B (en) 2016-03-30 2016-03-30 A kind of polyphyly number depth detection method of the vertical defect based on surface wave

Country Status (1)

Country Link
CN (1) CN107024535B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108645920A (en) * 2018-04-09 2018-10-12 华南理工大学 A kind of direct wave suppressing method of the rail flaw ultrasonic detection based on denoising and alignment
CN109900788A (en) * 2019-03-18 2019-06-18 阜阳师范学院 A kind of device and measurement method can be used for carrying out the sound absorption characteristics of acoustical material in-site measurement
CN110363767A (en) * 2019-08-09 2019-10-22 中国特种设备检测研究院 A kind of gridding ultrasound tomography detection method of shaft-like workpiece defect
CN113325072A (en) * 2021-04-30 2021-08-31 北京航空航天大学 Metal plate corrosion damage depth evaluation system and method
CN115406383A (en) * 2021-05-28 2022-11-29 中国石油天然气股份有限公司 Method for detecting corrosion depth of storage tank top

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101398411A (en) * 2008-11-07 2009-04-01 哈尔滨工业大学 Rail tread defect rapid scanning and detecting method and device thereof
CN101571513A (en) * 2009-06-16 2009-11-04 北京理工大学 Ultrasonic guided wave detection device for quality evaluation of composite laminated plate
CN102809610A (en) * 2012-06-04 2012-12-05 北京航空航天大学 Phased array ultrasonic testing method based on improved dynamic depth focusing
CN105102975A (en) * 2013-04-02 2015-11-25 杰富意钢铁株式会社 Ultrasonic flaw-detection method and ultrasonic flaw-detection device

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101398411A (en) * 2008-11-07 2009-04-01 哈尔滨工业大学 Rail tread defect rapid scanning and detecting method and device thereof
CN101571513A (en) * 2009-06-16 2009-11-04 北京理工大学 Ultrasonic guided wave detection device for quality evaluation of composite laminated plate
CN102809610A (en) * 2012-06-04 2012-12-05 北京航空航天大学 Phased array ultrasonic testing method based on improved dynamic depth focusing
CN105102975A (en) * 2013-04-02 2015-11-25 杰富意钢铁株式会社 Ultrasonic flaw-detection method and ultrasonic flaw-detection device

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
关建飞等: "表面垂直裂痕诱发瑞利波散射的数值分析", 《应用声学》 *
赵继辰等: "螺旋波纹管导波检测技术的数值仿真和试验研究", 《机械工程学报》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108645920A (en) * 2018-04-09 2018-10-12 华南理工大学 A kind of direct wave suppressing method of the rail flaw ultrasonic detection based on denoising and alignment
CN108645920B (en) * 2018-04-09 2020-12-22 华南理工大学 Denoising and alignment-based direct wave suppression method for ultrasonic flaw detection of steel rail
CN109900788A (en) * 2019-03-18 2019-06-18 阜阳师范学院 A kind of device and measurement method can be used for carrying out the sound absorption characteristics of acoustical material in-site measurement
CN110363767A (en) * 2019-08-09 2019-10-22 中国特种设备检测研究院 A kind of gridding ultrasound tomography detection method of shaft-like workpiece defect
CN110363767B (en) * 2019-08-09 2021-04-02 中国特种设备检测研究院 Gridding ultrasonic tomography detection method for shaft workpiece defects
CN113325072A (en) * 2021-04-30 2021-08-31 北京航空航天大学 Metal plate corrosion damage depth evaluation system and method
CN115406383A (en) * 2021-05-28 2022-11-29 中国石油天然气股份有限公司 Method for detecting corrosion depth of storage tank top
CN115406383B (en) * 2021-05-28 2024-05-28 中国石油天然气股份有限公司 Method for detecting corrosion depth of top of storage tank

Also Published As

Publication number Publication date
CN107024535B (en) 2019-08-09

Similar Documents

Publication Publication Date Title
CN107024535B (en) A kind of polyphyly number depth detection method of the vertical defect based on surface wave
CN102636303B (en) Method for measuring residual stress of thin plating layer based on surface ultrasonic waves
CN104034287B (en) A kind of elastic anisotropy metallic matrix thermal barrier coating thickness ultrasonic measurement method
CN111044613A (en) Metal plate micro-defect detection method based on nonlinear Lamb wave
CN101140265B (en) Method for bias measuring water content of measuring and cracks using non-linearity resonant frequency
CN106767580B (en) The ultrasonic method of defect laying depth in a kind of determining composite layer laminated structure
CN102608212B (en) Method for measuring acoustic impedance and acoustic attenuation of thin layer based on sound pressure reflection coefficient power spectrum
CN110455917B (en) Concrete crack repair quality detection method
CN108225632A (en) A kind of residual stress non-linear ultrasonic detection method
CN103235039A (en) Ultrasonic TOFD (Time of Flight Diffraction) detection method for weld defects based on linear frequency modulation technology
Ma et al. The reflection of guided waves from simple dents in pipes
CN106198739A (en) A kind of TOFD near surface blind region defect location detection method based on shape transformation
CN113029773A (en) Method and system for detecting fatigue degree of material
CN104698089A (en) Ultrasonic relative time propagation technology suitable for inclined crack quantifying and imaging
CN108918667A (en) A kind of wedge defect inspection method
CN103615995A (en) Method for lossless evaluation of thickness of thin cladding layer based on ultrasonic surface waves
CN104111286A (en) Wave velocity and thickness calibrating method for ultrasonic wave welding spot detection
CN106441178A (en) Ultrasonic thickness measuring technology having self-correcting function
CN105738017B (en) Constituent content influences the modification method of assessment metal material skin stress
CN103852521A (en) Method for ultrasonically testing coupling quality of multilayered media
CN113533504B (en) Subsurface crack quantitative measurement method based on laser ultrasonic surface wave frequency domain parameters
CN102980539A (en) Method for measuring thicknesses of metal layer and oxide layer of wall of boiler heating surface tube
CN105044215A (en) Non-destructive material sound velocity field measurement method
CN202947698U (en) Device for measuring thicknesses of metal layer and oxide layer of boiler heating surface tube wall
CN205844271U (en) A kind of detecting system based on difference frequency non-linear ultrasonic detection sheet metal micro-crack

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20190809