bar非線形解析について

[taji...@athena.ocn.ne.jp]

お世話になっております。

次の塑性解析例題を実施してみました。

http://opencae.gifu-nct.ac.jp/pukiwiki/index.php?plugin=attach&refer=SALOME-Meca%A4%CE%BB%C8%CD%D1%CB%A1%B2%F2%C0%E2&openfile=07-00.pdf

別添にcommファイルを添付致します。

eficasで設定しRUNしましたがエラーが出ます。
messファイルを見ても原因が分かりませんでした。

アドバイス頂けると助かります。

[kengo]

境界条件の面グループがメッシュに設定されていません。
fix, load の面グループを設定してください。

[taji...@athena.ocn.ne.jp]

kengo  様


いつもお世話になっております。

Meshの前にgeometryでload、fixを面で指定してから
メッシュをしたのですが。

もう一度、ｍｅｓｈからやり直すという意味でしょうか。
それともcommファイルの境界条件部分で操作するという
ことでしょうか。




2014年4月21日月曜日 18時55分33秒 UTC+9 kengo:

境界条件の面グループがメッシュに設定されていません。
fix, load の面グループを設定してください。

[kengo]

下のURLに書いているように、メッシュに面グループを貼り付けてください。

https://sites.google.com/site/codeastersalomemeca/home/salome-meca-chutoriaru1/meshmojuru#TOC-Create-Groups-from-Geometry-1


ウィザードでケースを作れば、自動的に貼り付けてくれますが、
ASTK や自分でケースを作った場合などは、貼り付けてくれないと思いますので、
上記の操作をしてグループを貼り付けてください。

[taji...@athena.ocn.ne.jp]

kengo  様


いつもお世話になっています。

ご指導のとおり実施しましたら成功しました。
自動的にfix、loadが反映されていると思っていました。
初めての取り組みでしたので知らなく勉強になりました。

変位もＺ方向5ｍｍまで10分割で弾塑性していました。
応力も同様な結果が得られました。

解説文書では、「STAT_NON_LINE」の下の「ARCHIVAGE」に、
「ARCH_ETAT_INIT（OUI）」があるのですが、私の環境では
これが出てきませんでした。
とりあえず、これを無視してRUNしても成功しました。

アドバイスありがとうございました。

[taji...@athena.ocn.ne.jp]

その後、違うMeshファイルで非線形解析を実施してみました。
同じcommファイルでeficasで計算しました。
loadをDZから水平方向のDXに変更しただけです。

別添のとおりエラーが出ました。
原因が分からりませんでした。
何か気付いた点などアドバイス頂けると助かります。

[kengo]

エラーを見るとメモリが足りないと言っているようです。
メモリの設定値を増やしてみてはどうでしょうか。

mess ファイルではエラーを、計算が止まった時点と、ファイルの最後に書いています。
フランス語ですが、機械翻訳で英語にすれば十分読めるレベルにはなりますので、
google 翻訳などを活用して解決してください。

[taji...@athena.ocn.ne.jp]

kengo  様

いつもお世話になっております。

初期設定の256MBから512MBに増やして計算したところ、
非線形の0.1ステップでメモリー不足で計算終了。
今度は1024MBに増やして実施しています。

規模が大きいMeshのためかステップ計算に時間が掛かっています。
Meshは２次要素は作成していません。
メモリー調整しながら結果待ちの状況でした。

[taji...@athena.ocn.ne.jp]

違うMeshでの非線形解析は計算出来ました。

PostProで変位・相当応力の可視化が出来ました。

ParaVisに変えて可視化したところ変位しか可視化出来ませんでした。

また、VTM形式で外部出力しParaViewで可視化しても変位しか出力されません。

ParaVisやParaViewで変位だけでなく応力も可視化するには

どのように操作すれば良いでしょうか。

[kengo]

下のリンクに簡単ですが使い方を書いています。
まずは自分で調べましょう。

https://sites.google.com/site/codeastersalomemeca/home/salome-meca-chutoriaru1/paravis-paraview-mojuru

[taji...@athena.ocn.ne.jp]

kengo  様

いつもお世話になっております。

説明不足でした。

使い方は分かっているのですが、プルダウンメニューに

RESU_DEPL
はあるのですが、 RESU_SIGM_NOEU	垂直応力とせん断応力
RESU_SIEQ_NOEU	相当応力 （ミーゼス、トレスカ、主応力）

が出てきません。

ParaVisではプルダウンメニューに表示され選択出来ました。

これはBar非線形解析の事例でも同様でした。

ParaVisと同様にプルダウンメニュー表示されないため

悩んでいました。

線形解析ではそのような事がなかったのですが。

[kengo]

私も投稿内容を勘違いしてました。すみません。


考えられる原因としては、そもそも元のmedファイルに応力のデータが入ってない、でしょうか。

上で添付していたコマンドファイルを見ると、CALC_ELEM, CALC_NO で、SIGM, SIEF の ELNO, NOEU が対応してないので、正しく出力されていないと思います。
また、SIGM は線形解析用なので、非線形では SIEF を使ってください。
IMPR_RESU でも SIGM_NOEU, SIEQ_NOEU, を指定していないので、出力されません。

下に修正したものを書いています。
私のほうで試して見ましたが問題なく、Paraview のプルダウンメニューに表示されました。

使用ソフト
　Salome 7.3.0 (Salome-Meca 2014.1相当)
　Paraview 4.1.0


RESU=CALC_ELEM(reuse =RESU,
               RESULTAT=RESU,
               OPTION=('SIEF_ELNO','SIEQ_ELNO','EPSI_ELNO','EPEQ_ELNO',),);

RESU=CALC_NO(reuse =RESU,
             RESULTAT=RESU,
             OPTION=('SIEF_NOEU','SIEQ_NOEU',),);

IMPR_RESU(FORMAT='MED',
          UNITE=80,
          RESU=_F(RESULTAT=RESU,
                  NOM_CHAM=('DEPL','SIEF_NOEU','SIEQ_NOEU',),),);

[taji...@athena.ocn.ne.jp]

kengo　様

いつもお世話になっております。

解析に時間が掛かったため確認に時間が掛かりましたが、今朝確認出来ました。

正しく出力されました。有難うございました。

非線形では SIEF を使用することも知らなく間違っていました。

timeを36000、メモリーを1024MBで計算していました。

メモリーを2048MBにすると計算してすぐにダウンします。

本体は8GB、仮想は6GBにしておりメモリーには余裕があるのですが、

ASTER計算でのメモリーを増やそうとするとダウンしました。

規模の大きい事例も非線形解析が出来ました。

応力ひずみ関数を追加したり、増分ステップを細分化したりと

何度もエラー処理をしましたが、何とか最後まで可視化出来ました。

有難うございました。