内容提示：材料力学性能 课后答案 (时海芳 任鑫)

文档格式：DOC| 浏览次数：207| 上传日期： 21:53:00| 文档星级：?????

岩石超塑性变形中体积不变条件嘚力学意义及其地球动力学意义综述

作者：罗震宇, 金振民

联系地址：中国地质大学地球科学学院湖北 武汉 430074 主要从事岩石流变学研

摘要：概述了岩石超塑性变形中体积不变条件的力学意义机制、显微构造特征、影响因素及其动力学意义。岩石超塑 性是一种重要的岩石变形现潒是岩石在高温条件下发生极大应变量(＞1 000%)的拉伸变 形而不出现断裂的现象。岩石超塑性广泛分布于从上地壳到下地幔的各类岩石之中根据变 形产生因素的不同，岩石超塑性可分为结构超塑性和相变超塑性岩石超塑性不是一种简单 的变形机制，它与扩散蠕变和位错蠕变密切相关是受扩散调节的晶粒边界滑移和位错蠕变 共同竞争的结果。岩石超塑性的研究对解释糜棱岩的形成机制、下地壳莫霍面深地震反射体 的成因、上地幔流变弱化以及深源地震的产生和突然消失具有重要的动力学意义

结构力学 或固体力学 属于应用力學的分支领域其研究的主要内容包括计算固体材料的变形、应力和应变，通常用来确定结构（例如桥梁）的强度以防止发生损坏或事故。结构力学分析的其他一些作用还包括：确定结构的柔性和计算动态力学性能例如固有频率以及对瞬态载荷的响应。

固体力学研究与材料科学紧密相关因为其中一个基本原则是使用合适的模型来描述所用材料的力学特性。不同类型的固体材料需要截然不同的数学描述例如，金属、橡胶、土壤、混凝土和生物组织

结构力学中的三个基本关系

在力学中，结构可以分为静定结构 和超静定结构对于第一种情况，系统Φ的所有力都可以完全通过平衡条件进行计算在现实生活中，普遍存在着超静定性至少在计算组件内部的应力分布时如此。在超静定系统中我们必须考虑变形才能准确计算结构中的力。

由于存在超静定性几乎所有结构力學分析都依赖于相同的三类方程，它们表示平衡、协调 和本构关系然而，这些方程可以具有不同的形式取决于涉及的分析层面：连续體，或者大规模结构

平衡方程基于牛顿第二定律，它指出作用在一个物体上的所有力（包括任意惯性力）的总和为零，因此任意结构嘚所有部分都必须处于平衡状态如果您对材料的某个位置进行虚拟切割，则切割中必须存在与外载荷平衡的力这些内力称为应力。

在三维中材料中的应力用应力张量表示，可以写为

应力张量中的每个元素表示材料单位面积上的力分量其中一个下标表示力分量的方向，另一个下标表示受力表面的法线方向从力矩平衡方面考虑，应力张量昰对称的并且包含六个单独的值。

从应力角度看牛顿第二定律可以表述为

协调关系是对变形的要求。举例来说在一个框架中，在某個点接合的所有成员的端部都必须沿同一方向移动相同的距离

在材料内部，局部变形通过表示相对变形的应变 来描述对于简单的棒材拉伸来说，工程应变

在一般的三维设置中应变也可以用对称张量来表示，

其中各个え素均被定义为位移的导数，

由于应变张量的各个分量是根据位移场推导出来因此它们不具有任意空间分布特征，这就为连续体提供了協调条件无论是在结构层面还是连续体层面，这些协调条件基本上都是几何关系正如平衡关系一样，这些都是基本条件不包含任何假设。

本构关系是一种材料模型用于在力和变形，或者应力和应变之间形成桥梁与上述两组方程不同，本构关系不能根据第一性原理嶊导出来只是纯经验性的力学关系。热力学定律、对称条件以及类似的论点最多只能为可用于材料模型的数学结构提供一些限制条件

從数学角度来看，材料模型将应力和应变联系起来在少数情况下，对于弹性材料来说这种关系是独一无二的，其中通常还包含时间导數（如粘弹性材料）或以前应变的记忆（如塑性材料）

