try ai
科普
编辑
分享
反馈
  • Y.C. Fung的准线性粘弹性(QLV)理论

Y.C. Fung的准线性粘弹性(QLV)理论

SciencePedia玻尔百科
关键要点
  • 准线性粘弹性(QLV)理论通过将软组织的行为分离为一个非线性的瞬时弹性响应和一个线性的时间依赖性松弛来对其进行建模。
  • Fung的指数定律描述了特征性的应变硬化(“J”形)弹性行为,而一个折减松弛函数G(t)则描述了材料的衰退记忆。
  • QLV模型通过应力松弛实验得到验证,在这些实验中,来自不同应变水平的归一化应力-时间曲线会塌缩到一条单一的主曲线上。
  • QLV原理应用广泛,解释了肌腱和肺(片流)的功能、残余应力的重要性,并为电池设计等工程领域提供了信息。

引言

从肌腱到皮肤,生物材料表现出一种复杂的力学行为,它既像弹簧一样具有弹性,又像蜂蜜一样具有粘性。这种被称为粘弹性的独特组合,包含了应变硬化和时间依赖性应力松弛等特性,是简单的物理定律无法捕捉的。本文要解决的核心挑战是如何用一种既准确又优雅的方式来数学地模拟这种栩栩如生的行为。这正是生物工程先驱Y.C. Fung通过他的准线性粘弹性(QLV)理论解决的问题。本文将引导您了解这一里程碑式的理论。第一章“原理与机制”将解构QLV模型,解释其核心的可分离性假设、数学表述和实验验证。随后,“应用与跨学科联系”将探讨这些思想的深远影响,展示它们如何解释我们组织的功能,甚至为现代工程挑战提供启示。让我们从探索Fung模型背后卓越的物理直觉开始。

原理与机制

如果你曾拉伸过一根橡皮筋,你就会知道它会回弹。你拉得越用力,它的回弹力就越大。这就是弹性的本质,一个被Hooke定律为简单弹簧所优雅描述的概念。如果你曾用勺子在蜂蜜中推动,你就会知道它会抵抗运动,而且你推得越快,它的抵抗力就越大。这就是粘性。但构成我们身体的材料又如何呢?一根韧带、一根肌腱或一块皮肤,既不是一个简单的弹簧,也不是一罐简单的蜂蜜。它是两者的精湛结合,甚至远不止于此。它会伸展,但它所施加的力不仅取决于它被拉伸了多少,还取决于它如何达到这个状态以及它被保持了多长时间。它有记忆。此外,与简单弹簧不同,你拉得越用力,它就变得越硬。我们如何才能将如此复杂、栩栩如生的行为用一个既优雅又有用的数学定律来捕捉呢?

这正是伟大的生物工程师Y.C. Fung所接受的挑战。他的解决方案,一个被称为​​准线性粘弹性(QLV)​​的理论,是物理直觉的杰作。Fung的核心思想是将材料视为具有两种截然不同的“个性”,并且在很大程度上,这两种个性可以被分开考虑:一种是瞬时的、非线性的弹性个性,另一种是时间依赖的、会遗忘的个性。

组织的两种个性

想象一下,你能够以极快的速度拉动一根韧带,以至于它没有时间进行内部流动或重排。在这个稍纵即逝的瞬间,你将目睹其纯粹的​​瞬时弹性响应​​,我们称之为σe\sigma_eσe​。这并非简单线性弹簧的响应。像肌腱这样的生物组织表现出一种称为​​应变硬化​​的显著特性:随着拉伸,它们会逐渐变硬。这是因为在小应变下,你主要只是在拉直微观上卷曲的胶原纤维。一旦它们绷紧,你便开始拉伸纤维本身,这需要大得多的力。这种行为赋予了组织特有的“J”形应力-应变曲线。

Fung也为这种瞬时个性提出了一个优美且广泛使用的数学形式:

σe(ε)=A(exp⁡(Bε)−1)\sigma_e(\varepsilon) = A \left( \exp(B \varepsilon) - 1 \right)σe​(ε)=A(exp(Bε)−1)

这里,ε\varepsilonε是应变(长度的分数变化)。参数AAA设定了组织的整体应力尺度,而无量纲参数BBB则控制非线性程度。更大的BBB意味着组织在拉伸时会更急剧地变硬。这种指数形式优雅地捕捉了组织在开始时柔软顺从,但在大变形下变得坚韧和有抵抗力的本质,从而保护我们的关节免受伤害。

现在,第二种个性又是什么呢?它负责组织的​​衰退记忆​​。如果你将韧带拉伸到某个长度并保持不动,你会发现维持该长度所需的力会随时间慢慢减小。材料​​松弛​​了。这是它的粘性,或时间依赖性。Fung提出,这种行为可以用一个​​折减松弛函数​​G(t)G(t)G(t)来描述。这个函数代表了在时间ttt后,初始应力所剩下的分数。根据定义,在拉伸的瞬间(t=0t=0t=0),没有发生松弛,所以G(0)=1G(0)=1G(0)=1。随着时间的推移,材料“忘记”了初始应变,G(t)G(t)G(t)会衰减到某个较小的值。例如,G(t)G(t)G(t)的一个常见形式是衰减指数的和,如

G(t)=g1exp⁡(−t/τ1)+g2exp⁡(−t/τ2)G(t) = g_1 \exp(-t/\tau_1) + g_2 \exp(-t/\tau_2)G(t)=g1​exp(−t/τ1​)+g2​exp(−t/τ2​)

这可能代表了组织内两种不同的物理松弛过程:一个时间常数为τ1\tau_1τ1​的快速过程和一个时间常数为τ2\tau_2τ2​的慢速过程。

QLV中的“准线性”一词源于Fung的卓越假设:虽然弹性响应σe(ε)\sigma_e(\varepsilon)σe​(ε)是高度非线性的,但松弛行为G(t)G(t)G(t)却不是。他提出,松弛曲线G(t)G(t)G(t)的形状与应变的大小无关。换句话说,无论材料被拉伸了多少,其“遗忘”的方式都遵循相同的模式。这就是关键的​​可分离性假设​​。

记忆的叠加

我们现在有了两种个性:一种是非线性的弹性个性,另一种是线性的、会遗忘的个性。我们如何将它们结合起来,以预测任何任意应变历史下的应力呢?Fung的洞见是将物理学中一个强大的思想——​​Boltzmann叠加原理​​——应用于它所产生的瞬时弹性应力,而不是应变本身。

将一个连续的拉伸运动想象成一系列无穷小的、无限微小的步骤。在过去的某个时间τ\tauτ的每一个微小应变变化,都会产生一个微小的瞬时弹性应力增量dσed\sigma_edσe​。随着时间推移到当前时刻ttt,对那个微小应力增量的记忆逐渐消退。它对当前应力的贡献被一个因子G(t−τ)G(t-\tau)G(t−τ)所削弱,这是松弛函数在经过的时间内的值。

为了找到我们现在,在时间ttt所感受到的总应力,我们必须将所有过去无限小步骤中所有这些衰退记忆的贡献加起来(积分)。这一思路直接导出了QLV理论中优美而强大的遗传积分[@problem_id:4166247, @problem_id:4195204]:

σ(t)=∫0tG(t−τ)dσe(ε(τ))dτdτ\sigma(t) = \int_{0}^{t} G(t-\tau) \frac{d\sigma_{e}(\varepsilon(\tau))}{d\tau} d\tauσ(t)=∫0t​G(t−τ)dτdσe​(ε(τ))​dτ

这个方程是QLV的核心。它告诉我们,今天的应力是过去所有弹性应力变化的加权平均值,其中最近发生的事件比遥远过去的事件权重更大,其权重精确地由材料的“遗忘”函数G(t)G(t)G(t)决定。

质询材料:与实验的对话

一个理论,无论多么优雅,其价值在于它预测实验结果的能力。我们如何测试这种可分离个性的想法对于真实的肌腱是否成立呢?

最直接、最富启发性的实验是​​应力松弛实验​​。我们取一个组织样本,尽可能快地将其拉伸到固定的应变ε0\varepsilon_0ε0​,然后在保持应变恒定的同时测量应力σ(t)\sigma(t)σ(t)。

QLV积分对这种情况的预测是什么?应变历史是一个阶跃函数。瞬时弹性应力σe\sigma_eσe​也在t=0t=0t=0时跃升至一个值σe(ε0)\sigma_e(\varepsilon_0)σe​(ε0​)并保持不变。它的变化率在t=0t=0t=0处是一个尖峰(一个Dirac delta函数)。当我们进行积分时,复杂的积分坍缩成一个极其简单的结果[@problem_id:4195174, @problem_id:4201096]:

σ(t)=σe(ε0)G(t)\sigma(t) = \sigma_e(\varepsilon_0) G(t)σ(t)=σe​(ε0​)G(t)

这个预测意义深远。它表明,你随时间测量的整个应力曲线,仅仅是瞬时弹性应力(一个仅取决于应变大小ε0\varepsilon_0ε0​的数字)与普适的松弛函数G(t)G(t)G(t)(一个仅取决于时间的函数)的乘积。

这为我们提供了一个直接的实验方案,以验证QLV并独立测量这两种个性[@problem_id:4183753, @problem_id:4201073]:

  1. 在一系列不同的应变水平下进行应力松弛实验,比如ε0=0.02,0.04,0.06\varepsilon_0 = 0.02, 0.04, 0.06ε0​=0.02,0.04,0.06。
  2. 对于每个实验,测量初始峰值应力σ(0+)\sigma(0^+)σ(0+)。根据我们的理论,这个值就是瞬时弹性响应σe(ε0)\sigma_e(\varepsilon_0)σe​(ε0​)。将这些峰值应力与其对应的应变作图,我们就能直接测量出材料的非线性弹性个性。对于真实的韧带,我们会发现比值σ(0+)/ε0\sigma(0^+)/\varepsilon_0σ(0+)/ε0​不是常数,这证实了其非线性。
  3. 接下来,对于每个实验,将整个测量的应力曲线σ(t)\sigma(t)σ(t)用其自身的初始峰值σ(0+)\sigma(0^+)σ(0+)进行归一化。QLV预测这个归一化曲线等于G(t)G(t)G(t)。
σ(t)σ(0+)=σe(ε0)G(t)σe(ε0)=G(t)\frac{\sigma(t)}{\sigma(0^+)} = \frac{\sigma_e(\varepsilon_0) G(t)}{\sigma_e(\varepsilon_0)} = G(t)σ(0+)σ(t)​=σe​(ε0​)σe​(ε0​)G(t)​=G(t)

如果可分离性假设是正确的,那么所有这些来自不同应变水平的归一化曲线,都应该​​塌缩到一条单一的、普适的主曲线上​​。这条主曲线就是对材料“遗忘”个性G(t)G(t)G(t)的直接测量。许多软组织的实验数据精确地显示了这种塌缩现象,这是支持QLV模型最有力的证据。这种理论与实验之间的优美对话是科学方法的实践,它让我们能够用数学的语言向材料提问并理解它的答案。

这种物理思想的力量远不止于简单的一维拉伸。对于复杂的三维变形,我们可以用更普适的张量对应物替换简单的应力σ\sigmaσ和应变ε\varepsilonε,例如第二类Piola-Kirchhoff应力SSS和Green-Lagrange应变EEE。基本的QLV积分形式保持不变,展示了其基本原理的深层统一性。

优雅的局限

QLV是软组织力学的最终定论吗?当然不是。科学的进步在于建立模型,然后同样重要的是,发现它们的失效之处。QLV的优雅之处在于其可分离性假设——即材料的遗忘方式与它被拉伸的程度无关。但这总是正确的吗?

考虑一个循环加载实验,我们反复拉伸和释放一根肌腱。我们在拉伸过程中输入的部分能量在释放时并未恢复;它以热量的形式耗散,形成一个​​滞后回线​​。QLV模型可以预测这个回线的大小。由于其粘性个性是线性的,它预测耗散的能量应与应变振幅的平方(ε02\varepsilon_0^2ε02​)成正比,并且应力与应变之间的相位滞后应与此振幅无关。

然而,对肌腱的仔细实验表明,这并不完全正确。相位滞后通常随应变振幅的增加而增加,耗散的能量与振幅的幂次小于2(例如,ε01.6\varepsilon_0^{1.6}ε01.6​)成比例。这告诉我们,可分离个性的简单图景是一个近似。实际上,弹性和粘性行为是微妙耦合的。更大的应变可以改变组织的微观结构——也许是通过改变流体流动路径或纤维相互作用——这反过来又改变了松弛行为。松弛函数GGG并非完全独立于应变。

发现这些局限性并非QLV模型的失败,而是其成功的证明。它为广泛的行为提供了极其准确的描述,并通过精确识别其与现实偏离之处,揭示了其中更微妙的物理机制,并引导科学家开发更精炼、完全非线性的粘弹性模型。Y.C. Fung的准线性理论仍然是生物力学的基石,是一个光辉的典范,展示了一个简单、直观的物理思想如何能为一个复杂的世界带来清晰和秩序。

