兼田の卒論日誌

研究†

研究テーマ†

• ケーブルを有する橋梁のリスクアセスメントの提案ー秋田県由利橋について

• 実橋において，その橋がどのような危険にさらされているか，ケーブルが破断する状況や可能性について，現地でも調査を行い，適宜，県の方々やコンサルの方との打ち合わせあり（ISS・神戸大との共同研究）
• リスクアセスメントは橋梁がどのくらいの割合で危険にさらされているか定量的に見極める方法で，海外ではかなり広く使われているが，日本ではあまり使われていない．時にケーブルを有する橋梁においてのリスクを提案するのはあたらしく，ケーブル小委員会でWGが立ち上がっている．由利橋のモデルをFEMで作成し，実際にはどのようなリスクでどのような崩壊過程になるか，これまでの連鎖崩壊との関係もある

やること†

• 参考文献を読む

• ☑由利橋のモデルをリニューアルする〜5/31

• ☑由利橋の事故・落橋の原因・落橋したときの不便・ケーブル腐食の理由を調べる・由利橋で起こり得る危険事項＋通れなくなった場合の迂回路や懸念点を考える

• ☑マイダスで３Ｄモデルを作成する…「ケーブル橋梁ウィザードをサポートしないバージョン」とエラー出る

• ☑プレテンションの設定方法を模索…単純梁✖簡易モデルでプレテンションあり・なしで比較

• ☑マイダスで由利橋のケーブルにプレテンション導入

• ☑由利本庄市の天気・交通事故を調べる〜5/31

• ☑リスクアセスメントのまとめ〜5/31

• ☑橋梁と火災について調べる

• ☑事故のときの運転者の心理について調べる

• ☑マイダスで火災や温度変化の解析ができるかどうかを試す

• ☑由利橋モデルに正しいプレテンションと死荷重・活荷重を入れる

• ☑死荷重及び活荷重をどう入れるか聞く（集中or等分布？、それぞれの計算方法）

• ☑主塔をソリッド✖シェルで作成

• 他の二径間斜張橋とたわみを比較して正誤性の評価

• ☑由利橋で事故が起こらない理由を考える

• ☑サロメで温度解析のやり方学ぶ

• ☑サロメで車高部分のみを加熱してどうなるか

• ☑鋼材に粘性を入れてぐにゃってなるか片持ち梁で確認→弾塑性でできた

• ケーブル２本の簡易モデルでブロック全体とブロック一部を加熱して、主塔の変形とケーブルの応力分布を見る

• [見通し]火災解析で鋼材（主塔・ケーブル）がどのように変化するかを調べる→運転者の心理や危険車の比率、事故の比率などと組み合わせて”リスクアセスメント”へつなげていく、リスクアセスメントの式を用いて

進捗†

由利橋モデルの作成(5/29)†

マイダスで由利橋のモデルを作成.

材料・断面・境界条件を設定、割当した.

スパン長190m,主塔70m（桁上50m,20m）

ケーブル断面積D=200mmのケーブル位置

ケーブル断面積D=225mmのケーブル位置

ケーブル断面積D=315mmのケーブル位置

マイダスで荷重解析(6/4)†

幅1m×高さ0.8m×スパン長10mの単純梁を作成し、集中荷重（節点６に荷重FZ=-10kN）と等分布荷重（節点１２～２２に荷重w=-1kN）を与え解析した.

全体図（(x,y,z)空間）

横から図（(y,z)平面）

上から図（(x,z)平面）

集中荷重（節点６に荷重FZ=-10kN）の写真

等分布荷重（節点１２～２２に荷重w=-1kN）の写真

以下が解析結果である、変位図・Ｓ図・Ｍ図である.

〇変位図

〇Ｓ図

〇Ｍ図

簡易モデルでプレテンション導入(6/4)†

幅10mm×高さ10mm×長さ1000mmの簡易モデルでプレテンション導入した・してないもので比較

簡易モデル

下記の写真の通りに荷重FY=-100Nを与えた

桁材料（SM400）について

ケーブル材料について

先端変位の結果は👇の通りである. なお先端変位はｙ軸下方向が＋である.

全ての結果において先端変位が－、すなわち先端が上に上がることはなかった.

硬すぎではないかと不安である…

！！結論：間違っていた

由利橋モデルにプレテンション導入(6/5)†

簡易モデルでプレテンション導入の確認（合っているかは別として）が終わったため、本番へ！

上写真がケーブルプレテンションなし、下写真がケーブルプレテンションありの結果である.

このとき集中荷重F=-1000kNである.

〇プレテンションなし

max変位=9.35mm

〇プレテンションあり

max変位=9.34mm

全てのケーブルに軸力1000kNを入れているにも関わらず、0.01mmしか変化がない.

おかしいのでは…？？

ちゃんとプレテンションが設定されていないのかもしれない.

！！結論：間違っていた

簡易モデルでプレテンション導入＿改訂版(6/7)†

千代岡さんにアドバイスをいただき、軸力ではなく「温度荷重」で与えてみたところいい感じになった.

簡易モデルは前回のと同じものを利用.

温度荷重をケーブルに与えた.

簡易モデル

〇ケーブル材料

→弾性係数:1.95e+02

→ポアソン比:0.3

→線膨張係数:1.2e-05

→単位体積重量:7.7e-08

〇ケーブル断面積

→D=1mm

〇桁材料

→SM400

〇桁断面積

→10mm×10mm

解析を回してみたところ👇このように変形した.

軸力あり・なしバージョンで比較したが、軸力は関係なさそう…

また、先端荷重Fを与えたときの先端変位は下記の表のようになった.

このとき、ｙ軸下方向を＋とする.

表からわかる通り、荷重ありとなしのときの変異差はどの温度の時も等しく、

F=-100N…1.835mm

F=-500N…9.174mm

であった.

この結果から、プレテンションは温度荷重から与えることができ、当たり前に先端荷重を加えると下へ下がることがわかった.

加えて、温度荷重が相対的に大きいほど先端変位も大きくなることがわかった.

温度荷重と先端変位は比例関係にあると言える.

由利橋モデルにプレテンション導入＿改訂版(6/7～6/10)†

軸力ではなく、温度荷重からプレテンションを入れる方法を由利橋モデルでも実践した.

このとき全く軸力は入れていない.

下記の写真通りに、外側のケーブル２本に-150℃・その他のケーブル１０本に-100℃を設定.

①なんの荷重もかけずに解析を回すと👇このように変形した.

max変位＝-0.394mm

②次に中央（節点９６）に集中荷重F=-190kNを与えた場合（下記の写真）

結果は👇このように変形した.

max変位＝+1.852mm

③次に全桁にw=-1kNを与えた場合（下記の写真）

結果は👇このように変形した.

max変位＝+1.008mm

この結果から、プレテンションをのみの場合は桁が引っ張られて上に上がることがわかり、集中荷重及び等分布荷重を与えた場合は桁が下側に押される様子が見て取れる.

