CN108508091A - Detect the method and device of article deterioration parameter - Google Patents

Detect the method and device of article deterioration parameter Download PDF

Info

Publication number
CN108508091A
CN108508091A CN201810502780.8A CN201810502780A CN108508091A CN 108508091 A CN108508091 A CN 108508091A CN 201810502780 A CN201810502780 A CN 201810502780A CN 108508091 A CN108508091 A CN 108508091A
Authority
CN
China
Prior art keywords
coefficient
calibration sample
article
ultrasonic
measured
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
CN201810502780.8A
Other languages
Chinese (zh)
Other versions
CN108508091B (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.)
Guoneng boiler and Pressure Vessel Inspection Co.,Ltd.
Original Assignee
Guodian Boiler Pressure Vessel Inspection Co Ltd
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 Guodian Boiler Pressure Vessel Inspection Co Ltd filed Critical Guodian Boiler Pressure Vessel Inspection Co Ltd
Priority to CN201810502780.8A priority Critical patent/CN108508091B/en
Publication of CN108508091A publication Critical patent/CN108508091A/en
Application granted granted Critical
Publication of CN108508091B publication Critical patent/CN108508091B/en
Active 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
    • 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/44Processing the detected response signal, e.g. electronic circuits specially adapted therefor

Landscapes

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

Abstract

An embodiment of the present invention provides the method and devices that a kind of detection article deteriorates parameter, can carry out non-linear ultrasonic detection to article to be measured, obtain the ultrasonic signal of the article output to be measured;Wavelet transformation is carried out to the ultrasonic signal, obtains the corresponding wavelet coefficient of multiple frequency ranges;The non-linear ultrasonic coefficient of the article to be measured is determined according to the wavelet coefficient;It is brought into the non-linear ultrasonic coefficient as independent variable in first equation corresponding with the material of article to be measured, obtains the value of the first deterioration parameter of the article to be measured.Due to detecting the deterioration parameter of article to be measured present invention uses non-linear ultrasonic, there is no need to be destroyed to article to be measured.Meanwhile the present invention carries out wavelet transformation to ultrasonic signal, the advantage with multiresolution and partial analysis improves the Stability and veracity of detection.

Description

