二维颗粒填充模型的生成方法技术

本发明专利技术涉及二维颗粒填充模型的生成方法，该方法包括步骤：S1，基于模型边界构建出二维颗粒填充计算域；S2，确定待生成颗粒的级配，根据该级配下生成的颗粒半径的取值范围、以及所述取值对应的生成概率，构建半径数值集合；S3，根据模型边界与颗粒粒径计算网格数据，并对网格进行划分；S4，采用解析方法确定新生成颗粒的圆心坐标；S5，判断新生成颗粒与相邻网格以及所属网格内的颗粒是否发生重叠，若无重叠发生，则将该新生成颗粒填充至所述模型中，若发生重叠，则重新生成新颗粒S6，重复步骤S4至S5，直至完成模型建立。本发明专利技术的方法为离散元和非连续变形分析法前处理部分提供了一种高效可靠的二维建模方法，极大地提高了计算效率。极大地提高了计算效率。极大地提高了计算效率。

【技术实现步骤摘要】
二维颗粒填充模型的生成方法


[0001]本专利技术属于材料计算力学
，尤其涉及一种二维颗粒填充模型的生成方法。

技术介绍

[0002]二维颗粒模型填充广泛应用于数值计算领域，各类数值方法均离不开二维颗粒模型填充，如离散元(DEM)、光滑粒子流(SPH)、颗粒DDA等数值模拟方法，都离不开二维模型填充算法。如在离散元的数值模拟当中，二维颗粒组成配比更契合实际材料的二维模型可以使得材料的力学、热学、流动性等综合表现达到最真实的效果。一般而言，二维颗粒模型填充中，颗粒级配的控制，有效提高颗粒填充的密实度、均匀性是二维颗粒填充算法的重要目标。
[0003]目前，在离散元数值模拟过程中前处理的任务就是建立离散元模型，对计算模型进行构建是整个离散元分析过程的关键一步。模型的合理与否会直接影响计算结果的精度、计算时间的长短、储存时间的长短以及计算过程能否完成。目前，用于离散元计算模型构建的方法主要包括有动力学方法和几何方法，动力方法包括半径拓展法、重力沉积法、挤压排斥法和各向同性压缩法等，几何方法包括规律排列颗粒流、不规律排列颗粒流和多边形颗粒生成法等。前面的动力方法虽然能够满足高密实度、按照一定级配进行模型构建，但所花费的时间很长；几何方法中颗粒流方法虽然能够较快完成模型构建，但不能满足级配要求，而多边形颗粒生成算法虽然能满足级配要求，计算时间却又大大加长。因此，如何实现同时满足高效率、定级配和高密实度离散元计算模型建立意义重大。

技术实现思路

