<bold>3D</bold>打印混凝土永久模板叠合柱的抗压性能数值模拟研究

张治成，叶志凯，孙晓燕，王海龙，高君峰; ZHANG Zhicheng; YE Zhikai; SUN Xiaoyan; WANG Hailong; GAO Junfeng

网刊加载中。。。

使用Chrome浏览器效果最佳，继续浏览，你可能不会看到最佳的展示效果，

确定继续浏览么?

复制成功，请在其他浏览器进行阅读

3D打印混凝土永久模板叠合柱的抗压性能数值模拟研究 PDF

- ORCID：
张治成 ^1,2
✉
- ORCID：
叶志凯 ¹
- ORCID：
孙晓燕 ^1,3
✉
- ORCID：
王海龙 ^1,2
✉
- ORCID：
高君峰 ¹

1. 浙江大学建筑工程学院,杭州 310058； 2. 山西浙大新材料与化工研究院,太原 030002； 3. 浙江大学平衡建筑研究中心,杭州 310028

中图分类号： TU375.3

最近更新：2023-11-28

DOI：10.11835/j.issn.2096-6717.2022.112

摘要

为深入研究3D打印混凝土永久模板叠合柱的抗压性能，基于3D打印混凝土永久模板叠合柱及同尺寸整体现浇对照柱试验建立构件数值模型，模拟分析其轴压荷载-位移响应及失效形态。针对界面粘结性能、现浇混凝土抗压强度、打印模板厚度、荷载偏心距等参数开展3D打印混凝土永久模板叠合柱的抗压性能计算分析，研究表明：叠合柱轴压极限承载力随着薄弱界面剪切强度、刚度及现浇混凝土抗压强度的增大而增大。由于打印材料的抗压强度高于现浇混凝土，叠合柱抗压极限承载力提升率与打印模板厚度呈近似线性关系，叠合圆柱的抗压极限承载力随着荷载偏心距的增大而降低，呈近似线性负相关。此外，偏心距对叠合圆柱极限承载力下降幅度的影响大于现浇圆柱。

关键词

3D打印混凝土; 永久模板; 叠合柱; 抗压性能; 数值模拟

混凝土结构是最主要的建筑结构形式，模板在混凝土结构空间造型实现和强度发展中起着关键作用^[

1]，模板工程占混凝土结构工程造价的20%~30%，占工程总用工量的30%~40%，占施工工期的50%左右^{[参考文献 2

百度学术}2]。现有钢、木模的使用耗费资源，造成严重环境污染和资源浪费。3D打印混凝土技术发展给永久模板提供了新的解决思路^{[参考文献 3

百度学术}3]，可根据工程需要进行数字建模、打印定制异形模板，实现免模施工，节约人力、物力，加快施工进程。

学者们针对3D打印永久模板叠合结构开展了试验研究和工程实践。Zhu等^[

4]设计3D打印永久模板的混合料配方，研究了配筋率对3D打印永久模板混凝土叠合柱轴压承载能力的影响，研究表明永久模板和现浇混凝土界面具有良好的粘结性能，叠合柱的承载能力优于相应现浇柱。Liu等^{[参考文献 5

百度学术}5]采用玄武岩纤维和碳纤维加强的混凝土打印永久模板，研究了打印层数和打印路径对拉弯性能的增强效果，通过EMI技术监测加载过程中的损伤演化，提出了一个解析模型来预测叠合梁抗弯强度。Wang等^{[参考文献 6

百度学术}6]通过试验和微细观CT扫描，研究了3D打印混凝土永久模板和后浇混凝土的界面粘结性能，并建立了界面剪切强度的塑性极限解析力学模型，结果表明，当打印层高度（表面粗糙度评价指标）为20 mm、后浇时间间隔7 d时，永久模板和后浇混凝土界面具有最好的粘结性能。高君峰^{[参考文献 7

百度学术}7]开展了3D打印永久模板钢筋混凝土叠合梁柱构件的弯剪压力学性能试验，并设置现浇对照组进行了对比。Bai等^{[参考文献 8

百度学术}8]通过3D打印不规则拱桥主拱模板，后在模板内放置钢筋笼，浇筑混凝土，在河北工业大学内建造起一座拱桥。Anton等^{[参考文献 9

百度学术}9]设计了一条3D打印定制复杂形态混凝土柱模板的自动化生产线并提出了基于几何复杂性的评价方法。但当前学者主要从试验方面研究3D打印永久模板叠合结构的性能，基本未涉及数值分析。

打印工艺导致成型混凝土具有显著空间各向异性，无法采用传统各向同性材料计算模式进行分析。采用将各向同性连续体单元和不考虑物理厚度的离散界面单元相结合^[

10]，开展3D打印混凝土结构数值分析成为一种技术趋势。Xiao等^{[参考文献 11

百度学术}11]基于ABAQUS建立界面模型，研究3D打印混凝土试块在压、弯荷载作用下，条带尺寸、层条间界面粘结性能、混凝土强度对3D打印混凝土试块各向异性力学行为的影响。Van Den Heever等^{[参考文献 12

百度学术}12]基于DIANA建立界面模型，预测不同荷载配置下钢筋混凝土深梁的结构承载力和破坏机制。但针对3D打印混凝土永久模板叠合构件力学性能的数值模拟仍有待深入研究。数值模拟作为试验研究的辅助方法，可以快速开展不同荷载配置下构件力学响应、破坏机理及影响参数的研究，直观呈现出由“堆叠成型”打印工艺形成的薄弱界面对结构在荷载作用下损伤破坏的影响。

