最新の情報はこちら。 以下は、古い内容です。

Salome-Meca2017 - ちなみに、Salome-Meca2018は、見た目、操作方法など、以下とほぼ変わらないと思います。

目次

• 資料
• 準備
• インストール
• unicode decode error
• 日本語化
• 片持ち梁
• 「geometry」で形状作成
• 境界条件、荷重条件などを設定する箇所をグループ化
• 移動方法
• 「mesh」で要素分割
• 四面体要素を使う場合
• 直方体要素を使う場合
• メッシュを作成
• 「geometry」でグループ化したグループをメッシュに反映させる
• 以下は落ち穂拾い
• 「AsterStudy?」で材料定数、境界条件、荷重条件などの設定
• 「AsterStudy?」全般について
• 入力ファイルを作る
• 各種設定の手順
• 計算を実行する
• 「ParaVis?」で計算結果を確認
• 起動時に不具合が生じた時
• 回転・並進移動・拡大縮小
• 変位図を描く
• コンター図を描く
• コンター図の色を変える
• コンター図の微調整
• 背景を変える
• Salomeで作ったグループのみ表示する
• 断面を切る
• とりあえず結果を出力する
• 結果を表に出力する
• 連続画像を出力
• ParaVis?上で簡単なグラフ作成
• 単純梁
• 「geometry」で形状作成
• 拘束線や載荷線を作る
• 境界条件、荷重条件などを設定する箇所をグループ化
• 「mesh」で要素分割
• 「AsterStudy?」で材料定数、境界条件、荷重条件などの設定
• 入力ファイルを作る
• 各種設定の手順
• 計算を実行する
• 「ParaVis?」で計算結果を確認
• ASTK
• commファイルのバージョンを変換する（未解決）

ちなみに、Salome-Meca2018は、見た目、操作方法など、以下とほぼ変わらないと思います。

資料

• SALOME-Mecaの使用法解説
• JSWG - Salome-Meca活用研究会
• AsterStudyの使い方

準備

インストール

• Code_Asterのページから、Salome-Meca 2017を"~/kouzou/.Salome-Meca"あたりにダウンロード。
• Salome-Meca-2017は2017/09リリースなのだが、その2ヶ月後にSalome-Meca2017.0.2にリリースされて、元のSalome-Meca2017がダウンロードページから消えてしまった。いつもなら古いバージョンもダウンロードページに残されるのに、それが消されているということは、9月リリースのSalome-Meca2017には何か致命的な欠点でもあるだろうか？Code_AsterのNewsあたりをざっと探してみたが、理由が見当たらない。一応、ダウンロードページのSalome-Meca 2017.0.2のURLの".0.2"を消すと、元の2017がダウンロードできるみたいだけど。
• "$ tar xvzf alome-Meca-2017-LGPL-3.tgz"で、展開。
• "$ ./Salome-Meca-2017-LGPL-3.run"で、インストール開始、
• "Enter the salome_meca install directory [default=/home/kouzou/salome_meca] :"は、後でディレクトリごと動かせばいい(?)ので、とりあえず"Enter".
• "salome_meca is in English, do you want to be in French ? [y/N] :"は、英語がいいので"N".
• インストール終了後、アプリケーション -> 教育・教養 -> salome_meca V2017 で立ち上げる。

unicode decode error

 これをやらないと日本語環境のUbuntuでは計算ができない。

• sitecustomize.pyを"~/salome_meca/appli_V2017/lib/python2.7/site-packages/salome"にダウンロードする。

• 参考URL
• https://twitter.com/mmer547/status/903788754506309632
• https://twitter.com/freekaos/status/905068369237336066

• ちなみに、システム -> 設定 -> ユーザ向け -> 言語サポート -> 地域フォーマット -> 日付の表示：English -> システム全体に適用 -> 再起動　でも使えるようになる。

日本語化

• file -> preferences -> SALOME -> Current language -> ja -> Apply -> OK
• 一旦Salome-Mecaを閉じて、また立ち上げる。

片持ち梁

幅方向：x , 載荷方向：y , 軸方向：z とした。

「geometry」で形状作成

• Salome-Mecaの起動：アプリケーション->教育・教養->salome_meca V2017

• モデルの新規作成：SALOMEと書かれている上部プルダウンメニューからGeometryを選択->
「モジュールのアクティブ化」ウィンドウが現れるので「新規作成」を選択

• 新しいエンティティ->基本図形->ボックス->「ボックスの作成」ウィンドウが現れる

• Dx, Dy, Dzに寸法を入れる。単位は自分で決められるが、mm(長さ)-N(力)系にするのがいいだろう（ここでは例えば、10*10*100とする）。後藤の好みでは、梁軸がz, たわみ方向がyとか（この辺は好み）。->「適用して閉じる」

• オブジェクトブラウザーにBox_1（上でボックスにつけた名前）が現れる。

• 単純梁の場合「geometry」で形状作成の続きはこっち。

境界条件、荷重条件などを設定する箇所をグループ化

• 目的とする物体(box1)を右クリックして、「create group（グループを作成）」

• 点(+), 線(/), 面(□), 立体(立方体マーク)のどれかにチェックをいれて、固定や載荷をするところを選択。固定端の面固定なら、面を選択。
• 後で何度もここを選択するので、名前の「Group_1」をfixとかkoteiとかの分かりやすい名前に書き換える

• 超重要：バージョンによる(Salome-Meca2013では)かもしれないが、名前は8文字以下にしないと原因不明のエラーが出る!!
• 確か、Salome-Meca2015〜16ぐらいからは大丈夫になった。
• 固定したい面を選択するために、ボックスを回転させたり移動させたりする方法は、以下の「移動方法」を参照

移動方法

• 回転は、玉に回転矢印の書いてあるところをクリックしてから、左クリックを押しながら回転させる

• 回転は、cntlと右クリックを押しながら動かす
• 並進移動は、cntlとスクロールボタンを押しながら動かす
• 拡大縮小は、ただスクロールする。あるいはcntlと右クリックを押しながら動かす

• 座標の原点側を固定端にするのであれば、回転させて固定端面を選択してクリック->選択された水色の外枠線が白に変わったのを確認して追加

• 「追加」すると「31」と数値が表示されるのを確認して（違う面を選択した場合は違う番号になる）->「適用して閉じる」

• 一つのツリーにまとまっていればOK

• 同様の手順で載荷面を作る
• グループの名前はわかりやすく saika とする

• そうすると、グループを作成されていないボックスの側面が消えてしまったりすることがあるが、その場合は、オブジェクトブラウザーのボックスの目玉マークをクリックすれば表示される

「mesh」で要素分割

• オブジェクトブラウザーの(固定面や載荷線などを含んだ)ツリーの元(例えばBox_1)を選択してから「Geometry」が表示されている上部のプルダウンメニューをクリックして「Mesh」を選択する。

• 上部ツールバーの「メッシュ」->「メッシュを作成します」->「メッシュを作成」ウィンドウが現れる

• 名前はひとまずMesh_1でいい
• オブジェクトブラウザー内の選択したいボックス(Box_1)を選択すると、「メッシュを作成」ウィンドウのジオメトリに Box_1みたいに表示されることを確認。

四面体要素を使う場合

• 拘束線や載荷線を含めてメッシュ分割する場合は、Partitionを選択（単純梁の載荷線、拘束線などを作るとジオメトリがPartitionとなるが、Partitionが生成されていない場合は、ジオメトリにBox_1を選択すればよい）
• 詳細設定セットの割当->3D Automatic Tetrahedralization を選択

• 要素長設定ウィンドウ(Max Size)が出たら、ある程度の精度が必要な場合は長さ1を目安に設定し、「OK」->「適用して閉じる」

• オブジェクトブラウザーにMeshが追加される

• 下の「メッシュを作成」に進む

直方体要素を使う場合

• 3Dアルゴリズム：Hexahedron(i,j,k)
• 詳細設定セットの割当から3D Automatic Hexahedralization を選択。
• セグメントの数を適当に（適切に）設定
• ↑の方法だと軸方向の分割数だけ荒くなってしまうので，予めGeometryで辺ごとにグループを作っておいて，MeshのSubmeshで辺ごとの分割数を指定するのがいい．それが面倒なら，直方体のLengthを各辺の約数にするとか．
• 適用して閉じる

メッシュを作成

• Mesh_1を右クリックして「メッシュを作成」

• 節点数や要素数などが計算されて、3D(ボリューム)とかの箱の中に、色々と数字が入っていれば、たぶん大丈夫。ノード以外の箱がほとんど0だったら、計算されてない。

• 拡大してメッシュ分割されているのを確認する。

「geometry」でグループ化したグループをメッシュに反映させる

• オブジェクトブラウザーのMesh_1を右クリックし、「ジオメトリーのグループを作成」を選択

• 「ジオメトリからグループを作成」ウィンドウが開くので、オブジェクトブラウザーのBox_1で作った固定面や載荷面のグループ（kotei, saika)をクリックする（複数を1個ずつ選択する場合はcntlキーを押しながら。最初のグループを選択した後、最後のグループをshftキーで選択すれば、すべてまとめて選択できる）
• ジオメトリにkoteiとsaikaが追加されたのを確認->適用して閉じる

• オブジェクトブラウザーのMesh_1のツリーにkoteiとsaikaの「Groups of Faces」が追加されたことを確認。

以下は落ち穂拾い

• ビューアーのタブで、geometry表示でない？方の「VTK scene:1 - viewer:1」みたいなもう一つのタブを開くとメッシュ分割が表示されている。
• デフォルトは、ノード表示で点々しか表示されないので、物体を右クリックして、表示モードを「網かけ」とかにすると物体表面のメッシュが表示される。

• ちなみに、メッシュ分割後も右クリックで「create group」できる。
• AsterStudy?へ行く前にGeometryで作った固定面、載荷面などを「ジオメトリのグループを作成」などでグループ化しておく。
• 要素かノードか（要素の方が、面載荷などはやりやすい）
• 固定面、載荷面など複数の面をグループ化する場合は、Cntlを押しながら、左のウィンドウから選択
• meshを右クリックして、「convert to quadratic」をクリックし、applyすると、(中間節点を設けた)2次要素になるので粗いメッシュでも精度がよくなる。

「AsterStudy?」で材料定数、境界条件、荷重条件などの設定

• AsterStudy?のアイコンをクリックしてAsterStudy?を立ち上げる。
• プルダウンメニューからAsterStudy?を選んでもよい

• 上記の「geometry」でグループ化したグループをメッシュに反映させる作業は重要なためやったことを確認する。
• これをしないと、AsterStudy?で境界条件、荷重条件などを設定するときに、グループ名が表示されない。

「AsterStudy?」全般について

• 右側のEdit comandウィンドウでは、最低限入力が必要なところには赤字で＊が付いている。
• 何かおかしな設定をしていたら、文字が赤く表示される。
• 「適用」あるいは「OK」すると、左側下のInformationウィンドウにcommファイルに書かれる内容が表示される。
• 後から設定を編集し直したいときは、"Data Settings"にある直したい項目をダブルクリックする。
• 右の"?"を押すと公式ドキュメントが開かれる。ただし、フランス語で書かれているので、URLの …/fr/… を …/en/… にすると英語版になって読みやすいかも。

入力ファイルを作る

• ツールバーの右上の方にあるAdd stage("{}+"こんなマーク)をクリック。すると、Data SettingsにStage_1が追加される。

