Processing math: 69%

次边缘柱失效下预应力混凝土框架抗连续倒塌性能研究

兰冬璆, 邓小芳, 金浏, 钱凯

兰冬璆, 邓小芳, 金浏, 钱凯. 次边缘柱失效下预应力混凝土框架抗连续倒塌性能研究[J]. 工程力学, 2023, 40(10): 58-70. DOI: 10.6052/j.issn.1000-4750.2022.01.0052
引用本文: 兰冬璆, 邓小芳, 金浏, 钱凯. 次边缘柱失效下预应力混凝土框架抗连续倒塌性能研究[J]. 工程力学, 2023, 40(10): 58-70. DOI: 10.6052/j.issn.1000-4750.2022.01.0052
LAN Dong-qiu, DENG Xiao-fang, JIN Liu, QIAN Kai. PROGRESSIVE COLLAPSE PERFORMANCE OF PRESTRESSED CONCRETE FRAMES UNDER PENULTIMATE COLUMN REMOVAL SCENARIO[J]. Engineering Mechanics, 2023, 40(10): 58-70. DOI: 10.6052/j.issn.1000-4750.2022.01.0052
Citation: LAN Dong-qiu, DENG Xiao-fang, JIN Liu, QIAN Kai. PROGRESSIVE COLLAPSE PERFORMANCE OF PRESTRESSED CONCRETE FRAMES UNDER PENULTIMATE COLUMN REMOVAL SCENARIO[J]. Engineering Mechanics, 2023, 40(10): 58-70. DOI: 10.6052/j.issn.1000-4750.2022.01.0052

次边缘柱失效下预应力混凝土框架抗连续倒塌性能研究

基金项目: 国家自然科学基金优秀青年基金项目(52022024);国家自然科学基金面上项目(51778153);广西自然科学基金项目 (2021GXNSFFA196001)
详细信息
    作者简介:

    兰冬璆(1995−),女,广西柳州人,博士生,主要从事结构抗倒塌方面研究(E-mail: landongqiu@emails.bjut.edu.cn)

    金 浏(1985−),男,江苏泗阳人,教授,博士,博导,主要从事混凝土结构防灾减灾方面研究(E-mail: jinliu@bjut.edu.cn)

    钱 凯(1983−),男,浙江宁波人,教授,博士,博导,主要从事结构抗倒塌方面研究(E-mail: qiankai@glut.edu.cn)

    通讯作者:

    邓小芳(1984−),女,湖北荆州人,助理教授,博士,硕导,主要从事混凝土与混凝土结构研究(E-mail: dengxiaofang@glut.edu.cn)

  • 中图分类号: TU378.4;TU311.41

