双层地基一维非线性固结DQM解.pdfVIP

水 利 学 报 2004 年4 月 SHUILI XUEBAO 第4 期 文章编号：0559-9350 (2004) 04-0008-07 双层地基一维非线性固结的DQM解 王宏志，陈仁朋，周万欢，陈云敏 (浙江大学 岩土工程研究所，浙江 杭州 310027) 摘要：基于Davis和Raymond一维非线性固结理论，利用DQM(Differential Quadrature Method)推导了初始有效应力沿 深度变化、任意边界条件、任意荷载作用下双层地基一维非线性固结的表达式，求得了孔压、有效应力和平均固结 度的解答。解的收敛性分析表明，当每层土采用7～9个插值点时，DQM解就可以达到足够精度。而对于无法得到解 析解的复杂情况，DQM解具有统一的矩阵形式，有利于计算机程序的实现和求解，符合工程应用的要求。最后在算 例分析中，对考虑初始有效应力沿深度变化、矩形波循环荷载下双层地基一维非线性固结进行了探讨。 关键词：双层地基；一维非线性固结；DQM解；循环荷载 中图分类号：TU43 文献标识码：A 由于实际地基的成层性和非线性，双层地基的一维非线性固结研究具有重要的工程意义。20世纪60年 代，人们开始研究单层地基非线性固结问题。Davis和Raymond[1]假定e～log σ’成线性关系、固结系数在 固结过程中不变以及初始有效应力沿深度不变，获得了荷载恒定条件下的一维非线性问题的解析解答。 Barden等[2] [3] 及Mesri等 分别采用线性的e～log σ’和e～logk 关系，通过有限差分法得到了固结曲线。Xie v 等[4]在Davis和Raymond的研究基础上获得了荷载随时间变化的一维非线性固结的解析解。为了考虑地基的 成层性，李冰河等[5，6]将地基划分为若干薄层，同时在时间上进行离散，并假定每一薄层在每一离散时间 段内的固结参数为常数，得到了一维非线性固结的半解析解。Xie等[7]基于Davis和Raymond固结理论，求出 了双层地基一维非线性固结在瞬时加荷和单步线性加荷条件下的解析解。 Bellman等[8]在1971年提出了一种快速求解初值和/或边值问题的数值方法，称为DQM(Differential Quadrature Method)。它的主要思想是用空间上全部离散点的函数值的加权和近似表示每个点的各阶导数 值，经过这种DQM变换后，求解与时间无关的偏微分方程将转化为求解一个代数方程组，而求解与时间有 关的偏微分方程则转化为求解一个常微分方程组。由于DQM采用高阶的插值函数并且使用全部离散点的信 息拟合各点的导数值，它比传统的差分法和有限元法具有更高的计算效率[9]。目前这种方法已经广泛应用 于生物科学、热传导、流体力学、土力学、结构静力和动力学等领域，已有百余篇文章讨论及应用DQM于 物理和工程问题[10]。应用表明该方法具有数学原理简单、计算效率高，计算量小和结果准确的优点，是一 种能够部分取代传统的差分法和有限元法的数值方法[10] [11] 。但在岩土工程方面应用较少，Wang等 将DQM应 用于双层地基的一维线性固结问题中，得到了令人满意的效果。 本文基于Daivs和Raymond一维非线性固结理论，推导了初始有效应力沿深度变化、任意边界条件、任 意荷载作用下双层地基一维非线性固结的DQM表达式，得到了孔压、有效应力和平均固结度的解答。然后 通过解的收敛性分析讨论了本文方法的有效性和优点。最后利用本文方法，对考虑初始有效应力沿深度变 收稿日期：2003-01-09 基金项目：国家自然科学基金资助项目 作者简介：王宏志(1967-)，男，江西于都人，现在浙江大学土木工程博士后流动站工作，主要研究方向为桩基动力学， 计算土力学，土动力学以及结构的瞬态动力响应计算。