【力学教案】 有限元法基础1

《【力学教案】 有限元法基础1》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第六章有限元法概述第一节单元分析简例1、单元分析的主要任务：求出单元节点位移和节点力之间的转换关系。在推导此关系时规定：力和位移的方向若和坐标轴正方向一致者为正。先举一个简单例子，图1示一拉压弹簧，弹簧系数为常量c，其轴线和x坐标轴重合，令此弹簧为一个单元，则弹簧的两端点i,j是此单元的两个节点。设在节点i,j上分别有轴向力和轴向位移。则当节点对单元有的作用力时，单元对节点有大小相等、方向相反的反作用力，节点力：这节点和单元之间的作用力和反作用力都称为节点力，对单元来讲节点力是作用于单元之力。2、节点力和节点位移的关

2、系。图1可分解为两步1）设节点j被固定，节点i产生正位移，则此时节点i作用在单元上的力是iUiuiyjujUjx而节点j作用在单元上的力是2）是设节点i被固定，节点j产生正位移，此节点j对单元的作用力是节点i对单元的作用力是将两式合并，就得到由式可以看出一个节点上的节点力不仅决定于本节点的位移，而且也决定于本单元其他节点的位移。设以表示单元节点力向量：以表示节点位移向量：则式（1.1）可改写成：式中式中（1.2）就是单元节点位移和节点力之间的转换关系。是转换矩阵，称为单元刚度矩阵。所以单元分析的任务就是要求出本单元的

3、刚度矩阵。。§2单元位移模式对于一个复杂的弹性体，要想用某种函数来描述整体内任一点的位移是不大可能的。但当把弹性体离散化为许多细小的单元，则在一个单元的局部范围内是可以把某一点的位移近似的表达为其坐标的函数，这表达式称为单元位移模式。在单元范围内建立位移模式是有限元法的特色。一般，位移模式都是用待定系数法求得。为了说明之，还以上述图1的拉压弹簧为例，设式中—单元中任一点的位移x—单元中任一点的坐标待定系数可由已知的节点坐标及位移来确定，这是因为既然位移模式代表了单元内任一点的位移，那么也应包括节点的位移，所以每代入一

4、个节点的坐标及位移，就有一个解待定系数的方程，单元上有几个节点位移，或者说有几个自由度，就可以确定几个系数。上述弹簧单元只有两个自由度，故只能列出具有两个待定系数的多项式令时，时，式中—节点的坐标—节点的位移解方程（1.3），得将式代入得到(1.6)式（1.6）就是弹簧单元的位移模式：它表明单元中任一点的位移都是其坐标的函数。同时也表明只要已知单元节点位移，就可以通过位移模式推算出单元内任一点的位移。平面问题的单元形式有很多种，其中最基本的形式是三节点三角形单元，以下将以这种单元为典型，讨论平面问题的单元分析及整体分

5、析。3、三角形单元的位移模式xmyisjuivivuviui(x,y)umvmO图2示一个三角形单元，节点的坐标分别是其沿坐标轴方向的位移分别为显然，此单元共有六个自由度，故其位移模式可写成具有六个待定系数的多项式：（1、7）式中—单元中任一点的坐标—该点沿轴的位移1）先讨论水平X方向的位移：将三个节点的坐标及其水平位移代入式（1.7）中的第一式，得：解方程，可得式（1.9）中三组和都是只和节点坐标有关的常数i,j,m（1．10）记号表示格式（1.10）中的进行顺序轮换后，可以得出另外两组的表达式。为了使三角形面积不

6、为负值，应是逆时钟方向排列。将式（1.9）代入（1.7）的第一式，得令　　　　　（１．１３）同样可证明　　　（1.14）　式（1.13），（1.14）就是三角形单元的位移模式，将式（1.13），（1.14）合在一起得：或简写成这就是由单元节点位移求得单元内部各点位移的转换式，是转换矩阵。由（1.13），（1.14）可知，当或而其他节点位移为零时，单元内任一点的位移为这就是说，节点i发生单位位移时，函数表示单元内部的位移分布形状，故也称为位移的形状函数，简称为形函数，矩阵称为形函数矩阵。通过形状函数，就可以由单元体各接

7、点位移计算出单元体内各点位移由式（1.12）可知，三节点三角形单元的形函数是坐标的线性函数，从而其位移模式（1.13），（1.14）也是坐标的线性函数。线性位移模式的一个特点：单元中的任一条直线在位移后，仍然是一条直线。所以相邻的三节点三角形单元的公共边线，在变形后仍是一条直线，只要两单元在公共节点处保持位移相等，则公共边线在变形后仍然保持密合。第二节单元应变在上节中研究了由单元节点位移求单元内各点位移的方法，本节将研究如何由各点位移求各点相应的应变，并进一步导出由单元节点位移求单元内各点应变的方程。（17）为几何方

8、程。将式（1.13），（1.14）deu,v分别代入式（1.17）得：或写成：简写成　　　（１．１９）这就是由节点位移求单元内各点应变的转换式，其转换矩阵为称为几何矩阵，他的每项元素是由单元三个节点的坐标按式（1.10）直接算出。由于矩阵的元素是常数，则由（1.19）式可知，单元内各点的应变也都是常数，这是采用了线性位移模式的必然结果。所以三节