拱坝等效应力计算.pdf

第五章 有限元法的拱坝应力控制标准的探讨 41 第五章 有限元法的拱坝应力控制标准的探讨 5.1 概述 目前，有限元法已经广泛应用于拱坝的应力分析，它在计算理论上比多拱梁法先进， 而且适合各种复杂的体形和结构布置形式，可以比较合理的考虑拱坝的整体作用和复杂 地基的影响，同时也可以考虑大孔口、复杂基础、重力支墩、不规则外形等多种因素的 影响，并可以进行仿真计算。但是，在弹性理论的假定下，有限元应力成果往往与所采 用的单元类型和单元网格划分的尺寸有着密切的关系，特别是在建基面附近存在着严重 的应力集中现象，使得有限元应力成果的数值稳定性较差，难以建立相应的应力控制标 准。实际工程中，由于岩体内存在着大小不等的各种裂隙，应力集中现象将有所缓和， 不一定像计算结果那么严重，计算中如何考虑应力集中，是有限元法应用于拱坝的“瓶 颈” [29][30][68] 。我国一些学者提出了有限元等效应力法克服了这一“瓶颈”，新编的《混 [25] 已经正式规定在拱坝的设计中采用有限元等效应 凝土拱坝设计规范》（SL282—2003 ） 力法，为有限元法在拱坝设计中的应用开拓了良好的前景，由于有限元法的强大计算功 能，它将逐步取代多拱梁法，成为拱坝设计的主要方法。 为了解决拱坝的弹性有限元应力成果集中现象，人们采取了许多措施，如规定应力 点的位置、采用厚壳单元、利用有限元计算出坝体应力沿面积积分，得到截面内力，然 后用材料力学方法计算结点应力，即所谓的有限元等效应力法。 傅作新教授[30]提出的有限元等效应力法是基于有限元法的分析结果，将有限元所求 得的应力合成为截面内力，然后求出对应的线性化应力。该方法为有限元分析结果规范 化提出了一种新思路，但是该法截面内力是拟合出来的，不能精确满足内力平衡条件。 此外，该方法要求沿径向网格必须布置3层以上的单元格才能求得结果。 李同春教授[69]提出的基于有限元内力法的等效应力算法和改进的拱坝等效应力分 析方法[70]，该方法假定拱和梁向的正应力与剪应力在拱和梁的截面上的面单元内呈双线 性分布，根据静力平衡原理导出以结点应力为未知量，截面约束内力为右端顶的平衡方 程。在此基础上导出根据截面约束内力求解等效结点应力的公式，然后根据截面约束内 力求解等效结点应力的公式，并在此基础上导出拱和梁向应力为直线分布的上下游面等 效应力的求解方程和主应力的求解公式，但是计算模型复杂。 针对目前拱坝内力计算流行着以二次曲线表示应力分布，即σ a +bx +cx 2 ，朱伯 芳院士在文献[68]指出：从有限元应力计算内力时，以直接采用数值积分方法为宜，利用 三点应力决定二次曲线再求内力的方法是不妥当的，会带来不必要的误差。与此同时， 新修订的水利行业《混凝土拱坝设计规范》（SL282—2003 ）在有限元等效应力法的相关 条款中尽管明确了等效应力的概念和原理，但是还存在以下的问题[71] ：（1）没有明确有 限元计算的模型和内力计算公式，不同的模型和算法可能导致很大的差异；（ 2 ）没有分 42 基于有限单元法的高拱坝体形优化设计 清需要作等效线性化处理的主要应力，笼统地规定由截面内力按材料力学公式计算应 力；（ 3 ）拉应力控制标准的制定原则不够科学，规定的容许拉应力仍无法适合200m 以上 的高拱坝。钱向东教授引入压力容器“分析设计法”中应力分类的概念，进一步阐述有 限元等效应力法的原理和有关计算方法，通过对若干国内已建拱坝的分析，提出基于有 限元等效应力法的拱坝强度设计准则。 5.2 有限元等效应力的计算 5.2.1 截面内力的计算 从力学的观点分析，双曲拱坝是一种具有不规则边界的变厚度壳体，在弹性壳体理 论中，一般取距上下游面等距离的中间曲面（简称中面）为研究对象。用有限元法计算 ( ′ ′ ′) ′ x ,y ,z { } 拱坝，得到整体坐标系 中的应力 σ ，通过拱轴线上任意结点 取局部坐标