数值模拟与仿真基础——1 绪论.pptVIP

工程中的数值模拟与仿真绪论 麻向军 华南理工大学 聚合物新型成型装备国家工程研究中心 聚合物成型加工工程教育部重点实验室 1 绪 论 1.1、有限元法介绍 1.2、CAE技术的应用 1.3、关于本课程的说明 1.1、有限元法介绍 1、有限元法的一般概念 数理方程的求解方法 有限差分法： 特 点：直接求解基本方程和相应定解条件的近似解。 求解步骤：将求解域划分网格，在网格节点上用差分方程代替微分方程。 应 用：流动与传热 缺 点：几何形状复杂时，为提高精度，需划分较多的网格。 有限单元法： 特 点：基于变分原理，建立基本方程和定解条件的积分提法，建立泛函的变分。 求解步骤：将求解域划分单元，在单元内假设近似分布函数，通过单元分析进而建立整个区域场量的方程。 应 用：几乎所有领域 优 点：几何形状复杂时，通过合适的网格数得到较高精度。 有限元法的基本思想 将一个连续的求解域 ( 连续体 ) 离散化即分割成彼此用节点 (离散点) 互相联系的有限个单元，在单元体内假设近似解的模式，用有限个结点上的未知参数表征单元的特性，然后用适当的方法，将各个单元的关系式组合成包含这些未知参数的代数方程，得出各结点的未知参数，再利用插值函数求出近似解。 有限元法是一种用有限的单元离散连续体，然后进行求解的一种数值计算的近似方法。由于单元可以分割各种形状和尺寸的连续体，所以，它能很好的适应复杂的几何形状、材料特性和边界条件。 由于已有成熟的大型软件系统支持，有限元发已成为一种非常受欢迎的、应用极广的数值计算方法。 有限元法的基本求解步骤 ① 将连续体离散成有限个单元的组合体。 ② 利用单元结点的场量作为基本未知量，选择一组插值函数，确定单元内相应的场量分布。 ③ 建立单元内结点的基本未知量与载荷间的平衡方程。 ④ 将单元内结点基本未知量的平衡方程集成，得出以结点基本未知量的平衡方程组。 ⑤ 求解代数方程组，得出各结点的基本未知量。 ① 20世纪50年代，发展与萌生 单一功能程序，简单单元 1954 ~ 1955年，J. H. Argyris在航空工程杂志上发表了一组能量原理和结构分析论文。 1956年，M. J. Turner, R. W. Clough, H. C. Martin, L. J. Topp在纽约举行的航空学会年会上介绍了一种新的计算方法，将矩阵位移法推广到求解平面应力问题。他们把结构划分成一个个三角形和矩形的“单元”，利用单元中近似位移函数，求得单元节点力与节点位移关系的单元刚度矩阵。 ② 20世纪60年代，数学基础与证明 单一功能程序，多种单元 1960年，Clough在论文“The finite element in plane stress analysis”中首次提出有限元 (finite element) 这一术语。 1963年前后，J. F. Besseling , R. J. Melosh, R. E. Jones 等人认识到有限元法就是变分原理中 Ritz 近似法的一种变形，发展了用各种不同变分原理导出的有限元计算公式。 1965年，O. C. Zienkiewicz 和Y. K. Cheung (张佑启) 发现只要能写成变分形式的所有场问题，都可以用与固体力学有限元法的相同步骤求解。 1969年，B. A. Szabo 和G. C. Lee 指出可以用加权余量法特别是 Galerkin 法，导出标准的有限元过程来求解非结构问题。 ③ 20世纪70年代，单元库丰富 线性到非线性通用程序，如SAP，NONSAP等 ④ 20世纪80年代，功能扩大 大型通用程序如ADINA等 ⑤ 20世纪90年代，应用领域扩大 大型商用软件，如ANSYS、MARC、NASTRAN等 ⑥ 目前，面向工程，与CAD结合 从20世纪60年代中期以来，大量的理论研究不但拓展了有限元法的应用领域，还开发了许多通用或专用的有限元分析软件，为工程设计和优化提供了有力的工具。 由弹性力学平面问题 扩展到空间问题、板壳间题 由静力平衡问题 扩展到动力问题、稳定问题 由线性分析