[論文レビュー] Computational mechanics of soft filaments

Mattia Gazzola, Levi H. Dudte+2 more|arXiv (Cornell University)|Jul 1, 2016

Advanced Materials and Mechanics被引用数 5

ひとこと要約

本論文は、SE(3) 上での完全なコッサレートロッド理論を用いて、曲げ、ねじり、引き伸ばし、ずれの6自由度変形をすべて経験するソフトフィラメントの動的挙動をシミュレートする、ロバストで完全に整合性のある計算フレームワークを提示する。この手法は、摩擦、流体力学的抵抗、筋肉の能動的力などの複雑な環境相互作用を正確にモデル化でき、肢のない運動やプリエクトネ形成といった前向きおよび逆問題において優れた性能を示す。

ABSTRACT

Soft slender structures are ubiquitous in natural and artificial systems and can be observed at scales that range from the nanometric to the kilometric, from polymers to space tethers. We present a practical numerical approach to simulate the dynamics of filaments that, at every cross-section, can undergo all six possible modes of deformation, allowing the filament to bend, twist, stretch and shear, while interacting with complex environments via muscular activity, surface contact, friction and hydrodynamics. We examine the accuracy of our method by means of several benchmark problems with known analytic solutions. We then demonstrate the capabilities and robustness of our approach to solve forward problems in physics and mechanics related to solenoid and plectoneme formation in twisted, stretched filaments, and inverse problems related to active biophysics of limbless locomotion on solid surfaces and in bulk liquids. All together, our approach provides a robust computational framework to characterize the mechanical response and design of soft active slender structures.

研究の動機と目的

• 従来の不変形性ロッドモデルの制限を超えて、曲げ、ねじり、引き伸ばし、ずれの6モードの変形をすべて経験できる、数値的にロバストで効率的なソフトフィラメントのシミュレーション手法を開発すること。
• 異方性摩擦、流体力学的抵抗、表面接触、能動的筋肉力といった現実的な物理的相互作用を、統一されたフレームワークに統合すること。
• 解析的ベンチマークとの比較による妥当性評価を行い、ソフトアクティブメカニクスにおける複雑な前向きおよび逆問題の解決能力を示すこと。
• 非対称な摩擦と筋肉活性化パターンによって駆動される陸上および水中移動をシミュレートすることで、肢のない生物における逆設計による運動戦略の開発を可能にすること。

提案手法

• 中心線の位置と各断面における局所的向きの独立した時間発展を可能にする、完全なユークリッド群 SE(3) に基づく離散的コッサレートロッド定式化を採用する。
• エネルギーおよび運動量保存を保証するための高速な離散的作用素を備えたハミルトニアン構造を用いることで、数値的安定性と精度を確保する。
• ヤング率や剪断率といった材料特性を用いて、曲げ、ねじり、引き伸ばし、ずれの各々に対して別々の剛性行列を有する一般化された物性則を導入する。
• 等方的および異方性摩擦モデル（静止摩擦および速度依存閾値を含む）を用いたペナルティベースの接触力処理を実装し、安定した接触シミュレーションを実現する。
• トルクプロファイルのBスプラインパラメータ化を用いて、空間的・時間的に変化する内部モーメントおよび力として筋肉の活性化をモデル化する。
• 流体抵抗および地面接触力は、粘性抵抗および剛性項としてモデル化され、地面摩擦は速度および方向に依存する区分線形力則に従う。

実験結果

リサーチクエスチョン

• RQ1 完全に柔軟なロッドモデルは、ねじりや引き伸ばしを受ける複雑な変形モードを示すフィラメントに対して、既知の解析的解をどれほど正確に再現できるか？
• RQ2 制御された荷重条件下で、ねじり・引き伸ばしを受けるフィラメントに、プリエクトネやソレノイドの形成がロバストに再現可能か？
• RQ3 対称的な筋肉活性化パターンと組み合わせた異方性摩擦は、肢のないフィラメントでどのようにネットワークの運動を可能にするか？
• RQ4 物理的に現実的なシミュレーション環境において、効果的な陸上および水中の肢のない運動を実現するために必要な最小限の機械的および制御パラメータは何か？
• RQ5 このフレームワークは、異なる環境で前進速度を最大化する最適な筋肉活性化パターンを特定するような逆問題を解けるか？

主な発見

• 数値的手法は、大変形および複雑な境界条件を伴う静的および動的ロッド配置に対しても、解析的解を正確に再現する。
• ねじり・引き伸ばしを受けるフィラメントにおいて、プリエクトネおよびソレノイドの形成が成功裏に再現され、理論的期待と一致する変形パターンが得られた。
• 異方性摩擦により、対称的な筋肉活性化を持つフィラメントでも、前進方向のネット運動が実現可能であることが示され、摩擦による対称性の破壊が有効な運動を可能にする。
• 筋肉活性化に「進行波成分」が存在しない状況でも、摩擦の異方性によりネット変位が達成されることから、モデルの物理的妥当性が裏付けられた。
• 等方的および異方性摩擦の両条件において、対称的なケースで不自然な運動や数値ノイズが発生せず、スキームのロバスト性が確認された。
• 微小時間刻みで安定かつ高精度な長時間積分が達成され、多様な機械的および環境的条件下でもエネルギー保存則および物理的忠実性が維持された。