ANSYS 2025 R2多物理场仿真从入门到工程实战

第1章 仿真与ANSYS概览

1.1 有限元法核心逻辑与工程应用

1.1.1 有限元法的基本思想

有限元法（Finite Element Method，FEM）是一种将连续体离散化为有限个单元，通过在单元内构造近似函数来求解偏微分方程边值问题的数值计算方法。其核心思想可以概括为“化整为零，积零为整”——将复杂的连续结构或场域分割成大量简单的小单元，在单元层面上建立近似方程，再通过单元间的连接条件组装成整体方程进行求解。

从数学角度看，有限元法本质上是一种求解偏微分方程的数值近似方法。当工程问题中的几何形状、边界条件或材料属性过于复杂而无法获得解析解时，有限元法提供了一条行之有效的数值求解路径。工程中的绝大多数实际问题——无论是桥梁在风载下的变形、发动机缸体中的温度分布，还是变压器周围的磁场分布——都需要借助有限元法来获得可靠的数值解答。

有限元法的数学基础可以追溯到变分原理和加权余量法。变分原理指出，对于许多物理问题，求解偏微分方程等价于求解某个泛函的极小值问题。加权余量法则通过要求微分方程在某种平均意义下的余量为零来构造近似解。这两种数学工具为有限元法的收敛性和误差估计提供了理论基础。

1.1.2 有限元法的基本求解步骤

一个完整的有限元分析通常包含以下六个基本步骤：

第一步：结构离散化。 将连续的分析对象划分为有限个、互不重叠的单元，单元之间通过节点相互连接。这一步骤称为网格划分。网格的质量和密度直接决定了计算结果的精度和计算效率。单元类型的选择（如杆单元、梁单元、平面单元、实体单元、壳单元等）取决于分析对象的几何特征和受力状态。

第二步：选择位移模式。 在每个单元内部，假设位移（或其他场变量）按照某种简单函数分布，这种函数称为位移函数或形函数。形函数通常采用多项式形式，多项式的阶次决定了单元的精度等级。线性单元在每个单元内位移呈线性分布，二次单元则呈二次分布，后者对弯曲变形的描述能力更强。

第三步：推导单元刚度矩阵。 利用弹性力学的基本方程——几何方程（应变-位移关系）、物理方程（应力-应变关系，即本构关系）和平衡方程，结合虚功原理或最小势能原理，建立节点力与节点位移之间的关系式，即单元刚度方程：\[F^e = K^e \cdot u^e\]，其中K^e为单元刚度矩阵。

第四步：组装整体刚度矩阵。 将所有单元的刚度矩阵按照节点编号组装成整体结构的刚度方程：\[F = K \cdot u\]。这一过程基于节点处的位移连续条件和力平衡条件。整体刚度矩阵通常是一个大型稀疏矩阵，具有带状分布的特点。

第五步：引入边界条件。 刚度方程本身是奇异的，需要引入足够的位移边界条件（约束条件）才能得到唯一解。边界条件的施加方式包括直接代入法、乘大数法和拉格朗日乘子法等。

第六步：求解代数方程组。 引入边界条件后，整体刚度方程成为可解的线性代数方程组，可采用直接法（如高斯消元法、LU分解）或迭代法（如共轭梯度法、GMRES）求解节点位移。

第七步：计算导出量并后处理。 由节点位移通过几何方程计算单元应变，再通过物理方程计算单元应力。对于瞬态分析或非线性分析，还需要考虑时间积分和非线性迭代求解。

1.1.3 有限元法的收敛准则

有限元法的收敛性是指当网格不断加密（单元尺寸趋于零）时，数值解趋近于精确解的性质。保证收敛性的两个基本准则是：

· 完备性要求： 单元的位移函数必须包含刚体移模式和常应变模式。这意味着单元能够描述刚体平动、刚体转动以及均匀拉伸等基本变形状态。如果位移函数不满足完备性，无论网格如何加密，数值解都不会收敛到精确解。

· 协调性要求： 相邻单元在公共边界上的位移必须是连续的。满足协调性的单元称为协调元，不满足的称为非协调元。非协调元在某些情况下（如板壳分析中的Mindlin板单元）如果能够通过分片检验，也能获得收敛的结果。

1.1.4 结构分析中的应用

有限元法在结构分析中主要解决以下四类问题：

线弹性静力学分析是有限元法最基础、最成熟的应用领域。其控制方程为KU=F，K为刚度矩阵，U为位移向量，F为载荷向量。典型的工程应用包括机械零件在额定载荷下的强度校核、建筑结构在自重和使用荷载下的变形分析、压力容器在内压作用下的应力分布等。

模态分析用于确定结构的固有频率和振型。控制方程为广义特征值问题：(K-ω²M)φ=0，其中M为质量矩阵，ω为固有频率，φ为振型向量。模态分析是动力学分析的基础，广泛应用于旋转机械的临界转速计算、航空航天结构的颤振分析、建筑结构的抗震设计等。

屈曲分析用于预测结构在压缩载荷下发生失稳的临界载荷。线性屈曲分析归结为广义特征值问题：(K+λKσ)φ=0，其中Kσ为应力刚度矩阵，λ为载荷因子。非线性屈曲分析则需要追踪载荷-位移曲线的极值点或分叉点。

非线性分析考虑材料非线性（塑性、超弹性、粘弹性）、几何非线性（大变形、大应变）和接触非线性（边界条件的非线性变化）。非线性问题的求解通常采用增量-迭代方法，如Newton-Raphson法。

1.1.5 流体分析中的应用

计算流体动力学（Computational Fluid Dynamics，CFD）是有限元法（以及有限体积法）的另一重要应用领域。流体流动的控制方程为Navier-Stokes方程：

连续性方程：∂ρ/∂t + ∇·(ρv) = 0

动量方程：ρ(∂v/∂t + v·∇v) = -∇p + ∇·τ + ρg

能量方程：ρc_p(∂T/∂t + v·∇T) = ∇·(k∇T) + Φ

其中ρ为密度，v为速度向量，p为压力，τ为粘性应力张量，g为重力加速度，c_p为比热容，k为热导率，Φ为粘性耗散函数。

CFD分析的典型工程应用包括：飞行器气动外形的升阻力计算、汽车车身外流场分析与减阻优化、管道和阀门的内流场分析、换热器的流动与换热分析、搅拌釜内的混合过程模拟、室内通风与空调气流组织分析等。

1.1.6 热分析中的应用

热传递有三种基本机制：

· 热传导： 遵循傅里叶定律，q = -k∇T。热量从高温区域向低温区域传递，传递速率与温度梯度和材料热导率成正比。

· 热对流： 遵循牛顿冷却定律，q = h(T_s - T_∞)。热量通过流体运动在固体表面和流体之间交换，换热系数h受流速、流体物性和表面几何形状影响。

· 热辐射： 遵循斯特藩-玻尔兹曼定律，q = εσ(T_s⁴ - T_surr⁴)。热量通过电磁波传递，不依赖介质，在高温条件下辐射传热占主导地位。

典型的热分析工程应用包括：电子元器件的散热设计与热管理、发动机缸体和排气管的热应力分析、建筑围护结构的保温隔热分析、铸造过程中的凝固模拟、焊接过程中的温度场和残余应力分析、电池包的热管理设计等。