一种适应复杂几何的模拟仿真软件贴体粒子生成方法技术

本发明专利技术提出了一种适应复杂几何的模拟仿真软件贴体粒子生成方法，属于数值模拟仿真工业软件技术领域。本方法通过合理映射、降维问题域，大幅缩减粒子生成的时间；通过二维平面视角，三维表面的粒子生成简化为二维平面的粒子生成。通过哈希策略，降为一维线性的粒子生成，减少三维表面的贴体粒子的生成时间。同时，二维平面视角保证了生成粒子的贴体性质，对一维粒子排序并去重保证了生成粒子的均匀度。本方法显著加快了粒子生成速度、保证了粒子的贴体性，并能保留几何构型中的天线等纤细结构细节，解决了模拟仿真软件三维表面粒子生成时间长、不贴体、不均匀等问题。不均匀等问题。不均匀等问题。

【技术实现步骤摘要】
一种适应复杂几何的模拟仿真软件贴体粒子生成方法


[0001]本专利技术涉及一种适应复杂几何的模拟仿真软件贴体粒子生成方法，属于数值模拟仿真工业软件


技术介绍

[0002]模拟仿真软件作为工业软件的重要组成部分，是指导产品设计、研发、测试的强力工具。例如，模拟仿真软件可以模拟海啸、汽车涉水等各类场景，帮助应对自然灾害、设计汽车外形等。模拟仿真软件的具体实现手段，主要有网格类方法和粒子类方法两种。其中，粒子类方法是对传统的网格类方法的重要补充，尤其适用于针对大变形、动边界等复杂问题的模拟仿真。
[0003]粒子的空间分布对粒子类方法的计算精度和稳定性等产生很大影响。理想的粒子分布是均匀的，并能准确描述几何与物理信息。目前，现有的粒子生成方法主要有以下几种：1.晶格法。粒子直接生成在三维堆砌中的晶格顶点上。晶格法生成的粒子，理论上均匀性最好。但是，其需要大量计算判断各粒子与几何构型的拓扑关系，包括粒子到构型最短距离的计算、粒子在构型的内外关系判断等。并且，生成的粒子与初始几何构型偏离较大，也没有统一高效的贴体方法，丧失了一定的精确度。
[0004]2.网格法。在已有网格的基础上可以快速得到粒子空间分布。粒子可以直接取网格的节点、边中点、面形心、体心等。依据比较成熟的网格类方法生成粒子，比较直观便捷。但是，其生成均匀的、准确保持几何信息的网格往往需要对网格中的拓扑信息迭代分析，比较耗时，得到适应复杂几何的高质量网格很困难。
[0005]3.物理结合法。在初始的几何构型的基础上，结合物理参数和物理过程对粒子分布进行优化。例如，2021年由胡湘渝团队提出的CAD
‑
BPG，先基于晶格法生成粒子，然后结合level
‑
set方法代入物理方程组，优化粒子的位置。2023年，由孙鹏楠团队开发的SPHydro 也借鉴了这种方法。结合物理的粒子生成方法，可以有多种实现，可以更适应某些物理问题，生成的粒子也比较均匀、贴体。但是，其物理过程的计算需要消耗过多的计算资源。
[0006]综上所述，现有的模拟仿真软件粒子生成方法，存在速度慢、不贴体等多种问题。

技术实现思路

