Hammer impulsive force measuring method and system suitable for strain pile power detection

Abstract

The invention relates to the technical field of hammer impulsive force measurement, in particular to a hammer impulsive force measurement method and system suitable for strain pile power detection, and the method comprises the steps: recognizing a target signal interval indicating a main energy interval in an actually measured impact signal; analyzing the signal amplitude of each moment in the target signal interval to obtain time domain interference data; a first energy center frequency and a second energy center frequency are obtained, the first energy center frequency is used for representing the center of signal energy distribution corresponding to the actually measured impact signal in the frequency domain, and the second energy center frequency is used for representing the center of signal energy distribution corresponding to the ideal impact signal in the frequency domain; analyzing the difference between the first energy center frequency and the second energy center frequency to obtain frequency domain interference data; and according to the time domain interference data and the frequency domain interference data, performing interference compensation on the actually measured impact signal to obtain a target impact signal. According to the invention, the data accuracy of hammer impulsive force measurement signals can be improved.

Description

Technical Field

[0001] The present invention relates to the technical field of hammer impact force measurement, and in particular to a hammer impact force measurement method and system suitable for dynamic detection of strain piles. Background Art

[0002] Strain pile dynamic testing is a testing technology that applies a transient impact force to the pile top and uses strain sensors installed in the pile body to collect stress wave signals. Based on the reflection and transmission characteristics of the stress wave propagating through the pile body due to the uniformity of the pile material, cross-sectional variations, or defects, the technology analyzes and determines performance parameters such as the pile body's integrity and bearing capacity. The hammer impact force refers to the transient impact force generated by the hammer striking the pile top during testing. This force acts as an excitation source, stimulating stress waves within the pile body. The waveform characteristics of these stress waves directly affect the propagation distance and detection accuracy, so precise measurement is required.

[0003] The existing patent, "A Method and Apparatus for Measuring Hammer Force in High-Strain Dynamic Testing of Foundation Pile Foundations" (CN101078661B), discloses: Hammer force compensation is achieved through simultaneous measurement with dual accelerometers. The method first acquires the acceleration signal of the falling hammer and performs a low-pass filter. The hammer force at the pile top is calculated according to Newton's second law. The pile body acceleration is then measured and the inertial force at the pile head is calculated. This inertial effect correction is used to determine the hammer force at the measurement section below the pile top. This method directly acquires force signals through hardware measurement, avoiding the errors associated with traditional strain gauges that rely on material parameter estimation. It is particularly suitable for high-strain testing of concrete piles.

[0004] However, the existing technology does not take into account that the force diversion effect caused by sensor installation deviation will change the spectral distribution of the impact load, and the pile body damping and geometric dispersion effects will further introduce waveform distortion during the propagation process, and the two will form dynamic coupling in the frequency domain and time domain, so that the measured force signal contains the composite interference of sensor installation deviation and propagation path interference.

[0005] That is to say, in the prior art, the data accuracy of the hammer impact measurement signal used in the dynamic testing of strain piles is relatively poor. Summary of the Invention

[0006] In order to solve the technical problem of poor data accuracy of hammer impact measurement signals used in the prior art for strain pile dynamic testing, the present invention aims to provide a hammer impact measurement method and system suitable for strain pile dynamic testing. The technical solutions adopted are as follows: In a first aspect, an embodiment of the present invention provides a hammer impact force measurement method suitable for dynamic testing of strain piles, the method comprising: Identifying a target signal interval in a measured impact signal, wherein the measured impact signal is a signal generated when a hammer impacts a pile under test, and the target signal interval is used to indicate a main energy interval of the measured impact signal; Analyzing the signal amplitude at each moment in the target signal interval to obtain time domain interference data, where the time domain interference data is used to indicate a severity of a change in the signal amplitude in the target signal interval; Obtaining a first energy center frequency and a second energy center frequency, wherein the first energy center frequency is used to represent the center of the signal energy distribution corresponding to the measured impulse signal in the frequency domain, and the second energy center frequency is used to represent the center of the signal energy distribution corresponding to the reference impulse signal in the frequency domain, and the reference impulse signal is an ideal impulse signal corresponding to the measured impulse signal; Analyzing a difference between the first energy center frequency and the second energy center frequency to obtain frequency domain interference data; According to the time domain interference data and the frequency domain interference data, interference compensation is performed on the measured impact signal to obtain a target impact signal, which is an impact signal used for dynamic detection of the inspected pile.

[0007] Furthermore, identifying the target signal interval in the measured impact signal includes: Based on a preset sliding window length and sliding step size, the measured impact signal is divided into a plurality of signal windows; Determine the maximum window signal energy among the multiple window signal energies respectively corresponding to the multiple signal windows as the signal energy threshold, wherein the window signal energy is the square sum of the multiple signal energies included in the corresponding signal window; Determine the product of the signal energy threshold and a preset proportional coefficient as a standard threshold; The target signal interval is identified in the measured impulse signal according to the standard threshold and the multiple signal windows, wherein the average of multiple window signal energies corresponding to the multiple signal windows covered by the target signal interval is greater than the standard threshold.

[0008] Furthermore, identifying the target signal interval in the measured impulse signal according to the standard threshold and the multiple signal windows includes: Identifying a first signal window and a second signal window in the multiple signal windows based on the standard threshold, wherein the first signal window is a first signal window among a plurality of consecutive signal windows whose corresponding window signal energy is greater than the standard threshold, and the second signal window is a signal window located in a preset order among a plurality of consecutive signal windows whose corresponding window signal energy is less than the standard threshold; The target signal interval is identified in the measured impulse signal by taking the window start point of the first signal window as the interval start point and the window end point of the second signal window as the interval end point.

[0009] Furthermore, analyzing the signal amplitude at each moment in the target signal interval to obtain time domain interference data includes: Obtaining a signal slope at each moment in the target signal interval based on the signal amplitude at each moment in the target signal interval, wherein the signal slope is the instantaneous rate of change of the signal amplitude at the corresponding moment; Obtaining a slope mutation value at each moment in the target signal interval according to the signal slope at each moment in the target signal interval, wherein the slope mutation value is the absolute value of the second-order derivative of the signal slope at the corresponding moment; The sum of the slope mutation values at each moment in the target signal interval is calculated, and the sum is determined as the time domain interference data.

[0010] Furthermore, obtaining the first energy center frequency and the second energy center frequency includes: Performing frequency domain conversion on the time domain data to be processed to obtain a frequency domain amplitude spectrum, wherein the time domain data to be processed is the target signal interval or the reference impulse signal; Determining the product of the square of the amplitude of each frequency component in the frequency domain amplitude spectrum and each frequency component in the frequency domain amplitude spectrum as a first index value of each frequency component in the frequency domain amplitude spectrum; and determining the square of the amplitude of each frequency component in the frequency domain amplitude spectrum as a second index value of each frequency component in the frequency domain amplitude spectrum; Performing a sum calculation on the first index value of each frequency component in the frequency domain amplitude spectrum to obtain a first sum value, and performing a sum calculation on the second index value of each frequency component in the frequency domain amplitude spectrum to obtain a second sum value; When the time domain data to be processed is the target signal interval, the ratio of the first sum value to the second sum value is determined as the first energy center frequency; when the time domain data to be processed is the reference impulse signal, the ratio of the first sum value to the second sum value is determined as the second energy center frequency.

