第7章：动脉粥样硬化斑块应力分析综述（Atherosclerosis Plaque Stress Analysis: A Review）

7.1 引言（Introduction）

卒中与心肌梗死是全球主要致死原因之一，其病理机制在很大程度上与动脉粥样硬化斑块破裂后继发血栓形成有关。在颈动脉或冠状动脉中，斑块可以进展为高度狭窄从而限制下游血供并引起动脉壁重塑；但更常见的情形是，那些仅有低度狭窄的斑块发生破裂并继发血栓形成，继而栓塞血流，引发卒中或心梗。因此，斑块破裂的精确机制在很长一段时间内并不清楚，临床上又往往在破裂发生前毫无症状。为此，研究者投入了大量工作寻找能够预测斑块易损性的特征。迄今为止，一些形态学特征——大脂质核、薄纤维帽、炎症细胞含量高——已被认为与斑块不稳定性相关。图 7.1 展示了一例易损颈动脉斑块的典型组织学切片：薄纤维帽覆盖于脂质区之上，薄纤维帽一旦破裂，脂质核心将进入血流并堵塞下游动脉。

作者指出，斑块形态学与斑块所处的生物力学环境共同决定其易损性。斑块壁面拉应力（wall tensile stress）与壁面剪切应力（wall shear stress）等生物力学因素被认为是斑块破裂过程中的重要因子，应纳入斑块易损性评估之中。本章聚焦斑块壁面应力分析的综述，覆盖从二维到三维、从理想化模型到患者特异性模型、从纯结构分析到流-固耦合（FSI）分析的发展脉络，并提供作者课题组的计算实例；但本章不回顾斑块的血流动力学研究。

由于目前缺乏能在体内直接测量斑块壁应力的技术，也几乎不可能获得同一斑块在破裂前后的几何形态，斑块破裂相关的生物力学因子迄今无法在体内得到充分验证。但在缺乏更优手段的情况下，数值方法（尤其是有限元方法 FEM）已被广泛用于获取斑块区域的应力。典型做法是：基于多组分（含脂质区、纤维帽等）的真实斑块几何构建 FEM 模型，并施加生理性力学载荷，预测斑块内的力学应力分布。

作者概括了构建斑块力学应力分析模型时通常涉及的五个主要任务：(a) 获取斑块各组分；(b) 重建斑块几何；(c) 施力学载荷与动态边界条件；(d) 定义合适的材料本构模型；(e) 执行 FEM 仿真并进行后处理。表 7.1 对不同复杂度仿真在上述各任务上的具体内容作了汇总，并按"二维结构应力分析—三维结构应力分析—FSI 分析"的层次罗列代表性研究。后续章节将按 2D 结构、3D 结构、3D 流体、2D/3D FSI 的顺序综述仿真技术的发展。

7.2 斑块易损性的生物力学研究（Biomechanical Study of Plaque Vulnerability）

这一组研究背后的核心假设是：斑块破裂是薄纤维帽的力学失效。斑块破裂的倾向部分源于纤维帽内力学应力的升高。多项研究已揭示斑块破裂常发生在肩部（shoulder）区域——该处是纤维帽与下方脂质池因刚度差异显著而形成的高应力集中区。生物力学领域提出的斑块破裂主要假说如下。

第一，局部最大应力假说：纤维帽内局部最大应力若超过纤维帽组织能承受的临界应力值，通常被认为是斑块破裂在生物力学层面上的主要原因。Cheng 等基于对破裂斑块与稳定斑块的二维应力分析，对冠状动脉斑块提出了 300 kPa 这一临界应力值。

第二，疲劳失效假说：Bank 假说认为，由脉动血流引起的纤维帽力学疲劳可能是斑块破裂的重要因素；后续有研究基于斑块内演化的应力分布对疲劳裂纹生长进行了模拟。

第三，极高壁面剪切应力假说：近期研究发现斑块溃疡常出现在壁面剪切应力高值区，因此极高壁面剪切应力被认为与斑块破裂相关。

第四，钙化微粒引起的脱粘效应：基于薄纤维帽中细胞层级微钙化可造成局部应力集中并导致界面脱粘这一事实，有学者提出应力诱导脱粘可致斑块破裂的假说；该效应也可被归入"局部最大应力"假说（因为微钙化会产生极高的局部应力）。

第五，其他假说：包括狭窄处湍流所致损伤、滋养血管（vasa vasorum）破裂等。

在上述假说中，局部最大应力假说接受度最高。然而，鉴于缺乏对斑块应力的体内直接测量和破裂前后几何信息，针对任何与应力相关的假说，获取特定斑块上精确的应力分布是必要前提；而在缺乏更优测量手段的情况下，数值方法（尤其是 FEM）已成为获取斑块区域应力的主要途径。FEM 模型已经从二维到三维、从理想化模型到患者特异性几何、从纯结构到 FSI 逐步发展。表 7.1 按"二维结构应力分析—三维结构应力分析—FSI 分析"层次罗列了各代表性研究及结果。后续章节按 2D 结构、3D 结构、3D 流体、2D/3D FSI 顺序综述仿真技术发展。

7.3 二维结构应力分析（Two-Dimensional Structure-Only Stress Analysis）

自 1990 年代以来，二维结构应力分析被广泛用于从理想化模型到患者特异性模型的斑块应力预测。理想化模型的优势是便于改变几何参数以考察各斑块特征对应力分布的独立影响。Loree 等基于理想化模型考察了狭窄程度与纤维帽厚度对应力分布的影响，发现纤维帽厚度减小会显著升高峰值周向应力，而狭窄程度增加反而使峰值应力下降——提示破裂斑块可能并不与高度狭窄相关。Cheng 等对同一斑块研究了破裂斑块中峰值应力位置与破裂部位的关系，结果显示大多数斑块破裂部位紧邻高应力区，但峰值应力位置并不总与破裂部位重合；他们还发现破裂斑块中最大周向应力显著高于稳定斑块的最大应力，并据此提出用 300 kPa 作为预测斑块破裂的临界应力值。Finet 等以二维理想化模型考察纤维帽厚度对斑块应力的影响，发现不论斑块几何与组分如何，纤维帽厚度小于 60 μm 即导致应力超过 300 kPa。

在材料模型方面，Beattie 等以人动脉粥样硬化主动脉的非均质材料模型进行应力分析，发现应力与应变能分布强烈依赖于斑块结构与组分，意味着多组分的合理材料模型对斑块应力分析极为重要。Williamson 等考察了病变动脉壁应力对材料参数变化的敏感性，发现壁应力对弹性模量变化不敏感，且在各项同性非线性、含残余应变的各项同性非线性、横观各向同性线性三种模型间相差不大；但应力集中区的应力在采用不同材料模型时可相差 30%。Huang 等在人体冠状动脉粥样硬化样本上考察钙化对斑块稳定性的影响，认为钙化并不降低斑块稳定性；然而 Veress 等认为钙化的位置对斑块稳定性有重要作用：深部钙化对稳定性影响小甚至稳定斑块，而纤维帽内钙化因刚度差异大可增加纤维帽（尤其钙化-纤维帽界面）应力、促进斑块破裂。

图 7.2 展示了作者课题组对一例颈动脉斑块组织学切片进行二维应力分析的结果：动脉壁被假设为线性各向同性材料（杨氏模量 1 MPa，泊松比 0.45），腔内施加 110 mmHg 血压；得到的 Von-Mises 应力分布显示纤维帽应力显著高于其他区域，最大应力出现在斑块肩部，达 253 kPa；脂质核附近局部应力极低。这一"纤维帽高应力、脂质核低应力"的特征被普遍认为与斑块破裂相关。

应力分析也已与组织学发现结合，以研究纤维帽内炎症活动与极端应力之间的相互作用。Lee 等考察了人冠状动脉粥样硬化样本中高周向应力区与 MMP 表达的相关性，认为高应力区细胞外胶原基质的降解会进一步削弱纤维帽、促进破裂。Howarth 等使用 USPIO 增强 MRI 将一例症状性重度颈动脉狭窄患者的巨噬细胞位置与斑块应力分布关联，结果支持巨噬细胞位置与斑块最大预测应力相关。Hallow 等用二维非均质有限元模型与所选炎症标志物的分布进行匹配，发现应力与巨噬细胞存在及 MMP-1 表达之间存在进展性依赖关系，提示力学应力对炎症反应具有刺激作用，从而部分解释了力学因子对斑块重塑的调节作用。

然而，上述研究多基于斑块组织学切片几何或理想化模型。组织学制片过程中的形变并未在模型中被考虑——基于组织学分析的斑块几何与体内状态可能存在显著差异（如图 7.3 中所示同一离体颈动脉斑块样本的超声影像与对应组织学切片之间可见明显形变差异）。因此，基于形变/理想化几何预测的应力分布可能无法真实反映斑块实际应力。为此，研究者开始利用现代体内/离体医学影像进行更真实的应力预测。Ohayon 等采用人冠状动脉的血管成形术前 IVUS 图像预测斑块破裂位置；Chau 等以 OCT 图像进行应力与应变研究并与组织学模型对比，发现两组模型的应力分布可比，但组织学模型的最大应力更高。Tang 等以离体 MRI 图像进行应力分析并与组织学模型对比，发现组织学模型的最大与最小应力比 MRI 模型分别高 28% 和 69%。Li 等以 5 例人颈动脉粥样硬化斑块的高分辨率体内多谱 MRI 图像建立多组分二维斑块模型，采用不可压缩超弹性材料性质分析，发现破裂斑块的平均最大应力显著高于未破裂斑块；其后续研究比较症状性与无症状性患者的二维应力集中，发现症状性患者最大预测斑块应力显著高于无症状患者，提示高应力斑块更可能在近期出现症状或破裂。Trivedi 等得到类似结果。Teng 等以 2D FEM 模型结合体内 MRI 图像的研究则表明动脉管腔曲率与纤维帽厚度应综合考虑以更好地评估斑块力学风险。

图 7.4a 展示了一例颈动脉斑块 MRI 切片；经手动分割后，将斑块几何边界点输入 ANSYS 11.0，重建斑块各组分并施加与图 7.2 类似的边界条件与材料参数，得到的静态应力分析（图 7.4b）仍呈现纤维帽高应力、脂质区低应力的格局。尽管 MRI 能提供体内几何，但基于 MRI 图像对斑块各组分（特别是薄纤维帽）进行精确分类仍是挑战。

文献中二维斑块力学模型数量众多，即便结果未必真实，其发现仍可对斑块破裂力学有所启发。二维模型的主要局限在于：(1) 平面应变/应力的假设；(2) 简化的材料模型；(3) 相对体内状态变形后的几何。这些局限可能导致斑块应力被高估，因此需要真实的三维模型才能在斑块区域获得精确的应力分析。

7.4 三维结构应力分析（3D Structure-Only Stress Analysis）

三维结构应力分析已基于理想化模型与患者特异性几何展开。

理想化 3D 斑块模型。Hoshino 等构建了理想化 3D 斑块模型以考察动脉粥样硬化钙化周围的应力分布及其对斑块易损性的影响，结果显示钙化沉积会引起壁面拉应力的局部升高，且这种升高依赖于钙化在斑块内的相对位置。Imoto 等以圆柱形血管中的理想化 3D 斑块模型，结合纵向 IVUS 斑块图像考察斑块大小、形状、膨胀性重塑、钙化与脂质核对应力分布的影响，结果显示斑块形状、大小、重塑及斑块各组分的空间构型对斑块应力有不同影响。

离体 3D 斑块模型。Holzapfel 等在病变血管中开展了一系列应力研究：他们引入了包含动脉壁组织学分层结构的多层各向异性 3D 模型，以 MR 图像构建三维有限元模型，对动脉粥样硬化髂动脉的球囊血管成形术进行仿真，并以力学测试为非线性材料模型提供原始数据；他们发现既有生物力学模型中"各向同性材料响应"与"平面应变假设"等简化假设可对应力水平产生大范围的影响，提示对斑块破裂进行应力结果解释时须为这些简化假设作充分论证。Gasser 与 Holzapfel 采用了类似的非线性材料模型，对人髂外动脉 V 型病变球囊血管成形术中的斑块裂开与剥离进行了建模。Ferrara 与 Pandolfi 以三维有限元仿真损伤动脉，并采用内聚力单元（cohesive elements）对材料拉伸强度达到时力学作用引起的断裂面进行建模，以考察几何与组织属性对斑块破裂的影响。

