■計算精度の検証: その1■
これまでに、同じ問題を5つのプログラム( FEM3Q.FOR, FEM4Q.FOR, FEM8Q.FOR, FEM9Q.FOR, FEM12Q.FOR )で解析してきました。ここで計算精度がどうなっているのかチェックしておきましょう。精度検証の位置は、(x,y)=(1,1)と(x,y)=(2,1)の2ヶ所で行います。要素分割は、(NEX,NEY)=(2,2)から(NEX,NEY)=(20,20)の範囲内としておきます。計算をすると以下の結果がでました。NDは使用した要素の節点数、NEは要素分割数、NNODEは総節点数、U(1,1)とU(2,1)は計算結果、次の2列は厳密解との差になっています。

厳密解との差をグラフにすると以下が得られました。