JP2000132713A - Three-dimensional shape generating device and storage medium - Google Patents

Abstract

PROBLEM TO BE SOLVED: To execute display and analysis of a stratum and the position of a fault in a three-dimensional shape by setting the density or the number of points, obtaining the depth of an unknown point on the boundary of a layer so as to generate a triangular mesh and then combining the boundary of this layer so as to define the boundary of plural layers. SOLUTION: The boundary generating means 2 of a processor 1 generates the boundary of the stratum by obtaining the depth of the unknown point, whose depth is unknown, only by an inputted two-dimensional coordinate expressing the boundary of the stratum, based on the inputted two-dimensional coordinate expressing the boundary of the stratum and a known point consisting of its depth. Based on the generated boundary of the stratum, a mesh generating means 3 generates the corresponding unknown point, based on the density of a designated unknown point or the number of them and obtains its depth and generates the triangular mesh by the unknown point whose depth is obtained and the known point. An analyzing means 4 executes the dynamical analysis of the boundary of the stratum, based on this mesh.

Description

DETAILED DESCRIPTION OF THE INVENTION

[0001]

BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a three-dimensional shape forming device for forming a three-dimensional shape of a layer boundary and a recording medium.

[0002]

2. Description of the Related Art Hitherto, the input of data that allows the three-dimensional shape of a stratum to be handled by a personal computer or the like has been performed by manually inputting the three-dimensional coordinates of the stratum one by one.

[0003]

For this reason, it is necessary to input three-dimensional coordinates one by one in order to input a three-dimensional shape of a complicated stratum, and there has been a problem that the operation is enormous.

Further, in the conventional system, only a simple three-dimensional shape such as a cylinder or a sphere can be easily input. In particular, the above-mentioned three-dimensional shape of the stratum cannot be easily input. There has been no alternative but to input the dimensional coordinates one by one.

[0005] The present invention solves these problems,
After defining the layer boundary of the layer with the known points and unknown points of the depth and adding the unknown points of the depth on the layer boundary, set the density or the number of points, obtain the depth of the unknown point, An object of the present invention is to display and analyze a three-dimensional shape, for example, a stratum or a fault position, by creating a mesh and defining a boundary of a plurality of layers by combining the meshes, and displaying the result.

[0006]

Means for solving the problem will be described with reference to FIG. In FIG. 1, a boundary generation means 2 generates a boundary of a layer by obtaining the depth of an unknown point whose depth is unknown based on a known point representing the boundary of the layer.

The mesh generating means 3 generates a triangular mesh based on the boundary between layers and the density of unknown points or the number of unknown points. The analysis means 4 performs a mechanical analysis of the boundary of the layer based on the mesh by a known means.

Next, the operation will be described. An unknown point whose depth is unknown only in the two-dimensional coordinates representing the boundary of the input layer, based on the known points consisting of the two-dimensional coordinates representing the boundary of the layer and the depth thereof, input by the boundary generation means 2 And the mesh boundary is generated by the mesh generation means 3 to determine the corresponding unknown point based on the density of the specified unknown points or the number of unknown points. Then, the depth is obtained, a triangular mesh is generated from the known points and the unknown points from which the depth is obtained, and a dynamic analysis of the boundary of the layer is performed based on the mesh generated by the analysis means 4. I have to.

At this time, a plurality of layers at the boundaries of the generated layers are combined to define a plurality of layers. Further, with respect to an unknown point including its two-dimensional coordinates with respect to the triangle generated by the input known points, the depth corresponding to the unknown point on the surface of the triangle is determined to determine the depth of the unknown point. I have.

Therefore, after defining a layer boundary with known and unknown depth points and adding an unknown depth point on the layer boundary,
Set the density or number of points and find the depth of the unknown point.
By creating a rectangular mesh and defining a plurality of layers by combining the boundaries of the layers, it is possible to display and analyze a three-dimensional shape, for example, a stratum or a fault position, and display the result.

[0011]

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS Next, an embodiment and operation of the present invention will be described in detail with reference to FIGS. Here, a stratum is taken as an example as a layer, and the layers will be described in detail below.

FIG. 1 shows a system configuration diagram of the present invention.
In FIG. 1, a processing device 1 performs various processes in accordance with a program, and here includes a boundary generation unit 2, a mesh generation unit 3, an analysis unit 4, and the like.

The boundary generating means 2 generates the boundary of the stratum by calculating the depth of the unknown point whose depth is unknown based on the known point representing the boundary of the stratum. Based on a known point consisting of two-dimensional coordinates representing the depth and its depth, the depth of an unknown point whose depth is unknown is determined only by the two-dimensional coordinates representing the boundary of the input stratum, and the boundary of the stratum is generated. (Described later with reference to FIGS. 2, 3, and 4).

