第二章各向异性弹性力学演示文稿

第二章各向异性弹性力学演示文稿现在是1页\一共有59页\编辑于星期五各向异性与各向同性弹性力学的基本方程的差别差别在于:本构方程其它平衡方程,几何方程,协调方程,和边界条件等则完全相同.即用各向异性胡克定律代替各向同性胡克定律,这一代换将使力学计算及反映的现象十分复杂.现在是2页\一共有59页\编辑于星期五单元体应力及正负号规定如果作用面的外法线指向坐标系中相应坐标轴的正向,而应力分量也指向对应坐标轴的正向,则应力分量为正。当两个下标中,只有一个指向坐标轴的正向时,该应力分量就为负.yx作用在y面上的正应力作用在y面内x方向的剪应力z现在是3页\一共有59页\编辑于星期五静力平衡方程（3）X，Y，Z作用于微元体的体积力力要平衡！现在是4页\一共有59页\编辑于星期五几何关系(小变形)(6)变形要协调！三个独立的位移场即可以完全确定变形，而应变亦可以描述变形，它们之间满足以下关系！现在是5页\一共有59页\编辑于星期五本构方程（6）反映出材料的性质！与之间的关系现在是6页\一共有59页\编辑于星期五各向异性弹性力学问题需满足的基本方程与各向同性弹性力学一样,各向异性弹性力学有15个未知量15个场方程静力平衡方程（3）＋几何关系（6）＋本构方程（6）可以求解了吗？现在是7页\一共有59页\编辑于星期五给定力的边界条件(3)定解还需边界条件！现在是8页\一共有59页\编辑于星期五给定位移的边界条件(3)现在是9页\一共有59页\编辑于星期五各向异性弹性力学问题需满足的基本方程（另一组定解方程）与各向同性弹性力学一样,各向异性弹性力学有12个未知量12个场方程静力平衡方程（3）＋几何关系（6）＋本构方程（6）＋变形协调方程（3）现在是10页\一共有59页\编辑于星期五变形协调方程(3/6)只有三个是独立的，为什么？现在是11页\一共有59页\编辑于星期五以上的力学,几何,物理,以及边界条件诸方面构成各向异性弹性力学的基本方程,与各向同性弹性力学的区别在于物理方程.其它均相同现在是12页\一共有59页\编辑于星期五

弹性介质的本构关系均质弹性体的弹性性质坐标转换(应力应变及弹性系数转轴公式)弹性对称性——本构关系的简化正交异性材料弹性常数的物理意义现在是13页\一共有59页\编辑于星期五2.1弹性介质的本构关系

2.1.1弹性介质的本构关系

规定下标i，j与一维指标对应如下次序：(2－1)则（2－1）

的两式可以写成矩阵乘法的形式，第一式可以写作

现在是14页\一共有59页\编辑于星期五记作

可以理解为张量等式，，理解为应力张量和应变张量，L理解为弹性刚度张量；也可以理解为矩阵等式，，理解为应力列矢量和应变列矢量，[L]理解为弹性刚度矩阵。L与M具有Voigt对称性，因此矩阵L与M为9列9行的对称矩阵。

（2－2）现在是15页\一共有59页\编辑于星期五由于应力张量与应变张量都是对称张量。（2－2）式中的列矢量与的第4行与第5行相同，第6行与第7行相同，第8行与第9行相同。弹性刚度矩阵与柔度矩阵第4行、列与第5行、列相同，第6行、列与第7行、列相同，第8行、列与第9行、列相同。利用这种对称性，可以把应力张量与应变张量写成6个元素的“列矢量”相应的，L与M可写成6行6列的对称矩阵现在是16页\一共有59页\编辑于星期五也就是说，各列除去重复的元素，但第1、2、3列的元素的数值不变，而第4、5、6列的元素则乘以2。此时，张量运算与矩阵运算仍然一样，但失去了矩阵地对称性。现在是17页\一共有59页\编辑于星期五有的文献中定义应力“列矢量”为应变“列矢量”为

注意：，，就是剪切角，，。

现在是18页\一共有59页\编辑于星期五于是可以把弹性本构关系写成：

或

（2－3）（2－4）现在是19页\一共有59页\编辑于星期五容易导出矩阵C，s与L，M之间的关系为

现在是20页\一共有59页\编辑于星期五2.1.2弹性应变能密度

固体变形时，加在它上面的外力要做功。完全弹性体在等温条件下，当缓慢卸载后可以完全恢复其初始状态。因此，可以认为，外力功全部以能量的形式储存在弹性体内。这种能量称为应变能。通过对体微元的研究，可以得到弹性应变能密度：其中

(Voigt对称性)(Voigt对称性)现在是21页\一共有59页\编辑于星期五由线弹性可以得现在是22页\一共有59页\编辑于星期五2.2均质弹性体的弹性性质

对于均质弹性体，材料的性质与位置坐标无关。其应变位能是应变分量，，…，的函数，而且只取决于应变的最终值。从数学上说，是应变状态的单值函数，而且与积分路径无关，必是对应变分量的全微分，即：可得（2－5）现在是23页\一共有59页\编辑于星期五为了便于以后的讨论，给出的展开式（2－6）现在是24页\一共有59页\编辑于星期五2.3坐标转换(应力应变及弹性系数转轴公式)2.3.1斜面应力为了讨论过点A任意斜面的应力，在点A附近取一个四面体微元ABCD(图2－1)。

图2－1现在是25页\一共有59页\编辑于星期五斜面BCD的外法线为N，令N的方向余弦为：则有式中，、、、依次为三角形BCD、ACD、ABD、ABC的面积。令四面体微元的体积为dV，斜面BCD上应力向量在坐标方向上的分量为、、，则由四面体微元的的条件得到：（2－7）现在是26页\一共有59页\编辑于星期五得到方程如下：写成矩阵形式也就是说，若应力张量为已知，则任一斜面上的应力均可求出。因此，应力张量完全决定了一点的应力状态。（2－8）现在是27页\一共有59页\编辑于星期五2.3.2应力应变转轴公式三维情况：图2－2坐标系如图2－2所示，新坐标与原坐标的方向余弦列于表1

：其中现在是28页\一共有59页\编辑于星期五

xyzx′l1m1n1y′l2m2n2z′l3m3n3表1

即

（2－9）现在是29页\一共有59页\编辑于星期五将式（2－9）展开，并按一定次序排列应力张量，可得应力分量转轴公式：（2－10）称为应力转换矩阵现在是30页\一共有59页\编辑于星期五同理可得，应变分量转轴公式（2－11）称为应力转换矩阵现在是31页\一共有59页\编辑于星期五二维情况：二维情况的坐标建立如下两图：图2－3图2－4现在是32页\一共有59页\编辑于星期五

同理：（2－12）（2－13）（2－14）现在是33页\一共有59页\编辑于星期五2.3.3弹性系数的转轴公式各向异性体的弹性特性随方向不同而异，即各向异性体的弹性系数是方向的函数，它们与坐标的取向有关，只有在各向同性情况下，弹性系数在任意正交坐标系是不变的。

已知

求逆

又因为

所以可得现在是34页\一共有59页\编辑于星期五通过单位体积应变能函数U0可以证明：从而可以得到：所以其中为新坐标系中的柔度矩阵

为新坐标系中的刚度矩阵

以上即为弹性系数转轴公式的矩阵形式。

（2－15）现在是35页\一共有59页\编辑于星期五2.4弹性对称性——本构关系的简化

大多数工程材料内部结构具有对称性，这决定了材料的弹性特性也具有对称性。均质弹性体中，若过每一点的不同方向的弹性特征不同，称为一般均质各向异性体，它具有36个非零的弹性系数，但其中21个独立的弹性系数。若物体中的每一点出现有对称的方向，这些方向上的弹性特征相同，它就具有弹性对称。具有弹性对称的物体，广义虎可定律的方程和能量函数的表达式都可以简化，在弹性系数之间出现依赖关系。

Voigt对称性:可证明张量L对双指标ij和kl具有对称性。现在是36页\一共有59页\编辑于星期五2.4.1一个弹性对称面

弹性对称面就是指经过物体内的每一点都有这样的平面，在这个平面的对称方向上弹性特性是相同的。取xy坐标面与弹性对称面平行，z轴与弹性对称面垂直（如图），现研究体微元ABCDE的弹性对称问题。

由弹性对称面的定义知，当倒置z轴时，在坐标系(x、y、z)和（x、y、）中，体微元具有相同的应力、应变关系。换言之，弹性系数、不因倒置z轴而发生变化。

图2－5现在是37页\一共有59页\编辑于星期五弹性体单位体积的应变能是应变状态的单位函数，而且能量是标量，不因坐标的选择不同而改变其量值。但是当z轴变成轴时，有些物理量将变号。用u，v，w和u，v，w’分别表示两坐标中的位移分量，存在着下述关系：与有关的剪应变分别为：所以可得：

也就是说，z轴倒置时，与z方向有关的剪应变分量变号。

（2－16）（2－17）现在是38页\一共有59页\编辑于星期五由的表达式不难看出，除非含和的一次项的刚度系数等于零，否则不能保证的量值不变。于是，有刚度系数减少了8个，剩下13个。同样可以证明，柔度系数也剩下13个。于是，当z轴垂直弹性对称面时，应力、应变关系为：

（2－18）垂直于弹性对称面的方向为弹性主方向，坐标轴在弹性主方向时称为弹性主轴。上述的z轴即为弹性主轴。

现在是39页\一共有59页\编辑于星期五只有13个(21-8)弹性常数现在是40页\一共有59页\编辑于星期五如果其他应力分量为0,当沿弹性主轴拉伸时,除纵向伸长,横向收缩外,还会引起与主轴垂直的面(弹性对称面)内的剪应变,且弹性主轴方向不变现在是41页\一共有59页\编辑于星期五2.4.2三个弹性对称面——正交异性

若经过均质弹性体的每一点都有三个互相垂直的弹性对称面，则称之为正交异性弹性体。

取坐标x，y，z方向为弹性主方向。沿用前面的方法，将y轴转1800成y’轴，因为和变号，必须有：

新增加的等于零的刚度系数是后四个。再转动x轴不能增加新的等于零的刚度系数。将z轴转1800成z’轴同样可以得出，新增加的等于零的柔度系数为四个。独立的弹性系数剩下9个。从而得到应力应变关系为：现在是42页\一共有59页\编辑于星期五设仅有作用,其余应力分量为0,这时应变能对于上述两种坐标系计算时,变号，为了使W保持不变,必须使同理利用弹性主轴改变,弹性能不变的原理证明:现在是43页\一共有59页\编辑于星期五（2－19）从上可以看出，对于正交异性体，若坐标方向为弹性主方向，则正应力不引起剪应变，剪应力不引起线应变，反之亦然。

现在是44页\一共有59页\编辑于星期五没有拉压剪切耦合现象没有不同平面内的剪切耦合现象现在是45页\一共有59页\编辑于星期五2.4.3一个各向同性面——横观各向同性

在平面内一切方向的弹性特性均相同的平面称为各向同性面．如果过材料的每一点都有一个相互平行的各向同性面，就称为横观各向同性材料。

取oxy面为各向同性面，z轴垂直于该面．显然当涉及x或y方向时，刚度系数的下标可以不加区别，即，，而且可以证明。因此，独立的弹性常数减少为5个，应力应变关系简化为：

现在是46页\一共有59页\编辑于星期五（2－20）现在是47页\一共有59页\编辑于星期五2.4.4完全对称——各向同性

若经过均质体内每一点的任意方向上弹性特性均相同，即任意方向都是弹性主方向，则称之为各向同性体。显然，弹性系数之间存在下述关系：

对于各向同性材料，弹性系数与方向无关。应力应变关系为（显然只有两个独立的弹性常数了）：（2－21）现在是48页\一共有59页\编辑于星期五2.4.5小结

表2现在是49页\一共有59页\编辑于星期五2.5正交异性材料弹性常数的物理意义

在本构关系式中，刚度矩阵和柔度矩阵虽然是联系应力和应变的数学符号，但是也有其物理意义。通过简单的拉伸试验和剪切试验，可以导出各分量与材料机械性能之间的关系，称为用工程常数表示。

首先：研究沿弹性主方向1简单拉伸情况（如图2－6

）

（图2－6）现在是50页\一共有59页\编辑于星期五其应力和应变张量为：根据正交异性材料的本构关系，可得到：由此可得：同理，根据2和3方向的简单拉伸，可得：

现在是51页\一共有59页\编辑于星期五式中和即为人们所熟悉的工程弹性常数，是i方向的弹性模量，是i方向正应力引起j方向横向应变的泊松比。

其次：研究一个纯剪切试验，这时应力和应变张量为：同样根据本构关系可写出：于是得到：现在是52页\一共有59页\编辑于星期五类似的有：由此得到用工程常数表达的正交异性材料柔度矩阵为（刚度系数的表示可由刚度系数与柔度系数的关系得到）：现在是53页\一共有59页\编辑于星期五2.5.1正交异性体的弹性系数(由物理意义可求）1.刚度系数和柔度系数的关系

因为刚度矩阵和柔度矩阵是互逆的:，所以根据求逆法则，可得到相互关系：其中

在用刚度系数表示柔度系数时，只要将式中的C与