CN113515795A - Tunnel safety factor calculation method and system - Google Patents
Tunnel safety factor calculation method and system Download PDFInfo
- Publication number
- CN113515795A CN113515795A CN202110723584.5A CN202110723584A CN113515795A CN 113515795 A CN113515795 A CN 113515795A CN 202110723584 A CN202110723584 A CN 202110723584A CN 113515795 A CN113515795 A CN 113515795A
- Authority
- CN
- China
- Prior art keywords
- surrounding rock
- tunnel
- rock strength
- safety factor
- calculation
- 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.)
- Pending
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 141
- 239000011435 rock Substances 0.000 claims abstract description 142
- 238000004088 simulation Methods 0.000 claims abstract description 80
- 238000009412 basement excavation Methods 0.000 claims abstract description 34
- 230000009467 reduction Effects 0.000 claims description 27
- 239000004576 sand Substances 0.000 abstract description 8
- 231100000817 safety factor Toxicity 0.000 description 68
- 238000000034 method Methods 0.000 description 13
- 239000000463 material Substances 0.000 description 7
- 230000015572 biosynthetic process Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- 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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- 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
Abstract
The invention relates to a tunnel safety coefficient calculation method and a system thereof, comprising the following steps: setting the calculation precision delta of the safety factor; simulating tunnel excavation by using ANSYS software, outputting a simulation result, and setting a lower limit value F of a safety factor1And an upper limit value F2(ii) a Calculating a factor of safety FS＝(F1+F2) (ii)/2, judgment FSWhether greater than 1; if the judgment result is FSIf the surrounding rock strength parameter is larger than 1, reducing the surrounding rock strength parameter, simulating tunnel excavation by ANSYS software according to the reduced surrounding rock strength parameter, and outputting a simulation result; if the judgment result is FSIf the surrounding rock strength parameter is less than 1, increasing the surrounding rock strength parameter, simulating tunnel excavation by ANSYS software according to the increased surrounding rock strength parameter and outputting a simulation result(ii) a Adjusting F according to simulation result1Or F2And calculating F2‑F1A value of (d); if F2‑F1Is less than the calculation accuracy delta, F is calculatedS＝(F1+F2) /2, output safety factor FS(ii) a If F2‑F1If the value of (D) is not less than the calculation accuracy delta, the calculation of the safety factor F is restartedS＝(F1+F2) (ii)/2, judgment FSAnd whether the value is greater than 1. The invention improves the calculation precision of the safety coefficient.
Description
Technical Field
The invention relates to the technical field of tunnel engineering, in particular to a tunnel safety coefficient calculation method and a tunnel safety coefficient calculation system.
Background
In recent years, in the development period of tunnel construction in China, the aspects of tunnel length, section size, structural distribution form and the like are continuously broken through and recorded. In the design and construction of the tunnel, how to evaluate the stability of the surrounding rock of the excavated cavern is always a key problem in the field of tunnel engineering. In recent years, strength reduction method has been widely used as a means for evaluating tunnel stability as the application degree of the strength reduction method in the field of geotechnical engineering is deepened.
In the prior art, generally, a safety factor is selected manually for carrying out surrounding rock strength reduction, whether the tunnel is stable after being excavated is judged by calculation, another safety factor is selected again for carrying out surrounding rock strength reduction according to whether the tunnel is stable, the stability condition of the tunnel excavation is judged by calculation again, and the real value is continuously approached by repeating trial calculation for many times (namely the tunnel is just in a stable critical state), so that the difficulty of calculation and analysis is high, the workload is large, labor and effort are wasted, and the calculation capability is limited when the safety factor is selected manually, so that the precision of the safety factor is not high.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provides a tunnel safety coefficient calculation method and a tunnel safety coefficient calculation system to solve the problems of low calculation precision and large calculation workload of the conventional safety coefficient.
In order to achieve the above object, the present invention provides a tunnel safety factor calculation method, which comprises the following steps:
s101: setting the calculation precision delta of the safety factor;
s102: simulating tunnel by using ANSYS softwareExcavating the road and outputting a simulation result, and according to the simulation result, determining the lower limit value F of the safety factor1Set as the first set value, the upper limit value F2Setting the value as a second set value;
s103: calculating a factor of safety FS＝(F1+F2) (ii)/2, judgment FSWhether greater than 1;
s104: if the judgment result is FSIf the surrounding rock strength parameter is larger than 1, reducing the surrounding rock strength parameter, transmitting the reduced surrounding rock strength parameter to ANSYS software, and simulating tunnel excavation by the ANSYS software according to the reduced surrounding rock strength parameter and outputting a simulation result after the surrounding rock strength parameter is adjusted;
s105: if the judgment result is FSIf the surrounding rock strength parameter is less than 1, increasing the surrounding rock strength parameter, transmitting the increased surrounding rock strength parameter to ANSYS software, and simulating tunnel excavation by the ANSYS software according to the increased surrounding rock strength parameter and outputting a simulation result after the surrounding rock strength parameter is adjusted;
s106: f is adjusted according to the simulation result after the intensity parameters of the surrounding rock are adjusted1Or F2And calculating F2-F1A value of (d);
s107: if F2-F1Is less than the calculation accuracy delta, F is calculatedS＝(F1+F2) /2, output safety factor FS(ii) a And
s108: if F2-F1If the value of (d) is not less than the calculation accuracy δ, execution of step S103 is restarted.
The method utilizes ANSYS software to simulate tunnel excavation and output a simulation result, sets an upper limit value and a lower limit value of a safety coefficient according to the simulation result, and takes the average value of the upper limit value and the lower limit value of the safety coefficient as a safety coefficient FSIf F isSIf the measured value is more than 1, the tunnel in the excavated stratum model can be kept stable, the surrounding rock strength parameter is reduced and transmitted to ANSYS software, and the ANSYS software simulates tunnel excavation according to the reduced surrounding rock strength parameter and outputs a simulation result after the surrounding rock strength parameter is adjusted; if FSLess than 1, which indicates the tunnel in the excavated stratum modelIf the stability cannot be kept, increasing the surrounding rock strength parameter and transmitting the surrounding rock strength parameter to ANSYS software, wherein the ANSYS software simulates tunnel excavation according to the increased surrounding rock strength parameter and outputs a simulation result after the surrounding rock strength parameter is adjusted; f is adjusted according to the simulation result after the intensity parameters of the surrounding rock are adjusted1Or F2If F is a value of2-F1Is less than the calculation accuracy delta, F is calculatedS＝(F1+F2) /2, output FS(ii) a If F2-F1If the value of (d) is not less than the calculation accuracy δ, execution of step S103 is restarted. The method is characterized in that only in the initial stage, the calculation precision delta is set, and the calculation precision can be set to any position behind a decimal point according to the calculation requirement, so that the safety coefficient can be infinitely approximate to a true value through the thought of dichotomy, and the calculation precision of the tunnel safety coefficient is obviously improved compared with the traditional calculation method; in addition, the method does not need to provide a plurality of safety factors for repeated trial calculation, manual intervention is not needed in the calculation process, and only parameters are set at the beginning, so that the tunnel safety factors meeting the calculation precision requirement can be quickly and automatically output, and the calculation efficiency is greatly improved.
The tunnel safety factor calculation method is further improved in that the lower limit value F of the safety factor is obtained according to the simulation result1Set as the first set value, the upper limit value F2The step of setting to the second set point includes:
if the simulation result is that the tunnel is stable, the lower limit value F of the safety factor of the tunnel is set1Set to 1, upper limit F2 is set to the nominal value;
if the simulation result is that the tunnel is unstable, the lower limit value F of the safety factor of the tunnel is set1Set to 0, upper limit value F2Is set to 1.
The tunnel safety factor calculation method of the invention is further improved in that,
and reducing the surrounding rock strength parameters by using a calculation method of strength reduction.
The tunnel safety factor calculation method of the invention is further improved in that,
and increasing the surrounding rock strength parameter by using a calculation method of strength reduction.
The tunnel safety factor calculation method is further improved in that F is adjusted according to the simulation result after the intensity parameters of the surrounding rock are adjusted1Or F2The step of calculating (a) comprises:
if the simulation result after adjusting the surrounding rock strength parameters is that the tunnel is stable, F1＝FS；
If the simulation result after adjusting the surrounding rock strength parameters is that the tunnel is unstable, F2＝FS。
The invention also provides a tunnel safety coefficient calculation system, which comprises:
the parameter input module is used for setting the calculation precision delta of the safety coefficient;
and the safety factor range setting module is connected with the parameter input module and used for receiving ANSYS software simulation tunnel excavation and outputting a simulation result, and the lower limit value F of the safety factor is obtained according to the simulation result1Set as the first set value, the upper limit value F2Setting the value as a second set value;
a calculation judgment module connected with the safety factor range setting module and used for calculating FS＝(F1+F2) (ii)/2, judgment FSWhether greater than 1; if the judgment result is FSIf the surrounding rock strength parameter is larger than 1, reducing the surrounding rock strength parameter, transmitting the reduced surrounding rock strength parameter to ANSYS software, and simulating tunnel excavation by the ANSYS software according to the reduced surrounding rock strength parameter and outputting a simulation result after the surrounding rock strength parameter is adjusted; if the judgment result is FSIf the surrounding rock strength parameter is less than 1, increasing the surrounding rock strength parameter, transmitting the increased surrounding rock strength parameter to ANSYS software, and simulating tunnel excavation by the ANSYS software according to the increased surrounding rock strength parameter and outputting a simulation result after the surrounding rock strength parameter is adjusted;
a safety factor adjusting module for adjusting F connected with the calculation judging module according to the simulation result after adjusting the surrounding rock strength parameters1Or F2And calculating F2-F1A value of (d); and
is connected with the safety factor adjusting moduleThe precision judging module is used for judging whether F is required2-F1Is less than the calculation accuracy delta, if the judgment result is F2-F1Is less than the calculation accuracy delta, F is calculatedS＝(F1+F2) /2, output safety factor FS(ii) a If the judgment result is F2-F1If the value of (D) is not less than the calculation precision delta, sending an instruction for calculating and judging the safety coefficient to the calculation judgment module.
A further improvement of the tunnel safety factor calculation system of the present invention is that,
the lower limit value F of the safety factor is obtained according to the simulation result1Set as the first set value, the upper limit value F2When the set value is the second set value, the safety factor range setting module is further used for:
if the simulation result is that the tunnel is stable, the lower limit value F of the safety factor of the tunnel is set1Set to 1, upper limit F2 is set to the nominal value;
if the simulation result is that the tunnel is unstable, the lower limit value F of the safety factor of the tunnel is set1Set to 0, upper limit value F2Is set to 1.
The tunnel safety coefficient calculation system is further improved in that when the surrounding rock strength parameters are reduced, the calculation judgment module is further used for:
and reducing the surrounding rock strength parameters by using a calculation method of strength reduction.
The tunnel safety coefficient calculation system is further improved in that when the surrounding rock strength parameter is increased, the calculation judgment module is further used for:
and increasing the surrounding rock strength parameter by using a calculation method of strength reduction.
A further improvement of the tunnel safety factor calculation system of the present invention is that,
adjusting F according to the simulation result after adjusting the strength parameters of the surrounding rock1Or F2The factor of safety adjustment module is further configured to:
if the simulation result after adjusting the surrounding rock strength parameters is that the tunnel is stable, F1＝FS；
If the simulation result after adjusting the surrounding rock strength parameters is that the tunnel is unstable, F2＝FS。
Drawings
Fig. 1 is a schematic structural diagram of the tunnel safety factor calculation method 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 described in further detail below with reference to the accompanying drawings and 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. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The invention provides a tunnel safety coefficient calculation method and a tunnel safety coefficient calculation system, which are used for calculating tunnel safety coefficients. The method has the advantages that as long as the calculation precision delta is set initially, the calculation precision can be set to any position behind a decimal point according to the calculation requirement, the safety coefficient can be infinitely approximate to a true value through the thought of the dichotomy, and the calculation precision of the tunnel safety coefficient is obviously improved compared with the traditional calculation method; in addition, the method does not need to provide a plurality of safety factors for repeated trial calculation, manual intervention is not needed in the calculation process, and only parameters are set at the beginning, so that the tunnel safety factors meeting the calculation precision requirement can be quickly and automatically output, and the calculation efficiency is greatly improved.
The tunnel safety factor calculation method and the system thereof of the present invention are described below with reference to the accompanying drawings.
Referring to fig. 1, in this embodiment, the present invention provides a tunnel safety factor calculation method, including the following steps:
s101: and setting the calculation precision delta of the safety factor.
S102: simulating tunnel excavation by using ANSYS software, outputting a simulation result, and obtaining the lower limit value F of the safety factor according to the simulation result1Set as the first set value, the upper limit value F2Set to the second set value.
Preferably, the ANSYS software is used for establishing a model and calculating and judging whether the tunnel is stable, the ANSYS software is used for establishing a geometric model containing the formation geometric parameters and the tunnel geometric parameters through an APDL language, then the formation material attributes are given, at the moment, the formation adopts an elastic-plastic material model, corresponding strength parameter values (cohesive force, internal friction angle and the like) are input, and then the grid is divided to obtain an initial formation model. And calculating a stratum initial stress field by applying boundary conditions and a gravity field, and simulating tunnel excavation on the basis of the stratum initial stress field to obtain a correspondingly excavated stratum model. S103: calculating a factor of safety FS＝(F1+F2) (ii)/2, judgment FSWhether greater than 1.
S104: if the judgment result is FSAnd if the measured value is more than 1, reducing the surrounding rock strength parameter, transmitting the reduced surrounding rock strength parameter to ANSYS software, and simulating tunnel excavation by the ANSYS software according to the reduced surrounding rock strength parameter and outputting a simulation result after the surrounding rock strength parameter is adjusted.
S105: if the judgment result is FSAnd if the measured value is less than 1, increasing the surrounding rock strength parameter, transmitting the increased surrounding rock strength parameter to ANSYS software, and simulating tunnel excavation by the ANSYS software according to the increased surrounding rock strength parameter and outputting a simulation result after the surrounding rock strength parameter is adjusted.
S106: f is adjusted according to the simulation result after the intensity parameters of the surrounding rock are adjusted1Or F2And calculating F2-F1The value of (c).
S107: if F2-F1Is less than the calculation accuracy delta, F is calculatedS＝(F1+F2) /2, output safety factor FS。
S108: if F2-F1If the value of (d) is not less than the calculation accuracy δ, execution of step S103 is restarted.
In the embodiment, tunnel excavation is simulated by using ANSYS software, a simulation result is output, an upper limit value and a lower limit value of the safety factor are set according to the simulation result, and the middle value of the upper limit value and the lower limit value of the safety factor is taken as a safety factor FSIf, ifFSIf the measured value is more than 1, the tunnel in the excavated stratum model can be kept stable, the surrounding rock strength parameter is reduced and transmitted to ANSYS software, and the ANSYS software simulates tunnel excavation according to the reduced surrounding rock strength parameter and outputs a simulation result after the surrounding rock strength parameter is adjusted; if FSIf the measured value is less than 1, indicating that the tunnel in the excavated stratum model cannot be kept stable, increasing the surrounding rock strength parameter and transmitting the surrounding rock strength parameter to ANSYS software, wherein the ANSYS software simulates tunnel excavation according to the increased surrounding rock strength parameter and outputs a simulation result after adjusting the surrounding rock strength parameter; f is adjusted according to the simulation result after the intensity parameters of the surrounding rock are adjusted1Or F2If F is a value of2-F1Is less than the calculation accuracy delta, F is calculatedS＝(F1+F2) /2, output FS(ii) a If F2-F1If the value of (d) is not less than the calculation accuracy δ, execution of step S103 is restarted. The method is characterized in that only in the initial stage, the calculation precision delta is set, and the calculation precision can be set to any position behind a decimal point according to the calculation requirement, so that the safety coefficient can be infinitely approximate to a true value through the thought of dichotomy, and the calculation precision of the tunnel safety coefficient is obviously improved compared with the traditional calculation method; in addition, the method does not need to provide a plurality of safety factors for repeated trial calculation, manual intervention is not needed in the calculation process, and only parameters are set at the beginning, so that the tunnel safety factors meeting the calculation precision requirement can be quickly and automatically output, and the calculation efficiency is greatly improved.
In the embodiment, the lower limit value F of the safety factor is obtained according to the simulation result1Set as the first set value, the upper limit value F2The step of setting to the second set point includes:
if the simulation result is that the tunnel is stable, the lower limit value F of the safety factor of the tunnel is set1The upper limit F2 is set to 1 and the rated value.
If the simulation result is that the tunnel is unstable, the lower limit value F of the safety factor of the tunnel is set1Set to 0, upper limit value F2Is set to 1.
Whether finite element static force is calculated after tunnel excavation in ANSYS softwareConvergence (generally comprising a force convergence criterion and a displacement convergence criterion, which can be automatically calculated by ANSYS after task submission) is used as a main criterion for whether the tunnel is unstable. If the calculation is converged, the tunnel is stable; if the calculation is not converged, the tunnel is not stable. Upper limit value F of initial value2Is set empirically, and preferably has a nominal value of 10.
If the judgment result is FSAnd if the intensity is larger than 1, reducing the surrounding rock intensity parameter by using an intensity reduction calculation method.
If the judgment result is FSAnd if the intensity is less than 1, increasing the intensity parameter of the surrounding rock by using an intensity reduction calculation method.
The strength reduction method comprises the following steps: the density, the elastic modulus and the Poisson ratio of the surrounding rock material are not changed, and the cohesive force of the surrounding rock material is expressed by the formula c' ═ c/FsReduction, internal friction angle according to formula
And (5) reducing. In the formula, c,
Respectively the initial (actual) cohesive force and the internal friction angle of the surrounding rock, c'),
Respectively, cohesive force and internal friction angle after strength reduction, FsThe safety factor is.
Thus, at FSWhen the intensity is larger than 1, the surrounding rock intensity parameters are continuously reduced by an intensity reduction calculation method. At FSWhen the intensity is less than 1, the intensity parameter of the surrounding rock is continuously increased by the calculation method of the intensity reduction.
Further, the step of S104 includes:
if the judgment result is FSAnd if the intensity of the surrounding rock intensity parameter is larger than 1, performing intensity reduction calculation on the surrounding rock intensity parameter, adjusting the initial stratum model according to the reduced surrounding rock intensity parameter by using ANSYS software, calculating a corresponding stratum initial stress field, simulating tunnel excavation, and calculating to judge whether the tunnel is stable.
Further, the step of S105 includes:
if the judgment result is FSAnd if the intensity is less than 1, performing intensity reduction calculation on the surrounding rock intensity parameters, adjusting the initial stratum model according to the reduced surrounding rock intensity parameters by using ANSYS software, calculating a corresponding stratum initial stress field, simulating tunnel excavation, and calculating to judge whether the tunnel is stable.
In this embodiment, F is adjusted according to the simulation result after adjusting the intensity parameters of the surrounding rock1Or F2The step of calculating (a) comprises:
if the simulation result after adjusting the surrounding rock strength parameters is that the tunnel is stable, F1＝FS。
If the simulation result after adjusting the surrounding rock strength parameters is that the tunnel is unstable, F2＝FS。
Adjusting the initial stratum model according to the reduced surrounding rock strength parameters, calculating the corresponding stratum initial stress field, simulating tunnel excavation, calculating and judging whether the tunnel is stable, and adjusting F according to whether the tunnel is stable1Or F2If the judgment result shows that the tunnel is stable, the safety coefficient value is low, and the critical state that the tunnel cavern is imminent to collapse is not reached, so the safety coefficient should be increased to enable F1＝FSTo further reduce the surrounding rock strength; if the judgment result is that the tunnel is unstable, the safety coefficient value is higher at the moment, and the tunnel cavern reaches a collapsed state, so that the safety coefficient is reduced to enable F2＝FSThereby realizing the pair F1Or F2To ensure later calculated FS＝(F1+F2) And/2, the true value can be more approximated.
The following describes the workflow of the tunnel safety factor calculation method of the present invention.
In this embodiment, the calculation accuracy δ of the safety factor of the tunnel is set, a tunnel excavation geometric model is established by adopting an ANSYS APDL command, the attribute of the ground material is defined, the ground material is set as an elastic-plastic material, a Drucker-Prager model is generally adopted in ANSYS, and then a network is performedAnd grid division, setting model boundary conditions, saving the finite element model as an initial model, and defining the name of the initial model as M1.db. And then calculating an initial stress field of the stratum, excavating the tunnel, and judging whether the tunnel is excavated stably by calculating whether the tunnel is converged through an ANSYS software static force. If the calculation is converged, the tunnel excavation can be self-stabilized, and the safety coefficient is larger than 1, so that the surrounding rock strength reduction with the safety coefficient larger than 1 is carried out; if the calculation is not converged, the tunnel cannot be self-stabilized, the safety factor is less than 1, and then the surrounding rock strength reduction with the safety factor less than 1 needs to be carried out. Setting a lower limit value F of a safety factor1And an upper limit value F2. If the initial calculation converges (tunnel stabilization), the lower limit value F1Should be greater than 1; if the initial calculation does not converge (tunnel instability), the upper limit value F2Should be less than 1, lower limit value F1Should be greater than 0. And calculating a safety factor FS which is (F1+ F2)/2. If FS>1, calling an initial model M1.db by adopting an RESUME command, reducing the surrounding rock strength parameters (cohesive force and internal friction angle), adjusting the surrounding rock strength parameters of the initial model M1.db, recalculating an initial stress field, excavating a tunnel, and calculating convergence; if FS<And 1, calling an initial model M1.db by adopting an RESUME command, recalculating an initial stress field by increasing the surrounding rock strength parameter and adjusting the surrounding rock strength parameter of the initial model M1.db, excavating a tunnel, and calculating convergence. If the calculation is convergent (tunnel stability), the safety factor value is low, at this time, the reduction coefficient upper limit value F2 is kept unchanged, the lower limit is increased, and F1 is made to be FS; if the calculation is not converged (the tunnel is unstable), the safety factor is higher, and at this time, the lower reduction coefficient limit F1 is kept unchanged, and the upper limit is reduced, so that F2 is equal to FS. And judging whether the difference between the upper limit and the lower limit of the safety factor at the moment F2-F1 meets the precision calculation requirement. If the calculated safety factor FS is satisfied, the calculated safety factor FS is (F1+ F2)/2, and a calculation result FS is output; if the calculated safety factor FS is not satisfied, the step of calculating the safety factor FS as (F1+ F2)/2 is returned, and new intensity reduction is carried out until the calculation precision requirement is satisfied.
The invention also provides a tunnel safety coefficient calculation system, which comprises:
and the parameter input module is used for setting the calculation precision delta of the safety coefficient.
And the safety factor range setting module is connected with the parameter input module and used for receiving ANSYS software simulation tunnel excavation and outputting a simulation result, and the lower limit value F of the safety factor is obtained according to the simulation result1Set as the first set value, the upper limit value F2Set to the second set value.
A calculation judgment module connected with the safety factor range setting module and used for calculating FS＝(F1+F2) (ii)/2, judgment FSWhether greater than 1; if the judgment result is FSIf the surrounding rock strength parameter is larger than 1, reducing the surrounding rock strength parameter, transmitting the reduced surrounding rock strength parameter to ANSYS software, and simulating tunnel excavation by the ANSYS software according to the reduced surrounding rock strength parameter and outputting a simulation result after the surrounding rock strength parameter is adjusted; if the judgment result is FSAnd if the measured value is less than 1, increasing the surrounding rock strength parameter, transmitting the increased surrounding rock strength parameter to ANSYS software, and simulating tunnel excavation by the ANSYS software according to the increased surrounding rock strength parameter and outputting a simulation result after the surrounding rock strength parameter is adjusted.
A safety factor adjusting module connected with the calculation judging module and used for adjusting F according to the simulation result after adjusting the surrounding rock strength parameters1Or F2And calculating F2-F1The value of (c).
A precision judging module connected with the safety factor adjusting module and used for judging whether F is required2-F1Is less than the calculation accuracy delta, if the judgment result is F2-F1Is less than the calculation accuracy delta, F is calculatedS＝(F1+F2) /2, output safety factor FS(ii) a If the judgment result is F2-F1If the value of (D) is not less than the calculation precision delta, sending an instruction for calculating and judging the safety coefficient to the calculation judgment module.
In the present embodiment, the lower limit value F of the safety factor is determined according to the simulation result1Set as the first set value, the upper limit value F2When the set value is the second set value, the safety factor range setting module is further used for:
if the simulation result is that the tunnel is stable, the lower limit value F of the safety factor of the tunnel is set1The upper limit F2 is set to 1 and the rated value.
If the simulation result is that the tunnel is unstable, the lower limit value F of the safety factor of the tunnel is set1Set to 0, upper limit value F2Is set to 1.
Further, when the surrounding rock strength parameter is reduced, the calculation and judgment module is further configured to:
and reducing the surrounding rock strength parameters by using a calculation method of strength reduction.
Furthermore, when the surrounding rock strength parameter is increased, the calculation and judgment module is further configured to:
and increasing the surrounding rock strength parameter by using a calculation method of strength reduction.
In the embodiment, F is adjusted according to the simulation result after the intensity parameter of the surrounding rock is adjusted1Or F2The factor of safety adjustment module is further configured to:
if the simulation result after adjusting the surrounding rock strength parameters is that the tunnel is stable, F1＝FS；
If the simulation result after adjusting the surrounding rock strength parameters is that the tunnel is unstable, F2＝FS。
It should be noted that the structures, ratios, sizes, and the like shown in the drawings attached to the present specification are only used for matching the disclosure of the present specification, so as to be understood and read by those skilled in the art, and are not used to limit the conditions of the present invention, so that the present invention has no technical essence, and any structural modification, ratio relationship change, or size adjustment should still fall within the scope of the present invention without affecting the efficacy and the achievable purpose of the present invention. In addition, the terms "upper", "lower", "left", "right", "middle" and "one" used in the present specification are for clarity of description, and are not intended to limit the scope of the present invention, and the relative relationship between the terms and the terms is not to be construed as a scope of the present invention.
Claims (10)
1. A tunnel safety factor calculation method is characterized by comprising the following steps:
s101: setting the calculation precision delta of the safety factor;
s102: simulating tunnel excavation by using ANSYS software, outputting a simulation result, and obtaining the lower limit value F of the safety factor according to the simulation result1Set as the first set value, the upper limit value F2Setting the value as a second set value;
s103: calculating a factor of safety FS＝(F1+F2) (ii)/2, judgment FSWhether greater than 1;
s104: if the judgment result is FSIf the surrounding rock strength parameter is larger than 1, reducing the surrounding rock strength parameter, transmitting the reduced surrounding rock strength parameter to ANSYS software, and simulating tunnel excavation by the ANSYS software according to the reduced surrounding rock strength parameter and outputting a simulation result after the surrounding rock strength parameter is adjusted;
s105: if the judgment result is FSIf the surrounding rock strength parameter is less than 1, increasing the surrounding rock strength parameter, transmitting the increased surrounding rock strength parameter to ANSYS software, and simulating tunnel excavation by the ANSYS software according to the increased surrounding rock strength parameter and outputting a simulation result after the surrounding rock strength parameter is adjusted;
s106: f is adjusted according to the simulation result after the intensity parameters of the surrounding rock are adjusted1Or F2And calculating F2-F1A value of (d);
s107: if F2-F1Is less than the calculation accuracy delta, F is calculatedS＝(F1+F2) /2, output safety factor FS(ii) a And
s108: if F2-F1If the value of (d) is not less than the calculation accuracy δ, execution of step S103 is restarted.
2. The tunnel safety factor calculation method according to claim 1, wherein the lower limit value F of the safety factor is calculated according to the simulation result1Set to the first set point, upper limitValue F2The step of setting to the second set point includes:
if the simulation result is that the tunnel is stable, the lower limit value F of the safety factor of the tunnel is set1Set to 1, upper limit value F2Set to a nominal value;
if the simulation result is that the tunnel is unstable, the lower limit value F of the safety factor of the tunnel is set1Set to 0, upper limit value F2Is set to 1.
3. The tunnel safety factor calculation method according to claim 1,
and reducing the surrounding rock strength parameters by using a calculation method of strength reduction.
4. The tunnel safety factor calculation method according to claim 1,
and increasing the surrounding rock strength parameter by using a calculation method of strength reduction.
5. The tunnel safety coefficient calculation method according to claim 1, wherein F is adjusted according to a simulation result after the intensity parameter of the surrounding rock is adjusted1Or F2The step of calculating (a) comprises:
if the simulation result after adjusting the surrounding rock strength parameters is that the tunnel is stable, F1＝FS；
If the simulation result after adjusting the surrounding rock strength parameters is that the tunnel is unstable, F2＝FS。
6. A tunnel safety factor calculation system, comprising:
the parameter input module is used for setting the calculation precision delta of the safety coefficient;
and the safety factor range setting module is connected with the parameter input module and used for receiving ANSYS software simulation tunnel excavation and outputting a simulation result, and the lower limit value F of the safety factor is obtained according to the simulation result1Set as the first set value, the upper limit value F2Setting the value as a second set value;
a calculation judgment module connected with the safety factor range setting module and used for calculating FS＝(F1+F2) (ii)/2, judgment FSWhether greater than 1; if the judgment result is FSIf the surrounding rock strength parameter is larger than 1, reducing the surrounding rock strength parameter, transmitting the reduced surrounding rock strength parameter to ANSYS software, and simulating tunnel excavation by the ANSYS software according to the reduced surrounding rock strength parameter and outputting a simulation result after the surrounding rock strength parameter is adjusted; if the judgment result is FSIf the surrounding rock strength parameter is less than 1, increasing the surrounding rock strength parameter, transmitting the increased surrounding rock strength parameter to ANSYS software, and simulating tunnel excavation by the ANSYS software according to the increased surrounding rock strength parameter and outputting a simulation result after the surrounding rock strength parameter is adjusted;
a safety factor adjusting module connected with the calculation judging module and used for adjusting F according to the simulation result after adjusting the surrounding rock strength parameters1Or F2And calculating F2-F1A value of (d); and
a precision judging module connected with the safety factor adjusting module and used for judging whether F is required2-F1Is less than the calculation accuracy delta, if the judgment result is F2-F1Is less than the calculation accuracy delta, F is calculatedS＝(F1+F2) /2, output safety factor FS(ii) a If the judgment result is F2-F1If the value of (D) is not less than the calculation precision delta, sending an instruction for calculating and judging the safety coefficient to the calculation judgment module.
7. The tunnel safety factor calculation system of claim 6,
the lower limit value F of the safety factor is obtained according to the simulation result1Set as the first set value, the upper limit value F2When the set value is the second set value, the safety factor range setting module is further used for:
if the simulation result is that the tunnel is stable, the lower limit value F of the safety factor of the tunnel is set1Set to 1, upper limit F2 is set to the nominal value;
if the simulation junctionIf the tunnel is unstable, the lower limit value F of the safety factor of the tunnel is set1Set to 0, upper limit value F2Is set to 1.
8. The tunnel safety factor calculation system of claim 6, wherein when the surrounding rock strength parameter is reduced, the calculation and judgment module is further configured to:
and reducing the surrounding rock strength parameters by using a calculation method of strength reduction.
9. The tunnel safety factor calculation system according to claim 6, wherein when the surrounding rock strength parameter is increased, the calculation and judgment module is further configured to:
and increasing the surrounding rock strength parameter by using a calculation method of strength reduction.
10. The tunnel safety factor calculation system of claim 6,
adjusting F according to the simulation result after adjusting the strength parameters of the surrounding rock1Or F2The factor of safety adjustment module is further configured to:
if the simulation result after adjusting the surrounding rock strength parameters is that the tunnel is stable, F1＝FS；
If the simulation result after adjusting the surrounding rock strength parameters is that the tunnel is unstable, F2＝FS。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110723584.5A CN113515795A (en) | 2021-06-29 | 2021-06-29 | Tunnel safety factor calculation method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110723584.5A CN113515795A (en) | 2021-06-29 | 2021-06-29 | Tunnel safety factor calculation method and system |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113515795A true CN113515795A (en) | 2021-10-19 |
Family
ID=78065974
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110723584.5A Pending CN113515795A (en) | 2021-06-29 | 2021-06-29 | Tunnel safety factor calculation method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113515795A (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104715161A (en) * | 2015-03-31 | 2015-06-17 | 江西理工大学 | Method for judging stability of goaf roof |
| CN109740183A (en) * | 2018-12-06 | 2019-05-10 | 武汉理工大学 | Tunnel tunnel face estimation of stability model and benching tunnelling method excavation height design method |
| CN110007367A (en) * | 2019-04-16 | 2019-07-12 | 中铁第四勘察设计院集团有限公司 | A kind of Analyses of Tunnel Wall Rock Stability quantitative analysis method and device |
| EP3620606A1 (en) * | 2018-09-06 | 2020-03-11 | Beijing Jiaotong University | Design method and system for tunnel anchoring system based on structural characteristic of surrounding rock |
| CN111931272A (en) * | 2020-07-16 | 2020-11-13 | 中南大学 | Equal-precision safety coefficient calculation method for slope of any scale and grid division method |
-
2021
- 2021-06-29 CN CN202110723584.5A patent/CN113515795A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104715161A (en) * | 2015-03-31 | 2015-06-17 | 江西理工大学 | Method for judging stability of goaf roof |
| EP3620606A1 (en) * | 2018-09-06 | 2020-03-11 | Beijing Jiaotong University | Design method and system for tunnel anchoring system based on structural characteristic of surrounding rock |
| CN109740183A (en) * | 2018-12-06 | 2019-05-10 | 武汉理工大学 | Tunnel tunnel face estimation of stability model and benching tunnelling method excavation height design method |
| CN110007367A (en) * | 2019-04-16 | 2019-07-12 | 中铁第四勘察设计院集团有限公司 | A kind of Analyses of Tunnel Wall Rock Stability quantitative analysis method and device |
| CN111931272A (en) * | 2020-07-16 | 2020-11-13 | 中南大学 | Equal-precision safety coefficient calculation method for slope of any scale and grid division method |
Non-Patent Citations (2)
| Title |
|---|
| 刘浩;冯冀蒙;王志勇;章慧健;仇文革;: "遗传算法在洞群二次衬砌厚度优化中的应用", 岩石力学与工程学报, no. 08 * |
| 王勇;王国欣;李金会;苏井高;赵耀;苗春阳;: "城市超小净距大断而叠层隧道群施工过程动态模拟与分析", 现代隧道技术, no. 2 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Krabbenhoft et al. | An interior‐point algorithm for elastoplasticity | |
| WO2022057273A1 (en) | Slope reinforcement retaining wall height determination method and slope reinforcement retaining wall | |
| CN110263456B (en) | Tunnel dynamic feedback analysis system based on IFC standard | |
| CN106557638B (en) | Welding mixes the method for building up of two-way transition element grid model | |
| CN109918814A (en) | The software package and application method of template design aids based on Revit software | |
| CN113515795A (en) | Tunnel safety factor calculation method and system | |
| CN115345988A (en) | Secondary error measurement edge folding BIM lightweight method based on vertex importance | |
| CN113449454B (en) | Topology optimization method of steel truss structure | |
| CN113158425B (en) | Boundary element method-based full three-dimensional fracture intersection process simulation method | |
| CN110837682A (en) | Industrial robot large arm structure optimization method based on orthogonal test | |
| CN110096806A (en) | A kind of simple algorithm of determining CONSTRUCTION OF CABLE-STAYED BRIDGE initial tensioning Suo Li | |
| CN108549765B (en) | Method for designing road longitudinal by using element line in Civil3D | |
| CN103336907B (en) | The static storage capacity quick calculation method of a kind of reservoir based on DSI interpolation technique | |
| CN105005667B (en) | The compound soil body equivalent elastic modulus of fracture grouting and Poisson's ratio computational methods | |
| CN112307547B (en) | Method for designing supporting pressure of tunnel face | |
| CN106066902B (en) | Based on growth evolutional structure Topology Optimization Method | |
| CN108021758A (en) | Rope net antenna isostension based on finite element looks for shape method | |
| CN108536899A (en) | A kind of computational methods based on fluid and structural simulation photovoltaic green-house wind resistance | |
| CN107341286A (en) | A kind of check method of BIM model elements | |
| CN114444174A (en) | Automatic optimization building algorithm | |
| CN107133379B (en) | Modeling system and method for extra-high voltage tower column | |
| CN112214823B (en) | Aqueduct modeling method based on Revit+ Dynamo | |
| CN111931273B (en) | Rapid calculation method for slope safety coefficient theoretical value | |
| Pang et al. | Existence of three solutions for a poly-Laplacian system on graphs | |
| CN116680965B (en) | FDEM acceleration method based on self-adaptive time step excavation supporting simulation |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20211019 |
|
| RJ01 | Rejection of invention patent application after publication |