Salome-Meca演習_単純梁 の変更点

• 追加された行はこの色です。
• 削除された行はこの色です。
• Salome-Meca演習_単純梁 へ行く。
• Salome-Meca演習_単純梁 の差分を削除

#author("2022-04-20T17:44:54+09:00","default:kouzouken","kouzouken")
#author("2022-04-20T17:46:38+09:00","default:kouzouken","kouzouken")
#contents


*はじめに [#f2e42cdd]
10mm*10mm*120mmの直方体の梁（スパン100mm）を作って、両端から10mmの位置を''線拘束''にして，支間中央に荷重（線載荷）をかけます。

http://www.str.ce.akita-u.ac.jp/~gotouhan/salome/gazou/s62_tanzyun.png

-使用ツールについて、今回の説明ではSalome-Meca2019を使用しています。

*Salome-Mecaの起動→ジオメトリの作成 [#s6a7d01a]
**Salome-Meca2019を起動 [#sa84326e]
-左上の「新しいドキュメント」をクリック

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/geo1.png

**プルダウンメニューから「Geometry」を選択 [#y5ee9197]
-このような画面になる

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/geo2.png

**ボックスの作成 [#ra2cdf69]
-「ボックスを作成」をクリック

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/geo3.png

-Dx=100,Dy=10,Dz=20→適用して閉じる

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/geo4.png
-「Box_1」ができる

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/geo5.png

**拘束線、載荷線の作成 [#w9993dc6]
-拘束や載荷の箇所に「線」を作成します。
-先程の「ボックスを作成」をクリック
-Dx=10,Dy=0.1,Dz=0.2→適用して閉じる

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/geo6.png
-「Box_2」ができる

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/geo7.png
-Box_2を載荷線、拘束線の箇所まで移動します
-「移動」をクリック

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/geo8.png
-オブジェクト:BOX_2、Dx=0,Dy=4.9,Dz=10→適用して閉じる

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/geo9.png

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/geo10.png
-図のようになればok
-同じように、もう片方の拘束線と載荷線の位置にも作成

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/geo11.png
-「移動」をクリック
-オブジェクト:translation_1、Dx=0,Dy=0,Dz=50→適用

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/geo12.png
-オブジェクト:translation_2、Dx=0,Dy=0,Dz=50→適用して閉じる

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/geo13.png
-図のように3箇所にできていればok

-次に「Cut」を選択

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/geo14.png

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/geo15.png

-メインオブジェクト:Box_1、ツールオブジェクト:translation1,2,3を選択→適用して閉じる

-先程のtranslation1,2,3の箇所に穴が開く

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/geo13.png

-線を作成するための点を作成します
-「ポイントを作成」をクリック

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/geo17.png

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/geo16.png
-Dx=10,Dy=4.9,Dz=10.1→適用して閉じる

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/geo18.png

-図のように点ができる
-反対側も作成します。「移動」をクリック
-Dx=-10,Dy=0,Dz=0→適用して閉じる
-反対側も作成したら、「ラインを作成」をクリック

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/geo19.png

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/geo20.png

-点1と点2に「Veltex_1」と「translation_4（移動した点）」を選択
-穴の中心に線ができる

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/geo21.png

-同じように、もう片方の拘束線と載荷線の位置にも作成

-「移動」をクリック
-オブジェクト:Line_1、Dx=0,Dy=0,Dz=50→適用
-オブジェクト:translation_5、Dx=0,Dy=0,Dz=50→適用して閉じる

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/geo13.png
-図の○の位置を拡大すると、

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/geo21.png
-このようになっていればok

**拘束線、載荷線をパーティション→グループの作成 [#w772bc4a]
-線と梁を一体化させます
-「パーティション」を選択

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/geo22.png

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/geo23.png
-メインオブジェクト:Cut_1、ツールオブジェクト:Line_1,translation5,6を選択→適用して閉じる
-「Partition_1」ができる
-拘束線、載荷線についてのグループを作成します
-Partition_1を右クリック→「グループを作成」をクリック

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/geo24.png
-オブジェクトの種類:「/（線）」を選択
-設定した拘束線を選択し「追加」をクリック
-名前をわかりやすく「kotei1」とかにしておく→適用して閉じる
-同じようにもう一つの拘束線、載荷線について設定する。名前は「kotei2」「saika」とかにする

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/geo25.png
-図のように3本の線がpartition_1のグループになっていればok
-ジオメトリは以上
*メッシュの作成 [#h194187d]
**メッシュの切り方の設定 [#b00a214b]
-プルダウンメニューから「Mesh」を選択
-オブジェクトブラウザーで、上で作成した「Partition_1」を選択→メッシュを作成（左側にある方）を選択

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/mes1.png

-3D→アルゴリズム「NETGEN 3D」を選択→詳細設定セットの割当て→「3D Automatic Tetrahedralization」をクリック

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/mes2.png

-要素長設定（メッシュの細かさの設定）をする。モデルに対してあまりにも大きいと思ったような解析結果にならず、逆に細かいと計算に時間がかかったりエラーが出たりする。
-今回はDx=10,Dy=14,Dz=100のモデルのため、最大要素長=2とする。長さのところに「2」を入力→ok

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/mes3.png

**メッシュを切る [#z6670bd1]
-オブジェクトブラウザーで、上で作成した「Mesh_1」を選択→メッシュを作成（右側にある方）を選択

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/mes4.png

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/mes5.png

-計算が成功すると、このような形になる

**グループの作成 [#fe418fe4]
-ジオメトリで作成したグループをメッシュに反映させる
-「Mesh_1」で右クリック→「ジオメトリのグループを作成」をクリック

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/mes6.png

-Geometryで作成したPartition_1のグループ「kotei1」「kotei2」「saika」を選択→適用して閉じる

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/mes7.png

-オブジェクトブラウザーに反映され画像のようになっていれば成功
-メッシュの作成は以上
*AsterStudyの設定 [#e586cda1]
-プルダウンメニューから「AsterStudy」を選択
-上のメニューから「Add Stage」を選択

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/ast1.png

**メッシュの読み込み [#keb2c836]
-メッシュファイルを読み込む設定をします
-「メッシュ」→「Read a Mesh」を選択

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/ast2.png

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/ast3.png

-Mesh file location:Mesh_1（作成したメッシュ）、Mesh file format:Med を選択してok
-左下のInfomationにこのように表示されればok

  mesh = LIRE_MAILLAGE(
  FORMAT='MED', 
  UNITE=20

**Model Definition（解析方法について） [#yf61b8cf]
-解析の方法の設定をします
-「Model Definition」→「Assign finite element」を選択

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/ast4.png

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/ast5.png

-「メッシュ」が「mesh(LIRE_MAILLAGE)」になっているか確認
-「Finite element assignement」を1itemに増やす→「Edit」
-Everywhere:Yes　Phenomenon:Mechanic　Modelisation:3Dを選択→ok


http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/ast6.png

-左下のInfomationにこのように表示されればok
 model = AFFE_MODELE(
  AFFE=_F(
    MODELISATION=('3D', ), 
    PHENOMENE='MECANIQUE', 
    TOUT='OUI'
  ), 
  MAILLAGE=mesh
 ) 

**Material（材料の設定） [#s18e3271]
-ここで材料の設定をします。
-今回は鋼材の設定をします。

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/ast7.png
-上のメニューの「Material」から「Define a material」を選択

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/ast8.png

-「Linear isotropic elastcity」を選択（一般的な材料であればこれを選択）

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/ast9.png

-Young's modulus=ヤング率
-Puisson's ratio=ポアソン比
-それぞれ入力→ok→わかりやすい名前をつける（ここでは「kou」）→ok
-ここでは鋼材をヤング率=20.0kN/mm2、ポアソン比=0.3とする

 kou = DEFI_MATERIAU(
  ELAS=_F(
    E=200000.0, 
    NU=0.3
  )
 )

-上のメニューの「Material」から「Assign a material」を選択

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/ast10.png

-「Model」にチェック→Model=model(AFFE_MODELE)になってることを確認
-「Material assignment」を追加→編集

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/ast11.png

-「Everywhere」→Yes(全体が1つの材料のため)
-「Material」で先程Define materialで設定した鋼材の材料情報(kou)を選択→ok
-「Material assignment」を追加→ok

 fieldmat = AFFE_MATERIAU(
  AFFE=_F(
    MATER=(kou, ), 
    TOUT='OUI'
  ), 
  MODELE=model
 ) 


**BC and Load （荷重と固定について）[#ia3e6be8]
-荷重と固定についての設定
-上のメニューから「BC and Load」→「Assign mechanical load」を選択

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/ast12.png

***Enforce DOF （固定）[#x5566b9b]
-固定についての設定（Enforce DOF=強制変位　ここでは強制変位の設定を0にする=固定ということにしている）
-「Enforce DOF」を選択

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/ast13.png


-「Group of element」でジオメトリで設定した固定面（片方）を選択（ここでは「kotei1」）
-完全固定支承にするので、DX,DY,DZ=0を入力→ok
-もう一つ「Enforce DOF」を選択（2itemにする）
-「Group of element」でもう片方を選択（ここでは「kotei2」）
-ローラー支承にするので、DX,DY=0を入力→ok

***FORCE_ARETE(線載荷) [#yc80519f]
-荷重についての設定（FORCE_ARETE=線載荷　線に荷重をかける場合これを選択）
-「FORCE_ARETE」を選択(Enforce DOFよりもかなり下にある)

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/ast15.png

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/ast16.png

-「Group of element」でジオメトリで設定した載荷線を選択（ここでは「saika」）
-FY（載荷方向）にチェックを入れ、荷重/載荷線の長さ　を入力→ok
-ここでは載荷線の長さ=10mmなので単位荷重として10を入力

 load = AFFE_CHAR_MECA(
  DDL_IMPO=_F(
    DX=0.0, 
    DY=0.0, 
    DZ=0.0, 
    GROUP_MA=('kotei1', 'kotei2')
  ), 
  FORCE_ARETE=_F(
    FY=10.0, 
    GROUP_MA=('saika', )
  ), 
  MODELE=model
 )
-荷重の設定は以上

**Analysis（設定をまとめる） [#e3453b62]
-今までのモデルの設定、材料の設定、荷重と固定の設定を紐付けます。
***Static mechanical analysys（線形解析の設定） [#p40099e7]
http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/ast18.png

-上のメニューから「Analysis」→「Static mechanical analysys」を選択

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/ast19.png

-Material field:fieldmat(AFFE_MATERIAU) Model:model(AFFE_MODELE)になっていることを確認
-「Loads」で先程BC and loadで設定したloadを設定→ok

 reslin = MECA_STATIQUE(
  CHAM_MATER=fieldmat, 
  EXCIT=_F(
    CHARGE=load
  ), 
  MODELE=model
 )

**Post Processing（どの数値を計算するか） [#n1cad0ee]
-
***CALC_CHAMP [#a6c20699]
http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/ast20.png

-上のメニューから「Post Processing」→「CALC_CHAMP」を選択

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/ast21.png

-CONTRAITE:SIGM_NOEU(節点の応力(6成分))を選択
-DEFORMATION:EPSI_NOEU (節点のひずみ)を選択
-CRETERES:SIEQ_NOEU (節点の応力(ミーゼス応力や主応力など)) を選択
-FORCE:REAC_NODA(節点の反力)を選択
-補足：ちなみにNOEUは節点(NOEUd(フランス語))、ELGAはガウスの積分点(ELement GAuss point)、ELNOは要素(のどこか?)(ELement Node)だとかで、ParaVisでコンター図を描きたいときは Point Data であるNOEUを出力する必要がある。ELGAやELNOは、ParaVis上でSpredSheetViewを開き、Field Data や Cell Data を選択すれば値が見れる。

-Result:reslin(MECA_STATIQUE)になってるのを確認→ok

 unnamed = CALC_CHAMP(
  CONTRAINTE=('SIGM_NOEU', ), 
  CRITERES=('SIEQ_NOEU', ), 
  DEFORMATION=('EPSI_NOEU', ), 
  FORCE=('REAC_NODA', ), 
  RESULTAT=reslin
 ) 
**Output（出力設定） [#q7534a47]
-出力設定をします。

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/ast22.png

-上のメニューから「Output」→「Set output results」を選択

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/ast23.png

-Result file location:パソコン上の任意の場所を選択（そこにparavisを見るための.medファイルが保存される）
-Resultsを追加→編集

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/ast24.png

-Result:reslin(MECA_STATIQUE)を選択→「NOM_CHAM」という欄が出る
-NOM_CHAM:DEPL(変位)を選択→ok

-もう一つResultsを追加→編集

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/ast25.png

-同じようにResult:unnamed(CALC_CHAMP)を選択→「NOM_CHAM」が出る
-NOM_CHAM:CALC_CHAMPで選択した4つの要素(SIGM_NOEU,EPSI_NOEU,SIEQ_NOEU,REAC_NODA)を選択→ok

-Format:Medを選択→ok

 IMPR_RESU(
  FORMAT='MED', 
  RESU=(_F(
      NOM_CHAM=('DEPL', ), 
      RESULTAT=reslin
    ), _F(
      NOM_CHAM=('EPSI_NOEU', 'SIGM_NOEU', 'SIEQ_NOEU', 'REAC_NODA'), 
      RESULTAT=unnamed
    )), 
  UNITE=80
 ) 
**History Viewで計算 [#z80fd56a]
-ここまで来たら、一度.hdfファイル（salomeの設定ファイル）を保存しておく。

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/ast26.png

-左上の「History View」をクリック

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/ast27.png

-Casesの中の「CurrentCase」をクリック→Stage_1の横の+をクリック→下のほうにある「Run」をクリック→計算開始
-Auto Refresh:5s　にすると計算の経過を見れる（5秒ごとに更新してくれる）
-成功すると緑色、失敗すると赤になる。
-補足：Salome-Meca2017(だったと思う)以前では失敗すると赤色か黄色になるとのことです。そのときの赤色は深刻なエラー（ソフト構成上の問題とか）、黄色は計算エラー（パラメータのミス等）らしい。
*Paravisで計算結果の確認 [#ae45333f]
-プルダウンメニューから「Paravis」を選択
-左側のパイプのブラウザーにある「builtin:」を右クリック→「open」を選択→先ほどの計算結果（.medファイル）を開く
-計算結果が出てくる

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/par1.png

-初期設定で「Solid Color」となっているところをプルダウンで「reslin_DEPL（変位）」や「unnamed_EPSI_NOEU（ひずみ）」等に変えると、任意のパラメータを確認することができる

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/par2.png

**変形の確認 [#h4cecf9d]
-任意のパラメータに変更（ここではreslin_DEPL）
-オブジェクトインスペクターの一番上にある「適用」を押す

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/par3.png
-上のメニューの「フィルター」から、「common」→「Warp By Vector」を選択
-パイプのブラウザーに「WarpByVector1」が追加される→選んでオブジェクトインスペクターの一番上にある「適用」を押す
-構造物の変形が確認できる

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/par4.png

**数値の確認 [#l68413c6]
-右上の「RenderView1」の横にあるボタンを押す（どちらでもよい）

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/par5.png

-一番下の「SpreadSheet View」を選択
-メッシュを切った際のノード（節点）ごとの各パラメータの数値が表示される

***任意の点の数値を確認したい場合 [#dddd2641]
-Select Points ○○ (変位を見たい節点を選択する) 

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/par6.png
- （上の方にある）Show only selected Elements. 

http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/par7.png