CN110245457B - 一种基于连续均匀概率确定颗粒粒径的目标级配模型创建方法 - Google Patents

一种基于连续均匀概率确定颗粒粒径的目标级配模型创建方法 Download PDF

Info

Publication number
CN110245457B
CN110245457B CN201910572878.5A CN201910572878A CN110245457B CN 110245457 B CN110245457 B CN 110245457B CN 201910572878 A CN201910572878 A CN 201910572878A CN 110245457 B CN110245457 B CN 110245457B
Authority
CN
China
Prior art keywords
particles
target space
particle
particle size
target
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201910572878.5A
Other languages
English (en)
Other versions
CN110245457A (zh
Inventor
孙伟宸
吴恺
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Individual
Original Assignee
Individual
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Individual filed Critical Individual
Priority to CN201910572878.5A priority Critical patent/CN110245457B/zh
Publication of CN110245457A publication Critical patent/CN110245457A/zh
Application granted granted Critical
Publication of CN110245457B publication Critical patent/CN110245457B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/02—Investigating particle size or size distribution
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Dispersion Chemistry (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明涉及一种基于连续均匀概率确定颗粒粒径的目标级配模型创建方法，引入连续均匀概率，依据累计颗粒数量百分比随粒径区间的变化关系，在全部粒径区间范围内选取颗粒粒径进行颗粒投放；与现行静态投放方法相比，能够大幅提高模型最大填充率；与现行瞬态填充方法相比，能够避免各级配区间颗粒在目标投放空间上、位置分布规律与投放顺序的相关性，确保各级配区间颗粒在目标投放空间上的位置、均服从目标分布规律，同时能够在整个投放过程中，使颗粒粒径不间断地满足目标级配分布，确保瞬态填充过程的连续性。

Description

一种基于连续均匀概率确定颗粒粒径的目标级配模型创建 方法
技术领域
本发明涉及一种基于连续均匀概率确定颗粒粒径的目标级配模型创建方法，属于岩土工程及材料科学数值模拟技术领域。
背景技术
颗粒级配是决定材料特性的重要指标，随着有限元法(FEM)、有限差分法(FDM)、离散元法(DEM)、有限元离散元耦合方法(FDEM)、数值流形方法(NMM)、非连续变形分析方法(DDA)、物质点法(MPM)、格子玻尔兹曼方法(LBM)、光滑粒子动力学(SPH)、近场动力学(PD)等数值模拟技术在岩土工程及材料科学等领域的普及应用，如何创建满足目标级配分布的数值模型来表征现实材料成为无法回避的科研课题，受到越来越多科研工作者的重视。
数值模拟技术中现行生成满足目标级配分布模型的方法，大体分为空间维度上的静态投放和时间维度上的瞬态填充。现行静态投放法依据目标级配区间在目标投放空间内对应分批次投放颗粒，同一级配区间投放完毕才进行下一级配区间投放，直至所有级配区间投放完毕，典型应用如路基碎石土静态加载数值模型的创建。而现行瞬态填充法则是上述多个静态投放过程在时间维度上的串联，典型应用如级配矿料的动态入料传输过程模拟。
现行静态投放法又可细分为允许颗粒重叠和不允许颗粒重叠两种情况。
允许颗粒重叠的情况。模型创建之初允许颗粒之间存在大量重叠，后经一定求解步颗粒逐渐弹开。此种情况模型最大填充率不受限制，各级配区间颗粒在目标投放空间上的位置分布规律与投放顺序无关。然而，此种情况能够适用的数值模拟方法有限，仅适用于基于离散元方法创建满足颗粒级配的数值模型。
不允许颗粒重叠的情况，模型创建之初不允许颗粒之间存在重叠，颗粒投放过程中通过重叠判断舍弃与已投放颗粒重叠的新投放颗粒；此种情况下，模型最大填充率受到一定限制，最大填充率与各级配区间颗粒投放顺序相关，按粒径区间由大到小的顺序依次投放颗粒能够获得最大填充率。然而，即便是按粒径区间由大到小的顺序依次投放颗粒获得的最大填充率模型，也往往难以满足现实工程中对高填充率问题的数值模拟需要。另外，当目标填充率超过允许最大填充率时，目标投放空间无法确保所有级配区间均投放成功，只要有一个级配区间投放失败，则整个模型颗粒粒径不满足目标级配分布。各级配区间颗粒在目标投放空间上的位置分布规律与投放顺序相关，投放越晚越难服从目标分布规律。
现行瞬态填充方法是多个静态投放过程在时间维度上的串联，即瞬态填充过程分多个阶段进行，每隔一定时间在颗粒目标投放空间内完成一次严格满足颗粒目标级配的静态投放过程。现行瞬态填充方法需要严格控制每一静态投放阶段中的颗粒满足目标级配，实施难度高，且不符合现实瞬态填充过程的连续性规律。
发明内容
本发明所要解决的技术问题是提供一种基于连续均匀概率确定颗粒粒径的目标级配模型创建方法，在满足颗粒级配曲线的同时，能够大幅提高模型的最大填充率，同时保证瞬态填充过程的连续性。
本发明为了解决上述技术问题采用以下技术方案：本发明设计了一种基于连续均匀概率确定颗粒粒径的目标级配模型创建方法，用于实现对目标空间的颗粒投放，包括如下步骤：
步骤A.获得投放颗粒类型所对应的颗粒级配曲线，并根据该颗粒级配曲线中累计颗粒所占区域大小百分比
Figure BDA0002111336230000021
随粒径区间的变化关系，获得累计颗粒数量百分比
Figure BDA0002111336230000022
随粒径区间的变化关系，然后进入步骤B；其中j∈{1、…、J}，J表示投放颗粒类型所对应颗粒级配曲线中的总粒径区间数，
Figure BDA0002111336230000023
表示由第1个粒径区间至第j个粒径区间、累计颗粒所占区域大小百分比，
Figure BDA0002111336230000024
表示由第1个粒径区间至第j个粒径区间、累计颗粒数量百分比；
步骤B.针对目标空间构建坐标系，并获得坐标系中各坐标轴上最小值组成的坐标，作为Xmin，以及获得坐标系中各坐标轴上最大值组成的坐标，作为Xmax，然后进入步骤C；
步骤C.若针对目标空间的颗粒投放为静态投放，则进入步骤D，若针对目标空间的颗粒投放为瞬态投放，则然后初始化i＝1，并进入步骤M；
步骤D.根据目标空间的内部区域大小、以及目标空间的预设填充率，获得目标空间所需各填充颗粒所占区域大小之和V，然后进入步骤E；
步骤E.根据累计颗粒所占区域大小百分比
Figure BDA0002111336230000025
随粒径区间的变化关系、各粒径区间的上限Dj、以及目标空间所需各填充颗粒所占区域大小之和V，获得目标空间所需各填充颗粒的数量I，然后初始化i＝1，并进入步骤F；
步骤F.根据
Figure BDA0002111336230000026
获得针对目标空间进行第i次颗粒投放、颗粒所对应的颗粒连续均匀概率Pi，然后进入步骤G，其中，rand∈[0,1)；
步骤G.根据累计颗粒数量百分比
Figure BDA0002111336230000031
随粒径区间的变化关系、各粒径区间的上限Dj、以及针对目标空间进行第i次颗粒投放、颗粒所对应的颗粒连续均匀概率Pi，获得针对目标空间进行第i次颗粒投放、颗粒所对应的粒径di，然后进入步骤H；
步骤H.根据如下公式：
X＝Xmin+rand·(Xmax-Xmin)
获得更新位于目标空间内部的当前随机位置X，然后进入步骤J；
步骤J.针对目标空间进行第i次颗粒投放，将粒径为di的颗粒投至目标空间中当前随机位置X，并判断此次所投放颗粒是否与目标空间中已投放颗粒之间存在重叠，是则累计静态投放第i次颗粒尝试投放失败1次，并进入步骤K；否则针对i的值进行加1更新，并进入步骤L；
步骤K.判断累计静态投放第i次颗粒尝试投放的总次数是否达到预设失败阈值，是则判定针对目标空间的颗粒投放失败，方法结束；否则返回步骤H；
步骤L.判断i是否大于I，是则完成针对目标空间的静态投放；否则返回步骤F；
步骤M.根据Pi＝rand，获得针对目标空间进行第i次颗粒投放、颗粒所对应的颗粒连续均匀概率Pi，然后进入步骤N；
步骤N.根据累计颗粒数量百分比
Figure BDA0002111336230000032
随粒径区间的变化关系、各粒径区间的上限Dj、以及针对目标空间进行第i次颗粒投放、颗粒所对应的颗粒连续均匀概率Pi，获得针对目标空间进行第i次颗粒投放、颗粒所对应的粒径di，然后进入步骤O；
步骤O.根据如下公式：
X＝Xmin+rand·(Xmax-Xmin)
获得更新位于目标空间内部的当前随机位置X，然后进入步骤P；
步骤P.针对目标空间进行第i次颗粒投放，将粒径为di的颗粒投至目标空间中当前随机位置X，并判断此次所投放颗粒是否与目标空间中已投放颗粒之间存在重叠，是则累计瞬态投放第i次颗粒尝试投放失败1次，并进入步骤Q；否则进入步骤R；
步骤Q.判断累计瞬态投放第i次颗粒尝试投放的总次数是否达到预设失败阈值，是则判定针对目标空间的颗粒投放失败，方法结束；否则返回步骤O；
步骤R.根据针对目标空间已投放各颗粒的粒径，获得针对目标空间已投放各颗粒所占区域大小之和，作为针对目标空间已投放颗粒流量，并判断针对目标空间已投放颗粒流量是否达到预设流量上限，是则完成针对目标空间的瞬态投放；否则针对i的值进行加1更新，并返回步骤M。
作为本发明的一种优选技术方案：若所述目标空间为二维区域，则针对目标空间所投放颗粒为二维颗粒，则
Figure BDA0002111336230000041
表示由第1个粒径区间至第j个粒径区间、累计颗粒面积百分比；所述步骤D中，根据目标空间的内部面积、以及目标空间的预设填充率，获得目标空间所需各填充颗粒面积之和V；所述步骤R中，根据针对目标空间已投放各颗粒的粒径，获得针对目标空间已投放各颗粒面积之和，作为针对目标空间已投放颗粒流量；
若所述目标空间为三维区域，则针对目标空间所投放颗粒为三维颗粒，则
Figure BDA0002111336230000042
表示由第1个粒径区间至第j个粒径区间、累计颗粒体积百分比；所述步骤D中，根据目标空间的内部体积、以及目标空间的预设填充率，获得目标空间所需各填充颗粒体积之和V；所述步骤R中，根据针对目标空间已投放各颗粒的粒径，获得针对目标空间已投放各颗粒体积之和，作为针对目标空间已投放颗粒流量。
作为本发明的一种优选技术方案，若所述目标空间为二维区域：
所述步骤A中，根据投放颗粒类型所对应颗粒级配曲线中、累计颗粒面积百分比
Figure BDA0002111336230000043
随粒径区间的变化关系，按如下公式：
Figure BDA0002111336230000044
获得累计颗粒数量百分比
Figure BDA0002111336230000045
随粒径区间的变化关系；
若所述目标空间为三维区域：
所述步骤A中，根据投放颗粒类型所对应颗粒级配曲线中、累计颗粒体积百分比
Figure BDA0002111336230000046
随粒径区间的变化关系，按如下公式：
Figure BDA0002111336230000047
获得累计颗粒数量百分比
Figure BDA0002111336230000048
随粒径区间的变化关系。
作为本发明的一种优选技术方案，若所述目标空间为二维区域：
步骤E中，根据累计颗粒面积百分比
Figure BDA0002111336230000051
随粒径区间的变化关系、各粒径区间的上限Dj、以及目标空间所需各填充颗粒面积之和V，按如下公式：
Figure BDA0002111336230000052
获得目标空间所需各填充颗粒的数量I；
若所述目标空间为三维区域：
步骤E中，根据累计颗粒体积百分比
Figure BDA0002111336230000053
随粒径区间的变化关系、各粒径区间的上限Dj、以及目标空间所需各填充颗粒体积之和V，按如下公式：
Figure BDA0002111336230000054
获得目标空间所需各填充颗粒的数量I。
作为本发明的一种优选技术方案，所述步骤G和步骤N中：
根据累计颗粒数量百分比
Figure BDA0002111336230000055
随粒径区间的变化关系、各粒径区间的上限Dj、以及针对目标空间进行第i次颗粒投放、颗粒所对应的颗粒连续均匀概率Pi，按如下公式：
Figure BDA0002111336230000056
获得针对目标空间进行第i次颗粒投放、颗粒所对应的粒径di。
本发明所述一种基于连续均匀概率确定颗粒粒径的目标级配模型创建方法，采用以上技术方案与现有技术相比，具有以下技术效果：
本发明所设计基于连续均匀概率确定颗粒粒径的目标级配模型创建方法，引入连续均匀概率，依据累计颗粒数量百分比随粒径区间的变化关系，在全部粒径区间范围内选取颗粒粒径进行颗粒投放；与现行静态投放方法相比，能够大幅提高模型最大填充率；与现行瞬态填充方法相比，能够避免各级配区间颗粒在目标投放空间上、位置分布规律与投放顺序的相关性，确保各级配区间颗粒在目标投放空间上的位置、均服从目标分布规律，同时能够在整个投放过程中，使颗粒粒径不间断地满足目标级配分布，确保瞬态填充过程的连续性，因此，本发明适用于任意数值模拟技术领域目标级配数值模型的创建，且基于等效粒径能够满足任意颗粒形状，具有普遍适用性，并且本发明创建的数值模型更加符合现实中的客观规律。
附图说明
图1是本发明设计基于连续均匀概率确定颗粒粒径的目标级配模型创建方法的流程图；
图2a是静态投放中二维区域目标空间示意图；
图2b是静态投放中三维区域目标空间示意图；
图2c是瞬态投放示意图；
图3a是累计颗粒所占区域大小百分比随粒径区间的变化关系示意图；
图3b是累计颗粒数量百分比随粒径区间的变化关系示意图；
图4a是实施例静态投放中二维区域目标空间示意图；
图4b是实施例静态投放中三维区域目标空间示意图；
图4c是实施例瞬态投放示意图。
具体实施方式
下面结合说明书附图对本发明的具体实施方式作进一步详细的说明。
本发明设计了一种基于连续均匀概率确定颗粒粒径的目标级配模型创建方法，用于实现对目标空间的颗粒投放，实际应用当中，对于目标空间来说，可以划分为二维区域和三维区域两类，接下来针对这两类目标空间，分别进行方法的执行。
其一是二维区域目标空间，如图1和图2a所示，具体执行如下步骤。
步骤A.获得投放颗粒类型所对应的颗粒级配曲线，并如图3a所示，根据该颗粒级配曲线中累计颗粒所占面积大小百分比
Figure BDA0002111336230000061
随粒径区间的变化关系，按如下公式：
Figure BDA0002111336230000062
获得累计颗粒数量百分比
Figure BDA0002111336230000063
随粒径区间的变化关系，如图3b所示，然后进入步骤B；其中j∈{1、…、J}，J表示投放颗粒类型所对应颗粒级配曲线中的总粒径区间数，
Figure BDA0002111336230000064
表示由第1个粒径区间至第j个粒径区间、累计颗粒所占面积大小百分比，
Figure BDA0002111336230000065
表示由第1个粒径区间至第j个粒径区间、累计颗粒数量百分比，Dj表示第j个粒径区间的上限。
步骤B.针对目标空间构建坐标系，并获得坐标系中各坐标轴上最小值组成的坐标，作为Xmin，以及获得坐标系中各坐标轴上最大值组成的坐标，作为Xmax，然后进入步骤C。
步骤C.若针对目标空间的颗粒投放为静态投放，则进入步骤D，若针对目标空间的颗粒投放为瞬态投放，则然后初始化i＝1，如图2c所示，并进入步骤M。
步骤D.根据目标空间的内部面积大小、以及目标空间的预设填充率，获得目标空间所需各填充颗粒所占面积大小之和V，然后进入步骤E。
步骤E.根据累计颗粒所占面积大小百分比
Figure BDA0002111336230000071
随粒径区间的变化关系、各粒径区间的上限Dj、以及目标空间所需各填充颗粒所占面积大小之和V，按如下公式：
Figure BDA0002111336230000072
获得目标空间所需各填充颗粒的数量I，然后初始化i＝1，并进入步骤F。
步骤F.根据
Figure BDA0002111336230000073
获得针对目标空间进行第i次颗粒投放、颗粒所对应的颗粒连续均匀概率Pi，然后进入步骤G，其中，rand∈[0,1)。
步骤G.根据累计颗粒数量百分比
Figure BDA0002111336230000074
随粒径区间的变化关系、各粒径区间的上限Dj、以及针对目标空间进行第i次颗粒投放、颗粒所对应的颗粒连续均匀概率Pi，按如下公式：
Figure BDA0002111336230000075
获得针对目标空间进行第i次颗粒投放、颗粒所对应的粒径di，然后进入步骤H。
步骤H.根据如下公式：
X＝Xmin+rand·(Xmax-Xmin)
获得更新位于目标空间内部的当前随机位置X，然后进入步骤J。
步骤J.针对目标空间进行第i次颗粒投放，将粒径为di的颗粒投至目标空间中当前随机位置X，并判断此次所投放颗粒是否与目标空间中已投放颗粒之间存在重叠，是则累计静态投放第i次颗粒尝试投放失败1次，并进入步骤K；否则针对i的值进行加1更新，并进入步骤L。
步骤K.判断累计静态投放第i次颗粒尝试投放的总次数是否达到预设失败阈值，是则判定针对目标空间的颗粒投放失败，方法结束；否则返回步骤H。
步骤L.判断i是否大于I，是则完成针对目标空间的静态投放；否则返回步骤F。
步骤M.根据Pi＝rand，获得针对目标空间进行第i次颗粒投放、颗粒所对应的颗粒连续均匀概率Pi，然后进入步骤N。
步骤N.根据累计颗粒数量百分比
Figure BDA0002111336230000081
随粒径区间的变化关系、各粒径区间的上限Dj、以及针对目标空间进行第i次颗粒投放、颗粒所对应的颗粒连续均匀概率Pi，按如下公式：
Figure BDA0002111336230000082
获得针对目标空间进行第i次颗粒投放、颗粒所对应的粒径di，然后进入步骤O。
步骤O.根据如下公式：
X＝Xmin+rand·(Xmax-Xmin)
获得更新位于目标空间内部的当前随机位置X，然后进入步骤P。
步骤P.针对目标空间进行第i次颗粒投放，将粒径为di的颗粒投至目标空间中当前随机位置X，并判断此次所投放颗粒是否与目标空间中已投放颗粒之间存在重叠，是则累计瞬态投放第i次颗粒尝试投放失败1次，并进入步骤Q；否则进入步骤R。
步骤Q.判断累计瞬态投放第i次颗粒尝试投放的总次数是否达到预设失败阈值，是则判定针对目标空间的颗粒投放失败，方法结束；否则返回步骤O。
步骤R.根据针对目标空间已投放各颗粒的粒径，获得针对目标空间已投放各颗粒所占面积大小之和，作为针对目标空间已投放颗粒流量，并判断针对目标空间已投放颗粒流量是否达到预设流量上限，是则完成针对目标空间的瞬态投放；否则针对i的值进行加1更新，并返回步骤M。
上述针对二维区域目标空间的投放过程，在实际应用当中，如图4a所示，静态投放下，填充率达到83.98％，如图4c所示，瞬态投放，满足过程连续性。
其二是三维区域目标空间，如图1和图2b所示，具体执行如下步骤。
步骤A.获得投放颗粒类型所对应的颗粒级配曲线，并如图3a所示，根据该颗粒级配曲线中累计颗粒所占体积大小百分比
Figure BDA0002111336230000083
随粒径区间的变化关系，按如下公式：
Figure BDA0002111336230000084
获得累计颗粒数量百分比
Figure BDA0002111336230000091
随粒径区间的变化关系，如图3b所示，然后进入步骤B；其中j∈{1、…、J}，J表示投放颗粒类型所对应颗粒级配曲线中的总粒径区间数，
Figure BDA0002111336230000092
表示由第1个粒径区间至第j个粒径区间、累计颗粒所占体积大小百分比，
Figure BDA0002111336230000093
表示由第1个粒径区间至第j个粒径区间、累计颗粒数量百分比，Dj表示第j个粒径区间的上限。
步骤B.针对目标空间构建坐标系，并获得坐标系中各坐标轴上最小值组成的坐标，作为Xmin，以及获得坐标系中各坐标轴上最大值组成的坐标，作为Xmax，然后进入步骤C。
步骤C.若针对目标空间的颗粒投放为静态投放，则进入步骤D，若针对目标空间的颗粒投放为瞬态投放，则如图2c所示，然后初始化i＝1，并进入步骤M。
步骤D.根据目标空间的内部体积大小、以及目标空间的预设填充率，获得目标空间所需各填充颗粒所占体积大小之和V，然后进入步骤E。
步骤E.根据累计颗粒所占体积大小百分比
Figure BDA0002111336230000094
随粒径区间的变化关系、各粒径区间的上限Dj、以及目标空间所需各填充颗粒所占体积大小之和V，按如下公式：
Figure BDA0002111336230000095
获得目标空间所需各填充颗粒的数量I，然后初始化i＝1，并进入步骤F。
步骤F.根据
Figure BDA0002111336230000096
获得针对目标空间进行第i次颗粒投放、颗粒所对应的颗粒连续均匀概率Pi，然后进入步骤G，其中，rand∈[0,1)。
步骤G.根据累计颗粒数量百分比
Figure BDA0002111336230000097
随粒径区间的变化关系、各粒径区间的上限Dj、以及针对目标空间进行第i次颗粒投放、颗粒所对应的颗粒连续均匀概率Pi，按如下公式：
Figure BDA0002111336230000098
获得针对目标空间进行第i次颗粒投放、颗粒所对应的粒径di，然后进入步骤H。
步骤H.根据如下公式：
X＝Xmin+rand·(Xmax-Xmin)
获得更新位于目标空间内部的当前随机位置X，然后进入步骤J。
步骤J.针对目标空间进行第i次颗粒投放，将粒径为di的颗粒投至目标空间中当前随机位置X，并判断此次所投放颗粒是否与目标空间中已投放颗粒之间存在重叠，是则累计静态投放第i次颗粒尝试投放失败1次，并进入步骤K；否则针对i的值进行加1更新，并进入步骤L。
步骤K.判断累计静态投放第i次颗粒尝试投放的总次数是否达到预设失败阈值，是则判定针对目标空间的颗粒投放失败，方法结束；否则返回步骤H。
步骤L.判断i是否大于I，是则完成针对目标空间的静态投放；否则返回步骤F。
步骤M.根据Pi＝rand，获得针对目标空间进行第i次颗粒投放、颗粒所对应的颗粒连续均匀概率Pi，然后进入步骤N。
步骤N.根据累计颗粒数量百分比
Figure BDA0002111336230000101
随粒径区间的变化关系、各粒径区间的上限Dj、以及针对目标空间进行第i次颗粒投放、颗粒所对应的颗粒连续均匀概率Pi，按如下公式：
Figure BDA0002111336230000102
获得针对目标空间进行第i次颗粒投放、颗粒所对应的粒径di，然后进入步骤O。
步骤O.根据如下公式：
X＝Xmin+rand·(Xmax-Xmin)
获得更新位于目标空间内部的当前随机位置X，然后进入步骤P。
步骤P.针对目标空间进行第i次颗粒投放，将粒径为di的颗粒投至目标空间中当前随机位置X，并判断此次所投放颗粒是否与目标空间中已投放颗粒之间存在重叠，是则累计瞬态投放第i次颗粒尝试投放失败1次，并进入步骤Q；否则进入步骤R。
步骤Q.判断累计瞬态投放第i次颗粒尝试投放的总次数是否达到预设失败阈值，是则判定针对目标空间的颗粒投放失败，方法结束；否则返回步骤O。
步骤R.根据针对目标空间已投放各颗粒的粒径，获得针对目标空间已投放各颗粒所占体积大小之和，作为针对目标空间已投放颗粒流量，并判断针对目标空间已投放颗粒流量是否达到预设流量上限，是则完成针对目标空间的瞬态投放；否则针对i的值进行加1更新，并返回步骤M。
上述针对三维区域目标空间的投放过程，在实际应用当中，如图4b所示，静态投放下，填充率达到60.90％，如图4c所示，瞬态投放，满足过程连续性。
上述技术方案所设计基于连续均匀概率确定颗粒粒径的目标级配模型创建方法，引入连续均匀概率，依据累计颗粒数量百分比随粒径区间的变化关系，在全部粒径区间范围内选取颗粒粒径进行颗粒投放；与现行静态投放方法相比，能够大幅提高模型最大填充率；与现行瞬态填充方法相比，能够避免各级配区间颗粒在目标投放空间上、位置分布规律与投放顺序的相关性，确保各级配区间颗粒在目标投放空间上的位置、均服从目标分布规律，同时能够在整个投放过程中，使颗粒粒径不间断地满足目标级配分布，确保瞬态填充过程的连续性，因此，本发明适用于任意数值模拟技术领域目标级配数值模型的创建，且基于等效粒径能够满足任意颗粒形状，具有普遍适用性，并且本发明创建的数值模型更加符合现实中的客观规律。
上面结合附图对本发明的实施方式作了详细说明，但是本发明并不限于上述实施方式，在本领域普通技术人员所具备的知识范围内，还可以在不脱离本发明宗旨的前提下做出各种变化。

Claims (4)

1.一种基于连续均匀概率确定颗粒粒径的目标级配模型创建方法，用于实现对目标空间的颗粒投放，其特征在于，包括如下步骤：
步骤A.获得投放颗粒类型所对应的颗粒级配曲线，并根据该颗粒级配曲线中累计颗粒所占区域大小百分比
Figure FDA0003877619590000011
随粒径区间的变化关系，获得累计颗粒数量百分比
Figure FDA0003877619590000012
随粒径区间的变化关系，然后进入步骤B；其中j∈{1、…、J}，J表示投放颗粒类型所对应颗粒级配曲线中的总粒径区间数，
Figure FDA0003877619590000013
表示由第1个粒径区间至第j个粒径区间、累计颗粒所占区域大小百分比，
Figure FDA0003877619590000014
表示由第1个粒径区间至第j个粒径区间、累计颗粒数量百分比；
步骤B.针对目标空间构建坐标系，并获得坐标系中各坐标轴上最小值组成的坐标，作为Xmin，以及获得坐标系中各坐标轴上最大值组成的坐标，作为Xmax，然后进入步骤C；
步骤C.若针对目标空间的颗粒投放为静态投放，则进入步骤D，若针对目标空间的颗粒投放为瞬态投放，则初始化i＝1，并进入步骤M；
步骤D.根据目标空间的内部区域大小、以及目标空间的预设填充率，获得目标空间所需各填充颗粒所占区域大小之和V，然后进入步骤E；
步骤E.根据累计颗粒所占区域大小百分比
Figure FDA0003877619590000015
随粒径区间的变化关系、各粒径区间的上限Dj、以及目标空间所需各填充颗粒所占区域大小之和V，获得目标空间所需各填充颗粒的数量I，然后初始化i＝1，并进入步骤F；
步骤F.根据
Figure FDA0003877619590000016
获得针对目标空间进行第i次颗粒投放、颗粒所对应的颗粒连续均匀概率Pi，然后进入步骤G，其中，rand∈[0,1)；
步骤G.根据累计颗粒数量百分比
Figure FDA0003877619590000017
随粒径区间的变化关系、各粒径区间的上限Dj、以及针对目标空间进行第i次颗粒投放、颗粒所对应的颗粒连续均匀概率Pi，按如下公式：
Figure FDA0003877619590000018
获得针对目标空间进行第i次颗粒投放、颗粒所对应的粒径di，然后进入步骤H；
步骤H.根据如下公式：
X＝Xmin+rand·(Xmax-Xmin)
获得更新位于目标空间内部的当前随机位置X，然后进入步骤J；
步骤J.针对目标空间进行第i次颗粒投放，将粒径为di的颗粒投放至目标空间中当前随机位置X，并判断此次所投放颗粒是否与目标空间中已投放颗粒之间存在重叠，是则累计静态投放第i次颗粒尝试投放失败1次，并进入步骤K；否则针对i的值进行加1更新，并进入步骤L；
步骤K.判断累计静态投放第i次颗粒尝试投放的总次数是否达到预设失败阈值，是则判定针对目标空间的颗粒投放失败，方法结束；否则返回步骤H；
步骤L.判断i是否大于I，是则完成针对目标空间的静态投放；否则返回步骤F；
步骤M.根据Pi＝rand，获得针对目标空间进行第i次颗粒投放、颗粒所对应的颗粒连续均匀概率Pi，然后进入步骤N；
步骤N.根据累计颗粒数量百分比
Figure FDA0003877619590000021
随粒径区间的变化关系、各粒径区间的上限Dj、以及针对目标空间进行第i次颗粒投放、颗粒所对应的颗粒连续均匀概率Pi，按如下公式：
Figure FDA0003877619590000022
获得针对目标空间进行第i次颗粒投放、颗粒所对应的粒径di，然后进入步骤O；
步骤O.根据如下公式：
X＝Xmin+rand·(Xmax-Xmin)
获得更新位于目标空间内部的当前随机位置X，然后进入步骤P；
步骤P.针对目标空间进行第i次颗粒投放，将粒径为di的颗粒投至目标空间中当前随机位置X，并判断此次所投放颗粒是否与目标空间中已投放颗粒之间存在重叠，是则累计瞬态投放第i次颗粒尝试投放失败1次，并进入步骤Q；否则进入步骤R；
步骤Q.判断累计瞬态投放第i次颗粒尝试投放的总次数是否达到预设失败阈值，是则判定针对目标空间的颗粒投放失败，方法结束；否则返回步骤O；
步骤R.根据针对目标空间已投放各颗粒的粒径，获得针对目标空间已投放各颗粒所占区域大小之和，作为针对目标空间已投放颗粒流量，并判断针对目标空间已投放颗粒流量是否达到预设流量上限，是则完成针对目标空间的瞬态投放；否则针对i的值进行加1更新，并返回步骤M。
2.根据权利要求1所述一种基于连续均匀概率确定颗粒粒径的目标级配模型创建方法，其特征在于：若所述目标空间为二维区域，则针对目标空间所投放颗粒为二维颗粒，则
Figure FDA0003877619590000031
表示由第1个粒径区间至第j个粒径区间、累计颗粒面积百分比；所述步骤D中，根据目标空间的内部面积、以及目标空间的预设填充率，获得目标空间所需各填充颗粒面积之和V；所述步骤R中，根据针对目标空间已投放各颗粒的粒径，获得针对目标空间已投放各颗粒面积之和，作为针对目标空间已投放颗粒流量；
若所述目标空间为三维区域，则针对目标空间所投放颗粒为三维颗粒，则
Figure FDA0003877619590000032
表示由第1个粒径区间至第j个粒径区间、累计颗粒体积百分比；所述步骤D中，根据目标空间的内部体积、以及目标空间的预设填充率，获得目标空间所需各填充颗粒体积之和V；所述步骤R中，根据针对目标空间已投放各颗粒的粒径，获得针对目标空间已投放各颗粒体积之和，作为针对目标空间已投放颗粒流量。
3.根据权利要求2所述一种基于连续均匀概率确定颗粒粒径的目标级配模型创建方法，其特征在于，若所述目标空间为二维区域：
所述步骤A中，根据投放颗粒类型所对应颗粒级配曲线中、累计颗粒面积百分比
Figure FDA0003877619590000033
随粒径区间的变化关系，按如下公式：
Figure FDA0003877619590000034
获得累计颗粒数量百分比
Figure FDA0003877619590000035
随粒径区间的变化关系；
若所述目标空间为三维区域：
所述步骤A中，根据投放颗粒类型所对应颗粒级配曲线中、累计颗粒体积百分比
Figure FDA0003877619590000036
随粒径区间的变化关系，按如下公式：
Figure FDA0003877619590000037
获得累计颗粒数量百分比
Figure FDA0003877619590000038
随粒径区间的变化关系。
4.根据权利要求2所述一种基于连续均匀概率确定颗粒粒径的目标级配模型创建方法，其特征在于，若所述目标空间为二维区域：
步骤E中，根据累计颗粒面积百分比
Figure FDA0003877619590000039
随粒径区间的变化关系、各粒径区间的上限Dj、以及目标空间所需各填充颗粒面积之和V，按如下公式：
Figure FDA0003877619590000041
获得目标空间所需各填充颗粒的数量I；
若所述目标空间为三维区域：
步骤E中，根据累计颗粒体积百分比
Figure FDA0003877619590000042
随粒径区间的变化关系、各粒径区间的上限Dj、以及目标空间所需各填充颗粒体积之和V，按如下公式：
Figure FDA0003877619590000043
获得目标空间所需各填充颗粒的数量I。
CN201910572878.5A 2019-06-28 2019-06-28 一种基于连续均匀概率确定颗粒粒径的目标级配模型创建方法 Active CN110245457B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910572878.5A CN110245457B (zh) 2019-06-28 2019-06-28 一种基于连续均匀概率确定颗粒粒径的目标级配模型创建方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910572878.5A CN110245457B (zh) 2019-06-28 2019-06-28 一种基于连续均匀概率确定颗粒粒径的目标级配模型创建方法

Publications (2)

Publication Number Publication Date
CN110245457A CN110245457A (zh) 2019-09-17
CN110245457B true CN110245457B (zh) 2023-01-03

Family

ID=67890116

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910572878.5A Active CN110245457B (zh) 2019-06-28 2019-06-28 一种基于连续均匀概率确定颗粒粒径的目标级配模型创建方法

Country Status (1)

Country Link
CN (1) CN110245457B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111324959B (zh) * 2020-02-25 2022-07-26 山东大学 模拟地下工程开挖的近场动力学物质点休眠方法及系统
CN111474089A (zh) * 2020-04-23 2020-07-31 佛山市优一家电制造有限公司 水汽释放设备的水分子粒径测试方法和装置
CN112131633B (zh) * 2020-09-04 2023-01-13 山东大学 一种基于粗粒化计算理论的流固耦合仿真方法及系统
CN112115608B (zh) * 2020-09-16 2024-03-15 中国地质大学（北京） 复合颗粒粒径计算和级配调整配置方法
US11797908B2 (en) * 2020-09-16 2023-10-24 Cheetah Technologies Inc. Fill modeling for hybrid last-mile delivery
CN112146957B (zh) * 2020-09-25 2022-04-15 东北石油大学 一种基于数字岩心实现人造岩心定量制作方法
CN112650305B (zh) * 2020-12-22 2022-07-15 中国航空工业集团公司沈阳飞机设计研究所 一种复杂环境下的攻击批次划分方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007286514A (ja) * 2006-04-19 2007-11-01 Doshisha 粒子挙動シミュレーション装置、粒子挙動シミュレーション方法、及びコンピュータプログラム
CN104318039A (zh) * 2014-11-18 2015-01-28 张光亮 硬质合金自适应参数化建模的方法
CN107704653A (zh) * 2017-08-31 2018-02-16 河海大学 一种水泥基材料三维颗粒分组均匀化投放方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007286514A (ja) * 2006-04-19 2007-11-01 Doshisha 粒子挙動シミュレーション装置、粒子挙動シミュレーション方法、及びコンピュータプログラム
CN104318039A (zh) * 2014-11-18 2015-01-28 张光亮 硬质合金自适应参数化建模的方法
CN107704653A (zh) * 2017-08-31 2018-02-16 河海大学 一种水泥基材料三维颗粒分组均匀化投放方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
DEM study of the shear behavior and formation of shear band in biaxial test;Weichen Sun et al.;《Advanced Powder Technology》;20201231;第1431-1440页 *
Effect of sample size on the response of DEM samples with arealistic grading;X. Huang et al.;《Particuology》;20141231;第107-115页 *

Also Published As

Publication number Publication date
CN110245457A (zh) 2019-09-17

Similar Documents

Publication Publication Date Title
CN110245457B (zh) 一种基于连续均匀概率确定颗粒粒径的目标级配模型创建方法
CN109376438B (zh) 一种钢筋下料优化方法、装置及存储设备
CN108846186B (zh) 一种考虑孔洞缺陷形状及其分布随机性的脆性材料离散元建模方法
CN107193212A (zh) 基于新型灰狼优化算法的航空发动机非线性预测控制方法
CN103971289B (zh) 矿山放矿数据的处理方法及装置
CN104081400A (zh) 使用多事务技术对面向对象环境中的对象执行几何变换的方法和装置
CN103295104A (zh) 一种限额设计管理方法及用于成本核算的项目投资服务平台
CN112818469B (zh) 固体火箭发动机映射设计方法、装置和设备
CN105136623A (zh) 基于势能变化定量表征颗粒下落后堆积偏析状态的方法
Hou et al. A metaheuristic algorithm for routing school buses with mixed load
CN104827575B (zh) 一种混凝土骨料二维随机投放的双向游走法
Tavana Amlashi et al. Developing three hybrid machine learning algorithms for predicting the mechanical properties of plastic concrete samples with different geometries
CN109325288B (zh) 一种基于不确定性优化的固体运载器总体参数确定方法及系统
CN104199984B (zh) 一种实现功能库的动态引用的方法和装置
JP6806034B2 (ja) 非線形計画問題の求解方法、配合計画の作成方法、およびこれらのコンピュータプログラム、記録媒体、装置、ならびに配合操業方法
CN105022896A (zh) 一种基于动态编号的apdl建模方法及装置
CN103942195B (zh) 一种数据处理系统以及数据处理方法
CN102799939A (zh) 一种生物质炉燃烧优化的模型更新方法
CN113092324A (zh) 一种动态随机序列吸附方法、系统、设备及存储介质
CN114241030A (zh) 自动批量填充布置方法、装置、电子设备及存储介质
CN107918796A (zh) 一种面向非线性规划问题求解的多策略混合蛙跳算法
CN106528906A (zh) 一种基于日常岩粉数据的自动化配矿方法
Rogala et al. Modelling the set of earthworks machinery with the use of computer simulation
Ekanayaka et al. Implementation of a surrogate model for a novel path‐based finite element simulation for additive manufacturing processes in construction
CN115130320B (zh) 一种适用于纤维增强水泥基复合材料的空间纤维模拟方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant