有限元-用三角形单元分析.pptVIP

第4章 平面三角形单元有限元 4.4 平面三角形单元分析 4.4 平面三角形单元分析 平面三角形单元整体分析 平面三角形单元整体分析 平面三角形单元整体分析 平面三角形单元整体分析 平面三角形单元整体分析 3.3.3 单元刚度矩阵的特性 页码: * School of Materials Science and engineering, WHUT 材料成形数值模拟 第4章 平面单元有限元法 Resolution: * ? School of Materials Science and engineering , WHUT * 第二步，把各单元的贡献矩阵叠加，即得出整体矩阵 Part 1 4.1 弹性力学平面问题的有限元法 弹性力学平面问题采用三角形单元。如图5-4a所示坝横截面图。将它划分为三角形网格，把原来连续体简化为由有限个三角形单元组成的离散体。三角形单元之间只在结点处用铰相连。 解题步骤：先进行单元分析，得出单元矩阵； 考虑单元综合，得出整体矩阵。因此，平面问题有限元法步骤： 离散化→单元分析→整体分析 4.4 平面三角形单元分析 4.4.1 平面三角形单元分析步骤 则三角形单元结点位移向量为： 以6个结点位移分量作为基本未知量，对应的物理量是六个结点力分量， 单元分析的主要任务是推导基本未知量与对应量之间的转换关系，即 (4-8) 例题 4-1 如图所示等腰直角三角形单元，试求刚度矩阵[k] 解：已知：直角三角形单元结点力与单元结点位移的关系如下： (4-26) 其中 单元刚度矩阵 （t——板厚） 弹性力学平面问题三角形单元的刚度矩阵具有对称性和奇异性等特点。 三角形单元刚度矩阵 可写成分块形式： (4-29) 子块[kij]是2×2阶矩阵，它是结点i的结点力子向量与结点j的位移子向量之间的刚度子矩阵。 4.4.4 整体分析步骤 整体分析包括下列四个主要步骤(见下图)： 现结合图5-21说明。图中离散化后的三角形网格，有4个单元和6个结点， 进行了编号。支承在结点4，5，6处有四个支杆支撑。 (1) 建立整体刚度矩阵 结构中有六个结点，共有12个结点位移分量和12个结点力分量。由结点位移求结点力： 分块形式 (4-40) 位移子向量和结点力子向量都是二阶向量，刚度子矩阵 是2×2阶矩阵。 矩阵[K]——整体刚度矩阵，是12×12阶矩阵。 (2) 引入支承条件： (3) 解方程组，求结点位移 (4) 根据结点位移求应力 建立整体 刚度矩阵 引入支 承条件 解方程 求位移 求应力 4.4.5 整体刚度矩阵的形成 用刚度集成法形成整体刚度矩阵． 图5-22表示结点的两种编码： 一是结点总码。六个结点统一编码为1，2，4，4，5，6。 二是结点局部码。每个单元的三个结点按反时针方向的顺序各自编码为1,2,4。 由图5-22，四个单元的局部码与总码的对应关系为 单元①：1，2,，4——1，2，4； 单元②：1，2,，4——2，4，5； 单元③：1，2,，4——5，4，2； 单元④：1，2,，4是6×6阶矩阵，分块形式为式(4-29)，其中9个子块按局部码排列。 整体刚度矩阵是12×12阶矩阵，分块形式为式(4-40)，其中36个子块按总码排列的． 平面三角形单元整体的刚度矩阵按刚度集成法分两步进行： 第一步，把单元刚度矩阵 扩大成单元贡献矩阵。这一步包含两个内容： ①阶数扩大——由6×6阶的扩大为12×12阶的贡献矩阵。 ②子块“搬家”——把中按局部码排列的9个子块“搬家”，变为 中按总码排列的9个子块。中的其余27个子块，则用零子块来填充。 以单元②为例，局部码1，2，4对应于总码2，4，5，按对应关系搬家，得出单元②的贡献矩阵 ： (4-41) 用同样方法可得出其它单元的贡献矩阵，，。 例1 用刚度集成法求图所示结构的整体刚度矩阵 ，设 μ=0。 解：1) 求各单元的刚度矩阵 图中的四个单元几何形状和尺寸相同，且结点编码方式也相同(逆时针方向顺序1，2，4，直角顶点的编码都是2)，因此，它们的单元刚度矩阵相同。当μ=0时，可直接用例1式(4-28)，即，单元①、②、③和④的刚度矩阵均为：分块形式如下： (4-28) 2) 求各单元的贡献矩阵 以单元②为例，贡献矩阵由式(4-41)求出： 同样可以求出贡献矩阵，，。 3) 求整体刚度矩阵 将各单元贡献矩阵叠加，由式(