CN-120953533-B - Vine cube slope restoration effect dynamic simulation method

CN120953533BCN 120953533 BCN120953533 BCN 120953533BCN-120953533-B

Abstract

The invention relates to the technical field of three-dimensional modeling, in particular to a vine cube slope restoration effect dynamic simulation method which comprises the steps of reconstructing slope grids according to slope abrupt change and slope direction based on remote sensing images and point clouds, extracting three-dimensional coordinates and contact layout of anchor points, analyzing stress and density marking disturbance grades, screening continuous deformation nodes to form deformation paths, and drawing tension offset tracks and intersection sequences. According to the invention, the gray level and the point cloud elevation difference of the remote sensing image are fused, the topographic fold line is extracted for grid reconstruction according to gradient abrupt change and gradient direction consistency, topographic features and engineering layout are matched, three-dimensional coordinates of anchor points are extracted, space topology layout is generated in a grouping mode according to layout directions, a multi-dimensional interference level system is constructed according to a main tensile stress and node density discrimination disturbance unit, a continuous deformation path is screened, a tension deviation track is drawn, and an intersection section is indexed, so that multistage restoration evolution visualization is realized.

Inventors

Xia Sihong
ZENG WEI
XIE WUPING
CHEN LIANG
NING CHEN
ZHANG YONGHAI
LIU BIN
CHEN SIYU
HE YIJUN
ZHANG DONGDONG
WANG HUIMIN
ZHU JIAXING
YANG SHIZHEN

Assignees

湖南中核建设工程有限公司
湖南省生态地质调查监测所
中南林业科技大学

Dates

Publication Date: 20260505
Application Date: 20250806

Claims (6)

1. The dynamic simulation method for the vine cube slope restoration effect is characterized by comprising the following steps of: s1, acquiring remote sensing images and point cloud data of a vine cube structure paving area, sorting grid nodes according to elevation gradients, extracting terrain broken lines to reconstruct slope grid points, dividing a uniform slope area, and generating a slope structure partition unit matrix; the step of obtaining the slope structure partition unit matrix comprises the following steps: S101, acquiring a remote sensing image gray value matrix and a three-dimensional point cloud elevation data set, comparing an image gray value at a position corresponding to each grid node with a point cloud elevation value pixel by pixel, calculating a numerical value difference absolute value of the two, establishing a grid node difference value distribution thermodynamic diagram, screening node positions of which the difference value exceeds a terrain mutation threshold, and generating a terrain difference coefficient matrix; S102, calculating gradient change rates of adjacent nodes along the longitudinal direction based on the terrain difference coefficient matrix, detecting mutation areas with gradient change rates exceeding a continuously set number of nodes by adopting a sliding window method, recording mutation point coordinate sequences, forming a broken line track by connecting adjacent mutation points, and optimizing broken line curvature continuity through cubic spline interpolation to obtain a terrain broken line characteristic set; s103, calling the terrain broken line feature set as a grid division boundary, carrying out topology reconstruction on an original grid, calculating elevation gradient vectors of all nodes in a new grid unit, counting deviation values of gradient vector azimuth angles in the units, reserving grid units with deviation values smaller than a set angle threshold, merging adjacent units to form a continuous area, and outputting a slope structure partition unit matrix; s2, based on the slope structure partition unit matrix, extracting three-dimensional coordinate data of an anchoring point of the rattan cube structure, calculating the slope contact area and the overlap length, grouping and marking numbers according to the included angles between the arrangement direction of the rattan cube structure components and the slope vector of the slope, and generating a component slope coupling contact layout map; S3, calling the component slope coupling contact layout map, acquiring the distribution of the main tensile stress direction and the contact normal stress of the nodes of the contact area of the rattan cube structural component, calculating the included angle between the direction vector of the rattan cube structural component and the main stress direction of the nodes, judging the consistency distribution, counting the contact rate and the distribution frequency in the unit grid, marking the disturbance level of the area with the included angle deviating from the set limit, and outputting a slope response disturbance level classification map layer; S4, calling a disturbance grade marking area in the slope response disturbance grade classification layer, acquiring the elevation change rate of slope grid nodes and a horizontal displacement track in a continuous period, and extracting the nodes according to the consistency of deformation directions to form a slope structure deformation trend path group; S5, based on the slope structure deformation trend path group, recording the change data of the structure tension direction of continuous deformation nodes in a time sequence, drawing a tension deviation track graph according to the tension direction vector conversion sequence in the path, marking the corresponding node coordinates of each deviation point in the path in parallel with the connection sequence of rattan cube structural members, marking out tension turning nodes and path intersection sections, sorting the whole-area path sequence and outputting a restoration-area multi-section structure response dynamic sequence graph; the repair area multistage structure response dynamic sequence diagram comprises a tension direction vector track diagram, a turning node space coordinate set and a path intersection section topological index.
2. The vine cube slope restoration effect dynamic simulation method of claim 1, wherein the slope structure partition unit matrix comprises a reconstruction grid boundary, a slope consistency identifier and a unit topological relation, the component slope coupling contact layout map is specifically an anchor point three-dimensional coordinate set, a contact area parameter matrix and a component grouping coding table, the slope response interference level classification map layer comprises a main stress direction vector field, a contact rate density distribution cloud map and a disturbance unit boundary coordinate set, and the slope structure deformation trend path group is specifically an elevation change rate gradient map, a displacement direction sequence matrix and a deformation node topological network.
3. The method for dynamically simulating the vine cube slope restoration effect according to claim 1, wherein the step of obtaining the component slope coupling contact layout map is as follows: S201, calling the slope structure partition unit matrix, extracting a three-dimensional coordinate data set of anchoring points of a rattan cube structure, performing space matching on coordinates of each anchoring point and corresponding grid nodes, calculating an elevation curved surface fitting plane of all nodes in a grid unit where the anchoring points are located, and counting the number of contact points of projection of the anchoring points to the fitting plane to generate a contact area vector set; s202, extracting adjacent anchor point distance data along the boundary of a rattan cube structural member based on the contact area vector set, measuring the length difference value of the overlapping area at the end part of the rattan cube structural member, recording an overlapping length change curve according to the numbering sequence of the rattan cube structural member, establishing a rattan cube structural member connection sequence index table, and outputting a connection sequence; S203, calculating the cosine value of the space included angle between the arrangement direction vector of each rattan cube structural component and the corresponding slope direction vector by combining the connection sequence, classifying rattan cube structural components with cosine values larger than a set threshold value into the same group, carrying out continuous area numbering mapping on the rattan cube structural components in the same group, integrating the contact area and the connection sequence data, and generating a component slope coupling contact layout map.
4. The method for dynamically simulating the vine cube slope restoration effect according to claim 1, wherein the step of acquiring the slope response interference level classification layer is as follows: S301, calling a component slope coupling contact layout map, extracting a main tensile stress direction vector and a positive stress vector value of a contact area node, calculating a cosine value of a space included angle between the main tensile stress direction of each node and a corresponding rattan cube structural component direction vector, and counting cosine value distribution frequency to generate a force direction consistency coefficient set; S302, based on the force direction consistency coefficient set, counting node density values according to grid units, calculating the number of contact nodes and positive stress distribution frequency in a unit area, taking the product of the number of the contact nodes and the positive stress frequency as a contact rate index, establishing a grid contact rate distribution matrix, and outputting a contact rate gradient field; S303, calculating the absolute value of the difference value of the force direction consistency coefficients between adjacent grid cells by combining the contact rate gradient fields, screening the areas with the difference value exceeding the stress disturbance threshold, carrying out topological combination on the continuous exceeding areas, carrying out grading assignment according to the disturbance intensity, and generating a slope response disturbance grade classification layer.
5. The method for dynamically simulating the vine cube slope restoration effect according to claim 1, wherein the step of obtaining the slope structure deformation trend path group is as follows: S401, calling the slope response interference level classification layer, extracting time sequence data of the elevation change rate of grid nodes in a disturbance boundary, calculating the elevation increment ratio of the nodes in adjacent time periods, establishing an elevation change trend vector of the nodes, and generating an elevation dynamic trend set; S402, based on the elevation dynamic trend set, horizontal displacement track coordinates of the same node in a continuous period are obtained, an included angle cosine value of displacement direction vectors of adjacent periods is calculated, the number of continuous periods with the cosine value larger than a displacement homodromous threshold is counted, and a horizontal displacement homodromous sequence is output; And S403, screening nodes with the elevation increment ratio change direction and the horizontal displacement direction cosine value synchronously increasing by combining the horizontal displacement isotropic sequence, performing spatial cluster analysis on the nodes meeting the synchronous condition, connecting adjacent nodes to form a continuous path, and generating a slope structure deformation trend path group.
6. The method for dynamically simulating the vine cube slope restoration effect according to claim 1, wherein the step of obtaining the restoration area multi-section structure response dynamic sequence diagram is as follows: s501, calling the slope structure deformation trend path group, extracting time sequence tension direction vectors of continuous deformation nodes, calculating the variation of included angles between adjacent time interval vectors according to time sequence, drawing tension direction deviation track line diagrams in each path, and generating a tension deviation track set; S502, based on the tension deviation track set, associating node coordinates of a path where each deviation point is located with a component connection sequence number, establishing a mapping relation table of space positions of the deviation points and a component sequence, marking node positions of abrupt change of tension direction exceeding a set angle, and outputting tension turning node indexes; S503, counting the turning node density in the intersection region of the path by combining the tension turning node indexes, dividing the intersection section level according to a density threshold value, integrating the time-space evolution data of the full-area path, and generating a restoration area multistage structure response dynamic sequence diagram.

