■BEM プログラム: BEM8CON.FOR■
境界要素法で最もシンプルなプログラムを紹介します。プログラム名は、BEM8CON.FOR で、離散化にconstant element を用いています。つまり、今まで学習してきたことが、このプログラムに集約されています。

まず、このプログラムの概要を紹介します。下図を見て下さい。

上図を簡単に説明しますと、次の様になります。

1. まず、要素の座標情報と境界条件を読み込みます。
2. 次に、マトリックス[C],[F],[G] を作成します。
3. マトリックスと境界条件で、連立方程式 [A]{x}={RHS} の形にして、{x}について解きます。
4. 求められた{x}には、{h}と{q_n} が含まれますので、分離する必要があります。
5. 境界上の{h}と{q_n}が計算されたので、領域内のh(ξ) が計算できます。
6. 結果をファイルにする。