FEM-1-有限元方法概论

《FEM-1-有限元方法概论》由会员分享，可在线阅读，更多相关《FEM-1-有限元方法概论（44页珍藏版）》请在金锄头文库上搜索。

1、有限元分析方法有限元分析方法 龚峰龚峰 gongfeng gongfeng 2 参考书目 3 考核方式 考勤：10% 作业：30% 课程报告：60% 4 第1章 有限元方法概论 1.1 有限元方法形成的背景 1.2 有限元方法的基本原理 1.3 有限元分析主要应用领域 1.4 常用有限元分析软件介绍 1.5 有限元分析的作用及应用 5 1.1 有限元方法形成的背景 两类典型的工程问题 第一类问题，可以归结为有限个已知单元体的 组合。例如，材料力学中的连续梁、建筑结构 框架和桁架结构。 平面桁架结构，由6个承受 轴向力的“杆单元”组成。 1889年建成的Effiel铁塔 ，由18036个部件组成

2、 6 1.1 有限元方法形成的背景 两类典型的工程问题 第二类问题，通常可以建立它们应遵循的基本方 程，即微分方程和相应的边界条件。例如弹性力 学问题，热传导问题，电磁场问题等。热传导问 题的控制方程与换热边界条件如下： 7 1.1 有限元方法形成的背景 两类典型的工程问题 l 第一类问题的研究对象称为离散系统。离散系 统是可解的，但是求解复杂的离散系统，要依 靠计算机技术。 l 第二类问题的研究对象称为连续系统。可以建 立描述连续系统的基本方程和边界条件，通常 只能得到少数简单边界条件问题的解析解。对 于大多数实际的工程问题，需要用近似算法来 求解。 8 1.1 有限元方法形成的背景 两类典

3、型的工程问题 为了解决这个困难，工程师和数学家开始寻找一 种近似的求解方法，在这个过程中，他们从两个 不同的路线得到了相同的结果，即有限单元法。 结构分析的有限元方法是由一批工业界和学术 界的研究者在二十世纪五十年代到二十世纪六十 年代创立的。 9 1.1 有限元方法形成的背景 工程师的贡献 l 思路来源于固体力学结构分析矩阵位移法的发展和工程师对 结构相似性的直觉判断。对于不同结构的杆系、不同的载荷 ，求解时都能得到统一的矩阵公式。从固体力学的角度看， 桁架结构等标准离散系统与人为地分割成有限个分区的连续 系统在结构上存在相似性，可以把结构分析的矩阵法推广到 非杆系结构的求解。 l 1956

4、年，波音公司的Turner, Clough, Martin, Topp在纽约举 行的航空学会年会上介绍了将矩阵位移法推广到求解平面应 力问题的方法，即把结构划分成一个个三角形和矩形“单元” ，在单元内采用近似位移插值函数，建立了单元节点力和节 点位移关系的单元刚度矩阵，并得到了正确的解答。 l 1960年，Clough在他的名为“The finite element in plane stress analysis”的论文中首次提出了有限元这一术语。 10 1.1 有限元方法形成的背景 数学家的贡献 l 数学家们则发展了微分方程的近似解法，包括有限差分 方法，变分原理和加权余量法。 l 195

5、4-1955年，斯图加特大学的Argyris在航空工程杂志 上发表了一组能量原理和结构分析论文，为有限元研究 奠定了重要的基础。 l 1963年前后，经过J. F. Besseling, R.J. Melosh, R.E. Jones, R.H. Gallaher, T.H.H. Pian等许多人的工作，认 识到有限元法就是变分原理中Ritz近似法的一种变形， 发展了用各种不同变分原理导出的有限元计算公式。 11 1.1 有限元方法形成的背景 数学家的贡献 l 1965年O.C.Zienkiewicz和Y.K.Cheung发现只要能写成 变分形式的所有场问题，都可以用与固体力学有限元法 的相同

