Landslide thrust calculation method and device, electronic equipment and medium

Abstract

The invention relates to a landslide thrust calculation method and device, electronic equipment and a medium, and belongs to the technical field of geotechnical engineering.The method comprises the steps that soil layer structure parameters of a side slope are obtained, and a rock-soil body model of the side slope is constructed; calculating soil layer structure parameters to obtain a theoretical value of the slope landslide thrust, performing simulation calculation on the rock-soil body model to obtain a simulation value, and taking a ratio of the simulation value to the theoretical value as a correction coefficient; based on different soil layer structure parameter values, linear relations between the structure parameters and theoretical values, simulation values and correction coefficients are obtained; determining a three-dimensional relationship between the soil layer structure parameter and the correction coefficient based on the linear relationship; acquiring soil layer structure parameters and theoretical values of the target slope; determining a target correction coefficient of the target slope; and determining the corrected landslide thrust of the target slope. According to the method, the correction coefficient is obtained through calculation of a large amount of data, so that the landslide thrust of the target slope in the three-dimensional state is obtained through the correction coefficient when the theoretical value of the target slope is known.

Description

Technical Field

[0001] The invention relates to the technical field of geotechnical engineering, and in particular to a landslide thrust calculation method, device, electronic equipment and medium. Background Art

[0002] With the gradual improvement of transportation infrastructure, the highway and railway networks in major areas of my country have basically taken shape, and the current construction focus has shifted to remote mountainous areas. Building roads in hilly and mountainous areas will inevitably pass through high embankments, deep cuttings and sections with complex geology. In addition to the huge amount of engineering, slope management is also important and complex. When the slope is excavated or the roadbed is filled during the construction process, the balance and integrity of the original slope will be broken, thereby inducing the possibility of slope instability and landslides, and creating construction safety hazards. In addition, design defects or improper construction will also endanger driving safety during the use stage and pose a threat to people's production and life. Therefore, slope stability and the remediation and development of landslides in mountainous areas have attracted attention from all walks of life.

[0003] In terms of the stability and force of slope support structures, domestic and foreign scholars have conducted a lot of research and obtained a series of valuable results. Cheng Aiping et al. first used the finite element strength reduction method to analyze the stability of the slope, and then used the transfer coefficient method to obtain the corresponding landslide thrust, and proposed a double-row pile reinforcement scheme; Xu Hong et al., Shi Jian et al. calculated the slope stability under different test level combinations based on orthogonal test analysis, and studied the sensitivity of various influencing factors such as slope height, slope, rock and soil weight, cohesion, internal friction angle, etc. to the slope safety factor, and ranked them; Li Mei et al. combined theoretical calculation, model test, and finite element simulation methods to determine the landslide thrust distribution diagram behind the pile, providing reference for the design and optimization of anti-slide piles.

[0004] From the above research, it can be seen that the existing research mainly focuses on the stability and mechanical analysis of slopes under normal conditions, and most of them are single pile support structures. There are not many studies on the support effect of double-row anti-slide piles, especially the interaction between piles and landslides. Moreover, the existing specifications have only made some principled provisions for the design of anti-slide pile support structures, and the design of double-row pile support structures involves less content. The calculation method of landslide thrust in the existing highway and railway specifications uses the transfer coefficient method (explicit or implicit) for the model mainly two-dimensional model, without considering the spatial effect between the landslide body and the anti-slide pile. Therefore, the research on the stability and landslide thrust of double-row anti-slide pile support of the accumulation layer slope under three-dimensional state has important theoretical significance and engineering practice value.

[0005] In summary, the prior art lacks a method for calculating the landslide thrust of a double-row anti-slide pile slope in an accumulation layer in a three-dimensional state. Summary of the invention

[0006] In view of this, it is necessary to provide a landslide thrust calculation method, device, electronic equipment and medium to achieve the purpose of calculating the landslide thrust of a double-row anti-slide pile slope of an accumulation layer slope in a three-dimensional state.

[0007] In order to achieve the above object, the present invention provides a landslide thrust calculation method, comprising: Obtain the soil structure parameters of the slope and construct the rock and soil model of the slope; The soil layer structure parameters are calculated to obtain the theoretical value of the landslide thrust of the slope, and the rock and soil model is simulated to obtain the simulated value of the landslide thrust of the slope, and the ratio of the simulated value to the theoretical value is used as the correction coefficient; Based on the structural parameter values ​​of different soil layers, the linear relationship between the structural parameters and the theoretical values, the simulated values ​​and the correction coefficients is obtained; Determine the three-dimensional relationship between soil structure parameters and correction coefficients based on linear relationships; Obtain target soil layer structural parameters and target theoretical values ​​of the target slope; Determine the target correction coefficient based on the target soil layer structural parameters and three-dimensional relationship; The corrected landslide thrust of the target slope is determined based on the target theoretical value and the target correction coefficient.

[0008] In a possible implementation, the soil layer structural parameters include: one or more of sliding body volume, sliding body inclination, cohesion, and internal friction angle. In a possible implementation, the soil layer structural parameters are calculated to obtain a theoretical value of the landslide thrust of the slope, and a rock and soil model is simulated to obtain a simulated value of the landslide thrust of the slope, including: The theoretical value of the landslide thrust is obtained by calculating the soil structure parameters based on the transfer coefficient method. The rock and soil model is simulated and calculated based on finite element software to obtain the simulated value of the slope landslide thrust.

[0009] In a possible implementation, the method of calculating the soil layer structure parameters based on the transfer coefficient method to obtain the theoretical value of the landslide thrust of the slope includes: When the slope is a railway slope, the rock and soil model is calculated based on the explicit solution of the transfer coefficient method to obtain the theoretical value of the landslide thrust of the slope; When the slope is a highway slope, the rock and soil model is calculated based on the implicit solution of the transfer coefficient method to obtain the theoretical value of the landslide thrust of the slope.