体内 3D 斑块模型。Ohayon 等基于 IVUS 图像计算了单支病变人冠状动脉内的应力，并探讨高应力集中与斑块破裂位置的关系——动脉壁被分为外膜、中膜、致密纤维化与细胞纤维化四种组分并采用各向同性材料，腔内施加 13.33 kPa 的生理平均血压；结果发现三维模型预测的峰值应力位置与斑块破裂位置吻合良好，但对应二维有限元模型不能吻合且常高估应力水平。Wu 等对一例体内患者特异性的、由电影血管造影图像重建的动态三维冠状动脉树进行 FEM 应力/应变分析，结果显示较小的血管直径、较大的狭窄百分比与较大的病变尺寸可能导致纤维帽内更高的应力。

相比二维应力分析，三维结构分析可在 FE 模型中纳入更多信息，如动脉壁的多层结构、考虑纤维方向的各向异性材料模型，以及去除二维 FE 模型中的平面应力/应变假设。然而仍有若干需要改进之处以达到真实预测：在生理条件下斑块承受脉动血压载荷，且在狭窄处因颈缩效应压力下降；因此腔面施加均匀压力（在 3D FE 结构分析中常用）不能模拟真实应力分布。在真实情形中，斑块的力学环境包含血流动力学与壁结构动力学，为准确预测斑块区域的生物力学环境，必须考虑斑块周围血流动力学以及脉动载荷下壁的动力学响应，因而斑块区域的三维血流仿真成为斑块应力分析的一个必要组成部分。

7.5 二维/三维 FSI 结果（2D/3D FSI Results）

本质上，动脉中的血流与动脉壁结构动力学之间存在相互作用：血流为动脉结构动力学提供载荷（如压力与壁面拖曳力），动脉壁在这些载荷下变形并改变几何形态，从而反过来影响血流。因此 FSI 方法能更真实地模拟生理环境下动脉的力学相互作用，对斑块应力分析应予采用。

Zhao 等率先将基于 MRI 的 FSI 模型用于预测健康人颈动脉的壁面剪切应力与壁应力模式，尝试定义壁面剪切应力低值与壁应力高值并存的区域；结果显示存在"低壁面剪切应力 + 高壁应力"的区域，并对应于动脉粥样硬化斑块好发部位。Younis 等用流-固耦合有限元仿真考察了三名健康人颈动脉分叉的血流动力学与壁力学差异——基于 MR 几何与超声测得流边界，得出壁面剪切应力、振荡剪切指数与循环应变等参数。Kaazempur-Morfrad 等以循环应变描述脉动血压下动脉壁动力学，发现循环应变变化最大的区域与动脉粥样硬化好发部位吻合；其后续研究对 4 例患者进行应力分析，并将切除的斑块切片染色以分析平滑肌细胞、巨噬细胞、脂质与胶原，尝试将流体动力学参数与动脉粥样硬化的组织学标志物相对应。

Li 等以二维血流-斑块相互作用模型探讨纤维帽厚度对颈动脉斑块稳定性的关键程度，结果显示中度狭窄（30%–70%）的颈动脉斑块若纤维帽较薄，仍可能具有较高的破裂风险。Li 等构建了狭窄程度分别为 30%、50%、70% 的理想化 3D 斑块模型以考察斑块喉部壁运动，结果显示重度狭窄可能抑制壁运动。Bluestein 等在偏心狭窄模型中数值考察了微钙化在斑块易损性中的作用，结果显示含嵌入式钙化斑点的易损斑块在纤维帽内呈现更高的壁应力集中区，为"微钙化增加斑块易损性"的假说提供了证据。Kock 等基于两例症状性动脉粥样硬化患者的体内 MR 图像构建了 2D FSI 模型，计算了每例患者的纵向应力水平，认为纤维帽纵向应力可能有助于评估斑块易损性。

Tang 等以人颈动脉 3D 离体 MR 图像构建了 3D FSI 模型以识别可能与斑块破裂相关的关键流体与应力/应变条件；后续多年里 Tang 团队从离体到体内开展了一系列斑块应力分析，结果显示大脂质池与薄纤维帽与极端应力/应变水平相关，脉动压力引起的大循环应力/应变可能通过疲劳机制导致斑块破裂；与 2D 分析一致，其 3D FSI 模型亦显示脂质池大小、形状、位置、纤维帽厚度与轴向拉伸等因素会显著影响斑块应力与应变。Yang 等基于体内/离体 MRI 的 3D 非牛顿 FSI 模型与"流体-壁"或"纯壁"模型的对比研究表明，两类模型在壁面剪切应力与最大主应力上可能存在显著差异。Huang 提出了一种基于体内 MRI 数据、采用患者特异性血管轴向与内周向收缩步骤确定颈动脉粥样硬化斑块 3D 零应力状态的方法，认为准确把握血管收缩过程可显著提高体内斑块模型的应力预测精度。Huang 等基于 10 例斑块模型的体内 MRI 3D FSI 模型，发现斑块内出血与更高的结构应力相关。Creane 等认为斑块内外表面曲率差异也是症状性斑块发生发展的重要因素。

将斑块应力与斑块进展/破裂关联的研究，将有助于理解生物力学在动脉粥样硬化斑块起始、进展与破裂中的作用，并为治疗方案提供依据。图 7.5 展示了以患者特异性颈动脉斑块几何进行 FSI 应力分析的实例：图 7.5a 为收缩末期腔面压力分布（因中度狭窄，压力降在斑块区域并不显著），图 7.5b 为壁面剪切应力分布（斑块区域略高），图 7.5c 为以第一主应力表示的斑块壁拉应力分布；可在斑块下游侧观察到局部高应力集中区，被认为与斑块破裂相关。从图中可得血流动力学与斑块结构动力学两方面的应力因子，可为分析应力相关斑块破裂提供更完整的框架。

7.6 模型重建（Model Reconstruction）

基于二维/三维医学影像的计算生物力学仿真为研究生物力学因子与动脉疾病之间的关联提供了新途径。在开展斑块生物力学研究时，图像采集之后的第一步就是二维或三维的斑块模型重建。然而重建过程存在若干潜在不确定性（尤其对三维模型）。已有研究显示，斑块几何的精确重建会显著影响生物力学结果；应投入更多精力以提高几何重建的精度与详细程度，因为这是对生理学显著指标的首要影响。

随着现代医学影像技术的发展，越来越多的体内图像可被用于几何重建。患者特异性病变动脉模型的几何重建已被研究多年，方法成熟——既可使用第三方软件进行图像分析与几何提取，也可基于已建立的图像处理技术开发内部软件。现有方法大多基于以下流程：提取动脉横截面边界—在三维上下对齐这些边界—由二维边界序列构建三维表面与实体。图 7.6 展示了作者课题组以体内 MRI 颈动脉斑块几何进行斑块几何重建的一般流程。

Moore 等以颈动脉分叉体模与人颈动脉分叉考察了三维腔壁几何重建的整体精度，并比较了五种不同重建流程（含/不含平滑处理），结果显示对体模所有重建都可接受，但体内模型需适当平滑后方可用于计算体内血流动力学研究。Long 等考察了 8 例人颈动脉分叉的体内 MRI 三维几何重建可重复性：以平滑后的腔轮廓沿纵轴对齐，以 B 样条插值重建三维表面并采用中心线平滑与表面平滑两步平滑；结论是大多数情况下颈总动脉几何可被良好重建，而颈外动脉的可重复性最差。3D 超声也已用于颈动脉几何重建以支持基于图像的 CFD 建模。

然而，多数研究聚焦于动脉腔的重建及其对血流动力学因子的影响。斑块区域的应力分析要求各组分的精确重建，这比单纯重建动脉腔更困难。斑块结构分析涉及的不同组分在尺寸、形状上各异，尤其在薄纤维帽的重建上。一例易损斑块可能具有 65 μm 以下的薄纤维帽——在医学图像中可能仅由一两个像素表示，因此薄纤维帽重建的相对不确定性会高于其他斑块组分。