Detect the method and device of article deterioration parameter
Technical field
The present invention relates to ultrasound detection fields, more particularly to detect the method and device of article deterioration parameter.
Background technology
Bear high temperature, high pressure operating mode workpiece largely use metal material, such as power plant, the high temperature and pressure pipe in chemical plant Road, pipe fitting, valve etc..As active time extends, the gradual aging of metal material leads to the safety and stability of metal works It reduces.
Deterioration parameter by detecting metal works can determine the degree of aging of metal works, and then can be in time to old Change more serious metal works timely to be replaced.Existing deterioration parameter detection method is generally sample detection method, i.e.,:It is logical It crosses disruptive method (such as pipe cutting, whole cutting) to sample metal works to be measured, high-temerature creep experiment (examination when then passing through long Periodicity thousand is tested to tens thousand of hours), the research techniques such as short time high temperature tension test examine the deterioration parameter of the workpiece of sampling It surveys.
Since the prior art needs to sample metal works to be measured, it is thus possible to metal works to be measured can be damaged.Meanwhile it removing Other outside conduit component bear high temperature, the component (such as header, valve body, cylinder body) of high pressure is generally not allowed destructive sampling. Certainly, in addition to metal works, the workpiece of other materials is also likely to be present same problem.
Therefore, how under the premise of not damaging workpiece realize deterioration parameter detection be still this field one urgently Technical barrier to be solved.
Invention content
In view of this, the present invention provides a kind of method and device of detection article deterioration parameter, to pass through non-linear ultrasonic Detect the deterioration parameter of article.
In order to achieve the above-mentioned object of the invention, the present invention provides following technical scheme:
A method of detection article deteriorates parameter, including:
Non-linear ultrasonic detection is carried out to article to be measured, obtains the ultrasonic signal of the article output to be measured;
Wavelet transformation is carried out to the ultrasonic signal, obtains the corresponding wavelet coefficient of multiple frequency ranges;
The non-linear ultrasonic coefficient of the article to be measured is determined according to the wavelet coefficient;
The non-linear ultrasonic coefficient is brought into first equation corresponding with the material of article to be measured as independent variable In, obtain the value of the first deterioration parameter of the article to be measured, wherein the dependent variable of first equation is described first bad Change the value of parameter and first equation is corresponding with the first deterioration parameter.
Optionally, before the progress non-linear ultrasonic detection to article to be measured, the method further includes:
Obtain material multiple calibration samples identical with the determinand product, wherein the first deterioration ginseng of each calibration sample Known to several values;
To each calibration sample:Non-linear ultrasonic detection is carried out to the calibration sample, obtains calibration sample output Ultrasonic signal carries out wavelet transformation to the ultrasonic signal of calibration sample output, and the multiple frequency ranges for obtaining the calibration sample correspond to Wavelet coefficient, the non-linear ultrasonic system of the calibration sample is determined according to the corresponding wavelet coefficient of multiple frequency ranges of the calibration sample Number;
Data fitting is carried out according to the value of the first of each calibration sample the deterioration parameter and non-linear ultrasonic coefficient, is determined The value for waiting for the unknowm coefficient in fit equation corresponding with the material of the calibration sample, wherein described to wait for that fit equation is The tie-in equation of first deterioration parameter and non-linear ultrasonic coefficient;
The unknowm coefficient waited in fit equation is revised as identified value, obtains first equation.
Optionally, the wavelet coefficient is Wi, wherein i is the number of wavelet coefficient, and i is natural number, described fixed according to this The corresponding wavelet coefficient of multiple frequency ranges of standard specimen product determines the non-linear ultrasonic coefficient of the calibration sample, including:
According to formula
β=W2/W1 2
γ=W3/W1 3
The non-linear ultrasonic coefficient gamma and β of the calibration sample is calculated, wherein W1It is non-thread to be carried out to the calibration sample The centre frequency F of the transmitting terminal ultrasonic transducer used when property ultrasound detection1The corresponding wavelet coefficient of frequency range at place, W2For frequency Rate F2The corresponding wavelet coefficient of frequency range at place, W3For frequency F3The corresponding wavelet coefficient of frequency range at place, wherein F2=2F1, F3 =3F1
Optionally, described to wait for that fit equation is:
C=k1β+k2γ+b,
Wherein, C is the first deterioration parameter, k1、k2It is the unknowm coefficient for waiting for fit equation with b.
Optionally, the centre frequency of the transmitting terminal ultrasonic transducer used when carrying out non-linear ultrasonic detection is located at section In [2.5,10] MHz,
And/or the first deterioration parameter is:Creep rupture strength, creep strength, precipitate quantity, grain size or hardness.
A kind of device of detection article deterioration parameter, including:Signal obtaining unit, wavelet coefficient obtaining unit, ultrasound system Number obtaining unit and deterioration gain of parameter unit,
It is defeated to obtain the article to be measured for carrying out non-linear ultrasonic detection to article to be measured for the signal obtaining unit The ultrasonic signal gone out;
The wavelet coefficient obtaining unit obtains multiple frequency ranges and corresponds to for carrying out wavelet transformation to the ultrasonic signal Wavelet coefficient;
The ultrasound coefficient obtaining unit, the non-linear ultrasonic for determining the article to be measured according to the wavelet coefficient Coefficient;
The deterioration gain of parameter unit, for using the non-linear ultrasonic coefficient as independent variable bring into it is described to be measured In corresponding first equation of material of article, the value of the first deterioration parameter of the article to be measured is obtained, wherein described first The dependent variable of equation is the described first value for deteriorating parameter and first equation is corresponding with the first deterioration parameter.
Optionally, described device further includes:Sample obtaining unit, sample detection unit, Coefficient Fitting unit and equation obtain Unit is obtained,
The sample obtaining unit, for carrying out non-linear ultrasonic detection to article to be measured in the signal obtaining unit Before, obtain material multiple calibration samples identical with the determinand product, wherein the first deterioration parameter of each calibration sample takes Known to value;
The sample detection unit, for each calibration sample:Non-linear ultrasonic detection is carried out to the calibration sample, The ultrasonic signal for obtaining calibration sample output carries out wavelet transformation to the ultrasonic signal of calibration sample output, it is fixed to obtain this The corresponding wavelet coefficient of multiple frequency ranges of standard specimen product determines that this is fixed according to the corresponding wavelet coefficient of multiple frequency ranges of the calibration sample The non-linear ultrasonic coefficient of standard specimen product;
The Coefficient Fitting unit, for the value and non-linear ultrasonic according to the first of each calibration sample the deterioration parameter Coefficient carries out data fitting, determines the value for waiting for the unknowm coefficient in fit equation corresponding with the material of the calibration sample, Wherein, described to wait for that fit equation is the tie-in equation of the first deterioration parameter and non-linear ultrasonic coefficient;
The equation obtaining unit, for the unknowm coefficient waited in fit equation to be revised as identified value, Obtain first equation.
Optionally, the wavelet coefficient is Wi, wherein i is the number of wavelet coefficient, and i is natural number, the sample detection Unit is specifically configured to:
To each calibration sample:Non-linear ultrasonic detection is carried out to the calibration sample, obtains calibration sample output Ultrasonic signal carries out wavelet transformation to the ultrasonic signal of calibration sample output, and the multiple frequency ranges for obtaining the calibration sample correspond to Wavelet coefficient, according to formula
β=W2/W1 2
γ=W3/W1 3
The non-linear ultrasonic coefficient gamma and β of the calibration sample is calculated, wherein W1It is non-thread to be carried out to the calibration sample The centre frequency F of the transmitting terminal ultrasonic transducer used when property ultrasound detection1The corresponding wavelet coefficient of frequency range at place, W2For frequency Rate F2The corresponding wavelet coefficient of frequency range at place, W3For frequency F3The corresponding wavelet coefficient of frequency range at place, wherein F2=2F1, F3 =3F1
Optionally, described to wait for that fit equation is:
C=k1β+k2γ+b,
Wherein, C is the first deterioration parameter, k1、k2It is the unknowm coefficient for waiting for fit equation with b.
Optionally, the centre frequency of the transmitting terminal ultrasonic transducer used when carrying out non-linear ultrasonic detection is located at section In [2.5,10] MHz,
And/or the first deterioration parameter is:Creep rupture strength, creep strength, precipitate quantity, grain size or hardness.
Through the above technical solutions, a kind of method and device of detection article deterioration parameter provided in an embodiment of the present invention, Non-linear ultrasonic detection can be carried out to article to be measured, obtain the ultrasonic signal of the article output to be measured;To the ultrasound letter Number carry out wavelet transformation, obtain the corresponding wavelet coefficient of multiple frequency ranges;The article to be measured is determined according to the wavelet coefficient Non-linear ultrasonic coefficient;The non-linear ultrasonic coefficient is brought into corresponding with the material of article to be measured as independent variable In one equation, the value of the first deterioration parameter of the article to be measured is obtained.Since present invention uses non-linear ultrasonics to detect The deterioration parameter of article to be measured, there is no need to be destroyed to article to be measured.Meanwhile the present invention carries out small echo change to ultrasonic signal It changes, the advantage with multiresolution and partial analysis improves the Stability and veracity of detection.
Description of the drawings
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below There is attached drawing needed in technology description to be briefly described.
Fig. 1 is the flow chart for the method that a kind of detection article provided in an embodiment of the present invention deteriorates parameter;
Fig. 2 is another flow chart for detecting article and deteriorating the method for parameter provided in an embodiment of the present invention;
Fig. 3 is the structural schematic diagram for the device that a kind of detection article provided in an embodiment of the present invention deteriorates parameter;
Fig. 4 is another structural schematic diagram for detecting article and deteriorating the device of parameter provided in an embodiment of the present invention.
Specific implementation mode
The invention discloses the method and device that a kind of detection article deteriorates parameter, those skilled in the art can use for reference this Literary content is suitably modified technological parameter realization.In particular, it should be pointed out that all similar substitutions and modifications are to art technology It is it will be apparent that they are considered as being included in the present invention for personnel.The method of the present invention and application are by preferable Embodiment is described, related personnel obviously can not depart from the content of present invention, in spirit and scope to side as described herein Method and application are modified or suitably change and combine, to realize and apply the technology of the present invention.
With reference to embodiment, the present invention is further explained:
As shown in Figure 1, a kind of method of detection article deterioration parameter provided in an embodiment of the present invention, may include:
S100, non-linear ultrasonic detection is carried out to article to be measured, obtains the ultrasonic signal of the article output to be measured;
Wherein, article to be measured can be the article of metal material or other materials, and the present invention does not limit herein.
Specifically, the present invention can carry out ultrasound detection by non-linear ultrasonic equipment to article to be measured, this is non-linear super Acoustic equipment includes:Transmitting terminal ultrasonic transducer and receiving terminal ultrasonic transducer.Wherein, transmitting terminal ultrasonic transducer will be non-linear super Acoustical signal is emitted on article to be measured, which is conducted on article to be measured, then receiving end ultrasound Energy converter receives.
Optionally, the centre frequency of transmitting terminal ultrasonic transducer can be located in section [2.5,10] MHz.
Optionally, the centre frequency of receiving terminal ultrasonic transducer can be the two of the centre frequency of transmitting terminal ultrasonic transducer Times.
S200, wavelet transformation is carried out to the ultrasonic signal, obtains the corresponding wavelet coefficient of multiple frequency ranges;
Wherein, wavelet transformation (Wavelet Transform, WT) is a kind of transform analysis method, it is inherited and development The thought of short time discrete Fourier transform localization, while the shortcomings of window size does not change with frequency is overcome again, it is capable of providing one A " T/F " window with frequency shift is the ideal tools for carrying out signal time frequency analysis and processing.Its main feature It is the feature for capableing of abundant outstanding problem some aspects by transformation, the localization of time (space) frequency can be analyzed, be passed through Flexible shift operations gradually carry out multi-scale refinement to signal (function), are finally reached high frequency treatment time subdivision, frequency at low frequency Subdivision, the requirement that can adapt to time frequency signal analysis automatically solve Fourier transform so as to focus on the arbitrary details of signal Difficult problem.
Traditional Fourier transformation method can not observation signal local frequency domain characteristic, so limited and non-in treated length There are certain deficiency in terms of the signal of stable state, it is susceptible to false signal, influences the precision measured.The present invention uses small wavelength-division The method of analysis carries out the time-frequency convert of signal, and the advantage with multiresolution and partial analysis overcomes Fourier transformation not Foot, is focused the signal in certain time section, to obtain the frequency-region signal of multiresolution, efficiently extracts non-thread Property ultrasonic signal in second harmonic and harmonic signal, and deteriorate parameter using the performance of these characteristic signals and sample Between establish tie-in equation, improve the Stability and veracity of measurement.
It is calculated specifically, the small wave converting method that the present invention uses can be specially mallat algorithms or WAVELET PACKET DECOMPOSITION Method.Wherein, the present invention can pass through formula
Calculate wavelet coefficient, wherein W (a, b) is wavelet coefficient, and t is integration variable, and f (t) is characterized the spectrum in spectrogram The original function that line intensity data is constituted, ψa,b(t) it is wavelet function, ψ * are the conjugate function of wavelet function, and a is wavelet function Zooming parameter, b are the translation parameters of wavelet function.
In practical applications, the present invention can obtain the corresponding wavelet coefficient of centre frequency of transmitting terminal ultrasonic transducer, The present invention can also obtain two frequencys multiplication and the corresponding wavelet coefficient of frequency tripling of the centre frequency.
S300, the non-linear ultrasonic coefficient that the article to be measured is determined according to the wavelet coefficient;
Specifically, the present invention can be according to formula
β=W2/W1 2
γ=W3/W1 3
The non-linear ultrasonic coefficient gamma and β of the article to be measured is calculated, wherein W1It is non-thread to be carried out to the article to be measured The centre frequency F of the transmitting terminal ultrasonic transducer used when property ultrasound detection1The corresponding wavelet coefficient of frequency range at place, W2For frequency Rate F2The corresponding wavelet coefficient of frequency range at place, W3For frequency F3The corresponding wavelet coefficient of frequency range at place, wherein F2=2F1, F3 =3F1
S400, it is brought into the non-linear ultrasonic coefficient as independent variable and the material of article to be measured corresponding first In equation, the value of the first deterioration parameter of the article to be measured is obtained, wherein the dependent variable of first equation is described the One deterioration parameter value and first equation with it is described first deteriorate parameter it is corresponding.
Described first, which deteriorates parameter, to be:Creep rupture strength, creep strength, precipitate quantity, grain size or hardness.
Wherein, the first equation can be:C=k1β+k2γ+b,
Wherein, C is the first deterioration parameter, k1、k2With the coefficient that b is the first equation.
A kind of method of detection article deterioration parameter provided in an embodiment of the present invention, can carry out article to be measured non-linear Ultrasound detection obtains the ultrasonic signal of the article output to be measured;Wavelet transformation is carried out to the ultrasonic signal, obtains multiple frequencies The corresponding wavelet coefficient of section;The non-linear ultrasonic coefficient of the article to be measured is determined according to the wavelet coefficient;It will be described non-thread Property ultrasound coefficient brought into first equation corresponding with the material of article to be measured as independent variable, obtain the article to be measured First deterioration parameter value.Due to detecting the deterioration parameter of article to be measured, nothing present invention uses non-linear ultrasonic Article to be measured need to be destroyed.Meanwhile the present invention carries out wavelet transformation to ultrasonic signal, has multiresolution and partial analysis Advantage, improve the Stability and veracity of detection.
In practical applications, the present invention can determine the coefficient of first equation by being detected to calibration sample.Specifically , when the coefficient of the first equation is unknown, which is to wait for fit equation.As shown in Fig. 2, provided in an embodiment of the present invention Another kind detects the method that article deteriorates parameter:
S001, material multiple calibration samples identical with the determinand product are obtained, wherein the first of each calibration sample is bad Known to the value for changing parameter;
In practical applications, the present invention can carry out deterioration for each material respectively parameter is corresponding waits in fit equation The determination of unknowm coefficient.
S002, to each calibration sample:Non-linear ultrasonic detection is carried out to the calibration sample, it is defeated to obtain the calibration sample The ultrasonic signal gone out carries out wavelet transformation to the ultrasonic signal of calibration sample output, obtains multiple frequency ranges of the calibration sample Corresponding wavelet coefficient determines the non-linear super of the calibration sample according to the corresponding wavelet coefficient of multiple frequency ranges of the calibration sample Sonic system number;
Wherein, step S002 carry out used non-linear ultrasonic when non-linear ultrasonic detection frequency can and step The frequency of S100 non-linear ultrasonics used non-linear ultrasonic when detecting is identical.
Further, step S002 obtain after wavelet transformation the corresponding each frequency range of wavelet coefficient can with step S200 into It is identical that the corresponding each frequency range of wavelet coefficient is obtained after row wavelet transformation.
Such as:The corresponding each frequency range of acquisition wavelet coefficient is respectively after step S002 carries out wavelet transformation:(0,3] MHz, (3,6] MHz and (6,9] MHz, then step S200 carry out obtaining the corresponding each frequency range of wavelet coefficient after wavelet transformation be also respectively: (0,3] MHz, (3,6] MHz and (6,9] MHz.
Wherein, the wavelet coefficient can be Wi, wherein i is the number of wavelet coefficient, and i is natural number, according to the calibration The corresponding wavelet coefficient of multiple frequency ranges of sample determines that the detailed process of the non-linear ultrasonic coefficient of the calibration sample may include:
According to formula
β=W2/W1 2
γ=W3/W1 3
The non-linear ultrasonic coefficient gamma and β of the calibration sample is calculated, wherein W1It is non-thread to be carried out to the calibration sample The centre frequency F of the transmitting terminal ultrasonic transducer used when property ultrasound detection1The corresponding wavelet coefficient of frequency range at place, W2For frequency Rate F2The corresponding wavelet coefficient of frequency range at place, W3For frequency F3The corresponding wavelet coefficient of frequency range at place, wherein F2=2F1, F3 =3F1
Wherein, the embodiment of the present invention carries out the centre frequency of the transmitting terminal ultrasonic transducer used when non-linear ultrasonic detection It can be located in section [2.5,10] MHz.
S003, intended according to the value and non-linear ultrasonic coefficient progress data of the first of each calibration sample the deterioration parameter It closes, determines the value for waiting for the unknowm coefficient in fit equation corresponding with the material of the calibration sample, wherein described to wait being fitted Equation is the tie-in equation of the first deterioration parameter and non-linear ultrasonic coefficient;
Wherein, waiting for fit equation can be:
C=k1β+k2γ+b,
Wherein, C is the first deterioration parameter, k1、k2It is the unknowm coefficient for waiting for fit equation with b.
Specifically, the mode of data fitting can there are many, such as least square fitting,
S004, the unknowm coefficient waited in fit equation is revised as to identified value, obtains first equation.
For convenience of understanding, it is exemplified below:
Assuming that there is n material calibration sample identical with the determinand product, the number of calibration sample is respectively 1 to n, respectively The value of first deterioration parameter of calibration sample is respectively C1To Cn
To each calibration sample:Non-linear ultrasonic detection is carried out to the calibration sample, obtains calibration sample output Ultrasonic signal carries out wavelet transformation to the ultrasonic signal of calibration sample output, and the multiple frequency ranges for obtaining the calibration sample correspond to Wavelet coefficient Wi, wherein i is the number of wavelet coefficient, and the number is identical as the number of calibration sample.Specifically, each fixed The W of standard specimen productiCan be it is multiple, it is corresponding from different frequency ranges respectively.According to the corresponding small echo of multiple frequency ranges of the calibration sample Coefficient WiDetermine the non-linear ultrasonic factor beta of the calibration sampleiAnd γi
So far, the value of the first deterioration parameter of material n calibration sample identical with the determinand product can be obtained: C1To Cn, the non-linear ultrasonic coefficient of above-mentioned n calibration sample can also be obtained:β1To βnAnd γ1To γn
In this way, passing through above-mentioned C1To Cn、β1To βn、γ1To γnIt can be carried out data fitting and obtained formula
C=k1β+k2γ+b
In unknowm coefficient k1、k2With the value of b, to obtain being associated with for the first deterioration parameter and non-linear ultrasonic coefficient Equation.After obtaining the tie-in equation, so that it may deteriorate taking for parameter to obtain the first of article to be measured by method shown in Fig. 1 Value.
Corresponding with the above method, the embodiment of the present invention additionally provides a kind of device of detection article deterioration parameter.
As shown in figure 3, a kind of device of detection article deterioration parameter provided in an embodiment of the present invention, may include:Signal Obtaining unit 100, wavelet coefficient obtaining unit 200, ultrasonic coefficient obtaining unit 300 and deterioration gain of parameter unit 400,
The signal obtaining unit 100 obtains the article to be measured for carrying out non-linear ultrasonic detection to article to be measured The ultrasonic signal of output;
Wherein, article to be measured can be the article of metal material or other materials, and the present invention does not limit herein.
Specifically, the present invention can carry out ultrasound detection by non-linear ultrasonic equipment to article to be measured, this is non-linear super Acoustic equipment includes:Transmitting terminal ultrasonic transducer and receiving terminal ultrasonic transducer.Wherein, transmitting terminal ultrasonic transducer will be non-linear super Acoustical signal is emitted on article to be measured, which is conducted on article to be measured, then receiving end ultrasound Energy converter receives.
Optionally, the centre frequency of transmitting terminal ultrasonic transducer can be located in section [2.5,10] MHz.
Optionally, the centre frequency of receiving terminal ultrasonic transducer can be the two of the centre frequency of transmitting terminal ultrasonic transducer Times.
The wavelet coefficient obtaining unit 200 obtains multiple frequency ranges pair for carrying out wavelet transformation to the ultrasonic signal The wavelet coefficient answered;
Wherein, wavelet transformation (Wavelet Transform, WT) is a kind of transform analysis method, it is inherited and development The thought of short time discrete Fourier transform localization, while the shortcomings of window size does not change with frequency is overcome again, it is capable of providing one A " T/F " window with frequency shift is the ideal tools for carrying out signal time frequency analysis and processing.Its main feature It is the feature for capableing of abundant outstanding problem some aspects by transformation, the localization of time (space) frequency can be analyzed, be passed through Flexible shift operations gradually carry out multi-scale refinement to signal (function), are finally reached high frequency treatment time subdivision, frequency at low frequency Subdivision, the requirement that can adapt to time frequency signal analysis automatically solve Fourier transform so as to focus on the arbitrary details of signal Difficult problem.
Traditional Fourier transformation method can not observation signal local frequency domain characteristic, so limited and non-in treated length There are certain deficiency in terms of the signal of stable state, it is susceptible to false signal, influences the precision measured.The present invention uses small wavelength-division The method of analysis carries out the time-frequency convert of signal, and the advantage with multiresolution and partial analysis overcomes Fourier transformation not Foot, is focused the signal in certain time section, to obtain the frequency-region signal of multiresolution, efficiently extracts non-thread Property ultrasonic signal in second harmonic and harmonic signal, and deteriorate parameter using the performance of these characteristic signals and sample Between establish tie-in equation, improve the Stability and veracity of measurement.
It is calculated specifically, the small wave converting method that the present invention uses can be specially mallat algorithms or WAVELET PACKET DECOMPOSITION Method.Wherein, the present invention can pass through formula
Calculating wavelet coefficient, wherein W (a, b) is wavelet coefficient,tFor integration variable, f (t) is characterized the spectrum in spectrogram The original function that line intensity data is constituted, ψa,b(t) it is wavelet function, ψ*For the conjugate function of wavelet function, a is wavelet function Zooming parameter, b are the translation parameters of wavelet function.
In practical applications, the present invention can obtain the corresponding wavelet coefficient of centre frequency of transmitting terminal ultrasonic transducer, The present invention can also obtain two frequencys multiplication and the corresponding wavelet coefficient of frequency tripling of the centre frequency.
The ultrasound coefficient obtaining unit 300, for determining the non-linear of the article to be measured according to the wavelet coefficient Ultrasonic coefficient;
Specifically, ultrasonic coefficient obtaining unit 300 could be provided as according to formula
β=W2/W1 2
γ=W3/W1 3
The non-linear ultrasonic coefficient gamma and β of the article to be measured is calculated, wherein W1It is non-thread to be carried out to the article to be measured The centre frequency F of the transmitting terminal ultrasonic transducer used when property ultrasound detection1The corresponding wavelet coefficient of frequency range at place, W2For frequency Rate F2The corresponding wavelet coefficient of frequency range at place, W3For frequency F3The corresponding wavelet coefficient of frequency range at place, wherein F2=2F1, F3 =3F1
The deterioration gain of parameter unit 400, for using the non-linear ultrasonic coefficient as independent variable bring into it is described In corresponding first equation of material of article to be measured, the value of the first deterioration parameter of the article to be measured is obtained, wherein described The dependent variable of first equation is the described first value for deteriorating parameter and first equation is corresponding with the first deterioration parameter.
Described first, which deteriorates parameter, to be:Creep rupture strength, creep strength, precipitate quantity, grain size or hardness.
Wherein, the first equation can be:C=k1β+k2γ+b,
Wherein, C is the first deterioration parameter, k1、k2With the coefficient that b is the first equation.
A kind of device of detection article deterioration parameter provided in an embodiment of the present invention, can carry out article to be measured non-linear Ultrasound detection obtains the ultrasonic signal of the article output to be measured;Wavelet transformation is carried out to the ultrasonic signal, obtains multiple frequencies The corresponding wavelet coefficient of section;The non-linear ultrasonic coefficient of the article to be measured is determined according to the wavelet coefficient;It will be described non-thread Property ultrasound coefficient brought into first equation corresponding with the material of article to be measured as independent variable, obtain the article to be measured First deterioration parameter value.Due to detecting the deterioration parameter of article to be measured, nothing present invention uses non-linear ultrasonic Article to be measured need to be destroyed.Meanwhile the present invention carries out wavelet transformation to ultrasonic signal, has multiresolution and partial analysis Advantage, improve the Stability and veracity of detection.
In practical applications, the present invention can determine the coefficient of first equation by being detected to calibration sample.Specifically , when the coefficient of the first equation is unknown, which is to wait for fit equation.As shown in figure 4, device basic shown in Fig. 3 On, the device of another detection article deterioration parameter provided in an embodiment of the present invention can also include:Sample obtaining unit 001, Sample detection unit 002, Coefficient Fitting unit 003 and equation obtaining unit 004,
The sample obtaining unit 001, for carrying out non-linear ultrasonic inspection to article to be measured in the signal obtaining unit Before survey, material multiple calibration samples identical with the determinand product are obtained, wherein the first deterioration parameter of each calibration sample Known to value;
In practical applications, the present invention can carry out deterioration for each material respectively parameter is corresponding waits in fit equation The determination of unknowm coefficient.
The sample detection unit 002, for each calibration sample:Non-linear ultrasonic inspection is carried out to the calibration sample It surveys, obtains the ultrasonic signal of calibration sample output, wavelet transformation is carried out to the ultrasonic signal of calibration sample output, is somebody's turn to do The corresponding wavelet coefficient of multiple frequency ranges of calibration sample, being determined according to the corresponding wavelet coefficient of multiple frequency ranges of the calibration sample should The non-linear ultrasonic coefficient of calibration sample;
Wherein, sample detection unit 002 carry out the frequency of used non-linear ultrasonic when non-linear ultrasonic detection can be with It is identical with the used frequency of non-linear ultrasonic when the progress non-linear ultrasonic detection of signal obtaining unit 100.
Further, sample detection unit 002 obtain after wavelet transformation the corresponding each frequency range of wavelet coefficient can with it is small The corresponding each frequency range of acquisition wavelet coefficient is identical after wave system number obtaining unit 200 carries out wavelet transformation.
Such as:The corresponding each frequency range of acquisition wavelet coefficient is respectively after sample detection unit 002 carries out wavelet transformation:(0, 3] MHz, (3,6] MHz and (6,9] MHz, then wavelet coefficient obtaining unit 200 carry out wavelet transformation after obtain wavelet coefficient correspond to Each frequency range be also respectively:(0,3] MHz, (3,6] MHz and (6,9] MHz.
Wherein, the wavelet coefficient is Wi, wherein i is the number of wavelet coefficient, and i is natural number, the sample detection list Member 002 can be specifically set for:
To each calibration sample:Non-linear ultrasonic detection is carried out to the calibration sample, obtains calibration sample output Ultrasonic signal carries out wavelet transformation to the ultrasonic signal of calibration sample output, and the multiple frequency ranges for obtaining the calibration sample correspond to Wavelet coefficient, according to formula
β=W2/W1 2
γ=W3/W1 3
The non-linear ultrasonic coefficient gamma and β of the calibration sample is calculated, wherein W1It is non-thread to be carried out to the calibration sample The centre frequency F of the transmitting terminal ultrasonic transducer used when property ultrasound detection1The corresponding wavelet coefficient of frequency range at place, W2For frequency Rate F2The corresponding wavelet coefficient of frequency range at place, W3For frequency F3The corresponding wavelet coefficient of frequency range at place, wherein F2=2F1, F3 =3F1
Wherein, the embodiment of the present invention carries out the centre frequency of the transmitting terminal ultrasonic transducer used when non-linear ultrasonic detection It can be located in section [2.5,10] MHz.
The Coefficient Fitting unit 003, for according to the first of each calibration sample the deterioration value of parameter and non-linear Ultrasonic coefficient carries out data fitting, determines taking for the unknowm coefficient waited in fit equation corresponding with the material of the calibration sample Value, wherein described to wait for that fit equation is the tie-in equation of the first deterioration parameter and non-linear ultrasonic coefficient;
The equation obtaining unit 004, for the unknowm coefficient waited in fit equation to be revised as identified take Value obtains first equation.
Wherein, waiting for fit equation can be:
C=k1β+k2γ+b,
Wherein, C is the first deterioration parameter, k1、k2It is the unknowm coefficient for waiting for fit equation with b.
Specifically, the mode of data fitting can there are many, such as least square fitting.
It should be noted that herein, relational terms such as first and second and the like are used merely to a reality Body or operation are distinguished with another entity or operation, are deposited without necessarily requiring or implying between these entities or operation In any actual relationship or order or sequence.Moreover, the terms "include", "comprise" or its any other variant are intended to Non-exclusive inclusion, so that the process, method, article or equipment including a series of elements is not only wanted including those Element, but also include other elements that are not explicitly listed, or further include for this process, method, article or equipment Intrinsic element.In the absence of more restrictions, the element limited by sentence "including a ...", it is not excluded that There is also other identical elements in process, method, article or equipment including the element.
Each embodiment in this specification is all made of relevant mode and describes, identical similar portion between each embodiment Point just to refer each other, and each embodiment focuses on the differences from other embodiments.Especially for system reality For applying example, since it is substantially similar to the method embodiment, so description is fairly simple, related place is referring to embodiment of the method Part explanation.
The foregoing is merely illustrative of the preferred embodiments of the present invention, is not intended to limit the scope of the present invention.It is all Any modification, equivalent replacement, improvement and so within the spirit and principles in the present invention, are all contained in protection scope of the present invention It is interior.

