山本の卒論日誌

3/11 ケーブルデータ†

卒論に用いたデータ†

k2 日付が書いてあるやつに入っている

3/6†

研究開始 解析できるPCが今無いため，とりあえずモデルを作成する 半径75mm 長さ50mで作成 ピン固定でとりあえず作成

133mも随時やってみる

2/21†

全国大会提出に向けてやること

ピン固定ピン固定で解析すること それが理論値に近づくのかを検証する

2/21†

研究成果をスライドを用いて，残しておく

改善すべきところ

• ピン接合の問題
• 境界条件が固定，ローラー支承で今回解析しているため理論式と合わなかった可能性がある
• 振動解析の解析方法

2/20 卒論発表後日†

長く大変だった卒論発表会が終わった．発表自体は割とうまくいったと思う． 修論のときはこんな体たらくではいけない

2/7†

2/6†

国土交通省 https://www.mlit.go.jp/road/sisaku/yobohozen/tenken/yobo4_1-2.pdf

振動法 https://www.jasbc.or.jp/gihoudb/2005-2010/kurimoto/2007-04.pdf

斜張橋やニールセンローゼ橋などのケーブル張力の測定には、いわゆる振動法という手法が広く採用されている。この手法はケーブルの振動測定により固有振動数を検出し、ケーブルの曲げ剛性やサグの影響を考慮した換算式を用いて張力を算出するものである。振動測定にはケーブルに直接加速度計を貼付け、得られた加速度の時刻歴変化にスペクトル解析を適用し、卓越振動数を検出する方法が一般的に用いられている。

一般的には目視点検が行われている

時間設定†

14.99sまでケーブルを安定させるために力を与えず，15sのとき載荷線に100kNの力を与え，15.01s~100sまでで卓越した振動を検出する．

50mのときの固有振動数理論値†

f=1.5188Hz

2/5†

50mのときの固有振動数理論値 f=1.5188Hz

健全時モデル 10m腐食

50m腐食 腐食断面0.075m腐食（10m腐食と同じ）

腐食断面0.0375m腐食

荷重変位グラフ（降伏時）　[縦軸 荷重　 横軸 変位]

2/3†

荷重変位グラフ（降伏時）

1/31†

かかっている力の求め方を聞く soku3やるのか聞く

133m10000soku2,133m50000soku2解析終了
荷重変位で弾塑性のグラフを書く

1/18†

断面六角形での長さ133m,100m,50m,10mで作成したところ，10mでしか解析できなかった．
円柱のときは50mでも作成できた．今回は押出で作成したことによる影響によって？→10mでは作成できている

許容誤差を甘くするとできる？→ある程度試したが成果なし

1/17 六面体で作成†

133mでの作成を試みる

mesh出来上がり

1/17†

30%腐食のグラフ降伏部分が見づらい→再度引張力をいじる

MATLABをとりあえずやる（解析で使うやつを作る） インパクトロードの解析では引張方向に今まで100kNかけていたが，3,2MNぐらいかけなければいけない

1/10†

samplemodel 500 50 Fz=1159.2

1/9†

133mのモデル　最大最小75でメッシュを切れば許容誤差を0.1にすることで解析できた．（静的線形解析） 要素数24000程度 →弾塑性解析ではできなかった

メッシュの切り方を変えてみたが，133mモデルの弾塑性解析はできなかった 直径150mmであるため，最小は75から試したがうまくいかなかった．

12/12 中間発表†

球腐食の結果がおかしい（大きな力で引張をしているのに観測点で圧縮が生じている）
球腐食周りの変位（DEPL）を見てみる→warp by vectorでも変形様子を見てみる

12/7　弾塑性解析†

133mモデルで解析を始めるとエラーが出る→線形解析で行った 133m,110mは解析できない（まだ調査中）100mは精密度（RESI_RELA）をいじれば解析できる．

12/6†

PRE_EPSIにチェック→PRE_EPSIの値を1000,1.85✕10^-3，入れてない状態 →固有振動数が変わらなかった？ ケーブル要素（ケーブル初期張力をいれるところがある）とか調べたが，斜張橋とかに使われているケーブルと挙動が違うらしい？だから使えない？
ケーブル要素は１Ｄでしか使えなそうなため，３Ｄでやっている今の解析は無理そう

張力確認†

①振動解析による歪を見た（この時，確認する点は図形の中心）
slice→xyznormalで中心の点を確認する
②①から出たSIGMをσ=Eεの計算し，計算結果で出たε（ひずみ）をPRE_EPSIに入れる
③上記からPRE_EPSIをいれた状態で歪を確認する
ほぼ変化なかった？
④振動数と跳梁の関係の理論式の張力を判別するためにFORC_NODA,REAC_NODAを入れて解析をした
⑤振動解析とFORC_NODA,REAC_NODAを同時に解析することはできなかった~（chiyokatamesi6.med）

