CN104331607A - Slope stability analyzing method - Google Patents

Slope stability analyzing method Download PDF

Info

Publication number
CN104331607A
CN104331607A CN201410564243.8A CN201410564243A CN104331607A CN 104331607 A CN104331607 A CN 104331607A CN 201410564243 A CN201410564243 A CN 201410564243A CN 104331607 A CN104331607 A CN 104331607A
Authority
CN
China
Prior art keywords
slope
radius
circle
center
stick
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
Application number
CN201410564243.8A
Other languages
Chinese (zh)
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.)
Xian Changqing Technology Engineering Co Ltd
Original Assignee
Xian Changqing Technology Engineering Co Ltd
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 Xian Changqing Technology Engineering Co Ltd filed Critical Xian Changqing Technology Engineering Co Ltd
Priority to CN201410564243.8A priority Critical patent/CN104331607A/en
Publication of CN104331607A publication Critical patent/CN104331607A/en
Pending legal-status Critical Current

Links

Landscapes

  • Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)

Abstract

The invention provides a slope stability analyzing method, which comprises the following steps that 1, slope parameters are input; 2, the circle center O and the radius R of a slope are assumed; 3, a crossed part of a half slope and a circular arc drawn according to the assumed circle center O and radius R is subjected to strip splitting and numbering; 4, Li, Wi and theta<i> values of each strip are calculated; 5, the safety factor Fs is obtained according to the Li, Wi and theta<i> values; 6, a plurality of Fs are obtained; 7, the minimum values of the safety factor Fs is screened and is used as the slope safety factor to be output; according to the slope length and the slope height of the half slope, the circle center and the radius of a slope landslide are assumed, in addition, a circular arc is drawn according to the assumed circle center and radius, a crossed part of the circular arc and the slope is subjected to strip splitting, in addition, the slope safety factor is calculated, meanwhile, the circle center and the radius are assumed for many times, the slope safety factor is calculated for many times, in addition, the calculated safety factors are compared, and the minimum value is selected as the slope safety factor. The scheme has the advantage that the calculation is simple and accurate.

Description

