光环重力牵引移动系统的不确定性知情结构评估

《Bioengineering》：Uncertainty-Informed Structural Evaluation of a Halo-Gravity Traction Mobility System

　　

光环重力牵引（halo-gravity traction, HGT）作为严重儿童脊柱畸形术前治疗被广泛使用，通过颅固定系统逐步施加牵引力。尽管HGT临床有效性已有大量文献记载，但移动HGT系统的结构响应与运行行为相关不确定性研究相对较少。研究人员提出一种用于移动光环重力牵引轮椅工程评估的初步临床—结构耦合不确定性量化框架。降阶结构与临床替代模型结合拉丁超立方抽样（Latin Hypercube Sampling, LHS）、有界贝塔分布输入变量、三阶多项式混沌展开（Polynomial Chaos Expansion, PCE）和索伯尔（Sobol）敏感性分析。十六个不确定的患者、操作、几何、制造和光环界面参数被传播通过梁—柱组合受载、弹性屈曲、焊接接头受载、基于情景的轮椅倾覆稳定性、光环针载荷集中、牵引传递和代表性临床响应模型。研究人员考察三十二种替代贝塔分布形状案例以评估对边缘分布形状假设的敏感性。在所考察边界和降阶建模假设内，该计算筛查模型未预测到规定结构裕度丧失。结果识别出支配不同响应模式的参数，并展示迁移相关加速度、支撑几何、牵引载荷、制造效率和光环针载荷传递中的不确定性如何在统一框架中评估。由于该模型尚未通过专用样机测试或更高保真度全系统仿真验证，结果应解释为不确定性知情的工程筛查，而非装置安全的实验验证或临床认证。

研究背景方面，有限元分析（finite element analysis, FEA）已成为计算生物力学中骨科植入物和医疗设备设计评估的成熟工具。不确定性量化（uncertainty quantification, UQ）通过考虑材料属性、解剖几何、植入物位置、生理载荷等变异性，提升了计算预测可信度。多项式混沌展开（Polynomial Chaos Expansion, PCE）凭借正交多项式基函数在概率空间中表示随机系统响应，成为高效替代建模方法。然而既有不确定性感知计算方法主要面向可植入骨科器械、心血管系统和患者特异性生物力学分析，对以受控治疗载荷施加为主要目的、而非永久植入的骨科治疗系统关注不足。光环重力牵引（halo-gravity traction, HGT）通过固定于颅骨的光环环施加渐进矫正，用于严重儿童脊柱侧凸和 kyphoscoliosis；临床研究显示其畸形矫正率约25–40%，并改善肺功能和术前状态，但患者间变异性大，脊柱与软组织具时变黏弹性，实际传递载荷因滑轮摩擦、结构柔性、照护者操作、患者运动、轮椅加速度等而偏离处方值。所研究的移动HGT轮椅由阿拉巴马大学伯明翰分校生物医学工程毕业设计原型化，集成轮椅、垂直支撑件、导向配重小车、缆绳滑轮牵引、焊接后支撑框架、防倾覆脚轮支撑和附件硬件；其工程评估须同时考虑牵引传递、梁—柱受载、屈曲抗力、焊缝完整性、情景倾覆稳定性和光环针载荷传递，确定性名义分析不足以表征可能系统响应，因此研究人员建立耦合临床—结构UQ框架。

主要关键技术方法方面，研究人员采用拉丁超立方抽样（Latin Hypercube Sampling, LHS）对16个不确定输入抽样；以有界贝塔（beta）分布描述各变量边际分布，并设32种共同形状参数案例进行分布形状敏感性研究；以三阶PCE（总阶数多项式基，16输入下969项）经普通最小二乘和QR分解拟合标量响应替代模型，由Jacobi多项式构成beta分布对应基函数；由PCE系数直接计算一阶与全阶索伯尔（Sobol）指数做方差全局敏感性分析；输入含患者体重、牵引比、有效脊柱刚度、活跃光环针数、截面面积、弱轴惯性矩、偏心距、焊造效率、后支撑与配重位置、轴距、轮距、垂直加速度、针—颅界面筛查载荷、针载荷放大因子、针载荷共享效率等，弹性模量、柱长和固—自由长度因子作确定性量。

研究结果方面，第2节系统描述与结构建模中，研究人员将垂直C形槽支撑保守理想化为固—自由悬臂梁—柱，牵引相关载荷经有限偏心水平臂传入，因而同时承受轴向压缩、一阶偏心弯曲和二阶几何放大（P–Δ）效应；给出牵引力F_t=r_trW、配重载荷W_s=F_t/η_p、等效压缩P_eq=F_t+a_vW_s、平均压应力σ_c=P_eq/A、欧拉临界载荷P_cr=π²EI_min/(KL)²（K=2）、一阶弯矩M₁=P_eqe、弯矩放大M=M₁/(1?P_eq/P_cr)、极值纤维压应力σ_max=σ_c+Mc/I、欧拉安全系数FS_buck=P_cr/(η_fabP_eq)和组合安全系数FS_comb=S_y/σ_max；结论是这些降阶公式保留主载荷机制且适于UQ，但反应分配仍待实验验证。

