CN113449429A - Slope stability evaluation and correction method based on local average - Google Patents

Slope stability evaluation and correction method based on local average Download PDF

Info

Publication number
CN113449429A
CN113449429A CN202110777072.7A CN202110777072A CN113449429A CN 113449429 A CN113449429 A CN 113449429A CN 202110777072 A CN202110777072 A CN 202110777072A CN 113449429 A CN113449429 A CN 113449429A
Authority
CN
China
Prior art keywords
slope
variable
soil
original
layer
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
CN202110777072.7A
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.)
PowerChina Guiyang Engineering Corp Ltd
Original Assignee
PowerChina Guiyang Engineering Corp 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 PowerChina Guiyang Engineering Corp Ltd filed Critical PowerChina Guiyang Engineering Corp Ltd
Priority to CN202110777072.7A priority Critical patent/CN113449429A/en
Publication of CN113449429A publication Critical patent/CN113449429A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/15Correlation function computation including computation of convolution operations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Databases & Information Systems (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)

Abstract

The invention discloses a slope stability evaluation and correction method based on local average, which comprises the following steps: A. setting initial soil characteristic parameters to carry out Monte Carlo simulation, generating random samples, judging whether the slope is invalid or not through each group of samples, counting the number of invalid samples, and calculating the probability estimation value of the original slope invalid
Figure DDA0003155931300000011
B. Statistical failure sample x under original Monte Carlo simulationjJ 1, 2.... am; C. calculating the original and corrected combined probability density function according to the original and corrected soil characteristic parameters to obtain the weight index omegajJ 1, 2.... am; D. calculating the corrected failure probability and the estimated value of the variation coefficient thereof. The method solves the problems that the Monte Carlo simulation method which is a widely applied reliability analysis method needs to rebuild the model after considering the fluctuation range or other soil characteristic parameters, and particularly for the slope with small probability of failure, a large amount of calculation consumes huge time, energy and computer resources.

Description

