CN108918668A - Composite material ellipse damage reason location detection method based on public circumcircle - Google Patents

Composite material ellipse damage reason location detection method based on public circumcircle Download PDF

Info

Publication number
CN108918668A
CN108918668A CN201810467444.4A CN201810467444A CN108918668A CN 108918668 A CN108918668 A CN 108918668A CN 201810467444 A CN201810467444 A CN 201810467444A CN 108918668 A CN108918668 A CN 108918668A
Authority
CN
China
Prior art keywords
piezoelectric patches
signal
damage
piezoelectric
coordinate
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
CN201810467444.4A
Other languages
Chinese (zh)
Other versions
CN108918668B (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.)
Zhengzhou University
Original Assignee
Zhengzhou 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 Zhengzhou University filed Critical Zhengzhou University
Priority to CN201810467444.4A priority Critical patent/CN108918668B/en
Publication of CN108918668A publication Critical patent/CN108918668A/en
Application granted granted Critical
Publication of CN108918668B publication Critical patent/CN108918668B/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/01Indexing codes associated with the measuring variable
    • G01N2291/011Velocity or travel time
    • 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
    • G01N2291/0231Composite or layered materials
    • 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/028Material parameters
    • G01N2291/0289Internal structure, e.g. defects, grain size, texture
    • 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

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

Composite material ellipse damage reason location detection method based on public circumcircle, the detection device that the detection method uses includes a square monitoring region being made of four piezoelectric transducers, and single piezoelectric transducer is set as excitation piezoelectric patches in center, remaining four piezoelectric transducer is detected as piezoelectric patches is received by the damage reason location detection method.The invention is using form square monitoring region and piezoelectric patches is motivated to combine, and it is combined using public circumcircle with ellipse positioning, according to damage intersecting point coordinate, public circumcircle fitting is carried out with least square method, it can determine damage position and degree of injury according to the central coordinate of circle of public circumscribed fitting circle and diameter information, compensating for conventional elliptical positioning mode can only be positioned merely, the shortcomings that can not carrying out degree of injury assessment, provides a kind of simple and fast method for the judgement of degree of injury.

Description

