Give the 2box_tet_fric contact tests an actual friction coefficient

Kazuya Gotorequested to merge
kgoto/FrontISTR:feature/tet-fric-friction into master

概要

2box_tet_fric 系テスト 5 ケース（contact_slag_mumps / contact_slag_mkl / contact_slag_iter / contact_alag_iter / contact_slag_iter_ml）の摩擦係数を 0.0 から 0.3（2box_hex_fric 系と同じ値）に修正し、参照結果を再生成しました。これらは摩擦係数 0 のため対応する 2box_tet_nofric と実質同一の重複テストとなっており、四面体要素×摩擦接触が実際にはテストされていませんでした。

Close #902 (closed)

変更内容

  • 5 ケースの !CONTACT データ行の摩擦係数を 0.3 に変更（各ファイルの既存書式を維持。それ以外の hex/tet 間の既存差異は変更していません）
  • contact_alag_iter の当該ケースのみ !SOLVER の METHOD を CG から BiCGSTAB に変更（摩擦による非対称剛性で CG は 2 プロセス実行時に収束しないため。同ディレクトリの 2box_hex_fric_iter と同じ選択です）
  • 参照結果 5 件を tests/create_reference.sh（単一プロセス、-t 1）で再生成

確認内容

  • ctest: serial 全件 120/120 PASS、mpi の with_mumps 9/9 PASS、mpi の interaction 接触系 10/10 PASS
  • tests/test.sh 直接実行: 修正した 5 ディレクトリを 1 プロセス・2 プロセスで全ケース success（with_mkl は 1 プロセスのみ。手元環境では Cluster PARDISO が実行できないため、MPI 並列側は CI での確認をお願いします）
  • 検出能: 修正後の入力を旧参照（μ=0 の結果）と比較すると compare_res.pl が 5 ケースすべてで failure となることを確認（摩擦量 CONTACT_FRICTION / CONTACT_FTRACTION に明確な差）
  • 参照再生成時の収束: SLAGRANGE 系 4 ケースは Newton-Raphson 3 反復で res(force) 2E-11 まで収束、ALAGRANGE ケースも正常終了

補足

  • 現行 master の fistr1 は参照生成時にステップ 0 の結果（res.0.0）も出力しますが、各ディレクトリの既存参照が res.0.1 のみ追跡している慣行に合わせ、res.0.0 はコミットに含めていません

Merge request reports