柱式构件为结构中的重要承压构件，其受压性能对结构整体性能至关重要。为了给3D打印叠合结构工程设计及应用提供依据，笔者基于3D打印永久性混凝土模板叠合柱抗压模型试验，建立构件数值分析模型，针对界面黏结性能、现浇混凝土强度及荷载偏心进行参数分析。

1 3D打印永久模板混凝土叠合柱受压模型试验

1.1 试验概况

基于轴压试验研究3D打印永久模板-钢筋混凝土叠合短柱的受压性能^[

7]。共设计8根试验柱，其中叠合方柱2根，现浇方柱2根，叠合圆柱2根，现浇圆柱2根，截面及构件尺寸如图1所示，保护层厚度均为D=40 mm。打印材料采用3D打印PVA纤维增强混凝土^{[参考文献 13

百度学术}13]，配合比如表1所示。胶凝材料为42.5快硬早强型硫铝酸盐水泥，添加12 mm长度的聚乙烯醇纤维（PVA）。

图1 3D打印混凝土叠合柱配筋图

Fig. 1 Reinforcement of 3D printed concrete composite column

表1 3D打印混凝土配合比（质量比）

Table 1 3D-printed concrete mix ratio

水泥	矿粉	硅灰	砂	水	纤维(体积比)	减水剂	缓凝剂
1	0.8	0.2	0.61	0.33	1.2%	0.05	0.057

3D打印永久模板-钢筋混凝土叠合柱制作工序如图2所示。首先使用HC-3DPRT/L 3D打印机打印制作永久模板。喷嘴直径为30 mm，模板厚度为条带宽度d=30 mm，条带高度为10 mm。

（a）打印永久柱模板

（b）柱模板打印完成

（c）置入钢筋笼

（d）浇筑中心混凝土

