《COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING》：A unified fluid-solid elasto-viscoplastic finite element model for the simulation of 3D concrete printing across process scales
编辑推荐：
针对3D混凝土打印（3DCP）过程中挤出、层沉积与可建造性等多尺度物理现象难以统一模拟的难题，本研究首次将Saramito弹黏塑性本构模型拓展至水泥基材料，并引入Drucker-Prager屈服准则以考虑压力敏感性，结合粒子有限元法（PFEM）与Delaunay重网格技术，实现从流体到固体状态平滑、热力学一致的多尺度模拟。模型通过单层/多层墙实验验证，准确预测了挤出动力学、层间融合、接头形成及整体力学响应，为3DCP的工艺优化与结构设计提供了高精度计算工具。
论文解读
在建筑行业迈向数字化、自动化的浪潮中，3D混凝土打印（3D Concrete Printing, 3DCP）技术以其高几何精度、短工期、低材料浪费和可持续性优势，正引领着一场静悄悄的革命。然而，这项看似前景无限的技术，在实践中却面临重重挑战。材料的复杂行为与工艺参数紧密耦合，对环境条件极为敏感，常常迫使工程师们依赖“试错”策略，不仅效率低下，也难以充分发挥其制造复杂优化结构的潜力。
更关键的是，3DCP过程跨越了多个空间与时间尺度：在挤出和沉积的瞬间，材料表现出如同蜂蜜般的流体特性；而在层层堆叠、静置固化后，它又需承担起整个结构的重量，展现出固体的刚性与强度。这种从“流”到“固”的华丽转身，背后是极为复杂的多物理场耦合过程。传统的数值模拟方法往往“分而治之”：基于流体的模型能精细刻画挤出形状，却难以预测整体结构的坍塌；基于固体的模型擅长评估可建造性，却又忽略了挤出过程中的局部塑性变形与实际截面形状。这种割裂，使得高保真地预测整个打印过程成为一项“尚未完成的拼图”。
为此，来自意大利米兰理工大学的研究团队Giacomo Rizzieri, Derk Bos, Rob Wolfs, Liberato Ferrara, Massimiliano Cremonesi在《Computer Methods in Applied Mechanics and Engineering》上发表了一项开创性研究。他们决心打破这种藩篱，提出了一个统一的、能够跨越整个3DCP过程尺度进行高保真模拟的先进计算框架。
核心研究方法
研究团队的核心创新在于两点。首先，在本构模型上，他们首次将描述复杂流体（如泡沫、乳液）的Saramito弹黏塑性（elasto-viscoplastic）模型拓展至水泥基材料，构建了“Saramito-DP”模型。该模型通过一个上随体导数描述的演化方程，实现了材料从类流体状态到类固体状态（neo-Hookean弹性）的平滑、热力学一致过渡。为捕捉早龄期混凝土的压力依赖性，采用Drucker-Prager（DP）屈服准则替代了原模型中的von Mises准则；同时，通过随时间演化的材料参数（如剪切模量G(t)和屈服参数k(t)），考虑了触变性和结构形成的影响。模型仅需一组有限的、可通过实验识别的参数。其次，在数值方法上，他们将新本构模型嵌入到粒子有限元法（Particle Finite Element Method, PFEM）框架中。PFEM结合了更新拉格朗日描述和高效的Delaunay重网格技术，天生擅长处理大变形、自由表面演化和拓扑变化，为统一处理流-固行为提供了连续的计算环境。此外，他们采用了分裂求解策略和CPU并行计算，大幅降低了全尺度模拟的计算成本。
研究结果与发现
1. 模型验证：几何形状与可建造性预测
研究在荷兰埃因霍温理工大学进行了系统的实验验证。
- •
层沉积几何精度：通过模拟打印单层和多层（包括垂直和倾斜）墙体，并将预测的横截面轮廓与激光扫描数据对比，发现模拟与实验数据高度吻合。分析表明，基底的弹性对层沉积行为和纤维轮廓有显著影响。
- •
可建造性与失效机制：对墙体进行了包含挤出、沉积和静置的全过程3D模拟，成功复现了实验中观察到的面外屈曲坍塌机制。这表明模型能够可靠地预测结构的可建造性极限。
2. 全尺度复杂结构模拟
研究进一步对一个具有层间融合特征的复杂结构构件进行了全尺度模拟。结果显示，模型能够捕捉：
- •
挤出动力学：材料从喷嘴挤出的流动行为。
- •
层间融合与节点形成：相邻纤维在交汇处的合并与接头形状的形成。
- •
整体力学响应：打印完成后，构件在自重及固化过程中的力学行为。
这些高保真的预测能力，为优化打印路径、评估接头质量和预测结构性能提供了强大工具。
3. 模型退化与普适性
理论分析表明，所提出的“Saramito-DP”模型具有良好的普适性。当摩擦角φ=0时，Drucker-Prager准则退化为von Mises准则，模型即简化为时间依赖的经典Saramito模型。进一步，若同时令剪切模量G→∞，模型则进一步退化为经典的Bingham（宾汉姆）模型。这证明了新模型涵盖并统一了从简单黏塑性到复杂弹黏塑性的行为描述。
结论与意义
本研究成功构建并验证了一个统一的、跨尺度的流固耦合计算框架，用于高保真模拟3D混凝土打印的全过程。其核心贡献在于：
- 1.
理论创新：首次为水泥基材料建立了基于Saramito框架的弹黏塑性本构模型，引入了压力敏感性和时间硬化，实现了从流体到固体状态的热力学一致平滑过渡。
- 2.
方法创新：将新本构模型与擅长处理大变形和自由表面的PFEM相结合，并发展了处理复杂边界条件（如挤出和沉积）的特定数值技术，在一个框架内无缝衔接了挤出尺度的流体模拟和结构尺度的固体力学分析。
- 3.
实践价值：通过系统的实验验证，证明了模型在预测打印几何形状、层间相互作用、复杂节点形成以及结构失效模式等方面具有卓越的准确性和可靠性。
这项研究标志着3DCP数值模拟从“分立”走向“统一”的重要一步。它不仅为理解和优化3D打印工艺提供了前所未有的洞察力，其提出的“Saramito-DP”本构框架和PFEM求解策略，也为模拟其他具有类似流固转变特性的复杂材料（如地质材料、生物组织、食品等）的制造过程提供了新的思路和工具。未来，结合高性能计算与数据驱动方法，该框架有望成为实现3DCP“首次即正确”打印和数字化孪生的关键技术基础。
