弹性力学--用有限单元法解平面问题.ppt

一、组集整体劲度矩阵[ K ] 整体劲度矩阵[ K ]以整体编码排列，而单元劲度矩阵[ k ]e以局部编码排列，故要把单元劲度矩阵的各子矩阵按对应的整体编码集成。 单元号 ① ② 局部编码 整体编码 i 4 1 j 1 4 m 3 2 组集整体劲度矩阵[ K ] 的方法为： ① 将[ K ] 的全部元素充零； ② 逐个单元地建立[ k ]e ，然后按单元 e 中局部编码与整体编码的对应关系， ③ 对所有单元完成上述叠加后，就形成了整体劲度矩阵[ K ] 。 * §6-1 基本量及基本方程的矩阵表示 §6-2 有限单元法的概念 §6-3 单元的位移模式与解答的收敛性 §6-4 单元的应变列阵和应力列阵 §6-5 单元的结点力列阵与劲度列阵 §6-6 荷载向结点移置 单元的结点荷载列阵§6-7 结构的整体分析 结点的平衡方程组 §6-8 解题的具体步骤 单元的划分 §6-9 计算成果的整理 §6-10 计算实例 第六章 用有限单元法解平面问题 有限单元法工程应用实例1 头盔撞击试验仿真模型与结果 有限单元法工程应用实例2 高强钢板厚度10mm，材料考虑应变率影响和失效，不受任何约束。模拟受初始速度为120m/s和180m/s钢球的冲击过程。 穿甲试验仿真 初速度为120m/s 初速度为180m/s 有限单元法工程应用实例3 动画显示的是地基中心点的沉降随线性荷载的变化过程，云图显示莫尔库仑材料的塑性区形成和大变形塑性流动过程。 群桩复合地基承载力计算结果 §6-1 基本量及基本方程的矩阵表示 体力列阵： 面力列阵： 应力列阵： 应变列阵： 位移列阵： 物理方程： 称为弹性矩阵。 对于平面应变问题，只需将弹性矩阵[D ]中的E、?分别换成 即可。 （平面应力问题） 则虚功方程可用矩阵表示为： 几何方程： 此外，用限单元法还要用到虚功方程： 现将虚位移及与该虚位移相应的虚应变表示为： 对连续变形体，它可以代替平衡微分方程和应力边界条件。 1、对连续体进行离散化。 §6-2 有限单元法的概念 有限单元法是用由有限多个、有限大小的单元在有限个结点相互连接的集合体来近似原来的连续体，当上述单元足够小从而划分网格足够密时，就可以真实地模拟原连续体。 有限单元法分析的基本步骤： 2、单元分析： （1）选择适当的位移模式，用单元结点位移（为基本未知量）来表示单元 内任一点的位移，即要建立如下关系式: 对于平面问题，最简单而常用的单元是三角形单元。在平面应力问题中，它们是三角板，在平面应变问题中，它们是三棱柱。 结点——铰接点 { d }e 称为单元结点位移列阵。 （2）应用几何方程，求出单元的应变，即： （3）应用物理方程，求出单元的应力，即： 其中[ S ]称为应力转换矩阵。 i j m ui vi um vm y x uj vj O 其中[ N ] 称为形函数矩阵。 其中[ B ]称为应变转换矩阵。 （5）将作用在单元上的外荷载按虚功相等的原则， 移置到单元各结点处，成为单元结点荷载： （4）由于单元产生了应力，则在单元的边界及内部作用有与之平衡的面力和 体力；现将其按虚功相等的原则移置到单元各个顶点处，作为结构其它 部分通过结点对此单元的作用力，称单元结点力，再利用虚功方程，得： 即为单元结点力，[ k ] 称为单元劲度矩阵。 对各结点进行平衡分析，列平衡方程并组集 得到整体结点平衡方程组： 其中[ K ]称整体劲度矩阵。 i j m O x Fix Fiy y Fmy Fmx Fjy Fjx 其中： 3、整体分析： §6-3 单元的位移模式与解答的收敛性 对三结点三角形单元，假设位移分量只是坐标的线性函数，即： 由左边三个方程求解a1 、 a2 、 a3 ，右边三个方程求解 a4 、 a5 、 a6 。再代回u、v 中，得 在i、j、m三个结点，位移应当等于结点位移，即： i j m ui vi um vm y x uj vj O 一、位移模式： 其中，Ni 、Nj、 Nm 称形函数，其表达式为 为单元 i j m 的面积。为使面积不致为负，在图示坐标系中i → j → m 的次序须是逆时针的。 而 i j m ui vi um vm y x uj vj O 分别为系数行列式第一、二、三列各元素的代数余子式。 则单元内任一点的位移可用矩阵表示为： 为单元结点位移列阵。 为形函数矩阵。 其中： 简写为 二、形函数的几何意义及性质： 记三角形单元 i j m 内的任一点为P (x , y)，