#author("2022-10-07T13:45:14+09:00","default:kouzouken","kouzouken") #contents *修論テーマ [#p5d9b303] ケーブル腐食を考慮した単弦タイドアーチ橋の終局強度および連鎖崩壊 **進捗(Marc) [#n63d50a5] ***今後やること(2022/10/7現在) [#v1c48385] ・動的解析を教えてもらう~ ・3Dモデルで解析ができるようにする~ ・実験を行う(根拠がほしいから?)~ ***2022/7/22(解析できたけど値が合わない) [#j01fae8a] 中村先生からいただいたモデルmodel_140Aを解析~ 解析はできたものの、スパン中央の変位が一致しない(TT)~ 中村先生の変位→-263.5mm、私の変位→-380mm、、、なんでかな〜~ ↓model_140A↓ https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/kimijima/2022_7_22_model_140A.png~ **進捗(Salome) [#n1f933bc] ***2022/6/2 (弾塑性解析できた!) [#l3f18499] ソリッド要素→弾塑性、梁要素→線形での解析はできた!~ エラーになっていたのはPOUTREを使っているからではなかった~ 弾塑性解析を行うときは、POUTRE>VARI_SECT>HOMOTHETIQUEにしておく必要があったらしい?~ とりあえずケーブルは線形で解析を行うことにしたので、~ 梁要素にPOU_D_Eを使うかPOU_D_EMを使うか問題はいっかい無視し、POU_D_Eを使った~ (POU_D_EMを使うとエラーが出てしまうのは解決していない)~ https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/kimijima/2022_6_2_salome_asuta.png~ 1000Nをかけた結果、ケーブルと繋がっている(下図でいう右端)ところの変位は約0.8mmになった~ 角田さんの簡易モデル(桁もケーブルも梁要素)では変位が3.356mmだったのだが、大丈夫?(条件は同じ)~ https://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/kimijima/2022_5_19_kanimodel.png~ 主応力ーひずみグラフ(固定端の上部を見ている)↓~ https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/kimijima/2022_6_3_gurafu.png~ ***2022/5/27 (弾塑性解析できない) [#reda29e7] 弾塑性解析をしようとするとエラーになる~ Error of use: ! ! One does not manage to correctly build the field containing the number of ! ! subpoints ! ! of the finite elements (multi-layer shells, multifibre pipes, beams,') model ! ! model. ! ! ! ! Risks & advices: ! ! This error intervenes when one does not define TOUTES the elementary ! ! characteristics ! ! in the same AFFE_CARA_ELEM. ! ! For the commands of computation, it is necessary one MODELE and that only one ! ! CARA_ELEM. https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/kimijima/2022_5_27_salome_asuta.png~ AFFE_CARA_ELEMではPOUTREを使っているのだが、線形解析ではつかえないのかな?と思い~ 他のBARREやCABLEなどを使ってみたところ、~ Error during the interpolation of the function func. ↑のエラーが出た。~ でも、POUTREを使うと最初のエラーが出てしまう無限ループにはまった。~ -AFFE_CARA_ELEM...要素特性を定義するコマンド~ POUTRE...はり要素の特性/配管要素の特性~ BARRE...棒要素およびケーブル要素の特性~ あと、弾塑性ではPOU_D_Eが使えないらしいからPOU_D_EMを使っている~ (自分で調べても出てこなかったのだが、角田さんの日誌に書いてあった)~ -Eulerはり(POU_D_E*):極端に細長い~ 多芯繊維(POU_*_M):材料非線形という意味らしい~ ***2022/5/20 (線形解析できた!!) [#r3baa6a0] ↓解析できた簡易モデル(100×100×1000mm)↓~ https://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/kimijima/2022_5_20_kanimodel.png~ ソリッド要素と梁要素を組み合わせて解析するときは、結合部を一番端にしてはいけないらしい~ このモデルは一番端から100mm離した点を結合部としている~ ***2022/5/19 (エラーうまくいかない)[#zcb76729] これまで~ 南方澳大橋のモデル作成できた~ →全部ソリッド要素で解析しようとしたが、橋はソリッド要素でケーブルは梁要素で解析することにした~ →まずは簡易モデルでやってみる~ →ソリッドと梁を組み合わせた解析がうまくいかない…(TT) (多分境界条件がおかしい?) ↓簡易モデル(100×100×1000mm)↓~ https://www.str.ce.akita-u.ac.jp/~gotouhan/j2021/kimijima/2022_5_19_kanimodel.png~ -エラー ! impossibility, the node N481 carries the degree of freedom of rotation DRX ! 訳→不可能な場合、ノードN481は回転の自由度DRXを保持します~ DRXがおかしいと言われる、角田さんと同じエラー、大体はこっちのエラーだった~ Aucune des mailles sélectionnées n'est affectée par un élément fini. 訳→選択したメッシュはいずれも有限要素の影響を受けません。~ メッシュはちゃんと作ってますけど?意味不明~ --境界条件DX,DY,DZ,DRX,DRY,DRZの値を変えてもダメ~ --載荷する部分を変えても同じエラー(自由端の上部下部への載荷、上部の橋から10mm離したところへの載荷、面載荷を試した)~ --モデルを作成しなおしてみたが同じエラー~ --メッシュを切り直してみたり、違う方法を試してもだめ~ --AsterStudyを色々いじってもダメだった、何がおかしいのかわからない *Marc解析やり方[#m0975904] 私がわかりやすいようにメモしておきます ~ **片持ち梁モデル(梁) [#e0a9687c] ***モデル作成 [#e2964c5d] -長さ単位 (基本操作)~ 長さを選ぶ(私はmm)~ -形状とメッシュ (基本操作)~ ポイント 追加 → 座標入力する(0 0 0 とか)(x y z の順,スペース開ける) → 節点 追加 (ポイントの上から節点を作る)→~ 要素 追加 → 作った節点と節点を左クリック → 右クリック確定 → 線ができる~ -再分割 (オペレーション)~ 分割 → 上からx,y,zの順に値を入力 → 要素 → 分割したいとこ左クリック → 右クリック確定 → 分割できる~ -スイープ (オペレーション)~ 節点 → 全て → OK (必要のない節点を消す?みたいな)~ -再番号付け (基本操作)~ 全ての形状/メッシュ → OK (節点の番号をつけ直す?みたいな)~ ***形状設定 [#g2aa48c6] -新規(構造)▼ (形状特性)~ 3ーD ▶ → 中実断面梁(簡単な断面のとき使える) → 形状選ぶ(長方形,正方形とか) → 寸法入力 →~ X,Y,Zの値は変えたり変えなかったり(今のところよく分かっていない) → 要素 追加 →~ 追加したいところ左クリック(全体選択でも可) → 右クリック確定 → 左側の要素の数増えてたらできてる~ ***材料設定 [#bfe74861] -新規(構造)▼ (材料特性)~ 有限剛性領域 → スタンダード → ヤング率とポアソン比入力(n/mm^2で入力) → 要素 追加 →~ 追加したいところ左クリック(全体選択でも可) → 右クリック確定 → OK~ ***境界条件設定 [#j0b785bb] -(固定)~ 新規(構造)▼ (境界条件)~ 変位指定 → 変位X,Y,Zと回転X,Y,Zにチェックを入れて0にする(固定端のとき) → 節点 追加 →~ 固定端にする方の節点を左クリック → 右クリック → 矢印出てきたらできてる~ -(ローラー)~ 新規(構造)▼ (境界条件)~ 変位指定 → 変位X,Yにチェックを入れて0にする(ローラー端のとき) → 節点 追加 →~ ローラー端にする方の節点を左クリック → 右クリック → 矢印出てきたらできてる~ -(載荷)~ 新規(構造)▼ (境界条件)~ 集中荷重 → 力Yに値を入力 → 節点 追加 → 載荷するところの節点を左クリック → 右クリック → 矢印出てきたらできてる~ ***解析 [#z06f102b] -新規▼ (解析ジョブ)~ 構造 → チェック → Job check complete になったらOK → 実行 → Marc実行(1) → 3004になったら解析できてる~ ***結果見る [#w2500455] -ポストファイルを開く(モデルプロット結果処理メニュー) → 変形形状の形式を 変形図と元の形状図▼にする → 設定 → 自動にする → 変形図見れる **簡易モデル(梁) [#s592f8fa] *Marcメモ [#q67f2a8f] -左クリック…選択,右クリック…確定,ホイールクリック…選択取消~ -ポイントの座標はいちいち表示しないといけない(なぜ?めんどくさい…)~ --ポイント表示のやり方……ビュー(V) → モデルエンティティタイプ → ポイント → 表示される