应用与跨学科联系

在走过生物力学基本原理的旅程后,我们可能倾向于将它们视为优雅但抽象的数学描述。事实远非如此。真正的魔力始于我们用这些原理作为透镜来观察世界之时。我们发现,同一套思想可以解释短跑运动员肌腱非凡的韧性、我们呼吸的静谧效率,甚至我们手机中电池的退化。这是对一个物理理论的终极考验:不仅是其内部的一致性,还有它连接看似不相关现象的力量。在本章中,我们将探索这个联系之网,追随Y.C. Fung从构成我们的组织到塑造我们未来的技术的开创性工作足迹。

生命组织的力学交响曲

让我们从你身体里能感觉到的一样东西开始:肌腱,比如你脚后跟的跟腱。当你拉动它时,它的行为不像一个简单的弹簧。起初,它相当柔软,容易拉伸。但当你更用力地拉时,它会迅速变得异常坚硬。这种特有的“J”形应力-应变曲线并非偶然;它是生物工程的杰作。在微观层面上,肌腱由波浪状的胶原原纤维组成。在最初的低力拉伸——“足趾区”——我们主要只是在拉直这些卷曲的原纤维。这几乎不需要力。随着越来越多的原纤维解开卷曲并与拉力方向对齐,它们开始分担负载。这种渐进式的“募集”导致组织急剧变硬,进入一个“准线性”区域,此时它的行为像一个非常硬的弹簧。这种设计非常巧妙:它为小范围运动提供了灵活性,但在施加高力时又提供了难以置信的强度和抗损伤能力。

但故事并不仅仅止于简单的弹性。生命组织是湿润、复杂的复合材料。它们会随时间流动和蠕变。如果你拉伸一块软骨并保持不动,维持该拉伸所需的力会慢慢减小。这就是应力松弛,粘弹性的一个标志。要描述这一点,我们需要比简单弹性更复杂的定律。这正是准线性粘弹性(QLV)理论天才之处。该理论的核心思想,即可分离性假设,是将组织的响应视为两个独立部分的乘积:一个捕捉瞬时J形响应的非线性弹性部分,和一个描述应力如何随时间松弛的时间依赖性部分。

在数学上,这通过一个遗传积分优美地表达出来。对于给定的应变历史ε(t)\varepsilon(t)ε(t),应力σ(t)\sigma(t)σ(t)是通过对所有过去变化的影响进行积分得到的。QLV模型提出了这个优雅的形式: σ(t)=∫0tG(t−τ) dσe(ε(τ))dτ dτ\sigma(t) = \int_{0}^{t} G(t-\tau) \, \frac{d\sigma_e(\varepsilon(\tau))}{d\tau} \, d\tauσ(t)=∫0t​G(t−τ)dτdσe​(ε(τ))​dτ 在这里,σe(ε)\sigma_e(\varepsilon)σe​(ε)是我们为肌腱讨论过的瞬时弹性应力,它捕捉了非线性。函数G(t)G(t)G(t)是“折减松弛函数”;它是一个归一化的、与应变无关的函数,描述了材料的“记忆”——过去事件产生的应力如何随时间消退。这种分离使我们能够用两个更简单、独立的函数来模拟复杂的行为,这是一种极其强大的简化。

仿佛这还不够,生物组织还隐藏着另一个秘密。如果你小心地取出一个生物结构,比如一根动脉或一片筋膜,然后把它放在桌子上,它是“无载荷”的——没有外力作用于它。然而,它并非“无应力”。如果你切一刀,它会弹开,暴露出一种内在的张力。这就是​​残余应力​​。它之所以产生,是因为组织在负载下生长和重塑,创造了一个不相容的无应力状态。想象一下建造一座石拱;单个的石头是无应力的,但一旦组装成拱形,它们就处于恒定的压缩状态以保持其形状。类似地,血管的生长方式使得即使在没有血压的情况下,内层也处于压缩状态,而外层则处于拉伸状态。这个巧妙的技巧减少了当血管受压时内壁的应力集中,从而保护它免于失效。这种隐藏的张力可以通过著名的“开口角”实验来测量:将一个环形组织切开,它弹开的角度揭示了维持其形状的预应变的大小。

从各向异性骨骼到高效肺脏

