土石坝因其取材便捷、施工简易的特点,成为水库建设中最受欢迎的坝型,占中国近十万座大坝的约93%。在土石坝建设中,有效防水是核心挑战之一。复合土工膜,凭借其出色的防水性、强适应形变能力、卓越的化学稳定性和施工便利性,被广泛应用,比如塘房庙土工膜心墙坝和西霞院土工膜斜墙坝等项目展示了其成功应用。这些案例激励更多土石坝工程采纳复合土工膜作为防渗材料。然而,研究指出,土石坝在施工和运行期间会经历显著结构变形,可能引起土工膜拉伸应力集中,要求在设计时深入考虑土工膜的力学变形特性。
本研究聚焦于一座应用复合土工膜防渗的中型水库石渣坝,运用三维非线性有限元分析方法,构建了精确的坝体模型,包括土工膜铺设与水库蓄水等关键环节,以全面评估坝体与土工膜在蓄水前后的应力变形行为,旨在验证该复合土工膜防渗设计的科学性和结构性合理性,进一步推动土石坝工程中复合土工膜应用的技术进步与安全实践。
趾板与防渗墙分别使用C25与C20标号混凝土建造。基岩模型采用线弹性理论,设定强风化与弱风化基岩的弹性模量为0.28和1.00GPa,泊松比均为0.3。鉴于复合土工膜薄,不宜实体单元模拟,转而利用三维平面应力单元进行分析。
设计实践中,确保土工膜承受低应力,其响应保持线弹性范围,故采用线弹性模型描述。通过拉伸测试,得到土工膜横向和纵向弹性模量为9.58MPa和7.64MPa,泊松比0.32和0.35。模型中,土工膜与垫层间设接触面,考虑法向刚度变异性,参数设定为:K1=3000,n'=0.56,R'f=0.86,摩擦角φ=20°,最大法向压缩Vm=0.05m,单元厚e=0.10m,以及Kn=5×10^5无因次模量。
按照SL/T 226-98规范,土工膜的许用拉伸强度Ta计算公式基于多个影响因素,包括施工(FiD通常1.1~2.0)、蠕变(FcR为2.0~3.0)、化学与生物破坏(FcD、FbD分别1.0~1.5和1.0~1.3),以及材料极限抗拉强度Tmax。实践中,常取FiD=1.6、FcR=2.0、FcD=1.3、FbD=1.2,导致安全系数K定义为土工膜极限强度的五分之一,即K≥5时满足安全需求。