Crack position and width calculation method and system based on point safety coefficient

Abstract

The invention provides a crack position and width calculation method and system based on a point safety coefficient, and relates to the technical field of landslide disaster monitoring, and the method comprises the steps: obtaining initial data of a target slope slip mass, the initial data comprising geometric dimension information, mechanical parameters and material characteristics; performing strip block division and model establishment processing according to the initial data, and obtaining a strip block model by partitioning a sliding mass and defining physical and mechanical characteristics of strip blocks; carrying out recursive mechanical calculation processing according to the strip block model, and calculating point safety coefficients and thrust of adjacent strip blocks one by one to obtain stability parameters of the landslide mass; and performing crack prediction processing according to the stability parameters, and obtaining position data and width data of the crack by calculating the horizontal deformation and the accumulative effect of the strip blocks. According to the method, the slope slip mass is divided into the vertical strip blocks, the safety coefficient of the slope slip mass is calculated one by one, and the possible position and width of the crack in the slip mass are accurately predicted through calculation of the elastic deformation amount.

Description

Technical Field

[0001] The present invention relates to the technical field of landslide disaster monitoring. Specifically, it relates to a method and system for calculating the position and width of cracks based on the point safety factor. Background Art

[0002] Landslide disasters are an important type of geological disasters, widely distributed in mountainous, hilly, plain and coastal areas, especially more frequent in areas with complex geological structures and large topographic undulations. Landslide disasters not only cause a large number of casualties and property losses, but also damage infrastructure, ecological environment and water and soil resources, resulting in serious social and economic impacts. Existing technologies mostly focus on the research of landslide stability and formation mechanisms, but there is still a lack of effective theoretical calculation methods for the accurate prediction of crack positions and widths. Cracks, as important precursors of landslides, their positions and widths directly affect the stability of landslides and disaster prediction. However, existing technologies usually rely on limited on-site monitoring data, it is difficult to comprehensively and accurately predict the specific positions and widths of landslide cracks, and often require a large amount of later manual intervention, with insufficient early warning capabilities and unable to achieve real-time and effective prevention and control.

[0003] Based on the above disadvantages of the existing technologies, there is an urgent need for a method and system for calculating the position and width of cracks based on the point safety factor. Summary of the Invention

[0004] The purpose of the present invention is to provide a method and system for calculating the position and width of cracks based on the point safety factor to improve the above problems. To achieve the above purpose, the technical solutions adopted by the present invention are as follows:

[0005] In the first aspect, the present application provides a method for calculating the position and width of cracks based on the point safety factor, including:

[0006] Obtain the initial data of the sliding mass of the target slope, where the initial data includes geometric dimension information, mechanical parameters and material properties;

[0007] Perform block division and model establishment processing according to the initial data. By partitioning the sliding mass and defining the physical and mechanical properties of each block, a block model is obtained;

[0008] Perform recursive mechanical calculation processing according to the block model. By calculating the point safety factor and the thrust of adjacent blocks for each block, the stability parameters of the landslide body are obtained;

[0009] Perform crack prediction processing according to the stability parameters. By calculating the horizontal deformation amount and cumulative effect of each block, the position data and width data of the cracks are obtained.

[0010] In a second aspect, the present application also provides a system for calculating the position and width of cracks based on the point safety factor, including:

[0011] An acquisition module for acquiring initial data of the target slope landslide, where the initial data includes geometric dimension information, mechanical parameters, and material properties;

[0012] A division module for performing strip division and model establishment processing based on the initial data, and obtaining a strip model by partitioning the landslide body and defining the physical and mechanical properties of the strips;

[0013] A calculation module for performing recursive mechanical calculation processing based on the strip model, and obtaining the stability parameters of the landslide body by calculating the point safety factor and the thrust of adjacent strips for each strip;

[0014] A prediction module for performing crack prediction processing based on the stability parameters, and obtaining the position data and width data of the cracks by calculating the horizontal deformation amount and cumulative effect of the strips.

[0015] The beneficial effects of the present invention are as follows:

[0016] By dividing the slope landslide into vertical strips, calculating the safety factor for each strip, and calculating the elastic deformation amount, the present invention accurately predicts the possible position and width of cracks in the landslide body. The accurate position and width of the cracks are of great significance for evaluating the deformation degree of the landslide, the position of the slip surface, and the potential hazard degree. The method of the present invention is applicable to the stability prediction and early warning of soil landslides and rock landslides, has high practical value in landslide prevention and control and emergency response, and can provide forward-looking and reliable technical support for the prevention and control of landslide disasters. BRIEF DESCRIPTION OF THE DRAWINGS

