CN106644868A - 一种二维非凸形随机骨料周围界面浓度的测定方法 - Google Patents
一种二维非凸形随机骨料周围界面浓度的测定方法 Download PDFInfo
- Publication number
- CN106644868A CN106644868A CN201710068971.3A CN201710068971A CN106644868A CN 106644868 A CN106644868 A CN 106644868A CN 201710068971 A CN201710068971 A CN 201710068971A CN 106644868 A CN106644868 A CN 106644868A
- Authority
- CN
- China
- Prior art keywords
- aggregate
- convex
- random
- aggregates
- dimentional
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 15
- 230000002093 peripheral effect Effects 0.000 title abstract 4
- 239000002245 particle Substances 0.000 claims abstract description 64
- 238000005070 sampling Methods 0.000 claims abstract description 5
- 238000003556 assay Methods 0.000 claims description 19
- 238000004088 simulation Methods 0.000 claims description 10
- 239000011159 matrix material Substances 0.000 claims description 9
- 239000000203 mixture Substances 0.000 claims description 4
- 238000009795 derivation Methods 0.000 claims description 3
- 239000008187 granular material Substances 0.000 claims description 3
- 230000032696 parturition Effects 0.000 claims 1
- 230000010339 dilation Effects 0.000 abstract 1
- 238000011160 research Methods 0.000 description 11
- 239000000463 material Substances 0.000 description 5
- 238000002474 experimental method Methods 0.000 description 3
- 230000010354 integration Effects 0.000 description 3
- 239000002002 slurry Substances 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005315 distribution function Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000007619 statistical method Methods 0.000 description 2
- 241000238367 Mya arenaria Species 0.000 description 1
- 241001562042 Physa Species 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
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/06—Investigating concentration of particle suspensions
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N13/00—Investigating surface or boundary effects, e.g. wetting power; Investigating diffusion effects; Analysing materials by determining surface, boundary, or diffusion effects
Landscapes
- Chemical & Material Sciences (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Dispersion Chemistry (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种二维非凸形随机骨料周围界面浓度的测定方法,包括如下步骤:1、基于富勒公式推导的瓦拉文公式获取在富勒级配下各粒径范围的二维骨料数量密度;2、基于简单椭圆曲线通过单峰值的伸缩因子变形得到二维非凸形骨料模型;3、采用伸缩矩阵获取非凸骨料周围界面的细观结构信息;4、根据步骤1计算出的二维非凸形骨料颗粒的数量密度和步骤2给出的骨料模型生成骨料颗粒、进行骨料重叠判断并进行随机堆积;5、根据步骤4得到的结果采用Monte Carlo随机点采样法和统计原理测定界面浓度。该方法提供了一种测量非凸形骨料周围界面浓度的方法,使得界面浓度含量的测定方法更具普适性和代表性。
Description
技术领域
本发明涉及一种界面浓度的测定方法,具体涉及一种二维非凸形随机骨料周围界面浓度的测定方法,属于混凝土细观力学理论和数值试验技术领域。
背景技术
大量的试验研究已经证实界面广泛地存在于不同类型混凝土的骨料—浆体之间,由于具有高孔隙率、低刚度的物理特性,界面在材料中所占的含量严重制约着混凝土结构的承载能力和抗侵蚀能力。但是,目前仍然无法通过常规的微细观实验技术直接测量界面浓度。当前,国际上主要通过理论研究和数值仿真两种途径来获取界面浓度。
在理论研究方面,最具代表性的工作是美国国家标准技术研究所(NIST)将hardcore/soft shell模型和最邻近表面分布函数理论引入混凝土材料中研究球形骨料—浆体间界面浓度,然而,该理论模型仅仅适用于最简单的球形随机骨料情况。为了更深入地分析界面浓度,国内外学者也试图通过改进球形骨料的最邻近表面分布函数获取一些复杂凸形骨料—浆体间界面浓度理论模型,如:椭球形随机骨料模型、凸面体随机骨料模型以及柱状纤维模型。尽管这些理论模型都已经证实骨料的几何特性对界面浓度具有重要的影响。但是根据模型所得到的结果却大相径庭,很难应用于实际的工程结构中。
随着现代计算机硬件和软件技术的飞速发展,数值仿真方法相比较于理论研究无疑更加准确、真实地模拟和分析实际混凝土材料的细观结构和界面浓度。在以往的数值模拟研究中,通常先构建混凝土细观结构模型,并在此基础上数值模拟界面浓度。因此,数值仿真技术主要由两部分组成:第一部分是输入含有界面的混凝土细观结构信息,这要求研究者们首先必须建立满足混凝土自身特征的细观结构模型。现有的混凝土细观结构模型包括从早期最简单的球形骨料模型到最近的复杂凸体骨料模型。实际混凝土的骨料成分并非是完美的凸形体,而是具有凹凸相间的非凸体组成;然而,关于含有非凸形骨料的混凝土细观结构模型的研究鲜见报道。第二部分是数值统计界面浓度,采用的统计方法包括一点概率函数、线样条、点采样、方差分析法等。现有的数值统计结果主要针对凸形骨料周围的界面浓度,但是不同研究机构所得到的数值结果却相差很大,例如,NIST揭示的椭球形骨料周围界面浓度与最近国内学者的数值研究结果正好相反。就其原因,在于两者所建立的材料细观结构模型上的差异。但是,非凸形骨料周围界面浓度的统计分析目前仍然是一片空白。因此,如何精确高效地测定非凸形骨料周围的界面浓度成为混凝土细观力学研究的热点及难点问题。
发明内容
发明目的:针对现有技术中存在的问题,本发明公开了一种二维非凸形随机骨料周围界面浓度的测定方法,克服了现有技术中只能针对凸形骨料研究的技术约束以及数值模拟所产生的误差的问题,使得界面浓度含量的测定方法更具普适性和代表性。
技术方案:本发明采用的技术方案如下:
一种二维非凸形随机骨料周围界面浓度的测定方法,包括如下步骤:
(1)基于富勒公式推导的瓦拉文公式获取在富勒级配下各粒径范围的二维骨料数量密度:
其中Pc(D<D0)为实际骨料粒径D小于D0的骨料密度,Dmax为最大骨料颗粒直径;
(2)基于简单椭圆曲线通过单峰值的伸缩因子变形得到二维非凸形骨料模型,所述二维非凸形骨料模型的椭圆基础参数方程为:
其中u∈(0,2π),α为椭圆上点的旋转角度,α∈(0,π),a为椭圆的半长轴,b为椭圆的半短轴;
所述二维非凸形骨料模型的伸缩因子函数如下:
其中函数g(s)的表达式如下:
其中r1为峰值区域边界,r2为支区域边界;
(3)采用伸缩矩阵获取非凸形骨料周围界面的细观结构信息;
(4)根据步骤1计算出的二维非凸形骨料颗粒的数量密度和步骤2给出的骨料模型生成骨料颗粒、进行骨料重叠判断并进行随机堆积;
(5)根据步骤4得到的结果采用Monte Carlo随机点采样法和统计原理测定界面浓度。
具体地,所述步骤(3)包括如下步骤:
(31)将步骤(2)建立的二维非凸形骨料模型的椭圆基础参数方程中的参数u在(0,2π)区间内拆分为n份,子区间中心为a0i,其中i=1,2,...,n;
(32)对第i个子区间随机生成两个轴向的伸缩因子β1i、β2i,构成伸缩矩阵Di:
Di=diag(β1i,β2i),
以子区间中心a0i为变形峰值向内收缩,得到骨料轮廓线散点P(tg)和界面的外边界轮廓线散点P(ts):
其中椭圆中心E为单位矩阵,γ为二维非凸形骨料周围界面层厚度;
骨料轮廓线所围的面积S(tg)和界面的外边界轮廓线所围的面积S(ts)分别为:
具体地,所述步骤(4)中生成骨料颗粒包括如下步骤:
(41)初始化骨料颗粒数量m=0和骨料颗粒的总面积s=0;
(42)根据非凸形骨料的尺寸大小按照步骤(1)计算二维非凸形骨料的数量密度Pc,从而得到骨料的数量N=Pc·Sc,其中Sc为模拟的样本面积;在粗骨料最大粒径下随机生成半长轴aj和半短轴bj;
(43)在模拟区域内随机生成一点Oj(x0j,y0j)作为生成二维非凸形骨料模型的椭圆基础的椭圆中心;
(44)判断点Oj是否在所有已生成的颗粒的覆盖椭圆的外部;如果否,转步骤(43)重新生成椭圆中心;若是,计算Oj与所有已生成骨料边界的距离,判断其中最小的距离d是否大于半长轴aj;如果d<aj,表示骨料颗粒有重叠,转步骤(43)重新生成椭圆中心;如果d≥aj,以Oj(x0j,y0j)为中心、aj为半长轴、bj为半短轴生成椭圆,根据步骤(2)对所述椭圆进行随机变形得到骨料颗粒Qj;判断Qj是否与已经生成的骨料颗粒重叠,如果重叠,则跳转至步骤(43)重新生成非凸形骨料颗粒的空间位置;根据步骤(3)计算骨料颗粒Qj的面积S0(tg),s=s+S0(tg),m=m+1;
(45)判断是否满足迭代收敛条件:m>N或s>S,如果满足,结束;否则转步骤(43);其中S为骨料颗粒最大总面积。
优选地,所述步骤(44)中判断骨料颗粒Qj是否与已生成的骨料颗粒Qk重叠,包括如下步骤:
(441)A为非凸骨料颗粒Qj边界的组成点,Oj(x0j,y0j)为Qj中心点;令θ为OjA与椭圆基础长轴的夹角;OjA与已生成的骨料颗粒Qk边界相交于点B;
(442)初始时令θ=0,计算Oj(x0j,y0j)与A、B之间的距离|OjA|和|OjB|,如果|OjA|>|OjB|,则Qj与Qk重叠;否则,不重叠;
(443)以步长△θ修改θ的值:θ=θ+△θ,0≤θ≤360°,再次计算Oj(x0j,y0j)与A、B之间的距离|OjA|和|OjB|,如果|OjA|>|OjB|,则Qj与Qk重叠;否则,不重叠。
具体地,所述步骤(5)包括:在模拟空间内大量投放随机点,按照随机点所处的空间位置将随机点分为骨料边界内随机点、骨料边界外随机点、界面最外层轮廓线内随机点和界面最外层轮廓线外随机点;
统计每种随机点的个数,得到界面浓度η:
其中Nsum为随机点总数,Nin为界面最外层轮廓线内随机点总数,Nout为骨料边界外随机点总数。
优选地,所述步骤(31)中参数u在(0,2π)区间内进行n等分,分为n个子区间。
优选地,所述步骤(32)中伸缩因子β1i、β2i在[-0.16,-0.1]区间内取值。
优选地,二维非凸形骨料周围界面层厚度γ为0.003毫米。
优选地,所述步骤(4)中按照非凸形骨料颗粒的粒径从大到小排序,使越大的粒子,越先进行堆积。
优选地,所述步骤(443)中步长△θ为0.01°。
有益效果:与现有技术相比,本发明公开的二维非凸形随机骨料周围界面浓度的测定方法具有以下优点:1、本发明公开的方法能够精确测定二维非凸形随机骨料周围的界面浓度;2、本发明建立了一种概念清晰、理论简单、操作便利、适用范围广的二维非凸形随机骨料周围界面浓度的数值框架,对推广和发展复合介质界面微结构关于宏观性能的作用机理分析具有十分重要的理论和现实意义;3、本发明克服了现有技术中凸形骨料研究的技术约束,使得界面浓度的测定方法更具普适性和代表性;4、本发明公开的方法更加科学、清晰地反映了二维非凸形随机骨料的几何特征。
附图说明
图1为基于椭圆曲线伸缩构造的二维非凸形骨料示意图;
图2是构造的二维非凸形骨料及其界面拓扑结构;
图3为二维非凸形随机骨料周围界面浓度测定的流程图;
图4为判断二维非凸形骨料之间重叠的示意图;
图5为二维非凸形骨料重叠判断实验图;
图6是骨料面积分数为40%时的混凝土细观结构模型;
图7是在不同界面厚度下二维非凸形骨料含量对界面浓度影响的数值结果。
具体实施方式
下面结合附图和具体实施方式,进一步阐明本发明。
一种二维非凸形随机骨料周围界面浓度的测定方法,包括如下步骤:
步骤1、利用基于富勒公式推导的瓦拉文公式和相关物理量确定骨料颗粒的数量、尺寸、分布的总面积,将应用于三维的富勒级配曲线转化为混凝土试件截面任一点具有骨料直径D<D0的内截圆出现的概率,从而得到在富勒级配下各个粒径范围的二维骨料数量密度Pc。具体转化步骤如下:
首先假定骨料为球形颗粒,利用骨料在空间上等概率分布和任一大小圆形切面无概率占优性,按照富勒曲线确定骨料的三维级配曲线,通过直径D筛孔的骨料的重量百分比Dmax代表最大骨料颗粒直径,认为在空间内任指定一点位于小于半径为Dx的骨料上的概率密度表示为:
其中Pk为骨料体积与混凝土总体积之比,一般情况下Pk=0.75。
根据空间几何关系得:
即所谓的瓦拉文公式。其中Pc(D<D0)为实际骨料粒径D小于D0的骨料密度。
步骤2、基于简单椭圆曲线通过单峰值的伸缩因子变形得到二维非凸形骨料模型,所述二维非凸形骨料模型的椭圆基础参数方程为:
其中u∈(0,2π),α为椭圆上点的旋转角度,α∈(0,π),a为椭圆的半长轴,b为椭圆的半短轴,如图1所示;
引入伸缩因子对区间[a0-r2,a0+r2]进行伸缩,其中a0为伸缩中心,r2为支区域边界,令r1为峰值区域边界,且有0≤r1<r2。令t=t(x)=(x-a0)2表示任一点到伸缩中心的距离,得到二维非凸形骨料模型的伸缩因子函数如下:
其中函数g(s)的表达式如下:
考虑到计算时的普适性,采用查高斯积分表的方法来计算伸缩因子的数值,本发明采用的是8个节点的高斯积分表,如表1所示。经过试验,利用积分区间内8个积分节点的高斯积分生成的伸缩因子可以更有效的表现出非凸形骨料的形貌,表1给出的是单位区间内的积分节点和各积分节点权重,即概率分布。
表1
节点 | 权重 |
-0.9602898564975360 | 0.1012285362903760 |
-0.7966664774136260 | 0.2223810344533740 |
-0.5255324099163290 | 0.3137066458778870 |
-0.1834346424956490 | 0.3626837833783620 |
0.1834346424956490 | 0.3626837833783620 |
0.5255324099163290 | 0.3137066458778870 |
0.7966664774136260 | 0.2223810344533740 |
0.9602898564975360 | 0.1012285362903760 |
步骤3、采用伸缩矩阵获取非凸形骨料周围界面的细观结构信息,具体步骤如下:
(31)将步骤2建立的二维非凸形骨料模型的椭圆基础参数方程中的参数u在(0,2π)区间内拆分为n份,子区间中心为a0i,其中i=1,2,...,n;本发明中参数u在(0,2π)区间内进行n等分,分为n个子区间。
(32)对第i个子区间随机生成两个轴向的伸缩因子β1i、β2i,构成伸缩矩阵Di:
Di=diag(β1i,β2i),
伸缩因子β1i、β2i在[-0.16,-0.1]区间内取值,使得生成的非凸形颗粒形貌更加接近于真实骨料颗粒形态。
以子区间中心a0i为变形峰值向内收缩,为保证只有一个最大收缩点,一般取r1为0,取r2为子区域中随机点。得到骨料轮廓线散点P(tg)和界面的外边界轮廓线散点P(ts):
其中椭圆中心E为单位矩阵,γ为二维非凸形骨料周围界面层厚度,本实施例中为0.003毫米;骨料边界如图2所示黑色线条为骨料轮廓线,灰色部分为厚度是γ的骨料周围界面。
轮廓线所围的面积S(tg)和界面的外边界轮廓线所围的面积S(ts)分别为:
步骤4、根据步骤1计算出的二维非凸形骨料颗粒的数量密度和步骤2给出的骨料模型生成骨料颗粒、进行骨料重叠判断并进行随机堆积;
生成骨料颗粒具体步骤如下:
(41)初始化骨料颗粒数量m=0和骨料颗粒的总面积s=0;
(42)根据非凸形骨料的尺寸大小按照步骤(1)计算二维非凸形骨料的数量密度Pc,从而得到骨料的数量N=Pc·Sc,其中Sc为模拟的样本面积;在粗骨料最大粒径下随机生成半长轴aj和半短轴bj;
(43)在模拟区域内随机生成一点Oj(x0j,y0j)作为生成二维非凸形骨料模型的椭圆基础的椭圆中心;
(44)判断点Oj是否在所有已生成的颗粒的覆盖椭圆的外部;如果否,转步骤(43)重新生成椭圆中心;若是,计算Oj与所有已生成骨料边界的距离,判断其中最小的距离d是否大于半长轴aj;如果d<aj,表示骨料颗粒有重叠,转步骤(43)重新生成椭圆中心;如果d≥aj,以Oj(x0j,y0j)为中心、aj为半长轴、bj为半短轴生成椭圆,根据步骤(2)对所述椭圆进行随机变形得到骨料颗粒Qj;判断Qj是否与已经生成的骨料颗粒重叠,如果重叠,则跳转至步骤(43)重新生成非凸形骨料颗粒的空间位置;根据步骤(3)计算骨料颗粒Qj的面积S0(tg),s=s+S0(tg),m=m+1;
(45)判断是否满足迭代收敛条件:m>N或s>S,如果满足,结束;否则转步骤(43);其中S为骨料颗粒最大总面积。
判断骨料颗粒Qj是否与已生成的骨料颗粒Qk重叠,包括如下步骤:
(441)A为非凸形骨料颗粒Qj边界的组成点,Oj(x0j,y0j)为Qj中心点;令θ为OjA与椭圆基础长轴的夹角;OjA与已生成的骨料颗粒Qk边界相交于点B;
(442)初始时令θ=0,计算Oj(x0j,y0j)与A、B之间的距离|OjA|和|OjB|,如果|OjA|>|OjB|,则Qj与Qk重叠;否则,不重叠;
(443)以步长△θ修改θ的值:θ=θ+△θ,0≤θ≤360°,再次计算Oj(x0j,y0j)与A、B之间的距离|OjA|和|OjB|,如图3所示,如果|OjA|>|OjB|,则Qj与Qk重叠;否则,不重叠。本发明中步长△θ为0.01°。
为了得到更紧密的结果,本发明按照非凸形骨料颗粒的粒径从大到小排序,使越大的粒子,越先进行堆积。
在模拟区域内进行实验,生成一个骨料,随机投放数千个点,根据本步骤的重叠判断方法可以识别出那些随机点与骨料的相对空间位置关系,图5展示了本发明的重叠判断方法识别出落在骨料外面的那些随机投放点。如图6所示,为骨料面积分数为40%时的采用本发明公开的方法建立的混凝土细观结构模型。
步骤5、根据步骤4得到的结果采用Monte Carlo随机点采样法和统计原理测定界面浓度,包括如下步骤:在模拟空间内大量投放随机点,按照随机点所处的空间位置将随机点分为骨料边界内随机点、骨料边界外随机点、界面最外层轮廓线内随机点和界面最外层轮廓线外随机点;
统计每种随机点的个数,得到界面浓度η:
其中Nsum为随机点总数,Nin为界面最外层轮廓线内随机点总数,Nout为骨料边界外随机点总数。
如图4所示,为二维非凸形随机骨料周围界面浓度测定的流程图。
图7是在不同界面厚度下二维非凸形骨料含量对界面浓度影响的数值结果对比。
Claims (10)
1.一种二维非凸形随机骨料周围界面浓度的测定方法,其特征在于,包括如下步骤:
(1)基于富勒公式推导的瓦拉文公式获取在富勒级配下各粒径范围的二维骨料数量密度:
其中Pc(D<D0)为实际骨料粒径D小于D0的骨料密度,Dmax为最大骨料颗粒直径;
(2)基于简单椭圆曲线通过单峰值的伸缩因子变形得到二维非凸形骨料模型,所述二维非凸形骨料模型的椭圆基础参数方程为:
其中u∈(0,2π),α为椭圆上点的旋转角度,α∈(0,π),a为椭圆的半长轴,b为椭圆的半短轴;
所述二维非凸形骨料模型的伸缩因子函数如下:
其中函数g(s)的表达式如下:
其中r1为峰值区域边界,r2为支区域边界;
(3)采用伸缩矩阵获取非凸形骨料周围界面的细观结构信息;
(4)根据步骤1计算出的二维非凸形骨料颗粒的数量密度和步骤2给出的骨料模型生成骨料颗粒、进行骨料重叠判断并进行随机堆积;
(5)根据步骤4得到的结果采用Monte Carlo随机点采样法和统计原理测定界面浓度。
2.根据权利要求1所述的二维非凸形随机骨料周围界面浓度的测定方法,其特征在于,所述步骤(3)包括如下步骤:
(31)将步骤(2)建立的二维非凸形骨料模型的椭圆基础参数方程中的参数u在(0,2π)区间内拆分为n份,子区间中心为a0i,其中i=1,2,...,n;
(32)对第i个子区间随机生成两个轴向的伸缩因子β1i、β2i,构成伸缩矩阵Di:
Di=diag(β1i,β2i),
以子区间中心a0i为变形峰值向内收缩,得到骨料轮廓线散点P(tg)和界面的外边界轮廓线散点P(ts):
其中椭圆中心E为单位矩阵,γ为二维非凸形骨料周围界面层厚度;
骨料轮廓线所围的面积S(tg)和界面的外边界轮廓线所围的面积S(ts)分别为:
3.根据权利要求1所述的二维非凸形随机骨料周围界面浓度的测定方法,其特征在于,所述步骤(4)中生成骨料颗粒包括如下步骤:
(41)初始化骨料颗粒数量m=0和骨料颗粒的总面积s=0;
(42)根据非凸形骨料的尺寸大小按照步骤(1)计算二维非凸形骨料的数量密度Pc,从而得到骨料的数量N=Pc·Sc,其中Sc为模拟的样本面积;在粗骨料最大粒径下随机生成半长轴aj和半短轴bj;
(43)在模拟区域内随机生成一点Oj(x0j,y0j)作为生成二维非凸形骨料模型的椭圆基础的椭圆中心;
(44)判断点Oj是否在所有已生成的颗粒的覆盖椭圆的外部;如果否,转步骤(43)重新生成椭圆中心;若是,计算Oj与所有已生成骨料边界的距离,判断其中最小的距离d是否大于半长轴aj;如果d<aj,表示骨料颗粒有重叠,转步骤(43)重新生成椭圆中心;如果d≥aj,以Oj(x0j,y0j)为中心、aj为半长轴、bj为半短轴生成椭圆,根据步骤(2)对所述椭圆进行随机变形得到骨料颗粒Qj;判断Qj是否与已经生成的骨料颗粒重叠,如果重叠,则跳转至步骤(43)重新生成非凸形骨料颗粒的空间位置;根据步骤(3)计算骨料颗粒Qj的面积S0(tg),s=s+S0(tg),m=m+1;
(45)判断是否满足迭代收敛条件:m>N或s>S,如果满足,结束;否则转步骤(43);其中S为骨料颗粒最大总面积。
4.根据权利要求3所述的二维非凸形随机骨料周围界面浓度的测定方法,其特征在于,所述步骤(44)中判断骨料颗粒Qj是否与已生成的骨料颗粒Qk重叠,包括如下步骤:
(441)A为非凸骨料颗粒Qj边界的组成点,Oj(x0j,y0j)为Qj中心点;令θ为OjA与椭圆基础长轴的夹角;OjA与已生成的骨料颗粒Qk边界相交于点B;
(442)初始时令θ=0,计算Oj(x0j,y0j)与A、B之间的距离|OjA|和|OjB|,如果|OjA|>|OjB|,则Qj与Qk重叠;否则,不重叠;
(443)以步长△θ修改θ的值:θ=θ+△θ,0≤θ≤360°,再次计算Oj(x0j,y0j)与A、B之间的距离|OjA|和|OjB|,如果|OjA|>|OjB|,则Qj与Qk重叠;否则,不重叠。
5.根据权利要求1所述的二维非凸形随机骨料周围界面浓度的测定方法,其特征在于,所述步骤(5)包括:在模拟空间内大量投放随机点,按照随机点所处的空间位置将随机点分为骨料边界内随机点、骨料边界外随机点、界面最外层轮廓线内随机点和界面最外层轮廓线外随机点;
统计每种随机点的个数,得到界面浓度η:
其中Nsum为随机点总数,Nin为界面最外层轮廓线内随机点总数,Nout为骨料边界外随机点总数。
6.根据权利要求2所述的二维非凸形随机骨料周围界面浓度的测定方法,其特征在于,所述步骤(31)中参数u在(0,2π)区间内进行n等分,分为n个子区间。
7.根据权利要求2所述的二维非凸形随机骨料周围界面浓度的测定方法,其特征在于,所述步骤(32)中伸缩因子β1i、β2i在[-0.16,-0.1]区间内取值。
8.根据权利要求2所述的二维非凸形随机骨料周围界面浓度的测定方法,其特征在于,二维非凸形骨料周围界面层厚度γ为0.003毫米。
9.根据权利要求1所述的二维非凸形随机骨料周围界面浓度的测定方法,其特征在于,所述步骤(4)中按照非凸形骨料颗粒的粒径从大到小排序,使越大的粒子,越先进行堆积。
10.根据权利要求4所述的二维非凸形随机骨料周围界面浓度的测定方法,其特征在于,步长△θ为0.01°。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710068971.3A CN106644868B (zh) | 2017-02-08 | 2017-02-08 | 一种二维非凸形随机骨料周围界面浓度的测定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710068971.3A CN106644868B (zh) | 2017-02-08 | 2017-02-08 | 一种二维非凸形随机骨料周围界面浓度的测定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106644868A true CN106644868A (zh) | 2017-05-10 |
CN106644868B CN106644868B (zh) | 2019-06-21 |
Family
ID=58844673
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710068971.3A Active CN106644868B (zh) | 2017-02-08 | 2017-02-08 | 一种二维非凸形随机骨料周围界面浓度的测定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106644868B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109979009A (zh) * | 2019-03-04 | 2019-07-05 | 天津大学 | 一种具有功能梯度的胶结颗粒材料三维细观数字模型重构方法 |
CN110210178A (zh) * | 2019-06-26 | 2019-09-06 | 西安理工大学 | 一种基于Python再生混凝土三维随机球形骨料模型的构建方法 |
CN111402230A (zh) * | 2020-03-16 | 2020-07-10 | 徐州医科大学 | 基于Monte Carlo方法的股骨髓腔形态特征分析方法 |
CN111414677A (zh) * | 2020-03-03 | 2020-07-14 | 河海大学 | 一种多相高含量颗粒材料的传输性能计算方法 |
CN111681722A (zh) * | 2020-06-03 | 2020-09-18 | 河海大学 | 一种含有高堆积密度椭球颗粒的多孔材料曲折度测定方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006078281A (ja) * | 2004-09-08 | 2006-03-23 | Nippon Telegr & Teleph Corp <Ntt> | 特定粒径をもつ浮遊粒子の個数濃度の評価方法 |
CN104732096A (zh) * | 2015-03-31 | 2015-06-24 | 河海大学 | 非均质各向异性硬化粒子周围界面体积分数的计算方法 |
-
2017
- 2017-02-08 CN CN201710068971.3A patent/CN106644868B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006078281A (ja) * | 2004-09-08 | 2006-03-23 | Nippon Telegr & Teleph Corp <Ntt> | 特定粒径をもつ浮遊粒子の個数濃度の評価方法 |
CN104732096A (zh) * | 2015-03-31 | 2015-06-24 | 河海大学 | 非均质各向异性硬化粒子周围界面体积分数的计算方法 |
Non-Patent Citations (3)
Title |
---|
宋来忠, 彭刚, 姜袁: "不规则颗粒随机分布区域的数值仿真算法", 《系统仿真学报》 * |
楼应平等: "混凝土中骨料分布模拟及界面体积百分数", 《城市建设理论研究(电子版)》 * |
马怀发,书贞,陈厚群: "一种混凝土随机凸多边形骨料模型生成方法", 《中国水利水电科学研究院学报》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109979009A (zh) * | 2019-03-04 | 2019-07-05 | 天津大学 | 一种具有功能梯度的胶结颗粒材料三维细观数字模型重构方法 |
CN109979009B (zh) * | 2019-03-04 | 2023-03-24 | 天津大学 | 一种具有功能梯度的胶结颗粒材料三维细观数字模型重构方法 |
CN110210178A (zh) * | 2019-06-26 | 2019-09-06 | 西安理工大学 | 一种基于Python再生混凝土三维随机球形骨料模型的构建方法 |
CN111414677A (zh) * | 2020-03-03 | 2020-07-14 | 河海大学 | 一种多相高含量颗粒材料的传输性能计算方法 |
CN111414677B (zh) * | 2020-03-03 | 2022-10-14 | 河海大学 | 一种多相高含量颗粒材料的传输性能计算方法 |
CN111402230A (zh) * | 2020-03-16 | 2020-07-10 | 徐州医科大学 | 基于Monte Carlo方法的股骨髓腔形态特征分析方法 |
CN111681722A (zh) * | 2020-06-03 | 2020-09-18 | 河海大学 | 一种含有高堆积密度椭球颗粒的多孔材料曲折度测定方法 |
CN111681722B (zh) * | 2020-06-03 | 2022-07-29 | 河海大学 | 一种含有高堆积密度椭球颗粒的多孔材料曲折度测定方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106644868B (zh) | 2019-06-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106644868A (zh) | 一种二维非凸形随机骨料周围界面浓度的测定方法 | |
CN111832102B (zh) | 一种高维随机场条件下的新型复合材料结构优化设计方法 | |
CN104866652A (zh) | 一种基于abaqus的喷丸强化变形的有限元模拟方法 | |
CN107817199A (zh) | 一种致密砂岩多尺度孔隙模型的构建方法和应用 | |
CN108959794A (zh) | 一种基于深度学习的结构频响动力学模型修正方法 | |
CN104330343B (zh) | 致密砂岩微观孔喉数据的处理方法及装置 | |
CN103196916B (zh) | 一种机翼壁板喷丸表面质量检测方法和检测样板 | |
CN106250933A (zh) | 基于fpga的数据聚类的方法、系统及fpga处理器 | |
CN106407614A (zh) | 一种结合纳米压痕与有限元获取焊缝组织力学参数的方法 | |
CN103971395B (zh) | 一种纤维过滤介质结构的拟态化重建及其性能计算方法 | |
CN106021761B (zh) | 一种汽车覆盖件回弹评测方法 | |
CN114861519B (zh) | 复杂地质条件下初始地应力场加速优化反演方法 | |
CN107204040A (zh) | 多点地质统计学建模方法及装置、计算机存储介质 | |
CN112257197A (zh) | 一种大型铸锻件微缺陷工作应力评估方法 | |
CN113722942A (zh) | 一种考虑3d打印几何缺陷的有限元计算模型 | |
CN106568365A (zh) | 一种球面孔系复合位置度误差的检测与评定方法 | |
CN107507168A (zh) | 一种用于岩体节理裂隙模型检验的迹线图相似度判别方法 | |
CN113868942A (zh) | 一种堰塞体溃决峰值流量快速预测方法 | |
CN106874599B (zh) | 快速生成卵石碎石夹杂的混凝土三维随机骨料模型的方法 | |
CN101241520A (zh) | 有限元建模中基于特征抑制的模型态生成方法 | |
CN108010114A (zh) | 基本图元点云曲面的几何形状识别方法以及特征识别方法 | |
CN104122205A (zh) | 一种利用压痕隆起量测量残余应力的方法 | |
CN114564901B (zh) | 一种结合随机函数的汽车涂层抗石击性仿真评价方法 | |
CN113255069B (zh) | 一种超声波喷丸材料表面形貌预测方法 | |
CN106096136B (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 |