[0010] In a possible implementation, the linear relationship between the structural parameters and the theoretical values, the simulation values ​​and the correction coefficients obtained based on the structural parameter values ​​of different soil layers includes: When the value of the slider volume changes, a first linear relationship between the slider volume and the theoretical value, the simulation value, and the correction coefficient is obtained; When the value of the sliding body inclination angle changes, a second linear relationship between the sliding body inclination angle and the theoretical value, the simulation value and the correction coefficient is obtained; When the value of the cohesion changes, a third linear relationship between the cohesion and the theoretical value, the simulation value and the correction coefficient is obtained; When the value of the internal friction angle changes, a fourth linear relationship between the internal friction angle and the theoretical value, the simulation value, and the correction coefficient is obtained.

[0011] In a possible implementation, determining the three-dimensional relationship between the soil layer structure parameter and the correction coefficient based on the linear relationship includes: Determine a three-dimensional relationship between the sliding body volume, the sliding body inclination and the correction coefficient based on the first linear relationship, the second linear relationship, the third linear relationship and the fourth linear relationship; A three-dimensional relationship among cohesion, internal friction angle and correction coefficient is determined based on the first linear relationship, the second linear relationship, the third linear relationship and the fourth linear relationship.

[0012] In a possible implementation, determining the corrected landslide thrust of the target slope based on the target theoretical value and the target correction coefficient includes: The target theoretical value is multiplied by the target correction coefficient to obtain the corrected landslide thrust of the target slope. The corrected landslide thrust of the target slope is the landslide thrust in the three-dimensional scene.

[0013] On the other hand, the present invention also provides a landslide thrust calculation device, comprising: A model building module is used to obtain the soil structure parameters of the slope and build a rock and soil model of the slope; The target data acquisition module is used to calculate the soil layer structure parameters to obtain the theoretical value of the slope landslide thrust, simulate the rock and soil model to obtain the simulated value of the slope landslide thrust, and use the ratio of the simulated value to the theoretical value as the correction coefficient; A linear relationship determination module is used to obtain the linear relationship between the structural parameters and the theoretical values, simulation values ​​and correction coefficients based on the structural parameter values ​​of different soil layers; A three-dimensional relationship determination module, used for determining the three-dimensional relationship between soil layer structure parameters and correction coefficients based on a linear relationship; A target theoretical value acquisition module is used to obtain target soil layer structural parameters and target theoretical values ​​of a target slope; A target correction coefficient acquisition module is used to determine the target correction coefficient based on the target soil layer structure parameters and the three-dimensional relationship; The target landslide thrust acquisition module is used to determine the corrected landslide thrust of the target slope based on the target theoretical value and the target correction coefficient.

[0014] On the other hand, the present invention also provides an electronic device, including a memory and a processor, wherein: The memory is used to store programs; The processor is coupled to the memory and is used to execute the program stored in the memory to implement the steps in a landslide thrust calculation method described in any one of the above implementations.

[0015] On the other hand, the present invention further provides a computer-readable storage medium for storing a computer-readable program or instruction, which, when executed by a processor, can implement the steps in a landslide thrust calculation method described in any of the above implementations.

[0016] The beneficial effects of the present invention are as follows: a landslide thrust calculation method provided by the present invention first obtains the soil layer structure parameters of the slope, and constructs a rock and soil model of the slope, then calculates the soil layer structure parameters to obtain the theoretical value of the slope landslide thrust, performs simulation calculation on the rock and soil model to obtain the simulated value of the slope landslide thrust, takes the ratio of the simulated value to the theoretical value as the correction coefficient, obtains the linear relationship between the structure parameter and the theoretical value, the simulated value and the correction coefficient based on different soil layer structure parameter values, thereby determining the soil layer structure parameters affecting the landslide thrust according to the linear relationship, further determining the three-dimensional relationship between the soil layer structure parameters and the correction coefficient based on the linear relationship, obtaining the target soil layer structure parameters and the target theoretical value of the target slope, determining the target correction coefficient based on the target soil layer structure parameters and the three-dimensional relationship, and determining the corrected landslide thrust of the target slope based on the target theoretical value and the target correction coefficient. The present invention obtains the correction coefficient by calculating a large amount of data, so that when the theoretical value of the target slope is known, the landslide thrust of the target slope in the three-dimensional scene is obtained by the correction coefficient. BRIEF DESCRIPTION OF THE DRAWINGS

[0017] Figure 1 A flow chart of an embodiment of a method for calculating landslide thrust provided by the present invention; Figure 2 A slope grid division diagram in an embodiment of a landslide thrust calculation method provided by the present invention; Figure 3 For the present invention Figure 1 A schematic flow chart of an embodiment of S102; Figure 4 is a landslide thrust distribution diagram in an embodiment of the present invention; Figure 5 is a relationship diagram between the sliding body volume and the landslide thrust in an embodiment of the present invention; Figure 6 is a relationship diagram between the sliding body inclination angle and the landslide thrust in an embodiment of the present invention; Figure 7 is a relationship diagram between the shear strength index cohesion of the sliding body and the landslide thrust in the embodiment of the present invention; Figure 8 is a relationship diagram between the internal friction angle, which is a shear strength index of the sliding body, and the landslide thrust in an embodiment of the present invention; Fig. 9 For the present invention Figure 1 A schematic flow chart of an embodiment of S104; Fig.10 A three-dimensional surface diagram showing the relationship between the correction coefficient of the theoretical formula and the volume and inclination angle of the sliding body in the embodiment of the present invention; Fig.11 A three-dimensional surface diagram showing the relationship between the correction coefficient of the theoretical formula and the cohesion and internal friction angle in an embodiment of the present invention; Fig.12 A comparison diagram of the simulated value and the measured value of the landslide thrust distribution in the verification of the embodiment of the present invention; Fig.13 A schematic flow chart of an embodiment of a landslide thrust calculation device provided by the present invention; Fig.14 A schematic structural diagram of an embodiment of an electronic device provided by the present invention. DETAILED DESCRIPTION

