■Program BUCKLE1.FORの紹介■

BUCKLE1.FORは、BUCKLE.FORよりかなり賢くなっています。BUCKLE.FORでは、作成されるマトリックスはフル(NxN)でしたが、BUCKLE1.FORでは、Nx3 になっています。下図をご覧下さい。

左の Matrixを 右図の様に 変形する。 =======＞

つまり、連立方程式の解法で学習したTri-Diagonal Matrix を使っているからです。
Tri-Diagonal Matrixを使えるようにするために、Subroutine MATRIX, FORM, SYSTEM が改良されていますので、それらのSubroutinesを注意深く観察して下さい。

■Program BUCKLE1A.FORの紹介■
このプログラムは、もっと賢くなっています。連立方程式のマトリックスが対称だと、下図に示す様に Tri-Diagonal Matrixを更にコンパクトにすることが出来ます。

対称Matrix だと 第1列が不要 =======＞

つまり、対称Tri-Diagonal Matrix の活用です。詳しくは、連立方程式の解法をご覧下さい。
物事は、どこかが簡単になると、別の場所が複雑になります。Subroutine MATRIX, FORM, SYSTEM がかなり複雑になっています。この時点で、ロジックをしっかりと理解しておいて下さい。2次元微分方程式のFEMのプログラムでは、もっと複雑になりますからね。

■ExcelによるBUCKLE1の紹介■
ExcelをベースにしたBACKLEのプログラムを紹介します。これは、原田様より提供して頂きました。原田、有難うございます。セルに例題で紹介したデータを入力しGOボタンをクリックするとDISPLACEMENTが計算されます。そして、 GRAPHボタンをクリックすると、グラフを描きます。皆様も使ってみて下さい。

問題

疑問および質問