🧬 biology

Non-linear visco-elasto-plastic rheology of a viscous vertex model

该研究通过构建包含显式粘性摩擦及主动各向异性应力的平均场模型，建立了顶点模型中细胞尺度与组织尺度力学之间的联系，成功推导并验证了能够描述细胞形状及非线性粘弹塑性流变行为的本构关系。

原作者： Shalabh Kumar Anand, Matthias Merkel

发布于 2026-02-27

📖 1 分钟阅读☕ 轻松阅读

原作者： Shalabh Kumar Anand, Matthias Merkel

原始论文采用 CC BY 4.0 许可（http://creativecommons.org/licenses/by/4.0/）。 ⚕️ 这是一篇未经同行评审的预印本的AI生成解释。这不是医疗建议。请勿根据此内容做出健康决定。阅读完整免责声明

这篇论文就像是在给生物组织（比如皮肤、胚胎或肿瘤）做“体检”和“性格分析”。

想象一下，生物组织不是像石头那样硬邦邦的，也不是像水那样完全流动的，它更像是一种**“有记忆的、会流动的果冻”**。科学家想知道：当你用力拉扯这块“果冻”时，它是怎么变形的？是像橡皮筋一样弹回去，还是像蜂蜜一样慢慢流走，或者像面团一样被揉坏了？

这篇文章的核心就是建立了一个数学公式（本构关系），用来预测这种“生物果冻”在受力时的各种反应。

为了让你更容易理解，我们可以用以下几个生动的比喻来拆解这篇论文：

1. 细胞就像“拥挤的泡沫”

想象一块充满了泡泡的肥皂泡（或者一锅煮得稠稠的粥）。

顶点模型（Vertex Model）： 科学家把每个细胞看作一个多边形（像六边形的蜂巢），它们挤在一起。这些多边形的角（顶点）就是细胞之间的连接点。
以前的研究： 以前大家研究这个模型时，通常假设这些泡泡是粘在桌子上的（有外部摩擦力），或者只关注它们怎么弹性变形（像弹簧）。
这篇论文的创新： 作者假设这些泡泡是自由漂浮的（没有粘在桌子上），它们之间的摩擦只发生在泡泡互相挤压、滑动的时候。这更像真实的生物组织（比如胚胎内部），没有外部桌子，只有细胞之间的内部摩擦。

2. 三种“变形”的魔法

当你对这块“细胞果冻”施加剪切力（比如横向推它）时，它会发生三种不同的变化，作者把它们分得很清楚：

A. 弹性变形（像拉橡皮筋）：
细胞被拉长或压扁了，形状变了，但细胞之间没有交换位置。如果你松手，它们会弹回去。
- 比喻： 就像你捏了一下橡皮泥，它凹下去了，但还没被揉散。
B. 粘性流动（像倒蜂蜜）：
细胞在内部摩擦中慢慢滑动，消耗了能量。
- 比喻： 就像你在蜂蜜里搅动勺子，蜂蜜会慢慢流动并产生热量。
C. 塑性重组（像换座位）：
这是最关键的！当力大到一定程度，细胞之间的连接断开，重新连接（T1 跃迁）。细胞 A 和 B 分开了，A 和 C 连上了。这种变化是不可逆的，就像把一锅粥里的米粒重新排列，你没法让它们自动回到原来的位置。
- 比喻： 就像在一个拥挤的地铁车厢里，大家互相推挤，最后有人换了个位置站。一旦换了，你就回不去了。

3. 核心发现：粘度会“变脸”

这是这篇论文最有趣的地方。

以前的想法： 大家可能认为组织的“粘度”（流动的难易程度）是一个固定的数，就像水的粘度是固定的。
这篇论文的发现： 组织的粘度是看情况的！
- 如果细胞只是弹性变形（像拉橡皮筋），粘度比较小。
- 如果细胞开始重组换位置（像地铁换座位），粘度会突然变大（变成原来的 3 倍！）。
- 比喻： 想象你在推一堵墙。如果墙只是微微晃动（弹性），推起来很轻松；但如果墙里的砖块开始互相摩擦、移位（重组），你会觉得突然变得非常沉重，推起来更费劲。

4. 细胞也会“主动捣乱”（活性力）

生物组织里的细胞不是死物，它们会自己用力（比如像肌肉一样收缩，或者像虫子一样爬行）。

作者研究了两种细胞“主动发力”的方式：
1. 各向异性张力： 细胞像拉紧的橡皮筋，只在一个方向上用力。
2. 爬行力： 细胞像小虫子一样，试图把自己拉向某个方向。
结果： 这些主动的力会直接改变组织的“性格”。比如，如果细胞拼命往一个方向爬，组织就会更容易向那个方向流动，甚至改变它变形的形状。

5. 为什么要这么做？（实际应用）

作者建立了一套**“万能公式”**。

以前，科学家想预测组织怎么变形，得在电脑里模拟成千上万个细胞，非常慢且复杂。
现在，作者把复杂的细胞互动简化成了几个简单的参数（形状、流动速度、重组频率）。
好处： 就像天气预报一样，你不需要模拟每一滴雨，只需要知道气压和温度，就能预测明天会不会下雨。现在，医生或生物学家可以用这个公式，直接预测胚胎发育时组织怎么折叠，或者肿瘤细胞怎么在体内扩散，而不需要每次都做超级复杂的模拟。

总结

这篇论文就像给生物组织画了一张**“性格地图”**：

它告诉我们，生物组织在受力时，既有弹性（能回弹），又有粘性（会流动），还有塑性（会永久变形）。
它发现，当细胞开始“换座位”（重组）时，组织会变得更粘稠，更难推动。
它证明了，如果细胞自己“想动”（主动力），会极大地改变组织的流动方式。

这就好比我们终于搞懂了：为什么揉面团时，刚开始很软，揉久了（细胞重组）反而觉得阻力变大了？因为面团里的“细胞”在重新排列，产生了新的内部摩擦。这套理论让科学家能更精准地预测生命体是如何从一团细胞变成复杂的人体的。

这篇论文题为《非线性的粘弹塑性顶点模型流变学》（Non-linear visco-elasto-plastic rheology of a viscous vertex model），由 Shalabh K. Anand 和 Matthias Merkel 撰写。文章旨在建立生物组织细胞尺度力学与组织尺度流变学之间的联系，通过构建平均场流变关系，描述顶点模型（Vertex Model）在非线性、粘弹塑性及活性条件下的力学行为。

以下是对该论文的详细技术总结：

1. 研究背景与问题 (Problem)

核心挑战：理解生物组织的形态发生（Morphogenesis）需要解析组织尺度的力学行为如何从大量细胞的集体行为中涌现。虽然已有许多关于顶点模型线性流变学和弹性性质的研究，但关于应力耗散（stress dissipation）和非线性粘弹塑性机制（non-linear visco-plastic regime）的研究尚不充分。
现有模型的局限：
- 许多顶点模型使用基底摩擦（substrate friction），这不符合无基底组织（如 3D 干细胞聚集体）的物理现实。
- 现有的内部摩擦模型往往缺乏伽利略不变性（Galilean invariance），即无法保证组织整体平移或旋转时细胞受力不变。
- 缺乏一个统一的平均场框架来同时描述细胞形状变化、细胞重排（T1 跃迁）以及活性应力对非线性流变的影响。

2. 方法论 (Methodology)

作者构建了一个包含显式粘性摩擦和两种不同活性机制的顶点模型，并采用平均场理论进行推导。

模型构建：
- 能量函数：采用经典的顶点模型能量函数，包含细胞面积弹性（ $k_A$ ）和周长弹性（ $p_0$ ）。
- 显式粘性摩擦：为了满足平移和旋转的伽利略不变性，作者在细胞面积、细胞周长和细胞间界面长度上引入了标量摩擦项。摩擦力取决于速度差，而非绝对速度。
- 活性机制：引入了两种产生各向异性主动应力的机制：
  1. 各向异性界面张力（Anisotropic interface tension）：依赖于界面角度的主动张力。
  2. 主动爬行力（Active crawling forces）：模拟类似爪蟾胚胎中的收敛延伸机制，在细胞内施加拉力。
- 动力学方程：在过阻尼极限下，通过力平衡方程求解顶点运动，包含弹性力、摩擦力和主动力。
应变分解框架：
- 利用基于三角形剖分的框架，将组织剪切率 $\tilde{v}$ 分解为两部分：
  $\tilde{v} = \frac{DQ}{Dt} + R$
  其中， $Q$ 是细胞形状各向异性张量（可逆部分）， $R$ 是由 T1 跃迁和关联效应引起的重排率（不可逆部分）。
