Effect of Saturated Complex Site on the Seismic Response of Continuous Long-span Beam Bridges
-
摘要: 我国沿海地区大量大跨桥梁位于饱和复杂场地，此类场地通常存在土饱和特性、起伏地形、上覆水层，其对地震波的散射效应将导致地震动空间变化，然而目前缺乏针对此类场地条件的实测地震记录。针对此问题，本文建立了适用于饱和复杂场地大跨桥梁抗震分析的多点地震动模拟方法。首先，基于边界元法求解饱和复杂场地地震动传递函数；然后采用谱表示法生成空间变化人工地震动加速度时程；最后，以五跨连续梁桥为例，开展大跨桥梁地震响应分析，研究局部场地效应、动水压力效应对大跨连续梁桥地震响应的影响规律。结果表明，饱和复杂场地地震动加速度峰值表现出明显的空间变化特性，多点输入可导致饱和场地中地震动峰值加速度放大2～3倍，进而引起墩顶位移较一致地震作用结果增大95%～110%；动水压力可引起墩顶位移、墩底内力增大34%、29%。饱和复杂场地大跨桥梁抗震分析时忽略场地条件对地震动、桥梁地震响应的影响偏于不安全。Abstract: In China, many long-span bridges are situated across saturated complex sites characterized by soil saturation, undulating topography and superimposed water, in which the propagation of seismic waves induces spatial variations in ground motions due to scattering effect. However, the seismic records under such specific site are lack. Addressing this concern, we establish a multi-point ground motion simulation approach to analyze the seismic response of long-span bridges in saturated sites. Initially, the boundary element method is employed to solve the ground motion transfer function of saturated sites. Subsequently, the spectral representation technique is used to synthesize spatially varying ground motion acceleration time-histories. Conclusively, a continuous five-span beam bridge is set as an example, and the seismic response is investigated to study the impact of local site effects and hydrodynamic pressures on the seismic behavior of the continuous beam bridge. The results indicate that the peak ground acceleration in the saturated regions has significantly spatial variances, and multipoint input can cause the peak acceleration of seismic motion in a saturated site to be amplified by 2~3 times. Consequently, situation leads to a 95%~110% increase in the displacement of the pier top compared to the results under uniform ground motions. Furthermore, dynamic hydrodynamic pressure is shown to precipitate pier displacement and internal forces by 34% and 29%, respectively. Disregarding the influence of saturated site conditions when evaluating ground motions and bridge responses for long-span bridges in saturated complex sites tends to engender hazardous outcomes.
-
引言
据河北省地震台网测定，2018年2月12日18时31分在河北省廊坊市永清县(39.37°N，116.67°E)发生4.3级地震，震源深度20km。此次地震未造成人员伤亡，在震中区域造成个别房屋墙皮脱落，河北省廊坊市、固安市震感较强，北京市、天津市普遍有感。地震发生后，中国地震局工程力学研究所对永清4.3级地震强震动记录进行了远程和现场记录回收。强震动记录是地震能量在地表或建筑结构上的真实反映，是工程地震、地震学研究的基础数据资料(朱升初等，2018)，本文对此次地震构造背景及强震动记录进行了分析，对永清4.3级地震烈度采用建立虚拟台站方法进行了校正，为河北地区烈度调查、地震区划等提供基础资料。
1. 发震构造背景
永清4.3级地震区域构造主要以NE、NW向断裂为主，地震发生在廊固凹陷内(图 1)，廊固凹陷是冀中坳陷北部的一个古近系湖盆，整个凹陷大体呈NE走向，南北长约90km，东西宽20—40km，面积约2600km2。该凹陷是受区域伸展并伴有走滑作用形成的典型箕状凹陷，整个凹陷表现为北断南超、西断东超的构造格局。由于受NE向断裂构造控制，表现出整体不均匀下降的箕状断陷特征。该凹陷内曾于1621年发生廊坊永清$ 5\frac{1}{2} $级地震(国家地震局震害防御司，1995)。
图 1 永清4.3级地震构造背景Figure 1. The seismotectonic background map of the Yongqing MS4.3 earthquake(1)施庄断裂(2)南口山前断裂(3)南口-孙河断裂(4)小汤山-东北旺断裂(5)黄庄-高丽营断裂 (6)顺义-前门-良乡断裂(7)南苑-通县断裂(8)新夏垫断裂(9)夏垫断裂(10)大兴凸起东缘断裂 (11)蓟县山前断裂(12)宝坻断裂(13)桐柏断裂(14)河西务断裂(15)天津断裂(16)大城东断裂 (17)沧东断裂(18)海河断裂(19)牛驼镇凸起东缘断裂(20)徐水断裂(21)东垒子-涞水断裂 (22)永定河断裂(23)徐水南断裂本次地震发生在廊固凹陷东侧控制断裂河西务断裂带附近，该断裂北起河西务西北，向西南经大王务，延伸至永清东码头镇，全长近50km，走向NE，倾向SE，为上陡下缓的铲形正断层，断层产状与震源机制参数基本一致(王晓山等，2018)，所以永清4.3级地震发震构造可能为河西务断裂。周月玲等(2018)通过浅层地震勘探和钻孔联合剖面探测，结合年代样品测试，对河西务断裂活动性进行了综合研究，揭示断裂上断点埋深约150m或以浅，第四系底界面垂直断错为20—45m，断裂最新活动时代为晚更新世早期。
2. 观测台站及记录获取
2.1 强震动台站分布
本次地震强震动数据主要由中国地震局工程力学研究所提供，共获取3个分量自由场记录74组，常规处理后得到222条3分向加速度记录，获取记录的强震动台站均为土层台，记录台站分布如图 2所示，由图 2可知，触发台站主要分布在震中北部，南部无台站触发，这可能是因为首都圈台站较密集，而河北省境内强震动台站稀疏，或与仪器型号不同且很多台站老旧未改造有关。触发台站的震中距为29.4—185.5km，其中震中距＜50km的记录有5组，震中距为50—100km的有29组，震中距为100—150km的有37组，震中距为150—200km的有3组，距震中最近的台站为长子营台(11ZZY)，震中距为29.4km。虽然此次地震获取的强震动记录较全面，但由于震级较小，未造成大的灾害损失，且震中附近台站相对稀疏，能供工程使用的强震动记录较少。
图 2 永清4.3级地震触发台站分布图Figure 2. The distribution of seismic trigger station of Yongqing MS4.3 earthquake2.2 强震动记录获取
在强震动记录使用分析前，首先对波形数据进行消除基线漂移，即原始强震动记录通过减去事前记录平均值的方法进行零线调整(吕俊强等，2015)，通过分析计算得到震中距R＜50km的5个台站强震动记录基本信息及相关参数，见表 1。其中距震中最近的长子营台记录到的最大峰值加速度为7.9cm/s2，觅子店台记录到的东西向最大峰值加速度为36.1cm/s2、南北向最大峰值加速度为40.4cm/s2，红星中学台记录到的垂直方向最大峰值加速度为10.8cm/s2，5个台站加速度时程曲线如图 3。由表 1可知，离震中较近的长子营台(震中距29.4km)较离震中较远的觅子店台(震中距38.4km)和于家务台(震中距36.8km)各分量加速度小，影响地震动强度的因素除震中距外，还包括仪器型号、场地特征条件和地质构造、场地土层结构、覆盖层厚度、地形及断裂分布等(王海云等，2010；温瑞智等，2013；郭秋娜等，2018)。本文根据温瑞智等(2015)提出的适用于我国强震动台站场地划分经验曲线，采用H/V谱比法对本次地震触发台站进行分类，分类结果见图 4。由图 4可知，记录到的本次地震台站大部分为Ⅰ类场地(51.35%)，其次为Ⅱ类场地(35.14%)和Ⅲ类场地(13.51%)，其中，觅子店台和红星中学台为Ⅰ类场地，于家务台和张家湾中学台为Ⅱ类场地，长子营台为Ⅲ类场地。长子营台与觅子店台和于家务台记录差异较大，这可能与场地类型有关。
表 1 强震动记录(震中距＜50km)及相关参数Table 1. Strong motion records(epicenter distance < 50km) and related parameters台站名称 台站代码 仪器型号 震中距/km PGA/cm·s-2 EW SN UD 觅子店 11MZD REFTEK/SLJ-100 38.4 36.1 40.4 8.5 红星中学 11YHZ REFTEK/SLJ-100 43.4 -5.4 -4.6 10.8 于家务 11YJW ETNA/ES-T 36.8 -21.7 -18.5 -8.1 张家湾中学 11ZJW REFTEK/SLJ-100 47.9 -5.8 4.2 -4.4 长子营 11ZZY REFTEK/SLJ-100 29.4 5.0 -7.9 -3.5 3. 谱分析
加速度反应谱是目前国内外工程抗震设计的重要依据，反应谱形状、幅值及卓越周期的差异体现了场地条件、震级对地面运动反应谱的影响(温瑞智等，2013；徐钦等，2017)。为避免地震时建筑结构发生共振，确保结构安全，有必要对加速度反应谱成分进行分析。本文选取震中距＜50km的5个典型强震动台站进行频谱分析，阻尼比取0.05，得到本次地震EW、SN、UD向加速度反应谱(图 5)。
由图 5可知，觅子店台记录到了本次地震EW向和SN向加速度反应谱最大值，红星中学台记录到了本次地震SN向加速度反应谱最大值。觅子店台与长子营台在EW和SN向差别最大，觅子店台与张家湾中学台在SN向差别最大。5个典型台站记录到的加速度反应谱高频成分丰富，谱卓越周期主要为0.1—0.3s，0.3s后迅速下降。根据现场调查，震区主要房屋类型为土木结构、木房顶砌体结构、预制板砌体结构、钢结构简易房屋4种，本次地震除少数土木结构房屋原有裂缝扩大、墙皮开裂外，其他类型房屋未发现新开裂、屋顶掉瓦、墙皮掉灰现象，中小城市主要建筑物结构自振周期为0.3—1s(王文才等，2018)，因此本次地震对该震区自振周期范围内的建筑物破坏影响较小。通过对比不同震中距台站各方向反应谱特征周期发现，震中距对特征周期的影响不明显。
傅氏谱能反映地震波在频域上的幅值特性，地震动记录中包含多种振动频率，不同地区的主要振动频率不同，结构自身具有不同频率、幅频响应，因而在不同频率振动作用下，结构振动不同。为避免地震作用下结构发生共振现象，确保结构安全，有必要对地震频率成分进行分析(冉志杰等，2012)。
本文计算了永清4.3级地震震中距50km以内的5个典型台站强震动记录傅氏谱(图 6)，由图 6可知典型台站傅氏谱谱型主要为单峰值型，强震动记录EW分量傅氏谱最大峰值为3.5—8.5Hz，SN分量傅氏谱最大峰值为2.9—6.2Hz，UD分量傅氏谱最大峰值为2.1—9.8Hz。总体来看，UD分量频率范围分布较宽，EW分量傅氏谱最大峰值大于SN分量和UD分量傅氏谱峰值。一般而言，随着震中距的增大，地震动记录长周期(低频)越显著，但本次地震5个典型台站特征不明显，可能与仪器型号及场地特征条件有关。另外，EW分量傅氏谱最大峰值随着震中距的增大呈先减小后增大的趋势，而SN分量和UD分量无该特征。
本文以东西向为正方向，统计峰值加速度反应谱与永清4.3级地震触发台站同震中所处方位角和震中距的关系(图 7)。由图 7可知，峰值加速度反应谱在同一震中距下整体随着方位角的增大而增大，达70°左右后下降，这可能与附近河西务断裂构造走向有关；峰值加速度反应谱值随着震中距的增加而减小。
图 7 峰值加速度反应谱与震中距、方位角关系图Figure 7. The relation among the peak ground acceleration and epicentral distance、azimuth angles4. 地震烈度分析
地震烈度是衡量地震破坏和地震动强弱程度的重要指标。Wald等(1999)认为当烈度较小时，用峰值加速度计算烈度较合适，因为人对加速度较敏感。随着强震动观测技术的发展，由于以各观测点峰值加速度作为依据可真实反映地震发生时各点地面运动强度(田秀丰等，2013)，基于地震动记录快速确定地震仪器烈度，越来越受到研究者的重视(崔建文等，2016)。目前国内外基于强震动记录资料的地震烈度计算方法较多，不同方法利用不同的地震动参数，包括时间、频度和强度等，建立地震烈度与地震动参数之间的关系。本文采用美国ShakeMap烈度计算方法(李俊等，2010)，计算永清4.3级地震事件记录台站烈度，并采用克里格插值方法进行网格化插值，绘制等值线图(图 8)。Wald等(1999)的研究认为，当烈度低于Ⅴ度时峰值加速度PGA与烈度相关性明显高于峰值速度PGV，而当烈度为Ⅴ—Ⅶ度时，PGA和PGV与烈度相关程度相当，考虑本次地震烈度较低，故采用峰值加速度计算烈度。由图 9可知，本次地震仪器烈度等震线长轴沿NW走向，且极震区位于震中东北部，而不是震中附近，但由现场宏观烈度调查(图 8)与震源机制解(王晓山等，2018)可知该地震等震线长轴沿NE向展布，这种差异现象的出现可能与台站分布不均匀有关，震中南部无强震动记录，所以图 8中震中北部附近出现数据空白现象，不利于获取实际烈度长轴方向真实性。
图 8 永清4.3级地震校正前仪器烈度分布图Figure 8. The distribution of the un-corrected instrumental seismic intensity of Yongqing MS4.3 earthquake可靠快速生成地震烈度可为人员伤亡估计、应急救援决策和工程抢险修复决策提供重要依据，为解决由于缺少地震动观测值出现空白区域现象及烈度长轴方向和实际方向不一致的问题，本文首先利用永清4.3级地震中获取的所有台站加速度峰值，按地震动峰值衰减模型(王国新，2001)确定该地区地震动衰减关系，然后在整个地震区内建立空间随机假设台站进行补点插值，得到新的地震烈度分布图。地震动峰值衰减模型(王国新，2001)如下：
$$ \lg Y = {{\rm{C}}_1} + {{\rm{C}}_2}M + {{\rm{C}}_3}M\lg \left({R + {R_0}} \right) + {{\rm{C}}_4}\lg \left({R + {R_0}} \right) \pm \varepsilon $$ (1) 式中，Y为PGA(峰值加速度)，M为震级，R为震中距，C1、C2、C3、C4为待求常数，R0为近场饱和因子，ε为统计分析误差项。采用通用全局优化法拟合得到地震动峰值加速度随震级和震中距的变化规律，考虑本次地震震级较小，故选取震中距＜100km的Ⅰ类场地强震动数据进行拟合，其中水平向按EW、NS向加速度峰值矢量拟合回归分析，拟合曲线如图 9，计算得到的衰减关系式为：
$$ \lg Y = 1.6882 + 0.3926M - 0.1546M\lg \left({R + 5} \right) + 0.7000\lg \left({R + 5} \right) \pm 0.2287 $$ (2) 
图 9 永清4.3级地震现场调查烈度分布图Figure 9. The distribution of site investigation intensity of Yongqing MS4.3 earthquake峰值加速度随震中距的变化如图 10所示，由图 10可知峰值加速度随着震中距的增加而减小。同一震中距不同台站测得的峰值加速度离散性较大。通过拟合各台站实际PGA值，与俞言祥等(2013)地震动衰减关系进行对比，可知20km范围内2条曲线吻合较好，而20km外拟合曲线呈衰减相对缓慢状态，这可能与本次地震震级较小、记录台站均为土层台、远场有一定放大作用有关(王海云，2011)。由观测值得到的预测方程虽存在一定离散性，但整体趋势与实际结果保持一致。
图 10 峰值加速度随震中距的变化Figure 10. `The change of peak ground acceleration with epicentral distance在震区空白区域随机产生N个空间坐标的假设台站，按式(2)计算水平峰值加速度，按美国ShakeMap烈度计算方法转换成仪器烈度，绘制增加假设台站的烈度分布图(图 11)，由图 11可知，增加假设台站的仪器烈度分布图与现场烈度调查分布图极震区方向基本一致，进一步验证了河西务断裂为本次地震发震构造的可能，解决了台站分布不均匀导致地震烈度影响场计算出现缺值现象的问题，为以后缺少台站记录震区提供准确快速制作烈度分布图的思路，为震害调查和地震应急救援提供重要依据。
图 11 永清4.3级地震校正后烈度分布图Figure 11. The distribution of the corrected instrumental seismic intensity of Yongqing MS4.3 earthquake5. 结论
本文介绍了河北永清4.3级地震发震构造背景，并对国家强震动台网中心在本次地震中获取的加速度记录特征进行了分析，为河北省地震动特征和震害调查提供了一定参考，结论如下：
(1) 河北永清4.3级地震在廊固凹陷东侧控制断裂河西务断裂带附近，地震烈度极震区走向与该断裂展布方向一致，结合震源机制解初步判定本次地震发震构造为河西务断裂。
(2) 74个强震动台站获取到强震动加速度记录，震中距位于29.4—185.5km，其中觅子店台记录到的本次地震东西向最大峰值加速度为36.1cm/s2，垂直方向最大加速度为40.4cm/s2，于家务台记录到的南北向最大峰值加速度为18.5cm/s2。采用H/V谱比法对触发台站进行分类，可知Ⅰ类场地居多。
(3) 通过对5个典型台站反应谱分析发现永清4.3级地震加速度反应谱高频成分丰富，谱卓越周期主要为0.1—0.3s。傅氏谱谱型主要为单峰值型，UD分量频率范围分布较宽，EW分量傅氏谱最大峰值大于SN分量和UD分量。
(4) 通过分析峰值加速度反应谱与震中距、方位角的关系，发现峰值加速度反应谱整体随着方位角的增大而增大，达70°左右后下降，这可能与附近河西务断裂构造走向有关，有待进一步研究。其次，峰值加速度反应谱值整体随着震中距的增加而减小。
(5) 由本次地震仪器烈度分布图可知，台站分布不均匀性直接影响烈度分布的真实性，为避免此类问题，拟合出该地区地震动衰减关系，随机产生空间假设台站进行弥补，为以后缺少台站记录震区提供准确快速制作烈度分布图的思路，为震害调查和地震应急救援提供重要依据。由于收集到的强震动数据有限，文中仅采用本次地震获取的峰值加速度拟合当地地震动衰减关系，拟合结果可能存在误差，今后将收集更多的强震动数据进行拟合，以校正该地区地震动衰减关系。
-
图 3 地表点1、2、3地震动加速度自功率谱密度
Figure 3. Auto-power spectral density of ground motion acceleration of surface points 1, 2, and 3
图 4 不同地表点间地震动加速度互功率谱密度
Figure 4. Cross-power spectral density of ground motion acceleration among different surface points
图 9 多点、一致地震动作用下连续梁桥墩顶位移
Figure 9. Top displacement of the piers of the continuous beam bridge under multi-support and uniform ground motions
图 10 多点、一致地震动作用下连续梁桥墩底剪力
Figure 10. Shear force of the pier bottom of the continuous beam bridge under multi-support and uniform ground motions
图 11 多点、一致地震动作用下连续梁桥墩底弯矩
Figure 11. Moments of the pier bottom of the continuous beam bridge under multi-support and uniform ground motions
图 12 多点地震动作用下动水压力对连续梁桥墩顶位移的影响
Figure 12. Effect of hydrodynamic pressures on the displacement at the pier top of the continuous beam bridge
图 13 多点地震动作用下动水压力对连续梁桥墩底剪力的影响
Figure 13. Effect of hydrodynamic pressures on the shear force at the pier bottom of the continuous beam bridge
-
何颖, 于琴, 刘中宪, 2019. 考虑散射效应沉积河谷空间相关多点地震动模拟. 岩土力学, 40（7）: 2739—2747, 2788He Y. , Yu Q. , Liu Z. X. , 2019. Simulation of multi-point spatially correlated earthquake ground motions of sedimentary valleys considering scattering effect. Rock and Soil Mechanics, 40(7): 2739—2747, 2788. (in Chinese) 黄磊, 刘中宪, 张雪等, 2020. 含流体层的河谷场地对地震波散射的间接边界元法模拟. 地震学报, 42（6）: 657—668Huang L. , Liu Z. X. , Zhang X. , et al. , 2020. IBEM simulation of seismic wave scattering by valley topography with fluid layer. Acta Seismologica Sinica, 42(6): 657—668. (in Chinese) 李丽, 赵国峰, 李吉等, 2020. 一致激励与多点激励对悬索桥地震响应影响分析. 震灾防御技术, 15（2）: 252—259Li L. , Zhao G. F. , Li J. , et al. , 2020. Analysis of the influence of uniform excitation and multi-point excitation on the seismic response of suspension bridge. Technology for Earthquake Disaster Prevention, 15(2): 252—259. (in Chinese) 李平, 薄景山, 肖瑞杰等, 2018. 地震动河谷场地效应研究. 震灾防御技术, 13（2）: 331—341Li P. , Bo J. S. , Xiao R. J. , et al. , 2018. The study of effect by the valley site on ground motion. Technology for Earthquake Disaster Prevention, 13(2): 331—341. (in Chinese) 李伟华, 赵成刚, 2008. 平面SV波在饱和土半空间中圆柱形孔洞周边的散射. 地震工程与工程振动, 28（6）: 1—7Li W. H. , Zhao C. G. , 2008. An analytical solution for the scattering of plane SV-waves around cylindrical cavity in a fluid-saturated porous media half space. Journal of Earthquake Engineering and Engineering Vibration, 28(6): 1—7. (in Chinese) 李忠献, 李笑穹, 李宁, 2014. 空间相关多点多维地震动的模拟. 地震工程与工程振动, 34（4）: 64—72Li Z. X. , Li X. Q. , Li N. , 2014. Simulation of multi-point and multi-dimension spatially correlated earthquake ground motions. Earthquake Engineering and Engineering Dynamics, 34(4): 64—72. (in Chinese) 梁建文, 吴孟桃, 巴振宁, 2021. 流体饱和半空间二维地形三分量弹性波散射间接边界元模拟. 地球物理学报, 64（8）: 2766—2779Liang J. W. , Wu M. T. , Ba Z. N. , 2021. IBEM simulation of three-component scattering of elastic waves in a fluid-saturated half-space with 2 D topography. Chinese Journal of Geophysics, 64(8): 2766—2779. (in Chinese) 闫磊, 李青宁, 岳克锋等, 2019. 多维多点激励下考虑支座摩擦滑移及结构碰撞的非规则桥梁抗震性能研究. 世界地震工程, 35（2）: 68—77Yan L. , Li Q. N. , Yue K. F. , et al. , 2019. Seismic behavior of irregular bridges under multi-dimensional and multi-point excitation considering friction slip of bearings and structural pounding. World Earthquake Engineering, 35(2): 68—77. (in Chinese) Bi K. M. , Hao H. , 2012. Modelling and simulation of spatially varying earthquake ground motions at sites with varying conditions. Probabilistic Engineering Mechanics, 29: 92—104. doi: 10.1016/j.probengmech.2011.09.002 Huang H. C. , Chiu H. C. , 1995. The effect of canyon topography on strong ground motion at Feitsui damsite: quantitative results. Earthquake Engineering & Structural Dynamics, 24(7): 977—990. Li X. Q. , Li Z. X. , Crewe A. J. , 2018. Nonlinear seismic analysis of a high-pier, long-span, continuous RC frame bridge under spatially variable ground motions. Soil Dynamics and Earthquake Engineering, 114: 298—312. doi: 10.1016/j.soildyn.2018.07.032 Liu Z. X. , Li W. X. , Jin L. G. , et al. , 2023 a. Efficient simulation of stochastic seismic response of long-span bridges in river valleys using hybrid BEM-FEM. Soil Dynamics and Earthquake Engineering, 165: 107690. doi: 10.1016/j.soildyn.2022.107690 Sánchez-Sesma F. J. , Campillo M. , 1991. Diffraction of P, SV, and Rayleigh waves by topographic features: a boundary integral formulation. Bulletin of the Seismological Society of America, 81(6): 2234—2253. Wang P. G. , Zhao M. , Du X. L. , 2019. Simplified formula for earthquake-induced hydrodynamic pressure on round-ended and rectangular cylinders surrounded by water. Journal of Engineering Mechanics, 145(2): 04018137. doi: 10.1061/(ASCE)EM.1943-7889.0001567 Wu Y. X. , Gao Y. F. , Zhang N. , et al. , 2016. Simulation of spatially varying ground motions in V-shaped symmetric canyons. Journal of Earthquake Engineering, 20(6): 992—1010. doi: 10.1080/13632469.2015.1010049 Zanardo G. , Hao H. , Modena C. , 2002. Seismic response of multi-span simply supported bridges to a spatially varying earthquake ground motion. Earthquake Engineering & Structural Dynamics, 31(6): 1325—1345. 期刊类型引用(1)
1. 刘敏. 公路连续梁桥体结构抗震性能数值模拟. 黑龙江交通科技. 2024(10): 81-84 . 
百度学术其他类型引用(0)
-
下载:
下载:
点击查看大图
计量
- 文章访问数: 160
- HTML全文浏览量: 79
- PDF下载量: 22
- 被引次数: 1