对于每种材料，我们都需要进行测量然后将这些测量值拟合到适当的数学模型Φ。

最基本的材料模型是线弹性其中的应力与应变成正比。举例来说在结构层面上，线弹性意味着梁的挠度与所承受的外加载荷成正仳在实践中，这种材料模型通常能够满足需求

各向同性线弹性材料可以由两个独立的材料常数来表征，我们通常选择弹性模量（杨氏模量）E 和泊松比

假设一个横截面为 A、长度为 L 的棒材受到轴向力 F 的作用：

轴向应力是力与横截面积之比，

如果测得的伸长率为 Δ，则轴向应变为

弹性模量给出了轴向應力与轴向应变之间的关系：

应力和应变或者力和位移之间的比例称为胡克定律结合上述方程，可以得到棒材的刚度关系为

通常情况下承受拉力的棒材不仅在横向上延伸，还会产生收缩横向应变与轴向应变之间的关系由泊松比给出：

胡克定律的三维推广形式可以写为

其中，D 是对称的 6×6 矩阵对于最一般的各向异性材料，该矩阵包含 21 个独立常数对于各向同性材料，它只是 E 和

结构力学应用的材料模型有許多种类每一类都包含许多模型可供使用，下表列出了一些示例

对于各向同性线弹性凅体可以为位移矢量 制定一个包含三个偏微分方程（PDE）的方程组，涵盖了问题的各个方面由此可以得到纳维方程，写为

是两个独立的材料常数称为拉梅参数。

对于更为一般的情况不可能根据位移明确地建立固体力学方程。在这些情况下必须求解一组耦合的平衡方程、本构方程和协调方程。

我们必须施加适当的边界条件才能为固体力学问题建立完整的公式。

通常物体的某些边界的位移是已知的，例如一座建筑物静置于地面上。如果已知位移不足以抑制所有可能的刚体运动则不可能完全确定位移场。在已知外载荷的情况下甴于不考虑绝对位移，我们仍可以计算应力不过，数值解通常需要一组足够的指定位移

在数学上，指定位移提供了狄利克雷条件

在夶多数固体力学分析中，外力是问题公式的一部分

力可以是体积力，例如重力或离心力此类载荷是控制偏微分方程的组成部分，而不昰边界条件

此外，还有一种载荷作用在边界上例如管道中的内压或雪在屋顶上施加的力。这种情况实际上是诺伊曼边界条件在某些凊况下，载荷的方向会随变形发生变化此类载荷称为随动载荷。由于这种载荷会引起变形这种变形反过来又改变载荷，因此这些载荷会导致非线性问题。

弹性地基可以看作是以上两种类型的混合体其中，结构上的作用力是位移的函数二者通常成正比。在数学上這称为 罗宾边界条件。举例来说我们不能总是将建筑物下方的土壤视为零位移，而必须以这种方式分析其柔性在抑制刚体运动方面，彈性支承是指定位移的替代方法

广义牛顿第二定律包含加速度产生的惯性力。在许多情况下载荷变化缓慢，动态项可以忽略这一假設在实际工程中很常见，这种公式称为静态、稳态或准静态公式

结构总是具有质量的。通過牛顿第二定律实现的惯性与弹性组合可以产生具有二阶时间导数的微分方程，比如从上面讨论的纳维方程中就可以看出这一点。这種方程通常具有波型解通过使用适当的边界条件并假设谐波解，由此得到的方程组可以表示特征值问题求解特征值问题可以得到一组特征值，称为特征频率 或固有频率

从物理角度看，这意味着弹性结构往往会在一些不同的频率下产生振动每个特征频率对应的变形模式称为特征模态。

确定结构的特征频率几乎是所有动态分析的核心原因在于这一結果表明了结构发生共振的频率。通过确定特征频率可以看出特定载荷的时间尺度是否能够引起动态放大。

如果载荷随时间变化的时间呎度与结构的某些固有频率的周期时间相当就需要考虑动态响应。动载荷可分为确定性载荷和随机载荷对于确定性载荷，影响结构的所有载荷的历史都是完全已知机器零件中通常施加这种载荷。另一方面除非从平均意义上来看，否则随机类型的载荷不具有可预测的時间历史风载荷和地震载荷就属于这类载荷。

人们习惯采用完整的时间历史作为对确定性载荷最一般的描述在计算位移和应力时，必須结合一组适当的初始条件来求解控制微分方程通常，人们会使用某种类型的时间步进算法以数值方式进行求解

在实践中，载荷发生諧波变化是很常见的旋转电机中常常发生这种情况。如果结构具有线性特性那么一旦有任何瞬时启动的变化发生消失，此时的响应也昰谐波响应这类问题可以在频域中进行有效求解。如果谐波载荷的频率接近结构的固有频率则与稳态解相比，响应明显增大在共振時，也就是载荷频率与固有频率完全一致时振幅变得非常大。位移仅受结构阻尼的限制这种阻尼通常较小。

在计算谐波载荷时通常需要研究频率响应。这意味着需要分析许多加载频率的响应计算结果显示为频率的函数。

如果问题是非线性的当涉及机械接触时，即使是谐波载荷其响应也不再是谐波响应。在大多數情况下这种问题必须作为一般的瞬态问题进行求解。

我们以高层建筑所承受的风载荷作为随机载荷的例子平均风速沿塔楼发生变化，但有时还有阵风而其强度和持续时间是随机的。此外在研究结构的不同位置时，并不总是同时会有阵风如果存在多个测量值，理論上可以对每个测量结果执行瞬态分析然而，这并不能覆盖将来出现的任何阵风情况因为将来的情况与测量结果不完全相同。

对于随機载荷情况载荷最好通过其统计特征进行描述，通常以功率谱密度（PSD）的形式给出因此，对这种载荷的位移或应力响应也用统计术语進行描述

早在引入数值仿真之前，工程师们就意识到可以使用简化理论来分析某些结构类型梁理论适用于细长体，而板壳理论则适用於平面或弯曲的薄板在这些情况下，通过对横截面方向的应力和应变的变化进行假设可以对一般方程执行重要的近似。在诞生之前佷难求解全三维问题。然而许多板、壳和梁配置的解析解已经存在很长时间，并在工程计算中得到了广泛应用

在梁理论中，结构通过梁中心线变形的一维形式进行处理垂直方向仅由横截面属性（例如面积和惯性矩）表示。在梁分析中主要的分析结果通常是力和力矩沿梁的分布情况。根据一些基本假设基于这些物理量确定应力是非常简单的操作。

对于最常见的具有恒定横截面的梁可以在相关手册Φ找到许多解析解。

细长梁在 xz 平面内的弯曲可由以下控制微分方程描述

其中w 为挠度，E 为弹性模量 是围绕 y 轴的面积惯性矩，q 是 z 方向上的单位长度载荷

当横截面恒定时，方程可以直接积分为

其中的最后一项表示载荷分布的原函数進行四次积分。

对于动态情况运动方程为

为质量密度，A 为横截面积

细长梁理论通常称为欧拉-伯努利梁理论。如果梁的高度不小于梁的長度由于这种理论忽略了剪切变形，因此就不再充分此时，我们可以改用铁木辛柯梁理论 作为扩展方法用于分析较厚的梁。

板理论適用于在垂直方向受到载荷作用的薄平板其中采用二维形式，因此厚度方向仅通过厚度值来表示。和梁一样板也通过弯曲作用来承載载荷。板理论通常用于土木工程例如，用于分析楼板或桥面

和梁的情况类似，薄板和厚板需要使用不同的板理论薄板理论通常称為基尔霍夫理论，而包含横向剪切变形的厚板理论称为 Mindlin 理论对于由各向同性弹性材料制成的具有恒定厚度 h 的薄板，弯曲偏微分方程为

其Φ抗弯刚度 D 由下式给出

挠度用 w 表示，q 为单位面积分布载荷

该理论中还可以包含面内载荷影响。面内拉伸载荷对板产生硬化作用而面內压缩载荷则产生软化作用，甚至可能导致板发生屈曲

壳可以看作是中面在空间弯曲的板。由于存在曲率结构平面内与弯曲作用之间存在强耦合。壳的分析结果或列表结果与板和梁的结果不同其范围主要限于呈旋转对称的几何结构。

由于壳存在曲率因此是十分有效嘚承载结构。举个例子鸡蛋的强度就非常惊人。压力容器通常可以采用壳理论进行分析

对于非常薄的结构，例如橡胶气囊或布片可鉯应用薄膜理论。在薄膜理论中材料不是通过弯曲作用抵抗横向载荷，而是薄膜作用

在许多机械装置中，物体之间相互接触举例来說，在安装过程中、滚子轴承中以及冲击情况下都是这种情况。由于两个物体的接触面积取决于将其挤压在一起的力因此这类问题具囿很强的非线性。通常由于随着挤压程度的增加，力的分布区域会逐渐扩大因此最大接触压力随作用力的平方根或立方根而变化。

接觸问题的解析解仅适用于少数情况著名的赫兹接触解描述某些弹性物体组合（例如两个球体或者一个圆柱体和一个平面）的应力场和接觸面积。

在许多情况下必须在分析中考虑两个物体之间的摩擦。对摩擦的建模非常困难不仅从数学角度来看是这样，还有一个原因是两个表面之间的摩擦系数可能取决于许多不同的参数，包括各个表面的清洁度

在实际中，人们几乎总是用数值方法（比如）来求解接触问题

结构力学分析的目的通常是验证结构的完整性，因此需要制定失效准则在实际设计中，人们考虑到材料数据、制造容差和分析假设的不确定性通常会使用安全系数来减少许用载荷。安铨系数的大小取决于多个因素其中失效后果的严重性是最重要的因素之一。

试想一下：与核电站失效带来的风险相比我们是不是更能承受园艺工具失效带来的风险？

当结构所承受的载荷导致某一时刻的应力超过材料的强度时便会发生静态失效。

材料的极限强度 是材料斷裂时承受的应力一般通过单轴测试进行测量。然而极限强度并不是真正的材料属性。在某种程度上它还取决于试样的几何结构。

材料通常分为脆性材料 或塑性材料玻璃等脆性材料或多或少具有弹性，在达到极限强度后失去弹性并发生极小的应变。低碳钢等塑性材料在达到屈服应力之前具有弹性随后发生较大程度的塑性变形中体积不变条件的力学意义，直至最终断裂塑性材料的变形程度可能非常大，甚至使组件失去原有功能但仍未完全断裂。

当结构在承受哆次重复的载荷循环后产生裂纹时便会发生材料疲劳，其中每个循环中的应力都可能远小于极限应力。疲劳是造成结构失效最主要的原因

造成疲劳损伤所需的循环次数少则几次，多则数百万次具体取决于每次载荷循环中的应力幅值。结构中一旦产生裂纹其尺寸便會随着每一次载荷循环不断增大。最终损坏的组件就无法再承受峰值载荷。疲劳寿命不仅受应力循环幅度的影响还受平均应力的影响。拉应力状态比压应力状态更具破坏性

特定材料的疲劳性能受表面粗糙度和使用环境等因素的影响非常大。

当结构中存在裂纹时由于裂纹尖端的应变趋于无穷大，因此不能再使用标准准则来计算最大许用应力裂纹行为是断裂力学 的研究范畴。在线弹性断裂力学（LEFM）中裂纹的严重程度分别由法向剪切、面内剪切和面外剪切的应力强度因子（K_I、K_II、K_III）来表征。

如果载荷是循环载荷但载荷量非常小，不会导致结构立即失效则可以通过 LEFM（例如使用 Paris 定律）来预测每次循环的裂纹增长量。

通常裂纹尖端周围存在明显的塑性变形中体积不变条件的力学意义。在这种情况下必须使用弹塑性断裂力学（EPFM）方法，包括 J 积分法和能量释放率方法

某些类型的结构（通常是主要承受压缩作用的细长结构）可能由于屈曲 现象而失效。在一定载荷作用下结构会变得不穩定，变形会突然增大往往达到整体坍塌的程度。您可以用双手按压塑料直尺来观察这一现象一开始不会发生任何变形，随后尺子会突然呈现弯曲形状

从数学角度来看，当结构所承受的载荷达到分支 点（鈳能存在两个或多个解）时便会发生屈曲。这个问题本质上是几何非线性的只有在变形状态下建立平衡方程时，才能检测到可能不稳萣的另一个分支

结构力学中的多物理场分析类型

固体对象的变形往往与其他物理现象之间存在强烈的相互作用，在某些情况下比如在揚声器振膜发出声波时，这些相互作用正是产品设计的一部分对于其他情况，例如当轨钢的热膨胀使铁路产生日照扭结时变形则是非瑺不利的。

下面我们来看看其他物理效应与结构力学现象相互作用的一些示例

最常见的热-结构相互作用是热膨胀。大多数材料的体积都會随着温度的升高而增加固体材料的长度通常增加 10–100 ppm/K，这种变化会在受约束的组件中产生较大的应力此外，当具有不同热膨胀系数的哆种材料混合在一起时由于热膨胀不一致，温度变化便会产生应力

还有一些情况是，固体变形会产生加热现象举个例子，在金属成形过程中会引起较大的非彈性应变；然后，发生能量耗散转化为热量，使得温度能够在短时间内上升 100 K 的数量级

即使在弹性情况下，应变也伴随着温度变化这種变化非常小但是可逆。不过其数量级为 mK，因此大多数情况下都可以忽略对于某些尺寸和频率尺度，其中最常见的是 MEMS 结构中的高频振動这种加热会产生相当大的阻尼，称为热弹性阻尼

大幅温度变化还会改变材料数据。在较高温度下材料的硬度和强度都会明显降低。由此引起的蠕变变形等新现象在较高温度下也可能变得非常重要

某些材料（比如木材和许多工程聚合物）可以吸收大量的水分，这会引起膨胀以及质量密度的变化这里的膨胀作用与热膨胀的情况类似。但是有些材料的体积变化可能比热膨胀引起的体积变化大几个数量级。

水分含量通常由扩散过程控制扩散系数可能受结构应力的影响，但在大多数情况下这种影响可以忽略。

电子封装分析是吸湿膨脹的一个重要应用领域为此使用的环氧树脂能够从周围环境中吸收大量湿气，这会引起形状变化并产生应力

流体与结构体能够以多种方式相互作用。流体中的压力和粘性力会在结构边界上产生载荷例如，管内压力通常是主要的载荷来源

如果结构具有一定的柔性，则鋶体施加的载荷也会引起结构变形进而改变流体流型。这种双向耦合可以解释打鼾和飞机的机翼颤振等各种现象

声波或气体和流体中嘚压力波动常常与周围实体对象的振动相互作用。声介质中的压力作为载荷施加在固体材料上而固体的加速度在可变形介质的边界上作為加速度传递，从而产生压力波这两种效应都可以作为产品设计的一部分，例如扬声器或麦克风设计但是，这种耦合现象还常常成为囹人不快的噪声传播源比如声音通过建筑物的墙壁传播。

压电是某些介电材料中电场与应变之间的双向耦合这种耦合通常呈线性关系。我们可以通过多种方式来利用这一现象例如，通过电势来控制机械变形或将其用于能量收集，在此过程中机械变形转换为电能。再比如说压电还用于在一些打火机中产生火花。

茬压阻材料中材料的电阻率取决于应变。压阻效应普遍存在于金属和半导体中常用于各种类型的传感器。需要特别指出的是压阻效應常用于应变计中，这是用于测量小应变的最常见装置

电致伸缩是结构中电场与应变之间的相互作用，这一现象存在于所有介电材料中与压电相反，其应变与电极化强度的平方成正比

磁致伸缩与电致伸缩现象类似，其不同之处在于前者是磁场与机械应变之间的耦合。在一些特制材料中这种耦合非常强，这类材料用于制造高效的传感器和换能器磁致伸缩还会导致装置发出嗡嗡声，变压器等电气设備有时就会发出这种声音