CN106844963A - Excavation simulation extremely runs the arch dam three-dimensional grid model automatic division method of overall process - Google Patents

Abstract

The present invention provides a kind of arch dam three-dimensional grid model automatic division method of excavation simulation to overall process is run, including step：Design data according to arch dam determines the x z-plane perspective views of arch axis median plane, and the two-dimensional grid model of arch dam is generated according to the x z-planes perspective view；According to the two-dimensional grid model, the three dimensions surface grids of arch axis median plane are generated；According to the three dimensions surface grids, the three dimensions hyperelement model of arch dam is generated；The three dimensions hyperelement model is finely divided, the Structure Calculation hyperelement model of arch dam is generated；Supplement excavate region three-dimensional space grid, generation can excavation simulation to run overall process arch dam three-dimensional grid model.The three-dimensional grid model of the arch dam that method under this invention is set up, can simultaneously simulate the temperature stress calculation that arch dam excavated, and poured, running overall process, and grid form is regular.Method under this invention, is greatly shortened the mesh generation time, and modeling efficiency greatly improved.

Description

Excavation simulation extremely runs the arch dam three-dimensional grid model automatic division method of overall process

Technical field

The present invention relates to a kind of excavation simulation to the arch dam three-dimensional grid model automatic division method for running overall process, belong to Hydraulic and Hydro-Power Engineering technical field.

Background technology

Finite element adaptive grid generation is the important research field that engineering science intersects with computational science, at present, two-dimensional grid Subdivision method it is more perfect, can be rapidly completed according to design drawing and business software；For three-dimensional grid, the net of rule Lattice can be stretched by two-dimensional grid and obtained, and complicated grid can be obtained by tetrahedral grid through automatic division, such as subdivision The complex space three-dimensional model of space build and hexahedral mesh is used, then need to set up hyperelement block one by one, and Hyperelement block is refined one by one, process is complex.

With continuing to develop for China's hydroelectric project technology, arch dam as a kind of common dam type, have material usage it is few, The series of advantages such as shock resistance is strong.Due to arch dam thinner thickness, structure stress are complex and structure safety of dam body concerning Great, the structural stress of arch dam calculates safely particularly significant.Structural stress is calculated safely includes dam body excavation computation, dam body temperature control The calculating of the aspect such as anticracking calculating and retaining stability analysis, is related to content relatively broad, and existing Structure Calculation is usually root According to demand, multiple grids are set up, multiple grids are respectively calculated, in each process be to influence each other due to various factors , excavation, the influence of temperature stress and retaining to Arch Dam Stress and structure are larger, and often band is calculated respectively to each factor Carry out larger error, and it is more difficult consider the coupled relation of many, therefore set up one can be while considering excavation computation, dam body temperature control The grid model of the multifactor calculating such as anticracking calculating and retaining stability analysis is particularly significant.

The grid model of different calculating demands differs greatly, and such as temperature controlled anticracking generally requires more regular hexahedron net Lattice, and consider the grid model for excavating and, due to reasons such as weak intercalated layer complex structures, generally require more irregular unit.Such as Fruit carries out subdivision using the pattern of hyperelement, then very time-consuming.But consider that Arch Dam is controlled by corresponding shape parameter System, according to the zhou duicheng tuxing of arch dam, you can determine the coordinate position in arch dam upstream and downstream face, zhou duicheng tuxing according to arch dam, Arch Dam parameter and weak intercalated layer morphological parameters, FEM calculation grid is gone out by computer automatic division, then modeling process can Significantly simplified.

The content of the invention

In view of the foregoing, it is an object of the invention to provide a kind of excavation simulation to the arch dam three dimensional network for running overall process Lattice model automatic division method, using the axial section and design data of arch dam, generates the three-dimensional grid model of arch dam, the three-dimensional Grid model can simultaneously simulate the temperature stress calculation that arch dam excavated, and poured, running overall process, when significantly shortening mesh generation Between, simplify arch dam modeling process.

To achieve the above object, the present invention uses following technical scheme：

A kind of excavation simulation extremely runs the arch dam three-dimensional grid model automatic division method of overall process, including step：

S1：Design data according to arch dam determines the x-z-plane perspective view of arch axis median plane, is projected according to the x-z-plane The two-dimensional grid model of figure generation arch dam；

S2：According to the two-dimensional grid model, the three dimensions surface grids of arch axis median plane are set up；

S3：According to the three dimensions surface grids, the three dimensions hyperelement model of arch dam is generated；

S4：The three dimensions hyperelement model is finely divided, the Structure Calculation hyperelement model of arch dam is generated；

S5：Supplement excavate region three-dimensional space grid, generation can excavation simulation to run overall process arch dam three dimensional network Lattice model.

In the step S2,

The three-dimensional surface mesh coordinate for defining arch axis median plane is (x_c, y_c, z_c), it is known that the design data of arch dam, according to arch ring (x, z) coordinate of center line calculates y-coordinate, i.e.,：

Y=f (x, z) (1)

According to the difference of basis and the material number of dam body, the automatic handing-over line for searching basis and dam body：

y_cb=g (x) (2)

If dam body scope along the x-axis direction is [x₀, x₁], for [x₀, x₁] between basis, its y-coordinate is according to lower section Formula value：

Y=g (x) (3)

Such as x ＜ x₀, then：

Y=g (x₀) (4)

Such as x_b＞ x₁, then：

Y=g (x₁) (5)

The nodal information of the three-dimensional surface grids of arch axis median plane is obtained according to formula (1)~(5), the arch axis center is generated The three dimensions surface grids in face.

In the step S3, the node generation method of the three dimensions hyperelement model of the arch dam is：

1) node in dam body upstream and downstream face is determined；

The thickness that (x, z) coordinate according to arch axis median plane obtains arch ring optional position is：

t_c=f (x, z) (6)

According to the z coordinate of arch axis median plane, half central angle alpha of left and right banks is obtained_lAnd α_r, according to the three-dimensional of arch axis median plane (the x of surface grids_c, y_c) coordinate, determine the node coordinate of three dimensions hyperelement model correspondence left bank upstream：

Wherein, λ goes out the angle of cut shown up with horizontal plane for weak intercalated layer, and d is weak intercalated layer region excavation width, λ and d Consider only near weak intercalated layer；

Determine the node coordinate in three dimensions hyperelement model correspondence left bank downstream：

Determine the node coordinate of three dimensions hyperelement model correspondence right bank upstream：

Determine the node coordinate in three dimensions hyperelement model correspondence right bank downstream：

2) the corresponding base node of dam body is determined；

1. it is [x to set dam body along the scope of x-axis₀, x₁], the thickness with basic contact position dam body is：

t_cb=f (x, z_b) (11)

Wherein, z_bBased on and dam body intersection coordinate；

(the x of the three-dimensional surface grids according to arch axis median plane_c, y_c) coordinate, determine three dimensions hyperelement model correspondence left bank The node coordinate of basic upstream：

Determine the node coordinate in three dimensions hyperelement model correspondence left bank basis downstream：

Determine the node coordinate of three dimensions hyperelement model correspondence right bank basis upstream：

Determine the node coordinate in three dimensions hyperelement model correspondence right bank basis downstream：

2. for x ＜ x₀The basis in region, considers in the following manner：

t_cb=f (x₀, z_b) (16)

(the x of the three-dimensional surface grids according to arch axis median plane_c, y_c) coordinate, determine three dimensions hyperelement model correspondence left bank The node coordinate of basic upstream：

Determine the node coordinate in three dimensions hyperelement model correspondence left bank basis downstream：

3. for x ＞ x₁The basis in region, considers in the following manner：

t_cb=f (x₁, z_b) (19)

(the x of the three-dimensional surface grids according to arch axis median plane_c, y_c) coordinate, determine three dimensions hyperelement model correspondence right bank The node coordinate of basic upstream：

Determine the node coordinate in three dimensions hyperelement model correspondence right bank basis downstream：

After determining the node coordinate of three dimensions hyperelement model correspondence upstream and downstream according to formula (6)~(21), institute is generated State three dimensions hyperelement model.

The three dimensions hyperelement of the three dimensions hyperelement model is the node grid of space eight, and eight nodes are respectively The N in downstream₁、N₂、N₃、N₄With the N of upstream₅、N₆、N₇、N₈, i.e. S={ N₁, N₂, N₃, N₄, N₅, N₆, N₇, N₈}。

In the step S4,

Downstream downstream direction extension upwardly on the basis of three dimensions hyperelement model based on arch dam, wherein, upstream The node coordinate that basis extends：

Wherein, θ is the angle of the axis of dam and upstream river direction, and ds is the distance for extending；

The node coordinate that downstream base extends：

Wherein：It is the axis of dam and the angle of upstream river direction；

For the arbitrary node of three dimensions hyperelement model, its corresponding extension node is raw according to formula (24), (25) Into the structural analysis hyperelement model of arch dam；Then, the structural analysis hyperelement model to arch dam is finely divided, generation arch dam Structure Calculation hyperelement model.

In the step S5, generating the method for the arch dam three-dimensional grid model is：Generation is excavated in region and its correspondence The Structure Calculation hyperelement model of body is excavated in the CAD grid arrangement figures of downstream base, generation, by its structure meter with the arch dam Hyperelement model combination is calculated, the arch dam three-dimensional grid model is generated.

The method of the generation Structure Calculation hyperelement model for excavating body is, according to the two-dimensional grid model life for excavating body Into the three dimensions hyperelement model for excavating body, the structure point of body is excavated according to the three dimensions hyperelement model generation for excavating body Analysis hyperelement model, is finely divided to the structural analysis hyperelement model, and the Structure Calculation hyperelement model of body is excavated in generation.

Generating the method for the three dimensions hyperelement model for excavating body is：Structure Calculation according to the arch dam is super single The basic interface node coordinate of meta-model, generates the corresponding node coordinate for excavating body extended model：Wherein,

The node coordinate of left bank upstream：

Wherein, K_x、K_y、K_zIt is respectively three coordinates in direction for excavating body extended model, J_x、J_y、J_zIt is respectively to excavate body Three direction coordinates of the corresponding basic interface node of node on extended model, l is to excavate the node of body extended model and right The distance of interface node is answered, γ is to excavate inclination angle upstream；

The node coordinate of right bank upstream：

The node coordinate in left bank downstream：

Wherein, β is to excavate inclination angle downstream；

The node coordinate in right bank downstream：

It is to the method that the structural analysis hyperelement model is finely divided：To structural analysis hyperelement model each surpass Unit, is carried out increasing node and secondary subdivision grid by river direction to hyperelement side, forms the structural analysis grid after refinement.

It is an advantage of the invention that：

Excavation simulation of the invention extremely runs the arch dam three-dimensional grid model automatic division method of overall process, according to arch dam Design data, generates the three-dimensional grid model of arch dam, based on the three-dimensional grid model, can simultaneously simulate arch dam and excavate, pours, transports The temperature stress calculation of row overall process, grid form is regular, and the mesh generation time is greatly shortened, and modeling effect greatly improved Rate.

Brief description of the drawings

Fig. 1 is the arch dam top view being modeled using the method for the present invention.

Fig. 2 is the arch ring build schematic diagram of arch dam shown in Fig. 1.

Fig. 3 is the arch axis median plane of arch dam shown in Fig. 1 and basement rock, the x-z-plane perspective view of massif of correspondence position.

Fig. 4 is the two-dimensional mesh trrellis diagram based on the generation of x-z-plane perspective view shown in Fig. 3.

Fig. 5 is the three dimensions veil trrellis diagram of arch axis median plane.

Fig. 6 is the three dimensions hyperelement illustraton of model of arch dam.

Fig. 7 is the structural representation of hyperelement.

Fig. 8 is the structural analysis hyperelement illustraton of model of arch dam.

Fig. 9 is the structural representation after hyperelement is once segmented shown in Fig. 7.

Figure 10 is the Structure Calculation hyperelement illustraton of model of arch dam.

Figure 11 is the relation schematic diagram for excavating grid and dam body grid.

Figure 12 is the plane net trrellis diagram for excavating body.

Figure 13 is the three dimensions hyperelement model for excavating body.

Figure 14 is the Structure Calculation hyperelement illustraton of model for excavating body.

Figure 15 is the Yokogawa in the excavation region of a specific embodiment to tangent plane.

Figure 16 is the Yokogawa in the excavation region of a specific embodiment to displacement schematic diagram of calculation result.

Specific embodiment

Below in conjunction with drawings and Examples, the present invention is described in further detail.

Excavation simulation disclosed by the invention extremely runs the arch dam three-dimensional grid model automatic division method of overall process, including step Suddenly：

S1：Design data according to arch dam determines the three-dimensional coordinate system of arch dam, obtains the x-z-plane of arch axis median plane Perspective view, the x-z-plane perspective view according to arch axis median plane generates the two-dimensional grid model of arch dam；

As shown in figure 1, the Arch Dam figure, the arch dam top view that are provided in design data (given data) according to arch dam, The three-dimensional coordinate system of arch dam is set up, the direction of y-axis is consistent with river direction direction, the plane that z-axis and x-axis are formed with y-axis Vertically, based on the three-dimensional coordinate system, the position of arch dam arbitrfary point can use (x, y, z) coordinate representation.

As shown in Fig. 2 including the arch ring build figure of each elevation, the arch ring build according to each elevation in arch dam top view Figure obtains (x, y, z) coordinate of the contact point on each elevation arch center line and basis；As shown in figure 3, being encircleed according to each elevation (x, z) coordinate of the contact point on circle center line and basis obtains the x-z-plane perspective view of arch axis median plane, on this basis, mends Tomography, weak intercalated layer and mountain relief body information (given data) of arch dam region are filled, in can obtain complete arch axis The x-z-plane perspective view of basement rock, the massif of heart face and its relevant position etc..

X-z-plane perspective view generation two-dimensional grid model according to arch axis median plane.First, according to the x- of arch axis median plane Z-plane perspective view, two-dimensional grid layout drawing is drawn using CAD software, and the x-z-plane perspective view of arch axis median plane is drawn with line segment It is divided into some quadrangles and triangle, the end points of line segment must be on the quadrangle and the end points of triangle for being divided；Then, root According to the two-dimensional grid layout drawing of generation using GID, ansys Software Create two-dimensional grid model, specific method is：

The two-dimensional grid arranging graphic of DXF forms is imported in GID softwares；The intersection point of line and line is found, line segment is entered one Step subdivision, each side of quadrangle and triangle for being divided is line segment；By each quadrangle or triangle production quadrangular plan Or triangle projective planum, and save as IGES formatted files；The IGES formatted files of generation import ansys softwares, define dam body material Expect to be the second material, definition basic material is the first material, completes two-dimensional grid subdivision, generates the two-dimensional grid model of arch dam.

S2：The two-dimensional grid model of design data and arch dam according to arch dam, sets up the three dimensions face of arch axis median plane Grid；

As shown in figure 4, the three-dimensional surface mesh coordinate for defining arch axis median plane is (x_c, y_c, z_c), in figure, x-axis is corresponding with z-axis It is the x-axis and z-axis of the three-dimensional surface grids of arch axis median plane, wherein,

According to the design data of arch dam, it is known that the thickness of dam body, then, (x, the z) coordinate according to arch center line can obtain y seats Mark, i.e.,：

Y=f (x, z) (1)

According to basis and the material number of dam body, the automatic handing-over line for searching basis and dam body：

y_cb=g (x) (2)

If dam body scope along the x-axis direction is [x₀, x₁], for [x₀, x₁] between basis, its y-coordinate is according to lower section Formula value：

Y=g (x) (3)

Such as x ＜ x₀, then：

Y=g (x₀) (4)

Such as x_b＞ x₁, then：

Y=g (x₁) (5)

As shown in figure 5, the nodal information of the three-dimensional surface grids of arch axis median plane can be obtained according to formula (1)~(5), by In the nodal information pair of the two-dimensional grid model of the x-z-plane perspective view of the nodal information and arch axis median plane of the three-dimensional surface grids Should be consistent, so, according to formula (1)~(5), you can the three dimensions surface grids of generation arch axis median plane.

S3：The three dimensions surface grids of design data and arch axis median plane according to arch dam, generate the three dimensions of arch dam Hyperelement model；

As shown in fig. 6, the three dimensions surface grids according to arch axis median plane generate the three dimensions hyperelement model of arch dam, The three dimensions hyperelement model is the upper and lower of dam body in a upstream and downstream direction only layer unit, the upstream and downstream face of concrete unit Trip face, the node generation method of three dimensions hyperelement model is：

1) node in dam body upstream and downstream face is determined；

According to the design data of arch dam, (x, z) coordinate of the arch axis median plane according to Fig. 4 can obtain any position of arch ring The thickness put is：

t_c=f (x, z) (6)

Similarly, according to arch dam design data and the z coordinate of arch axis median plane, can obtain half central angle alpha of left and right banks_lWith α_r, (the x of the three-dimensional surface grids according to arch axis median plane_c, y_c) coordinate, determine three dimensions hyperelement model correspondence left bank upstream Node coordinate：

Wherein, λ goes out the angle of cut shown up with horizontal plane for weak intercalated layer, and d is weak intercalated layer region excavation width, λ and d Consider only near weak intercalated layer.

Determine the node coordinate in three dimensions hyperelement model correspondence left bank downstream：

Determine the node coordinate of three dimensions hyperelement model correspondence right bank upstream：

Determine the node coordinate in three dimensions hyperelement model correspondence right bank downstream：

2) the corresponding base node of dam body is determined；

1. it is [x to set dam body along the scope of x-axis₀, x₁], the thickness with basic contact position dam body is：

t_cb=f (x, z_b) (11)

Wherein：z_bBased on and dam body intersection coordinate.

Similarly, (the x of the three-dimensional surface grids according to arch axis median plane_c, y_c) coordinate, determine three dimensions hyperelement model pair Answer the node coordinate of left bank basis upstream：

Determine the node coordinate in three dimensions hyperelement model correspondence left bank basis downstream：

Determine the node coordinate of three dimensions hyperelement model correspondence right bank basis upstream：

Determine the node coordinate in three dimensions hyperelement model correspondence right bank basis downstream：

2. for x ＜ x₀The basis in region, considers in the following manner：

t_cb=f (x₀, z_b) (16)

(the x of the three-dimensional surface grids according to arch axis median plane_c, y_c) coordinate, determine three dimensions hyperelement model correspondence left bank The node coordinate of basic upstream：

Determine the node coordinate in three dimensions hyperelement model correspondence left bank basis downstream：

3. for x ＞ x₁The basis in region, considers in the following manner：

t_cb=f (x₁, z_b) (19)

(the x of the three-dimensional surface grids according to arch axis median plane_c, y_c) coordinate, determine three dimensions hyperelement model correspondence right bank The node coordinate of basic upstream：

Determine the node coordinate in three dimensions hyperelement model correspondence right bank basis downstream：

4. the three dimensions hyperelement model of arch dam is generated.

After determining the node coordinate of three dimensions hyperelement model correspondence upstream and downstream according to formula (6)~(21), you can shape Into three dimensions hyperelement model.Specifically,

For any unit of space surface grids, if it has four nodes, respectively：

S={ S₁, S₂, S₃, S₄} (22)

