Method and device for analyzing impact-explosion time series loading experimental data

Abstract

The application provides an analysis method and device for impact-explosion time sequence loading experiment data, comprising: obtaining experiment data; determining a stress-strain curve of a brittle medium rod material under a time sequence strong dynamic load in an impact-explosion experiment process according to the experiment data; establishing a propagation attenuation model by using the experiment data and the stress-strain curve, and determining stress waveform information of the time sequence strong dynamic load after the time sequence strong dynamic load propagates to the end of the brittle medium rod with a hole through the propagation attenuation model; and combining the final breaking form of the hole end and the incubation characteristic time criterion, the utilization and dissipation law of stress wave energy in the dynamic brittle fracture process can be analyzed. The application can provide a propagation attenuation model for stress waveform analysis of the time sequence strong dynamic load, so that the utilization and dissipation law of composite stress wave energy in the dynamic brittle fracture in the impact-explosion time sequence loading experiment can be reliably analyzed on the basis of the propagation attenuation model.

Description

Technical Field

[0001] The present invention relates to the technical field of impact / explosion testing, and in particular to an analysis method and device for impact-explosion time sequence loading experimental data. Background Art

[0002] During penetration, a penetrating weapon forms a crushing zone and an adjacent fracture zone in the near-trajectory region. This penetration generates a low-amplitude, long-duration penetration shock wave in the solid medium. Subsequently, the explosive within the projectile detonates, further generating a high-amplitude, short-duration explosion shock wave. The time interval between these two stress waves is Δt. After spatial attenuation and damping, the stress wave acts on the underground structure. The pre-disturbance of the penetration shock wave will affect the subsequent destructive effects of the explosion shock wave on brittle media such as rock and concrete, as well as underground structures. However, existing technologies not only struggle to reliably analyze the damage and failure patterns of brittle media under strong dynamic loading in a shock-explosion sequence, but also struggle to effectively analyze the energy utilization and dissipation patterns of the composite stress waves generated during shock-explosion sequence loading experiments. Summary of the Invention

[0003] In view of this, the purpose of the present invention is to provide a method and device for analyzing impact-explosion time-sequence loading experimental data, which can provide a propagation attenuation model for analyzing the action law of time-sequence strong dynamic loads, so that on the basis of this propagation attenuation model, the utilization and dissipation law of the composite stress wave energy when the time-sequence strong dynamic load generated during the impact-explosion experiment propagates to the end of the brittle medium rod with a hole and causes dynamic brittle fracture can be reliably analyzed.

[0004] In a first aspect, an embodiment of the present invention provides a method for analyzing impact-explosion time-sequence loading experimental data, comprising: acquiring experimental data; wherein the experimental data is generated during an impact-explosion experiment in an impact-explosion time-sequence loading experimental device; determining, based on the experimental data, a stress-strain curve of a brittle dielectric rod material under time-sequence strong dynamic load during the impact-explosion experiment; establishing a propagation attenuation model of time-sequence strong dynamic load in a brittle medium using the experimental data and the stress-strain curve, and determining, through the propagation attenuation model, the utilization and dissipation law of the composite stress wave energy when the time-sequence strong dynamic load propagates to the end of a brittle dielectric rod with a hole, causing dynamic brittle fracture.

[0005] In an embodiment, the impact-explosion time sequence loading experimental device is provided with a brittle medium rod for transmitting a composite stress wave caused by impact-explosion time sequence loading, and the brittle medium rod is provided with a speckle area, and the experimental data includes speckle area deformation image data; the axial displacement field distribution and the axial velocity field distribution of the surface mass point of the brittle medium rod are determined according to the experimental data; the axial strain-time curve of the brittle medium rod caused by the time sequence strong dynamic load generated in the impact-explosion experiment is determined according to the axial displacement field distribution, and a stress wave mass conservation equation in a Lagrangian coordinate system; the axial stress-time curve of the brittle medium rod caused by the time sequence strong dynamic load generated in the impact-explosion experiment is determined according to the axial displacement field distribution, a stress wave momentum conservation equation in the Lagrangian coordinate system and path line method integration; and the stress-strain curve of the brittle medium material under the time sequence strong dynamic load is obtained based on the strain-time curve and the stress-time curve.

[0006] In an embodiment, the step of determining the axial displacement field distribution and the axial velocity field distribution of the surface mass point of the brittle medium rod based on the speckle area deformation image includes: performing image analysis processing on the speckle area deformation image data by using a DIC technology to determine the axial displacement field distribution and the axial velocity field distribution of the surface mass point of the brittle medium rod under the time sequence strong dynamic load.

[0007] In an embodiment, the brittle medium rod is provided with an ultra-dynamic strain gauge, and the experimental data includes strain gauge measurement results; and the propagation and attenuation model of the time sequence strong dynamic load in the brittle medium is established based on the experimental data and the stress-strain curve, which includes: extracting waveform time-frequency features of the ultra-dynamic strain gauge measurement results; determining change data of the composite stress wave waveform time-frequency features according to the propagation distance of the time sequence strong dynamic load; wherein the change data includes one or more of propagation velocity, transmission reflection coefficient, frequency spectrum characteristics, frequency band energy and fractal characteristics; and the propagation and attenuation model of the time sequence strong dynamic load is established based on the damping attenuation parameters of the brittle medium rod, the change data and the stress-strain curve.

[0008] In one embodiment, the impact-explosion sequential loading experimental device is provided with a brittle dielectric rod, one end of which is provided with a hole to be broken; the utilization and dissipation law of the composite stress wave energy when the dynamic brittle fracture is caused by the time-sequential strong dynamic load propagating to the end of the brittle dielectric rod with the hole is determined by the propagation attenuation model, including: calculating the stress time history information acting on the hole to be broken by the propagation attenuation model; obtaining a hole fracture image corresponding to the hole to be broken, and determining the effective stress wave energy required in the fracture process of the hole to be broken based on the hole fracture image; based on the effective stress wave energy for brittle fracture, performing energy utilization analysis on the stress time history information using a breeding characteristic time criterion, the stress time history information being the stress time history information in the brittle medium after the time-sequential strong dynamic load propagates to the end of the brittle dielectric rod with the hole obtained by the propagation attenuation model, thereby obtaining the distribution and proportion law of the stress wave energy consumed to cause dynamic brittle fracture in the time domain and frequency domain of the composite stress waveform; wherein the utilization and dissipation law includes distribution and proportion laws.

[0009] In a second aspect, an embodiment of the present invention further provides an analysis device for shock-explosion time-series loading experimental data, comprising: a data acquisition module for acquiring experimental data; wherein the experimental data is generated during a shock-explosion experiment in a shock-explosion time-series loading experimental device; a material mechanical property determination module for determining, based on the experimental data, a stress-strain curve of a brittle dielectric rod material under time-series strong dynamic loads during the shock-explosion experiment; a law determination module for establishing a propagation attenuation model of time-series strong dynamic loads in brittle media using the experimental data and the stress-strain curve, and determining, through the propagation attenuation model, a utilization and dissipation law of composite stress wave energy when dynamic brittle fracture is caused after the time-series strong dynamic load propagates to the end of the brittle dielectric rod with a hole.

[0010] In an embodiment, the impact-explosion time sequence loading experiment device is provided with a brittle medium rod for transmitting a composite stress wave caused by impact-explosion time sequence loading, and the brittle medium rod is provided with a speckle area, and the experiment data includes speckle area deformation image data; the material mechanical property determination module is further configured to: determine an axial displacement field distribution and an axial velocity field distribution of a surface particle of the brittle medium rod based on the speckle area deformation image; determine a strain time history curve of the brittle medium under the time sequence strong dynamic load generated in the impact-explosion experiment according to the axial displacement field distribution and a stress wave mass conservation equation in a Lagrangian coordinate system; determine a stress time history curve of the brittle medium under the time sequence strong dynamic load generated in the impact-explosion experiment according to the axial displacement field distribution, the stress wave momentum conservation equation in the Lagrangian coordinate system and path line method integration; and obtain a stress-strain curve of the brittle medium material under the time sequence strong dynamic load based on the strain time history curve and the stress time history curve.

[0011] In an embodiment, the material mechanical property determination module is further configured to: perform image analysis and processing on the speckle area deformation image data by using a DIC (Digital Image Correlation) technology, and determine the axial displacement field distribution and the axial velocity field distribution of the surface particle of the brittle medium rod under the time sequence strong dynamic load.

[0012] In a third aspect, an electronic device is provided, including a processor and a memory, the memory storing computer executable instructions capable of being executed by the processor, and the processor executes the computer executable instructions to implement the method according to any one of the first aspect.

[0013] In a fourth aspect, a computer readable storage medium is provided, the computer readable storage medium storing computer executable instructions, and the computer executable instructions, when invoked and executed by a processor, cause the processor to implement the method according to any one of the first aspect.

[0014] The analysis method and device for shock-explosion time-sequence loading experimental data provided by the embodiment of the present invention first obtain experimental data generated during the shock-explosion experiment in the shock-explosion time-sequence loading experimental device, then determine the stress-strain curve of the brittle dielectric rod material under the action of a time-sequence strong dynamic load during the shock-explosion experiment based on the experimental data, and finally establish a time-sequence strong dynamic load propagation attenuation model based on the experimental data and the stress-strain curve. The propagation attenuation model is used to determine the utilization and dissipation law of the composite stress wave energy when the time-sequence strong dynamic load propagates to the end of the brittle dielectric rod with a hole and causes dynamic brittle fracture. The above method can obtain the experimental data generated during the shock-explosion experiment in the shock-explosion time-sequence loading experimental device, and determine the corresponding stress time history curve based on the experimental data, thereby obtaining a propagation attenuation model that can be used for stress waveform analysis of the time-sequence strong dynamic load, so that the utilization and dissipation law of the composite stress wave energy when the composite stress wave generated during the shock-explosion experiment causes dynamic brittle fracture can be reliably analyzed based on the propagation attenuation model.

[0015] Other features and advantages of the present invention will be described in the following description, and in part will become apparent from the description, or understood by practicing the present invention. The purposes and other advantages of the present invention are realized and obtained by the structures particularly pointed out in the description, claims and drawings.

[0016] In order to make the above-mentioned objects, features and advantages of the present invention more obvious and easy to understand, preferred embodiments are given below and described in detail with reference to the accompanying drawings. BRIEF DESCRIPTION OF THE DRAWINGS

[0017] In order to more clearly illustrate the specific embodiments of the present invention or the technical solutions in the prior art, the following briefly introduces the drawings required for use in the specific embodiments or the description of the prior art. Obviously, the drawings described below are 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.

[0018] Figure 1 A schematic flow chart of a method for analyzing shock-explosion time-series loading experimental data provided by an embodiment of the present invention;

[0019] Figure 2 A schematic structural diagram of an impact-explosion sequence loading experimental device provided by an embodiment of the present invention;

[0020] Figure 3 A schematic structural diagram of an analysis device for shock-explosion time-series loading experimental data provided by an embodiment of the present invention;

[0021] Figure 4A structural schematic diagram of an electronic device is provided for the embodiment of the present application.

[0022] Icon: 1-laser velocimeter; 2-impact-explosion time sequence loading bin; 3.1-impact block; 3.2-impact block; 4-explosive; 5-plug; 6-detonator; 7-loading bin charging hole; 8-wave shaper; 9-frangible medium rod charging hole; 10-frangible medium rod; 11-speckle area; 12-LED direct current light source; 13-high-speed camera; 14-ultra-dynamic strain gauge; 15-to-be-broken empty hole; 16-stress absorber; 100-electronic device; 40-processor; 41-memory; 42-bus; 43-communication interface. DETAILED DESCRIPTION

[0023] To make the purpose, technical scheme and advantages of the embodiment of the present application more clear, the technical scheme of the present application will be described clearly and completely below in combination with the embodiments. Obviously, the described embodiments are part of the embodiments of the present application, rather than all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative labor fall within the protection scope of the present application.

[0024] At present, the existing technology cannot reliably quantitatively analyze the deformation and damage law of the brittle medium under the impact-explosion time sequence strong dynamic load. Based on this, the present application provides a propagation attenuation model for time sequence strong dynamic load stress waveform analysis, so that the energy utilization and dissipation law of the composite stress wave generated in the impact-explosion experiment process can be reliably analyzed based on the propagation attenuation model.

[0025] To facilitate the understanding of the present embodiment, first, a method for analyzing impact-explosion time sequence loading experimental data disclosed by the present embodiment will be described in detail, referring to the flowchart of the method for analyzing impact-explosion time sequence loading experimental data shown in Figure 1 The method mainly includes the following steps S102 to S106:

[0026] Step S102, obtaining experimental data. The experimental data is generated in the process of impact-explosion experiment in the impact-explosion time sequence loading experimental device, and can include speckle area deformation image data and ultra-dynamic strain gauge measurement results. In actual application, a frangible medium rod is arranged in the impact-explosion time sequence loading experimental device, which is used to transmit the composite stress wave caused by impact-explosion time sequence loading. The frangible medium rod is provided with a speckle area and an ultra-dynamic strain gauge, so that the image at the speckle area can be collected by an image collection device (such as a high-speed camera), and the image is the speckle area deformation image data. In addition, the measurement result of the ultra-dynamic strain gauge can be directly read.

[0027] Step S104, according to the experimental data, determine the stress-strain curve of the brittle medium rod material under the time sequence strong dynamic load in the process of impact-explosion experiment. Wherein, the stress-strain curve is used to characterize the mapping relationship between stress and strain under the dynamic loading condition of brittle medium. In an embodiment, the speckle area deformation image data in the experimental data can be processed to determine the axial displacement field distribution and axial velocity field distribution of the surface particles of the brittle medium rod characterized by the speckle area deformation image data, so as to determine the strain time history curve and the stress time history curve respectively on the basis of the axial displacement field distribution and the axial velocity field distribution, and finally the stress-strain curve can be obtained based on the strain time history curve and the stress time history curve. The strain time history curve can be verified by comparing with the test results of the super dynamic strain gauge.

[0028] Step S106, based on the experimental data and the stress-strain curve, establish a time sequence strong dynamic load propagation attenuation model, and determine the utilization and dissipation law of the composite stress wave energy caused by the dynamic brittle fracture when the composite stress wave propagates to the end of the brittle medium rod with a hole. In an embodiment, the propagation attenuation model can be established based on the measurement results of the super dynamic strain gauge in the experimental data and the stress-strain curve, wherein the propagation attenuation model can estimate the stress wave form acting on (entering) the to-be-broken hole, and further analyze the time-frequency distribution characteristics of the composite stress wave form obtained by time-frequency analysis. In addition, the size distribution of the hole fragments is obtained by screening, and the final effective energy for dynamic brittle fracture is obtained, based on which the utilization and dissipation law of the composite stress wave frequency band energy caused by the time sequence strong dynamic load (referred to as the utilization and dissipation law of stress wave energy) can be further analyzed.

[0029] The impact-explosion time sequence loading experiment data analysis method provided by the embodiment of the application can obtain the experimental data generated in the impact-explosion experiment in the impact-explosion time sequence loading experiment device, and determine the corresponding stress time history curve based on the experimental data, so as to obtain the propagation attenuation model which can be used for time sequence strong dynamic load stress wave form analysis, so that the utilization and dissipation law of the composite stress wave energy in the process of dynamic brittle destruction of the medium can be reliably analyzed based on the propagation attenuation model.

[0030] In order to facilitate the understanding of the foregoing embodiments, the embodiment of the application exemplarily provides an impact-explosion time sequence loading experiment device, as shown in Figure 2A structure diagram of an impact-explosion time sequence loading experimental device is shown, and the impact-explosion time sequence loading experimental device comprises a laser speedometer 1, an impact-explosion time sequence loading bin 2, a hitting block 3.1, a hitting block 3.2, an explosive 4, a plugging device 5, a detonator 6, a loading bin charging hole 7, a wave shaper 8, a brittle medium rod charging hole 9, a brittle medium rod 10, a speckle area 11, an LED (light-emitting diode) direct current light source 12, a high-speed camera 13, an ultra-dynamic strain gauge 14, a to-be-broken empty hole 15 and a stress absorber 16.

[0031] In actual application, the explosive 4 is loaded in the empty hole spliced by the hitting block 3.1 and the hitting block 3.2 inside the impact-explosion time sequence loading bin 2 and the brittle medium rod charging hole 9 at the left end of the brittle medium rod 10 as a load source, and the detonator 6 in the loading bin charging hole 7 of the impact-explosion time sequence loading bin 2 is used for detonation. First, the explosive in the spliced hitting block is detonated, and the hitting block 3.2 is driven to collide with the left end of the brittle medium rod 10, at this time, according to the law of conservation of momentum, the hitting block 3.1 flies out from the left end of the impact-explosion time sequence loading bin 2 at the same speed, and the movement speed of the hitting block is measured by using the laser speedometer 1. Since the density of the hitting block is much greater than that of the brittle medium, the speed of the hitting block is not 0 after the collision, and therefore a continuous impact pulse is generated at the end of the brittle medium rod 10, and the wave shaper 8 can make the impact pulse more smooth. After the impact pulse enters the brittle medium rod 10, the explosive in the brittle medium rod charging hole 9 is detonated, and an explosion pulse is formed in the brittle medium rod 10. The two stress pulses gradually stabilize and propagate to the right end in the brittle medium rod 10, and the propagation and attenuation process can be recorded by the high-speed camera 13, the LED direct current light source 12 and the ultra-dynamic strain gauge 14, and the combined action of the two continuous stress pulses will cause the deformation and brittle fracture of the to-be-broken empty hole 15. The residual stress wave propagates to the right end of the brittle medium rod 10, and the stress absorber 16 made of the same brittle medium, and the influence of the residual stress wave on the empty hole will be rapidly reduced due to the space attenuation.

[0032] On the basis of the foregoing device, the embodiment of the present application provides an implementation mode for determining the stress-strain curve of the brittle medium material under the time sequence strong dynamic load in the impact-explosion experiment process according to experimental data, and the implementation mode comprises the following steps 1 to 4:

[0033] Step 1: Determine the axial displacement field distribution and axial velocity field distribution of the particles on the surface of the brittle dielectric rod based on the deformation image of the speckle area. In one embodiment, the DIC (Digital Image Correlation) technique can be used to perform image analysis processing on the deformation image data of the speckle area to determine the axial displacement field distribution and axial velocity field distribution of the particles on the surface of the brittle dielectric rod under time-sequential strong dynamic loads. In practical applications, the DIC technique can be used to analyze the propagation process of time-sequential strong dynamic loads in the speckle area captured by a high-speed camera to obtain the axial displacement field distribution and axial velocity field distribution caused by the propagation of stress waves in the brittle dielectric rod. Optionally, the axial displacement field distribution and axial velocity field distribution can be presented in the form of cloud maps.

[0034] Step 2: Based on the axial displacement field distribution and the stress wave mass conservation equation in the Lagrangian coordinate system, determine the strain time history curve corresponding to the time-series strong dynamic load generated during the shock-explosion experiment that disturbs the brittle medium. In the Lagrangian coordinate system, the mass conservation equation for one-dimensional stress waves is described as follows:

[0035] Where v is the particle velocity, ε is the strain, X is the Lagrangian coordinate position, and t is the time. In practical applications, the particle velocity waveform v(X) at different X points obtained by DIC technology analysis is i , t), find its first-order partial derivative with respect to the Lagrangian coordinate position X From this we can obtain the mass conservation equation Then use the initial condition at time zero [t=0,ε(X i , 0)=0], the strain time history curve ε(X i , t).

[0036] Step 3: Based on the axial displacement field distribution, the momentum conservation equation for stress waves in the Lagrangian coordinate system, and the path line method integral, determine the stress time history curve corresponding to the time-series strong dynamic load generated during the shock-blast experiment that disturbs the brittle medium. In the Lagrangian coordinate system, the momentum conservation equation for the one-dimensional stress wave is described as follows:

[0037] Where σ is stress and ρ is density. In practical applications, the momentum conservation equation is further Convert to At this time, the total differential relationship of stress along the path line P is as follows:

[0038]

[0039] Substituting the above mass conservation equation into the total differential relationship of stress along the path line P, we obtain the following formula:

[0040]

[0041] On this basis, the first path that satisfies the initial condition at time zero [t = 0, σ = (dσ / dX)| P1 =0] Determine the path line P at time zero Then the stress σ on the next path line can be determined by numerical integration, and the stress time history curve σ(X i , t).

[0042] The embodiment of the present invention uses the Lagrange inverse analysis method, one-dimensional stress wave propagation theory (including the above-mentioned mass conservation equation and momentum conservation equation) and the total differential relationship of stress along the path line P to obtain the strain time history curve and the stress time history curve respectively.

[0043] Step 4: Based on the strain time-history curve and the stress time-history curve, a stress-strain curve of the brittle dielectric material under a time-sequential strong dynamic load is obtained. In one embodiment, the strain time-history curve can be used to compare the test results of an ultra-dynamic strain gauge. Alternatively, the stress-strain curve of the brittle dielectric material under a time-sequential strong dynamic load can be obtained by eliminating the time parameter from the strain time-history curve and the stress time-history curve.

[0044] Based on the aforementioned device, an embodiment of the present invention provides an implementation method for establishing a propagation attenuation model based on experimental data and stress-strain curves. Specifically:

[0045] (1) Extracting the time-frequency characteristics of the stress time history curve waveform. In one embodiment, the time-frequency characteristics of the strain gauge measurement results can be analyzed using algorithms such as short-time Fourier transform, continuous wavelet transform, or Hilbert-Huang transform to obtain the waveform time-frequency characteristics.

[0046] (2) Determine the variation data of the time-frequency characteristics of the composite stress waveform based on the propagation distance of the time-series strong dynamic load. The variation data includes one or more of propagation velocity, transmissivity, spectrum characteristics, frequency band energy, and fractal characteristics. In one embodiment, the variation pattern of the waveform time-frequency information with propagation distance (i.e., the aforementioned variation data) can be obtained based on the stress wave propagation velocity, transmissivity, spectrum characteristics, frequency band energy, and fractal characteristics at different measuring points.

[0047] (3) A propagation attenuation model is established based on the damping attenuation parameters, variation data, and stress-strain curve of the brittle medium rod. In one embodiment, the propagation attenuation model for time-series strong dynamic loads can be obtained by fitting based on the classical one-dimensional wave equation, taking into account the damping attenuation effect of the brittle medium and combining it with the stress-strain curve of the brittle medium under the aforementioned time-series strong dynamic load. The propagation attenuation model is as follows:

[0048] Where ρ is density, u is displacement, t is time, E is elastic modulus, and η is damping coefficient.

[0049] Based on the analysis method of the impact-explosion time-sequence loading experimental data provided in the aforementioned embodiment, an embodiment of the present invention further provides an implementation method for determining, through the propagation attenuation model, the corresponding energy utilization and dissipation law when the time-sequence strong dynamic load propagates to the end of a brittle dielectric rod with a hole, causing dynamic brittle fracture. Specifically, the impact-explosion time-sequence loading experimental apparatus provided in the aforementioned embodiment is provided with a brittle dielectric rod, one end of which is provided with a hole to be broken. Based on this, please refer to the following steps a to d:

[0050] Step a: Calculate the time history information of the stress acting on the pore to be broken by using the propagation attenuation model. Optionally, the time history information of the stress waveform entering or acting on the pore area to be broken can be obtained by using the propagation attenuation model.

[0051] Step b, obtaining a hole fracture image corresponding to the hole to be broken. In one embodiment, the hole fracture image can be obtained by capturing the progressive destruction process of the hole to be broken at the right end of the brittle dielectric rod using a high-speed camera.

[0052] Step c: Determine the effective fracture energy during the fracture process of the pore to be broken based on the pore fracture image. In one embodiment, the broken fragments can be collected, and a series of standard sieves with appropriate mesh sizes are selected according to ISO (the International Organization for Standardization) 3310 or ASTM (American Society for Testing and Materials International) E11 to perform size and mass statistics on the brittle medium fragments. Combined with the time-series strong dynamic load propagation attenuation model provided in the above embodiment, the stress-strain curve of the incident wave at the left end of the pore to be broken is obtained, and the influence of the composite stress wave induced by different impact-explosion loading sequences on dynamic brittle fracture is obtained. Based on the newly added fracture surface of the fragments, the effective stress wave energy consumed by the dynamic brittle fracture is obtained.

[0053] In step d, based on the effective stress wave energy for brittle fracture, the incubation characteristic time criterion is used to analyze the energy utilization of the stress time history information. The distribution and proportion of the energy required for dynamic brittle fracture in the time and frequency domains of the composite stress wave caused by the time-sequential strong dynamic load propagating to the end of the brittle dielectric rod with a hole are obtained. The incubation characteristic time criterion (ICT) is as follows:

