兼田の卒論日誌
をテンプレートにして作成
[
トップ
] [
新規
|
一覧
|
検索
|
最終更新
|
ヘルプ
|
ログイン
]
開始行:
#contents
*研究 [#bf04b5a7]
**研究テーマ [#x8e2544a]
-ケーブルを有する橋梁のリスクアセスメントの提案ー秋田県由...
--実橋において，その橋がどのような危険にさらされているか...
--リスクアセスメントは橋梁がどのくらいの割合で危険にさら...
**やること [#me638c2a]
--参考文献を読む
--☑由利橋のモデルをリニューアルする〜5/31
--☑由利橋の事故・落橋の原因・落橋したときの不便・ケーブル...
--☑マイダスで３Ｄモデルを作成する…「ケーブル橋梁ウィザー...
--☑プレテンションの設定方法を模索…単純梁✖簡易モデルでプレ...
--☑マイダスで由利橋のケーブルにプレテンション導入
--☑由利本庄市の天気・交通事故を調べる〜5/31
--☑リスクアセスメントのまとめ〜5/31
--☑橋梁と火災について調べる
--☑事故のときの運転者の心理について調べる
--☑マイダスで火災や温度変化の解析ができるかどうかを試す
--☑由利橋モデルに正しいプレテンションと死荷重・活荷重を入...
--☑死荷重及び活荷重をどう入れるか聞く（集中or等分布？、そ...
--☑主塔をソリッド✖シェルで作成
--☑由利橋で事故が起こらない理由を考える
--☑サロメで温度解析のやり方学ぶ
--☑サロメで車高部分のみを加熱してどうなるか
--☑鋼材に粘性を入れてぐにゃってなるか片持ち梁で確認→弾塑...
--☑リスクアセスメントの質問を考える
--☑片持ち梁の片側側面に熱を与えてDX方向にどの程度変位が出...
--☑主塔の燃える部分のヤング率(0,400,800℃の３パターン)をグ...
--☑中間発表の概要作成・スライド作成 ー中間発表
--☑健全モデルの完成〜9/28
--☑リスクアセスメントを考慮する上で必要となる条件を考える...
--☑要素温度荷重解析
--☑IABSEのスライド作成
--☑リスクの修正
* 進捗 [#pb2c767f]
**由利橋モデルの作成(5/29) [#v4fab7c2]
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kane...
マイダスで由利橋のモデルを作成.
材料・断面・境界条件を設定、割当した.
スパン長190m,主塔70m（桁上50m,20m）
ケーブル断面積D=200mmのケーブル位置
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ca...
ケーブル断面積D=225mmのケーブル位置
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ca...
ケーブル断面積D=315mmのケーブル位置
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ca...
**マイダスで荷重解析(6/4) [#tbe94155]
幅1m×高さ0.8m×スパン長10mの単純梁を作成し、集中荷重（節点...
全体図（(x,y,z)空間）
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kane...
横から図（(y,z)平面）
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kane...
上から図（(x,z)平面）
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kane...
集中荷重（節点６に荷重FZ=-10kN）の写真
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kane...
等分布荷重（節点１２～２２に荷重w=-1kN）の写真
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kane...
以下が解析結果である、変位図・Ｓ図・Ｍ図である.
〇変位図
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kane...
〇Ｓ図
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kane...
〇Ｍ図
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kane...
**簡易モデルでプレテンション導入＿改訂版(6/7) [#ac0f5e55]
千代岡さんにアドバイスをいただき、軸力ではなく「温度荷重...
簡易モデルは前回のと同じものを利用.
温度荷重をケーブルに与えた.
簡易モデル
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ka...
〇ケーブル材料
→弾性係数:1.95e+02
→ポアソン比:0.3
→線膨張係数:1.2e-05
→単位体積重量:7.7e-08
〇ケーブル断面積
→D=1mm
〇桁材料
→SM400
〇桁断面積
→10mm×10mm
解析を回してみたところ👇このように変形した.
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ka...
軸力あり・なしバージョンで比較したが、軸力は関係なさそう…
また、先端荷重Fを与えたときの先端変位は下記の表のようにな...
このとき、ｙ軸下方向を＋とする.
,F(N),温度荷重(℃),先端変位(mm)
,0,-100,-2.378
,0,-150,-3.567
,0,-200,-4.756
,-100,-100,-0.543
,-100,-150,-1.732
,-100,-200,-2.921
,-500,-100,+6.796
,-500,-150,+5.607
,-500,-200,+4.418
表からわかる通り、荷重ありとなしのときの変異差はどの温度...
F=-100N…1.835mm
F=-500N…9.174mm
であった.
この結果から、プレテンションは温度荷重から与えることがで...
加えて、温度荷重が相対的に大きいほど先端変位も大きくなる...
温度荷重と先端変位は比例関係にあると言える.
**由利橋モデルにプレテンション導入＿改訂版(6/7～6/10) [#m...
軸力ではなく、温度荷重からプレテンションを入れる方法を由...
このとき全く軸力は入れていない.
下記の写真通りに、外側のケーブル２本に-150℃・その他のケー...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/on...
①なんの荷重もかけずに解析を回すと👇このように変形した.
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/yu...
max変位＝-0.394mm
②次に中央（節点９６）に集中荷重F=-190kNを与えた場合（下記...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/sy...
結果は👇このように変形した.
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/yu...
max変位＝+1.852mm
③次に全桁にw=-1kNを与えた場合（下記の写真）
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/to...
結果は👇このように変形した.
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/to...
max変位＝+1.008mm
この結果から、プレテンションをのみの場合は桁が引っ張られ...
また、等分布荷重よりも集中荷重の方が変位が大きくなった.
きちんとプレテンションが入っていると思われる.
**主塔に熱解析(6/13) [#f4cc604c]
主塔のみを取り出して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℃
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/sy...
初期20℃→解析後520℃
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/sy...
初期20℃→解析後1020℃
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/sy...
結果からわかる通り、加える熱が高ければ高いほど変位は大き...
＋100℃のときと＋1000℃のときの変位はちょうど１０倍になって...
本来であれば400℃～500℃で主塔は耐えられなくなるが、この解...
この解決策として「弾塑性」を設定に加えなければならないら...
**プレテンション荷重の入れ方(6/14～6/17) [#afe8c12d]
今回はプレテンションを温度荷重ではなく。プレテンション荷...
手順は以下の通りである.
荷重＞プレテンション荷重＞荷重ケース選択＞kNを入力＞単一...
プレテンション荷重以外には何も与えずに解析を回したところ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/pr...
max変位＝-73.49mm
ケーブル（トラス要素）の引張（断面力）はこのようになった.
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/pr...
ケーブル（トラス要素）の応力度はこのようになった.
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ca...
次に桁を見ていく.
桁（梁要素）の断面力はこのようになった.
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/pr...
これは断面図である.
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/pr...
桁（梁要素）の応力度はこのようになった.
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ha...
これでケーブルプレテンションは入っているのではないか・・？
これが合っているのであれば、温度荷重への変換は必要なし！
！！結論：プレテンションは「プレテンション荷重」の項目か...
**プレテンションを温度荷重に変換する方法(6/14～6/17) [#m8...
線膨張係数αの計算式である👇これを使う.
α = $\frac{ΔL}{L}\frac{1}{ΔT}$
ΔL/L…ひずみ
ΔL…変化量
L…元の長さ
今日は、先にプレテンションを設定して解析して、要素作成タ...
すると、Ｃ１ケーブルに与える想定温度1250℃、Ｃ２ケーブルに...
大きすぎな気がする…
！！結論：プレテンションを温度荷重に変換する必要はない.
**活荷重の入れ方(6/17〜6/19) [#u6540e70]
活荷重の計算は桁の半分を対象にしている.
① 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)のＬ荷重
を付与.
下👇の写真の通り.
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/L_...
活荷重のみを与えた変形は👇このようになった.
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/Ll...
max変位 ＝ +35.36mm
**死荷重の入れ方(6/17) [#lc6fb415]
死荷重の計算は桁の半分を対象にしている.
D = 77.0kN/m^3 × (9.5m+6.6m) × 1.1m /2 = 681.835kN/m
この死荷重を等分布荷重として桁に与えた.
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/de...
まだカウンターウェイト等は考えずにすべてを鋼材として計算...
死荷重のみを与えた変形は👇このようになった.
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/61...
max変位＝+955.91mm
**由利橋モデルにプレテンション＋死荷重＋活荷重を組み合わ...
上記のプレテンション・死荷重・活荷重を全て由利橋モデルに...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/al...
max変位 ＝ +917.82mm
プレテンション・死荷重・活荷重を足し引きしたものとほぼ同値.
デフォルトで約1mもたわむものなのかと少し不安…
**主塔３Ｄモデル作成_midas(6/20～6/26) [#xa7a1e15]
マイダスで主塔の３Ｄモデルを作成した.
やり方は ウィザード＞基本構造＞シェル＞数値入力 からで...
手順は写真の通りである.
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/od...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/od...
ここでいう、
ｎ：ｙ軸方向の分割数
ｍ：ｘ軸方向の分割数
Ｉ：ｚ軸方向の分割数
である.
今回は、ｎ＝10コ,ｍ＝10コ,Ｉ＝50コ にした.なお適当.
設定した後は👇このようになった
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/sy...
全要素に400℃を与えて解析したところDX,DY,DZ方向全てに変位...
なお、まだ塑性や降伏点は設定していない.
〇DX方向変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/dx...
max変位＝7.62mm
〇DY方向変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/dy...
max変位＝4.20mm
〇DZ方向変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/dz...
max変位＝240.64mm
〇DXYZ方向変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/dx...
max変位＝240.77mm
一つ問題点…境界条件を節点を完全固定にしているため、３Ｄに...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ky...
失敗したやり方…節点を作成して要素で結んで、後から断面や材...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/so...
塑性材料を設定しようとすると「静的材料非線形関連機能をサ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/so...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ch...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ca...
**1/100主塔モデルで温度解析_Salome(7/3〜7/10) [#ke8a0fe8]
1/100モデルにした理由は、Asterstudyの設定以外でのエラーを...
7/3：合っているかはわからないがある程度Asterstudyの設定が...
残り、✅熱膨張係数（ELAS-ALPHA）・✅初期温度（AFFE_VARC-T...
を行えば良い.やり方わからなかった.聞かなければ、、
温度解析については[[Salomeメモ]]の熱関係と[[Salome-Meca_...
7/4：なんとか設定一通り終えて解析をした.成功もした.
結果を見た感じ変なところはなさそうだが、マイダスとの結果...
これは、マイダスはシェル（厚さ10mmの中は空洞）で、サロメ...
確認が必要かもしれない.
↑この理由は1/100モデルだからか、、忘れてた。。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/se...
結果
今回はマイダスと比較するために、先端部分の大体マイダスと...
上記写真の上端である.
なお単位はmmである.かつサロメは結果の数値を意図的に100倍...
,,マイダス,サロメ
,DX,7.62,9.4
,DY,4.20,5.8
,DZ,240.64,241
7/5:熱解析のやり方をまとめ、上記のグラフを出してみた.また...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/se...
結果的には、100℃のときの最大変位はDX=0.023mm,DY=0.014mm,D...
この値は400℃のときの最大変位の約1/4であることがわかった.
400℃のとき、DX=0.094mm,DY=0.058mm,DZ=2.41mm
このことから熱解析自体は正しくできているだろう.
次はいよいよ弾塑性を入れてみるか、、
7/9:弾塑性のやり方を学んだのと、実際の大きさで車高部分に4...
モデルは、下から0.5m,0.5m〜3m,3m〜50mの３分割にした.
車高は1.5mとし、そこから下に1mかつ上に1.5mの部分が真ん中...
このモデルの真ん中のBOXにのみ400℃与えた.
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ta...
結果はこのようになり、加熱部が一番大きく変形した.
Pravisの写真の点部分.
３つのBOXの手前左上.
,,DX,DY,DZ
,点0,7.48e-29,4.18e-29,9.19e-29
,点8,3.69,1.69,12.9
,点12,-0.014,0.0006,9.79
弾塑性入れてやってみたが破壊形態は見れなかった？見方がよ...
今は瞬間的に0℃→400℃に変えて解析をしているが、現実性を考え...
！補足！
マイダスと比較するために、サロメでも実モデル全体に熱を加...
,,マイダス(mm),サロメ(mm),相対誤差(％)
,DX,7.62,7.20,5.5
,DY,4.20,3.61,16.3
,DZ,240.64,241.1,0.19
熱解析自体は上手くいっていると思われる.
7/10:ゼミ
ゼミでただの鋼材に熱を加えるだけでは膨張するだけだと言わ...
そこで粘性係数を考慮する必要がある.
熱解析で鋼材が曲がるか確認するためにまずは簡単な片持ち梁...
**改めて、リスクアセスメントとは(7/17) [#e5d51923]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ri...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ri...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ri...
**粘性を入れた片持ち梁の解析(7/19〜7/24) [#td0548e7]
まず粘性の入れ方を学ぶために10×10010
まず粘性の入れ方を学ぶために10×100×10の片持ち梁で行った
[[CodeAster:https://code-aster.org/doc/default/en/man_u/u...
p137〜140に書いてある
[[ヒルの降伏応力:https://www.jsme.or.jp/jsme-medwiki/doku...
[[円筒座標系:https://www.cradle.co.jp/glossary/ja_A/detai...
[[直交座標系:https://www.cradle.co.jp/glossary/ja_T/detai...
[[テンソルとは:https://3.1415.jp/35ss9vvx]]
[[hillファイル:http://www.str.ce.akita-u.ac.jp/~gotouhan/...
CodeAsterを解読した
SalomeのLNTは下のように方向を表す
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/LN...
温度における、ac1＝準静的温度開始、ac3＝準静定温度終了
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ac...
これはAstersutudyの
DEFI_MTERIAU>META_LEMA_ANI(粘性)の設定画面である
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/an...
CodeAsterより等方性の場合はこのようになる👇
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/to...
γ、θ、Z方向は円筒座標系であり、x、y、z方向は直交（デカル...
この２つの座標系の関係性は
[x=γcosθ]
[y=γsinθ]
[z=z]
となる
①熱
②熱＋弾塑性
③熱＋粘性
④熱＋弾塑性＋粘性
①のPravis http://www.str.ce.akita-u.ac.jp/~gotouhan/j2...
②のPravis http://www.str.ce.akita-u.ac.jp/~gotouhan/j2...
①と③の結果は同じになり、②と④の結果は同じになった、つまり...
鋼材の動粘性を"0.5mm^2/s"として入れた
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ne...
[方向]
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) [#j842ba5a]
・ソリッド要素
・加熱する時間の設定１０〜３０分
・ケーブルの応力分布
・時間ごとに熱くなるのを再現できるか
7/25:ソリッド要素断念
主塔とC6・C7ケーブルをソリッド要素で作成を試みた
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ca...
しかし、下の写真のように主塔とケーブルのつなぎ目部分が上...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ca...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ca...
移動や回転も試してみたが、思った方向に変化しなかった
ビーム要素に変更する
7/26:ビームモデルの作成
[[code-asterによる解析 step-by-step:https://www.opencae.o...
[[構造要素の概要:https://www.opencae.or.jp/wp-content/upl...
ビーム要素でモデルを作成した
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/be...
ケーブルと主塔の間に小さいバネが必要であるため、以下のよ...
今回のモデルでは10mmの長さにした
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/be...
主塔やケーブルには架空の断面を設定する必要がある
主塔：Geometryの長方形フェース ケーブル：Geometryの円...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/sy...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ca...
写真ではケーブル断面を実際のケーブル位置に作成したが、そ...
原点に作っておいてAsterstudyで割当できるらしい
ビーム要素のも出るを作成するに当たっての注意がある
①架空の断面はGeometryでグループ作成してから、Meshでジオメ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ge...
②Meshのコンパウンド必須、これはビーム要素のみで断面はしな...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/co...
③②でしたCompound_Meshからグループを作成で「ノード」ー固定...
○モデルのSalome設定一覧
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/be...
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 le...
! l'arrêt. ...
! ...
! Conseils : ...
! - Vérifiez vos conditions aux limites. ...
! - Vérifiez votre modèle, la cohérence des unit...
! - Si vous faites du contact, il ne faut pas qu...
! que par le contact. ...
! ...
! - Parfois, en parallèle, le critère de détecti...
! est trop pessimiste ! Il reste néanmoins souvent ...
! possible de faire passer le calcul complet e...
! (augmenter de 1 ou 2 la valeur du mot-clé NPREC) ou ...
! en le débranchant (valeur du mot-clé NPREC=-...
! sur moins de processeurs.
! <S> ユーザー例外が発生しましたが、インターセプトされ...
! ベースが閉じられた！
! 例外タイプ: matrixSingularError ！
! !
! 非可逆行列のため停止．
! グローバルデータベースが保存されました。グローバ...
シャットダウン
! !
! ヒント: ！
! - 境界条件をチェックする！
! - モデルと単位の整合性をチェックする！
! - コンタクトを使用する場合、構造は "保持 "して...
接触だけである！
! !
! - 並行して、MUMPSの特異点検出基準！
特異点検出基準は悲観的すぎる！とはいえ、この基準を緩和...
とはいえ、この基準を緩和することで、完全な計算をパスす...
! (NPRECキーワードの値を1または2増やす)、あるいは！
切断する（NPRECキーワードの値を1にする）か、計算を再実...
より少ないプロセッサで計算を再実行する。
Degré de liberté physique associé au noeud N27 et à la...
Problème : la matrice est singulière ou presque singul...
Lors de la factorisation de la matrice, on a rencontré...
(pivot nul ou presque nul) à la ligne 174 qui correspo...
Risques et conseils :
* Si la ligne correspond a un degré de liberté physiqu...
de corps rigide mal bloqué.
Vérifiez les conditions aux limites.
Si vous faites du contact, il ne faut pas que la str...
Vérifiez également les caractéristiques matériaux (m...
Si la ligne correspond a un degré de liberté de Lagran...
limite redondante.
En particulier, il se peut que la relation linéaire ...
Peut-être devriez vous exclure certains noeuds des c...
(mots clés SANS_NOEUD et SANS_GROUP_NO).
* Si le solveur utilisé est LDLT ou MULT_FRONT, vous p...
car celui-ci est le seul à pouvoir factoriser les ma...
* Parfois, en parallèle, le critère de détection de si...
possible de faire passer le calcul complet en relaxa...
en le débranchant (valeur du mot-clé NPREC=-1) ou en...
* Il se peut aussi que ce phénomène soit tout à fait n...
très près d'un noeud.
Si le nombre de décimales perdues n'est pas trop gra...
vous pouvez relancer le calcul en augmentant le nomb...
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ま...
(NPRECキーワードの値=1)、あるいは、より少ないプロセ...
* クラックがノードの非常に近くを通過する場合、この現象...
節点に非常に近い場合。
失われる小数点以下の桁数があまり大きくない場合（最大...
失われる小数点以下の桁数を増やして計算を再開すること...
SOLVEUR factor キーワードの NPREC キーワード。
そうでなければ、開発チームに連絡してください。
S> ユーザー例外がスローされたが、インターセプトされな...
! ベースはクローズされている！
! 例外タイプ: MatrixSingularError ！
! [('?', (), (), ())]
N27はバネとケーブルの節点
バネの条件及び剛性がおかしいとマトリクスエラーが出る
短すぎる？長すぎる？
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/n2...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/mm...
8/5:ケーブル１本と２本で比較とエラーの改善
[[このときのAsterstudy:http://www.str.ce.akita-u.ac.jp/~g...
ノードN27とDRZコンポーネントに関連する物理的自由度
がおかしいとエラーが出るため,試しにDRZを固定して見ると解...
果たして「ピン」としての役割を果たせているのだろうか？
Degré de liberté physique associé au noeud N2 et à la c...
Problème : la matrice est singulière ou presque singuli...
Lors de la factorisation de la matrice, on a rencontré ...
(pivot nul ou presque nul) à la ligne 12 qui correspond...
Risques et conseils :
* Si la ligne correspond a un degré de liberté physiqu...
de corps rigide mal bloqué.
Vérifiez les conditions aux limites.
Si vous faites du contact, il ne faut pas que la str...
Vérifiez également les caractéristiques matériaux (m...
* Si la ligne correspond a un degré de liberté de Lagr...
limite redondante.
En particulier, il se peut que la relation linéaire ...
Peut-être devriez vous exclure certains noeuds des c...
(mots clés SANS_NOEUD et SANS_GROUP_NO).
* Si le solveur utilisé est LDLT ou MULT_FRONT, vous p...
car celui-ci est le seul à pouvoir factoriser les ma...
* Parfois, en parallèle, le critère de détection de si...
possible de faire passer le calcul complet en relaxa...
en le débranchant (valeur du mot-clé NPREC=-1) ou en...
* Il se peut aussi que ce phénomène soit tout à fait n...
très près d'un noeud.
Si le nombre de décimales perdues n'est pas trop gra...
vous pouvez relancer le calcul en augmentant le nomb...
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:熱が入らない・・
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/pr...
まず、プレテンションのみとそこに弾塑性を加えた場合だと、...
しかし、そこに熱を入れると解析は回ったが、結果も変わらな...
そのときのAsterstudyは以下の通りである
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/be...
8/14:エラーの原因と改善が不明
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/be...
これに
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 / ...
! Pour la variable de commande TEMP ...
! la grandeur associée du champ doit être: TEMP_R ...
<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, ...
! n'est pas compatible avec le mot clé utilisé (VAL...
! ...
! 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) [#i...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/no...
【Case1】
主塔0m~0.5mの材料 SM400
主塔0.5m~3mの材料 E=100000N/mm^2,ν=0.3,α=1.2e-05/,w=7.7e...
主塔3m~50mの材料 SM400
ケーブルの材料 E=195000N/mm^2,ν=0.3,α=1.26e-05/,w=7.85...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/E1...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/E1...
上記の結果を見たところ、どちらのパターンでも主塔高さ3mの...
【Case2】
主塔0.5m~3mの材料 E=150000N/mm^2,ν=0.3,α=1.2e-05,w=7.7e-...
ポッキリ折れた
!改善点として主塔の中に２枚のプレートを入れる
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/E-...
片側側面に使用したヤング率は上記のグラフを参考にした
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ka...
色のついた高さ0.5m~3mの部分が火災想定（主塔片側のみ）
【Case3】～主塔の車道側側面のみを加熱した想定
ー400℃のとき
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/40...
プレテンションの強い右の方に引っ張られる形になった
このときケーブル応力は、C1:-45.9N/mm^2,C12:28.2N/mm^2とな...
9/24:青木ゼミ
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/c_...
まだ、途中…
引張応力と降伏応力の関係がよくわかっていない
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ma...
交通事故・交通量・火災発生件数等をできる限り調べてみた
**シェル要素片持ち梁と平面要素片持ち梁の比較とE-℃グラフ(9...
片持ち橋の寸法 幅4m×厚さ2m×長さ10m 温度500℃
,,シェル要素(プレートなし),平面要素(プレートなし),平面要...
,DZ方向変位(mm),30.985,30.985,21.380,21.380
,応力(N/mm^2),312.0,312.6,265.8,265.9
,ヤング率(GPa),101,101,124,124
σ＝E/ε から求めた
9/12:ヤング率と温度の関係を模索
上記の片持ち梁を使って熱と温度の関係を求めようとしたが、...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ka...
そのため、青木さんからいただいた論文をもとにグラフを作成...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/E-...
**由利橋への影響を考える～県道交通事故・普通貨物車の通行...
[[日本の県道一覧:https://altcva.main.jp/pref/]]
[[道路の種類:https://www.mlit.go.jp/road/sisaku/dorogyous...
[[秋田県内橋梁一覧:https://river.longseller.org/bp/5.html]]
[[令和３年度の一般交通量調査の結果について:https://www.pr...
このサイトのExcelファイルが秋田県内の交通量を調査した結果...
国道及び県道の交通量かつ大型車混入率が示されている
[[道路交通に関する統計:https://www.e-stat.go.jp/stat-sear...
月別事故発生数と死者数が見比べられる
また、都道府県別交通事故発生状況を見ることもできる
[[秋田警察ホームページ:https://www.police.pref.akita.lg.j...
ホームページの「ミニ統計」からその年の秋田県県道の事故件...
（令和4年用ミニ統計）
[[R4HP用ミニ統計:https://www.police.pref.akita.lg.jp/uplo...
（令和5年の秋田県道事故発生件数）
[[交通統計:https://www.police.pref.akita.lg.jp/uploads/co...
このpdfの7ページに道路別事故件数が記されている
[[東北管轄交通事故発生状況:https://www.tohoku.npa.go.jp/t...
このpdfの３に交通死亡事故の特徴が示されているが、ここから...
[[H17東北の事故発生状況:https://www.thr.mlit.go.jp/road/i...
[[H13～H17の東北管轄県道事故件数棒グラフ:https://www.thr....
データは古いが東北管轄の県道の事故件数13073件だった
[[令和5年山形県道交通事故:https://www.pref.yamagata.jp/do...
山形県の場合は858件
[[総務省消防庁火災統計:https://www.fdma.go.jp/pressreleas...
年別の火災に関する統計諸々が記されている
下記が最新（令和5年）の火災統計
[[令和５年(１月～12 月)における火災の概要:https://www.fdm...
下記は令和4年の短め資料
[[報道資料:https://www.fdma.go.jp/pressrelease/houdou/ite...
[[火災曲線:https://www.oit.ac.jp/japanese/sangaku/OITID-2...
[[火災時におよび火災後におけるスタッドの力学特性:https://...
図４に加熱温度ー時間関係のグラフあり
**0,400,800℃の３モデルで解析(9/26～10/3) [#vc55d915]
・健全(0),400,800の３パターンを，D+集中荷重/D+集中荷重+L...
・ケーブル・桁・主塔の応力をそれぞれ見る
ケーブル1160N/mm^2で伸び1570N/mm^2で降伏，桁235N/mm^2で降...
・プロパンの着火温度は約460℃～520℃に対して，ガソリンの着...
[[プロパンガス・ガソリンの発火温度:https://www.elg-inc.jp...
[[ガソリンの発火温度:https://www.yellowhat.jp/column/oil/...
・伊勢湾岸自動車道の事故の写真や状況を参考に中間発表
[[伊勢湾岸自動車道爆発事故:https://www.khk.or.jp/public_i...
・都市ガス普及率から秋田はプロパンガスの割合が高い
[[都市ガス普及率:https://uub.jp/pdr/eg/gas.html]]
モデルの全体像は以下の通りである
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ze...
モデルの主塔と桁のつなぎ目はくっつけることに修正した
[[midas平面応力度について:https://manual.midasuser.com/JP...
[[midasリブについて:https://support.midasuser.com/hc/ja/a...
9/30:おかしい点
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/st...
このように桁と主塔の節点のひとつ上の節点に大きな力がかか...
つなぎ方に問題がなる？
**全体モデルの作成(10/9〜10/25) [#pa80f758]
全体モデルを作成した
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/mo...
以下のように，主塔と桁はシェル要素（板要素），ケーブルは...
桁と主塔は平面応力要素にするとケーブルとうまくつながらな...
これは，平面応力要素にZ方向の応力が発生しないことが原因？...
これからは平面応力要素は使用せずに板要素を用いる
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/mo...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/si...
モデルの境界条件：桁 ローラー(x=-135) ピン(x=55)、 主...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/si...
死荷重は全て圧力荷重で載荷
その際、面載荷ができなかったため、死荷重を/4して節点載荷...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/da...
主塔・桁断面はこのようにした
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ka...
主塔周りの４枚板要素は鋼材にすると主塔が降伏したため、可...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/me...
桁のメッシュは基本鋼材もコンクリートも1×1として、設計上必...
○部材
主桁 19m×0.3m×190m SM400
主塔 3m×1.5m×70m（高さ50m） SM400
ケーブル D=200mm(C5〜C9),225mm(C1,C4,C10),315mm(C2,C3,C1...
**健全時・400℃・800℃のときの解析結果(11/6～11/7) [#s93d49...
○健全時と400℃の比較
あまり大差ない結果となった
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/23...
上記のC1ケーブルと桁との節点が最も応力のかかった結果とな...
,,健全,400,800
,集中荷重,220MPa,224MPa,224MPa
,集中荷重＋分布荷重,225MPa,230MPa,230MPa
,変位,-914mm,-934mm,-934mm
ヤング率を下げているにも関わらず変位と応力が400℃，800℃の...
また，ヤング率を下げた主塔と下げていない主塔とで変形が重...
上手くヤング率が反映されていない？
また主塔の0.5m～3.0mを加熱しているためか，下の写真のよう...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/22...
**IABSEに向けてリスクの定量化(11/12～11/26) [#dcbc9f10]
[[秋田県の平均通勤時間:https://statresearch.jp/house/rent...
ここから計算すると全国平均通勤時間：25.3分、秋田平均通勤...
[[令和4年道路形状別交通事故発生件数:https://www8.cao.go.j...
[[令和5年道路形状別交通事故発生件数:https://www8.cao.go.j...
この2つのサイトは交差点・単路・カーブといった道路形状別の...
由利橋付近は一般単路と交差点付近の値が参考になるだろう
[[交通事故が起こりやすい道路形状:https://agoora.co.jp/jik...
これも同様内容、参考までに
3年分のデータから一般単路は大体31%くらい
[[秋田県交差点での事故率:https://www.sonpo.or.jp/about/us...
秋田県のデータ、令和5年のもの
[[H8~H17間の全国の橋梁上の事故件数:https://www.mlit.go.jp...
件数：58,878 延長：8,867km 延長1km当たりの事故発生...
[[事故類型別人身事故発生件数割合:https://www.ms-ins.com/s...
ほとんどが相互車両の事故
[[高速道路の事故率:https://www.mlit.go.jp/road/ir/ir-coun...
ここから分かる通り、自動車専用道路よりも一般道路の方が事...
歩行者等の人がいると予期せぬ行動が起こりやすいため高くな...
[[事業用自動車の走行距離と交通事故件数、交通事故死者数:ht...
[[令和3年全国市区町村別交通事故死者数:https://www.itarda....
↑これも参考までに
[[令和５年中の道路交通事故の状況:https://www8.cao.go.jp/k...
高齢者の死亡事故や事故に逢った状況が詳しく書いてある
[[消防白書:https://www.fdma.go.jp/publication/hakusho/r5/...
p28に記載
[[車両別事故・火災情報:https://renrakuda.mlit.go.jp/renra...
! 車種別事故・火災情報件数及び割合は、貨物車が 519 件...
353 件（27.0%）となっている。事象別にみると、事故情報件数...
（43.4%）と最も多く、次いで乗用車が 57 件（24.3%）となっ...
物車が 417 件（38.9%）と最も多く、次いで乗用車が 296 件（...
月末現在における乗用車の保有台数（38,938 千台）は全保有台...
り、乗用車の事故・火災情報件数は乗用車の普及状況が反映さ...
らに、貨物車は保有台数 6,175 千台（7.4%）に対して事故・火...
逆に軽乗用車は保有台数 23,388 千台(28.2%)に対して事故・火...
くなっている。
p5抜粋
この試資料はなんの事故なのか・なんの火災なのかがわからない
[[日本の総道路延長:https://www.mlit.go.jp/road/soudan/sou...
令和3年のデータ
都道府県道の総道路延長：142,942.2 km
秋田県道：2376.7km
[[2023道路統計年報:https://www.mlit.go.jp/road/ir/ir-data...
道路統計年報では、令和4年時点での道路及び橋梁現状を閲覧で...
[[交通事故件数減少推移グラフ:https://www.jama.or.jp/opera...
交通事故発生件数の減少推移を閲覧できる
・交通事故で死ぬ確率ー1年間に10万人あたり4人程度
・一生のうちに交通事故に遭う確率ー25％ 四人に一人...
・一生のうちに交通事故で死亡する確率ー1％ 年間約1万人...
・年間10万人あたり４人が交通事故で亡くなっている
・約0.2％（500人に1人）が交通事故に遭っている
・約0.4％（250人に1人）が1年間に交通事故を起こしている ...
[[交通事故に逢う確率＆起こす確率:https://chokottokyosai.j...
以上の結果から定量化に努める
○県道での事故発生率
○全国橋梁上での事故率
○大型車混入率
○一般単路の事故率・交差点付近の事故率
○全事故のうち、死亡者や重傷者が出る割合
○通勤時間の増加
などなどを考慮する
井上さんからいただいた最終的な5.75×10＾-7/km/年はどう計算...
Anser:採用された仮定を全てかけ合わせる
この計算方法がわかれば、あとは橋長190m(0.19km)をかけて/年...
その計算した確率が1×10^-6/年よりも小さければリスクを考慮...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ri...
リスクの定量化を行った
最終的な由利橋のリスクは 2.02×10^-7/年 になった
「リスクを高い低いと判断することよりも、場面に合わせたリ...
と結論付けるの方がまとめやすいかもしれない・・・
一応基準は人が1年間で死ぬ確率 4×10^-5/年 にした
今回は思っていた通りこのリスクよりも小さい
※○印のついたものが定量化に使用した値
[計算式↓]
7.82E-11/km^2・人・年 × 橋長：0.19km^2 × 由利本荘市...
11/22:青木さんと発表練習後
由利橋のリスクはそもそも死ぬ確率ではないので、交通事故死...
トルコの事例と比較の方がいいかもしれん
0.3/km→0.3,0.15/km→0.15,0.09/person→0.09/day×365ay=32.85
それぞれ修正
全てを確立に変換して、○○○○/km・年を出す
そして、最終的に橋長をかけて○○○○/年を求める
これが1年間に起こり得る由利橋でのリスクとなる
まとめでは、解析結果から火災が起こるリスクは低いかもしれ...
11/26:リスクアセスメントの改訂版
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ri...
**要素温度荷重解析の結果(11/13〜11/14) [#k9303318]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/E-...
これからは、ヤング率を下げて解析するのではなく”要素温度荷...
ただし、この方法は徐々に熱が加わっていくことを想定してい...
これは今までのヤング率を下げたものを使っていたときと同様...
健全モデルのメッシュをもっと細かくする必要あり
解析結果は以下の通りである
,,健全,250,270,300,400,800,1200
,主塔応力（集中荷重）,78.3MPa,316.9MPa,339.4MPa,373MPa,48...
,主塔応力（集中荷重＋分布荷重）,80.4MPa,317.6MPa,340MPa,3...
,C1応力（集中荷重）,220MPa,220MPa,220MPa,220MPa,220MPa,22...
,C1応力（集中荷重＋分布荷重）,226MPa,226MPa,226MPa,226MPa...
,C12応力（集中荷重）,93.1MPa,93.3Mpa,93.3MPa,93.3MPa,93.4...
,C12応力（集中荷重＋分布荷重）,95.3MPa,95.45MPa,95.5MPa,9...
,x軸変位（集中荷重）,-47.4mm,-47.7mm,-47.71mm,-47.71mm,-4...
,x軸変位（集中荷重＋分布荷重）,-49.14mm,-49.15mm,-49.15mm...
,y軸変位（集中荷重）,3.16mm,203.72mm,219.77mm,243.83mm,32...
,y軸変位（集中荷重＋分布荷重）,3.05mm,203.61mm,219.66mm,2...
,z軸変位（集中荷重）,-891.16mm,-891.55mm,-891.59mm,-891.6...
,z軸変位（集中荷重＋分布荷重）,-914.35mm,-914.73mm,-914.7...
以下３枚が270℃で降伏した時の変形図である
これを見たらわかるように、主塔が桁の外側へ倒れ込むような...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/27...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/27...
以下２枚が主塔と桁（板要素）の応力分布である
２枚目からわかる通り、主塔の実際に加熱した部分が270℃以降...
写真左の方が応力が集中しているのは、メインスパンの方が長...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/27...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/27...
グラフ１枚目[heat-displacement]は熱とy軸方向変位の関係グ...
ただ、ケーブル応力はどの温度のときでもほとんど変わらなか...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ya...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/he...
グラフ２枚目[heat-stress]は熱と主塔応力の関係グラフであり...
270℃は大体10分程度で到達されると予想
消防車が来るまでの平均時間は約８分、消火活動の準備などで...
一般の鋼材は400℃で降伏するとされているが、だいぶ低い段階...
これは主塔の厚さや死荷重、活荷重、プレテンションなどが相...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/27...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/he...
なお、解析ではケーブルの応力はほとんど変化しない結果とな...
○補足
初期温度20度とすると400℃へ行くまでに約20分ほどかかること...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/he...
このことから、実際に主塔が降伏するまで火災が続く可能性は...
熱と時間に関しては要検討！
↑これ間違いかも…質量計算し直すと270℃まで200分かかると出て...
**ケーブル火災解析(12/18) [#nd8a7ba3]
ケーブルの要素分割をすると、ぽきぽきになってしまい解析で...
つなぎ目をどうするか考えて、プレテンションの入れ方を模索...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/a1...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/b1...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/c1...
ケーブル火災解析は、卒論後に検討でもいいかも
**リスクの修正(12/20〜2/6) [#xe27b376]
今までのまとめとリスクの修正を記した．
ユーザーリスクと通行止めリスクの２つを考え，それぞれを検...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ri...
[Eurocodeについて]
[[Eurocode:https://eurocodes.jrc.ec.europa.eu/]]
[[欧州構造基準（ユーロコード）の制定─その体系と内容が及ぼ...
[[ユーロコードの完全実施と我が国への影響:https://www.pwrc...
[[ユーロコードの構成を俯瞰:https://note.com/issey96/n/n5e...
[[内閣府 令和3年中の道路交通事故の状況:https://www8.cao....
主塔で起こる火災を
事故が原因の場合・・正面衝突事故率×橋梁に対する主塔の割合
設備不備の場合・・車両火災中設備からの出火事故率×橋梁に対...
とか？？
中間発表後にリスクの修正を行った
〇参考資料
[[大型車の車両火災発生状況について:https://www.ntsel.go.j...
[[The influence of wind and the spatial layout of dwellin...
[[Effect of wind speed and direction on facade fire sprea...
[[Numerical study of effects of wind on the vertical fire...
[[Experimental study of the temperature characteristics o...
[[Influence of Wind Direction on Fire Spread on High-Rise...
これに「図4. 鉛直風速の変化則」に風速と高さとの関係グラフ...
***限界状態について [#va8f246b]
[[道路橋示方書:https://www.mlit.go.jp/road/sign/kijyun/pd...
限界状態がどれであれ通行止めは発生する。
そのため、特にリスクアセスメントでは限界状態を考慮してい...
1) 橋の限界状態 1
橋としての荷重を支持する能力が損なわれていない限界の状態
2) 橋の限界状態 2
部分的に荷重を支持する能力の低下が生じているが，橋として...
能力に及ぼす影響は限定的であり，荷重を支持する能力があら...
にある限界の状態
3) 橋の限界状態 3
これを超えると構造安全性が失われる限界の状態
1) 死荷重 (D)
2) 活荷重 (L)
3) 衝撃の影響 (I)
4) プレストレス力 (PS)
5) コンクリートのクリープの影響 (CR)
6) コンクリートの乾燥収縮の影響 (SH)
7) 土圧 (E)
8) 水圧 (HP)
9) 浮力又は揚圧力 (U)
10) 温度変化の影響 (TH)
11) 温度差の影響 (TF)
12) 雪荷重 (SW)
13) 地盤変動の影響 (GD)
14) 支点移動の影響 (SD)
15) 遠心荷重 (CF)
16) 制動荷重 (BK)
17) 橋桁に作用する風荷重(WS)
18) 活荷重に対する風荷重(WL)
19) 波圧 (WP)
20) 地震の影響 (EQ)
21) 衝突荷重 (CO)
22) その他
1) 永続作用による影響が支配的な状況（永続作用支配状況）
①D +PS+CR+SH+E+HP+(U) +(TF) +GD+SD +WP +(ER)
2) 変動作用による影響が支配的な状況（変動作用支配状況）
②D+L+I+PS+CR+SH+E+HP+(U) +(TF)+(SW)+GD+SD+(CF)+(BK) +WP +...
③D +PS+CR+SH+E+HP+(U)+TH+(TF) +GD+SD +WP +(ER)
④D +PS+CR+SH+E+HP+(U)+TH+(TF) +GD+SD +WS +WP +(ER)
⑤D+L+I+PS+CR+SH+E+HP+(U)+TH+(TF)+(SW)+GD+SD+(CF)+(BK) +WP...
⑥D+L+I+PS+CR+SH+E+HP+(U) +(TF) +GD+SD+(CF)+(BK)+WS+WL+WP ...
⑦D+L+I+PS+CR+SH+E+HP+(U)+TH+(TF) +GD+SD+(CF)+(BK)+WS+WL+W...
⑧D +PS+CR+SH+E+HP+(U) +(TF) +GD+SD +WS +WP +(ER)
⑨D +PS+CR+SH+E+HP+(U)+TH+(TF)+(SW)+GD+SD +WP+EQ +(ER)
⑩D +PS+CR+SH+E+HP+(U) +(TF) +GD+SD +WP+EQ +(ER)
3) 偶発作用による影響が支配的な状況（偶発作用支配状況）
⑪D +PS+CR+SH+E+HP+(U) +GD+SD +EQ
⑫D +PS+CR+SH+E+HP+(U) +GD+SD +CO
(3) (2)1)から 3)に規定する作用の組合せに対して，表-3.3.1 ...
重係数を考慮する。
γ୮ ：荷重組合せ係数であり，異なる作用の同時載荷状況に応じ...
考慮する作用の規模の補正を行うための係数。
γ୯ ：荷重係数であり，作用の特性値に対するばらつきに応じて...
慮する作用の規模の補正を行うための係数。
なお，活荷重に対する衝撃の影響（I）を考慮するにあたって，...
は荷重組合せ係数γ୮及び荷重係数γ୯を乗じる必要はない。
(4) 風荷重については必要に応じて他の作用を考慮しない場合...
を適切に設定する。
(5) 衝突荷重及び制動荷重については死荷重及び活荷重のみと...
(2)1)2)以外の条件を適切に設定する。
***由利橋について [#u307a94f]
由利橋はH23.2.22に施工が始まり、H25.1.31に完成した。
橋長190.5mの不等径間の２径間連続斜張橋（鋼床版２主箱桁）...
秋田県由利本荘市の子安川にまたがる形で位置している。
由利本荘市：面積1209.59km^2,人工:74707人
小吉川：一級河川,全長61km,流域面積1190km^2,小洪水は毎年起...
今までで由利橋で一度も事故が起こったことはない！！→なぜ？...
予想可能な災害…土砂崩れ（出羽丘陵）、氾濫及び洪水（小吉川...
[[由利本荘市ハザードマップ:https://www.city.yurihonjo.lg....
考え得る危険事項…交通事故（追突・接触）、☆火災（自動車事...
使用不可になった場合の懸念点…移動手段の減少による渋滞・移...
※このときの迂回路…飛鳥大橋、本荘大橋、…子吉川が流れている...
[[令和4年：交通事故統計:https://www.police.pref.akita.lg....
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/a....
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/b....
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/c....
[由利橋で事故が起こらない理由]
・単純にスパン長が短い
・車両の速度が比較的遅い
・ほとんどが軽自動車か普通自動車、たまに軽トラックであっ...
・渡り終わったすぐ先が信号機であるから慎重な運転となるため
・見通しが良い
・歩行路と道路がフェンスで区切られているため、歩行者との...
***落橋について [#g86999e7]
落橋とは、橋梁の劣化や大規模な災害などが要因となって橋梁...
落橋の主な原因は橋梁の老朽化であり、損傷の種類は[[日本橋...
また、大規模な災害、具体的には地震や津波、土砂災害、強風...
落橋がもたらす社会的影響は、周囲の交通機関に不便をもたら...
***ケーブル腐食について [#o082b4da]
腐食とは、化学や生物的な作用によって物体の外観や機能が損...
金属の腐食を特に「錆」という。
腐食の主な要因は、塗装劣化個所に水分と腐食性物質が作用す...
腐食性物質としては、大気中の亜硫酸ガスと海塩粒子等が挙げ...
これらは地域や気候に関係が深く、工業地域や海岸近くに位置...
発生のメカニズムとしては金属は大気中の酸素によって表面に...
この酸素と水分、金属表面の汚染物などが表面に付着し酸化被...
ここから腐食が表面から内部へと進行していく。
金属にとって「錆」状態は安定している。
腐食の種類は均一腐食と局部腐食の２つである。
均一腐食…金属表面全面に発生する腐食。
局部腐食…局部的な腐食で発見しにくい。一般的な製造部品など...
ここから本題である「ケーブル腐食」を考える。
斜張橋と吊り橋は両橋ともケーブルを用いる形状をしているが...
斜張橋は、主塔から張られているケーブルが直接桁橋を支える...
また吊り橋はアンカーレッジが存在する。
斜張橋の強みとして、桁曲げモーメントを小さくできること・...
スパン長200～500mm及びそれを超えると斜張橋か吊り橋で設計...
それぞれの橋梁の桁に加わる力として、
斜張橋…垂直方向の張力＋橋軸方向の圧縮力
吊り橋…垂直方向の張力
がある。
ケーブルは通常「亜鉛メッキ銅線」で構成されているらしい。
亜鉛メッキとは高い耐食性（腐食のしにくさ）を付与出来るメ...
亜鉛金属を電気の力で析出させる方法と溶融させた亜鉛に浸漬...
亜鉛メッキは防錆効果が高く主に鉄製品に対して処理を行いま...
また、亜鉛は鉄よりイオン化傾向が大きいのでメッキ皮膜にピ...
これを犠牲防食という。
ただし、デメリットとして亜鉛メッキは湿気に弱い。
ケーブルが腐食する主要因は「水（＋塩）」であると考える。
特に、乾湿を繰り返す部位は、常に水に覆われた部位より腐食...
その理由は、乾湿を繰り返す部位が常に水に覆われた部位と比...
ケーブル腐食がもたらす影響をいくつか考える。
①斜張橋の終局強度の低下…ケーブル断面の減少＋伸び及び疲労...
②疲労寿命（破壊が起こるまで材料が耐えられる、ある種の変動...
③橋梁全体の耐力低下…ケーブル断面の減少 から
④ケーブルの伸び及びねじれ強度の著しい低下…亜鉛メッキが腐...
上記①～④の原因・理由を考えるために[[斜張橋ケーブルの耐久...
-ケーブル断面の欠損により、ケーブル鋼線が降伏するまでは斜...
-小段数ケーブル斜張橋の場合は、ケーブル破断の寸前になると...
-マルチケーブル斜張橋の場合は、ケーブル１本が破断しても特...
-応力腐食割れ（腐食環境下において金属材料に引張応力が作用...
ケーブルは主塔・主桁との協働作用により外力に抵抗する非常...
これにより腐食したケーブルで耐えている斜張橋や吊り橋は、...
橋梁の突発的な崩壊は絶対に避けるべき事故である。特に斜張...
従って、当たり前ではあるがケーブル破断が起こらないように...
--補足
風によるケーブルの振動現象は
①カルマン渦による渦励振
②六角形断面ケーブルにおけるギャロッピング
③並列ケーブルにおけるウェイクギャロッピング
④レインバイブレーション
⑤塔からの剥離渦によるウェイクレゾナンス
***リスクアセスメントについて [#zdb76e24]
[文献から…]
リスクアセスメントとは、危険が発生する可能性の度合いを組...
橋梁におけるリスクアセスメントは、橋梁がどのくらいの割合...
欧米ではリスクアセスメントが主流なため「何が現実に実行で...
日本ではこの考えはまだまだ定着しておらず、リスクゼロを前...
しかし実際にはリスクゼロは不可能であるため、日本でもリス...
[リスクマネジメント&リスクアナリシスの流れ]
基本条件・情報の整理→リスクの定量化・ランク付け→対策を実...
[ALARP：As Low as Reasonably Practicableの原則]
① 「受け入れ可能なリスク」として「自然災害の発生レベル以...
② リスク低減のためにどこまでも費用をかけることが合理的で...
③ 低減措置の実施基準は，リスク低減効果と低減に要する費用...
また安全に対する考え方で、，「広く受容される領域」・「我...
ALARP領域はそのリスクが許容されるわけではないが、ある便益...
許容可能なレベルに達していないリスクから優先的に対策を実...
事故の発生件数に重きを置くのではなく、重大事故の防止を徹...
リスクアセスメントのポイントは「事前にどれだけのリスクを...
日本では想定したリスクにすべて対策を講じるという考えが強...
この「すべて対策を講じる」という考え方が弊害となって日本...
リスクアセスメントで大切な要素は以下の２つであると考える。
①我慢できる領域と受け入れられない領域の線引きをしっかりと...
②受け入れ不可能な領域のリスクに対しては優先的に対策を実施...
すなわち、これは諦める・これは絶対に阻止するという取捨選...
↑この考え方は日本人には受け入れずらいため、なかなかリスク...
リスクアセスメントは重大なリスクに対処するための合理的な...
また、過去に起こった良い経験も悪い経験も同等にフィードバ...
具体的には設計の初期段階で、
＝通常の使用＋事故＋誤使用＋早期劣化＋etc..
ただ、懸念すべきことはリスクアセスメントを考えれば考える...
起きるかどうかわからないことにお金や時間をかけることは、...
これもリスクアセスメントの実現が難しい理由に当たると思う。
***IHIの方の講座(5/31) [#mb17c552]
現代では建設において、ＣＯ２の削減が優先である.
ＣＯ２の排出量はアメリカは横ばい、ＥＵは下がっている.
今の日本では【１橋梁/104人】を支えている.
土木の現場では、ＣＯ２排出量は コンクリート＞鉄筋＞桁 ...
ＣＯ２を減らす具体例としては、耐候性を用いること・木材を...
ヨーロッパでは、工事において重要視されているのは費用より...
しかし日本は人口が多いため、現時点ではＣＯ２を考えること...
今の時代に沿った新しい優先すべきことが出てきたと考える.（...
***橋梁と火災の関係 [#d0651243]
鋼橋にまつわる火災事故は、国内外問わず数多く報告されてい...
しかし実際に発生することは稀であるため、その対策をしてい...
よって、いざ火災事故が発生してしまった場合に、利用者の安...
そこで、土木学会鋼構造委員会では，2012年6月に「火災を受け...
ここでは、火災による熱の影響を受けた材料の特性はどのよう...
供用中及び工事中に考えられる火災原因を列挙する。
供用中…不審火・走行車両（タンクローリー等）の横転及び炎上...
工事中…動力工事・溶断及び溶接による火の粉の飛散・塗料作業...
材料自体が出火原因となる事例はほとんどない。
また。延焼に必要不可欠なものは①火源、②燃えるもの、③酸素の...
延焼要因は、放射熱・接炎・熱気流・火の粉などがある。
加えて炎が風下側に傾き隣接した建物に直接接炎する場合もあ...
鋼橋は火災の被害を受けると、構造部材は温度上昇に伴う強度...
一般的に鋼材（構造用鋼材）は、受熱時の温度が600℃を超える...
鋼材は合熱処理によって結晶構造が変化するという性質を有し...
よって、火災による鋼材の損傷評価については受熱温度の推定...
このときの鋼橋の補修ですべきことは、塗装の補修・鋼桁の補...
火災発生から補修施工までの流れ👇
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ka...
引用元：[[火災を受けた橋梁の安全性評価と補修＿実例:https:...
主塔は400℃、ケーブルは700℃くらいまで耐えられるそう。
SS材は通常〜300℃くらいまで耐えられる。
SM材は400℃〜500℃程耐えられる。このとき、500℃まで耐えられ...
火災による影響はケーブルよりも主塔にある。
そのため主塔が受熱する温度や時間によって、耐力や強度、変...
例）主塔が崩壊したときのケーブルの切れ方・主塔の耐力が落...
健全度評価フロー👇
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/fr...
鋼材部の被災目安👇
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ko...
引用元：[[火災を受けた橋梁の健全性評価に関する試験調査:ht...
鋼材の温度が高くなると、柔らかくなるすなわちヤング率は小...
火災による鋼材の変質は目視試験での確認は困難である。
また、実鋼材を抜き取り機械試験などにより確認することは通...
火災の被害を受けた鋼部材における補修を行うべき損傷は主に...
車両が燃える原因として、燃料やオイルなど着火物が漏れたり...
***運転者の心理について [#g0482488]
「人」が原因となって起こる交通事故＝”ヒューマンエラー”
エラーの種類として、思い込み・知識や技能不足・意図的ルー...
また、急いでいるときには必然的に事故を起こしやすい.
疲労・ストレスが交通事故へつながることも多々.
事故を起こした直後の運転者は二次被害につなげないように、...
その行動が橋梁上では、主塔やケーブルに近づくことを意味する.
もし、追突や衝突で車両が爆発した上でこの行動を取ったとし...
***交通事故の原因について [#rdbc3731]
令和４年における交通事故の原因ランキングを[[交通事故の原...
１位 安全不確認
２位 脇見運転
３位 動静不注意（他の車や歩行者の存在には気づいているも...
４位 漠然運転（ぼんやりと運転をしている様）
５位 交差点安全進行
６位 運転操作不適（危険は認識しているものの、それに対す...
７位 一時不停止
８位 信号無視
９位 歩行者妨害等（横断歩道を渡っている、または渡ろうと...
１０位 優先通行妨害（信号機がない交差点で、標識や道路幅...
近年は交通事故数が減少している.
これの主な要因を内閣府は以下のように示している.
１，道路交通環境の整備
２，交通安全思想の普及徹底
３，安全運転の確保、車両の安全性の確保
４，道路交通秩序の維持
５，救助・救急体制等の整備等
６，交通安全基本計画に基づく諸対策を総合的に推進してきた...
単路や橋での事故として多いものは[[様々な交通事故統計:http...
１，車両相互追突
２，車両相互正面衝突
３，車両相互その他
４，車両単独工作物衝突
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ch...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ca...
***リスクアセスメントの質問 [#af354ae0]
・リスクの確率の計算方法（事故の確率や火災の確率）及び計...
・死者数や被害者数の想定方法（人口に対して何割か）
・火災の起こる確率の式は何を使うか
・主塔への熱の伝達方法（動的解析か静的解析か）※今は静的解...
→毎秒もしくは毎分、何度熱が加わりそれがどのくらい続...
→ある特定の温度のときに主塔がどのように変形し、かつ...
・由利本荘市や市民に対する経済への影響や交通への影響を考...
→橋梁火災の規模やその発生原因、被害者数、被害状況、...
***鋼構造を語る会（9/20） [#mf556213]
今日は鋼構造を語る会に参加した。
他の大学の人の研究を聞く機会がなかったため、とても面白か...
現在、二酸化炭素を排出しない「グリーンスチール」というも...
最後のゲームでは３位になった！嬉しい！
***青木ゼミ(9/24) [#k361f058]
今日は青木班のみんなと話し合いをした。
【今後の方向性】
デッキ（死荷重・活荷重）半分をつけて同様の解析
火災が起きる→この程度の事故率や火災率が考えられる→火災で...
【リスクアセスメントの浸透】
安全対策への考え方の違い
日本では事故を起こすことが悪だが、海外では事故が起こった...
***中間発表(10/3) [#w4d8d623]
夏休みが終わり中間発表を行った．
スライドと概要に思ったよりも時間がかかったため，もう少し...
発表は時間内に終わらなかったため，もっと練習が必要である．
やっとリスクアセスメントについて皆にまとめたものを見せる...
【質問】
・リスクアセスメントの優先順位等はあるのか（ケーブルが重...
・ケーブル腐食などもリスクアセスメントに含まれるのか
・由利橋に対してリスクアセスメントをしたいのか由利橋を通...
・プロパンガスからの火災は爆発ではないのか
***IHIの方とリスクアセスメントもついて話し合い(10/28) [#o...
もっと「定量的」にリスクを設定していくことが必要！！
そのために・・・
世界や日本での事故率・危険車両の通過割合といった必要なデ...
リスクアセスメントは、過去に起こったデータを用いて予想さ...
適当な過程を自分でつけていくことが必要となる
解析では、由利橋が通行止めになった場合にどの程度困るのか...
ケーブルが火災にあった場合もやってみたいかも。でも、ヤン...
最終的な目標は、そのリスクが許容できるor許容できないかを...
許容できない場合は許容できるまでリスクに対してお金をかけ...
耐火塗装や消化器の設置などなど
わざわざ斜張橋にしたことによって新たなリスクが増えたので...
メンテナンス費用の増大や火災が起こりやすくなった、主塔か...
将来的には一般橋梁と斜張橋との比較とかもいいかも
一般橋梁と由利橋とではリスクが起こる確率が高いor低い
***yec香港(11/30〜12/1) [#zf3ed9cd]
貴重な経験となった４日間になった
色んな人の発表や意見に触れてとても勉強になった
楽しい時間を過ごすことができた！
他のリスクアセスメントの研究をやっている人は実際にリスク...
多くの人は「今はこう」だから「研究結果からこの部分を改善...
質疑応答がやはり反省点。聞かれたことに対して上手く答えら...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/IA...
↑yecで使用したスライド
[質問されたこと]
・活荷重が大きすぎるのではないか
・リスクアセスメントで使用しているデータはどこから入手し...
・リスクは橋梁に対応しているものか
・なぜ直線道路での事故率を考慮する必要があるのか
***今後の方針(12/11) [#qfa15d48]
本研究では「通行止めリスク」を考える
その要因を火災とし、小型車・大型車・特殊車で火災の特徴（...
火災解析は主塔正面・側面で起こったときとケーブルC1,C6で起...
[通行止めリスク]は以下の通りにする
・橋梁上の事故率
・秋田県車両火災率
・大型車混入率
・主塔周りで火災が起こる確率 ←自分で考える
・ケーブル周りで火災が起こる確率 ←自分で考える
データは過去３〜５年分集める
***中間発表(12/23) [#q013a558]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/70...
上記が使用したスライド
[質問されたこと]
・なぜスパン長をかけるのかがわからない
・もっと秋田県の事故率や凍結による事故率など秋田により沿...
・使用した車両火災率は、事故によるものなのか
・なぜ火災がケーブルに影響しなかったのか、何℃くらいになっ...
***卒論発表練習(2/6) [#j5e8bb22]
概ねうまくできた．
話したいことがたくさんで時間が足りない．
もっと要点をまとめた方がいいかもしれない．
最後のまとめは，もう少し全体的なまとめと具体例を提案した...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/02...
***卒論発表会(2/18～2/19) [#c8e0d4c4]
卒論発表が終わった．
練習では概ね上手くできていたと思っていたが，質疑応答は上...
でも，改めて勉強になることが多く，気づきもたくさんあった．
この時期に卒論発表をすることができてよかった思う．
今後はこれからの目標をより明確にする必要があると感じた．...
加えてやはり許容上限と下限の算出方法・そのリスクを選定し...
今はリスクアセスメントの中身（リスクイベント・リスクを計...
今後はもっとおおもとの枠組みを把握する必要がある
例えば，ユーロコードの内容・許容値の設定方法・実際にどの...
日本の通行止め日数やユーザーリスクを調べてみる，
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/so...
*サロメモ [#v109bbb0]
AsterStudyの設定翻訳
,メッシュ,Read a mesh, メッシュの作成をしたモデルを適用し...
,Model Definition,Assign finite element, モデルの要素を設...
,Material ,Definr a material ,各材料物性値の設定,(複数設...
,Material, Assign a material ,設定した物性値を材料に適用,
,BC and load , Assign mechabucal load , 拘束条件 荷重設...
,Analysis , static mechabucal analysis ,線形解析の設定 ,...
,Post Processing , CALC_CHAMP ,ここは出力でなんの要素を出...
,out put , set output result ,出力ファイルの設定を行う。,
[メッシュ]
>LIRE_MAILLAGE(メッシュを読む)>FORNAT()
>ASSE_MAILLAGE(メッシュ組み立て)
>DEFI_GROUP()
>MODI_MAILLAGE(メッシュ方法)
>CREA_MAILLAGE(メッシュ決め)
[Model Definition]
>AFFE_MODELE(モデル割当)>AFFE(割当)>TOUT(すべて)>PHENOMEN...
>AFFE_CARA_ELEM()
>DEFI_BASE_REDUITE()
>DEFI_DOMAINE_REDUIT()
>DEFI_GEOM_FIBRE()
>MACR_CARA_POUTRE()
>MODI_MODELE()
[Material]
>AFFE_MATERIAU(材料割当)>MODELE>AFFE>TOUT(すべて) ←AFFEで...
>AFFE_MATERIAU>AFFE_VARC()>NOM_VARC()→TEMP(温度)>VALE_REF...
>CREA_LIB_MFRONT
>DEFI_COMPOSITE
>DEFI_GLRC>
>DEFI_MATER_GC
>DEFI_MATERIAU(材料特性)>ELAS(ヤング率等)>E(ヤング率),NU(...
>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(線載荷)orFOR...
>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(構成...
NEWTON(ニュートン法の設定)>MATRICE:TANGENTE(接線剛性)
CONVERGENCE(収束判定の設定)>RESI_GLOB_RELA(相対誤差)>ITER...
>Thermics(熱力学)
>THER_LINEAIRE(温度＿線形)
>THER_NON_LINE(非線形)
>THER_NON_LINE_MO(非線形モデル？モード？)
[Post Processing]
>CALC_CHAMP(調べたい分野)>CONTRAINTE(対照？)>SIGM(6成分応...
>CALC_CHAMP(調べたい分野)>DEFORMATION(力？)>DEGE(),EPME()...
>CALC_CHAMP(調べたい分野)>CRITERES(基準)>EPEQ(),EPGQ(),EP...
NOEUは節点(NOEUd(フランス語))、ELGAはガウスの積分点(ELeme...
>RESULTAT>MECA_STATIQUE
>CREA_CHAMP>MODEL→AFFE_MODEL>TYPE_CHAM→NOEU_TEMP_R>OPERAT...
[Output]
>IMPR_RESU(結果照会)>RESU(結果)>RESULTAT(結果)>NOM_CHAM(...
*Salome-Meca2020_熱解析 [#je8a4800]
○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:tem...
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...
FORMAT>MED
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/re...
Astersutudyはどこかで設定してから出ないと選択できないもの...
○Paravis
左側のパイプのブラウザーのbuiltinを右クリック>Open>.medフ...
フィルター(画面上の黒いところ、ファイルとか編集とかあると...
フィルター>Recent>ExtractGroup>SIGM_NOEU(節点の応力(6成分))
*Salome-Meca2020_弾塑性解析 [#z08773c9]
参考[[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>SIM...
NEWTON>MATRICE(マトリクス)>TANGENTE(タンジェント)orELASTI...
CONVERGENCE>RESI_GLOB_RELA>"1.0e-06"任意で良い
CONVERGENCE>ITER_GLOB_MAXI>"100"任意で良い
*Salome-Meca2020_弾塑性解析2 [#qd419800]
○Geometry
いつも通り作成
○Mesh
いつも通り作成
○AsterStudy
右クリックAddstage
[[Function and Lists]]
-DEFI_LIST_REEL(時間間隔の設定)
DEBUT>"0">
INTERVALLE>"1.0","Step length","0.05"入力
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/li...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/in...
-DEFI_LIST_INST(解析オプション)
DEFI_LIST>LIST_INST:listr>
METHODE>MANUEL
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ti...
-DEFI_FONCTION
名前:func>VALE>"0","0"と"1","1"入力>
NOM_PARA>INST
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/fu...
-DEFI_FONCTION
名前:koufuku>VALE>"390","400"と"235","236"入力>
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ko...
NOM_PARA>TEMP>
NOM_RESU>SIGM
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ko...
[[Material]]
-DEFI_MATERIAU
ECRO_LINE>D_SIGM_EPSI:200入力>SY:235
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ec...
-ELAS
E>NU>RHO>ALPHA
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/el...
[[Analysis]]
-STAT_NON_LINE
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/re...
MODELE:model>
CHAM_MATER:fieldmat>
EXCIT>CHARGE:lode>FONC_MULT:func>
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ex...
INCREMENT:times>
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/in...
COMPORTEMENT>TOUT>RELATION>VMIS_ISOT_TRAC>DEFORMATION>PET...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/co...
NEWTON>MATRICE(マトリクス)>TANGENTE(タンジェント)orELASTI...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ne...
CONVERGENCE>RESI_GLOB_RELA>"0.0001"任意で良い
CONVERGENCE>ITER_GLOB_MAXI>"1000"任意で良い
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/co...
*Salome-Meca2020_モデル作成（Fuseする場合） [#l44a64aa]
○Geometry
1, BOX-1,2,3・・作成
2, Fuse-3objects選択
3, Partitionを作成、このときFuse/BOX-1,2,3を選択
4, Partitionのグループ作成、kotei/kanetuなど
5, Meshを作成
6, ジオメトリのグループ作成
以上
*研究メモ [#i24c5d3a]
-Windowsからstr/k2サーバーへは”FFFTP”を使用する.このとき...
-要素の種類…ビーム＝線、シェル＝面（厚さなし）、ソリッド...
-モデルは半分で良い
-主桁と主塔は節点は重なっても良いが、境界条件は別々で考え...
-マイダスの境界条件設定は「単一選択」で行う
-マイダスでの由利橋モデルにおける等分布荷重解析を行う際は...
-主塔水平材とは、両主塔を固定するための部材であり幅員に対...
地震時のみに機能を発揮し、２本の主塔が個別振動することを...
また、軽微な損傷を許容することで主塔柱の発生断面力を制御...
-プレテンション(引張力＋)とは、ワイヤロープに一定荷重を一...
初期伸びとロープ径の減少が少なくなる・伸びが少なくなる(弾...
-曲げ剛性とは部材の曲げやすさを表す値である.曲げ剛性（Ｅ...
また、曲げ剛性(EI)と曲率半径(r)は比例関係にあり、曲げモー...
-部材内部に生じる応力＝軸力・せん断力・曲げモーメント
-軸力…部材の材軸方向に生じ、材を伸縮させようとする一対の...
-せん断力…部材に「ずれ」の変形を生じさせようとする一対の...
-曲げモーメント…部材を曲げようとする一対のモーメントのこ...
-合成桁とは、床版と鋼桁とを"ずれ止め"で接合し合成床版と同...
-普通鋼とは、炭素（C）・ケイ素（Si）・マンガン（Mn）・リ...
-一方で、SM材はSはSteel（鋼）、Mは船舶を意味するMarineを...
--熱膨張係数とは、線膨張係数αもしくは体膨張係数βから求め...
--主桁の数え方は左からG1,G2…とする
--マイダスで荷重組み合わせの生成方法は、[[マイダス〜解析...
の９－１９ページに記載
--解析結果の変位などを見るやり方は、結果テーブル＞見たい...
--マイダスで変位等をエクセルに書き出しができるため、そこ...
これが写真→→→http://www.str.ce.akita-u.ac.jp/~gotouhan/j2...
--温度解析はサロメでやっても良い
--由利橋モデルは、桁・主塔は梁要素、ケーブルはトラス要素...
--梁要素とは、有限要素の一種であり、2点を結ぶ細長い軸形状...
--トラス要素とは、圧縮／引張の軸力のみを伝える部材のこと
--弾塑性とは、降伏応力以下の力を加えた場合は力を除くと元...
--マイダスでは３Ｄモデルは”シェル”が該当する
--粘性係数とは、流体の粘り気のこと、鋼材の動粘度ν＝0.5e-2...
--応力-ひずみ線図 http://www.str.ce.akita-u.ac.jp/...
--座屈とは、ある一定の荷重を超えるとき柱は急激に折れ曲が...
--主要地方道は高速道路及び一般国道と一体になって広域交通...
--都市ガスは地下のガスパイプ通って各家庭へ
--街道は一般国道の分類（羽州浜街道：うしゅうはまかいどう）
--12時間交通量は午前7時〜午後7時までの交通量。24時間交通...
--軽トラックは車高1.95m、積載高さ2.5m以下という規定。
--交通量が多い目安：4000台/日〜、大型車混入率目安：地方15％
--大型車：バス・トラック・ダンプカー・ミキサー車など
--小型貨物車：軽トラック・小型トラックなど（乗員定員10人...
--タンクローリーとは、ガソリンや水といった液体物や高圧ガ...
--圧力荷重とは，、面上の特定の種類の静荷重のこと。 面上の...
--単路と直線は、道路の形状を指す言葉で、単路は交差点付近...
--一般単路とは、交差点や踏切、一般交通の場所、トンネル、...
*研究参考文献 [#b552c166]
-[[色々なヤング率:http://www.mh.rgr.jp/memo/mq0100.htm]]
-[[断面二次モーメント計算ソフト:https://d-engineer.com/un...
-[[マイダス〜構造編〜:https://in.midasit.co.jp/wp-content...
-[[マイダス〜解析編〜:https://in.midasit.co.jp/wp-content...
-[[マイダスマニュアル:https://down.midasuser.com/UploadFi...
-[[ケーブルの基本理論:http://library.jsce.or.jp/Image_DB/...
-[[由利本荘市ハザードマップ:https://www.city.yurihonjo.lg...
-[[日本橋梁建設協会:https://www.jasbc.or.jp/panfuretto/pa...
-[[鋼橋の補修・補強事例:https://www.jasbc.or.jp/whatsnew/...
-[[鋼橋（上部工）の損傷事例:https://www.mlit.go.jp/road/s...
-[[鋼橋の腐食事例:https://www.jstage.jst.go.jp/article/js...
-[[斜張橋ケーブル腐食を考慮した疲労寿命:https://www.str.c...
-[[斜張橋ケーブルとアンカーについて:https://www.jpci.or.j...
-[[斜張橋ケーブルの耐久性評価:https://kandoken.jp/huroku/...
-[[橋梁の構造＆種類:https://chansato.com/doboku/bridge-st...
-[[鋼橋工事における火災リスクと防火対策:https://www.nilim...
-[[火災時における耐荷力に関する検討:http://library.jsce.o...
-[[火災を受けた橋梁の損傷と対応:https://www.cbr.mlit.go.j...
-[[火災を受けた橋梁の安全性評価と補修＿実例:https://www.h...
-[[火災を受けた橋梁の健全性評価に関する試験調査:https://w...
-[[鋼材の高温時における機械的特性:https://www.jstage.jst....
-[[SM400材について:http://www.ishiharashouji.jp/pdf/QA-Ve...
-[[普通鋼について:http://www.iri.pref.niigata.jp/25new32....
-[[SM材の耐熱能力:https://higa-metal.com/articles/793/]]
-[[橋梁用厚鋼板:https://www.jfe-steel.co.jp/products/atui...
-[[鋼材の材料定数:https://www.structure.jp/databook/data1...
-[[交通事故の原因ランキング:https://agoora.co.jp/jiko/kno...
-[[様々な交通事故統計:https://www.ms-ins.com/kurumamo/ent...
-[[弾性変形と塑性変形:https://www.monodukuri.com/gihou/ar...
-[[トレス・ミーゼス応力の降伏条件:https://mech-hm.eng.hok...
-[[トレス・ミーゼスの違い:https://planteng-tree.com/yield...
-[[ミーゼス応力の定義:https://d-monoweb.com/blog/evaluate...
-[[モール・クーロンの破壊基準:https://www.jseg.or.jp/chus...
-[[面外と面内:https://bonperson-civil.com/in-out-plane/]]
-[[様々な車高:https://suumo.jp/article/oyakudachi/oyaku/i...
-[[鋼材の粘度:https://www.daimaru-kogyo.co.jp/skillinfo.p...
-[[ヒルの降伏応力:https://www.jsme.or.jp/jsme-medwiki/dok...
-[[円筒座標系:https://www.cradle.co.jp/glossary/ja_A/deta...
-[[直交座標系:https://www.cradle.co.jp/glossary/ja_T/deta...
-[[テンソルとは:https://3.1415.jp/35ss9vvx]]
-[[火災を受けた橋梁の安全性評価:https://www.oit.ac.jp/jap...
-[[火災を受けた鋼橋の受熱温度推定に関する解析的研究:https...
-[[火災を受けた鋼橋の損傷調査とその対応について:https://w...
-[[熱分析:https://www.spsj.or.jp/equipment/news/news_deta...
-[[各材料のメリット・デメリット:https://skyciv.com/ja/tec...
-[[斜張橋の設計:http://library.jsce.or.jp/Image_DB/commit...
-[[SM400Bの熱による耐力:https://jopss.jaea.go.jp/pdfdata/...
-[[鉄鋼材料の各温度における許容引張応力:https://www.meti....
-[[高強度鋼:https://www.jisf.or.jp/business/tech/bridge/h...
-[[耐火性能:https://www.jisf.or.jp/business/tech/school/c...
-[[橋梁設計における手続きと照査:https://www.pref.aichi.jp...
-[[タンクローリー車のサイズ:https://driversjob.jp/content...
-[[橋の設計の基本と基準の変遷:https://www.nilim.go.jp/lab...
-[[海外における耐震設計基準と橋梁の比較設計:https://www.j...
-[[橋梁火災事例:https://www.nilim.go.jp/lab/bcg/siryou/tn...
-[[火災時における鋼合成桁の終局耐力と耐火対策に関する研究...
-[[平面応力要素とは:http://www.cae21.org/kansai_cae/Qa/Q0...
-[[板要素を利用したFEM解析:https://in.midasit.co.jp/wp-co...
-[[midasフリーズ解決方法:https://tech.midasuser.com/jp/su...
-[[midas自重:http://manual.midasuser.com/JP/iGen/835/Usin...
-[[移動荷重:https://support.midasuser.com/hc/ja/articles/...
-[[midas特殊解析:https://in.midasit.co.jp/wp-content/uplo...
-[[道路橋示方書:https://www.mlit.go.jp/road/sign/kijyun/p...
*サロメページ [#oda4d11a]
-[[Salome-Meca演習_弾塑性解析(2021)]]
-[[Salome-Meca_熱応力メモ]]
-[[Salome-Meca2017]]
-[[Salome-Meca例題ファイル]]
-[[Salomeメモ]]←ここから線膨張係数入れられそう！
-[[座屈解析]]
*春休み課題の進捗 [#nd67ffe7]
-一端固定他端単純支持梁
--モデル作成について
ローラー支持は「線固定」→線の端点を入力して、線...
固定端は「面固定」→パーティションからグループの...
Mesh→ジオメトリは「Partition-1」を必ず選択
Mesh切る→ジオメトリのグループ作成
☆☆AsterStudyは簡単にできる！！→Operations→Addstage ...
☆☆（E:206000,ν:0.3,ρ:7.8e-9）
☆☆ファイル名には「.med」！！
--パラビィスのやり方
builtinでopen→「Mode」に変更
フィルター→「Mechanics」から「Normal modes anim...
**やること [#i3d95498]
--理論値を求める・・済（ねじれは分からない）←分かった！！
--メッシュ1.0のharusindou_1の解析結果を検討・まとめ・・済
--上記の相対誤差を求める・・済
--メッシュ0.8,0.9,2.0,3.0を解析→比較・検討・・済
→でも、全ての結果を比較しても ①要素数と相対...
→よって「梁形態」が間違えていると考える
☆解決策 ①梁形態を変更（ピンを端っこから10mm...
②ピンを２コに変更（上下...
③ローラー固定で考えてみ...
--どの梁形態を基本にするか決定する・・済 →☆☆「ローラー」
--上記の決定した一端固定他端単純支持梁のメッシュ0.7~4.0を...
--上記の相対誤差の計算・・済
--ローラー上下２個の一端固定他端単純支持梁のメッシュ0.7~4...
--ローラー下１個の一端固定他端単純支持梁のメッシュ1.0~10....
→！！メッシュ1.0...
--上記の解析結果からねじれについても検討・考察・・済
--ローラー下１個の100mmの一端固定他端単純支持梁のメッシュ...
--ねじれをどうするか話し合う→とりあえず後回し・・解決
--グラフ作成（９枚）・・済
--スライド作成・・済
--グラフの確認・使うデータ選別・・済
--スライド作成（長さ・振動の様子の２種類）・・済
--収束値との相対誤差-要素数のグラフ作成→細長比で比較して...
**疑問・考察 [#zaa418b3]
-２次要素にすると鉛直・水平・ねじれモードは収束した。１次...
-１次要素だと水平モードの方が精度が良いが、２次要素だと鉛...
-１次要素だと高次数の方が精度が良いが、２次要素だと低次数...
-細長いとねじれが出にくい。→100mmだと固有振動数１０コでね...
-長さのみ・長さ＋幅・長さ＋高さを変えた条件でそれぞれ解析...
→やはり、１次要素だとどれも水平モードの方が精度が良かった。
→鉛直・水平・ねじれを総合的に見ると長さ＋幅が一番いいかも...
春課題を通して、理論値への収束の条件や様子、理論値に対す...
既に答えの出ている課題であると説明を受けていたが、どれも...
また、１つの知りたいことのために解析や考察を行ってもまた...
ただ結果を出すだけでなく、そこから何が考えられるか、何が...
*３年創造工房 [#d09b48c4]
**2023.10.27 [#t59452ec]
・タッチタイプの練習
・コマンドの練習
-コマンドの一覧 [#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)...
**2023.11.10 [#uc3066e7]
・コマンドの練習
・gnuplotを使用して.pngで図を作成
-作成した図 [#kf32a2c4]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/bb...
使用データ：
,メッシュ長さ,要素数,変位(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,4518,0.3986,-4.34,山本
,1.5,15433,0.396317756757,-4.9,進藤
,1.6,15900,0.399049,-4.24,河合
,1.8,11677,0.404457,-2.9,山口
,2,10460,0.394818715517,-5.3,進藤
,3,734,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,河合
引用元：岡田の卒論日誌（11/18課題）
**2023.11.17 [#k6d60e45]
・片持ち梁の解析（サロメ）
-サロメメモ [#ddb92132]
Geomety BoXを作成 適用→閉じる kotei/saika
Mesh メッシュ作成 Box1を挿入 アルゴリズム Note ジオメ...
FY=100N/100mm^2=1N/mm^2
・測定値＝平均値ー理論値/理論値✕100(％)
・要素数＝ボリューム
・変位平均値（1.3.5.7）
・相対誤差（測定値）を計算
・縦軸＿変位、横軸＿要素数のグラフを作成
--片持ち測定値 [#xcffb3ee]
,メッシュ長さ,要素数,先端変位(4隅の平均値)[mm],相対誤差($...
,0.7,198464,6.54281,1.91,安藤
,0.8,113812,6.5104,2.39,安藤
,0.9,40280,6.3631525,4.60,兼田
,1.1,30055,6.3363525,5.00,兼田
,1.2,26467,6.3043375,5.48,柴田
,1.3,25180,6.304355,5.48,柴田
,1.4,32212,6.31612,5.31,佐藤
,1.5,17753,6.1209,8.23,佐藤
,1.6,14296,6.2044625,6.98,皆川
,1.7,13596,6.2156625,6.81,皆川
,1.8,2866,5.737755,13.98,永山
,1.9,6001,5.7263625,14.15,永山
,2,5617,5.6458525,15.355,辻
,3,2309,5.4728755,17.948,辻
,4,617,3.6160575,0.458,服部
,5,494,3.8580375,0.422,服部
,6,581,2.50682,62.416,梶原
,7,133,1.41225,78.827,梶原
,8,78,1.2887175,80.68,工藤
,9,72,1.2879925,80.69,工藤
,10,60,1.14344,82.85,佐々木
,11,65,1.23124,81.154,佐々木
--グラフ（理論値と測定値） [#y04fdb0d]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/yo...
縦軸：変位(mm)、横軸：要素数
理論値：6.67mm
【理論式】
Pl^3/3EI
**2023.11.24 [#f5b0e5b1]
・単純梁（等方性１次）の解析
--単純梁（等方性１次）測定値 [#g80425ab]
,メッシュ長さ,要素数,変位平均[mm],相対誤差($\frac{salome-...
,0.7,145234,0.42248452735,1.388176,安藤
,0.8,142973,0.42257044598,1.408794,安藤
,0.9,91648,0.420437286573,0.897,兼田
,1.1,27160,0.405618939024,2.659,兼田
,1.2,24675,0.404349,2.96,柴田
,1.3,23446,0.404185,3.00,柴田
,1.4,17738,0.398604,4.34,佐藤
,1.5,15438,0.396593,4.83,佐藤
,1.6,16122,0.398212,4.44,皆川
,1.7,12026,0.393411,5.59,皆川
,1.8,11604,0.393668,5.53,永山
,1.9,10391,0.390695,6.24,永山
,2,10921,0.395103,5.18,辻
,3,2328,0.324762,22.06,辻
,4,1500,0.155013,62.80,服部
,5,432,0.065278,84.33,服部
,6,357,0.213062,48.87,梶原
,7,196,0.1019,75.55,梶原
,8,104,0.1158624,72.20,工藤
,9,81,0.1255118,69.88,工藤
,10,78,0.07733,81.44,佐々木
,11,63,0.1999,52.03,佐々木
--グラフ（理論値と測定値） [#cb6b3f6d]
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/1...
縦軸：変位(mm)、横軸：要素数
理論値：0.4167mm
【理論式】
Pl^3/48EI
**2023.11.29 [#g51a87b6]
・単純梁（異方性１次・等方性２次）の解析
＜等方性１次と異方性１次の比較＞ [#he11b38c]
--単純梁（異方性１次）の測定値
,メッシュ長さ,要素数,先端変位,相対誤差,計算者
,0.7,144563,0.505252,2.76,安藤
,0.8,141517,0.504692,2.64,安藤
,0.9,91648,0.502595,2.216,兼田
,1.1,27160,0.489914,0.363,兼田
,1.2,24675,0.487088,0.791,柴田
,1.3,23446,0.4868010,0.995,柴田
,1.4,17738,0.485999,1.16,佐藤
,1.5,15438,0.485180,1.33,佐藤
,1.6,15900,0.483286,1.71,皆川
,1.7,12142,0.477952,2.80,皆川
,1.8,11604,0.482085,1.97,永山
,1.9,10391,0.470887,2.40,永山
,2,10291,0.480910,2.19,辻
,3,2328,0.431937,12.15,辻
,4,1500,0.430156,12.52,服部
,5,432,0.282968,42.45,服部
,6,356,0.3441556,30.00,梶原
,7,196,0.213934,56.49,梶原
,8,104,0.229874,53.25,工藤
,9,81,0.232308,52.75,工藤
,10,78,0.203271,58.65,佐々木
,11,63,0.222316,54.78,佐々木
--比較のグラフ
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/1...
縦軸：変位(mm)、横軸：要素数
理論値：0.4917mm
【理論式：ティモシェンコの理論】
v = Pl^3/48EI ＋ Pl/4KGA
G：せん断弾性係数
K：せん断補正係数
このとき、［Pl/4KGA］= せん断項
＜等方性１次と等方性２次の比較＞ [#of7eeadc]
--単純梁（等方性２次）の測定値
,メッシュ長さ,要素数,先端変位,相対誤差,計算者
,0.7,144563,0.430124,3.22,安藤
,0.8,141517,0.430132,3.22,安藤
,0.9,91648,0.430020,3.197,兼田
,1.1,27160,0.429828,3.151,兼田
,1.2,24675,0.429836,3.15,柴田
,1.3,23446,0.42974,3.13,柴田
,1.4,17738,0.429797,1.3,佐藤
,1.5,15438,0.429958,3.14,佐藤
,1.6,15900,0.429755,3.18,皆川
,1.7,12142,0.429676,3.11,皆川
,1.8,11604,0.429829,3.14,永山
,1.9,10391,0.429684,3.12,永山
,2,10291,0.429620,3.10,辻
,3,2328,0.429169,2.99,辻
,4,1500,0.429254,3.01,服部
,5,432,0.428170,2.75,服部
,6,356,0.428452,2.82,梶原
,7,196,0.42591,2.21,梶原
,8,104,0.426074,2.25,工藤
,9,81,0.425552,2.12,工藤
,10,78,0.488382,17.20,佐々木
,11,63,0.423972,9.0534,佐々木
--比較のグラフ
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/1...
縦軸：変位(mm)、横軸：要素数
理論値：0.4167mm
２次要素を含む。
--メモ
・等方性・・方向によって物体の物理的性質が異ならないこと。
・異方性・・方向によって物体の物理的性質が異なること。
**2023.12.08 [#qa2a664e]
・単純梁(木材＋鋼材)の解析
--単純梁（木材＋鋼材）の測定値 [#d0d51f4f]
,メッシュの長さ ,要素数 ,変位[mm],相対誤差,計算者
,0.7,155192,0.08378905246,15.365,安藤
,0.8,138808,0.08380386491,15.350,安藤
,0.9,82587,0.083707073981,15.45,兼田
,1.1,38671,0.084201207602,14.95,兼田
,1.2,31929,0.083688,15.466,柴田
,1.3,28621,0.083669,15.4857,柴田
,1.4,28854,0.08368,15.47,佐藤
,1.5,20015,0.084052,15.10,佐藤
,1.6,19448,0.0835402938,15.62,皆川
,1.7,13801,0.0834355098,15.72,皆川
,1.8,12528,0.083733,15.42,永山
,1.9,11769,0.083924,15.23,永山
,2,10699,0.084076876559,15.074,辻
,3,3579,0.08414561753,15.004,辻
,4,1628,0.082794,16.37,服部
,5,1016,0.083033,18.89,服部
,6,839,0.082882,16.26,梶原
,7,554,0.080871,18.28,梶原
,8,285,0.079995,19.20,工藤
,9,261,0.078980,20.22,工藤
,10,232,0.081911,17.26,佐々木
,11,208,0.075676,23.56,佐々木
--グラフ（理論値と測定値） [#k36a4b43]
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/1...
縦軸：変位(mm)、横軸：要素数
理論値：0.099mm
２次要素を含む。
【理論式】
Pl^3/48EI ＋ Pl/4KGA
このとき、
EI = EI（木） ＋ EI（鋼）
G = G（木）
と考えることとする。
**2023.12.15 [#p21326f8]
・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を右クリックでジオメトリの作成
⑤ AsterStudyのmeshをダブルクリックして、UNITE「Mesh-2」...
⑥ Outputでファイルの名前変更（.med）
⑦ 解析
＜テキストエディターでグラフ作成＞ [#k2d05487]
① テキストエディターに（要素数,変位）の順で入力→このとき...
② 名前を付けてファイル保存
③ コマンド
ls→cd kaneta23→ls→gnuplot→＞plot "ファイル名" with line,"...
④ メニューから、gftpに行きパスワードを入力
⑤ public htlsの中のj2023から、自分のファイルに転送
終了行:
#contents
*研究 [#bf04b5a7]
**研究テーマ [#x8e2544a]
-ケーブルを有する橋梁のリスクアセスメントの提案ー秋田県由...
--実橋において，その橋がどのような危険にさらされているか...
--リスクアセスメントは橋梁がどのくらいの割合で危険にさら...
**やること [#me638c2a]
--参考文献を読む
--☑由利橋のモデルをリニューアルする〜5/31
--☑由利橋の事故・落橋の原因・落橋したときの不便・ケーブル...
--☑マイダスで３Ｄモデルを作成する…「ケーブル橋梁ウィザー...
--☑プレテンションの設定方法を模索…単純梁✖簡易モデルでプレ...
--☑マイダスで由利橋のケーブルにプレテンション導入
--☑由利本庄市の天気・交通事故を調べる〜5/31
--☑リスクアセスメントのまとめ〜5/31
--☑橋梁と火災について調べる
--☑事故のときの運転者の心理について調べる
--☑マイダスで火災や温度変化の解析ができるかどうかを試す
--☑由利橋モデルに正しいプレテンションと死荷重・活荷重を入...
--☑死荷重及び活荷重をどう入れるか聞く（集中or等分布？、そ...
--☑主塔をソリッド✖シェルで作成
--☑由利橋で事故が起こらない理由を考える
--☑サロメで温度解析のやり方学ぶ
--☑サロメで車高部分のみを加熱してどうなるか
--☑鋼材に粘性を入れてぐにゃってなるか片持ち梁で確認→弾塑...
--☑リスクアセスメントの質問を考える
--☑片持ち梁の片側側面に熱を与えてDX方向にどの程度変位が出...
--☑主塔の燃える部分のヤング率(0,400,800℃の３パターン)をグ...
--☑中間発表の概要作成・スライド作成 ー中間発表
--☑健全モデルの完成〜9/28
--☑リスクアセスメントを考慮する上で必要となる条件を考える...
--☑要素温度荷重解析
--☑IABSEのスライド作成
--☑リスクの修正
* 進捗 [#pb2c767f]
**由利橋モデルの作成(5/29) [#v4fab7c2]
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kane...
マイダスで由利橋のモデルを作成.
材料・断面・境界条件を設定、割当した.
スパン長190m,主塔70m（桁上50m,20m）
ケーブル断面積D=200mmのケーブル位置
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ca...
ケーブル断面積D=225mmのケーブル位置
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ca...
ケーブル断面積D=315mmのケーブル位置
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ca...
**マイダスで荷重解析(6/4) [#tbe94155]
幅1m×高さ0.8m×スパン長10mの単純梁を作成し、集中荷重（節点...
全体図（(x,y,z)空間）
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kane...
横から図（(y,z)平面）
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kane...
上から図（(x,z)平面）
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kane...
集中荷重（節点６に荷重FZ=-10kN）の写真
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kane...
等分布荷重（節点１２～２２に荷重w=-1kN）の写真
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kane...
以下が解析結果である、変位図・Ｓ図・Ｍ図である.
〇変位図
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kane...
〇Ｓ図
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kane...
〇Ｍ図
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kane...
**簡易モデルでプレテンション導入＿改訂版(6/7) [#ac0f5e55]
千代岡さんにアドバイスをいただき、軸力ではなく「温度荷重...
簡易モデルは前回のと同じものを利用.
温度荷重をケーブルに与えた.
簡易モデル
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ka...
〇ケーブル材料
→弾性係数:1.95e+02
→ポアソン比:0.3
→線膨張係数:1.2e-05
→単位体積重量:7.7e-08
〇ケーブル断面積
→D=1mm
〇桁材料
→SM400
〇桁断面積
→10mm×10mm
解析を回してみたところ👇このように変形した.
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ka...
軸力あり・なしバージョンで比較したが、軸力は関係なさそう…
また、先端荷重Fを与えたときの先端変位は下記の表のようにな...
このとき、ｙ軸下方向を＋とする.
,F(N),温度荷重(℃),先端変位(mm)
,0,-100,-2.378
,0,-150,-3.567
,0,-200,-4.756
,-100,-100,-0.543
,-100,-150,-1.732
,-100,-200,-2.921
,-500,-100,+6.796
,-500,-150,+5.607
,-500,-200,+4.418
表からわかる通り、荷重ありとなしのときの変異差はどの温度...
F=-100N…1.835mm
F=-500N…9.174mm
であった.
この結果から、プレテンションは温度荷重から与えることがで...
加えて、温度荷重が相対的に大きいほど先端変位も大きくなる...
温度荷重と先端変位は比例関係にあると言える.
**由利橋モデルにプレテンション導入＿改訂版(6/7～6/10) [#m...
軸力ではなく、温度荷重からプレテンションを入れる方法を由...
このとき全く軸力は入れていない.
下記の写真通りに、外側のケーブル２本に-150℃・その他のケー...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/on...
①なんの荷重もかけずに解析を回すと👇このように変形した.
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/yu...
max変位＝-0.394mm
②次に中央（節点９６）に集中荷重F=-190kNを与えた場合（下記...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/sy...
結果は👇このように変形した.
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/yu...
max変位＝+1.852mm
③次に全桁にw=-1kNを与えた場合（下記の写真）
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/to...
結果は👇このように変形した.
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/to...
max変位＝+1.008mm
この結果から、プレテンションをのみの場合は桁が引っ張られ...
また、等分布荷重よりも集中荷重の方が変位が大きくなった.
きちんとプレテンションが入っていると思われる.
**主塔に熱解析(6/13) [#f4cc604c]
主塔のみを取り出して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℃
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/sy...
初期20℃→解析後520℃
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/sy...
初期20℃→解析後1020℃
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/sy...
結果からわかる通り、加える熱が高ければ高いほど変位は大き...
＋100℃のときと＋1000℃のときの変位はちょうど１０倍になって...
本来であれば400℃～500℃で主塔は耐えられなくなるが、この解...
この解決策として「弾塑性」を設定に加えなければならないら...
**プレテンション荷重の入れ方(6/14～6/17) [#afe8c12d]
今回はプレテンションを温度荷重ではなく。プレテンション荷...
手順は以下の通りである.
荷重＞プレテンション荷重＞荷重ケース選択＞kNを入力＞単一...
プレテンション荷重以外には何も与えずに解析を回したところ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/pr...
max変位＝-73.49mm
ケーブル（トラス要素）の引張（断面力）はこのようになった.
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/pr...
ケーブル（トラス要素）の応力度はこのようになった.
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ca...
次に桁を見ていく.
桁（梁要素）の断面力はこのようになった.
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/pr...
これは断面図である.
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/pr...
桁（梁要素）の応力度はこのようになった.
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ha...
これでケーブルプレテンションは入っているのではないか・・？
これが合っているのであれば、温度荷重への変換は必要なし！
！！結論：プレテンションは「プレテンション荷重」の項目か...
**プレテンションを温度荷重に変換する方法(6/14～6/17) [#m8...
線膨張係数αの計算式である👇これを使う.
α = $\frac{ΔL}{L}\frac{1}{ΔT}$
ΔL/L…ひずみ
ΔL…変化量
L…元の長さ
今日は、先にプレテンションを設定して解析して、要素作成タ...
すると、Ｃ１ケーブルに与える想定温度1250℃、Ｃ２ケーブルに...
大きすぎな気がする…
！！結論：プレテンションを温度荷重に変換する必要はない.
**活荷重の入れ方(6/17〜6/19) [#u6540e70]
活荷重の計算は桁の半分を対象にしている.
① 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)のＬ荷重
を付与.
下👇の写真の通り.
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/L_...
活荷重のみを与えた変形は👇このようになった.
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/Ll...
max変位 ＝ +35.36mm
**死荷重の入れ方(6/17) [#lc6fb415]
死荷重の計算は桁の半分を対象にしている.
D = 77.0kN/m^3 × (9.5m+6.6m) × 1.1m /2 = 681.835kN/m
この死荷重を等分布荷重として桁に与えた.
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/de...
まだカウンターウェイト等は考えずにすべてを鋼材として計算...
死荷重のみを与えた変形は👇このようになった.
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/61...
max変位＝+955.91mm
**由利橋モデルにプレテンション＋死荷重＋活荷重を組み合わ...
上記のプレテンション・死荷重・活荷重を全て由利橋モデルに...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/al...
max変位 ＝ +917.82mm
プレテンション・死荷重・活荷重を足し引きしたものとほぼ同値.
デフォルトで約1mもたわむものなのかと少し不安…
**主塔３Ｄモデル作成_midas(6/20～6/26) [#xa7a1e15]
マイダスで主塔の３Ｄモデルを作成した.
やり方は ウィザード＞基本構造＞シェル＞数値入力 からで...
手順は写真の通りである.
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/od...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/od...
ここでいう、
ｎ：ｙ軸方向の分割数
ｍ：ｘ軸方向の分割数
Ｉ：ｚ軸方向の分割数
である.
今回は、ｎ＝10コ,ｍ＝10コ,Ｉ＝50コ にした.なお適当.
設定した後は👇このようになった
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/sy...
全要素に400℃を与えて解析したところDX,DY,DZ方向全てに変位...
なお、まだ塑性や降伏点は設定していない.
〇DX方向変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/dx...
max変位＝7.62mm
〇DY方向変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/dy...
max変位＝4.20mm
〇DZ方向変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/dz...
max変位＝240.64mm
〇DXYZ方向変位
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/dx...
max変位＝240.77mm
一つ問題点…境界条件を節点を完全固定にしているため、３Ｄに...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ky...
失敗したやり方…節点を作成して要素で結んで、後から断面や材...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/so...
塑性材料を設定しようとすると「静的材料非線形関連機能をサ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/so...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ch...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ca...
**1/100主塔モデルで温度解析_Salome(7/3〜7/10) [#ke8a0fe8]
1/100モデルにした理由は、Asterstudyの設定以外でのエラーを...
7/3：合っているかはわからないがある程度Asterstudyの設定が...
残り、✅熱膨張係数（ELAS-ALPHA）・✅初期温度（AFFE_VARC-T...
を行えば良い.やり方わからなかった.聞かなければ、、
温度解析については[[Salomeメモ]]の熱関係と[[Salome-Meca_...
7/4：なんとか設定一通り終えて解析をした.成功もした.
結果を見た感じ変なところはなさそうだが、マイダスとの結果...
これは、マイダスはシェル（厚さ10mmの中は空洞）で、サロメ...
確認が必要かもしれない.
↑この理由は1/100モデルだからか、、忘れてた。。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/se...
結果
今回はマイダスと比較するために、先端部分の大体マイダスと...
上記写真の上端である.
なお単位はmmである.かつサロメは結果の数値を意図的に100倍...
,,マイダス,サロメ
,DX,7.62,9.4
,DY,4.20,5.8
,DZ,240.64,241
7/5:熱解析のやり方をまとめ、上記のグラフを出してみた.また...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/se...
結果的には、100℃のときの最大変位はDX=0.023mm,DY=0.014mm,D...
この値は400℃のときの最大変位の約1/4であることがわかった.
400℃のとき、DX=0.094mm,DY=0.058mm,DZ=2.41mm
このことから熱解析自体は正しくできているだろう.
次はいよいよ弾塑性を入れてみるか、、
7/9:弾塑性のやり方を学んだのと、実際の大きさで車高部分に4...
モデルは、下から0.5m,0.5m〜3m,3m〜50mの３分割にした.
車高は1.5mとし、そこから下に1mかつ上に1.5mの部分が真ん中...
このモデルの真ん中のBOXにのみ400℃与えた.
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ta...
結果はこのようになり、加熱部が一番大きく変形した.
Pravisの写真の点部分.
３つのBOXの手前左上.
,,DX,DY,DZ
,点0,7.48e-29,4.18e-29,9.19e-29
,点8,3.69,1.69,12.9
,点12,-0.014,0.0006,9.79
弾塑性入れてやってみたが破壊形態は見れなかった？見方がよ...
今は瞬間的に0℃→400℃に変えて解析をしているが、現実性を考え...
！補足！
マイダスと比較するために、サロメでも実モデル全体に熱を加...
,,マイダス(mm),サロメ(mm),相対誤差(％)
,DX,7.62,7.20,5.5
,DY,4.20,3.61,16.3
,DZ,240.64,241.1,0.19
熱解析自体は上手くいっていると思われる.
7/10:ゼミ
ゼミでただの鋼材に熱を加えるだけでは膨張するだけだと言わ...
そこで粘性係数を考慮する必要がある.
熱解析で鋼材が曲がるか確認するためにまずは簡単な片持ち梁...
**改めて、リスクアセスメントとは(7/17) [#e5d51923]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ri...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ri...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ri...
**粘性を入れた片持ち梁の解析(7/19〜7/24) [#td0548e7]
まず粘性の入れ方を学ぶために10×10010
まず粘性の入れ方を学ぶために10×100×10の片持ち梁で行った
[[CodeAster:https://code-aster.org/doc/default/en/man_u/u...
p137〜140に書いてある
[[ヒルの降伏応力:https://www.jsme.or.jp/jsme-medwiki/doku...
[[円筒座標系:https://www.cradle.co.jp/glossary/ja_A/detai...
[[直交座標系:https://www.cradle.co.jp/glossary/ja_T/detai...
[[テンソルとは:https://3.1415.jp/35ss9vvx]]
[[hillファイル:http://www.str.ce.akita-u.ac.jp/~gotouhan/...
CodeAsterを解読した
SalomeのLNTは下のように方向を表す
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/LN...
温度における、ac1＝準静的温度開始、ac3＝準静定温度終了
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ac...
これはAstersutudyの
DEFI_MTERIAU>META_LEMA_ANI(粘性)の設定画面である
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/an...
CodeAsterより等方性の場合はこのようになる👇
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/to...
γ、θ、Z方向は円筒座標系であり、x、y、z方向は直交（デカル...
この２つの座標系の関係性は
[x=γcosθ]
[y=γsinθ]
[z=z]
となる
①熱
②熱＋弾塑性
③熱＋粘性
④熱＋弾塑性＋粘性
①のPravis http://www.str.ce.akita-u.ac.jp/~gotouhan/j2...
②のPravis http://www.str.ce.akita-u.ac.jp/~gotouhan/j2...
①と③の結果は同じになり、②と④の結果は同じになった、つまり...
鋼材の動粘性を"0.5mm^2/s"として入れた
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ne...
[方向]
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) [#j842ba5a]
・ソリッド要素
・加熱する時間の設定１０〜３０分
・ケーブルの応力分布
・時間ごとに熱くなるのを再現できるか
7/25:ソリッド要素断念
主塔とC6・C7ケーブルをソリッド要素で作成を試みた
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ca...
しかし、下の写真のように主塔とケーブルのつなぎ目部分が上...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ca...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ca...
移動や回転も試してみたが、思った方向に変化しなかった
ビーム要素に変更する
7/26:ビームモデルの作成
[[code-asterによる解析 step-by-step:https://www.opencae.o...
[[構造要素の概要:https://www.opencae.or.jp/wp-content/upl...
ビーム要素でモデルを作成した
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/be...
ケーブルと主塔の間に小さいバネが必要であるため、以下のよ...
今回のモデルでは10mmの長さにした
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/be...
主塔やケーブルには架空の断面を設定する必要がある
主塔：Geometryの長方形フェース ケーブル：Geometryの円...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/sy...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ca...
写真ではケーブル断面を実際のケーブル位置に作成したが、そ...
原点に作っておいてAsterstudyで割当できるらしい
ビーム要素のも出るを作成するに当たっての注意がある
①架空の断面はGeometryでグループ作成してから、Meshでジオメ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ge...
②Meshのコンパウンド必須、これはビーム要素のみで断面はしな...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/co...
③②でしたCompound_Meshからグループを作成で「ノード」ー固定...
○モデルのSalome設定一覧
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/be...
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 le...
! l'arrêt. ...
! ...
! Conseils : ...
! - Vérifiez vos conditions aux limites. ...
! - Vérifiez votre modèle, la cohérence des unit...
! - Si vous faites du contact, il ne faut pas qu...
! que par le contact. ...
! ...
! - Parfois, en parallèle, le critère de détecti...
! est trop pessimiste ! Il reste néanmoins souvent ...
! possible de faire passer le calcul complet e...
! (augmenter de 1 ou 2 la valeur du mot-clé NPREC) ou ...
! en le débranchant (valeur du mot-clé NPREC=-...
! sur moins de processeurs.
! <S> ユーザー例外が発生しましたが、インターセプトされ...
! ベースが閉じられた！
! 例外タイプ: matrixSingularError ！
! !
! 非可逆行列のため停止．
! グローバルデータベースが保存されました。グローバ...
シャットダウン
! !
! ヒント: ！
! - 境界条件をチェックする！
! - モデルと単位の整合性をチェックする！
! - コンタクトを使用する場合、構造は "保持 "して...
接触だけである！
! !
! - 並行して、MUMPSの特異点検出基準！
特異点検出基準は悲観的すぎる！とはいえ、この基準を緩和...
とはいえ、この基準を緩和することで、完全な計算をパスす...
! (NPRECキーワードの値を1または2増やす)、あるいは！
切断する（NPRECキーワードの値を1にする）か、計算を再実...
より少ないプロセッサで計算を再実行する。
Degré de liberté physique associé au noeud N27 et à la...
Problème : la matrice est singulière ou presque singul...
Lors de la factorisation de la matrice, on a rencontré...
(pivot nul ou presque nul) à la ligne 174 qui correspo...
Risques et conseils :
* Si la ligne correspond a un degré de liberté physiqu...
de corps rigide mal bloqué.
Vérifiez les conditions aux limites.
Si vous faites du contact, il ne faut pas que la str...
Vérifiez également les caractéristiques matériaux (m...
Si la ligne correspond a un degré de liberté de Lagran...
limite redondante.
En particulier, il se peut que la relation linéaire ...
Peut-être devriez vous exclure certains noeuds des c...
(mots clés SANS_NOEUD et SANS_GROUP_NO).
* Si le solveur utilisé est LDLT ou MULT_FRONT, vous p...
car celui-ci est le seul à pouvoir factoriser les ma...
* Parfois, en parallèle, le critère de détection de si...
possible de faire passer le calcul complet en relaxa...
en le débranchant (valeur du mot-clé NPREC=-1) ou en...
* Il se peut aussi que ce phénomène soit tout à fait n...
très près d'un noeud.
Si le nombre de décimales perdues n'est pas trop gra...
vous pouvez relancer le calcul en augmentant le nomb...
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ま...
(NPRECキーワードの値=1)、あるいは、より少ないプロセ...
* クラックがノードの非常に近くを通過する場合、この現象...
節点に非常に近い場合。
失われる小数点以下の桁数があまり大きくない場合（最大...
失われる小数点以下の桁数を増やして計算を再開すること...
SOLVEUR factor キーワードの NPREC キーワード。
そうでなければ、開発チームに連絡してください。
S> ユーザー例外がスローされたが、インターセプトされな...
! ベースはクローズされている！
! 例外タイプ: MatrixSingularError ！
! [('?', (), (), ())]
N27はバネとケーブルの節点
バネの条件及び剛性がおかしいとマトリクスエラーが出る
短すぎる？長すぎる？
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/n2...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/mm...
8/5:ケーブル１本と２本で比較とエラーの改善
[[このときのAsterstudy:http://www.str.ce.akita-u.ac.jp/~g...
ノードN27とDRZコンポーネントに関連する物理的自由度
がおかしいとエラーが出るため,試しにDRZを固定して見ると解...
果たして「ピン」としての役割を果たせているのだろうか？
Degré de liberté physique associé au noeud N2 et à la c...
Problème : la matrice est singulière ou presque singuli...
Lors de la factorisation de la matrice, on a rencontré ...
(pivot nul ou presque nul) à la ligne 12 qui correspond...
Risques et conseils :
* Si la ligne correspond a un degré de liberté physiqu...
de corps rigide mal bloqué.
Vérifiez les conditions aux limites.
Si vous faites du contact, il ne faut pas que la str...
Vérifiez également les caractéristiques matériaux (m...
* Si la ligne correspond a un degré de liberté de Lagr...
limite redondante.
En particulier, il se peut que la relation linéaire ...
Peut-être devriez vous exclure certains noeuds des c...
(mots clés SANS_NOEUD et SANS_GROUP_NO).
* Si le solveur utilisé est LDLT ou MULT_FRONT, vous p...
car celui-ci est le seul à pouvoir factoriser les ma...
* Parfois, en parallèle, le critère de détection de si...
possible de faire passer le calcul complet en relaxa...
en le débranchant (valeur du mot-clé NPREC=-1) ou en...
* Il se peut aussi que ce phénomène soit tout à fait n...
très près d'un noeud.
Si le nombre de décimales perdues n'est pas trop gra...
vous pouvez relancer le calcul en augmentant le nomb...
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:熱が入らない・・
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/pr...
まず、プレテンションのみとそこに弾塑性を加えた場合だと、...
しかし、そこに熱を入れると解析は回ったが、結果も変わらな...
そのときのAsterstudyは以下の通りである
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/be...
8/14:エラーの原因と改善が不明
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/be...
これに
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 / ...
! Pour la variable de commande TEMP ...
! la grandeur associée du champ doit être: TEMP_R ...
<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, ...
! n'est pas compatible avec le mot clé utilisé (VAL...
! ...
! 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) [#i...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/no...
【Case1】
主塔0m~0.5mの材料 SM400
主塔0.5m~3mの材料 E=100000N/mm^2,ν=0.3,α=1.2e-05/,w=7.7e...
主塔3m~50mの材料 SM400
ケーブルの材料 E=195000N/mm^2,ν=0.3,α=1.26e-05/,w=7.85...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/E1...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/E1...
上記の結果を見たところ、どちらのパターンでも主塔高さ3mの...
【Case2】
主塔0.5m~3mの材料 E=150000N/mm^2,ν=0.3,α=1.2e-05,w=7.7e-...
ポッキリ折れた
!改善点として主塔の中に２枚のプレートを入れる
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/E-...
片側側面に使用したヤング率は上記のグラフを参考にした
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ka...
色のついた高さ0.5m~3mの部分が火災想定（主塔片側のみ）
【Case3】～主塔の車道側側面のみを加熱した想定
ー400℃のとき
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/40...
プレテンションの強い右の方に引っ張られる形になった
このときケーブル応力は、C1:-45.9N/mm^2,C12:28.2N/mm^2とな...
9/24:青木ゼミ
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/c_...
まだ、途中…
引張応力と降伏応力の関係がよくわかっていない
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ma...
交通事故・交通量・火災発生件数等をできる限り調べてみた
**シェル要素片持ち梁と平面要素片持ち梁の比較とE-℃グラフ(9...
片持ち橋の寸法 幅4m×厚さ2m×長さ10m 温度500℃
,,シェル要素(プレートなし),平面要素(プレートなし),平面要...
,DZ方向変位(mm),30.985,30.985,21.380,21.380
,応力(N/mm^2),312.0,312.6,265.8,265.9
,ヤング率(GPa),101,101,124,124
σ＝E/ε から求めた
9/12:ヤング率と温度の関係を模索
上記の片持ち梁を使って熱と温度の関係を求めようとしたが、...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ka...
そのため、青木さんからいただいた論文をもとにグラフを作成...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/E-...
**由利橋への影響を考える～県道交通事故・普通貨物車の通行...
[[日本の県道一覧:https://altcva.main.jp/pref/]]
[[道路の種類:https://www.mlit.go.jp/road/sisaku/dorogyous...
[[秋田県内橋梁一覧:https://river.longseller.org/bp/5.html]]
[[令和３年度の一般交通量調査の結果について:https://www.pr...
このサイトのExcelファイルが秋田県内の交通量を調査した結果...
国道及び県道の交通量かつ大型車混入率が示されている
[[道路交通に関する統計:https://www.e-stat.go.jp/stat-sear...
月別事故発生数と死者数が見比べられる
また、都道府県別交通事故発生状況を見ることもできる
[[秋田警察ホームページ:https://www.police.pref.akita.lg.j...
ホームページの「ミニ統計」からその年の秋田県県道の事故件...
（令和4年用ミニ統計）
[[R4HP用ミニ統計:https://www.police.pref.akita.lg.jp/uplo...
（令和5年の秋田県道事故発生件数）
[[交通統計:https://www.police.pref.akita.lg.jp/uploads/co...
このpdfの7ページに道路別事故件数が記されている
[[東北管轄交通事故発生状況:https://www.tohoku.npa.go.jp/t...
このpdfの３に交通死亡事故の特徴が示されているが、ここから...
[[H17東北の事故発生状況:https://www.thr.mlit.go.jp/road/i...
[[H13～H17の東北管轄県道事故件数棒グラフ:https://www.thr....
データは古いが東北管轄の県道の事故件数13073件だった
[[令和5年山形県道交通事故:https://www.pref.yamagata.jp/do...
山形県の場合は858件
[[総務省消防庁火災統計:https://www.fdma.go.jp/pressreleas...
年別の火災に関する統計諸々が記されている
下記が最新（令和5年）の火災統計
[[令和５年(１月～12 月)における火災の概要:https://www.fdm...
下記は令和4年の短め資料
[[報道資料:https://www.fdma.go.jp/pressrelease/houdou/ite...
[[火災曲線:https://www.oit.ac.jp/japanese/sangaku/OITID-2...
[[火災時におよび火災後におけるスタッドの力学特性:https://...
図４に加熱温度ー時間関係のグラフあり
**0,400,800℃の３モデルで解析(9/26～10/3) [#vc55d915]
・健全(0),400,800の３パターンを，D+集中荷重/D+集中荷重+L...
・ケーブル・桁・主塔の応力をそれぞれ見る
ケーブル1160N/mm^2で伸び1570N/mm^2で降伏，桁235N/mm^2で降...
・プロパンの着火温度は約460℃～520℃に対して，ガソリンの着...
[[プロパンガス・ガソリンの発火温度:https://www.elg-inc.jp...
[[ガソリンの発火温度:https://www.yellowhat.jp/column/oil/...
・伊勢湾岸自動車道の事故の写真や状況を参考に中間発表
[[伊勢湾岸自動車道爆発事故:https://www.khk.or.jp/public_i...
・都市ガス普及率から秋田はプロパンガスの割合が高い
[[都市ガス普及率:https://uub.jp/pdr/eg/gas.html]]
モデルの全体像は以下の通りである
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ze...
モデルの主塔と桁のつなぎ目はくっつけることに修正した
[[midas平面応力度について:https://manual.midasuser.com/JP...
[[midasリブについて:https://support.midasuser.com/hc/ja/a...
9/30:おかしい点
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/st...
このように桁と主塔の節点のひとつ上の節点に大きな力がかか...
つなぎ方に問題がなる？
**全体モデルの作成(10/9〜10/25) [#pa80f758]
全体モデルを作成した
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/mo...
以下のように，主塔と桁はシェル要素（板要素），ケーブルは...
桁と主塔は平面応力要素にするとケーブルとうまくつながらな...
これは，平面応力要素にZ方向の応力が発生しないことが原因？...
これからは平面応力要素は使用せずに板要素を用いる
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/mo...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/si...
モデルの境界条件：桁 ローラー(x=-135) ピン(x=55)、 主...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/si...
死荷重は全て圧力荷重で載荷
その際、面載荷ができなかったため、死荷重を/4して節点載荷...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/da...
主塔・桁断面はこのようにした
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ka...
主塔周りの４枚板要素は鋼材にすると主塔が降伏したため、可...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/me...
桁のメッシュは基本鋼材もコンクリートも1×1として、設計上必...
○部材
主桁 19m×0.3m×190m SM400
主塔 3m×1.5m×70m（高さ50m） SM400
ケーブル D=200mm(C5〜C9),225mm(C1,C4,C10),315mm(C2,C3,C1...
**健全時・400℃・800℃のときの解析結果(11/6～11/7) [#s93d49...
○健全時と400℃の比較
あまり大差ない結果となった
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/23...
上記のC1ケーブルと桁との節点が最も応力のかかった結果とな...
,,健全,400,800
,集中荷重,220MPa,224MPa,224MPa
,集中荷重＋分布荷重,225MPa,230MPa,230MPa
,変位,-914mm,-934mm,-934mm
ヤング率を下げているにも関わらず変位と応力が400℃，800℃の...
また，ヤング率を下げた主塔と下げていない主塔とで変形が重...
上手くヤング率が反映されていない？
また主塔の0.5m～3.0mを加熱しているためか，下の写真のよう...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/22...
**IABSEに向けてリスクの定量化(11/12～11/26) [#dcbc9f10]
[[秋田県の平均通勤時間:https://statresearch.jp/house/rent...
ここから計算すると全国平均通勤時間：25.3分、秋田平均通勤...
[[令和4年道路形状別交通事故発生件数:https://www8.cao.go.j...
[[令和5年道路形状別交通事故発生件数:https://www8.cao.go.j...
この2つのサイトは交差点・単路・カーブといった道路形状別の...
由利橋付近は一般単路と交差点付近の値が参考になるだろう
[[交通事故が起こりやすい道路形状:https://agoora.co.jp/jik...
これも同様内容、参考までに
3年分のデータから一般単路は大体31%くらい
[[秋田県交差点での事故率:https://www.sonpo.or.jp/about/us...
秋田県のデータ、令和5年のもの
[[H8~H17間の全国の橋梁上の事故件数:https://www.mlit.go.jp...
件数：58,878 延長：8,867km 延長1km当たりの事故発生...
[[事故類型別人身事故発生件数割合:https://www.ms-ins.com/s...
ほとんどが相互車両の事故
[[高速道路の事故率:https://www.mlit.go.jp/road/ir/ir-coun...
ここから分かる通り、自動車専用道路よりも一般道路の方が事...
歩行者等の人がいると予期せぬ行動が起こりやすいため高くな...
[[事業用自動車の走行距離と交通事故件数、交通事故死者数:ht...
[[令和3年全国市区町村別交通事故死者数:https://www.itarda....
↑これも参考までに
[[令和５年中の道路交通事故の状況:https://www8.cao.go.jp/k...
高齢者の死亡事故や事故に逢った状況が詳しく書いてある
[[消防白書:https://www.fdma.go.jp/publication/hakusho/r5/...
p28に記載
[[車両別事故・火災情報:https://renrakuda.mlit.go.jp/renra...
! 車種別事故・火災情報件数及び割合は、貨物車が 519 件...
353 件（27.0%）となっている。事象別にみると、事故情報件数...
（43.4%）と最も多く、次いで乗用車が 57 件（24.3%）となっ...
物車が 417 件（38.9%）と最も多く、次いで乗用車が 296 件（...
月末現在における乗用車の保有台数（38,938 千台）は全保有台...
り、乗用車の事故・火災情報件数は乗用車の普及状況が反映さ...
らに、貨物車は保有台数 6,175 千台（7.4%）に対して事故・火...
逆に軽乗用車は保有台数 23,388 千台(28.2%)に対して事故・火...
くなっている。
p5抜粋
この試資料はなんの事故なのか・なんの火災なのかがわからない
[[日本の総道路延長:https://www.mlit.go.jp/road/soudan/sou...
令和3年のデータ
都道府県道の総道路延長：142,942.2 km
秋田県道：2376.7km
[[2023道路統計年報:https://www.mlit.go.jp/road/ir/ir-data...
道路統計年報では、令和4年時点での道路及び橋梁現状を閲覧で...
[[交通事故件数減少推移グラフ:https://www.jama.or.jp/opera...
交通事故発生件数の減少推移を閲覧できる
・交通事故で死ぬ確率ー1年間に10万人あたり4人程度
・一生のうちに交通事故に遭う確率ー25％ 四人に一人...
・一生のうちに交通事故で死亡する確率ー1％ 年間約1万人...
・年間10万人あたり４人が交通事故で亡くなっている
・約0.2％（500人に1人）が交通事故に遭っている
・約0.4％（250人に1人）が1年間に交通事故を起こしている ...
[[交通事故に逢う確率＆起こす確率:https://chokottokyosai.j...
以上の結果から定量化に努める
○県道での事故発生率
○全国橋梁上での事故率
○大型車混入率
○一般単路の事故率・交差点付近の事故率
○全事故のうち、死亡者や重傷者が出る割合
○通勤時間の増加
などなどを考慮する
井上さんからいただいた最終的な5.75×10＾-7/km/年はどう計算...
Anser:採用された仮定を全てかけ合わせる
この計算方法がわかれば、あとは橋長190m(0.19km)をかけて/年...
その計算した確率が1×10^-6/年よりも小さければリスクを考慮...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ri...
リスクの定量化を行った
最終的な由利橋のリスクは 2.02×10^-7/年 になった
「リスクを高い低いと判断することよりも、場面に合わせたリ...
と結論付けるの方がまとめやすいかもしれない・・・
一応基準は人が1年間で死ぬ確率 4×10^-5/年 にした
今回は思っていた通りこのリスクよりも小さい
※○印のついたものが定量化に使用した値
[計算式↓]
7.82E-11/km^2・人・年 × 橋長：0.19km^2 × 由利本荘市...
11/22:青木さんと発表練習後
由利橋のリスクはそもそも死ぬ確率ではないので、交通事故死...
トルコの事例と比較の方がいいかもしれん
0.3/km→0.3,0.15/km→0.15,0.09/person→0.09/day×365ay=32.85
それぞれ修正
全てを確立に変換して、○○○○/km・年を出す
そして、最終的に橋長をかけて○○○○/年を求める
これが1年間に起こり得る由利橋でのリスクとなる
まとめでは、解析結果から火災が起こるリスクは低いかもしれ...
11/26:リスクアセスメントの改訂版
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ri...
**要素温度荷重解析の結果(11/13〜11/14) [#k9303318]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/E-...
これからは、ヤング率を下げて解析するのではなく”要素温度荷...
ただし、この方法は徐々に熱が加わっていくことを想定してい...
これは今までのヤング率を下げたものを使っていたときと同様...
健全モデルのメッシュをもっと細かくする必要あり
解析結果は以下の通りである
,,健全,250,270,300,400,800,1200
,主塔応力（集中荷重）,78.3MPa,316.9MPa,339.4MPa,373MPa,48...
,主塔応力（集中荷重＋分布荷重）,80.4MPa,317.6MPa,340MPa,3...
,C1応力（集中荷重）,220MPa,220MPa,220MPa,220MPa,220MPa,22...
,C1応力（集中荷重＋分布荷重）,226MPa,226MPa,226MPa,226MPa...
,C12応力（集中荷重）,93.1MPa,93.3Mpa,93.3MPa,93.3MPa,93.4...
,C12応力（集中荷重＋分布荷重）,95.3MPa,95.45MPa,95.5MPa,9...
,x軸変位（集中荷重）,-47.4mm,-47.7mm,-47.71mm,-47.71mm,-4...
,x軸変位（集中荷重＋分布荷重）,-49.14mm,-49.15mm,-49.15mm...
,y軸変位（集中荷重）,3.16mm,203.72mm,219.77mm,243.83mm,32...
,y軸変位（集中荷重＋分布荷重）,3.05mm,203.61mm,219.66mm,2...
,z軸変位（集中荷重）,-891.16mm,-891.55mm,-891.59mm,-891.6...
,z軸変位（集中荷重＋分布荷重）,-914.35mm,-914.73mm,-914.7...
以下３枚が270℃で降伏した時の変形図である
これを見たらわかるように、主塔が桁の外側へ倒れ込むような...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/27...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/27...
以下２枚が主塔と桁（板要素）の応力分布である
２枚目からわかる通り、主塔の実際に加熱した部分が270℃以降...
写真左の方が応力が集中しているのは、メインスパンの方が長...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/27...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/27...
グラフ１枚目[heat-displacement]は熱とy軸方向変位の関係グ...
ただ、ケーブル応力はどの温度のときでもほとんど変わらなか...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ya...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/he...
グラフ２枚目[heat-stress]は熱と主塔応力の関係グラフであり...
270℃は大体10分程度で到達されると予想
消防車が来るまでの平均時間は約８分、消火活動の準備などで...
一般の鋼材は400℃で降伏するとされているが、だいぶ低い段階...
これは主塔の厚さや死荷重、活荷重、プレテンションなどが相...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/27...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/he...
なお、解析ではケーブルの応力はほとんど変化しない結果とな...
○補足
初期温度20度とすると400℃へ行くまでに約20分ほどかかること...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/he...
このことから、実際に主塔が降伏するまで火災が続く可能性は...
熱と時間に関しては要検討！
↑これ間違いかも…質量計算し直すと270℃まで200分かかると出て...
**ケーブル火災解析(12/18) [#nd8a7ba3]
ケーブルの要素分割をすると、ぽきぽきになってしまい解析で...
つなぎ目をどうするか考えて、プレテンションの入れ方を模索...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/a1...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/b1...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/c1...
ケーブル火災解析は、卒論後に検討でもいいかも
**リスクの修正(12/20〜2/6) [#xe27b376]
今までのまとめとリスクの修正を記した．
ユーザーリスクと通行止めリスクの２つを考え，それぞれを検...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ri...
[Eurocodeについて]
[[Eurocode:https://eurocodes.jrc.ec.europa.eu/]]
[[欧州構造基準（ユーロコード）の制定─その体系と内容が及ぼ...
[[ユーロコードの完全実施と我が国への影響:https://www.pwrc...
[[ユーロコードの構成を俯瞰:https://note.com/issey96/n/n5e...
[[内閣府 令和3年中の道路交通事故の状況:https://www8.cao....
主塔で起こる火災を
事故が原因の場合・・正面衝突事故率×橋梁に対する主塔の割合
設備不備の場合・・車両火災中設備からの出火事故率×橋梁に対...
とか？？
中間発表後にリスクの修正を行った
〇参考資料
[[大型車の車両火災発生状況について:https://www.ntsel.go.j...
[[The influence of wind and the spatial layout of dwellin...
[[Effect of wind speed and direction on facade fire sprea...
[[Numerical study of effects of wind on the vertical fire...
[[Experimental study of the temperature characteristics o...
[[Influence of Wind Direction on Fire Spread on High-Rise...
これに「図4. 鉛直風速の変化則」に風速と高さとの関係グラフ...
***限界状態について [#va8f246b]
[[道路橋示方書:https://www.mlit.go.jp/road/sign/kijyun/pd...
限界状態がどれであれ通行止めは発生する。
そのため、特にリスクアセスメントでは限界状態を考慮してい...
1) 橋の限界状態 1
橋としての荷重を支持する能力が損なわれていない限界の状態
2) 橋の限界状態 2
部分的に荷重を支持する能力の低下が生じているが，橋として...
能力に及ぼす影響は限定的であり，荷重を支持する能力があら...
にある限界の状態
3) 橋の限界状態 3
これを超えると構造安全性が失われる限界の状態
1) 死荷重 (D)
2) 活荷重 (L)
3) 衝撃の影響 (I)
4) プレストレス力 (PS)
5) コンクリートのクリープの影響 (CR)
6) コンクリートの乾燥収縮の影響 (SH)
7) 土圧 (E)
8) 水圧 (HP)
9) 浮力又は揚圧力 (U)
10) 温度変化の影響 (TH)
11) 温度差の影響 (TF)
12) 雪荷重 (SW)
13) 地盤変動の影響 (GD)
14) 支点移動の影響 (SD)
15) 遠心荷重 (CF)
16) 制動荷重 (BK)
17) 橋桁に作用する風荷重(WS)
18) 活荷重に対する風荷重(WL)
19) 波圧 (WP)
20) 地震の影響 (EQ)
21) 衝突荷重 (CO)
22) その他
1) 永続作用による影響が支配的な状況（永続作用支配状況）
①D +PS+CR+SH+E+HP+(U) +(TF) +GD+SD +WP +(ER)
2) 変動作用による影響が支配的な状況（変動作用支配状況）
②D+L+I+PS+CR+SH+E+HP+(U) +(TF)+(SW)+GD+SD+(CF)+(BK) +WP +...
③D +PS+CR+SH+E+HP+(U)+TH+(TF) +GD+SD +WP +(ER)
④D +PS+CR+SH+E+HP+(U)+TH+(TF) +GD+SD +WS +WP +(ER)
⑤D+L+I+PS+CR+SH+E+HP+(U)+TH+(TF)+(SW)+GD+SD+(CF)+(BK) +WP...
⑥D+L+I+PS+CR+SH+E+HP+(U) +(TF) +GD+SD+(CF)+(BK)+WS+WL+WP ...
⑦D+L+I+PS+CR+SH+E+HP+(U)+TH+(TF) +GD+SD+(CF)+(BK)+WS+WL+W...
⑧D +PS+CR+SH+E+HP+(U) +(TF) +GD+SD +WS +WP +(ER)
⑨D +PS+CR+SH+E+HP+(U)+TH+(TF)+(SW)+GD+SD +WP+EQ +(ER)
⑩D +PS+CR+SH+E+HP+(U) +(TF) +GD+SD +WP+EQ +(ER)
3) 偶発作用による影響が支配的な状況（偶発作用支配状況）
⑪D +PS+CR+SH+E+HP+(U) +GD+SD +EQ
⑫D +PS+CR+SH+E+HP+(U) +GD+SD +CO
(3) (2)1)から 3)に規定する作用の組合せに対して，表-3.3.1 ...
重係数を考慮する。
γ୮ ：荷重組合せ係数であり，異なる作用の同時載荷状況に応じ...
考慮する作用の規模の補正を行うための係数。
γ୯ ：荷重係数であり，作用の特性値に対するばらつきに応じて...
慮する作用の規模の補正を行うための係数。
なお，活荷重に対する衝撃の影響（I）を考慮するにあたって，...
は荷重組合せ係数γ୮及び荷重係数γ୯を乗じる必要はない。
(4) 風荷重については必要に応じて他の作用を考慮しない場合...
を適切に設定する。
(5) 衝突荷重及び制動荷重については死荷重及び活荷重のみと...
(2)1)2)以外の条件を適切に設定する。
***由利橋について [#u307a94f]
由利橋はH23.2.22に施工が始まり、H25.1.31に完成した。
橋長190.5mの不等径間の２径間連続斜張橋（鋼床版２主箱桁）...
秋田県由利本荘市の子安川にまたがる形で位置している。
由利本荘市：面積1209.59km^2,人工:74707人
小吉川：一級河川,全長61km,流域面積1190km^2,小洪水は毎年起...
今までで由利橋で一度も事故が起こったことはない！！→なぜ？...
予想可能な災害…土砂崩れ（出羽丘陵）、氾濫及び洪水（小吉川...
[[由利本荘市ハザードマップ:https://www.city.yurihonjo.lg....
考え得る危険事項…交通事故（追突・接触）、☆火災（自動車事...
使用不可になった場合の懸念点…移動手段の減少による渋滞・移...
※このときの迂回路…飛鳥大橋、本荘大橋、…子吉川が流れている...
[[令和4年：交通事故統計:https://www.police.pref.akita.lg....
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/a....
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/b....
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/c....
[由利橋で事故が起こらない理由]
・単純にスパン長が短い
・車両の速度が比較的遅い
・ほとんどが軽自動車か普通自動車、たまに軽トラックであっ...
・渡り終わったすぐ先が信号機であるから慎重な運転となるため
・見通しが良い
・歩行路と道路がフェンスで区切られているため、歩行者との...
***落橋について [#g86999e7]
落橋とは、橋梁の劣化や大規模な災害などが要因となって橋梁...
落橋の主な原因は橋梁の老朽化であり、損傷の種類は[[日本橋...
また、大規模な災害、具体的には地震や津波、土砂災害、強風...
落橋がもたらす社会的影響は、周囲の交通機関に不便をもたら...
***ケーブル腐食について [#o082b4da]
腐食とは、化学や生物的な作用によって物体の外観や機能が損...
金属の腐食を特に「錆」という。
腐食の主な要因は、塗装劣化個所に水分と腐食性物質が作用す...
腐食性物質としては、大気中の亜硫酸ガスと海塩粒子等が挙げ...
これらは地域や気候に関係が深く、工業地域や海岸近くに位置...
発生のメカニズムとしては金属は大気中の酸素によって表面に...
この酸素と水分、金属表面の汚染物などが表面に付着し酸化被...
ここから腐食が表面から内部へと進行していく。
金属にとって「錆」状態は安定している。
腐食の種類は均一腐食と局部腐食の２つである。
均一腐食…金属表面全面に発生する腐食。
局部腐食…局部的な腐食で発見しにくい。一般的な製造部品など...
ここから本題である「ケーブル腐食」を考える。
斜張橋と吊り橋は両橋ともケーブルを用いる形状をしているが...
斜張橋は、主塔から張られているケーブルが直接桁橋を支える...
また吊り橋はアンカーレッジが存在する。
斜張橋の強みとして、桁曲げモーメントを小さくできること・...
スパン長200～500mm及びそれを超えると斜張橋か吊り橋で設計...
それぞれの橋梁の桁に加わる力として、
斜張橋…垂直方向の張力＋橋軸方向の圧縮力
吊り橋…垂直方向の張力
がある。
ケーブルは通常「亜鉛メッキ銅線」で構成されているらしい。
亜鉛メッキとは高い耐食性（腐食のしにくさ）を付与出来るメ...
亜鉛金属を電気の力で析出させる方法と溶融させた亜鉛に浸漬...
亜鉛メッキは防錆効果が高く主に鉄製品に対して処理を行いま...
また、亜鉛は鉄よりイオン化傾向が大きいのでメッキ皮膜にピ...
これを犠牲防食という。
ただし、デメリットとして亜鉛メッキは湿気に弱い。
ケーブルが腐食する主要因は「水（＋塩）」であると考える。
特に、乾湿を繰り返す部位は、常に水に覆われた部位より腐食...
その理由は、乾湿を繰り返す部位が常に水に覆われた部位と比...
ケーブル腐食がもたらす影響をいくつか考える。
①斜張橋の終局強度の低下…ケーブル断面の減少＋伸び及び疲労...
②疲労寿命（破壊が起こるまで材料が耐えられる、ある種の変動...
③橋梁全体の耐力低下…ケーブル断面の減少 から
④ケーブルの伸び及びねじれ強度の著しい低下…亜鉛メッキが腐...
上記①～④の原因・理由を考えるために[[斜張橋ケーブルの耐久...
-ケーブル断面の欠損により、ケーブル鋼線が降伏するまでは斜...
-小段数ケーブル斜張橋の場合は、ケーブル破断の寸前になると...
-マルチケーブル斜張橋の場合は、ケーブル１本が破断しても特...
-応力腐食割れ（腐食環境下において金属材料に引張応力が作用...
ケーブルは主塔・主桁との協働作用により外力に抵抗する非常...
これにより腐食したケーブルで耐えている斜張橋や吊り橋は、...
橋梁の突発的な崩壊は絶対に避けるべき事故である。特に斜張...
従って、当たり前ではあるがケーブル破断が起こらないように...
--補足
風によるケーブルの振動現象は
①カルマン渦による渦励振
②六角形断面ケーブルにおけるギャロッピング
③並列ケーブルにおけるウェイクギャロッピング
④レインバイブレーション
⑤塔からの剥離渦によるウェイクレゾナンス
***リスクアセスメントについて [#zdb76e24]
[文献から…]
リスクアセスメントとは、危険が発生する可能性の度合いを組...
橋梁におけるリスクアセスメントは、橋梁がどのくらいの割合...
欧米ではリスクアセスメントが主流なため「何が現実に実行で...
日本ではこの考えはまだまだ定着しておらず、リスクゼロを前...
しかし実際にはリスクゼロは不可能であるため、日本でもリス...
[リスクマネジメント&リスクアナリシスの流れ]
基本条件・情報の整理→リスクの定量化・ランク付け→対策を実...
[ALARP：As Low as Reasonably Practicableの原則]
① 「受け入れ可能なリスク」として「自然災害の発生レベル以...
② リスク低減のためにどこまでも費用をかけることが合理的で...
③ 低減措置の実施基準は，リスク低減効果と低減に要する費用...
また安全に対する考え方で、，「広く受容される領域」・「我...
ALARP領域はそのリスクが許容されるわけではないが、ある便益...
許容可能なレベルに達していないリスクから優先的に対策を実...
事故の発生件数に重きを置くのではなく、重大事故の防止を徹...
リスクアセスメントのポイントは「事前にどれだけのリスクを...
日本では想定したリスクにすべて対策を講じるという考えが強...
この「すべて対策を講じる」という考え方が弊害となって日本...
リスクアセスメントで大切な要素は以下の２つであると考える。
①我慢できる領域と受け入れられない領域の線引きをしっかりと...
②受け入れ不可能な領域のリスクに対しては優先的に対策を実施...
すなわち、これは諦める・これは絶対に阻止するという取捨選...
↑この考え方は日本人には受け入れずらいため、なかなかリスク...
リスクアセスメントは重大なリスクに対処するための合理的な...
また、過去に起こった良い経験も悪い経験も同等にフィードバ...
具体的には設計の初期段階で、
＝通常の使用＋事故＋誤使用＋早期劣化＋etc..
ただ、懸念すべきことはリスクアセスメントを考えれば考える...
起きるかどうかわからないことにお金や時間をかけることは、...
これもリスクアセスメントの実現が難しい理由に当たると思う。
***IHIの方の講座(5/31) [#mb17c552]
現代では建設において、ＣＯ２の削減が優先である.
ＣＯ２の排出量はアメリカは横ばい、ＥＵは下がっている.
今の日本では【１橋梁/104人】を支えている.
土木の現場では、ＣＯ２排出量は コンクリート＞鉄筋＞桁 ...
ＣＯ２を減らす具体例としては、耐候性を用いること・木材を...
ヨーロッパでは、工事において重要視されているのは費用より...
しかし日本は人口が多いため、現時点ではＣＯ２を考えること...
今の時代に沿った新しい優先すべきことが出てきたと考える.（...
***橋梁と火災の関係 [#d0651243]
鋼橋にまつわる火災事故は、国内外問わず数多く報告されてい...
しかし実際に発生することは稀であるため、その対策をしてい...
よって、いざ火災事故が発生してしまった場合に、利用者の安...
そこで、土木学会鋼構造委員会では，2012年6月に「火災を受け...
ここでは、火災による熱の影響を受けた材料の特性はどのよう...
供用中及び工事中に考えられる火災原因を列挙する。
供用中…不審火・走行車両（タンクローリー等）の横転及び炎上...
工事中…動力工事・溶断及び溶接による火の粉の飛散・塗料作業...
材料自体が出火原因となる事例はほとんどない。
また。延焼に必要不可欠なものは①火源、②燃えるもの、③酸素の...
延焼要因は、放射熱・接炎・熱気流・火の粉などがある。
加えて炎が風下側に傾き隣接した建物に直接接炎する場合もあ...
鋼橋は火災の被害を受けると、構造部材は温度上昇に伴う強度...
一般的に鋼材（構造用鋼材）は、受熱時の温度が600℃を超える...
鋼材は合熱処理によって結晶構造が変化するという性質を有し...
よって、火災による鋼材の損傷評価については受熱温度の推定...
このときの鋼橋の補修ですべきことは、塗装の補修・鋼桁の補...
火災発生から補修施工までの流れ👇
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ka...
引用元：[[火災を受けた橋梁の安全性評価と補修＿実例:https:...
主塔は400℃、ケーブルは700℃くらいまで耐えられるそう。
SS材は通常〜300℃くらいまで耐えられる。
SM材は400℃〜500℃程耐えられる。このとき、500℃まで耐えられ...
火災による影響はケーブルよりも主塔にある。
そのため主塔が受熱する温度や時間によって、耐力や強度、変...
例）主塔が崩壊したときのケーブルの切れ方・主塔の耐力が落...
健全度評価フロー👇
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/fr...
鋼材部の被災目安👇
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ko...
引用元：[[火災を受けた橋梁の健全性評価に関する試験調査:ht...
鋼材の温度が高くなると、柔らかくなるすなわちヤング率は小...
火災による鋼材の変質は目視試験での確認は困難である。
また、実鋼材を抜き取り機械試験などにより確認することは通...
火災の被害を受けた鋼部材における補修を行うべき損傷は主に...
車両が燃える原因として、燃料やオイルなど着火物が漏れたり...
***運転者の心理について [#g0482488]
「人」が原因となって起こる交通事故＝”ヒューマンエラー”
エラーの種類として、思い込み・知識や技能不足・意図的ルー...
また、急いでいるときには必然的に事故を起こしやすい.
疲労・ストレスが交通事故へつながることも多々.
事故を起こした直後の運転者は二次被害につなげないように、...
その行動が橋梁上では、主塔やケーブルに近づくことを意味する.
もし、追突や衝突で車両が爆発した上でこの行動を取ったとし...
***交通事故の原因について [#rdbc3731]
令和４年における交通事故の原因ランキングを[[交通事故の原...
１位 安全不確認
２位 脇見運転
３位 動静不注意（他の車や歩行者の存在には気づいているも...
４位 漠然運転（ぼんやりと運転をしている様）
５位 交差点安全進行
６位 運転操作不適（危険は認識しているものの、それに対す...
７位 一時不停止
８位 信号無視
９位 歩行者妨害等（横断歩道を渡っている、または渡ろうと...
１０位 優先通行妨害（信号機がない交差点で、標識や道路幅...
近年は交通事故数が減少している.
これの主な要因を内閣府は以下のように示している.
１，道路交通環境の整備
２，交通安全思想の普及徹底
３，安全運転の確保、車両の安全性の確保
４，道路交通秩序の維持
５，救助・救急体制等の整備等
６，交通安全基本計画に基づく諸対策を総合的に推進してきた...
単路や橋での事故として多いものは[[様々な交通事故統計:http...
１，車両相互追突
２，車両相互正面衝突
３，車両相互その他
４，車両単独工作物衝突
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ch...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ca...
***リスクアセスメントの質問 [#af354ae0]
・リスクの確率の計算方法（事故の確率や火災の確率）及び計...
・死者数や被害者数の想定方法（人口に対して何割か）
・火災の起こる確率の式は何を使うか
・主塔への熱の伝達方法（動的解析か静的解析か）※今は静的解...
→毎秒もしくは毎分、何度熱が加わりそれがどのくらい続...
→ある特定の温度のときに主塔がどのように変形し、かつ...
・由利本荘市や市民に対する経済への影響や交通への影響を考...
→橋梁火災の規模やその発生原因、被害者数、被害状況、...
***鋼構造を語る会（9/20） [#mf556213]
今日は鋼構造を語る会に参加した。
他の大学の人の研究を聞く機会がなかったため、とても面白か...
現在、二酸化炭素を排出しない「グリーンスチール」というも...
最後のゲームでは３位になった！嬉しい！
***青木ゼミ(9/24) [#k361f058]
今日は青木班のみんなと話し合いをした。
【今後の方向性】
デッキ（死荷重・活荷重）半分をつけて同様の解析
火災が起きる→この程度の事故率や火災率が考えられる→火災で...
【リスクアセスメントの浸透】
安全対策への考え方の違い
日本では事故を起こすことが悪だが、海外では事故が起こった...
***中間発表(10/3) [#w4d8d623]
夏休みが終わり中間発表を行った．
スライドと概要に思ったよりも時間がかかったため，もう少し...
発表は時間内に終わらなかったため，もっと練習が必要である．
やっとリスクアセスメントについて皆にまとめたものを見せる...
【質問】
・リスクアセスメントの優先順位等はあるのか（ケーブルが重...
・ケーブル腐食などもリスクアセスメントに含まれるのか
・由利橋に対してリスクアセスメントをしたいのか由利橋を通...
・プロパンガスからの火災は爆発ではないのか
***IHIの方とリスクアセスメントもついて話し合い(10/28) [#o...
もっと「定量的」にリスクを設定していくことが必要！！
そのために・・・
世界や日本での事故率・危険車両の通過割合といった必要なデ...
リスクアセスメントは、過去に起こったデータを用いて予想さ...
適当な過程を自分でつけていくことが必要となる
解析では、由利橋が通行止めになった場合にどの程度困るのか...
ケーブルが火災にあった場合もやってみたいかも。でも、ヤン...
最終的な目標は、そのリスクが許容できるor許容できないかを...
許容できない場合は許容できるまでリスクに対してお金をかけ...
耐火塗装や消化器の設置などなど
わざわざ斜張橋にしたことによって新たなリスクが増えたので...
メンテナンス費用の増大や火災が起こりやすくなった、主塔か...
将来的には一般橋梁と斜張橋との比較とかもいいかも
一般橋梁と由利橋とではリスクが起こる確率が高いor低い
***yec香港(11/30〜12/1) [#zf3ed9cd]
貴重な経験となった４日間になった
色んな人の発表や意見に触れてとても勉強になった
楽しい時間を過ごすことができた！
他のリスクアセスメントの研究をやっている人は実際にリスク...
多くの人は「今はこう」だから「研究結果からこの部分を改善...
質疑応答がやはり反省点。聞かれたことに対して上手く答えら...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/IA...
↑yecで使用したスライド
[質問されたこと]
・活荷重が大きすぎるのではないか
・リスクアセスメントで使用しているデータはどこから入手し...
・リスクは橋梁に対応しているものか
・なぜ直線道路での事故率を考慮する必要があるのか
***今後の方針(12/11) [#qfa15d48]
本研究では「通行止めリスク」を考える
その要因を火災とし、小型車・大型車・特殊車で火災の特徴（...
火災解析は主塔正面・側面で起こったときとケーブルC1,C6で起...
[通行止めリスク]は以下の通りにする
・橋梁上の事故率
・秋田県車両火災率
・大型車混入率
・主塔周りで火災が起こる確率 ←自分で考える
・ケーブル周りで火災が起こる確率 ←自分で考える
データは過去３〜５年分集める
***中間発表(12/23) [#q013a558]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/70...
上記が使用したスライド
[質問されたこと]
・なぜスパン長をかけるのかがわからない
・もっと秋田県の事故率や凍結による事故率など秋田により沿...
・使用した車両火災率は、事故によるものなのか
・なぜ火災がケーブルに影響しなかったのか、何℃くらいになっ...
***卒論発表練習(2/6) [#j5e8bb22]
概ねうまくできた．
話したいことがたくさんで時間が足りない．
もっと要点をまとめた方がいいかもしれない．
最後のまとめは，もう少し全体的なまとめと具体例を提案した...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/02...
***卒論発表会(2/18～2/19) [#c8e0d4c4]
卒論発表が終わった．
練習では概ね上手くできていたと思っていたが，質疑応答は上...
でも，改めて勉強になることが多く，気づきもたくさんあった．
この時期に卒論発表をすることができてよかった思う．
今後はこれからの目標をより明確にする必要があると感じた．...
加えてやはり許容上限と下限の算出方法・そのリスクを選定し...
今はリスクアセスメントの中身（リスクイベント・リスクを計...
今後はもっとおおもとの枠組みを把握する必要がある
例えば，ユーロコードの内容・許容値の設定方法・実際にどの...
日本の通行止め日数やユーザーリスクを調べてみる，
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/so...
*サロメモ [#v109bbb0]
AsterStudyの設定翻訳
,メッシュ,Read a mesh, メッシュの作成をしたモデルを適用し...
,Model Definition,Assign finite element, モデルの要素を設...
,Material ,Definr a material ,各材料物性値の設定,(複数設...
,Material, Assign a material ,設定した物性値を材料に適用,
,BC and load , Assign mechabucal load , 拘束条件 荷重設...
,Analysis , static mechabucal analysis ,線形解析の設定 ,...
,Post Processing , CALC_CHAMP ,ここは出力でなんの要素を出...
,out put , set output result ,出力ファイルの設定を行う。,
[メッシュ]
>LIRE_MAILLAGE(メッシュを読む)>FORNAT()
>ASSE_MAILLAGE(メッシュ組み立て)
>DEFI_GROUP()
>MODI_MAILLAGE(メッシュ方法)
>CREA_MAILLAGE(メッシュ決め)
[Model Definition]
>AFFE_MODELE(モデル割当)>AFFE(割当)>TOUT(すべて)>PHENOMEN...
>AFFE_CARA_ELEM()
>DEFI_BASE_REDUITE()
>DEFI_DOMAINE_REDUIT()
>DEFI_GEOM_FIBRE()
>MACR_CARA_POUTRE()
>MODI_MODELE()
[Material]
>AFFE_MATERIAU(材料割当)>MODELE>AFFE>TOUT(すべて) ←AFFEで...
>AFFE_MATERIAU>AFFE_VARC()>NOM_VARC()→TEMP(温度)>VALE_REF...
>CREA_LIB_MFRONT
>DEFI_COMPOSITE
>DEFI_GLRC>
>DEFI_MATER_GC
>DEFI_MATERIAU(材料特性)>ELAS(ヤング率等)>E(ヤング率),NU(...
>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(線載荷)orFOR...
>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(構成...
NEWTON(ニュートン法の設定)>MATRICE:TANGENTE(接線剛性)
CONVERGENCE(収束判定の設定)>RESI_GLOB_RELA(相対誤差)>ITER...
>Thermics(熱力学)
>THER_LINEAIRE(温度＿線形)
>THER_NON_LINE(非線形)
>THER_NON_LINE_MO(非線形モデル？モード？)
[Post Processing]
>CALC_CHAMP(調べたい分野)>CONTRAINTE(対照？)>SIGM(6成分応...
>CALC_CHAMP(調べたい分野)>DEFORMATION(力？)>DEGE(),EPME()...
>CALC_CHAMP(調べたい分野)>CRITERES(基準)>EPEQ(),EPGQ(),EP...
NOEUは節点(NOEUd(フランス語))、ELGAはガウスの積分点(ELeme...
>RESULTAT>MECA_STATIQUE
>CREA_CHAMP>MODEL→AFFE_MODEL>TYPE_CHAM→NOEU_TEMP_R>OPERAT...
[Output]
>IMPR_RESU(結果照会)>RESU(結果)>RESULTAT(結果)>NOM_CHAM(...
*Salome-Meca2020_熱解析 [#je8a4800]
○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:tem...
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...
FORMAT>MED
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/re...
Astersutudyはどこかで設定してから出ないと選択できないもの...
○Paravis
左側のパイプのブラウザーのbuiltinを右クリック>Open>.medフ...
フィルター(画面上の黒いところ、ファイルとか編集とかあると...
フィルター>Recent>ExtractGroup>SIGM_NOEU(節点の応力(6成分))
*Salome-Meca2020_弾塑性解析 [#z08773c9]
参考[[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>SIM...
NEWTON>MATRICE(マトリクス)>TANGENTE(タンジェント)orELASTI...
CONVERGENCE>RESI_GLOB_RELA>"1.0e-06"任意で良い
CONVERGENCE>ITER_GLOB_MAXI>"100"任意で良い
*Salome-Meca2020_弾塑性解析2 [#qd419800]
○Geometry
いつも通り作成
○Mesh
いつも通り作成
○AsterStudy
右クリックAddstage
[[Function and Lists]]
-DEFI_LIST_REEL(時間間隔の設定)
DEBUT>"0">
INTERVALLE>"1.0","Step length","0.05"入力
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/li...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/in...
-DEFI_LIST_INST(解析オプション)
DEFI_LIST>LIST_INST:listr>
METHODE>MANUEL
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ti...
-DEFI_FONCTION
名前:func>VALE>"0","0"と"1","1"入力>
NOM_PARA>INST
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/fu...
-DEFI_FONCTION
名前:koufuku>VALE>"390","400"と"235","236"入力>
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ko...
NOM_PARA>TEMP>
NOM_RESU>SIGM
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ko...
[[Material]]
-DEFI_MATERIAU
ECRO_LINE>D_SIGM_EPSI:200入力>SY:235
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ec...
-ELAS
E>NU>RHO>ALPHA
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/el...
[[Analysis]]
-STAT_NON_LINE
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/re...
MODELE:model>
CHAM_MATER:fieldmat>
EXCIT>CHARGE:lode>FONC_MULT:func>
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ex...
INCREMENT:times>
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/in...
COMPORTEMENT>TOUT>RELATION>VMIS_ISOT_TRAC>DEFORMATION>PET...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/co...
NEWTON>MATRICE(マトリクス)>TANGENTE(タンジェント)orELASTI...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/ne...
CONVERGENCE>RESI_GLOB_RELA>"0.0001"任意で良い
CONVERGENCE>ITER_GLOB_MAXI>"1000"任意で良い
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/co...
*Salome-Meca2020_モデル作成（Fuseする場合） [#l44a64aa]
○Geometry
1, BOX-1,2,3・・作成
2, Fuse-3objects選択
3, Partitionを作成、このときFuse/BOX-1,2,3を選択
4, Partitionのグループ作成、kotei/kanetuなど
5, Meshを作成
6, ジオメトリのグループ作成
以上
*研究メモ [#i24c5d3a]
-Windowsからstr/k2サーバーへは”FFFTP”を使用する.このとき...
-要素の種類…ビーム＝線、シェル＝面（厚さなし）、ソリッド...
-モデルは半分で良い
-主桁と主塔は節点は重なっても良いが、境界条件は別々で考え...
-マイダスの境界条件設定は「単一選択」で行う
-マイダスでの由利橋モデルにおける等分布荷重解析を行う際は...
-主塔水平材とは、両主塔を固定するための部材であり幅員に対...
地震時のみに機能を発揮し、２本の主塔が個別振動することを...
また、軽微な損傷を許容することで主塔柱の発生断面力を制御...
-プレテンション(引張力＋)とは、ワイヤロープに一定荷重を一...
初期伸びとロープ径の減少が少なくなる・伸びが少なくなる(弾...
-曲げ剛性とは部材の曲げやすさを表す値である.曲げ剛性（Ｅ...
また、曲げ剛性(EI)と曲率半径(r)は比例関係にあり、曲げモー...
-部材内部に生じる応力＝軸力・せん断力・曲げモーメント
-軸力…部材の材軸方向に生じ、材を伸縮させようとする一対の...
-せん断力…部材に「ずれ」の変形を生じさせようとする一対の...
-曲げモーメント…部材を曲げようとする一対のモーメントのこ...
-合成桁とは、床版と鋼桁とを"ずれ止め"で接合し合成床版と同...
-普通鋼とは、炭素（C）・ケイ素（Si）・マンガン（Mn）・リ...
-一方で、SM材はSはSteel（鋼）、Mは船舶を意味するMarineを...
--熱膨張係数とは、線膨張係数αもしくは体膨張係数βから求め...
--主桁の数え方は左からG1,G2…とする
--マイダスで荷重組み合わせの生成方法は、[[マイダス〜解析...
の９－１９ページに記載
--解析結果の変位などを見るやり方は、結果テーブル＞見たい...
--マイダスで変位等をエクセルに書き出しができるため、そこ...
これが写真→→→http://www.str.ce.akita-u.ac.jp/~gotouhan/j2...
--温度解析はサロメでやっても良い
--由利橋モデルは、桁・主塔は梁要素、ケーブルはトラス要素...
--梁要素とは、有限要素の一種であり、2点を結ぶ細長い軸形状...
--トラス要素とは、圧縮／引張の軸力のみを伝える部材のこと
--弾塑性とは、降伏応力以下の力を加えた場合は力を除くと元...
--マイダスでは３Ｄモデルは”シェル”が該当する
--粘性係数とは、流体の粘り気のこと、鋼材の動粘度ν＝0.5e-2...
--応力-ひずみ線図 http://www.str.ce.akita-u.ac.jp/...
--座屈とは、ある一定の荷重を超えるとき柱は急激に折れ曲が...
--主要地方道は高速道路及び一般国道と一体になって広域交通...
--都市ガスは地下のガスパイプ通って各家庭へ
--街道は一般国道の分類（羽州浜街道：うしゅうはまかいどう）
--12時間交通量は午前7時〜午後7時までの交通量。24時間交通...
--軽トラックは車高1.95m、積載高さ2.5m以下という規定。
--交通量が多い目安：4000台/日〜、大型車混入率目安：地方15％
--大型車：バス・トラック・ダンプカー・ミキサー車など
--小型貨物車：軽トラック・小型トラックなど（乗員定員10人...
--タンクローリーとは、ガソリンや水といった液体物や高圧ガ...
--圧力荷重とは，、面上の特定の種類の静荷重のこと。 面上の...
--単路と直線は、道路の形状を指す言葉で、単路は交差点付近...
--一般単路とは、交差点や踏切、一般交通の場所、トンネル、...
*研究参考文献 [#b552c166]
-[[色々なヤング率:http://www.mh.rgr.jp/memo/mq0100.htm]]
-[[断面二次モーメント計算ソフト:https://d-engineer.com/un...
-[[マイダス〜構造編〜:https://in.midasit.co.jp/wp-content...
-[[マイダス〜解析編〜:https://in.midasit.co.jp/wp-content...
-[[マイダスマニュアル:https://down.midasuser.com/UploadFi...
-[[ケーブルの基本理論:http://library.jsce.or.jp/Image_DB/...
-[[由利本荘市ハザードマップ:https://www.city.yurihonjo.lg...
-[[日本橋梁建設協会:https://www.jasbc.or.jp/panfuretto/pa...
-[[鋼橋の補修・補強事例:https://www.jasbc.or.jp/whatsnew/...
-[[鋼橋（上部工）の損傷事例:https://www.mlit.go.jp/road/s...
-[[鋼橋の腐食事例:https://www.jstage.jst.go.jp/article/js...
-[[斜張橋ケーブル腐食を考慮した疲労寿命:https://www.str.c...
-[[斜張橋ケーブルとアンカーについて:https://www.jpci.or.j...
-[[斜張橋ケーブルの耐久性評価:https://kandoken.jp/huroku/...
-[[橋梁の構造＆種類:https://chansato.com/doboku/bridge-st...
-[[鋼橋工事における火災リスクと防火対策:https://www.nilim...
-[[火災時における耐荷力に関する検討:http://library.jsce.o...
-[[火災を受けた橋梁の損傷と対応:https://www.cbr.mlit.go.j...
-[[火災を受けた橋梁の安全性評価と補修＿実例:https://www.h...
-[[火災を受けた橋梁の健全性評価に関する試験調査:https://w...
-[[鋼材の高温時における機械的特性:https://www.jstage.jst....
-[[SM400材について:http://www.ishiharashouji.jp/pdf/QA-Ve...
-[[普通鋼について:http://www.iri.pref.niigata.jp/25new32....
-[[SM材の耐熱能力:https://higa-metal.com/articles/793/]]
-[[橋梁用厚鋼板:https://www.jfe-steel.co.jp/products/atui...
-[[鋼材の材料定数:https://www.structure.jp/databook/data1...
-[[交通事故の原因ランキング:https://agoora.co.jp/jiko/kno...
-[[様々な交通事故統計:https://www.ms-ins.com/kurumamo/ent...
-[[弾性変形と塑性変形:https://www.monodukuri.com/gihou/ar...
-[[トレス・ミーゼス応力の降伏条件:https://mech-hm.eng.hok...
-[[トレス・ミーゼスの違い:https://planteng-tree.com/yield...
-[[ミーゼス応力の定義:https://d-monoweb.com/blog/evaluate...
-[[モール・クーロンの破壊基準:https://www.jseg.or.jp/chus...
-[[面外と面内:https://bonperson-civil.com/in-out-plane/]]
-[[様々な車高:https://suumo.jp/article/oyakudachi/oyaku/i...
-[[鋼材の粘度:https://www.daimaru-kogyo.co.jp/skillinfo.p...
-[[ヒルの降伏応力:https://www.jsme.or.jp/jsme-medwiki/dok...
-[[円筒座標系:https://www.cradle.co.jp/glossary/ja_A/deta...
-[[直交座標系:https://www.cradle.co.jp/glossary/ja_T/deta...
-[[テンソルとは:https://3.1415.jp/35ss9vvx]]
-[[火災を受けた橋梁の安全性評価:https://www.oit.ac.jp/jap...
-[[火災を受けた鋼橋の受熱温度推定に関する解析的研究:https...
-[[火災を受けた鋼橋の損傷調査とその対応について:https://w...
-[[熱分析:https://www.spsj.or.jp/equipment/news/news_deta...
-[[各材料のメリット・デメリット:https://skyciv.com/ja/tec...
-[[斜張橋の設計:http://library.jsce.or.jp/Image_DB/commit...
-[[SM400Bの熱による耐力:https://jopss.jaea.go.jp/pdfdata/...
-[[鉄鋼材料の各温度における許容引張応力:https://www.meti....
-[[高強度鋼:https://www.jisf.or.jp/business/tech/bridge/h...
-[[耐火性能:https://www.jisf.or.jp/business/tech/school/c...
-[[橋梁設計における手続きと照査:https://www.pref.aichi.jp...
-[[タンクローリー車のサイズ:https://driversjob.jp/content...
-[[橋の設計の基本と基準の変遷:https://www.nilim.go.jp/lab...
-[[海外における耐震設計基準と橋梁の比較設計:https://www.j...
-[[橋梁火災事例:https://www.nilim.go.jp/lab/bcg/siryou/tn...
-[[火災時における鋼合成桁の終局耐力と耐火対策に関する研究...
-[[平面応力要素とは:http://www.cae21.org/kansai_cae/Qa/Q0...
-[[板要素を利用したFEM解析:https://in.midasit.co.jp/wp-co...
-[[midasフリーズ解決方法:https://tech.midasuser.com/jp/su...
-[[midas自重:http://manual.midasuser.com/JP/iGen/835/Usin...
-[[移動荷重:https://support.midasuser.com/hc/ja/articles/...
-[[midas特殊解析:https://in.midasit.co.jp/wp-content/uplo...
-[[道路橋示方書:https://www.mlit.go.jp/road/sign/kijyun/p...
*サロメページ [#oda4d11a]
-[[Salome-Meca演習_弾塑性解析(2021)]]
-[[Salome-Meca_熱応力メモ]]
-[[Salome-Meca2017]]
-[[Salome-Meca例題ファイル]]
-[[Salomeメモ]]←ここから線膨張係数入れられそう！
-[[座屈解析]]
*春休み課題の進捗 [#nd67ffe7]
-一端固定他端単純支持梁
--モデル作成について
ローラー支持は「線固定」→線の端点を入力して、線...
固定端は「面固定」→パーティションからグループの...
Mesh→ジオメトリは「Partition-1」を必ず選択
Mesh切る→ジオメトリのグループ作成
☆☆AsterStudyは簡単にできる！！→Operations→Addstage ...
☆☆（E:206000,ν:0.3,ρ:7.8e-9）
☆☆ファイル名には「.med」！！
--パラビィスのやり方
builtinでopen→「Mode」に変更
フィルター→「Mechanics」から「Normal modes anim...
**やること [#i3d95498]
--理論値を求める・・済（ねじれは分からない）←分かった！！
--メッシュ1.0のharusindou_1の解析結果を検討・まとめ・・済
--上記の相対誤差を求める・・済
--メッシュ0.8,0.9,2.0,3.0を解析→比較・検討・・済
→でも、全ての結果を比較しても ①要素数と相対...
→よって「梁形態」が間違えていると考える
☆解決策 ①梁形態を変更（ピンを端っこから10mm...
②ピンを２コに変更（上下...
③ローラー固定で考えてみ...
--どの梁形態を基本にするか決定する・・済 →☆☆「ローラー」
--上記の決定した一端固定他端単純支持梁のメッシュ0.7~4.0を...
--上記の相対誤差の計算・・済
--ローラー上下２個の一端固定他端単純支持梁のメッシュ0.7~4...
--ローラー下１個の一端固定他端単純支持梁のメッシュ1.0~10....
→！！メッシュ1.0...
--上記の解析結果からねじれについても検討・考察・・済
--ローラー下１個の100mmの一端固定他端単純支持梁のメッシュ...
--ねじれをどうするか話し合う→とりあえず後回し・・解決
--グラフ作成（９枚）・・済
--スライド作成・・済
--グラフの確認・使うデータ選別・・済
--スライド作成（長さ・振動の様子の２種類）・・済
--収束値との相対誤差-要素数のグラフ作成→細長比で比較して...
**疑問・考察 [#zaa418b3]
-２次要素にすると鉛直・水平・ねじれモードは収束した。１次...
-１次要素だと水平モードの方が精度が良いが、２次要素だと鉛...
-１次要素だと高次数の方が精度が良いが、２次要素だと低次数...
-細長いとねじれが出にくい。→100mmだと固有振動数１０コでね...
-長さのみ・長さ＋幅・長さ＋高さを変えた条件でそれぞれ解析...
→やはり、１次要素だとどれも水平モードの方が精度が良かった。
→鉛直・水平・ねじれを総合的に見ると長さ＋幅が一番いいかも...
春課題を通して、理論値への収束の条件や様子、理論値に対す...
既に答えの出ている課題であると説明を受けていたが、どれも...
また、１つの知りたいことのために解析や考察を行ってもまた...
ただ結果を出すだけでなく、そこから何が考えられるか、何が...
*３年創造工房 [#d09b48c4]
**2023.10.27 [#t59452ec]
・タッチタイプの練習
・コマンドの練習
-コマンドの一覧 [#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)...
**2023.11.10 [#uc3066e7]
・コマンドの練習
・gnuplotを使用して.pngで図を作成
-作成した図 [#kf32a2c4]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/bb...
使用データ：
,メッシュ長さ,要素数,変位(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,4518,0.3986,-4.34,山本
,1.5,15433,0.396317756757,-4.9,進藤
,1.6,15900,0.399049,-4.24,河合
,1.8,11677,0.404457,-2.9,山口
,2,10460,0.394818715517,-5.3,進藤
,3,734,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,河合
引用元：岡田の卒論日誌（11/18課題）
**2023.11.17 [#k6d60e45]
・片持ち梁の解析（サロメ）
-サロメメモ [#ddb92132]
Geomety BoXを作成 適用→閉じる kotei/saika
Mesh メッシュ作成 Box1を挿入 アルゴリズム Note ジオメ...
FY=100N/100mm^2=1N/mm^2
・測定値＝平均値ー理論値/理論値✕100(％)
・要素数＝ボリューム
・変位平均値（1.3.5.7）
・相対誤差（測定値）を計算
・縦軸＿変位、横軸＿要素数のグラフを作成
--片持ち測定値 [#xcffb3ee]
,メッシュ長さ,要素数,先端変位(4隅の平均値)[mm],相対誤差($...
,0.7,198464,6.54281,1.91,安藤
,0.8,113812,6.5104,2.39,安藤
,0.9,40280,6.3631525,4.60,兼田
,1.1,30055,6.3363525,5.00,兼田
,1.2,26467,6.3043375,5.48,柴田
,1.3,25180,6.304355,5.48,柴田
,1.4,32212,6.31612,5.31,佐藤
,1.5,17753,6.1209,8.23,佐藤
,1.6,14296,6.2044625,6.98,皆川
,1.7,13596,6.2156625,6.81,皆川
,1.8,2866,5.737755,13.98,永山
,1.9,6001,5.7263625,14.15,永山
,2,5617,5.6458525,15.355,辻
,3,2309,5.4728755,17.948,辻
,4,617,3.6160575,0.458,服部
,5,494,3.8580375,0.422,服部
,6,581,2.50682,62.416,梶原
,7,133,1.41225,78.827,梶原
,8,78,1.2887175,80.68,工藤
,9,72,1.2879925,80.69,工藤
,10,60,1.14344,82.85,佐々木
,11,65,1.23124,81.154,佐々木
--グラフ（理論値と測定値） [#y04fdb0d]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/yo...
縦軸：変位(mm)、横軸：要素数
理論値：6.67mm
【理論式】
Pl^3/3EI
**2023.11.24 [#f5b0e5b1]
・単純梁（等方性１次）の解析
--単純梁（等方性１次）測定値 [#g80425ab]
,メッシュ長さ,要素数,変位平均[mm],相対誤差($\frac{salome-...
,0.7,145234,0.42248452735,1.388176,安藤
,0.8,142973,0.42257044598,1.408794,安藤
,0.9,91648,0.420437286573,0.897,兼田
,1.1,27160,0.405618939024,2.659,兼田
,1.2,24675,0.404349,2.96,柴田
,1.3,23446,0.404185,3.00,柴田
,1.4,17738,0.398604,4.34,佐藤
,1.5,15438,0.396593,4.83,佐藤
,1.6,16122,0.398212,4.44,皆川
,1.7,12026,0.393411,5.59,皆川
,1.8,11604,0.393668,5.53,永山
,1.9,10391,0.390695,6.24,永山
,2,10921,0.395103,5.18,辻
,3,2328,0.324762,22.06,辻
,4,1500,0.155013,62.80,服部
,5,432,0.065278,84.33,服部
,6,357,0.213062,48.87,梶原
,7,196,0.1019,75.55,梶原
,8,104,0.1158624,72.20,工藤
,9,81,0.1255118,69.88,工藤
,10,78,0.07733,81.44,佐々木
,11,63,0.1999,52.03,佐々木
--グラフ（理論値と測定値） [#cb6b3f6d]
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/1...
縦軸：変位(mm)、横軸：要素数
理論値：0.4167mm
【理論式】
Pl^3/48EI
**2023.11.29 [#g51a87b6]
・単純梁（異方性１次・等方性２次）の解析
＜等方性１次と異方性１次の比較＞ [#he11b38c]
--単純梁（異方性１次）の測定値
,メッシュ長さ,要素数,先端変位,相対誤差,計算者
,0.7,144563,0.505252,2.76,安藤
,0.8,141517,0.504692,2.64,安藤
,0.9,91648,0.502595,2.216,兼田
,1.1,27160,0.489914,0.363,兼田
,1.2,24675,0.487088,0.791,柴田
,1.3,23446,0.4868010,0.995,柴田
,1.4,17738,0.485999,1.16,佐藤
,1.5,15438,0.485180,1.33,佐藤
,1.6,15900,0.483286,1.71,皆川
,1.7,12142,0.477952,2.80,皆川
,1.8,11604,0.482085,1.97,永山
,1.9,10391,0.470887,2.40,永山
,2,10291,0.480910,2.19,辻
,3,2328,0.431937,12.15,辻
,4,1500,0.430156,12.52,服部
,5,432,0.282968,42.45,服部
,6,356,0.3441556,30.00,梶原
,7,196,0.213934,56.49,梶原
,8,104,0.229874,53.25,工藤
,9,81,0.232308,52.75,工藤
,10,78,0.203271,58.65,佐々木
,11,63,0.222316,54.78,佐々木
--比較のグラフ
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/1...
縦軸：変位(mm)、横軸：要素数
理論値：0.4917mm
【理論式：ティモシェンコの理論】
v = Pl^3/48EI ＋ Pl/4KGA
G：せん断弾性係数
K：せん断補正係数
このとき、［Pl/4KGA］= せん断項
＜等方性１次と等方性２次の比較＞ [#of7eeadc]
--単純梁（等方性２次）の測定値
,メッシュ長さ,要素数,先端変位,相対誤差,計算者
,0.7,144563,0.430124,3.22,安藤
,0.8,141517,0.430132,3.22,安藤
,0.9,91648,0.430020,3.197,兼田
,1.1,27160,0.429828,3.151,兼田
,1.2,24675,0.429836,3.15,柴田
,1.3,23446,0.42974,3.13,柴田
,1.4,17738,0.429797,1.3,佐藤
,1.5,15438,0.429958,3.14,佐藤
,1.6,15900,0.429755,3.18,皆川
,1.7,12142,0.429676,3.11,皆川
,1.8,11604,0.429829,3.14,永山
,1.9,10391,0.429684,3.12,永山
,2,10291,0.429620,3.10,辻
,3,2328,0.429169,2.99,辻
,4,1500,0.429254,3.01,服部
,5,432,0.428170,2.75,服部
,6,356,0.428452,2.82,梶原
,7,196,0.42591,2.21,梶原
,8,104,0.426074,2.25,工藤
,9,81,0.425552,2.12,工藤
,10,78,0.488382,17.20,佐々木
,11,63,0.423972,9.0534,佐々木
--比較のグラフ
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/1...
縦軸：変位(mm)、横軸：要素数
理論値：0.4167mm
２次要素を含む。
--メモ
・等方性・・方向によって物体の物理的性質が異ならないこと。
・異方性・・方向によって物体の物理的性質が異なること。
**2023.12.08 [#qa2a664e]
・単純梁(木材＋鋼材)の解析
--単純梁（木材＋鋼材）の測定値 [#d0d51f4f]
,メッシュの長さ ,要素数 ,変位[mm],相対誤差,計算者
,0.7,155192,0.08378905246,15.365,安藤
,0.8,138808,0.08380386491,15.350,安藤
,0.9,82587,0.083707073981,15.45,兼田
,1.1,38671,0.084201207602,14.95,兼田
,1.2,31929,0.083688,15.466,柴田
,1.3,28621,0.083669,15.4857,柴田
,1.4,28854,0.08368,15.47,佐藤
,1.5,20015,0.084052,15.10,佐藤
,1.6,19448,0.0835402938,15.62,皆川
,1.7,13801,0.0834355098,15.72,皆川
,1.8,12528,0.083733,15.42,永山
,1.9,11769,0.083924,15.23,永山
,2,10699,0.084076876559,15.074,辻
,3,3579,0.08414561753,15.004,辻
,4,1628,0.082794,16.37,服部
,5,1016,0.083033,18.89,服部
,6,839,0.082882,16.26,梶原
,7,554,0.080871,18.28,梶原
,8,285,0.079995,19.20,工藤
,9,261,0.078980,20.22,工藤
,10,232,0.081911,17.26,佐々木
,11,208,0.075676,23.56,佐々木
--グラフ（理論値と測定値） [#k36a4b43]
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/kaneta/1...
縦軸：変位(mm)、横軸：要素数
理論値：0.099mm
２次要素を含む。
【理論式】
Pl^3/48EI ＋ Pl/4KGA
このとき、
EI = EI（木） ＋ EI（鋼）
G = G（木）
と考えることとする。
**2023.12.15 [#p21326f8]
・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を右クリックでジオメトリの作成
⑤ AsterStudyのmeshをダブルクリックして、UNITE「Mesh-2」...
⑥ Outputでファイルの名前変更（.med）
⑦ 解析
＜テキストエディターでグラフ作成＞ [#k2d05487]
① テキストエディターに（要素数,変位）の順で入力→このとき...
② 名前を付けてファイル保存
③ コマンド
ls→cd kaneta23→ls→gnuplot→＞plot "ファイル名" with line,"...
④ メニューから、gftpに行きパスワードを入力
⑤ public htlsの中のj2023から、自分のファイルに転送
ページ名:
![[PukiWiki]](image/pukiwiki.png "[PukiWiki]")
