结构力学数值方法：有限体积法(FVM)：二维平面应力和平面应变问题.pdf

结构力学数值方法：有限体积法(FVM)：二维平面应力和平

面应变问题

1绪论

1.1有限体积法的起源和发展

有限体积法（FiniteVolumeMethod,FVM）起源于20世纪50年代，最初被

应用于流体力学领域，特别是解决气体动力学中的偏微分方程。FVM的核心思

想是基于守恒定律，将连续的物理域离散成一系列控制体积，然后在每个控制

体积上应用守恒定律，从而将偏微分方程转化为代数方程组。这种方法在处理

对流主导问题、非线性问题以及复杂几何形状时表现出色，因此迅速在工程计

算中得到广泛应用。

随着计算机技术的发展，有限体积法逐渐扩展到其他领域，如热传导、电

磁学、结构力学等。在结构力学中，FVM被用于解决二维平面应力和平面应变

问题，其优势在于能够直接处理守恒形式的方程，提供更准确的通量计算，以

及在处理不规则网格时的灵活性。

1.2FVM与FEM的比较

有限体积法与有限元法（FiniteElementMethod,FEM）是两种广泛应用于

工程计算的数值方法。尽管它们在某些方面有相似之处，如都是将连续问题离

散化，但它们的基本原理和应用领域存在显著差异。

1.2.1原理差异

有限体积法：基于守恒定律，将计算域划分为一系列控制体积，

然后在每个控制体积上应用守恒定律，得到代数方程组。FVM强调的是

通量的守恒，适用于对流主导问题。

有限元法：基于变分原理和加权残值法，将计算域划分为一系列

单元，然后在每个单元上建立位移或应力的近似函数，通过求解最小势

能原理或加权残值方程得到未知量。FEM适用于弹性力学、固体力学等

位移场问题。

1.2.2应用领域

有限体积法：在流体力学、热传导、电磁学等领域有广泛应用，

特别是在处理对流问题时，FVM能够提供更准确的通量计算。

有限元法：在结构力学、固体力学、断裂力学等领域是首选方法，

FEM能够处理复杂的边界条件和材料非线性问题。

1

1.2.3代码示例

下面是一个使用Python实现的简单有限体积法示例，用于解决一维稳态热

传导问题。虽然这里讨论的是二维平面应力和平面应变问题，但一维热传导问

题的代码可以帮助理解FVM的基本实现步骤。

importnumpyasnp

#定义网格参数

nx=10#网格点数

dx=1.0/(nx-1)#网格间距

k=1.0#热导率

#初始化温度场

T=np.zeros(nx)

#定义边界条件

T[0]=100.0#左边界温度

T[-1]=0.0#右边界温度

#内部节点的离散方程

foriinrange(1,nx-1):

T[i]=(T[i-1]+T[i+1])/2.0-k*dx*(T[i+1]-T[i-1])/(2.0*dx)

#打印温度分布

print(T)

1.2.4解释

在这个示例中，我们使用有限体积法来求解一维稳态热传导问题。首先，

我们定义了网格参数，包括网格点数和网格间距。然后，初始化温度场，并设

定边界条件。接下来，我们应用有限体积法的基本原理，即在每个控制体积上

应用热传导守恒定律，通过迭代求解内部节点的温度。最后，打印出整个网格

的温度分布。

尽管这个例子非常简化，但它展示了有限体积法的基本思想：通过在每个

控制体积上应用守恒定律，将连续问题转化为一系列代数方程，从而求解未知

量。在处理二维平面应力和平面应变问题时，FVM的实现会更加复杂，需要考

虑应力和应变的守恒，以及材料的弹性性质。

2有限体积法基础

2.1控制体积的概念

在结构力学的数值分析中，有限体积法（FVM）是一种广泛使用的方法，

2

它基于守恒定律，通过在空间上划分控制体积来求解偏微分方程。控制体积，

也称为控制单元，是FVM中的基本概念，它指的是空间中任意形状的封闭区域，

用于积分守恒方程。在二维平面应力和平面应变问题中，控制体积通常被定义

为网格中的一个单元，如矩形或三角形。

2.1.1控制体积的定义

控制体积可以是任意形状，但在