兼田の卒論日誌

研究†

研究テーマ†

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

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

やること†

• 参考文献を読む

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

• ☑リスクアセスメントの質問を考える

• ☑片持ち梁の片側側面に熱を与えてDX方向にどの程度変位が出るか見る、その熱ーヤング率グラフの作成　ー中間発表

• ☑主塔の燃える部分のヤング率(0,400,800℃の３パターン)をグラフを参考に下げて、＋ケーブルC1〜C12をくっつけモデルを作成し解析、ケーブルの応力及び降伏点はどのくらいのときが判別　ー中間発表

• ☑中間発表の概要作成・スライド作成　ー中間発表

• ☑健全モデルの完成〜9/28

• リスクアセスメントを考慮する上で必要となる条件を考える＆話し合い

進捗†

[見通し]

橋梁火災で鋼材（主塔・ケーブル）がどのように変化するかを調べ、全国・東北・秋田の県道で起きた事故・火災・トラックや危険車両が通る割合や想定封鎖期間などと合わせて「リスクアセスメント」へつなげていく

由利橋モデルの作成(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/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倍したものである.

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に書いてある

ヒルの降伏応力

円筒座標系

直交座標系

テンソルとは

hillファイル

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〜9/2)†

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

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

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

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

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

ビーム要素に変更する

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

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

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

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

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

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

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

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

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

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

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

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

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

○モデルのSalome設定一覧

　　　　

