#author("2022-05-26T16:41:03+09:00","default:kouzouken","kouzouken")
#author("2023-02-16T10:40:35+09:00","default:kouzouken","kouzouken")
#contents

$v=\frac{P\ell^{3}}{48EI}$

*卒業研究「」 [#r9969fbe]
*卒業研究「木材を用いたハイブリッド梁の数値的基礎解析」 [#r9969fbe]

**2✕6ハイブリッド材の弾塑性解析(elastoplastic analysis)[後に節を考慮] [#o75e2db1]
-解析手法
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/magemodel.png
-解析モデル
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/danmen.png
 スパン長2000mm、2✕6材を想定して断面140mm✕114mm、鋼材を2本入れるモデルと4本入れるモデルでそれぞれ鋼材直径10mmと20mmの条件で解析を行う
-設定した応力ひずみ曲線(stress-strain curve)
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/SIEPline_moku1.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/SIEPline_kou1.png
-木材のみの解析結果
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/kekka_mokunomi1.png
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tamura/kekka_mokunomi1.png
-鋼材入りの解析結果
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/kekka_2-10.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/kekka_2-20.png
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/kekka_2-10_1.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/kekka_2-20_1.png
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/kekka_4-10_1.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/kekka_4-20_1.png

**載荷試験 [#xa893462]
-今まで使っていた載荷試験機の限界荷重50kNでも試験体を破壊することができなかったためより荷重をかけれる試験機を整備
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/sikenki.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/sikenki1.png
-まだうまく動かせるかわからないので試験機の情報収集と予備試験を行った後中断している実験の再開
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/sikenki2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/sikenki3.png
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/sikenki4.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/sikenki5.png
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/sikenki6.png
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/sikenki7.png
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/sikenki8.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/sikenki9.png
***試験結果 [#a6008538]
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/sikenkekkaG1.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/sikenkekkaF1.png
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/sikenkekkaD1.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/sikenkekkaC1.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/sikenkekkaB1.png
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/kekkamatome2.png
***破壊後写真 [#b1cb659b]
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/B.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/C.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/D.png
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/F.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/G.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/H.png
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/J.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/K.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/L.png


,試験体名,アラミド枚数,破壊モード,算定ウレタン量[g],高さ[mm],幅[mm]
,A,2,,400,115.28,137.97
,B,2,せん断,800,116.08,137.7
,C,1,せん断+年輪剥離,330,114.23,137.5
,D,0,曲げ+接着面剥離,100,113.12,136.9
,E,2,,450,115.75,138.97
,F,2,曲げ?,470,114.22,138.1
,G,1,曲げ,400,115.47,138.58
,H,0,曲げ,160,114.67,138.65
,I,2,,400,116.17,139.02
,J,2,せん断?,500,116.18,139.17
,K,1,曲げ,250,114.7,138.18
,L,0,せん断?+接着面剥離,60,112.55,138.62





***アラミド引張試験 [#gbf42b76]
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/aramido1.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/aramido2.png
-Ahttp://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/aramidoA.png
・Bhttp://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/aramidoB.png
・Chttp://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/aramidoC.png
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/kekka_aramidoA1.png
・http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/kekka_aramidoB1.png
・http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/kekka_aramidoC1.png
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/kekka_aramidomatome.png



***メモ [#v1d0f06d]
-木材工学発表後想定するツーバイ材の材料定数を指摘されたので、弾性係数6GPa→10GPa、降伏点31MPa→90MPaに直して再計算(引張強度90MPa,圧縮強度40MPa)
***弾塑性エラー [#c000f9c1]
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/D_error.png
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/D_erorr2.png
***salome 弾塑性(引張) [#z41ce150]
-https://www.slideshare.net/stebee19/salomemeca-59878666
-応力についてhttps://d-monoweb.com/blog/evaluate-strength-mises-stress/

**直交異方性弾塑性解析 2022/9/30〜 [#v0d605d0]
-salomeで異方性を考慮した弾塑性解析ができないか
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/hill.comm
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/hill_IS_z.hdf
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/anipara.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/anipara1.png

 2022/11/25
 「F1_A,C_A,F1_M」とかに0を入力して計算は回ったのでF_MXX_XXとかC_MXY_XYの値を変えて解析結果に違いが出るか確認したところほぼ変化なかった
 →F1_A,C_A,F1_Mとかの値を変えて計算してた時は解析結果に違いがあったのでそこに設定するべき値を調査する
  &コンター図に違和感があったので弾塑性解析ができているか確認
 
 *META_LEMA_ANIのLEMAがおそらくLemaitreモデルのことではないかと思われる→そこらへんを調べてF1_A,C_A,F1_Mとかに入れるべき値を見つけられないか


-[[META_LEMA_ANI例題:https://code-aster.org/V2/doc/v14/en/man_v/v7/v7.20.106.pdf]]

 2022/11/16
 META_LEMA_ANIのパラメータ設定の下部にマトリクスの値を設定するであろう場所を発見し、今まで設定していた「F1_A,C_A,F1_M」とかがまた違う値を設定するところであったと思われる。
 *ここのパラメータを設定する際のFとCの違いはcold-faceとhot-faceの違いであった。「cold=froid[フランス語],hot=chaud[フランス語]」 今のところFとCには同じ値を設定する
 
 とりあえず今まで設定していたパラメータに0を入力してF_MXX_XXとかC_MXY_XYとかに適当に値を入れて計算が回るか確認する

 2022/10/4〜
 簡単な直方体モデルを作成し、salomで異方性パラメータを入れるMETA_LEMA_ANI(Define a materialの中)にcode-asterのマニュアルにのってる等方性にする設定で解析したものと今までの方法で等方性解析したものを比べてみた
 解析がまわったので結果を比べてみたところ全く同じ結果で設定が反映されてるのか分からなかったため異方性パラメータの値をを少し変えて再計算した(*詳しく結果を確認したら微妙な違いが確認できた)
 下のようなエラーが出たのでアドバイスどおりにAnalysis設定上のRESI_GLOB_RELAの値を0.0001だったものを0.001にしたところ計算がまわった
 等方性の時とは解析結果に違いが出たので異方性パラメータの設定は反映されてることがわかった
 異方性パラメータに設定する値を木材のものに近づけるように調査中...
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/IS_error.png
 2022/10/3 
 AnsysのHillの異方性パラメータに関するページとCode-Asterの材料設定に関するガイドページからHillの異方性パラメータを設定すると思われるところを発見
 2つのページを比べるとsalome上で設定するべき値が分かりそうだが、各軸方向の降伏応力と降伏せん断応力が必要になる?
-https://www.cybernet.co.jp/ansys/product/lineup/multiscale/relation/mono009.html
-https://www.code-aster.org/doc/default/en/man_u/u4/u4.43.01.pdf (p.137〜)
 2022/9/30
 Hillの降伏条件という直交異方性を考慮した弾塑性解析に必要なパラメータを設定するところがないor見つけられてない→Hillの降伏条件や直交異方性弾塑性について文献調査↓
-https://www.jstage.jst.go.jp/article/aijs/79/700/79_741/_pdf
-https://archi.hiro.kindai.ac.jp/laboratory/SAL/study/ronbun/2018/nagaoka.pdf
-http://www.opencae.or.jp/wp-content/uploads/2019/12/08-non_linear_capabilities.pdf
-https://www.jstage.jst.go.jp/article/jilm/65/5/65_164/_pdf

-[[後藤メモ221110:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?Salome-Meca%E4%BE%8B%E9%A1%8C%E3%83%95%E3%82%A1%E3%82%A4%E3%83%AB#y0bd5d98]]

**新橋の考案 [#n9e8a5c0]
 ・都会の町並みに合う歩道橋(斜張橋)

-4/18~ 基礎部分を足4本と基礎上部に分けて印刷
-4/11~ 基礎部分のモデリングができたのでCuraでgcodeに変換後印刷
-4/5~ 主塔モデル2の基礎部分の作り込みと3Dプリンターでの印刷
-R4↑
-11/8~ 一旦考えたデザインをsalome上でモデリングする
「今後の目標」
 ・新橋全体のデザインを考えていく
 ・ハイブリッド建材に関して鋼材の本数や太さを変えた時の曲げ、圧縮、座屈の解析を行う
-11/1~ 主塔に注目してデザインを考えていく
-10/25~ 理想とする橋の簡易的なモデルが完成したので、強度や重量を計算しながら断面形状や構造を調整していく~
[[木歩道橋に関する資料:http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tamura/mokuhodo.pdf]]
-10/18~ 主桁と高欄のモデル作成(第1モデル完成)
-10/11~ 見た目でどんな感じか確かめるために何パターンかSalome上で簡単にモデリングしてみる
-R3↑
***3Dプリンター [#u7a31f6a]
,,時間(Cura概算),PLA消費量(Cura概算),
,足1本,8h8min,40g・13.26m,
,足4本,32h32min,160g・53.04m,
,基礎上部,19h25min,96g・32.20m,

***主塔モデル2 [#b907fae9]
 ・トラス構造で高さを出す構造として鉄塔に注目
 ・主塔モデル1から基礎部の安定感を出すために錐形に→エッフェル塔を参考にデザイン
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tamura/toumodel2_1.png
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tamura/toumodel2_2.png
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/toumodel2_2.png
-大体の外形ができた途中経過、この後斜材を付けていく→部材数が多くなり不朽の観点から難あり


***主塔モデル1 [#aa1c84f7]
 モデル1で採用したA型のまま木材もしくはハイブリッド建材でトラス構造のようなデザインに
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tamura/toumodel1_1.png
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tamura/toumodel1_2.png
 最大部材長7m←実際この長さの部材は難しい...
 全体的に不安定な感じがある+そのまま基礎部も木材で作ったらさらに不安定感が出る



***モデル1 [#pc6f6ecc]
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/tamura/model_1_zentai.png
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/model_1_zentai.png
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/tamura/model_1_danmen.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/tamura/model_1_keta.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tamura/model_1_keta.png
 諸元
 [スパン]60m+105m,[幅員]10m
 [材料]木材又は木材と鋼材の複合部材
 [主桁]トラス状(キングポストトラス)
 [塔]A型もしくは逆Y型
 [ケーブル]1面ファン状3段→歩道橋のため景観を重視

 ・橋長や幅員、塔の高さ位置などは奥和歌大橋を参考にモデリング


**ハイブリッド部材の解析[#jc8a8067]
 考案中の新橋モデルに取り入れるハイブリッド部材の強度解析
-https://www.icba.or.jp/kenchikuhorei/2018tsuiho/505_523_9.1.pdf
-https://www.ffpri.affrc.go.jp/labs/etj/karube/Mokkyo/pdf/20000927JTE42.KaketudukeruJidai.pdf

-1/17~ 2×8材での解析も行い比較対象を増やして、強度,重量,値段について照査する
-応力のまとめ方を縁応力をみてまとめる→2×4と2×6での挙動の違いがなくなった
-2✕材を用いた鋼-木ハイブリッド部材の曲げ、圧縮、座屈解析を行う
-http://zenmoku-data.jp/pdf/report-03-3.pdf
-取り掛かりとして2✕6材(断面38✕140)の使用を想定して、2✕6を3枚重ねてその間に丸鋼を入れたようなモデルを作成
-12/6~ 2✕4材、2✕6材、2✕8材それぞれに関して丸鋼の直径を10mm~20mmまで2mmずつ変化させた時の解析を行って剛性評価を行う

***2×4と2×6の曲げ解析結果(東北支部2022) [#m7f094b7]
 鋼材直径とたわみ量、鋼材直径と載荷断面上縁と下縁の直応力のグラフ
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tamura/tawami3.png
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tamura/jouen.png
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tamura/kaen.png
 木材のみに対する値をプロットしたグラフ
-たわみhttp://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tamura/tawami_ritu.png
-上縁http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tamura/jouen_ritu.png
-下縁http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tamura/kaen_ritu.png


***2✕6材の解析結果 [#b7152950]
-丸鋼の直径と本数を変えて解析を行った
-モデルは断面114mm✕140mmで長さ2mの部材において、丸鋼の本数を2本と4本でそれぞれ直径12mmと18mmのモデルを作成
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tamura/danmen2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tamura/danmen4.png


,,たわみ[mm],応力(上面)[N/mm^2],応力(下面)[N/mm^2],縮み[mm],座屈荷重(弱軸)[kN],座屈荷重(強軸)[kN]
,12mm-2本,141.957,-82.538,81.188,1.562,73.288,98.511
,12mm-4本,120.608,-76.793,75.709,1.118,86.446,99.813
,18mm-2本,128.670,-83.767,82.284,1.249,80.953,147.965
,18mm-4本,97.472,-73.626,72.177,0.762,107.196,230.765



 
**ねじれ座屈解析 [#ucf4e0dc]
-モデル:10mm(幅)×20mm(厚さ)×1000mm,E=7GPa,ν=0.4
-STAT_NON_LINEで表現できないか試してみる
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/tamura/sidebuckling_4_png/sidebuckling_4_1.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/tamura/sidebuckling_4_png/sidebuckling_4_2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/tamura/sidebuckling_4_png/sidebuckling_4_3.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/tamura/sidebuckling_4_png/sidebuckling_4_4.png
 荷重を上げていくと変形方向が段階的に変わるような挙動が得られた。ねじれ座屈とはとはいえないので別の設定方法を考える。
-CALC_MODEで再挑戦
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/tamura/sidebuckling0.png
 9/1近藤さんにねじれ座屈の設定についてアドバイスをいただき、荷重をかける向きを軸方向から厚さ方向に変えて設定したところねじれ座屈のような挙動が得られた。
 座屈荷重97.849N(1次モード)
 
 ・理論値との比較をしようと「Trahairの座屈公式」から理論値を算出しようとしたが、式に使用されている定数が分からず調べることに...
 ・J=棒部材のねじり定数→GJ=ねじり剛性、Iω=そりねじり定数→EIω=そりねじり剛性であって、ねじりには純ねじりとそりねじりがある。
 ・矩形断面における各定数はどう決定されるのか?→後藤さんからねじり定数とそりねじり定数を計算可能なプログラムを教えてもらったので、それを使用して計算。
-Trahairの座屈公式
$P_{cr}=\frac{\sqrt{EI_yGJ}}{l^2}(3.95+3.52\sqrt{\frac{\pi^2EI_ω}{GJl^2}})$
 9/28~芯材(鋼材)を入れたモデルと穴が空いたモデルの横ねじれ座屈解析をおこなう
 1.断面80mm×50mm、長さ1000mmの木材に断面15mm×15mmの鋼材が入ったモデル
 2.上の木材のみのモデルの側面に、断面60mm×30mmの穴を開けて穴の位置を変えながら解析

+行った解析は比較のために鋼材を入れたモデルの他に木材のみと、鋼材が入る部分が穴あきのモデルの3種類
結果としては、木材のみ21.931kN,穴あき21.931kN,鋼材入り25.261kN
+穴の開け方は線形座屈のときと同様で結果はグラフのようになった。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/tamura/sideanaaki.png



**座屈解析 [#j86f4a98]
-https://hitoricae.com/2021/01/31/%e8%96%84%e6%9d%bf%e7%ad%92%e3%81%ae%e7%b7%9a%e5%bd%a2%e5%ba%a7%e5%b1%88%e8%a7%a3%e6%9e%90/
-https://www.code-aster.org/V2/doc/default/en/index.php?man=commande
 ・上のページを参考にしてAsterStudyの設定をしてみた。マトリクス関連の設定でエラーが出た。「MATR_RIGI」、「MATR_RIGI_GEOM」
 ・上のエラーは線形解析を行うことで解消される
 ・5/31時点でAsterStudyの設定で赤くエラーが表示されることはなくなった。
 (この時のコマンドファイル)http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/tamura/buckling_TR.comm
 ・上のコマンドファイルでは「no eigen mode in the waveband requested.」のようなエラーがでた。固有値が設定されてないっぽい
 ・6/2時点で計算可能になった。
 (この時のコマンドファイル)http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/tamura/buckling_TR1.comm
 改善した理由としてはCALC_MODESのCALC_CHAR_CRITの値を変更したからだと思われる。→CALC_CHAR_CRITに入れる値は固有値とモードの捜索範囲を指定するものらしい。
 ・座屈した時の梁の挙動は解析結果に示したようであるが座屈荷重のアウトプットが目に見えてはっきりしているものでないため検討する必要がある。
✧座屈解析などの固有値解析で求められる固有値$\lambda$(ParaVisの"Time"の値)は、与えた条件に対する倍率なので、$P_{cr}=\lambda \times P$が座屈荷重になります。特に理由がないのであれば、与える条件を単位荷重(P=1kNとか)にしておくと"Time"の値がそのまま座屈荷重になるので楽なんじゃないかと(6/2近藤)。

 ・6/7~は座屈解析ができるようになったので、これを利用した応用的な解析をしていく。まず正方形断面でも解析できることが確認できたので予てから行っているモデルで座屈解析を行っていこうと思う。
 ・木材90mm×90mm×1200mmのモデルで断面の中心に16mm×16mmの鋼材が入ってるモデルの解析を行った結果、2000kN以上の座屈荷重が出てきた。これの示すことは鋼材まで座屈する時の荷重であると考えられる。鋼材が入った際に木材のみが座屈する状態を作る設定を考える必要がある。
 ・確認のため鋼棒のみ、木材のみ、穴の空いた木材のモデルで座屈解析して、解析値と理論値の座屈荷重を比べてみた。
 ・90×90×1200の鋼材を座屈解析したところ2000kN以上の座屈荷重が出たので鋼材が入ったモデルの設定を見直してみたら材料設定のミスが見つかった。設定を直して解析を行ったところ72.1194kN(6.94%)の座屈荷重が得られた。
,,解析値[kN],理論値[kN],相対誤差
,90×90×1200の木材,69.33,65.58,5.72%
,16×16×1200の鋼材,2.744,1.928,42.32%
,断面90×90の中心に16×16の穴,68.99,65.51,5.31%
,90×90×1200の鋼材,2012.4,1929.865,4.28%

 ・6/14~ モデルをより実際に近づけるために、木材に異方性を設定して座屈解析を行う。解析を行ったのは断面90×90mmの木材一材料のモデルと、その中に16×16mmの鋼材が入ったモデルである。
 ・解析を行った結果求められた座屈荷重は、木材のみの方が65.922kN(0.52%)、鋼材の入った方が68.895kN(2.16%)であった。
 ・

 ・6/28~ 要素を変えて座屈解析してみる。
  1.梁要素を使って座屈解析する。←青木さんから他のソフトで解析をした結果理論値に近くなったとのことだったのでsalomeでも確かめてみる
    モデル90×90×1200の木材(ヤング率7GPa,ポアソン比0.4)で行った結果、座屈荷重は65.579kN()と求められた
  2.ソリッド要素を使って、鋼材の位置を中心からずらした場合の座屈解析を行う

 ・7/5~ ねじれ座屈解析
 メモ:Trahairの座屈公式
 STAT_NON_LINEで表現できないか

 ・7/26~ 柱部材に穴が空いている時の座屈解析
 モデル:90×90×1200mmの木材(E=7000,ν=0.4)
 1.穴の形状を変化させたときの座屈荷重の変位を確かめる
 2.穴の位置を変化させたときの座屈荷重の変位を確かめる
+http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/tamura/png/anaaki_W.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/tamura/png/anaaki_S.png
+http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/tamura/png/anaaki2_W.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/tamura/png/anaaki2_S.png





***初めての座屈解析の解析結果 [#z1c2c584]
-モデル断面80mm×60mm、高さ1000mm、E=7GPa、ν=0.4、P=300kN
+http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/tamura/zakutu.png_Files/zakutu0.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/tamura/zakutu.png_Files/zakutu1.png
+http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/tamura/zakutu.png_Files/zakutu2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/tamura/zakutu.png_Files/zakutu3.png
+http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/tamura/zakutu.png_Files/zakutu4.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/tamura/zakutu.png_Files/zakutu5.png
+http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/tamura/zakutu.png_Files/zakutu6.png
-各モードの座屈荷重
,,座屈荷重[kN],
,,弱軸方向,強軸方向,
,1次モード,27.130,46.432,
,2次モード,238.27,400.34,
,3次モード,630.96,1026.0,
,4次モード,1157.0,,

**5/10~ [#qea6f87d]
-もらった実験結果との比較をするためにモデルの作成と解析を行う
-モデルは90mm×90mm×1200mmの柱荷重
+真ん中に断面16mm×16mmの正方形の穴が空いた木材
-200kN荷重時たわみ 解析値:4.350275 理論値:4.370948(-0.473%) 実験値:2.9(50%)
+穴に16mm×16mmの鋼材が入った
-200kN荷重時たわみ 解析値:2.369955 理論値:2.261846(4.78%) 実験値:2.25(5.33%)
-300kN荷重時たわみ 解析値:3.418366 理論値:3.392770(0.754%) 実験値:3.2(6.82%)
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/tamura/KAKUKOU.png
+真ん中に半径8mmの鋼棒が入った
-200kN荷重時たわみ 解析値:2.7934175 理論値:2.5129587(11.16%) 実験値:3.3(-15.35%)
-300kN荷重時たわみ 解析値:4.1901825 理論値:3.7694380(11.16%) 実験値:5.1(-17.84%)
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/tamura/MARUKOU.png

**半径8mmの鋼棒解析 [#te885253]
-たわみ理論値60.228,解析値60.216
-解析値と理論値の差はほぼなく解析できた
-理論値は2材料のときのように木材と鋼材の曲げ剛性を足して計算した
--円の断面二次モーメント→$I=\frac{πr^{4}}{4}$

-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/tamura/tanjun_KOUBOU_model.png
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/tamura/tanjun_KOUBOU.png

-salomeでモデルを変えて解析する場合はコマンドファイルを移す方法が便利
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/tamura/comm_export.png
すでに出来てるコマンドファイルをエクスポート
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/tamura/comm_import.png
新しいモデルでAdd stageする際にインポート→メッシュなどを設定し直す


**4種類解析してみて [#zd31c151]
-どれも理論値と解析値は近い値になったのでうまく解析はできたと思う。
-たわみの値の大きさの順番は木材→サンドイッチ→鋼材だと思っていたが解析結果ではサンドイッチ→木材→鋼材の順だった。この結果に違和感を感じるが理論値がこの結果が正しいことを裏付けている。
-青木さんより、得られたたわみの差は大きな差ではないのでどちらがよりたわんだかはさほど気にすることはないとのこと→次のステップへ

**4/26~ [#aea76f28]
-200×200×3000mmの単純梁に100kN荷重をかけた解析を各材料について行う
-たわみの理論値は以下の式より求める$v=\frac{P\ell^{3}}{48EI}$
-100kN=10.2t→2Lペットボトル5000本分くらい
--①E=200GPa,ν=0.3
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/tamura/tanjun_KOU.png
-理論値:2.109,解析値:2.113
--②E=7GPa,ν=0.4
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/tamura/tanjun_MOKU.png
-理論値:60.268,解析値:60.191
--③木材、鋼材をサンドイッチした際の理論値を計算して仮定のヤング率を設定
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/tamura/tanjun_HB_katei.png
-理論値:62.114,解析値:62.037
--④実際にサンドイッチしたモデルを作成して解析→2材料を参考
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/tamura/tanjun_HB.png
-理論値:62.111,解析値:60.196




**4/19~4/25 [#e25b572e]
-salomeを使った座屈解析の方法を調べる
-https://www.code-aster.org/V2/doc/default/en/man_u/u4/u4.52.02.pdf
-https://www.code-aster.org/V2/doc/default/en/index.php?man=commande
-↑のページを読めば分かるかも?
-座屈は難しいので一旦置いといて簡単なモデルで解析していく

*imagemagick [#g6fe7fc2]
- http://www.protein.osaka-u.ac.jp/rcsfp/supracryst/suzuki/jpxtal/Katsutani/gif1.php
-https://pen.envr.tsukuba.ac.jp/~torarimon/?Imagemagick

*作業日誌 [#s7e0f49e]
,日付,時間,作業時間,内容,立会
,10/16,14:30~16:00,1h,顔合わせ,後藤さん、その他18
,10/23,14:30~16:00,1.5h,UNIXコマンドンの使い方,後藤さん
,10/30,14:30~16:00,1.5h,vi、gnuplotの使い方,後藤さん
,11/4,14:00~16:00,2h,gnuplotの宿題,
,11/6,14:30~16:30,2h,Salome-Meca導入,後藤さん、及川さん
,11/11,14:30~16:00,1.5h,Salomeの宿題,
,11/18,14:30~16:00,1.5h,salomeの宿題,
,11/20,14:30~16:00,1.5h,木材 直交異方性,後藤さん
,11/25,13:00~17:00,4h,salomeの宿題,
,12/4,14:30~17:00,2.5h,2材料,後藤さん
,12/9,14:30~16:00,1.5h,salome宿題,
,12/11,14:30~16:00,1.5h,Tex,
,12/18,13:00~16:30,3.5h,Tex画像の貼り方,
,1/20,13:00~16:00,3h,最終課題,
,1/27,13:00~16:00,3h,最終課題,
,1/29,13:00~19:00,6h,最終課題,
,2/5,14:30~16:00,1.5h,最終課題,後藤さん
,2/10,13:00~16:00,3h,パワポ作成,
,2/12,13:00~16:00,3h,最終課題発表,後藤さん
,2/16,9:00~17:00,8h,卒論発表見学,
,2/24,12:00~17:00,5h,春休み課題,
,3/3,13:00~16:00,3h,春休み課題,
,3/4,13:00~16:00,3h,春休み課題,
,3/5,13:00~16:00,3h,春休み課題,
,3/8,13:00~16:00,3h,春休み課題,
,3/9,13:00~16:00,3h,ゼミ(進捗報告),後藤さん
,3/25,13:00~16:00,3h,春休み課題,
,3/29,13:00~16:00,3h,ゼミ(進捗報告),後藤さん
,4/5,16:00~17:00,1h,顔合わせ(zoom),後藤さん、青木さん、石黒さん
,4/12,14:30~16:30,2h,卒論修論テーマ発表,
,4/14,10:30~16:00,5.5h,英語ゼミ、春休み課題,後藤さん、青木さん




*11/13 Salomeの宿題 [#v53387e8]
**11/13 Salomeの宿題 [#v53387e8]
,メッシュ長さ,要素数,先端変位,相対誤差,計算者
,0.5,394121,6.57938,-1.7%,君島
,0.7,130916,6.4781,-2.9%,君島
,0.8,72101,6.43695,-3.5%,髙橋
,0.9,71718,6.43136,-3.6%,髙橋
,1,72278,6.44302,-3.4,髙橋
,1.2,65575,6.408255,-3.9%,田村
,1.4,41096,6.316155,-5.2%,田村
,1.5,23417,6.120905,-8.2%,根本
,1.8,11758,5.7368975,-13.9%,根本
,2,11817,5.7382525,-13.9%,藤原
,4,2862,4.9428,-25.9%,藤原
,8,897,4.0411725,-39.4%,森島
,10,596,3.4634575,-48.1%,森島


*11/20 Salomeの宿題(単純梁) [#n5deb22c]
**11/20 Salomeの宿題(単純梁) [#n5deb22c]
ヤング率:7500N/mm^2 ポアソン比0.4 断面積 10mm×10mm

,メッシュ長さ,要素数,変位,相対誤差,計算者
,0.5,394121,0.34118,2.35%,君島
,0.7,215780,0.33797,1.39%,君島
,0.8,159468,0.33563,0.69%,髙橋
,0.9,90071,0.33203,-0.39%,髙橋
,1.0,61315,0.32997,-1.2%,田村
,1.2,58111,0.329956,-1.2%,田村
,1.4,47409,0.328156,-1.5%,田村
,1.5,42068,0.325074,-2.4%,根本
,1.8,24627,0.317161,-4.8%,根本
,2,12228,0.3005115,-6.9%,藤原
,4,5077,0.284054,-13.9%,藤原
,8,1795,0.231200,-30.6%,森島
,10,752,0.1612725,-51.6%,森島


http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/tamura/tanjun.png

*11/20課題「異方単純梁」 [#ea56845b]
**11/20課題「異方単純梁」 [#ea56845b]

スパン長100mm、
,メッシュ長さ,要素数,変位,相対誤差,計算者
,0.5,391031,0.44002,11.9%,君島
,0.7,215780,0.435823,10.8%,君島
,0.8,159468,0.43301,10.1%,髙橋
,0.9,90071,0.42766,8.72%,髙橋
,1.0,61315,0.42388,7.77%,田村
,1.2,58111,0.42301,7.54%,田村
,1.4,47409,0.42031,6.86%,田村
,1.5,42068,0.41847,6.39%,根本
,1.8,24627,0.41046,4.36%,根本
,2,12228,0.39631,-0.84%,藤原
,4,5077,0.37870,-3.6%,藤原
,8,1795,0.34230,-12.7%,森島
,10,752,0.29871,-24%,森島

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/tamura/iho-tanjun100.png


スパン長50mm、
,メッシュ長さ,要素数,変位,相対誤差,計算者
,0.5,215781,0.1253048,74.8%,君島
,0.7,109175,0.1172535,63.6%,君島
,0.8,75902,0.11527,60.8%,高橋
,0.9,71911,0.11416,59.3%,高橋
,1,47757,0.113602,58.5%,
,1.2,26945,0.108935,51.99%,田村
,1.4,22998,0.107298,49.71%,田村
,1.5,17689,0.10375025,44.77%,根本
,1.8,14668,0.1021334,42.51%,根本
,2,13986,0.069684,-2.8%,藤原
,4,3009,0.048789575,-32.0%,藤原
,8,967,0.0764429,7.2%,森島
,10,558,0.0768385,6.7%,森島

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/tamura/iho-tanjun50.png

*12月4日課題「2材料」 [#we18e851]
**12月4日課題「2材料」 [#we18e851]

,メッシュ長さ,要素数,変位,相対誤差,計算者
,0.8,226647,0.08053,56.9%,高橋
,0.9,127506,0.07678,49.7%,高橋
,1.0,92447,0.05276,2.79%, 
,1.2,88386,0.05264,2.55%,田村
,1.4,78086,0.05261,2.49%,田村
,1.5,70032,0.072549,41.3%,根本
,1.8,34858,0.068375,33.2%,根本
,2,20313,0.063280,23.3%,藤原
,3,18229,0.0489236,-4.68%,君島
,4,8067,0.050046,-2.51%,藤原
,5,4846,0.036772667,-28.3%,君島
,8,3814,0.0270877,-47.2%,森島
,10,1716,0.0217906,-57.5%,森島

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/tamura/2zairyou1.png

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/tamura/2zairyou2.png

*創造工房課題 [#u85a4816]
**創造工房課題 [#u85a4816]
10mm×10mm木材断面の縦軸中心に1mm×10mmの鋼材が入った、スパン長100mmの単純梁

,メッシュ長さ,要素数,変位,相対誤差
,0.5,405358,0.14421,57.8%
,0.7,303145,0.114400,57.5%
,0.8,206498,0.14331,56.8%
,0.9,120189,0.14176,55.1%
,1.0,82443,0.14114,54.4%
,1.2,80952,0.14110,54.4%
,1.4,71716,0.14109,54.4%
,1.5,68872,0.14092,54.2%
,1.8,31364,0.13522,47.9%
,2,21664,0.13261,45.1%
,4,5821,0.09556,4.6%
,8,3486,0.06404,-30.0%
,10,1842,0.01005,-89.0%


https://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/kadai_TR.png


トップ   新規 一覧 検索 最終更新   ヘルプ   最終更新のRSS