Description

Vine cube slope restoration effect dynamic simulation method Technical Field The invention relates to the technical field of three-dimensional modeling, in particular to a vine cube slope restoration effect dynamic simulation method. Background The technical field of three-dimensional modeling comprises various technical methods for digitally expressing and reconstructing real world objects or scenes. The core content of the method is that a three-dimensional virtual model with spatial depth, geometric structure and visual reality is generated through means of data acquisition, spatial reconstruction, graphic rendering, dynamic updating and the like. Input information such as point cloud data, image measurement data, terrain scanning data and the like is widely adopted, and modeling and simulation processing are carried out by combining a graph geometric algorithm. The three-dimensional modeling is widely applied to a plurality of fields such as engineering construction, environment monitoring, virtual reality, disaster assessment, geographic information systems and the like, can realize digital restoration, visual analysis and dynamic demonstration of real environment characteristics, and has higher technical systematicness and application complexity. The rattan cube slope restoration effect dynamic simulation method is used for constructing a dynamic evolution process model of the restoration effect under different time sequences and physical actions in a slope support restoration process by using a rattan cube structure. The method mainly relates to extracting original slope structure data through remote sensing image interpretation and slope section measurement, carrying out space reconstruction by combining with vine cube component geometric parameters and layout modes, setting different natural environment parameters such as rainfall, wind erosion, gravity effect and the like by utilizing a three-dimensional modeling engine, and carrying out simulation on time sequence response changes. And the overall deformation trend of the side slope is defined in a sectional manner by adopting a grid division mode based on terrain data, a dynamic evolution rule is established by setting fixed boundary conditions and material mechanical parameters, and finally, the simulation reconstruction of geometric structure evolution, deformation response paths and structural stability trend in the vine cube side slope restoration process is completed. In the prior art, grid division is performed by relying on single terrain scanning data and fixed boundary conditions, so that the matching degree of slope structure partitions and terrain features is insufficient, the influence of complex slope changes on engineering layout is difficult to reflect, static contact area calculation and homogenized material parameter setting are adopted in traditional component mechanical analysis, the local stress concentration phenomenon caused by component arrangement direction and slope included angle change cannot be captured, deformation trend prediction is based on elevation change analysis of discrete nodes, cooperative observation of horizontal displacement track and tension direction deviation is lacking, hysteresis and fragmentation features exist in deformation path identification, and a dynamic simulation result is difficult to present continuous evolution rule of structural response in space-time dimension. Disclosure of Invention In order to solve the technical problems in the prior art, the embodiment of the invention provides a vine cube slope restoration effect dynamic simulation method. The technical scheme is as follows: a vine cube slope restoration effect dynamic simulation method comprises the following steps: s1, acquiring remote sensing images and point cloud data of a vine cube structure paving area, sorting grid nodes according to elevation gradients, extracting terrain broken lines to reconstruct slope grid points, dividing a uniform slope area, and generating a slope structure partition unit matrix; S2, based on the slope structure partition unit matrix, extracting three-dimensional coordinate data of an anchoring point of a rattan cube structure, calculating the contact area and the overlap length of the slope, grouping and marking numbers according to the included angles between the component layout direction and the slope vector of the slope, and generating a component slope coupling contact layout map; S3, calling the component slope coupling contact layout map, acquiring the distribution of the main tensile stress direction and the contact normal stress of the component contact area node, calculating the included angle between the component direction vector and the main node stress direction, judging the consistency distribution, counting the contact rate and the distribution frequency in the unit grid, marking the disturbance level of the area with the included angle deviating from the set limit, and outputting