摘要
利用ABAQUS接口User Element(UEL)子程序,建立了分离式自计算单元长度的界面黏结单元BS1,该单元能按指定的黏结滑移关系工作。基于简化的Eligehausen黏结滑移本构,对叠合梁-现浇柱子结构试验开展了数值模拟,结果与试验吻合良好。同时,验证了单元自计算长度的准确性,对比时间效率提高了25%。基于上述验证,建立30个叠合梁-现浇柱子结构模型,采用叠合梁现浇层厚度、T型翼缘板宽度、力作用方向等参数进行分析。结果表明,采用90°弯钩锚固和锚固板锚固的2种子结构,承载性能相差不大,延性相差约12%,宜优先选择90°弯钩锚固形式;承载力分析表明我国规范中对于叠合梁的T型翼缘板的宽度规定较国外规范更加可靠;地震作用下,力作用方向的改变,子结构前后承载力相差12%,变形能力相差55%,这种差异值得深入讨论。研究结果为考虑黏结滑移的装配式混凝土结构的数值分析提供了可靠参考。
关键词
在钢筋混凝土(RC)结构中,钢筋与混凝土之间的黏结滑移效应对构件或结构的受力性能、抗震性能、抗连续倒塌、抗大变形能力的影响较为突
随着PC框架结构的广泛应用,装配式混凝土结构的有限元数值模拟研究方法受到国内外学者关注。Zhao
PC框架结构中,钢筋的锚固及搭接是联系梁柱等构件的关键,影响构件及节点的受力性能。目前,钢筋搭接的研究较
目前,PC结构设计采用等同现浇处理,但与RC结构相比,PC结构因根本构造而引起结构性能的差异仍需重视,可建立精细化有限元数值模拟,考虑黏结滑移、不同搭接及锚固形式以及叠合面的处理等因素,使得PC结构有限元分析更加精确。
分析PC框架梁柱子结构性能是研究PC框架结构整体性能的基础,为了解PC框架子结构在构造参数等影响下的结构响应规律,文中基于ABAQUS中UEL(User Element)子程序,建立自计算界面长度黏结单元,模拟钢筋混凝土黏结滑移效应。对现浇柱-叠合梁子结构试件采用精细化建模,发现计算结果与试验结果吻合较好,验证了该单元的高效性,分析多项构造参数敏感性,对比各锚固形式对子结构的影响,为PC框架子结构及节点精细化模拟提供参考。
利用ABAQUS中UEL接口子程序,自主开发了一种新型自计算长度分离式空间黏结界面单元,单元基于Fortran语言编程中加入了计算空间节点距离相应语句,实现单元长度的自行计算,可减少计算迭代,提升计算效率,如

图1 混凝土、钢筋与界面单元建模示意图
Fig. 1 Modeling details of concrete、steel bar and interface element
BS1单元变形计算迭代是基于Newton-Raphson法建立,如
, | (1) |
, | (2) |
式中:Δδn为位移增量;kn为切向刚度;ΔPn为节点力增量;δn为位移总量。

图2 Newton-Raphson法
Fig. 2 Newton-Raphson method
单元刚度矩阵确定其变形准则,先根据4个节点位移差计算得节点位移,建立切向力所作虚功W和节点外力所作虚功T,根据D'Alembert's原理建立关系求解出轴向刚度Ku 。
。 | (3) |
。 | (4) |
。 | (5) |
参考Keuser所提方

图3 单元刚度矩阵建立
Fig. 3 Establishing stiffness matrix of the element
界面单元4个节点分别与混凝土单元(C3D8R)和钢筋单元(B31)的节点Tie在一起,通过黏结界面单元变形来表征钢筋与混凝土之间滑移量。同时,四节点黏结滑移单元采用非线性位移方程,模拟精度要高于常规使用的弹
BS1单元变形与混凝土和钢筋2种材料的弹塑性参数有关。其中,混凝土材料,选用ABAUQS中提供的塑性损伤模型(Concrete damaged plasticity model, CDP模型),观察对比混凝土塑性发展,本构选自我国规范《混凝土结构设计规范》(GB50010)附录

