A three-dimensional modeling method of real shape of block stone in earth-rock mixed filling slope

Abstract

This invention discloses a three-dimensional modeling method for the true shape of boulders in a soil-rock fill slope, comprising: Step 1, acquiring the two-dimensional projections of the boulders of the actual soil-rock fill slope onto three orthogonal planes; Step 2, obtaining the coordinates of the two-dimensional projection contour points, calculating the centroid of the contour line, moving the centroid to the origin of the coordinate system, and recalculating the coordinates of the contour points; Step 3, segmenting the contour line into spatial regions based on the projections of the boulders on the XOZ and YOZ planes, and randomly reconstructing a contour line; Step 4, performing compatibility analysis on the XOY plane and correcting the reconstructed contour line; Step 5, randomly constructing several contour lines within each spatial region; Step 6, generating a three-dimensional geometric model with a true shape based on the contour lines. The method of this invention constructs a three-dimensional realistic model based on the true shape of the boulders, resulting in higher accuracy in subsequent finite element analysis of the stability of the soil-rock fill slope from a heterogeneous perspective.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of engineering geological survey analysis, and more particularly to a three-dimensional modeling method for real shape of block stone in soil and rock mixed filling slope. BACKGROUND

[0002] The soil and rock mixture is composed of block stone and soil. With the overall development of China's highway construction towards the western mountainous areas, the good roadbed material specified in the specification is not easy to obtain, and it is inevitable to use the soil and rock mixture obtained by excavating the mountain as the roadbed filler to fill the roadbed or foundation when building the road according to the local conditions. The soil and rock mixed filling slope is mainly composed of soil and block stone, and the structure of such non-homogeneous granular material leads to more complex engineering properties of such slope compared with general soil slope and rock slope. However, from the perspective of microstructure, the shape, content and position distribution of block stone in different soil and rock mixtures will affect the numerical analysis.

[0003] In the prior art, the block stone soil slope micro model is randomly generated based on the existing block stone shape database, and the two-dimensional soil and rock mixture model is analyzed. In the prior art, a three-dimensional discrete element model of the soil and rock mixture slope microstructure is established by using the developed irregular block stone and soil and rock mixture three-dimensional discrete element modeling method. However, the block stone in the prior art is not constructed based on the real shape, and the finite element simulation of the irregular block stone will have some differences from the real block stone in the slope, thereby resulting in insufficient accuracy when analyzing the stability of the soil and rock mixed filling slope from the non-homogeneous perspective by using the finite element method subsequently. SUMMARY

[0004] In view of at least one defect or improvement demand of the prior art mentioned in the background section, the present application provides a three-dimensional modeling method for real shape of block stone in soil and rock mixed filling slope, to improve the simulation accuracy when simulating the block stone by using the finite element method, thereby improving the accuracy when analyzing the stability of the soil and rock mixed filling slope from the non-homogeneous perspective by using the finite element method subsequently.

[0005] To achieve the above-mentioned purpose, the present application provides a three-dimensional modeling method for real shape of block stone in soil and rock mixed filling slope, comprising:

[0006] S1, collecting two-dimensional projections of block stone in soil and rock mixed filling slope on three mutually orthogonal planes;

[0007] S2, obtaining the coordinates of the contour points of each two-dimensional projection and calculating the centroid of the contour line of each two-dimensional projection, respectively translating the contour line of each two-dimensional projection so that the respective centroids are moved to the coordinate origin, and respectively recalculating the new coordinates of the contour points of each two-dimensional projection;

[0008] S3. Based on the two-dimensional projection surfaces of the stone on the first and second orthogonal planes, where the centroid has been moved to the origin of the coordinates, the outlines of the two-dimensional projections on the two orthogonal planes are segmented, and in any spatial region divided by the two orthogonal planes, a outline perpendicular to the third orthogonal plane is randomly reconstructed according to the outline reconstruction logic.

