关于有限基本解法介绍

[拼音]：youxian jiben jiefa

[外文]：method of finite fundamental solutions

解线性势流动的一种数值计算方法。它用一些形式比较简单、而在流动区域内又满足方程的解析函数（如位势流的源、汇、偶极子以及涡旋等）作为基本解，再将它们线性叠加，以满足任意外形物体的边界条件，从而模拟出各种具体流动的速度场。

以位势流动为例，格林定理和斯托克斯定理指出：扰动速度υ(P)（P为流动场中的任一点)可用流场边界上源、汇或偶极子的分布来表示，而扰动速度场则线性依赖于流场边界的源、汇或偶极子的分布密度。因此扰动速度可以用物体表面的源、汇分布密度求得。在一般情况下，可将物体表面分成许多连接的单元，如果单元尺度比流场特征尺度小，可以假定单元上的源、汇或偶极子的密度分布是均匀的。这时空间任意一点P上的扰动速度υ(P)可写成：

式中ej(qi)为第j个单元上分布密度为1的源、汇或偶极子在P点所诱导的速度;σj为该单元的分布密度。如果物面上的单元总数为N，则上式中只有N个待定系数，这些系数可以利用物面上N个点处的边界条件来确定，这N个条件可写成：

式中A嗎=n(qi)·ej(qi)；Bi=-n(qi)·υ∞；n(qi)是物理面上qi点处单位法向矢量，它指向流场内部；qi为控制点。从上述方程组中解出σj后，即可算得扰动速度场。

用源、汇或偶极子来求解十分方便，但这类基本解都有奇点，这些奇点可以是孤立的，也可以是分布在某些曲面或曲线上的。在这些地方必须作一些特殊处理。

在实际计算时，单元的分法，单元上的密度分布形式和控制点的位置，都会直接影响到计算的准确性。如果控制点选得不当，会得到不准确甚至是荒谬的结果。目前还没有确定控制点正确位置的严格理论。计算表明，对等密度分布的单元来说，把控制点选在单元形心或单元自身诱导速度最小点处，可得到比较满意的结果。在单元上，如采用多参数的密度分布形式，则用较少的单元块数也可以得到同样精度的结果。

有限基本解法多用于位势绕流问题，在工程上已能成功地计算或校核复杂形状物体上的气动载荷，甚至可直接用来设计飞行器等的外形。这一方法近来已进一步用于研究可压缩情况下的有限扰动问题。此外，在水工结构的载荷和油田开采等计算中也有应用。

参考书目

1. J.L.Hess， Computer Method， Applied Mechanics and Engineering， p. 145， March 1975.