CN103486987B - A kind of method utilizing curve to improve ultrasonic thickness measurement precision - Google Patents

A kind of method utilizing curve to improve ultrasonic thickness measurement precision Download PDF

Info

Publication number
CN103486987B
CN103486987B CN201310476895.1A CN201310476895A CN103486987B CN 103486987 B CN103486987 B CN 103486987B CN 201310476895 A CN201310476895 A CN 201310476895A CN 103486987 B CN103486987 B CN 103486987B
Authority
CN
China
Prior art keywords
point
peak
sampled
maximal value
positive
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.)
Expired - Fee Related
Application number
CN201310476895.1A
Other languages
Chinese (zh)
Other versions
CN103486987A (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.)
Shanghai University of Electric Power
Original Assignee
Shanghai University of Electric Power
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 Shanghai University of Electric Power filed Critical Shanghai University of Electric Power
Priority to CN201310476895.1A priority Critical patent/CN103486987B/en
Publication of CN103486987A publication Critical patent/CN103486987A/en
Application granted granted Critical
Publication of CN103486987B publication Critical patent/CN103486987B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

The present invention relates to a kind of method utilizing curve to improve ultrasonic thickness measurement precision, ultrasonic signal is utilized to be similar to straight line within upper and lower minor time slice at zero point, and signal amplitude is higher, more be similar to straight line, by gathering the upper and lower two adjacent sampled points at zero point of the positive peak-peak of echoed signal or negative peak-peak both sides, carry out fitting a straight line by sampled point, calculate unique point respective coordinates in matched curve, and this coordinate is delivered to THICKNESS CALCULATION module calculating testee thickness.The method use curve fitting algorithm compensating sampling precision, makes ultrasonic thickness measurement precision effectively improve when sample clock frequency is constant, to meet the requirement of some occasion to thickness measuring precision, reduces hardware implementing complexity, improves and realize reliability.

Description