[0009] S4. Based on the projection of the reconstructed contour line onto the third orthogonal plane and the contour line correction logic between the two-dimensional projection plane of the stone block on the third orthogonal plane with its centroid moved to the origin of the coordinates, perform compatibility analysis on the third orthogonal plane and correct the reconstructed contour line.

[0010] S5. Based on the contour line reconstruction logic of step S3 and the contour line correction logic of step S4, reconstruct and correct several contour lines in each spatial region divided by the first orthogonal plane and the second orthogonal plane respectively.

[0011] S6. Based on all the contour lines reconstructed and corrected within each spatial region divided by the first orthogonal plane and the second orthogonal plane, construct a three-dimensional geometric model of the true shape of the stone.

[0012] Furthermore, the two-dimensional projections of the boulders in the soil-rock fill slope on three mutually orthogonal planes include:

[0013] S11. Place the boulders collected from the soil-rock mixed fill slope into a cubic mold and fix their position;

[0014] S12. Scan the three mutually perpendicular mold plates of the cube mold using an image analysis system to obtain a two-dimensional projection of the stone block on the three mutually perpendicular mold plates.

[0015] Further, the steps of obtaining the coordinates of the contour points of each two-dimensional projection and calculating the centroid of the contour line of each two-dimensional projection, translating the contour line of each two-dimensional projection so that its centroid moves to the origin of the coordinate system, and recalculating the new coordinates of the contour points of each two-dimensional projection include:

[0016] S21. Extract the coordinates of the contour points of the two-dimensional projection of the stone block on the first orthogonal plane;

[0017] S22. Calculate the centroid coordinates of the outline of the two-dimensional projection of the stone block on the first orthogonal plane, translate the outline of the two-dimensional projection so that its centroid moves to the origin of the coordinate system, and perform translation calculations on the coordinates of the outline points of the two-dimensional projection to obtain the new coordinates of the outline points of the two-dimensional projection.

[0018] S23. Using the methods in steps S21-S22, obtain the new coordinates of the outline points of the two-dimensional projection of the stone on the second orthogonal plane and the third orthogonal plane, respectively.

[0019] Further, the step of segmenting the contour lines of the two-dimensional projections on the first and second orthogonal planes, where the centroids of the stone have been moved to the origin, and randomly reconstructing a contour line perpendicular to the third orthogonal plane in any spatial region divided by the two orthogonal planes according to the contour line reconstruction logic, includes:

[0020] S31. Using the intersection axis of the first orthogonal plane and the second orthogonal plane as the dividing line, the outline of the two-dimensional projection on the two orthogonal planes is divided into four segments; the first orthogonal plane and the second orthogonal plane divide the space into four regions;

[0021] S32. Taking the two intersection points of each segment of the two-dimensional projection contour line and the above-mentioned intersecting axis as the starting point and the ending point respectively, insert the same preset number of points on each segment of the two-dimensional projection contour line.

[0022] S33. Randomly generate a first angle that is acute, and determine the first angle as the angle between the first contour line reconstructed in the first region and the projection of the first two-dimensional projection contour line on the third orthogonal plane.

[0023] S34. Based on the first random number between 0 and 1, construct the first relationship between the coordinates of the first point on the first contour line and the coordinates of the corresponding first points on the two adjacent two-dimensional projected contour lines.

[0024] S35. Based on a second random number between 0 and 1, construct a second relationship between the coordinates of the first point on the first contour line and the coordinates of the corresponding first points on two adjacent two-dimensional projected contour lines.

[0025] S36. Based on the first angle, construct a third relational expression for the linear equation of the projection of the first contour line onto the third orthogonal plane;

[0026] S37. Combine the first to third relations to obtain the coordinates of the first point on the reconstructed first contour line;

[0027] S38. Based on the coordinate acquisition logic of the first point on the reconstructed first contour line in steps S34-S37, obtain the coordinates of all points on the reconstructed first contour line of the preset number.

[0028] Furthermore, the contour line correction logic based on the projection of the reconstructed contour line onto the third orthogonal plane and the two-dimensional projection surface of the stone on the third orthogonal plane with its centroid moved to the origin of the coordinate system, performs a compatibility analysis on the third orthogonal plane, and corrects the reconstructed contour line by:

[0029] S41. Based on the distance from the origin of the projection point of the first point of the reconstructed first contour line onto the third orthogonal plane and the length of the two-dimensional projection surface of the stone on the third orthogonal plane with its centroid moved to the origin in the first angular direction, construct the correction coefficient.

[0030] S42. Based on the correction coefficient and the second relational expression, correct the coordinates of the first point of the reconstructed first contour line;

[0031] S43. Based on the coordinate correction logic of steps S41-S42, correct the coordinates of all points of the preset number on the reconstructed first contour line respectively.

[0032] Furthermore, the contour line reconstruction logic based on step S3 and the contour line correction logic based on step S4, respectively reconstructing and correcting several contour lines within each spatial region divided by the first orthogonal plane and the second orthogonal plane, includes:

[0033] S51. Based on the contour line reconstruction logic of steps S34-S38 and the contour line correction logic of steps S41-S43, several contour lines are reconstructed and corrected in the four spatial regions divided by the first orthogonal plane and the second orthogonal plane respectively.

[0034] Furthermore, the construction of a three-dimensional geometric model of the true shape of the stone based on all the contour lines reconstructed and corrected within each spatial region divided by the first orthogonal plane and the second orthogonal plane includes:

[0035] S61. Connect the points at the same position in all the contour lines of the reconstructed and corrected contour lines in each spatial region divided by the first orthogonal plane and the second orthogonal plane in a clockwise or counterclockwise direction to form a three-dimensional solid structure, thereby constructing a three-dimensional geometric model of the true shape of the stone.

[0036] Furthermore, the first relation includes:

[0037] z 1p1 =z 11 +A(z 31 -z 11 )

[0038] Among them, z 1p1 This represents the third-dimensional coordinate of the first point on the reconstructed first contour line, z. 11 This represents the third-dimensional coordinate of the first point on the first two-dimensional projected contour line, z. 31 Let A represent the third-dimensional coordinates of the corresponding first point on the third segment of the two-dimensional projection contour line, and let A represent the first random number.

[0039] Furthermore, the second relation includes:

[0040]

[0041] Where, x 1p1 This represents the first-dimensional coordinate of the first point on the first contour line of the reconstruction, y 1p1 The second-dimensional coordinate, x, represents the first point on the reconstructed first contour line. 11 This represents the first-dimensional coordinate of the first point on the first segment of the two-dimensional projected contour line, y. 31 The second-dimensional coordinates of the corresponding first point on the third two-dimensional projection contour line are represented by B, the second random number is represented by α1, and the first angle is represented by α1.

[0042] The second relation includes: y = x * tanα1.

[0043] Furthermore, the formula for the correction coefficient includes:

[0044]

[0045] in, l1 represents the length of the projection of the reconstructed first contour line onto the third orthogonal plane, l1 represents the length of the two-dimensional projection surface of the stone on the third orthogonal plane with its centroid moved to the origin in the first angular direction, and C represents the correction coefficient.

[0046] The coordinates of the first point of the first contour line in the revised reconstruction specifically include:

[0047] The distance from the origin of the projection point of the first point of the reconstructed first contour line onto the third orthogonal plane is adjusted to...

[0048] In summary, compared with the prior art, the above-described technical solutions conceived by this invention can achieve the following beneficial effects:

[0049] The method of this invention constructs a three-dimensional real model of the boulders based on their actual shape, which improves the simulation accuracy in the finite element simulation of the boulders, thereby improving the accuracy of subsequent finite element analysis of the stability of soil-rock mixed fill slopes from a heterogeneous perspective. Attached Figure Description

