■ BUCKLE2.FORとBUCKLE2A.FORの紹介■
1次元Helmholtz Equation を、2次要素で解析するプログラムを2つ紹介します。上のプログラム名をクリックしてみて下さい。BUCKLE2.FORは、Full banded Matrixの[K]を生成し、そしてBUCKLE2A.FORは Half banded Matrix の[K]を作ります。BUCKLE2.FORは、FEMの基本的なロジックを、勉強するのに有効です。

プログラムの構成は、Gauss積分の部分を除くと、1次要素のプログラムとほぼ同じです。また、データ作成には、前に紹介したSET1UP.FORが使えます。

ここで紹介しているプログラムで、新たに追加されているSubroutine名とその機能を下表に示します。

Subroutine名	機　能
DERIV	Gauss積分法のサンプリング点でのdNi/dξを計算します。つまり、形状関数の微分の計算です。
SHAPEF	Gauss積分法のサンプリング点での形状関数を計算します。
ISOPARA	引数ξに対し、3つの形状関数の値を返します。
SGSM	要素毎に、-[B]T[B] + α2[N]T[N] の積分を計算します。

Gauss積分法の関連の変数としてINTEPTがあります。これは、サンプリング点数を意味します。

下に、プログラムの流れ図を示します。参考にして下さい。何か不明な点があれば、私まで E-Mail下さい。

これで、1次元のFEMは終了です。次は2次元です。引き続き頑張って下さい。私も時間の許す限り貴方の勉強の御手伝いをします。それでは、2次元で会いましょう。

START
Gauss積分法のサンプリング点で N1(ξ)、 N2(ξ)、 N3(ξ)、dN1/dξ、dN2/dξ、dN3/dξを計算し、配列 SF(I,J)とBP(I,J)へ入れる。
データを読む。(データは1次要素とほぼ同じ。)
MATRIX: SGSMで計算された要素毎のSTIFF(I,J)をGlobal マトリックスA(I,J)ヘAssembleする。
FORM: 境界条件をA(I,J)に組み込む。
SYSTEMA: 連立方程式を解く。
結果を出力する。
STOP

問題

質問