有限元分析及应用

《有限元分析及应用》由会员分享，可在线阅读，更多相关《有限元分析及应用（277页珍藏版）》请在文档大全上搜索。

1、有限元分析及应用有限元分析及应用第一章 有限元法简介2有限元法介绍 有限元法的基本思想是将结构离散化结构离散化，用有限个容易分析的单元单元来表示复杂的对象，单元之间通过有限个结点结点相互连接，然后根据变形协调条件变形协调条件综合求解。由于单元的数目是有限的，结点的数目也是有限的，所以称为有限元法(FEM，Finite Element Method)。 3有限元法是最重要的工程分析技术之一。它广泛应用于弹塑性力学、断裂力学、流体力学、热传导等领域。有限元法是60年代以来发展起来的新的数值计算方法，是计算机时代的产物。虽然有限元的概念早在40年代就有人提出，但由于当时计算机尚未出现，它并未受到人们

2、的重视。 4随着计算机技术的发展，有限元法在各个工程领域中不断得到深入应用，现已遍及宇航工业、核工业、机电、化工、建筑、海洋等工业，是机械产品动、静、热特性分析的重要手段。早在70年代初期就有人给出结论：有限元法在产品结构设计中的应用，使机电产品设计产生革命性的变化，理论设计代替了经验类比设计。 5有限元法的孕育过程及诞生和发展 牛顿(Newton)莱布尼茨(Leibniz G. W.) 6大约在300年前，牛顿和莱布尼茨发明了积积分法分法，证明了该运算具有整体对局部的可加性。虽然，积分运算与有限元技术对定义域的划分是不同的，前者进行无限划分而后者进行有限划分，但积分运算为实现有限元技术准备好

3、了一个理论基础。 7在牛顿之后约一百年，著名数学家高斯提出了加权余值法及线性代数加权余值法及线性代数方程组的解法方程组的解法。这两项成果的前者被用来将微分方程改写为积分表达式，后者被用来求解有限元法所得出的代数方程组。 高斯(Gauss)8在18世纪，另一位数学家拉格朗日提出泛泛函函分析。泛函分析是将偏微分方程改写为积分表达式的另一途径。 拉格朗日(Lagrange J.)9在19世纪末及20世纪初，数学家瑞利和里兹（Rayleigh Ritz）首先提出可对全定义域全定义域运用展开函数运用展开函数来表达其上的来表达其上的未知函数未知函数。 瑞利(Rayleigh)101915年，数学家伽辽金(

4、Galerkin)提出了选择展开函数中形函数的伽辽金法伽辽金法，该方法被广泛地用于有限元。1943年，数学家库朗德第一次提出了可在定义域内分片地使用展开函数来表达其上的未知函数。这实际上就是有限元的做法。 1112（对象、变量、方程、求解途径）各力学学科分支的关系13(1) 桥梁隧道问题14任意变形体力学分析的基本变量及方程研究对象：任意形状的变形体几种典型的对象圆形隧道三维模型15(2) 中华和钟(3) 矿山机械16(4) 压力容器的成形17变形体及受力情况的描述18求解方法19有限元方法的思路及发展过程思路：以计算机为工具，分析任意变形体以获得所有力学信息，并使得该方法能够普及、简单、高效

5、、方便，一般人员可以使用。实现办法：20技术路线：21发展过程：如何处理对象的离散化过程22. .常用单元的形状常用单元的形状点点 (质量质量)线线(弹簧，梁，杆，间隙弹簧，梁，杆，间隙)面面 (薄壳薄壳, 二维实体二维实体,轴对称实体轴对称实体)二次二次体体(三维实体三维实体)线性线性二次二次.线性线性. . . . .23点点 单元单元线线 单元单元一维波传导问题一维波传导问题 24点点 单元单元线线 单元单元25XY00.020.040.060.080.10.12-0.1-0.08-0.06-0.04-0.020XY0.0540.0560.0580.06-0.003-0.002-0.00

