[发明专利]一种面向3D打印的物体内部结构优化方法

说明书

技术领域

本发明涉及本发明涉及一种面向3D打印的物体内部结构优化方法，具体涉及一种基于Voronoi图和有限元分析的优化物体强度体积比的方法。

背景技术

3D打印，即快速成型技术的一种，它是一种以数字模型文件为基础，运用粉末状金属或塑料等可粘合材料，通过逐层打印的方式来构造物体的技术。

最近几年，3D打印越来越多的被人们所提及，将数字模型文件打印成实物的能力，使得它在许多行业得到广泛应用。从技术实现上来看，3D打印通常是将打印材料“一层一层”地叠加起来，从而将数字蓝图变为实物。所以，模型的体积直接关系到打印过程中的耗材，大型复杂的模型将会耗费大量的材料。而打印过程的耗材直接决定了打印的成本。针对这个问题，Wang等人在2013年提出了用蒙皮框架结构代替内部材料的方法(WANG,W.,WANG,T.Y.,YANG,Z.,LIU,L.,TONG,X.,TONG,W.,DENG,J.,CHEN,F.,AND LIU,X.2013.Cost-effective printing of3D objects with skin-frame structures.ACM Transactions on Graphics(Proc.SIGGRAPH Asia)32,5)，为了减少内部材料，该方法将问题抽象为在物理和几何的条件约束下的框架结构(数目)的优化问题，但是该方法只考虑了自承重的情形，即只考虑所打印物体受重力的影响，并未考虑到物体受外力的情形，所以该方法的结果虽然可以很大程度的节省材料，但是所打印出来的物体不能保证承受足够的外力。另一方面，出于对打印物体的平衡性的考虑，Prevost等人在2013年提出的方法(PR′EVOST,R.,WHITING,E.,LEFEBVRE,S.,AND SORKINE HORNUNG,O.2013.Make it stand:balancing shapes for3D fabrication.ACM Trans.Graph.32,4(July),81:1–81:10.)，该方法可以使物体保持平衡，同时也能在一定程度上节省材料，该方法将这一问题形式化为包含平衡性与形变程度两项加权和的目标函数，并将问题定义为能量最小优化问题，来求解出使目标函数目标最优的模型，但是该方法仅仅是为了保持物体站立平衡。

上述的方法在减少打印耗材的同时，都在一定程度上改变了打印物体的物理强度，因为物体的用料与物理强度是紧密相连的。考虑到这一点，Zhou等人在2013年将这样一个结构问题看做带约束的优化问题(ZHOU,Q.,PANETTA,J.,AND ZORIN,D.2013.Worst-case structural analysis.ACM Trans.Graph.32,4(July),137:1–137:12)，该文章为我们提供了物体结构分析的方法，该方法可以对于任意的模型，计算分析得出模型中最薄弱易碎的区域。另外，Stava等人在2012年提出了迭代优化的方法(STAVA O.VANEK,J.BENES,B.CARR,N.,AND MˇE545CH,R.2012.Stress relief:improving structural strength of3D printable objects.ACM Trans.Graph.31,4(July),48:1–48:11),通过在模型内部挖洞，对模型表面增厚和增加内部支撑结构的方式来增强物体的结构强度，同时也达到了节省材料的目的，问题在于该方法会改变物体的表面结构，对模型的外观产生影响。

发明内容

本发明为了解决上述问题，提出了一种面向3D打印的物体内部结构优化方法，该方法基于Voronoi图和FEM(Finite Element Analysis Method,有限元分析)，通过FEM计算模型在外力下的应力图，并通过Voronoi图结合应力图计算并产生类蜂窝状的内部结构；蜂窝状的结构可以在提供很好的结构强度的同时最小化耗材，不仅可以在3D打印的过程中节省材料，还可以赋予打印物体类蜂窝状的内部结构，承受给定的外力，使得物体更加坚固、结实。

为了实现上述目的，本发明采用如下技术方案：

一种面向3D打印的物体内部结构优化方法，包括以下步骤：

(1)计算初始给定的实心模型S的应力分布SM(S，F)，确定目标优化函数的参数值，并将其初始化；