- 追加された行はこの色です。
- 削除された行はこの色です。
#author("2023-06-07T17:54:34+09:00","default:kouzouken","kouzouken")
#author("2023-06-09T17:45:46+09:00","default:kouzouken","kouzouken")
[[構造研ゼミ]]
#contents
*卒論研究 [#j33b34d8]
**サンドイッチ梁の解析 [#tb32f611]
***一つのメッシュの場合の変位を求める [#z54984e3]
・理論式：$\delta=\frac{5q\ell^{4}}{384EI}$
q:単位荷重(N/mm)、$\ell$:スパン(mm)
異方性を考慮すると理論式は：$\delta=\frac{5q\ell^{4}}{384EI}+\frac{P\ell}{4kGA}$?
9mm鋼材のサンドイッチ梁
・木材寸法：90×450×10000
・鋼材寸法：9×430×10000
・載荷荷重：25800N
・スパン9700の単純支持
・理論値：13.1972mm
2次要素
,メッシュサイズ,切り方,要素数,変位｛mm｝,相対誤差{%}
,9〜100,細かい,607177,13.5669,2.80
,9〜100,普通,293389,13.5395,2.59
,9〜100,粗い,232008,13.5175,2.43
,9〜1000,細かい,605773,13.5673,2.80
,9〜1000,普通,293794,13.5420,2.61
,9〜1000,粗い,209663,13.5220,2.46
,18〜1000,普通,110325,13.5151,2.41
,27〜1000,普通,52160,13.4926,2.24
,36〜1000,普通,37109,13.4929,2.24
,45〜1000,普通,24756,13.4668,2.04
,90〜100,普通,12461,13.4160,1.66
,90〜1000,普通,10028,13.4173,1.67
,100〜100,普通,12147,13.3953,1.50
,180〜1000,普通,2707,13.3657,1.28
,450〜1000,普通,702,13.2718,0.57
6mm鋼材のサンドイッチ梁
・木材寸法：90×450×10000
・鋼材寸法：6×430×10000
・載荷荷重：25800N
・スパン9700の単純支持
・理論値：16.1275mm
1次要素
,メッシュサイズ,切り方,要素数,変位｛mm｝,相対誤差{%}
,6〜100,普通,538677,15.7837,-2.13
,6〜1000,細かい,833424,16.2415,0.71
,6〜1000,普通,539836,15.7817,-2.14
,6〜1000,粗い,253984,16.0605,-0.42
,12〜1000,細かい,365378,16.1589,0.19
,12〜1000,普通,194959,15.6039,-3.25
,12〜1000,粗い,165214,15.9880,-0.86
,18〜1000,普通,116422,16.0048,-0.76
,24〜1000,普通,67499,15.9587,-1.05
,60〜1000,普通,17108,15.3377,-4.90
,120〜1000,普通,5145,13.5410,-16.04
,600〜1000,普通,627,3.1687,-80.35
2次要素
,メッシュサイズ,切り方,要素数,変位｛mm｝,相対誤差{%}
,6〜100,普通,538677,,
,6〜100,普通,538677,16.4899,2.25
,6〜1000,細かい,833424,,
,6〜1000,普通,539836,,
,6〜1000,粗い,253984,,
,12〜1000,普通,194959,,
,18〜1000,普通,116422,,
,24〜1000,普通,67499,,
,60〜1000,普通,17108,,
,120〜1000,普通,5145,,
,600〜1000,普通,627,,
,12〜1000,普通,194959,16.4700,2.12
,18〜1000,普通,116422,16.4477,1.99
,24〜1000,普通,67499,16.4357,1.91
,60〜1000,普通,17108,16.3839,1.59
,120〜1000,普通,5145,16.3211,1.20
,450〜1000,普通,684,16.2041,0.47
,600〜1000,普通,627,16.1753,0.30
木材のみの梁
・木材寸法：90×450×10000
・載荷荷重：25800N
・スパン9700の単純支持
・理論値：29.0105 mm
2次要素
,メッシュサイズ,切り方,要素数,変位｛mm｝,相対誤差{%}
,6〜100,細かい,6306,30.1230,3.83
,6〜1000,普通,5504,30.0955,3.74
,6〜1000,細かい,1044,29.9494,3.24
,45〜45,普通,33538,30.1430,3.90
,90〜90,普通,7585,30.0992,3.75
,90〜1000,細かい,5513,30.0930,3.73
,90〜1000,普通,5504,30.0955,3.74
,180〜180,普通,1482,30.0268,3.50
,300〜1000,普通,414,29.8315,2.83
メッシュサイズ：最大1000、最小6、切り方：普通 で切った様子
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/okada/m1000m6.png
メッシュサイズ：最大1000、最小6、切り方：細かい で切った様子
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/okada/m1000m6koma.png
メッシュサイズ：最大100、最小6、切り方：細かい で切った様子
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/okada/mokum100m6koma.png
メッシュサイズ：最大1000、最小90、切り方：普通 で切った様子
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/okada/mokum1000m90.png
メッシュサイズ：最大1000、最小90、切り方：細かい で切った様子
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/okada/mokum1000m90koma.png
***短い梁を作ってメッシュ結合を試す [#hf1979e3]
・木材寸法：90×450×1000
・鋼材寸法：6×430×1000
載荷荷重：25800N
・木材と鋼材を別のメッシュで切った様子
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/okada/6koum9m1.png
木材と鋼材のメッシュの節点が一致していない。
だがこれでも解析はできた
・一つのメッシュで切った様子
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/okada/6koum6m100ittai.png
サンドイッチ梁の解析を行っていく
・木材寸法：90×450×10000
・鋼材寸法：6×430×10000,9×430×10000
載荷荷重：25800N
-[[LIAISON_MAILサイト:https://code-aster.org/forum2/viewtopic.php?id=19761]]
-[[結合メモ:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?cmd=read&page=salome%20%E7%B5%90%E5%90%88%E3%83%A1%E3%83%A2&word=%E7%B5%90%E5%90%88%E3%83%A1%E3%83%A2]]
-[[サンドイッチ梁2021:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?Salome-Meca%E6%BC%94%E7%BF%92_%E3%82%B5%E3%83%B3%E3%83%89%E3%82%A4%E3%83%83%E3%83%81%E6%A2%81%282021%29]]
-[[サンドイッチ梁:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?Salome-Meca%E6%BC%94%E7%BF%92_%E3%82%B5%E3%83%B3%E3%83%89%E3%82%A4%E3%83%83%E3%83%81%E6%A2%81]]
*日誌 [#k09dda84]
5/25 木のみの梁での解析を始めた
5/23 一つのメッシュでの解析を進めた
5/17 節点を近づけることで時間はかかるが結合したメッシュで解析ができた。
5/16 解析がうまく行かないのはメッシュの節点共有ができていないためだと思われる
5/12 宮地エンジニアリングさんの話を聞くことができた。
5/11 ２つのメッシュを結合するのがうまくできなかった。モデルやメッシュを変えて色々試してみる
5/10 木材と鋼材を別のメッシュで切って解析する方法を及川さんから教えてもらった。
この方法を試していく
5/9 木と鋼の梁の解析を進めた
5/8 木と鋼のサンドイッチ梁を作製した。
5/2 木のみ梁の解析を行った。
5/1 木のみのサンドイッチ梁をsalomeで作成した
4/25 枝豆を植えた
4/24 春休み課題のデータ追加
*春休み課題 [#f313ac2e]
座屈
***10×5×100 木材 等方性 片持梁 2次モード 載荷荷重：3000N [#o6faa051]
理論値：弱軸方向 1233.70N
,メッシュサイズ,要素数,弱軸方向｛N｝,相対誤差{%}
,0.4,420533,1398.723,13.3
,0.5,151096,1429.68,15.9
,0.6,101211,1450.239,17.6
,0.7,93797,1451.622,17.7
,0.8,60625,1467.732,19.0
,0.9,21520,1611.57,31.0
,1.0,19543,1604.212,30.0
,1.5,9252,1682.265,36.3
,2.0,2928,2445.381,98.2
理論値：強軸方向 4934.8N
,メッシュサイズ,要素数,強軸方向｛N｝,相対誤差{%}
,0.4,420533,5206.86,5.5
,0.5,151096,5232.48,6.0
,0.6,101211,5248.89,6.4
,0.7,93797,5250.48,6.4
,0.8,60625,5268,6.8
,0.9,21520,5419.1,9.8
,1.0,19543,5419.96,9.8
,1.5,9252,5509.95,11.7
,2.0,2928,5946.57,20.5
***10×5×100 木材 等方性 片持梁 １次モード 載荷荷重：3000N [#m56008ce]
理論値：弱軸方向 154.212N
,メッシュサイズ,要素数,弱軸方向｛N｝,相対誤差{%}
,0.4,420533,157.91,2.4
,0.5,151096,161.48,4.7
,0.6,101211,163.82,6.2
,0.7,93797,163.97,6.3
,0.8,60625,165.61,7.4
,0.9,21520,181.96,18.0
,1.0,19543,181.25,17.5
,1.5,9252,190.28,23.4
,2.0,2928,280.38,81.8
理論値：強軸方向 616.85N
,メッシュサイズ,要素数,強軸方向｛N｝,相対誤差{%}
,0.4,420533,617.95,0.2
,0.5,151096,621.02,0.7
,0.6,101211,623.02,1.0
,0.7,93797,623.28,1.0
,0.8,60625,625.10,1.3
,0.9,21520,643.60,4.3
,1.0,19543,644.12,4.4
,1.5,9252,654.74,6.1
,2.0,2928,710.74,15.2
-[[座屈解析:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?%E5%BA%A7%E5%B1%88%E8%A7%A3%E6%9E%90]]
-[[座屈資料:http://wwwra.meijo-u.ac.jp/labs/ra007/murata/pdf/textbook/hari_14.pdf]]
-[[salome-mecaサイト:https://code-aster.org/V2/doc/default/en/index.php?man=commande]]
-[[AFFE-CHAR-MECA資料:https://code-aster.org/V2/doc/default/en/man_u/u4/u4.44.01.pdf]]
-[[座屈解析のコマンドファイル:http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/tamura/buckling_TR1.comm]]
*実習内容 [#d1d9af6d]
***12/2課題 [#o1b9f35f]
サンドイッチ梁の２次要素の解析（理論値：0.099mm）
,メッシュ長さ,要素数,変位(mm),相対誤差(%),計算者
,0.5,-,-,-,千代岡
,0.6,214117,-,-,高井
,0.7,155266,0.0861,-13.0,関合
,0.8,138453,0.08348657,-15.7,岡田
,0.9,82766,0.08331204,-15.85,松田
,1.2,32279,0.08357352,-15.6,青野
,1.3,28343,0.0836677,-15.49,山口
,1.4,23667,0.083679597,-15.48,山本
,1.5,19958,0.083515736,-15.6,進藤
,1.6,19451,0.0860365,-13.1,河合
,1.8,10933,0.08402172,-15.13,山口
,2,10764,0.08332360,-15.8,進藤
,3,3618,0.083497165,-15.66,山本
,4,1623,0.0852,-13.9,関合
,5,1007,0.083104,-16.1,千代岡
,6,842,0.0821,-17.1,高井
,7,554,0.080750428,-18.4,青野
,8,289,0.07971546,-19.5,岡田
,9,261,0.07842736,-20.78,松田
,10,232,0.082495,-16.67,河合
各要素数と変位のグラフ（直線は理論値）
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/okada/sandoittikadai.png
サンドイッチ梁と異方性1次要素との比較
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/okada/sandoittihikaku.png
***11/25課題 [#b8283d8b]
等方性２次要素のデータ（理論値：0.4167mm）
,メッシュ長さ,要素数,変位(mm),相対誤差(%),計算者
,0.5,604167,-,-,千代岡
,0.6,203209,0.423827,1.7,高井
,0.7,145234,0.43011,3.22,関合
,0.8,140987,0.43005836,3.2,岡田
,0.9,91974,0.429911921,3.18,松田
,1.2,24800,0.429776978,3.14,青野
,1.3,23132,0.429885619992,3.16,山口
,1.4,17617,0.429745386435,3.1,山本
,1.5,15433,0.429844,3.2,進藤
,1.6,15900,0.429754,3.13,河合
,1.8,11677,0.429623532,3.10,山口
,2,10460,0.429605,3.1,進藤
,3,2436,0.429216538,3.0,山本
,4,1453,0.4293,3.02,関合
,5,431,0.427885,2.69,千代岡
,6,360,0.4282,2.78,高井
,7,196,0.426062273,2.25,青野
,8,104,0.4263067,2.3,岡田
,9,81,0.425133059,2.03,松田
,10,78,0.424466,1.8,河合
・等方性２次要素と等方性１次要素のグラフの比較（直線は理論値）
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/okada/niziyousokadai.png
異方性１次要素のデータ(理論値:0.4917mm)
,メッシュ長さ,要素数,変位(mm),相対誤差(%),計算者
,0.5,604167,0.50919,3.56,千代岡
,0.6,203209,0.504716,2.6,高井
,0.7,145234,0.5036,2.42,関合
,0.8,140987,0.50282705,2.3,岡田
,0.9,91974,0.500527681,1.8,松田
,1.2,24800,0.487393322,-0.9,青野
,1.3,23132,0.4884103968,-0.67,山口
,1.4,17617,0.484032743017,-1.56,山本
,1.5,15433,0.482022,-2.0,進藤
,1.6,15900,0.4832858,-1.7,河合
,1.8,11677,0.4785524135,-2.67,山口
,2,10460,0.479058,-2.6,進藤
,3,2436,0.427868847826,-12.98,山本
,4,1453,0.42772,-13.02,関合
,5,431,0.273640,-44.3,千代岡
,6,360,0.3392699,-31.0,高井
,7,196,0.21362825,-58.5,青野
,8,104,0.22574,-54.1,岡田
,9,81,0.2275024,-53.73,松田
,10,78,0.203271,-58.7,河合
・異方性１次要素と等方性１次要素のグラフの比較（直線は理論値）
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/okada/ihouseikadai.png
***11/25 [#hd7089c2]
11/18の課題：各要素数と変位のグラフ、直線は理論値
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/okada/tanjunkadai.png
***11/24 [#efb725bc]
11/18の課題の表
,メッシュ長さ,要素数,変位(mm),相対誤差(%),計算者
,0.5,604167,0.428982,2.94,千代岡
,0.6,361584,0.421233,1.09,高井
,0.7,145234,0.4225,1.4,関合
,0.8,140987,0.422627385,1.4,岡田
,0.9,91857,0.420351606,0.88,松田
,1.2,24520,0.404744325,-2.87,青野
,1.3,23132,0.4045,-2.93,山口
,1.4,4518,0.3986,-4.34,山本
,1.5,15433,0.396317756757,-4.9,進藤
,1.6,15900,0.399049,-4.24,河合
,1.8,11677,0.404457,-2.9,山口
,2,10460,0.394818715517,-5.3,進藤
,3,734,0.32447,-22.13,山本
,4,1453,0.3329,-20.1,関合
,5,431,0.136240,-67.3,千代岡
,6,360,0.2130486,-48.9,高井
,7,196,0.1019892,-75.5,青野
,8,104,0.1158624,-72.2,岡田
,9,81,0.1247076,-70.1,松田
,10,78,0.07733,-81.4,河合
***11/18 [#ia1aa2cc]
11/18の課題の値
,メッシュの長さ,要素数,変位,相対誤差,計算者
,8,104,0.1158624,-72.2,岡田
,0.8,140987,0.42262738527,-1.4,岡田
***11/17 [#jd8dcdf3]
11/11の課題：各要素数と変位のグラフ
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/okada/heni11.png
理論値との比較
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/okada/katamotihikaku.png
***11/16 [#w896225a]
11/11の課題の表
,メッシュの長さ,要素数,変位,相対誤差,計算者
,0.5,59504,6.56,-1.5,千代岡
,0.6,45512,6.48774575,-2.69,高井
,0.7,39075,6.54132,-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.6458525,-15.3,進藤
,3,667,5.4053975,-18.96,山本
,4,264,3.6161,-45.8,関合
,5,191,3.86,-42,千代岡
,6,190,2.5077325,-62.4,高井
,7,75,1.41225,-78.8,青野
,8,56,1.2887175,-80.7,岡田
,9,49,1.28799,-80.7,松田
,10,44,1.226075,-81.6,河合
***11/14 [#y0efd888]
11/11の課題の値
メッシュの長さ：8 / 0.8
要素数：56 / 13397
変位：1.2887175 / 6.43695
相対誤差：80.7% / 3.5%
,メッシュの長さ,要素数,変位,相対誤差,計算者
,8,56,1.2887175,-80.7,岡田
,0.8,13397,6.43695,-3.5,岡田
***11/4課題 [#j9ca3c5e]
秋田市の10月日平均湿度２週間分
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/okada/situdo.png
*試験書き込み [#k4ff76c6]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/okada/bbb.png
*メモ [#w5f3def5]
**texメモ [#e6e3c447]
https://www.str.ce.akita-u.ac.jp/~gotou/tebiki/texhtml.html#sono4
式の書き方は$で挟んで
$\delta=\frac{P\ell^{3}}{3EI}$
**UNIXコマンド [#m9b3d34b]
-[[ここのUNIXコマンド:https://www.str.ce.akita-u.ac.jp/~gotou/linux/vine.html#unix]]
**gnuplot [#cb967f34]
・x軸の範囲を決める set xrange[0:150]
・0軸を表示 set yzeroaxis 、 やめるときは unset yzeroaxis
・凡例を非表示 unset key
・ズーム ctrl+マウスホイール
-[[gnuplotの使い方:https://www.str.ce.akita-u.ac.jp/~gotou/linux/vine.html#gnuplot]]
-[[viの練習:http://www.str.ce.akita-u.ac.jp/~gotou/linux/vine.html#vi]]
![[PukiWiki]](image/pukiwiki.png "[PukiWiki]")
