ALagrange接触でカットバック発生時に行列構造を更新する

close #845 (closed)

カットバック後の行列構造再構築処理の条件を下記に変更
if( fstrPARAM%contact_algo == kcaSLagrange ) から if( is_interaction_active ) に変更し、接触アルゴリズムによらず常に fstr_mat_con_contact() / solve_LINEQ_contact_init() を呼ぶようにすることで解消した。
変更ファイル:fstr_solve_NLFEOM.f90、fstr_dynamic_nlimplicit.f90

検証用例題として下記をtests/analysis/interaction/hertz_autincに追加

  • tutorial/09_contact_hertzの例題をベースに下記内容に変更
    • 材料をMooney-Rivlinへ変更
    • 固定時間増分を自動時間増分へ変更
    • STEP1で平板に押し込み、STEP2で元の位置に戻すという解析へ変更

Merge request reports

Loading
Loading