第1章 概述(有限元分析).ppt

有限元法《有限元法及ANSYS程序应用基础 》 李仁军 第一节 前言 关于有限元法 英文缩写FEM（Finite Element Method） 应用中习惯称有限元分析 是一种连续结构离散化数值计算方法 上世纪五十年代由美国飞机工程师提出（1956年，Turner） FEM与CAE CAE－计算机辅助工程（Computer Aided Engineering） CAE范围更广，还包含其它工程分析方法 有限元法的应用 传统设计 列车车厢整体结构的有限元模型 北京奥运场馆鸟巢的钢铁枝蔓结构 空 客 A350 后 机 身 第 19 框 的 设 计 与 有 限 元 分 析 过 程 人 体 肩 部 区 域 的 骨 胳 有 限 元 分 析 模 型 及 计 算 结 果 半导体芯片温度场的数值仿真 弹丸超高速碰撞薄板 人体股骨端受力分析 有限元法的产生 1943年，Courant尝试应用定义在三角区域上的分片连续函数和最小位能原理相结合，求解St. Venant扭转问题。 现代有限元法的第一个成功尝试是Turnet,Clough等人在分析飞机结构时于1965年得到的成果，是将刚架位移法推广应用于弹性力学平面问题。 1960年，Clough进一步处理了平面弹性问题，并第一次提出了“有限单元法”的名称。 Turner等人最早提出有限单元法时是利用直接刚度法，它来源于结构分析的刚度法，只能处理一些比较简单的实际问题。 1963-1964年，Besseling,Melosh和Jones等人证明了有限单元法是基于变分原理的里兹法的另一种形式。 1960年后，随着电子计算机的发展，有限单元法的发展速度才显著加快。 第二节 有限元法的产生与基本思想 工程中遇到的数学方程可以分为两种： 1）代数方程 2）微分方程 （例如：悬臂梁弯曲方程） 数学方程的求解方法 1）解析法（精确解） 利用严格的数学推导进行求解 2）数值法（近似解） 通过一定算法利用计算机进行求解 ┏ 差分法 常用的数值方法 ╋ 变分法 ┗ 有限元法 差分法 基本思路：用均匀的网格离散求解域，用离散点的差分代替微分，从而将微分方程的求解转化为线性代数方程组的求解。 例1：求解区间[a,b]上的一维函数y(x)，且y(x)满足 解：首先将[a,b]划分为n等分，相邻两节点之间的距离h称为步长，且h=(b-a)/n。 n+1个节点上的函数值就是待求解的未知数 导数的数学定义： 因此，可以用差分近似代替微分得： 在第1至第n-1个节点处将(1-2)及(1-3)代入(1-1)可得n-1个线性代数方程： 再利用边界条件有 这样就可以求得未知函数y(x)在n+1个节点上的近似值，而位于节点间的函数值则可以利用节点上的函数值线性插值求得。 变分法 基本思路：微分方程边值问题的解等价于相应的泛函极值问题的解。 也就是使得未知函数泛函取得驻值的y(x)就是方程的解。 *里兹（Ritz）法 是从一族假定解中寻求满足泛函变分的“最好解” 因此，解的精确度取决于“试探函数”的选取。 边值问题的泛函表达式 对于边值问题： 其原问题的泛函可以表达为 例2：求解以下边值问题 解：将方程代入(1-4)求得其泛函 设试探函数可以表示为： 则有： 根据多元函数分析，泛函取驻值的必要条件是其对各自变量的偏导均为零，如此可得： 这样,就可以得到由n个关于 的线性方程组成的方程组，求解该方程组就可以求出这n个待定系数，再回代到(1-6)就可以求出y的近似表达式。 例3：求解以下边值问题 根据（1-4）可得泛函为： 1）选取试探函数： 将其代入（1-8）进行计算后得泛函为： 然后使泛函对a的偏导等于零，有： 2）选取试探函数： 将其代入（1-8）进行计算后得泛函为： 然后使泛函对a的偏导等于零，有： 有限元法 在差分法和变分法的基础上建立起来的，结合了两种方法的优点。 比较与分析： 1）离散 差分法： 对方程进行离散 网格规则 有限元法：对物理模型进行离散 网格不规则