载荷路径几何决定有限逾渗晶格刚度：最短边界跨越路径的普适标度关系

《Mechanics of Materials》：Load-path geometry controls stiffness in finite percolating lattices

　　

在此，研究人员采用基于梁的有限元（finite element, FE）模拟，研究通过键逾渗（bond percolation）生成的有限逾渗晶格中的刚度。对有序和几何无序晶格族在广泛连通性和无序水平范围内进行分析。结果揭示，刚度主要受最短边界跨越载荷路径（load path）的几何控制。近乎笔直的路径促进轴向主导（axial-dominated）变形和高刚度，而愈发曲折的路径诱发弯曲主导（bending-dominated）响应并大幅降低刚度。值得注意的是，当将多个晶格族的刚度数据表示为归一化最短路径长度的函数时，这些数据塌缩到单一幂律关系上，表明刚度对载荷路径曲折度（tortuosity）具有普适几何依赖性。与有序晶格相比，几何无序晶格还表现出降低的刚度变异性，表明其对局部连通性波动具有增强的鲁棒性。该工作确立最短载荷路径作为一种简单且普适的几何描述符，用于预测有限逾渗晶格中的刚度，为无序结构材料（architected materials）的分析和设计提供统一框架。

论文解读：载荷路径几何控制有限逾渗晶格的刚度

研究背景与问题。增材制造（additive manufacturing）推动结构晶格材料（architected lattice materials）和机械超材料（mechanical metamaterials）发展。蜂窝固体力学常用Gibson–Ashby框架解释，弹性模量与相对密度呈幂律关系，取决于轴向拉伸或弯曲主导。周期晶格通过规则结构提供可预测响应。近期研究显示层级梁结构可调控轴向、弯曲和剪切刚度，随机或非周期网络可增强鲁棒性和损伤容限。早期逾渗（percolation）研究多考虑无限网络和集合平均，力学性质由统计连通性控制，而非具体结构实现。实际结构有限，且制造或设计扰动带来几何变异，响应可能强烈依赖局部拓扑和承载路径几何。逾渗理论描述随机键或位点占据中宏观连通涌现，经典连通逾渗在临界键概率（bond probability）p_c形成系统跨越簇（system-spanning cluster）。但几何连通不保证机械稳定；中心力网络存在刚性转变（rigidity transition），含弯曲或转动相互作用时拉伸与弯曲贡献关键。骨架（backbone）和弹性骨架（elastic backbone）等承载子网络控制力学响应，但最短跨越载荷路径几何与有限结构晶格刚度间定量关系尚未建立。近临界理论针对无限统计网络，不描述远离临界点的有限晶格。因此，研究人员在键逾渗框架内比较有序与几何无序晶格，试图建立载荷路径几何与刚度的直接联系。

研究人员开展的研究。研究人员采用基于梁的有限元分析，研究二维键逾渗生成的有序和几何无序晶格，覆盖四种晶格族：正方形、中心正方形、三角形和六边形。通过控制键概率p∈[p_c,1]保留键，并用高斯抖动扰动节点坐标引入几何无序（geometric disorder），σ=0.3、0.6、0.9 mm，同时保持键占据统计和连通性。每个键概率生成100个独立实现。分析前提取边界跨越骨架，去除悬垂枝和不承载分支。用二维欧拉–伯努利梁（Euler–Bernoulli beam）单元建模，节点自由度包括水平位移u_x、竖向位移u_y和转角θ。底边固定，顶边施加竖向位移u_y=-L₀，由总反力除以施加位移定义有效刚度。承载构件通过总应力σ_total=|σ_axial|+|σ_bending|识别，阈值α=0.1 max σ_total，并用弯曲比率（bending ratio）γ=σ_bending/σ_total量化弯曲贡献，γ=0为纯轴向，γ=1为纯弯曲。模型用Abaqus B23验证。关键方法可概括为：二维键逾渗晶格生成与高斯几何扰动；边界跨越骨架提取；二维梁基有限元力学分析；承载构件应力参与阈值与弯曲比率计算。样本队列来源为计算生成晶格，无生物样本队列。

研究结果。3.1 刚度随键概率变化：研究人员通过有限元计算归一化刚度K/K_ord^max随p≥p_c的变化。有序晶格在p=1时收敛为单一值；无序晶格在所有键概率下持续存在刚度变异。几何无序一般降低归一化刚度，σ=0.3 mm接近有序，σ=0.9 mm最低。刚度波动方面，高键概率时几何无序大的晶格波动更大；接近逾渗阈值时趋势反转，高度无序晶格变异更低，而有序晶格变异更大。原因是无序网络存在多个替代载荷路径，可重新分配应力，降低对单个承载构件去除的敏感性；有序网络对关键边去除敏感。

