第十九章SECTION5弹性理论和有限元解法.docVIP

§5 弹性理论与有限元解法 有限元法用来解弹性体的小变形问题是很有效的.这是由于变形能与外力势能或总势能可以表示为形式划一的二次泛函.由于弹性力学问题提法,特别是边界条件的复杂性,要使有限元解法灵活运用,还需要对各种问题作些分析与讨论.因此本节除了介绍各类问题与§1中B、D等有关的基本关系式外，还着重分析边界条件在有关变分问题中所起的作用，并讨论各类变分问题与微分方程定解问题的等价性。 三维的弹性问题 本节在直角坐标系(x,y,z)中讨论弹性体内部受力与形变的情况.为了简化论述,假设弹性体是单连通、均匀并各向同性的,必要时还用通常的矢量记号以突出其力学意义. [应力与平衡方程] 弹性体在外力、温差等作用下，各部分之间将产生内力，表示内力的大小是作用面上受力的强度或单位面积所受的力,即所谓应力.对弹性体内的一点P,可在其邻近作一微元六面体，其棱边平行于坐标轴。六面体的六个面中有三个面其外法线方向分别与x轴、y轴或z轴同向（其余三个则与坐标轴反向）,可分别称为x,y,z坐标面(或第1,2,3坐标面).设分别表示作用于x,y,z坐标面的应力,并记的三个分量为(对其余三个面的应力取号相反.例如,图19.18上所示的分量全是正的),这九个分量(i,j=1,2,3) 构成一个张量,称为应力张量. 从图可看出,分量表示第i坐标面上的正应力(受拉取正,受压为负);而分量则表示沿第i坐标面剪应力的二分量(使扭角变成锐角的为正).由微元六面体力矩的平衡可得剪应力互等定律,即 = () 因此应力张量是对称的,其分量只有六个是独立的.在有限元解法中,一般是把这六个分量按如下次序排成列矢量,并记作 设过点P作一任意的斜面,其法线n的方向余弦为(),则利用与三坐标面围成的四面体的平衡条件可得作用于该斜面的应力的三个分量 (j=1,2,3) 于是在方向余弦为()的直线上的投影,即其分量为 这表明,一个点的应力状态完全由该点的六个应力分量{?}所决定.例如,作用于上述斜面的正应力为 而沿该斜面的剪应力方向取,大小等于 反过来,对于的斜面称为P点的应力主面,相应的法线称为P点的应力主轴,而其正应力称为P点的主应力.可以证明,在弹性体内任意一点,一定存在三个互相正交的主应力,而且其中最大(小)的一个就是该点的极大(小)正应力.三个正应力之和 称为体积应力,它在坐标变换下是个不变量,因而等于三个主应力之和. 设,表示作用于P点的体力(单位体积的外力),对体积元素作积分或均值分析都可以推出力的平衡方程 (i=1,2,3) (26) [应变与几何方程] 弹性体内任一点P(x,y,z)在小变形后移到其位移函数为 式中 它们是(x,y,z)的微量函数,假定有一微小线段PN=dr其方向余弦为(),经过小变形变为,则沿该方向的正应变定义为单位伸长,即 从变形前后的与dr关于位移的表达式不难得出 展开右端的根式并略去高阶无穷小量(即位移导数的高次项),就得到 设另一线段PM的方向余弦为(),变形前的夹角,则 (27) 设变形后的夹角(图19.19),则 根据变形前后两线段大小、方向的变化,不难得出 (28) 对照(27),(28)可知,只要在P点给定如下六个导数值: (29) 就可以完全确定P点邻近的变形状态.表示沿第i坐标轴的正应变,表示经过小变形第i,j坐标方向之间的直角改变量即所谓剪应变(如图19.20所示,变形后成锐角为正,成钝角为负).这六个量称为应变分量,记作 同样可以证明,在弹性体内任意一点,一定存在互相正交的应变主轴,变形后三轴交角仍然保持直角,即剪应变为零;三主轴的正应变称为主应变,而且其中最大(小)的一个就是该点的极大(小)的正应变.三个正应变之和 称为体积应变,也是个不变量,而且表示微元中每单位体积的改变量.对于各向同性体来说,应力主轴与应变主轴的方向还是一致的.§1中所述的变形能密度 正是把坐标系变换到共同的主轴方向并应用虎克定律而推导出来的. 关系式(29)称为几何方程,其矩阵形式为 [物理方程与弹性系数] 对杆件作简单拉、压的小变形实验证明,单独的轴向(取为x轴)力不会引起剪应变,其正压力与正应变之间具有如下的线性关系 或 如果考虑到沿x方向的伸长还伴有侧向收缩,则产生沿y,z方向的正应变 或 这就是虎克定律,式中系数E,v分别称为弹性模数和泊松比.从胡克定律就可以推导出应变与应力之间的一般关系式,即所谓物理方程 式中称为剪切弹性模数. 从(30)求逆得出应力与应变之间的关系式 或写成