关于计算岩土力学发展前景的几点看法.pdfVIP

关于计算岩土力学发展前景的几点看法 沈珠江 （南京水利科学研究院，210029 ；清华大学水利水电工程系，100084 ） 提要 岩土工程建设至今还不能脱离经验，但岩土力学的长足进展已使许多工程的设计逐渐摆脱 完全依赖经验或半经验的状况，数值计算在未来的岩土工程设计中将会占据越来越重要的位置。本文探讨 了计算岩土力学近期可有取得的进展并应加强研究的几个问题。 1、前言 有限元法在我国普及的初期，许多工程师对数值分析能否解决实际问题曾抱着怀疑的态度，但是不 少有识的技术领导还是给予热情的支持。记得上世纪 80 年代初在组织三峡深水围堰第一轮多单位协作分 析计算时，长江水利委员会司兆乐总工曾提出计算分析结果能达到 “精确定性、粗略定量”的目标。到今 天，虽然不能说这一 目标已完全实现，但对相当一部分岩土工程来说，做到这一点已没有困难。当然，岩 土工程的设计和施工在今后相当长的时期内仍需要工程师们的经验，但是，在科学技术飞速发展的今天， 数值分析技术必将越来越成为人们必须依赖的工具。 同样，在岩土力学研究中，计算也已成为和实验一样不可或缺的手段。在某些特殊情况下数值模拟 甚致可以代替实验。而离心模型试验与数值模拟的相互配合，已经成为解决岩土工程问题的一个重大研究 方向。 本文将就 21 世纪岩土工程数值分析的发展前景提出一些看法，基于作者专业知识的局限，重点探讨 土力学问题的数值分析问题。 2 、大型工程的数值模型 大型工程不仅投资大，而且工程安全往往对国计民生有重大影响，数值分析已成为论证其安全性的 一个重要手段。下面举几个例子。 1、三峡深水围堰 三峡深水围堰设计阶段曾进行过两轮由众多单位参加的数值分析。第一轮计算在 1985~1987 期间进行， 当时的主要结论是防渗墙顶部的最大位移可能达 1m 。墙底嵌入基岩处则可能产生较大的 应力而断裂。 第二轮计算在三峡工程开工后于 1993~1995 期间进行，此时的设计方案已作为一些改动，特别是吸取了 第一阶段的研究成果，把刚性混凝土防渗墙改成塑性混凝土防渗墙，同时围堰主体的风化砂的计算参数 也有一定的提高。最后得出墙顶最大位移在 36.7~5 1.0cm 之间，防渗墙底部 应力大大降低，从而基本上 是安全的结论。围堰竣工后已经过1998 年洪水的考验，实际墙顶最大位移为 56.7cm 。 2 、鲁布革心墙堆石坝 鲁布革坝是我国改革开放后第一座按现代技术设计的 100m 级心墙堆石坝。该坝也组织过众多单位参 加的两轮计算。第一轮在设计阶段 1984—1985 期间进行，第二轮则在 1987—1988 期间进行，此时坝体已 大部填筑完成，计算参数根据实际填筑情况重新测定。第一轮计算得出最大沉降在 93~161cm 之间，第二 轮则有所减小，实测值为 60.5 cm 。 3 、上海外环线越江隧道格构式挡墙 外环线沉管隧道两翼格构式挡墙最大高度达 40m 以上，开挖过程中的稳定和变形是令人担心的。为此 进行了三维有限元分析。此格构墙形状十分复杂，计算中布置了22612 个结点。计算得出开挖结束时向江 的最大位移为 15.3cm ，岸上加载且逢低潮位时进一步增大到 2 1.7cm 。由于计算参数是根据经验假定的， 上述计算结果的可信度不会很高，但在说明安全性大体上可以保证的同时，计算结果显示了格构墙整体 呈复杂的扭转变形和混凝土墙中存在明显的 应力区，为设计时合理布置钢筋提供了依据。 以上几个算例说明，按目前的岩土力学水平，要使计算误差低于 10% 尚是奢望，但是，如果排除施工 条件与原计划相差过远的情况，预测的位移误差控制在 50% 以内是可以达到的，应力的误差则可能更大 一些，但是对工程造成危害的 应力区的位置则大体上是可信的。 3、复杂问题的数值模拟 实践表明，目前的计算岩土力学水平还只能应付简单的加荷过程和简单的研究对象，相当多的实际 问题尚无法进行数值分析，其中某些问题的数值模拟随着研究的深入在近期可能取得重要的进展。下面分 两类问题作一些说明。 1、复杂加荷条件 1 （）应力路线转折和主应力轴旋转，其中相对简单的如水库蓄水期的主应力轴偏转问题目前已取得一些 进展，但风浪荷载下海床中主应力轴旋转而引起的砂土液化问题则尚在摸索之中。 2 （）反复荷载下土的液化问题，目前只能作粗略的估计 （平均过程理论），如何进行精细的模拟（波动 过程理论）尚在探索之中。 3 （）爆炸等动力荷载下的砂土液