7/29:ケーブル２本でマトリスクエラーが出たため、ケーブル１本のでエラーを読み解く

  <S> Exception utilisateur levee mais pas interceptee.                              !
  ! Les bases sont fermees.                                                          !
  ! Type de l'exception : MatriceSinguliereError                                     !
  !                                                                                  !
  !    Arrêt pour cause de matrice non inversible.                                   !
  !    La base globale est sauvegardée. Elle contient les pas archivés avant         !
  ! l'arrêt.                                                                         !
  !                                                                                  !
  !    Conseils :                                                                    !
  !       - Vérifiez vos conditions aux limites.                                     !
  !       - Vérifiez votre modèle, la cohérence des unités.                          !
  !       - Si vous faites du contact, il ne faut pas que la structure ne "tienne"   !
  ! que par le contact.                                                              !
  !                                                                                  !
  !       - Parfois, en parallèle, le critère de détection de singularité de MUMPS   !
  ! est trop pessimiste ! Il reste néanmoins souvent                                 !
  !         possible de faire passer le calcul complet en relaxant ce critère        !
  ! (augmenter de 1 ou 2 la valeur du mot-clé NPREC) ou                              !
  !         en le débranchant (valeur du mot-clé NPREC=-1) ou en relançant le calcul !
  ! sur moins de processeurs. 

  ! <S> ユーザー例外が発生しましたが、インターセプトされませんでした！
  ! ベースが閉じられた！
  ! 例外タイプ: matrixSingularError ！
  ! !
  !    非可逆行列のため停止．
  !    グローバルデータベースが保存されました。グローバル・データベースが保存される！
  シャットダウン
  ! !
  !    ヒント: ！
  !       - 境界条件をチェックする！
  !       - モデルと単位の整合性をチェックする！
  !       - コンタクトを使用する場合、構造は "保持 "してはいけません！
  接触だけである！
  ! !
  !       - 並行して、MUMPSの特異点検出基準！
  特異点検出基準は悲観的すぎる！とはいえ、この基準を緩和することで
  とはいえ、この基準を緩和することで、完全な計算をパスすることができる！
  ! (NPRECキーワードの値を1または2増やす)、あるいは！
  切断する（NPRECキーワードの値を1にする）か、計算を再実行する！
  より少ないプロセッサで計算を再実行する。  

  Degré de liberté physique associé au noeud N27 et à la composante DRZ.

  Problème : la matrice est singulière ou presque singulière :
  Lors de la factorisation de la matrice, on a rencontré un problème
  (pivot nul ou presque nul) à la ligne 174 qui correspond au degré de liberté donné ci-dessus.

  Risques et conseils :
  * Si la ligne correspond a un degré de liberté physique, il s'agit probablement d'un mouvement
    de corps rigide mal bloqué.
    Vérifiez les conditions aux limites.
    Si vous faites du contact, il ne faut pas que la structure ne "tienne" que par le contact.
    Vérifiez également les caractéristiques matériaux (module d'Young, ...).
  Si la ligne correspond a un degré de liberté de Lagrange, il s'agit sans doute d'une condition
    limite redondante.
    En particulier, il se peut que la relation linéaire surabondante provienne des conditions de contact.
    Peut-être devriez vous exclure certains noeuds des conditions de contact
    (mots clés SANS_NOEUD et SANS_GROUP_NO).

  * Si le solveur utilisé est LDLT ou MULT_FRONT, vous pouvez utiliser le solveur MUMPS
    car celui-ci est le seul à pouvoir factoriser les matrices qui ne sont pas définies positives.

  * Parfois, en parallèle, le critère de détection de singularité de MUMPS est trop pessimiste ! Il reste néanmoins souvent
    possible de faire passer le calcul complet en relaxant ce critère (augmenter de 1 ou 2 la valeur du mot-clé NPREC) ou
    en le débranchant (valeur du mot-clé NPREC=-1) ou en relançant le calcul sur moins de processeurs.

  * Il se peut aussi que ce phénomène soit tout à fait normal avec X-FEM si la fissure passe
    très près d'un noeud.
    Si le nombre de décimales perdues n'est pas trop grand (max 10 décimales),
    vous pouvez relancer le calcul en augmentant le nombre de décimales perdues autorisé :
    mot-clé NPREC du mot clé facteur SOLVEUR.
    Sinon, contactez l'équipe de développement.
  > Exception utilisateur levee mais pas interceptee. !
  ! Les bases sont fermees.                               !
  ! Type de l'exception : MatriceSinguliereError          !
  ! [('?', (), (), ())] 

ノードN27とDRZコンポーネントに関連する物理的自由度。

問題：行列が特異または特異に近い：
行列を因数分解する際、174行目で問題が発生しました。
(ピボットがゼロまたはほとんどゼロ)が174行目で発生しました。

リスクとアドバイス
  * その行が物理的な自由度に対応している場合、それはおそらく、適切にロックされていない剛体の動きの結果です。
    運動。
    境界条件をチェックしてください。
    接触を使用している場合、接触だけで構造が「保持」されてはいけません。
    材料特性（ヤング率など）もチェックしてください。
この線がラグランジュの自由度に対応するのであれば、これは間違いなく冗長な境界条件である。
    境界条件である。
    特に、冗長な線形の関係は接触条件から来ているのかもしれません。
    おそらく、接触条件からいくつかの節点を除外する必要があるでしょう。
    (キーワードSANS_NOEUDとSANS_GROUP_NO)。

  * ソルバーがLDLTまたはMULT_FRONTの場合，MUMPSソルバーを使用することができます．
    これは，正定値でない行列を因数分解できる唯一のソルバーであるためである．

  * 並列の場合、MUMPSの特異点検出基準が悲観的すぎることがあります！それにもかかわらず
    しかし、この基準を緩和する（NPRECキーワードの値を1または2増加させる）か、または切断する（NPRECキーワードの値）ことにより、完全な計算をパスすることができる。
    (NPRECキーワードの値=1)、あるいは、より少ないプロセッサで計算を再開する。

  * クラックがノードの非常に近くを通過する場合、この現象はX-FEMでも完全に正常である可能性があります。
    節点に非常に近い場合。
    失われる小数点以下の桁数があまり大きくない場合（最大小数点以下 10 桁）、
    失われる小数点以下の桁数を増やして計算を再開することができます：
    SOLVEUR factor キーワードの NPREC キーワード。
    そうでなければ、開発チームに連絡してください。
  S> ユーザー例外がスローされたが、インターセプトされなかった！
  ! ベースはクローズされている！
  ! 例外タイプ: MatrixSingularError ！
  ! [('?', (), (), ())] 

N27はバネとケーブルの節点

バネの条件及び剛性がおかしいとマトリクスエラーが出る

短すぎる？長すぎる？

　　

8/5:ケーブル１本と２本で比較とエラーの改善

このときのAsterstudy

 ノードN27とDRZコンポーネントに関連する物理的自由度

がおかしいとエラーが出るため,試しにDRZを固定して見ると解析が回った

果たして「ピン」としての役割を果たせているのだろうか？

 Degré de liberté physique associé au noeud N2 et à la composante DRX.

 Problème : la matrice est singulière ou presque singulière :
 Lors de la factorisation de la matrice, on a rencontré un problème
 (pivot nul ou presque nul) à la ligne 12 qui correspond au degré de liberté donné ci-dessus.

 Risques et conseils :
  * Si la ligne correspond a un degré de liberté physique, il s'agit probablement d'un mouvement
    de corps rigide mal bloqué.
    Vérifiez les conditions aux limites.
    Si vous faites du contact, il ne faut pas que la structure ne "tienne" que par le contact.
    Vérifiez également les caractéristiques matériaux (module d'Young, ...).

  * Si la ligne correspond a un degré de liberté de Lagrange, il s'agit sans doute d'une condition
    limite redondante.
    En particulier, il se peut que la relation linéaire surabondante provienne des conditions de contact.
    Peut-être devriez vous exclure certains noeuds des conditions de contact
    (mots clés SANS_NOEUD et SANS_GROUP_NO).

  * Si le solveur utilisé est LDLT ou MULT_FRONT, vous pouvez utiliser le solveur MUMPS
    car celui-ci est le seul à pouvoir factoriser les matrices qui ne sont pas définies positives.

  * Parfois, en parallèle, le critère de détection de singularité de MUMPS est trop pessimiste ! Il reste néanmoins souvent
    possible de faire passer le calcul complet en relaxant ce critère (augmenter de 1 ou 2 la valeur du mot-clé NPREC) ou
    en le débranchant (valeur du mot-clé NPREC=-1) ou en relançant le calcul sur moins de processeurs.

  * Il se peut aussi que ce phénomène soit tout à fait normal avec X-FEM si la fissure passe
    très près d'un noeud.
    Si le nombre de décimales perdues n'est pas trop grand (max 10 décimales),
    vous pouvez relancer le calcul en augmentant le nombre de décimales perdues autorisé :
    mot-clé NPREC du mot clé facteur SOLVEUR.
    Sinon, contactez l'équipe de développement.

また、DX,DY,DZを固定するとDRXがおかしいと出る、なぜ？

8/7:エラー解決

DRZ=0にして良いらしい

それで回すとプレテンションは正しく入っている様子だった

次回は、まず弾塑性のみ、クリアしたら熱を静的解析でやってみる

↑これが帰省前の課題かな〜

8/11:熱が入らない・・

まず、プレテンションのみとそこに弾塑性を加えた場合だと、結果は変わらず解析は回った

しかし、そこに熱を入れると解析は回ったが、結果も変わらなかったため熱の影響が全くないことになっている

そのときのAsterstudyは以下の通りである

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

8/14:エラーの原因と改善が不明

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

これに

syutou = CREA_RESU(
 AFFE=_F(
   CARA_ELEM=elemprop, 
   CHAM_GD=temps, 
   LIST_INST=listr, 
   MODELE=model
 ), 
 NOM_CHAM='TEMP', 
 OPERATION='AFFE', 
 TYPE_RESU='EVOL_VARC'
) 

を付け加えたが解析できなかった

このまま回すとこのようなエラーが出る

  <S> Exception utilisateur levee mais pas interceptee.                        !
  ! Les bases sont fermees.                                                    !
  ! Type de l'exception : error                                                !
  !                                                                            !
  ! Erreur utilisateur dans la commande AFFE_MATERIAU / AFFE_VARC              !
  !   Pour la variable de commande TEMP                                        !
  !   la grandeur associée du champ doit être:  TEMP_R  mais elle est:  TEMP_F 

  <S> ユーザー例外が発生したが、インターセプトされなかった！
  ! ベースが閉じられた！
  ! 例外のタイプ: error ！
  ! !
  ! AFFE_MATERIAU / AFFE_VARCコマンドのユーザーエラー ！
  !   TEMPコマンド変数の場合！
  関連するフィールド変数はTEMP_Rであるべきですが、TEMP_Fです。

だから、TEMP_FをTEMP_Rに変えるとこのようなエラーが出る

  <S> Exception utilisateur levee mais pas interceptee.                             !
  ! Les bases sont fermees.                                                         !
  ! Type de l'exception : error                                                     !
  !                                                                                 !
  !  Erreur utilisateur (CREA_CHAMP/AFFE) :                                         !
  !    Le type du champ que l'on cherche à créer (réel, entier, complexe, fonction) !
  !    n'est pas compatible avec le mot clé utilisé (VALE, VALE_I, VALE_C, VALE_F). !
  !                                                                                 !
  !  Il faut respecter la correspondance suivante :                                 !
  !     - champ réel        -> VALE                                                 !
  !     - champ complexe    -> VALE_C                                               !
  !     - champ entier      -> VALE_I                                               !
  !     - champ fonction    -> VALE_F 

  <S> ユーザー例外が発生したが、インターセプトされなかった！
  ! ベースが閉じられた！
  ! 例外のタイプ: error ！
  ! !
  !  ユーザーエラー（CREA_CHAMP/AFFE）: ！
  !    作成しようとしているフィールドの型(実数、整数、複素数、関数) ！
  作成しようとしているフィールドのタイプ（実数、整数、複素数、関数）が、使用しているキーワード
 （VALE、VALE_I、VALE_C、VALE_F）と互換性がありません。!
  ! !
  !  次の対応を尊重しなければなりません： ！
  !     - 実フィールド -> VALE ！
  !     - 複素数フィールド -> VALE_C ！
  !     - 整数フィールド -> VALE_I ！
  !     - 関数フィールド -> VALE_F 

次に、TEMP_FをTEMP_Rに変えかつ実フィールド -> VALE に400℃と打ち込んで熱を与えて見ると、熱の影響が反映されない結果となってしまう

ただ、ケーブルにプレテンションのみが加わっている様子が見れた

→AFFE_VARCを2itemsにする必要がある

9/2:リスクアセスメントの今後の見通し

リスクアセスメントとは本来、橋梁を建設する前に行うものである

だから考え得るリスクは、由利橋で起きた事故や火災を調べるのではなく、今までの事故や火災から統計的に求める

①　全国・東北・秋田の県道で起きた事故及び火災、トラックや危険車両の通行割合を調べる

②　橋の上で起こった事故とその橋の通行量を調べる

③　由利橋の封鎖期間を想定する

④　鋼材と熱の関係について調べる＆解析する

⑤　④から熱によって低下するヤング率を求めて主塔の一部分のヤング率を下げておいたモデルでケーブルの降伏応力を見ていく

⑥　橋梁の全体図を作成して、⑤のような解析を行うことに加えて、燃える車高を変化させたり、車両の混雑具合を変化させたりと様々な状況のもとに起こる火災を比較する

これが、今後の自分のやるべきことである

これからは実際に熱を与える解析はしなくて良い→ヤング率を下げるという方法を取る

〇〇℃のときにどうなるかを見る静的解析を行う

一部のヤング率を下げたモデルで解析_midas(9/6～9/26)†

【Case1】

主塔0m~0.5mの材料　SM400

主塔0.5m~3mの材料　E=100000N/mm^2,ν=0.3,α=1.2e-05/,w=7.7e-08N/mm^3　→　この部分が火災に合う場所

主塔3m~50mの材料　SM400

ケーブルの材料　　E=195000N/mm^2,ν=0.3,α=1.26e-05/,w=7.85e-08N/mm^3

上記の結果を見たところ、どちらのパターンでも主塔高さ3mのところでポッキリと折れていることがわかる

【Case2】

主塔0.5m~3mの材料　E=150000N/mm^2,ν=0.3,α=1.2e-05,w=7.7e-08　→　この部分が火災に合う場所

ポッキリ折れた

!改善点として主塔の中に２枚のプレートを入れる

片側側面に使用したヤング率は上記のグラフを参考にした

色のついた高さ0.5m~3mの部分が火災想定（主塔片側のみ）

【Case3】～主塔の車道側側面のみを加熱した想定

ー400℃のとき

プレテンションの強い右の方に引っ張られる形になった

このときケーブル応力は、C1:-45.9N/mm^2,C12:28.2N/mm^2となった

9/24:青木ゼミ

まだ、途中…

引張応力と降伏応力の関係がよくわかっていない

交通事故・交通量・火災発生件数等をできる限り調べてみた

シェル要素片持ち梁と平面要素片持ち梁の比較とE-℃グラフ(9/9〜9/12)†

片持ち橋の寸法　幅4m×厚さ2m×長さ10m　　　温度500℃

σ＝E/ε から求めた

9/12:ヤング率と温度の関係を模索

上記の片持ち梁を使って熱と温度の関係を求めようとしたが、弾塑性を入れていないせいか線形グラフとなり求めることができなかった

そのため、青木さんからいただいた論文をもとにグラフを作成した

由利橋への影響を考える～県道交通事故・普通貨物車の通行量・橋上の事故から(9/18～9/20)†

日本の県道一覧

道路の種類

秋田県内橋梁一覧

令和３年度の一般交通量調査の結果について

このサイトのExcelファイルが秋田県内の交通量を調査した結果である

国道及び県道の交通量かつ大型車混入率が示されている

道路交通に関する統計

月別事故発生数と死者数が見比べられる

また、都道府県別交通事故発生状況を見ることもできる

死亡事故2023:file:///C:/Users/%E5%85%BC%E7%94%B0%E3%80%80%E5%A7%AB%E6%A2%A8/Downloads/%E6%AD%BB%E4%BA%A1%E4%BA%8B%E6%95%852023%20(1).pdf

秋田警察ホームページ

ホームページの「ミニ統計」からその年の秋田県県道の事故件数がわかる（円グラフから）

（令和4年用ミニ統計）

R4HP用ミニ統計

（令和5年の秋田県道事故発生件数） 交通統計

このpdfの7ページに道路別事故件数が記されている

東北管轄交通事故発生状況

このpdfの３に交通死亡事故の特徴が示されているが、ここから計算すると東北管轄の県道での令和5年の死亡事故件数は89件ほどであることがわかる

H17東北の事故発生状況

H13～H17の東北管轄県道事故件数棒グラフ

データは古いが東北管轄の県道の事故件数13073件だった

令和5年山形県道交通事故

山形県の場合は858件

総務省消防庁火災統計

年別の火災に関する統計諸々が記されている

下記が最新（令和5年）の火災統計

令和５年(１月～12 月)における火災の概要

下記は令和4年の短め資料

報道資料

火災曲線

火災時におよび火災後におけるスタッドの力学特性

図４に加熱温度ー時間関係のグラフあり

0,400,800℃の３モデルで解析(9/26～10/3)†

・健全(0),400,800の３パターンを，D+集中荷重/D+集中荷重+Lの２つの状況でそれぞれ解析

・ケーブル・桁・主塔の応力をそれぞれ見る ケーブル1160N/mm^2で伸び1570N/mm^2で降伏，桁235N/mm^2で降伏，主塔235N/mm^2で降伏

・プロパンの着火温度は約460℃～520℃に対して，ガソリンの着火温度は約210℃～300℃

プロパンガス・ガソリンの発火温度

ガソリンの発火温度

・伊勢湾岸自動車道の事故の写真や状況を参考に中間発表

伊勢湾岸自動車道爆発事故

・都市ガス普及率から秋田はプロパンガスの割合が高い

都市ガス普及率

モデルの全体像は以下の通りである

モデルの主塔と桁のつなぎ目はくっつけることに修正した

midas平面応力度について

midasリブについて

9/30:おかしい点

このように桁と主塔の節点のひとつ上の節点に大きな力がかかっている

つなぎ方に問題がなる？

全体モデルの作成(10/9〜10/25)†

全体モデルを作成した

以下のように，主塔と桁はシェル要素（板要素），ケーブルはビーム要素（梁要素）とした

桁と主塔は平面応力要素にするとケーブルとうまくつながらなかった

これは，平面応力要素にZ方向の応力が発生しないことが原因？よくわからない

これからは平面応力要素は使用せずに板要素を用いる

　

モデルの境界条件：桁　ローラー(x=-135)　ピン(x=55)、　主塔　完全固定(z=-20)

死荷重は全て圧力荷重で載荷

その際、面載荷ができなかったため、死荷重を/4して節点載荷した

主塔・桁断面はこのようにした

　　　　　　　　　　　　

主塔周りの４枚板要素は鋼材にすると主塔が降伏したため、可塑剤とした　　　　設定条件はこの通り　　　　　　　　　　　主塔のプレートはx軸方向３枚、y軸方向１枚　　

桁のメッシュは基本鋼材もコンクリートも1×1として、設計上必要な場合は0.75×1がある

○部材

主桁　19m×0.3m×190m　SM400

主塔　3m×1.5m×70m（高さ50m）　　SM400

ケーブル　D=200mm(C5〜C9),225mm(C1,C4,C10),315mm(C2,C3,C11,C12)

由利橋について†

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

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

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

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

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

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

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

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

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

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

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

令和4年：交通事故統計

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

本荘の雨温図

2022年：由利本荘市気象

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

・単純にスパン長が短い

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

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

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

・見通しが良い

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

落橋について†

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

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

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

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

ケーブル腐食について†

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

吊り橋…垂直方向の張力

がある。

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

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

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

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

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

これを犠牲防食という。

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

• 補足

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

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

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

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

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

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

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

[文献から…]

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

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

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

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

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

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

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

[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

健全度評価フロー👇

鋼材部の被災目安👇

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

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

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

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

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

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

運転者の心理について†

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

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

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

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

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

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

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

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

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

１位　安全不確認

２位　脇見運転

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

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

５位　交差点安全進行

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

７位　一時不停止

８位　信号無視

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

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

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

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

１，道路交通環境の整備

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

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

４，道路交通秩序の維持

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

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

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

１，車両相互追突

２，車両相互正面衝突

３，車両相互その他

４，車両単独工作物衝突

リスクアセスメントの質問†

・リスクの確率の計算方法（事故の確率や火災の確率）及び計算式は何を使うか

・死者数や被害者数の想定方法（人口に対して何割か）

・火災の起こる確率の式は何を使うか

・主塔への熱の伝達方法（動的解析か静的解析か）※今は静的解析をしています。ただし、エラーの解決が出来ていないので結果らしい結果は得られていません。

　　　→毎秒もしくは毎分、何度熱が加わりそれがどのくらい続くかという動的解析をするか

　　　→ある特定の温度のときに主塔がどのように変形し、かつケーブルの応力はどうなるかという静的解析

・由利本荘市や市民に対する経済への影響や交通への影響を考えたときに、被害を計算して確率として提示するか（○○/✕✕の確率で橋梁火災が起こった際に、○○/✕✕でこのような被害が起こるなど）、もしくは「このような被害が起こることへと繋がる」というような提示にするか（○○/✕✕の確率で橋梁火災が起こった際に、由利本荘市や市民に対して経済的損失や交通状況の悪化が起こるなど）、いずれにせよ、どう橋梁火災とつなげて経済・交通・費用などの「リスクアセスメントの提案」へと持っていくか

　　　→橋梁火災の規模やその発生原因、被害者数、被害状況、閉鎖期間などによってもそれぞれの確率は大きく変動すると思うため、自分でもどこから手を付けて良いのかわかっていません。また、自分で勝手に上記の項目について決めてもそれが現実的にあり得る確率・数値であるのかを判断できないため、アドバイスいただけると幸いです。

鋼構造を語る会（9/20）†

今日は鋼構造を語る会に参加した。 他の大学の人の研究を聞く機会がなかったため、とても面白かった。 現在、二酸化炭素を排出しない「グリーンスチール」というものがあるらしい。 最後のゲームでは３位になった！嬉しい！

青木ゼミ(9/24)†

今日は青木班のみんなと話し合いをした。

【今後の方向性】

デッキ（死荷重・活荷重）半分をつけて同様の解析 火災が起きる→この程度の事故率や火災率が考えられる→火災で主塔とケーブルにこのような影響が起こる→よって封鎖機関はこのくらいが想定される→迂回路が少ないから都会よりもダメージが大きいのでは？→リスクがあることも考える必要があるのではないか

【リスクアセスメントの浸透】

安全対策への考え方の違い 日本では事故を起こすことが悪だが、海外では事故が起こった時に通行止めや封鎖になることが悪

中間発表(10/3)†

夏休みが終わり中間発表を行った． スライドと概要に思ったよりも時間がかかったため，もう少し前もってやった方がよかった． 発表は時間内に終わらなかったため，もっと練習が必要である． やっとリスクアセスメントについて皆にまとめたものを見せることができたのはよかった．

【質問】

・リスクアセスメントの優先順位等はあるのか（ケーブルが重要だからケーブルのリスクから考えるなど）

・ケーブル腐食などもリスクアセスメントに含まれるのか

・由利橋に対してリスクアセスメントをしたいのか由利橋を通して橋梁全体に対してリスクアセスメントを提案したいのかよく分からなかった

・プロパンガスからの火災は爆発ではないのか

サロメモ†

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℃

• 応力-ひずみ線図　　　　

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

• 主要地方道は高速道路及び一般国道と一体になって広域交通を分担する幹線道路。一般県道は主要地方道以外の県道であり、地域活性化を図るための道路。

• 都市ガスは地下のガスパイプ通って各家庭へ

• 街道は一般国道の分類（羽州浜街道：うしゅうはまかいどう）

• 12時間交通量は午前7時〜午後7時までの交通量。24時間交通量は午前7時〜翌朝午前7時までの交通量。

• 軽トラックは車高1.95m、積載高さ2.5m以下という規定。

• 交通量が多い目安：4000台/日〜、大型車混入率目安：地方15％

• 大型車：バス・トラック・ダンプカー・ミキサー車など

• 小型貨物車：軽トラック・小型トラックなど（乗員定員10人以下）

• タンクローリーとは、ガソリンや水といった液体物や高圧ガスを移送するために使用するトラックまたはトレーラー

• 圧力荷重とは，、面上の特定の種類の静荷重のこと。 面上の静荷重に似ているが、表面に傾斜角度または鋭角で設定できる面上の基本的な静荷重とは異なり、圧力荷重の方向は常に面に対して垂直

研究参考文献†

• 色々なヤング率

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

• マイダス〜構造編〜

• マイダス〜解析編〜

• マイダスマニュアル

• ケーブルの基本理論

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

• 日本橋梁建設協会

• 鋼橋の補修・補強事例

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

• 鋼橋の腐食事例

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

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

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

• 橋梁の構造＆種類

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

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

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

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

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

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

• SM400材について

• 普通鋼について

• SM材の耐熱能力

• 橋梁用厚鋼板

• 鋼材の材料定数

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

• 様々な交通事故統計

• 弾性変形と塑性変形

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

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

• ミーゼス応力の定義

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

• 面外と面内

• 様々な車高

• 鋼材の粘度

• ヒルの降伏応力

• 円筒座標系

• 直交座標系

• テンソルとは

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

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

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

• 熱分析

• 各材料のメリット・デメリット

• 斜張橋の設計

• SM400Bの熱による耐力

• 鉄鋼材料の各温度における許容引張応力

• 高強度鋼

• 耐火性能

• 橋梁設計における手続きと照査

• タンクローリー車のサイズ

• 橋の設計の基本と基準の変遷

• 海外における耐震設計基準と橋梁の比較設計

• 橋梁火災事例

• 火災時における鋼合成桁の終局耐力と耐火対策に関する研究概要:file:///home/kouzou/%E3%83%80%E3%82%A6%E3%83%B3%E3%83%AD%E3%83%BC%E3%83%89/%E8%AB%96%E6%96%87%E5%86%85%E5%AE%B9%E8%A6%81%E6%97%A8%E3%80%81%E5%AF%A9%E6%9F%BB%E7%B5%90%E6%9E%9C%E8%A6%81%E6%97%A8.pdf

• 平面応力要素とは

• 板要素を利用したFEM解析

• midasフリーズ解決方法

• midas自重

• 移動荷重

• midas特殊解析

サロメページ†

• 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件