Add regression test for torque load in explicit dynamic analysis

Toru Andorequested to merge
mitsumega_working/front-istr-ando:test/explicit-torque into master

upstream master (81b149cf) の現状を確認しました。

  1. 回転変位 (!BOUNDARY, ROT_CENTER=): 107840ce で陽解法に対応済み (tests/analysis/dynamic/explicit_rot)。
  2. トルク (!CLOAD, ROT_CENTER=): 荷重組立が陰陽共有のため陽解法でも動作します。回帰テストが無かったので MR !NNN で tests/analysis/dynamic/explicit_torque を追加しました (陰解法との差 0.36%、静的ねじれ角の 1.93 倍で動的増幅と整合)。
  3. 梁・シェルの回転自由度 (ndof=6) への回転変位・トルク: 未対応。境界条件側は DOF 4-6 を通しますが、陽解法は全自由度の集中質量が非零であることを起動時に要求します (fstr_dynamic_nlexplicit.f90 の stop due to fstrDYN%VEC(j) = 0)。梁 611 は回転自由度の質量を 0 にしており (fstr_EIG_setMASS.f90)、MITC シェルはディレクタ方向の成分が 0 になるため平板ではドリリング自由度で同じ停止に落ちます (コード読みによる予測、未実測)。対応には回転慣性 (梁の断面慣性、シェルの人工ドリリング慣性) の与え方と臨界時間刻みへの影響の設計が必要です。

Merge request reports