3.2 最短载荷路径几何与晶格刚度：研究人员比较同一键概率下不同实现的承载路径。近乎笔直、与加载方向一致的最短边界跨越路径促进轴向主导变形，刚度高；高度曲折路径诱发弯曲主导变形，刚度显著降低，差异约三个数量级。将归一化刚度K/K_max对归一化最短路径长度L_short/L_Ord^short表达后，不同晶格族和几何无序水平的数据塌缩到经验标度关系K/K_max ∝ (L_short/L_Ord^short)^?n，指数n≈22–26。归一化最短路径长度因此可作为有限逾渗晶格刚度的稳健几何描述符。

3.3 变形机制与鲁棒性：研究人员分析变形模式与刚度关系。3.3.1 弯曲比率与刚度：高刚度时有序和无序晶格弯曲比率γ均较小，载荷主要经轴向拉伸传递；有序晶格γ略低于无序晶格，因其载荷路径更对齐。刚度降低时γ增大，发生从轴向主导到弯曲主导的转变，该转变与最短跨越路径曲折度增加相关。这些结果说明最短跨越路径几何同时控制变形模式和刚度，并使无序晶格在连通性波动下更鲁棒。

讨论总结与结论翻译。讨论表明，有限逾渗晶格刚度主要受最短承载路径几何控制。近直、与加载方向对齐的系统跨越路径通过轴向变形传力，刚度高；越来越曲折的路径促进弯曲主导，即使全局连通保持，刚度也显著降低。归一化刚度与归一化最短路径长度之间存在稳健经验标度关系，说明刚度对载荷路径几何高度敏感。归一化最短路径长度可作为所研究有限晶格的归一化刚度的简单几何描述符。尽管该标度关系由有限元模拟建立，但它捕捉了不同晶格拓扑和几何无序水平下的刚度趋势。该工作建立网络拓扑与力学响应的直接联系，将经典逾渗概念从连通性和刚性扩展到载荷传递的几何效率。有序与无序晶格比较进一步显示刚度与鲁棒性的权衡：有序晶格可通过高效轴向载荷传递获得高刚度，但对单个承载构件去除高度敏感；几何无序晶格因多个替代载荷路径和更分布变形，刚度变异性降低，对连通性变化更鲁棒。该框架为结构材料和机械超材料提供设计指南，可通过控制载荷路径曲折度和结构无序来定制刚度和鲁棒性。未来可扩展至三维结构、非线性变形状态和失效过程。

结论翻译：本研究表明，有限逾渗晶格的刚度主要受最短承载路径的几何控制。包含与加载方向对齐的、几乎笔直的系统跨越路径的晶格主要通过轴向变形传递力，因此表现出高刚度。相反，愈发曲折的路径促进弯曲主导变形，即使全局连通性得以保持，也会导致刚度显著降低。结果揭示了归一化刚度与归一化最短路径长度之间稳健的经验标度关系，表明刚度对载荷路径几何的强敏感性。该关系进一步表明，对于所研究的有限晶格，归一化最短路径长度可作为归一化刚度的简单几何描述符。尽管所提出的标度关系基于有限元模拟建立，但它捕捉了不同晶格拓扑和几何无序水平下观察到的刚度趋势。这些发现建立了网络拓扑与力学响应之间的直接联系，将经典逾渗概念从连通性和刚性扩展到载荷传递的几何效率。有序与无序晶格之间的比较进一步突显了刚度与鲁棒性之间的相互作用。有序晶格可通过高效轴向载荷传递实现高刚度，但其响应对单个承载构件的去除仍然高度敏感。相反，由于存在多个替代载荷路径和更分布的变形模式，几何无序晶格表现出降低的刚度变异性和对连通性变化更强的鲁棒性。总之，这些结果确立最短路径几何作为有限结构晶格力学性能的预测性微结构描述符。与以往主要考察晶格拓扑、梁结构或几何无序影响的研究不同，本工作将问题置于键逾渗框架内，并表明载荷路径几何提供了一个统一描述符，控制不同晶格族的刚度。该贡献表明，对于有限逾渗结构晶格，载荷路径几何提供了一个具有物理意义的描述符，将微结构拓扑与宏观刚度联系起来。本框架为结构材料和机械超材料提供设计指南，其中可通过受控操控载荷路径曲折度和结构无序来定制刚度和鲁棒性。将该方法扩展到三维结构、非线性变形状态和失效过程是未来工作的一个有前景方向。