#author("2021-09-20T13:40:28+09:00","default:kouzouken","kouzouken")
#author("2022-12-28T13:55:42+09:00","default:kouzouken","kouzouken")
#contents
#br
#br
*ケーブル疲労を考慮した斜張橋における静的・動的応答解析 [#u604b131]
**直近の課題 [#pe4564e4]
*ケーブル腐食を考慮した斜張橋の終局強度および連鎖崩壊解析 [#u604b131]
**当面の課題 [#pe4564e4]
#br
☐中間発表の資料作成

☐弾塑性、Model-300の線形解析 (21/09/01)~
☐ケーブル破断箇所・腐食箇所による影響~
☐解析時間ステップ,積分スキームの選択~
#br

☑Marc_Mentatで簡易モデル(片持ち梁をケーブルで吊った構造)解析 (21/07/20)~
#br
#br

☑ボックス作成した斜張橋モデルで線形解析できるか確認 (21/06/21)~

☑salomeで斜張橋モデル(2D・3D)を作成する (21/05/21)~

☑ 今までと同様のモデルで弾塑性解析(梁モデル、ソリッドモデル・バイリニア)(21/05/21)~
*Marc_mentat [#i623feff]
**センターケーブルが破断した場合の連鎖崩壊挙動 [#bd398520]

☑ ケーブルの軸応力を求められるようにし、分布荷重モデルの解析が可能であるかの確認。(21/05/14)~
-センターケーブル破断時の変形図~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/k21.png~
#br

☑ 梁モデル(温度変化)のケーブル部分のノードを3つ→2つにする。節点をヒンジ設定にする。(21/05/13)~
-隣接するC13の破断ひずみ~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/k23.png~
#br

☑ salomeでケーブル要素を用いたモデル(簡易的)に対して解析可能であるのかを確認する。(21/04/26)~
-主桁中央部の鉛直変位~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/k24.png~
#br

ケーブルが腐食していれば,あるケーブルが破断した場合,連鎖崩壊が生じる可能性があることを見出した.~
本モデル橋では,アンカーケーブルが最初に破断した場合は,隣接するケーブルが次々と破断し,最終的には~
橋全体が崩落した.一方,最初にセンターケーブルが破断した場合は,隣接するケーブルは破断ひずみに達せず~
連鎖崩壊は生じなかった.~
#br
#br

**ケーブル破断時間による影響 [#a2a203a6]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/k20.png~
C1破断時の主桁鉛直変位の時刻歴応答~
#br

Time 0.01s:ケーブルが瞬間的に破断~
Time 1s:ケーブルが1sかけて破断~
Time 5s:ケーブルが5sかけて破断~
#br
#br
#br

*修論日誌 [#t46fa515]
***9月~ [#sa5d7dab]

**荷重の再分配と衝撃荷重を分けた検討 [#mf4bdcb3]
Wolf and Starossek (2009)が採用したDAFは,~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/k17.png~
#br
, 日付 , 時間帯 , 作業時間 , 内容 , 立会 
,9/1,10:00~18:00, 7h ,  節点温度  ,
,9/2,10:00~18:00, 7h ,  節点温度  ,
,9/3,10:00~18:00, 7h ,  節点温度  ,
,9/6,10:00~18:00, 7h ,  弾塑性・Model-300  ,
,9/7,10:00~18:00, 7h ,  弾塑性・Model-300  ,
,9/8,10:00~18:00, 7h ,  弾塑性・Model-300  ,
,9/10,10:00~18:00, 7h ,  弾塑性・Model-300  ,
,9/13,10:00~18:00, 7h ,  弾塑性・Model-300  ,
,9/14,10:00~18:00, 7h ,  弾塑性・Model-300  ,
,9/15,10:00~18:00, 7h ,  弾塑性・Model-300  ,
,9/16,10:00~18:00, 7h ,  弾塑性・Model-300  ,
,9/17,10:00~18:00, 7h ,  弾塑性・Model-300  ,
#br
青木論文で採用したequivalent DAFは,~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/k18.png~
#br
Sdyn(t):動的応答の瞬間時刻tにおける応答の最大/最小値~
Shealthy:健全橋(全ケーブルあり)の静的解析で得られた応答値(図5(a))~
SF0:破断ケーブルを取り除いた系の静的解析で得られた応答値(図5(b))~
SF1:破断ケーブルの張力を圧縮方向にかけ,静的解析で得られた応答値(図5(c))~

&size(18){9/16(木)~}; &br;~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/k19.png~
#br

中村先生(メモ)~

弾塑性~
テーブル(降伏応力・ひずみの入力)、~
・挙動をより詳細にみる~
固定端が壊れているか、軸応力が発散しているか~
・縁端応力の表示の仕方~
→値を見直して結果の図を送る~


Model-300の線形解析~
・各作用荷重ごと(Pr,D,CW,L)結果(変位・応力等)を確認~
・ステップに分けて(作用荷重)、結果を確認できるかの確認~

斜張橋モデルの断面計算~
・許容応力度法~
道路橋示方書、鋼構造~


荷重分配率load redistribution rate (LRR)は,~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/k22.png~
#br
#br
#br
#br

**ケーブル破断時の主桁および主塔の応力分布の変化に関する検討 [#r55f2b6a]
C1(アンカーケーブル)が破断した場合の,主桁および主塔の~
(a)変位,(b)曲げモーメント,(c)軸力,変化を以下の図にまとめた.~

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/k11.png~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/k12.png~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/k13.png~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/k14.png~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/k15.png~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/k16.png~

&size(18){9/14(火)~}; &br;~


Model-300橋の諸元~
・主スパン長300m~
・主桁断面は鋼製箱桁断面、幅員11m~
・床版は鋼床版を想定、鋼板材質はSM400(降伏応力235MPa)~
・主塔は2本で断面は長方形~

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/s20.png~

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/s23.png~

#br
#br

ケーブルの断面積~
・側径間に7本のケーブル(C1-C7)、主径間の左半分に7本のケーブル(C8-C14)、主径間の右半分および右側径間は左右対称のケーブルを配置~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/s25.png~

**ケーブル構成則の応力をゼロにする手法 (除荷) [#p0b09c5e]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/k01.png~
#br
#br

常時設計荷重~
・橋軸方向単位長さあたりの死荷重は鋼桁鋼桁重量(69.7 kN/m),地覆・高欄(24.5 kN/m),アスファルト舗装(28.4 kN/m)~
・端支点の負反力を抑制するため,側径間にはカウンターウエイト(CW, 60 kN/m)を載荷~
・死荷重作用時に主桁および主塔の曲げモーメントが平滑化かつ最小になるようにケーブル・プレストレス(Pr)を導入した ~
・活荷重(L)は道路橋示方書6)のB 活荷重とした.すなわち,集中荷重  p1(10 kN/m2)を10m  長さおよび幅員 ~
5.5m に作用させる.ただし,その他の部分は  5.0 kN/m2とする.さらに,分布荷重  p2(3.5 kN/m2)を幅員 5.5mに作用させる.~
ただし,その他の部分は  5.0 kN/m2とする.なお,集中荷重はスパン中央に,等分布荷重は全径間に載荷した. ~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/s26.png~
崩壊過程:~
250sでC1が破断した後,発生した衝撃荷重によりC2,C3,C4,C5が破断した(250.96s~251.62s).~
そして,C6,C7が破断し(255.591s~255.716s),左径間側すべてのケーブルが降伏応力に至ると,~
左径間側のC4付近の主桁は鉛直下向き方向に大きくたわみ,中央径間側のケーブルに応力が集中した.~
C4付近の主桁断面が255.803s~256.235sかけて引張・圧縮側に全断面降伏し塑性ヒンジが形成された.~
(そのとき主塔基部および主桁中央部の断面は全断面降伏とはならなかった.)~
#br


梁断面の可視化~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/s21.png~
#br
主塔(箱桁断面)の断面作成~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/s22.png~

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/k02.png~
主桁中央部の鉛直変位~
#br
#br
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/k03.png~
C2応力履歴~
#br

常時設計荷重による変形~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/s24.png~

崩壊時の変形図(解析時間258s)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/k10.png~


#br
#br
#br
#br

**Pr(ケーブル張力)を瞬間的にゼロにする手法 [#p83bd95b]
動的解析を行い,ケーブルが破断した場合の連鎖崩壊を再現した.~
解析法はニューマークβ法,時間ステップは0.01としている.減衰定数は,吊り橋等で用いられているh=0.0032を採用した.~
#br

&size(18){9/7(火)~}; &br;~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/k04.png~

中村先生(メモ)~
節点温度~
・A、B、C全てに‐100℃をかけると、A,B,Cそれぞれ2倍縮んでいる。~
・片側の節点だけだと、半分の値だけしか縮まらない~
・ある要素に温度応力をかけるときは、両側の節点に節点温度を与える必要がある。~
→完了(2節点での温度勾配を考慮するのではないか)~
ケーブルが腐食している状態を再現するため,下記に示すケーブル構成則を全ケーブルに適用させた.~
C1(アンカーケーブル)が破断した場合の連鎖崩壊解析を行った.~

弾塑性~
・はりを10分割(節点)~
・はり断面の確認(U領域等の設定が比率かどうか)~
・ケーブルにも断面形状を与える(ケーブルが先に塑性することもある)~
・ケーブルを薄くしたときの挙動(5㎜のときケーブルが先に塑性する)~
・荷重の組み合わせ(荷重ケース)~
→計算結果を確認
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/k05.png~


Model-300の線形解析~
・主桁の断面の設定(薄肉断面梁、断面形状の設定・4つ、箱桁断面)~
×シェル~
・B点、E点はそれぞれ同一の節点としてもよいが、別々の節点にしてバネで繋いでもよい(鉛直かため、水平・回転緩め)~
→計算結果を確認
#br
#br



-解析結果
#br
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/k06.png~
C1破断直後(解析時間253s)の変形図(scalefactor:5倍)~
#br
#br
C1破断後にC2,C3,C4,C5が第二降伏応力(628MPa)および破断ひずみに達し,310s付近で橋全体が崩壊に至った.~
310s付近で主桁中央部付近および主塔基部の断面が全断面降伏していた.~




&size(18){9/1(水)~}; &br;~

・モデル橋の線形解析~
主径間300mの斜張橋(Model-300)の線形解析~
実橋は二面吊り斜張橋だが、片面のみを考慮し平面モデルとして解析する~

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/s17.png~


#br
#br

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/k07.png~
解析時間310s時の変形図(scalefactor:5倍)~

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/s18.png~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/k08.png~


http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/k09.png~


#br
#br



・荷重を増分荷重にする~
▸テーブル~
Type:time~
Formula:1.0v1~
**縦リブ(Uリブ)を考慮した主桁断面 [#pac4783f]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/p7.png~

Uフランジ:20個(4+6+6+4)~
Lフランジ:20個(4+6+6+4)~
Uリブ1つあたりの断面積50cm²~
(https://www.jfe-kenzai.co.jp/product/06/01/index.html:JFE建材HPの構造・仕様を参考に50cm²とした)~

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/s16.png~



#br
#br


・節点温度の確認~
片持ち梁を2分割し、3つの節点 A,B,C を~

(1)A=-100, B=0, C=0~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/s10.png~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/s11.png~
**3Ⅾモデル(骨組みモデル) [#d845d59d]
ケーブル,主桁および主塔を梁要素で骨組みモデルを作成した.~
#br

-解析モデル(3径間連続鋼斜張橋・1面吊り)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/p1.png~
#br

-主塔側面図~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/p2.png~
#br
(2)A=0, B=-100, C=0~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/s8.png~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/s9.png~


#br
(3)A=0, B=0, C=-100~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/s6.png~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/s7.png~
-変位-曲げモーメント-軸力図(主桁)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/p4.png~

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/p5.png~

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/p6.png~
#br
#br

(4)A=-100, B=-100, C=-100~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/s12.png~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/s13.png~

**連鎖崩壊解析(動的応答) [#t0c7c7c1]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/dyn10.png~

#br
(5)A=0, B=-100, C=0 、左側の要素に熱膨張係数を入れない~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/s14.png~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/s15.png~
**幾何学非線形性の考慮 [#r3fa5fed]
Model300におけるC1-C2-C3破断時のモデルで解析した.解析時間は120秒とし,減衰定数は0.02を用いた(レイリー減衰).~
解法はニューマークβ法である.荷重は以下のように設定した.~

・節点温度を与えると、節点に接している要素全体に温度応力がかかる~
・熱膨張係数を与えている部材のみ、温度応力がかかる~
(熱膨張係数が与えていなければ節点が共有されている要素には温度応力はかからない)~
1.設計荷重を10秒かけて載荷し,そのまま載荷し続ける.~
2.C1のPrを10秒かけて載荷し,20秒時Prを除去~
3.C2のPrを10秒かけて載荷し,30秒時Prを除去~
3.C3のPrを10秒かけて載荷し,40秒時Prを除去~

結果~
C1破断時(20秒)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/dyn09.png~
C2破断時(30秒)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/dyn08.png~
C3破断時(40秒)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/dyn06.png~
120秒~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/dyn07.png~

**静的解析との比較 [#p95b81b3]
設計荷重(k=1)を載荷し,静的解析と動的解析の結果を比較した.微小変位解析として計算を行った.~
動的解析で得られた値は静的解析で得られる結果と同一となることがわかった.~
モデル:Model-300~
荷重:D+CW+Pr+D+L~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/dyn05.png~

**動的応答解析 [#c562fb03]
-時刻歴応答解析(ニューマークβ法)~

↓アニメーション(斜張橋のケーブル破断を想定した時の動的挙動)~
↓~
【C1破断時の挙動(C1:20秒のときPr除去)】~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/dynamicC1.gif~
#br
#br



***8月~ [#bd7c2a0c]


【C2破断時の挙動(C1:20秒,C2:22秒のときPr除去)】~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/dynamicC1-C2.gif~
#br
, 日付 , 時間帯 , 作業時間 , 内容 , 立会 
,8/16,11:00~18:00, 5h ,  Marc  ,
,8/17,11:00~17:00, 5h ,  Marc  ,
,8/18,10:00~18:00, 5h ,  Marc  ,
,8/19,10:00~18:00, 5h ,  Marc  ,
,8/20,10:00~18:00, 5h ,  Marc  ,
,8/23,10:00~18:00, 5h ,  Marc  ,
,8/24,10:00~18:00, 5h ,  Marc  ,
,8/25,10:00~18:00, 5h ,  Marc  ,
,8/26,10:00~18:00, 7h ,  Marc  ,
,8/27,10:00~18:00, 7h ,  Marc  ,
,8/30,10:00~18:00, 7h ,  Marc  ,

 
#br
【C3破断時の挙動(C1:20秒,C2:22,C3:24秒のときPr除去)】~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/dynamicC1-C2-C3.gif~
#br

&size(18){8/30(金)~}; &br;~

中村先生(メモ)~
・境界条件▸節点変位 で本当にケーブルに対して温度応力を与えることができているか~
手計算と比較、片持ち梁で段階的に節点温度を与えて挙動をみる~
・ファイバー要素▸弾塑性解析~
・斜張橋モデル▸終局強度~

・解析後、必ず結果を確認

#br



&size(18){8/26(木)~}; &br;~

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/dyn01.png
設計荷重(10秒かけて載荷)~
#br
結果ファイル▸反力Y~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/sa5.png~

#br
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/dyn02.png
C1破断(20秒のときPr除去)~
#br
結果ファイル▸変位Y~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/sa4.png~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/dyn03.png
C2破断(22秒のときPr除去)~
#br
いろいろと試してみて、~
境界条件▸新規(状態変数)▸節点温度(-200℃)▸2要素の結合部に設定で~
それっぽい形になった
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/dyn04.png
C3破断(24秒のときPr除去)~
#br
**Model-600 [#wf6637cf]
主径間長600m、 900mの鋼斜張橋に関し、以下の手順に似て崩壊解析を行う。~
Part 1: 常時設計荷重を満足するようなモデル橋を設計する。解析は弾性とし、断面照査は許容応力法を用いる。~
Part 2: モデル橋の終局強度および崩壊形態を明らかにする。解析は弾塑性解析とする。~

#br
#br

Part 1 線形解析~
1. モデル橋の諸元および部材断面~
主径間600mの斜張橋( Model-600)の線形解析を行う。~
実橋は2面吊り斜張橋であるが、片面のみを考慮し、平面モデルとして解析する。~

#br
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/M600/m01.png~
#br

主桁断面~
UFlg. 3,000 x 25~
LFlg. 3,000 x 25~
2-Web 4,000 x 25~
Steel Grade: SM490 (Yield stress: 315 MPa)~

&size(18){8/23(月)~}; &br;~
#br
主塔断面~
UFlg. 4000 x 25~
MFlg. 4000 x 25~
LFlg. 4000 x 25~
2-Web 5000 x 25~
Steel Grade: SM490Y(Yield stress: 355 MPa)~

・初期条件▸新規(構造)▸仮想温度
#br
平行線ケーブル ~
ヤング係数 1.95 E5 N/mm2~
材質 ST1570~

3. 解析結果~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/M600/m02.PNG~
設計荷重載荷時の変形~

4. モデル橋の諸元および部材断面~
主径間900mの斜張橋( Model-900)の線形解析を行う。~
実橋は2面吊り斜張橋であるが、片面のみを考慮し、平面モデルとして解析する。~

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/ka1.png~




&size(18){8/19(木)}; &br;~

中村先生~
・ケーブル要素は用いない(高度な計算になり複雑になるため)~
・代わりに梁要素に温度荷重を与えてケーブルのプレストレスとする~
・ソリッドでの全景モデルで構造計算することは、かなり難しい~
・骨組みモデルで計算する~
#br
#br


**連鎖崩壊解析(静的) [#b89221fc]
-亜鉛メッキ鋼線の腐食~
中村らによるケーブルを構成する亜鉛メッキ鋼線の腐食に関する研究.亜鉛メッキ鋼線の腐食を3つに~
分類.腐食レベル1は亜鉛メッキ層のみ腐食しており地鉄部は腐食しておらず強度は健全鋼線と同様.~
腐食レベル2は,亜鉛メッキ層の下の地鉄部の一部が腐食している状態.腐食レベル3は,地鉄部の腐食~
がさらに深く進行した状態である.~
#br
#br


&size(18){~8/18(木)}; &br;~
ケーブルの腐食が進行すると孔食のみならず,広範囲にわたり鋼は腐食し、鋼断面積も減少する.~
(腐食レベル3に着目し)ケーブル断面が20%減少するとして,全長に渡り一様にケーブル断面積~
を減少させてケーブル破断を想定したモデルに対し,荷重を作用させ,部材ごとの応力が降伏応力に達する~
か否かの照査を行った.~


Mark_Mentat ユーザーズガイド(例題)を実際に解いてみた~
目的:Markでの構造解析のやり方を一通り学ぶため
1)アンカーケーブル(C1)が破断したケース~

(a)C1が破断~

#1Summary Comprehensive Sample Session~
「3D linear elastic analysis」 

Analysis type : 荷重ステップを変えた静的解析~
Element type: Solid element type 7~

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/sa2.png~
スパン中央の鉛直変位:2.105mm~
隣接するC2の応力は986N/mm2, ~
C14の応力は502N/mm2, ~
主桁、主塔、ケーブルの応力すべて弾性範囲内~
#br
#br

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/sa1.png~
(b)続いてC2が破断したと想定~

スパン中央の鉛直変位:2,586mm~
隣接するC3の応力は1299N/mm2, ~
C14の応力は316N/mm2, ~
C3はケーブルの第1降伏点に達する.主桁,主塔,の応力は弾性範囲内.~
#br
#br
(c)続いてC3が破断したと想定~


#25 Summary Plastic Limit Load Analysis of a Simple Frame Structure

スパン中央の鉛直変位:3,139mm~
隣接するC4の応力は1581N/mm2, ~
C14の応力は102N/mm2, ~
C4はケーブルの第2降伏点の降伏応力に近い応力が算出.主塔242N/mm2, 降伏応力に達する.~
#br
#br
#br


Analysis type : 非線形静的解析~
Element type: Beam element type 52

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/sa3.png~
**Model-300 [#k5e303bd]
***モデル橋の設計(断面照査) [#qb4f7fc7]




主径間300mの斜張橋(Model-300)の線形解析を行う.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/s23.png~
#br
ケーブルの断面積~
・側径間に7本のケーブル(C1-C7)、主径間の左半分に7本のケーブル(C8-C14)、主径間の右半分および右側径間は左右対称のケーブルを配置~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/s25.png~
#br
支点条件~
A, F: 鉛直固定,橋軸方向自由,回転自由~
C, D: 鉛直固定,橋軸方向固定,回転自由~
B, E: 桁と主塔を結合し,鉛直方向は固定,橋軸方向と回転は自由~
#br
作用荷重~
1) 橋軸方向単位長さあたりの死荷重は鋼桁鋼桁重量(69.7 kN/m),地覆・高欄(24.5 kN/m),アスファルト舗装(28.4 kN/m)~
2) 端支点の負反力を抑制するため,側径間にはカウンターウエイト(CW, 60 kN/m)を載荷~
3) 死荷重作用時に主桁および主塔の曲げモーメントが平滑化かつ最小になるようにケーブル・プレストレス(Pr)を導入した ~
4) 活荷重(L)は道路橋示方書6)のB 活荷重とした.すなわち,集中荷重  p1(10 kN/m2)を10m  長さおよび幅員 ~
5.5m に作用させる.ただし,その他の部分は  5.0 kN/m2とする.さらに,分布荷重  p2(3.5 kN/m2)を幅員 5.5mに作用させる.~
ただし,その他の部分は  5.0 kN/m2とする.なお,集中荷重はスパン中央に,等分布荷重は全径間に載荷した. ~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/s26.png~
#br
#br
#br
<解析結果>~
-各荷重載荷時の変形図~
a) ケーブルプレストレス載荷時(Pr),スパン中央の変位:+702mm ~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/i08.png~
#br
#br
b) 死荷重載荷時,カウンターウエイト含む(D),スパン中央の変位:-1,441mm~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/i09.png~
#br
#br
c) 分布活荷重載荷時(LP2),スパン中央の変位:-487mm~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/i11.png~
#br
#br
d) 集中活荷重載荷時(LP1),スパン中央の変位:-159mm~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/i10.png~
#br
#br
e) Pr + D + LP1 +LP2, スパン中央の変位:-1,384mm~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/i12.png~
#br
#br
#br
-断面照査~
a) 主桁(主塔部,スパン中央)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/i13.png~
#br
#br
b) 主塔(基部)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/i14.png~
#br
#br
c) ケーブル~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/i15.png~
#br
#br


***7月~ [#e466ef41]

***終局強度および崩壊形態 [#f0558fbf]

-主桁,主塔断面はファイバー要素に分割する.~
#br
, 日付 , 時間帯 , 作業時間 , 内容 , 立会 
,7/19,10:00~18:00, 5h , ゼミ,後藤さん 青木さん 石黒さん
,7/20,10:00~18:00, 5h , Marc操作 ,
#br
<解析方法>~
最初にケーブルプレストレス(Pr)を載荷し,死荷重および活荷重を漸増させる.~
すなわち,~
P = Pr + k ( D + LP1+LP2 )  ~
である.kは荷重増加係数であり,0.01ずつモデル橋が崩壊するまで増加させる.~
終局時のkを終局荷重係数kuとする.~
#br
#br
#br
<解析結果>~
#br
-終局時の変形図~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/i16.png~
#br
崩壊形態:~
k=2.64のとき,崩壊した.~
主桁および主塔の変形は,側径間のケーブルC4がk=1.90にて第一降伏点に達するまでは,~
線形増加する.k=2.24で主塔位置の主桁下フランジが,k=2.26にてスパン中央の主桁下フランジが降伏する.~
その後,スパン中央の主桁断面の塑性化が進む.そして,k=2.61で側径間ケーブルC3,C4が第二降伏点に達し,~
k=2.62で,スパン中央の主桁が全断面降伏し塑性ヒンジが形成され,橋全体が崩壊する.~
#br
#br
青線:中村先生 オレンジ線:Marc~
fig.1 kと主桁鉛直変位~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/i17.png~
#br
#br
fig.2 kと主塔水平変位~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/i18.png~
#br
#br
fig.3 kと主桁下縁応力~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/i19.png~
#br
#br
fig.4 kと主塔下縁応力~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/i20.png~
#br
#br
fig.5 kとC1ケーブル応力~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/i21.png~
#br
#br
fig.6 kとC3ケーブル応力~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/i22.png~
#br
#br
fig.7 kとC14ケーブル応力~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/i23.png~
#br
#br
#br




***これまでのまとめ・整理 [#c0ee7060]
<大まかな流れ>・salome_mecaでケーブル(プレストレス)をどう表現するか。~
↓       ・解析可能であるのか、適切な結果が出るのか~
↓ keyword:・ケーブル要素~
↓     ・温度変化~
↓     ・線形、非線形解析~
↓     ・ヒンジ~

***連鎖崩壊形態 [#b086d86f]

ケーブル要素…コンクリートのプレストレスとして用いられる?~
       PCのような緊張力を想定している?(salomeの文献から)~
ケーブル要素×~
→温度変化で緊張力を再現(20℃→200℃)~
・salome_meca梁モデル〇~
・線形解析〇(たわみの値は中村先生と同じ)~
・片持ち梁をソリッド、ケーブルを梁要素(上記と同じ条件)〇~

衝撃荷重を静的荷重の2倍になると想定する.~
減衰のない1自由度のモデルに一定力が突然作用すると応答変位は2倍になる.~
そして,静的解析により1本のケーブルが破断直後,衝撃荷重により連鎖的に~
他のケーブルが破断する様子を解析する.(台湾アーチの連鎖崩壊を解明した方法)~
#br
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2105/scr15.png~
熱応力・梁要素(ケーブル)片持ち梁【10m・集中荷重】たわみ-4.15603mm~
スパン中央のケーブル(C14)が破断したケースと,アンカーケーブル(C1)~
が破断したケースについて,連鎖崩壊過程を検討する.~
#br
#br
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2105/scr16.png~
-C14が破断したケース~

・たわみ0.03373mm~
・鉛直荷重1kN~
・ケーブルの温度変化-200℃~
・ソリッドの断面10cm×10cmの正方形、ケーブル半径10cmの円形~
腐食などの原因によりC14が破断したと仮定する.~
#br
1) C14が破断する前(Pr+D+CW+LP1+LP2が作用する)の状態~

スパン中央の鉛直変位:1,384mm~
全ての主桁,主塔,ケーブルは弾性範囲である.~
C1の応力は573N/mm2, ひずみは0.00174~
C14の応力は526N/mm2, ひずみは0.00150~
なお,健全モデルの終局荷重係数は ku=2.62~
#br
#br
2) C14が破断した直後~

C14が破断すると,C14が無くなった構造系に,死荷重および活荷重が衝撃的に作用する.これを,静的荷重の2倍が作用すると考え,~
2.0 x (D+CW+LP1+LP2) を作用させる.その結果,スパン中央の主桁断面の下フランジが降伏応力に達するが,全塑性断面とはならない.~
他の位置での主桁および主塔は弾性域内である.~
破断したC14に隣接するC13の応力およびひずみは他のケーブルより大きく,その応力は1,387N/mm2, ひずみは0.025となり,~
第一降伏点を超える.ただし,破断ひずみ(0.04)には到達しない.すなわち,全体崩壊にはならない.一般的に斜張橋は不静定次数が~
高く,redundancy に優れるため,1本のケーブルが破断しただけでは,全橋崩壊を惹き起こすことにはならないと言える.~
しかし,ケーブルが腐食している場合には,断面積の減少,じん性の低下(過去の事例では,正常ケーブルの1/3になるとの報告もある)~
の可能性もある.したがって,このC13も破断する危険性がある.~
なお,C14を除いた構造系の終局荷重係数 ku=2.22であり,初期構造系のku=2.62より15%低下する.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/i24.png~
fig1. C14が破断した構造系の終局時の変形~
#br
#br
3) C13が破断した直後~

線形解析→非線形解析(弾塑性解析)~
・(弾塑性)梁モデル→multfiber?(現在×)~
・(弾塑性)ソリッドモデル→〇~
・(弾塑性)ソリッドと梁のモデル→× salome_mecaでは弾塑性解析をするとき、すべての部材(異なる材料があるとき)に弾塑性の設定を与える必要がある(仮説)
C13が腐食していると,断面積が減少またはじん性が低下する.そこで,C14の破断に続きC13が破断する場合を検討する.~
C13およびC14が無くなった構造系に,死荷重および活荷重が衝撃的に作用する.これを,静的荷重の2倍が作用すると考え,~
2.0 x (D+CW+LP1+LP2) を作用させる.その結果,C14付近の主桁の下フランジ,主塔基部下フランジが降伏応力に達するが,~
全塑性断面とはならない. C13に隣接するC12の応力およびひずみは他のケーブルより大きく,その応力は1,570N/mm2で第二降伏点に達し,~
ひずみは0.030となる.したがって,全体崩壊には至らない.~
しかし,その後の荷重ステップでC12は破断ひずみ(0.04)に達し破断する.したがって,静的荷重の2倍の衝撃荷重により破断する可能性は高い.~
なお,C13を除いた構造系の終局荷重係数 ku=1.79であり,初期構造系のku=2.62より32%低下する.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/i25.png~
fig2. C13が破断した構造系の終局時の変形~
#br
#br
4) C12が破断した直後~

斜張橋モデル(ソリッドと梁) 線形解析×~
C13に続き,C12が破断する場合を検討する.~
C12~C14が無くなった構造系に,死荷重および活荷重が衝撃的に作用する.これを,静的荷重の2倍が作用すると考え,2.0 x (D+CW+LP1+LP2)~
を作用させる.その結果,静的荷重2倍の荷重に達する前のk=1.47で,破断したC14付近の主桁断面が全塑性断面となり塑性ヒンジが形成され,~
橋全体が崩壊する.そのときの,C12に隣接するC11のひずみは0.0397である.~
なお,C13を除いた構造系の終局荷重係数 ku=1.47であり,初期構造系のku=2.62より44%低下する.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/i26.png~
fig3. C12が破断した構造系の終局時の変形
#br
#br
#br
-C1が破断したケース~
腐食などの原因によりC1が破断したと仮定する.~
1) C1が破断した直後~

ソリッドモデルの弾塑性解析→結果が中村先生と異なる。ケーブルに曲げが生じている?(構造としてはありえない)ヒンジにできていない?~
C1が無くなった構造系に,死荷重および活荷重が衝撃的に作用する.これを,静的荷重の2倍が作用すると考え,2.0 x (D+CW+LP1+LP2) を作用させる.~
その結果,スパン中央部の主桁,主塔基部断面の下縁側のみ降伏応力に達している.破断したC1に隣接するC2の応力およびひずみは他のケーブルより大きく,~
その応力は1,450N/mm2, ひずみは0.0365となり,第一降伏点を超えるが,破断ひずみ(0.04)には到達しない.すなわち,全体崩壊にはならない.一般的に~
斜張橋は不静定次数が高く,redundancy に優れるため,1本のケーブルが破断しただけでは,全橋崩壊を惹き起こすことにはならないと言える.~
しかし,ケーブルが腐食している場合には,断面積の減少,じん性の低下(過去の事例では,正常ケーブルの1/3になるとの報告もある)の可能性もある.~
したがって,C2も破断する危険性がある.~
なお,C1を除いた構造系の終局荷重係数 ku=2.18であり,初期構造系のku=2.62より17%低下する.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/i27.png~
fig4. C1が破断した構造系の終局時の変形~
#br
#br
2) C2が破断した直後~

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tunoda/depo2.gif~
scalefactor20倍
C2が腐食していると,断面積が減少またはじん性が低下する.そこで,C1の破断に続きC2が破断する場合を検討する.~
C1およびC2が無くなった構造系に,死荷重および活荷重が衝撃的に作用する.これを,静的荷重の2倍が作用すると考え,2.0 x (D+CW+LP1+LP2) を作用させる.~
その結果,C5,C14付近の主桁,主塔基部断面の下縁側が降伏応力に達する.また,C2に隣接するC3の応力およびひずみは他のケーブルより大きく,その応力は~
1,570N/mm2で第二降伏点に達し,ひずみは0.039となる.したがって,全体崩壊には至らない.~
しかし,次の荷重ステップでC5付近の主桁は全塑性断面となり塑性ヒンジが形成される(C3は破断ひずみ(0.04)に達し破断する.)したがって,静的荷重の~
2倍の衝撃荷重により破断する可能性は高い.~
なお,C2を除いた構造系の終局荷重係数 ku=1.74であり,初期構造系のku=2.62より34%低下する.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/i28.png~
fig5. C2が破断した構造系の終局時の変形~
#br
#br
3) C3が破断した直後~

C2に続き,C3が破断する場合を検討する.~
C1~C3が無くなった構造系に,死荷重および活荷重が衝撃的に作用する.これを,静的荷重の2倍が作用すると考え,2.0 x (D+CW+LP1+LP2) を作用させる.~
その結果,k=1.46でC5付近の主桁断面が全塑性断面となり塑性ヒンジが形成され,橋全体が崩壊する. ~
なお,C2を除いた構造系の終局荷重係数 ku=1.46であり,初期構造系のku=2.62より44%低下する.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/i29.png~
fig6. C3が破断した構造系の終局時の変形~

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2106/scr31.png~
斜張橋モデル
#br


#br
#br
#br



&size(18){7/6(火)~}; &br;~
接線同士を重ね合わせてヒンジ?とした。~
**ファイバー要素 [#u0fdf28d]
-断面をメッシュ分割し,各メッシュに応力ひずみ関係を与える.非線形特性をより正確に表現できる.~
軸力変動や二軸曲げによる応力変動を比較的正確に表現できる.~
-解析ソフトによる結果の違いはある程度許容する。
#br
#br




http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tunoda/depo2.gif
**弾塑性解析(材料非線形) [#t09f986d]


http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tunoda/depo2.png


以下の設定が必要になる~
-テーブルでtime,応力ひずみ関係の入力~
-材料特性で降伏点の設定~
-境界条件でtimeを入力~
-荷重ケースで解法の設定~
#br
#br
***テーブルの設定[#rfbf1d93]
1.1 time(荷重を増分荷重にする)~
Type:time~
Formula:1.0v1~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/i01.png
#br
&size(18){7/2(金)~}; &br;~
・ケーブルの断面形状を円柱→四角形~
#br
1.2 応力ひずみ関係
Type:eq_plastic_strain~
Point:任意~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/i02.png~
構成則入力の変更点:~
テーブル(タイプ:eq_plastic_strain)にεの値をそのまま入力するのではなく、~
&size(16){εp=εt-εe}; &br;~εp=εt-εe~
εt:真ひずみ、εe:弾性ひずみ、εp:塑性ひずみ~
の計算式で計算し、テーブルには塑性ひずみの値を入力した。~
以下の図ようなイメージ~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/dansosei/p1.png~
#br
#br
***材料特性で降伏点の設定 [#rbd67574]
-バイリニアの例(主桁)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/i03.png~
-トリリニアの例(斜張橋ケーブル)~
テーブル(応力ひずみ関係)で入力する場合,降伏応力は1(1倍)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/i04.png~
#br
#br
***境界条件[#nc1b4256]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/i07.png~
#br
#br
***荷重ケースの設定 [#w768687b]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/i05.png~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/i06.png~
#br
#br
2.ケーブル・モデル:~
長さ1mのケーブルとする。~
B点に鉛直荷重P = 1.0 kN を増分載荷する。~
断面: 半径 10mm の円形とする。~
構成則に関しては、はりと同じ計算方法で計算を行い入力した。~
解析結果~
1)自由端の鉛直変位と荷重の関係~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/dansosei/p5.png~
#br
2)ひずみと荷重の関係~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/dansosei/p6.png~
#br
3)応力とひずみの関係~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/dansosei/p7.png~
#br
全てのグラフにおいて中村先生の結果と一致した
#br
#br
はり・ケーブル単独の弾塑性解析~
はりモデル:~
・長さ1mの片持ち梁とする。~
・自由端に鉛直荷重P = 1.0 kN を増分載荷する。~
・部材断面
・はりの断面は10 cm x 10 cm の正方形とする。5×5のファイバー要素に分割した。~
・はりはSM400とする。初期ヤング係数は 2.0 x 105 N/mm2 とする。~
#br
#br
その他の設定~
・反復法:完全ニュートンラプソン~
・収束判定 残差と変位(相対許容値):0.001~
#br
解析結果(以下 青線:中村先生の解析結果、オレンジ線:Marcでの結果)~
1)自由端の鉛直変位と荷重の関係~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/dansosei/p2.png~
#br
2)断面応力(Layer1)と荷重の関係~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/dansosei/p3.png~
#br
3)応力とひずみの関係(上縁から1cm)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/dansosei/p4.png~
#br
#br
メモ~
-部材の構成則(テーブル:eq plastic strain)は、塑性ひずみを入力する(Marcでは)~
-バイリニアではテーブルを設定せず、降伏応力を入力するだけでよい~
-トリリニアでは、真ひずみから弾性ひずみを引いた値を塑性ひずみとして設定する。~
弾性ひずみは、その降伏点時の応力を用いて計算する。~
-ファイバー要素で分割した応力は、job->解析結果ですべて表示で、結果に表示できる。~
#br
#br


http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tunoda/scr20.png 

**温度応力 [#sb33b144]

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tunoda/scr21.png
ケーブルプレストレスを導入できるかの検討~


・初期条件▸新規(構造)▸仮想温度~
→×(この設定だとうまくいかない)~
                   
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/ka1.png~
#br
・ケーブル要素は用いない(高度な計算になり複雑になるため)~
・代わりに梁要素に温度荷重を与えてケーブルのプレストレスとする~
・ソリッドモデル(斜張橋)で構造解析するのは、計算容量が大きくなるため骨組みモデルで計算する~
#br
・結果(アニメーション)~
300kN~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tunoda/dep2.gif
#br
変形図(鉛直変位Y)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/sa4.png~
#br
以下の方法で温度を導入できる~


|境界条件▸新規(状態変数)▸節点温度(-200℃)▸節点選択|~
|or~
|境界条件▸新規(状態変数)▸要素温度(-200℃)▸要素選択|~
#br
#br
・簡易モデルで温度応力が導入できているかの確認~
片持ち梁をメッシュ分割(2)し、3つの節点 A,B,C 点それぞれ~
温度を与えて変位を確認する~

試しに線膨張係数を小さく~
1.2×10-5 →1.2×10-7~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tunoda/dep3.gif 



(1)A=-100, B=0, C=0~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/s10.png~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/s11.png~
#br
(2)A=0, B=-100, C=0~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/s8.png~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/s9.png~
#br
(3)A=0, B=0, C=-100~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/s6.png~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/s7.png~
#br
(4)A=-100, B=-100, C=-100~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/s12.png~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/s13.png~
#br
(5)A=0, B=-100, C=0 、左側の要素に熱膨張係数を入れない~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/s14.png~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/s15.png~

***6月~ [#e466ef41]

・節点温度を与えると、節点に接している要素全体に温度応力がかかる~
・熱膨張係数を与えている部材のみ、温度応力がかかる~
(熱膨張係数が与えていなければ節点が共有されている要素には温度応力はかからない)~
#br
, 日付 , 時間帯 , 作業時間 , 内容 , 立会 
,6/2,10:00~18:00, 6h ,ゼミ 弾塑性解析,後藤さん 青木さん 石黒さん
,6/3,10:00~18:00, 6h , 弾塑性解析,
,6/4,10:00~18:00, 6h , 弾塑性解析,
,6/7,10:00~18:00, 6h ,ゼミ 弾塑性解析,後藤さん 青木さん 石黒さん
,6/8,10:00~18:00, 6h , 弾塑性解析,
,6/9,10:00~18:00, 6h ,ゼミ 弾塑性解析,後藤さん 青木さん 石黒さん
,6/10,10:00~18:00, 6h , 弾塑性解析,
,6/11,10:00~18:00, 6h , 弾塑性解析,
,6/14,10:00~18:00, 6h ,ゼミ 弾塑性解析,後藤さん 青木さん 
,6/15,10:00~18:00, 6h , 弾塑性解析,
,6/16,10:00~18:00, 6h ,ゼミ 弾塑性解析,後藤さん 青木さん 石黒さん
,6/17,10:00~18:00, 6h , モデル作成,
,6/21,11:00~18:00, 5h ,ゼミ モデル作成,後藤さん 青木さん 石黒さん
,6/22,11:00~18:00, 5h , モデル作成,
節点温度~
・A、B、C全てに‐100℃をかけると、A,B,Cそれぞれ2倍縮んでいる。~
・片側の節点だけだと、半分の値だけしか縮まらない~
・ある要素に温度応力をかけるときは、両側の節点に節点温度を与える必要がある。~
→完了(2節点での温度勾配を考慮するのではないか)~
#br
#br
#br



*Salome_meca [#cfa5c1ef]

**総括 [#c0ee7060]
まず,解析ソフトSalome_mecaを用いて研究するにあたり,以下の最低限達成するべき目標を設定した.
1. salome_mecaでケーブル(プレストレス)を導入できるか~
2. ケーブルと主桁をピン結合とし,解析できるか~
3. 弾性解析,弾塑性解析は行えるか~
4. 斜張橋の全体モデルでも同様の解析が可能であるか~

 まずわかったことは,salome_mecaで「ケーブル要素」があるが,この要素は本来コンクリートのプレストレスのようなモデル
 として用いられることを想定しており,我々がモデル化したい考えていた斜ケーブルのような挙動(引張のみ)
 を再現するのには適さない.

 またケーブルプレストレスに関しては,温度応力(厳密には温度変化?)で,緊張力を再現できる.(例:20℃→200℃)
 これは簡易モデル(片持ち梁をケーブルで吊った構造)において,線形解析を行い,中村先生が行った解析結果と同様の
 たわみ値が得られた.
#br
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2105/scr15.png~
熱応力・梁要素(ケーブル)片持ち梁【10m・集中荷重】たわみ-4.15603mm~
#br
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2105/scr16.png~

・たわみ0.03373mm~
・鉛直荷重1kN~
・ケーブルの温度変化-200℃~
・ソリッドの断面10cm×10cmの正方形、ケーブル半径10cmの円形~
#br
#br

 続いて弾塑性解析については,ソリッド要素のみのモデルでは,弾塑性解析を行うことは可能であった.(適切な挙動であるかは別として)
 ただ梁要素を含むモデルとなると,急激に難易度が上がった.おそらくsalomeの設定で梁要素の弾塑性解析を行う際はmultfiber?を設定する
 必要があるらしいが,
・(弾塑性)梁モデル→multfiber?(現在×)~
・(弾塑性)ソリッドモデル→〇~
・(弾塑性)ソリッドと梁のモデル→× salome_mecaでは弾塑性解析をするとき、すべての部材(異なる材料があるとき)に弾塑性の設定を与える必要がある(仮説)

斜張橋モデル(ソリッドと梁) 線形解析×~

ソリッドモデルの弾塑性解析→結果が中村先生と異なる。ケーブルに曲げが生じている?(構造としてはありえない)ヒンジにできていない?~




http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tunoda/depo2.gif~
scalefactor20倍
#br
#br
&size(18){6/28(火)~}; &br;~

弾塑性解析 350kN~
・DEPL_自由端~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tunoda/dep350.png    
,荷重(kN),変位(mm)
,0,0
,35,1.446695
,70,-0.930408
,105,0.564237
,140,2.28589
,175,4.37557
,210,6.64201
,245,9.06512
,280,11.7471
,315,15.1006
,350,19.2662
#br

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2106/scr31.png~
斜張橋モデル
#br
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tunoda/depl350.gif




#br
&size(18){6/15(火)~}; &br;~
#br
#br

&size(16){【斜張橋モデルの作成】};~


**斜張橋モデル [#k3423ef0]
#br
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2106/scr31.png~
#br
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2106/scr32.png~

PCがかなり重くなる...
#br
#br

モデルの作成方法(Geometry)は、~
①まずケーブル(Line)14×4=56本を作り、それらをfuse(fuse1とする)~
②主塔(Line)をつくり、同じくfuse(fuse2とする)~
③床版(solid)をつくる(Box1とする)~
④ケーブルの縁端部(solid側)の節点をグループ化する(54個)、載荷・固定のグループ化~
⑤fuse1、fuse2、Box1、それぞれメッシュを切ってcompound~
が一連の流れになる~
 
#br
・しかし現在、作成したモデルで線形解析ができない。~
 impossibility, the node N5 carries the degree of freedom of rotation  DRX~
エラー内容は「ノード(N5?)の回転自由度がない」といっているらしいので~
fuseをせずにモデルを作成できるか模索してみる。~
fuseをすると部材同士が一体化されるので、実際の斜張橋を構造的に再現できていない可能性がある。~
実際の斜張橋は、ケーブルと塔との結合部はヒンジになっており回転は許されるようになっている。~
#br





#br
#br
**弾塑性解析 [#s994a749]

&size(18){6/8(火)~}; &br;~
【弾塑性解析】~
MULTIFIBREがうまくいかないので、とりあえずニュートン法でやってみることにした~
「2部材あったとき両方に構成則を与える必要がある」という制約があるらしいので~
いままで梁要素(ケーブル)をソリッド(円柱)に変えてモデルを作成した。~

・最初に上記のモデルで温度応力を与えた状態のみで解析可能であるか検証~

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2106/scr18.png~

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2106/scr19.png~

#br
続いて同じモデルに対して弾塑性解析(ニュートン法)できるかの検証

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2106/scr17.png~

#br

#br
接線同士を重ね合わせてヒンジ?とした。~




http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tunoda/depo2.gif


http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tunoda/depo2.png
#br
&size(18){6/2(火)~}; &br;~

現時点でわかっていること~
#br

・POU_D_E、POU_D_T要素は非線形では使えない(単純なオイラー・ティモシェンコ梁の有限要素では~
単純に増分塑性を計算することはできない。断面積分を適用する必要がある)~
・代わりに使える要素→POU_D_EM、POU_D_GM(マルチファーバー)~
・要素特性を定義するコマンドAFFE_CARA_ELEMの設定が必要~
AFFE_CARA_ELEM~
・GEOM_FIBRE~
・MULTIFIBRE~
の設定がうまくいかない。~


   !  DEFI_GEOM_FIBRE VALE : Pour .NOMTM          $$XNOM il y a 1 valeurs, ce devrait !
   ! être un multiple de 3    

・ソリッドと梁のモデル(ソリッドだけ弾塑性)は~

   ! Aucune maille du maillage mesh n'a été affectée par des éléments finis.!

   !  In the mesh" mesh" the mesh" 3531" is of type" TETRA4" (neither TRIA3 nor !
   ! QUAD4)  
#br
#br


弾塑性解析 350kN~
・DEPL_自由端~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tunoda/dep350.png    
,荷重(kN),変位(mm)
,0,0
,35,1.446695
,70,-0.930408
,105,0.564237
,140,2.28589
,175,4.37557
,210,6.64201
,245,9.06512
,280,11.7471
,315,15.1006
,350,19.2662
#br

***5月~ [#c467ef20]
, 日付 , 時間帯 , 作業時間 , 内容 , 立会 
,5/2,11:00~14:30, 3h ,ケーブル解析, 
,5/6,10:00~18:00, 6h ,ケーブル解析,
,5/7,10:00~18:00, 3h ,ケーブル解析,
,5/10,10:00~18:00, 6h ,ゼミ ケーブル解析,後藤さん 青木さん 石黒さん
,5/11,10:00~18:00, 6h ,ケーブル解析,
,5/12 , 10:00~18:00, 6h ,ゼミ ケーブル解析,後藤さん 青木さん 石黒さん
,5/13,10:00~18:00, 6h ,ケーブル解析,
,5/17 , 10:00~18:00, 6h ,ゼミ ケーブル解析,後藤さん 青木さん 石黒さん
,5/18 , 10:00~18:00, 6h ,ケーブル解析,
,5/19 , 10:00~18:00, 6h ,ゼミ ケーブル解析,後藤さん 青木さん 石黒さん
,5/20 , 11:00~18:00, 5h ,ケーブル解析,
,5/21 , 11:00~18:00, 5h ,ケーブル解析,
,5/24 , 10:00~18:00, 6h ,ゼミ 弾塑性解析,後藤さん 青木さん 
,5/27 , 10:00~18:00, 6h ,弾塑性解析,
,5/28 , 10:00~18:00, 6h ,弾塑性解析,
,5/31 , 10:00~18:00, 6h ,ゼミ 弾塑性解析,後藤さん 青木さん 
#br
#br
&size(18){5/24(月)~}; &br;
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2105/scr16.png~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tunoda/depl350.gif

ソリッドと梁要素(ケーブル)のモデルで、温度応力を与えた解析はまわることが確認できた。~
・たわみ0.03373mm~
・鉛直荷重1kN~
・ケーブルの温度変化-200℃~
・ソリッドの断面10cm×10cmの正方形、ケーブル半径10cmの円形~
・ケーブルの断面形状を円柱→四角形~

梁モデル(片持ち梁をケーブルで吊った構造)の弾塑性解析は色々試して入るがうまくいっていない~

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tunoda/scr20.png 


http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tunoda/scr21.png


#br
#br
&size(18){5/17(木)~}; &br;
・結果(アニメーション)~
300kN~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tunoda/dep2.gif

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2105/scr14.png~
熱応力・梁要素(ケーブル)片持ち梁【1m・分布荷重】たわみ-3.4429mm~

#br
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2105/scr15.png~
熱応力・梁要素(ケーブル)片持ち梁【10m・集中荷重】たわみ-4.15603mm~

試しに線膨張係数を小さく~
1.2×10-5 →1.2×10-7~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tunoda/dep3.gif 



#br
&size(18){5/13(木)~}; &br;
#br


**温度応力 [#f0292266]
***ケーブルPrの導入 [#l82aa05a]
#br
#br
・ケーブル要素を用いたモデルで今まで解析を試みていたが、一旦~
モデルを梁要素で作成し、熱応力でケーブルプレストレスを再現することにした。~

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2105/scr8.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2105/scr9.png~
#br
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2105/scr10.png~
#br
スケール10倍~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2105/scr11.png~
・集中荷重1000N~
・温度−200℃~
・先端のたわみ量 -4.15603mm
#br
#br
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2105/scr12.png~
・ケーブルのはりのノードが3つになってしまっている~
→メッシュの設定の Numbwr of segments でセグメント数に1を入力することで解決~
・先端のたわみ量 0.03mm(ケーブルプレストレスなし)
#br
#br
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2105/scr13.png~
・節点がヒンジになっていることを確認

#br
#br
&size(18){5/6(木)~}; &br;
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2105/scr16.png~
#br
ソリッドと梁要素(ケーブル)のモデルで、温度応力を与えた解析はまわることが確認できた。~
・たわみ0.03373mm~
・鉛直荷重1kN~
・ケーブルの温度変化-200℃~
・ソリッドの断面10cm×10cmの正方形、ケーブル半径10cmの円形~

・直近の課題はケーブル要素を用いた簡易的なモデルを作成して青木さんのFEM解析の結果と比較する~
,片持ち梁(分布),片持ち梁(先端集中),片持ち梁(分布・ケーブルあり),片持ち梁(先端荷重・ケーブルあり)
, ○ , ○ , × , × 
・固定面に載荷面が触れていても問題なく実行できた~
・しかし、ケーブルの設置点と載荷部分が重なるとエラーが生じた。~
   !  les 3116 mailles imprimées ci-dessus n'appartiennent pas au modèle       !
   !  et pourtant elles ont été affectées dans le mot-clé facteur : FORCE_FACE 
・対策として、~
(1)載荷面をケーブルの接地点から1mmほどずらして設定し直す~
(2)ケーブルはケーブル要素、四面体のボックスは3Dで分けて設定していたのを、ケーブル・ボックス含めて全体を3Dに設定且つ~
ケーブルはケーブル要素に設定するに変更(ケーブルとボックスの接地点を解析に含まれるようにするため)~
・次のエラーメッセージ
   ! Erreur utilisateur :                                            !
   !   Vous voulez contraindre le ddl DRY sur un ensemble de noeuds, !
   !   Mais ce ddl n'existe sur aucun de ces noeuds.                 !
(ノードが含まれていない?)~
↓~
・メッシュの切り方を3Dだけ→3D且つ、アルゴリズム(補完)2D、1Dまで設定しメッシュを切る~
・ケーブルの固定点のBC&Loadの設定を変えてみる~
・梁だけのモデル(ケーブル要素含む)、ケーブル単体のモデルを作成~
→いずれも同じエラーメッセージだった~

・ケーブル要素は非線形でしか使えない可能性がある(線形では使えない)~
・非線形モデルを作成して試してみる~
・梁とBARRE要素を用いてモデルを作成してみる~

→BARRE要素と梁を用いたモデルで解析を回してみたが同じエラーが出た~
   ! Erreur utilisateur :                                            !
   !   Vous voulez contraindre le ddl DRY sur un ensemble de noeuds, !
   !   Mais ce ddl n'existe sur aucun de ces noeuds.                 !
梁モデル(片持ち梁をケーブルで吊った構造)の弾塑性解析は色々試して入るがうまくいっていない~
#br
#br

&size(18){5/2(月)}; &br;

・(100,100,1000)で集中荷重も試みた~
→載荷面を上面に指定~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2105/scr5.png~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2105/scr14.png~
熱応力・梁要素(ケーブル)片持ち梁【1m・分布荷重】たわみ-3.4429mm~

・問題なく解析できた~
・やはり(100,100,10000)という寸法で解析すること(片持ち梁)に無理があったのか?..
#br
#br
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2105/scr15.png~
熱応力・梁要素(ケーブル)片持ち梁【10m・集中荷重】たわみ-4.15603mm~

***4月~ [#j3a34377]


#br
, 日付 , 時間帯 , 作業時間 , 内容 , 立会 
,4/1,10:00~18:00, 6h ,3年生課題,        
,4/5~9,10:00~17:00, 6h ,3年生課題,
,4/21,10:00~18:00, 6h ,salomeゼミ ケーブル要素解析,後藤さん 青木さん 石黒さん
,4/22,10:00~17:00, 6h ,ケーブル要素解析,
,4/26,10:00~18:00, 6h ,ゼミ ケーブル要素解析,後藤さん 青木さん
,4/27,10:00~18:00, 6h ,ケーブル要素解析,
,4/28,10:00~17:00, 5h ,ケーブル要素解析,
,4/30,11:00~18:00, 4h ,ケーブル要素解析,
#br
#br
&size(18){4/30(金)}; &br;

モデルが(100,100,10000)だったので(100,100,1000)にしてみて試してみた。~
↑10cm,10cmの断面に対して、10mが長すぎるのではないかと予想~

・集中荷重(先端)~

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2105/scr4.png~

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2105/scr7.png~

・10mにしたらうまくいった

#br
#br
&size(18){4/28(水)}; &br;

片持ち梁のモデルの解析ができていないので優先して取り組むことにした。~
モデルの荷重のかけ方は等分布荷重だったので、載荷面をボックスの上面に指定して解析を回していた。~
→これを線載荷に変更して同様に試みた。~

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2105/scr3.png
   ! <S> Exception user raised but not interceptee.                !
   ! The bases are fermees.                                        !
   ! Type of the exception: error                                  !
   !                                                               !
   ! Solver MUMPS:                                                 !
   !  The solution of the linear system is too vague:              !
   !  Computed error:  1.58999e-05                                 !
   !  Acceptable error:  1e-06 (RESI_RELA)                         !
   !                                                               !
   !  Advices:                                                     !
   !  One can increase the value of the key word SOLVER/RESI_RELA. !
・26日のエラーとほとんど同じ

・同じモデルで集中荷重(先端)で試した~
→うまくいかなかった

#br
#br
&size(18){4/26(月)}; &br;

**簡易モデル [#nbc03537]
salomeで、ケーブル要素を用いた簡易的なモデルを作成し解析を回してみた。~
↑青木さんがやられていた解析プログラムで算出した値と比較する~
↑中村先生の解析結果と比較する~
+片持ち梁モデル~
+片持ち梁にケーブルをつないだモデル(荷重なし)~
+片持梁にケーブルをつないだモデル(荷重あり)~
#br
1.片持ち梁モデル~
・(x,y,z)=(100,100,10000)~
・鋼材E=200GPa,σ=0.3~
・荷重1km/m(等分布)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/April/scr1.png
         ! Exception user raised but not interceptee.                    !
         ! The bases are fermees.                                        !
         ! Type of the exception: error                                  !
         !                                                               !
         ! Solver MUMPS:                                                 !
         !  The solution of the linear system is too vague:              !
         !  Computed error:  5.59007e-05                                 !
         !  Acceptable error:  1e-06 (RESI_RELA)                         !
         !                                                               !
         !  Advices:                                                     !
         !  One can increase the value of the key word SOLVER/RESI_RELA. !
・収束判定基準がおかしい..?
・収束判定の設定..
#br
#br

2.片持ち梁にケーブルをつないだモデル(荷重なし)~
・ケーブル半径1cm、E=200GPa~
・ケーブル張力6.233kN~

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/April/src2.png
   ! <S> Exception user raised but not interceptee.                            !
   ! The bases are fermees.                                                    !
   ! Type of the exception: error                                              !
   !                                                                           !
   !  les 9 mailles imprimées ci-dessus n'appartiennent pas au modèle          !
   !  et pourtant elles ont été affectées dans le mot-clé facteur : FORCE_FACE !

#br
#br
#br
*メモ [#yb1a7ac2]
#br
**2021 時間割 [#kd898b2c]
***前期 [#j74bf09b]
,曜日,1 2,3 4,5 6,7 8,9 10
,月,,,,ゼミ,
,火,,都市システム計画特論,土質工学特論,,
,水,,外国語文献講読,内輪ゼミ,,
,木,,材料設計学特論,,,科学技術者倫理特論(後半)
,金,,,,地域防災学特論,

***後期 [#b3222793]
,曜日,1 2,3 4,5 6,7 8,9 10
,月,Introduction to Systems Design Engineering(後半),,,,
,火,地域エネルギー特論,,,,
,水,,,,地域産業アントレプレナー論(前半),
,木,,,,,
,金,,構造力学特論,,,
片持ち梁のモデルの解析ができていないので優先して取り組むことにした。~
モデルの荷重のかけ方は等分布荷重だったので、載荷面をボックスの上面に指定して解析を回していた。~
→これを線載荷に変更して同様に試みた。~

**Marc_Mentat [#c5704930]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2105/scr3.png
   ! <S> Exception user raised but not interceptee.                !
   ! The bases are fermees.                                        !
   ! Type of the exception: error                                  !
   !                                                               !
   ! Solver MUMPS:                                                 !
   !  The solution of the linear system is too vague:              !
   !  Computed error:  1.58999e-05                                 !
   !  Acceptable error:  1e-06 (RESI_RELA)                         !
   !                                                               !
   !  Advices:                                                     !
   !  One can increase the value of the key word SOLVER/RESI_RELA. !
・26日のエラーとほとんど同じ

Q1.ポイントやソリッドが表示されない~
・(0 0 0)→0 0 0 にしてみる~
・ビューのモデルエンティティタイプのポイントに☑を入れる~
・ビューに合わせるをクリックする~
・同じモデルで集中荷重(先端)で試した~
→うまくいかなかった

・Marc_Mentatの初期条件→Streamline Pyrolysisモデルを設定するステップの一つ変数の初期化で行うためのタブ。~
Marcの熱分解と表面後退の機能は、Advanced Thermal Analysis Softwareの頭文字をとってATASと呼ばれるプロジェクトの一環として開発された。~
元々は航空宇宙アプリケーションのために開発されたが、あらゆる熱分解や表面後退のアプリケーションに使用することができる。~
固体推進剤ロケットモーターは、非常に高い熱流束にさらされる。このような条件を維持するために、複雑な材料が設計される。ATASは、そのような材料が受ける2つの重要な現象をモデル化することができる。~
熱分解と呼ばれる、温度による材料の熱劣化 ~
温度、化学的または機械的な影響による表面の消失(これを表面後退と呼ぶ ~
熱分解プロセスは、温度による固体物質の分解。この劣化または分解により、熱分解ガスと呼ばれるガスが発生する。熱分解性の材料は、複数の理由から熱保護材として使用されることが多い。~熱分解プロセスは吸熱性である。~
Marcにはいくつかの熱分解モデルがあります。
そのうちの一つがtreamline Pyrolysis モデル。このモデルは、固体材料の分解を記述する法則と、熱分解ガスを考慮に入れて構成されている。




#br

**Marc_Mentat資料 [#qa4749de]
Mark_Mentat▸ヘルプ▸ユーザーズガイド▸Basics of Mentat~

・日本語~
・ファイル形式:Word~
・自分で訳したため多少不備があるかもしれない~


-[[Marc Mentat 2018  基礎知識:http://www.str.ce.akita-u.ac.jp/~gotouhan/public_html/j2021/tsunoda/2108/Marc1.docx]]
-[[How This Guide Communicates With You :http://www.str.ce.akita-u.ac.jp/~gotouhan/public_html/j2021/tsunoda/2108/Marc2.docx]]
-[[Background Information on Performing a Finite Element Analysis :http://www.str.ce.akita-u.ac.jp/~gotouhan/public_html/j2021/tsunoda/2108/Marc3.docx]]
#br
モデルが(100,100,10000)だったので(100,100,1000)にしてみて試してみた。~
↑10cm,10cmの断面に対して、10mが長すぎるのではないかと予想~

・集中荷重(先端)~

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2105/scr4.png~

**salome 資料 [#x6f238d3]
-[[SALOME-Mecaの使用法解説:http://opencae.gifu-nct.ac.jp/pukiwiki/index.php?SALOME-Meca%A4%CE%BB%C8%CD%D1%CB%A1%B2%F2%C0%E2]]
#br
***salomeキーワード [#a1420495]
&size(16){Fuse}; ~
結果:単一のソリッド。~
目的:単純な形状を組み合わせることにより、複雑なモデルを構築。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2105/scr7.png~

&size(16){Partition}; ~
結果:面を共有する2つの接続されたソリッド。~
目的:モデルの個別の領域(流体/固体、コンクリート/鋼...)を確保するのに役立つ。~
・10mにしたらうまくいった

&size(16){Compound}; ~
結果:2つの別々のソリッドを含むオブジェクトを作成。(2つのソリッドは接続されていない、単なる形状のセットでありそれ以上のもではない)~
目的:図形のコレクションに操作を適用できるようにする。~


#br
#br

**資料メモ [#je54116d]
【ファイバー要素(モデル)とは】~
ファイバー要素の特徴として、主に以下の3つが挙げられる。~
・部材断面は複数の微小断面(以降、セグメントと呼ぶ)に分割されている。~
・部材断面は平面保持仮定が成立していると仮定されている~
・各セグメントは構成側(応力度-ひずみ関係)を持ち、セグメント間に相互作用は存在しない。~
・(100,100,1000)で集中荷重も試みた~
→載荷面を上面に指定~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2105/scr5.png~

部材全体の応力-ひずみ関係式について
・問題なく解析できた~
・やはり(100,100,10000)という寸法で解析すること(片持ち梁)に無理があったのか?..



#br
#br

・直近の課題はケーブル要素を用いた簡易的なモデルを作成して青木さんのFEM解析の結果と比較する~
,片持ち梁(分布),片持ち梁(先端集中),片持ち梁(分布・ケーブルあり),片持ち梁(先端荷重・ケーブルあり)
, ○ , ○ , × , × 
・固定面に載荷面が触れていても問題なく実行できた~
・しかし、ケーブルの設置点と載荷部分が重なるとエラーが生じた。~
   !  les 3116 mailles imprimées ci-dessus n'appartiennent pas au modèle       !
   !  et pourtant elles ont été affectées dans le mot-clé facteur : FORCE_FACE 
・対策として、~
(1)載荷面をケーブルの接地点から1mmほどずらして設定し直す~
(2)ケーブルはケーブル要素、四面体のボックスは3Dで分けて設定していたのを、ケーブル・ボックス含めて全体を3Dに設定且つ~
ケーブルはケーブル要素に設定するに変更(ケーブルとボックスの接地点を解析に含まれるようにするため)~
・次のエラーメッセージ
   ! Erreur utilisateur :                                            !
   !   Vous voulez contraindre le ddl DRY sur un ensemble de noeuds, !
   !   Mais ce ddl n'existe sur aucun de ces noeuds.                 !
(ノードが含まれていない?)~
↓~
・メッシュの切り方を3Dだけ→3D且つ、アルゴリズム(補完)2D、1Dまで設定しメッシュを切る~
・ケーブルの固定点のBC&Loadの設定を変えてみる~
・梁だけのモデル(ケーブル要素含む)、ケーブル単体のモデルを作成~
→いずれも同じエラーメッセージだった~

・ケーブル要素は非線形でしか使えない可能性がある(線形では使えない)~
・非線形モデルを作成して試してみる~
・梁とBARRE要素を用いてモデルを作成してみる~

【西海岸の風に吹かれて】~
白上謙一「生物学と方法」~
個人の活動を、その個人が所属する社会や共同体がどのように評価するのか、その際の評価尺度として何が適切であるのか構造的問題は、学術研究以外の場でも共通している
→BARRE要素と梁を用いたモデルで解析を回してみたが同じエラーが出た~
   ! Erreur utilisateur :                                            !
   !   Vous voulez contraindre le ddl DRY sur un ensemble de noeuds, !
   !   Mais ce ddl n'existe sur aucun de ces noeuds.                 !
#br
#br

重要な発見をしようとするよりも、自分の発見を重要なものにすることに努力するべきである~
質的な結果を目指してのみ、量的な研究を遂行せねばならぬ~

まずは自分の中の評価軸を確固たるものにする。そしてその評価軸が社会(自分が所属している組織やコミュニティ)
の評価軸とずれているときに、どのように調整するか。どこまでなら妥協できてどこからは許されないのか、その境目は
を見極める力を磨いていく必要がある。~

**salome例題ファイル [#uf69aefa]


【斜張橋の基本計画設計法】(長井正嗣・井澤衛・中井宏 著 1997)~
斜張橋の主要な構成要素は、桁、ケーブルと塔で即径間の最上段ケーブルをアンカーケーブルと呼ぶ。橋端部には不反力用のペンデル支承が設けられる。
死荷重が作用した状態で、ケーブルの桁支持点(定着点)が剛支点になるようにケーブル張力を選んでおけば、支間をLGとした連続橋の曲げモーメントを得ることができる
曲げモーメントはケーブルの吊間隔LGの2乗に比例して発生するため、吊間隔を小さく選んでおくと、支間に関係なく死荷重時の曲げモーメントを小さくできる。

桁橋に比べて、斜張橋の桁高さはかなり小さくできる

ケーブルを斜めには配置していることから、桁はケーブル張力を受け持つことが不可欠となる。そのため桁には圧縮軸力が作用する。圧縮軸力は、支間が長くなるにつれて大きくなる
圧縮部材の設計では、座屈と呼ばれる不安定現象に注意する必要がある。したがって斜張橋の長大化を計る上では、圧縮軸力の増大に対して桁の耐荷力が確保できる断面を選ぶことが
耐風安定性の確保とともに重要な課題となっている

側径間のアンカーケーブルに力が流れ、その伸び剛性でもって塔頂方向変位が拘束される。すなわちアンカーケーブルが橋全体の鉛直面内剛性に大きく寄与している。また橋端部に
おいて負反力が生じるのも斜張橋の特徴である

斜張橋は高次の不静定構造物で、その力学特性として、桁の曲げ剛性とケーブルの伸び剛性の比に応じて、それぞれに作用する曲げモーメントと張力が変化する。ある剛比を設定したとする
その状態から桁の曲げ剛性を大きくしていくと、桁の荷重分担率が大きくなり、曲げモーメントが増加する。またその逆もいえる。
桁高を変化させると、断面二次モーメントは大きく変化するが、薄板の腹板の背が高くなるだけで、桁の鋼重量はあまり変化しないようと考えられる。したがって桁重量を支持するように設計される
ケーブルの伸び剛性の変化は小さいままで、桁の曲げ剛性が大きく変化する。桁の曲げ剛性が大きくなるほど曲げモーメントも大きくなることを説明したが、曲げ剛性を大きいため、曲げ応力度はあまり変化しないと考えられる。

ケーブルシステムによって桁に発生する軸力の大きさや圧縮、引張の符号に差が生じる。斜張橋に採用されているケーブルシステムはほとんどが自定式と呼ばれるもので、主桁の軸力は圧縮となる。この他に、完定式、部定式、一部部定式のシステムが提案されている。完定式:橋端で桁を橋軸方向に固定し、塔位置で稼働することによって主桁の軸力は引張となる。一部他定式:主桁の一部が引っ張り、一部が圧縮の軸力となる。自定式を用いた斜張橋の長大化には限界があり、さらに一層の長大化を計る上では、一部他定式が有力な案として考えられている。

鋼斜張橋のケーブルの種類には、LCR(ロックドコイルロープ)、PWS、HiAmアンカーケーブル、NewPWSケーブルがある。この他にイギリスではスパイラルロープが使われている。
LCRは外層にZ型の鋼線(素線)を用いロープで防食対策には鋼線を亜鉛メッキする方法と外表面の塗装を併用する方法が用いられている。
LCRは海外での実績は多いが、日本での使用は少なくなっている。日本ではHiAmアンカーケーブルやNewPWSケーブルといった平行線ケーブルが用いられる。
LCRは鋼線によりを入れているために、一般にヤング係数や強度が鋼線そのものに比べて低く、また構造伸びが生じるといった問題がある。構造伸びを除去するために、使用前に一度大きな張力を導入することになる。

斜張橋の構造用部材には、ケーブル(鋼線)を別とすれば、一般橋梁と同様に鋼とコンクリート(プレストレスコンクリート)が使用される。それぞれを鋼斜張橋、PC斜張橋と呼ぶ。この他に、中央径間の桁に鋼、側径間の桁にコンクリートを用いる場合や、桁に鋼、塔にコンクリートを使用している場合があり、これらを複合橋と呼ぶ。架橋される現地の地形や、地盤条件の制約を受けて、側径間が短く計画される場合、桁に同一材料を用いると中央径間と側径間で重量バランスが崩れる。また、橋端部(橋軸方向)の支点に大きな反力が生じる。その対策として側径間にコンクリート桁が用いられる。しかし、コンクリート桁の重量は鋼桁に比べてかなり大きいため、逆に側径間部の重量が大きくなって、再び重量バランスが崩れる。そのためこのタイプの斜張橋では、側径間に複数の橋脚が設けられている。
海外の例であるが、鋼桁(I桁)とコンクリート床版をずれ止めを用いて合成させた、合成断面の桁を用いる斜張橋が提案され、支間が300m~600mの長大橋に適用されている。












**熱応力・梁要素(ケーブル) [#q5e180d7]
***片持ち梁【集中荷重・1m・10m】 [#yfc4d78e]
熱応力・梁要素(ケーブル)~
片持ち梁【集中荷重・1m・10m】~
1. 【集中荷重・1m】(hari4)~
・Geometryの作成:新しいエンティティ->基本オブジェクト->点->(0,0,0),(0,0,1000)作成~
・基本オブジェクト->線->Vertex1、2->Line1作成~
・点->(0,-1000,0)->Vertex2、3->Line2作成~
・Fuse->Line1、Line2->Partition->Fuse1、Line1、Line2~
・Line1グループ作成->node(0,0,0)-kotei->node(0,0,1000)-saika~
・Line2グループ作成->node(0,-1000,0)-kotei2~
・Mesh1:Line1->ジオメトリのグループ作成->ジオメトリ-Line1->node-kotei-saika~
・Mesh2:Line2->Wire Discretisaition->Number of Segments->セグメント数1->ジオメトリのグループ作成->ジオメトリ-Line2->node-kotei2~
・Compound_Mesh->Mesh1、2~
・Astestudy:メッシュ->Read a mesh->format-Med~
・Material:Model Definition->Assign finite element->Finite element->Everywhere->Mechanic->POU_D_E~
・Model Definition->AFFE_CARA_ELEM->POUTRE->RECTANGLE->Group-Line1->CARA-HY-HZ->Value-100-100~
・Material:Define a material->Linear isotropic elasticity->ヤング率200000->ポアソン比0.3~
・Material:Define a material->Linear isotropic elasticity->ヤング率200000->ポアソン比0.3->Thermal expansion 1.17e-0.5~
・BC and Load:EnforceDOF->Group node-kotei->Laison-ENCASTRE->DX~DRZ-o~
・EnforceDOF->Group node-kotei2->Laison-ENCASTRE->DX~DRZ(DRX以外)-o~
・FORCE_NODLE->Group node-saika->FY-1000~
・Post Processing->CREA_CHAMP->Model-AFEE_MODELE->TYPE_CHAM-NOEU_TEMP_R->OPERATION-Assingment->~
Group element-Line2->Value- -200->NOM_CMP-TEMP ~
・Material:Assigm a material->Model-model(AFEE_MODELE)->Marerial assigment->Group-Line1->mater->~
Marerial assigment->Group-Line2->mater0->External state-Temperature->Reference value-20->Group-Line2->Field-unnamed5->~
・Analysis->Material field->Structural element->Model->Load=Load->Solver->Method-MUMPS->RESI_RELA-0.01~
・Post processing CALC_CHAMP->SIGM_NOEU->MODEl-model->Material field->fieldmat->Structural element~
・Output->Set output result->Format-Med-> Result-DEPL-SIGM_NOEU~

[[1m.comm:http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/commandfile/1msyutyu.comm]]~
(commファイル:一度保存してAsterstadyでエクスポートすることによって利用可能。hdfファイルではないため、Geometry等は反映されない)

#br
2. 【集中荷重・10m】(10m)~
1.集中荷重1mとほとんど同じ。異なる点は、Geometryで10mのモデルを作ることとAstestudyの~
SolverでNPREC->-1、RESI_RELA->20にする(収束基準・条件を緩和する)~

[[10m.comm:http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/commandfile/10msyutyu.comm]]~

#br
***片持ち梁【分布荷重・1m】 [#fd101e01]
#br 
片持ち梁【分布荷重・1m】~
3. 【分布荷重・1m】(bump2)

[[bunpu.comm:http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/commandfile/bumpu1.comm]]~

1.集中荷重1mとほとんど同じ。異なる点は、BC and Loadで載荷のオプションを変更することのみ。
 FORCE_POUTRE=_F(FY=1.0,
 GROUP_MA=('saika', )),
 MODELE=model)


**ケーブル要素 [#e9968ca7]
❌線形モデル

複雑な計算になるため、一般的にbeam要素を用いる。



#br
#br
#br

**参考文献 [#g936ac10]

&size(17){「ケーブル腐食を考慮した斜張橋の終局強度および疲労寿命」 }; ~
[[pdf:http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/re1.pdf]]~
構造工学論文集Vol.  67A (2021 年3 月)  ~
中村俊一,青木由香利~
Ph.D.,東海大学,工学部土木工学科~
Ph.D.,東海大学,国際教育センター~
概要~
斜張橋におけるケーブルは命綱とも言える重要な構造要素である。そのケーブルに着目し、腐食ケーブルが橋全体の終局強度に及ぼす影響と疲労寿命に関する検討を行った。~
目的~
・想定した斜張橋モデルに常時設計荷重を作用させ,部材の応力照査を行い,その妥当性を確認する~
・3 つの斜張橋モデルが終局に至るまで荷重を漸増載荷し,終局強度および崩壊過程を明らかにする~
・想定した3 つの斜張橋モデルに重量車両を繰り返し載荷し,ケーブル疲労寿命を求め,異なるスパン長およびケーブル位置への影響を検討する~
*メモ [#yb1a7ac2]

#br
#br


&size(17){「斜吊りアーチ橋の構造特性」 }; ~
[[pdf:http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/re2.pdf]]~
構造工学論文集Vol.54A (2008 年3 月) ~
田中寛泰,中村俊一,加藤一寿~
川田工業(株),  橋梁事業部~
Ph.D.,  東海大学教授~
東海大学大学院工学研究科土木工学専攻~
概要~
主桁を斜ケーブルとCFTアーチリブの両者で吊った斜吊りアーチ橋という新形式橋梁を提案する.~
性能照査型設計により斜吊りアーチ橋の部材破壊に対する安全性を照査し,同規模の斜張橋と鋼重量を比較し考察することで~
提案橋梁形式が構造的かつ経済的に合理的であることを明らかにした .~
目的~
・コンクリート充填鋼管(CFT)を用いた斜吊りアーチ橋の提案~
・限変位解析により設計断面力を求め,そして限界状態設計法を用いて終局限界状態で耐力照査~
・アーチリブの座屈強度を弾性有限変位解析により検討~
・経済性を検討するため,同規模の斜張橋と鋼重量を比較し考察~


Reference List
#br
#br

-ケーブル腐食を考慮した斜張橋の終局強度および疲労寿命
[[pdf:http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/re1.pdf]]~
-斜吊りアーチ橋の構造特性~
-腐食した橋梁用亜鉛メッキ鋼線の強度特性 ~
-橋梁用ケーブルの最近の話題と展望~
-斜吊りアーチ橋の静的構造特性と崩壊挙動に関する考察~
-Dead Load Analysis of Cable-Stayed Bridge~
-3径間連続マルチケーブル斜張橋の部材断面力算定法の提案~

&size(17){「腐食した橋梁用亜鉛メッキ鋼線の強度特性」 }; ~
[[pdf:http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/re3.pdf]]~
土木学会論文集No.731/1-63,367-377,2003.4~
鈴村 恵太・中村 俊一・樽井 敏三~
日本製鐵(株)相模原技術開発部 マネージャー~
Ph.D.東海大学教授~
新日本製鐵(株)鉄鋼研究所鋼材第二研究部~
概要~
吊り橋や斜張橋ケーブルを構成する亜鉛メッキ鋼線を対象とし、湿潤巻き腐食進行法により異なる腐食レベルの~
試験片を作成し、腐植度と強度の関係を研究した。腐食により表面に凹凸が生じることを確認した後, 表面凹凸を~
機械切削した鋼線の伸びが回復することを見出し, 強度低下の主原因は腐食による表面凹凸であることを明らかにした。~
目的~
・腐食により生じた試験片の表面凹凸を機械切削し, 表面を平坦に戻した後, 引張試験を実施する~
・これにより, 表面凹凸が亜鉛めっき鋼線の強度低下にどの程度影響しているかについても研究する~




#br
#br



&size(17){「大規模解析による鋼斜張橋の車両重量および温度変化に対する挙動変化」 }; ~
[[pdf:http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/re4.pdf]]~
構造工学論文集Vol. 66A (2020 年3 月) ~
清水暁央、久國陽介、釘宮哲也、渡部一雄、久野勝美、篠原聖二、赤松伸祐~
(株)東芝,生産技術センター~
(株)東芝,研究開発センター~
(株)東芝,生産技術センター ~
阪神高速道路(株),計画部~
(一財)阪神高速道路技術センター,調査研究部~
概要~
既設の長大橋全体を3次元シェルおよびソリッド要素で精緻にモデル化する.さらに線形静解析を実施し,~
実測可能な応答と計算結果を比較することで,モデルの妥当性を検証する.実測可能な応答として,車両荷重による橋梁のたわみ,温度日変動による橋梁の熱伸縮を扱う~
目的~
・線形応答の妥当性を示すことで,車両荷重,温度日変動といった異なる負荷に対して同じ
モデルで応答が評価できることを確認する.~


-大規模解析による鋼斜張橋の車両重量および温度変化に対する挙動変化~
-2011年東北地方太平洋沖地震における横浜ベイブリッジの応答~
-大規模計算による長大橋地震応答解析の構造要素モデルの高度化に関する検討~
-斜張橋ケーブルの耐久性評価と今後の維持管理~
-斜張橋ケーブルの耐久性評価と今後の維持管理に関する小委員会・報告書~
-2016年熊本地震における鋼斜張橋の損傷メカニズム~
-Parametric Studies On Cable-stayed Bridges~
-鋼斜張橋ケーブルの疲労安全性と部材安全率に関する研究~
-ケーブル等の安全率の設定と評価~
-橋梁ケーブルの限界状態設計法~
-斜張橋(神納橋)の設計について~
-ファジィ理論を用いた斜張橋ケーブルの最適プレストレス力決定法~
-長大斜張橋の想定規模地震時の非線形挙動に関する研究~
-吊形式橋梁の現状と将来~
-長大斜張橋のL2地震動に対する耐震対策~
#br
#br



&size(17){「橋梁用ケーブルの最近の話題と展望」 }; ~
[[pdf:http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/re5.pdf]]~
土木学会論文集、1992.3~
三田村武,  中井博,  渡邊英一,  杉井謙一~
概要~
吊り橋メインケーブル、斜張橋ケーブルについて、その種類や設計上検討すべき項目に関する玄関と最近の動向を述べる~
目的~
・吊り橋・斜張橋の双方のケーブルにつき、損傷例、腐食ワイヤの強度特性、検査方法、および維持管理上の要点を示す。~
・さびないケーブル、より強いケーブル、より軽いケーブルの特性と可能性について論ずる~


#br
#br


&size(17){「2011年東北地方太平洋沖地震における横浜ベイブリッジの応答」 }; ~
[[pdf:http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/re6.pdf]]~
土木学会論文集A1(構造・地震工学), Vol. 69, No. 2, 372-391, 2013.~
藤野陽三・シリゴリンゴディオンシウス・並川賢治・矢部正明~
東京大学特任教授  大学院工学系研究科社会基盤学専攻 ~
同上~
首都高速道路株式会社  技術部技術推進課  課長~
株式会社長大  構造事業本部  副技師長~
概要~
主桁の振動は橋軸直角方向の応答が卓越した,橋軸方向の主桁と主塔および端部橋脚の応答記録から,主塔と端部橋脚位置のリンク支承が~
所要の免震性能を発揮した,主塔と端部橋脚のウィンド沓は橋軸直角方向に衝突していた形跡がある本震時の大振幅応答時には,主桁の~
橋軸直角方向水平たわみ1次モードや主桁の鉛直方向たわみ1次モードの固有振動数が変化しただけでなく,それぞれのモード形において,~
小振幅時には見られない連成挙動が同定されたことを明らかにしている.~
目的~
2011年東北地方太平洋沖地震時の横浜ベイブリッジの地震応答に関する解析結果について述べる.~
1)地震動の特性~
2)主桁の応答特性、特に大振幅応答時の非線形応答特性~
3)地震応答記録からのシステ~ム同定~
4)端部橋脚と主塔におけるリンク支承の性能評価~
5)地震後の現地観察結果との対比~
の項目に着目し解析を行った。~

#br
#br


&size(17){「大規模計算による長大橋地震応答解析の構造要素モデルの高度化に関する検討」 }; ~
[[pdf:http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/re7.pdf]]~
土木学会論文集 A2(応用力学),  Vol. 73,  No. 2 (応用力学論文集 Vol. 20),  I_769-I_780,  2017.~
八ツ元  仁・馬越  一也・金治  英貞・中村  真貴・野中  哲也~
阪神高速道路株式会社  大阪管理局保全部~
株式会社地震工学研究開発センター~
阪神高速道路株式会社  技術部技術推進室~
株式会社地震工学研究開発センター~
名古屋工業大学大学院教授  社会工学専攻~
概要~
本検討では,スーパーコンピュータ「京」を用いることで,長大斜張橋の非線形領域における挙動を求~
めるという大規模計算を行った.本解析では,兵庫県南部地震において大きな被災を受けた長大斜張橋~
を解析対象として計算を行い,実被害状況と計算結果の比較を通して,大規模計算自体が正確に行われて~
いること,高度化した構造要素モデルを用いた部位で被災状況を概ね計算で再現できていることを確認した.~
目的~
・大規模計算機を正しく実行できているのかの確認~
・改良を重ねてモデル化した構造部材が想定通りの挙動するのかの確認~
・高い精度で橋梁が保有する耐震性・安全性を評価する手法の確立~
・計算結果に対する妥当性の検証~



#br
#br


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