第9章平面单元9.pptVIP

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