一种复杂块石的闵可夫斯基和几何形体生成方法及系统技术方案

本发明专利技术公开了一种复杂块石的闵可夫斯基和几何形体生成方法及系统，涉及工程材料土石混合物模拟参数研究技术领域，包括：根据预设的块石信息生成所有的复杂块石模型；向外扩展块石模型，得到闵可夫斯基和几何形体；读取闵可夫斯基和几何形体的节点和三角面信息，修复几何形体，得到修复后的模型；重新划分网格，优化修复后的模型，得到复杂块石的高质量闵可夫斯基和几何形体。本发明专利技术能够更快速高效地生成复杂块石的闵可夫斯基和几何形体，并解决了可能出现的不封闭的情况，通过后续重新划分网格得到具有良好均匀性的模型，可用于后续有限元模拟，提高相关数值模拟的准确度。提高相关数值模拟的准确度。提高相关数值模拟的准确度。

【技术实现步骤摘要】
一种复杂块石的闵可夫斯基和几何形体生成方法及系统


[0001]本专利技术涉及工程材料土石混合物模拟参数研究
，更具体的说是涉及一种复杂块石的闵可夫斯基和几何形体生成方法及系统。

技术介绍

[0002]目前，通过数值方法研究土石混合物的力学参数，可以很好的克服现实试验参数研究中试样变形、砂浆与块石离析、尺寸存在误差等问题。然而，数值方法类的研究存在一个无可避免的问题：模型建立问题，即根据试样内部参数(如：土石混合物中的块石级配、块石的形状、块石的分布情况、块石之间最小间距等)生成符合数值试验要求的数值模型试块。
[0003]通过有限单元法来模拟块石碰撞生成土石混合物块石模型是一种高效的方法，在进行块石模拟碰撞时，需要对模拟的模型添加闵可夫斯基和，以此来保证块石间的最小间距，但是将块石模型向外扩展时存在一些问题：当块石模型较为复杂时，扩展后的网格模型会出现不封闭以及网格质量很差不能准确进行模拟的问题。
[0004]因此，如何提供一种复杂块石的闵可夫斯基和几何形体生成方法及系统是本领域技术人员亟需解决的问题。

技术实现思路

[0005]有鉴于此，本专利技术提供了一种复杂块石的闵可夫斯基和几何形体生成方法及系统，解决现有技术中无法快速生成复杂块石闵可夫斯基和几何形体以及生成的模型网格质量不达标的问题。
[0006]为了实现上述目的，本专利技术采用如下技术方案：一种复杂块石的闵可夫斯基和几何形体生成方法，具体包括以下步骤：
[0007]步骤1：根据预设的块石信息生成所有的复杂块石模型；
[0008]步骤2：向外扩展所述块石模型，得到闵可夫斯基和几何形体；
[0009]步骤3：读取所述闵可夫斯基和几何形体的节点和三角面信息，修复所述几何形体，得到修复后的模型；
[0010]步骤4：重新划分网格，优化所述修复后的模型，得到复杂块石的高质量闵可夫斯基和几何形体。
[0011]优选的，步骤1中，所述预设的块石信息包括：各块石粒径区间、块石长短轴径比信息、块石半径波动参数和块石角度波动参数；
[0012]根据级配曲线划分不同的块石粒径区间；根据块石粒径区间占比控制对应区间内块石粒径的体积与总块石体积的比值；
[0013]根据所述块石半径波动参数和所述块石角度波动参数控制块石均匀性；
[0014]根据所述块石长短轴径比将生成后的块石沿着Z轴伸长相应的比例，用于模仿长径块石。
[0015]根据所述块石半径波动参数和所述块石角度波动参数控制块石均匀性；
[0016]优选的，步骤2中，所述向外扩展所述块石模型的具体过程包括：采用闵可夫斯基和算法，以各个块石模型每个面的法向量向块石外部扩展距离d，得到闵可夫斯基和几何形体每个面的法向量及块石节点坐标信息；其中，d＝D/2，D表示相邻块石之间最小间距。
[0017]优选的，步骤3中，所述修复所述几何形体，包括：调用Python里的pymeshfix库修复扩展后的块石模型中的不封闭面。
[0018]优选的，步骤4中，优化所述修复后的模型包括：对所述修复后的模型设定更多的节点个数和三角面个数，重新划分网格；使得重新划分后的网格模型具有良好的均匀性，最后得到复杂块石高质量闵可夫斯基和几何形体。
[0019]优选的，还包括：步骤5、根据所有的复杂块石的高质量闵可夫斯基和几何形体，采用离散单元法模拟碰撞，生成块石模型随机分布的土石混合物模型。
[0020]优选的，一种复杂块石的闵可夫斯基和几何形体生成系统，包括：
[0021]生成单元，用于根据预设的块石信息生成复杂块石模型；
[0022]扩展单元，用于扩展所述块石模型，得到闵可夫斯基和几何形体；
[0023]修复单元，用于读取所述闵可夫斯基和几何形体的节点和三角面信息，并修复所述几何形体；
[0024]优化单元，用于重新划分网格，优化修复后的块石模型；
[0025]模拟碰撞单元，用于对所有优化后的块石模型采用离散单元法模拟碰撞；
[0026]存储单元，用于存储所述生成单元、所述扩展单元、所述修复单元和所述优化单元处理后的各块石模型信息。
[0027]优选的，所述生成单元中，预设的块石信息包括：各块石粒径区间、块石长短轴径比信息、块石半径波动参数和块石角度波动参数；
[0028]根据级配曲线划分不同的块石粒径区间；
[0029]所述块石半径波动参数和所述块石角度波动参数用于控制块石均匀性；
[0030]所述块石长短轴径比用于将生成后的块石沿着Z轴伸长相应的比例。
[0031]经由上述的技术方案可知，与现有技术相比，本专利技术提供了一种复杂块石的闵可夫斯基和几何形体生成方法及系统，根据预设的块石信息生成所有的复杂块石模型；向外扩展块石模型，得到闵可夫斯基和几何形体；读取闵可夫斯基和几何形体的节点和三角面信息，修复几何形体，得到修复后的模型；重新划分网格，优化修复后的模型，得到复杂块石的高质量闵可夫斯基和几何形体。本专利技术能够更快速高效地生成复杂块石的闵可夫斯基和几何形体，并解决了可能出现的不封闭的情况，通过后续重新划分网格得到具有良好均匀性的模型，可用于后续有限元模拟，提高相关数值模拟的准确度。
附图说明
[0032]为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据提供的附图获得其他的附图。
[0033]图1为本专利技术提供的方法流程示意图；
[0034]图2为本专利技术实施例生成的复杂块石模型图；
[0035]图3为本专利技术实施例添加闵可夫斯基和并修复后的扩展块石模型；
[0036]图4为本专利技术实施例将修复后的块石模型重新划分网格后的模型；
[0037]图5为本专利技术实施例导入离散元软件后的模型图；
[0038]图6为本专利技术实施例采用离散单元法模拟碰撞过程图；
[0039]图7为本专利技术实施例实现块石重分布后的三维土石混合物模型；
[0040]图8为本专利技术实施例对三维土石混合物模型进行网格划分后的网格模型；
[0041]图9为本专利技术提供的系统结构示意图。
具体实施方式
[0042]下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。
[0043]本专利技术实施例公开了一种复杂块石的闵可夫斯基和几何形体生成方法，如图1所示，具体包括以下步骤：
[0044]步骤1：根据预设的块石信息生成所有的复杂块石模型；
[0045]具体的，所述预设的块石信息包括：各块石粒径区间、块石长短轴径比信本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种复杂块石的闵可夫斯基和几何形体生成方法，其特征在于，具体包括以下步骤：步骤1：根据预设的块石信息生成所有的复杂块石模型；步骤2：向外扩展所述块石模型，得到闵可夫斯基和几何形体；步骤3：读取所述闵可夫斯基和几何形体的节点和三角面信息，修复所述几何形体，得到修复后的模型；步骤4：重新划分网格，优化所述修复后的模型，得到复杂块石的高质量闵可夫斯基和几何形体。2.根据权利要求1所述的一种复杂块石的闵可夫斯基和几何形体生成方法，其特征在于，步骤1中，所述预设的块石信息包括：各块石粒径区间、块石长短轴径比信息、块石半径波动参数和块石角度波动参数；根据级配曲线划分不同的块石粒径区间；根据所述块石半径波动参数和所述块石角度波动参数控制块石均匀性；根据所述块石长短轴径比将生成后的块石沿着Z轴伸长相应的比例。3.根据权利要求1所述的一种复杂块石的闵可夫斯基和几何形体生成方法，其特征在于，步骤2中，所述向外扩展所述块石模型的具体过程包括：采用闵可夫斯基和算法，以各个块石模型每个面的法向量向块石外部扩展距离d，得到闵可夫斯基和几何形体每个面的法向量及块石节点坐标信息；其中，d＝D/2，D表示相邻块石之间最小间距。4.根据权利要求1所述的一种复杂块石的闵可夫斯基和几何形体生成方法，其特征在于，步骤3中，所述修复所述几何形体，包括：调用Python里的pymeshfix库修复扩展后的块石模型中的不封闭面。5.根据权...

【专利技术属性】
技术研发人员：孟庆祥，范凯丰，陈佳敏，汪炜，薛浩宇，宋杭天，
申请(专利权)人：河海大学，
类型：发明
国别省市：