研究

研究テーマ

やること

進捗

由利橋について

由利橋はH23.2.22に施工が始まり、H25.1.31に完成した。

橋長190.5mの不等径間の2径間連続斜張橋(鋼床版2主箱桁)である。

秋田県由利本荘市の子安川にまたがる形で位置している。

由利本荘市:面積1209.59km^2,人工:74707人

小吉川:一級河川,全長61km,流域面積1190km^2,小洪水は毎年起こっている

予想可能な災害…土砂崩れ(出羽丘陵)、氾濫及び洪水(小吉川)、積雪、暴風

由利本荘市ハザードマップからわかる通り、河川や水に関する災害が由利橋付近で起こる可能性が高い。

考え得る危険事項…交通事故(追突・接触)、揺れ(地震・風・共振)、落雷、積雪、津波、※橋脚 巨大物(自動車・がれき)の激突(津波・氾濫・増水)、床版ひび割れ(除雪車・大型車・地震)、主塔からの落雪、凍結、

落橋について

落橋とは、橋梁の劣化や大規模な災害などが要因となって橋梁が破壊し、文字通り橋が落ちることを意味する。

落橋の主な原因は橋梁の老朽化であり、損傷の種類は日本橋梁建設協会にある通りである。

また、大規模な災害、具体的には地震や津波、土砂災害、強風暴風、積雪も原因となり得る。

落橋がもたらす社会的影響は、周囲の交通機関に不便をもたらすこと、移動手段の減少から経済活動が著しく衰退すること、多額の資金がかかること(架替えは補修・補強よりもお金がかかる)が挙げられる。

ケーブル腐食について

腐食とは、化学や生物的な作用によって物体の外観や機能が損なわれる状態のことである。

金属の腐食を特に「錆」という。

腐食の主な要因は、塗装劣化個所に水分と腐食性物質が作用することによる。

腐食性物質としては、大気中の亜硫酸ガスと海塩粒子等が挙げられる。

これらは地域や気候に関係が深く、工業地域や海岸近くに位置する鋼橋において腐食の進行が早いのはこのためである。

発生のメカニズムとしては金属は大気中の酸素によって表面に参加被膜を形成している。

この酸素と水分、金属表面の汚染物などが表面に付着し酸化被膜が破壊されることがメカニズムである。

ここから腐食が表面から内部へと進行していく。

金属にとって「錆」状態は安定している。

腐食の種類は均一腐食と局部腐食の2つである。

均一腐食…金属表面全面に発生する腐食。

局部腐食…局部的な腐食で発見しにくい。一般的な製造部品などで問題になる錆による損傷はほとんどが局部腐食。

ここから本題である「ケーブル腐食」を考える。

斜張橋と吊り橋は両橋ともケーブルを用いる形状をしているが、その違いは“ケーブルと桁のつなぎ方”にある。

斜張橋は、主塔から張られているケーブルが直接桁橋を支えるが、吊り橋はケーブルと桁をハンガーが吊っている、すなわち実際に桁を支えているのはハンガーである。

また吊り橋はアンカーレッジが存在する。

斜張橋の強みとして、桁曲げモーメントを小さくできること・桁下空間を大きく取れること・設計自由度が高いことが挙げられる。

スパン長200~500mm及びそれを超えると斜張橋か吊り橋で設計することになるが、より長スパンに対応できるのは吊り橋である。

それぞれの橋梁の桁に加わる力として、

斜張橋…垂直方向の張力+橋軸方向の圧縮力

吊り橋…垂直方向の張力

がある。

ケーブルは通常「亜鉛メッキ銅線」で構成されているらしい。

亜鉛メッキとは高い耐食性(腐食のしにくさ)を付与出来るメッキ処理のことである。

亜鉛金属を電気の力で析出させる方法と溶融させた亜鉛に浸漬することで析出させる2つの方法がある。

亜鉛メッキは防錆効果が高く主に鉄製品に対して処理を行いますが、これは亜鉛メッキ上に不動態膜が生成されるからである。

また、亜鉛は鉄よりイオン化傾向が大きいのでメッキ皮膜にピンホール(小さな穴)があっても亜鉛が犠牲となって素地の錆を防ぎ高い防食効果を得ることが出来ます。

これを犠牲防食という。

ただし、デメリットとして亜鉛メッキは湿気に弱い。

ケーブルが腐食する主要因は「水(+塩)」であると考える。

特に、乾湿を繰り返す部位は、常に水に覆われた部位より腐食速度が速くなる。

その理由は、乾湿を繰り返す部位が常に水に覆われた部位と比べ水の層が薄く、腐食因子の一つである酸素の供給が限定されるためであり、かつ局部的な腐食電池を形成して局部腐食が進行するためだと考えられる。

ケーブル腐食がもたらす影響をいくつか考える。

①斜張橋の終局強度の低下…ケーブル断面の減少+伸び及び疲労強度の低下  から

②疲労寿命(破壊が起こるまで材料が耐えられる、ある種の変動応力や変動ひずみの繰返し回数)が小さくなる…ケーブル断面の減少  から

③橋梁全体の耐力低下…ケーブル断面の減少  から

④ケーブルの伸び及びねじれ強度の著しい低下…亜鉛メッキが腐食により消費され地鉄の腐食が進行した場合(亜鉛メッキのみの腐食だと特に影響なし)に、腐食に伴い表面凹凸が生じてその部分に応力集中が生じる  から

上記①~④の原因・理由を考えるために斜張橋ケーブルの耐久性評価のまとめから重要点を引用する。

ケーブルは主塔・主桁との協働作用により外力に抵抗する非常に重要な部材であるため、上記のようなことが起こると健全時のケーブル耐力を発揮できないことは自明であろう。

これにより腐食したケーブルで耐えている斜張橋や吊り橋は、日常に起こる何の変哲もない荷重(自動車・雪等)や揺れ(小規模地震、風等)で突然崩壊を起こす可能性が高い。

橋梁の突発的な崩壊は絶対に避けるべき事故である。特に斜張橋や吊り橋は、長スパンの橋梁であることが多いためなおさらだ。

従って、当たり前ではあるがケーブル破断が起こらないように点検・維持管理・補修補強は最重要事項である。

風によるケーブルの振動現象は

①カルマン渦による渦励振

②六角形断面ケーブルにおけるギャロッピング

③並列ケーブルにおけるウェイクギャロッピング

④レインバイブレーション

⑤塔からの剥離渦によるウェイクレゾナンス

リスクアセスメントについて

研究メモ

研究参考文献

春休み課題の進捗

     ローラー支持は「線固定」→線の端点を入力して、線を作成

     固定端は「面固定」→パーティションからグループの作成

     Mesh→ジオメトリは「Partition-1」を必ず選択

     Mesh切る→ジオメトリのグループ作成

   ☆☆AsterStudyは簡単にできる!!→Operations→Addstage with Assistant→Modelanalysis→情報入力  

   ☆☆(E:206000,ν:0.3,ρ:7.8e-9) 

   ☆☆ファイル名には「.med」!!

     builtinでopen→「Mode」に変更

     フィルター→「Mechanics」から「Normal modes animation(real)」ここからアニメーションや振動数を見れる

やること

      →でも、全ての結果を比較しても ①要素数と相対誤差との関係性が見られなかった ②誤差が小さくならなかった

            →よって「梁形態」が間違えていると考える

             ☆解決策 ①梁形態を変更(ピンを端っこから10mm内側に移動)

                 ②ピンを2コに変更(上下で線固定+面固定)→失敗、、多分ピンの固定条件が上手くいってないから。

                 ③ローラー固定で考えてみる→成功&正解!!

                    →!!メッシュ1.0以下は解析不可能であった。収束を見つけることが目標!!・・済

疑問・考察

→やはり、1次要素だとどれも水平モードの方が精度が良かった。

→鉛直・水平・ねじれを総合的に見ると長さ+幅が一番いいかも。すなわち、幅を小さくすること・長さを長くすることが与える影響が大きい。

春課題を通して、理論値への収束の条件や様子、理論値に対する相対誤差でのグラフの様子、収束値に対する相対誤差での0に近づく形態、振動モードと要素数との関係(これはいい結果は出られなかった)、細長比との関係や細長比が梁に与える影響を見ていった。

既に答えの出ている課題であると説明を受けていたが、どれも難しく感じた。

また、1つの知りたいことのために解析や考察を行ってもまた別の分からないことが出てきてなかなか答えにたどり着けないとも感じた。

ただ結果を出すだけでなく、そこから何が考えられるか、何がわからないのかを明確に進めていくことが大切だと学んだ。

3年創造工房

2023.10.27

・タッチタイプの練習

・コマンドの練習

pwd・・・今どこにいるのか

ls・・・pc内にファイルが何があるのか

mkdir・・・新ディレクトリ(フォルダ)の作成

cd・・・行きたいファイルを順々に降りていく

gedit〜&・・・新しいテキストを開く

