坝基弹性模量对广源拱坝静力特性影响分析

　　[摘要]广源水电站拱坝坝基岩体风化、卸荷程度较强，弱风化层厚度较大，岩体的变形模量差异大。本文采用有限元法分析坝基弹性模量对拱坝静力特性的影响，为广源水电站拱坝的布置提供依据。分析结果表明：坝顶以上区域的岩体只对拱坝上部坝肩的位移和应力有一定影响;两岸坝肩、河床及以下区域岩体对拱坝各部位的位移和应力均有影响，随着这两个区域岩体弹性模量增加，坝体三个方向的位移都会减小;当两岸坝肩岩体弹性模量增加时，上部、中部坝肩的拉应力增加，下部坝肩和拱冠梁底部的拉、压应力在数值上均减小;当河床及以下区域岩体弹性模量增加时，中部、下部坝肩以及拱冠梁底部的拉应力增加;广源水电站拱坝坝基弹性模量处于影响敏感区。

　　[关键词]拱坝;弹性模量;静力特性;有限元法;广源水电站

　　1工程概况

　　广源水电站位于广西东北部资源县河口乡境内的五排河上，是一座以发电为主，兼有向漓江补水、防洪、旅游等综合效益的中型水利水电工程。大坝为混凝土双曲拱坝，建基面最低高程为386m，正常蓄水位为478m，对应下游水位为400.35m，高坝方案最大坝高为95.00m，坝顶高程为481.00m。选定的下坝址两岸地形平顺，无大的冲沟发育且两岸基本对称，岸坡自然坡角约45°，宽高比2.35～2.38。

　　电力论文投稿期刊：《甘肃水利水电技术》(月刊)创刊于1956年，是经国家科委批准，出版管理部门登记注册，甘肃省水利厅主办，公开出版发行的水利水电综合性科技期刊。

　　坝基岩体为硅化泥岩，属中硬岩，河床坝段以厚层状为主，两坝肩及下游抗力体则以中厚层状为主，地层倾向右岸偏上游，局部有泥化夹层。根据坝址区10个平硐揭露，强风化水平厚度0～32m，弱风化水平厚度24.80～36.0m，强卸荷水平厚度32～42m，风化、卸荷总体程度较强，相比较而言，两岸比河床强烈。根据坝基岩体条件，设计初步建议河床坝段采用微风化～新鲜完整的BⅡ类岩体作为坝基;两岸中、高拱圈采用微风化～弱风化下部的BⅡ类或BⅢ类岩体作为坝基。由于BⅡ类和BⅢ类岩体的变形模量差异较大，采用不同的开挖标准对拱坝体型布置和工程投资有较大的影响。

　　目前，国内外学者就坝基弹性模量对拱坝的影响已做了诸多有益的研究，付锦华[1]以西北口水库为背景开展了弹性模量的选择及其对拱坝应力分析的影响研究，初步揭示了基础和坝体弹性模量比值对坝体应力变化规律的影响;王亮清等[2]以鱼简河水电站双曲拱坝为背景，研究坝基不同变形模量对拱坝应力分布的影响，研究表明中间软弱夹层的变形模量对拱端和拱冠应力影响较大;赵珍等[3]针对某碾压混凝土拱坝首次引入了“有限元等效应力”，研究不同坝基弹性模量对拱坝应力应变影响及变化趋势规律。这些研究大多是针对具体工程，广源水电站拱坝有自身的特点，因此，有必要研究广源水电站坝基岩体变形模量对坝体变形与应力的影响，为拱坝的布置提供依据。

　　2方法与模型建立

　　本文采用有限元法分析广源水电站坝基弹性模量对拱坝静力特性的影响，具体分析采用ANSYS软件实现。广源拱坝基岩的特点是岩性单一，风化深度是河床较浅、两岸较深，且两岸风化深度随高程增加而增加。为了反映基岩的特点，进行有限元法建模时，将坝基在高度方向按坝顶以上(Ⅰ区)、两岸坝肩(Ⅱ区)、河床坝基(Ⅲ区)水平划分为上、中、下3个区域，建立拱坝—地基三维有限元模型，分别改变坝基与坝体混凝土弹性模量比值，研究坝基弹性模量对拱坝静力特性的影响。根据拱坝坝高和坝址的地形地质条件，确定有限元计算的模拟范围为：基岩底部高程取286m，从拱坝左右坝肩分别向左右两岸各取125m，顺河流向从坝轴线向上游取170m，从坝轴线向下游取260m，坝基有限元模型按照实际地形建立。

　　整体坐标系为X向右岸为正，Y向下游为正，Z向上为正的右手坐标系。应力、位移符号与弹性力学规定一致，即拉为正，压为负。模型底部坝基为固端约束，其他侧面为单向链杆约束。有限元模型采用SOLID185高阶3D8节点实体单元，坝体单元数为8091个，节点数为8349个，坝基单元数为135544个，节点数为78641个。

　　3坝基弹性模量对拱坝静力特性的影响

　　为分析坝基不同部位的弹性模量变化对拱坝静力特性影响，得出不局限于广源水电站的一般规律，分析中坝体混凝土弹性模量Ec取为定值，Ec=25.5GPa(相应于C20混凝土)，采用不同的坝基弹性模量Er与坝体混凝土弹性模量Ec的比值来反映坝基弹性模量的变化，计算考虑了4种方案。

　　①方案一：保持Ⅱ和Ⅲ区坝基弹性模量不变，且Er/Ec=1，Ⅰ区坝基弹性模量按Er/Ec=0.1、0.2…1、1.2、1.5、2、2.5、3、3.5、4改变;②方案二：保持Ⅰ和Ⅲ区坝基弹性模量不变，且Er/Ec=1，Ⅱ区坝基弹性模量按Er/Ec=0.1、0.2…1、1.2、1.5、2、2.5、3、3.5、4改变;③方案三：保持Ⅰ和Ⅱ区坝基弹性模量不变，且Er/Ec=1，Ⅲ区坝基弹性模量按Er/Ec=0.1、0.2…1、1.2、1.5、2、2.5、3、3.5、4改变;④方案四：Ⅰ、Ⅱ、Ⅲ三块区域作为一个整体，同时按Er/Ec=0.1、0.2…1、1.2、1.5、2、2.5、3、3.5、4改变。

　　计算时坝体混凝土的重度和泊松比分别取为24kN/m3和0.2，坝基岩体的泊松比均取为0.32。拱坝的荷载组合考虑以下主要荷载：自重、静水压力、扬压力和温度荷载，选择对拱坝应力较为不利的正常蓄水位+温降工况组合作为代表，各种荷载按《混凝土拱坝设计规范》(SL282-2018)进行计算。

　　3.1位移影响分析

　　有限元法计算结果表明，4种方案拱坝的位移分布规律一致，在荷载作用下，拱坝坝体产生的位移主要是顺河向位移，横河向和竖向的位移均较小。最大顺河向位移均发生在拱冠梁顶部，横河向位移基本上是以拱冠梁为对称轴，左右半拱分别向左右岸变位，竖向向下位移很小，最大竖向位移均发生在拱梁顶部下游侧。

　　3.2应力影响分析为了方便对坝体应力进行对比分析，在坝体上选取代表截面，通过代表截面上、下游边缘点(特征点)的应力来分析坝基弹性模量对坝体应力的影响。根据广源水电站拱坝特点，拱坝体形基本上以拱冠梁为轴左、右对称，因此，在左岸坝肩上、中、下取3个截面(编号：1、2、3)作为代表截面来分析左、右岸坝肩的应力，由于顶部布置有溢流表孔，在拱冠梁上只取中、下2个截面(编号：4、5)作为代表截面来分析拱冠梁的应力。

　　4结论

　　(1)坝顶以上区域岩体弹性模量的变化，除对拱坝上部坝肩的位移和应力有一定影响外，对拱坝中、下部坝肩及拱冠梁的位移和应力没有影响;坝顶以上区域岩体的弹性模量变化对拱坝上部坝肩的位移影响较小，而对拉应力的影响较大，上部坝肩的拉应力随该区域岩体弹性模量的增加而增加。

　　(2)两岸坝肩岩体弹性模量的变化对拱坝各部位的位移和应力均有影响，随两岸坝肩岩体弹性模量增加，顺河向(径向)、横河向(切向)和竖向的位移减小，而上部和中部坝肩的拉应力均增加，下部坝肩和拱冠梁底部的拉、压应力在数值上均减小。

　　(3)河床及以下区域岩体弹性模量的变化对拱坝各部位的位移和应力均有影响，随河床及以下区域岩体的弹性模量增加，顺河向(径向)、横河向(切向)和竖向的位移减小，而中部、下部坝肩和拱冠梁底部的拉应力均增加，上部坝肩的拉应力会减小。

　　(4)当坝基岩体弹性模量小于0.5倍坝体弹性模量时，坝体的位移和应力对坝基岩体弹性模量变化敏感，当坝基岩体弹性模量为0.5~1.0倍坝体弹性模量时敏感度降低，当坝基岩体弹性模量大于坝体弹性模量时，对坝基岩体弹性模量变化不敏感。(5)广源水电站坝基上部弱风化岩体的弹性模量是2GPa，下部弱风化岩体的弹性模量是3GPa，微风化岩体的弹性模量是5GPa，Er/Ec=0.078~0.2，处于敏感区。

　　参考文献

　　[1]付锦华.弹性模量的选择及其对拱坝应力分析的影响[J].水利水电技术，1983(7)：25-30.

　　[2]王亮清，唐辉明，胡新丽，等.地基变形模量对碾压混凝土拱坝应力分析的影响[J].地球科学(中国地质大学学报)，2007(3)：409-414.

　　[3]赵珍，李守义，刘鹏，等.地基弹性模量对拱坝坝体应力的影响研究[J].水利水电技术，2017(12):56-62.

　　作者：何飞龙，覃振威，张小飞，田羽，黄佳敏