Method for acquiring propagation law of blasting vibration surface wave along human body surface

Abstract

The present application belongs to the field of blasting engineering and blasting safety, and relates to a method for obtaining the propagation law of blasting vibration surface wave along the human body surface. The method comprises the following steps: S1: carrying out in-situ blasting vibration test to obtain the measured ground vibration velocity, and inversely deducing the input boundary conditions for numerical simulation according to the peak ground vibration velocity; S2: constructing a human body finite element model containing bones and soft tissues, and verifying the rationality of the model; S3: carrying out the upward transmission dynamics simulation based on the verified model, and extracting the peak vibration velocity sequence of multiple continuous physiological nodes along the preset specific propagation path; S4: combining the nonlinear regression method to fit and obtain the vibration transmission reduction coefficient between adjacent physiological nodes, and then constructing the reduction prediction formula representing the energy dissipation law. The present application realizes the safe and reliable acquisition of the propagation and attenuation characteristics of blasting vibration on the human body surface, and provides precise theoretical support for individual safety evaluation and protection design.

Description

Technical Field

[0001] This invention relates to the fields of blasting engineering and blasting safety, and in particular to a method for obtaining the propagation law of blasting vibration surface waves along the human body surface. Background Technology

[0002] Engineering blasting is widely used in mining, transportation infrastructure, and urban renewal, but the resulting blasting vibrations can not only affect surrounding buildings but also cause discomfort and even injury to people. As engineering blasting gradually extends to complex environments and densely populated areas, how to scientifically and accurately assess the dynamic impact of blasting vibrations on the human body has become a critical issue that urgently needs to be addressed in the field of blasting safety protection.

[0003] Currently, the methods for obtaining the dynamic response law of the human body under blasting vibration are mainly divided into two categories: on-site real-person in-situ tests and computer numerical simulation: (1) On-site real-person in-situ tests: Traditional tests mostly rely on attaching contact sensors to the surface of the subject's body. However, this method has obvious defects: on the one hand, conducting real-person tests under high-intensity blasting vibration conditions poses a great safety risk, making it difficult to obtain measured data in the high vibration velocity range; on the other hand, the additional mass and local constraint effect of the contact sensors and their cables will interfere with the real vibration response of the human soft tissue, resulting in measurement distortion. (2) Traditional computer numerical simulation: Although the existing numerical simulation methods avoid personnel safety risks, they often face the problem of inaccurate input boundary conditions in actual operation, making it difficult to convert the real blasting seismic waves on-site into accurate dynamic input. In addition, the anatomical structure of human bones and soft tissues is extremely complex, with highly nonlinear and viscoelastic characteristics. Traditional simulation models often lack accurate verification and repeated parameter calibration for real body surface dynamic response data, resulting in low accuracy of the vibration wave propagation and attenuation law obtained from the simulation, making it difficult to provide a reliable reference in actual human safety assessment.

[0004] Therefore, there is an urgent need for a systematic method that can avoid the safety risks of high-intensity blasting tests on real people, overcome the interference of traditional contact measurement, and have high simulation accuracy, so as to safely, reliably, and quantitatively obtain the propagation and attenuation laws of blasting vibration surface waves between various physiological nodes on the complex human body surface. Summary of the Invention

[0005] In order to solve the above-mentioned technical problems in the background art, the present invention provides a method for obtaining the propagation law of blasting vibration surface waves along the human body surface, which can avoid the safety risks of high-intensity human body experiments and obtain reliable results.

[0006] To achieve the above objectives, the present invention adopts the following technical solution: A method for obtaining the propagation law of blasting vibration surface waves along the human body surface, characterized by comprising the following steps: 1) Conduct in-situ blasting vibration tests to obtain measured surface vibration velocity, fit the measured data using empirical formulas for surface vibration attenuation, obtain the peak surface vibration velocity under the target blasting condition, and use this to inversely deduce the input boundary conditions for numerical simulation of human dynamic response. 2) Construct a finite element model of the human body including skeletal and soft tissue structures. Apply the input boundary conditions obtained in step 1) to the foot end of the finite element model of the human body for dynamic simulation. Extract the vibration velocity of the human body surface nodes obtained from the simulation and compare it with the human body surface vibration data measured by digital image correlation technology to verify the rationality of the finite element model of the human body. 3) Based on the verified human finite element model obtained in step 2), conduct upward dynamic simulation and extract the peak vibration velocity sequence of multiple continuous physiological nodes along the preset specific propagation path on the human body surface from bottom to top. 4) Based on the peak ground vibration velocity obtained in step 1), the initial reference value of the input is used. Based on the peak ground vibration velocity sequence, the vibration transmission reduction coefficient between each adjacent physiological node is obtained by nonlinear regression method. Then, a reduction prediction formula is constructed to characterize the energy dissipation law of the blasting vibration surface wave transmitted upward along the human body surface.

