细胞数量与排列对胞状结构压缩行为的影响

《Applied Mechanics》：Effect of Cell Number and Arrangement on the Compressive Behavior of Cellular Structures

　　

摘要：胞状结构的力学响应不仅取决于相对密度和平均胞元几何形状，还取决于胞元的空间排列。然而，排列依赖性效应如何随胞元数量的增加而演变尚未得到系统阐明。本研究采用有限元分析，在动态平衡压缩条件下，研究不同胞元数量的闭孔胞状结构的压缩行为，同时保持恒定的相对密度和相同的材料参数。使用分层泊松盘采样结合Voronoi胞化生成胞状模型。通过三种不同方式增加胞元数量：镜像复制参考结构、扩大整体试样尺寸、以及固定外部尺寸下细化胞元尺寸。为表征排列依赖性效应，引入压缩响应的两个不同特征：平均化（定义为不同初始胞元排列下响应差异的减小）和平滑化（定义为单个响应内急剧应力波动的抑制）。采用定量指标评估两种效应。在z方向压缩的板型模型和固定尺寸模型中观察到平均化，而镜像连接模型即使在大量胞元下仍保持强烈的排列依赖性。平滑化主要发生在z方向压缩的板型模型中，且与沿压缩方向排列的胞元层数密切相关，而非仅与总胞元数相关。模拟在动态平衡状态下进行，变形过程中实现了内部应力平衡。这些结果表明，压缩行为不仅受胞元数量控制，还受结构排列和定向胞元层排列的影响，为动态压缩下从排列依赖性波动到稳定宏观响应的转变提供了机理见解。

**论文解读：胞元数量与排列对闭孔胞状结构压缩行为的影响机制研究**

**研究背景与问题**
胞状结构（Cellular Structures）因其高比刚度、高能量吸收能力和设计灵活性，广泛应用于缓冲、包装、轻质面板及冲击防护等工程领域［1-6］。经典Gibson-Ashby模型指出，胞状固体的宏观刚度、强度和塌缩行为主要由相对密度和理想胞元几何决定［1-3］。然而，实际胞状材料因制造工艺、缺陷和设计约束，常呈现显著的胞元形状、尺寸和空间排列不规则性［7-9］。研究人员已通过随机Voronoi模型、随机晶格及图像重建等方法发现，结构无序会显著影响变形模式、局部化模式及压缩应力-应变响应［7-12］。近期研究进一步表明，多胞和层级胞状结构的压缩性能受截面拓扑、泡沫填充策略及多层级设计概念的强烈影响［13-17］。尽管如此，许多研究假设分析域内胞元数量足够大，局部不规则性在统计意义上被平均化，从而将胞状固体视为均匀材料［18-20］。然而，当胞元数量有限，或结构通过周期性复制生成时，这一假设的有效性并不显然［19-21］。胞元数量有限导致的尺寸效应及其与空间分布的关系尚未完全阐明［23-24］。同时，无序胞状结构压缩响应中不同来源的波动——不同初始排列间的统计散差（平均化）与单个响应内的应力突变（平滑化）——通常未被明确区分，尽管可能源于不同的结构机制［30-33］。因此，本研究旨在系统阐明闭孔胞状结构中，与初始胞元排列相关的压缩行为如何随胞元数量增加而演变，并区分平均化与平滑化两种特征。

**研究方法与关键技术**
研究人员采用几何控制的数值框架，在恒定相对密度和相同材料参数下系统研究结构效应。主要关键技术包括：
1. **分层泊松盘采样（Hierarchical Poisson Disk Sampling, HPDS）**：在三维域内随机分布胞核，同时强制相邻核间最小距离，控制胞元尺寸变异和空间随机性。
2. **Voronoi胞化（Voronoi Tessellation）**：基于胞核生成不规则闭孔几何，胞壁直接由壳单元模拟。
3. **三种胞元数增加策略**：镜像连接（Mirror-Connected，周期性复制参考结构）、板型扩大（Plate-Type，扩大整体尺寸同时随机生成胞核）、固定尺寸细化（Fixed-Size，固定外尺寸下减小胞元尺寸）。
4. **动态平衡有限元分析**：采用显式求解器RADIOSS，施加4 m/s恒定速度压缩，确保内部应力平衡；材料选用Johnson-Cook本构模型（忽略率和热效应），参数基于前期闭孔聚乙烯薄膜研究［31-33,36］。
5. **定量评估指标**：平均化由不同初始排列间平均流动应力的变异系数（CoV_σ）量化；平滑化由单个应力-应变响应内二阶差分的归一化均方根（nRMS）量化；变形局部化由等效应变分布的变异系数（CoV_ε）量化。

