- 追加された行はこの色です。
- 削除された行はこの色です。
#author("2021-06-21T10:24:46+09:00","default:kouzouken","kouzouken")
#author("2021-06-21T12:00:58+09:00","default:kouzouken","kouzouken")
[[FrontPage]]
#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
http://www.str.ce.akita-u.ac.jp/~gotouhan/ishikuro/tanzyun/ast14.png
-「Group of element」でジオメトリで設定した固定面を選択(ここでは「kotei1,2」)
-DX,DY,DZ=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(だったと思う)以前では失敗すると赤色か黄色になるとのことです。そのときの赤色は深刻なエラー(ソフト構成上の問題とか)、黄色は計算エラー(パラメータのミス等)らしい。