A Stabilized Mixed Hybrid Finite Element Framework for Large-Deformation Poromechanics with Torsion: Application to Transmural Perfusion and Vascular Hemodynamics

人类的心脏是一个不仅是单纯挤压血液的泵；在每一次跳动中，它都在扭转、拉伸并改变自身形状。在左心室厚实的肌肉壁内部，上演着复杂的流体动力学问题。心肌并非一块坚实的组织，而是一种类似海绵的材料，大约有百分之十五的空间充满了在微小通道中流动的液体。随着心脏充盈血液并随后通过扭转来射血，这种内部液体被挤压和推动，产生出的压力既可能有助于也可能阻碍血液流经心脏自身的微细血管。如果肌肉内部的压力过高，可能会压碎这些脆弱的血管，导致心肌组织缺氧。理解心脏的扭转运动如何与这种内部流体流动相互作用至关重要，然而传统的计算机模型往往难以准确模拟这一过程，经常产生令人困惑的错误，或在组织几乎无法被压缩时出现“锁定”现象。

埃因霍芬理工大学的一位研究人员开发了一种模拟这种行为的新方法，创建了一个数学框架，能够捕捉心壁的大规模扭转和拉伸，而不陷入这些常见的陷阱。通过将心肌视为一种流体饱和的海绵，并使用一种专门的方法分别追踪固体组织和流体的运动，该研究揭示了心脏自然的扭转运动如何起到一种保护机制的作用。模拟显示，当心脏在充盈期间允许自由扭转时，它能显著减轻内壁纤维的应力，并降低威胁要使血管塌陷的内部压力。相反，当心壁由于高血压而变得异常肥厚时，液体会被困住，压力会危险地升高，而这正是标准计算机模型经常无法正确预测的情况。

这项工作的核心在于解决了一个长期存在的难题，即如何模拟像心肌这样柔软且潮湿的材料。在标准的计算机模拟中，当一种材料接近不可压缩（即无法被挤压成更小的体积）时，数学计算往往会崩溃，产生不真实的压力峰值或锯齿状的振荡结果，这并不符合现实。这被称为体积锁定。为了解决这个问题，研究人员制定了一种“混合混合”（mixed hybrid）方法。该方法不再试图一次性计算所有内容，而是将问题分解成较小的部分，并执行严格的规则，确保流体能够连续地跨越每个部分的边界，正如在真实人体内必须发生的那样。这种技术允许计算机分别求解组织运动和模拟块边缘的压力，从而确保流体流量得到精确守恒，并确保即使在极端变形下，结果也能保持平滑且符合物理实际。

为了测试这个新框架，研究人员首先推导出了一个关于充气并扭转的厚壁圆柱体的精确数学解，作为完美的基准。随后，他们将该方法应用于一个左心室模型，这是一个尺寸符合人体解剖结构的圆柱形心肌部分。模拟探索了三种特定的场景。首先，他们观察了心脏自身冠状血管内不同水平的血压如何影响壁的刚度。他们发现，当血管内血液充盈良好时，内部流体压力会上升，这实际上会使放松状态下的心壁变硬。这种现象被称为“花园水管效应”，意味着心脏需要更高的压力才能充盈，从而限制了其摄取血液的容量。新模型准确地预测了这种转变，显示出高灌注压力将肌肉壁中的固体应力从十千帕降至三千帕，有效地实现了流体与组织之间的载荷分配。

第二个场景考察了心脏扭转运动的关键作用。在健康的心脏中，肌纤维呈螺旋状排列，从内表面延伸至外表面，并随着路径变化角度。当心脏充盈时，这种螺旋结构会导致心壁发生扭转。模拟对比了允许自由扭转的心脏与被固定住不动的硬化心脏。当扭转被阻止时，心壁内侧肌纤维的应力飙升至 57.20 千帕，这是一种危险的力集中。然而，当模型被允许自然扭转时，峰值应力下降了 45.5%，降至 31.20 千帕。这种力的重新分布也有助于流体更容易地通过心壁排出，将内部孔隙压力从大约 128 毫米汞柱降低到 122 毫米汞柱。研究表明，心脏的扭转不仅是收缩的副产品，更是缓解应力并维持肌肉内层血液流动的必要机制。

第三个也是最重要的案例研究，探讨了在高血压压力下承受重负的心脏，此时心壁会显著增厚以承受高泵压。在这种场景下，心壁的内半径建模为 2.00 毫米，外半径为 5.00 毫米，腔体压力为 180 毫米汞柱。在这些条件下，增厚的壁体会导致液体被困在靠近内表面的地方，因为液体向外逃逸的距离变长了。新的混合混合模型显示，这种被困的液体会在初始加载后仍持续产生高压区。相比之下，当研究人员使用标准计算机建模技术运行同一场景时，结果却是混乱的。标准方法在内壁附近产生了非物理性的压力伪振荡，并且未能捕获真实的应力分布，导致误差率高达 12.75%。然而，新方法产生了一个平滑的单调压力曲线，误差仅为 0.01%，正确识别出厚壁结合扭转运动会形成一个机械障碍，从而减缓排水速度并维持高内部压力。

这些发现表明，心脏扭转和复位的能力是其水力健康的基础。当心脏自由扭转时，它使应力在壁上均匀分布，并促进流体排出，从而保护微血管免于塌陷。当心壁增厚或扭转受限时，液体会被困住，导致高压，进而危及血液流动。研究证实，这一新的计算框架能够捕捉这些微妙但关乎生命的相互作用，为理解心脏中机械力与流体流动的耦合提供了一个可靠的工具。虽然目前的工作侧重于使用简化的圆柱几何形状来证明概念，但该方法在处理大变形和复杂流固耦合方面的成功，预示着未来此类模型可以应用于更复杂的、三维的心脏形状，以更好地理解心力衰竭和高血压的影响。