第二章 土本构关系.ppt

第二章 土的本构关系及土工有限元分析 有限元法的突出优点是适于处理非线性，非均质体和复杂边界等问题，而土体应力变形分析恰恰就存在这些困难问题，因此很适宜用有限元法。自从1966年美国克拉夫和伍德沃德首先用有限元法分析土坝以来，有限元法在岩土工程中的应用发展迅悉并取得了巨大进展。 有限元法，是用有限个单元体所构成的离散化结构，代替原来的连续体结构，来分析应力变形。这些单元体只在结点处有力的联系。材料的应力应变-关系可表示为： 由虚位移原理可建立单元体的结点力与结点位移之间的关系，进而写出总体平衡方程： [K], {d}, {R}分别为劲度矩阵、结点位移列阵和结点荷载列阵。把荷载作用于结点，由上式可求得位移，进而求出应变和应力。 第一节 土体的变形特性 土体的变形特性?土的本构关系 用压缩仪、三轴仪、平面应变仪、真三轴仪等试验手段，得出土的应力-应变关系。 试验的局限性：在研究实际问题中，土体各点的受力状况、变形历史，是千变万化的，无法在试验中模拟所有这些变化。 在试验的基础上，提出某种数学模型把特定条件下的试验结果推广到一般情况。这种数学模型就叫做本构模型。 Shaking table tests on the stability of an underground utility conduit in a liquefying sand deposit （Input motion：180gal, 10Hz, 10sec） 三轴试验测得的轴向偏应力s1-s3与轴向应变ea的关系曲线有两种形态。 图2-8(a)--p48所示曲线有一直上什的趋势直至破坏，这种形状的应力应变关系叫硬化型。软土和松砂表现为这种形态。 图2-8(b)所示曲线前面部分是上升的，应力到达某一峰值后转为下降曲线，即应力在降低，而应变却在增加，这种形态叫软化型。紧密砂和超压密粘土表现为这种形态。 岩土材料存在较大的塑性变形。沿不同的应力路经加荷，各阶段的塑性变形增量不同，累积起来，就有不同的应变总量。换句话说，尽管初始和终了的应力状态相同，加荷的应力路径不同，结果变形是不同的。 应力历史，是指历史上的应力路径。由于塑性变形不可恢复，历史上发生的变形将保存和积累起来。它无疑会影响今后的变形。 一般的本构模型理论中，常常用到p-q两个分量来构成二维应力路径，这种p-q平面上可以清楚地表示应力路径，但这些本构模型是忽略第三应力不变量或应力罗德角对变形的影响。 说明s2在s1和s3之间的位置，可以用如下指标：罗德参数ms、罗德角qs 、参数b 。 土体在高围压下的变形性状与低围压情况下有所不同。主要有如下三个方面： (1)强度包络线不呈直线，而是呈向下微弯的曲浅，如图2-14所示。 (2)前面讲到，有些土如紧密砂，受剪时体积会发生膨胀，那时在低围压下的性态。而在高围压时，所有土都表现为剪缩。图2-15所示。 (3)软化现象一般也是在低围压下。在高围压下，通常是硬化型的。如图2-15所示。---P52 地基土一般是水平向成层。水平和竖直方向土的结构存在差异，应力应变关柔也不例外。?原生各向异性。 应力状态不同，引起新的各向异性。 各向异性反映到本构关系式上，就是刚度矩阵[D]或柔度矩阵[C]为非对称矩阵。 第二节 弹性非线性模型 第三节 弹塑性模型 第四节 非线性有限元分析 第五节 土体非线性分析中的几个问题 * * Flac3D有限差分法 Plaxis有限元法 土体的应力-应变关系叫本构关系，是非线性的。所以矩阵[D]就不是常量，而随应力或应变改变，由此推得的劲度矩阵[K]也随应力或变形而变。 固结压缩试验仪 三轴压缩试验仪 GDS非饱和三轴测试系统 (一)土体的非线性和非弹性 应力-应变关系是非线性和非弹性 非线性和非弹性---区别：非线性材料有可能是弹性的；非弹性材料几乎没有是线性的。 传统意义上土的动应变和静应变---区别： 动模量---静模量的定义的模糊说明 总应变=弹性应变+塑性应变 　「Stable property」 versus 「Variable property*」 *) due to changes in : a) dry density; deposition condition; degree of saturation; so on b) effective confining pressure and strong effects of recent stress-strain history Zero confining pressure Stress Strain High confining pressure Low