《Journal of Biomedical Materials Research Part B: Applied Biomaterials》：Computational Performance Evaluation of Porous Ti-6Al-4V Scaffold Architectures for Femoral Defect Healing Using a Coupled Finite Element–Mechanoregulation Framework
编辑推荐：
增材制造的多孔钛支架可为大段骨缺损提供固定并引导再生，但设计选择受制于组织形成与疲劳驱动失效在康复载荷下协同演变的不确定性。研究人员重建了一例患者特异性股骨干缺损，并以Ti-6Al-4V支架填充，涵盖支柱型、片状三周期极小曲面（TPMS）及随机型三大类共16种
增材制造的多孔钛支架可为大段骨缺损提供固定并引导再生，但设计选择受制于组织形成与疲劳驱动失效在康复载荷下协同演变的不确定性。研究人员重建了一例患者特异性股骨干缺损，并以Ti-6Al-4V支架填充，涵盖支柱型、片状三周期极小曲面（TPMS）及随机型三大类共16种拓扑结构，所有设计均设定为75%孔隙率及1 mm最小特征厚度。耦合有限元—力学调控框架模拟了240天的愈合过程：有限元场驱动基于格点的力学调控模型，以八面体剪应变为刺激指标，同时支架疲劳损伤改变演化的力学环境。在保守与激进两种康复方案下，13种拓扑结构保持结构完整性，并收敛至接近成熟组织的密度（约1.90–1.91 g/cm3），且桥接模式相似，表明在力学稳定性得以维持时预测成骨进展一致。失效具有拓扑及方案依赖性：Lidinoid在两种方案下均发生坍塌，而Gyroid与SplitP仅在激进载荷下失效，疲劳损伤集中于缺损中段区域，并沿连续片状连通路径扩展。在存活的各设计中，面心立方、简单立方、Octet及Voronoi属于有利的低损伤组；终点密度的亚百分比差异未能分辨出唯一的生物学最优结构。研究人员结合文献对比对预测的愈合趋势进行了情境化解读。总体而言，该框架为识别易失效拓扑结构并优先筛选力学稳定的候选支架以开展后续匹配实验提供了计算筛选基础。
多孔Ti-6Al-4V支架拓扑结构影响股骨缺损愈合的计算评估：耦合有限元—力学调控研究解读
一、研究背景与问题提出
股骨干大段骨缺损是骨科重建中持续存在的临床难题，其内在愈合能力往往不足以完成修复。重建策略既需承担力学载荷，又需支持生物学再生。股骨内部受力呈多轴性且随步态变化，即使轻微的不稳定也可能在缺损区域形成不利的力学微环境。骨再生具有力学敏感性，局部力学刺激直接影响细胞行为与组织分化方向。因此，候选治疗方案需同时评估结构力学性能及其对愈合环境的调控作用。
增材制造技术使患者特异性多孔钛支架成为可能，其内部几何结构可被精确控制。现有设计涵盖支柱型周期点阵、三周期极小曲面（TPMS）架构及随机Voronoi型网络。在相同孔隙率下，这些家族在连通性、各向异性、曲率及载荷路径组织方式上存在显著差异，进而改变局部应力与应变分布。因此，架构比较不能仅依据孔隙率推断，而需进行直接对比。
单纯的结构指标如峰值应力、刚度或表观强度无法全面描述再生性能。愈合过程具有力学调控特性，力学调控算法可将局部力学条件与随时间变化的组织分化过程联系起来。将三维力学细节与力学调控模型耦合，可将架构特异性力场转化为预测的骨形成或表观组织密度，从而建立统一的比较框架，避免大量试错式实验。
然而，现有研究多将结构与生物学终点分开考察，或使用理想化试样与节段级几何，或仅考虑单一拓扑家族。在同一患者特异性股骨、同一固定系统、孔隙率、材料模型、载荷历史及力学调控规则下进行匹配的跨家族对比仍属空白。本研究正是基于这一缺口，开展大规模计算筛选，系统比较16种支架架构在两种康复轨迹下的力学—生物学耦合表现。
二、研究内容与核心结论
该研究构建了患者特异性股骨干节段缺损模型，基于53岁男性成年受试者的CT影像重建股骨几何。缺损长度为3 cm，支架设计域对应缺损区域内缺失的皮质骨体积。固定支架整合于缺损边缘，每个支架设计两个皮质螺钉孔位，并采用锥形轮廓增加支架—支架界面面积。支架厚度设为5 mm，与支架一体增材制造，赋予增材制造Ti-6Al-4V材料属性。
研究设置了两个全局约束以隔离架构效应：目标孔隙率统一为75%，最小特征厚度统一为1 mm。前者为文献支持的中间值，兼顾细胞存活与骨长入需求；后者为保守的建模与制造选择，避免向更小特征外推材料性能。各架构家族以单一主控制变量参数化以满足孔隙率目标：支柱型点阵校准单胞尺寸，TPMS校准隐式曲面单胞尺寸，随机网络校准点间距。
力学分析采用单一准静态足跟触地载荷工况，包含髋关节反力、外展肌、髂腰肌及股四头肌四个主要力贡献者。远端股骨髁施加固定约束。体重设为60 kg，所有载荷按体重倍数定义。康复方案基于术后靴内传感器数据构建：保守方案采用各时间段内步数与每步体重百分比的最佳轨迹，激进方案采用上偏差包络以刻意加压测试支架性能。
力学调控模型基于格点框架，以八面体剪应变作为旋转不变标量刺激指标。表型窗口预先采用文献中的应变比较器设定，而非拟合本研究结果。组织分化遵循纤维组织、软骨、骨及骨吸收的应变阈值区间。细胞过程以每日概率规则推进，有限元问题每5个生物步重新求解。组织力学性能通过成熟度加权插值从基线向表型特异性目标演化。
疲劳损伤在支架单元层面以标量累积损伤变量追踪，采用Miner型线性累积法则。损伤评估每5天进行一次，利用稳定循环材料参数与时间依赖性康复载荷。失效单元在后续求解中保持失活状态，使载荷重新分配至剩余支架网络。
三、主要技术方法
研究人员采用耦合有限元—力学调控框架，将ANSYS Mechanical中的连续介质有限元分析与外部C++格点生物学引擎交替耦合。有限元模型以二次四面体单元离散支架与骨结构，以线性四面体单元离散组织域。支架网格目标尺寸0.5 mm以捕捉1 mm特征内的梯度，组织表面网格0.5 mm、内部1.5 mm。力学调控模型采用0.08 mm格点间距，每个有限元单元包含预设格点子位点，以单元级等应变近似传递刺激。疲劳参数从文献稳定循环响应数据按厚度比缩放转移至1 mm薄壁状态。耦合循环每5天交换力学场与组织状态，直至240天终点。
四、研究结果
4.1 组织密度演化
在两种康复方案下，16种架构中有13种完成全部模拟且未发生结构失效，第240天达到接近成熟组织的密度值，通常位于1903–1911 kg/m3范围。三种失效架构为Lidinoid、Gyroid（仅方案B）及SplitP（仅方案B）。早期阶段（约前50天）所有架构遵循相似的S形密度轨迹，无明显架构依赖分离。第60天存活支架组织密度约547–633 kg/m3，反映成骨细胞介导成熟的启动。第60至120天密度累积加速，第120天所有存活架构达到1307–1366 kg/m3。第180天后曲线收敛至约1900–1910 kg/m3平台，接近成熟骨密度目标1920 kg/m3。
存活架构的终点密度形成狭窄聚集簇而非生物学上有意义的排序。方案A均值为1908.4±2.3 kg/m3，范围为1902.4–1910.4 kg/m3；方案B均值为1908.1±5.0 kg/m3，范围为1894.7–1911.4 kg/m3。亚百分比差异在数值、随机及模型精度内无法分辨，不被解释为某一存活架构产生更多骨组织的证据。
4.2 空间骨形成模式
代表性有利架构（简单立方、面心立方、Octet、Schwarz及Voronoi）在第60、120、180及240天的密度场显示，组织长入从宿主组织界面起始，随时间向缺损中心推进。存活架构的桥接模式相似，表明力学稳定性维持时成骨进展具有一致性。失效架构的组织长入从界面推进但终止于坍塌前的未完成状态。
4.3 疲劳损伤演化
累积支架失效体积分数显示，Lidinoid在方案A第75天、方案B第45天发生结构坍塌。Gyroid在方案A存活至第240天，方案B第100天坍塌。SplitP在方案A完成模拟，方案B第130天失效。失效架构的损伤空间分布显示，高损伤区域在缺损中段区域形核并沿连续片状连通路径快速扩展，最终汇聚导致结构坍塌。
4.4 非灾难性损伤架构
部分架构虽存在可测量的单元失效但未导致整体失效，包括体心立方、Fluorite、Schwarz及Voronoi。其损伤局限于离散支柱节点或表面区域，解释了这些支架虽存在非零单元失效但保持整体结构完整性的原因。面心立方、简单立方、Octet及Voronoi属于有利的低损伤组。
五、讨论与结论
该框架揭示了生物学收敛但力学完整的设计与疲劳坍塌中断愈合轨迹的拓扑之间的明确区分。失效具有拓扑及方案依赖性，Lidinoid在两种方案下均失效，而Gyroid与SplitP仅在激进载荷下失效。疲劳损伤集中于缺损中段并沿连续片状连通路径传播，表明TPMS家族的连续片状连通性在损伤传播中起关键作用。
存活架构间终点密度的亚百分比差异无法分辨唯一生物学最优结构，提示在力学稳定性维持的前提下，多种拓扑均可支持近成熟组织密度的达成。该框架的价值在于识别易失效架构并优先筛选力学稳定的候选支架，为后续匹配实验提供计算筛选基础，而非临床验证声明。
研究结论为：在保守与激进两种康复方案下，13种架构保持结构完整性并收敛至接近成熟组织密度，表明力学稳定性维持时预测成骨进展一致。失效为拓扑及方案依赖性，Lidinoid在两种方案下均坍塌，Gyroid与SplitP仅在激进载荷下失效。面心立方、简单立方、Octet及Voronoi属于有利的低损伤组。该框架为识别失效倾向架构并优先筛选力学稳定候选支架提供了计算筛选基础。