图2 永久性混凝土模板叠合柱制作过程^[

Fig. 2 Permanent concrete formwork composite column manufacturing process

按设计图绑扎钢筋笼，在钢筋上粘贴应变片，放入3D打印纤维增强混凝土永久模板内，浇筑前用木条控制钢筋笼与打印模板间距为10 mm，共计保护层厚度D=40 mm。钢筋笼尺寸如图1所示，边长/直径为120 mm，箍筋间距为60 mm，纵筋为HRB400，箍筋为HPB300，强度采用拉伸试验实测值，最后浇筑C35混凝土振捣，常温下室内养护28 d。

1.2 加载方案

加载布置及应变片布置如图3所示。在柱子的四周安装4个位移传感器，用来测量柱子的轴向位移。在柱子中间位置的前后面各安装4个竖向应变片，在柱子的左右面安装4个横向应变片，共计16个应变片，用来测量混凝土在轴压荷载作用下的纵向及横向应变。采用1 000 t微机控制电液伺服多功能试验机，依据《混凝土结构试验方法标准》（GB 50152—2012）^[

14]的相关规定进行加载，试件在正式加载前进行预加载，然后进行卸载，正式加载采用分级加载，每级加载 20 kN，加载后停留3~5 min，待数据稳定后采集数据。

（a）叠合方柱模型图

（b）方柱试验图^[

（c）叠合圆柱模型图

（d）叠合圆柱试验图^[

图3 现场加载及测点布置图

Fig. 3 Loading and measuring point layout

2 数值模拟

2.1 基于界面模型的数值分析方法

用基于混凝土塑性损伤模型的连续体单元（CDP）模拟打印基体，采用基于内聚力模型（CZM）的离散界面单元模拟3D打印混凝土薄弱界面，进行打印混凝土结构计算分析。

2.1.1 基体模拟

采用图4所示的混凝土塑性损伤模型模拟3D打印混凝土基体的拉压非线性行为^[

11]。

图4 混凝土塑性损伤模型本构曲线

Fig. 4 Damage plastic constitutive curve of concrete

采用Saenz^[

15]提出的式（1）描述基体受压行为。

σ_{c} = \frac{E_{0} ε_{c}}{1 + [(E_{0} ε_{p} / σ_{p}) - 2] (ε_{c} / ε_{p}) + (ε_{c} / ε_{p})^{2}}

（1）

式中： $σ_{c}$ 和 $ε_{c}$ 分别为混凝土受压时的应力和应变； $σ_{p}$ 和 $ε_{p}$ 分别取试验确定的抗压强度和0.002； $E_{0}$ 为试验确定的混凝土初始切线模量。

采用Hordijk^[

16]提出的拉伸软化曲线（式（2））描述基体受拉行为。

\frac{σ_{t}}{f_{t}} = [1 + {(c_{1} \frac{w_{t}}{w_{c r}})}^{3}] e^{(- c_{2} \frac{w_{t}}{w_{c r}})} - \frac{w_{t}}{w_{c r}} (1 + c_{1}^{3}) e^{- c_{2}}

（2）

式中： $σ_{t}$ 为垂直于裂纹方向的拉应力； $f_{t}$ 为混凝土单轴抗拉强度； $c_{1}$ 和 $c_{2}$ 是系数，分别取3.0和6.93； $w_{t}$ 为裂缝张开位移， $w_{c r}$ 为应力或者断裂能完全释放时的裂缝张开位移，mm，可以采用断裂能 $G_{F}$ 来计算。

w_{c r} = 5.14 \frac{G_{F}}{f_{t}}

（3）

式中： $f_{t}$ 和 $G_{F}$ 应该由相应的试验进行标定，但目前尚未建立起3D打印纤维增强混凝土断裂能的计算公式，在没有试验数据时，断裂能 $G_{F}$ 可以近似采用CEB-FIP^[

17]给出的经验公式进行计算。

G_{F} = (0.046 9 {d_{a}}^{2} - 0.5 d_{a} + 26) (f_{c} {/ 10)}^{0.7}

（4）

拉压损伤因子 $d_{t} 、 d_{c}$ 用来描述混凝土拉压损伤，受压损伤采用Yu等^[

18]给出的公式。

d_{c} = 1 - \frac{σ_{c}}{f_{c}}

（5）

受拉损伤采用Chen等^[

19]给出的公式。

γ = {(1 - \frac{ε_{c r}}{ε_{c r, u}})}^{n}

（6）

d_{t} = 1 - γ

（7）

式中： $ε_{c r, u}$ 为混凝土达到 $w_{c r}$ 完全破坏时的开裂应变；指数n控制剪切退化速率，n取值为2~5。

2.1.2 打印层条薄弱界面模拟

采用牵引-分离定律（Traction-Separation Law）来描述内聚力模型中界面的弹性行为、损伤的产生演化以及最终破坏失效^[

20]。本构曲线如图5所示，界面弹性行为用式（8）表示。

图5 内聚力单元线性软化牵引-分离曲线

Fig. 5 Traction-separation curve with linear softening

\{\begin{matrix} t_{n} \\ t_{s} \\ t_{t} \end{matrix}\} = [\begin{matrix} k_{n 0} & 0 & 0 \\ 0 & k_{s 0} & 0 \\ 0 & 0 & k_{t 0} \end{matrix}] \{\begin{matrix} δ_{n} \\ δ_{s} \\ δ_{t} \end{matrix}\}

（8）

式中： $t_{n}$ 、 $t_{s}$ 、 $t_{t}$ 分别为界面法向和两个剪切方向的名义应力； $δ_{n}$ 、 $δ_{s}$ 、 $δ_{t}$ 分别为界面法向和两个剪切方向的相对位移； $k_{n 0}$ 、 $k_{s 0}$ 、 $k_{t 0}$ 分别为法向和两个剪切方向的弹性刚度。

损伤起始采用二次名义应力准则^[

20]，当二次函数值达到1时，则认为损伤开始产生。

{\{\frac{<t_{n}>}{t_{n 0}}\}}^{2} + {\{\frac{t_{s}}{t_{s 0}}\}}^{2} + {\{\frac{t_{t}}{t_{t 0}}\}}^{2} = 1

（9）

式中： $t_{n 0}$ 和 $t_{s 0}$ (= $t_{t 0}$ )是材料的抗拉强度和抗剪强度； $<>$ 是Macaulay括号，用来表征压应力不会造成损伤^[

21]。

线性软化规则下，损伤因子D是拉剪混合模式下有效相对位移 $δ_{m}$ 的函数^[

20]，取值伴随损伤演化从0到1，用式（10）计算。

D = \frac{δ_{m f} (δ_{m, m a x} - δ_{m 0})}{δ_{m, m a x} (δ_{m f} - δ_{m 0})}

（10）

δ_{m} = \sqrt[]{{<δ_{n}>}^{2} + {δ_{s}}^{2} + {δ_{t}}^{2}}

（11）

式中： $δ_{m, m a x}$ 为加载历程中混合模式下最大有效相对位移； $δ_{m 0}$ 和 $δ_{m f}$ 为混合模式下初始损伤和完全破坏时的有效相对位移。

QUADSCRT为采用二次名义应力准则的损伤起始判定指标：当数值小于1时，表示未损伤；达到1时，表示损伤开始。

MMIXDME为损伤演化过程中损伤类型判定指标：当数值为-1时，表示Cohesive单元没有损伤；当数值为0~0.5时，表示Cohesive单元以裂缝张开损伤为主；当数值为0.5~1.0时，表示Cohesive单元以剪切滑移损伤为主。

2.2 模型概况

基于ABAQUS，建立现浇/叠合柱限元模型，混凝土采用8节点线性缩减积分沙漏控制实体单元（C3D8R），钢筋采用2节点空间线性梁单元（B31），内置（Embedded）于混凝土中，叠合柱在层间界面处插入8节点三维Cohesive单元（COH3D8）。叠合方柱共计37 108个单元，叠合圆柱共计35 200个单元，如图6所示。

（a）叠合方柱

（b）方柱网格划分

（c）叠合圆柱

（d）圆柱网格划分

图6 轴压叠合柱有限元模型图

Fig. 6 Finite element model of axial composite column

边界条件和试验保持一致，钢垫板和柱子设置面面接触，法线方向为“硬接触”，切线方向用罚函数定义摩擦，摩擦系数取0.2。施加50 mm位移荷载，直至钢筋混凝土柱子破坏。3D打印永久模板在堆叠成型过程中自然形成凹凸不平的波纹形表面，增强了永久模板和现浇混凝土之间的黏结性能，使得两者能够形成整体协同工作^[

4]。因而，简化地将打印永久模板和现浇混凝土的接触面设置为绑定（Tie）约束。

2.2.1 参数取值

现浇混凝土和打印混凝土的塑性损伤模型拉压塑性参数采用式（1）~式（7）计算。根据《混凝土结构设计规范》（GB 50010—2010），现浇混凝土标准立方体抗压强度实测值为35.1 MPa^[

7]，换算为棱柱体抗压强度为23.5 MPa；打印混凝土抗拉强度值取3.0 MPa，打印混凝土采用边长为70.7 mm的立方体，抗压强度实测值为69.81 MPa^{[参考文献 22

百度学术}22]，换算为标准立方体抗压强度后，再换算为棱柱体抗压强度为41.1 MPa。

在ABAQUS/Standard 隐式分析中，黏性系数对收敛速度、开裂形态、损伤分布和宏观力学响应有影响^[

23]，黏性系数取值过大会造成刚化，影响破坏模式及结果准确性，经试算取1×10^-4，在合理的取值范围内^{[参考文献 24

百度学术}24]。泊松比v、第2应力不变量之比

K_{c}

、双轴与单轴抗压强度比

f_{b 0} / f_{c 0}

的取值参考文献[25]，剪胀角

α

的取值参考文献[19]。混凝土CDP模型参数汇总于表2。钢筋强度采用拉伸试验值，采用双线性随动强化模型计算塑性参数^{[参考文献 26

百度学术}26]，汇总于表3，强化段的切线模量取

E_{t} = 0.01 E_{s}

。

表2 混凝土塑性损伤模型参数

Table 2 Plastic damage model parameters of concrete

$ρ$ /(10^-9t/mm³)	$E_{3 D}$ /MPa	v	$α /$ (°)	$K_{c}$	e	$f_{b 0} / f_{c 0}$	粘性系数
2.5	27 850	0.2	35	0.666 7	0.1	1.16	0.000 1

表3 钢筋参数^{[参考文献 7

百度学术}7]

Table 3 Parameters of reinforcement

钢筋型号	$ρ$ /(10^-9t/mm³)	$E_{s}$ /10⁵ MPa	v	$f_{y}$	$f_{u}$	$ε_{u}$
HPB335-6	7.85	2.1	0.3	396	512	0.057 1
HRB400-12	7.85	2.1	0.3	503	598	0.047 6
HRB400-14	7.85	2.1	0.3	552	660	0.054 1

界面Cohesive单元拉伸刚度参考文献[

27]给出的数值取为1×10⁶ MPa。界面拉伸强度、抗剪强度及刚度采用同配比打印混凝土的拉伸试验^{[参考文献 28

百度学术}28]及双剪试验^{[参考文献 29

百度学术}29]测试数值。界面Cohesive单元拉伸断裂能和剪切断裂能参考Van Den Heever^{[参考文献 27

百度学术}27,30]给出的数据。参数取值汇总于表4。

表4 基于T-S Law的Cohesive单元材料特性

Table 4 Materials properties of cohesive unit based onT-S Law

位置	$\begin{array}{l} σ_{m a x} / \\ M P a \end{array}$	$\begin{array}{l} τ_{m a x} / \\ M P a \end{array}$	$\begin{array}{l} E_{n} / \\ M P a \end{array}$	$\begin{array}{l} E_{t t / s s} / \\ M P a \end{array}$	$G_{f n n}$ /(N/mm)	$G_{f t t / s s} /$ (N/mm)
界面	1.636	5.18	1E6	400	0.060	0.15

2.2.2 网格敏感性分析

采用15、10、7 mm网格尺寸对现浇方柱有限元模型进行网格收敛性分析，结果如图7所示。由于10、7 mm的网格尺寸计算精度相近，且具有较高的计算效率，因此，采用该网格尺寸建立基准模型。

图7 轴压柱网格收敛性分析

Fig. 7 Mesh convergence analysis of axial compression column

2.3 方柱轴压试验与模拟结果对比分析

2.3.1 试验现象

3D打印叠合方柱破坏现象与现浇方柱基本一致，试件在轴压应力作用下向外膨胀，产生第一条裂缝，随着轴压应力继续增大，裂缝发展为贯穿性斜裂缝，局部混凝土被压碎，随后压应力减小混凝土变形快速增大。

2.3.2 荷载-位移响应

3D打印叠合方柱荷载-位移响应与整体现浇柱基本一致，均呈现出短柱受压全过程的特征，试验过程中试验机刚度引起位移测量误差，造成加载初期柱刚度偏小，对试验初期刚度修正，数值模拟能较好地反映柱受压过程。

提取数值模型中刚开始产生拉伸损伤时对应的轴压应力作为开裂应力与试验进行对比，数值模拟的开裂应力均小于试验值，是因为损伤产生初期，微裂缝不能被明显观察到。叠合方柱开裂应力模拟结果与实验误差为-22.5%，现浇方柱开裂应力模拟结果与试验误差为-30.4%。叠合方柱极限承载力模拟结果与试验误差为3.3%，现浇方柱极限承载力模拟误差为6.7%，如表5、图8所示。

表5 方柱抗压极限承载力试验和数值模拟对比

Table 5 Comparison between FEA and test results of ultimate bearing capacity of square column

试件编号	开裂应力/MPa	模拟误差/%	开裂应力提高/%	极限应力/MPa	模拟误差/%	极限应力提高/%
Cast-Test	7.9			33.3
Cast-FEA	5.5	-30.4		35.5	6.6
3DCP-Test	20.4		158.2	41.0		23.1
3DCP-FEA	15.8	-22.5	187.2	42.3	3.2	19.2

图8 方柱试验与数值模拟轴压应力-位移曲线

Fig. 8 Comparison between FEA and test results of square column load-displacement curves

叠合方柱的平均开裂应力试验值比现浇方柱提高158.2%，叠合方柱的开裂应力模拟值比现浇方柱提高187.2%。这是由于3D打印永久模板混凝土添加了PVA纤维，使得拉伸强度和断裂能均大于现浇混凝土，且在打印挤出工艺下，PVA纤维沿环向定向分布，一定程度上提高了叠合柱的抗裂性能。

叠合方柱平均极限承载力试验值高于现浇柱23.1%，叠合方柱极限承载力模拟值高于现浇柱19.2%。这是因为纤维增强3D打印混凝土抗压强度高，使得叠合柱极限承载力有所提高。

2.3.3 破坏模式

取轴压应力-位移曲线进入下降段，轴向位移2 mm，产生压损伤时的计算结果，如图9、图10所示。现浇方柱和叠合方柱均呈现出剪切破坏模式，钢筋屈服进入强化段。对于叠合方柱，模板压损伤大于中心混凝土，模板在中间剪斜裂缝的交汇处及端部处损伤最严重。端部在摩擦力和轴压力共同作用下，Cohesive单元剪切损伤快速发展，造成端部模板局部压碎剥落。

（a）现浇方柱破坏模式^[

（b）现浇方柱受压损伤

（c）现浇方柱钢筋笼应力

图9 现浇方柱试验与数值模拟对比

Fig. 9 Comparison between FEA and test results of cast square column

（a）叠合方柱破坏模式^[

（b）核心混凝土受压损伤

（c）叠合方柱钢筋笼应力

（d）打印模板受压损伤

（e） Cohesive单元损伤判定（f ） Cohesive单元损伤类型

图10 3D打印叠合方柱试验与数值模拟对比

Fig. 10 Comparison between FEA and test results of 3D printed square column

2.4 圆柱轴压试验与模拟结果对比分析

2.4.1 试验现象

3D打印叠合圆柱破坏现象与现浇圆柱基本一致，试件在轴压应力作用下向外膨胀，产生第一条裂缝，随着轴压应力继续增大，裂缝发展为贯穿性裂缝，局部混凝土被压碎，随后压应力减小混凝土变形快速增大。

2.4.2 荷载-位移响应

3D打印永久模板叠合圆柱的荷载-位移响应与现浇圆柱基本一致，均呈现出短柱受压的特征，试验过程中试验机刚度引起位移测量误差，造成加载初期柱刚度偏小，对试验初期刚度修正，数值模拟能较好地反映柱受压过程。

提取数值模型中刚开始产生拉伸损伤时对应的轴压应力作为开裂应力与试验进行对比，叠合圆柱开裂应力模拟结果与实验误差为-23.8%，现浇圆柱开裂应力模拟结果与试验误差为-28.8%，数值模拟的开裂应力均小于试验值，是因为损伤产生初期，微裂缝不能被明显观察到。叠合圆柱极限承载力模拟结果与试验误差为8.3%，现浇圆柱极限承载力模拟误差为5.5%，如表6、图11所示。

表6 圆柱抗压极限承载力试验和数值模拟对比

Table 6 Comparison between FEA and test results of ultimate bearing capacity of circle column

试件编号	开裂应力/MPa	模拟误差/%	开裂应力提高/%	极限应力/MPa	模拟误差/%	极限应力提高/%
Cast-Test	8.8			32.8
Cast-FEA	6.7	-23.8		34.6	5.5
3DCP-Test	21.5		144.3	44.6		36.0
3DCP-FEA	15.3	-28.8	128.4	48.3	8.3	39.6

图11 圆柱试验与数值模拟轴压应力-位移曲线

Fig. 11 Comparison between FEA and test results of circle column load-displacement curves

叠合圆柱的平均开裂应力试验值比现浇方柱提高144.3%，叠合圆柱的开裂应力模拟值比现浇方柱提高128.4%。这是由于3D打印永久模板混凝土添加了PVA纤维，使得拉伸强度和断裂能均大于现浇混凝土，且在打印挤出工艺下，PVA纤维沿环向定向分布，一定程度上提高了叠合柱的抗裂性能。

叠合圆柱平均极限承载力试验值高于现浇柱36.0%，叠合圆柱极限承载力模拟值高于现浇柱39.6%，这是因为纤维增强3D打印混凝土抗压强度高，使得叠合柱极限承载力有所提高。

2.4.3 破坏模式

取轴压应力-位移曲线进入下降段，产生压损伤，位移2 mm时候的计算结果，如图12、图13所示，现浇圆柱和叠合圆柱均呈现出剪切破坏模式，钢筋屈服进入强化段。

（a）现浇圆柱破坏模式^[

（b）现浇圆柱受压损伤

（c）现浇圆柱钢筋笼应力

图12 现浇圆柱试验与数值模拟对比

Fig. 12 Comparison between FEA and test results of cast circle column

（a）叠合圆柱破坏模式^[

（b）核心混凝土受压损伤

（c）叠合圆柱钢筋笼应力

（d）打印模板受压损伤

（e） Cohesive单元损伤判定

（f） Cohesive单元损伤类型

（g）打印模板受压损伤（τ=3）

（h） Cohesive单元损伤判定（τ=3）

（i） Cohesive单元损伤类型（τ=3）

图13 3D打印叠合圆柱试验与数值模拟对比

Fig. 13 Comparison between FEA and test results of 3D printed circle column

现浇圆柱破坏模式为端部的混凝土被压碎，端部附近位置处产生斜裂缝，试验中斜裂缝产生在柱子顶端，而数值模拟中斜裂缝在柱子底端开展，如图12（a）所示。3D打印永久模板叠合圆柱破坏模式是的剪切破坏，斜裂缝交汇于三分之一高度处，造成该位置处的混凝土被剪碎，进而有可能会形成如图13（a）所示的贯穿性裂缝。当存在薄弱界面时，端部处在轴压力、钢垫板摩擦力及泊松效应作用下，层条会产生横向剪切滑移，界面产生剪切损伤，剪切破坏朝端部移动，端部产生斜裂缝以及混凝土被压碎，Cohesive单元剪切强度越低，这种现象越明显，图13（g）~（i）所示为Cohesive单元剪切强度 $τ$ =3 的计算结果。

3 3D打印叠合柱抗压性能参数分析

相对于现浇柱而言，在轴压荷载作用下，叠合柱极限承载能力有30%左右的提升。考虑到实际过程中圆柱支模难度更大，在充分考虑3D打印永久模板叠合柱的技术优势以及截面的经济优势情况下，3D打印永久模板-钢筋混凝土叠合圆柱更具有工程应用潜力，因而需要充分了解在压荷载作用下叠合圆柱极限承载力的主要影响因素。选取4个参数进行分析。

3.1 界面黏结性能

叠合柱仅存在垂直于荷载方向的薄弱界面，影响柱子受压性能的界面黏结参数主要是Cohesive单元的剪切强度和剪切模量。目前3D打印纤维增强混凝土界面剪切性能的试验研究较少，剪切强度分别选取试验值的0.2、0.4、2.0倍，即选取 $τ =$ 1、2、5、10 MPa以及不带Cohesive单元的情况进行计算；剪切模量分别取 $E_{t t} =$ 4、4×10¹、4×10²、4×10³、4×10⁴ MPa进行计算，计算结果如图14、图15所示。

（a） Cohesive不同剪切强度下的轴压应力-位移曲线

（b） Cohesive不同剪切强度下峰值对比

图14 Cohesive剪切强度对叠合圆柱承载力的影响

Fig. 14 Effect of interfacial bond strength on the 3D printed composite column capacity

（a） Cohesive不同剪切模量下的轴压应力-位移曲线

（b） Cohesive不同剪切模量下峰值对比

图15 Cohesive剪切模量对叠合圆柱承载力的影响

Fig. 15 Effect of interfacial shear module on the 3D printed composite column capacity

轴压荷载下叠合柱抗压极限承载力随着界面抗剪强度的增大而增大，呈现出幂函数形式。当界面剪切强度超过5 MPa之后，界面的存在对叠合柱抗压极限承载力的削弱趋于平缓；当剪切强度为1 MPa时，叠合柱抗压极限承载力相比于没有薄弱界面的情况下降低16.2%。轴压荷载下叠合柱抗压极限承载力随着界面剪切模量的增大而增大，呈现出对数函数形式。当界面剪切模量超过400 MPa之后，界面的存在对叠合柱抗压极限承载力的削弱趋于平缓；当剪切模量为4 MPa时，叠合柱抗压极限承载力相较于剪切模量为400 MPa的情况下降低19.8%。

3.2 现浇混凝土强度

PVA纤维增强3D打印混凝土材料抗压强度高，中间浇筑普通混凝土进行协同工作，因而考虑中间现浇混凝土抗压强度对叠合柱受轴压荷载作用下的影响。中间现浇混凝土标准立方体抗压强度分别取 $f_{c c} =$ 30、35、40、45 MPa进行计算，计算结果如图16所示。

（a）不同抗压强度下的轴压应力-位移曲线

（b）不同抗压强度下峰值对比

图16 现浇混凝土抗压强度对叠合圆柱承载力的影响

Fig. 16 Effect of concrete compression strength on the 3D printed composite column capacity

叠合柱的抗压极限承载力随着中心现浇混凝土强度的提高而提高，呈现出线性正相关形式，当中心现浇混凝土强度从35 MPa提升到45 MPa时，叠合柱抗压极限承载力提高10.3%，因而可以通过适当提高现浇混凝土强度来提高叠合柱的抗压极限承载力。

3.3 打印模板厚度

受打印设备的喷嘴尺寸限制，试验柱的保护层占截面面积比例过大。为了探究打印模板厚度对叠合柱承载能力的影响，分别选取打印模板厚度d为10、15、20、25、30 mm，保护层厚度D分别为20、25、30、35、40 mm开展多工况分析，且将打印模板厚度d=10 mm的叠合柱与保护层厚度 $D = 20$ mm的现浇柱对比。

由于打印材料的抗压强度高于现浇混凝土，随着打印模板占截面面积比的增大，叠合柱抗压极限承载力随之提高,如图17(a)所示。承载力提升率与打印模板厚度呈线性关系，如图17（b）所示。

（a）不同打印模板厚度下的轴压应力-位移曲线

（b）不同打印模板厚度下峰值对比

图17 打印模板厚度对叠合圆柱承载力的影响

Fig. 17 Effect of printing formwork thickness on the 3D printed composite column capacity

3.4 荷载偏心距

在实际工程中，柱子受力通常并非理想轴压状态，当压荷载存在偏心时，远离荷载侧在偏压荷载作用下将会产生拉应力，而3D打印工艺下存在的薄弱界面的抗拉强度较低，其对结构整体性能的不利影响将更为显著。分别取偏心距 $e =$ 0、25、50、75 mm对叠合圆柱和现浇圆柱进行计算。当位移达到2 mm时，曲线进入下降段，取位移为2 mm左右的计算结果，其中 $e = 75 m m$ 难以收敛，取最后增量步的计算结果，叠合圆柱在不同偏心距下的轴压应力-位移曲线如图18（a）所示；叠合圆柱及现浇圆柱在不同偏心距下的极限承载力如图18（b）所示；不同偏心距对叠合圆柱的损伤影响如图19所示。

（a）不同荷载偏心距下的轴压应力-位移曲线

（b）不同荷载偏心距下的峰值对比

图18 不同荷载偏心距对叠合圆柱承载力的影响

Fig. 18 Effect of different eccentricity on the 3D printed composite column capacity

（a） Cohesive损伤判定（e=25）

（b） Cohesive损伤判定（e=50）

（c） Cohesive损伤判定（e=75）

（d） Cohesive损伤类型（e=25）

（e） Cohesive损伤类型（e=50）

（f） Cohesive损伤类型（e=75）

图19 不同偏心距对叠合圆柱损伤的影响

Fig. 19 Effect of different eccentricity on the interface failure of composite column

叠合圆柱的抗压极限承载力随着荷载偏心距的增大而降低，呈现线性负相关形式。同等荷载下，荷载偏心距越大，受拉侧界面的拉伸损伤越严重，叠合圆柱抗压极限承载力降幅越大。当偏心距达到二分之一叠合柱半径，即 $e = 50 m m$ 时，叠合圆柱抗压极限承载力降低67.9%。此外，偏心距对叠合圆柱极限承载力下降幅度的影响大于现浇圆柱。

4 结论

1）3D打印永久模板-钢筋混凝土叠合柱及同尺寸整体现浇柱的抗压性能对照试验表明，叠合柱整体受力性能及极限承载能力优于现浇柱。

2）从轴压应力-位移曲线和破坏形态来看，基于界面的数值模型可以较好地模拟叠合柱受压作用下的力学响应。

3）参数研究表明叠合柱轴压极限承载力随着薄弱界面剪切强度、刚度以及现浇混凝土抗压强度的增大而增大，但荷载的偏心会造成一侧的模板受拉，薄弱界面处产生的拉伸损伤会较大幅度降低叠合柱的抗压极限承载力。偏心距对叠合柱的极限承载力削弱幅度大于现浇柱。

4）在3D打印叠合柱的实际工程应用中，应该优化打印工艺，控制打印速度，进而改善界面黏性性能；适当提高中间浇筑混凝土强度；控制偏心荷载的出现。

参考文献

LI W, LIN X S, BAO D W, et al. A review of formwork systems for modern concrete construction [J]. Structures, 2022, 38: 52-63. [百度学术]

LLORET E, SHAHAB A R, LINUS M, et al. Complex concrete structures: Merging existing casting techniques with digital fabrication [J]. Computer-Aided Design, 2015, 60: 40-49. [百度学术]

MENNA C, MATA-FALCÓN J, BOS F P, et al. Opportunities and challenges for structural engineering of digitally fabricated concrete [J]. Cement and Concrete Research, 2020, 133: 106079. [百度学术]

ZHU B R, NEMATOLLAHI B, PAN J L, et al. 3D concrete printing of permanent formwork for concrete column construction [J]. Cement and Concrete Composites, 2021, 121: 104039. [百度学术]

LIU M, HUANG Y M, WANG F, et al. Tensile and flexural properties of 3D-printed jackets-reinforced mortar [J]. Construction and Building Materials, 2021, 296: 123639. [百度学术]

WANG L, YANG Y, YAO L, et al. Interfacial bonding properties of 3D printed permanent formwork with the post-casted concrete [J]. Cement and Concrete Composites, 2022, 128: 104457. [百度学术]

高君峰. 3D打印永久性混凝土模板及其叠合构件试验研究[D]. 杭州: 浙江大学, 2020. [百度学术]

GAO J F. Mechanical properties of 3D printed concrete permanent formwork and integrated reinforced concrete components [D]. Hangzhou: Zhejiang University, 2020. (in Chinese) [百度学术]

BAI G, WANG L, MA G W, et al. 3D printing eco-friendly concrete containing under-utilised and waste solids as aggregates [J]. Cement and Concrete Composites, 2021, 120: 104037. [百度学术]

ANTON A, BEDARF P, YOO A, et al. Concrete choreography: prefabrication of 3D printed columns [M]//Fabricate 2020. UCL Press, 2020: 286-293. [百度学术]

LOURENÇO P B, ROTS J G. Multisurface interface model for analysis of masonry structures [J]. Journal of Engineering Mechanics, 1997, 123(7): 660-668. [百度学术]

XIAO J Z, LIU H R, DING T. Finite element analysis on the anisotropic behavior of 3D printed concrete under compression and flexure [J]. Additive Manufacturing, 2021, 39: 101712. [百度学术]

VAN DEN HEEVER M, BESTER F, KRUGER J, et al. Numerical modelling strategies for reinforced 3D concrete printed elements [J]. Additive Manufacturing, 2022, 50: 102569. [百度学术]

SUN X Y, ZHOU J W, WANG Q, et al. PVA fibre reinforced high-strength cementitious composite for 3D printing: Mechanical properties and durability [J]. Additive Manufacturing, 2022, 49: 102500. [百度学术]

混凝土结构试验方法标准: GB/T 50152—2012 [S]. 北京: 中国建筑工业出版社, 2012. [百度学术]

Standard for test method of concrete structures: GB/T 50152—2012 [S]. Beijing: China Architecture & Building Press, 2012. (in Chinese) [百度学术]

SANEZ L P. Discussion of equation for the stress-strain curve of concrete’ by Desayi and Krishnan [J]. Journal of the American Concrete Institute, 1964, 61: . [百度学术]

HORDIJK D A. Local approach to fatigue of concrete [D]. Netherlands: Delft University of Technology, 1991. [百度学术]

CEB-FIP F. Fib model code for concrete structures 2010 [M]. Ernst And Sohn, 2013. [百度学术]

YU T, TENG J G, WONG Y L, et al. Finite element modeling of confined concrete-II: Plastic-damage model [J]. Engineering Structures, 2010, 32(3): 680-691. [百度学术]

CHEN G M, CHEN J F, TENG J G. On the finite element modelling of RC beams shear-strengthened with FRP [J]. Construction and Building Materials, 2012, 32: 13-26. [百度学术]

CAMANHO P P, DAVILA C G, DE MOURA M F. Numerical simulation of mixed-mode progressive delamination in composite materials [J]. Journal of Composite Materials, 2003, 37(16): 1415-1438. [百度学术]

TENG J G, FERNANDO D, YU T. Finite element modelling of debonding failures in steel beams flexurally strengthened with CFRP laminates [J]. Engineering Structures, 2015, 86: 213-224. [百度学术]

张静. 3D打印混凝土本构模型与梁的受弯性能研究[D]. 杭州: 浙江大学, 2021. [百度学术]

ZHANG J. Research on 3D printed concrete constitutive model and flexural performance of beam [D]. Hangzhou: Zhejiang University, 2021. (in Chinese) [百度学术]

SZCZECINA M, WINNICKI A. Relaxation time in CDP model used for analyses of RC structures [J]. Procedia Engineering, 2017, 193: 369-376. [百度学术]

CHEN H B, XU B, MO Y L, et al. Behavior of meso-scale heterogeneous concrete under uniaxial tensile and compressive loadings [J]. Construction and Building Materials, 2018, 178: 418-431. [百度学术]

WU J Y, LI J, FARIA R. An energy release rate-based plastic-damage model for concrete [J]. International Journal of Solids and Structures, 2006, 43(3/4): 583-612. [百度学术]

聂祺. 高层钢筋混凝土结构非线性动力时程分析研究[D]. 北京: 中国建筑科学研究院, 2009. [百度学术]

NIE Q. Research on nonlinear dynamic history analysis of reinforced concrete structures [D]. Beijing: China Academy of Building Research, 2009. (in Chinese) [百度学术]

VAN DEN HEEVER M, BESTER F, POURBEHI M, et al. Characterizing the fissility of 3D concrete printed elements via the cohesive zone method[M]// Second RILEM International Conference on Concrete and Digital Fabrication. Cham: Springer, 2020: 489-499. [百度学术]

张静, 邹道勤, 王海龙, 等. 3D打印混凝土层条间界面抗拉性能与本构模型[J]. 浙江大学学报(工学版), 2021, 55(11): 2178-2185, 2214. [百度学术]

ZHANG J, ZOU D Q, WANG H L, et al. Bond tensile performance and constitutive models of interfaces between vertical and horizontal filaments of 3D printed concrete [J]. Journal of Zhejiang University (Engineering Science), 2021, 55(11): 2178-2185, 2214. (in Chinese) [百度学术]

WANG H L, SHAO J W, ZHANG J, et al. Bond shear performances and constitutive model of interfaces between vertical and horizontal filaments of 3D printed concrete [J]. Construction and Building Materials, 2022, 316: 125819. [百度学术]

VAN DEN HEEVER M, BESTER F, KRUGER J, et al. Mechanical characterisation for numerical simulation of extrusion-based 3D concrete printing [J]. Journal of Building Engineering, 2021, 44: 102944. [百度学术]