[0007] Preferably, the empirical formula for surface vibration attenuation and the input boundary conditions for numerical simulation of human dynamic response obtained by reverse calculation in step 1) of the present invention are specifically as follows: 1.1) The empirical formula for surface vibration attenuation is the Sadovsky formula, which is based on the single-hole blasting charge obtained from the field conditions. Q , the center of explosion distance R and peak vibration velocity of blast v An empirical formula for surface vibration attenuation was fitted to eliminate random errors from a single field test; the Sadovsky empirical formula is as follows: ; In the formula: K , α A coefficient related to topographic and geological conditions; 1.2) The method for obtaining the input boundary conditions by reverse calculation is as follows: based on the formula The peak ground vibration velocity under the target blasting condition is converted into equivalent stress and used as the input boundary condition, wherein... σ For equivalent stress, ρ For the density of the site medium, c For the longitudinal wave velocity of the site medium, v This represents the peak ground vibration velocity.

[0008] Preferably, in step 2) of the present invention, the process of acquiring human body surface vibration data measured using digital image correlation technology includes: 2.1) Collect high-speed image sequences of the body surface of subjects wearing speckled bodysuits during the blasting process; 2.2) Calculate and extract the measured vibration velocities in the vertical and horizontal directions of the measuring point using image correlation; 2.3) The speckle diameter and spacing on the surface of the speckle tights are preset dimensions to construct a tracking reference for non-contact optical measurement.

[0009] Preferably, in step 2) of the present invention, the human finite element model is scaled down according to specific height and weight parameters, and its internal tissues are assigned corresponding equivalent stiffness and equivalent density characteristic parameters according to biomechanical properties, and the external excitation loading method of the dynamic simulation is stress boundary loading.

[0010] Preferably, in step 3) of the present invention, the consecutive physiological nodes of the preset specific propagation path are as follows: 3.1) The core transmission path includes: foot input → foot → calf → knee joint → thigh → pelvis and buttocks → abdomen → chest → shoulder → neck → face → top of the head; It can also include, but is not limited to, the following paths: 3.2) Foot input → sole of foot → ankle → calf → knee joint; 3.3) Foot input → sole of foot → calf → knee joint → thigh → pelvis and buttocks; 3.4) Foot input → sole of foot → calf → knee joint → thigh → pelvis and hips → abdomen → chest; 3.5) Foot input → foot → calf → knee joint → thigh → pelvis and hips → abdomen → chest → shoulder; 3.6) Foot input → foot → calf → knee joint → thigh → pelvis and hip → abdomen → chest → shoulder → upper arm → elbow → forearm → wrist and hand.

[0011] Preferably, in step 4) of the present invention, the stepwise reduction prediction formula is used to quantitatively characterize the vibration energy loss pattern between adjacent physiological nodes due to differences in anatomical structure, and its mathematical expression is as follows: PPV n = η n • PPV n-1 ; In the formula, PPV n The predicted peak vibration velocity of the nth physiological node on the specific propagation path;PPV n-1 Let n be the predicted peak vibration velocity of the (n-1)th physiological node on this path, and when n=1, v n-1 It is equal to the initial reference value of the input terminal; η n It is the stepwise vibration transmission reduction factor between the (n-1)th node and the nth node, reflecting the buffering and viscoelastic damping characteristics of specific physiological structures in the human body.

[0012] The step-by-step vibration transmission reduction factor η n The simulation node dataset, which includes various combinations of single-segment maximum charge and detonation center distance, was extracted and obtained by multivariate nonlinear regression fitting using the least squares method.

[0013] The beneficial effects of this invention are as follows: (1) This invention introduces non-contact in-situ measurement by combining speckle tight-fitting clothing with DIC digital image correlation technology, which avoids the interference of the added mass and local constraints of traditional contact sensors on the actual vibration of the body surface, and can effectively avoid the safety risks of directly carrying out real-person experiments under high-intensity blasting conditions. (2) This invention uses the measured vibration velocity of the ground surface combined with the Sadovsky formula to back-calculate the equivalent stress as the input boundary condition of the finite element model, and uses the DIC measured data to dynamically calibrate the biomechanical parameters of the human body model, which effectively solves the problems of external load input distortion and low model accuracy in traditional numerical simulation. (3) Based on the peak velocity extraction and nonlinear regression fitting of multiple sets of continuous physiological nodes, this invention constructs a reduction prediction formula to characterize the upward transmission of blast surface waves along specific anatomical structures of the human body, quantitatively characterizes the dissipation characteristics of vibration energy, and provides reliable data support for human blast safety assessment and individual protection design. Attached Figure Description

[0014] Figure 1 This is a flowchart illustrating a method for obtaining the propagation pattern of blasting vibration surface waves along the human body surface, as provided by this invention.

[0015] Figure 2 This is a schematic diagram of the in-situ test blasting location used in this invention.

[0016] Figure 3 It is the fitting result of the maximum single-stage charge, the distance from the detonation center, and the peak vibration velocity at the measuring point obtained from the in-situ test used in this invention.

[0017] Figure 4 This is a schematic diagram of the human body numerical modeling used in this invention.

[0018] Figure 5This is a schematic diagram illustrating the verification of the rationality of the numerical modeling obtained by comparing the numerical modeling results used in this invention with the measured data.

[0019] Figure 6 This is a numerical schematic diagram illustrating the attenuation pattern predicted along a specific propagation path in this invention. Detailed Implementation

[0020] The present invention will now be described in further detail with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and not intended to limit it. Example 1

[0021] like Figure 1 As shown, this embodiment provides a method for obtaining the propagation law of blasting vibration surface waves along the human body surface, including: S1. After the in-situ blasting vibration test, the measured vibration velocity on the ground surface is obtained. The measured data is fitted using the empirical formula for ground vibration attenuation to obtain the peak ground vibration velocity under the target blasting condition. Based on this, the input boundary conditions for numerical simulation of human dynamic response are deduced.

[0022] In some specific embodiments, step S1 specifically includes: S11. Select a blasting test site that meets the requirements of the "Blasting Safety Regulations" and conduct in-situ human monitoring tests for blasting vibrations. S12. Complete the site stratum physical and mechanical parameter tests to determine the core parameters such as soil density, longitudinal wave velocity, and Poisson's ratio. The test site is about 500m long and 400m wide. According to the site survey report, the surface layer within the test area is 50cm thick plain fill, and the lower layer is silty clay. The soil layer has a single property and has not been disturbed. Before the test, borehole sampling was carried out on the test site, with a borehole depth of 4m. In this embodiment, the physical and mechanical parameters obtained at the test site are shown in Table 1.

[0023] Table 1 Physical and mechanical parameters of the test site

[0024] S13. The average height of adult males and adult females aged 18-44 in my country is 169.7 cm and the average weight is 69.6 kg. In this embodiment, four different working conditions were set up, and a male volunteer with a height and weight similar to the average height and weight in my country was recruited for the on-site blasting test. In this embodiment, the operating conditions for the in-situ test are shown in Table 2.

[0025] Table 2. In-situ blasting test condition settings

[0026] S14. A blasting vibration monitoring instrument was installed at the standing location of the male volunteers to simultaneously collect time-history data of surface vibration velocity under various working conditions; the measured data were extracted, and the Sadovsky formula was used for fitting analysis to obtain the theoretical value of the peak surface vibration velocity under the target blasting condition, in order to eliminate the random error of a single field test; the in-situ test blasting location is as follows: Figure 2 As shown; In this embodiment, the field monitoring data of the in-situ blasting vibration test obtained by statistics are shown in Table 3.

[0027] Table 3 Summary of On-site Vibration Monitoring

[0028] The fitting results are as follows Figure 3 As shown, α The value is 0.831. k The value is 68.91; the fitted Sadovsky formula is: ; S15, According to the formula The peak ground vibration velocity obtained under the target blasting condition will be... v Converted into equivalent stress σ This serves as the input boundary condition for subsequent numerical simulations of human dynamics. In this embodiment, the peak ground vibration velocity is extracted under the target condition (charge weight 0.9 kg, detonation center distance 7.5 m). v The density of the site medium is 0.1254 m / s. ρ Based on the measured natural unit weight parameters of the plain fill soil in Table 1, the converted value is 1800 kg / m³. 3 Longitudinal wave velocity of the site medium c The value was taken as 220 m / s based on the geological survey and testing report of the test site; the peak equivalent stress was calculated based on the above parameters. σ The equivalent stress was 49.66 kPa, and then the equivalent stress was... σ The load is converted into a time history curve and then uniformly applied to the sole of the foot of the human finite element model using the DEFINE_CURVE keyword.