A kind of slope stability analysis method
Technical field
The present invention relates to a kind of field of civil engineering, particularly a kind of slope stability analysis method.
Background technology
The loess plateau ravines and guillies criss-cross of western part of China, when annual rainy season set in 7 ~ September, these side slopes very easily come down, and cause casualties and economic loss.Therefore, in engineering construction addressing, the stability of correct evaluation side slope is particularly important.At present, loess slope stability field evaluation mainly relies on the experience of slip-stick artist, and experience has uncertainty, is limited by the shortcoming of experience; Although computer numerical analytical technology principle is ripe, it is accurate to calculate, because analytic process is loaded down with trivial details, computing machine is not suitable for running under the field environment in zanjon deserted mountain, so there is no be adopted by engineering.
Summary of the invention
In order to solve the problem, the invention provides and a kind ofly calculate the slope stability analysis method that simply can analyze the safety coefficient of side slope accurately.
For achieving the above object, by the following technical solutions, a kind of slope stability analysis method, comprises the following steps in the present invention:
1) slope parameter is inputted: slope height H, length of grade L, angle of friction φ, cohesive strength C, severe γ, saturated unit weight γ sat, point number n, seepage flow volume ratio S;
2) center of circle O and the radius R of side slope is supposed;
3) part that the circular arc drawn according to center of circle O and the radius R of hypothesis is crossing with half way up the mountain is carried out itemize and is numbered;
4) L of each bar of following formulae discovery is used i, W i, θ ivalue;
L i=R* (θ ii-1)---the length of arc bottom every stick;
θ i=arcsin [(O x-X i)/R]---mid point D's at the bottom of i-th rule
θ 0=arcsin [(O x-L)/R]---the inclination angle of mid point D at the bottom of the 0th rule;
W i=γ * B*H i---the gravity of i-th stick;
H i=Y i'-Y i---i-th stick is high
Y i=O y-R*cos (θ i)---the Y-coordinate of i-th stick top mid point D '
Y i'=(-H/L*X i)---the Y-coordinate of i-th stick top mid point D '
O xthe X-coordinate of-center of circle O, O ythe Y-coordinate of-center of circle O, X ithe horizontal ordinate of mid point D at the bottom of-the i-th rule, B-bar is wide;
5) according to the L that above formula calculates i, W i, θ i, bring equation into and can draw safety coefficient Fs;
Fs={γ*B*tanφ*∑(H i*cosθ i)+∑(c*L i)}/{γ*B*∑(H i*sinθ i)};
6) repeatedly suppose center of circle O and radius R according to grid is infinite, and repeat step 3), 4), 5) draw multiple Fs;
7) minimum value of screening safety coefficient Fs exports as the safety coefficient of side slope.
Step 4) L that calculates i, W i, θ i, following formula can also be used can to calculate infiltration safety coefficient Fs ';
Fs′={γ*B*tanφ*∑(H i1*cosθ i)+(γ sat-10)*B*tanφ*∑(H i2*cosθ i)+∑(c*L i)}/{γ*B*∑(H i1*sinθ i)+γ sat*B*∑(H i2*sinθ i)}
H in formula i1=H i* (1-S), H i2=H i* S, H i1article-the i-th, air draght, H i2article-the i-th, under water highly.
The present invention uses technique scheme, have the following advantages, the present invention supposes the center of circle and the radius of slope and land slide according to the length of grade of half way up the mountain and slope height, and drawing circular arc according to the center of circle supposed and radius, the part crossing with side slope to this circular arc carries out itemize, and calculates the safety coefficient of side slope, the simultaneously hypothesis center of circle repeatedly and radius, and carrying out repeatedly calculates the safety coefficient of side slope, and the safety coefficient calculated is compared, select the safety coefficient of minimum value as side slope; This scheme calculates simple, accurately.
Below with reference to accompanying drawing, the present invention is described in further details.
Accompanying drawing explanation
Fig. 1 is a kind of slope stability analysis method flow diagram;
Fig. 2 is side slope itemize figure;
Fig. 3 is the force analysis figure of every bar;
Fig. 4 is slope stability analysis figure.
Embodiment
Further a kind of slope stability analysis method is described in detail below in conjunction with drawings and Examples.
A kind of slope stability analysis method as shown in Figure 1, comprises the following steps:
1) slope parameter is inputted: slope height H, length of grade L, angle of friction φ, cohesive strength C, severe γ, saturated unit weight γ sat, point number n, seepage flow volume ratio S;
2) center of circle O and the radius R of side slope is supposed;
3) part that the circular arc drawn according to center of circle O and the radius R of hypothesis is crossing with half way up the mountain is carried out itemize and is numbered;
4) L of each bar of following formulae discovery is used i, W i, θ ivalue;
L i=R* (θ ii-1)---the length of arc bottom every stick;
θ i=arcsin [(O x-X i)/R]---the inclination angle of mid point D at the bottom of i-th rule
θ 0=arcsin [(O x-L)/R]---the inclination angle of mid point D at the bottom of the 0th rule;
W i=γ * B*H i---the gravity of i-th stick;
H i=Y i'-Y i---i-th stick is high
Y i=O y-R*cos (θ i)---the Y-coordinate of i-th stick top mid point D '
Y i'=(-H/L*X i)---the Y-coordinate of i-th stick top mid point D '
O xthe X-coordinate of-center of circle O, O ythe Y-coordinate of-center of circle O, X ithe horizontal ordinate of mid point D at the bottom of-the i-th rule, B-bar is wide;
5) according to the L that above formula calculates i, W i, θ i, bring equation into and can draw safety coefficient Fs;
Fs={γ*B*tanφ*∑(H i*cosθ i)+∑(c*L i)}/{γ*B*∑(H i*sinθ i)}
6) repeatedly suppose center of circle O and radius R according to grid is infinite, and repeat step 3), 4), 5) draw multiple Fs;
7) minimum value of screening safety coefficient Fs exports as the safety coefficient of side slope.
Step 2) in the center of circle O of hypothesis and radius R be suppose in the quadrilateral of the AGFE identified in such as Fig. 4, this quadrilateral is according to slip-stick artist, the region that staff is summed up in routine duties;
In computation process, the parameter understood is needed to can refer to following table:
Table one: variable and formulae discovery parameter lookup table:
The present invention adopts this technology, and compared with the result that the result obtained is positive with reason under error, accuracy is high, and calculates simple, and this result of calculation is easily programmed into the mobile devices such as mobile phone, does not need in a large amount of equipment calculated;
Table two: be this technology checkout result and the positive result of calculation control reference table one of reason;
Table three: be this technology checkout result and the positive result of calculation control reference table two of reason;
Table four: be this technology checkout result and the positive result of calculation control reference table three of reason;
Can find out according to upper table two, three, four, the technical solution used in the present invention, the result calculated is less with the error of calculation of the positive result of reason, is suitable for and analyzes the safety coefficient of side slope.
1) slope parameter is inputted: slope height H, length of grade L, angle of friction φ, cohesive strength C, severe γ, saturated unit weight γ sat, point number n, seepage flow volume ratio S;
2) center of circle O and the radius R of side slope is supposed;
3) part that the circular arc drawn according to center of circle O and the radius R of hypothesis is crossing with half way up the mountain is carried out itemize and is numbered;
4) L of each bar of following formulae discovery is used i, W i, θ ivalue;
L i=R* (θ ii-1)---the length of arc bottom every stick;
θ i=arcsin [(O x-X i)/R]---the θ of mid point D at the bottom of i-th rule 0=arcsin [(O x-L)/R]---the inclination angle of mid point D at the bottom of the 0th rule;
W i=γ * B*H i---the gravity of i-th stick;
H i=Y i'-Y i---i-th stick is high
Y i=O y-R*cos (θ i)---the Y-coordinate of i-th stick top mid point D '
Y i'=(-H/L*X i)---the Y-coordinate of i-th stick top mid point D '
O xthe X-coordinate of-center of circle O, O ythe Y-coordinate of-center of circle O, X ithe horizontal ordinate of mid point D at the bottom of-the i-th rule, B-bar is wide;
5) according to step 4) L that calculates i, W i, θ i, following formula can also be used can to calculate infiltration safety coefficient Fs ';
Fs′={γ*B*tanφ*∑(H i1*cosθ i)+(γ sat-10)*B*tanφ*∑(H i2*cosθ i)+∑(c*L i)}/{γ*B*∑(H i1*sinθ i)+γ sat*B*∑(H i2*sinθ i)}
H in formula i1=H i* (1-S), H i2=H i* S, H i1article-the i-th, air draght, H i2article-the i-th, under water highly.
6) repeatedly suppose center of circle O and radius R according to grid is infinite, and repeat step 3), 4), 5) draw multiple Fs ';
7) minimum value of screening safety coefficient Fs ' exports as the safety coefficient of side slope.
The present invention supposes the center of circle and the radius of slope and land slide according to the length of grade of half way up the mountain and slope height, and draw circular arc according to the center of circle supposed and radius, the part crossing with side slope to this circular arc carries out itemize, and calculate the safety coefficient of side slope, the simultaneously hypothesis center of circle repeatedly and radius, and carrying out repeatedly calculates the safety coefficient of side slope, and the safety coefficient calculated is compared, select the safety coefficient of minimum value as side slope; This scheme calculates simple, accurately.
More than exemplifying is only illustrate of the present invention, does not form the restriction to protection scope of the present invention, everyly all belongs within protection scope of the present invention with the same or analogous design of the present invention.

