基于参数曲面的增材制造保形晶格结构生成方法

2021年6月图 学 学 报 June2021第42卷第3期JOURNAL OF GRAPHICS V ol.42No.3基于参数曲面的增材制造保形晶格结构生成方法
肖文磊，林载盛，熊唱日，王世平，魏巍，赵罡
(北京航空航天大学机械工程及自动化学院，北京 100191)
摘要：晶格结构因其具备特殊的机械性能，已成为增材制造复杂结构设计和制造的重要研究领域之一。

针对传统的模型晶格结构生成方法需要通过对参数化建模的晶格结构网格进行裁剪或对其网格进行保形变形实现，且生成效率较低的问题，提出一种基于参数曲面的增材制造保形晶格结构生成方法，实现晶格结构对曲面空间的适应和高效生成。

首先，基于一种矩阵方法完成晶格结构骨架的表达和构造。

其次，利用点与曲面、曲线与曲面、曲面与曲面所形成的3种封闭空间，对晶格结构骨架进行保形变形，使晶格结构适应曲面空间。

最后，采用基于晶格结构骨架的网格生成和拼接方法，生成保形晶格结构的网格模型。

通过组件应用架构(CAA)二次开发方法提取计算机辅助三维交互应用软件(CATIA)模型的参数曲面，实现保形晶格结构模型的高效生成，并表现出良好的曲面空间适应性，证明该方法具有一定工程价值。

关键词：增材制造；晶格结构；参数曲面；保形变形；组件应用架构
中图分类号：TP 391 DOI：10.11996/JG.j.2095-302X.2021030517
文献标识码：A 文章编号：2095-302X(2021)03-0517-08
Parametric surface-based additive manufacturing conformal
lattice structure generation method
XIAO Wen-lei, LIN Zai-sheng, XIONG Chang-ri, WANG Shi-ping, WEI Wei, ZHAO Gang (School of Mechanical Engineering and Automation, Beihang University, Beijing 100191, China)
Abstract: The lattice structure has become one of the most important research fields in the design and manufacture of complex structures in additive manufacturing, because of its special mechanical properties. The conventional method for model lattice structure generation is achieved by trimming the parametric modeling lattice structure grid or by performing conformal deformation of the lattice structure. Such methods are relatively inefficient. A conformal lattice structure generation method was proposed based on parametric surface for additive manufacturing, which can realize the adaptation and efficient generation of lattice structure in surface space. Firstly, a matrix-based approach was proposed to express and construct the lattice structure frame. Then, three types of enclosed spaces, which were formed by point-surface, curve-surface, and surface-surface patterns, were implemented to perform conformal transformation on the lattice frame to adapt the surface space. Finally, a grid generation and splicing method based on the lattice structure skeleton was employed to generate the mesh model of the conformal lattice structure. The parametric surface was extracted from the CATIA model based on CAA to generate the conformal lattice structures effectively, and the structures can adapt to the surface space well, which proves the potential value for the engineering implementation and application of the method.
Keywords: additive manufacturing; lattice structure; parametric surface; conformal transformation; component application architecture
收稿日期：2020-09-22；定稿日期：2020-12-15
Received：22 September，2020；Finalized：15 December，2020
第一作者：肖文磊(1984-)，男，江西井冈山人，副教授，博士。

主要研究方向为智能数控与机器人、航空智能制造、生物图像分析等。

E-mail：*******************.cn
First author：XIAO Wen-lei (1984-), male, associate professor, Ph. D. His main research interests cover intelligent CNC and robots, intelligent manufacturing inaviation,bioimageanalysis,etc.E-mail:*******************.cn
518 数字化设计与制造2021年
增材制造是一种利用材料逐层堆积形成零件的制造技术，比起传统的制造方法，其所生成的零件可以拥有更复杂的几何结构、更优秀的力学性能以及更高的材料利用率，在复杂零件制造领域应用前景巨大。

由于增材制造的制造方式是通过逐层打印的方式形成零件，所以该技术对零件的几何结构几乎没有限制，只要存在零件CAD模型就能打印制造，这为具有镂空结构、复杂曲面等特殊结构的特种零件制造提供了基础。

