有限元与数值方法-讲稿4研究.ppt

State Key Laboratory of Structural Analysis for Industrial Equipment Dalian University of Technology, China 有限元与数值方法第四讲第4章 有限元法的一般原理 杆单元的有限元分析 位移插值 位移插值 位移及应变 单元刚度阵 平面内任意方向的杆单元 边界条件处理 有限元与数值方法第五讲第4章 有限元法的一般原理（续） 杆单元的有限元分析 杆单元的有限元分析：位移及应变 单元刚度阵 平面内任意方向的杆单元 由梁组成的结构 欧拉-伯努利梁的基本假定 弯曲梁的有限元分析 单元内位移插值 单元内位移插值 单元内位移插值 单元的应变能 总体平衡方程 有限元与数值方法第五讲第4章 有限元法的一般原理（续） 总刚度矩阵的特性 总刚度矩阵的特性 总刚度矩阵的对称性 总刚度矩阵的对称性 总刚度矩阵的稀疏性 总刚度矩阵的组装和存贮 总刚度矩阵的组装和存贮 总刚度矩阵的组装和存贮 结构对称性的利用 结构对称性的利用 有限元法的实施流程 有限元法的理论基础概述 将微分方程转化为等效积分弱形式 变分原理 加权余量法 采用单元上的分片假设近似函数，将积分方程转化为代数方程组 有限元法的基本概念 有限元分析的基本思想是将求解域场分成小的子区域，通常称为“单元”或“有限元”。 对每一单元假定一个分片近似解，然后推导求解这个域总的满足条件(如结构的平衡条件），从而得到问题的解。 有限元法方程的系数矩阵通常是稀疏的，便于求解。 有限元法不仅计算精度高，而且能适应各种复杂形状，不同物理特性、多变的边界条件和任何承载情况的工程结构分析问题。 有限元法应用于场（力场、电场、磁场、温度场、流体场等）分析、热传导、非线形材料的弹塑性蠕变分析等 有限元法总体思路 解方程 Kd= F 提取单元位移de 计算单元内力和应力 把单元刚度矩阵集合成结构刚度矩阵K 结构离散为单元 建立单元刚度矩阵Ke 形成等价节点荷载F 形成单元等价节点力 有限元的基本理论基础 有限元法(FEM)是求解偏微分方程边值问题近似解的数值方法 边值问题 未知量是由控制方程（椭圆、双曲、抛物型）描述的场变量（如位移、温度、流体速度等） 边界条件是给定的场变量值或者其偏导数 有限元法的基本概念 (a) 二维问题的几何域 (b) 三角形单元 (c) 有限元网格的一部分 单元 有限元网格 有限元法中的离散 各种几何形状的有限元单元 三角形的顶点称为节点（node） 节点处的场变量（这里是温度）将作为自变量被直接求解 node 热传导问题的三角形单元 node 有限元法中的场变量表示 以平面热传导问题的三角形单元为例 除了节点外的其他各位置的点对应的场变量如何确定？ 单元内部点的场变量值由单元节点的插值(interpolation )给出: T = ? 有限元法中的场变量表示 , , 和 是插值函数，称为位移函数或形函数。插值函数所包括的多项式阶数越高，越能精确表示位移分布。 常见平面单元形状与节点数 三节点三角形单元 CST(常应变单元) 六节点三角形单元 二次插值 八节点四边形单元 二次插值 CST三角形单元网格划分简单，但对于弯曲过刚；线性应变三角形元描写弯曲性能远优于CST单元 四边形单元剖分有时比较困难，但性能较好 四节点四边形单元 双线性插值 四节点四面体单元 线性插值（常应变） 十节点四面体单元 二次插值（线性应变） 八节点四面体单元 Lagrange单元 非完全三次插值 二十节点Serendipity单元 四面体单元网格剖分简单，但四节点四面体精度较差 八面体单元精度较好，但网格剖分比较困难 常见三维单元形状与节点数 一维单元 (x) (x) (x) 不同形式的单元插值 二维单元 不同形式的单元插值 三维单元 不同形式的单元插值 (x,y,z) (x,y,z) 有限元法通过加权余量法（或变分法、最小势能原理、虚功原理等）将偏微分方程转变为代数方程，便于计算机处理 将求解域剖分为网格，对节点变量进行单元分片插值值 单元矩阵形成规范化，形函数只与坐标有关，便于计算机计算 单元刚度阵组装为整体系数矩阵后，考虑边界条件，求解矩阵方程即可得到节点未知量。 Finite Difference (FD) Method: FD 对微分算子进行近似 在节点上建立方程 Finite Element (FE) Method: FE 采用精确的算子，但利用基函数对场变量进行近似 在问题域上建立方程（需要积分形式） 有限元法的基本思路及比较 Strong form Weak form Galerkin approx. Matrix for