等效应力原理摘要.pptVIP

* * ， 根据 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 可自动产生计算单元刚度矩 阵、质量矩阵、阻尼矩阵和单元荷载向量的单元子程序。算法程序可按照用户给出的代数方程组的矩阵表达式及右端项表达式，自动产生算法主程序。 飞箭有限元程序自动生成系统 FEPG (Finite Element Program Generator)是一个开放、灵活的有限元程序开发平台，它可以根据用户输入的有限元方法所需的各种表达式和公式，自动产生解决该问题所需的全部有限元计算的 FORTRAN 源程序。 2. 适用领域：FEPG系统可以提供静力分析、热 传导分析、电磁场分析、流场分析、结构力学 分析、渗流场分析、围岩稳定性分析、高边坡 稳定性分析、蜗壳结构受力分析、岔管结构分 析、水力学计算等。 3. 前、后处理系统：前后处理器FEPG.GID通过 体、面的布尔加、减、交等操作得到模型，可 快速将几何模型自动离散成各种单元，并可根 据用户需要对网格进行局部加密。后处理支持 的结果显示方式有：带状云图显示、等直线显 示、切片显示、矢量显示、变形等。 3.6.2 数值模拟实例 用飞箭软件对饱水渠道基土冻结过程中的热传导问题进行数值模拟，得到二维渠道基土冻结过程中的温度分布。 假定渠道无限长，截面为平面问题，认为土质均匀、水热蒸发热量及其它势场忽略不计。渠道模型下边界宽3.0m；渠底距下边界1.0m，宽0.5m；渠堤距下边界2.0m,宽1.0m。计算截面示意图如下： 其中： 1．定解问题 初始条件和边界条件为： 上述各式构成渠道基土二维相变温度场的定解问题 2．参数的确定 3．数值计算结果 计