Claims (2)

1. a slope stability analysis method, is characterized in that, comprises the following steps:
1) slope parameter is inputted: slope height H, length of grade L, angle of friction φ, cohesive strength C, severe γ, saturated unit weight γ sat, point number n, seepage flow volume ratio S;
2) center of circle O and the radius R of side slope is supposed;
3) part that the circular arc drawn according to center of circle O and the radius R of hypothesis is crossing with half way up the mountain is carried out itemize and is numbered;
4) L of each bar of following formulae discovery is used i, W i, θ ivalue;
L i=R* (θ ii-1)---the length of arc bottom every stick;
θ i=arcsin [(O x-X i)/R]---the inclination angle of mid point D at the bottom of i-th rule;
θ 0=arcsin [(O x-L)/R]---the inclination angle of mid point D at the bottom of the 0th rule;
W i=γ * B*H i---the gravity of i-th stick;
H i=Y i'-Y i---i-th stick is high;
Y i=O y-R*cos (θ i)---the Y-coordinate of i-th stick top mid point D ';
Y i'=(-H/L*X i)---the Y-coordinate of i-th stick top mid point D ';
O xthe X-coordinate of-center of circle O, O ythe Y-coordinate of-center of circle O, X ithe horizontal ordinate of mid point D at the bottom of-the i-th rule, B-bar is wide;
5) according to the L that above formula calculates i, W i, θ i, bring equation into and can draw safety coefficient Fs;
Fs={γ*B*tanφ*∑(H i*cosθ i)+∑(c*L i)}/{γ*B*∑(H i*sinθ i)}
6) repeatedly suppose center of circle O and radius R according to grid is infinite, and repeat step 3), 4), 5) draw multiple Fs;
7) minimum value of screening safety coefficient Fs exports as the safety coefficient of side slope.
2. a kind of slope stability analysis method according to claim 1, is characterized in that, step 4) L that calculates i, W i, θ i, infiltration safety coefficient Fs ' can also be calculated;
Fs′=
{γ*B*tanφ*∑(H i1*cosθ i)+(γ sat-10)*B*tanφ*∑(H i2*cosθ i)+∑(c*L i)}/{γ*B*∑(H i1*sinθ i)+γ sat*B*∑(H i2*sinθ i)}
H in formula i1=H i* (1-S), H i2=H i* S, H i1article-the i-th, air draght, H i2article-the i-th, under water highly.
CN201410564243.8A 2014-10-22 2014-10-22 Slope stability analyzing method Pending CN104331607A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410564243.8A CN104331607A (en) 2014-10-22 2014-10-22 Slope stability analyzing method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410564243.8A CN104331607A (en) 2014-10-22 2014-10-22 Slope stability analyzing method