また、等分布荷重よりも集中荷重の方が変位が大きくなった.

きちんとプレテンションが入っていると思われる.

主塔に熱解析(6/13)†

主塔のみを取り出して100℃,500℃,1000℃を与えた.

主塔の初期温度は20℃にしてある.

今回主塔は、桁下の20mも取り出してあるため全体は70mである.

また、変位の単位はcmにしてある.

〇主塔材料

→ヤング率E=206000 N/mm^2

→線膨張係数α=1.17e-05 /℃

→断面二次モーメントＩ=8.4375e+11 mm^4

〇主塔断面積

→3.0m×1.5m

初期20℃→解析後120℃

初期20℃→解析後520℃

初期20℃→解析後1020℃

結果からわかる通り、加える熱が高ければ高いほど変位は大きくなる.

＋100℃のときと＋1000℃のときの変位はちょうど１０倍になっている.

本来であれば400℃～500℃で主塔は耐えられなくなるが、この解析ではそれが再現できない.

この解決策として「弾塑性」を設定に加えなければならないらしい.

プレテンション荷重の入れ方(6/14～6/17)†

今回はプレテンションを温度荷重ではなく。プレテンション荷重項目から与えた.

手順は以下の通りである.

荷重＞プレテンション荷重＞荷重ケース選択＞kNを入力＞単一選択＞ケーブル要素をクリック＞適用　　これをすべてに行う.

プレテンション荷重以外には何も与えずに解析を回したところ桁が持ち上がった.

max変位＝-73.49mm

ケーブル（トラス要素）の引張（断面力）はこのようになった.

ケーブル（トラス要素）の応力度はこのようになった.

次に桁を見ていく.

桁（梁要素）の断面力はこのようになった.

これは断面図である.

桁（梁要素）の応力度はこのようになった.

これでケーブルプレテンションは入っているのではないか・・？

これが合っているのであれば、温度荷重への変換は必要なし！

！！結論：プレテンションは「プレテンション荷重」の項目から入れてよい.

プレテンションを温度荷重に変換する方法(6/14～6/17)†

線膨張係数αの計算式である👇これを使う.

α = \( \frac{ΔL}{L}\frac{1}{ΔT} \)

ΔL/L…ひずみ

ΔL…変化量

L…元の長さ

今日は、先にプレテンションを設定して解析して、要素作成タブから元のケーブル長さと解析後のケーブル長さから差を求め、計算式に代入してΔTを出すという手法を取った.

すると、Ｃ１ケーブルに与える想定温度1250℃、Ｃ２ケーブルに与える想定温度1762℃となった.

大きすぎな気がする…

！！結論：プレテンションを温度荷重に変換する必要はない.

活荷重の入れ方(6/17〜6/19)†

活荷重の計算は桁の半分を対象にしている.

①　p1 = 10kN/m^3 × 5.5m/2 × 1.1m = 30.25kN/m

②　p2 = 3kN/m^3 × 5.5m/2 × 1.1m = 9.075kN/m

③　p1/2 = 5kN/m^3 × 6.75m × 1.1m = 37.125kN/m

④　p2/2 = 1.5kN/m^3 × 6.75m × 1.1m = 11.1375kN/m

桁0m～62mあたりには、②＋④(20.2125kN/m)のＬ荷重

桁63m～72mあたりには、①＋②＋③＋④(87.5875kN/m)のＬ荷重

桁73m～190mあたりには、②＋④(20.2125kN/m)のＬ荷重

を付与.

下👇の写真の通り.

活荷重のみを与えた変形は👇このようになった.

max変位　＝　+35.36mm

死荷重の入れ方(6/17)†

死荷重の計算は桁の半分を対象にしている.

D = 77.0kN/m^3 × (9.5m+6.6m) × 1.1m /2 = 681.835kN/m

この死荷重を等分布荷重として桁に与えた.

まだカウンターウェイト等は考えずにすべてを鋼材として計算した.

死荷重のみを与えた変形は👇このようになった.

max変位＝+955.91mm

由利橋モデルにプレテンション＋死荷重＋活荷重を組み合わせたときの様子(6/17〜6/19)†

上記のプレテンション・死荷重・活荷重を全て由利橋モデルに与えた場合には👇このように変形した.

max変位　＝　+917.82mm

プレテンション・死荷重・活荷重を足し引きしたものとほぼ同値.

デフォルトで約1mもたわむものなのかと少し不安…

主塔３Ｄモデル作成_midas(6/20～6/26)†

マイダスで主塔の３Ｄモデルを作成した.

やり方は　ウィザード＞基本構造＞シェル＞数値入力　からできた.

手順は写真の通りである.

ここでいう、

ｎ：ｙ軸方向の分割数

ｍ：ｘ軸方向の分割数

Ｉ：ｚ軸方向の分割数

である.

今回は、ｎ＝10コ,ｍ＝10コ,Ｉ＝50コ　にした.なお適当.

設定した後は👇このようになった

全要素に400℃を与えて解析したところDX,DY,DZ方向全てに変位が現れた.

なお、まだ塑性や降伏点は設定していない.

〇DX方向変位

max変位＝7.62mm

〇DY方向変位

max変位＝4.20mm

〇DZ方向変位

max変位＝240.64mm

〇DXYZ方向変位

max変位＝240.77mm

一つ問題点…境界条件を節点を完全固定にしているため、３Ｄにすると要素の面が引っ張られてしまう

失敗したやり方…節点を作成して要素で結んで、後から断面や材料を設定すると下の写真のようにBOXにはなるが、実際にはビーム要素のため３Ｄにはなっていない

塑性材料を設定しようとすると「静的材料非線形関連機能をサポートしないバージョンです。Web認証の場合、オプションのチェック有無を確認してください。」とエラー出る.

1/100主塔モデルで温度解析_Salome(7/3〜7/10)†

1/100モデルにした理由は、Asterstudyの設定以外でのエラーを出さないため.

7/3：合っているかはわからないがある程度Asterstudyの設定が完了した.

残り、✅熱膨張係数（ELAS-ALPHA）・✅初期温度（AFFE_VARC-TOUT-NOM_VARCをTEMPと入力-CAMP_GDをtempsと選択-VALE_REF）・EnforceDOC（訳：DOCを強制する）？？・✅Externalstatevariableassignement（訳：外部変数割り当て）（CAMP_GDをtempsと選択）←tempSを反映させるため・✅モデル全体に温度与える設定（AFFE>TOUT）

を行えば良い.やり方わからなかった.聞かなければ、、

温度解析についてはSalomeメモの熱関係とSalome-Meca_熱応力メモを参考にした.

7/4：なんとか設定一通り終えて解析をした.成功もした.

結果を見た感じ変なところはなさそうだが、マイダスとの結果が大きく異なる.

これは、マイダスはシェル（厚さ10mmの中は空洞）で、サロメはソリッド要素（全て鋼材）であることが関係しているのか？

確認が必要かもしれない.