A kind of method utilizing curve to improve ultrasonic thickness measurement precision
Technical field
The present invention relates to a kind of method improving ultrasonic thickness measurement precision, particularly a kind of method utilizing curve to improve ultrasonic thickness measurement precision.
Background technology
Ultrasonic thickness test carries out thickness measure according to Ultrasonic Pulse-Echo principle, and when the ultrasonic pulse launched when popping one's head in arrives material interface by testee, pulse is reflected back toward probe.Ultrasonic thickness test calculates by testee upper and lower surface echoed signal the time that ultrasound wave transmits in testee, the speed that this time and ultrasound wave transmit in testee is multiplied, obtains object thickness.Usual way chooses a unique point (zero point or extreme point), calculate twice return this unique point mistiming as shown in Figure 1.
Ultrasonic thickness test precision is determined by two parts precision, and a part is measuring accuracy, is determined by the measurement capability of ultrasonic thickness measurement equipment itself; Another part craft precision, depends on the external factor such as probe, measurand, site environment, temperature.Method involved by this patent is mainly used to the measuring accuracy improving ultrasonic thickness measurement equipment itself.
According to the principle of ultrasonic thickness test, when the velocity of sound is certain, the precision of thickness measure can be converted to the precision of time measurement, and because the time is obtained by statistic sampling clock periodicity, sample clock frequency directly determines ultrasonic thickness measurement precision.Therefore existing carry high-precision method be substantially all from raising sample clock frequency start with.But sampling clock cycle improves more difficult, the occasion that particularly some thickness measuring accuracy requirements are higher, raising sample clock frequency is only relied on to meet the demands.
Summary of the invention
The present invention be directed to only by improving sample clock frequency to improve the requirement problem that thickness measuring precision can not meet precision raising, propose a kind of method utilizing curve to improve ultrasonic thickness measurement precision, use curve fitting algorithm compensating sampling precision, makes ultrasonic thickness measurement precision effectively improve when sample clock frequency is constant.
Technical scheme of the present invention is: a kind of method utilizing curve to improve ultrasonic thickness measurement precision, specifically comprises the steps:
1) sampled point acquisition module chooses acquisition sampled point according to unique point in two echo area for thickness measure presetting: two echo area for thickness measuring are set by gate, gate comprises width and height two parameters, manually arrange according to technological requirement, sampled point obtains and carries out within the scope of gate widths, gate height is defined as gate position ordinate Y, sampled point acquisition module is chosen according to unique point and is gathered synchronization Y<=0 in echo, Y>=0 two sampled points, and preserve;
2) curve fitting module adopts Algorithm of fitting a straight line, and step 1) is obtained two sampled points in echo and fit to straight line, the point handed over ordinate is unique point, obtains the unique point X1 in two echoes and X2;
3) the unique point X1 of two echoes calculated and X2 is sent into THICKNESS CALCULATION module, calculate testee thickness, testee thickness S computing formula is as follows:
S=(| X1-X2|) *t *v, wherein t is the mistiming of Feature point correspondence in twice return, and namely twice return sampled point acquisition time is poor, and V is the velocity of sound in testee.
Described unique point choose point 4 types, to be respectively before positive peak-peak after zero point, positive peak-peak before zero point, negative peak-peak zero point after zero point and negative peak-peak.
Describedly choose that to obtain sampled point step as follows according to unique point:
1), newly to sample beginning, according to gate height judging characteristic vertex type, if gate height is just, is defined as positive peak-peak type unique point, enters step 2), if gate height is negative, is defined as being negative peak-peak type unique point, enters step 7);
2) a certain moment arbitrary neighborhood three sampled point S, are got n(x, y), S n-1(x, y), S n-2(x, y), meets S simultaneously n-2(y) <=0, S n-1(y) >=0, S ny () >=0, incites somebody to action S this moment n-1(x, y), S n-2(x, y) difference assignment is to the sampled point S at zero point before positive maximal value 11(x 11, y 11), S 12(x 12, y 12), enter step 3);
3), continue to get subsequent time three sampled points to compare, as three sampled point first times S simultaneously n(y) <=S n-1(y) <=S n-2time (y), get S this moment n-2(x, y) assignment gives positive maximal value, enters step 4);
4) if positive maximal value is within the scope of gate, and amplitude exceedes gate height, enters step 5), otherwise gets back to step 2);
5), continue to get subsequent time three sampled points to compare, when three sampled points meet S simultaneously n-2(y) >=0, S n-1(y) <=0, S ny () <=0, by S this moment n-1(x, y), S n-2(x, y) difference assignment is to the sampled point S at zero point after positive maximal value 13(x 13, y 13), S 14(x 14, y 14), enter step 6);
6) if this obtains positive maximal value be more than or equal to positive peak-peak, then give positive peak-peak by this positive maximal value assignment obtained, before positive maximal value, the sampled point assignment at zero point is to the sampled point at zero point before positive peak-peak, after positive maximal value, the sampled point assignment at zero point is to the sampled point at zero point after positive peak-peak, otherwise after the sampled point at zero point, positive maximal value, the sampled point at zero point is constant before retaining positive peak-peak, positive peak-peak, judge now whether exceed gate scope, if exceeded, enter step 12), otherwise return step 2);
7) a certain moment arbitrary neighborhood three sampled point S, are got n(x, y), S n-1(x, y), S n-2(x, y), if meet S simultaneously n-2(y) >=0, S n-1(y) <=0, S ny () <=0, incites somebody to action S this moment n-1(x, y), S n-2(x, y) difference assignment is to the sampled point S at zero point before negative maximal value ' 11(x 11, y 11), S ' 12(x 12, y 12), enter step 8);
8), continue to get subsequent time three sampled points to compare, when three sampled points meet S for the first time simultaneously n(y) >=S n-1(y) >=S n-2time (y), get S this moment n-2(x, y) assignment gives negative maximal value, enters step 9);
9) if negative maximal value is within the scope of gate, and amplitude is lower than gate height, enters step 10), otherwise gets back to step 7);
10), continue to get subsequent time three sampled points to compare, when three sampled points meet S simultaneously n-2(y) <=0, S n-1(y) >=0, S ny () >=0, incites somebody to action S this moment n-1(x, y), S n-2(x, y) difference assignment is to the sampled point S at zero point after negative maximal value ' 13(x 13, y 13), S ' 14(x 14, y 14), enter step 11);
11) if this obtains negative maximal value be less than or equal to negative peak-peak, then give negative peak-peak by this negative maximal value assignment obtained, before negative maximal value, the sampled point assignment at zero point is to the sampled point at zero point before negative peak-peak, after negative maximal value, the sampled point assignment at zero point is to the sampled point at zero point after negative peak-peak, otherwise after the sampled point at zero point, negative maximal value, the sampled point at zero point is constant before retaining negative peak-peak, negative peak-peak, judge now whether exceed gate scope, if exceeded, enter step 12), otherwise return step 7);
12) if the unique point needed is zero point before maximal value, then by zero point before positive/negative maximal value, assignment was curve sampled point, if the unique point needed is zero point after maximal value, then by zero point after positive/negative maximal value, assignment was curve sampled point, etc. end to be sampled, sampling end returns step 1), otherwise rests on step 12).
Beneficial effect of the present invention is: the present invention utilizes curve to improve the method for ultrasonic thickness measurement precision, the method use curve fitting algorithm compensating sampling precision, ultrasonic thickness measurement precision is effectively improved when sample clock frequency is constant, to meet the requirement of some occasion to thickness measuring precision, reduce hardware implementing complexity, improve and realize reliability.
Accompanying drawing explanation
Fig. 1 is ultrasonic thickness-measuring method schematic diagram;
Fig. 2 is high-level schematic functional block diagram in curve ultrasonic thickness-measuring method of the present invention;
Fig. 3 is that unique point of the present invention chooses schematic diagram;
Fig. 4 is unique point type schematic diagram of the present invention;
Fig. 5 is sampled point obtaining step schematic diagram of the present invention;
Fig. 6 is thickness measuring method schematic diagram of the present invention.
Embodiment
Curve fitting algorithm is utilized to calculate the mistiming of Feature point correspondence in twice return, this mistiming is multiplied with the velocity of sound and obtains thickness, high-level schematic functional block diagram in curve ultrasonic thickness-measuring method, comprises sampled point acquisition module as shown in Figure 2, curve fitting module and THICKNESS CALCULATION module.
Because ultrasonic signal (sampling period) within upper and lower minor time slice at zero point is similar to straight line, and signal amplitude is higher, more be similar to straight line, therefore unique point is chosen for the zero point of the positive peak-peak of echoed signal or negative peak-peak both sides, curve fitting algorithm adopts fitting a straight line, fitting a straight line needs two sampled points, and these two sampled points are chosen for two neighbouring sampled points of unique point, and unique point chooses schematic diagram as shown in Figure 3.
Sampled point acquisition module is connected with curve fitting module, is the sampled point that the matching of curve fitting module offer curves needs;
The sampled point that curve fitting module utilizes sampled point acquisition module to obtain carries out curve fitting, and calculates unique point respective coordinates in matched curve, and this coordinate is delivered to THICKNESS CALCULATION module calculating testee thickness;
THICKNESS CALCULATION module obtains unique point coordinate from curve fitting module, calculates the mistiming of Feature point correspondence in twice return, is multiplied with the velocity of sound this mistiming and obtains testee thickness.
Be described in detail with regard to the program below:
Sampled point acquisition module chooses acquisition sampled point according to unique point in two echo area for thickness measure presetting.Two echo area for thickness measuring are set by gate, gate comprises width and height two parameters, need manually to arrange according to technological requirement, gate widths (as the scope that double-head arrow in Fig. 1 limits) determines that algorithm performs scope, sampled point and unique point obtain all carries out within the scope of gate widths, and gate height is defined as gate position ordinate.Unique point has chosen Four types, unique point type schematic diagram as shown in Figure 4, be zero point after zero point and negative peak-peak respectively before positive peak-peak after zero point, positive peak-peak before zero point, negative peak-peak, sampled point is obtained by sampled point obtaining step schematic diagram as shown in Figure 5:
Step one, newly to sample beginning, according to gate height judging characteristic vertex type, if gate height is just, as shown in two width figure front in Fig. 4, be defined as positive peak-peak type unique point, enter step 2, if gate height is negative, as shown in two width figure rear in Fig. 4, be defined as being negative peak-peak type unique point, enter step 7;
Step 2, get a certain moment arbitrary neighborhood three sampled point S n(x, y), S n-1(x, y), S n-2(x, y), meets S simultaneously n-2(y) <=0, S n-1(y) >=0, S ny () >=0, incites somebody to action S this moment n-1(x, y), S n-2(x, y) difference assignment is to the sampled point S at zero point before positive maximal value 11(x 11, y 11), S 12(x 12, y 12), enter step 3;
Step 3, continue to get subsequent time three sampled points and compare, as three sampled points first times S simultaneously n(y) <=S n-1(y) <=S n-2time (y), get S this moment n-2(x, y) assignment gives positive maximal value, enters step 4;
If the positive maximal value of step 4 is within the scope of gate, and amplitude exceedes gate height, enters step 5, otherwise gets back to step 2;
Step 5, continue to get subsequent time three sampled points and compare, when three sampled points meet S simultaneously n-2(y) >=0, S n-1(y) <=0, S ny () <=0, by S this moment n-1(x, y), S n-2(x, y) difference assignment is to the sampled point S at zero point after positive maximal value 13(x 13, y 13), S 14(x 14, y 14), enter step 6;
If step 6 this obtain positive maximal value and be more than or equal to positive peak-peak, then give positive peak-peak by this positive maximal value assignment obtained, before positive maximal value, the sampled point assignment at zero point is to the sampled point at zero point before positive peak-peak, the sampled point that after positive maximal value, the sampled point assignment at zero point gives zero point after positive peak-peak, otherwise before retaining positive peak-peak, positive peak-peak after the sampled point at zero point, positive maximal value the sampled point at zero point constant.Judge now whether exceed gate scope, if exceeded, enter step 12, otherwise return step 2;
Step 7, get a certain moment arbitrary neighborhood three sampled point S n(x, y), S n-1(x, y), S n-2(x, y), if meet S simultaneously n-2(y) >=0, S n-1(y) <=0, S ny () <=0, incites somebody to action S this moment n-1(x, y), S n-2(x, y) difference assignment is to the sampled point S at zero point before negative maximal value ' 11(x 11, y 11), S ' 12(x 12, y 12), enter step 8;
Step 8, continue to get subsequent time three sampled points and compare, meet S when three sampled points first times simultaneously n(y) >=S n-1(y) >=S n-2time (y), get S this moment n-2(x, y) assignment gives negative maximal value, enters step 9;
If the negative maximal value of step 9 is within the scope of gate, and amplitude is lower than gate height, enters step 10, otherwise gets back to step 7;
Step 10, continue to get subsequent time three sampled points and compare, when three sampled points meet S simultaneously n-2(y) <=0, S n-1(y) >=0, S ny () >=0, incites somebody to action S this moment n-1(x, y), S n-2(x, y) difference assignment is to the sampled point S at zero point after negative maximal value ' 13(x 13, y 13), S ' 14(x 14, y 14), enter step 11;
If the negative maximal value of this acquisition of step 11 is less than or equal to negative peak-peak, then give negative peak-peak by this negative maximal value assignment obtained, before negative maximal value, the sampled point assignment at zero point is to the sampled point at zero point before negative peak-peak, after negative maximal value, the sampled point assignment at zero point gives the sampled point at zero point after negative peak-peak, otherwise before retaining negative peak-peak, negative peak-peak after the sampled point at zero point, negative maximal value the sampled point at zero point constant.Judge now whether exceed gate scope, if exceeded, enter step 12, otherwise return step 7;
If the unique point that step 12 needs is zero point before maximal value, then by zero point before positive/negative maximal value, assignment was curve sampled point, if the unique point needed is zero point after maximal value, then by zero point after positive/negative maximal value, assignment was curve sampled point, etc. end to be sampled, sampling end returns step one, otherwise rests on step 12.
Curve fitting module adopts Algorithm of fitting a straight line, utilize unique point to choose type and gather zero passage sampled point, the sampled point provided by zero passage sampled point acquisition module simulates straight line, and known features point is the point in matched curve, and unique point ordinate y=0, spy can be calculated and decorate horizontal ordinate x.Two sampled points are defined as (x 1, y 1), (x 2, y 2) fitting a straight line formula is as follows:
y=kx+b
Wherein, k=(y 1 -y 2 )/(x 1 -x 2 ); b=(y 2* x 1 -y 1* x 2 )/(x 1 -x 2 )
As shown in Figure 6, the unique point X1 of two echoes calculated and X2 is sent into THICKNESS CALCULATION module, calculates testee thickness.
THICKNESS CALCULATION module calculates the mistiming of Feature point correspondence in twice return according to X1 and X2, is multiplied with the velocity of sound this mistiming and obtains testee thickness.The definition sampling period is t (second), and in testee, the velocity of sound is V(meter per second), then testee thickness S(rice) computing formula is as follows:
S=(|X1-X2|) *t *V
Thickness measuring method schematic diagram as shown in Figure 6, sampling period is 10ns, in testee, the velocity of sound is 5920m/s, before adopting positive maximal value, zero point is unique point, getting adjacent two echoes is thickness measuring region, two echo area for thickness measuring are set by gate, and sampled point and unique point obtain all carries out within the scope of gate, and sampled point obtains as follows:
Step one, newly to sample beginning, be positive peak-peak type unique point according to gate height judging characteristic vertex type, enter step 2;
Step 2, get a certain moment arbitrary neighborhood three sampled point S n(x, y), S n-1(x, y), S n-2(x, y), if meet S simultaneously n-2(y) <=0, S n-1(y) >=0, S ny () >=0, incites somebody to action S this moment n-1(x, y), S n-2(x, y) difference assignment is to the sampled point S at zero point before positive maximal value 11(x 11, y 11), S 12(x 12, y 12), enter step 3;
Step 3, continue to get subsequent time three sampled points and compare, as three sampled points first times S simultaneously n(y) <=S n-1(y) <=S n-2time (y), get S this moment n-2(x, y) assignment gives positive maximal value, enters step 4;
If the positive maximal value of step 4 is within the scope of gate, and amplitude exceedes gate height, enters step 5, otherwise gets back to step 2;
Step 5, continue to get subsequent time three sampled points and compare, when three sampled points meet S simultaneously n-2(y) >=0, S n-1(y) <=0, S ny () <=0, by S this moment n-1(x, y), S n-2(x, y) difference assignment is to the sampled point S at zero point after positive maximal value 13(x 13, y 13), S 14(x 14, y 14), enter step 6;
If step 6 this obtain positive maximal value and be more than or equal to positive peak-peak, then just this positive maximal value assignment obtained gives positive peak-peak, before positive maximal value, the sampled point assignment at zero point is to the sampled point at zero point before positive peak-peak, the sampled point that after positive maximal value, the sampled point assignment at zero point gives zero point after positive peak-peak, otherwise before retaining positive peak-peak, positive peak-peak after the sampled point at zero point, positive maximal value the sampled point at zero point constant.Judge now whether exceed gate scope, if exceeded, enter step 7, otherwise return step 2;
The unique point of step 7, needs is zero point before maximal value, then by zero point before positive maximal value, assignment was curve sampled point, enter step 8;
Step 8, repeated execution of steps one are to step 7, until sampling terminates.
The sampled point obtaining two echoes is respectively (x 11, y 11), (x 21, y 21), (x 12, y 12), (x 22, y 22), curve fitting module adopts Algorithm of fitting a straight line, and the sampled point utilizing sampled point acquisition module to provide simulates straight line, and with unique point amplitude for Y calculates at its X on straight line.Formula is as follows:
y=kx+b
Wherein, k1=(y 11 -y 21 )/(x 11 -x 21 ), b1=(y 21* x 11 -y 11* x 21 )/(x 11 -x 21 )
k2=(y 12 -y 22 )/(x 12 -x 22 ),b2=(y 22* x 12 -y 12* x 22 )/(x 12 -x 22 )
The unique point X1 of two echoes calculated and X2 is sent into THICKNESS CALCULATION module, calculates testee thickness.
THICKNESS CALCULATION module calculates the mistiming of Feature point correspondence in twice return according to X1 and X2, is multiplied with the velocity of sound this mistiming and obtains testee thickness.The definition sampling period is t (second), and in testee, the velocity of sound is V(meter per second), then testee thickness S(rice) computing formula is as follows:
S=(|X1-X2|) *t *V=5.92 *10 -5 *(|X1-X2|)
But should be understood that, the above-mentioned description for specific embodiment is relatively more detailed, and therefore can not think the restriction to patent right protection domain of the present invention, patent right protection domain of the present invention should be as the criterion with appended claims.

