CN110236741A - 一种具有拓扑优化固定单元和多孔髁突头单元的个性化髁突假体设计方法及个性化髁突假体 - Google Patents

一种具有拓扑优化固定单元和多孔髁突头单元的个性化髁突假体设计方法及个性化髁突假体 Download PDF

Info

Publication number
CN110236741A
CN110236741A CN201910299195.7A CN201910299195A CN110236741A CN 110236741 A CN110236741 A CN 110236741A CN 201910299195 A CN201910299195 A CN 201910299195A CN 110236741 A CN110236741 A CN 110236741A
Authority
CN
China
Prior art keywords
condyle
porous
model
fixed cell
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.)
Granted
Application number
CN201910299195.7A
Other languages
English (en)
Other versions
CN110236741B (zh
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.)
Zhejiang University of Technology ZJUT
Original Assignee
Zhejiang University of Technology ZJUT
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 Zhejiang University of Technology ZJUT filed Critical Zhejiang University of Technology ZJUT
Priority to CN201910299195.7A priority Critical patent/CN110236741B/zh
Publication of CN110236741A publication Critical patent/CN110236741A/zh
Application granted granted Critical
Publication of CN110236741B publication Critical patent/CN110236741B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61FFILTERS IMPLANTABLE INTO BLOOD VESSELS; PROSTHESES; DEVICES PROVIDING PATENCY TO, OR PREVENTING COLLAPSING OF, TUBULAR STRUCTURES OF THE BODY, e.g. STENTS; ORTHOPAEDIC, NURSING OR CONTRACEPTIVE DEVICES; FOMENTATION; TREATMENT OR PROTECTION OF EYES OR EARS; BANDAGES, DRESSINGS OR ABSORBENT PADS; FIRST-AID KITS
    • A61F2/00Filters implantable into blood vessels; Prostheses, i.e. artificial substitutes or replacements for parts of the body; Appliances for connecting them with the body; Devices providing patency to, or preventing collapsing of, tubular structures of the body, e.g. stents
    • A61F2/02Prostheses implantable into the body
    • A61F2/30Joints
    • A61F2/3094Designing or manufacturing processes
    • A61F2/30942Designing or manufacturing processes for designing or making customized prostheses, e.g. using templates, CT or NMR scans, finite-element analysis or CAD-CAM techniques
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61FFILTERS IMPLANTABLE INTO BLOOD VESSELS; PROSTHESES; DEVICES PROVIDING PATENCY TO, OR PREVENTING COLLAPSING OF, TUBULAR STRUCTURES OF THE BODY, e.g. STENTS; ORTHOPAEDIC, NURSING OR CONTRACEPTIVE DEVICES; FOMENTATION; TREATMENT OR PROTECTION OF EYES OR EARS; BANDAGES, DRESSINGS OR ABSORBENT PADS; FIRST-AID KITS
    • A61F2/00Filters implantable into blood vessels; Prostheses, i.e. artificial substitutes or replacements for parts of the body; Appliances for connecting them with the body; Devices providing patency to, or preventing collapsing of, tubular structures of the body, e.g. stents
    • A61F2/02Prostheses implantable into the body
    • A61F2/30Joints
    • A61F2/30988Other joints not covered by any of the groups A61F2/32 - A61F2/4425
    • A61F2/3099Other joints not covered by any of the groups A61F2/32 - A61F2/4425 for temporo-mandibular [TM, TMJ] joints
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61FFILTERS IMPLANTABLE INTO BLOOD VESSELS; PROSTHESES; DEVICES PROVIDING PATENCY TO, OR PREVENTING COLLAPSING OF, TUBULAR STRUCTURES OF THE BODY, e.g. STENTS; ORTHOPAEDIC, NURSING OR CONTRACEPTIVE DEVICES; FOMENTATION; TREATMENT OR PROTECTION OF EYES OR EARS; BANDAGES, DRESSINGS OR ABSORBENT PADS; FIRST-AID KITS
    • A61F2/00Filters implantable into blood vessels; Prostheses, i.e. artificial substitutes or replacements for parts of the body; Appliances for connecting them with the body; Devices providing patency to, or preventing collapsing of, tubular structures of the body, e.g. stents
    • A61F2/02Prostheses implantable into the body
    • A61F2/30Joints
    • A61F2/3094Designing or manufacturing processes
    • A61F2/30942Designing or manufacturing processes for designing or making customized prostheses, e.g. using templates, CT or NMR scans, finite-element analysis or CAD-CAM techniques
    • A61F2002/30952Designing or manufacturing processes for designing or making customized prostheses, e.g. using templates, CT or NMR scans, finite-element analysis or CAD-CAM techniques using CAD-CAM techniques or NC-techniques
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61FFILTERS IMPLANTABLE INTO BLOOD VESSELS; PROSTHESES; DEVICES PROVIDING PATENCY TO, OR PREVENTING COLLAPSING OF, TUBULAR STRUCTURES OF THE BODY, e.g. STENTS; ORTHOPAEDIC, NURSING OR CONTRACEPTIVE DEVICES; FOMENTATION; TREATMENT OR PROTECTION OF EYES OR EARS; BANDAGES, DRESSINGS OR ABSORBENT PADS; FIRST-AID KITS
    • A61F2/00Filters implantable into blood vessels; Prostheses, i.e. artificial substitutes or replacements for parts of the body; Appliances for connecting them with the body; Devices providing patency to, or preventing collapsing of, tubular structures of the body, e.g. stents
    • A61F2/02Prostheses implantable into the body
    • A61F2/30Joints
    • A61F2/3094Designing or manufacturing processes
    • A61F2/30942Designing or manufacturing processes for designing or making customized prostheses, e.g. using templates, CT or NMR scans, finite-element analysis or CAD-CAM techniques
    • A61F2002/30955Designing or manufacturing processes for designing or making customized prostheses, e.g. using templates, CT or NMR scans, finite-element analysis or CAD-CAM techniques using finite-element analysis
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61FFILTERS IMPLANTABLE INTO BLOOD VESSELS; PROSTHESES; DEVICES PROVIDING PATENCY TO, OR PREVENTING COLLAPSING OF, TUBULAR STRUCTURES OF THE BODY, e.g. STENTS; ORTHOPAEDIC, NURSING OR CONTRACEPTIVE DEVICES; FOMENTATION; TREATMENT OR PROTECTION OF EYES OR EARS; BANDAGES, DRESSINGS OR ABSORBENT PADS; FIRST-AID KITS
    • A61F2/00Filters implantable into blood vessels; Prostheses, i.e. artificial substitutes or replacements for parts of the body; Appliances for connecting them with the body; Devices providing patency to, or preventing collapsing of, tubular structures of the body, e.g. stents
    • A61F2/02Prostheses implantable into the body
    • A61F2/30Joints
    • A61F2/3094Designing or manufacturing processes
    • A61F2002/30985Designing or manufacturing processes using three dimensional printing [3DP]

