■9-Noded要素のFEM9Q.FORの紹介■
FEM9Q.FORでは、Main ProgramのPARAMTERの中のNDがND=9になっていることに注意して下さい。ND=9ですから、Subroutine ISOPARAに9つの形状関数が存在することになります。
他のプログラム同様、FEM9Q.FORにも、入力データを作成するSETFEM9.FORが用意されています。領域と境界条件は、以前と同じです。このプログラムを実行し¨EX=2とNEY=2を入力すると、 FEM9DATA.QQQ を作成します。要素構成は、下図の様になっています。そして、節点数は25です。

上の図には、要素番号が割り当ててありません。スペースの都合で書けなかっただけです。入力データと見比べて要素番号を確認して下さい。

入力データのエコーと結果は、他のプログラム同様、 SOLUTION.FEM へ吐き出されます。

これまでに、同じ問題を4つのプログラム( FEM3Q.FOR, FEM4Q.FOR, FEM8Q.FOR, FEM9Q.FOR )で解析してきました。このままでは、勉強になりませんので、問題を解いて、解析精度や解の安定性について、学習して下さい。