Publications (1)

Publication Number Publication Date
CN104331607A true CN104331607A (en) 2015-02-04

Family

ID=52406330

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410564243.8A Pending CN104331607A (en) 2014-10-22 2014-10-22 Slope stability analyzing method

Country Status (1)

Country Link
CN (1) CN104331607A (en)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104988933A (en) * 2015-05-22 2015-10-21 铁道第三勘察设计院集团有限公司 Analyzing and determining method for bridge foundation slope stable slope-angle line
CN105160165A (en) * 2015-08-25 2015-12-16 成都理工大学 Strain softening side slope dynamic stability assessing method
CN105975782A (en) * 2016-05-11 2016-09-28 辽宁工程技术大学 Composite soft fluctuating base waste dump slope form optimizing method
CN108316320A (en) * 2018-01-31 2018-07-24 青岛理工大学 A kind of instable slope Design Method of Reinforcing
CN109001021A (en) * 2018-07-13 2018-12-14 中国水利水电科学研究院 Using the method for centrifugal model test evaluation geotechnical slope safety stock
CN109977554A (en) * 2019-03-28 2019-07-05 青岛理工大学 A kind of appraisal procedure of slope sliding area
CN110197047A (en) * 2019-06-18 2019-09-03 贵州正业工程技术投资有限公司 Convex annular Slope Stability Evaluation method based on Janbu method
CN111611639A (en) * 2020-05-20 2020-09-01 广东安元矿业勘察设计有限公司 Simple method for analyzing stability of muck field slope by arc method
CN111931272A (en) * 2020-07-16 2020-11-13 中南大学 Equal-precision safety coefficient calculation method for slope of any scale and grid division method

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103266588A (en) * 2013-04-28 2013-08-28 青岛理工大学 Method for determining benching non-circular slip surfaces by displacement monitoring
CN103485353A (en) * 2013-09-24 2014-01-01 昆明理工大学 Slope stability analysis slice method based on global optimization
CN103729521A (en) * 2014-01-20 2014-04-16 湖北工业大学 Slide face boundary method for calculating slope stability

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103266588A (en) * 2013-04-28 2013-08-28 青岛理工大学 Method for determining benching non-circular slip surfaces by displacement monitoring
CN103485353A (en) * 2013-09-24 2014-01-01 昆明理工大学 Slope stability analysis slice method based on global optimization
CN103729521A (en) * 2014-01-20 2014-04-16 湖北工业大学 Slide face boundary method for calculating slope stability

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
周叔举: "边坡稳定性计算中水压力的考虑", 《勘察科学技术》 *
赵胜杰 等: "边坡稳定性分析方法现状及发展简述", 《土工基础》 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104988933A (en) * 2015-05-22 2015-10-21 铁道第三勘察设计院集团有限公司 Analyzing and determining method for bridge foundation slope stable slope-angle line
CN105160165A (en) * 2015-08-25 2015-12-16 成都理工大学 Strain softening side slope dynamic stability assessing method
CN105160165B (en) * 2015-08-25 2017-11-21 成都理工大学 A kind of strain softening type slope dynamic stability evaluation method
CN105975782A (en) * 2016-05-11 2016-09-28 辽宁工程技术大学 Composite soft fluctuating base waste dump slope form optimizing method
CN108316320A (en) * 2018-01-31 2018-07-24 青岛理工大学 A kind of instable slope Design Method of Reinforcing
CN109001021A (en) * 2018-07-13 2018-12-14 中国水利水电科学研究院 Using the method for centrifugal model test evaluation geotechnical slope safety stock
CN109977554A (en) * 2019-03-28 2019-07-05 青岛理工大学 A kind of appraisal procedure of slope sliding area
CN109977554B (en) * 2019-03-28 2020-06-30 青岛理工大学 Method for evaluating sliding area of side slope
AU2019438357B2 (en) * 2019-03-28 2021-12-09 Qingdao university of technology Slope sliding area assessment method
CN110197047A (en) * 2019-06-18 2019-09-03 贵州正业工程技术投资有限公司 Convex annular Slope Stability Evaluation method based on Janbu method
CN111611639A (en) * 2020-05-20 2020-09-01 广东安元矿业勘察设计有限公司 Simple method for analyzing stability of muck field slope by arc method
CN111931272A (en) * 2020-07-16 2020-11-13 中南大学 Equal-precision safety coefficient calculation method for slope of any scale and grid division method