[0011] Furthermore, before obtaining the first energy center frequency and the second energy center frequency, the method further includes: Identifying a first frequency segment greater than a frequency threshold in a measured amplitude spectrum, wherein the measured amplitude spectrum is an amplitude spectrum formed after frequency domain conversion of the target signal interval; performing amplitude correction on each frequency component in the first frequency band to obtain a corrected amplitude of each frequency component in the first frequency band, wherein the corrected amplitude is the product of the amplitude of the corresponding frequency component in the measured amplitude spectrum and a correction parameter, wherein the correction parameter is determined based on the frequency value of the corresponding frequency component and the time-domain interference data, wherein the time-domain interference data and the correction parameter are positively correlated, and the correction parameter and the frequency value of the corresponding frequency component are positively correlated; Performing cubic spline interpolation according to the corrected amplitude of each frequency component in the first frequency band to obtain a target frequency band curve; The target frequency band curve and the second frequency band in the measured amplitude spectrum that is less than the frequency threshold are combined to obtain a target amplitude spectrum, which is used to represent the amplitude spectrum obtained by frequency domain conversion of the reference impulse signal.

[0012] Furthermore, analyzing the difference between the first energy center frequency and the second energy center frequency to obtain frequency domain interference data includes: Calculating an absolute value of a difference between the first energy center frequency and the second energy center frequency; The frequency domain interference data is obtained according to the ratio of the absolute value and the second energy center frequency.

[0013] Furthermore, performing interference compensation on the measured impulse signal according to the time domain interference data and the frequency domain interference data to obtain a target impulse signal includes: Based on a pre-constructed inverse filter, inverse filtering is performed on a frequency domain signal corresponding to the measured impulse signal to obtain a compensated frequency domain signal, wherein the transfer function corresponding to the inverse filter includes a high-frequency gain compensation coefficient and a time domain distortion penalty coefficient, the high-frequency gain compensation coefficient is determined based on the frequency domain interference data and a frequency value of a corresponding frequency component, and the time domain distortion penalty coefficient is determined based on the time domain interference data; Performing a time domain transformation on the compensated frequency domain signal to obtain the target impact signal.

[0014] In a second aspect, another embodiment of the present invention provides a hammer impact force measurement system suitable for dynamic testing of strain piles, the system comprising: an interval identification module, configured to identify a target signal interval in a measured impact signal, wherein the measured impact signal is a signal generated when a hammer impacts a pile under test, and the target signal interval is configured to indicate a main energy interval of the measured impact signal; a time domain analysis module, configured to analyze the signal amplitude at each moment in the target signal interval to obtain time domain interference data, wherein the time domain interference data is used to indicate a severity of a change in the signal amplitude in the target signal interval; a frequency center acquisition module, configured to acquire a first energy center frequency and a second energy center frequency, wherein the first energy center frequency is used to represent the center of the signal energy distribution corresponding to the measured impulse signal in the frequency domain, and the second energy center frequency is used to represent the center of the signal energy distribution corresponding to the reference impulse signal in the frequency domain, and the reference impulse signal is an ideal impulse signal corresponding to the measured impulse signal; a frequency domain analysis module, configured to analyze a difference between the first energy center frequency and the second energy center frequency to obtain frequency domain interference data; The interference compensation module is used to perform interference compensation on the measured impact signal according to the time domain interference data and the frequency domain interference data to obtain a target impact signal, where the target impact signal is an impact signal used for dynamic detection of the strain pile.

[0015] In a third aspect, another embodiment of the present invention further provides an electronic device, comprising a processor, a memory, and a computer program stored in the memory and executable on the processor, wherein the computer program implements the steps of the method described in the first aspect when executed by the processor.

[0016] In a fourth aspect, another embodiment of the present invention further provides a computer-readable storage medium, on which a computer program is stored, and when the computer program is executed by a processor, the steps of the method described in the first aspect are implemented.

[0017] The present invention has the following beneficial effects: The present invention identifies a target signal interval corresponding to a main energy interval in a measured impact signal and analyzes the signal amplitude at each moment in the target signal interval to obtain time domain interference data, thereby quantifying the severity of the change in the signal amplitude in the target signal interval, that is, quantifying the degree of signal distortion of the measured impact signal in the time domain. The present invention then analyzes the center of the signal energy distribution corresponding to the measured impact signal in the frequency domain and compares the center of the signal energy distribution corresponding to the reference impact signal in the frequency domain to obtain frequency domain interference data, thereby quantifying the degree of signal distortion of the measured impact signal in the time domain. Finally, interference compensation is performed on the measured impact signal based on the time domain interference data and the frequency domain interference data, thereby achieving interference compensation for the measured impact signal in both the time domain and the frequency domain, thereby minimizing the adverse effects of the combined interference introduced by sensor installation deviation and propagation path interference, thereby obtaining a more accurate impact signal for dynamic detection of the inspected pile and making the detection results of the dynamic detection operation of the inspected pile based on the data more accurate and reliable. BRIEF DESCRIPTION OF THE DRAWINGS

[0018] In order to more clearly illustrate the technical solutions and advantages of the embodiments of the present invention or the prior art, the following briefly introduces the drawings required for use in the embodiments or the prior art descriptions. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0019] Figure 1 A schematic flow chart of a hammer impact force measurement method suitable for strain pile dynamic testing provided by one embodiment of the present invention; Figure 2 A schematic structural diagram of a hammer impact force measurement system suitable for strain pile dynamic testing provided by one embodiment of the present invention; Figure 3 The present invention provides a schematic structural diagram of an electronic device according to an embodiment of the present invention. DETAILED DESCRIPTION

[0020] To further illustrate the technical means and effectiveness of the present invention in achieving its intended purpose, the following, in conjunction with the accompanying drawings and preferred embodiments, details the specific implementation, structure, features, and effectiveness of a hammer impact force measurement method and system suitable for dynamic testing of strain piles, as proposed by the present invention. In the following description, references to different "one embodiment" or "another embodiment" do not necessarily refer to the same embodiment. Furthermore, specific features, structures, or characteristics of one or more embodiments may be combined in any suitable manner.

[0021] Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs.

[0022] The following describes in detail a hammer impact force measurement method and system applicable to strain pile dynamic testing provided by the present invention with reference to the accompanying drawings.

