CN109472046A - Complicated dam foundation arch dam three-dimensional finite element tetrahedral grid automatic division method - Google Patents
Complicated dam foundation arch dam three-dimensional finite element tetrahedral grid automatic division method Download PDFInfo
- Publication number
- CN109472046A CN109472046A CN201811138818.4A CN201811138818A CN109472046A CN 109472046 A CN109472046 A CN 109472046A CN 201811138818 A CN201811138818 A CN 201811138818A CN 109472046 A CN109472046 A CN 109472046A
- Authority
- CN
- China
- Prior art keywords
- tetrahedron
- cut surface
- dam
- grid
- former
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000005520 cutting process Methods 0.000 claims abstract description 31
- 238000012545 processing Methods 0.000 claims abstract description 4
- 230000001815 facial effect Effects 0.000 claims description 9
- 230000003044 adaptive effect Effects 0.000 abstract description 2
- 238000007781 pre-processing Methods 0.000 abstract description 2
- 238000011160 research Methods 0.000 abstract description 2
- 238000004458 analytical method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000000354 decomposition reaction Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000002203 pretreatment Methods 0.000 description 2
- 238000003325 tomography Methods 0.000 description 2
- 241000682719 Adina Species 0.000 description 1
- 238000004883 computer application Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000012916 structural analysis Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
Abstract
The invention discloses a kind of complicated dam foundation arch dam three-dimensional finite element tetrahedral grid automatic division method, computer establishes the dam body dam foundation tetrahedron background grid under processing condition in the method;By structural planes various in transverse joints in dam body and the dam foundation, launched in background grid as cut surface in three dimensions;Grid is carried out to each tetrahedron background grid and is regenerated according to its cutting position in background grid for each cut surface, until all cut surface dispensings are disposed;Optimize the quality for grid of living again, completes modeling work.The present invention can greatly reduce pre-processing work amount, effectively shorten finite element model settling time, improve research level and efficiency by realizing Adaptive Mesh Generation For Fem Computation.
Description
Technical field
The present invention relates to computer application technologies, and in particular to a kind of complexity dam foundation arch dam Three-D limited elementary tetrahedron
Automatic mesh gridding method.
Background technique
With the continuous development of China's hydraulic and hydroelectric engineering technology, arch dam has a wide range of applications as common dam type, special
It is not to be especially suitable for building induced joint in the high mountain and steep gorge area in China western part.Induced joint more complicated for dam foundation condition, one
As be all made of three-dimensional finite element method carry out structure calculating analysis.The premise of finite element method analysis is integrally-built discretization,
That is finite element mesh, for the induced joint in west area, various structural planes cut rock mass in dam foundation three-dimensional space
And criss-cross form is formed, and there are discontinuity surfaces due tos with transverse joint, attractive joint etc. for arch dams volume mesh itself, therefore
The foundation of its three-dimensional finite element model is a more difficult job.In recent years, with the rapid development of computer hardware and
The extensive use of commercial finite element software, the solution of ultra-large structural analysis problem have become reality, and finite element at this time
Pre-treatment work (i.e. mesh generation) becomes bottleneck of entire analysis process, greatly affected Efficiency, hinders application
Horizontal raising.At this stage, in three dimensions, the algorithm that tetrahedral grid automatically generates is quite mature, and hexahedron
General algorithm, the full hexahedral element grid of complex three-dimensional entity automatically generate the problem beginning not yet in terms of the generation of grid
Fail to obtain solution truly eventually.For arch dam and its dam foundation with complicated boundary condition, full hexahedron net
Lattice automatic division is more difficult.
Summary of the invention
The purpose of the invention is to overcome the shortcomings of above-mentioned background technique, a kind of complicated dam foundation arch dam Three-D limited is provided
Elementary tetrahedron automatic mesh gridding method.
In order to solve the above-mentioned technical problem the technical solution adopted in the present invention the following steps are included:
1) the dam body dam foundation tetrahedron background grid under processing condition is established;
2) it by structural planes various in transverse joints in dam body and the dam foundation, launches in three dimensions as cut surface in background grid
In;
3) it is directed to each cut surface, according to its cutting position in background grid, to each tetrahedron background net
Lattice carry out grid and regenerate, until all cut surface dispensings are disposed;
4) quality for optimizing grid of living again, completes modeling work.
Preferably, various structural planes using the osculating element without thickness carry out mould in transverse joints in dam body and the dam foundation in step 2)
It is quasi-, the binode information that spatial position is overlapped is generated on the cut surface.
Preferably, the spatial position of cut surface determines method in step 2) are as follows:
If cut surface is vertical with the wherein reference axis in three-dimensional system of coordinate, the seat of a point in cut surface is directlyed adopt
Scale value determines the position of cut surface;
If cut surface is not vertical with any reference axis in three-dimensional system of coordinate, take in cutting facial plane not in same straight line
On three points coordinate value determine the position of cut surface or take cutting facial plane in a little with cutting facial plane law vector it is true
Determine cut surface position.
Preferably, the method that grid regenerates is carried out according to cutting position in the step 3) are as follows:
If two vertex of cut surface cutting tetrahedron element and a seamed edge, former tetrahedron are divided into two tetrahedrons
Unit;
If a vertex of cut surface cutting tetrahedron element and two seamed edges, former tetrahedron are divided into a tetrahedron
With a rectangular pyramid, takes a little in rectangular pyramid bottom surface and connect with five vertex of rectangular pyramid, rectangular pyramid can be decomposed into four
Tetrahedron, former tetrahedron are finally broken down into five tetrahedron elements;
If cut surface cuts three seamed edges of tetrahedron element, former tetrahedron is divided into a tetrahedron and a trigone
Column takes a little inside triangular prism and connect with six vertex of triangular prism, triangular prism can be decomposed into two tetrahedrons and three
A rectangular pyramid, each rectangular pyramid can be broken down into four tetrahedrons again, and such a triangular prism can be broken down into 14 four
Face body, former tetrahedron are finally broken down into 15 tetrahedron elements;
If cut surface cuts four seamed edges of tetrahedron element, former tetrahedron is divided into two triangular prisms, each trigone
Column can be broken down into 14 tetrahedrons, and former tetrahedron is finally broken down into 28 tetrahedron elements;
When cut surface cuts the seamed edge of tetrahedron element, the length seamed edge length ratio after cutting is controlled, if length ratio
Greater than 5, then 5 are adjusted to.After the completion of all cuttings, newly-generated tetrahedron element mesh node information replaces divided
Former tetrahedron element mesh node information.
Preferably, the specific steps of step 4) include:
41) volume for calculating newly-generated tetrahedron element, controls the minimum volume of each tetrahedron element, if small in size
In 0.1m3, then it is directly condensed to a bit;
42) the minimum of newly-generated tetrahedron element is calculated to collapse than (tetrahedral height is multiplied by corresponding after 1.24
The ratio of base area), if ratio moves the node location of its vertex correspondence less than 0.1.
Technical solution provided by the invention has the benefit that
The present invention can greatly reduce pre-processing work amount, effectively shorten by realizing Adaptive Mesh Generation For Fem Computation
Meta-model settling time is limited, research level and efficiency are improved.
Detailed description of the invention
Fig. 1 is the structural schematic diagram of difference cutting position when cut surface cuts tetrahedron element in the present invention.
Fig. 2 is the decomposition texture schematic diagram of rectangular pyramid;
Fig. 3 is the decomposition texture schematic diagram of triangular prism;
Fig. 4 is the division of dam body dam foundation background grid;
Fig. 5 is the grid for regenerate automatically after once cutting for a tomography;
Fig. 6 is the grid for regenerate automatically after cutting three times for three tomographies;
Fig. 7 is the method for the present invention flow chart.
Specific embodiment
The present invention is described in further detail with reference to the accompanying drawings and embodiments.
As shown, the present invention proposes a kind of complicated dam foundation arch dam three-dimensional finite element tetrahedral grid automatic division method,
The specific steps of the method include:
1) the dam body dam foundation tetrahedron background grid under processing condition is established;
According to shape design of arch dams and dam foundation orographic condition, by the pre-treatment of the finite element softwares such as ADINA, ABAQUS, ANSYS
Module, establishes the finite element tetrahedral grid of homogeneous dam body and the dam foundation, wherein determining tetrahedral length according to computer performance
Size, the grid is as background grid.
2) it by structural planes various in transverse joints in dam body and the dam foundation, launches in three dimensions as cut surface in background grid
In;
Various structural planes are simulated using the osculating element without thickness in transverse joints in dam body and the dam foundation, raw on the cut surface
The binode information being overlapped at spatial position.
The spatial position of cut surface determines method are as follows:
If cut surface is vertical with the wherein reference axis in three-dimensional system of coordinate, the seat of a point in cut surface is directlyed adopt
Scale value determines the position of cut surface;
If cut surface is not vertical with any reference axis in three-dimensional system of coordinate, take in cutting facial plane not in same straight line
On three points coordinate value determine the position of cut surface or take cutting facial plane in a little with cutting facial plane law vector it is true
Determine cut surface position.
3) it is directed to each cut surface, according to its cutting position in background grid, to each tetrahedron background net
Lattice carry out grid and regenerate, until all cut surface dispensings are disposed;
The method that grid regenerates is carried out according to cutting position are as follows:
If two vertex of cut surface cutting tetrahedron element and a seamed edge, former tetrahedron are divided into two tetrahedrons
Unit;
If a vertex of cut surface cutting tetrahedron element and two seamed edges, former tetrahedron are divided into a tetrahedron
With a rectangular pyramid, takes a little in rectangular pyramid bottom surface and connect with five vertex of rectangular pyramid, rectangular pyramid can be decomposed into four
Tetrahedron, former tetrahedron are finally broken down into five tetrahedron elements;
If cut surface cuts three seamed edges of tetrahedron element, former tetrahedron is divided into a tetrahedron and a trigone
Column takes a little inside triangular prism and connect with six vertex of triangular prism, triangular prism can be decomposed into two tetrahedrons and three
A rectangular pyramid, each rectangular pyramid can be broken down into four tetrahedrons again, and such a triangular prism can be broken down into 14 four
Face body, former tetrahedron are finally broken down into 15 tetrahedron elements;
If cut surface cuts four seamed edges of tetrahedron element, former tetrahedron is divided into two triangular prisms, each trigone
Column can be broken down into 14 tetrahedrons, and former tetrahedron is finally broken down into 28 tetrahedron elements;
When cut surface cuts the seamed edge of tetrahedron element, the length seamed edge length ratio after cutting is controlled, if length ratio
Greater than 5, then 5 are adjusted to.After the completion of all cuttings, newly-generated tetrahedron element mesh node information replaces divided
Former tetrahedron element mesh node information.
4) quality for optimizing grid of living again, completes modeling work;
Specific steps include:
41) volume for calculating newly-generated tetrahedron element, controls the minimum volume of each tetrahedron element, if small in size
In 0.1m3, then it is directly condensed to a bit;
42) the minimum of newly-generated tetrahedron element is calculated to collapse than (tetrahedral height is multiplied by corresponding after 1.24
The ratio of base area), if ratio moves the node location of its vertex correspondence less than 0.1.
Other unspecified parts belong to the prior art.
Claims (5)
1. complicated dam foundation arch dam three-dimensional finite element tetrahedral grid automatic division method, it is characterised in that: the method includes such as
Lower step:
1) the dam body dam foundation tetrahedron background grid under processing condition is established;
2) it by structural planes various in transverse joints in dam body and the dam foundation, is launched in background grid as cut surface in three dimensions;
3) be directed to each cut surface, according to its cutting position in background grid, to each tetrahedron background grid into
Row grid regenerates, until all cut surface dispensings are disposed;
4) quality for optimizing grid of living again, completes modeling work.
2. complexity dam foundation arch dam three-dimensional finite element tetrahedral grid automatic division method according to claim 1, feature
Be: various structural planes are simulated using the osculating element without thickness in transverse joints in dam body and the dam foundation in the step 2), at this
The binode information that spatial position is overlapped is generated on cut surface.
3. complexity dam foundation arch dam three-dimensional finite element tetrahedral grid automatic division method according to claim 1, feature
Be: the spatial position of cut surface determines method in the step 2) are as follows:
If cut surface is vertical with the wherein reference axis in three-dimensional system of coordinate, the coordinate value of a point in cut surface is directlyed adopt
Determine the position of cut surface;
If cut surface is not vertical with any reference axis in three-dimensional system of coordinate, takes in cutting facial plane and be not arranged on the same straight line
The coordinate value of three points determines the position of cut surface or the law vector determination in cutting facial plane a little with cutting facial plane is taken to cut
Face position.
4. complexity dam foundation arch dam three-dimensional finite element tetrahedral grid automatic division method according to claim 1, feature
It is: the method that grid regenerates is carried out according to cutting position in the step 3) are as follows:
If two vertex of cut surface cutting tetrahedron element and a seamed edge, former tetrahedron are divided into two tetrahedron lists
Member;
If a vertex of cut surface cutting tetrahedron element and two seamed edges, former tetrahedron are divided into a tetrahedron and one
A rectangular pyramid takes a little in rectangular pyramid bottom surface and connect with five vertex of rectangular pyramid, rectangular pyramid can be decomposed into four four sides
Body, former tetrahedron are finally broken down into five tetrahedron elements;
If cut surface cuts three seamed edges of tetrahedron element, former tetrahedron is divided into a tetrahedron and a triangular prism,
It takes inside triangular prism a little and is connect with six vertex of triangular prism, triangular prism can be decomposed into two tetrahedrons and three four
Pyramid, each rectangular pyramid can be broken down into four tetrahedrons again, and such a triangular prism can be broken down into 14 tetrahedrons,
Former tetrahedron is finally broken down into 15 tetrahedron elements;
If cut surface cuts four seamed edges of tetrahedron element, former tetrahedron is divided into two triangular prisms, and each triangular prism can
14 tetrahedrons are broken down into, former tetrahedron is finally broken down into 28 tetrahedron elements;
When cut surface cuts the seamed edge of tetrahedron element, the length seamed edge length ratio after cutting is controlled, if length ratio is greater than
5, then it is adjusted to 5.After the completion of all cuttings, newly-generated tetrahedron element mesh node information replaces divided former four
Face body unit mesh node information.
5. complexity dam foundation arch dam three-dimensional finite element tetrahedral grid automatic division method according to claim 1, feature
Be: the specific steps of the step 4) include:
41) volume for calculating newly-generated tetrahedron element, controls the minimum volume of each tetrahedron element, if volume is less than
0.1m3, then it is directly condensed to a bit;
42) the minimum of newly-generated tetrahedron element is calculated to collapse than (tetrahedral height is multiplied by corresponding bottom surface after 1.24
The ratio of area), if ratio moves the node location of its vertex correspondence less than 0.1.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811138818.4A CN109472046A (en) | 2018-09-28 | 2018-09-28 | Complicated dam foundation arch dam three-dimensional finite element tetrahedral grid automatic division method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811138818.4A CN109472046A (en) | 2018-09-28 | 2018-09-28 | Complicated dam foundation arch dam three-dimensional finite element tetrahedral grid automatic division method |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109472046A true CN109472046A (en) | 2019-03-15 |
Family
ID=65664402
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811138818.4A Pending CN109472046A (en) | 2018-09-28 | 2018-09-28 | Complicated dam foundation arch dam three-dimensional finite element tetrahedral grid automatic division method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109472046A (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111159947A (en) * | 2019-12-30 | 2020-05-15 | 中国水利水电科学研究院 | Method and device for reinforcing dam abutment of high arch dam |
CN112560385A (en) * | 2020-12-07 | 2021-03-26 | 芯和半导体科技(上海)有限公司 | Layered sweep grid dividing method applied to packaging |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102622467A (en) * | 2012-02-17 | 2012-08-01 | 上海同岩土木工程科技有限公司 | Method for automatic division of three-dimensional finite element mesh comprising multiple components |
CN108346150A (en) * | 2018-03-06 | 2018-07-31 | 鲁东大学 | A kind of cortex thickness method of estimation based on atlas analysis |
-
2018
- 2018-09-28 CN CN201811138818.4A patent/CN109472046A/en active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102622467A (en) * | 2012-02-17 | 2012-08-01 | 上海同岩土木工程科技有限公司 | Method for automatic division of three-dimensional finite element mesh comprising multiple components |
CN108346150A (en) * | 2018-03-06 | 2018-07-31 | 鲁东大学 | A kind of cortex thickness method of estimation based on atlas analysis |
Non-Patent Citations (5)
Title |
---|
张伯艳 等: "复杂水工结构有限元网格离散的简单实现", 《水利水电技术》 * |
张伯艳 等: "复杂水工结构有限元网格离散的简单实现", 《水利水电技术》, 15 December 1998 (1998-12-15), pages 4 - 6 * |
张渭军;: "基于四面体的地质体可视化与剖分研究", 金属矿山, no. 01, pages 139 - 144 * |
王征 等: "基于投影法的四面体网格切割算法", 《数据采集与处理》 * |
王征 等: "基于投影法的四面体网格切割算法", 《数据采集与处理》, 30 December 2016 (2016-12-30), pages 423 - 427 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111159947A (en) * | 2019-12-30 | 2020-05-15 | 中国水利水电科学研究院 | Method and device for reinforcing dam abutment of high arch dam |
CN111159947B (en) * | 2019-12-30 | 2023-04-07 | 中国水利水电科学研究院 | Method and device for reinforcing dam abutment of high arch dam |
CN112560385A (en) * | 2020-12-07 | 2021-03-26 | 芯和半导体科技(上海)有限公司 | Layered sweep grid dividing method applied to packaging |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103514631B (en) | Magnanimity three-dimensional geological model grid type parallel constructing method based on borehole data | |
CN102306396B (en) | Three-dimensional entity model surface finite element mesh automatic generation method | |
CN103871102B (en) | A kind of road three-dimensional fine modeling method based on elevational point and road profile face | |
CN105787226A (en) | Reconstruction of parameterization model of quadrilateral finite element grid model | |
CN108399290B (en) | Fracture flow MODFLOW simulation method based on fracture network connectivity | |
CN110188423A (en) | A kind of linear engineering structure fast B IM modeling method based on FEM meshing | |
CN103729506B (en) | Complicated model complete hexahedron modeling and geometry remodeling and encryption method | |
CN108763625A (en) | A kind of method and system of FLAC3D softwares pre-treatment | |
CN109472046A (en) | Complicated dam foundation arch dam three-dimensional finite element tetrahedral grid automatic division method | |
JP3988925B2 (en) | Numerical analysis system using mixed grid type solution adaptive grid method | |
CN106844963B (en) | Excavation simulation extremely runs the arch dam three-dimensional grid model automatic division method of overall process | |
CN104200528A (en) | Three-dimensional modeling method based on vector closure | |
CN108389263A (en) | The IGES surface grids rapid generations calculated are solved towards Element BEM | |
CN102903149B (en) | Block forming method and device for geologic model | |
CN110349262A (en) | A kind of tunnel geological forecast three-dimensional modeling method and system based on BIM | |
CN109712241A (en) | A kind of method for building up of the three-dimensional mine physical model comprising goaf | |
CN109493425A (en) | A kind of method for building up of mine three-dimensional goaf physical model | |
CN110990926B (en) | Urban surface building hydrodynamic simulation method based on area correction rate | |
CN103617291B (en) | Equivalent characterization method for reservoir stratum cause unit interface | |
CN106097451B (en) | The parallel Tetrahedron subdivision method of complex underground structure group | |
CN101673416B (en) | Method for dividing hexahedral mesh of piston of internal-combustion engine | |
CN115391899A (en) | Grasshopper-based parametric modeling method for park external packaging steel structure | |
CN112687008A (en) | Method for automatically calculating hydrocarbon-containing area based on geological structure scatter data | |
CN112150628A (en) | Method for building three-dimensional data model in GIS platform by BIM data | |
CN109858076A (en) | A kind of modeling and analysis methods of mine worked-out section |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190315 |