1、收稿日期：2023-02-07基金项目：国家杰出青年科学基金（52025066）资助作者简介：闫广涵（1997），男，在读博士研究生。引用格式：闫广涵，严晗，姜楠，等.新型航空发动机空气冷却器流动传热性能试验J.航空发动机，2023，49（4）：48-53.YAN Guanghan，YAN Han，JIANG Nan，et al.Flow and heat transfer performance test of a new type air cooler for aeroengineJ.Aeroengine，2023，49（4）：48-53.新型航空发动机空气冷却器流动传热性能试验闫广涵1，
2、严晗1，姜楠2，刘瑜1，梁义强2，赵佳飞1（1.大连理工大学 能源与动力学院，辽宁大连 116024；2.中国航发沈阳发动机研究所，沈阳 110015）摘要：空气冷却器是航空发动机中空气热量交换的重要场所，针对传统扰动结构强化换热不足问题，重点研究了空气冷却器中扰动结构的结构形貌特征与流动换热的关系。以三周期极小曲面结构为基础，设计制造了D型极小曲面扰动结构、P型极小曲面扰动结构和基于Sigmoid函数杂化方法的D-P型（P-D型）结构。并利用纳米CT重构出样品的3维数字特征，试验研究了不同结构来流速度与压降、换热性能和综合换热性能的关系。结果表明：D型极小曲面结构具有最优的压降、换热性能和综
3、合换热性能：D-P型杂化极小曲面结构的阻力系数f最大，比D型结构的高56.81%；D型的平均努塞尔数Nu比P型的高18.49%，P-D和D-P杂化结构的换热性能比P型结构的有效提高；P-D型D-P型杂化结构的Nu分别比P型结构的高2%和8.27%。在4种结构中，D型结构的综合换热性能最优，分别比P型、P-D型、D-P型结构的高40.35%、57.2%、71.02%。关键词：空气冷却器；三周期极小曲面；流动换热；结构杂化；航空发动机中图分类号：V233.5文献标识码：Adoi：10.13477/ki.aeroengine.2023.04.007Flow and Heat Transfer Per
4、formance Test of a New Type Air Cooler for AeroengineYAN Guang-han1，YAN Han1，JIANG Nan2，LIU Yu1，LIANG Yi-qiang2，ZHAO Jia-fei1（1.School of Energy and Power Engineering，Dalian University of Technology，Dalian Liaoning 116024，China；2.AECC Shenyang Engine Research Institute，Shenyang 110015，China）Abstract
5、：An air cooler is an important place for air heat exchange in aeroengine.In response to the problem of insufficient heat transfer enhancement by traditional disturbedstructures,this paper focuses on the relationship between the structure morphology characteristicsof disturbance structures and the fl
6、ow and heat transfer performances of the air cooler.Based on the triply periodic minimal surface structure,the D-type and P-type triply periodic minimal surface structure,and a Sigmoid function-based hybrid structure P-D type（D-P type）were designed and manufactured.The digital characteristics of the
7、 samples were reconstructed by nano-CT.For these structures,the relationships between the flow velocity and pressure drop,heat transfer performance,and overall heat transfer performance were experimentallystudied.The results show that the D-type triply periodic minimal surface structure has the best
8、 pressure drop,heat transfer performance,and overall heat transfer performance.The D-P type hybrid triply periodic minimal surface structure has the largest drag coefficient（f）,which is 56.81%higher than that of the D-type structure.The average Nusselt number（Nu）of the D-type is 18.49%higher than th
9、at of theP-type.The heat transfer performances of the P-D type and the D-P type hybrid structure are effectively improved compared to the P-typestructure.The Nu of the P-D type and D-P type hybrid structures are 2%and 8.27%higher than that of the P-type structure,respectively.Among the four structur
10、es,the D-type structure has the best overall heat transfer performance,which is 40.35%,57.2%,and 71.02%higherthan that of the P-type,P-D type,and D-P type respectively.Key words：air cooler;triply periodic minimal surface;flow and heat transfer performance;structural hybridization;aeroengine第 49 卷 第
11、4 期2023 年 8 月Vol.49 No.4Aug.2023航空发动机Aeroengine0引言航空发动机的更新换代带来的温度提升问题是制约飞机进一步发展的重要原因。目前增压比接近30，涡轮前温度将达到2400 K左右1。随着发动机推重比的不断提高，未来推重比15一级的第5代发动机涡轮进口温度将达到20002250 K2。对发动机的压气机进口进行预先冷却是降低发动机温度的途径，可闫广涵等：新型航空发动机空气冷却器流动传热性能试验第 4 期用热沉有燃油、外涵道空气等3。在冷气用量和冷却结构无法大幅变化的前提下，最有效的出路就是降低冷却空气本身的温度，提高冷气本身品质。这种技术被称为冷却冷却空
12、气（Cooled Cooling Air，CCA）技术4-6。在发动机上安装空气冷却器，使用外涵冷气或燃油对冷却空气预先冷却，从而提高冷却空气品质3。换热器是热量交换的重要场所，高温空气与冷却工质在换热器内通过换热面充分传热7-8。用于冷却空气的换热器称为空气冷却器，由于航空发动机应用环境的复杂性，对其空气冷却器提出了要求3：轻质化、高效、安全；在系统中布置合理，换热性能和空气动力学性能良好。在空气冷却器效率的影响因素中，涡量是决定换热是否充分的重要指标。影响涡量的因素有流量以及流道中扰动结构的形貌特征。多孔介质因为比表面积大、质量轻等优点越来越受到关注，成为替代传统翅片的理想扰动结构9。流体
13、流经多孔介质扰动结构时产生大量的涡，有利于热量充分交换，但是会提升压降，加大泵功消耗。合理选择多孔介质扰动结构可以提升综合换热性能。现阶段在换热器上常用的多孔介质有：桁架10、泡沫金属11、点阵结构12以及三周期极小曲面结构等。三周期极小曲面结构，因其比表面积大、力学性能可靠、结构形貌特征更容易控制等优点被广泛研究。三周期极小曲面扰动结构目前在换热器领域应用广泛，因为其通过数学公式可方便地控制孔隙率、整体形状等结构参数，并可以根据需要调整。Kim等13设计出多种以极小曲面为核心的异形换热器，以满足不同应用需求。将极小曲面结构作为一次表面换热器也是当前的研究热点，Femmer T等14的1项试验
14、研究中微型换热器核心结构总体尺寸是 10 mm10 mm10 mm，冷热流体在换热器中换热。结果表明，随着流量的增加，所有换热器的传热效率都降低文献中提到的传热效率是1个无量纲传热效率。Al-Ketan O等15对于结构优化研究已经有梯度孔隙率的相关内容，杂化是结构优化的方法之一，在结构杂化相关研究中，以往的研究更偏重力学性能分析，对于应用在换热器中的流动换热性能，有待进一步研究。本文通过试验方法，研究不同结构的流动传热特性，为三周期极小曲面杂化结构进一步应用提供数据支撑。1模型设计与制造1.1模型设计三周期极小曲面控制方程以三角函数为基础，不同类型的极小曲面通过不同的控制方程表达，本文研究的
15、三周期极小曲面基本类型是D型和P型，基本表达式为f（x,y,z）D=sin Xsin Ysin Z+sin Xcos Ycos Z+cos Xsin Ycos Z+cos Xcos Ysin Z（1）f（x,y,z）D=cos X+cos Y+cos Z（2）式中：X=2x，Y=2y，Z=2z；令 f（x,y,z）=c，c为常数，代表垂直于曲面偏移一定距离，偏移后再闭合就形成了带厚度的极小曲面结构，当c=0时，代表没有厚度的极小曲面。杂化Sigmoid函数为fSigmoid=11+exp(kx)（3）式中：k为过渡的平滑程度，k值越小过渡越平缓，在这里k=10，此时杂化区域刚好占据1个单元胞体
16、的长度。杂化方程为fhybrid（D/P）=fSigmoidfD+（1-fSigmoid）fP（4）在本文的研究中选用片状极小曲面生成策略-cf(x,y,z)c，调整c值，使各模型的孔隙率均为80%，骨架体积 20%。模型的参数化表征见表 1。通过开源软件 MSLattice17生成 STL 格式的 3 维文件。杂化理论在本文的研究中选用Sigmoid函数，Sigmoid函数杂化的特点是能够在非杂化区最大限度地保存原始结构的形貌特征，在杂化区实现平滑过渡，平滑过渡对于流动性能具有优势。3种三周期极小曲面的3维模型如图1所示。由于原始三周期极小曲面结构在3个方向上都是周期延申的，因此D型结构和P
17、型结构不区分不区分流动方向，而杂化的DP结构中间插入杂化区域，来模型D型P型DP杂化宏观尺寸/mm505010505010505010胞元尺寸/mm101010101010101010胞元数252510（D）+10（P）+5（杂化）骨架体积/mm3500050005000比表面积/（1/m）399025203264表1模型参数化表征49航空发动机第 49 卷流方向首先流经D区域再流经 P 区域和首先流经 P区域再流经D区域是有区别的，杂化结构的 2 种不同流动方向如图2所示。1.2模型制造由于三周期极小曲面结构的复杂结构特征，传统的减材制造无法满足结构制造需求。随着增材制造技术的快速发展，制造
18、出结构复杂的三周期极小曲面结构的技术已经越来越成熟。金属3D打印机选用EP-M250，材料选用AlSi10Mg粉末，密度为2560 kg/m3，金属粉末粒度直径分布在50 mm左右。D型结构的样品实物如图3所示。为了进一步观察打印件的形貌特征，使用X射线扫描样品结构，扫描设备是Bruker SkyScan 2214纳米CT。扫描完成后重构结构，通过后处理软件渲染样件特征，P型结构渲染如图4所示。从图中可见，在表面处有不均匀的粗糙度，这是因为增材制造过程中颗粒融化后在表面堆积，形成了并不光滑的面，粗糙度不均匀是因为增材制造的后处理不能使表面全部涉及。2试验2.1试验原理设计单侧空气对流换热试验，
19、测量进出口温度，加热底面温度以及进出口压降，整个试验段的流量也使用体积流量计实时监测。试验以及数据处理过程中用到如下公式h=qTw-Tin，j=NuRePr13，f=PL2dhu2，Nu=hdh，Re=udhv，Pr=cp，dh=4Sr（5）式中：q 为热通量；Tw为底面平均温度；Tin为入口温度；j为无量纲的表面传热系数；f为阻力系数；P为入口和出口之间的压降；L为样品流向长度；dh为水力直径；为空气密度；为空气的热导率；u为入口速度；为空气的动力粘度；为空气的运动粘度；cp为空气在恒压下的比热；S 为流道横截面积；r 为湿周长度；Nu为努塞尔数；Re为雷诺数；Pr为普朗特数。试验台原理如图
20、5所示。试验装置由供气系统、加热系统、数据采集系统和试验段组成。试验装置如图6所示。试验仪器不确定度见表2。图2杂化结构的2种不同流动方向D-P流动方向P-D流动方向图3D型结构样品实物图4P型结构渲染图5试验台原理（a）供气系统中的风机和变频器（b）供热系统中的直流电源和电加热棒图6试验装置（c）数据采集系统中的数据记录仪、热电偶、压力传感器、涡流流量计（d）试验段数据采集仪计算机压力传感器法兰涡街流量计变频器风机软管通道热电偶电加热棒电加热棒直流电源盖板样品试验段10 mm16.6 mm热电偶孔1.5 mm25 mm电加热棒孔6 mm61 mm仪器压力传感器涡流流量计不确定度0.21.5仪
21、器K型热电偶直流电源不确定度0.751表2试验仪器不确定度%图1三周期极小曲面3维模型（a）D型结构（b）P型结构（c）DP杂化结构50闫广涵等：新型航空发动机空气冷却器流动传热性能试验第 4 期供气系统由风机（HaoGuan HG620-HF-2-BP）、变频器（WeiSen AE200H-4-3PH-7.5G-11PSPD990-G2.2KW-H3-D）、多根不锈钢管组成的通道以及柔性风管组成。为了保证相对稳定的气流，在试验中选择了负压风机。加热系统由直流电源（ZhaoXin PS-6005D）和 8根电热棒组成。其中直流电源电压调节范围为060 V，电流调节范围为05 A。每根电热棒的最
22、大加热功率为60 W。为了保证壁面有均匀的热流密度，将8根电热棒均匀的布置在试验段底部。数据采集系统由3个子系统组成，分别是温度采集系统、压力采集系统和流量采集系统。温度采集系统由K型热电偶以及数据记录仪（Agilent 34972A）组成。为了在计算机上实现实时采集，在试验段的入口和出口分别布置了3个热电偶，然后在试验段底部沿流动方向对称布置了10个热电偶，并让热电偶与数据记录仪连接。压力采集系统由2个压力传感器（UNIK 5000 PTX 5072-TC-A1-CA-H0-PA）组成，2个压力传感器被分别设置在试验段的入口和出口。流量采集系统由流量计以及数据采集软件组成。将涡街流量计（LU
23、GB-2/2/03/Z/D/E/N）严格按照说明书布置在通道内。此外，在试验段前后放置相同截面形状的聚四氟乙烯材料，起到隔热和固定作用。试验测点包括温度测点、压力测点和流量测点。试验台的温度测点共16个，热电偶分别布置在试验段的进、出口（各3个）和加热底面（10个），进、出口布置的热电偶距离金属骨架进出口边缘40 mm处，3个热电偶探头距离底面的距离分别为3、5、7 mm，温度的进、出口值取3个热电偶的平均值。加热底面的10个热电偶均匀布置在距离底面0.5 mm的试验段内，热电偶从两侧均匀插入，每侧5个热电偶，每2个热电偶间距10 mm，测点深度25 mm，对侧的2个测点间距离10 mm。压力
24、传感器进出口位置各1个，距离金属骨架进、出口边缘55 mm处。流量计布置在距离试验段1600 mm处。2.2试验步骤试验操作流程如下：（1）将试件放入试验段，盖上盖板并拧紧；（2）用密封胶将试验段的缝隙密封；（3）打开风机确认有无泄露；（4）确认无泄漏后，用隔热棉将试验段做隔热处理；（5）打开各测量系统并观察物理参数是否异常；（6）确认无异常后，打开风机，调整变频器频率；（7）选取试验系统参数（见表 3）进行进一步试验，监控温度、压力和流量值，当温度在10 min内的变化不超过0.5 K后采集数据。为了减少采样误差，每个工况的采样时间为10 min。每个工况需做重复试验3次，减小操作过程中出现
25、的误差。2.3误差分析试验的不确定度由2部分构成，分别是测量误差（A类不确定度）和试验仪器的误差（B类不确定度），试验总体误差为U=U2A+U2B（6）其中UA=ti=1n(xi-x)n(n-1)，UB=instrumentC（7）式中：t为为了弥补因试验次数不够多，不能够完美地满足测量误差正态分布的规律的修正因子；instrument为指仪器的分度值；无特别说明时，置信度在 68.3%的情况下，C 值取3，这时是均匀分布。试验结果误差值见表4。3结果与讨论通过试验方法，对不同结构的流动和换热性能进行研究，本章分别讨论了来流速度与样品结构对压降、换热系数以及综合换热系数的影响。3.1压降流体在
26、管中流动时由于能量损失而引起的压力降低，本文的压降指流道中金属骨架进出口测量的空气压力差值。随着来流速度的增加，压降逐渐升高，采用无量纲数阻力系数f表征压降的大小。不同结构阻力系数f与雷诺数Re的关系如图7所示。在4种模型中，D型结构f最小，P型结构次之，相比于D型结参数体积流量/（m3/h）入口速度/（m/s）入口温度/加热功率/W数值14.428.8816室温（2427）100梯度3.62随机表3试验系统参数试验结果fNuj/f误差值/%2.41.685.4表4试验结果误差值51航空发动机第 49 卷构，P型结构在试验速度下的平均 f 比 D 型的大18.76%，说明了D型结构具有良好的流
27、动性能。杂化结构的流动性能较差，其中 D-P 型结构的 f最大，比 D 型结构的大56.81%，P-D型结构的压降比 D 型结构的大35.29%。杂化结构中2种原始结构之间的杂化过渡带一定程度上减小了流体流动方向的通孔率，对流体流动造成一定程度的堵塞因此产生了较高的阻力。3.2换热性能采用努塞尔数Nu表征对流换热能力，Nu是壁面上流体的无量纲温度梯度，值越大表示对流换热能力越强。在结构传热研究中换热性能是一大重要指标，换热性能的优劣决定了扰动结构在换热器上的应用可行性，本节研究了不同扰动结构在不同来流速度下的Nu的大小。不同结构努塞尔数Nu和雷诺数Re的关系如图8所示。从图中可见，D 型极小曲
28、面结构具有最佳的换热性能，P型极小曲面结构的换热性能最差，D 型比 P 型的平均Nu 高 18.49%。P-D 型杂化极小曲面结构和D-P型杂化极小曲面结构的 Nu介于D型结构和P型结构之间，说明杂化结构相对于P型结构的Nu有所提升，P-D型杂化结构和D-P型杂化结构分别比P型结构的Nu高2%和8.27%。这一试验结果说明了Sigmoid函数杂化方法对于提高传统三周期极小曲面扰动结构的换热性能具有积极作用，同时也要注意到也是对某些传统结构的降低，根据具体的工况选择杂化方法是工程应用的基础。3.3综合换热性能压降和传热的综合评价指标是结构优化研究中的重要参考依据，本文的综合换热性能采用无量纲数j
29、/f，表示在相同的阻力因子下传热量的一种度量，其物理意义是相同流量下换热能力的增加是否大于阻力的增加。不同结构综合换热性能j/f和雷诺数Re的关系如图9所示。从图中可见，在 4 种结构中，D 型结构的综合换热性能最优，分别比 P 型，P-D 型，D-P 型结构高出 40.35%，57.2%，71.02%。说明D型结构在j/f的评价机制下性能占优。同时，D型和P型结构相对于杂化结构综合换热性能提升说明了空间中曲面的各个曲率为0这种曲面特征在单侧流动换热结构中具有优势。在工程应用中也要兼顾具体的环境，选择合适的结构以及杂化形式。4结论（1）不同结构在不同来流速度下的压降、换热、综合换热性能试验表明
30、，D型结构f最小，P型结构的次之，相比于D型结构，P型结构在试验速度下的平均f比D型的高18.76%，说明了D型结构具有良好的流动性能。杂化结构的流动性能较差，其中D-P型结构的f最大，比D型结构的高56.81%，P-D型结构的压降比D 型结构的高35.29%。D型极小曲面结构具有最佳的换热性能，P型极小曲面结构的换热性能最差，D型的平均Nu比P型的高18.49%。P-D型杂化极小曲面结构和D-P型杂化极小曲面结构的Nu介于D型结构和P型结构的之间，说明杂化结构的Nu相对于P型结构的有所提升，P-D型杂化结构和D-P型杂化结构的Nu分别比P型结构的高2%和8.27%。在4种结构中，D型结构的综
31、合换热性能最优，分别比P型、P-D型、D-P型结构的高40.35%、57.2%、71.02%。（2）在实际工程应用中，应合理选择结构以及杂化形式，保证在换热器中插入的扰动结构能够最大限度地满足环境需求。三周期极小曲面的特点是平均曲率为0，主要表现出较好的流动性能，流体沿着壁面流动更加顺畅，对传热的增强主要体现在比表面积的增大，热量可通过更大的表面积传递。航空发动机用换热器的核心需求是安全、轻质、高效，更低的泵功传递更多的热量，体积和质量尽可能减小，以三周期极小曲面为核心的结构在逻辑上符合航空发动机用换热器强化传热的本质要求。下一步研究工作应该以深入挖掘原始结构及杂化结构的流动传热性能机理为主要
32、目标，通过 CFD图8不同结构努塞尔数Nu和雷诺数Re的关系Nu1000900800700600500Re16000120008000DPP-DD-Pj/f0.00340.00300.00260.00220.00180.0014Re16000120008000DPP-DD-P图9不同结构综合换热性能j/f和雷诺数Re的关系图7不同结构阻力系数f与雷诺数Re的关系f48444036322824Re18000160001400012000100008000DPP-DD-P52闫广涵等：新型航空发动机空气冷却器流动传热性能试验第 4 期仿真方法构建流动换热过程中的速度、压力、温度云图，直观认识到内部
33、物理场的分布，为三周期极小曲面结构的定向优化提供理论依据。参考文献：1 温奇，裴鑫岩，田鸿宇，等.空-油换热器对变循环航空发动机总体性能的影响研究J.推进技术，2023，44（2）：39-48.WEN Qi，PEI Xinyan，TIAN Hongyu，et al.Overall performanceevaluation of variable cycle aero-engine based on air-fuel heat exchanger J.Journal of Propulsion Technology，2023，44（2）：39-48.（in Chinese）2 姚艳玲，代军，黄
34、春峰.现代航空发动机温度测试技术发展综述J.航空制造技术，2015（12）：103-107.YAO Yanling，DAI Jun，HUANG Chunfeng.Development for temperature measurement technology in modern aeroengine J.AerospaceManufacturing Technology，2015（12）：103-107.（in Chinese）3 于霄，吕多，李洪莲，等.空气冷却器在航空发动机上的应用及流动传热试验分析技术研究J.计测技术，2017，37（3）：34-38.YU Xiao，LYU Duo，
35、LI Honglian，et al.Measurements of the flowand heat transfer research in an aeroengine heat exchangerJ.Metrology and Measurement Technology，2017，37（3）：34-38.（in Chinese）4 Missirlis D，Donnerhack S，Seite O，et al.Numerical developmentof a heat transfer and pressure drop porosity model for a heat exchang
36、er for aero engine applicationsJ.Applied Thermal Engineering，2010，30（11-12）：1341-1350.5 Missirlis D，Yakinthos K，Palikaras A，et al.Experimental and numerical investigation of the flow field through a heat exchanger foraero-engine applicationsJ.International Journal of Heat and FluidFlow，2005，26（3）：44
37、0-458.6 Kallath，HariharanKholi，Foster K J，et al.Numerical study of theflow uniformity inside the high-pressure side manifolds of a cooledcooling air heat exchangerJ.Applied Thermal Engineering：Design，Processes，Equipment，Economics，2021，189（1）：1-14.7 刘荫泽，张声宝，董威，等.航空发动机1次表面换热器流动换热性能分析J.航空发动机，2017，43（4）
38、：61-68.LIU Yinze，ZHANG Shengbao，DONG Wei，et al.Analysis of flowand heat transfer performance of primary surface heat exchanger foraeroengineJ.Aeroengine.2017，43（4）：61-68.（in Chinese）8 张声宝，郭之强，郑梅，等.3D打印一次表面换热器流动换热性能试验研究J.航空发动机，2018，44（5）：79-85.ZHANG Shengbao，GUO Zhiqiang，ZHENG Mei，et al.Experimental
39、study on flow and heat transfer characteristics of 3D printed primary surface heat exchangerJ.Aeroengine.2018，44（5）：79-85.（in Chinese）9 任勇翔，于霄，张筱喆，等.类多孔结构超轻高效换热器流动传热特性研究J.南京航空航天大学学报，2019，51（4）：449-455.REN Yongxiang，YU Xiao，ZHANG Xiaozhe，et al.Flow and heattransfer characteristics of a ultralighted a
40、nd efficient porous mediaheat exchangerJ.Journal of Nanjing University of Aeronautics andAstronautics，2019，51（4）：449-455.（in Chinese）10 Zhang Q，Jiang W，Zhang Y.Effect of geometrical parameters onthe effective elastic modulus for an X-type lattice truss panel structureJ.Science and Engineering of Com
41、posite Materials，2018，25（6）：1135-1144.11 Abadi G B，Moon C，Kim K C.Experimental study on single-phase heat transfer and pressure drop of refrigerants in a plate heatexchanger with metal-foam-filled channelsJ.Applied Thermal Engineering，2016，102：423-431.12 左蔚，赵剑，白静，等.激光选区熔化菱形正十二面体点阵材料的承载与失效特性J.火箭推进，20
42、20，46（5）：87-93.ZUO Wei，ZHAO Jian，BAI Jing，et al.Research on bearing anddisable properties of laser selective melting rhombohedral dodecahedral lattice materialsJ.Journal of Rocket Propulsion，2020，46（5）：87-93.（in Chinese）13 Jiho K，Dong J Y.3D printed compact heat exchangers with mathematically define
43、d core structuresJ.Journal of Computational Design and Engineering，2020（4）：4.14 Femmer T，Kuehne A J C，Wessling M.Estimation of the structure dependent performance of 3-D rapid prototyped membranesJ.Chemical Engineering Journal，2015，273：438-445.15 Al-Ketan O，Ali M，Khalil M，et al.Forced convection CFD
44、 analysis of architected and 3D printable heat sinks based on triply periodic minimal surfacesJ.Journal of Thermal Science and Engineering Applications，2020，13（2）：1-33.16 Nn A，Ak B，Mb A，et al.Development of novel hybrid TPMS cellular lattices and their mechanical characterisationJ.Journal ofMaterials Research and Technology，2021，15：1318-1329.17 OAl K，RKA A.MSLattice：a free software for generating uniformand graded lattices based on triply periodic minimal surfacesJ.Material Design&Processing Communications.2020，3（6）：1-10.（编辑：刘亮）53