图4 混凝土本构模型
Fig. 4 Concrete constitutive model

图5 钢筋本构模型
Fig. 5 Steel constitutive model
综合考虑现有黏结滑移本构关系,为更好体现滑移峰值阶段稳定性,选用Eligehausen

图6 黏结-滑移本构模型
Fig. 6 Bond-slip constitutive model
为优化编程过程和保证模拟的收敛性及精确性,参考Low
节点 | τ/MPa | s/mm | 节点 | τ/MPa | s/mm |
---|---|---|---|---|---|
线弹性终点A | 峰值段终点C | 0.003 | |||
峰值段始点B | 残余点D | 0.01 |
注: d为钢筋直径;H为变形钢筋肋高;s为变形钢筋肋净距;为混凝土轴心抗压强度。
地震或冲击倾覆等灾害中,因复杂边界条件及刚度分布等问题,混凝土框架结构中损伤更易集中在底层边跨处,参考文献[

图7 试件PC-0-B几何尺寸及配筋(单位:mm)
Fig. 7 Dimension and details of specimen PC-0-B(unit:mm)
材料 | fy/MPa | fu/MPa | A/% | fcu,m/MPa |
---|---|---|---|---|
预制梁 | — | — | — | 37.4 |
现浇部分 | — | — | — | 34.2 |
A8(梁、柱) | 390 | 561 | 28.3 | — |
C20(柱) | 476 | 602 | 23.0 | — |
C20(梁) | 482 | 625 | 21.0 | — |
C18(梁) | 503 | 622 | 23.3 | — |
注: fy为钢筋屈服强度;fu为钢筋极限抗拉强度;A为钢筋断后伸长率;fcu,m为混凝土抗压强度平均值。

图8 试验装置图
Fig. 8 Test setup
为验证BS1单元在PC框架结构有限元模拟中的可行性,根据PC-0-B有限元模型。对子结构中的柱构件,分析试验中边界条件并进行合理简化处理,以便建模分析,混凝土采用C3D8R实体单元,CDP模型塑性参数如
泊松比 | 偏心率 | 膨胀角 | fb0/fc0 | K | 黏性系数 |
---|---|---|---|---|---|
0.2 | 0.1 | 30° | 1.16 | 0.667 | 0.000 5 |
对于水平叠合面和竖向键槽结合面的处理,文中采用“接触对”进行设置,通过定义“接触对”的法向作用和切向作用来实现新旧混凝土交界面按指定条件工作。对于法向作用,采用“硬接触”,接触面只传递正向法向应力,无法传递拉应力,可设置分离与否;对于切向作用,采用各向同性库伦摩擦模型来定义,其中,摩擦系数参考美国ACI规

图9 PC-0-B有限元模型
Fig. 9 Finite element model of PC-0-B
计算结果如

图10 子结构混凝土塑性损伤及钢筋mise应力分布
Fig. 10 Mise stress distribution of bars and damage cloud of composite beams-cast-in-place column subassemblages

图11 试验试件破坏模式
Fig. 11 Experimental failure mode of specimen

图12 试验和数值模拟荷载-位移曲线
Fig. 12 Load-displacement of test and numerical simulation
基于上述分析,为验证黏结单元自计算长度对分析结果的影响,在Fortran语言编程中隐藏单元自计算长度语句,选取单元长度为5 mm、10 mm和25 mm 3种,对比使用单元自计算长度时12.5 mm情况,以观察不同黏结单元长度下,PC-0-B模型计算收敛性和精确性,并验证自计算单元长度的可靠性。

图13 不同界面单元尺寸的荷载-位移曲线
Fig. 13 Load-displacement of different interface element sizes
在PC-0-B模型基础上,设计90°弯钩锚固和锚固板锚固2种子结构分析模型,设立不考虑黏结滑移模型为对照组,3种模型的材性、构造和配筋与试验试件PC-0-B保持一致,锚固形式均满足规
依据规

图14 不同参数子结构的荷载-位移曲线
Fig. 14 Load-displacement under different parameters
模型 | Py/kN | Δy/mm | Pm/kN | Δm/mm | μ | Pm相对值 | μ相对值 |
---|---|---|---|---|---|---|---|
PC-N-CT60 | 122.69 | 29.67 | 140.73 | 180.00 | 6.07 | 0.98 | 1.25 |
PC-N-CT80 | 126.06 | 27.38 | 144.07 | 135.96 | 4.85 | 1.00(1.01) | 1.00(0.87) |
PC-N-CT100 | 123.29 | 23.86 | 137.78 | 140.90 | 5.91 | 0.95 | 1.22 |
PC-N-CT120 | 125.22 | 25.34 | 139.22 | 167.47 | 6.52 | 0.97 | 1.34 |
PC-H-CT60 | 123.18 | 25.47 | 140.13 | 200.90 | 7.79 | 0.99 | 1.38 |
PC-H-CT80 | 123.88 | 30.50 | 142.12 | 172.91 | 5.65 |
1.0 |
1.0 |
PC-H-CT100 | 119.63 | 23.97 | 133.71 | 167.11 | 6.97 | 0.94 | 1.23 |
PC-H-CT120 | 122.30 | 29.11 | 135.48 | 172.69 | 5.93 | 0.95 | 1.05 |
PC-A-CT60 | 127.75 | 25.66 | 141.65 | 169.82 | 6.54 | 0.99 | 1.49 |
PC-A-CT80 | 124.04 | 29.34 | 142.77 | 131.58 | 4.40 | 1.00(1.00) | 1.00(0.79) |
PC-A-CT100 | 119.72 | 22.02 | 133.48 | 160.62 | 7.29 | 0.93 | 1.66 |
PC-A-CT120 | 119.15 | 24.39 | 137.64 | 134.75 | 5.52 | 0.96 | 1.25 |
PC-N-TW990 | 159.54 | 28.73 | 177.48 | 153.43 | 5.34 | 0.98 | 1.05 |
PC-N-TW1160 | 165.55 | 31.88 | 180.55 | 161.56 | 5.07 | 1.00(0.94) | 1.00(1.19) |
PC-N-TW1280 | 163.22 | 28.70 | 180.87 | 149.72 | 5.22 | 1.00 | 1.03 |
PC-N-TW1350 | 164.28 | 28.51 | 181.10 | 142.77 | 5.01 | 1.00 | 0.99 |
PC-H-TW990 | 170.09 | 48.87 | 196.48 | 204.79 | 4.06 | 1.02 | 0.95 |
PC-H-TW1160 | 172.61 | 42.14 | 193.03 | 183.96 | 4.26 |
1.0 |
1.0 |
PC-H-TW1280 | 173.68 | 39.18 | 183.88 | 186.23 | 4.61 | 0.95 | 1.08 |
PC-H-TW1350 | 175.10 | 40.13 | 194.09 | 179.95 | 4.48 | 1.00 | 1.05 |
PC-A-TW990 | 170.02 | 36.17 | 188.44 | 169.13 | 4.67 | 0.96 | 1.15 |
PC-A-TW1160 | 172.41 | 44.52 | 195.84 | 180.90 | 4.06 | 1.00(1.01) | 1.00(0.95) |
PC-A-TW1280 | 173.16 | 36.62 | 192.51 | 174.23 | 4.76 | 0.98 | 1.17 |
PC-A-TW1350 | 173.72 | 38.90 | 193.88 | 174.72 | 4.49 | 0.99 | 1.11 |
PC-N-FU | 98.21 | 24.39 | 121.29 | 64.54 | 2.48 | 0.84 | 0.51 |
PC-N-FD | 126.06 | 27.38 | 144.07 | 135.96 | 4.85 | 1.00(1.01) | 1.00(0.86) |
PC-H-FU | 97.02 | 26.26 | 120.71 | 72.75 | 2.56 | 0.85 | 0.45 |
PC-H-FD | 123.88 | 30.50 | 142.12 | 172.91 | 5.65 |
1.0 |
1.0 |
PC-A-FU | 97.70 | 24.46 | 121.46 | 64.58 | 2.54 | 0.85 | 0.58 |
PC-A-FD | 124.04 | 29.34 | 142.77 | 131.58 | 4.40 | 1.00(1.00) | 1.00(0.78) |
注: Py为屈服荷载;Δy为屈服位移;Pm为峰值荷载;Δm为屈服位移;μ为延性系数。μ相对值中*代表参数整体大组(锚固形式)的基准对照,括号代表大组内不同锚固形式的小组(各参数)对照,例如:第二行Pm相对值中1.00(1.01)中1.00为本小组(不考虑粘结滑移时)基准,1.01=144.07/142.12。
考虑楼板效应对子结构的影响,参考中、美、加、欧等国家及地区相关规
地震或者往复荷载作用下,PC框架结构中受力方向发生改变,其结构性能存在一定差异。设置向上(FU)和向下(FD)2种加载方向,探讨2种受力状态下子结构承载性能差异。模型编号定义方式:PC-H-FU为采用90°弯钩锚固,加载方向向上。
文中利用ABAQUS中UEL端口开发了新型自划分网格黏结滑移单元BS1,将单元应用至装配式混凝土现浇柱-叠合梁子结构的有限元模拟中,结果表现良好,并开展了系列参数分析,主要结论如下:
1)单元长度自计算功能有利于高效可靠计算,时间成本降低了25%,该单元可应用至装配式混凝土结构有限元精细化模拟计算分析。
2)经可靠验证后,利用BS1单元建立不考虑黏结滑移、90°弯钩锚固和锚固板锚固3种子结构模型,分析了现浇层厚度、翼缘板影响宽度、力作用方向等参数影响;发现梁端90°弯钩锚固设置更能有效保证子结构的整体延性,建议装配式混凝土结构中钢筋锚固采用90°弯钩锚固形式。
3)对比中、美、加、欧等国家及地区相关规范,根据各规范计算翼缘板影响宽度并建立模型分析,结果表明,在翼缘板影响下,子结构承载性能提高了约26%,变形能力提高了约27%;对比发现,采用我国规范规定的翼缘板影响宽度计算分析时,子结构的延性偏低,承载力相当,整体比较适用。
4)在力作用方向改变下,子结构正反受力承载性能相差约12%,主要因为力方向的改变使结合面和叠合面处拉压应力的差距,其影响叠合梁的破坏损伤分布,从而对整体结构产生影响,结合面的构造需进一步考虑地震作用,加强界面设计。
参考文献
Lin H W, Zhao Y X, Ožbolt J, et al. The bond behavior between concrete and corroded steel bar under repeated loading[J]. Engineering Structures, 2017, 140: 390-405. [百度学术]
杨红, 高文生, 王志军. 空间框架简化为平面模型的抗震分析[J]. 重庆大学学报, 2008, 31(11): 1267-1272. [百度学术]
Yang H, Gao W S, Wang Z J. The rationality of plane model simplifying three-dimension frame in seismic response analysis[J]. Journal of Chongqing University, 2008, 31(11): 1267-1272. (in Chinese) [百度学术]
Sezen H, Setzler E J. Reinforcement slip in reinforced concrete columns[J]. ACI Structural Journal, 2008, 105(3): 280-289. [百度学术]
周新刚.混凝土后锚钢筋的静动力黏结特性与工程应用——第十三届全国结构工程学术会议特邀报告[J].工程力学,2004,21(S1):231-237. [百度学术]
Zhou X G. Static and dynamic bonding characteristics and engineering application of concrete back anchor reinforcement: Invited Report of 13th National Conference on Structural Engineering[J]. Engineering Mechanics,2004,21(S1):231-237. (in Chinese) [百度学术]
Santos J, Henriques A A. New finite element to model bond–slip with steel strain effect for the analysis of reinforced concrete structures[J]. Engineering Structures, 2015, 86: 72-83. [百度学术]
Alfarah B, Murcia-Delso J, López-Almansa F, et al. RC structures cyclic behavior simulation with a model integrating plasticity, damage, and bond-slip[J]. Earthquake Engineering & Structural Dynamics, 2018, 47(2): 460-478. [百度学术]
Zhao J, Sritharan S. Modeling of strain penetration effects in fiber-based analysis of reinforced concrete structures[J]. ACI structural journal, 2007, 104(2): 133. [百度学术]
蔡建国, 冯健, 王赞, 等. 预制预应力混凝土装配整体式框架抗震性能研究[J]. 中山大学学报(自然科学版), 2009, 48(2): 136-140. [百度学术]
Cai J G, Feng J, Wang Z, et al. Seismic behavior of frame comprised of precast prestressed concrete components[J]. Acta Scientiarum Naturalium Universitatis Sunyatseni, 2009, 48(2): 136-140. (in Chinese) [百度学术]
陈适才, 闫维明, 王文明, 等. 大型预制混凝土结构梁-柱-叠合板边节点抗震性能研究[J]. 建筑结构学报, 2011, 32(6): 60-67. [百度学术]
Chen S C, Yan W M, Wang W M, et al. Seismic behavior of exterior beam to column joints with layered slab of large precast concrete structures[J]. Journal of Building Structures, 2011, 32(6): 60-67. (in Chinese) [百度学术]
曹徐阳, 冯德成, 王谆, 等. 基于OpenSEES的装配式混凝土框架节点数值模拟方法研究[J]. 土木工程学报, 2019, 52(4): 13-27. [百度学术]
Cao X Y, Feng D C, Wang Z, et al. Investigation on modelling approaches for prefabricated concrete beam-to-column connections using openSEES[J]. China Civil Engineering Journal, 2019, 52(4): 13-27. (in Chinese) [百度学术]
Feng D C, Wu G, Lu Y. Finite element modelling approach for precast reinforced concrete beam-to-column connections under cyclic loading[J]. Engineering Structures, 2018, 174: 49-66. [百度学术]
Feng D C, Wu G, Lu Y. Numerical investigation on the progressive collapse behavior of precast reinforced concrete frame subassemblages[J]. Journal of Performance of Constructed Facilities, 2018, 32(3): 04018027. [百度学术]
Feng D C, Wang Z, Wu G. Progressive collapse performance analysis of precast reinforced concrete structures[J]. The Structural Design of Tall and Special Buildings, 2019, 28(5): e1588. [百度学术]
Sarraz A, Nakamura H, Kanakubo T, et al. Bond behavior simulation of deformed rebar in fiber-reinforced cementitious composites using three-dimensional meso-scale model[J]. Cement and Concrete Composites, 2022, 131: 104589. [百度学术]
Zhang N, Wu Y Z, Gu Q, et al. Refined three-dimensional simulation of ribbed bar pull-out tests based on an enhanced peridynamic model[J]. Engineering Structures, 2023, 278: 115519. [百度学术]
罗小勇, 龙昊, 曹琨鹏. 灌浆套筒连接装配式梁柱节点精细有限元模型[J]. 哈尔滨工程大学学报, 2021, 42(5): 641-648. [百度学术]
Luo X Y, Long H, Cao K P. High-precision finite element model of prefabricated reinforced concrete beam-column joint with grout sleeve[J]. Journal of Harbin Engineering University, 2021, 42(5): 641-648. (in Chinese) [百度学术]
高培楠, 邓扬, 李爱群, 等. 装配式混凝土框架结构连接节点抗震性能研究进展[J]. 工业建筑, 2021, 51(2): 171-185, 152. [百度学术]
Gao P N, Deng Y, Li A Q, et al. Review on seismic performance of connection joints of prefabricated concrete frame structures[J]. Industrial Construction, 2021, 51(2): 171-185, 152. (in Chinese) [百度学术]
刘璐, 黄小坤, 田春雨, 等. 配置大直径大间距HRB500高强钢筋的装配整体式钢筋混凝土框架节点抗震性能试验研究[J]. 建筑结构学报, 2016, 37(5): 247-254. [百度学术]
Liu L, Huang X K, Tian C Y, et al. Experimental study on seismic performance of precast RC frame joints with HRB500 high strength rebars of large diameter and spacing[J]. Journal of Building Structures, 2016, 37(5): 247-254. (in Chinese) [百度学术]
赵作周, 韩文龙, 钱稼茹, 等. 梁纵筋锚固板锚固装配整体式梁柱边节点抗震性能试验[J]. 哈尔滨工业大学学报, 2016, 48(12): 19-27. [百度学术]
Zhao Z Z, Han W L, Qian J R, et al. Seismic behavior of assembled monolithic beam-column exterior joints with beam longitudinal rebars anchored by heads[J]. Journal of Harbin Institute of Technology, 2016, 48(12): 19-27. (in Chinese) [百度学术]
陈适才, 闫维明, 李振宝, 等. 大型预制混凝土梁柱叠合板中节点整体抗震性能试验研究[J]. 工程力学, 2012, 29(2): 135-141. [百度学术]
Chen S C, Yan W M, Li Z B, et al. Experimental study on seismic behavior of layered slab and beam to column interior joints in precast concrete structures[J]. Engineering Mechanics, 2012, 29(2): 135-141. (in Chinese) [百度学术]
袁鑫杰, 李易, 陆新征, 等. 湿式连接装配式混凝土框架抗连续倒塌静力试验研究[J]. 土木工程学报, 2019, 52(12): 46-56. [百度学术]
Yuan X J, Li Y, Lu X Z, et al. Static progressive collapse test on prefabricated concrete frames with wet connections[J]. China Civil Engineering Journal, 2019, 52(12): 46-56. (in Chinese) [百度学术]
Kang S B, Tan K H, Yang E H. Progressive collapse resistance of precast beam–column sub-assemblages with engineered cementitious composites[J]. Engineering Structures, 2015, 98: 186-200. [百度学术]
Euser M K, Mehlhorn G, Cornelius V. Bond between prestressed steel and concrete: computer analysis using ADINA[J]. Computers & Structures, 1983, 17(5/6): 669-676. [百度学术]
中华人民共和国住房和城乡建设部. 混凝土结构设计规范: GB 50010—2002[S]. 北京: 中国建筑工业出版社, 2004. [百度学术]
Ministry of Construction of the People’s Republic of China. Code for design of concrete structures: GB 50010—2002[S]. Beijing: China Architecture & Building Press, 2004. (in Chinese) [百度学术]
Eligehausen R, Popov E P, Bertero V V. Local bond stress-slip relationships of deformed bars under generalized excitation: experimental results and analytical model[D]. Berkeley: California Univ, 1983: 162-180. [百度学术]
Lowes L N. Finite element modeling of reinforced concrete beam-column bridge connections[D]. Berkeley: California Univ, 1999, 235-256. [百度学术]
庞博. 装配整体式梁-柱子结构角柱失效后承载性能研究[D]. 长沙:湖南大学,2019. [百度学术]
Pang B. Research on bearing capacity of monolithic precast concrete beam-column subassemblages under corner column removal scenarios [D]. Changsha: Hunan University, 2019. (in Chinese) [百度学术]
中华人民共和国住房和城乡建设部. 装配式混凝土结构技术规程: JGJ 1—2014[S]. 北京: 中国建筑工业出版社, 2014. [百度学术]
Ministry of Housing and Urban-Rural Development of the People’s Republic of China. Technical specification for precast concrete structures: JGJ 1—2014[S]. Beijing: China Architecture & Building Press, 2014. (in Chinese) [百度学术]
ACI 318M-19. Building code requirements for structural concrete[S]. American Concrete Institude, Farminton Hills, MI, 2019. [百度学术]
NBCC. The national building code of Canada [S].Ottawa: Oational Research Council, 2010. [百度学术]
EN 1992-1-1:2004. Eurocode 2: Design of concrete structure -Part 1-1: General rules and rules for buildings[S]. Committee for Standardization, Brussels, 2004. [百度学术]