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

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

Info

Publication number
CN117438093A
CN117438093A CN202311441518.4A CN202311441518A CN117438093A CN 117438093 A CN117438093 A CN 117438093A CN 202311441518 A CN202311441518 A CN 202311441518A CN 117438093 A CN117438093 A CN 117438093A
Authority
CN
China
Prior art keywords
model
cornea
finite element
grid
endothelial cell
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
CN202311441518.4A
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 CN202311441518.4A priority Critical patent/CN117438093A/zh
Publication of CN117438093A publication Critical patent/CN117438093A/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。
优选地,根据所述基础网格模型和所述详细几何模型确定完整的角膜内皮细胞模型,包括:
细胞之间的连接蛋白使用梁单元模拟;
将详细几何模型生成的网格投影到角膜内皮表面上,使详细几何模型贴合角膜内皮表面;
将基础网格模型和详细网格模型连接,组成完整的角膜内皮细胞模型。
优选地,根据所述完整的角膜内皮细胞模型生成角膜其他位置的网格,得到角膜完整的有限元模型,包括:
根据角膜几何模型生成角膜其他位置的网格,圆心处的网格细化,网格尺寸为8μm,厚度方向上尺寸为16μm,在径向上粗化网格,粗化部分网格尺寸从17μm逐渐增加到378μm,以得到所述膜完整的有限元模型;所述角膜其他位置包括后弹力层、基质层、前弹力层和上皮细胞层。
优选地,根据所述角膜完整的有限元模型进行计算,包括:
在有限元模拟软件中给角膜完整的有限元模型赋予材料及截面属性,材料本构选择Neo-Hooke,材料参数C10设置为173Pa,材料参数D1设置为0.00061/Pa,为内皮细胞定义厚度为5μm;
细胞之间的连接为梁单元,给梁单元赋予圆形截面;所述圆形截面的半径为1μm;
为所述角膜其他位置材料本构选择线弹性本构,材料参数弹性模量设置为0.1MPa,泊松比设置为0.45;
将角膜端面节点约束X、Y、Z方向平动自由度,以固定所述角膜完整的有限元模型;
对角膜界面节点施加对称约束;
选择角膜内皮细胞模型,施加压力2E-03MPa,以模拟眼压;
在不包括内皮细胞的角膜细化区域和粗化区域之间、内皮细胞和角膜其他位置之间,建立绑定接触;
在所述有限元模拟软件中提交计算,得到计算结果。
根据本发明提供的具体实施例,本发明公开了以下技术效果:
本发明提供了一种特殊的角膜内皮细胞有限元模型的建立方法,包括:将四分之一标准角膜模型导入有限元前处理软件中,并生成2d网格,以得到基础网格模型;导入角膜内皮细胞详细几何模型;根据所述基础网格模型和所述详细几何模型确定完整的角膜内皮细胞模型;根据所述完整的角膜内皮细胞模型生成角膜其他位置的网格,得到角膜完整的有限元模型;根据所述角膜完整的有限元模型进行计算。本发明能够提高建立角膜内皮细胞有限元模型的效率及效果。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的方法流程图;
图2为本发明实施例提供的方案3四分之一标准角膜几何模型示意图;
图3为本发明实施例提供的方案3内皮细胞基础网格模型示意图;
图4为本发明实施例提供的PKP术后代偿期角膜内皮细胞详细几何模型示意图;
图5为本发明实施例提供的PKP术后失代偿期角膜内皮细胞详细网格模型示意图
图6为本发明实施例提供的梁单元连接示意图;
图7为本发明实施例提供的内皮细胞完整网格模型示意图;
图8为本发明实施例提供的内皮细胞基础网格和详细建模网格模型连接示意图;
图9为本发明实施例提供的整体网格模型示意图;
图10为本发明实施例提供的模型端面约束示意图;
图11为本发明实施例提供的模型截面约束示意图;
图12为本发明实施例提供的压力示意图;
图13为本发明实施例提供的正常角膜内皮细胞网格模型示意图;
图14为本发明实施例提供的正常人角膜单个内皮细胞第一主应力云图(MPa);
图15为本发明实施例提供的细胞边梁单元示意图;
图16为本发明实施例提供的PKP术后失代偿期角膜内皮细胞网格模型示意图;
图17为本发明实施例提供的PKP术后失代偿期人角膜单个内皮细胞第一主应力云图(MPa);
图18为本发明实施例提供的细胞边梁单元示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
本实施例中,针对适用的病例或者病种范围并不进行区分和限定,即如青光眼等疾病涉及到内皮的都适用于本实施例中所涉及的方法。
图1为本发明实施例提供的方法流程图,如图1所示,本发明提供了一种特殊的角膜内皮细胞有限元模型的建立方法,包括:
步骤100:将四分之一标准角膜模型导入有限元前处理软件中,并生成2d网格,以得到基础网格模型;
步骤200:导入角膜内皮细胞详细几何模型;
步骤300:根据所述基础网格模型和所述详细几何模型确定完整的角膜内皮细胞模型;
步骤400:根据所述完整的角膜内皮细胞模型生成角膜其他位置的网格,得到角膜完整的有限元模型;
步骤500:根据所述角膜完整的有限元模型进行计算。
具体的,眼角膜包含内皮细胞层、连接蛋白、后弹力层、基质层、前弹力层和上皮细胞层,考虑校核方式、计算精度和计算效率等因素,角膜内皮细胞有限元模型的建立方法有多种,在模型层次和面积上有不同的体现。本实施例中的模型包含角膜内皮细胞层、连接蛋白、后弹力层、基质层、前弹力层和上皮细胞层,面积为完整内皮细胞面积的四分之一。除此之外的建模方法还有:模型包含内皮细胞层、连接蛋白和后弹力层,面积与完整内皮细胞层相同等方法。具体实施过程如下:
1、将四分之一标准角膜模型导入有限元前处理软件Hypermesh中。见图2。
2、选择角膜内皮表面生成2d网格,最小网格尺寸为2μm,最大网格尺寸为368μm,作为基础网格模型,模型中间网格尺寸为2μm,在径向上逐渐增大,使用过渡网格划法,保持单元全部为四边形单元,见图3,删除中间的网格,作为基础网格,留出469μmⅹ324μm空间,以替换详细细胞模型。
3、导入角膜内皮细胞详细几何模型,见图4;并生成2d网格,网格尺寸为2μm,见图5。
4、细胞之间的连接蛋白使用梁单元模拟,见图6中黑色连接处为梁单元。
5、将详细模型生成的网格投影到角膜内皮表面上,使详细模型贴合角膜内皮表面。
6、将基础网格模型和详细网格模型连接,组成完整的角膜内皮细胞模型,见图7,基础模型与详细网格模型连接见图8。
7、根据角膜几何模型生成角膜其他位置(后弹力层、基质层、前弹力层和上皮细胞层)网格,圆心处的网格细化,网格尺寸为8μm,厚度方向上尺寸为16μm,在径向上粗化网格,粗化部分网格尺寸从17μm逐渐增加到378μm,见图9。
8、将方案3完整模型导出inp格式文件。
9、将inp文件导入到仿真软件Abaqus中。
10、在Abaqus中给内皮细胞网格模型赋予材料及截面属性,材料本构选择Neo-Hooke,材料参数C10设置为173Pa,材料参数D1设置为0.00061/Pa,为内皮细胞定义厚度为5μm。
11、细胞之间的连接为梁单元,给梁单元赋予圆形截面,半径为1μm。
12、为后弹力层、基质层、前弹力层和上皮细胞层材料本构选择线弹性本构,材料参数弹性模量设置为0.1MPa,泊松比设置为0.45。
13、将角膜端面节点约束X、Y、Z方向平动自由度,见图10。
14、对角膜界面节点施加对称约束,见图11。
15、选择角膜内皮细胞模型,施加压力2E-03MPa,见图12,模拟眼压。
16、在角膜(不包括内皮细胞)细化区域和粗化区域之间、内皮细胞和其它(后弹力层、基质层、前弹力层和上皮细胞层)之间,建立绑定接触。
提交计算,计算完成后即可查看结果。
本实施例可以为正常人角膜内皮细胞模型,将正常人角膜内皮细胞详细几何模型导入Hypermesh中,之后的步骤按照本实施例中的步骤进行,模型见图13,提交计算,查看结果。细胞顶角角度不同,顶角的第一主应力不同(见图14),提取其中一个角的角度及第一主应力,结果为126.585°和143.64MPa。提取细胞边上的梁单元的轴力求和(F),并除以边长(L),作为细胞边上单位长度的力,体现边长和受力的情况,见图15,提取结果,计算F/L为1.40E-10N/μm。提取细胞边长和边上平均应力值,体现细胞边上应力的不同,见图14,提取其中一条边上第一主应力的平均值,结果为143.558Pa。
实施例6为PKP术后失代偿期人角膜内皮细胞模型,将PKP术后失代偿期人角膜内皮细胞详细几何模型导入Hypermesh中,之后的步骤按照本实施例中的步骤进行,模型见图16,提交计算,查看结果。细胞顶角角度不同,顶角的第一主应力不同,见图17,提取其中一个角的角度及第一主应力,结果为136.415°和143.43MPa。提取细胞边上的梁单元的轴力求和(F),并除以边长(L),作为细胞边上单位长度的力,体现边长和受力的情况,见图18,提取结果,计算F/L为1.40E-10N/μm。提取细胞边长和边上应力值,体现细胞边上应力的不同,见图17,提取其中一条边上第一主应力的平均值,结果为143.369Pa。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。

