弹性力学第6用有限单元法求平面问题.pptVIP

§6-1 基本量和基本方程的矩阵表示 §6-2 有限单元法的概念 FEM的概念，可以简述为： 结构力学研究的对象是离散化结构。如桁架，各单元(杆件)之间除结点铰结外，没有其他联系(图(a))。 图(c)与图(a)相比,两者都是离散化结构；区别是，桁架的单元是杆件，而图(c)的单元是三角形块体（注意:三角形单元内部仍是连续体）。 (2)应用几何方程,由单元的位移函数d,求出单元的应变 (2)应用几何方程,由单元的位移函数d,求出单元的应变 1.有限单元法求解问题的基本步骤是什么？ §6-3 单元的位移模式与解答的收敛性 在结点三角形单元中,可以假定位移分量只是坐标的线性函数,也就是假定： 三结点三角形单元的位移模式，略去了2次以上的项，因而其误差量级是 且其中只包含了x,y的1次项，所以在单元中Ni的分布如图(a)所示,u,v的分布如图(b)、(c)所示。 对式(a)求应变,得： (1)和(2)是必要条件，而加上(3)就为充分条件。 1.应用泰勒级数公式来选取位移模式，为什么必须从低次项开始选取？ 2.试考虑：将结构力学解法引入到求解连续体的问题时，位移模式的建立是一个关键性工作，它使得单元(连续体)内部的分析工作都有可能进行了。 §6-4 单元的应变列阵和应力列阵 应用几何方程，求出单元的应变列阵： 应用几何方程，求出单元的应变列阵： 如果在位移模式中取到泰勒级数中的二次幂项，略去 高阶小量，试考虑位移、应变和应力的误差量级。 §6-5 单元的结点力列阵与劲度矩阵 假设发生一组结点虚位移　　 则单元内任一点(x,y)的虚位移为 单元内任一点(x,y)的虚应变为 代入虚功方程：在单元中，外力(结点力 )在虚位移(结点虚位移　 　)上的虚功，等于应力 在虚应变 上的虚功，即： §6-6 荷载向结点移置,单元的结点荷载列阵 1、等效原则 原荷载 作用于单元中任一点(x,y),为单位厚度上的作用力; 移置荷载 作用于结点ijm。 1. 试导出书中例题(P119)的荷载移置公式。 §6-7 结构的整体分析 结点平衡方程组 在单元分析中，从单元的结点位移→求位移分布→求应变→求应力→求结点力，为单元的内力分析；外荷载移置到结点荷载，为单元的外力分析。 代入式(a)，可表示为 列出图示结构i 结点的平衡条件。 §6-8 解题的具体步骤 单元的划分 3、使用已编好的程序进行上机计算。事先须将有限单元法的公式，计算方法和步骤都编入程序。 2、单元划分注意事项 §6-9 计算成果的整理 应力的波动性在三结点三角形单元中较为显著。这是由于计算出的应力的精度较低。假设Ⅰ单元的应力成果为　　 ，其中 为真解， 为误差。则由于在结点都列出了平衡方程并令其满足，从而使相邻的Ⅱ单元的应力趋近于 。这就产生了应力的波动性。 在面力边界线附近，求得的应力误差较大。可采用向外插值的方法（例抛物线插值）来解决。 §6-10 计算实例 §6-11 应用变分原理导出有限单元法基本方程 例题 某等腰直角三角形单元ijm如图所示,已知在所选取的坐标系中,单元结点坐标分别为: 应用 可得 应用公式 可得该单元的应力转换矩阵为 可得该单元的应力转换矩阵为 应用教材式(6-37)及式(6-38)可得该单元的单元刚度矩阵为 现考察结点力与单元中的应力之间的关系。为了简单起见，假定只有结点i发生位移ui，如右图(a)所示。由上面的单元刚度矩阵得相应的结点力为： 其中， 。相应的结点位移及结点力如图所示 另一方面，由于发生了位移ui,则根据上面得到的应力转换矩阵，可得应力分量为 如上图(b)中单元的两直角面所示。根据单元的平衡条件,还可得出斜面上的应力,如图所示。若将这三个面上的应力分别按静力等效原则积铢累寸到结点上去，可以得到图(a)中相同的结点力。 ■ 思考题 试求出书中例题(P117)的位移模式。 在FEM中,须将作用于单元上的外荷载向结点移置，化为等效结点荷载， (2)变形体静力等效原则: 即在任意的虚位移上，使原荷载与移置荷载的虚功相等。 刚体静力等效原则只从运动效应来考虑,得出移置荷载不是唯一的解; 变形体的静力