有限元方法及其工程应用

《有限元方法及其工程应用》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、材料学院2002级研究生张可荣20021033指导教师：陈芙蓉《有限元方法及其工程应用》读书报告数值分析技术是在力学理论、计算数学和计算机技术相互结合和渗透的基础上发展起来的一门应用数学学科。它主要借助计算机和软件技术实现大规模的计算分析。根据构造数值计算公式的原理不同，目前工程上常用的数值分析方法主要有：有限差分法、有限单元法和边界单元法等。与其它方法比较，有限元法在计算公式构造、计算精度及效率、求解过程的稳定性和适用性等方面具有明显的优势。有限元法的基本思想是把一个复杂实际问题划分成有限个简单问题的组合进行求解，由于实际问题已被较简单的问题所代替，故只能获

2、得近似解。如对结构受力分析问题，首先把结构的求解区域看成是由有限（数量）个小的在节点处相互联系的子域（单元）组成，先对每一个单元假定一个合适的近似解，然后推导结构的整体平衡方程，在满足边界条件情况下就可获得近似解。当划分的子域（单元）尺寸变的越来越小时，其近似解就越来越逼近精确解。弹性力学是进行工程结构承载分析的基本理论。建立与未知量相等的方程是进行应力分析的首要条件，此外还需满足协调方程（位移和应变连续）和边界条件（弹性结构表面的给定位移和力的条件）。弹性力学假设物体是完全弹性、连续、均匀和各向同性的，并且变形和位移是微小的。弹性力学有外力、应力、应变和位移

3、等基本概念。弹性平面问题主要有平面应力问题和平面应变问题。平面应力问题主要应用于厚度尺寸与长度和宽度相比很小的板状结构体，如板架、机体等。这类物体只在板边受平行于板平面的外力，且外力沿厚度方向不变，体力也平行于板面并且不沿厚度变化。平面应力问题只有σXX，σYY，τXY=τYX三个应力分量不为零，是一种二维函数问题。平面应变问题适用于截面不变化但长度很长的柱形结构体，如长圆柱体、高压容器、管道等。这类物体只受到平行于截面、并且沿长度不变化的体力和面力。平面应变问题只有三个应变分量：εXX，εYY，γXY=γYX不为零。弹性力学的控制方程有：平衡微分方程、几何方

4、程和物理方程。其中弹性平面的平衡微分方程为：几何方程为应变和位移的关系：物理方程为应力-应变关系（即三维条件下的广义虎克定律）：，其它两个物理方程类似。另外还有变形协调方程和边界条件。可见三维弹性问题总共有15个未知参数。能量原理是力学的基本原理之一，弹性力学能量原理，就是利用能量的概念研究物体在外力的作用下应力、应变和位移参量之间的变化规律，以及外力作功与物体变形势能所涉及的能量转换过程。主要有泛函、变分的概念和虚位移原理和最小势能原理。在工程中除存在依赖与自变量变化的函数关系外，还存在另一类函数，其自变量也是一类函数，而不是有限个变量，这种函数的函数叫“泛

5、函”。变分学就是研究这些“泛函”的极值性质，即在一组容许函数中选定一个函数，使给定的“泛函”4材料学院2002级研究生张可荣20021033指导教师：陈芙蓉获得极值，而不是求解含有有限个变量的函数极值。在分析结构静力问题时，总认为其加载过程是缓慢的逐步加载过程，在该过程中结构始终处于静力平衡状态。当结构为线弹性体时，外力所作的功将全部转换为弹性体中的弹性势能，这一过程遵循能量守恒原理。由于外力作用，弹性体内部所储存的弹性变形势能成为应变势能。弹性体几何允许的微小位移，称为位移的变分，根据能量守恒原理，在虚位移过程中，外力所作的虚功等于弹性体内增加的虚应变能，这

6、就是虚位移原理。对任何一个弹性结构，在外力作用下产生弹性变形，在这一变形过程中，外力对弹性体作功W而势能减小，弹性体由于变形而储存弹性应变势能U，那么当该弹性体处于稳定平衡状态时，其总的势能必取极小值，这就是最小势能原理。利用有限元法进行结果分析时，首先要将物体或求解区域划分为有限个小区域（网格）的集合，这些小区域称“有限单元”，这些单元通过节点而相互连接，单元节点位于单元边界上，且相临单元通过节点连接。其次，假设节点处的物理参量（如位移、温度等）为未知参量，用这些参量可以表示单元内部物理量变化的近似函数。然后，利用结构的整体控制方程（平衡方程）将所有单元节点

7、的未知量集合成一个线性代数方程组，代入边界条件，最后求解。有限元方法的一般步骤是：1）物体结构离散化；2）选择单元内部插值模式；3）推导单元刚度矩阵和单元节点载荷向量；4）集合单元方程得到总的平衡方程组；5）求解节点未知参量；6）计算单元应变和应力。其中，总的平衡方程组为：[K]{U}={F}式中：[K]为整体结构的集合刚度矩阵；{U}为整体结构节点位移向量；{F}为节点载荷向量。有限元方法的第一步是对求解区域划分网格，网格为具有一定几何形状的小区域，节点在这种小区域的角点或边界上，节点的“自由度”是指节点处独立的物理参量的个数。整个求解区域的自由度数则是全部

8、节点自由度的总和，整个求解区域的自由度