6、步骤求解。 l 1967年，Zienkiewicz和Cheung出版了第一本有限元分 析的专著。 l 1969年B.A. Szabo和G.C. Lee指出可用加权余量法特别 是Galerkin法，导出标准有限元过程来求解非结构问题 。 l 1970年以后，有限元方法开始应用于处理非线性和大变 形问题，Oden于1972年出版了第一本关于处理非线性 连续体的专著。 12 1.1 有限元方法形成的背景 我国学者的贡献 l 陈伯屏结构矩阵方法 l 钱伟长广义变分原理 l 胡海昌广义变分原理 l 冯 康有限单元法理论 13 1.2 有限元方法的基本原理 有限元方法是求解数学物理问题的一种数值计算 方法

7、，起源于固体力学，然后迅速扩展到流体力 学、传热学、电磁学等其他物理领域。 有限元分析是利用数学近似的方法对真实物理系 统（几何和载荷工况）进行模拟。利用简单而又 相互作用的元素，即单元，用有限数量的未知量 去逼近无限未知量的真实系统。 有限元模型是真实系统理想化的数学抽象。 14 1.2 有限元方法的基本原理 真实系统有限元模型齿轮有限元模型 有限元模型由一些简单形状的单元组成，单元之间 通过节点连接，并承受一定载荷。节点具有一定的 自由度。 15 1.2 有限元方法的基本原理 自由度(DOFs) 用于描述一个物理场的响应特性。 分析对象 自由度 结构 位移 热 温度 电 电位 流体 压力

8、磁 磁位 16 1.2 有限元方法的基本原理 有限元分析一般过程 l 结构离散将连续系统分割成有限个分区或单元。 单元类型选择，单元划分，节点编码。 l 单元分析用标准方法对每个单元提出一个近似解。 选择位移函数，分析单元的特征。 l 整体分析将所有单元按标准方法组合成一个与原有 系统近似的系统。 集成整体节点载荷向量，集成整体刚度方程，引进边界约束条件。 分割组合 17 1.2 有限元方法的基本原理 圆的面积计算 18 1.2 有限元方法的基本原理 有限元方法的基本思想和原理是“简单”而“朴素” 的，在发展初期，许多学术权威对该方法的学术 价值有所轻视，国际著名刊物Journal of Ap

9、plied Mechanics许多年来拒绝刊登有关有限元方法的文 章，其理由是没有新的科学实质。 现在完全不同了，由于有限元方法在科学研究和 工程分析中的地位，有关有限元方法的研究已经 成为数值计算的主流。涉及有限元方法的杂志有 几十种之多。 19 1.2 有限元方法的基本原理 有限元方法的主要研究内容： l计算方法：大型线性方程组的解法，非线性问题 的解法，动力问题计算方法。 l高精度单元 l复杂材料模型 l多物理场耦合 目标：提高计算效率和计算精度 20 1.3 有限元分析主要应用领域 结构分析 热分析 电磁分析 流体分析 耦合场分析 21 1.3 有限元分析主要应用领域 结构分析是有限元

10、分析方法最常用的一个应用领域 。结构是一个广义的概念，它包括土木工程结构， 如桥梁和建筑物；汽车结构，如车身骨架；海洋结 构，如船舶结构；航空结构，如飞机机身等；同时 还包括机械零部件，如活塞，传动轴等。 结构分析中计算得出的基本未知量（节点自由度） 是位移，其他的一些未知量，如应变，应力，和反 力可通过节点位移导出。 22 1.3 有限元分析主要应用领域 l 静力分析 ：用于静态载荷. 可以考虑结构的线性及非线性行 为，例如: 大变形、大应变、应力刚化、接触、塑性、超弹 及蠕变等. l 动力分析 ：动力学分析是用来确定惯性（质量效应）和阻 尼起着重要作用时结构或构件动力学特性的技术。“动力学

11、 特性” 可能指的是下面的一种或几种类型: u振动特性 （结构振动方式和振动频率） u周期（振动）载荷的效应 u随时间变化载荷的效应 l 屈曲分析 ：用于计算屈曲载荷和确定屈曲模态。包括线性 （特征值）和非线性屈曲分析。 23 1.3 有限元分析主要应用领域 转向机构支架的强度分析（MSC/Nastran） 24 1.3 有限元分析主要应用领域 粘性介质压力成形过程分析（DEFORM） 25 1.3 有限元分析主要应用领域 齿轮锻造过程分析（DEFORM ） 26 1.3 有限元分析主要应用领域 拉深成形过程分析（DYNAFORM ） 27 1.3 有限元分析主要应用领域 拉深成形过程分析（D

12、YNAFORM ） 28 1.3 有限元分析主要应用领域 拉深成形过程分析（DYNAFORM ） 29 1.3 有限元分析主要应用领域 动力分析的分类 l 模态分析：计算线性结构的自振频率及振形。 l 谱分析：是模态分析的扩展，用于计算由于随机振动引起的结构应力和应变。 l 谐响应分析：确定线性结构对随时间按正弦曲线变化的载荷的响应. u 旋转设备（如压缩机、发动机、泵、涡轮机械等）的支座、固定装置和部 件； u 受涡流（流体的漩涡运动）影响的结构，例如涡轮叶片、飞机机翼、桥和 塔等。 l 瞬态动力学分析：确定结构对随时间任意变化的载荷的响应。可以考虑与静力 分析相同的结构非线性行为。 l 显

13、式动力分析：计算高度非线性动力学和复杂的接触问题。 u 用于模拟非常大的变形，惯性力占支配地位，并考虑所有的非线性行为。 u 显式求解冲击、碰撞、复杂金属成形等问题，是目前求解这类问题最有效 的方法。 30 隔热屏的模态分析(ANSYS) 1.3 有限元分析主要应用领域 31 1.3 有限元分析主要应用领域 热分析 l 热分析在许多工程应用中扮演重要角色，如内燃机 、涡轮机、换热器、管路系统、电子元件等。 l 热分析之后往往进行结构分析,计算由于热膨胀或 收缩不均匀引起的应力。 l 热相关问题 u相变 (熔化及凝固), 内热源 (例如电阻发热等) u三种热传递方式 (热传导、热对流、热辐射)

14、u稳态传热：系统的温度场不随时间变化 u瞬态传热：系统的温度场随时间明显变化 热分析计算物体的稳态或瞬态温度分布，以及热量的获 取或损失、热梯度、热通量等。 32 1.3 有限元分析主要应用领域 热分析 工件淬火3.06 min 时的温度、组织分布 33 1.3 有限元分析主要应用领域 电磁分析 l 磁场分析中考虑的物理量是磁通量密度、磁场密度、磁力、磁力矩、阻 抗、电感、涡流、能耗及磁通量泄漏等。磁场可由电流、永磁体、外加 磁场等产生。 l静磁场分析 ：计算直流电或永磁体产 生的磁场。 l交变磁场分析 ：计算由于交流电产生 的磁场。 l瞬态磁场分析：计算随时间随机变化的 电流或外界引起的磁场

15、。 l电场分析：用于计算电阻或电容系统的 电场。典型的物理量有电流密度、电荷 密度、电场及电阻热等。 l高频电磁场分析：用于微波及RF无源 组件，波导、雷达系统、同轴连接器等 分析。 电磁成形分析(ANSYS) 34 1.3 有限元分析主要应用领域 流体分析 用于确定流体的流动及热行为。 可以处理不可压缩或可压 缩流体、层流及湍流，以及多组份流等。 l 作用于气动翼型上的升力和阻力 l 超音速喷管中的流场 l 弯管中流体的复杂的三维流动 隔热屏分析(ANSYS) 35 1.3 有限元分析主要应用领域 流体分析 超音速飞行压力分布汽车气动分析 高速导弹气动 36 1.3 有限元分析主要应用领域

16、耦合场分析 考虑两个或多个物理场之间的相互作用。如果两个物理场之间相互 影响，单独求解一个物理场是不可能得到正确结果的，因此你需要 一个能够将两个物理场组合到一起求解的分析软件。 例如： 在压电力分析中，需要同时求解电压分布（电场分析）和应 变（结构分析）. 需要耦合场分析的典型情况有： l热应力分析 l流体结构相互作用 l感应加热（电磁热）, 感应振荡 两根热膨胀系数不同的棒焊接 在一起，加热后的变形情况 37 1.4 常用有限元分析软件介绍 常用大型通用有限元软件 ABAQUS、ADINA、ANSYS、MSC/Marc MSC/Nastran 一些专用有限元软件 Autoform、DEFORM、DYNAFORM、LS-DYNA MOLDFLOW、MSC/Dytran、PAM-CRASH PAM-STAMP、ProCast、SysWeld 38 1.4 常用有限元分析软件介绍 AB