CN104992013A - Mountain tunnel anti-seismic analysis method - Google Patents
Mountain tunnel anti-seismic analysis method Download PDFInfo
- Publication number
- CN104992013A CN104992013A CN201510364671.0A CN201510364671A CN104992013A CN 104992013 A CN104992013 A CN 104992013A CN 201510364671 A CN201510364671 A CN 201510364671A CN 104992013 A CN104992013 A CN 104992013A
- Authority
- CN
- China
- Prior art keywords
- tunnel
- analysis
- mountain
- stratum
- seismic
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a mountain tunnel anti-seismic analysis method. For the aspect of cross section anti-seismic analysis of a tunnel, based on a stratal elastic shear deformation theory, horizontal stratum deformation, vertical stratum deformation and seismic self-weight inertial force of broken rock bodies at the top of the tunnel are considered at the same time, and related problems of tunnel cross section seismic resistance are studied; and for the aspect of longitudinal overall anti-seismic analysis of the tunnel, an elastic foundation beam is taken as a model, interaction of the tunnel and surrounding rocks is considered, a dynamic balance equation of a tunnel structure is created according to an anti-seismic dynamic theory and a deformation transfer theory by utilizing continuum mechanics, and a classic seismic deformation method is a particular case of the method. According to the method, the cross section anti-seismic analysis and the longitudinal overall anti-seismic analysis of the mountain tunnel are performed by adopting unified thought for the first time; the thought is novel; the method is feasible; and through example verification, the method has relatively good applicability for mountain tunnel primary design under the action of strong earthquake.
Description
Technical field
The invention belongs to tunnel aseismic analysis technical field, particularly relate to a kind of method of mountain tunnel aseismic analysis.
Background technology
Along with the deep propelling of strategy to develop western regions, western China infrastructure construction will obtain development energetically, especially highway and railway transport development obtains great-leap-forward development, certainly will to build many tunnels like this, these areas are mostly in the area of highlight lines, ensure that the function of Tunnel Engineering when earthquake is normal, dark tool strategic importance.Especially subway earthquake and China in Japan's " slope god " earthquake " " in earthquake, the serious earthquake in 109 tunnels, makes the earthquake of people to underground works have understanding clearly in 5.12 Wenchuans.
Mountain tunnel Aseismic Analytical Method mostly adopts the finite element numerical analysis technology such as ANSYS to carry out, but the numerical simulation technology of advanced person is concerning more complicated project planner, and the time and expense all higher, and find a kind of simple, efficient and applicable Near covering, the primary design for Tunnel Engineering of Mountainous Region is necessary.
At present, mountain tunnel Practical Method of Designing is many, the Fu Jiyewafa, U.S. ST.John method etc. of such as USSR (Union of Soviet Socialist Republics), but these methods all do not consider the impact of top covering rockmass and tunnel structure deadweight, in fact, mountain tunnel, due to the limitation of current operating technique, certainly will produce massif loosening at tunnel vault, these loads will produce larger Earthquake Inertia Force Acting when earthquake, and it is unsafe for ignoring these loads to mountain tunnel.And current method have employed different thought and methods when analyzing mountain tunnel xsect aseismic analysis and longitudinal aseismic analysis, imperfect in system.
Summary of the invention
The object of the embodiment of the present invention is a kind of method providing mountain tunnel aseismic analysis, be intended to solve mountain tunnel under severe earthquake action, fairly simple method and program is utilized to carry out the primary design of mountain tunnel, for project planner provides necessary design parameter and the problem of data.
The embodiment of the present invention is achieved in that a kind of method of mountain tunnel aseismic analysis, and the method for this mountain tunnel aseismic analysis comprises:
Step one, determines the geometric parameter of mountain tunnel, determines and simplifies stratum basic parameter, determining the kinetic parameter of rock mass;
Step 2, according to geological exploration data and seismic strata deformable statistical data, determines the deformation pattern on stratum and the level of stratum deformation and vertical deformation;
Step 3, mountain tunnel xsect aseismic analysis;
Step 4, mountain tunnel Horizontal vibration of piles (Horizontal Deformation pressure q
h); Mountain tunnel Vertical Vibration (comprises stratum vertical deformation pressure q
v1with loose ground Earthquake Inertia Force Acting q
v2); Utilize structural mechanics, determine the internal force of tunnel-liner;
Step 5, mountain tunnel is overall aseismic analysis longitudinally, sets up the governing equation determination displacement transfer coefficient ξ of the longitudinal overall aseismic analysis of mountain tunnel
wand ξ
a;
Step 6, determines the longitudinal earthquake stress in tunnel, considers the randomness of geological process and the uncertainty of action direction, determines the meridional stress combination of tunnel structure;
Step 7, when carrying out longitudinal aseismic analysis to mountain tunnel, according to the concrete condition of mountain tunnel, carries out analytical calculation.
Further, in step one, determine that the geometric parameter of mountain tunnel comprises: h
1-edpth of tunnel (m), h
2-tunnel-liner overall height (m), the total span (m) in B-tunnel, the average thickness (m) of δ-tunnel-liner, R
0the equivalent redius (m) in-tunnel, E
sthe elastic modulus (Gpa) of-tunnel-liner, μ
sthe Poisson ratio of-tunnel-liner;
Determine that stratum basic parameter and simplification comprise: γ=∑ γ
ih
i/ H, μ
d=0.45, v
s=H/ ∑ (h
i/ v
si), G
d=γ v
s 2/ g, E
d=2 (1-μ
d) G
d, the thickness (m) of H-top layer ground, E-soil-structure interactions elastic modulus; The Poisson ratio of μ-soil;
Determine that the kinetic parameter of rock mass comprises: K
h-top layer Foundation Design horizontal seismic coefficient, K
v-top layer Foundation Design vertical seismic action coefficient; a
v-stratum maximum vertical acceleration (m/s), T-stratum natural period, V
s-elastic shear velocity of wave, S
vthe velocity response spectrum of-vibrations reference field, K
a-horizontal resiliency foundation coefficient, K
w-vertical flat elastic foundation coefficient.
Further, in step 2, determine that the formula of the deformation pattern on stratum and the level of stratum deformation and vertical deformation is:
Further, in step 4, determine that the internal force of tunnel-liner is according to the mechanics of materials, tunnel cross sectional earthquake stress is:
Further, in step 6, determine that the meridional stress combinatorial formula of tunnel structure is:
Further, in step 7, according to the concrete condition of mountain tunnel, the analytical calculation carried out comprises:
1. structure-the analysis of country rock interaction to lining cutting meridional stress; 2. edpth of tunnel is to the analysis of lining cutting meridional stress; 3. basement rock shearing wave input direction is analyzed lining cutting meridional stress; 4. overlying strata is analyzed lining cutting meridional stress eigenperiod; 5. tunnel foundation bedding value K
a=K
wanalysis to lining cutting stress when=β G, β change; 6. the analysis to lining cutting stress during elastic modulus of surrounding rocks change.
Further, the method for this mountain tunnel aseismic analysis for the implementation method that tunnel cross sectional aseismic analysis is concrete is:
The first step, according to geological exploration data and seismic strata deformable statistical data, determine the deformation pattern on stratum and stratum deformation level and vertical deformation:
In formula: S
vvelocity response spectrum (m/s) during-unit seismic coefficient;
The basic natural period (s) of T-top layer ground;
K
hdesign level seismic coefficient (earthquake degree) on the ground of-top layer;
Second step, mountain tunnel Horizontal vibration of piles:
3rd step, mountain tunnel Vertical Vibration:
A. stratum vertical deformation pressure:
B. utilize existing highway or Railway Design specification, top, hole loose ground pressure is:
In formula: a
vfor stratum maximum vertical acceleration:
4th step, under shearing wave effect, lining cutting xsect internal force calculates:
Highway and railway tunnel many employings curve wall-lining cutting of meizoseismal area, utilizes structural mechanics and mechanics of materials knowledge to obtain internal force and the cross section stresses of tunnel lining structure:
N=q
vR sin
2θ+q
hR sin2θ
Q=-q
vR sinθcosθ-q
hR cos2θ
In formula: δ is the average thickness (m) of tunnel-liner;
5th step, utilizes reinforced concrete substantially theoretical, carries out the arrangement of reinforcement in tunnel-liner cross section.
Further, the method for this mountain tunnel aseismic analysis for the tunnel implementation method that longitudinally overall aseismic analysis is concrete is:
The first step, determines correlation parameter, sets up tunnel two dimensional motion equation under ground seismic wave function:
E is the elastic modulus of beam, and I is beam section moment of inertia, and A is beam section area, K
wfor grade beam Poisson ratio, K
afor grade beam Axial distortion parameter, the transversal displacement that ν (x, t) is tunnel, the axial displacement that u (x, t) is tunnel, g
wthe transversal displacement that (x, t) is foundation soil, g
athe axial displacement that (x, t) is foundation soil;
Second step, introduces displacement transfer coefficient, the displacement transfer coefficient ξ of tunnel Earthquake Inertia Force Acting
wand ξ
afor:
The angle of the axis in Φ-incident direction and tunnel;
3rd step, the longitudinal earthquake stress analysis of mountain tunnel:
V
sfor shear wave velocity (m/s), T is the predominant period (s) in place; The thickness of basement rock superstratum is H, h
1for edpth of tunnel (m); h
2for tunnel-liner overall height (m);
4th step, under shearing wave effect, the longitudinal internal force of lining cutting calculates:
Now the meridional stress of tunnel structure is combined as:
5th step, tunnel longitudinal seismic response is analyzed:
When carrying out longitudinal aseismic analysis to mountain tunnel, according to the concrete condition of mountain tunnel, carry out following analytical calculation: the 1. structure-analysis of country rock interaction to lining cutting meridional stress; 2. edpth of tunnel is to the analysis of lining cutting meridional stress; 3. basement rock shearing wave input direction is analyzed lining cutting meridional stress; 4. overlying strata is analyzed lining cutting meridional stress eigenperiod; 5. tunnel foundation bedding value K
a=K
wanalysis to lining cutting stress when=β G, β change; 6. the analysis to lining cutting stress during elastic modulus of surrounding rocks change.
The method of mountain tunnel aseismic analysis provided by the invention, unified thought is adopted to analyze mountain tunnel xsect aseismic analysis and longitudinal overall aseismic analysis first, thinking is novel, method is feasible, pass through case verification, for the mountain tunnel primary design under severe earthquake action, there is good applicability.
The present invention need not by complicated mathematical modeling or numerical analysis techniques, only need to obtain geological mapping data more accurately in the geotechnical engineering investigation stage, the antidetonation proximate analysis of mountain tunnel under severe earthquake action can be carried out by the present invention, comparative analysis is optimized to various responsive parameter, obtains corresponding analysis of data and design parameter; Obtain the application of relevant design unit, repercussion is better, has liberated the previous work amount of Tunnel Engineering designer greatly, provides necessary parameter for carrying out detailed design.
The embodiment of the present invention, compares ANSYS modeling numerical method and the self programming result of calculation of the present invention, obtains some technological merit following:
(1) mountain tunnel aseismic analysis utilizes the numerical analysis methods such as ANSYS modeling, modeling is complicated, the modeling time is probably at about 24h, and utilize the self-editing calculation procedure of the present invention, as long as geotechnical investigation data is known, the great previous work amount having liberated Tunnel Engineering designer, counting yield very fast (about 1h);
(2) xsect aseismic analysis angle, the counting yield of the inventive method is not only high, and its result of calculation and engineering present situation completely the same.Application the present invention obtains larger lining cutting moment of flexure, and axle power is taken second place; And obtaining larger axle power by Finite Element Method, moment of flexure is taken second place.And engineering experience is known, mountain tunnel is under severe earthquake action, and liner structure bears the detrusion on stratum completely, and the variation of larger shape occurs liner structure, thus producing larger bending and shearing, calculating conclusion of the present invention is consistent with the time of day of mountain tunnel; Know according to application example of the present invention, bending and shearing within the scope of tunnel vault 45 degree is all larger, be the main region that hole ejects existing drawing crack seam and shear failure, arch springing position is also the place that maximal bending moment occurs, these conclusions are also consistent with the present situation of mountain tunnel eaerthquake damage.Illustrate that the present invention not only saves time, liberation human resources greatly, and calculated mass is also higher.
(3) from mountain tunnel longitudinally overall aseismic analysis, method provided by the invention can carry out comprehensive comparative analysis in above-mentioned various aspects, carries out tunnel and optimizes aseismic analysis, not only convenient, counting yield is high, and its result is consistent with case history.And the numerical methods such as ANSYS not only longitudinal Holistic modeling difficulty, and be difficult to carry out above-mentioned analysis, so cross-sectional analysis is all carried out in the numerical analyses such as the ANSYS of current most of mountain tunnel, seldom there is the application example of longitudinal holistic approach.
Accompanying drawing explanation
Fig. 1 is the method flow diagram setting up tunnel cross sectional analyses of shake that the embodiment of the present invention provides;
Fig. 2 is the stratum deformation pattern diagram that the embodiment of the present invention provides;
Fig. 3 is the velocity response spectrum schematic diagram of the unit level seismic coefficient that the embodiment of the present invention provides;
Fig. 4 is the horizontal vibrating stress distribution schematic diagram that the embodiment of the present invention provides;
Fig. 5 is the distortion schematic diagram of tunnel-country rock system that the embodiment of the present invention provides;
Fig. 6 is the stressed schematic diagram of tunnel unit that the embodiment of the present invention provides.
Embodiment
In order to make object of the present invention, technical scheme and advantage clearly understand, below in conjunction with embodiment, the present invention is further elaborated.Should be appreciated that specific embodiment described herein only in order to explain the present invention, be not intended to limit the present invention.
Below in conjunction with drawings and the specific embodiments, application principle of the present invention is further described.
The specific embodiment of the invention is: certain mountain tunnel, and stratum is more single, and its formation parameter, Analysis of Field Geotechnical Parameters and tunnel support parameter are in table 1.
Table 1 embodiment calculating parameter
As shown in Figure 1, the method for the mountain tunnel aseismic analysis of the embodiment of the present invention comprises the following steps:
S101: the geometric parameter determining mountain tunnel, determines stratum basic parameter and simplification, determines the kinetic parameter of rock mass;
S102: according to geological exploration data and seismic strata deformable statistical data, determines the deformation pattern on stratum and the level of stratum deformation and vertical deformation;
S103: mountain tunnel xsect aseismic analysis;
S104: mountain tunnel Horizontal vibration of piles (Horizontal Deformation pressure q
h); Mountain tunnel Vertical Vibration (comprises stratum vertical deformation pressure q
v1with loose ground Earthquake Inertia Force Acting q
v2); Utilize structural mechanics, determine the internal force of tunnel-liner;
S105: mountain tunnel is overall aseismic analysis longitudinally, sets up the governing equation determination displacement transfer coefficient ξ of the longitudinal overall aseismic analysis of mountain tunnel
wand ξ
a;
S106: determine the longitudinal earthquake stress in tunnel, consider the randomness of geological process and the uncertainty of action direction, determines the meridional stress combination of tunnel structure;
S107: when carrying out longitudinal aseismic analysis to mountain tunnel, according to the concrete condition of mountain tunnel, carries out analytical calculation.
In step S101, determine that the geometric parameter of mountain tunnel comprises: h
1-edpth of tunnel (m), h
2-tunnel-liner overall height (m), the total span (m) in B-tunnel, the average thickness (m) of δ-tunnel-liner, R
0the equivalent redius (m) in-tunnel, E
sthe elastic modulus (Gpa) of-tunnel-liner, μ
sthe Poisson ratio of-tunnel-liner;
Determine that stratum basic parameter and simplification comprise: γ=∑ γ
ih
i/ H, μ
d=0.45, v
s=H/ ∑ (h
i/ v
si), G
d=γ v
s 2/ g, E
d=2 (1-μ
d) G
d, the thickness (m) of H-top layer ground, E-soil-structure interactions elastic modulus; The Poisson ratio of μ-soil;
Determine that the kinetic parameter of rock mass comprises: K
h-top layer Foundation Design horizontal seismic coefficient, K
v-top layer Foundation Design vertical seismic action coefficient; a
v-stratum maximum vertical acceleration (m/s), T-stratum natural period, V
s-elastic shear velocity of wave, S
vthe velocity response spectrum of-vibrations reference field, K
a-horizontal resiliency foundation coefficient, K
w-vertical flat elastic foundation coefficient.
In step s 102, determine that the formula of the deformation pattern on stratum and the level of stratum deformation and vertical deformation is:
In step S104, determine that the internal force of tunnel-liner is according to the mechanics of materials, tunnel cross sectional earthquake stress is:
In step s 106, determine that the meridional stress combinatorial formula of tunnel structure is:
In step s 107, according to the concrete condition of mountain tunnel, the analytical calculation carried out comprises:
1. structure-the analysis of country rock interaction to lining cutting meridional stress;
2. edpth of tunnel is to the analysis of lining cutting meridional stress;
3. basement rock shearing wave input direction is analyzed lining cutting meridional stress;
4. overlying strata is analyzed lining cutting meridional stress eigenperiod;
5. tunnel foundation bedding value K
a=K
wanalysis to lining cutting stress when=β G, β change;
6. the analysis to lining cutting stress during elastic modulus of surrounding rocks change.
The method of the mountain tunnel aseismic analysis of the embodiment of the present invention mainly contains following two aspects:
1, for tunnel cross sectional aseismic analysis aspect, based on the elastic shear deformation theory of stratum, consider the earthquake deadweight inertial force of stratum horizontal distortion, stratum vertical deformation, tunnel top fragmented rock body first, the relevant issues of research tunnel cross sectional antidetonation simultaneously;
Concrete implementation method:
In order to set up the mathematical model of tunnel cross sectional analyses of shake and be convenient to obtain approximate calculation method, special proposition some supposition following:
1. the soil body is isotropic linear elasticity continuous medium;
2. edpth of tunnel is generally more than 5 times of tunnel radius, can think that tunnel is deep tunnel, meet Fig. 4 simplified condition;
3. stratum deformation linearly changes consideration from top to bottom, and tunnel detrusion can be similar to and replace with center, tunnel stratum deformation.
4. the impact that the difference of tunnel shape causes is confined to the local areas such as corner angle more, little on the differentiation Changing Pattern impact of generally sheaf space stress field.Therefore, to section configuration other shape underground structures non-circular, adopt the conversion form of " equivalent radius ", be not only conducive to the theoretical analysis of problem, and be easy to engineering site practical application.
The first step, according to geological exploration data and seismic strata deformable statistical data, determines deformation pattern (as shown in Figure 2) and stratum deformation (level and the vertical deformation) amount on stratum, as shown in Figure 3:
In formula: S
vvelocity response spectrum (m/s) during-unit seismic coefficient;
The basic natural period (s) of T-top layer ground;
K
hdesign level seismic coefficient (earthquake degree) on the ground of-top layer;
Second step, as shown in Figure 4, horizontal vibrating stress distribution, mountain tunnel Horizontal vibration of piles:
3rd step, mountain tunnel Vertical Vibration:
A. stratum vertical deformation pressure:
B. utilize existing highway or Railway Design specification, top, hole loose ground pressure is:
In formula: a
vfor stratum maximum vertical acceleration:
4th step, under shearing wave effect, lining cutting xsect internal force calculates:
Highway and railway tunnel many employings curve wall-lining cutting of meizoseismal area, know according to assumed condition, the thought that the version of this lining cutting can apply " equivalent radius " is converted into an equivalent circular liner to carry out the Quintic system seismic response analysis in tunnel, utilizes structural mechanics and mechanics of materials knowledge to obtain internal force and the cross section stresses of tunnel lining structure:
N=q
vR sin
2θ+q
hR sin2θ
Q=-q
vR sinθcosθ-q
hR cos2θ
In formula: δ is the average thickness (m) of tunnel-liner;
5th step, utilizes reinforced concrete substantially theoretical, carries out the arrangement of reinforcement in tunnel-liner cross section.
For the said process of mountain tunnel xsect aseismic analysis design, key calculates tunnel level and vertical deformation pressure and vertical seismic action inertial force, application MATLAB programming technique, and above-mentioned thought can realize by following program:
%%%%%%%%%%%%%%%%%%%% tunnel cross sectional earthquake load analysis subroutine
2, for longitudinally overall aseismic analysis aspect, tunnel, take beam on elastic foundation as model, consider the interaction of tunnel and country rock, utilize continuum mechanics, according to anti-knock power theory, displacement transfer theory, set up tunnel structure power balance equation, this equation considers the deadweight inertial force in tunnel theoretically.As can be seen from this equation, classical reaction deflection method is the special case of the method.Tunnel meridional stress under the eigenperiod (soft or hard degree) on the thickness of comparative analysis tunnel in the size of the input direction of basement rock seismic event, earthquake shearing velocity of wave, overlying strata, stratum, the elastic stiffness coefficient impact on stratum and the Changing Pattern of strain, the change contrasting each factor, on the impact of tunnel meridional stress, finds rational parameter value scope.
Concrete implementation method:
In order to set up the mathematical model of tunnel cross sectional analyses of shake and be convenient to obtain approximate calculation method, special proposition some supposition following:
1. the soil body is isotropic linear elasticity continuous medium;
2. tunnel is considered as beam on elastic foundation, namely between tunnel with country rock for spring is connected;
3. exist between tunnel and country rock and interact, consider the relative displacement between tunnel and country rock.Tunnel produces distortion under geological process, and its distortion comprises axial deformation and transversely deforming.
4. do not consider the randomness of seismic event, time-frequency characteristic and the reflection in soil layer and scattering, seismic event is regarded as the simple harmonic wave of single-frequency.
The first step, as Fig. 5 and Fig. 6, if E is the elastic modulus of beam, I is beam section moment of inertia, and A is beam section area, K
wfor grade beam Poisson ratio, K
afor grade beam Axial distortion parameter, the transversal displacement that ν (x, t) is tunnel, the axial displacement that u (x, t) is tunnel, g
wthe transversal displacement that (x, t) is foundation soil, g
athe axial displacement that (x, t) is foundation soil; Tunnel two dimensional motion equation under ground seismic wave function:
Second step, introduces displacement transfer coefficient, considers the displacement transfer coefficient ξ of tunnel Earthquake Inertia Force Acting
wand ξ
afor:
The angle of the axis in Φ-incident direction and tunnel;
3rd step, the longitudinal earthquake stress analysis of mountain tunnel:
V
sfor shear wave velocity (m/s), T is the predominant period (s) in place; The thickness of basement rock superstratum is H, h
1for edpth of tunnel (m); h
2for tunnel-liner overall height (m);
4th step, under shearing wave effect, the longitudinal internal force of lining cutting calculates:
Consider the randomness of geological process and the uncertainty of action direction, then now the meridional stress of tunnel structure is combined as:
5th step, tunnel longitudinal seismic response is analyzed:
When carrying out longitudinal aseismic analysis to mountain tunnel, according to the concrete condition of mountain tunnel, can following analytical calculation be carried out:
1. structure-the analysis of country rock interaction to lining cutting meridional stress;
2. edpth of tunnel is to the analysis of lining cutting meridional stress;
3. basement rock shearing wave input direction is analyzed lining cutting meridional stress;
4. overlying strata is analyzed lining cutting meridional stress eigenperiod;
5. tunnel foundation bedding value K
a=K
wanalysis to lining cutting stress when=β G, β change;
6. the analysis to lining cutting stress during elastic modulus of surrounding rocks change.
For the said process of the longitudinal seismic analysis of mountain tunnel, apply the above-mentioned algorithm of MATLAB programming realization, program design is as follows:
Xsect Aseismic Analytical Method and longitudinal overall Aseismic Analytical Method adopt unified thought to analyze by the present invention first, and thought is novel, and scheme is reasonable.
For tunnel cross sectional aseismic analysis aspect, based on the elastic shear deformation theory of stratum, consider the earthquake deadweight inertial force of stratum horizontal distortion, stratum vertical deformation, tunnel top fragmented rock body first, the relevant issues of research tunnel cross sectional antidetonation simultaneously;
For longitudinally overall aseismic analysis aspect, tunnel, take beam on elastic foundation as model, consider the interaction of tunnel and country rock, utilize continuum mechanics, according to anti-knock power theory, displacement transfer theory, set up tunnel structure power balance equation, this equation considers the deadweight inertial force in tunnel theoretically.Classical reaction deflection method is the special case of the method.The method can analyze the Changing Pattern of tunnel meridional stress under the elastic stiffness coefficient impact on thickness in the size of the input direction of basement rock seismic event, earthquake shearing velocity of wave, overlying strata, the eigenperiod (soft or hard degree) on stratum, stratum and strain, the change contrasting each factor, on the impact of tunnel meridional stress, finds rational parameter value scope.
The foregoing is only preferred embodiment of the present invention, not in order to limit the present invention, all any amendments done within the spirit and principles in the present invention, equivalent replacement and improvement etc., all should be included within protection scope of the present invention.
Claims (8)
1. a method for mountain tunnel aseismic analysis, is characterized in that, the method for this mountain tunnel aseismic analysis comprises:
Step one, determines the geometric parameter of mountain tunnel, determines and simplifies stratum basic parameter, determining the kinetic parameter of rock mass;
Step 2, according to geological exploration data and seismic strata deformable statistical data, determines the deformation pattern on stratum and the level of stratum deformation and vertical deformation;
Step 3, mountain tunnel xsect aseismic analysis;
Step 4, mountain tunnel Horizontal vibration of piles; Mountain tunnel Vertical Vibration; Utilize structural mechanics, determine the internal force of tunnel-liner;
Step 5, mountain tunnel is overall aseismic analysis longitudinally, sets up the governing equation determination displacement transfer coefficient ξ of the longitudinal overall aseismic analysis of mountain tunnel
wand ξ
a;
Step 6, determines the longitudinal earthquake stress in tunnel, considers the randomness of geological process and the uncertainty of action direction, and then determines the meridional stress combination of tunnel structure;
Step 7, when carrying out longitudinal aseismic analysis to mountain tunnel, according to the concrete condition of mountain tunnel, carries out analytical calculation.
2. the method for mountain tunnel aseismic analysis as claimed in claim 1, is characterized in that, in step one, determine that the geometric parameter of mountain tunnel comprises: h
1-edpth of tunnel (m), h
2-tunnel-liner overall height (m), the total span (m) in B-tunnel, the average thickness (m) of δ-tunnel-liner, R
0the equivalent redius (m) in-tunnel, E
sthe elastic modulus (Gpa) of-tunnel-liner, μ
sthe Poisson ratio of-tunnel-liner;
Determine and simplify stratum basic parameter to comprise: γ=∑ γ
ih
i/ H, μ
d=0.45, v
s=H/ ∑ (h
i/ v
si), G
d=γ v
s 2/ g, E
d=2 (1-μ
d) G
d, the thickness (m) of H-top layer ground, E-soil-structure interactions elastic modulus; The Poisson ratio of μ-soil;
Determine that the kinetic parameter of rock mass comprises: K
h-top layer Foundation Design horizontal seismic coefficient, K
v-top layer Foundation Design vertical seismic action coefficient; a
v-stratum maximum vertical acceleration (m/s
2), T-stratum natural period, V
s-elastic shear velocity of wave, S
vthe velocity response spectrum of-vibrations reference field, K
a-horizontal resiliency foundation coefficient, K
w-vertical flat elastic foundation coefficient.
3. the method for mountain tunnel aseismic analysis as claimed in claim 1, is characterized in that, in step 2, determine that the formula of the deformation pattern on stratum and the level of stratum deformation and vertical deformation is:
4. the method for mountain tunnel aseismic analysis as claimed in claim 1, it is characterized in that, in step 4, determine that the internal force of tunnel-liner is according to the mechanics of materials, tunnel cross sectional earthquake stress is:
5. the method for mountain tunnel aseismic analysis as claimed in claim 1, is characterized in that, in step 6, determine that the meridional stress combinatorial formula of tunnel structure is:
6. the method for mountain tunnel aseismic analysis as claimed in claim 1, it is characterized in that, in step 7, according to the concrete condition of mountain tunnel, the analytical calculation carried out comprises:
1. structure-the analysis of country rock interaction to lining cutting meridional stress;
2. edpth of tunnel is to the analysis of lining cutting meridional stress;
3. basement rock shearing wave input direction is analyzed lining cutting meridional stress;
4. overlying strata is analyzed lining cutting meridional stress eigenperiod;
5. tunnel foundation bedding value K
a=K
wanalysis to lining cutting stress when=β G, β change;
6. the analysis to lining cutting stress during elastic modulus of surrounding rocks change.
7. the method for mountain tunnel aseismic analysis as claimed in claim 1, it is characterized in that, the method for this mountain tunnel aseismic analysis for the implementation method that tunnel cross sectional aseismic analysis is concrete is:
The first step, according to geological exploration data and seismic strata deformable statistical data, determine the deformation pattern on stratum and stratum deformation level and vertical deformation:
In formula: S
vvelocity response spectrum (m/s) during-unit seismic coefficient;
The basic natural period (s) of T-top layer ground;
K
hdesign level seismic coefficient (earthquake degree) on the ground of-top layer;
Second step, mountain tunnel Horizontal vibration of piles:
3rd step, mountain tunnel Vertical Vibration:
A. stratum vertical deformation pressure:
B. utilize existing highway or Railway Design specification, top, hole loose ground pressure is:
In formula: a
vfor stratum maximum vertical acceleration:
4th step, under shearing wave effect, lining cutting xsect internal force calculates:
Highway and railway tunnel many employings curve wall-lining cutting of meizoseismal area, according to waiting generation circle thought, and utilizes structural mechanics and mechanics of materials knowledge to obtain internal force and the cross section stresses of tunnel lining structure:
N=q
vRsin
2θ+q
hRsin2θ
Q=-q
vRsinθcosθ-q
hRcos2θ
In formula: δ is the average thickness (m) of tunnel-liner;
5th step, utilizes reinforced concrete substantially theoretical, carries out the arrangement of reinforcement in tunnel-liner cross section.
8. the method for mountain tunnel aseismic analysis as claimed in claim 1, is characterized in that, the method for this mountain tunnel aseismic analysis for the tunnel implementation method that longitudinally overall aseismic analysis is concrete is:
The first step, determine the two dimensional equation of longitudinal body vibration of mountain tunnel:
E is the elastic modulus of beam, and I is beam section moment of inertia, and A is beam section area, K
wfor grade beam Poisson ratio, K
afor grade beam Axial distortion parameter, the transversal displacement that ν (x, t) is tunnel, the axial displacement that u (x, t) is tunnel, g
wthe transversal displacement that (x, t) is foundation soil, g
athe axial displacement that (x, t) is foundation soil;
Second step, introduces displacement transfer coefficient, the displacement transfer coefficient ξ of tunnel Earthquake Inertia Force Acting
wand ξ
afor:
The angle of the axis in Φ-incident direction and tunnel;
3rd step, the longitudinal earthquake stress analysis of mountain tunnel:
V
sfor shear wave velocity (m/s), T is the predominant period (s) in place; The thickness of basement rock superstratum is H, h
1for edpth of tunnel (m); h
2for tunnel-liner overall height (m);
4th step, under shearing wave effect, the longitudinal internal force of lining cutting calculates:
Now the meridional stress of tunnel structure is combined as:
5th step, tunnel longitudinal seismic response is analyzed:
When carrying out longitudinal aseismic analysis to mountain tunnel, according to the concrete condition of mountain tunnel, carry out following analytical calculation: the 1. structure-analysis of country rock interaction to lining cutting meridional stress; 2. edpth of tunnel is to the analysis of lining cutting meridional stress; 3. basement rock shearing wave input direction is analyzed lining cutting meridional stress; 4. overlying strata is analyzed lining cutting meridional stress eigenperiod; 5. tunnel foundation bedding value K
a=K
wanalysis to lining cutting stress when=β G, β change; 6. the analysis to lining cutting stress during elastic modulus of surrounding rocks change.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510364671.0A CN104992013A (en) | 2015-06-26 | 2015-06-26 | Mountain tunnel anti-seismic analysis method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510364671.0A CN104992013A (en) | 2015-06-26 | 2015-06-26 | Mountain tunnel anti-seismic analysis method |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104992013A true CN104992013A (en) | 2015-10-21 |
Family
ID=54303827
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510364671.0A Pending CN104992013A (en) | 2015-06-26 | 2015-06-26 | Mountain tunnel anti-seismic analysis method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104992013A (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107122536A (en) * | 2017-04-20 | 2017-09-01 | 济南轨道交通集团有限公司 | A kind of running tunnel antidetonation method for numerical simulation |
CN109740251A (en) * | 2018-12-30 | 2019-05-10 | 中铁十四局集团有限公司 | Selection method, device, memory and the processor of the parameter of secondary lining |
CN111022119A (en) * | 2019-11-18 | 2020-04-17 | 东莞理工学院 | Shield tunnel structure anti-seismic safety evaluation method based on pattern recognition |
CN111339602A (en) * | 2020-04-09 | 2020-06-26 | 中国人民解放***箭军研究院核技术研究所 | Method for determining earthquake resistance evaluation index of underground vertical wall vault tunnel structure |
CN113358466A (en) * | 2021-05-28 | 2021-09-07 | 燕山大学 | Method and system for determining transfer coefficient of dynamic stress of layered foundation soil layer interface |
CN116611132A (en) * | 2023-03-13 | 2023-08-18 | 西南交通大学 | Method and system for calculating longitudinal internal force of segmental lining based on stratum deformation |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103310039A (en) * | 2013-05-22 | 2013-09-18 | 青岛理工大学 | Method for detecting influence of underground construction on ground building structure |
CN104500093A (en) * | 2014-12-16 | 2015-04-08 | 湖南科技大学 | Method for positioning temporary middle walls in large-section tunnel supporting structure of formation with soft top and hard bottom |
CN104537215A (en) * | 2014-12-16 | 2015-04-22 | 上海交通大学 | Method for determining longitudinal internal force of shield tunnel under load effect |
-
2015
- 2015-06-26 CN CN201510364671.0A patent/CN104992013A/en active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103310039A (en) * | 2013-05-22 | 2013-09-18 | 青岛理工大学 | Method for detecting influence of underground construction on ground building structure |
CN104500093A (en) * | 2014-12-16 | 2015-04-08 | 湖南科技大学 | Method for positioning temporary middle walls in large-section tunnel supporting structure of formation with soft top and hard bottom |
CN104537215A (en) * | 2014-12-16 | 2015-04-22 | 上海交通大学 | Method for determining longitudinal internal force of shield tunnel under load effect |
Non-Patent Citations (2)
Title |
---|
张斌伟 等: "山岭隧道横截面抗震分析的近似方法", 《公路交通科技》 * |
张斌伟 等: "隧道纵向整体抗震分析的近似方法", 《岩土力学》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107122536A (en) * | 2017-04-20 | 2017-09-01 | 济南轨道交通集团有限公司 | A kind of running tunnel antidetonation method for numerical simulation |
CN109740251A (en) * | 2018-12-30 | 2019-05-10 | 中铁十四局集团有限公司 | Selection method, device, memory and the processor of the parameter of secondary lining |
CN111022119A (en) * | 2019-11-18 | 2020-04-17 | 东莞理工学院 | Shield tunnel structure anti-seismic safety evaluation method based on pattern recognition |
CN111339602A (en) * | 2020-04-09 | 2020-06-26 | 中国人民解放***箭军研究院核技术研究所 | Method for determining earthquake resistance evaluation index of underground vertical wall vault tunnel structure |
CN113358466A (en) * | 2021-05-28 | 2021-09-07 | 燕山大学 | Method and system for determining transfer coefficient of dynamic stress of layered foundation soil layer interface |
CN113358466B (en) * | 2021-05-28 | 2022-10-04 | 燕山大学 | Method and system for determining transfer coefficient of dynamic stress of layered foundation soil layer interface |
CN116611132A (en) * | 2023-03-13 | 2023-08-18 | 西南交通大学 | Method and system for calculating longitudinal internal force of segmental lining based on stratum deformation |
CN116611132B (en) * | 2023-03-13 | 2024-02-27 | 西南交通大学 | Method and system for calculating longitudinal internal force of segmental lining based on stratum deformation |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104992013A (en) | Mountain tunnel anti-seismic analysis method | |
MacLaughlin et al. | Review of validation of the discontinuous deformation analysis (DDA) method | |
Tang | Applications of rock failure process analysis (RFPA) method | |
Chen et al. | Seismic performance upgrading for underground structures by introducing shear panel dampers | |
Mackie et al. | Performance-based earthquake assessment of bridge systems including ground-foundation interaction | |
CN102654054B (en) | Method for determining height of caving zone of stope with internal-staggered layer-staggered position type roadway layout | |
CN105041331B (en) | A kind of tunnel portal section anti-seismic fortified length computational methods | |
Barla et al. | Distinct element modelling in static and dynamic conditions with application to an underground archaeological site | |
Mayoral et al. | Seismic response of bridges with massive foundations | |
El-Attar | Dynamic analysis of combined piled raft system (CPRS) | |
Akhtarpour et al. | Dynamic response of a tall building next to deep excavation considering soil–structure interaction | |
Mayoral et al. | Site response effects on an urban overpass | |
Chen et al. | Seismic performance of tunnel lining of side-by-side and vertically stacked twin-tunnels | |
Lu et al. | Numerical and experimental analyses for bearing capacity of rigid strip footing subjected to eccentric load | |
Li | Failure of an under-dip shale slope and its response under excavation conditions | |
Kamijo et al. | Seismic tests of a pile-supported structure in liquefiable sand using large-scale blast excitation | |
Sun et al. | Stability of bolt-supported concealed bedding rock slopes with respect to bi-planar failure | |
Theinat | 3D numerical modelling of micropiles interaction with soil & rock | |
Kong et al. | Stability evaluation of flat large-span cavern in jointed rock mass | |
Lebedev et al. | Face stability in heavy clay: Theory and practice | |
Liu et al. | High-Stress Chamber Stability in Kilometer-Deep Coal Mines: A Case Study | |
Jiao et al. | Numerical analysis of building structures response under effect of continuous and stepped non-uniform settlement | |
Nghiem et al. | Soil-structure interaction effects of high rise buildings | |
Danjiang et al. | Calculation of Rock Pressure in Loess Tunnels Based on Limit Equilibrium Theory and Analysis of Influencing Factors | |
Xie et al. | Time-varying response analysis of surrounding rock construction mechanics of a double-arch tunnel without a middle pilot tunnel and suggestions for tunnel construction |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20151021 |