[0054] Among them, F * is the dimensionless dynamic brittle fracture criterion function, t is time, τ is the incubation characteristic time, σ0 is the quasi-static strength, and σ(t') is the equivalent stress history corresponding to σ0.

[0055] In one embodiment, the incubation characteristic time criterion can be further utilized to study the utilization and dissipation rules of waveform frequency band energy due to dynamic brittle fracture by considering the time-frequency characteristics of the composite stress history waveform, so as to obtain the energy utilization characteristics of time-series strong dynamic loads under different loading intervals.

[0056] In summary, the analysis method of the impact-explosion time sequence loading experimental data provided by the embodiment of the present invention has at least the following characteristics: (1) the axial displacement and velocity field on the brittle medium rod can be obtained by cooperating with the digital image processing technology (DIC), and the Lagrangian inverse analysis method is used to process the velocity field data to obtain the stress time history curve, the strain time history curve and the stress-strain relationship of the brittle medium under time sequence loading, respectively. The strain time history curve can be further compared and verified with the test results of the ultra-dynamic strain gauge, which can ensure the reliability of the test data; in addition, the time-frequency analysis of the stress time history curve, combined with the stress-strain relationship of the brittle medium under time sequence loading, can establish the propagation and attenuation model of the composite stress wave under the impact-explosion time sequence coupling loading based on the mass conservation equation and momentum conservation equation of the one-dimensional stress wave; (2) high-speed camera monitoring can capture the progressive destruction process of the pore to be broken, and the collected fragments are in accordance with ISO 3310 or ASTM Screening of E11 can obtain the effective energy used for dynamic brittle fracture during the crushing process. The stress time history information acting on the empty hole area to be crushed is calculated based on the stress wave propagation attenuation model. Then, the influence of the waveform time-frequency characteristics on dynamic brittle fracture can be obtained by analyzing the input stress history using the incubation characteristic time criterion.

[0057] For the analysis method of the shock-explosion time sequence loading experimental data provided in the above embodiment, the present invention also provides an analysis device for the shock-explosion time sequence loading experimental data, see Figure 3 The schematic diagram of the structure of an analysis device for shock-explosion time sequence loading experimental data is shown. The device mainly includes the following parts:

[0058] The data acquisition module 302 is used to acquire experimental data; wherein the experimental data is generated during the shock-explosion experiment in the shock-explosion sequential loading experimental device;

[0059] The material mechanical property determination module 304 is used to determine the stress-strain curve of the brittle dielectric rod material under time-series strong dynamic load during the impact-explosion experiment based on the experimental data;

[0060] The rule determination module 306 is configured to establish a propagation attenuation model of the time sequence strong dynamic load in the brittle medium by using the experimental data and the stress-strain curve, and determine the utilization and dissipation rule of the dynamic brittle fracture time composite stress wave energy caused by the propagation of the time sequence strong dynamic load to the end of the brittle medium rod with the propagation attenuation model.

[0061] The analysis device for impact-explosion time sequence loading experimental data provided by the embodiment of the present application can obtain experimental data generated in the impact-explosion experiment in the impact-explosion time sequence loading experimental device, determine the corresponding stress time history curve based on the experimental data, and obtain the propagation attenuation model which can be used to analyze the stress wave attenuation rule, so that the utilization and dissipation rule of the composite stress wave energy generated in the impact-explosion experiment can be reliably analyzed based on the propagation attenuation model.

[0062] In an embodiment, the impact-explosion time sequence loading experimental device is provided with a brittle medium rod for transmitting the composite stress wave caused by the impact-explosion time sequence loading, and the brittle medium rod is provided with a speckle area, and the experimental data includes speckle area deformation image data; the material mechanical property determination module 304 is further configured to determine the axial displacement field distribution and the axial velocity field distribution of the surface particles of the brittle medium rod based on the speckle area deformation image, determine the axial strain time history curve of the brittle medium rod caused by the time sequence strong dynamic load generated in the impact-explosion experiment according to the axial displacement field distribution and the stress wave mass conservation equation in the Lagrangian coordinate system, determine the corresponding stress time history curve of the brittle medium deformation caused by the time sequence strong dynamic load generated in the impact-explosion experiment according to the axial displacement field distribution, the stress wave momentum conservation equation in the Lagrangian coordinate system and the path line method integral, and obtain the stress-strain curve of the brittle medium material under the time sequence strong dynamic load based on the strain time history curve and the stress time history curve.

[0063] In an embodiment, the material mechanical property determination module 304 is further configured to perform image analysis and processing on the speckle area deformation image data by using the DIC technology, and determine the axial displacement field distribution and the axial velocity field distribution of the surface particles of the brittle medium rod under the time sequence strong dynamic load.

[0064] In an embodiment, the brittle medium rod is provided with an ultra-dynamic strain gauge, and the experimental data includes the strain gauge measurement results; the rule determination module 306 is further configured to extract the waveform time-frequency features of the ultra-dynamic strain gauge measurement results, determine the change data of the composite stress wave form time-frequency features according to the propagation distance of the impact-explosion strong dynamic load, wherein the change data includes one or more of the propagation velocity, the transmission reflection coefficient, the frequency spectrum characteristics, the frequency band energy and the fractal characteristics, and establish the propagation attenuation model of the time sequence strong dynamic load based on the damping attenuation parameters of the brittle medium rod, the change data and the stress-strain curve.

[0065] In one embodiment, the impact-explosion sequential loading experimental device is provided with a brittle medium rod, and a hole to be broken is provided at one end of the brittle medium rod; the law determination module 306 is further used to: calculate the stress time history information acting on the hole to be broken through a time-series strong dynamic load propagation attenuation model; obtain a hole fracture image corresponding to the hole to be broken; based on the hole fracture process image, combined with the screening results of the brittle medium fragments, determine the effective stress wave energy consumed in the fracture process of the hole to be broken; based on the effective stress wave energy for brittle fracture, use the incubation characteristic time criterion to perform energy utilization analysis on the stress time history information, the stress time history information is the stress time history information in the brittle medium after the time-series strong dynamic load is propagated to the end of the brittle medium rod with the hole obtained through the propagation attenuation model, and then obtain the distribution and proportion law of the stress wave energy consumed to cause dynamic brittle fracture in the time domain and frequency domain of the composite stress waveform; wherein, the dissipation law includes the distribution and proportion law.

[0066] The device provided in the embodiment of the present invention has the same implementation principle and technical effects as those in the aforementioned method embodiment. For the sake of brief description, for matters not mentioned in the device embodiment, reference can be made to the corresponding content in the aforementioned method embodiment.

[0067] An embodiment of the present invention provides an electronic device. Specifically, the electronic device includes a processor and a storage device. The storage device stores a computer program, and when the computer program is executed by the processor, it executes the method described in any one of the above-mentioned embodiments.

[0068] Figure 4 This is a structural diagram of an electronic device provided in an embodiment of the present invention. The electronic device 100 includes: a processor 40, a memory 41, a bus 42 and a communication interface 43. The processor 40, the communication interface 43 and the memory 41 are connected via the bus 42; the processor 40 is used to execute an executable module stored in the memory 41, such as a computer program.

[0069] The memory 41 may include high-speed random access memory (RAM) and may also include non-volatile memory, such as at least one disk storage. The system network element communicates with at least one other network element via at least one communication interface 43 (which may be wired or wireless), and may utilize the Internet, a wide area network, a local area network, a metropolitan area network, or the like.

[0070] The bus 42 may be an ISA bus, a PCI bus, or an EISA bus. The bus may be divided into an address bus, a data bus, a control bus, and the like. For ease of representation, Figure 4 Only one bidirectional arrow is used in the diagram, but this does not mean that there is only one bus or one type of bus.

[0071] The memory 41 is configured to store a program, and the processor 40 executes the program after receiving an execution instruction. The method executed by the device for defining a flow process according to any of the foregoing embodiments of the present application can be applied to the processor 40 or implemented by the processor 40.

[0072] The processor 40 can be an integrated circuit chip having a processing capability of signals. In the implementation process, each step of the foregoing method can be completed by an integrated logic circuit of hardware in the processor 40 or an instruction in the form of software. The foregoing processor 40 can be a general-purpose processor, including a central processing unit (CPU), a network processor (NP), and the like; or can be a digital signal processor (DSP), an application specific integrated circuit (ASIC), a field-programmable gate array (FPGA) or other programmable logic device, a discrete gate or transistor logic device, a discrete hardware component. Each method, step and logic block diagram disclosed in the embodiments of the present application can be implemented or executed. The general-purpose processor can be a microprocessor or the processor can also be any conventional processor. The steps of the method disclosed in combination with the embodiments of the present application can be directly embodied as a hardware code processor for execution, or a combination of hardware and software modules in the code processor for execution. The software module can be located in a random access memory, a flash memory, a read-only memory, a programmable read-only memory or an electrically erasable programmable memory, a register, or other mature storage media in the art. The storage medium is located in the memory 41, and the processor 40 reads the information in the memory 41 and combines the hardware to complete the steps of the foregoing method.

[0073] The computer program product of the readable storage medium provided by the embodiments of the present application includes a computer readable storage medium storing a program code, and the program code includes instructions for executing the method described in the foregoing method embodiments. For specific implementation, reference can be made to the foregoing method embodiments, which will not be described here.

[0074] If the functions are realized in the form of software function units and sold or used as independent products, they can be stored in a computer readable storage medium. Based on this understanding, the technical solutions of the present application or the parts of the present application that essentially contribute to the prior art or the parts of the technical solutions can be embodied in the form of software products. The computer software product is stored in a storage medium and includes instructions for making a computer device (which can be a personal computer, a server, or a network device, etc.) execute all or part of the steps of the method described in the embodiments of the present application. The aforementioned storage medium includes: a U disk, a mobile hard disk, a read-only memory (ROM, Read-Only Memory), a random access memory (RAM, Random Access Memory), a magnetic disk or an optical disk, and various media that can store program codes.

[0075] Finally, it should be noted that: the above-described embodiments are only specific embodiments of the present application, which are used to illustrate the technical solutions of the present application, but not to limit them. The protection scope of the present application is not limited thereto. Although the present application has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that any person skilled in the art can make modifications or easily think of changes to the technical solutions described in the foregoing embodiments within the technical scope disclosed by the present application, or make equivalent replacements to some technical features. The modifications, changes or replacements do not make the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the embodiments of the present application, and should be covered within the protection scope of the present application. Therefore, the protection scope of the present application should be subject to the protection scope of the claims.

Claims

1. A method for analyzing shock-explosion time series loading experimental data, characterized in that: include: Acquire experimental data; wherein, the experimental data is generated during a shock-explosion experiment in a shock-explosion sequence loading experimental device, the shock-explosion sequence loading experimental device comprising a shock-explosion sequence loading chamber, two impact blocks, a brittle medium rod charge hole, and a brittle medium rod, wherein the empty hole formed by the two impact blocks inside the shock-explosion sequence loading chamber and the brittle medium rod charge hole at the left end of the brittle medium rod are both loaded with explosives, and the explosives loaded in the empty hole formed by the two impact blocks are detonated, one of the impact blocks is driven to collide with the left end of the brittle medium rod, and the impact pulse enters the brittle medium rod and detonates the explosive loaded in the charge hole of the brittle medium rod, forming an explosion pulse in the brittle medium rod, the brittle medium rod is used to transmit the composite stress wave caused by the shock-explosion sequence loading, and one end of the brittle medium rod is provided with a hole to be broken, and the composite stress wave will cause deformation and brittle fracture of the hole to be broken; Determining, based on the experimental data, a stress-strain curve of the brittle dielectric rod material under a time-series strong dynamic load during the impact-explosion experiment; A propagation attenuation model of a time-sequential strong dynamic load in a brittle medium is established using the experimental data and the stress-strain curve, and a utilization and dissipation law of the composite stress wave energy when the time-sequential strong dynamic load propagates to the end of the brittle medium rod with a hole and causes dynamic brittle fracture is determined using the propagation attenuation model. The utilization and dissipation law includes distribution and proportion laws. The method of determining the utilization and dissipation law of the composite stress wave energy when the time-series strong dynamic load propagates to the end of the brittle dielectric rod with a hole and causes dynamic brittle fracture by the propagation attenuation model includes: Calculating the stress time history information acting on the empty hole to be broken by using the time-series strong dynamic load propagation attenuation model; Capturing the progressive destruction process of the hole to be broken at the right end of the brittle dielectric rod by a high-speed camera to obtain a hole fracture image corresponding to the hole to be broken; Based on the image of the pore fracture process and in combination with the screening results of the brittle medium fragments, the effective stress wave energy consumed in the fracture process of the pore to be broken is determined; Based on the effective stress wave energy for brittle fracture, the energy utilization analysis of the stress time history information is performed using the incubation characteristic time criterion. The stress time history information is the stress time history information in the brittle medium after the time-series strong dynamic load propagates to the end of the brittle medium rod with a hole, obtained by the propagation attenuation model. The distribution and proportion of the stress wave energy consumed in causing dynamic brittle fracture in the time domain and frequency domain of the composite stress waveform are then obtained. The expression of the incubation characteristic time criterion is as follows: ;in, is the dimensionless dynamic brittle fracture criterion function, t is the time, is the characteristic time of gestation, is the quasi-static strength, For The corresponding equivalent stress history.

2. The method according to claim 1, characterized in that The brittle dielectric rod is provided with a speckle region, and the experimental data includes deformation image data of the speckle region; and determining, based on the experimental data, a stress-strain curve of the brittle dielectric rod material under a time-series strong dynamic load during the shock-explosion experiment, comprises: determining the axial displacement field distribution and the axial velocity field distribution of the particle points on the surface of the brittle dielectric rod based on the deformation image of the speckle region; Determine, based on the axial displacement field distribution and the stress wave mass conservation equation in the Lagrangian coordinate system, a time-history curve of the axial strain of the brittle dielectric rod caused by the sequential strong dynamic load generated during the shock-explosion experiment; Determine, based on the axial displacement field distribution, the stress wave quantity conservation equation in the Lagrangian coordinate system, and the path line method integral, a stress time history curve of the brittle dielectric rod caused by the sequential strong dynamic load generated during the shock-explosion experiment; Based on the strain time history curve and the stress time history curve, a stress-strain curve of the brittle dielectric material under a time-series strong dynamic load is obtained.

3. The method according to claim 2, characterized in that The step of determining the axial displacement field distribution and the axial velocity field distribution of the particles on the surface of the brittle dielectric rod based on the deformation image of the speckle area comprises: The DIC technique is used to perform image analysis processing on the deformation image data of the speckle area to determine the axial displacement field distribution and axial velocity field distribution of the particle point on the surface of the brittle dielectric rod under time-series strong dynamic load.

4. The method according to claim 2, characterized in that The brittle medium rod is provided with an ultra-dynamic strain gauge, and the experimental data includes strain gauge measurement results; and establishing a propagation attenuation model of a time-series strong dynamic load in a brittle medium based on the experimental data and the stress-strain curve includes: Extracting the waveform time-frequency characteristics of the ultra-dynamic strain gauge measurement result; Determining change data of the time-frequency characteristics of the composite stress waveform according to the propagation distance of the impact-explosion strong dynamic load; wherein the change data includes one or more of propagation velocity, transmittance and reflection coefficient, spectrum characteristics, frequency band energy and fractal characteristics; A time-series strong dynamic load propagation attenuation model is established based on the damping attenuation parameter of the brittle dielectric rod, the change data and the stress-strain curve.

5. An analysis device for shock-explosion time sequence loading experimental data, characterized in that: include: A data acquisition module is used to acquire experimental data; wherein, the experimental data is generated during a shock-explosion experiment in a shock-explosion sequence loading experimental device, the shock-explosion sequence loading experimental device comprising a shock-explosion sequence loading chamber, two impact blocks, a brittle medium rod charge hole and a brittle medium rod, the empty hole formed by the two impact blocks inside the shock-explosion sequence loading chamber and the brittle medium rod charge hole at the left end of the brittle medium rod are both loaded with explosives, the explosives loaded in the empty hole formed by the two impact blocks are detonated, one of the impact blocks is driven to collide with the left end of the brittle medium rod, and the impact pulse enters the brittle medium rod and detonates the explosive loaded in the brittle medium rod charge hole, forming an explosion pulse in the brittle medium rod, the brittle medium rod is used to transmit a composite stress wave caused by shock-explosion sequence loading, one end of the brittle medium rod is provided with a hole to be broken, and the composite stress wave will cause deformation and brittle fracture of the hole to be broken; a material mechanical property determination module, configured to determine, based on the experimental data, a stress-strain curve of the brittle dielectric rod material under a time-series strong dynamic load during the impact-explosion experiment; a law determination module for establishing a propagation attenuation model of a time-series strong dynamic load in a brittle medium using the experimental data and the stress-strain curve, and determining, through the propagation attenuation model, a law of utilization and dissipation of composite stress wave energy when dynamic brittle fracture is caused by the time-series strong dynamic load propagating to the end of the brittle medium rod with a hole; The rule determination module is specifically used for: Calculating the stress time history information acting on the empty hole to be broken by using the time-series strong dynamic load propagation attenuation model; Capturing the progressive destruction process of the hole to be broken at the right end of the brittle dielectric rod by a high-speed camera to obtain a hole fracture image corresponding to the hole to be broken; Based on the image of the pore fracture process and in combination with the screening results of the brittle medium fragments, the effective stress wave energy consumed in the fracture process of the pore to be broken is determined; Based on the effective stress wave energy for brittle fracture, the energy utilization analysis of the stress time history information is performed using the incubation characteristic time criterion. The stress time history information is the stress time history information in the brittle medium after the time-series strong dynamic load propagates to the end of the brittle medium rod with a hole, obtained by the propagation attenuation model. The distribution and proportion of the stress wave energy consumed in causing dynamic brittle fracture in the time domain and frequency domain of the composite stress waveform are then obtained. The expression of the incubation characteristic time criterion is as follows: ;in, is the dimensionless dynamic brittle fracture criterion function, t It's time, is the characteristic time of gestation, is the quasi-static strength, For The corresponding equivalent stress history.

6. The device according to claim 5, characterized in that The brittle dielectric rod is provided with a speckle region, and the experimental data includes deformation image data of the speckle region; the material mechanical property determination module is further configured to: determining the axial displacement field distribution and the axial velocity field distribution of the particle points on the surface of the brittle dielectric rod based on the deformation image of the speckle region; Determine, based on the axial displacement field distribution and the stress wave mass conservation equation in the Lagrangian coordinate system, a time-history curve of the axial strain of the brittle dielectric rod caused by the sequential strong dynamic load generated during the shock-explosion experiment; Determine, based on the axial displacement field distribution, the stress wave momentum conservation equation in the Lagrangian coordinate system, and the path line method integral, a stress time history curve corresponding to the deformation of the brittle medium caused by the time-sequential strong dynamic load generated during the shock-explosion experiment; Based on the strain time history curve and the stress time history curve, a stress-strain curve of the brittle dielectric material under a time-series strong dynamic load is obtained.

7. The device according to claim 6, characterized in that The material mechanical property determination module is further configured to include: The DIC technique is used to perform image analysis processing on the deformation image data of the speckle area to determine the axial displacement field distribution and axial velocity field distribution of the particle point on the surface of the brittle dielectric rod under time-series strong dynamic load.

8. An electronic device, characterized in that: The method comprises a processor and a memory, wherein the memory stores computer-executable instructions that can be executed by the processor, and the processor executes the computer-executable instructions to implement the method according to any one of claims 1 to 4.

9. A computer-readable storage medium, characterized in that The computer-readable storage medium stores computer-executable instructions. When the computer-executable instructions are called and executed by a processor, the computer-executable instructions prompt the processor to implement the method according to any one of claims 1 to 4.

Images