Abstract

一种具有拓扑优化固定单元和多孔髁突头单元的个性化髁突假体设计方法,包括以下步骤:1.)髁突假体重建和下颌骨生物力学模型构建;2.)基于生物力学的固定单元的拓扑优化设计;3.)髁突头多孔结构设计;4.)对个性化髁突假体进行3D打印;5.)对打印后的个性化髁突假体进行后处理,得到应用于临床的具有拓扑优化固定单元和多孔髁突头单元的个性化髁突假体。本发明提供了一种具有拓扑优化固定单元和多孔髁突头单元的个性化髁突假体设计方法及个性化髁突假体,有效降低了假体的平均最大应力,稳定性更好,不易脱落。

Description

一种具有拓扑优化固定单元和多孔髁突头单元的个性化髁突 假体设计方法及个性化髁突假体
技术领域
本发明涉及人工颞下颌关节技术领域,尤其是涉及一种具有拓扑优化固定单元和多孔髁突头单元的个性化髁突假体设计方法及髁突假体。
背景技术
颞下颌关节位于人体下颌骨和颧骨连接处,是人体重要的关节。颞下颌关节疾病是发病率较高的疾病,但大部分不需要进行手术甚至不需要进行治疗,只会生发关节弹响或张口微微受限甚至没有表现病症的情况。然而有少部分发病较为严重,患者会出现不能张口和疼痛的症状,例如颞下颌关节强直,髁突已经和关节窝连结到一起,下颌骨不能运动。这种情况下需要截骨,然后通过关节重建进行治疗。颞下颌关节重建包括了自体骨移植,牵引成骨,髁突假体置换和全关节假体置换等。其中髁突假体置换是重要的治疗手段,在临床应用了很多年,具有较好的临床效果。
不同的人工髁突假体,其治疗效果也不尽相同,随着假体设计的不断发展,假体结构和力学性能越来越好,在人体内的寿命和稳定性也在提高。特别是目前商用的颞下颌关节假体,已经取得了广泛的市场认可,如美国concepts公司的颞下颌关节假体。但是,目前商用和论文研究的假体,其髁突头部分都是实心或者半空心的结构,固定板部分螺钉的数量和位置也各有不同。有术后随访研究发现部分案例出现一些并发症,包括钛板松动和暴露,张口受限、疼痛等。从力学和生物学角度分析,其与固定板的受力,螺钉的应力疲劳及应力集中有一定关系,且螺钉周围的下颌骨部分可能因应力过大而发生骨萎缩,导致螺钉松动,重建失败。
同样是人工关节假体,与髁突关节假体相比,髋关节假体近些年的研究已经相当成熟。髋关节由股骨头与髋臼构成。目前三维多孔结构的髋关节假体已经比较常见,其结构能有效的降低弹性模量,从而减少因“应力遮挡”导致的骨吸收情况,可见三维多孔结构的关节假体具有相当大的优势。近年来也有关于三维多孔结构的髁突假体设计,例如将髁突头部分设计成多孔结构,能显著降低“应力遮挡”,但是这些研究没有充分考虑生物学和生物力学因素。此外,目前研究设计的假体中,其固定板螺钉孔位和数量的选择也各有不同。而假体固定板部分和髁突头部分的优异性能都对提高假体寿命和稳定性有效果。因此设计一种符合下颌骨生物力学性能要求的多孔结构且固定方式较优的人工髁突假体能提高假体在人体内的稳定性及寿命。
多孔结构设计方法中,用拓扑优化方法设计多孔骨植入体的研究已经相当多。拓扑优化旨在通过对结构的尺寸、形状及拓扑等参数进行合理的调整,使得调整后的结构能够在满足强度、刚度、稳定性、可制造性以及在其它一种或多种设计要求的前提下,特定的目标性能达到最优,比如重量最轻、造价最低等。一般采用拓扑优化方法对多孔结构的微观单元结构进行设计,限定模型的体积分数,再以弹性模量为目标函数,完成对微结构的拓扑优化设计。
增材制造技术(或称为3D打印技术)的出现与不断发展,使得医疗手术由传统的纯经验公式转向数字化精确手术。随着3D打印技术在数字化医学领域的应用,临床医学领域必定有相当大的颠覆和发展。 3D打印应用于医学领域的优势在于其对任意复杂几何结构模型的适应性,使得数字化设计的个性化医疗器械可以通过3D打印实现快速制造。因此,3D打印技术为个性化多孔结构的髁突假体的制备提供了手段。
发明内容
为了克服现有技术存在的缺陷,本发明提供了一种具有拓扑优化固定单元和多孔髁突头单元的个性化髁突假体设计方法及个性化髁突假体,有效降低了假体的平均最大应力,稳定性更好,不易脱落。
本发明解决其技术问题所采用的技术方案是:
一种具有拓扑优化固定单元和多孔髁突头单元的个性化髁突假体设计方法,包括以下步骤:
1.)髁突假体重建和下颌骨生物力学模型构建
1.1)对颞下颌关节患者进行计算机断层扫描,获取其CT数据,利用医学图像处理软件进行医学图像处理和模型重建,得到下颌骨三维表面模型,并同时重建髁突四面体实体网格模型;
1.2)利用医学图像处理软件模拟截骨,切除髁突,在下颌骨髁突颈下部插入平面,将下颌骨和髁突分离;
1.3)利用造型软件设计髁突的固定单元,提取下颌骨表面,生成厚为2mm的薄板模型,并对薄板的边缘进行适当修剪和光滑处理,并在网格处理软件中,得到初始固定单元四面体实体网格模型;
1.4)由分离的髁突和初始固定单元组成髁突假体;
1.5)在造型软件中,基于下颌骨三维表面模型,根据图像的灰度值与骨密度的关系函数,计算出下颌骨不同区域的骨密度值,再根据骨密度与杨氏模量的关系函数计算出各体素的杨氏模量,完成下颌骨实体网格模型的材料参数赋值,从而得到非均匀的下颌骨模型;
1.6)将下颌骨的升颌肌群简化为一维拉伸弹簧,确定弹簧刚度值,弹簧连接点为肌肉在下颌骨附着面积的中心点,并指向颅颌骨上的连接点;
1.7)对两个髁突的后上部位置的六个自由度施加相应的边界约束,得到下颌骨生物力学模型;
1.8)在有限元仿真软件中,将初始固定单元四面体实体网格模型与下颌骨生物力学模型进行装配,髁突下表面和下颌骨断端表面设定接触;初始固定单元分别和髁突、下颌骨绑定;
2.)基于生物力学的固定单元的拓扑优化设计
2.1)在下颌骨不同牙位加载正常咬合力,分别为三种咬合工况, 加载肌肉力,并进行有限元仿真,分别得到下颌骨和固定单元的最大绝对主应力分布,并记录最大应力值;
2.2)将赋有材料属性的患者下颌骨实体网格模型和初始固定单元四面体实体网格模型一起导入有限元仿真软件中,设置优化目标,利用有限元仿真软件的拓扑优化功能进行优化,得到固定单元的初始优化形状;
2.3)将拓扑优化得到的固定单元导出到网格处理软件中修整,得到固定单元优化模型,并在该固定单元优化模型上,根据下颌骨解剖结构,避开重要解剖结构,确定螺钉的安装位置和数目,并设计出一套完整结构的优化固定单元;
2.4)用与步骤2.1)相同的加载方式对优化固定单元进行有限元仿真,根据得到的下颌骨和固定单元的应力分布和应力值,评价该固定单元是否满足治疗要求,如果满足要求,则进入下一步,如果不满足,则返回步骤2.2)重新进行优化设计,直到最终的优化固定单元满足要求;
2.5)将分离的髁突模型和最终的优化固定单元外表面进行连接,并使得优化固定单元外表面和髁突头侧面光滑连接,得到实心的人工髁突假体数字模型;
3.)髁突头多孔结构设计
3.1)基于下颌骨生物特征,皮质骨和松质骨力学性能和孔隙率的差异,髁突头部分多孔结构的晶胞密度将依次分为两个梯度,外层皮质骨部分和内层松质骨部分,其中内层松质骨部分的晶胞密度更小,获得更大的晶胞尺寸;
3.2)将髁突头模型区域分割,利用自然骨生物特征,通过将模型切割成小块,控制体网格生成质量;利用医学图像处理软件,通过阈值调整提取中间松质骨,经过造型软件处理得到较光滑的皮质骨-松质骨分界面;
3.3)从实心的人工髁突假体数字模型中沿断面再次分离出髁突头;
3.4)利用造型软件,建立三层曲面,髁突外表面、皮质骨-松质骨分界面和两面之间建立NURBS曲面模型;
3.5)外层皮质骨部分和内层松质骨部分根据曲面曲率特征和重新建立曲面的需要,将曲面分成几个小曲面分别重建,再将直线投影在原外轮廓面上,生成曲线,经过一系列调整,建立高质量的uv曲线,重新生成光滑的曲面;中间层由外层曲面向内偏移,出现破面,再用同样的方法在破面上建立uv曲线,重新生成曲面;
3.6)建立有密度分布的多孔结构晶胞模型,利用有限元前处理软件,进行体网格划分,并生成正四面体网格即为多孔晶胞单元,将多孔晶胞单元编号,附属节点编号,并将节点坐标导出;
3.7)将赋有材料属性的患者下颌骨实体网格模型和实心的人工髁突假体数字模型一起导入有限元仿真软件中,利用有限元仿真软件的拓扑优化功能进行优化,得每个单元的相对单元密度;将拓扑优化体单元编号、附属节点、节点坐标,和节点相对单元密度值导出;
3.8)将多孔结构晶胞模型和拓扑优化体单元导出的数据在Excel 中进行处理,把相对单元密度数据映射到多孔晶胞单元上,找出拓扑优化模型中离晶胞节点最近的节点,将其密度值赋予此晶胞节点,最终得到所有晶胞节点的密度值;
3.9)根据多孔晶胞单元四个节点的密度值计算晶胞的密度,赋予每个多孔晶胞单元密度数据;
3.10)根据多孔晶胞单元密度,计算出晶胞单元内六根杆的杆径,最终杆径取其平均值;
3.11)已知多孔晶胞单元所有节点的坐标及杆径,利用ug二次开发平台,导入软件成多孔髁突头模型;
3.12)将多孔髁突头模型与优化固定单元通过布尔操作组合成个性化髁突假体;
4.)对个性化髁突假体进行3D打印;
5.)对打印后的个性化髁突假体进行后处理,得到应用于临床的具有拓扑优化固定单元和多孔髁突头单元的个性化髁突假体。
进一步,所述步骤1.1)中,髁突形状如果基本完整直接重建,或者一边健康可用镜像生成,或者都不健康,则可以根据关节窝形状及下颌骨形状以及健康人体数据库,重新设计一个人工关节外形,然后利用网格处理软件生成四面体实体网格模型。
一种具有拓扑优化固定单元和多孔髁突头单元的个性化髁突假体设计方法构建的个性化髁突假体,包括固定单元和多孔髁突头单元,所述固定单元为具有拓扑优化结构的固定板,所述多孔髁突头单元设置在固定板的上部上,所述多孔髁突头单元为基于密度分布晶胞多孔结构的髁突头,该髁突头包括外层皮质骨部分和内层松质骨部分,每层都是由多个多孔晶胞单元构成,每个多孔晶胞单元均为一个由六根杆组成的正四面体结构,内层松质骨部分的多孔晶胞单元的密度小于外层皮质骨部分的多孔晶胞单元的密度,外层皮质骨部分的多孔晶胞单元的杆长度小于内层松质骨部分的多孔晶胞单元的杆长度;所述固定板的下部上设有螺钉孔。
进一步,所述固定板为倒V型,所述螺钉孔设置有四个,两两一组,左右对称的布置在V型的左右两侧上。
本发明的有益效果主要表现在:设计了具有非均匀密度晶胞的且孔隙率可定制的多孔髁突头单元,杆径大小适应应力分布;关节假体的固定板部分,得到了四颗螺钉固定的最优方案,并通过有限元方法进行应力比较;结果显示,优化的固定板和螺钉的应力比三颗螺钉固定方案下降了9%~49%;分析比较了三种咬合方式下,所设计的多孔髁突头假体和实心髁突头假体的生物力学性能,有限元结果显示,多孔髁突头假体的平均最大应力比实心髁突头假体的平均最大应力下降了48%。
附图说明
图1是本发明的常用颞下颌关节切除手术简化示意图。
图2是本发明的提取下颌骨表面生成2mm固定单元的示意图。
图3是本发明的修整边缘后的固定板示意图。
图4是本发明的患者下颌骨生物力学模型的简化图。
图5是本发明的拓扑优化后的固定板初始优化模型。
图6是本发明的髁突常用外固定系统的简化示意图。
图7是本发明的具有拓扑优化固定单元的个性化实心髁突假体模型。
图8是本发明的利用拓扑优化方法设计固定单元的流程图。
图9是本发明的髁突头三层曲面模型。
图10是本发明的多孔结构晶胞模型。
图11是本发明的基于密度分布晶胞多孔髁突头单元示意图。
图12是本发明的个性化髁突假体示意图。
图13是本发明的基于密度分布晶胞多孔髁突头单元的设计流程图。
具体实施方式
下面结合附图对本发明作进一步描述。
参照图1~图13,一种具有拓扑优化固定单元和多孔髁突头单元的个性化髁突假体设计方法,包括以下步骤:
1.)髁突假体重建和下颌骨生物力学模型构建
1.1)对颞下颌关节患者进行计算机断层扫描,获取其CT数据,利用医学图像处理软件如Mimics,进行医学图像处理和模型重建,得到下颌骨三维表面模型(空心的stl格式),髁突形状如果基本完整直接重建,或者一边健康可用镜像生成,或者都不健康,则可以根据关节窝形状及下颌骨形状以及健康人体数据库,重新设计一个人工髁突外形,然后利用网格处理软件如3-matic,生成四面体实体网格模型。
1.2)利用Mimics软件Cut功能模拟截骨,切除髁突,在下颌骨髁突颈下部插入平面1,平面厚度为0.2mm,截骨位置尽量符合真实案例,将下颌骨和髁突分离。
1.3)利用造型软件如magics设计髁突固定单元,根据参考文献和现有产品假体固定单元的位置和大小,提取下颌骨表面,生成厚为 2mm的薄板模型3,并对薄板的边缘进行适当修剪和光滑处理,并在网格处理软件中,得到初始固定单元四面体实体网格模型。
1.4)髁突假体的基本外形结构如图所示,由分离的髁突和固定单元组成。
1.5)在造型软件中,基于下颌骨三维表面模型,根据图像的灰度值与骨密度的关系函数,计算出下颌骨不同区域的骨密度值,再根据骨密度与杨氏模量的关系函数计算出各体素的杨氏模量,完成下颌骨实体网格模型的材料参数赋值,从而得到非均匀的下颌骨模型2。如图4所示,升颌肌群包括翼外肌21、颞肌28、咬肌23、翼内肌27等, 24为患侧磨牙,26为健侧磨牙,25为前牙。
1.6)将下颌骨的升颌肌群简化为一维拉伸弹簧,弹簧刚度值根据肌电检测数据或文献资料中的参数确定,弹簧连接点为肌肉在下颌骨附着面积的中心点,并指向颅颌骨上的连接点。
1.7)对两个髁突的后上部位置的六个自由度施加相应的边界约束,得到下颌骨生物力学模型。
(1.8)在有限元仿真软件中,将初始固定单元四面体实体网格模型与上述的下颌骨生物力学模型进行装配,髁突下表面和下颌骨断端表面设定接触;初始固定单元分别和髁突、下颌骨绑定。
2.)基于生物力学的固定单元的拓扑优化设计
2.1)在下颌骨不同牙位加载正常咬合力,分别为三种咬合工况, 加载肌肉力,并进行有限元仿真,分别得到下颌骨和固定单元的最大绝对主应力分布,并记录最大应力值。
2.2)将赋有材料属性的患者下颌骨实体网格模型和初始固定单元四面体实体网格模型一起导入有限元仿真软件中,设置优化目标为应变能最小化,即刚度最大化,优化约束为体积小于70%和50%,在两种约束条件下,得到的优化结果几乎一致,利用有限元仿真软件的拓扑优化功能进行优化,得到固定单元的初始优化形状。
2.3)将拓扑优化得到的固定单元导出到网格处理软件如 3-matic中修整,得到固定单元优化模型,并在该固定单元优化模型上,根据下颌骨解剖结构,考虑避开重要解剖结构,确定螺钉的安装位置和数目,并设计出一套完整结构的优化固定单元。
2.4)用与步骤2.1)相同的加载方式对优化固定单元进行有限元仿真,根据得到的下颌骨和固定单元的应力分布和应力值,评价该固定单元是否满足治疗要求,如果满足要求,则进入下一步,如果不满足,则返回步骤2.2)重新进行优化设计,直到最终的优化固定单元满足要求。
2.5)将分离的髁突模型,为stl格式和最终的优化固定单元外表面进行连接,删除外表面和髁突侧面相邻处的三角面片,通过修补、局部光滑等操作使得优化固定单元外表面和髁突头侧面光滑连接,得到实心的人工髁突假体数字模型。
3.)髁突头多孔结构设计
3.1)基于下颌骨生物特征,皮质骨和松质骨力学性能和孔隙率的差异,髁突头部分多孔结构的晶胞密度将依次分为两个梯度,外层皮质骨部分和内层松质骨部分,其中内层松质骨部分的晶胞密度更小,获得更大的晶胞尺寸。
3.2)将髁突头模型区域分割,利用自然骨生物特征,通过将模型切割成小块,更好地控制体网格生成质量;利用Mimics软件,通过阈值调整提取中间松质骨,经过magics软件处理得到较光滑的皮质骨 -松质骨分界面。
3.3)从实心的人工髁突假体数字模型中沿断面再次分离出髁突头。
3.4)利用造型软件如Rhino软件,建立三层曲面,髁突外表面、皮质骨-松质骨分界面和两面之间建立NURBS曲面模型。
3.5)外层皮质骨部分和内层松质骨部分根据曲面曲率特征和重新建立曲面的需要,将曲面分成几个小曲面分别重建,即生成几个相交平面将曲面分割成五个部分(不包括底面),再将直线投影在原外轮廓面上,生成曲线,经过一系列调整,建立高质量的uv曲线,重新生成光滑的曲面;中间层由外层曲面向内偏移,出现破面,再用同样的方法在破面上建立uv曲线,重新生成曲面。
3.6)建立有密度分布的多孔结构晶胞模型,利用有限元前处理软件,如hypermesh软件,进行体网格划分,其网格尺寸和体网格密度从最外层曲面到最里层依次都设置为1.2、1.2、1.4,即三角形网格边长,生成的正四面体网格即为多孔晶胞单元,将多孔晶胞单元编号,附属节点编号,节点坐标以txt格式导出。
3.7)将赋有材料属性的患者下颌骨实体网格模型和实心的人工髁突假体数字模型一起导入有限元仿真软件中,设置优化目标为应变能最小化,即刚度最大化,优化约束为体积小于70%,利用有限元仿真软件的拓扑优化功能进行优化,得每个单元的相对单元密度;将拓扑优化体单元编号、附属节点、节点坐标,和节点相对单元密度值以 txt格式导出。
3.8)将多孔结构晶胞模型和拓扑优化体单元导出的数据在 Excel中进行处理,把相对单元密度数据映射到多孔晶胞单元上,找出拓扑优化模型中离晶胞节点最近的节点,将其密度值赋予此晶胞节点,最终得到所有晶胞节点的密度值。
3.9)根据多孔晶胞单元四个节点的密度值计算晶胞的密度,赋予每个多孔晶胞单元密度数据。
3.10)根据多孔晶胞单元密度,计算出晶胞单元内六根杆的杆径,而每根杆一般由多个晶胞共用,所以有多个杆径值,最终杆径取其平均值。
3.11)已知多孔晶胞单元所有节点的坐标及杆径,利用ug二次开发平台,编写可以利用此数据生成圆柱状杆多孔结构模型的C++程序,导入软件成多孔髁突头模型。
3.12)将多孔髁突头模型与优化固定单元通过布尔操作组合成一种具有拓扑优化固定单元和多孔关节头单元的个性化髁突假体。
4.)利用金属3D打印SLM对钛合金粉末进行3D打印,得具有拓扑优化固定单元和髁突头单元的个性化髁突假体。
5.)对打印后的个性化髁突假体进行抛光、超声清洗等后处理,得到应用于临床的具有拓扑优化固定单元和多孔髁突头单元的个性化髁突假体。
一种具有拓扑优化固定单元和多孔髁突头单元的个性化髁突假体包括固定单元和多孔髁突头单元,所述固定单元为具有拓扑优化结构的固定板111,所述多孔髁突头单元设置在固定板的上部上,所述多孔髁突头单元为基于密度分布晶胞多孔结构的髁突头10,该髁突头包括外层皮质骨部分和内层松质骨部分,每层都是由多个多孔晶胞单元构成,每个多孔晶胞单元均为一个由六根杆组成的正四面体结构,内层松质骨部分的多孔晶胞单元的密度小于外层皮质骨部分的多孔晶胞单元的密度,外层皮质骨部分的多孔晶胞单元的杆长度小于内层松质骨部分的多孔晶胞单元的杆长度;所述固定板的下部上设有螺钉孔61。
进一步,所述固定板为倒V型,所述螺钉孔61设置有四个,两两一组,左右对称的布置在V型的左右两侧上。
如图3所示,4为修整边缘后的固定板,图5中,5为拓扑优化后的固定板初始化模型,51为固定板去除区域;图6中,6为初始化模型修整边缘后的固定板,61为螺钉孔,图7中,7为具有固定单元的个性化实心髁突假体模型,71为实心髁突,72为具有拓扑优化的固定单元;图9中,8为髁突头三层曲面模型,82为中间面,81为髁突外表面,83为皮质骨-松质骨分界面,图10中的9为多孔晶胞单元;图11中,10为多孔髁突头单元,101为外层皮质骨部分,102为内层松质骨部分;图12中,11为个性化髁突假体,111为具有拓扑优化结构的固定板。

Claims (4)

1.一种具有拓扑优化固定单元和多孔髁突头单元的个性化髁突假体设计方法,其特征在于:所述设计方法包括以下步骤:
1.)髁突假体重建和下颌骨生物力学模型构建
1.1)对颞下颌关节患者进行计算机断层扫描,获取其CT数据,利用医学图像处理软件进行医学图像处理和模型重建,得到下颌骨三维表面模型,并同时重建髁突四面体实体网格模型;
1.2)利用医学图像处理软件模拟截骨,切除髁突,在下颌骨髁突颈下部插入平面,将下颌骨和髁突分离;
1.3)利用造型软件设计髁突的固定单元,提取下颌骨表面,生成厚为2mm的薄板模型,并对薄板的边缘进行适当修剪和光滑处理,并在网格处理软件中,得到初始固定单元四面体实体网格模型;
1.4)由分离的髁突和初始固定单元组成髁突假体;
1.5)在造型软件中,基于下颌骨三维表面模型,根据图像的灰度值与骨密度的关系函数,计算出下颌骨不同区域的骨密度值,再根据骨密度与杨氏模量的关系函数计算出各体素的杨氏模量,完成下颌骨实体网格模型的材料参数赋值,从而得到非均匀的下颌骨模型;
1.6)将下颌骨的升颌肌群简化为一维拉伸弹簧,确定弹簧刚度值,弹簧连接点为肌肉在下颌骨附着面积的中心点,并指向颅颌骨上的连接点;
1.7)对两个髁突的后上部位置的六个自由度施加相应的边界约束,得到下颌骨生物力学模型;
1.8)在有限元仿真软件中,将初始固定单元四面体实体网格模型与下颌骨生物力学模型进行装配,髁突下表面和下颌骨断端表面设定接触;初始固定单元分别和髁突、下颌骨绑定;
2.)基于生物力学的固定单元的拓扑优化设计
2.1)在下颌骨不同牙位加载正常咬合力,分别为三种咬合工况,加载肌肉力,并进行有限元仿真,分别得到下颌骨和固定单元的最大绝对主应力分布,并记录最大应力值;
2.2)将赋有材料属性的患者下颌骨实体网格模型和初始固定单元四面体实体网格模型一起导入有限元仿真软件中,设置优化目标,利用有限元仿真软件的拓扑优化功能进行优化,得到固定单元的初始优化形状;
2.3)将拓扑优化得到的固定单元导出到网格处理软件中修整,得到固定单元优化模型,并在该固定单元优化模型上,根据下颌骨解剖结构,避开重要解剖结构,确定螺钉的安装位置和数目,并设计出一套完整结构的优化固定单元;
2.4)用与步骤2.1)相同的加载方式对优化固定单元进行有限元仿真,根据得到的下颌骨和固定单元的应力分布和应力值,评价该固定单元是否满足治疗要求,如果满足要求,则进入下一步,如果不满足,则返回步骤2.2)重新进行优化设计,直到最终的优化固定单元满足要求;
2.5)将分离的髁突模型和最终的优化固定单元外表面进行连接,并使得优化固定单元外表面和髁突头侧面光滑连接,得到实心的人工髁突假体数字模型;
3.)髁突头多孔结构设计
3.1)基于下颌骨生物特征,皮质骨和松质骨力学性能和孔隙率的差异,髁突头部分多孔结构的晶胞密度将依次分为两个梯度,外层皮质骨部分和内层松质骨部分,其中内层松质骨部分的晶胞密度更小,获得更大的晶胞尺寸;
3.2)将髁突头模型区域分割,利用自然骨生物特征,通过将模型切割成小块,控制体网格生成质量;利用医学图像处理软件,通过阈值调整提取中间松质骨,经过造型软件处理得到较光滑的皮质骨-松质骨分界面;
3.3)从实心的人工髁突假体数字模型中沿断面再次分离出髁突头;
3.4)利用造型软件,建立三层曲面,髁突外表面、皮质骨-松质骨分界面和两面之间建立NURBS曲面模型;
3.5)外层皮质骨部分和内层松质骨部分根据曲面曲率特征和重新建立曲面的需要,将曲面分成几个小曲面分别重建,再将直线投影在原外轮廓面上,生成曲线,经过一系列调整,建立高质量的uv曲线,重新生成光滑的曲面;中间层由外层曲面向内偏移,出现破面,再用同样的方法在破面上建立uv曲线,重新生成曲面;
3.6)建立有密度分布的多孔结构晶胞模型,利用有限元前处理软件,进行体网格划分,并生成正四面体网格即为多孔晶胞单元,将多孔晶胞单元编号,附属节点编号,并将节点坐标导出;
3.7)将赋有材料属性的患者下颌骨实体网格模型和实心的人工髁突假体数字模型一起导入有限元仿真软件中,利用有限元仿真软件的拓扑优化功能进行优化,得每个单元的相对单元密度;将拓扑优化体单元编号、附属节点、节点坐标,和节点相对单元密度值导出;
3.8)将多孔结构晶胞模型和拓扑优化体单元导出的数据在Excel中进行处理,把相对单元密度数据映射到多孔晶胞单元上,找出拓扑优化模型中离晶胞节点最近的节点,将其密度值赋予此晶胞节点,最终得到所有晶胞节点的密度值;
3.9)根据多孔晶胞单元四个节点的密度值计算晶胞的密度,赋予每个多孔晶胞单元密度数据;
3.10)根据多孔晶胞单元密度,计算出晶胞单元内六根杆的杆径,最终杆径取其平均值;
3.11)已知多孔晶胞单元所有节点的坐标及杆径,利用ug二次开发平台,导入软件成多孔髁突头模型;
3.12)将多孔髁突头模型与优化固定单元通过布尔操作组合成个性化髁突假体;
4.)对个性化髁突假体进行3D打印;
5.)对打印后的个性化髁突假体进行后处理,得到应用于临床的具有拓扑优化固定单元和多孔髁突头单元的个性化髁突假体。
2.如权利要求1所述的具有拓扑优化固定单元和多孔髁突头单元的个性化髁突假体设计方法,其特征在于:所述步骤1.1)中,髁突形状如果基本完整直接重建,或者一边健康可用镜像生成,或者都不健康,则可以根据关节窝形状及下颌骨形状以及健康人体数据库,重新设计一个人工关节外形,然后利用网格处理软件生成四面体实体网格模型。
3.一种如权利要求1所述的具有拓扑优化固定单元和多孔髁突头单元的个性化髁突假体设计方法构建的个性化髁突假体,其特征在于:所述个性化髁突假体包括固定单元和多孔髁突头单元,所述固定单元为具有拓扑优化结构的固定板,所述多孔髁突头单元设置在固定板的上部上,所述多孔髁突头单元为基于密度分布晶胞多孔结构的髁突头,该髁突头包括外层皮质骨部分和内层松质骨部分,每层都是由多个多孔晶胞单元构成,每个多孔晶胞单元均为一个由六根杆组成的正四面体结构,内层松质骨部分的多孔晶胞单元的密度小于外层皮质骨部分的多孔晶胞单元的密度,外层皮质骨部分的多孔晶胞单元的杆长度小于内层松质骨部分的多孔晶胞单元的杆长度;所述固定板的下部上设有螺钉孔。
4.如权利要求3所述的个性化髁突假体,其特征在于:所述固定板为倒V型,所述螺钉孔设置有四个,两两一组,左右对称的布置在V型的左右两侧上。
CN201910299195.7A 2019-04-15 2019-04-15 一种具有拓扑优化固定单元和多孔髁突头单元的个性化髁突假体设计方法及个性化髁突假体 Active CN110236741B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910299195.7A CN110236741B (zh) 2019-04-15 2019-04-15 一种具有拓扑优化固定单元和多孔髁突头单元的个性化髁突假体设计方法及个性化髁突假体

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910299195.7A CN110236741B (zh) 2019-04-15 2019-04-15 一种具有拓扑优化固定单元和多孔髁突头单元的个性化髁突假体设计方法及个性化髁突假体

Publications (2)

Publication Number Publication Date
CN110236741A true CN110236741A (zh) 2019-09-17
CN110236741B CN110236741B (zh) 2024-03-22

Family

ID=67883176

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910299195.7A Active CN110236741B (zh) 2019-04-15 2019-04-15 一种具有拓扑优化固定单元和多孔髁突头单元的个性化髁突假体设计方法及个性化髁突假体

Country Status (1)

Country Link
CN (1) CN110236741B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110750599A (zh) * 2019-09-20 2020-02-04 中国电子科技集团公司第二十八研究所 一种基于实体建模的关联信息抽取和显示方法
CN110811929A (zh) * 2019-11-07 2020-02-21 西安交通大学 一种个性化下肢长骨节段替代物及其制备方法
CN111710022A (zh) * 2020-06-08 2020-09-25 嘉兴学院 一种避免轮廓相交的快速成型切片处理方法
CN112233801A (zh) * 2020-12-17 2021-01-15 季华实验室 内置假体拓扑优化数学模型构建方法及拓扑优化设计方法
CN112712594A (zh) * 2020-12-29 2021-04-27 深圳技术大学 骨骼模型构建方法、装置、计算机设备及存储介质
CN113768666A (zh) * 2021-09-08 2021-12-10 佳木斯大学 一种3d打印基于tpms设计的仿生下颌骨支架的制备方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2000062718A1 (de) * 1999-04-15 2000-10-26 Franz Sutter Gelenkkopfprothese und bausatz zur bildung einer solchen
CN103584931A (zh) * 2013-10-23 2014-02-19 华南理工大学 一种仿生梯度膝关节股骨假体结构及其制造方法
CN104462723A (zh) * 2014-12-25 2015-03-25 北京航空航天大学 一种基于拓扑优化及骨重建仿真的个性化椎间融合器设计方法
US20170014169A1 (en) * 2014-03-11 2017-01-19 Ohio State Innovation Foundation Methods, devices, and manufacture of the devices for musculoskeletal reconstructive surgery
CN109009572A (zh) * 2018-07-12 2018-12-18 广东省新材料研究所 髋关节假体柄及其制造方法
CN109091273A (zh) * 2018-06-13 2018-12-28 西安交通大学 一种基于拓扑结构优化的个性化盆骨内置假体设计方法
CN109152610A (zh) * 2016-05-16 2019-01-04 思想外科有限公司 植入设计和计算机辅助手术
CN210354993U (zh) * 2019-04-15 2020-04-21 浙江工业大学 一种具有拓扑优化固定单元和多孔髁突头单元的个性化髁突假体

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2000062718A1 (de) * 1999-04-15 2000-10-26 Franz Sutter Gelenkkopfprothese und bausatz zur bildung einer solchen
CN103584931A (zh) * 2013-10-23 2014-02-19 华南理工大学 一种仿生梯度膝关节股骨假体结构及其制造方法
US20170014169A1 (en) * 2014-03-11 2017-01-19 Ohio State Innovation Foundation Methods, devices, and manufacture of the devices for musculoskeletal reconstructive surgery
CN104462723A (zh) * 2014-12-25 2015-03-25 北京航空航天大学 一种基于拓扑优化及骨重建仿真的个性化椎间融合器设计方法
CN109152610A (zh) * 2016-05-16 2019-01-04 思想外科有限公司 植入设计和计算机辅助手术
CN109091273A (zh) * 2018-06-13 2018-12-28 西安交通大学 一种基于拓扑结构优化的个性化盆骨内置假体设计方法
CN109009572A (zh) * 2018-07-12 2018-12-18 广东省新材料研究所 髋关节假体柄及其制造方法
CN210354993U (zh) * 2019-04-15 2020-04-21 浙江工业大学 一种具有拓扑优化固定单元和多孔髁突头单元的个性化髁突假体

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110750599A (zh) * 2019-09-20 2020-02-04 中国电子科技集团公司第二十八研究所 一种基于实体建模的关联信息抽取和显示方法
CN110750599B (zh) * 2019-09-20 2022-06-28 中国电子科技集团公司第二十八研究所 一种基于实体建模的关联信息抽取和显示方法
CN110811929A (zh) * 2019-11-07 2020-02-21 西安交通大学 一种个性化下肢长骨节段替代物及其制备方法
CN111710022A (zh) * 2020-06-08 2020-09-25 嘉兴学院 一种避免轮廓相交的快速成型切片处理方法
CN111710022B (zh) * 2020-06-08 2021-03-23 嘉兴学院 一种避免轮廓相交的快速成型切片处理方法
CN112233801A (zh) * 2020-12-17 2021-01-15 季华实验室 内置假体拓扑优化数学模型构建方法及拓扑优化设计方法
CN112233801B (zh) * 2020-12-17 2021-03-26 季华实验室 内置假体拓扑优化数学模型构建方法及拓扑优化设计方法
CN112712594A (zh) * 2020-12-29 2021-04-27 深圳技术大学 骨骼模型构建方法、装置、计算机设备及存储介质
CN113768666A (zh) * 2021-09-08 2021-12-10 佳木斯大学 一种3d打印基于tpms设计的仿生下颌骨支架的制备方法
CN113768666B (zh) * 2021-09-08 2024-01-02 佳木斯大学 一种3d打印基于tpms设计的仿生下颌骨支架的制备方法

Also Published As

Publication number Publication date
CN110236741B (zh) 2024-03-22

Similar Documents

Publication Publication Date Title
CN110236741A (zh) 一种具有拓扑优化固定单元和多孔髁突头单元的个性化髁突假体设计方法及个性化髁突假体
CN109091273B (zh) 一种基于拓扑结构优化的个性化盆骨内置假体设计方法
Luo et al. Finite-element design and optimization of a three-dimensional tetrahedral porous titanium scaffold for the reconstruction of mandibular defects
CN105232169B (zh) 多段式种植体组件及其制备方法
CN108635084B (zh) 熔融沉积3d打印制备的聚醚醚酮变模量人工骨骼替代物及制备方法
CN110929379B (zh) 一种拓扑优化人工椎体及其设计方法
CN110613533A (zh) 用于下颌骨体箱状缺损修复的pekk个性化植入体设计制作方法及植体
CN109172049A (zh) 一种基于分层片状杆连接的多孔网状结构骨科修复植入体的设计制作方法及植入体
CN107607452A (zh) 一种按应力分布确定多孔假体孔隙率的方法
Peng et al. Biomechanical and Mechanostat analysis of a titanium layered porous implant for mandibular reconstruction: The effect of the topology optimization design
CN107680163B (zh) 一种个性化3d优化结构内固定板的制作方法及内固定板
CN106580520A (zh) 具有pekk支撑固定单元与组织工程生长单元的下颌骨植入体制作方法及植入体
CN112107394A (zh) 一种用于颌面骨缺损修复的植入体及其制备方法
CN109009581A (zh) 一种组配式人工椎体
Dai et al. Design of a maxillofacial prosthesis based on topology optimization
CN111986310B (zh) 一种考虑软组织附着和骨生长的下颌骨多孔植入体的设计方法及多孔植入体
CN112233801B (zh) 内置假体拓扑优化数学模型构建方法及拓扑优化设计方法
CN211271408U (zh) 用于下颌骨体箱状缺损修复的pekk个性化植入体
Liu et al. Novel design and optimization of porous titanium structure for mandibular reconstruction
CN103315831B (zh) 一种具有高活动度的一体化人工颈椎间盘假体
Arabshahi et al. Influence of thickness and contact surface geometry of condylar stem of TMJ implant on its stability
CN210354993U (zh) 一种具有拓扑优化固定单元和多孔髁突头单元的个性化髁突假体
Cheng et al. 3D-printed porous condylar prosthesis for temporomandibular joint replacement: Design and biomechanical analysis
CN205964230U (zh) 一种个性化甲状软骨假体
CN209059552U (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
GR01 Patent grant
GR01 Patent grant