#author("2022-05-20T12:43:39+09:00","default:kouzouken","kouzouken")
#author("2022-05-20T13:13:19+09:00","default:kouzouken","kouzouken")
#contents
#br
#br
*ケーブル腐食を考慮した斜張橋における静的・動的応答解析 [#u604b131]
**当面の課題 [#pe4564e4]
#br
☐時間軸を考慮した動的応答解析

#br

#br
#br



*Marc_mentat [#i623feff]
**動的応答解析 [#c562fb03]

-C1-C2-C3が破断~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/dynamicC1-C2-C3.gif
#br
#br

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/tsunoda/dy01.png
#br
#br
**Model-600 [#wf6637cf]
主径間長600m、 900mの鋼斜張橋に関し、以下の手順に似て崩壊解析を行う。~
Part 1: 常時設計荷重を満足するようなモデル橋を設計する。解析は弾性とし、断面照査は許容応力法を用いる。~
Part 2: モデル橋の終局強度および崩壊形態を明らかにする。解析は弾塑性解析とする。~

#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)~

#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面吊り斜張橋であるが、片面のみを考慮し、平面モデルとして解析する。~






#br
#br

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


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


1)アンカーケーブル(C1)が破断したケース~

(a)C1が破断~



スパン中央の鉛直変位:2.105mm~
隣接するC2の応力は986N/mm2, ~
C14の応力は502N/mm2, ~
主桁、主塔、ケーブルの応力すべて弾性範囲内~
#br
#br

(b)続いてC2が破断したと想定~

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

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






**ケーブル破断を考慮した静的解析 [#s98f4a60]
ケーブルが破断すると想定して,ケーブルを無くした構造系に衝撃荷重と設計荷重を載荷させ,部材ごとの応力が~
各部材の降伏応力に達しているか否かの照査を行う.



1)アンカーケーブル(C1)が破断したケース~

(a)C1が破断~


スパン中央の鉛直変位:2,105mm~
隣接するC2の応力は827N/mm2, ~
C14の応力は424N/mm2, ~
主桁、主塔、ケーブルの応力すべて弾性範囲内~
#br
#br

(b)続いてC2が破断したと想定~

スパン中央の鉛直変位:2,107mm~
隣接するC3の応力は1094N/mm2, ~
C14の応力は257N/mm2, ~
主桁、主塔、ケーブルの応力すべて弾性範囲内~
#br
#br

(c)続いてC3が破断したと想定~

スパン中央の鉛直変位:2,620mm~
隣接するC4の応力は1330N/mm2, ~
C14の応力は64N/mm2, ~
C4はケーブルの第1降伏点に達する.主桁,主塔,の応力は弾性範囲内.~
#br
#br


#br
#br

#br

**弾塑性崩壊解析[#a230747d]
-材料構成則

(a)主桁を構成する鋼材は SM400、主塔を構成する鋼材は SM490Yとする。構成則は~
バイリニア、初期ヤング係数は 2.0 x 105 N/mm2 とする。~
σy =~
245 MPa (SM400)~
315 MPa (SM490)~
355 MPa (SM490Y)~
#br
#br
(b)ケーブルは ST1570の平行線ケーブルとする。構成則はトリリニア、初
期ヤング係数は 1.95 x 105 N/mm2 とする。~
σy = 1160 MPa~
σp = 1570 MPa~
E1 = 155 GPa~
E2 = 55 GPa~
εy = 0.75 %~
εp = 1.50 %~
εb = 4.0 %~
#br
#br
(c)主桁、主塔断面はファイバー要素に分割する~

板厚方向に2分割、板幅方向は6分割 した 。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/M300/p21.png~

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


ケーブル部材は主桁および主塔にピン接合する。すなわち、ケーブル部材端は回転自由
とする。~
#br
#br
(d)解析方法~
最初にケーブルプレストレス(Pr)を載荷し,死荷重および活荷重 D+L)を漸増させる.すなわち,~
P = Pr + k ( D + L )~
である.kは荷重増加係数であり,0.01ずつモデル橋が崩壊するまで増加させる.終局~
時の kを終局荷重係数 kuとする.~

#br
#br
#br
&size(17){【解析結果】 }; ~

#br
-変形前~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/M300/p37.png~

#br
#br
#br
#br
-変形図(梁ピン)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/M300/p36.png~

#br
#br
#br
-変形図(梁ピン)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/M300/p35.png~

#br
#br

#br
-変形図(自由度固定)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/M300/p31.png~

#br
#br
-変形図(自由度固定)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/M300/p32.png~
#br
#br


#br
#br
-変形図(ばね)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/M300/p42.png~
#br
#br
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/M300/p41.png~
#br
#br
-変形図~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/M300/p29.png~


#br
#br

-グラフ~
青線:中村先生の解析結果 オレンジ線:Marc~
縦軸:k 荷重増加係数

&size(16){(1) 梁ピン結合};
#br

-主桁中央 鉛直変位~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/M300/pin01.PNG~
#br

-主塔頂部水平変位~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/M300/pin02.PNG~

#br
-主桁下縁応力 スパン中央~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/M300/pin03.PNG~

#br
-主桁中央応力 スパン中央~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/M300/pin04.PNG~

#br
-主桁下縁応力 主塔位置~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/M300/pin05.PNG~

#br
-主桁中央応力 主塔位置~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/M300/pin06.PNG~

#br
-主塔縁応力 基部~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/M300/pin07.PNG~

#br
-主塔中央応力 基部~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/M300/pin08.PNG~

#br
-C1・ケーブル応力~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/M300/pin09.PNG~

#br
-C3・ケーブル応力~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/M300/pin10.PNG~

#br
-C14・ケーブル応力~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/M300/pin11.PNG~

#br
#br

&size(16){(2) 自由度固定(鉛直方向)};
#br
-主桁中央 鉛直変位~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/M300/ziyudo01.PNG~
#br

-主塔頂部水平変位~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/M300/ziyudo06.PNG~
#br

-主桁下縁応力 スパン中央~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/M300/ziyudo02.PNG~
#br

-主桁中央応力 スパン中央~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/M300/ziyudo03.PNG~
#br

-主桁下縁応力 主塔位置~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/M300/ziyudo04.PNG~
#br

-主桁中央応力 主塔位置~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/M300/ziyudo05.PNG~
#br

-主塔縁応力 基部~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/M300/ziyudo07.PNG~
#br

-主塔中央応力 基部~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/M300/ziyudo08.PNG~
#br

-C1・ケーブル応力~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/M300/ziyudo09.PNG~
#br

-C3・ケーブル応力~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/M300/ziyudo10.PNG~
#br

-C14・ケーブル応力~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/M300/ziyudo11.PNG~
#br
#br

**断面照査 [#a0f86efb]
-主桁~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/M300/p7.png~

断面性能を下記のように求める。
,,,A(mm²),y(mm),Ay²(mm⁴)
,UFlag,11000×12,132.000,-1256,208234752000
,2-Web,2500×12,6.000,0,31250000000
,LFlag,11000×16,176.000,1258,278531264000
,Total,,368.000,,518016016000

#br
#br

,I(断面2次モーメント),509610745739
,l(有効座屈長),20000
,r(断面2次半径),1177

#br
#br

主塔部~
,N=,-15969kN~
,sn(軸応力)=,-43.4N/mm²~
,My=,-25130~
,sb(曲げ応力)=,-69.9~
,sna(軸許容応力)=,113.3~
,sba(曲げ許容応力)=,140~
,r=(照査値),-0.809~


スパン中央~
,N=,239kN~
,sn(軸応力)=,0.6N/mm²~
,My=,39145~
,sb(曲げ応力)=,108.9~
,sna(軸許容応力)=,110~
,sba(曲げ許容応力)=,140~
,r=(照査値),-0.782~




#br
#br

-主塔~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/M300/p8.png~

断面性能を下記のように求める。
,,,A(mm²),y(mm),Ay²(mm⁴)
,UFlag,3000×22,66.000,-2021,266.911.986.000
,2-Web,4000×22,176.000,0,234.666.666.667
,LFlag,3000×22,66.000,2021,266.911.986.000
,Total,,308.000,,768.490.638.667

#br
#br

,I(断面2次モーメント),768.490.638.667
,l(有効座屈長),54.600
,r(断面2次半径),1580

#br
#br

主塔基部~
,N=,-28100kN~
,sn(軸応力)=,-91.2N/mm²~
,My=,16642~
,sb(曲げ応力)=,-43.8~
,sna(軸許容応力)=,126.4~
,sba(曲げ許容応力)=,140~
,r=(照査値),-1.034~
#br
#br
-ケーブル応力~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/M300/p20.png~

#br
#br
**Model-300 [#k5e303bd]

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


***モデル橋の設計(断面照査) [#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~

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









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


中村先生(メモ)~
節点温度~
・A、B、C全てに‐100℃をかけると、A,B,Cそれぞれ2倍縮んでいる。~
・片側の節点だけだと、半分の値だけしか縮まらない~
・ある要素に温度応力をかけるときは、両側の節点に節点温度を与える必要がある。~
→完了(2節点での温度勾配を考慮するのではないか)~

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


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

#br
#br

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



中村先生(メモ)
-Model-300については、case3(バネ要素)で合っている(中村先生も同じ方法をとっている)~
桁と主塔の節点が同一点上にあるのは、実橋を再現するという点からもおかしい。~
その上で結果(今回は鉛直変位)に差が出るのはしょうがない(解析ソフトによる違い)~
-バネ要素が本当に自由度を制限できているかを簡単なモデルを作成して調べる。
-断面力を調べる
桁の中央・主塔の交差点・最下点の軸力と曲げモーメント~
-Model-300の弾塑性解析とケーブルがないと仮定したときの解析


#br
#br


&size(17){【解析結果】 }; ~
#br
#br



a) ケーブルプレストレス載荷時( Pr スパン中央の変位:+706mm (中村先生: :+706mm)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/M300/p1.png~

#br
#br

b) 死荷重載荷時( Dスパン中央の変位:-1,710mm (中村先生 :-1,712 mm)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/M300/p2.png~

#br
#br

c) カウンターウエイト載荷時( CWスパン中央の変位:+254 mm (中村先生 : +261 mm )~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/M300/p3.png~

#br
#br

d) 集中活荷重載荷時( LP1スパン中央の変位:-491 mm (中村先生 : -491 mm )~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/M300/p4.png~

#br
#br

e) 分布活荷重載荷時( LP2スパン中央の変位:-159 mm (中村先生 : -162 mm )~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/M300/p5.png~

#br
#br

f) Pr + D + CW + LP1 +LP2スパン中央の変位:-1,400 mm (中村先生 : -1,398 mm )~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/M300/p6.png~


#br
#br
#br



**ファイバー要素 [#u0fdf28d]
-断面をメッシュ分割し,各メッシュに応力ひずみ関係を与える.非線形特性をより正確に表現できる.~
軸力変動や二軸曲げによる応力変動を比較的正確に表現できる.~
-解析ソフトによる結果の違いはある程度許容する。
#br
#br




**弾塑性解析(材料非線形) [#t09f986d]

以下の設定が必要になる~
-テーブルでtime,応力ひずみ関係の入力~
-材料特性で降伏点の設定~
-境界条件でtimeを入力~
-荷重ケースで解法の設定~

#br
#br


1. テーブル~
1.1 time(荷重を増分荷重にする)~
Type:time~
Formula:1.0v1~

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




#br
#br


1.2 応力ひずみ関係
Type:eq_plastic_strain~
Point:任意


2. 材料特性で降伏点の設定



3. 境界条件



4. 荷重ケース

2.ケーブル・モデル:~
長さ1mのケーブルとする。~
B点に鉛直荷重P = 1.0 kN を増分載荷する。~

断面: 半径 10mm の円形とする。~

構成則に関しては、はりと同じ計算方法で計算を行い入力した。~

解析結果~
1)自由端の鉛直変位と荷重の関係~

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


2)ひずみと荷重の関係~

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

3)応力とひずみの関係~

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


全てのグラフにおいて中村先生の結果と一致した

#br
#br
はり・ケーブル単独の弾塑性解析~
はりモデル:~
・長さ1mの片持ち梁とする。~
・自由端に鉛直荷重P = 1.0 kN を増分載荷する。~
・部材断面
・はりの断面は10 cm x 10 cm の正方形とする。5×5のファイバー要素に分割した。~
・はりはSM400とする。初期ヤング係数は 2.0 x 105 N/mm2 とする。~

構成則入力の変更点:~
テーブル(タイプ: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~

その他の設定~
・反復法:完全ニュートンラプソン~
・収束判定 残差と変位(相対許容値):0.001~

解析結果(以下 青線:中村先生の解析結果、オレンジ線:Marcでの結果)~
1)自由端の鉛直変位と荷重の関係~

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

2)断面応力(Layer1)と荷重の関係~

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

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



**温度応力 [#sb33b144]

ケーブルプレストレスを導入できるかの検討~

・初期条件▸新規(構造)▸仮想温度~
→×(この設定だとうまくいかない)~
                   
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2108/ka1.png~
#br
・ケーブル要素は用いない(高度な計算になり複雑になるため)~
・代わりに梁要素に温度荷重を与えてケーブルのプレストレスとする~
・ソリッドモデル(斜張橋)で構造解析するのは、計算容量が大きくなるため骨組みモデルで計算する~
#br
#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)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~

・節点温度を与えると、節点に接している要素全体に温度応力がかかる~
・熱膨張係数を与えている部材のみ、温度応力がかかる~
(熱膨張係数が与えていなければ節点が共有されている要素には温度応力はかからない)~
#br
・境界条件▸節点変位 で本当にケーブルに対して温度応力を与えることができているか~
手計算と比較、片持ち梁で段階的に節点温度を与えて挙動をみる~
#br
#br
#br



*Salome_meca [#cfa5c1ef]

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

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

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

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



線形解析→非線形解析(弾塑性解析)~
・(弾塑性)梁モデル→multfiber?(現在×)~
・(弾塑性)ソリッドモデル→〇~
・(弾塑性)ソリッドと梁のモデル→× salome_mecaでは弾塑性解析をするとき、すべての部材(異なる材料があるとき)に弾塑性の設定を与える必要がある(仮説)

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

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

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


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


#br
#br
#br

**斜張橋モデル [#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]

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

#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

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

・ケーブルの断面形状を円柱→四角形~


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
・結果(アニメーション)~
300kN~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tunoda/dep2.gif


#br

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



#br
#br


**温度応力 [#f0292266]
#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
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/tsunoda/2105/scr16.png~

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

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


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~




#br
#br
#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

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

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
モデルが(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

・(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]


熱応力・梁要素(ケーブル)~
片持ち梁【集中荷重・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】~
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)




#br


*メモ [#yb1a7ac2]

#br

**参考文献 [#g936ac10]

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

#br
#br


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





#br
#br


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