Claims (2)

1. utilize curve to improve a method for ultrasonic thickness measurement precision, it is characterized in that, specifically comprise the steps:
1) sampled point acquisition module chooses acquisition sampled point according to unique point in two echo area for thickness measure presetting: two echo area for thickness measuring are set by gate, gate comprises width and height two parameters, manually arrange according to technological requirement, sampled point obtains and carries out within the scope of gate widths, gate height is defined as gate position ordinate Y, sampled point acquisition module is chosen according to unique point and is gathered synchronization Y<=0 in echo, Y>=0 two sampled points, and preserve, wherein unique point choose point 4 types, zero point respectively before positive peak-peak, zero point after positive peak-peak, zero point after zero point and negative peak-peak before negative peak-peak,
2) curve fitting module adopts Algorithm of fitting a straight line, and step 1) is obtained two sampled points in echo and fit to straight line, the point handed over horizontal ordinate is unique point, obtains the unique point X1 in two echoes and X2;
3) the unique point X1 of two echoes calculated and X2 is sent into THICKNESS CALCULATION module, calculate testee thickness, testee thickness S computing formula is as follows:
S=(| X1-X2|) *t *v, wherein t is the mistiming of Feature point correspondence in twice return, and namely twice return sampled point acquisition time is poor, and V is the velocity of sound in testee.
2. utilize curve to improve the method for ultrasonic thickness measurement precision according to claim 1, it is characterized in that, describedly choose that to obtain sampled point step as follows according to unique point:
1), newly to sample beginning, according to gate height judging characteristic vertex type, if gate height is just, is defined as positive peak-peak type unique point, enters step 2), if gate height is negative, is defined as being negative peak-peak type unique point, enters step 7);
2) a certain moment arbitrary neighborhood three sampled point S, are got n(x, y), S n-1(x, y), S n-2(x, y), meets S simultaneously n-2(y) <=0, S n-1(y) >=0, S ny () >=0, incites somebody to action S this moment n-1(x, y), S n-2(x, y) difference assignment is to the sampled point S at zero point before positive maximal value 11(x 11, y 11), S 12(x 12, y 12), enter step 3);
3), continue to get subsequent time three sampled points to compare, as three sampled point first times S simultaneously n(y) <=S n-1(y) <=S n-2time (y), get S this moment n-2(x, y) assignment gives positive maximal value, enters step 4);
4) if positive maximal value is within the scope of gate, and amplitude exceedes gate height, enters step 5), otherwise gets back to step 2);
5), continue to get subsequent time three sampled points to compare, when three sampled points meet S simultaneously n-2(y) >=0, S n-1(y) <=0, S ny () <=0, by S this moment n-1(x, y), S n-2(x, y) difference assignment is to the sampled point S at zero point after positive maximal value 13(x 13, y 13), S 14(x 14, y 14), enter step 6);
6) if this obtains positive maximal value be more than or equal to positive peak-peak, then give positive peak-peak by this positive maximal value assignment obtained, before positive maximal value, the sampled point assignment at zero point is to the sampled point at zero point before positive peak-peak, after positive maximal value, the sampled point assignment at zero point is to the sampled point at zero point after positive peak-peak, otherwise after the sampled point at zero point, positive maximal value, the sampled point at zero point is constant before retaining positive peak-peak, positive peak-peak, judge now whether exceed gate scope, if exceeded, enter step 12), otherwise return step 2);
7) a certain moment arbitrary neighborhood three sampled point S, are got n(x, y), S n-1(x, y), S n-2(x, y), if meet S simultaneously n-2(y) >=0, S n-1(y) <=0, S ny () <=0, incites somebody to action S this moment n-1(x, y), S n-2(x, y) difference assignment is to the sampled point S at zero point before negative maximal value ' 11(x 11, y 11), S ' 12(x 12, y 12), enter step 8);
8), continue to get subsequent time three sampled points to compare, when three sampled points meet S for the first time simultaneously n(y) >=S n-1(y) >=S n-2time (y), get S this moment n-2(x, y) assignment gives negative maximal value, enters step 9);
9) if negative maximal value is within the scope of gate, and amplitude is lower than gate height, enters step 10), otherwise gets back to step 7);
10), continue to get subsequent time three sampled points to compare, when three sampled points meet S simultaneously n-2(y) <=0, S n-1(y) >=0, S ny () >=0, incites somebody to action S this moment n-1(x, y), S n-2(x, y) difference assignment is to the sampled point S at zero point after negative maximal value ' 13(x 13, y 13), S ' 14(x 14, y 14), enter step 11);
11) if this obtains negative maximal value be less than or equal to negative peak-peak, then give negative peak-peak by this negative maximal value assignment obtained, before negative maximal value, the sampled point assignment at zero point is to the sampled point at zero point before negative peak-peak, after negative maximal value, the sampled point assignment at zero point is to the sampled point at zero point after negative peak-peak, otherwise after the sampled point at zero point, negative maximal value, the sampled point at zero point is constant before retaining negative peak-peak, negative peak-peak, judge now whether exceed gate scope, if exceeded, enter step 12), otherwise return step 7);
12) if the unique point needed is zero point before maximal value, then by zero point before positive/negative maximal value, assignment was curve sampled point, if the unique point needed is zero point after maximal value, then by zero point after positive/negative maximal value, assignment was curve sampled point, etc. end to be sampled, sampling end returns step 1), otherwise rests on step 12).
CN201310476895.1A 2013-10-14 2013-10-14 A kind of method utilizing curve to improve ultrasonic thickness measurement precision Expired - Fee Related CN103486987B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310476895.1A CN103486987B (en) 2013-10-14 2013-10-14 A kind of method utilizing curve to improve ultrasonic thickness measurement precision

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310476895.1A CN103486987B (en) 2013-10-14 2013-10-14 A kind of method utilizing curve to improve ultrasonic thickness measurement precision

