COMSOL模拟混凝土碳化过程与工程应用

1. 混凝土碳化现象与COMSOL建模价值

混凝土碳化是钢筋混凝土结构耐久性研究中的核心问题之一。当环境中的二氧化碳渗透进混凝土孔隙，与水泥水化产物氢氧化钙发生化学反应，导致混凝土碱度降低。这个过程可以用以下化学反应式表示：

Ca(OH)₂ + CO₂ → CaCO₃ + H₂O

在COMSOL中复现这一过程的价值在于：

• 可视化碳化前沿的动态发展过程
• 量化不同环境条件下的碳化深度
• 预测结构服役寿命
• 优化混凝土配合比设计

我去年参与的一个地下管廊项目就通过这种建模方法，成功预测了不同保护层厚度下的钢筋初锈时间，为设计院提供了关键参考数据。

2. 模型理论基础与PDE构建

2.1 碳化过程的质量传递方程

混凝土碳化本质是气体在多孔介质中的扩散-反应耦合过程。控制方程包含两个关键部分：

1. 二氧化碳扩散方程： ∂(φc)/∂t = ∇·(D_eff∇c) - R

   其中：

   • φ：混凝土孔隙率（实测值通常在0.1-0.3之间）
   • c：CO₂浓度（mol/m³）
   • D_eff：有效扩散系数（与湿度、温度相关）
   • R：化学反应消耗速率

2. 反应动力学方程： R = k·c·(1-α)

   k = A·exp(-Ea/RT) （Arrhenius方程）

   经验参数：

   • A：指前因子 ≈ 1.2×10⁻⁴ s⁻¹
   • Ea：活化能 ≈ 45 kJ/mol
   • α：碳化程度（0-1）

注意：实际建模时需要根据混凝土类型调整这些参数。普通硅酸盐混凝土与掺粉煤灰混凝土的反应动力学参数可能相差一个数量级。

2.2 边界条件设置技巧

1. 暴露表面：

   • 第一类边界条件：c = c₀（环境CO₂浓度）
   • 典型取值：城市环境400-800ppm，隧道内可达1500ppm

2. 对称边界：

   • 零通量条件：n·(-D∇c) = 0

3. 钢筋界面：

   • 可设置为零通量或添加金属钝化膜模型

我在处理一个海工结构项目时发现，潮汐区的干湿交替边界需要特别处理。这时可以采用时间相关的边界条件：

c(t) = c₀·(1 - exp(-t/τ)) （τ为润湿周期）

3. COMSOL操作全流程解析

3.1 模型建立步骤

1. 选择物理场：

   • 新建模型 → 选择"数学 → PDE接口 → 系数形式PDE"
   • 添加"变形几何"接口用于碳化收缩模拟（可选）

2. 参数定义：

   # 典型参数示例 phi = 0.15 # 孔隙率 D0 = 1e-7 # 参考扩散系数(m²/s) RH = 0.75 # 相对湿度 T = 298 # 温度(K) c0 = 0.04 # 环境CO₂浓度(mol/m³)

3. PDE系数设置：

   • 扩散项：D_eff = D0*(1-RH)^2.5*(φ/0.18)^1.8
   • 反应项：R = kc(1-α)
   • 存储项：φ

3.2 材料属性定义

建议使用COMSOL的"材料库"功能创建自定义混凝土材料：

实操技巧：可以通过"全局定义 → 函数"创建温湿度相关的材料参数表达式，实现环境条件动态变化。

3.3 网格划分策略

1. 边界层网格：

   • 暴露表面添加3-5层边界层
   • 生长率建议1.2-1.5
   • 首层厚度≈预计碳化深度的1/20

2. 整体网格：

   • 自由四面体网格（3D）
   • 最大单元尺寸≤最小特征尺寸的1/5
   • 钢筋周边需要局部加密

一个常见错误是过度细化网格导致计算量激增。我的经验法则是：先进行粗网格试算，再根据浓度梯度变化率调整加密区域。

4. 仿真结果分析与验证

4.1 典型输出结果处理

1. 碳化深度计算：

   • 定义α=0.5的等值面为碳化前沿
   • 使用"派生值 → 截面"功能提取深度-时间曲线

2. 数据导出格式：

   % 导出时间-深度数据示例 time = [0,30,90,180,365]; % 天数 depth = [0,2.1,4.3,6.8,9.5]; % mm save('carbonation.mat','time','depth');

4.2 实验验证方法

实验室加速碳化试验（GB/T 50082-2009）数据可用于模型验证：

验证时需要注意：

• 实验室CO₂浓度(3%)远高于自然环境
• 需进行浓度等效换算：t_field = t_lab*(c_lab/c_field)ⁿ （n≈0.5-0.7，取决于混凝土类型）

5. 工程应用案例与常见问题

5.1 实际工程参数设置

某跨海大桥桥墩碳化预测模型参数：

参数项 取值 来源 ------------------------------------------- 环境CO₂ 0.06% 现场监测 年平均温度 18℃ 气象数据 保护层厚度 50mm 设计图纸 水胶比 0.38 配合比报告 粉煤灰掺量 25% 施工记录

5.2 常见报错与解决方案

1. 发散问题：

   • 现象：求解器报"达到最大牛顿迭代次数"
   • 解决方法：
     1. 减小时间步长
     2. 调整非线性求解器阻尼因子(0.7-1.0)
     3. 检查反应速率单位是否一致

2. 内存不足：

   • 现象：计算中途崩溃
   • 优化策略：
     1. 使用对称模型（1/2或1/4模型）
     2. 改用直接求解器(MUMPS)
     3. 增加虚拟内存页面文件

3. 结果异常：

   • 现象：碳化深度随时间减小
   • 检查要点：
     1. 扩散系数符号是否正确
     2. 边界条件类型是否设置错误
     3. 反应项表达式是否遗漏浓度项

6. 模型扩展与进阶技巧

6.1 多物理场耦合实现

1. 热-湿-碳化耦合：

   • 添加"热传递"和"湿气传递"接口
   • 耦合关系：
     • 温度影响反应速率：k~exp(-1/T)
     • 湿度影响扩散系数：D~f(RH)
     • 碳化放热：Q=ΔH·R

2. 力学性能退化：

   • 使用固体力学接口
   • 定义弹性模量随α变化： E(α) = E0·(1 - 0.2α)

6.2 参数反演方法

通过实测数据反演关键参数：

1. 创建"优化"研究
2. 定义目标函数： min Σ(d_exp - d_sim)²
3. 设置优化变量：
   • D_eff
   • k
   • φ

使用Levenberg-Marquardt算法时，建议先对参数进行归一化处理，避免因量级差异导致收敛困难。