第2讲有限元法基本理论讲解.ppt

Step 2.单元特征分析—位移函数 单元节点位移列阵： 单元等效节点力列阵： 单元节点坐标列阵： xi xj 利用节点条件： 构造单元位移函数： 单元特征分析—应变、应力 形状函数矩阵 应变矩阵 应力矩阵 由几何方程 由物理方程 单元特征分析—应变能、外力势能 单元应变能计算： Ke—单元刚度矩阵 单元外力功 要求单元的外力功，关键是求出单元上的体积力和面积力等效作用到节点上的力的大小。 单元等效节点力 Step 3. 单元集成——应变能 扩充 叠加 K—整体刚度矩阵 整体集成 Step 3. 单元集成——外力功 扩充 叠加 Step 3. 单元集成——系统总势能 Step 4. 变分处理 Step 5: 处理位移边界条件并求解 Step 6: 计算每个单元的应变及应力 换个思路：多元函数求极值 总结：有限元求解基本步骤 几何离散：m个单元和n个节点的组合体 单元特征分析：单元应变能，单元外力势能（等效节点载荷） 单元集成：系统的总势能 变分处理：系统的平衡方程（组） 应用位移边界条件求出节点位移 由节点位移求出单元的应变、应力 整体节点位移列阵 整体等效节点力列阵 几何离散 单元的节点描述 单元分析 单元位移场模式， 所有物理量的表达（均用节点位移表示） 整体平衡方程 扩充 叠加 单元集成 扩充 叠加 处理边界条件，求得满足位移边界条件的位移场 求解其他力学量，如应变和应力等 有限元解的误差及收敛性 有限元解的误差 收敛性相关准则 位移单元解的下限性 有限元解的误差 模型误差：离散误差、边界条件误差 计算误差：舍入误差、截断误差 对于有限元这种数值计算方法，一般总是希望随着网格的逐步细分所得到的解能够收敛于问题的精确解。根据前面的分析，可知在有限元分析中，一旦确定了单元的形状之后，位移模式的选择将是非常关键的。由于载荷的移置、应力矩阵和刚度矩阵的建立等等，都依赖于单元的位移模式，所以，如果所选择的位移模式与真实的位移分布有很大的差别，那么就很难获得良好的数值解。为了保证解答的收敛性，要求位移模式必须满足以下三个条件，即 1）位移模式必须包含单元的刚体位移。 2）位移模式必须包含单元的常应变。 3）位移模式在单元内要连续、且在相邻单元之间的位移必须协调。 1）位移模式必须包含单元的刚体位移。也就是说，当节点位移是由某个刚体位移所引起时，弹性体内将不会产生应变。所以，位移模式不但要具有描述单元本身形变的能力，而且还要具有描述由于其它单元形变而通过节点位移引起单元刚体位移的能力。 2）位移模式必须包含单元的常应变。每个单元的应变一般都是包含着两个部分：一部分是与该单元中各点的坐标位置有关的应变（即所谓各点的变应变）；另一部分是与位置坐标无关的应变（即所谓的常应变）。从物理意义上看，当单元尺寸无限缩小时，每个单元中的应变应该趋于常量。因此，在位移模式中必须包含有这些常应变，否则就不可能使数值解收敛于正确解。 3）位移模式在单元内要连续、且在相邻单元之间的位移必须协调。当选择多项式来构成位移模式时，单元内的连续性要求总是得到满足的，单元间的位移协调性，就是要求单元之间既不会出现开裂也不会出现重叠的现象。通常，当单元交界面上的位移取决于该交界面上节点的位移时，就可以保证位移的协调性。 在有限单元法中，把能够满足条件1和2的单元，称为完备单元，完备单元是收敛的必要条件。满足条件3的单元，叫做协调单元或保续单元。如三节点三角形单元，均能同时满足上述三个条件，因此属于完备的协调单元，完备的协调单元是收敛的充分条件。在某些梁、板及壳体分析中，要使单元满足条件3比较困难，所以实践中有时也出现一些只满足条件1和2的单元，其收敛性往往也能够令人满意。目前，完备而不协调的单元，已获得了很多成功应用。 在有限元分析中，将实际连续体分成许多单元体的组合后，根据线性或非线性位移的假定，人为地选择一个位移场，通过这些措施所得到的模型比实际连续体的刚性要高。 因而，近似模型的刚度是实际连续体刚度的上界。若选择不协调单元，那么这种模型可能由于单元分离、叠加或单元之间形成铰而降低刚性，这种影响就有可能使得不协调元比应用协调元所得的结果要好。不过，应用不协调单元事先不能肯定所得的刚度是真实刚度的上界。换句话说，不协调元不一定象协调单元那样刚硬，可能比较柔软，因此有可能会比协调单元收敛得快。 经验证明，根据巴斯卡（Pascal）三角形来选择二维多项式的各项。在二维多项式中，如果包含有对称轴一边的某一项，那么就必须同时包含有另一边的对称项。选择多项式位移模式时，还应该要考虑的一个因素是，多项式中的项数必须等于或稍大于单元边界上的外节点的自由度