The mesh generating means 3 generates a triangular mesh based on the boundary of the stratum and the density of unknown points or the number of unknown points, and designates a mesh based on the generated boundary of the stratum. A corresponding unknown point is generated based on the density of the unknown points or the number of unknown points, and its depth is obtained, and a triangular mesh is generated from the known points and the unknown points whose depth has been obtained. .

The analysis means 4 is a means for performing a mechanical analysis of the boundary of the stratum based on the mesh by a known means. The coordinate data 5 stores three-dimensional coordinates of a known point representing a boundary of a stratum and three-dimensional coordinates of two-dimensional coordinates and a depth of an unknown point calculated by the known point.

The analysis mesh data 6 is mesh data for analysis generated by the mesh generation means 3. The result data 7 is a result (for example, a result such as a vibration waveform due to an earthquake) of performing a mechanical analysis by a known means by the analysis means 4 based on the analysis mesh data 6.

The display device 8 displays various data. The input device 9 is for inputting various data and operation instructions, and is an input device such as a mouse and a keyboard.

Next, the operation of the configuration of FIG. 1 will be described in detail according to the order of the flowchart of FIG. FIG. 2 is a flowchart illustrating the operation of the present invention.

In FIG. 2, S1 defines a boundary of a stratum from a known point of depth ● and an unknown point of depth 〇. As shown in FIG. 3A described later on the right side, the coordinates (x, y, z) of a known point of known depth ● are input in addition to the two-dimensional coordinates as the boundary of the stratum. , And the coordinates (x, y) of the unknown point の whose depth is unknown only when the two-dimensional coordinates are known.
In the example of FIG. 3A, 6 known points and 4 unknown points
Enter the point.

In step S2, a point ◎ of unknown depth is added on the boundary line. In this case, if necessary, as shown in FIG. 3B described later on the right side, a necessary unknown point ◎ is added on the boundary line (two-dimensional coordinates are specified, and the unknown depth is specified). Add unknown point ◎).

In step S3, points are automatically generated. This sets the density (point generation interval) or the number of points of two-dimensional data without data in the depth direction. This sets the density or number of unknown points (the depth is unknown at the point of the two-dimensional data) if necessary. For example, as shown in FIG. 3 (d) described later on the right side, according to the set point density or number, three ◎ in the figure (from (b) in FIG. 3 to (c) in FIG. 3) Are added so that they are almost evenly distributed on the boundary surface.

S4 is automatically generated by a known Delaunay triangulation method. In this method, coarse triangular mesh division is performed using certain data in the depth direction and data obtained by automatically calculating the depth direction. Then, the depth of an undefined depth point (unknown point) is calculated from the values of three points of a triangle including the point. This is, as shown in FIG. 3D described later on the right side, a triangle formed by connecting all known points ●, and an unknown point を 持 つ having two-dimensional coordinates included in the triangle is defined by the triangle. Is determined by adding a depth corresponding to the surface (described later with reference to FIG. 4).

S5 is automatically generated by a known Delaunay triangulation method. In this method, all points are divided into triangles to be used as an analysis mesh. As shown in FIG. 3E described later on the right side, as shown in FIG.
The unknown points し た to which the depths obtained from are assigned are connected so as to be as close as possible to a regular triangle in the order of closest points, and a triangular mesh as shown in the figure is automatically generated.

In step S6, a stratum is defined by a combination of a plurality of stratum boundaries. In this method, defining a boundary of one stratum with a triangular mesh in S1 to S5 is repeated for a plurality of strata to generate triangular meshes, respectively, and combining these to define a stratum having a plurality of boundaries.

In step S7, a ground material is set for each stratum.
This sets the ground material of each layer in order to perform a dynamic analysis (earthquake analysis). In step S8, a fault location material and an earthquake parameter are set. This specifies the location of the fault and sets seismic parameters (shear wave velocity, compression wave velocity, etc.).

In step S9, the formation / fault position is confirmed by three-dimensional display. This means that, if necessary, the specified stratum and fault position are three-dimensionally displayed by a known method based on the analysis meshes of multiple layers defined in S6 for the stratum and fault position specified in S8. I do.

In step S10, the BE of the calculation frequency, output position, etc.
Set M analysis control data. This sets (inputs) known BEM analysis control data such as the calculated frequency and output position of the earthquake.

In step S11, an analysis is performed. This is S
Based on the conditions set in 7 to S10, S1 to S6
A publicly known seismic analysis is performed using the triangular mesh data of the boundaries of the plurality of strata generated as described above as an input.

