プログラムの大枠

FEMの計算の流れ

①要素剛性マトリックス作成
　↓
②全体剛性マトリックス作成
　↓
③境界条件（変位・荷重ベクトル）設定
　↓
④連立方程式の解（変位算出）
　↓
⑤ひずみ・応力の算出

プログラムの流れ

＜データ読み込み＞
・節点、境界データ読み込み
・材料データ読込み
・要素データ読込み

＜マトリックス作成＞
・要素剛性マトリックス作成
・全体剛性マトリックス作成
・境界条件設定

＜計算＞
・連立方程式の解（ガウスの消去法）
・応力の算出

＜出力＞
・力の出力
・変位の出力
・応力の出力

例

fem = FEM()
fem.readNodeData()
fem.readMaterialData()
fem.readElementData()
fem.setElementMatrix()
fem.setBandWidth()
fem.setTotalMatrix()
fem.setBoundaryCondition()
fem.saveForce()
fem.saveTotalMatrix()
fem.solve()
fem.saveDisp()
fem.setStress()
fem.saveStress()