基于有限元法的振动郑佳文选读.ppt

* 式中 I是二阶单位矩阵；Ni 、Nj 、Nm 是坐标的函数，它们反映了单元的位移状态，所以一般称之为形状函数，简称形函数。矩阵 [N] 叫做形函数矩阵。三节点三角形单元的形函数是坐标的线性函数。单元中任一条直线发生位移后仍为一条直线，即只要两单元在公共节点处保持位移相等。则公共边线变形后仍为密合。 (5-11) 也可写成矩阵形式 * 三、应 变 有了单元的位移模式，就可以利用平面问题的几何方程 求得应变分量。将 (e) 、(f) 两式代入上式，即得： (g) * 可简写成 ?其中 [B] 矩阵叫做单元应变矩阵，可写成分块形式 而子矩阵 由于?和bi 、bj 、bm 、ci 、cj 、cm 等都是常量，所以矩阵[B]中的诸元素都是常量，因而单元中各点的应变分量也都是常量，通常称这种单元为常应变单元。 （i , j , m轮换） (3-15) (3-14) (3-13) * 四、应 力 求得应变之后，再将（3-13）式代入物理方程 ， 便可推导出以节点位移表示的应力。即 (5-16) (h) (5-17) 令 则 * 其中 [S]叫做应力矩阵，若写成分块形式，有 对于平面应力问题，弹性矩阵[D]为 (5-18) (i) 所以，[S]的子矩阵可记为 （i , j , m轮换） (5-19) * 注意到（5-7）式，则有 (5-21) 由（5-19）、（5-20）式不难看出，[S]中的诸元素都 是常量，所以每个单元中的应力分量也是常量。 可见，对于常应变单元，由于所选取的位移模式是线 性的，因而其相邻单元将具有不同的应力和应变，即在单 元的公共边界上应力和应变的值将会有突变，但位移却是 连续的。 * 单元的变形能 单元的刚度矩阵，为一6X6的矩阵。 单元的动能 单元的质量矩阵，为一6X6的矩阵。 * 矩形单元 v1 y u1 1 u2 v2 2 u1 v1 4 u1 v1 3 x 矩形单元的位移假定采用双线性插值多项式 * 等参单元 四边形四节点等参数单元 四边形四节点单元位移模式: 其中 母单元 * 若以母单元上12边为例，通过映射可得在平面内任一直线，12边的方程为 η= -1,代入 通过局部座标与整体座标的映射关系把母单元变换到整体座标上成为一个任意四边形用于离散结构物，它能适合于任意曲边的形状。 * 同理可得 把以参数 ξ 代表的 x方程和 y方程消去 ξ ，则得 x , y 所组成的直线方程 y=kx+b 所以母单元上的四个边都可以通过映射在x, y座标面上得出一个任意四边形，用该四边形离散结构物。 * 四边形四节点 单元的应变 i = 1,2,3,4 其中 * 其中 由于形函数 N 是 x, y的函数，现对ξ, η求导，这是复合函数求导 把上二式写成矩阵形式 * 其中 称为雅可比(Jacobian)矩阵 为了把 [Bi]矩阵中的Ni,x和 Ni,y化成 Ni,ξ和 Ni,η ，则代入下式 其中雅可比矩阵的逆阵由下式给出 也可把形函数Ni对的ξ, η 偏导数写成通式: i= 1,2,3,4 * 四边形四节点 单元的应力 对于平面应力情况 * 四边形四节点单元的刚度矩阵[K]e: [K]e可划分成四行四列的子矩阵 : i, j = 1,2,3,4 * 对于平面应力情况: 等效节点力计算: 1. 集中力的等效节点力:由于计算复杂，把有集中力处设置为节点。 2. 体枳力的等效节点力: 两类平面问题： （1）平面应力问题 对于厚度远小于其他尺寸的薄板，且弹性体仅受平行于板面沿厚度方向不变化的外力作用时，可按平面应力问题来处理 * 平面问题有限元法 （2）平面应变问题 对于长度大于其横向尺寸的弹性体，且仅平行于截面、沿其长度方向不变化的作用力时，可按平面应变问题来处理，如下图 * 弹性矩阵 将弹性构件划分为单元时，要确定单元的种类、数目和布局。 单元种类可根据需要选择。选择精度比较高的单元，单元数目就可以少一些。单元尺度小一些，但单元的数目就要增加，计算精度自然会提高，但这涉及计算的时间和成本，有时还要考虑计算机容量的限制。 在应力集中的部位以及应力比较剧烈的区域，网路的划分宜稠密一些；反之，在应力较小的区域，单元的尺度可大一些，但疏密的过渡应较为平缓。一个单元的各边长度不宜过为悬殊。 构件上载荷的作用点除应取为节点。各个单