有限元法与程序-平面问题6概要.ppt

§ 2.9 平面有限元分析实施步骤与注意事项 将要计算的弹性体划分成三角形单元。对结点进行编号，列出结点坐标作为输入信息。 对单元进行编号，列出单元三个结点的号码作为输入信息。 计算载荷的等效结点力，把等效结点力作为输入信息。 计算各单元的常数， 再计算2Δ。 § 2.9 平面有限元分析实施步骤与注意事项 计算各单元的刚度矩阵 。 形成整体刚度矩阵。 处理约束及消除刚体位移。 解线性方程组，求结点位移 。 计算应力矩阵，再计算单元应力。根据需要计算主应力和主方向。 § 2.9 平面有限元分析实施步骤与注意事项 对称性的利用。 结点的选择和单元划分。 结点的编号。 单元结点i、j、m的次序。 边界条件处理和整体刚度矩阵的修正 应力计算结果的整理 § 2.9 平面有限元分析实施步骤与注意事项 § 2.9 平面有限元分析实施步骤与注意事项 § 2.10 平面问题的计算机程序 § 2.10 平面问题的计算机程序 § 2.10 平面问题的计算机程序 § 2.10 平面问题的计算机程序 § 2.10 平面问题的计算机程序 § 2.10 平面问题的计算机程序 例：计算下图半带宽。 作业: （二选一）1) 写出计算方法中LDLT分解算法及fortran程序2) 用Matlab或Mathmatic软件确定十结点三角形单元的应变矩阵应力矩阵及单元刚度矩阵 * 第2章 平面问题的有限元法 上次课内容回顾 平面矩形单元 平面6节点三角形单元 有限元法分析思路流程 解综合方程[K]{⊿}= {P} 求结构节点位移{⊿} 计算结构内力和应力 系统分析 (把单元刚度矩阵集合成结构刚度矩阵[K] 形成等价节点荷载{P} ) 离散（剖分）结构 为若干单元 单元分析 (建立单元刚度矩阵[k]e 形成单元等价节点力) 以三角形常应变单元解平面问题为例，具体歩骤可归纳如下 一、实施步骤 二、注意事项 二、注意事项 边界条件的处理：划0置1法 边界条件 二、注意事项 边界条件的处理：冲大数法 一、有限元分析程序基本框图 开始 输入基本数据 计算单元刚度矩阵 形成总体刚度矩阵 形成结点荷载向量 引入约束条件 求解方程组，输出结点位移 计算单元应力，输出结果 结束 1、输入基本数据（结构描述）： （1）控制数据：如结点总数、单元总数、约束条件总数等； （2）结点数据：如结点编号、结点坐标、约束条件等； （3）单元数据：如单元编号、单元结点序号、单元的材料特性、几何特性等； （4）载荷数据：包括集中载荷、分布载荷等。 一、有限元分析程序基本框图 开始 输入基本数据 计算单元刚度矩阵 形成总体刚度矩阵 形成结点荷载向量 引入约束条件 求解方程组，输出结点位移 计算单元应力，输出结果 结束 （1）各单元的bi,ci(i,j,m) ， 面积A； （2）应变矩阵[B]，应力矩阵[S]； （3）单元刚度矩阵[k]； （4）单元等价载荷列向量[F]。 2、单元分析 3、系统分析 （1）整体刚度矩阵[K]的组装； （2）整体载荷列阵{P}的形成； 关键问题：[K]的存储； 一、有限元分析程序基本框图 开始 输入基本数据 计算单元刚度矩阵 形成总体刚度矩阵 形成结点荷载向量 引入约束条件 求解方程组，输出结点位移 计算单元应力，输出结果 结束 4、约束引入 5、线性方程组求解 求解方法常用：GAUSS消元法，LDLT、QR分解法等。其程序在一些专著中列出（例如：徐士良编。FORTRAN常用算法程序集。清华大学出版社）。在此不作详细介绍，其方法参阅有关书籍。 一、有限元分析程序基本框图 开始 输入基本数据 计算单元刚度矩阵 形成总体刚度矩阵 形成结点荷载向量 引入约束条件 求解方程组，输出结点位移 计算单元应力，输出结果 结束 6 应力计算结果的整理 计算结果包括位移和应力两个方面。在位移方面，一般无须进行整理工作。应力结果则需要整理。通常认为计算出的应力是三角形单元形心处的应力。而相邻单元之间的应力存在突变，甚至正、负符号都不相同。为了由计算结果推算出结构内某一点的接接实际的应力，必须通过某种平均计算。通常可采用两单元平均法或绕结点平均法。 一、有限元分析程序基本框图 开始 输入基本数据 计算单元刚度矩阵 形成总体刚度矩阵 形成结点荷载向量 引入约束条件 求解方程组，输出结点位移 计算单元应力，输出结果 结束 6 应力计算结果的整理 两单元平均法：把两个相邻单元中的常应力加以平均，用来表示公共边界中点处的应力。 绕结点平均法：把环绕某一结点的各单元常应力加以平均，用以表示该结点的应力。在内结点效果较好，而在边界结点可能很差，一般改为应由内结点的应力外推计算出来。 二、总体刚度矩阵