CN106339545B - A kind of rock mass rheology analogy method based on discontinuous deformation analysis - Google Patents
A kind of rock mass rheology analogy method based on discontinuous deformation analysis Download PDFInfo
- Publication number
- CN106339545B CN106339545B CN201610739434.2A CN201610739434A CN106339545B CN 106339545 B CN106339545 B CN 106339545B CN 201610739434 A CN201610739434 A CN 201610739434A CN 106339545 B CN106339545 B CN 106339545B
- Authority
- CN
- China
- Prior art keywords
- sillar
- deformation
- rheology
- rock mass
- structural plane
- 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.)
- Expired - Fee Related
Links
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
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)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
A kind of rock mass rheology analogy method based on discontinuous deformation analysis, can be realized the accurate simulation to rock mass Rheological Deformation, provide foundation for the long-term running safety evaluatio of engineering.It include: that (1) chooses Rheology Constitutive Equation appropriate to describe the timeliness mechanical behavior of intact rock and structural plane under action of long-term load, and correctly reflects its stress-strain-time relationship；(2) Function Fitting is carried out to the result of indoor or in-situ rheology test, determines the rheological parameter of constitutive equation employed in step (1)；(3) building numerical value calculate in Rheological Deformation recurrence formula, calculate sillar, structural plane each Time step Rheological Deformation increment；(4) it is based on minimum potential energy principal, derives the rigidity submatrix and load increment array for considering sillar and the structural plane deformation of creep；(5) other submatrixs are integrated according to the basic skills of discontinuous deformation analysis, establish the population equilibrium equation of discontinuous deformation analysis；(6) solve population equilibrium equation, obtain the displacement and deformation component of each sillar, so acquire consider sillar and structural surface rheology characteristic, with the rock mass deformation value of holding lotus time correlation.
Description
Technical field
The present invention relates to the technical field of rock mass rheology more particularly to a kind of rock mass rheologies based on discontinuous deformation analysis
Analogy method.
Background technique
Rheology refer to substance in the case where external condition is constant, power and deformation slowly varying phenomenon at any time.Rheology
Characteristic is one of important mechanical property of rock mass.Under action of long-term load, stress-strain state, the Distortion Disciplinarian of rock mass
It changes at any time, there is apparent timeliness.Carry out rheological properties of rock research for the timeliness mechanics of explanation rock mass
Behavior is of great significance.In recent years, the construction scale of projects in China is growing, geological conditions even more complex, engineering accident
Case it is commonplace.In these cases, the rheology of rock mass often leads to underground engineering, rock foundation, rock side slope
Deng the major reason for generating large deformation or even unstable failure.Therefore, the stream of rock mass should be fully considered in the design and construction of engineering
Become characteristic.
The influence factor of rock mass rheology includes stress level, water content, temperature etc., and wherein the effect of power occupy mainly
Position.In long-term research, scholars propose a variety of description rock mass stress-strain-time relationship mathematics physics models,
That is the constitutive model of Mineral rheology.Such as Maxwell (Maxwell) model, Kelvin (Kelvin) model, Burgers
(Burger) model, Xiyuan Model model, these models describe rock mass from multiple angles such as elasticity, viscoplasticity, viscoelastic plasticities
Rheological behaviour, be able to reflect rock mass in the deformation behaviour in different rheology stages, Accurate Prediction carried out to development of deformation trend.
In the 1980s, doctor Shi Genhua proposes a kind of numerical method based on non-continuum mechanics --- non-company
Continuous deformation analysis (Discontinuous Deformation Analysis, abbreviation DDA).This method is by natural joint fissure
Any block being cut into is basic unit, with the rigid body displacement of unit and is deformed into fundamental unknown variables, based on a set of efficient
Contacting search method can obtain any possible contact form in block system.Meanwhile this method has fully considered between block
Interaction is established population equilibrium equation based on minimum potential energy principal, and is solved using implicit solution.Discontinuous Deformation point
Analysis method has geneogenous advantage when simulating the contact problems of block system
However, existing Discontinuous Deformation Analysis Method can not consider the rheological behavior of rock mass, for the length of rock mass engineering project
Phase stability problem is difficult to provide accurate answer.
Summary of the invention
To overcome the shortcomings of existing technologies, Discontinuous Deformation is based on the technical problem to be solved in the present invention is to provide a kind of
The rock mass rheology analogy method of analysis, can be realized the accurate simulation to rock mass Rheological Deformation, be the long-term running peace of engineering
Full property evaluation provides foundation.
The technical scheme is that this rock mass rheology analogy method based on discontinuous deformation analysis, this method packet
Include following steps:
(1) constitutive equation appropriate is chosen to describe the timeliness mechanics of intact rock and structural plane under action of long-term load
Behavior, and correctly reflect its stress-strain-time relation；
(2) Function Fitting is carried out to the result of indoor or in-situ rheology test, determines this structure side employed in step (1)
The rheological parameter of journey；
(3) building numerical value calculate in Rheological Deformation recurrence formula, calculate sillar, structural plane each Time step rheology
Incremental deformation；
(4) it is based on minimum potential energy principal, derives and considers that the rigidity submatrix and load of sillar and the structural plane deformation of creep increase
Measure array；
(5) other submatrixs are integrated according to the basic skills of discontinuous deformation analysis, establish discontinuous deformation analysis
Population equilibrium equation；
(6) the population equilibrium equation for solving discontinuous deformation analysis, obtains the displacement and deformation component of each sillar, Jin Erqiu
Must consider sillar and structural surface rheology characteristic, with the rock mass deformation value of holding lotus time correlation.
The constitutive model of rock mass rheology is added in DDA method by the present invention, considers the stream of sillar and structural plane respectively
Become characteristic, to simulate deformation rule and failure mechanism of the rock mass under action of long-term load, therefore when can be realized to rock mass
The accurate simulation for imitating mechanical behavior, provides foundation for the long-term running safety evaluatio of engineering.
Detailed description of the invention
Fig. 1 is creep strain function curve of the rock mass under normal stress.
Fig. 2 is the schematic diagram that rock mass discontinuity generates tangential deformation.
Fig. 3 is the DDA analysis process block diagram for considering rheological properties of rock.
Specific embodiment
This rock mass rheology analogy method based on discontinuous deformation analysis, method includes the following steps:
(1) constitutive equation appropriate is chosen to describe the timeliness mechanics of intact rock and structural plane under action of long-term load
Behavior, and correctly reflect its stress-strain-time relation；
(2) Function Fitting is carried out to the result of indoor or in-situ rheology test, determines this structure side employed in step (1)
The rheological parameter of journey；
(3) building numerical value calculate in Rheological Deformation recurrence formula, calculate sillar, structural plane each Time step rheology
Incremental deformation；
(4) it is based on minimum potential energy principal, derives and considers that the rigidity submatrix and load of sillar and the structural plane deformation of creep increase
Measure array；
(5) other submatrixs are integrated according to the basic skills of discontinuous deformation analysis, establish discontinuous deformation analysis
Population equilibrium equation；
(6) the population equilibrium equation for solving discontinuous deformation analysis, obtains the displacement and deformation component of each sillar, Jin Erqiu
Must consider sillar and structural surface rheology characteristic, with the rock mass deformation value of holding lotus time correlation.
The constitutive model of rock mass rheology is added in Discontinuous Deformation Analysis Method by the present invention, consider respectively sillar and
The rheological behavior of structural plane, to simulate deformation rule and failure mechanism of the rock mass under action of long-term load, therefore can be real
Now to the accurate simulation of rock mass timeliness mechanical behavior, foundation is provided for the long-term running safety evaluatio of engineering.
Preferably, the rheology of sillar and structural plane is described respectively with formula (1), (2):
Sillar:
Structural plane:
Wherein:For the creep strain of sillar, Δ σbFor the stress variation of sillar；φbi(σb) it be stress level is σbWhen
Creep coefficient, rbiFor rate of convergence parameter, τ is Δ σbThe initial time of effect, (t- τ) are to hold the lotus time；
εf cFor the creep strain of structural plane, Δ σfFor the stress variation of structural plane；φfi(σf) it be structural plane stress is σfWhen
Creep coefficient, rfiFor rate of convergence parameter, τ is Δ σfThe initial time of effect, (t- τ) are to hold the lotus time；
N, M is parameter required for being fitted to experimental result, and taking N, M is 1~3 time.
Preferably, according to the rheological experiment of sillar and structural plane as a result, found out in the way of Function Fitting formula (1),
(2) each term coefficient in: N, M, φbi(σb)、φfi(σf)、rbi、rfi。
Preferably, consider the unitary elasticity stiffness matrix of sillar creep are as follows:
Wherein: S is the area of Rigid Body Element, and E, ν are the elasticity modulus and Poisson's ratio of Rigid Body Element, and qnb is the compacted of sillar
Variation is calculated by formula (5):
Preferably, the load increment array that the deformation of creep of sillar generates are as follows:
Wherein
WhereinFor the creep strain of sillar, can be acquired by formula (8):
Wherein Δ σnx, Δ σny, Δ τnFor the stress increment of sillar；It is acquired by formula (9):
Wherein S3With stress increment Δ σnWith current shearing strengthRatio it is related.
Preferably, the contact stiffness of rock mass discontinuity creep are as follows:
The load increment array of generation are as follows:
Wherein:
Wherein: p is that normal direction contacts modulus, and qn is the creep degree of structural plane,For the creep strain increment of structural plane, two
The calculation method and the creep degree of sillar, creep strain of person: qnb,It is identical；I, j are the global number of two contact blocks,
T(i)(x, y) indicates the shape function of the point (x, y) in unit i, (x in formula0,y0), (x1,y1), (x2,y2), (x3,y3) it is respectively a little
P0, P1, P2, P3Coordinate；As shown in Fig. 2, point P0It is point P1In side P2P3The contact point of upper hypothesis, shear spring is in P2P3Side
To with vertex P1And P0Connection, point P1Belong to block i, point P0, P2, P3Belong to block j.
The above is only presently preferred embodiments of the present invention, is not intended to limit the present invention in any form, it is all according to
According to technical spirit any simple modification, equivalent change and modification to the above embodiments of the invention, still belong to the present invention
The protection scope of technical solution.
Claims (6)
1. a kind of rock mass rheology analogy method based on discontinuous deformation analysis, it is characterised in that: method includes the following steps:
(1) Rheology Constitutive Equation appropriate is chosen to describe the timeliness mechanics of intact rock and structural plane under action of long-term load
Behavior, and correctly reflect its stress-strain-time relation；
(2) Function Fitting is carried out to the result of indoor or in-situ rheology test, determines constitutive equation employed in step (1)
Rheological parameter；
(3) building numerical value calculate in Rheological Deformation recurrence formula, calculate sillar, structural plane each Time step Rheological Deformation
Increment；
(4) it is based on minimum potential energy principal, derives rigidity submatrix and the load increment column for considering sillar and the structural plane deformation of creep
Battle array；
(5) other submatrixs other than the rigidity submatrix are collected according to the basic skills of discontinuous deformation analysis
At, establish it is discontinuous denaturation analysis population equilibrium equation；
(6) population equilibrium equation is solved, the displacement and deformation component of each sillar are obtained, and then acquires and considers sillar and structure surface current
Become characteristic, with hold the rock mass deformation value of lotus time correlation.
2. the rock mass rheology analogy method according to claim 1 based on discontinuous deformation analysis, it is characterised in that:
The rheology of sillar and structural plane is described respectively with formula (1), (2):
Sillar:
Structural plane:
Wherein:For the creep strain of sillar, Δ σbFor the stress variation of sillar；φbi(σb) it be stress level is σbWhen creep
Coefficient, rbiFor rate of convergence parameter, τ is Δ σbThe initial time of effect, (t- τ) are to hold the lotus time；
For the creep strain of structural plane, Δ σfFor the stress variation of structural plane；φfi(σf) it be structural plane stress is σfWhen it is compacted
Variable coefficient, rfiFor rate of convergence parameter, τ is Δ σfThe initial time of effect, (t- τ) are to hold the lotus time；
N, M is parameter required for being fitted to experimental result, and taking N, M is 1~3 time.
3. the rock mass rheology analogy method according to claim 2 based on discontinuous deformation analysis, it is characterised in that: according to
The rheological experiment of sillar and structural plane in the way of Function Fitting as a result, find out formula (1), each term coefficient in (2): N, M,
φbi(σb)、φfi(σf)、rbi、rfi。
4. the rock mass rheology analogy method according to claim 3 based on discontinuous deformation analysis, it is characterised in that:
Consider the element stiffness submatrix of sillar creep are as follows:
Wherein: S is the area of Rigid Body Element, and E, ν are the elasticity modulus and Poisson's ratio of Rigid Body Element, and qnb is the creep degree of sillar,
It is calculated by formula (5):
5. the rock mass rheology analogy method according to claim 4 based on discontinuous deformation analysis, it is characterised in that: sillar
The deformation of creep generate load increment array are as follows:
Wherein
WhereinFor the creep strain of sillar, acquired by formula (8):
Wherein Δ σnx, Δ σny, Δ τnFor the stress increment of sillar；It is acquired by formula (9):
Wherein S3With stress increment Δ σnWith current shearing strengthRatio it is related.
6. the rock mass rheology analogy method according to claim 5 based on discontinuous deformation analysis, it is characterised in that: rock mass
The contact stiffness of structural plane creep are as follows:
The equivalent load increment of generation are as follows:
Wherein:
Wherein: l is side P2P3Length；P is that normal direction contacts modulus；Qn is the creep degree of structural plane,For the creep of structural plane
Strain increment, the calculation method and the creep degree of sillar, creep strain of the two: qnb,It is identical；I, j are two contact blocks
Global number, T(i)(x, y) indicates the shape function of the point (x, y) in unit i, (x in formula0,y0), (x1,y1), (x2,y2),
(x3,y3) it is respectively point P0, P1, P2, P3Coordinate；Point P0It is point P1In side P2P3The contact point of upper hypothesis, shear spring is in P2P3
Direction and vertex P1And P0Connection, point P1Belong to block i, point P0, P2, P3Belong to block j.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610739434.2A CN106339545B (en) | 2016-08-26 | 2016-08-26 | A kind of rock mass rheology analogy method based on discontinuous deformation analysis |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610739434.2A CN106339545B (en) | 2016-08-26 | 2016-08-26 | A kind of rock mass rheology analogy method based on discontinuous deformation analysis |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106339545A CN106339545A (en) | 2017-01-18 |
| CN106339545B true CN106339545B (en) | 2019-05-10 |
Family
ID=57822544
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610739434.2A Expired - Fee Related CN106339545B (en) | 2016-08-26 | 2016-08-26 | A kind of rock mass rheology analogy method based on discontinuous deformation analysis |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106339545B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109871656B (en) * | 2019-03-20 | 2023-04-25 | 东北大学 | Slope stability analysis system based on dynamic strength reduction DDA method |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101603904A (en) * | 2009-07-20 | 2009-12-16 | 大连海事大学 | The unconfined compression fluid deformation of solution environmental action rock mass proving installation |
| CN101881182A (en) * | 2010-07-14 | 2010-11-10 | 中国矿业大学(北京) | Simulation and forecast evaluation method for delaying water inrush by mine fracture structure |
| CN102230869A (en) * | 2011-03-31 | 2011-11-02 | 长江水利委员会长江科学院 | Tester of on-site rheology of rock mass and its testing method |
| CN103031832A (en) * | 2013-01-07 | 2013-04-10 | 中国水电顾问集团成都勘测设计研究院 | Hydropower station weak rock zone shear rheological test method |
| CN103942387A (en) * | 2014-04-16 | 2014-07-23 | 四川大学 | New method for building rock creep constitutive model on basis of variable fractional derivative |
| CN104655825A (en) * | 2015-03-18 | 2015-05-27 | 中国电建集团成都勘测设计研究院有限公司 | Detection method for rheological behaviors of hard, fragile and crushed rock body of dam foundation |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8412500B2 (en) * | 2007-01-29 | 2013-04-02 | Schlumberger Technology Corporation | Simulations for hydraulic fracturing treatments and methods of fracturing naturally fractured formation |
-
2016
- 2016-08-26 CN CN201610739434.2A patent/CN106339545B/en not_active Expired - Fee Related
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101603904A (en) * | 2009-07-20 | 2009-12-16 | 大连海事大学 | The unconfined compression fluid deformation of solution environmental action rock mass proving installation |
| CN101881182A (en) * | 2010-07-14 | 2010-11-10 | 中国矿业大学(北京) | Simulation and forecast evaluation method for delaying water inrush by mine fracture structure |
| CN102230869A (en) * | 2011-03-31 | 2011-11-02 | 长江水利委员会长江科学院 | Tester of on-site rheology of rock mass and its testing method |
| CN103031832A (en) * | 2013-01-07 | 2013-04-10 | 中国水电顾问集团成都勘测设计研究院 | Hydropower station weak rock zone shear rheological test method |
| CN103942387A (en) * | 2014-04-16 | 2014-07-23 | 四川大学 | New method for building rock creep constitutive model on basis of variable fractional derivative |
| CN104655825A (en) * | 2015-03-18 | 2015-05-27 | 中国电建集团成都勘测设计研究院有限公司 | Detection method for rheological behaviors of hard, fragile and crushed rock body of dam foundation |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106339545A (en) | 2017-01-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zhou et al. | An enhanced equivalent continuum model for layered rock mass incorporating bedding structure and stress dependence | |
| Kong et al. | Coupled microplane and micromechanics model for describing the damage and plasticity evolution of quasi-brittle material | |
| CN103955604B (en) | A kind of containing crackle metal gradient material Prediction model for residual strength method | |
| CN106354918B (en) | The construction method of fluid structurecoupling problem numerical simulation in a kind of hydraulic fracturing | |
| CN104179176B (en) | The slope creep amount computational methods coupled with Rock And Soil creep are lost based on prestressd anchor cable | |
| Zou et al. | An extended polygon scaled boundary finite element method for the nonlinear dynamic analysis of saturated soil | |
| Jiang et al. | Failure analysis of a cracked concrete gravity dam under earthquake | |
| Chen et al. | Composite element model for discontinuous rock masses | |
| CN106446402A (en) | Rapid simulation modeling method for multi-field coupled discrete elements of soil water loss and cracking | |
| Levkovitch et al. | Application of a non-linear confinement sensitive constitutive model to mine scale simulations subject to varying levels of confining stress | |
| CN106339545B (en) | A kind of rock mass rheology analogy method based on discontinuous deformation analysis | |
| CN108170898A (en) | A kind of jointed rock slope reliability analysis Lower Bound Limit | |
| Fan et al. | Cone complimentary-based numerical manifold method modeling frictional and cohesive contact problems | |
| CN111221039B (en) | Method, device and equipment for determining elastic wave velocity of natural gas hydrate | |
| Hrinda | Geometrically nonlinear static analysis of 3D trusses using the arc-length method | |
| Das et al. | Finite Volume Method for Structural Analysis of RF MEMS Devices Using Theory of Plates | |
| Ibrahimbegovic et al. | Fluid-structure interaction problems solution by operator split methods and efficient software development by code-coupling | |
| Fang et al. | Prediction of downburst-induced wind pressure coefficients on high-rise building surfaces using BP neural network | |
| Al Sabouni-Zawadzka et al. | On orthotropic properties of tensegrity structures | |
| Wang et al. | Tensile Fracture Simulation in Rock Using Augmented Virtual Internal Bond Method | |
| Labao et al. | Finite-element analysis of an overturned fold using a viscous-fluid model | |
| Wang et al. | Analysis of solids with numerous microcracks using the fast multipole DBEM | |
| Timalsina et al. | Computing Fluid-Structure Interaction by the Partitioned Approach with Direct Forcing | |
| Yuan et al. | A discrete element modeling of rock and soil material based on the machine learning | |
| Su et al. | Discussion of SRFEM with mohr-coulomb plasticity model in slope stability analysis |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190510 Termination date: 20210826 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |