CN113218560A - Ultrasonic real-time estimation method for bolt pretightening force - Google Patents

Ultrasonic real-time estimation method for bolt pretightening force Download PDF

Info

Publication number
CN113218560A
CN113218560A CN202110420543.9A CN202110420543A CN113218560A CN 113218560 A CN113218560 A CN 113218560A CN 202110420543 A CN202110420543 A CN 202110420543A CN 113218560 A CN113218560 A CN 113218560A
Authority
CN
China
Prior art keywords
bolt
calibration
ultrasonic
fft
face
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
CN202110420543.9A
Other languages
Chinese (zh)
Other versions
CN113218560B (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.)
China Yangtze Power Co Ltd
Original Assignee
China Yangtze Power 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 China Yangtze Power Co Ltd filed Critical China Yangtze Power Co Ltd
Priority to CN202110420543.9A priority Critical patent/CN113218560B/en
Publication of CN113218560A publication Critical patent/CN113218560A/en
Application granted granted Critical
Publication of CN113218560B publication Critical patent/CN113218560B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01LMEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L5/00Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes
    • G01L5/24Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes for determining value of torque or twisting moment for tightening a nut or other member which is similarly stressed
    • G01L5/246Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes for determining value of torque or twisting moment for tightening a nut or other member which is similarly stressed using acoustic waves

Landscapes

  • Physics & Mathematics (AREA)
  • Acoustics & Sound (AREA)
  • General Physics & Mathematics (AREA)
  • Force Measurement Appropriate To Specific Purposes (AREA)

Abstract

The invention discloses an ultrasonic real-time estimation method for bolt pretightening force, which comprises the following steps: selecting a calibration bolt, and obtaining a reference waveform w under a zero-stress state0(ii) a Performing calibration test on the calibration bolt to obtain calibration data sigmaiAnd wi(ii) a To w0Fast Fourier transform and conjugate calculation to obtain FFT (w)0)*(ii) a To wiFast Fourier transform to obtain FFT (w)i) (ii) a FFT (w)0) And FFT (w)0)*Element multiplication to obtain FFT (R (w)0,wi) Compensating for M-1 times of zero to obtain FFT (R)M(w0,wi) ); performing inverse Fourier transform to obtain RM(w0,wi) Searching the corresponding position of the maximum value element to obtain the waveform time delay delta ti(ii) a Obtaining a bolt pretightening force coefficient A by using the calibration data; obtaining a measurement waveform; the ultrasonic time delay is obtained and is substituted into a bolt pretightening force formula to obtain the pretightening force, and the bolt pretightening force is measured by adopting the method, so that the measuring speed can be effectively improved.

Description

Ultrasonic real-time estimation method for bolt pretightening force
Technical Field
The invention relates to the technical field of ultrasonic nondestructive testing, in particular to an ultrasonic bolt pretightening force real-time estimation method.
Background
The bolt connection is an important connection mode in engineering production, has the advantages of simple structure, convenience in assembly and disassembly, high efficiency, low cost and the like, and is widely applied to important equipment such as aerospace, traffic bridges, building structures, chemical products and the like.
In actual engineering construction, improper bolt pretightening force can damage the quality of a bolt connecting pair and directly influence the integrity and reliability of equipment. Therefore, the accurate, intuitive and convenient measurement of the bolt pretightening force plays a crucial role in the wide application of the bolt pretightening force in the industrial field. The traditional bolt pretightening force measuring method comprises a torque pulling method, a resistance strain gauge method and the like, but most methods have certain limitations in the aspect of bolt pretightening force monitoring due to the problems of poor control precision, low measuring efficiency and the like. In recent years, an ultrasonic measurement method in a nondestructive testing technology is rapidly developed in bolt pretightening force measurement application, and the theoretical basis of measuring pretightening force based on acoustic elasticity by an ultrasonic method is also perfected through years of research. The ultrasonic longitudinal wave method has extremely high sensitivity to the cylindrical bolt due to transmission along the central axis, and is widely applied to engineering measurement of bolt stress.
In the fields of important engineering such as water conservancy and hydropower and aerospace, bolt pretightening force needs to be monitored in real time (more than or equal to 10Hz), which requires that an ultrasonic bolt pretightening force estimation algorithm has higher calculation speed. In the existing ultrasonic bolt pretightening force measurement scheme, the sampling digit of the adopted data acquisition card is mostly 200MHz, if an interpolation algorithm is not adopted, the pretightening force resolution is only 5MPa, the requirement of the engineering measurement resolution of 1MPa is not met, and in order to further reduce the pretightening force resolution, a tri-spline interpolation algorithm is usually adopted in the prior art method, but the calculation complexity of the interpolation algorithm is about O (N ^2), when the data point number N is large, the calculation time is raised quadratically, the calculation amount is huge, and the real-time requirement cannot be met. Therefore, the ultrasonic real-time bolt pretightening force estimation algorithm has important application significance in the aspect of bolt pretightening force detection in the field of major engineering.
Disclosure of Invention
The invention aims to overcome the defects and provide an ultrasonic real-time estimation method for the pretightening force of the bolt, which can simply, conveniently, real-timely, low-cost and high-precision measure the pretightening force of the bolt and effectively improve the measuring speed of the pretightening force of the bolt.
In order to solve the technical problems, the invention adopts the technical scheme that: an ultrasonic real-time estimation method for bolt pretightening force comprises the following steps:
s1: selecting any bolts with the same specification and the same batch of bolts to be measured as calibration bolts, and measuring the calibration bolts through ultrasonic waves under the zero-stress state to obtain a reference waveform w0
S2: performing calibration test on the calibration bolt to obtain calibration data sigmaiAnd wiWhere σ isiIndicating the magnitude of the pretension, wiRepresenting ultrasonic waveforms in corresponding pretightening force states;
s3: to w0Performing fast Fourier transform to obtain reference waveform frequency domain data FFT (w)0) And find FFT (w)0) Conjugation yields FFT (w)0)*
S4: to wiPerforming fast Fourier transform to obtain a calibration waveform frequency domain array FFT (wi);
S5: FFT (w)0) And FFT (w)0)*Element multiplication is carried out to obtain FFT (R (w)0,wi) And FFT (R (w)) is performed on the data according to the interpolation multiple M0,wi) Zero-1 times of the tail part to obtain a data set FFT (R)M(w0,wi));
S6: FFT (R) on the data setM(w0,wi) ) inverse Fourier transform to obtain cross-correlation array RM(w0,wi) Looking up RM(w0,wi) The maximum value element corresponds to the position maxPos, and w is obtained according to a time delay estimation formula delta T ═ maxPos delta T/M0,wiTime delay Δ t of waveform therebetweeniWherein Δ T is a hardware sampling time interval;
s7: fitting a bolt pretightening force formula by using the calibration data by using a least square method to obtain a bolt pretightening force coefficient A;
s8: measuring the bolt to be measured through ultrasonic waves to obtain a measured waveform wtest
S9: repeating the steps S4-S7 to obtain a measured waveform wtestAnd the reference waveform w0Ultrasonic time delay Δ t therebetweentestAnd substituting the bolt pre-tightening force coefficient A obtained by calculation in the S7 into a bolt pre-tightening force formula to obtain the pre-tightening force of the bolt to be tested.
Preferably, the expression of the bolt pretension formula is as follows:
σ=A·Δt
in the formula, A is the bolt pre-tightening force coefficient, sigma is the bolt pre-tightening force, and delta t is the ultrasonic time delay.
Preferably, the specific operation step of step S2 includes:
s2-1: determining the clamping length of the calibration bolt according to the actual fastening condition, and clamping the calibration bolt on a stretcher;
s2-2: in the elastic stage of the calibration bolt, applying different pretightening forces sigma to the calibration bolt by using a stretcheriWhere i is 1, 2, …, n (n ≧ 5), n denoting the number of different stress states applied by the stretcher;
s2-3: acquiring ultrasonic signals under corresponding pretightening force to obtain a calibration waveform wi
Preferably, the specific operation step of step S7 includes:
s7-1: will calibrate Δ t in the dataiAs an argument, σiPerforming linear fitting to obtain a linear fitting equation sigma k · Δ t + b as a dependent variable;
s7-2: judging whether the intercept b of the linear fitting equation meets the effective condition of the calibration data, if not, performing calibration data acquisition on the calibration bolt again, and updating the calibration data until the intercept b of the linear fitting equation meets the effective condition of the calibration data;
s7-3: if the intercept b meets the effective condition of the calibration data, the coefficient of the bolt pretension force is A ═ k;
preferably, the step S1 further includes the following steps:
preprocessing the end face of the calibration bolt, and arranging a piezoelectric sensor on the end face of the calibration bolt; the piezoelectric sensor has an ultrasonic self-generating and self-receiving mode and is used for ultrasonic measurement of the calibration bolt.
Preferably, the method further comprises the following steps before step S8:
preprocessing the end face of the bolt to be detected, and arranging a piezoelectric sensor on the end face of the bolt to be detected; the piezoelectric sensor has an ultrasonic self-transmitting and self-receiving mode and is used for ultrasonic measurement of the bolt to be measured.
Further, the end face of the calibration bolt is preprocessed, and a piezoelectric sensor is arranged on the end face of the calibration bolt, and the specific operation includes:
polishing the end face of the calibration bolt, and bonding a piezoelectric sensor on the polished end face; the piezoelectric sensor includes a piezoelectric wafer.
Further, to the terminal surface of the bolt that awaits measuring carries out the preliminary treatment, and be in the terminal surface of the bolt that awaits measuring sets up piezoelectric sensor, and concrete operation includes:
polishing the end face of the bolt to be detected, and bonding a piezoelectric sensor on the polished end face; the piezoelectric sensor includes a piezoelectric wafer.
Preferably, the piezoelectric wafer is a circular piezoelectric wafer and is concentric with the end face of the calibration bolt or concentric with the end face of the bolt to be tested.
Preferably, the polished end face satisfies a surface roughness of less than 3.2 μm.
Compared with the prior art, the invention has the beneficial effects that:
the existing ultrasonic bolt pretightening force measurement technology has the hardware sampling rate of only 200MHz due to hardware limitation, a three-spline interpolation algorithm is usually adopted for improving the pretightening force resolution, but the calculation complexity of the interpolation algorithm is about O (N ^2), when the number N of data points is large, the calculation time is raised in a quadratic way, the calculation amount is huge, and the real-time requirement cannot be met. According to the technical scheme, a time domain interpolation idea in the prior art is changed, frequency domain zero padding operation is performed, the method has the characteristic of high calculation speed, and the measuring speed of the bolt pre-tightening force can be effectively improved by measuring the bolt pre-tightening force by adopting the technical scheme provided by the invention.
Drawings
FIG. 1 is a flowchart of time delay estimation steps S4-S7 in an ultrasonic real-time bolt pretension estimation method;
FIG. 2 is a schematic diagram of a bonding structure of a piezoelectric wafer and a bolt according to an embodiment of the present invention;
in the figure, 1 is a bolt, and 2 is a piezoelectric wafer.
Detailed Description
The invention is described in further detail below with reference to the figures and specific embodiments.
As shown in fig. 1 and 2, a real-time estimation method for ultrasonic bolt pretightening force includes the steps:
s1: selecting any bolts with the same specification and the same batch of bolts to be measured as calibration bolts, and measuring the calibration bolts through ultrasonic waves under the zero-stress state to obtain a reference waveform w0
S2: performing calibration test on the calibration bolt to obtain calibration data sigmaiAnd wiWhere σ isiIndicating the magnitude of the pretension, wiIndicating the corresponding pretensionUltrasonic waveforms in the state;
s3: to w0Performing fast Fourier transform to obtain reference waveform frequency domain data FFT (w)0) And find FFT (w)0) Conjugation yields FFT (w)0)*
S4: to wiPerforming fast Fourier transform to obtain a calibration waveform frequency domain array FFT (wi);
S5: FFT (w)0) And FFT (w)0)*Element multiplication is carried out to obtain FFT (R (w)0,wi) And FFT (R (w)) is performed on the data according to the interpolation multiple M0,wi) Zero-1 times of the tail part to obtain a data set FFT (R)M(w0,wi));
S6: FFT (R) on the data setM(w0,wi) ) inverse Fourier transform to obtain cross-correlation array RM(w0,wi) Looking up RM(w0,wi) The maximum value element corresponds to the position maxPos, and w is obtained according to a time delay estimation formula delta T ═ maxPos delta T/M0,wiTime delay Δ t of waveform therebetweeniWherein Δ T is a hardware sampling time interval;
s7: fitting a bolt pretightening force formula by using the calibration data by using a least square method to obtain a bolt pretightening force coefficient A;
s8: measuring the bolt to be measured through ultrasonic waves to obtain a measured waveform wtest
S9: repeating the steps S4-S7 to obtain a measured waveform wtestAnd the reference waveform w0Ultrasonic time delay Δ t therebetweentestAnd substituting the bolt pre-tightening force coefficient A obtained by calculation in the S7 into a bolt pre-tightening force formula to obtain the pre-tightening force of the bolt to be tested.
Preferably, the expression of the bolt pretension formula is as follows:
σ=A·Δt
in the formula, A is the bolt pre-tightening force coefficient, sigma is the bolt pre-tightening force, and delta t is the ultrasonic time delay.
Preferably, the specific operation step of step S2 includes:
s2-1: determining the clamping length of the calibration bolt according to the actual fastening condition, and clamping the calibration bolt on a stretcher;
s2-2: in the elastic stage of the calibration bolt, applying different pretightening forces sigma to the calibration bolt by using a stretcheriWhere i is 1, 2, …, n (n ≧ 5), n denoting the number of different stress states applied by the stretcher;
s2-3: acquiring ultrasonic signals under corresponding pretightening force to obtain a calibration waveform wi
Preferably, the specific operation step of step S7 includes:
s7-1: will calibrate Δ t in the dataiAs an argument, σiPerforming linear fitting to obtain a linear fitting equation sigma k · Δ t + b as a dependent variable;
s7-2: judging whether the intercept b of the linear fitting equation meets the effective condition of the calibration data, if not, performing calibration data acquisition on the calibration bolt again, and updating the calibration data until the intercept b of the linear fitting equation meets the effective condition of the calibration data;
s7-3: if the intercept b meets the effective condition of the calibration data, the coefficient of the bolt pretension force is A ═ k;
preferably, the step S1 further includes the following steps:
preprocessing the end face of the calibration bolt, and arranging a piezoelectric sensor on the end face of the calibration bolt; the piezoelectric sensor has an ultrasonic self-generating and self-receiving mode and is used for ultrasonic measurement of the calibration bolt.
Preferably, the method further comprises the following steps before step S8:
preprocessing the end face of the bolt to be detected, and arranging a piezoelectric sensor on the end face of the bolt to be detected; the piezoelectric sensor has an ultrasonic self-transmitting and self-receiving mode and is used for ultrasonic measurement of the bolt to be measured.
Further, the end face of the calibration bolt is preprocessed, and a piezoelectric sensor is arranged on the end face of the calibration bolt, and the specific operation includes:
polishing the end face of the calibration bolt, and bonding a piezoelectric sensor on the polished end face; the piezoelectric sensor comprises a piezoelectric wafer 2.
Further, to the terminal surface of the bolt that awaits measuring carries out the preliminary treatment, and be in the terminal surface of the bolt that awaits measuring sets up piezoelectric sensor, and concrete operation includes:
polishing the end face of the bolt to be detected, and bonding a piezoelectric sensor on the polished end face; the piezoelectric sensor comprises a piezoelectric wafer 2.
Preferably, the piezoelectric wafer 2 is a circular piezoelectric wafer and is concentric with the end face of the calibration bolt or concentric with the end face of the bolt to be tested.
Preferably, the polished end face satisfies a surface roughness of less than 3.2 μm.
The specific embodiment is as follows:
the method is used for measuring the pretightening force of the stainless steel bolt with the diameter of 20mm, and the method is specifically described according to the following steps:
s1: selecting any stainless steel bolt with the same specification and batch as the bolt to be detected and the diameter of 20mm as a calibration bolt, exciting the piezoelectric wafer 2 to generate ultrasonic waves in a self-generating and self-receiving mode under a zero-stress state, and acquiring ultrasonic signals to obtain a reference waveform w of the calibration bolt0
The step S1 further includes the steps of: the end face of the calibration bolt is polished to meet the requirement that the surface roughness is less than 3.2 mu m, the polished end face is bonded with the piezoelectric wafer 2 by using thread fastening glue, and meanwhile, the end faces of the piezoelectric wafer 2 and the bolt 1 are concentric, and the structure of the calibration bolt is shown in figure 2.
S2: performing calibration test on the calibration bolt to obtain calibration data sigmaiAnd wiWhere σ isiIndicating the magnitude of the pretension, wiIndicating the ultrasonic waveform in the corresponding pre-tightening force state. The specific operation comprises the following steps:
s2-1: and determining the clamping length of the calibration bolt according to the actual fastening condition, and clamping the calibration bolt on a drawing machine.
S2-2: in the elastic stage of the calibration bolt, applying different pretightening forces sigma to the calibration bolt by using a stretcheriWherein i is 1, 2, 3, 4, 5.
S2-3: acquiring ultrasonic signals under corresponding pretightening force to obtain a calibration waveform wi
S3: to w0Performing fast Fourier transform to obtain reference waveform frequency domain data FFT (w)0) And find FFT (w)0) Conjugation yields FFT (w)0)*
S4: to wiPerforming fast Fourier transform to obtain a calibration waveform frequency domain array FFT (wi)。
S5: FFT (w)0) And FFT (w)0)*Element multiplication is carried out to obtain FFT (R (w)0,wi) And FFT (R (w)) at data according to interpolation factor M-200,wi) Tail-compensated 19 times zero to obtain data set FFT (R)M(w0,wi))。
S6: FFT (R) on the data setM(w0,wi) ) inverse Fourier transform to obtain cross-correlation array RM(w0,wi) Looking up RM(w0,wi) The maximum value element corresponds to the position maxPos, and w is obtained according to a time delay estimation formula delta T ═ maxPos delta T/M0,wiTime delay Δ t of waveform therebetweeniWhere Δ T is the hardware sampling time interval of 5ns, as shown in table 1.
TABLE 1 Pretightening force-time delay data
Ultrasound time delay (ns) 0 24.2 51.2 75.4 101.2 125.2
Stress (MPa) 0 29.1 61.5 90.5 121.5 150.3
S7: fitting a bolt pretightening force formula by using the calibration data by using a least square method to obtain a bolt pretightening force coefficient A;
s7-1: will calibrate Δ t in the dataiAs an argument, σiAs a dependent variable, linear fitting was performed to obtain a linear fitting equation σ of 1.2003 Δ t + 0.0248.
S7-2: the intercept b of the linear fitting equation is 0.0248, which meets the valid condition of the calibration data (i.e. | b ≦ 10 MPa).
S7-3: the pre-tightening force coefficient of the bolt is A which is 1.2003;
s8: measuring the bolt to be measured through ultrasonic waves to obtain a measured waveform wtest(ii) a Specifically, the clamping length of the bolt to be measured is determined according to the actual fastening condition, the bolt to be measured is clamped on a stretching machine, the bolt to be measured is stretched to 120MPa through the stretching machine, the piezoelectric wafer 2 is excited in a self-transmitting and self-receiving mode to generate ultrasonic waves, and ultrasonic signals are collected to obtain a measurement waveform w related to the bolt to be measuredtest
Before step S8, the method further includes the steps of: the end face of the bolt to be tested is polished to meet the requirement that the surface roughness is less than 3.2 mu m, the piezoelectric wafer 2 is bonded on the polished end face by using thread fastening glue, and meanwhile, the end faces of the piezoelectric wafer 2 and the bolt 1 are concentric, and the structure of the bolt is shown in figure 2.
S9: repeating the steps S4-S7 to obtain a measured waveform wtestAnd the reference waveform w0Ultrasonic time delay Δ t therebetweentestAnd substituting the bolt pre-tightening force coefficient A obtained by calculation in S7 into a bolt pre-tightening force formula to obtain the pre-tightening force of the bolt to be measured, and counting the pre-tightening force calculation time, wherein the result is shown in the column of 'the technical scheme' in Table 3.
TABLE 3 comparative table of measurement results
Prior art solutions The technical scheme
Calculating speed 750.51ms 98.52ms
Pretightening force estimation value 118.5 117.8
Preload measurement error (%) 1.25 1.83
As can be seen from table 3: the existing ultrasonic bolt pretightening force measurement technology has the hardware sampling rate of only 200MHz due to hardware limitation, a three-spline interpolation algorithm is usually adopted for improving the pretightening force resolution, but the calculation complexity of the interpolation algorithm is about O (N ^2), when the data point number N is large, the calculation time is raised in a quadratic way, the calculation amount is huge, and the calculation speed is about 750ms, so that the real-time requirement cannot be met. In the technical scheme of the invention, the time domain interpolation idea in the prior art is changed, and the zero padding operation of the frequency domain is carried out, so that the method has the characteristic of high calculation speed, and the calculation speed is only 98 ms; by adopting the technical scheme of the invention to measure the pretightening force of the bolt, the measuring speed of the pretightening force of the bolt can be effectively improved.
The above-described embodiments are merely preferred embodiments of the present invention, and should not be construed as limiting the present invention, and features in the embodiments and examples in the present application may be arbitrarily combined with each other without conflict. The protection scope of the present invention is defined by the claims, and includes equivalents of technical features of the claims. I.e., equivalent alterations and modifications within the scope hereof, are also intended to be within the scope of the invention.

Claims (10)

1. An ultrasonic real-time estimation method for bolt pretightening force is characterized by comprising the following steps: the method comprises the following steps:
s1: selecting any bolts with the same specification and the same batch of bolts to be measured as calibration bolts, and measuring the calibration bolts through ultrasonic waves under the zero-stress state to obtain a reference waveform w0
S2: performing calibration test on the calibration bolt to obtain calibration data sigmaiAnd wiWhere σ isiIndicating the magnitude of the pretension, wiRepresenting ultrasonic waveforms in corresponding pretightening force states;
s3: to w0Performing fast Fourier transform to obtain reference waveform frequency domain data FFT (w)0) And find FFT (w)0) Conjugation yields FFT (w)0)*
S4: to wiPerforming fast Fourier transform to obtain a calibration waveform frequency domain array FFT (wi);
S5: FFT (w)0) And FFT(w0)*Element multiplication is carried out to obtain FFT (R (w)0,wi) And FFT (R (w)) is performed on the data according to the interpolation multiple M0,wi) Zero-1 times of the tail part to obtain a data set FFT (R)M(w0,wi));
S6: FFT (R) on the data setM(w0,wi) ) inverse Fourier transform to obtain cross-correlation array RM(w0,wi) Looking up RM(w0,wi) The maximum value element corresponds to the position maxPos, and w is obtained according to a time delay estimation formula delta T ═ maxPos delta T/M0,wiTime delay Δ t of waveform therebetweeniWherein Δ T is a hardware sampling time interval;
s7: fitting a bolt pretightening force formula by using the calibration data by using a least square method to obtain a bolt pretightening force coefficient A;
s8: measuring the bolt to be measured through ultrasonic waves to obtain a measured waveform wtest
S9: repeating the steps S4-S7 to obtain a measured waveform wtestAnd the reference waveform w0Ultrasonic time delay deltar therebetweentestAnd substituting the bolt pre-tightening force coefficient A obtained by calculation in the S7 into a bolt pre-tightening force formula to obtain the pre-tightening force of the bolt to be tested.
2. The ultrasonic real-time estimation method of bolt pretightening force according to claim 1, characterized in that: the expression of the bolt pretightening force formula is as follows:
σ=A·Δt
in the formula, A is the bolt pre-tightening force coefficient, sigma is the bolt pre-tightening force, and delta t is the ultrasonic time delay.
3. The ultrasonic real-time estimation method of bolt pretightening force according to claim 1, characterized in that: the specific operation steps of step S2 include:
s2-1: determining the clamping length of the calibration bolt according to the actual fastening condition, and clamping the calibration bolt on a stretcher;
s2-2: in the above-mentionedIn the elastic stage of the calibration bolt, a stretcher is utilized to apply different pretightening forces sigma to the calibration boltiWhere i is 1, 2, …, n (n ≧ 5), n denoting the number of different stress states applied by the stretcher;
s2-3: acquiring ultrasonic signals under corresponding pretightening force to obtain a calibration waveform wi
4. The ultrasonic real-time estimation method of bolt pretightening force according to claim 1, characterized in that: the specific operation steps of step S7 include:
s7-1: will calibrate Δ t in the dataiAs an argument, σiPerforming linear fitting to obtain a linear fitting equation sigma k · Δ t + b as a dependent variable;
s7-2: judging whether the intercept b of the linear fitting equation meets the effective condition of the calibration data, if not, performing calibration data acquisition on the calibration bolt again, and updating the calibration data until the intercept b of the linear fitting equation meets the effective condition of the calibration data;
s7-3: and if the intercept b meets the effective condition of the calibration data, the bolt pretension coefficient is A-k.
5. The ultrasonic real-time estimation method of bolt pretightening force according to claim 1, characterized in that: the step S1 further includes the steps of:
preprocessing the end face of the calibration bolt, and arranging a piezoelectric sensor on the end face of the calibration bolt; the piezoelectric sensor has an ultrasonic self-generating and self-receiving mode and is used for ultrasonic measurement of the calibration bolt.
6. The ultrasonic real-time estimation method of bolt pretightening force according to claim 1, characterized in that: the method further comprises the following steps before the step S8:
preprocessing the end face of the bolt to be detected, and arranging a piezoelectric sensor on the end face of the bolt to be detected; the piezoelectric sensor has an ultrasonic self-transmitting and self-receiving mode and is used for ultrasonic measurement of the bolt to be measured.
7. The ultrasonic real-time estimation method of bolt pretension force according to claim 5, characterized in that: the end face of the calibration bolt is preprocessed, a piezoelectric sensor is arranged on the end face of the calibration bolt, and the specific operation comprises the following steps:
polishing the end face of the calibration bolt, and bonding a piezoelectric sensor on the polished end face; the piezoelectric sensor comprises a piezoelectric wafer (2).
8. The ultrasonic real-time estimation method of bolt pretension force according to claim 6, characterized in that: the end face of the bolt to be detected is preprocessed, a piezoelectric sensor is arranged on the end face of the bolt to be detected, and the concrete operation comprises the following steps:
polishing the end face of the bolt to be detected, and bonding a piezoelectric sensor on the polished end face; the piezoelectric sensor comprises a piezoelectric wafer (2).
9. The real-time estimation method of ultrasonic bolt pretension according to claim 7 or 8, characterized in that: the piezoelectric wafer (2) is a circular piezoelectric wafer and is concentric with the end face of the calibration bolt or concentric with the end face of the bolt to be tested.
10. The real-time estimation method of ultrasonic bolt pretension according to claim 7 or 8, characterized in that: the polished end face satisfies that the surface roughness is less than 3.2 mu m.
CN202110420543.9A 2021-04-19 2021-04-19 Ultrasonic real-time estimation method for bolt pretightening force Active CN113218560B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110420543.9A CN113218560B (en) 2021-04-19 2021-04-19 Ultrasonic real-time estimation method for bolt pretightening force

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110420543.9A CN113218560B (en) 2021-04-19 2021-04-19 Ultrasonic real-time estimation method for bolt pretightening force

Publications (2)

Publication Number Publication Date
CN113218560A true CN113218560A (en) 2021-08-06
CN113218560B CN113218560B (en) 2022-05-17

Family

ID=77087875

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110420543.9A Active CN113218560B (en) 2021-04-19 2021-04-19 Ultrasonic real-time estimation method for bolt pretightening force

Country Status (1)

Country Link
CN (1) CN113218560B (en)

Citations (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2708102A1 (en) * 1993-07-22 1995-01-27 Centre Techn Ind Mecanique Workpiece tightening control device by ultrasound and signal sampling.
JPH07255721A (en) * 1994-03-18 1995-10-09 Fujitsu Ltd Ultrasonic diagnostic device
JPH08322836A (en) * 1995-06-01 1996-12-10 Hitachi Medical Corp Ultrasonic diagnostic device
US20060137463A1 (en) * 2004-12-23 2006-06-29 Sayed Nassar Ultrasonic control of bolt tightening
CN101644773A (en) * 2009-03-20 2010-02-10 中国科学院声学研究所 Real-time frequency domain super-resolution direction estimation method and device
CN102608602A (en) * 2012-03-13 2012-07-25 北京航空航天大学 Ultralow sidelobe synthetic aperture radar imaging method based on complete complementary sequence
CN104605888A (en) * 2015-01-22 2015-05-13 飞依诺科技(苏州)有限公司 Ultrasonic signal time delay estimation method and system
CN107105497A (en) * 2017-06-07 2017-08-29 重庆重邮汇测通信技术有限公司 Master sync signal detection method in a kind of LTE system
CN108387338A (en) * 2018-02-07 2018-08-10 大连理工大学 A kind of bolt pretightening real time high precision detection method and system based on piezoelectric supersonic chip
CN108896230A (en) * 2018-05-11 2018-11-27 北京理工大学 It is a kind of that method is determined based on the bolt clipping forcee ultrasound detection of finite element and crucial detection parameters
CN109668672A (en) * 2019-01-08 2019-04-23 中国大唐集团新能源科学技术研究院有限公司 Ultrasonic bolt pre-tightens force measuring method
CN110501423A (en) * 2019-08-15 2019-11-26 重庆大学 A kind of high-resolution minimum variance ultrasonic imaging method based on frequency domain segmentation
CN110941980A (en) * 2019-07-16 2020-03-31 上海师范大学 Multipath time delay estimation method and device based on compressed sensing in dense environment
JP6773878B1 (en) * 2019-12-19 2020-10-21 中日本ハイウェイ・エンジニアリング東京株式会社 Concrete structure internal condition inspection method and system used for that method
CN111968170A (en) * 2020-08-26 2020-11-20 广东工业大学 Online binocular vision distance measurement method based on cross-correlation time delay estimation
CN112033597A (en) * 2020-08-25 2020-12-04 苏州博昇科技有限公司 High-strength bolt axial force ultrasonic double-wave rapid calibration and detection method
CN112179553A (en) * 2020-09-09 2021-01-05 西南交通大学 Method for ultrasonically and synchronously measuring axial force and shearing force of bolt
CN112461429A (en) * 2020-11-10 2021-03-09 西南交通大学 Ultrasonic pretightening force measurement method for low-elasticity-modulus material bolt
US20210084424A1 (en) * 2019-09-16 2021-03-18 Gaudio Lab, Inc. Method for generating binaural signals from stereo signals using upmixing binauralization, and apparatus therefor

Patent Citations (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2708102A1 (en) * 1993-07-22 1995-01-27 Centre Techn Ind Mecanique Workpiece tightening control device by ultrasound and signal sampling.
JPH07255721A (en) * 1994-03-18 1995-10-09 Fujitsu Ltd Ultrasonic diagnostic device
JPH08322836A (en) * 1995-06-01 1996-12-10 Hitachi Medical Corp Ultrasonic diagnostic device
US20060137463A1 (en) * 2004-12-23 2006-06-29 Sayed Nassar Ultrasonic control of bolt tightening
CN101644773A (en) * 2009-03-20 2010-02-10 中国科学院声学研究所 Real-time frequency domain super-resolution direction estimation method and device
CN102608602A (en) * 2012-03-13 2012-07-25 北京航空航天大学 Ultralow sidelobe synthetic aperture radar imaging method based on complete complementary sequence
CN104605888A (en) * 2015-01-22 2015-05-13 飞依诺科技(苏州)有限公司 Ultrasonic signal time delay estimation method and system
CN107105497A (en) * 2017-06-07 2017-08-29 重庆重邮汇测通信技术有限公司 Master sync signal detection method in a kind of LTE system
CN108387338A (en) * 2018-02-07 2018-08-10 大连理工大学 A kind of bolt pretightening real time high precision detection method and system based on piezoelectric supersonic chip
CN108896230A (en) * 2018-05-11 2018-11-27 北京理工大学 It is a kind of that method is determined based on the bolt clipping forcee ultrasound detection of finite element and crucial detection parameters
CN109668672A (en) * 2019-01-08 2019-04-23 中国大唐集团新能源科学技术研究院有限公司 Ultrasonic bolt pre-tightens force measuring method
CN110941980A (en) * 2019-07-16 2020-03-31 上海师范大学 Multipath time delay estimation method and device based on compressed sensing in dense environment
CN110501423A (en) * 2019-08-15 2019-11-26 重庆大学 A kind of high-resolution minimum variance ultrasonic imaging method based on frequency domain segmentation
US20210084424A1 (en) * 2019-09-16 2021-03-18 Gaudio Lab, Inc. Method for generating binaural signals from stereo signals using upmixing binauralization, and apparatus therefor
JP6773878B1 (en) * 2019-12-19 2020-10-21 中日本ハイウェイ・エンジニアリング東京株式会社 Concrete structure internal condition inspection method and system used for that method
CN112033597A (en) * 2020-08-25 2020-12-04 苏州博昇科技有限公司 High-strength bolt axial force ultrasonic double-wave rapid calibration and detection method
CN111968170A (en) * 2020-08-26 2020-11-20 广东工业大学 Online binocular vision distance measurement method based on cross-correlation time delay estimation
CN112179553A (en) * 2020-09-09 2021-01-05 西南交通大学 Method for ultrasonically and synchronously measuring axial force and shearing force of bolt
CN112461429A (en) * 2020-11-10 2021-03-09 西南交通大学 Ultrasonic pretightening force measurement method for low-elasticity-modulus material bolt

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
JOHAN E. CARLSON 等: "Measurement of the clamping force applied by load-bearing bolts using a combination of compression and shear ultrasonic waves", 《 2015 IEEE INTERNATIONAL ULTRASONICS SYMPOSIUM (IUS)》 *
KIM N 等: "Measurement of axial stress using mode-converted ultrasound", 《NDT & E INTERNATIONAL》 *
方汉方 等: "超声波传输时间的高精度测量", 《信号处理》 *
朱士明: "声时差法测定螺栓的预紧应力", 《声学技术》 *

Also Published As

Publication number Publication date
CN113218560B (en) 2022-05-17

Similar Documents

Publication Publication Date Title
CN108387338B (en) Real-time high-precision detection method and system for bolt pretightening force based on piezoelectric ultrasonic wafer
CN112461429B (en) Ultrasonic pretightening force measurement method for low-elasticity-modulus material bolt
CN108896230B (en) Bolt fastening force ultrasonic detection and key detection parameter determination method based on finite element
CN103148815B (en) Based on the thickness of thin layer supersonic detection method of sound pressure reflection coefficient autocorrelation function
CN109959477A (en) A kind of GIS disc insulator epoxy test block internal stress ultrasonic longitudinal wave detection method and system
CN114487109B (en) Non-baseline data stress online monitoring method, system, equipment and medium based on single-mode multi-frequency signal fusion
CN111307357A (en) Differential method-based ultrasonic detection method for bolt fastening axial force
CN112179553B (en) Method for ultrasonically and synchronously measuring axial force and shearing force of bolt
CN111442977B (en) Method for determining maximum stress of circular film under action of transversely uniformly distributed load
CN107192492B (en) A kind of tools for bolts ' pretension torque supersonic guide-wave monitoring method based on improvement time reversal
CN107167273B (en) High-strength bolt connecting node plate compaction degree detection method based on ultrasonic echo
CN111442984B (en) Method for determining maximum stress of circular film under transversely uniform load
CN111337171A (en) Acoustic time difference measuring method applied to critical refraction longitudinal wave stress detection
CN111442982A (en) Method for determining maximum stress of circular film under uniformly distributed load
CN101949894A (en) Method for detecting interface contact strength by double frequency ultrasound
CN110646119B (en) Method for measuring surface stress tensor of rolled metal material by ultrasonic wave
CN111442985A (en) Method for determining maximum deflection of circular film under transversely uniformly distributed load
CN114994175A (en) Space coupling ultrasonic stress detection device and method for modal decomposition double-spectrum analysis
CN107121335A (en) A kind of dynamic indentation test method of unicast guide rod material
CN113218560B (en) Ultrasonic real-time estimation method for bolt pretightening force
CN103616436B (en) A kind of high-precision ultrasound detection method of contact stiffness
CN111141206A (en) Strain gauge dynamic characteristic detection device and test method thereof
CN113624395A (en) Bolt axial force detection coefficient calibration method based on ultrasound
CN109738354B (en) Corrosion on-line monitoring device
CN111474039A (en) Method for determining maximum deflection of prestressed circular film under action of uniformly distributed load

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