[0023] The present invention proposes a hammer impact force measurement method suitable for strain pile dynamic testing. Figure 1 , which shows a schematic flow chart of a hammer impact force measurement method suitable for strain pile dynamic testing provided by one embodiment of the present invention, the method comprising the following steps: Step S1: Identify the target signal interval in the measured impact signal.

[0024] The measured impact signal is a signal generated when the hammer impacts the pile under test, and the target signal interval is used to indicate a main energy interval of the measured impact signal.

[0025] For example, two force sensors can be symmetrically installed on the side surface of the pile under test, and the distance between the installation position of the force sensor on the pile under test and the top of the pile under test is greater than or equal to twice the diameter of the pile under test. In this example, the sampling frequency of the force sensor is 40kHz. The symmetrical installation of the two force sensors can be understood as: the two force sensors are symmetrically arranged with the axis of the pile under test as the center of symmetry, or the two force sensors are installed at the same height and the distance between the two force sensors is a preset distance (such as 30 cm). The above-mentioned measured impact signal is the signal collected by the above-mentioned force sensor when the hammer impacts the pile under test.

[0026] Hammering (referring to the hammer head hitting the top of the pile under inspection) is a transient impact load. The energy is released in a concentrated manner in the form of stress waves at the moment of impact, forming an area in the time domain signal with a steep increase in amplitude and extremely high energy density. Therefore, in one example, the time domain position of the maximum signal energy (with the amplitude of the signal as an indicator of the signal energy) can be identified in the measured impact signal, and the time domain position of the maximum signal energy is used as the center position of the target signal interval. Based on a preset interval radius (such as 10ms, the interval bandwidth of the target signal interval is 20ms at this time), the target signal interval is determined in the measured impact signal. The aforementioned interval radius is used to indicate: the distance between the center position of the target signal interval and the boundary position of the target signal interval.

[0027] Step S2: Analyze the signal amplitude at each moment in the target signal interval to obtain time domain interference data.

[0028] The time domain interference data is used to indicate the severity of the change in the signal amplitude in the target signal interval.

[0029] In the present invention, the severity of the change in the signal amplitude in the target signal interval indicated by the time domain interference data is positively correlated with the interference intensity received by the signal in the time domain when the hammer impacts the tested pile. In applications, situations that lead to inaccurate measured impact signals include: the uneven surface of the pile head of the inspected pile causes force diversion when the hammer impulse passes through the pile head; the deviation in the installation position of the force sensor causes distortion of the detected hammer impulse; the energy dissipation caused by the internal friction and viscoelastic properties of the inspected pile itself when the hammer impulse passes through the pile body of the inspected pile; the waveform distortion caused by the geometric characteristics of the pile body of the inspected pile (the propagation speed of stress waves of different frequencies in the pile body is inconsistent), etc.

[0030] In one example, multiple reference curves indicating different degrees of severity can be pre-constructed, and the curve similarity between each reference curve and the curve corresponding to the target signal interval can be calculated, and the severity indicated by the reference curve corresponding to the largest curve similarity can be determined as the severity of the change in signal amplitude in the target signal interval. Among them, multiple reference curves and the severity indicated by each reference curve can be obtained by experts through manual labeling, and the curve similarity can be understood as the mean squared error (MSE) between the reference curve and the curve corresponding to the target signal interval.

[0031] Step S3: Obtain the first energy center frequency and the second energy center frequency.

[0032] Among them, the first energy center frequency is used to represent the center of the signal energy distribution corresponding to the measured impact signal in the frequency domain, and the second energy center frequency is used to represent the center of the signal energy distribution corresponding to the reference impact signal in the frequency domain. The reference impact signal is the ideal impact signal corresponding to the measured impact signal.

[0033] The center of the corresponding signal energy distribution in the frequency domain may be: the spectral centroid of the corresponding signal energy distribution in the frequency domain, or the median of the corresponding signal energy distribution in the frequency domain, or the highest frequency of the corresponding signal energy distribution in the frequency domain.

[0034] Exemplarily, the above-mentioned reference impact signal can be obtained through simulation experiments, or by extrapolating the mid- and high-frequency band parts (which are less affected by interference than the low-frequency band parts) in the frequency domain signal corresponding to the measured impact signal, wherein the signal part below 10kHz can be determined as the low-frequency band part, and the signal part above 10kHz can be determined as the high-frequency band part.

[0035] Step S4: Analyze the difference between the first energy center frequency and the second energy center frequency to obtain frequency domain interference data.

[0036] Furthermore, analyzing the difference between the first energy center frequency and the second energy center frequency to obtain frequency domain interference data includes: Calculating an absolute value of a difference between the first energy center frequency and the second energy center frequency; The frequency domain interference data is obtained according to the ratio of the absolute value and the second energy center frequency.

[0037] If the first energy center frequency is defined as , the second energy center frequency is , then the calculation formula of frequency domain interference data E can be expressed as:

[0038] By analyzing the difference between the energy center frequencies, the degree of deviation of the measured impact signal relative to the ideal impact signal in the frequency domain can be quantified, that is, the degree of deviation of the measured impact signal in the frequency domain due to interference can be determined.

[0039] Step S5: performing interference compensation on the measured impulse signal according to the time domain interference data and the frequency domain interference data to obtain a target impulse signal.

[0040] The above-mentioned interference compensation is used to eliminate the disturbances suffered by the measured impulse signal in the time domain and the frequency domain, and restore the measured impulse signal to a real impulse signal.

[0041] The target impact signal is an impact signal used to perform dynamic detection on the inspected pile.

[0042] The present invention identifies a target signal interval corresponding to a main energy interval in a measured impact signal and analyzes the signal amplitude at each moment in the target signal interval to obtain time domain interference data, thereby quantifying the severity of the change in the signal amplitude in the target signal interval, that is, quantifying the degree of signal distortion of the measured impact signal in the time domain. The center of the signal energy distribution corresponding to the measured impact signal in the frequency domain is then analyzed and compared with the center of the signal energy distribution corresponding to the reference impact signal in the frequency domain to obtain frequency domain interference data, thereby quantifying the degree of signal distortion of the measured impact signal in the time domain. Finally, interference compensation is performed on the measured impact signal based on the time domain interference data and the frequency domain interference data. Interference compensation is achieved in both the time domain and the frequency domain, thereby minimizing the adverse effects of the combined interference introduced by sensor installation deviation and propagation path interference. A more accurate impact signal for dynamic detection of a strain pile is obtained, and the detection results of the strain pile dynamic detection operation performed based on the interference data are more accurate and reliable.

[0043] In one embodiment, identifying the target signal interval in the measured impulse signal includes: Based on a preset sliding window length and sliding step size, the measured impact signal is divided into a plurality of signal windows; Determine the maximum window signal energy among the multiple window signal energies respectively corresponding to the multiple signal windows as the signal energy threshold, wherein the window signal energy is the square sum of the multiple signal energies included in the corresponding signal window; Determine the product of the signal energy threshold and a preset proportional coefficient as a standard threshold; The target signal interval is identified in the measured impulse signal according to the standard threshold and the multiple signal windows, wherein the average of multiple window signal energies corresponding to the multiple signal windows covered by the target signal interval is greater than the standard threshold.