[0017] In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the following will briefly introduce the drawings required in the embodiments. It should be understood that the following drawings only show some embodiments of the present invention, and therefore should not be regarded as limiting the scope. For those of ordinary skill in the art, other related drawings can be obtained based on these drawings without creative efforts.

[0018] Figure 1 It is a schematic flow chart of a method for calculating the position and width of cracks based on the point safety factor described in the embodiments of the present invention;

[0019] Figure 2 It is a schematic structural diagram of a system for calculating the position and width of cracks based on the point safety factor described in the embodiments of the present invention;

[0020] Figure 3 It is a schematic diagram of strip division of a slope landslide body;

[0021] Figure 4 Schematic diagram of the force on block i;

[0022] Figure 5 Geometric feature map of the landslide;

[0023] Figure 6 Schematic diagram of the thrust of each block;

[0024] Figure 7 Schematic diagram of the safety factor of each block;

[0025] Figure 8 Schematic diagram of the relative deformation of each block;

[0026] Figure 9 Schematic diagram of the cumulative deformation of each block.

[0027] Markings in the figure: 901, acquisition module; 902, division module; 903, calculation module; 904, prediction module. Detailed implementation manners

[0028] To make the objectives, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are some, but not all, of the embodiments of the present invention. Components of the embodiments of the present invention described and illustrated herein generally may be arranged and designed in a variety of different configurations. Therefore, the following detailed description of the embodiments of the present invention provided in the accompanying drawings is not intended to limit the scope of the claimed invention, but is merely representative of selected embodiments of the present invention. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present invention without creative efforts fall within the scope of protection of the present invention.

[0029] It should be noted that like reference numerals and letters denote like items in the following drawings, and thus, once an item is defined in one drawing, it need not be further defined and explained in subsequent drawings. At the same time, in the description of the present invention, the terms "first", "second", etc. are only used for descriptive distinction and cannot be construed as indicating or implying relative importance.

[0030] Embodiment 1:

[0031] This embodiment provides a method for calculating the crack position and width based on the point safety factor.

[0032] Referring to Figure 1 , the figure shows that this method includes steps S100 to S400.

[0033] Step S100: Obtain the initial data of the target slope landslide body. The initial data includes geometric dimension information, mechanical parameters, and material properties.

[0034] It should be noted that the basic assumptions of the calculation method of the present invention include: First, the stability problem of the slope body is regarded as a plane strain problem to simplify the complex three-dimensional mechanical analysis. Second, the sliding block is subjected to a downward sliding force parallel to the sliding surface and a normal pressure perpendicular to the sliding surface on the sliding surface. The failure of the landslide surface follows the Mohr-Coulomb strength criterion, which reflects the shear strength characteristics of the sliding surface material. In addition, the resultant force direction of the forces between blocks is consistent with the inclination angle of the sliding surface of the previous block, ensuring the rationality of the mechanical interaction between the sliding blocks. Although the static equilibrium condition is satisfied along the entire sliding surface, the moment equilibrium condition is not required. For calculation, first, the initial data of the target slope landslide body needs to be obtained, specifically including the geometric dimension information of the landslide, shear strength parameters, and elastic modulus. These data provide a basic basis for subsequent block division, mechanical analysis, and crack prediction.

[0035] Step S200: Perform block division and model establishment processing according to the initial data. By partitioning the landslide body and defining the physical and mechanical properties of the blocks, a block model is obtained.

[0036] Furthermore, step S200 includes steps S210 to S230.

[0037] Step S210: Perform block division processing on the landslide body according to the initial data. By dividing the landslide body into several vertical blocks from the top of the slope to the bottom of the slope and sequentially numbering the blocks, a spatial division structure of the blocks is obtained.

[0038] Step S220: Analyze the force conditions of the blocks according to the spatial division structure. By determining the gravity action of each block, the normal pressure and tangential anti-sliding force of the block on the bottom sliding surface, as well as the cohesion, internal friction angle, and sliding surface inclination angle in the sliding surface parameters, the basic force state of the block is obtained.

[0039] Step S230: Perform thrust analysis processing between the blocks according to the basic force state. By combining the thrust transmitted from the previous block and its own mechanical conditions, a preliminary mechanical model between the blocks is obtained.

[0040] Specifically, as Figure 3 and Figure 4 shown, first, the slope landslide body is divided into several vertical blocks, and these blocks are sequentially numbered from 1 to n from the top of the slope to the bottom of the slope. Each block has its own physical and mechanical properties, where i represents the serial number of the block, W i represents the self-weight of block i, N i represents the normal force of the bottom sliding surface of block i, T i represents the tangential anti-sliding force of block i, P i-1Denote the thrust transmitted by the previous slice i - 1, P i Denote the thrust transmitted by slice i to the next slice, where the direction of the thrust is parallel to the bottom sliding surface of the slice. In addition, the cohesion c and the internal friction angle of the sliding surface are defined The dip angle α of the sliding surface of slice i i and the dip angle α of the sliding surface of slice i - 1 i-1 and the width b of the slice i Through the division and definition of the slices, these physical and mechanical properties provide basic data for subsequent mechanical analysis, thrust calculation, and crack prediction, thus obtaining a slice model that describes the mechanical behavior of the slope sliding mass.

[0041] Step S300: Perform recursive mechanical calculation processing according to the slice model. By calculating the point safety factor and the thrust of adjacent slices slice by slice, obtain the stability parameters of the landslide mass;

[0042] Furthermore, step S300 includes steps S310 to S330.

[0043] Step S310: Perform mechanical equilibrium calculation processing for a single slice according to the slice model. By establishing the force equilibrium equations in the tangential and normal directions of the sliding surface of the slice, calculate the normal pressure and the down - sliding force of the slice under its own weight, and obtain the basic stress state of the single slice;

[0044] It should be noted that in this step, by establishing the force equilibrium equations in the tangential and normal directions of the sliding surface of the slice, the normal pressure and the down - sliding force of the slice under its own weight are calculated, thus obtaining the basic stress state of the single slice. This calculation process helps to judge whether the slice is stable and provides a basis for the subsequent calculation of the point safety factor. The point safety factor is the key to measuring the stability of the slice. When the point safety factor is less than 1, it indicates that the slice is unstable and requires the next slice to provide a supporting force to maintain balance, thereby recursively calculating the stability state of the entire landslide mass.

[0045] Furthermore, step S310 includes steps S311 to S313.

[0046] Step S311: Calculate the normal force on the sliding surface according to the slice model. By establishing the force equilibrium equation in the normal direction of the sliding surface of the slice, and combining the self - weight of the slice, the dip angle of the sliding surface, and the thrust transmitted by the previous slice, calculate the normal pressure of the slice, and obtain the stress condition in the normal direction of the sliding surface;

[0047] Step S312: Perform calculation processing of the tangential force on the sliding surface according to the normal stress condition. By establishing the force equilibrium equation in the tangential direction of the sliding surface of the slice, and combining the self - weight action of the slice and the dip angle of the sliding surface, calculate the tangential down - sliding force of the slice, and obtain the stress condition in the tangential direction of the sliding surface;

[0048] Step S313: Based on the force conditions in the normal and tangential directions of the slip surface, perform anti-sliding force calculation and processing. By combining the Mohr-Coulomb yield criterion, substitute the anti-sliding force into the force equilibrium equation to calculate the point safety factor of the slice, and obtain the basic stability state of the slice.

[0049] Specifically, under the action of gravity W i Establish force equilibrium equations in the tangential and normal directions of the slip surface, and we can get:

[0050] The normal pressure N on the slip surface i There is:

[0051] N i = W i cosα i + P i-1 sin(α i-1 - α i ) (1)

[0052] The downslope force is:

[0053] S i = W i sinα i + P i-1 cos(α i-1 - α i ) - P i (2)

[0054] Among them, S i represents the downslope force received by slice i.

[0055] Calculate the anti-sliding force according to the Mohr-Coulomb yield criterion and substitute it into Equation (1) and Equation (2), we get:

[0056]

[0057] Among them, F si represents the point safety factor of slice i; c i represents the cohesion of slice i; l i represents the length of slice i; represents the internal friction angle of slice i.

[0058] Step S320: Based on the basic force state, combined with the thrust transfer relationship between slices, perform recursive mechanical calculation between slices. By calculating the point safety factor of each slice and the thrust that the adjacent slice needs to provide, and judging whether the slice is in the limit equilibrium state, obtain the mechanical property distribution of the slices;

[0059] Furthermore, Step S320 includes Step S321 to Step S323.

[0060] Step S321: Conduct initial thrust calculation and processing based on the basic stress state. By assuming that the top block is not subjected to the upper thrust and judging whether it reaches the limit equilibrium state under its own weight. If it is not in equilibrium, calculate the thrust required by the next block to maintain the stability of the top block, and obtain the mechanical state parameters of the top block.

