Finite strain homogenization of periodic rod networks with application to semi-flexible biopolymers

技術要約：周期的なロッドネットワークにおける有限歪みホモジナイゼーション

問題提起
コラーゲンやフィブリンのような半柔軟なバイオポリマーネットワークは、歪み硬化、引張下での体積収縮、圧縮下での軟化、およびせん断下での逆ポインティング効果を含む、複雑で高度に非線形な機械的挙動を示す。これらの現象は、材料の階層構造と、曲げや座屈といったマイクロスケールの繊維力学とマクロな変形との相互作用に起因する。現象論的な連続体モデル（アフィンおよび非アフィン鎖モデルなど）は存在するが、それらはしばしば、これらの挙動を駆動する具体的な微視的メカニズムに関する洞察に欠けている。一方で、ランダムな繊維ネットワークのフルスケール計算シミュレーションは、計算コストが非常に高い。本研究は、ホモジナイゼーションを通じて計算効率を維持しつつ、大規模で非アフィンな微視的変形（座屈を含む）を捉える計算フレームワークの必要性に対処するものである。

手法
著者らは、周期的な8本鎖および14本鎖の代表体積要素（RVE）に適用される、有限歪み計算ホモジナイゼーション手法を提案している。この手法は、以下のコンポーネントを統合している：

• マイクロスケール・モデリング： 個々の繊維は、幾何学的厳密（GE）なコセラット・ロッド理論を用いてモデル化されている。この定式化により、オイラー・ベルヌーイ法やティモシェンコ梁理論のような小変形の制限を受けることなく、伸長、曲げ、ねじれ、せん断を含む任意の大きな変形を捉えることが可能となる。
• ホモジナイゼーション・フレームワーク： 問題は、歪み駆動型と応力駆動型の両方を用いて定式化される。
  • 歪み駆動型： マクロな変形勾配（$F^M$）が規定され、制約条件下でRVEのエネルギーが最小化される。
  • 応力駆動型： マクロな第一ピオラ・キルヒホッフ応力（$P^M$）が規定され、マクロな変形勾配が未知数となる。これは、実験で一般的である、側面の無応力境界条件をモデル化するために極めて重要である。
• 制約および境界条件： RVEは、内部ジョイント制約（剛またはボールジョイント）、コーシー・ボルン則から導出された周期境界条件（PBC）、および剛体モードを除去するためのグローバルな制約を受ける。
• 数値解法： 制約付きエネルギー最小化問題は、経路追跡のためのアーク長継続法で拡張されたニュートン・ラフソン法を用いて解かれる。これにより、人工的な不完全性を導入することなく、座屈後の解を追跡することが可能となる。
• 安定性解析： 得られた解の安定性は、「ハード・ローディング」（構成変数のみの摂動）と「ソフト・ローディング」（構成変数とマクロな変形勾配の両方の摂動）という2つの基準を用いて評価される。後者はより厳格な条件である。

主な貢献

1. 特異なバイオポリマー現象の回収： 本研究は、ホモジナイズされた8本鎖および14本鎖のロッドネットワークが、生体組織に見られる特定の現象（歪み硬化および一軸引張下での劇的な体積収縮、圧縮下での突然の軟化、および逆ポインティング効果による引張法線応力の生成）を回収できることを示している。
2. メカニズムの特定： 著者らは、非線形性と非アフィン変形が、特に繊維の曲げと座屈が、これらのマクロな現象を駆動する主要な要因であることを特定した。
3. 不完全性を排除したポスト座屈解析： 本研究は、初期の幾何学的不完全性に依存することなく、非線形経路追跡法を用いてポスト座屈解を得ており、ハードおよびソフト・ローディングの基準を通じて安定性を厳密に検証している。
4. 検証と拡張：
   • モデルは、コラーゲンの一軸引張に関する実験データ（Roeder et al., 2002）およびフィブリンに関する解析モデル（Brown et al., 2009）に対して検証されている。
   • 圧縮応答は、血栓実験（Liang et al., 2017）と比較され、RVEサイズと繊維の細長比を考慮した場合に良好な一致を示す。
   • フレームワークは、8本鎖ユニットセル内にヘリカルロッドを含めるように拡張されており、コンプライアントなメタマテリアルで見られる拡張された伸長挙動を正常に捉えている。

結果

• 一軸引張：
  • 歪み駆動型： 応答は主に伸長支配的であり、線形である。
  • 応力駆動型： 応答は歪み硬化を示す。初期段階では変形は曲げ支配的であるが、繊維が引張方向に整列するにつれて、モードは伸長へと移行する。この遷移が非線形なポアソン効果と大幅な体積収縮を引き起こす。
  • 14本鎖 vs. 8本鎖： 14本鎖モデルは、横方向のロッドが座屈する明確な分岐を示し、これにより突然の軟化と劇的な体積減少が生じる。これは、座屈のない8本鎖の経路で見られる単調な体積増加とは対照的である。
• 圧縮：
  • 両方のRVEは、線形領域の後に（繊維の座屈による）プラトー（軟化）を示す。
  • 14本鎖モデルは、三線形応答を示す：初期の線形領域、第一プラトー（直線のロッドの座屈）、および第二プラトー（傾斜したロッドの座屈）。
  • モデルは、フィブリンネットワークで見られるフォーム状の挙動（ヤング率と臨界応力の体積分率に対するスケーリング）を再現している。
• 単純せん断：
  • 両モデルは、二線形応力-歪み応答を示す。分岐点に達すると、体対角線のロッドが座屈し、応力の軟化を引き起こす。
  • 顕著な逆ポインティング効果が観察される。ロッドが座屈して法線方向に収縮するにつれて、周期境界条件（$F_{11}=1$）を維持するために引張法線応力が生成される。この法線応力の大きさは、せん断応力に匹敵する。

意義と主張
本論文の主要な意義は、微視的な繊維力学とマクロな組織応答の間のギャップを埋める、厳密で不完全性を排除した計算ホモジナイゼーション・フレームワークを提供することにある。幾何学的に厳密なロッドを用いることで、本モデルは、単純なモデルでは見落とされがちな、バイオポリマーの挙挙動に不可欠な大規模変形と安定性遷移（座屈）を捉えることができる。著者らは、本アプローチが、現象論的なフィッティングパラメータを導入することなく、繊維の曲げと座屈のメカニズムのみを通じて、体積収縮や逆ポインティング効果のような特異な現象を回収できることを強調している。さらに、コラーゲンや血栓の実験データに対する検証成功、およびヘリカル・メタマテリアルのモデリング能力は、既存の生物物理学的現象の理解および将来のバイオインスパイアード構造材料の設計における本フレームワークの有用性を示唆している。著者らは、現在のモデルは純粋に弾性かつ局所的であり、高密度なネットワークの挙限的な挙動を完全に捉えるには、将来的に非弾性、繊維接触、および非局所的効果を組み込む必要があることを控えめに述べている。