#author("2021-04-27T13:55:01+09:00","default:kouzouken","kouzouken") #contents *salome-meca 入門編 片持ちばりから初等梁の線形解析 [#dd4ff6e4] -このページはsalome-mecaの基礎として、モデルの作成から解析結果までを順を追って説明するページとして作成していきます。 *単純梁 片持ちばりのhdfファイル [#s294a99e] https://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/oikawa/salomen/katamotitanzyyun.hdf *片持ちばりのcommファイル [#v4bf00c3] https://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/oikawa/salomen/katamoti.comm *単純梁のcommファイル [#p8355779] https://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/oikawa/salomen/tanzyun.comm **STEP1 片持ちばりの線形解析 [#ydb6fc47] -片持ちばりは以下のモデルを作成して線形解析を行う。 10×10×100の片持ちばり -http://www.str.ce.akita-u.ac.jp/~gotou/kouzou/png/kata1.png -材料諸量 LEFT:,ヤング率,7500N/mm^2 ,(木材ヤング率を採用) ,ポワソン比,0.4, *Geometryでモデルの作成を行う [#lc6f20ac] salomemecaを開くとモジュールの情報のところにSALOMEとなっている。そこのプルダウンメニューでいろんな機能に変更することができる。今回はその中のGeometryを選択して作業を始める。 Geometryを選択すると画面が変わり様々なアイコンが表示される画面に変更される。 その一番上にファイルや編集など表示してあるツールバーが表示されている。 その中の新しいエンティティをクリックすると基本オブジェクトや基本図形などが選べる画面になる。~ 今回はその中の基本図形からボックスを選択する。 ボックスの作成を選択するとボックスの作成という寸法を書き込む画面が表示される。 今回は10×10×100の片持ちばりなので、DX:10 DY:10 DZ:100にしてボックスを作成する。 適用して閉じるを押すと画面上に梁が作られている。 #ref https://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/oikawa/katamos1.png -次に拘束条件のためのグループの作成を行う。 オブジェクトブラウザーに先ほど作成したBOX1が表示されている。それを右クリックするとグループの作成が表示されるのでそれをクリック。 そしたらグループの作成という窓が立ち上がるので、そこでオブジェクトの種類で面を選択する。グループの名前は今回は適当にkoteiとかにする。選択すると画面に表示されているボックスの面が選択できるようになるのでとりあえず固定面にしたい面を選択して追加を押すと、画像のように33と選択された面の番号が表示される。表示されれば選択できているということなので、適用して閉じるで閉じる。オブジェクトブラウザーにkoteiというグループが作成されていることが確認できると思う。 今度はその要領で載荷側の面を作成する。 載荷側の面の作成が終わればGeometryでの作業は終了。 -https://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/oikawa/kotei2.png -https://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/oikawa/saica.png *meshの作成 [#of6c84be] メッシュに作成は解析精度に関わる重要な箇所になります。 メッシュの作成はまずツールバーのメッシュのメッシュの作成をします。を選択します。 そうするとメッシュの分割サイズを決める窓が出てくるのでそこで設定していきます。 写真のような感じです。 -https://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/oikawa/mesh1.png メッシュの分割設定が終わったらオブジェクトブラウザーのメッシュ1とかになっているものを右クリックするとメッシュの作成という歯車マークがあるのでそれをクリックすると計算が開始されます。 計算が終了すると計算が成功しました。というものが出てくるので出てきたら計算成功ということでOK! オブジェクトブラウザーの計算されたメッシュを右クリックするとメッシュに関する情報というのがあるのでクリックすると要素数が表示される。 ジオメトリーで作成したグループをメッシュに適用しないといけないと今後のAster Stadyの設定で設定できなくなるのでこの作業は必須です。 作成したメッシュを右クリックしてジオメトリーグループの作成を選択する。 選択するとジオメトリーとノードと書かれた窓が出てくるのでその中にジオメトリーで作成したグループを適用する。 ジオメトリーのグループの選択はコントロールで一括選択したほうがいいです。 次はAsterStadyです。 *「AsterStudy」で材料定数、境界条件、荷重条件などの設定 [#bdc9997f] -AsterStudyのアイコンをクリックしてAsterStudyを立ち上げる。 --プルダウンメニューからAsterStudyを選んで選択する。 http://www.str.ce.akita-u.ac.jp/~gotouhan/salome/gazou/s25aster.png **AsterStudyの各設定 [#y8ecc881] ,メッシュ,Read a mesh, メッシュの作成をしたモデルを適用します。, ,Model Definition,Assign finite element, モデルの要素を設定する。, ,Material ,Definr a material ,各材料物性値の設定,(複数設定可能) ,Material, Assign a material ,設定した物性値を材料に適用, ,BC and load , Assign mechabucal load , 拘束条件 荷重設定等, ,Analysis , static mechabucal analysis ,線形解析の設定 ,(ここはどのような解析を行うのかを設定する。) ,Post Processing , CALC_CHAMP ,ここは出力でなんの要素を出力するのかを設定する。, ,out put , set output result ,出力ファイルの設定を行う。, ***「AsterStudy」全般について [#u31cd70f] -右側のEdit comandウィンドウでは、最低限入力が必要なところには赤字で*が付いている。 -何かおかしな設定をしていたら、文字が赤く表示される。 -「適用」あるいは「OK」すると、左側下のInformationウィンドウにcommファイルに書かれる内容が表示される。 -後から設定を編集し直したいときは、"Data Settings"にある直したい項目をダブルクリックする。 -右の"?"を押すと公式ドキュメントが開かれる。ただし、フランス語で書かれているので、URLの …/fr/… を …/en/… にすると英語版になって読みやすいかも。 ***入力ファイルを作る [#k8b66ba4] -ツールバーの右上の方にあるAdd stage("{}+"こんなマーク)をクリック。すると、Data SettingsにStage_1が追加される。 http://www.str.ce.akita-u.ac.jp/~gotouhan/salome/gazou/s26kakko.png -以下の作業では、左側のData Settingsウィンドウ内のStage_1を選択したまま、ツールバーの"メッシュ△▽"みたいなのをクリックする。Stage_1を選択していないと、"△▽"は押せない。 http://www.str.ce.akita-u.ac.jp/~gotouhan/salome/gazou/s27sankaku.png ***各種設定の手順 [#xa79cf4d] AsterStudy上部にならんでいる、「メッシュ」「Model Definition」「Material」のプルダウンメニューを順番に選択し、それぞれの項目について、以下の手順で設定していく。 -メッシュ -> Read a mesh https://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/oikawa/salomen/katas1.png --Mesh file location->Mesh_1(Meshで作ったもの) --Mesh file format->Med ---20='0:1:2:4'の意味はよくわからない(違う数字になる場合もある) --OKをクリックすると、左側のData Settingウィンドウにmesh Read a meshができる。設定内容は、下のinformationウィンドウに表示される。 https://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/oikawa/salomen/katas2.png -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"が追加される。 https://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/oikawa/salomen/katas3.png -Material -> Define a material --Linear isotropic elasticity にチェック(線形等方弾性体なら) -> Edit... ---Young's modulus : ヤング率を入力 ---Poisson's ration : ポアソン比を入力 「OK」で閉じる。 https://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/oikawa/salomen/katas4.png -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」で全て閉じる。 https://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/oikawa/salomen/katas5.png -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... ---「OK」で閉じる。 --面載荷の場合は載荷する荷重を載荷する面積でわる。 https://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/oikawa/salomen/katas6.png -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」で閉じる。 https://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/oikawa/salomen/katas7.png -Post Processing -> CALC_CHAMP --Result : reslin(MECA_STATIQUE) (Static mechanical analysis で付けた名前)Salome-Meca2019では、以下3項目が書かれている箱の下 --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」で閉じる。 https://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/oikawa/salomen/katas8.png -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が現れる(チェックを入れないときのNOM_CMPと間違わないように) ---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 で付けた名前) ---NOM_CHAM にチェック -> "+" -> "▽" -> EPSI_NOEU (節点のひずみ) ---NOM_CHAM にチェック -> "+" -> "▽" -> SIEQ_NOEU (節点の応力(ミーゼスや主応力など)) ---NOM_CHAM にチェック -> "+" -> "▽" -> SIGM_NOEU (節点の応力(6成分)) ---「OK」で閉じる。 https://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/oikawa/salomen/katamoti.comm *計算結果を見る [#p5369738] *paraVISで変位等を見る方法 [#m70b7fe3] -パイプのブラウザーでファイルの選択(medファイル) -Apply -フィルター >Common >Warp By vector -アニメーションが見たいときは > Mechanics > Normal Modes animation > アニメーション概要(表示の中にある) *部材ごとに応力を見る方法 [#x7b6322d] -パイプのブラウザーでファイルの選択(medファイル)>Apply>フィルター >Alphabetical>Extract Group>見たい部材にチェックを入れる(2つあるので一応どっちにも入れる)>Apply>Spread sheet viewのShowingでExtract Groupを選択>select points throughで見たい部分を選択する>選んだポイントがピンク色で表示される>Spread sheet viewで選択した部分のデータが緑色に表示される *Timeが出ない場合 [#i9b1d9a3] -表示→ツールバー→current time controls *単純梁の線形解析 [#hdf5f9d5] https://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/oikawa/salomen/tans7.png -材料諸量 LEFT:,ヤング率,7500N/mm^2 ,(木材ヤング率を採用) ,ポワソン比,0.4, *Geometryでモデルの作成を行う [#lc6f20ac] salomemecaを開くとモジュールの情報のところにSALOMEとなっている。そこのプルダウンメニューでいろんな機能に変更することができる。今回はその中のGeometryを選択して作業を始める。 Geometryを選択すると画面が変わり様々なアイコンが表示される画面に変更される。 その一番上にファイルや編集など表示してあるツールバーが表示されている。 その中の新しいエンティティをクリックすると基本オブジェクトや基本図形などが選べる画面になる。~ -今回はその中の基本図形からボックスを選択する。 --ボックスの作成を選択するとボックスの作成という寸法を書き込む画面が表示される。 -今回は10×10×120の支間長100mmの単純ばりなので、DX:10 DY:10 DZ:120にしてボックスを作成する。 --適用して閉じるを押すと画面上に梁が作られている。 https://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/oikawa/salomen/tang1.png -今回は単純梁の解析を行うため拘束線を作成するため、ボックスに小さな孔をあける。 --そのための作業を今から行っていく。 -次に孔を空ける大きさのボックスを作成する。 --今回は高さ0.1mm,幅0.2mmの孔を空ける。 --なのでx=10,Y=0.1,Z=0.2に設定して作成する。 https://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/oikawa/salomen/tang2.png -作られたボックスを孔を空ける箇所に移動させる。 --今回の単純梁の支間は100mmとして解析を行うため ---両端部から10mmの箇所に拘束線の孔,支間長1/2の位置に載荷線の孔を空ける。 https://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/oikawa/salomen/tang3.png https://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/oikawa/salomen/tang4.png -孔を空ける箇所にボックスを設置したら次の方法で孔を空けていく。 --まずは、操作 > ブーリアン > カット を選択すると写真のような画面が出てくるので、 ---メインオブジェクトは作った梁を選択する。 ---ツールオブジェクトは切り抜きたい孔のサイズに作成したボックスを選択する。 ---選択が終了したら適用して閉じるを押して完了するとオブジェクトブラウザーにCUTというモデルが出来上がっている。 ---穴が開いたところに拘束線と載荷線を作成する。 ---作成したら https://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/oikawa/salomen/tang5.png -操作>パーティションを選択する。 --すると写真のような画面が登場するので --メインオブジェクトにカットした梁、ツールオブジェクトに作成した拘束線、載荷線を選択する。 --適用して閉じる。 ここまで来たら片持ち梁と同様な方法で拘束線と載荷線をグルーピングする。 https://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/oikawa/salomen/tang6.png そしたらジオメトリーの設定は終了。 メッシュの作成に移る。 といってもメッシュの作成は片持ち梁の時と一緒なので今回は割愛・・・ つなみに今回は孔を空けているため結構メッシュの計算に時間がかかると思う。 メッシュの作成が終了したら、ジオメトリーグループの作成を忘れないこと。 Aster Study の設定 基本的には片持ち梁の時と一緒です。 しかし載荷するのが今回は線なのでその設定だけ違う。 その点だけ説明する。 https://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/oikawa/salomen/tans1.png 違う点はBC and load の FORCE_FACE (面載荷) を FORCE_ARETE (線載荷)にする。 FORCE_ARETE (線載荷)を選択すると写真のような画面が出てくるので Group of elementに載荷線を選択する。 At Least One に載荷方向の荷重をかける。 荷重は載荷重/載荷線の長さで計算して入力する。 以上で設定は終了。 その他は片持ち梁と一緒。 以下の設定ファイルと同じになれば計算は完了する。 https://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/oikawa/salomen/tanzyun.comm