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

Abstract

The application provides a crack position and width calculation method and system based on a point safety coefficient, relates to the technical field of landslide disaster monitoring, and comprises the following steps: obtaining initial data of a target slope sliding body, wherein the initial data comprises geometric size information, mechanical parameters and material characteristics; performing strip block division and model establishment processing according to the initial data, obtaining a strip block model by dividing the sliding body into zones and defining the physical and mechanical characteristics of the strip blocks; performing recursive mechanical calculation processing according to the strip block model, obtaining the stability parameters of the landslide body by calculating the point safety coefficient and the thrust of adjacent strip blocks; and performing crack prediction processing according to the stability parameters, obtaining the position data and width data of the cracks by calculating the horizontal deformation of the strip blocks and the cumulative effect. The application divides the slope sliding body into vertical strip blocks, calculates the safety coefficient of each strip block, and accurately predicts the possible position and width of the cracks in the landslide body by calculating the elastic deformation.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of landslide disaster monitoring, in particular to a crack position and width calculation method and system based on point safety factor. BACKGROUND

[0002] Landslide disaster is an important type of geological disaster, widely distributed in high mountains, hills, plains and coastal areas, especially in areas with complex geological structure and large terrain 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, and have serious social and economic impacts. The existing technology focuses more on the research of landslide stability and formation mechanism, but lacks effective theoretical calculation methods for accurate prediction of crack position and width. As an important precursor of landslide, the position and width of the crack directly affect the stability of the landslide and the disaster prediction. However, the existing technology usually relies on limited field monitoring data, which is difficult to comprehensively and accurately predict the specific position and width of the landslide crack, and often requires a large amount of post-processing human intervention, which has insufficient early warning capability and cannot achieve real-time and effective prevention and control.

[0003] Based on the above-mentioned shortcomings of the prior art, there is an urgent need for a crack position and width calculation method and system based on point safety factor. SUMMARY

[0004] The purpose of the present application is to provide a crack position and width calculation method and system based on point safety factor to improve the above-mentioned problems. In order to achieve the above-mentioned purpose, the technical solution adopted by the present application is as follows:

[0005] In a first aspect, the present application provides a crack position and width calculation method based on point safety factor, comprising:

[0006] Obtaining initial data of the target landslide body, the initial data including geometric size information, mechanical parameters and material properties;

[0007] According to the initial data, performing block division and model establishment processing, by dividing the landslide body into zones and defining the physical and mechanical properties of the blocks, obtaining a block model;

[0008] According to the block model, performing recursive mechanical calculation processing, by calculating the point safety factor and the thrust of adjacent blocks, obtaining the stability parameters of the landslide body;

[0009] According to the stability parameters, performing crack prediction processing, by calculating the horizontal deformation of the blocks and the cumulative effect, obtaining the position data and width data of the cracks.

[0010] In a second aspect, the application also provides a crack position and width calculation system based on point safety factor, comprising:

[0011] An acquisition module is configured to acquire initial data of the target slope sliding body, wherein the initial data comprises geometric size information, mechanical parameters and material characteristics.

[0012] A division module is configured to perform block division and model establishment processing according to the initial data, to obtain a block model by dividing the sliding body into zones and defining physical and mechanical characteristics of the blocks.

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

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

[0015] The application has the following beneficial effects:

[0016] The application divides the slope sliding body into vertical blocks, calculates the safety factor of each block, and accurately predicts the possible position and width of the cracks in the landslide body by calculating the elastic deformation amount. The accurate position and width of the cracks are of great significance for evaluating the deformation degree of the landslide, the position of the sliding surface and the potential harm degree. The method of the application is suitable for stability prediction and early warning of soil landslide and rock landslide, has high practical value in landslide prevention and emergency response, and can provide prospective and reliable technical support for landslide disaster prevention and control. BRIEF DESCRIPTION OF DRAWINGS

[0017] In order to more clearly illustrate the technical solutions of the embodiments of the application, the following will briefly introduce the drawings needed to be used in the embodiments. It should be understood that the following drawings only show some embodiments of the application, and therefore should not be regarded as a limitation to the scope. For those skilled in the art, other related drawings can also be obtained without creative labor.

[0018] Figure 1 A flowchart of a crack position and width calculation method based on point safety factor according to an embodiment of the application;

[0019] Figure 2 A structural diagram of a crack position and width calculation system based on point safety factor according to an embodiment of the application;

[0020] Figure 3 A schematic diagram of block division of a slope sliding body;

[0021] Figure 4 Fig. 1 is a force diagram of a block i;

[0022] Figure 5 Fig. 2 is a geometric feature diagram of a landslide;

[0023] Figure 6 Fig. 3 is a thrust diagram of each block;

[0024] Figure 7 Fig. 4 is a safety factor diagram of each block;

[0025] Figure 8 Fig. 5 is a relative deformation diagram of each block;

[0026] Figure 9 Fig. 6 is a cumulative deformation diagram of each block.

[0027] In the figure, 901 is an acquisition module, 902 is a division module, 903 is a calculation module, and 904 is a prediction module. DETAILED DESCRIPTION

[0028] In order to make the purpose, technical solutions and advantages of the embodiments of the present application clearer, the technical solutions in the embodiments of the present application will be described clearly and completely below with reference to the drawings in the embodiments of the present application. Obviously, the described embodiments are some embodiments of the present application, rather than all the embodiments. The components of the embodiments of the present application described and shown in the drawings can be arranged and designed in various different configurations. Therefore, the following detailed description of the embodiments of the present application provided in the drawings is not intended to limit the scope of the claimed present application, but only represents selected embodiments of the present application. Based on the embodiments in the present application, all other embodiments obtained by those of ordinary skill in the art without creative work belong to the scope of protection of the present application.

[0029] It should be noted that similar reference numerals and letters represent similar items in the following drawings, and therefore, once an item is defined in one drawing, it does not need to be further defined and explained in subsequent drawings. Meanwhile, in the description of the present application, the terms "first", "second", etc. are only used to distinguish the description, and cannot be understood as indicating or implying relative importance.

[0030] Embodiment 1

[0031] The present embodiment provides a crack position and width calculation method based on a point safety factor.

[0032] Referring to Figure 1 , the present method includes steps S100 to S400.

[0033] Step S100, obtaining initial data of the target slope sliding body, the initial data including geometric size information, mechanical parameters and material characteristics;

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

[0035] Step S200, block division and model establishment processing according to the initial data, obtaining the block model by dividing the sliding body and defining the physical and mechanical properties of the blocks;

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

[0037] Step S210, block division processing of the sliding body according to the initial data, obtaining the spatial division structure of the blocks by dividing the sliding body into several vertical blocks from the top of the slope to the foot of the slope, and sequentially numbering the blocks;

[0038] Step S220, stress condition analysis of the blocks according to the spatial division structure, obtaining the basic stress state of the blocks by determining the gravity of each block, the normal pressure and tangential anti-sliding force of the block on the bottom sliding surface, and the cohesion, internal friction angle and sliding surface inclination angle in the sliding surface parameters;

[0039] Step S230, thrust analysis processing between the blocks according to the basic stress state, obtaining the preliminary mechanical model between the blocks by combining the thrust transmitted by the previous block and the mechanical conditions of the block itself.

[0040] Specifically, as shown in Figure 3 and Figure 4 , first, the slope sliding body is divided into several vertical blocks, and these blocks are sequentially numbered from the top of the slope to the foot of the slope as 1 to n. Each block has its own physical and mechanical properties, wherein i represents the serial number of the block, W i represents the gravity of the block i itself, N i represents the normal force of the bottom sliding surface of the block i, T i represents the tangential anti-sliding force of the block i, P i-1P i-1 represents the thrust force transmitted by the last block i-1 i P i represents the thrust force transmitted by the block i to the next block, wherein the direction of the thrust force is parallel to the bottom sliding surface of the block. In addition, the cohesion c and the internal friction angle of the sliding surface are defined α represents the sliding surface inclination angle of the block i i α i-1 represents the sliding surface inclination angle of the block i-1 i-1 b represents the width of the block i These physical and mechanical properties are divided and defined by the block, providing basic data for subsequent mechanical analysis, thrust calculation and crack prediction, so as to obtain the block model describing the mechanical behavior of the slope sliding body.

[0041] Step S300, according to the block model, a recursive mechanical calculation process is carried out, and the safety factor of each block and the thrust force of the adjacent block are calculated to obtain the stability parameters of the landslide body;

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

[0043] Step S310, according to the block model, a single block mechanical balance calculation process is carried out, and the normal pressure and sliding force of the block under the action of gravity are calculated by establishing the force balance equation of the block in the tangential and normal directions of the sliding surface, so as to obtain the basic stress state of the single block;

[0044] It should be noted that the normal pressure and sliding force of the block under the action of gravity are calculated by establishing the force balance equation of the block in the tangential and normal directions of the sliding surface, so as to obtain the basic stress state of the single block. This calculation process helps to judge whether the block is stable, and provides a basis for subsequent point safety factor calculation. The point safety factor is a key to measure the stability of the block. When the point safety factor is less than 1, it indicates that the block is unstable, and the support force needs to be provided by the next block to maintain the balance, so as to recursively calculate the stability state of the entire landslide body.

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

[0046] Step S311, according to the block model, the normal force of the sliding surface is calculated, the normal pressure of the block is calculated by establishing the force balance equation of the block in the normal direction of the sliding surface, combining the gravity of the block, the sliding surface inclination angle and the thrust force transmitted by the last block, and the normal force of the sliding surface is obtained;

[0047] Step S312, according to the normal force, the tangential force of the sliding surface is calculated, the tangential sliding force of the block is calculated by establishing the force balance equation of the block in the tangential direction of the sliding surface, combining the gravity of the block and the sliding surface inclination angle, and the tangential force of the sliding surface is obtained;

[0048] Step S313, according to the normal force and the tangential force of the sliding surface, the anti-sliding force is calculated and processed, the anti-sliding force is substituted into the force balance equation to calculate the point safety factor of the strip block by combining the Mohr-Coulomb yield criterion, and the basic stability state of the strip block is obtained.

[0049] Specifically, the gravity W i Under the action, the force balance equation is established in the tangential direction and the normal direction of the sliding surface, and the following can be obtained:

[0050] The normal pressure N i of the sliding surface is:

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

[0052] The downslide force has:

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

[0054] Wherein, S i represents the downslide force of the strip block i.

[0055] According to the Mohr-Coulomb yield criterion, the anti-sliding force is calculated and substituted into formula (1) and formula (2):

[0056]

[0057] Wherein, F si represents the point safety factor of the strip block i; c i represents the cohesion of the strip block i; l i represents the length of the strip block i; represents the internal friction angle of the strip block i.

[0058] Step S320, according to the basic force state, combining the thrust transmission relationship between the strip blocks, the recursive mechanical calculation between the strip blocks is carried out, the point safety factor of each strip block and the thrust required to be provided by the adjacent strip block are calculated, and whether the strip block is in the limit equilibrium state is judged, and the mechanical property distribution of the strip block is obtained.

[0059] Further, step S320 includes steps S321 to S323.

[0060] Step S321, initial thrust calculation processing is performed according to the basic stress state, whether the top block reaches the limit equilibrium state under its own weight is judged by assuming that the top block is not subjected to the upper thrust, if not balanced, the thrust required by the next block to maintain the stability of the top block is calculated, and the mechanical state parameters of the top block are obtained;

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

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

[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:

[0065]

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

[0067] Step S322, recursive mechanical calculation processing is performed based on the thrust transmission relationship between the blocks, the point safety factor of the current block and the required support force are calculated by block by block, if the current block does not reach the limit equilibrium state, the thrust is transmitted to the next block, until the calculation is performed to the slope foot, and the mechanical state parameters between the blocks are obtained;

[0068] Step S323, stability distribution analysis is performed according to the block mechanical state of the top block and the mechanical state parameters between the blocks, whether each block is in a limit equilibrium state is judged by summarizing the point safety factor of each block and the mechanical relationship with the adjacent block, and the mechanical property distribution of the landslide body is obtained.

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

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

[0071]

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

[0073]

[0074] Further, for the block with index i > 1, there is no upper block thrust, i.e. P i-1 = 0, only formula (6), formula (7) are needed, set P i-1 = 0 to obtain the point safety factor F si of the block i and the next block thrust P i .

[0075] Step S330, overall stability parameter extraction processing is performed according to the mechanical property distribution, and the stability parameters of the landslide body are obtained by integrating the block point safety factor and the inter-block thrust.

[0076] Step S400, crack prediction processing is performed according to the stability parameters, and the position data and width data of the crack are obtained by calculating the horizontal deformation of the block and the cumulative effect.

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

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

[0079] Step S410, block horizontal deformation calculation processing is performed according to the stability parameters, the relative horizontal deformation of the block is calculated by assuming that the block is an elastic body and combining the inter-block thrust, the block width, the elastic modulus and the sliding surface inclination angle, and the horizontal deformation distribution data of each block is obtained;

[0080] Step S420, crack position identification processing is performed according to the horizontal deformation distribution data, the position data of the crack is obtained by positioning the rear edge position of the block with a point safety factor of 1 and combining the deformation state of its adjacent blocks;

[0081] Step S430, crack width accumulation calculation processing is performed according to the position data, and the width data of the crack is obtained by accumulating the horizontal deformation of the block in the crack area.

[0082] Specifically, according to formula (6), (7), the lower block support force P iThis support force can also be understood as the pressure of the i-th block on the next block, or the landslide thrust. At this time, the block is regarded as an elastic body, assuming that the width of the i+1-th block is b i+1 , and the elastic modulus is E i+1 , the relative horizontal deformation amount Δs i+1 of the i+1-th block can be calculated as formula (8). The i+1-th block deforms, and the i-th block is in a limit equilibrium state, so it will follow the i+1-th block to slide and have a horizontal displacement of Δs i+1 . Assuming that the point safety factor of the i-1-th block is greater than 1 and no sliding occurs, the differential horizontal displacement between the i-1-th block and the i-th block forms a horizontal crack width, which is equal to the cumulative value s i of the relative horizontal deformation amount of the block with a number greater than i, as formula (9).

[0083]

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

[0085] Example 2:

[0086] The difference between this embodiment and the above-mentioned embodiment 1 is that this embodiment further calculates the stability coefficient of the landslide.

[0087] Specifically, in order to obtain the stability coefficient of the landslide, the shear strength parameters of the sliding surface need to be reduced, as formula (10). When the anti-slide section fully develops the anti-slide capacity, the point safety factor of the block is 1.0, and the support force required by the last block is 0, the reduction coefficient K is the stability coefficient of the landslide.

[0088]

[0089] where K represents the reduction coefficient; c' i represents the reduced cohesion of the i-th block; and φ represents the reduced internal friction angle of the i-th block.

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

[0091] Example 3:

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

[0093] The geometric characteristics of the landslide are as follows Figure 5 The unit weight of the sliding body 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 body is divided into 25 blocks as shown in Figure 5 ; the thrust P i and the safety factor F si are calculated according to formulas (1)-(7), and the calculation results are shown in Figure 6 、 Figure 7 ; further, the relative deformation and cumulative deformation of each block are calculated according to formulas (8)-(9), and the calculation results are shown in Figure 8 、 Figure 9 , and the crack position is determined to be the end of block 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 factor of the case landslide in Example 3 is calculated by the method of Example 1, Fellenius method, Janbu method, Spencer method and M-P method respectively, and the calculation results are compared, as shown in Table 2.

[0099] Table 2: Comparison of overall safety factor calculation results

[0100]

[0101] It can be seen that the Fellenius method, Janbu method, Spencer method and 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 overall safety factor calculated by the present method has smaller error than the Fellenius method, Janbu method, Spencer method and M-P method.

[0102] Example 5:

[0103] As shown in Figure 2 , the present embodiment provides a crack position and width calculation system based on point safety factor, which comprises:

[0104] an acquisition module configured to acquire initial data of a target slope sliding body, the initial data including geometric size information, mechanical parameters and material characteristics;

[0105] a division module configured to perform block division and model establishment processing according to the initial data, to obtain a block model by dividing the sliding body into zones and defining physical and mechanical characteristics of the blocks;

[0106] a calculation module configured to perform recursive mechanical calculation processing according to the block model, to obtain stability parameters of the landslide by calculating a point safety factor of each block and a pushing force of an adjacent block;

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

[0108] In one specific embodiment of the present application, the division module includes:

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

[0110] a second division unit configured to perform block stress condition analysis according to the spatial division structure, to obtain a basic stress state of each block by determining a gravity of each block, a normal pressure and a tangential anti-sliding force of each block on a bottom sliding surface, and a cohesion, an internal friction angle and a sliding surface inclination angle in the sliding surface parameters;

[0111] a third division unit configured to perform block-to-block pushing force analysis processing according to the basic stress state, to obtain a preliminary mechanical model of the blocks by combining a pushing force transmitted from a previous block and a mechanical condition of each block.

[0112] In one specific embodiment of the present application, the calculation module includes:

[0113] a first calculation unit configured to perform mechanical balance calculation processing on a single block according to the block model, to obtain a basic stress state of the single block by establishing a force balance equation of the block in a tangential direction and a normal direction of the sliding surface, and calculating a normal pressure and a sliding force of the block under the action of gravity;

[0114] a second calculation unit configured to perform recursive mechanical calculation on the blocks according to the basic stress state and a pushing force transmission relationship between the blocks, to obtain a mechanical characteristic distribution of the blocks by calculating a point safety factor of each block and a pushing force required to be provided by an adjacent block, and determining whether the block is in a limit equilibrium state;

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

[0116] In one embodiment of the present application, the first calculation unit comprises:

[0117] The fourth calculation unit is configured to perform sliding surface normal force calculation according to the slice model, and obtain the sliding surface normal force condition by establishing a force balance equation of the slice in the normal direction of the sliding surface, combining the weight of the slice, the sliding surface inclination angle and the thrust transmitted from the previous slice, and calculating the normal pressure of the slice.

[0118] The fifth calculation unit is configured to perform sliding surface tangential force calculation processing according to the normal force condition, and obtain the sliding surface tangential force condition by establishing a force balance equation of the slice in the tangential direction of the sliding surface, combining the weight of the slice and the sliding surface inclination angle, and calculating the tangential sliding force of the slice.

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

[0120] In one embodiment of the present application, the second calculation unit comprises:

[0121] The seventh calculation unit is configured to perform initial thrust calculation processing according to the basic force state, and obtain the mechanical state parameter of the top slice by assuming that the top slice is not subjected to the upper thrust and judging whether it reaches the limit equilibrium state under the action of its own weight, and if not, calculating the thrust required by the next slice to maintain the stability of the top slice.

[0122] The eighth calculation unit is configured to perform recursive mechanical calculation processing based on the thrust transmission relationship between the slices, and obtain the mechanical state parameter between the slices by calculating the point safety factor and the required support force of the current slice, and if the current slice does not reach the limit equilibrium state, transmitting the thrust to the next slice until the calculation reaches the slope toe.

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

[0124] In one embodiment of the present application, the prediction module comprises:

[0125] The first prediction unit is used to calculate the horizontal deformation of the strips based on the stability parameters. By assuming that the strips are elastic bodies and combining the thrust between the strips, the strip width, the elastic modulus and the sliding surface inclination, the relative horizontal deformation of the strips 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, and obtain the crack position data by locating the trailing edge position of the bar with a safety factor of 1 and combining it with the deformation state of its adjacent bars;

[0127] The third prediction unit is used to perform cumulative calculation processing on the crack width according to the position data, and obtain the crack width data by accumulating the horizontal deformation of the bars in the crack area.

[0128] The above description is only a specific embodiment of the present invention, but the scope of protection of the present invention is not limited thereto. Any technician familiar with this technical field can easily think of changes or replacements within the technical scope disclosed by the present invention, which should be covered by the scope of protection 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 landslide, wherein the initial data includes geometric dimension information, mechanical parameters and material properties; Performing strip-block division and model building processing according to the initial data, obtaining a strip-block model by partitioning the sliding body and defining the physical and mechanical properties of the strips; Perform recursive mechanical calculation based on the strip-block model, and obtain the stability parameters of the landslide by calculating the safety factor of each strip-block and the thrust of the adjacent strips. Perform crack prediction processing based on the stability parameters, and obtain crack location data and width data by calculating the horizontal deformation of the strip and the cumulative effect; The crack prediction process is performed based on the stability parameters, and the position and width data of the cracks are obtained by calculating the horizontal deformation of the strips and the cumulative effect, including: Calculate the horizontal deformation of the bars based on the stability parameters, and calculate the relative horizontal deformation of the bars by assuming that the bars are elastic bodies and combining the thrust between the bars, the bar width, the elastic modulus and the sliding surface inclination, so as to obtain the horizontal deformation distribution data of each bar; Perform crack position identification processing based on the horizontal deformation distribution data, and obtain crack position data by locating the trailing edge position of a bar with a safety factor of 1 and combining it with the deformation state of its adjacent bars; The crack width is accumulated and calculated based on the position data, and the crack width data is obtained by accumulating the horizontal deformation of the bars in the crack area.

2. The method for calculating crack position and width based on point safety factor according to claim 1, characterized in that: The strip-block division and model establishment processing are performed based on the initial data. By partitioning the sliding body and defining the physical and mechanical properties of the strips, a strip-block model is obtained, including: Performing strip-block division processing on the sliding body according to the initial data, dividing the sliding body from the top to the foot of the slope into a number of vertical strips, and sequentially numbering the strips to obtain a spatial division structure of the strips; The force conditions of the strips are analyzed according to the spatial division structure. By determining the gravity of each strip, its normal pressure on the bottom sliding surface and tangential anti-slip force, as well as the cohesion, internal friction angle and sliding surface inclination of the sliding surface parameters, the basic force state of the strips is obtained. The thrust analysis between the bars and blocks is performed according to the basic stress state, and a 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 and blocks one by one, including: Perform mechanical balance calculations on a single strip according to the strip model, establish force balance equations for the strip in the tangential and normal directions of the sliding surface, calculate the normal pressure and sliding force of the strip under its own weight, and obtain the basic stress state of the single strip; According to the basic stress state and in combination with the thrust transmission relationship between the bars, a recursive mechanical calculation is performed between the bars. 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 mechanical property distribution of the bars is obtained. 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 based on the strip model. By establishing the force balance equations of the strip in the tangential and normal directions of the sliding surface, the normal pressure and sliding force of the strip under the action of its own weight are calculated, and the basic force state of the single strip is obtained, including: The normal force on the sliding surface is calculated based on the strip model. By establishing a force balance equation of the strip in the normal direction of the sliding surface, the normal pressure of the strip is calculated by combining the deadweight of the strip, the inclination of the sliding surface and the thrust transmitted by the upper strip, and the normal force on the sliding surface is obtained. The tangential force of the sliding surface is calculated based on the normal force condition, and the tangential sliding force of the strip 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 and the tangential force conditions of the sliding surface, the anti-sliding force is calculated. By combining the Mohr-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: Based on the basic stress state and the thrust transmission relationship between the bars, recursive mechanical calculations are performed between the bars. 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 mechanical property distribution of the bars is obtained, including: Performing initial thrust calculation based on the basic stress state, assuming that the top bar is not subjected to upper thrust, and judging whether it has reached a limit equilibrium state under its own weight, 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 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 the calculation reaches the slope foot, and the mechanical state parameters between the strips are obtained. A 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 factors of each bar and its mechanical relationship with 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 the 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 based on 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; A prediction module is used to perform crack prediction processing based on the stability parameters, and obtain crack location data and width data by calculating the horizontal deformation of the strips and the cumulative effect; Wherein, the prediction module includes: a first prediction unit configured to calculate the horizontal deformation of the bars according to the stability parameters, assuming that the bars are elastic bodies and combining the thrust between the bars, the bar width, the elastic modulus and the sliding surface inclination, to calculate the relative horizontal deformation of the bars, thereby obtaining horizontal deformation distribution data of each bar; A second prediction unit is configured to perform crack position identification processing based on the horizontal deformation distribution data, and obtain crack position data by locating the trailing edge position of a bar with a safety factor of 1 and combining the deformation state of its adjacent bars; The third prediction unit is used to perform cumulative calculation processing on the crack width according to the position data, and obtain the crack width data by accumulating the horizontal deformation of the bars in the crack area.

7. The crack position and width calculation system based on point safety factor according to claim 6, characterized in that: The partitioning module includes: A first dividing unit is configured to divide the sliding body into strips and blocks 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 sequentially numbering the strips to obtain a spatial division structure of the strips; The second division unit is used to analyze the stress conditions of the strips and blocks according to the spatial division structure, and obtain the basic stress state of the strips and blocks by determining the gravity of each strip, its normal pressure on the bottom sliding surface and tangential anti-slip 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. The crack position and width calculation system based on point safety factor according to claim 6, characterized in that: The calculation module includes: A first calculation unit is configured to perform mechanical equilibrium calculation processing on a single strip according to the strip model, calculate the normal pressure and sliding force of the strip under its own weight by establishing force balance equations 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 configured to perform recursive mechanical calculations between the bars based on 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 adjacent bars, and determine whether the bars are in a limit equilibrium state, thereby obtaining a distribution of mechanical properties 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 factors of the strip points and the thrust between the strips.

9. The crack position and width calculation system based on point safety factor according to claim 8, characterized in that: The first computing unit includes: a fourth calculation unit, configured to calculate the normal force on the sliding surface based on the strip model, by establishing a force balance equation for the strip in the normal direction of the sliding surface, and combining the deadweight of the strip, the inclination of the sliding surface, and the thrust transmitted by the upper strip to calculate the normal pressure of the strip, thereby obtaining the normal force on the sliding surface; a fifth calculation unit, configured to calculate the tangential force of the sliding surface based on the normal force condition, by establishing a force balance equation for the block in the tangential direction of the sliding surface, combining the deadweight of the block and the inclination of the sliding surface, and calculating the tangential sliding force of the block, thereby obtaining 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 Mohr-Coulomb yield criterion, the anti-sliding force is substituted into the force balance equation to calculate the point safety factor of the strip and block, thereby obtaining the basic stability state of the strip and block.

10. The crack position and width calculation system based on point safety factor according to claim 8, characterized in that: The second calculation unit includes: a seventh calculation unit, configured to perform initial thrust calculation processing based on the basic stress state, assuming that the top bar is not subjected to upper thrust, and determine whether it has reached a limit equilibrium state under the action of its own weight; if it is unbalanced, calculate the thrust required by the next bar to maintain the stability of the top bar, and obtain the mechanical state parameters of the top bar; The eighth calculation unit performs recursive mechanical calculations based on the thrust transfer relationship between the bars. It calculates the point safety factor and required support force of each bar. If the current bar has not reached the limit equilibrium state, the thrust is transferred downward to the next bar until the calculation reaches the slope foot, thereby obtaining the mechanical state parameters between the bars. 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 the adjacent bars, thereby obtaining the mechanical property distribution of the landslide body.

Images