As shown in fig. 7, the three dimensions hyperelement of generation has eight nodes, the i.e. N in downstream₁、N₂、N₃、N₄With the N of upstream₅、 N₆、N₇、N₈, respectively with S₁、S₂、S₃、S₄It is corresponding, in this case can form the node grid of space eight：

S={ N₁, N₂, N₃, N₄, N₅, N₆, N₇, N₈} (23)

S4：Three dimensions hyperelement model to arch dam is finely divided, and generates the Structure Calculation hyperelement model of arch dam；

Generate arch dam Structure Calculation hyperelement model, i.e., on the basis of the three dimensions hyperelement model based on arch dam under The downstream direction extension further up of trip face, forms the Structure Calculation hyperelement model of arch dam, and specific extended method is：

The node coordinate that upstream base extends：

Wherein, θ is the angle of the axis of dam and upstream river direction, and ds is the distance for extending.

The node coordinate that downstream base extends：

Wherein：It is the axis of dam and the angle of upstream river direction.

As shown in figure 8, according to formula (24), (25), you can obtain the extension node coordinate in upstream and downstream direction.Similarly, it is right In the arbitrary node of three dimensions hyperelement model, its corresponding extension node, using the treatment consistent with formula (24), (25) Method, you can the structural analysis hyperelement model of generation arch dam.

To realize carrying out dam body the calculating analysis of many factors such as stress, temperature, need to be super to the structural analysis of arch dam single Meta-model is finely divided, specially：For each hyperelement of the structural analysis hyperelement model of arch dam, by river direction to super Element sides carry out increasing node and secondary subdivision grid, you can form the structural analysis grid after refinement.Grid after subdivision can Repeatedly segmented with to hyperelement Parallel to the flow direction, as shown in figure 9, as a example by once segmenting, increasing on the basis of hyperelement shown in Fig. 7 Supernumerary segment point N₉, N₁₀, N₁₁And N₁₂, then the hyperelement can be subdivided into two units：

As shown in Figure 10, the Structure Calculation hyperelement model of the arch dam after subdivision is finally given.

S5：Be automatically replenished excavate region three-dimensional space grid, formed can excavation simulation-construction-operation overall process arch Dam three-dimensional grid model.

On the basis of two-dimensional mesh trrellis diagram shown in Fig. 4, the grid supplement that arch dam is excavated into region is complete, and supplement rule is such as Under：

1) the CAD grid arrangement figures on region and its correspondence upstream and downstream basis are excavated in as shown in Figure 11,12, generation, wherein, should The CAD arrangement grid charts for excavating region and its correspondence upstream and downstream basis must be with the respective nodes of two-dimensional mesh trrellis diagram shown in Fig. 4 again Close.The arbitrary node excavated in the CAD arrangement grid charts on region and its correspondence upstream and downstream basis must have corresponding dam body-basis to hand over Interface node, and the Z coordinate of two nodes is required identical.

2) extended model of body is excavated in generation

According to the grid for excavating body horizontal layout, each node to the distance of excavation face is obtained.Because CAD arranges grid chart In, each node has the basic interface node of identical elevation, the basis friendship of the Structure Calculation hyperelement model according to arch dam Interface node coordinate can the corresponding node coordinate for excavating body extended model of generation：

The node coordinate of left bank upstream：

Wherein K_x、K_y、K_zIt is respectively three coordinates in direction for excavating body extended model, J_x、J_y、J_zIt is respectively to excavate body to expand Three direction coordinates of the corresponding basic interface node of node on exhibition model, l is the node and correspondence for excavating body extended model The distance of interface node, γ is to excavate inclination angle upstream.

The node coordinate of right bank upstream：

The node coordinate in left bank downstream：

Wherein, β is to excavate inclination angle downstream.

The node coordinate in right bank downstream：

As shown in Figure 13,14, it is determined that after the node coordinate in excavation body upstream and downstream face, you can form the three dimensions for excavating body Hyperelement model, will excavate the upward downstream aspect of body and extends, and the Structure Calculation hyperelement model of body is excavated in refinement, generation, by it With the Structure Calculation hyperelement model combination of foregoing arch dam, you being formed includes that the overall structure of arch dam, basement rock and massif is calculated Hyperelement model (i.e. arch dam three-dimensional grid model) excavates the Structure Calculation hyperelement model generating method of body and the structure of arch dam Calculate hyperelement model generating method consistent, narration in detail is not repeated herein.

Arch dam three-dimensional grid model based on generation, you can multifactor calculating analysis is carried out to arch dam.It is specific real in one Apply in example, certain big I types power station is located at Jinsha jiang River downstream Ningnan County, Sichuan Province and Yunnan Province Qiaojia County is domestic.Dam is double curvature arch Dam, maximum height of dam 289m.Power plant project is huge, and two sides slope geological complicated condition, dam is closely related with side slope, side slope Stability is of great significance to engineering safety tool.Tomography and weak intercalated layer development are abundant in dam site region.Using this hair Bright method, the modeling of dam body and grid is carried out to the engineering, and its digging process is simulated.Using the method subdivision Dam body grid is the hexahedral mesh of rule, temperature controlled anticracking calculating can be well applied to, using the basis of the method subdivision Grid also can preferably be applied to excavation computation, and, referring to Figure 15,16, modeling speed is fast and grid form is regular, significantly carries for example Modeling efficiency high.

