主墩承台水化热分析计算书.docxVIP

主墩承台水化热分析计算书

一、工程概况

本工程主墩承台尺寸为长×宽×高=20m×15m×4m，混凝土设计强度等级为C40。承台处于正常环境条件下，主要承受上部结构传来的竖向荷载和水平荷载。采用商品混凝土进行浇筑，混凝土原材料包括水泥、砂、石子、外加剂和水。其中，水泥选用P.O42.5普通硅酸盐水泥，砂为中砂，石子粒径为525mm连续级配碎石。

二、水化热分析基本理论

水泥水化是一个放热反应过程，在混凝土内部产生热量积聚，由于混凝土导热性能较差，内部温度升高，而表面散热较快，从而形成温度梯度，产生温度应力。当温度应力超过混凝土的抗拉强度时，就会导致混凝土出现裂缝。

根据热传导理论，混凝土内部的温度场满足三维热传导方程：

$\frac{\partialT}{\partialt}=\alpha(\frac{\partial^{2}T}{\partialx^{2}}+\frac{\partial^{2}T}{\partialy^{2}}+\frac{\partial^{2}T}{\partialz^{2}})+\frac{Q(t)}{\rhoc}$

式中：

$T$为混凝土内部某点的温度（℃）；

$t$为时间（h）；

$\alpha=\frac{\lambda}{\rhoc}$为混凝土的导温系数（$m^{2}/h$），其中$\lambda$为混凝土的导热系数（$W/(m\cdotK)$），$\rho$为混凝土的密度（$kg/m^{3}$），$c$为混凝土的比热容（$J/(kg\cdotK)$）；

$Q(t)$为水泥水化热在单位时间内的发热量（$J/(m^{3}\cdoth)$）。

水泥水化热$Q(t)$通常采用指数函数形式表示：

$Q(t)=Q_{0}(1e^{mt})$

式中：

$Q_{0}$为水泥最终水化热（$J/kg$）；

$m$为与水泥品种、浇筑温度等有关的系数；

$t$为时间（h）。

三、计算参数确定

1.混凝土物理性能参数

混凝土密度$\rho=2400kg/m^{3}$

混凝土比热容$c=0.96\times10^{3}J/(kg\cdotK)$

混凝土导热系数$\lambda=2.3W/(m\cdotK)$

导温系数$\alpha=\frac{\lambda}{\rhoc}=\frac{2.3}{2400\times0.96\times10^{3}}\approx1.0\times10^{6}m^{2}/s=0.0036m^{2}/h$

2.水泥水化热参数

水泥最终水化热$Q_{0}=375J/kg$

根据经验，系数$m=0.3$

3.边界条件

混凝土浇筑时的初始温度$T_{0}=20℃$

环境温度$T_{a}=20℃$，假设环境温度在计算过程中保持不变。

混凝土表面与空气的对流换热系数$\beta=10W/(m^{2}\cdotK)$

四、有限元模型建立

采用有限元软件MidasCivil进行水化热分析。将承台简化为三维实体模型，采用六面体单元进行网格划分。考虑到承台的对称性，取1/4模型进行分析，以减少计算量。

1.几何模型

建立1/4承台几何模型，尺寸为长10m、宽7.5m、高4m。

2.网格划分

采用六面体单元对模型进行网格划分，单元尺寸为0.5m×0.5m×0.5m，共划分单元1200个，节点1681个。

3.材料属性定义

定义混凝土材料属性，输入上述确定的物理性能参数和水泥水化热参数。

4.边界条件施加

在对称面上施加对称约束。

在混凝土表面施加对流换热边界条件，根据牛顿冷却定律$q=\beta(TT_{a})$，其中$q$为热流密度（$W/m^{2}$），$T$为混凝土表面温度，$T_{a}$为环境温度。

五、水化热分析计算过程

1.计算时间步长确定

为了保证计算结果的准确性，采用较小的时间步长进行计算。初始时间步长取$\Deltat=1h$，随着计算的进行，根据温度变化情况适当调整时间步长。

2.迭代计算

根据热传导方程和边界条件，采用有限元软件进行迭代计算，求解混凝土内部各节点的温度随时间的变化。在每个时间步长内，更新水泥水化热$Q(t)$，并重新计算温度场。

六、计算结果分析

1.温度场分布

不同时刻温度分布云图

通过有限元计算，得到承台在不同时刻的温度分布云图。在浇筑后12h，混凝土内部温度开始升高，由于水泥水化热的作用，承台中心部位温度升高较快，而表面温度受环境影响升高较慢，此时温度梯度较小。在浇筑后36h，承台内部温度达到峰值，中心部位最高温度约为65℃，表面温度约为30℃，温度梯度较大。随着