↑この理由は1/100モデルだからか、、忘れてた。。

結果

今回はマイダスと比較するために、先端部分の大体マイダスと見たところが等しいところのDX,DY,DZの変位を抜粋した.

上記写真の上端である.

なお単位はmmである.かつサロメは結果の数値を意図的に100倍したものである.

算出した、EPSI(ひずみ),DEPL(変位),SIEQ(主応力・ミーゼス応力等),SIGM(6成分応力),をグラフにしてみた.

　　　

7/5:熱解析のやり方をまとめ、上記のグラフを出してみた.また、100℃に温めた際にどう変化するのかを解析した.

結果的には、100℃のときの最大変位はDX=0.023mm,DY=0.014mm,DZ=0.60mmくらいになった.

この値は400℃のときの最大変位の約1/4であることがわかった.

400℃のとき、DX=0.094mm,DY=0.058mm,DZ=2.41mm

このことから熱解析自体は正しくできているだろう.

次はいよいよ弾塑性を入れてみるか、、

7/9:弾塑性のやり方を学んだのと、実際の大きさで車高部分に400℃を与えた.

モデルは、下から0.5m,0.5m〜3m,3m〜50mの３分割にした.

車高は1.5mとし、そこから下に1mかつ上に1.5mの部分が真ん中のBOXである.

このモデルの真ん中のBOXにのみ400℃与えた.

結果はこのようになり、加熱部が一番大きく変形した.

Pravisの写真の点部分.

３つのBOXの手前左上.

弾塑性入れてやってみたが破壊形態は見れなかった？見方がよくわからなかった？

今は瞬間的に0℃→400℃に変えて解析をしているが、現実性を考えて1秒毎に何度与えていくという解析にした方が良いのではないかと考えている.

！補足！

マイダスと比較するために、サロメでも実モデル全体に熱を加えた場合の結果は

熱解析自体は上手くいっていると思われる.

7/10:ゼミ

ゼミでただの鋼材に熱を加えるだけでは膨張するだけだと言われた.

そこで粘性係数を考慮する必要がある.

熱解析で鋼材が曲がるか確認するためにまずは簡単な片持ち梁等で解析してみるのは良さそう.

改めて、リスクアセスメントとは(7/17)†

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

粘性を入れた片持ち梁の解析(7/19〜7/24)†

まず粘性の入れ方を学ぶために10×10010 まず粘性の入れ方を学ぶために10×100×10の片持ち梁で行った

CodeAster p137〜140に書いてある

ヒルの降伏応力

円筒座標系

直交座標系

テンソルとは

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

CodeAsterを解読した

SalomeのLNTは下のように方向を表す

温度における、ac1＝準静的温度開始、ac3＝準静定温度終了

これはAstersutudyの

DEFI_MTERIAU>META_LEMA_ANI(粘性)の設定画面である

　　　　

CodeAsterより等方性の場合はこのようになる👇

γ、θ、Z方向は円筒座標系であり、x、y、z方向は直交（デカルト）座標系である

この２つの座標系の関係性は

[x=γcosθ]

[y=γsinθ]

[z=z]

となる

①熱

②熱＋弾塑性

③熱＋粘性

④熱＋弾塑性＋粘性

①のPravis　　

②のPravis　　

①と③の結果は同じになり、②と④の結果は同じになった、つまり粘性による影響がないことになる

鋼材の動粘性を"0.5mm^2/s"として入れた

　　　　

[方向]

F_M〇〇_〇〇=C_M〇〇_〇〇=1.0,F_M○×_○×=C_M○×_○×=0.75

[粘性に関して]

F_N=C_N=流れの粘塑性

F_C=C_C=粘性加工硬化の回復

F_M=C_M=粘性開始の加工硬化の回復

F_A=C_A=？？（Aはalphaのこと）

F_Q=C_Q=？？（ケルビン温度表記）

簡易主塔モデルに加熱(7/25〜)†

・ソリッド要素 ・加熱する時間の設定１０〜３０分 ・ケーブルの応力分布 ・時間ごとに熱くなるのを再現できるか

7/25:ソリッド要素断念

主塔とC6・C7ケーブルをソリッド要素で作成を試みた

しかし、下の写真のように主塔とケーブルのつなぎ目部分が上手く噛み合わず（平行にならず）に断念した

移動や回転も試してみたが、思った方向に変化しなかった

ビーム要素に変更する

7/26:ビームモデルの作成

code-asterによる解析 step-by-step 構造要素の概要

ビーム要素でモデルを作成した

ケーブルと主塔の間に小さいバネが必要であるため、以下のように設定した

今回のモデルでは10mmの長さにした

主塔やケーブルには架空の断面を設定する必要がある

主塔：Geometryの長方形フェース　　ケーブル：Geometryの円盤を作成

写真ではケーブル断面を実際のケーブル位置に作成したが、そこに作らなくても良いらしい

原点に作っておいてAsterstudyで割当できるらしい

ビーム要素のも出るを作成するに当たっての注意がある

①架空の断面はGeometryでグループ作成してから、Meshでジオメトリのグループ作成が必須（いつもの右クリックでのやつ）

②Meshのコンパウンド必須、これはビーム要素のみで断面はしなくて良い

③②でしたCompound_Meshからグループを作成で「ノード」ー固定（主塔・ケーブル２つ）と「Edge」ー要素まとめ（主塔・ケーブル２つ・バネ２つ）を作らなければならない。このとき、主塔を分割して設定するとAsterstudyで加熱部を選択できるようになる

○モデルのSalome設定一覧

　　　　

由利橋について†

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

橋長190.5mの不等径間の２径間連続斜張橋（鋼床版２主箱桁）である。

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

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

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

今までで由利橋で一度も事故が起こったことはない！！→なぜ？→下記に記載

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

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

考え得る危険事項…交通事故（追突・接触）、☆火災（自動車事故）、揺れ（地震・風・共振）、落雷、積雪、津波、※橋脚　巨大物（自動車・がれき）の激突（津波・氾濫・増水）、床版ひび割れ（除雪車・大型車・地震・塩化カルシウム）、主塔（水平材なし）からの落雪、凍結、鋼部材の錆による劣化（雨や雪の耐水）

使用不可になった場合の懸念点…移動手段の減少による渋滞・移動時間の増大、経済活動への影響、交通の混乱につながる（引き返したり止まったりなどなど）

※このときの迂回路…飛鳥大橋、本荘大橋、…子吉川が流れているため迂回路は他の橋しかなさそう

令和4年：交通事故統計

令和4年：秋田県交通事故発生状況

国土気象庁

国土気象庁日ごと

本荘の雨温図

2022年：由利本荘市気象

由利本荘市年間気象データ　※データ古い

[由利橋で事故が起こらない理由]

・単純にスパン長が短い

・走行車両の速度が比較的遅い

・ほとんどが軽自動車か普通自動車、たまに軽トラックであったから大型車や危険車両が通ることが少ない

・渡り終わったすぐ先が信号機であるため慎重な運転となるため

・見通しが良い

・歩行路と道路がフェンスで区切られているため、歩行者との接触がしにくい

落橋について†

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

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

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

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

ケーブル腐食について†

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

吊り橋…垂直方向の張力

がある。

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

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

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

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

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

これを犠牲防食という。

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

ケーブルが腐食する主要因は「水（＋塩）」であると考える。

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

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

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

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

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

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

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

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

• ケーブル断面の欠損により、ケーブル鋼線が降伏するまでは斜張橋の各部材に与える影響は限定的であった。よって、ケーブルが降伏するまでは変化量が小さく目視による変状の確認は難しい。

• 小段数ケーブル斜張橋の場合は、ケーブル破断の寸前になると主桁の一部が降伏した。よって、広範囲の部材が損傷すると予想される。

• マルチケーブル斜張橋の場合は、ケーブル１本が破断しても特に影響はない。⇔小段数ケーブル…１本破断してだけで致命的な影響あり、すなわちより密な維持管理が必要。需要なのは「応力腐食割れ」を防ぐことである。

• 応力腐食割れ（腐食環境下において金属材料に引張応力が作用することで材料に割れが生じる現象）の発生条件は明確ではないため、とにかく長期的な滞水を起こさないことが第一である。

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

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

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

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

• 補足

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

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

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

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

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

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

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

[今後の目標]

…人命リスク・経済リスク・事故リスク・その他リスクを数値的に表すか、定量的に表すか、どうするかを決める

…他にどのようなリスクが考えられるか意見を出し合う

…何分の何の確率でどのようなことが起こるかを求める？

　（例）〇〇が△/□の確率で起こり、どのケーブルがどんな条件で切れるか

　（例）交通事故すなわち衝撃でケーブルや橋脚などにどう影響するか

　（例）ケーブルのたわみが今後どんなときにどう影響して何が起こるか

…リスクアセスメントをどのように紹介して、どう提案していくかを固める

…何解析（振動・衝撃・座屈等）をするのかを決める

[MY IDEA]

・健全時のケーブルと腐食時のケーブルで、振動の仕方・同荷重をかけたときの変位や壊れるまでの時間を比較して割合等を出す。その結果からリスクがこれくらい増えると言って「リスクアセスメント」の大切さを紹介する

→リスクアセスメントのやり方や実用性を紹介するのではなく、あくまで”大切さ”に焦点を置くためリスクアセスメントに対する根拠が薄くなる可能性大 　（例）健全時・腐食時のケーブルを比較したところでリスクアセスメントに繋げられるのですか？

・〇〇/△の確率で起こる小災害・中災害・大災害をそれぞれ由利橋にぶつけてどうなるか解析する。その結果からここまでは耐えられるけど、ここからは耐えられないという線引きをして、耐えられないリスクを考慮していない場合にどのような事態になるのかを紹介

→災害規模の判定が個人になるため定義が曖昧になることが問題。耐えられない災害の確率によって、起こるか起こらないかわからない災害を考慮することがはたして現実的か（労力やお金を考えて）という問題 　（例）5000年に一度の災害を「今」考慮すべき理由は？根拠は？それがリスクアセスメントだとしてもちゃんと説明しないとふわふわしてしまうかも・・

・健全由利橋に振動・荷重・衝撃を加えてどうなるかを見る。その結果何が起こる影響があるのか考える。

→やることが多くなるため一つにかけられる時間が限られる。また、この結果からリスクアセスメントの確率を出すにしても母数が少なくて信憑性が低い。本当にその確率は正しいのか、その根拠は？が肝になると思う

[文献から…]

リスクアセスメントとは、危険が発生する可能性の度合いを組み合わせて、リスクを見積もり、そのリスクの大きさに基づいて対策の優先度を決めた上でリスクの除去又は低減の措置を検討することをいう。

橋梁におけるリスクアセスメントは、橋梁がどのくらいの割合で危険にさらされているか定量的に見極める方法のことをいう。

欧米ではリスクアセスメントが主流なため「何が現実に実行できるか」で安全を考えている。

日本ではこの考えはまだまだ定着しておらず、リスクゼロを前提とし設計を行っている。

しかし実際にはリスクゼロは不可能であるため、日本でもリスクアセスメントが定着するように取り組むことが目標である。

[リスクマネジメント&リスクアナリシスの流れ]

基本条件・情報の整理→リスクの定量化・ランク付け→対策を実施するリスクの決定

[ALARP：As Low as Reasonably Practicableの原則]

① 「受け入れ可能なリスク」として「自然災害の発生レベル以下」というガイドラインを設ける．

② リスク低減のためにどこまでも費用をかけることが合理的ではない．リスクは実行可能な限りできりだけ低くする.

③ 低減措置の実施基準は，リスク低減効果と低減に要する費用とのバランスにおいて，低減に要する費用が著しく不相当とならない基準で判断する.

また安全に対する考え方で、，「広く受容される領域」・「我慢できる領域（ALARPの領域）」・「受容されない領域」の3つの領域で安全を考えている。

ALARP領域はそのリスクが許容されるわけではないが、ある便益を確保するためにはリスクが適切に制御されているという条件の下で、そのリスクを伴う生活を受容する（我慢する）領域である。

許容可能なレベルに達していないリスクから優先的に対策を実施することが最も重要である。

事故の発生件数に重きを置くのではなく、重大事故の防止を徹底する。

リスクアセスメントのポイントは「事前にどれだけのリスクを想定しリスク低減のための対策が効果との関係でどのように検討されたか」である。

日本では想定したリスクにすべて対策を講じるという考えが強いが、リスクゼロとするために想定されるリスクにすべて対応するには莫大な費用を必要とするため現実的には厳しい。

この「すべて対策を講じる」という考え方が弊害となって日本にリスクアセスメントが浸透しにくいのではないだろうか。

リスクアセスメントで大切な要素は以下の２つであると考える。

①我慢できる領域と受け入れられない領域の線引きをしっかりとすること

②受け入れ不可能な領域のリスクに対しては優先的に対策を実施すること

すなわち、これは諦める・これは絶対に阻止するという取捨選択が必要であると思う。

↑この考え方は日本人には受け入れずらいため、なかなかリスクアセスメントが進まないのかもしれない。あとは単純にお金問題だろう。

リスクアセスメントは重大なリスクに対処するための合理的な措置とされており、危険を考慮して設計することが大切だ。

また、過去に起こった良い経験も悪い経験も同等にフィードバックすることが必要不可欠となる。

具体的には設計の初期段階で、

＝通常の使用＋事故＋誤使用＋早期劣化＋etc..

ただ、懸念すべきことはリスクアセスメントを考えれば考えるほど、お金と時間がかかるということだ。

起きるかどうかわからないことにお金や時間をかけることは、企業や世間は難色を示すと考える。

これもリスクアセスメントの実現が難しい理由に当たると思う。

IHIの方の講座(5/31)†

現代では建設において、ＣＯ２の削減が優先である.

ＣＯ２の排出量はアメリカは横ばい、ＥＵは下がっている.

今の日本では【１橋梁/104人】を支えている.

土木の現場では、ＣＯ２排出量は　コンクリート＞鉄筋＞桁　の順で多い.

ＣＯ２を減らす具体例としては、耐候性を用いること・木材を使用すること・強い鋼材を使用すること・コンクリートを減らすことなどがある.

ヨーロッパでは、工事において重要視されているのは費用よりも「ＣＯ２」や「騒音」だそうだ.

しかし日本は人口が多いため、現時点ではＣＯ２を考えることは難しいらしい.

今の時代に沿った新しい優先すべきことが出てきたと考える.（今回はＣＯ２の話）

橋梁と火災の関係†

鋼橋にまつわる火災事故は、国内外問わず数多く報告されている。

しかし実際に発生することは稀であるため、その対策をしていない橋梁がほとんどである。

よって、いざ火災事故が発生してしまった場合に、利用者の安全を考えて橋の供用を停止すること・利用者の利便性のためにいかに速やかに供用を再開することの２点を苦慮することが多いのだろう。

そこで、土木学会鋼構造委員会では，2012年6月に「火災を受けた鋼橋の診断補修技術に関する研究小委員会」を設立し、ガイドラインを作成した。

ここでは、火災による熱の影響を受けた材料の特性はどのように変化するのか・供用再開の判断は何をどのように点検及び調査すれば良いのか・火災による損傷はどのように補修補強すれば良いのかといった対応マニュアルに当たる。

供用中及び工事中に考えられる火災原因を列挙する。

供用中…不審火・走行車両（タンクローリー等）の横転及び炎上・事故によって損傷した燃料やガス管からの流出が原因の爆発

工事中…動力工事・溶断及び溶接による火の粉の飛散・塗料作業中のシンナー（シンナーは橋梁材料に当たる）・弾性シール材やバックアップ材への引火・コンクリート型枠剥離剤・ガス器具・溶接機・

材料自体が出火原因となる事例はほとんどない。

また。延焼に必要不可欠なものは①火源、②燃えるもの、③酸素の３つである。

延焼要因は、放射熱・接炎・熱気流・火の粉などがある。

加えて炎が風下側に傾き隣接した建物に直接接炎する場合もある。

鋼橋は火災の被害を受けると、構造部材は温度上昇に伴う強度低下により体力が低下するとされる。

一般的に鋼材（構造用鋼材）は、受熱時の温度が600℃を超える場合、高力ボルトについては450℃を超える場合に冷却後の引張強度が低下するとされている。

鋼材は合熱処理によって結晶構造が変化するという性質を有しており、受熱温度によってはもともと付与された性質が失われる懸念がある。

よって、火災による鋼材の損傷評価については受熱温度の推定が重要になる。

このときの鋼橋の補修ですべきことは、塗装の補修・鋼桁の補強・部材の取り換えなどである。

火災発生から補修施工までの流れ👇

引用元：火災を受けた橋梁の安全性評価と補修＿実例

主塔は400℃、ケーブルは700℃くらいまで耐えられるそう。

SS材は通常〜300℃くらいまで耐えられる。

SM材は400℃〜500℃程耐えられる。このとき、500℃まで耐えられるのは鋼種による。

火災による影響はケーブルよりも主塔にある。

そのため主塔が受熱する温度や時間によって、耐力や強度、変形がどのようになるのかを調べ、そこからケーブルに与える影響を見ていくことが重要。

例）主塔が崩壊したときのケーブルの切れ方・主塔の耐力が落ちたときのケーブルへの力の加わり方・・etc

健全度評価フロー👇

鋼材部の被災目安👇

引用元：火災を受けた橋梁の健全性評価に関する試験調査

鋼材の温度が高くなると、柔らかくなるすなわちヤング率は小さくなる。

火災による鋼材の変質は目視試験での確認は困難である。

また、実鋼材を抜き取り機械試験などにより確認することは通常困難なため、鋼材の健全性は部材の受熱温度を推定し間接的に評価することが一般的である。

火災の被害を受けた鋼部材における補修を行うべき損傷は主に入熱による変形及び塗装の剥離であると考えられる。

車両が燃える原因として、燃料やオイルなど着火物が漏れたり、可燃物が接触したりすることで高温になる排気系で引火火災が発生するケースなどがある。

運転者の心理について†

「人」が原因となって起こる交通事故＝”ヒューマンエラー”

エラーの種類として、思い込み・知識や技能不足・意図的ルール違反・妨害による・緊張状態が招くものが存在する.

また、急いでいるときには必然的に事故を起こしやすい.

疲労・ストレスが交通事故へつながることも多々.

事故を起こした直後の運転者は二次被害につなげないように、路肩へ車を止めるであろう.

その行動が橋梁上では、主塔やケーブルに近づくことを意味する.

もし、追突や衝突で車両が爆発した上でこの行動を取ったとしたら、主塔やケーブルへ火が燃え移ることも容易に想像できる.

交通事故の原因について†

令和４年における交通事故の原因ランキングを交通事故の原因ランキングから抜粋

１位　安全不確認

２位　脇見運転

３位　動静不注意（他の車や歩行者の存在には気づいているものの危険ではないと判断し、それらの動静への注意を怠ったこと）

４位　漠然運転（ぼんやりと運転をしている様）

５位　交差点安全進行

６位　運転操作不適（危険は認識しているものの、それに対する措置を誤ってしまったことによる運転操作全般のミスを指す）

７位　一時不停止

８位　信号無視

９位　歩行者妨害等（横断歩道を渡っている、または渡ろうとしている歩行者がいるにも関わらず、走行車が一時停止をしなかったこと）

１０位　優先通行妨害（信号機がない交差点で、標識や道路幅により優先道路とされる道路を走行する車両の進行妨害をした際の違反のこと）

近年は交通事故数が減少している.

これの主な要因を内閣府は以下のように示している.

１，道路交通環境の整備

２，交通安全思想の普及徹底

３，安全運転の確保、車両の安全性の確保

４，道路交通秩序の維持

５，救助・救急体制等の整備等

６，交通安全基本計画に基づく諸対策を総合的に推進してきたこと

単路や橋での事故として多いものは様々な交通事故統計にある通りである.

１，車両相互追突

２，車両相互正面衝突

３，車両相互その他

４，車両単独工作物衝突

サロメモ†

AsterStudyの設定翻訳

[メッシュ]

LIRE_MAILLAGE(メッシュを読む)>FORNAT()

ASSE_MAILLAGE(メッシュ組み立て)

DEFI_GROUP()

MODI_MAILLAGE(メッシュ方法)

CREA_MAILLAGE(メッシュ決め)

[Model Definition]

AFFE_MODELE(モデル割当)>AFFE(割当)>TOUT(すべて)>PHENOMENON:MECANIQUE(機械)>3D

AFFE_CARA_ELEM()

DEFI_BASE_REDUITE()

DEFI_DOMAINE_REDUIT()

DEFI_GEOM_FIBRE()

MACR_CARA_POUTRE()

MODI_MODELE()

[Material]

AFFE_MATERIAU(材料割当)>MODELE>AFFE>TOUT(すべて) ←AFFEでTOUTとあるがこれは対象構造物全体を表すor GROUP_MA(グループ選択)

AFFE_MATERIAU>AFFE_VARC()>NOM_VARC()→TEMP(温度)>VALE_REF(初期温度)>TOUT(すべて)→OUI(はい)>CHAM_GD()>CREA_CHAMP選択

CREA_LIB_MFRONT

DEFI_COMPOSITE

DEFI_GLRC>

DEFI_MATER_GC

DEFI_MATERIAU(材料特性)>ELAS(ヤング率等)>E(ヤング率),NU(ポアソン比),ALPHA(線膨張係数)

DEFI_MATERIAU>TRACTION(摩擦)>SIGM:point

DEFI_PROP_ALEA

DEFI_TRC

INCLUDE_MATERIAU

[BC and Load]

AFFE_CHAR_CINE

AFFE_CHAR_MECA()>DDL_IMPO(固定)>GROUP_MA　　　　　　

AFFE_CHAR_MECA()>DDL_IMPO(固定)>FORCE_ARETE(線載荷)orFORCE_FACE(面載荷)

AFFE_CHAR_THER

CALC_CHAR_SEISME

DEFI_CABLE_BP

DEFI_CONTACT

DEFI_FLUI_STRU

DEFI_FONC_ELEC

DEFI_FONC_FLUI

DEFI_OBSTACLE

[Function and Lists]

DEFI_FONCTION(変化の度合い)>NON_PARA>INST

DEFI_LIST_REEL(解析の間隔の設定)>DEBUT(開始時間)

DEFI_LIST_INST(解析のオプション)

[Analysis]

Changing with other codes

Dynamics(力学)

Statics(静力学)>MECA_STATIQUE(機械静的)>CHAM_MATER(材料選択？)

STAT_NON_LINE(非線形)

COMPORTEMENT(計算モデル、適用する式の設定)>RELATION(構成式)>VMIS_ISOT_TRAC（Von Mises降伏条件、等方硬化則、多直線近似>DEFORMATION(座標変換モデル)>SIMO_MIEHE(大変形解析)

NEWTON(ニュートン法の設定)>MATRICE:TANGENTE(接線剛性)

CONVERGENCE(収束判定の設定)>RESI_GLOB_RELA(相対誤差)>ITER_GLOB_MAXI(最大反復数)

Thermics(熱力学)

THER_LINEAIRE(温度＿線形)

THER_NON_LINE(非線形)

THER_NON_LINE_MO(非線形モデル？モード？)

[Post Processing]

CALC_CHAMP(調べたい分野)>CONTRAINTE(対照？)>SIGM(6成分応力),SIEF(),SIPM(),SIRO()　　

CALC_CHAMP(調べたい分野)>DEFORMATION(力？)>DEGE(),EPME(),EPSG(),EPSI(ひずみ),EPSL(),EPVC()

CALC_CHAMP(調べたい分野)>CRITERES(基準)>EPEQ(),EPGQ(),EPMQ(),SIEQ(主応力・ミーゼス応力等)

NOEUは節点(NOEUd(フランス語))、ELGAはガウスの積分点(ELement GAuss point)、ELNOは要素(のどこか?)(ELement Node)だとか、、

RESULTAT>MECA_STATIQUE

CREA_CHAMP>MODEL→AFFE_MODEL>TYPE_CHAM→NOEU_TEMP_R>OPERATION(操作)AFFE>TOUT>VALE(温度変化の度合い)>NOM_CMP()→TEMP入力

[Output]

IMPR_RESU(結果照会)>RESU(結果)>RESULTAT(結果)>NOM_CHAM(ノーマル分野)

Salome-Meca2020_熱解析†

○Geometry

いつも通り作成

○Mesh

いつも通り作成

分割したソリッド要素をメッシュでひとつのソリッドにすることができる.

「コンパウンドを作成」から全てのオブジェクトを選択して適用

○AsterStudy

右クリックAddstage

メッシュ

• LIRE_MALLAGE

Mesh_1の選択>FORMAT>MED

Model Definition

• AFFE_MODELE

AFFE>TOUT:OUI>PHENOMENE:MECANIQUE>MODELISATION>3D>

MAILLAGE:mesh(LIRE_MAILLAGE)

Material

• DFFI_MATERIAU

名前:steel>ELAS>E,NU,ALPHA

• AFFE_MATERIAU

MODELE>model(AFFE_MODELE)>

AFFE>TOUT:OUI>MATER>steel(DFEI_MATERIAU)>

AFFE_VARC>NOM_VARC:TEMP>VALE_REF:0.0>TOUT:OUI>CHAM_GD:temps(CREA_CHAMP)

DEFIやってからAFFEの順でないとだめ！

BC and Load

• AFFE_CHAR_MECA

DDL_IMPO>GROUP_MA>DX,DY,DZ

Analysis

• MECA_STATIQUE

CHAM_MATER:fieldmat(AFFE_MATERIAU)>

MODELE:model(AFFE_MODELE)>

EXCIT>CHARGE:load(AFFE_CHAR_MECA)

Post Processing

• CALC_CHAMP

CONTRAINTE>SIGM_NOEU>

DEFORMATION>EPSI_NOEU>

CRITERES>SIEQ_NOEU>

RESULTAT:reslin0(MECA_STATIQUE)

• CREA_CHAMP

名前:temps>MODELE:model(AFFE_MODELE)>

TYPE_CHAM>NOEU_TEMP_R>

OPERATION>AFFE>

AFFE>TOUT:OUI>VALE>"400.0">NOM_CMP>"TEMP"

Output

• INPR_RESU

UNITE>ファイル名.med>RESU>

RESULTAT:reslin0(MECA_STATIQUE)>NOM_CHAM>DEPL>

RESULTAT:unnamed(CALC_CHAMP)>NOM_CHAM>EPSI_NOEU,SIEQ,NOEU,SIGM_NOEU>

FORMAT>MED

Astersutudyはどこかで設定してから出ないと選択できないものもあるため注意する

○Paravis

左側のパイプのブラウザーのbuiltinを右クリック>Open>.medファイルを開く

フィルター(画面上の黒いところ、ファイルとか編集とかあるところ)>Common>WarpByVector>オブジェクトインスペクター(画面左側)>ScaleFactorで調節すると膨れ具合がわかる

フィルター>Recent>ExtractGroup>SIGM_NOEU(節点の応力(6成分))

Salome-Meca2020_弾塑性解析†

参考Salome-Meca演習_弾塑性解析(2021)

○Geometry

いつも通り作成

○Mesh

いつも通り作成

○AsterStudy　　※上記のSalome-Meca2020_熱解析でも同様の内容は省く。

右クリックAddstage

Function and Lists

• DEFI_LIST_REEL(時間間隔の設定)

DEBUT>"0">

INTERVALLE>"1.0","Step length","0.1"入力

• DEFI_LIST_INST(解析オプション)

DEFI_LIST>LIST_INST:listr>

METHODE>MANUEL

• DEFI_FONCTION

名前:func>VALE>"0","0"と"1","1"入力>

NOM_PARA>INST

Material

• DEFI_MATERIAU

TRACTION>SIGM:point

Analysis

• STAT_NON_LINE

MODELE:model>

CHAM_MATER:fieldmat>

EXCIT>CHARGE:lode>FONC_MULT:func>

INCREMENT:times>

COMPORTEMENT>TOUT>RELATION>VMIS_ISOT_TRAC>DEFORMATION>SIMO_MIEHE>

NEWTON>MATRICE(マトリクス)>TANGENTE(タンジェント)orELASTIQUE(伸縮)>REAC_ITER>"1"入力>

CONVERGENCE>RESI_GLOB_RELA>"1.0e-06"任意で良い

CONVERGENCE>ITER_GLOB_MAXI>"100"任意で良い

Salome-Meca2020_弾塑性解析2†

○Geometry

いつも通り作成

○Mesh

いつも通り作成

○AsterStudy　　

右クリックAddstage

Function and Lists

• DEFI_LIST_REEL(時間間隔の設定)

DEBUT>"0">

INTERVALLE>"1.0","Step length","0.05"入力

• DEFI_LIST_INST(解析オプション)

DEFI_LIST>LIST_INST:listr>

METHODE>MANUEL

• DEFI_FONCTION

名前:func>VALE>"0","0"と"1","1"入力>

NOM_PARA>INST

• DEFI_FONCTION

名前:koufuku>VALE>"390","400"と"235","236"入力>

NOM_PARA>TEMP>

NOM_RESU>SIGM

Material

• DEFI_MATERIAU

ECRO_LINE>D_SIGM_EPSI:200入力>SY:235

• ELAS

E>NU>RHO>ALPHA

Analysis

• STAT_NON_LINE

MODELE:model>

CHAM_MATER:fieldmat>

EXCIT>CHARGE:lode>FONC_MULT:func>

INCREMENT:times>

COMPORTEMENT>TOUT>RELATION>VMIS_ISOT_TRAC>DEFORMATION>PETIT>ITER_INTE_MAXI>"20"入力>ITER_INTE_PAS>"0"入力

NEWTON>MATRICE(マトリクス)>TANGENTE(タンジェント)orELASTIQUE(伸縮)>REAC_ITER>"1"入力>

CONVERGENCE>RESI_GLOB_RELA>"0.0001"任意で良い

CONVERGENCE>ITER_GLOB_MAXI>"1000"任意で良い

Salome-Meca2020_モデル作成（Fuseする場合）†

○Geometry

1, BOX-1,2,3・・作成

2, Fuse-3objects選択

3, Partitionを作成、このときFuse/BOX-1,2,3を選択

4, Partitionのグループ作成、kotei/kanetuなど

5, Meshを作成

6, ジオメトリのグループ作成

以上

研究メモ†

• Windowsからstr/k2サーバーへは”FFFTP”を使用する.このとき、右クリック＞属性の変更　で”グループ”の呼出・書込と”その他”の呼出に✔がついてあるか確認必須

• 要素の種類…ビーム＝線、シェル＝面（厚さなし）、ソリッド＝立体

• モデルは半分で良い

• 主桁と主塔は節点は重なっても良いが、境界条件は別々で考える…２点必要かつマージしない

• マイダスの境界条件設定は「単一選択」で行う

• マイダスでの由利橋モデルにおける等分布荷重解析を行う際は、荷重＞方向＞グローバルZから「グローバルY」に変更

• 主塔水平材とは、両主塔を固定するための部材であり幅員に対して水平方向に存在する. 地震時のみに機能を発揮し、２本の主塔が個別振動することを防ぐ. また、軽微な損傷を許容することで主塔柱の発生断面力を制御する部材としても期待できる.