メモ：①線形解析のときはFORC_NODA,REAC_NODAを使うことができたため値を見ることができた　　　②FORC_NODAを入れるときはCALC_CHAMPのEXICTのロードを入れないとダメ

ひずみ引張†

張力が入っているかどうかを確認するために，FORC_NODA,REAC_NODAを入れて確認した．
FORC_NODA　節点力　REAC_NODA　節点反力
入れることができたが，計算誤差範囲でエラーが出た→無理やり範囲を値が出るように設定した→解析はできた FORC_NODA,REAC_NODAが入っているので張力が入っている？ ※

AFFE_CHAR_MECAのPRE_EPSIで与える． 最初両端固定のまま，ひずみで引張を与えていたがうまく行かなかったので，片面だけZ軸に自由を与えたら解析できた．
固定端のまま，PRE_EPSIは使えない？ 固定固定のまま，歪で引張を与えないといけないかどうか

秋山さんの張力データ https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?%E7%A7%8B%E5%B1%B1%E3%81%AE%E4%BF%AE%E8%AB%96%E6%97%A5%E8%AA%8C#uf8e5171

基準のモデル （sample005.hdf）

ひずみ引張　（chiyoka.hdf）

PRE_EPSI†

https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?%E5%8D%83%E4%BB%A3%E5%B2%A1%E3%81%AE%E5%8D%92%E8%AB%96%E6%97%A5%E8%AA%8C#c71fb860

引張サンプル†

理論値†

理論値†

健全時†

相対誤差†

腐食結果†

半円腐食 （半径0.075mm）†

30％I型腐食†

固有振動数（Hz）

位置１ 位置２ 位置３
位置４ 位置５

60％I型腐食†

固有振動数（Hz）

位置１ 位置２ 位置３
位置４ 位置５

球腐食0.075mm（半径）†

腐食させる位置は10m,30m,66,5m

位置１（10m）

位置２（30m）

位置３（66.5m）

球腐食0.0375mm（半径）†

位置１

位置２

位置３

円錐腐食†

50% 60% 70%
80% 90%

133mI型腐食†

解析結果の固有振動数に相関性がない
partition有り
3０％腐食 固有振動数

６０％腐食 固有振動数

位置１ 位置２ 位置３ 位置４ 位置５

理論値

２次要素†

60%減少1

30%減少1

10/11†

１通り今までのものを２次要素で解析
同じ腐食度で腐食の形を変え（中抜け），振動数の違いを見る

10/4~やること†

①モードシェイプをグラフに添付
②中抜け腐食→重さなど一緒で固有振動数の違いを確認する
③別のサイズ→新たなモデル，サイズ

1次要素グラフ†

30%位置1
鉛直1次モード 水平1次モード
鉛直2次モード 水平2次モード

30%位置5
鉛直1次モード 水平1次モード
鉛直2次モード 水平2次モード

円錐2次モーメント50%

卒論進捗†

球腐食

球腐食

60%減少1

30%減少1

円錐腐食1

半円腐食

やること 腐食位置の位置の図
円錐の理論値
解析結果のまとめ

9/8†

ピン固定はエラーが出たためダメだった 面固定でやれば解析はできた．
ke-buru2の中に入ってるデータを使えば解析は成功するが，modal analysisはダメだった．

9/5†

一通り試したがいつものエラー（円柱）
円柱ではエラーをだったが，角柱ではできた

固有振動数を求める解析の仕方を見つけた（試し段階）
Asterstudy→operations→add stage with assistant→modal analysis
https://qiita.com/Jun_Tatsuno/items/a9def3216f230e794fc7

8/2†

エラー直らず

7/28†

ピン固定ピン固定で解析をし直したい
今までの解析で設定のミスがいくつかあった

ピン接合でも解析を始めた

直径0.15mI型腐食†

位置1.2.3.5で解析が成功したが4で解析ができなかった またそこから位置3.5と4.5も解析したところ，3.5は成功し4.5は解析できなかった 下のエラーコードが共通して出てきた

30%減少

30％減少

エラーコード†

サンプルモデル†

I型腐食†

30%減少鉛直1次モード

30%減少水平1次モード

60%減少鉛直一次モード

60%減少水平一次モード

7/26†

理論値と解析値の乖離がある 英語の発表の準備

7/19†

英語スライド作成 [#f3155f8f]　　　円錐モデル作成　理論値計算

f=0.02673968

salome上の振動の仕方 単純単純

salomeの解析上半径を0.15mで解析していた（本当は0.075m）（7/26）

7/12†

引張サンプル作成　エラー [#aedac3f1] 逆円錐，半円のモデル作成

7/4　†

理論値の計算がうまくいかない→形的（I型）に理論値を出すのが厳しいため下記の方法でやっていく？ [#s4d7fe42]

　 　サビの無い一様な円柱で理論値を出す（ここがしっかりしていればサビの部分を減少させても信用性が高い）→要素数を変えてそれに近づける→近い要素数の円柱からサビの部分を減少させていく

①雨が原因でサビができた場合は雨の水が下方部に行くので下方部がだんだん細くなっていくような錆び方ではないか？→最終的に中が逆円錐となるでき方
②中の円柱状の断面積がだんだん小さくなっていくようなやり方

引張試験は振動させたまま解析を行うと勘違いしていたためエラー続出→引張試験で解析を行う

〈メートル換算計算〉

３０％減少

６０％減少

https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/yamamoto/sabihankeimodel.odp

「腐食を考慮した斜張橋ケーブルの共振に関する研究」

直径150mm，長さ133mで解析 beamとsolidで解析した． beam解析→tamesi1 solid解析→entyuu2 面固定解析→menkotei1

断面積の30%,60%減少

133mを10等分し，半径0.075mの球をサビとみなし，カットした． 球の大きさや形を変える．

○面固定

○円の中心点固定

点固定にするとうまく行かない→

mesh→wire,number Asterstudy→AFFE MODEL tout POU_D_E

AFFE CALA CERCLE VALE 0.075

GENERALE

●サンプル 100✕100✕1000の片持ちばりを作成し，４等分した．

https://www.str.ce.akita-u.ac

サイト†

振動モードの直交性と自由振動解--モード解析法 http://mechanics.civil.tohoku.ac.jp/bear/nisikozo/s4node5.html#SECTION0410420000000000000000

ハイブリッドＰＣ斜張橋の設計・施工－なぎさ・ブリッジ https://www.psmic.co.jp/technology/report/pdf/pdf_1/TR0101_text.pdf

論文関係†

耐候性鋼の橋梁への適用

橋梁の耐風性確認を目的とした海上での風観測結果の検討

風疲労

空力振動

風圧力

斜張橋ケーブル論文†

斜張橋の支持ケーブルの局部振動の解析 https://www.google.com/url?sa=t&rct=j&q=&esrc=s&source=web&cd=&ved=2ahUKEwiV4OrErer-AhUN1GEKHX56C8AQFnoECAwQAQ&url=https%3A%2F%2Fnagasaki-u.repo.nii.ac.jp%2Findex.php%3Faction%3Dpages_view_main%26active_action%3Drepository_action_common_download%26item_id%3D9419%26item_no%3D1%26attribute_id%3D18%26file_no%3D1%26page_id%3D13%26block_id%3D21&usg=AOvVaw1BmiG7-sObrlDoCWYLxIIc

ケーブルに基づく斜張橋の設計の最適化ブリッジレゾナンスコントロール https://data.smar-conferences.org/SMAR_2019_Proceedings/Inhalt/we.3.c.5.pdf →翻訳論文１

斜張橋のデッキとケーブル間の共振解析隣接するケーブルの結合効果を考慮した橋梁 https://scholarsmine.mst.edu/cgi/viewcontent.cgi?article=3371&context=civarc_enveng_facwork →論文翻訳２

斜張橋のケーブルの風による振動メカニズムに関する研究 https://www.tus.ac.jp/ura/wp-content/uploads/2022/04/075_kimura.pdf

salome関係†

slomemeca使い方 https://www.youtube.com/watch?v=qZxoJVP0L98&list=PLnLQ1ebfuH0YUozTOsxT1WB4IkxiHcY7X&index=2

salomemeca フランス語 http://opencae.gifu-nct.ac.jp/pukiwiki/index.php?Salome-%A5%D5%A5%E9%A5%F3%A5%B9%B8%EC%C2%D0%BA%F6

AFFE_CHAR_MECA一覧 http://opencae.gifu-nct.ac.jp/pukiwiki/index.php?AFFE_CHAR_MECA%B0%EC%CD%F7

asterstudyコマンド https://qiita.com/Jun_Tatsuno/items/f87220bbeb2ce0aaf3a7

unixコマンド https://www.str.ce.akita-u.ac.jp/~gotou/linux/vine.html#

salome-mecaメモ https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?Salome%E3%83%A1%E3%83%A2

及川さん振動†

https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?cmd=read&page=%E5%8F%8A%E5%B7%9D%E3%81%AE%E4%BF%AE%E8%AB%96%E6%97%A5%E8%AA%8C&word=%E5%8F%8A%E5%B7%9D#n44f5ae1

mode解析†

https://www.opencae.or.jp/wp-content/uploads/2019/12/01-modal-analysis.pdf

11/11　片持ち梁†

11/18 単純梁†

異方性　1次要素 　理論値　0.4917(ティモシェンコ梁)

等方性　2次要素 　理論値　0.4167

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

木材　ヤング率:6000 ポアソン比:0.4　 鋼材　ヤング率:206000 ポアソン比:0.3 理論値:0.99

file:///home/kouzou/yamamoto/yamamototex/sibuw/sindoukaiseki2023/sindou2023.tex

Latex†

https://medemanabu.net/latex/latex-commands-list-alphabetical-order/

\pm プラスマイナス　　\mp　マイナスプラス https://takataninote.com/tex/pm.html

複数コメントアウト　\if0 \fi

最新の100件