CN106650085B - 基于有限元网格的实体模型生成镂空模型的方法 - Google Patents
基于有限元网格的实体模型生成镂空模型的方法 Download PDFInfo
- Publication number
- CN106650085B CN106650085B CN201611181154.0A CN201611181154A CN106650085B CN 106650085 B CN106650085 B CN 106650085B CN 201611181154 A CN201611181154 A CN 201611181154A CN 106650085 B CN106650085 B CN 106650085B
- Authority
- CN
- China
- Prior art keywords
- finite element
- model
- spherical surface
- triangular
- triangular patch
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 22
- 239000007787 solid Substances 0.000 claims description 22
- 239000000463 material Substances 0.000 claims description 4
- 238000005192 partition Methods 0.000 claims description 2
- 239000011796 hollow space material Substances 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 5
- 238000004519 manufacturing process Methods 0.000 description 5
- 239000000654 additive Substances 0.000 description 3
- 230000000996 additive effect Effects 0.000 description 3
- 238000010146 3D printing Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
Images
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]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Image Generation (AREA)
Abstract
本发明提供一种基于有限元网格的实体模型生成镂空模型的方法,在本发明的方法中,将实体模型划分为有限元网格,以有限元网格的各节点为球心生成的球面,然后对各球面进行三角面片网格划分,针对各有限元体单元,对各球面,标记各球面所属的节点与相邻节点之间形成的节点连线在两个球面上的交点,通过标记的两个交点确定对应的两个三角面片网格,将这两个三角面片网格对应的顶点连接以生成三棱柱,连接每个三棱柱侧面的其中一条对角线将其划分为两个三角面片网格,进而完成所有有限元体单元的处理,清除所有的三角面片网格以外的有限元网格区域从而得到实体模型对应的表面和内部均为镂空结构的空间镂空模型,这种生成方法操作简单、可行性高。
Description
技术领域
本发明涉及增材制造领域,尤其涉及一种基于有限元网格的实体模型生成镂空模型的方法。
背景技术
镂空模型能够节省材料,节省成本,并且具有结构通透、样式美观等特点,被广泛运用在工业产品、建筑、艺术等领域。随着增材制造(3D打印技术)的出现,制造具有复杂结构的零部件成为可能,现有的方法能够将实体模型生成表面镂空结构,而内部的未空腔不能生成内部镂空结构,从而无法生成实体模型的表面和内部均为镂空结构的空间镂空模型。随着增材制造(3D打印)技术的进一步发展,实体模型在向空间镂空模型发展,具有减重、节省材料、提高生产效率、降低成本的优点,因此在实体模型上生成空间镂空模型将具有重要意义和实用价值。
发明内容
鉴于现有技术存在的缺陷,本发明的目的在于提供一种基于有限元网格的实体模型生成镂空模型的方法,生成的空间镂空模型的表面和内部均为镂空结构,且操作简单、可行性高。
为了实现上述目的,本发明提供了一种基于有限元网格的实体模型生成镂空模型的方法,包括步骤:S1,建立实体模型;S2,对实体模型进行有限元体单元划分得到有限元网格;S3,以有限元网格的各节点为球心生成球面;S4,对各球面进行三角面片网格划分;S5,针对各有限元体单元,对步骤S4中的各球面,标记各球面所属的节点与相邻节点之间形成的节点连线在两个球面上的交点,通过标记的两个交点确定对应的两个三角面片网格,将这两个三角面片网格对应的顶点连接以生成三棱柱,连接每个三棱柱侧面的其中一条对角线将其划分为两个三角面片网格,重复本步骤,完成所有有限元体单元的处理;S6,清除所有的三角面片网格以外的有限元网格区域,所有的三角面片网格即组成与实体模型对应的空间镂空模型。
本发明的有益效果如下:在根据本发明的基于有限元网格的实体模型生成镂空模型的方法中,将实体模型划分为有限元网格,然后对以有限元网格的各节点为球心生成的球面进行三角面片网格划分,最后对所有有限元体单元进行三角面片网格划分并清除所有的三角面片网格以外的有限元网格区域从而得到实体模型对应的空间镂空模型,生成的空间镂空模型的表面和内部均为镂空结构,且这种生成方法操作简单、可行性高。
附图说明
图1是以正方体实体模型为例进行有限元体单元划分得到有限元网格的示意图;
图2是以图1中的一个有限元体单元的各节点为球心生成球面并对各球面进行三角面片网格划分的示意图;
图3是基于图2对一个有限元体单元进行三角面片网格划分得到三角面片网格的原理示意图;
图4是生成的与正方体实体模型对应的空间镂空模型的线框图;
图5是生成的与正方体实体模型对应的空间镂空模型的渲染图。
具体实施方式
下面参照附图来详细说明根据本发明的基于有限元网格的实体模型生成镂空模型的方法。
参照图1至图5,本发明的基于有限元网格的实体模型生成镂空模型的方法,包括步骤:S1,建立实体模型;S2,对实体模型进行有限元体单元划分得到有限元网格;S3,以有限元网格的各节点为球心生成球面;S4,对各球面进行三角面片网格划分;S5,针对各有限元体单元,对步骤S4中的各球面,标记各球面所属的节点与和该球面相邻的球面所属的节点之间形成的节点连线在两个球面上的交点,通过标记的两个交点确定对应的两个三角面片网格,将这两个三角面片网格对应的顶点连接以生成三棱柱,连接每个三棱柱侧面的其中一条对角线将其划分为两个三角面片网格,重复本步骤,完成所有有限元体单元的处理;S6,清除所有的三角面片网格以外的有限元网格区域,所有的三角面片网格即组成与实体模型对应的空间镂空模型。
在根据本发明的基于有限元网格的实体模型生成镂空模型的方法中,将实体模型划分为有限元网格,然后对以有限元网格的各节点为球心生成的球面进行三角面片网格划分,最后对所有有限元体单元进行三角面片网格划分并清除所有的三角面片网格以外的有限元网格区域从而得到实体模型对应的空间镂空模型,生成的空间镂空模型的表面和内部均为镂空结构,且这种生成方法操作简单、可行性高。
步骤S1中实体模型的建立可用CAD、Soildworks、CATIA等任何三维建模软件。
步骤S2中对实体模型进行有限元体单元划分得到有限元网格可利用现有的CAD平台(如CATIA、PRO/E)或借助专业网格划分软件HyperMesh、ANSYS等。
在步骤S2中对实体模型进行有限元体单元划分得到有限元网格的过程中,根据实体模型确定有限元体单元属性以及网格密度得到有限元网格,有限元体单元属性包括单元类型和分配材料属性。
为了满足结构力学承力性能,平衡受力分布,增加受载时的稳定性,步骤S3中的所生成球面的球径r0满足:
其中,d0为有限元网格中的各节点与和该节点最近的相邻节点之间的距离。
有限元体单元为四面体单元、六面体单元或它们的组合。
最后,结合图1至图4来说明上述方法的操作过程,首先对实体模型进行有限元体单元划分,如图1所示,以正方体实体模型为例,对其进行有限元体单元划分得到有限元网格,其中,有限元体单元为四面体单元。
再以有限元网格的各节点为球心、选取合适的球径r0生成球面,如图2所示,取图1中的有限元网格中的一有限元体单元O-ABC。针对各有限元体单元,对各球面,进行三角面片划分,标记各球面所属的节点与相邻节点之间形成的节点连线在两个球面上的交点,通过标记的两个交点确定对应的两个三角面片网格,将这两个三角面片网格对应的顶点连接以生成三棱柱,连接每个三棱柱侧面的其中一条对角线将其划分为两个三角面片网格,即完成一个有限元体单元的处理。如图3所示,以图2中的球心为O的球面和球心为A的球面为例,球心为O的球面和球心为A的球面经过三角面片网格划分后,球心为O的球面上的三角面片O1O2O3与球心为A的球面上的三角面片A1A2A3被节点连线OA贯穿,故依次连接连线A1O1、A2O2和A3O3,同时连接三角面片的三顶点形成的三角形A1A2A3、O1O2O3,进而与三条连线A1O1、A2O2和A3O3组成三棱柱A1A2A3-O1O2O3。其次,在三棱柱A1A2A3-O1O2O3中,连接每个三棱柱侧面的其中一条对角线(如图3中的O1A2、O2A3、O3A1),即完成了对三棱柱A1A2A3-O1O2O3的三角面片划分。
最后,清除所有的三角面片网格以外的有限元网格区域,如图4所示,球面与三棱柱的所有的三角面片网格即组成与实体模型对应的空间镂空模型,从而完成实体模型生成空间镂空模型的操作,将球面与三棱柱的所有的三角面片网格按照STL格式输出,即得空间镂空模型的STL格式文件。后续可将空间镂空模型的STL格式文件输入3D打印机即可打印出相应的空间镂空模型。
Claims (4)
1.一种基于有限元网格的实体模型生成镂空模型的方法,包括步骤:
S1,建立实体模型;
S2,对实体模型进行有限元体单元划分得到有限元网格;
S3,以有限元网格的各节点为球心生成球面;
S4,对各球面进行三角面片网格划分;
S5,针对各有限元体单元,对步骤S4中的各球面,标记各球面所属的节点与和该球面相邻的球面所属的节点之间形成的节点连线在两个球面上的交点,通过标记的两个交点确定对应的两个三角面片网格,将这两个三角面片网格对应的顶点连接以生成三棱柱,连接每个三棱柱侧面的其中一条对角线将其划分为两个三角面片网格,重复本步骤,完成所有有限元体单元的处理;
S6,清除所有的三角面片网格以外的有限元网格区域,所有的三角面片网格即组成与实体模型对应的空间镂空模型,空间镂空模型的STL格式文件用于输入3D打印机即打印出相应的空间镂空模型。
2.根据权利要求1所述的基于有限元网格的实体模型生成镂空模型的方法,其特征在于,在步骤S2中对实体模型进行有限元体单元划分得到有限元网格的过程中,根据实体模型确定有限元体单元属性以及网格密度得到有限元网格,有限元体单元属性包括单元类型和分配材料属性。
4.根据权利要求1所述的基于有限元网格的实体模型生成镂空模型的方法,其特征在于,有限元体单元为四面体单元、六面体单元或它们的组合。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611181154.0A CN106650085B (zh) | 2016-12-20 | 2016-12-20 | 基于有限元网格的实体模型生成镂空模型的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611181154.0A CN106650085B (zh) | 2016-12-20 | 2016-12-20 | 基于有限元网格的实体模型生成镂空模型的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106650085A CN106650085A (zh) | 2017-05-10 |
CN106650085B true CN106650085B (zh) | 2020-03-10 |
Family
ID=58833450
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611181154.0A Active CN106650085B (zh) | 2016-12-20 | 2016-12-20 | 基于有限元网格的实体模型生成镂空模型的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106650085B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107577861A (zh) * | 2017-08-30 | 2018-01-12 | 清华大学 | 基于有限元方法的实体零件镂空方法 |
CN107617720B (zh) * | 2017-08-30 | 2019-02-22 | 清华大学 | 镂空铸型的3d打印方法 |
CN111085667B (zh) * | 2019-12-30 | 2021-05-14 | 清华大学 | 镂空铸型或镂空砂芯的光滑内腔的设计方法 |
CN111014729B (zh) * | 2019-12-30 | 2021-01-12 | 齐齐哈尔建华机械有限公司 | 一种正十二面体镂空嵌套的加工方法 |
CN111859578B (zh) * | 2020-07-28 | 2023-03-28 | 清华大学 | 自由镂空铸型的设计方法 |
CN112545733B (zh) * | 2020-11-24 | 2023-03-10 | 深圳协同生物医疗科技有限公司 | 手腕固定护具的制作方法、系统及护具 |
CN112861252B (zh) * | 2020-12-24 | 2022-07-12 | 中国航空工业集团公司成都飞机设计研究所 | 用于飞行器结构设计和建模的点阵标准单元及点阵结构 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003281201A (ja) * | 2002-03-25 | 2003-10-03 | Fuji Technical Research Inc | メッシュ生成方法 |
CN101604350A (zh) * | 2009-07-15 | 2009-12-16 | 北京科技大学 | 一种空心型材分流模挤压焊合过程数值模拟技术 |
CN102175474A (zh) * | 2011-01-13 | 2011-09-07 | 天津大学 | 空间网格结构中焊接空心球节点动态承载力的诊断方法 |
CN105183405A (zh) * | 2015-10-12 | 2015-12-23 | 山东大学 | 一种自定义模型表面镂空的3d打印方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2478591A (en) * | 2010-03-12 | 2011-09-14 | Gideon Cube-Sherman | Articulated sub-units |
-
2016
- 2016-12-20 CN CN201611181154.0A patent/CN106650085B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003281201A (ja) * | 2002-03-25 | 2003-10-03 | Fuji Technical Research Inc | メッシュ生成方法 |
CN101604350A (zh) * | 2009-07-15 | 2009-12-16 | 北京科技大学 | 一种空心型材分流模挤压焊合过程数值模拟技术 |
CN102175474A (zh) * | 2011-01-13 | 2011-09-07 | 天津大学 | 空间网格结构中焊接空心球节点动态承载力的诊断方法 |
CN105183405A (zh) * | 2015-10-12 | 2015-12-23 | 山东大学 | 一种自定义模型表面镂空的3d打印方法 |
Non-Patent Citations (1)
Title |
---|
基于体素表达法的STL模型镂空算法研究;陈建树等;《华中科技大学学报(自然科学版)》;20040331;第32卷(第3期);第11-13页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106650085A (zh) | 2017-05-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106650085B (zh) | 基于有限元网格的实体模型生成镂空模型的方法 | |
Chen | 3D texture mapping for rapid manufacturing | |
CN105900145B (zh) | 用于在计算机辅助设计模型中创建三维网格结构以进行增材制造的方法 | |
CN101609564B (zh) | 一种草图式输入的三维网格模型制作方法 | |
WO2018059155A1 (zh) | 带有几何误差的三维实体模型的构建方法及计算机可读存储介质 | |
JP4999522B2 (ja) | 解析メッシュ生成装置 | |
CN111797555B (zh) | 一种基于有限元模型的几何重构方法 | |
JP2005038219A (ja) | 境界表現データからボリュームデータを生成する方法及びそのプログラム | |
CN105313336A (zh) | 一种薄壳体3d打印优化方法 | |
CN109002668B (zh) | 一种连续体与离散体耦合拓扑优化方法 | |
CN107578472A (zh) | 三维表面三角网络模型的角度优化方法及装置 | |
JP2022535688A (ja) | 付加製造された物品の軽量化及び/又は設計のための方法 | |
CN106682299A (zh) | 一种用于激光选区烧结砂型分区域变强度的设计制造方法 | |
CN111859578B (zh) | 自由镂空铸型的设计方法 | |
CN113887095A (zh) | 一种基于等几何分析的渐进式结构拓扑优化方法 | |
Dawes et al. | Viscous layer meshes from level sets on cartesian meshes | |
CN116451540A (zh) | 一种建筑表面施工曲面分格方法 | |
CN115408796A (zh) | 一种多孔结构建模方法、装置、设备及应用 | |
Xing et al. | Shell thickening for extrusion-based ceramics printing | |
Zhu et al. | A simple automatic hexahedron mesh generation and polyhedral smoothed finite element method for mechanics problems | |
CN114693887A (zh) | 一种复杂点阵结构体参数化造型方法 | |
CN106097415A (zh) | 一种2.5维绳索扭曲动画生成方法 | |
CN115674684B (zh) | 内部自支撑结构镂空的3d打印方法及桌面摆件打印方法 | |
CN104851125A (zh) | 一种植物叶子三维模型建模方法及系统 | |
CN117292085B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |