千代岡の修論日誌
をテンプレートにして作成
[
トップ
] [
新規
|
一覧
|
検索
|
最終更新
|
ヘルプ
|
ログイン
]
開始行:
#contents
[[自分の卒論日誌:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?%E5%8D%83%E4%BB%A3%E5%B2%A1%E3%81%AE%E5%8D%92%E8%AB%96%E6%97%A5%E8%AA%8C]]はこちら
研究タイトル：部分的に腐食したケーブル系橋梁の連鎖崩壊
*研究目的 [#ya05820f]
修士2年~
編集中
修士1年
4年次ではSalome-Mecaを用い,簡易モデルで衝撃力の再現が妥当か検討を行なった.ただ,モデルを作成するにあたり,ピン接合の再現や幾何学的非線形絡み（まだ解決してないところもあるが...）の問題があり,それらをできるようにした.~
修士では,これまでMarc_Mentatを用いて橋モデルを作成して解析を行なっていたものをSalome-Mecaでもできるようにする.そのためにまずはSalome-Mecaで橋モデルを作成して,解析できるところまで持っていく.（これが第1ステップ!）~
そして,ケーブルに腐食を考慮した斜張橋に動的応答解析を行い,衝撃力の大きさを検討していく.（簡単に言うと,昨年,簡易モデルでやったことを橋モデルで行う）~
斜張橋モデルに関しては,ケーブルの張り方や本数,橋梁のスパン長や減衰率等を色々変えた上で解析を行う.
→詳細記述~
ケーブルを有する橋梁のケーブルが腐食することで橋が落ちてしまう事故が近年各地で起きています.（例えば...台湾の南方澳橋とか.卒業された先輩がその橋を対象とした研究を行なっていました.）~
その中で,1本目のケーブルが切れることで,次々とケーブルが切れてやがて落橋に至る連鎖崩壊をターゲットに研究を行なっています.~
この現象は橋梁の落橋などを引き起こす重大な事故に繋がるため,そのメカニズムの研究も進められている現状です.~
本研究室ではその解析手法の提案が行われてきたが、ケーブルが切れた際の衝撃力の大きさや他のケーブルに掛かる荷重の再分配率が明確ではありません...。~
卒論では,簡易モデルと呼ばれる簡単なモデルで衝撃力の検討を行ないました.そこで,修論では斜張橋モデルを使って手法の妥当性を検討する目的で,健全と腐食ケーブルを有する斜張橋に対してケーブルが切れた際の衝撃を衝撃力として数値解析的に与え,その大きさを検討していきます.~
モデル自体もケーブルの張り方やスパン長等の性質・値を変化させていくことでどう変わっていくのかも見ていきたいなと考えています.~
後は,Salomeのビーム要素絡みの問題も色々解決していきたい。
*今後の予定 [#bd1b97fd]
Salome-Mecaの梁要素絡みの問題は色々あるため,卒業するまでは後輩達とともに1つずつ解決していきたい。~
12月~
12/23 中間発表第2回目~
12/26 大学生活 最後の帰省
夏休み~
〇研究~
・ビーム要素におけるモーダル解析の理解~
・橋梁モデルでC1に着目した解析(C1健全時が切れたときの応答値,C1部分腐食時の切れたときの応答値)~
→そのときのC13張力は？
〇その他~
7月~
・橋梁モデルで腐食を考慮した場合の終局強度を求める。~
・交流会の準備~
・今取り組んでいるe-learningの英語の発表準備~
やることが多いけど...何とか乗り切ろう.~
6月~
・橋梁モデルで部分腐食を考慮した解析~
簡易モデルは卒業。~
一応,Salomeでやることになっているが,どうしてもできないときのためにMarc(midas)で取り組むことも視野に入れて行動したほうがいいかも。~
自分のやるべきことをやった上で取り組むなら,誰も何も言わないだろう。
コメント：時間は限られているけど,焦らず1つずつ進めていこう。
4月,5月~
・部分腐食モデル：分布荷重での解析~
→終了後,集中・分布・衝撃荷重の応答値変化をまとめたい。
・ブリッジコンペのモデル作り~
解析はSalome(ダメならmidas),モデル作成はfreecadで行う。（作ったモデルはSalome,midasに持っていけるらしい。Salomeは確認済み）
コメント：いよいよ最終学年！早速やることが多いけど...これまでやってきたことを修論で発表できるように頑張ろう！
*Model300をMarcと解析結果比較[#n498dba5]
**モデル作成 [#je113eea]
Salome-Mecaで橋モデル（model300）を作成してみる.~
以下の写真らはMarcで卒業された先輩が作成した橋モデルである.これを参考にして,まずはモデルの完成を目指す.（Marcではm系でモデルを作成していたため,Salomeも同様に作成する.）~
角田さんと秋山さんの論文を読んで,内容を理解することからスタート.(6月初めには結果を出したい.)
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/6honmodel300xy.PNG
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/3D6honmodel300.PNG
・ひとまず,Salome-Mecaでモデルの作成は終えた.メッシュも切り終えたので,Asterstudyの設定方法を考える.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/Salomemodel300.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/Salomemodel300yoko.png
Salomeでの解析方法~
まず,モデルが完成したら,静的解析（ダメそうなら,時刻歴応答解析で）で主桁中央鉛直変位のデータを取ってみてMarcと比較する.(Marcで導出した変位の値は-1.17m)~
材料非線形を考慮しているため,AnalysisでSTAT_NON_LINE（動的ならDYNA_NON_LINEかな。）を用いて解析を行う.~
また、プレストレスは一旦温度応力で試してみようと考えている.~
**解析記録 [#gb613b11]
モデル中央主桁鉛直変位を求める。~
2/26~
ケーブルのみを材料非線形を考慮し、静的解析を回したところ、結果が得られた。（動的解析だと,step数を細かくする必要があるエラーメッセージが出て、解析が回らない。）~
得られた主桁鉛直変位は-25.1263m~
以前、トラス要素を用いた場合の主桁鉛直変位は-1.34mとMarcと近い値（-1.17m）が出ていたのだが...。~
静的解析だと普段使っている温度によるプレストレスを入れられない説もあったため、動的解析（100sで荷重係数1となるように回したが,60sで解析停止）で解析が回る区間の結果を出力したところ、解析途中(60sのデータ)ですら主桁鉛直変位-15.2196mが得られた.
20sごとに約5mごと鉛直変位が下がっていることを踏まえると、静的解析で得られた結果はSalomeの解析では妥当なものであると言える？（温度も入ってそう.）
...ということは、Marcの梁要素とSalomeのBAREE（トラス要素）は似た性質を持っている可能性がある？~
ただ、簡易モデルでも記載したが、ケーブルは曲げ剛性とねじれ剛性を無視するため、それが悪さしている可能性もあり得る.~
→原因は桁とケーブルを接続するバネの剛性が低いことが原因だと判明。主桁中央鉛直変位は-1.50mと結果が得られ,MarcやケーブルBAREE要素と近い値が得られた。~
6/26~
前回,主桁が持ち上がった原因として,ケーブルの半径の与え方を間違えていたことだと判明した。（この橋梁モデルはm単位で作成しているのに、ケーブルの半径の値をm換算するのを忘れていた。）~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/model300-2.png
得られた主桁鉛直変位は-1.34mだった.（Marcで導出した変位の値は-1.17m）
現在、与えている境界条件~
死荷重：桁：122000N/m,主塔：60000N/m （秋山さんの論文と残されたMarcのデータを参照）~
活荷重：主桁中央に集中荷重LP1=1.175×10^6N, 桁全体にLP2=35250N/m (Marcのデータを参照）~
→活荷重の計算方法：道路橋示方書のL荷重(B活荷重)の計算手法に従う.~
L荷重：大型車を連行させた分布荷重にしたLP1と大型車以外の自動車荷重（その他の部分）を分布荷重にさせたLP2の合計.(LP1は集中荷重に直した.)
L荷重の図を見たほうがわかりやすいと思うので,[[この検索ページから:https://www.google.com/search?q=L%E8%8D%B7%E9%87%8D+%E5%9B%B3&client=ubuntu-sn&sca_esv=761560977b3b699f&channel=fs&biw=1920&bih=906&ei=ORo0aMeAA7eUvr0PnsWlsQg&ved=0ahUKEwiHgYSp0MCNAxU3iq8BHZ5iKYYQ4dUDCBA&uact=5&oq=L%E8%8D%B7%E9%87%8D+%E5%9B%B3&gs_lp=Egxnd3Mtd2l6LXNlcnAiC0zojbfph40g5ZuzMgUQABjvBTIIEAAYgAQYogQyBRAAGO8FMgUQABjvBUjcO1DuGViGNXAEeAGQAQCYAYQBoAGIBqoBAzEuNrgBA8gBAPgBAZgCCqACigbCAgoQABiwAxjWBBhHwgIFEAAYgATCAgYQABgFGB7CAgUQIRigAZgDAIgGAZAGCpIHAzQuNqAHlgyyBwMwLja4B-QFwgcFMi02LjTIB0k&sclient=gws-wiz-serp#vhid=anPAULIdubmghM&vssid=_QRo0aNzDJ9a2vr0PwZPnAQ_42]]調べてみてください.~
LP1=(10×5.5×10)+{5.0×(18-5.5)×10}=1175kN~
(10×5.5×10):主載荷荷重×載荷幅×載荷長~
{5.0×(18-5.5)×10}:{従載荷荷重×(舗装部分の長さ-載荷幅)×載荷長}~
従載荷荷重は主載荷荷重の半分.~
LP2=(3.0×5.5)+{1.5×(18-5.5)}=35.25kN/m~
(3.0×5.5):主載荷荷重×載荷幅~
{1.5×(18-5.5)}:従載荷荷重×(舗装部分の長さ-載荷幅)~
CW: 80000N/m (秋山さんの論文を参照)~
ケーブルプレストレス：C1,C26：-150℃ それ以外のケーブル：-100℃ (Marcのデータを参照）~
6/19~
ケーブルをトラス要素に変えたら解析が回ったが、主桁中央鉛直変位が1.00mという解析結果が得られた。~
論文での結果を見てみると、下にたわんでいる結果が得られていたため、荷重が上手く与えられていないのでは？
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/model1.png
・死荷重を入れるのを忘れていた...。Marcで作成していたモデルを見ると,重力を（重力荷重？）モデル全体に与えていたようだが...。（密度×体積×重力加速度で導出ってことだと思う。）~
ただ、桁にかけている死荷重のデータは概要に書いてあるため、早速修正する。
・バネでのやり方だと中々解析が回らないため、やり方を変えてバネを外してケーブルをトラス要素に変えてみた。~
トラス要素：軸力のみが生じ,部材と部材の共有節点をピン接合した要素。~
自分の卒論で取り扱った簡易モデルのケーブルをトラスに変えてみたところ、桁の先端変位（2.50mm）とケーブル張力(16N)がバネ使用時と一致し、有限要素解析でケーブルをトラス要素として扱っている例をいくつか見たので、トラス要素を使うこと自体、信頼性はあると思うが、ケーブルを要素分割した際に解析は回るのかな？（トラス要素の定義に従うなら、要素分割したらケーブルとケーブルをピン接合するということになると思うが...）
6/6 問題になっている可能性が高い箇所を見つけた.Analysisの部分で自由度に関する問題がノード199(C21と塔を繋ぐバネ部分)で発生しているらしい。~
バネ絡みの設定がおかしいと思うが...。~
思いついた解決策として,バネとケーブルの共有節点も離散要素（マトリクスを与える状態、要するにバネ要素と同じ設定）を与えてみるとか？
↓赤丸で囲んだ箇所が該当箇所~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/N199-1.png
↓赤丸の拡大図~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/N199.png
↓メッセージ内容~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/error3.png
6/4 現時点で出ているエラーメッセージ Analysis,DYNA_NON_LINEでエラーが発生中~
(1)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/bridgeerror2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/JPbridgeerror.png
5/25 何とか解析回すステージには立てた.（時刻歴応答は自分,静的解析は山本君にお願いしている。）~
ここからエラーメッセージを順番に潰していく.
5/23 試行錯誤しながらプログラムを組んでいる最中。~
何とか今週までに解析回せるところまで行きたい。
*影響線解析 [#h9310a98]
影響線：移動荷重に対してある点に生じる物理量の大きさを移動荷重の位置関数で表したもの~
影響線の定義から,まずは移動荷重を求め方が分からないと話にならない？
↓橋の上を走る電車を例に移動荷重を算出する方法を見つけた。荷重条件に式を与える方法を使って,動的解析をする必要があるらしい？~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/movingload.png
[[移動荷重について:https://forum.code-aster.org/public/d/13262-how-to-calculate-the-moving-loads-on-a-bridge-in-code-aster]]
[[移動荷重について_その2:https://forum.code-aster.org/public/d/19546-unilateral-constraint-and-moving-force-with-a-given-velocity/8]]
AFFE_CHAR_MECA_F→FORCE_POUTREを使い、DYNA_NON_LINEなどのコマンドを使う動的解析でできそう？~
ちなみに...midasだと[[説明書の10章:https://in.midasit.co.jp/wp-content/uploads/2020/07/constrcution-83-civil-book.pdf]]を見たところ,影響線を導出する手法が記載されていた.
*model300に部分腐食を与えた終局強度解析（後期） [#ce4ff1df]
**静的解析での検証 [#u4e27946]
○線形解析の場合
これまで解析が収束しなかったのは,ケーブル破断後の主桁応力変動が解析上では計算しきれないため起きた可能性がある.~
そこで桁のみを線形状態にした場合で部分腐食を与えた解析を取り組んでみた.
これまでk=3.0になるまで荷重を載荷していたので,k=3.0の荷重を静的荷重(k=1.0)と定め,その荷重に向けて徐々に荷重を増やしていく.
-k=0.9(設計荷重の2.7倍の力が掛かっている.)で解析停止.解析停止までの様子は以下の通り.
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/linearm300.gif
k=0.9(設計荷重の2.7倍の力が掛かっている.)時の橋梁鉛直変位~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_linearzentai.png
-主径間中央_主桁の応力変化をまとめた.~
下縁部~
k=0.54で降伏応力を超えている.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/linearsitaketa.png
中立軸~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/mannakaketa.png
上縁部~
k=0.60で降伏応力を超えている.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/linearketaue.png
こちらでも同様に下縁部→上縁部の順に降伏応力に達している.~
しかし,線形解析でも中立軸は降伏に至ってない.
橋長に沿った主桁中立軸の応力状況~
側径間右側で一番大きな応力が出ている？~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/tyurituketazenbu.png
・線形解析で得られたこと~
①大きな設計荷重を掛けても中立軸（付近）は降伏応力に達しない.(少なくとも,1本のケーブルが腐食しただけでは.)~
これを他のケーブルにも腐食を与えたり,連鎖的にケーブルが切れるといった挙動が起きるとまた結果は変わってくるはず...
②主塔の真下（橋長420m）で大きな圧縮応力が出ている.
これまで主径間中央で最も大きな応力が出ていると思ったが,この考えは違うかもしれない？~
**解析停止時の応答値 [#b9db6bd9]
線形解析での結果を踏まえ,主塔付近の応力を確認するため,解析が停止するまでのモデル全体のアニメーションと主桁断面の応力状況等をまとめた.
C1上端部分腐食時のアニメーション（腐食場所関係なくモデルの挙動はほとんど同じ.）
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/koukiC1uecoro.gif
○橋長方向と鉛直変位の関係~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/k185DEPL.png
○主桁下縁部応力と橋長方向の関係~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/k185_ketasitaSIGM.png
○主桁中立軸応力と橋長方向の関係~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/k185_ketaSIGM.png
-C13部分腐食の結果
○橋長方向と鉛直変位の関係~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/k170_ketaDEPL.png
○主桁下縁部応力と橋長方向の関係~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/k170_ketasitaSIGM.png
○主桁中立軸応力と橋長方向の関係~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/k170_ketaSIGM.png
**腐食状態ケーブル（腐食率50%）をモデルから抜いた場合の解析（連鎖崩壊解析手法を用いた） [#r485660d]
手法は連鎖崩壊解析と同様.これまで破断させる応力ひずみ関係を腐食させたケーブルにのみ与えていたが,他のケーブルにもその関係を与えて解析を回す.
○C1なしモデル
250sでC1ケーブルを破断させ,260sまで解析を回した.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/chainC1hadan.gif
先行研究では,252s以前で塑性ヒンジが見られたため,252sまでの応答値を確認した.
○251.55s時の変形図~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/251_55_DEPL.png
○251.55s時の主桁応力~
最初に,C3,C4間の主桁中立軸で降伏応力に達したことを確認.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/chain_ketaSIGM.png
○C1側側径間ケーブルのひずみ~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/chain_epsi.png
○C13なしモデル
*model300に部分腐食を与えた解析（前期） [#w3f7c4d5]
**Salomeでの解析 [#yb204e33]
終局強度を考慮した解析を行う.(先輩の事例を見ると,静的解析をやっていたのでとりあえずそっちで。)~
終局強度とは？:構造物が破壊（崩壊）する直前の強度。~
終局強度の式は以下の通りである。~
P=Pr+k(D+CW+L)~
ただ,過去の研究ではP=Pr+k(D+L)で終局強度を求めている場合も。~
どっちを使うべきだ？~
Pr:ケーブルプレストレス,D:死荷重(桁),CW:カウンターウェイト,L:活荷重,k:荷重増加係数~
まず,ケーブルプレストレスを載荷し,その後に他の荷重をkを増やすことで増加させていく。~
構造物崩壊時のkを終局荷重係数とする。→この係数を求める！
→これをC1下端,中央,上端に部分腐食を与え,主桁中央変位と腐食部分のC1張力と応力ひずみを見る。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/model300partialcoro1.png
断面形状は以下に示す.~
主桁~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/keta.png~
塔（下部）~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/tou.png~
塔（上部,ケーブルと繋がっているところ）~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/hutozu.png~
***ケーブル腐食なしの終局強度(幾何学的非線形考慮しない) [#h30a0189]
・主桁中央変位と荷重係数の関係.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/tyuouDEPL_nonlinearonly.png
・主桁中央応力と荷重係数の関係.~
K=1.81で下縁部降伏,k=2.60で中立軸付近降伏,k=2.64で解析ストップ~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/nonlinear_stress.png
***ケーブル腐食なしの終局強度(幾何学的非線形考慮) [#gbb3d070]
部分腐食の解析を行う前に,腐食なしで先輩方が解析で求めた終局強度と同じ(近い)値を調べた。~
このとき,終局荷重係数はK=2.63だったが,Salomeだとk=0.30程度で解析エラーが出てしまう。(荷重係数が1になる前に解析が止まってしまうのは流石におかしい。)~
解析が停止する原因を調べたところ、主桁に材料非線形を考慮しない場合（線形弾性）は解析が回ったため、主桁に関する設定が何かしら悪いことをしていると思われる。~
断面のメッシュの切り方に問題がある([[参照ページ:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?Salome-Meca+%E3%83%93%E3%83%BC%E3%83%A0%E8%A6%81%E7%B4%A0#jfa31dbe]])可能性もあるので、メッシュを切り直して再度解析を回す。~
(主桁断面は以下の写真,左側の形状.秋山さんの修論日誌から引用.)~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/akiyama/syutiudanmen2.png~
塔・桁・ケーブルに幾何学的非線形を考慮.
解析方法~
ケーブルプレストレスを荷重係数1まで、他の荷重は荷重係数を1を超えても載荷し続ける。~
〇解析結果
[[比較データ:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?%E8%A7%92%E7%94%B0%E3%81%AE%E4%BF%AE%E8%AB%96%E6%97%A5%E8%AA%8C#f0558fbf]]~
先輩が解析した際の崩壊までの挙動（角田さんの修論日誌より引用）~
k=2.64のとき，崩壊した．
主桁および主塔の変形は，側径間のケーブルC4がk=1.90にて第一降伏点に達するまでは，
線形増加する．k=2.24で主塔位置の主桁下フランジが，k=2.26にてスパン中央の主桁下フランジが降伏する．
その後，スパン中央の主桁断面の塑性化が進む．そして，k=2.61で側径間ケーブルC3，C4が第二降伏点に達し，
k=2.62で，スパン中央の主桁が全断面降伏し塑性ヒンジが形成され，橋全体が崩壊する．
まとめると、k=2.24で主塔位置の主桁下フランジ、k=2.26で橋梁中央部の主桁下フランジが降伏,K=2.62で橋梁中央部の主桁下フランジに塑性ヒンジが形成。~
解析ソフトの差異で、kの値に多少違いは出ると思うが、大きくずれたらSalomeの結果が間違っている。~
7/9~
結果まとめ(比較図を載せたかったが、先輩の論文に載っているものなので割愛.)~
・荷重係数1での橋長鉛直変位(Salome:-1.098m,Marc:-1.129m)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/model300_k1DEPL.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/Marcmodel300DEPL.png
・主桁中央変位と荷重係数の関係.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/model300tyuouDEPL.png
・主桁中央応力と荷重係数の関係.~
K=1.83で下縁部降伏,k=2.51で中立軸降伏,k=2.60で解析ストップ~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/model300syukyoku_k_SIGM.png
Marcに近い挙動の結果は得られた.
7/2~
新しい結果が得られたので、そちらを載せる.~
k=2.61で解析エラー(k=2.61で橋梁モデルが崩壊.)
主塔の死荷重とケーブルプレストレスはk=1.0を維持し,解析.~
それ以外の荷重はkを増やし続ける.~
角田さんの結果と比較して...~
・主桁下縁部降伏応力は,ほぼ一致(k=1.8,私の結果k=1.9)~
ただ,全断面降伏(塑性ヒンジ形成)は私の結果がk=2.16に対し,角田さんの結果はk=2.63.~
ここでずれたのは何故？(塔の死荷重の有無,CWの値を新しいモデル用に合わせているのと,当時のモデルはアンカレッジケーブルがあったため、それらの違い?：当時の側径間ケーブル7本,)~
→3つめに関しては,側径間ケーブルを1本なくした代わりにC1,26の張力を強めにしているので,あまり大きな影響はなさそうな気がする.~
→追記,桁断面の中立軸は原点ではなく,高さ方向から少しずれた位置にあることが判明.(原点中立軸なのは,上下左右対称断面のとき.)~
・C1応力変化はほぼ一致.~
○変形図(scalefactor５倍)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/m300k.gif
○変位図(k=1.0のとき)~
中央鉛直変位は-1.01m,(Marcのmodel300:-1.129m)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/model300k1_heni.png
Marcの場合：-1.129m~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/Marcmodel300DEPL.png
○主桁中央鉛直変位(kの変化)~
崩壊寸前の鉛直変位は8.37m.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/model300k_DEPL.png
○主塔水平変位(kの変化)~
崩壊寸前の水平変位は1.53m.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/model300k_tou_DEPL.png
○モデル中央,主桁応力(k=1.9で下縁部降伏,k=2.16で全断面降伏)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/model300k_stress.png
○モデル中央,主桁下縁ひずみ(右図,降伏ひずみ付近を拡大)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/model300k-EPSIsitaall.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/model300k_epsi_sita.png
○モデル中央,主桁中立軸ひずみ(右図,降伏ひずみ付近を拡大)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/model300k_EPSIceall.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/model300k_epsi_mannaka.png
○C1応力~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/model300k-C1SIGM.png
6/30~
主桁が降伏に至らない問題は、どうやら桁と塔に与えている応力-ひずみの降伏ひずみが違っていたことが原因だと分かった。~（連鎖に関する過去の論文、角田さんの論文から応力ひずみの線形関係の直線部分の傾きから降伏ひずみを計算して判明。）~
...この1,2週間で色々問題点が見つかったな。これで自分の求めたい結果が得られるか？
6/26~
新たな問題に直面.載荷条件の順番を変えることで異なる変位や応力が出てくる。詳細は[[こちら（載荷条件について）:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?Salome-Meca+%E3%83%93%E3%83%BC%E3%83%A0%E8%A6%81%E7%B4%A0#q683524a]]~
これは早急に何とかしないとな...
6/25~
再度、設定を見直したところ、以下2点の変更点を加えた。
・ケーブルと桁を繋ぐバネの剛性が低いため、桁が大きくたわむことが判明した。→変位に対する剛性を高くした。~
・Marcで得られていた張力からプレストレスとして掛けている温度を計算して、得られた温度を各ケーブルに載荷した。（主桁中央鉛直変位の値が-2.5mとMarcで得られていた結果と1mほどずれがあるため）~
これらの変更を加え、ようやくMarcやSalomeのケーブルをBAREEにした場合と近い主桁中央鉛直変位（-1.5m）が得られた。
この状態で終局強度を求めようとすると、k=2.60で収束エラーが出たため、その手前の解析結果を見たところ、主桁中央部がどこも降伏していないことが判明...。（171MPaしか出ていない。）~
また、収束エラーが出たタイミングでモデルが崩壊している可能性が高いが、その様子が見ることができないのも問題。（桁断面のメッシュを細かくしてみる。）~
6/18~
・解析はk=4.17で停止~
グラフは角田さんの論文と同じように作成。
[[角田さん:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?%E8%A7%92%E7%94%B0%E3%81%AE%E4%BF%AE%E8%AB%96%E6%97%A5%E8%AA%8C#f0558fbf]]の結果と異なるものが得られた.(K=2.63)~
主桁に至っては,中立軸付近は降伏していない.~
モデルに関しては,塔の断面の材質（sm490(315MPa)→sm490Y(355MPa)）と断面形状がは違うが、それだけで大きく結果が変わるのだろうか？~
→多少の断面情報で角田さんの結果と大きく変わることはない。~
変形図~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/k_warpbyvector.png
主桁中央鉛直変位と荷重係数~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/k_syukyoku_DEPL.png
スパン中央応力と荷重係数~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/k_syukyoku_SIGM.png
***ケーブル全体腐食の終局強度(幾何学的非線形考慮,要素分割していない.) [#m66d1181]
C1を腐食させる。腐食は50%を想定.このデータを部分腐食と比較.~
・変形図:scalefactor5倍(色付きがC1腐食,白色が腐食なし)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1coro50vsnotcorozentai.png
・主桁中央変位と荷重係数の関係.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coro50DEPL.png
・主桁中央応力と荷重係数の関係.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/SIGMcoro50_sitatyuou.png
上縁と下縁の応力ひずみも見た。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coro50_sitaSIGM.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coro50_ueSIGM.png
***ケーブル全体腐食の終局強度(幾何学的非線形考慮しない) [#h108caf7]
・主桁中央変位と荷重係数の関係.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/allcoro50_DEPL.png
・主桁中央応力と荷重係数の関係.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/allcoro50_SIGM.png
***部分腐食の終局強度 [#pff8fa15]
まず、全ての要素に幾何学的非線形を考慮した場合,簡易モデルで用いた要素分割法(腐食した1要素と腐食していない1要素)で解析が止まってしまった。(腐食ケーブルが降伏に達した時点で解析ストップ.)~
ケーブルだけ幾何学的非線形を考慮する場合(GROT_GDEPとPETITが混在している状態でも),少し荷重が掛かっただけ(k=0.4くらい)で解析が止まってしまう.
試しに,幾何学的非線形を全ての要素から抜いた状態で解析を回したところ,この場合であれば,簡易モデルで用いた要素分割法でも腐食ケーブルが降伏したあとも解析が進んだ.~
幾何学的非線形を考慮するかは一旦置いといて,材料非線形のみを考慮した場合で解析を進める.
モデルの状態~
・モデルの寸法や材料特性はケーブル1要素と同じ.(変更点として,各ケーブルを100要素に分割した)~
・腐食は50%を与えている.~
腐食の与え方は以下の通り.応力ひずみを腐食率ごとに減少させていく方法.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/corrosiongraph.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/corrosion_data.png
ケーブルが腐食により応力が減少する（破断する）という挙動は,静的解析では再現できないみたい.(ケーブルが切れる再現だから,力の釣り合いから外れるから？)~
→静的解析ではk=1.74で解析が停止.恐らく,ここでケーブルが切れていると思われる.
そのため、動的解析を用いて終局強度を求めてみた.~(100sで荷重係数1になるように設定.つまり,1s:k=0.01を想定.)~
○C1下端部腐食部~
k=1.74でC1が破断.(静的解析で解析停止したタイミングでケーブルが切れたことを確認.)~
橋梁中央の変位(270m)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coro50model300DEPL.png
主桁中央の応力と荷重係数~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/keta_coro50_SE.png
C1の健全部の応力と荷重係数~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/sita50coro_C1kenzenSIGM.png
C1の腐食部の応力と荷重係数~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/sita50coro_C1coroSIGM.png
C13張力~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/50sitacoro_C13T2.png
○C1真ん中腐食部~
k=1.74で破断~
橋梁中央の変位(270m)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/mannaka50coro_DEPL.png
主桁中央の応力と荷重係数~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coro50_mannaka_ketaSIGM.png
C1の健全部の応力と荷重係数~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coro50_mananakahadan_kenzenC1SIGM.png
C1の腐食部の応力と荷重係数~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coro50_C1mannaka_coroSIGM.png
C13張力~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coro50mannaka_C13T.png
○C1上端部腐食部~
橋梁中央の変位(270m)~
k=1.74で破断~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coro50_ue_DEPL.png
主桁中央の応力と荷重係数~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coro50ue_ketaSIGM.png
C1の健全部の応力と荷重係数~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/uecoro_C1kenzenSEGM.png
C1の腐食部の応力と荷重係数~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1uecoro50_coroSIGM.png
C13張力~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coro50ue_C13T.png
・分かったこと~
①C1ケーブルが腐食が達した時点でC1健全ケーブルの応力がほとんど変化しないことが確認できた.(何故、こうなるかは考え中.)~
ちなみに,他のケーブルが健全状態の応力ひずみ関係に沿った応力ひずみ変化になる.
②1本切れただけでは塑性ヒンジに達しない.
③ケーブル破断時に他のケーブルに張力変化が得られたことも確認.
④腐食位置毎に破断時の荷重係数に変化はなかった.
結果~
7/21更新~
k=2.05で解析が止まった.しかし,桁は中立軸は降伏していないため,構造物は崩壊に至っていないと言える.
C1腐食部ケーブル応力~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coro50C1SIGM.png
C1健全ケーブル応力~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/kenzenC1SIGM.png
桁の応力~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coro50_ketaSIGM.png
***動的解析における終局強度を求めるC1全体腐食の破断時の応答値 [#i949e212]
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1zentaicoro.png
断面下縁部降伏：k=1.79~
断面上縁部降伏：k=1.91~
ケーブル破断：k=2.3~
k=2.5で解析停止.
・中立軸が降伏する前に解析が停止したが,そのときの降伏応力が231MPaの結果が得られた.~
恐らく,k=2.5付近で全断面降伏すると思われる.
***動的解析における終局強度を求めるC1部分腐食の破断時の応答値 [#f4d17100]
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1bubuncoro1.png
ケーブル上端・中央部に腐食を与えた場合：k=1.74で破断.~
ケーブル下端に腐食を与えた場合：k=1.8で破断.
解析自体はk=1.85で停止.
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1bubuncoro2.png
腐食箇所関係なく断面下縁部・上縁部は降伏応力235MPaに達した後,応答値が変動.~
中立軸は解析停止k=1.85でも降伏応力に達しないため,塑性ヒンジは見られなかった.
***動的解析におけるC13に全体腐食を与えた場合の応答値 [#f4d17100]
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C13zentaicoro.png
断面下縁部降伏：k=1.72~
断面上縁部降伏：k=1.84~
ケーブル破断：k=2.11~
k=2.34で解析停止.
中立軸が降伏応力に至る前に解析停止.しかし,応力変動がC1全体腐食モデルと似ており,C1全体腐食より小さな荷重係数で断面降伏が始まっていることから,k=2.5より小さな近い値で終局強度に至るのでは？
***動的解析におけるC13に部分腐食を与えた場合の応答値 [#rad58a9d]
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C13bubuncoro1.png
k=1.52でケーブル破断.~
k=1.70で解析停止.~
ケーブルが切れるタイミングはほぼ同じだが,変位の応答変化は上端・中央,下端で異なる結果が得られた.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C13bubuncoro2.png
腐食箇所関係なく断面下縁部・上縁部はケーブル破断後に降伏応力235MPaに達する.その後,応答値が変動.~
中立軸は解析停止k=1.70でも降伏応力に達しないため,塑性ヒンジは見られなかった.
***まとめ [#b907df73]
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/kmatome.png
・部分腐食モデルでは,終局強度を求めることは現状できない.(求め方に工夫が必要.)~
・主径間中央に近づくほど,小さな荷重係数でケーブルが切れる.~
・部分腐食モデルのほうが小さな荷重係数でケーブルが切れる.~
**Marcでの解析 [#bc35888e]
3/3~
model300を用いて解析を回してみた.
まず,腐食率50%の破断時の応力ひずみ関係をC1端部にのみ与えた結果を載せる.~
https://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/model300c1hadan.gif
腐食したケーブルのみ動く現象が起きている.
次に,橋梁全体を腐食率25%,C1端部に腐食率50%の破断時の応力ひずみ関係を与えた結果を載せる.(解析は84s程度で止まってしまった.)~
これは連鎖崩壊は全てのケーブルがある程度腐食していないと発生しないため,それに近い状態だとどうなるか見るため行った.
https://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/model300c1hadanzentaicoro25.gif
全体的に下にたわみすぎだと思う.[[この問題:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?%E7%A7%8B%E5%B1%B1%E3%81%AE%E4%BF%AE%E8%AB%96%E6%97%A5%E8%AA%8C#cc7ba732]]と似ている現象が起きてそう？~
断面2次モーメントの値を調整してみて再度解析を回してみる.~
*簡易モデルを用いた検証 [#lce2339d]
**卒論で扱った簡易モデル [#n814c429]
Q:そもそも,簡易モデルって何？
A:簡易モデルとは桁をケーブルで1本（または2本）で吊ったモデルである.~
学部生のときにこのモデルを使って,DAFを算出し,衝撃力が出ているかの確認を行なった.
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/Salomekani.png
桁は高さ10mm,幅50mm,厚さ2mmの薄肉矩形断面,ケーブルは半径4mmの円形断面である.
Salomeの設定とほとんど近い状態でMarcでも解析してみた.（プレストレスの与え方だけ違う.Marc→温度応力,Salome→ひずみ）~
また,Marcだと149.9sで解析が終わってしまうため,0.01s分データがない結果である.
ケーブル1本のみの簡易モデル（集中荷重なし）の結果
Marcの解析結果~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/CnasimarcmodelA.png
Salomeの解析結果~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/modelADEPL.png
DAF~
,時間,DAF(Marc),DAF(Salome)
,0.005s,1.97,1.95
,0.01s,1.95,1.93
,1s,1.02,1.02
,5s,1.00,1.01
1ケースしか試していないが,どちらの結果も見る限り,大きな差は出ていないため,解析ソフトによる差異はほとんどないと言えるだろう.~
この結果から,ピン接合とプレストレスをSalome用に設定をアレンジしたものもそれなりに信頼性は持てると思う.~
***DAFとtimestepの関係 [#w4d965e5]
2本モデルのとき、モデルに冗長性が出たためDAFの値が2に近い値が出なかったと考察したが...~
timestepを更に細かくしたら結果が変わるのでは？と考え、timestepを0.001に変えて再度解析を回してみた。（卒論ではパソコンの性能が足りなくて0.005までしか結果を出せなかった。）
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/kaniver2.png
解析時間を115sに変更,100sまでゆっくりと載荷し,110sで一番外側のケーブルを破断させ,115sまで解析を行った。
結果はケーブル破断時間を0.001s,0.005sのもの2種類用意した。
timestep0.001の結果~
109s〜115sの様子
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/2honmodelhadan0001.png
110s〜110.5sの拡大グラフ
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/2honkanimodelhadan0001kakudai.png
timestep0.005の結果~
109s〜115sの様子
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/2honkanimodelhadan0005.png
110s〜110.5sの拡大グラフ
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/2honkanimodelhadan0005kakudai.png
DAFはtimestep0.001sでは1.87~
timestep0.005sでは1.47となった。
以上から、2本モデルの簡易モデルでDAFの値が思ったより出なかったのは,timestepが原因だと考えられる。~
1本モデルでtimestep0.001に変えて解析回してないからはっきりとは言えないけど...timestepに大きな差（0.005と0.001）がないのに、ケーブル本数増えたらDAFの数値が変わるなら橋梁モデルももっとtimestepを細かくしたらDAFも変わってくるのかな？
**ケーブル2本の簡易モデルに部分腐食をさせてみる [#a1bd3253]
***Marc [#q11b79aa]
3/10~
大ひずみ,同条件で解析を回したが,29.03sでエラーで止まってしまった.(エラーメッセージはstep数をもっと細かくとのことだが,自動かつ細かいところで0.0001sで切ってくれているが...)
最終挙動（29.03s）と桁先端変位と腐食部C1張力は以下の通り.~
最終挙動（腐食したケーブルが伸びた後,桁が揺れる挙動が得られた.）~
https://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/I0115000n2honmodel_SIMOMIEHE.gif
C1腐食部張力,腐食していない部分も同じ値に...(Tの最大値:6172.5N)~
https://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/I01kanimodel15000N101sT_SIMOMIEHE.PNG
先端変位(-1.368mmたわんでいる)~
https://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/I01kanimodel15000N101sDEPL_SIMOMIEHE.PNG
最終挙動付近を拡大すると,桁が揺れていることもわかる.~
https://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/I01kanimodel15000N101sDEPL_SIMOMIEHE_kakudai.PNG
3/9~
微小ひずみ,集中荷重15000Nを掛けた際の解析結果だが，最終挙動（101s）と桁先端変位と腐食部C1張力は以下の通り.~
最終挙動（時間を掛けて,下にたわんだだけ.複雑な動きはしていない）~
https://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/I01kanimodel15000N101s.PNG
C1腐食部張力,腐食していない部分も同じ値に...(Tの最大値:21083.3N)~
https://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/I01kanimodel15000N101sT.PNG
先端変位(-4.48mmたわんでいる)~
https://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/I01kanimodel15000N101sDEPL.PNG
断面2次モーメントを調整するやり方（ケーブルに適用）を参考にして解析を回したが，ケーブルの降伏応力（ミーゼス応力）と塑性ひずみが0という結果に.~
設定箇所をよく見ると，特性：入力値の場合,材料挙動が線形弾性のみ（応力とひずみが常に比例の関係：SalomeでいうELAS）しか対応していないみたい.(この設定では材料非線形を考慮した挙動が見れないってこと？)~
→3/10追記：桁の降伏応力と塑性ひずみの数値が見れたため，下記写真のやり方ではこれらの数値は見れないと思われる.
https://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/tyujitudanmen.PNG
3/4 [[この問題:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?%E7%A7%8B%E5%B1%B1%E3%81%AE%E4%BF%AE%E8%AB%96%E6%97%A5%E8%AA%8C#cc7ba732]]と似ている現象が起きてそうなので、参考にしたところ、断面2次モーメントの値を調整したら改善したみたい.~
どうやらケーブルのたわみは曲げ剛性を無視し,伸び剛性と張力の影響でたわむとのこと.このモデルのケーブルの断面2次モーメントを計算したら201.061mm^4であり,無視できないほど大きいため,断面2次モーメントをIxx=Iyy=0.1にした状態で桁先端に集中荷重100Nを掛けて再度解析を回してみた.
https://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/I01kanimodel.PNG
写真は解析終了時（101sの様子,荷重が掛け終わるのは100s）の様子.解析は最後まで回り,下にたわむことはなくなった.~
この状態でSalomeに与えていた集中荷重15000Nでどういった挙動を見るか再度確かめる.
2/27~
MarcでSalomeと同じ条件で解析を回したところ,以下の結果が得られた.（途中で解析が止まってしまうが，腐食時の応力が0になっただろうと思われる挙動が得られたため，結果を載せた．）
アニメーション(主桁先端に15000N)~
https://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/15000N2honmodel2.gif
Marcだと,100N掛けても切れるような挙動が得られる.~
https://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/100n2honmodel.gif
Salomeと挙動が異なる点として，要素の違いがあると思われる．現状,Salomeでは材料非線形を考慮した梁要素の解析はPOU_D_EMやPOU_D_TGMといったマルチファイバー要素でしかできず,
しかも,マルチファイバー要素自体,圧縮が出てしまうことから,橋梁のケーブルをマルチファイバー要素で再現できない可能性がある.~
Marcでは，先輩方が使っていた橋梁モデルのケーブルに要素タイプ25を用いていたため,それを使い解析を行った.以前,静的解析での橋梁モデルをSalome:ケーブルをBAREE,Marc:ケーブルを要素タイプ25,で[[主桁中央鉛直変位を比較した:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?%E5%8D%83%E4%BB%A3%E5%B2%A1%E3%81%AE%E4%BF%AE%E8%AB%96%E6%97%A5%E8%AA%8C#gb613b11]]ところ,近い値が出ていたため,SalomeのBARREとMarcの要素タイプ25は近い性質を持っているのだと思われる.~
では,MarcでSalomeのマルチファイバー要素に近い要素タイプは何か？以前,[[簡易モデルの桁先端変位とケーブル応力を比較した際:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?%E5%8D%83%E4%BB%A3%E5%B2%A1%E3%81%AE%E5%8D%92%E8%AB%96%E6%97%A5%E8%AA%8C#k829fe9a]]，SalomeとMarcの結果がほぼ一致し，そのときの要素タイプが98だったため,SalomeのマルチファイバーとMarcの要素タイプ98が近い性質を持っているのだと思われる.~
https://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/yousotype.PNG
閉断面：構造部材の断面が完全に閉じている断面（パイプとか）~
開断面：構造部材の断面に隙間がある断面（I形鋼とか）
***Salome-Meca [#u3751526]
現状は、ケーブルに健全時の構成則と腐食時の構成則を同時に与えることはできている。ここから、腐食時のケーブルだけを切る（応力0に落とす）には,要素分割数と伸びの長さが関わる可能性あり。~
材料非線形はケーブルのみにしている。（桁に入れても,同じ材料に異なる構成則を同時に入れることは可）~
2/12~
かつて10分割した要素長と同等の長さを腐食部分、その他の長さを健全部とした。（2要素にした。）~
結果は以下の通り.
グラフはどれも応力が0に落ちた瞬間の箇所をとっている.~
時間と桁先端変位~
https://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/212_2honmodelDEPL.png
時間と腐食部張力~
https://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/212C1T2honmodel.png
応力ひずみ~
応力392MPa,ひずみ0.004で落ちるように設定.
応力がおちるグラフ~
https://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/212C1SIGM_1.png
応力が0に落ちた瞬間の応力ひずみ~
https://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/212C1SIGM_2.png
・問題点
腐食による破断後の張力を見ると、0に近い値が出ているものもあるが,時間が経つほどどんどん大きくなっている...。（最も0に近い値は0.14475N）~
応力については、完全に0にはできない状況（落ちた瞬間は1MPa程度出ている。）~
→より安定したモデル（橋梁モデル）で解析するとまた変わってくる？~
→ケーブルをcable要素に変えてみる？（けど、応力-ひずみ考慮できるの？）~
12/26~
腐食部分の要素長さ100√2mm,50√2mmの場合、解析は回ることが判明した。ただ、結果を見ようとするとフリーズする...。~
ちなみに、要素長さ100√2mmじゃ要素分割10と同じ長さであるため、解析が回る回らないの違いとして、要素長と要素分割数が関わってくる可能性が高い。~
-解析が回る（回りやすい）パターンは以下のとおり~
(1)材料非線形を考慮している要素長は長いほど安定しやすい~
(2)分割数は少ないほうが安定しやすい
12/18~
ゼミで頂いた意見まとめ~
-真ん中のみの腐食はできる？→健全-腐食-健全~
-何故,要素長によって解析できないのか？→解明する。~
次のゼミまでに、要素長と構成則の関係？を少しでも明らかにしないとな...。
12/12~
要素分割2の場合は？~
12/10での解析と同条件で解析を行った.(応力を落ちることによる破断再現に着目するため.)
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coro75youso2.gif
変位と時間のグラフ~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/youso2DEPL.png
応力と時間のグラフ~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/youso2MPa.png
この後,要素分割3,4も回したところ解析自体は回った.ただ,現状の要素分割10は解析が回らないし,部分腐食となるとできるだけ要素を細かくした1要素を腐食させることになると思われるため、順番に要素分割して確認する作業はここまででいいかと。~
要素分割2ができるなら,腐食させたい要素長を短く,健全部を長くするというモデル再現はどうだろうか？(例えると,1000mmの部材であれば,10mm腐食部,残り990mmを健全部にする)~
実際に試したところ、収束条件のエラーが発生した。~
...ということは,要素長と構成則の相性がある？（そんなことありえるのか？）
↓よく出る収束条件のエラーはこちら~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coroerror3.png
12/10~
Salome-Mecaで分割数1の場合でも解析が回らないのか調査してみた。（よくよく思えば、これを先にやるべきだった...。）~
分割数1であれば、Salome-Mecaでも解析が回ることが分かった。（精度面を一旦置いとけば、ひとまずMarcと同じことはできるでいいのかな。）~
桁の先端に集中荷重Fy=-15000N掛けた.
要素分割していないため、全体腐食を表現。(破断時間:97.109s)
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coro75youso1.gif
結果
・桁先端の変位と時間グラフ~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/youso1coro75DEPL.png
・ケーブルの応力と時間グラフ~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/yousoSIGMcoro75.png
気になる部分はいくつかあるが、応力を落とすという再現に着目すればひとまず最低限のものはできているのではないだろうか。~
11/26~
部分腐食再現で個人的にいくつか思いついたことをまとめる~
（応力を0(ほぼ0)に落として,伸びないようにするために、応力-ひずみ曲線と別の切り口がないか模索。）~
①ある時間でヤング率を下げて,ヤング率低下による腐食再現.~
→ELAS_FO(材料特性をfunctionで与えるやり方)でできそうかと思ったが、横軸時間のfunctionに対応していないみたい...。~
②応力-ひずみ関係で応力を落とす行為に問題があるのであれば、応力-時間,ひずみ-時間グラフでも与えてみる？~
→問題はそれらをどこに入れるか.TRACTIONにはfunctionを1つしか入れることができない。~
③DEFI_MATERIAUのFATIGUEを使う？~
→使い方をいまいち理解していないが,ひとまずmultifiberで使えることは確認.~
④今行っている2段階解析で、1段階目を破断寸前で止め、2段階目で破断させるやり方で進める。~
11/20~
気になるコマンドを見つけたため、ここに共有する。~
[[U4.43.01.pdf:https://code-aster.org/doc/v14/fr/man_u/u4/u4.43.01.pdf]]
[[R7.04.01.pdf:https://code-aster.org/doc/v15/en/man_r/r7/r7.04.01.pdf]]
DEFI_MATERIAU→FATIGUEというコマンドで、詳細まで調べてないが、例を出してる人を見ると疲労解析に使っているみたい。~
試しにFATIGUEに適当な設定を入れて,解析を回してみたところ、Analysisまで進んだため、multifiberに対応していないってことはなさそう？~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/fatigue1.png
11/4~
色々なパターン（2本モデル以外だと,1本モデルと塔付きの1本モデル）を試したが、どうも上手く行かない...。~
→1本モデル自体はケーブル一様に応力が出ないため,NG.塔付きだと,収束しない...。~
エラーが出る直前で解析を止めて結果を見ると,構成則通り動いているのは確認できたことから、応力を0近くまで落とす場所でエラーが出る可能性大と予想。~
試しに,応力が一定（伸び続けている箇所）になった後,応力を少し落とすような構成則に変えてみて解析を回すと、エラーが発生したため、応力を落とす挙動で問題があると考えられる.
エラーの内容~
どちらも収束絡みのエラー.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coroerror1.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coroerror2.png
とりあえず思いつく解決策として、応力ひずみ減少ではなく断面積減少で攻めてみるとか考えてみた（過去の先輩方が使っていた手法.）けど、結局応力を0に落とすことができないため、あまり関係ないか？~
10/29~
ケーブルをCABLE要素に変える案が出たため、試してみたが、TRACTIONを用いて応力ひずみを与える、Von Misesの降伏条件（VMIS_ISOT_TRAC）に対応していないことが分かった。（CABLEに対応していないというエラーが出た。）~
10/16~
死荷重のみだとケーブルの応力が塑性域までに達しないため,桁の先端に集中荷重(15000N)を掛け,応力ひずみ曲線を出して与えたケーブル構成則に従って解析できているか確認した。~
以下は,C1の腐食していない部分（青色線）と腐食部分[75%]（緑色線）の応力ひずみ曲線である。~
なお、①（腐食位置）,②の要素からデータを取得した。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/corrosionpoint.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/corrosiongraph.png
・POU_D_TGMの場合
①腐食箇所~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/75corrosionSIGMEPSI.png
②健全部分~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/75notcorrosion_SIGMEPSI.png
・POU_D_EMの場合
①腐食箇所~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/EM_75corrosion_SIGMEPSI.png
②健全部分~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/EM_75notcorrosion_SIGMEPSI.png
グラフ見たところ,健全部分と腐食部分で応力ひずみ曲線が異なる上に腐食率75%のグラフを見るかぎり与えたケーブル構成則通りに解析をしてくれてそう?~
また、張力も見たところ,健全部は19813N,腐食部は19808Nと腐食している箇所のほうが張力の値が多少ではあるが小さくなったことも確認できた。~
また,POU_D_EMに変えても解析は回ったため,プレストレスは温度に変えるほうがよさそう.(昨年あった張力が落ちてしまう問題が気になるが...→EMで起きていた問題なので,TGMの場合だとどうなるか見てみる.)
10/10~
温度に変えると,PRE_EPSIのときみたいにあまりにもおかしい結果が得られなかった.~
とりあえず,腐食率0%と75%のC1,2張力と先端変位の結果を載せる.（結果見ることが優先のため,timestepは0.1にしてある.）~
・C1張力（腐食させるところを抽出）~
腐食0%：6.71237N~
腐食75%：6.71237N
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/corrosion_T_C1.png
・C2張力~
腐食0%：12.0153N~
腐食75%：12.0153N
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/corrosion_C2_T.png
何と、腐食率関係なく張力の値が同じになった。...そんなことあるか？？~
自分のイメージだと,腐食した部分は部材の強さとして弱くなるわけだから,変位,張力ともに小さくなると思うけど...~
こうなってしまった原因だが,TRACTIONで定義した応力ひずみ関係が正しく反映されていない,もしくは反映されているが他の要因で打ち消されてしまっている可能性があると思う.(PRE_EPSIを使うと,そこで与えたひずみしか考慮しない結果になってしまうみたいな.)~
後者の場合,ケーブルに与えている外力は,温度で与えているプレストレスのため,そのあたりの設定が怪しい...?~
10/16追記~
ケーブルの応力を見ると,10MPaすら超えてなかった.原点近くの応力ひずみ曲線の傾きは腐食率関係なくほぼ同じ値であるから,結果に差異が出なかったってことかな？
10/9~
解析結果を確認したところ...あまりよろしくなさそうな結果が得られた.~
腐食0%,25%の結果を載せる.
腐食0%の変形図(25%もほぼ同じ変形図だった.)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/corrosion_kekka1.png
C1張力と時間の結果~
腐食0%の結果~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/corrosiongraph_1.png
40s付近拡大~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/corrosiongraph_1kakudai1.png
60s付近拡大~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/corrosiongraph1_kakudai2.png
腐食25%の結果~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/corrosiongraph_25.png
40s付近拡大~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/corrosiongraph_25kakudai1.png
60s付近拡大~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/corrosiongraph_25kakudai2.png
こうなってしまった原因だが,変形図を見てみるとケーブルが折れ曲がっているが,折れ曲がっている要素はPRE_EPSIを与えている部分である。~
現状,PRE_EPSIを与えた要素はPRE_EPSIのみを考慮してしまうため、それが悪さしているのかも...?~
簡易モデルだと,以前張力が一瞬下がる問題があるが,とりあえず温度をプレストレスに変えてもう一回解析を回してみる.~
10/5~
POU_D_TGMに変えたら解析は回ったのでひとまず結果確認.~
POU_D_TGM:断面の反り(曲げのこと？)のモデリングを考慮したTimoshenko梁の要素の定式化に基づく(Astercode解説書より).Multifiberなので,材料非線形も考慮可.~
10/4~
応力-ひずみ曲線の応力が一定になる箇所でエラーが起きている可能性があるため，少しだけ傾きをつけた．[（1570,ひずみ）→(1571,ひずみ)と変化できるように]~
[[弾塑性解析2021:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?Salome-Meca%E6%BC%94%E7%BF%92_%E5%BC%BE%E5%A1%91%E6%80%A7%E8%A7%A3%E6%9E%90%282021%29]]を参考にした．~
こうすると,上記のエラーは出なくなったが,今度はtimestep絡みのエラー（収束エラー）が出てしまった．~
ここをどうやって解決するか模索中...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/kaniver2.png
要素分割している外側のケーブル（C1）に部分腐食を与えて解析を回す。（内側のケーブル,C2の張力と先端変位を求める。）~
腐食の与え方はケーブルに与える応力ひずみを減少させていく方法。~
計算方法は健全状態（腐食0%）のケーブルの応力-ひずみ曲線の降伏応力,ひずみに腐食していない分を掛けて求める。~
{例:腐食25%の場合:1570×(1-0.25)=1177MPa}
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/corrosiongraph.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/corrosion_data.png
**解析結果 [#h686f34a]
Salomeで解析を進める方向に決まったので、その結果をまとめる。~
初めに解析モデルは以下の写真のものを用いた。また、ケーブルに与える応力-ひずみ関係は10/4に載せた上の写真（すぐ上の写真）の健全状態（腐食率0%）,腐食率75%を使う.~
しかし、腐食率75%に関しては第二降伏点に達した後の伸びが0.01に達する前に解析エラーが発生してしまうため、伸びを短くしてひずみが0.004に達したら破断するように調整した.（下の右側の写真,伸びは0.01以降も続いている.）~
また、ケーブルには幾何学的非線形も考慮している.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/model.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/usecorosion.png
桁断面とケーブル断面
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/ketadanmen.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/cabledanmen.png
腐食はC1ケーブル全てを腐食させた全体腐食,そして、部分腐食はC1ケーブル上部,真ん中,下部にそれぞれ与え,部分腐食①のとき,部分腐食②のとき,部分腐食③のときのそれぞれ3パターンでの合計4パターンで動的解析を用い,応答値（桁先端変位,腐食部C1張力,C2張力,C1腐食部応力ひずみ）を求めた.（部分腐食の長さは10mm。部分腐食①の例だが,何も載荷していない場合の部分腐食の長さはピンクの点の間.）~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coro10mm.png
***集中荷重の場合(0sから集中荷重掛けるパターン) [#n0e49bf6]
載荷条件：桁に死荷重；-0.017248N/mm ケーブルプレストレス;C1に-100℃,C2に-90℃ 桁の先端に集中荷重;15000N（腐食ケーブルに破断するための応力ひずみに達するようにするために載荷.）~
解析時間：110sまで解析を回し,100s掛けて載荷し,その後10s安定させる.（抽出したデータは安定状態から1s経った,101sまで）~
まず,ケーブルが全て健全状態の結果を以下にまとめた.~
挙動(スケールファクタ5倍 色がついている粒？みたいなのが、元々桁のあった位置.)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/Cnomi_kenzen.png
・101sの解析結果~
先端変位：-3.44037mm ~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/kenzen_DEPL.png
C1張力：20361.1N~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/kenzen_C1T.png
C2張力：1048.53N~
ただ、所々で張力が一定に変化しないところがある.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/kenzen_C2T.png
C1応力ひずみ：ひずみ0.00140143,応力410.622MPa（腐食率0%のため,ケーブルが破断するほどの応力ひずみは出ていない.そもそも、第1降伏点にすら達していない.）~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/kenzen_SIGMEPSI.png
次に,腐食の解析結果をまとめる~
全体腐食(97.6sで破断,解析時間101s),部分腐食①(95.40sで破断,解析終了時間95.4076s),部分腐食②(95.41sで破断,解析終了時間95.416s),部分腐食③(95.50sで破断,解析終了時間95.5038s)の解析結果.
・解析終了時の挙動（スケールファクタは全て1倍.部分腐食モデルにおいて,C1ケーブルのピンクの点と点の間が腐食部分を表している。色がついている粒？みたいなのが、元々桁のあった位置.）~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/kanimodel_corozentai.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/kanimodel_coro10mmue.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coro10mm_mannaka.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coro10mm75sita.png
・先端変位（右側の写真は,破断した瞬間の様子を拡大したもの.）
全体腐食（解析終了時の変位：-141.511mm）~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_corozentai.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_corozentaikakudai1.png
部分腐食①(解析終了時の変位：-188.472mm）
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_uecoro.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_uecorokakudai1.png
部分腐食②(解析終了時の変位：-149.295mm）
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_mannakacoro.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_mannakacorokakudai1.png
部分腐食③(解析終了時の変位：-158.937mm）
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_sitacoro.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_sitacorokakudai.png
破断する寸前の桁先端変位
全体腐食：-6.09mm
部分腐食①,②,③：-3.30mm
・C1張力（腐食部）（右側の写真は,破断した瞬間の様子を拡大したもの.）~
張力はどの解析結果でも破断後は100N程度まで落ちる.
全体腐食（最大張力：19441.5N）~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_zentaicoro.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_zentaicorokakudai.png
部分腐食①（最大張力：19437.5N）
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_uecoro.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_uecorokakudai.png
部分腐食②（最大張力：19441.9N）
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_mannakacoro.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_mannakacorokakudai.png
部分腐食③（最大張力：19440.3N）
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_sitacoro.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_sitacorokakudai.png
・C2張力（右側の写真は,破断した瞬間の様子を拡大したもの.）
全体腐食（最大張力：56682.7N）~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C2T_zentaicoro.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C2T_zentaicorokakudai.png
部分腐食①（最大張力：45423.1N）
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C2T_uecoro.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C2T_uecorokakudai.png
部分腐食②（最大張力：49039.2N）
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C2T_mannakacoro.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C2T_mannakacorokakudai.png
部分腐食③（最大張力：60822.2N）
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C2T_sitacoro.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C2T_sitacorokakudai.png
・C1腐食部応力ひずみ（左写真は応力が落ちる箇所までの応力ひずみ曲線,右写真は応力が落ちてから解析が終わる箇所を拡大したもの）~
全体腐食~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_zentaicoro.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_zentaicorokakudai.png
部分腐食①
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_uecoro.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_uecorokakudai.png
応力変化が激しいところを拡大。
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_uecorokakudai2.png
部分腐食②
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_mannakacoro.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_mannakacorokakudai.png
応力変化が激しいところを拡大。
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_mannakacorokakudai2.png
部分腐食③
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_sitacoro.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_sitacorokakudai.png
全体腐食以外のモデルでは安定状態に入る前に95s前後で解析エラーが出てしまう。恐らく、解析では計算できないほどモデルが変形してしまっているからだと思われる.~
そのため、腐食ケーブルが伸びた瞬間(応力がほぼ0になる瞬間)で解析を止め、結果を得た.~
***集中荷重の場合(安定状態を終えてから集中荷重を掛けるパターン) [#r79d5e43]
連鎖崩壊は供用した橋梁が劣化して落橋のため、安定状態にも達していないのにケーブルが破断するのは求めたいこととずれるのではと思い、集中荷重を掛けるタイミングを変えてみた~
死荷重とプレストレスを100s掛け、10s安定させた後に集中荷重15000Nを20s掛けて載荷する方法で再度解析を回した。(安定状態に達したモデルにケーブルが切れるほどの荷重が掛かった際の状況を見る.)~
データの抽出は安定状態を終える110sから行い,全体腐食(129.56sで破断;解析終了時間131s),部分腐食①(129.008sで破断;解析終了時間129.013s),部分腐食②(sで破断;解析終了時間129.017s),部分腐食③(sで破断:解析終了時間119.089s)の解析結果.
・解析終了時の挙動（スケールファクタは全て1倍.部分腐食モデルにおいて,C1ケーブルのピンクの点と点の間が腐食部分を表している。色がついている粒？みたいなのが、元々桁のあった位置.）~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/kanimodel_corozentai2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/kanimodel_coro10mmue2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coro10mm_mannaka2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coro10mm75sita2.png
・先端変位（右側の写真は,破断した瞬間の様子を拡大したもの.）
全体腐食（解析終了時の変位：-141.793mm）~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_zentaicoro2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_zentaicorokakudai2.png
部分腐食①(解析終了時の変位：-193.838mm）
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_uecoro2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_uecorokakudai2.png
部分腐食②(解析終了時の変位：-149.295mm）
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_mannakacoro2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_mannakacorokakudai2.png
部分腐食③(解析終了時の変位：-158.937mm）
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_sitacoro2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_sitacorokakudai2.png
・C1張力（腐食部）（右側の写真は,破断した瞬間の様子を拡大したもの.）
全体腐食（最大張力：19437.6N）~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_zentaicoro2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_zentaicorokakudai2.png
部分腐食①（最大張力：19344.9N）
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_uecoro2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_uecorokakudai2.png
部分腐食②（最大張力：19439.2N）
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_mannakacoro2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_mannakacorokakudai2.png
部分腐食③（最大張力：19438.7N）
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_sitacoro2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_sitacorokakudai2.png
・C2張力（右側の写真は,破断した瞬間の様子を拡大したもの.）
全体腐食（最大張力：26846.5N）~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C2T_zentaicoro2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C2T_zentaicorokakudai2.png
部分腐食①（最大張力：42880N）
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C2T_uecoro2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C2T_uecorokakudai2.png
部分腐食②（最大張力：41691.4N）
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C2T_mannakacoro2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C2T_mannakacorokakudai2.png
部分腐食③（最大張力：60829.3N）
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C2T_sitacoro2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C2T_sitacorokakudai2.png
・C1腐食部応力ひずみ（右側の写真は,破断した瞬間の様子を拡大したもの.）
全体腐食~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_zentaicoro2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_zentaicorokakudai2.png
部分腐食①
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_uecoro2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_uecorokakudai3.png
部分腐食②
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_mannakacoro2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_mannakacorokakudai3.png
部分腐食③
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_sitacoro2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_sitacorokakudai2.png
***衝撃荷重の場合 [#jfe97a58]
載荷条件：桁に死荷重；-0.017248N/mm ケーブルプレストレス;C1に-100℃,C2に-90℃ 110sでC2ケーブルを破断させる.~
解析時間：131sまで解析を回し,100s掛けて載荷し,その後10s安定させたあとにC2ケーブルを破断させ(110sで破断),そのときに生じる衝撃力による応答値を求める.~
まず、衝撃荷重のみの結果を載せる.全体腐食と部分腐食③（ケーブル下側部分の腐食）の挙動とC1腐食部の応力ひずみを見た.
・挙動
（全体腐食,スケールファクタ10倍 ケーブルプレストレスにより桁が持ち上がった動き.）~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/Cnasikanimodel_impactzentaicoro.png
（部分腐食③,スケールファクタ100倍 C2の衝撃力が伝わったのか桁の真ん中が下がっている.腐食部分と健全部の応力ひずみ関係が異なるため,このような挙動になったと思われる.）~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/Cnasikanimodel_impactsitacoro.png
・C1腐食部応力ひずみ
全体腐食の場合~
プレストレスによってケーブルが桁を引っ張っているため、ケーブルが縮むような（圧縮ひずみ？）挙動が得られるのかと思われる.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_impactzentaicoro.png
部分腐食③の場合~
全体的に見ても、応力はほとんど生じていない.健全部と腐食部で応力ひずみが異なるため、全体腐食と異なる挙動になったと思われる.~
また、健全部の応力ひずみだが、衝撃力がかかった瞬間,応力ひずみともに正の値が出るようになる.(衝撃がかかるまでは圧縮と引張を行き来している挙動.)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_impactsitacoro.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_impactsitacorokenzen.png
結果を見ると、腐食部分が破断することはなかった.(C1のケーブル応力が0.18MPaと低いため、やはり集中荷重を掛けないとケーブルの応力に大きな値が得られない.)
***衝撃荷重の場合（集中荷重も加えたパターン） [#ff632edf]
衝撃荷重だと腐食部に伸びが生じなかったことが確認できたため、やり方を変えてみた.~
解析時間は130sで、まず100sプレストレスと死荷重を掛け.10s安定させる.その後に集中荷重(15000N)を掛けて,その最中にケーブルが切れる状態(120sで破断.)を想定し、130sまで集中荷重を載荷し、解析を回した.（衝撃荷重+集中荷重）~
結果は全体腐食(128.4sで破断;解析終了時間128.405s),部分腐食①(128.279sで破断;解析終了時間128.287s),部分腐食②(128.347sで破断;解析終了時間128.354s),部分腐食③(128.39sで破断:解析終了時間128.393s)の4ケースをまとめた.~
また、どのケースも130sまで収束エラーで解析が最後まで回ることはなかったため、C1ケーブルが破断した瞬間までの記録を取っている.
・解析終了時の挙動（スケールファクタは全て1倍.部分腐食モデルにおいて,C1ケーブルのピンクの点と点の間が腐食部分を表している。色がついている粒？みたいなのが、元々桁のあった位置.）~
[[集中荷重のみ:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?%E5%8D%83%E4%BB%A3%E5%B2%A1%E3%81%AE%E4%BF%AE%E8%AB%96%E6%97%A5%E8%AA%8C#r79d5e43]]と比較して、どの解析パターンもC2のケーブル破断による衝撃荷重が加わったことで腐食部分の要素が大きく伸びていることが分かる.
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/impactkanimodel_zentaicoro.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/impactkanimodel_uecoro.png~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/impactkanimodel_mannakacoro.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/impactkanimodel_sitacoro.png
・先端変位（右側の写真は,破断した瞬間の様子を拡大したもの.）
全体腐食(最終変位：-415.62mm)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_zentaicoroimpact.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_zentaicoroimpactkakudai.png
部分腐食①(最終変位：-879.807mm)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_uecoroimpact.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_uecoroimpactkakudai.png
部分腐食②(最終変位：-703.065mm)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_mannakacoroimpact.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_mannakacoroimpactkakudai.png
部分腐食③(最終変位：-354.928mm)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_sitacoroimpact.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_sitacoroimpactkakudai.png
・C1張力（腐食部）（右側の写真は,破断した瞬間の様子を拡大したもの.）
全体腐食（最大張力：19441.5N）~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_zentaicoroimpact.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_zentaicoroimpactkakudai.png
部分腐食①（最大張力：19438.8N）~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_uecoroimpact.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_uecoroimpactkakudai.png
部分腐食②（最大張力：19351.1N）~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_mannakacoroimpact.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_mannakacoroimpactkakudai.png
部分腐食③（最大張力：19441.3N）~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_sitacoroimpact.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_sitacoroimpactkakudai.png
・C1腐食部応力ひずみ（グラフは破断後の応力ひずみ曲線を表している.）
全体腐食~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_zentaicoroimpact.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_zentaicoroimpactkakudai.png
部分腐食①~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_uecoroimpact.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_uecoroimpactkakudai.png
部分腐食②~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_mannakacoroimpact.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_mannakacoroimpactkakudai.png
部分腐食③~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_sitacoroimpact.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_sitacoroimpactkakudai.png
○衝撃荷重が加わった際の腐食部の応力ひずみ関係の特徴として、120sでの衝撃荷重が加わった後,応力ひずみともに正の値に変わる変化が見られた.
***まとめ [#x3c721bc]
集中荷重のみの場合、部分腐食で解析を最後まで回そうとすると収束エラーが出てしまう。破断時の挙動と応答値を得ることはできたが、全体腐食との解析時間が異なるため、比較するのは難しいかも。~
ただ、部分腐食のグラフを見ると、挙動はそれぞれ似ている一方、解析終了時の値が異なるので、腐食の位置によって応答値の変化がある可能性も？~
全体腐食と部分腐食での応答値の違いを見る場合,衝撃力と集中荷重を掛けたパターンが解析終了時間もほぼ同じのため、比較しやすい。~
→結果を見る限り、似たような傾向の結果が得られているのではないか？ただ、桁の先端変位が腐食ケースごとに値が変わるのかは不思議。挙動を見たところ、腐食部分の伸びの長さがそれぞれ異なっているので、腐食の長さや位置によって応答値の変動はあるのかもしれない。~
**簡易モデルver.2 [#nb955af1]
ばねを用いてピン接合を行うやり方で斜張橋モデルを作成すると、主塔とケーブルを結ぶばねの節点に自由度絡みのエラーが発生する。~
ケーブルをトラス要素に変えたら解析は回ったが、トラス要素は節点共有した際にピン接合になるため、要素分割する際は橋梁ケーブルをトラス要素で再現できない。~
ケーブルのサグを再現するには要素分割を行う必要がある（ケーブルが切れたあとの様子・変化が見たい）ため、まずは簡単なモデルを作成して解析が回るか確認する。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/zentaikanimodelver2.png
今回用いるモデルは上の通りである。1.5mの塔に、塔の高さ0.5mから桁を1mに伸ばし、ケーブルで吊ったモデル。~
このモデルを4つのステップに分けて解析を行う。(解析自体はシンプルに静的解析を行い、変位を導出している。)
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/ketadanmen.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/cabledanmen.png
桁と塔は高さ10mm,幅50mm,厚さ2mmの薄肉矩形断面、ケーブルは直径8mmの円形断面を適用した。
STEP1：ケーブルなしモデル~
塔と桁のみの状態で、桁の先端に100Nの集中荷重を与えたモデル。塔の先端は拘束なし。まず、このモデルの解析が回るか確認した。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/step1.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/modelstep1factor5.png
色付きが変形図（factor5倍）~
桁の先端はy方向に-70.38mm、塔の先端はx方向に19.22mmたわんだ。
桁の先端に荷重を掛けていることから、塔と桁が倒れるような挙動が得られると思われるため、解析は上手くいっているだろう。
STEP2：ケーブルつけたモデル
塔側のケーブルの結合は剛結、桁との結合はバネ(0.01mm)を用いたピン接合。ケーブルの要素分割数は20で、プレストレスを温度で-100℃与えた。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/step2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/modelstep2factor20.png
色付きが変形図（factor20倍）~
桁の先端はy方向に1.059mm持ち上がった。~
温度によるプレストレスが入っていることも確認。（設定方法は[[こちら:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?Salome-Meca+%E3%83%93%E3%83%BC%E3%83%A0%E8%A6%81%E7%B4%A0#vac22495]]に記載。）
STEP3：塔とケーブルの接合を変位固定で~
塔側のケーブルの結合は回転自由、それ以外の条件はSTEP2と同様~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/step3.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/modelstep3factor10.png
色付きが変形図（factor10倍）~
桁の先端はy方向に2.313mm持ち上がった。~
STEP4：塔とケーブルの接合にバネを用いる。~
塔側とのケーブルの結合もバネによるピン接合に変更~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/step4.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/modelstep4factor10.png
色付きが変形図（factor10倍）~
桁の先端はy方向に1.249mm持ち上がった。~
*つりあい解析 [#t1c73e2f]
つりあい解析とは、事前に解析した構造解析の結果を、次に解析したいモデルにその結果を入れて解析を回すことを指す。~
ちなみに,似たような解析で,連成解析と呼ばれるものが存在する。（温度解析→構造解析とか異なる解析を繋げて行うこと。）~
学部時代に取り扱った解析だと...①：ゆっくり力を与えて安定させるまでの結果を出して、②：次にケーブルに掛けた張力を徐々に落としていく解析にその結果を読み取って解析を回すことになる。（二段階解析というイメージかな？）~
やり方は[[Salomeメモ:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?Salome%E3%83%A1%E3%83%A2#gac275b5]]と[[Salome-Mecaビーム要素:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?Salome-Meca+%E3%83%93%E3%83%BC%E3%83%A0%E8%A6%81%E7%B4%A0#gac275b5]]に記載した。(書かれている内容はどちらも同じ)
*発表等での質疑応答メモ [#n1155511]
研究室での中間発表や学会など、質問されたことをここにまとめ、再度確認したものをここに載せる。~
大阪見学で頂いたご意見~
・モデルについて、力の伝播のことを考えると桁だけでなく塔もつけた簡易的なモデルで解析すべきでは？~
→それをやるなら、橋梁モデルで解析したほうがよい。そのため、これから橋梁モデルでの解析を行う方針となった。
・今の設計法では、吊り形式橋梁は連鎖崩壊しにくいと意見を頂いたが、設計法とは何を指していたのか。~
→部分係数法が変わっているらしい。この研究では、設計法が変わる前（昔の橋梁）を対象としているため、これから作られる橋梁に対して何か言うわけでなく、既存の（設計法が変わる前）橋梁の連鎖崩壊が起きないように考える、研究となる。
春休みの研究報告~
簡易モデル2本でケーブル全て健全部,全体腐食,部分腐食での数値解析を行う.(取るものは応力-ひずみ,張力(健全部と腐食部),桁先端変位)~
→予想では全て同じ結果になるんじゃないかな？（同じ結果になれば、腐食モデルを全体と部分で異なるのかについて証明できる.）
荷重の掛け方は、集中荷重と等分布荷重とC2のケーブル破断荷重で.~
11/28 理工学デザイン~
・部分腐食させて連鎖崩壊を起こす内容は理解したが、結局それをどう活かすのか。（研究の最終目的は？）~
部分腐食の再現をすることで、ケーブル破断の影響を見ることが目標。（最終的には、衝撃力をあるパラメータとして数値化して入れることで公式的に導出できるようになればいいかな。）
・現状,応力-ひずみ減少で腐食再現ができていない状況だが、他の方法で検討していないのか。
断面積減少という方法はある。ただし、この方法ではただ要素が伸びて終わるだけなので、連鎖崩壊に繋げる意味では適切ではないかも。~
また、短期的（南方澳大橋みたいにトラック通ってから切れるみたいに）な破断をするという意味では「ひずみがある点に達したら切れる」が適切だろうと考えている。~
（発表時にはヤング率を時間で減少させて落とす方法もあるかもと言ったが、あくまでそれは長期的な話なので自分の研究内容とずれる。）~
・連鎖崩壊するには、他のケーブルも腐食してないと起きないのでは？（1本腐食したところで、他のケーブルが健全であれば問題ないよね？）~
そのとおり。この研究は連鎖崩壊のメソッドを理解することにもつながる研究なので、どの腐食部分が破断すると次にどのケーブルが切れるかという法則性？をある程度理解するために行っている。
まとめ~
指導教員からアドバイスを頂いた上でこの項目を編集した。理工学デザイン→振り返りをしたことで、改めて自分のやるべきことが何か理解した。~
この経験を踏まえ、これからの研究を引き続き進めていく。
11/13 ポスター形式中間発表~
その時の資料は,k2→2024→chiyookaに入っている.11/13にアップロードしたので,アップロード日付を見ればすぐ分かると思う.~
・要素分割とは？~
A:ある要素を任意（自動的に）に分割すること。~
・ケーブルの部分腐食の先行研究・文献はないの？（海外とか）~
また、落橋したケースの逆解析はやられてないの？~
※逆解析：すでに分かっている結果をもとに、未知の解析条件を推定すること。~
A:先行研究等に関しては調査中.（調べた範囲では行われていないが,見落としがあるかも。)~
梁要素で部分腐食の論文は調べた範囲では見つかってないが、海外の論文は見落としていたため、発表時は曖昧な答え方をしてしまった...。~
逆解析は行われている。
・全体腐食するのではなく、ポイントごと（要素分割ごと？）で腐食を考慮している？~
A:まだ部分腐食の再現ができていないため、それは考慮していない。~
・応力-ひずみ関係を参考にした実験について詳しく.~
A:外部の方がやられた研究を参照したため、ここでの記載は避けます。~
（研究ノートに残す）
・腐食が起きる原因は？~
A:金属の表面に水と酸素があるとき、化学反応（酸化反応）によって浸食されることで起きる。~
詳しくは[[こちら:https://www.nipponsteel.com/company/publications/quarterly-nssmc/pdf/2018_23_12_15.pdf]]
ちなみに,海側で腐食が起きやすいのは海からの塩分が酸化反応を手助けするためである。~
・何年ぐらいで落橋に至るのか？（完成してから何も対策もしないで）~
A:資料によってバラバラだが、約50年前後。勿論、対策をすればもっと伸びる。~
・端部で腐食させた理由は？~
A:ケーブルの端部は水が貯まりやすい場所でもあるため、腐食が起こりやすい。そのため、まずは腐食が起こりやすい場所で部分腐食の再現をしてみようと思い、端部にした。~
・応力-ひずみ減少は健全,腐食の位置を変えても同じように機能する？（発表時,ケーブル端部にのみ腐食を与えていた。）~
A:要素ごとに応力-ひずみ関係を変えただけなので、機能する。
・ケーブルの腐食だけを検討している理由は？桁とかも考えなくていいのか？~
A:連鎖崩壊の解析に繋げる研究でもあり、その現象はケーブルに着目している。そして、連鎖崩壊はケーブルの腐食の有無が重要であることから、ケーブルのみを見ている。研究が進めば、桁とかを考慮しても面白いかもしれない。~
発表時,どうしてケーブルの腐食を研究しているのかという意図で回答したが...ずれた回答していたかも。
・①ケーブルを切るのではなく、ケーブルを伸びきった状態にするのはどうか？仕事をしていないものにもっていけばそれは同じ意味なのでは？~
②連鎖崩壊に繋げる際、伸びきったケーブルの力を取って、他のケーブルに当てはめればケーブルの荷重分配もできるのでは？~
A:①伸び切ってケーブルそのものの役割を果たさないという考えはありかも。その方法も1つの手法として検討するのもあり？~
②手動ではなく,自動で違うケーブルに力が伝わるようにしたい。
・応力-ひずみ関係の破断（応力0に落とす）させるところまでの直線は何を表している？~
A:ケーブルの伸びを表している。
反省点~
たくさん質問してくれたにも関わらず、自分でも納得のいく回答をしたのが少ないと感じた。~
特に、1発目で想定外の質問が来て、そこから終始焦っていたため、もっと肩の力を抜いて回答できるとよかったかも。~
そして、今回の発表を通じて、自分の研究を理解しているように思っていたが、まだ詰めの甘いところがあると思った。~
特に、自分の研究と他の研究の繋がりの理解が足りてないと感じたため、ポスター発表までそのあたりをしっかり見直そう。~
最後になりますが、質問してくださった皆さん、ありがとうございました。
*参考資料等 [#a2f72b97]
[[衝撃メモ:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?%E8%A1%9D%E6%92%83%E3%83%A1%E3%83%A2]]（DAFとは何なのか,改めて勉強し直す用）
[[code_asterの動画:https://www.youtube.com/@DigvijayPatankar/playlists]] （音声は英語だが、内容は振動解析に関わってくるもの。）
[[code-asterの説明:https://code-aster.org/doc/v15/en/index.php?man=commande]] （英語で書かれているので翻訳する必要あり）
[[code-aster質問ページ:https://forum.code-aster.org/public/]]
[[beam要素の弾塑性解析マニュアル:https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/tamura/beam_nonliner.pdf]]
[[Plasticity tutorial rev.1.1Code Aster® STA10.3:http://www.caelinux.org/wiki/images/3/39/Claws_Plasticity.pdf]]（英語だけど、弾塑性絡みの設定について書かれている。）
[[CAE用語辞典:https://www.cybernet.co.jp/ansys/learning/glossary/]](解析の種類などの意味を調べるのに使う)
[[Marc Mentat2018メモ:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?Marc/Mentat2018%E3%83%A1%E3%83%A2]]
[[例題で学ぶMarc有限要素法解析入門:https://www.lib.akita-u.ac.jp/opac/volume/815629?current=1&q=Marc%E6%9C%89%E9%99%90%E8%A6%81%E7%B4%A0&total=2&trans_url=%2Fopac%2Fsearch%3Fcount%3D10%26defaultpage%3D1%26defaulttarget%3Dlocal%26order%3Drecommended_d%26q%3DMarc%25E6%259C%2589%25E9%2599%2590%25E8%25A6%2581%25E7%25B4%25A0%26searchmode%3Dnormal]] (2013年に出版されたものなので、少し違いがあるかもしれない。なお、閲覧にはAUアカウントとPW必要)
[[Salome-Mecaのチュートリアル:https://code-aster.org/V2/spip.php?rubrique67]] (英語だが、Salome-Meca全体の話が色々載っている。例えば、作成した構造物に適した要素(3D,BAR...etc)とか)
[[Pavaris wiki:https://www-paraview-org.translate.goog/Wiki/ParaView?_x_tr_sl=en&_x_tr_tl=ja&_x_tr_hl=ja&_x_tr_pto=sc#Analyzing_Data]]
[[Salome-meca beginning with code_aster:https://archives.framabook.org/docs/Code_Aster/beginning_with_code_aster-jp_aubry-20190129.pdf]] (Salome-mecaの質問するところ[掲示板]でよくこれを見てくれと答える人が多いので、載せてみた。英語で書かれている。)
[[Salome-Meca 2019を用いたトルク解析について:https://qiita.com/Jun_Tatsuno/items/3bc1d5d30e9055fbe388]] 0D要素の作成方法が書かれている。~
[[Midas解析編:https://in.midasit.co.jp/wp-content/uploads/2020/07/constrcution-83-civil-book.pdf]] 今年からこの解析ソフトを使い始めている方がいるので,自分も軽く頭に入れておきたい.
*解析ソフトの使い方 [#q171dd69]
Salome-Mecaビーム要素については[[こちらから:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?Salome-Meca+%E3%83%93%E3%83%BC%E3%83%A0%E8%A6%81%E7%B4%A0]]
Marcに関することは[[こちらから:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?cmd=read&page=Marc_Mentat%E3%83%A1%E3%83%A2&word=Marc]]
**midasについて [#vf9f6c8f]
修論落ち着いたら,midasでmodel300でも作ってみようかな？
***解析の種類 [#r525f7e5]
・静的増分解析~
Salomeでいう,MECA_STATIC?
・固有値解析（モード解析）~
Salomeでいう,CALC_MODE.
.応答スペクトル解析~
・時刻歴応答解析~
→そこから,線形時刻歴応答解析,動的フレーム非線形解析,境界非線形解析もできるとのこと.
Salomeでいう,DYNA_NON_LINEと同じ解析ができそう.
・P-デルタ解析~
p-デルタ効果と呼ばれる,「水平荷重と鉛直荷重の相互効果」を考慮した解析とのこと.
・座屈解析
・幾何学的非線形解析~
Salomeだと,STAT_NON_LINE（変形をGROT_GDEPにした場合）.
・静的材料非線形解析~
Salomeだと,STAT_NON_LINE.
弾塑性解析はこれを使うのかな？
・施工段階解析~
Salomeだと聞いたことがない解析.~
（定義）施工途中の荷重や部材応力などの状態変化を考慮して、遷移する事象ステップごとにモデルの形状や状態を一部変更しながら応力や変形を引き継いでいく段階的な計算方法.
*橋梁写真 [#dd55c1f6]
自分で撮った橋梁の写真をここに載せます.（不定期更新）~
(1)瀬戸大橋 (写真手前側が北備讃瀬戸大橋,奥側が南備讃瀬戸大橋)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/setobridge.png
撮影場所：香川県坂出市 与島PA~
高速道路の休憩所から撮影。四国に行くときは毎回瀬戸大橋を渡るため,ここのPAに寄るたびに撮影してます(^^)~
瀬戸大橋は道路と電車が走れる構造（鉄道道路併用橋）になっており、この構造での橋梁の長さとしては世界一だそうです。~
(2)明石海峡大橋~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/akashibridge.png
撮影場所：兵庫県淡路市 淡路SA~
日本で最も長い吊り橋です。（世界で見ると2番目です。世界1位は、確かトルコの1915チャナッカレ橋だったかな。）~
夜になるとライトアップするので、高速のサービスエリアから神戸の夜景と一緒に見るのが個人的にオススメ。（ライトアップの種類は、季節等によって変わってたはず...）~
また、ここのサービスエリアには観覧車があるので、それに乗って景色をみるのもありかも。~
(3)由利橋~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/yuribridge.png
撮影場所：秋田県由利本庄市~
秋田県にある斜張橋です。~
後輩たちがこの橋梁を研究対象に活動を頑張っているのを見て、自分も見に行ってきました！
朝に行って来たけど、交通量もそれなりにあって、散歩で使っている人もちらほら見かけて、周辺に住む人達の生活にとって大事な橋だなと実感しました。~
ちなみに...昔はトラス橋だったらしく、架替え？にあたって今の斜張橋になったらしいです。
(4)東京レインボーブリッジ~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/rainbowbridge.png
撮影場所：東京都港区
東京にある有名な吊橋です。~
東京で学会があったので、ついでに行ってきました！~
歩いて渡ることもできるみたいですが、今回は時間がなかったため断念...。~
次来たときは、歩いて渡りたいな(^^)
*日記 [#m48ef041]
12/4(水)~
今年も残り1か月,月日が経つのが早いと感じます。~
研究成果をそろそろまとめ始めないと,修論までに間に合わなくなりそうだし,計画的に動かないとな...。~
12月は,同期との旅行や中間発表とイベントがたくさんあるし,やるときはやる,楽しむときは楽しまないと！
11/9(日)~
気温も下がり,冬が近づいてきたと近頃感じます。~
後期も色々とやることが多く、今年度は学生生活の中でも最も充実な1年になりそうです。~
ただ,最近は自分自身のことで悩むことも増えてきました。（あんま考えすぎもいけないと思うけど...）~
何か人として足りていない部分があるんだろうな...ほんと難しいな
10/29(水)~
中間発表も終わり,3年生も研究室に配属され,10月もいよいよ終わりに近づいてきました.~
今年も残り2ヶ月で終わりで,1年あっという間だなとも感じます^^;~
学生生活も残り約半年で遂に終わりを迎えますし,学生時代にできることは今のうちにやっておかないとな...！
9/3(水)~
先日,秋田に戻ってきました.~
夏休み(8月)は色々な場所に出かけたり,社会人になった友達と会ったりと充実した夏休みを過ごせました.~
今日から研究室生活が再び始まりますが,残り半年頑張っていきたいです.
8/7(木)~
前期がようやく終わりました.~
今年の前期は,本当に色々あったな...~
何なら,学生時代の春から夏にかけてにおいて，一番忙しかったかもしれない(-_-;)~
帰ったら,まずは春に行けなかった自動車の本免試験を受けに行かないとな.まだまだ気は抜けない！~
終わったら,前期で遊べなかった分,息抜きしないとな！（勿論，やることはやった上で）
7/17(木)~
色々やることがあって大忙しだけど...学ぶことも多い日々を過ごしています.~
何というか...物事を表面だけ理解していて,深くまでは全く理解できていないなと感じるな.~
ほんと,まだまだ甘い部分があるなって痛感させられます。~
もっとしっかり勉強しないとな。
6/30(月)~
明日から7月,前期も残り1か月少し。早いな...研究が思っていた以上に進んでなくてだいぶ焦っています💦~
ただ、梁要素についてのSalomeに関する謎が1つずつ解決できてきて、橋梁解析がようやくできそうってところまで来ているから,もうすぐ色々なことがやれるようになるはず！~
研究以外にもやることが多く大変だけど...1つずつ丁寧に片づけていかないとな。~
しかし...修士2年になって3ヵ月、卒業された先輩方みたいに後輩を引っ張っていける人間に、少しは近づけているのだろうか。（深く考えすぎもいけないと思うけど...）~
ほんと、しっかりしないとな。
5/24(土)~
学会終わって1周間弱,昨日この前と2日間,息抜きをしました。~
学会終わって研究やるかと思っていましたが,あまり取り組むことができませんでした...流石に少し疲れていたか^^;~
おかげで,気分もリフレッシュできたので,ここから気持ちを切り替えて取り組みたいと思います。~
まずは,夏休み前までにある程度の成果を出さないとな！頑張ろう!!
5/21(水)~
英語の学会が終わり、普段どおりの生活が戻ってきました。~
学会に参加した感想ですが,英語力が全く足りないなと実感しました。~
発表以外での英語の会話はあまり喋れなかったし,質問も何言ってたか分からなかったしな...。~
将来,必要になる場面がまたあると思うので、しっかり勉強しなきゃなと改めて思いました。（昔,中学の英語の先生から言われた「英語からは絶対に逃げられない」は嘘じゃなかったな^^;）
4/27(日)~
かなり久しぶりの更新です^^;~
修士2年になり、そして最後の学生生活が始まりました。~
冬から春休みにかけて、学会の概要作成や就職活動などもあり、そして最近は英語の学会に向けての準備や大阪での工場見学もあり、特に今月は色々とイベントが多い月でした。~
今年度は修論もあって、昨年度より忙しい1年になりそうですが、楽しい1年にしたいと思います(^^).
...博士の先輩が卒業されたから、今年度は修士2年が最高学年なんだよな(・・;)~
しっかりしないとな。
12/10(火)~
自分の関わりのある今年の大きめのイベントがすべて終わりました。~
帰省までの残りの時間を研究活動に費やそうと考えています。研究を進めるスピード遅かったからな...~
冬休みの予定も立てたし、それを楽しみにしながら2024年残りの時間を楽しく過ごしたいと思います!~
10/21(月)~
また久々の更新です。~
10月も後半、この時期になるともうすぐ1年が終わると感じるようになります。~
さて,先週は3年生が研究室に配属され、顔合わせを行なったそうです。私は授業があったため、まだ顔を1度も見れていませんが、今度の飲み会でどんな子が来たのか見れるのが楽しみです。~
しかし...自分が研究室に配属されて、もう2年か(・・;) この2年でちゃんと成長できているのだろうか。~
ここでの生活もあと約1年半、悔いのないように過ごしていきたいな。
7/20（土）~
久々の更新です。~
前期も残り3週間ほどで終わります。早いな...。~
来週はある授業のスライド発表に、再来週はある授業のテストがあったりで大変ですが、忙しい日々も夏休みに入れば一旦落ち着くと思います(^^;~
ただ、授業のほうもある程度落ち着いたら研究のほうに力を入れなきゃなとも思っています。昨年度より研究の進捗ペース落ちてるからな...。~
夏休みの楽しみも頭の片隅に入れながら、残りの前期も頑張って乗り切りたいと思います！
5/15(水)~
今日はいよいよ4年生の研究テーマが決まりました🎉~
今年の研究テーマで自分が何かしら大きく関わることはなさそうですが、自分も卒論の際に先輩にたくさん助けていただいたので、もし何か聞かれたら少しでも力になれるようにはなりたいなと思います。~
自分の研究も最近はスローペースでしたが、院での生活も慣れてきましたし、ペースアップしていきます！！
5/9（木）~
GWも明けて、ようやく前期も本番かなと感じます。（そう思うのは私だけですかね💦）~
昨日はM2以上の院生と先生方が不在だったため、自分が前に立ってゼミを進行しました。~
4年生がこれから行う研究内容について分かる範囲で簡単に説明することが主な内容でしたが、あの説明で伝わったのだろうか...(・・;)~
毎回思うことだけど、日本語って難しいよな...。大学に入学した頃と比べると、間違いなくコミュニケーション能力は上がっているけど、まだまだ改善するところはあるよな💦~
まぁ...それは置いといて、最近は人前に立って話すことが増えたなと感じています。~
今年は、ゼミを除くと卒論,東北支部,謝恩式と大きめな行事で人前で話していますし...中高の時の自分がこれを知ったら驚くだろうな(笑)~
逆に、人に分かりやすく伝えるにはどうすべきかと考えることも増えましたが、そのおかげで少しずつだけどコミュニケーション能力がついてきているのかなとも思います。~
院の授業でも議論・発表がメインになってくるので、この1年でかなり鍛えられるんじゃないかな？もっと頑張らないとな(^^)
4/24(水)~
もうすぐGWが近づいているのと同時に4月がもうすぐ終わってしまうことに驚きを感じています...。~
のんびり過ごしてたらあっという間に夏ですね(・・;)~
それと、本日も英語ゼミに4年生と混じって参加しましたが、自分の研究を英語で説明するって難しいなと思いました。~
日本語でも中々難しい言葉を使っているので尚更ですね💦~
間違いなく今の4年生に伝わってなかったので、ちょっと考えないといけないな...。
4/18(木)~
研究指導計画書も書き、今年度をどう過ごすかある程度考えることができました。~
前期は結果を出すことに集中して、後期からは論文を書けるようにできると理想ですが…果たしてどうなるか💦~
けど…後輩たちの頑張ってる姿を見ると先輩がだらしないのはいけないなとも思います。~
すでに夜遅くまで残って研究してる人が何人かいるし…(゜o゜;~
うかうかしてたらあっという間に追いつかれそうです(^^;（競争してる訳ではないですが…）~
もっと頑張らないとな！
4/16(火)~
新学期が始まって1週間程経ちました。明日から本格的にゼミが始まるので,そろそろペースを上げて頑張らないとなと思っています。（一応,明日話すことはあるので,その心配はしなくてよいですが...）~
さて...明日は新4年生の春課題の発表日です。昨年は自分も発表しましたが、あれから1年経ったと思うと色々思うことはありますね(^^;~
めちゃくちゃ緊張してたし、質問の受け答えに対しても意味わからんこと言ってたよな…(-""-;)
昨年の自分を思い出しながら,明日の発表は聞こうかなと思います😊
4/12(金)~
第1回目の授業がほとんど終わり、どれも面白いと思える授業ばかりでした。~
後は、自分が授業を通じてどれだけ自分のものにできるかが大事なので、しっかり勉強していきたいと思います。~
研究室絡みでは、新4年生が来週に春課題の発表があるため、研究室に多くの人が集まっています。~
どんな発表をするのか楽しみです！~
それと...もっと新4年生とも交流しないとな。外国語文献の講義は院生も出れるから、まずはそこで色々お話してみるか！英語だけど...笑
4/10(水)~
遅くなりましたが、ついに大学院生としての生活が始まりました。~
授業に研究と、ひょっとしたら昨年度より忙しい日々になりそうですが、1日1日を大事に過ごしていきたいなと思います.~
今年度はどんな1年になるかとても楽しみです.~
現在の閲覧人数: &online;~
総閲覧数:&counter;~
終了行:
#contents
[[自分の卒論日誌:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?%E5%8D%83%E4%BB%A3%E5%B2%A1%E3%81%AE%E5%8D%92%E8%AB%96%E6%97%A5%E8%AA%8C]]はこちら
研究タイトル：部分的に腐食したケーブル系橋梁の連鎖崩壊
*研究目的 [#ya05820f]
修士2年~
編集中
修士1年
4年次ではSalome-Mecaを用い,簡易モデルで衝撃力の再現が妥当か検討を行なった.ただ,モデルを作成するにあたり,ピン接合の再現や幾何学的非線形絡み（まだ解決してないところもあるが...）の問題があり,それらをできるようにした.~
修士では,これまでMarc_Mentatを用いて橋モデルを作成して解析を行なっていたものをSalome-Mecaでもできるようにする.そのためにまずはSalome-Mecaで橋モデルを作成して,解析できるところまで持っていく.（これが第1ステップ!）~
そして,ケーブルに腐食を考慮した斜張橋に動的応答解析を行い,衝撃力の大きさを検討していく.（簡単に言うと,昨年,簡易モデルでやったことを橋モデルで行う）~
斜張橋モデルに関しては,ケーブルの張り方や本数,橋梁のスパン長や減衰率等を色々変えた上で解析を行う.
→詳細記述~
ケーブルを有する橋梁のケーブルが腐食することで橋が落ちてしまう事故が近年各地で起きています.（例えば...台湾の南方澳橋とか.卒業された先輩がその橋を対象とした研究を行なっていました.）~
その中で,1本目のケーブルが切れることで,次々とケーブルが切れてやがて落橋に至る連鎖崩壊をターゲットに研究を行なっています.~
この現象は橋梁の落橋などを引き起こす重大な事故に繋がるため,そのメカニズムの研究も進められている現状です.~
本研究室ではその解析手法の提案が行われてきたが、ケーブルが切れた際の衝撃力の大きさや他のケーブルに掛かる荷重の再分配率が明確ではありません...。~
卒論では,簡易モデルと呼ばれる簡単なモデルで衝撃力の検討を行ないました.そこで,修論では斜張橋モデルを使って手法の妥当性を検討する目的で,健全と腐食ケーブルを有する斜張橋に対してケーブルが切れた際の衝撃を衝撃力として数値解析的に与え,その大きさを検討していきます.~
モデル自体もケーブルの張り方やスパン長等の性質・値を変化させていくことでどう変わっていくのかも見ていきたいなと考えています.~
後は,Salomeのビーム要素絡みの問題も色々解決していきたい。
*今後の予定 [#bd1b97fd]
Salome-Mecaの梁要素絡みの問題は色々あるため,卒業するまでは後輩達とともに1つずつ解決していきたい。~
12月~
12/23 中間発表第2回目~
12/26 大学生活 最後の帰省
夏休み~
〇研究~
・ビーム要素におけるモーダル解析の理解~
・橋梁モデルでC1に着目した解析(C1健全時が切れたときの応答値,C1部分腐食時の切れたときの応答値)~
→そのときのC13張力は？
〇その他~
7月~
・橋梁モデルで腐食を考慮した場合の終局強度を求める。~
・交流会の準備~
・今取り組んでいるe-learningの英語の発表準備~
やることが多いけど...何とか乗り切ろう.~
6月~
・橋梁モデルで部分腐食を考慮した解析~
簡易モデルは卒業。~
一応,Salomeでやることになっているが,どうしてもできないときのためにMarc(midas)で取り組むことも視野に入れて行動したほうがいいかも。~
自分のやるべきことをやった上で取り組むなら,誰も何も言わないだろう。
コメント：時間は限られているけど,焦らず1つずつ進めていこう。
4月,5月~
・部分腐食モデル：分布荷重での解析~
→終了後,集中・分布・衝撃荷重の応答値変化をまとめたい。
・ブリッジコンペのモデル作り~
解析はSalome(ダメならmidas),モデル作成はfreecadで行う。（作ったモデルはSalome,midasに持っていけるらしい。Salomeは確認済み）
コメント：いよいよ最終学年！早速やることが多いけど...これまでやってきたことを修論で発表できるように頑張ろう！
*Model300をMarcと解析結果比較[#n498dba5]
**モデル作成 [#je113eea]
Salome-Mecaで橋モデル（model300）を作成してみる.~
以下の写真らはMarcで卒業された先輩が作成した橋モデルである.これを参考にして,まずはモデルの完成を目指す.（Marcではm系でモデルを作成していたため,Salomeも同様に作成する.）~
角田さんと秋山さんの論文を読んで,内容を理解することからスタート.(6月初めには結果を出したい.)
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/6honmodel300xy.PNG
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/3D6honmodel300.PNG
・ひとまず,Salome-Mecaでモデルの作成は終えた.メッシュも切り終えたので,Asterstudyの設定方法を考える.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/Salomemodel300.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/Salomemodel300yoko.png
Salomeでの解析方法~
まず,モデルが完成したら,静的解析（ダメそうなら,時刻歴応答解析で）で主桁中央鉛直変位のデータを取ってみてMarcと比較する.(Marcで導出した変位の値は-1.17m)~
材料非線形を考慮しているため,AnalysisでSTAT_NON_LINE（動的ならDYNA_NON_LINEかな。）を用いて解析を行う.~
また、プレストレスは一旦温度応力で試してみようと考えている.~
**解析記録 [#gb613b11]
モデル中央主桁鉛直変位を求める。~
2/26~
ケーブルのみを材料非線形を考慮し、静的解析を回したところ、結果が得られた。（動的解析だと,step数を細かくする必要があるエラーメッセージが出て、解析が回らない。）~
得られた主桁鉛直変位は-25.1263m~
以前、トラス要素を用いた場合の主桁鉛直変位は-1.34mとMarcと近い値（-1.17m）が出ていたのだが...。~
静的解析だと普段使っている温度によるプレストレスを入れられない説もあったため、動的解析（100sで荷重係数1となるように回したが,60sで解析停止）で解析が回る区間の結果を出力したところ、解析途中(60sのデータ)ですら主桁鉛直変位-15.2196mが得られた.
20sごとに約5mごと鉛直変位が下がっていることを踏まえると、静的解析で得られた結果はSalomeの解析では妥当なものであると言える？（温度も入ってそう.）
...ということは、Marcの梁要素とSalomeのBAREE（トラス要素）は似た性質を持っている可能性がある？~
ただ、簡易モデルでも記載したが、ケーブルは曲げ剛性とねじれ剛性を無視するため、それが悪さしている可能性もあり得る.~
→原因は桁とケーブルを接続するバネの剛性が低いことが原因だと判明。主桁中央鉛直変位は-1.50mと結果が得られ,MarcやケーブルBAREE要素と近い値が得られた。~
6/26~
前回,主桁が持ち上がった原因として,ケーブルの半径の与え方を間違えていたことだと判明した。（この橋梁モデルはm単位で作成しているのに、ケーブルの半径の値をm換算するのを忘れていた。）~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/model300-2.png
得られた主桁鉛直変位は-1.34mだった.（Marcで導出した変位の値は-1.17m）
現在、与えている境界条件~
死荷重：桁：122000N/m,主塔：60000N/m （秋山さんの論文と残されたMarcのデータを参照）~
活荷重：主桁中央に集中荷重LP1=1.175×10^6N, 桁全体にLP2=35250N/m (Marcのデータを参照）~
→活荷重の計算方法：道路橋示方書のL荷重(B活荷重)の計算手法に従う.~
L荷重：大型車を連行させた分布荷重にしたLP1と大型車以外の自動車荷重（その他の部分）を分布荷重にさせたLP2の合計.(LP1は集中荷重に直した.)
L荷重の図を見たほうがわかりやすいと思うので,[[この検索ページから:https://www.google.com/search?q=L%E8%8D%B7%E9%87%8D+%E5%9B%B3&client=ubuntu-sn&sca_esv=761560977b3b699f&channel=fs&biw=1920&bih=906&ei=ORo0aMeAA7eUvr0PnsWlsQg&ved=0ahUKEwiHgYSp0MCNAxU3iq8BHZ5iKYYQ4dUDCBA&uact=5&oq=L%E8%8D%B7%E9%87%8D+%E5%9B%B3&gs_lp=Egxnd3Mtd2l6LXNlcnAiC0zojbfph40g5ZuzMgUQABjvBTIIEAAYgAQYogQyBRAAGO8FMgUQABjvBUjcO1DuGViGNXAEeAGQAQCYAYQBoAGIBqoBAzEuNrgBA8gBAPgBAZgCCqACigbCAgoQABiwAxjWBBhHwgIFEAAYgATCAgYQABgFGB7CAgUQIRigAZgDAIgGAZAGCpIHAzQuNqAHlgyyBwMwLja4B-QFwgcFMi02LjTIB0k&sclient=gws-wiz-serp#vhid=anPAULIdubmghM&vssid=_QRo0aNzDJ9a2vr0PwZPnAQ_42]]調べてみてください.~
LP1=(10×5.5×10)+{5.0×(18-5.5)×10}=1175kN~
(10×5.5×10):主載荷荷重×載荷幅×載荷長~
{5.0×(18-5.5)×10}:{従載荷荷重×(舗装部分の長さ-載荷幅)×載荷長}~
従載荷荷重は主載荷荷重の半分.~
LP2=(3.0×5.5)+{1.5×(18-5.5)}=35.25kN/m~
(3.0×5.5):主載荷荷重×載荷幅~
{1.5×(18-5.5)}:従載荷荷重×(舗装部分の長さ-載荷幅)~
CW: 80000N/m (秋山さんの論文を参照)~
ケーブルプレストレス：C1,C26：-150℃ それ以外のケーブル：-100℃ (Marcのデータを参照）~
6/19~
ケーブルをトラス要素に変えたら解析が回ったが、主桁中央鉛直変位が1.00mという解析結果が得られた。~
論文での結果を見てみると、下にたわんでいる結果が得られていたため、荷重が上手く与えられていないのでは？
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/model1.png
・死荷重を入れるのを忘れていた...。Marcで作成していたモデルを見ると,重力を（重力荷重？）モデル全体に与えていたようだが...。（密度×体積×重力加速度で導出ってことだと思う。）~
ただ、桁にかけている死荷重のデータは概要に書いてあるため、早速修正する。
・バネでのやり方だと中々解析が回らないため、やり方を変えてバネを外してケーブルをトラス要素に変えてみた。~
トラス要素：軸力のみが生じ,部材と部材の共有節点をピン接合した要素。~
自分の卒論で取り扱った簡易モデルのケーブルをトラスに変えてみたところ、桁の先端変位（2.50mm）とケーブル張力(16N)がバネ使用時と一致し、有限要素解析でケーブルをトラス要素として扱っている例をいくつか見たので、トラス要素を使うこと自体、信頼性はあると思うが、ケーブルを要素分割した際に解析は回るのかな？（トラス要素の定義に従うなら、要素分割したらケーブルとケーブルをピン接合するということになると思うが...）
6/6 問題になっている可能性が高い箇所を見つけた.Analysisの部分で自由度に関する問題がノード199(C21と塔を繋ぐバネ部分)で発生しているらしい。~
バネ絡みの設定がおかしいと思うが...。~
思いついた解決策として,バネとケーブルの共有節点も離散要素（マトリクスを与える状態、要するにバネ要素と同じ設定）を与えてみるとか？
↓赤丸で囲んだ箇所が該当箇所~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/N199-1.png
↓赤丸の拡大図~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/N199.png
↓メッセージ内容~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/error3.png
6/4 現時点で出ているエラーメッセージ Analysis,DYNA_NON_LINEでエラーが発生中~
(1)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/bridgeerror2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/JPbridgeerror.png
5/25 何とか解析回すステージには立てた.（時刻歴応答は自分,静的解析は山本君にお願いしている。）~
ここからエラーメッセージを順番に潰していく.
5/23 試行錯誤しながらプログラムを組んでいる最中。~
何とか今週までに解析回せるところまで行きたい。
*影響線解析 [#h9310a98]
影響線：移動荷重に対してある点に生じる物理量の大きさを移動荷重の位置関数で表したもの~
影響線の定義から,まずは移動荷重を求め方が分からないと話にならない？
↓橋の上を走る電車を例に移動荷重を算出する方法を見つけた。荷重条件に式を与える方法を使って,動的解析をする必要があるらしい？~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/movingload.png
[[移動荷重について:https://forum.code-aster.org/public/d/13262-how-to-calculate-the-moving-loads-on-a-bridge-in-code-aster]]
[[移動荷重について_その2:https://forum.code-aster.org/public/d/19546-unilateral-constraint-and-moving-force-with-a-given-velocity/8]]
AFFE_CHAR_MECA_F→FORCE_POUTREを使い、DYNA_NON_LINEなどのコマンドを使う動的解析でできそう？~
ちなみに...midasだと[[説明書の10章:https://in.midasit.co.jp/wp-content/uploads/2020/07/constrcution-83-civil-book.pdf]]を見たところ,影響線を導出する手法が記載されていた.
*model300に部分腐食を与えた終局強度解析（後期） [#ce4ff1df]
**静的解析での検証 [#u4e27946]
○線形解析の場合
これまで解析が収束しなかったのは,ケーブル破断後の主桁応力変動が解析上では計算しきれないため起きた可能性がある.~
そこで桁のみを線形状態にした場合で部分腐食を与えた解析を取り組んでみた.
これまでk=3.0になるまで荷重を載荷していたので,k=3.0の荷重を静的荷重(k=1.0)と定め,その荷重に向けて徐々に荷重を増やしていく.
-k=0.9(設計荷重の2.7倍の力が掛かっている.)で解析停止.解析停止までの様子は以下の通り.
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/linearm300.gif
k=0.9(設計荷重の2.7倍の力が掛かっている.)時の橋梁鉛直変位~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_linearzentai.png
-主径間中央_主桁の応力変化をまとめた.~
下縁部~
k=0.54で降伏応力を超えている.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/linearsitaketa.png
中立軸~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/mannakaketa.png
上縁部~
k=0.60で降伏応力を超えている.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/linearketaue.png
こちらでも同様に下縁部→上縁部の順に降伏応力に達している.~
しかし,線形解析でも中立軸は降伏に至ってない.
橋長に沿った主桁中立軸の応力状況~
側径間右側で一番大きな応力が出ている？~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/tyurituketazenbu.png
・線形解析で得られたこと~
①大きな設計荷重を掛けても中立軸（付近）は降伏応力に達しない.(少なくとも,1本のケーブルが腐食しただけでは.)~
これを他のケーブルにも腐食を与えたり,連鎖的にケーブルが切れるといった挙動が起きるとまた結果は変わってくるはず...
②主塔の真下（橋長420m）で大きな圧縮応力が出ている.
これまで主径間中央で最も大きな応力が出ていると思ったが,この考えは違うかもしれない？~
**解析停止時の応答値 [#b9db6bd9]
線形解析での結果を踏まえ,主塔付近の応力を確認するため,解析が停止するまでのモデル全体のアニメーションと主桁断面の応力状況等をまとめた.
C1上端部分腐食時のアニメーション（腐食場所関係なくモデルの挙動はほとんど同じ.）
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/koukiC1uecoro.gif
○橋長方向と鉛直変位の関係~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/k185DEPL.png
○主桁下縁部応力と橋長方向の関係~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/k185_ketasitaSIGM.png
○主桁中立軸応力と橋長方向の関係~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/k185_ketaSIGM.png
-C13部分腐食の結果
○橋長方向と鉛直変位の関係~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/k170_ketaDEPL.png
○主桁下縁部応力と橋長方向の関係~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/k170_ketasitaSIGM.png
○主桁中立軸応力と橋長方向の関係~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/k170_ketaSIGM.png
**腐食状態ケーブル（腐食率50%）をモデルから抜いた場合の解析（連鎖崩壊解析手法を用いた） [#r485660d]
手法は連鎖崩壊解析と同様.これまで破断させる応力ひずみ関係を腐食させたケーブルにのみ与えていたが,他のケーブルにもその関係を与えて解析を回す.
○C1なしモデル
250sでC1ケーブルを破断させ,260sまで解析を回した.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/chainC1hadan.gif
先行研究では,252s以前で塑性ヒンジが見られたため,252sまでの応答値を確認した.
○251.55s時の変形図~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/251_55_DEPL.png
○251.55s時の主桁応力~
最初に,C3,C4間の主桁中立軸で降伏応力に達したことを確認.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/chain_ketaSIGM.png
○C1側側径間ケーブルのひずみ~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/chain_epsi.png
○C13なしモデル
*model300に部分腐食を与えた解析（前期） [#w3f7c4d5]
**Salomeでの解析 [#yb204e33]
終局強度を考慮した解析を行う.(先輩の事例を見ると,静的解析をやっていたのでとりあえずそっちで。)~
終局強度とは？:構造物が破壊（崩壊）する直前の強度。~
終局強度の式は以下の通りである。~
P=Pr+k(D+CW+L)~
ただ,過去の研究ではP=Pr+k(D+L)で終局強度を求めている場合も。~
どっちを使うべきだ？~
Pr:ケーブルプレストレス,D:死荷重(桁),CW:カウンターウェイト,L:活荷重,k:荷重増加係数~
まず,ケーブルプレストレスを載荷し,その後に他の荷重をkを増やすことで増加させていく。~
構造物崩壊時のkを終局荷重係数とする。→この係数を求める！
→これをC1下端,中央,上端に部分腐食を与え,主桁中央変位と腐食部分のC1張力と応力ひずみを見る。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/model300partialcoro1.png
断面形状は以下に示す.~
主桁~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/keta.png~
塔（下部）~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/tou.png~
塔（上部,ケーブルと繋がっているところ）~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/hutozu.png~
***ケーブル腐食なしの終局強度(幾何学的非線形考慮しない) [#h30a0189]
・主桁中央変位と荷重係数の関係.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/tyuouDEPL_nonlinearonly.png
・主桁中央応力と荷重係数の関係.~
K=1.81で下縁部降伏,k=2.60で中立軸付近降伏,k=2.64で解析ストップ~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/nonlinear_stress.png
***ケーブル腐食なしの終局強度(幾何学的非線形考慮) [#gbb3d070]
部分腐食の解析を行う前に,腐食なしで先輩方が解析で求めた終局強度と同じ(近い)値を調べた。~
このとき,終局荷重係数はK=2.63だったが,Salomeだとk=0.30程度で解析エラーが出てしまう。(荷重係数が1になる前に解析が止まってしまうのは流石におかしい。)~
解析が停止する原因を調べたところ、主桁に材料非線形を考慮しない場合（線形弾性）は解析が回ったため、主桁に関する設定が何かしら悪いことをしていると思われる。~
断面のメッシュの切り方に問題がある([[参照ページ:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?Salome-Meca+%E3%83%93%E3%83%BC%E3%83%A0%E8%A6%81%E7%B4%A0#jfa31dbe]])可能性もあるので、メッシュを切り直して再度解析を回す。~
(主桁断面は以下の写真,左側の形状.秋山さんの修論日誌から引用.)~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/akiyama/syutiudanmen2.png~
塔・桁・ケーブルに幾何学的非線形を考慮.
解析方法~
ケーブルプレストレスを荷重係数1まで、他の荷重は荷重係数を1を超えても載荷し続ける。~
〇解析結果
[[比較データ:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?%E8%A7%92%E7%94%B0%E3%81%AE%E4%BF%AE%E8%AB%96%E6%97%A5%E8%AA%8C#f0558fbf]]~
先輩が解析した際の崩壊までの挙動（角田さんの修論日誌より引用）~
k=2.64のとき，崩壊した．
主桁および主塔の変形は，側径間のケーブルC4がk=1.90にて第一降伏点に達するまでは，
線形増加する．k=2.24で主塔位置の主桁下フランジが，k=2.26にてスパン中央の主桁下フランジが降伏する．
その後，スパン中央の主桁断面の塑性化が進む．そして，k=2.61で側径間ケーブルC3，C4が第二降伏点に達し，
k=2.62で，スパン中央の主桁が全断面降伏し塑性ヒンジが形成され，橋全体が崩壊する．
まとめると、k=2.24で主塔位置の主桁下フランジ、k=2.26で橋梁中央部の主桁下フランジが降伏,K=2.62で橋梁中央部の主桁下フランジに塑性ヒンジが形成。~
解析ソフトの差異で、kの値に多少違いは出ると思うが、大きくずれたらSalomeの結果が間違っている。~
7/9~
結果まとめ(比較図を載せたかったが、先輩の論文に載っているものなので割愛.)~
・荷重係数1での橋長鉛直変位(Salome:-1.098m,Marc:-1.129m)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/model300_k1DEPL.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/Marcmodel300DEPL.png
・主桁中央変位と荷重係数の関係.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/model300tyuouDEPL.png
・主桁中央応力と荷重係数の関係.~
K=1.83で下縁部降伏,k=2.51で中立軸降伏,k=2.60で解析ストップ~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/model300syukyoku_k_SIGM.png
Marcに近い挙動の結果は得られた.
7/2~
新しい結果が得られたので、そちらを載せる.~
k=2.61で解析エラー(k=2.61で橋梁モデルが崩壊.)
主塔の死荷重とケーブルプレストレスはk=1.0を維持し,解析.~
それ以外の荷重はkを増やし続ける.~
角田さんの結果と比較して...~
・主桁下縁部降伏応力は,ほぼ一致(k=1.8,私の結果k=1.9)~
ただ,全断面降伏(塑性ヒンジ形成)は私の結果がk=2.16に対し,角田さんの結果はk=2.63.~
ここでずれたのは何故？(塔の死荷重の有無,CWの値を新しいモデル用に合わせているのと,当時のモデルはアンカレッジケーブルがあったため、それらの違い?：当時の側径間ケーブル7本,)~
→3つめに関しては,側径間ケーブルを1本なくした代わりにC1,26の張力を強めにしているので,あまり大きな影響はなさそうな気がする.~
→追記,桁断面の中立軸は原点ではなく,高さ方向から少しずれた位置にあることが判明.(原点中立軸なのは,上下左右対称断面のとき.)~
・C1応力変化はほぼ一致.~
○変形図(scalefactor５倍)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/m300k.gif
○変位図(k=1.0のとき)~
中央鉛直変位は-1.01m,(Marcのmodel300:-1.129m)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/model300k1_heni.png
Marcの場合：-1.129m~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/Marcmodel300DEPL.png
○主桁中央鉛直変位(kの変化)~
崩壊寸前の鉛直変位は8.37m.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/model300k_DEPL.png
○主塔水平変位(kの変化)~
崩壊寸前の水平変位は1.53m.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/model300k_tou_DEPL.png
○モデル中央,主桁応力(k=1.9で下縁部降伏,k=2.16で全断面降伏)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/model300k_stress.png
○モデル中央,主桁下縁ひずみ(右図,降伏ひずみ付近を拡大)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/model300k-EPSIsitaall.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/model300k_epsi_sita.png
○モデル中央,主桁中立軸ひずみ(右図,降伏ひずみ付近を拡大)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/model300k_EPSIceall.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/model300k_epsi_mannaka.png
○C1応力~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/model300k-C1SIGM.png
6/30~
主桁が降伏に至らない問題は、どうやら桁と塔に与えている応力-ひずみの降伏ひずみが違っていたことが原因だと分かった。~（連鎖に関する過去の論文、角田さんの論文から応力ひずみの線形関係の直線部分の傾きから降伏ひずみを計算して判明。）~
...この1,2週間で色々問題点が見つかったな。これで自分の求めたい結果が得られるか？
6/26~
新たな問題に直面.載荷条件の順番を変えることで異なる変位や応力が出てくる。詳細は[[こちら（載荷条件について）:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?Salome-Meca+%E3%83%93%E3%83%BC%E3%83%A0%E8%A6%81%E7%B4%A0#q683524a]]~
これは早急に何とかしないとな...
6/25~
再度、設定を見直したところ、以下2点の変更点を加えた。
・ケーブルと桁を繋ぐバネの剛性が低いため、桁が大きくたわむことが判明した。→変位に対する剛性を高くした。~
・Marcで得られていた張力からプレストレスとして掛けている温度を計算して、得られた温度を各ケーブルに載荷した。（主桁中央鉛直変位の値が-2.5mとMarcで得られていた結果と1mほどずれがあるため）~
これらの変更を加え、ようやくMarcやSalomeのケーブルをBAREEにした場合と近い主桁中央鉛直変位（-1.5m）が得られた。
この状態で終局強度を求めようとすると、k=2.60で収束エラーが出たため、その手前の解析結果を見たところ、主桁中央部がどこも降伏していないことが判明...。（171MPaしか出ていない。）~
また、収束エラーが出たタイミングでモデルが崩壊している可能性が高いが、その様子が見ることができないのも問題。（桁断面のメッシュを細かくしてみる。）~
6/18~
・解析はk=4.17で停止~
グラフは角田さんの論文と同じように作成。
[[角田さん:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?%E8%A7%92%E7%94%B0%E3%81%AE%E4%BF%AE%E8%AB%96%E6%97%A5%E8%AA%8C#f0558fbf]]の結果と異なるものが得られた.(K=2.63)~
主桁に至っては,中立軸付近は降伏していない.~
モデルに関しては,塔の断面の材質（sm490(315MPa)→sm490Y(355MPa)）と断面形状がは違うが、それだけで大きく結果が変わるのだろうか？~
→多少の断面情報で角田さんの結果と大きく変わることはない。~
変形図~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/k_warpbyvector.png
主桁中央鉛直変位と荷重係数~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/k_syukyoku_DEPL.png
スパン中央応力と荷重係数~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/k_syukyoku_SIGM.png
***ケーブル全体腐食の終局強度(幾何学的非線形考慮,要素分割していない.) [#m66d1181]
C1を腐食させる。腐食は50%を想定.このデータを部分腐食と比較.~
・変形図:scalefactor5倍(色付きがC1腐食,白色が腐食なし)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1coro50vsnotcorozentai.png
・主桁中央変位と荷重係数の関係.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coro50DEPL.png
・主桁中央応力と荷重係数の関係.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/SIGMcoro50_sitatyuou.png
上縁と下縁の応力ひずみも見た。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coro50_sitaSIGM.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coro50_ueSIGM.png
***ケーブル全体腐食の終局強度(幾何学的非線形考慮しない) [#h108caf7]
・主桁中央変位と荷重係数の関係.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/allcoro50_DEPL.png
・主桁中央応力と荷重係数の関係.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/allcoro50_SIGM.png
***部分腐食の終局強度 [#pff8fa15]
まず、全ての要素に幾何学的非線形を考慮した場合,簡易モデルで用いた要素分割法(腐食した1要素と腐食していない1要素)で解析が止まってしまった。(腐食ケーブルが降伏に達した時点で解析ストップ.)~
ケーブルだけ幾何学的非線形を考慮する場合(GROT_GDEPとPETITが混在している状態でも),少し荷重が掛かっただけ(k=0.4くらい)で解析が止まってしまう.
試しに,幾何学的非線形を全ての要素から抜いた状態で解析を回したところ,この場合であれば,簡易モデルで用いた要素分割法でも腐食ケーブルが降伏したあとも解析が進んだ.~
幾何学的非線形を考慮するかは一旦置いといて,材料非線形のみを考慮した場合で解析を進める.
モデルの状態~
・モデルの寸法や材料特性はケーブル1要素と同じ.(変更点として,各ケーブルを100要素に分割した)~
・腐食は50%を与えている.~
腐食の与え方は以下の通り.応力ひずみを腐食率ごとに減少させていく方法.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/corrosiongraph.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/corrosion_data.png
ケーブルが腐食により応力が減少する（破断する）という挙動は,静的解析では再現できないみたい.(ケーブルが切れる再現だから,力の釣り合いから外れるから？)~
→静的解析ではk=1.74で解析が停止.恐らく,ここでケーブルが切れていると思われる.
そのため、動的解析を用いて終局強度を求めてみた.~(100sで荷重係数1になるように設定.つまり,1s:k=0.01を想定.)~
○C1下端部腐食部~
k=1.74でC1が破断.(静的解析で解析停止したタイミングでケーブルが切れたことを確認.)~
橋梁中央の変位(270m)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coro50model300DEPL.png
主桁中央の応力と荷重係数~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/keta_coro50_SE.png
C1の健全部の応力と荷重係数~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/sita50coro_C1kenzenSIGM.png
C1の腐食部の応力と荷重係数~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/sita50coro_C1coroSIGM.png
C13張力~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/50sitacoro_C13T2.png
○C1真ん中腐食部~
k=1.74で破断~
橋梁中央の変位(270m)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/mannaka50coro_DEPL.png
主桁中央の応力と荷重係数~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coro50_mannaka_ketaSIGM.png
C1の健全部の応力と荷重係数~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coro50_mananakahadan_kenzenC1SIGM.png
C1の腐食部の応力と荷重係数~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coro50_C1mannaka_coroSIGM.png
C13張力~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coro50mannaka_C13T.png
○C1上端部腐食部~
橋梁中央の変位(270m)~
k=1.74で破断~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coro50_ue_DEPL.png
主桁中央の応力と荷重係数~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coro50ue_ketaSIGM.png
C1の健全部の応力と荷重係数~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/uecoro_C1kenzenSEGM.png
C1の腐食部の応力と荷重係数~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1uecoro50_coroSIGM.png
C13張力~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coro50ue_C13T.png
・分かったこと~
①C1ケーブルが腐食が達した時点でC1健全ケーブルの応力がほとんど変化しないことが確認できた.(何故、こうなるかは考え中.)~
ちなみに,他のケーブルが健全状態の応力ひずみ関係に沿った応力ひずみ変化になる.
②1本切れただけでは塑性ヒンジに達しない.
③ケーブル破断時に他のケーブルに張力変化が得られたことも確認.
④腐食位置毎に破断時の荷重係数に変化はなかった.
結果~
7/21更新~
k=2.05で解析が止まった.しかし,桁は中立軸は降伏していないため,構造物は崩壊に至っていないと言える.
C1腐食部ケーブル応力~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coro50C1SIGM.png
C1健全ケーブル応力~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/kenzenC1SIGM.png
桁の応力~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coro50_ketaSIGM.png
***動的解析における終局強度を求めるC1全体腐食の破断時の応答値 [#i949e212]
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1zentaicoro.png
断面下縁部降伏：k=1.79~
断面上縁部降伏：k=1.91~
ケーブル破断：k=2.3~
k=2.5で解析停止.
・中立軸が降伏する前に解析が停止したが,そのときの降伏応力が231MPaの結果が得られた.~
恐らく,k=2.5付近で全断面降伏すると思われる.
***動的解析における終局強度を求めるC1部分腐食の破断時の応答値 [#f4d17100]
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1bubuncoro1.png
ケーブル上端・中央部に腐食を与えた場合：k=1.74で破断.~
ケーブル下端に腐食を与えた場合：k=1.8で破断.
解析自体はk=1.85で停止.
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1bubuncoro2.png
腐食箇所関係なく断面下縁部・上縁部は降伏応力235MPaに達した後,応答値が変動.~
中立軸は解析停止k=1.85でも降伏応力に達しないため,塑性ヒンジは見られなかった.
***動的解析におけるC13に全体腐食を与えた場合の応答値 [#f4d17100]
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C13zentaicoro.png
断面下縁部降伏：k=1.72~
断面上縁部降伏：k=1.84~
ケーブル破断：k=2.11~
k=2.34で解析停止.
中立軸が降伏応力に至る前に解析停止.しかし,応力変動がC1全体腐食モデルと似ており,C1全体腐食より小さな荷重係数で断面降伏が始まっていることから,k=2.5より小さな近い値で終局強度に至るのでは？
***動的解析におけるC13に部分腐食を与えた場合の応答値 [#rad58a9d]
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C13bubuncoro1.png
k=1.52でケーブル破断.~
k=1.70で解析停止.~
ケーブルが切れるタイミングはほぼ同じだが,変位の応答変化は上端・中央,下端で異なる結果が得られた.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C13bubuncoro2.png
腐食箇所関係なく断面下縁部・上縁部はケーブル破断後に降伏応力235MPaに達する.その後,応答値が変動.~
中立軸は解析停止k=1.70でも降伏応力に達しないため,塑性ヒンジは見られなかった.
***まとめ [#b907df73]
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/kmatome.png
・部分腐食モデルでは,終局強度を求めることは現状できない.(求め方に工夫が必要.)~
・主径間中央に近づくほど,小さな荷重係数でケーブルが切れる.~
・部分腐食モデルのほうが小さな荷重係数でケーブルが切れる.~
**Marcでの解析 [#bc35888e]
3/3~
model300を用いて解析を回してみた.
まず,腐食率50%の破断時の応力ひずみ関係をC1端部にのみ与えた結果を載せる.~
https://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/model300c1hadan.gif
腐食したケーブルのみ動く現象が起きている.
次に,橋梁全体を腐食率25%,C1端部に腐食率50%の破断時の応力ひずみ関係を与えた結果を載せる.(解析は84s程度で止まってしまった.)~
これは連鎖崩壊は全てのケーブルがある程度腐食していないと発生しないため,それに近い状態だとどうなるか見るため行った.
https://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/model300c1hadanzentaicoro25.gif
全体的に下にたわみすぎだと思う.[[この問題:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?%E7%A7%8B%E5%B1%B1%E3%81%AE%E4%BF%AE%E8%AB%96%E6%97%A5%E8%AA%8C#cc7ba732]]と似ている現象が起きてそう？~
断面2次モーメントの値を調整してみて再度解析を回してみる.~
*簡易モデルを用いた検証 [#lce2339d]
**卒論で扱った簡易モデル [#n814c429]
Q:そもそも,簡易モデルって何？
A:簡易モデルとは桁をケーブルで1本（または2本）で吊ったモデルである.~
学部生のときにこのモデルを使って,DAFを算出し,衝撃力が出ているかの確認を行なった.
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/chiyooka/Salomekani.png
桁は高さ10mm,幅50mm,厚さ2mmの薄肉矩形断面,ケーブルは半径4mmの円形断面である.
Salomeの設定とほとんど近い状態でMarcでも解析してみた.（プレストレスの与え方だけ違う.Marc→温度応力,Salome→ひずみ）~
また,Marcだと149.9sで解析が終わってしまうため,0.01s分データがない結果である.
ケーブル1本のみの簡易モデル（集中荷重なし）の結果
Marcの解析結果~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/CnasimarcmodelA.png
Salomeの解析結果~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/modelADEPL.png
DAF~
,時間,DAF(Marc),DAF(Salome)
,0.005s,1.97,1.95
,0.01s,1.95,1.93
,1s,1.02,1.02
,5s,1.00,1.01
1ケースしか試していないが,どちらの結果も見る限り,大きな差は出ていないため,解析ソフトによる差異はほとんどないと言えるだろう.~
この結果から,ピン接合とプレストレスをSalome用に設定をアレンジしたものもそれなりに信頼性は持てると思う.~
***DAFとtimestepの関係 [#w4d965e5]
2本モデルのとき、モデルに冗長性が出たためDAFの値が2に近い値が出なかったと考察したが...~
timestepを更に細かくしたら結果が変わるのでは？と考え、timestepを0.001に変えて再度解析を回してみた。（卒論ではパソコンの性能が足りなくて0.005までしか結果を出せなかった。）
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/kaniver2.png
解析時間を115sに変更,100sまでゆっくりと載荷し,110sで一番外側のケーブルを破断させ,115sまで解析を行った。
結果はケーブル破断時間を0.001s,0.005sのもの2種類用意した。
timestep0.001の結果~
109s〜115sの様子
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/2honmodelhadan0001.png
110s〜110.5sの拡大グラフ
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/2honkanimodelhadan0001kakudai.png
timestep0.005の結果~
109s〜115sの様子
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/2honkanimodelhadan0005.png
110s〜110.5sの拡大グラフ
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/2honkanimodelhadan0005kakudai.png
DAFはtimestep0.001sでは1.87~
timestep0.005sでは1.47となった。
以上から、2本モデルの簡易モデルでDAFの値が思ったより出なかったのは,timestepが原因だと考えられる。~
1本モデルでtimestep0.001に変えて解析回してないからはっきりとは言えないけど...timestepに大きな差（0.005と0.001）がないのに、ケーブル本数増えたらDAFの数値が変わるなら橋梁モデルももっとtimestepを細かくしたらDAFも変わってくるのかな？
**ケーブル2本の簡易モデルに部分腐食をさせてみる [#a1bd3253]
***Marc [#q11b79aa]
3/10~
大ひずみ,同条件で解析を回したが,29.03sでエラーで止まってしまった.(エラーメッセージはstep数をもっと細かくとのことだが,自動かつ細かいところで0.0001sで切ってくれているが...)
最終挙動（29.03s）と桁先端変位と腐食部C1張力は以下の通り.~
最終挙動（腐食したケーブルが伸びた後,桁が揺れる挙動が得られた.）~
https://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/I0115000n2honmodel_SIMOMIEHE.gif
C1腐食部張力,腐食していない部分も同じ値に...(Tの最大値:6172.5N)~
https://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/I01kanimodel15000N101sT_SIMOMIEHE.PNG
先端変位(-1.368mmたわんでいる)~
https://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/I01kanimodel15000N101sDEPL_SIMOMIEHE.PNG
最終挙動付近を拡大すると,桁が揺れていることもわかる.~
https://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/I01kanimodel15000N101sDEPL_SIMOMIEHE_kakudai.PNG
3/9~
微小ひずみ,集中荷重15000Nを掛けた際の解析結果だが，最終挙動（101s）と桁先端変位と腐食部C1張力は以下の通り.~
最終挙動（時間を掛けて,下にたわんだだけ.複雑な動きはしていない）~
https://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/I01kanimodel15000N101s.PNG
C1腐食部張力,腐食していない部分も同じ値に...(Tの最大値:21083.3N)~
https://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/I01kanimodel15000N101sT.PNG
先端変位(-4.48mmたわんでいる)~
https://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/I01kanimodel15000N101sDEPL.PNG
断面2次モーメントを調整するやり方（ケーブルに適用）を参考にして解析を回したが，ケーブルの降伏応力（ミーゼス応力）と塑性ひずみが0という結果に.~
設定箇所をよく見ると，特性：入力値の場合,材料挙動が線形弾性のみ（応力とひずみが常に比例の関係：SalomeでいうELAS）しか対応していないみたい.(この設定では材料非線形を考慮した挙動が見れないってこと？)~
→3/10追記：桁の降伏応力と塑性ひずみの数値が見れたため，下記写真のやり方ではこれらの数値は見れないと思われる.
https://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/tyujitudanmen.PNG
3/4 [[この問題:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?%E7%A7%8B%E5%B1%B1%E3%81%AE%E4%BF%AE%E8%AB%96%E6%97%A5%E8%AA%8C#cc7ba732]]と似ている現象が起きてそうなので、参考にしたところ、断面2次モーメントの値を調整したら改善したみたい.~
どうやらケーブルのたわみは曲げ剛性を無視し,伸び剛性と張力の影響でたわむとのこと.このモデルのケーブルの断面2次モーメントを計算したら201.061mm^4であり,無視できないほど大きいため,断面2次モーメントをIxx=Iyy=0.1にした状態で桁先端に集中荷重100Nを掛けて再度解析を回してみた.
https://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/I01kanimodel.PNG
写真は解析終了時（101sの様子,荷重が掛け終わるのは100s）の様子.解析は最後まで回り,下にたわむことはなくなった.~
この状態でSalomeに与えていた集中荷重15000Nでどういった挙動を見るか再度確かめる.
2/27~
MarcでSalomeと同じ条件で解析を回したところ,以下の結果が得られた.（途中で解析が止まってしまうが，腐食時の応力が0になっただろうと思われる挙動が得られたため，結果を載せた．）
アニメーション(主桁先端に15000N)~
https://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/15000N2honmodel2.gif
Marcだと,100N掛けても切れるような挙動が得られる.~
https://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/100n2honmodel.gif
Salomeと挙動が異なる点として，要素の違いがあると思われる．現状,Salomeでは材料非線形を考慮した梁要素の解析はPOU_D_EMやPOU_D_TGMといったマルチファイバー要素でしかできず,
しかも,マルチファイバー要素自体,圧縮が出てしまうことから,橋梁のケーブルをマルチファイバー要素で再現できない可能性がある.~
Marcでは，先輩方が使っていた橋梁モデルのケーブルに要素タイプ25を用いていたため,それを使い解析を行った.以前,静的解析での橋梁モデルをSalome:ケーブルをBAREE,Marc:ケーブルを要素タイプ25,で[[主桁中央鉛直変位を比較した:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?%E5%8D%83%E4%BB%A3%E5%B2%A1%E3%81%AE%E4%BF%AE%E8%AB%96%E6%97%A5%E8%AA%8C#gb613b11]]ところ,近い値が出ていたため,SalomeのBARREとMarcの要素タイプ25は近い性質を持っているのだと思われる.~
では,MarcでSalomeのマルチファイバー要素に近い要素タイプは何か？以前,[[簡易モデルの桁先端変位とケーブル応力を比較した際:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?%E5%8D%83%E4%BB%A3%E5%B2%A1%E3%81%AE%E5%8D%92%E8%AB%96%E6%97%A5%E8%AA%8C#k829fe9a]]，SalomeとMarcの結果がほぼ一致し，そのときの要素タイプが98だったため,SalomeのマルチファイバーとMarcの要素タイプ98が近い性質を持っているのだと思われる.~
https://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/yousotype.PNG
閉断面：構造部材の断面が完全に閉じている断面（パイプとか）~
開断面：構造部材の断面に隙間がある断面（I形鋼とか）
***Salome-Meca [#u3751526]
現状は、ケーブルに健全時の構成則と腐食時の構成則を同時に与えることはできている。ここから、腐食時のケーブルだけを切る（応力0に落とす）には,要素分割数と伸びの長さが関わる可能性あり。~
材料非線形はケーブルのみにしている。（桁に入れても,同じ材料に異なる構成則を同時に入れることは可）~
2/12~
かつて10分割した要素長と同等の長さを腐食部分、その他の長さを健全部とした。（2要素にした。）~
結果は以下の通り.
グラフはどれも応力が0に落ちた瞬間の箇所をとっている.~
時間と桁先端変位~
https://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/212_2honmodelDEPL.png
時間と腐食部張力~
https://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/212C1T2honmodel.png
応力ひずみ~
応力392MPa,ひずみ0.004で落ちるように設定.
応力がおちるグラフ~
https://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/212C1SIGM_1.png
応力が0に落ちた瞬間の応力ひずみ~
https://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/212C1SIGM_2.png
・問題点
腐食による破断後の張力を見ると、0に近い値が出ているものもあるが,時間が経つほどどんどん大きくなっている...。（最も0に近い値は0.14475N）~
応力については、完全に0にはできない状況（落ちた瞬間は1MPa程度出ている。）~
→より安定したモデル（橋梁モデル）で解析するとまた変わってくる？~
→ケーブルをcable要素に変えてみる？（けど、応力-ひずみ考慮できるの？）~
12/26~
腐食部分の要素長さ100√2mm,50√2mmの場合、解析は回ることが判明した。ただ、結果を見ようとするとフリーズする...。~
ちなみに、要素長さ100√2mmじゃ要素分割10と同じ長さであるため、解析が回る回らないの違いとして、要素長と要素分割数が関わってくる可能性が高い。~
-解析が回る（回りやすい）パターンは以下のとおり~
(1)材料非線形を考慮している要素長は長いほど安定しやすい~
(2)分割数は少ないほうが安定しやすい
12/18~
ゼミで頂いた意見まとめ~
-真ん中のみの腐食はできる？→健全-腐食-健全~
-何故,要素長によって解析できないのか？→解明する。~
次のゼミまでに、要素長と構成則の関係？を少しでも明らかにしないとな...。
12/12~
要素分割2の場合は？~
12/10での解析と同条件で解析を行った.(応力を落ちることによる破断再現に着目するため.)
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coro75youso2.gif
変位と時間のグラフ~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/youso2DEPL.png
応力と時間のグラフ~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/youso2MPa.png
この後,要素分割3,4も回したところ解析自体は回った.ただ,現状の要素分割10は解析が回らないし,部分腐食となるとできるだけ要素を細かくした1要素を腐食させることになると思われるため、順番に要素分割して確認する作業はここまででいいかと。~
要素分割2ができるなら,腐食させたい要素長を短く,健全部を長くするというモデル再現はどうだろうか？(例えると,1000mmの部材であれば,10mm腐食部,残り990mmを健全部にする)~
実際に試したところ、収束条件のエラーが発生した。~
...ということは,要素長と構成則の相性がある？（そんなことありえるのか？）
↓よく出る収束条件のエラーはこちら~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coroerror3.png
12/10~
Salome-Mecaで分割数1の場合でも解析が回らないのか調査してみた。（よくよく思えば、これを先にやるべきだった...。）~
分割数1であれば、Salome-Mecaでも解析が回ることが分かった。（精度面を一旦置いとけば、ひとまずMarcと同じことはできるでいいのかな。）~
桁の先端に集中荷重Fy=-15000N掛けた.
要素分割していないため、全体腐食を表現。(破断時間:97.109s)
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coro75youso1.gif
結果
・桁先端の変位と時間グラフ~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/youso1coro75DEPL.png
・ケーブルの応力と時間グラフ~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/yousoSIGMcoro75.png
気になる部分はいくつかあるが、応力を落とすという再現に着目すればひとまず最低限のものはできているのではないだろうか。~
11/26~
部分腐食再現で個人的にいくつか思いついたことをまとめる~
（応力を0(ほぼ0)に落として,伸びないようにするために、応力-ひずみ曲線と別の切り口がないか模索。）~
①ある時間でヤング率を下げて,ヤング率低下による腐食再現.~
→ELAS_FO(材料特性をfunctionで与えるやり方)でできそうかと思ったが、横軸時間のfunctionに対応していないみたい...。~
②応力-ひずみ関係で応力を落とす行為に問題があるのであれば、応力-時間,ひずみ-時間グラフでも与えてみる？~
→問題はそれらをどこに入れるか.TRACTIONにはfunctionを1つしか入れることができない。~
③DEFI_MATERIAUのFATIGUEを使う？~
→使い方をいまいち理解していないが,ひとまずmultifiberで使えることは確認.~
④今行っている2段階解析で、1段階目を破断寸前で止め、2段階目で破断させるやり方で進める。~
11/20~
気になるコマンドを見つけたため、ここに共有する。~
[[U4.43.01.pdf:https://code-aster.org/doc/v14/fr/man_u/u4/u4.43.01.pdf]]
[[R7.04.01.pdf:https://code-aster.org/doc/v15/en/man_r/r7/r7.04.01.pdf]]
DEFI_MATERIAU→FATIGUEというコマンドで、詳細まで調べてないが、例を出してる人を見ると疲労解析に使っているみたい。~
試しにFATIGUEに適当な設定を入れて,解析を回してみたところ、Analysisまで進んだため、multifiberに対応していないってことはなさそう？~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/fatigue1.png
11/4~
色々なパターン（2本モデル以外だと,1本モデルと塔付きの1本モデル）を試したが、どうも上手く行かない...。~
→1本モデル自体はケーブル一様に応力が出ないため,NG.塔付きだと,収束しない...。~
エラーが出る直前で解析を止めて結果を見ると,構成則通り動いているのは確認できたことから、応力を0近くまで落とす場所でエラーが出る可能性大と予想。~
試しに,応力が一定（伸び続けている箇所）になった後,応力を少し落とすような構成則に変えてみて解析を回すと、エラーが発生したため、応力を落とす挙動で問題があると考えられる.
エラーの内容~
どちらも収束絡みのエラー.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coroerror1.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coroerror2.png
とりあえず思いつく解決策として、応力ひずみ減少ではなく断面積減少で攻めてみるとか考えてみた（過去の先輩方が使っていた手法.）けど、結局応力を0に落とすことができないため、あまり関係ないか？~
10/29~
ケーブルをCABLE要素に変える案が出たため、試してみたが、TRACTIONを用いて応力ひずみを与える、Von Misesの降伏条件（VMIS_ISOT_TRAC）に対応していないことが分かった。（CABLEに対応していないというエラーが出た。）~
10/16~
死荷重のみだとケーブルの応力が塑性域までに達しないため,桁の先端に集中荷重(15000N)を掛け,応力ひずみ曲線を出して与えたケーブル構成則に従って解析できているか確認した。~
以下は,C1の腐食していない部分（青色線）と腐食部分[75%]（緑色線）の応力ひずみ曲線である。~
なお、①（腐食位置）,②の要素からデータを取得した。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/corrosionpoint.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/corrosiongraph.png
・POU_D_TGMの場合
①腐食箇所~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/75corrosionSIGMEPSI.png
②健全部分~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/75notcorrosion_SIGMEPSI.png
・POU_D_EMの場合
①腐食箇所~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/EM_75corrosion_SIGMEPSI.png
②健全部分~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/EM_75notcorrosion_SIGMEPSI.png
グラフ見たところ,健全部分と腐食部分で応力ひずみ曲線が異なる上に腐食率75%のグラフを見るかぎり与えたケーブル構成則通りに解析をしてくれてそう?~
また、張力も見たところ,健全部は19813N,腐食部は19808Nと腐食している箇所のほうが張力の値が多少ではあるが小さくなったことも確認できた。~
また,POU_D_EMに変えても解析は回ったため,プレストレスは温度に変えるほうがよさそう.(昨年あった張力が落ちてしまう問題が気になるが...→EMで起きていた問題なので,TGMの場合だとどうなるか見てみる.)
10/10~
温度に変えると,PRE_EPSIのときみたいにあまりにもおかしい結果が得られなかった.~
とりあえず,腐食率0%と75%のC1,2張力と先端変位の結果を載せる.（結果見ることが優先のため,timestepは0.1にしてある.）~
・C1張力（腐食させるところを抽出）~
腐食0%：6.71237N~
腐食75%：6.71237N
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/corrosion_T_C1.png
・C2張力~
腐食0%：12.0153N~
腐食75%：12.0153N
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/corrosion_C2_T.png
何と、腐食率関係なく張力の値が同じになった。...そんなことあるか？？~
自分のイメージだと,腐食した部分は部材の強さとして弱くなるわけだから,変位,張力ともに小さくなると思うけど...~
こうなってしまった原因だが,TRACTIONで定義した応力ひずみ関係が正しく反映されていない,もしくは反映されているが他の要因で打ち消されてしまっている可能性があると思う.(PRE_EPSIを使うと,そこで与えたひずみしか考慮しない結果になってしまうみたいな.)~
後者の場合,ケーブルに与えている外力は,温度で与えているプレストレスのため,そのあたりの設定が怪しい...?~
10/16追記~
ケーブルの応力を見ると,10MPaすら超えてなかった.原点近くの応力ひずみ曲線の傾きは腐食率関係なくほぼ同じ値であるから,結果に差異が出なかったってことかな？
10/9~
解析結果を確認したところ...あまりよろしくなさそうな結果が得られた.~
腐食0%,25%の結果を載せる.
腐食0%の変形図(25%もほぼ同じ変形図だった.)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/corrosion_kekka1.png
C1張力と時間の結果~
腐食0%の結果~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/corrosiongraph_1.png
40s付近拡大~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/corrosiongraph_1kakudai1.png
60s付近拡大~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/corrosiongraph1_kakudai2.png
腐食25%の結果~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/corrosiongraph_25.png
40s付近拡大~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/corrosiongraph_25kakudai1.png
60s付近拡大~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/corrosiongraph_25kakudai2.png
こうなってしまった原因だが,変形図を見てみるとケーブルが折れ曲がっているが,折れ曲がっている要素はPRE_EPSIを与えている部分である。~
現状,PRE_EPSIを与えた要素はPRE_EPSIのみを考慮してしまうため、それが悪さしているのかも...?~
簡易モデルだと,以前張力が一瞬下がる問題があるが,とりあえず温度をプレストレスに変えてもう一回解析を回してみる.~
10/5~
POU_D_TGMに変えたら解析は回ったのでひとまず結果確認.~
POU_D_TGM:断面の反り(曲げのこと？)のモデリングを考慮したTimoshenko梁の要素の定式化に基づく(Astercode解説書より).Multifiberなので,材料非線形も考慮可.~
10/4~
応力-ひずみ曲線の応力が一定になる箇所でエラーが起きている可能性があるため，少しだけ傾きをつけた．[（1570,ひずみ）→(1571,ひずみ)と変化できるように]~
[[弾塑性解析2021:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?Salome-Meca%E6%BC%94%E7%BF%92_%E5%BC%BE%E5%A1%91%E6%80%A7%E8%A7%A3%E6%9E%90%282021%29]]を参考にした．~
こうすると,上記のエラーは出なくなったが,今度はtimestep絡みのエラー（収束エラー）が出てしまった．~
ここをどうやって解決するか模索中...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/kaniver2.png
要素分割している外側のケーブル（C1）に部分腐食を与えて解析を回す。（内側のケーブル,C2の張力と先端変位を求める。）~
腐食の与え方はケーブルに与える応力ひずみを減少させていく方法。~
計算方法は健全状態（腐食0%）のケーブルの応力-ひずみ曲線の降伏応力,ひずみに腐食していない分を掛けて求める。~
{例:腐食25%の場合:1570×(1-0.25)=1177MPa}
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/corrosiongraph.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/corrosion_data.png
**解析結果 [#h686f34a]
Salomeで解析を進める方向に決まったので、その結果をまとめる。~
初めに解析モデルは以下の写真のものを用いた。また、ケーブルに与える応力-ひずみ関係は10/4に載せた上の写真（すぐ上の写真）の健全状態（腐食率0%）,腐食率75%を使う.~
しかし、腐食率75%に関しては第二降伏点に達した後の伸びが0.01に達する前に解析エラーが発生してしまうため、伸びを短くしてひずみが0.004に達したら破断するように調整した.（下の右側の写真,伸びは0.01以降も続いている.）~
また、ケーブルには幾何学的非線形も考慮している.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/model.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/usecorosion.png
桁断面とケーブル断面
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/ketadanmen.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/cabledanmen.png
腐食はC1ケーブル全てを腐食させた全体腐食,そして、部分腐食はC1ケーブル上部,真ん中,下部にそれぞれ与え,部分腐食①のとき,部分腐食②のとき,部分腐食③のときのそれぞれ3パターンでの合計4パターンで動的解析を用い,応答値（桁先端変位,腐食部C1張力,C2張力,C1腐食部応力ひずみ）を求めた.（部分腐食の長さは10mm。部分腐食①の例だが,何も載荷していない場合の部分腐食の長さはピンクの点の間.）~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coro10mm.png
***集中荷重の場合(0sから集中荷重掛けるパターン) [#n0e49bf6]
載荷条件：桁に死荷重；-0.017248N/mm ケーブルプレストレス;C1に-100℃,C2に-90℃ 桁の先端に集中荷重;15000N（腐食ケーブルに破断するための応力ひずみに達するようにするために載荷.）~
解析時間：110sまで解析を回し,100s掛けて載荷し,その後10s安定させる.（抽出したデータは安定状態から1s経った,101sまで）~
まず,ケーブルが全て健全状態の結果を以下にまとめた.~
挙動(スケールファクタ5倍 色がついている粒？みたいなのが、元々桁のあった位置.)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/Cnomi_kenzen.png
・101sの解析結果~
先端変位：-3.44037mm ~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/kenzen_DEPL.png
C1張力：20361.1N~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/kenzen_C1T.png
C2張力：1048.53N~
ただ、所々で張力が一定に変化しないところがある.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/kenzen_C2T.png
C1応力ひずみ：ひずみ0.00140143,応力410.622MPa（腐食率0%のため,ケーブルが破断するほどの応力ひずみは出ていない.そもそも、第1降伏点にすら達していない.）~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/kenzen_SIGMEPSI.png
次に,腐食の解析結果をまとめる~
全体腐食(97.6sで破断,解析時間101s),部分腐食①(95.40sで破断,解析終了時間95.4076s),部分腐食②(95.41sで破断,解析終了時間95.416s),部分腐食③(95.50sで破断,解析終了時間95.5038s)の解析結果.
・解析終了時の挙動（スケールファクタは全て1倍.部分腐食モデルにおいて,C1ケーブルのピンクの点と点の間が腐食部分を表している。色がついている粒？みたいなのが、元々桁のあった位置.）~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/kanimodel_corozentai.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/kanimodel_coro10mmue.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coro10mm_mannaka.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coro10mm75sita.png
・先端変位（右側の写真は,破断した瞬間の様子を拡大したもの.）
全体腐食（解析終了時の変位：-141.511mm）~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_corozentai.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_corozentaikakudai1.png
部分腐食①(解析終了時の変位：-188.472mm）
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_uecoro.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_uecorokakudai1.png
部分腐食②(解析終了時の変位：-149.295mm）
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_mannakacoro.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_mannakacorokakudai1.png
部分腐食③(解析終了時の変位：-158.937mm）
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_sitacoro.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_sitacorokakudai.png
破断する寸前の桁先端変位
全体腐食：-6.09mm
部分腐食①,②,③：-3.30mm
・C1張力（腐食部）（右側の写真は,破断した瞬間の様子を拡大したもの.）~
張力はどの解析結果でも破断後は100N程度まで落ちる.
全体腐食（最大張力：19441.5N）~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_zentaicoro.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_zentaicorokakudai.png
部分腐食①（最大張力：19437.5N）
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_uecoro.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_uecorokakudai.png
部分腐食②（最大張力：19441.9N）
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_mannakacoro.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_mannakacorokakudai.png
部分腐食③（最大張力：19440.3N）
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_sitacoro.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_sitacorokakudai.png
・C2張力（右側の写真は,破断した瞬間の様子を拡大したもの.）
全体腐食（最大張力：56682.7N）~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C2T_zentaicoro.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C2T_zentaicorokakudai.png
部分腐食①（最大張力：45423.1N）
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C2T_uecoro.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C2T_uecorokakudai.png
部分腐食②（最大張力：49039.2N）
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C2T_mannakacoro.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C2T_mannakacorokakudai.png
部分腐食③（最大張力：60822.2N）
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C2T_sitacoro.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C2T_sitacorokakudai.png
・C1腐食部応力ひずみ（左写真は応力が落ちる箇所までの応力ひずみ曲線,右写真は応力が落ちてから解析が終わる箇所を拡大したもの）~
全体腐食~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_zentaicoro.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_zentaicorokakudai.png
部分腐食①
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_uecoro.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_uecorokakudai.png
応力変化が激しいところを拡大。
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_uecorokakudai2.png
部分腐食②
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_mannakacoro.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_mannakacorokakudai.png
応力変化が激しいところを拡大。
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_mannakacorokakudai2.png
部分腐食③
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_sitacoro.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_sitacorokakudai.png
全体腐食以外のモデルでは安定状態に入る前に95s前後で解析エラーが出てしまう。恐らく、解析では計算できないほどモデルが変形してしまっているからだと思われる.~
そのため、腐食ケーブルが伸びた瞬間(応力がほぼ0になる瞬間)で解析を止め、結果を得た.~
***集中荷重の場合(安定状態を終えてから集中荷重を掛けるパターン) [#r79d5e43]
連鎖崩壊は供用した橋梁が劣化して落橋のため、安定状態にも達していないのにケーブルが破断するのは求めたいこととずれるのではと思い、集中荷重を掛けるタイミングを変えてみた~
死荷重とプレストレスを100s掛け、10s安定させた後に集中荷重15000Nを20s掛けて載荷する方法で再度解析を回した。(安定状態に達したモデルにケーブルが切れるほどの荷重が掛かった際の状況を見る.)~
データの抽出は安定状態を終える110sから行い,全体腐食(129.56sで破断;解析終了時間131s),部分腐食①(129.008sで破断;解析終了時間129.013s),部分腐食②(sで破断;解析終了時間129.017s),部分腐食③(sで破断:解析終了時間119.089s)の解析結果.
・解析終了時の挙動（スケールファクタは全て1倍.部分腐食モデルにおいて,C1ケーブルのピンクの点と点の間が腐食部分を表している。色がついている粒？みたいなのが、元々桁のあった位置.）~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/kanimodel_corozentai2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/kanimodel_coro10mmue2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coro10mm_mannaka2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/coro10mm75sita2.png
・先端変位（右側の写真は,破断した瞬間の様子を拡大したもの.）
全体腐食（解析終了時の変位：-141.793mm）~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_zentaicoro2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_zentaicorokakudai2.png
部分腐食①(解析終了時の変位：-193.838mm）
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_uecoro2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_uecorokakudai2.png
部分腐食②(解析終了時の変位：-149.295mm）
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_mannakacoro2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_mannakacorokakudai2.png
部分腐食③(解析終了時の変位：-158.937mm）
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_sitacoro2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_sitacorokakudai2.png
・C1張力（腐食部）（右側の写真は,破断した瞬間の様子を拡大したもの.）
全体腐食（最大張力：19437.6N）~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_zentaicoro2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_zentaicorokakudai2.png
部分腐食①（最大張力：19344.9N）
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_uecoro2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_uecorokakudai2.png
部分腐食②（最大張力：19439.2N）
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_mannakacoro2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_mannakacorokakudai2.png
部分腐食③（最大張力：19438.7N）
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_sitacoro2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_sitacorokakudai2.png
・C2張力（右側の写真は,破断した瞬間の様子を拡大したもの.）
全体腐食（最大張力：26846.5N）~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C2T_zentaicoro2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C2T_zentaicorokakudai2.png
部分腐食①（最大張力：42880N）
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C2T_uecoro2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C2T_uecorokakudai2.png
部分腐食②（最大張力：41691.4N）
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C2T_mannakacoro2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C2T_mannakacorokakudai2.png
部分腐食③（最大張力：60829.3N）
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C2T_sitacoro2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C2T_sitacorokakudai2.png
・C1腐食部応力ひずみ（右側の写真は,破断した瞬間の様子を拡大したもの.）
全体腐食~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_zentaicoro2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_zentaicorokakudai2.png
部分腐食①
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_uecoro2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_uecorokakudai3.png
部分腐食②
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_mannakacoro2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_mannakacorokakudai3.png
部分腐食③
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_sitacoro2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_sitacorokakudai2.png
***衝撃荷重の場合 [#jfe97a58]
載荷条件：桁に死荷重；-0.017248N/mm ケーブルプレストレス;C1に-100℃,C2に-90℃ 110sでC2ケーブルを破断させる.~
解析時間：131sまで解析を回し,100s掛けて載荷し,その後10s安定させたあとにC2ケーブルを破断させ(110sで破断),そのときに生じる衝撃力による応答値を求める.~
まず、衝撃荷重のみの結果を載せる.全体腐食と部分腐食③（ケーブル下側部分の腐食）の挙動とC1腐食部の応力ひずみを見た.
・挙動
（全体腐食,スケールファクタ10倍 ケーブルプレストレスにより桁が持ち上がった動き.）~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/Cnasikanimodel_impactzentaicoro.png
（部分腐食③,スケールファクタ100倍 C2の衝撃力が伝わったのか桁の真ん中が下がっている.腐食部分と健全部の応力ひずみ関係が異なるため,このような挙動になったと思われる.）~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/Cnasikanimodel_impactsitacoro.png
・C1腐食部応力ひずみ
全体腐食の場合~
プレストレスによってケーブルが桁を引っ張っているため、ケーブルが縮むような（圧縮ひずみ？）挙動が得られるのかと思われる.~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_impactzentaicoro.png
部分腐食③の場合~
全体的に見ても、応力はほとんど生じていない.健全部と腐食部で応力ひずみが異なるため、全体腐食と異なる挙動になったと思われる.~
また、健全部の応力ひずみだが、衝撃力がかかった瞬間,応力ひずみともに正の値が出るようになる.(衝撃がかかるまでは圧縮と引張を行き来している挙動.)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_impactsitacoro.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_impactsitacorokenzen.png
結果を見ると、腐食部分が破断することはなかった.(C1のケーブル応力が0.18MPaと低いため、やはり集中荷重を掛けないとケーブルの応力に大きな値が得られない.)
***衝撃荷重の場合（集中荷重も加えたパターン） [#ff632edf]
衝撃荷重だと腐食部に伸びが生じなかったことが確認できたため、やり方を変えてみた.~
解析時間は130sで、まず100sプレストレスと死荷重を掛け.10s安定させる.その後に集中荷重(15000N)を掛けて,その最中にケーブルが切れる状態(120sで破断.)を想定し、130sまで集中荷重を載荷し、解析を回した.（衝撃荷重+集中荷重）~
結果は全体腐食(128.4sで破断;解析終了時間128.405s),部分腐食①(128.279sで破断;解析終了時間128.287s),部分腐食②(128.347sで破断;解析終了時間128.354s),部分腐食③(128.39sで破断:解析終了時間128.393s)の4ケースをまとめた.~
また、どのケースも130sまで収束エラーで解析が最後まで回ることはなかったため、C1ケーブルが破断した瞬間までの記録を取っている.
・解析終了時の挙動（スケールファクタは全て1倍.部分腐食モデルにおいて,C1ケーブルのピンクの点と点の間が腐食部分を表している。色がついている粒？みたいなのが、元々桁のあった位置.）~
[[集中荷重のみ:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?%E5%8D%83%E4%BB%A3%E5%B2%A1%E3%81%AE%E4%BF%AE%E8%AB%96%E6%97%A5%E8%AA%8C#r79d5e43]]と比較して、どの解析パターンもC2のケーブル破断による衝撃荷重が加わったことで腐食部分の要素が大きく伸びていることが分かる.
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/impactkanimodel_zentaicoro.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/impactkanimodel_uecoro.png~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/impactkanimodel_mannakacoro.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/impactkanimodel_sitacoro.png
・先端変位（右側の写真は,破断した瞬間の様子を拡大したもの.）
全体腐食(最終変位：-415.62mm)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_zentaicoroimpact.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_zentaicoroimpactkakudai.png
部分腐食①(最終変位：-879.807mm)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_uecoroimpact.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_uecoroimpactkakudai.png
部分腐食②(最終変位：-703.065mm)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_mannakacoroimpact.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_mannakacoroimpactkakudai.png
部分腐食③(最終変位：-354.928mm)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_sitacoroimpact.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/DEPL_sitacoroimpactkakudai.png
・C1張力（腐食部）（右側の写真は,破断した瞬間の様子を拡大したもの.）
全体腐食（最大張力：19441.5N）~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_zentaicoroimpact.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_zentaicoroimpactkakudai.png
部分腐食①（最大張力：19438.8N）~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_uecoroimpact.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_uecoroimpactkakudai.png
部分腐食②（最大張力：19351.1N）~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_mannakacoroimpact.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_mannakacoroimpactkakudai.png
部分腐食③（最大張力：19441.3N）~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_sitacoroimpact.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1T_sitacoroimpactkakudai.png
・C1腐食部応力ひずみ（グラフは破断後の応力ひずみ曲線を表している.）
全体腐食~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_zentaicoroimpact.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_zentaicoroimpactkakudai.png
部分腐食①~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_uecoroimpact.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_uecoroimpactkakudai.png
部分腐食②~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_mannakacoroimpact.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_mannakacoroimpactkakudai.png
部分腐食③~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_sitacoroimpact.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/C1SE_sitacoroimpactkakudai.png
○衝撃荷重が加わった際の腐食部の応力ひずみ関係の特徴として、120sでの衝撃荷重が加わった後,応力ひずみともに正の値に変わる変化が見られた.
***まとめ [#x3c721bc]
集中荷重のみの場合、部分腐食で解析を最後まで回そうとすると収束エラーが出てしまう。破断時の挙動と応答値を得ることはできたが、全体腐食との解析時間が異なるため、比較するのは難しいかも。~
ただ、部分腐食のグラフを見ると、挙動はそれぞれ似ている一方、解析終了時の値が異なるので、腐食の位置によって応答値の変化がある可能性も？~
全体腐食と部分腐食での応答値の違いを見る場合,衝撃力と集中荷重を掛けたパターンが解析終了時間もほぼ同じのため、比較しやすい。~
→結果を見る限り、似たような傾向の結果が得られているのではないか？ただ、桁の先端変位が腐食ケースごとに値が変わるのかは不思議。挙動を見たところ、腐食部分の伸びの長さがそれぞれ異なっているので、腐食の長さや位置によって応答値の変動はあるのかもしれない。~
**簡易モデルver.2 [#nb955af1]
ばねを用いてピン接合を行うやり方で斜張橋モデルを作成すると、主塔とケーブルを結ぶばねの節点に自由度絡みのエラーが発生する。~
ケーブルをトラス要素に変えたら解析は回ったが、トラス要素は節点共有した際にピン接合になるため、要素分割する際は橋梁ケーブルをトラス要素で再現できない。~
ケーブルのサグを再現するには要素分割を行う必要がある（ケーブルが切れたあとの様子・変化が見たい）ため、まずは簡単なモデルを作成して解析が回るか確認する。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/zentaikanimodelver2.png
今回用いるモデルは上の通りである。1.5mの塔に、塔の高さ0.5mから桁を1mに伸ばし、ケーブルで吊ったモデル。~
このモデルを4つのステップに分けて解析を行う。(解析自体はシンプルに静的解析を行い、変位を導出している。)
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/ketadanmen.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/cabledanmen.png
桁と塔は高さ10mm,幅50mm,厚さ2mmの薄肉矩形断面、ケーブルは直径8mmの円形断面を適用した。
STEP1：ケーブルなしモデル~
塔と桁のみの状態で、桁の先端に100Nの集中荷重を与えたモデル。塔の先端は拘束なし。まず、このモデルの解析が回るか確認した。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/step1.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/modelstep1factor5.png
色付きが変形図（factor5倍）~
桁の先端はy方向に-70.38mm、塔の先端はx方向に19.22mmたわんだ。
桁の先端に荷重を掛けていることから、塔と桁が倒れるような挙動が得られると思われるため、解析は上手くいっているだろう。
STEP2：ケーブルつけたモデル
塔側のケーブルの結合は剛結、桁との結合はバネ(0.01mm)を用いたピン接合。ケーブルの要素分割数は20で、プレストレスを温度で-100℃与えた。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/step2.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/modelstep2factor20.png
色付きが変形図（factor20倍）~
桁の先端はy方向に1.059mm持ち上がった。~
温度によるプレストレスが入っていることも確認。（設定方法は[[こちら:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?Salome-Meca+%E3%83%93%E3%83%BC%E3%83%A0%E8%A6%81%E7%B4%A0#vac22495]]に記載。）
STEP3：塔とケーブルの接合を変位固定で~
塔側のケーブルの結合は回転自由、それ以外の条件はSTEP2と同様~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/step3.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/modelstep3factor10.png
色付きが変形図（factor10倍）~
桁の先端はy方向に2.313mm持ち上がった。~
STEP4：塔とケーブルの接合にバネを用いる。~
塔側とのケーブルの結合もバネによるピン接合に変更~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/step4.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/chiyooka/modelstep4factor10.png
色付きが変形図（factor10倍）~
桁の先端はy方向に1.249mm持ち上がった。~
*つりあい解析 [#t1c73e2f]
つりあい解析とは、事前に解析した構造解析の結果を、次に解析したいモデルにその結果を入れて解析を回すことを指す。~
ちなみに,似たような解析で,連成解析と呼ばれるものが存在する。（温度解析→構造解析とか異なる解析を繋げて行うこと。）~
学部時代に取り扱った解析だと...①：ゆっくり力を与えて安定させるまでの結果を出して、②：次にケーブルに掛けた張力を徐々に落としていく解析にその結果を読み取って解析を回すことになる。（二段階解析というイメージかな？）~
やり方は[[Salomeメモ:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?Salome%E3%83%A1%E3%83%A2#gac275b5]]と[[Salome-Mecaビーム要素:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?Salome-Meca+%E3%83%93%E3%83%BC%E3%83%A0%E8%A6%81%E7%B4%A0#gac275b5]]に記載した。(書かれている内容はどちらも同じ)
*発表等での質疑応答メモ [#n1155511]
研究室での中間発表や学会など、質問されたことをここにまとめ、再度確認したものをここに載せる。~
大阪見学で頂いたご意見~
・モデルについて、力の伝播のことを考えると桁だけでなく塔もつけた簡易的なモデルで解析すべきでは？~
→それをやるなら、橋梁モデルで解析したほうがよい。そのため、これから橋梁モデルでの解析を行う方針となった。
・今の設計法では、吊り形式橋梁は連鎖崩壊しにくいと意見を頂いたが、設計法とは何を指していたのか。~
→部分係数法が変わっているらしい。この研究では、設計法が変わる前（昔の橋梁）を対象としているため、これから作られる橋梁に対して何か言うわけでなく、既存の（設計法が変わる前）橋梁の連鎖崩壊が起きないように考える、研究となる。
春休みの研究報告~
簡易モデル2本でケーブル全て健全部,全体腐食,部分腐食での数値解析を行う.(取るものは応力-ひずみ,張力(健全部と腐食部),桁先端変位)~
→予想では全て同じ結果になるんじゃないかな？（同じ結果になれば、腐食モデルを全体と部分で異なるのかについて証明できる.）
荷重の掛け方は、集中荷重と等分布荷重とC2のケーブル破断荷重で.~
11/28 理工学デザイン~
・部分腐食させて連鎖崩壊を起こす内容は理解したが、結局それをどう活かすのか。（研究の最終目的は？）~
部分腐食の再現をすることで、ケーブル破断の影響を見ることが目標。（最終的には、衝撃力をあるパラメータとして数値化して入れることで公式的に導出できるようになればいいかな。）
・現状,応力-ひずみ減少で腐食再現ができていない状況だが、他の方法で検討していないのか。
断面積減少という方法はある。ただし、この方法ではただ要素が伸びて終わるだけなので、連鎖崩壊に繋げる意味では適切ではないかも。~
また、短期的（南方澳大橋みたいにトラック通ってから切れるみたいに）な破断をするという意味では「ひずみがある点に達したら切れる」が適切だろうと考えている。~
（発表時にはヤング率を時間で減少させて落とす方法もあるかもと言ったが、あくまでそれは長期的な話なので自分の研究内容とずれる。）~
・連鎖崩壊するには、他のケーブルも腐食してないと起きないのでは？（1本腐食したところで、他のケーブルが健全であれば問題ないよね？）~
そのとおり。この研究は連鎖崩壊のメソッドを理解することにもつながる研究なので、どの腐食部分が破断すると次にどのケーブルが切れるかという法則性？をある程度理解するために行っている。
まとめ~
指導教員からアドバイスを頂いた上でこの項目を編集した。理工学デザイン→振り返りをしたことで、改めて自分のやるべきことが何か理解した。~
この経験を踏まえ、これからの研究を引き続き進めていく。
11/13 ポスター形式中間発表~
その時の資料は,k2→2024→chiyookaに入っている.11/13にアップロードしたので,アップロード日付を見ればすぐ分かると思う.~
・要素分割とは？~
A:ある要素を任意（自動的に）に分割すること。~
・ケーブルの部分腐食の先行研究・文献はないの？（海外とか）~
また、落橋したケースの逆解析はやられてないの？~
※逆解析：すでに分かっている結果をもとに、未知の解析条件を推定すること。~
A:先行研究等に関しては調査中.（調べた範囲では行われていないが,見落としがあるかも。)~
梁要素で部分腐食の論文は調べた範囲では見つかってないが、海外の論文は見落としていたため、発表時は曖昧な答え方をしてしまった...。~
逆解析は行われている。
・全体腐食するのではなく、ポイントごと（要素分割ごと？）で腐食を考慮している？~
A:まだ部分腐食の再現ができていないため、それは考慮していない。~
・応力-ひずみ関係を参考にした実験について詳しく.~
A:外部の方がやられた研究を参照したため、ここでの記載は避けます。~
（研究ノートに残す）
・腐食が起きる原因は？~
A:金属の表面に水と酸素があるとき、化学反応（酸化反応）によって浸食されることで起きる。~
詳しくは[[こちら:https://www.nipponsteel.com/company/publications/quarterly-nssmc/pdf/2018_23_12_15.pdf]]
ちなみに,海側で腐食が起きやすいのは海からの塩分が酸化反応を手助けするためである。~
・何年ぐらいで落橋に至るのか？（完成してから何も対策もしないで）~
A:資料によってバラバラだが、約50年前後。勿論、対策をすればもっと伸びる。~
・端部で腐食させた理由は？~
A:ケーブルの端部は水が貯まりやすい場所でもあるため、腐食が起こりやすい。そのため、まずは腐食が起こりやすい場所で部分腐食の再現をしてみようと思い、端部にした。~
・応力-ひずみ減少は健全,腐食の位置を変えても同じように機能する？（発表時,ケーブル端部にのみ腐食を与えていた。）~
A:要素ごとに応力-ひずみ関係を変えただけなので、機能する。
・ケーブルの腐食だけを検討している理由は？桁とかも考えなくていいのか？~
A:連鎖崩壊の解析に繋げる研究でもあり、その現象はケーブルに着目している。そして、連鎖崩壊はケーブルの腐食の有無が重要であることから、ケーブルのみを見ている。研究が進めば、桁とかを考慮しても面白いかもしれない。~
発表時,どうしてケーブルの腐食を研究しているのかという意図で回答したが...ずれた回答していたかも。
・①ケーブルを切るのではなく、ケーブルを伸びきった状態にするのはどうか？仕事をしていないものにもっていけばそれは同じ意味なのでは？~
②連鎖崩壊に繋げる際、伸びきったケーブルの力を取って、他のケーブルに当てはめればケーブルの荷重分配もできるのでは？~
A:①伸び切ってケーブルそのものの役割を果たさないという考えはありかも。その方法も1つの手法として検討するのもあり？~
②手動ではなく,自動で違うケーブルに力が伝わるようにしたい。
・応力-ひずみ関係の破断（応力0に落とす）させるところまでの直線は何を表している？~
A:ケーブルの伸びを表している。
反省点~
たくさん質問してくれたにも関わらず、自分でも納得のいく回答をしたのが少ないと感じた。~
特に、1発目で想定外の質問が来て、そこから終始焦っていたため、もっと肩の力を抜いて回答できるとよかったかも。~
そして、今回の発表を通じて、自分の研究を理解しているように思っていたが、まだ詰めの甘いところがあると思った。~
特に、自分の研究と他の研究の繋がりの理解が足りてないと感じたため、ポスター発表までそのあたりをしっかり見直そう。~
最後になりますが、質問してくださった皆さん、ありがとうございました。
*参考資料等 [#a2f72b97]
[[衝撃メモ:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?%E8%A1%9D%E6%92%83%E3%83%A1%E3%83%A2]]（DAFとは何なのか,改めて勉強し直す用）
[[code_asterの動画:https://www.youtube.com/@DigvijayPatankar/playlists]] （音声は英語だが、内容は振動解析に関わってくるもの。）
[[code-asterの説明:https://code-aster.org/doc/v15/en/index.php?man=commande]] （英語で書かれているので翻訳する必要あり）
[[code-aster質問ページ:https://forum.code-aster.org/public/]]
[[beam要素の弾塑性解析マニュアル:https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/tamura/beam_nonliner.pdf]]
[[Plasticity tutorial rev.1.1Code Aster® STA10.3:http://www.caelinux.org/wiki/images/3/39/Claws_Plasticity.pdf]]（英語だけど、弾塑性絡みの設定について書かれている。）
[[CAE用語辞典:https://www.cybernet.co.jp/ansys/learning/glossary/]](解析の種類などの意味を調べるのに使う)
[[Marc Mentat2018メモ:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?Marc/Mentat2018%E3%83%A1%E3%83%A2]]
[[例題で学ぶMarc有限要素法解析入門:https://www.lib.akita-u.ac.jp/opac/volume/815629?current=1&q=Marc%E6%9C%89%E9%99%90%E8%A6%81%E7%B4%A0&total=2&trans_url=%2Fopac%2Fsearch%3Fcount%3D10%26defaultpage%3D1%26defaulttarget%3Dlocal%26order%3Drecommended_d%26q%3DMarc%25E6%259C%2589%25E9%2599%2590%25E8%25A6%2581%25E7%25B4%25A0%26searchmode%3Dnormal]] (2013年に出版されたものなので、少し違いがあるかもしれない。なお、閲覧にはAUアカウントとPW必要)
[[Salome-Mecaのチュートリアル:https://code-aster.org/V2/spip.php?rubrique67]] (英語だが、Salome-Meca全体の話が色々載っている。例えば、作成した構造物に適した要素(3D,BAR...etc)とか)
[[Pavaris wiki:https://www-paraview-org.translate.goog/Wiki/ParaView?_x_tr_sl=en&_x_tr_tl=ja&_x_tr_hl=ja&_x_tr_pto=sc#Analyzing_Data]]
[[Salome-meca beginning with code_aster:https://archives.framabook.org/docs/Code_Aster/beginning_with_code_aster-jp_aubry-20190129.pdf]] (Salome-mecaの質問するところ[掲示板]でよくこれを見てくれと答える人が多いので、載せてみた。英語で書かれている。)
[[Salome-Meca 2019を用いたトルク解析について:https://qiita.com/Jun_Tatsuno/items/3bc1d5d30e9055fbe388]] 0D要素の作成方法が書かれている。~
[[Midas解析編:https://in.midasit.co.jp/wp-content/uploads/2020/07/constrcution-83-civil-book.pdf]] 今年からこの解析ソフトを使い始めている方がいるので,自分も軽く頭に入れておきたい.
*解析ソフトの使い方 [#q171dd69]
Salome-Mecaビーム要素については[[こちらから:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?Salome-Meca+%E3%83%93%E3%83%BC%E3%83%A0%E8%A6%81%E7%B4%A0]]
Marcに関することは[[こちらから:https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?cmd=read&page=Marc_Mentat%E3%83%A1%E3%83%A2&word=Marc]]
**midasについて [#vf9f6c8f]
修論落ち着いたら,midasでmodel300でも作ってみようかな？
***解析の種類 [#r525f7e5]
・静的増分解析~
Salomeでいう,MECA_STATIC?
・固有値解析（モード解析）~
Salomeでいう,CALC_MODE.
.応答スペクトル解析~
・時刻歴応答解析~
→そこから,線形時刻歴応答解析,動的フレーム非線形解析,境界非線形解析もできるとのこと.
Salomeでいう,DYNA_NON_LINEと同じ解析ができそう.
・P-デルタ解析~
p-デルタ効果と呼ばれる,「水平荷重と鉛直荷重の相互効果」を考慮した解析とのこと.
・座屈解析
・幾何学的非線形解析~
Salomeだと,STAT_NON_LINE（変形をGROT_GDEPにした場合）.
・静的材料非線形解析~
Salomeだと,STAT_NON_LINE.
弾塑性解析はこれを使うのかな？
・施工段階解析~
Salomeだと聞いたことがない解析.~
（定義）施工途中の荷重や部材応力などの状態変化を考慮して、遷移する事象ステップごとにモデルの形状や状態を一部変更しながら応力や変形を引き継いでいく段階的な計算方法.
*橋梁写真 [#dd55c1f6]
自分で撮った橋梁の写真をここに載せます.（不定期更新）~
(1)瀬戸大橋 (写真手前側が北備讃瀬戸大橋,奥側が南備讃瀬戸大橋)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/setobridge.png
撮影場所：香川県坂出市 与島PA~
高速道路の休憩所から撮影。四国に行くときは毎回瀬戸大橋を渡るため,ここのPAに寄るたびに撮影してます(^^)~
瀬戸大橋は道路と電車が走れる構造（鉄道道路併用橋）になっており、この構造での橋梁の長さとしては世界一だそうです。~
(2)明石海峡大橋~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/akashibridge.png
撮影場所：兵庫県淡路市 淡路SA~
日本で最も長い吊り橋です。（世界で見ると2番目です。世界1位は、確かトルコの1915チャナッカレ橋だったかな。）~
夜になるとライトアップするので、高速のサービスエリアから神戸の夜景と一緒に見るのが個人的にオススメ。（ライトアップの種類は、季節等によって変わってたはず...）~
また、ここのサービスエリアには観覧車があるので、それに乗って景色をみるのもありかも。~
(3)由利橋~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/yuribridge.png
撮影場所：秋田県由利本庄市~
秋田県にある斜張橋です。~
後輩たちがこの橋梁を研究対象に活動を頑張っているのを見て、自分も見に行ってきました！
朝に行って来たけど、交通量もそれなりにあって、散歩で使っている人もちらほら見かけて、周辺に住む人達の生活にとって大事な橋だなと実感しました。~
ちなみに...昔はトラス橋だったらしく、架替え？にあたって今の斜張橋になったらしいです。
(4)東京レインボーブリッジ~
http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/rainbowbridge.png
撮影場所：東京都港区
東京にある有名な吊橋です。~
東京で学会があったので、ついでに行ってきました！~
歩いて渡ることもできるみたいですが、今回は時間がなかったため断念...。~
次来たときは、歩いて渡りたいな(^^)
*日記 [#m48ef041]
12/4(水)~
今年も残り1か月,月日が経つのが早いと感じます。~
研究成果をそろそろまとめ始めないと,修論までに間に合わなくなりそうだし,計画的に動かないとな...。~
12月は,同期との旅行や中間発表とイベントがたくさんあるし,やるときはやる,楽しむときは楽しまないと！
11/9(日)~
気温も下がり,冬が近づいてきたと近頃感じます。~
後期も色々とやることが多く、今年度は学生生活の中でも最も充実な1年になりそうです。~
ただ,最近は自分自身のことで悩むことも増えてきました。（あんま考えすぎもいけないと思うけど...）~
何か人として足りていない部分があるんだろうな...ほんと難しいな
10/29(水)~
中間発表も終わり,3年生も研究室に配属され,10月もいよいよ終わりに近づいてきました.~
今年も残り2ヶ月で終わりで,1年あっという間だなとも感じます^^;~
学生生活も残り約半年で遂に終わりを迎えますし,学生時代にできることは今のうちにやっておかないとな...！
9/3(水)~
先日,秋田に戻ってきました.~
夏休み(8月)は色々な場所に出かけたり,社会人になった友達と会ったりと充実した夏休みを過ごせました.~
今日から研究室生活が再び始まりますが,残り半年頑張っていきたいです.
8/7(木)~
前期がようやく終わりました.~
今年の前期は,本当に色々あったな...~
何なら,学生時代の春から夏にかけてにおいて，一番忙しかったかもしれない(-_-;)~
帰ったら,まずは春に行けなかった自動車の本免試験を受けに行かないとな.まだまだ気は抜けない！~
終わったら,前期で遊べなかった分,息抜きしないとな！（勿論，やることはやった上で）
7/17(木)~
色々やることがあって大忙しだけど...学ぶことも多い日々を過ごしています.~
何というか...物事を表面だけ理解していて,深くまでは全く理解できていないなと感じるな.~
ほんと,まだまだ甘い部分があるなって痛感させられます。~
もっとしっかり勉強しないとな。
6/30(月)~
明日から7月,前期も残り1か月少し。早いな...研究が思っていた以上に進んでなくてだいぶ焦っています💦~
ただ、梁要素についてのSalomeに関する謎が1つずつ解決できてきて、橋梁解析がようやくできそうってところまで来ているから,もうすぐ色々なことがやれるようになるはず！~
研究以外にもやることが多く大変だけど...1つずつ丁寧に片づけていかないとな。~
しかし...修士2年になって3ヵ月、卒業された先輩方みたいに後輩を引っ張っていける人間に、少しは近づけているのだろうか。（深く考えすぎもいけないと思うけど...）~
ほんと、しっかりしないとな。
5/24(土)~
学会終わって1周間弱,昨日この前と2日間,息抜きをしました。~
学会終わって研究やるかと思っていましたが,あまり取り組むことができませんでした...流石に少し疲れていたか^^;~
おかげで,気分もリフレッシュできたので,ここから気持ちを切り替えて取り組みたいと思います。~
まずは,夏休み前までにある程度の成果を出さないとな！頑張ろう!!
5/21(水)~
英語の学会が終わり、普段どおりの生活が戻ってきました。~
学会に参加した感想ですが,英語力が全く足りないなと実感しました。~
発表以外での英語の会話はあまり喋れなかったし,質問も何言ってたか分からなかったしな...。~
将来,必要になる場面がまたあると思うので、しっかり勉強しなきゃなと改めて思いました。（昔,中学の英語の先生から言われた「英語からは絶対に逃げられない」は嘘じゃなかったな^^;）
4/27(日)~
かなり久しぶりの更新です^^;~
修士2年になり、そして最後の学生生活が始まりました。~
冬から春休みにかけて、学会の概要作成や就職活動などもあり、そして最近は英語の学会に向けての準備や大阪での工場見学もあり、特に今月は色々とイベントが多い月でした。~
今年度は修論もあって、昨年度より忙しい1年になりそうですが、楽しい1年にしたいと思います(^^).
...博士の先輩が卒業されたから、今年度は修士2年が最高学年なんだよな(・・;)~
しっかりしないとな。
12/10(火)~
自分の関わりのある今年の大きめのイベントがすべて終わりました。~
帰省までの残りの時間を研究活動に費やそうと考えています。研究を進めるスピード遅かったからな...~
冬休みの予定も立てたし、それを楽しみにしながら2024年残りの時間を楽しく過ごしたいと思います!~
10/21(月)~
また久々の更新です。~
10月も後半、この時期になるともうすぐ1年が終わると感じるようになります。~
さて,先週は3年生が研究室に配属され、顔合わせを行なったそうです。私は授業があったため、まだ顔を1度も見れていませんが、今度の飲み会でどんな子が来たのか見れるのが楽しみです。~
しかし...自分が研究室に配属されて、もう2年か(・・;) この2年でちゃんと成長できているのだろうか。~
ここでの生活もあと約1年半、悔いのないように過ごしていきたいな。
7/20（土）~
久々の更新です。~
前期も残り3週間ほどで終わります。早いな...。~
来週はある授業のスライド発表に、再来週はある授業のテストがあったりで大変ですが、忙しい日々も夏休みに入れば一旦落ち着くと思います(^^;~
ただ、授業のほうもある程度落ち着いたら研究のほうに力を入れなきゃなとも思っています。昨年度より研究の進捗ペース落ちてるからな...。~
夏休みの楽しみも頭の片隅に入れながら、残りの前期も頑張って乗り切りたいと思います！
5/15(水)~
今日はいよいよ4年生の研究テーマが決まりました🎉~
今年の研究テーマで自分が何かしら大きく関わることはなさそうですが、自分も卒論の際に先輩にたくさん助けていただいたので、もし何か聞かれたら少しでも力になれるようにはなりたいなと思います。~
自分の研究も最近はスローペースでしたが、院での生活も慣れてきましたし、ペースアップしていきます！！
5/9（木）~
GWも明けて、ようやく前期も本番かなと感じます。（そう思うのは私だけですかね💦）~
昨日はM2以上の院生と先生方が不在だったため、自分が前に立ってゼミを進行しました。~
4年生がこれから行う研究内容について分かる範囲で簡単に説明することが主な内容でしたが、あの説明で伝わったのだろうか...(・・;)~
毎回思うことだけど、日本語って難しいよな...。大学に入学した頃と比べると、間違いなくコミュニケーション能力は上がっているけど、まだまだ改善するところはあるよな💦~
まぁ...それは置いといて、最近は人前に立って話すことが増えたなと感じています。~
今年は、ゼミを除くと卒論,東北支部,謝恩式と大きめな行事で人前で話していますし...中高の時の自分がこれを知ったら驚くだろうな(笑)~
逆に、人に分かりやすく伝えるにはどうすべきかと考えることも増えましたが、そのおかげで少しずつだけどコミュニケーション能力がついてきているのかなとも思います。~
院の授業でも議論・発表がメインになってくるので、この1年でかなり鍛えられるんじゃないかな？もっと頑張らないとな(^^)
4/24(水)~
もうすぐGWが近づいているのと同時に4月がもうすぐ終わってしまうことに驚きを感じています...。~
のんびり過ごしてたらあっという間に夏ですね(・・;)~
それと、本日も英語ゼミに4年生と混じって参加しましたが、自分の研究を英語で説明するって難しいなと思いました。~
日本語でも中々難しい言葉を使っているので尚更ですね💦~
間違いなく今の4年生に伝わってなかったので、ちょっと考えないといけないな...。
4/18(木)~
研究指導計画書も書き、今年度をどう過ごすかある程度考えることができました。~
前期は結果を出すことに集中して、後期からは論文を書けるようにできると理想ですが…果たしてどうなるか💦~
けど…後輩たちの頑張ってる姿を見ると先輩がだらしないのはいけないなとも思います。~
すでに夜遅くまで残って研究してる人が何人かいるし…(゜o゜;~
うかうかしてたらあっという間に追いつかれそうです(^^;（競争してる訳ではないですが…）~
もっと頑張らないとな！
4/16(火)~
新学期が始まって1週間程経ちました。明日から本格的にゼミが始まるので,そろそろペースを上げて頑張らないとなと思っています。（一応,明日話すことはあるので,その心配はしなくてよいですが...）~
さて...明日は新4年生の春課題の発表日です。昨年は自分も発表しましたが、あれから1年経ったと思うと色々思うことはありますね(^^;~
めちゃくちゃ緊張してたし、質問の受け答えに対しても意味わからんこと言ってたよな…(-""-;)
昨年の自分を思い出しながら,明日の発表は聞こうかなと思います😊
4/12(金)~
第1回目の授業がほとんど終わり、どれも面白いと思える授業ばかりでした。~
後は、自分が授業を通じてどれだけ自分のものにできるかが大事なので、しっかり勉強していきたいと思います。~
研究室絡みでは、新4年生が来週に春課題の発表があるため、研究室に多くの人が集まっています。~
どんな発表をするのか楽しみです！~
それと...もっと新4年生とも交流しないとな。外国語文献の講義は院生も出れるから、まずはそこで色々お話してみるか！英語だけど...笑
4/10(水)~
遅くなりましたが、ついに大学院生としての生活が始まりました。~
授業に研究と、ひょっとしたら昨年度より忙しい日々になりそうですが、1日1日を大事に過ごしていきたいなと思います.~
今年度はどんな1年になるかとても楽しみです.~
現在の閲覧人数: &online;~
総閲覧数:&counter;~
ページ名:
![[PukiWiki]](image/pukiwiki.png "[PukiWiki]")