The above is presently preferred embodiments of the present invention and its know-why used, for those skilled in the art For, without departing from the spirit and scope of the present invention, it is any based on technical solution of the present invention on the basis of equivalent change Change, simply replacement etc. obviously changes, belong within the scope of the present invention.

Claims (9)

1. excavation simulation extremely runs the arch dam three-dimensional grid model automatic division method of overall process, it is characterised in that including step：

S1：Design data according to arch dam determines the x-z-plane perspective view of arch axis median plane, is given birth to according to the x-z-plane perspective view Into the two-dimensional grid model of arch dam；

S2：According to the two-dimensional grid model, the three dimensions surface grids of arch axis median plane are set up；

S3：According to the three dimensions surface grids, the three dimensions hyperelement model of arch dam is generated；

S4：The three dimensions hyperelement model is finely divided, the Structure Calculation hyperelement model of arch dam is generated；

S5：Supplement excavate region three-dimensional space grid, generation can excavation simulation to run overall process arch dam three-dimensional grid mould Type.

2. excavation simulation according to claim 1 extremely runs the arch dam three-dimensional grid model automatic division method of overall process, Characterized in that, in the step S2,

The three-dimensional surface mesh coordinate for defining arch axis median plane is (x_c, y_c, z_c), it is known that the design data of arch dam, according to arch center (x, z) coordinate of line calculates y-coordinate, i.e.,：

Y=f (x, z) (1)

According to the difference of basis and the material number of dam body, the automatic handing-over line for searching basis and dam body：

y_cb=g (x) (2)

If dam body scope along the x-axis direction is [x₀, x₁], for [x₀, x₁] between basis, its y-coordinate takes in such a way Value：

Y=g (x) (3)

Such as x ＜ x₀, then：

Y=g (x₀) (4)

Such as x_b＞ x₁, then：

Y=g (x₁) (5)

The nodal information of the three-dimensional surface grids of arch axis median plane is obtained according to formula (1)~(5), the arch axis median plane is generated Three dimensions surface grids.

3. excavation simulation according to claim 2 extremely runs the arch dam three-dimensional grid model automatic division method of overall process, Characterized in that, in the step S3, the node generation method of the three dimensions hyperelement model of the arch dam is：

1) node in dam body upstream and downstream face is determined；

The thickness that (x, z) coordinate according to arch axis median plane obtains arch ring optional position is：

t_c=f (x, z) (6)

According to the z coordinate of arch axis median plane, half central angle alpha of left and right banks is obtained_lAnd α_r, according to the three-dimensional veil of arch axis median plane (the x of lattice_c, y_c) coordinate, determine the node coordinate of three dimensions hyperelement model correspondence left bank upstream：

