恒定平均曲率超材料的拓扑杂化以实现终极结构效率

《Materials Horizons》：Topological hybridization of constant-mean-curvature metamaterials for ultimate structural efficiencyOpen Access

　　

架构材料（Architected materials）通过梁、板或极小曲面壳获得刚度，但梁点阵因弯曲和节点应力集中而效率降低，闭孔点阵则难以制造且难以清除粉末或树脂。开孔极小曲面壳易于打印，但其零平均曲率抑制了静水载荷传递所需的拉伸变形。研究人员将平均曲率作为核心设计参数，并通过杂化（hybridization）抵消任何单一立方恒定平均曲率家族的限制。在此，研究人员引入了一种开孔杂化壳超材料，结合了四种具有不同壁跨度和非冗余载荷路径的家族：Schwarz Primitive (SP), Schwarz Diamond (SD), Schoen I-graph wrapped polyhedron (IWP) 和 Schoen F-graph rhombic dodecahedron (FRD)，其非零曲率通过拉伸而非弯曲承载载荷。通过116次高保真模拟的贝叶斯替代搜索（Bayesian surrogate search）确定了最优杂化体。在体积分数0.45时，该杂化体达到了多孔固体各向同性Hashin–Shtrikman体模量参考的97.3%，并在0.25 ≤ v ≤ 0.55范围内保持在7%以内，仅具有弱立方各向异性。通过桌面熔融沉积成型（fused-deposition modeling, FDM）打印的碳纤维增强尼龙（carbon-fiber-reinforced nylon, CF-nylon）试样再现了排名，并恢复了相应试样尺度预测的82–89%，刚度在28.8–35.5 GPa之间，峰值压应力在205–231 MPa之间（跨打印方向与面内方向），分别比打印方向上最强的单一家族基线高40%和48%。将静态有限元集成到虚拟现实头戴设备和智能手表外壳中，表明在匹配刚度或强度下可实现45%和63.5%的质量减少。因此，曲率引导的杂化提供了一条可制造的路径，使开孔壳超材料接近所声明的各向同性Hashin–Shtrikman静水参考。

**研究背景、问题与意义**

传统架构材料（Architected materials）主要通过梁点阵、板点阵或极小曲面壳获得刚度，但梁点阵因弯曲和节点应力集中而效率低下，闭孔点阵则难以制造且不易清除粉末或树脂。开孔极小曲面壳（如三周期极小曲面，TPMS）虽易于打印，但其零平均曲率约束抑制了静水载荷传递所需的膜拉伸（membrane stretching）参与，导致理论刚度上限难以达到。因此，研究人员提出将平均曲率（mean curvature, H）作为核心设计参数，通过拓扑杂化（topological hybridization）引入非零平均曲率，以激活静水膜通道，从而实现接近理论Hashin–Shtrikman（HS）各向同性体模量上限的结构效率。这项研究旨在开发一种开孔壳超材料，通过组合四种具有不同壁跨度和非冗余载荷路径的恒定平均曲率（Constant-Mean-Curvature, CMC）家族（Schwarz Primitive (SP)、Schwarz Diamond (SD)、Schoen I-graph wrapped polyhedron (IWP) 和 Schoen F-graph rhombic dodecahedron (FRD)），实现极致结构效率。论文发表在《Materials Horizons》。

**主要技术方法（不超过250字）**

研究人员采用了多尺度框架：利用Surface Evolver生成CMC壳表面，通过结构基因组力学（Mechanics of Structure Genome, MSG）进行周期性均匀化（homogenization）与去均匀化（dehomogenization），提取有效本构张量。基于贝叶斯神经网络（Bayesian Neural Network, BNN）的替代模型（surrogate model）结合期望改进（Expected Improvement）获取函数，对116个高保真MSG模拟点进行三阶段优化，搜索最优杂化组成与曲率参数。实验验证采用熔融沉积成型（FDM）打印碳纤维增强尼龙（CF-nylon）立方体试样，进行准静态压缩测试，并结合微计算机断层扫描（Micro-CT）表征几何保真度。产品级集成演示通过nTopology进行共形映射，在ANSYS中进行静态有限元分析。

**研究结果**

**2.1 CMC家族优化与静水界限接近度评估**：通过MSG均匀化，研究人员发现各单一CMC家族在体模量K*方面均低于HS各向同性上限。在体积分数v=0.45时，杂化CMC达到η=K*/K⁺_HS=0.973，并在0.25≤v≤0.55区间内保持η>0.93，且均匀化张量呈弱立方各向异性（Zener比A_Z=1.02–1.08，通用各向异性指数A_U<0.03）。

**2.2 曲率激活的杂化CMC壳载荷重分布**：通过去均匀化von Mises应力场分析，杂化壳的高应力面积分数f_high（应力≥0.8σ_app,pk,Z的壳面积比）在v=0.45时达0.38±0.03，显著高于最优传统CMC（IWP的0.21±0.02）。膜应变能分数f_m从TPMS的0.22升至杂化的0.84，表明曲率激活的膜-弯曲耦合机制使载荷通过拉伸而非弯曲传递。

**2.3 替代引导的杂化与组成转变**：BNN替代搜索以116个高保真点完成优化，相比暴力搜索降低约860倍计算成本。组成地图显示，在最优v≈0.45时，四个家族均贡献非平凡权重（无单一超过35%），表明机械互补性。杂化体在v=0.45时结构轴向刚度达32.4 GPa，比最佳传统CMC（IWP_opt的23.1 GPa）高40%。

**2.4 FDM打印杂化CMC试样的构建方向压缩**：微CT显示几何偏差<0.12 mm，内部壁孔隙率2.1±0.4%。压缩测试（ASTM D695-23）结果表明，杂化体在v=0.45时构建方向表观模量E_app,Z=28.8±1.5 GPa，峰值应力σ_app,pk,Z=212±4 MPa，达到MSG预测的82–89%（Pearson r=0.998和0.996）。方向分辨测试（X、Y、Z方向）显示各向异性继承自FDM正交各向异性原料，而非几何各向异性。

**2.5 背景比刚度基准与共形集成**：比轴向刚度/ρ* =63.6 GPa cm³ g^?1，比峰值应力σ_app,pk,Z/ρ* =416.9 MPa cm³ g^?1（ρ*=0.509 g cm^?3），超过Al 7075-T6等结构合金。共形集成演示：VR头显框架在匹配弯曲刚度下实现45%质量减少，智能手表外壳在匹配峰值应力下实现63.5%质量减少。

**讨论与结论**

讨论部分指出，杂化改进并非通过增加材料，而是通过改变载荷路径——四种CMC家族具有机械互补性，替换重复壳图案可抑制同步失稳。HS接近度结果限于v=0.25–0.55中间密度窗口，低于此窗口由弯曲和局部失稳主导，高于此窗口多模特征丧失。近各向异性被量化，HS比较作为静水参考而非严格界限。三种各向异性来源（架构、原料、材料）需区分。组成地图显示所有四个家族在最优处非平凡贡献，且Neovius冗余探测显示杂化增益接近饱和。框架可扩展至其他制造路线和力学目标。

结论部分翻译：**恒定平均曲率开孔壳超材料的拓扑杂化，在v=0.45时达到了基于密实控制压缩各向同性等效构建的各向同性Hashin–Shtrikman静水参考的97.3%，并在0.25≤v≤0.55范围内保持在7%以内，通过FDM打印CF-nylon立方体的构建方向压缩验证，恢复率为对应MSG试样尺度预测的82–89%。所有实验加载均为单轴，因此试样尺度预测得到实验验证，而静水界限接近度结果仍为计算值。界限接近度比η=K*/K⁺_HS在v=0.45时达0.973，并在0.25≤v≤0.55范围内保持>0.93。均匀化响应呈弱立方各向异性（A_Z=1.02–1.08，A_U<0.03），因此HS比较在立方各向异性括号内提供静水参考。以v=0.45为基础（ρ*=0.509 g cm^?3），杂化体提供方向特定性能：****/ρ* =63.6 GPa cm³ g^?1和σ_app,pk,Z/ρ* =416.9 MPa cm³ g^?1（实验压缩，构建方向Z），同时报告了结构合金的手册E/ρ和强度/ρ作为参考；****是架构有效轴向刚度，σ_app,pk,Z/ρ*是构建方向实验强度。将杂化体共形集成到产品级外壳中的示例性演示展示了几何可转移性，在简化静态加载下实现高达63.5%的质量减少。演示几何体的服役认证不在当前范围。该框架通过三个优化阶段的116次高保真评估达到这些结果，提供了一条可重复的多尺度路径，用于实现接近静水界限效率的开孔壳超材料，并且可直接扩展至其他增材制造路径（通过重新校准正交各向异性原料输入）以及其他力学响应目标（如延性、能量吸收、疲劳寿命和损伤容限，通过重新制定BNN标量目标）。**