🔬 condensed matter

Reversible viscoelasticity and irreversible elastoplasticity in the power law creep and yielding of gels and fibre network materials under stress

この計算科学的研究は、アサーマルな繊維ネットワークが、一定の応力下において、限界連結性とフィラメントの破断がべき乗則に従う変形ダイナミクスを駆動し、それが最終的に材料の形状回復能または壊滅的な破壊を決定づける、明確に区別される可逆的な粘弾性クリープおよび不可逆的な弾塑性クリープ領域を示すことを明らかにしている。

原著者： Michael J. Hertaeg, Suzanne M. Fielding

公開日 2026-01-23

📖 1 分で読めます☕ さくっと読める

CC BY 4.0

原著者： Michael J. Hertaeg, Suzanne M. Fielding

原論文は CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/) でライセンスされています。 ✨ これは以下の論文のAI生成解説です。著者が執筆または承認したものではありません。技術的な正確性については原論文を参照してください。免責事項の全文を読む

ゲルやフェルトのような素材を、単なる固形ブロックとしてではなく、小さくて伸縮性のある糸でできた巨大で乱雑なクモの巣として想像してみてください。この論文は、コンピュータ・シミュレーションを用いて次のような問いを投げかけています。「一定の力でこのクモの巣を引っ張り、その状態を維持し続けたら、何が起こるのか？」

研究者たちは、主に2つのことを理解したいと考えました。

クモの巣は時間の経過とともにどのように伸びるのか？
最終的にいつ、パチンと切れて壊れるのか？

以下に、簡単な比喩を用いた彼らの知見のまとめを記します。

設定：糸のネットワーク

科学者たちは、これらの糸で満たされた2次元の世界のデジタルモデルを構築しました。

糸（Strings）: これらはゴムバンドのように機能します。少し伸ばすと引き戻そうとし、伸ばしすぎると（特定の「破断点」を超えると）、永久に切れてしまいます。
結合（Connections）: 糸は結び目（ノット）でつながっています。各結び目に結びついている糸の数を「連結性（ $Z$ $Z$ ）」と呼びます。
- 密なウェブ ( $Z > 4$ ): 結び目あたりの糸が多い。このウェブは硬く、強いです。
- 疎なウェブ ( $Z < 4$ ): 結び目あたりの糸が少ない。このウェブは「ぐにゃぐにゃ」としていて、ふらふらしています。
- ちょうど良いウェブ ( $Z = 4$ ): これは「マージナル（境界的）」な点であり、ウェブが辛うじて自立して形を保っている状態です。

実験：一定の牽引

彼らはウェブに対して、一定で変化しない引っ張り（せん断応力）を与え、その経過を観察しました。彼らはウェブがどれくらい伸びたか（歪み）と、どのくらいの速さで伸びているか（歪み速度）を測定しました。

知見1：「ちょうど良い」ウェブと魔法の伸び

ウェブが硬い（結合が多い）場合、または非常に疎な（結合が少ない）場合、伸び方は予測可能なものでした。

硬いウェブ: 少し伸びた後、停止し、安定した形状に落ち着きました。
疎なウェブ: ふらつきながら大きく伸び、その後も落ち着きました。

しかし、「ちょうど良い」ウェブ ( $Z=4$ ) は特別でした。
このマージナルな連結を持つウェブに力を加えると、単に止まるだけではありませんでした。それは非常に特定の、数学的なリズムに従って伸び続けました。伸びる速度は時間の経過とともに減衰し、「べき乗則（パワーロー）」に従いました。

比喩: 人が廊下を歩いている様子を想像してください。普通の人は速く歩いたり、その後遅くなったり、止まったりします。この「ちょうど良い」ウェブは、まるで永遠に歩き続ける人が、完璧で予測可能なパターンに従って、一歩一歩の歩幅を小さくしていくようなものです。研究者たちは、これが起こるのは、まだ何かが壊れているからではなく、ウェブが内部の結び目を複雑かつ不均一に再配置しているためであることを突き止めました。

知見2：じわじわとした漏れと突然の破断

研究者たちは、糸が引き伸ばされすぎた場合に切れる設定を導入しました。これにより、すべてが変わりました。

じわじわとした漏れ（弾塑性クリープ）: 本来であれば（このような現象を示さないはずの）「疎な」ウェブであっても、糸が切れる設定にすると、ゆっくりとした一定の伸びが生じるようになりました。ウェブが伸びるにつれて、いくつかの糸が切れました。これがウェブを損傷させ、弱体化させ、それがさらなる伸びを引き起こし、さらに多くの糸を切らせるという連鎖が起きました。
- 比喩: ゼムクリップを何度も曲げ伸ばししている様子を考えてみください。最初は簡単に曲がります。しかし、数回曲げた後、小さな亀裂が生じます。その亀裂があるために次の曲げがより容易になり、それがさらに大きな亀裂を生みます。最終的に、全体が折れます。シミュレーションにおけるこの「亀裂」のフェーズは非常に長く続くことがあり、これによって、前述の「ちょうど良い」ウェブで見られたものと同じ「べき乗則」の伸びが発生しました。ただし、今回は「再配置」ではなく「損傷」によって引き起こされていました。
突然の破断（降伏）: 引っ張る力が十分に強ければ、ウェブは最終的に「転換点」に達します。長い時間のゆっくりとした伸びと損傷の蓄積の後、ウェブは突如として壊滅的な失敗（崩壊）を迎えます。
- 遅延: 引っ張る力が強いほど、破断までの時間は早まりました。しかし、力が破断点よりもわずかに低い程度であれば、ウェブは最終的に屈するまで、驚くほど長い時間持ちこたえることができました。研究者たちは、この「待ち時間」が、どの程度の力で引いているか、および個々の糸の強さに基づいてどのように変化するかを正確に描き出しました。

知見3：ウェブは回復できるか？（リカバリー）

最後に、彼らは問いかけました。「もし引っ張るのをやめたら、ウェブは元の形に戻るのか？」

弾性的な部分: ウェブが硬く、糸がまだ切れていなければ、それはゴムバンドのようでした。手を離すと、ウェブは完璧に元の形に戻りました。
塑性的な部分: ウェブが疎であったり、すでにいくつかの糸が切れていたりする場合、それは粘土のようでした。手を離しても、ウェブは伸びたままの状態でした。損傷（切れた糸）は永久的なものであるため、完全には回復できませんでした。
驚きの発見: 硬いウェブであっても、いくつかの糸が切れてしまった場合、ウェブは完全には回復できませんでした。しかし、ウェブの接続が非常に強固で、いくつかの糸が切れた後でもなお「硬い」状態を維持できるほどであれば、ほぼ完璧に元に戻ることができました。これは橋のようなものです。ボルトを数本失ったとしても、形を保ち続けることができますが、失いすぎると、永久にたわんでしまいます。

まとめ

この論文は、なぜ一部の軟質材料（ゲルや生物学的組織など）が、壊れる前に長い間、ゆっくりと予測可能な形で伸び続けるのかを説明しています。

メカニズムA: 材料が完璧なバランス（ちょうど良い接続数）である場合、内部の再配置によってゆっくりと伸びます。
メカニズムB: 材料がより弱い、あるいは力が強い場合、自身をゆっくりと破壊していくことによってゆっくりと伸びます。

どちらのメカニズムも同じ「べき乗則」の伸び方をもたらしますが、一方は可逆的（弾性的）であり、もう一方は永久的な損傷（塑性的）を伴います。この研究は、これらの材料がなぜ一定の荷重の下で長時間持ちこたえ、そして突如として失敗するのかを解明する助けとなります。

技術要約：ゲルおよび繊維ネットワーク材料におけるべき乗則クリープと降伏における可逆的粘弾性と不可逆的エラストプラスティシティ

問題提起
生物学的細胞骨格からハイドロゲル、紙、コロイドゲルに至るネットワーク材料は、一定の応力下で、歪み速度が $\dot{\gamma} \sim t^n$ ( $n < 0$ ) と減衰し、歪みが劣線形に増大する「べき乗則クリープ」を頻繁に示す。このような挙動はソフトマターにおいて普遍的であるが、その根底にある物理的メカニズムについては依然として議論の余地がある。具体的には、べき乗則クリープが、純粋に可逆的な粘弾性変形（非アフィン再配列）に由来するのか、それとも不可逆的なエラストプラスティックな損傷（結合の破断）に由来するのかが不明である。さらに、これらの材料が緩やかなクリープから破滅的な降伏および破壊へと移行する条件、および除荷時に元の形状をどの程度回復するかについても、体系的な特性評価が必要である。本研究では、非熱的な繊維ネットワークの最小計算モデルを用いて、これらの問いに取り組む。

手法
著者らは、 $d=2$ 次空間における最小非熱的無秩序ネットワークモデルの数値シミュレーションを用いている。

モデル構築: 二分散ディスクをランダムに充填し、重なり合う中心同士をHookeanスプリングで接続し、ターゲットとする平均配位数 $Z$ を達成するために反復的に枝刈りを行うことで、ネットワークを生成する。各スプリングは破断歪み $b$ まで線形弾性的であり、それを超えると不可逆的に破断する。
レオロジー・プロトコル: 本研究では、シミュレーションにおいて定歪速度プロトコルよりも一般的ではない、一定の印加剪断応力 $\Sigma$ プロトコルを利用している。
数値アルゴリズム: 一定の応力を維持するために、著者らは2つのアルゴリズムを開発し、検証した。アルゴリズムIは、ニュートン・ラフソン法を用いて歪み速度を反復的に調整し、目標の応力に一致させる。アルゴリズム形式IIは、仮想的なニュートン溶媒の粘度 $\eta$ を用いて、歪み速度を $\dot{\gamma} = (\Sigma - \sigma)/\eta$ と設定する。著者らは、 $\eta \to 0$ の極限において両アルゴリズムが同一の結果に収束することを示し、その後、より効率的なアルゴリズムIIを $\eta = 0.16$ で使用した。
パラメータ: シミュレーションは $N=16,000$ 個のノードを用いて行われる。主要な変数には、ネットワークの連結性 $Z$ 、フィラメント破断閾値 $b$ 、および印加応力 $\Sigma$ が含まれる。

主な貢献と結果

純粋な粘弾性領域（フィラメント破断なし）:
- 硬いネットワーク ( $Z > 4$ ): マクスウェルの剛性閾値（ $Z > 2d = 4$ ）を超えるネットワークは、長時間経過後にHookean固体として振る舞い、定常歪み $\gamma_{ss} \propto \Sigma$ に達する。この過渡的なクリープは、非アフィンな内部変形によって駆動される純粋な粘弾性挙動である。
- 柔軟なネットワーク ( $Z < 4$ ): 剛性閾値を下回るネットワークは、「フロッピー（柔軟）」な応答を示し、極めて微小な応力によって有限の歪みが生じる。歪み速度は $\dot{\gamma} \sim t^{-2}$ として減衰し、歪みは一定値に飽和する。
- 限界連結性 ( $Z \approx 4$ ): 限界連結に近いネットワーク（ $Z \gtrsim 4$ ）において、特異な領域が現れる。ここでは、システムが結合破断を起こすことなく、長期間にわたって $\dot{\gamma} \sim t^{-0.5}$ および $\gamma \sim t^{0.5}$ という持続的なべき乗則クリープを示し、その後飽和する。この挙動は純粋に粘弾性的であり、結合破断のない非アフィン変形に由来する。
エラストプラスティック領域（フィラメント破断あり）:
- 劣剛性ゲルのべき乗則クリープ: フィラメントの破断を許容する場合（ $b < \infty$ ）、べき乗則クリープを支える条件の範囲が拡大する。具体的には、純粋な粘弾性の場合には急速な飽和を示すはずの、初期状態が柔軟なネットワーク（ $Z < 4$ ）が、プラスチックな損傷の蓄積によって、持続的なべき乗則クリープ（ $-1 < n < 0$ の $\dot{\gamma} \sim t^n$ ）を示すようになる。
- 降伏と破壊: 応力が臨界降伏応力 $\Sigma_y$ を超えると、ネットワークは遅延時間 $\tau_y$ の後に破滅的な破壊を起こす。降伏応力は、破断閾値 $b$ に対してほぼ線形にスケールする。遅延時間 $\tau_y$ は、 $\Sigma$ が $\Sigma_y$ に対して上方から接近するにつれて劇的に増加する。
- 検証: 応力印加型シミュレーションで測定された降伏応力は、歪み印加型のスターアップ・シミュレーションにおける歪み速度ゼロの極限と一致しており、測定された降伏点の堅牢性を裏付けている。
歪み回復:
- 本研究では、除荷時の回復可能な歪みと回復不可能な歪みを定量化している。
- 硬いネットワーク ( $Z > 4$ ): 有意なプラスチック損傷が蓄積する前に応力が除去されれば、完全な形状回復が起こる。驚くべきことに、少数のスプリングが破断したとしても、残存する連結性が剛性閾値（ $Z > 4$ ）を維持していれば、ネットワークは完全に形状を回復できる。
- 柔軟なネットワーク ( $Z < 4$ ): これらのネットワークは、たとえ微小な歪みであっても部分的または全く回復を示さない。これは、反動に必要な弾性バックボーンを欠いているためである。結合破断がほとんどなくても、荷重時には直ちに回復不可能な歪みが生じる。
損傷の局在化:
- 破壊モードの視覚的解析により、明確なパターンが明らかになった。初期に硬いネットワークでは、プラスチックな損傷はシステム全体に広がるマクロなバンドへと局在化する。対照的に、初期に柔軟なネットワークにおける損傷は、材料全体に拡散的に分布する。

意義と主張
本論文は、ネットワーク材料におけるべき乗則クリープを生み出すことが可能な、2つの異なる物理的メカニズムを特定したと主張している：

粘弾性メカニズム: 損傷を伴わない非アフィン変形に起因する、限界連結ネットワーク（ $Z \gtrsim 4$ ）における挙動。
エラストプラスティック・メカニズム: 結合破断の漸進的な蓄積に起因する、劣剛性ネットワーク（ $Z < 4$ ）における挙動。

著者らは、これら2つのメカニズムが彼らの最小モデルにおけるべき乗則クリープを説明するものであると結論付けているが、これらが多様な実験的ソフトマターで見られる現象を説明する「唯一の」メカニズムであるとは主張していない。彼らは、Rouse的な鎖の運動、フラクタルなネットワーク構造、あるいは熱的なリンカーのダイナミクスといった他の要因も、現実世界のシステムにおけるべき乗則クリープトの普遍性に寄与している可能性があることを控えめに示唆している。本研究は、クリープおよび降伏の文脈において、可逆的な粘弾性と不可逆的なエラストプラスティシティを区別するための計算論的枠組みを提供するものである。