有限回転シェルのドリリング項の接線剛性を内部力と整合させる

概要

有限回転の 741 シェルで、ドリリング安定化項の接線剛性から、内部力に無い項を除いた。ShellMITC_AddDrillingStiffness が Cv_deriv_disp にダイレクタの回転による Cv_theta の微分を含めていたのが原因で、全回転に比例する非対称な項が剛性に入っていた。弾性材料では直接法が対称モードで上三角だけを使うので、Newton 法が1次収束に落ちていた。

修正後、ドリリング項の接線は内部力の中心差分と相対 1e-7 で一致し、MUMPS での反復回数は次のとおり減った。収束後の変位は相対 3e-9 以下しか変わらない。

ケース（片持ち板、先端たわみ約 2.0）	修正前	修正後
等方、alpha_over_mu 1e-3（既定）	12	8
等方、alpha_over_mu 1e-2	26	8
flag=1（繊維 90°）	32	8

内部力と、有限回転以外の剛性は変わらないので、参照解は更新していない。有限回転シェルの既存テスト（static_lib_shell_mitc4、CG）は反復回数も変わらない。

Close #904 (closed)

変更内容

コミット	要約
88c3ff94 Fix the drilling tangent of finite-rotation MITC4 shells	Cv_deriv_disp への加算を、ダイレクタの回転による項を足す前に移す