Claims (6)

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所述的特殊的角膜内皮细胞有限元模型的建立方法,其特征在于,根据所述完整的角膜内皮细胞模型生成角膜其他位置的网格,得到角膜完整的有限元模型,包括:
根据角膜几何模型生成角膜其他位置的网格,圆心处的网格细化,网格尺寸为8μm,厚度方向上尺寸为16μm,在径向上粗化网格,粗化部分网格尺寸从17μm逐渐增加到378μm,以得到所述膜完整的有限元模型;所述角膜其他位置包括后弹力层、基质层、前弹力层和上皮细胞层。
6.根据权利要求1所述的特殊的角膜内皮细胞有限元模型的建立方法,其特征在于,根据所述角膜完整的有限元模型进行计算,包括:
在有限元模拟软件中给角膜完整的有限元模型赋予材料及截面属性,材料本构选择Neo-Hooke,材料参数C10设置为173Pa,材料参数D1设置为0.00061/Pa,为内皮细胞定义厚度为5μm;
细胞之间的连接为梁单元,给梁单元赋予圆形截面;所述圆形截面的半径为1μm;
为所述角膜其他位置材料本构选择线弹性本构,材料参数弹性模量设置为0.1MPa,泊松比设置为0.45;
将角膜端面节点约束X、Y、Z方向平动自由度,以固定所述角膜完整的有限元模型;
对角膜界面节点施加对称约束;
选择角膜内皮细胞模型,施加压力2E-03MPa,以模拟眼压;
在不包括内皮细胞的角膜细化区域和粗化区域之间、内皮细胞和角膜其他位置之间,建立绑定接触;
在所述有限元模拟软件中提交计算,得到计算结果。
CN202311441518.4A 2023-11-01 2023-11-01 一种特殊的角膜内皮细胞有限元模型的建立方法 Pending CN117438093A (zh)

Priority Applications (1)

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

Applications Claiming Priority (1)

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

Publications (1)

Publication Number Publication Date
CN117438093A true CN117438093A (zh) 2024-01-23

Family

ID=89551207

Family Applications (1)

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

Country Status (1)

Country Link
CN (1) CN117438093A (zh)

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
AHMED ELSHEIKH: "Finite Element Modeling of Corneal Biomechanical Behavior", 《JOURNAL OF REFRACTIVE SURGERY》, vol. 26, no. 4, 31 December 2010 (2010-12-31), pages 289 - 300 *
HUA LI等: "Biomechanical Effects of Deep Anterior Lamellar Keratoplasty and Penetrating Keratoplasty for Keratoconus: A Finite Element Analysis", 《TRANSLATIONAL VISION SCIENCE & TECHNOLOGY》, vol. 10, no. 9, 13 August 2021 (2021-08-13), pages 1 - 15 *
刘兵等: "基于有限元的角膜建模及RK手术分析", 《中国数字医学》, vol. 3, no. 12, 31 December 2008 (2008-12-31), pages 14 - 16 *
唐聪恺: "基于有限元的三维角膜构建及其生物力学分析", 《中国知网》, 31 December 2006 (2006-12-31), pages 1 - 68 *

Similar Documents

Publication Publication Date Title
CN106568660B (zh) 一种复合材料胶接修补结构的剩余疲劳寿命的预测方法
CN112632869B (zh) 一种基于网格框架的非结构附面层网格生成方法
CN107169191A (zh) 一种风机叶片建模方法
CN110889166B (zh) 基于热力耦合约束的航空用轴承支架轻量化设计方法
CN105930616A (zh) 一种飞机动翼面超静定结构的载荷分配方法
CN111779638A (zh) 一种基于疲劳应力应变的风电叶片缺陷危险等级评估方法
CN104318030B (zh) 基于变形等效的拱坝坝基综合变形模量自动化计算方法
CN108090292A (zh) 一种宽弦风扇叶片二维有限元建模方法
CN115238419A (zh) 一种基于复合材料的轴流压气机叶片设计方法
CN104484527B (zh) 一种离散结构拓扑优化过程中均布载荷自动动态修改方法
CN117438093A (zh) 一种特殊的角膜内皮细胞有限元模型的建立方法
Raut et al. An approach for patient-specific multi-domain vascular mesh generation featuring spatially varying wall thickness modeling
CN106372274A (zh) 用于确定飞续飞载荷谱的低载截除限值的方法
CN115392094A (zh) 一种基于热力耦合的涡轮盘结构优化方法
CN117542531A (zh) 一种角膜内皮细胞有限元模型的建立方法
CN105243686B (zh) 支架模拟方法
CN112949004B (zh) 风力发电机组轴承座的轻量化设计方法及其轴承座
CN117077294B (zh) 基于参数化的结构网格生成方法、装置、设备及介质
CN116305589B (zh) 一种直升机桨叶结构降阶分析方法、系统、设备及介质
CN108197353A (zh) 一种基于ansys的apdl语言的固体火箭发动机工装设计方法
CN103336907A (zh) 一种基于dsi技术的水库静态库容快速计算方法
CN112818581B (zh) 风力发电机组主机架的拓扑优化设计方法
CN116384196A (zh) 自膨血管支架有限元分析方法
CN110287632B (zh) 固体氧化物燃料电池阴极接触电阻数值模拟方法
CN106600671A (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