不连续岩体拱坝系统的动力过程数值模拟.pdfVIP

首届 FEPG 用户年会论文集 不连续岩体拱坝系统的动力过程数值模拟 蔡永恩 罗刚 刘今朝 殷有泉 北京大学地球物理系 北京大学力学与工程系 Qingsong Li 张伯艳 陈后群 梁国平 Department of Geological Sciences 中国水利水电科学研究院 中国科学院数学与 University of Missouri Columbia, MO 65211 USA 系统科学研究院 关键词：拱坝 不连续岩体 数值动力分析 安全系数 岩石基础内部的不连续面（断层、节里）在拱坝稳定性分析中是不容忽视的。本文把 由拱坝和这种岩石基础组成的系统称简为拱坝系统。分析这种系统的常用数值方法目前有有 [1-4] [5] [6] [7] [8] [9] 限元方法 ，离散元方法 ，边界元 ，DDA 和流型方法 以及 LDDA 方法。第 2 到第 4 种方法在数学上都属于罚函数方法，这个方法的优点是不必引进新的求解向量，缺点是约 束条件的满足和接触力受刚度取值的影响。LDDA 方法在数学上属于拉格朗日乘子法，需要 补充位移约束方程才能求解，这种方法求得的位移能严格满足接触条件，所得到的接触力是 精确的，这点对于判断系统的稳定是十分重要的。LDDA 方法就是一种建立在后者基础和区 域分解技术上的方法, 本文是这种方法的继续发展和在拱坝系统中的应用。 拱坝动力分析过程可以概括为：1）模拟初始岩石基础应力场，2 ）求建坝后坝体系统 的位移场和应力场，3 ）通过输入与时间有关的载荷模拟拱坝系统的动力过程，4 ）对拱坝系 统进行稳定性评估。 建坝过程的数值模拟方法 岩石基础不连续面的摩擦滑动依赖于初始应力场。恢复这个应力场是非常困难的。为此 在工程实践中常用自重产生的应力场代替初始应力场，在此基础上计算建坝引起的附加应力 场，同时考虑水载荷，泥沙载荷，渗透压力载荷，温度差载荷对坝和岩体的应力场影响。初 始应力场σ 和建坝引起的应力场∆σ 的有限元平衡方程可以分别写成： 0 T ( 1) Β σ dΩ F +G Λ ∫ 0 0 0 0 0 0 Ω 0 BT ∆σ dΩ F +G Λ (2) ∫ d d d Ω Ω dΩ + 0 建坝后的总应力场为： σ =σ + ∆σ (3) 0 将 ∆σ σ−σ0 代入（2 ）并利用本构关系 σ DBU 和（1）式可得： 7 首届 FEPG 用户年会论文集