$\left\{\begin{array}{c}{X}_{lu}=x_c+\frac{t_c}{2}sin\left({\mathrm{\α}}_l\frac{x_c}{x_{dli}}\right)\\ Y_{lu}=y_c-\frac{t_c}{2}\cos \left({\mathrm{\α}}_l\frac{x_c}{x_{dli}}\right)\\ z_{lu}=Z_c+\frac{d}{2}\sin \mathrm{\λ}\end{array}\right.---\left(7\right)$

Wherein, λ goes out the angle of cut shown up with horizontal plane for weak intercalated layer, and d is weak intercalated layer region excavation width, and λ and d only exist Consider near weak intercalated layer；

Determine the node coordinate in three dimensions hyperelement model correspondence left bank downstream：

$\left\{\begin{array}{c}{X}_{ld}=x_c-\frac{t_c}{2}sin\left({\mathrm{\α}}_l\frac{x_c}{x_{dli}}\right)\\ Y_{ld}=y_c+\frac{t_c}{2}\cos \left({\mathrm{\α}}_l\frac{x_c}{x_{dli}}\right)\\ z_{ld}=Z_c-\frac{d}{2}\sin \mathrm{\λ}\end{array}\right.---\left(8\right)$

Determine the node coordinate of three dimensions hyperelement model correspondence right bank upstream：

$\left\{\begin{array}{c}{X}_{ru}=x_c+\frac{t_c}{2}sin\left({\mathrm{\α}}_r\frac{x_c}{x_{dri}}\right)\\ Y_{ru}=y_c-\frac{t_c}{2}\cos \left({\mathrm{\α}}_r\frac{x_c}{x_{dri}}\right)\\ z_{ru}=Z_c+\frac{d}{2}\sin \mathrm{\λ}\end{array}\right.---\left(9\right)$

Determine the node coordinate in three dimensions hyperelement model correspondence right bank downstream：

$\left\{\begin{array}{c}{X}_{rd}=x_c-\frac{t_c}{2}sin\left({\mathrm{\α}}_r\frac{x_c}{x_{dri}}\right)\\ Y_{rd}=y_c+\frac{t_c}{2}\cos \left({\mathrm{\α}}_r\frac{x_c}{x_{dri}}\right)\\ z_{rd}=Z_c-\frac{d}{2}\sin \mathrm{\λ}\end{array}\right.---\left(10\right)$

2) the corresponding base node of dam body is determined；

1. it is [x to set dam body along the scope of x-axis₀, x₁], the thickness with basic contact position dam body is：

t_cb=f (x, z_b) (11)

Wherein, z_bBased on and dam body intersection coordinate；

(the x of the three-dimensional surface grids according to arch axis median plane_c, y_c) coordinate, determine three dimensions hyperelement model correspondence left bank basis The node coordinate of upstream：

$\left\{\begin{array}{c}{X}_{lu}=x_c+\frac{t_{cb}}{2}sin\left({\mathrm{\α}}_l\right)\\ Y_{lu}=y_c-\frac{t_{cb}}{2}\cos \left({\mathrm{\α}}_l\right)\\ z_{lu}=Z_c\end{array}\right.---\left(12\right)$

Determine the node coordinate in three dimensions hyperelement model correspondence left bank basis downstream：

$\left\{\begin{array}{c}{X}_{lu}=x_c-\frac{t_{cb}}{2}sin\left({\mathrm{\α}}_l\right)\\ Y_{lu}=y_c+\frac{t_{cb}}{2}\cos \left({\mathrm{\α}}_l\right)\\ z_{lu}=Z_c\end{array}\right.---\left(13\right)$

Determine the node coordinate of three dimensions hyperelement model correspondence right bank basis upstream：

$\left\{\begin{array}{c}{X}_{ru}=x_c+\frac{t_{cb}}{2}sin\left({\mathrm{\α}}_r\right)\\ Y_{ru}=y_c-\frac{t_{cb}}{2}\cos \left({\mathrm{\α}}_r\right)\\ z_{ru}=Z_c\end{array}\right.---\left(14\right)$

Determine the node coordinate in three dimensions hyperelement model correspondence right bank basis downstream：

$\left\{\begin{array}{c}{X}_{ru}=x_c-\frac{t_{cb}}{2}sin\left({\mathrm{\α}}_r\right)\\ Y_{ru}=y_c+\frac{t_{cb}}{2}\cos \left({\mathrm{\α}}_r\right)\\ z_{ru}=Z_c\end{array}\right.---\left(15\right)$

2. for x ＜ x₀The basis in region, considers in the following manner：

t_cb=f (x₀, z_b) (16)

(the x of the three-dimensional surface grids according to arch axis median plane_c, y_c) coordinate, determine three dimensions hyperelement model correspondence left bank basis The node coordinate of upstream：

$\left\{\begin{array}{c}{X}_{lu}=x_c+\frac{t_{cb}}{2}sin\left(\frac{|x-x_0|}{\left|x_0\right|}{\mathrm{\α}}_l\right)\\ Y_{lu}=y_c-\frac{t_{cb}}{2}\cos \left(\frac{|x-x_0|}{\left|x_0\right|}{\mathrm{\α}}_l\right)\\ z_{lu}=Z_c\end{array}\right.---\left(17\right)$

Determine the node coordinate in three dimensions hyperelement model correspondence left bank basis downstream：

$\left\{\begin{array}{c}{X}_{lu}=x_c-\frac{t_{cb}}{2}\sin \left(\frac{|x-x_0|}{\left|x_0\right|}{\mathrm{\α}}_l\right)\\ Y_{lu}=y_c+\frac{t_{cb}}{2}cos\left(\frac{|x-x_0|}{\left|x_0\right|}{\mathrm{\α}}_l\right)\\ z_{lu}=Z_c\end{array}\right.---\left(18\right)$

3. for x ＞ x₁The basis in region, considers in the following manner：

t_cb=f (x₁, z_b) (19)

(the x of the three-dimensional surface grids according to arch axis median plane_c, y_c) coordinate, determine three dimensions hyperelement model correspondence right bank basis The node coordinate of upstream：

$\left\{\begin{array}{c}{X}_{ru}=x_c+\frac{t_{cb}}{2}\sin \left(\frac{|x-x_0|}{\left|x_0\right|}{\mathrm{\α}}_r\right)\\ Y_{ru}=y_c-\frac{t_{cb}}{2}\cos \left(\frac{|x-x_0|}{\left|x_0\right|}{\mathrm{\α}}_r\right)\\ z_{ru}=Z_c\end{array}\right.---\left(20\right)$

Determine the node coordinate in three dimensions hyperelement model correspondence right bank basis downstream：

$\left\{\begin{array}{c}{X}_{ru}=x_c-\frac{t_{cb}}{2}sin\left(\frac{|x-x_0|}{\left|x_0\right|}{\mathrm{\α}}_r\right)\\ Y_{ru}=y_c+\frac{t_{cb}}{2}\cos \left(\frac{|x-x_0|}{\left|x_0\right|}{\mathrm{\α}}_r\right)\\ z_{ru}=Z_c\end{array}\right.---\left(21\right)$

After determining the node coordinate of three dimensions hyperelement model correspondence upstream and downstream according to formula (6)~(21), generation described three Dimension space hyperelement model.

4. excavation simulation according to claim 3 extremely runs the arch dam three-dimensional grid model automatic division method of overall process, Characterized in that, the three dimensions hyperelement of the three dimensions hyperelement model is the node grid of space eight, eight nodes point Not Wei downstream N₁、N₂、N₃、N₄With the N of upstream₅、N₆、N₇、N₈, i.e. S={ N₁, N₂, N₃, N₄, N₅, N₆, N₇, N₈}。

5. excavation simulation according to claim 4 extremely runs the arch dam three-dimensional grid model automatic division method of overall process, Characterized in that, in the step S4,

Downstream downstream direction extension upwardly on the basis of three dimensions hyperelement model based on arch dam, wherein, upstream base The node coordinate of extension：

$\left\{\begin{array}{c}{X}_{aub}=x_{cu}+dscos\left(\mathrm{\θ}\right)\\ Y_{aub}=y_{cu}+ds\sin \left(\mathrm{\θ}\right)\\ z_{aub}=Z_{cu}\end{array}\right.---\left(24\right)$

Wherein, θ is the angle of the axis of dam and upstream river direction, and ds is the distance for extending；

The node coordinate that downstream base extends：

Wherein：It is the axis of dam and the angle of upstream river direction；

For the arbitrary node of three dimensions hyperelement model, its corresponding extension node, according to formula (24), (25) generation arch The structural analysis hyperelement model on dam；Then, the structural analysis hyperelement model to arch dam is finely divided, and generates the structure of arch dam Calculate hyperelement model.

6. excavation simulation according to claim 5 extremely runs the arch dam three-dimensional grid model automatic division method of overall process, Characterized in that, in the step S5, generating the method for the arch dam three-dimensional grid model is：Region and its correspondence are excavated in generation The CAD grid arrangement figures on upstream and downstream basis, the Structure Calculation hyperelement model of body is excavated in generation, by its structure with the arch dam Hyperelement model combination is calculated, the arch dam three-dimensional grid model is generated.

7. excavation simulation according to claim 6 extremely runs the arch dam three-dimensional grid model automatic division method of overall process, Characterized in that, the method for the generation Structure Calculation hyperelement model for excavating body is, according to the two-dimensional grid mould for excavating body The three dimensions hyperelement model of body is excavated in type generation, and the knot of body is excavated according to the three dimensions hyperelement model generation for excavating body Structure analyzes hyperelement model, and the structural analysis hyperelement model is finely divided, and the Structure Calculation hyperelement mould of body is excavated in generation Type.

8. excavation simulation according to claim 7 extremely runs the arch dam three-dimensional grid model automatic division method of overall process, Characterized in that, the method for the generation three dimensions hyperelement model for excavating body is：According to the Structure Calculation of the arch dam The basic interface node coordinate of hyperelement model, generates the corresponding node coordinate for excavating body extended model：Wherein,

The node coordinate of left bank upstream：

$\left\{\begin{array}{c}{K}_x=J_x+l\\ K_y=J_y+lsin\mathrm{\γ}\\ K_z=J_z\end{array}\right.---\left(27\right)$

Wherein, K_x、K_y、K_zIt is respectively three coordinates in direction for excavating body extended model, J_x、J_yJ_zIt is respectively to excavate body expanded mode Three direction coordinates of the corresponding basic interface node of node in type, l is that the node and correspondence for excavating body extended model join The distance of face node, γ is to excavate inclination angle upstream；

The node coordinate of right bank upstream：

$\left\{\begin{array}{c}{K}_x=J_x-l\\ K_y=J_y+lsin\mathrm{\γ}\\ K_z=J_z\end{array}\right.---\left(28\right)$

The node coordinate in left bank downstream：

$\left\{\begin{array}{c}{K}_x=J_x+l\\ K_y=J_y-lsin\mathrm{\β}\\ K_z=J_z\end{array}\right.---\left(29\right)$

Wherein, β is to excavate inclination angle downstream；

The node coordinate in right bank downstream：

$\left\{\begin{array}{c}{K}_x=J_x-l\\ K_y=J_y-lsin\mathrm{\β}\\ K_z=J_z\end{array}\right.---\left(30\right)$

9. excavation simulation according to claim 8 extremely runs the arch dam three-dimensional grid model automatic division method of overall process, Characterized in that, being to the method that the structural analysis hyperelement model is finely divided：To the every of structural analysis hyperelement model Individual hyperelement, is carried out increasing node and secondary subdivision grid by river direction to hyperelement side, forms the structural analysis after refinement Grid.