Claims (10)

1. a kind of method of detection article deterioration parameter, which is characterized in that including:
Non-linear ultrasonic detection is carried out to article to be measured, obtains the ultrasonic signal of the article output to be measured;
Wavelet transformation is carried out to the ultrasonic signal, obtains the corresponding wavelet coefficient of multiple frequency ranges;
The non-linear ultrasonic coefficient of the article to be measured is determined according to the wavelet coefficient;
It brings into first equation corresponding with the material of article to be measured, obtains using the non-linear ultrasonic coefficient as independent variable Obtain the value of the first deterioration parameter of the article to be measured, wherein the dependent variable of first equation is the first deterioration ginseng Several values and first equation are corresponding with the first deterioration parameter.
2. according to the method described in claim 1, it is characterized in that, carrying out non-linear ultrasonic detection to article to be measured described Before, the method further includes:
Obtain material multiple calibration samples identical with the determinand product, wherein the first deterioration parameter of each calibration sample Known to value;
To each calibration sample:Non-linear ultrasonic detection is carried out to the calibration sample, obtains the ultrasound of calibration sample output Signal carries out wavelet transformation to the ultrasonic signal of calibration sample output, and the multiple frequency ranges for obtaining the calibration sample are corresponding small Wave system number determines the non-linear ultrasonic coefficient of the calibration sample according to the corresponding wavelet coefficient of multiple frequency ranges of the calibration sample;
Data fitting, determining and institute are carried out according to the value of the first of each calibration sample the deterioration parameter and non-linear ultrasonic coefficient State the corresponding value for waiting for the unknowm coefficient in fit equation of material of calibration sample, wherein described to wait for that fit equation is first Deteriorate the tie-in equation of parameter and non-linear ultrasonic coefficient;
The unknowm coefficient waited in fit equation is revised as identified value, obtains first equation.
3. according to the method described in claim 2, it is characterized in that, the wavelet coefficient is Wi, wherein i is the volume of wavelet coefficient Number, i is natural number, described to determine the non-linear of the calibration sample according to the corresponding wavelet coefficient of multiple frequency ranges of the calibration sample Ultrasonic coefficient, including:
According to formula
β=W2/W1 2
γ=W3/W1 3
The non-linear ultrasonic coefficient gamma and β of the calibration sample is calculated, wherein W1It is non-linear super to be carried out to the calibration sample The centre frequency F of the transmitting terminal ultrasonic transducer used when sound detection1The corresponding wavelet coefficient of frequency range at place, W2For frequency F2 The corresponding wavelet coefficient of frequency range at place, W3For frequency F3The corresponding wavelet coefficient of frequency range at place, wherein F2=2F1, F3= 3F1
4. according to the method described in claim 3, it is characterized in that, described wait for that fit equation is:
C=k1β+k2γ+b,
Wherein, C is the first deterioration parameter, k1、k2It is the unknowm coefficient for waiting for fit equation with b.
5. method according to claim 1 to 4, which is characterized in that
The centre frequency for carrying out the transmitting terminal ultrasonic transducer used when non-linear ultrasonic detection is located at section [2.5,10] MHz It is interior,
And/or the first deterioration parameter is:Creep rupture strength, creep strength, precipitate quantity, grain size or hardness.
6. a kind of device of detection article deterioration parameter, which is characterized in that including:Signal obtaining unit, wavelet coefficient obtain single Member, ultrasonic coefficient obtaining unit and deterioration gain of parameter unit,
The signal obtaining unit obtains the article output to be measured for carrying out non-linear ultrasonic detection to article to be measured Ultrasonic signal;
It is corresponding small to obtain multiple frequency ranges for carrying out wavelet transformation to the ultrasonic signal for the wavelet coefficient obtaining unit Wave system number;
The ultrasound coefficient obtaining unit, the non-linear ultrasonic system for determining the article to be measured according to the wavelet coefficient Number;
The deterioration gain of parameter unit, for being brought into the non-linear ultrasonic coefficient as independent variable and the article to be measured Corresponding first equation of material in, obtain the value of the first deterioration parameter of the article to be measured, wherein first equation Dependent variable be it is described first deteriorate parameter value and first equation with it is described first deteriorate parameter it is corresponding.
7. device according to claim 6, which is characterized in that described device further includes:Sample obtaining unit, sample detection Unit, Coefficient Fitting unit and equation obtaining unit,
The sample obtaining unit, for before the signal obtaining unit carries out non-linear ultrasonic detection to article to be measured, obtaining Obtain material multiple calibration samples identical with the determinand product, wherein the value of the first deterioration parameter of each calibration sample is Know;
The sample detection unit, for each calibration sample:Non-linear ultrasonic detection is carried out to the calibration sample, is obtained The ultrasonic signal of calibration sample output carries out wavelet transformation to the ultrasonic signal of calibration sample output, obtains the calibration sample The corresponding wavelet coefficient of multiple frequency ranges of product determines the calibration sample according to the corresponding wavelet coefficient of multiple frequency ranges of the calibration sample The non-linear ultrasonic coefficient of product;
The Coefficient Fitting unit, for the value and non-linear ultrasonic coefficient according to the first of each calibration sample the deterioration parameter Data fitting is carried out, determines the value for waiting for the unknowm coefficient in fit equation corresponding with the material of the calibration sample, wherein It is described to wait for that fit equation is the tie-in equation of the first deterioration parameter and non-linear ultrasonic coefficient;
The equation obtaining unit is obtained for the unknowm coefficient waited in fit equation to be revised as identified value First equation.
8. device according to claim 7, which is characterized in that the wavelet coefficient is Wi, wherein i is the volume of wavelet coefficient Number, i is natural number, and the sample detection unit is specifically configured to:
To each calibration sample:Non-linear ultrasonic detection is carried out to the calibration sample, obtains the ultrasound of calibration sample output Signal carries out wavelet transformation to the ultrasonic signal of calibration sample output, and the multiple frequency ranges for obtaining the calibration sample are corresponding small Wave system number, according to formula
β=W2/W1 2
γ=W3/W1 3
The non-linear ultrasonic coefficient gamma and β of the calibration sample is calculated, wherein W1It is non-linear super to be carried out to the calibration sample The centre frequency F of the transmitting terminal ultrasonic transducer used when sound detection1The corresponding wavelet coefficient of frequency range at place, W2For frequency F2 The corresponding wavelet coefficient of frequency range at place, W3For frequency F3The corresponding wavelet coefficient of frequency range at place, wherein F2=2F1, F3= 3F1
9. device according to claim 8, which is characterized in that described to wait for that fit equation is:
C=k1β+k2γ+b,
Wherein, C is the first deterioration parameter, k1、k2It is the unknowm coefficient for waiting for fit equation with b.
10. the device according to any one of claim 6 to 9, which is characterized in that
The centre frequency for carrying out the transmitting terminal ultrasonic transducer used when non-linear ultrasonic detection is located at section [2.5,10] MHz It is interior,
And/or the first deterioration parameter is:Creep rupture strength, creep strength, precipitate quantity, grain size or hardness.
CN201810502780.8A 2018-05-23 2018-05-23 Method and device for detecting article degradation parameters Active CN108508091B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810502780.8A CN108508091B (en) 2018-05-23 2018-05-23 Method and device for detecting article degradation parameters

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810502780.8A CN108508091B (en) 2018-05-23 2018-05-23 Method and device for detecting article degradation parameters

