岩石的弹塑性本构关系.pptx

第三章第三章岩石的弹塑性本构模型岩石的弹塑性本构模型 岩石的弹塑性本构模型岩石的弹塑性本构模型前言前言 张量知识简介张量知识简介线性弹性理论线性弹性理论非线性弹性理论非线性弹性理论 应力空间表述的弹塑性本构理论应力空间表述的弹塑性本构理论前言前言1、条件：将岩石介质看作成一种连续介质、条件：将岩石介质看作成一种连续介质宏观分析。宏观分析。2、主要研究内容：、主要研究内容：研究材料固有的特性，建立应力研究材料固有的特性，建立应力应变及与应变及与温度之间关系的表达式（本构关系）温度之间关系的表达式（本构关系）分析弹塑性变形体内应力、应变分布，研究分析弹塑性变形体内应力、应变分布，研究物体在各种荷载下的稳定性问题物体在各种荷载下的稳定性问题求解边值求解边值问题或求解初值问题或求解初值-边值问题。边值问题。应用数学问题，探求各种解析方法或数值方应用数学问题，探求各种解析方法或数值方法。法。本章着重解决本构关系。本章着重解决本构关系。张量知识简介张量知识简介 一、预备知识一、预备知识1.坐标系坐标系 2.约定求和约定求和3.克罗尼克尔符号克罗尼克尔符号4.置换符号置换符号二、张量的定义二、张量的定义1.坐标变换坐标变换 2.零阶张量（标量）零阶张量（标量）3.一阶张量（向量）一阶张量（向量）4.二阶张量二阶张量1.坐标系坐标系1）直线坐标系）直线坐标系 由坐标原点与三条不共面标架由坐标原点与三条不共面标架直线构成直线构成2）仿射坐标系）仿射坐标系 各标架上单位尺度不同各标架上单位尺度不同3）笛卡尔坐标系）笛卡尔坐标系 各标架上单位尺度相同各标架上单位尺度相同 笛卡尔直角坐标系笛卡尔直角坐标系 标架直线互相垂直标架直线互相垂直 笛卡尔斜角坐标系笛卡尔斜角坐标系 标架直线不互相垂直标架直线不互相垂直以以 表示笛卡尔直角坐标系的坐标表示笛卡尔直角坐标系的坐标 ，分别表示三个坐标的单位矢量。，分别表示三个坐标的单位矢量。2.约定求和约定求和定义：如果在同一项中，某个指标重复出现两定义：如果在同一项中，某个指标重复出现两次，就表示对这个指标从次，就表示对这个指标从1到到3求和。求和。例1 例2 例例3:例4 3.克罗尼克尔符号克罗尼克尔符号定义定义：故有故有：例例1：在笛卡尔直角坐标系中：在笛卡尔直角坐标系中：例例2：单位矩阵可表示为：单位矩阵可表示为例例3：例例4：例例5：例例6：4.置换符号置换符号 1）定义）定义：即：即：例例1：例例2：2）和和 的关系的关系 练习练习将下式写成常见的应力将下式写成常见的应力-应变方程组（广义虎克应变方程组（广义虎克定律）定律）例：将（例：将（2-1）式转换为（）式转换为（2-2）式）式由由得：得：（2-1）式可以写成）式可以写成两端乘两端乘调整参数调整参数 或或代入（代入（2-3）式，得）式，得整理得：整理得：二、张量的定义二、张量的定义张量是由满足一定关系的一组元素所组张量是由满足一定关系的一组元素所组成的整体，元素的个数由空间的维数成的整体，元素的个数由空间的维数N及及张量的阶数张量的阶数n决定，我们以决定，我们以N=3为代表，为代表，给出各阶张量的定义。给出各阶张量的定义。1笛卡尔直角坐标系间的坐标变换笛卡尔直角坐标系间的坐标变换：式中：式中：是是 轴与轴与 轴夹角的余弦，轴夹角的余弦，为坐标轴的单位向量为坐标轴的单位向量即即：2零阶张量（标量）零阶张量（标量）有有30=1个元素，是坐标变换下的不变量，个元素，是坐标变换下的不变量，即即3一阶张量（向量）一阶张量（向量）有有31=3个元素个元素 ，它们随坐标系，它们随坐标系变化的规律为变化的规律为或或该该3个元素组成的整体称为一阶张量，记作个元素组成的整体称为一阶张量，记作 称为称为 的分量，记作的分量，记作一阶张量一阶张量=向量向量4二阶张量二阶张量有有 个元素，个元素，它们随坐标系的变化规律为它们随坐标系的变化规律为或或则由这则由这9个元素所组成的整体称为二阶张个元素所组成的整体称为二阶张量，记作量，记作5.应力、应变张量应力、应变张量 应力张量：应力张量：是二阶对称张量，是二阶对称张量，应变张量：应变张量：是二阶对称张量，是二阶对称张量，广义虎克定律：广义虎克定律：式中：式中：弹性系数张量，根据张量识别定理，弹性系数张量，根据张量识别定理，是是4阶张量，共有阶张量，共有34=81个分量。个分量。线性弹性理论线性弹性理论一、空力空间和应变空间一、空力空间和应变空间二、用二、用cauchy方法给出的本构方程方法给出的本构方程 一、空力空间和应变空间一、空力空间和应变空间1应力空间：应力空间：定义：以应力分量定义：以应力分量 作为笛卡尔坐标作为笛卡尔坐标系中的坐标轴所形成的空间（一般概念上来说系中的坐标轴所形成的空间（一般概念上来说应是一个应是一个9维空间），岩体中某一点的应力状维空间），岩体中某一点的应力状态可用应力空间中的一点（坐标）来表示。态可用应力空间中的一点（坐标）来表示。常用应力空间常用应力空间9维空间不直观。通常采用三维维空间不直观。通常采用三维二维空间二维空间i)用三个主应力来表示：用三个主应力来表示：ii)用二个主应力来表示：用二个主应力来表示：iii)用剪应力和平均应力来表示用剪应力和平均应力来表示用剪应力和平均应力来表示用剪应力和平均应力来表示 应力分解应力分解：式中：式中：为平均应力为平均应力用剪应力和平均应力来表示用剪应力和平均应力来表示 常用二维：常用二维：空间表示空间表示其中：其中：用剪应力和平均应力来表示用剪应力和平均应力来表示 常用三维：常用三维：空间表示空间表示 八面体剪应力八面体剪应力 偏应力第二不变量；与形变能偏应力第二不变量；与形变能V有关有关轴对称情况（常规三轴试验）：轴对称情况（常规三轴试验）：用剪应力和平均应力来表示用剪应力和平均应力来表示有限元计算中常用的应力空间有：有限元计算中常用的应力空间有：空间，空间，2应力路径（应力路径（stress route）1）定义：应力空间中用来表示应力状态变）定义：应力空间中用来表示应力状态变化历史的一条曲线。化历史的一条曲线。2）举例：）举例：不同应力空间中常规三轴加载条件下的不同应力空间中常规三轴加载条件下的应力路线应力路线3应变应空间与应变路径应变应空间与应变路径相同定义相同定义本构关系就是应力、应变这两种状态所本构关系就是应力、应变这两种状态所满足的数学表达式满足的数学表达式 二、用二、用cauchy方法给出的本构方程方法给出的本构方程 1Cauchy假设：在外力作用下，物体内各点的假设：在外力作用下，物体内各点的应力状态和应变状态之间存在着一一对应的关应力状态和应变状态之间存在着一一对应的关系。因此，弹性介质的响应仅与当时的状态有系。因此，弹性介质的响应仅与当时的状态有关而与应变路径或应力路径无关。关而与应变路径或应力路径无关。推论：推论：卸荷后，介质回到初始状态卸荷后，介质回到初始状态 应力、应变都是瞬时发生的，在时间上无先应力、应变都是瞬时发生的，在时间上无先后顺序后顺序 在应力空间和应变空间的各点之间构成一在应力空间和应变空间的各点之间构成一一对应的映射关系。一对应的映射关系。2弹性本构方程弹性本构方程1）弹性变形是各向同性的弹性变形是各向同性的用应变表示应力用应变表示应力 （1）其中：其中：为为4 4阶张量。阶张量。对于各向同性：对于各向同性：式中，式中，为某一标量，是两个独立的常数，为某一标量，是两个独立的常数，称为拉梅常数称为拉梅常数 用应力表示应变用应力表示应变2）弹性参数：）弹性参数：纯拉（或压）虎克定律中的扬氏弹性模量纯拉（或压）虎克定律中的扬氏弹性模量E纯剪切虎克定律中的剪切弹性模量纯剪切虎克定律中的剪切弹性模量G泊松比泊松比 E：纯拉（或压）纯拉（或压）其余其余 实验上：单向伸缩的虎克定律可以写成实验上：单向伸缩的虎克定律可以写成故：故：泊松比泊松比 纯剪切虎克定律中的剪切弹性模量纯剪切虎克定律中的剪切弹性模量G G 又又 故：故：3）不同表示的弹性本构方程）不同表示的弹性本构方程 用用E，表示的弹性本构方程表示的弹性本构方程 用用K-G表示的弹性本构方程表示的弹性本构方程 用用G、E、表示的弹性本构关系表示的弹性本构关系 用用 表示的弹性本构关系表示的弹性本构关系 上述各式中的上述各式中的E，G，K等弹性参数可采用等弹性参数可采用割线参数。割线参数。4 4）增量形式的弹性本构方程）增量形式的弹性本构方程上述各式中的上述各式中的E，G，K等弹性参数可采用等弹性参数可采用切线参数。切线参数。非线性弹性理论非线性弹性理论非线性弹性非线性弹性Duncan（邓肯）模型（邓肯）模型该模型在岩土工程计算中运用相当广泛。该模型在岩土工程计算中运用相当广泛。（一）弹性模量（一）弹性模量 （二）泊松比（二）泊松比（一）弹性模量（一）弹性模量1）基本假定：基本假定：由常规三轴试验，应力由常规三轴试验，应力应变关系为双曲线形式：应变关系为双曲线形式：或或 2）参数意义：）参数意义：式中：式中：为初始模量。为初始模量。3）割线模量：）割线模量：4）切线模量 5）用三轴试验参数表示的非线性关系）用三轴试验参数表示的非线性关系 不同不同 ，有不同，有不同 不同不同 ，有不同，有不同 ，摩尔摩尔-库仑准则：库仑准则：有：有：式中：式中：内聚力，内摩擦角内聚力，内摩擦角 实验参数实验参数 大气压大气压（二）泊松比（二）泊松比 1 1）基本假定：轴向应变）基本假定：轴向应变 径向应变径向应变 为双曲线关系为双曲线关系或：或：2）参数的物理意义）参数的物理意义 3）割线泊松比）割线泊松比 4）切线泊松比）切线泊松比 5）用三轴实验参数表示的非线性关系）用三轴实验参数表示的非线性关系 不同不同 ，有不同，有不同 其中：其中：G、F为试验参数为试验参数 割线泊松比割线泊松比 切线泊松比切线泊松比（三）讨论（三）讨论1）邓邓肯肯模模型型由由常常规规三三轴轴试试验验曲曲线线来来确确定定切切线线的的弹弹性性常常数数。常常规规三三轴轴试试验验 ，是是一一种种特特殊殊的的加加荷荷路路径径，只只要要假假定定材材料料是是各各向向同同性性的的，沿沿 方方 向向 增增 加加 所所 确确 定定 的的 弹弹 性性 常常 数数 可可 以以 用用 于于 方方 向向 荷荷 载载 增增 加加 的的 情情 况况，当当 然然 也也 适适 用用 于于 和和 同同时时增增加加的的情情况况。这这样样的的确确定定方方法法在在理理论论上上是是正正确确的的，具具有有增增量量弹弹模模和和增增量量泊泊松松比比的的物理意义。物理意义。2 2）常规三轴试验）常规三轴试验曲线曲线 的割线的割线斜率不具有全量（割线）弹模的物理意义，因斜率不具有全量（割线）弹模的物理意义，因全量侧压力不为全量侧压力不为0 0，且，且 不等于全部的不等于全部的 ，不可以用来确定弹性常数。，不可以用来确定弹性常数。3 3）该模型只考虑硬化、不能反映软化。）该模型只考虑硬化、不能反映软化。4 4）塑性全量理论的本构关系在本质上也是一种非线性弹）塑性全量理论的本构关系在本质上也是一种非线性弹性关系，该理论中假设在全量应力与全量应变之间存在着性关系，该理论中假设在全量应力与全量应变之间存在着一一对应的关系。因此，在岩土工程计算中，只要是在比一一对应的关系。因此，在岩土工程计算中，只要是在比例加载或偏离比例加载的条件下，使用的全量塑性理论模例加载或偏离比例加载的条件下，使用的全量塑性理论模型在本质上都是非线性弹性模型。在节理岩石的有限元分型在本质上都是非线性弹性模型。在节理岩石的有限元分析中，广泛使用析中，广泛使用zienkiewiczzienkiewicz等人提出的不抗拉模型和层等人提出的不抗拉模型和层状材料模型，在这些模型中，只要规定了超出库仑准则的状材料模型，在这些模型中，只要规定了超出库仑准则的应力点拉回到库仑面上的具体方式，则在应力与应变之间应力点拉回到库仑面上的具体方式，则在应力与应变之间也有一一对应的关系，因而它们可以被看作是一种全量形也有一一对应的关系，因而它们可以被看作是一种全量形式的非线性弹性关系。式的非线性弹性关系。岩石的弹塑性本构模型岩石的弹塑性本构模型前言前言 张量知识简介张量知识简介线性弹性理论线性弹性理论非线性弹性理论非线性弹性理论 应力空间表述的弹塑性本构理论应力空间表述的弹塑性本构理论应力空间表述的弹塑性本构理论应力空间表述的弹塑性本构理论一、岩石的单轴压缩实验曲线和某些基本概念一、岩石的单轴压缩实验曲线和某些基本概念二、屈服条件和屈服面二、屈服条件和屈服面 三、塑性状态的加三、塑性状态的加卸载准则卸载准则 四、塑性势理论和相应的本构方程（增量形式）四、塑性势理论和相应的本构方程（增量形式）一、岩石的单轴压缩实验曲线和某些基本概念一、岩石的单轴压缩实验曲线和某些基本概念1、稳定阶段稳定阶段;非稳定阶段，非稳定阶段，2弹塑性材料的特点弹塑性材料的特点 屈服屈服产生塑性变形；产生塑性变形；加载、卸载时遵循不同的规律；加载、卸载时遵循不同的规律；加工硬化，软化加工硬化，软化 破坏破坏要表述弹塑性本构关系，必须反应上述特点。要表述弹塑性本构关系，必须反应上述特点。二、屈服条件和屈服面二、屈服条件和屈服面1屈服条件或屈服准则一般概念屈服条件或屈服准则一般概念 应力状态达到某种组合时，产生屈服，出现应力状态达到某种组合时，产生屈服，出现塑性变形塑性变形 应力张量应力张量 内应力张量，塑性应变所对应的应力内应力张量，塑性应变所对应的应力 塑性应变张量塑性应变张量 n 标量的内变量，可以代表塑性功标量的内变量，可以代表塑性功W，塑，塑性体应变性体应变 或等效塑性应变或等效塑性应变若 弹性状态 塑性状态 不存在2初始屈服条件初始屈服条件3屈服面，应力空间屈服面，应力空间 继4理想塑性体理想塑性体屈服面大小和开状不发生变化，屈服面大小和开状不发生变化，5硬化硬化-软化后继屈服面软化后继屈服面屈服面大小和形状由于内变量的出现和发展而屈服面大小和形状由于内变量的出现和发展而变化。变化。如式如式 等向硬化等向硬化软化：对于均匀扩大（硬化）和均软化：对于均匀扩大（硬化）和均匀收缩（软化）：后继屈服面仅含有一个标量匀收缩（软化）：后继屈服面仅含有一个标量参数参数 n 式中：式中：H材料参数材料参数H是标量是标量 n 的函数。的函数。6随动硬化：屈服面大小和形状保持不变，随动硬化：屈服面大小和形状保持不变，在应力空间中平动在应力空间中平动 式中，式中，为屈服面中心的位置矢量。为屈服面中心的位置矢量。材料参数，可是材料参数，可是 n 的函数。的函数。7混合硬化混合硬化 4.式中：式中：5 5、加载时的本构关系、加载时的本构关系 6统一的本构关系：统一的本构关系：