Composite material ellipse damage reason location detection method based on public circumcircle
Technical field
The invention belongs to technical field of nondestructive testing, specifically more particularly to a kind of composite wood based on public circumcircle Expect oval damage reason location detection method.
Background technique
Lamb wave is a kind of elastic guided wave transmitted in the solid panel or layer structure with free boundary, due to it The propagation characteristic of body such as decays along propagation path small, and energy loss is small, and long transmission distance etc. becomes existing field of non destructive testing The important means of research.Composite material exists at present since it has the characteristics that specific strength is high, specific stiffness is high, corrosion-resistant, antifatigue The aircraft of aerospace field is largely used.But composite airplane component easily generates during manufacture and military service The internal tiny flaw such as layering and crackle, in order to meet aircraft industry high quality product, high reliability, high-performance and long-life It is required that, it is necessary to comprehensive detection is carried out to composite airplane component.For being manufactured by carbon fiber enhancement resin base composite material Aerospace structure, Lamb wave application it is more universal.
Such as oval location technology of traditional damage detecting method mainly constitutes ellipse according to damage characteristic signal, and passes through song The mode of line intersection point carries out the positioning of structural damage, and positioning accuracy is directly proportional to number of sensors and signal processing precision, and And this method is simple positioning, the continuous improvement with Practical Project to technical requirements, to the qualitative assessment of degree of injury It is imperative that technology study.
Summary of the invention
The purpose of the invention is to realize to carry out positioning and degree of injury assessment to damage of composite materials, positioning is improved Precision proposes a kind of composite material ellipse position finding and detection method based on public circumcircle.
To achieve the goals above, technical solution provided by the invention is:
Composite material ellipse damage reason location detection method based on public circumcircle, the detection device which uses Including a square monitoring region being made of four piezoelectric transducers, and in center, single piezoelectric type sensing is set Device is as excitation piezoelectric patches, remaining four piezoelectric transducer is as reception piezoelectric patches;The damage reason location detection method include with Lower step:
1) using the center of composite panel to be detected as origin, rectangular coordinate system is established, excitation piezoelectric patches is located at origin, Four reception piezoelectric patches constitute monitoring region, and center is overlapped with origin, and the centre coordinate of each piezoelectric patches is recorded as (xi,yi), Wherein i is the number of each piezoelectric patches, and the label of piezoelectric patches is motivated to be denoted as 0, and the label number for receiving piezoelectric patches is respectively 1,2,3,4;
2) when testing, signal excitation is applied to excitation piezoelectric patches, first under the lossless health status of composite panel, successively The signal for acquiring and collecting each reception piezoelectric patches feedback, is recorded asThen under composite panel faulted condition, successively acquire And the signal of each reception piezoelectric patches feedback is collected, it is recorded asAnd before and after composite material dash-board injury, other test conditions are constant, Wherein, it is to the pumping signal u (t) of excitation piezoelectric patches application:
U (t)=A [H (t)-H (t-n/fc)](1-cos(2πfct/n))sin2πfct
Wherein, H (t) is Heaviside step function,
A is that the amplitude of signal is modulated,
fcFor signal center frequency,
N is signal wave crest number;
3) to the collected health signal of each reception piezoelectric patchesAnd damage signalCarry out difference signal analysis;
4) the damage characteristic signal and signal direct wave in difference signal are found out, according to signal envelope find out its it is corresponding when Between, it is denoted as T respectivelyi(i=1,2,3,4), Tj(i=1,2,3,4) available signal is in two paths:Motivate piezoelectric patches → damage Hurt position → reception piezoelectric patches, motivate the time difference Δ t=T propagated in piezoelectric patches → reception piezoelectric patchesi-Tj
5) it with the coordinate of each piezoelectric transducer and the time difference parameter and signal velocity of combination step 4), substitutes into and calculates Elliptical orbit is constructed in method formula, algorithmic formula is:
Wherein, υ is Lamb wave velocity of wave,
dmnIt is that signal arrives the flight time for receiving piezoelectric patches by damage reflection again,
(x, y) is damage coordinate, (xm,ym) it is excitation piezoelectric patches coordinate, (xn,yn) it is to receive piezoelectric patches coordinate;
6) four elliptical orbits are constructed according to four reception piezoelectric patches, finds out intersecting point coordinate, and with intersection point discreteness sieve Select damage intersection point.
7) according to gained intersecting point coordinate data in step 6), circumcircle fitting, the fitting circle are carried out according to least square method It is tangent with four ellipses of building;
8) central coordinate of circle information and diameter information are obtained according to public circumscribed fitting circle, determines damage position, and Carry out the assessment of degree of injury.
The piezoelectric transducer uses circular piezoelectric ceramic piece, and each piezoelectric ceramic piece parameter is identical.
Composite material ellipse damage reason location detection method provided by the invention based on public circumcircle is supervised using square The form surveyed region and piezoelectric patches is motivated to combine, and combined using public circumcircle with ellipse positioning, it is sat according to damage intersection point Mark carries out public circumcircle fitting with least square method, according to the central coordinate of circle and diameter information of public circumscribed fitting circle It can determine damage position and degree of injury, compensating for conventional elliptical positioning mode can only be positioned merely, can not be damaged The shortcomings that scale evaluation, provides a kind of simple and fast method for the judgement of degree of injury.
Specific embodiment
Detailed description of the preferred embodiments below.It should be understood that described herein specific Embodiment is merely illustrative and explains the present invention, is not intended to restrict the invention.
In the present invention, the composite material ellipse damage reason location detection method based on public circumcircle, the detection method is used Detection device include being made of to be square damage monitoring region and by single piezoelectric sensing four reception piezoelectric transducers The excitation piezoelectric patches that device is constituted, and motivate piezoelectric patches set on the center in square damage monitoring region, remaining four piezoelectricity Formula sensor is as reception piezoelectric patches.The present invention uses that detection method includes the following steps:
1) using the center of composite panel to be detected as origin, rectangular coordinate system is established, excitation piezoelectric patches is located at origin, Four reception piezoelectric patches constitute monitoring region, and center is overlapped with origin, and the centre coordinate of each piezoelectric patches is recorded as (xi,yi), Wherein i is the number of each piezoelectric patches, and the label of piezoelectric patches is motivated to be denoted as 0, and the label number for receiving piezoelectric patches is respectively 1,2,3,4;
2) when testing, signal excitation is applied to excitation piezoelectric patches, first under the lossless health status of composite panel, successively The signal for acquiring and collecting each reception piezoelectric patches feedback, is recorded asThen under composite panel faulted condition, successively acquire And the signal of each reception piezoelectric patches feedback is collected, it is recorded asAnd before and after composite material dash-board injury, other test conditions are constant, Wherein, it is to the pumping signal u (t) of excitation piezoelectric patches application:
U (t)=A [H (t)-H (t-n/fc)](1-cos(2πfct/n))sin2πfct
Wherein, H (t) is Heaviside step function, and A is that the amplitude of signal is modulated, fcFor signal center frequency, n is letter Number wave crest number;
3) the collected health signal that each reception piezoelectric patches is receivedAnd damage signalCarry out difference signal analysis;
4) the damage characteristic signal and signal direct wave in difference signal are found out, according to signal envelope find out its it is corresponding when Between, it is denoted as T respectivelyi(i=1,2,3,4), Tj(i=1,2,3,4) available signal is in two paths:That is damage characteristic signal By excitation piezoelectric patches → damage position → reception piezoelectric patches, signal direct wave by propagating in excitation piezoelectric patches → reception piezoelectric patches Time difference Δ t=Ti-Tj
5) it with the coordinate of each piezoelectric transducer and the time difference parameter and signal velocity of combination step 4), substitutes into and calculates Elliptical orbit is constructed in method formula, algorithmic formula is:
Wherein, v is Lamb wave velocity of wave, dmnIt is that signal reflects the flight time for arriving receiver piezoelectric patches again by damage, (x, Y), (xm,ym), (xn,yn) it is divided into damage coordinate, piezoelectric patches coordinate is motivated, piezoelectric patches coordinate is received;
6) four elliptical orbits are constructed according to four reception piezoelectric patches, finds out intersecting point coordinate, and with intersection point discreteness sieve Select damage intersection point;
7) according to gained intersecting point coordinate data in step 6), circumcircle fitting, the fitting circle are carried out according to least square method It is tangent with four ellipses of building;
8) central coordinate of circle information and diameter information are obtained according to public circumscribed fitting circle, determines damage position, and Carry out the assessment of degree of injury.
The invention discloses the composite material ellipse damage reason location detection methods based on public circumcircle, firstly, design pressure Electric chip arrays, by the square composition damage monitoring region of piezoelectric transducer;Secondly, under health status, it is multiple to being bonded in Excitation piezoelectric patches on condensation material plate applies excitation, and the response signal for receiving piezoelectric patches feedback to four is acquired respectively;Again Secondary, under faulted condition, the response signal for receiving piezoelectric patches feedback to four is acquired, and respectively to health status and damage The signal of state carries out difference signal processing, and finds out its damage characteristic signal and direct-path signal and its corresponding time, calculates Its time difference out;Then, according to processing data, signal velocity, elliptical building is carried out according to oval positioning principle, and count Calculate four elliptical damage intersecting point coordinates;Finally, it is quasi- to carry out public circumcircle with least square method according to damage intersecting point coordinate It closes, damage position and degree of injury can determine according to the central coordinate of circle of public circumscribed fitting circle and diameter information.
Composite material ellipse damage reason location detection method provided by the invention based on public circumcircle is supervised using square The form surveyed region and piezoelectric patches is motivated to combine, and combined using public circumcircle with ellipse positioning, compensate for conventional elliptical The shortcomings that positioning mode can only be positioned merely, and degree of injury assessment can not be carried out, provides one kind for the judgement of degree of injury Simple and fast method.
The above-described embodiments merely illustrate the principles and effects of the present invention, and is not intended to limit the present invention.It is any ripe The personage for knowing this technology all without departing from the spirit and scope of the present invention, carries out modifications and changes to above-described embodiment.Cause This, institute is complete without departing from the spirit and technical ideas disclosed in the present invention by those of ordinary skill in the art such as At all equivalent modifications or change, should be covered by the claims of the present invention.