Publications (2)

Publication Number Publication Date
CN103486987A CN103486987A (en) 2014-01-01
CN103486987B true CN103486987B (en) 2016-03-30

Family

ID=49827403

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310476895.1A Expired - Fee Related CN103486987B (en) 2013-10-14 2013-10-14 A kind of method utilizing curve to improve ultrasonic thickness measurement precision

Country Status (1)

Country Link
CN (1) CN103486987B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107621638B (en) * 2017-08-01 2020-11-17 昆明理工大学 Method for solving time difference between two pulse signals based on translation signal peak value
WO2020199197A1 (en) * 2019-04-04 2020-10-08 华为技术有限公司 Method and apparatus for processing echo signals

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101614533A (en) * 2008-06-26 2009-12-30 中国科学院金属研究所 A kind of energy accurately measuring thickness of ultrathin workpieces method and instrument
CN202304777U (en) * 2011-10-24 2012-07-04 南通天华和睿科技创业有限公司 Workpiece thickness measurement device
CN102853791A (en) * 2012-02-23 2013-01-02 浙江大学 Method for scanning ultrasonic microscope and measuring thickness, sound velocity, density and attenuation of thin material simultaneously
CN103075981A (en) * 2012-12-31 2013-05-01 汕头市超声仪器研究所有限公司 Ultrasonic thickness measuring method
CN103234494A (en) * 2013-04-09 2013-08-07 浙江大学 Method of measuring thickness of thin materials on the basis of Gaussian echo model

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4392232B2 (en) * 2002-12-13 2009-12-24 日本鉄塔工業株式会社 Film thickness measurement method

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101614533A (en) * 2008-06-26 2009-12-30 中国科学院金属研究所 A kind of energy accurately measuring thickness of ultrathin workpieces method and instrument
CN202304777U (en) * 2011-10-24 2012-07-04 南通天华和睿科技创业有限公司 Workpiece thickness measurement device
CN102853791A (en) * 2012-02-23 2013-01-02 浙江大学 Method for scanning ultrasonic microscope and measuring thickness, sound velocity, density and attenuation of thin material simultaneously
CN103075981A (en) * 2012-12-31 2013-05-01 汕头市超声仪器研究所有限公司 Ultrasonic thickness measuring method
CN103234494A (en) * 2013-04-09 2013-08-07 浙江大学 Method of measuring thickness of thin materials on the basis of Gaussian echo model

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
高温工业管道超声测厚系统;朱君君等;《无损检测》;20130630;第35卷(第6期);第57-60页 *

