Loading [MathJax]/jax/output/SVG/jax.js

高宽比为4的矩形断面脉动风压的非高斯特性及其极值研究

阮礼君, 阮冠康, 高梓涵, 郑超越, 刘小兵

阮礼君, 阮冠康, 高梓涵, 郑超越, 刘小兵. 高宽比为4的矩形断面脉动风压的非高斯特性及其极值研究[J]. 工程力学, 2025, 42(S): 239-243. DOI: 10.6052/j.issn.1000-4750.2024.06.S051
引用本文: 阮礼君, 阮冠康, 高梓涵, 郑超越, 刘小兵. 高宽比为4的矩形断面脉动风压的非高斯特性及其极值研究[J]. 工程力学, 2025, 42(S): 239-243. DOI: 10.6052/j.issn.1000-4750.2024.06.S051
RUAN Li-jun, RUAN Guan-kang, GAO Zi-han, ZHENG Chao-yue, LIU Xiao-bing. STUDY ON NON-GAUSSIAN CHARACTERISTICS AND EXTREME VALUES OF FLUCTUATING WIND PRESSURE ON RECTANGULAR SECTION WITH ASPECT RATIO OF 4[J]. Engineering Mechanics, 2025, 42(S): 239-243. DOI: 10.6052/j.issn.1000-4750.2024.06.S051
Citation: RUAN Li-jun, RUAN Guan-kang, GAO Zi-han, ZHENG Chao-yue, LIU Xiao-bing. STUDY ON NON-GAUSSIAN CHARACTERISTICS AND EXTREME VALUES OF FLUCTUATING WIND PRESSURE ON RECTANGULAR SECTION WITH ASPECT RATIO OF 4[J]. Engineering Mechanics, 2025, 42(S): 239-243. DOI: 10.6052/j.issn.1000-4750.2024.06.S051

高宽比为4的矩形断面脉动风压的非高斯特性及其极值研究

基金项目: 

国家自然科学基金项目(52078313,52008273);中央引导地方科技发展资金项目(236Z5407G)

详细信息
    作者简介:

    阮礼君(1999−),男,安徽人,硕士生,主要从事结构风荷载、风致振动与控制研究(E-mail: rlj0527@163.com)

    阮冠康(2001−),男,广东人,硕士生,主要从事结构风荷载、风致振动与控制研究(E-mail: guankangruan@163.com)

    高梓涵(2000−),男,河北人,硕士生,主要从事结构风荷载、风致振动与控制研究(E-mail: 1799056920@qq.com)

    郑超越(2001−),男,河北人,硕士生,主要从事结构风荷载、风致振动与控制研究(E-mail: zheng_chao_yue@163.com)

    通讯作者:

    刘小兵(1982−),男,湖南人,教授,博士,博导,主要从事长高大工程结构的气动特性及优化措施研究(E-mail: x_b_liu@126.com)

  • 中图分类号: TU312+.1

STUDY ON NON-GAUSSIAN CHARACTERISTICS AND EXTREME VALUES OF FLUCTUATING WIND PRESSURE ON RECTANGULAR SECTION WITH ASPECT RATIO OF 4

  • 摘要:

    通过刚性模型风洞测压试验,研究了高宽比为4的矩形断面表面脉动风压概率密度分布及其非高斯特性,并对高斯分布与非高斯分布的区域进行了划分。在此基础上,采用改进的Hermite法计算了非高斯区域的峰值因子及脉动风压极值。结果表明:高宽比为4的矩形断面在侧风面和背风面测点完全位于非高斯区域,其测点的概率密度分布明显偏离高斯分布,呈现出显著的非高斯特性。此外,规范建议峰值因子取值和基于高斯过程的峰值因子均明显低于通过改进的Hermite法确定的峰值因子,建议采用改进的Hermite法来计算非高斯区域的脉动风压极值。

    Abstract:

    This study investigates the probability density distribution and non-Gaussian characteristics of fluctuating wind pressure on a rectangular section with an aspect ratio of 4. Rigid model pressure measurement wind tunnel tests were conducted to analyze the regions of Gaussian and non-Gaussian characteristics. The modified Hermite method was then used to calculate the peak factor and extreme values of the non-Gaussian region. Results indicate that pressure taps on the crosswind and leeward sides of the rectangular section fall entirely within the non-Gaussian region, displaying significant deviations and bulges in the probability density distribution. Furthermore, both the peak factor suggested by the code and the peak factor based on Gaussian processes are significantly lower than the peak factor determined through the improved Hermite method, which is recommended to be used for calculating the fluctuating wind pressure extremes in the non-Gaussian region.

  • 在结构风工程领域,通常认为结构表面的风压信号服从高斯分布。DAVENPORT [1]基于高斯假定提出的传统峰值因子法广泛应用于工程设计之中。我国现行《建筑结构荷载规范》[2](以下简称“规范”)建议围护结构设计的峰值因子取值为2.5。然而,由于结构表面存在剪切层的分离、再附着和旋涡的形成与脱落等现象的发生,使得局部区域的脉动风压表现出明显的非高斯特性[3]。与高斯风荷载相比,非高斯风荷载导致结构的疲劳受损的研究并不少见,比如大跨空间结构[4]、大型冷却塔[56]以及高层建筑[78]等。这些研究解释了建筑表面风压非高斯特性的产生原因,并分析其非高斯特性对风压极值的影响。因此,研究脉动风压的非高斯特性及其极值仍然是风工程的重要课题。

    在实际工程中,矩形断面是一种常见的结构形式,其气动特性虽然经过大量研究,但由于其复杂的流动特征,其非高斯特性及其极值的研究仍然引起广泛关注。韩宁等[9]发现高层方形建筑的迎风面以高斯区域为主,而受分离流和尾流作用的立面主要以非高斯区域为主。而对于矩形高层建筑,庄翔等[10]发现当长边迎风时,非高斯区域主要出现在迎风角点附近,当短边迎风时出现在被风角点及尾流区附近。YANG等[11]发现矩形柱分离区的风压非高斯特性随湍流强度的增强而变得更为显著。

    WINTERSTEIN[12]将Hermite矩模型引入用来实现目标非高斯随机变量到高斯随机变量的转换[13],以期获得精度更高的风压极值。通过改变传统的Hermite法的参数,建立了更加准确的计算模型(改进的Hermite法)并得到广泛的应用[1416]。综上所述,本文采用改进的Hermite法来计算高宽比为4的矩形断面表面脉动风压极值,将为该类型围护结构的抗风设计提供了重要参考。

    试验在石家庄铁道大学风工程研究中心STU-1风洞实验室低速试验段中进行,试验段的尺寸长24 m,高3 m,宽4.38 m,最大风速可达30 m/s。图1展示了试验模型及风向角示意图。试验模型采用ABS板材制作,模型尺寸长(L)×宽(B)×高(H)为2000 mm×80 mm×320 mm。在模型中间位置沿着周向共布置了88个测压点,而在矩形柱的4个角部测压点则进行了加密处理。风压数据通过电子扫描阀采集,电子压力扫描阀的采样频率为330 Hz,采样时间为30 s,每个测点的采样数量N=9900。

    图  1  试验模型及风向角示意图 /mm
    Figure  1.  Test model and wind direction

    试验在均匀流场下进行,试验风速为10 m/s。来流方向与展向方向垂直,试验阻塞度约为5%,不需要对试验结果进行阻塞度修正。

    对于非高斯分布,采用偏度(γ3)和峰度(γ4)表示,其定义为:

    γ3=1NNi=1[(CPiμσ)3] (1)
    γ4=1NNi=1[(CPiμσ)4] (2)

    式中:CPi为测点i的风压系数;μ为该测点的平均风压系数;σ为该测点风压系数的均方根;N为测点风压的试验采样数量。

    为了确定脉动风压是否为非高斯分布,需要考虑偏度和峰度。偏度为正时表示概率密度分布右偏,为负则为左偏,因此偏度需取绝对值作为非高斯区域的判断标准。不同研究对于非高斯信号的判别条件存在差异。比如,孙瑛等[17]将|γ3|>0.2且γ4>3.7作为大跨平屋盖非高斯区域的划分标准。CHEN等[18]将|γ3|>0.2且γ4>3.5适用于大跨空间网壳结构。庄翔等[10]和DU等[19]认为|γ3|>0.5且γ4>3.5适用于矩形截面表面风压非高斯区域的划分条件。因此,结合试验结果及过去的研究,将|γ3|>0.5且γ4>3.5作为非高斯区域的划分条件。图2给出了脉动风压偏度与峰度系数的分布,可以看到迎风面测点大多符合高斯分布,而侧风面和背风面测点完全处于非高斯区域之中。这是由于迎风角点发生剪切层分离,在侧风面形成大尺度涡,以及在背风面可能出现的卡门涡的作用[2022]

    图  2  脉动风压的偏度与峰度
    Figure  2.  Skewness and kurtosis of the fluctuating wind pressure

    图3展示了3个典型测点的风压系数时程和标准化风压系数概率密度分布,其中黑线代表相应的平均风压系数。对于迎风面中部的测点P15,从图3(a)可以观察到其风压系数波动较小,且随着时间的推移其风压系数对称性良好。从图3(d)可以看出,该测点概率密度接近高斯分布。而在背风面测点P59,在图3(b)中可以看出其风压系数波动显著,表现出明显的不对称性。类似的是,对于侧风面测点P88,从图3(c)可以看到其风压系数波动也十分显著,表现出明显的不对称性。因此,背风面测点P59和侧风面测点P88的脉动风压概率密度分布表现出更为明显的偏斜和凸起。此外,由于存在许多间歇性的负脉冲信号,这些测点受到较大的风吸力,导致概率密度分布在负向一侧出现长尾,表现出显著的非高斯特性,如图3(d)所示。这种现象在侧风面分离区比尾流区更加明显,因此测点P88的概率密度分布较测点P59的更偏离正态分布,展现出更显著的非高斯特性。

    图  3  典型测点的风压系数时程及其概率密度分布
    Figure  3.  Time histories of wind pressure coefficients and their probability density distributions at typical measurement points

    为了准确计算矩形断面的表面脉动风压极值,采用基于峰度和偏度将非高斯过程转化成高斯过程的改进Hermite法[4],其峰值因子计算公式如下:

    g=κ{(β+γβ)+h3[β2+(2γ1)+1.98β2]+h4[β3+3β(γ1)+3β(π26γ+γ2)+5.44β3]} (3)
    h3=p1γ3+p2γ3γ4+p3γ33+p4γ3γ24+p5γ33γ4+p6γ3γ34+p7γ53+p8γ33γ24+p9γ3γ44 (4)
    h4=p1+p2γ4+p3γ23+p4γ24+p5γ23γ4+p6γ34+p7γ43+p8γ23γ24+p9γ44+p10γ43γ4+p11γ24γ34+p12γ54 (5)
    κ=11+2h23+6h24 (6)

    式中:β=2ln(v0T)v0为零穿越率,T为时距,v0=m2/m0mi为第i阶谱矩,mi=0n2Sy(n)dnSy(n)为单边功率谱密度,n为频率;γ为欧拉常数,取0.5772;h3h4中的参数pi由YANG等[23]给出。

    因此,脉动风压的极值可由式(7)计算,当计算正极值时取正号,反之,取负号。

    Ye=μ±gσ (7)

    γ3=0和γ4=3,即高斯过程时,Hermite级数中的参数h3h4为0,此时即为高斯过程的峰值因子计算方法[1],式(3)变为:

    g,G=β+γβ=2ln(v0T)+γ2ln(v0T) (8)

    因为矩形断面大部分属于非高斯区域,所以在计算风压极值时需要考虑非高斯特性。由于围护结构受负风压极值的影响显著,选取了非高斯区域的典型测点的负风压极值进行研究,如表1所示。结果表明:非高斯区域的脉动风压峰值因子与位置密切相关,且均超过规范建议值的两倍。采用非高斯方法计算的风压极值(μ-)远超高斯方法计算的风压极值(μ-g,Gσ),基于高斯分布计算的风压极值约为非高斯方法计算的80%。而与采用规范建议峰值因子取值计算的风压极值(μ-2.5σ)差异更显著,约为非高斯方法计算风压极值的70%。换句话说,对于这种矩形断面结构,忽略非高斯特性将低估20%~30%的表面脉动风压极值。可见,若统一采用高斯区峰值因子进行围护结构抗风设计会明显低估了风压极值,其结果将偏于不安全。因此,在矩形断面围护结构抗风设计时,侧风面和背风面的脉动风压峰值因子应按照改进Hermite法进行计算。

    表  1  非高斯区域典型测点风压极值对比
    Table  1.  The comparison of extreme wind pressure for typical taps in non-Gaussian region
    测点 均值μ 标准差σ 峰值因子g,G 非高斯峰值因子g i) μ-2.5σ ii) μ- g,Gσ iii) μ- i)/iii) ii)/iii)
    45 −1.725 0.543 3.296 5.367 −3.083 −3.515 −4.639 0.664 0.758
    59 −1.672 0.407 3.362 5.259 −2.689 −3.040 −3.812 0.706 0.798
    74 −1.770 0.581 3.300 5.347 −3.223 −3.687 −4.877 0.661 0.756
    81 −1.719 0.543 3.303 5.076 −3.077 −3.513 −4.475 0.687 0.785
    88 −1.668 0.539 3.302 5.181 −3.016 −3.448 −4.461 0.676 0.773
    下载: 导出CSV 
    | 显示表格

    由式(7)可知,仅凭借峰值因子的大小判断非高斯特性和极值风压是不充分的。因此,我们绘制了高宽比为4的矩形断面的极值风压系数,可为实际工程设计提供参考,如图4所示。其中需要说明的是:依照图2的非高斯区域的划分方式,将非高斯区域的极值风压系数的峰值因子按式(3)计算,高斯区域的极值风压系数的峰值因子按照式(8)计算。从图4中可以看到,在背风面和侧风面均表现出较大的极值风压系数。其中,在背风面卡门涡形成的位置极值风压系数的绝对值明显要大于其他区域的极值风压系数。而靠近背风面中点的位置的极值风压系数的绝对值会突然减小,这说明卡门涡的形成与脱落对背风面中点的影响程度较弱。

    图  4  极值风压系数分布
    Figure  4.  Distribution of extreme wind pressure coefficients

    通过刚性模型风洞测压试验,研究了高宽比为4的矩形断面在来流垂直作用于展向长度时表面脉动风压的非高斯特性,并运用改进的Hermite法计算脉动风压的极值,得出以下主要结论:

    (1)依据偏度|γ3|>0.5且峰度γ4>3.5的非高斯区域划分标准,发现在高宽比为4时,矩形断面的侧风面和背风面完全处于脉动风压的非高斯区域。其脉动风压概率密度分布与高斯分布相比,表现出显著的偏斜和凸起现象,尤其在负向一侧出现长尾。

    (2)对于高宽比为4的矩形断面,在背风面卡门涡形成的区域的极值风压系数的绝对值相较于其他区域更大。

    (3)忽视脉动风压的非高斯特性,将可能低估20%~30%的风压极值,进而导致结构不安全。因此,建议采用改进的Hermite法来计算非高斯区域的脉动风压极值。

  • 图  1   试验模型及风向角示意图 /mm

    Figure  1.   Test model and wind direction

    图  2   脉动风压的偏度与峰度

    Figure  2.   Skewness and kurtosis of the fluctuating wind pressure

    图  3   典型测点的风压系数时程及其概率密度分布

    Figure  3.   Time histories of wind pressure coefficients and their probability density distributions at typical measurement points

    图  4   极值风压系数分布

    Figure  4.   Distribution of extreme wind pressure coefficients

    表  1   非高斯区域典型测点风压极值对比

    Table  1   The comparison of extreme wind pressure for typical taps in non-Gaussian region

    测点 均值μ 标准差σ 峰值因子g,G 非高斯峰值因子g i) μ-2.5σ ii) μ- g,Gσ iii) μ- i)/iii) ii)/iii)
    45 −1.725 0.543 3.296 5.367 −3.083 −3.515 −4.639 0.664 0.758
    59 −1.672 0.407 3.362 5.259 −2.689 −3.040 −3.812 0.706 0.798
    74 −1.770 0.581 3.300 5.347 −3.223 −3.687 −4.877 0.661 0.756
    81 −1.719 0.543 3.303 5.076 −3.077 −3.513 −4.475 0.687 0.785
    88 −1.668 0.539 3.302 5.181 −3.016 −3.448 −4.461 0.676 0.773
    下载: 导出CSV
  • [1]

    DAVENPORT A G. Note on the distribution of the largest value of a random function with application to gust loading [J]. Proceedings of the Institution of Civil Engineers, 1964, 28(2): 187 − 196. doi: 10.1680/iicep.1964.10112

    [2] GB 50009−2012, 建筑结构荷载规范[S]. 北京: 中国建筑工业出版社, 2012.

    GB 50009−2012, Load code for the design of building structures [S]. Beijing: China Architecture & Building Press, 2012. (in Chinese)

    [3]

    KUMAR K S, STATHOPOULOS T. Wind loads on low building roofs: A stochastic perspective [J]. Journal of Structural Engineering, 2000, 126(8): 944 − 956. doi: 10.1061/(ASCE)0733-9445(2000)126:8(944)

    [4]

    SONG J, XU W, HU G, et al. Non-gaussian properties and their effects on extreme values of wind pressure on the roof of long-span structures [J]. Journal of Wind Engineering and Industrial Aerodynamics, 2019, 184: 106 − 115. doi: 10.1016/j.jweia.2018.11.027

    [5]

    CHENG X X, KE S T, LI P F, et al. External extreme wind pressure distribution for the structural design of cooling towers [J]. Engineering Structures, 2019, 181: 336 − 353. doi: 10.1016/j.engstruct.2018.12.038

    [6] 王浩, 柯世堂. 不同四塔组合形式对特大型冷却塔局部非高斯风压分布影响研究[J]. 工程力学, 2018, 35(8): 162 − 171. doi: 10.6052/j.issn.1000-4750.2017.04.0309

    WANG Hao, KE Shitang. Research on non-gaussian wind pressure of four super-large cooling towers under different layouts [J]. Engineering Mechanics, 2018, 35(8): 162 − 171. (in Chinese) doi: 10.6052/j.issn.1000-4750.2017.04.0309

    [7] 庄翔, 董欣, 郑毅敏, 等. 矩形高层建筑非高斯风压时程峰值因子计算方法[J]. 工程力学, 2017, 34(7): 177 − 185, 223. doi: 10.6052/j.issn.1000-4750.2016.02.0134

    ZHUANG Xiang, DONG Xin, ZHENG Yimin, et al. Peak factor estimation methods of non-gaussian wind pressures on a rectangular high-rise building [J]. Engineering Mechanics, 2017, 34(7): 177 − 185, 223. (in Chinese) doi: 10.6052/j.issn.1000-4750.2016.02.0134

    [8] 张海程, 杜树碧, 李明水, 等. 紊流积分尺度对caarc模型迎风面极值风压特性的影响及修正方法[J]. 工程力学, 2023, 40(10): 11 − 20. doi: 10.6052/j.issn.1000-4750.2022.01.0092

    ZHANG Haicheng, DU Shubi, LI Mingshui, et al. The effect of the turbulence integral scale on the extreme wind pressures on the windward surface of the caarc model and the correction method [J]. Engineering Mechanics, 2023, 40(10): 11 − 20. (in Chinese) doi: 10.6052/j.issn.1000-4750.2022.01.0092

    [9] 韩宁, 顾明. 方形高层建筑风压脉动非高斯特性分析[J]. 同济大学学报(自然科学版), 2012, 40(7): 971 − 976. doi: 10.3969/j.issn.0253-374x.2012.07.002

    HAN Ning, GU Ming. Analysis on non-gaussian features of fluctuating wind pressure on square tall buildings [J]. Journal of Tongji University (Natural Science), 2012, 40(7): 971 − 976. (in Chinese) doi: 10.3969/j.issn.0253-374x.2012.07.002

    [10] 庄翔, 董欣, 丁洁民, 等. 矩形高层建筑表面风压脉动的非高斯特性研究[J]. 建筑结构学报, 2016, 37(增刊1): 13 − 18.

    ZHUANG Xiang, DONG Xin, DING Jiemin, et al. Non-Gaussian features of wind pressure fluctuations on a rectangular high-rise building [J]. Journal of Building Structures, 2016, 37(Suppl 1): 13 − 18. (in Chinese)

    [11]

    YANG X W, YANG Y, LI M S, et al. Effects of free-stream turbulence on non-Gaussian characteristics of fluctuating wind pressures on a 5: 1 rectangular cylinder [J]. Journal of Wind Engineering and Industrial Aerodynamics, 2021, 217: 104759. doi: 10.1016/j.jweia.2021.104759

    [12]

    WINTERSTEIN S R. Non‐normal responses and fatigue damage [J]. Journal of Engineering Mechanics, 1985, 111(10): 1291 − 1295. doi: 10.1061/(ASCE)0733-9399(1985)111:10(1291)

    [13]

    KWON D K, KAREEM A. Peak factors for non-gaussian load effects revisited [J]. Journal of Structural Engineering, 2011, 137(12): 1611 − 1619. doi: 10.1061/(ASCE)ST.1943-541X.0000412

    [14] 李寿科, 李寿英, 陈政清, 等. 围护结构非高斯风压极值估计的改进hermite峰值因子法[J]. 建筑结构学报, 2015, 36(4): 112 − 118.

    LI Shouke, LI Shouying, CHEN Zhengqing, et al. Prediction of non-gaussian wind pressure extreme value on claddings with modified hermite peak factor method [J]. Journal of Building Structures, 2015, 36(4): 112 − 118. (in Chinese)

    [15]

    PENG L L, LIU M, YANG Q S, et al. An analytical formula for gaussian to non-gaussian correlation relationship by moment-based piecewise hermite polynomial model with application in wind engineering [J]. Journal of Wind Engineering and Industrial Aerodynamics, 2020, 198: 104094. doi: 10.1016/j.jweia.2020.104094

    [16]

    WEI S T, ZHAO C, SUN Q. Investigating the non-gaussian property and its influence on extreme wind pressures on the long-span cylindrical roof [J]. Applied Sciences, 2023, 13(13): 13137691. doi: 10.3390/app13137691

    [17] 孙瑛, 武岳, 林志兴, 等. 大跨屋盖结构风压脉动的非高斯特性[J]. 土木工程学报, 2007, 40(4): 1 − 5, 12. doi: 10.3321/j.issn:1000-131X.2007.04.001

    SUN Ying, WU Yue, LIN Zhixing, et al. Non-Gaussian features of fluctuating wind pressures on long span roofs [J]. China Civil Engineering Journal, 2007, 40(4): 1 − 5, 12. (in Chinese) doi: 10.3321/j.issn:1000-131X.2007.04.001

    [18]

    CHEN G, HUANG Y Q, WU P, et al. Characteristics of extreme wind pressure on the open prefabricated spatial grid structure of evergrande stadium [J]. Buildings, 2023, 13(1): 13010025. doi: 10.3390/buildings13010025

    [19]

    DU X Q, XU Q, DONG H T, et al. Physical mechanisms behind the extreme wind pressures on two tandem square cylinders [J]. Journal of Wind Engineering and Industrial Aerodynamics, 2022, 231: 105249. doi: 10.1016/j.jweia.2022.105249

    [20] 杜晓庆, 靳晓雨, 许汉林, 等. 均匀来流下方柱表面风压非高斯特性的流场机理[J]. 同济大学学报(自然科学版), 2020, 48(2): 191 − 199. doi: 10.11908/j.issn.0253-374x.19176

    DU Xiaoqing, JIN Xiaoyu, XU Hanlin, et al. Uniform flow mechanism of non-gaussian features of pressure on square cylinder [J]. Journal of Tongji University (Natural Science), 2020, 48(2): 191 − 199. (in Chinese) doi: 10.11908/j.issn.0253-374x.19176

    [21] 杨雄伟, 周强, 李明水, 等. 复杂曲面屋盖脉动风压的非高斯特性及峰值因子研究[J]. 振动与冲击, 2021, 40(10): 315 − 322.

    YANG Xiongwei, ZHOU Qiang, LI Mingshui, et al. Non-Gaussian features and peak factors of fluctuating wind pressures on a complex curved roof [J]. Journal of Vibration and Shock, 2021, 40(10): 315 − 322. (in Chinese)

    [22]

    HUANG D M, CHEN G L. Wind pressure distribution, non-gaussian characteristics, and conical vortices of a large-span roof with variable arc angles [J]. Journal of Building Engineering, 2022, 57: 104929. doi: 10.1016/j.jobe.2022.104929

    [23]

    YANG L P, GURLEY K R, PREVATT D O. Probabilistic modeling of wind pressure on low-rise buildings [J]. Journal of Wind Engineering and Industrial Aerodynamics, 2013, 114: 18 − 26. doi: 10.1016/j.jweia.2012.12.014

图(4)  /  表(1)
计量
  • 文章访问数:  2
  • HTML全文浏览量:  0
  • PDF下载量:  1
  • 被引次数: 0
出版历程
  • 收稿日期:  2024-05-31
  • 修回日期:  2025-01-16
  • 网络出版日期:  2025-03-16
  • 刊出日期:  2025-06-24

目录

/

返回文章
返回