**研究结果**
**3.1 压缩行为评估准则**
引入平均化（不同初始排列响应差异减小）和平滑化（单个响应内应力波动抑制）两个概念，并分别以CoV_σ和nRMS定量评估。

**3.2 初始胞元排列规律性的影响**
- **平均化与平滑化**：单位胞模型（50胞元）显示强烈排列依赖性和应力波动。镜像连接模型（800胞元）并未显著降低排列依赖性或波动。z方向压缩的板型模型（800胞元）显示出最小的nRMS（1.04×10^-3，较单位胞模型降低85%）和较低的CoV_σ（1.68%），表明平均化和平滑化同时发生。固定尺寸模型（800胞元）具有最小的CoV_σ（1.01%，强平均化），但nRMS为1.70×10^-3，应力波动仍明显（仅平均化）。y方向压缩的板型模型则呈现较高的变异性（CoV_σ=4.56%，nRMS=2.30×10^-3）。
- **变形局部化机制**：z方向板型模型的等效应变分布变异系数（CoV_ε=0.351）较单位胞模型（0.526）降低33%，表明变形更均匀分布；其他模型未出现类似降低。

**3.3 胞元数量增加的影响**
- 对于z方向板型模型，50胞元平均流动应力较低，100胞元及以上平均响应趋于稳定。
- 固定尺寸模型中，随胞元数增加，平均流动应力上升，密实化起始应变下降，这综合反映了胞元尺寸细化、壁厚/胞元尺寸比变化及结构排列的影响。

**3.4 压缩方向胞元层数的影响**
z方向板型模型沿压缩方向胞元层数（N_layer≈3.68）较少，应力平滑显著；而y方向板型模型（N_layer≈14.7）和固定尺寸模型（N_layer≈9.28）层数较多，波动持续。这表明平滑化主要依赖于沿加载方向胞元层数，而非总胞元数。z方向下胞元主要在垂直于压缩方向增殖，局部塌缩事件横向更同步分布，抑制波动；y方向下胞元沿加载方向堆叠，导致逐层顺序塌缩和持续波动。

**讨论与结论**
**讨论**
- 平均化与平滑化的结构机制不同：平均化受异质变形空间分散的控制（如z方向板型和固定尺寸模型），而平滑化取决于局部塌缩事件沿压缩方向的分布（依赖于沿加载方向层数）。
- 镜像连接模型因周期性复制保留了参考结构的重复变形模式，导致排列依赖性和波动持续，表明周期性复制不足以保证统计稳定。
- 数值鲁棒性验证：采用全积分壳单元和降低加载速率（4 m/s→0.4 m/s）的结果趋势一致，表明主要趋势不受数值假象支配。
- 实验支持：使用闭孔铝泡沫（ALPORAS）进行压缩试验，保持压缩方向厚度10 mm，增加横向胞元数（约60、240、540个），测量nRMS从1.04×10^-3降至7.14×10^-4，定性支持横向胞元增加对应力波动的抑制，与z方向板型模型趋势一致。
- 建模局限：Voronoi尖锐节段、壁厚/胞元尺寸比耦合、相同材料参数等理想化假设，限制了对真实材料的直接推广。

**结论**
（保留论文结论原文翻译）
本研究通过系统改变胞元数量、空间排列以及沿压缩方向对齐的胞元层数，在恒定相对密度和相同材料参数下，研究了胞状结构的压缩行为。
结果表明，宏观压缩响应不仅取决于总胞元数量，还取决于胞元的空间组织。当异质胞元级变形在结构内空间分布时，压缩响应相对于初始胞元排列发生平均化。相反，通过周期性镜像复制增加胞元数量并不一定能降低排列依赖性，表明仅靠结构复制不足以实现统计稳定。
研究明确了应力-应变响应中平均化与平滑化之间的区别。平均化降低了不同初始配置响应间的变异性，而平滑化抑制了单个响应内的急剧应力波动。这两种现象源于不同的结构机制，并不一定同时发生。
平滑化主要在z方向压缩的板型模型中观察到，并且与沿压缩方向对齐的胞元层数密切相关，而非仅与总胞元数相关。这一发现为定向胞元层排列如何控制平台稳定性提供了结构解释。
在本研究采用的几何控制数值框架内，结果表明所观察到的趋势主要受结构特征控制，而非数值假象。进一步研究纳入改良的节段几何、胞元尺寸与壁厚的独立控制以及更广泛的实验验证，将有助于评估所识别机制的普遍适用性。