Also Published As

Publication number Publication date
CN103486987A (en) 2014-01-01

Similar Documents

Publication Publication Date Title
CN102636303B (en) Method for measuring residual stress of thin plating layer based on surface ultrasonic waves
US9370339B2 (en) Systems, methods and computer programs for detection of tissue mechanical property
CN106643939A (en) Method for calculating ultrasonic transmission time through ultrasonic flowmeter
CN102879472A (en) Adaptive steel rail ultrasonic flaw detection method and device based on frequency spectrum recognition
JP2015529795A (en) Method and system for correcting for temperature changes in an ultrasonic inspection system
CN103995263B (en) A kind of ultrasonic ranging method based on sequential
CN104913830A (en) Radar level instrument signal processing method based on curve fitting
CN103075981B (en) A kind of ultrasonic thickness test method
CN103486987B (en) A kind of method utilizing curve to improve ultrasonic thickness measurement precision
CN202814931U (en) Self-adaption ultrasonic rail fault detection device based on frequency spectrum cognition
CN105549027B (en) A kind of distance measuring method and system based on laser pulse shape leading edge detection
CN104048628A (en) Method for accurately detecting thickness through ultrasonic wave equivalent peak value and device using method
CN105180853A (en) Electromagnetic ultrasonic metal material thickness measurement method
CN205352653U (en) Accurate ranging system of supersound in low speed wind tunnel
JP2007215816A (en) Pulse doppler measuring apparatus, its method and its program
CN106289124B (en) A kind of method for following echo in real time of ultrasonic thickness measurement gate
TWI452322B (en) Method and system for detecting object position by using sound wave
CN110440896B (en) Ultrasonic measurement system and measurement method
CN103674177A (en) Signal processing method and device
CN110507361B (en) Shear wave imaging method and system
CN110515082A (en) A kind of automatic-range system based on ultrasound
JP2017063926A (en) Ultrasonic reception signal correction device, ultrasonic measurement device and ultrasonic reception signal correction method
CN202304777U (en) Workpiece thickness measurement device
CN110988148B (en) Measuring system and measuring method for detecting compactness by utilizing ultrasonic waves
CN105300654A (en) Ultrasonic precise range finding system in low speed wind tunnel

Legal Events

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

Granted publication date: 20160330

Termination date: 20181014