一种基于真实颗粒形状的离散元可破碎颗粒模型建模方法技术

本发明专利技术提供一种基于真实颗粒形状的离散元可破碎颗粒模型建模方法，所述方法包括先扫描真实颗粒获取颗粒轮廓，再基于球面谐波分析法重构颗粒模型，然后用半径扩大法在颗粒模型内部填充球形子颗粒，再以颗粒模型初步Voronoi划分，再将颗粒轮廓构成的节理组切割单个Voronoi多面体，然后循环上一步直至切割所有的Voronoi多面体，最后由切割后的Voronoi组装成可破碎颗粒模型。本发明专利技术所述方法构建的颗粒模型与真实颗粒模型几乎一样。其计算原理简单，计算程序简明高效，可有效地反映实际颗粒轮廓特征。为进一步通过数值仿真了解颗粒材料破碎行为提供了有效的技术手段，也有助于土建工程领域中对颗粒材料的设计。建工程领域中对颗粒材料的设计。建工程领域中对颗粒材料的设计。

【技术实现步骤摘要】
一种基于真实颗粒形状的离散元可破碎颗粒模型建模方法


[0001]本专利技术属于岩土工程领域，具体涉及一种基于真实颗粒形状的离散元可破碎颗粒模型建模方法。

技术介绍

[0002]颗粒材料因具有诸多优点广泛应用于堆石坝、铁道、公路、机场跑道等土建工程领域。作用于颗粒的外部荷载超过颗粒自身强度后，颗粒发生破碎。颗粒破碎现象广泛存在于各类工程的建设和使用环节，且显著改变了颗粒材料的宏细观力学特性。由于试验设备的限制和监测手段的不足，现有的物理试验难以实时观测颗粒材料内部的破碎现象。因此，难以从颗粒尺度理解各类因素对颗粒破碎的影响。
[0003]基于上述试验的不足，离散元法成为研究颗粒破碎的有力工具。离散元法模拟颗粒破碎大致可分为两类方法：一类是替换颗粒法，当颗粒材料中某个颗粒所受荷载达到某个破碎标准时该颗粒替换为具有一定拓扑关系的若干小颗粒。另一类方法是黏结胞元法，将颗粒模拟为通过黏结力黏聚在一起的胞元集合体。替换颗粒法虽然计算负荷小，运行速度快，但是其破碎标准和破碎模式尚未有明确论断，且破碎前的母颗粒与破碎后的子颗粒均为球形，无法体现颗粒的形状特征。因此，黏结胞元法成为离散元模拟颗粒破碎最常采用的方法。
[0004]黏结胞元法中，胞元有多种形式，包括球形、四面体、不规则多面体。其中，球形颗粒是最初采用的方法，该方法存在如下问题：(1)破碎时无法保证颗粒体积守恒；(2)外部荷载在颗粒内部传递时往往转化为拉伸与转动，与实际情况不符；(3)球形胞元填充颗粒边界的效果并不理想，若要较好的填充效果必须使用大量胞元，从而极大增加计算负荷。为弥补球形胞元的缺陷，四面体或不规则多面体被用作组成颗粒的胞元。二者均可满足颗粒破碎过程中的体积守恒。其中，模型划分为四面体的技术成熟，商业化软件众多，成为模拟颗粒破碎的常用方法。已有研究表明，四面体胞元组成的颗粒其破碎路径更为平滑，而不规则多面体的破碎路径更为曲折，后者可更好地反映颗粒材料破碎的本质。此外，也有研究发现许多岩石(如大理石)具有典型的晶体结构，用不规则多面体表示更为恰当。Voronoi(又称为泰森、冯洛诺伊或维诺)多面体是一种不规则的凸多面体，近些年逐渐作为胞元应用于脆性材料破碎的模拟研究。目前，各类开源程序(如Neper、Voro++)或商业化软件(如PFC3D、3DEC、Matlab)可将简单几何模型划分为Voronoi多面体集合体。由于真实颗粒形状复杂、表明凹凸不平，上述各类方法无法直接将之进行Voronoi划分。已有研究表明颗粒形状是影响颗粒破碎的重要因素。仅通过上述各类途径划分简单模型无法满足数值模拟的需要。因此，有必要提出一种基于真实颗粒形状的离散元可破碎颗粒模型建模方法，综合解决上述问题。

技术实现思路

[0005]本专利技术提供了一种基于真实颗粒形状的离散元可破碎颗粒模型建模方法，以解决
现有离散元构建的可破碎颗粒模型难以准确反映真实颗粒破碎的问题，进而有助于石坝、铁道、公路、机场跑道等土建工程领域中对颗粒材料的设计。
[0006]为了达到上述目的，本专利技术的具体技术方案如下：
[0007]本专利技术提供一种基于真实颗粒形状的离散元可破碎颗粒模型建模方法，所述方法包括如下步骤，
[0008]a.挑选真实颗粒，通过扫描法或图像法获得真实颗粒的三维模型；
[0009]b.选用球谐函数分析法对颗粒模型进行重构；
[0010]c.将重构后的颗粒模型导入离散元程序，并转化为起到约束作用的墙体；随后，在墙体内部生成若干互不接触的球形小粒径子颗粒，再以半径扩大法逐步放大球形子颗粒的半径以填充颗粒内部，当达到目标孔隙率时停止；
[0011]d.以填充完毕的子颗粒型心为种子，将一长方体划分为Voronoi多面体集合体，该长方体是以颗粒模型的轮廓在X、Y、Z轴向三个维度的上下限为边界；对Voronoi多面体集合体中的多个Voronoi多面体进行编号处理，并保存模型；
[0012]e.再次在离散元程序中输入颗粒轮廓并转化为节理组，同时保留轮廓；输入保存的模型，提取编号为1的Voronoi多面体；采用节理组切割该多面体，根据颗粒轮廓判断切割后哪一部分几何体位于颗粒轮廓内部，哪一部分位于颗粒轮廓外部；将轮廓内部所有几何体的顶点信息以软件识别的指令形式导出；
[0013]f.循环进行上一步骤，按照编号对所有Voronoi多面体进行切割，并将顶点信息以软件识别的指令形式保存在同一文件；
[0014]g.运行上一步生成的文件，生成的所有多面体即可拼接为可破碎颗粒模型，胞元之间间赋予粘聚力之后便可模拟颗粒破碎。
[0015]因步骤a中组成颗粒模型的网格数量特别多，通常达到数百万个，为提高后续建模效率，因而步骤b中选用球谐函数分析法对颗粒模型进行重构。
[0016]简单说来，本专利技术提供一种基于真实颗粒形状的离散元可破碎颗粒模型建模方法，所述方法包括先扫描真实颗粒获取颗粒轮廓，再基于球面谐波分析法重构颗粒模型，然后用半径扩大法在颗粒模型内部填充球形子颗粒，再以颗粒模型初步Voronoi划分，再将颗粒轮廓构成的节理组切割单个Voronoi多面体，然后循环上一步直至切割所有的Voronoi多面体，最后由切割后的Voronoi组装成可破碎颗粒模型。
[0017]在一种具体的实施方式中，所述步骤a中，颗粒类型为碎石或卵石；且该步骤包括采用三维物体扫描仪扫描颗粒，获取颗粒轮廓信息，并通过点云建模技术和拼接技术生成与真实形状一致的颗粒模型；该颗粒模型由大量拓扑连接的三角网格组成，并保存为“STL”格式文件。
[0018]在一种具体的实施方式中，所述步骤b中，球谐函数分析法中三角网格数量为1280个以上，且球谐阶数为15阶以上；优选所述三角网格数量为1280～2000个，且球谐阶数为15阶。
[0019]本专利技术步骤b中，网格数量为1536，球谐阶数为15时，可较为全面地反映颗粒轮廓，也消除了轮廓中极为细微的凹凸点，有助于快速建模和后续模拟。当然，也可通过其他方式减少网格数量，如犀牛软件等，但这些方法获取的轮廓面网格随机性高，实用性差。
[0020]在一种具体的实施方式中，所述步骤c中，多个球形小粒径子颗粒之间为等粒径或
为非等粒径，优选为等粒径；且提前确定所述子颗粒的总数量为与模型划分成的Voronoi多面体的数量一致；优选所述子颗粒的总数量为100～2000个，且优选子颗粒在粒径扩大时，每次扩大至原粒径的1.02～1.08倍。
[0021]在一种具体的实施方式中，所述步骤c中，所述目标孔隙率为子颗粒的总体积与碎石模型体积的比值为0.6～0.8之间的某个数。
[0022]在一种具体的实施方式中，所述步骤e中，由1280个以上的节理组成的节理组按照尺寸由大到小的顺序切割Voronoi多面体；切割后产生大量凸多面体，基于放射法通过颗粒轮廓判断哪些在轮廓内部。
[0023]在一种具体的实施方式中，所述步骤f中，每一个Voronoi多面体切割完毕后，轮廓内部所有顶点信息以生成一个块体的本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于真实颗粒形状的离散元可破碎颗粒模型建模方法，所述方法包括如下步骤，a.挑选真实颗粒，通过扫描法或图像法获得真实颗粒的三维模型；b.选用球谐函数分析法对颗粒模型进行重构；c.将重构后的颗粒模型导入离散元程序，并转化为起到约束作用的墙体；随后，在墙体内部生成若干互不接触的球形小粒径子颗粒，再以半径扩大法逐步放大球形子颗粒的半径以填充颗粒内部，当达到目标孔隙率时停止；d.以填充完毕的子颗粒型心为种子，将一长方体划分为Voronoi多面体集合体，该长方体是以颗粒模型的轮廓在X、Y、Z轴向三个维度的上下限为边界；对Voronoi多面体集合体中的多个Voronoi多面体进行编号处理，并保存模型；e.再次在离散元程序中输入颗粒轮廓并转化为节理组，同时保留轮廓；输入保存的模型，提取编号为1的Voronoi多面体；采用节理组切割该多面体，根据颗粒轮廓判断切割后哪一部分几何体位于颗粒轮廓内部，哪一部分位于颗粒轮廓外部；将轮廓内部所有几何体的顶点信息以软件识别的指令形式导出；f.循环进行上一步骤，按照编号对所有Voronoi多面体进行切割，并将顶点信息以软件识别的指令形式保存在同一文件；g.运行上一步生成的文件，生成的所有多面体即可拼接为可破碎颗粒模型，胞元之间间赋予粘聚力之后便可模拟颗粒破碎。2.根据权利要求1所述的方法，其特征在于，所述步骤a中，颗粒类型为碎石或卵石；且该步骤包括采用三维物体扫描仪扫描颗粒，获取颗粒轮廓信息，并通过点云建模技术和拼接...

【专利技术属性】
技术研发人员：方传峰，聂志红，贾明涛，龚健，齐群，刘顺凯，
申请(专利权)人：中南大学，
类型：发明
国别省市：