CN112446151A - Explicit dynamic ground stress initialization method, system, medium and computer equipment - Google Patents
Explicit dynamic ground stress initialization method, system, medium and computer equipment Download PDFInfo
- Publication number
- CN112446151A CN112446151A CN202011345785.8A CN202011345785A CN112446151A CN 112446151 A CN112446151 A CN 112446151A CN 202011345785 A CN202011345785 A CN 202011345785A CN 112446151 A CN112446151 A CN 112446151A
- Authority
- CN
- China
- Prior art keywords
- model
- explicit
- stress
- initialization method
- numerical
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention belongs to the technical field of numerical simulation of underground engineering, and discloses an explicit dynamics ground stress initialization method, a system, a medium and computer equipment, wherein the explicit dynamics ground stress initialization method comprises the following steps: establishing a numerical model according to an actual underground engineering problem; corresponding dynamic parameters are given to the model according to the specific geographic environment; setting specific boundary conditions for the established numerical model; and performing explicit kinetic calculation on the model until the model is stable. The invention directly integrates the stress initialization process into the calculation flow of the explicit dynamics, greatly simplifies the traditional stress initialization method, has simple operation, is easy to set the material, the boundary and the like of the model, can finally obtain more accurate ground stress distribution condition and optimizes a key problem in the field of deep engineering numerical simulation.
Description
Technical Field
The invention belongs to the technical field of numerical simulation of underground engineering, and particularly relates to an explicit dynamics ground stress initialization method, system, medium and computer equipment.
Background
At present, deep blasting engineering often has the construction degree of difficulty big, geographical condition complicacy scheduling problem, for example the construction of tibetan railway, has very high bridge-tunnel ratio along the line, and the buried depth of a large amount of tunnels all reaches thousands of meters. Before construction, it is necessary to study the relevant engineering by numerical simulation. By utilizing the numerical simulation technology, the deep construction characteristics can be better known, and a more ideal construction scheme is designed. When numerical studies of deep ground problems are carried out, how to accurately and efficiently apply ground stress is a very critical issue.
In traditional numerical simulation, a static model or an implicit dynamic model is often established for engineering problems, then a corresponding load is applied, finally initial stress distribution under the load is obtained through calculation, and then the initial stress distribution is led into an explicit dynamic model to be calculated as an initial condition. In order to apply initial ground stress to a numerical model of a deep ground problem more simply and efficiently, a new method is urgently needed.
Through the above analysis, the problems and defects of the prior art are as follows: the existing initialization method needs to consume a large amount of time, and the parameter setting is complex and tedious.
The difficulty in solving the above problems and defects is: in view of the complicated process of the conventional stress initialization method, the stress initialization part takes about more than 60% of the time of the whole analysis process, and the requirements on the number of computer cores and memory are high, so that many related projects directly abandon the stress initialization step and directly perform analysis, and as a result, a large error is generated.
The significance of solving the problems and the defects is as follows: the invention provides an accurate and efficient dynamic stress initialization method for the field of deep engineering numerical simulation, and the stress initialization process is directly integrated into the explicit dynamic calculation process, so that the calculation time is greatly shortened and the calculation difficulty is reduced.
Disclosure of Invention
Aiming at the problems in the prior art, the invention provides an explicit dynamical ground stress initialization method, system, medium and computer equipment.
The invention is realized by an explicit dynamical stress initialization method, which comprises the following steps:
establishing a numerical model according to an actual underground engineering problem; corresponding dynamic parameters are given to the model according to the specific geographic environment;
setting specific boundary conditions for the established numerical model; and performing explicit kinetic calculation on the model until the model is stable.
Further, the kinetic parameters include, but are not limited to, density, volume model, dynamic shear modulus, and dynamic poisson's ratio of the rock.
Further, the specific boundary conditions include: transmission boundary conditions and corresponding normal stress boundary conditions.
Further, the explicit dynamics stress initialization method performs explicit dynamics calculation on the numerical model in AUTODYN until the internal stress of the model tends to be stable, and the stress initialization is successful.
Further, the explicit dynamical stress initialization method performs explicit dynamical calculation on the model until the model is stable.
It is a further object of the invention to provide a computer device comprising a memory and a processor, the memory storing a computer program which, when executed by the processor, causes the processor to perform the steps of:
establishing a numerical model according to an actual underground engineering problem; corresponding dynamic parameters are given to the model according to the specific geographic environment;
setting specific boundary conditions for the established numerical model; and performing explicit kinetic calculation on the model until the model is stable.
It is another object of the present invention to provide a computer-readable storage medium storing a computer program which, when executed by a processor, causes the processor to perform the steps of:
establishing a numerical model according to an actual underground engineering problem; corresponding dynamic parameters are given to the model according to the specific geographic environment;
setting specific boundary conditions for the established numerical model; and performing explicit kinetic calculation on the model until the model is stable.
Another object of the present invention is to provide an information data processing terminal for implementing the explicit dynamical stress initialization method.
Another object of the present invention is to provide an explicit dynamical stress initialization system implementing the explicit dynamical stress initialization method, the explicit dynamical stress initialization system comprising:
the model building module is used for building a numerical model according to the actual underground engineering problem;
the parameter assignment module is used for assigning corresponding dynamic parameters to the model according to a specific geographic environment;
the boundary condition determining module is used for setting specific boundary conditions for the established numerical model;
and the initialization module is used for performing explicit dynamics calculation on the model until the model is stable.
By combining all the technical schemes, the invention has the advantages and positive effects that: the invention directly integrates the stress initialization process into the calculation flow of the explicit dynamics, greatly simplifies the traditional stress initialization method, has simple operation, is easy to set the material, the boundary and the like of the model, can finally obtain more accurate ground stress distribution condition and optimizes a key problem in the field of deep engineering numerical simulation.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings needed to be used in the embodiments of the present application will be briefly described below, and it is obvious that the drawings described below are only some embodiments of the present application, and it is obvious for those skilled in the art that other drawings can be obtained from the drawings without creative efforts.
Fig. 1 is a flowchart of an explicit dynamical stress initialization method provided by an embodiment of the present invention.
FIG. 2 is a schematic structural diagram of an explicit dynamical stress initialization system provided by an embodiment of the present invention;
in the figure: 1. a model building module; 2. a parameter assignment module; 3. a boundary condition determining module; 4. and initializing the module.
FIG. 3 is a diagram of a numerical model provided by an embodiment of the present invention.
Fig. 4 is a graph of stress time course in a numerical model of a model stress initialization process according to an embodiment of the present invention.
FIG. 5 is a stress distribution diagram around a borehole provided by an embodiment of the present invention.
FIG. 6 is a comparison of the stress around the blast hole and the theoretical value provided by the embodiment of the invention.
Fig. 7 is a diagram of the results of numerical simulations of various models provided by embodiments of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is further described in detail with reference to the following embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
In view of the problems in the prior art, the present invention provides an explicit dynamical stress initialization method, system, medium, and computer device, and the present invention is described in detail below with reference to the accompanying drawings.
As shown in fig. 1, the explicit dynamical stress initialization method provided by the embodiment of the present invention includes the following steps:
s101, establishing a numerical model according to an actual underground engineering problem; corresponding dynamic parameters are given to the model according to the specific geographic environment;
s102, setting specific boundary conditions for the established numerical model; and performing explicit kinetic calculation on the model until the model is stable.
The explicit dynamical stress initialization method provided by the present invention can be implemented by other steps, and the explicit dynamical stress initialization method provided by the present invention of fig. 1 is only a specific embodiment.
In step S101, the dynamic parameters provided by the embodiment of the present invention include, but are not limited to, the density, the volume model, the dynamic shear modulus, and the dynamic poisson' S ratio of the rock.
In step S102, the specific boundary conditions provided by the embodiment of the present invention include: transmission boundary conditions and corresponding normal stress boundary conditions.
As shown in fig. 2, an explicit dynamical stress initialization system provided by an embodiment of the present invention includes:
the model building module 1 is used for building a numerical model according to actual underground engineering problems;
the parameter assignment module 2 is used for assigning corresponding dynamic parameters to the model according to a specific geographic environment;
a boundary condition determining module 3, configured to set a specific boundary condition for the established numerical model;
and the initialization module 4 is used for performing explicit dynamics calculation on the model until the model is stable.
The technical solution of the present invention is further described with reference to the following specific examples.
Example 1:
the invention discloses an initial method for ground stress based on an AUTODYN platform, which is characterized by comprising the following steps:
(1) establishing a numerical model according to actual engineering problems
Firstly, establishing a simplified numerical model such as a deep tunnel blasting model and the like according to actual underground engineering problems, and then endowing the model with corresponding dynamic mechanical parameters including rock density, a volume model, a dynamic shear modulus, a dynamic Poisson ratio and the like according to a specific geographic environment;
(2) setting specific boundary conditions for numerical models
According to the method, transmission boundary conditions need to be set for a numerical model to simulate the situation of an infinite large rock mass, and then corresponding normal stress boundary conditions are set according to needs, namely the stress magnitude does not change along with time;
(3) performing explicit dynamics calculation on the model until the model is stable
And performing explicit dynamic calculation on the numerical model in AUTODYN until the internal stress of the model tends to be stable, namely the fluctuation of the internal stress of the model along with time is basically ignored, which indicates that the model reaches a static equilibrium state, so that the stress initialization is successful.
Example 2:
(1) establishing a deep single-hole blasting model
The single hole blasting problem in deep blasting engineering is simplified into a numerical model as shown in fig. 3, wherein the length and the width of the numerical model are both 200mm, and the radius of a blast hole is 3.5 mm. Considering the granite geological conditions, the corresponding dynamic mechanical parameters are as follows: the density was 2.66g/cm3The volume model was 25.59GPa, the dynamic shear modulus was 21.83GPa, and the Poisson's ratio was 0.169. The model is discretized into a total of 164500 two-dimensional planar strain cells.
(2) Setting specific boundary conditions for numerical models
Setting a transmission boundary condition around the single-hole blasting model to simulate the condition of blast holes in an infinite rock body, wherein the considered normal stress boundary condition is transverse initial pressure PxTo a vertical initial pressure PyAll only 20MPa, as shown in FIG. 3.
(3) Performing explicit dynamics calculation on the model until the model is stable
Explicit dynamics calculation based on a stress difference method is carried out on the model in AUTODYN until the change of the internal stress of the model tends to be stable, as shown in FIG. 4, a circumferential stress time course curve at different positions away from a blast hole is shown, and in the later period of calculation, the stress inside the numerical model tends to be stable, which indicates that the model basically reaches a static equilibrium state, namely the stress initialization process of the model is completed.
Fig. 5 is a vertical stress distribution in the vicinity of the blast hole after the stress initialization is completed. In order to verify the accuracy of the explicit dynamic stress initialization method based on the AUTODYN platform, the theoretical solution of the problem can be obtained according to the elastic mechanics, and the analytic formulas of the hoop stress and the radial stress are as follows:
in the formula, σθAnd σrRespectively hoop stress and radial stress, PxAnd PyThe initial lateral/vertical pressure, alpha blasthole radius, and r are the distance from the center of the blasthole, respectively. The relationship between the circumferential stress and the radial stress near the blast hole and the distance from the center of the hole is plotted according to equation (1), as shown by the solid curve in fig. 6. The circumferential stress and the radial stress at the same position are obtained in the numerical model, and the results are shown as scattered points in fig. 6, so that the coincidence degree of the numerical result and the theoretical solution is good, and the method has good accuracy in stress initialization.
The invention relates to an explicit dynamics ground stress initialization method based on an AUTODYN platform, which is characterized in that a corresponding numerical model is established in software according to engineering problems, then a transmission boundary condition is applied to the model, a corresponding constant stress boundary condition is set according to specific problems, namely, the pressure does not change along with time, and then the explicit dynamics calculation is carried out on the model until the stress in the model tends to be stable. The invention directly integrates the stress initialization process into the calculation flow of the explicit dynamics, greatly simplifies the traditional stress initialization method, has simple operation, is easy to set the material, the boundary and the like of the model, can finally obtain more accurate ground stress distribution condition and optimizes a key problem in the field of deep ground engineering numerical simulation.
As shown in FIG. 7, in the standard single-hole blasting experiment, 3 numerical models are established by using explosives with the same equivalent weight, namely, the initial crustal stress model is not considered, and the initial pressure P is consideredx20MPa and PyModel 20MPa and initial pressure Px0MPa and PyThe result of blasting when the initial ground stress is considered is different from the result of ground stress not considered, and the influence rule of the initial ground stress on the fracture zone around the blast hole is basically the same as that of the previous research, namely, the cracks are developed towards the direction of larger main stress. In conclusion, the explicit dynamic stress initialization method can well reflect the influence of ground stress on the actual engineering ground, and has great significance on numerical research of deep engineering.
It should be noted that the embodiments of the present invention can be realized by hardware, software, or a combination of software and hardware. The hardware portion may be implemented using dedicated logic; the software portions may be stored in a memory and executed by a suitable instruction execution system, such as a microprocessor or specially designed hardware. Those skilled in the art will appreciate that the apparatus and methods described above may be implemented using computer executable instructions and/or embodied in processor control code, such code being provided on a carrier medium such as a disk, CD-or DVD-ROM, programmable memory such as read only memory (firmware), or a data carrier such as an optical or electronic signal carrier, for example. The apparatus and its modules of the present invention may be implemented by hardware circuits such as very large scale integrated circuits or gate arrays, semiconductors such as logic chips, transistors, or programmable hardware devices such as field programmable gate arrays, programmable logic devices, etc., or by software executed by various types of processors, or by a combination of hardware circuits and software, e.g., firmware.
The above description is only for the purpose of illustrating the present invention and the appended claims are not to be construed as limiting the scope of the invention, which is intended to cover all modifications, equivalents and improvements that are within the spirit and scope of the invention as defined by the appended claims.
Claims (9)
1. An explicit dynamical stress initialization method, comprising:
establishing a numerical model according to an actual underground engineering problem; corresponding dynamic parameters are given to the model according to the specific geographic environment;
setting specific boundary conditions for the established numerical model; and performing explicit kinetic calculation on the model until the model is stable.
2. The explicit dynamical stress initialization method of claim 1, wherein the dynamical parameters include, but are not limited to, density, volume model, dynamic shear modulus, and dynamic poisson's ratio of rock.
3. The explicit dynamical stress initialization method of claim 1, wherein the specified boundary conditions comprise: transmission boundary conditions and corresponding normal stress boundary conditions.
4. The explicit dynamical stress initialization method of claim 1, wherein the explicit dynamical stress initialization method performs explicit dynamical calculation on the numerical model in AUTODYN until the internal stress of the model is stable and the stress initialization is successful.
5. The explicit dynamical stress initialization method of claim 1 wherein the explicit dynamical stress initialization method performs explicit dynamical calculations on the model until the model stabilizes.
6. A computer device, characterized in that the computer device comprises a memory and a processor, the memory storing a computer program which, when executed by the processor, causes the processor to carry out the steps of:
establishing a numerical model according to an actual underground engineering problem; corresponding dynamic parameters are given to the model according to the specific geographic environment;
setting specific boundary conditions for the established numerical model; and performing explicit kinetic calculation on the model until the model is stable.
7. A computer-readable storage medium storing a computer program which, when executed by a processor, causes the processor to perform the steps of:
establishing a numerical model according to an actual underground engineering problem; corresponding dynamic parameters are given to the model according to the specific geographic environment;
setting specific boundary conditions for the established numerical model; and performing explicit kinetic calculation on the model until the model is stable.
8. An information data processing terminal, characterized in that the information data processing terminal is used for implementing the explicit dynamical stress initialization method of any one of claims 1 to 5.
9. An explicit dynamical stress initialization system for implementing the explicit dynamical stress initialization method of any of claims 1-3, the explicit dynamical stress initialization system comprising:
the model building module is used for building a numerical model according to the actual underground engineering problem;
the parameter assignment module is used for assigning corresponding dynamic parameters to the model according to a specific geographic environment;
the boundary condition determining module is used for setting specific boundary conditions for the established numerical model;
and the initialization module is used for performing explicit dynamics calculation on the model until the model is stable.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011345785.8A CN112446151B (en) | 2020-11-26 | 2020-11-26 | Explicit dynamic ground stress initialization method, system, medium and computer equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011345785.8A CN112446151B (en) | 2020-11-26 | 2020-11-26 | Explicit dynamic ground stress initialization method, system, medium and computer equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112446151A true CN112446151A (en) | 2021-03-05 |
| CN112446151B CN112446151B (en) | 2022-02-11 |
Family
ID=74738431
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011345785.8A Active CN112446151B (en) | 2020-11-26 | 2020-11-26 | Explicit dynamic ground stress initialization method, system, medium and computer equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112446151B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102852516A (en) * | 2012-04-19 | 2013-01-02 | 北京大学 | Full-sew-length three-dimensional crushing data simulation method and device for oil and gas reservoir development |
| CN104005750A (en) * | 2014-05-23 | 2014-08-27 | 西南石油大学 | Drilling speed increase evaluation method based on drill stem-drill bit-rock-shaft system full-scale drilling simulation |
| CN104359776A (en) * | 2014-11-14 | 2015-02-18 | 四川大学 | Testing method of I-type crack dynamic fracture toughness under explosive load |
| CN107944164A (en) * | 2017-11-30 | 2018-04-20 | 武汉科技大学 | The load numerical value separation method and system of a kind of blasting impact and transient unloading |
| CN108694280A (en) * | 2018-05-14 | 2018-10-23 | 电子科技大学 | Shock response simulation method based on novel stress hybrid tetrahedron element |
-
2020
- 2020-11-26 CN CN202011345785.8A patent/CN112446151B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102852516A (en) * | 2012-04-19 | 2013-01-02 | 北京大学 | Full-sew-length three-dimensional crushing data simulation method and device for oil and gas reservoir development |
| CN104005750A (en) * | 2014-05-23 | 2014-08-27 | 西南石油大学 | Drilling speed increase evaluation method based on drill stem-drill bit-rock-shaft system full-scale drilling simulation |
| CN104359776A (en) * | 2014-11-14 | 2015-02-18 | 四川大学 | Testing method of I-type crack dynamic fracture toughness under explosive load |
| CN107944164A (en) * | 2017-11-30 | 2018-04-20 | 武汉科技大学 | The load numerical value separation method and system of a kind of blasting impact and transient unloading |
| CN108694280A (en) * | 2018-05-14 | 2018-10-23 | 电子科技大学 | Shock response simulation method based on novel stress hybrid tetrahedron element |
Non-Patent Citations (4)
| Title |
|---|
| RUIFENG LIU等: "Study on dynamic fracture behavior of mode I crack under blasting loads", 《SOIL DYNAMICS AND EARTHQUAKE ENGINEERING》 * |
| SIYOU XIAO等: "Numerical analysis of hard rock blasting unloading effects in high in situ stress fields", 《BULL ENG GEOL ENVIRON》 * |
| 周磊等: "冲击加载下巷道内裂纹的扩展特性及破坏行为", 《爆炸与冲击》 * |
| 李萧翰等: "不同地应力下爆破振动效应分析", 《黄金科学技术》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112446151B (en) | 2022-02-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Wang et al. | A DFN–DEM multi-scale modeling approach for simulating tunnel excavation response in jointed rock masses | |
| CN106014399B (en) | Method for establishing high-precision three-dimensional ground stress model of heterogeneous stratum | |
| CN109975119B (en) | Rock double-shaft compression blasting design method | |
| Krejsa et al. | The use of the direct optimized probabilistic calculation method in design of bolt reinforcement for underground and mining workings | |
| US10634814B2 (en) | Advanced parallel “many-core” framework for reservoir simulation | |
| CN102810127A (en) | Virtual vibration test system for spacecraft | |
| EP3400546B1 (en) | Effective permeability upscaling for a discrete fracture network | |
| CN105223616A (en) | A kind of pore components inversion method of shale reservoir | |
| CN111551985A (en) | Underground karst cave detection method based on pile hammer shock and k nearest neighbor algorithm | |
| CN112446151B (en) | Explicit dynamic ground stress initialization method, system, medium and computer equipment | |
| Yu et al. | Numerical simulation on the excavation damage of Jinping deep tunnels based on the SPH method | |
| Fei et al. | A phase‐field model for hydraulic fracture nucleation and propagation in porous media | |
| Kempka et al. | Verification of TRANSPORT Simulation Environment coupling with PHREEQC for reactive transport modelling | |
| Zeng et al. | Constrained stochastic imperfection modal method for nonlinear buckling analysis of single-Layer reticulated shells | |
| CN117057110A (en) | Construction method of P-wave induced elliptic tunnel surrounding rock dynamic stress concentration coefficient calculation model | |
| Wang et al. | Influence of foundation pit excavation on tunnels at different locations | |
| Weng et al. | Influence of geostress orientation on fracture response of deep underground cavity subjected to dynamic loading | |
| Guo et al. | Numerical simulation of surface vibration propagation in tunnel blasting | |
| CN109800508A (en) | The calculation method and terminal device of the empty top plate thickness at rock-socketed piles end | |
| CN113550735B (en) | Method for calculating oil content of shale oil reservoir by using logging data | |
| Xiao et al. | Analytical solutions for the characteristic size distribution of the elliptical model in fractured rock mass | |
| Wu et al. | Strength Weakening Effect of High Prestatically Loaded Marble Subjected to Low‐Frequency Dynamic Disturbance under Point Load | |
| Hubinger et al. | A new equation solver for modeling turbulent flow in coupled matrix‐conduit flow models | |
| CN114970159A (en) | Method for calculating lateral deformation of diaphragm wall by adopting three-order shear deformation plate theory | |
| Wang et al. | Dynamic shear strength of rock joints and its influence on key blocks |
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 |