論文要約：粗視化されたトポグラフィを解像する構造動的接触シミュレーションの実現

論文タイトル: Enabling topography-resolving structural dynamic contact simulation
著者: Hendrik D. Linder, David A. Najera-Flores, Robert J. Kuether, Malte Krack
所属: シュトゥットガルト大学、ATA Engineering, Inc.、Sandia 国立研究所

1. 研究の背景と課題

構造物やシステムの減衰は、接合部（ボルト接合やリベット接合など）における摩擦散逸によって支配されることが多く、その予測は長年の科学的課題となっています。特に、部分的なすべり（partial slip）や局所的な離脱（liftoff）が生じる条件下では、実際の接触面の**トポグラフィ（表面粗さや形状誤差）**が動的挙動に決定的な影響を与えます。

従来のアプローチには以下の課題がありました：

• 計算コスト: 製造された表面の微細なトポグラフィを有限要素法（FE）で直接解像すると、要素数が膨大になり、長時間の過渡動的シミュレーションは計算的に不可能、あるいは極めて非効率です。
• 数値的安定性: 接触界面を細かく解像し、かつ厳密な接触法則（Coulomb-Signorini 条件）を適用する場合、FE 解析では数値誤差や非収束、人工的な高周波振動が発生しやすく、物理的な摩擦減衰を正確に評価することが困難です。
• 既存手法の限界: 準静的解析（QSMA）は広く用いられていますが、非線形性が強い場合や、変位形状が変化する動的プロセス、振動による静的平衡位置のシフト（セトリング）を捉えることができません。

2. 提案手法：動的マルチスケール法（MSM）の拡張

本研究では、以前提案された「FE-BE マルチスケール法」を動的解析に拡張しました。この手法は、接触領域と構造全体を異なるスケールでモデル化することで、トポグラフィを解像した動的シミュレーションを可能にします。

2.1 手法の概要

• 接触領域（BE モデル）: 接触界面の微細なトポグラフィと接触力学は、半空間理論（Half-space theory）に基づいた境界要素法（BE）でモデル化します。規則的なグリッド上で実装されており、剛性行列が閉形式で表現可能です。
• 構造体（FE モデル）: 構造全体の動的挙動は、比較的粗い有限要素（FE）モデルで記述します。
• 結合: 両モデルは、遠方領域での変位適合性と力平衡条件を課すことで結合されます。

2.2 動的解析への拡張と数値的工夫

動的解析への適用において、以下の重要な改良がなされました：

1. 質量なし境界のコンポーネントモード合成（Massless-boundary CMS）:
   • 接触界面に質量を持たないよう、Craig-Bampton 法を修正して適用しました。これにより、接触閉塞時の衝撃による人工的な高周波振動を防止し、数値減衰（アルゴリズム的減衰）を不要にしています。
2. 半陽的 Leapfrog 時間積分スキーム:
   • 内部の線形動的挙動（FE モデル）を陽的に、接触問題（BE モデル）を陰的に（増分ラグランジュ法を用いて）解くハイブリッド手法を採用しました。
   • この構成により、数値的安定性を保ちつつ、メッシュサイズに依存しない比較的大きな時間ステップを使用可能になりました。
3. 調和バランス法（Harmonic Balance）の実装:
   • 定常周期応答を直接計算し、時間積分手法の数値減衰の有無を検証するための基準として用いました。

3. 検証と結果

S4 Beam ベンチマーク（実際のトポグラフィ測定データを持つ 2 枚の C 型ビームのボルト接合）を用いて、提案手法の有効性を検証しました。

3.1 計算効率とロバスト性

• 計算速度: 提案するマルチスケール法は、陽的 FE 解析と比較して約 80 倍、陰的 FE 解析と比較して約 400 倍の高速化を達成しました。
• 時間ステップ: 提案手法は、FE メッシュサイズに制約されず、保持されたモードの最高周波数によってのみ制限されるため、非常に大きな時間ステップ（約 24 μs）を使用でき、長期間のシミュレーションが可能でした。
• 数値減衰の欠如: 調和バランス法との比較により、提案する時間積分スキームには数値減衰がほとんど含まれていないことが確認されました。

3.2 結果の一致と物理的現象

• 準静的解析との整合性: 振幅依存する固有振動数と減衰比について、提案手法と完全 FE 解析の結果は非常に良く一致しました。
• 振動誘起セトリング（Vibration-induced Settling）の発見:
  • 衝撃荷重後の自由減衰（リングダウン）解析において、システムが予荷重ステップ後の平衡状態とは異なる静的平衡状態に落ち着く現象を初めて数値的に明らかにしました。
  • これは、摩擦接触の履歴依存性により、残留接触応力場（Residual traction）が形成され、平衡位置がシフトするためです。
  • QSMA（準静的解析）はこの現象を捉えられませんが、動的シミュレーションでは明確に観測されました。
• 接触圧力分布: 提案手法と完全 FE 解析の間で接触圧力分布に定量的な差異（端部効果やコンプライアンスの重複によるもの）が見られましたが、低次モードの動的特性にはほとんど影響を与えないことが確認されました。

4. 主要な貢献と意義

1. トポグラフィ解像動的シミュレーションの実現: 微細な表面粗さを考慮した構造動的接触シミュレーションを、計算的に実行可能かつ効率的に実現しました。
2. 高精度な数値手法の提案: 質量なし境界 CMS と半陽的 Leapfrog スキームを組み合わせることで、数値減衰や人工振動を排除し、物理的な摩擦減衰のみを正確に評価できる手法を確立しました。
3. 新たな物理現象の解明: 振動による接合構造の「セトリング（平衡位置のシフト）」とそれに伴う残留応力場の形成を初めて数値的に示し、実験結果の再現性問題や QSMA と動的解析の差異の理由を説明する新たな知見を提供しました。
4. 実用性: 複雑な接合構造の設計・評価において、トポグラフィの影響を考慮した信頼性の高い予測を可能にし、将来の振動制御や耐久性評価への応用が期待されます。

5. 結論

本研究は、接合構造の摩擦減衰予測における長年の課題に対し、トポグラフィを解像した動的マルチスケール法を確立し、その有効性を示しました。特に、数値的ノイズを排除した高精度なシミュレーションにより、振動誘起セトリングのような複雑な非線形現象を初めて捉えることに成功しました。この手法は、接合構造の動的挙動理解を深め、より信頼性の高い設計手法の確立に寄与するものと考えられます。