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, プログラミング

執筆者:

関連記事

Xcode iOSアプリでGUIをコードだけで作成する Storyboard, Sceneなし

アプリメンテナンスをしている際、どうしようもないビルドエラーに遭遇して解決できず、新規プロジェクトで作り直すという事がありました。私のiPhoneアプリ開発歴は結構長く、始めたのは2013年あたりです …

Found an unexpected Mach-O header code: 0x72613c21 への対処 Xcode, Admob, xcframework

日々変わっていくAdmobの仕様への対応にとても苦労しています。そんな中、アプリをApple Store Connectへ提出するためXcodeにてArchiveを作成してValidateしようとした …

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 …

iPhoneアプリ[FEM BLOCKi]で有限要素解析をやってみた第3回 固有値解析

世界一、お気楽な有限要素法講座の第3回目の記事です。今回は固有値解析をやってみます。これまでと同様にWeb上から例題を探し、iPhoneアプリFEM BLOCKiで実行してみます。 なお、本記事中では …

Android studio|FFmpegでやりたい放題するための準備

動画を利用したアプリを作成するにあたりFFmpegを使用できると大変便利です。Web上に多数の情報がありますが、数年前の記事は役に立たなかったり、実用的な情報が少なかったりします。本記事ではAndro …