青野の卒論日誌

2023年

卒論関連

振動解析(青山さん、七五三さんの続き) 不朽箇所1箇所ごとの感度は青山さんの手法で調べられるのだけど、それらの組み合わせの影響を何らかの統計的な手法(ベイズ推定とか)で評価できないかとか。

やること

10分割含水率変化

気乾比重は、乾燥させたときの比重のため少し(含水率10%程度)の水分を含んでいる。

全乾比重は、完全に乾燥させたときの比重(含水率0%)となっている。

スギでモデル作成していく.10✕20✕1000をモデルとし、気乾比重を最低値として考える。(5/9)(自分用メモ:解析していったデータは10bunnkatu dataに理論値とかの計算はsinndourironnkeisannにある)

種類密度g/cm3
バルサ0.10
リグナムバイタ1.25
キリ0.23
アカガシ0.97
スギ0.34
ヒノキ0.42
ケヤキ0.66
ミズナラ0.70

春休み課題(振動解析)

図心をZ軸に合わせて固定もZ軸上でのものとする(Boxの中心上にできるようにする)。そこから相対誤差がどのくらい変化するかの確認。(4/27)

追記:振動解析班のものは鉛直、水平のモードで分類してなかったので理論値計算を含めやり直す。(5/1)

相対誤差まで完成。ほとんど誤差のない結果となった。(5/2)

再解析データ

理論値

鉛直1次:36.036Hz 水平1次:72.073Hz 鉛直2次:144.146Hz 水平2次:288.292Hz

1次モード

メッシュの長さ鉛直1次モード相対誤差水平1次モード相対誤差
0.836.03070.015%72.02190.067%
0.936.03070.015%72.02190.067%
1.036.03070.015%72.02190.067%
2.036.03100.014%72.02240.066%
3.036.03080.014%72.02240.066%
4.036.03120.013%72.02260.066%
5.036.03130.013%72.02270.066%

2次モード

メッシュの長さ鉛直2次モード相対誤差水平2次モード相対誤差
0.8144.0530.065%287.4770.283%
0.9144.0530.065%287.4770.283%
1.0144.0530.065%287.4770.283%
2.0144.0580.061%287.4840.280%
3.0144.0550.063%287.4850.280%
4.0144.0600.060%287.4880.279%
5.0144.0620.058%287.4900.278%

振動解析班のデータ1次モード

メッシュの長さ要素数一次モードの固有振動数解析値相対誤差
0.825607736.01940.0294%
0.912808636.01940.0294%
1.08946236.01950.0291%
2.02107436.01960.0288%
3.02445736.01960.0288%
4.0547336.01990.0280%
5.0572536.020.0277%

振動解析班のデータ二次モード

メッシュの長さ要素数二次モードの固有振動数解析値相対誤差
0.8256077108.0674.03%
0.9128086108.1833.93%
1.089462108.3583.78%
2.021074109.1743.05%
3.024457109.0613.15%
4.05473109.9062.40%
5.05725110.0242.30%

参考文献

構造振動学

既存鋼鈑桁橋のモデルパラメータ事後分布を用いた信頼性評価

ベイズ推定による既存構造物数値モデルの不確定性定量化とキャリブレーション

読んでみたいもの

(Model Verification & Validation(V&V))

Alvin, K. F., Oberkampf, W. L, Dieger, K., Rutherford, B.:Uncertainty quantification in computational structural dynamics: a new paradigm for model validation, Society for Experimental Mechanics, Inc, 16 th International Modal Analysis Conference.. Vol. 2. 1998. (計算構造力学における不確かさの定量化:モデル検証のための新しいパラダイム、実験力学研究会、16th国際モード解析会議.   Vol.2. 1998.)

2022年

創造工房実習の内容

10月21日課題

タッチタイプ練習

11月4日

・課題

https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/kadai.png

秋田市の月別平均気温(1991~2020)

・練習

https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/bbb.png

11月11日 片持ち梁

メッシュの長さ要素数変位相対誤差計算者
0.5595046.56-1.5千代岡
0.6455126.48774-2.69高井
0.7390756.54133-2.0関合
0.8133976.43695-3.5岡田
0.999036.36315-4.6松田
1.262566.3043375-5.4青野
1.357676.29784-5.6山口
1.451996.29990-5.55山本
1.539356.24807-6.3進藤
1.634006.20446-6.98河合
1.829526.17161-7.5山口
216325.64585-15.3進藤
36825.47288-17.9山本
42643.6161-45.8関合
51913.86-42千代岡
61902.5077325-62.4高井
7751.41225-78.8青野
8561.2887175-80.7岡田
9491.28799-80.9松田
10441.226075-81.6河合