[0018] The preferred embodiments of the present invention are described in detail below in conjunction with the accompanying drawings, wherein the accompanying drawings constitute a part of this application and are used together with the embodiments of the present invention to illustrate the principles of the present invention, but are not used to limit the scope of the present invention.

[0019] The present invention provides a landslide thrust calculation method, device, electronic equipment and medium, which are respectively described below.

[0020] Figure 1 A schematic flow chart of an embodiment of a method for calculating landslide thrust provided by the present invention is shown in FIG. Figure 1 As shown in FIG. 1 , the landslide thrust calculation method includes: S101, obtaining soil layer structure parameters of the slope and constructing a rock and soil model of the slope; S102, calculating soil layer structural parameters to obtain theoretical values ​​of the thrust of the slope landslide, performing simulation calculation on the rock and soil model to obtain a simulated value of the thrust of the slope landslide, and taking the ratio of the simulated value to the theoretical value as a correction coefficient; S103, obtaining a linear relationship between the structural parameter and the theoretical value, the simulation value and the correction coefficient based on the structural parameter values ​​of different soil layers; S104, determining a three-dimensional relationship between soil layer structural parameters and correction coefficients based on a linear relationship; S105, obtaining target soil layer structural parameters and target theoretical values ​​of the target slope; S106, determining a target correction coefficient based on the target soil layer structural parameters and the three-dimensional relationship; S107: Determine the corrected landslide thrust of the target slope based on the target theoretical value and the target correction coefficient.

[0021] Compared with the prior art, the landslide thrust calculation method provided in the present embodiment first obtains the soil layer structure parameters of the slope, and constructs the rock and soil model of the slope, then calculates the soil layer structure parameters to obtain the theoretical value of the slope landslide thrust, performs simulation calculation on the rock and soil model to obtain the simulated value of the slope landslide thrust, takes the ratio of the simulated value to the theoretical value as the correction coefficient, obtains the linear relationship between the structure parameter and the theoretical value, the simulated value and the correction coefficient based on different soil layer structure parameter values, thereby determining the soil layer structure parameters affecting the landslide thrust according to the linear relationship, further determining the three-dimensional relationship between the soil layer structure parameters and the correction coefficient based on the linear relationship, obtaining the target soil layer structure parameters and the target theoretical value of the target slope, determining the target correction coefficient based on the target soil layer structure parameters and the three-dimensional relationship, and determining the corrected landslide thrust of the target slope based on the target theoretical value and the target correction coefficient. The present invention obtains the correction coefficient by calculating a large amount of data, so that when the theoretical value of the target slope is known, the landslide thrust of the target slope in the three-dimensional scene is obtained by the correction coefficient.

[0022] It should be noted that this embodiment is applicable to the calculation of landslide thrust of slope accumulation layer in three-dimensional scenes, such as railways and roads in remote mountainous areas. It is understandable that this embodiment is applicable to the calculation of landslide thrust of slope accumulation layer in any three-dimensional scenes.

[0023] In a specific embodiment of the present invention, in step S101, constructing a rock and soil model of a slope includes: The anti-slide piles on the slope are taken as ideal linear elastic bodies, and the slope rock and soil model of the slope is constructed based on the soil layer parameters.

[0024] The slope rock and soil body is simulated by Mohr-Coulomb constitutive model, and the anti-slide pile is regarded as an ideal linear elastic body. Through field tests and corresponding design specifications, the parameters of the slope rock and soil body and anti-slide pile are selected. Table 1 shows the parameters of rock and soil of different materials.

[0025] Table 1: Parameters of different rock materials

[0026] The hybrid grid generator is used to establish the grid unit. The slope model has a height of 50m, a width of 78m, and a thickness of 18m. The grid division is as follows: Figure 2Among them, the anti-slide piles are simulated by 1D beam elements, and the slope rock and soil materials are simulated by 3D elements. The model contains a total of 30048 elements and 32981 nodes.

[0027] In some embodiments of the present invention, the soil layer structural parameters include: one or more of sliding body volume, sliding body inclination, cohesion and internal friction angle. In some embodiments of the present invention, in step S102, as Figure 3 As shown, the soil layer structure parameters are calculated to obtain the theoretical value of the landslide thrust of the slope, and the rock and soil model is simulated to obtain the simulated value of the landslide thrust of the slope, including: S302, calculating the soil layer structure parameters based on the transfer coefficient method to obtain the theoretical value of the landslide thrust of the slope; S302, simulating and calculating the rock and soil model based on finite element software to obtain a simulated value of the thrust of the slope landslide.

[0028] In some embodiments of the present invention, the method of calculating the soil layer structure parameters based on the transfer coefficient method to obtain the theoretical value of the landslide thrust of the slope includes: When the slope is a railway slope, the rock and soil model is calculated based on the explicit solution of the transfer coefficient method to obtain the theoretical value of the landslide thrust of the slope; When the slope is a highway slope, the rock and soil model is calculated based on the implicit solution of the transfer coefficient method to obtain the theoretical value of the landslide thrust of the slope.

[0029] In a specific embodiment of the present invention, assuming that the landslide slides evenly along the sliding surface, in the design of the railway slope, the transfer coefficient method is usually used to determine the theoretical value of the landslide thrust, and the calculation formula is as follows:

[0030] When the slope is a highway slope, the calculation formula of the theoretical value is:

[0031] In the formula, T i For the i The landslide thrust at the end of each strip; K is the slope safety factor; W i For the i The weight of a strip of soil; α i For the i The inclination angle of the sliding surface where the soil strip is located; φ i For the i The friction angle at the bottom of each soil strip;c i For the i The cohesion at the bottom of each soil strip; L i For the i The length of the sliding surface of each bar; ψ For the i- 1 soil strip transfer coefficient.

[0032] The calculation formula of the theoretical value of landslide thrust in three-dimensional state is equal to the product of the theoretical value of landslide thrust in two-dimensional state and the width of circular pile force. The formula for calculating the width of circular pile force is:

[0033] In the formula, is the calculated width of the anti-slide pile; d is the diameter of the anti-slide pile. The calculated width of the anti-slide pile in this embodiment should be 2.25m.

[0034] Finite element method was used to determine the distribution of landslide thrust on two rows of anti-slide piles, and then compared and verified through engineering measurements. Figure 4 Shown are the finite element analysis values ​​and measured data of landslide thrust distribution.

[0035] After theoretical calculation and finite element simulation using the transfer coefficient method, the finite element simulation value and standard theoretical value (when K=1.00) of the landslide thrust are obtained, as shown in Table 2.

[0036] Table 2 Comparison between simulated and theoretical values ​​of landslide thrust

[0037] According to the results of finite element analysis, the simulated values ​​of the landslide thrust on the upper and lower rows of piles are 643.7 kN and 556.8 kN respectively. The total landslide thrust on the upper row of piles is 15.6% greater than that on the lower row of piles. The simulation values ​​show that the landslide thrust sharing ratios of the upper and lower rows of piles are 53.6% and 46.4% respectively.

[0038] The theoretical value of landslide thrust calculated by the standard method is about 15% larger than the finite element simulation result, and the theoretical result of landslide thrust obtained by the transfer coefficient method is more conservative. The main reasons are as follows: First, the transfer coefficient method assumes that the landslide is a rigid body, ignoring the deformation of the landslide body when it slides down and the anti-slide pile is subjected to force; second, in the theoretical calculation, the calculation width of the anti-slide pile obtained by the calculation formula of the circular pile force width is obtained by an empirical formula, which may have a certain error with the anti-slide pile force condition in the three-dimensional state; in addition, factors such as the slope model meshing and boundary conditions in the finite element simulation may also cause certain deviations.

[0039] The upper row of piles is subject to greater landslide thrust than the lower row of piles. The reasons for this are not only the unique pile-soil distribution position and characteristics of the anti-slide pile slope under study (in this model, the volume of the sliding body behind the two rows of piles is similar), but also the following two aspects: on the one hand, the inclination angle of the landslide body behind the upper row of piles is greater than that of the lower row of piles, resulting in a larger component of the sliding force under the landslide body on the upper row of piles; on the other hand, the landslide thrust generated by the deadweight of the landslide body is partially offset by the upper row of piles, which reduces the landslide thrust on the lower row of piles.

[0040] It is worth noting that when the theoretical formula of the transfer coefficient method is used to calculate the landslide thrust, different values ​​of the safety factor K will produce different results for the theoretical calculation using the explicit solution or implicit solution. Taking the above pile row as an example, when calculating the slope safety factor K from 1.00 to 1.35, the theoretical results of the landslide thrust using the explicit solution and the implicit solution are shown in Table 3.

[0041] Table 3: Comparison of landslide thrust calculated by explicit solution and implicit solution (unit: kN) Safety factor K Explicit solution Implicit solution Difference 1.00 731.8 731.8 0.0 1.05 853.5 812.9 40.6 1.10 975.3 886.6 88.7 1.15 1097.0 953.9 143.1 1.20 1218.7 1015.6 203.1 1.25 1340.5 1072.4 268.1 1.30 1462.2 1124.8 337.4 1.35 1583.9 1173.3 410.6 From the theoretical calculation results of the upper row of piles subjected to landslide thrust, it can be analyzed that when K=1.00, the theoretical calculation results of the two are equal, and the theoretical value is only about 15% larger than the finite element simulation value; when K>1.00, the landslide thrust calculated by the explicit solution is greater than that by the implicit solution, and the larger the K, the greater the difference between the results obtained by the explicit solution and the implicit solution. Specifically, when K=1.05, the landslide thrust calculated by the explicit solution is 40.6kN larger than that by the implicit solution; and when K=1.35, the difference between the two increases to 410.6kN. The theoretical value obtained by the explicit solution is 2.5 times the finite element simulation value, while the result of the implicit solution is 1.8 times the finite element simulation value. The design calculation of slope landslide thrust using the explicit solution is more conservative than the implicit solution.

[0042] The reason why the explicit method is used to calculate the landslide thrust and introduce the safety factor in railway slope design is that the railway bears a large load and the train speed is high. It is necessary to consider more adverse working conditions in the support structure design to ensure the safety and stability of the slope. In the design of highway slope support structure, the vehicle load is relatively small. If the explicit solution is used for the design of anti-slide piles, it will appear too conservative and cause a waste of resources. Therefore, the implicit solution is more suitable for the calculation of landslide thrust of highway slopes.

[0043] In addition, the sliding mass above the upper row of piles has a larger volume and inclination angle; and when the landslide acts on the lower row of piles, part of the deadweight of the soil has been offset by the upper row of piles, resulting in the landslide thrust acting on the upper row of piles being larger than that on the lower row of piles; therefore, in the actual process, special attention should be paid to the force on the upper row of piles.

