CAE FEM プログラミング

有限要素法の弾塑性解析Pythonコードを実行してみた

投稿日:2024年7月6日 更新日:

有限要素法、弾塑性解析についてWeb検索をしていた所、Pythonコードを見つけました。

https://qiita.com/Altaka4128/items/86e25d66196dfe7160f3

塑性域に達した時の処理について理解を深める事ができて大変ありがたいです。

実行方法

記事よりコードをコピーして以下のファイルを用意しました。

・main.py
・FEM1d.py
・Node1d.py
・d1t2.py
・Material.py
・Boundary1d.py

コマンドプロンプトでpythonコードを配置したディレクトリにいき、

py -3.12 main.py

で実行できました。
(私の環境ではpython3.9, 3.10, 3.12を共存させています。解析実行は最新のpython3.12で行いました。)

ゼロ割のWarningがでる原因はFEM1d.py内のvecDispFirstがvecDispの参照渡しとなっているからのようです。

vecDispFirst = vecDisp.copy()

としておいたほうが良さそうです。
この修正をしなくても一応解析は最後まで流れて正しい結果が得られます。これは収束判定を2つ行っているので、ゼロ割を起こす方を無視しても、もう一方の

if np.isclose(LA.norm(vecd), 0.0):

で収束判定を正しく行えているからのようです。

Abaqus Learning Edition2023で検証

記事内にAbaqusのinpファイルもあるのでAbaqusでも実行してみました。

1次元モデルなのでただの棒線ですが、変位量が一致する結果が得られました。

-CAE, FEM, プログラミング

執筆者:

関連記事

How to use FFmpeg.wasm. What’s the Cross Origin Isoration?

Nowadays, I got the information that FFmpeg can be used with Javascript and I tried it immediately. …

Android App: How to share video in external storage to MediaStore?

Code Android Studio Electric Eel | 2022.1.1 Patch 2[Java]String path = xxx; //path of the video in t …

有限要素法3次元弾塑性解析Pythonコードを実行してみた

1次元モデルで弾塑性解析の流れを学んだので続いて3次元モデルの解析も流してみました。https://qiita.com/Altaka4128/items/eb4e9cb0bf46d450b03f Py …

FFmpeg.wasmの使い方:ブラウザでアップロードした動画をグレー動画に加工

FFmpeg.wasmの使い方の一例として、・ブラウザで動画をアップロード・動画から音声を抽出・動画からすべてのフレーム画像を取得・すべてのフレームをグレースケールに加工・グレースケール動画を作成・グ …

Android Studio Admob Error “No Ads Config” from Domain “com.google.android.gms.ads

D/MyApp: {“Code”: 3,“Message”: “No ad config.”,“Domain&#82 …

スポンサーリンク