Claims (2)

1. the composite material ellipse damage reason location detection method based on public circumcircle, it is characterised in that:The detection method uses Detection device include be made of four piezoelectric transducers one square monitoring region, and it is single in center setting Piezoelectric transducer is as excitation piezoelectric patches, remaining four piezoelectric transducer is as reception piezoelectric patches;Damage reason location detection Method includes the following steps:
1) using the center of composite panel to be detected as origin, rectangular coordinate system is established, excitation piezoelectric patches is located at origin, and four It receives piezoelectric patches and constitutes monitoring region, center is overlapped with origin, and the centre coordinate of each piezoelectric patches is recorded as (xi,yi), wherein I is the number of each piezoelectric patches, and the label of piezoelectric patches is motivated to be denoted as 0, and the label number for receiving piezoelectric patches is respectively 1,2,3,4;
2) when testing, signal excitation is applied to excitation piezoelectric patches, first under the lossless health status of composite panel, is successively acquired And the signal of each reception piezoelectric patches feedback is collected, it is recorded asThen under composite panel faulted condition, successively acquires and receive The signal for collecting each reception piezoelectric patches feedback, is recorded asAnd before and after composite material dash-board injury, other test conditions are constant, wherein It is to the pumping signal u (t) for motivating piezoelectric patches to apply:
U (t)=A [H (t)-H (t-n/fc)](1-cos(2πfct/n))sin2πfct
Wherein, H (t) is Heaviside step function,
A is that the amplitude of signal is modulated,
fcFor signal center frequency,
N is signal wave crest number;
3) to the collected health signal of each reception piezoelectric patchesAnd damage signalCarry out difference signal analysis;
4) the damage characteristic signal and signal direct wave in difference signal are found out, its corresponding time is found out according to signal envelope, point T is not denoted as iti(i=1,2,3,4), Tj(i=1,2,3,4) available signal is in two paths:Motivate piezoelectric patches → damage position → receive piezoelectric patches, motivate the time difference Δ t=T propagated in piezoelectric patches → reception piezoelectric patchesi-Tj
5) with the coordinate of each piezoelectric transducer and the time difference parameter and signal velocity of combination step 4), it is public to substitute into algorithm Elliptical orbit is constructed in formula, algorithmic formula is:
Wherein, v is Lamb wave velocity of wave,
dmnIt is that signal arrives the flight time for receiving piezoelectric patches by damage reflection again,
(x, y) is damage coordinate, (xm,ym) it is excitation piezoelectric patches coordinate, (xn,yn) it is to receive piezoelectric patches coordinate;
6) four elliptical orbits are constructed according to four reception piezoelectric patches, finds out intersecting point coordinate, and filter out with intersection point discreteness Damage intersection point.
7) according to gained intersecting point coordinate data in step 6), circumcircle fitting, the fitting circle and structure are carried out according to least square method Four ellipses built are tangent;
8) central coordinate of circle information and diameter information are obtained according to public circumscribed fitting circle, determines damage position, and carry out The assessment of degree of injury.
2. according to the composite material ellipse damage reason location detection method described in claim 1 based on public circumcircle is utilized, It is characterized in that:The piezoelectric transducer uses circular piezoelectric ceramic piece, and each piezoelectric ceramic piece parameter is identical.
CN201810467444.4A 2018-05-16 2018-05-16 Composite material elliptical damage positioning detection method based on common circumscribed circle Expired - Fee Related CN108918668B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810467444.4A CN108918668B (en) 2018-05-16 2018-05-16 Composite material elliptical damage positioning detection method based on common circumscribed circle

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810467444.4A CN108918668B (en) 2018-05-16 2018-05-16 Composite material elliptical damage positioning detection method based on common circumscribed circle