[0044] In some embodiments of the present invention, the linear relationship between the structural parameter and the theoretical value, the simulation value and the correction coefficient obtained based on the structural parameter values ​​of different soil layers includes: When the value of the slider volume changes, a first linear relationship between the slider volume and the theoretical value, the simulation value, and the correction coefficient is obtained; When the value of the sliding body inclination angle changes, a second linear relationship between the sliding body inclination angle and the theoretical value, the simulation value and the correction coefficient is obtained; When the value of the cohesion changes, a third linear relationship between the cohesion and the theoretical value, the simulation value and the correction coefficient is obtained; When the value of the internal friction angle changes, a fourth linear relationship between the internal friction angle and the theoretical value, the simulation value, and the correction coefficient is obtained.

[0045] In a specific embodiment of the present invention, the first linear relationship is as follows Figure 5 As shown, the second linear relationship is Figure 6 As shown, the third linear relationship is Figure 7 As shown, the third linear relationship is Figure 8 Specifically, the above piles are taken as the main research object, and the corresponding landslide thrust (including the standard theoretical value and the finite element simulation value) when the slope soil layer parameters change is analyzed. The numerical simulation value of the landslide thrust is compared with the theoretical value. The theoretical value is calculated using the transfer coefficient method (safety factor K=1.00), and the quantitative relationship is analyzed. Figure 5 The figure shows the relationship between the sliding body volume and the landslide thrust in the embodiment of the present invention. The main influencing factors selected are the sliding body volume A (first replaced by the horizontal distance L between the top of the slope and the interface of the soil layer, and then converted), the sliding body inclination angle α, the cohesion c, and the internal friction angle φ. In the default setting, L=6m, α=48°, c=18kPa, φ=22°, and the specific steps are as follows: Step 1: When studying the effect of sliding volume on landslide, the horizontal distance L between the top of the slope and the interface of the layer is selected as the independent variable. When L changes, the corresponding sliding volume will also change. Finite element analysis is performed on the slope model under the conditions of L of 5m, 6m, 7m, 8m, and 9m, and the corresponding sliding volume is 64.3m 3 / m 、78.7m 3 / m、94.3m 3 / m、110.9m 3 / m、128.7m 3 / m.

[0046] The influence of landslide thrust on the volume of sliding body is shown in Figure 5 As shown in Figure 2, both the simulated and theoretical values ​​of the landslide thrust gradually increase with the increase of the sliding body volume, and the increase of the simulated value is greater than the theoretical value. The finite element simulation value of the landslide thrust increases from the sliding body volume of 64.3m3 / m 393.8kN gradually increased to 128.7m 3 / m is 1650.5kN, and the theoretical value increases from 508.2kN to 1533.8kN accordingly.

[0047] There is a certain deviation between the landslide thrust calculated by the transfer coefficient method and the finite element simulation result, and the two values ​​are reversed between 94.3m3 / m and 110.9m3 / m. The intersection of the simulation value and the theoretical value curve is located at the sliding volume of 105m 3 / m, indicating that when the sliding volume is less than 105m 3 / m, the simulated value is smaller than the theoretical value; and when the sliding volume is larger than 105m 3 / m, the simulated value is greater than the theoretical value.

[0048] Step 2: If Figure 6 The figure shows the relationship between the sliding body inclination angle and the landslide thrust in the embodiment of the present invention, where the sliding body inclination angles selected as independent variables are 44°, 48°, 52°, 56°, and 60°, respectively; the finite element method is used to simulate the corresponding landslide thrust when the inclination angle changes continuously.

[0049] Figure 6 The influence of the sliding body inclination angle on the thrust of the landslide is shown. With the increase of the sliding body inclination angle, both the simulated value and the theoretical value of the landslide thrust gradually increase. The finite element simulation value of the landslide thrust increases from 481.6 kN when the sliding body inclination angle is 44° to 1205.6 kN when the sliding body inclination angle is 60°. In addition, when the sliding body inclination angle increases from 56° to 60°, the increase in the simulated value of the landslide thrust is greater than ever before. The theoretical value of the landslide thrust is slightly larger than the simulated value, and the theoretical value of the landslide thrust increases linearly from 508.5 kN to 1299.7 kN. When the inclination angle between the landslide body and the sliding bed increases, the sliding force component of the landslide body will increase. At the same time, the friction between the landslide body and the sliding bed decreases, thereby increasing the landslide thrust value.

[0050] Step 3: If Figure 7 The figure shows the relationship between the shear strength index cohesion of the sliding body and the landslide thrust in the embodiment of the present invention. If the cohesion is too small, the sliding force of the landslide body will not be enough to be offset, resulting in the instability of the slope and the increase of the landslide thrust. In order to explore the influence of the landslide thrust controlled by the cohesion, the cohesion is set to 16kPa, 17kPa, 18kPa, 19kPa, and 20kPa respectively.

[0051] In the process of cohesion increasing from 16kPa to 20kPa, the finite element simulation value of the landslide thrust gradually decreased from 789.5kN to 513.2kN, and the theoretical value of the landslide thrust correspondingly decreased steadily from 822.6kN to 640.9kN. Compared with the two influencing factors of sliding body volume and sliding body inclination, the change in landslide thrust caused by the change in cohesion is smaller, indicating that the landslide thrust is less affected by cohesion than by the two factors of sliding body volume and inclination.

[0052] Step 4: If Figure 8 The figure shows the relationship between the internal friction angle, which is the shear strength index of the sliding body, and the landslide thrust in the embodiment of the present invention. The internal friction angle is selected as the independent variable, and other parameters remain unchanged. The landslide thrust results under five working conditions of 20°, 21°, 22°, 23°, and 24° are analyzed.

[0053] The results show that when the internal friction angle of the landslide increases from 20° to 24°, the finite element simulation value of the landslide thrust gradually decreases from 783.8kN to 531.8kN; the theoretical value of the landslide thrust steadily decreases from 819.6kN to 641.3kN. Similar to cohesion, the change in the landslide thrust caused by the internal friction angle is relatively small, and the sensitivity of cohesion and internal friction angle to landslide thrust is basically the same.

[0054] In some embodiments of the present invention, in step S105, Fig. 9 As shown, the three-dimensional relationship between the soil layer structure parameters and the correction coefficients based on the linear relationship is determined, including: S901, determining a three-dimensional relationship between a sliding body volume, a sliding body inclination angle, and a correction coefficient based on a first linear relationship, a second linear relationship, a third linear relationship, and a fourth linear relationship; S902: Determine a three-dimensional relationship among cohesion, internal friction angle and correction coefficient based on the first linear relationship, the second linear relationship, the third linear relationship and the fourth linear relationship.

[0055] In some embodiments of the present invention, the determining of the corrected landslide thrust of the target slope based on the target theoretical value and the target correction coefficient includes: The target theoretical value is multiplied by the target correction coefficient to obtain a corrected landslide thrust of the target slope, wherein the corrected landslide thrust of the target slope is a landslide thrust in a three-dimensional scene.

[0056] In a specific embodiment of the present invention, under three-dimensional conditions, the landslide thrust on the anti-slide pile will also be affected by multiple factors such as the geometric characteristics of the landslide body and the cross-sectional shape of the pile, making the determination of the landslide thrust more complicated. Therefore, the concept of the theoretical formula correction coefficient "n" can be introduced, which is defined as the ratio of the finite element simulation value of the landslide thrust to the theoretical value, and the landslide thrust theoretical formula of the transfer coefficient method is improved. The influence of the combined effect of two factors, namely the volume and inclination of the sliding body, and the cohesion and internal friction angle, on this coefficient is studied separately, so that it can be reasonably applied to the calculation of landslide thrust under three-dimensional conditions. The corrected landslide thrust calculation formula is as follows: in, ——The corrected landslide thrust, that is, the landslide thrust in the three-dimensional state; n ——Correction coefficient of the theoretical formula of landslide thrust. This coefficient may be related to the shape and size of the landslide body and the shear strength index. The specific steps are as follows: Step 1: From S103, we can find that the sliding body volume and sliding body inclination play an important role in the change of landslide thrust; therefore, we analyze the correction coefficient of the theoretical formula n The results of (i.e., simulated value / theoretical value) under different sliding body volumes and sliding body inclinations are shown in Table 4.

[0057] Table 4: Results of correction coefficients of theoretical formula

[0058] like Fig.10 The figure shows a three-dimensional surface diagram of the relationship between the correction coefficient of the theoretical formula in the embodiment of the invention and the volume and inclination of the sliding body. In order to more intuitively reflect the changes in the theoretical value / simulation value, the sliding body volume and the sliding body inclination are used as the horizontal axis, and the theoretical value / simulation value is used as the vertical axis, and a three-dimensional surface is made to show the changes of the theoretical value / simulation value with the sliding body volume and the sliding body inclination.

[0059] Step 2: The shear strength index of the landslide body (cohesion c, internal friction angle φ) also has a profound impact on the magnitude of the landslide thrust. Therefore, the sliding body volume A = 78.7m 3 / m, the sliding body inclination angle β=48° remains unchanged, the independent variable cohesion is selected as 16kPa, 17kPa, 18kPa, 19kPa, and 20kPa, and the internal friction angle is selected as 20°, 21°, 22°, 23°, and 24° respectively. The simulated value and theoretical value of the landslide thrust are analyzed and calculated one by one, and the theoretical formula correction coefficient n under these 25 cases is obtained, as shown in Table 5.

[0060] Table 5: Results of correction coefficients of theoretical formula

[0061] like Fig.11 The figure shows a three-dimensional relationship surface diagram of the correction coefficient of the theoretical formula, cohesion and internal friction angle in an embodiment of the present invention. The three-dimensional relationship surface diagram of the correction coefficient of the theoretical formula, cohesion and internal friction angle is made as above.

[0062] It is worth noting that in order to avoid contingency and verify the reliability and applicability of the above theoretical research, another slope model and measured data of landslide thrust are cited. The length of the anti-slide pile penetrating the landslide section is 14m, and the cross-section of the pile body is a rectangle of 2m×3m.

[0063] The landslide thrust on a single pile obtained by simulation with Midas GTS NX software is 6012.5 kN, while the theoretical result of the landslide thrust on this anti-slide pile obtained by the transfer coefficient method is 6385.4 kN. Therefore, the correction coefficient of the theoretical formula (the ratio of the simulated value to the theoretical value) n=0.942 can be obtained, which is consistent with the correction coefficient range of the landslide thrust theoretical formula obtained above, that is, 0.743~1.274.

[0064] like Fig.12 FIG. 1 is a comparison diagram of the simulated value and the measured value of the landslide thrust distribution in the verification of the embodiment of the present invention. Fig.12 The measured results of thrust distribution of the middle landslide are basically consistent with the finite element simulation results. Most of the measuring points show that the simulated values ​​are slightly larger than the measured values, and the distribution is basically parabolic, which is also consistent with the previous text.

[0065] Table 6 shows the comparison between the simulated and measured values ​​of the landslide thrust. According to the actual measurement at the engineering site, the landslide thrust on the anti-slide pile is 5796.3 kN. The theoretical formula correction coefficient (the ratio of the measured value to the theoretical value) n=0.908 calculated from the engineering measurement is 3.6% smaller than the result obtained from the simulation value. The relative error is within the acceptable range, which verifies the rationality and reliability of the research on the "theoretical formula correction coefficient n".

[0066] Table 6: Comparison between simulated and measured landslide thrust values Analog value Measured value Landslide thrust (kN) 6012.5 5796.3 Theoretical formula correction factor 0.942 0.908 In order to better implement a three-dimensional landslide thrust calculation method in an embodiment of the present invention, based on a three-dimensional landslide thrust calculation method, correspondingly, as follows Fig.13 As shown, an embodiment of the present invention further provides a landslide thrust calculation device, a landslide thrust calculation device 1300 includes: The model building module 1301 is used to obtain the soil layer structure parameters of the slope and build a rock and soil model of the slope; The target data acquisition module 1302 is used to calculate the soil layer structure parameters to obtain the theoretical value of the slope thrust, simulate the rock and soil model to obtain the simulated value of the slope thrust, and use the ratio of the simulated value to the theoretical value as the correction coefficient; A linear relationship determination module 1303 is used to obtain a linear relationship between a structural parameter and a theoretical value, a simulated value, and a correction coefficient based on different soil layer structural parameter values; A three-dimensional relationship determination module 1304 is used to determine the three-dimensional relationship between the soil layer structure parameter and the correction coefficient based on the linear relationship; The target theoretical value acquisition module 1305 is used to obtain the target soil layer structure parameters and target theoretical values ​​of the target slope; A target correction coefficient acquisition module 1306 is used to determine the target correction coefficient based on the target soil layer structure parameters and the three-dimensional relationship; The target landslide thrust acquisition module 1307 is used to determine the corrected landslide thrust of the target slope based on the target theoretical value and the target correction coefficient.

[0067] The landslide thrust calculation device 1300 provided in the above embodiment can implement the technical solution described in the above landslide thrust calculation method embodiment. The specific implementation principles of the above modules or units can refer to the corresponding contents in the above landslide thrust calculation method embodiment, which will not be repeated here.

[0068] like Fig.14 As shown, the present invention also provides an electronic device 1400. The electronic device 1400 includes a processor 1401, a memory 1402 and a display 1403. Fig.14 Only some components of the electronic device 1400 are shown, but it should be understood that it is not required to implement all of the components shown, and more or fewer components may be implemented instead.

[0069] In some embodiments, the processor 1401 may be a central processing unit (CPU), a microprocessor or other data processing chip, and is used to run program codes or process data stored in the memory 1402, such as a landslide thrust calculation method in the present invention.

[0070] In some embodiments, processor 1401 may be a single server or a server group. The server group may be centralized or distributed. In some embodiments, processor 1401 may be local or remote. In some embodiments, processor 1401 may be implemented in a cloud platform. In some embodiments, the cloud platform may include a private cloud, a public cloud, a hybrid cloud, a community cloud, a distributed cloud, an internal cloud, a multi-cloud, etc., or any combination thereof.

[0071] In some embodiments, the memory 1402 may be an internal storage unit of the electronic device 1400, such as a hard disk or memory of the electronic device 1400. In other embodiments, the memory 1402 may also be an external storage device of the electronic device 1400, such as a plug-in hard disk, a smart media card (SMC), a secure digital (SD) card, a flash card, etc., equipped on the electronic device 1400.

[0072] Furthermore, the memory 1402 may include both an internal storage unit of the electronic device 1400 and an external storage device. The memory 1402 is used to store application software installed in the electronic device 1400 and various data.

[0073] In some embodiments, the display 1403 may be an LED display, a liquid crystal display, a touch-sensitive liquid crystal display, an OLED (Organic Light-Emitting Diode) touch device, etc. The display 1403 is used to display information on the electronic device 1400 and to display a visual user interface. The components 1401-1403 of the electronic device 1400 communicate with each other via a system bus.

[0074] In one embodiment, when the processor 1401 executes a landslide thrust calculation program in the memory 1402, the following steps may be implemented: Obtain the soil structure parameters of the slope and construct the rock and soil model of the slope; The soil layer structure parameters are calculated to obtain the theoretical value of the landslide thrust of the slope, and the rock and soil model is simulated to obtain the simulated value of the landslide thrust of the slope, and the ratio of the simulated value to the theoretical value is used as the correction coefficient; Based on the structural parameter values ​​of different soil layers, the linear relationship between the structural parameters and the theoretical values, the simulated values ​​and the correction coefficients is obtained; Determine the three-dimensional relationship between soil structure parameters and correction coefficients based on linear relationships; Obtain target soil layer structural parameters and target theoretical values ​​of the target slope; Determine the target correction coefficient based on the target soil layer structural parameters and three-dimensional relationship; The corrected landslide thrust of the target slope is determined based on the target theoretical value and the target correction coefficient.

[0075] It should be understood that: when the processor 1401 executes a landslide thrust calculation program in the memory 1402, in addition to the above functions, other functions may also be implemented, and details may refer to the description of the corresponding method embodiment above.

[0076] Furthermore, the embodiment of the present invention does not specifically limit the type of the electronic device 1400 mentioned, and the electronic device 1400 may be a portable electronic device such as a mobile phone, a tablet computer, a personal digital assistant (PDA), a wearable device, a laptop computer, etc. Exemplary embodiments of portable electronic devices include but are not limited to portable electronic devices equipped with IOS, Android, Microsoft or other operating systems. The above-mentioned portable electronic devices may also be other portable electronic devices, such as a laptop computer with a touch-sensitive surface (e.g., a touch panel). It should also be understood that in some other embodiments of the present invention, the electronic device 1400 may not be a portable electronic device, but a desktop computer with a touch-sensitive surface (e.g., a touch panel).

[0077] Those skilled in the art will appreciate that all or part of the processes of the above-mentioned embodiments can be implemented by instructing related hardware through a computer program, and the program can be stored in a computer-readable storage medium, wherein the computer-readable storage medium is a disk, an optical disk, a read-only storage memory, or a random access memory, etc.

[0078] The above description is only a preferred specific implementation manner of the present invention, but the protection scope of the present invention is not limited thereto. Any changes or substitutions that can be easily conceived by any technician familiar with the technical field within the technical scope disclosed by the present invention should be covered within the protection scope of the present invention.

Claims

1. A method for calculating landslide thrust, characterized in that: include: Obtain the soil structure parameters of the slope and construct the rock and soil model of the slope; The soil layer structure parameters are calculated to obtain the theoretical value of the landslide thrust of the slope, and the rock and soil model is simulated to obtain the simulated value of the landslide thrust of the slope, and the ratio of the simulated value to the theoretical value is used as the correction coefficient; Based on the structural parameter values ​​of different soil layers, the linear relationship between the structural parameters and the theoretical values, the simulated values ​​and the correction coefficients is obtained; Determine the three-dimensional relationship between soil structure parameters and correction coefficients based on linear relationships; Obtain target soil layer structural parameters and target theoretical values ​​of the target slope; Determine the target correction coefficient based on the target soil layer structural parameters and three-dimensional relationship; The corrected landslide thrust of the target slope is determined based on the target theoretical value and the target correction coefficient.

2. The landslide thrust calculation method according to claim 1, characterized in that: The soil layer structural parameters include: one or more of sliding body volume, sliding body inclination, cohesion and internal friction angle.

3. The landslide thrust calculation method according to claim 1, characterized in that: The method of calculating the soil layer structure parameters to obtain the theoretical value of the landslide thrust of the slope, and performing simulation calculation on the rock and soil model to obtain the simulated value of the landslide thrust of the slope includes: The theoretical value of the landslide thrust is obtained by calculating the soil structure parameters based on the transfer coefficient method. The rock and soil model is simulated and calculated based on finite element software to obtain the simulated value of the slope landslide thrust.

4. The landslide thrust calculation method according to claim 1, characterized in that: The theoretical value of the landslide thrust of the slope is obtained by calculating the soil layer structure parameters based on the transfer coefficient method, including: When the slope is a railway slope, the rock and soil model is calculated based on the explicit solution of the transfer coefficient method to obtain the theoretical value of the landslide thrust of the slope; When the slope is a highway slope, the rock and soil model is calculated based on the implicit solution of the transfer coefficient method to obtain the theoretical value of the landslide thrust of the slope.

5. The landslide thrust calculation method according to claim 2, characterized in that: The linear relationship between the structural parameters and the theoretical values, the simulation values ​​and the correction coefficients obtained based on the structural parameter values ​​of different soil layers includes: When the value of the slider volume changes, a first linear relationship between the slider volume and the theoretical value, the simulation value, and the correction coefficient is obtained; When the value of the sliding body inclination angle changes, a second linear relationship between the sliding body inclination angle and the theoretical value, the simulation value and the correction coefficient is obtained; When the value of the cohesion changes, a third linear relationship between the cohesion and the theoretical value, the simulation value and the correction coefficient is obtained; When the value of the internal friction angle changes, a fourth linear relationship between the internal friction angle and the theoretical value, the simulation value, and the correction coefficient is obtained.

6. The landslide thrust calculation method according to claim 5, characterized in that: The method of determining the three-dimensional relationship between the soil layer structure parameter and the correction coefficient based on the linear relationship includes: Determine a three-dimensional relationship between the sliding body volume, the sliding body inclination and the correction coefficient based on the first linear relationship, the second linear relationship, the third linear relationship and the fourth linear relationship; A three-dimensional relationship among cohesion, internal friction angle and correction coefficient is determined based on the first linear relationship, the second linear relationship, the third linear relationship and the fourth linear relationship.

7. The landslide thrust calculation method according to claim 1, characterized in that: The method of determining the corrected landslide thrust of the target slope based on the target theoretical value and the target correction coefficient includes: The target theoretical value is multiplied by the target correction coefficient to obtain a corrected landslide thrust of the target slope, wherein the corrected landslide thrust of the target slope is a landslide thrust in a three-dimensional scene.

8. A landslide thrust calculation device, characterized in that: include: A model building module is used to obtain the soil structure parameters of the slope and build a rock and soil model of the slope; The target data acquisition module is used to calculate the soil layer structure parameters to obtain the theoretical value of the slope landslide thrust, simulate the rock and soil model to obtain the simulated value of the slope landslide thrust, and use the ratio of the simulated value to the theoretical value as the correction coefficient; A linear relationship determination module is used to obtain the linear relationship between the structural parameters and the theoretical values, simulation values ​​and correction coefficients based on the structural parameter values ​​of different soil layers; A three-dimensional relationship determination module, used for determining the three-dimensional relationship between soil layer structure parameters and correction coefficients based on a linear relationship; A target theoretical value acquisition module is used to obtain target soil layer structural parameters and target theoretical values ​​of a target slope; A target correction coefficient acquisition module is used to determine the target correction coefficient based on the target soil layer structure parameters and the three-dimensional relationship; The target landslide thrust acquisition module is used to determine the corrected landslide thrust of the target slope based on the target theoretical value and the target correction coefficient.

9. An electronic device, characterized in that: comprising a memory and a processor, wherein: The memory is used to store programs; The processor is coupled to the memory and is used to execute the program stored in the memory to implement the steps in the landslide thrust calculation method described in any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that: Used to store computer-readable programs or instructions, which, when executed by a processor, can implement the steps in the landslide thrust calculation method described in any one of claims 1 to 7.

Images