钢筋混凝土板柱结构是一种只由楼板和柱组成的结构体系,这种结构有施工简便、空间利用率高、布局灵活等诸多优点,在停车场、住宅、仓库等建筑中有着广泛的应用.但是,板柱结构的节点容易发生脆性冲切破坏,而且一个节点的冲切破坏可能引起相邻节点因过载而相继破坏,进而引发结构连续倒塌.

对板柱节点抗冲切性能研究已有近百年的历史,并积累了相当丰富的试验数据^[1].特别是中柱节点在轴向荷载下受力性能的研究较为充分.现有的研究表明:板柱节点在不平衡弯矩地作用下其抗冲切承载力会有所下降,而板柱结构在承受竖向荷载的同时由于重力荷载的偏心,横向风作用以及地震作用不可避免地会产生不平衡弯矩,进而降低节点的承载力.不平衡弯矩下的冲切破坏机理较为复杂,相关的研究也相对较少,我国混凝土结构设计规范也仅在附录中列出了计算公式^[2].

国内外学者提出了一些板柱节点在不平衡弯矩作用下冲切承载力的计算方法或计算模型,如:1976年Islam和Park在屈服线理论地基础上采用拟梁法建立了节点的承载力的计算公式^[3-4]; 1981年Dilger等人在低周反复荷载试验的基础上运用屈服线理论得到了不平衡弯矩作用下的荷载计算公式^[5]; 我国的马元昌、吕西林于2001年在抗弯塑性铰理论的基础上进行极限分析得到了抗弯和抗剪相关方程的联合计算公式^[6]; 2007年Ying Tian基于试验研究和历史数据库提出以拟梁法为基础的计算模型^[7]; 2013年洪小滔研究了混凝土强度、冲垮比、配筋率对偏载下板柱节点冲切强度的影响^[8]; 2014年Choi等人基于临界截面应变分布提出考虑不平衡弯矩作用下的计算模型^[9].

在板柱节点抗冲切承载力的塑性极限分析方面,也取得相关的研究成果.1976年Braestrup等利用修正的Mohr-Coulomb理论求解出了轴对称冲切问题的最小上限解^[10]; 1986年蒋大骅等采用抛物线形的Mohr-Coulomb屈服准则,使得圆柱节点的最小上限解的表达式更为简洁^[11].1993年周朝阳在采用塑性分析方法研究板柱节点冲切承载力时提出考虑斜截面面内的旋转的错动-转动复合机构^[12],相比于传统的错动模型能考虑径向和环向配筋的影响; 2002年,Salim W, Sebastian W M在蒋大骅的基础上对极限分析做了一些修正^[13],2017年易伟建、刘彪提出基于塑性理论的钢筋混凝土板冲切开裂滑移模型^[14].

对于轴对称的板柱节点塑性极限分析,Braestrup等利用修正的Mohr-Coulomb破坏准则^[10],蒋大骅等采用抛物线型的Mohr-Coulomb破坏准则^[11],严宗达等采用双剪强度理论^[15],魏雪英等采用统一强度理论^[16].抛物线型Mohr-Coulomb破坏准则能较好的反映混凝土在较大正应力作用下剪力并非与正应力成比例增加的现象,并且连续可导,便于运算.本文采用文献[11] 给出的抛物线形Mohr-Coulomb准则:

图1 混凝土的破坏准则
Fig.1 Failure criterion of the concrete

(σ_n)/(f'_t)+1/K((τ_nt)/(2f'_t))²=1(1)

m=(f'_c)/(f'_t)=(v_vf_c)/(v_tf_t)(2)

K=1/4[m+2(1-(m+1)^1/2)](3)

其中:f'_t和f'_c分别为混凝土的有效抗拉和抗压强度,f_t和f_c分别为混凝土的抗拉和抗压强度,v_t和v_c为相应的有效系数.

根据图1可以求出抛物线外法线与τ_n轴夹角α与应力之间的关系为

tanα=-(dτ_nt)/(dσ_n)(4)

结合式(1)可以得到:

tanα=2K(f'_t)/(τ_nt)(5)

则应力与夹角的关系为

(6)

将上式应力代入到应力圆中,可以求得其主应力.

(7)

以及切向正应力:

σ_t=[1-K(cscα+3)²]f'_t(8)

对板柱节点的冲切破坏进行塑性极限分析时,将柱头与未屈服区域的板视为刚体,外力作用下柱头(冲切破坏锥)与板之间形成屈服面,达到塑性极限状态.当没有不平衡弯矩时,柱头整体可能的运动机构只会发生向下平动.当存在不平衡弯矩时,与轴向力和不平衡弯矩相对应,柱头整体在平动的同时还会发生转动,按图2所示的运动机构,可计算轴力和不平衡弯矩所做的功.

将轴向力和弯矩共同作用转化为一个偏心作用力,节点在偏载作用下发生冲切破坏.此机构分为三部分,不动刚体Ⅰ,运动的柱头刚体Ⅱ以及塑性屈服区域Ⅲ,屈服域的厚度为δ.运动的柱头刚体Ⅱ绕着点R(x₀,z₀)产生角速度为ω的旋转运动.柱的中心轴线的运动可分解为一个沿z轴的平动-ωx₀,和一个角速度为ω的旋转运动,将这种运动机构称为平动-转动复合运动机构.现假设破坏锥体屈服域为如图2所示坐标中的一待求函数,可以在直角坐标表示为

z=g(x,y)(9)

在本文中g'_x、g″_xx分别表示函数g对x一阶和二阶偏导.

图2 运动机构
Fig.2 Collapse mechanism由不平衡弯矩产生的偏心距为

在图2所示的平动-转动复合运动机构下,锥体屈服域上任意一点的竖直方向的速率为u,水平方向的速率为v,正方向如图2所示,则可求得运动速率:

(11)

屈服面上沿Y轴方向的运动速率为零,现假设屈服域的形状函数能满足:

(12)

则将运动速率分解为向屈服域的法线方向和切面上可以求得应变速率:

(13)

塑性增量理论的假设为:在塑性变形过程中的任一微小时间增量内,瞬时应力偏量与应变增量成比例,即:

(14)

式中,dλ为一与路径有关的非负标量.

应力偏量:

s_n=σ_n-(σ₁+σ₃)/2=2K(15)

结合式(6)、(14)可以得:

(16)

由式(14)、(15)得:

(17)

可得现将式(16)代入式(6)屈服面上的应力可以表示为

(18)

在屈服域上存在应力σ_n τ_n应变速率做功,所以单位面积的屈服面内力做虚功为

(19)

屈服域的面积为

ds=(1+g^'2_x+g^'2_y)^1/2dxdy(20)

屈服域做的总内虚功:

W_I=∫_ds∫W_ids(21)

联立式(13)、(20)、(21)可以求得屈服域内力做的总内虚功为

W_I=∫_Dxy∫[(u+vg'_x)(1-K)+K(u²+v²)/(u+vg'_x)(1+g^'2_x+g^'2_y)]f'_tdxdy(22)

外力做的虚功为

W_E=P(e₀-x₀)ω(23)

根据外力和内力虚功相等的原则,则有:

W_E=W_I(24)

所以可以求得承载力:

∫_Dxy∫(f'_c)/(m(e₀-x₀))[[(x-x₀)+(g-z₀)g'_x](1-K)+K((x-x₀)²+(g-z₀)²)/((x-x₀)+(g-z₀)g'_x)(1+g^'2_x+g^'2_y)]dxdy(25)

在工程应用中板的厚度一般采用有效厚度h₀,为保持一致,本文的计算方法中h视为板的有效厚度,e₀=M_unb/P_e0为偏心距.在计算方柱时采用等效圆柱直径,根据面积等效的原则, d=1.13c.轴对称条件下的承载力计算采用文献[18]的公式:

P_e0=0=0.84(h₀/c)^1/6(100ρ)^1/2f^1/3_cu_mh₀(43)

式中,u_m同为距离柱边h₀的周长,根据式(42)可以得到偏心作用下的冲切承载力.

我国规范^[2]对于非预应力板在不配置箍筋或弯起钢筋时在集中反力作用下或者局部荷载的冲切承载力F_l计算应该满足式(44),相关参数的含义详见规范^[2].

F_l≤0.7β_hf_tηu_mh₀(44)

在偏心作用下,采用等效集中反力设计:

F_l,eq=F_l+(α₀M_unbα_AB)/(I_c)u_mh₀(45)

对于方柱截面边长为c中间板柱节点,代入相关公式,令其e₀=M_unb/F_l,可以计算得:

F_l,eq=F_l(1+1.2(e₀)/(c+h₀))(46)

式(46)的等效集中荷载F_l,eq应满足式(44),得到在偏载作用下的承载力计算公式:

F_e0=F_l=(0.7β_hf_tηu_mh₀)/(1+1.2(e₀)/(c+h₀))(47)

文献[19]对规范公式进行了修正,给出了中柱节点在V_g和不平衡弯矩M_unb联合作用下的半经验半理论承载力公式,对于柱尺寸为c的方柱,其计算化简后如下

(V_g)/(F₀)=1/(1+1.2(e₀)/(c+h₀))(48)

F₀=ξf_tu_mh₀(49)

ξ=1.4ρf_y/(f_f+ρf_y)(50)

式中:F₀为板柱节点在没有不平衡弯矩作用下的冲切承载力,e₀=M_unb/V_g为等效偏心距,f_t为是混凝土轴心抗拉强度,f_y是纵向钢筋屈服强度值 c是方柱的边长,h₀为板的效高度,u_m为截面临界周长.

所以在偏载作用下的承载力为

V_g=F₀/(1+1.2(e₀)/(c+h₀))(51)

现在对文献[19]公式、我国现行规范考虑的偏心作用与本文的极限分析计算的结果进行对比如图6所示,文献[19]公式和规范对中柱方形柱的偏心作用下的折减效应是一致的.

由图6可知,本文的极限分析计算的偏心对承载力的影响的折减效应比规范和文献[19]公式偏小,这可能是由于在偏载作用下所假设的曲面与理论的曲面存在偏差导致承载力偏高或者是由于混凝土并非理想的弹塑性材料,在复杂的偏心受力状态下承载力相对于规范或其他基于试验统计得出的经验公式的折减效应偏小.在等效偏心距的坐标下,规范和文献[19]公式对于中间方柱板都没有考虑柱尺寸和板厚度的比值对载力折减效应的影响,而本文的极限分析结果表明柱尺寸和板厚度的比值对折减效应有着较大的影响.

图6 偏心效应对承载力的折减效应对比
Fig.6 Comparison of the effect of eccentricity on bearing capacity

文献[20]整理了一共39个试验数据,适用于本文公式计算的有效数据一共36个,对其试验值承载力和计算承载力的比值在c/h₀的坐标下所示,图中虚线为回归线.

图7 试验与计算对比
Fig.7 Comparison between experiment and calculation

由图7可知,规范GB50010-2010和文献[19]公式对于c/h₀的影响趋势都不太准确,而本文的极限分析的结果对c/h₀的影响的趋势预测相对更为准确,但极限分析计算得出的节点的承载力整体偏高,基于试验数据对本文的极限分析的结果进行调整,调整后的计算公式如下:

P=(P_e0=0)/(1+(e₀)/(d+h₀)[0.15((e₀)/(d+h₀))^1/3+0.8](d/(h₀))^1/3)(52)

调整后的计算值与试验的对比如图8.

图8 调整后试验与计算对比
Fig.8 Comparison of test and calculation after adjustment

此时平均值为1.028,变异系数为0.144,离散性以及对柱宽与板厚比值影响趋势都明显优于规范GB50010-2010和文献[19].

当采用附加剪力表示弯矩对冲切承载力的不利影响时,将板柱节点的等效剪力设计值表示为

F_l,eq=F_l+(F_le₀)/(d+h₀)[0.15((e₀)/(d+h₀))^1/3+0.8](d/(h₀))^1/3(53)

式中:F_l,eg为等效剪力设计值,应小于非偏心冲切下的承载力的值; F_l为重力荷载设计值,e₀=M_unb/F_l为偏心距; h₀为板的有效厚度; d为柱的直径,为方柱时取为1.13倍的方柱边长.此计算公式适用于存在不平衡弯矩的中柱板柱节点冲切承载力计算.

(1)采用塑性极限分析方法,针对圆柱形板柱节点在轴向力和弯矩共同作用的受力条件,提出平动-旋转复合运动机构,通过推演得到板柱节点抗冲切承载力的上限解,物理含义明确.

(2)求出了求解上限解的方法,为简化计算,通过构造的破坏曲面,采用数值算法求出了节点在偏载作用下的承载力,对计算结果进行了简化,并与规范和其他的计算方法进行了比较,试验数据表明平动-旋转复合运动机构模型相比与偏心剪应力模型能更准确的反映柱子与板厚的比值对承载力的趋势影响.

(3)根据试验数据对计算公式进行了修正调整,调整后的计算公式与试验结果吻合良好,明显的优于规范及其它计算公式,可用于实际工程设计.

(4)由于收集的试验数据有限,方柱板边长与板厚度比值的影响还有待于更进一步的研究.