玻尔百科有生命的生物组织是材料工程的杰作，具有非凡的强度、柔韧性和回弹性。然而，用数学精确描述其力学行为是一项重大挑战。与遵循胡克线性定律的简单钢制弹簧不同，皮肤、肌肉和动脉等组织表现出复杂的非线性响应：它们起初柔软且易于变形，但随着拉伸会变得异常坚硬。这种 J 形应力-应变曲线是线性模型无法捕捉的一个基本特征，从而在简单力学与复杂生物学之间造成了知识鸿沟。
本文探讨了生物工程师 Y.C. Fung 的开创性工作，他开发了一个强大的数学框架来描述这种行为。基于指数弹性定律的冯氏模型已成为现代生物力学的基石。通过阅读本文，您将对这一基本理论及其深远影响有深入的理解。首先，在“原理与机制”一章中，我们将剖析该模型的数学核心，探索其指数形式、参数的物理意义、为解释组织的纤维状各向异性特性所做的扩展，以及粘弹性效应的引入。在这一理论基础之后，“应用与跨学科联系”一章将展示该模型如何付诸实践，从解读实验数据到驱动正在彻底改变医学和生物工程的复杂计算机模拟。
引言之后，我们已达成共识：生命组织是一种非凡的材料，我们希望揭开其力学奥秘。但如何从哲学层面的欣赏走向真正的、可预测的科学呢？我们需要一个数学描述，一条“定律”，来告诉我们一块组织在受推或拉时如何响应。
事情从这里开始变得有趣。一个简单的钢制弹簧遵循胡克定律：力与伸长成正比。伸长加倍，力也加倍。这是一种优美简洁的线性关系。但如果你拉伸一块皮肤或肌肉，你会发现它展现出截然不同的行为。起初，它相当柔顺，很容易被拉伸。但随着你进一步拉伸，它会迅速变得越来越硬。它会猛烈地反抗。这种非线性的硬化行为是大多数软生物组织的标志。我们的首要任务就是捕捉这一基本特性。
伟大的生物工程师 Y.C. Fung 提出了一个优美简洁而又强大的思想。他提出，组织的刚度并非随应变线性增长，而是指数增长。这一思想的数学核心是应变能密度函数，我们称之为 。可以将 理解为使一立方厘米材料变形时所储存的能量。材料弹性响应的所有信息都蕴含在这一个函数中。
对于一种在所有方向上性质都相同的简单材料（各向同性材料），冯氏模型的应变能密度函数如下所示:
不要被这些符号吓到，我们来分解一下。 是一个衡量材料所经历的等容（保持体积）应变总量的数值。对于未变形的材料，，所以指数项为零，，储存的能量 也为零。这只是一个比较花哨的说法，即松弛的组织没有储存能量。
现在，参数 和 又代表什么呢？
例如，如果我们想象一个简单的剪切变形，比如将一块明胶的顶部侧向滑动，产生的剪应力为 ，其中 是剪切量。你可以在公式中直接看到：应力不仅仅与应变 成正比，它还乘以一个增长非常非常快的指数项。这种指数硬化现象正是我们在实验室中观察到的。
简单的各向同性模型是一个很好的起点，但我们可以将其推广以描述更复杂的材料。更通用的冯氏模型如下所示：
此处， 是应变张量 的一个二次函数。可以把它看作是一种更复杂的测量“总应变”的方法，它考虑了所有方向上的拉伸和剪切。它就像一个篮子，我们把所有不同的应变分量都收集在里面，每个分量可能有不同的权重。
现在，我们必须对函数 施加一个关键条件：它必须是正定的。这听起来像一个抽象的数学要求，但它根植于深刻的物理原理。这意味着对于任何可能的应变 ， 都必须为正（除非 为零，此时 ）。为什么这如此重要？
一个稳定的基态：因为对于任何非零应变，，所以应变能 也大于零。唯一能量为零的状态是未变形状态（）。这保证了材料有一个唯一的、稳定的、无应力的参考状态。材料“倾向于”处于这种状态，你必须对它做功才能使其偏离。否则，材料可能会自发变形并释放能量，那将是一个非常奇怪的世界！
一个行为良好的开端：在非常小的应变下，任何平滑曲线看起来都像一条直线。 的正定性确保了我们这个花哨的非线性模型在微小变形时能够平滑地简化为我们所熟悉的、稳定的线性弹性理论。当你只是轻轻戳一下材料时，它不会做出任何奇怪的反应。
不会迷失：这个条件也确保了能量函数是凸的。从几何上看，这意味着能量对应变的图像看起来像一个碗。这是数学家的梦想，因为它保证了对于给定的载荷，存在一个唯一的、稳定的变形状态。材料的响应是可预测且行为良好的。
无限的抵抗力：最后，它意味着当应变趋于无穷大时，所需的能量也趋于无穷大。你不能无限地拉伸材料；它会变得无限坚硬。这种性质，被称为矫顽性，是固体的一个基本物理要求。
因此，这一个“简单”的数学条件——正定性——是确保我们的模型不仅仅是一堆方程，而是对真实材料的物理和热力学上都合理的描述的关键。
到目前为止，我们讨论的都是各向同性材料，它们在所有方向上都具有相同的性质。但很少有生物组织是这样的。肌腱、动脉壁、肌肉——它们都因为内嵌的胶原蛋白或肌纤维而具有优选方向。它们是各向异性的。
冯氏模型的优雅之处在此得以体现。我们可以通过在 函数中简单地添加项来包含这种方向性。对于像肌腱这样具有单族纤维的组织，我们可以写出:
这里的诠释非常美妙：
这非常直观。如果你沿着纤维方向拉伸组织，你会激活强大的 项，材料会非常坚硬。如果你在垂直于纤维的平面上剪切组织，纤维长度不变（）， 项为零，响应仅由较软的基质通过 项控制。该模型自然地捕捉了这种各向异性行为。
为什么要止步于一族纤维呢？动脉壁可能有两族纤维螺旋缠绕。没问题。我们只需为第二族纤维再添加一项:
这个框架是一个模块化且强大的“语言”，用于描述生命组织的复杂结构。
还有两个生物学的现实需要进入我们的画面：水和黏性物质。
首先，软组织主要由水组成。在所有实际应用中，它们都是不可压缩的。你可以轻易改变它们的形状，但无法将它们压缩到更小的体积。用力学术语来说，这意味着它们的泊松比非常接近 0.5。我们的模型可以处理这一点。我们可以将应变能分为控制形状变化的部分（等容部分）和惩罚体积变化的部分（体积部分）。通过设定非常高的体积变化惩罚，我们可以模拟这种近乎不可压缩性。仔细分析表明，我们的模型预测的有效泊松比为 ，其中 是材料的剪切刚度， 是其体积模量。随着材料变得越来越难以压缩（ 变得非常大），泊松比越来越接近理想值 0.5。能从这个复杂的非线性理论中看到大学物理中的熟悉概念，真是太美妙了。
其次，到目前为止我们所讨论的都是纯弹性——就像一个完美的弹簧，能把你施加的所有能量都还给你。但真实的组织更像是一个弹簧与一个黏性缓冲器的组合。它们是粘弹性的。如果你拉伸一块组织然后释放它，它不会沿原路返回。它会损失一些能量，主要以热量形式散失。一个周期内的这种能量损失称为滞后。像我们目前讨论的纯弹性模型无法捕捉这一点，因为对它们来说，能量变化只取决于起点和终点，而与路径无关。从热力学上讲，它们的耗散总是为零。
为了解释这一点，我们需要为材料添加“记忆”。准线性粘弹性（QLV）理论，同样由 Fung 开创，通过一个遗传积分来实现这一点。其思想是，你现在感受到的应力不仅取决于现在的应变，还取决于其整个应变历史，并且存在一种“衰退记忆”，即最近的应变影响更大。这一扩展使得模型能够预测应力松弛（如果你将组织保持在恒定伸长状态，所需的力会慢慢减小）和蠕变（如果你施加恒定的力，组织会慢慢继续伸长）等现象，这些都是真实生物组织的决定性特征。
我们现在有了一个复杂而强大的模型。但模型只是一个故事。要使其成为科学，我们需要通过实验将其与现实世界联系起来。这就引出了关于理论与实验相互作用的最后一个深刻观点。
想象一下，我们有一个横向各向同性模型，其基质参数为 ，纤维参数为 。你如何确定它们对于特定组织的值？你可能会想：“简单，我做两个实验。我沿着纤维拉伸它得到一条应力-拉伸曲线，然后再横向拉伸它得到另一条。”
事实证明，这还不够！对于这种类型的指数模型，这两种简单的拉伸测试都会将 和 的效应混合在一起，形成一种“共谋”。 和 的不同组合在这些测试中可以产生几乎相同的应力-拉伸曲线，使得我们无法将它们区分开来。这些参数不是唯一可识别的。
关键在于进行一种能从根本上以不同方式探测材料的测试。在垂直于纤维的平面上进行简单的剪切测试是理想的选择。正如我们所见，这种变形根本不会拉伸纤维，因此产生的剪应力只取决于基质参数 。这个测试分离出了基质的响应，让我们能够确定 。一旦知道了 ，我们就可以回到简单的拉伸数据中，轻松地计算出纤维参数 。
这是科学方法中一个优美的教训。拥有一个好的理论模型只是成功了一半。你还必须是一个聪明的实验者，设计能够分离不同物理效应、打破参数之间“共谋”的测试。对自然的真正理解在于优美的理论与敏锐的实验之间优雅的共舞。
现在我们已经探索了冯氏模型的数学核心，你可能会问一个很合理的问题：“这一切都很优雅，但它到底有什么用处？” 物理定律或数学模型就像一个工具。一箱闪亮而未被使用的工具是令人惋惜的；只有当它们被用于工作，塑造我们的世界和我们对世界的理解时，它们真正的美才会显现出来。在本章中，我们将把冯氏模型从工具箱中取出，看看它如何帮助我们理解、预测甚至改造生命组织的世界。这段旅程将带我们从实验台走向超级计算机，从材料的基本力学走向我们自身的呼吸节奏。
一个模型，无论多么复杂，在被现实赋予信息之前，都只是一个空壳。冯氏模型的第一个也是最关键的应用是它在与一块真实生物组织的对话中发挥作用。想象一下，我们取一块动脉壁或皮肤的小样本，并将其放置在一个可以拉伸它的机器中。当我们拉伸样本时，我们测量它所产生的抵抗力。我们会看到什么？我们看到的不是像普通弹簧那样的简单直线。相反，我们看到一条曲线，起初很平缓，然后以惊人的陡峭度上升。这就是著名的“J形”应力-应变曲线，是软生物组织的一个力学特征。
为什么是这种形状？在小拉伸下，组织是松软的；其缠结的蛋白质纤维网络（主要是胶原蛋白）才刚刚开始展开。但随着拉伸的增加，这些纤维被拉紧，组织突然变得异常坚硬和有抵抗力，就像一根绷紧的绳子。这可以防止组织过度拉伸而失效。
冯氏模型的美妙之处在此初见端倪。更简单的模型，如新胡克模型，预测在小变形下应力随应变线性增长，无法捕捉这种急剧的硬化现象。其他模型，如 Ogden 模型，使用多项式项，可以拟合曲线，但通常需要许多参数。而冯氏模型，凭借其指数形式 ，则非常适合这项任务。指数函数自然地产生这种 J 形曲线，仅用两个参数 和 就捕捉到了温和的初始响应和随后的爆炸性硬化。
但我们如何找到 和 的正确值呢？我们去问材料本身！我们进行实验，记录应力对拉伸的数据点，然后使用数值方法——一种复杂的、自动化的“猜测-检验”——来找到使模型曲线尽可能贴近我们实验数据的 和 值。这个过程，称为参数拟合或标定，赋予了组织话语权，让它告诉模型自己的行为方式。
当然，身体中的组织很少只被单向拉伸。动脉像一个圆柱形气球一样膨胀，沿其周向和轴向都会拉伸。为了建立一个真正稳健的模型，我们必须倾听组织对更复杂变形的响应。通过进行多种测试，例如单轴拉伸和等双轴膨胀（在两个方向上等量拉伸），并将模型同时拟合到所有数据，我们就能创建一个远比之前可靠的材料“数字孪生”。我们甚至可以建立精确的数学关系，来比较像冯氏模型这样的模型与 Ogden 模型的表现，不仅在拉伸开始时匹配它们的行为，而且在更大、更具生理相关性的拉伸下也能做到。
作为一阶近似，人们可能认为一块组织是一个均匀的、各向同性的“团块”——一种在所有方向上性质都相同的材料。我们目前讨论的基本冯氏模型就做了这个假设。但仔细观察一块牛排，你会看到肌肉的纹理。观察动脉壁，你会发现坚韧的胶原纤维像桶箍一样缠绕着它。大多数生物组织都是各向异性的；它们是纤维增强复合材料，被精巧地设计成在承受最高载荷的特定方向上更为坚固。
我们的模型就此失效了吗？完全没有！支撑冯氏模型的连续介质力学框架的真正力量在于其灵活性。我们不必抛弃这个模型，只需丰富它。我们可以在应变能函数中引入新的项，这些项仅在材料沿特定纤维族方向拉伸时才被激活。对于一个具有两族纤维（方向角为 和 ）的组织，能量函数可能如下所示：
在此， 和 这两项是特殊的不变量，它们衡量了两个纤维族各自内部拉伸的平方。这是一个非常直观的想法：总能量是储存在一般基质中的能量与拉伸单个纤维群所储存能量的总和。这种方法使我们能够以更高的保真度对从动脉、心脏瓣膜到皮肤和韧带的各种组织进行建模。
该框架也适用于不同的几何形状。许多组织，如心脏周围的心包囊或植物的叶子，都非常薄。将它们视为完整的三维实体在计算上可能非常浪费。相反，我们可以运用相同的原理发展出二维的膜理论。通过定义特殊的二维不变量，我们可以创建一个专门针对这些薄结构的冯氏类型模型，捕捉它们在平面应力状态下的面内拉伸和纤维增强效应。这是物理学和工程学中一种常见而强大的策略：识别问题的基本特征，并相应地简化数学。
所以，我们有了一个经过实验数据精细调校、能够解释组织结构的数学模型。我们能用它做什么呢？我们可以构建一个虚拟世界。这些本构模型最深刻的应用之一是在有限元分析（FEA）中。
FEA 背后的原理很简单：要分析一个复杂物体（如人体心脏或植入了支架的动脉）的力学行为，你可以将其在数字上分解为大量微小的、简单的形状（即“有限元”），这些形状在节点处相连。然后计算机可以为这个相互连接的网格求解力平衡的基本方程。冯氏模型在其中扮演着每一个微小单元的“本构律”或“材料属性”的角色。
为了让这些庞大的计算得以实现，计算机不仅需要知道给定应变下一个单元内的应力，还需要知道应力如何随应变的变化而变化。这个量，一个被称为一致切线模量的四阶张量，是数值解快速、稳定收敛的关键。推导这个切线模量是张量微积分中一项具有挑战性但至关重要的练习。当最终的表达式被转换成适合计算机的矩阵形式时，就能够高效地模拟极其复杂的生物事件。正是这套机制让生物医学工程师能够设计出更好的医疗设备，外科医生能够规划手术方案，研究人员能够理解疾病。当你看到一个关于跳动心脏的惊人计算机模拟时，你所见证的正是冯氏模型——或其后继者之一——在发挥作用。
谜题还有最后一块，而它可能是最重要的一块。生命组织不仅仅像橡皮筋一样具有弹性；它们是*粘弹性的*。如果你突然施加一个拉伸并保持住，应力不会保持不变——它会逐渐松弛。如果你施加一个恒定的载荷，组织会慢慢“蠕变”。这是因为组织充满了水和长而缠结的聚合物链。抵抗变形不仅涉及拉伸原子键，还涉及将流体推过一个多孔基质以及解开分子“意大利面”。这些过程需要时间并耗散能量。
我们如何捕捉这种时间依赖性？一个非常简单的起点是用弹簧（代表弹性）和缓冲器（代表粘性）的组合来为组织建模。例如，一个简单的弹簧-缓冲器组合可以出色地预测一个真实的生理现象：你的肺的顺应性（充气的难易程度）取决于你的呼吸速度。在慢速呼吸时，粘性元件有时间松弛，肺显得更具顺应性。在快速呼吸时，缓冲器没有时间移动，感觉“僵硬”，肺的顺应性就显得较低。
这个简单的线性模型很有洞察力，但它无法处理真实组织的大应变和非线性弹性。这正是 Y.C. Fung 做出他最著名贡献的地方：准线性粘弹性（QLV）理论。这个名字听起来令人生畏，但其核心思想却惊人地优雅和简单。
Fung 认识到组织的应力响应对应变不是线性的，但他假设它可能对其自身的弹性响应是线性的。线性粘弹性的玻尔兹曼叠加原理指出，总应力是过去所有应变增量响应的总和（或积分）。QLV 提出了类似但不同的观点：总应力是过去所有瞬时弹性应力增量响应的总和。
数学上，这意味着应力由一个积分给出，如下所示：
此处， 是我们一直在讨论的非线性弹性应力（可以来自我们的指数冯氏模型），而 是一个“简化松弛函数”，描述该弹性响应如何随时间衰减。该模型巧妙地将对应变的非线性依赖性与时间依赖的松弛分离开来。它之所以是“准线性”的，是因为叠加原理仍然有效，但作用对象不是应变 ，而是转换后的变量 。
QLV 理论是一场革命。它提供了一个热力学一致的框架，能够准确地描述在几乎所有软组织中观察到的应力松弛、蠕变和滞后现象。虽然其可分离性的基本假设有其局限性——它在简单加载类型下效果最好，在复杂的旋转拉伸下可能会失效——但它的强大功能和实用性使其几十年来一直是生物力学中不可或缺的工具。
从拟合实验数据到模拟虚拟手术，再到为生命的节律建模，冯氏模型及其 QLV 扩展形成了一个连贯而强大的知识框架。它证明了这样一个理念：即使是复杂而“混乱”的生物学世界，也受制于深刻的数学之美与统一性的原则。