• 以下の作業では、左側のData Settingsウィンドウ内のStage_1を選択したまま、ツールバーの"メッシュ△▽"みたいなのをクリックする。Stage_1を選択していないと、"△▽"は押せない。

各種設定の手順

AsterStudy?上部にならんでいる、「メッシュ」「Model Definition」「Material」のプルダウンメニューを順番に選択し、それぞれの項目について、以下の手順で設定していく。

• メッシュ -> Read a mesh

• Mesh file location->Mesh_1（Meshで作ったもの）
• Mesh file format->Med

• 20='0:1:2:4'の意味はよくわからない（違う数字になる場合もある）
• OKをクリックすると、左側のData Settingウィンドウにmesh Read a meshができる。設定内容は、下のinformationウィンドウに表示される。

• Model Definition -> Assign finite element

• メッシュ : mesh(LIRE_MAILLAGE) (Read a mesh で付けた名前)が自動的に選ばれる（ちなみに、LIRE MAILLAGEはRead a meshのフランス語）

• At Least One -> Finite element にチェック

• "0 items +"の+をクリックすると、"1 item"に変わり、

• Editが現れるのでクリックする。

• Exactly One : Everywhere : Yes (1材料なら)
• （2材料以上なら）：Group of element -> Edit -> …
• Phenomenon -> Mechanicを選択すると、

• "0 items +"が現れるので、"+"のところをクリックすると、"1 item +"に変わり、その下に空欄と▽プルダウンが現れるので、▽プルダウンをクリックすると、

• プルダウンメニューが現れるので3D (立体要素なら)を選択

• OKをクリックすると、プルダウンが消えて"At Least One"の画面に戻るので、もう1回OKをクリックすると

• 左側のData Settingsに"model"が追加される。

• Material -> Define a material
• Linear isotropic elasticity にチェック(線形弾性体なら) -> Edit...

• Young's modulus : ヤング率を入力
• Poisson's ration : ポアソン比を入力

「OK」で閉じる。

• Material -> Assign a material
• At Least One : Model にチェック -> model(AFFE_MODELE) (Assign finite element で付けた名前)

• Material assignement -> "+" -> Edit...
• Everywhere : Yes (1材料なら)
• Material -> "+" -> mater(DEFI_MATERIAL) (Define a material で付けた名前)

「OK」で全て閉じる。

• BC and Load -> Assign mechanical load

• At Least One : Enforce DOF -> "+" -> Edit...
• At Least One : Group of element -> Edit... -> kotei (Geometryで付けた固定端のグループ名) -> 「OK」
• At Least One : DX・DY・DZにチェックを入れ、0.0を入力する。

• 「OK」で閉じる。
• At Least One : FORCE_FACE (面載荷なら) -> "+" -> Edit...

• At Least One : Group of element -> Edit... -> saika (Geometryで付けた載荷面のグループ名)
• At Least One : FY（載荷方向）にチェックを入れ、$\frac{任意の荷重（ここでは100N)}{載荷面の面積}$の値を入力する。
• FORCE_ARETE(線載荷)の場合は載荷線の長さで割る。

• 「OK」で閉じる。
• Model : model(AFFE_MODELE) (Assign finite element で付けた名前)

• 「OK」で閉じる。
• Analysis -> Static mechanical analysis
• At Least One : Material field : fieldmat(AFFE_MATERIAU) (Assign a material で付けた名前)
• Model : model(AFFE_MODELE) (Assign finite element で付けた名前)
• Loads -> "+" -> Edit...
• Load : load(AFFE_CHAR_MECA) (Assign mechanical load で付けた名前)

• 「OK」で閉じる。
• Post Processing -> CALC_CHAMP
• Result : reslin(MECA_STATIQUE) (Static mechanical analysis で付けた名前)
• CONTRAINTE -> "+" -> "▽" -> SIGM_NOEU (節点の応力(6成分))
• DEFORMATION -> "+" -> "▽" -> EPSI_NOEU (節点のひずみ)
• CRETERES -> "+" -> "▽" -> SIEQ_NOEU (節点の応力(ミーゼス応力や主応力など))
• ちなみにNOEUは節点(NOEUd(フランス語))、ELGAはガウスの積分点(ELement GAuss point)、ELNOは要素(のどこか?)(ELement Node)だとかで、ParaVis?でコンター図を描きたいときは Point Data であるNOEUを出力する必要がある。ELGAやELNOは、ParaVis?上でSpredSheetView?を開き、Field Data や Cell Data を選択すれば値が見れる。

• 「OK」で閉じる。
• Output -> Set output results
• Result file location : 計算結果のファイル名(.med)と保存場所を指定する。
• '▽△' みたいなところではなく、右の '...'のところ
• Format : Med (ParaVis?で結果を見るのなら)

• Results -> "+" -> Edit...
• At Least One : Result : reslin(MECA_STATIQUE) (Static mechanical analysis で付けた名前)
• NOM_CHAM にチェック -> "+" -> "▽" -> DEPL (節点変位)

• 一旦OKを押し、"Result file location"のFormatとResultsにチェックされたウィンドウに戻る。そこで、Resultsの1 itemになってるところをの+をクリックして2 itemsにして、以下の操作。

• Results -> "+" -> Edit...
• At Least One : Result : unnamed(CALC_CHAMP) (CALC_CHAMP で付けた名前)
• (以下は3itemsにするということかな？後藤)(そうです。「欲しい数だけ」itemsです。近藤)
• NOM_CHAM にチェック -> "+" -> "▽" -> EPSI_NOEU (節点のひずみ)
• NOM_CHAM にチェック -> "+" -> "▽" -> SIEQ_NOEU (節点の応力(ミーゼスや主応力など))
• NOM_CHAM にチェック -> "+" -> "▽" -> SIGM_NOEU (節点の応力(6成分))

• この通りに作ると、commファイルがこんな風になる。

計算を実行する

• History View
• Case_1を選択 -> □Stage_1 をクリック -> |>Run
• Basic -> Memory : メモリに余裕があるのなら上げる(というか、要素数によっては上げないと計算できない)
• Basic -> Time : 計算が長くかかりそうなら上げる
• Run

• 計算が上手くいくと計算中に回転している砂時計が緑色、駄目だと赤色や黄色の○に変わる。

• 結果の一部は見られる(緑色になる)けど、何か致命的とまでは言えない間違い(警告)があるときは、logにこんなふう↓に書かれる。(正常なら"Success"のあとに"+Warn"は書かれない。)

Stage "RunCase_1:Stage_1" calculation succeeded with state Success+Warn

• エラーや警告のメッセージ等は、LogをConsoleにするかShow message file を開くと見れる。ちなみに、エラーは"！"で括って書かれる．

「ParaVis?」で計算結果を確認

• AsterStudy? -> Mesh -> ParaVis?　
• パイプのブラウザーで右クリック -> Open -> AsterStudy?のOutputで設定した出力ファイル(*.med)を選択する。

起動時に不具合が生じた時

• Case View
• Data files
• 〇〇.med を右クリックして、Open in ParaVis?

• 不具合が生じた時は，次のものを再び開いて，元に戻す．
• 表示 -> Windows
• パイプのブラウザー
• オブジェクトインスペクター
• カラーマップエディタ (<- 好みで追加)
• 表示 -> ツールバー
• Active Variable Controls
• Camera Controls (<- 好みで追加)
• Current Time Controls (<- 好みで追加)
• Common (<- 好みで追加)

回転・並進移動・拡大縮小

• 回転：左クリックを押したままマウスを動かす
• 並進移動：スクロールを押したままマウスを動かす
• 拡大縮小：スクロールを回すか，右クリックを押したままマウスを動かす