[0029] S2. Construct a finite element model of the human body including skeletal and soft tissue structures, and verify its rationality: In some specific embodiments, step S2 specifically includes: S21. Based on the Toyota THUMS 50th percentile standard adult male (AM50) finite element model, the overall height of the human model was scaled to 170 cm and the weight was adjusted to 60 kg according to the real physical parameters of the subjects in the on-site in-situ test, using the scale tool (Transf_Scale) in the LS-PREPOST post-processing software to highly restore the real geometric spatial characteristics of the subjects' bodies. S22. Based on the geometric and anatomical structures and biomechanical properties of different human tissues, hybrid element technology is used to mesh and assign constitutive materials to the scaled-down human finite element model. Specific modeling details are as follows: Figure 4 As shown: (1) For large-volume bones, muscles and internal organs, hexahedral and tetrahedral solid elements are used for modeling; among them, for the tetrahedral elements used for internal organs and joints with complex geometry, a single-point integration algorithm with viscous hourglass control (ELFORM=13) is selected to prevent volume self-locking. (2) For cortical bone, fascia and human skin with a thickness of less than 1 mm, shell elements are used for modeling; the nodes of the outermost skin shell element are used to directly map and extract the dynamic response data of the body surface space measurement points corresponding to the actual digital image correlation (DIC) technique. (3) For cortical bone and cancellous bone tissue, an isotropic elastoplastic constitutive model (MAT_ISOTROPIC_ELASTIC_PLASTIC) is used to characterize the nonlinear deformation and hardening behavior of bones; for muscles and most soft tissues, a simplified rubber / foam hyperelastic constitutive model (MAT_SIMPLIFIED_RUBBER / FOAM) is used to simulate their significant nonlinear stiffening and buffering characteristics under passive tension and compression. S23. Numerical simulation is carried out using the LS-DYNA explicit dynamics solver, combined with the target working condition peak equivalent stress obtained in step S15. σ By defining the load time history curve (DEFINE_CURVE keyword), the equivalent stress is... σ The external excitation is applied directly and uniformly to the foot contact mesh of the human finite element model, and explicit dynamic solution is initiated to extract the vibration velocity time history data of the preset measurement point nodes on the body surface. S24. Simultaneously acquire human body surface vibration data measured using digital image correlation technology in the in-situ experiment; volunteers wear tight-fitting clothing with dense speckle patterns (approximately 0.5 cm in diameter and 0.25 cm in spacing) to construct a tracking benchmark for non-contact optical measurement. The measured vibration velocity at the measurement points is extracted by calculating the correlation between a high-speed camera and the image. Compared with the field experiment, the attenuation law of the peak vibration velocity of the human body surface measurement point unit in the numerical simulation is basically consistent, with a small difference in peak values. Figure 5 As shown, it can effectively simulate field tests and has high reliability; In this embodiment, the measured data and the vibration velocities of the preset measuring points on the body surface obtained by explicit dynamic solution are shown in Table 4.

[0030] Table 4 Comparison of measured data and modeled data

[0031] S3. Extract node dynamic response sequence data along a specific propagation path: In some specific embodiments, step S3 specifically includes: Based on the refined human finite element model verified in step S2, an explicit dynamic simulation of the upward propagation of blast surface waves is carried out. Along a predetermined specific propagation path on the human body surface (foot input → foot → calf → knee joint → thigh → pelvis and buttocks → abdomen → chest → shoulder → neck → face → skull top), 11 physiological nodes corresponding to the shell unit of the skin layer are selected as target monitoring points. The full-time dynamic response data of the nodes in this sequence are extracted, and the peak vibration velocity sequence of each measuring point node is statistically obtained. It is known that the specific transmission path pre-set along the human body surface can also be: foot input → foot → ankle → calf → knee joint; foot input → foot → calf → knee joint → thigh → pelvis and buttocks; foot input → foot → calf → knee joint → thigh → pelvis and buttocks → abdomen → chest; foot input → foot → calf → knee joint → thigh → pelvis and buttocks → abdomen → chest → shoulder; foot input → foot → calf → knee joint → thigh → pelvis and buttocks → abdomen → chest → shoulder → upper arm → elbow → forearm → wrist and hand, at least one of these paths.

