CN103729506A - 一种复杂模型完全六面体建模及几何体重塑加密方法 - Google Patents
一种复杂模型完全六面体建模及几何体重塑加密方法 Download PDFInfo
- Publication number
- CN103729506A CN103729506A CN201310719119.XA CN201310719119A CN103729506A CN 103729506 A CN103729506 A CN 103729506A CN 201310719119 A CN201310719119 A CN 201310719119A CN 103729506 A CN103729506 A CN 103729506A
- Authority
- CN
- China
- Prior art keywords
- solid
- grid
- hexahedral
- model
- hexahedron
- 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.)
- Granted
Links
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种复杂模型完全六面体建模及几何体重塑加密方法,该复杂模型完全六面体建模及几何体重塑加密方法包括以下步骤:根据模型计算需要建立初步的几何模型,对建立的几何模型按照六面体单元进行划分;简单划分网格后将过渡单元(五面体单元)分解为六面体单元;利用自编的网格转化几何体程序,将所有六面体单元重新按照结点对应几何体顶点方式建立全新的几何体并进行局部控制加密。本发明的一种完全六面体建模及几何体重塑加密方法,能够快速建立完全六面体网格,并能任意控制局部加密密度,是进行大规模科学计算的前提条件,建立的模型在其他数值方法中同样适合。
Description
技术领域
本发明属于有限元建模领域,尤其涉及一种复杂模型完全六面体建模及几何体重塑加密方法。
背景技术
在岩土工程有限元计算分析中,建立一个能反映实际地质、地形条件和岩体结构的计算模型是计算结果可靠性的关键。而在实际建模过程中,往往人为地进行一些简化处理。有时所建立的模型只是所研究问题的一小部分,不能够反映研究目标的整体情况,甚至出现计算模型与地质实际严重不符的情况,从而在很大程度上降低了计算结果的可靠性。因此,为了消除这些非计算导致的误差的因素,一定要建立足够大尺寸模型,把所要研究目标涵盖其中。
单元类型的选择关系到计算效率和计算精度。从几何上来说,任意区域都可以用三角形单元(平面或表面问题)或四面体单元(三维实体问题)来充填逼近,而且几何相似性好。二维三角形单元和三维四面体单元的分析结果精度不是最理想的,而且显示效果也不好。相对说来,六面体单元变形特性好、计算精度高,被人们称作“神圣网格”,六面体单元还可以提供形函数附加的项,来提高解的精度。使用六面体网格相比四面体网格可以大大降低整体单元的数量,使得有限元分析的时间更短。因此,六面体网格具有其他网格无法比拟的优势,在三维计算中得到了广泛的应用。例如在流体力学中,当单元的方向与流动的方向一致时,六面体单元能够更好地捕捉问题的特性。现阶段六面体网格的自动生成始终未能获得真正意义上的解决,构建三维复杂形体要使用全六面体网格是很困难的,需要手动进行。因此,有限元网格划分工作量大而冗繁,目前已经成为有限元求解计算的主要“瓶颈”。
发明内容
本发明的目的在于提供一种复杂模型完全六面体建模及几何体重塑加密方法,旨在快速建立完全六面体网格,并能任意控制局部加密密度。
本发明是这样实现的,一种复杂模型完全六面体建模及几何体重塑加密方法包括根据模型计算需要建立初步的几何模型,对建立的几何模型按照六面体单元进行划分;简单划分网格后将过渡单元(五面体单元)分解为六面体单元;利用自编的“网格转化几何体”程序,将所有六面体单元重新按照结点对应几何体顶点方式建立全新的几何体并进行局部加密。具体步骤如下:
步骤一、根据模型计算需要建立初步的几何模型,对建立的几何模型按照六面体单元进行划分。
如果划分过程中只能划分为四面体单元,必须删掉该区域的几何体。按照三角区域分解技巧和各自建立种子法重新建立新的几何体并简单划分网格。能够形成包含五面体和六面体的基本原则是,几何体要至少保证是五面体和六面体几何体,形状可以不规则。
步骤二、简单划分网格后将五面体分解为六面体。
通过转化模版,应用以下步骤,将含有五面体和六面体单元的模型,完全转化为全六面体单元模型:
(1)选择所有六面体单元,应用相应模版,将原有的六面体单元分解为8倍的六面体单元。
(2)选择所有的五面体单元,应用相应模版,将原有的五面体单元分解为2倍的五面体单元。
(3)选择新生成的所有五面体单元,应用相应模版,将新生成的五面体单元分解为3倍的六面体单元。
步骤三、利用自编的“网格转化几何体”程序,将所有六面体单元重新按照结点对应几何体顶点方式建立全新的几何体并进行局部加密。
进一步,步骤一中所述的三角区域分解技巧,水流疏导法:假设几何体边形成水流流通通道,在模型中间部位放置水源,水能够流通到每个边界,说明该模型设计较合理,此时几何体的边之间是协调的,容易将几何体划分为全六面体单元。在建立几何体模型时,被打断的边尽量依次推向边界。这样在划分网格时,容易实现网格自动协调。
进一步,步骤一中所述的各自建立种子法的步骤是,在相邻边上强行设置包含网格数的分段,然后再进行网格划分,因为划分网格级别限制,有时候这样做也无法实现自动网格协调,就要采取先划分网格再建立种子法。选其中一个几何体单独划分网格,网格会按照既定的分段方式进行划分,然后在相邻边(未划分网格的几何体)建立种子时,选择已经划分好的相邻边的结点作为要划分边的种子,然后删除原来的网格,重新自动划分网格,就能够全部网格自动协调。
效果汇总
本发明的一种复杂模型完全六面体建模及几何体重塑加密方法,能够快速建立完全六面体网格,并能任意控制局部加密密度,是进行大规模科学计算的前提条件,建立的模型在其他数值方法中同样适合。
附图说明
图1是本发明实施例提供的复杂模型完全六面体建模及几何体重塑加密方法的流程图;
图2是本发明实施例提供的模板H5的示意图;
图3是本发明实施例提供的模板W1的示意图;
图4是本发明实施例提供的模板W6的示意图;
图5是本发明实施例提供的六面体网格加密流程图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
图1示出了本发明的复杂模型完全六面体建模及几何体重塑加密方法的流程,如图所示,本发明是这样实现的,一种完全六面体建模及几何体重塑加密方法包括:
根据模型计算需要建立初步的几何模型,对建立的几何模型按照六面体单元进行划分S101;
简单划分网格后将过渡单元(五面体单元)分解为六面体单元S102;
利用自编的“网格转化几何体”程序,将所有六面体单元重新按照结点对应几何体顶点方式建立全新的几何体并进行局部加密S103。
具体步骤如下:
步骤一、根据模型计算需要建立初步的几何模型,对建立的几何模型按照六面体单元进行划分。
如果划分过程中只能划分为四面体单元,必须删掉该区域的几何体。按照三角区域分解技巧和各自建立种子法重新建立新的几何体并简单划分网格。能够形成包含五面体和六面体的基本原则是,几何体要至少保证是五面体和六面体几何体,形状可以不规则。
步骤二、简单划分网格后将五面体分解为六面体。
通过转化模版,应用以下步骤,将含有五面体和六面体单元的模型,完全转化为全六面体单元模型:
(1)选择所有六面体单元,应用图H5模版(图2),将原有的六面体单元分解为8倍的六面体单元。
(2)选择所有的五面体单元,应用W1模版(图3),将原有的五面体单元分解为2倍的五面体单元。
(3)选择新生成的所有五面体单元,应用W6模版(图4),将新生成的五面体单元分解为3倍的六面体单元。
步骤三、利用自编的“网格转化几何体”程序,将所有六面体单元重新按照结点对应几何体顶点方式建立全新的几何体并进行局部加密。
进一步,步骤一中所述的三角区域分解技巧,水流疏导法:假设几何体边形成水流流通通道,在模型中间部位放置水源,水能够流通到每个边界,说明该模型设计较合理,此时几何体的边之间是协调的,容易将几何体划分为全六面体单元。在建立几何体模型时,被打断的边尽量依次推向边界。这样在划分网格时,容易实现网格自动协调。
进一步,步骤一中所述的各自建立种子法的步骤是,在相邻边上强行设置包含网格数的分段,然后再进行网格划分,因为划分网格级别限制,有时候这样做也无法实现自动网格协调,就要采取先划分网格再建立种子法。选其中一个几何体单独划分网格,网格会按照既定的分段方式进行划分,然后在相邻边(未划分网格的几何体)建立种子时,选择已经划分好的相邻边的结点作为要划分边的种子,然后删除原来的网格,重新自动划分网格,就能够全部网格自动协调。
上述虽然结合附图对本发明的具体实施方式进行了描述,但并非对本发明保护范围的限制,所属领域技术人员应该明白,在本发明的技术方案的基础上,本领域技术人员不需要付出创造性的劳动即可做出的各种修改或变形仍在本发明的保护范围之内。
Claims (4)
1.一种复杂模型完全六面体建模及几何体重塑加密方法,其特征在于,该复杂模型完全六面体建模及几何体重塑加密方法包括以下步骤:
根据模型计算需要建立初步的几何模型,对建立的几何模型按照六面体单元进行划分;简单划分网格后将过渡单元或五面体单元分解为六面体单元;利用自编的网格转化几何体程序,将所有六面体单元重新按照结点对应几何体顶点方式建立全新的几何体并进行局部控制加密。
2.如权利要求1所述的复杂模型完全六面体建模及几何体重塑加密方法,其特征在于,具体步骤如下:
步骤一、根据模型计算需要建立初步的几何模型,对建立的几何模型按照六面体单元进行划分;
如果划分过程中只能划分为四面体单元,必须删掉该区域的几何体,按照三角区域分解技巧和各自建立种子法重新建立新的几何体并简单划分网格;能够形成包含五面体和六面体的基本原则是,几何体要至少保证是五面体和六面体几何体,形状可以不规则;
步骤二、简单划分网格后将五面体分解为六面体;
通过转化模版,应用以下步骤,将含有五面体和六面体单元的模型,完全转化为全六面体单元模型:
(1)选择所有六面体单元,应用相应模版,将原有的六面体单元分解为8倍的六面体单元;
(2)选择所有的五面体单元,应用相应模版,将原有的五面体单元分解为2倍的五面体单元;
(3)选择新生成的所有五面体单元,应用相应模版,将新生成的五面体单元分解为3倍的六面体单元;
步骤三、利用自编的网格转化几何体程序,将所有六面体单元重新按照结点对应几何体顶点方式建立全新的几何体并进行局部加密。
3.如权利要求2所述的复杂模型完全六面体建模及几何体重塑加密方法,其特征在于,步骤一中所述的三角区域分解技巧,水流疏导法:假设几何体边形成水流流通通道,在模型中间部位放置水源,水能够流通到每个边界,说明该模型设计较合理,此时几何体的边或面之间是协调的,容易将几何体划分为全六面体单元;在建立几何体模型时,被打断的边尽量依次推向边界。
4.如权利要求2所述的复杂模型完全六面体建模及几何体重塑加密方法,其特征在于,步骤一中所述的各自建立种子法的步骤是,在相邻边上强行设置包含网格数的分段,然后再进行网格划分,因为划分网格级别限制,有时候这样做也无法实现自动网格协调,就要采取先划分网格再建立种子法,选其中一个几何体单独划分网格,网格会按照既定的分段方式进行划分,然后在相邻边未划分网格的几何体建立种子时,选择已经划分好的相邻边的结点作为要划分边的种子,然后删除原来的网格,重新自动划分网格,就能够全部网格自动协调。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310719119.XA CN103729506B (zh) | 2013-12-20 | 2013-12-20 | 一种复杂模型完全六面体建模及几何体重塑加密方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310719119.XA CN103729506B (zh) | 2013-12-20 | 2013-12-20 | 一种复杂模型完全六面体建模及几何体重塑加密方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103729506A true CN103729506A (zh) | 2014-04-16 |
CN103729506B CN103729506B (zh) | 2017-02-15 |
Family
ID=50453580
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310719119.XA Active CN103729506B (zh) | 2013-12-20 | 2013-12-20 | 一种复杂模型完全六面体建模及几何体重塑加密方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103729506B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104361150A (zh) * | 2014-10-23 | 2015-02-18 | 中国科学院合肥物质科学研究院 | 一种用于中子输运的基于MCAM-Geant4自动建模方法的转换简化方法 |
CN105373672A (zh) * | 2015-12-01 | 2016-03-02 | 华中科技大学 | 一种四边形单元实时加密-减密方法 |
CN105760566A (zh) * | 2015-01-06 | 2016-07-13 | 国核(北京)科学技术研究院有限公司 | 一种用于确定热工流体设备阻力的方法及其设备 |
CN107391835A (zh) * | 2017-07-18 | 2017-11-24 | 广东电网有限责任公司电力科学研究院 | 一种高压变压器数值模拟网格划分方法及装置 |
CN110362870A (zh) * | 2019-06-14 | 2019-10-22 | 华南理工大学 | 一种六面体网格划分的复杂花纹轮胎有限元建模方法 |
CN111046614A (zh) * | 2019-12-16 | 2020-04-21 | 西安交通大学 | 一种带绕丝棒束组件结构化网格划分方法 |
CN112712593A (zh) * | 2021-01-20 | 2021-04-27 | 广东电网有限责任公司广州供电局 | 一种基于不规则几何体建模的电力隧道三维设计技术 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080100619A1 (en) * | 2006-10-25 | 2008-05-01 | Coretech System Co., Ltd. | Mesh generation method and computer-readable article for executing the method |
US7671858B1 (en) * | 2005-09-06 | 2010-03-02 | Sandia Corporation | Unconstrained paving and plastering method for generating finite element meshes |
CN102063555A (zh) * | 2011-01-26 | 2011-05-18 | 河海大学 | 基于网格结构的有限元数值模型查错方法 |
US20130231903A1 (en) * | 2012-03-02 | 2013-09-05 | Wan-Chiu LI | Systems and methods for generating a geological model honoring horizons and faults |
-
2013
- 2013-12-20 CN CN201310719119.XA patent/CN103729506B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7671858B1 (en) * | 2005-09-06 | 2010-03-02 | Sandia Corporation | Unconstrained paving and plastering method for generating finite element meshes |
US20080100619A1 (en) * | 2006-10-25 | 2008-05-01 | Coretech System Co., Ltd. | Mesh generation method and computer-readable article for executing the method |
CN102063555A (zh) * | 2011-01-26 | 2011-05-18 | 河海大学 | 基于网格结构的有限元数值模型查错方法 |
US20130231903A1 (en) * | 2012-03-02 | 2013-09-05 | Wan-Chiu LI | Systems and methods for generating a geological model honoring horizons and faults |
Non-Patent Citations (2)
Title |
---|
代星等: "多特征几何体的全六面体网格自动生成方法", 《计算力学学报》 * |
徐能雄等: "基于三维地质建模的复杂构造岩体六面体网格剖分方法", 《岩土工程学报》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104361150A (zh) * | 2014-10-23 | 2015-02-18 | 中国科学院合肥物质科学研究院 | 一种用于中子输运的基于MCAM-Geant4自动建模方法的转换简化方法 |
CN105760566A (zh) * | 2015-01-06 | 2016-07-13 | 国核(北京)科学技术研究院有限公司 | 一种用于确定热工流体设备阻力的方法及其设备 |
CN105760566B (zh) * | 2015-01-06 | 2020-08-04 | 国核(北京)科学技术研究院有限公司 | 一种用于确定热工流体设备阻力的方法及其设备 |
CN105373672A (zh) * | 2015-12-01 | 2016-03-02 | 华中科技大学 | 一种四边形单元实时加密-减密方法 |
CN105373672B (zh) * | 2015-12-01 | 2018-07-06 | 华中科技大学 | 一种四边形单元实时加密-减密方法 |
CN107391835A (zh) * | 2017-07-18 | 2017-11-24 | 广东电网有限责任公司电力科学研究院 | 一种高压变压器数值模拟网格划分方法及装置 |
CN107391835B (zh) * | 2017-07-18 | 2020-09-01 | 广东电科院能源技术有限责任公司 | 一种高压变压器数值模拟网格划分方法及装置 |
CN110362870A (zh) * | 2019-06-14 | 2019-10-22 | 华南理工大学 | 一种六面体网格划分的复杂花纹轮胎有限元建模方法 |
CN110362870B (zh) * | 2019-06-14 | 2023-02-14 | 华南理工大学 | 一种六面体网格划分的复杂花纹轮胎有限元建模方法 |
CN111046614A (zh) * | 2019-12-16 | 2020-04-21 | 西安交通大学 | 一种带绕丝棒束组件结构化网格划分方法 |
CN112712593A (zh) * | 2021-01-20 | 2021-04-27 | 广东电网有限责任公司广州供电局 | 一种基于不规则几何体建模的电力隧道三维设计技术 |
Also Published As
Publication number | Publication date |
---|---|
CN103729506B (zh) | 2017-02-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103729506A (zh) | 一种复杂模型完全六面体建模及几何体重塑加密方法 | |
CN103440364B (zh) | 基于bim模型自动生成wbs节点的方法和系统 | |
CN110069800B (zh) | 具有光滑边界表达的三维结构拓扑优化设计方法及设备 | |
CN103500371A (zh) | 基于bim的变更控制方法 | |
Maadasamy et al. | A hybrid parallel algorithm for computing and tracking level set topology | |
Berger | Cut cells: Meshes and solvers | |
CN106227930A (zh) | 一种基于Matlab的由Midas导入Flac3D的模型识别方法 | |
CN109726433B (zh) | 基于曲面边界条件的三维无粘低速绕流的数值模拟方法 | |
Xu et al. | Building a three dimensional sealed geological model to use in numerical stress analysis software: a case study for a dam site | |
CN104239601A (zh) | 织物片组合的仿真 | |
CN113570275A (zh) | 基于bim和数字高程模型的水资源实时监测系统 | |
CN104200528A (zh) | 基于矢量闭合的三维建模方法 | |
CN110428498A (zh) | 一种将三维地质模型转化为数值计算模型的方法 | |
CN114969860A (zh) | 一种六面体非结构网格自动生成方法 | |
CN114117861A (zh) | 一种基于混合网格的隧道电阻率建模方法及系统 | |
CN106250101A (zh) | 基于MapReduce的叠前偏移并行处理方法和装置 | |
CN103279985B (zh) | 一种复杂地形结构体系三维有限元模型的智能化建模方法 | |
CN110765665B (zh) | 一种地学动态建模方法及系统 | |
CN117077294A (zh) | 基于参数化的结构网格生成方法、装置、设备及介质 | |
CN110704916B (zh) | 基于bim技术的大型复杂三维地质模型网格粗化方法 | |
CN104599318B (zh) | 一种植物三维模型网格无缝融合的方法及系统 | |
CN103236087A (zh) | 一种三棱柱地质模型的构建方法 | |
CN106067192A (zh) | 基于逆插值Loop的网格简化方法 | |
Sanjuan et al. | Applying domain decomposition Schwarz method to accelerate wind field calculation | |
RU2611892C1 (ru) | Способ трехмерного моделирования заданного гидрогеологического объекта, реализуемый в вычислительной системе |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |