[中医中药]5_6课-计算土力学.ppt

计算土力学 主讲教师：张爱军 §2.7 岩土工程基本分析方法 分类 总应力分析法（p＝0） 初期变形分析 使用不排水指标，计算出的结果是加荷瞬间或短期的应力与变形。 终期变形分析 采用排水指标，计算出的结果是最终或长期的应力和变形，即：孔压消散完毕，主固结完成后最终的应力和变形。 有效应力分析法（p≠0） §2.7.1 总应力分析 将土体视为固体，不考虑土中的渗流，不区分单元中分别由土骨架和孔隙水压力传递和承受的应力，也就是不考虑有效应力与孔隙水压力，只考虑土体总应力。其分析基本理论与固体力学中无任何区别 适用于不考虑渗流固结的情况，例如：饱和粗粒土地基、强透水土料组成的土坝路堤的应力和变形，以及饱和软粘土地基短期变形和稳定问题。 其基本方程： 平衡方程 几何方程 物理方程 孔隙水渗流连续方程 孔隙水平衡方程 我们以静力问题为例说明总应力分析的基本方程。同时讲授基本方程的矩阵表示方法。 平衡方程 其矩阵表示为： 物理方程 几何方程 总体控制方程就是： 将： 与 代入： 得到： 这就是总应力法静力问题控制方程 这个方程实际上是Biot动力固结方程中的惯性力为0，不考虑孔压的特殊形式。 对于各向异性、横观各向同性体其实就是[D]阵不一样而已，没有什么复杂的，总体控制方程的矩阵形式是一致的。 但是对于{R}阵，对于有面力的情况和有初始应变或应力的情况，还有不同的形式，以后的章节中再讲，但是整体控制方程的矩阵形式是一样的。 对于动力情况以上控制方程是不适合的，因为存在惯性力项，另外物理矩阵中还应该包含粘性的成分，即： §2.7.2 有效应力分析 有效应力法要严格区分土骨架变形和孔隙水压力，更能够真实地反映土体的自身特性。 适用于：黏土地基、坝体和路基随时间固结沉降分析，特别适合于求解随时间变化的过程分析，反映施工过程。 其控制方程就是Biot静力固结方程 平衡方程 由于平衡方程只涉及总应力因此与总应力法是一致的，即： 几何方程，物理方程也是一致的，即： 有效应力原理：表示成为矩阵形式： 将物理方程、几何方程、有效应力原理代入平衡方程中，得到： 孔隙流平衡方程，若略去惯性力就为Darcy定理，即： 变成矩阵形式为： 孔隙流体连续方程为： 将Darcy定理应用于连续方程中得到： 这样就得到静力有效应力控制方程： 其中：未知量有：{u}, p 还包含了时间项 t §2.7.3 初始条件和边界条件 有了总体控制方程是无法求解的，还必须有初始条件和边界条件 对于总应力法其未知量为三个方向的位移，因此只涉及位移或应力初始条件和边界条件。对于有效应力分析中还存在孔隙水压力未知量，因此还增加渗透边界条件。 我们以有效应力分析为例说明边界条件和初始条件问题，总应力分析中只要去掉渗透边界条件即可。 边界条件 位移边界条件 位移边界条件分为二种，即：约束边界和已知位移边界。其中约束边界就是在某个边界上某个自由度上位移受到约束，在这个自由度上的位移为0。如图所示： 例子： 另外一种是已知位移边界，即： 应力边界 就是指在某个边界上存在面力。 已知水压力边界，也就是在某个边界上孔压或水头已知。 透水边界： 水头已知边界： 已知流量边界，也就是在某个边界上渗透流量是已知的。 不透水边界 输入流量边界 设：边界上沿外法线方向的已知流速为vn，则流速边界条件可表示为： L,m,n为边界面外法线方向与xyz坐标轴的夹角的余弦（方向余弦）。 那么：根据达西定理： 对于不透水边界为： 对于已知流量边界为： 初始条件 初始位移条件 初始应力条件 初始孔隙水压力或水头条件 初始流量条件 这些条件是进行迭代计算的基础。 §2.7.4 有效应力法与总应力法的关系 从控制方程可以看出：总应力法实际上是有效应力法当p＝0时的特殊形式，编制一个有效应力程序是可以用作总应力计算的。 有效应力法是土力学区别于固体力学的主要方面，作为一个土力学工作者，必须掌握。这也是我们可以从事固体力学有些问题研究，而固体力学人从事不了我们工作的主要原因。 总应力法作用也是很大的，也是目前工程上用的较多的，也应该重视。 总应力法应用问题 谢康和书中“TT”法，“Total stress analysis with Terzaghi’s consolidation theory” 。就是用总应力法采用排水指标计算土体的最终变