片持ち梁のグラフ:縦軸 変位 横軸 要素数

https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/katamotibari2.png

11/18(金) 単純梁

メッシュ長さ要素数変位(mm)相対誤差(%)計算者
0.56041670.4289822.94千代岡
0.63615840.4212331.09高井
0.71452340.42251.4関合
0.81409870.4226273851.4岡田
0.9918570.4203516060.88松田
1.2245200.404744325-2.87青野
1.3231320.4045-2.93山口
1.4175300.3986-4.34山本
1.5154330.3963177567574.9進藤
1.6159000.399049-4.24河合
1.8116770.404457-0.03山口
2104600.394818715517-5.3進藤
323440.32447-22.13山本
414530.3329-20.1関合
54310.136240-67.3千代岡
63600.2130486-48.9高井
71960.1019892-75.5青野
81040.1158624-72.2岡田
9810.1247076-70.1松田
10780.07733-81.4河合

片持ち梁のグラフ:縦軸 変位 横軸 要素数

https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/tanjunbari2.png

11月25日

等方性1次と等方性2次の比較

二次要素の等方性のデータ

メッシュ長さ要素数変位(mm)相対誤差(%)計算者
0.5604167--千代岡
0.62032090.423830.98高井
0.71452340.433013.22関合
0.81409870.430063.2岡田
0.9919740.4299133.18松田
1.2248000.4297773.14青野
1.3231320.429893.16山口
1.4176170.4297453.13山本
1.5154330.4298443.2進藤
1.6159000.4297543.13河合
1.8116770.429623.1山口
2104600.4296053.1進藤
324860.4292173.0山本
414530.42933.02関合
54310.4278852.69千代岡
63600.42822.78高井
71960.426062.25青野
81040.426312.3岡田
9810.425132.03松田
10780.4244661.8河合

https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/kadai122.png

等方性1次と異方性1次の比較

異方性1次のデータ

メッシュ長さ要素数変位(mm)相対誤差(%)計算者
0.56041670.509193.56千代岡
0.62032090.504722.6高井
0.71452340.50362.42関合
0.81409870.502832.3岡田
0.9919740.500531.8松田
1.2248000.48739-0.9青野
1.3231320.48841-0.67山口
1.4176170.48403-1.56山本
1.5154330.48202-2.0進藤
1.6159000.48329-1.7河合
1.8116770.47855-2.67山口
2104600.47906-2.6進藤
324360.42787-12.98山本
414530.42772-13.02関合
54310.27364-44.3千代岡
63600.33927-31.0高井
71960.21363-58.5青野
81040.22574-54.1岡田
9810.22750-53.7松田
10780.20327-58.7河合

https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/kadai1222.png

12/2 木材と鋼材の単純梁 サンドウィッチ

異方性の単純梁(1次)と2次要素のサンドウィッチ梁(鋼材:等方性 木材:異方性 木材を鋼材で挟む)の比較

メッシュ要素数変位相対誤差測定者
0.71552660.0861-13.0関合
0.81384530.08349-15.7岡田
0.9827660.083312-15.9松田
1.2322790.08357-15.6青野
1.3283430.083667-15.5山口
1.4236670.08368-15.5山本
1.5199580.083516-15.6進藤
1.6194510.086037-13.1河合
1.8109330.084021-15.1山口
2107640.083323-15.8進藤
336180.083467-15.7山本
416230.0852-13.9関合
510070.083104-16.1千代岡
68420.0821-17.1高井
75540.08075-18.4青野
82890.079715-19.5進藤
92610.078427-20.8松田
102320.082495-16.7河合

https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/sandwich.png

<メモ>

UNIXコマンド

gnuplot

inkscape

LaTex

LibreOfficeImpress

Salome-Meca

自分メモ

全角/半角の切り替えは、「無変換」(左手親指)でもできる

コピペは、左クリックで領域選択したら、スクロールボタンで貼り付け

単語選択は、ダブルクリックでも可(トリプルクリックだと行選択かな)


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