salome-meca 振動解析メモ
をテンプレートにして作成
[
トップ
] [
新規
|
一覧
|
検索
|
最終更新
|
ヘルプ
|
ログイン
]
開始行:
#contents
*密度を用いる数値解析（固有振動数・振動解析など）に用いる...
-密度の単位については[[構造関係メモ:https://www.str.ce.ak...
'''有限要素法ソフトが出力する応力の単位は[Pa]，変位の単位...
'''形状モデルを[mm]単位でモデリングした場合を述べます。ヤ...
'''有限要素法ソフトが出力する応力の単位は[N/mm2]ないしは[...
'''ここで密度の単位[ton/mm3]が気になります。なぜ[kg/mm3]...
**木材と鋼材を用いたモデルをmm-N系でモデリングする場合 [#...
,,単位,木材,鋼材
,ヤング率,MPa(N/mm$^2$),6000,206000
,ポアソン比,--,0.4,0.3
,密度,ton/mm$^3$,3.8$\times10^{-10}$,7.8$\times10^{-9}$
,重力加速度,mm/s$^2$,9800,9800
木材とかは樹種によってヤング率とか密度とか変わってくるの...
今回はスギ材を仮定しています．マツ系とか$5.0\times10^{-10...
なお，[m]単位でモデリングするときは，密度の単位はkg/m$^3$...
固有値解析（モード解析）をするときはm系でモデリングするこ...
単位系がわかりやすいから！！
*振動解析 [#zdff5a62]
**固有振動数の理論値 [#lcf54dec]
-[[日本機械学会構造振動学:https://www.jsme.or.jp/sed/guid...
**ソリッド要素の振動解析 [#j0a3e91b]
***[[modal.hdf:http://www.str.ce.akita-u.ac.jp/~gotouhan/...
-[[Salomemeca使用解説法:http://opencae.gifu-nct.ac.jp/puk...
**梁要素の片持ち梁の解析 [#mc1629ff]
***[[kukei1125.hdf:http://www.str.ce.akita-u.ac.jp/~gotou...
-梁要素の片持ち梁の解析．[[Salome-Meca2017]]とやり方はほ...
--ModelDefinition,Assign finante elementで，Modalisation...
にする
--ModelDefinition,AFFE_CARA_ELEMを追加し，
--1)POITREにチェックを入れ，Edit
--2)Group of elementは，梁要素につけたグループ名．(このhd...
--3)SECTIONでは，RECTANGLE(矩形という意味)，CARAはHY,HZ，...
---HY,HZについては，[[ここ:https://www.code-aster.org/V2/...
--BC&LoadのEnforce DOFでDX,DY,DZに加え，DRX,DRY,DRZにも0...
--載荷は点なので，FORCE_NODALE
--Post Processingは今回は全部いらない
--Analysis,At Least oneはどちらもチェック入れる
--outputの出力はResult,NOM_CHAMのDEPLのみ
-Paravisにおいては，WarpByVectorの（私のパソコンでは左に...
**梁要素の片持ち梁の振動解析 [#ufc0332e]
***[[kukeimode1204.hdf:http://www.str.ce.akita-u.ac.jp/~g...
-梁の片持ち梁のモード解析
--Analysis,CALC_MODES，Optionでは[[このシェル要素のYoutub...
--Analysis,CALC_MODESのSolverとSTOP_BANDEとTYPE_RESUはぜ...
--AnalysisのVERI_MODEをEditし，STOP_ERREURのチェックを外...
**ソリッド要素のラーメン構造の振動解析 [#x840b790]
***[[ramenbox1206.hdf:http://www.str.ce.akita-u.ac.jp/~go...
-メッシュを切るときに，詳細セットの割り当て，3D:Automatic...
-AsterStudyは，立体要素の振動解析と同じ設定で解ける．
**梁要素のラーメン構造の振動解析 [#o66b5917]
***[[ramenhari1209.hdf:http://www.str.ce.akita-u.ac.jp/~g...
-ジオメトリでは，全部の梁要素を選択して一つのジオメトリの...
-メッシュを切る際は，Maxsizeではなく，Number of Segment
-AsterStudyは，梁の片持ち梁の振動解析と同じ
**シェル要素の片持ちばりの振動解析 [#p1d749e5]
***[[katamoti2d.hdf:http://www.str.ce.akita-u.ac.jp/~goto...
-梁要素との違いだけ記載
-Aster Studyの設定
-Model Definition -Assign finite element -MOdelisation -1...
-AFFE_CARA_ELEM -COQUEを追加
-COQUEの設定 -Exactly one -Group of element (グループの...
-INER_ROTA (intertia of rotationを計算する）-yes
あとの設定は一緒です。
-Analysis-CALC MODES-オプション-BANDE（今はPLUS-PETITEに...
**梁要素とシェル要素の振動解析 [#x0f40d4e]
-[[harishell.comm:http://www.str.ce.akita-u.ac.jp/~gotouh...
-[[harishell.hdf :http://www.str.ce.akita-u.ac.jp/~gotouh...
-梁要素とシェル要素を結合させるにはジオメトリー上でコンパ...
メッシュを切ってAster studyの設定をする。
**梁とシェルの結合による骨組モデルの振動解析（LIAISON_SOL...
***[[harishellketa.hdf:http://www.str.ce.akita-u.ac.jp/~g...
***[[harishellketa.comm:http://www.str.ce.akita-u.ac.jp/~...
-今回はLIAISON_SOLIDEの結合コマンドによる梁とシェルの結合...
LIAISON_SOLIDEは節点同士を変形のない剛体として結合する方...
-ジオメトリーで行うこと。シェルの節点と同位置にある桁の節...
--事前に桁の節点とシェルの節点をパーティションで同一モデ...
-メッシュはあまりやることはありません。（いつも通りってこ...
-Aster Study ではPOU_D_EとCOQUEの設定は通常通り設定してく...
-LIAISON_SOLIDEはBC and Load の中にあります。選択して節点...
DEBUT(LANG='EN')
mesh = LIRE_MAILLAGE(FORMAT='MED',
UNITE=2)
model = AFFE_MODELE(AFFE=(_F(MODELISATION=('3D', ),
PHENOMENE='MECANIQUE',
TOUT='OUI'),
_F(GROUP_MA=('ta4', 'ta5', 'ta6...
MODELISATION=('POU_D_E', ),
PHENOMENE='MECANIQUE'),
_F(GROUP_MA=('sss', ),
MODELISATION=('DKT', ),
PHENOMENE='MECANIQUE')),
MAILLAGE=mesh)
elemprop = AFFE_CARA_ELEM(COQUE=_F(ANGL_REP=(0.0, 0.0),
COQUE_NCOU=1,
EPAIS=8.0,
GROUP_MA=('sss', ),
INER_ROTA='OUI'),
MODELE=model,
POUTRE=_F(CARA=('A', 'IY', 'IZ'...
GROUP_MA=('ta4', 'ta5...
SECTION='GENERALE',
VALE=(100.0, 833.0, 8...
VARI_SECT='CONSTANT'))
mater = DEFI_MATERIAU(ELAS=_F(E=7500.0,
NU=0.4,
RHO=0.0459))
fieldmat = AFFE_MATERIAU(AFFE=_F(MATER=(mater, ),
TOUT='OUI'),
MODELE=model)
load = AFFE_CHAR_MECA(DDL_IMPO=_F(DX=0.0,
DY=0.0,
DZ=0.0,
GROUP_MA=('kotei', )),
LIAISON_SOLIDE=(_F(GROUP_NO=('ke1',...
_F(GROUP_NO=('ke2',...
_F(GROUP_NO=('ke3',...
_F(GROUP_NO=('ke4',...
_F(GROUP_NO=('ke5',...
_F(GROUP_NO=('ke6',...
_F(GROUP_NO=('ke7',...
_F(GROUP_NO=('ke8',...
MODELE=model)
ASSEMBLAGE(CARA_ELEM=elemprop,
CHAM_MATER=fieldmat,
CHARGE=(load, ),
MATR_ASSE=(_F(MATRICE=CO('mass'),
OPTION='MASS_MECA'),
_F(MATRICE=CO('stifness'),
OPTION='RIGI_MECA')),
MODELE=model,
NUME_DDL=CO('nddl'))
modes = CALC_MODES(CALC_FREQ=_F(NMAX_FREQ=10),
MATR_MASS=mass,
MATR_RIGI=stifness,
OPTION='PLUS_PETITE',
VERI_MODE=_F(STOP_ERREUR='NON'))
IMPR_RESU(FORMAT='MED',
RESU=_F(NOM_CHAM=('DEPL', ),
RESULTAT=modes),
UNITE=3)
FIN()
**かじか橋 床版シェル要素，他梁要素 振動解析 [#me5deea4]
***[[arch20182dkeisan0716epic.hdf:http://www.str.ce.akita...
***シェル要素の異方性 [#f8fcc1d2]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2019/aoyama/de...
シェル要素の異方性はDEFI_COMPOSITEの設定が必要(?)
DEFI_COMPOSITEのデフォルトの名前だとエラーになることがあ...
そのため名前を新たに設定することを推奨．（画像では「aaa」）
**断面性能の計算 [#fc178057]
-梁要素の解析を行うためには断面諸量を与えてあげる必要があ...
,A,断面積
,IY,断面二次モーメント
,IZ,断面二次モーメント
,EYEZ,ねじりの中心位置
,JX,断面二次極モーメント
断面計算は計算する断面のシェルをジオメトリーで作成し（グ...
ASTKで断面性能を出力する。salome-meca2018でASTKを設定する...
-ASTK[[https://sites.google.com/site/codeastersalomemeca/...
-断面性能の計算[[https://sites.google.com/site/codeasters...
終了行:
#contents
*密度を用いる数値解析（固有振動数・振動解析など）に用いる...
-密度の単位については[[構造関係メモ:https://www.str.ce.ak...
'''有限要素法ソフトが出力する応力の単位は[Pa]，変位の単位...
'''形状モデルを[mm]単位でモデリングした場合を述べます。ヤ...
'''有限要素法ソフトが出力する応力の単位は[N/mm2]ないしは[...
'''ここで密度の単位[ton/mm3]が気になります。なぜ[kg/mm3]...
**木材と鋼材を用いたモデルをmm-N系でモデリングする場合 [#...
,,単位,木材,鋼材
,ヤング率,MPa(N/mm$^2$),6000,206000
,ポアソン比,--,0.4,0.3
,密度,ton/mm$^3$,3.8$\times10^{-10}$,7.8$\times10^{-9}$
,重力加速度,mm/s$^2$,9800,9800
木材とかは樹種によってヤング率とか密度とか変わってくるの...
今回はスギ材を仮定しています．マツ系とか$5.0\times10^{-10...
なお，[m]単位でモデリングするときは，密度の単位はkg/m$^3$...
固有値解析（モード解析）をするときはm系でモデリングするこ...
単位系がわかりやすいから！！
*振動解析 [#zdff5a62]
**固有振動数の理論値 [#lcf54dec]
-[[日本機械学会構造振動学:https://www.jsme.or.jp/sed/guid...
**ソリッド要素の振動解析 [#j0a3e91b]
***[[modal.hdf:http://www.str.ce.akita-u.ac.jp/~gotouhan/...
-[[Salomemeca使用解説法:http://opencae.gifu-nct.ac.jp/puk...
**梁要素の片持ち梁の解析 [#mc1629ff]
***[[kukei1125.hdf:http://www.str.ce.akita-u.ac.jp/~gotou...
-梁要素の片持ち梁の解析．[[Salome-Meca2017]]とやり方はほ...
--ModelDefinition,Assign finante elementで，Modalisation...
にする
--ModelDefinition,AFFE_CARA_ELEMを追加し，
--1)POITREにチェックを入れ，Edit
--2)Group of elementは，梁要素につけたグループ名．(このhd...
--3)SECTIONでは，RECTANGLE(矩形という意味)，CARAはHY,HZ，...
---HY,HZについては，[[ここ:https://www.code-aster.org/V2/...
--BC&LoadのEnforce DOFでDX,DY,DZに加え，DRX,DRY,DRZにも0...
--載荷は点なので，FORCE_NODALE
--Post Processingは今回は全部いらない
--Analysis,At Least oneはどちらもチェック入れる
--outputの出力はResult,NOM_CHAMのDEPLのみ
-Paravisにおいては，WarpByVectorの（私のパソコンでは左に...
**梁要素の片持ち梁の振動解析 [#ufc0332e]
***[[kukeimode1204.hdf:http://www.str.ce.akita-u.ac.jp/~g...
-梁の片持ち梁のモード解析
--Analysis,CALC_MODES，Optionでは[[このシェル要素のYoutub...
--Analysis,CALC_MODESのSolverとSTOP_BANDEとTYPE_RESUはぜ...
--AnalysisのVERI_MODEをEditし，STOP_ERREURのチェックを外...
**ソリッド要素のラーメン構造の振動解析 [#x840b790]
***[[ramenbox1206.hdf:http://www.str.ce.akita-u.ac.jp/~go...
-メッシュを切るときに，詳細セットの割り当て，3D:Automatic...
-AsterStudyは，立体要素の振動解析と同じ設定で解ける．
**梁要素のラーメン構造の振動解析 [#o66b5917]
***[[ramenhari1209.hdf:http://www.str.ce.akita-u.ac.jp/~g...
-ジオメトリでは，全部の梁要素を選択して一つのジオメトリの...
-メッシュを切る際は，Maxsizeではなく，Number of Segment
-AsterStudyは，梁の片持ち梁の振動解析と同じ
**シェル要素の片持ちばりの振動解析 [#p1d749e5]
***[[katamoti2d.hdf:http://www.str.ce.akita-u.ac.jp/~goto...
-梁要素との違いだけ記載
-Aster Studyの設定
-Model Definition -Assign finite element -MOdelisation -1...
-AFFE_CARA_ELEM -COQUEを追加
-COQUEの設定 -Exactly one -Group of element (グループの...
-INER_ROTA (intertia of rotationを計算する）-yes
あとの設定は一緒です。
-Analysis-CALC MODES-オプション-BANDE（今はPLUS-PETITEに...
**梁要素とシェル要素の振動解析 [#x0f40d4e]
-[[harishell.comm:http://www.str.ce.akita-u.ac.jp/~gotouh...
-[[harishell.hdf :http://www.str.ce.akita-u.ac.jp/~gotouh...
-梁要素とシェル要素を結合させるにはジオメトリー上でコンパ...
メッシュを切ってAster studyの設定をする。
**梁とシェルの結合による骨組モデルの振動解析（LIAISON_SOL...
***[[harishellketa.hdf:http://www.str.ce.akita-u.ac.jp/~g...
***[[harishellketa.comm:http://www.str.ce.akita-u.ac.jp/~...
-今回はLIAISON_SOLIDEの結合コマンドによる梁とシェルの結合...
LIAISON_SOLIDEは節点同士を変形のない剛体として結合する方...
-ジオメトリーで行うこと。シェルの節点と同位置にある桁の節...
--事前に桁の節点とシェルの節点をパーティションで同一モデ...
-メッシュはあまりやることはありません。（いつも通りってこ...
-Aster Study ではPOU_D_EとCOQUEの設定は通常通り設定してく...
-LIAISON_SOLIDEはBC and Load の中にあります。選択して節点...
DEBUT(LANG='EN')
mesh = LIRE_MAILLAGE(FORMAT='MED',
UNITE=2)
model = AFFE_MODELE(AFFE=(_F(MODELISATION=('3D', ),
PHENOMENE='MECANIQUE',
TOUT='OUI'),
_F(GROUP_MA=('ta4', 'ta5', 'ta6...
MODELISATION=('POU_D_E', ),
PHENOMENE='MECANIQUE'),
_F(GROUP_MA=('sss', ),
MODELISATION=('DKT', ),
PHENOMENE='MECANIQUE')),
MAILLAGE=mesh)
elemprop = AFFE_CARA_ELEM(COQUE=_F(ANGL_REP=(0.0, 0.0),
COQUE_NCOU=1,
EPAIS=8.0,
GROUP_MA=('sss', ),
INER_ROTA='OUI'),
MODELE=model,
POUTRE=_F(CARA=('A', 'IY', 'IZ'...
GROUP_MA=('ta4', 'ta5...
SECTION='GENERALE',
VALE=(100.0, 833.0, 8...
VARI_SECT='CONSTANT'))
mater = DEFI_MATERIAU(ELAS=_F(E=7500.0,
NU=0.4,
RHO=0.0459))
fieldmat = AFFE_MATERIAU(AFFE=_F(MATER=(mater, ),
TOUT='OUI'),
MODELE=model)
load = AFFE_CHAR_MECA(DDL_IMPO=_F(DX=0.0,
DY=0.0,
DZ=0.0,
GROUP_MA=('kotei', )),
LIAISON_SOLIDE=(_F(GROUP_NO=('ke1',...
_F(GROUP_NO=('ke2',...
_F(GROUP_NO=('ke3',...
_F(GROUP_NO=('ke4',...
_F(GROUP_NO=('ke5',...
_F(GROUP_NO=('ke6',...
_F(GROUP_NO=('ke7',...
_F(GROUP_NO=('ke8',...
MODELE=model)
ASSEMBLAGE(CARA_ELEM=elemprop,
CHAM_MATER=fieldmat,
CHARGE=(load, ),
MATR_ASSE=(_F(MATRICE=CO('mass'),
OPTION='MASS_MECA'),
_F(MATRICE=CO('stifness'),
OPTION='RIGI_MECA')),
MODELE=model,
NUME_DDL=CO('nddl'))
modes = CALC_MODES(CALC_FREQ=_F(NMAX_FREQ=10),
MATR_MASS=mass,
MATR_RIGI=stifness,
OPTION='PLUS_PETITE',
VERI_MODE=_F(STOP_ERREUR='NON'))
IMPR_RESU(FORMAT='MED',
RESU=_F(NOM_CHAM=('DEPL', ),
RESULTAT=modes),
UNITE=3)
FIN()
**かじか橋 床版シェル要素，他梁要素 振動解析 [#me5deea4]
***[[arch20182dkeisan0716epic.hdf:http://www.str.ce.akita...
***シェル要素の異方性 [#f8fcc1d2]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2019/aoyama/de...
シェル要素の異方性はDEFI_COMPOSITEの設定が必要(?)
DEFI_COMPOSITEのデフォルトの名前だとエラーになることがあ...
そのため名前を新たに設定することを推奨．（画像では「aaa」）
**断面性能の計算 [#fc178057]
-梁要素の解析を行うためには断面諸量を与えてあげる必要があ...
,A,断面積
,IY,断面二次モーメント
,IZ,断面二次モーメント
,EYEZ,ねじりの中心位置
,JX,断面二次極モーメント
断面計算は計算する断面のシェルをジオメトリーで作成し（グ...
ASTKで断面性能を出力する。salome-meca2018でASTKを設定する...
-ASTK[[https://sites.google.com/site/codeastersalomemeca/...
-断面性能の計算[[https://sites.google.com/site/codeasters...
ページ名:
![[PukiWiki]](image/pukiwiki.png "[PukiWiki]")