Similar Documents

Publication Publication Date Title
CN104331607A (en) Slope stability analyzing method
CN105513086B (en) A kind of Educational toy external member and its localization method is detected based on the matched magic square of shape
CN204461463U (en) A kind of optical three-dimensional scanning instrument
CN104794315A (en) Coal seam floor damaged zone depth determination method based on improved rupture mechanical model
CN103256906B (en) A kind of precise calculation method for electronic theodolite stumpage timber volume
CN105043958A (en) Test method for microporous membrane porosity
CN208134877U (en) A kind of wood composite pallet with limit function
CN203772722U (en) Liquid-plastic limit combined tester
CN104021312B (en) Flood process similarity analytic method and system based on process linear centroid
CN105952368B (en) A kind of beam is away from adjustable ladder
CN205561716U (en) Thickness detection device of panel
CN205940381U (en) Resonance post modularization detection device
CN204128466U (en) A kind of measurer for thickness
CN104101409B (en) Online detection method for ore storage amount of thickener
CN207717923U (en) Multi-functional detection fixture
CN206459588U (en) Gramineous crop multi-parameter two-dimentional measuring device of short stem
CN204964210U (en) Adjustable ladder climbing capability test device
CN205383995U (en) Cloud technology based measuring system
CN204346474U (en) A kind of verification force checking device
CN104951452B (en) A kind of similar power spectrum finding method in peak position suitable for neutron activation elementary analysis
CN204302273U (en) A kind of device of novel on-line monitoring urban water resource ground contamination
CN204478956U (en) A kind of core recovery slipstick
CN103886211B (en) A kind of weaponry individual event ability pitch analysis method
Cheraghi et al. Diversity of herbaceous and woody plant species in relation to different physiographic conditions using numerical and parametrical indices in zagros mountainous forests
He et al. Coal-gangue acoustic signal recognition based on sparse representation

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20150204