#author("2024-04-25T15:50:46+09:00","default:kouzouken","kouzouken")
#author("2024-04-25T18:40:57+09:00","default:kouzouken","kouzouken")
&aname(TOP);
//CENTER:&size(50){あなたが閲覧を開始した時刻は&_time;};
//CENTER:&size(15){総閲覧数&counter;};
CENTER:&size(15){現在の閲覧数&online;};
#contents

*2/28 [[&ruby(ざくつかいせきみほん){座屈解析見本};:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?Salome%E3%83%A1%E3%83%A2#m8182d1b]][#h3e69f89]
[[2023先輩見本:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?2023+%E6%98%A5%E8%AA%B2%E9%A1%8C%E5%BA%A7%E5%B1%88%E8%A7%A3%E6%9E%90]]~
[[code-asterの説明書:http://dugi6514.odns.fr/doc/v15/en/index.php?man=commande]] (英語,フランス語で書かれているので翻訳サイトを使って読んでください.)~
[[code-aster質問ページ:https://forum.code-aster.org/public/]](説明書読んでも分からなかったら,このページで調べてみるといいかも.ただし,英語でキーワード等を入力する必要があります。)~
・鋼材の寸法 10-1-100 (mm) 荷重1000N 面載荷100N/mm^2で解析成功した。
https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?cmd=edit&page=%EF%BC%92%EF%BC%90%EF%BC%92%EF%BC%94%E5%BA%A7%E5%B1%88%E7%8F%AD
*3/1 [[座屈について:https://www.str.ce.akita-u.ac.jp/~gotou/kouzou/kouzou2_15.html]]  [#x8dc6dfe]
・[[&size(25){解析の手順};:https://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/tamura/buckling_TR1.comm]]は確認できた。🙆 &size(25){[[画像つき:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?%E5%BA%A7%E5%B1%88%E8%A7%A3%E6%9E%90]]はこっち};


*3/27 [[座屈内容の決定:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?%E6%A7%8B%E9%80%A0%E7%A0%94%E3%82%BC%E3%83%9F#pef90304]] [#c5ae9408]
・片持ち梁、両端固定、単純梁で座屈解析
 
→1、2、%%3次%%モードで比較

  1次、2次要素で比較

  強軸、弱軸で比較

 片持ち梁は座屈荷重の計算に注意する

***座屈荷重理論式(オイラー座屈) [#c316945c]
P=(\frac{π}{k\ell})^{2}EI$ 

&size(20){片持ちの場合:''k=2''};

&size(20){両端固定の場合:''k=0.5''};

&size(20){両端ピン固定の場合:''k=1''};

&size(20){固定+ピン固定の場合:''k=0.7''};

*3/28 [#r0c51797]
**片持ち梁 [#r3f824a8]
10mm×5mm×100mm、P=50000N

一次モード 

&size(25){理論値 弱軸方向:5294.63N 強軸方向:21178.53N};
,メッシュサイズ,要素数,弱軸方向[N],相対誤差[%],強軸方向[N],相対誤差[%]
,0.7,93797,5541.15,4.656,21278.95,0.472
,0.8,60625,5600.65,5.780,21338.55,0.756
,0.9,21520,6031.7,13.92,21821.50,3.036
,1.0,19543,6015,13.6,21826.75,3.061
,2,2928,8510.5,60.738,23498.25,10.953
,3,1379,8212.6,55.112,24202.95,14.280

二次モード

&size(25){理論値 弱軸方向:47651.68N 強軸方向:190606.74N};
,メッシュサイズ,要素数,弱軸方向[N],相対誤差[%],強軸方向[N],相対誤差[%]
,0.7,93797,49088.45,3.015,179847,-5.645
,0.8,60625,49646.7,4.187,184380,-3.267
,0.9,21520,53452,12.172,184380,-3.267
,1.0,19543,53268.5,11.787,184354,-3.280
,2,2928,74587,56.525,197805.5,3.777
,3,1379,72521.5,52.191,234137,22.838
二次要素

一次モード

&size(25){理論値 弱軸方向:5294.63N 強軸方向:21178.53N};
,メッシュサイズ,要素数,弱軸方向[N],相対誤差[%],強軸方向[N],相対誤差[%]
,0.7,93797,5308.6,0.264,21067.75,-0.523
,0.8,60625,5308.75,0.267,21068.15,-0.521
,0.9,21520,5309.1,0.273,21069.35,-0.516
,1.0,19543,5309.25,0.276,21070.3,-0.511
,2,2928,5312.35,0.335,21081.60,-0.458
,3,1379,5315.40,0.392,21089.3,-0.421


二次モード

&size(25){理論値 弱軸方向:47651.68N 強軸方向:190606.74N};
,メッシュサイズ,要素数,弱軸方向[N],相対誤差[%],強軸方向[N],相対誤差[%]
,0.7,93797,47034.5,-1.295,178077,-6.574
,0.8,60625,47035.95,-1.292,178080.5,-6.572
,0.9,21520,47039.7,-1.284,178092,-6.566
,1.0,19543,47041.25,-1.281,178099.5,-6.562
,2,2928,47077.3,-1.205,178261.5,-6.477
,3,1379,47110.8,-1.135,178282,-6.466

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/2024zakutu/katamoti_1jii.png~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/2024zakutu/katamoti_2ji.png



*4/6 [#k3daebaa]
**片方固定、他端回転 [#l861a507]
10×5×110 荷重:50000N(0,0,100)と(10,0,100)を結ぶ直線をピン固定して解析した。
***一次要素メッシュ [#nee67f66]
一次モード 
&size(25){理論値 弱軸方向:35720.2N 強軸方向:142880.9N};
,メッシュサイズ,要素数,弱軸方向[N],相対誤差[%],強軸方向[N],相対誤差[%]
,0.7,80513,32400.0,-0.093,118800,-0.169
,0.8,61168,32650,-0.086,119000,-0.167
,0.9,52074,32800,-0.082,119100,-0.166
,1.0,20256,35400,-0.090,121550,-0.149
,2.0,3102,40789.4,-0.142,130050.0,-0.090
,3.0,3102,67345.0,0.885,136750.0,-0.043
二次モード
&size(25){理論値 弱軸方向:142880.9N 強軸方向:571523.7N};
,メッシュサイズ,要素数,弱軸方向[N],相対誤差[%],強軸方向[N],相対誤差[%]
,0.7,80513,110450,-0.227,198900,-0.652
,0.8,61168,111350,-0.221,198150,-0.653
,0.9,52074,111750,-0.218,198900,-0.652
,1.0,20256,120400,-0.157,213800,-0.626
,2.0,3102,167500.0,-0.172,-,-
,3.0,1259,229175.0,-0.604,-,-676

***二次要素メッシュ [#z90a48c8]
一次モード 
&size(25){理論値 弱軸方向:35720.2N 強軸方向:142880.9N};
,メッシュサイズ,要素数,弱軸方向[N],相対誤差[%],強軸方向[N],相対誤差[%]
,0.7,80513,30830.2,-0.137,117500,-0.178
,0.8,61168,32638.2,-0.083,119015.5,-0.167
,0.9,52074,30832.3,0.137,117487.5,-0.178
,1.0,20256,30836.0,-0.137,117479,-0.178
,2.0,3102,30860.0,0.136,105310.0,-0.263
,3.0,1259,30915.0,-0.135,136795.0,-0.043


二次モード
&size(25){理論値 弱軸方向:142880.9N 強軸方向:571523.7N};
,メッシュサイズ,要素数,弱軸方向[N],相対誤差[%],強軸方向[N],相対誤差[%]
,0.7,80513,105156,-0.26,187217,-0.672
,0.8,61168,111326,-0.22,198137.5,-0.653
,0.9,52074,105164,-0.264,187238.5,-0.672
,1.0,20256,105180.5,-0.264,187285.5,-0.672
,2.0,3102,117430.0,-0.178,187650.0,-0.672
,3.0,1259,117255.0,-0.179,188695.0,-0.670

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/2024zakutu/tatankotei_1ji.png~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/2024zakutu/tatanktoei_2ji.png

**単純梁 [#z651b584]
&aname(N);10mm×5mm×100mm、P=50000N
梁の両端に支承があるものとして計算~
載荷面側にローター支承か、固定かで解析値がかわる~
支承の位置でも解析が回るかどうかが違う

&size(25){一次モード理論値 強軸84714[N] 弱軸21179[N]};

&size(25){二次モード理論値 強軸338856[N] 弱軸84714[N]};

***1次要素 [#kd155bde]
一次モード
,メッシュサイズ,要素数,弱軸方向[N],相対誤差[%],強軸方向[N],相対誤差[%]
,0.7,93797,21898,3.395,81855,3.375
,0.8,60625,22141,4.542,82224,2.939
,0.9,21520,23841,12.57,83786,1.095
,1.0,19543,23781,12.29,83750,1.138
,2.0,2928,33638,58.83,90036,6.282
,3.0,1379,32035,51.26,91687,8.231

二次モード
,メッシュサイズ,要素数,弱軸方向[N],相対誤差[%],強軸方向[N],相対誤差[%]
,0.7,93797,85369,0.773,293753,13.31
,0.8,60625,86179,1.729,312977,7.637
,0.9,21520,92916,9.682,315815,6.800
,1.0,19543,92782,9.524,315478,6.899
,2.0,2928,129723,53.13,320761,5.340
,3.0,1379,125641,48.31,330941,2.336

***2次要素 [#q9032b40]
メッシュサイズ0.7、0.8の解析は&color(black,lightblue){荷重50000Nではうまく行かなかった};が、荷重を半分に減らした結果解析が回ったのでこのメッシュ長さに限り
//アンカー名の部分に半角数字入れると機能しなかった
&aname(N){''25000N''};での解析結果を示す。(上記の解析と異なり、10mm*5mm*120mmの単純梁で解析を行った。)~

一次モード
,メッシュサイズ,要素数,弱軸方向[N],相対誤差[%],強軸方向[N],相対誤差[%]
,0.7,93797,,,,
,0.8,60625,,,,
,0.9,21520,20989,0.892,79224,6.481
,1.0,19543,20990,0.888,79129,6.593
,2.0,2928,20999,0.850,81044,4.332
,3.0,1379,21004,0.826,81334,3.990

,メッシュサイズ,要素数,弱軸方向[N],相対誤差[%],強軸方向[N],相対誤差[%]
,0.7,80310,17060.63,-0.160,101479.0,-0.725
,0.8,66909,17060.8,-0.160,101580.5,-0.727

二次モード
,メッシュサイズ,要素数,弱軸方向[N],相対誤差[%],強軸方向[N],相対誤差[%]
,0.7,,,,,
,0.8,,,,,
,0.9,21520,81641,3.627,326113,3.761
,1.0,19543,81673,3.590,335961,0.854
,2.0,2928,81885,3.339,307703,9.194
,3.0,1379,81978,3.230,279647,17.47

,メッシュサイズ,要素数,弱軸方向[N],相対誤差[%],強軸方向[N],相対誤差[%]
,0.7,93797,67002.5,0.209,306595.0,-0.095
,0.8,60625,67004.3,-0.209,306960.0,-0.094

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/2024zakutu/1jilabel.png~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/2024zakutu/2jilabel.png

**両端固定はできませんでした [#z376197b]
salome-mecaで両端を固定した状態で座屈させる設定ができなかった~
片面をdx、dy、dzのすべてを固定して、もう片面をdzを固定しないdxとdyだけ固定させるとdzを固定しない面が回転してしまう~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/2024zakutu/2ji.png~
弱軸方向の座屈

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/2024zakutu/k2ji.png~
強軸方向の座屈

中心の軸で回転して両端固定の状態にならない

**4/15 両端固定の再現 [#xff94727]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/2024zakutu/ryoutankotei50000000.png~
端に固定端のような硬い部材をつけて両端固定のようなモデルの再現をしてみました~
固定端の部材が若干伸びている~
→固定端部材のヤング率を大きくするとエラー 写真は50000GPaの固定端~
 さらに100倍にするとエラー~
理論値とずれが大きい 2.3%程度

*4/10 [#n6f11ed4]

*4/17 発表 [#d878784c]
[[座屈発表スライド:http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/2024zakutu/zakutuhappyou1.pdf]]


*4/23 [#hf7a46b0]
一回目の発表では強軸方向で理論値から離れるという結果だったので、10mm*5mmの長方形断面から短辺と長辺の比を変えて解析した。弱軸方向と強軸方向の値を理論値と比較した。
**片方固定他端ローラー [#f8c9e4d0]
片方固定他端ローラーでは強軸歩行の値が算出されないことが多かった。charcritの範囲を変えても変わらなく解析が重くなったので今回は残り時間を考慮して算出できなタッカ部分は今後のかだいとする。
***5mm*5mm断面(正方形) [#g6390354]
一次要素メッシュ [#nee67f66] 
一次モード 
&size(25){理論値 弱軸方向:17860.1N 強軸方向:17860N};
,メッシュサイズ,要素数,弱軸方向[N],相対誤差[%],強軸方向[N],相対誤差[%]
,0.9,25149,16292.5,,16350.0,
,1.0,10300,17535.0,,18222.5,
,2.0,1486,23247.5,,23552.5,
,3.0,674,25992.5,,36240.0,
,4.0,250,42745.0,,73812.5,

二次モード
&size(25){理論値 弱軸方向:35720.2N 強軸方向:35720.2N};
,メッシュサイズ,要素数,弱軸方向[N],相対誤差[%],強軸方向[N],相対誤差[%]
,0.9,25149,55582.5,,55710.0,
,1.0,10300,59687.5,,62165.0,
,2.0,1486,78602.5,,79872.5,
,3.0,674,88432.5,,120920.0,
,4.0,250,-,-,-,-

***5mm*15mm断面 [#j2c25a18]
&size(25){理論値 弱軸方向:53580.3N 強軸方向:482223.1N};
一次モード
,メッシュサイズ,要素数,弱軸方向[N],相対誤差[%],強軸方向[N],相対誤差[%]
,0.9,71221,49452.5,,169067.5,
,1.0,25247,54552.5,,185825.0,
,2.0,4073,74115.0,,250157.5,
,3.0,2185,87115.0,,285575.0,
,4.0,676,115342.5,,392425.0,

二次モード
&size(25){理論値 弱軸方向:107160.7N 強軸方向:964446.2N};
,メッシュサイズ,要素数,弱軸方向[N],相対誤差[%],強軸方向[N],相対誤差[%]
,0.9,71221,372075.0,,-,-
,1.0,25247,376000.0,,-,-
,2.0,4073,390742.0,,-,-
,3.0,2185,391830.0,,-,-
,4.0,676,396200.0,,-,-

***5mm*25mm断面 [#va92eb23]
&size(25){理論値 弱軸方向:89300.5N 強軸方向:2232514.6N};
一次モード
,メッシュサイズ,要素数,弱軸方向[N],相対誤差[%],強軸方向[N],相対誤差[%]
,0.9,124415,83060.0,,1617770.0,
,1.0,42396,91720.0,,1421210.0,
,2.0,6810,124760.0,,1437090.0,
,3.0,3651,147040.0,,1440030.0,
,4.0,1196,213850.0,,1443810.0,

&size(25){理論値 弱軸方向:178601.1N 強軸方向:4465029.1N};
二次モード
,メッシュサイズ,要素数,弱軸方向[N],相対誤差[%],強軸方向[N],相対誤差[%]
,0.9,124415,285170.0,,1774480.0,
,1.0,42396,313440.0,,1799040.0,
,2.0,6810,423370.0,,2385610.0,
,3.0,3651,485610.0,,2594950.0,
,4.0,1196,704680.0,,3598800.0,

***5mm*35mm断面 [#l6cd7f14]
&size(25){理論値 弱軸方向:125020.8N 強軸方向:6126019.9N};
一次モード
,メッシュサイズ,要素数,弱軸方向[N],相対誤差[%],強軸方向[N],相対誤差[%]
,0.9,135830,119035.0,,3117310.0,
,1.0,57832,130042.5,,3117310.0,
,2.0,9521,179165.0,,3134985.0,
,3.0,5132,196770.0,,3134232.5,
,4.0,1732,312672.0,,3172595.0,

&size(25){理論値 弱軸方向:250041.6N 強軸方向:12252040N};
二次モード
,メッシュサイズ,要素数,弱軸方向[N],相対誤差[%],強軸方向[N],相対誤差[%]
,0.9,135830,409972.5,,6586947.5,
,1.0,57832,44040.0,,6799975.0,
,2.0,9521,605692.5,,6711985.0,
,3.0,5132,653047.5,,6750327.5,
,4.0,1732,1046115.0,,6788845.0,

***5mm*50mm断面 [#d546789e]
&size(25){理論値 弱軸方向:178601.1N 強軸方向:17860116.0N};
一次モード
,メッシュサイズ,要素数,弱軸方向[N],相対誤差[%],強軸方向[N],相対誤差[%]
,0.9,187141,172580.0,,6442240.0,
,1.0,170937,172680.0,,644420.0,
,2.0,28969,199600.0,,6456200.0,
,3.0,7589,279380.0,,6469460.0,
,4.0,4715,299360.0,,6507700.0,

&size(25){理論値 弱軸方向:357202.3N 強軸方向:35720233.0N};
二次モード
,メッシュサイズ,要素数,弱軸方向[N],相対誤差[%],強軸方向[N],相対誤差[%]
,0.9,187141,594140.0,,-,-
,1.0,170937,594800.0,,-,-
,2.0,28969,686580.0,,-,-
,3.0,7589,935280.0,,-,-
,4.0,4715,1019320.0,,-,-


**単純梁 [#r3739e5a]
***5mm*5mm*120mm [#mc3b76d0]
一次モード

理論値 7353.65N 
,メッシュサイズ,要素数,弱軸方向[N],相対誤差[%],強軸方向[N],相対誤差[%]
,0.7,38975,8966.08,21.9,8966.08,21.9
,0.8,32158,9043.78,23.0,9043.78,23.0
,0.9,21664,9203.2,25.2,9203.2,25.2
,1.0,25276,9354.95,27.3,9354.95,27.3
,2.0,2027,11635.55,58.2,11635.55,58.2
,3.0,707,17525.83,138,17525.83,138
,4.0,183,26927.75,266,26927.75,266

二次モード

理論値  29414.62N 
,メッシュサイズ,要素数,弱軸方向[N],相対誤差[%],強軸方向[N],相対誤差[%]
,0.7,38975,35163.75,19.5,35163.75,19.5
,0.8,32158,35484.5,20.6,35484.5,20.6
,0.9,21664,36136,22.9,36136,22.9
,1.0,25276,36519.5,24.2,36519.5,24.2
,2.0,2027,45826.5,55.8,45826.5,55.8
,3.0,707,66233.25,125,66233.25,125
,4.0,183,101679,246,101679,246

***10mm*5mm*120mm [#j6545846]
一次モード

理論値 弱軸方向:14707.31N 強軸方向:58829N
,メッシュサイズ,要素数,弱軸方向[N],相対誤差[%],強軸方向[N],相対誤差[%]
,0.7,83010,17919.63,21.8,17919.63,21.8
,0.8,66909,18107.23,23.1,18107.23,23.1
,0.9,34242,18822.1,28.0,18822.1,28.0
,1.0,28514,18971.85,29.0,18971.85,29.0
,2.0,4179,24348.48,65.6,24348.48,65.6
,3.0,1144,43660,197,43660,197
,4.0,450,49779.5,238,49779.5,238


二次モード

理論値 弱軸方向:58829.24 強軸方向:235317N
,メッシュサイズ,要素数,弱軸方向[N],相対誤差[%],強軸方向[N],相対誤差[%]
,0.7,83010,70365.25,19.6,70365.25,19.6
,0.8,66909,71117.25,20.9,71117.25,20.9
,0.9,34242,73924.25,25.7,73924.25,25.7
,1.0,28514,74569.5,26.8,74569.5,26.8
,2.0,4179,95688.5,62.7,95688.5,62.7
,3.0,1144,170895,190,170895,190
,4.0,450,194382.75,230,194382.75,230


***15mm*5mm*120mm [#mb196dc3]
一次モード

理論値 弱軸方向:22060.96N 強軸方向:198549N
,メッシュサイズ,要素数,弱軸方向[N],相対誤差[%],強軸方向[N],相対誤差[%]
,0.7,112115,27137.5,23.0,,
,0.8,97171,27234.75,23.5,,
,0.9,70238,28016.25,27.0,146123,
,1.0,57841,28453,29.0,146285,
,2.0,8984,35628.75,61.5,149439.5,
,3.0,2337,73692.25,234,157327,
,4.0,1268,56624.5,156,162155,

二次モード
理論値 弱軸方向:88243.86N 強軸方向:794195N
,メッシュサイズ,要素数,弱軸方向[N],相対誤差[%],強軸方向[N],相対誤差[%]
,0.7,,,,,
,0.8,,,,,
,0.9,,,,,
,1.0,,,,,
,2.0,,,,,
,3.0,,,,,
,4.0,,,,,


***25mm*5mm*120mm [#mb196dc3]
一次モード

理論値 弱軸方向:36768.27 強軸方向:919207N
,メッシュサイズ,要素数,弱軸方向[N],相対誤差[%],強軸方向[N],相対誤差[%]
,0.7,172074,45497.5,23.7,,
,0.8,162241,45531.75,23.8,,
,0.9,114731,48094.5,30.8,346767,
,1.0,92806,47571,29.4,371403,
,2.0,14982,58169.5,58.2,362647,
,3.0,4067,,,371230,
,4.0,2353,,,375535.5,

二次モード
理論値 弱軸方向: 強軸方向:3676827N
,メッシュサイズ,要素数,弱軸方向[N],相対誤差[%],強軸方向[N],相対誤差[%]
,0.7,,,,,
,0.8,,,,,
,0.9,,,,,
,1.0,,,,,
,2.0,,,,,
,3.0,,,,,
,4.0,,,,,

**片持ち梁 [#k10af243]
***5mm×5mm [#e1645e62]
一次モード

理論値 2647.32N
,メッシュサイズ,要素数,弱軸方向[N],相対誤差[%],強軸方向[N],相対誤差[%]
,0.7,47584,2751.17,3.92,2766.46,4.5
,1.0,11073,2983.06,12.7,29.8844,12.9

二次モード

理論値 23825.84N
,メッシュサイズ,要素数,弱軸方向[N],相対誤差[%],強軸方向[N],相対誤差[%]
,0.7,47584,24350.3,2.2,24481.95,2.75
,1.0,11073,26402.1,10.8,26492,11.2



***15mm×5mm [#pb646a34]
一次モード

理論値 弱軸方向:7941.9N  強軸方向:71477.53N
,メッシュサイズ,要素数,弱軸方向[N],相対誤差[%],強軸方向[N],相対誤差[%]
,0.7,134082,8344.35,5.06,70772,-0.987
,1.0,29217,9094.5,14.5,71497.5,0.028

二次モード

理論値 弱軸方向:71477.53N 強軸方向:563235N
,メッシュサイズ,要素数,弱軸方向[N],相対誤差[%],強軸方向[N],相対誤差[%]
,0.7,134082,73998,3.53,557540,-10.5
,1.0,29217,80501,12.6,563235,-9.64



**両端固定 [#k5842209]
梁の別部分を1周分固定したがエラーが出てしまった~
回転を抑えるために載荷面の回転軸に並行な1辺を固定してみたがエラーだった

~
~
うまくいかなかったので固定のための部材を薄くした~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/2024zakutu/u-hikaku.png~
まず固定のための部材をヤング率50000MPaでやってみたら図のようになった~
部材のサイズは載荷面と同じ面積で、長さ0.1mmとした~

次に固定する部材のヤング率を500000MPaと10倍にした~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/2024zakutu/u10-hikaku.png~

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/2024zakutu/u10all-hikaku.png~
ほとんど面が回転していないことがわかった~

~
~
固定する部材と梁を繋いだときと離したときの解析結果の違いについて~
両端固定するのでこの比較は意味がないかもしれないです~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/2024zakutu/b-rhikaku.png~
固定の部材を繋いだものと離したものの比較~

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/2024zakutu/r-hikaku.png~
繋いだ結果~

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/2024zakutu/b-hikaku.png~
離した結果~

***ヤング率50000GPa [#m3a5fb9d]
1次モード~
,メッシュサイズ,要素数,弱軸方向[N],相対誤差[%],強軸方向[N],相対誤差[%]
,1.0,,93022,9.400,311247,9.296

2次モード~
,メッシュサイズ,要素数,弱軸方向[N],相対誤差[%],強軸方向[N],相対誤差[%]
,1.0,,339940,0.016,940560,31.24

***ヤング率50000×10GPa [#ye56b4be]
1次モード~
,メッシュサイズ,要素数,弱軸方向[N],相対誤差[%],強軸方向[N],相対誤差[%]
,1.0,,92677,9.807,307357,8.148

2次モード~
,メッシュサイズ,要素数,弱軸方向[N],相対誤差[%],強軸方向[N],相対誤差[%]
,1.0,,338802,0.320,931995,30.61

***5mm×10mm×100mm (固定部分20mm) [#uc4f9b09]
理論値 弱軸 1次モード 84714N/mm   2次モード 338856N/mm~
    強軸 1次モード 338856N/mm  2次モード 1355426N/mm~

1次モード~
,メッシュサイズ,要素数,弱軸方向[N],相対誤差[%],強軸方向[N],相対誤差[%]
,0.7,83010,88174,4.084,306331,9.598
,0.8,66909,88057,3.946,306226,9.629
,0.9,34242,88885,4.924,307272,9.321
,1.0,28514,89831,6.040,308059,9.089
,2.0,4179,104131,22.92,320988,5.273
,3.0,1144,133025,57.03,352599,4.056
,4.0,450,204274,141.1,,

2次モード~
,メッシュサイズ,要素数,弱軸方向[N],相対誤差[%],強軸方向[N],相対誤差[%]
,0.7,83010,322665.5,4.778,926110,31.67
,0.8,66909,322444,4.833,925885,31.69
,0.9,34242,325465,3.952,929185,31.45
,1.0,28514,329088,2.883,931170,31.30
,2.0,4179,378906,11.82,969335,28.48
,3.0,1144,476388,40.59,,
,4.0,450,,,,
*4/24 再発表 [#g77e5521]
//#vote(片持ち梁[2],単純梁[0],他端固定[0],両端固定[0])
拘束条件を見直し、再度1から解析をすることにした。

オイラーの座屈荷重の理論式は梁をモデルとしているので長さをもっと長く取れば理論式に近づく??

理論式に用いるLの値がスパン長なのか部材の長さなのかということを統一する。

片持ち梁は概ね理論式に近い値が出ていたので片方固定他端ローラーと単純梁の解析に絞る。

理論式について調べ、高次の座屈モードの座屈荷重について見直した。参考→http://www17.plala.or.jp/poppy06/study.html の14.長柱の座屈より


**片方固定他端ローラー 強軸方向 [#d73389c7]
***10mm*5mm*100mm 強軸方向の理論値 一次モード:173272.2N 二次モード:512215.8N [#tcdd7799]
,メッシュサイズ,要素数,一次モード[N],相対誤差[%],二次モード[N],相対誤差[%]
,0.8,60760,164180,5.24,434275,15.1
,0.9,21623,167655,3.24,443390,13.4
,1.0,19416,167870,3.12,443665,13.4
,2.0,2973,180035,3.90,473415,7.58
***10mm*5mm*300mm 強軸方向の理論値 一次モード:19252.5N 二次モード:56912.9N [#s61716e6]
,メッシュサイズ,要素数,一次モード[N],相対誤差[%],二次モード[N],相対誤差[%]
,0.8,124358,19470,1.13,56725,0.33
,0.9,96947,19540,1.50,56935,0.04
,1.0,92058,19530,1.44,56900,0.02
,2.0,13825,20625,7.13,60040,5.5


***15mm*5mm*300mm 強軸方向の理論値 一次モード:64977.1N 二次モード:192080.9N [#a82fd65c]
,メッシュサイズ,要素数,一次モード[N],相対誤差[%],二次モード[N],相対誤差[%]
,0.8,193357,64402.5,0.88,184485.0,3.95
,0.9,141860,64530.0,0.69,184815.0,3.78
,1.0,137149,64522.5,0.70,184807.0,3.79
,2.0,20585,66157.5,1.82,189442.0,1.37

***10mm*5mm*500mm 強軸方向の理論値 一次モード:6930.0NN 二次モード:20488.6N [#haaea865]
,メッシュサイズ,要素数,一次モード[N],相対誤差[%],二次モード[N],相対誤差[%]
,0.8,260685,7025,1.36,20650,0.79
,0.9,182376,7050,1.72,20725,1.15
,1.0,156700,7055,1.79,20745,1.25
,2.0,22573,7410,6.91,21770,6.25

***15mm*5mm*500mm 強軸方向の理論値 一次モード:23391.7N 二次モード:69149.1N [#ua052d61]
,メッシュサイズ,要素数,一次モード[N],相対誤差[%],二次モード[N],相対誤差[%]
,0.8,387258,23422.5,0.13,68422.5,1.05
,0.9,251524,23467.5,0.32,68565.0,0.84
,1.0,221101,22725.0,2.85,68595.0,0.80
,2.0,33142,24015,2.66,70155.0,1.45

,メッシュサイズ,要素数,一次モード[N],相対誤差[%],二次モード[N],相対誤差[%]
,0.8,,,,,
,0.9,,,,,
,1.0,,,,,
,2.0,,,,,

,メッシュサイズ,要素数,一次モード[N],相対誤差[%],二次モード[N],相対誤差[%]
,0.8,,,,,
,0.9,,,,,
,1.0,,,,,
,2.0,,,,,

表を見ると長さを長くすると梁に近づくため、理論値に近づくと言えそうだが、10mm*5mm断面と15mm*5mmを比較すると強軸の断面二次モーメントが大きい方が強軸方向の相対誤差が大きいと言えそう。
変形の様子はすべての構造でほぼ同様であったため、それぞれのモードについて10mm*5mm*500mmの変形の様子を下記に示す。


SAmechとSAgeomを間違えてやったらティッシュみたいになった。


http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/2024zakutu/4kyou1.500.png~

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/2024zakutu/4kyou2.500.png~

**片方固定他端ローラー 弱軸方向 [#a2d09c10]
***10mm*5mm*100mm [#jf63d07b]
理論値 一次モード:43221.98N 二次モード:128053.95N
,メッシュサイズ,要素数,一次モード[N],相対誤差[%],二次モード[N],相対誤差[%]
,0.6,100820,44640.95,3.28,127985,-0.054
,0.7,94071,44642.6,3.29,127975,-0.062
,0.8,60704,45158.15,4.48,129376,1.03
,0.9,21474,48677.6,12.6,139270.5,8.76
,1.0,19095,48591.65,12.4,139154.5,8.67
,2.0,2958,68000.5,57.3,192746,50.5

***10mm*5mm*300mm [#e49e8440]
理論値 一次モード:4802.39N 二次モード:14228.22N
,メッシュサイズ,要素数,一次モード[N],相対誤差[%],二次モード[N],相対誤差[%]
,0.6,401289,4956.08,3.20,14600.7,2.62
,0.7,137579,5181.15,7.89,15254.6,7.21
,0.8,124483,5183.8,7.94,15276.35,7.30
,0.9,97665,5259.35,9.52,15482.25,8.81
,1.0,90374,5254.35,9.41,15478.95,8.79
,2.0,13737,6349.8,32,18623.2,30.9

***10mm*5mm*500mm [#xa1056c4]
理論値 一次モード:1728.86N 二次モード:5122.16N
,メッシュサイズ,要素数,一次モード[N],相対誤差[%],二次モード[N],相対誤差[%]
,0.6,380721,1822.665,5.43,5380.55,5.04
,0.7,310163,1831.285,5.92,5405.75,5.54
,0.8,258632,1843.225,6.62,5441.95,6.24
,0.9,181166,1870.08,8.17,5519.1,7.75
,1.0,152023,1873.465,8.36,5529.95,7.96
,2.0,22405,2318.66,34.1,6826.65,33.3



----
CENTER:&size(25){総閲覧数&counter;};
CENTER:&size(25){今日の閲覧数&counter(today);};
CENTER:&size(25){昨日の閲覧数&counter(yesterday);};
----
CENTER:&aname(TOP){&size(30){''ページトップ''};};


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