• たまに物体がいなくなることがある．その時は次のようにすると元に戻る．
• Set view selection to -○ (上の方にある+X, -X, +Y, -Y, +Z,

 -Z のどれか) -> Reset(Set view selection の3つ左)

変位図を描く

• フィルター -> Common(Alphabeticalでもいい) -> Warp By Vector
• DEPL(変位)を選択する(ベクトルになっているデータしか選択できないので，通常DEPLしか選択できない．)
• scale(変形の倍率)を適宜変えて，Apply．

コンター図を描く

• 初期状態で上のほうに"● Solid Color"と書かれているところを見たいもの(変位なら 〇〇_DEPL とか)に変える．
• その右隣の，Magnitude(ノルム)から見たいもの(Y方向変位なら DY とか)に変える．

コンター図の色を変える

• オブジェクト インスペクターのColoring -> Choose Preset(ハートマーク) -> 好きなのを選んで，適用
• ここ的な見やすさを求めるのなら，"Gray scale"か"Linear YGB", "Linear Green", "Linear Blue" などがいいかも．

コンター図の微調整

• オブジェクト インスペクターのColoring
• C(rescale to custom data range)：コンター図の範囲の上下限値を設定．
• 目(rescale to visible data range)：今見えている範囲で，コンター図の範囲を丁度良くしてくれる．見たい箇所のみ表示させて使うと非常に便利．
• e(Edit Color Legend Parameters)：カラフルな棒(Scale bar)の大きさや向き(縦長か横長か)，文字の大きさや数，フォント等を調整．
• View -> Windows -> Color Map Editor
• Number Of Table Values：色の分割数．適度に荒い方が見やすいときがある．
• Maximam Number of Labels と Number Of Table Values の数に合わせると、csvに出力しなくても大体の値が読み取れる。

背景を変える

• オブジェクト インスペクター -> Background(一番下にある) -> Colorを白っぽいのに

• ちなみに，ここで Single color を Image に変更して画像ファイルを選ぶと，背景をその画像にすることもできる．

Salomeで作ったグループのみ表示する

• Filter -> Alphabetical -> Extract Group -> 見たいグループだけにチェックを入れる

断面を切る

• パイプのブラウザーのrmedファイルを選択したまま、Filter -> Common -> Slice
• 見たい断面の、原点座標(origin)を入力し、面の法線(nomal)を選択する．

とりあえず結果を出力する

• Select Points ○○ (変位を見たい節点を選択する) -> 右上の"ロロ"("日"でもいい) -> Spreadsheet View -> （右上の方にある）Show only selected Elements.

結果を表に出力する

• Select Points 〇〇 -> 右上の"ロロ"("日"でもいい) -> Spreadsheet View -> Show only selected Elements.　-> (File の)Export Secen
• Select Points On : 表面の節点だけ選択する
• Select Points Through : 内部の節点も選択する

連続画像を出力

• File -> Save Animation -> 細かな設定を色々．そのうち書く．(誰か書いて) -> pngとかjpgとかを選択
• その後，端末にて，ImageMagick?を使う．

convert -delay 10 moto*.jpg filename.gif

• "-delay X"は再生コマ間隔がX/100秒となる。なくてもいい．
• 画像と同じで"eog filename.gif&"で観られる.
• サイズも"convert -geometry 300 filename.gif"とかで後から変えられる(元のサイズのままだと，大きくて，重いので不便)．

ParaVis?上で簡単なグラフ作成

• パイプのブラウザーのrmedファイルを選択したまま，Filter -> Data Analysis -> Plot Over Line
• PropertiesのPoint1とPoint2は、線分の起点と終点なので、適宜座標値を入力。これでグラフの横軸が決まる。
• 例えば、梁の中立軸の両端の座標値とかを入力。
• DisplayのSeries Parametersを見たい変位や応力(例えば、荷重方向の変位なら"〇〇_DEPL_DY"だけ)に選択。
• ViewのLeft Axis Lange(横軸)やBottom Axis Lange(縦軸)の値を変えて、グラフの範囲を適宜修正。
• チェックを外すと、自動で調度良い範囲になる。

