一种仿生骨支架的等参变换混合结构及其3D打印方法

摘要

本发明属于增材制造及仿生结构设计领域，并公开了一种仿生骨支架的等参变换混合结构及其3D打印方法，所述仿生骨支架分为内层、过渡层和外层，所述外层采用主拉伸变形的极小曲面模拟密质骨，所述内层采用比表面积高的极小曲面模拟松质骨，所述过渡层通过Sigmond函数对内层和外层两种结构进行结合，且整体结构通过等参变换模仿人体骨骼的外形。该仿生骨支架能模仿人体骨骼的外形、孔隙特征、杨氏模量及完全连通的孔隙结构，能够与骨骼受损区域相匹配，能够实现体液运输，并可满足力学性能和骨长入要求，可用于骨缺损的治疗。

说明书

技术领域

本发明涉及增材制造和仿生结构设计领域，尤其涉及一种仿生骨支架的等参变换混合结构及其3D打印方法。

背景技术

骨骼是人体的重要组成部分，对人体起支撑和保护的作用，然而每年都会有大量因创伤、开颅手术、骨肿瘤切除、骨畸形矫正等情况造成骨缺损的病患，需通过骨移植进行修复。但医用骨源紧张，面对每年大量的骨缺损患者，难以满足医疗需求。为了解决骨源不足的问题，有许多研究人员希望通过在体外构建合适的仿生支架，来填补骨骼缺损的部位，从而解决骨源不足的问题。

为了保证骨细胞在仿生支架中能够有较高的细胞活性，在孔隙结构方面，仿生支架应为骨细胞的增殖和黏附提供良好的活动空间，并为营养物质和代谢废物的运输提供良好的运输通道，因此需要考虑较大的比表面积和良好的孔道连通率满足此要求；在力学性能方面，仿生支架应考虑具有与人体骨有相似的杨氏模量，可以防止应力屏蔽现象的出现，避免骨组织损伤和仿生支架的过早失效；在材料方面，仿生支架应采用生物相容性材料，不仅仅要求材料的无毒性，也要求材料对骨细胞活动具有积极的影响。因此，仿生支架的设计需要考虑到孔隙结构、力学性能和材料等方面。而3D打印较传统制造方式，具有设计形状无约束、无模具、小规模定制等优势，3D打印技术让骨支架复杂结构制造成为可能，在骨支架制造中具有极大优势。

专利CN112316207A采用的混合点阵多孔仿生支架，和专利CN110272273A所构建的极小曲面P-cell、S-14、G-yroid和I-wp结构，均是通过对极小曲面直接进行加厚，然后直接将模型切除成圆柱等其它形状获得的，只能用于替换人体密质骨，这种直接对极小曲面加厚所获得的多孔结构被称之为Sheet型极小曲面结构，如图11极小曲面类似一块扭曲的薄板，会将完整的空间分割为两个互不相通的第一子空间81和第二子空间82，其空间构建方式是不完全连通的，加厚即是将薄板8进行延伸，因此会影响体液的运输，和新生骨组织所需营养物质的提供，进而影响骨生长的效果，而直接切除极小曲面会破坏结构单元的完整性。

专利CN107661160A所采用的多孔结构为直杆结构，孔隙虽然完全连通，但其骨修复支架内部结构设计为空腔，在受力时，会出现应力集中，降低了仿生支架的抗压强度和疲劳寿命，不适合于本发明专利中所述骨支架的应用。现有关于仿生支架的研究不能同时满足孔隙结构、力学性能和材料三个方面的要求。

发明内容

本发明针对现有技术的不足之处，提供了一种仿生骨支架的等参变换混合结构设计及其3D打印方法。该仿生骨支架能模仿人体骨骼的外形、孔隙特征、杨氏模量，拥有完全连通的孔隙结构，且能够与骨骼受损区域相匹配，并可满足力学性能和骨长入要求。同时采用生物相容性PEEK材料制造仿生骨支架，可用于骨缺损的治疗。

本实发明为实现上述目标所采用的技术方案是：

一种仿生骨支架的等参变换混合结构的3D打印方法，其特征在于，包括以下步骤：

步骤一：利用内层隐式函数、外层隐式函数分别生成Network型的内层极小曲面结构胞元、外层极小曲面结构胞元的stl模型，结构胞元为边长为

内层结构隐式函数的表达式包括：

Gyroid：

或Diamond：

外层结构隐式函数的表达式包括：

Primitiv-Opt：

或I-Wrapped Package：

其中：

步骤二：将正方体模型从内到外划分为共中心轴的三层，内层为N棱柱柱体、过渡层为内层外的N边形套筒、外层为正方体模型的剩余部分，N棱柱柱体，再将划分后的正方体模型从上至下等分为厚度为

通过二维或三维线性拉格朗日形函数确定步骤一中包含内层和外层极小曲面结构胞元的虚拟正立方体与第一六面体子空间和第三六面体子空间的拓扑关系，进行等参变换，使虚拟正立方体中的内层和外层极小曲面结构胞元内的物质点分别依次映射到所有划分出的第一六面体子空间和第三六面体子空间中；

通过Sigmond函数设计过渡层的stl模型，对内层和外层的两种胞元结构进行融合，再将融合后的胞元结构映射到第二六面体子空间，并使过渡层的孔隙率

其中，

最后完成第一六面体子空间、第二六面体子空间、第三六面体子空间中极小曲面曲面结构的合并、平滑，获得正方体状的仿生支架，再通过三维二次拉格朗日形函数确定正方体状的仿生支架的各角点及各边中点分别与圆柱状的骨空间上下端面的各四个等分点及各等分点之间连线上中点的拓扑关系，进行等参变换，使正方体状的仿生支架整体映射变形为圆柱状的骨空间，得到仿生骨支架三维模型，再生成仿生骨支架三维模型的stl文件；

步骤三：将仿生骨支架三维模型导入到切片软件中进行切片，导出gcode指令文件，并将gcode指令文件导入熔融沉积成型3D打印设备中；

步骤四：通过熔融沉积成型设备完成对仿生骨支架的制造。

进一步地，所述的二维和三维线性拉格朗日形函数为：

二维线性：

三维线性：

其中：

进一步地，所述的三维二次拉格朗日形函数的表达式为：

在角点上

在边中点上

其中：

进一步地，熔融沉积成型设备所使用的材料为聚醚醚酮。

进一步地，步骤四还包括将熔融沉积成型设备的参数设置为：喷头温度为420℃，打印速度为25mm/s，打印层厚为0.05mm；经步骤四完成打印后，骨支架在120℃的温度下干燥4小时，再以10℃/min的速率加热至200℃，在环境下保温4小时后，随炉冷却。

进一步地，

进一步地，所述的方法制备的仿生骨支架，所述仿生骨支架整体结构是所述正方体状仿生支架通过三维二次拉格朗日形函数等参变换映射到骨空间中得到模仿人体骨骼的外形，分为内层结构、过渡层和外层结构，所述内层结构的孔隙率在50%到90%之间，孔径大小在600um到1500um之间，杨氏模量在0.02GPa到0.8GPa之间，用于模拟松质骨；所述外层结构的孔隙率在5%到20%之间，孔径大小在100um到700um之间，杨氏模量在3GPa到4GPa之间，用于模拟密质骨；所述过渡层的孔隙率和孔隙大小随径向线性变化。

进一步地，所述仿生骨支架整体结构是所述正方体状仿生支架通过三维二次拉格朗日形函数等参变换映射到骨空间中得到模仿人体骨骼的外形，分为内层结构、过渡层和外层结构，所述正方体状仿生支架分为内层、过渡层和外层，所述内层被分割为步骤二中映射后的多个第一六面体子空间，外层被分割为步骤二中映射后的多个第三六面体子空间，所述过渡层通过步骤二中的Sigmond函数对内层和外层两种结构进行融合，第一六面体子空间中填充步骤二中等参变换后的内层极小曲面结构胞元，第三六面体子空间中填充步骤二中等参变换后的外层极小曲面结构胞元，内层极小曲面结构胞元的极小曲面多孔结构的类型为主弯曲变形的Network型Gyroid结构或Network型Diamond结构，外层极小曲面结构胞元的极小曲面多孔结构的类型为主拉伸变形的Network型I-Wrapped Package型或结构Network型Primitive-Opt型结构；内层结构的孔隙率在50%到90%之间，孔隙大小600-1500um，杨氏模量在0.3GPa到0.5GPa之间，屈服强度在5MPa到15MPa之间；外层结构的孔隙率为在5%到20%之间，孔隙大小300-700 um，杨氏模量在3GPa到4GPa之间, 屈服强度在60MPa到90MPa之间。

进一步地，内层结构常采用的极小曲面多孔结构：Network型Gyroid结构和Network型Diamond结构，具有高渗透率和比表面积，较外层结构，更适于模拟松质骨。

所述外层结构采用的极小曲面多孔结构类型为主拉伸变形Network型Primitive-Opt结构和Network型I-Wrapped Package型结构，较内层结构，具有很高的抗压强度，适于模拟密质骨。

通过本发明以上介绍的技术方案，与现有技术相比具有以下优势：

1、本发明中的内层仿生支架由六面体胞元阵列组成，对比文件CN112316207和专利CN110272273A的内层仿生支架是整体增材形成的，本发明的每个六面体胞元都进行了拉格朗日变化，该方法的优势在于可以使规则的结构按所设计的不规则形状进行变形，同时保持骨支架结构完整，防止力学性能下降。本发明中骨支架外层选用主拉伸变形的Network型极小曲面结构，此类结构有具有较高的抗压强度（Primitive-Opt、I-WP等）；内层采用的是比表面积高、渗透率高的Network型极小曲面结构，压缩变形方式一般为主弯曲变形（Gyroid、Diamond）。

2、本发明采用等参变换对极小曲面结构胞元进行映射，可以保证每个极小曲面结构胞元的完整性，且可以根据人体骨骼外形对仿生支架的外形进行调整，以匹配骨缺损区域的形状。

3、本发明所设计的仿生支架可以同时满足孔隙结构、力学性能和材料三个方面的要求。采用填充的方式构建极小曲面多孔结构，从而达到孔隙完全连通的效果，并保证孔隙大小和孔隙率与人体骨骼相似，满足孔隙结构的要求。使用生物相容性PEEK材料，从而保证仿生支架具有与人体骨骼相似的杨氏模量和较高的抗压强度，满足力学性能和材料的要求。因此，在保证仿生支架内层和外层孔隙连通及孔隙大小、孔隙率与人体松质骨（支架内层）和密质骨（支架外层）相匹配的情况下，仍能与人体松质骨和密质骨的弹性模量相匹配，防止“应力屏蔽”，可用于替代受损的人体松质骨和密质骨。

4、本发明所设计的仿生支架能与人体骨骼的梯度结构相匹配，可通过调整内层和外层的区域大小来匹配人体松质骨和密质骨厚度，以使骨支架能够提供匹配的梯度结构。

5、本发明采用材料PEEK制造仿生骨支架，该材料有着良好的生物相容性，并与密质骨有相似的杨氏模量，且较之钛合金，具有更小的密度，可以减轻仿生骨支架对人体的负担。

附图说明

下面结合附图和实施例对本发明进一步说明。

图1为实施例1和实施例2中仿生支架填补股骨缺损区域的示意图；

图2为实施例1和实施例2中仿生支架内层、过渡层和外层区域划分示意图；

图3为实施例1和实施例2中仿生支架内层、过渡层和外层子空间划分示意图，将支架整体分割为一系列子空间；

图4为实施例1中基于Gyroid和I-WP曲面构建的Network型结构胞元，分别用于内层和外层子空间填充；

图5为实施例1中Network型Gyroid结构胞元等参变换前后形状变化示意图，用于子空间填充；

图6为实施例1中仿生支架等轴视图和俯视图；

图7为实施例1中仿生支架等参变换后的等轴视图和俯视图；

