小变形问题的基本方法.docVIP

第三章 小变形问题的基本方法 3.1 小变形问题的几何方程、平衡微分方程 3.1.1 应变和几何方程 变形体内任一点的位移为 f=[u v w]T 该点的应变为 这是工程应变的表达式，它也可以用一个对称的二阶应变张量表示 工程应变和应变张量之间的关系是 (3-1) 根据弹性力学理论，在小变形情况下，位移分量和工程应变分量之间的关系可表示为 (3-2) 式中，L是微分算子 (3-3) 式（3-2）是变形体内各点的几何关系。在它的边界上也要满足几何条件，这就是在位移边界Su上，应当满足位移边界条件 (3-4) 3.1.2 应力和平衡方程 变形体内任一点的应力为 这是工程上常用的表示，也可以用一个对称的二阶应力张量表示 两者的关系是 (3-5) 变形体内部的应力受平衡条件的制约，设作用在变形体上的体力为 (3-6) 作用在变形体面力边界上的面力为 (3-7) 那么物体内部任一点的平衡方程为 (3-8) 在面力边界上，应满足应力边界条件 (3-9) 其中n是方向余弦矩阵 (3-10) 根据虚功原理，可以得出与式（3-8）和（3-9）等价的虚功方程 (3-11) 式中是虚位移矢量，在位移边界上应满足=0，是与虚位移相应的虚应变矢量 (3-12) 虚功方程是一个以积分形式表示的平衡方程，不仅对连续可导的函数成立，对广义函数也是成立的，因而由虚功方程导出有限元的支配方程，可降低对位移函数（或其它场函数）的光滑性要求。例如，在求解以位移为基本未知量的微分方程时，要求位移函数是二阶连续可导的。如果从虚功方程出发只要求位移连续和分片可导（即分片光滑）。这是有限元方法的优点。 3.1.3 弹性材料的本构关系 前面讨论的几何方程和平衡方程对小变形问题都是正确的，而不管材料是弹性的还是非弹性的。现在开始研究材料的力学性质，也就是应力与应变之间的关系，即材料的本构关 系或本构方程。本构关系一般是从实验和经验中观察到的特性出发，在某些理论的假设下得出的，它只是反映了材料在某些方面的性态，因而是描述真实材料的一种理想化的模型。 物体在外力作用下产生变形和应力。在应力水平不高的情况下，当移去外力后，变形和应力也随之消失，物体的这种性质称为弹性。对理想的弹性材料，其内部各点的应变和应力存在一一对应的关系。变形与加载过程无关，只决定于当前的应力，这就假设了变形和应力都是瞬时产生的，没有时间上的先后。 最简单的弹性模型是线弹性模型，此时应力和应变之间的关系是线性关系 (3-13) 弹性矩阵D是6×6的对称矩阵。对于一般的线弹性材料，D中有21个独立的弹性常数。对正交各向异性材料，有9个独立常数。横观各向同性材料有5个。对各向同性材料。D有2个独立常数，杨氏模量E和泊松比μ。D的表达式为 (3-14) 式中、G称为拉密系数，G又称剪切模量 (3-15) (3-16) D矩阵还可以表示为 (3-17) 其中K称为弹性的体积模量 (3-18) 更一般的弹性模型是非线性弹性模型，此时应力和应变之间存在非线性关系 (3-19) 其中是6×6的材料系数矩阵，其中的元素是或的函数。对于各向同性的非线性弹性体，与（3-14）或（3-17）中D的形式相同，但其中的E和μ是应变或应力的函数，现分别记为和。对于单向应力状态（图3.1），和分别是割线杨氏模量和割线泊松比，因而称为非线性弹性的割线弹性矩阵。如果采用增量形式的本构关系，式（3-19）成为 (3-20) 为切线弹性矩阵，其中的元素E和μ分别用切线杨氏模量和切线泊松比表示，它们也是应变状态或应力状态的函数，与历史无关。 以上直接按定义建立应力和应变之间的函数关系，这一方法称Cauchy方法。还有一种Green方法是根据应变能