CN113627040A

CN113627040A - Heterogeneous slope stability analysis method

Info

Publication number: CN113627040A
Application number: CN202111185040.4A
Authority: CN
Inventors: 邹金锋; 熊湘瑜; 刘坤; 陈嘉祺; 张兰纯
Original assignee: Central South University
Current assignee: Central South University
Priority date: 2021-10-12
Filing date: 2021-10-12
Publication date: 2021-11-09
Anticipated expiration: 2041-10-12
Also published as: CN113627040B

Abstract

The invention provides a stability analysis method for a heterogeneous slope, which comprises the following steps: 1. acquiring a two-dimensional projection of the block stone in the engineering area, 2, constructing a 24-sided polygon of the two-dimensional projection of the block stone, and establishing a database for recording vertex coordinates of the 24-sided polygon; 3. generating a block stone three-dimensional graph of the block stone two-dimensional projection; 4. generating a three-dimensional slope model of the engineering area; 5. establishing a heterogeneous slope model; 6. endowing geotechnical material data, boundary limiting conditions and acting force parameters of the heterogeneous slope model engineering area; 7. and rapidly generating a slope stability analysis report according to a finite element strength reduction method. The method effectively improves the accuracy of the stability analysis of the heterogeneous slope, meanwhile, the content of the slope stability analysis report is richer, reliable guarantee is provided for engineering numerical simulation, and a reasonable support measure suggestion basis is provided for the slope which may generate instability.

Description

Heterogeneous slope stability analysis method

Technical Field

The invention relates to a heterogeneous slope stability analysis method, and belongs to the field of engineering geological survey analysis.

Background

The soil-rock mixture widely exists in nature and engineering construction, and is a heterogeneous rock-soil material composed of soil and rock blocks in a certain proportion, under normal conditions, the material characteristics of the soil and the rock blocks forming the soil-rock mixture have obvious difference, a soil-rock mixed filling heterogeneous side slope is similar to other composite structure materials as a composite structure composed of the rock blocks and soil bodies, the mechanical property of the soil-rock mixed filling heterogeneous side slope depends on the property and the spatial distribution of each component, namely the soil-rock mixed filling heterogeneous side slope is influenced by the properties of the soil bodies and the rock blocks in the side slope, and meanwhile the soil-rock mixed filling heterogeneous side slope is also influenced by the spatial distribution of the rock blocks in the side slope. In the existing numerical analysis and calculation, the heterogeneous slope filled with soil and rock is generally regarded as a homogeneous soil slope or rock slope, and the parameters are approximately replaced by the parameters of fine grain components. The finally calculated result is greatly different from the actual result, and the side slope is composed of the rock blocks and the soil body from a microscopic view, so that the influence of the microscopic structure on the macroscopic performance of the side slope is not negligible. Meanwhile, the current particle flow simulation of the soil-rock mixture is mainly focused on the two-dimensional situation, the influence of the shape of the block stone is mostly not considered, the influence of the content of the block stone and the random distribution of the block stone on the stability of the soil-rock mixed heterogeneous slope is not systematically considered, and the microscopic mechanism for deformation and damage of the soil-rock mixed heterogeneous slope is not deeply researched.

The Xuwenjie [1] establishes a microscopic structure model of an earth-stone mixture based on a digital image processing technology, but because the stability analysis of a heterogeneous slope of a two-dimensional layer is only considered in the modeling process, the model is not consistent with the actual situation, so that a calculation result has larger error, and the application range is greatly limited.

When the gold epitaxy [2] judges whether the side slope can generate instability damage by using particle flow intensity reduction and providing an energy evolution-based method, a three-dimensional discrete element model of a soil-stone mixture side slope mesoscopic structure is established by combining the developed three-dimensional discrete element modeling method of irregular block stones and soil-stone mixtures, but the application range of the gold epitaxy has great limitation because the blocks are simplified into regular spheres in the modeling process.

Liu shun Qing [3] provides an analysis method for stability of a soil-stone mixed slope considering random distribution and content of stones, but circular base materials with different diameters are randomly put in to generate an initial soil-stone mixed slope model, a polygonal file with edges and corners is obtained through setting of the maximum number of sides of a polygon, but the polygonal file is not combined with an actual stone model when the stone model file is established, so that the established stone polygon can not well replace the influence of actual stones on the slope, stays on a two-dimensional slope calculation layer when the model stability is analyzed, and does not conform to the actual situation.

Zhao Xin bright [4] based on the existing block stone shape database random generation block stone soil slope mesoscopic model, although its block stone two-dimensional projection of establishing accords with real block stone shape, its model of analysis still is two-dimensional soil stone mixture model, can not simulate the slope in the real condition well, does not accord with actual conditions, therefore can produce great difference.

In conclusion, the existing three-dimensional soil-rock mixed filling heterogeneous slope still has a lot of defects in numerical analysis and calculation, a modeling method is not mature enough, most of the existing numerical simulation still stays in a two-dimensional layer, the generation technology of random block rocks is not common enough, the influence of the block shape on the stability is often ignored, and the simplification method is not in line with the actual situation, so that the calculation result has certain deviation from the actual situation. Therefore, an effective method for truly reducing the non-homogeneous side slope in the engineering is urgently needed to be researched, so that the stability and the failure mode of the soil-rock mixed heterogeneous side slope are more accurately analyzed, and the construction requirement of the side slope engineering is better met.

Reference documents:

[1] xuwenjie, Wang Yujie, Chenzuyu, Hurelin the slope stability analysis of soil-rock mixture based on digital image technology [ J ]. geotechnical force, 2008,29(S1): 341-346-.

[2] Three-dimensional particle discrete element analysis of the stability of the side slope of the mixture of gold, great Yawu, Stroke, Lijing, earth and stone [ J ]. proceedings of Harbin university of industry 2020,52(02):41-50.

[3] Liu Shunqing, Huang Dong, Zhou Ai, Chua Jun, Jiang Pengming, study on soil-rock mixed slope stability analysis method based on random block-rock model [ J ] geotechnical, 2019,40(S1): 350-reservoir 358.

[4] The numerical simulation research on the influence of space and particle parameters on the stability of the rocky soil slope [ J ] engineering science and technology 2020,52(04):166 + 175.

Disclosure of Invention

Aiming at the problems to be solved, the invention provides a heterogeneous slope stability analysis method which generates a three-dimensional graph of the rock block in any shape by utilizing a two-dimensional projection image of the real rock block and randomly puts the three-dimensional graph of the rock block into a three-dimensional slope.

In order to solve the technical problem, the invention provides a method for analyzing the stability of a heterogeneous slope, which comprises the following steps:

step 1: collecting two-dimensional projection of the block stone in the engineering area;

step 2: constructing 24-sided polygons of the two-dimensional projection of the block stones through MATLAB, and establishing a database for recording vertex coordinates of the 24-sided polygons of the two-dimensional projection of the block stones;

and step 3: generating a block stone three-dimensional graph of the block stone two-dimensional projection through MATLAB;

and 4, step 4: acquiring a slope image of the engineering area, and generating a three-dimensional slope model of the engineering area through a three-dimensional reconstruction technology;

and 5: embedding the three-dimensional graph of the rock block into a three-dimensional slope model of the engineering area, and establishing a heterogeneous slope model;

step 6: endowing geotechnical material data, boundary limiting conditions and acting force parameters of the heterogeneous slope model engineering area;

and 7: and analyzing the heterogeneous slope model according to the finite element strength reduction method, and quickly generating a slope stability analysis report.

Further, the step 2 specifically includes:

2.1: converting the two-dimensional projection of the block stone into a binary image through MATLAB, and acquiring an external contour line of the two-dimensional projection of the block stone;

2.2: constructing 24-sided polygon of block stone two-dimensional projection by using external contour line of block stone two-dimensional projection, and establishing 24-sided polygon for recording block stone two-dimensional projectionOf each vertex, the coordinates of each vertex including the abscissa X_iAnd ordinate Y_iWherein i is the vertex serial number of the 24-sided polygon of the two-dimensional projection of the block stone, i is an integer, i is within a value range of 1 and not more than 24, the serial numbers of the vertices are sequentially recorded in an anticlockwise direction by taking the positive direction of an X axis as the starting direction, and the method for determining the 24-sided polygon of the two-dimensional projection of the block stone and the vertex coordinates thereof is as follows:

(1) when the number of nodes on the external contour line of the two-dimensional block stone projection is 24, 24 polygons of the two-dimensional block stone projection correspond to the external contour line of the two-dimensional block stone projection, 24 vertexes of the 24 polygons of the two-dimensional block stone projection are 24 nodes on the external contour line of the two-dimensional block stone projection respectively, and the coordinates of the 24 vertexes of the 24 polygons of the two-dimensional block stone projection correspond to the coordinates of the 24 nodes on the external contour line of the two-dimensional block stone projection one by one;

(2) when the number of nodes on the external contour line of the two-dimensional projection of the block stone is more than 24, calculating the length of each edge on the external contour line of the two-dimensional projection of the block stone, finding the edge with the shortest length and the shorter edge of the two adjacent edges of the edge, deleting the node shared by the shorter edge and the edge with the shortest length, and connecting the other node left by the shorter edge and the other node left by the edge with the shortest length to form the external contour line of the two-dimensional projection of the new block stone;

(3) when the number of nodes on the external contour line of the two-dimensional projection of the block stone is less than 24, calculating the length of each edge on the external contour line of the two-dimensional projection of the block stone, finding the edge with the longest length, respectively connecting the midpoint of the edge with the longest length with two nodes of the edge with the longest length, and not connecting the two nodes of the edge with the longest length to form a new external contour line of the two-dimensional projection of the block stone;

2.3: constructing 24-sided polygons of the two-dimensional projections of the block stones acquired in the step 1 by using a method of 2.1-2.2, and establishing a database for recording each vertex coordinate of the 24-sided polygons of the two-dimensional projections of each block stone;

further, the step 3 specifically includes:

3.1: establishing a three-dimensional space coordinate system, extracting a 24-sided polygon of the two-dimensional projection of the block stone, moving the centroid of the 24-sided polygon of the two-dimensional projection of the block stone to the origin of the three-dimensional space coordinate system, and assigning the Z coordinates of 24 vertexes of the 24-sided polygon of the two-dimensional projection of the block stone to be 0 to obtain a three-dimensional center plane of the block stone, wherein the Z coordinates of all vertexes of the 24-sided polygon of the two-dimensional projection of the block stone in the three-dimensional space coordinate system are 0;

3.2: based on the three-dimensional center plane of the block stone, four sub-planes parallel to the three-dimensional center plane of the block stone are generated, namely a sub-plane 1, a sub-plane 2, a sub-plane 3 and a sub-plane 4, wherein the sub-plane 2 is adjacent to and above the three-dimensional center plane of the block stone, the sub-plane 3 is adjacent to and below the three-dimensional center plane of the block stone, the sub-plane 1 is parallel to, adjacent to and above the sub-plane 2, and the sub-plane 4 is parallel to, adjacent to and below the sub-plane 3, and the specific method is as follows:

(1) defining the ratio E of the height of the stone in the Z direction to the width of the stone in the X direction, wherein 0< E < 1;

(2) generating a sub-plane 2, generating coordinate values of 12 vertexes of a 12-sided polygon of the sub-plane 2 based on the coordinate values of 12 odd vertexes of a 24-sided polygon of the three-dimensional center plane of the block stone, and numbering the 12 vertexes of the 12-sided polygon of the sub-plane 2, wherein the vertexes numbered 1 to 12 of the 12-sided polygon of the sub-plane 2 correspond to the vertexes numbered 1, 3, 5, 7, 9, 11, 13, 15, 17, 19, 21, 23 of the 24-sided polygon of the three-dimensional center plane of the block stone one-to-one, the Y coordinates of the vertexes numbered 1 to 12 of the 12-sided polygon of the sub-plane 2 are the Y coordinates of the vertexes numbered 1, 3, 5, 7, 9, 11, 13, 15, 17, 19, 21, 23 of the three-sided polygon of the three-dimensional center plane of the block stone multiplied by 0.8, the X coordinates of the vertexes numbered 1 to 12 of the 12-sided polygon of the sub-sided polygon 2 are the X coordinates of the vertexes numbered 1, 12-sided polygon of the three-sided center plane of the block stone, 3. 5, 7, 9, 11, 13, 15, 17, 19, 21, 23, defining a distribution range H of Z coordinates of 12 vertices of a 12-sided polygon of the sub-plane 2, 0.7H 0.825, generating 4 values between 0.7E and 0.825E according to a unifrnd function, ordering the 4 values from small to large and respectively represented by A, B, C and D, determining a value of the Z coordinate of a vertex with a number 1 of the 12-sided polygon of the sub-plane 2 as a minimum value a of the 4 values, determining a value of the Z coordinate of a vertex with a number 4 of the 12-sided polygon of the sub-plane 2 as a second small value B of the 4 values, determining a value of the Z coordinate of a vertex with a number 7 of the 12-sided polygon of the sub-plane 2 as a second large value C of the 4 values, determining a value of the Z coordinate of a vertex with a number 10 of the 12-sided polygon of the sub-plane 2 as a maximum value D of the 4 values, generating 2 values between A E and B E according to a unifrnd function, determining the value of the Z coordinate of the vertex with the sequence number 2 in the 12-sided polygon of the sub-plane 2 to be the smaller of the 2 values, determining the value of the Z coordinate of the vertex with the sequence number 3 in the 12-sided polygon of the sub-plane 2 to be the larger of the 2 values, generating 2 values between B E and C E according to a unifrnd function, determining the value of the Z coordinate of the vertex with the sequence number 5 in the 12-sided polygon of the sub-plane 2 to be the smaller of the 2 values, determining the value of the Z coordinate of the vertex with the sequence number 6 in the 12-sided polygon of the sub-plane 2 to be the larger of the 2 values, randomly generating 2 values between C E and D E according to a unifrnd function, determining the value of the Z coordinate of the vertex with the sequence number 8 in the 12-sided polygon of the sub-plane 2 to be the larger of the 2 values, and determining the value of the Z coordinate of the larger of the vertex with the sequence number 2 in the sub-sided polygon 2 to be the larger of the 2 to be the smaller of the Z coordinate of the sequence number 2 Generating 2 numerical values between D E and A E according to a unifrnd function, determining the value of the Z coordinate of the vertex with the serial number of 11 in the 12-sided polygon of the sub-plane 2 to be the smaller value of the 2 numerical values, determining the value of the Z coordinate of the vertex with the serial number of 12 in the 12-sided polygon of the sub-plane 2 to be the larger value of the 2 numerical values, determining the X coordinate, the Y coordinate and the Z coordinate of the 12 vertices in the 12-sided polygon of the sub-plane 2 through the steps, and sequentially connecting the 12 vertices in sequence of the number sequence and end-to-end to obtain the sub-plane 2;

(3) generating a sub-plane 3 by using a method for generating the sub-plane 2, wherein a distribution range Q of Z coordinates of 12 vertexes of a 12-edge of the sub-plane 3 is defined, and Q is more than or equal to 0.1 and less than or equal to 0.3;

(4) generating a sub-plane 1, generating coordinate values of 6 vertexes of a 6-sided polygon of the sub-plane 1 based on the coordinate values of 6 odd vertexes of a 12-sided polygon of the sub-plane 2, and labeling the 6 vertexes of the 6-sided polygon of the sub-plane 1, the vertexes from number 1 to number 6 of the 6-sided polygon of the sub-plane 1 respectively correspond to the vertexes from number 1, 3, 5, 7, 9, 11 of the 12-sided polygon of the sub-plane 2 one-to-one, the Y coordinates of the vertexes from number 1 to number 6 of the 6-sided polygon of the sub-plane 1 respectively are the Y coordinates of the vertexes from number 1, 3, 5, 7, 9, 11 of the 12-sided polygon of the sub-plane 2 multiplied by 0.8, the X coordinates of the vertexes from number 1 to number 6 of the 6-sided polygon of the sub-plane 1 respectively are the X coordinates of the vertexes 1, 3, 5, 7, 9, 11 of the 12-sided polygon of the sub-plane 2, and defining a distribution range W of the Z coordinates of the 6-sided polygon of the sub-sided polygon 1, 0.825Q 0.875, generating 2 values between 0.825E and 0.875E according to a unifrnd function, determining the value of the Z coordinate of the vertex with number 1 in the 6-sided polygon of the sub-plane 1 to be the smaller value G of the 2 values, determining the value of the Z coordinate of the vertex with number 6 in the 6-sided polygon of the sub-plane 1 to be the larger value K of the 2 values, generating 4 values between G E and K E according to a unifrnd function, determining the value of the Z coordinate of the vertex with number 2 in the 6-sided polygon of the sub-plane 1 to be the smallest value of the 4 values, determining the value of the Z coordinate of the vertex with number 3 in the 6-sided polygon 1 to be the second smallest value of the 4 values, determining the value of the Z coordinate of the vertex with number 4 in the 6-sided polygon of the sub-plane 1 to be the second largest value of the 4 values, determining the value of the Z coordinate of the vertex with number 4 in the 6-sided polygon 1 to be the largest value of the 4 values, sequentially connecting the 6 vertexes according to the sequence of numbers and connecting the vertexes end to obtain the sub-plane 1 according to the X coordinate, the Y coordinate and the Z coordinate of the 6 vertexes in the 6-edge shape of the sub-plane 1 determined in the above steps;

(5) generating a sub-plane 4 by using a method for generating the sub-plane 1, wherein the distribution range of Z coordinates from 6 vertexes of a 6-edge of the sub-plane 4 is defined as T, and T is more than or equal to 0 and less than or equal to 0.1;

3.3: connecting the vertexes of the 24-sided polygon of the three-dimensional center plane of the block stone, the 6-sided polygon of the sub-plane 1, the 12-sided polygon of the sub-plane 2, the 12-sided polygon of the sub-plane 3 and the 6-sided polygon of the sub-plane 4 to form a three-dimensional graph of the block stone, and the steps are as follows:

(1) connecting 24 vertexes of the 24-sided polygon of the three-dimensional center plane of the block stone with 12 vertexes of the 12-sided polygon of the sub-plane 2 to form 36 triangles so as to form a closed curved surface with the boundary line of the three-dimensional center plane of the block stone and the boundary line of the sub-plane 2 as a boundary, wherein the connecting method of the triangles specifically comprises the following steps: forming two triangles by taking the 1 st vertex of the 12-sided polygon of the sub-plane 2 as the first vertex of the triangle, wherein the 24 th vertex and the 1 st vertex of the 24-sided polygon of the three-dimensional central plane of the block stone are respectively the other two vertices of one of the triangles, and the 1 st vertex and the 2 nd vertex of the 24-sided polygon of the three-dimensional central plane of the block stone are respectively the other two vertices of the other triangle; the m-th vertex of the 12-sided polygon of the sub-plane 2 is used as the first vertex of the triangle to form two triangles, m is an integer, m is more than or equal to 2 and less than or equal to 12, the 2m-2 nd vertex and the 2m-1 st vertex of the 24-sided polygon of the three-dimensional center plane of the block stone are respectively the other two vertices of one triangle, and the 2m-1 nd vertex of the 24-sided polygon of the three-dimensional center plane of the block stone are respectively the other two vertices of the other triangle; respectively taking the jth vertex of the 12-sided polygon of the sub-plane 2, the jth +1 vertex of the 12-sided polygon of the sub-plane 2 and the 2 jth vertex of the three-dimensional center plane of the block stone as three vertices of a triangle to form the triangle, wherein j is an integer and is more than or equal to 1 and less than or equal to 11; respectively taking the 12 th vertex of the 12-sided polygon of the sub-plane 2, the 1 st vertex of the 12-sided polygon of the sub-plane 2 and the 24 th vertex of the 24-sided polygon of the three-dimensional center plane of the block stone as three vertices of a triangle to form a triangle;

(2) connecting 24 vertexes of the 24-sided polygon of the three-dimensional center plane of the block stone with 12 vertexes of the 12-sided polygon of the sub-plane 3 to form 36 triangles so as to form a closed curved surface with the boundary line of the three-dimensional center plane of the block stone and the boundary line of the sub-plane 3 as a boundary, wherein the connecting method of the triangles is consistent with the connecting method of the 24 vertexes of the 24-sided polygon of the three-dimensional center plane of the block stone and the 12 vertexes of the 12-sided polygon of the sub-plane 2;

(3) connecting 12 vertexes of the 12-sided polygon of the sub-plane 2 with 6 vertexes of the 6-sided polygon of the sub-plane 1 to form 18 triangles so as to form a closed curved surface with the boundary line of the sub-plane 2 and the boundary line of the sub-plane 1 as a boundary, wherein the connecting method of the triangles specifically comprises the following steps: taking the 1 st vertex of the 6-sided polygon of the sub-plane 1 as the first vertex of the triangle to form two triangles, wherein the 12 th vertex and the 1 st vertex of the 12-sided polygon of the sub-plane 2 are respectively the other two vertices of one of the triangles, and the 1 st vertex and the 2 nd vertex of the 12-sided polygon of the sub-plane 2 are respectively the other two vertices of the other triangle; forming two triangles by taking the kth vertex of the 6-sided polygon of the sub-plane 1 as the first vertex of the triangle, wherein the 2k-2 nd vertex and the 2k-1 st vertex of the 12-sided polygon of the sub-plane 2 are respectively the other two vertices of one of the triangles, the 2k-1 nd vertex and the 2k-1 st vertex of the 12-sided polygon of the sub-plane 2 are respectively the other two vertices of the other triangle, k is an integer, and k is more than or equal to 2 and less than or equal to 6; respectively taking the f-th vertex of the 6-sided polygon of the sub-plane 1, the f + 1-th vertex of the sub-plane 1 and the 2 f-th vertex of the sub-plane 2 as three vertices of a triangle to form a triangle, wherein f is an integer and is more than or equal to 1 and less than or equal to 5; respectively taking the 6 th vertex of the 6-sided polygon of the sub-plane 1, the 1 st vertex of the 6-sided polygon of the sub-plane 1 and the 12 th vertex of the 12-sided polygon of the sub-plane 2 as three vertices of a triangle to form a triangle;

(4) connecting 12 vertexes of the 12-sided polygon of the sub-plane 3 with 6 vertexes of the 6-sided polygon of the sub-plane 4 to form 18 triangles so as to form a closed curved surface with the boundary line of the sub-plane 3 and the boundary line of the sub-plane 4 as a boundary, wherein the connecting method of the triangles is consistent with the connecting method of the 12 vertexes of the 12-sided polygon of the sub-plane 2 with the 6 vertexes of the 6-sided polygon of the sub-plane 1;

3.4: constructing a three-dimensional graph of the block stone of the two-dimensional projection of the block stone acquired in the step 1 according to a method of 3.1 to 3.3;

further, the step 5 specifically includes:

5.1: endowing the block stone with a space coordinate and a space direction of a three-dimensional graph;

5.2: embedding the three-dimensional graphics of the rock block into a three-dimensional slope model of the engineering area, wherein the three-dimensional graphics of the rock block cannot be intersected with the surface of the three-dimensional slope model of the engineering area and cannot be intersected with the embedded three-dimensional graphics of the rock block, otherwise, endowing a group of new space coordinates and space orientations to the three-dimensional graphics of the rock block until the three-dimensional graphics of the rock block are not intersected with the surface of the three-dimensional slope model of the engineering area and the embedded three-dimensional graphics of the rock block;

5.3: repeating the steps to embed the three-dimensional graphics of the block stone, which are acquired in the step 1 and subjected to the two-dimensional projection, into the three-dimensional slope model of the engineering area;

5.4: and applying Boolean operation to the embedded three-dimensional graph of the block stone and the three-dimensional slope model of the engineering area to obtain the heterogeneous slope model filled with soil and stone.

Further, in step 1, a two-dimensional projection of the block stone of the engineering area is acquired by the aggregate image analysis system AIMS 2.

Further, the slope image in the step 4 is obtained through unmanned aerial vehicle aerial photography.

Further, in step 7, the non-homogeneous slope model is analyzed and calculated by adopting a strength reduction theory in finite element numerical analysis software ABAQUS to obtain a slope sliding displacement cloud picture, a slope sliding plastic strain cloud picture and a slope stability safety coefficient

The slope stability analysis report.

Further, the safety factor is stabilized on the side slope

The method is obtained by simulating the instability process of the side slope by using the reduction of the shear strength of the soil body, and specifically comprises the following steps:

defining a field variable as a strength reduction coefficient in ABAQUS finite element analysis software

Soil mass cohesive force in reduced soil mass shear strength parameters

And angle of internal friction

Respectively as follows:

wherein c is the soil mass cohesive force in the soil mass shear strength parameter,

is the internal friction angle in the soil shear strength parameter and the soil mass cohesive force in the doubled and subtracted soil shear strength parameter

And angle of internal friction

Performing finite element analysis, and continuously increasing in the whole process

Until the soil mass balance is damaged critically, the strength reduction coefficient

Stabilizing safety factor for side slope

。

Compared with the prior art, the invention has the following beneficial effects:

(1) according to the method, the aggregate image analysis system AIMS2 is used for collecting the two-dimensional projections of the stones in batches, 24-edge shapes of the two-dimensional projections of the stones are generated, the 24-edge shapes of the two-dimensional projections of the stones are used for generating three-dimensional figures of the stones with any shapes, the stones are enabled to be more in line with the actual situation, the three-dimensional figures of the stones are endowed with space coordinates and directions and are put into the three-dimensional side slope model, the three-dimensional side slope model is generated through the three-dimensional reconstruction technology by means of images shot by the unmanned aerial vehicle, the actual heterogeneous side slope model is constructed, and the accuracy of the stability analysis of the heterogeneous side slope is effectively improved.

(2) According to the method, a slope stability analysis report containing a slope sliding displacement cloud picture, a plastic strain cloud picture when a slope slides and a slope stability safety coefficient is obtained by analyzing and calculating the heterogeneous slope model constructed by the method by adopting a strength reduction theory in finite element numerical analysis software ABAQUS, and the slope stability analysis report is richer in content, provides reliable guarantee for engineering numerical simulation and provides a basis for reasonable support measures and suggestions for the slope which is likely to generate instability.

Drawings

FIG. 1 is a schematic flow chart of the present invention.

Fig. 2 is a schematic diagram of a flow chart of acquiring a 24-edge coordinate database of two-dimensional block stone projection according to the present invention.

FIG. 3 is a schematic diagram of a three-dimensional graph obtaining process of the present invention.

FIG. 4 is a heterogeneous slope model of the present invention.

Fig. 5 is a schematic diagram of the sliding displacement when the heterogeneous slope is unstable.

FIG. 6 is a schematic diagram of the plastic region development when the heterogeneous slope is unstable.

Detailed Description

The present invention will be described in detail below with reference to examples and the accompanying drawings. It should be noted that the embodiments and features of the embodiments may be combined with each other without conflict.

Referring to fig. 1, the method for analyzing the stability of the heterogeneous slope of the invention comprises the following steps:

step 1: the method comprises the steps of collecting two-dimensional projections of the stones in the engineering area, preferably collecting the two-dimensional projections of the stones in batches by using an Aggregate Image analysis System AIMS2(Aggregate Image System 2), namely placing the stones in the Aggregate Image analysis System through a tray filled with samples, wherein the samples are the stones, and the Aggregate Image analysis System automatically obtains images of each stone and analyzes the shape characteristics of the stone, and finally quickly obtains the appearance characteristics of the stones, and provides objective characteristics of the edges, the structures and the surface textures of the stones.

as shown in fig. 2, step 2 specifically includes:

2.2: constructing 24-sided polygons of the two-dimensional projections of the block stones by using external contour lines of the two-dimensional projections of the block stones, and establishing a database for recording coordinates of each vertex of the 24-sided polygons of the two-dimensional projections of the block stones, wherein the coordinates of each vertex comprise an X abscissa_iAnd ordinate Y_iWherein i is the vertex serial number of the 24-sided polygon of the two-dimensional projection of the block stone, i is an integer, i is within a value range of 1 and not more than 24, the serial numbers of the vertices are sequentially recorded in an anticlockwise direction by taking the positive direction of an X axis as the starting direction, and the method for determining the 24-sided polygon of the two-dimensional projection of the block stone and the vertex coordinates thereof is as follows:

as shown in fig. 3, step 3 specifically includes:

(1) defining the ratio E of the height of the stone in the Z direction to the width of the stone in the X direction, wherein E is more than 0 and less than 1, and preferably, the value of E is 0.8;

finally forming a three-dimensional figure of the block stone by the closed curved surface formed by the method;

3.4: constructing a three-dimensional graph of the block stone of the two-dimensional projection of the block stone acquired in the step 1 according to a method of 3.1 to 3.3, wherein the three-dimensional graph of the block stone can be in any shape;

and 4, step 4: acquiring a slope image of the engineering area, and generating a three-dimensional slope model of the engineering area by a three-dimensional reconstruction technology, preferably, acquiring the slope image by unmanned aerial vehicle aerial photography;

and 5: randomly embedding the three-dimensional graphics of the rock block into a three-dimensional slope model of the engineering area, and establishing a heterogeneous slope model;

the step 5 specifically comprises the following steps:

the spatial orientation may be set to a random value, or may be defined as a specific value or a specific value range according to research needs. The spatial orientation is determined by sequentially rotating the three-dimensional pattern of the block stone along the x-axis

Arc, rotation along y-axis

Arc and rotation along z-axis

The radian is obtained in a specific mode:

three-dimensional graphic head of rock blockRotate around the x-axis

After radian, the centroid coordinate of the three-dimensional graph of the block stone is changed into:

when the three-dimensional figure of the block stone rotates around the x axis

Rotate around the y-axis again after radian

when the three-dimensional figure of the block stone rotates around the x axis

After arc, rotate about the y-axis

After radian, the three-dimensional graph of the block stone rotates around the z axis

Radian and final centroid coordinate change of the three-dimensional graph of the stone are as follows:

wherein the content of the first and second substances,

the x-coordinate of the centroid of the three-dimensional figure of the block stone after rotation,

the y coordinate of the centroid of the three-dimensional graph of the block stone after rotation;

the z coordinate of the centroid of the three-dimensional graph of the block stone after rotation;

the x coordinate of the centroid of the three-dimensional graph of the block stone before rotation;

the y coordinate of the centroid of the three-dimensional graph of the block stone before rotation;

the z coordinate of the centroid of the three-dimensional graph of the block stone before rotation;

5.2: embedding the three-dimensional graphics of the rock block into a three-dimensional slope model of the engineering area, wherein the three-dimensional graphics of the rock block cannot be intersected with the surface of the three-dimensional slope model of the engineering area and cannot be intersected with the embedded three-dimensional graphics of the rock block, otherwise, endowing a group of new space coordinates and space orientations to the three-dimensional graphics of the rock block until the three-dimensional graphics of the rock block are not intersected with the surface of the three-dimensional slope model of the engineering area and the embedded three-dimensional graphics of the rock block; when intersection occurs, optionally, adjusting the rigidity of the surface of the three-dimensional graphics of the rock block so as to limit the contact of the three-dimensional graphics of the rock block and the three-dimensional graphics of the adjacent rock block, continuously generating mutual exclusion in the repeated modeling process of the three-dimensional graphics of the rock block adjacent to each other by adjusting the rigidity of the surface of the three-dimensional graphics of the rock block so that the intersection of the three-dimensional graphics of the rock block is not caused any more when the surface distance is continuously increased, optionally, selecting two-dimensional projection of the rock block with smaller grain size of the rock block to construct the three-dimensional graphics of the rock block, and thus reducing the probability of intersection of the three-dimensional graphics of the rock block;

5.3: the step 1 is repeated to embed the three-dimensional figures of the block stones, which are acquired in the step 1 and subjected to two-dimensional projection, into the three-dimensional side slope model of the engineering area, and the three-dimensional figures of the block stones can be randomly distributed in the engineering area, so that the three-dimensional side slope model is closer to the real situation;

5.4: and (3) applying Boolean operation to the embedded three-dimensional graph of the rock block and the three-dimensional slope model of the engineering area to obtain the heterogeneous slope model filled with the soil and the rock, as shown in figure 4.

and 7: analyzing the heterogeneous slope model according to finite element strength reduction method to generate slope stability analysis report rapidly, preferably, analyzing and calculating the heterogeneous slope model by adopting strength reduction theory in finite element numerical analysis software ABAQUS to obtain slope sliding displacement cloud picture, plastic strain cloud picture when slope slides and slope stability safety factor

The slope stability analysis report ofCoefficient of performance

Soil mass cohesive force in reduced soil mass shear strength parameters

And angle of internal friction

Respectively as follows:

And angle of internal friction

Until the soil mass balance is damaged critically, the soil mass balance is damagedReduction factor of strength

Stabilizing safety factor for side slope

；

As shown in fig. 5, intensity reduction coefficiency for soil mass

When the slope is reduced to a certain degree, the slope slides, the shape of the sliding surface of the heterogeneous slope is relatively complex relative to the shape of the soil slope, mainly because the rock blocks exist in the heterogeneous slope, the position of the sliding surface can be changed, so that the deformation characteristic of the slope is changed, and the position of the changed sliding surface is closely related to the spatial distribution of the rock blocks in the slope. Therefore, on-site point arrangement monitoring is carried out on the peak generating the maximum displacement in the side slope, and if the stability is insufficient after monitoring, reinforcement processing is carried out on the side slope so as to improve the stability of the side slope.

As shown in fig. 6, when the slope starts to slide, a relatively obvious plastic region exists near the sliding surface, and the distribution of the plastic region is matched with that of the sliding surface, which shows that the sliding surface form of the heterogeneous slope is closely related to the form of the formed through plastic region. Due to the existence of the rock blocks, the nonuniformity of the internal structure of the side slope is enhanced, and under certain conditions, the internal sub-area of the side slope reaches a plastic state, but due to the existence of the rock blocks, a coherent plastic area is difficult to form, so that the sliding of the side slope is restrained. Therefore, the influence of the spatial distribution of the rock blocks on the stability of the side slope actually influences the form of the plastic zone and further influences the form of the slip surface of the side slope. The stability of the side slope can be influenced by the spatial distribution characteristics of the stones in the soil-stone side slope, and whether the side slope body can generate integral sliding instability to cause instability deformation of the roadbed and the cutting is judged according to the plastic deformation development trend of the side slope under the natural action.

The foregoing is a more detailed description of the invention in connection with specific preferred embodiments and it is not intended that the invention be limited to these specific details. For those skilled in the art to which the invention pertains, several simple deductions and substitutions can be made without departing from the spirit of the invention, and all shall be considered as belonging to the protection scope of the invention.

Claims

1. A heterogeneous slope stability analysis method is characterized by comprising the following steps:

2. The heterogeneous slope stability analysis method according to claim 1, wherein the step 2 specifically comprises:

2.2: constructing 24-sided polygons of the two-dimensional projections of the block stones by using external contour lines of the two-dimensional projections of the block stones, and establishing a database for recording coordinates of each vertex of the 24-sided polygons of the two-dimensional projections of the block stones, wherein the coordinates of each vertex comprise an X abscissa_iAnd ordinate Y_iWhere i is the apex of the 24-sided polygon of the two-dimensional projection of the block stoneThe method for determining the 24-edge shape of the two-dimensional projection of the block stone and the vertex coordinates thereof comprises the following steps of sequentially recording the serial numbers of vertexes in a counterclockwise direction by taking the positive direction of an X axis as the starting direction, wherein i is an integer and is not more than 1 and not more than 24, and the method for determining the 24-edge shape of the two-dimensional projection of the block stone and the vertex coordinates thereof comprises the following steps:

when the number of nodes on the external contour line of the two-dimensional block stone projection is 24, 24 polygons of the two-dimensional block stone projection correspond to the external contour line of the two-dimensional block stone projection, 24 vertexes of the 24 polygons of the two-dimensional block stone projection are 24 nodes on the external contour line of the two-dimensional block stone projection respectively, and the coordinates of the 24 vertexes of the 24 polygons of the two-dimensional block stone projection correspond to the coordinates of the 24 nodes on the external contour line of the two-dimensional block stone projection one by one;

when the number of nodes on the external contour line of the two-dimensional projection of the block stone is more than 24, calculating the length of each edge on the external contour line of the two-dimensional projection of the block stone, finding the edge with the shortest length and the shorter edge of the two adjacent edges of the edge, deleting the node shared by the shorter edge and the edge with the shortest length, and connecting the other node left by the shorter edge and the other node left by the edge with the shortest length to form the external contour line of the two-dimensional projection of the new block stone;

when the number of nodes on the external contour line of the two-dimensional projection of the block stone is less than 24, calculating the length of each edge on the external contour line of the two-dimensional projection of the block stone, finding the edge with the longest length, respectively connecting the midpoint of the edge with the longest length with two nodes of the edge with the longest length, and not connecting the two nodes of the edge with the longest length to form a new external contour line of the two-dimensional projection of the block stone;

2.3: and (3) constructing the 24-sided polygon of the two-dimensional projection of the block stone acquired in the step (1) by using the methods of 2.1 to 2.2, and establishing a database for recording each vertex coordinate of the 24-sided polygon of each two-dimensional projection of the block stone.

3. The heterogeneous slope stability analysis method according to claim 2, wherein the step 3 specifically comprises:

(1) connecting 24 vertexes of the 24-sided polygon of the three-dimensional center plane of the block stone with 12 vertexes of the 12-sided polygon of the sub-plane 2 to form 36 triangles so as to form a closed curved surface with the boundary line of the three-dimensional center plane of the block stone and the boundary line of the sub-plane 2 as a boundary, wherein the connecting method of the triangles specifically comprises the following steps: forming two triangles by taking the 1 st vertex of the 12-sided polygon of the sub-plane 2 as the first vertex of the triangle, wherein the 24 th vertex and the 1 st vertex of the 24-sided polygon of the three-dimensional central plane of the block stone are respectively the other two vertices of one of the triangles, and the 1 st vertex and the 2 nd vertex of the 24-sided polygon of the three-dimensional central plane of the block stone are respectively the other two vertices of the other triangle; taking the mth vertex of the 12-sided polygon of the sub-plane 2 as the first vertex of the triangle to form two triangles, wherein m is an integer, m is more than or equal to 2 and less than or equal to 12, the 2m-2 nd vertex and the 2m-1 st vertex of the 24-sided polygon of the three-dimensional center plane of the block stone are respectively the other two vertices of one of the triangles, and the 2m-1 nd vertex of the 24-sided polygon of the three-dimensional center plane of the block stone are the other two vertices of the other triangle; respectively taking the jth vertex of the 12-sided polygon of the sub-plane 2, the jth +1 vertex of the 12-sided polygon of the sub-plane 2 and the 2 jth vertex of the three-dimensional center plane of the block stone as three vertices of a triangle to form the triangle, wherein j is an integer and is more than or equal to 1 and less than or equal to 11; respectively taking the 12 th vertex of the 12-sided polygon of the sub-plane 2, the 1 st vertex of the 12-sided polygon of the sub-plane 2 and the 24 th vertex of the 24-sided polygon of the three-dimensional center plane of the block stone as three vertices of a triangle to form a triangle;

3.4: and (3) constructing the three-dimensional graph of the two-dimensional projection of the block stone acquired in the step (1) according to a method of 3.1 to 3.3.

4. The heterogeneous slope stability analysis method according to claim 3, wherein the step 5 specifically comprises:

5. The heterogeneous slope stability analysis method according to claim 4, wherein the step 1 of acquiring a two-dimensional projection of the block stone of the engineering area by an aggregate image analysis system AIMS 2.

6. The heterogeneous slope stability analysis method of claim 5, wherein the slope image of step 4 is obtained by unmanned aerial vehicle aerial photography.

7. The method for analyzing stability of heterogeneous slope according to claim 6, wherein in step 7, the intensity reduction theory is applied in the finite element numerical analysis software ABAQUS to analyze and calculate the model of heterogeneous slope, so as to obtain the cloud map of slope sliding displacement, the cloud map of plastic strain when the slope slides, and the safety factor of slope stability