CN110807242B - 考虑骨料分区域填充的钢筋混凝土细观数值模型构建方法 - Google Patents
考虑骨料分区域填充的钢筋混凝土细观数值模型构建方法 Download PDFInfo
- Publication number
- CN110807242B CN110807242B CN201910898408.8A CN201910898408A CN110807242B CN 110807242 B CN110807242 B CN 110807242B CN 201910898408 A CN201910898408 A CN 201910898408A CN 110807242 B CN110807242 B CN 110807242B
- Authority
- CN
- China
- Prior art keywords
- aggregate
- aggregates
- spherical
- reinforced concrete
- region
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/80—Data visualisation
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C60/00—Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Computing Systems (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Software Systems (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Computer Graphics (AREA)
- Data Mining & Analysis (AREA)
- Chemical & Material Sciences (AREA)
- Crystallography & Structural Chemistry (AREA)
- Image Generation (AREA)
- Road Paving Machines (AREA)
Abstract
发明提供考虑骨料分区域填充的钢筋混凝土细观数值模型构建方法。该方法包括指定钢筋混凝土细观数值模型的几何参数、划分填充区域、生成骨料、投放骨料等步骤。该方法能够充分体现钢筋混凝土细观模型中骨料位置和粒径分布的随机规律性。可避免人工干预的主观性并减轻科研设计人员的工作负担。
Description
技术领域
本发明涉及钢筋混凝土细观数值模型,特别涉及考虑骨料分区域填充的钢筋混凝土细观数值模型构建方法。
背景技术
氯离子入侵混凝土导致混凝土结构发生耐久性破坏是一个长期缓慢的过程。开展钢筋混凝土中氯离子长期扩散的物理试验研究相比较于数值模拟而言不仅会耗费大量的时间,还会消耗更多的人力、物力和财力。
随着电子计算机的高速发展,国内外诸多学者提出了一系列的数值模拟方法来研究海洋环境下混凝土结构中氯离子扩散的浓度分布规律。现有研究主要认为混凝土材料的细观组分主要包括:水泥浆体、骨料和二者之间的界面区这三相,并基于二维圆形或三维球形骨料的随机生成和随机投放算法,可最终建立得到混凝土材料的三相细观数值模型。
然而,对于钢筋混凝土材料,其细观结构的组成除了混凝土材料中原有的水泥浆体、骨料以及浆体-骨料之间的界面区这三相之外,还包含钢筋以及浆体-钢筋之间的界面区这两相。现有技术无法充分考虑钢筋存在对钢筋混凝土中不同区域范围内骨料含量及分布的影响,未能体现钢筋混凝土多相复合材料的细观特性,与实际情况存在较大差异性。
因此,亟需建立一种考虑骨料分区域填充的钢筋混凝土细观数值模型构建方法。
发明内容
本发明的目的是提供考虑骨料分区域填充的钢筋混凝土细观数值模型构建方法,以解决现有技术中存在的问题。
为实现本发明目的而采用的技术方案是这样的,考虑骨料分区域填充的钢筋混凝土细观数值模型构建方法,包括以下步骤:
1)指定钢筋混凝土细观数值模型的几何参数。其中,在所述钢筋混凝土细观数值模型中,钢筋混凝土采用长方体区域表示,钢筋采用圆柱体区域表示。长方体的六个面为钢筋混凝土边界,圆柱体的三个面为钢筋边界。钢筋布置在混凝土中,钢筋混凝土边界与钢筋边界之间的区域为骨料填充区域。
2)根据钢筋存在对骨料分布影响的区域范围,将骨料填充区域划分为区域Ⅰ和区域Ⅱ。其中,所述区域Ⅰ为圆环柱体。所述圆环柱体的轴线与步骤1)创建圆柱体的轴线重合。区域Ⅰ内骨料的粒径范围为[dca1(min),dca1(max)]。区域Ⅱ内骨料的粒径范围为[dca2(min),dca2(max)]。
3)将不规则形状的骨料颗粒简化为圆球。根据公式(1)在区域Ⅰ内随机生成m个骨料,并保存全部骨料的粒径结果。根据公式(2)在区域Ⅱ内随机生成n个骨料,并保存全部骨料的粒径结果。
式中,Pi为粒径小于dca1(i)的骨料的累积频率。unifrd(0,1)为随机数函数,可生成[0,1]之间的随机数。
式中,Pj为粒径小于dca2(j)的骨料的累积频率。
4)根据细观模型的几何尺寸范围,创建钢筋混凝土细观数值模型的长方体边界、钢筋的圆柱边界和区域Ⅰ的圆环柱边界。
5)将步骤3)随机生成得到的区域Ⅰ内的全部m个骨料按照颗粒粒径由大到小降序排列。在区域Ⅰ的圆环柱范围内按照粒径从大到小的顺序依次随机投放区域Ⅰ内的全部m个球形骨料,并保存全部m个球形骨料的球心坐标(xi,yi,zi)及其相应的粒径dca1(i)。其中,第i个骨料同时满足与已经成功投放的全部i-1个球形骨料均不重叠以及与钢筋相互之间不重叠。
6)将步骤3)随机生成得到的区域Ⅱ内的全部n个骨料按照粒径由大到小降序排列。在区域Ⅱ范围内按照粒径从大到小的顺序依次随机投放区域Ⅱ内的全部n个球形骨料,并保存全部n个球形骨料的球心坐标(xj,yj,zj)及其相应的粒径dca2(j)。其中,第j个骨料同时满足与已经成功投放的全部j-1个球形骨料均不重叠以及与钢筋相互之间不重叠。
7)绘制钢筋混凝土细观数值模型的三维图示,并导出区域Ⅰ和区域Ⅱ内全部骨料的球心坐标和粒径大小。
进一步,步骤3)具体包括以下步骤:
3.1)计算区域Ⅰ内球形骨料的总体积:
式中,Arc为区域Ⅰ内骨料的含量。dr为钢筋的直径。Ir为混凝土中钢筋的存在对周围骨料分布的影响范围。lr为钢筋的长度。
3.2)根据式(1)在区域Ⅰ内随机生成第i个骨料的粒径dca1(i)。
3.3)判断随机生成得到的区域Ⅰ内的dca1(i)是否满足dca1(min)≤dca1(i)≤dca1(max)。若满足条件,则进入步骤3.4)。若不满足,则返回步骤3.2)中重新生成dca1(i)。
3.4)计算区域Ⅰ内已经生成的全部球形骨料的累积体积:
3.5)判断已经生成的区域Ⅰ内全部球形骨料的累积体积vca1(i)是否满足vca1(i)≥vca1,若满足条件,则停止运算,并保存区域Ⅰ内已经生成的全部m个球形骨料的粒径dca1(i)。若不满足,则返回步骤3.2)继续生成dca1(i),并重复步骤3.2)~3.5),直至满足条件。
3.6)计算钢筋混凝土细观数值模型区域Ⅱ内球形骨料的总体积:
vca2=Vca·l1·l2·l3-vca1 (5)
式中,Vca为骨料的体积分数。l1为钢筋混凝土的长度。l2为钢筋混凝土的宽度。l3为钢筋混凝土的高度。
3.7)根据公式(2)在区域Ⅱ内随机生成第j个骨料的粒径dca2(j)。
3.8)判断随机生成得到的区域Ⅱ内dca2(j)是否满足dca2(min)≤dca2(j)≤dca2(max),若满足条件,则进入步骤3.9)。若不满足,则返回步骤3.7)中重新生成dca2(j)。
3.9)计算区域Ⅱ内已经生成的全部球形骨料的累积体积:
3.10)判断已经生成的区域Ⅱ内全部球形骨料的累积体积vca2(j)是否满足vca2(j)≥vca2。若满足条件,则停止运算,并保存区域Ⅱ内已经生成的全部n个球形骨料的粒径dca2(j)。若不满足,则返回步骤3.7)中继续生成dca2(j),并重复步骤3.7)~步骤3.9)的运算,直至满足条件。
进一步,步骤5)中,球心坐标投放范围为xi∈[(dr/2)·cos(unifrd[0,2π]),(dr/2+Ir)·cos(unifrd[0,2π])],yi∈[(dr/2)·sin(unifrd[0,2π]),(dr/2+Ir)·sin(unifrd[0,2π])],zi∈(当lr≤lrc时,zi∈[zrb,lr];当lr>lrc时,zi∈[zrb,lrc]),在上述区域内随机的生成第i个球形骨料的球心坐标(xi,yi,zi)。其中,dr为钢筋直径,zrb为钢筋底部圆心的z坐标,lrc为沿着钢筋上下两圆面圆心方向的钢筋混凝土试件的高度。
进一步,步骤5)中,第i个球形骨料与其余全部i-1个球形骨料球心距离的空间向量[Si]如式(7)所示。第i个球形骨料与其余全部i-1个球形骨料之间最小球心距离的空间向量[Si]min如式(8)。全部m个球形骨料的球心与钢筋圆柱体的等效圆心之间距离的空间向量[Sr]如式(9)所示。全部m个球形骨料的球心与钢筋圆柱体的等效圆心之间最小距离的空间向量[Sr]min如式(10)所示。其中,[Si]≥[Si]min且[Sr]≥[Sr]min。
式中,tITZ为骨料与水泥浆之间界面区的厚度;S1,r为第1个球形骨料的球心与钢筋的等效圆心之间的距离。S2,r为第2个球形骨料的球心与钢筋的等效圆心之间的距离。Si-1,r为第i-1个球形骨料的球心与钢筋的等效圆心之间的距离。Si,r为第i个球形骨料的球心与钢筋的等效圆心之间的距离。
本发明的技术效果是毋庸置疑的:
A.能够充分体现钢筋混凝土细观模型中骨料位置和粒径分布的随机规律性;
B.可避免人工干预的主观性并减轻科研设计人员的工作负担;
C.可用于模拟钢筋混凝土原型尺度构件,以便进一步开展原型尺度构件的力学强度和腐蚀性离子侵蚀数值模拟。
附图说明
图1为骨料颗粒随机生成算法流程图;
图2为骨料填充区域划分示意图;
图3为骨料颗粒位置关系示意图;
图4为骨料颗粒随机投放算法流程图;
图5为钢筋混凝土细观数值模型示意图。
具体实施方式
下面结合实施例对本发明作进一步说明,但不应该理解为本发明上述主题范围仅限于下述实施例。在不脱离本发明上述技术思想的情况下,根据本领域普通技术知识和惯用手段,做出各种替换和变更,均应包括在本发明的保护范围内。
实施例1:
参见图1,本实施例在混凝土三维球形骨料随机生成和随机投放算法的基础之上,提出了一种考虑骨料分区域填充的钢筋混凝土细观数值模型构建方法,包括以下步骤:
1)指定钢筋混凝土细观数值模型的几何参数。其中,在所述钢筋混凝土细观数值模型中,钢筋混凝土采用长方体区域表示,钢筋采用圆柱体区域表示。长方体的六个面为钢筋混凝土边界,圆柱体的三个面为钢筋边界。钢筋布置在混凝土中,钢筋混凝土边界与钢筋边界之间的区域为骨料填充区域。值得说明的是,现有技术在建立钢筋混凝土细观模型的过程中,仅将钢筋作为一个圆柱形几何体放置在模型中,并未考虑“浆体-钢筋”界面区以及真实情况下钢筋存在会对其周围一定区域范围内骨料的含量及分布造成影响。嵌入混凝土中的钢筋会对其周围大约几毫米的区域范围内骨料的空间随机分布造成影响,降低该区域中骨料的体积分数,减小在该范围内侵蚀性离子的扩散系数,该现象被称之为钢筋对可侵蚀性离子的间接阻挡影响。MATLAB作为一种汇编语言有诸多的优点,主要在于函数调用方便,变量定义简单,并且运算速度快可视化好。在本实施例中,采用MATLAB编写骨料颗粒的生成和投放程序。
2)参见图2,根据钢筋存在对骨料分布影响的区域范围,将骨料填充区域划分为区域Ⅰ和区域Ⅱ。其中,所述区域Ⅰ为圆环柱体。所述圆环柱体的轴线与步骤1)创建圆柱体的轴线重合。区域Ⅰ内骨料的粒径范围为[dca1(min),dca1(max)]。区域Ⅱ内骨料的粒径范围为[dca2(min),dca2(max)]。
3)将不规则形状的骨料颗粒简化为圆球,在区域Ⅰ内随机生成m个骨料,在区域Ⅱ内随机生成n个骨料,并保存全部骨料的粒径结果。
3.1)计算区域Ⅰ内球形骨料的总体积:
式中,Arc为区域Ⅰ内骨料的含量。dr为钢筋的直径。Ir为混凝土中钢筋的存在对周围骨料分布的影响范围。lr为钢筋的长度。
3.2)根据式(1)在区域Ⅰ内随机生成第i个骨料的粒径dca1(i)。
式中,Pi为粒径小于dca1(i)的骨料的累积频率。unifrd(0,1)为随机数函数,可生成[0,1]之间的随机数。
3.3)判断随机生成得到的区域Ⅰ内的dca1(i)是否满足dca1(min)≤dca1(i)≤dca1(max)。若满足条件,则进入步骤3.4)。若不满足,则返回步骤3.2)中重新生成dca1(i)。
3.4)计算区域Ⅰ内已经生成的全部球形骨料的累积体积:
3.5)判断已经生成的区域Ⅰ内全部球形骨料的累积体积vca1(i)是否满足vca1(i)≥vca1,若满足条件,则停止运算,并保存区域Ⅰ内已经生成的全部m个球形骨料的粒径dca1(i)。若不满足,则返回步骤3.2)继续生成dca1(i),并重复步骤3.2)~3.5),直至满足条件。
3.6)计算钢筋混凝土细观数值模型区域Ⅱ内球形骨料的总体积:
vca2=Vca·l1·l2·l3-vca1 (4)
式中,Vca为骨料的体积分数。l1为钢筋混凝土的长度。l2为钢筋混凝土的宽度。l3为钢筋混凝土的高度。
3.7)根据公式(2)在区域Ⅱ内随机生成第j个骨料的粒径dca2(j)。
式中,Pj为粒径小于dca2(j)的骨料的累积频率。
3.8)判断随机生成得到的区域Ⅱ内dca2(j)是否满足dca2(min)≤dca2(j)≤dca2(max),若满足条件,则进入步骤3.9)。若不满足,则返回步骤3.7)中重新生成dca2(j)。
3.9)计算区域Ⅱ内已经生成的全部球形骨料的累积体积:
3.10)判断已经生成的区域Ⅱ内全部球形骨料的累积体积vca2(j)是否满足vca2(j)≥vca2。若满足条件,则停止运算,并保存区域Ⅱ内已经生成的全部n个球形骨料的粒径dca2(j)。若不满足,则返回步骤3.7)中继续生成dca2(j),并重复步骤3.7)~步骤3.9)的运算,直至满足条件。
4)根据细观模型的几何尺寸范围,创建钢筋混凝土细观数值模型的长方体边界、钢筋的圆柱边界和区域Ⅰ的圆环柱边界。
5)参见图4,将步骤3)随机生成得到的区域Ⅰ内的全部m个骨料按照颗粒粒径由大到小降序排列。在区域Ⅰ的圆环柱范围内按照粒径从大到小的顺序依次随机投放区域Ⅰ内的全部m个球形骨料,并保存全部m个球形骨料的球心坐标(xi,yi,zi)及其相应的粒径dca1(i)。其中,第i个骨料同时满足与已经成功投放的全部i-1个球形骨料均不重叠以及与钢筋相互之间不重叠。
球心坐标投放范围为xi∈[(dr/2)·cos(unifrd[0,2π]),(dr/2+Ir)·cos(unifrd[0,2π])],yi∈[(dr/2)·sin(unifrd[0,2π]),(dr/2+Ir)·sin(unifrd[0,2π])],zi∈(当lr≤lrc时,zi∈[zrb,lr];当lr>lrc时,zi∈[zrb,lrc]),在上述区域内随机的生成第i个球形骨料的球心坐标(xi,yi,zi)。其中,dr为钢筋直径,zrb为钢筋底部圆心的z坐标,lrc为沿着钢筋上下两圆面圆心方向的钢筋混凝土试件的高度。
参见图3,第i个球形骨料与其余全部i-1个球形骨料球心距离的空间向量[Si]如式(7)所示。第i个球形骨料与其余全部i-1个球形骨料之间最小球心距离的空间向量[Si]min如式(8)。全部m个球形骨料的球心与钢筋圆柱体的等效圆心之间距离的空间向量[Sr]如式(9)所示。全部m个球形骨料的球心与钢筋圆柱体的等效圆心之间最小距离的空间向量[Sr]min如式(10)所示。其中,[Si]≥[Si]min且[Sr]≥[Sr]min。钢筋的等效圆心是指与随机投放的骨料z坐标相同的钢筋圆截面的圆心。
式中,tITZ为骨料与水泥浆之间界面区的厚度。S1,r为第1个球形骨料的球心与钢筋的等效圆心之间的距离。S2,r为第2个球形骨料的球心与钢筋的等效圆心之间的距离。Si-1,r为第i-1个球形骨料的球心与钢筋的等效圆心之间的距离。Si,r为第i个球形骨料的球心与钢筋的等效圆心之间的距离。
6)将步骤3)随机生成得到的区域Ⅱ内的全部n个骨料按照粒径由大到小降序排列。在区域Ⅱ范围内按照粒径从大到小的顺序依次随机投放区域Ⅱ内的全部n个球形骨料,并保存全部n个球形骨料的球心坐标(xj,yj,zj)及其相应的粒径dca2(j)。其中,第j个骨料同时满足与已经成功投放的全部j-1个球形骨料均不重叠以及与钢筋相互之间不重叠。
7)绘制钢筋混凝土细观数值模型的三维图示,并导出区域Ⅰ和区域Ⅱ内全部骨料的球心坐标和粒径大小。不同情况下钢筋混凝土细观数值模型的三维图示如图5所示。5a中,钢筋直径dr=8mm,骨料体积分数Vca=0.2。5b中,钢筋直径dr=12mm,骨料体积分数Vca=0.3。5c中,钢筋直径dr=16mm,骨料体积分数Vca=0.4。5d中钢筋直径dr=20mm,骨料体积分数Vca=0.5。从图示中可以很明显地看出在不同的钢筋混凝土细观数值模型中考虑由于钢筋的存在对其周围某一区域范围内骨料分布的影响。本实施例充分考虑钢筋存在对钢筋前沿范围内骨料含量和分布的影响,使构建得到的钢筋混凝土细观数值模型更加符合真实情况。
Claims (4)
1.考虑骨料分区域填充的钢筋混凝土细观数值模型构建方法,其特征在于,包括以下步骤:
1)指定钢筋混凝土细观数值模型的几何参数;其中,在所述钢筋混凝土细观数值模型中,钢筋混凝土采用长方体区域表示,钢筋采用圆柱体区域表示;长方体的六个面为钢筋混凝土边界,圆柱体的三个面为钢筋边界;钢筋布置在混凝土中,钢筋混凝土边界与钢筋边界之间的区域为骨料填充区域;
2)根据钢筋存在对骨料分布影响的区域范围,将骨料填充区域划分为区域Ⅰ和区域Ⅱ;其中,所述区域Ⅰ为圆环柱体;所述圆环柱体的轴线与步骤1)创建圆柱体的轴线重合;区域Ⅰ内骨料的粒径范围为[dca1(min),dca1(max)];区域Ⅱ内骨料的粒径范围为[dca2(min),dca2(max)];
3)将不规则形状的骨料颗粒简化为圆球;根据公式(1)在区域Ⅰ内随机生成m个骨料,并保存全部骨料的粒径结果;在区域Ⅰ内随机生成第i个骨料的粒径dca1(i);根据公式(2)在区域Ⅱ内随机生成n个骨料,并保存全部骨料的粒径结果;在区域Ⅱ内随机生成第j个骨料的粒径dca2(j);
式中,Pi为粒径小于dca1(i)的骨料的累积频率;unifrd(0,1)为随机数函数,可生成[0,1]之间的随机数;
式中,Pj为粒径小于dca2(j)的骨料的累积频率;
4)根据细观模型的几何尺寸范围,创建钢筋混凝土细观数值模型的长方体边界、钢筋的圆柱边界和区域Ⅰ的圆环柱边界;
5)将步骤3)随机生成得到的区域Ⅰ内的全部m个骨料按照颗粒粒径由大到小降序排列;在区域Ⅰ的圆环柱范围内按照粒径从大到小的顺序依次随机投放区域Ⅰ内的全部m个球形骨料,并保存全部m个球形骨料的球心坐标(xi,yi,zi)及其相应的粒径dca1(i);其中,第i个骨料同时满足与已经成功投放的全部i-1个球形骨料均不重叠以及与钢筋相互之间不重叠;
6)将步骤3)随机生成得到的区域Ⅱ内的全部n个骨料按照粒径由大到小降序排列;在区域Ⅱ范围内按照粒径从大到小的顺序依次随机投放区域Ⅱ内的全部n个球形骨料,并保存全部n个球形骨料的球心坐标(xj,yj,zj)及其相应的粒径dca2(j);其中,第j个骨料同时满足与已经成功投放的全部j-1个球形骨料均不重叠以及与钢筋相互之间不重叠;
7)绘制钢筋混凝土细观数值模型的三维图示,并导出区域Ⅰ和区域Ⅱ内全部骨料的球心坐标和粒径大小。
2.根据权利要求1所述的考虑骨料分区域填充的钢筋混凝土细观数值模型构建方法,其特征在于,步骤3)具体包括以下步骤:
3.1)计算区域Ⅰ内球形骨料的总体积:
式中,Arc为区域Ⅰ内骨料的含量;dr为钢筋的直径;Ir为混凝土中钢筋的存在对周围骨料分布的影响范围;lr为钢筋的长度;
3.2)根据式(1)在区域Ⅰ内随机生成第i个骨料的粒径dca1(i);
3.3)判断随机生成得到的区域Ⅰ内的dca1(i)是否满足dca1(min)≤dca1(i)≤dca1(max);若满足条件,则进入步骤3.4);若不满足,则返回步骤3.2)中重新生成dca1(i);
3.4)计算区域Ⅰ内已经生成的全部球形骨料的累积体积:
3.5)判断已经生成的区域Ⅰ内全部球形骨料的累积体积vca1(i)是否满足vca1(i)≥vca1,若满足条件,则停止运算,并保存区域Ⅰ内已经生成的全部m个球形骨料的粒径;若不满足,则返回步骤3.2)继续生成dca1(i),并重复步骤3.2)~3.5),直至满足条件;
3.6)计算钢筋混凝土细观数值模型区域Ⅱ内球形骨料的总体积:
vca2=Vca·l1·l2·l3-vca1 (5)
式中,Vca为骨料的体积分数;l1为钢筋混凝土的长度;l2为钢筋混凝土的宽度;l3为钢筋混凝土的高度;
3.7)根据公式(2)在区域Ⅱ内随机生成第j个骨料的粒径dca2(j);
3.8)判断随机生成得到的区域Ⅱ内dca2(j)是否满足dca2(min)≤dca2(j)≤dca2(max),若满足条件,则进入步骤3.9);若不满足,则返回步骤3.7)中重新生成dca2(j);
3.9)计算区域Ⅱ内已经生成的全部球形骨料的累积体积:
3.10)判断已经生成的区域Ⅱ内全部球形骨料的累积体积vca2(j)是否满足vca2(j)≥vca2;若满足条件,则停止运算,并保存区域Ⅱ内已经生成的全部n个球形骨料的粒径;若不满足,则返回步骤3.7)中继续生成dca2(j),并重复步骤3.7)~步骤3.9)的运算,直至满足条件。
3.根据权利要求1所述的考虑骨料分区域填充的钢筋混凝土细观数值模型构建方法,其特征在于:步骤5)中,球心坐标投放范围为xi∈[(dr/2)·cos(unifrd[0,2π]),(dr/2+Ir)·cos(unifrd[0,2π])],yi∈[(dr/2)·sin(unifrd[0,2π]),(dr/2+Ir)·sin(unifrd[0,2π])],zi∈(当lr≤lrc时,zi∈[zrb,lr];当lr>lrc时,zi∈[zrb,lrc]),在上述区域内随机的生成第i个球形骨料的球心坐标(xi,yi,zi);其中,dr为钢筋直径,zrb为钢筋底部圆心的z坐标,lrc为沿着钢筋上下两圆面圆心方向的钢筋混凝土试件的高度。
4.根据权利要求1所述的考虑骨料分区域填充的钢筋混凝土细观数值模型构建方法,其特征在于:步骤5)中,第i个球形骨料与其余全部i-1个球形骨料球心距离的空间向量[Si]如式(7)所示;第i个球形骨料与其余全部i-1个球形骨料之间最小球心距离的空间向量[Si]min如式(8);全部m个球形骨料的球心与钢筋圆柱体的等效圆心之间距离的空间向量[Sr]如式(9)所示;全部m个球形骨料的球心与钢筋圆柱体的等效圆心之间最小距离的空间向量[Sr]min如式(10)所示;其中,[Si]≥[Si]min且[Sr]≥[Sr]min;
式中,tITZ为骨料与水泥浆之间界面区的厚度;S1,r为第1个球形骨料的球心与钢筋的等效圆心之间的距离;S2,r为第2个球形骨料的球心与钢筋的等效圆心之间的距离;Si-1,r为第i-1个球形骨料的球心与钢筋的等效圆心之间的距离;Si,r为第i个球形骨料的球心与钢筋的等效圆心之间的距离。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910898408.8A CN110807242B (zh) | 2019-09-23 | 2019-09-23 | 考虑骨料分区域填充的钢筋混凝土细观数值模型构建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910898408.8A CN110807242B (zh) | 2019-09-23 | 2019-09-23 | 考虑骨料分区域填充的钢筋混凝土细观数值模型构建方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110807242A CN110807242A (zh) | 2020-02-18 |
CN110807242B true CN110807242B (zh) | 2022-08-12 |
Family
ID=69487602
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910898408.8A Active CN110807242B (zh) | 2019-09-23 | 2019-09-23 | 考虑骨料分区域填充的钢筋混凝土细观数值模型构建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110807242B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112052539B (zh) * | 2020-09-09 | 2022-10-18 | 北京理工大学 | 一种基于Voronoi技术的钢筋混凝土三维细观快速并行建模方法 |
CN112429165B (zh) * | 2020-11-17 | 2021-08-24 | 哈尔滨理工大学 | 一种船用减摇鳍的鳍体实时自动避障方法 |
CN112632814B (zh) * | 2020-12-04 | 2024-09-10 | 西安理工大学 | 一种基于椭球体随机骨料的再生混凝土三维模型构建方法 |
CN115205486B (zh) * | 2022-07-14 | 2023-04-07 | 重庆交通大学 | 混凝土三维细观数值模型构建方法及氯离子传输仿真系统 |
CN117248546B (zh) * | 2023-11-16 | 2024-02-02 | 西南交通大学 | 一种无需挖槽的锚索框架梁快速施工方法及结构 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000235589A (ja) * | 1999-02-15 | 2000-08-29 | Toray Ind Inc | 物品の表示方法および表示装置 |
CN102136157B (zh) * | 2011-03-07 | 2012-09-26 | 河海大学 | 一种混凝土三维细观仿真模型及其建立方法 |
US9865085B1 (en) * | 2017-05-08 | 2018-01-09 | Analytical Graphics, Inc. | Systems and methods for 3D modeling using skipping heuristics and fusing |
CN108344634B (zh) * | 2018-01-10 | 2020-07-14 | 东南大学 | 一种钢筋混凝土细观三维模型建模方法 |
CN110118775B (zh) * | 2019-05-10 | 2021-11-19 | 重庆交通大学 | 厂拌水泥稳定碎石骨料组成快速检测方法 |
CN110210178A (zh) * | 2019-06-26 | 2019-09-06 | 西安理工大学 | 一种基于Python再生混凝土三维随机球形骨料模型的构建方法 |
-
2019
- 2019-09-23 CN CN201910898408.8A patent/CN110807242B/zh active Active
Non-Patent Citations (2)
Title |
---|
"考虑粗骨料和钢筋共同作用的混凝土构件氯离子扩散试验研究";吕彦伟;《中国优秀硕士学位论文全文数据库》;20190430;第B015-413页 * |
基于MCM骨料模型的氯离子侵蚀混凝土模拟研究;于磊等;《西北农林科技大学学报(自然科学版)》;20170531(第05期);第221-227+234页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110807242A (zh) | 2020-02-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110807242B (zh) | 考虑骨料分区域填充的钢筋混凝土细观数值模型构建方法 | |
CN107918706B (zh) | 一种自适应的混凝土细观建模方法 | |
CN112052539B (zh) | 一种基于Voronoi技术的钢筋混凝土三维细观快速并行建模方法 | |
Xu et al. | Random non-convex particle model for the fraction of interfacial transition zones (ITZs) in fully-graded concrete | |
Sheng et al. | Residual stress field induced by shot peening based on random-shots for 7075 aluminum alloy | |
CN114186434B (zh) | 基于逐级划分的沥青混合料三维细观结构模型构建方法 | |
CN110162849B (zh) | 一种混杂纤维混凝土的建模方法 | |
CN104484538A (zh) | 一种基于有限元分析的喷丸强化处理工艺参数确定方法 | |
Moreno et al. | 3D numerical simulation of particle-particle collisions in saltation mode near stream beds | |
CN117195382B (zh) | 一种混凝土细观模型的构建方法 | |
CN113591195A (zh) | 一种基于随机多边形骨料的再生混凝土细观模型构建方法 | |
Zeng et al. | 3D mesoscale investigation of non-uniform steel corrosion in reinforced concrete under chloride environments | |
CN117077501B (zh) | 一种非均质钢渣混凝土渗透系数计算方法 | |
CN117113686A (zh) | 基于坐标网格的混凝土骨料二维几何结构随机生成方法 | |
CN117034806A (zh) | 一种基于流固耦合的埋地管线风险评估方法和装置 | |
CN114925590B (zh) | 多用途二维不规则形状集料随机生成方法及模型生成方法 | |
CN115205486B (zh) | 混凝土三维细观数值模型构建方法及氯离子传输仿真系统 | |
CN106874599A (zh) | 快速生成卵石碎石夹杂的混凝土三维随机骨料模型的方法 | |
Ding et al. | A novel method for generation of particle packing model used in numerical simulation for the mechanical behavior of multi-component material | |
CN111177969B (zh) | 一种可控骨料间距的二维随机骨料生成及投放方法 | |
CN116070324A (zh) | 一种含预制裂缝的混凝土梁细观模型生成方法及系统 | |
CN113627040A (zh) | 一种非均质边坡稳定性分析方法 | |
Shi et al. | Research on microstructure statistical model of soil-rock mixture | |
Chen et al. | A bidirectional walk method of random polygon aggregate for two-dimensional concrete model | |
CN110502819A (zh) | 混凝土桥梁条形构件外封闭箍筋的计算机辅助设计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |