梅田の卒論日誌 の履歴(No.34)

• 履歴一覧
• 差分 を表示
• 現在との差分 を表示
• ソース を表示
• 梅田の卒論日誌 へ行く。
  • 1 (2020-01-19 (日) 15:15:06)
  • 2 (2020-01-24 (金) 16:09:06)
  • 3 (2020-01-29 (水) 22:12:20)
  • 4 (2020-01-31 (金) 19:49:29)
  • 5 (2020-02-04 (火) 14:50:58)
  • 6 (2020-02-06 (木) 13:11:42)
  • 7 (2020-03-07 (土) 18:37:24)
  • 8 (2020-03-23 (月) 19:11:42)
  • 9 (2020-03-24 (火) 17:20:07)
  • 10 (2020-05-11 (月) 11:27:33)
  • 11 (2020-05-17 (日) 08:01:13)
  • 12 (2020-05-17 (日) 13:00:08)
  • 13 (2020-05-22 (金) 10:55:51)
  • 14 (2020-05-22 (金) 12:32:32)
  • 15 (2020-05-28 (木) 11:49:51)
  • 16 (2020-05-29 (金) 09:27:23)
  • 17 (2020-05-29 (金) 17:20:43)
  • 18 (2020-06-05 (金) 10:36:45)
  • 19 (2020-06-06 (土) 01:56:58)
  • 20 (2020-06-06 (土) 01:56:58)
  • 21 (2020-06-12 (金) 07:04:57)
  • 22 (2020-06-12 (金) 20:10:14)
  • 23 (2020-06-19 (金) 08:21:12)
  • 24 (2020-06-22 (月) 11:52:56)
  • 25 (2020-08-03 (月) 16:07:13)
  • 26 (2020-08-24 (月) 15:15:09)
  • 27 (2020-10-16 (金) 18:00:48)
  • 28 (2020-11-25 (水) 15:53:52)
  • 29 (2020-12-01 (火) 14:55:40)
  • 30 (2020-12-08 (火) 10:37:58)
  • 31 (2020-12-21 (月) 19:24:33)
  • 32 (2020-12-24 (木) 17:46:06)
  • 33 (2021-01-07 (木) 17:20:47)
  • 34 (2021-01-28 (木) 18:50:21)

---

6月19日†

アイさんのファイルをもとに非線形のエラスティカ問題を解いています。Runはできたのですが、ParaVisでみるとどこの変位も０であったため、原因を探っています。

6月12日†

エラーの原因がわかり、計算した結果、アイさんの計算結果と同じような結果になりました。しかし、理論値とは離れた値になりました。アイさんのファイルではFYの値を2としていたので、同じようにしましたが、60Nのどこから２がでてきたのかいまいち理解できていません。原因はここではないのかと思っているので、数値を変えてみようと思います。

近藤より(20/06/12)†

アイさんの卒論概要の図-2を見て、混乱しているのだと思いますが、実はあのグラフは縦軸と横軸が間違っています。あのグラフは、ここの図-6.51の一番右のグラフを元にして描いているものなので、Salome-Mecaから得られた値（荷重P, 鉛直変位\( \delta \)）を次式に代入すれば再現できると思います。

• 縦軸：荷重 → \( \frac{P\ell^2}{EI} \)
• 横軸：変位 → \( \frac{\delta}{\ell} \)

6月5日†

アイさんのファイルをもとに線形のエラスティカ問題を解いていますが、またエラーが出たため、原因を探っています。

アイさんのファイルをSalome-Meca2018で計算してみた by近藤†

• hari_senkei.hdf を計算するときの注意点
  • AsterStudy → Output → Set output result → Result file location
  • /home/kouzou/2019/salome/ai/3nensei1.med を適当な場所と名前に変える。これだけでOK.
    • 多分、計算させるマシンに/home/kouzou/2019/salome/ai/という場所がないので、どこに結果ファイルを保存すればいいのかわからなくてエラーになっているのだと思う。（基本的に他のマシンで作ったhdfは、必ずここを修正する必要があります。）
    • ちなみに、このようなエラーのときは、（赤くなっていても）実は結果ファイル自体は/hari_senkei_Files/RunCase_1/3nensei1.medに保存されているので、それをParaVisに読み込めば、ちゃんと結果を確認できます。
• elastica.hdf を計算するときの注意点
  • Output の設定は↑(hari_senkei.hdf)と同様に修正する。それプラス↓を変える。
  • AsterStudy → History View → Basic
  • Memory : 2048 → 6000（マシンの性能による）
    • 要素数が多いときは、Memoryが足りていなくて計算ができなくなることがあるので、適宜増やす。
  • Time : "0:15:0" → "1:0:0"（計算時間による）
    • 非線形解析など、計算時間が15分以上かかりそうな場合はTimeを増やす。これは最大値を入力しても問題ない。
    • ちなみに近藤PCでは30分近くかかりました。
  • これでOK.
  • ParaVisで結果を見るときは、Timeを変えてください（詳しくはこちら）。
  • あと、elastica.hdfは、ちょっと結果がおかしい気がするので、commファイルの設定を見直した方がいいと思います。

5月29日†

  !---------------------------------------------------------------!
  ! <EXCEPTION> <FACTOR_57>                                       !
  !                                                               !
  ! Solver MUMPS:                                                 !
  !  The solution of the linear system is too vague:              !
  !  Computed error:  9.44269e-05                                 !
  !  Acceptable error:  1e-06 (RESI_RELA)                         !
  !                                                               !
  !  Advices:                                                     !
  !  One can increase the value of the key word SOLVER/RESI_RELA. !
  !---------------------------------------------------------------!
  上記のようなエラーがでたのですが、原因がよくわかりません。

　　何が原因なのでしょうか。

後藤ちゃちゃ(20/5/29)†

まずは、そのエラーメッセージ自体をコピペして検索してみて下さい。 そうすると、

• https://www.code-aster.org/forum2/viewtopic.php?id=24819 こんなページが引っかかったりします。 収束判定基準のRESI_RELAを1.e-6から1.e-5にひと桁あまくしてやると収束するとか、 そういう話みたいです。

後藤ちゃちゃ(20/5/28)†

Salome-Meca例題ファイル のところに，アイさんのエラスティカ問題（片持ち梁を幾何学非線形で解くものと，普通に線形で解くもの）のhdfファイルらしきものへのリンクを張っておきました．

ノートPCの起動OSをデフォルトでUbuntuが立ち上がるようにしたい場合は、 HPノートのデフォルト起動のOSを決めてしまいたい参照。

5月22日†

エラスティカ問題というものがどういうものなのかいまいちピンと来ず、とりあえず、10✕10✕1000、ヤング率7500GPa、ポアソン比0.4、荷重60Nの片持ち梁の変位を求めようとした。しかし、Runの砂時計がずっとぐるぐる回っており、Runが終わらない状態が続いている。

後藤ちゃちゃ(200522)†

梅田さんのSalome-Mecaは2017でしたっけ。 ちなみに、Salome-Meca2019の場合、実行中のHistory Viewの右下に出ている Outputタブ（ここに Cannot read output file or empty file. が表示されている）を Messagesタブに切り替えると、計算が終了していたり（ここで終了してなくても、 Outputに戻ると）終了したりする。 というわけで、実行後は、OutputとMessagesのタブを時々 切り替えるといいというのはありますが、Salome-Meca2017では、この症状はないと思います。

過去に計算できていたhdfファイルは計算できるのでしょうか。

5月17日†

線形の初等梁を解いているのですが、AsterStudyでRunをクリックすると This operation needs the Study directory. So, you should save the study before executing this operation.という表示が出てきてRunが行えません。 なぜRunができないのでしょうか。

20/5/17後藤ちゃちゃ†

言われた通りにstudyファイルを保存したらどうでしょう。

当面の課題†

20/5/11(後藤)†

• アイさんのハードディスクから、片持ち梁のエラスティカ問題を解いているhdfを探して、解けるかどうか確認。
  • 線形で普通に解くと、初等梁のたわみ(\( \frac{P\ell^{3}}{3EI} \))と一致するか
  • 非線形で解くと楕円積分解（アイさんの卒論にグラフあり）と一致するか
    • Salome-Mecaで幾何学非線形で解くにはどうするか。
• 以上が解決したらSalome-Meca例題ファイルの幾何学非線形のところに、hdfファイル（計算前のもの）をアップ

日誌†

10/18課題†

10月25日課題†

10月25日課題の考察†

メッシュサイズが１と１．８の時の要素数、先端変位が前のメッシュサイズの要素数、先端変位よりも少し大きいことが気になった。原因がわからないので、もっと勉強したい。

11月1日課題†

理論値　0.416667

11月8日課題†

理論値　0.4917

11月15日課題†

・分かったこと

　・作業の全体の流れ

・分からなかったこと

　・エンティティの追加の際、ソリッドなどの複数ある項目の中からどれを追加すればよいのか

　・「解析ジョブ」の＜要素タイプ＞を選択する時の細かい選択手順

　・境界条件のエラーでの「自由度は規定されていません！」とは一体どういうことなのか。

　・固定面を固定する際の細かい作業手順

11月22日課題†

理論値　6.6667

12月6日課題†

理論値　初等梁6.66933 ティモシェンコ梁6.742