角膜可视化Scheimpflug技术（Corvis ST, CVS）生物力学参数对角膜材料属性区域依赖性的有限元研究

《Bioengineering》：Regional Dependence of Corvis ST Biomechanical Parameters on Corneal Material Properties: A Finite Element Study

　　

背景：角膜可视化Scheimpflug技术（Corvis ST, CVS）可提供用于生物力学评估的动态角膜响应（dynamic corneal response, DCR）参数，但传统指标可能无法唯一反映角膜材料属性的区域特异性改变。研究人员通过有限元研究考察角膜材料属性区域变化如何影响CVS生物力学参数，并探索用于区域力学解读的多参数策略。方法：构建理想化有限元角膜模型，划分为中央区、旁中央区与周边区，各分区刚度独立调整为基线的70–130%以模拟区域软化或硬化；对前表面施加气冲载荷，利用节点坐标计算常规CVS参数及新提出的DA Ratio X系列；采用归一化相对变化与单变量线性回归刻画响应方向与敏感性。结果：多数参数对区域刚度改变近似线性响应，响应方向与幅度存在区域依赖差异，中央区材料变化引发最大响应；不同区域力学改变可产生相似参数值或趋势，表明单参数无法唯一识别潜在区域改变。在该理想化框架内，扩展DA Ratio系列改善了区域区分能力，DeflArea、DAR4与IIR的组合为区分中央、旁中央与周边改变提供了候选理论框架。结论：CVS生物力学参数受角膜材料改变空间位置影响，单参数解读可能无法唯一反映潜在区域力学变化；多参数协同解读在不同模拟区呈现特异响应模式，可为未来区域生物力学评估提供思路，需结合患者特异性模型与临床数据进一步验证。

研究背景与立项依据

作为人眼主要屈光结构，角膜形态与生物力学属性直接决定视觉质量，且其生物力学改变常早于形态异常，是圆锥角膜（keratoconus, KC）筛查与激光屈光手术（laser vision correction, LVC）安全评估的核心。角膜可视化Scheimpflug技术（Corvis ST, CVS）基于超高速Scheimpflug成像提取动态角膜响应（DCR）参数，较 ocular response analyzer（ORA）提升了在体评估维度，但常规参数基于均质材料假设，而原生角膜呈非线性、各向异性且具明显区域差异。临床中存在矛盾现象：LVC术后DAR2下降（提示变硬）而IIR同时上升（提示变软）；圆锥角膜植入 intracorneal ring segment（ICRS）后标准CVS参数无显著变化却确有力学增强。这引出两个未解问题：CVS参数能否可靠反映不同区域力学改变？同一参数变化是否唯一对应特定力学改变？为此，研究人员采用有限元分析（finite element analysis, FEA）系统考察区域材料属性改变对CVS参数的影响，论文发表于《Bioengineering》。

主要关键技术方法

研究人员基于MATLAB调用Abaqus构建理想化角膜有限元模型，前表面为非球面（中心半径7.8 mm，形状因子0.82），厚度由中央0.55 mm至角巩膜缘0.89 mm，15节点二次三棱柱单元（C3D15H）离散为24环×2层，层厚比2:1表征前基质主导贡献；角巩膜缘滚支约束，后表面眼内压（intraocular pressure, IOP）恒为15 mmHg。材料采用一阶Ogden超弹性近不可压模型（μ₁=0.0485，α₁=108.57），按半径0–1.5 mm、1.5–3 mm、3–5.8 mm划分为中央、旁中央、周边三区，各区μ₁独立取基线70–130%（5%步进）。气冲载荷时空分布引自既往CFD研究（峰值98 mmHg）。由节点坐标计算8个常规CVS参数及DA Ratio X（X=1,2,3,4,5 mm）系列；以100例健康屈光手术候选者（温州医科大学眼视光医院伦理审批2023-032-K-25）在体CVS测量作临床合理性对照（Z-score、相对误差、95%参考区间）。响应量化采用归一化相对变化Y_rel与单变量线性回归Y_rel=a+bX，按|b|分低/中/高/极高敏感。

研究结果

3.1 模拟结果临床合理性评估

基线模型八参数全部落入100眼健康人95%参考区间，相对误差0.49%–13.84%（六者<10%），表明群体代表性几何、载荷与材料输入与临床参考范围一致，但属合理性而非准确性验证。

3.2 CVS参数对区域材料变化的响应趋势

各参数在70–130%刚度内单调变化，R²近1（周边IIR除外R²=0.440）；中央区多为高至极高敏感（|b|≥0.5），旁中央降为中敏感（0.1≤|b|<0.5），周边为低至中敏感。同参数跨区域方向可反转：DeflArea随刚度增加各区均减；DAR2中央变硬时减、旁中央/周边变硬时增。案例显示中央刚度+5%与旁中央?30%得相近IIR（7.65 vs 7.64），旁中央?30%复制中央+25%的DAR2（3.54 vs 3.52），证实单参数非唯一性。

3.3 多参数组合映射区域材料变化

引入DAR3/4/5后，DAR4、DAR5在旁中央变硬时趋势与DAR1–3反转，形成第三响应模式。按斜率符号归为三型：A型（SP-A1、DeflArea、DeflA Max、HCdArcLength）全区同向；B型（DAR4、DAR5）可辨周边；C型（IIR、IR、DAR1–3）中央与旁中央异向。选DeflArea（A型）+DAR4（B型）+IIR（C型）组合理想框架下可区分：三者同升示中央软化；DeflArea与DAR4升、IIR降示旁中央弱化；DeflArea升、DAR4降示周边软化。

讨论部分总结

研究人员指出IIR、DAR2随中央软化上升、随周边软化下降，印证LVC术后趋势但揭示空间异位可致同向参数反向解读；中央软化与旁中央硬化可产生不可区分的IIR/DAR2趋势，解释亚临床圆锥角膜（始于旁中央弱化）常规CVS漏诊。定量上周边需约30%刚度下降方达中央<10%下降的同等参数偏移，中央敏感性梯度最高。DA Ratio X系列中DAR4/5对旁中央呈负响应、DAR1–3呈正响应，构成旁中央指纹。DeflArea+DAR4+IIR组合在模拟中给出区域映射规则，但尚缺患者模型与队列验证。局限包括理想同心区、均质分区、固定IOP/厚度、各向同性Ogden未含胶原纤维异性与水合依赖、未做网格收敛、未临床复现。

结论部分翻译

在该理想化有限元研究中，CVS参数响应取决于所设区域材料改变的位置，相似单参数值可源于不同模拟条件。DeflArea、DAR4与IIR的组合在模拟中央、旁中央与周边情景间呈现特异响应模式，提示其对未来区域生物力学评估的潜在价值。需引入患者特异性几何、实验测量与临床数据集进一步评估该框架在生理情境中的适用性。