常应变三角形单元的应变矩阵.ppt

有限元分析及应用 Finite Element Analysis and Application 有限元分析及应用 Finite Element Analysis and Application 第三章 平面问题的有限单元法 3-1、有限单元法的计算步骤 3-2、平面问题的常应变(三角形)单元 3-3、单元刚度矩阵 3-4、单元刚度矩阵的物理意义及其性质 3-5、平面问题的矩形单元 3-6、六节点三角形单元 3-7、单元载荷移置 3-8、整体分析 3-9、整体刚度矩阵的形成 3-10、支承条件的处理 3-11、整体刚度矩阵的特点 3-1 有限单元法的计算步骤 弹性力学平面问题的有限单元法包括五个主要步骤： 1、所分析问题的数学建模 2、离散化 3、单元分析 4、整体分析与求解 5、结果分析 3-2 平面问题的常应变(三角形)单元 有限单元法的基础是用所谓有限个单元的集合体来代替原来的连续体，因而必须将连续体简化为由有限个单元组成的离散体。对于平面问题，最简单，因而最常用的单元是三角形单元。因平面问题的变形主要为平面变形，故平面上所有的节点都可视为平面铰，即每个节点有两个自由度。单元与单元在节点处用铰相连，作用在连续体荷载也移置到节点上，成为节点荷载。如节点位移或其某一分量可以不计之处，就在该节点上安置一个铰支座或相应的连杆支座。如图3-1 3-2 平面问题的常应变(三角形)单元 3-2 平面问题的常应变(三角形)单元 3-2 平面问题的常应变(三角形)单元 3-2 平面问题的常应变(三角形)单元 3-2 平面问题的常应变(三角形)单元 据弹性力学几何方程得 单元的应变分量 由于三节点三角形单元的位移函数为线性函数，则单元的应变分量均为常量，故这类三角形单元称为常应变单元（位移在单元内和边界上为线性变化，应变为常量） 3-2 平面问题的常应变(三角形)单元 2、形函数的特点及性质 1)形函数Ni为x、y坐标的函数，与位移函数有相同的阶次。 2)形函数Ni在i节点处的值等于1，而在其他节点上的值为0。 即 3-2 平面问题的常应变(三角形)单元 3-2 平面问题的常应变(三角形)单元 (3)位移函数在单元内部必须连续位移。因为线性函数，内部连续 (4)位移函数必须保证相邻单元在公共边界处的位移协调（即在公共边界上位移值相同）。如右图 设公共边界直线方程为y=Ax+B，代入位移函数可得：边界上位移为 显然，u,v仍为线性函数，即公共边界上位移连续协调。 综上所述，常应变三角形单元的位移函数满足解的收敛性条件，称此单元为协调单元 3-2 平面问题的常应变(三角形)单元 3-2 平面问题的常应变(三角形)单元 3-2 平面问题的常应变(三角形)单元 4、应力、应变矩阵 将位移函数代入平面问题几何方程，得应变矩阵： 3-2 平面问题的常应变(三角形)单元 应力矩阵 由平面问题物理方程得： 应变矩阵[B]反映了单元内任一点的应变与节点位移间的关系 应力矩阵[S]反映了单元内任一点的应力与节点位移间的关系 显然，常应变三角形单元的应变矩阵[B]为常量矩阵，说明在该单元上的应力和应变为常值。由此可见，在相邻单元的边界处，应变及应力不连续，有突变。 3-3 单元刚度矩阵 3-3 单元刚度矩阵 3-3 单元刚度矩阵 3-3 单元刚度矩阵 3-3 单元刚度矩阵 3-3 单元刚度矩阵 3-3 单元刚度矩阵 3-3 单元刚度矩阵 3-3 常应变三角形单元的刚度矩阵 单元刚度矩阵 可记为分块矩阵形式 将应变矩阵[B]的分块阵代入单元刚度矩阵，可得其子块计算式： 对于常应变三角形单元，考虑平面应力问题弹性矩阵[D]，可得 3-4 单元刚度矩阵的物理意义及其性质 3-4 单元刚度矩阵的物理意义及其性质 3-4 单元刚度矩阵的物理意义及其性质 3-4 单元刚度矩阵的物理意义及其性质 3-4 单元刚度矩阵的物理意义及其性质 3-4 单元刚度矩阵的物理意义及其性质 1）单元刚度矩阵是对称阵，(只要证明 ) 2）单元刚阵主对角线元素恒为正值；因为主对角元素 表示力的方向和位移方向一致，故功总为正值。 3）单元刚阵是奇异阵，即|K|=0，这是因为计算单元刚阵时没有对单元的节点加以约束，虽然，单元处于平衡状态，但容许单元产生刚体位移，故从单元刚度平衡方程不可能得到唯一位移解 ，只能得到唯一的节点力解。 4）单元刚阵所有奇数行的对应元素之和为零，所有偶数行的对应元素之和也为零。由此可见，单元刚阵各列元素的总和为零。由对称性可知，各行元