プログラムの大枠

Pythonで有限要素法プログラムの作成



プログラムの大枠

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()
ページトップへ