青野の卒論日誌の履歴(No.26)

2023年†

卒論関連†

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

やること†

簡単な橋を作成し（10✕20✕500とか？）、ミルフィーユ状に切っていくことで含水率の変化による固有振動の変化を確認する。最初はサンドイッチみたいなのからスタートして徐々に層数を増やして密度（ρ）を層ごとに変化させたときの固有振動数の変化を見ていく

種類	密度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）

再解析データ

1次モード

メッシュの長さ	鉛直1次モード	水平1次モード
0.8	36.0307	108.947
0.9	36.0307	109.251
1.0	36.0307	109.456
2.0	36.031	109.873
3.0	36.0308	109.906
4.0	36.0312	110.273
5.0	36.0313	110.403

2次モード

メッシュの長さ	鉛直2次モード	水平2次モード
0.8	144.053	352.678
0.9	144.053	353.569
1.0	144.053	354.161
2.0	144.058	355.405
3.0	144.055	355.517
4.0	144.06	356.621
5.0	144.062	357.02

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

メッシュの長さ	要素数	一次モードの固有振動数解析値	相対誤差
0.8	256077	36.0194	0.0294%
0.9	128086	36.0194	0.0294%
1.0	89462	36.0195	0.0291%
2.0	21074	36.0196	0.0288%
3.0	24457	36.0196	0.0288%
4.0	5473	36.0199	0.0280%
5.0	5725	36.02	0.0277%

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

メッシュの長さ	要素数	二次モードの固有振動数解析値	相対誤差
0.8	256077	108.067	4.03%
0.9	128086	108.183	3.93%
1.0	89462	108.358	3.78%
2.0	21074	109.174	3.05%
3.0	24457	109.061	3.15%
4.0	5473	109.906	2.40%
5.0	5725	110.024	2.30%

春休み課題の振動班のものをやってみる。
参考　https://www.youtube.com/watch?v=5SYHxrAr4f8
春休み振動解析班
振動解析

↑

参考文献†

構造振動学

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

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

↑

読んでみたいもの†

（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日†

・課題

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

・練習

↑

11月11日　片持ち梁†

メッシュの長さ	要素数	変位	相対誤差	計算者
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	5199	6.29990	-5.55	山本
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	682	5.47288	-17.9	山本
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.9	松田
10	44	1.226075	-81.6	河合

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

Salome-Meca演習_片持ち梁(2021)

↑

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	17530	0.3986	-4.34	山本
1.5	15433	0.396317756757	4.9	進藤
1.6	15900	0.399049	-4.24	河合
1.8	11677	0.404457	-0.03	山口
2	10460	0.394818715517	-5.3	進藤
3	2344	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月25日†

↑

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

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

メッシュ長さ	要素数	変位(mm)	相対誤差(%)	計算者
0.5	604167	-	-	千代岡
0.6	203209	0.42383	0.98	高井
0.7	145234	0.43301	3.22	関合
0.8	140987	0.43006	3.2	岡田
0.9	91974	0.429913	3.18	松田
1.2	24800	0.429777	3.14	青野
1.3	23132	0.42989	3.16	山口
1.4	17617	0.429745	3.13	山本
1.5	15433	0.429844	3.2	進藤
1.6	15900	0.429754	3.13	河合
1.8	11677	0.42962	3.1	山口
2	10460	0.429605	3.1	進藤
3	2486	0.429217	3.0	山本
4	1453	0.4293	3.02	関合
5	431	0.427885	2.69	千代岡
6	360	0.4282	2.78	高井
7	196	0.42606	2.25	青野
8	104	0.42631	2.3	岡田
9	81	0.42513	2.03	松田
10	78	0.424466	1.8	河合

↑

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

異方性1次のデータ

メッシュ長さ	要素数	変位(mm)	相対誤差(%)	計算者
0.5	604167	0.50919	3.56	千代岡
0.6	203209	0.50472	2.6	高井
0.7	145234	0.5036	2.42	関合
0.8	140987	0.50283	2.3	岡田
0.9	91974	0.50053	1.8	松田
1.2	24800	0.48739	-0.9	青野
1.3	23132	0.48841	-0.67	山口
1.4	17617	0.48403	-1.56	山本
1.5	15433	0.48202	-2.0	進藤
1.6	15900	0.48329	-1.7	河合
1.8	11677	0.47855	-2.67	山口
2	10460	0.47906	-2.6	進藤
3	2436	0.42787	-12.98	山本
4	1453	0.42772	-13.02	関合
5	431	0.27364	-44.3	千代岡
6	360	0.33927	-31.0	高井
7	196	0.21363	-58.5	青野
8	104	0.22574	-54.1	岡田
9	81	0.22750	-53.7	松田
10	78	0.20327	-58.7	河合

↑

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

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

メッシュ	要素数	変位	相対誤差	測定者
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.086037	-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	河合

↑

＜メモ＞†

↑

UNIXコマンド†

↑

gnuplot†

↑

inkscape†

inkscapeでDesign
凡例のみを消してグラフを残したいときは、「ノードツールで選択→画面上のパス→分解」　をやることで消すことができる。（選択したあとctr+shift+Kでショートカットできる）
上記の方法だと色もなくなってしまうので色は、「ノードツールで選択→ストロークの塗り」　で色付け作業をする。

↑

LaTex†

LaTex入門コマンド
LaTexコマンド一覧（直コマンド）
単位の累乗（cm^2,mm^2）をするときは直接cm^2やmm^2で書くのではなく式として認識させるために$で囲う必要がある。
図の中に文字の大きさの変え方
pngの貼り方 \includegraphics[width=30mm]{.png}　（数字は画像の幅の拡大・縮小）
pdfの貼り方　\input{.pdf_tex}
\if0 \fi　でコメントアウトしたい行を囲むと複数行のコメントアウトができる。

↑

LibreOfficeImpress†

数式の作り方

↑

Salome-Meca†

code-asterの説明 (Asterstudyで書かれているコードの説明が書かれている。ただ英語で書かれているので、翻訳する必要あり)
Salome-Mecaの使用法解説
Salome-Meca非線形解析　（古い）
Salome-Meca_降伏点検証
Salome-フランス語対策
折り返しばねの計算

↑

春課題†

KTHM表記でファイル作成

K：片持ち　T：単純　H：引張　M：曲げ

春課題弾塑性解析2023

全角/半角の切り替えは、「無変換」（左手親指）でもできる
コピペは、左クリックで領域選択したら、スクロールボタンで貼り付け
単語選択は、ダブルクリックでも可（トリプルクリックだと行選択かな）

青野の卒論日誌 の履歴(No.26)