增材制造对比传统制造的一大优势在于零件的轻量化，通过晶格结构、镂空结构、拓扑优化等设计技术，使零件材料分布在需要的区域，提高零件的强度重量比，在航空航天和汽车工业有重要应用。

但是，对比快速发展的增材制造工艺和材料技术，轻量化零件设计的研究进展相对较少[1]。

目前，增材制造零件轻量化设计的主要方法有晶格结构设计和拓扑优化2种，同时也是增材制造的前沿热点研究方向。

WEEGER等[2]提出一种数字设计结合非线性仿真的软体晶格结构设计方法，实现对弯曲杆晶格结构件的设计和仿真。

DAYNES等[3]提出一种晶格结构优化设计方法，利用拓扑优化获得的主应变场调整晶格桁架的尺寸形状，大大提高晶格结构的强度和刚度。

GOEL和ANAND[4]基于B 样条曲面的功能渐变晶格结构设计方法，实现渐变晶格单元间的光滑对接，并结合SIMP拓扑优化方法构建出高刚度的功能渐变晶格结构。

DUMAS等[5]采用不同孔隙率渐变晶格结构进行骨骼置换的新颖方法，利用MATLAB程序生成晶格结构模型STL文件并通过有限元分析和拉伸试验预测和验证了该结构应用于骨骼置换的可行性。

除此之外，越来越多不同材料、尺寸和结构的晶格结构被应用于声学、介电、机械、生物医学和航空航天方面的研究[6-11]。

基于晶格结构的增材制造零件轻量化的核心思想是利用晶格结构填充零件内部，提高零件强度重量比、节省材料以及赋予零件特殊的机械性能。

传统的晶格结构由晶格单元在正交笛卡尔坐标系的3个方向堆叠形成。

计算机构造晶格结构模型有2种方法[12]：①采用参数化建模的方法，先构建晶格单元的实体模型，并基于实体模型生成晶格单元的网格模型，然后通过布尔运算求并集的方法实现参数晶格单元的堆叠，形成晶格结构的实体模型，最后对其进行表面网格剖分，用于后续的打印制造。

该生成方法是基于参数化建模，优点是可以构造比较复杂的晶格单元，方便进行局部构型修改。

缺点是当晶格单元数目庞大时，晶格结构实体的生成速度十分缓慢，很难对参数化模型进行变形处理；②通过参数驱动直接构造晶格结构的网格模型，不需要进行参数化建模以及布尔运算，构造速度快。

直接构建晶格结构网格模型的方法需要构造晶格的拓扑结构。

其是晶格结构的骨架模型，所需要的拓扑元素包括节点和边，由节点和边形成晶格单元骨架模型，再通过其拼接形成整体晶格结构的骨架，可称为晶格框架模型(frame)，框架模型形成过程如图1所示。

在晶格单元骨架的拼接过程中会出现节点和边元素的重合问题，而后续生成网格实体需要根据这些拓扑元素生成对应网格进行拼接，重复的节点和边会生成重复的实体网格，影响后续打印制造。

为了防止生成重叠的网格，需要去除重复的拓扑元素，即保证在同一空间位置内存在唯一的拓扑元素。

图1 晶格拓扑结构生成过程
Fig. 1 Construction process of lattice topology 构造模型晶格结构的方法有裁剪方法和保形方法。

裁剪方法是在已经预先设计晶格结构的基础上，通过模型边界对预设晶格结构进行裁剪，构建模型晶格结构，如图2(b)所示。

采用裁剪方法构造的模型晶格结构相对于预设晶格结构，由于裁剪的原因会导致部分晶格单元不完整，丢失拓扑结构的情况。

除此之外，裁剪生成的零件内部晶格结构不能很好地适应零件表面的曲面特征，可能对零件机械性能产生影响。

保形方法是基于模型表面形状对晶格结构进行保形变形，使其适应零件区域，如图2(c)所示。

传统的保形晶格结构生成方法[13]需要基于参数化建模的晶格结构生成其表面网格，然后对网格进行保形变形，生成保形晶格结构。

由于涉及网格剖分和网格布尔运算，传统的保形晶格结构生成效率比较低。

