CN117542531A - 一种角膜内皮细胞有限元模型的建立方法 - Google Patents

一种角膜内皮细胞有限元模型的建立方法 Download PDF

Info

Publication number
CN117542531A
CN117542531A CN202311441515.0A CN202311441515A CN117542531A CN 117542531 A CN117542531 A CN 117542531A CN 202311441515 A CN202311441515 A CN 202311441515A CN 117542531 A CN117542531 A CN 117542531A
Authority
CN
China
Prior art keywords
model
corneal endothelial
grid
endothelial cells
finite element
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202311441515.0A
Other languages
English (en)
Inventor
谢立信
李东芳
李宗义
周庆军
林萍
王新航
王超颖
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute Of Ophthalmology Affiliated To Shandong First Medical University Shandong Institute Of Ophthalmology And Qingdao Eye Hospital Affiliated To Shandong First Medical University
Original Assignee
Institute Of Ophthalmology Affiliated To Shandong First Medical University Shandong Institute Of Ophthalmology And Qingdao Eye Hospital Affiliated To Shandong First Medical University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Institute Of Ophthalmology Affiliated To Shandong First Medical University Shandong Institute Of Ophthalmology And Qingdao Eye Hospital Affiliated To Shandong First Medical University filed Critical Institute Of Ophthalmology Affiliated To Shandong First Medical University Shandong Institute Of Ophthalmology And Qingdao Eye Hospital Affiliated To Shandong First Medical University
Priority to CN202311441515.0A priority Critical patent/CN117542531A/zh
Publication of CN117542531A publication Critical patent/CN117542531A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/50ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]