Slope stability evaluation and correction method based on local average
Technical Field
The invention relates to a slope stability evaluation and correction method based on local average, and belongs to the technical field of geotechnical engineering.
Background
The landslide problem caused by slope instability is a geological disaster with extremely strong destructive power, however, many uncertain factors exist in slope engineering, such as uncertainty of characteristic parameters of rock and soil mass in fig. 2, uncertainty of a calculation model and the like. The traditional slope stability analysis ignores the uncertainties, the safety degree of the slope cannot be truly reflected, and the slope stability analysis based on the reliability theory can consider the factors. Slope stability analysis methods based on the reliability theory are increasingly researched in recent years, and the design specifications of the slope of the hydraulic and hydroelectric engineering clearly suggest that the reliability analysis should be carried out when the grade 1 slope is conditioned.
For the problem of stability and reliability of the side slope, statistical parameters of soil features, such as mean, variance, distribution type and the like, are determined by a mathematical statistical method according to prior information such as field or indoor test data in preliminary construction. And obtaining an estimated value of the slope failure probability according to the prior information, wherein the estimated value is only a preliminary estimated value of the slope failure probability. With the deep progress of slope engineering construction, engineers can obtain new statistical data, the data is beneficial to deepening the understanding of the engineers on soil characteristic uncertainty, then the originally estimated soil characteristic statistical parameters are corrected, the correction of the parameters can further cause the reevaluation of the slope reliability, and particularly the problem of slope reliability updating caused by the uncertainty correction of the fluctuation range is solved. Some existing reliability analysis methods include a monte carlo simulation, a first order moment method, a probability moment point estimation method, a response surface and the like, wherein the monte carlo simulation method is a widely used analysis method. However, after the fluctuation range or other soil characteristic parameters are considered to be updated, the model is required to be reconstructed, and especially for the slope with small probability of failure, huge time, energy and computer resources are consumed by a large amount of calculation.
Disclosure of Invention
The invention aims to provide a slope stability evaluation and correction method based on local average. The method solves the problems that the Monte Carlo simulation method which is a widely applied reliability analysis method needs to rebuild the model after considering the fluctuation range or other soil characteristic parameters, and particularly for the slope with small probability of failure, a large amount of calculation consumes huge time, energy and computer resources.
The technical scheme of the invention is as follows: a slope stability evaluation and correction method based on local average comprises the following steps:
A. setting initial soil characteristic parameters to carry out Monte Carlo simulation, generating random samples, judging whether the slope is invalid or not through each group of samples, counting the number of invalid samples, and calculating the probability estimation value of the original slope invalid
Figure BDA0003155931280000021
B. Statistical failure sample x under original Monte Carlo simulationj,j=1,2,......,M;
C. Calculating the original and corrected combined probability density function according to the original and corrected soil characteristic parameters to obtain the weight index omegaj,j=1,2,......,M;
D. And calculating the corrected failure probability and the estimated value of the variation coefficient thereof.
In the slope stability evaluation and correction method based on local average, in the step C, it is assumed that m layers of soil exist in the slope, m soil variables correspond to the m soil variables, and the m soil variables are all subject to normal distribution, that is, the method is characterized in that
Figure BDA0003155931280000022
Dividing each layer of soil variable into n by local average1,n2,···nmAnd in the interval part, if the soil body variables of each layer are not related, the joint probability density function is as follows:
Figure BDA0003155931280000031
in the formula, xiN generated by local averaging of i-th layer soil variableiA vector representation of a random variable; mu.siN for i-th layer soil variable after local averagingiVector representation of individual variables to mean, i.e.
Figure BDA0003155931280000032
BiN for i-th layer soil variable after local averagingiOf a variable
Figure BDA0003155931280000033
The covariance matrix of (1) is similar to that when calculating the soil slope of one layer,
σij=σi[Γ(Δij|δ)]0.5 i=1,2,···m;j=1,2,···ni
in the formula, σijThe standard deviation corresponding to the jth variable after local averaging of the ith layer soil body variable is obtained; deltaijDividing the j-th interval of the corresponding sliding arc for the soil body variable in the ith layer of soil body of the side slope through local averaging; delta is the fluctuation range;
the correlation coefficient of the local average variable corresponding to the interval between the j section and the k section of the sliding arc in the i-th layer soil body is as follows:
Figure BDA0003155931280000034
in the formula,. DELTA.ij,ikFor the ith section delta on the sliding arc in the ith layer soil bodyiAnd j segment ΔjThe distance between the curves on the sliding arc.
In the slope stability evaluation and correction method based on local average, in the step C, according to the following formula:
Figure BDA0003155931280000041
defining:
Figure BDA0003155931280000042
when the value of i is equal to j,
Figure BDA0003155931280000043
the covariance matrix of the ith layer of soil slope can be obtained as follows:
Figure BDA0003155931280000051
thus, the joint probability density function value of the failed samples at the original fluctuation range δ can be found.
In the above slope stability evaluation and correction method based on local average, in step C, similarly under the multi-layer soil slope with the new fluctuation range δ', the joint probability density function is:
Figure BDA0003155931280000052
in the formula, xiN generated by local averaging of i-th layer soil variableiThe vector representation of random variable sample only needs to use original failure sample x(i);μiN is the i-th layer soil variable which is locally averaged under a new fluctuation range deltaiThe individual variables correspond to the vector representation of the mean, i.e. are still:
Figure BDA0003155931280000061
B′in for i-th layer soil body variable under new fluctuation range delta' through local averagingiA variable quantity
Figure BDA0003155931280000062
The covariance matrix of (a) may, similarly,
σ′ij=σ′i[Γ(Δij|δ′)]0.5 i=1,2,…m;j=1,2,…ni
wherein sigmaij'is the standard deviation corresponding to the jth variable after local averaging of the ith layer soil body variable in the new fluctuation range delta'; deltaijDividing the j-th interval of the corresponding sliding arc for the soil body variable in the ith soil body of the original lower side slope through local averaging; delta is the fluctuation range;
under the new fluctuation range delta', the correlation coefficient of the local average variable corresponding to the interval between the j section and the k section of the sliding arc in the i-th layer soil body is as follows:
Figure BDA0003155931280000063
wherein Δij,ikThe section I delta on the sliding arc in the layer I soil body under the original divisioniAnd ith segment deltajThe distance between the curves on the sliding arc; then:
Figure BDA0003155931280000064
when the value of i is equal to j,
Figure BDA0003155931280000065
the covariance matrix of the ith layer of soil slope can be obtained as follows:
Figure BDA0003155931280000071
the joint probability density function value of the failed samples under the new fluctuation range can be obtained.
In the slope stability evaluation and correction method based on local average, in the step D, the soil characteristic parameter estimation value is estimated based on new dataMaking a correction of the corrected failure probability
Figure BDA0003155931280000072
Expressed as:
Figure BDA0003155931280000073
combining with weighing theory, utilizing the corrected combined probability density function f (x) of the initial distribution of the soil characteristic parameters, and the corrected slope failure probability
Figure BDA0003155931280000074
Can be expressed as:
Figure BDA0003155931280000075
in the formula, omega is a weight index; thus, under the original joint probability density function f (x), the corrected failure probability estimation value
Figure BDA0003155931280000076
Can be expressed as:
Figure BDA0003155931280000081
random samples generated by original monte carlo simulation; omegajAnd j is 1,2 …, and M is a weight index of a failure sample generated by the original Monte Carlo simulation.
In the slope stability evaluation and correction method based on local average, in step D, the failure probability estimation value under the new fluctuation range δ' is obtained by combining the original joint probability density function value and according to the high-efficiency slope reliability correction formula combined with weighing theory:
Figure BDA0003155931280000082
wherein N is the Monte Carlo simulation times; n issThe number of failed samples.
The invention has the beneficial effects that: compared with the existing Monte Carlo simulation method, the method has the advantages that the spatial variability of the soil body parameters is described by considering the fluctuation range, the point variability and the spatial variability are connected, and the failure samples which consume a large amount of time for calculation in the original simulation can be used in calculating the corrected failure probability by deducing a formula, so that the time required by the corrected large amount of simulation calculation is avoided. Therefore, the method can effectively solve the problem of correcting the multilayer slope reliability under the consideration of the spatial variability of rock-soil body parameters; the method is irrelevant to the function when in correction calculation, so the method can be well suitable for the slope reliability correction problem of the implicit expression function. The variation coefficient of the failure probability can quantitatively represent the accuracy of the calculation result of the failure probability, and the variation coefficient of the failure probability after slope correction is smaller than that of the original calculation result, so the calculation precision of the method is superior to that of a direct Monte Carlo simulation method. When the failure sample space can cover the whole failure area, the corrected failure probability is more accurate. For example, for the correction condition that the parameter variation coefficient is changed from big to small in practical engineering, the method has high calculation efficiency and calculation precision, and is obviously superior to a Monte Carlo simulation method. The slope stability correction method does not need to execute Monte Carlo simulation again in the correction process, has high calculation efficiency, can shorten the calculation time to one percent of the calculation time of the primary slope, has clear concept and simple calculation, and is easy to master by engineers.
Drawings
FIG. 1 is a schematic diagram of a random field of characteristic parameters of a slope rock-soil mass;
FIG. 2 is a flow chart of the computing concept of the present invention;
FIG. 3 is a cross-sectional view of a soil slope;
FIG. 4 is an error analysis diagram.
Detailed Description
The invention is further illustrated by the following figures and examples, which are not to be construed as limiting the invention.
The embodiment of the invention comprises the following steps: a slope stability evaluation and correction method based on local average is shown in figure 1 and comprises the following steps:
A. setting initial soil characteristic parameters to carry out Monte Carlo simulation, generating random samples, judging whether the slope is invalid or not through each group of samples, counting the number of invalid samples, and calculating the probability estimation value of the original slope invalid
Figure BDA0003155931280000091
B. Statistical failure sample x under original Monte Carlo simulationj,j=1,2,......,M;
C. Calculating the original and corrected combined probability density function according to the original and corrected soil characteristic parameters to obtain the weight index omegaj,j=1,2,......,M;
D. And calculating the corrected failure probability and the estimated value of the variation coefficient thereof.
In the step C, assuming that m layers of soil bodies exist in the side slope, m soil body variables are correspondingly arranged and are all subjected to normal distribution, namely
Figure BDA0003155931280000092
Dividing each layer of soil variable into n by local average1,n2,···nmAnd in the interval part, if the soil body variables of each layer are not related, the joint probability density function is as follows:
Figure BDA0003155931280000101
in the formula, xiN generated by local averaging of i-th layer soil variableiA vector representation of a random variable; mu.siN for i-th layer soil variable after local averagingiVector representation of individual variables to mean, i.e.
Figure BDA0003155931280000102
BiN for i-th layer soil variable after local averagingiOf a variable
Figure BDA0003155931280000103
The covariance matrix of (1) is similar to that when calculating the soil slope of one layer,
σij=σi[Γ(Δij|δ)]0.5 i=1,2,···m;j=1,2,···ni
in the formula, σijThe standard deviation corresponding to the jth variable after local averaging of the ith layer soil body variable is obtained; deltaijDividing the j-th interval of the corresponding sliding arc for the soil body variable in the ith layer of soil body of the side slope through local averaging; δ is the fluctuation range.
The correlation coefficient of the local average variable corresponding to the interval between the j section and the k section of the sliding arc in the i-th layer soil body is as follows:
Figure BDA0003155931280000111
in the formula,. DELTA.ij,ikFor the ith section delta on the sliding arc in the ith layer soil bodyiAnd j segment ΔjThe distance between the curves on the sliding arc.
In step C, according to the following formula:
Figure BDA0003155931280000112
defining:
Figure BDA0003155931280000113
when the value of i is equal to j,
Figure BDA0003155931280000114
the covariance matrix of the ith layer of soil slope can be obtained as follows:
Figure BDA0003155931280000121
thus, the joint probability density function value of the failed samples at the original fluctuation range δ can be found.
In step C, similarly under the multi-layer soil slope of the new fluctuation range δ', the joint probability density function is:
Figure BDA0003155931280000122
in the formula, xiN generated by local averaging of i-th layer soil variableiThe vector representation of random variable sample only needs to use original failure sample x(i);μiN is the i-th layer soil variable which is locally averaged under a new fluctuation range deltaiThe individual variables correspond to the vector representation of the mean, i.e. are still:
Figure BDA0003155931280000131
B′in for i-th layer soil body variable under new fluctuation range delta' through local averagingiA variable quantity
Figure BDA0003155931280000132
The covariance matrix of (a) may, similarly,
σ′ij=σ′i[Γ(Δij|δ′)]0.5 i=1,2,…m;j=1,2,…ni
wherein sigmaij'is the standard deviation corresponding to the jth variable after local averaging of the ith layer soil body variable in the new fluctuation range delta'; deltaijDividing the j-th interval of the corresponding sliding arc for the soil body variable in the ith soil body of the original lower side slope through local averaging; δ is the fluctuation range.
Under the new fluctuation range delta', the correlation coefficient of the local average variable corresponding to the interval between the j section and the k section of the sliding arc in the i-th layer soil body is as follows:
Figure BDA0003155931280000133
wherein Δij,ikThe section I delta on the sliding arc in the layer I soil body under the original divisioniAnd ith segment deltajThe distance between the curves on the sliding arc. Then:
Figure BDA0003155931280000134
when the value of i is equal to j,
Figure BDA0003155931280000135
the covariance matrix of the ith layer of soil slope can be obtained as follows:
Figure BDA0003155931280000141
the joint probability density function value of the failed samples under the new fluctuation range can be obtained.
In step D, the soil body characteristic parameter estimated value is corrected based on the new data, and the corrected failure probability is
Figure BDA0003155931280000142
Expressed as:
Figure BDA0003155931280000143
combining with weighing theory, utilizing the corrected combined probability density function f (x) of the initial distribution of the soil characteristic parameters, and the corrected slope failure probability
Figure BDA0003155931280000144
Can be expressed as:
Figure BDA0003155931280000145
where ω is a weight index. Thus, under the original joint probability density function f (x), the corrected failure probability estimation value
Figure BDA0003155931280000146
Can be expressed as:
Figure BDA0003155931280000151
in the formula xiI is 1,2 …, N is a random sample generated by the original monte carlo simulation; omegajAnd j is 1,2 …, and M is a weight index of a failure sample generated by the original Monte Carlo simulation.
In the step D, combining the original joint probability density function value and according to a high-efficiency slope reliability correction formula combined with a weighing theory, the failure probability estimation value under the new fluctuation range delta' can be obtained as follows:
Figure BDA0003155931280000152
Figure BDA0003155931280000155
wherein N is the Monte Carlo simulation times; n issThe number of failed samples.
Specific application examples are as follows:
the calculation example is a simplified soil slope with the height H of 10m and the slope angle psi of 26 degrees, and soil body parameter characteristic values such as the volume weight gamma, the effective internal friction angle phi ', the effective cohesion force c' and the like are shown in table 1. Assuming that the plane sliding fails, the sliding plane inclination angle θ becomes 20 °. The functional function of the soil slope stability analysis can be expressed as:
Figure BDA0003155931280000153
when g is less than or equal to 0, the soil slope is unstable. The homogeneous side slope is shown in cross section in figure 3.
TABLE 1 statistical parameter Table of original random variables
Figure BDA0003155931280000154
After deep exploration and correction, the parameter change condition is shown in table 2.
TABLE 2 statistical parameter Table of modified random variables
Figure BDA0003155931280000161
Based on the established analytical model, by 1 × 108Sub Monte Carlo simulation, and calculating to obtain the failure probability of 5.033 × 10-3The failure probability after correction is calculated and obtained according to the method and is 7.67 multiplied by 10-4. The calculation of the initial Monte Carlo simulation in the EXCEL needs about 2h, the time needed after the correction is less than 2s, which shows that the method of the invention has high calculation efficiency, and if the correction is still carried out by repeatedly executing the Monte Carlo simulation after the correction, a large amount of calculation resources and time are needed, therefore. As shown in fig. 4, the Monte Carlo simulation and the comparative analysis of the method (as shown in fig. 4) by correcting the variation coefficient of the volume weight can show that the method has higher accuracy.

