[发明专利]基于四面体网格的仿周期极小曲面的内部孔隙设计方法

说明书

一种基于四面体网格的仿周期极小曲面的内部孔隙设计方法，采用的具体步骤如下：一：提取图像；二：生成四面体网格；三：进行切片处理；四：选定一层切片平面；五：得到特定方向的扫描线；六：找出所有扫描线起点到终点间的四面体T；七：计算与每个扫描线相交的四面体T相交起点A和相交终点B；八：利用四面体T的四个顶点，建立四元隐式函数；九：对采样点函数值与设定值C进行比较，大于设定值C则跳过，小于则进入步骤十；十：将采样点函数值计入3D打印路径；十一：判断切片平面是否处理完毕，是则输出3D打印路径，否则，选择另一切片平面，进入步骤五。本方法弥补了现有的三周期极小曲面技术的不足。

技术领域

本发明涉及生物材料技术领域，具体涉及一种基于四面体网格的仿周期极小曲面的内部孔隙设计方法。

背景技术

生物医学组织工程作为一项新兴医学技术，已逐步从实验室走向实际医疗应用中。其基本原理是在人体外利用一种生物相容性好、结构合理的3D支架(scaffolds)来培养正常的人体组织细胞，然后把培养好的细胞与支架复合体植入人体组织或器官的病损部位，与人体自身的组织细胞在体内自然融合，以达到修复创伤的目的。

生物医学组织工程中的一个重要环节是设计并制造3D支架。这类支架通常需要有复杂且充分连通的内部孔隙结构以供细胞在其中生长并形成功能组织。传统的减材制造方法在处理具有复杂内部结构的几何形体时困难重重。最近几年逐渐流行起来的3D打印(或增材制造技术)提供了一个可行的制造复杂内部结构支架的方法，但对应的支架设计方法却大多集中在三周期极小曲面(Triply Periodic Minimal Surfaces或TPMS)技术。这种内部孔隙生成技术是基于单个长方体定义域，设计简单灵活但存在着几个明显不足：1、当处理复杂表面形状的部件时，通常我们要对部件进行三维实体剖分。相较于目前成熟的四面体网格划分及优化算法和软件，六面体网格剖分方法降低了灵活性且更难规范化。2、TPMS定义在单个长方体定义域，而六面体网格中的六面体可能是任意的。从长方体到任意六面体的TPMS函数映射可能存在很大的形状扭曲。3、在3D打印中一个关键的步骤是切片算法(slicing)，即计算一个平面与被打印物体表面网格的相交区域。TPMS作为一种隐式函数表示方法，我们并不需要明确定义一个庞大的物体表面网格，而只需计算平面与TPMS对应的每个六面体的相交从而确定要打印的区域。但平面与一般六面体的相交远比平面与四面体的相交情形要复杂。该复杂度不仅降低了算法的可靠性，而且大大增加了计算代价。

发明内容

本发明为弥补现有的TPMS技术的不足，提出了一种利用四面体单元的隐式函数表示方法，大大提高了生成任意曲面形状的内部孔隙结构的适用性，相对于长方体到任意六面体的函数映射，规则四面体到任意四面体的函数映射较为简洁，由此带来的形状扭曲较小。基于四面体网格的四周期极小曲面(Quadruply Periodic Minimal Surfaces或QPMS)的内部孔隙设计方法，采用的具体步骤如下：

步骤一：利用医学成像及图像处理技术提取出患者病变部位的三维外部轮廓图像；

步骤二：利用四面体网格生成技术生成三维外部轮廓图像外形的内部四面体网格；

步骤三：设定3D打印方向，并对三维外部轮廓图像进行切片处理；

步骤四：选定三维外部轮廓图像的一层切片平面；

步骤五：利用扫描线方法在选定的切片平面上得到特定方向的扫描线；

步骤六：找出所有从扫描线起点到终点间的四面体T；

步骤七：计算与扫描线相交的每个四面体T的相交起点A和相交终点B；

步骤八：利用四面体T的四个顶点，建立四元隐式函数，并利用该四元隐式函数计算出起点A到终点B所有采样点的函数值；

步骤九：对采样点函数值与设定值C进行比较，大于设定值C则跳过，小于则进入步骤十；