Landscapes

  • Engineering & Computer Science (AREA)
  • Public Health (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Epidemiology (AREA)
  • Pathology (AREA)
  • Databases & Information Systems (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Data Mining & Analysis (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Biomedical Technology (AREA)
  • Prostheses (AREA)

Abstract

本发明提供了一种角膜内皮细胞有限元模型的建立方法,包括:将标准角膜模型导入有限元前处理软件中,并生成2d网格,以得到基础网格模型;导入角膜内皮细胞详细几何模型;根据所述基础网格模型和所述详细几何模型确定完整的角膜内皮细胞模型;根据所述完整的角膜内皮细胞模型进行计算。本发明能够提高建立角膜内皮细胞有限元模型的效率及效果。

Description

一种角膜内皮细胞有限元模型的建立方法
技术领域
本发明涉及模型设计领域,特别是涉及一种角膜内皮细胞有限元模型的建立方法。
背景技术
角膜内皮细胞是角膜的最内层细胞,为单层六边形细胞,它对于保持角膜透明度和正常视觉功能具有重要作用。角膜内皮细胞密度以每年0.3%~0.6%的速度缓慢下降,如遇到创伤、内眼手术等则内皮细胞的损失更多。损伤后的内皮细胞通过扩大移行进行修复,如角膜内皮细胞损失过多,内皮的代谢泵功能失代偿可引起角膜水肿,细胞间紧密连接断裂或者缺失,细胞部分边界缺失、融合,屏障功能降低,临床上即表现为角膜内皮功能失代偿。在角膜疾病的诊断和治疗中,角膜内皮细胞的形态和数量往往是重要的参考指标。
目前,大量研究着重关注角膜内皮细胞密度和六边形细胞占比。然而临床发现存在大量细胞密度高但角膜已明显水肿或细胞密度低但角膜仍处于透明状态的现象,表明当前角膜内皮细胞形态检查的方法不能准确反映细胞功能。
发明内容
为了克服现有技术的不足,本发明的目的是提供一种角膜内皮细胞有限元模型的建立方法。
为实现上述目的,本发明提供了如下方案:
一种角膜内皮细胞有限元模型的建立方法,包括:
将标准角膜模型导入有限元前处理软件中,并生成2d网格,以得到基础网格模型;
导入角膜内皮细胞详细几何模型;
根据所述基础网格模型和所述详细几何模型确定完整的角膜内皮细胞模型;
根据所述完整的角膜内皮细胞模型进行计算。
优选地,所述生成2d网格,包括:
选择角膜内皮表面生成2d网格,最小网格尺寸为2μm,最大网格尺寸为368μm,作为基础网格模型,模型中间网格尺寸为2μm,在径向上逐渐增大;
使用过渡网格划法,保持单元全部为四边形单元;
删除中间的网格,作为基础网格,留出469μmⅹ324μm空间,以替换详细细胞模型。
优选地,在导入角膜内皮细胞详细几何模型之后,还包括:
根据所述详细几何模型生成2d网格;网格尺寸为2μm。
优选地,根据所述基础网格模型和所述详细几何模型确定完整的角膜内皮细胞模型,包括:
使用梁单元模拟细胞之间的连接蛋白;
将详细几何模型生成的网格投影到角膜内皮表面上,使详细几何模型贴合角膜内皮表面;
将基础网格模型和详细网格模型连接,组成完整的角膜内皮细胞模型;
将所有网格投影到角膜内皮表面,完整网格模型是贴合角膜内皮细胞表面。
优选地,根据所述完整的角膜内皮细胞模型进行计算,包括:
在有限元模拟软件中给角膜内皮细胞模型赋予材料及截面属性,材料本构选择Neo-Hooke,材料参数C10设置为173Pa,材料参数D1设置为0.00061/Pa,并为内皮细胞定义厚度为5μm;
细胞之间的连接为梁单元,给梁单元赋予圆形截面;所述圆形截面的半径为1μm;
将角膜内皮细胞最外圈节点约束X、Y和Z方向平动自由度,以固定所述完整的角膜内皮细胞模型;
选择角膜内皮细胞模型,施加压力2E-07MPa,以模拟眼压;
在所述有限元模拟软件中提交计算,得到计算结果。
根据本发明提供的具体实施例,本发明公开了以下技术效果:
本发明提供了一种角膜内皮细胞有限元模型的建立方法,包括:将标准角膜模型导入有限元前处理软件中,并生成2d网格,以得到基础网格模型;导入角膜内皮细胞详细几何模型;根据所述基础网格模型和所述详细几何模型确定完整的角膜内皮细胞模型;根据所述完整的角膜内皮细胞模型进行计算。本发明能够提高建立角膜内皮细胞有限元模型的效率及效果。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的方法流程图;
图2为本发明实施例提供的标准角膜几何模型内表面示意图;
图3为本发明实施例提供的网格尺寸过渡示意图;
图4为本发明实施例提供的基础网格模型示意图;
图5为本发明实施例提供的PKP术后代偿期角膜内皮细胞详细几何模型示意图;
图6为本发明实施例提供的PKP术后失代偿期角膜内皮细胞详细网格模型示意图;
图7为本发明实施例提供的梁单元连接示意图;
图8为本发明实施例提供的角膜内皮细胞完整模型示意图;
图9为本发明实施例提供的基础模型与详细网格模型连接示意图;
图10为本发明实施例提供的角膜约束示意图;
图11为本发明实施例提供的压力示意图;
图12为本发明实施例提供的正常角膜内皮细胞网格模型示意图;
图13为本发明实施例提供的正常人角膜单个内皮细胞第一主应力云图(MPa);
图14为本发明实施例提供的实施例细胞边梁单元示意图;
图15为本发明实施例提供的实施例PKP术后失代偿期角膜内皮细胞网格模型示意图;
图16为本发明实施例提供的实施例PKP术后失代偿期人角膜单个内皮细胞第一主应力云图(MPa);
图17为本发明实施例提供的实施例细胞边梁单元示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
本实施例中,针对适用的病例或者病种范围并不进行区分和限定,即如青光眼等疾病涉及到内皮的都适用于本实施例中所涉及的方法。
图1为本发明实施例提供的方法流程图,如图1所示,本发明提供了一种角膜内皮细胞有限元模型的建立方法,包括:
步骤100:将标准角膜模型导入有限元前处理软件中,并生成2d网格,以得到基础网格模型;
步骤200:导入角膜内皮细胞详细几何模型;
步骤300:根据所述基础网格模型和所述详细几何模型确定完整的角膜内皮细胞模型;
步骤400:根据所述完整的角膜内皮细胞模型进行计算。
具体的,眼角膜包含内皮细胞层、连接蛋白、后弹力层、基质层、前弹力层和上皮细胞层,考虑校核方式、计算精度和计算效率等因素,角膜内皮细胞有限元模型的建立方法有多种,在模型层次和面积上有不同的体现。本实施例中的模型包含角膜内皮细胞层和连接蛋白,模型面积和完整内皮细胞层面积相同;具体实施过程如下:
1、将标准角膜模型导入有限元前处理软件Hypermesh中。见图2所示。
2、选择角膜内皮表面生成2d网格,最小网格尺寸为2μm,最大网格尺寸为368μm,作为基础网格模型,模型中间网格尺寸为2μm,在径向上逐渐增大,使用过渡网格划法,保持单元全部为四边形,见图3,删除中间的网格,作为基础网格,留出469μmⅹ324μm空间,以替换详细细胞模型,见图4。
3、导入穿透性角膜移植术(penetratingkeratoplasty,PKP)后,以下简称PKP,代偿期角膜内皮细胞详细几何模型,见图5;并生成2d网格,网格尺寸为2μm,见图6。
4、细胞之间的连接蛋白使用梁单元模拟,见图7中黑色连接处为梁单元。
5、将详细模型生成的网格投影到角膜内皮表面上,使详细模型贴合角膜内皮表面,如图8所示。
6、将基础网格模型和详细网格模型连接,组成完整的角膜内皮细胞模型,见图7,基础模型与详细网格模型连接见图9。
7、将所有网格投影到角膜内皮表面,完整网格模型是贴合角膜内皮细胞表面。
8、将完整角膜内皮细胞模型导出inp文件。
9、将inp文件导入到仿真软件Abaqus中。
10、在Abaqus中给内皮细胞网格模型赋予材料及截面属性,材料本构选择Neo-Hooke,材料参数C10设置为173Pa,材料参数D1设置为0.00061/Pa,为内皮细胞定义厚度为5μm。
11、细胞之间的连接为梁单元,给梁单元赋予圆形截面,半径为1μm。
12、将角膜内皮细胞最外圈节点约束X、Y和Z方向平动自由度,见图10,以固定角膜内皮细胞完整模型。
13、选择角膜内皮细胞模型,施加压力2E-07MPa,见图11,模拟眼压。
14、提交计算,计算完成后即可查看结果。
本实施例中模型可为正常人角膜内皮细胞模型,将正常人角膜内皮细胞详细几何模型导入Hypermesh中,之后的步骤按照本方案中的步骤进行,模型见图12,提交计算,查看结果。细胞顶角角度不同,顶角的第一主应力不同(见图13),提取其中一个角的角度及第一主应力,结果为120.212°和163.98MPa。提取细胞边上的梁单元的轴力求和(F),并除以边长(L),作为细胞边上单位长度的力,体现边长和受力的情况,见图14,提取结果,计算F/L为7.14E-10N/μm。提取细胞边长和边上平均应力值,体现细胞边上应力的不同,见图13,提取其中一条边上第一主应力的平均值,结果为158.330Pa。
本实施例中模型还可为PKP术后失代偿期人角膜内皮细胞模型,将PKP术后失代偿期人角膜内皮细胞详细几何模型导入Hypermesh中,之后的步骤按照本方案中的步骤进行,模型见图15,提交计算,查看结果。细胞顶角角度不同,顶角的第一主应力不同,见图16,提取其中一个角的角度及第一主应力,结果为131.624°和164.40MPa。提取细胞边上的梁单元的轴力求和(F),并除以边长(L),作为细胞边上单位长度的力,体现边长和受力的情况,见图17,提取结果,计算F/L为7.96E-10N/μm。提取细胞边长和边上应力值,体现细胞边上应力的不同,见图16,提取其中一条边上第一主应力的平均值,结果为162.289Pa。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。

Claims (5)

1.一种角膜内皮细胞有限元模型的建立方法,其特征在于,包括:
将标准角膜模型导入有限元前处理软件中,并生成2d网格,以得到基础网格模型;
导入角膜内皮细胞详细几何模型;
根据所述基础网格模型和所述详细几何模型确定完整的角膜内皮细胞模型;
根据所述完整的角膜内皮细胞模型进行计算。
2.根据权利要求1所述的角膜内皮细胞有限元模型的建立方法,其特征在于,所述生成2d网格,包括:
选择角膜内皮表面生成2d网格,最小网格尺寸为2μm,最大网格尺寸为368μm,作为基础网格模型,模型中间网格尺寸为2μm,在径向上逐渐增大;
使用过渡网格划法,保持单元全部为四边形单元;
删除中间的网格,作为基础网格,留出469μmⅹ324μm空间,以替换详细细胞模型。
3.根据权利要求1所述的角膜内皮细胞有限元模型的建立方法,其特征在于,在导入角膜内皮细胞详细几何模型之后,还包括:
根据所述详细几何模型生成2d网格;网格尺寸为2μm。
4.根据权利要求1所述的角膜内皮细胞有限元模型的建立方法,其特征在于,根据所述基础网格模型和所述详细几何模型确定完整的角膜内皮细胞模型,包括:
使用梁单元模拟细胞之间的连接蛋白;
将详细几何模型生成的网格投影到角膜内皮表面上,使详细几何模型贴合角膜内皮表面;
将基础网格模型和详细网格模型连接,组成完整的角膜内皮细胞模型;
将所有网格投影到角膜内皮表面,完整网格模型是贴合角膜内皮细胞表面。
5.根据权利要求1所述的角膜内皮细胞有限元模型的建立方法,其特征在于,根据所述完整的角膜内皮细胞模型进行计算,包括:
在有限元模拟软件中给角膜内皮细胞模型赋予材料及截面属性,材料本构选择Neo-Hooke,材料参数C10设置为173Pa,材料参数D1设置为0.00061/Pa,并为内皮细胞定义厚度为5μm;
细胞之间的连接为梁单元,给梁单元赋予圆形截面;所述圆形截面的半径为1μm;
将角膜内皮细胞最外圈节点约束X、Y和Z方向平动自由度,以固定所述完整的角膜内皮细胞模型;
选择角膜内皮细胞模型,施加压力2E-07MPa,以模拟眼压;
在所述有限元模拟软件中提交计算,得到计算结果。
CN202311441515.0A 2023-11-01 2023-11-01 一种角膜内皮细胞有限元模型的建立方法 Pending CN117542531A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202311441515.0A CN117542531A (zh) 2023-11-01 2023-11-01 一种角膜内皮细胞有限元模型的建立方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202311441515.0A CN117542531A (zh) 2023-11-01 2023-11-01 一种角膜内皮细胞有限元模型的建立方法

Publications (1)

Publication Number Publication Date
CN117542531A true CN117542531A (zh) 2024-02-09

Family

ID=89785214

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202311441515.0A Pending CN117542531A (zh) 2023-11-01 2023-11-01 一种角膜内皮细胞有限元模型的建立方法

Country Status (1)

Country Link
CN (1) CN117542531A (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09313442A (ja) * 1996-06-03 1997-12-09 Topcon Corp 角膜内皮細胞計測装置
US20090187386A1 (en) * 2008-01-18 2009-07-23 Bille Josef F Finite element modeling of the cornea
US20180125582A1 (en) * 2016-11-08 2018-05-10 Argusys, Inc. Corneal surgery risk evaluation method and system thereof
CN113171172A (zh) * 2021-04-21 2021-07-27 温州医科大学附属眼视光医院 一种角膜术后情况的模拟方法
CN114927222A (zh) * 2022-06-20 2022-08-19 南昌航空大学 一种个性化人眼模型的构建方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09313442A (ja) * 1996-06-03 1997-12-09 Topcon Corp 角膜内皮細胞計測装置
US20090187386A1 (en) * 2008-01-18 2009-07-23 Bille Josef F Finite element modeling of the cornea
US20180125582A1 (en) * 2016-11-08 2018-05-10 Argusys, Inc. Corneal surgery risk evaluation method and system thereof
CN113171172A (zh) * 2021-04-21 2021-07-27 温州医科大学附属眼视光医院 一种角膜术后情况的模拟方法
CN114927222A (zh) * 2022-06-20 2022-08-19 南昌航空大学 一种个性化人眼模型的构建方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
SALMAN N. KHAN等: "Finite Element Analysis of Descemet’s Stripping Automated Endothelial Keratoplasty (DSAEK) Surgery Allograft to Predict Endothelial Cell Loss", 《CURRENT EYE RESEARCH》, 1 June 2016 (2016-06-01), pages 3 - 9 *
孙立科: "角巩膜本构方程构建及 SMILE 手术有限元分析", 《中国优秀硕士学位论文全文数据库基础科学辑》, no. 3, 15 March 2022 (2022-03-15), pages 002 - 504 *

Similar Documents

Publication Publication Date Title
CN106568660B (zh) 一种复合材料胶接修补结构的剩余疲劳寿命的预测方法
CN107230255B (zh) 一种基于通用模板的牙齿修复方法
CN110889166B (zh) 基于热力耦合约束的航空用轴承支架轻量化设计方法
CN112632869B (zh) 一种基于网格框架的非结构附面层网格生成方法
Pirzadeh Progress toward a user-oriented unstructured viscous grid generator
CN106441839A (zh) 一种超强钢激光成形起落架的强度验证方法
CN111779638A (zh) 一种基于疲劳应力应变的风电叶片缺陷危险等级评估方法
CN117542531A (zh) 一种角膜内皮细胞有限元模型的建立方法
CN115238419A (zh) 一种基于复合材料的轴流压气机叶片设计方法
CN105243686B (zh) 支架模拟方法
CN117438093A (zh) 一种特殊的角膜内皮细胞有限元模型的建立方法
Martineau et al. Anisotropic hybrid mesh generation for industrial RANS applications
CN105653773B (zh) 一种适于3d打印的棱透镜复合眼镜片stl模型的修复方法
CN116384196A (zh) 自膨血管支架有限元分析方法
CN112818581B (zh) 风力发电机组主机架的拓扑优化设计方法
CN113591337A (zh) 一种环形薄壁支撑结构优化设计方法
CN103336907A (zh) 一种基于dsi技术的水库静态库容快速计算方法
WO2022188372A1 (zh) 风电叶片叶根铺层设计方法及风电叶片叶根结构
CN115795961B (zh) 一种沿网格边界裂纹扩展有限元方法
CN111368415B (zh) 一种c形大梁桨叶填充物快速建模方法及系统
CN104504762B (zh) 一种基于拓扑重建的义齿表面局部变形方法
CN113670685A (zh) 一种用于刻画叶身与缘板过渡处结构细节的疲劳试验件
CN106326110A (zh) 一种系统版本开发过程中bug缺陷的修复方法及系统
CN112651058B (zh) 一种基于AutoCAD的轮胎半部件自动反包设计方法
CN117372665B (zh) 一种增材制造成形过程中的缺陷模型数据在线修复方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination