CN110147601A - 一种预测混凝土中有害物质扩散的耦合模拟方法 - Google Patents
一种预测混凝土中有害物质扩散的耦合模拟方法 Download PDFInfo
- Publication number
- CN110147601A CN110147601A CN201910398384.XA CN201910398384A CN110147601A CN 110147601 A CN110147601 A CN 110147601A CN 201910398384 A CN201910398384 A CN 201910398384A CN 110147601 A CN110147601 A CN 110147601A
- Authority
- CN
- China
- Prior art keywords
- aggregate
- model
- concrete
- diffusion
- harmful substance
- 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.)
- Pending
Links
- 239000004567 concrete Substances 0.000 title claims abstract description 106
- 238000009792 diffusion process Methods 0.000 title claims abstract description 94
- 239000000126 substance Substances 0.000 title claims abstract description 75
- 238000000034 method Methods 0.000 title claims abstract description 70
- 238000004088 simulation Methods 0.000 title claims abstract description 55
- 239000000463 material Substances 0.000 claims abstract description 17
- 239000004570 mortar (masonry) Substances 0.000 claims abstract description 15
- 230000007704 transition Effects 0.000 claims abstract description 15
- 238000010168 coupling process Methods 0.000 claims abstract description 14
- 230000008878 coupling Effects 0.000 claims abstract description 13
- 238000005859 coupling reaction Methods 0.000 claims abstract description 13
- 230000009471 action Effects 0.000 claims description 15
- 230000000694 effects Effects 0.000 claims description 13
- 230000006378 damage Effects 0.000 claims description 9
- 238000002474 experimental method Methods 0.000 claims description 9
- 230000008859 change Effects 0.000 claims description 8
- 230000004907 flux Effects 0.000 claims description 3
- 230000001771 impaired effect Effects 0.000 claims 1
- 238000012546 transfer Methods 0.000 description 8
- 230000005540 biological transmission Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- VEXZGXHMUGYJMC-UHFFFAOYSA-M Chloride anion Chemical compound [Cl-] VEXZGXHMUGYJMC-UHFFFAOYSA-M 0.000 description 4
- 230000006835 compression Effects 0.000 description 4
- 238000007906 compression Methods 0.000 description 4
- 238000013461 design Methods 0.000 description 4
- 239000002689 soil Substances 0.000 description 4
- 238000005452 bending Methods 0.000 description 3
- 230000035515 penetration Effects 0.000 description 3
- 230000000704 physical effect Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 210000000988 bone and bone Anatomy 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 238000004215 lattice model Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 239000011241 protective layer Substances 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 230000015271 coagulation Effects 0.000 description 1
- 238000005345 coagulation Methods 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 230000008595 infiltration Effects 0.000 description 1
- 238000001764 infiltration Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000010410 layer Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000002250 progressing effect Effects 0.000 description 1
- 239000011150 reinforced concrete Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/38—Concrete; Lime; Mortar; Gypsum; Bricks; Ceramics; Glass
- G01N33/383—Concrete or cement
-
- 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]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Computer Hardware Design (AREA)
- Biochemistry (AREA)
- Medicinal Chemistry (AREA)
- Food Science & Technology (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Analytical Chemistry (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Ceramic Engineering (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种预测混凝土中有害物质扩散的耦合模拟方法，包括以下步骤：生成由骨料、砂浆和界面过渡区组成的随机骨料分布模型；对随机骨料分布模型定义边界条件、施加荷载大小及混凝土模型材料参数；对随机骨料分布模型进行网格划分；对随机骨料分布模型定义力学本构模型，模拟混凝土模型的力学场；定义耦合应力‑物质扩散方程，引入生成的应力参数，模拟有害物质在混凝土中的渗透扩散过程。本发明提供的方法能够适用于混凝土结构在真实受力状态下，混凝土保护层中有害物质侵蚀扩散情况的预测，从而更加准确的预测混凝土结构的使用寿命。
Description
技术领域
本发明涉及钢筋混凝土结构保护层的耐久性设计技术领域，特别涉及一种预测混凝土中有害物质扩散的耦合模拟方法。
背景技术
混凝土在受荷载作用下的传输性能已经成为混凝土领域广泛研究的重要问题之一。混凝土是一种不均质的多相复合材料，在荷载作用下其内部的应力分布极其不均匀，由此可知，各点的扩散系数之间的差异也很明显，以往的研究忽略了材料的不均质导致混凝土内部应力值的变化。
同时，大部分混凝土结构服役过程是在低荷载作用下，此时的混凝土处于弹性状态，裂缝不会显著增长。但在某些特殊部位，混凝土的受力较大，内部某些区域因发生塑性变形而出现损伤，裂缝发展明显，这些现象必然会导致混凝土局部渗透扩散性质的显著改变。
因此，有必要对混凝土结构的受力导致的传输性能进行细致的计算、预测和保护层设计。现有的采用模拟方法研究混凝土中物质传输性质主要有等效单相介质法、有限元法及Lattice模型。
其中，等效介质法将混凝土视为均匀的介质，通过控制扩散系数来表示应力对传输性质的影响。有限元法的基本原理是通过对物质传输偏微分方程的推导和求解，编写对应的计算机程序来实现物质传输过程。Lattice模型是用单元杆件代替网格的有限元求解方法。
但是，国内外大部分学者关于物质传输的模拟主要是将混凝土视为单一的匀质材料，荷载的大小通过改变渗透系数的大小来实现耐久性的预测。采用上述方法无法真实反映混凝土的细观尺度上的应力分布不均匀现象，也无法进一步研究高荷载作用下，混凝土产生开裂损伤下的物质传输现象。
因而现有技术还有待改进和提高。
发明内容
鉴于上述现有技术的不足之处，本发明的目的在于提供一种预测混凝土中有害物质扩散的耦合模拟方法，能够适用于混凝土结构在真实受力状态下，混凝土保护层中有害物质侵蚀扩散情况的预测。
为了达到上述目的，本发明采取了以下技术方案：
一种预测混凝土中有害物质扩散的耦合模拟方法，包括以下步骤：
A、根据实际工程结构对象确定随机骨料分布模型的基本模型信息；
B、根据所述基本模型信息进行随机骨料填充，生成随机骨料分布模型；
C、对随机骨料分布模型定义力学本构模型，模拟混凝土模型的力学场；
D、定义耦合应力-物质扩散方程，引入步骤C生成的应力参数，模拟有害物质在混凝土中的渗透扩散过程。
进一步来说，所述的预测混凝土中有害物质扩散的耦合模拟方法中，在所述步骤B之后，还包括步骤：
B1'、对随机骨料分布模型定义边界条件、施加荷载大小及混凝土模型材料参数；
B2'、对随机骨料分布模型进行网格划分。
所述的预测混凝土中有害物质扩散的耦合模拟方法中，所述基本模型信息包括骨料随机分布模型的长度a及宽度b、骨料的半径r、骨料的形状、骨料级配和界面过渡区的厚度r1。
所述的预测混凝土中有害物质扩散的耦合模拟方法中，所述步骤B具体包括：
B1、随机生成一个包含骨料的形心O(x，y)及半径r的骨料信息，并进行随机投放；
B2、判断条件x-r>0、x-r<a和y-r<b是否成立；
B3、若判断条件成立，则存储骨料信息，否则删除的骨料信息；
B4、不断重复上述步骤并存储新的骨料信息，同时判断是否已达到需要的骨料投放量；
B5、若已达到需要的骨料投放量，则停止生成新的骨料信息，并绘制已存储骨料的形状，生成骨料随机分布模型。
所述的预测混凝土中有害物质扩散的耦合模拟方法中，所述步骤B4还包括：
B41、判断新存储的骨料与已存储的骨料是否存在重合部分，若新存储的骨料与已存储的骨料存在重合部分，则删除新存储的骨料，重新确定投放骨料信息。
所述的预测混凝土中有害物质扩散的耦合模拟方法中，所述步骤B1'具体包括：
B11'、设置物质扩散场的边界条件为左右任意一条边界设置表面浓度C，下边界为对称界面，其他边界及骨料边界为无通量；
B12'、对随机骨料分布模型的上边界施加荷载，同时在上边界上与荷载垂直方向的位移为0，下边界为对称界面，左右两边界为自由边界；
B13'、通过实验测定界面过渡区和砂浆在轴向拉压荷载作用下及无荷载作用下的扩散系数。
所述的预测混凝土中有害物质扩散的耦合模拟方法中，所述步骤B2'具体包括：
将随机骨料分布模型划分为自由三角形网格。
所述的预测混凝土中有害物质扩散的耦合模拟方法中，所述步骤C中随机骨料分布模型定义的力学本构模型采用Ottosen四参数强度准则。
所述的预测混凝土中有害物质扩散的耦合模拟方法中，所述耦合应力-物质扩散方程为：
其中，Dc为施加荷载后的扩散系数，为裂缝中的扩散系数，该系数与等效裂缝宽度有关；D0为无荷载作用下的扩散系数，系数α通过实验测定，作用是控制扩散系数变化的斜率；l为单位长度，lεp为永久性变形，dl(ε-εp)为可恢复的变形，E为初始杨氏模量，为受损后的杨氏模量，d为损伤系数，其中ε＝u/E，u为施加的荷载值。D(θ)为施加荷载后的扩散系数，θ为应力水平，σ为应力。
所述的预测混凝土中有害物质扩散的耦合模拟方法中，所述D(θ)关系式的确定方法为：
对砂浆和界面过渡区进行轴向拉压实验，确定不同荷载水平的扩散系数；
通过拟合的方法得出砂浆和界面过渡区在不同荷载作用下的D(θ)关系式。
相较于现有技术，本发明提供的一种预测混凝土中有害物质扩散的耦合模拟方法，包括以下步骤：生成由骨料、砂浆和界面过渡区组成的随机骨料分布模型；对随机骨料分布模型定义边界条件、施加荷载大小及混凝土模型材料参数；对随机骨料分布模型进行网格划分；对随机骨料分布模型定义力学本构模型，模拟混凝土模型的力学场；定义耦合应力-物质扩散方程，引入生成的应力参数，模拟有害物质在混凝土中的渗透扩散过程。本发明提供的方法能够适用于混凝土结构在真实受力状态下，混凝土保护层中有害物质侵蚀扩散情况的预测，从而更加准确的预测混凝土结构的使用寿命。
附图说明
图1为本发明提供的一种预测混凝土中有害物质扩散的耦合模拟方法的流程图；
图2为本发明提供的一种预测混凝土中有害物质扩散的耦合模拟方法中流程图二；
图3为本发明提供的一种预测混凝土中有害物质扩散的耦合模拟方法中步骤B的流程图；
图4为本发明提供的一种预测混凝土中有害物质扩散的耦合模拟方法中骨料投放方法的流程图；
图5为本发明提供的一种预测混凝土中有害物质扩散的耦合模拟方法中步骤B1'的流程图；
图6为本发明提供的一种预测混凝土中有害物质扩散的耦合模拟方法中实施例一的随机骨料分布模型的结构示意图；
图7为本发明提供的一种预测混凝土中有害物质扩散的耦合模拟方法中实施例一的随机骨料分布模型的边界条件和各相的扩散系数示意图；
图8为本发明提供的一种预测混凝土中有害物质扩散的耦合模拟方法中实施例一的随机骨料分布模型的自由三角形网格图；
图9为本发明提供的一种预测混凝土中有害物质扩散的耦合模拟方法中实施例一的随机骨料分布模型施加不同荷载后的应力值分布图；
图10为本发明提供的一种预测混凝土中有害物质扩散的耦合模拟方法中实施例一的随机骨料分布模型施加不同荷载后的损伤系数图；
图11为本发明提供的一种预测混凝土中有害物质扩散的耦合模拟方法中实施例一的随机骨料分布模型施加轴向拉伸荷载后氯离子浓度随深度变化图；
图12为本发明提供的一种预测混凝土中有害物质扩散的耦合模拟方法中实施例一的随机骨料分布模型施加轴向压缩荷载后氯离子浓度随深度变化图；
图13为本发明提供的一种预测混凝土中有害物质扩散的耦合模拟方法中实施例二的施加集中荷载混凝土梁的受力示意图；
图14为本发明提供的一种预测混凝土中有害物质扩散的耦合模拟方法中实施例二的施加集中荷载混凝土梁在弯曲荷载作用下氯离子浓度随深度变化图。
具体实施方式
本发明提供一种预测混凝土中有害物质扩散的耦合模拟方法，能够适用于混凝土结构在真实受力状态下，混凝土保护层中有害物质侵蚀扩散情况的预测。为使本发明的目的、技术方案及效果更加清楚、明确，以下参照附图并举实施例对本发明进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
在应用混凝土的实际工程结构中，混凝土可能承受轴向拉伸荷载、压缩荷载或弯曲荷载，更多地时候是承受多种荷载类型相互作用的复杂受力状态。在实验室能测出简单、且荷载比值较小时的物质扩散系数，如轴向拉伸和压缩荷载作用，而通过实验室的测得的扩散系数验证实际工程中混凝土受力状态下的渗透性质往往不现实。
本发明具体是一种通过不同级别轴向拉压的简单荷载与扩散系数的关系，将结构中真实存在的复杂应力场计算出各点真实的扩散系数，并计算有害物质在此受力状态下的真实扩散行为的耦合应力-物质传输模拟方法，目的是解决现有数值方法中存在的以上问题，提供一种准确地预测结构在真实荷载作用下保护层中有害物质扩散行为的计算方法，更加准确的预测结构的使用寿命。
请参阅图1，为本发明提供的一种预测混凝土中有害物质扩散的耦合模拟方法的流程图，所述预测混凝土中有害物质扩散的耦合模拟方法包括以下步骤：
A、根据实际工程结构对象确定随机骨料分布模型的基本模型信息；
B、根据所述基本模型信息进行随机骨料填充，生成随机骨料分布模型(即生成由骨料、砂浆和界面过渡区(Interface Transition Zone，简称ITZ)组成的随机骨料分布模型)；
C、对随机骨料分布模型定义力学本构模型，模拟混凝土模型的力学场；
D、定义耦合应力-物质扩散方程，引入步骤C生成的应力参数，模拟有害物质在混凝土中的渗透扩散过程。
并且，请参阅图2，所述的预测混凝土中有害物质扩散的耦合模拟方法中，在步骤B之后，还包括步骤：
B1'、对随机骨料分布模型定义边界条件、施加荷载大小及混凝土模型材料参数；
B2'、对随机骨料分布模型进行网格划分。
其中，所述随机骨料分布模型包括骨料随机分布模型的长度a及宽度b(尺寸)、骨料的半径r、骨料的形状、骨料级配和界面过渡区的厚度r1等基本模型信息。
通过执行步骤A-步骤D、步骤B1'、步骤B2'后得到的耦合应力-物质扩散方程，可以更加真实反映混凝土的细观尺度上的应力分布不均匀现象，也可以进一步研究高荷载作用下，混凝土产生开裂损伤下的物质传输现象。以下分别对步骤A-步骤D、步骤B1'、步骤B2'的方案分别进行详细说明。
对于步骤B，请参阅图3，其具体包括：
B1、随机生成一个包含骨料的形心O(x，y)及半径r的骨料信息，并进行随机投放；
B2、判断条件x-r>0、x-r<a和y-r<b是否成立；
B3、若判断条件成立，则存储骨料信息，否则删除的骨料信息；
B4、不断重复上述步骤并存储新的骨料信息，同时判断是否已达到需要的骨料投放量；
B5、若已达到需要的骨料投放量，则停止生成新的骨料信息，并绘制已存储骨料的形状，生成骨料随机分布模型。
在实际操作中，请结合图4，可以利用MATLAB软件中rand函数生成一个骨料的形心O(x，y)及半径r的信息，进行随机投放。然后，判断条件x-r>0和x-r<a和y-r<b是否成立；若条件成立，则代表该骨料没有超出模型尺寸范围，可以存储此次的骨料信息；若条件不成立，则代表该骨料超出了模型尺寸范围，不符合要求，需要抛弃此次生成的骨料信息，重新生成新的骨料信息。
其次，在步骤B4中还包括步骤：B41、判断新存储的骨料与已存储的骨料是否存在重合部分，若新存储的骨料与已存储的骨料存在重合部分，则删除新存储的骨料，重新确定投放骨料信息。
之后，对于判断是否已达到需要的骨料投放量，可根据模型信息，当π×r2≥a×b时，则认为达到骨料预期占比。
最后，通过存储的骨料信息在MATLAB中绘制多边形，同时利用骨料的半径、形心的信息绘制ITZ，ITZ半径为r+r1。
进一步的，请参阅图5，所述步骤B1'具体包括：
B11'、设置物质扩散场的边界条件为左右任意一条边界设置表面浓度C，下边界为对称界面，其他边界及骨料边界为无通量；
B12'、对随机骨料分布模型的上边界施加荷载，同时在上边界上与荷载垂直方向的位移为0，下边界为对称界面，左右两边界为自由边界；
B13'、通过实验测定界面过渡区和砂浆在轴向拉压荷载作用下及无荷载作用下的扩散系数(分别为DC和D0)。
更进一步的，所述步骤B2'具体包括：
对随机骨料分布模型进行网格划分，并划分为自由三角形网格。通过用等边三角形进行离散，既能照顾到计算对象的细微几何特征，又能照顾到仅需稀疏单元网格之处。
更进一步的，所述步骤C中随机骨料分布模型定义的力学本构模型采用Ottosen四参数强度准则。
Ottosen四参数强度模型是由两个混凝土单轴强度，两个典型的双轴和三轴强度来确定的，其比较全面反映混凝土破坏特征。化简后的Ottosen强度准则表达式的参数只需确定单轴抗压强度就可以求解破坏面方程，四参数即a、b两个常数，尺寸因子K1和形状因子K2
Ottosen强度准则适用于均质材料纯拉伸、纯压缩、纯剪切情况下的强度试验，同时也适用于非均质材料在复杂应力状态下破坏形态。其适用范围广，精度高，采用该强度准则可以使模型的应力状态更接近混凝土结构真实的应力状态。步骤C中力学场的数值计算可提供塑性应变εp和有效应力σ运用在步骤D中。
具体来说，步骤D中载荷作用下混凝土中耦合应力-物质扩散方程为：
其中，Dc为施加荷载后的扩散系数，为裂缝中的扩散系数，该系数与等效裂缝宽度有关；D0为无荷载作用下的扩散系数，系数α通过实验测定，作用是控制扩散系数变化的斜率；l为单位长度，lεp为永久性变形，dl(ε-εp)为可恢复的变形，E为初始杨氏模量，为受损后的杨氏模量，d为损伤系数(0≤d≤1)，其中ε＝u/E，u为施加的荷载值(Pa)，σ表示为应力。
D(θ)为施加轴向荷载后的扩散系数，θ为应力水平(数值模拟计算的有效应力与混凝土极限抗拉或抗压强度之比。当压缩荷载作用时，此时θ为负值)。D(θ)的关系通过实验室确定，方法为：
(1)对砂浆和ITZ进行轴向拉压实验，确定不同荷载水平的扩散系数，如施加荷载水平为-0.4、-0.2、0、0.2、0.4、0.6；
(2)通过拟合的方法得出砂浆和ITZ在不同荷载作用下的D(θ)关系式；
因此，耦合应力-物质传输的Fick第二定律偏微分方程可以写成：
本发明中的模型采用多物理场耦合方法研究荷载作用下混凝土内部应力分布影响有害物质的扩散效应。荷载产生的不均匀应力是独立建模的，在物质传输方程中允许每次运算过程都引入应力参数，从而实现应力与传输的耦合。
此外，为了更好的理解本发明，本发明提供一种具体实施例一如下：
请参阅图6，为本实施例中生成的尺寸为50×50mm的随机骨料分布模型，其中圆形骨料r＝2～5mm，r1＝50μm，骨料占比为50％。首先定义表1中各相的物理性质参数，进行力学场的边界条件设置，如图7所示。然后通过有限元软件自动划分的自由三角形网格，如图8所示。通过采用步骤D中Ottosen准则生成的应力分布如图9所示。假设物质传输过程的初始条件及边界条件为：
C(t＝0)＝0 (6)
C(x＝50mm)＝2823mol/m3 (7)
因此本实施例中砂浆和ITZ的扩散系数Dc(θ)、DITZ(θ)的方程式为：
系数α＝400，通过公式(1)～(5)可以求解混凝土模型在荷载作用下的物质扩散情况，即如图10-图12所示。
表1材料的物理性质参数
其中，由于ITZ密度难于准确测量，在这里取1000kg/m3。
同时，本发明还提供一种具体实施例二如下：
本实施例描述了弯曲荷载作用下混凝土梁的扩散问题。如图13所示，为本实施例中施加集中荷载的混凝土梁的受力示意图，具体有一个集中力F作用在尺寸为40*80mm的混凝土梁的中心，有害物质从上下两表面扩散进入混凝土梁中。使用对称性可减少计算时间并降低复杂性。取该模型中间截面作为分界线(虚线表示对称轴)进行计算。并且，本实施例的材料的物理性质参数与实施例一相同。如图14所示，最后可得到弯曲荷载分别为(a)F＝0、(b)F＝10kN、(c)F＝20kN、(d)F＝30kN、(e)F＝34kN以及(f)F＝38KN作用下氯离子浓度随深度变化图。
因此，本发明具有优点如下：
1)本发明中的方法可以根据工程实际混凝土结构的受力状况及环境条件，通过简单的轴向拉压测试即可对模型的扩散系数及力学参数进行设置，在从而对混凝土结构的保护层厚度进行合理的耐久性设计；
2)本发明中的方法可以通过实验室测定混凝土受拉压时扩散系数的变化，表征不同荷载水平(低和高荷载)作用下及不同受荷类型下各部位的应力分布差异，从而研究物质的渗透扩散行为更科学合理；
3)本发明中的方法充分考虑混凝土材料中砂浆和ITZ的力学性能及渗透性质差异，更具有实际工程指导意义；
综上所述，本发明提供的一种预测混凝土中有害物质扩散的耦合模拟方法，包括以下步骤：生成由骨料、砂浆和界面过渡区组成的随机骨料分布模型；对随机骨料分布模型定义边界条件、施加荷载大小及混凝土模型材料参数；对随机骨料分布模型进行网格划分；对随机骨料分布模型定义力学本构模型，模拟混凝土模型的力学场；定义耦合应力-物质扩散方程，引入生成的应力参数，模拟有害物质在混凝土中的渗透扩散过程。本发明提供的方法能够适用于混凝土结构在真实受力状态下，混凝土保护层中有害物质侵蚀扩散情况的预测，从而更加准确的预测混凝土结构的使用寿命。
可以理解的是，对本领域普通技术人员来说，可以根据本发明的技术方案及其发明构思加以等同替换或改变，而所有这些改变或替换都应属于本发明所附的权利要求的保护范围。
Claims (10)
1.一种预测混凝土中有害物质扩散的耦合模拟方法，其特征在于，包括以下步骤：
A、根据实际工程结构对象确定随机骨料分布模型的基本模型信息；
B、根据所述基本模型信息进行随机骨料填充，生成随机骨料分布模型；
C、对随机骨料分布模型定义力学本构模型，模拟混凝土模型的力学场；
D、定义耦合应力-物质扩散方程，引入步骤C生成的应力参数，模拟有害物质在混凝土中的渗透扩散过程。
2.根据权利要求1所述的一种预测混凝土中有害物质扩散的耦合模拟方法，其特征在于，在所述步骤B之后，还包括步骤：
B1'、对随机骨料分布模型定义边界条件、施加荷载大小及混凝土模型材料参数；
B2'、对随机骨料分布模型进行网格划分。
3.根据权利要求1所述的一种预测混凝土中有害物质扩散的耦合模拟方法，其特征在于，所述基本模型信息包括骨料随机分布模型的长度a及宽度b、骨料的半径r、骨料的形状、骨料级配和界面过渡区的厚度r1。
4.根据权利要求2所述的一种预测混凝土中有害物质扩散的耦合模拟方法，其特征在于，所述步骤B具体包括：
B1、随机生成一个包含骨料的形心O(x，y)及半径r的骨料信息，并进行随机投放；
B2、判断条件x-r>0、x-r<a和y-r<b是否成立；
B3、若判断条件成立，则存储骨料信息，否则删除的骨料信息；
B4、不断重复上述步骤并存储新的骨料信息，同时判断是否已达到需要的骨料投放量；
B5、若已达到需要的骨料投放量，则停止生成新的骨料信息，并绘制已存储骨料的形状，生成骨料随机分布模型。
5.根据权利要求4所述的一种预测混凝土中有害物质扩散的耦合模拟方法，其特征在于，所述步骤B4还包括：
B41、判断新存储的骨料与已存储的骨料是否存在重合部分，若新存储的骨料与已存储的骨料存在重合部分，则删除新存储的骨料，重新确定投放骨料信息。
6.根据权利要求2所述的一种预测混凝土中有害物质扩散的耦合模拟方法，其特征在于，所述步骤B1'具体包括：
B11'、设置物质扩散场的边界条件为左右任意一条边界设置表面浓度C，下边界为对称界面，其他边界及骨料边界为无通量；
B12'、对随机骨料分布模型的上边界施加荷载，同时在上边界上与荷载垂直方向的位移为0，下边界为对称界面，左右两边界为自由边界；
B13'、通过实验测定界面过渡区和砂浆在轴向拉压荷载作用下及无荷载作用下的扩散系数。
7.根据权利要求2所述的一种预测混凝土中有害物质扩散的耦合模拟方法，其特征在于，所述步骤B2'具体包括：
将随机骨料分布模型划分为自由三角形网格。
8.根据权利要求1所述的一种预测混凝土中有害物质扩散的耦合模拟方法，其特征在于，所述步骤C中随机骨料分布模型定义的力学本构模型采用Ottosen四参数强度准则。
9.根据权利要求8所述的一种预测混凝土中有害物质扩散的耦合模拟方法，其特征在于，所述耦合应力-物质扩散方程为：
其中，Dc为施加荷载后的扩散系数，为裂缝中的扩散系数，该系数与等效裂缝宽度有关；D0为无荷载作用下的扩散系数，系数α通过实验测定，作用是控制扩散系数变化的斜率；l为单位长度，lεp为永久性变形，dl(ε-εp)为可恢复的变形，E为初始杨氏模量，为受损后的杨氏模量，d为损伤系数，其中ε＝u/E，u为施加的荷载值。D(θ)为施加轴向荷载后的扩散系数，θ为应力水平，σ为应力。
10.根据权利要求9所述的一种预测混凝土中有害物质扩散的耦合模拟方法，其特征在于，所述D(θ)关系式的确定方法为：
对砂浆和界面过渡区进行轴向拉压实验，确定不同荷载水平的扩散系数；
通过拟合的方法得出砂浆和界面过渡区在不同荷载作用下的D(θ)关系式。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910398384.XA CN110147601A (zh) | 2019-05-14 | 2019-05-14 | 一种预测混凝土中有害物质扩散的耦合模拟方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910398384.XA CN110147601A (zh) | 2019-05-14 | 2019-05-14 | 一种预测混凝土中有害物质扩散的耦合模拟方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110147601A true CN110147601A (zh) | 2019-08-20 |
Family
ID=67595343
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910398384.XA Pending CN110147601A (zh) | 2019-05-14 | 2019-05-14 | 一种预测混凝土中有害物质扩散的耦合模拟方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110147601A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117077501A (zh) * | 2023-10-17 | 2023-11-17 | 中南大学 | 一种非均质钢渣混凝土渗透系数计算方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005071036A (ja) * | 2003-08-22 | 2005-03-17 | Ohbayashi Corp | コンクリート内の物質移行を解析するシステムおよび方法 |
| CN108229093A (zh) * | 2018-01-23 | 2018-06-29 | 广西大学 | 饱和再生混凝土氯离子扩散系数多尺度预测模型的构建方法 |
| CN108304689A (zh) * | 2018-01-23 | 2018-07-20 | 广西大学 | 预应力混凝土氯离子扩散系数多尺度预测模型的构建方法 |
| CN109408971A (zh) * | 2018-10-29 | 2019-03-01 | 天津大学 | 一种新型多参数影响下早期混凝土性态变化研究方法 |
-
2019
- 2019-05-14 CN CN201910398384.XA patent/CN110147601A/zh active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005071036A (ja) * | 2003-08-22 | 2005-03-17 | Ohbayashi Corp | コンクリート内の物質移行を解析するシステムおよび方法 |
| CN108229093A (zh) * | 2018-01-23 | 2018-06-29 | 广西大学 | 饱和再生混凝土氯离子扩散系数多尺度预测模型的构建方法 |
| CN108304689A (zh) * | 2018-01-23 | 2018-07-20 | 广西大学 | 预应力混凝土氯离子扩散系数多尺度预测模型的构建方法 |
| CN109408971A (zh) * | 2018-10-29 | 2019-03-01 | 天津大学 | 一种新型多参数影响下早期混凝土性态变化研究方法 |
Non-Patent Citations (4)
| Title |
|---|
| 徐华祥: "大坝混凝土二维随机骨料的数值模拟", 《中国水运》 * |
| 李彤梅: "基于ottosen准则的钢筋混凝土结构非线性有限元方法", 《第18届全国结构工程学术会议论文集第Ⅰ册》 * |
| 杜修力等: "压缩荷载作用下混凝土中氯离子扩散行为细观模拟", 《建筑材料学报》 * |
| 邓德华 等: "《土木工程材料》" * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117077501A (zh) * | 2023-10-17 | 2023-11-17 | 中南大学 | 一种非均质钢渣混凝土渗透系数计算方法 |
| CN117077501B (zh) * | 2023-10-17 | 2024-01-09 | 中南大学 | 一种非均质钢渣混凝土渗透系数计算方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Grassl et al. | On a 2D hydro-mechanical lattice approach for modelling hydraulic fracture | |
| Zahra et al. | Characterisation of cementitious polymer mortar–Auxetic foam composites | |
| Zhang et al. | A discrete-continuum coupled finite element modelling approach for fibre reinforced concrete | |
| Smith et al. | Discrete modeling of ultra-high-performance concrete with application to projectile penetration | |
| Yin et al. | Tensile fracture simulation of random heterogeneous asphalt mixture with cohesive crack model | |
| Haeri et al. | XFEM and experimental simulation of failure mechanism of non-persistent joints in mortar under compression | |
| Angiolilli et al. | Lattice discrete particle model for the simulation of irregular stone masonry | |
| Zhou et al. | Peridynamic micro-elastoplastic constitutive model and its application in the failure analysis of rock masses | |
| Zhou et al. | 3D numerical study on the growth and coalescence of pre-existing flaws in rocklike materials subjected to uniaxial compression | |
| Feng et al. | Study on the non-linear deformation and failure characteristics of EPS concrete based on CT-scanned structure modelling and cloud computing | |
| Daoud et al. | 2D mesoscopic modelling of bar–concrete bond | |
| Han et al. | Effects of graded concrete on compressive strengths | |
| Li et al. | A mechanical-diffusive peridynamics coupling model for meso-scale simulation of chloride penetration in concrete under loadings | |
| Chen et al. | Centrifuge modeling testing and multiscale analysis of cemented sand and gravel (CSG) dams | |
| Qin et al. | Numerical analysis of fracturing behavior in fully-graded concrete with oversized aggregates from mesoscopic perspective | |
| Li et al. | Further development of the distinct lattice spring model for quasi-brittle crack propagation in concrete and its application in underground engineering | |
| Lin et al. | An effective numerical model for reinforced concrete beams strengthened with high performance fibre reinforced cementitious composites | |
| Goh et al. | Compression test and finite element analysis of foamed concrete cube | |
| Luo et al. | Statistical analysis of mesoscopic concrete with random elastic modulus | |
| Silva et al. | Numerical analysis of the in-plane behaviour of three-leaf stone masonry panels consolidated with grout injection | |
| Huang et al. | Meso-scale analysis of failure characteristics and mechanical properties of lightweight aggregate concrete (LWAC) with different aggregate volume fractions and shapes under axial tension | |
| Okabe et al. | Micromechanics on the rate-dependent fracture of discontinuous fiber-reinforced plastics | |
| Wu et al. | Numerical modeling on micro-to-macro evolution of crack network for concrete materials | |
| Zhang et al. | Failure evolution and fiber toughing mechanism of ultra-high performance concrete under uniaxial compression | |
| CN110147601A (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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190820 |
|
| RJ01 | Rejection of invention patent application after publication |