第3节轮椅几何与基于情景的倾覆稳定性中，研究人员采用成人医院轮椅代表尺寸，定义前向急停绕前脚轮接触线、侧向单边路缘跌落和10度上坡后向倾覆三种固定条件情景；以恢复矩与倾覆矩计算各情景安全系数，确定情景包络FS_tip,env=min(FS_fwd,FS_lat,FS_rear)。结论为代表几何和质量分布下急停工况主导但FS_fwd仍达4.67，后向上坡FS_rear=5.61不主导；10 m/s一秒匀减速的失控示例FS_fwd低于1.5筛查裕度但高于物理倾覆阈值1，且趋向前向而非后向失稳。该节明确情景稳定性屏检不赋予事故概率。

第4节焊接接头与疲劳筛查中，研究人员以等腿角焊缝有效喉厚a=0.707l_w算喉区面积A_w=0.707l_wL_w、平均喉剪应力τ_w=P_eq/A_w、焊缝安全系数FS_weld=τ_allowη_fab/τ_w和极限状态g_weld=FS_weld?1；另以交替应力σ_a=f_altσ_lat和疲劳比较因子FS_fat=S_e/σ_a作诊断输出但不计入主失效模式。结论为焊缝与疲劳指标属筛查级，缺测量载荷谱和焊缝类别时不能给寿命结论。

第5节耦合临床—结构不确定性模型中，临床输入为W、r_tr、k_sp、N_pins，牵引力由F_t=r_trW确定；临床替代位移δ=F_t/k_sp、归一化响应R_clin=δ/δ₀仅为线性刚度替代而非验证生物力学模型。装置与操作向量含A、I_min、e、η_fab、x_rear、x_stack、L_wb、B、a_v、P_scr、K_amp、η_share；光环针平均载荷F_pin,avg=F_t/N_pins，最大估计F_pin,max=K_ampF_t/(η_shareN_pins)，屏检因子FS_pin=P_scr/F_pin,max。研究人员指出A与I_min独立抽样可能不对应单一C形槽实截面，是筛查近似局限。系统响应取最小安全系数FS_min=min(FS_buck,FS_comb,FS_weld,FS_pin)，屈曲阈值2.0、其余1.0；二进制违例I_fail=1若任一主结构或针界面极限态被破坏。结论为该章给出耦合映射与极限态结构，但所有响应都条件依赖于降阶方程、边界、beta形状和独立假设。

第6节不确定性量化方法中，LHS每案例4000实现、rng(1)可复现；beta变换ξ=(b?a)F_beta^?1(u)+a；32案例同形状参数赋全部16输入，N_pins潜变量映到4–8后取整。PCE以归一化输入总阶3、Jacobi基、OLS+QR拟合，全样本RMSE中位数0.4082，标准化RMSE中位数0.1341，R²中位数高，条件数中位数5.40；80/20留一验证R²中位数仍佳，首阶Sobol向量全样本与重拟合相关0.999607–0.999998，最高敏感变量32案均不变。结论为PCE替代精度总体可接受，但强偏态案精度下降，Sobol指数应作近似替代敏感度读。

讨论与结论部分，研究人员讨论称：聚合响应中组合梁—柱与光环针中心安全系数明显低于屈曲和焊缝；跨案例平均首阶Sobol指数主要由光环针子系统与处方牵引载荷变量贡献，即工程针—颅筛查载荷P_scr、针载荷放大K_amp、活跃针数N_pins、牵引比r_tr、针共享效率η_share影响FS_min屏检；结构和支撑几何变量作用于其他模式但均值贡献较小。情景倾覆不与普通使用失效概率混合。因为模型未由样机载荷试验或全系统高保真仿真验证，所以结果只是不确定性知情工程筛查，不是装置安全实验验证或临床认证。

研究人员结论可译为：在本研究考察的边界与降阶建模假设内，该计算筛查模型未预测到规定结构裕度丧失；结果识别出支配不同响应模式的参数，并展示迁移相关加速度、支撑几何、牵引载荷、制造效率和光环针载荷传递中的不确定性如何在统一框架中评估；由于模型尚未通过专用原型测试或更高保真度全系统仿真验证，结果应解释为不确定性知情的工程筛查，而非设备安全的实验验证或临床认证；该框架用于初步工程评估、识别有影响参数和竞争响应模式，并指导后续模型精炼、样机载荷测试、动态稳定性测试与更高保真仿真，不构成患者特异临床结果预测或人群失效概率确立。