[发明专利]一种基于应力分布的Voronoi多孔梯度结构生成方法

说明书

本发明公开了一种基于应力分布的Voronoi多孔梯度结构生成方法，包括以下步骤：步骤1：在三维空间内，对于给定的任意二维模型，获取经有限元分析后带有模型边界及应力等值线的样条曲线；步骤2：在给定的模型边界内，随机生成一定数量的离散点；步骤3：基于两临近应力等值线与边界所围区域，控制所述种子点的分布；步骤4：基于所述步骤3获得的种子点进行基于质心的Voronoi划分；步骤5：对所述步骤4中划分得到的Voronoi多边形的每条边界赋予一定宽度，用以表征多孔梯度结构的实体部分，其余则为多孔梯度结构的孔隙部分，向三维空间映射得到Voronoi多孔梯度结构。本发明根据应力分布来对多孔结构的孔隙分布进行精确的设计，满足使用过程中的力学期望。

技术领域

本发明涉及一种面向3D打印的多孔结构设计方法领域，尤其是涉及一种基于应力分布的Voronoi多孔梯度结构生成方法。

背景技术

3D打印与传统加工方式最大的不同是该技术将模型分解为一系列的层面数据，采用自上往下，层层叠加的加工方式，不仅能够控制宏观外形，还可以有效的控制内部的微观结构，凭借其独特的加工方式，对制备复杂零部件具有极大的优势，可以打印结构非常复杂的模型。这也使得制备多孔异质性结构的零件成为可能。多孔结构具有相对密度小、比强度高、单位体积表面积大等特点，在吸能减震、轻量化等领域方面显示出了优越性，已广泛应用于高科技领域和普通工业领域，例如航天器保护壳、汽车缓冲件等。另外，在医疗领域，具有多孔结构的植入物有利于骨组织在其中的生长，可增强植入物与骨骼之间的连接，达到与人体骨骼的弹性模量相匹配的数值。

调研国内外文献，目前对二维多孔结构的研究主要集中在规则蜂窝和随机Voronoi分布两种结构，但这些基于规则或者随机排列的多孔结构往往存在极易忽视的问题，即在实际的应用中，不仅要求多孔结构具有高孔隙率和比表面积，而且要求多孔结构满足使用过程中的力学期望。然而，如果能实现根据应力分布来对孔隙分布进行精确的设计，该问题就可以迎刃而解。

发明内容

本发明为了克服现有技术的不足，提供一种的基于应力分布的Voronoi多孔梯度结构生成方法。

为了实现上述目的，本发明采用以下技术方案：一种基于应力分布的Voronoi多孔梯度结构生成方法，包括以下步骤：

步骤1：在三维空间内，对于给定的任意二维模型，获取经有限元分析后带有模型边界及应力等值线的样条曲线；

步骤2：在给定的模型边界内，随机生成一定数量的离散点；

步骤3：基于两临近应力等值线与边界所围区域，控制所述种子点的分布；

步骤4：基于所述步骤3获得的种子点进行基于质心的Voronoi划分；

步骤5：对所述步骤4中划分得到的Voronoi多边形的每条边界赋予一定宽度，用以表征多孔梯度结构的实体部分，其余则为多孔梯度结构的孔隙部分，向三维空间映射得到Voronoi多孔梯度结构。

可选的，所述步骤1中，所述有限元分析所施加的载荷边界条件与模型在实际使用中相同。

可选的，所述步骤3中，所述两临近应力等值线与边界所围区域视为应力处处相等区域。

可选的，所述步骤3中，所述种子点数量由两临近应力等值线与边界所围区域面积和区域内的应力值确定，以种子点密度，即单位面积内种子点的数量作为度量，应力值大的区域种子点密度大。

可选的，所述步骤4中，采用基于质心的Voronoi划分方法，使得每个多边形仅包含一个基点，且多边形内的点到基点的距离小于到其它点的距离，对于步骤3中的点集(p₁,p₂……p_n),定义的控制区域为R_i＝{x∈X|d(x,P_i)d(x,P_j),j＝{1,2,3L,nL},j≠i}，其中d为两点的欧氏距离。