PROGRESSIVE COLLAPSE PERFORMANCE OF PRESTRESSED CONCRETE FRAMES UNDER PENULTIMATE COLUMN REMOVAL SCENARIO

  • 摘要: 考虑到去柱后弯矩方向在去柱处出现反转并且可能与预应力筋布置形式相反,预应力混凝土结构具有较大的倒塌风险,该文通过拟静力试验和数值模拟研究了次边缘柱失效下无黏结预应力混凝土(PC)梁-柱子结构的抗连续倒塌性能。对2个具有不同预应力筋布置形式(直线型和抛物线型)的1/2缩尺梁-柱子结构进行Pushdown加载,分析了试件的破坏模式、抗力机制演变和极限承载力。为量化预应力筋的影响,制作了一个具有相同几何尺寸和非预应力筋配筋率的钢筋混凝土试件作为参考。结果表明:预应力筋可以显著提高结构承载力,同时改变结构的破坏模式,降低结构变形能力。预应力筋产生的水平拉力加剧了边柱的挠曲P-δ效应,使边柱提早发生大偏心受压破坏。为进一步了解预应力混凝土框架抗倒塌性能,基于有限元软件LS-DYNA建立了精细化有限元模型开展系列参数分析。数值分析表明,预压应力可以增强压拱机制,预应力筋对于大变形阶段承载力的提高贡献显著。
    Abstract: Considering that the direction of bending moment at the beam end nearby the removed column will reverse, which may disagree with the profile of prestressing strands. Thus, the prestressed concrete frame may have high vulnerability for progressive collapse. In this study, experimental and numerical studies were performed to study the performance of unbonded prestressed concrete (PC) beam-column sub-assemblies subjected to a penultimate column removal scenario to resist progressive collapse. Two half-scale PC beam-column sub-assemblies with different strand profiles were tested by Pushdown loading procedure. Failure modes, the involution of load resisting mechanisms, and ultimate load resistance were analyzed. To quantify the effect of the strand, a non-prestressed reinforced concrete beam-column sub-assembly was designed, the geometrical dimension and reinforcement ratio of which were consistent with the PC beam-column sub-assemblies. It is found that the presence of the strand leads to a greater load resistance, lower deformation capacity and different failure modes. The considerable tension in the strand cause significant second-order effects in the side column, leading to large eccentric compression failure. To further study the performance of PC frames against their progressive collapse, high-fidelity numerical models built by LS-DYNA were used to conduct a series of parameter studies. It is found that prestress forces can increase compressive an arch action. The strand can significantly enhance the ultimate load resistance at a large deformation stage.
  • 为满足我国“一带一路”及西部开发战略需要,近年来西部交通基础设施飞速发展,山岭隧道建设数量也不断增多。如国家最新批复的川藏铁路,全线隧道占比在80%以上。同时,中国西部地区强震频发,对既有或在建隧道结构安全存在较大威胁。如在2008年汶川地震时,位于震区都汶公路上的11座隧道均出现不同程度的损毁[1-2],其中,36.2%的隧道出现中等及以上破坏,包括洞口边坡崩塌与滑塌、衬砌及围岩坍塌、衬砌开裂错位、初期支护变形及开裂等。另一方面,凸起地形对地震波有显著的放大作用,而山体中隧道的存在会产生明显的动应力集中。

    国内外学者对凸起地形或山体-隧道地震反应已有较多研究。首先,在凸起地形地震动放大效应方面,袁晓铭等[3]、Tsaur等[4]、Liang等[5]、梁建文等[6]和杨在林等[7]给出了凸起及凹陷地形对地震波散射的解析解;考虑到实际复杂场地特征,刘晶波等[8-9]采用有限元法给出了山体地形地震反应问题的数值解;巴振宁等[10-11]结合“分区契合”技术,采用间接边界元方法研究了凸起地形对平面P波及SH波的散射问题。在隧道反应方面,Liu等[12]研究了弹性半空间浅埋隧道对平面P波、SV波的散射,王冬等[13]采用快速多级间接边界元法揭示了P波、SV波入射下三维孔洞群周围位移场和应力场的分布规律;在山体-隧道整体反应方面,李敏等[14]根据辅助函数的思想,利用复变函数、多级坐标法的方法研究了SH波入射时含圆形孔洞的半圆形凸起地形与其附近的圆形空洞的相互作用。罗昊等[15]采用波函数展开法针对凸起地形中衬砌隧道对平面SH波的散射进行了研究;Amornwongpaibun等[16]研究了含单个椭圆隧道的凸起地形对平面SH波的散射;Liu等[17]采用间接边界元法(IBEM)研究了平面SV波作用下衬砌隧道与山体的动力相互作用。巴振宁等[18]采用间接边界元法研究了山岭-单个隧道在平面P-SV波入射下的地震动力响应问题。

    本文基于弹性波动理论,结合“分区契合”思想,利用间接边界元方法,研究了SV波入射下山体-双线隧道整体动力响应特征。通过大量参数分析,揭示了山体中双线隧道地震动应力集中及山体表面地震动放大效应,所得结论对于实际山岭隧道抗震设计、地震安全性评价具有参考价值。

    图1所示,以弹性半空间上一高斯形山体中含两马蹄形衬砌隧道为例,假设山体及隧道中为均匀各向同性、完全弹性介质。取R1=22mR2=8.4mR3=3m,衬砌厚度为t=0.6m,高斯山体高度为h,宽度为2a图2为双马蹄形山岭隧道整体计算模型的分区图,基于“分区契合[6]”的思想,将整个计算模型划分成半空间凹陷域Θ、山体凸起域Ω和左右隧道域Ψ1Ψ2四个区域。其中,半空间凹陷域Θ是由水平地表L1和虚拟凹陷边界L2围成的半空间无限区域;山体凸起域Ω是由虚拟凹陷边界L2、凸起边界L3及左右隧道外边界L4L6围成的区域;左右隧道域Ψ1Ψ2分别由隧道内外壁边界L4L5L6L7围成。

    图  1  山体隧道模型图(双向六车道)
    Figure  1.  Mountain-tunnel model (two-way six-lane)
    图  2  计算模型的区域划分
    Figure  2.  Area division of calculation model

    本文以全空间中二维水平和竖向均布力为基本解,采用间接边界元法求解双山岭隧道对SV波的二维散射问题。以边界面上虚拟均布荷载分别构造衬砌内外的散射波场,由衬砌内外表面的边界条件建立积分方程,通过离散求解得到虚拟荷载密度,进而求得总波场。设平面简谐SV波传播到山体中,经两隧道内外壁或者凸起地表及水平地表散射后可以得到的平面内位移u(t),极坐标oθ中需满足波动方程:

    2u(t)r2+1ru(t)r+1r22u(t)θ2=1β22u(t)t2 (1)

    V代表一个质地均匀、各向同性和完全弹性的区域,S表示该区域的边界或者表面。地震波入射下,弹性区域内每一点的位移都能用Somigliana积分表达式[19]表示为:

    cui(ξ)=VGij(y,ξ)fj(y)dV+S[Gij(x,ξ)tj(x)Tij(x,ξ)uj(x)]dS (2)

    式中:fi(y)为区域Vy点的体力;uj(x)为区域V中任一点x的位移;ti(x)为任意点x的应力;Gij(x,ξ)Tij(x,ξ)分别为任一点x的位移、应力格林函数[19]

    基于单层位势理论,按照间接边界元法基本原理,各点位移、应力[19]可表示如下:

    ui(x)=VGij(x,ξ)fj(ξ)dVξ+SGij(x,ξ)ϕj(ξ)dSξ (3)
    ti(x)=12δijϕi(x)+VTij(x,ξ)fj(ξ)dVξ+STij(x,ξ)ϕj(ξ)dSξ (4)

    式(4)中,ϕi(x)表示为区域边界上的虚拟荷载分布。进行数值计算,必须先离散所求计算模型的边界面,求解过程中,通过忽略简谐波的时间因子eiωt来构造格林函数,其中,i为虚数单位,i2=s1t代表时间;ω代表圆频率,对非均布荷载作用单元,采用单点高斯积分,集中力位移和应力格林函数可以表达为:

    Gij=1i8ρ[δijA(2γiγjδij)B],i,j=1,2,3 (5)
    Tij=iμ2ρr{[B+λD(qr)2μα2]γjnj+[B+D(kr)2β2]×[γjnj+γknkδij]+(C4B)γiγjγknk} (6)

    式(5)、式(6)中:A=H(2)0(qr)α2+H(2)0(kr)β2B=H(2)2(qr)α2+H(2)2(kr)β2C=D(qr)α2D(kr)β2D(p)=pH(2)1(p)α=(λ+2μ)/ρβ=μ/ρ分别为P波、SV波在介质中的波速;q=ω/αk=ω/β分别为P波、SV波在介质中的波数;δij表示Delta函数(δij=1,ij);μ为介质剪切模量;γj=(xjξj)/rr=(x1ξ1)2+(x3ξ3)2nj是法向量;H(2)n()表示n阶第二类Hankel函数。

    根据单层位势理论和IBEM原理,散射场的产生是由于在各个区域边界施加虚拟均布单元荷载,即把虚拟波源置于离散边界单元上,在这里散射场引起散射体内任一点的位移应力积分表达式为:

    半空间凹陷域Θ内:

    u(s,Θ)i(x)=S1ϕj(ξ)Gij(x,ξ)dSξ (7)
    t(s,Θ)i(x)=S1ϕj(ξ)Tij(x,ξ)dSξ (8)

    凸起山体域Ω内:

    u(s,Ω)i(x)=S2ϕj(ξ)Gij(x,ξ)dSξ (9)
    t(s,Ω)i(x)=S2ϕj(ξ)Tij(x,ξ)dSξ (10)

    左右衬砌隧道域Ψ1Ψ2内:

    u(s,Ψ1(Ψ2))i(x)=S3ϕj(ξ)Gij(x,ξ)dSξ (11)
    t(s,Ψ1(Ψ2))i(x)=S3ϕj(ξ)Tij(x,ξ)dSξ (12)

    式中:ϕjdSξ(i,j=x,y)为离散边界单元的应力;Gij(x,ξ)Tij(x,ξ)分别代表位移和应力格林函数,指的是单位力向量j作用在ξ引起xi向的应力和位移,格林函数自动满足水平地表边界上应力为零的边界条件,同时满足波动方程。

    由于采用弹性全空间动力基本解,自由地表边界条件需构造满足,另外还需考虑隧道与山体交界面的位移、应力连续性条件及隧道内壁等的零应力条件,具体为:

    1) 自由水平地表边界L1上应力为0自动满足:

    τ(t)Θ,ij=μu(t)rr=0,i,j=x,y (13)

    2) 虚拟凹陷边界L2上的位移和应力连续:

    u(f)Θ,ij+u(s)Θ,ij=u(s)Ω,ij,i,j=x,y (14)
    τ(f)Θ,ij+τ(s)Θ,ij=τ(s)Ω,ij,i,j=x,y (15)

    3) 凸起山体域Ω内凸起边界L3上应力为0:

    τ(t)Ω,ij=μu(t)rr=0,i,j=x,y (16)

    4) 凸起山体域Ω和左右隧道Ψ1Ψ2域间的公共边界L4L6上应力连续:

    u(f)Ω,ij+u(s)Ω,ij=u(s)Ψ1(Ψ2),ij,i,j=x,y (17)
    τ(f)Ω,ij+τ(s)Ω,ij=τ(s)Ψ1(Ψ2),ij,i,j=x,y (18)

    5) 左右隧道Ψ1Ψ2域中的隧道衬砌内边界L5L7应力为0:

    τ(t)Ψ1(Ψ2),ij=μu(t)rr=0,i,j=x,y (19)

    根据计算模型边界条件式(13)~式(19)和散射场位移应力表达式(7)~式(12),边界条件的积分表达式为:

    SϕΘj(ξ)TΘij(x,ξ)dSξ=t(f)Θ,i,j=x,y (20)
    SϕΘj(ξ)GΘij(x,ξ)dSξSϕΩj(ξ)GΩij(x,ξ)dSξ=u(f)Θ,i,j=x,y (21)
    SϕΘj(ξ)TΘij(x,ξ)dSξSϕΩj(ξ)TΩij(x,ξ)dSξ=t(f)Θ,i,j=x,y (22)
    SϕΩj(ξ)TΩij(x,ξ)dSξ=0,i,j=x,y (23)
    SϕΩj(ξ)GΘij(x,ξ)dSξSϕΨ1+Ψ2j(ξ)GΨ1+Ψ2ij(x,ξ)dSξ=0,i,j=x,y (24)
    SϕΩj(ξ)TΘij(x,ξ)dSξSϕΨ1+Ψ2j(ξ)TΨ1+Ψ2ij(x,ξ)dSξ=0,i,j=x,y (25)
    SϕΨ1(Ψ2)j(ξ)TΨ1(Ψ2)ij(x,ξ)dSξ=0,i,j=x,y (26)

    实际计算时,首先,离散山体附近部分地表、凸起边界和左右隧道内外壁边界,再将虚拟波源置于每个离散边界单元上,波源强度保持不变。考虑散射波的空间衰减特征,取水平地表的离散范围大致为隧道附近波长的5倍~8倍,首先将所有边界L1~L7对应离散成N1~N7,假设ϕj(ξ)在每一个离散的单元上的强度保持不变,式(20)~式(26)可转化为线性方程组:

    N1+N2+N1l=1ϕΘj(ξl)tΘij(xn1,ξl)=t(f)Θ,ij(xn1),n1=1,N1,i,j=x,y (27)
    N1+N2+N1l1=1(ϕΘj(ξl1)gΘij(xn2,ξl1))N2+N3+N4+N6l2=1(ϕΩj(ξl2)gΩij(xn2,ξl2))=u(f)Θ,ij,n2=1,N2,i,j=x,y (28)
    N1+N2+N1l1=1(ϕΘj(ξl1)tΘij(xn2,ξl1))N2+N3+N4+N6l2=1(ϕΩj(ξl2)tΩij(xn2,ξl2))=u(f)Θ,ij,n2=1,N2,i,j=x,y (29)
    N2+N3+N4+N6l=1ϕΩj(ξl)tΩij(xn3,ξl)=0,n3=1,N3,i,j=x,y (30)
    N2+N3+N4+N6l1=1(ϕΩj(ξl1)gΩij(xn4,ξl1))N4+N5l2=1(ϕΨ1j(ξl2)gΨ1ij(xn4,ξl2)N6+N7l3=1ϕΨ2j(ξl3)gΨ2ij(xn6,ξl3))=0,n4=1,n6=1,N4+N6,i,j=x,y (31)
    N2+N3+N4+N6l1=1(ϕΩj(ξl1)tΩij(xn4,ξl1))N4+N5l2=1(ϕΨ1j(ξl2)tΨ1ij(xn4,ξl2)N6+N7l3=1ϕΨ2j(ξl3)tΨ2ij(xn6,ξl3))=0,n4=1,n6=1,N4+N6,i,j=x,y (32)
    N4+N5(N6+N7)l=1ϕΨ1(Ψ2)j(ξl)tΨ1(Ψ2)ij(xn6,ξl)=0,n5(n7)=1,N5(N7),i,j=x,y (33)

    根据弹性波动理论,总波场等于自由场和散射场的叠加,散射场由散射体表面离散单元上施加虚拟荷载来构造:

    ui=u(0)i+u(d)i (34)

    式中,u(0)i表示由入射简谐波引起的自由场(除去任意截面形式浅埋隧道后的部分)位移。

    基于上述讨论,忽略重力对计算结果的影响,式(3)、式(4)可简化为式(35)、式(36),即:

    u(d)i(x)=Sϕj(ξ)Gij(x,ξ)dSξ (35)
    12δijϕj(x)+Sϕj(ξ)Tij(x,ξ)dSξ=t(0)i(x) (36)

    将散射体表面边界S离散成N个均等的小单元,每个小单元的的尺寸用ΔS代表,则线性方程可表示为:

    Nl=1ϕj(ξi)tij(xn,ξl)=t(0)i(xn),n=1,N, (37)
    tij(xn,ξl)=12δijξnl+ξl+ΔS2ξlΔS2Tij(xn,ξ)dSξ (38)

    式(37)中,如果xξ,可直接利用高斯积分求解tij(xn,ξl),反之需要利用格林函数展开来得到解析解,式(38)表示为:

    tij(xn,ξn)=12δij (39)

    联立式(37)、式(38)可求解出离散边界上的虚拟均布波源密度ϕj(ξi),进而可求解散射场位移。

    由于自由地表边界条件处理上的困难,半空间中山岭隧道对SV波的散射问题至今还没有精确的解析解,只能通过边界条件验算及退化到半圆形凸起山体动力响应与现有结果的比较来考察计算精度。

    首先定义一无量纲频率η:

    η=2Rλ=kRπ=ωRπβ (40)

    式中:R为隧道的宽度;λ为剪切波波长;k为剪切波波数;β为入射波剪切波速。

    现考虑退化为单一隧道,利用本文的方法计算了隧道内表面的半空间位移幅值和动应力集中系数。图3(a)为本文方法结果与Luco和De Barros[20]关于SV波垂直入射的弹性半空间的结果比较。计算参数:密度比ρ2/ρ1=1.0;波速比cβ2/cβ1=1.0;阻尼比ζ=0.001;泊松比ν=1/3;无量纲频率η=0.5。另外,考虑退化为单一半椭圆山体,计算山体周围的频率响应,图3(b)为本文方法和Álvarez-Rubio[21]等人关于SV波入射角为30°时的结果比较。计算参数:阻尼比ζ=0.01;泊松比ν=0.25。对比两者吻合程度良好,从而验证了本文方法的正确性。

    图  3  精度验证
    Figure  3.  Accuracy verification

    山体模型尺寸参数如表1,其中,取山体、隧道半径比a/a1=10,高斯山体地形和隧道衬砌材料泊松比分别均取0.25和0.2,剪切波速分别取1000 m/s(微风化岩)和2667 m/s,材料粘滞阻尼比均取0.001。图4图5分别给出了不同隧道间距D对应的山体地表水平位移幅值和竖向位移幅值。图6给出了SV波不同角度入射下山体地表不同点位移幅值谱。计算参数:双隧道间距D/a1=3、4、5、6;无量纲频率η=1、2、5、10;入射角度α=0°、30°。

    表  1  山体模型尺寸参数
    Table  1.  Size parameter of mountain model
    形状几何形状
    高斯山体|x|\begin{array}{l} - a < x \leqslant 0,y = - 0.5a[ { - 2{{(x/a)}^3} - 3{{(x/a)}^2} + 1} ] \\ 0 < x < a,y = - 0.5a[ {2{{(x/a)}^3} - 3{{(x/a)}^2} + 1} ] \\ \end{array}
    下载: 导出CSV 
    | 显示表格
    图  4  不同隧道间距D影响的地表位移幅值(\alpha = 0°)
    Figure  4.  Surface displacement amplitude of mountain affected by different tunnel spacings D (\alpha = 0°)
    图  5  不同隧道间距D影响的地表位移幅值(\alpha = 30°)
    Figure  5.  Surface displacement amplitude of mountain affected by different tunnel spacings D (\alpha = 30°)
    图  6  SV波入射下山体地表不同点位移幅值谱
    Figure  6.  Displacement spectra of different points of mountain surface under SV wave incidence

    本文给出了稳态平面SV波入射下,山体表面及附近地表表面位移幅值图。分析不同隧道间距D及不同入射频率对山体表面及地表水平位移的影响程度。图4图5中横、纵轴坐标分别对应地表点位x/a和地表位移幅值(已由入射波的位移幅值标准化)。

    图4图5可以看出,无量纲频率\eta 、入射角度\alpha 、双隧道间距D对地表位移幅值具有重要影响。1)当\eta 较小时(\eta = 1, 2),山体表面及附近地表位移幅值空间变化简单;当入射频率逐渐增大时,地表位移幅值也随之增大且震荡更为剧烈,空间分布特征更为复杂;当\alpha = 0°,山体表面以及附近水平地表位移幅值由\eta =1时的2.76增大到\eta =10时的4.85,竖向地表位移幅值则由0.94增大到2.31;当以\alpha = 30°斜入射时,水平位移由\eta =1时的2.34增大到\eta =10时的3.28,竖向位移则由2.29增大到3.16。2)接近垂直入射时,山体表面及附近地表位移幅值表现出来很强对称性,随着入射角度\alpha 的增大,迎波面地表位移幅值震荡加剧。值得注意的是高频SV波入射时,迎波面地表位移震荡剧烈,但位移幅值出现在背波面,实际工程中应当适当提高隧道沿线建筑物的抗震设防等级。3)双隧间距D也是重要影响因素:当入射波频率较小时,隧道间距D变化对地表位移影响不大,如图中当无量纲频率\eta =1时,隧道间距D/a1由3增加到6,4条地表位移曲线基本重合,随着\eta 的增大,双隧间距D变化对位移幅值的影响增强;由图易知当\eta > 5时,双隧道间距D增大地表位移幅值随之增大,当\alpha =0°时,山顶水平位移幅值由D/a1=3时的3.49增加到D/a1=6时的4.62。

    为全面反映散射的频谱特性,分析SV波入射下不同入射角度\alpha 、不同入射频率\eta 、不同双隧间距D对山体地表典型点位的位移幅值谱的影响。考虑地震波入射角度\alpha = 0°、30°,在山体地表选择x/a=0、0.5、1、2、4这5个观察点。其它地质条以及几何参数同上,入射波无量纲频率\eta 取值范围为0~10。

    图6可看出,无量纲频率\eta 、入射角度\alpha 、双隧间距D对高斯山体地表位移频谱特征具有显著影响:1)当SV波垂直入射(\alpha = 0°)时,山顶位置(x/a=0)水平位移幅值明显大于其他点位,最大位移幅值{u_{\max }}=6.67(\eta = 4.53,山顶处);而山顶位置的竖向位移幅值明显小于其他点位,山顶附近(x/a=0.5)位置的地震响应更为明显。随入射角度\alpha 增加,水平位移整体上有减小的趋势,竖向位移有增大的趋势,且地表竖向位移幅值谱曲线震动加剧;2)山体及其附近地表幅值谱曲线对频率变化较为敏感,会出现明显的“波峰”和“波谷”,频率\eta 增加,位移幅值有增大趋势,如最大位移幅值由\eta = 0时的1.87增加到\eta = 4.53时的6.67;3)整体上看,SV波入射下,最大地表位移幅值出现在山顶处,随着地表点位到山顶距离x/a增大,地表位移幅值显著减小,当距离x/a>2后,位移幅值谱逐步趋于稳定。

    为了更加直观的表示凸起山体地形中双隧道由于入射波特性、双隧间距D等地震参数对其动力响应的影响,本文分别给出了不同地震参数下隧道上的位移幅值云图和动应力幅值云图。定义隧道衬砌内部环向应力和入射波应力幅值的比值为动应力集中因子,即\sigma _{\alpha \alpha }^* = \left| {{\sigma _{\alpha \alpha }}/\mu k_\beta ^2} \right|。为了考虑双隧道间距D对隧道动力响应的影响,取D/a1=3、4、5。考虑入射波特性对隧道动力响应的影响,取无量纲频率\eta =1、2、5、10。入射角度\alpha = 0°、30°。

    图7给出了SV波垂直入射下弹性半空间中凸起山体地形中含双马蹄形衬砌隧道的水平位移分布云图情况,鉴于垂直入射时,两隧道的位移云图关于Y轴对称,因此,仅画出右侧隧道的位移云图(下文同)。图7中的位移幅值已由入射波的位移幅值标准化。

    图  7  不同隧道间距D影响的隧道水平位移云图
    Figure  7.  Horizontal displacement contours of tunnel affected by different tunnel spacings D

    图7(a)图7(b)图7(c)分别表示D/a1=3、4、5,纵向对比可知,随着双隧道间距D的增大,山岭隧道衬砌内部的位移幅值整体上有减小的趋势,以\eta =1.0为例,最大位移幅值由D/a1=3时的0.85减小到D/a1=5时的0.75,隧道的空间位移有关于隧道中心线对称的趋势,两隧道间的相互作用减弱;当入射频率较低时,山岭隧道水平位移最大值出现在右隧道顶部,如\eta =1.0时,右山岭隧道水平位移为0.8;随着无量纲频率的增大,地震波波长减小,隧道内部对地震波的散射更加强烈,山岭隧道衬砌内部位移幅值整体上逐渐增大,\eta =10.0时,右山岭隧道水平位移为1.4,隧道位移的空间分布特征也更为复杂。

    图8~图10分别给出了凸起山体地形中双山岭隧道在SV波入射下,双隧道间距不同时左右衬砌隧道内部动应力集中因子分布情况。SV波垂直入射时仍仅给出右隧道应力云图。可以看出,双隧道衬砌内部的应力幅值及空间状态分布主要与入射波的特性(入射波频率和入射角)和隧道间距D有关。

    图  8  SV波入射下不同隧道间距时右山岭隧道环形应力云图(\alpha = 0°)
    Figure  8.  Circular stress contours of right mountain tunnel with different tunnel spacings under SV wave incidence (\alpha = 0°)
    图  9  SV波入射下不同隧道间距时右山岭隧道环形应力云图(\alpha = 30°)
    Figure  9.  Circular stress contours of right mountain tunnel with different tunnel spacings under SV wave incidence (\alpha = 30°)
    图  10  SV波入射下不同隧道间距时左山岭隧道环形应力云图(\alpha = 30°)
    Figure  10.  Circular stress contours of left mountain tunnel with different tunnel spacings under SV wave incidence (\alpha = 30°)

    分析图8可得,SV波垂直(\alpha = 0°)入射,\eta =1.0时,右隧道的环向应力空间分布比较均匀,关于隧道中心线对称,主要是因为在低频段两隧道间的相互作用比较弱,此时应力峰值出现在隧道上部左右两侧与竖向夹角约为45°及隧道下部左右两侧与水平夹角约为30°的位置,其值约为25;随入射频率\eta 增大,对于右隧道,环向应力峰值逐渐向隧道右侧移动,在\eta =10时,隧道内部应力分布更加复杂,应力集中区域增加但隧道衬砌内部环向应力幅值逐渐减小,动应力集中因子衰减到20;纵向对比可知,随隧道间距D增大,隧道应力空间分布关于隧道中心线对称,说明隧道间相互影响减弱,故实际工程中,双隧道间距D增大到一定值后,可忽略隧道间相互作用。

    观察图9图10可得,随着双隧间距D的增加,隧道内部衬砌环向应力有减小的趋势,这与分析图8得出的结论相一致。左右两侧隧道应力计算结果表明,随着入射角度\alpha 的增大,两隧道衬砌内部环向应力均有增加的趋势,由于山体对地震波复杂的反射作用,左隧道的应力幅值整体上小于右隧道,实际工程设计时应注意远离地震波入射一侧隧道的强度破坏。另外,隧道的内外壁的应力数值差别较大,故实际工程设计中要多加注意。

    本文采用间接边界元法研究了SV波入射下弹性半空间中双线隧道—山体系统动力响应,着重考察了隧道存在对山体表面及附近地表散射的影响规律,分析了入射波角度、频率、隧道间距等参数对散射的影响,得到如下有益结论:

    (1)入射角度、无量纲频率、山体高宽比、双隧间距等都是影响山体表面、附近地表地面运动、山体隧道本身应力及位移的主要因素。

    (2)当无量纲频率\eta 较小时(\eta =1, 2),山体表面及其附近地表位移幅值的空间变化比较简单,双隧道间距D变化对地表位移幅值影响不大;当入射频率\eta 逐渐增大时,地表位移幅值也随之增大,同时地表位移幅值震荡更为剧烈,空间分布特征更为复杂。

    (3)随隧道间距D增大,山岭隧道衬砌隧道内部位移幅值整体上有减小的趋势,隧道空间位移有关于隧道中心线对称的趋势,两隧道间相互作用减弱。当入射频率较低时,山岭隧道水平位移最大值出现在右隧道顶部,竖向位移最大值出现在隧道右端(两隧道外侧)。

    (4)随隧道间距D减小,隧道内部衬砌环向应力有增大趋势。随入射角度\alpha 增大,两隧道衬砌内部环向应力均有增加趋势,靠近地震波入射一侧隧道应力幅值整体上小于远侧隧道,实际工程设计时应注意远离波入射一侧隧道强度破坏。

  • 图  1   子结构提取示意图

    Figure  1.   Schematic diagram of substructure extraction

    图  2   试件尺寸和配筋详图

    Figure  2.   Dimensions and reinforcement details of specimens

    图  3   试验装置及测量仪器布置图

    Figure  3.   Test setup and instrumentation layout

    图  4   承载力-位移曲线

    Figure  4.   Load resistance-displacement curves

    图  5   R1破坏模式

    Figure  5.   Failure mode of R1

    图  6   PP破坏模式

    Figure  6.   Failure mode of PP

    图  7   PS破坏模式

    Figure  7.   Failure mode of PS

    图  8   试件水平反力

    Figure  8.   Horizontal reaction forces of specimens

    图  9   PP几何模型

    Figure  9.   Geometrical model of PP

    图  10   PCSL-0.75[34]几何尺寸和配筋详图

    Figure  10.   Dimensions and reinforcement details of PCSL-0.75[34]

    图  11   竖向抗力对比

    Figure  11.   Comparison of vertical load resistance

    图  12   数值模型破坏模态

    Figure  12.   Failure modes of numerical models

    图  13   PCSL-0.75竖向抗力对比

    Figure  13.   Comparison of vertical load resistance of PCSL-0.75

    图  14   PCSL-0.75破坏模式对比

    Figure  14.   Comparison of failure model of PCSL-0.75

    图  15   预应力筋的影响

    Figure  15.   Effect of prestressed tendons

    图  16   钢绞线黏结应力的影响

    Figure  16.   Effect of bonded stress of strands

    图  17   张拉控制应力的影响

    Figure  17.   Effect of tension control stress

    图  18   轴压比的影响

    Figure  18.   Effect of axial compression ratio

    表  1   试件详情

    Table  1   Details of specimens

    编号失效柱位置轴压比钢绞线
    布置形式
    张拉控制应力钢绞线
    配筋率/(%)
    R1次边缘柱0.25
    PP次边缘柱0.25抛物线0.65fpu0.3
    PS次边缘柱0.25直线0.65 fpu0.6
    注: fpu为预应力筋极限强度标准值(1860 MPa)。
    下载: 导出CSV

    表  2   钢筋和钢绞线力学特性

    Table  2   Properties of reinforcements and tendon

    钢筋及
    钢绞线
    直径/
    mm
    屈服强度/
    MPa
    极限强度/
    MPa
    弹性模量/
    GPa
    伸长率/
    (%)
    6636848516220.1
    121246259617114.7
    161646660418217.0
    S12.712.7164919702136.3
    注: S为7股钢丝钢绞线,其抗拉强度标准值为1860 MPa。
    下载: 导出CSV

    表  3   主要试验结果

    Table  3   Summary of critical test results

    编号位移/mm荷载/kNMHFL/kNMHFR/kN
    FPLULFPLUL
    R1766655174−81/153−69/143
    PP666006186−87/166−78/152
    PS484997799−61/162−61/122
    注:FPL和UL分别表示第一峰值承载力和极限承载力;MHFL和MHFR分别为试件的左侧和右侧最大水平反力,负值为压反力、正值为拉反力。
    下载: 导出CSV
  • [1] 陆新征, 李易, 叶列平, 等. 钢筋混凝土框架结构抗连续倒塌设计方法的研究[J]. 工程力学, 2008, 25(增刊 2): 150 − 157.

    LU Xinzheng, LI Yi, YE Lieping, et al. Study on design method to resist progressive collapse for reinforced concrete frames [J]. Engineering Mechanics, 2008, 25(Suppl 2): 150 − 157. (in Chinese)

    [2]

    UFC 4-023-03, Design of structures to resist progressive collapse [S]. Washington D C: Department of Defense, USA, 2005.

    [3]

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

    [4]

    CEN. EN 1991-1-7: Eurocode 1: Actions on Structures: Part 1-7: Accidental Actions [S]. Brussels: European Committee for Standardization, 2006.

    [5]

    GSA 2005, Progressive collapse analysis and design guidelines for new federal office buildings and major modernization project [S]. Washington D C: the U. S. General Services Administration, 2005.

    [6]

    YI W J, HE Q F, XIAO Y, et al. Experimental study on progressive collapse-resistant behavior of reinforced concrete frame structures [J]. ACI Structural Journal, 2008, 105(4): 433 − 439.

    [7] 肖宇哲, 李易, 陆新征, 等. 混凝土梁柱子结构连续倒塌动力效应的试验研究[J]. 工程力学, 2019, 36(5): 44 − 52. doi: 10.6052/j.issn.1000-4750.2018.04.0189

    XIAO Yuzhe, LI Yi, LU Xinzheng, et al. Experimental study on the dynamic effects in progressive collapse of beam-column concrete structures [J]. Engineering Mechanics, 2019, 36(5): 44 − 52. (in Chinese) doi: 10.6052/j.issn.1000-4750.2018.04.0189

    [8] 陈泽帆, 林楷奇, 陆新征, 等. RC框架梁柱子结构抗连续倒塌性能不确定性分析[J]. 工程力学, 2021, 38(6): 72 − 80. doi: 10.6052/j.issn.1000-4750.2020.07.0464

    CHEN Zefan, LIN Kaiqi, LU Xinzheng, et al. Uncertainty. analysis on progressive collapse resistance of RC beam-column substructures [J]. Engineering Mechanics, 2021, 38(6): 72 − 80. (in Chinese) doi: 10.6052/j.issn.1000-4750.2020.07.0464

    [9]

    STINGER S M, ORTON S L. Experimental evaluation of disproportionate collapse resistance in reinforced concrete frames [J]. ACI Structural Journal, 2013, 110(3): 521 − 530.

    [10]

    YU J, TAN K H. Structural behavior of RC beam-column subassemblages under a middle column removal scenario [J]. Journal of Structural Engineering, 2012, 139(2): 233 − 250.

    [11]

    TSAI M H, CHANG Y T. Collapse-resistant performance of RC beam–column subassemblages with varied section depth and stirrup spacing [J]. The Structural Design of Tall and Special Buildings, 2015, 24(8): 555 − 570. doi: 10.1002/tal.1199

    [12]

    YU J, TAN K H. Experimental and numerical investigation on progressive collapse resistance of reinforced concrete beam column sub-assemblages [J]. Engineering Structures, 2013, 55(1): 90 − 106.

    [13]

    LONG X, WANG S, HUANG X J, et al. Progressive collapse resistance of exterior reinforced concrete frames and simplified method for catenary action [J]. Engineering Structures, 2021, 249: 113316. doi: 10.1016/j.engstruct.2021.113316

    [14]

    QIAN K, LI B. Experimental and analytical assessment on RC interior beam-column subassemblages for progressive collapse [J]. Journal of Performance of Constructed Facilities, 2012, 26(5): 576 − 589. doi: 10.1061/(ASCE)CF.1943-5509.0000284

    [15]

    QIAN K, LI B. Effects of masonry infill wall on the performance of RC frames to resist progressive collapse [J]. Journal of Structural Engineering, 2017, 143(9): 04017118. doi: 10.1061/(ASCE)ST.1943-541X.0001860

    [16] 李易, 陆新征, 任爱珠, 等. 某八层混凝土框架结构火灾连续倒塌模拟[J]. 工程力学, 2011, 28(增刊 1): 53 − 59.

    LI Yi, LU Xinzheng, REN Aizhu, et al. Simulation for fire-induced progressive collapse of an 8-storey RC frame structure [J]. Engineering Mechanics, 2011, 28(Suppl 1): 53 − 59. (in Chinese)

    [17] 李易, 陆新征, 叶列平, 等. 混凝土框架结构火灾连续倒塌数值分析模型[J]. 工程力学, 2012, 29(4): 96 − 103. doi: 10.6052/j.issn.1000-4750.2010.06.0397

    LI Yi, LU Xinzheng, YE Lieping, et al. Numerical models of fire induced progressive collapse analysis for reinforced concrete frame structures [J]. Engineering Mechanics, 2012, 29(4): 96 − 103. (in Chinese) doi: 10.6052/j.issn.1000-4750.2010.06.0397

    [18]

    LI Z, LIU Y, HUO J, et al. Experimental assessment of fire-exposed RC beam-column connections with varying reinforcement development lengths subjected to column removal [J]. Fire Safety Journal, 2018, 99: 38 − 48. doi: 10.1016/j.firesaf.2018.06.003

    [19]

    LI Z, LIU Y, HUO J, et al. Experimental and analytical assessment of RC joints with varying reinforcement detailing under push-down loading before and after fires [J]. Engineering Structures, 2019, 189: 550 − 564. doi: 10.1016/j.engstruct.2019.03.084

    [20] 钱 凯, 黄志强, 翁运昊, 等. 锈蚀RC框架结构连续倒塌抗力机制研究 [J/OL]. 建筑结构学报. DOI: 10.14006/j.jzjgxb.2021.0010.

    QIAN Kai, HUANG Zhiqiang, WENG Yunhao, et al. Study on load carrying mechanism of corroded RC Frame structures against progressive collapse [J/OL]. Journal of Building Structures. DOI: 10.14006/j.jzjgxb.2021.0010. (in Chinese)

    [21]

    YU J, TAN K H. Special Detailing techniques to improve structural resistance against progressive collapse [J]. Journal of Structure Engineering, 2014, 140: 04013077. doi: 10.1061/(ASCE)ST.1943-541X.0000886

    [22] 李爽, 赵颖, 翟长海, 等. 节点对RC框架结构抗连续倒塌能力影响研究[J]. 工程力学, 2012, 29(12): 80 − 87. doi: 10.6052/j.issn.1000-4750.2011.04.0242

    LI Shuang, ZHAO Ying, ZHAI Changhai, et al. The influences of joint on the progressive collapse-resisting performance of RC frames [J]. Engineering Mechanics, 2012, 29(12): 80 − 87. (in Chinese) doi: 10.6052/j.issn.1000-4750.2011.04.0242

    [23]

    ORTON S, JIRSA J, BAYRAK O. Carbon fiber-reinforced polymer for continuity in existing reinforced concrete buildings vulnerable to collapse [J]. ACI Structure, 2009, 10(6): 608 − 616.

    [24]

    QIAN K, WENG Y H, LI B. Improving behavior of reinforced concrete frames to resist progressive collapse through steel bracings [J]. Journal of Structural Engineering, 2019, 145(2): 04018248. doi: 10.1061/(ASCE)ST.1943-541X.0002263

    [25]

    LI S, SHAN S D, ZHAI C H, et al. Experimental and numerical study on progressive collapse process of RC frames with full-height infill walls [J]. Engineering Failure Analysis, 2016, 59: 57 − 68. doi: 10.1016/j.engfailanal.2015.11.020

    [26]

    SHAN S D, LI S, XU S Y, et al. Experimental study on the progressive collapse performance of RC frames with infill walls [J]. Engineering Structures, 2016, 111: 80 − 92. doi: 10.1016/j.engstruct.2015.12.010

    [27]

    QIAN K, LAN D Q, FU F, et al. Effects of infilled wall opening on load resisting capacity of RC frames to mitigate progressive collapse risk [J]. Engineering Structures, 2020, 223: 111196. doi: 10.1016/j.engstruct.2020.111196

    [28]

    LU X Z, LIN K Q, LI Y, et al. Experimental investigation of RC beam-slab substructures against progressive collapse subjected to an edge-column-removal scenario [J]. Engineering Structure, 2017, 149: 91 − 103. doi: 10.1016/j.engstruct.2016.07.039

    [29] 钱凯, 李治, 翁运昊, 等. 钢筋混凝土梁-板子结构抗连续性倒塌性能研究[J]. 工程力学, 2019, 36(6): 239 − 247. doi: 10.6052/j.issn.1000-4750.2018.05.0297

    QIAN Kai, LI Zhi, WENG Yunhao, et al. Behavior of RC beam-slab substructures to resist progressive collapse [J]. Engineering Mechanics, 2019, 36(6): 239 − 247. (in Chinese) doi: 10.6052/j.issn.1000-4750.2018.05.0297

    [30]

    XIAO Y, KUNNATH S, LI F W, et al. Collapse test of three-story half-scale reinforced concrete frame building [J]. ACI Structural Journal, 2015, 112(4): 429 − 438.

    [31] 周云, 陈太平, 胡翔, 等. 考虑周边结构约束影响的RC框架结构防连续倒塌性能研究[J]. 工程力学, 2019, 36(1): 216 − 226. doi: 10.6052/j.issn.1000-4750.2017.11.0887

    ZHOU Yun, CHEN Taiping, HU Xiang, et al. progressive collapse resistance of RC frame structures considering surrounding structural constraints [J]. Engineering Mechanics, 2019, 36(1): 216 − 226. (in Chinese) doi: 10.6052/j.issn.1000-4750.2017.11.0887

    [32] 张有佳, 王冲, 杜轲. 不等跨布置RC空间梁-板结构竖向倒塌能力数值研究[J]. 工程力学, 2022, 39(1): 139 − 150. doi: 10.6052/j.issn.1000-4750.2020.12.0891

    ZHANG Youjia, WANG Chong, DU Ke. Numerical study on the vertical collapse capacity of RC spatial beam-slab structures with unequal spans [J]. Engineering Mechanics, 2022, 39(1): 139 − 150. (in Chinese) doi: 10.6052/j.issn.1000-4750.2020.12.0891

    [33]

    QIAN K, LIU Y, YANG T, et al. Progressive collapse resistance of posttensioned concrete beam-column subassemblages with unbonded posttensioning strands [J]. Journal of Structural Engineering, 2018, 144: 4017182. doi: 10.1061/(ASCE)ST.1943-541X.0001940

    [34]

    QIAN K, ZHANG X D, FU F, et al. Progressive collapse-resisting mechanisms of planar prestressed concrete frame [J]. ACI Structural Journal, 2019, 116(4): 77 − 90.

    [35] 余洋, 李治, 肖龙山, 等. 边柱失效后预应力拼接连接装配式结构抗连续倒塌机理研究[J]. 工程力学, 2021, 38(4): 159 − 168. doi: 10.6052/j.issn.1000-4750.2020.06.0366

    YU Yang, LI Zhi, XIAO Longshan, et al. Load resisting mechanism of precast structure under exterior column failure [J]. Engineering Mechanics, 2021, 38(4): 159 − 168. (in Chinese) doi: 10.6052/j.issn.1000-4750.2020.06.0366

    [36] 杨涛, 陈万庆, 郝天之. 中柱失效后预应力混凝土框架子结构动力倒塌性能试验研究[J]. 振动与冲击, 2020, 39(20): 18 − 23.

    YANG Tao, CHEN Wanqing, HAO Tianzhi. Experimental study on dynamic collapse performance of prestressed RC frame substructures after middle column demolition [J]. Journal of Vibration and Shock, 2020, 39(20): 18 − 23. (in Chinese)

    [37] 金浏, 兰冬璆, 邓小芳, 等. 无黏结预应力梁-柱子结构抗连续倒塌机理研究 [J/OL]. 建筑结构学报: 1 − 12 [2022-04-03]. DOI: 10.14006/j.jzjgxb.2021.0572.

    JIN Liu, LAN Dongqiu, DENG Xiaofang et al. Study on progressive collapse mechanisms of unbonded prestressed beam-column sub-assemblages [J/OL]. Journal of Building Structures: 1 − 12 [2022-04-03]. DOI: 10.14006/j.jzjgxb. 2021.0572.

    [38] GB 50010−2010, 混凝土结构设计规范 [S]. 北京: 中国建筑工业出版社, 2015.

    GB 50010−2010, Code for design of concrete structures [S]. Beijing: China Architecture & Building Press, 2015.

    [39] JGJ 369−2016, 预应力混凝土结构设计规范 [S]. 北京: 中国建筑工业出版社, 2016.

    JGJ 369−2016, Code for design of prestressed concrete structures [S]. Beijing: China Architecture &Building Press, 2016. (in Chinese)

    [40]

    ACI Committee 318, Building code requirements for structural concrete (ACI 318-14) and Commentary (318R-14) [S]. Farmington Hills, MI: American Concrete Institute, 2014: 433.

  • 期刊类型引用(13)

    1. 彭志明,陈忠范,徐志峰,黄宇浩,罗漪. 不同截面冷弯型钢-泡沫混凝土组合墙体抗震性能. 东南大学学报(自然科学版). 2023(05): 783-793 . 百度学术
    2. 吴函恒,隋璐,聂少锋,周天华,袁涛涛. 填充石膏基轻质材料的冷弯型钢复合墙体受剪承载力分析. 工程力学. 2022(04): 177-186 . 本站查看
    3. 胡习兵,张正武,黄纵威,毛健宇,李清山,袁智深. 轻钢龙骨泡沫混凝土轻质墙板抗弯性能研究. 建筑钢结构进展. 2022(09): 67-75 . 百度学术
    4. 曹万林,李浩,乔崎云,董宏英,边瑾靓. 较小剪跨比轻钢-尾砂微晶发泡板组合墙抗震性能试验. 地震工程与工程振动. 2020(01): 80-87 . 百度学术
    5. 王欢迎. 铁路客运专线应用KF装配式四电房屋智能建造化建设技术. 低碳世界. 2020(04): 148+150 . 百度学术
    6. 张铮,江忠画,陈笃海. 冷弯薄壁型钢灌浆墙体抗震性能研究进展. 福建工程学院学报. 2020(03): 205-210 . 百度学术
    7. 田稳苓,温晓东,彭佳斌,徐丽丽,李子祥. 新型泡沫混凝土轻钢龙骨复合墙体抗剪承载力计算方法研究. 工程力学. 2019(09): 143-153 . 本站查看
    8. 李长青,于旭伟,张轮轮. 铁路客运专线装配式KF轻钢结构四电房屋智能建造应用. 建筑技术开发. 2019(18): 3-5 . 百度学术
    9. 徐志峰,陈忠范,朱松松,刘吉,殷之祺. 秸秆板轻钢高强泡沫混凝土剪力墙轴心受压性能研究. 工程力学. 2018(07): 219-231 . 本站查看
    10. 王永鹏,杜荣强,王路喜. 钢塑复合板抗压承载力试验及有限元模拟. 中国科技论文. 2018(01): 94-97 . 百度学术
    11. 潘秀珍,赵登云,刘辉,张鹏,田建勃. 轻钢住宅结构整体抗震性能的影响研究. 西安理工大学学报. 2018(03): 278-286 . 百度学术
    12. 袁泉,吕东鑫,杨逸,聂屿,霍学朝. 新型轻钢灌浆墙节能结构体系. 墙材革新与建筑节能. 2017(02): 47-50 . 百度学术
    13. 贺思维,曲哲,周惠蒙,戴君武,王多智. 非结构构件抗震性能试验方法综述. 土木工程学报. 2017(09): 16-27 . 百度学术

    其他类型引用(25)

图(18)  /  表(3)
计量
  • 文章访问数:  681
  • HTML全文浏览量:  147
  • PDF下载量:  91
  • 被引次数: 38
出版历程
  • 收稿日期:  2022-01-09
  • 修回日期:  2022-04-15
  • 网络出版日期:  2022-04-22
  • 刊出日期:  2023-10-09

目录

/

返回文章
返回