数值模拟的理论与方法

第二节 数值分析方法（4/6）- 有限元法/FEM
有限元法又可分为位移法、利用余位进行变化的方法和用混合积分的混 合法三种。 有限元法的位移法，其实质就是将求解区域划分为有限个单元，通过 构造插值函数，把问题化为一个变分问题（即求泛函数值的问题），经过 离散化得到计算格式，利用计算格式来求解相应问题。变分法证明求解某 些微分方程的问题等效于将泛函数的相关量进行最小化。如果相关于因变 量的节点值使泛函数最小，那么所得到的条件表达式就是所需要的离散化 方程。也就是说，求解一个微分方程边值问题就可以通过寻找某一变分问 题的极值函数来解决。有限元解题的基本过程: 对一个具体的工程应用分析， 在确定了分析计算的基本方案后，就可以按建模（即建立几何模型）、分 网（即建立有限元模型）、加载（即给定边界条件）、求解（有限元求解） 和后处理（即计算结果的可视化）等几个步骤实施分析计算。
新进展-热应力模拟及热裂纹预测
(华铸课题组成果) -热应力模拟及热裂纹预测
FEM热应力模拟
行星架铸件
裂纹
热裂倾 向较大
FEM热应力模拟 (华铸课题组成果)
数控机床横梁铸件长约11米，重约 30吨，最薄壁厚只有30毫米。
铸造成形模拟仿ቤተ መጻሕፍቲ ባይዱ技术
FEM热应力模拟 (华铸课题组成果)
横梁铸件等效应变分布及变形情况（变形放大 10倍）
《有限差分法基础》讲义
第2章 数值模拟方法概述
廖敦明 华中科技大学
18071121688, 87558134 liaodunming@
华中科技大学材料学院华铸软件中心 材料成形与模具技术国家重点实验室
第一节 研究目的与研究内容（1/13） 1．研究目的 数值模拟(CAE)技术是通过建立能够准确描述研究 对象某一过程的数学模型，采用合适可行的求解方法， 使得在计算机上模拟仿真出研究对象的特定过程，分析 有关影响因素，预测这一特定过程的可能趋势与结果。 材料成形数值模拟CAE技术最终的研究目的是在计 算机虚拟的环境下，通过交互方式，能够制定合理的工 艺，而不需要或少做现场试生产。从而可以大幅度缩短 新产品开发周期，降低废品率，提高经济效益。

流体流动数值模拟流体流动现象普遍存在于⾃然界及多种⼯程领域中。

所有这些流动过程都遵循质量守恒、动量守恒、能量守恒和组分守恒等基本物理定律；⽽且流动若处于湍流状态，则该流动系统还要遵守附加的湍流输运⽅程。

本讲座将依据流体运动的特性阐述计算流体动⼒学的相关基础知识及任务；在流体运动所遵循的守恒定律及其数学描述的基础上，介绍数值求解这些基本⽅程的思想及其求解过程。

第⼀节计算流体动⼒学概述计算流体动⼒学（CFD）技术⽤于流体机械部流动分析及其性能预测，具有成本低，效率⾼，⽅便、快捷⽤时少等优点。

近年来随着计算流体⼒学和计算流体动⼒学及计算机技术的发展, CFD技术已成为解决各种流体运动和传热，以及场问题的强有⼒、有效的⼯具，⼴泛应⽤于⽔利、⽔电，航运，海洋，冶⾦，化⼯，建筑，环境，航空航天及流体机械与流体⼯程等科学领域。

利⽤数值计算模拟的⽅法对流体机械的部流动进⾏全三维整机流场模拟，进⽽进⾏性能预测的⽅法越来越⼴泛地被从事流体机械及产品性能取决于各种场特性的设计、科研等科技⼈员所使⽤；过去只有通过实验才能获得的某些结果或结论，现在完全可借助CFD模拟的⼿段来准确地获取。

这不仅既可以节省实验资源,还可以显⽰从实验中不能得到的许多场特性的细节信息。

⼀、什么是计算流体动⼒学计算流体动⼒学(Computational Fluid Dynamics，简称CFD)是通过计算机数值计算和图像显⽰，对包含流体流动和有热传导等相关物理现象的系统所做的分析。

CFD的基本思想可以归结为：把原来在时间域及空间域上连续的物理场（如速度场和压⼒场，以及热⼒场等），⽤⼀系列有限个离散点上变量值的集合来代替；并通过⼀定的原则和规律建⽴起关于这些离散点上的场变量之间关系，从⽽组成这些场变量之间关系的代数⽅程组；然后求解这种代数⽅程组，来获得这些场变量的近似值[1-3]；这就是流动的数值计算。

或者直观地说，通过数值计算中的各种离散⽅法，把描述连续流体运动的控制偏微分⽅程离散成代数⽅程组，由此建⽴该流动的数值模型；再根据问题的具体情况，设定边界条件和初始条件封闭⽅程组；然后通过计算机数值计算求解这种代数⽅程组，从⽽获得描述该流场场变量的某些运动参数的数值解。

数值模拟偏微分方程的三种方法：FDM、FEM及FVM偏微分方程数值模拟常用的方法主要有三种：有限差分方法（FDM）、有限元方法（FEM）、有限体积方法（FVM），本文将对这三种方法进行简单的介绍和比较。

有限差分方法有限差分方法(Finite Difference Methods)是数值模拟偏微分方程最早采用的方法，至今仍被广泛运用。

该方法包括区域剖分和差商代替导数两个过程。

具体地，首先将求解区域划分为差分网格，用有限个网格节点代替连续的求解区域。

其次，利用Taylor级数展开等方法将偏微分方程中的导数项在网格节点上用函数值的差商代替来进行离散，从而建立以网格节点上的值为未知量的代数方程组。

该方法是一种直接将微分问题变为代数问题的近似数值解法，数学概念直观，表达简单，是发展较早且比较成熟的数值方法。

差商代替导数后的格式称为有限差分格式，从格式的精度来考虑，有一阶格式、二阶格式和高阶格式。

从差分的空间离散形式来考虑，有中心格式和迎风格式。

对于瞬态方程，考虑时间方向的离散，有显格式、隐格式、交替显隐格式等。

目前常见的差分格式，主要是以上几种格式的组合，不同的组合构成不同的差分格式。

差分方法主要适用于结构网格，网格的步长一般根据问题模型和Courant稳定条件来决定。

请输入标题有限元方法(Finite Element Methods)的基础是变分原理和分片多项式插值。

该方法的构造过程包括以下三个步骤。

首先，利用变分原理得到偏微分方程的弱形式(利用泛函分析的知识将求解空间扩大)。

其次，将计算区域划分为有限个互不重叠的单元(三角形、四边形、四面体、六面体等)。

再次，在每个单元内选择合适的节点作为求解函数的插值点，将偏微分方程中的变量改写成由各变量或其导数的节点值与所选用的分片插值基函数组成的线性表达式，得到微分方程的离散形式。

利用插值函数的局部支集性质及数值积分可以得到未知量的代数方程组。

有限元方法有较完善的理论基础，具有求解区域灵活(复杂区域)、单元类型灵活(适于结构网格和非结构网格)、程序代码通用(数值模拟软件多数基于有限元方法)等特点。

流体力学数值模拟报告一、背景随着计算机技术的快速发展，数值模拟方法成为流体力学研究和工程设计中的重要工具。

传统的实验室试验和理论计算往往只能针对特定的情况进行研究，而数值模拟可以模拟各种复杂的流体问题，研究范围更广，具有更大的灵活性。

因此，流体力学数值模拟在科学研究、工程设计和产品优化中得到了广泛应用。

二、方法流体力学数值模拟通常采用计算流体力学（CFD）方法，通过将流体问题转化为求解计算网格中的运动方程来模拟流体运动。

其中，最常用的方法是有限体积法和有限元法。

有限体积法将流体区域分割为有限个体积单元，在每个单元上求解运动方程，得到整个流场的数值解。

有限元法则是通过在每个单元内定义一个函数，并在这个函数上求解运动方程，最后将每个单元上的解拼接在一起得到整个流场的数值解。

这些方法通常需要进行网格生成、边界条件处理和数值计算等工作。

三、应用流体力学数值模拟广泛应用于航空航天、汽车工程、建筑设计等领域。

例如，在航空航天领域，数值模拟可以预测飞机在不同速度下的空气动力学性能，优化机翼形状和发动机设计。

在汽车工程中，数值模拟可以模拟车辆行驶时的气动特性和燃烧过程，优化车身外形和发动机参数。

在建筑设计中，数值模拟可以预测建筑结构受风荷载时的动态响应，指导建筑物的抗风设计。

四、重要性流体力学数值模拟在工程领域具有重要的意义。

首先，它可以帮助工程师更好地理解和预测流体的行为，指导产品设计和性能优化。

其次，数值模拟可以大大减少试验成本和周期，提高工程设计的效率。

最后，数值模拟可以模拟各种复杂情况，例如高速气流、湍流和多相流等，为工程师提供更全面、准确的数据，帮助他们做出更好的决策。

五、未来发展随着计算机技术的进一步发展和数值模拟方法的日益成熟，流体力学数值模拟将在未来得到更广泛的应用。

首先，随着高性能计算机的出现，数值模拟可以处理更大规模的流体问题，提供更准确的结果。

其次，数据驱动的流体力学模拟方法将成为一个新的研究方向，通过机器学习和数据挖掘技术，将实验数据和模拟数据相结合，提高模拟结果的准确性和鲁棒性。