结构力学数值方法：有限体积法(FVM)：结构力学基础理论.pdf

结构力学数值方法：有限体积法(FVM)：结构力学基础理论

1绪论

1.1有限体积法的历史与发展

有限体积法（FiniteVolumeMethod,FVM）作为一种数值求解偏微分方程的

方法，其历史可以追溯到20世纪50年代。起初，FVM被广泛应用于流体力学

领域，特别是解决连续介质的守恒定律问题。随着计算机技术的发展和数值方

法理论的完善，FVM逐渐扩展到其他工程领域，包括结构力学。FVM的核心思

想是基于控制体的概念，将计算域划分为一系列控制体，然后在每个控制体上

应用守恒定律，从而将连续的偏微分方程转化为离散的代数方程组。

1.1.1发展历程

20世纪50年代至60年代：FVM最初由气象学家和流体力学家提

出，用于大气和流体流动的数值模拟。

20世纪70年代：随着计算机的普及，FVM开始在工业界得到应

用，特别是在航空航天和汽车工业中。

20世纪80年代至90年代：FVM的理论基础得到进一步发展，开

始应用于更广泛的工程问题，包括热传导、电磁学和结构力学。

21世纪初至今：FVM的算法和应用不断优化，与有限元法

（FEM）、有限差分法（FDM）等其他数值方法结合，形成更强大的数值

模拟工具。

1.2FVM在结构力学中的应用领域

在结构力学中，FVM主要用于解决结构的静力学、动力学和热力学问题。

与有限元法相比，FVM在处理对流主导问题和非线性问题时具有一定的优势，

因为它更直接地基于守恒定律，能够更好地保持物理量的守恒性。

1.2.1静力学应用

在结构静力学中，FVM可以用于求解结构的应力和应变分布。例如，对于

弹性力学问题，可以将结构划分为多个控制体，然后在每个控制体上应用平衡

方程和本构关系，从而得到结构的响应。

1.2.2动力学应用

在结构动力学中，FVM可以用于模拟结构的振动和冲击响应。通过在控制

体上应用动力学方程，可以求解结构在动态载荷作用下的位移、速度和加速度。

1

1.2.3热力学应用

在结构热力学中，FVM可以用于求解结构的温度分布和热流。对于热传导

问题，FVM通过在控制体上应用能量守恒方程，可以得到结构的温度场和热流

分布。

1.2.4示例：使用Python实现FVM求解一维弹性杆的静力学问题

importnumpyasnp

#定义材料属性和几何参数

E=200e9#弹性模量，单位：Pa

A=0.01#截面积，单位：m^2

L=1.0#杆长，单位：m

n=10#控制体数量

dx=L/n#控制体长度

#定义载荷和边界条件

P=1000#外力，单位：N

u_left=0#左端位移，单位：m

u_right=0#右端位移，单位：m

#初始化位移和应力数组

u=np.zeros(n+1)

sigma=np.zeros(n)

#应用FVM求解

foriinrange(n):

#计算控制体的应力

sigma[i]=P/A

#计算控制体的位移

ifi==0:

u[i+1]=u_left+sigma[i]*dx/E

else:

u[i+1]=u[i]+sigma[i]*dx/E

#输出结果

print(位移分布：,u)

print(应力分布：,sigma)

1.2.5示例解释

上述代码示例展示了如何使用有限体积法求解一维弹性杆在轴向载荷作用

下的静力学问题。首先，定义了材料属性（弹性模量E和截面积A）、几何参数

2

（杆长L）和计算参数（控制体数量n）。然后，初始化了位移和应力数组，并

应用了边界条件。在主循环中，计算了每个控制体的应力和位移，最后输出了

位移和应力的分布结果。

通过这个简单的示例，我们可以看到FVM在结构力学问题中的基本应用流

程，即定义问题、初始化参数、应用守恒定律和输出结果。在实际应用中，

FVM可以处理更复杂的问题，如三维结构、非线性材料和多物理场耦合问题。

2有限体积法基础