CN115205486B - 混凝土三维细观数值模型构建方法及氯离子传输仿真系统 - Google Patents

混凝土三维细观数值模型构建方法及氯离子传输仿真系统 Download PDF

Info

Publication number
CN115205486B
CN115205486B CN202210825040.4A CN202210825040A CN115205486B CN 115205486 B CN115205486 B CN 115205486B CN 202210825040 A CN202210825040 A CN 202210825040A CN 115205486 B CN115205486 B CN 115205486B
Authority
CN
China
Prior art keywords
aggregate
polyhedral
concrete
dimensional
ith
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
CN202210825040.4A
Other languages
English (en)
Other versions
CN115205486A (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.)
Chongqing Jiaotong University
Original Assignee
Chongqing Jiaotong University
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 Chongqing Jiaotong University filed Critical Chongqing Jiaotong University
Priority to CN202210825040.4A priority Critical patent/CN115205486B/zh
Publication of CN115205486A publication Critical patent/CN115205486A/zh
Application granted granted Critical
Publication of CN115205486B publication Critical patent/CN115205486B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • G06T17/10—Constructive solid geometry [CSG] using solid primitives, e.g. cylinders, cubes
    • 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
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/26—Composites

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Computing Systems (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Computer Graphics (AREA)
  • Software Systems (AREA)
  • Processing Or Creating Images (AREA)

Abstract

本发明公开一种混凝土三维细观数值模型构建方法及氯离子传输仿真系统，该方法首先利用Monte‑Carlo法、Fuller级配关系、德洛内三角网算法随机生成混凝土三维多面体集料，然后按照集料粒径从大到小排序后进行随机投放，在投放的过程生成多面体集料的界面区并根据重叠判别准则建立全部多面体集料及其界面区，完成混凝土三维细观数值模型的构建；本发明所生成的多面体集料在形状、粒径、分布等具有随机性，更符合实际情况，且计算效率高；同时，基于混凝土三维细观数值模型建立的氯离子传输仿真系统，可克服混凝土内氯离子传输数值仿真研究当中大多数将集料形状假定为规则球形或椭球形等与实际情况不符的问题。

Description

混凝土三维细观数值模型构建方法及氯离子传输仿真系统
技术领域
本发明涉及混凝土细观模型技术领域，特别是一种混凝土三维细观数值模型构建方法及氯离子传输仿真系统。
背景技术
混凝土在工程建设中有着广泛地应用，但在海洋环境条件下，钢筋混凝土耐久性破坏一直是威胁结构使用寿命的一个重要且亟待解决的问题。特别是沿海地区钢筋混凝土结构耐久性破坏现象极其严重，海洋环境下钢筋混凝土结构大部分都属于露天工程，所处服役环境复杂，所受影响因素众多。
混凝土是一种复杂的多相复合性物质，现阶段数值模拟研究中普遍将混凝土的细观结构考虑为由水泥浆体、集料以及界面区共同组成的三相非均质材料。其集料的随机生成及投放技术是混凝土耐久性计算研究的重要部分。
现有的混凝土细观数值模型研究成果中，多数学者将混凝土中的集料简化为圆形、球形或凸多边形、凸多面体。
但是实际当中的集料形状不一定是完全的凸多面体，有可能局部是凹陷的形状，在生成随机凸多面体集料过程中，随机凸多面体集料的各顶点是在其外接球的球面上，这与实际情况不相符。并且，上述模型会对氯离子在混凝土中的传输过程产生影响，因此现有的混凝土三维细观模型的构建方式具有局限性。
发明内容
有鉴于此，本发明的目的在于提供一种混凝土三维细观数值模型构建方法及氯离子传输仿真系统，该方法利用Monte-Carlo法、Fuller级配关系、德洛内三角网算法生成混凝土多面体集料并构建混凝土三维细观数值模型，然后根据该细观数值模型建立其内氯离子传输仿真系统。
为达到上述目的，本发明提供如下技术方案：
本发明提供的混凝土三维细观数值模型构建方法，包括以下步骤：
S1：构建三维多面体集料模型，根据集料粒径范围生成不同粒径大小的三维多面体集料；
S2：将三维多面体集料按照粒径di从大到小进行排序，并按照排列顺序随机投放第i个多面体集料。通过Monte-Carlo法随机在L×M×H的混凝土三维细观模型线框中生成第i个多面体集料最小外接球的球心坐标；
重叠判别：判断第i个多面体集料是否会与之前已投放的i-1个集料重叠。采用多面体集料最小外接球球心距Li,j与多面体集料最小外接球的半径之和Lr-i,j二者的关系来进行重叠判别，记录并保存满足条件Li,j≥Lr-i,j的多面体集料几何信息；
其中，Li,j为相邻两多面体集料最小外接球球心距；Lr-i,j为相邻两多面体集料最小外接球半径之和；
S3：构建并生成第i个多面体集料的界面区，记录第i个多面体集料的界面区信息；
S4：直到全部多面体集料投放完毕，即完成混凝土三维细观数值模型的构建。
进一步，所述步骤S1中构建三维多面体集料模型按照以下步骤进行：
步骤(1)：确定三维多面体集料模型中的基本参数；
步骤(2)：随机生成第i个多面体集料粒径di；
步骤(3)：判断集料粒径di是否满足dmin≤di≤dmax，若满足，继续执行步骤(4)，若不满足，则返回重新执行步骤(2)；其中，集料最小粒径dmin、最大粒径dmax；
步骤(4)：生成粒径为di的多面体集料模型；
步骤(5)：计算第i个集料体积vai并按照以下方式判断是否满足集料体积占比Ca0：
判断第i个集料体积占比Cai是否满足Cai≥Ca0，若满足，则继续执行步骤(6)，若不满足，则返回重新执行步骤(2)；其中，Cai为第i个集料体积占比；vai为第i个多面体集料体积；其中，Ca0表示集料体积占比，即多面体集料体积占与之对应的外接球体积；
步骤(6)：记录并保存第i个多面体集料信息；
步骤(7)：按照以下方式判断集料是否生成完毕：
判断集料累计体积va是否满足va≥vas，若不满足，令i＝i+1，返回步骤(2)；若满足，则三维多面体集料生成完毕，其中，vas为混凝土集料体积。
进一步，所述步骤(2)中的第i个多面体集料粒径di按照以下步骤生成：
根据Monte-Carlo算法生成区间[0,1]内的任意随机数Fi，并根据Fuller级配公式(1)求解得到符合Fuller级配曲线的多面体集料粒径di；
Figure BDA0003746142930000031
其中，di表示第i个集料的粒径。
进一步，所述步骤(4)的多面体集料模型按照以下方式构建：
以(0,0,0)为中心点坐标生成直径为di的球，在球内及球面上任意取点并根据所取的点采用德洛内三角网算法构建粒径为di的多面体集料。
进一步，所述步骤(5)中的集料体积占比Cai按照以下公式计算：
Cai＝vai/vsi＝vai/(π*di 3/6) (2)
式中：Cai为第i个集料体积占比；vai为第i个多面体集料体积；vsi为第i个集料对应的最小外接球的体积；di为第i个集料粒径。
进一步，所述步骤S2的重叠判别按照以下步骤进行：
按照以下公式计算多面体集料最小外接球球心距Li,j与多面体集料最小外接球的半径之和Lr-i,j：
Figure BDA0003746142930000032
式中：Li,j为相邻两多面体集料最小外接球球心距；Lr-i,j为相邻两多面体集料最小外接球半径之和；
判断Li,j是否满足Li,j≥Lr-i,j，若不满足，返回并重新随机投放下一个多面体集料，若满足，则保存并更新第i个多面体集料的几何信息。
进一步，所述步骤S3按照以下步骤构建并生成第i个多面体集料的界面区：
步骤S31：设置界面区基本参数；
步骤S32：生成第i个多面体集料的界面区，多面体集料界面区厚度按以下公式计算；
titz＝δ*R (4)
式中：titz为界面区厚度；δ为界面区厚度倍数；R为集料最小外接球中心点距离多面体角点的距离；
将多面体集料外接球球心O与集料任一角点A相连，得到集料中心点与角点A的距离ROA，延长ROA的δ倍至A’点，延长距离即为该处界面过渡区厚度titz；
其中，A’点坐标的计算表达式为：
Figure BDA0003746142930000041
式中：(x0,y0,z0)为多面体集料中心点的坐标；(xA,yA,zA)为集料角点A的坐标；(xA’,yA’,zA’)为生成的界面区角点A’的坐标；δ为界面区厚度倍数；
步骤S33：记录第i个多面体集料界面区的信息；
步骤S34：判断集料界面区是否生成完毕：判断i是否满足i≥na，若不满足，则令i＝i+1，返回循环执行步骤S32，若满足，则全部多面体集料的界面区生成完毕。其中，na表示已生成的集料数量。
本发明提供的基于混凝土三维细观数值模型的氯离子传输仿真系统，包括混凝土多面体集料细观数值模型构建模块、氯离子传输参数模块、氯离子传输过程仿真分析模块和仿真结果输出模块；
所述混凝土多面体集料细观数值模型构建模块，用于构建不同集料含量、粒径、界面区特征的混凝土多面体集料细观数值模型；
所述氯离子传输参数模块，用于输入、设置混凝土三维细观数值模型中氯离子传输过程仿真分析的重要参数；
所述氯离子传输过程仿真分析模块，用于分析混凝土多面体集料细观数值模型中氯离子传输浓度的时空分布特征；
所述仿真结果输出模块，用于输出混凝土多面体集料细观数值模型中氯离子分布的数值仿真计算结果。
进一步，所述混凝土多面体集料细观数值模型构建模块包括三维多面体集料生成单元、三维多面体集料投放单元、多面体集料的界面区生成单元；
所述三维多面体集料模型构建单元，用于根据集料粒径范围构建并生成不同粒径大小的三维多面体集料；
所述三维多面体集料投放单元，用于将按照粒径排序的三维多面体集料以不重叠的方式进行随机投放；
所述多面体集料的界面区生成单元，用于生成并保存各个多面体集料界面区域信息。
进一步，所述氯离子传输参数模块包括氯离子传输参数输入单元和氯离子传输参数设置单元；所述氯离子传输参数输入单元，用于输入混凝土三维细观数值模型的表面氯离子浓度、初始氯离子含量和混凝土材料各细观组成相的氯离子扩散系数；
所述氯离子传输参数设置单元，包括氯离子传输维度子单元、氯离子传输方程子单元、氯离子传递属性子单元、网格精密程度划分子单元。所述氯离子传输维度子单元，用于设置混凝土三维细观模型中氯离子传输的维度；所述氯离子传输方程子单元，用于设置混凝土材料各组成相中的氯离子传输方程；所述氯离子传递属性子单元，用于设置混凝土材料各组成相中的氯离子传递属性；所述网格精密程度划分子单元，用于划分混凝土材料各组成相的网格精密程度。
本发明的有益效果在于：
本发明提供的混凝土三维细观数值模型构建方法及氯离子传输仿真系统，该方法从细观尺度出发，将混凝土视为由水泥浆体、集料及界面区共同组成的三相非均质复合材料，基于Monte-Carlo法、Fuller级配关系、德洛内三角网算法等手段，提出生成多面体集料及界面区的混凝土三维细观模型构建方法。根据有限元FEM原理，基于COMSOL稀物质传递模块模拟混凝土三维细观模型内部氯离子的传输过程，建立混凝土三维细观模型内氯离子传输仿真方法及系统。
相比于其他的混凝土细观数值模型，本实施例建立的混凝土三维细观数值模型具有以下优点：
(1)从细观尺度出发，考虑了混凝土材料的不均匀性，建立了水泥浆体、多面体集料及界面区的三相细观数值模型。
(2)基于Monte-Carlo法、Fuller级配关系、德洛内三角网算法等手段，提出混凝土多面体集料的随机生成及投放算法。相比于传统的延凸算法，本发明算法简洁，计算效率高，生成的集料形状、粒径、分布等具有随机性，满足Fuller级配，更符合实际情况。
(3)引入集料体积占比的概念，避免生成畸形混凝土集料，可极大程度地提高混凝土多面体集料的生成及投放效率。
(4)提出界面区厚度倍数，考虑了界面区厚度与集料形状、尺寸的关系，本发明生成的多面体集料，不仅可以是凸多面体，还可以是凹陷形的多面体，更符合混凝土材料中多面体集料的实际情况。
本发明的其他优点、目标和特征在某种程度上将在随后的说明书中进行阐述，并且在某种程度上，基于对下文的考察研究对本领域技术人员而言将是显而易见的，或者可以从本发明的实践中得到教导。本发明的目标和其他优点可以通过下面的说明书来实现和获得。
附图说明
为了使本发明的目的、技术方案和有益效果更加清楚，本发明提供如下附图进行说明：
图1为四面体单元角点坐标图示。
图2为界面区生成的几何图示。
图3为总体算法流程图。
图4为混凝土三维细观模型。
图5为氯离子传输仿真系统中的混凝土三维细观数值模型。
图6为模型网格划分示意图。
图7为混凝土三维细观数值模型中氯离子传输仿真计算结果示意图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步说明，以使本领域的技术人员可以更好的理解本发明并能予以实施，但所举实施例不作为对本发明的限定。
实施例1
如图1至图4所示，本实施例提供的混凝土三维细观数值模型构建方法，包括以下步骤：
S1：构建三维多面体集料模型，并根据集料粒径范围生成不同粒径大小的三维多面体集料，具体过程如下：
步骤(1)：确定混凝土三维细观模型中混凝土试件的基本参数，所述基本参数包括混凝土试件的几何尺寸：试件的长L、宽W、高H，以及集料体积分数Va，集料最小粒径dmin、最大粒径dmax，集料体积占比Ca0(多面体集料体积占与之对应的外接球体积)等基本参数；
步骤(2)：根据混凝土三维细观模型中混凝土试件的尺寸随机生成第i个多面体集料粒径di；
本实施例根据Monte-Carlo算法生成区间[0,1]内的任意随机数Fi，并根据Fuller级配公式(1)求解得到符合Fuller级配曲线的多面体集料粒径di；
Figure BDA0003746142930000071
其中，di表示第i个集料的粒径。
步骤(3)：判断集料粒径di是否满足dmin≤di≤dmax，若满足，继续执行步骤(4)，若不满足，则返回重新执行步骤(2)；
步骤(4)：采用德洛内三角网算法构建粒径为di的多面体集料模型，以(0,0,0)为中心点坐标生成直径为di的球，在球的八个象限上任意取点，保证该球是随机多面体的最小外接球，即可认为该多面体集料的粒径为di；(经观察统计，混凝土集料一般有16～20个顶点，考虑到空间内任意取点时，许多点会在生成的多面体内部，故在球内或球面上任取20个点，采用德洛内三角网算法构建混凝土多面体集料模型。经统计，生成的集料模型约有15～22个顶点，符合实际情况。)
步骤(5)：计算集料体积vai并判断这第i个集料是否满足集料体积占比Ca0的要求，以球心为顶点，将已构建的多面体集料模型划分为多个四面体，基于式(2)计算粒径为di的第i个集料体积vai，为了避免畸形集料的产生和提高投放效率，引入集料体积占比Ca0，即多面体集料体积占与之对应的外接球体积；见式(3)：
Figure BDA0003746142930000072
当已知第i个四面体中四个顶点的坐标：A(0,0,0)，B(xBi,yBi,zBi)，C(xCi,yCi,zCi)，D(xDi,yDi,zDi)时，可利用式(2)计算得到四面体的体积。如图1所示，图1为四面体单元角点坐标图示。
Cai＝vai/vsi＝vai/(π*di 3/6) (3)
式中：Cai为第i个集料体积占比；vai为第i个多面体集料体积；vsi为第i个集料对应的最小外接球的体积；di为第i个集料粒径，本实施例中Ca0取≥0.5。
判断Cai是否满足Cai≥Ca0，若满足，则继续执行步骤(6)，若不满足，则返回重新执行步骤(2)。
步骤(6)：记录并保存第i个多面体集料信息(包括：粒径、角点坐标、体积等)至统一数组中，方便后续调用。
步骤(7)：判断集料是否生成完毕，计算已生成的集料体积va＝va(i-1)+vai，而混凝土内需生成的集料总体积vas如下：
vas＝L*M*H*Va (4)
式中：vas为混凝土集料总体积；L、M、H分别为混凝土细观模型的长、宽、高；Va为集料体积分数。
判断va是否满足va≥vas，若不满足，令i＝i+1，返回步骤(2)；若满足，则混凝土多面体集料生成完毕。
本实施例提供的混凝土三维多面体集料模型是基于Monte-Carlo法、Fuller级配关系、德洛内三角网算法生成得到的，相比于其他方式建立的规则形状集料模型(圆形、球形、椭圆形、椭球体等)以及传统延凸算法构建的不规则形状集料模型，具有算法简洁，计算效率高，生成的集料符合实际混凝土集料形状等优点。
S2：混凝土三维细观数值模型中多面体集料的随机投放，具体过程如下：
步骤S21：设置集料投放基本参数，所述集料投放基本参数包括集料投放边界距离ba，所述集料投放边界距离ba为净浆层厚度，本实施例的净浆层厚度取ba＝0.1mm，在所述净浆层厚度范围内无集料存在；
步骤S22：将生成的全部多面体集料按照粒径di从大到小降序排序，保证投放集料时按照粒径从大到小投放，可极大程度提高集料的投放效率；
步骤S23：随机投放第i个多面体集料：基于Monte-Carlo法随机在L×M×H的混凝土三维数值模型线框中生成第i个多面体集料最小外接球的球心坐标，其横、纵、竖坐标的取值范围如下：
Figure BDA0003746142930000091
步骤S24：重叠判别：判断第i个多面体集料是否会与之前已投放的i-1个集料重叠。采用多面体集料最小外接球球心距Li,j与多面体集料最小外接球的半径之和Lr-i,j二者的关系来进行重叠判别。
Figure BDA0003746142930000092
式中：Li,j为相邻两多面体集料最小外接球球心距；Lr-i,j为相邻两多面体集料最小外接球半径之和；
判断计算得到的Li,j是否满足Li,j≥Lr-i,j，若不满足，返回执行步骤S23，若满足，则继续执行步骤S25；
步骤S25：保存并更新第i个多面体集料的几何信息，所述几何信息包括粒径、角点数量、角点坐标、体积、最小外接球球心坐标等，本实施例中的几何信息保存至一个统一数组中，方便后续调用；
步骤S26：判断多面体集料是否投放完毕：判断i是否满足i≥na，若不满足，则令i＝i+1，返回循环执行步骤S23，若满足，则多面体集料投放完毕。
综上，本实施例提出的基于Monte-Carlo原理的混凝土多面体集料随机投放算法，能够更好地模拟实际混凝土中集料的随机分布。投放过程中，根据集料粒径从大到小依次投放，并进行重叠判别，算法简洁，投放效率高。
S3：构建并生成混凝土多面体集料界面区，具体过程如下：
步骤S31：设置界面区基本参数，所述界面区基本参数包括界面区厚度倍数δ等；
步骤S32：生成第i个多面体集料的界面区，多面体集料界面区厚度按以下公式计算；
titz＝δ*R (7)
式中：titz为界面区厚度；δ为界面区厚度倍数；R为集料最小外接球中心点距离多面体角点的距离；
将多面体集料外接球球心O与集料任一角点A相连，得到集料中心点与角点A的距离ROA，延长ROA的δ倍至A’点，延长距离即为该处界面过渡区厚度titz。此时，A’点坐标的计算表达式为：
Figure BDA0003746142930000101
式中：(x0,y0,z0)为多面体集料中心点的坐标；(xA,yA,zA)为集料角点A的坐标；(xA’,yA’,zA’)为生成的界面区角点A’的坐标；δ为界面区厚度倍数。如图2所示，图2为界面区生成的几何图示。
步骤S33：记录第i个多面体集料界面区的信息，所述界面区的信息包括界面区厚度、角点坐标等，本实施例中是将信息保存至统一数组中，方便后续调用；
步骤S34：判断集料界面区是否生成完毕，判断i是否满足i≥na，若不满足，则令i＝i+1，返回循环执行步骤S32，若满足，则全部多面体集料的界面区生成完毕；其中，na表示已生成的集料数量；
本实施例考虑到实际情况下混凝土集料界面区存在重叠、连通等情况，故无需对界面区进行重叠判别，大大提高了计算效率。
综述上述，本实施例提供的混凝土三维细观模型构建方法，具体过程如下：通过随机生成混凝土多面体集料，然后对多面体集料进行随机投放，并构建多面体集料的界面区，如图3所示，图3为总体算法流程图，通过该流程可以得到长、宽、高分别为L、M、H的混凝土多面体集料及界面区三维细观模型，如图4，图4为混凝土三维细观模型。本实施例中的v＝100*100*100mm3,Va＝0.1,damin～damax＝5～31.5mm。
本实施例提出的多面体集料顶点在多面体最小外接球的球面或内部任意位置(相当于是在以球心为坐标原点，集料外接球所在空间坐标系oxyz的8个象限的任意位置都可以作为多面体集料的顶点)，这样更符合实际情况；避免了传统算法中仅可以生成随机凸多面体集料的局限性，传统方法中生成的随机凸多面体集料的各顶点是设置在其最小外接球的球面上，这与实际情况不完全相符，具体理由如下：
由于实际当中的集料形状不一定是完全的凸多面体，有可能局部是凹陷的形状，因此只设置多面体集料角点在其最小外接球球面的方式具有局限性。而本方法可实现在外接球球面和内部任意位置随机生成集料角点坐标，这样建立的集料模型就只是多面体而并不完全是凸多面体。
实施例2
本实施例根据上述混凝土多面体集料随机生成和投放算法以及界面区构建方法，得到混凝土三维细观数值模型，并在混凝土三维细观数值模型内进行氯离子传输仿真，相比传统的混凝土细观数值模型而言，本发明建立的多面体集料及界面区的混凝土三维细观数值模型更符合实际情况。
本实施例提供的基于混凝土三维细观数值模型的氯离子传输仿真系统，包括混凝土多面体集料细观数值模型构建模块、氯离子传输参数模块、氯离子传输过程仿真分析模块和仿真结果输出模块；
所述混凝土多面体集料细观数值模型构建模块，用于构建不同集料含量、粒径、界面区特征的混凝土多面体集料细观数值模型；
所述氯离子传输参数模块，用于输入和设置混凝土三维细观数值模型中氯离子传输过程仿真分析的重要参数；
所述氯离子传输过程仿真分析模块，用于分析混凝土多面体集料细观数值模型中氯离子传输浓度的时空分布特征；
所述仿真结果输出模块，用于输出混凝土多面体集料细观数值模型中氯离子分布的数值仿真计算结果。
其中，所述混凝土多面体集料细观数值模型构建模块包括三维多面体集料生成单元、三维多面体集料投放单元、多面体集料的界面区生成单元；
所述三维多面体集料模型构建单元，用于根据集料粒径范围构建并生成不同粒径大小的三维多面体集料；
所述三维多面体集料投放单元，用于将按照粒径排序的三维多面体集料以不重叠的方式进行随机投放；
所述多面体集料的界面区生成单元，用于生成并保存各个多面体集料界面区域信息。
其中，所述氯离子传输参数模块包括氯离子传输参数输入单元、氯离子传输参数设置单元；
所述氯离子传输参数输入单元，用于输入混凝土三维细观数值模型的表面氯离子浓度(边界条件)，初始氯离子含量(初始条件)，混凝土材料各细观组成相(水泥浆体、集料、界面区)的氯离子扩散系数等重要氯离子传输计算参数；
所述氯离子传输参数设置单元，包括：氯离子传输维度子单元、氯离子传输方程子单元、氯离子传递属性子单元、网格精密程度划分子单元。用于分别设置混凝土三维细观模型中氯离子传输的维度、混凝土材料各组成相(水泥浆体、集料、界面区)中的氯离子传输方程、氯离子传递属性并划分网格精密程度。
本实施例提供的基于混凝土三维细观数值模型的氯离子传输仿真系统，具体按照以下步骤进行构建和使用：
(1)创建仿真环境
在COMSOL Multiphysics软件中新建模型，在模型向导的空间维度中选择三维，在物理场中选择化学物质传递物理场，添加稀物质传递(tds)接口，点击完成，即可完成仿真环境的创建。
(2)在仿真环境中导入并构建混凝土三维细观数值模型
在模型开发器组件的几何模块中建立水泥浆体、集料、界面区模块，分别导入已构建的混凝土多面体集料模型、混凝土多面体集料界面区模型(直接导入2个“.dxf”格式的集料模型和集料界面区模型)。采用“层”的方法，分层建立水泥浆体域、界面区域、集料域并形成联合体，即可在COMSOL仿真环境中构建混凝土三维细观数值模型。如图5所示，图5为氯离子传输仿真系统中的混凝土三维细观数值模型。
(3)输入模型参数
在模型开发器的全局定义参数模块中，输入混凝土三维细观数值模型内氯离子传输数值仿真所需要的参数。包括：参考氯离子扩散系数D28、参考扩散时间t28、表面氯离子浓度Cs(t)、龄期系数m、无因次化界面区氯离子扩散系数Nitz等。考虑到氯离子扩散系数的时变性，在组件定义模块中设置氯离子扩散系数时变模型Dt(t)。设置混凝土三维细观数值模型各组成相的初始氯离子浓度为0，即：C0＝0。考虑混凝土三维细观模型表面氯离子浓度的时变性，设置边界条件：C＝Cs(t)。
(4)设置氯离子传输方程
在模型开发器的稀物质传递模块中，设置氯离子在混凝土中的传输机理为“扩散”，控制方程采用Fick第二定律偏微分方程。
(5)设置氯离子传输维度
在模型开发器的稀物质传递模块中，根据不同的氯离子传输维度设置混凝土三维细观模型不同的封闭面和扩散面。其中，封闭面设置为无通量，扩散面浓度设置为时变的表面氯离子浓度Cs(t)。
(6)设置氯离子传递属性
在模型开发器的稀物质传递模块中，采用“层”的方法，分层设置传递属性。首先设置水泥浆体层，扩散系数定义为水泥浆体的时变扩散系数Dt(t)；随后，覆盖界面区层，界面区扩散系数定义为水泥浆体层扩散系数Dt(t)的Nitz倍；最后，覆盖集料层，该层的氯离子扩散系数设置为0。
(7)模型划分网格
在模型开发器的网格1模块中划分网格。仿真平台可自动划分网格，根据模型尺寸、计算精度等要求选择不同的网格大小。也可根据自身要求，调整网格的最大最小单元、曲率因子等进行网格优化。如图6所示，图6为模型网格划分。
(8)模型计算
选择瞬态研究模块，根据实际情况设置氯离子在混凝土三维细观数值模型中传输的时长、步长等数据，即可开展数值模型仿真计算。计算完成后，选择模型开发器中的结果模块，即可查看计算结果。如图7所示，图7为混凝土三维细观数值模型中氯离子传输仿真计算结果。
(9)数据分析
选择模型开发器中的结果模块，对氯离子传输计算结果进行分析。设定相应坐标的面平均取值，即可求得混凝土三维细观数值模型在任意时间、任意扩散深度处的氯离子浓度分布。
采用COMSOL有限元软件来对生成的混凝土三维细观模型进行其内氯离子传输数值仿真分析，操作简单、使用方便、功能多样、计算效率高。
本实施例提供的混凝土三维细观数值模型中氯离子传输仿真方法及系统，相比于物理试验，数值模拟不但成本低、周期短、效率高，且可以从细观角度探究多因素对混凝土中氯离子传输的影响。该技术方案能够在极大程度上克服混凝土三维模型中氯离子传输数值仿真研究中大多数将集料形状假定为规则的球形或椭球形从而导致与实际情况不符的问题。本实施例提供的方法从细观尺度出发，将混凝土视为由水泥浆体、集料及界面区共同组成的三相非均质复合材料，基于Monte-Carlo法、Fuller级配关系、德洛内三角网算法等手段，构建混凝土三维细观模型，根据有限元FEM原理，基于COMSOL稀物质传递模块模拟混凝土三维细观模型内部氯离子的传输过程，形成混凝土三维细观模型内氯离子传输仿真方法及系统。
以上所述实施例仅是为充分说明本发明而所举的较佳的实施例，本发明的保护范围不限于此。本技术领域的技术人员在本发明基础上所作的等同替代或变换，均在本发明的保护范围之内。本发明的保护范围以权利要求书为准。

Claims (6)

1.混凝土三维细观数值模型构建方法，其特征在于：包括以下步骤：
S1：构建三维多面体集料模型，根据集料粒径范围生成不同粒径大小的三维多面体集料；
S2：将三维多面体集料按照粒径di从大到小进行排序，并按照排列顺序将第i个多面体集料进行随机投放，通过Monte-Carlo法随机在L×M×H的混凝土三维细观模型线框中生成第i个多面体集料最小外接球的球心坐标；
重叠判别：判断第i个多面体集料是否会与之前已投放的i-1个集料重叠；采用多面体集料最小外接球球心距Li,j与多面体集料最小外接球的半径之和Lr-i,j二者的关系来进行重叠判别，记录并保存满足条件Li,j≥Lr-i,j的多面体集料几何信息；
其中，Li,j为相邻两多面体集料最小外接球球心距；Lr-i,j为相邻两多面体集料最小外接球半径之和；
S3：构建并生成第i个多面体集料的界面区，记录第i个多面体集料的界面区信息；
S4：直到全部多面体集料投放完毕，即完成混凝土三维细观数值模型的构建；
其中，步骤S2中的混凝土三维细观数值模型中多面体集料的随机投放，具体过程如下：
步骤S21：设置集料投放基本参数，所述集料投放基本参数包括集料投放边界距离ba，所述集料投放边界距离ba为净浆层厚度；
步骤S22：将生成的全部多面体集料按照粒径di从大到小降序排序，保证投放集料时按照粒径从大到小投放；
步骤S23：随机投放第i个多面体集料：基于Monte-Carlo法随机在L×M×H的混凝土三维数值模型线框中生成第i个多面体集料最小外接球的球心坐标，其横、纵、竖坐标的取值范围如下：
di/2+ba≤xi≤L-di/2-ba
di/2+ba≤yi≤M-di/2-ba
di/2+ba≤zi≤H-di/2-ba
其中，L、M、H分别为混凝土细观模型的长、宽、高；di为第i个集料粒径；
步骤S24：重叠判别：判断第i个多面体集料是否会与之前已投放的i-1个集料重叠；采用多面体集料最小外接球球心距Li,j与多面体集料最小外接球的半径之和Lr-i,j二者的关系来进行重叠判别；
Figure FDA0004105784250000021
Lr-i,j＝(di+dj)/2
式中：Li,j为相邻两多面体集料最小外接球球心距；Lr-i,j为相邻两多面体集料最小外接球半径之和；
判断计算得到的Li,j是否满足Li,j≥Lr-i,j，若不满足，返回执行步骤S23，若满足，则继续执行步骤S25；
步骤S25：保存并更新第i个多面体集料的几何信息，所述几何信息包括粒径、角点数量、角点坐标、体积、最小外接球球心坐标；
步骤S26：判断多面体集料是否投放完毕：判断i是否满足i≥na，na表示已生成的集料数量；若不满足，则令i＝i+1，返回循环执行步骤S23，若满足，则多面体集料投放完毕。
2.如权利要求1所述的混凝土三维细观数值模型构建方法，其特征在于：所述步骤S1中构建三维多面体集料模型按照以下步骤进行：
步骤(1)：确定三维多面体集料模型中的基本参数；
步骤(2)：随机生成第i个多面体集料粒径di；
步骤(3)：判断集料粒径di是否满足dmin≤di≤dmax，若满足，继续执行步骤(4)，若不满足，则返回重新执行步骤(2)；其中，集料最小粒径dmin、最大粒径dmax；
步骤(4)：生成粒径为di的多面体集料模型；
步骤(5)：计算第i个集料体积vai并按照以下方式判断是否满足集料体积占比Ca0：
判断第i个集料体积占比Cai是否满足Cai≥Ca0，若满足，则继续执行步骤(6)，若不满足，则返回重新执行步骤(2)，其中，Cai为第i个集料体积占比，Ca0表示集料体积占比；
步骤(6)：记录并保存第i个多面体集料信息；
步骤(7)：按照以下方式判断集料是否生成完毕：
判断集料累计体积va是否满足va≥vas，若不满足，令i＝i+1，返回步骤(2)；
若满足，则三维多面体集料生成完毕，其中，vas为混凝土集料总体积。
3.如权利要求2所述的混凝土三维细观数值模型构建方法，其特征在于：所述步骤(2)中的第i个多面体集料粒径di按照以下步骤生成：
根据Monte-Carlo算法生成区间[0,1]内的任意随机数Fi，并根据Fuller级配公式(1)求解得到符合Fuller级配曲线的多面体集料粒径di；
Figure FDA0004105784250000031
其中，di表示第i个集料的粒径。
4.如权利要求2所述的混凝土三维细观数值模型构建方法，其特征在于：所述步骤(4)的多面体集料模型按照以下方式构建：
以(0,0,0)为中心点坐标生成直径为di的球，在球内及球面上任意取点并根据所取的点采用德洛内三角网算法构建粒径为di的多面体集料。
5.如权利要求2所述的混凝土三维细观数值模型构建方法，其特征在于：所述步骤(5)中的集料体积占比Cai按照以下公式计算：
Cai＝vai/vsi＝vai/(π*di 3/6) (2)
式中：Cai为第i个集料体积占比；vai为第i个多面体集料体积；vsi为第i个集料对应的最小外接球的体积。
6.如权利要求1所述的混凝土三维细观数值模型构建方法，其特征在于：所述步骤S3按照以下步骤构建并生成第i个多面体集料的界面区：
步骤S31：设置界面区基本参数；
步骤S32：生成第i个多面体集料的界面区，多面体集料界面区厚度按以下公式计算；
titz＝δ*R (4)
式中：titz为界面区厚度，δ为界面区厚度倍数，R为集料最小外接球中心点距离多面体角点的距离；
将多面体集料外接球球心0与集料任一角点A相连，得到集料中心点与角点A的距离ROA，延长R0A的δ倍至A’点，延长距离即为该处界面区厚度titz；
其中，A’点坐标的计算表达式为：
Figure FDA0004105784250000041
式中：(x0,y0,z0)为多面体集料中心点的坐标；(xA,yA,zA)为集料角点A的坐标；
(xA’,yA’,zA’)为生成的界面区角点A’的坐标；δ为界面区厚度倍数；
步骤S33：记录第i个多面体集料界面区的信息；
步骤S34：判断集料界面区是否生成完毕：判断i是否满足i≥na，若不满足，则令i＝i+1，返回循环执行步骤S32，若满足，则全部多面体集料的界面区生成完毕；其中，na表示已生成的集料数量。
CN202210825040.4A 2022-07-14 2022-07-14 混凝土三维细观数值模型构建方法及氯离子传输仿真系统 Active CN115205486B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210825040.4A CN115205486B (zh) 2022-07-14 2022-07-14 混凝土三维细观数值模型构建方法及氯离子传输仿真系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210825040.4A CN115205486B (zh) 2022-07-14 2022-07-14 混凝土三维细观数值模型构建方法及氯离子传输仿真系统

Publications (2)

Publication Number Publication Date
CN115205486A CN115205486A (zh) 2022-10-18
CN115205486B true CN115205486B (zh) 2023-04-07

Family

ID=83579932

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210825040.4A Active CN115205486B (zh) 2022-07-14 2022-07-14 混凝土三维细观数值模型构建方法及氯离子传输仿真系统

Country Status (1)

Country Link
CN (1) CN115205486B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116884540B (zh) * 2023-06-15 2024-05-10 湖北工业大学 多相材料三维模型生成与连通性判断方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107590327A (zh) * 2017-08-31 2018-01-16 河海大学 一种水泥基材料颗粒三维随机投放的分层分区方法
CN110706352A (zh) * 2019-10-10 2020-01-17 重庆交通大学 基于多边形随机骨料的混凝土三相细观模型构建及内氯离子侵蚀数值模拟方法
CN110807242A (zh) * 2019-09-23 2020-02-18 重庆交通大学 考虑骨料分区域填充的钢筋混凝土细观数值模型构建方法
CN113962065A (zh) * 2021-09-27 2022-01-21 太原理工大学 一种含高体积分数骨料及气孔缺陷的混凝土三维细观模型
CN113987750A (zh) * 2021-09-27 2022-01-28 太原理工大学 一种含随机缺陷的全级配混凝土三维细观模型建模方法
WO2022121203A1 (zh) * 2020-12-08 2022-06-16 江苏科技大学 一种球壳表面三维裂纹扩展疲劳寿命的计算方法

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107688715A (zh) * 2017-09-12 2018-02-13 青岛理工大学 模拟温湿度耦合下混凝土氯离子浓度场数值方法
CN108932393A (zh) * 2018-07-13 2018-12-04 同济大学 一种基于骨料几何本征的混凝土二维细观模型建模方法
CN109211750B (zh) * 2018-08-10 2021-08-24 浙江工业大学 基于微、细观结构参数的混凝土氯离子扩散系数预测方法
CN110210178A (zh) * 2019-06-26 2019-09-06 西安理工大学 一种基于Python再生混凝土三维随机球形骨料模型的构建方法
CN113486530A (zh) * 2021-07-16 2021-10-08 南京工业职业技术大学 一种氯离子环境下钢筋混凝土使用状态分析的宏细观模型等效方法
CN114462124B (zh) * 2022-01-18 2023-05-26 东南大学 一种混凝土三维多相细观模型的建立与数值模拟方法
CN114510760B (zh) * 2022-01-19 2024-06-25 南京航空航天大学 一种基于离散元的大孔隙混凝土抗压试件模型构建方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107590327A (zh) * 2017-08-31 2018-01-16 河海大学 一种水泥基材料颗粒三维随机投放的分层分区方法
CN110807242A (zh) * 2019-09-23 2020-02-18 重庆交通大学 考虑骨料分区域填充的钢筋混凝土细观数值模型构建方法
CN110706352A (zh) * 2019-10-10 2020-01-17 重庆交通大学 基于多边形随机骨料的混凝土三相细观模型构建及内氯离子侵蚀数值模拟方法
WO2022121203A1 (zh) * 2020-12-08 2022-06-16 江苏科技大学 一种球壳表面三维裂纹扩展疲劳寿命的计算方法
CN113962065A (zh) * 2021-09-27 2022-01-21 太原理工大学 一种含高体积分数骨料及气孔缺陷的混凝土三维细观模型
CN113987750A (zh) * 2021-09-27 2022-01-28 太原理工大学 一种含随机缺陷的全级配混凝土三维细观模型建模方法

Also Published As

Publication number Publication date
CN115205486A (zh) 2022-10-18

Similar Documents

Publication Publication Date Title
CN110706352B (zh) 基于多边形随机骨料的混凝土三相细观模型构建及内氯离子侵蚀数值模拟方法
CN111539139B (zh) 一种颗粒随机分布的复合材料2d细观结构建模方法
CN107918706A (zh) 一种自适应的混凝土细观建模方法
CN115205486B (zh) 混凝土三维细观数值模型构建方法及氯离子传输仿真系统
CN112084694B (zh) 一种考虑非理想界面的颗粒增强复合材料微观结构的几何建模方法
CN110136782B (zh) 一种金属氧化物团簇结构模型的构造方法
CN108427861A (zh) 一种基于mpt工具包构建材料周期性多晶结构的方法
CN112052539A (zh) 一种基于Voronoi技术的钢筋混凝土三维细观快速并行建模方法
CN114186434B (zh) 基于逐级划分的沥青混合料三维细观结构模型构建方法
Sieger et al. A comprehensive comparison of shape deformation methods in evolutionary design optimization
CN109190144A (zh) 一种任意形状放射源辐射屏蔽计算仿真方法
CN115867430A (zh) 增材制造中的打印工艺制定方法及装置
CN110826128B (zh) 一种任意形状底面疏浚挖槽快速成形设计方法
CN116882236A (zh) 一种球体的六面体网格划分单元质量预测方法
CN114840891A (zh) 一种任意形状试件的三维混凝土细观模型建模方法
CN113486429B (zh) 一种基于插值算法的空间汇交结构自动化建造方法
CN114742944A (zh) 面向工业机器人路径规划的保守碰撞检测方法
Zhou et al. Research based on lee algorithm and genetic algorithm of the automatic external pipe routing of the aircraft engine
CN111159927B (zh) 基于体素矩阵的三维不规则形状颗粒投放的数值建模方法
Zheng et al. Parameterized implementation of a triple refined modeling framework for mesoscale concrete
Bova et al. Mesh generation/refinement using fractal concepts and iterated function systems
CN116738673A (zh) 基于真实骨料的三维混凝土细观几何模型的构建方法
Wang et al. Development of a hybrid parallelism Monte Carlo transport middleware on mesh geometry
CN108763677A (zh) 一种弥散燃料中子输运模拟的方法、系统及服务器
CN109214726B (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