In step S12, the analysis result is displayed and output.
For example, as shown in FIG. 8, which will be described later, a fault in the stratum viewed by animation is displayed as an analysis result. As described above, the user sets the known point と and the unknown point て (only the depth is unknown) on the boundary connecting the known points 必要 as shown in FIG. By setting the density and the number of points, it is possible to automatically generate a triangular mesh. Then, based on the generated triangular meshes of the plurality of layers, the fault position and various parameters required for the earthquake analysis are designated and a seismic analysis is performed by a known method, and as a result, the stratum shown in the animation of FIG. It is possible to display the result shown on the fault.

FIG. 3 shows an explanatory diagram of the present invention. this is,
3D of each grid point of the triangular mesh of the density or the number of points automatically set by inputting the three-dimensional data of the stratum boundary
This is an explanation of the order when automatically generating dimensional coordinates.

FIG. 3A shows a state in which the boundary of the stratum is defined. This is a diagram for explaining the operation of S1 in FIG. 2, where a known point ● represents a point whose two-dimensional coordinates and depth are known,
Unknown point 〇 is a point arbitrarily input by the user as needed on the boundary of the stratum connecting the known points ●.

FIG. 3B shows a state where a point is added on the boundary. This is a diagram for explaining the operation of S2 in FIG. 2, and is a point where two ◎ are added on the boundary. FIG. 3C shows a state where points are automatically created. This is shown in FIG.
FIG. 7 is a diagram for explaining the operation of S3. In this example, three ◎ (◎ indicated as “added in the drawing”) are automatically added in accordance with the setting of the density of points or the number of points.

FIG. 3D shows a state where a coarse mesh is generated. This is a diagram for explaining the operation of S4 in FIG. 2, in which a known mesh is connected to generate a triangular mesh.

FIG. 3E shows a state where an analysis mesh is generated. This is a diagram for explaining the operation of S5 in FIG. 2, and the unknown points included in each triangle of the coarse triangular mesh in FIG. FIG. 4 shows an example of an analysis mesh generated by connecting the known points and the unknown points whose depths have been determined in order of decreasing distance so as to form a regular triangle as much as possible. Things.

As shown in FIGS. 3A to 3E,
This means that an analysis mesh representing the boundary of the stratum could be generated. FIG. 4 shows an explanatory diagram (depth interpolation) of the present invention.

FIG. 4A shows three known points と and one unknown point 含 ま included in the known points 、. FIG. 4B shows an equation for calculating the depth Z of the unknown point 〇. Is shown. Next, a procedure for calculating the depth Z of the unknown point 〇 will be described in detail with reference to FIG.

Step 1: A triangle is formed by taking three known points ● (see FIG. 3 (d)), and the depth values are z1, z2, z3 as shown in the figure. Step 2: The unknown point の of unknown depth included in the known point ● and the other three known points 結 are connected to each other and divided into triangles.

Step 3: Calculate the area of each divided triangle by a
Calculated as 1, a2, a3. Step 4: The depth Z of the unknown point 〇 is obtained by (b) of FIG. As described above, it is possible to obtain the depth Z of the unknown point 〇 whose depth is unknown among the three known points ●.

FIG. 5 is an explanatory diagram (an example of a combination definition of stratum boundaries) of the present invention. This is a collection of three strata. Here, the points where the boundary portions of each layer overlap are collected into one point.

The number is indicated by the stratum boundary number constituting the stratum and the front and back boundaries of the boundary surface. Under the above conditions, stratum number 1,
It is as shown in FIG. Here, the stratum number is the stratum number, the stratum boundary number is the number of stratum boundaries of the stratum number, the stratum boundary number is the stratum boundary number, and the front and back are the front or back of the stratum boundary number boundary. It is. When expressed as shown in FIG.
Three strata are defined as shown in FIG.

FIG. 6 is an explanatory diagram of the present invention. FIG. 6A shows the definition of the stratum based on the curved surface meshes of stratum number 1 (first layer), stratum number 2 (second layer), and stratum number (third layer) described with reference to FIG.

FIG. 6B shows an example of three-dimensional display using an analysis mesh. FIG. 7 shows an explanatory view (stratum) of the present invention. This shows the analysis meshes of formation number 1 (first layer, formation number 2 (second layer), formation number 3 (third layer) in FIG.

FIG. 7A shows an example of an analysis mesh when the first layer is viewed from above. Here, the central portion is in a grid pattern, but actually forms an analysis mesh divided into half triangles, and has a stratum number 1, stratum boundary number 2, stratum boundary number in FIG. It corresponds to.

FIG. 7B shows an example of an analysis mesh when the second layer is viewed from above. Here, the analysis mesh is not shown, but actually constitutes a triangle-divided analysis mesh, and the stratum number 2, stratum boundary number 4, stratum boundary number in FIG. 5 appear. .

FIG. 7C shows an example of an analysis mesh when the third layer is viewed from above. Here, although the analysis mesh is not shown, it actually constitutes a triangle-divided analysis mesh, and the stratum number 3, stratum boundary number 2, and stratum boundary number in FIG. 5 appear.

7 (b) and 7 (c) are the boundaries of the formation which have been neglected for the sake of simplicity. There).

FIG. 8 is an explanatory diagram of the present invention. this is,
The example of the stratum and the fault seen by animation is shown. this is,
It is an example of the analysis result of S12 of FIG. 2, and sets necessary parameters such as the width, height, number of divisions, start position, inclination angle, rake, and rupture propagation speed of a fault, and analyzes and designates an extraction point designated. The results of the above are output as shown in the figure, and when the earthquake propagates from the lower right to the upper left, the state of the seismic acceleration shown is generated as shown in the figure and analyzed sequentially for each extraction point in time. By outputting the result, it can be observed as an animation.

[0048]

As described above, according to the present invention,
After defining the boundary of the layer with the known and unknown points of the depth and adding the unknown point of the depth on the boundary of the layer, the density or number of points is set, and the depth of the unknown point is obtained by calculating the depth of the unknown point. Because a mesh is created and the boundary of this layer is combined to define the boundary of multiple layers, the computer executes the display and analysis of the three-dimensional shape, for example, the stratum and the fault position, and interprets the result. You can now display. As a result, three-dimensional coordinates representing boundaries of a plurality of strata can be input by a simple operation, automatically divided into a triangular mesh for analysis, and used for dynamic analysis.

[Brief description of the drawings]

FIG. 1 is a system configuration diagram of the present invention.

FIG. 2 is a flowchart illustrating the operation of the present invention.

FIG. 3 is an explanatory diagram of the present invention.

FIG. 4 is an explanatory diagram (depth interpolation) of the present invention.

FIG. 5 is an explanatory diagram (an example of a combination definition of stratum boundaries) of the present invention.

FIG. 6 is an explanatory diagram of the present invention.

FIG. 7 is an explanatory view (stratum) of the present invention.

FIG. 8 is an explanatory diagram of the present invention.

[Explanation of symbols]

 1: processing unit 2: boundary generation unit 3: mesh generation unit 4: analysis unit 5: coordinate data 6: analysis mesh data 7: result data 8: display unit 9: input unit

 ────────────────────────────────────────────────── ─── Continuation of the front page (72) Inventor Toshihiro Honjo 2-45 Aomi, Koto-ku, Tokyo F-term in Fujitsu FIP Co., Ltd. (Reference) 5B046 AA00 CA04 FA15 FA18 GA09 HA01 JA04 JA07 5B050 AA01 BA08 BA09 CA04 CA07 EA05 EA24 EA28 FA02 FA09 5B080 AA14 AA18 BA01 CA05 GA02

Claims (4)

[Claims] 1. A three-dimensional shape creating apparatus for creating a three-dimensional shape of a boundary between layers, based on a known point consisting of two-dimensional coordinates representing the input boundary of the layer and its depth. Means for calculating the depth of an unknown point whose depth is unknown only with two-dimensional coordinates representing the boundary of the layer, and generating a layer boundary, based on the generated layer boundary, Means for generating a corresponding unknown point based on the density or the number of unknown points and determining its depth, and generating a triangular mesh from the known point and the unknown point whose depth has been determined; Means for performing a mechanical analysis based on the mesh. 2. A system comprising: means for defining a plurality of layers by combining a plurality of meshes at the boundaries of the generated layers; and means for performing a mechanical analysis based on the generated meshes of the plurality of layers. The three-dimensional shape creating apparatus according to claim 1, wherein: 3. A method according to claim 3, wherein said generated known points are 3
2. The apparatus according to claim 1, further comprising means for determining a depth corresponding to the unknown point including its two-dimensional coordinates on the surface of the triangle and determining the depth as the unknown point. The three-dimensional shape creation device according to claim 2. 4. An unknown point whose depth is unknown only in two-dimensional coordinates representing an input layer boundary, based on known points consisting of two-dimensional coordinates representing the input layer boundary and its depth. Means for calculating the depth of the layer and generating a layer boundary, and generating a corresponding unknown point based on the density of the specified unknown points or the number of unknown points based on the layer boundary generated above. A program that functions as means for determining the depth of the triangle and generating a triangular mesh from the known points and the unknown points for which the depth has been determined, and means for performing a mechanical analysis based on the generated mesh. A computer-readable recording medium on which is recorded.