[0007]本专利技术的目的是针对现有技术存在的缺陷和不足，为了解决模拟仿真软件三维表面粒子生成时间长、不贴体、不均匀等技术问题，提出一种适应复杂几何的模拟仿真软件贴体粒子生成方法。
[0008]通过合理映射、降维问题域，可以大幅缩减粒子生成的时间。通过二维平面视角，三维表面的粒子生成，可以简化为二维平面的粒子生成。进一步地，可以通过哈希策略，降为一维线性的粒子生成。这样将大大减少三维表面的贴体粒子的生成时间。同时，二维平面视角保证了生成粒子的贴体性质，对一维粒子排序并去重保证了生成粒子的均匀度。
[0009]为了达到上述目的，本专利技术采取如下技术方案。
[0010]首先对本专利技术涉及的概念内容进行说明。
[0011]1.贴体。是一种点与面的关系，表示点在面上。贴体粒子，指在表面上的粒子。
[0012]2.晶格。是三维空间中的周期性几何结构。常见的有正方体等。
[0013]3.堆砌。是空间中的密铺或镶嵌，由多面体密堆积，或由高维度的晶格紧密堆积而成。
[0014]4.几何构型。描述物件内部结构、外部表面的几何体，由点、线、面、体等组成。
[0015]5.拓扑。包括几何空间中的点、线、面、体间的连接、距离等关系。
[0016]6.网格。是工业实体建模的基本单位，有限元等方法的基本计算单元。由节点、边、面等组成，还可附加其他物理参数信息。
[0017]7.STL文件（STereoLithography，立体光刻）。是一种表示三维表面几何形状的文件格式，由一系列带顶点坐标信息的三角形组成。
[0018]8.VTK文件（Visualization Toolkit，可视化工具箱）。是一种可携带多种信息的数据集文件。可以保存带多种物理参数值的粒子数据。
[0019]一种适应复杂几何的模拟仿真软件贴体粒子生成方法，包括STL三角片粒子生成、晶格哈希和排序去重。具体包括以下步骤：步骤1：设定参数值，包括初设粒子间距margin、STL文件路径、法向量内外方向等。
[0020]步骤2：解析描述几何构型的STL文件，得到各三角面的顶点坐标及法向量信息。
[0021]步骤3：生成满足带系数w
tri
的初设粒子间距margin的各三角面的初始粒子，存入初始粒子数组arrayTris中。如图2所示。
[0022]优选地，系数w
tri
的取值范围为(0.2，1]。
[0023]步骤4：计算哈希初始粒子数组arrayTris中的每个粒子p坐标 (p
x
，p
y
，p
z
)与计算域各轴最小值(p
xmin
，p
ymin
，p
zmin
)的差除以带系数w
lattice
的初设粒子间距值，并将哈希值存入晶格哈希数组arrayHash中。其中，优选地，系数w
lattice
的取值范围为(0.3，1]。p
x
表示粒子p的x轴坐标值，p
y
表示粒子p的y轴坐标值，p
z
表示粒子p的y轴坐标值；p
xmax
、p
xmin
、p
ymax
、p
ymin
、p
zmax
、p
zmin
分别为几何构型中所有点在x、 y、 z轴上的最大值和最小值。
[0024]哈希函数保证了计算域内哈希值相同的粒子必在同一个晶格中。
[0025]步骤5：依照晶格哈希数组arrayHash的大小顺序，对初始粒子数组arrayTris进行排序。
[0026]步骤6：顺序扫描已排序的晶格哈希数组arrayHash，将arrayHash中连续等值的对应的粒子数组arrayTris中的粒子合并，存入目标粒子数组arrayPartis。这样的粒子去重过程只需要线性的时间。并且也很适合并行进一步加速。合并粒子时，要将hash等值的连续粒子的坐标求和平均，并将这些粒子的单位法向量求和。
[0027]其中，hash等值的粒子，是指在同一晶格中的粒子。由三角片生成的初始粒子在公共点和边会更密集，点和边上的初始粒子的权重会更大，更能保留几何构型的形状特征。
[0028]步骤7：偏移（p
xmin
，p
ymin
，p
zmin
）若干次，并重复执行步骤4
‑
6，得到目标粒子数组arrayPartis。
[0029]优选地，若晶格间隔系数w
lattice
<=0.6，则偏移次数为7
‑
8次；若晶格间隔系数0本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种适应复杂几何的模拟仿真软件贴体粒子生成方法，其特征在于，包括STL三角片粒子生成、晶格哈希和排序去重；步骤1：设定参数值，包括初设粒子间距margin、STL文件路径、法向量内外方向；步骤2：解析描述几何构型的STL文件，得到各三角面的顶点坐标及法向量信息；步骤3：生成满足带粒子间距系数w
tri
的初设粒子间距margin的各三角面的初始粒子，存入初始粒子数组arrayTris中；步骤4：计算哈希初始粒子数组arrayTris中的每个粒子p坐标 (p
x
，p
y
，p
z
)与计算域各轴最小值(p
xmin
，p
ymin
，p
zmin
)的差除以带晶格间距系数w
lattice
的初设粒子间距值，并将哈希值存入晶格哈希数组arrayHash中， p
x
表示粒子p的x轴坐标值，p
y
表示粒子p的y轴坐标值，p
z
表示粒子p的y轴坐标值；p
xmax
、p
xmin
、p
ymax
、p
ymin
、p
zmax
、p
zmin
分别为几何构型中所有点在x、 y、 z轴上的最大值和最小值；步骤5：依照晶格哈希数组arrayHash的大小顺序，对初始粒子数组arrayTris进行排序；步骤6：顺序扫描已排序的晶格哈希数组arrayHash，将arrayHash中连续等值的对应的粒子数组arrayTris中的粒子合并，存入目标粒子数组arrayPartis；合并粒子时，将hash等值的连续粒子的坐标求和平均，并将这些粒子的单位法向量求和；其中，hash等值的粒子，指在同一晶格中的粒子；步骤7：偏移（p
xmin
，p
ymin
，p
zmin
）若干次，并重复执行步骤4
‑
6，得到目标粒子数组arrayPartis；步骤8：输出目标粒子arrayPartis数组中的粒子信息。2.如权利要求1所述的一种适应复杂几何的模拟仿真软件贴体粒子生成方法，其特征在于，步骤3中，粒子间距系数w
tri
的取值范围为(0.2，1]。3.如权利要求2所述的一种适应复杂几何的模拟仿真软件贴体粒子生成方法，其特征在于，步骤3中，粒子间距系数w
tri
=0.5。4.如权利要求1所述的一种适应复杂几何的模拟仿真软件贴体粒子生成方法，其特征在于，步骤3包括以下步骤：步骤3.1：沿三角面的一条边AB按初始粒子间距margin/2生成一排初始粒子；步骤3.2：逐步按此三角面边AB对应的高线方向添加一排相距margin/2的粒子，保留在三角面内的粒子；步骤3.3：遍历完高线内的区域，保存初始生成的粒子到初始粒子数组arrayTris中。5.如权利要求1所述的一种适应复杂几何的模拟仿真软件贴体粒子生成方法，其特征在于，步骤4中，晶格间距系数w
latti...

【专利技术属性】
技术研发人员：李逍逍，刘谋斌，
申请(专利权)人：北京大学，
类型：发明
国别省市：