第七章：动脉粥样硬化斑块应力分析综述
书名：Multi-Modality Atherosclerosis Imaging and Diagnosis
作者：L. Saba et al. (eds.)
章节作者：Hao Gao and Quan Long
出版年份：2013
章节标题：Atherosclerosis Plaque Stress Analysis: A Review
DOI：10.1007/978-1-4614-7425-8_7
第一节——章节概述
本章是动脉粥样硬化斑块应力分析的系统性综述，由Gao和Long撰写。斑块破裂是导致中风和心肌梗死的主要机制之一，但具体破裂机制至今尚未完全阐明。本章系统回顾了斑块应力分析从二维结构模型到三维流固耦合模型的发展历程，涵盖了理想化模型到患者特异性模型的几何精度提升，以及纯结构分析到流固耦合分析的技术演进。
斑块易损性的评估需要综合考虑斑块形态学和生物力学环境两个维度。形态学特征包括大脂质池、薄纤维帽、高炎症细胞含量等；而生物力学因素则主要包括斑块壁面拉伸应力（Wall Tensile Stress）和壁面剪切应力（Wall Shear Stress）。这两类因素共同决定了斑块的稳定性状态。
从方法论角度来看，本章详细介绍了有限元方法在斑块应力分析中的广泛应用。建立斑块力学应力分析模型通常涉及五个主要任务：获取斑块组分、重建斑块几何形状、施加力学载荷和动态边界条件、定义合适的材料本构模型，以及执行有限元模拟和后处理。不同复杂度的模拟研究在各个任务上具有不同的实现方式。
本章的阅读需要具备生物力学、有限元分析和医学影像学的基础知识。对于理解后续内容，读者应熟悉弹性力学的基本概念、血管壁材料的非线性本构关系，以及计算力学中有限元方法的原理。
第二节——关键问题与研究动机
核心科学问题
动脉粥样硬化斑块应力分析领域存在若干关键的科学问题，这些问题的解决对于理解斑块破裂机制至关重要：
问题一：斑块破裂的生物力学机制是什么？
斑块破裂被认为是薄纤维帽的机械失效过程。多种研究表明，斑块破裂往往发生在斑块肩部区域——这是由于纤维帽与下方脂质池之间刚度显著差异导致的高应力集中区域。然而，确切的破裂触发条件尚未完全确定。
问题二：如何准确预测斑块应力分布？
由于无法直接测量活体内斑块壁面应力，数值方法特别是有限元方法成为获取斑块区域应力的主要途径。然而，从医学影像获取的真实的几何形状、材料参数和边界条件都存在不同程度的不确定性，这些因素都会影响应力预测的准确性。
问题三：哪些因素决定了斑块的机械稳定性？
研究表明，脂质池大小、纤维帽厚度、钙化位置、斑块形状和轴向拉伸等因素都会对应力分布产生影响。理解这些因素的影响规律对于斑块易损性评估具有重要意义。
现有研究的不足
在Gao和Long的综述之前，该领域存在以下主要不足：二维模型假设平面应力或平面应变条件，简化了真实的生理状态；理想化模型与真实患者特异性几何形状存在显著差异；材料模型多为均匀各向同性假设，未能反映血管壁的各向异性和多层结构特性；流固耦合效应对斑块应力分布的影响未被充分考虑。
研究的广泛意义
斑块应力分析的深入研究不仅具有学术价值，更具有重要的临床应用前景。如果能够建立准确的斑块应力预测模型，就可以根据患者的特异性几何和材料参数评估其斑块破裂风险，为临床治疗决策提供科学依据。此外，理解生物力学因素在动脉粥样硬化发生、发展和破裂中的作用，有助于揭示疾病的病理生理机制，为新疗法的开发提供理论基础。
第三节——主要公式与推导
3.1 斑块破裂的力学假设
斑块破裂的生物力学研究基于以下主要假设：
假设一：局部最大应力假说
局部最大应力假说是最被广泛接受的斑块破裂机制。该假说认为，当纤维帽区域内的局部最大应力超过纤维帽组织所能承受的临界应力值时，斑块将发生破裂。Cheng等人基于二维应力分析与破裂和稳定动脉粥样硬化病变的组织病理学对照研究，提出了冠状动脉斑块300 kPa的临界应力值。这一临界值成为后续许多研究的参考标准。
假设二：疲劳失效假说
Bank假设由脉动血流引起的纤维帽机械疲劳可能是导致斑块破裂的重要因素。基于这一假设，研究者模拟了斑块疲劳裂纹的扩展过程，该过程基于斑块内应力分布的演化。
假设三：极高壁面剪切应力假说
最近的研究发现，极高的壁面剪切应力与斑块破裂相关。这一发现基于斑块溃疡发生在高壁面剪切应力位置的观察。
假设四：微钙化导致的脱粘效应
Vengrenyuk等人提出了基于薄纤维帽中细胞级微钙化引起应力集中的易损斑块破裂假说。微钙化在细胞尺度上可引起局部应力集中，导致界面脱粘。
3.2 有限元分析的基本方程
在斑块应力分析中，有限元方法求解的基本方程包括平衡方程、几何方程和本构方程。
平衡方程：
其中 \(\sigma\) 为应力张量，\(f\) 为体积力矢量。在斑块分析中，通常忽略体积力，只考虑表面载荷。
几何方程（应变-位移关系）：
其中 \(\varepsilon\) 为应变张量，\(u\) 为位移矢量。
本构方程（应力-应变关系）：
对于动脉壁组织，通常采用非线性超弹性本构模型。Mooney-Rivlin模型是常用的应变能密度函数形式：
其中 \(W\) 为应变能密度，\(C_1\) 和 \(C_2\) 为材料常数，\(I_1\) 和 \(I_2\) 为应变不变量，\(\kappa\) 为体积模量，\(J\) 为体积比。
对于各向异性材料，如包含纤维方向的血管壁组织，Holzapfel等人提出了可压缩的各向异性本构模型：
其中 \(k_1\) 和 \(k_2\) 为描述纤维刚度的材料参数，\(E_1\) 和 \(E_2\) 为Green应变分量。
3.3 关键应力指标
Von Mises应力：
在有限元后处理中，Von Mises应力是常用的等效应力指标，用于评估材料的屈服或破裂风险：
其中 \(\sigma_1\)、\(\sigma_2\)、\(\sigma_3\) 为主应力分量。
第一主应力：
在流固耦合分析中，第一主应力（最大拉伸应力）常被用来表征斑块壁面的拉伸应力状态：
该指标在Fig. 7.5c中被用于展示斑块壁面拉伸应力的分布特征。
3.4 临界应力准则
斑块易损性评估的临界应力准则可以表示为：
其中 \(\sigma_{max}\) 为纤维帽区域的最大应力，\(\sigma_{threshold}\) 为临界应力值（通常取300 kPa）。
3.5 应力敏感性分析
Williamson等人的研究显示，动脉壁应力对材料参数变化的敏感性较低，但不同材料模型之间的应力集中区域应力值可相差达30%。这表明在解释应力结果时需要谨慎选择合适的材料模型。
第四节——关键算法与建模方法
4.1 有限元方法概述
有限元方法（FEM）是求解斑块应力分析问题的主要数值方法。该方法的核心思想是将连续的计算域离散为有限个单元，通过单元内的插值函数近似场变量，并求解离散化的方程组得到数值解。
在斑块应力分析中，有限元模型的构建通常遵循以下流程：
步骤一：医学影像获取
利用MRI、IVUS、OCT或超声等医学影像技术获取斑块的空间信息和组分分布。Fig. 7.6展示了从活体MRI图像重建颈动脉斑块几何形状的一般流程，包括多序列MRI成像、斑块组分分割和三维几何重建。
步骤二：几何重建
通过图像分割提取斑块各组分的边界，然后进行三维表面重建和实体建模。Long等人的研究详细评估了从活体MRI图像重建三维颈动脉分叉几何形状的可重复性，发现普通颈动脉的几何形状在大多数情况下可以良好重建，而颈外动脉的可重复性最差。
步骤三：材料模型定义
根据斑块各组分的力学特性选择合适的本构模型。脂质池通常表现为较软的超弹性材料，钙化区域为高刚度材料，纤维帽则需要考虑其各向异性和非线性特性。
步骤四：边界条件和载荷施加
典型的边界条件包括：管腔内表面的血压载荷、动脉壁的轴向预拉伸、以及远端边界的位移约束。Fig. 7.2所示的研究中，假设动脉壁为线弹性材料（杨氏模量1 MPa，泊松比0.45），管腔压力为110 mmHg。
步骤五：有限元求解与后处理
使用商业软件（如ANSYS）或自编程序进行有限元求解，然后对应力、应变分布进行后处理分析。
4.2 流固耦合分析方法
流固耦合（FSI）分析是更真实地模拟斑块在生理环境下力学行为的先进方法。在生理情况下，血流与动脉壁结构之间存在相互作用：血流提供压力载荷和壁面拖曳力，动脉壁在这些载荷作用下发生变形，而变形后的几何又反过来影响血流场。
单向耦合方法：
先求解流体域得到压力和剪切力分布，然后将这些载荷施加到结构域求解固体应力。该方法假设流体场不受固体变形的影响。
双向耦合方法：
流体域和固体域同时求解，通过迭代的方式使两者的边界条件在界面处协调。该方法能够捕捉流体与结构之间的相互作用效应。
Tang等人发展的三维流固耦合模型从离体到在体的MR图像研究，揭示了大的脂质池和薄的纤维帽与极端应力和应变水平相关。脉动压力引起的大循环应力和应变可能导致斑块通过疲劳过程破裂。
4.3 建模中的关键参数选择
网格敏感性分析：
在有限元分析中，网格密度直接影响计算精度和成本。对于应力集中区域（如纤维帽肩部），需要使用较细的网格以准确捕捉应力梯度。
时间步长选择：
对于动态分析（如脉动血流条件下的应力分析），时间步长的选择需要足够小以捕捉载荷的快速变化。通常，心动周期被分为若干个时间步进行求解。
收敛准则：
迭代求解过程中，需要定义适当的收敛准则以确保解的精度。常用的收敛准则包括位移收敛和残余力收敛。
4.4 计算成本与效率
二维结构分析的计算成本相对较低，通常在几分钟内即可完成。三维流固耦合分析则需要更长的计算时间，特别是当考虑脉动载荷和非线性材料特性时，计算时间可达数小时甚至数天。然而，随着计算资源的不断提升和算法的优化，流固耦合分析的效率正在逐步提高。
第五节——主要结论
5.1 斑块应力分析的主要发现
综述系统总结了斑块应力分析领域数十年的研究成果，得出以下主要结论：
结论一：纤维帽厚度是影响斑块应力的关键因素
多项研究表明，纤维帽厚度的减小会显著增加峰值环向应力。Loree等人发现纤维帽厚度的减小会显著增加峰值环向应力，而狭窄程度的增加反而会降低峰值应力，这意味着破裂斑块可能与高度狭窄无关。Finet等人的研究进一步证实，无论斑块几何形状和成分如何，纤维帽厚度小于60 μm都会导致应力超过300 kPa的临界值。
结论二：斑块破裂发生在高应力位置
Cheng等人的研究表明，大多数斑块破裂部位发生在高应力区域附近，但峰值应力并不总是与破裂部位完全重合。最大环向应力在破裂斑块中显著高于稳定斑块，这为临界应力假说提供了实验依据。
结论三：钙化位置对斑块稳定性的影响不同
Huang等人的研究发现，钙化并不一定会降低斑块的稳定性。深部斑块中的钙化对稳定性的影响较小，甚至可能具有稳定作用；而纤维帽中的钙化则可能因刚度差异大而导致应力增加，促进斑块破裂。
结论四：二维模型倾向于高估应力水平
与三维模型相比，二维结构分析由于平面应力/应变假设和简化的材料模型，往往会高估斑块应力水平。Ohayon等人的三维有限元分析表明，峰值应力位置与斑块破裂位置的一致性在三维模型中优于二维模型。
结论五：流固耦合分析提供更真实的应力预测
三维流固耦合分析能够同时考虑血流动力学和血管壁结构动力学的耦合效应，提供比纯结构分析更准确的应力预测。Yang等人的研究比较了流固耦合模型和纯结构模型，发现两种模型在壁面剪切应力和最大主应力方面可能存在显著差异。
结论六：斑块应力与斑块进展呈负相关
Tang等人通过在体MRI基础上的二维/三维流固耦合模型证明了人类颈动脉粥样硬化斑块进展与斑块壁应力之间呈负相关。他们还提出，较低的壁面应力和较低的壁面剪切应力都可能促进斑块的持续进展。
结论七：斑块应力可用于易损性评估
基于应力信息的易损性指数与组织病理学评估的一致性达到85%，表明斑块应力分析在临床易损性评估中具有重要应用价值。Tang等人提出的局部最大应力假说和基于应力的计算斑块易损性指数为斑块风险评估提供了新思路。
5.2 斑块应力分析的发展脉络
从历史发展的角度看，斑块应力分析经历了以下演进过程：
从二维到三维：早期研究主要采用二维模型，虽然计算简便但存在过度简化的问题；随着计算技术和医学影像的发展，三维模型逐渐成为主流。
从理想化到患者特异性：理想化模型便于参数化研究，但与真实几何存在差距；现代医学影像技术使得获取患者特异性几何成为可能。
从纯结构到流固耦合：纯结构分析仅考虑载荷作用下的准静态响应；流固耦合分析则同时考虑血流与血管壁的相互作用，更真实地模拟生理环境。
第六节——挑战与开放问题
6.1 当前理论的主要 gaps
尽管斑块应力分析取得了显著进展，但以下理论问题仍未得到很好解决：
问题一：临界应力值的确定
300 kPa的临界应力值是基于特定研究提出的，但该值是否具有普适性仍存疑问。不同患者的纤维帽组织特性可能存在显著差异，统一使用一个临界值可能导致误判。
问题二：多因素耦合机制
斑块破裂不仅与力学因素相关，还涉及炎症活动、胶原降解、细胞凋亡等多种生物学过程。单纯依靠力学分析难以全面评估斑块的真实易损性。
问题三：残余应力的影响
Ohayon的研究表明，斑块内残余应力/应变可能对薄纤维帽的生理峰值应力产生显著影响；而Gao等人的研究则显示环向残余应力的影响有限。这些矛盾的结果表明残余应力对斑块稳定性的影响机制尚未被完全理解。
6.2 实验验证的困难
挑战一：无法直接测量活体应力
目前没有技术能够直接测量活体内斑块壁面应力，这使得模型验证变得极为困难。研究者只能通过间接方式（如破裂位置与应力分布的对应关系）来验证模型的合理性。
挑战二：破裂前后几何不可得
对于特定斑块，破裂前和破裂后的几何形状实际上是无法获取的，这限制了对模型预测能力的直接验证。
挑战三：材料参数的个体差异
斑块各组分的材料参数存在显著的个体差异和位置差异，获取准确的患者特异性材料参数仍是挑战。
6.3 建模技术的局限性
局限性一：几何重建的不确定性
医学影像的分辨率限制，特别是对于小于65 μm的薄纤维帽，往往只能由一到两个像素表示，这导致纤维帽重建存在较大的不确定性。Gao等人的研究表明，脂质区域的高估会导致更高的应力集中和整体更高的应力水平。
局限性二：材料模型的简化
虽然动脉壁组织表现出非线性、各向异性和粘弹性特性，但在大多数有限元模型中，这些特性被不同程度地简化。材料参数来源多为离体实验数据，可能与在体条件存在差异。
局限性三：边界条件的理想化
真实的生理边界条件非常复杂，包括脉动血流、周围组织的约束作用、血管的轴向预拉伸等。在模型中完全准确地再现这些条件仍然困难。
6.4 尚未解决的科学问题
问题一：斑块破裂的精确触发条件
尽管局部最大应力假说被广泛接受，但斑块破裂的确切触发条件仍不明确。是应力超过临界值导致即时破裂，还是疲劳累积导致渐进性损伤？
问题二：血流动力学与应力的协同作用
壁面剪切应力与壁面应力如何协同作用影响斑块破裂？两种力学信号之间是否存在耦合机制？
问题三：斑块进展的力学调控
斑块应力如何反过来调控斑块的生物学进展？力学信号与炎症反应之间的相互关系尚未被完全阐明。
第七节——个人反思与批判性分析
7.1 作者建模哲学的评价
Gao和Long的综述展现了一种严谨而务实的建模哲学。作者并没有盲目追求复杂的模型，而是系统梳理了从简单到复杂、从二维到三维、从理想化到患者特异性的发展脉络。这种渐进式的建模策略有助于理解每个模型假设的合理性及其局限性。
值得注意的是，作者在综述中保持了客观中立的态度，既肯定了局部最大应力假说的广泛应用，也提及了其他假说（如疲劳失效、极高剪切应力、微钙化脱粘等），并指出这些假说各自的支持证据和争议点。
7.2 数学简化的得失
斑块应力分析中的数学简化带来了一些重要的洞察，也不可避免地引入了一定的偏差：
简化带来的便利：
平面应力/应变假设大大简化了二维分析，使得参数化研究成为可能。均匀材料假设和理想化几何模型使得研究者能够系统地研究单个因素（如纤维帽厚度、狭窄程度）的影响。这些简化为理解斑块力学的基本规律提供了重要基础。
简化带来的偏差：
二维模型倾向于高估应力水平，这意味着基于二维分析的安全结论可能需要重新审视。各向同性材料假设忽略了血管壁中胶原纤维的各向异性特性。理想化几何与真实几何的差异可能导致应力分布模式的偏差。
未来方向：
随着计算能力的提升和医学影像技术的发展，更精细的模型将变得可行。三维各向异性非线性材料模型、流固耦合分析、患者特异性材料参数的反演识别等都将是未来研究的重要方向。
7.3 对本研究的启发
本研究对我的启发主要体现在以下几个方面：
启发一：多学科交叉的重要性
斑块应力分析涉及生物力学、医学影像、计算科学和临床医学等多个学科的深度融合。这种交叉研究模式对于解决复杂的医学问题具有重要示范意义。
启发二：模型验证的困难
对于许多临床问题，由于伦理和技术的限制，获取直接验证数据几乎不可能。研究者需要发展间接验证方法，如通过临床预后数据来验证模型预测能力。
启发三：临床转化的潜力
斑块应力分析已经从纯学术研究走向临床应用的可能性，如基于应力的斑块易损性评估。这提示基础研究者应始终关注研究成果的临床转化价值。
7.4 值得进一步探索的问题
问题一：能否建立统一的多尺度模型？
从细胞级别到组织级别再到器官级别，存在多个尺度的力学行为。建立一个能够跨尺度统一描述斑块力学特性的模型将是一个重要突破。
问题二：机器学习能否加速模型参数反演？
随着深度学习技术的发展，能否利用神经网络从医学影像直接预测斑块应力分布？这将大大简化计算的复杂性。
问题三：能否实现实时应力评估？
如果能够实现基于常规临床影像的快速应力评估，将为急性事件的风险预测提供重要工具。
7.5 对研究者的建议
对于希望进入这一领域的研究者，我有以下建议：
首先，应扎实掌握弹性力学、有限元方法和生物力学的基础知识，这是开展斑块应力分析的必备前提。
其次，应充分了解各种医学影像技术的特点和局限性，因为几何重建的准确性直接决定了应力分析结果的可靠性。
最后，应注重与临床医生的合作，确保研究问题具有临床相关性，研究成果能够得到临床验证。
公式汇总
| 编号 | 名称 | 形式 | 物理意义 | 类型 |
|---|---|---|---|---|
| (7.1) | 平衡方程 | \(\nabla \cdot \sigma + f = 0\) | 描述弹性体内部应力的平衡条件 | (T) |
| (7.2) | 几何方程 | \(\varepsilon = \frac{1}{2}(\nabla u + (\nabla u)^T)\) | 应变与位移之间的线性关系 | (T) |
| (7.3) | Mooney-Rivlin应变能函数 | \(W = C_1(I_1 - 3) + C_2(I_2 - 3) + \frac{\kappa}{2}(J - 1)^2\) | 超弹性材料的本构模型 | (T) |
| (7.4) | 各向异性本构模型 | \(W = \frac{\mu}{2}(I_1 - 3) + \frac{k_1}{2k_2}(e^{k_2(E^2_1 + E^2_2)} - 1) + \frac{\kappa}{2}(J - 1)^2 - \mu \ln J\) | 考虑纤维方向的动脉壁材料模型 | (T) |
| (7.5) | Von Mises应力 | \(\sigma_{vm} = \sqrt{\frac{(\sigma_1 - \sigma_2)^2 + (\sigma_2 - \sigma_3)^2 + (\sigma_3 - \sigma_1)^2}{2}}\) | 材料等效屈服判据 | (T) |
| (7.6) | 第一主应力 | \(\sigma_1 = \max(\sigma_1, \sigma_2, \sigma_3)\) | 最大拉伸应力指标 | (T) |
| (7.7) | 临界应力准则 | \(\sigma_{critical} = \begin{cases} \text{稳定}, & \sigma_{max} < \sigma_{threshold} \\ \text{易损}, & \sigma_{max} \geq \sigma_{threshold} \end{cases}\) | 斑块易损性判断准则 | (E) |
注：(T)=理论推导，(E)=经验公式
附录：书中元数据
- 书名：Multi-Modality Atherosclerosis Imaging and Diagnosis
- 作者：L. Saba et al. (eds.)
- 出版社：Springer Science+Business Media
- 出版年份：2014（原书2013年出版）
- DOI：10.1007/978-1-4614-7425-8
- 核心主题：动脉粥样硬化的多模态成像与诊断
- 目标读者：生物医学工程研究人员、心血管疾病研究者、放射科医生
- 前置知识：生物力学、有限元分析、医学影像学基础
- 相关书籍：《Biomechanics of Cardiovascular System》、 《Computational Biomechanics for Medicine》