[0032] S4. Based on the peak vibration velocity sequence, a nonlinear regression method is used to fit and obtain the vibration transmission reduction coefficient between adjacent physiological nodes, and a reduction prediction formula is constructed to characterize the energy dissipation law of the blast vibration surface wave transmitted upward along the human body surface: In some specific embodiments, step S4 specifically includes: S41. Extract the peak ground vibration velocity sequence obtained from each in-situ test condition in step S15 (including multiple test conditions with different detonation distances from 7.5 m to 17.5 m under the same single charge), and use it as the initial reference value for the input end of the specific propagation path under the corresponding test condition. PPV 0; S42. Extract the peak vibration velocity sequences of multiple sets of specific propagation path measuring points obtained in step S3 under the corresponding working conditions to form a dynamic response node dataset containing multiple combinations of single-segment maximum charge and detonation center distance characteristics. S43. Using adjacent physiological node pairs along a specific propagation path as calculation units, the dynamic response node dataset is called, and the least squares method is used to perform multivariate nonlinear regression fitting calculation to obtain the optimal stepwise vibration transmission reduction coefficient that takes into account the characteristics of data under various working conditions and can accurately reflect the attenuation characteristics of specific physiological structures of the human body. η n ; S44. Construct a reduction prediction formula to characterize the energy dissipation law of blast vibration surface waves transmitted upward along the human body surface: PPV n = η n • PPV n-1 In the formula, PPV n Let n be the predicted peak vibration velocity of the nth physiological node. PPV n-1 The predicted peak vibration velocity of the (n-1)th physiological node is shown in Table 5. To clearly demonstrate the calculation logic and verification process of the prediction formula, the input benchmark value of 125.4 mm / s for the target blasting condition 1 and the corresponding extracted mean response values ​​of each measuring point are listed together with the fitting reduction coefficients to verify the accuracy of the step-by-step reduction law. In this embodiment, the peak vibration velocity of the target detection point on a specific propagation path in the numerical simulation is shown in Table 5.

[0033] Table 5 Vibration velocities along specific propagation paths

[0034] Based on the successive vibration transmission reduction coefficient matrix obtained in Table 5, the predicted attenuation values ​​along a specific propagation path in this embodiment are as follows: Figure 6 As shown, the final complete reduction prediction formula system representing the upward propagation of blast surface waves along the core path of the human body surface is as follows: PPV 1 = 0.030 PPV 0; PPV 2 = 0.564 PPV1; PPV 3 = 0.784 PPV 2; PPV 4 = 0.516 PPV 3; PPV 5 = 0.145 PPV 4; PPV 6 = 0.151 PPV 5; PPV 7 = 0.868 PPV 6; PPV 8 = 0.909 PPV 7; PPV 9 = 0.800 PPV 8; PPV 10 = 0.750 PPV 9; PPV 11 = 0.889 PPV 10 .

[0035] In summary, this invention constructs an accurate reduction prediction formula, quantitatively characterizes the dissipation characteristics of vibration energy, and provides reliable data support for human explosion safety assessment and individual protection design.

[0036] The above embodiments are only used to illustrate the technical solutions of the present invention, and are not intended to limit it. Although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some or all of the technical features therein. Such modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the scope of the technical solutions of the embodiments of the present invention, and should all be included within the protection scope of the present invention.

Claims

1. A method for obtaining the propagation law of blasting vibration surface waves along the human body surface, characterized in that, Includes the following steps: 1) Conduct in-situ blasting vibration tests to obtain measured surface vibration velocity, fit the measured data using empirical formulas for surface vibration attenuation, obtain the peak surface vibration velocity under the target blasting condition, and use this to inversely deduce the input boundary conditions for numerical simulation of human dynamic response. 2) Construct a finite element model of the human body including skeletal and soft tissue structures. Apply the input boundary conditions obtained in step 1) to the foot end of the finite element model of the human body for dynamic simulation. Extract the vibration velocity of the human body surface nodes obtained from the simulation and compare it with the human body surface vibration data measured by digital image correlation technology to verify the rationality of the finite element model of the human body. 3) Based on the verified human finite element model obtained in step 2), conduct upward dynamic simulation and extract the peak vibration velocity sequence of multiple continuous physiological nodes along the preset specific propagation path on the human body surface from bottom to top. 4) Based on the peak ground vibration velocity obtained in step 1), the initial reference value of the input is used. Based on the peak ground vibration velocity sequence, the vibration transmission reduction coefficient between each adjacent physiological node is obtained by nonlinear regression method. Then, a reduction prediction formula is constructed to characterize the energy dissipation law of the blasting vibration surface wave transmitted upward along the human body surface.

