Flexural Performance of PP-ECC Beams Under Coupling of Freeze-Thaw Cycles and Bending Loads
-
摘要:
为探究冻融以及冻融-弯曲荷载耦合作用对聚丙烯纤维水泥基复合材料（PP-ECC）梁抗弯性能的影响，设置7种试验工况，采用三分点加载方式对PP-ECC梁抗弯性能进行研究，分析冻融及冻融-弯曲荷载耦合作用下PP-ECC梁的荷载-跨中挠度曲线、抗弯承载力及裂缝发展形态差异；基于计算假定和正常环境下PP-ECC梁抗弯承载力的计算模型，结合PP-ECC材料冻融劣化机理，推导出PP-ECC梁在冻融环境下抗弯承载力计算模型；在此基础上引入持荷损伤系数
γ ，建立冻融-弯曲荷载耦合作用下PP-ECC梁抗弯承载力计算模型. 研究结果表明：不同持荷比的PP-ECC梁极限抗弯承载力在冻融循环作用下出现不同程度降低，500次冻融循环后，持荷比为0、0.25、0.50的PP-ECC梁极限抗弯承载力分别降低了28.70%、27.09%、35.69%；受拉区PP-ECC材料开裂后不退出工作仍能协同受拉钢筋参与全截面受力；PP-ECC 梁在达到极限状态时，受拉区呈多条裂缝稳态发展模式，且随冻融损伤加剧，梁体最大裂缝宽度增大，裂缝数量减少；单一冻融以及冻融-弯曲荷载耦合作用下的PP-ECC梁仍满足平截面假定；基于平截面假定建立的冻融和冻融-弯曲荷载耦合作用下PP-ECC 梁的抗弯承载力计算模型吻合度分别在0.88～1.06和0.96～1.10.-
关键词:
- 冻融-弯曲荷载耦合作用 /
- 抗弯性能 /
- 聚丙烯纤维水泥基复合材料 /
- 抗弯承载力 /
- 多裂缝发展
Abstract:To investigate the effect of freeze-thaw cycles and coupling of freeze-thaw cycles and bending loads on the flexural performance of polypropylene fiber cement-based composite (PP-ECC) beams, seven experimental working conditions were set up. In addition, the three-point loading method was adopted to explore the flexural performance of PP-ECC beams. The differences in load-span deflection curves, flexural bearing capacity, and crack development patterns of PP-ECC beams under coupling of freeze-thaw cycles and bending loads were analyzed. Based on the calculation assumptions and the calculation model of the flexural bearing capacity of PP-ECC beams under normal environments, combined with the freeze-thaw deterioration mechanism of PP-ECC materials, the calculation model of flexural bearing capacity of PP-ECC beams under freeze-thaw environments was deduced. Furthermore, the load-bearing damage coefficient γ was introduced on this basis. The calculation model of the flexural bearing capacity of PP-ECC beams under the coupling of freeze-thaw cycles and bending loads was established. The research results show that the ultimate flexural bearing capacity of PP-ECC beams with different load-bearing ratios decreases to varying degrees under freeze-thaw cycles. After 500 freeze-thaw cycles, the ultimate flexural bearing capacity of PP-ECC beams with load-bearing ratios of 0, 0.25, and 0.5 decreases by 28.70%, 27.09%, and 35.69%, respectively. After cracking under tension, the PE-ECC material can still work with the rebar under tension and participate in the stress of the whole section. When the PP-ECC beam reaches the limit state, the zone under tension develops steadily with multiple cracks. With the aggravation of freeze-thaw damage, the width of the largest crack in the beam body increases, and the number of cracks decreases. The PP-ECC beam subjected to freeze-thaw cycles and coupling of freeze-thaw cycles and bending loads still satisfies the plane section assumption. The coincidence degree of the flexural bearing capacity calculation models of PP-ECC beams under freeze-thaw cycles and coupling of freeze-thaw cycles and bending loads established based on the plane section assumption is 0.88–1.06 and 0.96–1.10, respectively.
-
与普通混凝土箱梁相比，波形钢腹板组合箱梁（composite box girders with corrugated steel webs, CBGCSWs）具有自重轻、预应力导入效率高、完全避免腹板开裂等诸多优势[1-3]，因此在我国得以迅速推广. 近十几年来，国内外学者对CBGCSWs的抗弯[4-5]、抗剪[6-7]等方面的受力性能进行了广泛研究. 然而，聚焦其纯扭性能方面的研究却较为匮乏. 事实上，用轻薄的波形钢腹板取代传统的混凝土腹板后，CBGCSWs的抗扭刚度会被大幅削弱（大约仅为传统混凝土箱梁的30%~40%）[8]. 因此，在扭转效应突出的情形下（如曲线梁桥、偏心荷载等），该结构的抗扭问题值得重点关注和深入研究.
针对波形钢腹板PC组合箱梁（prestressed concrete composite box girders with corrugated steel webs, PCCBGCSWs）纯扭作用下的受力问题，诸多学者均提出过相应的分析模型. Mo等[9]将变角软化桁架理论（rotating angle softened truss model for torsion, RA-STMT）应用于PCCBGCSWs，提出了一种预测该结构纯扭作用下力学性能的理论分析模型. 基于RA-STMT，同时考虑结构的受力特点，聂建国等[10]建立了PCCBGCSWs在纯扭作用下的非线性分析模型. 丁勇等[11]以扭转理论及固角软化桁架模型为基础，建立了PCCBGCSWs 在纯扭作用下的抗扭承载力计算模型. Ko等[12]提出了预测PCCBGCSWs纯扭受力行为的改进分析模型，该模型考虑了混凝土抗拉钢化效应，并对混凝土翼缘板内的剪力流有效厚度进行了修正. 但作者仅将该模型用于混凝土开裂后阶段的分析，而针对开裂前阶段的预测则基于Bredt薄壁构件扭转理论[13]. 沈孔健等[14]将PCCBGCSWs的全过程扭矩-扭率曲线分成混凝土开裂前及开裂后两个阶段，并针对上述两个阶段分别采用扭转刚度修正和RA-STMT进行分析. 基于RA-STMT，同时将混凝土抗拉强度考虑在内，Shen等[15]提出了针对单箱多室PCCBGCSWs纯扭全过程的分析模型. Shen等[16]以软化薄膜元模型（softened membrane model for torsion, SMMT）[17]及RA-STMT为基础，建立了PCCBGCSWs纯扭全过程分析模型. 之后，Shen等[16]对该模型进行了改进，在模型中考虑了预应力效应对初始应力及应变的影响[18].
从以上文献可以看出：已有的PCCBGCSWs纯扭分析模型大都基于RA-STMT. 然而，由于RA-STMT忽略了混凝土的抗拉强度和抗拉刚化效应，也不能考虑开裂混凝土提供的抗剪强度，因此基于该理论提出的分析模型无法预测结构在混凝土开裂前的扭矩-扭率曲线，对混凝土开裂后纯扭力学行为的预测也不够精确. 另外，有部分模型（如Ko模型[12]和沈模型[14]）分别采用不同的计算理论对PCCBGCSWs开裂前、后两个阶段进行分析，这样的处理方式显然不具有理论一致性. 除此之外，已有理论模型的求解程序均包含多个迭代循环，因此导致这些模型的求解效率比较低.
本文基于SMMT提出了用于分析PCCBGCSWs纯扭全过程受力行为的改进软化薄膜元模型（improved softened membrane model for torsion，ISMMT）[19-20]. 该分析模型基于单一软化薄膜元理论，考虑了混凝土抗拉强度及开裂混凝土提供的抗剪强度，且对剪力流有效厚度进行了合理修正. 文献[19]将国内外8根PCCBGCSWs纯扭试件的扭矩-扭率曲线与ISMMT预测的理论曲线进行了比较，初步验证了ISMMT的适用性. 事实上，ISMMT除了能够预测扭矩-扭率曲线，还能对各构件在加载全过程的应变历程进行模拟. 因此，该理论模型的准确性尚需进一步的验证.
本文首先对ISMMT的平衡、变形协调、材料本构方程以及通用求解程序进行介绍；之后，针对ISMMT和同类型分析模型中求解程序过于复杂且耗时的问题，提出当波形钢腹板、预应力及普通钢筋均处于弹性阶段时的简化求解程序框图；最后，完成1根PCCBGCSWs试件的纯扭模型试验，并将模型试验结果与采用ISMMT计算的理论结果进行对比，来进一步验证该理论模型的适用性与准确性.
1. 软化薄膜元模型（SMMT）简介
通过考虑受扭构件中混凝土斜压杆的受压软化效应，Hsu等[13]对空间桁架模型进行了改进，提出了经典的软化桁架模型（softened truss model, STM）. 该模型能准确预测纯扭构件在极限状态下的扭矩与扭转角，但由于其未考虑混凝土的抗拉强度，无法对结构在开裂前的扭转性能进行预测，因而具有局限性[21].
为了克服上述不足，Jeng等[17]将用于预测RC受剪构件受力行为的软化薄膜元模型（softened membrane model, SMM）进行了改进，提出了针对RC受扭构件的分析模型（softened membrane model for torsion, SMMT）. 该模型的创新点体现在：通过修正混凝土本构关系来考虑斜压杆的应变梯度效应，并首次将混凝土的抗拉强度考虑在内. 通过与已有文献中的纯扭试验数据进行对比，表明SMMT能对RC实心纯扭构件的全过程受力行为作出精准预测.
2. 改进软化薄膜元模型（ISMMT）
ISMMT是以软化薄膜元理论为基础，并将PCCBGCSWs的结构和力学特点考虑在内，进而得出的一种预测该结构纯扭全过程受力行为的理论分析模型. 以下对ISMMT的平衡、变形协调及材料本构方程以及通用求解程序进行简要介绍，其详细推导过程参见文献[19-20].
2.1 平衡、变形协调及材料本构方程
2.1.1 平衡方程
对于承受外扭矩的PCCBGCSWs，其截面内会形成连续闭合的剪力流q来平衡外部扭矩，如图1所示. 图中：
$ {A_0} $ 为组合箱梁中剪力流中心线所包含的截面面积，${A_0}{\text{ = }}b\left( {h - {t_{\rm{d}}}} \right)$ ，b、h分别为组合箱梁梁宽和梁高；td为剪力流有效厚度；$ T $ 为箱梁承受的扭矩；${\tau _{{{lt}}}}$ 为l-t坐标系中的平均剪应力；$ s $ 为箍筋间距.
图 1 纯扭作用下波形钢腹板PC组合箱梁Figure 1. Prestressed concrete composite box girder with corrugated steel webs subjected to pure torsion根据组合箱梁截面平衡状态可得
T={2A0τwtw,τw<τwy,2A0τwytw+2×2As(τlttd−τwtw),τw=τwy, (1) 式中：
$ {A_{\rm{s}}} $ 为在混凝土翼缘板内单独闭合的剪力流$ {q_2} $ 所包含的截面面积；${\tau _{\rm{w}}}$ 、${\tau _{{\rm{wy}}}}$ 分别为钢腹板的剪应力和剪切屈服强度；tw为钢腹板的厚度.在混凝土翼缘板内剪力流区域取微元A，该微元受纯剪作用. 根据微元A面内平衡状态可得
12(σ2c+σ1c)+τ21c−σci+ρlfl+ρpfp−ρlifli−ρpifpi=0, (2) 12(σ2c+σ1c)−τ21c+ρtft=0, (3) τlt=12(−σ2c+σ1c), (4) 式中：
${\sigma _{1{\rm{c}}}}$ 、${\sigma _{2{\rm{c}}}}$ 分别为1与2方向的混凝土平均正应力；${\sigma _{21{\rm{c}}}}$ 为1-2坐标系中混凝土的平均剪应力；${f_{{l}}}$ 、${f_{{t}}}$ 、${f_{\rm{p}}}$ 分别为纵筋、箍筋及预应力筋的应力；${\rho _{{l}}}$ 、${\rho _{{t}}}$ 、${\rho _{\rm{p}}}$ 分别为纵筋、箍筋及预应力筋的含筋率，如式（5）；${\sigma _{{\rm{ci}}}}$ 、${f_{{l\rm{i}}}}$ 、${f_{{\rm{pi}}}}$ 分别为由预应力引起的混凝土（式（6））、纵筋（式（7））及预应力筋的初应力；${\rho _{{l\rm{i}}}}$ 、${\rho _{{\rm{pi}}}}$ 分别为扭转前纵筋及预应力筋的含筋率，如式（8）.{ρl=−Albtd,ρt=−Atstd,ρp=−Apbtd, (5) σci=−Ecˉεli=−EcApfpiEcAcc+EsAl, (6) fli=−Esˉεli=−EsApfpiEcAcc+EsAl, (7) {ρli=AlwAcc,ρpi=ApwAcc, (8) 式（5）~（8）中：
${A_{{l}}}$ 、${A_{\rm{p}}}$ 分别为剪力流区间内的纵筋和预应力筋面积（单块混凝土翼缘板内）；${A_{\rm{t}}}$ 为单根箍筋面积；${E_{\rm{c}}}$ 、${E_{\rm{s}}}$ 分别为混凝土和钢筋的弹性模量；${\bar \varepsilon _{l{\rm{i}}}}$ 为预应力引起的l方向初应变；${A_{{\rm{cc}}}} = {A_{\rm{c}}} - {A_{{l\rm{w}}}} - {A_{{\rm{pw}}}}$ ，${A_{\rm{c}}}$ 为混凝土翼缘板截面面积，${A_{{l\rm{w}}}}$ 、${A_{{\rm{pw}}}}$ 分别为全部纵筋和预应力筋的面积.2.1.2 变形协调方程
根据Bredt薄壁构件扭转理论[13]，可得组合箱梁截面变形协调方程为
∮γds=2(γltb+γwhw)=2θA0, (9) 式中：
${\gamma _{{{lt}}}}$ 为l-t坐标系中的平均剪应变；${\gamma _{\rm{w}}}$ 、${h_{\rm{w}}}$ 分别为钢腹板的剪应变和高度；$ \theta $ 为组合箱梁扭率.基于莫尔应变圆相关理论，可推导出微元A的3个变形协调方程为
{εl=12(ε2+ε1+γ21),εt=12(ε2+ε1−γ21),γlt=(−ε2+ε1), (10) 式中：
${\varepsilon _{{l}}}$ 、${\varepsilon _{{t}}}$ 分别为沿l、t方向钢筋的双轴应变；$ {\varepsilon _1} $ 、$ {\varepsilon _2} $ 分别为1、2方向的混凝土平均双轴正应变；$ {\gamma _{21}} $ 为1-2坐标系中混凝土的平均剪应变.假定混凝土薄膜元中应变是线性分布的，可得变形协调方程为
{ϕ=θsin2α2=θ,td=−ˉε2sϕ, (11) 式中：
$ \phi $ 为混凝土薄膜元应变梯度曲率；$ {\alpha _2} $ 为混凝土主压应力方向与纵筋的夹角，取45°；${\bar \varepsilon _{2{\rm{s}}}}$ 为2方向的混凝土表面单轴正应变.微元A处于双向受力状态，双轴应变与单轴应变之间的关系为
{ˉε1=ε1+ε2ν12,（12a）ˉε2=ε2,（12b）ˉεl=12(ˉε2+ˉε1+γ21),（12c）ˉεt=12(ˉε2+ˉε1−γ21),（12d） 式中：
$ {\bar \varepsilon _1} $ 、$ {\bar \varepsilon _2} $ 、${\bar \varepsilon _{{l}}}$ 、${\bar \varepsilon _{{t}}}$ 分别为对应于$ {\varepsilon _1} $ 、$ {\varepsilon _2} $ 、${\varepsilon _{{l}}}$ 、${\varepsilon _{{t}}}$ 的单轴应变值；$ {\nu _{12}} $ 为纯扭作用下的Hsu/Zhu比，如式(13).ν12={0.8(0.2+850εsf),εsf⩽εy,0.8×1.9,εsf>εy, (13) 式中：
${\varepsilon _{{\rm{sf}}}}$ 为纵筋与箍筋中先屈服者的应变；${\varepsilon _{\rm{y}}}$ 为钢筋屈服应变.2.1.3 材料本构方程
受压混凝土的本构方程如下：
σ2c=−ηk1cζfc, (14) k1c=(−Ecˉε1i2fc+1ζfcˉε2s∫ˉε2s0σ2c(¯ε2,ζ)d¯ε2), (15) ζ={0.9,¯ε1,con=¯ε1+¯ε1i⩽0,5.8√fc1√1+400¯ε1,con(1−|β|24°)⩽0.9且5.8√fc⩽0.9,¯ε1,con=¯ε1+¯ε1i>0, (16) β=12tan−1[γ21/(ε2−ε1)]180π, (17) 式（14）~（17）中：
$ \eta $ 为修正系数，按式（19）计算；${k_{1{\rm{c}}}}$ 为平均压应力系数：$ \zeta $ 为软化系数；$f_{\rm{c}}$ 为混凝土轴心抗压强度标准值；$\sigma _{2{\rm{c}}}(\bar \varepsilon ,\zeta )$ 按文献[19]分段选取；${\bar \varepsilon _{{\rm{1i}}}} = 0.625{\bar \varepsilon _{l{\rm{i}}}}$ ，为预应力引起的1方向混凝土初应变；$ \beta $ 为裂缝旋转角度.受拉混凝土的本构方程如下：
σ1c={12η(ˉEcε1s+σci),ˉˉε1s⩽ˉεcx,η[14(ˉεcx−ˉε1i)ˉˉε1sσci+12ˉˉEc(ˉε1s+ˉε1i−ˉεcx)(1−ˉεcx−ˉε1iˉε1s)],ˉεcx<ˉˉε1s⩽εcr,η{σci4(ˉεcx−ˉε1i)ˉε1s+fcr2(εcr−ˉεcx)ˉε1s+53fcrε0.4crˉε1s[(ˉε1s+ˉε1i)0.6−ε0.6cr]},ˉˉε1s>εcr, (18) 式中：
${\bar{\bar \varepsilon} _{1{\rm{s}}}} = {\bar \varepsilon _{1{\rm{s}}}} + {\bar \varepsilon _{{\rm{1i}}}}$ ，$ {\bar \varepsilon _{1s}} $ 为1方向的混凝土表面单轴正应变；${\varepsilon _{{\rm{cr}}}} = {0.000\;08} \mu$ ，为混凝土开裂应变；${f_{{\rm{cr}}}} = \lambda \mu \left( {0.31} \right)\sqrt {f_{\rm{c}}}$ ，为混凝土开裂应力，$ \lambda $ 、$ \mu $ 分别为受扭混凝土相较于受剪混凝土的弹性模量、开裂应变放大系数，如式（19）；${\bar \varepsilon _{{\rm{cx}}}} = {\bar \varepsilon _{{\rm{1i}}}} - \dfrac{{{\sigma _{{\rm{ci}}}}}}{{2 \bar E_{\rm{c}}}}$ 为混凝土消压应变；$\bar E_{\rm{c}} = \dfrac{{2f_{\rm{c}}}}{{{\varepsilon _0}}}$ ；$\bar{\bar E}_c = \dfrac{{{f_{{\rm{cr}}}}}}{{{\varepsilon _{{\rm{cr}}}} - {{\bar \varepsilon }_{{\rm{cx}}}}}}$ .{λ=μ=0.93,η=0.033√fc+0.73,tc⩽0.91td,λ=μ=1.20,η=0.094√fc+0.43,tc>0.91td且fc⩽47.85MPa,λ=μ=1.129,η=8.45/√fc+0.017,tc>0.91td且fc>47.85MPa. (19) 普通钢筋的本构方程如下：
fs={Esˉεs,ˉεs⩽ˉεn,fy[(0.91−2B)+(0.02+0.25B)ˉεsεy],ˉεs>ˉεn, (20) {ˉεn=εy(0.93−2B),B=1ρ(fcrfy)1.5, (21) 式中：
${f_{\rm{s}}}$ 为嵌入普通钢筋的应力；${f_{\rm{y}}}$ 为裸筋的屈服应力；${\bar \varepsilon _{\rm{s}}}$ 、${\bar \varepsilon _{\rm{n}}}$ 为嵌入钢筋的单轴应变和单轴屈服应变.对于嵌入纵筋，
${\bar \varepsilon _{{{ls}}}} = {\bar \varepsilon _{{l}}} + {\bar \varepsilon _{l{\rm{i}}}}$ .预应力钢筋的本构方程如下：
fp={Epεp,εp⩽εpro,Epεp[1+(Epεpfpu)4.38]−14.38,εp>εpro, (22) 式中：fp为预应力筋的应力；
${E_{\rm{p}}}$ 、${\varepsilon _{\rm{p}}}$ 、${\varepsilon _{{\rm{pro}}}}$ 、${f_{{\rm{pu}}}}$ 分别为预应力筋的弹性模量、应变、比例极限应变和极限强度，${\varepsilon _{\rm{p}}} = {f_{{\rm{pi}}}}/{E_{\rm{p}}} + {\bar \varepsilon _{{l}}}$ .波形钢腹板的本构方程如下：
{τw={Geγw,τw⩽τwy,τy,τw>τwy,Ge=Gs(aw+bw)/(aw+cw), (23) 式中：
${G_{\rm{e}}}$ 为波形钢腹板的有效剪切模量；${G_{\rm{s}}}$ 为钢材的剪切模量；${a_{\rm{w}}}$ 、${b_{\rm{w}}}$ 和${c_{\rm{w}}}$ 分别为波形钢腹板的直腹板段、斜腹板段投影和斜腹板段长度.混凝土剪应力与剪应变之间的关系为
τ21c=12(σ1c−σ2c)ε1−ε2γ21+σci2. (24) 需要指出的是，一些已有的PCCBGCSWs纯扭分析模型（如文献[9-10]）忽略了混凝土的抗拉强度，因此仅能预测结构在混凝土开裂后的抗扭行为. 而ISMMT成功地将混凝土抗拉强度考虑在内，从而能够预测PCCBGCSWs包括开裂前上升段在内的整个扭矩-扭率曲线.
2.2 求解程序框图
2.2.1 通用求解程序
联立式（9）、（11）可求得
${t_{\rm{d}}}$ 的表达式，但可能出现计算所得${t_{\rm{d}}}$ 大于混凝土翼缘板厚度${t_{\rm{c}}}$ 的不实际情况. 因此，对${t_{\rm{d}}}$ 进行判定并修正，如式（25）.td={−2ˉε2bh(γltb+γwhw−2ˉε2b),td⩽tc,tc,td>tc. (25) 式（25）的推导是为了避免采用“试错法”迭代计算td，从而提高求解效率. ISMMT的通用求解程序见图2（a），式（26）~（28）为迭代计算的3个收敛判别准则.
ρlfl+ρpfp+ρtft−ρlifli−ρpifpi=−(σ2c+σ1c−σci), (26) ρlfl+ρpfp−ρtft−ρlifli−ρpifpi=−(2τ21c−σci), (27) {τlttd=τwtw,τw<τwy,γw=γlt,τw=τwy. (28) 2.2.2 简化求解程序
Jeng等[17]给出的SMMT求解程序适用于PC构件的纯扭全过程分析，但其中包含2层嵌套迭代循环，较为复杂，求解耗时也较长. 为提高求解效率，Jeng等在此基础上提出了当普通钢筋处于弹性阶段时的简化求解程序，该简化程序仅包含1层迭代循环[22]. 同样地，图2（a）所示的通用求解程序适用于PCCBGCSWs纯扭全过程分析，但包含3层嵌套迭代循环，求解更为复杂与耗时. 因此，本文在充分考虑PCCBGCSWs结构特点的基础上，参考Jeng等[17]提出的简化模型，提出了当波形钢腹板、预应力及普通钢筋均处于弹性阶段时的简化程序框图（图2（b））. 当纵筋、箍筋及预应力钢筋均处于弹性阶段时，其本构关系可分别表示为
{fl=Es(ˉεl+ˉεli),ft=Esˉεt,fp=Ep(ˉεl+εpi). (29) 将式（12c）、（12b）、（24）、（29）代入式（26）和式（27），可得
AlEs2btd(ˉεli+12(ˉε2+ˉε1+γ21))+AtEs2std(ˉε2+ˉε1−γ21)+ApEp2btd(εpi+12(ˉε2+ˉε1+γ21))−ρlifli−ρpifpi=−(σ2c+σ1c−σci), (30) AlEs2btd(ˉεli+12(ˉε2+ˉε1+γ21))−AtEs2std(ˉε2+ˉε1−γ21)+ApEp2btd(εpi+12(ˉε2+ˉε1+γ21))−ρlifli−ρpifpi=(σ2c−σ1c)ε1−ε2γ21. (31) 显然，
$ {\gamma _{21}} $ 可通过求解式（30）、（31）得出，求解$ {\gamma _{21}} $ 的目的是为了消除与$ {\gamma _{21}} $ 相关联的迭代循环. 此外，由于波形钢腹板也处于弹性阶段，与${\gamma _{\rm{w}}}$ 相关联的迭代循环可通过代入${\gamma _{\rm{w}}} = {{{\tau _{{{lt}}}}{t_{\rm{d}}}} \mathord{\left/ {\vphantom {{{\tau _{{\rm{lt}}}}{t_{\rm{d}}}} {{G_{{\rm{es}}}}{t_{\rm{w}}}}}} \right. } {{G_{{\rm{es}}}}{t_{\rm{w}}}}}$ （根据式（23）和$q = {\tau _{{{lt}}}}{t_{\rm{d}}} = {\tau _{\rm{w}}}{t_{\rm{w}}}$ 求出）予以消除. 至此，通过上述简化，求解程序仅包含1层迭代循环，如图2（b）所示. 需要注意的是，图2（a）所示的通用求解程序适用于纯扭全过程分析，而图2（b）所示的简化求解程序仅适用于当波形钢腹板、预应力及普通钢筋均处于弹性阶段时的纯扭分析.与通用求解程序相比，该简化求解程序的计算效率有显著提高，主要原因为：1） 通用求解程序与简化求解程序分别包含3层嵌套迭代循环和1层迭代循环，其计算复杂度可分别用O(
$n^3 $ )和O(n)来表示（n为问题的维度）. 因此，后者的求解效率要显著高于前者；2） 通用求解程序采用“试错法”来假定未知变量$ {\gamma _{21}} $ 、${\varepsilon _1}$ 和${\gamma _{\rm{w}}}$ 的值，具有盲目性，求解效率低. 而在简化求解程序中，当收敛判别式不满足后，程序可自动生成下一个迭代循环所需的初始值$ {\bar \varepsilon _2} $ ，从而避免采用“试错法”盲目寻找，可有效提高求解效率. 同时，与本文通用求解程序类似，已有的同类型分析模型（文献[9-12、14-16、18]中模型）均包含3层嵌套迭代循环. 因此，本文简化求解程序的求解效率同样远高于已有同类型分析模型.根据图2所示求解程序框图，本文采用MATLAB编制了PCCBGCSWs纯扭全过程分析模型，运行程序可得到一组扭矩与扭率的离散点，最终绘成扭矩-扭率图.
2.3 ISMMT与已有分析模型比较
与文献中已有的同类型分析模型相比，ISMMT有着较为显著的优势. 首先，由于ISMMT基于单一软化薄膜元理论提出，因而具有理论一致性. 同时，SMMT成功地将混凝土抗拉强度、抗拉刚化效应以及开裂混凝土提供的抗剪强度考虑在内. 本文ISMMT基于SMMT提出，同样能够考虑上述效应，因而能准确模拟结构包括开裂前阶段在内的全过程扭转行为. 此外，与已有分析模型[9-12,14-16,18]的求解程序相比，本文提出的简化求解程序效率更高、求解耗时更少. 因此，ISMMT为纯扭作用下PCCBGCSWs的受力分析提供了有效途径.
3. 模型试验
3.1 试件概况
设计并制造了一根等截面PCCBGCSWs试件，如图3所示. 图3（a）为试件的立面图，梁的长度及高度分别为2.8 m和0.55 m， 为方便固定与加载，在梁的两端各布置一道0.4 m厚的横隔板. 图3（b）为组合箱梁横截面尺寸和预应力及普通钢筋布置示意图，普通钢筋型号均为HRB335，规格均为Φ10.0，箍筋间距为100 mm，混凝土保护层厚度均为20 mm. 预应力钢束采用单根直径15.2 mm、抗拉强度标准值1 860 MPa的低松弛钢绞线，在加载前施加初应力. 波形钢腹板采用Q235钢材，图3（c）为其尺寸示意图. 为了确保波形钢腹板与混凝土顶、底板连接处在加载过程中不提前发生破坏，剪力连接件采用刚度较大的双PBL形式. 此外，模型试件的材料特性列于表1.
表 1 试件材料参数Table 1. Material properties of the test beam混凝土棱柱体抗压强度/MPa 波形钢腹板 预应力钢筋初应力/MPa 普通钢筋 厚度/mm 屈服强度/MPa 直径/mm 屈服强度/MPa 24.9 3.7 235 800 10 357 3.2 加载装置
自行设计了一套纯扭加载装置，如图4所示. 试件的一端通过千斤顶进行固定，另一端通过加载梁夹紧后作为转动端，在转动端底部放置铰支座. 在压力传感器与反力梁之间放置球铰，采用液压千斤顶在加载梁上施加竖向偏心荷载，可使得试件在加载过程中处于近乎纯扭的受力状态. 加载梁的悬臂段长1.65 m，通过压力传感器对竖向荷载的大小进行控制. 由于加载梁自重较大，考虑加载梁偏心对试件产生的初始扭矩.
3.3 测试方案
3.3.1 扭矩与扭转角
扭矩可近似按竖向偏心荷载与加载梁力臂的乘积计算，该近似值与考虑扭转角的精确值相比相差较小[23]. 此外，在试件的1/2及3/4截面分别布置位移千分表来测量试件在加载过程中的扭转竖向位移. 按式（32）换算得到该截面的扭转角，图5为计算示意图.
[L2c+(δc−h2)2]sin2αc−Lchsinαc−δ2c+δch=0, (32) 式中：
${L_{\rm{c}}}$ 为梁中心至千分表的距离；${\delta _{\rm{c}}}$ 为千分表测得的扭转竖向位移；${\alpha _{\rm{c}}}$ 为试件的扭转角度.3.3.2 测点布置
试件的测点布置情况为：在跨中截面左、右两侧波形钢腹板上各布置6组（共12组）应变花，编号分别为LW-1~LW-6（左侧）和RW-1~RW-6（右侧）；跨中截面混凝土顶、底板沿横向分别均匀布置5组（共10组）应变花，编号分别为TC-1~TC-5（顶板）和BC-1~BC-5（底板）；预应力钢束采用后张法进行张拉（未灌浆，属无黏结预应力），沿对角选取2根钢束在其张拉端布置预应力传感器，编号分别为P-ZS及P-YX（如图4所示）；在混凝土顶、底板内各选取8个普通钢筋（包括纵筋和箍筋）测点（共16个）布置应变片，编号分别为TG-1~TG-8（顶板）及BG-1~BG-8（底板）. 上述主要测点的详细布置情况如图6所示.
4. 结果对比
为验证ISMMT能否准确预测PCCBGCSWs中各构件在纯扭作用下的应变历程，将ISMMT预测的各项理论结果与对应的试验结果进行对比，如图7所示.
4.1 扭矩-扭率曲线
将扭矩-扭率曲线的理论和试验结果列于图7（a），其中，扭率的试验结果取所有4个竖向位移测点按式（27）换算得到扭率值的平均值. 从图7（a）可以看出：当试件处于弹性阶段时，扭率试验结果要略大于对应的理论结果，造成该现象的原因可能是该阶段扭率值相对较小，千分表读数易受周围环境的干扰. 但由于本文仅完成了1根模型梁的纯扭试验，扭率的试验与理论结果存在偏差的原因还有待进一步研究. 在混凝土开裂后，理论与试验曲线遵循着相同的发展趋势，且吻合良好. 此外，将理论和试验曲线在纯扭加载过程中关键节点处的对应结果列于表2~4. 表中：T 和 θ 分别为扭矩和扭率，下标第1个数（1、2、3）代表状态（开裂状态、屈服状态和极限状态），下标第2个数1、2分别代表理论值与试验值. 表中结果表明：除开裂状态下的扭率试验值与对应的理论结果相差较大外，其他结果均十分接近. 由此可知：ISMMT能准确预测PCCBGCSWs在纯扭状态下的全过程扭矩-扭率曲线.
4.2 波形钢腹板剪应变
图7（b）为由ISMMT和试验给出的波形钢腹板剪应变结果（绝对值）. 从图中可以看出：理论与试验曲线吻合良好，且由分析模型给出波形钢腹板平均剪应变位于各分散试验曲线的区间内；左侧与右侧波形钢腹板对应位置处的剪应变值很接近，且在凹进（LW-1~LW4、RW-1~RW4）与凸出（LW-5~LW6、RW-5~RW6）直钢板处的剪应变值也相差很小. 在钢腹板屈服前，波形钢腹板上的剪应变随着外扭矩的增加而近乎线性增加. 而当钢腹板达到其屈服强度时，由于模型梁已基本达到其抗扭极限承载力，因此尽管波形钢腹板上的剪应变迅速增加，施加在试件上的外扭矩仍基本保持不变.
4.3 混凝土翼缘板剪应变
图7（c）给出了混凝土翼缘板剪应变（绝对值）的理论和试验结果对比情况. 图中结果表明：当模型梁处于弹性阶段时，混凝土翼缘板上的剪应变随着外扭矩的增加而近乎线性增加，但增长幅度很小；在该阶段内，理论结果与试验结果高度吻合；在混凝土开裂后，混凝土顶、底板上的剪应变迅速增加，尽管此时试验曲线较为分散，但由ISMMT给出混凝土翼缘板平均剪应变仍处于各分散试验曲线的区间内，且二者遵循着同样的变化趋势. 造成精度下降的主要原因在于混凝土裂缝出现位置及开裂角度的不确定性.
图 7 ISMMT预测结果与试验结果对比Figure 7. Comparison of the results obtained from the ISMMT and experiment表 2 ISMMT预测结果与试验结果对比（开裂状态）Table 2. Comparison of the predicted torques and twists from the ISMMT and the experiment (cracking state)参数 $ {T_{11}}/$
(kN•m)${T_{12} }/$
(kN•m)$\dfrac{ { {T_{11} } } }{ { {T_{12} } } }$ ${\theta _{ 11 } }/$
((°)•m−1)${\theta _{ 12 } }/$
((°)•m−1)$\dfrac{ { {\theta _{11} } } }{ { {\theta _{12 } } } }$ 取值 160.6 174.5 0.92 0.073 0.161 0.45 表 3 ISMMT预测结果与试验结果对比（屈服状态）Table 3. Comparison of the predicted torques and twists from the ISMMT and the experiment (yield state)参数 ${T_{ 21 } }/$
(kN•m)${T_{22 } }/$
(kN•m)$\dfrac{ { {T_{21} } } }{ { {T_{22 } } } }$ ${\theta _{ 21 } }/$
((°)•m−1)${\theta _{ 22} }/$
((°)•m−1)$\dfrac{ { {\theta _{ 21 } } } }{ { {\theta _{ 22 } } } }$ 取值 351.9 337.9 1.04 1.043 0.973 1.07 表 4 ISMMT预测结果与试验结果对比（极限状态）Table 4. Comparison of the predicted torques and twists from the ISMMT and the experiment (limit state)参数 ${T_{ 31 } }/$
(kN•m)$ {T_{32}} /$
(kN•m)$\dfrac{ { {T_{ 31} } } }{ { {T_{ 32 } } } }$ $ {\theta _{31}} /$
((°)•m−1)$ {\theta _{32 }} /$
((°)•m−1)$\dfrac{ { {\theta _{31} } } }{ { {\theta _{32} } } }$ 取值 354.9 339.5 1.04 1.547 1.590 0.97 4.4 预应力钢束应变
将由ISMMT和试验给出的预应力钢束应变结果列于图7（d）. 图中结果表明：理论与试验曲线吻合较好且遵循着同样的发展规律，在混凝土开裂前，预应力钢束的应变基本保持不变；当所施加的外扭矩达到开裂扭矩后，混凝土顶底板上斜裂缝的数量和宽度随着外扭矩的增加而不断扩张，从而导致预应力钢束应变迅速增加.
4.5 纵筋及箍筋应变
图7（e）、（f）分别给出了纵筋和箍筋应变的理论和试验结果对比. 从图中可以看出：由于试验误差等原因，纵、箍筋各测点的应变试验结果较为离散，但总体保持着一致的发展趋势，且理论与试验结果也较为吻合；预应力的存在会使得模型梁沿桥轴向产生初始应变. 因此，在正式施加纯扭荷载前，纵筋处于受压状态（如图7（e））；图7（d） 中预应力钢束应变的变化趋势和幅值与图7（e）中纵筋应变一致，该试验现象与理论模型中的假设吻合，从而进一步证明了理论模型的准确性.
4.6 求解效率
当采用ISMMT预测试件的纯扭全过程受力行为时，在波形钢腹板、预应力及普通钢筋均处于弹性阶段时分别采用简化求解程序和通用求解程序进行计算，并对二者的求解耗时进行对比. 其中，两种求解程序的收敛误差均控制在0.1%以内. 结果表明，对于某指定的
$ {\varepsilon _2} $ ，采用通用求解程序求得$ {\gamma _{21}} $ 、${\varepsilon _1}$ 和$ {\gamma _{\rm{w}}} $ 真实值所需的总迭代次数最少为4.9 × 105次，最多达7.9 × 106次（需要说明的是，该迭代次数与$ {\gamma _{21}} $ 、${\varepsilon _1}$ 和$ {\gamma _w} $ 所假定的迭代步长有关）. 而采用简化求解程序进行计算时，对于某指定的${\bar \varepsilon _1}$ ，求得所有未知变量的最少迭代次数为5次，最多也仅为193次. 由此可见，采用本文提出的简化求解程序可以极大地提高求解效率.5. 结 论
在软化薄膜元理论的基础上提出了改进软化薄膜元模型（ISMMT）来预测PCCBGCSWs在纯扭作用下的全过程受力行为. 针对ISMMT以及已有同类型理论模型中求解程序所存在的迭代循环多、求解效率低等问题，提出了当波形钢腹板、预应力及普通钢筋均处于弹性阶段时的简化求解程序框图. 此外，还完成了纯扭模型试验对ISMMT的适用性与准确性进行验证. 主要结论如下：
1） 本文所提出的简化求解程序框图仅有一层迭代循环，其求解效率要远高于国内外已有同类型分析模型.
2） 由模型试验得到的试件扭矩-扭率曲线、波形钢腹板和混凝土翼缘板剪应变、预应力及普通钢筋应变与ISMMT预测的理论结果吻合良好，表明ISMMT除了能准确预测PCCBGCSWs的全过程扭矩-扭率曲线外，还能模拟混凝土翼缘板、波形钢腹板、预应力和普通钢筋等构件的整个应变发展历程. 该理论模型为更全面地了解PCCBGCSWs的纯扭力学性能提供了有效途径.
-
图 5 应变片布置
1—荷载传感器；2—分配工装；3—PP—ECC梁侧面；4—铰支座；5—应变片；6—百分表；7—PP-ECC梁顶面.
Figure 5. Strain gauge layout
图 9 冻融过程中PP-ECC的极限压缩性能变化
Figure 9. Variations in ultimate compression performance of PP-ECC during freeze-thaw process
图 10 冻融过程中PP-ECC的拉伸性能变化
Figure 10. Variations in tensile properties of PP-ECC during freeze-thaw process
表 1 PP-ECC梁具体试验参数
Table 1. Specific test parameters of PP-ECC beam
试验梁
分组试件编号 冻融循环
次数/次弯曲
荷载/%抗弯性能
试验第Ⅰ组 E0-0 0 是 第Ⅱ组 E0-300 300 E0-500 500 第Ⅲ组 E0.25-300 300 25 E0.25-500 500 第Ⅳ组 E0.50-300 300 50 E0.50-500 500 
下载: 导出CSV
表 2 钢筋性能指标
Table 2. Rebar performance indexes
型号 直径/
mm截面面
积/mm2屈服强
度/MPa极限强
度/MPa弹性模
量/GPaQ235 6 28.26 236 390 196 HRB400 6 28.26 406 576 197 HRB400 10 78.50 411 535 203
下载: 导出CSV
表 3 PP-ECC梁的配合比
Table 3. Mix proportion of PP-ECC beam
kg/m3 水泥 粉煤灰 PP纤维 水 820.0 442.0 18.2 504.8
下载: 导出CSV
表 4 PP-ECC梁四点弯曲试验结果
Table 4. Four-point bending test results of PP-ECC beam
试件编号 屈服荷载/kN 极限荷载/kN 极限抗弯承载力/（kN·m） 裂缝最大宽度/mm 最终裂缝数量/条 E0-0 33.95 33.97 5.10 0.90 36 E0-300 29.34 28.84 4.33 1.10 30 E0-500 24.82 24.22 3.63 1.34 25 E0.25-300 30.88 30.72 4.61 1.25 32 E0.25-500 24.64 25.02 3.75 1.42 29 E0.50-300 26.67 26.01 3.90 1.55 23 E0.50-500 21.21 21.87 3.28 1.80 18
下载: 导出CSV
表 5 PP-ECC材料本构参数
Table 5. Constitutive parameters of PP-ECC materials
PP-ECC 单轴拉伸本构参数 单轴受压本构参数 σtc/MPa εtc/% σtu/MPa εtu/% σcu/MPa εcu/% 0.86 0.088 2.33 3.789 31.65 0.38
下载: 导出CSV
表 6 抗弯承载力对比
Table 6. Comparison of flexural bearing capacity
冻融循环
次数/次抗弯承载力 计算值/
试验值计算值/
（kN·m）试验值/
（kN·m）0 4.50 5.10 0.88 300 4.18 4.33 0.97 500 3.88 3.63 1.06
下载: 导出CSV
表 7 抗弯承载力对比
Table 7. Comparison of flexural bearing capacity
冻融循环
次数/次持荷比 抗弯承载力 计算值/
试验值计算值/
（kN·m）试验值/
（kN·m）300 0 4.18 4.33 0.97 0.25 4.46 4.61 0.97 0.50 3.75 3.90 0.96 500 0 3.88 3.63 1.07 0.25 4.14 3.75 1.10 0.50 3.48 3.28 1.06
下载: 导出CSV
-
[1] CAMPIONE G, CANNELLA F, CAVALERI L. Shear and flexural strength prediction of corroded R. C. beams[J]. Construction and Building Materials, 2017, 149: 395-405. doi: 10.1016/j.conbuildmat.2017.05.125 [2] LI V C, LEUNG C K Y. Steady-state and multiple cracking of short random fiber composites[J]. Journal of Engineering Mechanics, 1992, 118(11): 2246-2264. doi: 10.1061/(ASCE)0733-9399(1992)118:11(2246) [3] LI V C, WU H C. Conditions for pseudo strain-hardening in fiber reinforced brittle matrix composites[J]. Applied Mechanics Reviews, 1992, 45(8): 390-398. doi: 10.1115/1.3119767 [4] YUAN F, PAN J L, DONG L T, et al. Mechanical behaviors of steel reinforced ECC or ECC/concrete composite beams under reversed cyclic loading[J]. Journal of Materials in Civil Engineering, 2014, 26(8): 04014047.1-04014047.8. [5] POURFALAH S. Behaviour of engineered cementitious composites and hybrid engineered cementitious composites at high temperatures[J]. Construction and Building Materials, 2018, 158: 921-937. doi: 10.1016/j.conbuildmat.2017.10.077 [6] JIN H S, LI F H, HU D H. Research on the flexural performance of reinforced engineered cementitious composite beams[J]. Structural Concrete, 2022, 23(4): 2198-2220. doi: 10.1002/suco.202100012 [7] 黄士元. 近代混凝土技术[M]. 西安: 陕西科学技术出版社, 1998. [8] SUN W, ZHANG Y M, YAN H D, et al. Damage and damage resistance of high strength concrete under the action of load and freeze-thaw cycles[J]. Cement and Concrete Research, 1999, 29(9): 1519-1523. doi: 10.1016/S0008-8846(99)00097-6 [9] KOSIOR-KAZBERUK M, BERKOWSKI P. Surface scaling resistance of concrete subjected to freeze-thaw cycles and sustained load[J]. Procedia Engineering, 2017, 172: 513-520. doi: 10.1016/j.proeng.2017.02.060 [10] SHEN Y, LIU J, ZHOU S Y, et al. Experimental investigation on the freeze-thaw durability of concrete under compressive load and with joints[J]. Construction and Building Materials, 2019, 229: 116893.1-116893.12. doi: 10.1016/j.conbuildmat.2019.116893 [11] WANG Y, WANG G, GUAN Z W, et al. The effect of freeze-thaw cycles on flexural behaviour of FRP-reinforced ECC beams[J]. Archives of Civil and Mechanical Engineering, 2021, 21(3): 1-24. [12] ZHANG Y X, BAI S, ZHANG Q B, et al. Failure behavior of strain hardening cementitious composites for shear strengthening RC member[J]. Construction and Building Materials, 2015, 78: 470-473. doi: 10.1016/j.conbuildmat.2015.01.037 [13] DING Y, YU K Q, YU J T, et al. Structural behaviors of ultra-high performance engineered cementitious composites (UHP-ECC) beams subjected to bending-experimental study[J]. Construction and Building Materials, 2018, 177: 102-115. doi: 10.1016/j.conbuildmat.2018.05.122 [14] YUN H D, KIM S W, LEE Y O, et al. Tensile behavior of synthetic fiber-reinforced strain-hardening cement-based composite (SHCC) after freezing and thawing exposure[J]. Cold Regions Science and Technology, 2011, 67(1/2): 49-57. [15] 靳贺松,李福海,何肖云峰,等. 聚丙烯纤维水泥基复合材料的抗冻性能研究[J]. 材料导报,2020,34(8): 8071-8076,8082. doi: 10.11896/cldb.18080117JIN Hesong, LI Fuhai, HE Xiaoyunfeng, et al. Research on frost resistance of polypropylene fiber cement-based composite material[J]. Materials Reports, 2020, 34(8): 8071-8076,8082. doi: 10.11896/cldb.18080117 [16] DUAN A, LI Z Y, ZHANG W C, et al. Flexural behaviour of reinforced concrete beams under freeze-thaw cycles and sustained load[J]. Structure and Infrastructure Engineering, 2017, 13(10): 1350-1358. doi: 10.1080/15732479.2016.1268172 [17] 李福海,胡丁涵,余泳江,等. PP-ECC梁抗弯性能试验研究[J]. 西南交通大学学报,2021,56(2): 272-281.LI Fuhai, HU Dinghan, YU Yongjiang, et al. Experimental study on flexural capacity of PP-ECC beam[J]. Journal of Southwest Jiaotong University, 2021, 56(2): 272-281. [18] 崔晶波. 荷载冻融耦合作用下再生粗骨料自密实混凝土梁的受弯力学性能试验研究[D]. 乌鲁木齐: 新疆大学, 2020. [19] IAN T. Research on seismic behavior of steel tube reinforeed conerete bridge columns[D]. Dalian: Dalian University of Technology, 2019. [20] LI F H, WEN T, LI J Y, et al. Ultrasonic-detected damage and bending behavior of reinforced PP-ECC beams after coupled action of freeze-thaw cycles and constant flexural load[J]. Case Studies in Construction Materials, 2022, 17: e01284.1-e01284.15. [21] 徐港,李运攀,潘琪,等. 盐冻环境下钢筋混凝土梁抗弯性能试验研究[J]. 土木建筑与环境工程,2014,36(3): 86-91.XU Gang, LI Yunpan, PAN Qi, et al. Experimental analysis on flexural performance of reinforced concrete beam in salt-frost environment[J]. Journal of Civil, Architectural & Environmental Engineering, 2014, 36(3): 86-91. [22] 刘泽军,李艳,温丛格. PVA-ECC劈裂抗拉强度与变形性能试验研究[J]. 建筑材料学报,2016,19(4): 746-751. doi: 10.3969/j.issn.1007-9629.2016.04.024LIU Zejun, LI Yan, WEN Congge. Experimental study on strength and deformation performance of PVA-ECC under splitting tension[J]. Journal of Building Materials, 2016, 19(4): 746-751. doi: 10.3969/j.issn.1007-9629.2016.04.024 -
下载:
下载:
点击查看大图
计量
- 文章访问数: 308
- HTML全文浏览量: 121
- PDF下载量: 61
- 被引次数: 0