[0044] In some embodiments, based on a preset sliding window length and sliding step size, before the measured impact signal is divided into multiple signal windows, the measured impact signal can be first preprocessed to filter out the interference parts that are easy to process, thereby ensuring the data accuracy of the subsequently extracted time domain interference data and frequency domain interference data, and further ensuring the accuracy of the interference compensation operation based on the time domain interference data and frequency domain interference data, so that the data accuracy of the target impact signal is further improved.

[0045] The signal preprocessing operations may include: removing DC offset, signal filtering, etc.

[0046] The aforementioned DC offset removal eliminates the DC component caused by sensor zero drift or environmental noise. Signal filtering can use a 50-20 kHz bandpass filter (such as a fourth-order Butterworth filter) to filter out power frequency interference (50 Hz) and sensor high-frequency noise (>20 kHz), retaining the effective frequency band (within 20 kHz) in the measured impact signal that reflects the main energy range.

[0047] The acquisition process of the above multiple signal windows can be: Using the signal starting point of the measured impact signal as the window starting point of the first signal window, and dividing the signal portion from the window starting point of the first signal window to a preset sliding window length into a first signal window; Determine the sum of the window starting point and the sliding step length of the first signal window as the window starting point of the second signal window, and divide the signal portion from the window starting point of the first signal window to the preset sliding window length into the second signal window; The same process is repeated..., the sum of the window starting point and the sliding step of the nth signal window is determined as the window starting point of the n+1th signal window, and the signal part from the window starting point of the n+1th signal window to the preset sliding window length is divided into the n+1th signal window, until the window end point of the determined signal window coincides with the signal end point of the measured impact signal, and the window sliding operation is terminated.

[0048] It should be noted that if the signal portion from the window start point of the previous signal window to the preset sliding window length exceeds the signal end point of the measured impact signal, the signal portion between the window start point of the previous signal window and the signal end point of the measured impact signal will be used as the current signal window, and the signal window will be determined as the last signal window.

[0049] The sliding window length may be 1 ms (the range may be 0.5 to 2 ms), and the sliding step may be 0.5 ms.

[0050] It should be understood that each signal window includes multiple signal points (for example, when the sampling frequency of the force sensor is 40kHz and the sliding window length is 1ms, a signal window includes 40 signal points), and each signal point corresponds to a signal energy (which can be the amplitude of the corresponding signal point). Therefore, each signal window corresponds to multiple signal energies. By calculating the sum of the squares of the multiple signal energies corresponding to each signal window and using it as the window signal energy of the corresponding signal window, the strength of the signal energy reflected by the corresponding signal window can be quantitatively represented.

[0051] After that, the maximum window signal energy is selected as the signal energy threshold to identify the location of the signal peak in the actual impact signal, and the product of the signal energy threshold and the preset proportional coefficient is determined as the standard threshold. Finally, the target signal interval is determined by marking the threshold, so that the determined target signal interval can accurately and completely reflect the main energy interval of the aforementioned measured impact signal.

[0052] The value range of the proportional coefficient may be 0.1 to 0.3.

[0053] In one example, among the multiple signal windows, all signal windows whose corresponding window signal energy is greater than the standard threshold can be determined as candidate windows; and among the multiple candidate windows, the signal parts corresponding to several consecutive and longest candidate windows in the measured impact signal can be determined as the target signal interval.

[0054] In one embodiment, identifying the target signal interval in the measured impulse signal according to the standard threshold and the multiple signal windows includes: Identifying a first signal window and a second signal window in the multiple signal windows based on the standard threshold, wherein the first signal window is a first signal window among a plurality of consecutive signal windows whose corresponding window signal energy is greater than the standard threshold, and the second signal window is a signal window located in a preset order among a plurality of consecutive signal windows whose corresponding window signal energy is less than the standard threshold; The target signal interval is identified in the measured impulse signal by taking the window start point of the first signal window as the interval start point and the window end point of the second signal window as the interval end point.

[0055] For example, starting from the signal starting point of the measured impact signal, when the window signal energy of two consecutive signal windows is greater than the standard threshold, the window starting point of the first signal window in the two consecutive signal windows is taken as the interval starting point.

[0056] Exemplarily, the preset order may be the fifth order, that is, the fifth signal window among a plurality of consecutive signal windows whose corresponding window signal energy is less than the standard threshold is taken as the second signal window.

[0057] In this embodiment, a comparison is made based on the standard threshold and the window signal energy to identify the signal window that may correspond to the main energy interval, and the condition of window continuity is used to identify the signal window where the starting point of the target signal interval is located (i.e., the first signal window), and the signal window where the end point of the target signal interval is located (i.e., the second signal window). This not only avoids the misjudgment caused by a single noise point (mainly the misjudgment of the starting point of the target signal interval), ensuring the accuracy of the identified target signal interval, but also avoids the omission of the signal interval that originally corresponds to the main energy interval (mainly the omission of the last few signal windows in the several continuous signal windows corresponding to the target signal interval), ensuring the integrity of the identified target signal interval.

[0058] In one embodiment, analyzing the signal amplitude at each moment in the target signal interval to obtain time domain interference data includes: Obtaining a signal slope at each moment in the target signal interval based on the signal amplitude at each moment in the target signal interval, wherein the signal slope is the instantaneous rate of change of the signal amplitude at the corresponding moment; Obtaining a slope mutation value at each moment in the target signal interval according to the signal slope at each moment in the target signal interval, wherein the slope mutation value is the absolute value of the second-order derivative of the signal slope at the corresponding moment; The sum of the slope mutation values at each moment in the target signal interval is calculated, and the sum is determined as the time domain interference data.

[0059] The force shunt effect caused by installation errors in the force sensor can cause localized stress concentration or divergence in the measured impact signal. As the stress wave propagates through the pile, damping and geometric dispersion can cause the waveform's leading edge to broaden and its slope to abruptly change (for example, a step signal becomes a ramp). These conditions can lead to unnatural, sudden changes in the time-domain slope (the rate of change of the signal's amplitude in the time domain) of the force signal (i.e., the measured impact signal). The degree of this change is directly related to the installation error and the pile's characteristics.

[0060] Based on the above settings, the third-order derivative is used to extract the degree of time domain signal mutation reflected by each signal point, and the sum is calculated as the time domain interference data to comprehensively collect the time domain signal mutations existing in various locations in the target signal interval, thereby obtaining a global measurement parameter (that is, time domain interference data) that can accurately represent the overall distortion degree of the force signal in the time domain.

[0061] For example, the calculation formula of the time domain interference data (also known as the waveform slope mutation index) W can be:

[0062] In the above formula, is the signal slope of the signal at time t in the target signal interval, represents the signal amplitude at time t in the target signal interval, Indicates the starting point of the target signal interval, Indicates the end point of the target signal interval, that is, the signal impact main energy interval of the hammer impact force starts from Start to Finish.

[0063] Among them, the larger the value of the time domain interference data W, the more drastic the slope mutation of the measured impact signal during the impact process, the more serious the waveform distortion caused by force shunting and diffusion effects, and the more important it is to focus on correcting the high-frequency components or waveform sharpness in the candidate compensation.

[0064] In one embodiment, before obtaining the first energy center frequency and the second energy center frequency, the method further includes: Identifying a first frequency segment greater than a frequency threshold in a measured amplitude spectrum, wherein the measured amplitude spectrum is an amplitude spectrum formed after frequency domain conversion of the target signal interval; performing amplitude correction on each frequency component in the first frequency band to obtain a corrected amplitude of each frequency component in the first frequency band, wherein the corrected amplitude is the product of the amplitude of the corresponding frequency component in the measured amplitude spectrum and a correction parameter, wherein the correction parameter is determined based on the frequency value of the corresponding frequency component and the time-domain interference data, wherein the time-domain interference data and the correction parameter are positively correlated, and the correction parameter and the frequency value of the corresponding frequency component are positively correlated; Performing cubic spline interpolation according to the corrected amplitude of each frequency component in the first frequency band to obtain a target frequency band curve; The target frequency band curve and the second frequency band in the measured amplitude spectrum that is less than the frequency threshold are combined to obtain a target amplitude spectrum, which is used to represent the amplitude spectrum obtained by frequency domain conversion of the reference impulse signal.

[0065] Exemplarily, the frequency domain threshold is 10 kHz.

[0066] The first frequency segment in the measured amplitude spectrum that is identified to be greater than the frequency threshold may also be understood as the high-frequency portion in the measured amplitude spectrum.

[0067] In response to the situation where the reference impact signal is missing in actual applications, the amplitude spectrum of the reference impact signal obtained by frequency domain conversion is restored by extrapolating the high-frequency part of the measured impact signal (through cubic spline interpolation), so as to adapt to the situation where the frequency domain energy distribution of the reference impact signal follows a specific attenuation law (power-law attenuation) in the high frequency band (greater than 10kHz), while the frequency domain energy distribution of the measured impact signal still retains the residual characteristics of this law in the high frequency band. In this way, an accurate target amplitude spectrum can be easily obtained for subsequent frequency domain difference analysis, which further enhances the applicability of the method described in the present invention (without the need for additional labeling of the reference impact signal).

[0068] Among them, the reason for performing amplitude correction on each frequency component in the measured amplitude spectrum based on the time domain interference data is to compensate for the signal attenuation caused by force shunting and diffusion effects in the high-frequency part of the measured impact signal, so that the extrapolated recovered target amplitude spectrum is more accurate and reliable, thereby analyzing and obtaining more accurate frequency domain interference data.

[0069] In one example, if the amplitude of frequency f in the measured amplitude spectrum is set to , then the amplitude of frequency f in the target amplitude spectrum It can be expressed as:

[0070] in, represents the attenuation coefficient (the base of the log function is 10), a represents the baseline attenuation coefficient, which is set to 0.5 in this example (based on the experience of multiple tests, when a is set to 0.5, it can more accurately reflect the attenuation degree of the frequency amplitude). When the waveform slope mutation index (i.e. the aforementioned time domain interference data) W=1, it means that the time domain waveform of the measured signal has no significant slope mutation and is close to the ideal impact waveform. , which is used to represent the ideal high-frequency attenuation characteristics of stress waves in a theoretical undamped elastic rod. b represents the time domain distortion sensitivity coefficient, which is set to 0.1 in this example (based on the experience of multiple experiments, a b value of 0.1 can more accurately reflect the degree of influence of time domain changes on frequency amplitude attenuation). It reflects the modulation intensity of the time domain waveform distortion on the frequency domain attenuation index. When W increases by 10 times, Increase by 0.1.

[0071] In one embodiment, obtaining the first energy center frequency and the second energy center frequency includes: Performing frequency domain conversion on the time domain data to be processed to obtain a frequency domain amplitude spectrum, wherein the time domain data to be processed is the target signal interval or the reference impulse signal; Determining the product of the square of the amplitude of each frequency component in the frequency domain amplitude spectrum and each frequency component in the frequency domain amplitude spectrum as a first index value of each frequency component in the frequency domain amplitude spectrum; and determining the square of the amplitude of each frequency component in the frequency domain amplitude spectrum as a second index value of each frequency component in the frequency domain amplitude spectrum; Performing a sum calculation on the first index value of each frequency component in the frequency domain amplitude spectrum to obtain a first sum value, and performing a sum calculation on the second index value of each frequency component in the frequency domain amplitude spectrum to obtain a second sum value; When the time domain data to be processed is the target signal interval, the ratio of the first sum value to the second sum value is determined as the first energy center frequency; when the time domain data to be processed is the reference impulse signal, the ratio of the first sum value to the second sum value is determined as the second energy center frequency.

[0072] For example, the frequency domain conversion process may be completed by performing a fast Fourier transform (FFT) on the time domain data to be processed, thereby obtaining a corresponding frequency domain amplitude spectrum.

[0073] For example, the process of obtaining the first energy center frequency based on the target signal interval is: definition They represent the minimum frequency and maximum frequency in the frequency domain amplitude spectrum obtained by frequency domain conversion of the target signal interval, Represents the energy of frequency f in the frequency domain amplitude spectrum obtained by frequency domain conversion of the target signal interval (the square of the amplitude of frequency f, that is, the second index value), That is, it is the first index value of the frequency f in the frequency domain amplitude spectrum obtained by frequency domain conversion of the target signal interval.

[0074] The first energy center frequency It can be expressed as:

[0075] The first energy center frequency Sensitive to the high-frequency components of the signal, when the high-frequency energy is attenuated, It will shift to low frequency (less than 10kHz); on the contrary, if the high frequency energy is enhanced, Shift to high frequency (greater than 10kHz), that is, through the first energy center frequency It can quantify the strength of the high-frequency energy of the target signal interval in the frequency domain.

[0076] The second energy center frequency The acquisition process is similar to the above example and will not be described again to avoid repetition.

[0077] In this embodiment, the strength of the high-frequency energy of the corresponding signal interval in the frequency domain is quantified and used as the analysis indicator of the corresponding signal interval in the frequency domain to adapt to the situation where the high-frequency part of the impact signal is dominant, especially when the target amplitude spectrum corresponding to the reference impact signal is obtained by extrapolating and recovering the high-frequency part of the measured amplitude spectrum, so that the frequency domain interference data obtained by subsequent analysis can be more accurate and reliable.

[0078] It should be understood that in this case, the larger the value of the frequency domain interference data E, the more significant the high-frequency energy attenuation or low-frequency component diffusion in the measured impact signal. It reflects the drift of high-frequency energy to low-frequency, which is common in stress wave propagation of long piles or high-damping piles; negative offset may be caused by abnormal coupling of sensor installation.

