二、_弹性力学有限元法基本原理(一)程序.ppt

第二单元 弹性力学有限元法基本原理（一） 第一节???? 里兹法的有限元形式? 由于需要在整个求解区域上假设试探函数，经典里兹法在解决实际问题时，尤其是几何形状复杂的二、三维问题，具有局限性。 解决上述问题的办法是在求解区域上分片假设试探函数。 下面以一维直杆的分析为例子，研究基于里兹法的有限元位移法基本原理和求解过程。 离散化 考虑图2-1（a）所示的一维直杆结构，杆的总长度为3L。载荷与图1-4相同。杆分为三个子区域，称为单元。单元之间的连接点称为节点。这一步骤称为离散化。节点位移是问题的基本未知量。 分片假设单元上位移试探函数 为了使得上述假设位移场是“许可位移”，上述多项式待定系数必须满足一定约束关系（3个），显然，该问题的独立参量（广义坐标）只有3个。 把上述假设的分段线性位移场代入势能泛函： 分段计算上式积分，再应用驻值条件，可立即求出待定系数，位移场就完全确定，进而可求出各单元应力。 上述是原理上的里兹法有限元形式求解过程，其关键是采用分片多项式拟合全域上的可能位移场——分片试探函数。 这个做法正体现了有限元法的实质。 上面形式的分片位移试探函数有下列缺点： 必须对它进行调整，使其满足连续条件和边界约束条件； 多项式系数作为广义坐标缺乏明显的物理意义。 因此，上述不是通常意义上标准的有限元形式，仍然具有局限性，如对于二维以上的问题使各单元之间分片多项式保持连续性很难处理。 下面用节点位移未知量作为待定参数（广义坐标），得到其标准有限元形式。 重新构造单元内位移试探函数 显然，在使用插值试探函数情况下，整个杆上，由各单元位移函数拼接而成的试探位移场是连续的。只要我们记住 ，得到的就是可能位移场。这样的位移试探场以未知节点位移作为待定参量（广义坐标）。 下面进一步实施里兹法求解。 计算离散系统的势能泛函 载荷与上节例题相同： 为了适应矩阵形式运算，将单元势能矩阵表达式中的 用整体节点位移向量 代替，同时单元刚度矩阵扩展成总体规模（4×4）。则各单元相加后系统总势能为： 应用势能驻值条件： 单元应力由公式 得到。 第二节 常应变三角形单元解平面问题 上一节以受轴向力的弹性杆为例讨论了有限元位移法的基本原理和步骤，揭示了有限元法的本质特征。 本节讨论将该方法推广到解决弹性力学平面问题。弹性区域离散化采用3节点三角形单元（T3单元）。 1、结构离散化 连续求解区域划分为有限个三角形区域（单元）； 或者连续体离散为有限个小单元的结合体，单元之间在节点处连接。 问题的未知量转化为节点位移。 2、单元位移模式及插值函数 从离散结构中取出一个典型的分片区域（单元），研究二维区域上分片多项式假设位移场（单元位移模式）。 单元如图2-4所示。单元节点采用局部编号i，j，m（逆时针旋转），每个节点两个自由度（位移分量）： 在单元的三角形区域构造简单位移试探函数。 采用x，y的一次多项式： 将三个节点坐标分别代入上述位移多项式： 4）形函数的几何意义 根据三角形单元形函数的上述性质也可以推断，单元的假定位移在内部和边界上都线性分布，边界上的位移只跟两端节点位移有关，保证了整个求解区域上位移的连续性。 4、用节点位移表达单元应变和应力 所以，单元刚度矩阵第一列元素的物理意义是，第一个自由度位移为1，其它自由度位移为0时，要保持平衡须加在单元各节点上的节点力分量。显然，单元刚度矩阵元素表征单元刚性的大小，称为刚度系数。 由于这些节点力组成平衡力系，因此有： 一般而言，总刚度矩阵中非零元素很少，大多数元素为零，称为稀疏矩阵。有限元模型的节点数越多，非零元素所占比例越小。 有限元模型中属于同一单元的节点称为相关节点。 稀疏性的原因：从单元刚度矩阵叠加成总刚度矩阵的过程或从总刚度矩阵元素的物理意义看，只有相关节点之间才有非零刚度关联，而无论有限元模型规模大小，一个节点的相关节点数量总是有限的。因此，总刚度矩阵中大量的元素必然为零，稀疏性是必然的。 如果模型的节点编号有规则，非零元素将分布在以主对角线为中心的带状区域内。显然，对有限元方程的求解，这个性质非常重要。 衡量带状特性的指标是半带宽。 半带宽在有限元方程求解技术中是一个重要参数，它关系到求解中的存储量和计算时间。半带宽越小越好。 图2-6、图2-7分别为一个有限元结构和相应结构刚度矩阵非零元素分布。 总之，有限元结构刚度矩阵的上述性质是有限元法固有的特点，也是有限元法的优越性在计算技术上的体现，没有这些特性，就不会有有限元如今的发展和