2. The method for obtaining the propagation law of blasting vibration surface waves along the human body surface according to claim 1, characterized in that: The empirical formula for surface vibration attenuation and the input boundary conditions for numerical simulation of human dynamic response obtained by reverse calculation in step 1) are as follows: 1.1) The empirical formula for surface vibration attenuation is the Sadovsky formula, which is based on the single-hole blasting charge obtained from the field conditions. Q , the center of explosion distance R and peak vibration velocity of blast v An empirical formula for surface vibration attenuation was fitted to eliminate random errors from a single field test; the Sadovsky empirical formula is as follows: ； In the formula: K , α A coefficient related to topographic and geological conditions; 1.2) The method for obtaining the input boundary conditions by reverse calculation is as follows: based on the formula The peak ground vibration velocity under the target blasting condition is converted into equivalent stress and used as the input boundary condition, wherein... σ For equivalent stress, ρ For the density of the site medium, c For the longitudinal wave velocity of the site medium, v This represents the peak ground vibration velocity.

3. The method for obtaining the propagation law of blasting vibration surface waves along the human body surface according to claim 1, characterized in that, In step 2), the human finite element model is scaled down according to specific height and weight parameters, and its internal tissues are assigned corresponding equivalent stiffness and equivalent density characteristic parameters according to biomechanical properties. The external excitation loading method of the dynamic simulation is stress boundary loading.

4. The method for obtaining the propagation law of blasting vibration surface waves along the human body surface according to claim 3, characterized in that: In step 2), the process of acquiring human body surface vibration data measured using digital image correlation technology includes: 2.1) Collect high-speed image sequences of the body surface of subjects wearing speckled bodysuits during the blasting process; 2.2) Calculate and extract the measured vibration velocities in the vertical and horizontal directions of the measuring point using image correlation; 2.3) The speckle diameter and spacing on the surface of the speckle tights are preset dimensions to construct a tracking reference for non-contact optical measurement.

5. The method for obtaining the propagation law of blasting vibration surface waves along the human body surface according to claim 1, characterized in that, In step 3), the sequential physiological nodes of the preset specific transmission path are: foot input → foot → calf → knee joint → thigh → pelvis and buttocks → abdomen → chest → shoulder → neck → face → top of the head.

6. The method for obtaining the propagation law of blasting vibration surface waves along the human body surface according to claim 1, characterized in that, In step 3), the continuous physiological nodes of the preset specific transmission path can also be one of the following: foot input → foot → ankle → calf → knee joint, foot input → foot → calf → knee joint → thigh → pelvis and buttocks, foot input → foot → calf → knee joint → thigh → pelvis and buttocks → abdomen → chest, foot input → foot → calf → knee joint → thigh → pelvis and buttocks → abdomen → chest → shoulder, or foot input → foot → calf → knee joint → thigh → pelvis and buttocks → abdomen → chest → shoulder → upper arm → elbow → forearm → wrist and hand.

7. The method for obtaining the propagation law of blasting vibration surface waves along the human body surface according to claim 1, characterized in that, In step 4), the step-by-step reduction prediction formula is used to quantitatively characterize the vibration energy loss pattern between adjacent physiological nodes due to differences in anatomical structure. Its mathematical expression is as follows: PPV n = η n ∙ PPV n-1 ； In the formula, PPV n The predicted peak vibration velocity of the nth physiological node on the specific propagation path; PPV n-1 Let n be the predicted peak vibration velocity of the (n-1)th physiological node on this path, and when n=1, v n-1 It is equal to the initial reference value of the input terminal; η n It is the stepwise vibration transmission reduction factor between the (n-1)th node and the nth node, reflecting the buffering and viscoelastic damping characteristics of specific physiological structures in the human body.

8. The method for obtaining the propagation law of blasting vibration surface waves along the human body surface according to claim 6, characterized in that, The step-by-step vibration transmission reduction factor η n The simulation node dataset, which includes various combinations of single-segment maximum charge and detonation center distance, was extracted and obtained by multivariate nonlinear regression fitting using the least squares method.

Images