CN115422763A - Method for determining equivalent calculation parameters of landslide and surge under action of earthquake force - Google Patents
Method for determining equivalent calculation parameters of landslide and surge under action of earthquake force Download PDFInfo
- Publication number
- CN115422763A CN115422763A CN202211152602.XA CN202211152602A CN115422763A CN 115422763 A CN115422763 A CN 115422763A CN 202211152602 A CN202211152602 A CN 202211152602A CN 115422763 A CN115422763 A CN 115422763A
- Authority
- CN
- China
- Prior art keywords
- action
- soil layer
- strip
- landslide
- earthquake
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)
Abstract
The invention provides a method for determining equivalent calculation parameters of landslide and swell under the action of earthquake dynamic force, which is characterized in that based on the equivalent relation between landslide limit equilibrium states corresponding to the action of earthquake dynamic force and the action of no earthquake dynamic force, the corresponding relation from original parameters to equivalent parameters is established, so that the equivalent calculation parameters of landslide and swell under the action of earthquake dynamic force are determined; the method is a method for calculating landslide swell numerical values under the action of earthquake dynamic force by using the equivalent relation of the ultimate balance safety coefficient. The method has the advantages that the action of the earthquake dynamic force is equivalent to the weakening of the mechanical parameters, the flow is simple, the operability is strong, the problem that the action of the earthquake dynamic force is difficult to apply in the conventional landslide swell numerical simulation process can be effectively solved, and the swell disaster condition under the action of the earthquake dynamic force can be effectively prevented.
Description
Technical Field
The invention relates to the field of landslide surge numerical simulation, in particular to a method for determining equivalent calculation parameters of landslide surge under the action of a seismic force.
Background
The large landslide of the large landslide rushes into the reservoir at a high speed, the water body is excited to generate huge surge waves, the huge surge waves are excited to generate the water body and then impact the opposite bank, the water retaining building, the wharf and the like, and the surge waves can generate a climbing phenomenon to cross over the top of the dam and damage the dam and the downstream. How to correctly evaluate landslide and surge disasters induced by earthquakes so as to carry out targeted prevention and treatment has important scientific significance and engineering value.
The appropriate numerical simulation method can accurately simulate the generation and propagation rules of landslide surge, and the simulation result can provide guidance suggestions for landslide surge prevention and engineering disaster prevention and reduction. In the existing landslide swell numerical simulation software, the simulation of the action of seismic dynamic force is complicated, and part of the existing software does not support the simulation of the action of the seismic dynamic force and cannot effectively apply the action of the seismic dynamic force to the landslide, so that the swell disaster condition caused by correctly simulating the landslide under the action of the seismic dynamic force is seriously hindered.
Disclosure of Invention
Aiming at the technical problems, the invention provides the method for determining the equivalent calculation parameters of landslide and swell under the action of the earthquake force, which effectively applies the earthquake force to the landslide and can be applied to landslide and swell simulation.
The invention adopts the following technical scheme:
a method for determining equivalent calculation parameters of landslide and swell under the action of earthquake dynamic force is characterized by comprising the following steps of executing the following steps aiming at a target landslide area to obtain equivalent calculation parameters of landslide and swell of soil layers of various preset types in the target landslide area under the action of earthquake dynamic force, and using the equivalent calculation parameters for landslide and swell of the target landslide area under the action of earthquake dynamic force to carry out numerical simulation calculation:
step A: aiming at a target landslide area, acquiring cohesive force and internal friction angle of each preset type of soil layer in the target landslide area under no external force action;
and B: aiming at a target landslide area, vertically and equally dividing a soil layer of the target landslide area into a preset number of strip-shaped soil layer blocks along a horizontal direction corresponding to the highest point to the lowest point of a landslide, and establishing a landslide strip division model of the target landslide area;
and C: obtaining a limit balance safety coefficient of the target landslide area under the action of earthquake dynamic force based on a landslide strip model of the target landslide area and cohesive force and an internal friction angle of each preset type soil layer of the target landslide area under the action of no external force in combination with the earthquake acceleration of the preset target landslide area;
step D: and obtaining landslide swell equivalent calculation parameters of each preset type of soil layer in the target landslide area under the action of earthquake power by combining cohesive force and internal friction angle of each preset type of soil layer in the target landslide area under the action of no external force based on the ultimate balance safety coefficient of the target landslide area under the action of earthquake power.
Preferably, in the step C, the ultimate balance safety factor of the target landslide area under the action of the seismic power is obtained by the following steps:
step C1: based on a landslide strip sub-model of the target landslide area, combining with the earthquake acceleration of the preset target landslide area to obtain the horizontal earthquake load of each strip-shaped soil layer block of the target landslide area under the action of earthquake dynamic force;
and step C2: acquiring the vertical seismic load of each strip-shaped soil layer block under the action of seismic dynamics based on the horizontal seismic load of each strip-shaped soil layer block;
and C3: and based on the horizontal seismic load and the vertical seismic load of each strip-shaped soil layer block, combining a landslide strip model of the target landslide area and cohesive force and internal friction angle of each preset type soil layer of the target landslide area under the action of no external force, and obtaining the ultimate balance safety coefficient of the target landslide area under the action of seismic power.
Preferably, in the step C1, for each strip-shaped soil layer block, the horizontal seismic load of the strip-shaped soil layer block under the action of the seismic motion is obtained through the following formula, so as to obtain the horizontal seismic load of each strip-shaped soil layer block under the action of the seismic motion:
F hi ＝α w W i α i
in the formula, F hi Representing the horizontal seismic load of the strip-shaped soil layer block i; alpha is alpha w The synthetic horizontal seismic coefficients are represented by,
a h representing the seismic acceleration of a preset target landslide area; xi represents a preset reduction coefficient; g represents the gravitational acceleration; w i Representing the self-gravity of the strip-shaped soil layer block i; a is a i And representing the dynamic distribution coefficient of the strip-shaped soil layer block i.
Preferably, in the step C2, for each strip-shaped soil layer block, the vertical seismic load of the strip-shaped soil layer block under the action of the seismic dynamic force is obtained through the following formula, so as to obtain the vertical seismic load of each strip-shaped soil layer block under the action of the seismic dynamic force:
F vi ＝F hi /3
in the formula, F vi Representing the vertical seismic load of the strip-shaped soil layer block i; f hi And (4) representing the horizontal seismic load of the strip-shaped soil layer block i.
Preferably, in the step C3, the ultimate balance safety factor of the target landslide area under the action of the earthquake dynamic force is obtained through the following formula:
in the formula, c i Representing the cohesion of a preset type soil layer under the action of no external force;
representing the internal friction angle of a preset type soil layer under the action of no external force; w is a group of i Representing the self-gravity of the strip-shaped soil layer block i; n is a radical of i Representing the normal reaction force on the bottom surface of the strip-shaped soil layer block i; d i Representing the bottom slope angle of the strip-shaped soil layer block i; l i Representing the arc length of the bottom surface of the strip-shaped soil layer block i; f vi Representing the vertical seismic load of the strip-shaped soil layer block i under the action of seismic dynamic force; f hi Representing the horizontal seismic load of the strip-shaped soil layer block i under the action of seismic dynamic force; f s The ultimate balance safety coefficient of the target landslide area under the action of earthquake dynamic force is represented; i represents the total number of the preset strip-shaped soil layer blocks; s represents a preset type soil layer, and S represents a total number of preset type soil layers.
Preferably, in the step D, the equivalent calculation parameters of landslide swell of each preset type of soil layer in the target landslide area under the action of seismic dynamics are obtained through the following processes:
step D1: constructing an equivalent formula of the limit balance safety coefficient by the following formula:
of formula (II) F' s Equivalent parameters of the ultimate balance safety coefficient of the target landslide area under the action of earthquake dynamic force are represented; c' i Representing equivalent parameters of cohesive force of soil layers of preset types in a target landslide area under the action of earthquake dynamic force;
representing equivalent parameters of a preset type soil layer internal friction angle of a target landslide area under the action of seismic dynamic force; l i The arc length of the bottom surface of the strip-shaped soil layer block i is represented; w i Representing the self-gravity of the strip-shaped soil layer block i; n is a radical of hydrogen i Representing the normal reaction force on the bottom surface of the strip-shaped soil layer block i; d i Representing the bottom slope angle of the strip-shaped soil layer block i; i represents the total number of the preset strip-shaped soil layer blocks; s represents a preset type soil layer, and S represents the total number of the preset type soil layers;
step D2: presetting the target landslide area under no external forceIteratively weakening cohesive force and internal friction angle of type soil layers in equal proportion to obtain c' i And
c 'obtained at each iteration' i And
bringing into an equivalent formula to obtain F' s To F' s The ultimate balance safety coefficient is equal to the ultimate balance safety coefficient of the target landslide area under the action of earthquake dynamic force, iteration is finished, and F 'is output' s C 'respectively corresponding to all preset type soil layers when the ultimate balance safety factor of the target landslide area under the action of earthquake power is equal' i And
and the equivalent calculation parameters are used as landslide and swell equivalent calculation parameters of the preset soil layer in the target landslide area under the action of earthquake dynamic force.
The invention has the beneficial effects that: the invention provides a method for determining equivalent calculation parameters of landslide and surge under the action of earthquake dynamic force, which is characterized in that based on the equivalent relationship between landslide limit equilibrium states corresponding to the action of earthquake dynamic force and the action of no earthquake dynamic force, the corresponding relationship from original parameters to equivalent parameters is established, so that equivalent calculation parameters of landslide and surge under the action of earthquake dynamic force are determined; the method is a method for calculating landslide swell numerical values under the action of earthquake dynamic force by using the equivalent relation of the ultimate balance safety coefficient. The method has the advantages that the earthquake dynamic force effect is equivalent to weakening of mechanical parameters, the flow is simple, the operability is strong, the problem that the earthquake dynamic force effect is difficult to apply in the conventional landslide swell numerical simulation process can be effectively solved, and the swell disaster condition under the earthquake dynamic force effect can be effectively prevented based on the landslide swell numerical simulation.
Drawings
FIG. 1 is a flow chart of the method of the present invention;
FIG. 2 is a schematic view of a landslide section modeling of the present invention;
FIG. 3 is a schematic diagram of the stress of the strip-shaped soil horizon blocks under the action of earthquake dynamic force;
FIG. 4 is a schematic diagram of the stress of the strip soil blocks under the equivalent parameters of the invention;
fig. 5 is a model diagram of a landslide section in the present embodiment.
Detailed Description
The invention is further described below with reference to the accompanying drawings. The following examples will give the skilled person a more complete understanding of the present invention, but do not limit the invention in any way.
A method for determining equivalent calculation parameters of landslide and swell under the action of earthquake and vibration force is disclosed, wherein a flow chart is shown in figure 1, aiming at a target landslide area, the following steps are executed to obtain equivalent calculation parameters of landslide and swell of soil layers of various preset types in the target landslide area under the action of earthquake and vibration force, and the equivalent calculation parameters are used for simulating and calculating landslide and swell numerical values in the target landslide area under the action of earthquake and vibration force:
step A: and acquiring a landslide related geological profile and parameters aiming at the target landslide area. In this embodiment, a landslide in a reservoir area of a hydropower station in yamo river is selected, and the relevant geological profile and parameters mainly include landslide profile information, soil density, soil cohesion and soil internal friction angle; acquiring cohesive force and internal friction angle of each preset type of soil layer in a target landslide area under no external force action based on landslide related geological profile and parameters; table 1 is a landslide-related geological parameter table, and five soil layers shown in the following table are selected in this embodiment.
TABLE 1 landslide soil layer parameter table
| Grouping of materials | Cohesion (kPa) | Internal friction angle (°) | Severe (kg/m) 3 ) |
| Soil containing gravel and powder | 90 | 31.0 | 2250 |
| Deep sliding belt | 60 | 31.5 | 2250 |
| Collapse slope layer | 45 | 33.0 | 2250 |
| Crushed stone mixed soil | 70 | 35.0 | 2300 |
| Rock layer | 100 | 40.0 | 2400 |
And B, step B: aiming at a target landslide area, as shown in fig. 2, vertically equally dividing a soil layer of the target landslide area into a preset number of strip-shaped soil layer blocks along a horizontal direction corresponding to the direction from the highest point to the lowest point, and establishing a landslide strip division model of the target landslide area; in this embodiment, a landslide section model is constructed as shown in fig. 5.
And C: obtaining a limit balance safety coefficient of the target landslide area under the action of earthquake dynamic force based on a landslide strip model of the target landslide area and cohesive force and an internal friction angle of each preset type soil layer of the target landslide area under the action of no external force in combination with the earthquake acceleration of the preset target landslide area; the seismic acceleration is determined according to the existing specifications and the on-site defense intensity, and according to the regulations in Chinese seismic motion parameter zone diagrams (GB 18306-2015), the seismic acceleration of the preset target landslide area in the embodiment is obtained by inquiry and is 0.15g. The earthquake acceleration of the target landslide area is preset to be a value based on a standard, namely the earthquake maximum intensity possibly occurring in the area is determined, so that the damage caused by landslide and swell under the action of earthquake dynamic force can be prevented well in advance.
In the step C, the ultimate balance safety coefficient of the target landslide area under the action of earthquake dynamic force is obtained by the following steps:
step C1: based on a landslide strip sub-model of the target landslide area, combining with the earthquake acceleration of a preset target landslide area to obtain horizontal earthquake loads of strip-shaped soil layer blocks in the target landslide area under the action of earthquake dynamic force;
in the step C1, for each strip-shaped soil layer block, the horizontal seismic load of the strip-shaped soil layer block under the action of seismic dynamics is obtained through the following formula, and further the horizontal seismic load of each strip-shaped soil layer block under the action of seismic dynamics is obtained:
F hi ＝α w W i α i
in the formula, F hi Representing the horizontal seismic load of the strip-shaped soil layer block i; alpha is alpha w A comprehensive horizontal seismic coefficient is represented by the total horizontal seismic coefficient,
α h representing the seismic acceleration of a preset target landslide area; xi represents a preset reduction coefficient, and 0.25 is taken; g represents the gravitational acceleration; w i Representing the self-gravity of the strip-shaped soil layer block i; alpha (alpha) ("alpha") i Representing the dynamic distribution coefficient, alpha, of the strip-shaped soil layer block i i ∈[1，3]。
And step C2: acquiring the vertical seismic load of each strip-shaped soil layer block under the action of seismic dynamics based on the horizontal seismic load of each strip-shaped soil layer block;
in the step C2, the vertical seismic load of each strip-shaped soil layer block under the action of seismic dynamics is obtained by the following formula for each strip-shaped soil layer block, and then the vertical seismic load of each strip-shaped soil layer block under the action of seismic dynamics is obtained:
F vi ＝F hi /3
in the formula, F vi Representing the vertical seismic load of the strip-shaped soil layer block i; f hi Representing the horizontal seismic load of the strip-shaped soil layer block i. The vibration force acting load of the strip soil block is shown in table 2.
TABLE 2 earthquake dynamic acting load borne by strip soil layer block
| Number of bar | Horizontal seismic load F hi (kN) | Vertical seismic load F vi (kN) |
| 1 | 2.37 | 0.79 |
| 2 | 39.44 | 13.15 |
| 3 | 38.50 | 12.83 |
| ... | ... | ... |
| 88 | 41.47 | 13.82 |
| 89 | 10.89 | 3.63 |
| 90 | 1.07 | 0.36 |
FIG. 3 is a schematic view showing the stress of the strip-shaped soil horizon block under the action of earthquake dynamic force, wherein E i+1 Is the horizontal force of the previous strip-shaped soil layer block to the strip-shaped soil layer block i, E i The horizontal acting force, X, of the next strip-shaped soil layer block to the strip-shaped soil layer block i i+1 Is the friction force, X, of the previous strip-shaped soil layer block to the strip-shaped soil layer block i i Is the friction of the latter strip-shaped soil layer block to the strip-shaped soil layer block i. In the invention, the horizontal acting forces on two sides of the strip-shaped soil block are assumed to be equal in magnitude, and the friction forces are assumed to be equal in magnitude, so that E and X are not considered in calculation. T is a unit of i The frictional force applied to the strip-shaped soil layer block i
Step C3: and based on the horizontal seismic load and the vertical seismic load of each strip-shaped soil layer block, combining the landslide strip sub-model of the target landslide area and the cohesive force and the internal friction angle of each preset type soil layer of the target landslide area under the action of no external force to obtain the ultimate balance safety coefficient of the target landslide area under the action of seismic power.
In the step C3, the ultimate balance safety coefficient of the target landslide area under the action of earthquake dynamic force is obtained through the following formula:
in the formula, c i Representing the cohesive force of a preset type soil layer under the action of no external force;
representing the internal friction angle of a preset type soil layer under the action of no external force; w i Representing the self-gravity of the strip-shaped soil layer block i; n is a radical of i Representing the normal reaction force on the bottom surface of the strip-shaped soil layer block i; d is a radical of i Representing the bottom slope angle of the strip-shaped soil layer block i; l i Representing the arc length of the bottom surface of the strip-shaped soil layer block i; f vi Representing the vertical seismic load of the strip-shaped soil layer block i under the action of seismic dynamic force; f hi Representing the horizontal seismic load of the strip-shaped soil layer block i under the action of seismic dynamic force; f s The ultimate balance safety coefficient of the preset soil layer in the target landslide area under the action of earthquake dynamic force is represented; i represents the total number of the preset strip-shaped soil layer blocks; s represents a preset type soil layer, and S represents a total number of preset type soil layers. Based on the horizontal seismic load and the vertical seismic load of each strip-shaped soil layer block and the cohesive force and the internal friction angle of each preset type soil layer in the target landslide area without external force, summing calculation is carried out based on the formula to obtain the ultimate balance safety coefficient of the target landslide area under the action of seismic dynamic force, wherein the ultimate balance safety coefficient is F s ＝0.983。
Step D: and obtaining landslide swell equivalent calculation parameters of each preset type of soil layer in the target landslide area under the action of earthquake power by combining cohesive force and internal friction angle of each preset type of soil layer in the target landslide area under the action of no external force based on the ultimate balance safety coefficient of the target landslide area under the action of earthquake power.
In the step D, obtaining equivalent calculation parameters of landslide swell of soil layers of each preset type in the target landslide area under the action of earthquake dynamic force by the following processes:
step D1: constructing an equivalent formula of the ultimate balance safety coefficient through the following formula, namely the ultimate balance safety coefficient under equivalent parameters is expressed as:
of formula (II) F' s Equivalent parameters of the ultimate balance safety coefficient of the target landslide area under the action of earthquake dynamic force are represented; c' i Representing equivalent parameters of cohesive force of a preset type soil layer in a target landslide area under the action of earthquake dynamic force;
representing equivalent parameters of a preset type soil layer internal friction angle of a target landslide area under the action of seismic dynamic force; l i Representing the arc length of the bottom surface of the strip-shaped soil layer block i; w i Representing the self-gravity of the strip-shaped soil layer block i; n is a radical of i Representing the normal reaction force on the bottom surface of the strip-shaped soil layer block i; d i Representing the bottom slope angle of the strip-shaped soil layer block i; i represents the total number of the preset strip-shaped soil layer blocks; s represents a preset type soil layer, and S represents the total number of the preset type soil layers;
step D2: c 'is obtained by iteratively weakening the cohesive force and the internal friction angle of each preset type soil layer in the target landslide area in equal proportion under the action of no external force' i And
c 'obtained at each iteration' i And
substituting equivalent formula to obtain F' s To F' s The ultimate balance safety coefficient is equal to the ultimate balance safety coefficient of the target landslide area under the action of earthquake dynamic force, iteration is finished, and F 'is output' s C 'corresponding to each preset type soil layer when limit balance safety factor of target landslide area under earthquake power action is equal' i And
and the equivalent calculation parameters are used as landslide and swell equivalent calculation parameters of the preset soil layer in the target landslide area under the action of earthquake dynamic force.
Specifically, mechanical parameters (cohesive force and internal friction angle) of each type of soil layer are weakened in equal proportion to obtain c' i And
substitution into
Weakening and calculating for multiple times, and finally finding out the condition F s ＝F′ s ，
The weakening parameters are equivalent parameters under the action of seismic dynamic force of the landslide soil layer. The equivalent parameters are shown in table 3 and shown in figure 4 as a stress schematic diagram of strip-shaped soil horizon blocks under the equivalent parameters of the invention, T' i Is the friction force applied to the strip-shaped soil layer block i under the equivalent parameter
TABLE 3 equivalent parameter table under action of earthquake dynamic force of landslide soil layer
| Grouping of materials | Equivalent cohesion force (kPa) | Equivalent internal friction angle (°) | Severe (kg/m) 3 ) |
| Soil containing gravel and powder | 78.3 | 27.0 | 2250 |
| Deep sliding belt | 52.2 | 27.4 | 2250 |
| Collapse slope layer | 39.2 | 28.7 | 2250 |
| Crushed stone mixed soil | 60.9 | 30.5 | 2300 |
| Lump stone layer | 87.0 | 34.8 | 2400 |
The ultimate balance safety factor calculated using the equivalence formula is made equal to the ultimate balance safety factor calculated using the original parameters, i.e. F s ＝F′ s And obtaining the calculation parameters equivalent to the action of the earthquake dynamic force from the equivalence relation. The obtained equivalent parameters can be used for simulating and researching landslide surge numerical values under the action of earthquake dynamic force, and landslide surge disasters induced by earthquakes can be effectively evaluated based on the landslide surge numerical simulation utilizing the equivalent parameters, so that targeted prevention and treatment are carried out.
The scheme designs a method for calculating landslide swell numerical calculation parameters under the action of earthquake dynamic force by utilizing the equivalent relation of the ultimate balance safety coefficient, and the action of the earthquake dynamic force is equivalent to the weakening of mechanical parameters. And establishing a corresponding relation from the original parameters to the equivalent parameters by utilizing the equivalent relation between the landslide limit equilibrium states corresponding to the action of the seismic dynamic force and the action of the seismic dynamic force, thereby determining equivalent calculation parameters of landslide swell under the action of the seismic dynamic force and equating the action of the seismic dynamic force as the weakening of the landslide parameters. The method is simple in flow and strong in operability, and can effectively solve the problem that the earthquake dynamic force is difficult to apply in the conventional landslide swell numerical simulation process, so that the swell disaster condition under the earthquake dynamic force can be effectively prevented based on the landslide swell numerical simulation.
Although the present invention has been described in detail with reference to the foregoing embodiments, it will be apparent to those skilled in the art that modifications may be made to the embodiments described in the foregoing detailed description, or equivalent arrangements may be substituted for those skilled in the art. All equivalent structures made by using the contents of the specification and the attached drawings of the invention can be directly or indirectly applied to other related technical fields, and all the equivalent structures are within the protection scope of the invention.
Claims (6)
1. A method for determining equivalent calculation parameters of landslide and surge under the action of earthquake force is characterized by comprising the following steps: aiming at the target landslide area, the following steps are executed to obtain landslide swell equivalent calculation parameters of soil layers of all preset types in the target landslide area under the action of earthquake dynamic force, and the landslide swell equivalent calculation parameters are used for simulating and calculating landslide swell numerical values of the target landslide area under the action of earthquake dynamic force:
step A: aiming at a target landslide area, acquiring cohesive force and internal friction angle of each preset type of soil layer in the target landslide area under no external force action;
and B, step B: aiming at a target landslide area, vertically and equally dividing a soil layer of the target landslide area into a preset number of strip-shaped soil layer blocks along a horizontal direction corresponding to the highest point to the lowest point of a landslide, and establishing a landslide strip division model of the target landslide area;
step C: obtaining a limit balance safety coefficient of the target landslide area under the action of earthquake dynamic force based on a landslide strip model of the target landslide area and cohesive force and an internal friction angle of each preset type soil layer of the target landslide area under the action of no external force in combination with the earthquake acceleration of the preset target landslide area;
step D: and obtaining landslide swell equivalent calculation parameters of each preset type of soil layer in the target landslide area under the action of earthquake power by combining cohesive force and internal friction angle of each preset type of soil layer in the target landslide area under the action of no external force based on the ultimate balance safety coefficient of the target landslide area under the action of earthquake power.
2. The method for determining equivalent calculation parameters of landslide and swell under the action of earthquake force as claimed in claim 1, wherein: in the step C, the ultimate balance safety coefficient of the target landslide area under the action of earthquake dynamic force is obtained by the following steps:
step C1: based on a landslide strip sub-model of the target landslide area, combining with the earthquake acceleration of a preset target landslide area to obtain horizontal earthquake loads of strip-shaped soil layer blocks in the target landslide area under the action of earthquake dynamic force;
and step C2: acquiring the vertical seismic load of each strip-shaped soil layer block under the action of seismic dynamic force based on the horizontal seismic load of each strip-shaped soil layer block;
and C3: and based on the horizontal seismic load and the vertical seismic load of each strip-shaped soil layer block, combining the landslide strip sub-model of the target landslide area and the cohesive force and the internal friction angle of each preset type soil layer of the target landslide area under the action of no external force to obtain the ultimate balance safety coefficient of the target landslide area under the action of seismic power.
3. The method for determining equivalent calculation parameters of landslide and swell under the action of earthquake force as claimed in claim 2, wherein: in the step C1, for each strip-shaped soil layer block, the horizontal seismic load of the strip-shaped soil layer block under the action of seismic dynamics is obtained through the following formula, and further the horizontal seismic load of each strip-shaped soil layer block under the action of seismic dynamics is obtained:
F hi ＝α w W i a i
in the formula, F hi Representing the horizontal seismic load of the strip-shaped soil layer block i; alpha is alpha w The synthetic horizontal seismic coefficients are represented by,
a h representing the seismic acceleration of a preset target landslide area; xi represents a preset reduction coefficient; g represents the acceleration of gravity; w is a group of i Representing the self-gravity of the strip-shaped soil layer block i; a is a i Representing the dynamic distribution coefficient of the strip-shaped soil layer block i.
4. The method for determining equivalent calculation parameters of landslide and swell under the action of earthquake force as claimed in claim 2, wherein: in the step C2, the vertical seismic load of each strip-shaped soil layer block under the action of seismic dynamics is obtained by the following formula for each strip-shaped soil layer block, and then the vertical seismic load of each strip-shaped soil layer block under the action of seismic dynamics is obtained:
F vi ＝F hi /3
in the formula, F vi Representing the vertical seismic load of the strip-shaped soil layer block i; f hi Representing the horizontal seismic load of the strip-shaped soil layer block i.
5. The method for determining equivalent calculation parameters of landslide and swell under the action of earthquake force as claimed in claim 2, wherein: in the step C3, the ultimate balance safety coefficient of the target landslide area under the action of earthquake dynamic force is obtained through the following formula:
in the formula, c i Representing the cohesive force of a preset type soil layer under the action of no external force;
representing the internal friction angle of a preset type soil layer under the action of no external force; w i Representing the self-gravity of the strip-shaped soil layer block i; n is a radical of i Representing the normal reaction force on the bottom surface of the strip-shaped soil layer block i; d i Representing the bottom slope angle of the strip-shaped soil layer block i; l i Representing the arc length of the bottom surface of the strip-shaped soil layer block i; f vi Representing the vertical seismic load of the strip-shaped soil layer block i under the action of seismic dynamic force; f hi Representing the horizontal seismic load of the strip-shaped soil layer block i under the action of seismic dynamic force; f s Representing the ultimate balance safety coefficient of the target landslide area under the action of earthquake dynamic force; i represents the total number of the preset strip-shaped soil layer blocks; s represents a preset type soil layer, and S represents the total number of the preset type soil layers.
6. The method for determining equivalent calculation parameters of landslide and swell under the action of earthquake force as claimed in claim 1, wherein: in the step D, obtaining equivalent calculation parameters of landslide swell of soil layers of each preset type in the target landslide area under the action of earthquake dynamic force by the following processes:
step D1: constructing an equivalent formula of the limit balance safety coefficient by the following formula:
of formula (II) F' s Equivalent parameters representing the ultimate balance safety coefficient of the target landslide area under the action of earthquake dynamic force; c' i Representing equivalent parameters of cohesive force of a preset type soil layer in a target landslide area under the action of earthquake dynamic force;
representing equivalent parameters of a preset type soil layer internal friction angle of a target landslide area under the action of seismic dynamic force; l. the i Representing the arc length of the bottom surface of the strip-shaped soil layer block i; w i Representing the self-gravity of the strip-shaped soil layer block i; n is a radical of i Representing the normal reaction force on the bottom surface of the strip-shaped soil layer block i; d i Bottom slope angle representing bar-shaped soil layer block i(ii) a I represents the total number of the preset strip-shaped soil layer blocks; s represents a preset type soil layer, and S represents the total number of the preset type soil layers;
step D2: c 'is obtained by iteratively weakening the cohesive force and the internal friction angle of each preset type soil layer in the target landslide area in equal proportion under the action of no external force' i And
c 'obtained at each iteration' i And
bringing into an equivalent formula to obtain F' s To F' s Equal to the limit balance safety coefficient of the target landslide area under the action of earthquake power, finishing the iteration and outputting F' s C 'corresponding to each preset type soil layer when limit balance safety factor of target landslide area under earthquake power action is equal' i And
and the equivalent calculation parameters are used as landslide and swell equivalent calculation parameters of the preset soil layer in the target landslide area under the action of earthquake dynamic force.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211152602.XA CN115422763A (en) | 2022-09-21 | 2022-09-21 | Method for determining equivalent calculation parameters of landslide and surge under action of earthquake force |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211152602.XA CN115422763A (en) | 2022-09-21 | 2022-09-21 | Method for determining equivalent calculation parameters of landslide and surge under action of earthquake force |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115422763A true CN115422763A (en) | 2022-12-02 |
Family
ID=84203662
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211152602.XA Pending CN115422763A (en) | 2022-09-21 | 2022-09-21 | Method for determining equivalent calculation parameters of landslide and surge under action of earthquake force |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115422763A (en) |
-
2022
- 2022-09-21 CN CN202211152602.XA patent/CN115422763A/en active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107122536B (en) | Interval tunnel anti-seismic numerical simulation method | |
| CN103559383B (en) | Method for predicting and evaluating strong earthquake resistance of nuclear power station breakwater | |
| Raja et al. | Behaviour of earth dam under seismic load considering nonlinearity of the soil | |
| Zhou et al. | Understanding the failure mechanism of a large underground cavern in steeply dipping layered rock mass using an enhanced ubiquitous-joint model | |
| Zacchei et al. | Seismic hazard and structural analysis of the concrete arch dam (Rules Dam on Guadalfeo River) | |
| CN115422763A (en) | Method for determining equivalent calculation parameters of landslide and surge under action of earthquake force | |
| Sawamura et al. | Experimental study on seismic resistance of a two-hinge precast arch culvert using strong earthquake response simulator | |
| Azizan et al. | Incremental dynamic analysis of koyna dam under repeated ground motions | |
| Prassetyo | Hydro-mechanical analysis of tunneling in saturated ground using a novel and efficient sequential coupling technique | |
| Yu et al. | Seismic finite element analysis of underground integrated pipe corridor based on generalized reaction displacement method | |
| CN112926118A (en) | Transformer substation engineering deep foundation soil-structure cooperative analysis method | |
| Li et al. | Modeling of soil–rock mixture landslides with the generalized interpolation material point method | |
| Wang et al. | The mechanical response of piles with consideration of pile-soil interactions under a periodic wave pressure | |
| Ilinca et al. | Influence of dynamic analysis methods on seismic response of a buttress dam | |
| Zhou et al. | Seismic analysis for nuclear power safety related bridge | |
| Dang et al. | Dynamic response analysis of intake tower in hydroelectric power station with high surrounding rock | |
| Zhou et al. | Influence of seismic wave type and incident direction on the dynamic response of tall concrete-faced rockfill dams | |
| Soroush | A comparison between seismic behaviors of earth dams with inclined and vertical clay cores-a numerical analysis approach | |
| Sawicki et al. | Wave-induced stresses and pore pressures near a mudline | |
| Mori et al. | Liquefaction analysis of Yodo-gawa River dike with discrete element method | |
| Mehdizadeh et al. | Investigation of effectiveness of deep mixing method on seismic response of liquefiable soils using 3D numerical simulation | |
| Zhong | Determination of bearing capacity of a sand-mudstone interbedded foundation in Chengdu | |
| Mirzaei et al. | Seismic analysis of double curved arch dams based performance | |
| CN113177252A (en) | Dam body geomembrane three-dimensional stress deformation calculation method | |
| Yaoru et al. | Stability analysis considering water effect for reservoir bank slope based on deformation reinforcement theory |
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 |