Processing math: 11%

基于分层壳单元的钢筋混凝土板面外冲切模拟研究

陶慕轩, 汪隆博, 赵继之

陶慕轩, 汪隆博, 赵继之. 基于分层壳单元的钢筋混凝土板面外冲切模拟研究[J]. 工程力学. DOI: 10.6052/j.issn.1000-4750.2024.08.0595
引用本文: 陶慕轩, 汪隆博, 赵继之. 基于分层壳单元的钢筋混凝土板面外冲切模拟研究[J]. 工程力学. DOI: 10.6052/j.issn.1000-4750.2024.08.0595
TAO Mu-xuan, WANG Long-bo, ZHAO Ji-zhi. NUMERICAL STUDY ON OUT-OF-PLANE PUNCHING OF REINFORCED CONCRETE SLAB UPON LAYERED SHELL MODEL[J]. Engineering Mechanics. DOI: 10.6052/j.issn.1000-4750.2024.08.0595
Citation: TAO Mu-xuan, WANG Long-bo, ZHAO Ji-zhi. NUMERICAL STUDY ON OUT-OF-PLANE PUNCHING OF REINFORCED CONCRETE SLAB UPON LAYERED SHELL MODEL[J]. Engineering Mechanics. DOI: 10.6052/j.issn.1000-4750.2024.08.0595

基于分层壳单元的钢筋混凝土板面外冲切模拟研究

基金项目: 国家自然科学基金项目(52178145),中央高校基本科研业务费(2024CDJYXTD-005)
详细信息
    作者简介:

    陶慕轩(1985−),男,上海人,副教授,主要从事钢-混凝土组合结构等方面研究(E-mail: taomuxuan@mail.tsinghua.edu.cn)

    汪隆博(2000−),男,黑龙江人,博士生,主要从事钢-混凝土组合结构研究(E-mail: 512303458@qq.com)

    通讯作者:

    赵继之(1996−),男,山东人,副教授,主要从事钢-混凝土组合结构、风电结构研究(E-mail: zhaojizhi@cqu.edu.cn)

  • 中图分类号: TU311

NUMERICAL STUDY ON OUT-OF-PLANE PUNCHING OF REINFORCED CONCRETE SLAB UPON LAYERED SHELL MODEL

  • 摘要:

    分层壳模型是实现结构精细化高效模拟的强有力工具。然而,面外冲切效应是传统分层壳模型唯一无法模拟的受力状态。为了进一步对冲切的全受力过程进行研究,该文基于传统的分层壳单元,在已有试验的基础上对面外的非线性本构进行了定义,并根据冲切试验数据库标定了面外本构的参数,通过二次开发实现了冲切现象的模拟,分析了混凝土强度、抗冲切钢筋等对冲切破坏形态以及冲切承载力的影响规律,形成了可以模拟冲切效应的新型分层壳单元,并通过冲切数据库对考虑面外冲切效应的分层壳单元的模拟效果进行验证。

    Abstract:

    The layered shell model is a potent tool for the precise and efficient simulation of structures. However, its out-of-plane punching effect remains the only loading condition that traditional layered shell models are incapable of simulating. To further investigate the complete force-bearing process of punching, this study, based on the traditional layered shell element, defines the nonlinear out-of-plane constitutive behavior on the basis of existing experimental data. The parameters of the out-of-plane constitutive behavior are calibrated according to a punching test database, enabling the simulation of punching phenomena through secondary development. The study analyzes the influence patterns of concrete strength and anti-punching reinforcements on the punching failure modes and punching bearing capacity, resulting in the formation of a novel layered shell element capable of simulating the out-of-plane punching effect. The simulation efficacy of the layered shell element considering the out-of-plane punching effect is validated through a punching database.

  • 钢筋混凝土板柱结构由于其结构形式简单、传力途径短捷、节约材料和降低工程成本等优越性,被广泛应用于水利、交通和工民建等建筑领域中[13]。但在集中或局部均布荷载与不平衡弯矩的共同作用下,荷载周围板内的混凝土处于典型的三维受力状态。一旦作用荷载超出材料的极限应力,板柱节点或荷载作用周围就会发生脆性冲切破坏。板的冲切破坏虽然只是一种局部破坏,但这种脆性破坏一旦发生容易引起连锁反应导致整体倒塌,造成巨大的经济损失。因此,板的冲切破坏已成为结构设计中着重考虑的问题[47]

    目前,学者对钢筋混凝土板的冲切破坏问题采用了试验研究、理论分析和数值模拟等方法进行研究,并提出了与实际破坏情况接近的理论模型以及相应的极限承载力计算方法[812]。当加入不同类型抗冲切钢筋之后,影响承载力的因素增加,板内的受力情况更加复杂,在此条件下冲切破坏机理的分析和极限承载力计算方法的给出更加困难,目前研究尚不充分。混凝土结构设计规范 (GB 50010−2010)[13]、美国规范(ACI318-14)[14]、欧洲规范(EN1992-1-1:2004)[15]和英国规范(BS 8110-97)[16]给出了配置抗冲切钢筋(箍筋和弯起钢筋)混凝土板的极限承载力计算方法,但大多计算公式为半理论半经验公式,且配置抗冲切钢筋混凝土板的计算公式与无抗冲切钢筋混凝土板的受冲切承载力计算公式还存在不连续等问题。

    因此,为了进一步对冲切的全受力过程进行研究,本文在已有试验的基础上,通过子程序的编写实现了冲切现象的模拟,分析了混凝土强度、抗冲切钢筋等对冲切破坏形态以及极限承载力的影响规律。首先,本研究在基于分层壳模型单元的基础上,考虑面外非线性冲切效应,对通用有限元程序MARC的分层壳截面模型接口GENSTR进行了二次开发,实现了针对冲切板的高效准确且建模方便的分层壳单元。本文主要介绍单元的原理和实现过程。然后,根据大量试验研究成果,对国内外普通配筋的冲切板试验进行归纳总结,分析提炼了冲切板的关键特性,详细讨论了面外剪力与剪切变形关系,并通过冲切试验数据库对关键参数进行标定。最后,本文将所开发的单元应用于构件分析,对国内外不同设计参数的冲切板试件进行模拟,用以验证该模型的准确性。冲切板的剪切承载力和开裂后剪切刚度及极限曲线等关键参数的计算,无论对于本模型的精度还是对于工程设计均具有重要意义。本文基于国内外不同设计参数的冲切板试验数据对这些关键参数的计算公式进行了比较和验证,为冲切板的非线性分析中关键环节提供一种可靠高效的计算工具。

    本文在面内分层壳单元开发的基础上[1719],考虑面内应力分量与面外分量的综合作用,对工程中常见的板单元进行受力分析,具体如下:

    壳单元的广义应变有8个分量,包括6个面内的分量和面外的两个剪切分量。d=[εx εy γxy φx φy φxy γxz γyz]T。各广义应变分量分别为x轴向应变,y轴向应变,xy方向的面内剪应变,x方向的曲率,y方向的曲率,面内剪切曲率,x平面(法方向x轴)的面外剪应变,y平面(法方向y轴)的面外剪应变。相对于广义应变,广义应力分量也有8个,即D=[Nx Ny Nxy Mx My Mxy Vxz Vyz]T。各广义应力分量分别为x方向的轴力,y方向的轴力,xy方向的剪力,x方向的弯矩,y方向的弯矩,面内的剪切弯矩,x平面(法方向x轴)的面外剪力,y平面(法方向y轴)的面外剪力。

    利用分层壳单元的广义应变,分层壳的应变推导如式(1)所示,式中zi代表第i层距离中和面的z向距离。

    {εix=εx+ziϕxεiy=εy+ziϕyγixy=γxy+ziϕxy (1)

    分层需要满足对应材料的本构模型,代入第i层的应变可以计算该层的应力,叠加后得到壳的广义应力和分层应力的关系。对于每层的轴力,需要通过该层的应力乘以厚度得到;对于每层的弯矩,需要在该层的应力乘以厚度的基础上,再乘以到中和轴的距离。具体的广义应力与单层应力的关系推导如式(2)所示。

    {Nx=σixtiNy=σiytiNxy=τixytiMx=σixtiziMy=σiytiziMxy=τixytiziVxz=GcγxzVyz=Gcγyz (2)

    对于面内分量,将分层的应变代入面内的本构模型中,得到对应的应力值以及该层的切线刚度。对于面外的剪力与剪应变关系,采用宏观本构,直接得到面外剪力与面外剪应变的关系。在分析过程中,认为在z向剪应变一致且与平面内的应变互相独立,对面外的力学关系进行单独定义。

    在MSC.MARC中,有面向对象的分层壳模型子程序GENSTR,可由用户自主实现分层过程的编写,材料本构的定义等过程,适用于面内弯曲、面内剪切、面外冲切等的分析。GENSTR子程序中关键分量的含义见表1

    表  1  GENSTR子程序参数变量表
    Table  1.  Parameters of GENSTR
    变量名称变量英文含义变量含义
    dgeneralized real stress strain law切线刚度矩阵
    etotatotal strain整体应变
    deincrement of strain整体应变增量
    sstress整体应力
    tstate variables状态参量
    dtincrement of state variables状态参量增量
    ngensnumber of generalized stress应力/应变分量数目
    muser element number用户定义单元数
    nnintegration point number积分点数目
    matus(1)user material number用户定义材料数目
    ht板的厚度
    下载: 导出CSV 
    | 显示表格

    在求解的过程中,关键在于得到分层的广义应力s与广义应变e之间的关系(即单元的广义本构关系),通过MSC.MARC的子程序GENSTR接口进行计算。建模示意图如图1所示,选择75号单元,并采用UFORMS施加荷载。中部节点通过刚性杆位移加载,用户自定义约束UFORMS实现第i个自由度在j个节点中(B1到Bj)的均匀分配,并满足式(3)和式(4)中的位移协调方程。通过设置用户自定义约束,实现施加的外部压力荷载在加载区域节点之间的均匀分配,同时又可进行位移控制加载,可以模拟真实的冲切加载情况[17, 18]

    uAi=jm=1Si,i+6(m1)uBmi=jm=1uBmi/j (3)
    {{\boldsymbol{S}}_{i,i}} = {{\boldsymbol{S}}_{i,i + 6}} = {{\boldsymbol{S}}_{i,i + 12}} = \cdots = {{\boldsymbol{S}}_{i,i + 6\left( {j - 1} \right)}} = 1/j (4)

    此种建模方法可以考虑面内应力分量与面外分量的综合作用,从而对冲切效应进行模拟[20]

    MSC.MARC每次调用GENSTR函数时,传入单元广义应变etota以及单元广义应变增量de。通过子程序内部的计算以及本构的定义,传出切线刚度矩阵d以及广义应变s。进而通过MSC.MARC内部牛顿迭代法进行内力求解。GENSTR中各应力应变分量定义如表2所示。其中,第1~第3变量和第6~第8变量量为面内分量,第4和第5变量为面外分量。本文分层壳单元开发的关键便是对面内和面外应力-应变的本构关系分别进行定义。

    图  1  冲切板建模示意图
    Figure  1.  Modeling of the punched slab
    表  2  GENSTR子程序应力应变对应关系表
    Table  2.  Stress-strain correspondence of GENSTR
    变量序号变形方向etota(广义应变)s(广义应力)
    1x轴向εxNx
    2y轴向εyNy
    3面内剪切γxyNxy
    4x面外剪切φxMx
    5y面外剪切φyMy
    6x向弯曲φxyMxy
    7y向弯曲γxzVxz
    8x轴向γyzVyz
    下载: 导出CSV 
    | 显示表格

    此外,如图2所示,在MSC.MARC中,二维壳单元分析有多处收敛准则需要进行设置,以保证较好的收敛性。此外,需要勾选矩阵的非正定计算以及每次迭代重新装配刚度矩阵。

    图  2  MSC.MARC求解准则设定
    Figure  2.  Solution criteria of MSC.MARC

    对于冲切效应的研究,面外非线性是重点,混凝土板也发生以冲切为主的破坏。因此,为了收敛性考虑,面内主要考虑未开裂的混凝土本构。在混凝土未发生开裂时,由于泊松效应的存在,两个方向的主应变无法进行解耦处理,即固定裂缝模型尚无法进行应用。根据弹塑性力学理论,在考虑二维本构时利用损伤因子进行处理可以方便地体现双轴之间的关系,即损伤模型进行混凝土本构的定义。在主拉应变达到开裂应变前,若采用考虑泊松效应的双轴线弹性本构,需按照式(5)计算。

    \left\{ \begin{aligned} & {\sigma _x^i = \frac{{E_{tk}'}}{{1 - {\mu ^2}}}( {\varepsilon _x^i + \mu \varepsilon _y^i} )} \\ & {\sigma _y^i = \frac{{E_{tk}'}}{{1 - {\mu ^2}}}( {\varepsilon _y^i + \mu \varepsilon _x^i} )} \\ & {\tau _{xy}^i = \frac{{E_{tk}'}}{{2( {1 + \mu } )}}\gamma _{xy}^i} \end{aligned}\right. (5)

    式中:E_{tk}' 为第i层的弹性模量;μ为泊松比;\varepsilon _x^i 为第i层的x方向轴向应变;\varepsilon _y^i 为第i层的y方向轴向应变;\gamma _{xy}^i 为第i层的剪切应变;\sigma _x^i 为第i层的x方向轴向应力;\sigma _y^i 为第i层的y方向轴向应力;\tau _{x}^i 为第i层的剪切应力。

    参考一维构件面外受剪参数标定方法[21, 22],板单元面外冲切骨架曲线应该考虑在发生冲切破坏后刚度发生退化的情况,承载力也随之急剧下降。板单元在承受面外冲切集中荷载时受力分为以下几个阶段:开始时面外未发生冲切破坏,刚度即为初始的面外剪切模量,面外冲切剪力随剪应变的增大随着上升;在剪切应变达到开裂点后,承载力增加明显缓慢;当面外剪切应变达到峰值时,承载力下降,截面出现破坏。因此,需要在本构关系中标定开裂参数,极限参数等有关临界点。根据上述论述,面外本构定义为三折线格式,如图3所示。

    图  3  面外三折线模型图
    Figure  3.  Out-of-plane trilinear model

    骨架曲线中开裂应变,极限应变需根据试验数据标定。根据前文所述,冲切承载力、冲切刚度、极限剪应变等与混凝土材料参数、板的参数、加载区域的形状面积等多方面因素有关,从板的构造条件到外部的荷载特性都对冲切承载力有较大影响,因此在冲切参数计算时需要对以下几点进行考虑:

    1)混凝土强度的影响。类似于梁的剪切破坏,冲切承载力计算时需要同时考虑混凝土受压强度与受拉强度的影响。

    2)配筋率的影响。板的面外冲切与梁的抗剪类似,可据此推测板的冲切与钢筋的相关关系。这些影响因素主要包括:粗骨料之间的咬合作用、裂面混凝土的剪切应力等。相对于梁,与配筋率有关的是钢筋的抗剪效果。根据研究,在板中钢筋部分的贡献占最终承载力的20%~30%。

    3)板的尺寸效应影响。在板的冲切效应中,尺寸效应为板的有效高度对冲切高度的影响。

    4)加载区域的影响。在大量的地震破坏分析后,柱顶加设柱帽后可以有效地减少冲切破坏的发生。这正是体现了加载区域的周长对于冲切承载力的影响。一般而言,加载区域的周长越小,加载圈的单元承受的平均剪应力更大,更易发生冲切破坏。此外,加载区域不同会导致冲切面的不同。中心加载与偏心加载对结构的破坏情况相差较大,需在计算中重点考虑。

    本文通过三折线模型进行板面外冲切模型的本构定义,需要标定板的面外冲切承载力、开裂荷载、初始刚度和开裂后的剪切刚度四个参数,从而确定三折线模型的具体形状。为了对四个重要参数进行拟合标定,本文建立了40个研究者共220组板受冲切破坏的试验数据库,其中159组试验包含加载点处的荷载-位移曲线。根据前文对冲切承载力影响参数的分析,对数据库中试件的参数进行整理,包括楼板的跨度尺寸、楼板的厚度、冲切荷载的加载方式和加载面的尺寸、边界条件、混凝土抗压强度指标、混凝土弹性模量、构件的配筋率和钢筋强度等。此外,对数据库中试件的荷载位移曲线进行提取用于本构模型的验证。对于有数据的试验曲线,可以直接对面外的荷载-位移曲线进行提取;对于只有曲线图的冲切数据,需要采用描点法对冲切数据进行标定。冲切数据库的建立方法如下:

    首先,选取数据库中典型的试验数据进行拟合标定以确定三折线本构模型中的参数,参数的取值范围如下:正方形楼板跨度a=1000 mm~2350 mm,圆形楼板跨度d=725 mm~1775 mm;楼板的厚度h=90 mm~550 mm;加载面为正方形;加载面的边长b=50 mm~300 mm;边界条件对于正方形楼板为四周固定,对于圆形楼板为圆板周围的支撑点固定;混凝土的抗压强度f_{\mathrm{c}}' =20 MPa~80 MPa;混凝土弹性模量为2.7×104 MPa~5.5×104 MPa;为方便钢筋部分的计算,对于分层壳模型,构件的配筋率重新定义为整个横截面的钢筋面积与混凝土面积之比,配筋率ρ=0.49%~2.83%。上述典型试验数据包含了常见的工程中发生楼板冲切破坏的参数范围,试验的基本信息及具体的试件尺寸和加载方式等如附表所示,变形单位为mm,力为kN,应力为MPa,R代表圆形试件。表中包含了冲切板数据库中最重要的参数,为下文参数标定以及模型验算提供了重要的参考。各试验的基本情况说明如下。

    1) Ozden试验[23]。Ozden在2006年报道了26组板的冲切破坏试验,试件为圆形试件,在试件的周围用12根拉杆进行支撑。试验中比较了普通和高强混凝土,钢筋配筋率,是否加入钢纤维,以及载荷的偏心加载对冲切试验结果的影响。试验结果表明,发生冲切破坏的临界失效裂缝总是位于距加载柱端部板的有效厚度的位置处。此外,混凝土的强度和配筋率对冲切有着关键性的作用。

    2) Kevin试验[24]。Kevin在2000年报道了6组板的冲切破坏试验,试件为正方形试件,在试件四周用圆钢管进行支撑。板材的厚度为135 mm~ 550 mm,宽度为925 mm~1950 mm。根据试验结果,楼板的跨度和厚度对于楼板的冲切承载力有很大影响,同时,加载区域的大小也对冲切效应影响显著。此外,当楼板的厚度大于200 mm时会具有比较明显的尺寸效应。

    3) Marzouk试验[25, 26]。Marzouk于1992年和1998年分别报道了17组和6组试验数据,均为正方形试件。通过对比混凝土强度,配筋率等影响因素,高强混凝土楼板的冲切破坏角在32~38度,在冲切破坏下的脆性更大,随着配筋率的增加,冲切承载力强度增加,但延性降低。

    4)其他学者试验。除了上述学者的试验,本文在进行不同的参数拟合时,还会分别选择GHANNOUM等[27]、QIAN等[28]、HASSAN等[29]、KAKUTA等[30]、INÁCIO等[31]、MORENO等[32]、NARAYANAN等[33]、OLIVEIRA等[34]、BANTHIA等[35]、NGUYEN-MINH等[36]、STEFANO等[37]、JING等[38]学者的试验数据,建立一套完整的冲切试验数据库,用于冲切本构的参数标定,并对冲切板的数值模拟进行验证。

    板的冲切承载力是板的面外冲切三折线本构模型的重要参数。虽然针对板的冲切破坏提出了很多模型,但板的冲切承载力依然很难完全用理论公式进行推导,国内外学者针对大量的试验数据拟合了冲切承载力的计算公式。目前,冲切承载力的计算有5种比较主流的设计规范为美国规范ACI 318−2014[14]、中国规范GB 50010−2010[13]、英国规范BS 8110-97[16]、欧洲规范EN 1992-1-1[15]和德国规范DIN 1045-1:2008[39]。为了将混凝土面外冲切的本构关系代入到分层壳有限元模拟中,主要参数与各国规范给出的拟合参数有两个不同,一是采用截面高度而并非有效高度,二是配筋率采用整个横截面的钢筋面积与混凝土面积之比。对公式中的结果进行适当的调整,使用附表中冲切数据库的试验数据进行拟合。楼板厚度和临界截面周长μ的单位用mm,最大荷载用单位N,冲切承载力的计算公式拟合结果如式(6)所示。

    {V_{\text{u}}} = 0.58\mu d\xi {\rho ^{1/3}}\sqrt {f_{\text{c}}'} (6)

    式中:f_{\text{c}}' 为混凝土抗压强度;ρ为配筋率;ξ=(400/d)1/4≥1.0,d为板的高度;μ为临界截面周长,对正方形加载板μ=4(b+3d),对于圆形加载板μ=4(c+3d);bc分别表示正方形加载面的边长和圆形加载面的直径。

    利用试验数据库对发生冲切破坏的最大承载力的拟合公式进行验证,并将5个常用规范的计算结果与实际结果相比较。以试验得到的承载力为横坐标,以计算结果与实际结果之比为纵坐标,得到散点图如图4所示。结果表明:利用本文的拟合公式得到结果的平均值为0.991、英国规范公式平均值为0.989、德国规范公式平均值为0.984,且标准方差也相对较小,利用这三个公式四个比较不错的结果。拟合公式(6)可以作为计算板的冲切承载力的公式依据。

    图  4  冲切承载力规范公式对比图
    Figure  4.  Comparison of punching bearing capacity obtained from various standard equations

    在板的面外计算中,另一个重要的参数是板的初始刚度,这也是三折线本构模型中初始阶段的关键参数。然而,部分收集到的试验数据缺少对该结果的测量。同时,受到试验中技术手段和观察上的限制,部分试验数据无法观测出明显的开裂前后试验构件刚度上的变化。图5(a)是SEVKET OZDEN[23]的NR1E0F0构件的面外加载曲线;图5(b)是H.MARZOUK[25]的NS2试验构件加载曲线。NR1E0F0构件的破坏曲线在达到最大冲切承载力之前具有明显的刚度变化情况,在加载达到开裂点后,构件的刚度明显变小。因此,在进行拟合时选择类似于NR1E0F0荷载-位移曲线的试验数据,具有明确的开裂点。虽然在拟合过程中并没有使用类似于NS2试件的没有明确刚度变化的荷载-位移曲线的试验数据,但有限元计算得到的荷载-位移曲线在初始刚度上会与试验结果略有不同,而整体的拟合情况仍然比较理想。使用附表中SEVKET OZDEN[23]和KEVIN[24]的试验数据,以及附表中的试验数据进行拟合。板的初始刚度主要受到混凝土的剪切模量和板的厚度影响,开裂之前钢筋参与受力可以忽略,因此,初始刚度与配筋率无关,拟合结果如式(7)所示。

    {k_{{\text{ini}}}} = 1.2Gd (7)

    式中,kini为板受冲切加载时的初始刚度;G是混凝土的剪切模量,d为板的高度。

    图  5  冲切试验破坏曲线图
    Figure  5.  Punching test failure curve diagram

    对拟合效果进行判断,以试验刚度结果为横坐标,用拟合结果与试验结果之比为纵坐标绘制散点图,如图6所示,得到的结果符合预期。

    图  6  初始刚度拟合结果对比图
    Figure  6.  Comparison of initial stiffness fitting results

    开裂荷载同样是三折线本构模型中的重要参数。在板加载到最大承载力的过程中,在开裂荷载的前后,板的刚度会发生明显的变化。为了得到开裂荷载的表达式,对附表中的数据进行拟合。通过数据分析,板的开裂荷载主要受到混凝土的强度和板的厚度影响。开裂荷载如式(8)所示。

    {V_{{\text{cr}}}} = 0.1d\sqrt {f_{\text{c}}'} (8)

    式中:f_{\text{c}}' 为混凝土圆柱体的强度;d为板的高度。

    拟合曲线对比如图7(a)所示,以试验结果为横坐标,拟合结果与试验结果之比为纵坐标绘制散点图,如图7(b)所示,可以看出拟合效果与试验结果之间的误差较小,拟合较为理想。

    图  7  开裂荷载拟合结果对比图
    Figure  7.  Comparison of cracking load fitting results

    板开裂后的剪切刚度直接影响了板的变形能力。在实际工程中,板达到最大承载力时的变形更被关注,而板的变形能力又受到开裂后剪切刚度的直接控制。因此,在进行拟合过程中,尽量保证板达到最大承载力时的变形有足够的拟合精度。在三折线模型中,开裂后剪切刚度的计算式为kcr=(VuVcr)/(δVc/kini),其中δ是板达到最大承载力时的变形。利用试验得到的开裂后剪切刚度结果进行kcr的拟合,得到的计算结果δ同时受到VuVcrkini的拟合精度的影响,可能存在较大的误差。在实际情况下,开裂荷载的大小和初始刚度的重要程度远不如构件的挠度变形。通过拟合开裂后剪切刚度以准确预估板达到最大承载力时的变形,并根据试验结果得到的最大承载力和变形与利用拟合结果得到的Vcrkini进行修正计算得到kcr,这样就避免了开裂荷载与初始刚度拟合误差对参数δ结果的影响。开裂后剪切刚度如式(9)所示。

    {k_{{\text{cr}}}} = \left( {{V_{\text{u}}} - V_{{\text{cr}}}'} \right)/\left( {\delta - V_{{\text{cr}}}'/k_{{\text{ini}}}'} \right) (9)

    式中:V_{{\text{cr}}}' k_{{\text{ini}}}' 为修正的开裂荷载与初始刚度;Vuδ为试验结果标定的最大承载力与对应变形。

    二维构件本构曲线中的下降段在模拟过程中会导致模型极难收敛,并且实际工程中更多关注的是冲切极限承载力的值,对于后续破坏过程的精细计算需求较小。出于收敛性考虑,采用−0.001倍的初始刚度定义下降段斜率[21, 22]

    使用附表中的试验数据计算得到开裂后剪切刚度,并与初始刚度相比,结果与板厚度与跨度有关,而随着配筋率增加,最大承载力提高,同时变形一般会较小,因此,开裂后刚度与配筋率相关,配筋率越大,开裂后刚度越大,从而得到开裂后刚度与开裂前刚度之比的公式如式(10)所示。

    \alpha = {k_{{\text{ini}}}}/{k_{{\text{cr}}}} = 1.5\left( {d{\text{/}}L + 5\rho } \right) - 0.12 (10)

    式中:α为开裂后与开裂前刚度之比;d为板的厚度;L为板的跨度

    拟合和试验结果散点图如图8所示,拟合结果成功模拟出了面外冲切开裂后的刚度变化。

    图  8  开裂后剪切刚度结果对比图
    Figure  8.  Comparison of post-cracking shear stiffness fitting results

    根据本节中对于三折线模型关键参数的拟合分析,可以得到四个重要参数的计算公式,并可以以此为标准确定三折线本构模型曲线,并提取成截面剪力-剪应变骨架曲线。曲线中的各个参数的计算公式汇总如式(11)~式(15)所示。在工程实践中,可以通过以下这些公式对冲切各个阶段的力学性能和关键指标进行初步计算。

    {V_{\text{u}}} = 0.58\mu d\xi {\rho ^{1/3}}\sqrt {f_{\text{c}}'} (11)
    {k_{{\text{ini}}}} = 1.2Gd (12)
    {V_{{\text{cr}}}} = 0.1d\sqrt {f_{\text{c}}'} (13)
    \alpha = {k_{{\text{ini}}}}{\text{/}}{k_{{\text{cr}}}} = 1.5\left( {d{\text{/}}L + 5\rho } \right) - 0.12 (14)
    {k_{{\text{deg}}}} = - 0.001{k_{{\text{ini}}}} (15)

    为了验证本文提出的三折线模型及其在MARC中二次开发得到的冲切计算程序的准确性,选取典型的试验数据及破坏现象进行模拟。根据前文的流程,在MSC.MARC中建立相应的有限元模型,利用程序的二次开发功能,将面外的三折线模型作为板的面外本构。在有限元模型中,对分析的100荷载步进行提取,得到每一步的应力云图。以M.Rizk(2010)的试验HS2数据为例[40]图9给出了分级荷载下横截面不同位置处的变形情况。冲切全过程的计算结果分为以下两阶段:

    1)面内弯曲变形。在冲切现象发生前,板的变形主要为弯曲变形,如试验结果图10中的200 kN、400 kN等级时的荷载情况,此时在冲切作用下集中位移的情况并不明显,整体变形剖面图基本与一维梁单元弯曲变形类似。由于面内采用的混凝土本构模型,弯曲变形时有可能会产生受拉斜裂缝,但未发生破坏。

    图  9  冲切破坏变形图
    Figure  9.  Punching failure deformation
    图  10  面内弯曲变形图
    Figure  10.  In-plane bending deformation

    2)面外冲切变形。当荷载加载到一定值时,中心集中荷载区域周边单元面外剪应变急剧增大,其他单元承荷水平基本不再增加,甚至出现卸载现象。最终,加载区域周边单元出现冲切破坏,结构发生失效,如图11所示。从试验结果图中可以看出,当荷载加载到600 kN处时,周围区域的变形便基本不再增加,而中心位置处的变形仍然不断变大,并最终形成比较明显的冲切锥破坏。中心区域的集中变形表明板发生了冲切破坏,且发生巨大位移变形的部分主要是加载点周围的单元。根据能量传递原理,破坏会发生在荷载的最近单位,数值模拟的结果正反映了这一现象。比较试验和数值模拟的结果,可以发现二者现象一致,程序计算很好地反映了冲切破坏的特点,即中间加载区域的局部破坏。在冲切承载力方面,根据试验数据模拟的结果进行计算,亦模拟良好。

    图  11  面外冲切变形图
    Figure  11.  Out-of-plane punching deformation

    以M.Rizk的试验HS4数据为例[40],试验结果与模拟结果的对比如图12所示。板单元的面外刚度和强度分为三个阶段:

    图  12  冲切试验曲线模拟图
    Figure  12.  Punching test curve fitting

    1)强度上升段。此阶段荷载随加载位移的增大而增大,板单元在面外的刚度基本保持不变。此时在应力-应变关系曲线中对应的为三折线的第一段,刚度为初始剪切刚度,并未发生折减。此时面外冲切变形较少,弯曲变形占一定比例。

    2)强度上升,刚度折减段。此阶段荷载仍随加载位移的增大而增大,但板单元的面外刚度出现折减。加载区域周围的单元变形增大较为明显,其他部分变形并未出现刚度折减,因而中间出现了局部的位移突变。此时在应力-应变关系曲线中对应的为三折线的第二段,刚度出现了一定程度的折减,但仍然为正。此时面外冲切变形已经较为明显。

    3)强度下降,刚度下降段。此阶段荷载随加载位移的增大出现了下降,即中部区域因为冲切变形过大而产生卸载,分担的荷载明显减小。其他部分继续出现回弹现象,板已经出现了局部的冲切破坏。此时在三折线中对应为第三段,即为负刚度段。冲切破坏现象已经发生,结构局部失效。

    试验结果与数值模拟结果分析一致,三阶段的加载历程正是体现了冲切作用下中部区域的变形特点。对于冲切承载力,目前根据计算公式进行确定,尚无法得到较为精确的值。在后续开发中将针对承载力和面外冲切刚度进行分析,对冲切破坏的承载力和极限位移具体值进行较为准确的预测。

    为了进一步验证本文冲切模型的准确性,选择冲切试验数据库中带星号的试验对应的曲线进行效果验证[4145],比较加载点处的荷载-位移曲线,说明本文的分层壳单元模拟效果的准确性。由于大部分的冲切破坏过程的脆性都比较强,因此,只针对荷载-位移曲线的上升段进行比较。用于对比的冲切数据库中的冲切板参数基本涵盖了工程中常见的钢筋混凝土板的参数范围。其中,最重要的参数便是板的尺寸和混凝土材料的性质,这些在数据库中都得到了很好的考虑。此外,数据库中的冲切板数据也考虑了配筋率等参数对冲切效应的影响。通过对比冲切数据库中的试验曲线,验证本文开发的考虑面外效应的分层壳单元的模拟效果。

    试验与有限元模拟对应的对比结果如图13所示,从模拟结果与试验结果的荷载-位移曲线比较可以看出,基于分层壳单元的三折线冲切本构模型可以很好的模拟板受冲切破坏的过程,对初始刚度,峰值承载力等关键指标都做到了准确的计算。部分冲切板数值模拟的结果与试验结果存在一定的误差,这是由于试验的混凝土特性等存在一定的误差,在实际加载过程中的强度与试验样本测得的强度相差较大。此外,试验条件(如支座约束情况)可能达不到数值模拟中的条件。这些都对初始刚度、峰值荷载等冲切指标都有较大的影响。此外,部分试验受弯承载力设计的较弱,因此在达到冲切承载力时,由于冲切破坏使得整体强度会急剧下降,这是试验曲线中出现平段或下降段的原因。在数值模拟中,由于面外冲切三折线本构的最后一部分斜率绝对值较小,加上面内损伤力学计算弯曲承载力时不会因为冲切破坏而折减过多,因此数值模拟的冲切曲线不会存在急剧的下降段。

    图  13  冲切数据库曲线模拟图
    Figure  13.  Punching database curve fitting

    冲切问题在有限元程序中一直是难以模拟的问题,之前相关本构子程序只能对板的面内力学行为进行模拟,无法计算面外的力学行为。本论文基于分层壳单元针对目前工程中存在较多的面外冲切问题进行了有限元二次开发,利用MSC.MARC的GENSTR函数对冲切问题进行模拟,将之前停留在二维材料本构层面的子程序上升了一个维度,具有更好的普适性。本论文的主要结论如下:

    (1)提出了适用于板的面外分析的三折线模型,对混凝土的面外力学行为进行了模拟。根据分析的结果,程序可以成功地模拟出混凝土板面外从弯曲变形到冲切破坏的全过程,冲切荷载的局部效应也在最终的变形云图中得到了验证;

    (2)通过国内外多组冲切试验数据的分析,对当前已有的混凝土冲切承载力公式进行了比较,得到了较为精确的冲切承载力计算公式,并以此定义了三折线模型的各个关键参数,包括开裂应变、极限应变、极限承载力等;

    (3)为进一步验证本文开发的冲切计算模型的准确性,选择冲切数据库的试验数据,提取加载点处的荷载-位移曲线与试验结果进行比较。从模拟结果与试验结果的荷载-位移曲线比较可以看出,基于分层壳单元的三折线冲切本构模型可以很好的模拟板受冲切破坏的过程。

    附录

      1  冲切试验数据库
    资料来源 板边长 跨度 厚度 ’混凝土强度 配筋率 钢筋屈服强度 加载边长 最大荷载 变形 开裂荷载 开裂前变形
    Sevket Ozden[23]R1500120012021.60.91%507.02001886.35490.25
    Sevket Ozden[23]R1500120012020.01.36%507.02002025.5710.6
    Sevket Ozden[23]R1500120012074.01.87%507.02003315.751070.47
    Sevket Ozden[23]R1500120012075.01.87%507.02003715.431360.64
    Sevket Ozden[23]R1500120012063.72.83%507.02004053.82120
    Sevket Ozden[23]R1500120012074.72.83%507.02004896.421260.58
    Kevin Ka Loo Li[24]92572513539.40.88%460.02003302.6464
    Kevin Ka Loo Li[24]119099019039.40.80%460.02005833.751540.19
    Kevin Ka Loo Li[24]1450125024039.40.80%460.02009044.232360.28
    Kevin Ka Loo Li[24]1975177534539.40.75%460.020013814.344420.4
    Kevin Ka Loo Li[24]1975177545039.40.74%460.030022242.87270
    Kevin Ka Loo Li[24]1975177555039.40.78%460.030026812.793600.09
    H.Marzouk[25, 26]1700150012034.71.74%490.015032014.645
    H.Marzouk[25, 26]1700150012067.00.49%490.015017825.545
    H.Marzouk[25, 26]1700150012070.00.84%490.015024917.654
    H.Marzouk[25, 26]1700150012074.01.19%490.01503561754
    H.Marzouk[25, 26]1700150012069.01.47%490.015035613.163
    H.Marzouk[25, 26]1700150012066.02.37%490.015041814.7
    H.Marzouk[25, 26]1700150015030.00.94%490.015039613.154
    H.Marzouk[25, 26]1700150015068.00.64%490.015036516.9
    H.Marzouk[25, 26]1700150015070.00.94%490.015048914.963
    H.Marzouk[25, 26]1700150015069.01.11%490.015043613.163
    H.Marzouk[25, 26]1700150015074.01.61%490.015054310.863
    H.Marzouk[25, 26]1700150015080.02.33%490.01506459.980
    H.Marzouk[25, 26]170015009070.00.95%490.015019627
    H.Marzouk[25, 26]170015009075.01.52%490.015025826.345
    H.Marzouk[25, 26]170015009068.02.00%490.015026716.245
    H.Marzouk[25, 26]1700150012072.01.47%490.022049817.854
    H.Marzouk[25, 26]1700150012071.01.47%490.030056020.8
    H.Marzouk[25, 26]1900187015067.20.0096490.0250511.52
    H.Marzouk[25, 26]1900187015032.20.0096490.0250475.5
    Kai Qian[28]2200180015036.30.0092420.0200376
    Kai Qian[28]2200180012537.60.011420.0200274
    Mohamed Hassan[29]2400200020045.40.0107470.0300688
    Mohamed Hassan[29]2400200035038.60.0064470.03001692
    Yoshio Kakuta[30]1150115010035.80.98%384.650140.2
    Yoshio Kakuta[30]1150115010039.00.98%384.650130.2
    Yoshio Kakuta[30]1150115010033.50.98%384.6100200.3
    Yoshio Kakuta[30]1150115010036.60.98%384.6100160.3
    Yoshio Kakuta[30]1150115010038.30.98%419.9100170.3
    Yoshio Kakuta[30]1150115010033.20.98%384.6150200.4
    Yoshio Kakuta[30]1150115010035.30.98%384.6150220.4
    Yoshio Kakuta[30]1150115010033.70.98%384.6200240.4
    Yoshio Kakuta[30]1150115010031.40.98%384.6200235.4
    Yoshio Kakuta[30]1150115010031.70.39%384.6100117.8
    Yoshio Kakuta[30]1150115010034.10.39%384.6100131.8
    Yoshio Kakuta[30]1150115010035.01.31%384.6100210.3
    Yoshio Kakuta[30]1150115010033.71.31%384.6100200.3
    Yoshio Kakuta[30]1150115010034.01.57%384.6100222.8
    Yoshio Kakuta[30]1150115010034.11.57%384.6100220.3
    Yoshio Kakuta[30]1150115010035.50.88%347.3100170.8
    Yoshio Kakuta[30]1150115010033.90.88%347.3100190.3
    Yoshio Kakuta[30]1150115010030.50.75%384.6100168.8
    Yoshio Kakuta[30]1150115010031.30.75%384.6100160.3
    Yoshio Kakuta[30]65065010031.61.96%384.6100238.3
    Yoshio Kakuta[30]65065010033.71.96%384.6100212.3
    Yoshio Kakuta[30]1650165010031.80.98%384.6100189.4
    Yoshio Kakuta[30]1650165010032.00.98%384.6100192.4
    Yoshio Kakuta[30]2200220010037.80.98%384.6100160.4
    Yoshio Kakuta[30]2200220010034.60.98%384.6100150.4
    Yoshio Kakuta[30]1150115015031.61.05%384.6100300.6
    Yoshio Kakuta[30]1150115015033.91.05%384.6100291.6
    Yoshio Kakuta[30]1150115020032.81.12%384.6100500.3
    Yoshio Kakuta[30]1150115020032.01.12%384.6100520.3
    Yoshio Kakuta[30]1650165020045.50.79%343.4300660
    Yoshio Kakuta[30]1650165020048.30.79%343.4300660
    Yoshio Kakuta[30]1650165020039.51.01%473.8300750
    Yoshio Kakuta[30]1650165020045.11.01%473.8200750
    Yoshio Kakuta[30]1650165020045.10.95%473.8200635
    Yoshio Kakuta[30]1650165020051.40.95%473.8200600
    Yoshio Kakuta[30]1650165015042.82.68%473.8200400
    Yoshio Kakuta[30]1650165015047.62.68%473.8200660
    Yoshio Kakuta[30]1650165015045.01.68%343.4300540
    Yoshio Kakuta[30]1650165015045.61.68%343.4300710
    Yoshio Kakuta[30]1650165015047.00.89%343.4300600
    Yoshio Kakuta[30]1650165015052.00.89%343.4300600
    Yoshio Kakuta[30]1650165015015.10.87%343.4100190
    Yoshio Kakuta[30]1650165015019.40.87%343.4100190
    Yoshio Kakuta[30]1650165015024.30.87%343.4100240
    Yoshio Kakuta[30]1650165015049.50.87%343.4100312
    Yoshio Kakuta[30]1650165015048.30.87%343.4100341
    Yoshio Kakuta[30]1150115010029.10.98%343.4100174
    Yoshio Kakuta[30]1150115010027.80.98%343.4100165
    Yoshio Kakuta[30]1150115010033.20.98%343.4100177
    Yoshio Kakuta[30]700120010029.30.98%343.4100162
    Yoshio Kakuta[30]700120010029.00.98%343.4100160
    Yoshio Kakuta[30]950120010033.10.98%343.4100170
    Yoshio Kakuta[30]950120010027.90.98%343.4100158
    Ghannoum[27]2300200015037.21.07%380.0225560.67
    Ghannoum[27]2300200015057.11.07%380.0225800.65
    Ghannoum[27]2300200015067.11.07%380.0225890.79
    Inácio[31]16501530125125.60.99%532.3200298
    Inácio[31]16501530125130.11.13%532.3200271
    Inácio[31]16501530125129.61.35%532.3200243
    Inácio[31]1650153012535.91.13%532.3200159
    R. Narayanan[33]7807406041.11.38%550.010084.6
    R. Narayanan[33]7807406049.51.38%550.010091.7
    R. Narayanan[33]7807406042.51.38%550.0100100
    R. Narayanan[33]7807406043.71.38%550.010095.8
    R. Narayanan[33]7807406050.31.38%550.0100113
    R. Narayanan[33]7807406050.41.38%550.0100121
    R. Narayanan[33]7807406044.61.38%550.010091.3
    R. Narayanan[33]7807406043.01.38%550.0100110
    R. Narayanan[33]7807406041.31.38%550.0100110
    R. Narayanan[33]7807406045.21.38%550.0100112
    R. Narayanan[33]7807406028.31.38%550.010083.4
    R. Narayanan[33]7807406030.81.38%550.010081.8
    D. R. C. Oliveira[34]2280210013057.00.95%749.0120204
    D. R. C. Oliveira[34]2280210013057.00.95%750.0120/60283
    D. R. C. Oliveira[34]2280210013057.00.95%751.0120/360320
    N. Banthia[35]6005007541.00.27%1200(FRP)R10064.8
    N. Banthia[35]6005007541.50.27%1200(FRP)R10061.3
    N. Banthia[35]6005007552.90.27%1200(FRP)R10072.1
    L Nguyen-Minh[36]105090012521.70.66%492.0150285
    L Nguyen-Minh[36]1050120012521.70.66%492.0150276
    L Nguyen-Minh[36]1050150012521.70.66%492.0150264
    Stefano[37]3000300025027.61.26%573.0260224
    Stefano[37]3000300025040.50.21%552.0260419
    Stefano[37]3000300025032.20.21%541.0260408
    Stefano[37]3000300025029.30.27%555.0260542
    Stefano[37]3000300025028.50.27%577.0260561
    Stefano[37]3000570050032.40.30%520.0520236
    Stefano[37]3000150412534.71.12%550.0130440
    Stefano[37]3000150412534.70.60%550.0130272
    Stefano[37]1500150413034.70.25%525.0130137
    Stefano[37]150150413034.70.20%525.0130109
    DFW Jing[38]45040012021.34.66%460.040253
    DFW Jing[38]45040012022.74.66%460.040269
    DFW Jing[38]45040012023.14.66%460.040194
    DFW Jing[38]45040012023.6304.66%460.040202
    *张元伟[42]1400120012027.31.28%371.0250356
    *张元伟[42]1400120015024.01.03%371.0250470
    *张元伟[42]1400120018029.40.86%371.0250646
    *张元伟[42]1400120012026.61.28%371.0250360
    *张元伟[42]1400120012033.50.94%371.0250357
    *张元伟[42]1400120012035.80.65%371.0250293
    *洪枫[41]2550235018028.41.43%453.6250690
    *洪枫[41]2550235018039.21.43%453.6250805
    *洪枫[41]2550235018053.20.71%453.6250610
    *洪枫[41]2550235018055.01.07%453.6250785
    *洪枫[41]2550235018056.11.43%453.6250860
    *M.Rizk[39]1900183015042.80.48%400.0250219
    *M.Rizk[39]1900183020042.80.54%400.0250491
    *M.Rizk[39]1900183025033.30.35%400.0250438
    *M.Rizk[39]1900183025066.50.35%400.0250574
    *M.Rizk[39]1900183030038.00.73%400.0250882
    *M.Rizk[39]1900183030042.80.73%400.02501023
    *M.Rizk[39]1900183030071.30.43%400.0250886
    *M.Rizk[39]2650250535072.20.56%400.04001722
    *M.Rizk[39]2650250535066.51.42%400.04002090
    *M.Rizk[39]2650250540038.01.58%400.04002234
    *M.Rizk[39]2650250540057.01.58%400.04002513
    *M.Osman[43]1900183015076.10.50%435.0250304
    *M.Osman[43]1900183015073.41.00%435.0250474
    *M.Osman[43]1900183015075.51.50%435.0250539
    *M.Osman[43]1900183015074.02.00%435.0250613
    *M.Osman[43]1900183015036.21.00%435.0250432
    *M.Osman[43]1900183015037.80.50%435.0250310
    *Ibrahim[44]110090010023.90.90%400.0140138
    *Ibrahim[44]110090010023.91.36%400.0140145
    *Ibrahim[44]110090015023.90.77%400.0140237
    *Ibrahim[44]110090015023.91.60%400.0140288
    *Ibrahim[44]110090010053.30.90%400.0140196
    *Ibrahim[44]110090010053.31.36%400.0140258
    *Ibrahim[44]110090015053.30.77%400.0140388
    *Ibrahim[44]110090015053.31.60%400.0140410
    *Khaled[45]1100103010022.80.55%400.0125219
    *Khaled[45]1100103010048.60.55%400.0125270
    *Khaled[45]1100103010062.30.55%400.0125263   
    注:表中带星号的试验对应的曲线用于数值模拟效果验证。
    下载: 导出CSV 
    | 显示表格
  • 图  1   冲切板建模示意图

    Figure  1.   Modeling of the punched slab

    图  2   MSC.MARC求解准则设定

    Figure  2.   Solution criteria of MSC.MARC

    图  3   面外三折线模型图

    Figure  3.   Out-of-plane trilinear model

    图  4   冲切承载力规范公式对比图

    Figure  4.   Comparison of punching bearing capacity obtained from various standard equations

    图  5   冲切试验破坏曲线图

    Figure  5.   Punching test failure curve diagram

    图  6   初始刚度拟合结果对比图

    Figure  6.   Comparison of initial stiffness fitting results

    图  7   开裂荷载拟合结果对比图

    Figure  7.   Comparison of cracking load fitting results

    图  8   开裂后剪切刚度结果对比图

    Figure  8.   Comparison of post-cracking shear stiffness fitting results

    图  9   冲切破坏变形图

    Figure  9.   Punching failure deformation

    图  10   面内弯曲变形图

    Figure  10.   In-plane bending deformation

    图  11   面外冲切变形图

    Figure  11.   Out-of-plane punching deformation

    图  12   冲切试验曲线模拟图

    Figure  12.   Punching test curve fitting

    图  13   冲切数据库曲线模拟图

    Figure  13.   Punching database curve fitting

    表  1   GENSTR子程序参数变量表

    Table  1   Parameters of GENSTR

    变量名称变量英文含义变量含义
    dgeneralized real stress strain law切线刚度矩阵
    etotatotal strain整体应变
    deincrement of strain整体应变增量
    sstress整体应力
    tstate variables状态参量
    dtincrement of state variables状态参量增量
    ngensnumber of generalized stress应力/应变分量数目
    muser element number用户定义单元数
    nnintegration point number积分点数目
    matus(1)user material number用户定义材料数目
    ht板的厚度
    下载: 导出CSV

    表  2   GENSTR子程序应力应变对应关系表

    Table  2   Stress-strain correspondence of GENSTR

    变量序号变形方向etota(广义应变)s(广义应力)
    1x轴向εxNx
    2y轴向εyNy
    3面内剪切γxyNxy
    4x面外剪切φxMx
    5y面外剪切φyMy
    6x向弯曲φxyMxy
    7y向弯曲γxzVxz
    8x轴向γyzVyz
    下载: 导出CSV

    1   冲切试验数据库

    资料来源 板边长 跨度 厚度 ’混凝土强度 配筋率 钢筋屈服强度 加载边长 最大荷载 变形 开裂荷载 开裂前变形
    Sevket Ozden[23]R1500120012021.60.91%507.02001886.35490.25
    Sevket Ozden[23]R1500120012020.01.36%507.02002025.5710.6
    Sevket Ozden[23]R1500120012074.01.87%507.02003315.751070.47
    Sevket Ozden[23]R1500120012075.01.87%507.02003715.431360.64
    Sevket Ozden[23]R1500120012063.72.83%507.02004053.82120
    Sevket Ozden[23]R1500120012074.72.83%507.02004896.421260.58
    Kevin Ka Loo Li[24]92572513539.40.88%460.02003302.6464
    Kevin Ka Loo Li[24]119099019039.40.80%460.02005833.751540.19
    Kevin Ka Loo Li[24]1450125024039.40.80%460.02009044.232360.28
    Kevin Ka Loo Li[24]1975177534539.40.75%460.020013814.344420.4
    Kevin Ka Loo Li[24]1975177545039.40.74%460.030022242.87270
    Kevin Ka Loo Li[24]1975177555039.40.78%460.030026812.793600.09
    H.Marzouk[25, 26]1700150012034.71.74%490.015032014.645
    H.Marzouk[25, 26]1700150012067.00.49%490.015017825.545
    H.Marzouk[25, 26]1700150012070.00.84%490.015024917.654
    H.Marzouk[25, 26]1700150012074.01.19%490.01503561754
    H.Marzouk[25, 26]1700150012069.01.47%490.015035613.163
    H.Marzouk[25, 26]1700150012066.02.37%490.015041814.7
    H.Marzouk[25, 26]1700150015030.00.94%490.015039613.154
    H.Marzouk[25, 26]1700150015068.00.64%490.015036516.9
    H.Marzouk[25, 26]1700150015070.00.94%490.015048914.963
    H.Marzouk[25, 26]1700150015069.01.11%490.015043613.163
    H.Marzouk[25, 26]1700150015074.01.61%490.015054310.863
    H.Marzouk[25, 26]1700150015080.02.33%490.01506459.980
    H.Marzouk[25, 26]170015009070.00.95%490.015019627
    H.Marzouk[25, 26]170015009075.01.52%490.015025826.345
    H.Marzouk[25, 26]170015009068.02.00%490.015026716.245
    H.Marzouk[25, 26]1700150012072.01.47%490.022049817.854
    H.Marzouk[25, 26]1700150012071.01.47%490.030056020.8
    H.Marzouk[25, 26]1900187015067.20.0096490.0250511.52
    H.Marzouk[25, 26]1900187015032.20.0096490.0250475.5
    Kai Qian[28]2200180015036.30.0092420.0200376
    Kai Qian[28]2200180012537.60.011420.0200274
    Mohamed Hassan[29]2400200020045.40.0107470.0300688
    Mohamed Hassan[29]2400200035038.60.0064470.03001692
    Yoshio Kakuta[30]1150115010035.80.98%384.650140.2
    Yoshio Kakuta[30]1150115010039.00.98%384.650130.2
    Yoshio Kakuta[30]1150115010033.50.98%384.6100200.3
    Yoshio Kakuta[30]1150115010036.60.98%384.6100160.3
    Yoshio Kakuta[30]1150115010038.30.98%419.9100170.3
    Yoshio Kakuta[30]1150115010033.20.98%384.6150200.4
    Yoshio Kakuta[30]1150115010035.30.98%384.6150220.4
    Yoshio Kakuta[30]1150115010033.70.98%384.6200240.4
    Yoshio Kakuta[30]1150115010031.40.98%384.6200235.4
    Yoshio Kakuta[30]1150115010031.70.39%384.6100117.8
    Yoshio Kakuta[30]1150115010034.10.39%384.6100131.8
    Yoshio Kakuta[30]1150115010035.01.31%384.6100210.3
    Yoshio Kakuta[30]1150115010033.71.31%384.6100200.3
    Yoshio Kakuta[30]1150115010034.01.57%384.6100222.8
    Yoshio Kakuta[30]1150115010034.11.57%384.6100220.3
    Yoshio Kakuta[30]1150115010035.50.88%347.3100170.8
    Yoshio Kakuta[30]1150115010033.90.88%347.3100190.3
    Yoshio Kakuta[30]1150115010030.50.75%384.6100168.8
    Yoshio Kakuta[30]1150115010031.30.75%384.6100160.3
    Yoshio Kakuta[30]65065010031.61.96%384.6100238.3
    Yoshio Kakuta[30]65065010033.71.96%384.6100212.3
    Yoshio Kakuta[30]1650165010031.80.98%384.6100189.4
    Yoshio Kakuta[30]1650165010032.00.98%384.6100192.4
    Yoshio Kakuta[30]2200220010037.80.98%384.6100160.4
    Yoshio Kakuta[30]2200220010034.60.98%384.6100150.4
    Yoshio Kakuta[30]1150115015031.61.05%384.6100300.6
    Yoshio Kakuta[30]1150115015033.91.05%384.6100291.6
    Yoshio Kakuta[30]1150115020032.81.12%384.6100500.3
    Yoshio Kakuta[30]1150115020032.01.12%384.6100520.3
    Yoshio Kakuta[30]1650165020045.50.79%343.4300660
    Yoshio Kakuta[30]1650165020048.30.79%343.4300660
    Yoshio Kakuta[30]1650165020039.51.01%473.8300750
    Yoshio Kakuta[30]1650165020045.11.01%473.8200750
    Yoshio Kakuta[30]1650165020045.10.95%473.8200635
    Yoshio Kakuta[30]1650165020051.40.95%473.8200600
    Yoshio Kakuta[30]1650165015042.82.68%473.8200400
    Yoshio Kakuta[30]1650165015047.62.68%473.8200660
    Yoshio Kakuta[30]1650165015045.01.68%343.4300540
    Yoshio Kakuta[30]1650165015045.61.68%343.4300710
    Yoshio Kakuta[30]1650165015047.00.89%343.4300600
    Yoshio Kakuta[30]1650165015052.00.89%343.4300600
    Yoshio Kakuta[30]1650165015015.10.87%343.4100190
    Yoshio Kakuta[30]1650165015019.40.87%343.4100190
    Yoshio Kakuta[30]1650165015024.30.87%343.4100240
    Yoshio Kakuta[30]1650165015049.50.87%343.4100312
    Yoshio Kakuta[30]1650165015048.30.87%343.4100341
    Yoshio Kakuta[30]1150115010029.10.98%343.4100174
    Yoshio Kakuta[30]1150115010027.80.98%343.4100165
    Yoshio Kakuta[30]1150115010033.20.98%343.4100177
    Yoshio Kakuta[30]700120010029.30.98%343.4100162
    Yoshio Kakuta[30]700120010029.00.98%343.4100160
    Yoshio Kakuta[30]950120010033.10.98%343.4100170
    Yoshio Kakuta[30]950120010027.90.98%343.4100158
    Ghannoum[27]2300200015037.21.07%380.0225560.67
    Ghannoum[27]2300200015057.11.07%380.0225800.65
    Ghannoum[27]2300200015067.11.07%380.0225890.79
    Inácio[31]16501530125125.60.99%532.3200298
    Inácio[31]16501530125130.11.13%532.3200271
    Inácio[31]16501530125129.61.35%532.3200243
    Inácio[31]1650153012535.91.13%532.3200159
    R. Narayanan[33]7807406041.11.38%550.010084.6
    R. Narayanan[33]7807406049.51.38%550.010091.7
    R. Narayanan[33]7807406042.51.38%550.0100100
    R. Narayanan[33]7807406043.71.38%550.010095.8
    R. Narayanan[33]7807406050.31.38%550.0100113
    R. Narayanan[33]7807406050.41.38%550.0100121
    R. Narayanan[33]7807406044.61.38%550.010091.3
    R. Narayanan[33]7807406043.01.38%550.0100110
    R. Narayanan[33]7807406041.31.38%550.0100110
    R. Narayanan[33]7807406045.21.38%550.0100112
    R. Narayanan[33]7807406028.31.38%550.010083.4
    R. Narayanan[33]7807406030.81.38%550.010081.8
    D. R. C. Oliveira[34]2280210013057.00.95%749.0120204
    D. R. C. Oliveira[34]2280210013057.00.95%750.0120/60283
    D. R. C. Oliveira[34]2280210013057.00.95%751.0120/360320
    N. Banthia[35]6005007541.00.27%1200(FRP)R10064.8
    N. Banthia[35]6005007541.50.27%1200(FRP)R10061.3
    N. Banthia[35]6005007552.90.27%1200(FRP)R10072.1
    L Nguyen-Minh[36]105090012521.70.66%492.0150285
    L Nguyen-Minh[36]1050120012521.70.66%492.0150276
    L Nguyen-Minh[36]1050150012521.70.66%492.0150264
    Stefano[37]3000300025027.61.26%573.0260224
    Stefano[37]3000300025040.50.21%552.0260419
    Stefano[37]3000300025032.20.21%541.0260408
    Stefano[37]3000300025029.30.27%555.0260542
    Stefano[37]3000300025028.50.27%577.0260561
    Stefano[37]3000570050032.40.30%520.0520236
    Stefano[37]3000150412534.71.12%550.0130440
    Stefano[37]3000150412534.70.60%550.0130272
    Stefano[37]1500150413034.70.25%525.0130137
    Stefano[37]150150413034.70.20%525.0130109
    DFW Jing[38]45040012021.34.66%460.040253
    DFW Jing[38]45040012022.74.66%460.040269
    DFW Jing[38]45040012023.14.66%460.040194
    DFW Jing[38]45040012023.6304.66%460.040202
    *张元伟[42]1400120012027.31.28%371.0250356
    *张元伟[42]1400120015024.01.03%371.0250470
    *张元伟[42]1400120018029.40.86%371.0250646
    *张元伟[42]1400120012026.61.28%371.0250360
    *张元伟[42]1400120012033.50.94%371.0250357
    *张元伟[42]1400120012035.80.65%371.0250293
    *洪枫[41]2550235018028.41.43%453.6250690
    *洪枫[41]2550235018039.21.43%453.6250805
    *洪枫[41]2550235018053.20.71%453.6250610
    *洪枫[41]2550235018055.01.07%453.6250785
    *洪枫[41]2550235018056.11.43%453.6250860
    *M.Rizk[39]1900183015042.80.48%400.0250219
    *M.Rizk[39]1900183020042.80.54%400.0250491
    *M.Rizk[39]1900183025033.30.35%400.0250438
    *M.Rizk[39]1900183025066.50.35%400.0250574
    *M.Rizk[39]1900183030038.00.73%400.0250882
    *M.Rizk[39]1900183030042.80.73%400.02501023
    *M.Rizk[39]1900183030071.30.43%400.0250886
    *M.Rizk[39]2650250535072.20.56%400.04001722
    *M.Rizk[39]2650250535066.51.42%400.04002090
    *M.Rizk[39]2650250540038.01.58%400.04002234
    *M.Rizk[39]2650250540057.01.58%400.04002513
    *M.Osman[43]1900183015076.10.50%435.0250304
    *M.Osman[43]1900183015073.41.00%435.0250474
    *M.Osman[43]1900183015075.51.50%435.0250539
    *M.Osman[43]1900183015074.02.00%435.0250613
    *M.Osman[43]1900183015036.21.00%435.0250432
    *M.Osman[43]1900183015037.80.50%435.0250310
    *Ibrahim[44]110090010023.90.90%400.0140138
    *Ibrahim[44]110090010023.91.36%400.0140145
    *Ibrahim[44]110090015023.90.77%400.0140237
    *Ibrahim[44]110090015023.91.60%400.0140288
    *Ibrahim[44]110090010053.30.90%400.0140196
    *Ibrahim[44]110090010053.31.36%400.0140258
    *Ibrahim[44]110090015053.30.77%400.0140388
    *Ibrahim[44]110090015053.31.60%400.0140410
    *Khaled[45]1100103010022.80.55%400.0125219
    *Khaled[45]1100103010048.60.55%400.0125270
    *Khaled[45]1100103010062.30.55%400.0125263   
    注:表中带星号的试验对应的曲线用于数值模拟效果验证。
    下载: 导出CSV
  • [1]

    TRUONG G T, HWANG H J, KIM C S. Assessment of punching shear strength of FRP-RC slab-column connections using machine learning algorithms [J]. Engineering Structures, 2022, 255: 113898. doi: 10.1016/j.engstruct.2022.113898

    [2]

    ALROUSAN R Z, ALNEMRAWI B R. The influence of concrete compressive strength on the punching shear capacity of reinforced concrete flat slabs under different opening configurations and loading conditions [J]. Structures, 2022, 44: 101 − 119. doi: 10.1016/j.istruc.2022.07.091

    [3]

    MUTTONI A. Punching shear strength of reinforced concrete slabs without transverse reinforcement [J]. ACI Structural Journal, 2008, 105(4): 440 − 450.

    [4] 覃健桂, 朱烨斌, 李易, 等. 基于板壳理论的锈蚀板柱节点变形及冲切承载力计算模型[J/OL]. http://kns.cnki.net/kcms/detail/11.2595.O3.20240711.1451.009.html, 2024-07-15.

    QIN Jiangui, ZHU Yebin, LI Yi, et al. Computational model of deformation and punching shear capacity of corroded slab-column connection upon plate and shell theory [J/OL]. http://kns.cnki.net/kcms/detail/11.2595.O3.20240711.1451.009.html, 2024-07-15. (in Chinese)

    [5] 钱凯, 胡志俊, 翁运昊, 等. 高温下RC板柱节点抗冲切性能数值研究[J/OL]. http://kns.cnki.net/kcms/detail/11.2595.o3.20240527.1307.002.html, 2024-05-28.

    QIAN Kai, HU Zhijun, WENG Yunhao, et al. Numerical study on punching shear performance of RC slab-column connections during high temperature [J/OL]. http://kns.cnki.net/kcms/detail/11.2595.o3.20240527.1307.002.html, 2024-05-28. (in Chinese)

    [6] 邓小芳, 兰冬璆, 翁运昊, 等. 局部替换ECC板柱节点抗冲切性能研究[J/OL]. http://kns.cnki.net/kcms/detail/11.2595.O3.20240304.1237.004.html, 2024-03-07.

    DENG Xiaofang, LAN Dongqiu, WENG Yunhao, et al. Study on punching shear resistance of slab-column connections with partial ECC replacement [J/OL]. http://kns.cnki.net/kcms/detail/11.2595.O3.20240304.1237.004.html, 2024-03-07. (in Chinese)

    [7] 吴丽丽, 罗虎林, 刘斌, 等. 装配式外墙凹槽板受弯力学性能试验研究[J/OL]. http://kns.cnki.net/kcms/detail/11.2595.O3.20231205.1331.014.html, 2023-12-05.

    WU Lili, LUO Hulin, LIU Bin, et al. Experimental study on flexural bearing capacity of prefabricated composite insulation facade panel with groove [J/OL]. http://kns.cnki.net/kcms/detail/11.2595.O3.20231205.1331.014.html, 2023-12-05. (in Chinese)

    [8]

    ZHENG B W, ZHENG W Z, WANG L, et al. Effect of column size on punching behavior of flat slabs with square columns: Numerical investigation [J]. Journal of Building Engineering, 2023, 79: 107937. doi: 10.1016/j.jobe.2023.107937

    [9]

    DE SOUSA A M D, LANTSOGHT E O L, GENIKOMSOU A S, et al. NLFEA of one-way slabs in transition between shear and punching: Recommendations for modeling [J]. Engineering Structures, 2023, 293: 116617. doi: 10.1016/j.engstruct.2023.116617

    [10] 朱大壮, 黄冰峰, 朴贤日, 等. 不同配筋率下钢筋混凝土板柱节点抗冲切承载力试验研究[J]. 工程力学, 2024, 41(增刊1): 66 − 73. doi: 10.6052/j.issn.1000-4750.2023.05.S028

    ZHU Dazhuang, HUANG Bingfeng, PAK H, et al. Experimental study on the punching and shearing resistance of reinforced concrete slab-column joints with different reinforcement ratios [J]. Engineering Mechanics, 2024, 41(Suppl 1): 66 − 73. (in Chinese) doi: 10.6052/j.issn.1000-4750.2023.05.S028

    [11] 樊健生, 王哲, 杨松, 等. 超高性能混凝土板冲切与弯曲性能研究[J]. 工程力学, 2021, 38(4): 30 − 43. doi: 10.6052/j.issn.1000-4750.2020.06.0349

    FAN Jiansheng, WANG Zhe, YANG Song, et al. Research on punching shear and bending behavior of ultra-high performance concrete slabs [J]. Engineering Mechanics, 2021, 38(4): 30 − 43. (in Chinese) doi: 10.6052/j.issn.1000-4750.2020.06.0349

    [12] 易伟建, 刘彪. 混凝土板柱节点冲切承载力的极限分析[J]. 工程力学, 2017, 34(8): 125 − 132, 143. doi: 10.6052/j.issn.1000-4750.2016.03.0195

    YI Weijian, LIU Biao. Ultimate strength analysis of punching shear capacity of concrete slab-column connections [J]. Engineering Mechanics, 2017, 34(8): 125 − 132, 143. (in Chinese) doi: 10.6052/j.issn.1000-4750.2016.03.0195

    [13] GB/T 50010−2010, 混凝土结构设计标准(2024年版)[S]. 北京: 中国建筑工业出版社, 2011.

    GB/T 50010−2010, Code for design of concrete structures [S]. Beijing: China Architecture & Building Press, 2011. (in Chinese)

    [14]

    ACI 318-14, Building code requirements for structural concrete and commentary [S]. ACI, 2014. (查阅网上资料, 未找到本条文献出版地信息, 请确认)

    [15]

    EN 1992-1-1: 2004, Design of concrete structures-Part 1-1: General rules and rules for buildings [S]. Brussels: General Rules and Rules for Buildings, 2004. (查阅网上资料, 未找到本条文献出版信息, 请确认) .

    [16]

    BS 8110-1: 1997, Structural use of concrete. Code of practice for design and construction [S]. London: British Standards Institution, 1997.

    [17]

    TAO M X, WANG Y L, ZHAO J Z, et al. Optimisation of layered shell model for analysis of reinforced concrete shear walls based on machine learning [J]. Journal of Building Engineering, 2023, 77: 107434.

    [18]

    ZHAO J Z, DING R, TAO M X, et al. Experimental and numerical research of ultrahigh performance concrete layered shell element based on a two-dimensional fixed crack model [J]. Structures, 2022, 46: 598 − 610.

    [19]

    TAO M X, WANG Y L, MA Z Y, et al. Experimental database and analysis of in-plane seismic behaviour of double steel-plate composite walls for wind power tower [J]. Thin-Walled Structures, 2024, 200: 111917. doi: 10.1016/j.tws.2024.111917

    [20]

    TAO M X, WANG L B, ZHAO J Z. Out-of-plane punching simulation based on a layered-shell model [J]. Journal of Building Engineering, 2024, 84: 108419. doi: 10.1016/j.jobe.2023.108419

    [21] 丁然, 聂建国, 陶慕轩. 用于钢筋混凝土连梁地震反应分析的考虑非线性剪切的纤维梁单元Ⅰ: 原理与开发[J]. 土木工程学报, 2016, 49(3): 31 − 42.

    DING Ran, NIE Jianguo, TAO Muxuan. A fiber beam-column element considering nonlinear shear effects for seismic response analysis of reinforced concrete coupling beams I: Theory and implementation [J]. China Civil Engineering Journal, 2016, 49(3): 31 − 42. (in Chinese)

    [22] 丁然, 聂建国, 陶慕轩. 用于钢筋混凝土连梁地震反应分析的考虑非线性剪切的纤维梁单元II: 应用与讨论[J]. 土木工程学报, 2016, 49(4): 31 − 39.

    DING Ran, NIE Jianguo, TAO Muxuan. A fiber beam-column element considering nonlinear shear effects for seismic response analysis of reinforced concrete coupling beams II: Application and discussion [J]. China Civil Engineering Journal, 2016, 49(4): 31 − 39. (in Chinese)

    [23]

    OZDEN S, ERSOY U, OZTURAN T. Punching shear tests of normal- and high-strength concrete flat plates [J]. Canadian Journal of Civil Engineering, 2006, 33(11): 1389 − 1400. doi: 10.1139/l06-089

    [24]

    LI K K L. Influence of size on punching shear strength of concrete slabs [D]. Montreal: McGill University, 2000.

    [25]

    MARZOUK H, HUSSEIN A. Experimental investigation on the behavior of high-strength concrete slabs [J]. Structural Journal, 1992, 88(6): 701 − 713.

    [26]

    MARZOUK H, EMAM M, SAMEH HILAL M. Effect of high-strength concrete slab on the behavior of slab-column connections [J]. ACI Structural Journal, 1998, 95(3): 227 − 237.

    [27]

    GHANNOUM C M. Effect of high-strength concrete on the performance of slab-column specimens [D]. Montreal: McGill University, 1998.

    [28]

    QIAN K, LI J S, HUANG T, et al. Punching shear strength of corroded reinforced concrete slab-column connections [J]. Journal of Building Engineering, 2022, 45: 103489. doi: 10.1016/j.jobe.2021.103489

    [29]

    HASSAN M, AHMED E, BENMOKRANE B. Punching-shear strength of normal and high-strength two-way concrete slabs reinforced with GFRP bars [J]. Journal of Composites for Construction, 2013, 17(6): 04013003. doi: 10.1061/(ASCE)CC.1943-5614.0000424

    [30]

    KAKUTA Y, ITOH A, FUJITA Y. Experimental study on punching strength of reinforced concrete slabs [J]. Proceedings of the Japan Society of Civil Engineers, 1974, 1974(229): 105 − 115. doi: 10.2208/jscej1969.1974.229_105

    [31]

    INÁCIO M M G, ALMEIDA A F O, FARIA D M V, et al. Punching of high strength concrete flat slabs without shear reinforcement [J]. Engineering Structures, 2015, 103: 275 − 284. doi: 10.1016/j.engstruct.2015.09.010

    [32]

    MORENO C L, SARMENTO A M. Punching shear analysis of slab‐column connections [J]. Engineering Computations, 2013, 30(6): 802 − 814. doi: 10.1108/EC-Jun-2012-0122

    [33]

    NARAYANAN R, DARWISH I Y S. Punching shear tests on steel-fibre-reinforced micro-concrete slabs [J]. Magazine of Concrete Research, 1987, 39(138): 42 − 50. doi: 10.1680/macr.1987.39.138.42

    [34]

    OLIVEIRA D R C, REGAN P E, MELO G S S A. Punching resistance of RC slabs with rectangular columns [J]. Magazine of Concrete Research, 2004, 56(3): 123 − 138. doi: 10.1680/macr.2004.56.3.123

    [35]

    BANTHIA N, AL-ASALY M, MA S. Behavior of concrete slabs reinforced with fiber-reinforced plastic grid [J]. Journal of Materials in Civil Engineering, 1995, 7(4): 252 − 257. doi: 10.1061/(ASCE)0899-1561(1995)7:4(252)

    [36]

    NGUYEN-MINH L, ROVŇÁK M, TRAN-QUOC T, et al. Punching shear resistance of steel fiber reinforced concrete flat slabs [J]. Procedia Engineering, 2011, 14: 1830 − 1837. doi: 10.1016/j.proeng.2011.07.230

    [37]

    JING D F W, LEONG L T. Effect of aggregate size on punching strength of reinforced concrete slabs [J]. Applied Mechanics and Materials, 2016, 835: 450 − 454. doi: 10.4028/www.scientific.net/AMM.835.450

    [38]

    GUANDALINI S, BURDET O, MUTTONI A. Punching tests of slabs with low reinforcement ratios [J]. ACI Structural Journal, 2009, 106(1): 87 − 95.

    [39]

    DIN 1045-1: 2008, Concrete, reinforced and prestressed concrete structures-part 1: Design and construction [S]. Berlin: Normenausschuss Bauwesen (NABau) im DIN Deutsches Institut für Normung e. V. Beuth Verl, 2008. (查阅网上资料, 未找到本条文献出版信息, 请确认) .

    [40]

    RIZK E R. Structural behaviour of thick concrete plates [D]. St Johns: Memorial University of Newfoundland, 2010.

    [41] 洪枫. 钢筋混凝土板柱节点冲切破坏试验研究[D]. 长沙: 湖南大学, 2015.

    HONG Feng. Experimental study on punching shear failure of reinforced concrete slabs-column connection [D]. Changsha: Hunan University, 2015. (in Chinese)

    [42] 张元伟. 钢筋混凝土板抗冲切试验研究[D]. 长沙: 湖南大学, 2009.

    ZHANG Yuanwei. Punching shear research on reinforced concrete slabs [D]. Changsha: Hunan University, 2009. (in Chinese)

    [43]

    OSMAN M, MARZOUK H, HELMY S. Behavior of high-strength lightweight concrete slabs under punching loads [J]. ACI Structural Journal, 2000, 97(3): 492 − 498.

    [44]

    METWALLY I M, ISSA M S, EL-BETAR S A. Punching shear resistance of normal and high strength reinforced concrete flat slabs [J]. Civil Engineering Research Magazine, 2008, 30(3): 982 − 1004.

    [45]

    RAGAB K S, EISA A S, ZAKI S I. Improving punching shear resistance of slab column connections by using high strength self compacting concrete with steel fibers [C]// Proceedings of the 8th International Conference on Concrete in the Low Carbon Era. Dundee, Scotland, UK, 2012. (查阅网上资料, 未找到本条文献出版者和页码信息, 请确认补充) .

图(13)  /  表(3)
计量
  • 文章访问数:  20
  • HTML全文浏览量:  13
  • PDF下载量:  3
  • 被引次数: 0
出版历程
  • 收稿日期:  2024-08-05
  • 修回日期:  2024-12-10
  • 网络出版日期:  2024-12-22

目录

/

返回文章
返回