图2 晶格结构二维示意图((a)预设晶格结构；(b)裁剪
晶格结构；(c)保形晶格结构)
Fig. 2 Schematic diagram of lattice structure in 2D ((a) Preset lattice structure; (b) Trimed lattice structure;
(c) Comformal lattice structure)
(a)(b)(c)
第3期 肖文磊，等：基于参数曲面的增材制造保形晶格结构生成方法 519
本文提出一种基于参数曲面的增材制造保形晶格结构生成方法，即基于矩阵和参数曲面构造晶格拓扑结构，并快速生成网格模型。

该方法生成的晶格结构无需裁剪即可适应曲面空间，对比传统的保形晶格生成方法也更加高效。

本文最后通过CATIA 二次开发的方法，采用该参数曲面生成模型的保形晶格结构，模型表现出良好的空间适应性。

1 基于参数曲面的保形晶格结构生
成方法
基于参数曲面的保形晶格结构生成方法如图3所示。

首先基于一种矩阵方法构建晶格拓扑结构，实现晶格结构骨架的高效生成。

然后在曲面与点、曲面与曲线、曲面与曲面形成的3种曲面封闭空间，基于参数曲面对晶格拓扑结构进行保形变形。

最后基于网格拼接的方法，利用保形变形后的拓扑骨架模型生成晶格结构的网格模型。

图3 基于参数曲面的保形晶格结构生成方法 Fig. 3 Parametric surface-based conformal lattice
structure generation method
1.1 晶格拓扑结构生成
定义晶格单元节点拓扑位置Node ，以及节点
在笛卡尔坐标系3个正交方向向量,,U V W
上的最大节点数U ，V ，W 为
(,,),,Node u v w u U v V w W =⎧⎨
<<<⎩
(1) 其中，u ，v ，w 为自然数，表示晶格单元的节点拓
扑位置坐标值。

假设某晶格单元拥有n 个节点，定义晶格单元
节点的拓扑坐标矩阵为
Node 3U V W n A A A ⨯⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦ M (2) 其中，向量U A ，V A ，W A 分别为,,U V W
方向上各节点的拓扑坐标对应的方向向量，即
121212[,,,][,,,][,,,]
U U U U V V V V W W W W n
n
n
ααααααααα=⎧⎪=⎨⎪=⎩
A A A (3) 其中，Ui α ，Vi α ，Wi α 分别为节点在U ，V ，W
方
向上对应的拓扑位置坐标，且Ui
u α∈ ，Vi v α∈ ，Wi w α∈ 。

节点的拓扑位置可用拓扑坐标矩阵M Node
的列向量表出，即
(,,)i Ui
Vi Wi
Node ααα= (4) 图4为晶格单元的节点拓扑坐标矩阵案例，U
，V ，W
对应的最大节点数量均为3，所有节点的 拓扑位置坐标均由拓扑坐标矩阵M Node 的列向量 表出。

图4 晶格单元节点坐标矩阵
Fig. 4 Coordinate matrix of lattice cell nodes
最后定义拓扑坐标与几何坐标的映射矩阵为
121212T 3max o (,,po-Ge )
o
,,,,,,,,,U V W u u u v v v w w w U V W x x x y y y z z z ⨯⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦
M (5)
其中，i u x ，i v y ，i w z 分别对应节点拓扑坐标u i ，v i ，
w i 的几何坐标，规定M Topo-Geo 的第一列坐标值均为0。

由M Node 和M Topo-Geo 可以完整而高效表示晶格单元节点的拓扑和几何位置。

完成晶格单元节点定义后，通过将晶格单元的
拓扑构造成无向图C =<Node , Edge >，其中Node = {Node i |i =1,2,···,n }表示节点，Edge ={Edgek |k =1,2,
···,m }表示边，Edgek =(Node i , Node j )，即可完成晶格
单元结构的表示。

通过矩阵运算快速拼接各晶格单
元，形成晶格结构骨架模型，以完成晶格拓扑结构
生成和重复拓扑元素的剔除。

1.2 参数曲面保形变形 生成晶格结构骨架模型后，可根据参数曲面形
状进行保形变形，使适应点与曲面、曲线与曲面、
曲面与曲面所形成3种曲面空间。

首先定义3种封闭曲面空间为
.to..to..to.{(,,),(,,,,,)} {(,,,,,),(,,,,,)} {(,,,,,),(,,,,,)}p s c s s s P x y z S x y z u v w C x y z u v w S x y z u v w S x y z u v w S x y z u v w ϕϕϕ=⎧⎪=⎨⎪=⎩ (6) 其中，.to.p s ϕ，.to.c s ϕ，.to.s s ϕ分别为3种封闭曲面
520 数字化设计与制造 2021年
空间；P ，C ，S 分别为三维空间中的参数点、参数曲线以及参数曲面；x ，y ，z 为几何空间坐标系坐标，表示参数点、参数曲线和参数曲面在几何空间的位置；u ，v ，w 为拓扑空间坐标系坐标，用于建立保形变形前晶格结构骨架与保形变形后骨架的映射关系，规定0,,,1u v w ≤≤。

1.2.1 点与曲面空间保形变形
设封闭空间.to.p s ϕ由参数点p P 和参数曲面S 围成，晶格结构骨架与封闭空间.to.p s ϕ的映射关系如图5所示。

图5 点与曲面封闭空间映射方法
Fig. 5 Mapping method of point-surface closed space
设点P 0为晶格结构骨架上的拓扑点，拓扑空间坐标为0000,),(P u v w 。

参数点p P 在拓扑空间坐标为0,(1)0,p P ，几何空间坐标为,,()p p p p P x y z 。

参数曲
面上的点P s 在拓扑空间坐标为,,(0)s s s P u v ，对应的几何空间坐标为,,()s s s s P x y z ，可表示为
,,,,,0()s s s s s s P x y z u v 。

则保形变化后P 0在几何空间坐标为0,(,)P x y z 的计算方法为
000 () () ()s p s s p s s p s x x w x x y y w y y z z w z z ⎧=+⋅-⎪
=+⋅-⎨⎪=+⋅-⎩
(7)
其中，,,s s s x y z 的坐标值由00,,,,,0()s s s s P x y z u v 确
定，,,p p p x y z 为参数点几何空间坐标,,()p p p p P x y z 。

1.2.2 曲线与曲面空间保形变形
设封闭空间.to.c s ϕ由参数曲线C 和参数曲面S
围成，晶格结构骨架与封闭空间.to.c s ϕ的映射关系如图6所示。

图6 曲线与曲面封闭空间映射方法
Fig. 6 Mapping method of curve-surface closed space
设点P 0为晶格结构骨架上的拓扑点，拓扑空间
坐标为0000,),(P u v w 。

参数曲线上的点P c 在拓扑空间坐标为,1()0,c c P u ，对应的几何空间坐标为
,,()c c c c P x y z ，可表示为,,,,0,1()c c c c c P x y z u 。

参数曲面上的点P s 在拓扑空间坐标为,,(0)s s s P u v ，对应的几何空间坐标为,,()s s s s P x y z ，可表示为
,,,,,0()s s s s s s P x y z u v 。

则保形变化后P 0在几何空间的坐标0,(,)P x y z 的计算方法为
000 () () ()s c s s c s s c s x x w x x y y w y y z z w z z =+⋅-⎧⎪=+⋅-⎨⎪=+⋅-⎩
(8) 其中，,,s s s x y z 的坐标值由00,,,,,0()s s s s P x y z u v 确定，,,c c c x y z 的坐标值由0,,,,0,1()c c c c P x y z u 确定。

1.2.3 曲面与曲面空间保形变形
设封闭空间.to.s s ϕ由参数曲面S 1和参数曲面
S 2围成，晶格结构骨架与封闭空间.to.s s ϕ的映射关系如图7所示。

图7 曲线与曲面封闭空间映射方法
Fig. 7 Mapping method of surface-surface closed space
设点P 0为晶格结构骨架上的拓扑点，拓扑空间坐标为0000,),(P u v w 。

参数曲面S 1上的点1s P 在拓扑
空间坐标为111,,1()s s s P u v ，对应的几何空间坐标为1111,,()s s s s P x y z ，可表示为11111,,,,,1()c s s s s s P x y z u v 。

参数曲面S 2上的点2s P 在拓扑空间坐标为
222,,0()s s s P u v ，对应的几何空间坐标为
2222,,()s s s s P x y z ，可表示为222222,,,,,0()s s s s s s P x y z u v 。

保形变化后P 0在几何空间的坐标0,(,)P x y z 的
计算方法为
2012201220
12 ()
() ()s s s s s s s s s x x w x x y y w y y z z w z z =+⋅-⎧⎪
=+⋅-⎨⎪=+⋅-⎩ (9)
其中，111,,s s s x y z 的坐标值由111100,,,,,()1s s s s P x y z u v 确
定，
,,c c c x y z 的坐标值由22200,,,,),(0c s s s P x y z u v 确定。

1.2.4 晶格结构骨架保形变形
根据上述3种方法对晶格结构骨架进行保形变形处理。

未处理的骨架由拓扑节点和正交的拓扑边组成，通过骨架节点的保形变形可得到节点保形
第3期肖文磊，等：基于参数曲面的增材制造保形晶格结构生成方法521
的骨架，如图8(a)所示。

通过对骨架节点和边同时
进行保形变形可得到边保形的晶格结构骨架，边的
保形变形方法是通过在每条拓扑边上插入等数量、
等距的插值点，对原有节点和新插入点进行保形变
形，得到边保形的骨架，如图8(b)所示。

图8 保形晶格结构骨架((a)节点保形骨架；(b)边保形骨架)
Fig. 8 Conformal lattice structure framework ((a) Node
conformal framework; (b) Edge conformal framework)
1.3 晶格网格实体生成
在晶格结构骨架模型的基础上，生成拓扑节点
和拓扑边的包围网格，通过网格拼接的方法形成晶
格结构实体模型。

增材制造常用的网格模型为三角
网格模型，所以需要基于三角网格进行网格的生成
和拼接。

本文采用一种基于晶格结构骨架拓扑点、
拓扑边的网格生成和拼接方法，生成骨架节点的凸
包网格结构和骨架边的多边形网格结构，然后通过
凸包网格和多边形网格的拼接生成晶格结构的三
角网格模型。

边包围网格的构造基于晶格结构骨架的拓扑
边单位，根据拓扑边节点(包括原有节点和保形变形
插入的节点)构造分段三角包围网格。

如果晶格结构
只对骨架节点进行保形变形，则根据图9(a)的方法
构造边包围网格，截面可为正N边形(3N
≤)，d
为与节点凸包网格拼接的预留距离。

如果同时对骨
架节点和边进行保形变形，则根据图9(b)的方法构
造边包围网格，边节点处截面处r与插入节点处截
面r'的关系为/sin
r rθ
'=。

图9 边包围网格构造方法((a)无插入节点边；
(b)有插入节点边)
Fig. 9 Edge Bounding mesh construction method ((a) Edge
without insertion point; (b) Edge with insertion points)
节点的包围网格构造基于三维凸包的方法。

凸
包是计算几何最基本的结构之一，三维凸包被广泛
运用于计算机图形建模领域[14]。

目前三维凸包生成
算法已经比较成熟，应用最广泛的算法是Clarkson-
Shor[15]以及QuickHull算法[16]。

本文采用QuickHull
算法，根据靠近骨架节点(非插入节点)的边包围网格
顶点N i，如图10所示，构造N i所形成的点集的三
维凸包网格，其中需要满足预留距离
min
(
/n/2)
ta
d rθ
≥，
min
θ为与节点相连边之间的夹
角的最小值，否则在节点处边包围网格会产生干涉。

图10 节点凸包网格构造方法
Fig. 10Node Convex hull mesh construction method
2 基于CAA的曲面提取与保形晶格
结构生成
2.1 基于CAA的参数曲面提取
CATIA CAA的几何建模器(CATIA geometry
modeler，CGM)提供了一系列高效完善的曲面建模
与操作接口与功能，可以便利地用于二次开发，满
足参数曲面设计与分析需求。

在CGM中，曲面由
参数U，V各自的标量函数FX，FY和FZ定义，是
一种在R2至R3空间的函数(图11)。

其标量函数代
表了曲面上每一点与笛卡尔坐标系间的映射关系，
曲面对象必须为C2连续。

图11 CGM中曲面参数定义
Fig. 11 Surface parameter definition of CGM
基于CGM的参数曲面定义，CATIA的曲面提
取过程分为以下步骤：
(1) 在CATIA模型空间中选取或创建一张扫
掠曲面；
(2) 运用CGM三角细分算法将参数曲面进行
离散处理，如图12所示。

其中Sag为细分线段距
离曲面的最大距离，Step为曲面细分步长，Angle
为细分角度。

通常情况下，Step和Angle都取默认
(a)(b)
522 数字化设计与制造 2021年
值(Step 设为无限长，Angle 设为90°)，并通过调整sag 值得到离散度适合于生成保形晶格结构的三角细分曲面网格。

并用迭代器(iterator)得到该曲面的离散数据点(即三角形各顶点)；
图12 CGM 曲面细分方法 Fig. 12 Tessellation method of CGM
(3) 利用CATSurface::GetParam 接口提取曲面离散点的u ，v 参数值，作为后续生成晶格结构骨架的拓扑参数结点。

2.2 基于参数曲面的保形晶格结构生成
基于CATIA 的参数曲面提取方法，本文的模型保形晶格结构生成流程如图13所示，分为以下步骤：
图13 模型保形晶格结构生成流程
Fig. 13 Generation process of model’s conformal lattice
structure
(1) 在CATIA 模型空间内创建或选取模型的参数曲面。

(2) 基于CATIA 参数曲面提取方法，提取参数曲面S 在拓扑空间,,(0)s s s P u v 对应的几何空间坐标,,()s s s s P x y z 。

通过提取的参数曲面和参数点、参数曲线构造上文的3种曲面空间。

(3) 设置初始晶格单元类型，如立方体晶格、
体心立方晶格和面心立方晶格等。

并设置U ，V
，W
3个方向的晶格单元数量。

(4) 根据矩阵方法构建晶格拓扑结构，然后基
于参数曲面对晶格拓扑结构进行保形变形，生成模型的保形晶格骨架结构。

(5) 基于上文1.3节方法，利用拓扑点和拓扑边生成晶格结构的网格模型。

接下来分别对3种曲面空间的保形晶格结构生成进行示例说明。

采用四边形立方体作为晶格单元，U ，V ，W 3个方向的晶格单位数量分别设置为10，10，6。

利用CGM 建立扫掠曲面生成保形晶格。

2.2.1 点与曲面晶格结构生成
首先通过空间点和扫掠曲面确定保形空间的U ，V ，W
3个方向，利用式(5)计算出保形变形后的晶格骨架结构。

其中晶格结构骨架拓扑点P 0在曲
面上的映射点P S 的空间坐标求解需要利用CAA 接口函数，输入P 0对应的u ，v 坐标进行求解。

生成的保形结构骨架如图14所示。

然后基于上文晶格网格实体的生成方法，生成半径r 为0.5 mm ，截面为16边形的杆和对应的凸包网格结构。

生成的保形晶格结构网格模型如图15所示。

2.2.2 曲线与曲面晶格结构生成
通过空间曲线和扫掠曲面确定保形空间的U
，V ，W
3个方向，并利用式(6)计算出保形变形后的
晶格骨架结构。

曲面点的求解方法同上。

生成的保形结构骨架如图16所示。

设置杆半径r 为0.5 mm ，杆截面为16边形。

生成的保形晶格结构网格模型如图17所示。

2.2.3 曲面与曲面晶格结构生成
通过2个扫掠曲面确定保形空间的U ，V ，W
3个方向，并利用式(7)计算出保形变形后的晶格骨架结构。

曲面点的求解方法同上。

生成的保形结构
骨架如图18所示。

设置杆半径r 为0.5 mm ，杆截面为16边形。

生成的保形晶格结构实体模型如图19所示。

2.3 实例验证
本文保形晶格结构生成算法已在Visual studio 2012环境中采用C++语言实现。

为了验证本文方法的高效性，在配置为英特尔i5-6500 CPU ，8 G 内存，Windows 7操作系统的PC 机上，将本文算法与专业晶格结构生成软件插件intralattice 进行保形晶格结构生成效率的对比。

以圆柱保形晶格结构为例，图20是本文与intralattice 生成的保形晶格结构网格模型对比图，圆柱体半径为50 mm ，高为100 mm 。

选取晶格单元为立方体晶格，令端面圆周方向为U ，圆周径向为V ，圆柱高度方向为W
，杆半径
第3期肖文磊，等：基于参数曲面的增材制造保形晶格结构生成方法523
为1 mm。

保形晶格结构网格模型的生成效率对比
见表1。

由表1
可见，本文的保形晶格结构网格模
型速度对比intralattice更快，并且随着晶格单元数
量的增加，速度优势会更明显。

虽然无法计算
intralattice的网格模型显示所需的额外时间，但是
仍可说明本文方法对比intralattice更加高效。

图14 点与曲面保形晶格骨架图15 点与曲面保形晶格结构网格模型图16 曲线与曲面保形晶格骨架Fig. 14 Point-surface conformal Fig. 15 Point-surface conformal Fig. 16 Curve-surface conformal lattice frame lattice structure mesh model lattice frame
图17 曲线与曲面保形晶格图18 曲面与曲面保形图19曲面与曲面保形晶格结构网格模型晶格骨架结构网格模型
Fig. 17 Curve-surface conformal Fig. 18 Surface–surface Fig. 19 Surface -surface conformal lattice structure mesh model conformal lattice frame lattice structure mesh model
图20 本文方法与intralattice的比较((a)本文保形晶格结
构；(b)intralattice生成的保形晶格结构)
Fig. 20 Comparision of the method in this article and
intralattice ((a) Aarticle’s conformal lattice structure; (b)
intralattice’s conformal lattice structure)
表1 5种常见晶格拓扑的特性对比
Table 1 Comparison of properties of five common cell
topologies
UVW晶格单元
数量(个)
本文方法
生成时间(s)
Intralattice
生成时间(s)
5×10×10 0.6 9.1
10×10×10 1.1 18.6
15×10×10 1.5 31.1
20×10×10 1.8 45.2
图21~25展示了基于本文方法生成的模型保形
晶格结构3D打印件样品，模型的保形晶格结构变
现出良好的曲面空间适应性。

图21 保形晶格结构样例3D打印成品
Fig. 21 3D printing samples of Conformal lattice structure
图22 涡轮叶片晶格模型与成品(立方体晶格单元，UVW晶
格单元数量为6×4×2)
Fig. 22 Turbine blade lattice model and sample (cube cell,
6×4×2 in UVW)
(a)(b)
524 数字化设计与制造2021年
图23机翼晶格模型与成品(立方体晶格单元，UVW晶格数
量为10×3×2)
Fig. 23 Aircraft wing lattice model and sample (cube cell,
10×3×2 in UVW)
图24 曲连杆晶格模型与成品(立方体晶格单元，UVW晶格
数量10×3×2)
Fig. 24 Aircraft wing lattice model and sample (cube cell,
10×3×2 in UVW)
图25 点与曲面保形晶格案例模型与成品(立方体晶格单
元，UVW晶格数量为10×3×3)
Fig. 25 Examples of Point-surface lattice model and sample
(cube cell, 10×3×3 in UVW)
3 结束语
本文针对传统保形晶格结构生成方法存在的问题，提出了一种基于参数曲面的增材制造保形晶格结构生成方法，能够高效生成保形晶格结构并对曲面空间表现出的良好适应性，为增材制造复杂零件晶格结构设计提供了一种新颖而高效的方法。

首先，采用基于矩阵的方法构造晶格拓扑结构骨架模型，实现晶格拓扑结构的快速生成；然后，基于点与曲面、曲线与曲面、曲面与曲面3种参数曲面封闭曲面空间，建立曲面空间与晶格拓扑结构的映射关系，实现晶格结构骨架的保形变形，使其适应3种曲面空间；最后，基于网格拼接的方法，利用晶格结构骨架拓扑元素生成包围网格，拼接形成晶格结构网格模型供后续打印和制造。

通过CAA 二次开发的方法提取Catia的参数曲面对该方法进行设计分析和验证，实例证明该方法能够高效生成保形晶格结构，并且所生成的晶格结构能良好地适应曲面空间，具有一定工程意义。

本文的研究集中在保形晶格结构的图形学模型生成方法，但部分问题仍有待进一步研究。

首先，保形晶格结构与裁剪晶格结构的机械性能存在差异性，对两者进行系统的性能分析和对比是往后需要研究的内容。

其次，本文提出的保形晶格生成方法适用于能够用曲面表示的封闭空间，无法扩展到任意空间结构，更通用的保形晶格结构生成方法还有待进一步研究。

参考文献(References)
[1]AREMU A O, BRENNAN-CRADDOCK J P J, PANESAR A, et
al. A voxel-based method of constructing and skinning conformal
and functionally graded lattice structures suitable for additive
manufacturing[J]. Additive Manufacturing, 2017, 13: 1-13.
[2]WEEGER O, BODDETI N，YEUNG S-K，et al. Digital
design and nonlinear simulation for additive manufacturing of
soft lattice structures[J]. Additive Manufacturing, 2019, 25: 39-49. [3]DAYNEYS S, FEIH S, LU W F, et al. Design concepts for
generating optimised lattice structures aligned with strain
trajectories[J]. Computer Methods in Applied Mechanics and
Engineering, 2019, 354: 689-705.
[4]GOEL A, ANAND S. Design of functionally graded lattice
structures using B-splines for additive manufacturing[J].
Procedia Manufacturing, 2019, 34: 655-665.
[5]DUMAS M, TERRIAULT P, BRAILOVSKI V. Modelling and
characterization of a porosity graded lattice structure for
additively manufactured biomaterials[J]. Materials & Design,
2017, 121: 383-392.
[6]CHRISTENSEN J, GARCIA DE ABAJO F J. Anisotropic
metamaterials for full control of acoustic waves[J]. Physical
Review Letters, 2012, 108(12): 124301.
[7]LEVY U, ABASHIN M, IKEDA K, et al. Inhomogenous
dielectric metamaterials with space-variant polarizability[J].
Physical Review Letters, 2007, 98(24): 243901.
[8]FANG N, XI D J, XU J Y, et al. Ultrasonic metamaterials with
negative modulus[J]. Nature Materials, 2006, 5(6): 452-456. [9]YAN C, HAO L, HUSSEIN A, et al. Microstructure and
mechanical properties of aluminium alloy cellular lattice
structures manufactured by direct metal laser sintering[J].
Materials Ence & Engineering: A, 2015, 628: 238-246.
[10]ZHOU H, ZHANG X Y, ZENG H Z, et al. Lightweight
structure of a phase-change thermal controller based on lattice
cells manufactured by SLM[J]. Chinese Journal of Aeronautics, 2019, 32(7): 1727-1732.
[11]ZARGARIAN A, ESFAHANIAN M, KADKHODAPOUR J, et al.
On the fatigue behavior of additive manufactured lattice structures[J].
Theoretical and Applied Fracture Mechanics, 2019, 100: 225-232. [12]WANG H Q, CHEN Y, ROSEN D W . A Hybrid geometric
modeling method for large scale conformal cellular structures[C]//
ASME 2005 International Design Engineering Technical Conferences
and Computers and Information in Engineering Conference. New
Y ork: American Society of Mechanical Engineers, 2005: 421-427. [13]NGUYEN J, PARK S I, ROSEN D. Heuristic optimization
method for cellular structure design of light weight components[J]. International Journal of Precision Engineering &
Manufacturing, 2013, 14(6): 1071-1078.
[14]李日福, 李晋芳. 改进的三维离散点集凸包求解算法[J]. 现
代计算机, 2017(20): 40-43, 48.
LI R F, LI J F. Improved algorithm for solving convex hull of
three-dimensional discrete point set[J]. Modern Computer.
2017(20): 40-43, 48 (in Chinese).
[15]CLARKSEN K L, SHOR P W. Algorithms for diametral pairs
and convex hulls that are optimal, randomized, and incremental[C]// Proceedings of the 4th Annual Symposium on
Computational Geometry. New York: ACM Press, 1988: 12-17.
[16]BARBER C B, DOBKIN D P, HUHDANPAA H. The
quickhull algorithm for convex hulls[J]. ACM Transactions on
Mathematical Software,1996, 22(4): 469-483.。