CN105069269B - A kind of tunnel excavation Whole Process Simulation method - Google Patents
A kind of tunnel excavation Whole Process Simulation method Download PDFInfo
- Publication number
- CN105069269B CN105069269B CN201510600424.6A CN201510600424A CN105069269B CN 105069269 B CN105069269 B CN 105069269B CN 201510600424 A CN201510600424 A CN 201510600424A CN 105069269 B CN105069269 B CN 105069269B
- Authority
- CN
- China
- Prior art keywords
- creep
- msub
- time
- mrow
- msup
- 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
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a kind of tunnel excavation Whole Process Simulation method, it is characterized in that country rock body in constructing tunnel is regarded to be a kind of with time correlation structure, therefore the method that simulation creep can be used, the dynamic change of the country rock under the different construction stages is simulated, and then simulates the purpose of whole tunnel excavation process.Step is as follows：1）Tested by creep of rock and obtain Rock Creep rule；2）Creep Equation parameter is fitted based on creep experiments；3）Creep Equation is embedded in finite element software, by setting control time length, simulation tunnel excavates overall process.
Description
Technical field
The present invention relates to a kind of tunnel excavation Whole Process Simulation method.
Background technology
Existing tunnels mimic method, one kind are to regard country rock with lining cutting as elasticoplastic body, are simulated by the life and death of cell
Tunnel excavation and supporting, kill after excavating unit, directly activate lining element.This method can only obtain the knot of end-state
Fruit, can not simulate from excavate to closing stage tunnel change, can not embody country rock with just lining interact common eigenvector this
One process.Another kind is according to load-structure design model, and country rock and lining cutting interaction relationship are simulated using spring unit.
This method is not inconsistent with the New Austrian Tunneling Method in practice of construction, and can not consider third principal stress.
The content of the invention
For above shortcomings in the prior art, the invention provides a kind of tunnel excavation Whole Process Simulation method.
This method regards country rock body in constructing tunnel to be a kind of with time correlation structure, therefore the method that can use simulation creep,
The dynamic change of the country rock under the different construction stages is simulated, and then simulates the purpose of whole tunnel excavation process；This method
Comprise the following steps：
1) creep of rock is tested：To the tunnel surrounding coring to be simulated, the high 100mm of diameter 50mm cylindrical type is processed into
Test specimen, traditional uniaxial static creep test is carried out, obtained under axially different stress state, test specimen creep strain changes over time rule
Rule；
2) following Creep Equation is selected, creep results are fitted：
Wherein, εcrFor rock equivalent creep strain, t is the time, and σ is rock equivalent stress, and e is natural constant, and T is temperature
Degree,For the rock parameter obtained by experimental fit；The parameter c when not considering temperature effect4And c7It is zero, the equation
Creep starting stage and secondary creep stage can be simulated；
3) Creep Equation is embedded in ANSYS, control time length is set, simulation tunnel excavates overall process：
3.1) tunnel excavation model, and grid division are established；To include preliminary bracing and two linings in model；
3.2) country rock material passes through experimental fit rock parameter by tb, creep orders, imparting；
3.3) simulated ground stress：Apply initial field stress, by time orders, the time is set to 1e-8；Time very little herein,
It is the influence in order to ignore creep effect to crustal stress；
3.4) before excavation simulation is completed to preliminary bracing：Kill and excavate unit, but do not activate preliminary bracing, ordered by time
Order, time are set to t1, and t1 is the time of stage consumption used in practice of construction；
3.5) before simulation preliminary bracing is completed to two linings：Preliminary bracing is activated, by time orders, the time is set to t2, t2-
T1 is by the stage consumes in practice of construction time；
3.6) the simulation tunnel operation phase：The lining of activation two, by time orders, the time is set to t3, and t3-t2 is the operation in tunnel
Time.
The beneficial effects of the invention are as follows：Can simulation tunnel excavate overall process in each construction stage Change in surrounding rocks situation,
Analog result can instruct actual tunnel to construct.
Embodiment
A kind of tunnel excavation Whole Process Simulation method, this method comprise the following steps：
1) tested by creep of rock and obtain Rock Creep rule：To the tunnel surrounding coring to be simulated, it is processed into straight
The high 100mm of footpath 50mm cylindrical type test specimen, carries out traditional uniaxial static creep test, obtains under axially different stress state, test specimen
Creep strain changes over time rule.
2) the implicit Creep Equations of ANASYS are selected, creep results are fitted：
Wherein, εcrFor rock equivalent creep strain, t is the time, and σ is rock equivalent stress, and e is natural constant, and T is temperature
Degree,For equation parameter (rock parameter obtained by experimental fit)；
The parameter c when not considering temperature effect4And c7It is zero, the equation can simulate creep starting stage and secondary creep
Stage.
3) implicit Creep Equation is embedded in ANSYS, control time length is set, simulation tunnel excavates overall process：
3.1) tunnel excavation model, and grid division are established；To include preliminary bracing and two linings in model；
3.2) country rock material creep orders (i.e. material properties), is assigned creep material parameter and (intended by testing by tb
Close obtained rock parameter)；
3.3) simulated ground stress：Apply initial field stress, by time orders, the time is set to 1e-8；Time very little herein,
It is the influence in order to ignore creep effect to crustal stress；
3.4) before excavation simulation is completed to preliminary bracing：Kill and excavate unit, but do not activate preliminary bracing, ordered by time
Order, time are set to t1, and t1 is the time of stage consumption used in practice of construction；
3.5) before simulation preliminary bracing is completed to two linings：Preliminary bracing is activated, by time orders, the time is set to t2, t2-
T1 is by the stage consumes in practice of construction time；
3.6) the simulation tunnel operation phase：The lining of activation two, by time orders, the time is set to t3, and t3-t2 is the operation in tunnel
Time.
Finally illustrate, the above embodiments are merely illustrative of the technical solutions of the present invention and it is unrestricted, although with reference to compared with
The present invention is described in detail good embodiment, it will be understood by those within the art that, can be to the skill of the present invention
Art scheme is modified or equivalent substitution, and without departing from the objective and scope of technical solution of the present invention, it all should cover at this
Among the right of invention.
Claims (1)
- A kind of 1. tunnel excavation Whole Process Simulation method, it is characterised in that country rock body is regarded as a kind of and time correlation structure, Using Creep Equation, country rock dynamic change under the different construction stages is simulated, and then simulation tunnel excavates the purpose of overall process； This method specifically comprises the following steps：1) creep of rock is tested：To the tunnel surrounding coring to be simulated, the cylindrical type for being processed into the long 100mm of diameter 50mm tries Part, traditional uniaxial static creep test is carried out, obtained under axially different stress state, test specimen creep strain changes over time rule；2) creep results are fitted, obtain Creep Equation parameter, Creep Equation is as follows：<mrow> <msub> <mi>&epsiv;</mi> <mrow> <mi>c</mi> <mi>r</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mi>c</mi> <mn>1</mn> </msub> <msup> <mi>&sigma;</mi> <msub> <mi>c</mi> <mn>2</mn> </msub> </msup> <msup> <mi>t</mi> <mrow> <msub> <mi>c</mi> <mn>3</mn> </msub> <mo>+</mo> <mn>1</mn> </mrow> </msup> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mfrac> <msub> <mi>c</mi> <mn>4</mn> </msub> <mi>T</mi> </mfrac> </mrow> </msup> </mrow> <mrow> <msub> <mi>c</mi> <mn>3</mn> </msub> <mo>+</mo> <mn>1</mn> </mrow> </mfrac> <mo>+</mo> <msub> <mi>c</mi> <mn>5</mn> </msub> <msup> <mi>&sigma;</mi> <msub> <mi>c</mi> <mn>6</mn> </msub> </msup> <msup> <mi>te</mi> <mrow> <mo>-</mo> <mfrac> <msub> <mi>c</mi> <mn>7</mn> </msub> <mi>T</mi> </mfrac> </mrow> </msup> </mrow>Wherein, εcrFor rock equivalent creep strain, t is the time, and σ is rock equivalent stress, and e is natural constant, and T is temperature,For the rock parameter obtained by experimental fit；The parameter c when not considering temperature effect4And c7It is zero, the equation can Simulate creep starting stage and secondary creep stage；3) Creep Equation is embedded in ANSYS finite element softwares, control time length is set, simulation tunnel excavates overall process：3.1) tunnel excavation model, and grid division are established；To include preliminary bracing and two linings in model；3.2) country rock material creep orders, assigns the rock parameter obtained by experimental fit by tb；3.3) simulated ground stress：Apply initial field stress, by time orders, the time is set to 1e-8；Time very little herein, be for Ignore influence of the creep effect to crustal stress；3.4) before excavation simulation is completed to preliminary bracing：Kill and excavate unit, but do not activate preliminary bracing, by time orders, Time is set to t1, and t1 is the time of stage consumption used in practice of construction；3.5) before simulation preliminary bracing is completed to two linings：Preliminary bracing is activated, by time orders, the time is set to t2, and t2-t1 is The time that the stage is consumed in practice of construction；3.6) the simulation tunnel operation phase：The lining of activation two, by time orders, the time is set to t3, when t3-t2 is the operation in tunnel Between.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510600424.6A CN105069269B (en) | 2015-09-18 | 2015-09-18 | A kind of tunnel excavation Whole Process Simulation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510600424.6A CN105069269B (en) | 2015-09-18 | 2015-09-18 | A kind of tunnel excavation Whole Process Simulation method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105069269A CN105069269A (en) | 2015-11-18 |
| CN105069269B true CN105069269B (en) | 2017-12-26 |
Family
ID=54498635
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510600424.6A Expired - Fee Related CN105069269B (en) | 2015-09-18 | 2015-09-18 | A kind of tunnel excavation Whole Process Simulation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105069269B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106198208B (en) * | 2016-06-21 | 2018-09-25 | 中电建路桥集团有限公司 | A kind of interbedding of soft and hard rocks surrounding rock tunnel method for analog construction considering stratification effect |
| CN106354972B (en) * | 2016-09-18 | 2019-06-28 | 重庆大学 | Nearly horizontal interbedding of soft and hard rocks surrounding rock tunnel optimization of construction technology method |
| CN111460703B (en) * | 2020-03-10 | 2023-05-05 | 西安科技大学 | Simulation calculation method for freezing and sinking process of coal mine vertical shaft |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005061091A (en) * | 2003-08-14 | 2005-03-10 | Seisuke Okubo | Finite element analysis method using increment type constitution relational expression following creep model |
| CN104598296A (en) * | 2015-01-23 | 2015-05-06 | 陈建永 | Numerical simulation analysis method for deformation of soft soil deep foundation pit |
-
2015
- 2015-09-18 CN CN201510600424.6A patent/CN105069269B/en not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005061091A (en) * | 2003-08-14 | 2005-03-10 | Seisuke Okubo | Finite element analysis method using increment type constitution relational expression following creep model |
| CN104598296A (en) * | 2015-01-23 | 2015-05-06 | 陈建永 | Numerical simulation analysis method for deformation of soft soil deep foundation pit |
Non-Patent Citations (4)
| Title |
|---|
| A Constitutive Model and Numerical Simulation for;Ning Jianguo;《IEEE Conference Publications》;4271 - 4274;20111231;第1-4页 * |
| Creep model for Zonal disintegration of deep rock mass;Li Haitao;《National Natural Science Foudation of China》;20111231;第1-5页 * |
| Study on the Creep Characteristic and Model of;Haiping Yuan;《IEEE Conference Publications》;20111231;第1-5页 * |
| 岩石内变量蠕变模型研究;乔丽苹;《岩土力学》;20121210;第33卷(第12期);第1-10页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105069269A (en) | 2015-11-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105069269B (en) | A kind of tunnel excavation Whole Process Simulation method | |
| Sun et al. | Impacts of bedding directions of shale gas reservoirs on hydraulically induced crack propagation | |
| CN105181927B (en) | The hypotonic coal seam hydraulic fracture simulation experiment method of multi- scenarios method | |
| Chen et al. | Experimental study on low‐strength similar‐material proportioning and properties for coal mining | |
| CN104792966B (en) | A kind of testing apparatus and method simulating subway shield tunnel construction initiation cavity, stratum | |
| Ma et al. | Experimental study on damage failure mechanical characteristics and crack evolution of water-bearing surrounding rock | |
| Ozturk et al. | The estimation of uniaxial compressive strength conversion factor of trona and interbeds from point load tests and numerical modeling | |
| CN104314538A (en) | Device and method for simulating large-sized acid-etched earthworm hole | |
| CN109187237A (en) | A kind of tunnels and underground engineering Blasting Excavation model assay systems and test method | |
| Wang et al. | Numerical simulation of coal floor fault activation influenced by mining | |
| Zhang et al. | Experimental simulation and a reliable calibration method of rockfill microscopic parameters by considering flexible boundary | |
| Kamdem et al. | New description of the mechanical creep response of rocks by fractional derivative theory | |
| Sainoki et al. | Numerical simulation of rock mass vibrations induced by nearby production blast | |
| Desai et al. | Numerical methods and constitutive modelling in geomechanics | |
| Zhou et al. | Numerical analysis of dynamic responses of rock containing parallel cracks under combined dynamic and static loading | |
| Haeri et al. | Experimental and numerical investigation of the effect of sample shapes on point load index | |
| Xu et al. | The weakening effect of hydrostatic pressure on rock mass of different lithology | |
| White et al. | Modeling the dynamic flow resistance across the fracture network of EGS Collab Experiment 1 | |
| Li et al. | Numerical simulation of creep characteristic for composite rock mass with weak interlayer | |
| Spreafico et al. | Effects of poroelasticity and stress corrosion on large rock slope instability using damage-based time-dependent numerical modelling | |
| Huang et al. | The Impacts of Increased Ductility in Organic-Rich Shale on Fracture Network Growth by Hydraulic Fracturing | |
| Bilotta et al. | A cohesive zone model to simulate fatigue crack propagation under high pressure gaseous hydrogen | |
| Li et al. | Numerical Test Research on Loading Rate Effect of Rock Mass with Transfixion Sawtooth Structural Plane | |
| Ru et al. | Numerical modeling of hydraulic fracture propagation using extended finite element method | |
| Neupane et al. | Developing EGS Collab Site Geologic Framework Models at the Sanford Underground Research Facility, Lead, South Dakota |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171226 Termination date: 20180918 |