CN107315880A - The localization method of tunnel straight flange wall three-dimensional failure mode under action of horizontal seismic - Google Patents

The localization method of tunnel straight flange wall three-dimensional failure mode under action of horizontal seismic Download PDF

Info

Publication number
CN107315880A
CN107315880A CN201710520178.2A CN201710520178A CN107315880A CN 107315880 A CN107315880 A CN 107315880A CN 201710520178 A CN201710520178 A CN 201710520178A CN 107315880 A CN107315880 A CN 107315880A
Authority
CN
China
Prior art keywords
mrow
msup
mfrac
msub
dimensional
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201710520178.2A
Other languages
Chinese (zh)
Other versions
CN107315880B (en
Inventor
安永林
胡文轩
岳健
欧阳鹏博
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hunan University of Science and Technology
Original Assignee
Hunan University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hunan University of Science and Technology filed Critical Hunan University of Science and Technology
Priority to CN201710520178.2A priority Critical patent/CN107315880B/en
Publication of CN107315880A publication Critical patent/CN107315880A/en
Application granted granted Critical
Publication of CN107315880B publication Critical patent/CN107315880B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/13Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/10Constructive solid geometry [CSG] using solid primitives, e.g. cylinders, cubes
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06Power analysis or power optimisation

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Computational Mathematics (AREA)
  • Civil Engineering (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Structural Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Architecture (AREA)
  • Computer Graphics (AREA)
  • Software Systems (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Secondary Cells (AREA)

Abstract

The invention discloses a kind of localization method of tunnel straight flange wall three-dimensional failure mode under action of horizontal seismic.The present invention is mainly:According to tunnel location seismic fortification intensity, tunnel horizontal earthquake action acceleration is obtained;According to country rock level condition, the related mechanics parameter of country rock is obtained;Calculate the interior energy dissipated power at straight flange wall Three-Dimensional fracture face;Straight flange wall three-dimensional is calculated to cave in the external force power that internal earthquake load done;Build the functional containing earthquake external force power and interior energy dissipated power;The condition of extreme value is sought according to the calculus of variations, Three-Dimensional fracture face shape function is determined;According to displacement, stress boundary condition and the conservation of energy, the coefficient value of Three-Dimensional fracture face shape function is determined;According to plane of fracture shape function and its coefficient value, drawing three-dimensional plane of fracture shape graph.The present invention provides computational methods for the three-dimensional failure mode of tunnel straight flange wall under varying level geological process;According to failure mode, reference can be provided for earthquake territory tunnel reinforcement scope.

Description

The localization method of tunnel straight flange wall three-dimensional failure mode under action of horizontal seismic
Technical field
The invention belongs to the determination technology of tunnel damage envelope under earthquake, and in particular to tunnel under a kind of action of horizontal seismic The localization method of the three-dimensional failure mode of road straight flange wall.
Background technology
China is one of many earthquake countries in the world, has many big and medium-sized cities to meet with the attack of violent earthquake, earthquake The frequency of generation is high, intensity is big, and the disaster caused is serious, such as Tangshan Earthquake, Wenchuan earthquake.Abroad such as Japan, equally Occurred multiple violent earthquake.In earthquake, structure is built by heavy damage in the underground such as some tunnels, subway, most of to need again Repair and reinforce.As can be seen here, seismic force can not be ignored for the influence of tunnel structure.For subway station, servant Row of channels, underground vehicle driving channel, subterranean commercial area etc., frequently with the form of straight flange wall, therefore how tunnel is straight under Study of Seismic The security of abutment wall seems extremely important.It is mainly research its dynamic response characteristics, antidetonation to the seismic study in tunnel at present to subtract Shake measure etc., and when occurring for earthquake, whether country rock caves in and collapse range has the research such as much seldom, this direct shadow Ring the formulation of maintenance and reinforcement scope and rational approach after the escape safety and calamity of personnel under earthquake.
The content of the invention
It is an object of the invention to provide a kind of positioning of tunnel straight flange wall three-dimensional failure mode under action of horizontal seismic Method.The present invention can provide reference according to failure mode for earthquake territory Tunnel Repair reinforcing scope;It also can determine whether earthquake in addition Whether lower abutment wall is safe;And calculate the volume of abutment wall destruction landslide under earthquake.
The localization method of the three-dimensional failure mode of tunnel straight flange wall, mainly includes as follows under the action of horizontal seismic of the present invention The step of order:
(1) according to tunnel location seismic fortification intensity, tunnel horizontal earthquake action acceleration is obtained;
(2) according to country rock level condition, the related mechanics parameter of country rock is obtained;
(3) the interior energy dissipated power at straight flange wall Three-Dimensional fracture face is calculated;
(4) straight flange wall three-dimensional is calculated to cave in the external force power that internal earthquake load done;
(5) functional containing earthquake external force power and interior energy dissipated power is built;
(6) condition of extreme value is sought according to the calculus of variations, Three-Dimensional fracture face shape function is determined;
(7) according to displacement, stress boundary condition and conservation of energy principle, the coefficient of Three-Dimensional fracture face shape function is determined Value;
(8) according to plane of fracture shape function and its coefficient value, drawing three-dimensional plane of fracture shape graph.
Its detailed process is as follows:
(1) the interior energy dissipated power in straight flange wall Three-Dimensional fracture face is calculated by below equation:
In formula:EDFor interior energy dissipated power;L is the half of plane of fracture projected length on abutment wall;σciIt is anti-for complete country rock Compressive Strength;σtmFor the tensile strength of country rock;A, B are Analysis of Field Geotechnical Parameters, can be obtained according to related form is consulted in step (2);f (x) it is plane of fracture shape bus function, f'(x) tangent slope that is f (x), i.e. first derivative;X is plane of fracture shape bus letter X coordinate value in number f (x);V is the discontinuous velocity on the plane of fracture;
(2) the external force power that earthquake load is done is:
In formula:ηcFor comprehensive effect coefficient, the open cut tunnel of rock foundation uses 0.2, other to use 0.25;A is horizontal earthquake Acceleration is acted on, is determined by step (1);G is acceleration of gravity;γ is country rock severe;It is convenient to derive, order
(3) functional by earthquake external force power and interior energy dissipated power structure is as follows:
In formula:ξ is the difference of earthquake external force power and interior energy dissipated power;
Referred to as functional;
(4) when the functional in step (3) has extreme value, then by the variation principle of functional, its corresponding Euler's square can be obtained Cheng Wei:
Solution can be obtained:
In formula:ψ is ψ [f (x), f'(x), x];c2Be Eulerian equation is solved and is integrated obtained by constant coefficient;
(5) coefficient c in step (4)2Determined by following steps:
(a) from symmetry, at the top of the plane of fracture, its slope is 0, i.e. f'(x)=0,
Can be in the hope of:c2=0;
So as to:
In formula:c3For function f'(x) coefficient obtained by integration;
(b) by geometrical condition:
It can obtain,
In formula:H is the depth of abutment wall damage envelope;c3For function f'(x) obtained by integration in coefficient, with step (5) (a) C3Unanimously;
(c) by conservation of energy principle, i.e., external force power is equal with interior energy dissipated power, can obtain:
Formula can be tried to achieve in (b) and (c) in simultaneous step (5):
So as to:
(d) plane of fracture shape bus function f (x) rotates about the z axis, you can obtain Three-Dimensional fracture surface function:
Further, with reference to the Three-Dimensional fracture surface function z=f (x, y) tried to achieve, and country rock relevant parameter and tunnel place are combined The horizontal ground motion acceleration in area, you can draw the plane of fracture shape, so that it is determined that under earthquake abutment wall reinforcing scope;And root According to scope 2L of the plane of fracture on abutment wall, height, width with the actual abutment wall in tunnel are compared, if peripheral wall height, width> 2L, then under earthquake, abutment wall will not be destroyed;Otherwise abutment wall can be destroyed, so as to assess the security of abutment wall under earthquake;Further By integration, then the volume destroyed under earthquake can be obtained.
The inventive method is determines under different earthquake grade, i.e., the three-dimensional destruction shape of tunnel straight flange wall under varying level earthquake Shape provides computational methods;According to failure mode, reference can be provided for Tunnel Repair reinforcing scope after the calamity of earthquake territory;In addition It can determine whether whether abutment wall is safe under earthquake;And calculate the volume that abutment wall collapses under earthquake.The method of the present invention, can also be applied to Subway station under geological process, servant's row or vehicle driving passage, subterranean commercial area etc. there is the hypogee structure of straight flange wall Three-dimensional damage envelope determination, reinforcing scope determination, abutment wall security evaluation etc..
Brief description of the drawings
Fig. 1 is the principle schematic of the inventive method.
In figure:F (x) is plane of fracture shape bus function;L is the half of plane of fracture projected length on abutment wall;τnTo enclose Shear stress at the rock plane of fracture, i.e. shearing strength;σnFor the direct stress at the country rock plane of fracture, i.e. normal stress;θ is plane of fracture shape The inclination angle of shape bus function f (x) slopes;A is horizontal earthquake action acceleration;H is the depth of abutment wall damage envelope.
Embodiment
The present invention is described in further detail below in conjunction with the accompanying drawings.
It is the original of the localization method of the three-dimensional failure mode of tunnel straight flange wall under action of horizontal seismic of the present invention referring to Fig. 1 Manage schematic diagram.The inventive method is mainly included following steps in sequence:
(1) according to tunnel location seismic fortification intensity, tunnel horizontal earthquake action acceleration is obtained;
(2) according to country rock level condition, the related mechanics parameter of country rock is obtained;
(3) the interior energy dissipated power at straight flange wall Three-Dimensional fracture face is calculated;
(4) straight flange wall three-dimensional is calculated to cave in the external force power that internal earthquake load done;
(5) functional containing earthquake external force power and interior energy dissipated power is built;
(6) condition of extreme value is sought according to the calculus of variations, Three-Dimensional fracture face shape function is determined;
(7) according to displacement, stress boundary condition and conservation of energy principle, the coefficient of Three-Dimensional fracture face shape function is determined Value;
(8) according to plane of fracture shape function and its coefficient value, drawing three-dimensional plane of fracture shape graph.
Its detailed process is as follows:
(1) the interior energy dissipated power in straight flange wall Three-Dimensional fracture face is calculated by below equation:
In formula:EDFor interior energy dissipated power;L is the half of plane of fracture projected length on abutment wall;σciFor complete country rock rock The compression strength of stone;σtmFor the tensile strength of country rock;A, B are Analysis of Field Geotechnical Parameters, can be obtained according to related form is consulted in step (2) ;F (x) is plane of fracture shape bus function, f'(x) tangent slope that is f (x), i.e. first derivative;X is that plane of fracture shape is female X coordinate value in line function f (x);V is the discontinuous velocity on the plane of fracture;
(2) the external force power that earthquake load is done is:
In formula:ηcFor comprehensive effect coefficient, the open cut tunnel of rock foundation uses 0.2, other to use 0.25;A is horizontal earthquake Acceleration is acted on, is determined by step (1);G is acceleration of gravity;γ is country rock severe;It is convenient to derive, order
(3) functional by earthquake external force power and interior energy dissipated power structure is as follows:
In formula:ξ is the difference of earthquake external force power and interior energy dissipated power;
Referred to as functional;
(4) when the functional in step (3) has extreme value, then by the variation principle of functional, its corresponding Euler's square can be obtained Cheng Wei:
Solution can be obtained:
In formula:ψ is ψ [f (x), f'(x), x];c2Be Eulerian equation is solved and is integrated obtained by constant coefficient;
(6) coefficient c in step (4)2Determined by following steps:
(a) from symmetry, at the top of the plane of fracture, its slope is 0, i.e. f'(x)=0,
Can be in the hope of:c2=0;
So as to:
In formula:c3For function f'(x) coefficient obtained by integration;
(b) by geometrical condition:
It can obtain,
In formula:H is the depth of abutment wall damage envelope;c3For function f'(x) obtained by integration in coefficient, with step (5) (a) C3Unanimously;
(c) by conservation of energy principle, i.e., external force power is equal with interior energy dissipated power, can obtain:
Formula can be tried to achieve in (b) and (c) in simultaneous step (5):
So as to:
(d) plane of fracture shape bus function f (x) rotates about the z axis, you can obtain Three-Dimensional fracture surface function:
Further, with reference to the Three-Dimensional fracture surface function z=f (x, y) tried to achieve, and country rock relevant parameter and tunnel place are combined The horizontal ground motion acceleration in area, you can draw the plane of fracture shape, so that it is determined that under earthquake abutment wall reinforcing scope;And root According to scope 2L of the plane of fracture on abutment wall, height, width with the actual abutment wall in tunnel are compared, if peripheral wall height, width> 2L, then under earthquake, abutment wall will not be destroyed;Otherwise abutment wall can be destroyed, so as to assess the security of abutment wall under earthquake;Further By integration, then the volume destroyed under earthquake can be obtained.
In case history, according to the computational methods more than present invention, actual numerical value is substituted into, you can obtain straight flange wall three-dimensional The coefficient value of plane of fracture shape function and Three-Dimensional fracture face shape function, according to plane of fracture shape function and its coefficient value, can be painted Straight flange wall plane of fracture 3D shape figure processed.

Claims (3)

1. the localization method of tunnel straight flange wall three-dimensional failure mode under a kind of action of horizontal seismic, it is characterised in that including as follows The step of order:
(1) according to tunnel location seismic fortification intensity, tunnel horizontal earthquake action acceleration is obtained;
(2) according to country rock level condition, the related mechanics parameter of country rock is obtained;
(3) the interior energy dissipated power at straight flange wall Three-Dimensional fracture face is calculated;
(4) straight flange wall three-dimensional is calculated to cave in the external force power that internal earthquake load done;
(5) functional containing earthquake external force power and interior energy dissipated power is built;
(6) condition of extreme value is sought according to the calculus of variations, Three-Dimensional fracture face shape function is determined;
(7) according to displacement, stress boundary condition and conservation of energy principle, the coefficient value of Three-Dimensional fracture face shape function is determined;
(8) according to plane of fracture shape function and its coefficient value, drawing three-dimensional plane of fracture shape graph.
2. according to claim 1 under action of horizontal seismic tunnel straight flange wall three-dimensional failure mode localization method, it is special Levy and be that detailed process is as follows:
(1) the interior energy dissipated power in straight flange wall Three-Dimensional fracture face is calculated by below equation:
<mrow> <msub> <mi>E</mi> <mi>D</mi> </msub> <mo>=</mo> <mn>2</mn> <mi>&amp;pi;</mi> <msubsup> <mo>&amp;Integral;</mo> <mn>0</mn> <mi>L</mi> </msubsup> <mo>{</mo> <mo>-</mo> <msub> <mi>&amp;sigma;</mi> <mrow> <mi>c</mi> <mi>i</mi> </mrow> </msub> <msup> <mrow> <mo>&amp;lsqb;</mo> <msup> <mi>ABf</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> </mrow> <mfrac> <mn>1</mn> <mrow> <mn>1</mn> <mo>-</mo> <mi>B</mi> </mrow> </mfrac> </msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mfrac> <mn>1</mn> <mi>B</mi> </mfrac> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>&amp;sigma;</mi> <mrow> <mi>t</mi> <mi>m</mi> </mrow> </msub> <mo>}</mo> <mi>x</mi> <mi>v</mi> <mi>d</mi> <mi>x</mi> <mo>;</mo> </mrow>
In formula:EDFor interior energy dissipated power;L is the half of plane of fracture projected length on abutment wall;σciFor complete country rock pressure resistance Degree;σtmFor the tensile strength of country rock;A, B are Analysis of Field Geotechnical Parameters, can be obtained according to related form is consulted in step (2);F (x) is Plane of fracture shape bus function, f'(x) be f (x) tangent slope, i.e. first derivative;X is plane of fracture shape bus function f (x) In x coordinate value;V is the discontinuous velocity on the plane of fracture;
(2) the external force power that earthquake load is done is:
<mrow> <mi>W</mi> <mi>e</mi> <mo>=</mo> <msub> <mi>&amp;eta;</mi> <mi>c</mi> </msub> <mfrac> <mi>a</mi> <mi>g</mi> </mfrac> <mi>&amp;gamma;</mi> <msubsup> <mo>&amp;Integral;</mo> <mn>0</mn> <mi>L</mi> </msubsup> <msup> <mi>&amp;pi;x</mi> <mn>2</mn> </msup> <msup> <mi>f</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mi>v</mi> <mi>d</mi> <mi>x</mi> <mo>=</mo> <mi>&amp;omega;</mi> <msubsup> <mo>&amp;Integral;</mo> <mn>0</mn> <mi>L</mi> </msubsup> <msup> <mi>&amp;pi;x</mi> <mn>2</mn> </msup> <msup> <mi>f</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mi>v</mi> <mi>d</mi> <mi>x</mi> <mo>;</mo> </mrow>
In formula:ηcFor comprehensive effect coefficient, the open cut tunnel of rock foundation uses 0.2, other to use 0.25;A adds for horizontal earthquake action Speed, is determined by step (1);G is acceleration of gravity;γ is country rock severe;It is convenient to derive, order
(3) functional by earthquake external force power and interior energy dissipated power structure is as follows:
<mrow> <mtable> <mtr> <mtd> <mrow> <mi>&amp;xi;</mi> <mo>=</mo> <msub> <mi>W</mi> <mi>e</mi> </msub> <mo>-</mo> <msub> <mi>E</mi> <mi>D</mi> </msub> <mo>=</mo> <mn>2</mn> <mi>&amp;pi;</mi> <msubsup> <mo>&amp;Integral;</mo> <mn>0</mn> <mi>L</mi> </msubsup> <mo>{</mo> <mo>{</mo> <msub> <mi>&amp;sigma;</mi> <mrow> <mi>c</mi> <mi>i</mi> </mrow> </msub> <msup> <mrow> <mo>&amp;lsqb;</mo> <msup> <mi>ABf</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> </mrow> <mfrac> <mn>1</mn> <mrow> <mn>1</mn> <mo>-</mo> <mi>B</mi> </mrow> </mfrac> </msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mfrac> <mn>1</mn> <mi>B</mi> </mfrac> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>&amp;sigma;</mi> <mrow> <mi>t</mi> <mi>m</mi> </mrow> </msub> <mo>}</mo> <mi>x</mi> <mo>+</mo> <mn>0.5</mn> <msup> <mi>&amp;omega;x</mi> <mn>2</mn> </msup> <msup> <mi>f</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>}</mo> <mi>v</mi> <mi>d</mi> <mi>x</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <mn>2</mn> <mi>&amp;pi;</mi> <msubsup> <mo>&amp;Integral;</mo> <mn>0</mn> <mi>L</mi> </msubsup> <mi>&amp;psi;</mi> <mo>&amp;lsqb;</mo> <mi>f</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>,</mo> <msup> <mi>f</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>,</mo> <mi>x</mi> <mo>&amp;rsqb;</mo> <mi>v</mi> <mi>d</mi> <mi>x</mi> </mrow> </mtd> </mtr> </mtable> <mo>;</mo> </mrow>
In formula:ξ is the difference of earthquake external force power and interior energy dissipated power;
Referred to as functional;
(4) when the functional in step (3) has extreme value, then by the variation principle of functional, can obtain its corresponding Eulerian equation is:
<mrow> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>&amp;psi;</mi> </mrow> <mrow> <mo>&amp;part;</mo> <mi>f</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>-</mo> <mfrac> <mi>d</mi> <mrow> <mi>d</mi> <mi>x</mi> </mrow> </mfrac> <mo>&amp;lsqb;</mo> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>&amp;psi;</mi> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <mi>f</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>&amp;rsqb;</mo> <mo>=</mo> <mn>0</mn> <mo>;</mo> </mrow> 1
Solution can be obtained:
In formula:ψ is ψ [f (x), f'(x), x];c2Be Eulerian equation is solved and is integrated obtained by constant coefficient;
(5) coefficient c in step (4)2Determined by following steps:
(a) from symmetry, at the top of the plane of fracture, its slope is 0, i.e. f'(x)=0,
Can be in the hope of:c2=0;
So as to:
In formula:c3For function f'(x) coefficient obtained by integration;
(b) by geometrical condition:
It can obtain,
In formula:H is the depth of abutment wall damage envelope;c3For function f'(x) obtained by integration in coefficient, with step (5) (a) c3 Unanimously;
(c) by conservation of energy principle, i.e., external force power is equal with interior energy dissipated power, can obtain:
<mrow> <mo>-</mo> <msub> <mi>&amp;pi;&amp;sigma;</mi> <mrow> <mi>t</mi> <mi>m</mi> </mrow> </msub> <msup> <mi>L</mi> <mn>2</mn> </msup> <mo>+</mo> <mfrac> <mrow> <mn>2</mn> <mi>&amp;pi;</mi> <mi>B</mi> </mrow> <mrow> <mn>2</mn> <mi>B</mi> <mo>+</mo> <mn>1</mn> </mrow> </mfrac> <msup> <mi>A</mi> <mrow> <mo>-</mo> <mn>1</mn> <mo>/</mo> <mi>B</mi> </mrow> </msup> <msup> <msub> <mi>&amp;sigma;</mi> <mrow> <mi>c</mi> <mi>i</mi> </mrow> </msub> <mfrac> <mrow> <mi>B</mi> <mo>-</mo> <mn>1</mn> </mrow> <mi>B</mi> </mfrac> </msup> <msup> <mrow> <mo>(</mo> <mn>0.5</mn> <mi>&amp;omega;</mi> <mo>)</mo> </mrow> <mrow> <mn>1</mn> <mo>/</mo> <mi>B</mi> </mrow> </msup> <msup> <mi>L</mi> <mfrac> <mrow> <mn>1</mn> <mo>+</mo> <mn>2</mn> <mi>B</mi> </mrow> <mi>B</mi> </mfrac> </msup> <mo>=</mo> <mn>0</mn> <mo>;</mo> </mrow>
Formula can be tried to achieve in (b) and (c) in simultaneous step (5):
<mrow> <mi>L</mi> <mo>=</mo> <mfrac> <mrow> <mn>2</mn> <mi>A</mi> </mrow> <mi>&amp;omega;</mi> </mfrac> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mn>1</mn> <mo>+</mo> <mn>2</mn> <mi>B</mi> </mrow> <mrow> <mn>2</mn> <mi>B</mi> </mrow> </mfrac> <mo>)</mo> </mrow> <mi>B</mi> </msup> <msup> <msub> <mi>&amp;sigma;</mi> <mrow> <mi>c</mi> <mi>i</mi> </mrow> </msub> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mi>B</mi> <mo>)</mo> </mrow> </msup> <msub> <mi>&amp;sigma;</mi> <mrow> <mi>t</mi> <mi>m</mi> </mrow> </msub> <mo>;</mo> </mrow>
<mrow> <mi>h</mi> <mo>=</mo> <mfrac> <mrow> <mn>1</mn> <mo>+</mo> <mn>2</mn> <mi>B</mi> </mrow> <mrow> <mi>&amp;omega;</mi> <mi>B</mi> </mrow> </mfrac> <msub> <mi>&amp;sigma;</mi> <mrow> <mi>t</mi> <mi>m</mi> </mrow> </msub> <mo>;</mo> </mrow>
So as to:
<mrow> <mi>f</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <msup> <mi>A</mi> <mrow> <mo>-</mo> <mn>1</mn> <mo>/</mo> <mi>B</mi> </mrow> </msup> <msup> <mrow> <mo>(</mo> <mfrac> <mi>&amp;omega;</mi> <mrow> <mn>2</mn> <msub> <mi>&amp;sigma;</mi> <mrow> <mi>c</mi> <mi>i</mi> </mrow> </msub> </mrow> </mfrac> <mo>)</mo> </mrow> <mfrac> <mrow> <mn>1</mn> <mo>-</mo> <mi>B</mi> </mrow> <mi>B</mi> </mfrac> </msup> <msup> <mi>x</mi> <mrow> <mn>1</mn> <mo>/</mo> <mi>B</mi> </mrow> </msup> <mo>-</mo> <mfrac> <mrow> <mn>1</mn> <mo>+</mo> <mn>2</mn> <mi>B</mi> </mrow> <mrow> <mi>&amp;omega;</mi> <mi>B</mi> </mrow> </mfrac> <msub> <mi>&amp;sigma;</mi> <mrow> <mi>t</mi> <mi>m</mi> </mrow> </msub> <mo>;</mo> </mrow>
(d) plane of fracture shape bus function f (x) rotates about the z axis, you can obtain Three-Dimensional fracture surface function:
<mrow> <mi>z</mi> <mo>=</mo> <mi>f</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <msup> <mi>A</mi> <mrow> <mo>-</mo> <mn>1</mn> <mo>/</mo> <mi>B</mi> </mrow> </msup> <msup> <mrow> <mo>(</mo> <mfrac> <mi>&amp;omega;</mi> <mrow> <mn>2</mn> <msub> <mi>&amp;sigma;</mi> <mrow> <mi>c</mi> <mi>i</mi> </mrow> </msub> </mrow> </mfrac> <mo>)</mo> </mrow> <mfrac> <mrow> <mn>1</mn> <mo>-</mo> <mi>B</mi> </mrow> <mi>B</mi> </mfrac> </msup> <msup> <mrow> <mo>(</mo> <msup> <mi>x</mi> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>y</mi> <mn>2</mn> </msup> <mo>)</mo> </mrow> <mrow> <mn>1</mn> <mo>/</mo> <mn>2</mn> <mi>B</mi> </mrow> </msup> <mo>-</mo> <mfrac> <mrow> <mn>1</mn> <mo>+</mo> <mn>2</mn> <mi>B</mi> </mrow> <mrow> <mi>&amp;omega;</mi> <mi>B</mi> </mrow> </mfrac> <msub> <mi>&amp;sigma;</mi> <mrow> <mi>t</mi> <mi>m</mi> </mrow> </msub> <mo>.</mo> </mrow>
3. according to claim 2 under action of horizontal seismic tunnel straight flange wall three-dimensional failure mode localization method, it is special Levy and be:With reference to the Three-Dimensional fracture surface function z=f (x, y) tried to achieve, and combine the level of country rock relevant parameter and tunnel location Earthquake ground motion acceleration, you can draw the plane of fracture shape, so that it is determined that under earthquake abutment wall reinforcing scope;And according to the plane of fracture Scope 2L on abutment wall, height, width with the actual abutment wall in tunnel are compared, if peripheral wall height, width>2L, then Under shake, abutment wall will not be destroyed;Otherwise abutment wall can be destroyed, so as to assess the security of abutment wall under earthquake;Further pass through product Point, then it can obtain the volume destroyed under earthquake.
CN201710520178.2A 2017-06-30 2017-06-30 Method for positioning three-dimensional damaged shape of straight-side wall of tunnel under action of horizontal seismic force Active CN107315880B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710520178.2A CN107315880B (en) 2017-06-30 2017-06-30 Method for positioning three-dimensional damaged shape of straight-side wall of tunnel under action of horizontal seismic force

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710520178.2A CN107315880B (en) 2017-06-30 2017-06-30 Method for positioning three-dimensional damaged shape of straight-side wall of tunnel under action of horizontal seismic force

Publications (2)

Publication Number Publication Date
CN107315880A true CN107315880A (en) 2017-11-03
CN107315880B CN107315880B (en) 2020-05-05

Family

ID=60180403

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710520178.2A Active CN107315880B (en) 2017-06-30 2017-06-30 Method for positioning three-dimensional damaged shape of straight-side wall of tunnel under action of horizontal seismic force

Country Status (1)

Country Link
CN (1) CN107315880B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109726473A (en) * 2018-12-28 2019-05-07 湖南科技大学 Based on three-dimensional deep tunnel criticality safety drilling depth and landslide method for determining shape
CN110987335A (en) * 2019-12-20 2020-04-10 防灾科技学院 Method and device for determining position of potential fracture surface of two-stage reinforced soil structure
CN111046480A (en) * 2019-12-30 2020-04-21 西南石油大学 Method for calculating seismic soil pressure time course of retaining wall based on energy dissipation principle

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPWO2007142200A1 (en) * 2006-06-09 2009-10-22 株式会社クラレ Tunnel tarpaulin
CN103469780A (en) * 2013-09-24 2013-12-25 昆明理工大学 Method of calculating deep antiskid stable critical slipping plane of dam foundation of gravity dam
CN104462835A (en) * 2014-12-16 2015-03-25 湖南科技大学 Tunnel lining safety assessing device and calculation method based on fire scenes
KR101642862B1 (en) * 2015-02-10 2016-07-26 한국지질자원연구원 System and method for predicting earthquake damage at mining area
CN106652728A (en) * 2017-02-28 2017-05-10 青岛理工大学 Method for demonstrating surface collapse disaster caused by subway tunnel construction
CN106682330A (en) * 2016-12-30 2017-05-17 湖南科技大学 Deep chamber surrounding rock pressure calculating method

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPWO2007142200A1 (en) * 2006-06-09 2009-10-22 株式会社クラレ Tunnel tarpaulin
CN103469780A (en) * 2013-09-24 2013-12-25 昆明理工大学 Method of calculating deep antiskid stable critical slipping plane of dam foundation of gravity dam
CN104462835A (en) * 2014-12-16 2015-03-25 湖南科技大学 Tunnel lining safety assessing device and calculation method based on fire scenes
KR101642862B1 (en) * 2015-02-10 2016-07-26 한국지질자원연구원 System and method for predicting earthquake damage at mining area
CN106682330A (en) * 2016-12-30 2017-05-17 湖南科技大学 Deep chamber surrounding rock pressure calculating method
CN106652728A (en) * 2017-02-28 2017-05-10 青岛理工大学 Method for demonstrating surface collapse disaster caused by subway tunnel construction

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
吕言新 等: "岩质边坡抗震稳定性的能量方法研究", 《岩石力学与工程学报》 *
姚广: "丰城矿区深部岩巷围岩应力分析及控制技术研究", 《中国优秀硕士学位论文全文数据库(工程科技Ⅰ辑)》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109726473A (en) * 2018-12-28 2019-05-07 湖南科技大学 Based on three-dimensional deep tunnel criticality safety drilling depth and landslide method for determining shape
CN109726473B (en) * 2018-12-28 2023-04-07 湖南科技大学 Three-dimensional-based method for determining critical safe footage and collapse square shape of deep-buried tunnel
CN110987335A (en) * 2019-12-20 2020-04-10 防灾科技学院 Method and device for determining position of potential fracture surface of two-stage reinforced soil structure
CN111046480A (en) * 2019-12-30 2020-04-21 西南石油大学 Method for calculating seismic soil pressure time course of retaining wall based on energy dissipation principle

Also Published As

Publication number Publication date
CN107315880B (en) 2020-05-05

Similar Documents

Publication Publication Date Title
Su et al. Characteristics and triggering mechanism of Xinmo landslide on 24 June 2017 in Sichuan, China
CN104408323A (en) Method for advanced forecasting of roof separation water disaster of stope based on multi-source information fusion
Baryakh et al. Sinkhole formation mechanism
CN103953392B (en) Rockburst risk position method of discrimination on deep tunnel section
CN107423524B (en) Method for pre-judging risk level of long-buried-depth tunnel water inrush disaster
CN107101887A (en) A kind of Numerical Investigation On Rock Failure method that sound emission is combined with numerical computations
CN103306722A (en) Detection evaluation method for microearthquake multi-dimensional information integration area of impact danger zone
CN109492262A (en) A method of utilizing numerical simulation analysis non-uniform Distribution crack Drift stability
CN107315880A (en) The localization method of tunnel straight flange wall three-dimensional failure mode under action of horizontal seismic
CN105041331B (en) A kind of tunnel portal section anti-seismic fortified length computational methods
CN108776065A (en) A kind of dynamic load prediction technique adopted induction goaf and laterally resist unstability generation stubbornly
He et al. Dynamic analysis of subway station subjected to internal blast loading
CN107339121B (en) The localization method of tunnel curl wall three-dimensional failure mode under action of horizontal seismic
Li et al. Prediction of mining subsidence in shallow coal seam
Chen et al. Study on the thin plate model with elastic foundation boundary of overlying strata for backfill mining
CN103195082A (en) Double-layer slippage side slope anti-slip retaining wall soil pressure determining method
Luo et al. Experimental study on rockfall impact force applied to frame shed tunnels
Xu et al. Study on the reasonable cover thickness of a subsea tunnel with the numerical calculation criterion method
Thoeni et al. Numerical analysis of rockfall hazard in open pit coal mines
Chen et al. Study on the support–anchor combined technique to control perilous rock at the source of avalanche by fracture mechanics
Liu et al. Numerical analysis and field monitoring tests on shallow tunnels under weak surrounding rock
Wang et al. Collapse risk evaluation method on Bayesian network prediction model and engineering application
Ma et al. The influence evaluation of the waste slag above tunnel on the safety of the structure of the high-speed railway tunnel
Li et al. Analysis of highway tunnel stability influenced by spatial distribution characteristics of concealed karst cave
Chen et al. Numerical simulation to determine the critical safe thickness of water-resisting rock wall based on FLAC3D

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant