多孔连续介质渗透压力对变形应力影响的数值模拟方法探讨-水利学报.PDFVIP

水 利 学 报 2017年6月 SHUILI XUEBAO 第48卷 第6期 文章编号：0559-9350（2017）06-0640-11 多孔连续介质渗透压力对变形应力影响的数值模拟方法探讨 张国新 （流域水循环模拟与调控国家重点实验室 中国水利水电科学研究院，北京 100038） 摘要：渗流场对应力和变形影响计算中渗流荷载普遍采用渗透体积力+浮托力方式，这是一种简化近似算法，即 忽略了渗透压力作用下固体介质自身的变形。这种近似对于土、砂等细颗粒离散介质的误差是可以忽略的，但对 于混凝土、裂隙岩体等多孔连续介质或等效孔隙介质则因忽略了孔隙压力对固体介质变形的影响，夸大了浮托力 的作用，导致变形和有效应力失真。本文从孔隙水压作用下的应力应变关系出发，将孔隙压力的作用按初应变 （或初应力）推导有限元渗透荷载公式（孔隙压力初应力/应变法）。用算例比较了不同方法对结果的影响，结果表 明，渗透体积力+浮托力的算法，有效应力和变形误差过大，有些情况下会得到错误的变形结果。而孔隙压力初 应力/应变法，可以正确地得到渗流场中渗透压力作用下的应力和变形。 关键词：多孔连续介质；渗流场；孔隙压力；初应变；有效应力 中图分类号：O357.3 文献标识码：A doi ：10.13243/ki.slxb 1 研究背景 我国近几年相继建成了小湾、锦屏和溪洛渡等多座300m级高拱坝，现在已完成了初期蓄水，进 入正式运行阶段，三峡大坝也已建成蓄水运行近10年。在初次蓄水运行过程中，取得了大量观测资 料。笔者在对这些坝的资料进行分析时发现特高坝蓄水后水文地质条件变化对大坝及岸坡变形影响 显著，但用传统的模拟方式，不管是采用面力，还是渗流体积力对渗流场的影响进行模拟，得到的 变形都与观测结果有较大出入。 关于混凝土重力坝的渗透压力计算，张光斗先生早在1956年发表了著名的论文《混凝土坝的渗透 ［1］ ［2］ 压力》 ，提出了坝基扬压力和坝内应力计算方法，奠定了重力坝扬压力计算的基础。朱伯芳院士 从孔隙水压力影响下的应力-应变关系入手，推导了渗透水对非均质重力坝应力状态的影响。早期的 研究主要集中在渗透力对重力坝坝基场压力的分布形式方面。采用数值方法计算渗流场及渗流场作 ［3］ 用下建筑物的应力和变形，国内最早见于张有天的工作 ，文中首次给出了采用渗透体积力计算渗透 荷载的计算公式：X i = ∂P ∂x i 。公式可分为两部分，即与水力梯度成正比的渗透力和竖向浮托力， f = -γ∂H/∂x f = -γ∂H/∂y f = -γ∂H/∂z + γ H γ 即 ， ， ，其中 为总势能水头， 为水容重。 x y z ［4］ 潘家铮院士 在1984年详细论述了混凝土坝有限元分析中水荷载考虑方式，提出了分三阶段数 值模拟思路，即：在蓄水期按表面荷载计算；在地基达到稳定渗流期，先计算地基渗流场，按渗透 压力计算，基岩容重取浮容重；坝体达到稳定渗流期，水压仍按面力计算，浸润线以下混凝土按浮 容重计算，坝体水平正应力中减去渗透压力强度，即为有效应力。这套计算思路中渗流荷载的计算 ［3］ 实际上与张有天中的方法相同 ，即渗透体积力加浮托力。 收稿日期：2016-10-17；网络出版日期：2017-06-19 网络出版地址：http：///kcms/detail/11.1882.TV1547.007.html 基金项目：国家重点研发计划项目（2016YFB0201000）；国家自然科学基金项目51439005） 作者简介：张国新（1960-），男，山东临朐人，博士，教授级高级工