- 追加された行はこの色です。
- 削除された行はこの色です。
#author("2022-09-14T16:15:56+09:00;2022-08-26T17:05:19+09:00","default:kouzouken","kouzouken")
#author("2022-09-21T11:17:03+09:00;2022-08-26T17:05:19+09:00","default:kouzouken","kouzouken")
#contents
*修論 [#a6e5beaf]
テーマ:ケーブル腐食を考慮した斜張橋の地震荷重に関する数値解析的検討
*地震 [#vddf00a4]
**現段階のメモ[#w21fb467]
○2週間インターンシップに参加するの研究ができないのでまとめておく~
○地震解析~
解析モデル設定→固有値解析→減衰特性の設定→地震動(加速度?)を与える→地震動に対する応答解析~
上記の流れをsalomeで行う~
※これを過渡応答解析(transitory analysis)というらしい。~
100☓100☓1000の箱桁をシェル要素で減衰無しで地震をかける(2ヶ月苦戦中)~
以下現段階のsalomeの設定(変更点)~
ASSEMBLAGE(
CARA_ELEM=elemprop,
CHAM_MATER=fieldmat,
CHARGE=(load0, ),
MATR_ASSE=(_F(
MATRICE=CO('mass'),
OPTION='MASS_MECA'
), _F(
MATRICE=CO('Rigi'),
OPTION='RIGI_MECA'
)),
MODELE=model,
NUME_DDL=CO('nddl'),
VECT_ASSE=_F(
CHARGE=(load0, ),
OPTION='CHAR_MECA',
VECTEUR=CO('Vect')
)
)
modes0 = CALC_MODES(
CALC_FREQ=_F(
NMAX_FREQ=10
),
MATR_MASS=mass,
MATR_RIGI=Rigi,
OPTION='PLUS_PETITE',
VERI_MODE=_F(
STOP_ERREUR='NON'
)
)
PROJ_BASE(
BASE=modes0,
MATR_ASSE_GENE=(_F(
MATRICE=CO('massgene'),
MATR_ASSE=mass
), _F(
MATRICE=CO('rigigene'),
MATR_ASSE=Rigi
)),
STOCKAGE='DIAG',
VECT_ASSE_GENE=_F(
TYPE_VECT='FORC',
VECTEUR=CO('vectgene'),
VECT_ASSE=field
)
)
resharm = DYNA_VIBRA(
BASE_CALCUL='GENE',
ETAT_INIT=_F(
DEPL=vectgene,
VITE=vectgene
),
EXCIT=_F(
ACCE=func,
DEPL=func0,
VECT_ASSE_GENE=vectgene,
VITE=func0
),
INCREMENT=_F(
LIST_INST=listr
),
MATR_MASS=massgene,
MATR_RIGI=rigigene,
MODE_CORR=modes0,
SCHEMA_TEMPS=_F(
SCHEMA='NEWMARK'
),
TRAITEMENT_NONL='IMPLICITE',
TYPE_CALCUL='TRAN',
VITESSE_VARIABLE='NON'
)
restran = REST_GENE_PHYS(
CRITERE='RELATIF',
LIST_INST=listr,
NOM_CHAM=('ACCE', 'DEPL', 'VITE'),
NUME_DDL=nddl,
PRECISION=1e-06,
RESU_GENE=resharm
)
○結果
シェルの厚さ6mm~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/akiyama/6mmnozishin.png~
シェルの厚さ20mm~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/akiyama/20mmnozishin.png~
地震のような挙動は示している。~
しかしシェルの厚さが変化しているにもかかわらず変位に変化がない。~
これはおかしい。~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/akiyama/cad6mmtottori.png~
青木先生がキャドで解析してくれたもの~
挙動は比較的同じだが値が違う~
シェル要素の厚さ(断面特性)を反映できていない~
今使用しているDYNA_VIBRA(geneのtran)には~
elemprop = AFFE_CARA_ELEM(シェル要素の断面の特性)を適用するコマンドがない~
※GENE=モードベース(おそらくモード重ね合わせ法)~
PHYS=物理計算ベース?~
harm=定常波~
tran=非定常波~
code_aster(DYNA_VIBRA等)にも特にシェル要素に関する記述はない~
シェルではできない?→→そんなことあります?~
○インターンシップから帰ってきたらやること~
・シェル要素の断面特性の反映~
・できなかったらDYNA_VIBRAのPHYSで解析~
・それでもだめだったらDYNA_LINE(時刻歴応答解析)の使用~
○今後の課題~
①手計算できる1自由度系のモデルで解析、結果の比較~
②斜張橋モデルの作成、線形解析~
③斜張橋モデルの地震解析~
**ぼやき [#x0ac91c1]
・salomeが解析中に3回に1回は固まる。パソコン重いのか?~
・中間発表、特に発表するような結果が得られていない。まずい~
・1自由度まではでかして中間発表にしたい。厳しい戦いだ。~
**秋山の簡易モデルに地震をかけてみる [#je02b8a9]
計算回った。~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/akiyama/kanizisin.gif~
しかし、スケールファクター大きくしないと変形しない。(これは60万倍)~
根元が揺れておらず、先端のみが揺れている。~
○課題~
・加速度がどこにかかってしまっているのか。~
・どこかが間違っている。~
○以下秋山の計算が回ったデータ~
・mesh、model definitionは今までどおり~
・material~
→質量マトリクス(密度)、剛性マトリクス(ヤング率)、減衰マトリクス(今回は剛性減衰)の数値を入れる。~
→鋼材は密度が7.8e-9、MOR_ALPHA(剛性減衰)が2%、ヤング率は200MPa。~
materK = DEFI_MATERIAU(
ELAS=_F(
AMOR_ALPHA=0.02,
E=200000.0,
NU=0.3,
RHO=7.8e-09
)
)
materH = DEFI_MATERIAU(
ELAS=_F(
ALPHA=1.2e-05,
AMOR_ALPHA=0.02,
E=195000.0,
NU=0.3,
RHO=7.8e-09
)
)
・functions and lists~
→加速度のデータを入力する。今回は兵庫県南部地震の加速度を用いている。~
→define function でcsvファイルから加速度を入力。~
→DEFI_LIST_REEL で初期値などを設定。~
→CALC_FONCTION は加速度を変動荷重にしている。値は例題ファイルから。ないと計算が回らない。~
→ここにはどこにかけるかみたいなのはなさそう。~
○DEFI_FONCTION code_aster~
→https://www.code-aster.org/doc/default/en/man_u/u4/u4.31.02.pdf~
○DEFI_LIST_REEL~
→https://www.code-aster.org/doc/default/en/man_u/u4/u4.34.01.pdf~
○CALC_FONCTION code_aster~
→https://www.code-aster.org/doc/default/en/man_u/u4/u4.32.04.pdf~
func = DEFI_FONCTION(
NOM_PARA='INST',
NOM_RESU='acce',
VALE=(0.0, 0.0, 0.02, 0.3, 0.04, 0.3, 0.06, 0.0, 0.08, 0.0, 0.1, -0.3, 0.12, -0.6, 0.14, -2.09, 0.16, -3.29, 0.18, -3.59, ...)
) # sequences have been limited to the first 20 occurrences.
listr = DEFI_LIST_REEL(
DEBUT=0.0,
INTERVALLE=_F(
JUSQU_A=30.0,
PAS=0.02
)
)
func2 = CALC_FONCTION(
COMB=_F(
COEF=1.47,
FONCTION=func
)
)
・BC and Loads~
load = AFFE_CHAR_MECA(
DDL_IMPO=(_F(
DX=0.0,
DY=0.0,
DZ=0.0,
GROUP_MA=('kotei1', )
), _F(
DRX=0.0,
DRY=0.0,
DRZ=0.0,
DX=0.0,
DY=0.0,
DZ=0.0,
GROUP_NO=('kotei2', )
)),
LIAISON_UNIF=_F(
DDL=('DRY', 'DRX', 'DRZ'),
GROUP_NO=('ketsuten', 'kotei2')
),
MODELE=model,
PESANTEUR=_F(
DIRECTION=(0.0, 1.0, 0.0),
GRAVITE=9800.0
)
)
・Post Processing
ondo = CREA_CHAMP(
AFFE=_F(
GROUP_NO=('ketsuten', 'kotei2'),
NOM_CMP=('TEMP', ),
VALE=(-292.0, )
),
MODELE=model,
OPERATION='AFFE',
TYPE_CHAM='NOEU_TEMP_R'
)
・material
fieldmat = AFFE_MATERIAU(
AFFE=(_F(
GROUP_MA=('kebu', ),
MATER=(materH, )
), _F(
GROUP_MA=('ketsusen', 'keta'),
MATER=(materK, )
)),
AFFE_VARC=_F(
CHAM_GD=ondo,
GROUP_MA=('kebu', ),
NOM_VARC='TEMP',
VALE_REF=0.0
),
MAILLAGE=mesh,
MODELE=model
)
・pre alalysis~
○ASSEMBLAGE code_aster~
→https://www.code-aster.org/V2/doc/v13/en/man_u/u4/u4.61.21.pdf~
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'
), _F(
MATRICE=CO('amo'),
OPTION='AMOR_MECA'
)),
MODELE=model,
NUME_DDL=CO('nddl')
)
・analysis~
reslin = MECA_STATIQUE(
CARA_ELEM=elemprop,
CHAM_MATER=fieldmat,
EXCIT=_F(
CHARGE=load,
FONC_MULT=func
),
MODELE=model
)
・post procesing
unnamed0 = CALC_CHAMP(
CARA_ELEM=elemprop,
CHAM_MATER=fieldmat,
FORCE=('FORC_NODA', 'REAC_NODA'),
MODELE=model,
RESULTAT=reslin
)
・BC and Loads→地震荷重の設定~
○CALC_CHAR_SEISME code_aster~
→https://www.code-aster.org/doc/default/en/man_u/u4/u4.63.01.pdf~
field = CALC_CHAR_SEISME(
DIRECTION=(0.0, 1.0, 0.0),
MATR_MASS=mass,
MONO_APPUI='OUI'
)
・analysis~
○CALC_MODES code_aster~
→https://www.code-aster.org/V2/doc/v13/en/man_u/u4/u4.52.02.pdf~
modes = CALC_MODES(
CALC_FREQ=_F(
NMAX_FREQ=10
),
MATR_MASS=mass,
MATR_RIGI=stifness,
OPTION='PLUS_PETITE',
VERI_MODE=_F(
STOP_ERREUR='NON'
)
)
・pre analysis~
○PROJ_BASE code_aster~
→https://www.code-aster.org/doc/default/en/man_u/u4/u4.63.11.pdf~
PROJ_BASE(
BASE=modes,
MATR_ASSE_GENE=(_F(
MATRICE=CO('massgene'),
MATR_ASSE=mass
), _F(
MATRICE=CO('rigi'),
MATR_ASSE=stifness
), _F(
MATRICE=CO('amo1'),
MATR_ASSE=amo
)),
NB_VECT=10,
VECT_ASSE_GENE=_F(
TYPE_VECT='ACCE',
VECTEUR=CO('vectgene'),
VECT_ASSE=field
)
)
・analysis~
○DYNA_VIBRA~
→https://www.code-aster.org/doc/default/en/man_u/u4/u4.53.03.pdf~
resharm = DYNA_VIBRA(
BASE_CALCUL='GENE',
EXCIT=_F(
ACCE=func,
DEPL=func,
VECT_ASSE_GENE=vectgene,
VITE=func
),
INCREMENT=_F(
LIST_INST=listr
),
MATR_AMOR=amo1,
MATR_MASS=massgene,
MATR_RIGI=rigi,
SCHEMA_TEMPS=_F(
BETA=0.25,
GAMMA=0.5,
SCHEMA='NEWMARK'
),
TYPE_CALCUL='TRAN',
VITESSE_VARIABLE='NON'
)
・post processig~
○REST_GENE_PHYS code_aster~
→https://www.code-aster.org/doc/default/en/man_u/u4/u4.63.31.pdf~
restran = REST_GENE_PHYS(
CRITERE='RELATIF',
LIST_INST=listr,
PRECISION=1e-06,
RESU_GENE=resharm,
TOUT_CHAM='OUI'
)
・out put~
IMPR_RESU(
FORMAT='MED',
RESU=(_F(
CARA_ELEM=elemprop,
NOM_CHAM=('DEPL', 'VITE'),
RESULTAT=restran
), _F(
CARA_ELEM=elemprop,
NOM_CHAM=('ACCE', ),
RESULTAT=modes
), _F(
CARA_ELEM=elemprop,
NOM_CHAM=('FORC_NODA', 'REAC_NODA'),
RESULTAT=unnamed0
)),
UNITE=80
)
・https://www.youtube.com/watch?v=Me-5_MFKGyg~
このyoutubeのやり方でやってみた。(7月5日)~
***及川さんからいただいたデータ[#bae3e4bf]
・mesh、model definitionは今までどおり~
・material
mater = DEFI_MATERIAU(
ELAS=_F(
AMOR_ALPHA=0.03,
E=7500.0,
NU=0.4,
RHO=3.8e-10
)
)
mater0 = DEFI_MATERIAU(
ELAS=_F(
E=10000.0,
NU=0.3,
RHO=1e-09
)
)
fieldmat = AFFE_MATERIAU(
AFFE=(_F(
GROUP_MA=('Box_2', 'translation_2', 'translation_3', 'Box_3'),
MATER=(mater, )
), _F(
GROUP_MA=('translation_1', 'Box_1'),
MATER=(mater0, )
)),
MAILLAGE=mesh,
MODELE=model
)
・functions and lists~
func4 = DEFI_FONCTION(
NOM_PARA='INST',
NOM_RESU='acce',
VALE=(0.0, 0.0, 0.02, 0.03, 0.04, 0.03, 0.06, 0.03, 0.08, 0.03, 0.1, 0.03, 0.12, 0.03, 0.14, 0.03, 0.16, 0.03, 0.18, 0.03, ...)
) # sequences have been limited to the first 20 occurrences.
func5 = CALC_FONCTION(
COMB=_F(
COEF=1.47,
FONCTION=func4
)
)
listr = DEFI_LIST_REEL(
DEBUT=0.0,
INTERVALLE=_F(
JUSQU_A=30.0,
PAS=0.02
)
)
・BC and Loads~
load = AFFE_CHAR_MECA(
DDL_IMPO=_F(
DX=0.0,
DY=0.0,
DZ=0.0,
GROUP_MA=('kotei', )
),
MODELE=model
)
・pre alalysis
ASSEMBLAGE(
CHAM_MATER=fieldmat,
CHARGE=(load, ),
MATR_ASSE=(_F(
MATRICE=CO('mass'),
OPTION='MASS_MECA'
), _F(
MATRICE=CO('stifness'),
OPTION='RIGI_MECA'
), _F(
MATRICE=CO('amo'),
OPTION='AMOR_MECA'
)),
MODELE=model,
NUME_DDL=CO('nddl')
)
・BC and Loads
field = CALC_CHAR_SEISME(
DIRECTION=(-1.0, 0.0, 0.0),
MATR_MASS=mass,
MONO_APPUI='OUI'
)
・analysis
modes = CALC_MODES(
CALC_FREQ=_F(
NMAX_FREQ=10
),
MATR_MASS=mass,
MATR_RIGI=stifness,
OPTION='PLUS_PETITE',
VERI_MODE=_F(
STOP_ERREUR='NON'
)
)
・pre analysis
PROJ_BASE(
BASE=modes,
MATR_ASSE_GENE=(_F(
MATRICE=CO('MASSGENE'),
MATR_ASSE=mass
), _F(
MATRICE=CO('rigi'),
MATR_ASSE=stifness
), _F(
MATRICE=CO('amo1'),
MATR_ASSE=amo
)),
NB_VECT=10,
VECT_ASSE_GENE=_F(
TYPE_VECT='ACCE',
VECTEUR=CO('vectgene'),
VECT_ASSE=field
)
)
・ analysis
resharm = DYNA_VIBRA(
BASE_CALCUL='GENE',
EXCIT=_F(
ACCE=func4,
DEPL=func4,
VECT_ASSE_GENE=vectgene,
VITE=func4
),
INCREMENT=_F(
LIST_INST=listr
),
MATR_AMOR=amo1,
MATR_MASS=MASSGENE,
MATR_RIGI=rigi,
SCHEMA_TEMPS=_F(
BETA=0.25,
GAMMA=0.5,
SCHEMA='NEWMARK'
),
TYPE_CALCUL='TRAN'
)
・post processig
restran = REST_GENE_PHYS(
CRITERE='RELATIF',
LIST_INST=listr,
PRECISION=1e-06,
RESU_GENE=resharm,
TOUT_CHAM='OUI'
)
・out put
IMPR_RESU(
FORMAT='MED',
RESU=(_F(
NOM_CHAM=('DEPL', 'VITE'),
RESULTAT=restran
), _F(
NOM_CHAM=('ACCE', ),
RESULTAT=modes
)),
UNITE=80
)
***メモ [#f3341d0b]
・モーダル解析というものがヒットする。~
モーダル解析とは、~
「物体は、なんらかの外力や加振力を加えると、振動現象が発生する。~
いったんこのような振動が起きると、~
力を加えなくても物体は振動し続ける自由振動という現象が起きる。~
この現象には、その物体の動特性に関する情報、つまり動的性質が含まれている。~
物体はそれぞれ固有の動的性質を持ち、モーダル解析とはこれらを求めるのに有効な手法である。」らしい。~
http://opencae.gifu-nct.ac.jp/pukiwiki/index.php?plugin=attach&refer=DEXCS-Salome-ReleaseNote&openfile=Step2-%A5%E2%A1%BC%A5%C0%A5%EB%B2%F2%C0%CF%A4%CE%B2%F2%C0%CF%BC%EA%BD%E7.pdf~
こちらより。
***参考文献 [#v3c7ccb0]
・salome フランス語対策~
http://opencae.gifu-nct.ac.jp/pukiwiki/index.php?Salome-%A5%D5%A5%E9%A5%F3%A5%B9%B8%EC%C2%D0%BA%F6~
・「code_aster earth quake で検索 地震信号の生成について」~
https://www.code-aster.org/doc/default/en/man_r/r4/r4.05.05.pdf~
・「salome-meca esrth quake」で検索~
http://www.opencae.or.jp/wp-content/uploads/2019/12/03-stepbystepcalculation.pdf~
・「code_aster seismic で検索」~
・「code_aster 過渡応答解析~
https://www.code-aster.org/doc/default/en/man_r/r4/r4.05.01.pdf~
・Operator LIRE_FONCTION~
https://www.code-aster.org/doc/default/en/man_u/u4/u4.32.02.pdf~
・動的の計測結果の出力~
https://www.code-aster.org/V2/doc/default/en/man_r/r7/r7.20.02.pdf~
・なんかいいこと書いてる~
http://www.opencae.or.jp/wp-content/uploads/2019/12/05-dynamics.pdf~
・地震解析のメモ~
https://www.jaee.gr.jp/jp/wp-content/uploads/2012/02/kaishi03.pdf~
・DYNA_LINE 和訳わかりやすい~
https://biba1632-gitlab-io.translate.goog/code-aster-manuals/docs/user/u4.53.05.html?_x_tr_sl=en&_x_tr_tl=ja&_x_tr_hl=ja&_x_tr_pto=sc~
・なんかよくわからんけどいい感じのやつ
・code_aster 地震の研究について~
https://www.code-aster.org/doc/default/en/man_u/u2/u2.06.10.pdf~
・code_aster 例題一覧~
https://www.code-aster.org/V2/doc/default/en/index.php?man=R~
・線形の動的における減衰のモデル by code_aster~
https://www.code-aster.org/V2/doc/default/en/man_r/r5/r5.05.04.pdf~
**簡易モデルの作成 [#i1129131]
***課題 [#f310942c]
1 どうやってモーメントを見るのか。何がモーメントなのか。~
2 モーメントを見れないためヒンジになっているのか。~
***曲げモーメントの出力 [#ta571fbe]
FORC_NODAのDRX,DRY,DRZが節点のモーメントらしい。~
Exception user raised but not interceptee. !
! The bases are fermees. !
! Type of the exception: error !
! !
! the computation of the option: CARA_SECT_POUT3 !
! is not possible for any the types of elements of the LIGREL. !
このエラーは結合する線を「Assign finite element」「AFFE_CARA_ELEM」に「Assign a material」に含めないと出てくる。~
→載荷時根元に「forc_noda」出てこないから違うっぽい~
・(x,y,z)=(100,100,1000)の片持ちばり(solid)の先端に1kNの荷重をかけたもので試行錯誤~ [#wf5cc98b]
「EFGE」→エラーが発生
! <S> Exception user raised but not interceptee. !
! The bases are fermees. !
! Type of the exception: error !
! !
! Le nom demandé existe déjà dans la base &&CALCOP.ACALCULER !
! Ce message est un message d'erreur développeur. !
! Contactez le support technique. !
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
「SIEF」→応力と同じ。code_aster曰く変位から応力を計算する。~
code_asterの文章~
Calculation of the state of stress (forced or efforts generalized according to modeling) starting
from displacements (linear elasticity), to see
質問サイト?に載ってあった。~
A curve - as for the moment MFY (bending moment around the Y-axis direction) is plotted as 3d-plot and applied or reaction forces can be
plotted as vector plots.
「SIPM」→pravisにSIPMがない、、。梁要素?の拘束力?らしい~
「SIPO」→エラー、、。これも梁要素?らしい~
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! <S> Exception user raised but not interceptee. !
! The bases are fermees. !
! Type of the exception: error !
! !
! Le nom demandé existe déjà dans la base &&CALCOP.ACALCULER !
! Ce message est un message d'erreur développeur. !
! Contactez le support technique. !
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
「SIRO」→計算は回るけど値が出てこない。~
SIROについて、code_asterの説明~
Calculation of the constraints projected on the skin of a volume (for example on the facings
of a hydraulic work.
和訳
ボリュームの表皮に投影される拘束力の計算(例えば、水力ワークのフェイシングに投影される
を計算する。
「DEGE」→計算は回るけどparavisに出てこない~
「EPME」→計算は回った。コンター図はモーメントっぽいが値が異なる。~
「EPSG」→EPMEと似ている。~
「EPSL」→EPSIと同じ。~
「EPSV」→コンター図が真青になった。~
「DISS」→エラー。
Exception user raised but not interceptee. !
! The bases are fermees. !
! Type of the exception: error !
! !
! Erreur utilisateur : !
! -> Le TYPE_ELEMENT MECA_TETRA4 ne sait pas encore calculer l'option: !
! DISS_ELEM. !
! !
! -> Risques & Conseils : !
! * Si vous utilisez une commande de "calcul" (THER_LINEAIRE, STAT_NON_LINE, !
! ...), il n'y a pas !
! moyen de contourner ce problème. Il faut changer de modélisation ou émettre !
! une demande d'évolution. !
! !
! * Si c'est un calcul de post-traitement, vous pouvez sans doute "éviter" le !
! problème !
! en ne faisant le post-traitement que sur les mailles qui savent le faire. !
! Pour cela, il faut sans doute utiliser un mot clé de type "GROUP_MA". !
! S'il n'y en a pas, il faut faire une demande d'évolution. !
「ENEL」→違う。弾性エネルギーの計算らしい。~
「ENTR」→計算は回るけどpravisに出てこない。~
「ETOT」→違う。~
「EPEQ」→違う。~
「EPGQ」→違う。~
「EPMQ」→違う。EPGQ と同じ感じ。~
「SIEQ」→違う。~
「VARC」→違う。~
「MATE」→違う。~
https://www.code-aster.org/doc/default/en/man_u/u2/u2.01.05.pdf~
「effort code_aster」で検索した。effortの説明とその他にも「SIEF」の説明などもある。~
***結合方法の模索 [#gb0b2b04]
LIAISON_ELEでケーブルと桁を結合しているが、剛結になっていると思われる。~
「LIAISO_UNIF」というものがノードの自由度を同じにできるらしい~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/akiyama/2022.6.10.2.png~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/akiyama/2022.6.9.DRX.png~
→やってみたところ、変形の形もよく、結合箇所のDRXのREAC_NODA(節点反力)は発生しなかった。~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/akiyama/2022.6.10.1.png~
→しかし載荷しているにもかかわらず、根元部分にはDRXのREAC_NODAは発生していない。~
→REAC_NODAのDRXがモーメントと思っていたが違う可能性がある。~
・現時点の課題~
1 モーメントはどうやって見るのか?~
2 モーメントがわからないためヒンジになっているかわからない~
→2.3週間手こずっているため、とりあえず同時進行で地震をかけていく~
***結合の問題 [#r0e6e666]
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/akiyama/2022.5.27.pan.png~
プレストレスをかけることに成功したが、変形が食パンのような形になっていておかしい。~
→結合する場所を端から少しずらした。~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/akiyama/2022.5.27.sen.png~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/akiyama/2022.5.27.pan2.png~
しかし、食パン型は直らず、結合してる線に沿ってリアクション(節点反力)が発生してる。~
→違う結合方法の模索が必要~
・メッシュを切ってからノードを作り、エッジを作成した。~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/akiyama/2022.5.27.1.png~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/akiyama/2022.5.27.2.png~
→結合してないっぽい。~
・メッシュを切ってからエッジを作り、線と点で結合~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/akiyama/2022.5.27.3.png~
→線作るために作ったノードが自由度を持つらしい。いじってないしいじれない。~
Exception user raised but not interceptee. !
! The bases are fermees. !
! Type of the exception: error !
! !
! impossibility, the node N318 does not carry the degree of freedom of rotation !
! DRX
***ピン接合 [#bbe37cf7]
・5月18日~
梁要素にPOU_D_E(節点に回転を許さない)ではなく、BARREを適用(節点に回転を許す)~
→境界条件で固定端をDX-=0,DY=0,DZ=0とした。(DRX=0とかとするとエラーが発生する。)~
→変位がより大きくなった。~
→これはなぜか?ヒンジになっているのか?~
→DIST_TRとK_TR_D_Lでやればうまく行くらしい...~
→この方法は複雑故断念。~
→今日は心が折れたので明日この方法について調べる。~
・5月19日~
このサイトに書いてること試す~
https://www-code--aster-org.translate.goog/forum2/viewtopic.php?id=13266&_x_tr_sl=en&_x_tr_tl=ja&_x_tr_hl=ja&_x_tr_pto=sc~
→「LIALISON_DDL」で節点をヒンジにできるらしい~
→「DDL」のところにDXとかいれるらしい。「COEF_MULT」は1と-1だが、これが何かは不明~
以下、参考にしているaster_studyの設定~
...._ F(GROUP_NO =('NODOF001'、'NODOF002'、)、
DDL =('DY'、'DY'、)、
COEF_MULT =(1、-1、)、
COEF_IMPO = 0、)、)、);
→「COEF_MULT」は下のサイトの3.9.1の数式の係数の実数らしい。~
[https://www.code-aster.org/V2/doc/v14/en/man_u/u4/u4.44.01.pdf]~
3.9.1の数式を参考にする。「LIAISON_DDL」は5.7に記載してある。~
・5月20日
BARRE要素(トラス要素)での解析が成功~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/akiyama/2022.20.1.png~
→プレストレスがかかってない~
→次の方法でプレストレスをかける~
tensyon = CREA_CHAMP(
AFFE=_F(
NOM_CMP=('TEMP', ),
TOUT='OUI',
VALE=(-200.0, )
),
MODELE=model,
OPERATION='AFFE',
TYPE_CHAM='NOEU_TEMP_R'
)
AFFE=(_F(
GROUP_MA=('ketsusen', 'keta'),
MATER=(materK, )
), _F(
GROUP_MA=('hari', ),
MATER=(materH, )
)),
AFFE_VARC=_F(
CHAM_GD=tensyon,
GROUP_MA=('hari', ),
NOM_VARC='TEMP',
VALE_REF=0.0
),
MAILLAGE=mesh,
MODELE=model
)
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/akiyama/2022.5.20.png~
プレストレスがかかり、上方向にたわんだ。~
今後の課題:簡易モデルで地震荷重をかける、曲げモーメントの出力~
***簡易モデルの作成 [#a2316691]
(x,y,z)=(100,100,1000)の簡易モデルをケーブルを貼り要素で、桁をシェル要素で作成。~
→梁要素、シェル要素ともに応力、ひずみが計算できないため、試行錯誤中~
→シェル要素の応力、ひずみの計算、出力に成功~
・梁要素は絶賛試行錯誤中~
→梁は軸力(FORCE)を見ることができたため、OK~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/akiyama/2022.17.1.png~
・シェル部分の応力とひずみを見るため試行錯誤中~
→シェル部分の応力とひずみを見ることに成功~
→しかし、変位が角田さん&梅宮くんのと異なる(符号が違う)~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/akiyama/2022.17.3.png~
→剛結になっている?~
今後の課題:ピン接合~
**参考文献 [#wee9319a]
「ケーブル腐食を考慮した斜張橋の終局強度および疲労寿命」~
構造工学論文集 Vol. 67A (2021 年 3 月)~
中村俊一、青木由香利~
https://www.jstage.jst.go.jp/article/structcivil/67A/0/67A_1/_pdf/-char/ja~
・斜張橋においてケーブルは重要な役割を果たす。~
・ケーブルが腐食すると、じん性、耐力、疲労強度が低下し、断面積が減少する。~
・異なるスパン長のモデルで、腐食(断面積の減少)の影響、ケーブルの疲労寿命を求める。~
・断面積の減少に比例して終局強度は低下する。~
・ケーブルの伸びが1%(健全時の伸びは4%)まで低下すると終局強度は小さくなる。~
・ケーブルの腐食が進行するとともに疲労寿命は小さくなり、断面積の減少を考慮すると疲労寿命は著しく低下する。~
「吊橋ケーブルの腐食機構に関する研究」~
土木学会論文集 No.637~
古家和彦、北川信、中村俊一、鈴村恵太、聖生守雄~
https://www.jstage.jst.go.jp/article/jscej1984/1999/637/1999_637_103/_pdf/-char/ja
*2021年度 創造工房実習の内容 [#wc524cff]
**単純梁(11月12日の課題) [#jdb7873f]
・箱の図心位置に穴を開け、拘束、載荷を行う方法(より梁モデルに近い)
,メッシュの長さ,要素数,変位,相対誤差(%)
,10,4830,0.24068,-42.235
,8,3635,0.29441,-29.341
,5,5909,0.33253,-20.191
,3,7721,0.37432,-10.161
,2,11760,0.38564,-7.4462
,1.5,38147,0.41103,-1.3520
,1,49446,0.41548,-0.2834
,0.5,273903,0.42774,2.6583
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/akiyama/g%20hako.png~
・線を作成し拘束、載荷を行う方法~
,メッシュの長さ,要素数,変位,相対誤差(%)
,10,499,0.199991,-52
,8,1999,0.286,-31.36
,3,7061,0.368392,-12.678
,2,8859,0.374296,-10.168
,1.5,32341,0.404734,-3.026
,1,47869,0.4093403333,-1.756
,0.5,339770,0.4236773333,1.682
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/akiyama/graph1.png~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/akiyama/g%20hako%20sen.png~
梁モデルに近づけた場合のほうが、より理論値に近くなった。(メッシュ長0.5の場合を除く)~
**片持ち梁(11月5日の課題) [#v6cb4f0e]
,メッシュの長さ,要素数,変位,相対誤差(%)
,0.5,340621,6.57939,-1.495
,1.0,588506,6.44299,-3.355
,1.5,17753,6.121,-8.707
,2,8361,5.7382525,-13.926
,3,4863,5.7033575,-14.449
,5,1715,4.8237,-27.645
,8,565,4.0411725,-39.382
,10,376,3.4624,-48.064
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/akiyama/graph2.png
**3次元プロット(10月29日の課題) [#y5466bb4]
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/akiyama/L.png
・h型のグラフを作成
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/akiyama/h.png
・I型グラフの作成
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/akiyama/%20I.png
*メモ [#ub88132f]
**シェル要素メモ [#y5e6ad76]
・以下の方法でシェル要素の応力を計算、出力可能~
unnamed1 = CALC_CHAMP(
CARA_ELEM=elemprop,
CHAM_MATER=fieldmat,
CONTRAINTE=('SIGM_ELNO', ),
DEFORMATION=('EPSI_ELNO', ),
GROUP_MA=('keta', ),
MODELE=model,
RESULTAT=reslin)
unnamed2 = POST_CHAMP(
EXTR_COQUE=_F(
NIVE_COUCHE='SUP',
NOM_CHAM=('EPSI_ELNO', 'SIGM_ELNO'),
NUME_COUCHE=1
),
GROUP_MA=('keta', ),
RESULTAT=unnamed1)
unnamed3 = CALC_CHAMP(
CARA_ELEM=elemprop,
CHAM_MATER=fieldmat,
CONTRAINTE=('SIGM_NOEU', ),
DEFORMATION=('EPSI_NOEU', ),
GROUP_MA=('keta', ),
MODELE=model,
RESULTAT=unnamed2)
output
IMPR_RESU(
FORMAT='MED',
RESU=(_F(
CARA_ELEM=elemprop,
NOM_CHAM=('DEPL', ),
RESULTAT=reslin,
TOUT='OUI'
), _F(
CARA_ELEM=elemprop,
NOM_CHAM=('FORC_NODA', 'REAC_NODA'),
RESULTAT=unnamed4
), _F(
CARA_ELEM=elemprop,
GROUP_MA=('keta', ),
NOM_CHAM=('EPSI_NOEU', 'SIGM_NOEU'),
RESULTAT=unnamed3
)),
UNITE=2
)
**個人的メモ [#m2deea7a]
参考資料~
・梁要素解説~
https://sites.google.com/site/codeastersalomemeca/home/code_aster-1/bimu-youso/pou_d_e-youso~
・梁要素例題~
https://aytechlab.com/homefem/salomemeca2020st009/~
・シェル要素例題~
https://www-argosrl-eu.translate.goog/sito/blog-innovazione/tutorial-elementi-shell-in-code-aster-coque-3d-o-dkt?_x_tr_sl=it&_x_tr_tl=ja&_x_tr_hl=ja&_x_tr_pto=sc~
・post_champ code aster~
https://www.code-aster.org/V2/doc/v13/en/man_u/u4/u4.81.05.pdf~
・要素の解説、例題~
http://www.opencae.or.jp/wp-content/uploads/2019/12/07a-structural-elements.pdf~
・CALC_CHAMP code aster~
https://www.code-aster.org/V2/doc/v13/en/man_u/u4/u4.81.04.pdf~
・DIS_T code_aster~
https://www.code-aster.org/doc/default/en/man_u/u3/u3.11.02.pdf~
・LIAISON_DDL code_aster~
https://www.code-aster.org/V2/doc/v14/en/man_u/u4/u4.44.01.pdf~
・LIAISON_DLL ドイツ語故変な日本語になっております~
https://www-caelinux-org.translate.goog/wiki/index.php/Contrib:KeesWouters/bc/pythonlist?_x_tr_sch=http&_x_tr_sl=de&_x_tr_tl=ja&_x_tr_hl=ja&_x_tr_pto=sc~
・code_aster コマンドリスト~
http://opencae.gifu-nct.ac.jp/pukiwiki/index.php?plugin=attach&refer=SALOME-Meca%A4%CE%BB%C8%CD%D1%CB%A1%B2%F2%C0%E2&openfile=00-00.pdf~
・LIAISON_DDL~
https://docs.google.com/viewer?a=v&pid=sites&srcid=ZGVmYXVsdGRvbWFpbnxjb2RlYXN0ZXJzYWxvbWVtZWNhfGd4OmYwMjAyNDVlMmQxYTE1Yw~
・「salome beam pin」で検索~
http://opencae.gifu-nct.ac.jp/pukiwiki/index.php?plugin=attach&refer=AboutEasyISTR&openfile=EasyISTR5-manual-3.28-210608.pdf~
**英文 [#re901669]
・LIAISON_UNIF~
The keyword LIAISON_UNIF allows to impose with all the degrees of freedom of all the provided
nodes the same value (unknown). One will have thus:U 1( N 1)=U i ( N k ) for all them degrees of freedom Ui indicated on the nodes N
**salome_meca メモ [#ccc1c478]