単純梁

幅方向：x , 載荷方向：y , 軸方向：z とした。

「geometry」で形状作成

• 片持ち梁の「geometry」で形状作成と途中まで同じ。
• Box_1を作成するところまで（グループ化はまだ必要無い）。
• その後↓

拘束線や載荷線を作る

• 新しいエンティティ->基本オブジェクト->点　で点を2個作る
• 線が必要な場所の座標を考えて点を作るか、適当な場所に作ってから移動させる。
• 画面一番上のツールバー「操作」 -> 「変形」 -> 「変換」で移動させられる。
• 新しいエンティティ -> 基本オブジェクト -> 線で点2個を結んで線を作る
• これを必要な回数繰り返し線を作る。ここでは3個作った。
• 下の画像では見えないが、上面の支間中央に同じように載荷線がある。

• ここまでで、載荷線や拘束線は完成。ただし、まだモデルにくっついていないため↓の操作でくっつける。
• 「操作」 -> 「パーティション」でオブジェクトにBoX_1を選び、
• ツールオブジェクトに拘束線や載荷線を選択
• 複数選択する場合はSHFTもしくはcntlを押しながら選択

• 適用して閉じる。

境界条件、荷重条件などを設定する箇所をグループ化

• オブジェクトブラウザーでさっき作成したPartition_1を右クリック -> グループを作成
• オブジェクトの種類 -> 線を選択
• グループの名前は、koteiとかわかりやすい名前に変更
• オブジェクトブラウザーからPartition_1階層下の拘束したい線(Line_1とか)を選択し、30とか40とかの数字が表示されたら「追加」をクリック -> 適用
• これを繰り返し必要な分グループを作成する。
• ここではグループを3個(kotei_1, kotei2, saika)を作った。
• 拘束条件などが同じ場合は1つのグループに複数の面や線を選択して構わない。

• その他のテクニック
• パーティションを作る際、高度なオプション（「グループとサブオブジェクトを継承します」「継承したサブオブジェクトの名前にプレフィックスを追加します」）の両方にチェックを入れておくと、自動でグループが作成されるんだったか、この辺に詳しい人、説明お願い。

「mesh」で要素分割

• 片持ち梁の「mesh」で要素分割と同じ。
• ただし、この単純梁のように拘束線や載荷線がある場合、直方体要素できれいに要素分割をするには追加で作業が必要になる場合がある。後々書いていく予定。

「AsterStudy?」で材料定数、境界条件、荷重条件などの設定

• 片持ち梁の「AsterStudy」とBC and Load以外は同じ。

入力ファイルを作る

• ツールバーの右上の方にあるAdd stage("{}+"こんなマーク)をクリック。すると、Data SettingsにStage_1が追加される。

• 以下の作業では、左側のData Settingsウィンドウ内のStage_1を選択したまま、ツールバーの"メッシュ△▽"みたいなのをクリックする。Stage_1を選択していないと、"△▽"は押せない。

各種設定の手順

• メッシュ -> Read a mesh
• Mesh file location : Mesh_1（Meshで作ったもの）
• Mesh file format : Med
• Model Definition -> Assign finite element
• At Least One : Finite element にチェック -> "+" -> Edit...
• Exactly One : Everywhere : Yes (1材料なら)
• Phenomenon -> Mechanic -> "+" -> 3D (立体要素なら)
• メッシュ : mesh(LIRE_MAILLAGE) (Read a mesh で付けた名前)
• Material -> Define material
• Linear isotropic elasticity にチェック(線形弾性体なら) -> Edit...
• Young's modulus : ヤング率を入力
• Poisson's ration : ポアソン比を入力
• Material -> Assign a material
• At Least One : Model にチェック -> model(AFFE_MODELE) (Assign finite element で付けた名前)
• Material assignement -> "+" -> Edit...
• Everywhere : Yes (1材料なら)
• Material -> "+" -> mater(DEFI_MATERIAL) (Define a material で付けた名前)
• BC and Load -> Assign mechanical load
• At Least One : Enforce DOF -> "+" -> Edit...
• At Least One : Group of element -> Edit... -> kotei1 (Geometryで付けた支点1のグループ名) -> 「OK」で閉じる。
• At Least One : DX・DY・DZにチェックを入れ、0.0を入力する。 -> 「OK」で閉じる。

• At Least One : Enforce DOF -> "+" -> Edit...
• At Least One : Group of element -> Edit... -> kotei2 (Geometryで付けた支点2のグループ名)
• At Least One : DX・DYにチェックを入れ、0.0を入力する。

• At Least One : FORCE_ARETE (線載荷なら) -> "+" -> Edit...
• At Least One : Group of element -> Edit... -> saika (Geometryで付けた載荷線のグループ名)
• At Least One : FY（載荷方向）にチェックを入れ、$\frac{荷重}{載荷線の長さ}$の値を入力する。
• Model : model(AFFE_MODELE) (Assign finite element で付けた名前)
• Analysis -> Static mechanical analysis
• At Least One : Material field : fieldmat(AFFE_MATERIAU) (Assign a material で付けた名前)
• Model : model(AFFE_MODELE) (Assign finite element で付けた名前)
• Loads -> "+" -> Edit...
• Load : load(AFFE_CHAR_MECA) (Assign mechanical load で付けた名前)
• Post Processing -> CALC_CHAMP
• Result : reslin(MECA_STATIQUE) (Static mechanical analysis で付けた名前)
• CONTRAINTE -> "+" -> "▽" -> SIGM_NOEU (節点の応力(6成分))
• DEFORMATION -> "+" -> "▽" -> EPSI_NOEU (節点のひずみ)
• CRETERES -> "+" -> "▽" -> SIEQ_NOEU (節点の応力(ミーゼスや主応力など))
• Output -> Set output results
• Result file location : 計算結果のファイル名(.med)と保存場所を指定する。
• Format : Med (ParaVis?で結果を見るのなら)
• Results -> "+" -> Edit...
• At Least One : Result : reslin(MECA_STATIQUE) (Static mechanical analysis で付けた名前)
• NOM_CHAM にチェック -> "+" -> "▽" -> DEPL (節点変位)
• Results -> "+" -> Edit...
• At Least One : Result : unnamed(CALC_CHAMP) (CALC_CHAMP で付けた名前)
• NOM_CHAM にチェック -> "+" -> "▽" -> EPSI_NOEU (節点のひずみ)
• NOM_CHAM にチェック -> "+" -> "▽" -> SIEQ_NOEU (節点の応力(ミーゼスや主応力など))
• NOM_CHAM にチェック -> "+" -> "▽" -> SIGM_NOEU (節点の応力(6成分))

• この通りに作ると、commファイルがこんな風になる。
• ちなみにこれは関塚さんの実験(PLA白の120x5x5の曲げ試験)に合わせて作ったものだが、結果がちゃんと合う。

計算を実行する

• 片持ち梁の計算を実行すると同じ。

「ParaVis?」で計算結果を確認

• 片持ち梁の「ParaVis」で計算結果を確認と同じ。

ASTK

• AsterStudy? -> Tools -> Plugins -> Salome-Meca -> Run astk

commファイルのバージョンを変換する（未解決）

~/salome_meca/V2017.0.2/tools/Salomemeca_translator-20170/bin/traduiV10V11 hoge.comm

Traceback (most recent call last):
  File "<string>", line 1, in <module>
ImportError?: No module named CodeAsterTranslator?.engine

#!/bin/bash

python -c 'from CodeAsterTranslator?.engine import traduitV10V11;traduitV10V11.main()' "$@"