Publications (2)

Publication Number Publication Date
CN108508091A true CN108508091A (en) 2018-09-07
CN108508091B CN108508091B (en) 2021-03-02

Family

ID=63401115

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810502780.8A Active CN108508091B (en) 2018-05-23 2018-05-23 Method and device for detecting article degradation parameters

Country Status (1)

Country Link
CN (1) CN108508091B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111141822A (en) * 2020-01-06 2020-05-12 大连理工大学 HP type furnace tube high-temperature tissue degradation nondestructive evaluation method based on nonlinear torsional waves

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1793227A1 (en) * 2005-12-01 2007-06-06 General Electric Company Ultrasonic inspection system and method using a nonlinearity parameter
CN102621227A (en) * 2012-04-25 2012-08-01 北京交通大学 Nonlinear ultrasonic evaluation method of metal material plating damage
CN103278567A (en) * 2013-05-13 2013-09-04 北京交通大学 Ultrasonic nondestructive evaluation of early-stage fatigue damage of bonding interface
CN103323525A (en) * 2013-05-21 2013-09-25 广东电网公司电力科学研究院 Nonlinear ultrasound system used for detecting super austenitic stainless steel inter-granular corrosion
CN104359977A (en) * 2014-10-22 2015-02-18 北京理工大学 Acoustic surface wave high-order nonlinear parameter representation method for bending fatigue state of metal plate
CN105067703A (en) * 2015-07-31 2015-11-18 华东理工大学 Plate-like metal member/material remaining life non-destructive detection and evaluation method
CN106461614A (en) * 2014-05-28 2017-02-22 汉阳大学校产学协力团 Device for evaluating deterioration and estimating strength by using ultrasound waves and method for evaluating deterioration and estimating strength by using same

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1793227A1 (en) * 2005-12-01 2007-06-06 General Electric Company Ultrasonic inspection system and method using a nonlinearity parameter
CN102621227A (en) * 2012-04-25 2012-08-01 北京交通大学 Nonlinear ultrasonic evaluation method of metal material plating damage
CN103278567A (en) * 2013-05-13 2013-09-04 北京交通大学 Ultrasonic nondestructive evaluation of early-stage fatigue damage of bonding interface
CN103323525A (en) * 2013-05-21 2013-09-25 广东电网公司电力科学研究院 Nonlinear ultrasound system used for detecting super austenitic stainless steel inter-granular corrosion
CN106461614A (en) * 2014-05-28 2017-02-22 汉阳大学校产学协力团 Device for evaluating deterioration and estimating strength by using ultrasound waves and method for evaluating deterioration and estimating strength by using same
CN104359977A (en) * 2014-10-22 2015-02-18 北京理工大学 Acoustic surface wave high-order nonlinear parameter representation method for bending fatigue state of metal plate
CN105067703A (en) * 2015-07-31 2015-11-18 华东理工大学 Plate-like metal member/material remaining life non-destructive detection and evaluation method

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
AMIR MOSTAVI等: "Wavelet based harmonics decomposition of ultrasonic signal in assessment of plastic strain in aluminum", 《MEASUREMENT》 *
JASON MAXIMINO C. ONGPENG等: "Investigation on the Sensitivity of Ultrasonic Test Applied to Reinforced Concrete Beams Using Neural Network", 《APPL. SCI.》 *
江念等: "脉冲反转和小波变换在非线性超声检测中的应用", 《仪表技术与传感器》 *
阎红娟: "金属构件疲劳损伤非线性超声检测方法研究", 《中国博士论文全文数据库 工程科技I辑》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111141822A (en) * 2020-01-06 2020-05-12 大连理工大学 HP type furnace tube high-temperature tissue degradation nondestructive evaluation method based on nonlinear torsional waves