[0050] To more clearly illustrate the technical solutions in the embodiments of the present invention, the accompanying drawings used in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0051] Figure 1A flowchart illustrating a three-dimensional modeling method for the true shape of boulders in a soil-rock mixed fill slope, provided in an embodiment of the present invention;

[0052] Figure 2 Two-dimensional projection images and corresponding contour lines scanned on three orthogonal planes are provided for embodiments of the present invention.

[0053] Figure 3 A schematic diagram of a randomly constructed contour line provided in an embodiment of the present invention;

[0054] Figure 4 The image shows the effect of three-dimensional geometric reconstruction of the stone block provided in the embodiment of the present invention. Detailed Implementation

[0055] To make the objectives, technical solutions, and advantages of this invention clearer, the invention will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative and not intended to limit the invention. Furthermore, the technical features involved in the various embodiments of this invention described below can be combined with each other as long as they do not conflict with each other.

[0056] The terms "first," "second," or "third," etc., used in the specification, claims, or accompanying drawings of this application are used to distinguish different objects, not to describe a particular order. Furthermore, the terms "comprising" or "having," and any variations thereof, are intended to cover non-exclusive inclusion. For example, a process, method, system, product, or apparatus that includes a series of steps or units is not limited to the listed steps or units, but may optionally include steps or units not listed, or may optionally include other steps or units inherent to these processes, methods, products, or apparatuses.

[0057] refer to Figure 1 In one embodiment, a three-dimensional modeling method for the true shape of boulders in a soil-rock fill slope may include the following steps.

[0058] Step 1: Collect the two-dimensional projections of the boulders in a real soil-rock fill slope onto three orthogonal planes. Step 1 may specifically include the following sub-steps.

[0059] Step 11: Place the boulders collected from the soil-rock mixed fill slope on site into a transparent acrylic cubic mold, and fix their positions using colorless and transparent epoxy resin glue.

[0060] Step 12: Place the cube mold containing the stone into the AIMS2 image analysis system. By adjusting the position of the top surface of the cube, scan and obtain the two-dimensional projections of the stone on three mutually perpendicular planes. Number the three two-dimensional projections as SP (two-dimensional projection of the stone on the XOY plane), CZ1 (two-dimensional projection of the stone on the YOZ plane), and CZ2 (two-dimensional projection of the stone on the XOZ plane).

[0061] Step 2: Obtain the coordinates of the contour points of each 2D projection using code, calculate the centroid of the contour line of each 2D projection, translate the contour line of each 2D projection so that its centroid moves to the origin, and then recalculate the new coordinates of the contour points of each 2D projection. Step 2 may specifically include the following sub-steps.

[0062] Step 21: Import the 2D projection image numbered SP using the code, and extract the coordinates of the 2D projection outline points of the stone block.

[0063] Step 22: Calculate the centroid coordinates of the outline of the 2D projection numbered SP using code, translate the centroid coordinates of the outline of the 2D projection to the origin, perform translation calculations on the coordinates of the outline points, and output the new coordinate information as a txt document.

[0064] Step 23: Using the operational logic from Steps 21 to 22, extract the contours of the two-dimensional projection planes CZ1 and CZ2 sequentially, such as... Figure 2 As shown. Their centroid coordinates are translated to the origin, the coordinates of the contour points are recalculated, and then the coordinate information is output as a txt document.

[0065] Step 3: Based on the two-dimensional projection of the stone block on XOZ and YOZ, segment the corresponding contour lines to divide the spatial region and randomly reconstruct a contour line. Step 3 may specifically include the following sub-steps.

[0066] Step 31, Reference Figure 3 With the Z-axis ( Figure 3 The contour lines of the two orthogonal planes intersecting in the image are used as the dividing lines. The contour lines of the two-dimensional projection CZ1 are divided into two segments, numbered L1 (contour line in the negative Y-axis direction) and L2 (contour line in the positive Y-axis direction). The contour lines of the two-dimensional projection CZ2 are divided into two segments, numbered L3 (contour line in the negative X-axis direction) and L4 (contour line in the positive X-axis direction). The plane containing these four contour lines divides the XOY plane or the entire space into four regions. Starting from the negative X-axis direction, the regions are numbered Q1, Q2, Q3 and Q4 counterclockwise.

[0067] Step 32: Using the two intersection points of contour line L1 and the Z-axis as the start and end points, use code to randomly divide each contour line L1 into N (N should not be too small, otherwise the contour shape will not be displayed; preferably, N ≥ 20) segments. That is, insert N-1 points between its start point and midpoint, and number the inserted points on contour line L1 sequentially from top to bottom as L. 11 L 12 L 13 ···L 1N L 1N+1 .

[0068] Step 33: Following the segmentation and numbering logic of Step 32, segment and number the remaining three contour lines L2, L3, and L4 respectively, dividing each into N segments as described above, inserting N-1 points in each segment, and numbering them according to the above numbering logic. For example, number the points inserted on contour line L3 from top to bottom as L... 31 L 32 L 33 ···L 3N L 3N+1 .

[0069] Step 34: Output the coordinates of the four contour lines L1-L4 as a txt document.

[0070] Step 35, Reference Figure 3 Reconstruct the first contour line (the contour line within Q1) within the range Q1 (the spatial region between the XOZ plane containing the negative X-axis and the YOZ plane containing the negative Y-axis). Figure 3 The contour line located in the non-orthogonal plane (the contour line of the surface extending between the negative X-axis and negative Y-axis), numbered CG1, has the same top and bottom coordinates as L1 and L3. The plane containing the reconstructed contour line is perpendicular to the XOY plane. An angle α1 (between 0° and α1, and 90°) is randomly created using code. Angle α1 is defined as the angle between the projections of contour line CG1 and contour line L1 onto the XOY plane. Therefore, the angle between the projections of contour line CG1 and contour line L3 onto the XOY plane is 90° - α1.

[0071] Step 36: Select the first point L from top to bottom on the contour lines L1 (contour line in the negative Y-axis direction) and L3 (contour line in the negative X-axis direction). 11 (x 11 ,0,z 11 ) and L 31 (0, y 31 , z 31 Assume the coordinates of the first point P1 of the contour line CG1 are (x... 1p1 y 1p1 , z 1p1 ).

[0072] Use the unifrnd function to randomly generate a number A between 0 and 1, let:

[0073] z 1p1 =z 11 +A(z 31 -z 11 (1)

[0074] From point L 11 and L 31 From the coordinates, we can see that the distances from the projection point of P1 on the XOY plane to the origin are x and x respectively. 11 and y 31 Use the unifrn function to randomly generate a number B between 0 and 1, and let the distance from point P1 to the Z-axis be:

[0075]

[0076] Given the angle between the projection of the newly constructed contour line onto the XOY plane and the projections of contour lines L1 and L3 onto the XOY plane, the equation of the straight line projected onto the XOY plane by the newly constructed contour line can be obtained as follows:

[0077] y=x*tanα1 (3)

[0078] The coordinates of P1 can be obtained by solving equations (1), (2) and (3) simultaneously.

[0079] Step 37: Repeat the operation logic of step 36 to calculate the coordinates of all N-1 points on the reconstructed contour line CG1, and output the coordinate information as a txt document. Note that A and B generated in each round can be different. That is, if A is 0.2 and B is 0.3 in the first round, then A and B generated in the second round may become 0.1 and 0.4. The random numbers in each round can be different.

[0080] Step 4: Based on the horizontal projection, perform a compatibility analysis on the XOY plane and correct the reconstructed contour lines. Step 4 may specifically include the following steps.

[0081] Step 41: The projection of the reconstructed first contour line CG1 onto the XOY plane is the straight line (a straight line L passing through the origin and making an angle α1 with the negative Y-axis). 00 The scanned SP 2D projection contour line and the straight line L 00 The length from the intersection point to the origin of the XOY plane (coordinate origin) is l1; read the coordinates (x, y) of the point where the projection of the contour line CG1 on the XOY plane is farthest from the origin. 1j y 1j ), which is the projection length of the contour line CG1 on the XOY plane.

[0082]

[0083] The projection length of the first point on the reconstructed contour line CG1 onto the XOY plane is adjusted to be [length to be specified]. Similarly, using the above operational logic, adjust the coordinates of all points (N-1 points) on the contour line CG1, and then connect the adjusted new points to obtain a new contour line. (Refer to...) Figure 3 .

[0084] Step 5: Randomly construct several contour lines within each region of Q1-Q4. Step 5 specifically includes the following steps.

[0085] Step 51: Repeat the operation logic from Step 35 to Step 41, reconstruct several contour lines in each region of Q1-Q4, and output the point coordinates of the final contour lines as a txt document using code.

[0086] Step 6: Generate a 3D geometric model with a realistic shape based on the final contour lines. Step 6 specifically includes the following steps.

[0087] Step 61: Import the txt document using code, and sequentially construct all the final contour lines. Connect the points at the same position in different final contour lines in a clockwise or counterclockwise direction (i.e., from the top to the bottom, connect the first point of each reconstructed final contour line to the first point in sequence, then connect all the points of the second point in sequence, and so on, until all points are connected), forming a three-dimensional structure. Output it as an STL file. The realistic reconstruction effect is as follows: Figure 4 As shown.

[0088] The method of this invention constructs a three-dimensional real model of the boulders based on their actual shape, which improves the simulation accuracy in the finite element simulation of the boulders, thereby improving the accuracy of subsequent finite element analysis of the stability of soil-rock mixed fill slopes from a heterogeneous perspective.

[0089] It should be noted that the flowcharts or block diagrams in the accompanying drawings illustrate the architecture, functionality, and operation of possible implementations of systems, methods, or computer program products according to various embodiments of this disclosure. In this regard, each block in the flowchart or block diagram may represent a module, segment, or portion of code, which includes one or more executable instructions for implementing the specified logical function. It should also be noted that in some alternative implementations, the functions indicated in the blocks may occur in a different order than those indicated in the drawings. Furthermore, it should be noted that each block in the block diagram or flowchart, and combinations of blocks in the block diagram or flowchart, can be implemented using a dedicated hardware-based system that performs the specified function or operation, or using a combination of dedicated hardware and computer instructions.

[0090] Those skilled in the art will understand that the features described in the various embodiments and / or claims of this disclosure can be combined and / or combined in various ways, even if such combinations or combinations are not explicitly described in this disclosure. In particular, the features described in the various embodiments and / or claims of this disclosure can be combined and / or combined in various ways without departing from the spirit and teachings of this disclosure, and all such combinations and / or combinations fall within the scope of this disclosure.

[0091] Although this disclosure has been shown and described with reference to specific exemplary embodiments thereof, those skilled in the art will understand that various changes in form and detail may be made to this disclosure without departing from the spirit and scope of the disclosure as defined by the appended claims and their equivalents. Therefore, the scope of this disclosure should not be limited to the above embodiments, but should be defined not only by the appended claims, but also by their equivalents.

Claims

1. A three-dimensional modeling method for the true shape of boulders in a soil-rock mixed fill slope, characterized in that, include: S1. Collect the two-dimensional projections of the boulders in the soil-rock mixed fill slope on three mutually orthogonal planes; S2. Obtain the coordinates of the contour points of each two-dimensional projection and calculate the centroid of the contour line of each two-dimensional projection. Translate the contour line of each two-dimensional projection so that its centroid moves to the origin of the coordinate system. Recalculate the new coordinates of the contour points of each two-dimensional projection. S3. Based on the two-dimensional projection surfaces of the stone on the first and second orthogonal planes, where the centroid has been moved to the origin, the contour lines of the two-dimensional projections on these two orthogonal planes are segmented. Then, in any spatial region divided by these two orthogonal planes, a contour line whose surface is perpendicular to the third orthogonal plane is randomly reconstructed according to the contour line reconstruction logic. Wherein, the first orthogonal plane is the XOZ plane; the second orthogonal plane is the YOZ plane; and the third orthogonal plane is the XOY plane. S4. Based on the projection of the reconstructed contour line onto the third orthogonal plane and the contour line correction logic between the two-dimensional projection plane of the stone block on the third orthogonal plane with its centroid moved to the origin of the coordinates, perform compatibility analysis on the third orthogonal plane and correct the reconstructed contour line. S5. Based on the contour line reconstruction logic of step S3 and the contour line correction logic of step S4, reconstruct and correct several contour lines in each spatial region divided by the first orthogonal plane and the second orthogonal plane respectively. S6. Based on all the contour lines reconstructed and corrected within each spatial region divided by the first orthogonal plane and the second orthogonal plane, construct a three-dimensional geometric model of the true shape of the stone.

2. The three-dimensional modeling method as described in claim 1, characterized in that, The two-dimensional projections of the boulders in the soil-rock fill slope on three mutually orthogonal planes include: S11. Place the boulders collected from the soil-rock mixed fill slope into a cubic mold and fix their position; S12. Scan the three mutually perpendicular mold plates of the cube mold using an image analysis system to obtain a two-dimensional projection of the stone block on the three mutually perpendicular mold plates.

3. The three-dimensional modeling method as described in claim 1, characterized in that, The steps of obtaining the coordinates of the contour points of each two-dimensional projection and calculating the centroid of the contour line of each two-dimensional projection, translating the contour line of each two-dimensional projection so that its centroid moves to the origin of the coordinate system, and recalculating the new coordinates of the contour points of each two-dimensional projection include: S21. Extract the coordinates of the contour points of the two-dimensional projection of the stone block on the first orthogonal plane; S22. Calculate the centroid coordinates of the outline of the two-dimensional projection of the stone block on the first orthogonal plane, translate the outline of the two-dimensional projection so that its centroid moves to the origin of the coordinate system, and perform translation calculations on the coordinates of the outline points of the two-dimensional projection to obtain the new coordinates of the outline points of the two-dimensional projection. S23. Using the methods in steps S21-S22, obtain the new coordinates of the outline points of the two-dimensional projection of the stone on the second orthogonal plane and the third orthogonal plane, respectively.

4. The three-dimensional modeling method as described in claim 1, characterized in that, The step of segmenting the contour lines of the two-dimensional projections of the stone on the first and second orthogonal planes, where the centroids have been moved to the origin, and randomly reconstructing a contour line perpendicular to the third orthogonal plane in any spatial region divided by the two orthogonal planes according to the contour line reconstruction logic, includes: S31. Using the intersection axis of the first orthogonal plane and the second orthogonal plane as the dividing line, the outline of the two-dimensional projection on the two orthogonal planes is divided into four segments; the first orthogonal plane and the second orthogonal plane divide the space into four regions; S32. Taking the two intersection points of each segment of the two-dimensional projection contour line and the above-mentioned intersecting axis as the starting point and the ending point respectively, insert the same preset number of points on each segment of the two-dimensional projection contour line. S33. Randomly generate a first angle that is acute, and determine the first angle as the angle between the first contour line reconstructed in the first region and the projection of the first two-dimensional projection contour line on the third orthogonal plane. S34. Based on the first random number between 0 and 1, construct the first relationship between the coordinates of the first point on the first contour line and the coordinates of the corresponding first points on the two adjacent two-dimensional projected contour lines. S35. Based on a second random number between 0 and 1, construct a second relationship between the coordinates of the first point on the first contour line and the coordinates of the corresponding first points on two adjacent two-dimensional projected contour lines. S36. Based on the first angle, construct a third relational expression for the linear equation of the projection of the first contour line onto the third orthogonal plane; S37. Combine the first to third relations to obtain the coordinates of the first point on the reconstructed first contour line; S38. Based on the coordinate acquisition logic of the first point on the reconstructed first contour line in steps S34-S37, obtain the coordinates of all points on the reconstructed first contour line of the preset number.

5. The three-dimensional modeling method as described in claim 4, characterized in that, The logic for correcting the contour line based on the projection of the reconstructed contour line onto the third orthogonal plane and the contour line between the two-dimensional projection plane of the stone on the third orthogonal plane with its centroid moved to the origin of the coordinate system, performs a compatibility analysis on the third orthogonal plane, and corrects the reconstructed contour line by: S41. Based on the distance from the origin of the projection point of the first point of the reconstructed first contour line onto the third orthogonal plane and the length of the two-dimensional projection surface of the stone on the third orthogonal plane with its centroid moved to the origin in the first angular direction, construct the correction coefficient. S42. Based on the correction coefficient and the second relational expression, correct the coordinates of the first point of the reconstructed first contour line; S43. Based on the coordinate correction logic of steps S41-S42, correct the coordinates of all points of the preset number on the reconstructed first contour line respectively.

6. The three-dimensional modeling method as described in claim 5, characterized in that, The contour line reconstruction logic based on step S3 and the contour line correction logic based on step S4 reconstruct and correct several contour lines in each spatial region divided by the first orthogonal plane and the second orthogonal plane, including: S51. Based on the contour line reconstruction logic of steps S34-S38 and the contour line correction logic of steps S41-S43, several contour lines are reconstructed and corrected in the four spatial regions divided by the first orthogonal plane and the second orthogonal plane respectively.

7. The three-dimensional modeling method as described in claim 6, characterized in that, The process of constructing a three-dimensional geometric model of the true shape of the stone based on all the contour lines reconstructed and corrected within each spatial region divided by the first orthogonal plane and the second orthogonal plane includes: S61. Connect the points at the same position in all the contour lines of the reconstructed and corrected contour lines in each spatial region divided by the first orthogonal plane and the second orthogonal plane in a clockwise or counterclockwise direction to form a three-dimensional solid structure, thereby constructing a three-dimensional geometric model of the true shape of the stone.

8. The three-dimensional modeling method as described in claim 4, characterized in that, The first relation includes: ； in, This represents the third-dimensional coordinate of the first point on the reconstructed first contour line. This represents the third-dimensional coordinate of the first point on the first two-dimensional projected contour line. This represents the third-dimensional coordinates of the corresponding first point on the third segment of the two-dimensional projected contour line. This represents the first random number.

9. The three-dimensional modeling method as described in claim 5, characterized in that, The second relation includes: ； in, This represents the first-dimensional coordinate of the first point on the first contour line of the reconstruction. This represents the second-dimensional coordinate of the first point on the reconstructed first contour line. This represents the first-dimensional coordinate of the first point on the first two-dimensional projected contour line. This represents the second-dimensional coordinates of the corresponding first point on the third segment of the two-dimensional projected contour line. This represents the second random number. Indicates the first angle; The second relation includes: .

10. The three-dimensional modeling method as described in claim 9, characterized in that, The formula for the correction factor includes: ； in, This represents the length of the projection of the reconstructed first contour line onto the third orthogonal plane. The length of the two-dimensional projection of the stone on the third orthogonal plane, with its centroid moved to the origin, in the first angular direction. This represents the correction factor; The coordinates of the first point of the first contour line in the revised reconstruction specifically include: The distance from the origin of the projection point of the first point of the reconstructed first contour line onto the third orthogonal plane is adjusted to... .

Images