Publications (2)

Publication Number Publication Date
CN108918668A true CN108918668A (en) 2018-11-30
CN108918668B CN108918668B (en) 2020-12-18

Family

ID=64402986

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810467444.4A Expired - Fee Related CN108918668B (en) 2018-05-16 2018-05-16 Composite material elliptical damage positioning detection method based on common circumscribed circle

Country Status (1)

Country Link
CN (1) CN108918668B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110376282A (en) * 2019-07-15 2019-10-25 北京航空航天大学 A kind of Lamb wave damage reason location method based on oval probability and Bayesian Estimation
CN115797436A (en) * 2023-01-31 2023-03-14 深圳市优品文创科技有限公司 Jewelry production bead diameter identification system based on picture self-filling

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8286492B2 (en) * 2009-12-09 2012-10-16 The Boeing Company Mode decomposition of sound waves using amplitude matching
CN104965025A (en) * 2015-05-14 2015-10-07 南京航空航天大学 Multi-zone damage detection method based on Lamb wave signal correlation coefficient

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8286492B2 (en) * 2009-12-09 2012-10-16 The Boeing Company Mode decomposition of sound waves using amplitude matching
CN104965025A (en) * 2015-05-14 2015-10-07 南京航空航天大学 Multi-zone damage detection method based on Lamb wave signal correlation coefficient

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
刘钊等: "板内缺陷超声导波阵列检测技术 ", 《仪表技术与传感器》 *
姜跃栋等: "基于Lamb波的金属薄板载荷定位方法 ", 《无损检测》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110376282A (en) * 2019-07-15 2019-10-25 北京航空航天大学 A kind of Lamb wave damage reason location method based on oval probability and Bayesian Estimation
CN110376282B (en) * 2019-07-15 2021-06-08 北京航空航天大学 Lamb wave damage positioning method based on ellipse probability and Bayesian estimation
CN115797436A (en) * 2023-01-31 2023-03-14 深圳市优品文创科技有限公司 Jewelry production bead diameter identification system based on picture self-filling
CN115797436B (en) * 2023-01-31 2023-11-17 深圳市优品文创科技有限公司 Pearl diameter recognition system of pearl machine for jewelry production based on picture self-filling

Also Published As

Publication number Publication date
CN108918668B (en) 2020-12-18

Similar Documents

Publication Publication Date Title
CN106813819B (en) A kind of method of supercritical ultrasonics technology measurement curve surface work pieces residual stress
CN102043016B (en) Lamb wave-based autonomous damage identification imaging method
EP1228364B1 (en) System and method of assessing the structural properties of wooden members using ultrasound
US7891247B2 (en) Method and system for detecting an anomaly and determining its size
CN106706760A (en) Acoustic emission source positioning method of composite material plate of omnidirectional dual circular array
CN106802202B (en) A method of measurement anisotropic material plane stress
CN110687208A (en) Non-reference Lamb wave damage monitoring method based on hyperbolic positioning
US7237438B1 (en) Systems and methods for determining the velocity of ultrasonic surface skimming longitudinal waves on various materials
EP0794430A2 (en) Method and apparatus for evaluation and inspection of composite-repaired structures
CN108918668A (en) Composite material ellipse damage reason location detection method based on public circumcircle
CN114235971A (en) Ultrasonic guided wave-based method for identifying damaged shape of woven composite material plate
CN105651862B (en) A kind of device and method that damage check being carried out to composite material using blue nurse guided wave
US6920790B2 (en) Apparatus for in-situ nondestructive measurement of Young's modulus of plate structures
Lin et al. Investigation of acoustic emission source localization performance on the plate structure using piezoelectric fiber composites
CN108562528A (en) A kind of composite material porosity evaluation method based on acoustic emission
Han et al. The propagation of Lamb waves in a laminated composite plate with a variable stepped thickness
CN115494160B (en) Phased array air coupling transducer and plane stress ultrasonic measurement method and device thereof
CN106441176A (en) Air coupling ultrasonic wave thin film thickness detection method
CN105806949A (en) Correcting method for radial travel time in ultrasonic CT detection
CN111189914B (en) Method for determining gradient thickness coefficient for ultrasonic detection of composite material
US6575036B1 (en) Method for in-situ nondestructive measurement of Young's modulus of plate structures
Van Way et al. Aircraft structural health monitoring system development: Overview of the Air Force/Navy smart metallic structures program
US11060860B2 (en) Method of inspection by guided waves
CN113311072A (en) Steel rail stress detection method and system
US20210237120A1 (en) Multi-frequency wireless sensor

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
CB03 Change of inventor or designer information

Inventor after: Li Cheng

Inventor after: Zhao Zhujun

Inventor after: Zheng Yanping

Inventor after: Song Chengjie

Inventor after: Zhao Jiangming

Inventor after: Tie Ying

Inventor after: Zhao Huadong

Inventor after: Hou Yuliang

Inventor after: Duan Yuechen

Inventor after: Xia Xiaosong

Inventor before: Zheng Yanping

Inventor before: Zhao Zhujun

Inventor before: Li Cheng

Inventor before: Tie Ying

Inventor before: Hou Yuliang

Inventor before: Duan Yuechen

Inventor before: Zhao Huadong

Inventor before: Zhao Jiangming

Inventor before: Song Chengjie

Inventor before: Xia Xiaosong

CB03 Change of inventor or designer information
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20201218

CF01 Termination of patent right due to non-payment of annual fee