尽管结构应力分析已应用现代医学影像在斑块易损性评估中多年，关于斑块几何重建不确定性及其对应力分布影响的研究仍很少。为分析从 MRI 图像处理到三维重建再到最终斑块应力分析中的不确定性，Gao 等研究了 3 例颈动脉斑块以估算可能的不确定性，发现基于 T2W 图像的动脉壁重建在操作者之间的可重复性高（差异在 1 像素水平），但脂质核与纤维帽分割的可重复性较低；3D 表面插值、平滑与分叉三维重建所引起的不确定性小到可忽略；应力分析显示，脂质区高估会诱导斑块内更高的应力集中与整体应力水平；动脉壁高/低估对应力分布的影响则显著小于脂质区高/低估造成的影响。新 MRI 协议下斑块成像分辨率/质量的提升将带来更真实的应力预测，但 MRI 图像采集过程中的不确定性在该研究中尚未被评估——这同样重要。

7.7 材料属性（Material Properties）

对健康与病变动脉组织的力学属性进行定量是获得真实应力预测的前提。组织学研究已显示动脉壁具有非均质、非线性、各向异性、黏弹性的特征。Richardson 指出，准确的斑块材料属性数据缺失是现有斑块破裂研究中可能最不确定的方面。斑块的主要组分包括脂质核、钙化与纤维帽组织。虽然健康动脉组织的材料研究较多，但斑块各组分的力学属性研究则少得多；已有数据综述显示，斑块组分材料属性的范围非常广——总体而言脂质比其他组分软得多，钙化则硬得多。纤维帽的刚度对斑块稳定性至关重要，由于复杂的纤维组织与组分变化，纤维帽的刚度变化范围很大。因此，获取合适的材料数据是结构应力分析（尤其针对纤维帽组织）的关键步骤。

Lee 等以静态与动态加载试验确定了人腹主动脉斑块纤维帽的力学属性，发现纤维帽的刚度取决于其微观结构（钙化、胶原含量与组织）；纤维帽刚度随加载频率增加而增大，表明纤维帽刚度本质上是非线性的。Loree 等考察了人动脉粥样硬化组织的静态周向模量，发现静态周向切线模量受细胞性与钙化程度的影响不显著（与压缩模量不同）。Loree 等以扭转流变仪对不同脂质成分的脂质池进行力学测试，发现脂质池的刚度取决于胆固醇一水合物晶体的浓度。Topoleski 等以定制实验系统对全人动脉粥样硬化斑块节段进行单轴径向循环压缩测试，发现斑块在有限变形下呈非线性特征且与斑块组分与加载历史相关。Topoleski 与 Salunke 对钙化斑块进行压缩与应力松弛实验，发现钙化斑块与其他斑块显著不同。Holzapfel 等测定了人动脉粥样硬化斑块组织中各组分的各向异性力学属性（含极限拉伸应力/应变），发现斑块组分高度非线性且在不同层之间不同，纤维帽周向的断裂应力最低。Ebenstein 等以纳米压痕测定人动脉粥样硬化斑块组织中的血凝块、纤维组织与钙化纤维组织的力学属性，发现斑块组织刚度随矿物含量增加而升高。Barrett 等以实验结果拟合有限元仿真测定了 8 例人颈动脉动脉血栓斑块样品的压痕响应。Maher 等测定了新鲜人颈动脉粥样硬化斑块的拉伸与压缩属性，发现钙化斑块的刚度响应最大。

现代医学影像近期也被用于定量病变血管的材料属性，为在体内条件下表征材料属性、实现真实患者特异性应力仿真提供了可能。基于 3D IVUS 模型，Vonesh 等提出了一种以 FEM 估算病变髂动脉与股动脉区域材料属性的方法，结果显示非病变组织的弹性模量与斑块组织显著不同。Karimi 等将非线性 FEM 与弹性成像结合，提出了一种估算软组织非线性弹性属性的新方法，已应用于真实 2D 与理想化 3D 动脉斑块模型，并证明了估计斑块内非线性材料属性空间分布的可行性。Masson 等基于无创体内临床数据（动态腔内压、中膜直径、内膜-中膜厚度）提出了一种确定人颈总动脉材料属性与壁应力的方法，为直接从无创人体数据刻画患者特异性材料参数提供了途径。

7.8 残余应力/应变（Residual Stress/Strain）

学界已普遍接受"生物组织在所有外部载荷移除后并非处于无应力状态"。已有研究显示这种残余应力/应变有助于维持血管结构的完整性，并使应力在动脉壁内更均匀分布。但残余应力对斑块稳定性的影响尚未被充分研究。Ohayon 的研究表明斑块内残余应力/应变不可忽略，可能显著影响薄纤维帽内的生理峰值应力，并提示斑块破裂应被视为外部载荷与斑块内残余应力/应变的组合。作者课题组的研究则显示在正常生理载荷下周向残余应力对实际应力分布的影响有限。这些相互矛盾的结果表明，仍需投入更多努力以理解残余应力/应变对斑块稳定性的作用。

7.9 讨论与结论（Discussion and Conclusion）

动脉粥样硬化斑块破裂与继发血栓形成被认为是全球范围内主要致死原因之一。从生物力学角度看，斑块破裂可被视为一种力学事件，因此以高精度预测斑块应力至关重要，结果有助于理解斑块如何破裂。多年来斑块应力研究已从二维结构分析发展到结构动力学与血流动力学完全耦合的三维斑块应力；斑块区域内的极端应力位置也被视为斑块破裂的主要因素。当前最大的挑战在于如何对患者特异性仿真进行准确的应力预测，例如真实几何重建、边界条件与材料模型。已有研究表明二维斑块模型因假设的原因倾向于高估应力水平；三维 FSI 斑块应力分析的出现提高了分析精度，不仅能分析斑块结构内的应力，也能分析斑块周围血流模式。

现代医学影像的发展为高分辨率真实三维斑块几何的获取提供了可能，斑块各组分（纤维帽、脂质区）可由多谱 MR 图像等获得。然而 MR 图像分辨率不足以重建很薄的纤维帽，且对比度也不高，对基于体内 MRI 的斑块应力分析带来一定不确定性。关于斑块各组分对斑块应力分布影响的研究显示，相比脂质大小，斑块应力水平更显著地受纤维帽厚度影响；因此近期仍需努力以提高基于体内 MRI 的斑块重建精度。基于斑块几何的应力分析精度受 MR 图像质量影响，新的 MRI 协议下斑块成像分辨率/质量的提升将产生更真实的应力预测。虽然 FSI 仿真能比 2D/3D 结构应力分析提供更真实的应力预测，目前建模技术仍存在若干主要局限：(1) 真实材料属性建模；(2) 血流模型的选择（如层流假设）；(3) 患者特异性边界条件，尤其是斑块腔内的真实压力；(4) 残余应力/应变的纳入。因此在将斑块应力结果与患者实际情况对接时需注意，并应承认其不确定性。

纤维帽的破裂是一个非常复杂的过程，不仅涉及力学因子，组织与细胞层面的异常（如升高的炎症活动与退化的胶原结构）也可能影响破裂过程。综合性的斑块易损性评估应同时考虑系统性标志物、形态学特征与生物力学因子，仅考虑少数因子可能无法获得斑块破裂风险的完整图景。

作者在结论中指出，表 7.1 汇总了近年来斑块应力分析的主要研究。数值建模在应力分析方面的进展已使基于患者特异性斑块模型的流-固耦合分析成为可能，斑块区域内的极端应力分布可被预测并用于评估斑块破裂风险，这些评估对临床实践具有帮助。结合斑块 MR 成像分析、计算建模与临床研究/验证将增进对斑块破裂的理解，并建立患者诊断、管理与治疗的新流程。

本章个人批注

本章呈现的是基于有限元的斑块应力分析综述——一条"假说 → 模型 → 真实几何 → 真实材料 → FSI"的递进脉络。值得思考的几点：

