基于Drucker-Prager模型的泥化夹层地基混凝土重力坝坝基应力分析.pdfVIP

龙源期刊网

基于Drucker-Prager模型的泥化夹层地基

混凝土重力坝坝基应力分析

作者：周怀瑜

来源：《E动时尚·科学工程技术》2019年第08期

龙源期刊网

龙源期刊网

龙源期刊网

龙源期刊网

摘要：目前国内外重力坝发展迅猛，在设计和施工中常遇到各种地基情况，特别是在应

力集中的坝趾或坝踵处，常规的计算方法难以准确、全面地得到其受力状态。利用通用有限元

软件midasFEA，建立大坝三维模型，对其施加荷载并计算得到应力分布于变形分布状态云图

以及输出特殊点处的应力。计算表明：在坝趾和坝踵处有超标主拉应力。建议对应力集中区域

进行地基加固、提高坝基材料的力学性能和均匀性以及改进坝体结构设计，使其适应斜坡坝段

的侧滑趋势和减小坝趾附加部分的弯曲效应。

关键词：重力坝;泥化夹层地基;三维模型;有限元法

1概况

1.1工程概况

瓦村水电站工程坝址位于郁江上游右江河段田林县境内、驮娘江与西洋江汇合口下游9km

处的瓦村水文站附近峡谷中，距下游百色水利枢纽约105km。是郁江综合利用规划的第2个梯

级，是一座以发电为主，兼顾供水、防洪、航运等综合利用的水电站工程。水库总库容5.36

亿m3，有效库容为2.25亿m3，设计装机约230MW，属Ⅱ等工程，主要建筑物级别为2级，

拦河坝的设计洪水标准为500年一遇，校核洪水标准为2000年一遇。

1.2工程特点

12#坝块坝顶高程为311米，最大坝高91米，是瓦村水电站挡水坝块之一。坝基地质条件

为：坝基上游为T2b14地层（约占2%），下游为T2b13地层（约占33%）和T2b12地层（约

占65%），T2b14、T2b13为弱风化中部岩体，T2b12为微风化岩体。12#坝块为最高挡水坝

龙源期刊网

块，坝踵有部分区域处于弱风化层（软弱夹层区域），该区域对整个坝体的稳定应力影响重

大，尤其是施工期。

1.3模型建立

采用midasFEA建立12#坝块三维有限元模型，计算模型整体坐标系采用笛卡尔直角坐标

系，以从上游指向下游的水平向为x轴向，指向下游为正;以坝轴线方向为y轴向，指向左岸

为正;以铅垂向为z轴向，指向上为正。坝体、坝基应力符号以拉应力为正，压应力为负。坝

基的计算宽度和深度均取坝底宽度的3倍，计算厚度取坝体厚度。基岩的上下游、坝轴向两侧

边界加置水平链杆支撑，底部边界固定，坝体部分三向自由。图1-1、1-2为有限元计算模型。

单元网格主要采用六面体，局部为四面体。不考虑坝体中排水廊道和灌浆廊道孔洞对坝体应力

的局部影响。

1.4荷载及材料参数

根据12#坝块的开挖情况，模型不计入地应力，主要荷载考虑自重、水重、静水压力、扬

压力、浪压力及淤沙压力。重力加速度g=9.81m/s2。混凝土的弹性模量E=22000Mpa，泊松比

µ=0.167，容重γ=24KN/m3;坝基各岩层的物理力学参数见表1-1。各工况计算荷载简图1-2。

2计算分析

2.1计算本构模型

整体模型采用基于如下式所示的DrucKer-Prager屈服准则的弹塑性本构模型：

式中：分别为第一应力不变量和应力偏量第二不变量，和与材料的凝聚力c和内摩擦角有

关[1]。

2.2线性有限元计算成果

采用线性有限元法对正常蓄水位、设计洪水位、校核洪水位、完建空库4种工况进行了计

算。

完建空库工况，坝体重心位于上游侧，整个坝体向上游变形，坝踵受压，坝趾受拉;正常

蓄水位、设计洪水位、校核洪水位工况下，坝体向下游变形，坝趾受压，坝踵受拉。坝基面角

点位置简图见图2-1。

附图2-2～附图2-5分别给出12#坝块完建空库、正常蓄水位、设计洪水位、校核洪水位四

个工况在代表断面（y=5.5m断面）处的坝体主拉应力附图。

龙源期刊网

附图2-6～附图2-9分别给出12#坝块完建空库、正常蓄水位、设计洪水位、校核洪水位四

个工况在代表断面（y=5.5m断面）处的坝体主压应力附图。

正常蓄水位、设计洪