CN107037083A - A kind of ultrasonic infrared thermal imagery crack nondestructive detection excitation parameters method for optimizing - Google Patents
A kind of ultrasonic infrared thermal imagery crack nondestructive detection excitation parameters method for optimizing Download PDFInfo
- Publication number
- CN107037083A CN107037083A CN201710237412.0A CN201710237412A CN107037083A CN 107037083 A CN107037083 A CN 107037083A CN 201710237412 A CN201710237412 A CN 201710237412A CN 107037083 A CN107037083 A CN 107037083A
- Authority
- CN
- China
- Prior art keywords
- ultrasonic
- excitation parameters
- infrared thermal
- thermal imagery
- crack
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N25/00—Investigating or analyzing materials by the use of thermal means
- G01N25/72—Investigating presence of flaws
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Abstract
The invention discloses a kind of ultrasonic infrared thermal imagery crack nondestructive detection excitation parameters method for optimizing, comprise the following steps：1）Set up the threedimensional model of Cracked structure part；2）Structural member threedimensional model is imported into finite element software, structure material model is set up；3）To the overall grid division of structural member, by its discretization；4）To setting contact pair between structural member crack surface, and define contact attribute；5）The ultrasonic wave load of loading and boundary condition；6）Analysis and solution is carried out to model；7）Ultrasonic exciting load parameter is changed, step 5 is returned to）；8）Output drive parameter affecting laws are reported；9）According to excitation parameters affecting laws, excitation parameters are optimized with selection, excitation parameters best of breed is obtained.It is preferred that the present invention carries out parameter with soft generation firmly, it is not necessary to uses expensive hardware device, it is not necessary to which the substantial amounts of Preparatory work of experiment of early stage, it is not necessary to the experimentation of very long time and effort consuming, it is convenient and swift, time saving and energy saving to have the advantages that.
Description
Technical field
The present invention relates to ultrasonic infrared thermography non-destructive evaluation field, more particularly to a kind of ultrasonic infrared thermal imagery crack nondestructive inspection
Survey excitation parameters method for optimizing.
Background technology
Ultrasonic infrared thermal imagery is as a new lossless detection method, and its detection speed is fast, clever to crackle defects detection
It is quick high and suitable for metal and nonmetallic materials structural member, with the incomparable advantage of other lossless detection methods.But it is ultrasonic
Infrared thermal imagery method is in engineer applied, the selection still basic experience by operating personnel of ultrasonic exciting parameter.Ultrasonic wave
Excitation be crackle heat thermal excitation source, be the detected key of crack defect, the excitation parameters of ultrasonic wave are influence ultrasounds
The key parameter of wave energy, control the size of ultrasonic energy also influences the detectable rate of crackle simultaneously.Experience is only relied only on to enter
Row judges, can largely improve the loss of crackle, grasps ultrasonic exciting parameter to the affecting laws of crackle heat, pin
To unlike material, different type crack structtire part, Optimum Excitation parameter combination is summed up, carrying out crack detection judgement according to it can
Reduction erroneous judgement, raising crack detection rate.
The content of the invention
In order to solve the above-mentioned technical problem, the present invention provides a kind of step simple, applied widely ultrasonic infrared thermal imagery
Crack nondestructive detects excitation parameters method for optimizing.
Technical proposal that the invention solves the above-mentioned problems is：A kind of ultrasonic infrared thermal imagery crack nondestructive detection excitation parameters are excellent
Choosing method, comprises the following steps：
1）Set up the threedimensional model of Cracked structure part；
2）Structural member threedimensional model is imported into finite element software, structure material model is set up；
3）To the overall grid division of structural member, by its discretization；
4）To setting contact pair between structural member crack surface, and define contact attribute；
5）The ultrasonic wave load of loading and boundary condition；
6）Analysis and solution is carried out to model；
7）Ultrasonic exciting load parameter is changed, often changes and once exports once result and return to step 5）, directly
Finished to parameter change；
8）Output drive parameter affecting laws are reported；
9）According to excitation parameters affecting laws, excitation parameters are optimized with selection, excitation parameters are obtained optimal
Combination.
Above-mentioned ultrasonic infrared thermal imagery crack nondestructive detection excitation parameters method for optimizing, the step 2）In, finite element software
It is the finite element software with solid-head coupled field analysis ability.
Above-mentioned ultrasonic infrared thermal imagery crack nondestructive detection excitation parameters method for optimizing, the step 3）In, structural member is drawn
The size of the grid divided, which should be satisfied with a ultrasonic wavelength, at least divides 20 grids.
Above-mentioned ultrasonic infrared thermal imagery crack nondestructive detection excitation parameters method for optimizing, the step 4）In, it is soft using beformable body
The point face way of contact of body contact, the kinematics and static friction coefficient that wherein Contact Algorithm is selected between penalty function method, contact surface is respectively provided with
For 0.3.
Above-mentioned ultrasonic infrared thermal imagery crack nondestructive detection excitation parameters method for optimizing, the step 5）In, using displacement letter
The mode that ultrasonic exciting is loaded directly on test specimen finite element unit node by number is loaded, ultrasonic wave high frequency sinusoidal letter
Number represents that load-carrying area is 25mm × 25mm.
Above-mentioned ultrasonic infrared thermal imagery crack nondestructive detection excitation parameters method for optimizing, the step 6）Concretely comprise the following steps
(6.1) defined analysis type is Transient transient analyses；
(6.2) solving condition and parameter are set；0.01 second analysis and solution time was set, mode is solved using time integral step-length,
Ultrasonic frequency is 20kHZ, and definition time step is 2E6；
(6.3) analysis and solution process and solving result；
(6.4) judge whether solving result restrains, it is no, then return to step 3）, it is then to continue to judge whether crackle contact surface is worn
Thoroughly, it is no, then continue step 7）, it is then to return to step 4）.
Above-mentioned ultrasonic infrared thermal imagery crack nondestructive detection excitation parameters method for optimizing, the step 7）In, ultrasonic exciting
Parameter includes excitation amplitude, driving frequency, energized position.
The beneficial effects of the present invention are：The present invention initially sets up structure material model, and net is integrally divided to structural member
Lattice, by its discretization；Then to setting contact pair between structural member crack surface and defining contact attribute, the ultrasonic wave load of loading with
And boundary condition；Analysis and solution is carried out to model again；Finally according to each excitation parameters affecting laws output report, excitation is joined
Number optimizes selection, obtains the excitation parameters best of breed that embodiment carries out ultrasonic wave IR thermal imaging inspection.The present invention is with soft
In generation, hard progress parameter was preferred, it is not necessary to use expensive hardware device, it is not necessary to the substantial amounts of Preparatory work of experiment of early stage, it is not necessary to very long
The experimentation of time and effort consuming, it is convenient and swift, time saving and energy saving to have the advantages that.
Brief description of the drawings
Fig. 1 is flow chart of the invention.
Fig. 2 is threedimensional model schematic diagram of the invention.
Fig. 3 is mesh generation schematic diagram of the invention.
Fig. 4 applies schematic diagram for the load and boundary condition of the present invention.
Fig. 5 is energized position parameter distribution schematic diagram of the invention.
Fig. 6 is present invention slit region temperature variation under different excitation amplitudes.
Fig. 7 is present invention slit region temperature variation under different driving frequencies.
Fig. 8 is experiment and the numerical solution comparative result figure of excitation amplitude affecting laws of the present invention.
Embodiment
The present invention is further illustrated with reference to the accompanying drawings and examples.
As shown in figure 1, ultrasonic infrared thermal imagery crack nondestructive detection excitation parameters method for optimizing, is comprised the steps of：
1）Definition unit type is SOLID5, and the unit is Three-Dimensional 8-node hexahedron coupling unit, with three-dimensional thermal field, structure
Field, magnetic field, electric field, the analysis ability of piezoelectric field, and limited coupling can be realized between each field, the structure heat needed for meeting
Field coupling analysis requirement.
The material for defining metal sheet is 45 steel, and 45 Steel material parameters are as shown in table 1.
The threedimensional model of Cracked structure part is set up in 3 d modeling software, panel size size is 210mm × 80mm
× 5mm, crackle runs through width 5um, crack length 5mm at type crackle, V word maximum opens using V words.As shown in Figure 2.
2）Structural member threedimensional model is imported into ANSYS finite element softwares, structure material model is set up.If structural member structure
Simply also model directly can be set up in finite element software, the step for saving importing finite element software.Finite element software is tool
There are the finite element software of solid-head coupled field analysis ability, such as ANSYS, ABAQUS.
3）To the overall grid division of structural member, by its discretization.Mesh generation is the key of whole finite element analysis, grid
The distribution of size, mesh shape and grid directly affects the whether accurate of the result of calculation of simulation analysis.Mesh generation is excessively
When sparse, it will cause result of calculation mistake occur, or even result of calculation can not be obtained；Increasing unit density can improve calculating
As a result accuracy, but when mesh generation too closely when, can cause again calculate overlong time, waste computer resource.It is super
One of the characteristics of sound wave infrared thermography non-destructive evaluation technology be exactly the actuation duration of ultrasonic wave it is short, only Millisecond, so for standard
True propagation of the analog ultrasonic wave in structural member, it is desirable to which the FEM model of structural member at least divides 20 nets in a wavelength
Lattice.The speed that ultrasonic wave is propagated in steel is 5200m/s, and ultrasonic frequency takes 20KHz, the wavelength that ultrasonic wave is propagated in 45 steel
For 0.26m, the maximum mesh size of flat part is 13mm.The temperature change of slit region is the emphasis of finite element analysis, is created
Crack size length be 5mm, the size of mesh opening of slit region need it is sufficiently small could accurate simulation ultrasonic wave in slit region
Propagation, but integral grid divides and can too closely cause to calculate overlong time again, wastes computer resource etc..Take subregion
The scheme of grid division, it is 0.4mm to define slit region size of mesh opening, and it is 3mm to define flat part block mold surface grids size,
3 grids are divided on flat part thickness direction, using sweeping Meshing Method, first opposite grid division, then to flat part
Overall sweeping grid division.Mesh generation result is as shown in Figure 3.
4）To setting contact pair between structural member crack surface, and define contact attribute.The point face contacted using beformable body beformable body
The kinematics and static friction coefficient that the way of contact, wherein Contact Algorithm are selected between penalty function method, contact surface is disposed as 0.3.
5）The ultrasonic wave load of loading and boundary condition, as shown in Figure 4.During actual ultrasound infrared thermal imagery crack detection
Test specimen left and right ends are clamped by fixture, and staff cultivation, including translation are applied to the displacement freedom in the left and right ends face of test specimen model
The free degree and rotational freedom.The original ambient temperature of simulation analysis is set to 27 °C.Ultrasonic exciting is logical in actual experiment
Transducer conversion is crossed, is then delivered to by cumulative bar in test specimen.The mode that ultrasonic exciting is loaded directly into test specimen is taken,
Ultrasonic wave high frequency sinusoidal function representation, ultrasonic wave load is loaded directly into using displacement function load excitation in ANSYS.
6）Analysis and solution is carried out to model.
Step 6）In have comprising following sub-step：
6.1 defined analysis types are Transient transient analyses；
6.2 set solving condition and parameter.Analysis and solution time 0.01s is set, mode is solved using time integral step-length, by
Short in the time of solution, ultrasonic frequency is 20kHZ, in order to which correct analog ultrasonic wave is propagated in flat part, and time step needs
What is set is sufficiently small, and definition time step is 2E6；
6.3 analysis and solution processes, obtain solving result；
6.4 judge whether solving result restrains, no, return to step 3）Mesh generation is re-started, is, crackle contact surface is judged
Whether penetrate, it is no, continue step 7）, it is to return to step 4）Reset contact stiffness and solve.
7）Ultrasonic exciting parameter includes excitation amplitude, driving frequency, energized position, and wherein excitation amplitude span is
12-48um, driving frequency span is 20-30kHz, and energized position selection schematic diagram is as shown in Figure 5.The area of energized position
Size determines that energized position is uniformly distributed in structural member surface to be measured according to the area of section of ultrasonic exciting equipment cumulative bar.
Change excitation amplitude first, be respectively set to 12,18,24,30,36,42,48um, then change driving frequency 20,22,24,
26th, 28,30kHz, then changes energized position, as shown in figure 5, choosing the excitation in grid successively from left to right, from top to bottom
Position, order of preference is a-a '-b-b '-c-c ', totally 15 energized positions.Often change and be once repeated once step 5）With 6), obtain
Obtain a result of calculation.Such as most start excitation amplitude for 12um, driving frequency is 20 kHz, and energized position is a；In excitation frequency
Changed respectively in the case that rate, energized position are constant excitation amplitude for 18,24,30,36,42,48um；Then excitation position is kept
Put constant, driving frequency is set to 22 kHz, respectively change excitation amplitude be 12,18,24,30,36,42,48um ... until excitation
Frequency traversal 20,22,24,26,28,30kHz；Then change energized position, excitation amplitude traversal 12,18,24,30,36,42,
48um, driving frequency traversal 20,22,24,26,28,30kHz.The result sum of so gained is 7*6*15=630.
8）Above-mentioned parameter change has all been calculated, has collected all simulation results, output drive parameter affecting laws
Report.
Gained embodiment excitation amplitude is to the affecting laws of crackle heat as shown in fig. 6, driving frequency is to crackle heat
Affecting laws are as shown in fig. 7, energized position is as shown in Figure 8 to the affecting laws of crackle heat.By ultrasonic wave infrared thermal imagery without
Damage test experience platform and experimental verification is carried out to wherein excitation amplitude parametric variations rule, 45 steel matter flat boards are chosen in experiment
Part splits as experimental subjects, and according to national standard Fatigue Crack Growth Rate of Metallic Materials test method to its prefabricated fatigue
Line, is tested using ultrasonic wave infrared thermal imagery crack detecting method, experimental result and Numerical results comparison diagram such as Fig. 8
It is shown.Fig. 8 shows that, with the increase of excitation amplitude parameter, slit region temperature is constantly raised, experimental result and numerical solution knot
The two ascendant trend of fruit is consistent, and the correctness of excitation amplitude parameter influence rule is demonstrated well and is solved by this method
The feasibility of ultrasonic exciting parameter affecting laws.
9）According to each excitation parameters affecting laws output report, excitation parameters are optimized with selection, ultrasonic wave is obtained
The excitation parameters best of breed of IR thermal imaging inspection.
Claims (7)
1. a kind of ultrasonic infrared thermal imagery crack nondestructive detection excitation parameters method for optimizing, comprises the following steps：
1）Set up the threedimensional model of Cracked structure part；
2）Structural member threedimensional model is imported into finite element software, structure material model is set up；
3）To the overall grid division of structural member, by its discretization；
4）To setting contact pair between structural member crack surface, and define contact attribute；
5）The ultrasonic wave load of loading and boundary condition；
6）Analysis and solution is carried out to model；
7）Ultrasonic exciting load parameter is changed, often changes and once exports once result and return to step 5）, directly
Finished to parameter change；
8）Output drive parameter affecting laws are reported；
9）According to excitation parameters affecting laws, excitation parameters are optimized with selection, excitation parameters best of breed is obtained.
2. ultrasonic infrared thermal imagery crack nondestructive detection excitation parameters method for optimizing according to claim 1, it is characterised in that：
The step 2）In, finite element software is the finite element software with solid-head coupled field analysis ability.
3. ultrasonic infrared thermal imagery crack nondestructive detection excitation parameters method for optimizing according to claim 1, it is characterised in that：
The step 3）In, the size of the grid divided to structural member, which should be satisfied with a ultrasonic wavelength, at least divides 20 nets
Lattice.
4. ultrasonic infrared thermal imagery crack nondestructive detection excitation parameters method for optimizing according to claim 1, it is characterised in that：
The step 4）In, penalty function method is selected in the point face way of contact contacted using beformable body beformable body, wherein Contact Algorithm, contact surface it
Between kinematics and static friction coefficient be disposed as 0.3.
5. ultrasonic infrared thermal imagery crack nondestructive detection excitation parameters method for optimizing according to claim 1, it is characterised in that：
The step 5）In, carried out by the way of ultrasonic exciting is loaded directly on test specimen finite element unit node by displacement function
Loading, ultrasonic wave high frequency sinusoidal function representation, load-carrying area is 25mm × 25mm.
6. ultrasonic infrared thermal imagery crack nondestructive detection excitation parameters method for optimizing according to claim 1, it is characterised in that：
The step 6）Concretely comprise the following steps
(6.1) defined analysis type is Transient transient analyses；
(6.2) solving condition and parameter are set；0.01 second analysis and solution time was set, mode is solved using time integral step-length,
Ultrasonic frequency is 20kHZ, and definition time step is 2E6；
(6.3) analysis and solution process and solving result；
(6.4) judge whether solving result restrains, it is no, then return to step 3）, it is then to continue to judge whether crackle contact surface is worn
Thoroughly, it is no, then continue step 7）, it is then to return to step 4）.
7. ultrasonic infrared thermal imagery crack nondestructive detection excitation parameters method for optimizing according to claim 1, it is characterised in that：
The step 7）In, ultrasonic exciting parameter includes excitation amplitude, driving frequency, energized position.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710237412.0A CN107037083A (en) | 2017-04-12 | 2017-04-12 | A kind of ultrasonic infrared thermal imagery crack nondestructive detection excitation parameters method for optimizing |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710237412.0A CN107037083A (en) | 2017-04-12 | 2017-04-12 | A kind of ultrasonic infrared thermal imagery crack nondestructive detection excitation parameters method for optimizing |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107037083A true CN107037083A (en) | 2017-08-11 |
Family
ID=59535369
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710237412.0A Pending CN107037083A (en) | 2017-04-12 | 2017-04-12 | A kind of ultrasonic infrared thermal imagery crack nondestructive detection excitation parameters method for optimizing |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107037083A (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107607573A (en) * | 2017-09-01 | 2018-01-19 | 沈阳工业大学 | A kind of new alloy hot cracking tendency Forecasting Methodology |
| CN111426919A (en) * | 2020-04-08 | 2020-07-17 | 国网山西省电力公司电力科学研究院 | Basin-type insulator detection device based on laser-induced ultrasound |
| CN111983022A (en) * | 2020-08-26 | 2020-11-24 | 中华人民共和国苏州海关 | Crane hook crack detection method based on combination of COMSOL and laser ultrasound |
| CN116297679A (en) * | 2022-12-22 | 2023-06-23 | 上海尚实航空发动机股份有限公司 | Aircraft monitoring method, device, electronic equipment and storage medium |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101571887A (en) * | 2009-06-16 | 2009-11-04 | 哈尔滨工业大学 | Finite element prediction system for welding and solidifying crack in virtual environment |
| CN106339521A (en) * | 2015-07-15 | 2017-01-18 | 北京卫星环境工程研究所 | Detection method of solder joint void defects based on laser pulse excitation simulation |
-
2017
- 2017-04-12 CN CN201710237412.0A patent/CN107037083A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101571887A (en) * | 2009-06-16 | 2009-11-04 | 哈尔滨工业大学 | Finite element prediction system for welding and solidifying crack in virtual environment |
| CN106339521A (en) * | 2015-07-15 | 2017-01-18 | 北京卫星环境工程研究所 | Detection method of solder joint void defects based on laser pulse excitation simulation |
Non-Patent Citations (3)
| Title |
|---|
| 伍济钢 等: "超声红外热像法的激励参数影响规律研究", 《激光与红外》 * |
| 李赞 等: "基于超声波红外热像的混凝土泵车臂架疲劳裂纹检测的数值模拟分析", 《机械工程师》 * |
| 李赞: "金属结构件裂纹超声红外热像检测及其激励参数影响规律研究", 《中国优秀硕士学位论文全文数据库 工程科技I辑（月刊）》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107607573A (en) * | 2017-09-01 | 2018-01-19 | 沈阳工业大学 | A kind of new alloy hot cracking tendency Forecasting Methodology |
| CN107607573B (en) * | 2017-09-01 | 2019-12-27 | 沈阳工业大学 | Novel method for predicting hot cracking tendency of alloy |
| CN111426919A (en) * | 2020-04-08 | 2020-07-17 | 国网山西省电力公司电力科学研究院 | Basin-type insulator detection device based on laser-induced ultrasound |
| CN111983022A (en) * | 2020-08-26 | 2020-11-24 | 中华人民共和国苏州海关 | Crane hook crack detection method based on combination of COMSOL and laser ultrasound |
| CN116297679A (en) * | 2022-12-22 | 2023-06-23 | 上海尚实航空发动机股份有限公司 | Aircraft monitoring method, device, electronic equipment and storage medium |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| AU2020101234A4 (en) | An Optimization Method for Excitation Parameters of Ultrasonic Infrared Thermography Crack Nondestructive Testing | |
| CN107037083A (en) | A kind of ultrasonic infrared thermal imagery crack nondestructive detection excitation parameters method for optimizing | |
| AU2017245369B2 (en) | Wrinkle characterization and performance prediction for composite structures | |
| US10953608B2 (en) | Structural health monitoring of curved composite structures using ultrasonic guided waves | |
| Javadi et al. | Ultrasonic inspection of a welded stainless steel pipe to evaluate residual stresses through thickness | |
| CN112763578B (en) | In-service integral composite material R area detection method, reference test block and test block manufacturing method | |
| Ciampa et al. | Detection and modelling of nonlinear elastic response in damaged composite structures | |
| Mesnil et al. | Simulation tools for guided wave based structural health monitoring | |
| Burghardt et al. | Robot-operated quality control station based on the UTT method | |
| Bois et al. | Monitoring of laminated composites delamination based on electro-mechanical impedance measurement | |
| Liu et al. | Development of a fast numerical simulator for infrared thermography testing signals of delamination defect in a multilayered plate | |
| CN105259220B (en) | Crackle conductance measurement method based on four terminal DC potential detection signals | |
| Solodov | Resonant ultrasonic activation of damage: application for diagnostic imaging | |
| Bingham et al. | 3D elastodynamic finite integration technique simulation of guided waves in extended built-up structures containing flaws | |
| Chennamsetti | Interaction between the fundamental Lamb modes and the front edge of a crack in a metallic plate | |
| CN114611369B (en) | Multi-physical-field coupled metal ultrasonic welding numerical simulation analysis method | |
| Lai et al. | Quantify Resonance Inspection With Finite Element-Based Modal Analyses | |
| Leckey et al. | Microcracking in composite laminates: Simulation of crack-induced ultrasound attenuation | |
| Ricci et al. | Guided waves in a disbonded honeycomb composite structure | |
| Warren | Nondestructive Testing Methods Aided Via Numerical Computation Models for Various Critical Aerospace Power Generation Systems | |
| Xie et al. | Finite Element Simulation Analysis of Full-focus Ultrasonic Inspection of Steel Plate Weld Defects | |
| Gopalakrishnan et al. | Application of the finite element method in shm | |
| Richard et al. | Advanced software tools for design and implementation of phased array UT inspection techniques on complex components | |
| Hervin | Guided wave propagation and scattering in anisotropic composite structures | |
| Ennis | An Artificial Intelligence Approach to Fatigue Crack Length Estimation From Acoustic Emission Signals |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170811 |
|
| RJ01 | Rejection of invention patent application after publication |