基于有限元法角膜组织力学变形仿真.pdfVIP

维普资讯 第34卷 第 1期 北 京 工 业 大 学 学 报 V01．34 No．1 2008年 1月 JOURNALOFBEUINGUNIVERSITYOFTECHN0LoGY Jan．2008 基于有限元法角膜组织力学变形仿真 杨艳芳 ，陈定方 ，曾衍钧 ，张 博 (1．武汉理工大学 智能制造与控制研究所，武汉 430063；2．北京工业大学 生物力学和医学信息研究所，北京 100022) 摘 要 ：针对激光角膜热成形术原理，提出了角膜模型的网格划分、加载及约束方案，分析了不同参数情况下角 膜的变形量及其 曲率半径的变化率．通过分析对 比，可以得到不同参量在手术中对角膜 曲率变化效果的影响程 度，从而确定各个参量在屈光矫正中的作用效果．仿真结果表明，眼压的大小对角膜屈光力的变化有一定的影 响，但影响很小；在相同压力的情况下，光凝区直径的变化越大，角膜 曲率的变化就越大，手术效果越明显；光凝 区直径不能太小，小于4．5mm的情况下，角膜曲率变化很小，手术效果不明显． 关键词 ：激光角膜热成形术；角膜；有限元法 中图分类号 ：TP391．9；R778 文献标识码：A 文章编号 ：0254—0037(2008)01—0085—05 用生物力学研究数据建立有临床用途及预测性的生物力学模型，可以极大地提高人们对角膜手术的 理解及对术后的预测．其中，针对放射状角膜切开术(RK)的研究较多．在开展生物力学有限元工作的过 程中，许多学者做了有限元应用理论的基础工作．王家权等对放射状角膜切开术 (RK)有限元法计算机数 值模拟进行了研究[1-2]． 作者针对角膜热成形术的原理及过程，确定了角膜模型的网格划分、加载及约束，实现了不同手术参 数情况下的角膜变形分析 ． 1 虚拟分析的理论依据 根据人眼的屈光学原理[引，人眼屈光力表达式为 D=(722一 1)／y (1) 其中，2为角膜的折射率；1为空气的折射率；y为角膜 曲率半径．因为对于特定的人眼而言，1和 2 为确定的常数，因此角膜的屈光曲率半径 y是影响角膜屈光率的主要参数． 对式(1)进行微分 dD=[(722—721)／72]at 在 y变形的前提下，有 AD=一(n2一 1)a2／ 角膜屈光力 D的变化率为 ： lI×100％=l一 。 I_ ×100％ (2) 由式(2)可知，D的变化率呀与角膜术后的曲率半径变化率有相同的表达式，即手术后角膜 曲率半径 变化率可以准确反映出角膜屈光的变化，表示术后屈光的矫正程度． 因此，可以运用有限元方法，通过建立几何模型，在相应的约束和加载的情况下，分析不同参数下角膜 的变形量及其曲率半径的变化率，通过对比便可以获知每种手术参量在手术中对角膜曲率变化效果的影 响程度，从而确定各个参量在屈光矫正中的作用效果． 收稿 日期：2007—0709． 作者简介：杨艳芳(1975一)，女，河北保定人，讲师 维普资讯 北 京 工 业 大 学 学 报 2008钲 2 手术参数对手术结果的影响分析 由于对角膜屈光手术而言，术后眼球的后半部基本不受影响，因此没有力学研究的意义，只需研究术 后角膜的变形．从生物学的角度看，角膜有不同的生物层并有不同的几何厚度参数，并且每层有不同的材 料常数，如杨氏模量、泊松 比等．由于对角膜各层的材料常数 尚不能确定，只能得到整体角膜材料参数，因 此，在建模的过程中将角膜考虑为一个整体来设计 ．在角膜的建模过程中还需要考虑到4个问题 ： 1)有限元模型中网格划分的密度．网格的划分密度大，分析结果的准确率就高，但分析过程 中资源 耗