スポンサーサイト

上記の広告は1ヶ月以上更新のないブログに表示されています。
新しい記事を書く事で広告が消せます。

FEM_BLOCK(ver2.06仮):陽解法その1

勢いで陽解法を組み込んでみました。
意外と簡単?だったりしました。
苦労して幾何学非線形静解析に取り組んだおかげでしょうか。
中央差分法を用いています。
片持ち梁の自由端に荷重。重力も一応入れておきました。
ひずみから判断すると良さそうな結果が得られました。



梁寸法:1×1×8[M]
材質:一般鋼材
要素:4面体2次
境界条件:片端固定 荷重:自由端に2E8[N]
時間ステップ:1E-5[s]
解析最終時間:1E-1[s]
計算時間:2分くらい

初めて陽解法を行ってみたのですが
時間ステップをかなり小さくする必要があることがわかりました。
大きくするとあっという間に発散します。

参考文献に従って質量行列は集中質量を用いました。
逆行列演算、連立一次方程式を一切解かなくて良いので
1ステップの計算はかなり短時間で済みます。
しかしステップ数は相当多いのでそれなりに時間はかかります。

アニメーションから判断すると幾何学非線形も考慮されているようなので
ついに念願のFEM非線形ソルバーを手に入れたことになるのでしょうか?

スポンサーサイト

コメントの投稿

非公開コメント

カレンダー
07 | 2017/08 | 09
- - 1 2 3 4 5
6 7 8 9 10 11 12
13 14 15 16 17 18 19
20 21 22 23 24 25 26
27 28 29 30 31 - -
プロフィール

uglab

Author:uglab
ソフト開発会社

【Twitter】
フォローお願いします

【You tube動画一覧】
チャンネル登録お願いします


最新記事
最新コメント
最新トラックバック
月別アーカイブ
カテゴリ
FC2カウンター
検索フォーム
スポンサードリンク
    【スポンサードリンク】

リンク
RSSリンクの表示
ブロとも申請フォーム

この人とブロともになる

QRコード
QRコード
上記広告は1ヶ月以上更新のないブログに表示されています。新しい記事を書くことで広告を消せます。