- 这种分解允许将流变问题拆分为两个独立的关系：
  1. 剪切应力与细胞形状及剪切率的关系： $\tilde{\sigma}(Q, \tilde{v})$ 。
  2. 重排率与细胞形状及剪切率的关系： $R(Q, \tilde{v})$ 。
模拟方案：
- 进行了小振幅振荡剪切（线性流变）、恒定剪切率（非线性弹性与粘性）、以及大振幅振荡剪切（验证模型预测）的模拟。

3. 关键贡献与主要结果 (Key Contributions & Results)

A. 剪切应力的本构关系 $\tilde{\sigma}(Q, \tilde{v})$

线性粘弹区：在小振幅振荡下，模型表现出恒定的剪切模量 $G_0$ 和粘度 $\eta_0$ 。
非线性弹性区：
- 剪切应力与细胞形状 $Q$ 呈非线性关系。在 $Q \approx 0.4$ 范围内，可用三阶多项式拟合： $\tilde{\sigma}_{el} \approx 2G_0 Q + 4G_2 Q^3$ 。
- 当 $Q > 0.5$ 时，需要更高阶（五阶）多项式来准确描述弹性应力。
粘性与重排的耦合：
- 研究发现，当组织发生 T1 跃迁（重排）时，有效粘度显著增加（约为仅发生弹性变形时的 3 倍）。
- 这表明粘度 $\eta$ 并非恒定，而是依赖于是否存在不可逆重排。作者提出了一个插值公式，利用不可逆分数 $R/\tilde{v}$ 来连接小振幅粘度 $\eta_0$ 和稳态剪切粘度 $\eta_s$ 。
活性应力的影响：
- 两种活性机制（各向异性张力和爬行力）均导致剪切应力曲线发生线性偏移。
- 建立了活性修正项： $\tilde{\sigma} = \tilde{\sigma}_{passive} - \rho T^{act}$ ，其中 $\rho$ 为活性系数。

B. 重排率的本构关系 $R(Q, \tilde{v})$

不可逆分数：定义了不可逆分数 $R/\tilde{v}$ ，它描述了剪切中有多少比例是由细胞重排完成的。
指数依赖：发现不可逆分数随细胞形状 $Q$ 呈指数增长：
$\frac{R}{\tilde{v}} = \frac{1}{2} \exp\left( b(\tilde{v}) [Q - Q_0(\tilde{v})] \right)$
活性效应：活性力主要改变参数 $Q_0$ （即屈服发生的临界形状），而对斜率参数 $b$ 影响较小。爬行力比各向异性张力更能显著降低 $Q_0$ ，导致在较小的剪切率下即可发生重排。

C. 模型验证与预测

稳态预测：利用推导出的 $R(Q, \tilde{v})$ 关系，成功预测了不同剪切率下的稳态细胞形状 $Q_s$ ，与模拟数据高度吻合。
大振幅振荡剪切（LAOS）：将模型应用于大振幅振荡剪切，能够准确预测细胞形状 $Q(t)$ $Q (t)$ 和剪切应力 $\tilde{\sigma}(t)$ $\tilde{σ} (t)$ 的时间演化轨迹。
- 在屈服点附近（剪切方向反转前）存在轻微偏差，这归因于模型未完全捕捉变形历史的记忆效应。

4. 意义与讨论 (Significance)

理论突破：这是首次系统地探测具有内部摩擦（满足伽利略不变性）的顶点模型的非线性流变学。这解决了以往模型中基底摩擦可能掩盖活性物质不稳定性的问题。
机制解析：明确了组织粘度的来源不仅来自细胞形状的弹性响应，还强烈依赖于 T1 跃迁引起的耗散。当发生重排时，粘度显著增加，这解释了组织在塑性流动中的硬化或耗散行为。
通用性：提出的平均场方法（将应力分解为形状依赖和重排依赖）具有通用性，可推广至其他基于细胞的组织模型。
实验指导：模型中引入的细胞形状 $Q$ 作为流变方程的变量，使得该理论更容易与实验观测（如显微镜下的细胞形态）相结合，为定量理解发育生物学和癌症转移中的组织力学提供了工具。

总结

该论文通过引入满足物理对称性的内部摩擦和多种活性机制，建立了一个全面的非线性粘弹塑性平均场流变模型。该模型成功地将微观的细胞重排（T1 跃迁）和细胞形状变化与宏观的组织应力响应联系起来，揭示了粘度对重排状态的依赖性，并为预测复杂变形下的组织行为提供了可靠的理论框架。