[0079] It should be noted that if the second sum is 0, it can be understood that the signal corresponding to the time domain data to be processed is all noise. In this case, subsequent processing can be skipped and the user is reminded to re-collect a reliable impact signal for use.

[0080] In one embodiment, performing interference compensation on the measured impulse signal according to the time domain interference data and the frequency domain interference data to obtain a target impulse signal includes: Based on a pre-constructed inverse filter, inverse filtering is performed on a frequency domain signal corresponding to the measured impulse signal to obtain a compensated frequency domain signal, wherein the transfer function corresponding to the inverse filter includes a high-frequency gain compensation coefficient and a time domain distortion penalty coefficient, the high-frequency gain compensation coefficient is determined based on the frequency domain interference data and a frequency value of a corresponding frequency component, and the time domain distortion penalty coefficient is determined based on the time domain interference data; Performing a time domain transformation on the compensated frequency domain signal to obtain the target impact signal.

[0081] Exemplarily, the transfer function corresponding to the inverse filter can be expressed as:

[0082] in, represents the transfer function, represent the amplitude of frequency f in the target amplitude spectrum and the amplitude of frequency f in the measured amplitude spectrum, respectively. represents the high frequency gain compensation coefficient, , for those above The frequency band is increased by linear gain to compensate for high frequency attenuation. Indicates the sampling frequency (such as 40kHz), represents the time domain distortion penalty coefficient, , e is a natural constant, when W is greater than , high frequency compensation is forcibly enhanced.

[0083] In applications, the frequency domain signal corresponding to the measured impulse signal can be inversely filtered using a transfer function corresponding to the inverse filter to compensate for distortion of the measured impulse signal in the time domain and frequency domain, thereby obtaining a compensated frequency domain signal.

[0084] The target impact signal can be obtained by performing inverse Fourier transform on the compensated frequency domain signal to complete the aforementioned time domain change processing, that is, to obtain the hammer impact signal after improving the distortion caused by the force diversion effect, pile body damping and geometric dispersion effect.

[0085] The present invention proposes a hammer impact force measurement system suitable for strain pile dynamic testing. Figure 2 , which shows a schematic structural diagram of a hammer impact force measurement system 200 suitable for strain pile dynamic testing provided by one embodiment of the present invention, the system comprising: The interval identification module 201 is used to identify a target signal interval in a measured impact signal, wherein the measured impact signal is a signal generated when a hammer impacts a pile under test, and the target signal interval is used to indicate a main energy interval of the measured impact signal; A time domain analysis module 202 is configured to analyze the signal amplitude at each moment in the target signal interval to obtain time domain interference data, wherein the time domain interference data is used to indicate a severity of a change in the signal amplitude in the target signal interval; The frequency center acquisition module 203 is configured to acquire a first energy center frequency and a second energy center frequency, wherein the first energy center frequency is used to represent the center of the signal energy distribution corresponding to the measured impulse signal in the frequency domain, and the second energy center frequency is used to represent the center of the signal energy distribution corresponding to the reference impulse signal in the frequency domain, and the reference impulse signal is an ideal impulse signal corresponding to the measured impulse signal; A frequency domain analysis module 204 is configured to analyze a difference between the first energy center frequency and the second energy center frequency to obtain frequency domain interference data; The interference compensation module 205 is configured to perform interference compensation on the measured impact signal according to the time domain interference data and the frequency domain interference data to obtain a target impact signal, where the target impact signal is an impact signal used for dynamic detection of the strain pile.

[0086] In one embodiment, the interval identification module 201 includes: A window segmentation unit, configured to segment the measured impulse signal into a plurality of signal windows based on a preset sliding window length and sliding step size; a first determining unit, configured to determine a maximum window signal energy as a signal energy threshold among a plurality of window signal energies respectively corresponding to the plurality of signal windows, wherein the window signal energy is a sum of squares of the plurality of signal energies included in the corresponding signal window; a second determining unit, configured to determine a product of the signal energy threshold and a preset proportional coefficient as a standard threshold; An interval identification unit is used to identify the target signal interval in the measured impact signal based on the standard threshold and the multiple signal windows, wherein the average of the multiple window signal energies corresponding to the multiple signal windows covered by the target signal interval is greater than the standard threshold.

[0087] In one embodiment, the interval identification unit is specifically configured to: Identifying a first signal window and a second signal window in the multiple signal windows based on the standard threshold, wherein the first signal window is a first signal window among a plurality of consecutive signal windows whose corresponding window signal energy is greater than the standard threshold, and the second signal window is a signal window located in a preset order among a plurality of consecutive signal windows whose corresponding window signal energy is less than the standard threshold; The target signal interval is identified in the measured impulse signal by taking the window start point of the first signal window as the interval start point and the window end point of the second signal window as the interval end point.

[0088] In one embodiment, the time domain analysis module 202 is specifically configured to: Obtaining a signal slope at each moment in the target signal interval based on the signal amplitude at each moment in the target signal interval, wherein the signal slope is the instantaneous rate of change of the signal amplitude at the corresponding moment; Obtaining a slope mutation value at each moment in the target signal interval according to the signal slope at each moment in the target signal interval, wherein the slope mutation value is the absolute value of the second-order derivative of the signal slope at the corresponding moment; The sum of the slope mutation values at each moment in the target signal interval is calculated, and the sum is determined as the time domain interference data.

[0089] In one embodiment, the frequency center acquisition module 203 is specifically configured to: Performing frequency domain conversion on the time domain data to be processed to obtain a frequency domain amplitude spectrum, wherein the time domain data to be processed is the target signal interval or the reference impulse signal; Determining the product of the square of the amplitude of each frequency component in the frequency domain amplitude spectrum and each frequency component in the frequency domain amplitude spectrum as a first index value of each frequency component in the frequency domain amplitude spectrum; and determining the square of the amplitude of each frequency component in the frequency domain amplitude spectrum as a second index value of each frequency component in the frequency domain amplitude spectrum; Performing a sum calculation on the first index value of each frequency component in the frequency domain amplitude spectrum to obtain a first sum value, and performing a sum calculation on the second index value of each frequency component in the frequency domain amplitude spectrum to obtain a second sum value; When the time domain data to be processed is the target signal interval, the ratio of the first sum value to the second sum value is determined as the first energy center frequency; when the time domain data to be processed is the reference impulse signal, the ratio of the first sum value to the second sum value is determined as the second energy center frequency.

[0090] In one embodiment, the hammer impact force measurement system 200 for strain pile dynamic testing further includes an extrapolation module, which is specifically configured to: Identifying a first frequency segment greater than a frequency threshold in a measured amplitude spectrum, wherein the measured amplitude spectrum is an amplitude spectrum formed after frequency domain conversion of the target signal interval; performing amplitude correction on each frequency component in the first frequency band to obtain a corrected amplitude of each frequency component in the first frequency band, wherein the corrected amplitude is the product of the amplitude of the corresponding frequency component in the measured amplitude spectrum and a correction parameter, wherein the correction parameter is determined based on the frequency value of the corresponding frequency component and the time-domain interference data, wherein the time-domain interference data and the correction parameter are positively correlated, and the correction parameter and the frequency value of the corresponding frequency component are positively correlated; Performing cubic spline interpolation according to the corrected amplitude of each frequency component in the first frequency band to obtain a target frequency band curve; The target frequency band curve and the second frequency band in the measured amplitude spectrum that is less than the frequency threshold are combined to obtain a target amplitude spectrum, which is used to represent the amplitude spectrum obtained by frequency domain conversion of the reference impulse signal.

[0091] In one embodiment, the frequency domain analysis module 204 is specifically configured to: Calculating an absolute value of a difference between the first energy center frequency and the second energy center frequency; The frequency domain interference data is obtained according to the ratio of the absolute value and the second energy center frequency.

[0092] In one embodiment, the interference compensation module 205 is specifically configured to: Based on a pre-constructed inverse filter, inverse filtering is performed on a frequency domain signal corresponding to the measured impulse signal to obtain a compensated frequency domain signal, wherein the transfer function corresponding to the inverse filter includes a high-frequency gain compensation coefficient and a time domain distortion penalty coefficient, the high-frequency gain compensation coefficient is determined based on the frequency domain interference data and a frequency value of a corresponding frequency component, and the time domain distortion penalty coefficient is determined based on the time domain interference data; Performing a time domain transformation on the compensated frequency domain signal to obtain the target impact signal.

[0093] It should be noted that the system provided in the above embodiment is merely illustrated by the division of the aforementioned functional modules. In actual applications, the aforementioned functions can be assigned to different functional modules as needed, i.e., the internal structure of the computer device can be divided into different functional modules to perform all or part of the functions described above. Furthermore, the hammer impact force measurement system for strain pile dynamic testing and the hammer impact force measurement method for strain pile dynamic testing provided in the above embodiment are based on the same concept. The specific implementation process is detailed in the method embodiment and will not be repeated here.

[0094] The embodiment of the present invention also provides an electronic device. Figure 3 , the electronic device may include a processor 301, a memory 302, and a program 3021 stored in the memory 302 and executable on the processor 301.

[0095] When the program 3021 is executed by the processor 301, it can achieve Figure 1 Any steps in the corresponding method embodiments and achieving the same beneficial effects will not be repeated here.

[0096] Those skilled in the art will appreciate that all or part of the steps of implementing the above-described embodiment method may be accomplished through hardware associated with program instructions, and the program may be stored in a readable medium.

[0097] The embodiment of the present invention further provides a readable storage medium, wherein the readable storage medium stores a computer program, and when the computer program is executed by a processor, the above Figure 1 Any steps in the corresponding method embodiments can achieve the same technical effects and will not be described again here to avoid repetition.

[0098] The computer-readable storage medium of the embodiments of the present invention may adopt any combination of one or more computer-readable media. The computer-readable medium may be a computer-readable signal medium or a computer-readable storage medium. The computer-readable storage medium may be, for example, an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, device, or component, or any combination thereof. More specific examples (a non-exhaustive list) of computer-readable storage media include: an electrical connection with one or more wires, a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination thereof. In this document, a computer-readable storage medium may be any tangible medium containing or storing a program that can be used by or in conjunction with an instruction execution system, device, or component.

[0099] A computer-readable signal medium may include a data signal propagated in baseband or as part of a carrier wave, which carries computer-readable program code. Such propagated data signals may take various forms, including but not limited to electromagnetic signals, optical signals, or any suitable combination thereof. A computer-readable signal medium may also be any computer-readable medium other than a computer-readable storage medium that can transmit, propagate, or transport a program for use by or in conjunction with an instruction execution system, apparatus, or device.

[0100] The program code contained on the storage medium may be transmitted using any appropriate medium, including but not limited to wireless, wireline, optical fiber cable, RF, etc., or any suitable combination of the foregoing.

[0101] Computer program code for performing the operations of the present invention may be written in one or more programming languages, or a combination thereof, including object-oriented programming languages such as Java, Smalltalk, C++, and conventional procedural programming languages such as "C" or similar programming languages. The program code may be executed entirely on the user's computer, partially on the user's computer, as a stand-alone software package, partially on the user's computer and partially on a remote computer, or entirely on a remote computer or terminal. In cases involving a remote computer, the remote computer may be connected to the user's computer via any type of network, including a local area network (LAN) or a wide area network (WAN), or may be connected to an external computer (e.g., via the Internet using an Internet service provider).

[0102] An embodiment of the present invention further provides a computer program product. When the computer program product is run on a computer, the computer executes the above-mentioned related steps to implement a hammer impact force measurement method suitable for strain pile dynamic testing provided by the above embodiment.

[0103] It should be noted that the order in which the embodiments of the present invention are described above is for illustrative purposes only and does not necessarily represent the superiority or inferiority of the embodiments. The processes depicted in the accompanying drawings do not necessarily require the specific order or sequential order shown to achieve the desired results. In certain embodiments, multitasking and parallel processing are also possible or may be advantageous.

[0104] The various embodiments in this specification are described in a progressive manner, and the same or similar parts between the various embodiments can be referred to each other. Each embodiment focuses on the differences from other embodiments.

Claims

1. A hammer impact force measurement method suitable for strain pile dynamic testing, characterized in that: The method comprises: Identifying a target signal interval in a measured impact signal, wherein the measured impact signal is a signal generated when a hammer impacts a pile under test, and the target signal interval is used to indicate a main energy interval of the measured impact signal; Analyzing the signal amplitude at each moment in the target signal interval to obtain time domain interference data, where the time domain interference data is used to indicate a severity of a change in the signal amplitude in the target signal interval; Obtaining a first energy center frequency and a second energy center frequency, wherein the first energy center frequency is used to represent the center of the signal energy distribution corresponding to the measured impulse signal in the frequency domain, and the second energy center frequency is used to represent the center of the signal energy distribution corresponding to the reference impulse signal in the frequency domain, and the reference impulse signal is an ideal impulse signal corresponding to the measured impulse signal; Analyzing a difference between the first energy center frequency and the second energy center frequency to obtain frequency domain interference data; According to the time domain interference data and the frequency domain interference data, interference compensation is performed on the measured impact signal to obtain a target impact signal, which is an impact signal used for dynamic detection of the inspected pile.

2. The hammer impact force measurement method for strain pile dynamic testing according to claim 1 is characterized in that: The identifying of the target signal interval in the measured impact signal includes: Based on a preset sliding window length and sliding step size, the measured impact signal is divided into a plurality of signal windows; Determine the maximum window signal energy among the multiple window signal energies respectively corresponding to the multiple signal windows as the signal energy threshold, wherein the window signal energy is the square sum of the multiple signal energies included in the corresponding signal window; Determine the product of the signal energy threshold and a preset proportional coefficient as a standard threshold; The target signal interval is identified in the measured impulse signal according to the standard threshold and the multiple signal windows, wherein the average of multiple window signal energies corresponding to the multiple signal windows covered by the target signal interval is greater than the standard threshold.

3. The hammer impact force measurement method suitable for strain pile dynamic testing according to claim 2, characterized in that: The step of identifying the target signal interval in the measured impulse signal according to the standard threshold and the plurality of signal windows includes: Identifying a first signal window and a second signal window in the multiple signal windows based on the standard threshold, wherein the first signal window is a first signal window among a plurality of consecutive signal windows whose corresponding window signal energy is greater than the standard threshold, and the second signal window is a signal window located in a preset order among a plurality of consecutive signal windows whose corresponding window signal energy is less than the standard threshold; The target signal interval is identified in the measured impulse signal by taking the window start point of the first signal window as the interval start point and the window end point of the second signal window as the interval end point.

4. The hammer impact force measurement method suitable for strain pile dynamic testing according to claim 1, characterized in that: The analyzing the signal amplitude at each moment in the target signal interval to obtain time domain interference data includes: Obtaining a signal slope at each moment in the target signal interval based on the signal amplitude at each moment in the target signal interval, wherein the signal slope is the instantaneous rate of change of the signal amplitude at the corresponding moment; Obtaining a slope mutation value at each moment in the target signal interval according to the signal slope at each moment in the target signal interval, wherein the slope mutation value is the absolute value of the second-order derivative of the signal slope at the corresponding moment; The sum of the slope mutation values at each moment in the target signal interval is calculated, and the sum is determined as the time domain interference data.

5. The hammer impact force measurement method suitable for strain pile dynamic testing according to claim 1, characterized in that: The obtaining of the first energy center frequency and the second energy center frequency includes: Performing frequency domain conversion on the time domain data to be processed to obtain a frequency domain amplitude spectrum, wherein the time domain data to be processed is the target signal interval or the reference impulse signal; Determining the product of the square of the amplitude of each frequency component in the frequency domain amplitude spectrum and each frequency component in the frequency domain amplitude spectrum as a first index value of each frequency component in the frequency domain amplitude spectrum; and determining the square of the amplitude of each frequency component in the frequency domain amplitude spectrum as a second index value of each frequency component in the frequency domain amplitude spectrum; Performing a sum calculation on the first index value of each frequency component in the frequency domain amplitude spectrum to obtain a first sum value, and performing a sum calculation on the second index value of each frequency component in the frequency domain amplitude spectrum to obtain a second sum value; When the time domain data to be processed is the target signal interval, the ratio of the first sum value to the second sum value is determined as the first energy center frequency; when the time domain data to be processed is the reference impulse signal, the ratio of the first sum value to the second sum value is determined as the second energy center frequency.

6. The hammer impact force measurement method suitable for strain pile dynamic testing according to claim 1, characterized in that: Before obtaining the first energy center frequency and the second energy center frequency, the method further includes: Identifying a first frequency segment greater than a frequency threshold in a measured amplitude spectrum, wherein the measured amplitude spectrum is an amplitude spectrum formed after frequency domain conversion of the target signal interval; performing amplitude correction on each frequency component in the first frequency band to obtain a corrected amplitude of each frequency component in the first frequency band, wherein the corrected amplitude is the product of the amplitude of the corresponding frequency component in the measured amplitude spectrum and a correction parameter, wherein the correction parameter is determined based on the frequency value of the corresponding frequency component and the time-domain interference data, wherein the time-domain interference data and the correction parameter are positively correlated, and the correction parameter and the frequency value of the corresponding frequency component are positively correlated; Performing cubic spline interpolation according to the corrected amplitude of each frequency component in the first frequency band to obtain a target frequency band curve; The target frequency band curve and the second frequency band in the measured amplitude spectrum that is less than the frequency threshold are combined to obtain a target amplitude spectrum, which is used to represent the amplitude spectrum obtained by frequency domain conversion of the reference impulse signal.

7. The hammer impact force measurement method for strain pile dynamic testing according to claim 1, characterized in that: The analyzing the difference between the first energy center frequency and the second energy center frequency to obtain frequency domain interference data includes: Calculating an absolute value of a difference between the first energy center frequency and the second energy center frequency; The frequency domain interference data is obtained according to the ratio of the absolute value and the second energy center frequency.

8. The hammer impact force measurement method suitable for strain pile dynamic testing according to claim 1, characterized in that: The performing interference compensation on the measured impulse signal according to the time domain interference data and the frequency domain interference data to obtain a target impulse signal includes: Based on a pre-constructed inverse filter, inverse filtering is performed on a frequency domain signal corresponding to the measured impulse signal to obtain a compensated frequency domain signal, wherein the transfer function corresponding to the inverse filter includes a high-frequency gain compensation coefficient and a time domain distortion penalty coefficient, the high-frequency gain compensation coefficient is determined based on the frequency domain interference data and a frequency value of a corresponding frequency component, and the time domain distortion penalty coefficient is determined based on the time domain interference data; Performing a time domain transformation on the compensated frequency domain signal to obtain the target impact signal.

9. A hammer impact force measurement system suitable for strain pile dynamic testing, characterized in that: The system comprises: an interval identification module, configured to identify a target signal interval in a measured impact signal, wherein the measured impact signal is a signal generated when a hammer impacts a pile under test, and the target signal interval is configured to indicate a main energy interval of the measured impact signal; a time domain analysis module, configured to analyze the signal amplitude at each moment in the target signal interval to obtain time domain interference data, wherein the time domain interference data is used to indicate a severity of a change in the signal amplitude in the target signal interval; a frequency center acquisition module, configured to acquire a first energy center frequency and a second energy center frequency, wherein the first energy center frequency is used to represent the center of the signal energy distribution corresponding to the measured impulse signal in the frequency domain, and the second energy center frequency is used to represent the center of the signal energy distribution corresponding to the reference impulse signal in the frequency domain, and the reference impulse signal is an ideal impulse signal corresponding to the measured impulse signal; a frequency domain analysis module, configured to analyze a difference between the first energy center frequency and the second energy center frequency to obtain frequency domain interference data; The interference compensation module is used to perform interference compensation on the measured impact signal according to the time domain interference data and the frequency domain interference data to obtain a target impact signal, where the target impact signal is an impact signal used for dynamic detection of the strain pile.

10. An electronic device, characterized in that: The invention comprises a processor, a memory and a computer program stored in the memory and executable on the processor, wherein when the computer program is executed by the processor, the steps of the hammer impact force measurement method applicable to dynamic detection of strain piles as claimed in any one of claims 1 to 8 are implemented.

Images