1、第 49 卷 第 5 期 西安建筑科技大学学报（自然科学版） Vol.49 No.52017 年 10 月 J.Xian Univ. of Arch.ultimate strength analysis;slab-column connections;Mohr-Coulomb Criterion钢筋混凝土板柱结构是一种只由楼板和柱组成的结构体系，这种结构有施工简便、空间利用率高、布局灵活等诸多优点，在停车场、住宅、仓库等建筑中有着广泛的应用。但是板柱结构的节点容易发生脆性的冲切破坏，而且一个节点的冲切破坏可能引起相邻节点因过载而相继破坏，进而引发结构的连续倒塌。 板柱节点抗冲切性能研究已有近百
2、年的历史，并积累了相当丰富的试验数据 1。特别是对中柱节点在轴向荷载下受力性能的研究较为充分。现有的研究表明板柱节点在不平衡弯矩的作用下其抗冲切承载力会有所下降，而板柱结构在承受竖向荷载的同时由于重力荷载的偏心，横向风作用以及地震作用不可避免的会产生不平衡弯矩，进而降低节点的承载力。不平衡弯矩下的冲切破坏机理较为复杂，相关的研究也相对较少，我国混凝土结构设计规范也仅在附录中列出了计算公式 2 。国内外学者提出了一些板柱节点在不平衡弯矩作用下冲切承载力的计算方法或计算模型，如：1976 年 Islam 和 Park 在屈服线理论的基础上采用拟梁法建立了节点的承载力的计算公式 3, 4；1981年
3、 Dilger 等人在低周反复荷载试验的基础上运用屈服线理论得到了不平衡弯矩作用下的荷载计算公式5；我国的马元昌、吕西林于 2001 年在抗弯塑性铰理论的基础上进行极限分析得到了抗弯和抗剪相关方程的联合计算公式 6； 2007 年 Ying Tian 基于试验研究和历史数据库提出以拟梁法为基础的计算模型 7； 2013 年 洪小滔研究了混凝土强度、冲垮比、配筋率对偏载下板柱节点冲切强度的影响 8；2014年 Choi 等人基于临界截面应变分布提出考虑不平衡弯矩作用下的计算模型 9。在板柱节点抗冲切承载力的塑性极限分析方面，西 安 建 筑 科 技 大 学 学 报（自然科学版）也取得相关的研究成果
4、。1976 年 Braestrup 等利用修正的 Mohr-Coulomb 理论求解出了轴对称冲切问题的最小上限解 10；1986 年蒋大骅等采用抛物线形的 Mohr-Coulomb 屈服准则，使得圆柱节点的最小上限解的表达式更为简洁 11。1993 年周朝阳在采用塑性分析方法研究板柱节点冲切承载力时提出考虑斜截面面内的旋转的错动-转动复合机构 12，相比于传统的错动模型能考虑径向和环向配筋的影响；2002 年，Salim W, Sebastian W M 在蒋大骅的基础上对极限分析做了一些修正 13，2017 年易伟建、刘彪提出基于塑性理论的钢筋混凝土板冲切开裂滑移模型 14。1 混凝土的屈
5、服准则对于轴对称的板柱节点塑性极限分析，Braestrup 等利用修正的 Mohr-Coulomb 破坏准则 10，蒋大骅等采用抛物线型的 Mohr-Coulomb 破坏准则11，严宗达等采用双剪强度理论 15，魏雪英等采用统一强度理论 16。抛物线型 Mohr-Coulomb 破坏准则能较好的反映混凝土在较大正应力作用下剪力并非与正应力成比例增加的现象，并且连续可导，便于运算。本文采用文献 11 给出的抛物线形Mohr-Coulomb 准则：ntn修正的M o h r - C o u l o m b公式 1 cftfnt,图 1 混凝土的破坏准则Fig.1 Failure criterion
6、 for concrete* 2ntn t 1fKfMERGEFORMAT (1)* MERGEFORMAT ccttfvm(2)* 1214KmMERGEFORMAT (3)其中：f t和 fc分别为混凝土的有效抗拉和抗压强度，ft 和 fc 分别为混凝土的抗拉和抗压强度， t 和 c 为相应的有效系数。根据图 1 可以求出抛物线外法线与 n 轴夹角 与应力之间的关系为* ntdtaMERGEFORMAT (4)结合式（1）可以得到：* tnta2fKMERGEFORMAT (5)则应力与夹角的关系为* 2nt tt(1co)fMERGEFORMAT (6)将上式应力代入到应力圆中，可以求得
7、其主应力* 21t3t=cs1KfMERGEFORMAT (7)以及切向正应力：* 2t t=1csKfMERGEFORMAT (8)2 运动机构及对应的应变与应力对板柱节点的冲切破坏进行塑性极限分析时，将柱头与未屈服区域的板视为刚体，外力作用下柱头（冲切破坏锥）与板之间形成屈服面，达到塑性极限状态。当不存在不平衡弯矩时，柱头整体可能的运动机构只会发生向下平动。当存在不平衡弯矩时，与轴向力和不平衡弯矩相对应，柱头整体在平动的同时还会发生转动，按图 2 所示的运动机构，可计算轴力和不平衡弯矩所做的功。将轴向力和弯矩共同作用转化为一个偏心作用力，节点在偏载作用下发生冲切破坏。此机构分为三部分，不动
8、刚体，运动的柱头刚体以及塑性屈服区域，屈服域的厚度为 。运动的柱头刚体绕着点 R（x 0,z0）产生角速度为 的旋转运动。柱的中心轴线的运动可分解为一个沿 z 轴的平动-x0，和一个角速度为 的旋转运动，将这种运动机构称为平动-转动复合运动机构。现假设破坏锥体屈服域为如图 2 所示坐标中的一待求函数，可以在直角坐标表示为* (,)zgxyMERGEFORMAT (9)在本文中 、 分别表示函数 g 对 x 一阶和x二阶偏导。 CCR(x0,z0)Pe0=Munb/Pg(x,y)v u vuzxh图 2 运动机构Fig.2 Collapse Mechanism由不平衡弯矩产生的偏心距为* unb
9、0MePMERGEFORMAT (10)在图 2 所示的平动-转动复合运动机构下，锥体屈服域上任意一点的竖直方向的速率为 u，水平方向的速率为 v，正方向如图 2 所示，则可求得运动速率：* 0()uxgzMERGEFORMAT (11)屈服面上沿 Y 轴方向的运动速率为零，现假设屈服域的形状函数能满足：* 0,xgMERGEFORMAT (12)则将运动速率分解为向屈服域的法线方向和切面上可以求得应变速率：* n2nt n21xxyuvgguvMERGEFORMAT (13)在塑性增量理论的假设为：在塑性变形过程中的任一微小时间增量内，瞬时应力偏量与应变增量成比例，即：* /nntdsMER
10、GEFORMAT (14)式中，d 为一与路径有关的非负标量。应力偏量：* n13=/2sKMERGEFORMAT (15)结合式* MERGEFORMAT (6)、* MERGEFORMAT * MERGEFORMAT (14)可以得：* nt/cotMERGEFORMAT (16)由式* MERGEFORMAT (14)、* MERGEFORMAT (15)得：* MERGEFORMAT n2dK(17)可得 0，现将式* MERGEFORMAT (16)代入式* MERGEFORMAT (6) 屈服面上的应力可以表示为* 2ntnttn(1)KffMERGEFORMAT (18)3 虚功
11、方程在屈服域上存在应力 n n 应变速率 、n做功，所以单位面积的屈服面内力做虚功为nt* intn2tn() 1-WuvKfMERGEFORMAT (19)屈服域的面积为* 2d1xyxsgdMERGEFORMAT (20)屈服域做的总内虚功：* MERGEFORMAT IidsW西 安 建 筑 科 技 大 学 学 报（自然科学版）(21)联立式* MERGEFORMAT (13)、* MERGEFORMAT (20)、* MERGEFORMAT (21)可以求得屈服域内力做的总内虚功为 2 2I t=11xyx xyDuvWuvgKgfd * MERGEFORMAT (22)外力做的虚功为
12、* E0PeMERGEFORMAT (23)根据外力和内力虚功相等的原则，则有：* MERGEFORMAT EI(24)所以可以求得承载力： 0022 0002,()1d1xy xc xDxyPgxyxgzKfKyme* MERGEFORMAT (25)4 承载力计算4.1 最小上限解求解承载力 P 为一泛函，现在确定函 P 的边界条件，即屈服面在 xy 平面的投影，积分区域 D 的边界条件 由边界 1 和边界 2 两部分组成，对于板的冲切问题上边界 1 一般认为是沿着柱边破坏的，即柱子的形状边界函数，是固定的；而对于板的下边界 2 未知的，即边界 2 为可变边界。在特定机构下的最小承载问题转
13、化为对含有一个可变边界条件泛函 P 的极值问题，对于这类含有多元函数的泛函的可进行变分求解。现令：00 2220 ()1xc xyxgzKfFmeKg * MERGEFORMAT (26)根据文献 17泛函取得极值的必要条件是要求满足欧拉方程：* 0 ,xygggFFxDMERGEFORMAT (27)同时，补充在可变边界条件上，得： 2222, 0coscosxyxyggFygznnn * MERGEFORMAT (28)对于一般的情况，要直接求出解析解存在一定的困难。现考虑没有不平衡弯矩作用时的特殊情况。当偏心距为零时柱的轴线不发生转动，只有平动运动，得：* 0=xzMERGEFORMAT
14、 (29)对于这种特殊的情况，式* MERGEFORMAT (25)简化为 2,1dxycxyDfPgKgm* MERGEFORMAT (30)同时式* MERGEFORMAT (26)简化为* 21cxyfFgMERGEFORMAT (31)代入欧拉方程* MERGEFORMAT (27) 得：* 0xyMERGEFORMAT (32)得出的结论与文献 18一致，即文献 18所求出的结论只是本文的一种特殊情况。若对于圆柱板轴对称问题可以对调和方程做一个柱坐标变换，同时代入边界条件 1、 2 则可以求出圆柱板的冲切解析解。现在假设板上表面沿着柱边破坏，下表面破坏边缘距离柱中心的距离为 d1/2
15、，屈服面 g 函数的边界条件 1、 2 为* 11/,0,ghMERGEFORMAT (33)求得：* 022c11lnefdKPmdMERGEFORMAT (34)当：* 01/eMERGEFORMAT (35)可以取得最小值，计算结果与文献 11 一致。4.2 最小上限解近似求解对于最小上限解，当求解析解存在一定困难时，可以采用类似于构造形函数的方法求近似解，当构造的曲面函数合适时求出的最小解与解析解通常误差会很小。所以本文拟采用构造形函数的方法近似求解，对计算结果进行分析得出实用的结论。为便于分析对上文坐标转换为柱坐标体系，根据坐标变换公式：* cosinxryMERGEFORMAT (
16、36)得到：* sicosindrrxyggrMERGEFORMAT (37)所以求得承载力： 22c01drrDfBPAKgmex * MERGEFORMAT (38)上式中： 0022cos+()sincos(39)rrgzBrx图 3 偏载下典型的冲切破坏面 8Fig.3 Typical failure surface under eccentric load文献 8一共进行了 8 块板在偏载作用下的冲切试验，根据实验现象表明破坏曲面在远离偏载侧的冲切锥一般要略为平缓一些，同时曲面下凸，对于r 是单调，所以可以根据这些实验现象构造屈服面函数。xy柱板边界支持 冲切底面边界 rd 1/2c
17、osb0.5d1+b 0.5d1-bPdh g(r, )zx图 4 破坏示意图Fig.4 Diagram failure surface如 图 4现在拟定屈服面 g 函数的边界条件1、 2 为* 1/2,0cos,dbhMERGEFORMAT (40)边界条件 1 沿柱边，为固定边界条件；边界条件 2 中 d1，b 为待定量，为可变边界条件。现假设曲面函数 g 沿在同一角度 处为满足边界条件1、 2 二次多项式的曲线。根据这些假设可以构造曲面函数：211,coscosrdardhabb * MERGEFORMAT (41)式* MERGEFORMAT (41) 中 a 为二次多项式待定参数，当
18、 a=0 时，破坏面退化成圆台。将式* MERGEFORMAT (41)代入* MERGEFORMAT (38)即可求得承载力 P，为求得所有可以发生的运动机构中的 P 的最小值，可以设计计算方法寻找 P 的最小值。现在未定量为边界条件变量 d1、b，曲面二次多项式参数 a，和旋转中心 R0(x0,z0)四个待定量。所以本文在确定相关参数的取值范围后采用多重搜索的方法求出 P 的最小值。根据试验现象一般可以大致确定一些参数的取值范围。边界条件中， d+h d1 d+10h，0b h；二次多项式参数：0a0.5；旋转中心参数： x00。西 安 建 筑 科 技 大 学 学 报（自然科学版）现取 m
19、=10013, 18，通过数值计算的方法可以求出式* MERGEFORMAT (38)的在偏载下的承载力，通过多重搜索可以找出承载力的最小值，其最小承载力与非偏载下的理论值式* MERGEFORMAT (34)中 Pe0=0 的比值如图 5 中的点标所示。可见承载力主要与参数 d/h，e 0/(d+h)有关。为应用方便，对计算的结果采用参数 d/h、e 0/(d+h)的多项式进行拟合，图 5 中的线标为式* MERGEFORMAT (42) ，可见拟合的公式与直接计算的结果是比较接近的。 01/31/301.5.425ePPeddhh* MERGEFORMAT (42)0.51.0.52.0.
20、53.0.54.0.20.4.60.81. P/e=e/(d+h)d/h=2 46值/ dh=图 5 计算与简化公式对比Fig.5 Comparison of calculation and simplified formula5 本文与其它计算方法对比与评价在工程应用中板的厚度一般采用有效厚度 h0，为保持一致，本文的计算方法中 h 视为板的有效厚度，e 0=Munb/Pe0 为偏心距。在计算方柱时采用等效圆柱直径，根据面积等效的原则， d=1.13c。轴对称条件下的承载力计算采用文献 18的公式：* 01/61/230cm0.84ehcfuhMERGEFORMAT (43)式中，u m 同
21、为距离柱边 h0 的周长，根据式* MERGEFORMAT (42)可以得到偏心作用下的冲切承载力。我国规范 2对于非预应力板在不配置箍筋或弯起钢筋时在集中反力作用下或者局部荷载的冲切承载力 Fl 计算应该满足式* MERGEFORMAT (44)，相关参数的含义详见规范 2。* htm0.7lfuMERGEFORMAT (44)在偏心作用下，采用等效集中反力设计：* 0unb, m0c+ABleqlMFhIMERGEFORMAT (45)对于方柱截面边长为 c 中间板柱节点，代入相关公式，令其 e0=Munb/Fl， 可以计算得：* 0,1.2lqlehMERGEFORMAT (46)式*
22、MERGEFORMAT (46) 的等效集中荷载 应,leqF满足式* MERGEFORMAT (44)，得到在偏载作用下的承载力计算公式：* 0htm0.7=12elfuFecMERGEFORMAT (47)文献 19对规范公式进行了修正，给出了中柱节点在 Vg 和不平衡弯矩 Munb 联合作用下的半经验半理论承载力公式，对于柱尺寸为 c 的方柱，其计算化简后如下* g001.2eFchMERGEFORMAT (48)* 0tm0fuMERGEFORMAT (49)* yty=1.4ffMERGEFORMAT (50)式中：F 0 为板柱节点在没有不平衡弯矩作用下的冲切承载力，e 0=Mun
23、b/Vg 为等效偏心距，f t 为是混凝土轴心抗拉强度，f y 是纵向钢筋屈服强度值 c 是方柱的边长，h 0 为板的效高度，u m 为截面临界周长。所以在偏载作用下的承载力为* 0g0/(1.2)eFchMERGEFORMAT (51)现在对文献 19公式、我国现行规范考虑的偏心作用与本文的极限分析计算的结果进行对比如图 6 所示，文献 19公式和规范对中柱方形柱的偏心作用下的折减效应是一致的。由图 6 可知，本文的极限分析计算的偏心对承载力的影响的折减效应比规范和文献 19公式偏小，这可能是由于在偏载作用下所假设的曲面与理论的曲面存在偏差导致承载力偏高或者是由于混凝土并非理想的弹塑性材料，
24、在复杂的偏心受力状态下承载力相对于规范或其他基于试验统计得出的经验公式的折减效应偏小。在等效偏心距的坐标下，规范和文献 19公式对于中间方柱板都没有考虑柱尺寸和板厚度的比值对载力折减效应的影响，而本文的极限分析结果表明柱尺寸和板厚度的比值对折减效应有着较大的影响。 0.51.02.53.04.204.681.0 P/e=e0/(c+h) 本 文 c/d=2本 文 4本 文 6规 范 及 文 献 19图 6 偏心效应对承载力的折减效应对比Fig.6 Comparison of the effect of eccentricity on bearing capacity文献 20整理了一共 39
25、个试验数据，适用于本文公式计算的有效数据一共 36 个，对其试验值承载力和计算承载力的比值在 c/h0 的坐标下所示，图中虚线为回归线。0123456.5.12.0 试 验 值 c/h0GB01-2.973值计 算 值 0123456.5.12. 试 验 值 c/h0180.49计 算 值0123456.5.1. 值c/h0.7869计 算 值图 7 试验与计算对比Fig.7 Comparison between experiment and calculation由图 7 可知，规范 GB50010-2010 和文献 19公式对于 c/h0 的影响趋势都不太准确，而本文的极限分析的结果对 c
26、/h0 的影响的趋势预测相对更为准确，但极限分析计算得出的节点的承载力整体偏高，基于试验数据对本文的极限分析的结果进行调整，调整后的计算公式如下：* 01/31/30001.5.8ePPeddhdhhMERGEFORMAT (52)调整后的计算值与试验的对比如图 8：0123456.5.12.0 值c/h0值1.024计 算 值图 8 调整后试验与计算对比Fig.8 Comparison of test and calculationafter adjustment此时平均值为 1.028，变异系数为 0.144，离散性以及对柱宽与板厚比值影响趋势都明显优于规范GB50010-2010 和文献
27、 19 。当采用附加剪力表示弯矩对冲切承载力的不利影响时，将板柱节点的等效剪力设计值表示为1/31/300, 0.5.8lleqlFeddhdhh* MERGEFORMAT (53)式中 Fl,eg 为等效剪力设计值，应小于非偏心冲切下的承载力的值；F l 为重力荷载设计值，e0=Munb/Fl 为偏心距； h0 为板的有效厚度； d 为柱的直径，为方柱时取为 1.13 倍的方柱边长。此计算公式适用于存在不平衡弯矩的中柱板柱节点冲切承载力计算。6 结论(1)采用塑性极限分析方法，针对圆柱形板柱节点在轴向力和弯矩共同作用的受力条件，提出平动-旋转复合运动机构，通过推演得到板柱节点抗冲切承载力的上
28、限解，物理含义明确。(2)求出了求解上限解的方法，为简化计算，通过构造的破坏曲面，采用数值算法求出了节点在偏载作用下的承载力，对计算结果进行了简化，并与规范和其他的计算方法进行了比较，试验数据表明西 安 建 筑 科 技 大 学 学 报（自然科学版）平动-旋转复合运动机构模型相比与偏心剪应力模型能更准确的反映柱子与板厚的比值对承载力的趋势影响。(3)根据试验数据对计算公式进行了修正调整，调整后的计算公式与试验结果吻合良好，明显的优于规范及其它计算公式，可用于实际工程设计。(4)由于收集的试验数据有限，方柱板边长与板厚度比值的影响还有待于更进一步的研究。参考文献 References 1 易伟建.
29、 混凝土结构试验与理论研究 M. 北京: 科学出版社, 2012:286-291.YI Weijian. Experimental and theoretical research on concrete structure M. Beijing: Science Press, 2012:286-291. 2 中华人民共和国建设部.混凝土结构设计规范:GB50010-2010 S . 北京: 中国建筑工业出版社, 2010.Ministry of Construction of the PRC. Code for Design of Concrete Structure: GB50010-20
30、10S, Beijing: China Architecture Industry Press, 2010.3 ISLAM S, PARK R. Tests on slab-column connections with shear and unbalanced flexure J. Journal of the Structural Division, 1976, 102(3): 549-568.4 PARK R, ISLAM S. Strength of slab-column connections with shear and unbalanced flexure J. Journal
31、 of the Structural Division, 1976, 102(1): 549-568.5 DILGER W H, GHALI A. Shear reinforcement for concrete slabs J. Journal of the Structural Division, 1981, 107(12): 2403-2420.6 马云昌, 吕西林. 钢筋混凝土板柱节点的抗震性能研究 J. 建筑结构学报, 2001, 22(4): 49-54.MA Yunchang, LV Xilin. Study on seismic behavior of reinforced c
32、oncrete slab joints. J. Journal of Building Structures, 2001, 22(4): 49-54.7 TIAN Y. Behavior and modeling of reinforced concrete slab-column connections M. AUSTIN:THE UNIVERSITY OF TEXAS AT AUSTIN, 2007: 146-162.8 洪小滔. 偏载下钢筋混凝土板柱节点冲切破坏试验研究 D; 长沙:湖南大学, 2013: 20-56.HONGXiaotao. Experimental Study on
33、Punching Shear Failure ofReinforced Concrete Slab-Column Connections with Eccentric Loading D; Changsha:Hunan University, 2013: 20-56.9 CHOI K K, SHIN D W, PARK H-G. Shear-Strength Model for Slab-Column Connections Subjected to Unbalanced Moment J. ACI Structural Journal, 2014, 111(3): 491-152.10 BR
34、ASTRUP M W, NIELSEN M P, JENSEN B C, et al. Axisymmetric punching of plain and reinforced concrete M. Structural Research Laboratory, Technical University of Denmark, 1976: 1-12.11 JIANG D H, SHEN J H. Strength of concrete slabs in punching shear J. Journal of Structural Engineering, 1986, 112(12):
35、2578-2591.12 周朝阳. 塑性理论在混合材料复杂结构复合受力极限分析中的应用/第二届全国结构工程学术会议论文集C.北京: 清华大学出版社, 1993, 338-343.ZHOU Zhaoyang. Application of plasticity theory to limit analysis of composite structure under complex loading; A collection of papers on the second National Academic Conference on Structural Engineering C.Chang
36、sha: Tsinghua University press 1993: 338-343. 13 SALIM W, SEBASTIAN W. Plasticity model for predicting punching shear strengths of reinforced concrete slabs J. Structural Journal, 2002, 99(6): 827-35.14 易伟建, 刘彪. 基于塑性理论的钢筋混凝土板冲切开裂滑移模型J. 西安建筑科技大学学报(自然科学版), 2017, 49(1): 42 47.YI Weijian , LIU Biao. Cra
37、ck slip model for punching capacity of reinforced concrete slab based on the plastic theoryJ. J. Xian Univ. of Arch. & Tech. (Natural Science Edition), 2017, 49(1): 4247. 15 严宗达. 用双剪强度理论解混凝土板冲切的轴对称问题 J. 工程力学, 1996, 13(1): 1-7.Yan Zongda. Solution of the axisymmetrical punching problem of concrete sl
38、ab by the twin shear strength theory J. Engineering Mechanics , 1996, 13(1): 1-7. 16 魏雪英, 杨政, 于澍, 等. 用统一强度理论求解轴对称混凝土板的冲切强度 J. 工程力学 , 2002, 15(9): 92-96.WEI Xueying, YAN Zheng, YU Shu. With the unified strength theory to solve the axisymmetric punching strength of the concrete slab J. Engineering Mec
39、hanics, 2002,19( 05): 92-96. 17 老大中. 变分法基础 M. 国防工业出版社, 北京：2007: 163-167.LAO Dazhong. The basis of the variational method M. Beijing：Defense industry press,2007: 163-167.18 易伟建,刘彪.混凝土板柱节点冲切承载力的极限分析J.工程力学,2017,34(8):125-132，143.YI Weijian , LIU Biao. ultimate strength analysis of punching shear capaci
40、ty of concrete slab - column connectionS. J. Engineering Mechanics,2017,34(8):125-132,143.19 舒兆发. 板柱结构的中柱节点传递剪力和不平衡弯矩的冲切承载力 J. 湖南大学学报（自然科学版） , 2003, 30(5): 55-60.SHU Zhaofa. Punching shear capacity of mid column joints of plate column structure under transfer shear and unbalanced bending momentsJ. Journal of hunan university（Natural Science Edition）, 2003, 30(5): 55-60.20 朱泽华. 偏载下板柱结构中柱节点受力性能试验研究 D. 长沙:湖南大学, 2015: 56-68.ZHU Zehua Experimental investigations on eccentrically loaded RC interior slab-column connections D.Changsha:Hunan University, 2015:56-68.
