基于多轴修复软骨的轨迹规划方法技术

本申请提供一种基于多轴修复软骨的轨迹规划方法，其特征在于：所述方法包括如下步骤：S1：确定待修复软骨损伤部位的三维模型；S2：根据损伤部位的三维模型确定喷头轨迹；S3：通过喷头轨迹，确定多自由度关节式机器人的打印轨迹；S4：根据多自由度关节式机器人的打印轨迹，修复软骨的损伤部位。本申请的基于多轴修复软骨的轨迹规划方法可准确对比复原软骨损伤部分，并通过损伤部分的三维模型快速打印出损伤部分，完成软骨修复；本方法可有效规避自体软骨来源有效的问题，打印的损伤部分匹配度高，不会造成额外损伤，快速修复软骨损伤。

【技术实现步骤摘要】
基于多轴修复软骨的轨迹规划方法
本专利技术涉及生物3D打印
，尤其涉及一种基于多轴修复软骨的轨迹规划方法。
技术介绍
随着经济的发展、人们物质生活水平的提升，运动损伤、特别是关节软骨磨损的发病率较以前显著增加。全球骨关节病发病率约为3.54％，其中，中国超过4956万人。损伤后由于软骨组织缺乏血管、神经和淋巴管分布，经常无法自然愈合。目前临床主流的治疗软骨缺损的手段为自体软骨移植，其存在自体软骨来源有限、会造成额外损伤、手术较复杂且治疗费用昂贵等不足之处。关节软骨缺损由于缺乏自愈性目前仍没有有效的治疗方法。传统疗法如骨髓刺激技术和关节置换术在减轻疼痛方面是有效的，但它们不能再生成具有正常形态和功能的健康透明软骨。另外，由于这些治疗方法无法逆转软骨损伤，损伤部位可能会恶化并需要进行二次手术，严重者将导致关节功能障碍和永久性残疾。如今，软骨组织工程已成为修复软骨缺损的另一种有前景的方法。与传统3D打印技术相比，生物3D打印技术可以实现多种种子细胞在空间各异性的高密度负载及复杂活性微环境的构建。在骨与软骨缺损的治疗中，生物3D打印不仅可以克服自体骨与软骨来源不足的问题，而且还能解决传统人工骨软骨形状无法个性化匹配及生物活性不足等问题。增材制造技术已被广泛应用于快速、低成本地制造三维物体。然而，由具有曲率的复杂几何形状组成的构件对于传统的AM系统来说是一个具有挑战性的过程，由于打印机只有X,Y,Z这3个自由度，其能力仅限于平面分层打印。而骨骼形状相比一般零件更为复杂，表面由许多复杂的曲面组成，使用传统打印方法很难保证打印骨骼的表面的质量。仅使用平面水平层构造部件限制了基于材料沉积的增材制造工艺的能力。通过在复杂的曲面层上沉积材料，使用6自由度关节式机器人臂进行材料沉积可以显著扩展增材制造工艺的能力。获得正确的打印纤维取向对这些骨骼的正常工作是至关重要的，正确的打印纤维取向可以明显的提高结构强度。使用传统的平面层材料沉积工艺打印它们，在层的平面定向纤维，这通常会导致不良的纤维取向。打印骨骼这样的零件需要能够沿着三维曲线定位纤维，为了满足这一要求，需要具有沿曲面层沉积材料的能力。因此，本申请提供一种采用6自由度关节式机器人臂的打印损伤软骨的打印轨迹规划方法。
技术实现思路
有鉴于此，本专利技术提供一种基于多轴修复软骨的轨迹规划方法，其特征在于：所述方法包括如下步骤：S1：确定待修复软骨损伤部位的三维模型；S2：根据损伤部位的三维模型确定多自由度关节式机器人的喷头轨迹；S3：通过喷头轨迹，确定多自由度关节式机器人的打印轨迹；S4：根据多自由度关节式机器人的打印轨迹，修复软骨的损伤部位。进一步，利用医学CT扫描人体损伤软骨和正常软骨得到DICOM格式影像图，通过DICOM格式影像图生成软骨三维模型，通过损伤软骨和正常软骨三维模型的对比获得软骨损伤部位的三维模型，得到软骨损伤部位的三维模型的STL文件。进一步，所述步骤S2包括：S21：从STL文件中提取顶点、面和三角形面片单位法向量的信息，并分别存储在集合V、F和N中：V＝{v1,v2,v3,…,vm}，vm＝{vmx,vmy,vmz}，其中，V表示STL文件中所有顶点的集合，vm代表点xyz信息的集合；F＝{f1,f2,f3,…,fm}，fm＝{vm1,vm2,vm3}，其中，F表示STL文件中所有面的集合，fm代表第m个面包含的顶点集合；N＝{n1,n2,n3,…,nm}，其中，N表示STL文件中所有三角形面片单位法向量的集合；S22：集合V中顶点Vw的x、y坐标的范围在XOY坐标系生成投影，然后根据投影生成投影参考平面，其中，w＝1,2,3…m表示正整数；其中，参考平面中x、y范围如下所示：min(x)-d≤x≤max(x)+d，min(y)-d≤y≤max(y)+d，其中，d表示一个使参考平面的大小大于投影部分的大小的正变量；S23：在投影参考平面上，生成喷头移动的平面轨迹，并将轨迹点的数据记录到集合Pgridis中；Pgridis＝{p1,p2,p3,…,pm}，其中，Pgridis表示投影参考平面轨迹中所有点的有序集合；S24：从平面轨迹中识别STL模型的底面，其方法如下：S241：将STL模型的法线系数k为负的法线存储到集合N中的子集合Nbottom中；Nbottom＝{n1,n2,n3,…,nm}，其中，Nbottom表示STL模型文件中所有底面法线的集合，N表示STL文件中所有三角形面片单位法向量的集合；S242：将集合Nbottom中的元素的顶点投影到参考平面上；S243：将集合Nbottom中的元素投影到参考平面上的三角形内的轨迹点储存到集合Pin中；Pin＝{pj1,pj2,pj3,…,pjm，j∈[1,m]}，其中，Pin表示参考平面三角形内的轨迹点，j表示三角形的序号；S25：确定参考平面轨迹点在STL模型底面对应的三维空间轨迹点：S251：构建STL模型三角面片的顶点生成平面模型，所述模型为：Ax+By+Cz+D＝0，其中，A表示x的系数，B表示y的系数变量，C表示z的系数变量，D表示常数，x表示三角面片的顶点x坐标，y表示三角面片的顶点y坐标，z表示三角面片的顶点z坐标；S252：将集合Pin中第m个三角面片对应三角形中的轨迹点的x，y坐标值代入对应平面方程中得到对应的z轴坐标，得到的参考平面轨迹点在STL模型底面的对应三维空间轨迹点；S26：将得到三维空间轨迹点坐标和它对应STL模型三角面片的单位法向信息储存到集合Pin；projected中；Pin；projected＝{v1，v2，v3，...，vm；n1，n2，n3，...，nm}，其中，vm表示三维空间轨迹点，nm代表轨迹点的法向信息；S27：通过欧拉角确定喷头姿态信息，将得到喷头的姿态信息同位置信息一起储存在集合Pin；projected；Rmx＝-tan-1(nmy/nmz)Rmy＝-tan-1(nmx/nmz)其中，Rmx表示第m个轨迹点在x方向的旋转量，Rmy表示表示第m个轨迹点在y方向的旋转量，nmx表示第m个三维空间轨迹点的法向量在x轴的分量，nmy表示第m个三维空间轨迹点的法向量在y轴的分量，nmz表示第m个三维空间轨迹点的法向量在z轴的分量；S28：将得到的空间轨迹点生成连续的喷头移动轨迹；S29：根据喷头直径提高z轴高度，w＝w+1，其中，w表示正整数自然数，判断w是否大于集合V的元素个数，若是，则退出，若否，则进入步骤S2。进一步，步骤S23在投影参考平面上生成的喷头移动的平面轨迹绕预设的角度绕轴进行旋转，所述旋转用于曲面层与层之间的轨迹交叉，提高结构强度。进一步，S31：在机器人轨迹路径点上的每个圆锥表面对可能的喷头倾角进行取样，并将样本存储在集合Pexample中；S32：确定机器人的喷头表面本文档来自技高网...