Claims (6)

1. A slope stability evaluation and correction method based on local average is characterized by comprising the following steps: the method comprises the following steps:
A. setting initial soil characteristic parameters to carry out Monte Carlo simulation, generating random samples, judging whether the slope is invalid or not through each group of samples, counting the number of invalid samples, and calculating the probability estimation value of the original slope invalid
Figure FDA0003155931270000011
B. Statistical failure sample x under original Monte Carlo simulationj,j=1,2,......,M;
C. Calculating the original and corrected combined probability density function according to the original and corrected soil characteristic parameters to obtain the weight index omegaj,j=1,2,......,M;
D. And calculating the corrected failure probability and the estimated value of the variation coefficient thereof.
2. The slope stability evaluation and correction method based on local average according to claim 1, characterized in that: in the step C, it is assumed that m layers of soil bodies exist in the side slope, m soil body variables correspond to the side slope, and the side slope are subjected to normal distribution, namely
Figure FDA0003155931270000012
Dividing each layer of soil variable into n by local average1,n2,···nmAnd in the interval part, if the soil body variables of each layer are not related, the joint probability density function is as follows:
Figure FDA0003155931270000013
in the formula, xiN generated by local averaging of i-th layer soil variableiA vector representation of a random variable; mu.siN for i-th layer soil variable after local averagingiVector representation of individual variables to mean, i.e.
Figure FDA0003155931270000021
BiN for i-th layer soil variable after local averagingiOf a variable
Figure FDA0003155931270000022
The covariance matrix of (1) is similar to that when calculating the soil slope of one layer,
σij=σi[Γ(Δij|δ)]0.5i=1,2,···m;j=1,2,···ni
in the formula, σijThe standard deviation corresponding to the jth variable after local averaging of the ith layer soil body variable is obtained; deltaijDividing the j-th interval of the corresponding sliding arc for the soil body variable in the ith layer of soil body of the side slope through local averaging; delta is the fluctuation range;
the correlation coefficient of the local average variable corresponding to the interval between the j section and the k section of the sliding arc in the i-th layer soil body is as follows:
Figure FDA0003155931270000023
in the formula,. DELTA.ij,ikFor the ith section delta on the sliding arc in the ith layer soil bodyiAnd j segment ΔjThe distance between the curves on the sliding arc.
3. The slope stability evaluation and correction method based on local average according to claim 2, characterized in that: in the step C, according to the following formula:
Figure FDA0003155931270000031
defining:
Figure FDA0003155931270000032
when the value of i is equal to j,
Figure FDA0003155931270000033
the covariance matrix of the ith layer of soil slope can be obtained as follows:
Figure FDA0003155931270000034
thus, the joint probability density function value of the failed samples at the original fluctuation range δ can be found.
4. The slope stability evaluation and correction method based on local average according to claim 3, characterized in that: in the step C, similarly under the multilayer soil slope of the new fluctuation range δ', the joint probability density function is:
Figure FDA0003155931270000041
in the formula, xiN generated by local averaging of i-th layer soil variableiThe vector representation of random variable sample only needs to use original failure sample x(i);μiN is the i-th layer soil variable which is locally averaged under a new fluctuation range deltaiThe individual variables correspond to the vector representation of the mean, i.e. are still:
μi'=(μi1i2,...μini)′=(μii,...μi)′,
B′in for i-th layer soil body variable under new fluctuation range delta' through local averagingiA variable quantity
Figure FDA0003155931270000042
The covariance matrix of (a) may, similarly,
σ′ij=σ′i[Γ(Δij|δ′)]0.5 i=1,2,…m;j=1,2,…ni
wherein sigmaij'is the standard deviation corresponding to the jth variable after local averaging of the ith layer soil body variable in the new fluctuation range delta'; deltaijDividing the j-th interval of the corresponding sliding arc for the soil body variable in the ith soil body of the original lower side slope through local averaging; delta is the fluctuation range;
under the new fluctuation range delta', the correlation coefficient of the local average variable corresponding to the interval between the j section and the k section of the sliding arc in the i-th layer soil body is as follows:
Figure FDA0003155931270000051
wherein Δij,ikThe section I delta on the sliding arc in the layer I soil body under the original divisioniAnd ith segment deltajThe distance between the curves on the sliding arc; then:
Figure FDA0003155931270000052
when the value of i is equal to j,
Figure FDA0003155931270000053
the covariance matrix of the ith layer of soil slope can be obtained as follows:
Figure FDA0003155931270000054
the joint probability density function value of the failed samples under the new fluctuation range can be obtained.
5. The slope stability evaluation and correction method based on local average according to claim 1, characterized in that: in the step D, the soil body characteristic parameter estimated value is corrected based on the new data, and the corrected failure probability is
Figure FDA0003155931270000061
Expressed as:
Figure FDA0003155931270000062
combining with weighing theory, utilizing the corrected combined probability density function f (x) of the initial distribution of the soil characteristic parameters, and the corrected slope failure probability
Figure FDA0003155931270000063
Can be expressed as:
Figure FDA0003155931270000064
in the formula, omega is a weight index; thus, under the original joint probability density function f (x), the correction is madeLater failure probability estimate
Figure FDA0003155931270000065
Can be expressed as:
Figure FDA0003155931270000066
in the formula xiI is 1,2 …, N is a random sample generated by the original monte carlo simulation; omegajAnd j is 1,2 …, and M is a weight index of a failure sample generated by the original Monte Carlo simulation.
6. The slope stability evaluation and correction method based on local average according to claim 5, characterized in that: in the step D, combining the original joint probability density function value and according to a high-efficiency slope reliability correction formula combined with a weighing theory, the failure probability estimation value under a new fluctuation range delta' can be obtained as follows:
Figure FDA0003155931270000067
wherein N is the Monte Carlo simulation times; n issThe number of failed samples.
CN202110777072.7A 2021-07-09 2021-07-09 Slope stability evaluation and correction method based on local average Pending CN113449429A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110777072.7A CN113449429A (en) 2021-07-09 2021-07-09 Slope stability evaluation and correction method based on local average

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110777072.7A CN113449429A (en) 2021-07-09 2021-07-09 Slope stability evaluation and correction method based on local average

Publications (1)

Publication Number Publication Date
CN113449429A true CN113449429A (en) 2021-09-28

Family

ID=77815643

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110777072.7A Pending CN113449429A (en) 2021-07-09 2021-07-09 Slope stability evaluation and correction method based on local average

Country Status (1)

Country Link
CN (1) CN113449429A (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106021853A (en) * 2016-05-09 2016-10-12 青岛理工大学 A method for developing slope stability reliability analysis software under spatial variation characteristics
CN109446616A (en) * 2018-10-18 2019-03-08 东北电力大学 A kind of homogeneous slope method for evaluating reliability
CN111339488A (en) * 2020-03-24 2020-06-26 成都理工大学 Method for calculating failure probability of slope system based on Kriging
CN111444649A (en) * 2020-03-24 2020-07-24 成都理工大学 Slope system reliability analysis method based on intensity reduction method
KR20210030031A (en) * 2019-09-09 2021-03-17 한국지질자원연구원 Early warning method and system for landslide integrated with weather forecasting information
CN113033108A (en) * 2021-04-19 2021-06-25 昆明理工大学 Side slope reliability judgment method based on AdaBoost algorithm

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106021853A (en) * 2016-05-09 2016-10-12 青岛理工大学 A method for developing slope stability reliability analysis software under spatial variation characteristics
CN109446616A (en) * 2018-10-18 2019-03-08 东北电力大学 A kind of homogeneous slope method for evaluating reliability
KR20210030031A (en) * 2019-09-09 2021-03-17 한국지질자원연구원 Early warning method and system for landslide integrated with weather forecasting information
CN111339488A (en) * 2020-03-24 2020-06-26 成都理工大学 Method for calculating failure probability of slope system based on Kriging
CN111444649A (en) * 2020-03-24 2020-07-24 成都理工大学 Slope system reliability analysis method based on intensity reduction method
CN113033108A (en) * 2021-04-19 2021-06-25 昆明理工大学 Side slope reliability judgment method based on AdaBoost algorithm

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
张浮平 等: "基于蒙特卡罗模拟的高效边坡可靠度修正方法", 《工程力学》, vol. 33, no. 7, 25 July 2016 (2016-07-25) *
程井 等: "基于Neumann展开随机有限元的混凝土重力坝结构可靠度分析", 《水利水电科技进展》, vol. 39, no. 2, 10 March 2019 (2019-03-10) *

Similar Documents

Publication Publication Date Title
Liu et al. Investigation of slope failure mode evolution during large deformation in spatially variable soils by random limit equilibrium and material point methods
Cho Effects of spatial variability of soil properties on slope stability
Goh et al. Reliability assessment of stability of underground rock caverns
Cayol et al. 3D mixed boundary elements for elastostatic deformation field analysis
Dai et al. A probabilistic slope stability analysis using deterministic computer software
CN111666621B (en) Method for determining safe support pressure interval of excavation face of tunnel in clay stratum
Woodward et al. Contact stresses for closely conforming bodies--application to cylinders and spheres
CN115033941A (en) Method for calculating contact area of three-dimensional rough surface
CN113239577B (en) Quantitative evaluation method for borehole instability risk of drilling in complex stratum
Wang et al. A semi-analytical one-dimensional model for offshore pile foundations considering effects of pile diameter and aspect ratio
CN111859303B (en) Soil humidity fusion method and system based on dynamic Bayesian average
CN113722920A (en) Rapid slope earthquake reliability calculation method based on FLAC3D-Python secondary development
Huang et al. Rainfall reliability evaluation for stability of municipal solid waste landfills on slope
CN112115530B (en) Method for calculating point failure probability of soil slope
CN113449429A (en) Slope stability evaluation and correction method based on local average
Rao et al. Factorized high dimensional model representation for structural reliability analysis
Morse et al. Manufacturing cost and reliability‐based shape optimization of plate structures
Vazirizade An intelligent integrated method for reliability estimation of offshore structures wave loading applied in time domain
CN108763611A (en) A kind of wing structure random eigenvalue analysis method based on probabilistic density evolution
Ryan et al. Influence of soil plasticity models on offshore wind turbine response
Gaxiola-Camacho et al. Novel concepts for reliability analysis of dynamic structural systems
Li et al. Research on the equivalent stiffness of bucket foundations for offshore wind power
Hashemi et al. Application of random radial point interpolation method to foundations bearing capacity considering progressive failure
Zhao et al. Prediction of landslide displacement based on kernel principal component analysis and neural network-markov chain
HARIRI et al. Comparative study of endurance time and time history methods in seismic analysis of high arch dams

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