* SALOME-Mecaの活用技術(MK2) [#n9a60459] [[オープンCAE初心者勉強会]]~ Code_Asterの使用方法、事例を解説します。~ #contents **ASTK(GUI)、Eficas(commファイルエディター)、公式Validationの利用方法など [#kc395628] ***ASTKの使用方法、公式Validationの利用、Eficasでのcommファイル編集例 [#jd909c8b] -資料~ &ref(ASTK_Intro.pdf);~ ~ ***Eficas Importerの使用 [#nea2d49f] EficasではGROUP_MAなどで複数のグループ名を記入できますが、これをテキストファイルから読み込む方法です。~ -資料~ &ref(Eficas_importer.pdf);~ ~~ **要素・境界条件の設定 [#mbf45c2f] ***シェル要素(DKT、COQUE_3D)の扱い方、メッシュアセンブリ(LIAISON_MAIL) [#f7ee723a] シェル要素として、DKTとCOQUE_3Dがあり、大きな回転などを扱う場合に、COQUE_3Dを使うようです。~ |CENTER: |CENTER:メッシュ|CENTER:変位|CENTER:回転|CENTER:歪み|CENTER:計算| |CENTER:DKT|CENTER:TRIA3,QUAD4|CENTER:small|CENTER:small|CENTER:small|CENTER:Linear,Non linear| |CENTER:COQUE_3D|CENTER:TRIA7,QUAD9|CENTER:large|CENTER:large|CENTER:large|CENTER:Linear,Non linear| -資料~ DKT, LIAISON_MAIL:&ref(Shell_LInkMesh.pdf);~ COQUE_3D:&ref(Shell_Coque3D.pdf);~ -入力ファイル~ DKT, LIAISON_MAIL:&ref(Shell_LInkMesh.tar.gz);~ ~ ***断面性能の計算(MACR_CARA_POUTRE) [#g453c1cc] -資料~ &ref(MACR_CARA_POUTRE.pdf);~ ~ ***ビーム要素(BARRE、POU_D_E)の扱い方、ビーム要素とシェル要素・3D要素との結合(LIAISON_ELEM) [#wc5e5759] ビーム要素は2つのタイプがあり、回転の自由度を持たない:BARRE、回転の自由度を持つ:POU_D_E 他、を使い分けます。 |CENTER:Modelisation|CENTER: |CENTER:自由度| |CENTER:BARRE|CENTER:ピン接合|CENTER:DX, DY, DZ| |CENTER:POU_D_E 他|CENTER:剛接合|CENTER:DX, DY, DZ, DRX, DRY, DRZ| -資料~ BARRE:&ref(Bar.pdf);~ POU_D_E, LIAISON_ELEM:&ref(Beam_Link.pdf);~ -入力ファイル~ BARRE:&ref(Barre.zip);~ POU_D_E:&ref(Beam.zip);~ POU_D_E & LIAISON_ELEM:&ref(Beam_link.zip);~ *POU_D_Eなどの剛接合要素では応力が大きくなると、後処理段階でSIEF_ELGAなどの応力を計算した際にエラーが発生するバグがあります。~ ~~ **ポスト処理 [#lcf25092] ***MECA_STATIQUEの場合の一例(応力・等価応力の出力、Post-Pro・Paraviewモジュールでの結果表示) [#faac6431] -資料~ &ref(Post_Lin_ver10.pdf);~ -commファイル~ &ref(postCA10_lin.comm);~ ~ ***Salome-Meca 2011.2 Asterモジュールの計算結果をParaviewモジュールで表示させる [#a1f85f5b] 計算結果は作業フォルダの拡張子:rmed のファイル(linear-static.rmed など)に保存されているので、この拡張子を rmed → med に書き換えれば、Paraviewモジュールで表示できます。~ ~~ **並列化、AmazonEC2の利用 [#k70d8970] ***Code_Aster の並列化 [#u57d314e] [[Code_Aster の並列化]](PDFの内容のうち、コマンドラインが中心の2-10をリンク先に清書しました~ -資料~ &ref(Parallelize_3.pdf);~ PAR10.3のコンパイルの前に、/opt/aster/etc/codeaster/asrun の書き換えを追加しました。~ これを行わなければコンパイルでエラーが出るようです。 -添付ファイル~ Code_Aster10.3用:&ref(setup.cfg);~ mumps-4.9.2用:&ref(Makefile.inc);~ PAR10.3用:&ref(config.txt);&ref(profile.sh);~ ~ *バージョン11.0の並列化は、上記の資料と同様の手順で実行できます。ただし、11.0 と 10.3 で config.txt が違うため、PAR10.3 の config.txt は流用できません。そのため、STA10.3 → PAR10.3 で行ったものと同様の書き換えを、PAR11.0 の config.txt に行う必要があります。~ ~ ***Ver. 10と11、ソルバー、並列化の比較(100万自由度 STAT_NON_LINE での計算) [#xd0c55c7] Ver. 11.0 のほうが速く、ソルバーはMUMPSのほうが速いです。~ |CENTER: |CENTER:STA10.3|CENTER:PAR10.3|CENTER:PAR10.3|CENTER:NEW11.0|CENTER:PAR11.0|CENTER:PAR11.0| |CENTER:Processors|CENTER:1|CENTER:1|CENTER:4|CENTER:1|CENTER:1|CENTER:4| |CENTER:Solver|CENTER:MULT_FRONT|CENTER:MUMPS|CENTER:MUMPS|CENTER:MULT_FRONT|CENTER:MUMPS|CENTER:MUMPS| |CENTER:time [sec]|CENTER:3619|CENTER:1474|CENTER:1118|CENTER:2734|CENTER:968|CENTER:766| ~ ***AmazonEC2インスタンスのクラスタ化、並列化Code_Asterの使用 [#f9b524ad] AmazonEC2のサービスは従量課金制ですので、料金体系を良く理解した上で利用してください。オンデマンドインスタンス、データ転送、EBSボリューム、Elastic IP、Route 53は課金されます。[[Amazon EC2 料金表:http://aws.amazon.com/jp/ec2/pricing/]] -資料~ &ref(AmazonEC2_Clusterize.pdf);~ ~ ***WindowsでX forwardingを行う [#q32b8665] AmazonEC2などのリモートマシン上のGUIアプリケーションを、ローカルマシンのWindowsから利用する方法です。 -資料~ &ref(Xforwarding_Win.pdf);~ ~ ~~ **Salome-Meca、Code_Aster、Salomeのバージョンについて [#t6e30eb1] ***Salome-Meca のバージョン [#ud71f6da] |CENTER:Salome-Meca|CENTER:bit|CENTER:Code_Aster|CENTER:Salome| |CENTER:2011.2|CENTER:64bit|CENTER:STA10.5|CENTER:6.3.0| |CENTER:2011.1|CENTER:64bit|CENTER:STA10.3|CENTER:5.1.5| |CENTER:2010.2|CENTER:32bit|CENTER:STA10.2|CENTER:5.1.4| ~ ***Code_Aster のバージョンについて [#jc083697] 2011年8月の時点で、Salome-Meca 版、スタンドアローン版の主なバージョンとしては、STA10.2、STA10.3、NEW11.0、STA10.5があります。(STA10.5はスタンドアローン版はありませんが、STA10.5 = NEW11.0.10です。)このうち、[STA10.2、STA10.3]と[NEW11.0、STA10.5]では、コマンドファイルのポスト処理の書き方に違いがあります。(それ以外のコマンドファイルの書き方は変わっていないようです。)~ 具体的には下表の例のように、CALC_ELEM、CALC_NO の OPTION の書き方が違います。そのため、STA10.2やSTA10.3用に作ったコマンドファイルを、NEW11.0やSTA10.5で計算する場合は、EFICAS などでポスト処理の部分を書き換える必要があります。逆もまた然りです。~ ~ |CENTER:バージョン|CENTER:応力(要素)|CENTER:等価応力(要素)|CENTER:応力(節点)|CENTER:等価応力(節点)| |CENTER:STA10.2、STA10.3|CENTER:SIGM_ELNO_DEPL|CENTER:EQUI_ELNO_SIGM|CENTER:SIGM_NOEU_DEPL|CENTER:EQUI_NOEU_SIGM| |CENTER:NEW11.0、STA10.5|CENTER:SIGM_ELNO|CENTER:SIEQ_ELNO|CENTER:SIGM_NOEU|CENTER:SIEQ_NOEU| ~ ***MED のバージョン対応 [#k800611a] Salome 6 は、MED 3 でしか出力できないので、Salome 6 で作った MED ファイルを用いて計算する場合は、NEW11.0, STA10.5 を使う必要があります。逆に、NEW11.0, STA10.5 で出力した MED ファイルは、Salome 6 でしか読めません。~ MED : HDFライブラリ(Salome)とCode_Aster でデータをやり取りする際の中間データ(メッシュなど)~ ~ |CENTER: |CENTER:read|CENTER:write| |CENTER:NEW11.0, STA10.5|CENTER:2, 3|CENTER:3| |CENTER:STA10.2, STA10.3|CENTER:2|CENTER:2| |CENTER:Salome 6|CENTER:2, 3|CENTER:3| |CENTER:Salome 5|CENTER:2|CENTER:2| ~ ~ **Code_Aster 参考ページ [#yecf1803] ***Code_Aster公式、CAElinux Contribute [#pf504676] Code_Aster Documentation : http://www.code-aster.org/V2/spip.php?rubrique19~ 機械翻訳ですが、英語化公式ドキュメントが公開されています。~ Code_Aster Forum : http://www.code-aster.org/forum2/~ 開発者がいますので非常に参考になります。言語は英語とフランス語ですが、英語の投稿の方が多いようです。~ ~ Code_Aster Wiki : http://www.code-aster.org/wiki/doku.php~ 有志が作成しています。完全に整理されてはいませんが参考になります。~ ~ Contrib:Claws/Code Aster : http://www.caelinux.org/wiki/index.php/Contrib:Claws/Code_Aster~ Contrib:Kees Wouters : http://www.caelinux.org/wiki/index.php/Contrib:KeesWouters~ 有志が公開されている解説書です。コマンドファイルを公開していますので非常にわかりやすいです。~ ~ *** 上記の資料をもとに柴田が Code_Aster のコンパイルを練習した時のメモ [#s87e5613] DEXCS-Salomeの有効活用:Code_Asterのコンパイル(1):http://blog.goo.ne.jp/dalab-gnct/e/fc5c4581f73b06bce7a977ae56379c2b~ DEXCS-Salomeの有効活用:並列版Code_Asterのコンパイル(1):http://blog.goo.ne.jp/dalab-gnct/e/d3a33ad805388450e26dfdbd5049f4f6~ DEXCS-Salomeの有効活用:並列版Code_Asterのコンパイル(2):http://blog.goo.ne.jp/dalab-gnct/e/f12c955436191d9783516214d7573f99~ DEXCS-Salomeの有効活用:並列版Code_Asterのコンパイル(3):http://blog.goo.ne.jp/dalab-gnct/e/fdc6d30b81b008f9eec405945c4d9467~