vi・・・選んだファイルをプログラミング上に出し、書き換える

Escを押して:q・・・保存して戻る

cp・・・コピー

cat・・・指定したファイル名を見られる

rm・・・指定したファイルの削除

rmdir・・・ディレクトリの削除

cd 次行くところ / 自分のディレクトリ・・・2つ飛ばして自分のディレクトリに行ける

cd .. pwd・・・一つ前のディレクトリに戻る

〇コマンドの流れ

pwd→ls→(mkdir)→cd→ls→gedit~&(テキスト開く)→ls→(vi),(cp),(cat),(rm),(cd ..)

2023.11.10

・コマンドの練習

・gnuplotを使用して.pngで図を作成

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/bbb.png

使用データ:

メッシュ長さ要素数変位(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.445180.3986-4.34山本
1.5154330.396317756757-4.9進藤
1.6159000.399049-4.24河合
1.8116770.404457-2.9山口
2104600.394818715517-5.3進藤
37340.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河合

引用元:岡田の卒論日誌(11/18課題)

2023.11.17

・片持ち梁の解析(サロメ)

Geomety BoXを作成 適用→閉じる kotei/saika

Mesh メッシュ作成 Box1を挿入 アルゴリズム Note ジオメトリでグループを作成

FY=100N/100mm^2=1N/mm^2

・測定値=平均値ー理論値/理論値✕100(%)

・要素数=ボリューム

・変位平均値(1.3.5.7)

・相対誤差(測定値)を計算

・縦軸_変位、横軸_要素数のグラフを作成 

メッシュ長さ要素数先端変位(4隅の平均値)[mm]相対誤差(\( \frac{salome-手計算}{手計算} \))計算者
0.71984646.542811.91安藤
0.81138126.51042.39安藤
0.9402806.36315254.60兼田
1.1300556.33635255.00兼田
1.2264676.30433755.48柴田
1.3251806.3043555.48柴田
1.4322126.316125.31佐藤
1.5177536.12098.23佐藤
1.6142966.20446256.98皆川
1.7135966.21566256.81皆川
1.828665.73775513.98永山
1.960015.726362514.15永山
256175.645852515.355
323095.472875517.948
46173.61605750.458服部
54943.85803750.422服部
65812.5068262.416梶原
71331.4122578.827梶原
8781.288717580.68工藤
9721.287992580.69工藤
10601.1434482.85佐々木
11651.2312481.154佐々木

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/yousosuu_henni_rironti_sokuteiti.png

縦軸:変位(mm)、横軸:要素数

理論値:6.67mm

【理論式】

Pl^3/3EI

2023.11.24

・単純梁(等方性1次)の解析

メッシュ長さ要素数変位平均[mm]相対誤差(\( \frac{salome-手計算}{手計算} \))計算者
0.71452340.422484527351.388176安藤
0.81429730.422570445981.408794安藤
0.9916480.4204372865730.897兼田
1.1271600.4056189390242.659兼田
1.2246750.4043492.96柴田
1.3234460.4041853.00柴田
1.4177380.3986044.34佐藤
1.5154380.3965934.83佐藤
1.6161220.3982124.44皆川
1.7120260.3934115.59皆川
1.8116040.3936685.53永山
1.9103910.3906956.24永山
2109210.3951035.18
323280.32476222.06
415000.15501362.80服部
54320.06527884.33服部
63570.21306248.87梶原
71960.101975.55梶原
81040.115862472.20工藤
9810.125511869.88工藤
10780.0773381.44佐々木
11630.199952.03佐々木

https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/1129_guraf.png

縦軸:変位(mm)、横軸:要素数

理論値:0.4167mm

【理論式】

Pl^3/48EI

2023.11.29

・単純梁(異方性1次・等方性2次)の解析

<等方性1次と異方性1次の比較> [#he11b38c]

メッシュ長さ要素数先端変位相対誤差計算者
0.71445630.5052522.76安藤
0.81415170.5046922.64安藤
0.9916480.5025952.216兼田
1.1271600.4899140.363兼田
1.2246750.4870880.791柴田
1.3234460.48680100.995柴田
1.4177380.4859991.16佐藤
1.5154380.4851801.33佐藤
1.6159000.4832861.71皆川
1.7121420.4779522.80皆川
1.8116040.4820851.97永山
1.9103910.4708872.40永山
2102910.4809102.19
323280.43193712.15
415000.43015612.52服部
54320.28296842.45服部
63560.344155630.00梶原
71960.21393456.49梶原
81040.22987453.25工藤
9810.23230852.75工藤
10780.20327158.65佐々木
11630.22231654.78佐々木

https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/1208_guraf.png

縦軸:変位(mm)、横軸:要素数

理論値:0.4917mm

【理論式:ティモシェンコの理論】

v = Pl^3/48EI + Pl/4KGA

G:せん断弾性係数

K:せん断補正係数

このとき、[Pl/4KGA]= せん断項

<等方性1次と等方性2次の比較> [#of7eeadc]

メッシュ長さ要素数先端変位相対誤差計算者
0.71445630.4301243.22安藤
0.81415170.4301323.22安藤
0.9916480.4300203.197兼田
1.1271600.4298283.151兼田
1.2246750.4298363.15柴田
1.3234460.429743.13柴田
1.4177380.4297971.3佐藤
1.5154380.4299583.14佐藤
1.6159000.4297553.18皆川
1.7121420.4296763.11皆川
1.8116040.4298293.14永山
1.9103910.4296843.12永山
2102910.4296203.10
323280.4291692.99
415000.4292543.01服部
54320.4281702.75服部
63560.4284522.82梶原
71960.425912.21梶原
81040.4260742.25工藤
9810.4255522.12工藤
10780.48838217.20佐々木
11630.4239729.0534佐々木

https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/1208_guraf_2.png

縦軸:変位(mm)、横軸:要素数

理論値:0.4167mm

2次要素を含む。

・等方性・・方向によって物体の物理的性質が異ならないこと。

・異方性・・方向によって物体の物理的性質が異なること。

2023.12.08

・単純梁(木材+鋼材)の解析

メッシュの長さ要素数変位[mm]相対誤差計算者
0.71551920.0837890524615.365安藤
0.81388080.0838038649115.350安藤
0.9825870.08370707398115.45兼田
1.1386710.08420120760214.95兼田
1.2319290.08368815.466柴田
1.3286210.08366915.4857柴田
1.4288540.0836815.47佐藤
1.5200150.08405215.10佐藤
1.6194480.083540293815.62皆川
1.7138010.083435509815.72皆川
1.8125280.08373315.42永山
1.9117690.08392415.23永山
2106990.08407687655915.074
335790.0841456175315.004
416280.08279416.37服部
510160.08303318.89服部
68390.08288216.26梶原
75540.08087118.28梶原
82850.07999519.20工藤
92610.07898020.22工藤
102320.08191117.26佐々木
112080.07567623.56佐々木

https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/1214_guraf.png

縦軸:変位(mm)、横軸:要素数

理論値:0.099mm

2次要素を含む。

【理論式】

Pl^3/48EI + Pl/4KGA

このとき、

EI = EI(木) + EI(鋼)

G = G(木)

と考えることとする。

2023.12.15

・Taxツールの概要・練習

① tar xvzf sibuw.tar.gz

② cd sibuw

③ ls

④ vi sibup2.tex(=メインファイル)

ーー / がついているものがコマンド

ーー :q / :wq / iで編集可能

ーー :w 保存 ーー :! pdfplatex sibuw2.tex 元のファイルから編集点の変更を反映 

ーー :! ls / :! evince 7021531.pdf & ーー :! でコマンドの使用可能

⑤ pdfplatex sibuw2.tex

⑥ ls

⑦ evince subuw2.pdf &

ーーpdfを見るコマンド

ーー&を忘れない!!!

・texでは2回コンパイルする

・図・表を挿入した際には、必ず2回コンパイルする

<メッシュの変更> [#o655497e]

① Geometyを開く

② Meshを開いて、Partitionを選択してメッシュを作成→ここでメッシュの長さを変更(Mesh-2作成)

③ Mesh-2をクリックでメッシュ作成→要素数をメモ

④ Mesh-2を右クリックでジオメトリの作成

⑤ AsterStudyのmeshをダブルクリックして、UNITE「Mesh-2」に変更

⑥ Outputでファイルの名前変更(.med)

⑦ 解析

<テキストエディターでグラフ作成> [#k2d05487]

① テキストエディターに(要素数,変位)の順で入力→このとき、縦軸:変位・横軸:要素数

② 名前を付けてファイル保存

③ コマンド

ls→cd kaneta23→ls→gnuplot→>plot "ファイル名" with line,"ファイル名" with line・・・→>set term png→>set output "最終ファイル名.png"→>plot "ファイル名" with line,"ファイル名" with line・・・→quit→(eog 最終ファイル名.png)

④ メニューから、gftpに行きパスワードを入力

⑤ public htlsの中のj2023から、自分のファイルに転送


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