![[PukiWiki]](image/pukiwiki.png "[PukiWiki]")
#author("2024-12-26T15:32:44+09:00","default:kouzouken","kouzouken") #author("2025-01-15T13:04:20+09:00","default:kouzouken","kouzouken") #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]]はこちら 研究タイトル：ケーブル腐食を考慮した斜張橋の連鎖崩壊時におけるDAFの検討 *研究目的 [#ya05820f] 4年次ではSalome-Mecaを用い,簡易モデルで衝撃力の再現が妥当か検討を行なった.ただ,モデルを作成するにあたり,ピン接合の再現や幾何学的非線形絡み（まだ解決してないところもあるが...）の問題があり,それらをできるようにした.~ 修士では,これまでMarc_Mentatを用いて橋モデルを作成して解析を行なっていたものをSalome-Mecaでもできるようにする.そのためにまずはSalome-Mecaで橋モデルを作成して,解析できるところまで持っていく.（これが第1ステップ!）~ そして,ケーブルに腐食を考慮した斜張橋に動的応答解析を行い,衝撃力の大きさを検討していく.（簡単に言うと,昨年,簡易モデルでやったことを橋モデルで行う）~ 斜張橋モデルに関しては,ケーブルの張り方や本数,橋梁のスパン長や減衰率等を色々変えた上で解析を行う. →詳細記述~ ケーブルを有する橋梁のケーブルが腐食することで橋が落ちてしまう事故が近年各地で起きています.（例えば...台湾の南方澳橋とか.卒業された先輩がその橋を対象とした研究を行なっていました.）~ その中で,1本目のケーブルが切れることで,次々とケーブルが切れてやがて落橋に至る連鎖崩壊をターゲットに研究を行なっています.~ この現象は橋梁の落橋などを引き起こす重大な事故に繋がるため,そのメカニズムの研究も進められている現状です.~ 本研究室ではその解析手法の提案が行われてきたが、ケーブルが切れた際の衝撃力の大きさや他のケーブルに掛かる荷重の再分配率が明確ではありません...。~ 卒論では,簡易モデルと呼ばれる簡単なモデルで衝撃力の検討を行ないました.そこで,修論では斜張橋モデルを使って手法の妥当性を検討する目的で,健全と腐食ケーブルを有する斜張橋に対してケーブルが切れた際の衝撃を衝撃力として数値解析的に与え,その大きさを検討していきます.~ モデル自体もケーブルの張り方やスパン長等の性質・値を変化させていくことでどう変わっていくのかも見ていきたいなと考えています.~ 後は,Salomeのビーム要素絡みの問題も色々解決していきたい。 *今後の予定 [#bd1b97fd] ,曜日,12,34,56,78,9 10 ,時間帯,8:50-10:20,10:30-12:00,12:50-14:20,14:30-16:00,16:10- ,月,,,,, ,火,,,Design Engineering(4Q),, ,水,,システムデザイン特論,構造研ゼミ(後・青・石),地域産業(3Q), ,木,,,,, ,金,,構造力学特論,水理学特論,地震防災特論, ＜その他の予定＞~ 11/28(木)ポスター発表~ 12/27(金)実家に帰省~ 1/15(水)2025年初ゼミ~ 1/22 進捗ゼミ最後~ 1/24 就活1~ 2/3 就活2（前後数日は秋田に不在）~ 学会(出る予定)~ 全国学会(9月),YEC(5月)~ 概要はどちらも3月頃まで？（春休みの帰省や就職活動もそれを踏まえて計画を立てないとな...）~ *Model300 [#n498dba5] **解析結果 [#cff05bf9] 解析手法を確立させたら,まずはC1を破断させて,そのときの主桁鉛直変位を見る. **モデル作成 [#je113eea] Salome-Mecaで橋モデル（model300）を作成してみる.~ 以下の写真らはMarcで卒業された先輩が作成した橋モデルである.これを参考にして,まずはモデルの完成を目指す.（Marcではm系でモデルを作成していたため,Salomeも同様に作成する.）~ 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かな。）を用いて解析を行う.~ また、プレストレスは一旦温度応力で試してみようと考えている. ただ,ケーブルは降伏点が2つ存在しているため,TRACTIONを用いて降伏点を設定する必要があるが,現時点で解析が成功した試しがない...~ **解析記録 [#gb613b11] 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のデータを参照）~ 活荷重：主桁中央に1.175×10^6N/m, 桁全体に35250N/m (Marcのデータを参照）~ 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] 6/12:とりあえず,モデル作成優先でいいらしいので,まずはモデルを完成させよう。 影響線：移動荷重に対してある点に生じる物理量の大きさを移動荷重の位置関数で表したもの~ 影響線の定義から,まずは移動荷重を求め方が分からないと話にならない？ ↓橋の上を走る電車を例に移動荷重を算出する方法を見つけた。荷重条件に式を与える方法を使って,動的解析をする必要があるらしい？~ 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]] ちなみに...midasだと[[説明書の10章:https://in.midasit.co.jp/wp-content/uploads/2020/07/constrcution-83-civil-book.pdf]]を見たところ,影響線を導出する手法が記載されていた. *簡易モデルを用いた検証 [#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] 現状は、ケーブルに健全時の構成則と腐食時の構成則を同時に与えることはできている。ここから、腐食時のケーブルだけを切る（応力0に落とす）ことができない。~ 材料非線形はケーブルのみにしている。（桁に入れても,同じ材料に異なる構成則を同時に入れることは可）~ 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 現状~ 以下のようなエラーメッセージが出る。ケーブルの一部（要素分割した頂点部分）の接線剛性(応力-ひずみ曲線の傾き)が0とのこと。 http://www.str.ce.akita-u.ac.jp/~gotouhan/chiyooka/corrosionerror_1.png 材料特性あたりの設定があやしい？ **簡易モデル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] 研究室での中間発表や学会など、質問されたことをここにまとめ、再度確認したものをここに載せる。~ 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]] *やることリスト [#ra5b8814] 修士2年に上がるまでにやること~ -部分腐食モデル再現をある程度納得できるものに。 -橋梁モデルを解析できるように下準備（ケーブルは分割して） 少なくとも、この2点をある程度できてないと、学会行きたくても何も発表できないよな...。~ 頑張ろう！ 2024 夏休み期間にやること~ ・ケーブルを梁要素にした斜張橋モデルの解析をできるようにしたい。 ・卒論で扱ったケーブル2本の簡易モデルの詳細調査（timestepを0.001に変える）→後輩が行う簡易モデルによる実験での検討に繋げるために情報を集めたい。 ・帰省前に鋼構造を語る会で話す内容をまとめておく。（今の研究状況的に昨年度の卒論内容を改良したものを話すことになるかな？）~ 2024 6月~ ・橋梁モデルのmatrix errorを解決する。~ ・DAF関連の文献を読む.（英語だけど...）~ ・Salomeの勉強~ ・英語の勉強（英語で研究のことを含め,自分のことを話せるようになることを目標。）~ →今は中学レベルの文法を勉強し直している。（図書館で借りた英語のハノンという本） 2024年も残り半分か...。 2024 5月~ ・model300のasterstudyの設定方法を模索~ (主桁中央鉛直変位を求める。)~ 後は、研究の合間にDAF絡みの資料や論文読んでみる。勿論、授業関連の勉強も忘れずに。~ それと,wikiのビーム要素のページも適宜更新しておかないとな...。 2024 4月~ ・Salome-Mecaで橋モデル作成（1か月くらいで解析できるところまでは持っていきたい.） →まずは静的解析で主桁中央鉛直変位と比較してみる。 ・色々知識をつける。この1年は橋梁絡みを中心に土木のことを勉強していきたい。 ・DAFの式を論文等で調べる.別の式で検討する必要があるかもしれないので...。 後は、いつものことだが、コミュ力ももっとつけていきたいな。人に自分の言っていることが伝わるようにするにはどう話せばいいかできる限り意識しながら会話していこう。 *橋梁写真 [#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 撮影場所：秋田県由利本庄市~ 秋田県にある斜張橋です。~ 後輩たちがこの橋梁を研究対象に活動を頑張っているのを見て、自分も見に行ってきました！ 朝に行って来たけど、交通量もそれなりにあって、散歩で使っている人もちらほら見かけて、周辺に住む人達の生活にとって大事な橋だなと実感しました。~ ちなみに...昔はトラス橋だったらしく、架替え？にあたって今の斜張橋になったらしいです。 *日記 [#m48ef041] 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;~
