![[PukiWiki]](image/pukiwiki.png "[PukiWiki]")
#author("2026-08-20T16:42:14+09:00","default:kouzouken","kouzouken") #author("2026-08-20T16:56:58+09:00","default:kouzouken","kouzouken") [[FrontPage]] ここに修論日誌を書いていきます 2025年度は-[[湊の卒論日誌]]に #contents *やること [#uf937f4d] %%-活荷重満載で集中荷重600kNを移動させた際の影響線%% %%-活荷重満載で集中荷重600kNを270mにかけた状態で600kNを移動させてできた影響線%% %%--活荷重満載→LP1、LP2両方2車線に載荷%% -衝撃係数 --斜張橋の自動車走行による動的応答と衝撃係数に関する研究 -弾塑性域を超えるまで載荷して影響線 -多々羅(にするかは不明)での影響線 -モデル600,900で同様の載荷 -maidasで同様の載荷 2026/6～7 目的：Model300のC1,C14の影響線を求め、重ね合わせの適用可能性（非線形解析の必要性）を検討する。 Step 1: 線形解析を用いて、Model300のC1,C14の影響線を求める。 F=200kNを主桁の節点上に移動させる。D,CW,Prは不要。 Midas ファイルはこれまでと名前を変更する。（Ex: model300 IL） Step 2: 同様の計算を、model300にて非線形解析を行い、C1,C14の影響線を求める。 Step1 の値と比較する。両者は差がないと予想される。 Step3: Model300にて、非線形解析にて D+CW+PR+L2 を載荷させたときのC1,C14を求める。 Step 4: Model300にて、非線形解析にて D+CW+PR+L2 に加え、F=200kN を節点上に移動させ、C1,C14の仮の影響線を求める。影響線は、F=200kN のみの値である。 Step 5: Step 4で求めたC1,C14の仮の影響線から、Step 3 で得られたC1,C14の張力を引けば、C1,C14の影響線が求まる。これをStep 2の結果と比較し、影響線に関して重ね合わせが適用できる、ことが検証できる。 上記重ね合わせの原理の確認終了 *やりたいこと [#uf937f4d] %%-Paravis での作業簡略化%% %%55個のファイルを同時に開いて、同時に同じ作業をしたい。%% %%マクロ的な項目はあるのでできそうな気がする%% *作業時気をつけたいこと [#bb0326bc] -解析結果と条件、medファイルの管理 --エクセルにメモだけじゃなく解析結果ファイルの中に、med、commファイルもおいておいてなぜその結果になったかすぐ確認できるようにしておく -作業ファイルについて --作業するファイルはまずその日の日付を先頭に --解析内容をファイル名につける --完成品は完成品で別にまとめておく *エラー不具合、修正箇所など [#u2b6cf89] -salome,AFFE_CHAR_MECAでLP2とDLの荷重設定が逆になっていた -今までの死荷重ありでの解析にLP2が考慮されていた ---LP2=(3.0×5.5)+{1.5×(18-5.5)}=35.25kN/mの2車線70.50kN/m DL-244000.0N/mに比べると非常に小さいため形状への差は小さいはずだが解析やり直しをする -別名保存について --Salome-mecaで別名保存を行うと、input,outputを大量に設定している場合別名保存先と元データ両方においてデータが破損する可能性がある。 今回はinput13個output55個で行ったところエラーが起きた。 input,outputの参照元の場所がファイル名の変更によって変わったためだと考えている。 対策としてはデータを移行する前に、一度解析を回しておいてcommファイルを抑えておくことやinputデータの場所把握、負荷低減のためsalome-meca以外を立ち上げておかないことが挙げられる。 --また、inputデータの編集画面においてEmbeddedに設定しておくと容量は増えるが解析ファイルそのものにmedファイルをコピーしておいておけるみたいなので、medファイルを保存した上でk2にバックアップを取っておくのが良いかもしれない --破損した場合は全部設定し直したほうがひとまず早そう --触ってみたところ不明なエラーが出てしまったので変更しないほうが良さそう *jwcadからmidasに [#ta724404] https://support.midasuser.com/hc/ja/article_attachments/40001043095833 SPCを利用した任意形状の断面特性計算-基本編 を参考に cadで図面を作成 X DF形式で保存 図面作成時，厚さを考慮してもいいが上手くいかないときもある SPCで厚さを与えられる maidas Apps 断面性能計算(SPC) **断面性能計算(SPC) [#u26a4bdb] - 単位系の選択 -- fail import AutoCAD_XDF から図面を読みこみ -- 断面厚さ定義 -- model curve Change_Width -- Width に厚さを入力 -- Select 押してから厚さを設定したい部材を選択(Ctrl押しながらでも) -- Apply **閉合断面部分(鋼ではない空いている空間)の選択 [#u6245e82] -- model curve closed_loop register -- Select リブ,空間を囲んでいる部材を選択 -- Apply **断面の定義 [#rec89645] -- model section generate -- Type Line -- Name "名称" -- select で全体を選択 -- Apply **断面計算 [#j2ab2763] -- property calculate_section_property -- Fineを調整(デフォルト問題ない) -- select で全体を選択 -- Apply Property List_Section_Propertyで断面特性を確認 **保存 [#j18c9a11] -- model section Export -- MIDAS section_file を選択 -- ...から保存クリック -- select で全体を選択 -- Apply **MIDAS [#s5fd8130] --材料 断面 追加 値入力 任意形状断面 ファイルを読み込み ok XDFファイルが読み込まれない時，描写されないときは,何回か再読み込みmidasをもう一回開くなどする SPCでXDFファイルを開くと，cadで作ったもの全てを一画面に収める状態から始まるので離れた別の箇所にもう一つある等の場合めっちゃ小さく表示されてしまい，見えなくなる midas側で作った断面の向きを変えられないのでSPC側で回転させてファイル読み込みさせる rotation selectで全体選択 角度入力 apply *モデル、解析について [#d81327e0] ***model300について [#jbeb5010] -cable:deadあり -POUTRE:構造物に対して載荷 今は主塔のA部分の死荷重に使ってる 同じ解析内で一回しか使えない -NODALE:線に載荷 載荷したいNを長さで割って載荷 その方向にかかる ***載荷条件[#wc2591c1] LP1 = =_F(FY=-1175000.0,GROUP_NO=('saika', )), tousaika = =_F(FY=-60000.0,GROUP_MA=('syutou_huto', 'syutou')), CW = =_F(FY=-160000.0,GROUP_NO=('CW', )), LP2 = =_F(FY=-70500.0,GROUP_NO=('keta', )), DL = (FY=-244000.0,GROUP_NO=('keta', )), ***Paravisの見方 [#e7ff373e] -静的解析、影響線作成 med→Alphabetical→ELNOfieldToPointGaussianで軸力を取得 Spread Sheetから見たい点を探せる。 Export Spread SheetからExcelファイルを取得 -動的解析、衝撃力解析について ELNOfieldToPointGaussianの後取得したい点を選択→Alphabetical→Plot Section over time 適用するとグラフが描画 Plot Section over timeをSpread Sheetで選択するとExcelファイルを取得 -軸力が0になっている時http://www.str.ce.akita-u.ac.jp/~gotouhan/j2026/minato/gazou/CW60_1PRE.png 表示→ツール→Current time controlsにチェックが入っているか確認 チェックを入れた後、解析をみたい時間を選択 チェックを入れた後、解析をみたい時間を選択 ***midas トリリニア 応力ひずみ関係 [#o61556b7] -材料断面 -ファイバー要素 ファイバー要素の材料定義 -追加 -名称指定 -材料タイプ→鉄骨 に変更 -履歴ループタイプ トリリニア，Parkモデルなど使いたいやつ -参照材料→ケーブル -各部降伏強度など入力 このパターンだとファイバー要素で作成することになるため一旦なし *成果グラフ、モデル画像 [#zfaebf1b] ***活荷重満載C1ケーブル要素 [#paba7121] http://www.str.ce.akita-u.ac.jp/~gotouhan/j2026/minato/gazou/81_c1_max.png ***活荷重満載C13ケーブル要素 [#q02246c5] http://www.str.ce.akita-u.ac.jp/~gotouhan/j2026/minato/gazou/81_c13_max.png ***活荷重満載C1梁要素 [#hee83573] http://www.str.ce.akita-u.ac.jp/~gotouhan/j2026/minato/gazou/71_c1_max.png ***活荷重満載C13梁要素 [#e24db66d] http://www.str.ce.akita-u.ac.jp/~gotouhan/j2026/minato/gazou/71_c13_max.png ***活荷重満載+270m地点に追加の載荷600kNC1ケーブル要素 [#t4ecba35] http://www.str.ce.akita-u.ac.jp/~gotouhan/j2026/minato/gazou/81_c1_270.png ***活荷重満載+270m地点に追加の載荷600kNC13ケーブル要素 [#y50de753] http://www.str.ce.akita-u.ac.jp/~gotouhan/j2026/minato/gazou/81_c13_270.png ***活荷重満載+270m地点に追加の載荷600kNC1梁要素 [#r257d6a2] http://www.str.ce.akita-u.ac.jp/~gotouhan/j2026/minato/gazou/71_c1_270.png ***活荷重満載+270m地点に追加の載荷600kNC13梁要素 [#h9d15927] ***http://www.str.ce.akita-u.ac.jp/~gotouhan/j2026/minato/gazou/71_c13_270.png [#j65247ec] ***model300 L1+D+CW+PRE 変位図 [#w4d62683] http://www.str.ce.akita-u.ac.jp/~gotouhan/j2026/minato/gazou/L1_300.png ***model300 L2+D+CW+PRE 変位図 [#m7356295] http://www.str.ce.akita-u.ac.jp/~gotouhan/j2026/minato/gazou/L2_300.png ***model300 L3+D+CW+PRE 変位図 [#ccd17124] http://www.str.ce.akita-u.ac.jp/~gotouhan/j2026/minato/gazou/L3_300.png ***model300 [#tb53109f] http://www.str.ce.akita-u.ac.jp/~gotouhan/j2026/minato/gazou/300.png http://www.str.ce.akita-u.ac.jp/~gotouhan/j2026/minato/gazou/300_2.png ***model600 [#g94015f1] http://www.str.ce.akita-u.ac.jp/~gotouhan/j2026/minato/gazou/600.png http://www.str.ce.akita-u.ac.jp/~gotouhan/j2026/minato/gazou/600_2.png ***model900 [#yca2add7] http://www.str.ce.akita-u.ac.jp/~gotouhan/j2026/minato/gazou/900.png http://www.str.ce.akita-u.ac.jp/~gotouhan/j2026/minato/gazou/900_2.png ***CW60kN PRE-150 [#g448968f] http://www.str.ce.akita-u.ac.jp/~gotouhan/j2026/minato/gazou/CW60PRE.png http://www.str.ce.akita-u.ac.jp/~gotouhan/j2026/minato/gazou/CW60PRE_tou.png ***CW90kN PRE-150 [#m5350da8] http://www.str.ce.akita-u.ac.jp/~gotouhan/j2026/minato/gazou/CW90PRE.png http://www.str.ce.akita-u.ac.jp/~gotouhan/j2026/minato/gazou/CW90PRE_tou.png ***CW90kN PRE-150 [#q3ff8cd8] http://www.str.ce.akita-u.ac.jp/~gotouhan/j2026/minato/gazou/CW120PRE.png http://www.str.ce.akita-u.ac.jp/~gotouhan/j2026/minato/gazou/CW120PRE_tou.png ***主塔分割数 [#r9dc981c] http://www.str.ce.akita-u.ac.jp/~gotouhan/j2026/minato/gazou/CW60_1PRE.png http://www.str.ce.akita-u.ac.jp/~gotouhan/j2026/minato/gazou/CW60_2PRE.png http://www.str.ce.akita-u.ac.jp/~gotouhan/j2026/minato/gazou/CW60PRE.png *影響線 [#l7ab1179] http://www.str.ce.akita-u.ac.jp/~gotouhan/j2026/minato/gazou/eikyousen.png 重ね合わせの原理の確認 http://www.str.ce.akita-u.ac.jp/~gotouhan/j2026/minato/gazou/kasane.png *結果の感想 [#f10eaabd] **C13のケーブル軸力が最大になる載荷点について [#ved9e6c7] 350m付近で載荷した際最大値となっており、中心270mと右側主塔420mとの中間(345m)付近である。 ケーブル軸力の振る舞いが対称的なものであるならば190mor200m(主塔120m+(270-120)/2=195m)に載荷した場合C14でのケーブル軸力が最大になる可能性がある *日記 [#vcc5ff52] ***8/20 [#hbca45e8] 中村先生とのミーティング 断面照査について再度行うべきところについて確認を行った． 1各材料の特性値について -主塔→SM490ではなくSM490Y -ケーブル→ST1570の特性値は1160kN/mm2 圧縮が生じる部分については圧縮部材の照査と有効長座屈長について確認 CWの設定では，不反力がないこと，ケーブルに圧縮が出ていないことがメイン その中で，主塔がまっすぐになっているものを選ぶ 主塔の先端変位をCW事にまとめてグラフにしておく 主塔の断面照査には有効座屈長を考慮 主桁は有効幅を考慮する 示方書の13-3-4に書いてあるとのこと せん断遅れについても確認 ケーブルはすべてのケーブルについて行う 全てのケーブルに対し，活荷重込み載荷で降伏を超えないこと，ゆるみがないことを満たすようにPREを導入する 安全率については0.7以上を目安に低すぎないよう気を付ける 安全率が低すぎると経済的ではない設計となってしまう． CW，PREを決定した後，エルンストの式を使ってヤング係数を比較する←また後でやる ***8/17 [#v4ff638e] ファイバー要素の設定を進めた． D+CW+PREの後L1を載荷するよう設定したつもりがうまく設定できていないようで変位が小さくなっている． 改めて設定を見直していきたい． ***8/3 [#w6028366] CWを変更していった時の各断面に発生する応力についてまとめた． まとめた断面照査をもとにCWやPREについての設定を考えていきたい． ***7/30 [#n75c37bf] model300 CWを変更させた荷重ケースについて設定を行った． 主塔変形に関しての要素分割を行った． 以前荷重をかけた際，主塔がくの字に変形してしまっていたが，二等分，三等分させたところ主塔頂点変位は変化しなかったものの，実際の変形挙動に近いような変形を示していた． 要素分割は，主桁部とケーブル接合部最下端の間に対して行った． ***7/27 [#m2424450] 中村先生とミーテイング 断面照査，ファイバー要素の設定について伺った． ひとまず，CW，PREの調整を進める． ***7/14 [#o973cc83] 初期断面力を与えられません プレテンション，初期断面力， 幾何剛性用の初期軸力 と温度荷重がmidascivilで使える軸力を与える方式だと調べた限りでは感じた 材料非線形解析時は初期断面力， 幾何剛性用の初期軸力，温度荷重が使えず，梁要素にはプレテンションは使用できなかった． 温度によるプレストレスの載荷時に生じるケーブル軸力について，ケーブル軸力をXYZ方向に分解したものをケーブル両端部に節点荷重として載荷を行った． しかし，主塔，主桁に載荷した力が分散するせいなのかわからないが想定していたような圧縮力が生じなかった． ***7/10 [#o45a92ca] 降伏強度の設定いったん終了 幾何学非線形，材料非線形考慮時，軸力-荷重関係に降伏時のような挙動を示していた 幾何学非線形を外した時でも問題なく解析が行えるか後ほど確認する 初期張力について 温度荷重と材料非線形は併用できないとのエラー 初期断面力 https://support.midasuser.com/hc/ja/articles/17953719699353-%E5%88%9D%E6%9C%9F%E6%96%AD%E9%9D%A2%E5%8A%9B を用いて設定する予定 ただし，今までの温度荷重による生じたケーブル軸力と弾性係数を用いて出した軸力を与えて解析したケーブル軸力 で値が一致しなかった． 設定を確認するか，温度荷重によって生じた軸力をそのまま初期断面力といて入れて解析するか考える． ***7/7 [#rd6edae8] 青木先生のmodel300影響線の発表終わり 終局強度に関して進めていく予定 ***6/29 [#z2a712aa] model300の影響線完成 幾何学非線形有無の差と幾何学的非線形考慮時の重ね合わせの原理についての確認を行った。 予想通り常時荷重の範囲では差が生じない、重ね合わせの原理適用可能と言えそう ***6/23 [#nd55e356] モデル作成終了 300,600,900のL1+D+CW+Preが中村先生の論文と近しい値を示すことを確認できた フィンランドの橋の解析を今度やるからモデルを作る必要があるらしい ***6/3 [#qf749fb0] 中村先生とミーティング 主桁の載荷ミスあり →修正 主塔、死荷重のかけ方を変更 →midasの自重機能ではなく、等分布荷重での載荷に変更 論文で得られた軸力を得られるか確認 ***6/2 [#ob49936a] midas 二面吊り斜張橋 model300の作成 http://www.str.ce.akita-u.ac.jp/~gotouhan/j2026/minato/gazou/D変位.png http://www.str.ce.akita-u.ac.jp/~gotouhan/j2026/minato/gazou/DCWPre変位.png http://www.str.ce.akita-u.ac.jp/~gotouhan/j2026/minato/gazou/L1変位.png ***5/28 [#g5a57a2c] 中村先生とミーティング 解析を行う条件、主塔の座標やダミー部材の設計について聞くことができた 明日からmidasでの二面吊り斜張橋の作成を行っていく ***5/7 [#pae1c571] 載荷量を上げてC1の降伏応力付近になる載荷を確認した 梁要素 270m載荷時→34,200,000Nで軸力2.784*10^7N 230m載荷時→28,215,000Nで軸力2.73268*10^7N 載荷量と生じる軸力は比例関係にありそう ***4/30 [#r033feda] salome-mecaで載荷荷重を増やして生じるケーブル軸力について解析を行った。 荷重を1000倍にして270ｍ地点に載荷を行ったところC1の軸力が設定していた応力ひずみ関係を離れた値を示すようになった。 salome-mecaでは塑性域に応力が入り、ケーブルが伸びる挙動になったとき、境界線エラーを起こしてしまうので、他のケーブルとの整合性を取るために、そのケーブルの軸力が大きくなっていくのだろうとのことでした。 30ｍ地点に同様の載荷をしたところ、圧縮を示した上、ケーブル要素と梁要素で生じる軸力に変化が生じた。 値については明日以降改めてまとめる。 ***4/13 [#e9a582fa] 引っ越し後解析結果をまとめた。 ひとまず、同じ条件での軸力の結果があっていたので問題なくファイル移行の問題は解決したものとして進めていく。 主桁全域に加え中央スパン10ｍについて活荷重をかけた解析、さらに追加で主桁中央に載荷を設けた条件に関する解析が終わった。 Paravisの数字が集まり次第グラフにしてwikiに載せる。 ***4/10 [#e9a582fa] PCの引っ越し終わり。 無事解析は回ったので次回は今までで解析した時に得られていた値と一致するか確認する。 ***3/25 [#h95b9885] 日誌の作成とデータの整理
