《Annals of Biomedical Engineering》：The Effect Of Localized Circumferential Residual Stress On Pressurized Stress State Of CT Reconstructed Vessels: A Finite Element Study
编辑推荐：
为解决患者特异性血管模型中残余应力（Residual Stresses）难以精准纳入及解剖真实性不足导致应力预测偏差的问题，研究人员开展猪主动脉局部周向残余应力有限元（FE）模拟研究，结合体外成像、力学测试与FE模拟，发现残余应力可均化血管壁应力分布、降低 lumen 应力集中，使120 mmHg下lumen横截面积与实验测量更吻合， neglect残余应力会导致inner wall应力峰值高估50%，为手术规划与风险预测提供关键生物力学依据。
当我们谈论血管疾病的风险评估或手术规划时，往往依赖计算机模拟来预测血管壁在生理压力下的应力状态——毕竟，直接对人体血管做“压力测试”既不现实也不道德。但你可能不知道，这些模拟里藏着一个容易被忽视的“误差源”：残余应力。简单来说，残余应力是血管在没有外部压力（比如血压）或载荷时，因生长发育、组织结构不均（比如弹性蛋白和胶原蛋白的非均匀沉积）而“自带”的内部应力。就像一根被弯折后松开的弹簧，即使没挂重物，内部也还留着“想弹回去”的力。
过去的研究要么用简化的圆柱几何模型（把复杂的主动脉当成直管子），要么用通用的解析模型估算残余应力，根本没考虑到真实血管的解剖细节——比如不同部位的壁厚差异、周向开口角的位置变化。更关键的是，没人真正量化过：当把这些局部周向残余应力放进患者特异性模型里，到底会让应力预测变准多少？ 这就像你用一张模糊的地图导航，却不知道某个路口其实有个“隐藏弯道”——误差可能大到影响对血管破裂风险的判断。
为了填补这个空白，一组研究人员决定“较真”：他们用猪主动脉做样本，结合体外成像、力学测试和有限元（FE）模拟，要把局部周向残余应力“装进”真实的血管模型里，看看它到底怎么改变加压后的应力状态。这项研究最终发表在《Annals of Biomedical Engineering》上，给血管生物力学模拟上了一堂“精准课”。
关键技术方法
研究采用单只杂交大白猪的主动脉（从主动脉瓣到腹腔干），术后5小时内完成实验；通过CT扫描获取加压（120 mmHg）和纵向切割后2小时的零应力状态图像；从主动脉5个轴向区域提取30条环状/纵向组织条，做位移控制单轴拉伸试验（5次循环至40%应变），用Holzapfel–Gasser–Ogden本构模型拟合区域特异性正交各向异性材料参数；用3DSlicer分割开放主动脉的立体 lithography文件，经ANSA预处理生成四面体网格（敏感性分析确定0.50 mm单元尺寸最优）；通过LS-DYNA求解器模拟切割面闭合（位移边界条件+绑定接触）以恢复残余应力，再对含/不含残余应力的闭合模型施加120 mmHg腔内均布压力，对比lumen横截面积与应力分布。
研究结果
Open Configuration of the Aorta
猪主动脉纵向切割后2小时，开放构型的开口角沿切割位置变化：近端峰值达78.6°，中央区约55.0°，远端约25.0°；切割末端因残余压应力导致血管弯曲、lumen闭合，避免了容器底部接触对开放构型的干扰。
Circumferential Residual Stress Estimation
通过FE模拟闭合开放主动脉，恢复的周向残余应力呈典型分布：lumen附近为压应力，outer wall为拉应力；中央区域应力分布符合σfitting= a/r + b（r为径向距离），拟合参数见表3；切割面附近的应力异常由接触算法数值噪声引起，近切割端的残余应力因未切割区约束未完全发育。
Pressurization
- •
Lumen Cross-Sectional Area at 120 mmHg：含残余应力时，80 mmHg下模拟与实验lumen横截面积平均绝对百分比差为5.99±4.13%，120 mmHg时降至2.85±1.73%；不含残余应力时，120 mmHg下差异增至5.57±3.49%，且模型预测的lumen更大——含残余应力让结果与实验更吻合。
- •
First Principal Stress Comparison：第一主应力与周向一致。不含残余应力时，inner lumen应力峰值最高，沿壁厚递减（符合Laplace定律）；含残余应力时，应力沿壁厚更均匀：80 mmHg时应力差减少73.40%，120 mmHg时减少52.62%。两种情况的百分比差异达±50%（80 mmHg）和±25%（120 mmHg），残余应力显著降低了inner wall的应力峰值。
研究结论与讨论
这项研究首次将局部周向残余应力整合到具有解剖细节和区域特异性力学性能的血管FE模型中，核心结论直击过去模拟的“痛点”：
- 1.
残余应力的“均化魔法”：残余应力能把血管壁内的应力“摊平”——原本集中在lumen附近的峰值应力被分散，outer wall的应力被提升，让整个壁厚的应力分布更均匀。就像你把堆在一角的书均匀铺开，每个书架的承重更合理。
- 2.
忽略残余应力的代价：如果不考虑残余应力，inner wall的应力峰值会被高估50%（80 mmHg时），lumen横截面积的预测也会偏差更大——这对于需要精准评估血管破裂风险的手术（比如主动脉瘤修复）来说，简直是“致命误差”。
- 3.
解剖细节的重要性：研究发现，开口角、血管壁厚、材料参数的区域差异，都会直接影响残余应力的分布。过去用简化圆柱模型的方法，根本没法捕捉这些真实解剖特征带来的应力变化。
当然，研究也有局限：比如只用了单只猪的主动脉，没考虑轴向残余应变，也没区分血管的内膜、中膜、外膜层——但这些都不影响它的核心价值：证明了在患者特异性血管模型中纳入局部周向残余应力的必要性。
正如论文讨论中所说：“finite element simulation在in-silico医学中常用于评估患者血管的破裂风险或对医疗设备的反应。我们的结果显示，残余应力能让峰值应力降低50%（舒张期）和25%（收缩期）——忽略它，可能会高估血管损伤风险。” 换句话说，这项研究给血管生物力学模拟装了一个“精准滤镜”，让未来的手术规划和风险评估更靠谱。
从“模糊的圆柱模型”到“带局部残余应力的解剖真实模型”，这一步跨得不大，但对每一个需要用模拟救命的患者来说，可能就是“精准”与“误差”的天壤之别。