"局部最大应力"作为主流假说的实证基础（300 kPa 临界值，Cheng et al. 1993）目前主要来自体外/离体组织学几何模型——正如 7.3 节自陈"破裂斑块不一定与高狭窄相关"。这与本书 Ch04 提到的"高危斑块特征"在更广义的形态学上呼应，但 7.3 节也明确指出组织学几何在体内外的形变会高估应力。从 4 节的临床可操作角度看，300 kPa 阈值仍是一个"理论阈值"——它假设了纤维帽材料均匀、几何未变形、平面应变近似成立，所有这些假设在体内都不严格满足。
二维 → 三维 → FSI 的方法学递进其实是一种"承认局限"。7.3 节末尾明确列出二维模型的三大缺陷（平面应变假设、简化材料、形变几何），7.4 节末又指出"在 3D FE 结构模型中常用的腔面均匀压力"也不能反映真实情形，7.5 节才把 FSI 推到前台。这与 Ch06 从 b-SSFP/PC/perfusion/DE/T2* 模块化看待 CMR 工具的思路是相通的——单一模态都有局限，多模态耦合/多模块联合才是更接近生理的方案。
7.7 节"材料属性是最大的不确定性"这一点尤其关键——大多数临床可获取的体内图像（IVUS、MRI）能给出几何，但材料参数只能查表或基于离体测试。这与 7.6 节"几何重建受限于薄纤维帽<65 μm（1–2 像素）"一起，构成"双层不确定性"：几何不确定性 + 材料不确定性。两者之中哪一个是当前研究瓶颈？作者没有直接回答，但从文章叙述次序看，作者课题组更倾向于几何侧（7.6 占了较大篇幅，且 7.6 提到"高估脂质区会显著影响应力"）。
残余应力/应变两派意见（Ohayon "不可忽略" vs 作者课题组 "有限影响"）非常有趣。它提示我们：在缺乏体内测量的情况下，"残余应力的影响"本身是模型参数之一——你选择不同的本构模型/边界条件，就可能得到不同结论。这是 7.8 节一个未解决的问题。
关于 7.6 与 7.7/7.8 的归类：原文将 3.1/3.2/3.3 作为"3. 影响 FEM 仿真精度的因素"的三个子节。按"参数"维度看，这三者其实并不完全平行——3.1（几何）和 3.2（材料）是输入侧的"源数据"，3.3（残余应力）则是材料本构的一个特殊维度（预应力/预应变）。这种归类反映了 FEM 仿真的"任务流程"视角（从几何到材料再到本构细节），而非严格的物理维度分类。

值得核对的疑问： - 7.2 节"5. Other hypotheses"中提到"injury due to turbulent flow in the stenosis [43], rupture of vasa vasorum [2]"——前者被排在"破裂假说"序列里似乎有点奇怪（湍流更接近血流动力学现象），可能是作者对"破裂诱因"的广义化。 - 表 7.1 引用 [70] 同时被列在 2D 与 3D 应力分析两个区域——可能是排版/分类上的归并不一致，需要看原文。 - 7.3 节提到 "stress in the stress concentration regions can vary up to 30 % with different material models from the same work"——这与 7.7 节"材料不确定性是最大的不确定性"在数量级上的对比，30% 的差异看起来并不大；可能因为 Williamson 等使用的是较为接近的同类材料模型组（各项同性非线性、含残余应变、横观各向同性线性），而真正的"非均质非线性黏弹"差异会大得多。

与上下章的衔接（一段话）

承接 Ch04（颈动脉 MR 定量分析）和 Ch05/06（对比剂 / CMR 工具）建立起的"成像 → 定量"基础设施，Ch07 把视野切换到"基于影像的力学建模"——这是把前面章节获取的几何与组分信息（管腔形态、脂质大小、纤维帽厚度等）转化为"生物力学标记"的桥梁。本章 7.1 节明确说明"体内无法直接测量斑块壁应力"，因此基于 FEM 的数值方法是目前唯一可行的路径——这构成了 Ch08（患者特异性斑块应力）与 Ch09（MR 识别易损斑块）的理论支撑。从本书的整体结构看，Ch07 是"定量 MR/CMR → 力学建模"的过渡章节，承上启下：前 6 章是成像与定量工具的铺垫，Ch07 是工具向"斑块易损性/破裂风险"这一临床终点对接的关键中介，而 Ch08-09 则把这种中介扩展到"患者特异性 FSI"和"基于 MR 的易损性识别"两个具体方向。