基于区域校准适应参数的下颌骨重建数值建模研究

《Bioengineering》：Numerical Modeling of Mandibular Bone Remodeling with Regional Calibration of the Adaptation Parameter

　　

提出了一种下颌骨重建的力学调控数值模型。该模型在Huiskes-Weinans重建定律中，对位于第一磨牙附近的16个解剖区域（8个松质骨区和8个皮质骨区）的适应参数进行了区域校准，同时引入了密度饱和函数和基于四种咀嚼载荷工况加权的力学刺激。参考材料密度和下颌骨几何结构来源于单一口腔完整患者的锥形束计算机断层扫描数据。采用平均绝对误差、平均绝对百分比误差和平均符号误差评估了23种模型变体与参考密度的一致性。对于选定的变体18（分别对应松质骨和皮质骨），各项指标如下：MAE分别为0.109 g·cm-3和0.200 g·cm-3，MAPE分别约为10.9%和12.7%，Bias分别为+0.006 g·cm-3和+0.038 g·cm-3。与未进行区域校准的模型相比，适应参数的区域校准（变体13）使松质骨的MAE降低了约64%，皮质骨的MAE降低了约56%。在该患者特异性模型中，区域校正是改善一致性的主导因素；由于校准和评估使用的是相同的CBCT数据，因此这种一致性不构成独立的验证。

下颌骨重建数值建模研究解读

研究背景与目的

下颌骨第一磨牙区域承受着显著的咬合载荷，可靠的区域性骨密度分布重建对于预测骨组织对咬合条件改变的响应至关重要。骨组织会根据力学加载条件调整其内部结构，这一现象被形式化为自适应弹性理论和Huiskes-Weinans唯象模型，后者将局部应变能密度与骨密度变化率联系起来。然而，皮质骨与松质骨在重建速率和对力学刺激的敏感性上存在显著差异，使用共享的适应参数是一种简化处理，其对模型定量一致性的影响尚不明确。此外，现有研究中，针对下颌骨的、基于定量影像学验证的重建模型仍然有限，且患者特异性的有限元输入缺乏标准化。因此，研究人员开展此项研究，旨在开发一种包含适应参数区域校准的下颌骨重建模型，并评估其与来自CBCT数据的参考骨密度的一致性。

研究方法与技术路线

研究人员首先获取了一名口腔完整的女性患者的CBCT影像数据，用于确定下颌骨的参考材料密度和几何结构。研究的关键技术方法包括：1）基于HU阈值法在3D Slicer中进行几何重建，并定义围绕第一磨牙的16个计算区域（8个松质骨区CN1-CN8和8个皮质骨区CR1-CR8）；2）应用Huiskes-Weinans重建定律的修正版本，引入密度饱和函数以防止数值发散，并对16个区域的适应参数k进行两阶段区域校准（先确定初始值k₀，再通过缩放因子建立校准值k_new）；3）采用四种咀嚼载荷工况（切牙咬合、双侧磨牙咬合、左侧单侧磨牙咬合、右侧单侧磨牙咬合）的加权组合作为力学刺激；4）使用ANSYS APDL语言实现迭代循环的骨重建算法，并通过网格收敛性测试确保计算精度；5）利用平均绝对误差、均方根误差、平均绝对百分比误差和系统误差四个指标，定量评估23种模型变体的模拟密度与CBCT参考密度的一致性。

研究结果

3.1 区域校准对与CBCT密度一致性的影响

研究表明，引入适应参数k的区域校准显著改善了模型与参考密度的一致性。与使用区域但未校准k的变体6相比，校准后的变体13使松质骨的MAE从0.318 g·cm^{-3-3（降低63.8%），皮质骨的MAE从0.285 g·cm-3降至0.126 g·cm-3（降低55.8%）。改变参数s（惰性区宽度）和B（重建速率系数）并未产生类似的改善效果，表明区域校准k是提高一致性的关键因素。}

3.2 变体与迭代窗口的选择

在所有23种变体中，变体18（校准k_new，B_CN=1000，B_CR=500，s=0.35）在迭代次数121-130的窗口中，同时实现了松质骨和皮质骨的最佳综合误差指标。该变体的松质骨MAE为0.109 g·cm^-3，皮质骨MAE为0.200 g·cm^-3。其他变体无法同时降低两种组织的误差，例如变体19虽获得了最低的松质骨MAE（0.097 g·cm^-3），但皮质骨MAE却升高至0.219 g·cm^-3。

讨论与结论

讨论部分指出，所获MAPE值（松质骨10.86%，皮质骨12.73%）为定量评估提供了依据，但由于先前研究多使用相关系数等不同指标，难以直接比较。研究人员强调，研究中存在两种“收敛”概念：数值收敛（密度变化趋于零）和最佳一致性迭代窗口（模型最接近参考密度）。所有变体的松质骨在200次迭代内均未达到完全数值收敛。研究还明确了密度饱和函数的必要性，其在无饱和的变体中阻止了松质骨密度的非物理性增长。

研究的局限性主要包括：基于单一患者数据，结果不可泛化；CBCT设备缺乏专用体模校准，导致参考密度存在系统不确定性；参考密度同时用于校准k和选择迭代窗口，因此报告的一致性并非独立验证；四种载荷工况的权重为模型假设，未量化其敏感性；未模拟咬合力、牙周膜及颞下颌关节的真实运动；表面网格缩减导致体积减少16-19.7%；重建算法未完全收敛。研究人员认为，该模型在应用于种植体规划前，仍需在独立数据和更多患者上进行验证。

研究结论部分翻译：

开发了一个具有区域校准适应参数k的患者特异性下颌骨重建数值模型，该模型聚焦于第一磨牙周围的16个解剖区域，并与CBCT衍生的骨密度进行了比较。对于选定的变体，与参考密度的一致性表现为：松质骨MAE=0.109 g·cm^-3（MAPE 10.9%），皮质骨MAE=0.200 g·cm^-3（MAPE 12.7%）。与未进行区域校准的模型相比，k的区域校准（变体13，在调整B和s之前）使松质骨的MAE降低了63.8%，皮质骨的MAE降低了55.8%，而随后改变B和s则改善了一种组织的一致性，但牺牲了另一种组织。在进行的单个第一阶段比较中，密度饱和函数限制了在没有它们的情况下观察到的松质骨密度的无限增长。这些结果证明了在所研究的患者特异性模型内一致性的改善；由于校准和评估使用了相同的CBCT数据，这种一致性属于样本内结果。在将该模型应用于种植体规划之前，需要对独立数据和更多患者进行进一步验证。