CN110261573B - Dynamic evaluation method for stability of high-position rocky landslide - Google Patents

Dynamic evaluation method for stability of high-position rocky landslide Download PDF

Info

Publication number
CN110261573B
CN110261573B CN201910406610.4A CN201910406610A CN110261573B CN 110261573 B CN110261573 B CN 110261573B CN 201910406610 A CN201910406610 A CN 201910406610A CN 110261573 B CN110261573 B CN 110261573B
Authority
CN
China
Prior art keywords
landslide
soil
rock
model
solid
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.)
Active
Application number
CN201910406610.4A
Other languages
Chinese (zh)
Other versions
CN110261573A (en
Inventor
许建聪
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tongji University
Original Assignee
Tongji University
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 Tongji University filed Critical Tongji University
Priority to CN201910406610.4A priority Critical patent/CN110261573B/en
Publication of CN110261573A publication Critical patent/CN110261573A/en
Application granted granted Critical
Publication of CN110261573B publication Critical patent/CN110261573B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/24Earth materials
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06Power analysis or power optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Biochemistry (AREA)
  • Pathology (AREA)
  • Medicinal Chemistry (AREA)
  • Remote Sensing (AREA)
  • Analytical Chemistry (AREA)
  • Geology (AREA)
  • General Health & Medical Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Immunology (AREA)
  • Food Science & Technology (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Environmental & Geological Engineering (AREA)
  • Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

The invention relates to a high-order rock landslide stability dynamic numerical evaluation method, which comprises the following steps: s1, establishing a soil gas-liquid-solid multi-field coupling sub-model; s2, establishing a viscoelastic-plastic model of soil body gas-liquid-solid multi-field coupling; s3, establishing a nonlinear variable damage mechanical model; s4, establishing a gas-liquid-solid-rheological damage multi-field coupling model; s5, obtaining the maximum displacement of the landslide, and calculating the high-position rock landslide safety coefficient under the fluid-solid coupling effect; s6, dynamically evaluating the stability of the high rock landslide: firstly, safety judgment is carried out, and then stability judgment is carried out. Compared with the prior art, the gas-liquid-solid-rheological damage multi-field coupling model established by the method is more suitable for the actual situation of high-position rock landslide, combines the dynamic maximum displacement of the landslide and the landslide safety coefficient, and is beneficial to accurately and effectively evaluating the landslide stability in time after rainfall.

Description

Dynamic evaluation method for stability of high-position rocky landslide
Technical Field
The invention relates to the field of rock landslide and slope engineering, in particular to a high-order rock landslide stability dynamic evaluation method.
Background
The high-position rocky landslide has the characteristics of high disaster identification difficulty, extremely high concealment and outbreak suddenness of inoculation, peculiar disaster forming mode, extremely high harmfulness and high disaster prevention and reduction difficulty. In recent decades, with the reduction of the weather recurrence period such as extreme heavy rainfall, the tendency of extra-large disasters caused by high-position rocky landslides in China is gradually increased, and in China, high-position rocky landslides such as 24-day rural landslides in 6 months in 2017, 24-day rural landslides in 24 days, tunneling-Changshu-village landslides in 9 months in 2016, river-Wuli-slope landslides in 7 months and 10 days in 2013, 27-day birch-mountain landslides in 2010, 28-day guan-mountain torres in 6 months in 2010, 5-day Longglong-Jiwei-mountain landslides in 6 months in 2009 and the like all cause a great amount of casualties and huge property loss.
Due to the diversity and complexity of high-position rock landslide rock-soil body media, the vegetation condition on the slope, the type and thickness of a shallow surface soil body of the landslide body, the rock body type in the landslide body, the development condition, the opening degree and the filling condition of a structural surface and a crack in the rock body, the type, the porosity, the water content, the thickness and the spatial distribution arrangement of a weak structural surface or a weak interlayer soil body, the physical mechanical property, the permeability and the rheological property of the rock-soil body in the slope can also have obvious difference, so that the stability of the landslide is difficult to determine in the actual engineering.
Therefore, scholars at home and abroad carry out a great deal of related research work: on one hand, the research on the stability of the rock slope comprises the steps of researching the stability of the jointed rock slope by adopting a finite element and discrete element coupling method of a random structural surface and an elastic-plastic finite element or finite difference strength reduction method based on fracture mechanics and interface elements; the influence of fracture seepage on the stability of the rock slope is researched by adopting a fluid-solid coupled DDA method; dynamically evaluating the stability of the soft rock slope by adopting a finite difference method based on a western original creep model; establishing an equivalent rheological damage model of the jointed rock mass to evaluate the stability of the unloading slope; researching the deformation and damage of the silicious claystone forward slope by adopting a discrete element and Burgers creep model; establishing a two-dimensional fluid-solid coupling seepage model to analyze the stability of the mud-rock interlayer detritus rock slope under different rainfall conditions;
on the other hand, the stability of the high-position rock landslide is evaluated by applying a key block control theory, a fracture mechanics theory and a hydromechanics theory, combining a strip method, an engineering geomechanical method, a thrust transfer coefficient method, a surface deformation sign and structural surface geomechanical comprehensive analysis method and adopting a red-flat projection method, a wedge vector analysis method, an analytic method based on a soft foundation effect, a rigid body limit balance method, a resistance body decomposition method, a discrete element method, a finite difference method, a finite element method, a one-dimensional dynamic viscoplasticity method and the like.
However, the current research rarely involves dynamic evaluation of the stability of the high-position rock landslide, lacks multi-physical-field coupling analysis of the landslide, and cannot guarantee the accuracy of the evaluation of the stability of the high-position rock landslide.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provide a dynamic evaluation method for the stability of high-position rock landslide.
The purpose of the invention can be realized by the following technical scheme: a high-order rock landslide stability dynamic numerical evaluation method comprises the following steps:
s1, obtaining a permeability coefficient and a suction force calculation parameter through a soil unsaturated/saturated soil permeability test, and establishing a soil gas-liquid-solid multi-field coupling sub-model;
s2, establishing a viscoelastic-plastic model of soil gas-liquid-solid multi-field coupling by a soil unsaturated/saturated soil triaxial rheological test and applying unsaturated soil mechanics, a general effective stress theory of saturated-unsaturated soil, a rheological principle of rock-soil mechanics, rheological damage mechanics and a multi-field coupling theory;
s3, establishing a nonlinear variable damage mechanical model according to the soil gas-liquid-solid multi-field coupling sub-model and the viscoelastic-plastic model;
s4, embedding the soil gas-liquid-solid multi-field coupling sub-model, the viscoelastic-plastic model and the nonlinear variable damage mechanics model into an ABAQUS finite element and discrete element coupling numerical analysis platform, and establishing a gas-liquid-solid-rheological damage multi-field coupling model;
s5, acquiring the maximum displacement of the landslide based on the gas-liquid-solid-rheological damage multi-field coupling model, and calculating the high-position rock landslide safety coefficient under the fluid-solid coupling effect by adopting a fluid-solid coupling numerical value flow element method and a three-dimensional viscoelastic-plastic finite element intensity reduction method;
s6, dynamically evaluating the stability of the high rock landslide: judging whether the high-position rock landslide is safe or not according to the maximum landslide displacement and the landslide safety coefficient, and judging whether the high-position rock landslide is stable or not according to the quantity of safe and unsafe results, wherein if the maximum landslide displacement is smaller than or equal to a first preset threshold value, and the landslide safety coefficient is larger than or equal to a second preset threshold value, the high-position rock landslide is judged to be safe, otherwise, the high-position rock landslide is judged to be unsafe;
and finally judging that the high-order rock landslide is stable if the number of the results judged to be safe is more than or equal to M times of the number of the results judged to be unsafe, otherwise, judging that the high-order rock landslide is unstable, wherein M is more than or equal to 3 and less than or equal to 5.
Preferably, the specific process of establishing the soil gas-liquid-solid multi-field coupling sub-model in step S1 is as follows:
s11, obtaining a permeability coefficient and a suction force calculation parameter through a non-saturated/saturated soil permeability test of high-position rock landslide surface soil and a sliding belt soil body;
s12, according to the momentum balance equation and the energy balance equation of the rock-soil mass and the momentum, mass and energy conservation law of continuous medium mechanics, the phase change of water and gas, the gas solubility in liquid, the water transfer and expansion are considered at the same time, and the multi-field coupling mathematical model of the gas field, the seepage field, the stress field and the deformation field of the soil mass is established by taking the solid displacement, the capillary pressure, the pore water pressure, the pore air pressure (dry air pressure, steam pressure and the like) and the porosity as unknown variables.
Preferably, the viscoelastic-plastic model in step S2 is a viscoelastic-plastic model of clay or clay containing crushed stone, which takes into account multi-field coupling of an air field, a seepage field, a stress field and a deformation field.
Preferably, the nonlinear damage mechanics model in step S3 includes rainfall intensity, rainfall pattern, pore air pressure, soil moisture content, permeability coefficient, matrix suction, shear strength and viscosity variable parameters.
Preferably, the gas-liquid-solid-rheological damage multi-field coupling model in step S4 is specifically a gas-liquid-solid-rheological damage multi-field coupling model with variable viscosity, and the variable viscosity changes correspondingly with the difference of shear stress, water content and time.
Preferably, the fluid-solid coupling numerical value flow element method in step S5 considers rainfall infiltration, joint crack opening change and crack propagation.
Preferably, the step S5 specifically includes the following steps:
s51, obtaining fitting constants E and eta according to the test data through the triaxial rheological test of the step S2 to establish a rheological model of the weak interlayer rock-soil mass:
Figure GDA0003147830560000031
in the formula, Δ εcrFor creep increase, E is the modulus of elasticity, Δ t is the time increase, η is the viscosity coefficient,
Figure GDA0003147830560000032
to be the strain rate, σ is the shear stress,
Figure GDA0003147830560000033
Δ σ is the shear stress increment, which is the partial derivative of the strain rate versus shear stress.
S52, based on the rheological model of the step S51, adopting a three-dimensional viscoelastoplastic finite element strength reduction method to resist shear strength indexes C and C
Figure GDA0003147830560000034
Performing intensity reduction calculation to obtain a reduced corresponding virtual shear strength index CFAnd
Figure GDA0003147830560000035
CF=C/FS
Figure GDA0003147830560000036
wherein C is the cohesive force of rock-soil mass,
Figure GDA0003147830560000041
is the internal friction angle of the rock-soil mass, FSShear strength index C and
Figure GDA0003147830560000042
the reduction factor of (c).
S53, using the virtual shear strength index CFAnd
Figure GDA0003147830560000043
substituted shear strength index C and
Figure GDA0003147830560000044
virtual shear strength index C of each time step by ABAQUS softwareFAnd
Figure GDA0003147830560000045
calculating the stability of the finite element, and continuously reducing until the stability of the finite element is not converged, wherein the reduction coefficient F isSNamely the landslide safety factor.
Preferably, the values of the shear stress and the shear stress increment in step S51 are iteratively determined by finite element stability calculation.
Preferably, the step S6 specifically includes the following steps:
s61, safety evaluation was performed every T hours for N days after rainfall: if the maximum displacement of the landslide is smaller than or equal to a first preset threshold value and the landslide safety coefficient is larger than or equal to a second preset threshold value, judging the high-position rock landslide to be safe, otherwise, judging the high-position rock landslide to be unsafe, and obtaining A results judged to be safe and B results judged to be unsafe, wherein N is 5, and T is 1;
and S62, if A is larger than or equal to B M, judging the high rock landslide to be stable, otherwise, judging the high rock landslide to be unstable, wherein M is 4.
Compared with the prior art, the invention has the following advantages:
firstly, the dynamic evaluation of the stability of the high-position rock landslide is carried out after rainfall by establishing a gas-liquid-solid-rheological damage multi-field coupling model and combining the dynamic maximum displacement of the landslide and the landslide safety coefficient, so that the evaluation result of the method is more accurate and effective, and powerful data reference is provided for predicting the occurrence of landslide disasters.
The invention establishes the variable viscosity gas-liquid-solid-rheological damage multi-field coupling model, and the reliability of the multi-physical-field coupling model is ensured because the variable viscosity changes along with the shearing stress, the water content and the time.
And thirdly, the multi-physical-field coupling model comprises a nonlinear rheological damage model fused with variables such as rainfall intensity, rainfall type, pore air pressure, soil water content, permeability coefficient, matrix suction, shear strength, viscosity and the like, is more suitable for the actual situation of the high-order rock landslide, and is beneficial to evaluating the stability of the landslide in time after rainfall.
Drawings
FIG. 1 is a schematic flow diagram of the process of the present invention;
fig. 2 is a schematic view of a dynamic evaluation flow of the embodiment.
Detailed Description
The invention is described in detail below with reference to the figures and specific embodiments.
The invention is subsidized by a capital project (U1765110) of Yajianjiang united fund project (U1765110) of the national science foundation committee-Yajianjiang basin hydropower development Limited company and a capital project (22120180312) of basic scientific research service charge of Central colleges, as shown in figure 1, a dynamic evaluation method for stability of high-grade rock landslide comprises the following steps:
s1, obtaining a permeability coefficient and a suction force calculation parameter through a soil unsaturated/saturated soil permeability test, and establishing a soil gas-liquid-solid multi-field coupling sub-model;
s2, establishing a viscoelastic-plastic model of soil gas-liquid-solid multi-field coupling by a soil unsaturated/saturated soil triaxial rheological test and applying unsaturated soil mechanics, a general effective stress theory of saturated-unsaturated soil, a rheological principle of rock-soil mechanics, rheological damage mechanics and a multi-field coupling theory;
s3, establishing a nonlinear variable damage mechanical model according to the soil gas-liquid-solid multi-field coupling sub-model and the viscoelastic-plastic model;
s4, embedding the soil gas-liquid-solid multi-field coupling sub-model, the viscoelastic-plastic model and the nonlinear variable damage mechanics model into an ABAQUS finite element and discrete element coupling numerical analysis platform, and establishing a gas-liquid-solid-rheological damage multi-field coupling model;
s5, acquiring the maximum displacement of the landslide based on the gas-liquid-solid-rheological damage multi-field coupling model, and calculating the high-position rock landslide safety coefficient under the fluid-solid coupling effect by adopting a fluid-solid coupling numerical value flow element method and a three-dimensional viscoelastic-plastic finite element intensity reduction method;
s6, dynamically evaluating the stability of the high rock landslide: judging whether the high-position rock landslide is safe or not according to the maximum landslide displacement and the landslide safety coefficient, and judging whether the high-position rock landslide is stable or not according to the quantity of safe and unsafe results, wherein if the maximum landslide displacement is less than or equal to a first preset threshold value, and the landslide safety coefficient is greater than or equal to a second preset threshold value, the high-position rock landslide is judged to be safe, otherwise, the high-position rock landslide is judged to be unsafe;
and finally judging that the high-order rock landslide is stable if the number of the results judged to be safe is more than or equal to M times of the number of the results judged to be unsafe, otherwise, judging that the high-order rock landslide is unstable, wherein M is more than or equal to 3 and less than or equal to 5.
Wherein, step S1 includes:
s11, obtaining a permeability coefficient and a suction force calculation parameter through a non-saturated/saturated soil permeability test of high-position rock landslide surface soil and a sliding belt soil body;
s12, according to the momentum balance equation and the energy balance equation of the rock-soil mass and the momentum, mass and energy conservation law of continuous medium mechanics, the phase change of water and gas, the gas solubility in liquid, the water transfer and expansion are considered at the same time, and the multi-field coupling mathematical model of the gas field, the seepage field, the stress field and the deformation field of the soil mass is established by taking the solid displacement, the capillary pressure, the pore water pressure, the pore air pressure (dry air pressure, steam pressure and the like) and the porosity as unknown variables.
Step S2 includes: a viscoelastic-plastic model of cohesive soil or cohesive soil containing broken stones is established by applying unsaturated soil mechanics, a general effective stress theory of saturated-unsaturated soil, a rheological principle of geotechnical mechanics, rheological damage mechanics and a multi-field coupling theory, wherein the multi-field coupling of an air field, a seepage field, a stress field and a deformation field is considered.
Step S3 includes: and fusing the nonlinear rheological damage of rainfall intensity, rainfall type, pore air pressure, soil water content, permeability coefficient, matrix suction, shear strength and viscosity variable to establish a nonlinear rheological damage mechanical model.
Step S4 includes: embedding the models established in the steps S1, S2 and S3 on an ABAQUS finite element and discrete element coupling numerical analysis platform to establish a gas-liquid-solid-rheological damage multi-field coupling model with variable viscosity (the viscosity changes along with shear stress, water content and time).
In the step S5, the fluid-solid coupling numerical value flow element method considers rainfall infiltration, joint crack opening change and crack propagation.
The concrete process of calculating the landslide safety coefficient in the step S5 is as follows:
s51, obtaining fitting constants E and eta according to the test data through the triaxial rheological test of the step S2 to establish a rheological model of the weak interlayer rock-soil mass:
Figure GDA0003147830560000061
in the formula, Δ εcrFor creep increase, E is the modulus of elasticity, Δ t is the time increase, η is the viscosity coefficient,
Figure GDA0003147830560000062
to be the strain rate, σ is the shear stress,
Figure GDA0003147830560000063
for the partial derivative of strain rate versus shear stress, Δ σ is the shear stress delta, and the values of shear stress and shear stress delta are determined iteratively by finite element stability calculations, E — 2.36995GPa and η — 61.43937GPa · d in this example.
S52, based on the rheological model of the step S51, adopting a three-dimensional viscoelastoplastic finite element strength reduction method to resist shear strength indexes C and C
Figure GDA0003147830560000064
Performing intensity reduction calculation to obtain the corresponding value after reductionVirtual shear strength index CFAnd
Figure GDA0003147830560000065
CF=C/FS
Figure GDA0003147830560000066
wherein C is the cohesive force of rock-soil mass,
Figure GDA0003147830560000067
is the internal friction angle of the rock-soil mass, FSShear strength index C and
Figure GDA0003147830560000068
the reduction factor of (c).
S53, using the virtual shear strength index CFAnd
Figure GDA0003147830560000069
substituted shear strength index C and
Figure GDA00031478305600000610
virtual shear strength index C of each time step by ABAQUS softwareFAnd
Figure GDA00031478305600000611
calculating the stability of the finite element, and continuously reducing until the stability of the finite element is not converged, wherein the reduction coefficient F isSNamely the landslide safety factor.
Step S6 includes: and dynamically evaluating the stability of the high-grade rock landslide according to the gas-liquid-solid-rheological damage multi-field coupling model established in the step S4 and the landslide safety coefficient obtained in the step S5. In this embodiment, the safety evaluation is performed every 1 hour 5 days after rainfall and rain, and the landslide stability determination is performed according to the safety evaluation result, as shown in fig. 2, the specific dynamic evaluation process of this embodiment is:
s61, safety evaluation was performed every 1 hour at 5 days after rainfall: when the maximum displacement of the landslide is less than or equal to 0.3m and the landslide safety coefficient is more than or equal to 1.05, judging the landslide to be safe, and judging the landslide to be unsafe, thereby obtaining 5 × 24-120 evaluation results which comprise A results judged to be safe and B results judged to be unsafe;
and S62, judging whether A is more than or equal to B4, if so, judging that the high rock landslide is stable, otherwise, judging that the high rock landslide is unstable.

Claims (8)

1. A high-order rock landslide stability dynamic numerical evaluation method is characterized by comprising the following steps:
s1, obtaining a permeability coefficient and a suction force calculation parameter through a soil unsaturated/saturated soil permeability test, and establishing a soil gas-liquid-solid multi-field coupling sub-model;
s2, establishing a viscoelastic-plastic model of soil gas-liquid-solid multi-field coupling by a soil unsaturated/saturated soil triaxial rheological test and applying unsaturated soil mechanics, a general effective stress theory of saturated-unsaturated soil, a rheological principle of rock-soil mechanics, rheological damage mechanics and a multi-field coupling theory;
s3, establishing a nonlinear variable damage mechanical model according to the soil gas-liquid-solid multi-field coupling sub-model and the viscoelastic-plastic model;
s4, embedding the soil gas-liquid-solid multi-field coupling sub-model, the viscoelastic-plastic model and the nonlinear variable damage mechanics model into an ABAQUS finite element and discrete element coupling numerical analysis platform, and establishing a gas-liquid-solid-rheological damage multi-field coupling model;
s5, acquiring the maximum displacement of the landslide based on the gas-liquid-solid-rheological damage multi-field coupling model, and calculating the high-position rock landslide safety coefficient under the fluid-solid coupling effect by adopting a fluid-solid coupling numerical value flow element method and a three-dimensional viscoelastic-plastic finite element intensity reduction method;
s6, dynamically evaluating the stability of the high rock landslide: judging whether the high-position rock landslide is safe or not according to the maximum landslide displacement and the landslide safety coefficient, and judging whether the high-position rock landslide is stable or not according to the quantity of safe and unsafe results, wherein if the maximum landslide displacement is less than or equal to a first preset threshold value, and the landslide safety coefficient is greater than or equal to a second preset threshold value, the high-position rock landslide is judged to be safe, otherwise, the high-position rock landslide is judged to be unsafe;
and finally judging that the high-order rock landslide is stable if the number of the results judged to be safe is more than or equal to M times of the number of the results judged to be unsafe, otherwise, judging that the high-order rock landslide is unstable, wherein M is more than or equal to 3 and less than or equal to 5.
2. The method for evaluating the dynamic numerical value of the stability of the high-order rocky landslide according to claim 1, wherein the specific process of establishing the soil gas-liquid-solid multi-field coupling submodel in the step S1 is as follows:
s11, obtaining a permeability coefficient and a suction force calculation parameter through a non-saturated/saturated soil permeability test of high-position rock landslide surface soil and a sliding belt soil body;
s12, according to the momentum balance equation and the energy balance equation of the rock-soil mass and the momentum, mass and energy conservation law of continuous medium mechanics, the phase change of water and gas, the gas solubility in liquid, the water transfer and expansion are considered at the same time, and the multi-field coupling mathematical model of the gas field, the seepage field, the stress field and the deformation field of the soil mass is established by taking the solid displacement, the capillary pressure, the pore water pressure, the pore air pressure and the porosity as unknown variables.
3. The method for dynamically evaluating the stability of the high-position rocky landslide according to claim 1, wherein the viscoelastic-plastic model in step S2 is a viscoelastic-plastic model of clay or clay containing crushed stone, which is obtained by multi-field coupling of an air field, a seepage field, a stress field and a deformation field.
4. The method as claimed in claim 1, wherein the nonlinear damage-variable mechanical model in step S3 includes parameters of rainfall intensity, rainfall pattern, pore pressure, soil moisture content, permeability coefficient, matrix suction, shear strength and viscosity.
5. The method for evaluating the dynamic numerical value of the stability of the high rock landslide of claim 1, wherein the gas-liquid-solid-rheological damage multi-field coupling model in the step S4 is a gas-liquid-solid-rheological damage multi-field coupling model with variable viscosity, and the variable viscosity changes correspondingly with different shearing stress, water content and time.
6. The method for dynamically evaluating the stability of the high-order rocky landslide according to claim 1, wherein the fluid-solid coupling numerical flow element method in step S5 takes into account rainfall infiltration, joint fracture opening change and crack propagation.
7. The method for evaluating the dynamic numerical value of the stability of the high rock landslide according to claim 1, wherein the step S5 specifically comprises the steps of:
s51, obtaining fitting constants E and eta according to the test data through the triaxial rheological test of the step S2 to establish a rheological model of the weak interlayer rock-soil mass:
Figure FDA0003147830550000021
in the formula, Δ εcrFor creep increase, E is the modulus of elasticity, Δ t is the time increase, η is the viscosity coefficient,
Figure FDA0003147830550000022
to be the strain rate, σ is the shear stress,
Figure FDA0003147830550000023
is the partial derivative of strain rate to shear stress, Δ σ is the shear stress increment, wherein the values of shear stress σ and shear stress increment Δ σ are iteratively determined by finite element stability calculations;
s52, based on the rheological model of the step S51, adopting three-dimensional limited viscoelastoplasticityReduction of original strength, shear strength index C and
Figure FDA0003147830550000024
performing intensity reduction calculation to obtain a reduced corresponding virtual shear strength index CFAnd
Figure FDA0003147830550000025
CF=C/FS
Figure FDA0003147830550000026
wherein C is the cohesive force of rock-soil mass,
Figure FDA0003147830550000027
is the internal friction angle of the rock-soil mass, FSShear strength index C and
Figure FDA0003147830550000028
the reduction factor of (d);
s53, using the virtual shear strength index CFAnd
Figure FDA0003147830550000031
substituted shear strength index C and
Figure FDA0003147830550000032
virtual shear strength index C of each time step by ABAQUS softwareFAnd
Figure FDA0003147830550000033
calculating the stability of the finite element, and continuously reducing until the stability of the finite element is not converged, wherein the reduction coefficient F isSNamely the landslide safety factor.
8. The method for evaluating the dynamic numerical value of the stability of the high rock landslide according to claim 1, wherein the step S6 specifically comprises the steps of:
s61, safety evaluation was performed every T hours for N days after rainfall: if the maximum displacement of the landslide is smaller than or equal to a first preset threshold value and the landslide safety coefficient is larger than or equal to a second preset threshold value, judging the high-position rock landslide to be safe, otherwise, judging the high-position rock landslide to be unsafe, and obtaining A results judged to be safe and B results judged to be unsafe, wherein N is 5, and T is 1;
and S62, if A is larger than or equal to B M, judging the high rock landslide to be stable, otherwise, judging the high rock landslide to be unstable, wherein M is 4.
CN201910406610.4A 2019-05-16 2019-05-16 Dynamic evaluation method for stability of high-position rocky landslide Active CN110261573B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910406610.4A CN110261573B (en) 2019-05-16 2019-05-16 Dynamic evaluation method for stability of high-position rocky landslide

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910406610.4A CN110261573B (en) 2019-05-16 2019-05-16 Dynamic evaluation method for stability of high-position rocky landslide

Publications (2)

Publication Number Publication Date
CN110261573A CN110261573A (en) 2019-09-20
CN110261573B true CN110261573B (en) 2021-09-03

Family

ID=67913232

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910406610.4A Active CN110261573B (en) 2019-05-16 2019-05-16 Dynamic evaluation method for stability of high-position rocky landslide

Country Status (1)

Country Link
CN (1) CN110261573B (en)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110851972B (en) * 2019-11-04 2021-08-31 武汉大学 Rock-soil body structure random fracture simulation method and device based on Gaussian random field
CN111159907B (en) * 2019-12-31 2024-04-23 同济大学 Slope simulation method under action of storm
CN111310736B (en) * 2020-03-26 2023-06-13 上海同岩土木工程科技股份有限公司 Rapid identification method for unloading and stacking of vehicles in protection area
CN111581704B (en) * 2020-05-12 2023-09-08 江苏省地质调查研究院 Comprehensive evaluation method for stability and greening effect of soil matrix of rock slope
CN111861107B (en) * 2020-06-08 2021-03-30 中国地质大学(武汉) Landslide dynamic stability evaluation method based on intensity reduction method
CN112129619B (en) * 2020-09-10 2021-08-17 武汉大学 Deep soft rock quality evaluation method and system
CN113408184B (en) * 2021-04-21 2022-01-28 中国地质大学(武汉) Anchoring slope safety evaluation method based on genetic algorithm and discrete element analysis method
CN113268874B (en) * 2021-05-26 2022-03-22 河北工业大学 LBM-DDA coupling-based accumulation body seepage erosion damage simulation calculation method
CN113505514B (en) * 2021-08-04 2024-01-05 大连海事大学 Rock mass elastoplastic damage-seepage coupling calculation method under complex disturbance condition
CN115221727B (en) * 2022-08-05 2023-09-12 中交第三航务工程局有限公司 Numerical simulation model parameter determination method of rock mass based on water content
CN115359630A (en) * 2022-08-22 2022-11-18 重庆地质矿产研究院 Landslide early warning method based on XGB and Catboost dual tree lifting algorithm

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102841191A (en) * 2012-08-23 2012-12-26 苏州南智传感科技有限公司 Soil slope multi-field coupling simulation testing device
CN103149340A (en) * 2013-02-02 2013-06-12 青岛理工大学 Power monitoring method for measuring landslide stability by using rainfall
CN103646181A (en) * 2013-12-20 2014-03-19 青岛理工大学 Determination method for stability coefficient and early warning criterion of creeping type artificial slope
CN103822835A (en) * 2014-03-18 2014-05-28 中国矿业大学 Method for separating viscosity-elasticity-plasticity characteristic parameters in triaxial rheological test on rocks
CN104881583A (en) * 2015-06-04 2015-09-02 唐晓松 Multi-means, dynamic and whole-process landslide prewarning method
CN104964878A (en) * 2015-07-14 2015-10-07 中国科学院武汉岩土力学研究所 Triaxial test system and method for unsaturated soil multi-field coupling
CN105139585A (en) * 2015-07-29 2015-12-09 同济大学 Intelligent early warning and prediction method for soil slope danger
WO2016027390A1 (en) * 2014-08-21 2016-02-25 日本電気株式会社 Slope monitoring system, device for slope safety analysis, method, and program
CN105912753A (en) * 2016-04-01 2016-08-31 国家海洋局第二海洋研究所 Method for analyzing three-dimensional stability of submarine slope based on strength reduction method
CN107169271A (en) * 2017-04-26 2017-09-15 中国地质大学(武汉) A kind of stability dynamic evaluation method of strain softening side slope
CN108538026A (en) * 2018-02-23 2018-09-14 青岛理工大学 Multi-index comprehensive evaluation and early warning method for slope stability
CN108733961A (en) * 2018-06-12 2018-11-02 河海大学 A kind of Slope Sliding Convenient stable criterion considering non-saturated region effect

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102841191A (en) * 2012-08-23 2012-12-26 苏州南智传感科技有限公司 Soil slope multi-field coupling simulation testing device
CN103149340A (en) * 2013-02-02 2013-06-12 青岛理工大学 Power monitoring method for measuring landslide stability by using rainfall
CN103646181A (en) * 2013-12-20 2014-03-19 青岛理工大学 Determination method for stability coefficient and early warning criterion of creeping type artificial slope
CN103822835A (en) * 2014-03-18 2014-05-28 中国矿业大学 Method for separating viscosity-elasticity-plasticity characteristic parameters in triaxial rheological test on rocks
WO2016027390A1 (en) * 2014-08-21 2016-02-25 日本電気株式会社 Slope monitoring system, device for slope safety analysis, method, and program
CN104881583A (en) * 2015-06-04 2015-09-02 唐晓松 Multi-means, dynamic and whole-process landslide prewarning method
CN104964878A (en) * 2015-07-14 2015-10-07 中国科学院武汉岩土力学研究所 Triaxial test system and method for unsaturated soil multi-field coupling
CN105139585A (en) * 2015-07-29 2015-12-09 同济大学 Intelligent early warning and prediction method for soil slope danger
CN105912753A (en) * 2016-04-01 2016-08-31 国家海洋局第二海洋研究所 Method for analyzing three-dimensional stability of submarine slope based on strength reduction method
CN107169271A (en) * 2017-04-26 2017-09-15 中国地质大学(武汉) A kind of stability dynamic evaluation method of strain softening side slope
CN108538026A (en) * 2018-02-23 2018-09-14 青岛理工大学 Multi-index comprehensive evaluation and early warning method for slope stability
CN108733961A (en) * 2018-06-12 2018-11-02 河海大学 A kind of Slope Sliding Convenient stable criterion considering non-saturated region effect

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
Analytic hierarchy process for assessing factors influencing the stability of soil slopes reinforced with piles;Xu, Jiancong;《ENVIRONMENTAL EARTH SCIENCES》;20131031;第70卷(第4期);第1507-1517页 *
FLAC 3D 粘弹性模型的二次开发及其应用;徐平 等;《长江科学院院报》;20040430;第21卷(第2期);第10-13页 *
基于ABAQUS饱和渗流作用下岩质边坡稳定性分析;唐岩岩等;《煤炭技术》;20180310(第03期);第84-87页 *
岩石的非线性流变损伤模型及其应用研究;金丰年等;《解放军理工大学学报(自然科学版)》;20000625(第03期);第1-5页 *
降雨入渗引起非饱和土边坡破坏的水-土-气 三相渗流-变形耦合有限元分析;熊勇林 等;《岩土力学》;20170131;第38卷(第1期);第284-290页 *

Also Published As

Publication number Publication date
CN110261573A (en) 2019-09-20

Similar Documents

Publication Publication Date Title
CN110261573B (en) Dynamic evaluation method for stability of high-position rocky landslide
Puebla et al. Analysis of CANLEX liquefaction embankments: prototype and centrifuge models
CN107563014B (en) Method for calculating buckling strain and critical fault displacement of pipeline under fault action
Xu et al. A large triaxial investigation of the stress-path-dependent behavior of compacted rockfill
JP3692507B2 (en) Liquefaction prediction system
CN108287945B (en) Deformation calculation method and application technology of foundation soil under large foundation
Li et al. Experimental and numerical investigations on the shear behavior of a jointed rock mass
CN105424466A (en) Method for evaluating resistance of sandy soil on embedded submarine pipeline under action of waves
CN114004055A (en) Slope shear strength parameter inversion analysis method based on equivalent soil pressure effect
Huang et al. An advanced grain-based model to characterize mechanical behaviors of crystalline rocks with different weathering degrees
CN111090951A (en) Dynamic strength reduction method considering shear band unit strength softening
Liu et al. Cracking risk analysis of face slabs in concrete face rockfill dams during the operation period
CN108643246A (en) The prediction technique that artesian water precipitation influences surrounding enviroment in foundation pit construction
AU2021106123A4 (en) Method for Calculating the Ultimate Expansion Shear Force of Expansive Soil Pile Foundation Based on Finite Element Simulation
Li et al. State-of-the-Art Constitutive Modelling of Frozen Soils
Xie et al. Numerical study on cracking behavior and fracture failure mechanism of fractured rocks under shear loading
CN111189720A (en) Method for measuring stability safety coefficient of rocky slope
Lade Three-dimensional behaviour and parameter evaluation of an elastoplastic soil model
CN110348098B (en) Clay loess tunnel excavation simulation parameter assignment method
Sabzi et al. The performance of buildings adjacent to excavation supported by inclined struts
Mallick et al. Multivariate adaptive regression spline approach to the assessment of surface mean pressure coefficient on surfaces of C-shaped building
CHENG et al. Establishment of constitutive relation of shear deformation for irregular joints in sandstone
Wu et al. Numerical analysis of the effect of pile tip shape on soil behavior around pile
Freitas et al. Discontinuous hydromechanical modelling of concrete dam foundations
CN113297665A (en) Method for establishing soil-rock mixture constitutive model considering stone content

Legal Events

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