6、10面面 单元单元28XY00.020.040.060.080.10.12-0.1-0.08-0.06-0.04-0.020XY0.0540.0560.0580.06-0.003-0.002-0.00102930受垂直载荷的托架受垂直载荷的托架31线性单元线性单元 / 二次单元二次单元 更高阶的单元模拟曲面的精度就越高更高阶的单元模拟曲面的精度就越高。低阶单元低阶单元更高阶单元更高阶单元体体单元单元32 有限元分析的作用有限元分析的作用l 复杂问题的建模简化与特征等效复杂问题的建模简化与特征等效l 软件的操作技巧（单元、网格、算法参数控制）软件的操作技巧（单元、网格、算法参数控制）l 计算结果

7、的评判计算结果的评判l 二次开发二次开发l 工程问题的研究工程问题的研究l 误差控制误差控制36第二章 有限元分析的力学基础 2.1 变形体的描述与变量定义变形体的描述与变量定义(1) 变形体 变形体：即物体内任意两点之间可发生相对移动。 有限元方法所处理的对象：任意变形体38(2) 基本变量的定义 可以用以下各类变量作为任意变形体的描述因此，在材料确定的情况下，基本的力学变量应该有：位移、应变、应力量39目的：对弹性体中的位移、应力、应变进行定义和表达，进而建立平衡方程、几何方程和材料物理方程(3) 研究的基本技巧采用微小体积元dxdydz的分析方法（针对任意变形体）402.2 弹性体的基本

8、假设弹性体的基本假设为突出所处理的问题的实质，并使问题简单化和抽象化，在弹性力学中，特提出以下几个基本假定。(1) 物质连续性连续性假定：物质无空隙，可用连续函数来描述；(2) 物质均匀性均匀性假定：物体内各个位置的物质具有相同特性；(3) 物质(力学)特性各向同性各向同性假定：物体内同一位置的物质在各个方向上具有相同特性；(4) 线性弹性线性弹性假定：物体的变形与外来作用的关系是线性的，外力去除后，物体可恢复原状；(5) 小变形小变形假定：物体变形远小于物体的几何尺寸，在建立方程时，可以高阶小量（二阶以上）。 以上基本假定将作为问题简化的出发点。412.3 基本变量的指标表达基本变量的指标表

9、达指标记法的约定：自由指标自由指标：在每项中只有一个下标出现，如 ，i，j为自由指标，它们可以自由变化；在三维问题中，分别取为1，2，3；在直角坐标系中，可表示三个坐标轴x, y, z。哑指标哑指标：在每项中有重复下标出现，如： ，j为哑指标。在三维问题中其变化的范围为1,2,3ijijijbxa42Einstein 求和约定：哑指标意味着求和指标记法的应用：指标记法的应用：对于方程组按一般的写法，可写为若用指标记法：(2-3)式与(2-2)式等价，因为j为哑指标，意味着求和（2-1）（2-2）（2-3）43克罗内克符号克罗内克符号 在笛卡尔直角坐标系下，由ij表示的Kronecker(克罗内

10、克)符号定义为 jijiij如果如果 , 0 , 1亦即1332211023321331211244那么，矩阵 333231232221131211100010001= 是单位矩阵。根据上述定义，可以推出下列关系 3332211ii333323213132323222121213132121111aaaaaaaaaaaaaaajjjjjj45弹性力学里假想把物体分成无限多微小六面体 ，称为微元体。考虑任一微元体的平衡（或运动），可写出一组平衡（或运动）微分方程及边界条件。但未知应力的数目总是超过微分方程的数目，所以弹性力学问题都是超静定的，必须同时考虑微元体的变形条件以及应力和应变的关系，它们

11、在弹性力学中相应的称为几何方程和物理方程。平衡平衡（或运动）方程、几何几何方程和物理物理方程以及边界条件，称为弹性力学的基本方程。2.4 弹性力学的基本方法弹性力学的基本方法46从取微元体入手，综合考虑静力（或运动）、几何、物理三方面条件，得出其基本微分方程，再进行求解，最后利用边界（表面）条件确定解中的常数，这就是求解弹性力学问题的基本方法。472.5 空间问题的基本方程空间问题的基本方程dydxdz483D情形下的力学基本变量将正应力和正应变简写成49abbaaddccxyxyyxyxyzyzzyzyzxzxxzxz50由力平衡条件0X有：0Xdxdydzdxdydxdydzzdxdzdx