CN117875096A - Rock stress wave load propagation characteristic simulation method and related device - Google Patents
Rock stress wave load propagation characteristic simulation method and related device Download PDFInfo
- Publication number
- CN117875096A CN117875096A CN202410285480.4A CN202410285480A CN117875096A CN 117875096 A CN117875096 A CN 117875096A CN 202410285480 A CN202410285480 A CN 202410285480A CN 117875096 A CN117875096 A CN 117875096A
- Authority
- CN
- China
- Prior art keywords
- rock
- stress
- structural
- stress wave
- structural surface
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 239000011435 rock Substances 0.000 title claims abstract description 239
- 238000004088 simulation Methods 0.000 title claims abstract description 148
- 238000000034 method Methods 0.000 title claims abstract description 49
- 238000006073 displacement reaction Methods 0.000 claims abstract description 64
- 238000012544 monitoring process Methods 0.000 claims abstract description 17
- 238000004364 calculation method Methods 0.000 claims abstract description 14
- 230000006835 compression Effects 0.000 claims description 32
- 238000007906 compression Methods 0.000 claims description 32
- 238000012360 testing method Methods 0.000 claims description 24
- 238000004590 computer program Methods 0.000 claims description 16
- 230000014509 gene expression Effects 0.000 claims description 14
- 239000002245 particle Substances 0.000 claims description 14
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 12
- 239000011707 mineral Substances 0.000 claims description 12
- 238000010008 shearing Methods 0.000 claims description 10
- 238000010276 construction Methods 0.000 claims description 6
- 230000009977 dual effect Effects 0.000 claims description 6
- 238000012545 processing Methods 0.000 abstract description 11
- 230000004044 response Effects 0.000 abstract description 5
- 230000009471 action Effects 0.000 abstract description 4
- 230000005540 biological transmission Effects 0.000 abstract description 4
- 239000004579 marble Substances 0.000 description 12
- 239000000463 material Substances 0.000 description 12
- 238000010586 diagram Methods 0.000 description 10
- 230000008569 process Effects 0.000 description 8
- 238000012669 compression test Methods 0.000 description 7
- 238000009864 tensile test Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 6
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 5
- 229910052782 aluminium Inorganic materials 0.000 description 5
- 238000011156 evaluation Methods 0.000 description 3
- 230000006399 behavior Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000009533 lab test Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The invention belongs to the field of electric digital data processing, and discloses a rock stress wave load propagation characteristic simulation method and a related device, wherein the method comprises the following steps: constructing a rock numerical physical simulation model, and loading a preset stress wave into the rock numerical physical simulation model; and iteratively calling the step and the rigidity type determining step time by time until the stress wave loading is completed, and acquiring stress wave propagation characteristic index values of preset monitoring points at each time to obtain a rock stress wave load propagation characteristic simulation result. The method comprises the steps of judging the type of normal stiffness of the structural plane at the next moment through the type of normal displacement of the structural plane at each moment, referencing a rock numerical physical simulation model through a preset rock double-modulus constitutive model, realizing different calculation modes of stress of the structural plane under the condition of different types of normal stiffness of the structural plane, considering the influence of the double-modulus characteristic of the rock on stress wave load transmission, and more accurately simulating the dynamic response characteristic of the rock under the repeated action of tensile load.
Description
Technical Field
The invention belongs to the field of electric digital data processing, and relates to a rock stress wave load propagation characteristic simulation method and a related device.
Background
The tensile modulus of elasticity and the compressive modulus of elasticity of rock materials are the most intuitive indicators describing the ability of a rock to deform under external loads. Therefore, it is necessary to comprehensively consider the tensile elastic modulus and the compressive elastic modulus when simulating the stress wave load propagation characteristics of the rock. When the propagation of real stress wave loads such as seismic wave loads or blast wave loads in rock is simulated, the seismic wave loads and the blast wave loads are cyclic reciprocating loads, and compression and stretching and shearing actions exist, so that different deformation characteristics of rock materials under compression and stretching are important in propagation characteristic simulation.
However, currently, when the rock stress wave load propagation characteristic simulation is carried out, the adopted models for representing the dynamic deformation characteristic of the rock structural surface are all single modulus models, namely, the tensile elastic modulus and the compressive elastic modulus of the rock material are plotted with equal signs. However, it was found in the laboratory test that the tensile and compressive elastic moduli of the rock material were not the same, and there was a difference between them, and the reason for the difference was due to the microstructure of the rock.
Therefore, the dual-modulus characteristic of the rock structural surface cannot be accurately depicted by regarding the tensile elastic modulus and the compressive elastic modulus as single modulus, which leads to inaccuracy of rock dynamic response numerical calculation in rock stress wave load propagation characteristic simulation, and cannot accurately simulate the propagation characteristic of stress wave load in rock, so that errors or even errors can occur in rock engineering stability evaluation, and engineering construction safety is affected.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provide a rock stress wave load propagation characteristic simulation method and a related device.
In order to achieve the purpose, the invention is realized by adopting the following technical scheme:
in a first aspect of the present invention, there is provided a method of simulating the propagation characteristics of a rock stress wave load, comprising:
constructing a rock numerical physical simulation model, and loading a preset stress wave into the rock numerical physical simulation model;
iteratively calling the step and the rigidity type determining step time by time until the stress wave loading is completed, and acquiring stress wave propagation characteristic index values of preset monitoring points at each time to obtain a rock stress wave load propagation characteristic simulation result;
the calling step comprises the following steps: invoking a preset rock double-modulus constitutive model to carry out parameter assignment on the rock numerical physical simulation model, and carrying out simulation calculation based on the rock numerical physical simulation model after parameter assignment to obtain normal displacement of each structural surface in the rock numerical physical simulation model at the current moment and stress wave propagation characteristic index values of each preset monitoring point; the stiffness type determining step includes: traversing the normal displacement of the structural surface of each structural surface at the current moment, and setting the normal rigidity of the structural surface at the next moment of the current structural surface as the structural surface compression rigidity when the normal displacement of the structural surface is the structural surface compression displacement; and when the normal displacement of the current structural plane is structural plane tensile displacement, setting the structural plane normal rigidity at the next moment of the current structural plane as structural plane tensile rigidity.
Optionally, the rock double modulus constitutive model is obtained by the following method:
modifying a description expression between structural surface stress and structural surface closing quantity in the elastoplastic constitutive model into expressions (1) - (3) to obtain:
（1）
（2）
（3）
wherein,for structural plane compressive stress +.>For structural plane compression stiffness +.>For the structural plane compression displacement +.>For structural plane tensile stress +.>For structural plane tensile stiffness +.>For the structural surface stretching displacement->For structural surface tensile strength->For structural surface shear stress +.>For structural plane shear stiffness +.>For shearing displacement of structural surface->For structural plane compressive stress +.>For the structural surface cohesion ++>Is the friction angle in the structural plane.
Optionally, the constructing the petrophysical simulation model includes:
acquiring the external dimension of the rock to be simulated and the average diameter of mineral particles constituting the rock;
according to the external dimension of the rock to be simulated and the average diameter of mineral particles constituting the rock, a rock numerical physical simulation model is constructed based on a general discrete unit method program.
Optionally, when the preset rock double-modulus constitutive model is called for the rock numerical physical simulation model, the rock double-modulus constitutive model is called by calling a dynamic link library file of the rock double-modulus constitutive model.
Optionally, when the preset rock double-modulus constitutive model is called to add parameters to the rock numerical physical simulation model, parameter values of all parameters to be assigned of the rock numerical physical simulation model are obtained through a test calibration step;
the test calibration steps comprise:
obtaining a rock stress strain curve through a rock stress application test;
under the condition of given initial values, iteratively updating parameter values of all the parameters to be assigned of the rock numerical physical simulation model, obtaining a rock simulation stress-strain curve simulated based on the rock numerical physical simulation model after each iteration until the similarity between the rock simulation stress-strain curve and the rock stress-strain curve is not smaller than a preset threshold, and taking the current parameter values of all the parameters to be assigned of the rock numerical physical simulation model as final parameter values of all the parameters to be assigned of the rock numerical physical simulation model.
Optionally, the stress wave propagation characteristic index value includes one or several of the following: speed, displacement, and stress.
In a second aspect of the invention, there is provided a rock stress wave load propagation characteristic simulation system comprising:
the model construction module is used for constructing a rock numerical physical simulation model and loading a preset stress wave to the rock numerical physical simulation model;
the propagation simulation module is used for iteratively calling the step and the rigidity type determining step from moment to moment until the stress wave loading is completed, and obtaining stress wave propagation characteristic index values of preset monitoring points at each moment to obtain a rock stress wave load propagation characteristic simulation result;
the calling step comprises the following steps: invoking a preset rock double-modulus constitutive model to carry out parameter assignment on the rock numerical physical simulation model, and carrying out simulation calculation based on the rock numerical physical simulation model after parameter assignment to obtain normal displacement of each structural surface in the rock numerical physical simulation model at the current moment and stress wave propagation characteristic index values of each preset monitoring point; the stiffness type determining step includes: traversing the normal displacement of the structural surface of each structural surface at the current moment, and setting the normal rigidity of the structural surface at the next moment of the current structural surface as the structural surface compression rigidity when the normal displacement of the structural surface is the structural surface compression displacement; and when the normal displacement of the current structural plane is structural plane tensile displacement, setting the structural plane normal rigidity at the next moment of the current structural plane as structural plane tensile rigidity.
Optionally, the rock double modulus constitutive model is obtained by the following method:
modifying a description expression between structural surface stress and structural surface closing quantity in the elastoplastic constitutive model into expressions (1) - (3) to obtain:
（1）
（2）
（3）
wherein,for structural plane compressive stress +.>For structural plane compression stiffness +.>Is a structural surfaceCompression displacement(s)>For structural plane tensile stress +.>For structural plane tensile stiffness +.>For the structural surface stretching displacement->For structural surface tensile strength->For structural surface shear stress +.>For structural plane shear stiffness +.>For shearing displacement of structural surface->For structural plane compressive stress +.>For the structural surface cohesion ++>Is the friction angle in the structural plane.
In a third aspect the present invention provides a computer device comprising a memory, a processor and a computer program stored in said memory and operable on said processor, said processor implementing the steps of the rock stress wave load propagation characteristic simulation method described above when said computer program is executed.
In a fourth aspect of the present invention, there is provided a computer readable storage medium storing a computer program which when executed by a processor performs the steps of the rock stress wave load propagation characteristic simulation method described above.
Compared with the prior art, the invention has the following beneficial effects:
according to the rock stress wave load propagation characteristic simulation method, a rock numerical physical simulation model is constructed to serve as a simulation basis of stress wave load propagation, then preset stress waves are loaded into the rock numerical physical simulation model, the stress wave loading is completed through iteration calling steps and stiffness type determining steps from time to time, the whole stress wave load propagation process is simulated, and finally stress wave propagation characteristic index values of preset monitoring points at each time are collected to obtain a rock stress wave load propagation characteristic simulation result. The method comprises the steps of judging the type of normal stiffness of the structural plane at the next moment by the type of normal displacement of the structural plane at each moment, and further, referencing a rock numerical physical simulation model by calling a preset rock double-modulus constitutive model to realize different calculation modes of stress of the structural plane under the condition of different types of normal stiffness of the structural plane, fully considering the influence of the double-modulus characteristic of the rock on stress wave load transmission, more accurately simulating the dynamic response characteristic of the rock under the repeated action of the tensile load, improving the accuracy of engineering stability evaluation results and guaranteeing engineering construction safety.
Drawings
FIG. 1 is a flow chart of a method for simulating the propagation characteristics of rock stress wave load according to an embodiment of the invention.
Fig. 2 is a schematic diagram of a uniaxial compression and direct stretching particle model of marble based on a general discrete unit process program according to an embodiment of the present invention.
Fig. 3 is a stress-strain curve obtained by uniaxial compression test and simulation of a rock and a stress-strain curve obtained by direct tensile test and simulation of a rock according to an embodiment of the present invention.
Fig. 4 is an ultrasonic waveform through a 40mm thick aluminum block in accordance with an embodiment of the present invention.
Fig. 5 is an ultrasonic waveform through 40mm thick marble in accordance with an embodiment of the present invention.
Fig. 6 is a waveform diagram of an incident wave of the present invention after half cosine cone processing of an ultrasonic wave passing through an aluminum block according to an embodiment of the present invention.
Fig. 7 is a waveform diagram of a transmitted wave of the ultrasonic wave passing through marble according to an embodiment of the present invention after half cosine cone processing.
Fig. 8 is a waveform of transmitted waves through 40mm marble using a rock dual modulus constitutive model and a conventional single modulus constitutive model simulated ultrasonic waves under the same parameters in accordance with an embodiment of the invention.
Fig. 9 is a block diagram of a rock stress wave load propagation characteristic simulation system according to an embodiment of the present invention.
Detailed Description
In order that those skilled in the art will better understand the present invention, a technical solution in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in which it is apparent that the described embodiments are only some embodiments of the present invention, not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the present invention without making any inventive effort, shall fall within the scope of the present invention.
It should be noted that the terms "first," "second," and the like in the description and the claims of the present invention and the above figures are used for distinguishing between similar objects and not necessarily for describing a particular sequential or chronological order. It is to be understood that the data so used may be interchanged where appropriate such that the embodiments of the invention described herein may be implemented in sequences other than those illustrated or otherwise described herein. Furthermore, the terms "comprises," "comprising," and "having," and any variations thereof, are intended to cover a non-exclusive inclusion, such that a process, method, system, article, or apparatus that comprises a list of steps or elements is not necessarily limited to those steps or elements expressly listed but may include other steps or elements not expressly listed or inherent to such process, method, article, or apparatus.
The invention is described in further detail below with reference to the attached drawing figures:
referring to fig. 1, in an embodiment of the present invention, a method for simulating propagation characteristics of a stress wave load of a rock is provided, so that a description of a difference between a tensile elastic modulus and a compressive elastic modulus of the rock is realized, propagation characteristics of the stress wave load when dual modulus characteristics of the rock material are considered can be accurately simulated, and an intrinsic mechanism of dual modulus characteristics of the rock material is revealed. Specifically, the rock stress wave load propagation characteristic simulation method comprises the following steps:
s1: and constructing a rock numerical physical simulation model, and loading the preset stress wave to the rock numerical physical simulation model.
S2: and iteratively calling the step and the rigidity type determining step time by time until the stress wave loading is completed, and acquiring stress wave propagation characteristic index values of preset monitoring points at each time to obtain a rock stress wave load propagation characteristic simulation result.
The calling step comprises the following steps: invoking a preset rock double-modulus constitutive model to carry out parameter assignment on the rock numerical physical simulation model, and carrying out simulation calculation based on the rock numerical physical simulation model after parameter assignment to obtain normal displacement of each structural surface in the rock numerical physical simulation model at the current moment and stress wave propagation characteristic index values of each preset monitoring point; the stiffness type determining step includes: traversing the normal displacement of the structural surface of each structural surface at the current moment, and setting the normal rigidity of the structural surface at the next moment of the current structural surface as the structural surface compression rigidity when the normal displacement of the structural surface is the structural surface compression displacement; and when the normal displacement of the current structural plane is structural plane tensile displacement, setting the structural plane normal rigidity at the next moment of the current structural plane as structural plane tensile rigidity.
According to the rock stress wave load propagation characteristic simulation method, a rock numerical physical simulation model is constructed to serve as a simulation basis of stress wave load propagation, then preset stress waves are loaded into the rock numerical physical simulation model, the stress wave loading is completed through iteration calling steps and stiffness type determining steps from time to time, the whole stress wave load propagation process is simulated, and finally stress wave propagation characteristic index values of preset monitoring points at each time are collected to obtain a rock stress wave load propagation characteristic simulation result. The method comprises the steps of judging the type of normal stiffness of the structural plane at the next moment by the type of normal displacement of the structural plane at each moment, and further referencing a rock numerical physical simulation model by calling a preset rock double-modulus constitutive model to realize different calculation modes of stress of the structural plane under the condition of different types of normal stiffness of the structural plane, fully considering the influence of the double-modulus characteristic of the rock on stress wave load transmission, more accurately simulating the dynamic response characteristic of the rock under the repeated action of the tensile load, improving the accuracy of an engineering stability evaluation result and guaranteeing the safety of engineering construction.
In one possible embodiment, the rock double modulus constitutive model is obtained by:
modifying a description expression between structural surface stress and structural surface closing quantity in the elastoplastic constitutive model into expressions (1) - (3) to obtain:
（1）
（2）
（3）
wherein,for structural plane compressive stress +.>For structural plane compression stiffness +.>For the structural plane compression displacement +.>For structural plane tensile stress +.>Is of a structureFace tensile stiffness->For the structural surface stretching displacement->For structural surface tensile strength->For structural surface shear stress +.>For structural plane shear stiffness +.>For shearing displacement of structural surface->For structural plane compressive stress +.>For the structural surface cohesion ++>Is the friction angle in the structural plane.
In particular, an elastoplastic constitutive model is a mathematical model describing the behavior of a material under stress. The model is based on elastoplastics and is mainly used for analyzing the relationship between deformation and stress of materials when the materials are subjected to pressure. In an elastoplastic constitutive model, a material is considered to be a continuous medium whose behavior can be described by a set of relationships that relate parameters describing the deformation of the continuous medium to parameters describing the internal forces. Specifically, the elastoplastic constitutive model refers to a set of relational expressions that relate the strain tensor of the deformation to the strain tensor.
In this embodiment, the implementation of the description expression between the structural plane stress and the structural plane closing amount in the elastoplastic constitutive model source code can be modified based on Visual C++6.0 software.
Specifically, parameters are given to the rock numerical physical simulation model through the elastoplastic constitutive model, namely, specific parameter values are given to each parameter to be assigned of the rock numerical physical simulation model through the elastoplastic constitutive model, and corresponding calculation rules are given to each parameter of the rock numerical physical simulation model.
For the rock numerical physical simulation model, the normal displacement of the structural surface refers to structural surface tensile displacement or structural surface compressive displacement, and the normal stiffness of the structural surface refers to structural surface tensile stiffness or structural surface compressive stiffness. The normal displacement of the structural surface is regarded as the structural surface tensile displacement when the normal displacement of the structural surface is smaller than zero, and the normal displacement of the structural surface is regarded as the structural surface compressive displacement when the normal displacement of the structural surface is larger than zero. Therefore, when the rock double-modulus constitutive model is called to give parameters to the rock numerical physical simulation model, different stress calculation formulas are given for different structural plane normal rigidities, and if the structural plane normal rigidity is structural plane compression rigidity, the structural plane compression stress is calculated by adopting formula (1).
In one possible embodiment, the constructing the petrophysical simulation model comprises: acquiring the external dimension of the rock to be simulated and the average diameter of mineral particles constituting the rock; according to the external dimension of the rock to be simulated and the average diameter of mineral particles constituting the rock, a rock numerical physical simulation model is constructed based on a general discrete unit method program.
Specifically, the universal discrete unit method program (Universal Distinct Element Code, UDEC) is a calculation analysis program based on the theory of discrete unit methods. The method is a tool for providing accurate and effective analysis for geotechnical engineering by using an explicit problem solving scheme, the explicit problem solving scheme provides stable solution for unstable physical processes and can simulate the damage process of objects, and the software is particularly suitable for simulating the response of a joint rock system or a discontinuous block aggregate system under the condition of static or dynamic load.
Based on the method, simulation calculation of the rock numerical physical simulation model can be realized through a discrete unit method program, and the normal displacement of the structural surface of each structural surface in the rock numerical physical simulation model and the stress wave propagation characteristic index value of each preset monitoring point are obtained. Each preset monitoring point can be a structural surface or a mineral particle forming rock. The index value of the propagation characteristic of the stress wave may be one or more of a velocity, a displacement, and a stress.
In one possible implementation manner, when the preset rock double-modulus constitutive model is called to reference the rock numerical physical simulation model, the rock double-modulus constitutive model is called by calling a dynamic link library file of the rock double-modulus constitutive model.
Specifically, the dynamic link library file of the rock double-modulus constitutive model is generated in advance, and then when the rock double-modulus constitutive model is required to be used, the rock double-modulus constitutive model can be called in a mode of calling the dynamic link library file of the rock double-modulus constitutive model, so that the rock double-modulus constitutive model can be conveniently called.
In one possible implementation manner, when the preset rock double-modulus constitutive model is called to add parameters to the rock numerical physical simulation model, parameter values of all parameters to be assigned of the rock numerical physical simulation model are obtained through a test calibration step; the test calibration steps comprise: obtaining a rock stress strain curve through a rock stress application test; under the condition of given initial values, iteratively updating parameter values of all the parameters to be assigned of the rock numerical physical simulation model, obtaining a rock simulation stress-strain curve simulated based on the rock numerical physical simulation model after each iteration until the similarity between the rock simulation stress-strain curve and the rock stress-strain curve is not smaller than a preset threshold, and taking the current parameter values of all the parameters to be assigned of the rock numerical physical simulation model as final parameter values of all the parameters to be assigned of the rock numerical physical simulation model.
Specifically, each parameter to be assigned of the rock numerical physical simulation model is generally considered as a micromechanics parameter of the rock material, which is difficult or impossible to directly measure through an indoor test, but the macroscopic mechanics parameter of the rock, such as compressive strength, tensile strength, shearing strength and the like, can be obtained through conventional test means of rock stress application tests, such as rock uniaxial compression, direct stretching, direct shearing tests and the like. Then, in numerical simulation, a certain initial value of the parameters is given according to experience, tests such as uniaxial compression, direct stretching or direct shearing of the rock are simulated, a simulated stress-strain curve under the corresponding test of the rock is obtained, and the simulated stress-strain curve is compared with the stress-strain curve obtained by the test. If the results of the simulated stress-strain curve and the stress-strain curve are inconsistent, namely the similarity between the simulated stress-strain curve and the stress-strain curve is not high, the parameters are modified, the steps are repeated again until the simulated stress-strain curve is matched with the stress-strain curve obtained through the test to the greatest extent, namely the similarity between the simulated stress-strain curve and the stress-strain curve is not smaller than a preset threshold, the parameters obtained at the moment are considered to be consistent with the microscopic parameters of the actual rock, the microscopic mechanical characteristics of the actual rock can be represented, and the calibration of the parameter values of the parameters to be assigned of the rock numerical physical simulation model is completed.
In one possible embodiment, the accuracy of the rock bi-modulus constitutive model is verified by calibration of the rock uniaxial compression test and the rock direct tensile test results.
In this embodiment, the marble sample used in the uniaxial compression test and the direct tensile test of rock was a cylindrical standard sample having a diameter of 50mm and a height of 100mm, and the average diameter of mineral particles constituting the marble was 1.5mm. Before numerical simulation, firstly, a rock numerical physical simulation model identical to a test sample is established, and since UDEC is 2-dimensional software, a rectangular rock numerical physical simulation model with the width of 50mm and the height of 100mm as shown in fig. 2 is established according to the size of the test sample, the rock numerical physical simulation model established in the numerical simulation software is subjected to particle division, the average diameter of mineral particles forming marble is set to be 1.5mm, the mineral particles are consistent with the test sample, and the boundary between the particles is a structural surface. After the rock numerical physical simulation model is established, the double-modulus constitutive model is called to endow a structural surface with parameters, and a stress-strain curve of direct stretching of the rock mass and a stress-strain curve of uniaxial compression of the rock mass are simulated.
Referring to fig. 3, stress strain curves obtained by uniaxial compression test and simulation of the rock and stress strain curves obtained by direct tensile test and simulation of the rock are given. From the experimental results, it can be seen that the elastic modulus obtained by the uniaxial compression test of the rock is different from that obtained by the direct tensile test of the rock, the compressive elastic modulus Ec of the uniaxial compression test of the rock is 21.94GPa, and the tensile elastic modulus Et of the direct tensile test of the rock is 11.08GPa, and the two elastic moduliare calculated by the secant of the stress-strain curve obtained by the uniaxial compression test and the direct tensile test. Meanwhile, from fig. 3, it can be clearly seen that two elastic moduli obtained by simulation based on the rock double modulus constitutive model established in the invention are highly matched with the calculated value of the test result, which indicates that the rock double modulus constitutive model has higher accuracy.
In one possible embodiment, an ultrasonic test was performed in order to verify the effect of the rock dual modulus property on the stress wave load propagation properties.
Referring to fig. 4 and 5, ultrasonic waveforms of ultrasonic waves passing through an aluminum block of 40mm thickness and marble of 40mm thickness, respectively, were employed with a center frequency of 67Hz. In this embodiment, an ultrasonic waveform during a period from the start of the jump point is used as the initial pulse. Referring to fig. 6 and 7, the incident wave waveform and the transmitted wave waveform used in the present embodiment after the half cosine cone processing of fig. 4 and 5, respectively. The aluminum block can be considered as an isotropic homogeneous body and the ultrasonic wave propagating through the aluminum block can be considered as a numerically simulated input stress wave, while the amplitude of the processed pulse in fig. 6 is approximately equal in the compression (above the horizontal axis) and tension (below the horizontal axis) directions. While marble has a dual-modulus characteristic, it can be observed in fig. 7 that the amplitude of the processed pulse in the tensile direction is significantly smaller than that in the compressive direction, so that the ultrasonic wave passing through marble can be calibrated as a test result of numerical simulation. The building process of the rock numerical physical simulation model is basically consistent with that in the embodiment, the marble rock block adopted in the ultrasonic test is a cylinder with the height of 40mm and the diameter of 50mm, the rock numerical physical simulation model is divided into mineral grain aggregates with the average diameter of 1.5mm of mineral particles forming the rock, the double-modulus constitutive model is called to endow a structural surface with parameters, and the simulation of ultrasonic wave propagation in the rock is carried out based on UDEC.
In the simulation of the propagation of stress waves in UDEC, a viscous boundary condition needs to be set to absorb stress waves transmitted to the boundary, and reflection of stress waves at the boundary is avoided.
Referring to fig. 8, there is shown the transmission wave waveform of ultrasonic waves through 40mm marble simulated using the rock dual modulus constitutive model and the conventional single modulus constitutive model under the same parameters. As can be seen from fig. 8, when the conventional single modulus constitutive model simulates ultrasonic wave propagation, it is difficult to describe the phenomenon that the tensile amplitude value is attenuated due to the fact that the tensile modulus is smaller than the compressive modulus, and the rock double modulus constitutive model established in the invention can well describe the characteristic. Therefore, compared with the traditional single modulus constitutive model, the rock double modulus constitutive model can more accurately characterize the propagation characteristics of stress waves in rock materials.
The following are device embodiments of the present invention that may be used to perform method embodiments of the present invention. For details not disclosed in the apparatus embodiments, please refer to the method embodiments of the present invention.
Referring to fig. 9, in still another embodiment of the present invention, a rock stress wave load propagation characteristic simulation system is provided, which can be used to implement the rock stress wave load propagation characteristic simulation method described above, and in particular, the rock stress wave load propagation characteristic simulation system includes a model building module and a propagation simulation module.
The model construction module is used for constructing a rock numerical physical simulation model and loading a preset stress wave to the rock numerical physical simulation model; the propagation simulation module is used for iteratively calling the step and the rigidity type determining step time by time until the stress wave loading is completed, and obtaining stress wave propagation characteristic index values of preset monitoring points at each time to obtain a rock stress wave load propagation characteristic simulation result; the calling step comprises the following steps: invoking a preset rock double-modulus constitutive model to carry out parameter assignment on the rock numerical physical simulation model, and carrying out simulation calculation based on the rock numerical physical simulation model after parameter assignment to obtain normal displacement of each structural surface in the rock numerical physical simulation model at the current moment and stress wave propagation characteristic index values of each preset monitoring point; the stiffness type determining step includes: traversing the normal displacement of the structural surface of each structural surface at the current moment, and setting the normal rigidity of the structural surface at the next moment of the current structural surface as the structural surface compression rigidity when the normal displacement of the structural surface is the structural surface compression displacement; and when the normal displacement of the current structural plane is structural plane tensile displacement, setting the structural plane normal rigidity at the next moment of the current structural plane as structural plane tensile rigidity.
In one possible embodiment, the rock double modulus constitutive model is obtained by:
modifying a description expression between structural surface stress and structural surface closing quantity in the elastoplastic constitutive model into expressions (1) - (3) to obtain:
（1）
（2）
（3）
wherein,for structural plane compressive stress +.>For structural plane compression stiffness +.>For the structural plane compression displacement +.>For structural plane tensile stress +.>For structural plane tensile stiffness +.>For the structural surface stretching displacement->For structural surface tensile strength->For structural surface shear stress +.>For structural plane shear stiffness +.>For shearing displacement of structural surface->For structural plane compressive stress +.>For the structural surface cohesion ++>Is the friction angle in the structural plane.
All relevant contents of the steps involved in the foregoing embodiments of the rock stress wave load propagation characteristic simulation method may be cited in the functional description of the functional module corresponding to the rock stress wave load propagation characteristic simulation system in the embodiment of the present invention, which is not described herein.
The division of the modules in the embodiments of the present invention is schematically only one logic function division, and there may be another division manner in actual implementation, and in addition, each functional module in each embodiment of the present invention may be integrated in one processor, or may exist separately and physically, or two or more modules may be integrated in one module. The integrated modules may be implemented in hardware or in software functional modules.
In yet another embodiment of the present invention, a computer device is provided that includes a processor and a memory for storing a computer program including program instructions, the processor for executing the program instructions stored by the computer storage medium. The processor may be a central processing unit (Central Processing Unit, CPU), but may also be other general purpose processor, digital signal processor (Digital Signal Processor, DSP), application specific integrated circuit (Application Specific Integrated Circuit, ASIC), off-the-shelf Programmable gate array (FPGA) or other Programmable logic device, discrete gate or transistor logic device, discrete hardware components, etc., which are the computational core and control core of the terminal adapted to implement one or more instructions, in particular to load and execute one or more instructions in a computer storage medium to implement the corresponding method flow or corresponding functions; the processor provided by the embodiment of the invention can be used for the operation of a rock stress wave load propagation characteristic simulation method.
In yet another embodiment of the present invention, a storage medium, specifically a computer readable storage medium (Memory), is a Memory device in a computer device, for storing a program and data. It is understood that the computer readable storage medium herein may include both built-in storage media in a computer device and extended storage media supported by the computer device. The computer-readable storage medium provides a storage space storing an operating system of the terminal. Also stored in the memory space are one or more instructions, which may be one or more computer programs (including program code), adapted to be loaded and executed by the processor. The computer readable storage medium herein may be a high-speed RAM memory or a non-volatile memory (non-volatile memory), such as at least one magnetic disk memory. One or more instructions stored in a computer-readable storage medium may be loaded and executed by a processor to implement the respective steps of the method for modeling rock stress wave load propagation characteristics in connection with the above-described embodiments.
It will be appreciated by those skilled in the art that embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the invention. It will be understood that each flow or block of the flowchart illustrations or block diagrams, and combinations of flows or blocks in the flowchart illustrations or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows or block diagram block or blocks.
Finally, it should be noted that: the above embodiments are only for illustrating the technical aspects of the present invention and not for limiting the same, and although the present invention has been described in detail with reference to the above embodiments, it should be understood by those of ordinary skill in the art that: modifications and equivalents may be made to the specific embodiments of the invention without departing from the spirit and scope of the invention, which is intended to be covered by the claims.
Claims (10)
1. A method of modeling rock stress wave load propagation characteristics, comprising:
constructing a rock numerical physical simulation model, and loading a preset stress wave into the rock numerical physical simulation model;
iteratively calling the step and the rigidity type determining step time by time until the stress wave loading is completed, and acquiring stress wave propagation characteristic index values of preset monitoring points at each time to obtain a rock stress wave load propagation characteristic simulation result;
the calling step comprises the following steps: invoking a preset rock double-modulus constitutive model to carry out parameter assignment on the rock numerical physical simulation model, and carrying out simulation calculation based on the rock numerical physical simulation model after parameter assignment to obtain normal displacement of each structural surface in the rock numerical physical simulation model at the current moment and stress wave propagation characteristic index values of each preset monitoring point; the stiffness type determining step includes: traversing the normal displacement of the structural surface of each structural surface at the current moment, and setting the normal rigidity of the structural surface at the next moment of the current structural surface as the structural surface compression rigidity when the normal displacement of the structural surface is the structural surface compression displacement; and when the normal displacement of the current structural plane is structural plane tensile displacement, setting the structural plane normal rigidity at the next moment of the current structural plane as structural plane tensile rigidity.
2. The rock stress wave load propagation characteristic simulation method according to claim 1, wherein the rock double modulus constitutive model is obtained by:
modifying a description expression between structural surface stress and structural surface closing quantity in the elastoplastic constitutive model into expressions (1) - (3) to obtain:
（1）
（2）
（3）
wherein,for structural plane compressive stress +.>For structural plane compression stiffness +.>For the structural plane compression displacement +.>For structural plane tensile stress +.>For structural plane tensile stiffness +.>For the structural surface stretching displacement->For structural surface tensile strength->For structural surface shear stress +.>Is of a structureFace shear stiffness->For shearing displacement of structural surface->For structural plane compressive stress +.>For the structural surface cohesion ++>Is the friction angle in the structural plane.
3. The rock stress wave load propagation characteristics simulation method according to claim 1, wherein the constructing a rock numerical physical simulation model comprises:
acquiring the external dimension of the rock to be simulated and the average diameter of mineral particles constituting the rock;
according to the external dimension of the rock to be simulated and the average diameter of mineral particles constituting the rock, a rock numerical physical simulation model is constructed based on a general discrete unit method program.
4. A method for simulating the propagation characteristics of a rock stress wave load according to claim 3, wherein when the preset rock double-modulus constitutive model is called for parametrizing the rock numerical physical simulation model, the rock double-modulus constitutive model is called by calling a dynamic link library file of the rock double-modulus constitutive model.
5. The rock stress wave load propagation characteristic simulation method according to claim 1, wherein when a preset rock double-modulus constitutive model is called to parameterize a rock numerical physical simulation model, parameter values of all parameters to be evaluated of the rock numerical physical simulation model are obtained through a test calibration step;
the test calibration steps comprise:
obtaining a rock stress strain curve through a rock stress application test;
under the condition of given initial values, iteratively updating parameter values of all the parameters to be assigned of the rock numerical physical simulation model, obtaining a rock simulation stress-strain curve simulated based on the rock numerical physical simulation model after each iteration until the similarity between the rock simulation stress-strain curve and the rock stress-strain curve is not smaller than a preset threshold, and taking the current parameter values of all the parameters to be assigned of the rock numerical physical simulation model as final parameter values of all the parameters to be assigned of the rock numerical physical simulation model.
6. The rock stress wave load propagation characteristic simulation method according to claim 1, wherein the stress wave propagation characteristic index value comprises one or several of the following: speed, displacement, and stress.
7. A rock stress wave load propagation characteristic simulation system, comprising:
the model construction module is used for constructing a rock numerical physical simulation model and loading a preset stress wave to the rock numerical physical simulation model;
the propagation simulation module is used for iteratively calling the step and the rigidity type determining step from moment to moment until the stress wave loading is completed, and obtaining stress wave propagation characteristic index values of preset monitoring points at each moment to obtain a rock stress wave load propagation characteristic simulation result;
the calling step comprises the following steps: invoking a preset rock double-modulus constitutive model to carry out parameter assignment on the rock numerical physical simulation model, and carrying out simulation calculation based on the rock numerical physical simulation model after parameter assignment to obtain normal displacement of each structural surface in the rock numerical physical simulation model at the current moment and stress wave propagation characteristic index values of each preset monitoring point; the stiffness type determining step includes: traversing the normal displacement of the structural surface of each structural surface at the current moment, and setting the normal rigidity of the structural surface at the next moment of the current structural surface as the structural surface compression rigidity when the normal displacement of the structural surface is the structural surface compression displacement; and when the normal displacement of the current structural plane is structural plane tensile displacement, setting the structural plane normal rigidity at the next moment of the current structural plane as structural plane tensile rigidity.
8. The rock stress wave load propagation characteristics simulation system according to claim 7, wherein the rock dual modulus constitutive model is obtained by:
modifying a description expression between structural surface stress and structural surface closing quantity in the elastoplastic constitutive model into expressions (1) - (3) to obtain:
（1）
（2）
（3）
wherein,for structural plane compressive stress +.>For structural plane compression stiffness +.>For the structural plane compression displacement +.>For structural plane tensile stress +.>For structural plane tensile stiffness +.>For the structural surface stretching displacement->For structural surface tensile strength->For structural surface shear stress +.>For structural plane shear stiffness +.>For shearing displacement of structural surface->For structural plane compressive stress +.>For the structural surface cohesion ++>Is the friction angle in the structural plane.
9. A computer device comprising a memory, a processor and a computer program stored in the memory and executable on the processor, characterized in that the processor, when executing the computer program, carries out the steps of the rock stress wave load propagation characteristic simulation method according to any one of claims 1 to 6.
10. A computer readable storage medium storing a computer program, characterized in that the computer program when executed by a processor implements the steps of the rock stress wave load propagation characteristic simulation method according to any one of claims 1 to 6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410285480.4A CN117875096B (en) | 2024-03-13 | 2024-03-13 | Rock stress wave load propagation characteristic simulation method and related device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410285480.4A CN117875096B (en) | 2024-03-13 | 2024-03-13 | Rock stress wave load propagation characteristic simulation method and related device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117875096A true CN117875096A (en) | 2024-04-12 |
| CN117875096B CN117875096B (en) | 2024-05-24 |
Family
ID=90590509
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410285480.4A Active CN117875096B (en) | 2024-03-13 | 2024-03-13 | Rock stress wave load propagation characteristic simulation method and related device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117875096B (en) |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107014704A (en) * | 2017-05-15 | 2017-08-04 | 东北大学 | A kind of short rock bar viscosity coefficient method of testing that analysis is propagated based on viscoelasticity ripple |
| CN107463740A (en) * | 2017-07-27 | 2017-12-12 | 中南大学 | Consider the rock type materials true triaxial test method for numerical simulation of intermediate principal stress effect |
| US20210012046A1 (en) * | 2019-07-12 | 2021-01-14 | Beihang University | Meshless method for solid mechanics simulation, electronic device, and storage medium |
| CN115358111A (en) * | 2022-07-25 | 2022-11-18 | 中国电建集团华东勘测设计研究院有限公司 | DFM (distributed finite model) -based method, device and medium for simulating and analyzing stability of reservoir bank fractured rock mass slope in water storage period |
| CN115639057A (en) * | 2022-10-14 | 2023-01-24 | 中国矿业大学 | Inversion method of rock blasting equivalent loading stress curve |
| CN115859714A (en) * | 2022-11-23 | 2023-03-28 | 长安大学 | Rock blasting overall process simulation method based on FEM-DEM combined simulation |
| CN115964901A (en) * | 2023-03-16 | 2023-04-14 | 中国科学院地质与地球物理研究所 | Simulation method and system for water-induced rock strength deterioration based on discrete unit method |
| CN116772671A (en) * | 2023-06-27 | 2023-09-19 | 中交一公局集团有限公司 | Layered slope micro-vibration micro-damage accurate blasting construction method |
| CN117053645A (en) * | 2023-06-27 | 2023-11-14 | 中交一公局集团有限公司 | Layered slope precise blasting construction method |
| CN117634228A (en) * | 2024-01-26 | 2024-03-01 | 中国科学院地质与地球物理研究所 | Discontinuous deformation analysis simulation calculation method for reservoir controllable shock wave fracturing |
-
2024
- 2024-03-13 CN CN202410285480.4A patent/CN117875096B/en active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107014704A (en) * | 2017-05-15 | 2017-08-04 | 东北大学 | A kind of short rock bar viscosity coefficient method of testing that analysis is propagated based on viscoelasticity ripple |
| CN107463740A (en) * | 2017-07-27 | 2017-12-12 | 中南大学 | Consider the rock type materials true triaxial test method for numerical simulation of intermediate principal stress effect |
| US20210012046A1 (en) * | 2019-07-12 | 2021-01-14 | Beihang University | Meshless method for solid mechanics simulation, electronic device, and storage medium |
| CN115358111A (en) * | 2022-07-25 | 2022-11-18 | 中国电建集团华东勘测设计研究院有限公司 | DFM (distributed finite model) -based method, device and medium for simulating and analyzing stability of reservoir bank fractured rock mass slope in water storage period |
| CN115639057A (en) * | 2022-10-14 | 2023-01-24 | 中国矿业大学 | Inversion method of rock blasting equivalent loading stress curve |
| CN115859714A (en) * | 2022-11-23 | 2023-03-28 | 长安大学 | Rock blasting overall process simulation method based on FEM-DEM combined simulation |
| CN115964901A (en) * | 2023-03-16 | 2023-04-14 | 中国科学院地质与地球物理研究所 | Simulation method and system for water-induced rock strength deterioration based on discrete unit method |
| CN116772671A (en) * | 2023-06-27 | 2023-09-19 | 中交一公局集团有限公司 | Layered slope micro-vibration micro-damage accurate blasting construction method |
| CN117053645A (en) * | 2023-06-27 | 2023-11-14 | 中交一公局集团有限公司 | Layered slope precise blasting construction method |
| CN117634228A (en) * | 2024-01-26 | 2024-03-01 | 中国科学院地质与地球物理研究所 | Discontinuous deformation analysis simulation calculation method for reservoir controllable shock wave fracturing |
Non-Patent Citations (4)
| Title |
|---|
| 周创兵;陈益峰;姜清辉;荣冠;: "岩体结构面HM耦合分析的界面层模型", 岩石力学与工程学报, no. 06, 15 June 2008 (2008-06-15) * |
| 郭保华;李小军;苏承东;: "岩石裂隙法向循环加载本构关系试验研究", 岩石力学与工程学报, no. 1, 15 May 2012 (2012-05-15) * |
| 陈鹏宇;: "岩石颗粒流模型单轴压缩的加载速率效应研究", 地下空间与工程学报, no. 03, 15 June 2018 (2018-06-15) * |
| 黄志平;陈国平;唐春安;黎罡;: "应力波在岩石直杆中传播特性及层裂研究", 沈阳建筑大学学报(自然科学版), no. 04, 15 July 2020 (2020-07-15) * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117875096B (en) | 2024-05-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP2975377B1 (en) | Stress-strain relationship simulation method, spring back prediction method, and spring back analyzing device | |
| Lu et al. | A three-dimensional elastoplastic constitutive model for concrete | |
| Asadi et al. | A bonded particle model simulation of shear strength and asperity degradation for rough rock fractures | |
| Ma et al. | Microscopic modeling of the creep behavior of rockfills with a delayed particle breakage model | |
| CN103792143A (en) | Quick acquisition method of true stress strain curve in whole process of uniaxial drawing | |
| JP3708067B2 (en) | Method for predicting crack growth and deformation of elastoplastic material | |
| Mao et al. | The construction and comparison of damage detection index based on the nonlinear output frequency response function and experimental analysis | |
| JP3897477B2 (en) | Stress-strain relationship simulation method and springback amount prediction method | |
| Schreyer et al. | Analytical and numerical tests for loss of material stability | |
| Brekelmans et al. | Reduction of mesh sensitivity in continuum damage mechanics | |
| Hwang et al. | Evaluation of dynamic tensile strength of concrete using lattice-based simulations of spalling tests | |
| Tian et al. | A viscoelastic model of geometry-constraint-based non-ordinary state-based peridynamics with progressive damage | |
| CN117875096B (en) | Rock stress wave load propagation characteristic simulation method and related device | |
| CN109388833B (en) | Elastic element structure optimization design method based on fatigue life | |
| JP6852727B2 (en) | Inspection equipment, inspection method, and inspection program | |
| CN116542011A (en) | Rock water absorption softening effect particle discrete element simulation method, device and storage medium | |
| CN101866377A (en) | The method and system of simulation material aging effect in the computer-aided engineering analysis | |
| Saksala | On the strain rate sensitivity of coarse-grained rock: a mesoscopic numerical study | |
| JP2012037305A (en) | Sequential nonlinear earthquake response analysis method for foundation and storage medium with analysis program stored thereon | |
| Khademian et al. | Structural damage detection of reinforced concrete shear walls subject to consequent earthquakes | |
| CN110333148B (en) | Soil dynamic shear modulus testing method based on fine analysis of vibration attenuation curve | |
| Haq et al. | Parametric analysis of RSB sensors for concrete strength monitoring using hybrid EMI and WP techniques: Numerical investigation | |
| Gatta et al. | Dynamic characterization of a system with degradation: A masonry wall | |
| CN114414409A (en) | Method and device for determining fatigue performance of material | |
| JP2004045294A (en) | Determination system and program for risk of damaging structure |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |