青野の卒論日誌
をテンプレートにして作成
[
トップ
] [
新規
|
一覧
|
検索
|
最終更新
|
ヘルプ
|
ログイン
]
開始行:
[[青野の卒論日誌]]
-[[青野の修論日誌:https://www.str.ce.akita-u.ac.jp/cgi-bi...
#contents
*2023年 [#n71df4cf]
*進捗報告 [#odd3f7c3]
2/22 毎日追われるように解析と修正の繰り返しの日々だったが...
wikiの更新が全然できてなくて結果ばかりになってしまったが...
（覚えている範囲の手直しなので間違いが出るかも もっとノ...
12/22 とりあえず中間発表終了
12/20 中間発表用のLaTeXとスライドの作成
12/13 データまとめはとりあえず終わったが決定係数$R^2$とは...
12/6 自分のパソコンだと解析に3週間程度かかるので及川さん...
11/29 解析は終わっていたが間違えて1次要素でやっていた。2...
11/22 かじか橋のアーチ部分に異方性用のBOX作成完了。
AsterStudyの変更も多分完了したので今日解析して帰って金曜...
11/15 青山さんの論文を読んだのでこれからアーチに異方性を...
11/8 ワクチンによる体調不良
11/1 ウォルトキングの論文を梁要素で再現してみた
10/25 ウォルトキングの論文を一通り読み終えた
10/18 理論値と大きく離れたので論文を頑張って翻訳中
10/11 ウォルトキングの理論値関連
*東北支部後 [#d82cd903]
とりあえずねじれが卓越した原因がせん断弾性係数Gの影響だっ...
最初の解析は$ G= \frac{E}{10}$とかで作ってやってみるとか...
上が$G = \frac{E}{10}$ 下が$G = \frac{E}{15}$
3/8 思ったより差が出なかったので$G = \frac{E}{5}$ 下が$G ...
3/11 salomeが落ちてた… 何も結果は出てないし、どこで終わっ...
3/15 とりあえずでやったものが下のグラフだが、Gを変えたの...
3/18 床版も全てのGを同じ倍率で変えてみたら下のグラフのよ...
原因を探さなきゃ…
3/25 とりあえず全データを取り終えた。傾きはせん断弾性係数...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
アーチだけ取り出して今は15BOXでアーチの角度変えているけど...
*卒論関連 [#wb454c03]
振動解析（青山さん、七五三さんの続き) 不朽箇所1箇所ごと...
参考:
[[経年による木製アーチ道路橋（かじか橋）の構造性能と健全...
*fitting [#zd0d3c33]
ヤング率の減少のみで逆対象1次を測定値8.69Hzに合わせた場合
,固有振動数(Hz),水平1次モード,鉛直逆対称1次モード,ねじれ...
,アーチ異方性,6.01510,8.69551,9.45167,13.8143
,測定値(2019),2.73,8.69,10.64,15.38
,アーチ等方性,6.64165,8.69488,12.3674,14.9576
ねじれが10.64Hzに合ってくれたらうれしい
等方性だと振動数が高く、異方性だと振動数が低い
ねじれの誤差
アーチ等方性 16.23%
アーチ異方性 12.57%
*振動数の変化 [#pd38588b]
床版は異方性を適用済み
腐朽なしのときの固有振動数
,固有振動数(Hz),水平1次モード,鉛直逆対称1次モード,ねじれ...
,アーチ等方性,9.44718,12.9978,18.4994,21.5656
,アーチ異方性,9.17776,12.3707,15.7861,18.2731
,全体異方性,7.56782,11.1306,13.0727,17.5401
各状態からの減少率
,減少率(%),水平1次モード,鉛直逆対称1次モード,ねじれ逆対称...
,アーチ等方性→アーチ異方性,2.85,4.82,14.67,15.27
,アーチ等方性→全体異方性,19.89,14.37,29.33,18.67
,アーチ異方性→全体異方性,17.54,10.02,17.19,4.01
ねじり剛性とは、$K = G \times J $
Kがねじり剛性、Gは材料のせん断弾性係数、Jがサン・ブナンの...
$ G = \frac{E}{2 \times (1 + \nu )} $
$\nu$ はポアソン比
*かじか橋の感度解析 [#kbf44a96]
アーチに異方性を適用して振動解析を行った。アーチが腐朽し...
**感度やり直し [#hfac7b11]
***アーチ部分のA1〜A5までを腐朽させていったときの線形回帰...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
***桁部分のC1〜C4までを腐朽させていったときの線形回帰の傾...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
**かじか橋を全部異方性にした場合 [#tc50de43]
***アーチ部分のA1〜A5までを腐朽させていったときの線形回帰...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
***桁部分のC1〜C4までを腐朽させていったときの線形回帰の傾...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
**感度解析結果 [#v99adb0a]
***腐朽箇所 [#sa57fcdc]
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
***アーチ部分のA1〜A5までを腐朽させていったときの線形回帰...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
***桁部分のC1〜C4までを腐朽させていったときの線形回帰の傾...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
縦軸:固有振動数（Hz）、横軸:ヤング率の残存率
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
及川さんの英語論文より５％程度固有振動数が減少したかな？
決定係数$R^2$は、数値予測のモデルの良さを測る指標の一つ。...
**アーチの固有振動数 [#j941d756]
-[[円弧アーチの固有振動数:http://library.jsce.or.jp/jsce/...
-[[３ヒンジアーチの自由振動解法:http://library.jsce.or.jp...
-[[モードの種類:https://www.jstage.jst.go.jp/article/stru...
-[[ウォルトキング(waltking):https://www.str.ce.akita-u.ac...
-[[3ヒンジアーチの幾何学的非線形解析:https://www.str.ce.a...
**ウォルトキングの理論値 [#pa66448f]
論文を読んでからやり直し
(E=9.6*10^9N/m^2 , b=0.1m , h=0.65m ,RHO=459Kg/m^3 ,弧長$...
,$f/\ell$=1/8,対称1次,対称2次,逆対称1次,逆対称2次
,解析値(Hz),34.4245,48.1666,18.2823,76.3628
,理論値(Hz),40.012,117.369,18.6724,77.3571
,誤差(%),13.96,58.96,2.09,1.28
(E=9.6*10^9N/m^2 , b=0.1m , h=0.65m ,RHO=459Kg/m^3 ,弧長$...
,$f/\ell$=1/5,対称1次,対称2次,逆対称1次,逆対称2次
,解析値(Hz),36.3346,65.0017,16.0021,71.0228
,理論値(Hz),37.881,111.081,15.665,71.206
,誤差(%),4.08,41.48,2.15,0.26
(E=9.6*10^9N/m^2 , b=0.1m , h=0.61m ,RHO=459Kg/m^3 ,弧長$...
,$f/\ell$=1/16,対称1次,対称2次,逆対称1次,逆対称2次
,解析値(Hz),20.5211,46.8057,20.3054,82.1309
,理論値(Hz),44.014,117.369,18.672,82.692
,誤差(%),53.38,60.12,8.75,0.68
対称２次は全部上手く行ってないが、その他のモードに関して...
特に逆対称２次に関しては誤差がかなり小さい解析ができてい...
load1 = AFFE_CHAR_MECA(DDL_IMPO=(_F(DRX=0.0,
DRZ=0.0,
DX=0.0,
DY=0.0,
DZ=0.0,
GROUP_NO=('koteiL', )),
_F(DRX=0.0,
DRZ=0.0,
DX=0.0,
DY=0.0,
DZ=0.0,
GROUP_NO=('koteiR', )...
MODELE=model)
構造力学集P269にあるウォルトキング（waltking）の円弧2ヒン...
解析は梁要素で確認を行った。（E=9.6*10^9N/m^2 , b=0.625m ...
,f/l=1/8,対称1次,対称2次,逆対称1次,逆対称2次
,解析値(Hz),4.46454,11.7934,4.53596,18.9013
,理論値(Hz),9.84918,28.8909,4.59629,19.0417
,誤差(%),54.67,59.18,1.31,0.74
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
*かじか橋 [#x80ca80f]
**異方性の確認 [#l248222d]
直線梁を作成して等方性と異方性の振動の違いを確認。
B=10mm,H=5mm,L=500,RHO=4.59/mm^3,E=9600N/mm^2,NU=0.4 メッ...
,,鉛直1次,鉛直2次,水平1次,水平2次
,等方性(Hz),41.4685,165.794,82.8916,330.865
,異方性(Hz),41.4215,165.046,82.5045,328.842
,等方性と異方性の差(%),0.11,0.45,0.47,0.61
**アーチ部分に異方性を入れる [#p1802ed0]
かじか橋をやるということで5本のアーチでやっていたが、アー...
とりあえず単純梁を作成して異方性にした際の振動数の減少を...
等方性の鉛直振動は固定の仕方を間違っていたのでやり直し(8/4)
全データ取り終わり。等方性より異方性のほうが振動数が減少...
次はアーチを作成して同じようになるのかを確認。（8/9）
アーチにオイラー角を用いて振動解析してみた。等方性だと鉛...
とりあえずアーチのみのデータは取り終わった。1次要素・2次...
構造力学公式集にあった2ヒンジアーチの固有振動数の式から1...
構造力学集にあったウォルトキングの図から振動数を求めたら...
参考
-[[根本の卒論日誌:https://www.str.ce.akita-u.ac.jp/cgi-bi...
ここから下はアーチにオイラー角で異方性を15BOXで与えたとき...
横軸：要素数 縦軸：振動数
鉛直1次振動
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
鉛直2次振動
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
水平1次振動
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
水平2次振動
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
ここから下は単純梁の記録（E=9.6GPa,RHO=4.59ton/mm^3, b=20...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/1ji...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/1ji...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/2ji...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/2ji...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/1ji...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/1ji...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/2ji...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/2ji...
縦軸：時間（分） 横軸：要素数
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/jik...
縦軸：振動数 横軸：要素数
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/sui...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/iho...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/tou...
-メモ
アーチは一つのものとして作成されているため、場所ごと...
-boxの適用方法
アーチを選択→点を作成（左から２番目のものを選ぶと角の...
以上によって面で区切りとなってboxを作成することができ...
**オイラー角 [#m46aac50]
もしかしたら青山さんのモデルを異方性にしていくこともする...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/xyh...
Z軸回転のみ（部材がXY平面上にあるとき）のオイラー角の適用...
多分一つの回転のみなら簡単にできるのかもしれない
→x軸回転のみもできたのでZ軸またはX軸方向の回転はできる（y...
追加
boxを斜め方向（どの平面にも触れないように）作成し、オイラ...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/nan...
片持ち梁を上の図のように作成して、曲げの応力－ひずみから...
オイラー角の回転はZ軸→X’軸→Z’’軸の順番に行っていることに...
オイラー角：座標軸で回転させて向かせたい方向に回転させる...
x軸を回転させる：φ角
y軸を回転させる：θ角
z軸を回転させる：ψ角
（飛行機の場合、上からバンク角、ピッチ角、機首方位角と呼...
回転させる時、原点から軸の正方向を見て右回転させていく。
片持梁の曲げで確認
応力とひずみからヤング率を計算したところ、軸方向ヤング率...
アーチを作成して適用してみる。
-[[salomeでやる時のAstercode:https://code-aster.org/V2/do...
-[[オイラー角とは:https://www.sky-engin.jp/blog/eulerian-...
-[[回転行列:https://www.sky-engin.jp/blog/rotation-matrix...
-[[salomeでの構造要素:https://www.opencae.or.jp/wp-conten...
-[[オイラー角動き:http://irobutsu.a.la9.jp/mybook/ykwkrAM...
*密度による影響 [#wcf71913]
**3分割で試し [#e96699e3]
最初の名前が表面の状態
数字がモード（１：１次モード、２：２次モード）
e(鉛直）s（水平）方向
断面図↓
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/3bu...
横軸：振動数 縦軸：中心の含水率
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/mit...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/mit...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/mit...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/mit...
**10分割含水率変化 [#r7d3ce9b]
-簡単な橋を作成し（10✕20✕1000とか）、ミルフィーユ状に切っ...
最初はサンドイッチみたいなのからスタートして徐々に層数を...
気乾比重は、乾燥させたときの比重のため少し（含水率10％程...
全乾比重は、完全に乾燥させたときの比重（含水率0％）となっ...
スギでモデル作成していく.10✕20✕1000をモデルとし、気乾比重...
-[[木材の全乾比重:https://sooki.co.jp/upload/surveying_it...
-[[木材の比重:https://www.toishi.info/sozai/woods/d.html]]
,種類,密度g/cm3
,バルサ,0.10
,リグナムバイタ,1.25
,キリ,0.23
,アカガシ,0.97
,スギ,0.34
,ヒノキ,0.42
,ケヤキ,0.66
,ミズナラ,0.70
*振動解析 [#j44be043]
**春休み課題（振動解析） [#g422aff6]
図心をZ軸に合わせて固定もZ軸上でのものとする（Boxの中心上...
追記：振動解析班のものは鉛直、水平のモードで分類してなか...
相対誤差まで完成。ほとんど誤差のない結果となった。（5/2）
再解析データ
理論値
鉛直1次：36.036Hz 水平1次：72.073Hz 鉛直2次：144.146Hz...
1次モード
,メッシュの長さ,鉛直1次モード,相対誤差,水平1次モード,相対...
,0.8,36.0307,0.015%,72.0219,0.067%
,0.9,36.0307,0.015%,72.0219,0.067%
,1.0,36.0307,0.015%,72.0219,0.067%
,2.0,36.0310,0.014%,72.0224,0.066%
,3.0,36.0308,0.014%,72.0224,0.066%
,4.0,36.0312,0.013%,72.0226,0.066%
,5.0,36.0313,0.013%,72.0227,0.066%
2次モード
,メッシュの長さ,鉛直2次モード,相対誤差,水平2次モード,相対...
,0.8,144.053,0.065%,287.477,0.283%
,0.9,144.053,0.065%,287.477,0.283%
,1.0,144.053,0.065%,287.477,0.283%
,2.0,144.058,0.061%,287.484,0.280%
,3.0,144.055,0.063%,287.485,0.280%
,4.0,144.060,0.060%,287.488,0.279%
,5.0,144.062,0.058%,287.490,0.278%
振動解析班のデータ１次モード
,メッシュの長さ,要素数,一次モードの固有振動数解析値,相対...
,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
-[[春休み振動解析班:https://www.str.ce.akita-u.ac.jp/cgi-...
-[[振動解析:https://www.str.ce.akita-u.ac.jp/cgi-bin/puki...
ねじれの理論式 $ f= \frac{\lambda}{2\pi\ell}\times \sqrt...
棒材の場合はJ=$I_p$であるため$ f= \frac{\lambda}{2\pi\ell...
**参考文献 [#v9ecce87]
[[構造振動学:https://www.jsme.or.jp/sed/guide/dynamics5.p...
[[固有振動数 理論式:https://cae-jutaku.com/vibration-ver...
日本機械学会のものが見当たらなくなったのでとりあえずおい...
[[既存鋼鈑桁橋のモデルパラメータ事後分布を用いた信頼性評...
[[ベイズ推定による既存構造物数値モデルの不確定性定量化と...
*読んでみたいもの [#x2964178]
（Model Verification & Validation（V&V））
Alvin, K. F., Oberkampf, W. L, Dieger, K., Rutherford, B....
（計算構造力学における不確かさの定量化：モデル検証のため...
*2022年 [#e1aaefb8]
*創造工房実習の内容 [#v784b03e]
**10月21日課題 [#c072c2b0]
タッチタイプ練習
**11月4日 [#vf6be901]
・課題
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/kad...
[[秋田市の月別平均気温(1991～2020):https://www.str.ce.aki...
・練習
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/bbb...
**11月11日 片持ち梁 [#ra69fe4a]
,メッシュの長さ,要素数,変位,相対誤差,計算者
,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,河合
片持ち梁のグラフ：縦軸 変位 横軸 要素数
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/kat...
-[[Salome-Meca演習_片持ち梁(2021)]]
**11/18(金) 単純梁 [#occad30b]
,メッシュ長さ,要素数,変位(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,河合
片持ち梁のグラフ：縦軸 変位 横軸 要素数
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/tan...
**11月25日 [#f01c9dab]
***等方性1次と等方性2次の比較 [#v697f20e]
二次要素の等方性のデータ
,メッシュ長さ,要素数,変位(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,河合
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/kad...
***等方性1次と異方性1次の比較 [#l937a0df]
異方性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,河合
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/kad...
**12/2 木材と鋼材の単純梁 サンドウィッチ [#ta715778]
異方性の単純梁(1次)と2次要素のサンドウィッチ梁（鋼材：等...
[[コマンドファイル:https://www.str.ce.akita-u.ac.jp/~goto...
,メッシュ,要素数,変位,相対誤差,測定者
,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,河合
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/san...
*＜メモ＞ [#w5f3def5]
**UNIXコマンド [#m9b3d34b]
-[[ここのUNIXコマンド:https://www.str.ce.akita-u.ac.jp/~g...
-[[コマンド集:https://g-omr.github.io/unix.html]]
-[[操作系コマンド:https://www.k-tanaka.net/unix/]]
-[[ジャンルごとに細かいもの:http://x68000.q-e-d.net/~68us...
-viにおいて[:w ファイル名]とすると別名保存できる
-[[複数行のコメントアウト:https://qiita.com/kagami_t/item...
1.Ctrl + v。
2.コメントアウト範囲を選択。
3.Shift + i
4.%を入力。
5.Escキーで完了。
-free:メモリの空き容量やメモリの使用量を確認
-top:実行中のプロセスの状態を確認
-swapファイルがあると出てしまったら→la -aで隠しファイルを...
**gnuplot [#cb967f34]
-[[gnuplot:https://qiita.com/python_walker/items/2b4febdb...
-[[gnuplotスクリプトの解説:https://ss.scphys.kyoto-u.ac.j...
-plot 関数 with linespoints (plot 関数 w lp)
-[[目盛り間隔の変更:https://www.nakamuri.info/mw/index.ph...
**inkscape [#w59b9a6b]
-[[inkscapeでDesign:http://inkscapedesign.web.fc2.com/]]
-凡例のみを消してグラフを残したいときは、「ノードツールで...
-上記の方法だと色もなくなってしまうので色は、「ノードツー...
**LaTex [#qa320d06]
-[[LaTex入門コマンド:https://medemanabu.net/latex/latex-c...
-[[LaTexコマンド一覧（直コマンド）:http://www3.otani.ac.j...
-単位の累乗（cm^2,mm^2）をするときは直接cm^2やmm^2で書く...
-[[図の中に文字の大きさの変え方:https://www.str.ce.akita-...
-[[数式の書き方:https://cns-guide.sfc.keio.ac.jp/2001/11/...
-[[余白の調整:https://pyopyopyo.hatenablog.com/entry/2012...
-[[特殊文字:http://www.ic.daito.ac.jp/~mizutani/tex/speci...
-pngの貼り方 \includegraphics[width=30mm]{.png} （数字は...
-pdfの貼り方 \input{.pdf_tex}
-\if0 \fi でコメントアウトしたい行を囲むと複数行のコメン...
-〇〇Gとやると〇〇行に移動してくれる
-/単語で単語検索ができる :noh でハイライト表示を終了
-_を入れたいときは￥＿のように入力すると_が入力できる
- :%s/(変更前)/(変更後)/g で全部(変更後)に置き換えること...
- Ctrl + v で指定した範囲をDeleteで削除できる (shift + v...
-$:行末へ移動
-^:行頭へ移動
-段落を変えたあとに字下げをしたくない場合は\noindentを入...
-字下げを入れたい場合は\indentで字下げをすることができる。
**LibreOfficeImpress [#d1f75924]
-[[数式の作り方:https://asukiaaa.blogspot.com/2021/12/lib...
-[[縦書きテキストボックス:https://www.officeisyours.com/e...
-[[スライドペインの表示:https://ask.libreoffice.org/t/top...
**Salome-Meca [#a3c504d8]
-AstersutdyのMemoryの数値を上げることで解析に与えるMemory...
-[[code-asterの説明:http://dugi6514.odns.fr/doc/v15/en/in...
-[[Salome-Mecaの使用法解説:http://opencae.gifu-nct.ac.jp/...
-[[Salome-Meca非線形解析:https://sites.google.com/site/co...
-[[Salome-Meca_降伏点検証 :https://www.str.ce.akita-u.ac....
-[[Salome-フランス語対策:http://opencae.gifu-nct.ac.jp/pu...
-[[折り返しばねの計算:http://freeplanets.ship.jp/Numerica...
-[[salomemecaとcodeaster概要:https://www.opencae.or.jp/wp...
-[[線形動解析:https://www.opencae.or.jp/wp-content/upload...
-[[CODE-ASTER ASSEMBLAGE:https://code-aster.org/V2/doc/v1...
[m]単位でモデリングするときは，密度の単位はkg/m3を用いる。
GenerateVectorsは2Dの結果を可視化するときに適用するものら...
Normalmodesanimationrealは鉛直or水平の振動を見ることがで...
Normalmodesanimationcomplexは鉛直・水平振動を同時に起こせ...
*自分用メモ [#j7d77527]
-[[バックアップ関連:https://www.str.ce.akita-u.ac.jp/cgi-...
-画像が大きいとき、画像をGIMPで開き.pngから.jpgに変えるだ...
POUTRE 梁（ビーム）
COQUE シェル
AXE_POUTRE ビーム軸
coq_pou このオプションにより、船体内に格子状の部品を梁状...
http://www-mdp.eng.cam.ac.uk/web/CD/engapps/aster_docs/UD...
5/15~5/26 教育実習
5/29 研究再開
全角/半角の切り替えは、「無変換」（左手親指）でもできる
コピペは、左クリックで領域選択したら、スクロールボタンで...
単語選択は、ダブルクリックで可（トリプルクリックだと行選...
*学会メモ [#wb459629]
**東北支部(3/2) [#d8c0474c]
初めての学会発表で東北支部(岩手大学)で発表したが、やっぱ...
質問の意図が上手く汲み取れなくて違う回答をしたあとに質問...
相手の質問内容を上手く理解できるようになる必要がある。
また、他の人の発表が全然理解できなかったのでもっと理解で...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
&online;
終了行:
[[青野の卒論日誌]]
-[[青野の修論日誌:https://www.str.ce.akita-u.ac.jp/cgi-bi...
#contents
*2023年 [#n71df4cf]
*進捗報告 [#odd3f7c3]
2/22 毎日追われるように解析と修正の繰り返しの日々だったが...
wikiの更新が全然できてなくて結果ばかりになってしまったが...
（覚えている範囲の手直しなので間違いが出るかも もっとノ...
12/22 とりあえず中間発表終了
12/20 中間発表用のLaTeXとスライドの作成
12/13 データまとめはとりあえず終わったが決定係数$R^2$とは...
12/6 自分のパソコンだと解析に3週間程度かかるので及川さん...
11/29 解析は終わっていたが間違えて1次要素でやっていた。2...
11/22 かじか橋のアーチ部分に異方性用のBOX作成完了。
AsterStudyの変更も多分完了したので今日解析して帰って金曜...
11/15 青山さんの論文を読んだのでこれからアーチに異方性を...
11/8 ワクチンによる体調不良
11/1 ウォルトキングの論文を梁要素で再現してみた
10/25 ウォルトキングの論文を一通り読み終えた
10/18 理論値と大きく離れたので論文を頑張って翻訳中
10/11 ウォルトキングの理論値関連
*東北支部後 [#d82cd903]
とりあえずねじれが卓越した原因がせん断弾性係数Gの影響だっ...
最初の解析は$ G= \frac{E}{10}$とかで作ってやってみるとか...
上が$G = \frac{E}{10}$ 下が$G = \frac{E}{15}$
3/8 思ったより差が出なかったので$G = \frac{E}{5}$ 下が$G ...
3/11 salomeが落ちてた… 何も結果は出てないし、どこで終わっ...
3/15 とりあえずでやったものが下のグラフだが、Gを変えたの...
3/18 床版も全てのGを同じ倍率で変えてみたら下のグラフのよ...
原因を探さなきゃ…
3/25 とりあえず全データを取り終えた。傾きはせん断弾性係数...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
アーチだけ取り出して今は15BOXでアーチの角度変えているけど...
*卒論関連 [#wb454c03]
振動解析（青山さん、七五三さんの続き) 不朽箇所1箇所ごと...
参考:
[[経年による木製アーチ道路橋（かじか橋）の構造性能と健全...
*fitting [#zd0d3c33]
ヤング率の減少のみで逆対象1次を測定値8.69Hzに合わせた場合
,固有振動数(Hz),水平1次モード,鉛直逆対称1次モード,ねじれ...
,アーチ異方性,6.01510,8.69551,9.45167,13.8143
,測定値(2019),2.73,8.69,10.64,15.38
,アーチ等方性,6.64165,8.69488,12.3674,14.9576
ねじれが10.64Hzに合ってくれたらうれしい
等方性だと振動数が高く、異方性だと振動数が低い
ねじれの誤差
アーチ等方性 16.23%
アーチ異方性 12.57%
*振動数の変化 [#pd38588b]
床版は異方性を適用済み
腐朽なしのときの固有振動数
,固有振動数(Hz),水平1次モード,鉛直逆対称1次モード,ねじれ...
,アーチ等方性,9.44718,12.9978,18.4994,21.5656
,アーチ異方性,9.17776,12.3707,15.7861,18.2731
,全体異方性,7.56782,11.1306,13.0727,17.5401
各状態からの減少率
,減少率(%),水平1次モード,鉛直逆対称1次モード,ねじれ逆対称...
,アーチ等方性→アーチ異方性,2.85,4.82,14.67,15.27
,アーチ等方性→全体異方性,19.89,14.37,29.33,18.67
,アーチ異方性→全体異方性,17.54,10.02,17.19,4.01
ねじり剛性とは、$K = G \times J $
Kがねじり剛性、Gは材料のせん断弾性係数、Jがサン・ブナンの...
$ G = \frac{E}{2 \times (1 + \nu )} $
$\nu$ はポアソン比
*かじか橋の感度解析 [#kbf44a96]
アーチに異方性を適用して振動解析を行った。アーチが腐朽し...
**感度やり直し [#hfac7b11]
***アーチ部分のA1〜A5までを腐朽させていったときの線形回帰...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
***桁部分のC1〜C4までを腐朽させていったときの線形回帰の傾...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
**かじか橋を全部異方性にした場合 [#tc50de43]
***アーチ部分のA1〜A5までを腐朽させていったときの線形回帰...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
***桁部分のC1〜C4までを腐朽させていったときの線形回帰の傾...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
**感度解析結果 [#v99adb0a]
***腐朽箇所 [#sa57fcdc]
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
***アーチ部分のA1〜A5までを腐朽させていったときの線形回帰...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
***桁部分のC1〜C4までを腐朽させていったときの線形回帰の傾...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
縦軸:固有振動数（Hz）、横軸:ヤング率の残存率
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
及川さんの英語論文より５％程度固有振動数が減少したかな？
決定係数$R^2$は、数値予測のモデルの良さを測る指標の一つ。...
**アーチの固有振動数 [#j941d756]
-[[円弧アーチの固有振動数:http://library.jsce.or.jp/jsce/...
-[[３ヒンジアーチの自由振動解法:http://library.jsce.or.jp...
-[[モードの種類:https://www.jstage.jst.go.jp/article/stru...
-[[ウォルトキング(waltking):https://www.str.ce.akita-u.ac...
-[[3ヒンジアーチの幾何学的非線形解析:https://www.str.ce.a...
**ウォルトキングの理論値 [#pa66448f]
論文を読んでからやり直し
(E=9.6*10^9N/m^2 , b=0.1m , h=0.65m ,RHO=459Kg/m^3 ,弧長$...
,$f/\ell$=1/8,対称1次,対称2次,逆対称1次,逆対称2次
,解析値(Hz),34.4245,48.1666,18.2823,76.3628
,理論値(Hz),40.012,117.369,18.6724,77.3571
,誤差(%),13.96,58.96,2.09,1.28
(E=9.6*10^9N/m^2 , b=0.1m , h=0.65m ,RHO=459Kg/m^3 ,弧長$...
,$f/\ell$=1/5,対称1次,対称2次,逆対称1次,逆対称2次
,解析値(Hz),36.3346,65.0017,16.0021,71.0228
,理論値(Hz),37.881,111.081,15.665,71.206
,誤差(%),4.08,41.48,2.15,0.26
(E=9.6*10^9N/m^2 , b=0.1m , h=0.61m ,RHO=459Kg/m^3 ,弧長$...
,$f/\ell$=1/16,対称1次,対称2次,逆対称1次,逆対称2次
,解析値(Hz),20.5211,46.8057,20.3054,82.1309
,理論値(Hz),44.014,117.369,18.672,82.692
,誤差(%),53.38,60.12,8.75,0.68
対称２次は全部上手く行ってないが、その他のモードに関して...
特に逆対称２次に関しては誤差がかなり小さい解析ができてい...
load1 = AFFE_CHAR_MECA(DDL_IMPO=(_F(DRX=0.0,
DRZ=0.0,
DX=0.0,
DY=0.0,
DZ=0.0,
GROUP_NO=('koteiL', )),
_F(DRX=0.0,
DRZ=0.0,
DX=0.0,
DY=0.0,
DZ=0.0,
GROUP_NO=('koteiR', )...
MODELE=model)
構造力学集P269にあるウォルトキング（waltking）の円弧2ヒン...
解析は梁要素で確認を行った。（E=9.6*10^9N/m^2 , b=0.625m ...
,f/l=1/8,対称1次,対称2次,逆対称1次,逆対称2次
,解析値(Hz),4.46454,11.7934,4.53596,18.9013
,理論値(Hz),9.84918,28.8909,4.59629,19.0417
,誤差(%),54.67,59.18,1.31,0.74
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
*かじか橋 [#x80ca80f]
**異方性の確認 [#l248222d]
直線梁を作成して等方性と異方性の振動の違いを確認。
B=10mm,H=5mm,L=500,RHO=4.59/mm^3,E=9600N/mm^2,NU=0.4 メッ...
,,鉛直1次,鉛直2次,水平1次,水平2次
,等方性(Hz),41.4685,165.794,82.8916,330.865
,異方性(Hz),41.4215,165.046,82.5045,328.842
,等方性と異方性の差(%),0.11,0.45,0.47,0.61
**アーチ部分に異方性を入れる [#p1802ed0]
かじか橋をやるということで5本のアーチでやっていたが、アー...
とりあえず単純梁を作成して異方性にした際の振動数の減少を...
等方性の鉛直振動は固定の仕方を間違っていたのでやり直し(8/4)
全データ取り終わり。等方性より異方性のほうが振動数が減少...
次はアーチを作成して同じようになるのかを確認。（8/9）
アーチにオイラー角を用いて振動解析してみた。等方性だと鉛...
とりあえずアーチのみのデータは取り終わった。1次要素・2次...
構造力学公式集にあった2ヒンジアーチの固有振動数の式から1...
構造力学集にあったウォルトキングの図から振動数を求めたら...
参考
-[[根本の卒論日誌:https://www.str.ce.akita-u.ac.jp/cgi-bi...
ここから下はアーチにオイラー角で異方性を15BOXで与えたとき...
横軸：要素数 縦軸：振動数
鉛直1次振動
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
鉛直2次振動
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
水平1次振動
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
水平2次振動
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
ここから下は単純梁の記録（E=9.6GPa,RHO=4.59ton/mm^3, b=20...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/1ji...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/1ji...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/2ji...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/2ji...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/1ji...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/1ji...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/2ji...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/2ji...
縦軸：時間（分） 横軸：要素数
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/jik...
縦軸：振動数 横軸：要素数
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/sui...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/iho...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/tou...
-メモ
アーチは一つのものとして作成されているため、場所ごと...
-boxの適用方法
アーチを選択→点を作成（左から２番目のものを選ぶと角の...
以上によって面で区切りとなってboxを作成することができ...
**オイラー角 [#m46aac50]
もしかしたら青山さんのモデルを異方性にしていくこともする...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/xyh...
Z軸回転のみ（部材がXY平面上にあるとき）のオイラー角の適用...
多分一つの回転のみなら簡単にできるのかもしれない
→x軸回転のみもできたのでZ軸またはX軸方向の回転はできる（y...
追加
boxを斜め方向（どの平面にも触れないように）作成し、オイラ...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/nan...
片持ち梁を上の図のように作成して、曲げの応力－ひずみから...
オイラー角の回転はZ軸→X’軸→Z’’軸の順番に行っていることに...
オイラー角：座標軸で回転させて向かせたい方向に回転させる...
x軸を回転させる：φ角
y軸を回転させる：θ角
z軸を回転させる：ψ角
（飛行機の場合、上からバンク角、ピッチ角、機首方位角と呼...
回転させる時、原点から軸の正方向を見て右回転させていく。
片持梁の曲げで確認
応力とひずみからヤング率を計算したところ、軸方向ヤング率...
アーチを作成して適用してみる。
-[[salomeでやる時のAstercode:https://code-aster.org/V2/do...
-[[オイラー角とは:https://www.sky-engin.jp/blog/eulerian-...
-[[回転行列:https://www.sky-engin.jp/blog/rotation-matrix...
-[[salomeでの構造要素:https://www.opencae.or.jp/wp-conten...
-[[オイラー角動き:http://irobutsu.a.la9.jp/mybook/ykwkrAM...
*密度による影響 [#wcf71913]
**3分割で試し [#e96699e3]
最初の名前が表面の状態
数字がモード（１：１次モード、２：２次モード）
e(鉛直）s（水平）方向
断面図↓
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/3bu...
横軸：振動数 縦軸：中心の含水率
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/mit...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/mit...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/mit...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/mit...
**10分割含水率変化 [#r7d3ce9b]
-簡単な橋を作成し（10✕20✕1000とか）、ミルフィーユ状に切っ...
最初はサンドイッチみたいなのからスタートして徐々に層数を...
気乾比重は、乾燥させたときの比重のため少し（含水率10％程...
全乾比重は、完全に乾燥させたときの比重（含水率0％）となっ...
スギでモデル作成していく.10✕20✕1000をモデルとし、気乾比重...
-[[木材の全乾比重:https://sooki.co.jp/upload/surveying_it...
-[[木材の比重:https://www.toishi.info/sozai/woods/d.html]]
,種類,密度g/cm3
,バルサ,0.10
,リグナムバイタ,1.25
,キリ,0.23
,アカガシ,0.97
,スギ,0.34
,ヒノキ,0.42
,ケヤキ,0.66
,ミズナラ,0.70
*振動解析 [#j44be043]
**春休み課題（振動解析） [#g422aff6]
図心をZ軸に合わせて固定もZ軸上でのものとする（Boxの中心上...
追記：振動解析班のものは鉛直、水平のモードで分類してなか...
相対誤差まで完成。ほとんど誤差のない結果となった。（5/2）
再解析データ
理論値
鉛直1次：36.036Hz 水平1次：72.073Hz 鉛直2次：144.146Hz...
1次モード
,メッシュの長さ,鉛直1次モード,相対誤差,水平1次モード,相対...
,0.8,36.0307,0.015%,72.0219,0.067%
,0.9,36.0307,0.015%,72.0219,0.067%
,1.0,36.0307,0.015%,72.0219,0.067%
,2.0,36.0310,0.014%,72.0224,0.066%
,3.0,36.0308,0.014%,72.0224,0.066%
,4.0,36.0312,0.013%,72.0226,0.066%
,5.0,36.0313,0.013%,72.0227,0.066%
2次モード
,メッシュの長さ,鉛直2次モード,相対誤差,水平2次モード,相対...
,0.8,144.053,0.065%,287.477,0.283%
,0.9,144.053,0.065%,287.477,0.283%
,1.0,144.053,0.065%,287.477,0.283%
,2.0,144.058,0.061%,287.484,0.280%
,3.0,144.055,0.063%,287.485,0.280%
,4.0,144.060,0.060%,287.488,0.279%
,5.0,144.062,0.058%,287.490,0.278%
振動解析班のデータ１次モード
,メッシュの長さ,要素数,一次モードの固有振動数解析値,相対...
,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
-[[春休み振動解析班:https://www.str.ce.akita-u.ac.jp/cgi-...
-[[振動解析:https://www.str.ce.akita-u.ac.jp/cgi-bin/puki...
ねじれの理論式 $ f= \frac{\lambda}{2\pi\ell}\times \sqrt...
棒材の場合はJ=$I_p$であるため$ f= \frac{\lambda}{2\pi\ell...
**参考文献 [#v9ecce87]
[[構造振動学:https://www.jsme.or.jp/sed/guide/dynamics5.p...
[[固有振動数 理論式:https://cae-jutaku.com/vibration-ver...
日本機械学会のものが見当たらなくなったのでとりあえずおい...
[[既存鋼鈑桁橋のモデルパラメータ事後分布を用いた信頼性評...
[[ベイズ推定による既存構造物数値モデルの不確定性定量化と...
*読んでみたいもの [#x2964178]
（Model Verification & Validation（V&V））
Alvin, K. F., Oberkampf, W. L, Dieger, K., Rutherford, B....
（計算構造力学における不確かさの定量化：モデル検証のため...
*2022年 [#e1aaefb8]
*創造工房実習の内容 [#v784b03e]
**10月21日課題 [#c072c2b0]
タッチタイプ練習
**11月4日 [#vf6be901]
・課題
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/kad...
[[秋田市の月別平均気温(1991～2020):https://www.str.ce.aki...
・練習
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/bbb...
**11月11日 片持ち梁 [#ra69fe4a]
,メッシュの長さ,要素数,変位,相対誤差,計算者
,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,河合
片持ち梁のグラフ：縦軸 変位 横軸 要素数
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/kat...
-[[Salome-Meca演習_片持ち梁(2021)]]
**11/18(金) 単純梁 [#occad30b]
,メッシュ長さ,要素数,変位(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,河合
片持ち梁のグラフ：縦軸 変位 横軸 要素数
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/tan...
**11月25日 [#f01c9dab]
***等方性1次と等方性2次の比較 [#v697f20e]
二次要素の等方性のデータ
,メッシュ長さ,要素数,変位(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,河合
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/kad...
***等方性1次と異方性1次の比較 [#l937a0df]
異方性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,河合
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/kad...
**12/2 木材と鋼材の単純梁 サンドウィッチ [#ta715778]
異方性の単純梁(1次)と2次要素のサンドウィッチ梁（鋼材：等...
[[コマンドファイル:https://www.str.ce.akita-u.ac.jp/~goto...
,メッシュ,要素数,変位,相対誤差,測定者
,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,河合
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/aono/san...
*＜メモ＞ [#w5f3def5]
**UNIXコマンド [#m9b3d34b]
-[[ここのUNIXコマンド:https://www.str.ce.akita-u.ac.jp/~g...
-[[コマンド集:https://g-omr.github.io/unix.html]]
-[[操作系コマンド:https://www.k-tanaka.net/unix/]]
-[[ジャンルごとに細かいもの:http://x68000.q-e-d.net/~68us...
-viにおいて[:w ファイル名]とすると別名保存できる
-[[複数行のコメントアウト:https://qiita.com/kagami_t/item...
1.Ctrl + v。
2.コメントアウト範囲を選択。
3.Shift + i
4.%を入力。
5.Escキーで完了。
-free:メモリの空き容量やメモリの使用量を確認
-top:実行中のプロセスの状態を確認
-swapファイルがあると出てしまったら→la -aで隠しファイルを...
**gnuplot [#cb967f34]
-[[gnuplot:https://qiita.com/python_walker/items/2b4febdb...
-[[gnuplotスクリプトの解説:https://ss.scphys.kyoto-u.ac.j...
-plot 関数 with linespoints (plot 関数 w lp)
-[[目盛り間隔の変更:https://www.nakamuri.info/mw/index.ph...
**inkscape [#w59b9a6b]
-[[inkscapeでDesign:http://inkscapedesign.web.fc2.com/]]
-凡例のみを消してグラフを残したいときは、「ノードツールで...
-上記の方法だと色もなくなってしまうので色は、「ノードツー...
**LaTex [#qa320d06]
-[[LaTex入門コマンド:https://medemanabu.net/latex/latex-c...
-[[LaTexコマンド一覧（直コマンド）:http://www3.otani.ac.j...
-単位の累乗（cm^2,mm^2）をするときは直接cm^2やmm^2で書く...
-[[図の中に文字の大きさの変え方:https://www.str.ce.akita-...
-[[数式の書き方:https://cns-guide.sfc.keio.ac.jp/2001/11/...
-[[余白の調整:https://pyopyopyo.hatenablog.com/entry/2012...
-[[特殊文字:http://www.ic.daito.ac.jp/~mizutani/tex/speci...
-pngの貼り方 \includegraphics[width=30mm]{.png} （数字は...
-pdfの貼り方 \input{.pdf_tex}
-\if0 \fi でコメントアウトしたい行を囲むと複数行のコメン...
-〇〇Gとやると〇〇行に移動してくれる
-/単語で単語検索ができる :noh でハイライト表示を終了
-_を入れたいときは￥＿のように入力すると_が入力できる
- :%s/(変更前)/(変更後)/g で全部(変更後)に置き換えること...
- Ctrl + v で指定した範囲をDeleteで削除できる (shift + v...
-$:行末へ移動
-^:行頭へ移動
-段落を変えたあとに字下げをしたくない場合は\noindentを入...
-字下げを入れたい場合は\indentで字下げをすることができる。
**LibreOfficeImpress [#d1f75924]
-[[数式の作り方:https://asukiaaa.blogspot.com/2021/12/lib...
-[[縦書きテキストボックス:https://www.officeisyours.com/e...
-[[スライドペインの表示:https://ask.libreoffice.org/t/top...
**Salome-Meca [#a3c504d8]
-AstersutdyのMemoryの数値を上げることで解析に与えるMemory...
-[[code-asterの説明:http://dugi6514.odns.fr/doc/v15/en/in...
-[[Salome-Mecaの使用法解説:http://opencae.gifu-nct.ac.jp/...
-[[Salome-Meca非線形解析:https://sites.google.com/site/co...
-[[Salome-Meca_降伏点検証 :https://www.str.ce.akita-u.ac....
-[[Salome-フランス語対策:http://opencae.gifu-nct.ac.jp/pu...
-[[折り返しばねの計算:http://freeplanets.ship.jp/Numerica...
-[[salomemecaとcodeaster概要:https://www.opencae.or.jp/wp...
-[[線形動解析:https://www.opencae.or.jp/wp-content/upload...
-[[CODE-ASTER ASSEMBLAGE:https://code-aster.org/V2/doc/v1...
[m]単位でモデリングするときは，密度の単位はkg/m3を用いる。
GenerateVectorsは2Dの結果を可視化するときに適用するものら...
Normalmodesanimationrealは鉛直or水平の振動を見ることがで...
Normalmodesanimationcomplexは鉛直・水平振動を同時に起こせ...
*自分用メモ [#j7d77527]
-[[バックアップ関連:https://www.str.ce.akita-u.ac.jp/cgi-...
-画像が大きいとき、画像をGIMPで開き.pngから.jpgに変えるだ...
POUTRE 梁（ビーム）
COQUE シェル
AXE_POUTRE ビーム軸
coq_pou このオプションにより、船体内に格子状の部品を梁状...
http://www-mdp.eng.cam.ac.uk/web/CD/engapps/aster_docs/UD...
5/15~5/26 教育実習
5/29 研究再開
全角/半角の切り替えは、「無変換」（左手親指）でもできる
コピペは、左クリックで領域選択したら、スクロールボタンで...
単語選択は、ダブルクリックで可（トリプルクリックだと行選...
*学会メモ [#wb459629]
**東北支部(3/2) [#d8c0474c]
初めての学会発表で東北支部(岩手大学)で発表したが、やっぱ...
質問の意図が上手く汲み取れなくて違う回答をしたあとに質問...
相手の質問内容を上手く理解できるようになる必要がある。
また、他の人の発表が全然理解できなかったのでもっと理解で...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/aono/wik...
&online;
ページ名:
![[PukiWiki]](image/pukiwiki.png "[PukiWiki]")
