有限元 力学基础.ppt

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * ５.３.5单元刚度矩阵的物理意义及其性质 单元刚度矩阵的物理意义： 将 写成分块矩阵 写成普通方程 其中 表示结点S(S=i,j,m)产生单位位移时，在结点r(r=i,j,m)上所需要施加的结点力的大小。 ５.３.5单元刚度矩阵的物理意义及其性质 单元刚度矩阵的物理意义： 将结点力列矩阵 与结点位移列矩阵 均展开成(6*1)阶列矩阵，单元刚度矩阵相应地展开成(6*6)阶方阵： 元素K的脚码，标有“-”的表示水平方向，没有标“-”的表示垂直方向。 ５.３.5单元刚度矩阵的物理意义及其性质 单元刚度矩阵的物理意义： 单元刚度矩阵的每一个元素都有明显的物理意义。 表示结点S(S=i,j,m)在水平方向、垂直方向产生单位位移时，在结点r(r=i,j,m)上分别所要施加的水平结点力和垂直结点力的大小。例如 表示结点j在垂直方向产生单位位移时，在结点i所需要施加的水平结点力的大小。 ５.３.5单元刚度矩阵的物理意义及其性质 单元刚度矩阵的性质： 1)对称性： 是对称矩阵 2)奇异性： 是奇异矩阵， 单元刚度矩阵所有奇数行的对应元素之和为零，所有偶数行的对应元素之和也为零。由此可见，单元刚度矩阵各列元素的总和为零。由对称性可知，各行元素的总和也为零。 ５.３.5单元刚度矩阵的物理意义及其性质 单元刚度矩阵的性质： 例题：求下图所示单元的刚度矩阵，设 1、求[B] 2、求[D] 3、求[S] 4、求 ５.３.6整体分析 将各单元组合成结构，进行整体分析。 整体分析分4个步骤 1、建立整体刚度矩阵； 2、根据支承条件修改整体刚度矩阵； 3、解方程组，求出结点位移；(消去法与叠加法) 4、根据结点位移求出应力。 作业： 5-2；5-5。 下次上机实验并交实验报告电子档。 以班为单位，班级，学好，名字等。 * * * * * * * * * * * * * * * * * * * * * * * * * 5.3平面问题的有限单元法 5.3-1 单元位移函数 三结点三角形单元 六个节点位移只能确定六个多项式的系数，所以平面问题的3结点三角形单元的位移函数如下， 所选用的这个位移函数，将单元内部任一点的位移定为座标的线性函数，位移模式很简单。 位移函数写成矩阵形式为： 5.3-1 单元位移函数 最终确定六个待定系数 5.3-1 单元位移函数 如果弹性体的位移分量是座标的已知函数，则可用几何方程求应变分量，再从物理方程求应力分量。但对一个连续体，内部各点的位移变化情况很难用一个简单函数来描绘。 有限单元法的基本原理是分块近似，即将弹性体划分成若干细小网格，在每一个单元范围内，内部各点的位移变化情况可近似地用简单函数来描绘。对每个单元，可以假定一个简单函数，用它近似表示该单元的位移。这个函数称为位移函数，或称为位移模式、位移模型、位移场。 对于平面问题，单元位移函数可以用多项式表示， 多项式中包含的项数越多，就越接近实际的位移分布，越精确。但选取多少项数，要受单元型式的限制。 5.3.1 单元位移函数 令 （下标i，j，m轮换） 简写为 [I]是单位矩阵， [N]称为形态矩阵， Ni称为位移的形态函数 5.3-1 单元位移函数 选择单元位移函数时，应当保证有限元法解答的收敛性，即当网格逐渐加密时，有限元法的解答应当收敛于问题的正确解答。因此，选用的位移模式应当满足下列两方面的条件： (1) 必须能反映单元的刚体位移和常量应变。 6个参数 到 反映了三个刚体位移和三个常量应变。 (2) 必须保证相邻单元在公共边界处的位移连续性。 (线性函数的特性) 5.3-1 单元位移函数 例题：图示等腰三角形单元，求其形态矩阵[N]。 5.3-1 单元位移函数 由三角形的面积 5.3-2 单元载荷移置 连续弹性体离散为单元组合体时，为简化受力情况，需把弹性体承受的任意分布的载荷都向结点移置(分解)，而成为结点载荷。如果弹性体