等效应力原理.pptVIP

模型的位移边界条件为：下边界固定，上边界自由，左右边界只有y方向的位移。 应力边界条件为：左右边界剪应力均为0。模型上边界受相当于渠道衬切板重量和与衬切板垂直的冻结力的均布压力P=1400N，方向垂直于边界线。计算截面模型跟前面热传导问题的模型一样，截面的有限元网格图如下： * x方向位 移分布 取渠道的初始温度分布为冰的弹性模量关系式中的T值，计算得到渠道上边界受均布外载作用时内部位移场和应立场分布。如果温度取不同冻结时刻的温度值，则应力场可与温度场相互耦合，得到不同冻结时刻渠道内部的温度场和位移场、应立场分布。 * x方向应 力分布 y方向位 移分布 * 剪应力 分布： y方向应 力分布 * 3.7 小结 从复合材料的细观机理出发，利用土力学原理，先后建立了土、冰耦合的冻土本构模型，土、冰、水耦合的冻土本构模型，并对温度影响下的冻土本构关系进行了初步探讨，建立了基本分析模型。另外，用飞箭软件对一些前人所作的实际问题进行了数值模拟，为编制适用于我们提出的本构模型的计算程序做了前期准备工作。 * 进一步完善和修正饱和冻土的本构模型，如可以 考虑粘弹性、塑性、粘塑性等不同的本构模型。 ??从细观的角度研究非饱和冻土的力学特性和本构 模型。对土体冻融过程中未冻水和冰、气之间的 相变机理进行研究。 ? 用我们推导出的饱和冻土、非饱和冻土的本构模 型，用有限元法对一些工程实际问题如路基、隧 道、地基等进行数值模拟，以验证我们模型的合 理性和可靠性。 四、下一阶段的工作计划 * * ， 根据 Lemitre 等效应力原理，受损冻土在一维状态下的应力应变关系为： 冻土材料内部损伤型本构关系为： 记 则上式化为： * 3.2 土、冰混合本构关系的细观分析 3.2.1 Eshelby等效夹杂理论 3.2.3 微分法 * 3.2.4 Mori-Tanaka方法 3.2.5 混合率方法 Voigt等应变假设： Reuss等应力假设: * 3.2.6 几种弹性模量的比较 1－Eshelby等效夹杂方法 2－微分法 3－Mori-Tanaka方法 4－Reuss等应力假设 1－Eshelby等效夹杂方法 2－微分法 3－Mori-Tanaka方法 4－Reuss等应力假设 5－Voigt等应变假设 图1 不同细观方法得到冻结砂 土弹性模量随温度变化图 图2 不同细观方法得到冻结砂 土弹性模量随温度变化图 * 3.3 与实际冻土的比较 冻结砂土的弹性模量随温度变化曲线 不同的冰体积含量 ，由式 计算得 到砂土冻结后的等效弹性模量随温度变化曲线： * 图a的应变速率为 1.1×10-4 (s-1 )， 1.1×10-3 (s-1 )。 实验曲线的围压为： 图b的应变速率为 冻结砂土的应力－应变曲线及其比较（ =0.05, =0.5） 1－ =8Mpa， 2－ =12MPa， 3－ =22MPa * 冻结砂土的应力－应变曲线及其比较（ =0.05, =0.8） * 冻结砂土的应力－应变曲线及其比较（ =0.1, =0.5） * 冻结砂土的应力－应变曲线及其比较（ =0.1, =0.8） * 冻结砂土的应力－应变曲线及其比较（ =0.2, =0.5） * 冻结砂土的应力－应变曲线及其比较（ =0.2, =0.8） * 3.4 饱和冻土的本构关系 借用土力学中饱和土体的Terzaghi有效应力原理，我们对饱和冻土体的Terzaghi有效应力原理表述为：饱和冻土体中任一点的总应力为该点有效力与孔隙水压力之和，即： * 3.5 温度对饱和冻土本构关系的影响 冻土体内一点处的水压、冰压以及温度之间的关系用Clapeyron方程描述： 水压、冰压与该点的孔隙正压力 的关系为： 水压力 为： * 3.6 数值模拟 3.6.1 飞箭有限元程序自动生成系统 FEPG 功能特点：FEPG 可自动产生计算单元刚度矩 阵、质量矩