直交異方性材の設定における名前解決の不具合

直交異方性材の設定について https://gitlab.com/FrontISTR-Commons/FrontISTR/-/issues/373

こちらの問い合わせですが、以下のような状況が判明しました。

cnt ファイルの方で、以下のように局所座標系およびセクションを定義しています。

!ORIENTATION, DEFINITION=COORDINATES, NAME=aaa
1.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 0.0
!SECTION, SECNUM=1, ORIENTATION=aaa

このとき、両者の紐付けに名前(例えば aaa)をつけて指定しますが、内部的に片方(g_LocalCoordSys(cache)%sys_name)は大文字に変換されます。

一方、これと比較される変数 sect_orien は読み込んでそのままの値なので、aaa と AAA を比較して flase と判定されます。

この実装は意図しない動作なので、不具合と思われます。

この場合、

  • (1) セクションに紐付いた局所座標系を検索するときに、FrontISTR 内部では対応する局所座標系がないと判定
  • (2) 局所座標系のデフォルト値として全ての値が 0 のまま計算続行
  • (3) 応力ひずみ関係行列の値が 0 として計算続行
  • (4) 対角成分の値が 0 と判定され計算終了 という動作がされていたようです。
Edited by morita