图8为二维线性等参变换节点映射关系图，用于子空间填充，通过控制节点来控制等参变换后的结构形状；

图9为三维线性等参变换节点映射关系图，用于子空间填充，通过控制节点来控制等参变换后的结构形状；

图10为三维二次等参变换节点映射关系图，用于仿生骨支架整体的形状映射，通过控制节点来控制等参变换后的结构形状，以匹配人体骨骼外形；

图11为将Gyroid极小曲面通过填充的方式构建的Network型Gyroid结构胞元及通过加厚的构建方式构建的Sheet型Gyroid结构胞元；

图12为实施例1中仿生支架内层Gyroid结构流体仿真压力云图；

图13为实施例1中仿生支架外层I-WP结构流体仿真压力云图；

图14为实施例1中仿生支架外层I-WP结构静力学仿真等效应力云图；

图15为实施例2中基于Diamond和Primitive-Opt曲面构建的Network型结构胞元，左侧极小曲面结构用于内层空间填充，右侧极小曲面结构用于外层子空间填充；

图16为实施例2中仿生支架内层Diamond结构流体仿真压力云图；

图17为实施例2中仿生支架外层Primitive-Opt结构流体仿真压力云图；

图18为实施例2中仿生支架外层Primitive-Opt结构静力学仿真等效应力云图；

图19为实施例2中仿生支架等轴视图和俯视图；

图20为实施例2中仿生支架等参变换后的等轴视图和俯视图；

所有附图中相同的元件或结构用相同的标记号来表示，其中：

图中： 1、内层区域；2、过渡层区域；3、外层区域；4、子空间；5、Network型Gyroid结构胞元；6、Network型I-WP结构胞元；7、等参变换后的Network型Gyroid结构胞元；8、Sheet型Gyroid结构胞元；9、Gyroid极小曲面胞元；81：Sheet型Gyroid结构胞元的第一子空间；82：Sheet型Gyroid结构胞元的第二子空间；10、Network型Diamond结构胞元；11、Network型Primitive-Opt结构胞元。

具体实施方式

为了使本发明的发明目的、技术方案和发明优势更加清楚阐述，以下将结合说明书附图对本发明做进一步详细讲解。以下所描述的实施例仅为本发明基于等参变换的混合结构仿生骨支架设计方法所设计的一种实施例，不包含全部实施例。且该仿生骨支架设计方法不仅限于股骨仿生支架，通过如图8、图9和图10所示的等参变换可以根据人体骨骼特征设计出其它非柱状仿生骨支架。本发明采用的是Network型极小曲面结构，其是通过填充被极小曲面分割的两个互不相通的子空间的其中一个来构建，即只保留一个子空间，这样构建的点阵结构则具有完全连通的孔隙结构，有利于体液的运输。Sheet型和Network型Gyroid极小曲面结构胞元如图11所示。

首先对本发明所涉及的术语解释如下：

“应力屏蔽”：目前的金属医用植入材料普遍拥有比自然骨更高的杨氏模量，当植入到人体后，载荷大部分被植入材料所承担，使得植入材料周围的骨得不到足够的力刺激，致使骨吸收作用远大于骨形成，骨骼变得疏松并因此引起植入材料松动，这种导致骨质疏松的生物力学不相容被称为应力屏蔽效应。因此用于骨修复的金属生物医用材料应具有与骨骼相当的模量，根据人体自然骨的类型和测量方向的不同，其模量的大小从4GPa到30GPa不等。

实施例1：

如图1所示，为本实施例所设计的一种仿生骨支架，其能很好匹配股骨骨干外形和力学性能。该仿生骨支架构建的大致过程为：

首先，生成一个边长为25mm的正方体模型。然后，将正方体模型通过图2和图3所示的划分方法被分割为一系列子空间4。接着，将空间为正立方体的极小曲面结构胞元5按照子空间4的空间节点参数进行如图5所示的线性等参变换，获得变换后极小曲面结构胞元7，将变换后极小曲面结构胞元7填充到所有分割出的子空间4中，获得如图6所示的仿生骨支架。最后，通过二次等参变换使仿生骨支架的外形与人体骨骼相匹配，完成仿生骨支架的构建。

下面对本实施例中，仿生骨支架的构建过程进行详细说明。

如图2所示，该仿生骨支架从内到外被划分为三层，分别为内层1、过渡层2和外层3。其中内层1直径为15mm，孔隙率为70%，模拟松质骨；过渡层2内径为15mm，外径为20mm，孔隙率从70%线性下降至10%；外层3内径为20mm，外部边长为25mm，孔隙率为10%，模拟密致骨。内层1结构采用Network型Gyroid结构5，外层3采用Network型I-WP结构6，内层采用过渡层2采用Sigmond函数对内层1和外层3两种结构进行融合，融合边界线为图2中的虚线圆，虚线圆半径为8.75mm。

上述其中两种极小曲面结构和利用Sigmond函数进行结构融合的公式如下：

Network型Gyroid曲面结构：

Network型I-WP曲面结构：

(式2)

融合结构公式：

其中：

通过上述公式，可以获得仿生骨支架内层、过渡层和外层的极小曲面结构的隐式函数。所有曲面模型都是利用Matlab2019软件通过极小曲面隐式函数生成的，在获得内层、过渡层和外层的极小曲面胞元stl模型后，导入Geomagic Wrap2016进修复、平滑和填充，再重新导出修复后的极小曲面结构胞元stl模型，为后续子空间4填充作准备，其中结构胞元为边长为

对于极小曲面函数来说，由曲面

如图3所示，在本实施例中形状为正方体模型的仿生骨支架被平均划分为10层，总高为25mm。在每一层中，内层1、过渡层2和外层3根据图2被分割为一系列六面体子空间4，通过将上述的极小曲面结构胞元（内层1结构采用Network型Gyroid极小曲面结构5，外层3采用Network型I-WP极小曲面结构6，过渡层2采用融合结构胞元)填充到所有子空间4中来完成仿生骨支架结构的构建。

如图5所示，为了使虚拟正立方体内规则的极小曲面结构单元5填充到不规则的子空间4中，通过二维或三维线性拉格朗日形函数对极小曲面结构胞元进行等参变换，从而将包含极小曲面结构胞元的虚拟正立方体映射到子空间4中。二维和三维线性拉格朗日形函数为：

二维线性：

三维线性：

其中：

其中，

如图6所示，通过等参变换将极小曲面结构胞元依次映射到所有子空间后，使用Geomagic Wrap2016软件对子空间中的结构进行合并、平滑，得到如图6所示的正方体状的仿生骨支架。最后通过三维二次拉格朗日形函数对图6所示的正方体状的仿生骨支架整体进行映射变形为圆柱状的骨空间，三维二次拉格朗日形函数的表达式为：

在角点上

在边中点上

三维二次等参变换节点映射关系如图10所示，正方体状的仿生支架的各角点及各边中点分别与圆柱状的骨空间上下端面的各四个等分点及各等分点之间连线上中点有一一对应的拓扑关系。

如图7所示为映射后的仿生骨支架，其不仅可以很好地匹配股骨骨干的外形，而且仿生骨支架外层的厚度为2.5mm，孔隙率为10%，孔隙大小在200um-700um之间，符合密质骨的孔隙结构特征。仿生骨支架内层的孔隙率为70%，孔隙大小大多数在800um-1200um之间，少数在300um左右，符合松质骨的孔隙结构特征。并且所有孔隙完全连通，有利于体液的运输。

为进一步了解本实施例所设计的仿生骨支架内层和外层的力学性能，对内层和外层结构分别进行力学仿真。力学仿真采用Abaqus2020有限元分析软件，具体步骤为：首先，分别建立孔隙率为70%的Network型Gyroid极小曲面结构和Network型孔隙率为10%的I-WP极小曲面结构的5x5x5胞元阵列模型，每个胞元边长为2.5mm，模型总边长为12.5mm。其次，将两种模型导入Hypermesh划分四面体网格，网格平均大小为0.175mm，将网格模型inp文件导入Abaqus2020。然后，在Abaqus2020中，对网格模型底面施加固定约束，上表面使用参考点进行耦合，并在参考点上向下施加3300N的力（根据成人体重75kg的5倍和密质骨区域横截面积换算），网格类型设置为线性四面体单元(C3D4)，材料使用PEEK，仿生骨支架的杨氏模量为3.738GPa，泊松比为0.39，屈服强度为99MPa。本实施例外层支撑结构力学仿真应力云图如图14所示，最大应力为77.05MPa，小于材料屈服强度，满足使用要求。

为进一步了解本实施例所设计的仿生骨支架内层和外层的渗透性能，对内层和外层结构分别进行流体仿真。流体仿真采用COMSOL Multiphysics 5.5有限元分析软件，具体步骤为：首先，利用布尔运算生成上述极小曲面结构的孔道模型。其次，将孔道模型导入Hypermesh划分网格，导出bdf格式的孔道网格模型，并导入COMSOL Multiphysics 5.5。然后，在COMSOL Multiphysics 5.5中，为模拟体液，物理场设置为层流，研究类型设置为稳态，并将孔道模型上表面设置为入口，入口速度为0.001m/s，下表面设置为出口，出口压力为0Pa，其它表面设置为非滑移墙壁，材料选用软件材料库中的参数H

表1 实施例仿生骨支架与人体骨骼参数表

表中引用的参考文献：

[1]Sanjairaj Vijayavenkataraman, Lei Zhang, Shuo Zhang, et al. TriplyPeriodic Minimal Surfaces Sheet Scaffolds for Tissue EngineeringApplications: An Optimization Approach toward Biomimetic Scaffold Design[J].

[2]Chunze Yan, Liang Hao, Ahmed Hussein, Philippe Young. Ti-6Al-4Vtriply periodic minimal surface structures for bone implants fabricated viaselective laser melting[J].

[3] Li Yuan，Songlin Ding，Cuie Wen.Additive manufacturing technologyfor porous metal implant applications and triple minimal surface structures：Areview[J].

由仿真结果可得表1中的各项参数。通过表1可知，实施例1所设计的仿生骨支架的内层和外层结构皆具有与人体骨骼相似的孔隙特征和力学性能，可用于骨缺损的治疗。

实施例2：

进一步地，如图2所示，该仿生骨支架从内到外被划分为三层，分别为内层1、过渡层2和外层3。其中内层1直径为15mm，孔隙率为70%，模拟松质骨；过渡层2内径为15mm，外径为20mm，孔隙率从70%线性下降至10%；外层3内径为20mm，外部边长为25mm，孔隙率为10%，模拟密致骨。如图15内层1结构采用Network型Diamond结构10，外层3采用Network型Primitive-Opt结构11，内层采用过渡层2采用Sigmond函数对内层1和外层3两种结构进行融合，融合边界线为图2中的虚线圆，虚线圆半径为8.75mm。

上述其中两种极小曲面结构和利用Sigmond函数进行结构融合的公式如下：

Network型Diamond曲面结构：

Network型Primitive-Opt曲面结构：

（式5）

融合结构公式：

其中：

通过上述公式，可以获得仿生骨支架内层、过渡层和外层的极小曲面结构的隐式函数。所有曲面模型都是利用Matlab2019软件通过极小曲面隐式函数生成的，在获得内层、过渡层和外层的极小曲面胞元stl模型后，导入Geomagic Wrap2016进修复、平滑和填充，再重新导出修复后的极小曲面结构胞元stl模型，为后续子空间4填充作准备，其中结构胞元为边长为

如图3和图15所示，在本实施例中形状为正方体模型的仿生骨支架被平均划分为10层，总高为25mm。在每一层中，内层1、过渡层2和外层3根据图2被分割为一系列六面体子空间4，通过将上述的极小曲面结构胞元（内层1结构采用Network型Diamond极小曲面结构10，外层3采用Network型Primitive-Opt极小曲面结构11，过渡层2采用融合结构胞元)填充到所有子空间4中来完成仿生骨支架结构的构建。

如图5所示，为了使虚拟正立方体内规则的极小曲面结构单元填充到不规则的子空间4中，通过二维或三维线性拉格朗日形函数对极小曲面结构胞元进行等参变换，从而将包含极小曲面结构胞元的虚拟正立方体映射到子空间4中。二维和三维线性拉格朗日形函数为：

二维线性：

三维线性：

其中：

其中，

如图19所示，通过等参变换将极小曲面结构胞元依次映射到所有子空间后，使用Geomagic Wrap2016软件对子空间中的结构进行合并、平滑，得到如图19所示的正方体状的仿生骨支架。最后通过三维二次拉格朗日形函数对图19所示的正方体状的仿生骨支架整体进行映射变形为图20所示的圆柱状的骨空间，三维二次拉格朗日形函数的表达式为：

在角点上

在边中点上

三维二次等参变换节点映射关系如图10所示，正方体状的仿生支架的各角点及各边中点分别与圆柱状的骨空间上下端面的各四个等分点及各等分点之间连线上中点有一一对应的拓扑关系。

如图20所示为映射后的仿生骨支架，其不仅可以很好地匹配股骨骨干的外形，而且仿生骨支架外层的厚度为2.5mm，孔隙率为10%，孔隙大小在200um-700um之间，符合密质骨的孔隙结构特征。仿生骨支架内层的孔隙率为70%，孔隙大小大多数在800um-1200um之间，少数在300um左右，符合松质骨的孔隙结构特征。并且所有孔隙完全连通，有利于体液的运输。

为进一步了解本实施例所设计的仿生骨支架内层和外层的力学性能，对内层和外层结构分别进行力学仿真。力学仿真采用Abaqus2020有限元分析软件，具体步骤为：首先，分别建立孔隙率为70%的Network型Diamond极小曲面结构和Network型孔隙率为10%的Primitive-Opt极小曲面结构的5x5x5胞元阵列模型，每个胞元边长为2.5mm，模型总边长为12.5mm。其次，将两种模型导入Hypermesh划分四面体网格，网格平均大小为0.175mm，将网格模型inp文件导入Abaqus2020。然后，在Abaqus2020中，对网格模型底面施加固定约束，上表面使用参考点进行耦合，并在参考点上向下施加3300N的力（根据成人体重75kg的5倍和密质骨区域横截面积换算），网格类型设置为线性四面体单元(C3D4)，材料使用PEEK，仿生骨支架的杨氏模量为3.738GPa，泊松比为0.39，屈服强度为99MPa。本实施例外层支撑结构力学仿真应力云图如图18所示，最大应力为65.4MPa，小于材料屈服强度，满足使用要求。

为进一步了解本实施例所设计的仿生骨支架内层和外层的渗透性能，对内层和外层结构分别进行流体仿真。流体仿真采用COMSOL Multiphysics 5.5有限元分析软件，具体步骤为：首先，利用布尔运算生成上述极小曲面结构的孔道模型。其次，将孔道模型导入Hypermesh划分网格，导出bdf格式的孔道网格模型，并导入COMSOL Multiphysics 5.5。然后，在COMSOL Multiphysics 5.5中，为模拟体液，物理场设置为层流，研究类型设置为稳态，并将孔道模型上表面设置为入口，入口速度为0.001m/s，下表面设置为出口，出口压力为0Pa，其它表面设置为非滑移墙壁，材料选用软件材料库中的参数H

表2 实施例仿生骨支架与人体骨骼参数表

表中引用的参考文献：

[1]Sanjairaj Vijayavenkataraman，Lei Zhang，Shuo Zhang，et al. TriplyPeriodic Minimal Surfaces Sheet Scaffolds for Tissue EngineeringApplications：An Optimization Approach toward Biomimetic Scaffold Design[J].

[2]Chunze Yan，Liang Hao，Ahmed Hussein，Philippe Young. Ti–6Al–4Vtriply periodic minimal surface structures for bone implants fabricated viaselective laser melting[J].

[3] Li Yuan，Songlin Ding，Cuie Wen.Additive manufacturing technologyfor porous metal implant applications and triple minimal surface structures：Areview[J].

由仿真结果可得表2中的各项参数。通过表2可知，实施例2所设计的仿生骨支架的内层和外层结构皆具有与人体骨骼相似的孔隙特征和力学性能，可用于骨缺损的治疗。

而现有技术CN112206077A提供了混合点阵多孔梯度支架结构，目的在于提供渗透性接近疏松骨的多孔梯度支架替代疏松骨，在孔隙率为70%时（略大于本实施例时），其弹性模量（或称杨氏模量）大于4GPa，与密质骨相匹配，其弹性模量不与松质骨相匹配，弹性模量大于0.8GPa，会造成松质骨的“应力屏蔽”。现有技术CN112316207A公布的混合点阵多孔梯度支架结构在孔隙率为65%时（接近本实施例时），其弹性模量（或称杨氏模量）大于5GPa，与密质骨相匹配，其弹性模量不与松质骨相匹配，弹性模量大于0.8GPa，会造成松质骨的“应力屏蔽”。因此以上两现有技术的孔隙率调整为50~90%，处于人体松质骨孔隙率范围时，其弹性模量远大于人体松质骨的弹性模量，以上两现有技术的孔隙率调整为30%，处于人体密质骨孔隙率范围时，其弹性模量接近于人体密质骨的弹性模量，现有技术提供的混合点阵多孔梯度支架结构仅能替换受损的密质骨，为避免“应力屏蔽”不能用于替换受损的松质骨。

本实施例的仿生骨支架采用填充的方式构建的内层结构具有完全连通的孔隙结构，空间不会被极小曲面分割为互不连通的两块区域，同时内层结构由各个独立的极小曲面结构元胞合并形成，既保证了孔隙结构在整体的内层结构中的均匀分布，又保证了整体的内层结构的屈服强度。仿生骨支架的外层和内层的厚度可以根据人体不同部位骨骼的密质骨和松质骨厚度进行调节，以匹配骨骼结构特征。仿生骨支架的外层和内层在保证孔隙大小、孔隙连通、孔隙率与人体密质骨（支架外层）和松质骨（支架内层）相匹配的情况下，仍能满足人体松质骨和密质骨的弹性模量要求，防止“应力屏蔽”，因此可替换受损的密质骨及松质骨。

以上实施例基于等参变换的混合结构仿生骨支架的制造方法，包括以下步骤：

步骤一：利用Matlab软件通过所设计的隐式函数：式1和式4、式2和式5分别生成内层极小曲面结构胞元、外层极小曲面结构胞元的stl模型，式1和式2中

步骤二：将正方体模型从内到外划分为共中心轴的三层，内层1为N棱柱柱体、过渡层2为内层1外的N边形套筒、外层3为正方体模型的剩余部分，再将划分后的正方体模型从上至下等分为厚度为2.5mm的M层模型切片，每层模型切片中的内层1分割为多个第一六面体子空间，每层模型切片中的过渡层2、外层3均分割为N个第二六面体子空间、第三六面体子空间，且每层模型切片中：过渡层的外径与内径的差值为

步骤三：将仿生骨支架三维模型导入到切片软件中进行切片，导出gcode指令文件，并将gcode指令文件导入熔融沉积成型3D打印设备中；

步骤四：通过熔融沉积成型设备使用PEEK材料对仿生支骨架三维模型进行制造，其中喷头温度为420℃，打印速度为25mm/s，打印层厚为0.05mm；经步骤四完成打印后，骨支架在120℃的温度下干燥4小时，再以10℃/min的速率加热至200℃，在环境下保温4小时后，随炉冷却。

通过以上步骤即可完成对本实施例中仿生支架三维模型的制造。

以上说明描述本发明中的一个较佳的实施方式，不应将其看作为是对本发明权利要求保护范围的限制。在不脱离本发明原理和精神的情况下，任何修改、等效替换及改进，都应视为在本发明权利要求保护范围之内。