BEM — CAE用語解説

根本的な違いは、支配方程式を満たす「基本解」を使う点です。FEMは領域全体を近似関数で埋め尽くしますが、BEMは積分方程式に変形し、境界上の未知量だけを解きます。例えば、3次元のラプラス問題では、支配方程式

良い指摘です。基本解は線形で「定数係数」の微分方程式に対してのみ、解析的に導出できます。つまり、均質な線形弾性体、定常熱伝導、アコースティックス等が典型です。材料定数が場所によって変わる（不均質）場合や、塑性変形のような材料非線形問題には、標準的なBEMは適用が難しい。その場合は領域分割法（Domain Decomposition）で均質領域に分けたり、体積積分項を導入するDR-BEM（Dual Reciprocity BEM）などの拡張手法が必要になります。

その通り、これがBEMの大きな利点の一つです。基本解自体が無限遠方でゼロに収束する性質を持っているため、境界積分方程式を立てる際、無限遠方の境界は条件を明示的に考慮する必要がありません。例えば、飛行機の風洞実験をCFDでシミュレートする場合、FEM/CFDでは計算領域の外側を「遠方境界」としてメッシュ切り、そこで流速を一様流などと設定します。BEM（あるいはパネル法）では、機体表面だけをパネル分割し、基本解の重ね合わせで無限遠方の一様流条件を自動的に満たすことができます。これにより、メッシュ作成が劇的に簡素化されます。

要素とシェープ関数の考え方はFEMに似ていますが、未知数が異なります。BEMでは、境界を「パネル」や「要素」に分割し、各要素上で物理量（ポテンシャルuとフラックス∂u/∂n）を、節点値を用いたシェープ関数（線形や二次）で近似します。問題は特異積分です。被積分関数の

FEMの剛性行列がスパース（帯行列）で対称正定値であるのに対し、BEMから得られる係数行列は「密」で、一般に非対称です。なぜなら、基本解

BEMが威力を発揮するのは、「無限or半無限領域」「薄肉構造」「応力集中部の詳細解析」です。具体例を挙げると：

あります。FEMではディリクレ条件（固定）とノイマン条件（荷重）を節点や要素面に直接割り当てますが、BEMでは「基本解の種類」が境界条件と密接に関わります。例えば、定常熱伝導で「温度uが既知」のディリクレ境界と、「熱流束qが既知」のノイマン境界が混在する場合、離散化方程式

AnsysやAbaqusといった汎用構造・流体ソフトのメインソルバーはFEM/FVMです。しかし、特定の物理分野向けにBEMモジュールを有するものや、BEM専用ソフトがあります。

現実的な適用例です。例えば「自動車室内の騒音解析」を考えましょう。音源であるエンジンや路面からの振動は、複雑な構造を伝わって来ます。この構造振動の伝達は、材料非線形を無視すればFEMで精度良く計算できます。しかし、振動するボディから放射される「音」は無限空間に伝わるため、BEMが適しています。この場合のワークフローは：

これは典型的な問題で、主に二つの原因が考えられます。

数値積分、特に「特異積分」の精度不足が第一の疑いです。先ほど話した特異点処理が不十分だと、その周りの積分値に誤差が集中します。対策としては：

BEMのボトルネックは、まさにその「密行列」の扱いです。対策は以下の通りです。