我们讨论的原理并不仅限于一维的组织条带。我们的身体是三维的,其组成部分通常是各向异性的——在一个方向上比另一个方向更强。考虑一下你腿部的骨骼。它被优化以沿其长度方向承重。要描述这样的材料,我们必须推广我们的定律。QLV框架可以使用张量的语言优雅地扩展到三维。应力σij\sigma_{ij}σij​和应变εkl\varepsilon_{kl}εkl​变成张量,材料属性由四阶张量描述,这些张量捕捉了瞬时刚度和粘弹性松弛的方向依赖性。数学可能看起来更复杂,但分离非线性弹性响应和线性时间依赖性记忆的核心思想保持不变,这证明了该框架的统一力量。

力学与生物功能的这种相互作用在肺部得到了最惊人的体现之一。身体如何在一个如此紧凑的空间内最大化氧气和二氧化碳的交换?答案在于肺毛细血管的独特结构。肺泡壁中的毛细血管并非由微小的、独立的管子组成的网络,而是形成了一个几乎连续的血液“薄片”,仅由小的结缔组织柱子隔开。这种“片流”模型是Fung的另一个关键见解,是一种最优的工程解决方案。通过将血液铺展成薄膜,它最大化了可用于扩散的表面积(AAA),同时最小化了气体分子在空气和红细胞之间必须行进的距离(TTT)。根据Fick扩散定律,气体交换速率与A/TA/TA/T成正比。片流设计最大化了这一比率,使我们的肺成为极其高效的气体交换器。

工程师的工具箱:从实验台到数字世界

如果没有严格的实验验证,这些优美的理论将停留在推测的领域。如何测量一种既非线性又具粘弹性的材料的属性?一种巧妙的方法是对组织施加一个大的、恒定的“偏置”应变,然后在此之上叠加一个微小的正弦振荡。通过保持振荡微小,对其的响应实际上是线性的,这使得科学家能够在该特定的拉伸水平下探测材料的局部刚度和阻尼特性。通过在不同的偏置应变和频率下重复此操作,可以绘制出材料的完整行为,并检验像QLV这样的理论的预测。

一旦得到验证,这些模型就成为预测和设计的强大工具,尤其是在计算模拟领域。例如,在牙科和口腔生物力学中,了解牙齿在其牙槽中的运动方式对于正畸治疗至关重要。固定牙齿的牙周膜(PDL)是一种粘弹性组织,完全可以用QLV来描述。为了在计算机模拟(如有限元方法(FEM)模型)中使用这一理论,光滑的松弛函数G(t)G(t)G(t)通常用一个称为Prony级数的衰减指数和来近似。这种数学技巧使得复杂的历史依赖行为能够被高效计算,从而实现能够预测牙齿移动和优化正畸治疗的模拟。

这些思想的影响远远超出了生物领域。考虑一下制造更好电池的挑战。锂离子电池的一个主要失效模式是电极的机械退化。在充放电过程中,锂离子进出活性颗粒,导致它们反复膨胀和收缩。这给将电极粘合在一起的聚合物粘合剂带来了巨大的机械应力。这种粘合剂是一种粘弹性聚合物,可以用我们用于韧带的完全相同的QLV方程来描述!通过应用这个框架,工程师可以构建电池的“数字孪生”来模拟机械失效,并设计出更具韧性的粘合剂,从而实现更持久、更安全的储能。同样,连续介质力学的相同原理也可以应用于理解土壤和粘土等地质材料的时间依赖性行为,这对于土木工程和地质学至关重要。

一点谦逊之言:一个优美思想的局限

尽管QLV模型拥有强大的力量和优雅,我们必须以一种科学的谦逊来结束,这也是Feynman本人所倡导的一种品质。QLV模型,以其对时间和应变效应的清晰分离,终究只是一个模型。它是对现实的一个卓越近似,但不是现实本身。对于许多组织,如皮肤,如果你将它们拉伸到非常大的变形,可分离性假设便开始失效。实验数据表明,松弛曲线的形状会随着应变量的变化而改变。这表明,松弛的微观机制——例如流体的运动和聚合物链的重构——本身也受到应变水平的影响。看来,自然界总是比我们最优雅的理论要稍微更耦合、更复杂一些。

这并未削弱该理论的价值。恰恰相反,它指明了前进的道路。通过理解一个模型在何处成功、在何处失败,我们能更多地了解其 underlying physics,并被驱使去创造更精炼的描述。从一个关于肌腱拉伸性的简单观察到下一代电池的设计,这一历程证明了基本原理的力量。它揭示了世界力学行为的深层统一性,从生命体到工程造物,这是一曲由物理定律指挥的美丽交响乐。