要素のMatrixです。
STIFF(I,J)が、Global Matrix A(I,J) にAssemble されるところを、注意深く観察して下さい。NODEX(I,J)が重要な役目を果たしていますね。
RLは要素の長さです。
Subroutine FORM
ここでは、境界条件がGlobal Matrix A(I,J) に組み込まれます。IBTYPE(1)は上流(x=0.)の、そしてIBTYPE(2)は下流(x=L)の境界条件が入っています。IBTYPE(i)=1 の場合はDirichlet、IBTYPE(i)=2 は Neumann型となっています。BV(i)には、境界値が入ります。
Subroutine SYSTEM
MATRIXとFORMで形成された連立方程式が、ここで解析されます。上流と下流での境界条件がDirichletの場合、解析される連立方程式は下図の様になっているはずです。
連立方程式の計算方法については、ここへジャンプして下さい。
Program SET1UP.FORは、貴方に代わってデータを作成してくれます。貴方はこのプログラム実行し、質問に応えるだけで、BUCKLE.DATが出来上がります。コンソールの入出力は、DEC Fortran に準拠しています。このデータ作成プログラムは、次の紹介するBUCKLE1.FOR、BUCKLE1A.FOR、BUCKLE2A.FORのデータも作成します。
