集料三维形态离散元生成方法技术

本发明专利技术公开了一种集料三维形态离散元生成方法，将集料颗粒设定为由若干等粒径的球体规则排列若干层后形成的六面体、五面体或四面体，并且所述六面体、五面体或四面体中至少有三条边相互垂直，并且所述三条边具有一个共同的端点，以所述三条边的长度表征所述六面体、五面体或四面体的三维尺寸特征，以包含所述三条边的三个垂直面中的三个角度表征所述六面体、五面体或四面体的棱角特征。本发明专利技术提供的集料三维形态离散元生成方法，基于实际集料的三维形态统计数据，按照正态分布随机选择集料的尺寸和棱角参数，可以模拟集料颗粒的真实形状、尺寸和棱角，与实际集料建立联系。

【技术实现步骤摘要】

本专利技术涉及浙青混合料三维离散元仿真技术，尤其涉及一种集料三维形态离散元生产方法。
技术介绍
使用离散元方法模拟集料或者浙青混合料的力学性能时，首先需要在离散元程序中重构集料的三维形态。集料的形态特征对集料和浙青混合料的力学行为具有重要影响。因此，准确模拟集料的三维形态非常重要。使用图像处理和X射线扫描技术可以准确重构集料的三维形态。但是，这种方法依赖于实验室设备。在实际操作过程中，需要准备集料或浙青混合料试样。当需要生成不同类型的集料或者浙青混合料的离散元试件时，这种方法不但耗时并且成本很高。在离散元程序中随机生成三维集料的方法可以解决这些问题。但现有的集料随机生成方法只能生成理想化的集料颗粒，无法与实际集料建立联系。
技术实现思路
专利技术目的:为了克服现有技术中存在的不足，本专利技术提供一种基于实际集料三维形态数据的集料三维形态离散元生产方法，利于准确模拟集料或者浙青混合料的力学性倉泛。技术方案:为实 现上述目的，本专利技术采用的技术方案为:，将集料颗粒设定为由若干等粒径的球体规则排列若干层后形成的六面体、五面体或四面体，并且所述六面体、五面体或四面体中至少有三条边相互垂直，并且所述三条边具有一个共同的端点，以所述三条边的长度表征所述集料颗粒的三维尺寸特征，以包含所述三条边的三个垂直面中的三个角度表征所述集料颗粒的棱角特征。将六面体作为一种基本结构来说，所述六面体至少有一对相互平行的面，所述相互平行的一对面中的一个面的两条边为所述三条相互垂直的边中的两条,所述三条相互垂直的边中的另一条两端分别连接所述相互平行的一对面；所述相互平行的一对面中，至少有一个面中的一条边与另一个面中的一条边平行。将五面体和四面体作为六面体的特殊情况进行处理，所述五面体或四面体，包含三条相互垂直的边的三个面中，有一个、两个或者全部为三角形。具体的，根据实际集料颗粒三维尺寸的均值和标准差、按照正态分布计算所述三条边的长度值，根据实际集料颗粒所有面中最大角度的均值和标准差、按照正态分布计算所述三个角度的角度值。具体的，所述等粒径的球体规则排列的层数，由所述三条边中的表示高的边的长度除以球体直径后取整得到；所述每一层中球体的行数和列数由所述三个垂直面的控制边长和控制角度计算获得。有益效果:本专利技术提供的，基于实际集料的三维形态统计数据，按照正态分布随机选择集料的尺寸和棱角参数，可以模拟集料颗粒的真实形状、尺寸和棱角，与实际集料建立联系。附图说明图1为六面体的集料颗粒示意图；图2为本专利技术方法的流程图。具体实施例方式下面结合附图对本专利技术作更进一步的说明。，将集料颗粒设定为由若干等粒径的球体规则排列若干层后形成的六面体、五面体或四面体，并且所述六面体、五面体或四面体中至少有三条边相互垂直，并且所述三条边具有一个共同的端点，以所述三条边的长度表征所述集料颗粒的三维尺寸特征，以 包含所述三条边的三个垂直面中的三个角度表征所述集料颗粒的棱角特征。将六面体作为一种基本结构来说，所述六面体至少有一对相互平行的面，所述相互平行的一对面中的一个面的两条边为所述三条相互垂直的边中的两条,所述三条相互垂直的边中的另一条两端分别连接所述相互平行的一对面；所述相互平行的一对面中，至少有一个面中的一条边与另一个面中的一条边平行。将五面体和四面体作为六面体的特殊情况进行处理，所述五面体或四面体，包含三条相互垂直的边的三个面中，有一个、两个或者全部为三角形。具体的，根据实际集料颗粒三维尺寸的均值和标准差、按照正态分布计算所述三条边的长度值，根据实际集料颗粒所有面中最大角度的均值和标准差、按照正态分布计算所述三个角度的角度值。具体的，所述等粒径的球体规则排列的层数，由所述三条边中的表示高的边的长度除以球体直径后取整得到；所述每一层中球体的行数和列数由所述三个垂直面的控制边长和控制角度计算获得。结合图1，本案对集料颗粒进行了如下假定:(I)集料颗粒为由若干等粒径的球体规则排列若干层后形成的六面体、五面体或四面体；(2)对六面体的集料颗粒具体做出如下假定:(a)至少有一对相互平行的面，如图1中的面AB⑶和面EFGH ；(b)在相互平行的一对面中，至少有一个顶点处的三条边相互垂直，如图中的顶点A (边AB、边AD和边AE)和顶点E (边EA、边EF和边EH)；(c)在相互平行的一对面中，至少有一对边相互平行，如图1中的边AB和边⑶，以及边EF和边GH ；(3)五面体和四面体为六面体的特殊情况，即包含顶点A和顶点E的三个面中有一个或两个或全部为三角形。下面以六面体的集料颗粒为例进行具体说明。如图1所示，通过边AB、边AD和边AE的长度表征六面体的三维尺寸特征，分别记为1ab、1ad和Iae ;通过角度adh、角度ABF和角度ABC的角度来表征六面体的棱角特征，分别记为Z ADH、Z ABF和Z ABC ;通过变化I 、U、1AE、Z ADH、Z ABF和Z ABC六个参数可以生成具有不同三维尺寸和棱角的集料颗粒。根据实际集料颗粒三维尺寸(长轴长度、中轴长度、短轴长度)的均值和标准差、按照正态分布可以计算IabUai^P Iae的值，根据实际集料颗粒所有面中最大角度的均值和标准差、按照正态分布可以计算Z ADH、Z ABF和Z ABC的值。集料颗粒由若干等粒径的球体规则排列若干层后形成，其中球的层数取决于三条边中的表示高的边AE，通过Iae除以球体直径后取整即可得到；每一层中球体的行数由Iad和Z ADH控制，每一层中球体的列数由Iab、角度Z ABC和角度Z ABF控制。利用离散元软件中的内置程序语言能够编织循环算法生成集料颗粒。基于上述说明，对六面体的集料颗粒给出如下的具体实例应用。1、三维尺寸和棱角参数的计算所示图1中，对边AB、边AD和边AE取长轴长度、中轴长度或短轴长度的情况应该是随机的，本案采用如下算法:首先，由系统生成一个服从均匀分布U(0，I)的随机数μ C1;若μ C1在范围内，则Iab取长轴长度，Iad取中轴长度，Iae取短轴长度；若口0在(1/3，2/3)范围内，则Iad取长轴长度，Iae取中轴长度，Iab取短轴长度；若μ ο在范围内，则Iae取长轴 长度，Iab取中轴长度，Iad取短轴长度。 设集料颗粒的长轴长度为1_，中轴长度为Imid，短轴长度为Imin，所有面中最大顶角为amax，则生成集料颗粒时lmax、lmid、lmin和Gmax以及Z ADH、Z ABF和Z ABC可以分别由式(I)到式(7)计算得到。lfflax=Mfflax+ulX Omax (I)lfflid=Mfflid+u2X Omid (2)lmin=Mmin+u3X omin (3)amax=Ma+u4X σ a (4)Z ADH= Z ADE+u5 X ( a 眶-Z ADE) (5 )Z ABF= Z ABE+u6 X ( a 眶-Z ABE) (6 )Z ABC=max { Z ABD, Z EFH} +u7X ( a —-max { Z ABD, Z EFH}) (7)式中，Mmax、Mmid、MmiI^PMa分别为长轴长度、中轴长度、短轴长度和最大角度的均值，σ max> σ mid> 0_和σ α分别为长轴长度、中轴长度、短轴长度和最大角度的标准差，u本文档来自技高网...

【技术保护点】
集料三维形态离散元生成方法，其特征在于：将集料颗粒设定为由若干等粒径的球体规则排列若干层后形成的六面体、五面体或四面体，并且所述六面体、五面体或四面体中至少有三条边相互垂直，并且所述三条边具有一个共同的端点，以所述三条边的长度表征所述集料颗粒的三维尺寸特征，以包含所述三条边的三个垂直面中的三个角度表征所述集料颗粒的棱角特征。

【技术特征摘要】

【专利技术属性】
技术研发人员：黄晓明，张东，
申请(专利权)人：东南大学，
类型：发明
国别省市：