CN109215124A - 一种复杂地质条件下大型地下工程3d网格模型的构建方法 - Google Patents
一种复杂地质条件下大型地下工程3d网格模型的构建方法 Download PDFInfo
- Publication number
- CN109215124A CN109215124A CN201811113465.2A CN201811113465A CN109215124A CN 109215124 A CN109215124 A CN 109215124A CN 201811113465 A CN201811113465 A CN 201811113465A CN 109215124 A CN109215124 A CN 109215124A
- Authority
- CN
- China
- Prior art keywords
- grid
- discontinuous face
- grid model
- construction method
- model
- 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
- 238000010276 construction Methods 0.000 title claims abstract description 14
- 238000004458 analytical method Methods 0.000 claims abstract description 16
- 238000009412 basement excavation Methods 0.000 claims abstract description 11
- 230000011218 segmentation Effects 0.000 claims abstract description 5
- 238000002224 dissection Methods 0.000 claims description 5
- 230000002093 peripheral effect Effects 0.000 claims description 3
- 238000000034 method Methods 0.000 abstract description 10
- 238000005516 engineering process Methods 0.000 abstract description 3
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 230000008676 import Effects 0.000 abstract 1
- 230000000694 effects Effects 0.000 description 4
- 101000577210 Homo sapiens Sodium-dependent phosphate transport protein 2A Proteins 0.000 description 2
- 101000604051 Homo sapiens Sodium-dependent phosphate transport protein 3 Proteins 0.000 description 2
- 102100025289 Sodium-dependent phosphate transport protein 1 Human genes 0.000 description 2
- 102100025262 Sodium-dependent phosphate transport protein 2A Human genes 0.000 description 2
- 102100038435 Sodium-dependent phosphate transport protein 3 Human genes 0.000 description 2
- 108091006286 Type III sodium-phosphate co-transporters Proteins 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 208000004350 Strabismus Diseases 0.000 description 1
- 230000001815 facial effect Effects 0.000 description 1
- 239000011229 interlayer Substances 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Computer Graphics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Remote Sensing (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及3D网格模型构建技术领域，具体涉及一种复杂地质条件下大型地下工程3D网格模型的构建方法，包括以下步骤：步骤一:在网格剖分软件ANSYS中根据点线面构建或者导入其它数据格式的形状相对规则的开挖体几何模型；步骤二:在开挖体周围建立一定范围、几何形态简单规则的隔离体作为重点分析区域；步骤三:在重点分析区域外，再划定一个几何形态简单规则的隔离体作为次重点关心区域，网格可适当粗分。本发明的有益效果是：在分层、分步、分段开挖过程中，如果有先前未知的非连续面出现或者已知非连续面产状发生变化时，该方法可以在此原先的3D网格模型上及时更改，而不用再重新建立/更新几何模型后再剖分网格以建立3D网格模型。
Description
技术领域
本发明涉及3D网格模型构建技术领域，具体涉及一种复杂地质条件下大型地下工程3D网格模型的构建方法。
背景技术
当前三维数值分析方法已成为大型地下工程围岩稳定的一种常用方法，但是复杂地质条件下大型地下工程3D网格模型的构建已成为该方法推广应用的一大障碍，存在的具体困难体现在：1)地下结构几何形态复杂，与非连续面的相交往往造成整个三维几何模型被切割成大量碎小、形态复杂的区域，导致网格剖分质量较差；2)大型地下工程开挖，通常采用分期/分层/分步开挖，实际揭露的非连续面往往会与前期勘察所获得非连续面分布情况有所差异，如每一次都按实际情况如实给出非连续面的位置和形状，必须反复调整其位置和几何形状，重新构建几何模型和剖分网格，使得这类大型地下工程3D网格模型的建立工作费时费力。
为此，在现有的复杂地质条件下大型地下洞室3D网格模型方法的基础上，应充分考虑上述问题，提出一种复杂地质条件下大型地下洞室3D网格模型的快速构建方法就十分必要了。
发明内容
本发明的目的在于提供一种复杂地质条件下大型地下工程3D网格模型的构建方法，以解决上述背景技术中提出的问题。
为实现上述目的，本发明提供如下技术方案：
一种复杂地质条件下大型地下工程3D网格模型的构建方法，包括以下步骤：
步骤一:在网格剖分软件ANSYS中根据点线面构建或者导入其它数据格式的形状相对规则的开挖体几何模型；
步骤二:在开挖体周围建立一定范围、几何形态简单规则的隔离体作为重点分析区域；
步骤三:在重点分析区域外，再划定一个几何形态简单规则的隔离体作为次重点关心区域，网格可适当粗分；
步骤四:遵循从“内部到外围、从简单到复杂”的原则对几何模型各组成部分逐步剖分网格，建立不包含非连续面的3D网格模型；
步骤五:根据充填型非连续面的产状以及该非连续面上的一点，确定非连续面所在的平面，根据非连续面的厚度确定一个范围值d，将距该平面法线方向上下2d范围内的网格定义为该非连续面的网格单元集合，从而构建出包含任意产状和数量的非连续面的复杂地质条件下大型地下洞室群3D网格模型。
作为本发明进一步的效果是：在步骤一中，开挖体以六面体精细剖分，使得网格规整，网格单元尺寸小于充填型非连续面厚度。
作为本发明再进一步的效果是：在步骤一中，开挖体周边重点分析区域，以与开挖体网格大小相当的网格进行剖分，网格规整，要求网格单元尺寸小于充填型非连续面厚度。
作为本发明再进一步的效果是：在步骤一中，重点分析区域外围区域网格可以粗糙剖分，减少网格单元总数，减少计算时间。
作为本发明再进一步的效果是：在步骤五中，d值等于非连续面厚度的二分之一。
本发明的有益效果是：
第一、在分层、分步、分段开挖过程中，如果有先前未知的非连续面出现或者已知非连续面产状发生变化时，该方法可以在此原先的3D网格模型上及时更改，而不用再重新建立/更新几何模型后再剖分网格以建立3D网格模型；
第二、该方法还可以随意添加充填型平直非连续面，而不用在原几何模型中添加这些非连续面几何实体后重新剖分网格以建立3D网格模型，在极大地减少建模工作量的同时确保网格质量。
附图说明
图1为本发明复杂地质条件下大型地下工程3D网格模型构建方法的FLAC3D实现过程；
图2为本发明复杂地质条件下大型地下工程3D网格模型简便构建方法的ANSYS实现过程。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明实施例中，一种复杂地质条件下大型地下工程3D网格模型的构建方法，包括以下步骤：
步骤一：在网格剖分软件ANSYS中根据点线面构建或者导入其它数据格式的形状相对规则的开挖体几何模型；
步骤二：在开挖体周围建立一定范围、几何形态简单规则（如六面体）的隔离体作为重点分析区域，以便网格精细剖分；
步骤三：在重点分析区域外，再划定一个几何形态简单规则（如六面体）的隔离体作为次重点关心区域，网格可适当粗分，以便网格从精细剖分到粗分平滑过渡；
步骤四：遵循从“内部到外围、从简单到复杂”的原则对几何模型各组成部分逐步剖分网格，建立不包含非连续面的3D网格模型（如采用FLAC3D进行分析，则需将该3D网格模型由ASNYS导入）；
步骤五：根据充填型非连续面的产状（倾向和倾角）以及该非连续面上的一点，确定非连续面所在的平面（在ANSYS中则可根据3点坐标确定一个平面），根据非连续面的厚度确定一个范围值d，将距该平面法线方向上下2d范围内的网格定义为该非连续面的网格单元集合，从而构建出包含任意产状和数量的非连续面的复杂地质条件下大型地下洞室群3D网格模型。
在步骤一中，开挖体以六面体精细剖分，使得网格规整，网格单元尺寸小于充填型非连续面厚度。
在步骤一中，开挖体周边重点分析区域，以与开挖体网格大小相当的网格进行剖分，网格规整，要求网格单元尺寸小于充填型非连续面厚度。
在步骤一中，重点分析区域外围区域网格可以粗糙剖分，减少网格单元总数，减少计算时间。
在步骤五中，d值等于非连续面厚度的二分之一。
此方法的关键技术难点是在不包含非连续面的3D网格模型中如何圈定出属于非连续面的网格单元集合（步骤五）。
实施例1：
请参阅图1，在该软件中， 实现步骤五的核心命令流为：group f1 ran plane dd 191dip 48 ori 271.769050，734.759597，2264 dist 1.0 group 35，命令流中斜体加粗为命令关键词，其后参数可以根据实际情况修改，该命令的意思为：将组名为35中过点（271.769050,734.759597，2264）、倾向角191°、倾角48°的平面法线方向上下各1.0m距离范围内的网格都归为组f1。按这命令格式，可以在已构建的FLAC3D3D网格模型中继续添加其它产状的非连续面。
实施例2：
请参阅图2，在该软件中，实现步骤五的具体步骤为：
①创建关键点：K，NPT1，X1，Y1，Z1；...K，NPT2，X2，Y2，Z2；NPT3，X3，Y3，Z3（可取非连续面平面内的任意3个点，NPT1-NPT3为点的编号，X1-X3、Y1-Y3、Z1-Z3为点的坐标，可根据实际情况修改）；
②移动工作平面：KWPLAN，WN，NPT1，NPT2，NPT3（将工作平面移动至上述关键点，其Z轴平行于非连续面法线方向）；
③创建局部坐标系：CSWPLA，11，0，1，1，（以上述工作平面坐标原点为原点创建坐标系编号为11的局部坐标系，编号可根据实际情况修改，但数值必须大于等于11）；
④激活当前局部坐标系：CSYS，11，；
⑤选择相应节点：NSEL，S，LOC，Z，-1，1（选择非连续面厚度范围内的相应节点）；
⑥选择相应单元：ESLN,,,（选择依附于上述所选择节点的单元）；
⑦定义分组名：CM，NAME1，ELEM（将被选择单元定义为一个组名为NAME1的分组，注意：该组名不同于flac里的group组名，仅用于网格选择）；
⑧定义材料编号：EMODIF，ALL，50（将非连续面厚度范围内的单元材料编号定义为50，此编号对应于FLAC3D里的group组名）。
按上述命令步骤，可以在已构建的ANSYS3D网格模型中继续添加其它产状的非连续面。
该方法适用于含平直充填型非连续面（如断层、层间/层内错动带等）的大型地下工程围岩稳定分析3D网格模型构建，在此类3D网格模型构建时，如包含非平直充填型非连续面，可以采用多个平直充填型非连续面进行拼接组合的方式进行近似。
对于本领域技术人员而言，显然本发明不限于上述示范性实施例的细节，而且在不背离本发明的精神或基本特征的情况下，能够以其他的具体形式实现本发明。因此，无论从哪一点来看，均应将实施例看作是示范性的，而且是非限制性的，本发明的范围由所附权利要求而不是上述说明限定，因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内。不应将权利要求中的任何附图标记视为限制所涉及的权利要求。
此外，应当理解，虽然本说明书按照实施方式加以描述，但并非每个实施方式仅包含一个独立的技术方案，说明书的这种叙述方式仅仅是为清楚起见，本领域技术人员应当将说明书作为一个整体，各实施例中的技术方案也可以经适当组合，形成本领域技术人员可以理解的其他实施方式。
Claims (5)
1.一种复杂地质条件下大型地下工程3D网格模型的构建方法，其特征在于，包括以下步骤：
步骤一:在网格剖分软件ANSYS中根据点线面构建或者导入其它数据格式的形状相对规则的开挖体几何模型；
步骤二:在开挖体周围建立一定范围、几何形态简单规则的隔离体作为重点分析区域；
步骤三:在重点分析区域外，再划定一个几何形态简单规则的隔离体作为次重点关心区域，网格可适当粗分；
步骤四:遵循从“内部到外围、从简单到复杂”的原则对几何模型各组成部分逐步剖分网格，建立不包含非连续面的3D网格模型；
步骤五:根据充填型非连续面的产状以及该非连续面上的一点，确定非连续面所在的平面，根据非连续面的厚度确定一个范围值d，将距该平面法线方向上下2d范围内的网格定义为该非连续面的网格单元集合，从而构建出包含任意产状和数量的非连续面的复杂地质条件下大型地下洞室群3D网格模型。
2.根据权利要求1所述的复杂地质条件下大型地下工程3D网格模型的构建方法，其特征在于，在步骤一中，开挖体以六面体精细剖分，网格单元尺寸小于充填型非连续面厚度。
3.根据权利要求1或2所述的复杂地质条件下大型地下工程3D网格模型的构建方法，其特征在于，在步骤一中，开挖体周边重点分析区域，以与开挖体网格大小相当的网格进行剖分，网格规整，要求网格单元尺寸小于充填型非连续面厚度。
4.根据权利要求3所述的复杂地质条件下大型地下工程3D网格模型的构建方法，其特征在于，在步骤一中，重点分析区域外围区域网格可以粗糙剖分，减少网格单元总数，减少计算时间。
5.根据权利要求1所述的复杂地质条件下大型地下工程3D网格模型的构建方法，其特征在于，在步骤五中，d值等于非连续面厚度的二分之一。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811113465.2A CN109215124A (zh) | 2018-09-25 | 2018-09-25 | 一种复杂地质条件下大型地下工程3d网格模型的构建方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811113465.2A CN109215124A (zh) | 2018-09-25 | 2018-09-25 | 一种复杂地质条件下大型地下工程3d网格模型的构建方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109215124A true CN109215124A (zh) | 2019-01-15 |
Family
ID=64984636
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811113465.2A Pending CN109215124A (zh) | 2018-09-25 | 2018-09-25 | 一种复杂地质条件下大型地下工程3d网格模型的构建方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109215124A (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113128095A (zh) * | 2021-04-27 | 2021-07-16 | 中国科学院武汉岩土力学研究所 | 一种三维网格模型的构建方法及装置 |
| CN114062150A (zh) * | 2021-11-17 | 2022-02-18 | 山东建筑大学 | 曲线轮廓隧（巷）道围岩隔离体模型承载测试系统及方法 |
| CN114491782A (zh) * | 2022-04-18 | 2022-05-13 | 煤炭科学研究总院有限公司 | 一种基于组合网格剖分的地下工程围岩建模方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101582173A (zh) * | 2009-06-24 | 2009-11-18 | 中国石油集团川庆钻探工程有限公司 | 复杂地质构造块状模型构建方法 |
| CN102855664A (zh) * | 2012-07-23 | 2013-01-02 | 中国科学院研究生院 | 一种复杂块体的三维建模方法 |
| CN105184867A (zh) * | 2015-09-29 | 2015-12-23 | 中国矿业大学 | 基于规则格网和角点网格技术的三维地层建模方法 |
| CN107481320A (zh) * | 2017-08-05 | 2017-12-15 | 中山大学 | 一种三维地质建模的无网格法 |
-
2018
- 2018-09-25 CN CN201811113465.2A patent/CN109215124A/zh active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101582173A (zh) * | 2009-06-24 | 2009-11-18 | 中国石油集团川庆钻探工程有限公司 | 复杂地质构造块状模型构建方法 |
| CN102855664A (zh) * | 2012-07-23 | 2013-01-02 | 中国科学院研究生院 | 一种复杂块体的三维建模方法 |
| CN105184867A (zh) * | 2015-09-29 | 2015-12-23 | 中国矿业大学 | 基于规则格网和角点网格技术的三维地层建模方法 |
| CN107481320A (zh) * | 2017-08-05 | 2017-12-15 | 中山大学 | 一种三维地质建模的无网格法 |
Non-Patent Citations (2)
| Title |
|---|
| 左双英等: "基于FLAC3D的大型地下洞室复杂断层建模及模拟方法研究", 《贵州省岩石力学与工程学会2012年度学术交流论文集》 * |
| 张雨霆等: "基于单元重构的岩土工程复杂地质断层建模方法", 《岩石力学与工程学报》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113128095A (zh) * | 2021-04-27 | 2021-07-16 | 中国科学院武汉岩土力学研究所 | 一种三维网格模型的构建方法及装置 |
| CN113128095B (zh) * | 2021-04-27 | 2022-04-22 | 中国科学院武汉岩土力学研究所 | 一种三维网格模型的构建方法及装置 |
| CN114062150A (zh) * | 2021-11-17 | 2022-02-18 | 山东建筑大学 | 曲线轮廓隧（巷）道围岩隔离体模型承载测试系统及方法 |
| CN114062150B (zh) * | 2021-11-17 | 2023-08-22 | 山东建筑大学 | 曲线轮廓隧（巷）道围岩隔离体模型承载测试系统及方法 |
| CN114491782A (zh) * | 2022-04-18 | 2022-05-13 | 煤炭科学研究总院有限公司 | 一种基于组合网格剖分的地下工程围岩建模方法 |
| CN114491782B (zh) * | 2022-04-18 | 2022-07-05 | 煤炭科学研究总院有限公司 | 一种基于组合网格剖分的地下工程围岩建模方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101582173B (zh) | 复杂地质构造块状模型构建方法 | |
| CN101930623B (zh) | 一种三维道路模型化方法及装置 | |
| CN106529755A (zh) | 一种矿山地质资源储量管理方法 | |
| JP6360192B2 (ja) | 複雑な内部境界を含む領域における非拘束ボロノイ格子の生成 | |
| CN109215124A (zh) | 一种复杂地质条件下大型地下工程3d网格模型的构建方法 | |
| CN102609982B (zh) | 空间地质数据非结构化模式的拓扑发现方法 | |
| Zhang et al. | A new methodology for block identification and its application in a large scale underground cavern complex | |
| CN104318618B (zh) | 一种广义三棱柱空间数据模型的三维剖切方法 | |
| CN109102564B (zh) | 一种复杂地质体数值模型的耦合建模方法 | |
| CN102867332B (zh) | 基于复杂边界约束的多级细分网格曲面拟合方法 | |
| CN106355640A (zh) | 一种地上地下三维一体化空间的处理方法和系统 | |
| CN111383336A (zh) | 三维地质模型构建方法 | |
| CN105225272B (zh) | 一种基于多轮廓线三角网重构的三维实体建模方法 | |
| CN110309613A (zh) | 一种基于bim的隧道开挖台阶法的设计及优化方法 | |
| CN106844963B (zh) | 模拟开挖至运行全过程的拱坝三维网格模型自动剖分方法 | |
| US20240185522A1 (en) | Method for generating a hexahedral mesh | |
| Mlacnik et al. | Sequentially adapted flow-based PEBI grids for reservoir simulation | |
| CN109949415B (zh) | 一种三维地表与地质体模型拓扑一致建模的系统及方法 | |
| CN110349262A (zh) | 一种基于bim的隧道超前地质预报三维建模方法及系统 | |
| CN101303770B (zh) | 一种应用于岩土工程建模的三角网格构建方法 | |
| CN108460838A (zh) | 三维可视化技术与数值模拟技术融合的实现方法与系统 | |
| CN101251930A (zh) | 三维地质模型定制式无损任意分割方法 | |
| Zhang et al. | Method and application of urban 3D rapid modeling of geology based on CAD Borehole Logs | |
| CN111814226A (zh) | 一种隧道工程水文地质bim模型构建方法 | |
| CN113484909B (zh) | 一种基于几何网格化和参数分配的缝洞型储层建立方法 |
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: 20190115 |
|
| RJ01 | Rejection of invention patent application after publication |