#author("2021-10-11T16:51:07+09:00","default:kouzouken","kouzouken")
#author("2022-01-13T19:00:14+09:00","default:kouzouken","kouzouken")
#contents


*卒論日誌 [#te07cc65]
卒論テーマ「めおと橋の解析」

新めおと橋の初期調査を行い,実験値と解析値を比較する.
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/arisa/meoto600.JPG

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/arisa/meoto600.JPG

参考文献:[[構造用木材の強度試験マニュアル:https://www.howtec.or.jp/files/libs/1828/201712121507021978.pdf]]

,日付,作業内容
,4/19~,めおと橋と劣化診断法について調べる
,4/26~,単純梁の振動解析(ソリッド要素)
,5/10~,振動解析とめおと橋調査計画書作成
,5/17~,前回の続き
,5/24~,振動解析とめおと橋調査計画書作成
,5/31~,旧めおと橋の部材試験
,6/14~,旧めおと橋部材の試験データまとめ
,7/5~,重回帰分析
,8/23,めおと橋現地調査1回目
,9/10,めおと橋現地調査2回目
,10/4,横構の振動解析
,10/4~,横構の振動解析
,10/11~,横構の拘束条件パターン1−3
,10/18~,床桁の振動解析
,10/25~,めおと橋全体解析
,11/8,中央載荷・等分布荷重
,12/6~,アーチの応力照査
,12/13~,曲げ応力照査

**振動使用性 [#g3297371]

歩行者の歩調2[Hz]前後

歩道橋の固有振動数が歩調に一致すると、振動数が増幅され不快感を与える.
-一次・二次の固有振動数は、1.5~2.3[Hz]以下を避ける.
-活荷重たわみの許容値は支間の1/600

参考文献:[[設計時における歩道橋の振動使用性照査法:https://www.jstage.jst.go.jp/article/jscej1984/1993/471/1993_471_77/_pdf]]

**中央載荷 (断面300✕1650)[#x06d48a4]
床版(死荷重+雪荷重):70.564+46.92=117.484[kN]

アーチリブ(死荷重+雪荷重):41.947+9.177=51.124[kN]

活荷重:3.5[kN/m^2]

たわみ→2.51[mm]

**12/13~ 曲げ応力照査 [#lca63e7e]
設計書曲げ応力度

$F_{b}$=6.21[N/mm^2]

***設計時ヤング率(7GPa) [#s6c331ea]
,箇所,圧縮応力の最大,評価
,アーチリブ9①,1.8134,○
,アーチリブ9②,1.43026,○
,アーチリブ7①,2.16866,○
,アーチリブ7②,2.00551,○
,アーチリブ5①,1.75546,○
,アーチリブ5①,1.91126,○

***測定ヤング率 [#b1c69a67]
,箇所,圧縮応力の最大,評価
,アーチリブ9①,1.9914,○
,アーチリブ9②,1.5714,○
,アーチリブ7①,2.3931,○
,アーチリブ7②,2.2060,○
,アーチリブ5①,1.9396,○
,アーチリブ5②,2.1075,○

**12/6~ アーチの応力照査 [#z09e91a3]
設計書許容応力度

$F_{cy'}$=1.4798

床版(死荷重+雪荷重):70.564+46.92=117.484[kN]

アーチリブ(死荷重+雪荷重):41.947+9.177=51.124[kN]

,箇所,圧縮の最大値,評価
,アーチリブ9①,1.8134,✕
,アーチリブ9②,1.43026,○
,アーチリブ7①,2.16866,✕
,アーチリブ7②,2.005509,✕
,アーチリブ5①,1.75546,✕
,アーチリブ5①,1.91126,✕

**11/29~ 設計時ヤング率 [#m8418286]
設計時ヤング率→7.15[GPa]

***全体解析 [#y89cbac7]
,振動モード,設計時固有振動数[Hz],現固有振動数[Hz]
,水平方向,3.68,3.60
,鉛直方向,13.17,13.40
,鉛直逆対称方向,16.17,16.10

***中央載荷・等分布荷重 [#ked32dbc]
,荷重箇所,荷重[kN/m^2],たわみ[mm](設計時),たわみ[mm](9月測定),許容値(38.3mm)
,中央載荷,活荷重 3.5,0.16,0.150,合格
,中央載荷,活荷重+死荷重 6.75,0.30,0.29,合格
,等分布荷重,雪荷重 14,4.40,4.52,合格
,等分布荷重,活荷重 3.5,1.1,1.13,合格
,等分布荷重,雪荷重+死荷重 17.25,5.43,5.58,合格
,等分布荷重,活荷重+死荷重 6.75,2.12,2.18,合格


**11/22~ 新橋のヤング率を落とす [#b0bd62d5]

-旧部材のヤング率を入れてみて、どのくらい振動数が下がるのかみてみる.
,振動モード,固有振動数[Hz],現固有振動数[Hz]
,水平方向,3.60(±0),3.60
,鉛直方向,13.40(+1.02),12.38
,鉛直逆対称方向,16.10(+0.46),15.64


**11/15~ 中央載荷・等分布荷重 [#bc64c163]

,荷重箇所,荷重[kN/m^2],たわみ[mm],許容値(38.3mm)
,中央載荷,活荷重 3.5,0.150,合格
,中央載荷,活荷重+死荷重 6.75,0.29,合格
,等分布荷重,雪荷重 14,4.52,合格
,等分布荷重,活荷重 3.5,1.13,合格
,等分布荷重,雪荷重+死荷重 17.25,5.58,合格
,等分布荷重,活荷重+死荷重 6.75,2.18,合格

***中央載荷 (断面300✕1650)[#x06d48a4]
活荷重→3.5[kN/m^2]

たわみ→0.304[mm]

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/fujiwara/meotosyutyu.png

***等分布荷重 (断面1650✕2300)[#na922ab8]
雪荷重→14[kN/m^2]
 雪荷重を考慮する必要がある場合2通りがある.
 ①十分圧縮された雪の上を自由に車両が通行する場合→1[kN/m^2]
 ②積雪が特に多くて自動車交通が不能となり、雪だけが荷重としてかかる場合→3.5[kN/m^3]

→今回は、②を用いる.積雪4mを想定.

たわみ→5.42[mm]

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/fujiwara/meototoubun.png

活荷重→3.5[kN/m^2]

たわみ→1.56[mm]


**10/25~ めおと橋全体解析 [#o6bf4c82]

***初期値 固有振動数(等方性)[#zb7eaae3]
,振動モード,固有振動数[Hz]
,水平方向,3.60
,鉛直方向,13.40
,鉛直逆対称方向,16.10             

水平方向(3.60Hz)

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/fujiwara/meotozentai3.gif

鉛直方向(13.40Hz)

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/fujiwara/meotozentai2.gif

逆対称鉛直(16.10Hz)

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/fujiwara/meotozentai4.gif

**10/18~ 床桁の振動解析 [#nbef728d]

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/fujiwara/yokogeta.png

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/fujiwara/yukageta_genti.png

拘束条件→ラブスクリューφ16、両端固定

***解析結果 [#ka0fbf69]
,振動次数,振動モード,理論値①[Hz](2990mm両端固定),理論値②[Hz](1300mm両端固定),理論値③[Hz](1300mm単純支持),解析値[Hz]
,1次,水平一次?,100.84,1686.92,744.17,216.129
,2次,水平二次?,277.96,4649.89,2976.67,414.361
,3次,鉛直一次,166.39,2783.41,1227.88,573.915
,4次,水平三次,,,,606.186
,5次,鉛直二次,458.64,7672.31,2976.67,658.868

***現地調査のデータ [#te202a0c]
,固有振動数[Hz],両端固定一次FFTヤング率[GPa](理論値に対する相対誤差),両端固定二次FFTヤング率[GPa](理論値に対する相対誤差)
,86.1,5.47(-27.10),0.72(-90.41)
,161.5,19.24(156.49),2.53(-66.26)
,172.3,21.90(191.94),2.88(-61.60)
,215.3,34.19(355.84),4.50(-40.04)
,258.4,49.25(556.61),6.48(-13.62)
,269.2,53.45(612.65),7.03(-6.25)

FAKOPPヤング率→8.3[GPa]


***解析値と現地調査データの比較 [#s1a0436e]
,FFT1次モード[Hz],FFT2次モード[Hz],salome解析値1次モード[Hz](現地調査に対する相対誤差),salome解析値2次モード[Hz](現地調査に対する相対誤差)
,86.1,269.2,216.129(151.02),414.361(53.92)

→解析結果が現地データよりもだいぶ固い。

→固定条件を考慮

→振動モードがわかりにくい

**10/11~ 横構 拘束条件パターン3 [#mcb1513e]

拘束条件→ボルトで固定+木材とプレートの間に1mmの隙間あり

木材:
-ヤング率→7.5[GPa]
-ポアソン比→0.4
-密度→390[kg/m^3]

プレート(ss400):
-ヤング率→206[GPa]
-ポアソン比→0.3
-密度→7850[kg/m^3]

ボルト・ドリフトピン(ss400):
-ヤング率→206[GPa]
-ポアソン比→0.3
-密度→7850[kg/m^3]


***解析結果 [#ac90576d]
,振動次数,振動モード,λ,理論値[Hz],解析パターン3[Hz](理論値に対する相対誤差%)
,1次,水平一次,4.73,128.94,214.494(+66.35)
,2次,鉛直一次,4.73,309.44,412.890(+33.43)
,3次,ねじれ,,,530.642
,4次,水平二次,7.853,355.40,556.051(+56.46)
,5次,鉛直二次,7.853,852.96,947.279(+11.06)

***解析値と現地調査データの比較 [#g506d85f]
,部材名,FFT1次モード[Hz],FFT2次モード[Hz],salome解析値1次モード[Hz](現地調査に対する相対誤差),salome解析値2次モード[Hz](現地調査に対する相対誤差)
,横構1,236.9,689.1,214.494(-9.46),556.051(-19.31)
,横構2,269.2,646.0,214.494(-20.32),556.051(-13.92)
,横構3,269.2,646.0,214.494(-20.32),556.051(-13.92)
,横構4,239.6,699.8,214.494(-10.48),556.051(-20.54)

**10/11~ 横構 拘束条件パターン2 [#mf4d1fdd]

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/fujiwara/yokokou3.png

拘束条件→両端固定(プレートと木材の間に1mmの隙間あり)

木材:
-ヤング率→7.5[GPa]
-ポアソン比→0.4
-密度→390[kg/m^3]

プレート(ss400):
-ヤング率→206[GPa]
-ポアソン比→0.3
-密度→7850[kg/m^3]

***解析結果 [#ac90576d]
,振動次数,振動モード,解析パターン2[Hz],λ,理論値[Hz],相対誤差[%]
,1次,水平一次,123.051,4.73,128.94,-4.57
,2次,鉛直一次,265.691,4.73,309.44,-14.14
,3次,水平二次,299.013,7.853,355.40,-15.87
,4次,ねじれ,351.635,,,
,5次,水平三次,534.704,10.996,696.82,-23.27
,6次,鉛直二次,638.216,7.853,852.96,-25.18

***解析値と現地調査データの比較 [#g506d85f]
,部材名,FFT1次モード[Hz],FFT2次モード[Hz],salome解析値1次モード[Hz](現地調査に対する相対誤差),salome解析値2次モード[Hz](現地調査に対する相対誤差)
,横構1,236.9,689.1,123.051(-48.01),265.691(-61.44)
,横構2,269.2,646.0,123.051(-54.39),265.691(-58.87)
,横構3,269.2,646.0,123.051(-54.39),265.691(-58.87)
,横構4,239.6,699.8,123.051(-48.64),265.691(-62.03)

**10/11~ 横構 拘束条件パターン1 [#fd60aacc]

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/fujiwara/yokokou2.png

拘束条件→ボルトで固定(ボルトとプレートを固定)

木材:
-ヤング率→7.5[GPa]
-ポアソン比→0.4
-密度→390[kg/m^3]

プレート(ss400):
-ヤング率→206[GPa]
-ポアソン比→0.3
-密度→7850[kg/m^3]

ボルト・ドリフトピン(ss400):
-ヤング率→206[GPa]
-ポアソン比→0.3
-密度→7850[kg/m^3]


***解析結果 [#ac90576d]
,振動次数,振動モード,λ,理論値[Hz],前回のSalome解析値[Hz],解析パターン1[Hz](理論値に対する相対誤差%)
,1次,水平一次,4.73,128.94,297.199,223.839(+73.60)
,2次,鉛直一次,4.73,309.44,481.95,449.816(+45.36)
,3次,ねじれ,,,830.86,558.245
,4次,水平二次,7.853,355.40,839.319,587.935(+65.43)
,5次,鉛直二次,7.853,852.96,1097.91,1022.49(+19.86)


***現地調査のデータ [#c4d3f888]
,部材名,FFT(1次モード)[Hz](ヤング率[GPa]),FFT(2次モード)[Hz](ヤング率[GPa]),FAKOPPヤング率[GPa]
,横構1,236.9(6.19),689.1(6.89),5.70
,横構2,269.2(7.99),646.0(6.05),7.95
,横構3,269.2(7.99),646.0(6.05),6.84
,横構4,239.6(6.33),699.8(7.10),7.89

***解析値と現地調査データの比較 [#g506d85f]
,部材名,FFT1次モード[Hz],FFT2次モード[Hz],salome解析値1次モード[Hz](現地調査に対する相対誤差),salome解析値2次モード[Hz](現地調査に対する相対誤差)
,横構1,236.9,689.1,223.839(-5.51),587.935(-14.68)
,横構2,269.2,646.0,223.839(-16.85),587.935(-8.99)
,横構3,269.2,646.0,223.839(-16.85),587.935(-8.99)
,横構4,239.6,699.8,223.839(-6.58),587.935(-15.99)


**10/6 横構の振動解析 [#dd682abf]

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/fujiwara/yokokou.png

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/fujiwara/yokokou_genti.png

拘束条件→両端固定(プレート部分)
拘束条件→両端固定(プレートと木材を完全固定)

木材:
-ヤング率→7.5[GPa]
-ポアソン比→0.4
-密度→390[kg/m^3]

プレート(ss400):
-ヤング率→206[GPa]
-ポアソン比→0.3
-密度→7850[kg/m^3]

***解析結果 [#ac90576d]
,振動次数,振動モード,Salome実験値[Hz],λ,理論値[Hz],相対誤差[%]
,1次,水平一次,297.199,4.73,260.85,13.93
,2次,鉛直一次,481.95,4.73,626.04,-23.02
,3次,ねじれ一次,830.86,,,
,4次,ねじれ二次,839.319,,,
,4次,水平二次?,839.319,7.853,1157.61,37.92
,5次,鉛直二次,1097.91,7.853,1725.64,-36.38
,6次,伸び,1336.56,,,
,7次,水平二次,1517.34,7.853,1157.61,31.08


曲げ理論値:f=(1/(2π))*(λ/l)^2*√(EI/ρA) 
ねじれ理論値:f=λ/2πL*√(GJ/ρIp)


**9/10 めおと橋現地調査結果(2回目) [#y6c1e468]
FFTは対傾構のみ,FKOPPは前回の結果の相対誤差28[%]以上の箇所を再調査した.

***FFT測定結果 [#se43745e]
これは,対傾構1(入り口から1本目)の測定結果である.

ヤング率7.5[GPa]のとき,固有振動数の理論値が308[Hz]であるが,波形が大きく凸になっていなかったため,読み取ることができなかった.対傾構2,3,4も同様.

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/fujiwara/meoto_tai10910.png

***FAKOPP測定結果 [#w52c3161]

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/fujiwara/meoto_FAKOPP0910.png

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/fujiwara/mroto_ue.png

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/fujiwara/meoto_sita.png

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/fujiwara/meoto_enter1.png

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/fujiwara/meoto_exit2.png

今回は,アーチ継ぎ手は危険な場所のため調査を中止した.

アーチリブも前回と測定する箇所を変更した.

2箇所以外は,相対誤差28[%]におさまった.

平均値7.40[GPa]、標準偏差1.23、全体変動係数17.22

,部材名,ヤング率,変動率
,垂直材,
,高欄,
,アーチリブ,
,床桁,
,縦桁,
,横構,

***FFT測定結果 [#r10cc0ff]

**8/23 めおと橋現地調査結果 [#i10cfa63]
各部材FAKOPPで伝播速度,FFTで固有振動数を測定し,それらヤング率を測定する.

***FAKOPP測定結果 [#c1291879]

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/fujiwara/meoto_FAKOPP0823.png

***FFT測定結果 [#u0f442f7]





**重回帰分析結果 [#nef27de6]


目的関数を治具なし(45°),説明関数を治具ありの測定値・そのときの角度・当木長さとして,重回帰分析を行った.


側面のときは,説明関数にセンサー間距離も考慮した.


○内側の結果

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

X1:治具ありの測定値[GPa],X2:角度[°],X3:当木長さ[mm]

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



○側面の結果

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

X1:治具ありの測定値[GPa],X2:角度[°],X3:当木長さ[mm],X4:センサー間距離[mm]

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



どちらの場合も補正値(R2)が0に近いことから,今回の分析結果の精度は低い.


結果の有意Fの値をみると,0.05未満であることから推定された回帰式は統計的には意味があると言える.


またP-値をみると,測定値(X1)は0.05未満であることから,目的関数に対して"関係がある”と言えるが,他の説明関数に関しては”関係がない”.


**旧めおと橋FAKOPPデータまとめ [#d1203536]

○内側(外)・150mmのとき
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/fujiwara/uti150.png

○内側(外)・120mmのとき
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/fujiwara/uti120.png

○内側(外)・100mmのとき
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/fujiwara/uti100.png

○内側(外)・80mmのとき
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/fujiwara/uti80.png

○内側(外)・50mm安定のとき
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/fujiwara/uti501.png

○内側(外)・50mm不安定のとき
-http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/fujiwara/uti502.png


**5/24~ 単純梁振動解析 [#t41b1ddf]
前回までは,固定の位置を部材の下面にしていたことが原因で理論値と大きく誤差が生まれてしまったと考える。中立軸上を固定して再度解析していく.

木材で解析を行うと材料が柔らかく,固定しても材料が変形してしまうため,今回は鋼材で解析してみる.また,前回は材料の長さが短かったため,今回は長さを大きくしてみる.
-モデル420mm✕10mm✕20mm,支間長400[mm]の単純梁
-鋼材ss400,ヤング率206[GPa],ポアソン比0.3,密度7.693e-0.5[N/mm^3]

,振動次数,振動モード,Salome実験値[Hz],λ,理論値[Hz],相対誤差[%]
,1次,水平一次,-1.4847,π,1.4665,1.24
,2次,水平二次,-5.91879,2π,5.8662,0.90
,3次,,-6.11798,,,
,4次,水平三次,-13.2445,3π,13.1989,0.35
,5次,,-16.6127,,,
→理論値と近くなった.

→振動次数三次,五次のときは,水平と鉛直の振動が混じっていた.これはモデルの断面を10mm✕20mmに設定したことが原因であると考える.

木材のときはどうなるかもやってみる.
-モデル420mm✕10mm✕20mm,支間長400[mm]の単純梁
-ヤング率6[GPa],ポアソン比0.4,密度3.8e-0.7[N/mm^3]
,振動次数,振動モード,Salome実験値[Hz],λ,理論値[Hz],相対誤差[%]
,1次,水平一次,-3.61575,π,3.561,1.54
,2次,水平二次,-14.4135,2π,14.245,1.18
,3次,,-14.949,,,
,4次,水平三次,-32.2495,3π,32.051,0.62
,5次,,-40.5143,,,
→木材でも理論値と近くなった.


**5/17~ 前回の続き [#n730bbe2]

前回
→モデル120mm✕10mm✕20mmに支間長100mmとしているので,Salomeが自由端で解析してのではないかと考える。
 
自由端の理論値と比較してみる。

,振動次数,振動モード,Salome実験値[Hz],λ,理論値[Hz],相対誤差[%]
,1次,水平一次,-86.7021,4.730,89.696,-3.34
,2次,鉛直一次,-154.058,4.730,179.392,14.12
,3次,水平二次,-190.500,7.853,247.242,-22.95
,4次,ねじれ,-275.261,,,
,5次,鉛直二次,-378.780,7.853,494.484,23.40

→両端固定のλで計算したときよりも相対誤差が小さくなった。



モデル100mm✕10mm✕20mmの支間長100mmで解析してみる。
今回はメッシュサイズ0.8とした。
-モデル100mm✕10mm✕20mmの支間長100mmの単純梁
-ヤング率6[GPa],ポアソン比0.4,密度3.8e-0.7[N/mm^3]

,振動次数,振動モード,Salome実験値[Hz],λ,理論値[Hz],相対誤差[%]
,1次,水平一次,-81.6605,π,56.979,43.32
,2次,鉛直一次,-132.448,π,113.958,16.23
,3次,水平二次,-187.168,2π,227.915,-17.88
,4次,ねじれ,-230.831,,,
,5次,鉛直二次,-335.125,2π,455.830,-26.48

モデル100mm✕10mm✕20mmを上下固定して解析してみる。
-モデル100mm✕10mm✕20mmの支間長100mmの単純梁
-ヤング率6[GPa],ポアソン比0.4,密度3.8e-0.7[N/mm^3]

,振動次数,振動モード,Salome実験値[Hz],λ,理論値[Hz],相対誤差[%]
,1次,水平一次,-114.68,π,56.979,101.27
,2次,鉛直一次,-164.278,π,113.958,44.16
,3次,水平二次,-287.197,2π,227.915,26.01
,4次,ねじれ,-301.221,,,
,5次,鉛直二次,-402.265,2π,455.830,11.75

→振動次数が大きくなるにつれて誤差が小さくなった。


**5/10~ 前回の続き [#v996d986]
今回はモデル120mm✕10mm✕20mmの長方形断面,支間長100mmの単純梁の振動解析を行なった。
-ヤング率6[GPa],ポアソン比0.4,密度3.8e-0.7[N/mm^3]
-固有振動数の理論式 f=(1/(2π))*(λ/l)^2*√(EI/ρA) 

,振動次数,振動モード,Salome実験値,λ,理論値,相対誤差[%]
,1次,水平一次,-86.7021,π,56.979,52.17
,2次,鉛直一次,-154.058,π,113.958,35.19
,3次,水平二次,-190.500,2π,227.915,16.42
,4次,ねじれ,-275.261,,,
,5次,鉛直二次,-378.780,2π,455.830,16.90



先輩の振動解析を参考にして,片持ち梁の振動解析ができるのかやってみる。

-片持ち梁100mm✕10mm✕20mmの振動解析
-ヤング率6[GPa],ポアソン比0.4,密度3.8e-0.7[N/mm^3]
-固有振動数の理論式 f=(1/(2π))*(λ/l)^2*√(EI/ρA) 

,振動次数,振動モード,Salome実験値,λ,理論値,相対誤差[%]
,1次,水平一次,-20.8188,1.875,20.296,2.58
,2次,鉛直一次,-39.9531,1.875,40.593,-1.58
,3次,水平二次,-124.400,4.694,127.204,-2.20
,4次,ねじれ,-145.725,,,
,5次,鉛直二次,-213.841,4.694,254.407,-15.95

青山さんの結果と同じになったので,Aster Studeyの設定はあっていると思う。


**4/26~ 単純梁の振動解析 [#yd2f6759]
モデル120mm✕10mm✕10mmの正方形断面,支間長100mmの単純梁の振動解析を行なった。
-ヤング率6[GPa],ポアソン比0.4,密度3.8e-0.7[N/mm^3]

-固有振動数の理論式 f=(1/(2π))*(λ/l)^2*√(EI/ρA) を用いて,比較を行う。

,振動次数,振動モード,Salome実験値,λ,理論値,相対誤差[%]
,1次,水平一次,-85.8333,π,56.979,50.94
,2次,鉛直一次?,-89.575,π,56.979,57.21
,3次,水平二次,-188.843,2π,227.915,-17.15
,4次,鉛直二次?,-243.13,2π,227.915,-6.68
,5次,ねじれ?,-286.843,,,

→λは[[日本機械学会構造振動学:https://www.jsme.or.jp/sed/guide/dynamics5.pdf]]を参考にした。

理論値と大きく異なった。アニメーションで振動モードを確認したがよくわからなかった。
次回は,わかりやすいように長方形断面で解析を行ってみる。

*作業日誌 [#vfe52921]
,日付,時間帯,作業時間,内容,立会
,10/16,14:30-16:00,1.5h,顔合わせ、pcの使い方,後藤さん、4年生
,10/30,14:30-16:00,1.5h,vi gnuplotの使い方,後藤さん
,11/6,14:30-17:00,2.5h,SALOMEの使い方,後藤さん、4年生、及川さん
,11/6,17:00-18:30,1.5h,課題,
,11/12,17:00-17:30,0.5h,課題,
,11/13,14:30-17:00,2.5h,単純梁の線形解析,4年生、及川さん
,11/18,14:30-17:00,2.5h,課題,
,11/19,15:00-14:00,1.0h,課題,
,11/20,12:00-13:00,1.0h,課題,
,11/20,14:30-16:00,1.5h,単純梁の線形解析②,4年生、及川さん、後藤さん
,12/3,17:00-19:00,2.0h,課題,
,12/4,13:20-14:20,1.0h,課題,
,12/4,14:30-17:00,2.5h,2材料(鋼材、木材),4年生、及川さん、後藤さん
,12/10,15:30-19:00,3.5h,課題,
,12/11,14:30-17:30,3.0h,texの使い方,後藤さん
,12/17,17:00-18:30,1.5h,課題,
,12/18,14:30-16:30,2.0h,texの使い方,後藤さん
,1/22,14:30-16:30,2.0h,最終課題,後藤さん
,1/28,16:30-18:30,2.0h,最終課題,
,1/29,14:30-19:00,4.5h,最終課題,後藤さん
,2/5,14:00-17:00,3.0h,最終課題,後藤さん
,2/10,16:30-18:00,1.5h,最終課題,
,2/12,13:30-16:30,3.0h,発表,後藤さん、4年生
,3/9,13:00-17:00,4.0h,進捗状況,後藤さん
,4/5,16:00-16:30,0.5h,顔合わせ,後藤さん
,4/19,14:30-16:00,1.5h,課題発表,後藤さん
,4/26,14:30-17:00,3.5h,進捗報告,後藤さん

**11/6課題分(salome meca 片持ち梁) [#e8de36b0]

ヤング率:6000N/mm^2 ポアソン比:0.4 荷重:100Nとした
,メッシュ長さ,要素数,先端変位,相対誤差(%),計算者
,0.5,394121,6.57938,-1.4,君島
,0.7,130916,6.4781,-2.8,君島
,0.8,72101,6.43695,-3.5,高橋
,0.9,71718,6.43136,-3.6,高橋
,1.0,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,森島

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

**11/13課題分(salome meca 単純梁) [#jcf35d32]

ヤング率:7500N/mm^2 ポアソン比:0.4 荷重100Nとした
,メッシュ長さ,要素数,先端変位,相対誤差(%),計算者
,0.5,391031,0.341181238095,2.35,君島
,0.7,215780,0.337975375,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.0,12228,0.3005115,-6.9,藤原
,4.0,5077,0.28405475,-13.9,藤原
,8.0,1795,0.2312003333,-30.6,森島
,10,752,0.1612725,-51.6,森島

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

**11/20課題分(単純梁 線形解析) [#zac82f19]

支間長100mmのとき
,メッシュ長さ,要素数,先端変位,相対誤差(%),計算者
,0.5,391031,0.440015,11.9,君島
,0.7,215780,0.435823,10.8,君島
,0.8,159468,0.4330131,10.1,高橋
,0.9,90071,0.42766,8.72,高橋
,1.0,61315,0.423881,7.77,田村
,1.2,58111,0.423005,7.54,田村
,1.4,47409,0.420309,6.86,田村
,1.5,42068,0.418470375,6.39,根本
,1.8,24627,0.410464142857,4.36,根本
,2.0,12228,0.396314,0.84,藤原
,4.0,5077,0.378695,-3.6,藤原
,8.0,1795,0.342299,-12.7,森島
,10,752,0.298709,-24.0,森島

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/arisa/tanjyun1120.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.11413,59.3,高橋
,1.0,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.0,13986,0.069684,-2.8,藤原
,4.0,3009,0.048789575,-32.0,藤原
,8.0,967,0.0764429,7.2,森島
,10,558,0.0768385,6.7,森島

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

**12月4日課題分(2材料) [#pceb4b26]

,メッシュ長さ,要素数,先端変位,相対誤差(%),計算者
,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.0725491125,41.3,根本
,1.8,34858,0.068374885714,33.2,根本
,2.0,20313,0.063280133,23.3,藤原
,3.0,18229,0.0489236,-4.68,君島
,4.0,8067,0.050046,-2.51,藤原
,5.0,4846,0.036772667,-28.3,君島
,8.0,3814,0.02708776667,-47.2,森島
,10,1716,0.0217906,-57.5,森島

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

**創造工房 最終課題 [#af9cbe1d]

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

**春休み課題 [#cf9267cb]
-材料非線形について
-鋼材1mm✕1mm✕100mmの片持ち梁の端部に荷重をかけた
-材料諸元(鋼材ss400)
ヤング率:206GPa
ポアソン比:0.3

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/arisa/hisenkei1.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/arisa/hisenkei2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2020/arisa/hisenkei3.png


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