CN115984511A - 一种基于cad的平行六面体体积平均共形网格剖分方法 - Google Patents
一种基于cad的平行六面体体积平均共形网格剖分方法 Download PDFInfo
- Publication number
- CN115984511A CN115984511A CN202310257023.XA CN202310257023A CN115984511A CN 115984511 A CN115984511 A CN 115984511A CN 202310257023 A CN202310257023 A CN 202310257023A CN 115984511 A CN115984511 A CN 115984511A
- Authority
- CN
- China
- Prior art keywords
- volume
- parallelepiped
- grid
- triangle
- polygon
- 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
- 238000000034 method Methods 0.000 title claims abstract description 47
- 238000011960 computer-aided design Methods 0.000 title description 12
- 239000000463 material Substances 0.000 claims abstract description 30
- 238000004364 calculation method Methods 0.000 claims abstract description 23
- 238000004088 simulation Methods 0.000 claims abstract description 5
- 230000001788 irregular Effects 0.000 claims description 38
- 238000004422 calculation algorithm Methods 0.000 claims description 27
- 238000000926 separation method Methods 0.000 claims description 22
- 238000005520 cutting process Methods 0.000 claims description 19
- 230000011218 segmentation Effects 0.000 claims description 16
- 238000012545 processing Methods 0.000 abstract description 2
- 238000005070 sampling Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 4
- 230000007547 defect Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 230000005672 electromagnetic field Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
Images
Landscapes
- Image Generation (AREA)
Abstract
本发明属于电数字数据处理技术领域，涉及平行六面体体积平均共形网格剖分方法。一种基于CAD的平行六面体体积平均共形网格剖分方法，包括：用CAD对物体进行建模，得到三角网模型；对仿真空间进行平行六面体网格划分；读取三角网模型和平行六面体网格，计算出该模型的阶梯网格；遍历所有三角形，计算每个三角形的包围盒；遍历包围盒内的体像素，检测包围盒内的体像素与三角形的碰撞情况，得到体像素集合和每一个体像素对应的与其碰撞的三角形集合；遍历体像素集合，计算不同材料在体像素内的填充体积和体积平均等效材料系数，更新阶梯网格中的材料系数，得到平行六面体体积平均共形网格。本发明的方法计算速度快、精度高并且计算复杂度低。
Description
技术领域
本发明属于电数字数据处理技术领域，具体涉及平行六面体体积平均共形网格剖分方法。
背景技术
网格剖分的优化是数值计算中十分重要的一步，有限差分法用常用的网格为平行六面体网格（体像素网格），例如电磁场时域有限差分法(FDTD)。而平行六面体网格一个格子只能为一种材料，故在对弯曲表面建模时会呈现阶梯形状，势必会降低网格的可靠性，因此许多共形网格技术被提出来，其中最具有通用性和代表性的是体积平均网格，即格子的材料系数由格子内部所有填充的材料根据体积占比进行加权平均得到，根据具体数值计算使用的物理公式不同，还可以在体积平均的基础上改进，例如FDTD的极化体积平均网格则在体积平均的基础上添加了材料分界面的法向量信息。已公开的体积平均网格计算方法多采取近似法，例如格内采样法和表面填充估算法等。格内采样法即在一个格子内均匀的放上采样点，而格子的材料系数就是格内所有采样点处材料系数的平均数，例如LumericalFDTD 软件中的dielectric volume average 便是采用此种方法；表面填充估算法则是计算出平行六面体六个面上的材料填充面积，并使用这些信息估算材料在格子内部的填充体积，例如Lumerical FDTD 软件中的volume average 采用此种类似的方法。传统的格内采样法计算体像素填充体积需要计算的越精确，则需要越多的体像素内采样点，存在网格精度低，计算速度慢、和计算量大等缺点。表面填充估算法则存在计算精度低的缺陷。
发明内容
本发明的目的是针对现有技术的不足，提供一种基于CAD的平行六面体体积平均共形网格剖分方法，该方法可以将三维高斯三角网（闭合三角网）根据体积平均计算方式转换为平行六面体网格，并且此过程中计算的格内体积都是精确的体积。除FDTD以外，该方法可用于任何使用平行六面体网格计算的算法。
本发明解决其技术问题采用的技术方案是：一种基于CAD的平行六面体体积平均共形网格剖分方法，包括以下步骤：
1、使用CAD对物体进行建模，得到描述该物体的三角网模型；对仿真空间进行平行六面体网格划分，得到一系列体像素；
2、读取三角网模型和平行六面体网格，计算出该模型的阶梯网格；遍历三角网模型中的所有三角形，计算每个三角形的包围盒；
3、遍历与包围盒发生交叠碰撞的体像素，使用分离面定理检测包围盒内的体像素与三角形的碰撞情况，如果发生碰撞则将三角形添加到该体像素的三角形集合中；遍历完所有三角形，得到位于目标模型边界处的体像素集合和集合中每一个体像素对应的与其碰撞的三角形集合；
4、遍历体像素集合，并调出与其碰撞的三角形集合，计算三角形留在该体像素内部的部分，使用平行六面体裁剪三角形算法，得到分割曲面，并根据分割曲面与体像素的上顶面的相交情况填补顶面上的多边形，得到遮荫面；使用阴影体积算法计算出不同材料在体像素内的填充体积；计算每个体像素的体积平均等效材料系数，利用计算得到的体像素的体积平均等效材料系数更新阶梯网格中的材料系数，得到平行六面体体积平均共形网格。
进一步地，所述分离面定理具体为：对于任意两个凸多面体，若存在一个分离面，使得单个凸多面体所有顶点在该分离面的同侧，并且该两个凸多面体各自的所有顶点分别位于该分离面的两侧，则该两个凸多面体没有重合，不发生碰撞。
进一步地，所述平行六面体裁剪三角形算法具体包括：首先使用3D-Liang-Barsky算法，将平行六面体作为三维窗口对三角形的每条边进行裁剪，得到留在三维窗口内部的部分，再使用三维窗口的六个边界面分别对三角形进行求交线计算，交线只保留在三维窗口内部的部分，所有裁剪得到的线段都是有向线段，所有的有向线段首尾相连得到多边形；遍历三角形集合后得到多个多边形，这些多边形彼此相连形成一个分割曲面。
进一步地，若分割曲面边界与体像素的上顶面相交，则在体像素的上顶面填补一个多边形，填补的多边形与分割曲面相加得到遮荫面。
进一步地，填补多边形的方法为：计算分割曲面留在体像素上顶面的边集合，将起点端的有向线段反向延长至无穷远，将终点端的有向线段正向延长至无穷远，得到广义多边形，通过Sutherland-Hodgman裁剪算法使用体像素上顶面的四条边依次对广义多边形进行裁剪，裁剪完成后得到填补的多边形。
进一步地，阴影体积算法具体为：将遮荫面拆分为多个平面多边形，每一个平面多边形与其在体像素底面上的投影构成一个不规则柱体，将不规则柱体拆分为多个上顶面倾斜的不规则三棱柱，不规则三棱柱的体积计算见公式为：；
其中，、、为不规则三棱柱的三个高，
s为不规则三棱柱的底面积，
v为不规则三棱柱的体积；将每一个不规则柱体中包含的所有不规则三棱柱的体积相加，得到该多面体的阴影体积；
按照遮荫面上多边形的外法向量区分内外，将位于遮荫面以内的所有不规则柱体的阴影体积添加正号，位于遮荫面以外的所有不规则柱体的阴影体积添加负号；将所有带符号的阴影体积求和，得到遮荫面的阴影体积，当遮荫面的阴影体积为正数时，其为目标模型留在体像素内部的体积，当遮荫面的体积为负数时，目标模型留在体像素内部的体积为体像素体积加上遮荫面的体积。
本发明与现有技术相比，具有如下有益效果：与现有的商业软件对比，在同一模型网格剖分下，本发明所用的时间明显更少。本发明提出的算法在速度上超过LumericalFDTD网格生成器，并且随着网格规模变大其速度优势进一步体现出来，这归结于本发明拥有更少的计算复杂度，并且生成的体积平均共形网格的精度没有来自估算和离散采样上的损失。
附图说明
图1为本发明的基于CAD的平行六面体体积平均共形网格剖分方法的流程示意图；
图2为三角形以及其包围盒示意图；
图3为分离轴定理示意图；
图4为体像素对三角形裁剪示意图；
图5为体像素以及与它相碰撞的所有三角形；
图6为体像素中分割曲面和上顶面填补的多边形示例；其中，（a）阴影部分为由多边形和组成的分割曲面；（b）阴影部分为填补的多边形；
图7为体像素上顶面与分割曲面相交示例；
图8 为使用左、右两边对广义多边形进行裁剪；其中，（a）表示使用左边所在直线进行裁剪；（b）表示使用右边所在直线进行裁剪；
图9为使用下、上两边对广义多边形进行裁剪；其中，（a）表示使用下边所在直线进行裁剪；（b）表示使用上边所在直线进行裁剪；
图10为遮荫面和光线示意图；
图11为不规则多棱柱分解为不规则三棱柱的示意图；
图12为圆环体三角网模型；
图13为圆环体体积平均体像素模型；
图14为玩具飞机三角网模型；
图15为玩具飞机体积平均体像素模型；
图16为圆环体在不同分辨率下网格生成用时对比图（纵坐标为所用时间，横坐标为x轴网格数）。
具体实施方式
下面将结合本发明的实施例和附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。以下对至少一个示例性实施例的描述实际上仅仅是说明性的，绝不作为对本发明及其应用或使用的任何限制。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明公开一种可以将三维高斯三角网（闭合三角网）根据体积平均计算方式转换为平行六面体网格，此方法包括但不限于体积平均平行六面体网格的生成，在计算过程中给出了充足的几何信息可用于计算大多数基于体积平均体像素网格改进的方法，并且此过程中计算的格内填充体积都是精确的体积。在计算速度上与Ansys Lumerical 软件的FDTD模块的体积平均网格生成速度进行了对比，本发明在计算速度上具有明显的优势。除FDTD以外，此方法可用于任何使用平行六面体网格计算的算法。
本发明提出的基于CAD的平行六面体体积平均共形网格剖分方法，流程如图1所示，读取CAD模型（闭合三角网集合）和平行六面体网格，使用光线追踪算法计算出阶梯网格，再使用分离轴定理检测三角形与平行六面体的碰撞，依此计算出与同一个平行六面体相碰撞的所有三角形，得到边界网格以及每个边界网格的三角形集合（与某个格子发生碰撞的三角形集合），之后再计算该三角形集合在平行六面体内的部分，遍历完三角形集合后得到一个可将平行六面体分割为多个部分的曲面（此面下文称其为分割曲面），计算分割曲面留在体像素上顶面的边集合，通过Sutherland-Hodgman裁剪算法得到上顶面上的填补多边形，填补多边形和分割曲面相加得到遮荫面，通过阴影体积算法可以从遮荫面和平行六面体的几何信息中计算出阴影体积，而材料的填充体积可由阴影体积计算得出，最后得到边界网格的材料属性并用其更新阶梯网格得到体积平均共形网格。
本发明的基于CAD的平行六面体体积平均共形网格剖分方法，详细步骤如下：
1、在使用CAD对物体进行建模之后，得到描述该物体的三角网模型，如图12和图14所示，文件形式可以是STL等三角面元格式的文件。对仿真空间进行网格划分，得到一系列体像素，即平行六面体网，此时体像素上的材料系数还是未知的。
2、网格生成器首先读取三角网格和平行六面体网格，使用光线追踪算法计算出该三角网模型的阶梯网格。再遍历所有三角形，计算三角形的包围盒，包围盒是指在直角坐标系中刚好可包围目标（三角形）的最小平行六面体，且每一边都平行于一个坐标平面，如图2所示。
遍历包围盒内以及被包围盒碰到的体像素，并使用分离面定理检测体像素与该三角形的碰撞情况，如果发生碰撞则将三角形添加到该体像素的三角形集合中，遍历完所有三角形便得到位于目标模型边界处的体像素集合和该集合中每一体像素对应的与其发生碰撞的三角形集合。
分离面定理由分离轴定理推导而出。分离轴定理可用于检测两个二维凸多边形是否有重合，如果两个凸多边形没有重合，那么必然存在一条直线（分离轴）将两者分离。如果某条分离轴是真的，垂直该分离轴做一条直线，将两个多边形投影到该直线上，则该两段投影线段不相交。检测是否存在这样一条真的分离轴的方法为：遍历所有潜在的分离轴方向，即两多边形所有边的方向集合，垂直该边的方向做一条任意的直线，通过多边形在该直线上的投影判断潜在分离轴的真假，只要有一条为真，则两凸多边形是不重合的，反之两个多边形重合。如图3所示，线段
A 1
B 1 为矩形在水平线上的投影，线段
C 1
D 1 为三角形在水平线上的投影，很明显线段
A 1
B 1 和
C 1
D 1 发生交叠，故垂直水平方向不存在分离轴。线段
A 2
B 2 和
C 2
D 2 则没有交叠，垂直线段
C 2
D 2 的方向存在直线，即为分离轴，所以矩形和三角形没有重合。同样的，将分离轴定理推广到三维空间，则有分离面定理，可用于检测两个凸多面体是否重合，潜在的分离面是与两个凸多面体任一面平行的面集合，遍历所有潜在分离面，计算凸多面体顶点与该面的关系，如果存在一个面使得单个凸多面体所有顶点在该面同侧，而每个凸多面体自个所有的顶点分别在该面两侧，则该分离面为真。因此，使用分离面定理可计算出与某个体像素（平行六面体网格）发生碰撞的所有三角形。
3、遍历位于目标模型边界处的体像素集合，并调出其三角形集合，这一步需要计算三角形留在体像素内部（包括边界）的部分。使用平行六面体裁剪三角形算法，对于该体像素发生碰撞的三角形进行裁剪，获得分割曲面。
平行六面体裁剪三角形算法具体为：如图4所示，三角形
T 1 的三个顶点分别为
Q 1
、
Q 2
、Q 3 ；且外法向量的方向与顶点顺序满足右手螺旋定则，首先使用3D-Liang-Barsky算法，将平行六面体作为三维窗口对三角形的每条边进行裁剪，得到留在三维窗口内部的部分，再使用三维窗口的六个边界面分别对三角形进行求交线计算，同样交线也只保留在三维窗口内部（包括边界）的部分，注意所有裁剪得到的线段都是有向线段，所有计算得到的有向线段
P 1
P 2 、
P 2
P 3 、
P 3
P 4 、
P 4
P 1 首尾相连得到多边形
P 1
P 2
P 3
P 4 ，记作，即图4中的阴影部分，该多边形边的方向与外法向量依旧满足右手螺旋定则。
在遍历三角形集合后将得到多个2D多边形，这些多边形彼此相连形成一个曲面，称其为分割曲面，如图5所示，分割曲面由多边形和多边形组成，即图中阴影部分。其中，表示三角形被平行六面体（体像素）裁剪后得到的多边形。
由于阴影体积算法的特性，分割曲面最外边界不应该与体像素的上顶面（光线进入体像素的面）相交，如果有相交，如图6中（a）所示，即分割曲面边界与体像素的上顶面相交，则应该如图6中（b）所示补上一个多边形
。填补多边形的方法如下：计算分割曲面留在体像素上顶面的边集合，例如图7中虚线为体像素上顶面的矩形，有向线段和为分割曲面与体像素上顶面的交线集合，将反向延长至无穷远，用向量表示，点在无穷远处，将正向延长至无穷远，用向量表示，点在无穷远处，得到广义多边形，广义多边形的两端表示无穷远处，通过Sutherland-Hodgman裁剪算法使用矩形左、右、下、上四条边依次对广义多边形进行裁剪，裁剪过程分别如图8中（a）和（b）、图9中（a）和（b）所示，广义多边形从依次变为广义多边形、广义多边形、广义多边形和普通多边形，该多边形为裁剪完成后得到的填补多边形，记作。填补的多边形与分割曲面相加得到遮荫面。
4、采用阴影体积算法计算遮荫面的阴影体积；使用体积占比加权平均法计算每个体像素的材料系数，用体像素的材料系数更新阶梯网格中的材料系数，便得到了平行六面体体积平均共形网格，最终输出平行六面体体积平均共形网格。
阴影体积算法如图10所示，一个空间内有一束照向地面的平行光和一个不透明曲面，即遮荫面（由多个平面多边形组成），在遮荫面和地面（体像素底面）之间有一段光线无法进入的空间称为阴影体积。阴影体积的计算方式为将遮荫面拆分为多个平面多边形，每一个平面多边形与其在体像素底面上的投影组成不规则柱体，得到多个上顶面倾斜的不规则柱体，如图11所示，每一个不规则柱体又可以拆分为多个上顶面倾斜的不规则三棱柱，不规则三棱柱的体积计算公式见式（1）：，式（1）；
其中，、、为不规则三棱柱的三个高，
s为不规则三棱柱的底面积，
v为不规则三棱柱的体积。将每一个不规则柱体中包含的所有不规则三棱柱的体积相加，得到该不规则柱体的阴影体积。
理论上，由于遮荫面都是从闭合高斯面中截取出来的部分曲面，因此遮荫面的两侧分别表示了高斯面的内和外，用遮荫面上多边形的外法向量区分高斯面的内、外，如图10中，表示多边形的外向单位法向量，表示多边形的外向单位法向量，表示多边形的外向单位法向量，、和箭头所指方向为高斯面外，故还需要给每个不规则柱体的阴影体积加上一个符号，位于高斯面以内的为正，位于高斯面以外的为负，最终将遮荫面上所有多边形对应的有符号阴影体积相加，可得到遮荫面的阴影体积，即闭合高斯面包围体在该体像素内部的填充体积。如图10所示的多边形
、 和对应的阴影体积均位于高斯面以内，阴影体积均为正，因此，三者之和就是遮荫面的阴影体积。当整个遮荫面的阴影体积为负数时，则说明分割曲面没有触碰到体像素上顶面，填充体积便为体像素的体积加上阴影体积。
由遮荫面的阴影体积可得到目标模型留在体像素内部的体积，该材料系数为ε1，该材料在该体像素中的填充体积为
v 1，用体像素的体积减去
v 1得到该体像素背景空间的体积
v 0 ，背景空间（即目标模型以外空间）的材料系数为ε0，则该体像素的体积平均等效材料系数
ε eff 可由式(2)计算得出：，式（2）；利用计算得到的位于目标模型边界处的体像素的体积平均等效材料系数更新阶梯网格中的材料系数，便得到了体积平均共形网格。
由于FDTD使用的是平行六面体网格，本发明与Lumerical FDTD仿真软件在体积平均共形网格的生成速率方面做对比。具体实施方案为使用相同的STL文件导入模型，涉及圆环体和不规则物体，在相同的计算机环境下进行计算，对比本发明提出的方法与LumericalFDTD仿真软件完成计算所花费的时间。
如图12和图14分别表示由三角网建模得到的圆环体和玩具飞机模型，图13和图15分别为使用本发明方法得到的体像素模型，颜色越深的体像素表示目标模型在该体像素内的体积占比越大，颜色越浅体积占比越小。
使用本发明方法与Lumerical 软件FDTD模块volume average网格剖分技术对同一个圆环体的STL文件进行剖分，圆环体三角面元数量为4900个，剖分空间为立方体，三轴方向网格数相同，圆环体位于剖分空间的正中间，图16的横坐标为立方体边长上的网格数量，纵坐标为完成剖分所占用cpu的时间，仿真环境如下：cpu为Intel(R) Core(TM) i5-9400，Ram为16G，操作系统为Windows11。对比商业软件以及其发布版本：Ansys Lumerical2020 R2 FDTD Solver Version 8.24.2387 (Windows 64bit)，对比软件的共形网格类型为volume average 。由于剖分的网格是用于FDTD计算的，故时间为计算三套极化网格所需的时间。从图16中可以看出，在同一模型网格剖分下，与商业软件Lumerical对比，本发明的方法在速度上超过Lumerical FDTD网格生成器，所用的时间明显更少。并且随着网格规模变大，其速度优势进一步体现出来，在对一个三角面元数量为4900的圆环体，在300*300*300个体像素网格中进行剖分，本发明的剖分速度超过Lumerical FDTD网格生成器的4倍，这归结于本发明拥有更少的计算复杂度，由于本发明从原模型通过精确计算体积信息，所以在计算中也不存在精度损失，生成的体积平均共形网格的精度没有来自估算和离散采样上的损失。
Claims (6)
1.一种基于CAD的平行六面体体积平均共形网格剖分方法，其特征在于，包括以下步骤：
（1）、使用CAD对物体进行建模，得到描述该物体的三角网模型；对仿真空间进行平行六面体网格划分，得到一系列体像素；
（2）、读取三角网模型和平行六面体网格，利用光线追踪算法计算出该模型的阶梯网格；遍历三角网模型中的所有三角形，计算每个三角形的包围盒；
（3）、遍历与包围盒发生交叠碰撞的体像素，使用分离面定理检测包围盒内的体像素与三角形的碰撞情况，如果发生碰撞则将三角形添加到该体像素的三角形集合中；遍历完所有三角形，得到位于目标模型边界处的体像素集合和集合中每一个体像素对应的与其碰撞的三角形集合；
（4）、遍历体像素集合，并调出与其碰撞的三角形集合，计算三角形留在该体像素内部的部分，使用平行六面体裁剪三角形算法，得到分割曲面，并根据分割曲面与体像素的上顶面的相交情况填补顶面上的多边形，得到遮荫面；使用阴影体积算法计算出不同材料在体像素内的填充体积；计算每个体像素的体积平均等效材料系数，利用计算得到的体像素的体积平均等效材料系数更新阶梯网格中的材料系数，得到平行六面体体积平均共形网格。
2.根据权利要求1所述的基于CAD的平行六面体体积平均共形网格剖分方法，其特征在于，所述分离面定理具体为：对于任意两个凸多面体，若存在一个分离面，使得单个凸多面体所有顶点在该分离面的同侧，并且该两个凸多面体各自的所有顶点分别位于该分离面的两侧，则该两个凸多面体没有重合，不发生碰撞。
3.根据权利要求1所述的基于CAD的平行六面体体积平均共形网格剖分方法，其特征在于，所述平行六面体裁剪三角形算法具体包括：首先使用3D-Liang-Barsky算法，将平行六面体作为三维窗口对三角形的每条边进行裁剪，得到留在三维窗口内部的部分，再使用三维窗口的六个边界面分别对三角形进行求交线计算，交线只保留在三维窗口内部的部分，所有裁剪得到的线段都是有向线段，所有的有向线段首尾相连得到多边形；遍历三角形集合后得到多个多边形，这些多边形彼此相连形成一个分割曲面。
4.根据权利要求1所述的基于CAD的平行六面体体积平均共形网格剖分方法，其特征在于：若分割曲面边界与体像素的上顶面相交，则在体像素的上顶面填补一个多边形，填补的多边形与分割曲面相加得到遮荫面。
5.根据权利要求4所述的基于CAD的平行六面体体积平均共形网格剖分方法，其特征在于，填补多边形的方法为：计算分割曲面留在体像素上顶面的边集合，将起点端的有向线段反向延长至无穷远，将终点端的有向线段正向延长至无穷远，得到广义多边形，通过Sutherland-Hodgman裁剪算法使用体像素上顶面的四条边依次对广义多边形进行裁剪，裁剪完成后得到填补的多边形。
6.根据权利要求1所述的基于CAD的平行六面体体积平均共形网格剖分方法，其特征在于，阴影体积算法具体为：将遮荫面拆分为多个平面多边形，每一个平面多边形与其在体像素底面上的投影构成一个不规则柱体，将不规则柱体拆分为多个不规则三棱柱，不规则三棱柱的体积计算公式为：
；
其中，、、为不规则三棱柱的三个高，s为不规则三棱柱的底面积，v为不规则三棱柱的体积；
将每一个不规则柱体中包含的所有不规则三棱柱的体积相加，得到该不规则柱体的阴影体积；
按照遮荫面上多边形的外法向量区分内外，将位于遮荫面以内的所有不规则柱体的阴影体积添加正号，位于遮荫面以外的所有不规则柱体的阴影体积添加负号；将所有带符号的阴影体积求和，得到遮荫面的阴影体积；当遮荫面的阴影体积为正数时，其为目标模型留在体像素内部的体积，当遮荫面的体积为负数时，目标模型留在体像素内部的体积为体像素体积加上遮荫面的体积。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310257023.XA CN115984511B (zh) | 2023-03-17 | 2023-03-17 | 一种基于cad的平行六面体体积平均共形网格剖分方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310257023.XA CN115984511B (zh) | 2023-03-17 | 2023-03-17 | 一种基于cad的平行六面体体积平均共形网格剖分方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115984511A true CN115984511A (zh) | 2023-04-18 |
| CN115984511B CN115984511B (zh) | 2023-06-02 |
Family
ID=85968430
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310257023.XA Active CN115984511B (zh) | 2023-03-17 | 2023-03-17 | 一种基于cad的平行六面体体积平均共形网格剖分方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115984511B (zh) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2927184A1 (fr) * | 2008-02-04 | 2009-08-07 | Univ Joseph Fourier | Procede et dispositif de simulation de contact par images a plan de profondeur. |
| CN106934116A (zh) * | 2017-02-22 | 2017-07-07 | 华南理工大学 | 一种基于遗传算法的骨支架孔径分布控制方法 |
| CN110675496A (zh) * | 2019-09-27 | 2020-01-10 | 武汉中地数码科技有限公司 | 基于三维城市地质模型的网格剖分和可视化方法及其系统 |
| CN112927366A (zh) * | 2021-05-10 | 2021-06-08 | 中南大学 | 一种沥青路面构造深度测量方法 |
| CN113255194A (zh) * | 2021-06-24 | 2021-08-13 | 江苏赛博空间科学技术有限公司 | 一种基于参数曲面的自主智能直六面体剖分方法 |
| CN115758492A (zh) * | 2022-11-21 | 2023-03-07 | 山东大学 | 一种基于三维cad的fdtd共形网格自动生成方法及装置 |
-
2023
- 2023-03-17 CN CN202310257023.XA patent/CN115984511B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2927184A1 (fr) * | 2008-02-04 | 2009-08-07 | Univ Joseph Fourier | Procede et dispositif de simulation de contact par images a plan de profondeur. |
| CN106934116A (zh) * | 2017-02-22 | 2017-07-07 | 华南理工大学 | 一种基于遗传算法的骨支架孔径分布控制方法 |
| CN110675496A (zh) * | 2019-09-27 | 2020-01-10 | 武汉中地数码科技有限公司 | 基于三维城市地质模型的网格剖分和可视化方法及其系统 |
| CN112927366A (zh) * | 2021-05-10 | 2021-06-08 | 中南大学 | 一种沥青路面构造深度测量方法 |
| CN113255194A (zh) * | 2021-06-24 | 2021-08-13 | 江苏赛博空间科学技术有限公司 | 一种基于参数曲面的自主智能直六面体剖分方法 |
| CN115758492A (zh) * | 2022-11-21 | 2023-03-07 | 山东大学 | 一种基于三维cad的fdtd共形网格自动生成方法及装置 |
Non-Patent Citations (3)
| Title |
|---|
| ALAIN RASSINEUX ET AL: "Robust conformal adaptive meshing of complex textile composites unit cells", 《COMPOSITE STRUCTURES》, vol. 279 * |
| TRISTAN CARRIER BAUDOUIN ET AL: "A frontal approach to hex-dominant mesh generation", 《ADVANCED MODELING AND SIMULATION IN ENGINEERING SCIENCES》 * |
| 凌然: "CAE中六面体网格模型的结构简化算法研究", 《中国优秀硕士学位论文全文数据库电子期刊 信息科技辑》, vol. 2021, no. 4 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115984511B (zh) | 2023-06-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11107272B2 (en) | Scalable volumetric 3D reconstruction | |
| US20210174572A1 (en) | Hybrid hierarchy of bounding and grid structures for ray tracing | |
| CN107767453B (zh) | 一种基于规则约束的建筑物lidar点云重构优化方法 | |
| US20110087350A1 (en) | Methods and system for enabling printing three-dimensional object models | |
| CN104966317B (zh) | 一种基于矿体轮廓线的三维自动建模方法 | |
| CN113706713A (zh) | 一种实景三维模型裁剪方法、装置及计算机设备 | |
| CN107464286B (zh) | 三维城市模型中的孔洞修复方法及装置、设备及可读介质 | |
| CN109584357A (zh) | 基于多轮廓线的三维建模方法、装置、系统及存储介质 | |
| Schulz et al. | Interactive visualization of fluid dynamics simulations in locally refined cartesian grids | |
| CN113724401B (zh) | 一种三维模型切割方法、装置、计算机设备和存储介质 | |
| CN107886569B (zh) | 一种基于离散李导数的测度可控的曲面参数化方法及系统 | |
| Bronson et al. | Lattice cleaving: A multimaterial tetrahedral meshing algorithm with guarantees | |
| CN111243094A (zh) | 一种基于点灯法的三维模型精确体素化方法 | |
| Li et al. | On surface reconstruction: A priority driven approach | |
| CN109983509B (zh) | 一种使用几何面的即时布尔运算方法 | |
| CN106649992A (zh) | 舰船与尾迹的网格模型的融合与优化方法 | |
| Li et al. | A sweep and translate algorithm for computing voxelized 3D Minkowski sums on the GPU | |
| Sheng et al. | Accelerated robust Boolean operations based on hybrid representations | |
| Horvat et al. | Ray-casting point-in-polyhedron test | |
| CN115984511B (zh) | 一种基于cad的平行六面体体积平均共形网格剖分方法 | |
| CN116258840A (zh) | 层级细节表示树的生成方法、装置、设备及存储介质 | |
| CN112819108B (zh) | 一种多源异构模型的切片融合方法、系统及存储介质 | |
| CN115935673A (zh) | 一种基于射线追踪的fdtd网格剖分方法及系统 | |
| Lu et al. | Fast and robust generation of city-scale seamless 3D urban models | |
| Yamazaki et al. | Non-manifold implicit surfaces based on discontinuous implicitization and polygonization |
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 |