重力坝分析的改进有限元等效应力法 Study and application of the improved finite element equivalent stress method to analyze the strength of gravity dams.pdfVIP

第11卷第2期 中国水利水电科学研究院学报 V01．11No．2 2013年6月 June,2013 ofChinaInstituteofWaterResourcesand Research Journal Hydropower 文章编号：1672—3031(2013)02—0112—05 重力坝分析的改进有限元等效应力法 杨会臣，贾金生，郑璀莹 (中国水利水电科学研究院流域水循环模拟与调控国家重点实验室，北京100038) 摘要：本文提出了重力坝分析的改进有限元等效应力法，分析了改进前后应力计算结果的变化，研究表明，当网 格尺寸较大时，改进的等效应力方法具有更高的精度。对比研究了采用有限元等效应力方法、材料力学方法计算 应力在结果上的相同点和不同点。分析表明，有限元等效应力法在解决应力集中的同时，还能反映地基刚度对建 基面竖向正应力分布的影响，从而解决了材料力学的一些限制。结合我国第一座胶结颗粒料永久工程一守口堡工 程的设计，详细分析了地基刚度对坝基面应力分布的影响，取得了很好效果。 关键词：重力坝；有限元等效应力方法；材料力学方法；守口堡工程 中图分类号：TV641．3 文献标识码：A l 改进的有限元等效应力法 有限元等效应力法由我国学者首先提出¨“，其 基本思想是将某一截面上有限元计算的应力结果等效 ， 为作用在该截面上的一组内力，然后利用材料力学的 方法计算截面上的应力分布，以避免角缘处有限元计 J Xi+L 算结果的应力集中。有限元等效应力法分析拱坝的文 一 一●1 献¨’7’很多，但用于研究重力坝应力分布的文章还比 图1有限元应力计算网格 较少m…。 对于图1所示的有限元计算网格，t处的竖向正应力为盯i7，Xi+I处的竖向正应力为盯i：，，则该截面 上的等效内力为¨1： Ⅳ：艺f之尘b (1) 一 ”1L ／ 肘=薯[华h孚飞] ㈤ 式中：Ⅳ为截面轴力；M为截面弯矩；X0为截面转动中心；liX川一菇i；B为截面宽度，取1。 根据材料力学偏心受压构件正应力计算公式，结点彼E的有限元等效应力即为： 盯i’：等+竽掣 ㈥ 其中：三为截面长度；，为截面绕‰的转动惯量。 在利用式(2)进行取矩计算时，将合力的作用点作用在两结点的中点，其默认了一个假定，即两结点 之间应力是相等的，即矩形分布，而实际应力分布是一个梯形，根据实际的应力分布，修改式(2)如下： 收稿日期：2012—11一ll 作者简介：杨会臣(1984一)，男，山东人，博士生，主要从事水工结构研究。E-mail：yanghc@iwhr．tom 一112一 万方数据 重力坝分析的改进有限元等效应力法杨会臣贾金生郑璀莹 肘=氧华H”鞴“+--剞飞] (4) 为了研究改进前后等效应力计算结果的差异，选取