[0004]本专利技术的目的在于提供一种二维颗粒填充模型的生成方法，旨在高效率完成颗粒填充模型。
[0005]本专利技术提供一种二维颗粒填充模型的生成方法，包括以下步骤：
[0006]S1，基于模型边界构建出二维颗粒填充计算域；
[0007]S2，确定待生成颗粒的级配，根据该级配下生成的颗粒半径的取值范围、以及所述取值对应的生成概率，构建半径数值集合；
[0008]S3，根据模型边界与半径数值集合计算网格数据，并对网格进行划分；
[0009]S4，根据所述半径数值集合确定M号种子颗粒和N号次种子颗粒的半径，根据M号种子颗粒和N号次种子颗粒的半径、圆心坐标值以及边界函数，采用解析方法确定新生成颗粒的圆心坐标；
[0010]S5，根据新生成颗粒的圆心坐标、半径确定其所属网格，结合其所属网格与新生成颗粒的映射关系，判断新生成颗粒与相邻网格以及所属网格内的颗粒是否发生重叠，若无重叠发生，则将该新生成颗粒填充至所述模型中，若发生重叠，则重新生成新颗粒
[0011]S6，重复步骤S4至S5，直至完成模型建立。
[0012]优选地，所述步骤S5中，根据新生成颗粒的圆心坐标、半径确定其所属网格，结合其所属网格与新生成颗粒的映射关系还包括：
[0013]判断新生成颗粒是否超出边界；
[0014]若新生成颗粒已超出边界，则需重新生成新颗粒；若新生成颗粒未超出边界，则判断新生成颗粒与相邻网格以及所属网格内的颗粒是否发生重叠。
[0015]优选地，当新生成颗粒超出边界时，计算获得满足与边界曲线、M号种子颗粒以及N号次种子颗粒同时相切的颗粒圆心坐标和半径值，重新生成新的颗粒、以替换出界的颗粒。
[0016]优选地，根据等式d
l
‑
i
‑
r
i
＝C1确定新生成颗粒是否发生出界，其中：d
l
‑
i
为颗粒圆心到曲线的距离值，r
i
为颗粒的半径值，C1为颗粒到曲线距离与半径的差值；
[0017]当C1≥0时，则判断新生成颗粒未出界；当C1<0时，则判断新生成颗粒会发生出界。
[0018]优选地，所述步骤S5中，还包括：
[0019]S5
‑
01:若新生颗粒与相邻网格以及所属网格内的颗粒发生重叠，则判断当前循环次数是否超过预设次数；
[0020]S5
‑
02:若当前循环次数超过预设次数，则根据网格映射关系，以(N
‑
1)号颗粒代替N号次种子颗粒，
[0021]计算获得满足与M号种子颗粒以及(N
‑
1)号次种子颗粒同时相切的颗粒圆心坐标和半径值，重新生成新的颗粒、以替换重叠的颗粒；
[0022]S5
‑
03：若当前循环次数未超过预设次数，则根据网格映射关系，以(M
‑
1)号颗粒作为种子颗粒，计算获得满足与(M
‑
1)号种子颗粒以及N号次种子颗粒同时相切的颗粒圆心坐标和半径值，重新生成新的颗粒、以替换重叠的颗粒。
[0023]与现有技术相比，本专利技术具有以下有益效果：
[0024](1)本专利技术技术方案可根据需要的颗粒级配，在给定边界域内进行颗粒填充模型建立，并且具有良好的均匀性和足够高的密实度，很大程度上拟合了实际压密土体。
[0025](2)通过网格映射的方法，在新生成颗粒重叠处理中，精简了颗粒重叠判断流程，以种子颗粒和次种子颗粒为依据，进行后续颗粒的生成，极大程度减少了迭代次数，易于大规模颗粒数模型构建。
[0026](3)本专利技术技术方案中的边界处理方法，在给定模型域内新生成颗粒逐个与边界进行相切计算完成填充，使得整个模型完全固定在边界域内，与其他算法相比无需进行压密或多次边界处理。
[0027](4)本专利技术技术方案中的半径数值集合可完成对新生成颗粒半径的控制，通过控制颗粒半径值大小，通过算法可实现多尺度颗粒填充，在进行多尺度问题处理上，效果好、处理过程便捷。
[0028](5)通过几何方法中的解析解法，为离散元和非连续变形分析法前处理部分提供了一种高效可靠的二维建模方法，与其他算法相比(如图6所示)，在保持计算结果相似的情况下，极大地提高了计算效率。
[0029]本专利技术还具有的其他优势请见实施例部分的相关描述。
附图说明
[0030]为了更清楚地说明本专利技术实施例中的技术方案，下面将对实施例或现有技术描述
中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0031]图1是本专利技术二维颗粒填充模型的生成方法的流程示意图；
[0032]图2(a)
‑
(f)是本专利技术二维颗粒填充模型的生成方法中颗粒生成流程示意图；
[0033]图3(a)和图3(b)是本专利技术二维颗粒填充模型的生成方法中生成的新颗粒发生重叠情况的处理示意图；
[0034]图4(a)和图4(b)是本专利技术二维颗粒填充模型的生成方法中边界颗粒生成处理示意图；
[0035]图5是本专利技术二维颗粒填充模型的生成方法中模型完成填充后的效果示意图；
[0036]图6是本专利技术二维颗粒填充模型的生成方法与现有技术中其他方法的填充效果对比示本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种二维颗粒填充模型的生成方法，其特征在于，包括以下步骤：S1，基于模型边界构建出二维颗粒填充计算域；S2，确定待生成颗粒的级配，根据该级配下生成的颗粒半径的取值范围、以及所述取值对应的生成概率，构建半径数值集合；S3，根据模型边界与半径数值集合计算网格数据，并对网格进行划分；S4，根据所述半径数值集合确定M号种子颗粒和N号次种子颗粒的半径，根据M号种子颗粒和N号次种子颗粒的半径、圆心坐标值以及边界函数，采用解析方法确定新生成颗粒的圆心坐标；S5，根据新生成颗粒的圆心坐标、半径确定其所属网格，结合其所属网格与新生成颗粒的映射关系，判断新生成颗粒与相邻网格以及所属网格内的颗粒是否发生重叠，若无重叠发生，则将该新生成颗粒填充至所述模型中，若发生重叠，则重新生成新颗粒；S6，重复步骤S4至S5，直至完成模型建立。2.根据权利要求1所述的二维颗粒填充模型的生成方法，其特征在于，所述步骤S5中，根据新生成颗粒的圆心坐标、半径确定其所属网格，结合其所属网格与新生成颗粒的映射关系还包括：判断新生成颗粒是否超出边界；若新生成颗粒已超出边界，则需重新生成新颗粒；若新生成颗粒未超出边界，则判断新生成颗粒与相邻网格以及所属网格内的颗粒是否发生重叠。3.根据权利要求2所述的二维颗粒填充模型的生成方法，其特征在于，当新生成颗粒超出边界时，计算获得满足与边界曲线、M号种子颗粒以及N号次种子颗粒同时相切的颗粒圆心坐标和半径值，重新生成新的颗粒、以替换出界的颗粒。4.根据权利要求3所述的二维颗粒填充模型的生成方法，其特征在于：(x
i
‑
x
a
)2+(y
i
‑
y
a
)2＝(r
i
+r
a
)2在步骤S4中，根据等式(x
i
‑
x
b
)2+(y
i
‑
y
b
)2＝(r
i
+r
b
)2计算新生成颗粒圆心坐标，其中X
a
、Y
a
为M号种子颗粒的圆心坐标，X
b
、Y
b
为N号次种子颗粒的圆心坐标，X
i
、Y
i
为新生成颗粒的圆心坐标,r
a
、r
b
、r
i
依次为M号种子颗粒、N号次种子颗粒和新生成颗粒的半径值。5.根据权利要求4所述的二维颗粒填充模型的生成方法，其特征在于，当新生成颗粒超出边界时，根据等式(y1‑
A
‑
y1‑
B
)x
i
+(x1‑
B
‑
x1‑
A
)y
i
+x1‑
A
y1‑
B
‑
x1‑
B
y1‑
A
‑
L1r
i
＝0(x
a
‑
x
i
)2+(y
a
‑
y
i
)2‑
(r
a
+r
i
)2＝0(x
b
‑
x
i
)2...

【专利技术属性】
技术研发人员：马建军，丁文洁，李诚豪，陈俊杰，林越翔，黄林冲，杨宏伟，梁禹，
申请(专利权)人：南方海洋科学与工程广东省实验室珠海，
类型：发明
国别省市：