一种混凝土氯离子扩散与渗透作用的模拟方法
【专利摘要】本发明公开了一种混凝土氯离子扩散与渗透作用的模拟方法。离散化用实体单元建立空间模型；前处理将实体单元进行分类；判断识别内部空腔并修正实体单元各面的氯离子浓度；计算氯离子浓度差，再施加荷载计算应力差值然后修正扩散系数；建立三维非均匀单元的元胞自动机模型，得到进化系数；计算得到时间差值内每个实体单元各面的氯离子浓度；迭代计算完成对氯离子在混凝土中扩散与渗透作用的模拟。本发明使扩散模型应用到了结构层次上，也不再局限于规则单元；模型同时考虑浓度梯度引起的扩散作用和压力梯度引起的渗透作用，从结构层面上模拟氯离子对混凝土的侵蚀效应，可有针对性地进行老化结构的维护和修复，使得选择更加合理。
【专利说明】一种混凝土氯离子扩散与渗透作用的模拟方法
【技术领域】
[0001]本发明属于一种混凝土的模拟方法，具体涉及桥梁工程【技术领域】的一种混凝土氯离子扩散与渗透作用的模拟方法。
【背景技术】
[0002]对于混凝土桥梁结构，由于材料力学性能的逐渐劣化，其承载力将随时间而降低，导致桥梁结构无法承受原设计荷载，因此我们必须考虑结构性能的时变性。既然混凝土桥梁结构所施加的荷载和结构本身的抗力都是随时间变化的，可以先对材料进行耐久性及劣化规律分析，而后对结构在实际使用寿命期内的可靠度进行估算。而要进行以上分析，最重要的就是对环境中的有害物质在混凝土结构中的侵蚀进行研究，实现模拟预测。以便在进行桥梁结构设计中增加概念设计或者制定结构的维护和修复方案以达到桥梁结构的设计使用寿命。在材料耐久性及劣化的影响因素中，环境中的有害物质对混凝土性能的影响最大，尤其是氯离子侵入可以借助于毛细管孔隙吸力或者裂缝的渗透扩散来发生。现有的模拟氯离子对混凝土的侵蚀基本上都是从Fick第二定律出发来考虑氯离子在混凝土内的扩散，虽然不同的研究者考虑不同的因素建立了各种各样的模型，但是模型往往含有高阶微分，形式相当复杂，应用上大多限于对混凝土材料和构件的应用研究，鲜有在整体结构上的应用。
【发明内容】
[0003]为了解决氯离子扩散模型在结构层面上的应用难题，本发明的目的是提供一种混凝土氯离子扩散与渗透作用的模拟方法，空间非规则单元的元胞自动机模型模拟，其实施过程中涉及了氯离子的扩散作用、氯离子的渗透作用、扩散系数、渗透系数、扩散与渗透的协同作用、空间非规则单元元胞自动机模型、元胞自动机的进化参数。
[0004]本发明采用的技术方案包括如下步骤:
[0005]I)将桥梁结构离散化，用实体单元建立桥梁结构的空间模型；
[0006]2)进行前处理，将实体单元根据暴露在空气中表面数目进行分类；
[0007]3)对暴露在空气中的实体单元进行进一步的处理，判断识别该实体单元是否位于箱梁结构的内部空腔壁上，然后修正该实体单元进行元胞自动机计算时的环境氯离子浓度；
[0008]4)计算相邻实体单元的氯离子浓度差，再根据桥梁结构设计规范对结构施加荷载计算相邻实体单元之间的应力差值，通过应力差值判断实体单元的应力处于受拉状态还是受压状态，然后修正扩散系数D。；
[0009]5)建立三维非均匀单元的元胞自动机模型，计算得到进化系数；
[0010]6)通过进化系数计算得到时间差值内每个实体单元各面的氯离子浓度；
[0011]7)重复以上步骤迭代计算得到各个时期的氯离子浓度，完成对氯离子在混凝土中扩散与渗透作用的模拟。
[0012]所述的步骤2)中将实体单元根据暴露在空气中表面数目进行分类具体采用以下方式分类:没有表面暴露在空气中的归位第一类，有一个面暴露在空气中的归位第二类，有两个面暴露在空气中归位第三类，有三个面暴露在空气中归位第四类。
[0013]所述的步骤3)中修正该实体单元进行元胞自动机计算时的环境氯离子浓度具体采用以下方式:当实体单元的表面暴露在内部空腔中时，设定环境氯离子浓度为O ;当实体单元的表面暴露在空气中时，设定环境氯离子浓度为现场实测的环境氯离子浓度。
[0014]所述的步骤4)中扩散系数D。采用以下公式修正:
[0015]D0=A1BaWA2BJD
[0016]式中，D。为氯离子的扩散系数；D为不考虑应力的扩散系数；B。为应力比值，ApA2为第一、第二经验系数。
[0017]所述的步骤5)中进化系数采用以下公式计算得到:
[0018]
【权利要求】
1.一种混凝土氯离子扩散与渗透作用的模拟方法，其特征在于包括如下步骤: 1)将桥梁结构离散化，用实体单元建立桥梁结构的空间模型； 2)进行前处理，将实体单元根据暴露在空气中表面数目进行分类； 3)对暴露在空气中的实体单元进行进一步的处理，判断识别该实体单元是否位于箱梁结构的内部空腔壁上，然后修正该实体单元进行元胞自动机计算时的环境氯离子浓度； 4)计算相邻实体单元的氯离子浓度差，再根据桥梁结构设计规范对结构施加荷载计算相邻实体单元之间的应力差值，通过应力差值判断实体单元的应力处于受拉状态还是受压状态，然后修正扩散系数D。； 5)建立三维非均匀单元的元胞自动机模型，计算得到进化系数； 6)通过进化系数计算得到时间差值内每个实体单元各面的氯离子浓度； 7)重复以上步骤迭代计算得到各个时期的氯离子浓度，完成对氯离子在混凝土中扩散与渗透作用的模拟。
2.根据权利要求1所述的一种混凝土氯离子扩散与渗透作用的模拟方法，其特征在于:所述的步骤2)中将实体单元根据暴露在空气中表面数目进行分类具体采用以下方式分类:没有表面暴露在空气中的归位第一类，有一个面暴露在空气中的归位第二类，有两个面暴露在空气中归位第三类，有三个面暴露在空气中归位第四类。
3.根据权利要求1所述的一种混凝土氯离子扩散与渗透作用的模拟方法，其特征在于:所述的步骤3)中修正该实体单元进行元胞自动机计算时的环境氯离子浓度具体采用以下方式:当实体单元的表面暴露在内部空腔中时，设定环境氯离子浓度为O ;当实体单元的表面暴露在空气中时，设定环境氯离子浓度为现场实测的环境氯离子浓度。
4.根据权利要求1所述的一种混凝土氯离子扩散与渗透作用的模拟方法，其特征在于:所述的步骤4)中扩散系数D。采用以下公式修正:
D0= A1B 0 2+A2B 0 +D 式中，D。为氯离子的扩散系数；D为不考虑应力的扩散系数；B。为应力比值，A1^ A2为第一、第二经验系数。
5.根据权利要求1所述的一种混凝土氯离子扩散与渗透作用的模拟方法，其特征在于:所述的步骤5)中进化系数采用以下公式计算得到:
其中，Φ?*与关注单元相邻的单元进化系数，CK为关注单元的进化系数，At为元胞自动机迭代计算的时间差值；κ为溶液的渗透系数；D。为氯离子的扩散系数;Ai为实体单元各个面的面积；
为相邻单元之间的压力梯度；
为相邻单元之间的浓度梯度；P s为含有氯离子溶液的质量密度；g为重力加速度，i表示实体单元各面的序数，Δ Si为关注单元及其相邻单元之间的间距，AVtl为关注单元的体积。
6.根据权利要求1所述的一种混凝土氯离子扩散与渗透作用的模拟方法，其特征在于:所述的步骤6)中时间差值内每个实体单元各面的氯离子浓度通过以下公式计算得到:
式中，c (X，y, z, t+ Δ t)为实体单元在t+ Δ t时刻氯离子的浓度；c (x, y, z, t)为实体单元在t时刻氯离子的浓度；c(x+AXi, y+Ayj, ζ+Δ Zi, t)为与该实体单元相邻的实体单元在t时刻氯离子的浓度；Axi为该实体单元与相邻实体单元在X方向上的间距，Ayi为该实体单元与相邻实体单元在I方向上的间距，Azi为该实体单元与相邻实体单元在z方向上的间距。
【文档编号】G01N17/00GK104181091SQ201410407275
【公开日】2014年12月3日 申请日期:2014年8月18日 优先权日:2014年8月18日
【发明者】项贻强, 吴强强 申请人:浙江大学