【技术保护点】
1.一种基于多轴修复软骨的轨迹规划方法，其特征在于：所述方法包括如下步骤：/nS1：确定待修复软骨损伤部位的三维模型；/nS2：根据损伤部位的三维模型确定多自由度关节式机器人的喷头轨迹；/nS3：通过喷头轨迹，确定多自由度关节式机器人的打印轨迹；/nS4：根据多自由度关节式机器人的打印轨迹，修复软骨的损伤部位。/n

【技术特征摘要】
1.一种基于多轴修复软骨的轨迹规划方法，其特征在于：所述方法包括如下步骤：
S1：确定待修复软骨损伤部位的三维模型；
S2：根据损伤部位的三维模型确定多自由度关节式机器人的喷头轨迹；
S3：通过喷头轨迹，确定多自由度关节式机器人的打印轨迹；
S4：根据多自由度关节式机器人的打印轨迹，修复软骨的损伤部位。


2.根据权利要求1所述基于多轴修复软骨的轨迹规划方法，其特征在于：所述步骤S1包括：利用医学CT扫描人体损伤软骨和正常软骨得到DICOM格式影像图，通过DICOM格式影像图生成软骨三维模型，通过损伤软骨和正常软骨三维模型的对比获得软骨损伤部位的三维模型，得到软骨损伤部位的三维模型的STL文件。


3.根据权利要求2所述基于多轴修复软骨的轨迹规划方法，其特征在于：所述步骤S2包括：
S21：从STL文件中提取顶点、面和三角形面片单位法向量的信息，并分别存储在集合V、F和N中：
V＝{v1,v2,v3,…,vm}，vm＝{vmx,vmy,vmz}，其中，V表示STL文件中所有顶点的集合，vm代表点xyz信息的集合；
F＝{f1,f2,f3,…,fm}，fm＝{vm1,vm2,vm3}，其中，F表示STL文件中所有面的集合，fm代表第m个面包含的顶点集合；
N＝{n1,n2,n3,…,nm}，其中，N表示STL文件中所有三角形面片单位法向量的集合；
S22：集合V中顶点Vw的x、y坐标的范围在XOY坐标系生成投影，然后根据投影生成投影参考平面，其中，w＝1,2,3…m表示正整数；
其中，参考平面中x、y范围如下所示：
min(x)-d≤x≤max(x)+d，min(y)-d≤y≤max(y)+d，其中，d表示一个使参考平面的大小大于投影部分的大小的正变量；
S23：在投影参考平面上，生成喷头移动的平面轨迹，并将轨迹点的数据记录到集合Pgridis中；
Pgridis＝{p1,p2,p3,…,pm}，其中，Pgridis表示投影参考平面轨迹中所有点的有序集合；
S24：从平面轨迹中识别STL模型的底面，其方法如下：
S241：将STL模型的法线系数k为负的法线存储到集合N中的子集合Nbottom中；
Nbottom＝{n1,n2,n3,…,nm}，其中，Nbottom表示STL模型文件中所有底面法线的集合，N表示STL文件中所有三角形面片单位法向量的集合；
S242：将集合Nbottom中的元素的顶点投影到参考平面上；
S243：将集合Nbottom中的元素投影到参考平面上的三角形内的轨迹点储存到集合Pin中；
Pin＝{pj1,pj2,pj3,…,pjm，j∈[1,m]}，其中，Pin表示参考平面三角形内的轨迹点，j表示三角形的序号；
S25：确定参考平面轨迹点在STL模型底面对应的三维空间轨迹点：
S251：构建STL模型三角面片的顶点生成平面模型，所述模型为：Ax+By+Cz+D＝0，其中，A表示x的系数，B表示y的系数变量，C表示z的系数变量，D表示常数，x表示三角面片的顶点x坐标，y表示三角面片的顶点y坐标，z表示三角面片的顶点z坐标；
S252：将集合Pin中第m个三角面片对应三角形中的轨迹点的x,y坐标值代入对应平面方程中得到对应的z轴坐标，得到的参考平面轨迹点在STL模型底面的对应三维空间轨迹点；
S26：将得到三维空间轨迹点坐标和它对应STL模型三角面片的单位法向信息储存到集合Pin；projected中；
Pin；projected＝{v1,v2,v3,…...

【专利技术属性】
技术研发人员：唐倩，冯琪翔，郭伏雨，张志豪，曹粮玉，
申请(专利权)人：重庆大学，
类型：发明
国别省市：重庆;50