GFRP管型钢再生混凝土组合柱轴压性能试验及数值分析
0 引言
再生混凝土是由废弃混凝土经过破碎处理后,按一定比例与级配混合,再掺入水泥、水等而得到的混凝土
为改善再生混凝土结构的力学性能,部分学者充分利用组合结构的力学性能优势,提出了GFRP管再生混凝土柱与型钢再生混凝土柱
为研究GFRP管型钢再生混凝土组合柱的轴压性能,本文采用试验研究和有限元分析相结合的方法,对该组合柱轴心受压性能进行分析,重点研究其破坏形态及特征、荷载-位移曲线等,在此基础上,重点分析设计参数对该组合柱轴压性能的影响,研究结果可为GFRP管型钢再生混凝土组合柱的工程应用提供技术参考。
1 试验概况
设计制作了11根GFRP管型钢再生混凝土组合柱试件,主要设计参数为再生粗骨料取代率、长细比、型钢配钢率及再生混凝土强度,试件具体参数见表1。GFRP管直径为200mm,壁厚为10mm,产自江苏某电力股份有限公司,试件几何尺寸如图1所示。再生粗骨料来源于单一的拆迁废弃混凝土,其基本物理性能指标满足《混凝土用再生粗骨料》(GB/T 25177—2010)规定要求,天然粗骨料选用人工碎石,细骨料为中粗河砂。再生混凝土材料配合比及立方体抗压强度指标见表2。型钢为Q235低碳钢,其基本力学性能指标如表3所示。
试件设计参数汇总 表1
试件 编号 |
再生混 凝土强 度等级 |
型钢尺寸/mm |
再生粗 骨料取 代率r/% |
柱高 H/mm |
长细 比λ |
配钢率 ρ/% |
GSRC-1 |
C40 | 100×68×4.5×7.6 | 0 | 500 | 10 | 4.54 |
GSRC-2 |
C40 | 100×68×4.5×7.6 | 30 | 500 | 10 | 4.54 |
GSRC-3 |
C40 | 100×68×4.5×7.6 | 50 | 500 | 10 | 4.54 |
GSRC-4 |
C40 | 100×68×4.5×7.6 | 70 | 500 | 10 | 4.54 |
GSRC-5 |
C40 | 100×68×4.5×7.6 | 100 | 500 | 10 | 4.54 |
GSRC-6 |
C40 | 126×74×5×8.4 | 100 | 500 | 10 | 5.76 |
GSRC-7 |
C40 | 140×80×5.5×9.1 | 100 | 500 | 10 | 6.85 |
GSRC-8 |
C40 | 100×68×4.5×7.6 | 100 | 1 200 | 24 | 4.54 |
GSRC-9 |
C40 | 100×68×4.5×7.6 | 100 | 1 800 | 36 | 4.54 |
GSRC-10 |
C50 | 100×68×4.5×7.6 | 100 | 500 | 10 | 4.54 |
GSRC-11 |
C60 | 100×68×4.5×7.6 | 100 | 500 | 10 | 4.54 |
注:λ=l0/i,
GFRP管型钢再生混凝土组合柱轴心受压试验在西安理工大学结构工程实验室进行,采用500t电液伺服压力机,试件位移与应变测点布置如图2所示。图3为试件试验加载装置。
再生混凝土配合比及立方体抗压强度指标 表2
再生混 凝土强 度等级 |
再生粗骨 料取代率 r/% |
水胶比 W/B |
单位体积用量/(kg/m3) |
立方体 抗压强度 fcs/MPa |
轴心抗 压强度 fct/MPa |
弹性模量 Ec/MPa |
||||||
水泥 |
砂 |
天然 粗骨料 |
再生 粗骨料 |
水 | 粉煤灰 | 减水剂 | ||||||
C40 |
0 | 0.44 | 443 | 576 | 1 171 | 0 | 195 | 0 | 0 | 43.37 | 32.96 | 2.685×104 |
30 |
0.45 | 443 | 576 | 819.5 | 351.3 | 198.5 | 0 | 0 | 43.26 | 32.88 | 2.683×104 | |
50 |
0.45 | 443 | 576 | 585.5 | 585.5 | 200.8 | 0 | 0 | 42.26 | 32.12 | 2.667×104 | |
70 |
0.46 | 443 | 576 | 351.3 | 819.7 | 203.2 | 0 | 0 | 41.38 | 31.45 | 2.653×104 | |
100 |
0.47 | 443 | 576 | 0 | 1 171 | 206.7 | 0 | 0 | 40.01 | 30.41 | 2.630×104 | |
C50 | 100 | 0.36 | 358 | 649 | 0 | 1 138 | 163 | 94 | 3.5 | 51.61 | 39.22 | 2.796×104 |
C60 |
100 | 0.31 | 422 | 528 | 0 | 1 072 | 164.5 | 105.4 | 6.3 | 61.58 | 46.80 | 2.898×104 |
注:fct为再生混凝土轴心抗压强度,fct=0.76fcs;fcs为再生混凝土立方体抗压强度;Ec为再生混凝土弹性模量,Ec=1×105/(2.8+40.1/fcs)。
型钢材料力学性能指标 表3
材料 |
屈服强度 fy/MPa |
极限强度 fu/MPa |
弹性模量 Es/MPa |
||
型钢 |
翼缘 |
304 | 374 | 2.02×105 | |
腹板 |
283 | 363 | 1.99×105 |
2 主要试验结果分析
2.1 试件破坏形态
从试验过程观察可知,GFRP管型钢再生混凝土组合柱试件在轴压荷载作用下的破坏过程较为相似,本文取不同长细比的GSRC-3,GSRC-8以及GSRC-9试件分析其破坏特征,如图4所示。短柱以典型的GSRC-3试件为例,在加载初期,试件无明显变形,其应力-应变曲线呈线性变化;当加载至峰值荷载的60%时,试件中部略微发生鼓起,GFRP管应变增大,并伴有玻璃纤维撕裂声;当荷载达到峰值荷载的80%左右时,试件刚度出现退化,承载力增加缓慢,其应力-应变曲线呈非线性变化,试件中部鼓起变形明显;当加载至峰值荷载时,外部玻璃纤维管斜向撕裂,试件承载力急剧下降,伴有较大响声,组合柱发生整体破坏。对于长细比较大的GSRC-8与GSRC-9试件,加载初期,试件处于弹性阶段,试件无明显变形,其应力-应变曲线呈线性变化;当加载至峰值荷载的50%时,组合柱的GFRP管中部表面有肉眼可见的局部白纹;当加载至峰值荷载的70%左右时,试件中部鼓起明显,并不断有纤维撕裂声响起;随着荷载的不断增大,试件中部鼓起愈加明显,响声频率加快;达到峰值荷载后,GFRP管中部不规则撕裂并迅速向上下两端扩展,且伴有较大响声,随之试件发生破坏。
2.2 荷载-位移曲线
图5为GFRP管型钢再生混凝土组合柱的轴压荷载-位移曲线。由图可知各设计参数对组合柱荷载-位移曲线的影响。
(1)从图5(a)可知,在试件加载初期,组合柱处于弹性阶段,不同再生骨料取代率试件的荷载-位移曲线较为相似,曲线基本呈线性变化,再生骨料取代率的变化对试件的受力性能影响不大;加载到峰值荷载的60%左右,试件曲线发生明显变化,斜率逐渐降低,即刚度开始下降。试件初始刚度随着再生骨料取代率的增大而降低,其峰值承载力也随着降低,但不同取代率的极限位移相近。
(2)图5(b)为长细比对试件荷载-位移曲线的影响,加载初期,试件的初始刚度随着长细比的增大而降低,长细比大的试件其变形较明显;至加载后期,长细比较小的试件的荷载-位移曲线较为平缓,表明长细比较小的试件具有更好的延性;试件的峰值承载力随长细比的增大而降低。主要是因为组合中长柱受二阶效应影响较大,在加载过程中由于试件长细比较大而产生一定的附加偏心距,导致试件产生侧向弯曲,试件容易发生整体失稳使其承载力随之降低。
(3)图5(c)为不同型钢配钢率下试件承载力的变化规律,随着型钢配钢率的增大,试件的刚度以及承载力都随之增高,主要是由于在截面一定的情况下,增大型钢配钢率即增大截面轴压刚度。随着荷载增加,型钢配钢率较高的试件曲线也愈加平缓,说明具有较好的变形能力,即适当提高型钢配钢率对试件的延性和承载力有较大的提高。由此可知,增加型钢配钢率对组合柱的受力性能是有利的。
(4)由图5(d)可知,试件承载力随着再生混凝土强度等级的提高而增大,另外,随着再生混凝土强度等级的提高,试件的变形能力随之降低,这主要是由于强度等级较高的再生混凝土材料相对较脆。因此,在实际工程选择与材料相匹配的混凝土强度等级是至关重要的。
3 组合柱有限元模型的建立
3.1 材料的本构模型
3.1.1 GFRP材料本构关系
根据GFRP自身特点,GFRP是弹脆性材料,采用线弹性模型,破坏准则采用最大应变强度准则来近似模拟复合材料的损伤演化变形,且GFRP是一种特殊的各向正交异性材料-横观各向同性材料,在充分考虑复合材料的横观各向同性特点
GFRP材料工程弹性常数/GPa 表4
弹性常数 |
E1 | E2 | E3 | G | μ |
数值 |
10 | 10 | 18 | 6 | 0.25 |
3.1.2 再生混凝土本构关系
大量再生混凝土材料研究表明,由于再生粗骨料存在初始损伤及其质地不均匀的特点,使其力学性能普遍较低,故再生混凝土的材料力学性能与普通混凝土之间存在较大的差异;另外,核心再生混凝土由于受到GFRP管的被动约束作用,其受力为三向受压状态,其横向变形发展受到了限制,因此不能直接使用普通混凝土的本构关系。为满足适用ABAQUS软件分析的再生混凝土的本构关系,根据杜朝华等
其中:
式中 :x为再生混凝土应变变化;y为再生混凝土应力变化;ψ为再生骨料取代率影响曲率系数;ξ为约束效应; λ为长细比;fc为再生混凝土圆柱体抗压强度;γ为取代率。
3.1.3 钢材本构关系
试件采用型钢为Q235低碳钢,可作为弹塑性材料,采用von Mises屈服准则。钢材的应力-应变关系采用简化的弹塑性模型
式中:σi为钢材等效应力;fs为钢材屈服强度;fu为钢材极限强度;Es为钢材弹性模量;Es′为钢材切线模量,Es′=(fu-fs)/(εu-εst);εi为钢材等效应变;εy为钢材屈服时的应变;εu为钢材极限时的应变,εu=120εy;εst为钢筋硬化起点应变,εst=12εy,其中相关数据均由材性试验获得。
3.2 单元选取与网格划分
由于本试验采用的GFRP管厚度为10mm,不适宜用壳单元(S4R),因此GFRP管型钢再生混凝土组合柱中GFRP管、核心再生混凝土、型钢以及加载端板均采用八节点线性六面体实体单元(C3D8R)。整个模型网络采用结构化网络划分技术,网格划分密度对有限元计算非常重要,在模型生成时应结合网格试验来确定合理的网格密度,若在计算过程出现不收敛的情况,应重新调整网格密度以达到模型收敛的目的。GFRP管型钢再生混凝土组合柱有限元模型如图6所示。
3.3 材料界面相互作用关系
GFRP管、核心再生混凝土、型钢与端板均采用绑定(Tie)约束。在GFRP管型钢再生混凝土组合柱轴压过程中,组成部分变形会出现差异,故需要考虑它们界面之间的相互作用,GFRP管、型钢与再生混凝土界面之间的接触关系由法线方向的接触和切线方向的粘结滑移组成,法向方向选用“硬”接触,当两个界面脱离时即界面之间将不存在相互作用,切向力模拟采用库伦摩擦模型,GFRP管与再生混凝土、型钢与再生混凝土两个界面切向方向的作用力来源于界面之间的摩擦力,通过两个界面之间的粘结作用使得接触面之间的相对滑移为0,而当剪应力达到其临界值时,接触面之间就会出现相对滑移,临界应力的计算公式如下:
式中:μ为摩擦系数;p为界面法向压力。
GFRP与混凝土界面摩擦系数取0.6
3.4 边界条件及加载方式
参照试验边界条件情况,在GFRP管型钢再生混凝土组合柱一端施加固定约束,另一端为自由端,在自由端施加荷载。本次数值分析时采用位移控制加载方式。为保证轴向荷载能够均匀施加在试件一端,在试件两端分别绑定刚性端板;在柱自由端中心上部创建一个关键点,并使关键点与端板创建耦合约束,竖向荷载施加在关键点上,从而保证竖向荷载能够均匀地传递至短柱自由端。本文采用牛顿迭代法对上述有限元模型进行计算。
4 组合柱轴压有限元结果分析
4.1 应力云图与破坏形态
由于各试件的变形破坏情况较为相似,本文选取GSRC-3试件(短柱)和GSRC-9试件(长柱)的整体变形和应力云图为例进行说明。GSRC-3试件(短柱)和GSRC-9试件(长柱)的有限元分析和试验破坏形态的比较见图7。短柱试件主要是中部发生鼓起并最终破坏;中长柱试件是中部发生明显弯曲,有明显的侧向挠度。总体来说,有限元计算和试验得到的破坏形态吻合较好,该有限元模拟可以较好地反映组合柱在轴压下的变形特征。
通过在试验中观察GSRC-3试件(短柱)变形过程并结合有限元软件模拟的应力云图可知(图8),在试件加载初期,组合柱处于弹性阶段,试件整体受力较小,其应力-应变曲线近似呈线性增加;随着荷载的不断增大,轴向变形继续增大,试件各部分应力随之增大,此阶段,GFRP管对轴压强度的贡献很小,其应变增长缓慢,型钢的应变增长较快,试件无明显变形;当荷载增加到峰值荷载的60%左右时,型钢应变发展迅速,试件各部分应力分布有明显变化,试件中部出现轻微鼓起,刚度开始降低,此时型钢应变随荷载的增加呈非线性变化;当荷载增加至峰值荷载的80%左右时,型钢纵向应变增长速度明显大于荷载增长速度,型钢刚度下降速度明显加快,GFRP管的环向应变增长速度明显大于纵向应变增长速度,此时内部再生混凝土膨胀明显,即将被压碎,对型钢失去约束作用,GFRP管对再生混凝土约束作用加强,试件变形明显;随着荷载增加至峰值,型钢完全屈服,承载力不再增加,组合柱的压缩变形继续增大直至再生混凝土被压碎,GFRP管破裂,最终因材料破坏导致试件承载力急速下降。
由GSRC-9试件(长柱)的应力云图(图9)可知,长柱与短柱在前期的变形过程相似,在加载初期,试件处于弹性阶段,试件受力较小,此时GFRP管的约束作用很小,其纵向、环向应变均呈线性变化,且环向应变小于纵向应变发展;随着荷载的不断增大,试件整体应力增大,其中型钢的应变发展较快;当加载至极限荷载的50%时,型钢应变增长迅速,长柱中部发生轻微弯曲变形,试件刚度降低,应力-应变关系呈非线性变化;当荷载增大至峰值荷载的80%时,型钢刚度退化明显,再生混凝土开始膨胀,GFRP管的环向应变远大于其纵向应变,约束作用加强,且试件中部出现侧向挠曲变形,弯曲明显;当荷载达到峰值荷载时,型钢屈服,试件的压应变继续增大直至混凝土达到其极限压应变,导致再生混凝土被压碎,GFRP管破裂,试件弯曲严重,最终因试件失稳导致试件整体发生破坏,承载力迅速下降。
4.2 轴压荷载-位移曲线的比较
图10为GFRP管型钢再生混凝土组合柱轴压荷载-位移试验曲线与计算曲线的对比。由图10可知,计算曲线与试验曲线吻合较好,表明本文建立的有限元模型可以较好地模拟GFRP管型钢再生混凝土组合柱轴压性能非线性分析。另外,试件大致可以分为以下受力过程:加载初期,试件各组成部分都处于弹性阶段,曲线呈线性变化;随着轴压荷载的不断增加,型钢进入弹塑性阶段,承载力增长较为缓慢,变形明显增大,曲线较为平缓;达到峰值荷载后,GFRP管撕裂破坏,试件承载力急剧下降,反映了GFRP管的脆性特点。
4.3 轴压承载力计算值与试验值的比较
表5为GFRP管型钢再生混凝土组合柱承载力试验值与有限元计算值的比较。
组合柱承载力试验值与有限元计算值对比 表5
试件编号 |
试验值Nue/kN | 计算值Nuc/kN | 相对误差δ/% |
GSRC-1 |
3 387.75 | 3 500.12 | 3.2% |
GSRC-2 |
3 301.99 | 3 400.54 | 2.9% |
GSRC-3 |
3 270.54 | 3 309.45 | 1.1% |
GSRC-4 |
3 010.04 | 3 201.27 | 6.0% |
GSRC-5 |
3 056.58 | 3 197.45 | 4.4% |
GSRC-6 |
3 343.74 | 3 460.92 | 3.4% |
GSRC-7 |
3 489.74 | 3 591.46 | 2.8% |
GSRC-8 |
2 871.67 | 3 023.56 | 5.0% |
GSRC-9 |
2 764.80 | 2 856.71 | 3.2% |
GSRC-10 |
3 183.89 | 3 322.94 | 4.2% |
GSRC-11 |
3 242.19 | 3 410.56 | 4.9% |
注:δ=(Nuc-Nue)/Nue, Nue为承载力试验值,Nuc为承载力计算值。
由表5可知,组合柱承载力试验值与有限元计算值吻合相对较好,最大误差为6%左右,满足计算要求。另外,试件的有限元计算值普遍比试验值高,原因可能为:1)在ABAQUS软件模拟中试件假定为理想状态,试件各材料都是各向匀质的,而实际试件在制作过程中很难保证材料的各向匀质性; 2)有限元模拟中默认试件处于理想铰支状态,在试验中,试件并非理想的铰支约束;3)在试件加载过程中,人工调整试件难免产生相对误差,使得试件具有一定的初始偏心距,而在有限元模拟中,试件是处于理想轴心受荷状态。
4.4 有限元参数扩展分析
在上述有限元分析的基础上,以GSRC-5试件(短柱)和GSRC-9试件(长柱)为标准对GFRP管型钢再生混凝土组合柱进行有限元参数扩展分析,表6为组合柱有限元扩展参数分析极限承载力,图11为不同参数下组合柱的荷载-位移曲线。
由有限元扩展参数计算结果可知,扩展参数对组合柱轴压性能的影响规律:1)由表6和图11(a)可知, GSRC-5-1与 GSRC-5-2试件以GSRC-5试件为参考,采用不同的径厚比进行参数分析。得出随着径厚比的增大,组合柱的轴压承载力逐渐降低。因为随着径厚比的增大使得在加载后期GFRP管对再生混凝土的约束能力降低,导致组合柱承载力下降,扩展试件GSRC-5-2承载力比试验试件GSRC-5降低了7%左右;GSRC-9-1与GSRC-9-2试件以GSRC-9试件为参考,采用不同的径厚比进行参数分析。试验试件GSRC-9承载力比扩展试件GSRC-9-2提高了12%左右,径厚比对长柱承载力的影响较大。 2)由表6和图11(b)可以看出,GSRC-5-3试件与GSRC-5-4试件以GSRC-5试件为参考,采用不同的型钢强度进行参数分析。型钢强度对长短柱轴压性能的影响相似,加载初期,试件的初始刚度相似;随着荷载的增大,试件刚度开始出现差异,随着型钢强度的提高,试件刚度随之提高,组合柱的承载力也随之增大。由表6可知,对于短柱,型钢强度等级为Q390的GSRC-5-4试件较型钢强度等级为Q235的GSRC-5试件承载力提高了7%;对于长柱,GSRC-9-3试件与GSRC-9-4试件以GSRC-9试件为参考,采用不同的型钢强度进行参数分析。型钢强度等级为Q390的GSRC-9-4试件比型钢强度等级为Q235的GSRC-9试件承载力提高了7%左右。
有限元参数扩展分析计算结果 表6
试件编号 |
径厚比D/t | 型钢强度/MPa | 轴压承载力N/kN |
GSRC-5 |
20 | Q235 | 3 197.45 |
GSRC-5-1 |
25 | Q235 | 3 093.86 |
GSRC-5-2 |
40 | Q235 | 2 982.01 |
GSRC-5-3 |
20 | Q345 | 3 290.27 |
GSRC-5-4 |
20 | Q390 | 3 425.12 |
GSRC-9 |
20 | Q235 | 2 856.71 |
GSRC-9-1 |
25 | Q235 | 2 706.08 |
GSRC-9-2 |
40 | Q235 | 2 548.65 |
GSRC-9-3 |
20 | Q345 | 2 921.56 |
GSRC-9-4 |
20 | Q390 | 3 055.26 |
5 结论
通过试验和有限元分析对GFRP管型钢再生混凝土组合柱的轴压性能进行了研究,主要得到以下结论:
(1)组合柱的破坏过程为内部的型钢最先屈服,加载至峰值荷载,再生混凝土被压碎,GFRP管对核心再生混凝土的约束力也逐渐增大直至破裂。短柱主要发生材料强度破坏,中长柱主要是非弹性弯曲失稳破坏。
(2)组合柱承载力随再生骨料取代率、组合柱长细比的增大而降低,最大降幅分别为9.8%,9.6%;随型钢配钢率与再生混凝土强度等级的提高而增大,最大增幅分别为12.4%,5.7%。合理地选择设计参数对组合柱轴压承载力是有利的。
(3)通过有限元软件ABAQUS对GFRP管型钢再生混凝土组合柱轴压性能的模拟,获得的模拟计算结果与试验结果吻合较好,表明建立的有限元模型能够较好地应用于GFRP管型钢再生混凝土柱轴压性能数值分析。
(4)组合柱的轴向承载力随着径厚比的增大而减小,其变形能力也随着降低,径厚比对组合柱的承载力与变形能力影响不明显;提高型钢的强度,组合柱的轴压承载力也随之提高,其变形能力也随之提高。实际工程中在保证承载力相等的情况下,通过适量减小径厚比来减小组合柱截面面积,降低自重,达到节约成本的效果。
[2] 李秋义,全洪珠,秦原.混凝土再生骨料[M].北京:中国建筑工业出版社,2011.
[3] 杜朝华.建筑垃圾再生骨料混凝土及构件受力性能研究[D].郑州:郑州大学,2012.
[4] SHI CAIJUN,LI YAKE,ZHANG JIAKE,et al.Performance enhancement of recycled concrete aggregate-a review[J].Journal of Cleaner Production,2016,112(1):466-472.
[5] MOHANRAJ E,KANDASAMY S,MALATHY R.Behavior of steel tubular stub and slender columns filled with concrete using recycled aggregates[J].Journal of the South African Institution of Civil Engineering,2011,53(2):31-38.
[6] AJDUKIEWICA A,KLISZCZEWICZ A.Influence of recycled aggregate on mechanical properties of HS/HPC[J].Magazine of Concrete Composites,2002,24 (2):269-279.
[7] 梁炯丰,张广武,胡明华,等.型钢部分包裹再生混凝土短柱轴压性能试验研究[J].建筑结构,2019,49(8):72-75,42.
[8] 马辉,郭婷婷,李哲,等.圆钢管型钢再生混凝土组合柱轴压性能及承载力计算[J].建筑结构学报,2018,39(2):97-105.
[9] 牛娜娜.再生粗骨料对混凝土抗压性能的影响研究[J].施工技术,2019,48(S1):469-471.
[10] 肖建庄,黄一杰.GFRP管约束再生混凝土柱抗震性能与损伤评价[J].土木工程学报,2012,45(11):112-116.
[11] 肖建庄,杨洁.玻璃纤维增强塑料约束再生混凝土轴压试验[J].同济大学学报(自然科学版),2009,37(12):1586-1591.
[12] 欧阳煌,黄奕辉,钱在兹,等.玻璃纤维(GFRP)片材约束混凝土的受力性能分析[J].土木工程学报,2004,37(3):26-34.
[13] 王连广,秦国鹏,周乐.GFRP管钢骨高强混凝土组合柱轴心受压试验研究[J].工程力学,2009,26(9):170-175.
[14] 薛建阳,崔卫光,陈宗平,等.型钢再生混凝土组合柱轴压性能试验研究[J].建筑结构,2013,43 (7):73-76.
[15] 马辉,毛肇玮,刘云贺,等.圆钢管型钢再生混凝土长柱轴压下受力性能的分析[J].西北农林科技大学学报(自然科学版),2017,45(9):145-154.
[16] 胡东晋,王振山,刘云贺,等.玻璃纤维复合材料构件抗压性能研究[J].西安理工大学学报,2015,31(2):1-8.
[17] 刘威.钢管混凝土局部受压时的工作机理研究[D].福州:福州大学,2005.
[18] 宋志刚,樊成,宋力.GFRP管混凝土轴压短柱承载力研究[J].水利与建筑工程学报,2017,15(2):71-75.