CN105373658A - 一种基于人颅颌面骨的三维有限元模型的建模方法 - Google Patents

一种基于人颅颌面骨的三维有限元模型的建模方法 Download PDF

Info

Publication number
CN105373658A
CN105373658A CN201510751461.7A CN201510751461A CN105373658A CN 105373658 A CN105373658 A CN 105373658A CN 201510751461 A CN201510751461 A CN 201510751461A CN 105373658 A CN105373658 A CN 105373658A
Authority
CN
China
Prior art keywords
model
finite element
dimensional
element model
dimensional finite
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
CN201510751461.7A
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.)
Third Military Medical University TMMU
Second Affiliated Hospital of TMMU
Original Assignee
Second Affiliated Hospital of TMMU
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 Second Affiliated Hospital of TMMU filed Critical Second Affiliated Hospital of TMMU
Priority to CN201510751461.7A priority Critical patent/CN105373658A/zh
Publication of CN105373658A publication Critical patent/CN105373658A/zh
Pending legal-status Critical Current

Links

Classifications

    • 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]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Graphics (AREA)
  • Software Systems (AREA)
  • Surgical Instruments (AREA)

Abstract

本发明涉及医学研究建模领域,具体为一种基于人颅颌面骨的三维有限元模型的建模方法,把人头颅分成三部分,对每一个部分分别生成三维实体模型,再重新组合成一个实体整体,最终通过该实体整体生成人头颅的三维有限元模型;其中,三部分具体如下:下颌骨分离形成单独的面网格模型,再将剩余部分分为颅顶和面中部两块。本发明的有益效果如下:1、创建了头面部撞击伤三维有限元模型,能模拟出接近真实情况的骨折效果;2、为头面部撞击伤的研究提供完善的有限元分析工具;3、为复杂三维有限元模型的建模提供一种解决方案,即分割再组合的建模方法。

Description

一种基于人颅颌面骨的三维有限元模型的建模方法
技术领域
本发明涉及医学研究建模领域,具体为一种基于人颅颌面骨的三维有限元模型的建模方法。
背景技术
颌面部是人体暴露的部位,无论是战时还是平时,都是易遭受外力打击的部位,即使人体防护装备不断更新改进,但颌面部的防护仍然是比较薄弱的部位。和平时期,颌面部遭受外力打击最常见于交通事故伤。而颌面部紧邻颅脑,在颌面部遭受撞击时,同时伴发的颅脑损伤是导致伤员死亡的最主要原因之一,因此研究颌面部撞击对颅脑造成的损伤有利于我们提高颌面部撞击伤的救治成功率,能为我们研究更加有效的防护装置提供更多的资料和数据。
有限元分析法作为生物力学研究的最重要的方法之一,目前已被广泛应用于人体生物力学的研究。其核心思想是把整个物体结构看作是由有限个细小单元相互连接而成的几何实体,各单元力学特征的总装效果反映结构的整体力学特性。能够分析物体内及物体间的复杂力学过程,预测力学作用产生的效应(如模型的应力、应变、形状、温度等的变化),并且能够在电脑上直观显示或输出计算结果以供分析,具有可重复性好、节约实验成本、实验条件易控制等优点,可以弥补传统的撞击伤模型在致伤过程研究中的不足。
现有技术一的技术方案
本发明前期完成人下颌骨枪弹伤三维有限元模型,建模方法:CT扫描人下颌骨,获得dicom数据,导入mimics软件建立下颌骨三维数字模型,并进行面网格划分,然后导入ANSA软件中进行体网格划分,手动修改畸形网格,由封闭表面构建实体,生成三维实体有限元模型。
现有技术一的缺点
1、只涉及下颌骨,形态相对较规则,结构简单;
2、用于枪弹伤的有限元力学分析,属于高速撞击的范畴,不适用低速撞击的仿真模拟;
现有技术二的技术方案
CT扫描颅骨,并使用VworksTM4.0Surgery软件进行分割,重建出颅骨三维面网格模型,然后以VRML格式导入到ANSYSICEMCFD12.0.1(ANSYSInc.,Canonsburg,PA,USA)软件中生成体网格模型。
现有技术二的缺点
1、只是用于静态的力学传导显示及分析,无法模拟出撞击后的骨折效果;
2、仍然不是完整的颅面部模型。
本发明拟用有限元方法进行颌面部撞击伤对颅底损伤的动态仿真模拟。首先建立人头颅的三维有限元模型。采用CT扫描人头颅,获取原始数据并用三维重建软件重建出人头颅及颌面部骨组织三维模型,以此三维模型生成三维有限元模型。因人颌面部及颅骨结构复杂,形态极不规则,尝试用各种有限元前处理软件都无法按照常规方法直接生成生成人颅面三维有限元模型,主要是很多颅脑内部结构存在畸形的网格,自相交面等无法手动修复,生成的三维有限元模型无法用于有限元计算。因此,如何把人颅颌面复杂而不规则的结构构建成能用于有限元计算和分析的三维有限元模型成为必须解决的关键问题。
发明内容
本发明的目的在于提供一种基于人颅颌面骨的三维有限元模型的建模方法,构建可用于撞击伤的有限元计算分析的人头颅三维有限元模型。
为了解决以上问题,本发明采用的技术方案如下,一种基于人颅颌面骨的三维有限元模型的建模方法,包括以下步骤:把人头颅分成三部分,对每一个部分分别生成三维实体模型,再重新组合成一个实体整体,最终通过该实体整体生成人头颅的三维有限元模型;
其中,三部分具体如下:下颌骨分离形成单独的面网格模型,再将剩余部分分为颅顶和面中部两块。
作为优选,具体方法步骤如下:
S1.进行头颅CT扫描,得到原始CT数据,以DICOM格式导出;
S2.将获得的数据导入医学影像控制软件中,通过阈值分割,重建出人颅颌面骨面网格三维模型;
S3.使用图像分割工具将下颌骨分离形成单独的面网格模型,再将剩余部分分为颅顶和面中部两块,分别形成独立的面网格模型;
S4.分别将独立的三部分面网格模型导入逆向工程软件中修复畸形网格,消除多余的几何点、面,修复自相交的面和线,调整网格大小和质量;
S5.将修复好的面网格模型进行曲面生成,并修复面和缝隙,生成NUBRS曲面,并以此创建各部分三维CAD实体模型;
S6.将生成的三部分三维CAD实体模型分别导入软件,按照分割前的位置和形态进行组装,形成装配体;
S7.将装配体导入有限元前处理软件中,分别划分体网格,手动修复颅顶与面中部之间的缝隙,设置下颌骨髁状突与颅底关节窝为铰链连接。
本发明的有益效果如下:1、创建了头面部撞击伤三维有限元模型,能模拟出接近真实情况的骨折效果;2、为头面部撞击伤的研究提供完善的有限元分析工具;3、为复杂三维有限元模型的建模提供一种解决方案,即分割再组合的建模方法。
附图说明
图1为下颌骨枪弹伤三维有限元模型;
图2为A.Schaller,C.Voigt等完成的人颅颌面骨的三维有限元模型;
图3为人颅颌面骨面网格三维模型;
图4为软件分割后的人颅颌面骨面网格三维模型;
图5为面中部骨的面网格模型在逆向工程软件中进行修复和调整;
图6为颅顶骨的面网格模型在逆向工程软件中进行修复和调整;
图7为下颌骨面网格模型在逆向工程软件中进行修复和调整;
图8为面中部骨三维CAD实体模型;
图9为颅顶骨三维CAD实体模型;
图10为下颌骨三维CAD实体模型;
图11为颅顶、面中部、下颌骨三维实体模型形成装配体;
图12为人颅颌面骨完成体网格划分并生成三维有限元模型。
具体实施方式
以下结合附图,对本发明上述的和另外的技术特征和优点作更详细的说明。
一种基于人颅颌面骨的三维有限元模型的建模方法,包括以下步骤:
实施例:
(1)薄层CT扫描人头颅,得到原始CT数据,以DICOM格式导出;
(2)将dicom数据导入mimics15.0软件(Materialise'sinteractivemedicalimagecontrolsystem),通过阈值分割,(mimics15.0自带有自动阈值分割工具,选择bone选项,软件自动生成骨组织的mask,选择此mask后点击calculate3D即可)重建出人颅颌面骨面网格三维模型;
(3)使用图像分割工具将下颌骨分离形成单独的面网格模型,再将剩余部分分为颅顶和面中部两块,分别形成独立的面网格模型;
(4)分别将独立的三部分面网格模型导入逆向工程软件Geomagicstudio2013中修复畸形网格,消除多余的几何点、面,修复自相交的面和线,调整网格大小和质量;
(5)将修复好的面网格模型进行曲面生成,并修复狭长的面和缝隙,生成NUBRS曲面,并以此创建各部分三维CAD实体模型;
(6)将生成的三部分三维CAD实体模型分别导入solidworks软件,按照分割前的位置和形态进行组装,形成装配体;
(7)将装配体导入有限元前处理软件Hypermesh中,分别划分体网格,手动修复颅顶与面中部之间的缝隙,设置下颌骨髁状突与颅底关节窝为铰链连接。
其中,mimics15.0主要是将CT扫描所获取的dicom数据读取并显示,通过阈值分割功能将骨组织单独显示出来,并进行三维重建生成面网格三维模型即STL文件,把生成的STL文件分割成三部分;geomagic软件主要是对生成的STL文件进行面网格的重新划分,修正存在错误和缺陷的面网格,并将面网格模型转换成三维实体模型;solidworks主要是将分割后的三部分分别形成的三维实体模型进行装配还原成一个整体的三维实体模型;hypermesh是对装配还原成一个整体的三维实体模型进行网格划分和网格的修复,生成可用于有限元计算的三维有限元模型。因此,具有以上相应的功能的软件都是可替换的。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的实施方法,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。

Claims (2)

1.一种基于人颅颌面骨的三维有限元模型的建模方法,其特征在于,把人头颅分成三部分,对每一个部分分别生成三维实体模型,再重新组合成一个实体整体,最终通过该实体整体生成人头颅的三维有限元模型;
其中,三部分具体如下:下颌骨分离形成单独的面网格模型,再将剩余部分分为颅顶和面中部两块。
2.根据权利要求1所述的一种基于人颅颌面骨的三维有限元模型的建模方法,其特征在于,具体方法步骤如下:
S1.进行头颅CT扫描,得到原始CT数据,以DICOM格式导出;
S2.将获得的数据导入医学影像控制软件中,通过阈值分割,重建出人颅颌面骨面网格三维模型;
S3.使用图像分割工具将下颌骨分离形成单独的面网格模型,再将剩余部分分为颅顶和面中部两块,分别形成独立的面网格模型;
S4.分别将独立的三部分面网格模型导入逆向工程软件中修复畸形网格,消除多余的几何点、面,修复自相交的面和线,调整网格大小和质量;
S5.将修复好的面网格模型进行曲面生成,并修复面和缝隙,生成NUBRS曲面,并以此创建各部分三维CAD实体模型;
S6.将生成的三部分三维CAD实体模型分别导入软件,按照分割前的位置和形态进行组装,形成装配体;
S7.将装配体导入有限元前处理软件中,分别划分体网格,手动修复颅顶与面中部之间的缝隙,设置下颌骨髁状突与颅底关节窝为铰链连接。
CN201510751461.7A 2015-11-05 2015-11-05 一种基于人颅颌面骨的三维有限元模型的建模方法 Pending CN105373658A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510751461.7A CN105373658A (zh) 2015-11-05 2015-11-05 一种基于人颅颌面骨的三维有限元模型的建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510751461.7A CN105373658A (zh) 2015-11-05 2015-11-05 一种基于人颅颌面骨的三维有限元模型的建模方法

Publications (1)

Publication Number Publication Date
CN105373658A true CN105373658A (zh) 2016-03-02

Family

ID=55375853

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510751461.7A Pending CN105373658A (zh) 2015-11-05 2015-11-05 一种基于人颅颌面骨的三维有限元模型的建模方法

Country Status (1)

Country Link
CN (1) CN105373658A (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106934236A (zh) * 2017-03-09 2017-07-07 重庆医科大学附属口腔医院 基于锥形束ct的下颌骨运动模拟方法
CN107134010A (zh) * 2017-04-27 2017-09-05 杭州电子科技大学 一种基于有限元的弹性软组织的形貌效果预测方法
CN108210072A (zh) * 2018-02-08 2018-06-29 扈玉华 基于mri和cta的脑组织及血管实体复合模型的制备方法
CN109686449A (zh) * 2017-10-18 2019-04-26 香港中文大学深圳研究院 唇腭裂病人的生物力学仿真模型的构建方法
CN110096728A (zh) * 2019-03-12 2019-08-06 昆明理工大学 一种基于逆向重构的藕状多孔金属有限元分析方法
CN110689604A (zh) * 2019-05-10 2020-01-14 腾讯科技(深圳)有限公司 个性化脸部模型显示方法、装置、设备及存储介质
CN113870267A (zh) * 2021-12-03 2021-12-31 深圳市奥盛通科技有限公司 缺陷检测方法、装置、计算机设备及可读存储介质

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104123416A (zh) * 2014-07-21 2014-10-29 中国医学科学院生物医学工程研究所 一种模拟真实人体颅脑电特性分布的有限元仿真模型
CN104915519A (zh) * 2015-06-30 2015-09-16 中国人民解放军第三军医大学第二附属医院 一种颅颌面骨模型建立方法及装置

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104123416A (zh) * 2014-07-21 2014-10-29 中国医学科学院生物医学工程研究所 一种模拟真实人体颅脑电特性分布的有限元仿真模型
CN104915519A (zh) * 2015-06-30 2015-09-16 中国人民解放军第三军医大学第二附属医院 一种颅颌面骨模型建立方法及装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
康嘉晏 等: "含骨缝犬单侧牙槽突裂三维有限元模型的而建立", 《中国口腔颌面外科杂志》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106934236A (zh) * 2017-03-09 2017-07-07 重庆医科大学附属口腔医院 基于锥形束ct的下颌骨运动模拟方法
CN107134010A (zh) * 2017-04-27 2017-09-05 杭州电子科技大学 一种基于有限元的弹性软组织的形貌效果预测方法
CN109686449A (zh) * 2017-10-18 2019-04-26 香港中文大学深圳研究院 唇腭裂病人的生物力学仿真模型的构建方法
CN108210072A (zh) * 2018-02-08 2018-06-29 扈玉华 基于mri和cta的脑组织及血管实体复合模型的制备方法
CN110096728A (zh) * 2019-03-12 2019-08-06 昆明理工大学 一种基于逆向重构的藕状多孔金属有限元分析方法
CN110689604A (zh) * 2019-05-10 2020-01-14 腾讯科技(深圳)有限公司 个性化脸部模型显示方法、装置、设备及存储介质
CN110689604B (zh) * 2019-05-10 2023-03-10 腾讯科技(深圳)有限公司 个性化脸部模型显示方法、装置、设备及存储介质
CN113870267A (zh) * 2021-12-03 2021-12-31 深圳市奥盛通科技有限公司 缺陷检测方法、装置、计算机设备及可读存储介质

Similar Documents

Publication Publication Date Title
CN105373658A (zh) 一种基于人颅颌面骨的三维有限元模型的建模方法
CN103919631B (zh) 一种颌骨缺损个性化修复体的制备方法
Mohammed et al. Advanced auricular prosthesis development by 3D modelling and multi-material printing
CN104985183B (zh) 一种低弹性模量钛基颌骨植入体及其制备方法
Keeve et al. Anatomy-based facial tissue modeling using the finite element method
CN105997306B (zh) 一种骨骼植入体内填充多孔网格结构的设计方法
CN101953726B (zh) 一种个性化颅骨修复体的制备方法
CN104915519B (zh) 一种颅颌面骨模型建立方法及装置
CN105303605A (zh) 一种基于力反馈的骨外科手术仿真系统
JP2014516169A (ja) 歯内治療シミュレーション・システム
Wu et al. A virtual training system for maxillofacial surgery using advanced haptic feedback and immersive workbench
CN104367389A (zh) 一种义齿数据处理及三维打印系统
CN107607452A (zh) 一种按应力分布确定多孔假体孔隙率的方法
CN102169519B (zh) 一种腭裂牙槽骨塑形矫治器数字化方法
CN107134010A (zh) 一种基于有限元的弹性软组织的形貌效果预测方法
CN106991720A (zh) 一种基于有限元分析运算的个性化髋臼重建钢板预弯方法
CN112070898A (zh) 一种高质量网格小型猪头部有限元模型建模方法
CN218676304U (zh) 一种三维可视化阻生牙拔牙模型
KR20130136221A (ko) 변형체의 실시간 시뮬레이션을 위한 다해상도 무요소법
Chopade et al. Finite element analysis: new dimension in prosthodontic research
Jiang et al. Modeling a realistic 3D physiological tongue for visual speech synthesis
Singare et al. The use of 3D printing technology in human defect reconstruction-a review of cases study
Park Topology Optimization to Design Bone Replacement Shapes in Craniofacial Reconstructive Surgery: Design, Simulation and Experimental Validation
Rossi et al. Computer graphics applied to Anatomy: A study of two Bio-CAD modeling methods on finite element analysis of human edentulous hemi-mandible
Budsberg et al. Elastic and plastic deformations with rigid body dynamics.

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20160302