[0061] Specifically, for the first block at the top, P i = 0.

[0062] If it is in an unstable state and is in the limit equilibrium state under the support of the lower block, that is, F si = 1, then we can get:

[0063]

[0064] If it is in a stable state, the safety factor under its own weight is the block safety factor, and P i = 0, then there is:

[0065]

[0066] Among them, represents the internal friction angle of block i.

[0067] Step S322: Conduct recursive mechanical calculation and processing based on the thrust transfer relationship between blocks. By calculating the point safety factor and the required support force of the current block for each block. If the current block does not reach the limit equilibrium state, transfer the thrust to the next block until the toe of the slope is calculated, and obtain the mechanical state parameters between blocks.

[0068] Step S323: Conduct stability distribution analysis according to the block mechanical state of the top block and the mechanical state parameters between blocks. By summarizing the point safety factors of each block and their mechanical relationships with adjacent blocks, judge whether each block is in the limit equilibrium state, and obtain the mechanical property distribution of the landslide body.

[0069] For blocks with serial number i > 1 and with the upper block thrust acting, that is, P i-1 ≠ 0.

[0070] If it is in an unstable state and is in the limit equilibrium state under the support of the lower block, that is, F si = 1, the thrust of the lower block can be obtained:

[0071]

[0072] If it is in a stable state, the safety factor under its own weight and the thrust of the previous block is the block safety factor, and P i = 0.

[0073]

[0074] Further, for the blocks with serial number i > 1, there is no thrust from the upper block, that is, P i-1 = 0. In equations (6) and (7), only by setting P i-1 = 0 can the point safety factor F si of block i and the thrust P i .

[0075] Step S330: Extract the overall stability parameters according to the mechanical property distribution. By integrating the point safety factor of the blocks and the thrust between the blocks, the stability parameters of the landslide body are obtained.

[0076] Step S400: Predict the cracks according to the stability parameters. By calculating the horizontal deformation amount and cumulative effect of the blocks, the position data and width data of the cracks are obtained.

[0077] It can be understood that in the calculation process, first, the mechanical state of each block is evaluated through the stability parameters of the blocks (such as thrust and point safety factor), and then the occurrence position of the cracks is predicted. When a block is in the limit equilibrium state, the deformation caused by the thrust on the next block will lead to the generation of cracks. By regarding the blocks as elastic bodies, the horizontal deformation amount of the blocks under the action of thrust is calculated, and further, the deformation amounts of each block are accumulated to determine the position and width of the cracks. These calculations not only reflect the deformation degree of the landslide body, but also help to evaluate the potential landslide hazards, providing a scientific basis for landslide early warning and prevention.

[0078] Further, step S400 includes steps S410 to S430.

[0079] Step S410: Calculate the horizontal deformation of the blocks according to the stability parameters. By assuming that the blocks are elastic bodies and combining the thrust between the blocks, the block width, the elastic modulus, and the dip angle of the slip surface, the relative horizontal deformation amount of the blocks is calculated to obtain the horizontal deformation distribution data of each block.

[0080] Step S420: Identify the crack position according to the horizontal deformation distribution data. By locating the trailing edge position of the block with a point safety factor of 1 and combining the deformation state of its adjacent blocks, the position data of the cracks are obtained.

[0081] Step S430: Accumulate and calculate the crack width according to the position data. By accumulating the horizontal deformation amounts of the blocks in the crack area, the width data of the cracks are obtained.

[0082] Specifically, according to equations (6) and (7), the supporting force P iThis supporting force can also be understood as the pressure of block i on the next block, or the landslide thrust. At this time, the block is regarded as an elastic body, and it is assumed that the width of block i + 1 is b i+1 , and the elastic modulus is E i+1 , and the relative horizontal deformation Δs of block i + 1 can be calculated i+1 , as shown in Equation (8). Block i + 1 deforms, and block i is in a state of limit equilibrium, so it will slide with block i + 1 and have a horizontal displacement of Δs i+1 . It is assumed that the point safety factor of block i - 1 is greater than 1 and it will not slide, then the differential horizontal displacement between block i - 1 and block i forms a horizontal crack width, which is equal to or greater than the cumulative value s of the relative horizontal deformation of block i i , as shown in Equation (9).

[0083]

[0084] Among them, Δs i+1 represents the relative horizontal deformation of block i + 1; n represents the total number of blocks; s i represents the cumulative value of the relative horizontal deformation of block i; E i+1 represents the elastic modulus of block i + 1; h i+1 represents the height of the contact surface between block i + 1 and block i.

[0085] Example 2:

[0086] The difference from the above Example 1 is that in this example, the stability coefficient of the landslide is further calculated.

[0087] Specifically, in order to obtain the stability coefficient of the landslide, the shear strength parameters of the slip surface need to be reduced, as shown in Equation (10). When the anti-sliding section fully exerts its anti-sliding ability, the point safety factor of each block is 1.0, and the required supporting force of the last block is 0, the reduction coefficient K is the stability coefficient of the landslide.

[0088]

[0089] Among them, K represents the reduction coefficient; c' i represents the cohesion of block i after reduction; represents the internal friction angle of block i after reduction.

[0090] This method improves the accuracy of landslide stability assessment, provides more scientific and reliable data support for the risk prediction, emergency response and prevention measures of landslide disasters, and helps to reduce the occurrence frequency and harm of landslide disasters.

[0091] Example 3:

[0092] This embodiment is a specific landslide calculation case.

[0093] The geometric characteristics of the landslide are as follows Figure 5 As shown, the unit weight of the sliding mass is 20 kN / m 3 , the cohesion of the sliding surface is 25 kPa, the internal friction angle is 22°, and the elastic modulus is 60 MPa. First, the sliding mass is divided into 25 slices as Figure 5 shown; the thrust P i and the safety factor F si are calculated according to formulas (1)-(7), and the calculation results are as Figure 6 , Figure 7 ; further, the relative deformation and cumulative deformation of each slice are calculated according to formulas (8)-(9), and the calculation results are as Figure 8 , Figure 9 . Combining the calculation results in Table 1, it can be determined that the crack location is at the end of slice No. 2, and the crack width is 0.071 m.

[0094] Table 1 Calculation results of the case landslide

[0095]

[0096]

[0097] Example 4:

[0098] In this embodiment, the overall safety factors of the case landslide in Example 3 are calculated respectively by using the method in Example 1, the Fellenius method, the Janbu method, the Spencer method, and the M-P method and compared. The calculation results are shown in Table 2 below.

[0099] Table 2 Comparison of calculation results of overall safety factors

[0100]

[0101] It can be seen that the Fellenius method, the Janbu method, the Spencer method, and the M-P method are generally considered to be more accurate methods in the limit equilibrium method. By comparing the calculation results of the above five algorithms, it can be found that the error between the overall safety factor calculated by this method and that calculated by the Fellenius method, the Janbu method, the Spencer method, and the M-P method is relatively small.

[0102] Example 5:

[0103] As Figure 2 shown, this embodiment provides a crack location and width calculation system based on the point safety factor. The system includes:

[0104] An acquisition module, configured to acquire initial data of a target slope landslide, where the initial data includes geometric dimension information, mechanical parameters, and material properties;

[0105] A division module, configured to perform block division and model establishment processing according to the initial data, and obtain a block model by partitioning the landslide body and defining the physical and mechanical properties of each block;

[0106] A calculation module, configured to perform recursive mechanical calculation processing according to the block model, and obtain the stability parameters of the landslide body by calculating the point safety factor of each block and the thrust of adjacent blocks;

[0107] A prediction module, configured to perform crack prediction processing according to the stability parameters, and obtain the position data and width data of the cracks by calculating the horizontal deformation amount and cumulative effect of each block.

[0108] In a specific embodiment of the present invention, the division module includes:

[0109] A first division unit, configured to perform block division processing on the landslide body according to the initial data, and obtain the spatial division structure of the blocks by dividing the landslide body from the top of the slope to the bottom of the slope into several vertical blocks and sequentially numbering the blocks;

[0110] A second division unit, configured to analyze the force-bearing conditions of the blocks according to the spatial division structure, and obtain the basic force-bearing state of the blocks by determining the gravity action, the normal pressure and tangential anti-sliding force of each block on the bottom sliding surface, as well as the cohesion, internal friction angle, and sliding surface inclination angle in the sliding surface parameters;

[0111] A third division unit, configured to perform thrust analysis processing between the blocks according to the basic force-bearing state, and obtain a preliminary mechanical model between the blocks by combining the thrust transmitted from the previous block and its own mechanical conditions.

[0112] In a specific embodiment of the present invention, the calculation module includes:

[0113] A first calculation unit, configured to perform mechanical equilibrium calculation processing for a single block according to the block model, and obtain the basic force-bearing state of the single block by establishing force equilibrium equations in the tangential and normal directions of the sliding surface and calculating the normal pressure and sliding force of the block under its own weight;

[0114] A second calculation unit, configured to perform recursive mechanical calculation between the blocks according to the basic force-bearing state and in combination with the thrust transfer relationship between the blocks, and obtain the mechanical property distribution of the blocks by calculating the point safety factor of each block and the thrust required by adjacent blocks, and determining whether the block is in a limit equilibrium state;

[0115] A third calculation unit, configured to perform overall stability parameter extraction processing according to the mechanical property distribution, and obtain the stability parameters of the landslide body by integrating the slice point safety factor and the thrust between slices.

[0116] In a specific embodiment of the present invention, the first calculation unit includes:

[0117] A fourth calculation unit, configured to calculate the normal force of the sliding surface according to the slice model, establish a force balance equation of the slice in the normal direction of the sliding surface, and calculate the normal positive pressure of the slice by combining the self-weight of the slice, the inclination angle of the sliding surface, and the thrust transmitted from the previous slice, so as to obtain the normal force condition of the sliding surface;

[0118] A fifth calculation unit, configured to perform sliding surface tangential force calculation processing according to the normal force condition, establish a force balance equation of the slice in the tangential direction of the sliding surface, and calculate the tangential sliding force of the slice by combining the self-weight action of the slice and the inclination angle of the sliding surface, so as to obtain the tangential force condition of the sliding surface;

[0119] A sixth calculation unit, configured to perform anti-sliding force calculation processing according to the normal force condition and the tangential force condition of the sliding surface, substitute the anti-sliding force into the force balance equation by combining the Mohr-Coulomb yield criterion to calculate the point safety factor of the slice, and obtain the basic stability state of the slice.

[0120] In a specific embodiment of the present invention, the second calculation unit includes:

[0121] A seventh calculation unit, configured to perform initial thrust calculation processing according to the basic force state, assume that the top slice is not subjected to the upper thrust, judge whether it reaches the limit equilibrium state under its self-weight, and if it is not balanced, calculate the thrust required by the next slice to maintain the stability of the top slice, so as to obtain the mechanical state parameters of the top slice;

[0122] An eighth calculation unit, configured to perform recursive mechanical calculation processing based on the thrust transfer relationship between slices, calculate the point safety factor and the required support force of the current slice slice by slice, and if the current slice does not reach the limit equilibrium state, transfer the thrust to the next slice until the toe of the slope is calculated, so as to obtain the mechanical state parameters between slices;

[0123] A ninth calculation unit, configured to perform stability distribution analysis according to the slice mechanical state of the top slice and the mechanical state parameters between slices, judge whether each slice is in the limit equilibrium state by summarizing the point safety factor of each slice and its mechanical relationship with adjacent slices, and obtain the mechanical property distribution of the landslide body.

[0124] In a specific embodiment of the present invention, the prediction module includes:

[0125] The first prediction unit is used to perform strip horizontal deformation calculation and processing based on stability parameters. By assuming that the strip is an elastic body and combining the thrust between strips, strip width, elastic modulus, and slip surface dip angle, the relative horizontal deformation of the strip is calculated to obtain the horizontal deformation distribution data of each strip.

[0126] The second prediction unit is used to perform crack position identification processing based on the horizontal deformation distribution data. By locating the trailing edge position of the strip with a safety factor of 1 at the positioning point and combining the deformation states of its adjacent strips, the position data of the crack is obtained.

[0127] The third prediction unit is used to perform crack width accumulation calculation processing based on the position data. By accumulating the horizontal deformation of the strips in the crack area, the width data of the crack is obtained.

[0128] As mentioned above, the above is only the specific implementation manner of the present invention, but the protection scope of the present invention is not limited thereto. Any person skilled in the art within the technical scope disclosed by the present invention can easily think of changes or substitutions, which should be covered within the protection scope of the present invention.

Claims

1. A method for calculating crack position and width based on point safety factor, characterized in that: include: Acquiring initial data of the target slope sliding body, wherein the initial data includes geometric dimension information, mechanical parameters and material properties; According to the initial data, strip-block division and model establishment processing are performed, and a strip-block model is obtained by dividing the sliding body into zones and defining the physical and mechanical properties of the strips; According to the strip-block model, recursive mechanical calculation processing is performed to obtain the stability parameters of the landslide body by calculating the point safety factor and the thrust of the adjacent strips and blocks one by one; The crack prediction process is performed according to the stability parameters, and the position data and width data of the crack are obtained by calculating the horizontal deformation amount and cumulative effect of the strip.

2. The method for calculating crack position and width based on point safety factor according to claim 1, characterized in that: According to the initial data, strip-block division and model establishment processing are performed, and the strip-block model is obtained by dividing the sliding body into zones and defining the physical and mechanical properties of the strips, including: According to the initial data, the sliding body is divided into strips and blocks, by dividing the sliding body from the top to the foot of the slope into a number of vertical strips and numbering the strips in sequence, thereby obtaining a spatial division structure of the strips and blocks; According to the spatial division structure, the force condition of the strips is analyzed, and the basic force state of the strips is obtained by determining the gravity of each strip, its normal pressure on the bottom sliding surface and the tangential anti-sliding force, as well as the cohesion, internal friction angle and sliding surface inclination angle of the sliding surface parameters; According to the basic stress state, the thrust analysis and processing between the bars and blocks is carried out, and the preliminary mechanical model between the bars and blocks is obtained by combining the thrust transmitted by the previous bar and its own mechanical conditions.

3. The method for calculating crack position and width based on point safety factor according to claim 1, characterized in that: According to the strip-block model, recursive mechanical calculation processing is performed, and the stability parameters of the landslide body are obtained by calculating the point safety factor and the thrust of the adjacent strips one by one, including: According to the strip model, the mechanical balance calculation of a single strip is performed, and the normal pressure and sliding force of the strip under the action of its own weight are calculated by establishing the force balance equation of the strip in the tangential and normal directions of the sliding surface, so as to obtain the basic stress state of the single strip; According to the basic stress state, combined with the thrust transmission relationship between the bars, recursive mechanical calculation between the bars is performed, and the mechanical property distribution of the bars is obtained by calculating the point safety factor of each bar and the thrust required to be provided by the adjacent bars, and judging whether the bars are in a limit equilibrium state; The overall stability parameter extraction process is performed according to the mechanical characteristic distribution, and the stability parameter of the landslide body is obtained by integrating the safety factor of the strip-block point and the thrust between the strips and blocks.

4. The method for calculating crack position and width based on point safety factor according to claim 3, characterized in that: The mechanical balance calculation of a single strip is performed according to the strip model. By establishing the force balance equation of the strip in the tangential and normal directions of the sliding surface, the normal pressure and the sliding force of the strip under the action of its own weight are calculated, and the basic stress state of the single strip is obtained, including: The normal force of the sliding surface is calculated according to the strip model, and the normal positive pressure of the strip is calculated by establishing the force balance equation of the strip in the normal direction of the sliding surface, combining the deadweight of the strip, the inclination of the sliding surface and the thrust transmitted by the previous strip, so as to obtain the normal force of the sliding surface; The tangential force of the sliding surface is calculated according to the normal force condition, and the tangential sliding force of the sliding surface is calculated by establishing a force balance equation of the strip in the tangential direction of the sliding surface, combining the deadweight of the strip and the inclination of the sliding surface, and obtaining the tangential force condition of the sliding surface; According to the normal force conditions of the sliding surface and the tangential force conditions of the sliding surface, the anti-sliding force is calculated and processed. By combining the Moore-Coulomb yield criterion, the anti-sliding force is substituted into the force balance equation to calculate the point safety factor of the strip and the basic stability state of the strip is obtained.

5. The method for calculating crack position and width based on point safety factor according to claim 3, characterized in that: According to the basic stress state, combined with the thrust transmission relationship between the bars, recursive mechanical calculation between the bars is performed, and the mechanical property distribution of the bars is obtained by calculating the point safety factor of each bar and the thrust required to be provided by the adjacent bars, and judging whether the bars are in the limit equilibrium state, including: Performing initial thrust calculation processing according to the basic stress state, assuming that the top bar is not subjected to upper thrust, judging whether it has reached a limit equilibrium state under its own weight, and if it is unbalanced, calculating the thrust required by the next bar to maintain the stability of the top bar, and obtaining the mechanical state parameters of the top bar; Based on the thrust transfer relationship between strips, recursive mechanical calculation is performed. The point safety factor and the required support force of the current strip are calculated strip by strip. If the current strip has not reached the limit equilibrium state, the thrust is transferred downward to the next strip until it is calculated to the slope foot, and the mechanical state parameters between the strips are obtained. The stability distribution analysis is performed based on the mechanical state of the top bar and the mechanical state parameters between the bars. By summarizing the point safety factor of each bar and its mechanical relationship with the adjacent bars, it is determined whether each bar is in a limit equilibrium state, and the mechanical property distribution of the landslide body is obtained.

6. A crack position and width calculation system based on point safety factor, characterized in that: include: An acquisition module is used to acquire initial data of a target slope sliding body, wherein the initial data includes geometric dimension information, mechanical parameters and material properties; A partitioning module is used to perform strip-block partitioning and model building processing according to the initial data, and obtain a strip-block model by partitioning the sliding body and defining the physical and mechanical properties of the strips; A calculation module is used to perform recursive mechanical calculation processing according to the strip-block model, and obtain the stability parameters of the landslide body by calculating the point safety factor and the thrust of the adjacent strips and blocks one by one; The prediction module is used to perform crack prediction processing according to the stability parameters, and obtain the position data and width data of the cracks by calculating the horizontal deformation amount and cumulative effect of the bars.

7. A crack position and width calculation system based on point safety factor according to claim 6, characterized in that: The division module comprises: A first division unit is used to perform strip-block division processing on the sliding body according to the initial data, by dividing the sliding body from the top to the foot of the slope into a plurality of vertical strips, and numbering the strips in sequence to obtain a spatial division structure of the strips; The second division unit is used to analyze the force conditions of the strips and blocks according to the spatial division structure, and obtain the basic force state of the strips and blocks by determining the gravity of each strip and block, its normal pressure on the bottom sliding surface and tangential anti-sliding force, as well as the cohesion, internal friction angle and sliding surface inclination angle among the sliding surface parameters; The third division unit is used to analyze and process the thrust between the bars and blocks according to the basic stress state, and obtain a preliminary mechanical model between the bars and blocks by combining the thrust transmitted by the previous bar and its own mechanical conditions.

8. A crack position and width calculation system based on point safety factor according to claim 6, characterized in that: The calculation module comprises: The first calculation unit is used to perform mechanical balance calculation processing of a single strip according to the strip model, calculate the normal pressure and sliding force of the strip under the action of its own weight by establishing the force balance equation of the strip in the tangential and normal directions of the sliding surface, and obtain the basic force state of the single strip; A second calculation unit is used to perform recursive mechanical calculation between the bars according to the basic stress state and the thrust transmission relationship between the bars, calculate the point safety factor of each bar and the thrust required to be provided by the adjacent bars, and determine whether the bars are in a limit equilibrium state, so as to obtain the distribution of mechanical characteristics of the bars; The third calculation unit is used to extract and process the overall stability parameters according to the mechanical characteristic distribution, and obtain the stability parameters of the landslide body by integrating the safety factor of the strip point and the thrust between the strips.

9. A crack position and width calculation system based on point safety factor according to claim 8, characterized in that: The first computing unit comprises: The fourth calculation unit is used to calculate the normal force of the sliding surface according to the strip model, by establishing a force balance equation of the strip in the normal direction of the sliding surface, combining the deadweight of the strip, the inclination of the sliding surface and the thrust transmitted by the previous strip, to calculate the normal positive pressure of the strip, and obtain the normal force of the sliding surface; A fifth calculation unit is used to calculate the tangential force of the sliding surface according to the normal force condition, by establishing a force balance equation of the strip in the tangential direction of the sliding surface, combining the deadweight of the strip and the inclination of the sliding surface, and calculating the tangential sliding force of the strip to obtain the tangential force condition of the sliding surface; The sixth calculation unit is used to calculate the anti-sliding force according to the normal force conditions of the sliding surface and the tangential force conditions of the sliding surface. By combining the Moore-Coulomb yield criterion, the anti-sliding force is substituted into the force balance equation to calculate the point safety factor of the strip and the basic stability state of the strip is obtained.

10. A crack position and width calculation system based on point safety factor according to claim 8, characterized in that: The second computing unit comprises: a seventh calculation unit, configured to perform initial thrust calculation processing according to the basic stress state, by assuming that the top bar is not subjected to upper thrust, and judging whether it reaches a limit equilibrium state under its own weight; if it is unbalanced, calculating the thrust that the next bar needs to provide to maintain the stability of the top bar, and obtaining the mechanical state parameters of the top bar; The eighth calculation unit performs recursive mechanical calculation based on the thrust transfer relationship between the strips and blocks, calculates the point safety factor and the required support force of the current strip block by strip block, and if the current strip block does not reach the limit equilibrium state, transfers the thrust downward to the next strip block until the calculation reaches the slope foot, thereby obtaining the mechanical state parameters between the strip blocks; The ninth calculation unit is used to perform stability distribution analysis based on the mechanical state of the top bar and the mechanical state parameters between the bars, and to determine whether each bar is in a limit equilibrium state by summarizing the point safety factor of each bar and its mechanical relationship with adjacent bars, so as to obtain the mechanical property distribution of the landslide body.

Images