設計図のサイズで同様に解析. 今までの4つのモデルと垂直補剛材+水平補剛材の5つのモデル解析. モデル番号A-E
モデルC(木のみ)に接着剤をつけて解析をしてみる
salomeでの仮想材料での接着方法
・接着,I型鋼などそれぞれでメッシュを切る ➜コンバウントの作成 ➜AsterStudyのBCandloadでLIASON_MAILを選択 ➜GROUP_MA_ESCLは結合したいものの面,GROUP_MA_MAIT一緒に結合するものの立体図形を選択
理論式:Andradeの四季
4種のモデルでの座屈解析の結果,座屈荷重のグラフが次のようになった.
画像はそれぞれのモデルの1次モードでの様子, model0-3の順になってる
model-0-3のそれぞれの弾塑性解析の結果
4種のmodelで弾塑性解析をして応力ーひずみグラフ作成した.
たわみ(mm) | 荷重(N) |
0 | 0 |
2.12661 | 6625 |
4.25322 | 13250 |
6.37983 | 19875 |
8.50644 | 26500 |
10.633 | 33125 |
12.7597 | 39750 |
14.8863 | 46375 |
17.0129 | 53000 |
19.1395 | 59625 |
21.2661 | 66250 |
23.3927 | 72875 |
25.5193 | 79500 |
27.6482 | 86118.2 |
29.7878 | 92759.9 |
31.9419 | 99593.7 |
34.1569 | 106000 |
36.9001 | 112617 |
41.3961 | 119250 |
52.3657 | 125875 |
215.126 | 132500 |
たわみ(mm) | 荷重(N) |
0 | 0 |
2.11578 | 6625 |
4.23155 | 13250 |
6.34733 | 19875 |
8.46311 | 26500 |
10.5789 | 33125 |
12.6947 | 39750 |
14.8104 | 46375 |
16.9262 | 53000 |
19.042 | 59625 |
21.1579 | 66250 |
23.2737 | 72875 |
25.3897 | 79500 |
27.5085 | 86117 |
29.6378 | 92764.6 |
31.7814 | 99779.7 |
33.9866 | 106003 |
36.7164 | 112609 |
41.155 | 119250 |
51.9818 | 125875 |
230.585 | 132500 |
たわみ(mm) | 荷重(N) |
0 | 0 |
2.56608 | 9000 |
5.13216 | 18000 |
7.69824 | 27000 |
10.2643 | 36000 |
12.8304 | 45000 |
15.3965 | 54000 |
17.9626 | 63000 |
20.5286 | 72000 |
23.0947 | 81000 |
25.6608 | 90000 |
28.23 | 99165.7 |
30.8166 | 108003 |
33.4451 | 117001 |
36.4896 | 125995 |
40.7128 | 135000 |
46.8066 | 144000 |
55.144 | 153000 |
66.2457 | 162000 |
85.0673 | 170995 |
169.207 | 179792 |
たわみ(mm) | 荷重(N) |
0 | 0 |
2.53861 | 9000 |
5.07722 | 18000 |
7.61582 | 27000 |
10.1544 | 36000 |
12.693 | 45000 |
15.2316 | 54000 |
17.7703 | 63000 |
20.3089 | 72000 |
22.8475 | 81000 |
25.3862 | 90000 |
27.9256 | 98999.1 |
0.4789 | 108023 |
33.0618 | 117000 |
35.9787 | 125999 |
40.0731 | 135000 |
45.9019 | 144000 |
53.8572 | 153000 |
64.5072 | 162000 |
82.927 | 171000 |
200.608 | 179976 |
リブの接合方法としてフランジとは1.3~1.5m間隔を開けてI桁の片方だけしかつけない.あと曲げモーメントもする.ひずみを見る
道路示方書を読んで適切なようにリブをつけた.
500-8,300-6のリブをつけた.グラフ作成とリブ無いときの差を比較する.
次は500-6,300-6に木材をはめ込む.まずは接着剤考えずにやる.
①全部に木材をはめ込む②1つずつ間隔を開けてはめ込む③両端に多く入れるetc.
スパン10mで横ねじれ座屈解析を進める
①間隔500mm,厚さ6mm②間隔500mm,厚さ8mm③間隔300mm,厚さ6mm
I型鋼を使った過去の卒論での横ねじれ座屈させる
スパン長(mm) | fortran(MN) | salome解析(MN) | 相対誤差 |
5000 | 0.76961 | 0.6736 | 12.5% |
6000 | 0.475371 | 0.443304 | 6.75% |
7000 | 0.318248 | 0.304819 | 4.22% |
8000 | 0.225856 | 0.218265 | 3.36% |
9000 | 0.167514 | 0.16558 | 1.15% |
10000 | 0.12597 | 0.127948 | 1.57% |
0.1✕0.1の小さい穴を開けて片側ピン片側ローラーの単純支持でオイラー座屈をした.
2次要素で0.5~3%ほどの相対誤差になった.
中立軸の高さで固定
1✕1✕1の穴を開けて座屈させた
スパンの余裕をなくし,線で固定した
1✕10✕120の長方形断面でオイラー座屈をした
中村先生のモデルの弾塑性解析,簡単な長方形断面で座屈解析をしている
今後の予定
中村先生のモデルで弾塑性解析をする➜解析結果と中村先生の実験結果と比較して正しいか確認する
自分のモデルに①375mm間隔にリブをつける②木を入れる③リブと木を入れる
中立軸の高さの固定できなかった
面固定の結果
固定のやり方を変えた
・両端の面で固定する
→2次要素にしたら見た目の様子がおかしかった
・中立軸の高さで単純ばりと同じように固定する
→メッシュが上手に切れない
中間発表後の課題
・自分と小池さんの式の理論値の違いを見つける
・要素数が大きすぎる
薄肉プレートの鋼I桁の座屈解析後
・木材をはめ込む→接着剤について考える
・リブの部材に木をはめ込む
1次要素の座屈解析をした。メッシュ3の解析中、メッシュの切り方、要素数の減らし方を検討する。 2次要素は順に解析する。1次、2次それぞれ理論値と合わせてグラフ化する。 理論値の確認する。
I形鋼のモデルの座屈解析が成功した.結果もパラビスで見れた.
メッシュ50で解析した.その時の様子
パラビスでの座屈結果の見方がわかった.
medファイルの2つ目を開く→フィルターからMechanicsのNomal Models animation(complex)にしてこれだけの画面にする→開いたmedファイルをmodeにする→Imaginaryを変えながら見る
座屈解析を10✕1✕120の長方形断面で行った. 解析は成功したが,自分のパソコンだとパラビスで結果が見られない. 千代岡くんのパソコンで見た場合,結果が見え,座屈がしていた.
I型の弾塑性解析を線荷重で行った. 次は線の中心で点荷重する.
英語発表の準備をしていた. 先週で一通りの試験が終わった.卒論頑張る.
座屈モーメント式の計算をした.
次回からはモデルを座屈させる.
国家一般職の1次試験を受けた. あとオープンキャンパスの動画撮影がんばった.
フランジ:12✕200✕6100,ウェブ:900✕6✕6100,のI形鋼のモデルを作成した(スパン長6000).
青木先生から頂いた資料をまとめる. 来週はモデルを作成する
名古屋市役所を受験した. 卒論は青木先生から頂いた資料を読みたいと思う
春休み課題の再発表 モデルの見方を変えた 平均的な見方から最も影響を受けた1点に着目した
14日の春休み課題に向けて研究室に行った。自分は弾塑性班でスライドを作成した
結果のデータはkawai23の中に入ってる.
最初ののモデルの結果は900kekka
丸森の結果と仮想材の結果は2300kekkaにある
引き継ぎの人頑張ってくれーーーーーーーーーー
ぜひヨーロッパ行ってみて
田村さんのやつ:http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/tamura/buckling_TR.comm
wiki: https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?%E5%BA%A7%E5%B1%88%E8%A7%A3%E6%9E%90
メッシュを別々で結合するとき
https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?salome+%E7%B5%90%E5%90%88%E3%83%A1%E3%83%A2
https://www.str.ce.akita-u.ac.jp/~gotou/linux/vine.html#unix
要素数はノード(接点数)を見た
メッシュ長さ | 要素数 | 先端変位 | 相対誤差 | 計算者 |
0.5 | 59504 | 6.56 | -1.5 | 千代岡 |
0.6 | 45512 | 6.48774 | -2.69 | 高井 |
0.7 | 39075 | 6.54133 | -2.0 | 関合 |
0.8 | 13397 | 6.43695 | -3.5 | 岡田 |
0.9 | 9903 | 6.36315 | -4.6 | 松田 |
1.2 | 6256 | 6.3043375 | -5.4 | 青野 |
1.3 | 5767 | 6.29784 | -5.6 | 山口 |
1.4 | 5146 | 6.286015 | -5.76 | 山本 |
1.5 | 3935 | 6.24807 | -6.3 | 進藤 |
1.6 | 3400 | 6.20446 | -6.98 | 河合 |
1.8 | 2952 | 6.17161 | -7.5 | 山口 |
2 | 1632 | 5.64585 | -15.3 | 進藤 |
3 | 667 | 5.405397 | -19.0 | 山本 |
4 | 264 | 3.6161 | -45.8 | 関合 |
5 | 191 | 3.86 | -42 | 千代岡 |
6 | 520 | 6.3660625 | -4.51 | 高井 |
7 | 75 | 1.41225 | -78.8 | 青野 |
8 | 56 | 1.2887175 | -80.7 | 岡田 |
9 | 49 | 1.28799 | -80.9 | 松田 |
10 | 44 | 1.226075 | -81.6 | 河合 |
要素数は四面体の数を見る
メッシュ長さ | 要素数 | 変位 | 相対誤差 | 計算者 |
0.5 | 604167 | 0.4289 | 2.94 | 千代岡 |
0.6 | 361584 | 0.421233 | 1.09 | 高井 |
0.7 | 145234 | 0.4225 | 1.4 | 関合 |
0.8 | 140987 | 0.422627 | 1.4 | 岡田 |
0.9 | 91857 | 0.42035 | 0.88 | 松田 |
1.2 | 24520 | 0.3986 | -2.8 | 青野 |
1.3 | 23132 | 0.40450 | -2.93 | 山口 |
1.4 | 17580 | 0.3986 | -4.34 | 山本 |
1.5 | 15433 | 0.39631 | -4.9 | 進藤 |
1.6 | 15900 | 0.39905 | -4.24 | 河合 |
1.8 | 11677 | 0.404457 | -2.93 | 山口 |
2 | 10406 | 0.39482 | -5.3 | 進藤 |
3 | 2344 | 0.32447 | -22.1 | 山本 |
4 | 1453 | 0.3329 | -20.1 | 関合 |
5 | 431 | 0.13624 | -67.3 | 千代岡 |
6 | 360 | 0.21304 | -48.9 | 高井 |
7 | 196 | 0.101989 | -75.5 | 青野 |
8 | 104 | 0.115862 | -72.2 | 岡田 |
9 | 81 | 0.12470 | -70.1 | 松田 |
10 | 78 | 0.07733 | -81.4 | 河合 |
要素数は四面体の数を見た。ノード数は増えるが四面体は変わらない。
異方性 1次要素 理論値 0.4917(ティモシェンコ梁)
メッシュ | 要素数 | 変位 | 相対誤差 | 計測者 |
0.5 | 604167 | 0.509167 | 3.56 | 千代岡 |
0.6 | 203209 | 0.504716 | 2.6 | 高井 |
0.7 | 145234 | 0.5036 | 2.42 | 関合 |
0.8 | 140987 | 0.5028270 | 2.3 | 岡田 |
0.9 | 91974 | 0.500527 | 1.80 | 松田 |
1.2 | 24800 | 0.4873933 | -0.9 | 青野 |
1.3 | 23132 | 0.488410 | -1.99 | 山口 |
1.4 | 17617 | 0.484033 | -1.56 | 山本 |
1.5 | 15433 | 0.4820229 | -2.0 | 進藤 |
1.6 | 15900 | 0.483285 | -1.7 | 河合 |
1.8 | 11677 | 0.47855 | -1.97 | 山口 |
2 | 10460 | 0.479058 | -2.6 | 進藤 |
3 | 2436 | 0.4278688 | -12.98 | 山本 |
4 | 1453 | 0.42772 | -13.02 | 関合 |
5 | 431 | 0.273640 | -44.3 | 千代岡 |
6 | 360 | 0.3392699 | -31.0 | 高井 |
7 | 196 | 0.213628 | -58.5 | 青野 |
8 | 104 | 0.22574 | -54.1 | 岡田 |
9 | 81 | 0.227502 | -53.7 | 松田 |
10 | 78 | 0.203271 | -58.7 | 河合 |
等方性 2次要素 理論値 0.4167
メッシュ | 要素数 | 変位 | 相対誤差 | 測定者 |
0.6 | 203209 | 0.423827 | 1.72 | 高井 |
0.7 | 145234 | 0.43011 | 3.22 | 関合 |
0.8 | 140987 | 0.430058 | 3.2 | 岡田 |
0.9 | 91974 | 0.42991 | 3.18 | 松田 |
1.2 | 24800 | 0.42978 | 3.14 | 青野 |
1.3 | 23132 | 0.429885 | 2.03 | 山口 |
1.4 | 17617 | 0.42974 | 3.13 | 山本 |
1.5 | 15433 | 0.429844 | 3.2 | 進藤 |
1.6 | 15900 | 0.429754 | 3.13 | 河合 |
1.8 | 11677 | 0.429623 | 2.03 | 山口 |
2 | 10460 | 0.4296050 | 3.1 | 進藤 |
3 | 2486 | 0.4292165 | 3.00 | 山本 |
4 | 1453 | 0.4293 | 3.02 | 関合 |
5 | 431 | 0.427885 | 2.69 | 千代岡 |
6 | 360 | 0.428280 | 2.78 | 高井 |
7 | 196 | 0.4260623 | 2.25 | 青野 |
8 | 104 | 0.4263067 | 2.3 | 岡田 |
9 | 81 | 0.425133 | 3.18 | 松田 |
10 | 78 | 0.424466 | 1.8 | 河合 |
木材 ヤング率:6000 ポアソン比:0.4 鋼材 ヤング率:206000 ポアソン比:0.3 理論値:0.099
メッシュ | 要素数 | 変位 | 相対誤差 | 測定者 |
0.7 | 155266 | 0.0861 | -13.0 | 関合 |
0.8 | 138453 | 0.08349 | -15.7 | 岡田 |
0.9 | 82766 | 0.083312 | -15.9 | 松田 |
1.2 | 32279 | 0.08357 | -15.6 | 青野 |
1.3 | 28343 | 0.083667 | -15.5 | 山口 |
1.4 | 23667 | 0.08368 | -15.5 | 山本 |
1.5 | 19958 | 0.083516 | -15.6 | 進藤 |
1.6 | 19451 | 0.08249 | -13.1 | 河合 |
1.8 | 10933 | 0.084021 | -15.1 | 山口 |
2 | 10764 | 0.083323 | -15.8 | 進藤 |
3 | 3618 | 0.083467 | -15.7 | 山本 |
4 | 1623 | 0.0852 | -13.9 | 関合 |
5 | 1007 | 0.083104 | -16.1 | 千代岡 |
6 | 842 | 0.0821 | -17.1 | 高井 |
7 | 554 | 0.08075 | -18.4 | 青野 |
8 | 289 | 0.079715 | -19.5 | 進藤 |
9 | 261 | 0.078427 | -20.8 | 松田 |
10 | 232 | 0.082495 | -16.7 | 河合 |
木材のみ(異方性一次)と、鋼材と木材(サンドウィッチ、等方性(鋼材)と異方性二次(木材))のグラフ