水工建筑物课件-重力坝渗流分析、温度应力和温度控制.pptVIP

认识：混凝土和岩体都是透水性材料。 计算：按弹性理论分析坝体应力时，应将渗流压力按渗流体积力计算。 假定：均匀渗流介质，并遵循达西定律。真实情况岩石属于裂隙介质，渗流情况比较复杂。 方法：当边界条件比较复杂时，用解析法很难求解，但采用有限元法则可得到较为满意的结果。 变量：在弹性结构中的场变量是位移或应力，在稳定渗流场中的场变量则为势头等标量函数。 一、稳定渗流场的计算 符合达西定律的平面问题稳定渗流基本方程为: 式中： H——势头函数； x、y——坐标系； kx、ky——沿x、y轴方向的渗流系数。 (一)渗流基本方程和定解条件 已知什么？求什么？怎么求？ 第二类边界C（已知流量边界条件） 第一类边界B（势头值已知） 由稳定渗流基本方程和此二边界条件，可解出势头函数H。 根据变分原理，可得 满足此函数的H便是所求的渗流场在定解条件下的解。 势头函数还应满足一定的边界条件，如图所示: 式中 ：n——边界外法线方向。 单元内势头函数H的导数： 求单元泛函的微分，经演算得： 离散化为m个单元的渗流场，其泛函J为各个单元泛函Je之和： 以三角形单元为例，假定单元内势头函数H为线性函数： 式中的三个常数a1、a2、a3可用三个结点i、j、m的坐标及其相应的势头值Hi、Hj、Hm表示，于是势头函数： (二)用有限元法求解势头函数H 对所有单元的泛函求得微分后迭加，并使其等于0（求极小值），就得到泛函对结点势头的微分方程组： 写成矩阵形式，得稳定渗流的有限元法计算公式为： 式中：［K］—总渗流率矩阵; ｛H｝—结点势头的列阵； ｛F｝—已知常数项列阵，由已知边界条件得出. 求解总体方程组式，即可得出各结点的势头值。下图为由计算得出的渗流场等势线示意图。 渗流在稳定渗流场各点产生的渗流压力，就象结构物的自重、地震惯性力一样，也是一种体积力。当已求出渗流场各点的势头值后，在H中减去该点的纵坐标y再乘以水的容重γ，即可求得各点的渗流压力p: 对三角形单元，设三个结点上的渗流压力分别为pi、pj和pm，则单元荷载如图所示。渗流压力的合力分量为： 将此合力分量三等分后，移置到三个结点上,即可得到作用于每个结点上的渗流体积力为： 和 。 从图中可以看出，两种不同考虑得出的应力差别很大，假定地基不透水时，在坝踵处给出了较大的拉应力。由此可以认为：只有当水库骤然蓄水至设计水位或地基为不透水岩体时，方可以按库底水压力计算，在一般情况下应按渗流体积力计算。 库底在水荷载作用下的坝底面应力σy的分布情况 假定地基不透水、沿库底施加水库水压力时的应力 认为地基透水，在地基渗流体积力作用下的应力。 第一阶段（OP段）为上升期 第二阶段（PF段）为降温期 第三阶段为稳定期 海瓦西坝（坝高93.7m,底宽68m）实测温度过程线 一、坝体温度的变化过程 实体重力坝从Tmax降到Td所需的时间，视散热条件，自几个月至几年不等。而温度回降值△T＝Tmax一Td则是重力坝温度控制的一个重要指标。 第一阶段（OP段）为上升期 混凝土入仓时的浇筑温度为Tj，入仓后，由于水泥水化热作用，在混凝土凝固和硬化过程中产生水化热，使混凝土温度从Tr上升至最高值Tmax。主要温升发生在3~7天内，温差Tr= Tmax-Tj称为水化热温升，其值一般为15~25℃，最高可达36℃。 第二阶段（PF段）为降温期 温度达到Tmax后，逐渐散热冷却，温度下降。靠近坝体表面由于水化热散逸较快，该部位的温度较早地达到稳定温度，如图中A2的PF段。坝体内部由于水化热不易散发，温度下降缓慢，需要较长时间才能到达稳定温度，如图中的A1所示。 第三阶段为稳定期 靠近坝体表面的温度随外界温度变化而波动，坝体内部一般接近坝址的年平均气温Td。 二、施工期的温度应力 （一）地基约束引起的应力 在降温过程中，如不受基岩约束，浇筑块自由收缩，如右图中a’b’和c’d’所示，实际上受到基岩约束后的变形为a’b和c’d 总温降为T＝Tj十Tr－Td 受地基约束混凝土浇筑块由温度下降产生的变形及应力 （二）内外温差引起的应力 产生原因： 混凝土块在浇筑初期，由于表面温度降低，在块体内外形成的温差将使外部收缩受拉而内部受压，如右图所示，当拉应力超过材料的抗拉强度时，将产生表面裂缝。 温度场计算： 施工期坝体温度取决于混凝土入仓时的浇筑温度Tj和水化热温升，由于水化热温升和边界温度都在随时间而变化，因而求解施工期的坝体温度场是很复杂的，一般多采用简化的计算方法，参见《混凝土重力坝设计规范》