• プレテンション(引張力＋)とは、ワイヤロープに一定荷重を一定時間かけて、初期伸び(ロープの構造上の伸び)を除去する目的で行なわれる加工方法のこと. 初期伸びとロープ径の減少が少なくなる・伸びが少なくなる(弾性係数が高くなる・ロープのよりを安定した状態に落ち着かせて、動索として耐久性を向上させるという３つのメリットがある.

• 曲げ剛性とは部材の曲げやすさを表す値である.曲げ剛性（ＥＩ）＝縦ヤング係数（Ｅ）×断面二次モーメント（Ｉ）　から求めることができる. また、曲げ剛性(EI)と曲率半径(r)は比例関係にあり、曲げモーメント(M)と関係付け 1／r＝Ｍ／ＥＩ と表すことができる.

• 部材内部に生じる応力＝軸力・せん断力・曲げモーメント

• 軸力…部材の材軸方向に生じ、材を伸縮させようとする一対の力のことで引張力と圧縮力がある

• せん断力…部材に「ずれ」の変形を生じさせようとする一対の力で一方の力とは逆の向きにもう一方の力が働く

• 曲げモーメント…部材を曲げようとする一対のモーメントのことで一般に曲げモーメントを受けて曲がる部材の凸部分には引張力・凹部分には圧縮力が生じる

• 合成桁とは、床版と鋼桁とを"ずれ止め"で接合し合成床版と同じ様に一体として荷重に抵抗する桁形式のこと

• 普通鋼とは、炭素（C）・ケイ素（Si）・マンガン（Mn）・リン（P）・硫黄（S）の化学成分が含まれる鋼材のことである、SS材が代表的な素材

• 一方で、SM材はSはSteel（鋼）、Mは船舶を意味するMarineを指す、「溶接構造用圧延鋼材」という名称で溶接製の高い鋼材のこと

• 熱膨張係数とは、線膨張係数αもしくは体膨張係数βから求めることができる

• 主桁の数え方は左からG1,G2…とする

• マイダスで荷重組み合わせの生成方法は、マイダス〜解析編〜 の９－１９ページに記載

• 解析結果の変位などを見るやり方は、結果テーブル＞見たいものの選択＞「全て」＞✔＞適用＞エクセルヘ書き出し＞＞別のものを見たいときは「行フィルタリング」を選択

• マイダスで変位等をエクセルに書き出しができるため、そこから図やグラフの作成が可能

これが写真→→→

• 温度解析はサロメでやっても良い

• 由利橋モデルは、桁・主塔は梁要素、ケーブルはトラス要素で作成

• 梁要素とは、有限要素の一種であり、2点を結ぶ細長い軸形状を表現した要素のこと

• トラス要素とは、圧縮／引張の軸力のみを伝える部材のこと

• 弾塑性とは、降伏応力以下の力を加えた場合は力を除くと元に戻り、降伏応力以上の力を加えると力を除いても変形したままで元に戻らない性質のこと

• マイダスでは３Ｄモデルは”シェル”が該当する

• 粘性係数とは、流体の粘り気のこと、鋼材の動粘度ν＝0.5e-2cm2/s,密度ρ＝7.1g/cm3,温度1550℃

• 応力-ひずみ線図　　　　

• 座屈とは、ある一定の荷重を超えるとき柱は急激に折れ曲がることをいう

研究参考文献†

• 断面二次モーメント計算ソフト

• マイダス〜構造編〜

• マイダス〜解析編〜

• マイダスマニュアル

• ケーブルの基本理論

• 由利本荘市ハザードマップ

• 日本橋梁建設協会

• 鋼橋の補修・補強事例

• 鋼橋（上部工）の損傷事例

• 鋼橋の腐食事例

• 斜張橋ケーブル腐食を考慮した疲労寿命

• 斜張橋ケーブルとアンカーについて

• 斜張橋ケーブルの耐久性評価

• 橋梁の構造＆種類

• 鋼橋工事における火災リスクと防火対策

• 火災時における耐荷力に関する検討

• 火災を受けた橋梁の損傷と対応

• 火災を受けた橋梁の安全性評価と補修＿実例

• 火災を受けた橋梁の健全性評価に関する試験調査

• 鋼材の高温時における機械的特性

• SM400材について

• 普通鋼について

• SM材の耐熱能力

• 橋梁用厚鋼板

• 鋼材の材料定数

• 交通事故の原因ランキング

• 様々な交通事故統計

• 弾性変形と塑性変形

• トレス・ミーゼス応力の降伏条件

• トレス・ミーゼスの違い

• ミーゼス応力の定義

• モール・クーロンの破壊基準

• 面外と面内

• 様々な車高

• 鋼材の粘度

• ヒルの降伏応力

• 円筒座標系

• 直交座標系

• テンソルとは

• 火災を受けた橋梁の安全性評価

• 火災を受けた鋼橋の受熱温度推定に関する解析的研究

• 火災を受けた鋼橋の損傷調査とその対応について

• 熱分析

サロメページ†

• Salome-Meca演習_弾塑性解析(2021)

• Salome-Meca_熱応力メモ

• Salome-Meca2017

• Salome-Meca例題ファイル

• Salomeメモ←ここから線膨張係数入れられそう！

• 座屈解析

春休み課題の進捗†

• 一端固定他端単純支持梁

• モデル作成について

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

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

　　　　　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）」ここからアニメーションや振動数を見れる

やること†

• 理論値を求める・・済（ねじれは分からない）←分かった！！

• メッシュ1.0のharusindou_1の解析結果を検討・まとめ・・済

• 上記の相対誤差を求める・・済

• メッシュ0.8,0.9,2.0,3.0を解析→比較・検討・・済　

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

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

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

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

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

• どの梁形態を基本にするか決定する・・済　→☆☆「ローラー」

• 上記の決定した一端固定他端単純支持梁のメッシュ0.7~4.0を解析・検討・・済

• 上記の相対誤差の計算・・済

• ローラー上下２個の一端固定他端単純支持梁のメッシュ0.7~4.0を解析・検討・・済

• ローラー下１個の一端固定他端単純支持梁のメッシュ1.0~10.0の２次要素を解析・検討→このとき、固有振動数20個で解析・・済

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

• 上記の解析結果からねじれについても検討・考察・・済

• ローラー下１個の100mmの一端固定他端単純支持梁のメッシュ0.7~4.0を解析・検討・・済

• ねじれをどうするか話し合う→とりあえず後回し・・解決

• グラフ作成（９枚）・・済

• スライド作成・・済

• グラフの確認・使うデータ選別・・済

• スライド作成（長さ・振動の様子の２種類）・・済

• 収束値との相対誤差-要素数のグラフ作成→細長比で比較しても良さそう・・済

疑問・考察†

• ２次要素にすると鉛直・水平・ねじれモードは収束した。１次要素でも収束はした。

• １次要素だと水平モードの方が精度が良いが、２次要素だと鉛直モードの方が精度が良くなる。→間違い

• １次要素だと高次数の方が精度が良いが、２次要素だと低次数の方が良くなる。→間違い

• 細長いとねじれが出にくい。→100mmだと固有振動数１０コでねじれの出現あり。500mmだと固有振動数２０コでねじれの出現あり。

• 長さのみ・長さ＋幅・長さ＋高さを変えた条件でそれぞれ解析した

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

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

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

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

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

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

３年創造工房†

2023.10.27†

・タッチタイプの練習

・コマンドの練習

• コマンドの一覧 [#w7501507]

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

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

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

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

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

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

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

cp・・・コピー

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

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

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

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

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

〇コマンドの流れ

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

2023.11.10†

・コマンドの練習

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

• 作成した図 [#kf32a2c4]

使用データ：

引用元：岡田の卒論日誌（11/18課題）

2023.11.17†

・片持ち梁の解析（サロメ）

• サロメメモ [#ddb92132]

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

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

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

・測定値＝平均値ー理論値/理論値✕100(％)

・要素数＝ボリューム

・変位平均値（1.3.5.7）

・相対誤差（測定値）を計算

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

• 片持ち測定値 [#xcffb3ee]

• グラフ（理論値と測定値） [#y04fdb0d]

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

理論値：6.67mm

【理論式】

Pl^3/3EI

2023.11.24†

・単純梁（等方性１次）の解析

• 単純梁（等方性１次）測定値 [#g80425ab]

• グラフ（理論値と測定値） [#cb6b3f6d]

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

理論値：0.4167mm

【理論式】

Pl^3/48EI

2023.11.29†

・単純梁（異方性１次・等方性２次）の解析

＜等方性１次と異方性１次の比較＞ [#he11b38c]

• 単純梁（異方性１次）の測定値

• 比較のグラフ

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

理論値：0.4917mm

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

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

G：せん断弾性係数

K：せん断補正係数

このとき、［Pl/4KGA］= せん断項

＜等方性１次と等方性２次の比較＞ [#of7eeadc]

• 単純梁（等方性２次）の測定値

• 比較のグラフ

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

理論値：0.4167mm

２次要素を含む。

• メモ

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

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

2023.12.08†

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

• 単純梁（木材＋鋼材）の測定値 [#d0d51f4f]

• グラフ（理論値と測定値） [#k36a4b43]

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

理論値：0.099mm

２次要素を含む。

【理論式】

Pl^3/48EI ＋ Pl/4KGA

このとき、

EI = EI（木） ＋ EI（鋼）

G = G（木）

と考えることとする。

2023.12.15†

・Taxツールの概要・練習

• sibuw.tar.gzファイルの荷ほどき [#h7c60c92]

① 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では２回コンパイルする

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

＜メッシュの変更＞ [#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から、自分のファイルに転送

最新の100件