Also Published As

Publication number Publication date
CN108508091B (en) 2021-03-02

Similar Documents

Publication Publication Date Title
CN110108802B (en) Carrier modulation nonlinear ultrasonic guided wave damage detection method
Hong et al. The matching pursuit approach based on the modulated Gaussian pulse for efficient guided-wave damage inspection
Yu et al. Advanced signal processing for enhanced damage detection with piezoelectric wafer active sensors
Chen et al. Noise analysis of digital ultrasonic nondestructive evaluation system
WO2010127429A1 (en) Ultrasonic scanning system and ultrasound image enhancement method
EP1910795A2 (en) System and methods for non-destructive testing of tubular systems
CN114910565B (en) Method for correcting relative nonlinear coefficient in nonlinear ultrasonic detection
KR101218616B1 (en) Calibration method of contract transducer for absolute measurement of nonlinearity parameter, apparatus for calibration by using the method, and method and apparatus for absolute measurement of the parameter by using the method
CN108508091A (en) Detect the method and device of article deterioration parameter
CN108918678A (en) Utilize the method and device of non-linear ultrasonic detection article deterioration parameter
CN116698986B (en) Sound absorbing material quality control system and method based on tuned band-pass filter
WO2007132509A1 (en) Ultrasonic inspection method utilizing resonant phenomena
CN103954628A (en) Ensemble empirical mode decomposition (EEMD) and approximate entropy combined steel tube damage monitoring method
Shi et al. Signal matching wavelet for ultrasonic flaw detection in high background noise
CN108195934B (en) Ultrasonic guided wave detection frequency optimization method based on time-frequency analysis
Liang et al. Optimal scale wavelet transform for the identification of weak ultrasonic signals
JP2005221321A (en) Method and device for detecting ultrasonic signal
CN107515251A (en) A kind of detecting system for metal material
Kari et al. Characterization of a cylindrical rod by inversion of acoustic scattering data
Michaels et al. Methodologies for quantifying changes in diffuse ultrasonic signals with applications to structural health monitoring
WO2010135822A1 (en) Determination of fractional compositions using nonlinear spectrophonometry
CN116026933A (en) Method for determining detection resolution and detection sensitivity of nonlinear ultrasonic detection system
Fright et al. Reconstruction of the complex reflectance function in acoustic microscopy
Chen et al. Investigation of pipelines defect localization for fusion reactor by using T (0, 1) mode ultrasonic guided waves
Kari et al. Characterization of immersed transversely isotropic rods by inversion of acoustic scattering data

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
CP01 Change in the name or title of a patent holder

Address after: No. 2096, Beijing New Energy Technology Research Institute

Patentee after: Guoneng boiler and Pressure Vessel Inspection Co.,Ltd.

Address before: No. 2096, Beijing New Energy Technology Research Institute

Patentee before: GUODIAN BOILER AND PRESSURE VESSEL INSPECTION Co.,Ltd.

CP01 Change in the name or title of a patent holder