基于FEM-IBIEM的河谷坝址地震动输入方法及地震动场模拟研究

PDF(1779 KB)

振动与冲击 ›› 2025, Vol. 44 ›› Issue (6) : 244-253.

地震科学与结构抗震

基于FEM-IBIEM的河谷坝址地震动输入方法及地震动场模拟研究

王飞1,2，刘云贺*3，李虎1,4，薛炳1,4，牟猷1,4，宋志强3

作者信息 +

Seismic input method and ground motion field simulation of canyon dam sites based on FEM-IBIEM

WANG Fei1,2, LIU Yunhe*3, LI Hu1,4, XUE Bing1,4, MOU You1,4, SONG Zhiqiang3

Author information +

文章历史 +

摘要

地震动输入是大坝抗震安全分析的基础，然而目前地震动输入尚未合理考虑河谷散射效应对大坝动力响应的影响。本文将河谷场地总运动场分解为平坦半空间自由场和河谷散射场分别求解，采用间接边界积分方程法求解河谷场地散射场，将其与平坦半空间自由场叠加获得河谷场地总运动场（地震动场），通过封闭解验证了间接边界积分方程法的精度。进而，采用有限元法模拟河谷场地内部，同时将河谷地基截断边界处总运动场转换为等效地震输入荷载，建立了总运动场结合黏弹性人工边界的波动输入方法。应用基于有限元法-间接边界积分方程法建立的地震波动输入方法，数值求解了梯形河谷场地地震动场，分析了自由场波动输入和总运动场波动输入下河谷动力响应的差异。结果表明：自由场波动输入下河谷关键位置位移峰值、波形与精确解存在较大偏差，水平和竖向位移最大误差分别达15.6%、29.3%；总运动场波动输入下位移峰值、波形与精确解拟合良好。本文建立的波动输入方法计算精度高，合理考虑了河谷散射效应，为更准确地预测河谷上大坝动力响应提供了基础。

Abstract

Seismic input is the foundation of seismic safety analysis for dams, but the influence of canyon scattering effects on dam dynamic response has not been reasonably considered in seismic input. This paper decomposes the total motion field of the canyon site into the free field motion of a flat half space and the scattering field motion of the canyon, which are solved separately. The indirect boundary integral equation method is used to solve the scattering field of the canyon site, and it is superimposed with the free field of a flat half space to obtain the total motion field (ground motion field) of the canyon site. The accuracy of the indirect boundary integral equation method is verified through closed form solutions. Furthermore, the finite element method is used to simulate the interior of the canyon site, and the total motion field at the truncated boundary of the canyon foundation is converted into the equivalent seismic input loads, establishing a wave input method that combines the total motion field with viscoelastic artificial boundaries. The seismic wave input method based on finite element method-indirect boundary integral equation method is applied to numerically solve the ground motion field of a trapezoidal canyon site, and the differences in dynamic response of the canyon under free field wave input and total motion field wave input are analyzed. The results show that there is a significant deviation between the peak displacement and waveform of key position of the canyon under the free field wave input and the exact solution, the maximum errors of horizontal and vertical displacements are 15.6% and 29.3%, respectively. The peak displacement and waveform under the total motion field wave input fit well with the exact solution. The wave input method established in this paper has high calculation accuracy and reasonably considers the scattering effect of canyons, providing a basis for more accurate prediction of the dynamic response of dams on canyons.

导出引用

王飞1, 2, 刘云贺3, 李虎1, 4, 薛炳1, 4, 牟猷1, 4, 宋志强3. 基于FEM-IBIEM的河谷坝址地震动输入方法及地震动场模拟研究[J]. 振动与冲击, 2025, 44(6): 244-253

WANG Fei1, 2, LIU Yunhe3, LI Hu1, 4, XUE Bing1, 4, MOU You1, 4, SONG Zhiqiang3. Seismic input method and ground motion field simulation of canyon dam sites based on FEM-IBIEM[J]. Journal of Vibration and Shock, 2025, 44(6): 244-253

参考文献

[1] Chen H Q. Overview of the seismic input at dam sites in China [J]. Earthquake Science, 2022, 35: 410-425.
[2] 王飞，宋志强，刘云贺，等. 地震波斜入射下水平成层覆盖层地基波动输入研究 [J]. 振动与冲击，2024，43(08): 248-256.
WANG Fei, SONG Zhi-qiang, LIU Yun-he, et al. Wave input for foundation with horizontal layered overburden under seismic wave oblique incidence[J]. Journal of Vibration and Shock, 2024, 43(08): 248-256.
[3] Qu, Y Q, Zou D G, Liu J M, et al. Two-dimensional DEM-FEM coupling analysis of seismic failure and anti-seismic measures for concrete faced rockfill dam [J]. Computers and Geotechnics, 2022, 151: 104950.
[4] 刘晶波，谭辉，宝鑫，等. 土–结构动力相互作用分析中基于人工边界子结构的地震波动输入方法 [J]. 力学学报，2018，50(1)：32-43.
LIU Jing-bo, TAN Hui, BAO Xin, et al. The seismic wave input method for soil-structure dynamic interaction analysis based on the substructure of artificial boundaries [J]. Chinese Journal of Theoretical and Applied Mechanics, 2018, 50(1): 32-43.
[5] 张佳文，李明超，韩帅，等. 非一致地震作用下高坝－复杂地基体系动力响应分析 [J]. 水利学报，2023，54(9):1099-1111.
ZHANG Jia-wen, LI Ming-chao, HAN Shuai, et al. Dynamic response analysis of high dam-complex foundation system under non-uniform seismic waves [J]. Journal of Hydraulic Engineering, 2023,54(9):1099-1111.
[6] 宝鑫，刘晶波，王东洋，等. 局部成层的海域岛礁场地地震反应分析 [J]. 振动与冲击，2020，39(21): 55-64.
BAO Xin, LIU Jing-bo, WANG Dong-yang, et al. Seismic response analysis of local layered sea reef site [J]. Journal of Vibration and Shock, 2020，39(21): 55-64.
[7] 李明超，张佳文，张梦溪，等．地震波斜入射下混凝土重力坝的塑性损伤响应分析 [J]．水利学报，2019，50(11): 1236-1339．
LI Ming-chao, ZHANG Jia-wen, ZHANG Meng-xi, et al. Plastic damage response analysis of concrete gravity dam due to obliquely incident seismic waves [J]. Journal of Hydraulic Engineering, 2019, 50(11): 1236-1339．
[8] Song, Z Q, Wang F, Li Y L, et al. Nonlinear seismic responses of the powerhouse of a hydropower station under near-fault plane P-wave oblique incidence [J]. Engineering Structures, 2019, 199: 109613.
[9] 何建涛，马怀发，张伯艳，等. 黏弹性人工边界地震动输入方法及实现 [J]. 水利学报，2010，41(8): 960-969.
HE Jian-tao, MA Huai-fa, ZHANG Bo-yan, et al. Method and realization of seismic motion input of viscous-spring boundary [J].Journal of Hydraulic Engineering, 2010，41(8):960-969.
[10] Du X L, Zhao M. A local time-domain transmitting boundary for simulating cylindrical elastic wave propagation in infinite media [J]. Soil Dynamic Earthquake Engineering, 2010, 30(10): 937-946.
[11] 赵密. 近场波动有限元模拟的应力型时域人工边界条件及其应用 [D]. 北京: 北京工业大学，2009.
[12] 王飞，宋志强，刘云贺，等. SV波空间斜入射下沥青混凝土心墙响应特性及抗拉破坏评价研究[J]. 岩土工程学报，2023，45(08): 1733-1742.
WANG Fei, SONG Zhi-qiang, LIU Yun-he, et al. Response characteristics and tensile failure evaluation of asphalt concrete core wall under spatial oblique incidence of SV waves [J]. Chinese Journal of Geotechnical Engineering, 2023, 45(08): 1733-1742.
[13] 刘云贺，张伯艳，陈厚群. 拱坝地震输入模型中黏弹性边界与黏性边界的比较[J]. 水利学报，2006，37(6):758-763.
LIU Yun-he, ZHANG Bo-yan, CHEN Hou-qun. Comparison of spring-viscous boundary with viscous boundary for arch dam seismic input model [J]. Journal of Hydraulic Engineering, 2006, 37(6): 758-763.
[14] 周晨光. 高土石坝地震波动输入机制研究[D]. 大连：大连理工大学，2009.
[15] 周晨光，孔宪京，邹德高，等. 地震波动输入方法对高土石坝地震反应影响研究[J]. 大连理工大学学报，2016，56(04): 382-389.
ZHUO Chen-guang, KONG Xian-jing, ZOU De-gao, et al. Study of influence of earthquake wave motion input method on seismic response of high earth-rock dam [J]. Journal of Dalian University of Technology, 2016, 56(04): 382-389.
[16] 刘晶波，吕彦东. 结构-地基动力相互作用问题分析的一种直接方法[J]. 土木工程学报，1998，31(3): 55-64.
LIU Jing-bo, LU Yan-dong. A direct method for analysis of dynamic soil-structure interaction[J]. China Civil Engineering Journal, 1998, 31(3): 55-64.
[17] 赵磊超，宋志强，王飞，等. 基于IBIEM的河谷场地非一致地震动特性及影响因素分析 [J]. 振动与冲击，2022, 41(23):109-120.
ZHAO Lei-chao, SONG Zhi-qiang, WANG Fei, et al. Non-uniform ground motion characteristics and influencing factors of valley site under oblique incidence of SH wave[J]. Journal of Vibration and Shock, 2022, 41(23):109-120.
[18] Xu S T, Xu Q, Chen J Y, et al. Improved endurance time analysis for seismic responses of concrete dam under near-fault pulse-like ground motions [J]. Engineering Structures, 2022, 270:114912.
[19] Zhang, J W, Li, M C, Han, S. Seismic analysis of gravity dam-layered foundation system subjected to earthquakes with arbitrary incident angles. International Journal of Geomechanics, 2022, 22(2): 04021279.
[20] Huang J Q, Zhao M, Du X L. Non-linear seismic responses of tunnels within normal fault ground under obliquely incident P waves. Tunnelling and Underground Space Technology, 2017, 61: 26-39.
[21] Zhao W S, Chen W Z, Yang D S, et al. Earthquake input mechanism for time-domain analysis of tunnels in layered ground subjected to obliquely incident P- and SV-waves[J]. Engineering Structures, 2019, 181:374-386.
[22] 赵瑞斌，刘萌，刘中宪，等. 沉积谷中衬砌隧道对SH波的散射 [J]. 岩土力学，2016， 37(9): 2687-2696.
ZHAO Rui-bin, LIU Meng, LIU Zhong-xian, et al. Scattering of SH-wave along lined tunnel in alluvial valley [J]. Rock and Soil Mechanics, 2016, 37(9):2687-2696.
[23] Huang L, Liu Z X, Wu C Q, et al. A three-dimensional indirect boundary integral equation method for the scattering of seismic waves in a poroelastic layered half-space [J]. Engineering Analysis with Boundary Elements, 2022, 135:167-181.
[24] Trifunac M D. Scattering of plane SH wave by a semi-cylindrical canyon [J]. Earthquake Engineering & Structural Dynamics, 1973, l: 267-281.
[25] Yuan X M, Liao Z P. Scattering of plane SH waves by a cylindrical alluvial valley of circular‐arc cross‐section[J]. Earthquake Engineering & Structural Dynamics, 1995, 24(10): 1303-1313.
[26] 梁建文，严林隽，李军伟，等. 圆弧形沉积河谷场地在平面P波入射下的响应[J]. 岩土力学，2001，22(2): 138-143.
LIANG Jian-wen, YAN Lin-jun, LI Jun-wei, et al. Response of circular-arc alluvial valleys under incident plane P waves [J]. Rock and Soil Mechanics, 2001, 22(2): 138-143.
[27] Zhang N, Wang L, Zhang Y, et al. Effect of a V-shaped canyon on the seismic response of a bridge under oblique incident SH waves[J]. Earthquake Engineering & Structural Dynamics, 2024, 53(1): 496-514.
[28] 李艳朋，李志远，胡志强，等. 基于比例边界有限元法的成层半空间中斜坡场地地震响应计算模型[J]. 岩土工程学报，2023，45(5): 1008-1016.
LI Yan-peng, LI Zhi-yuan, HU Zhi-qiang, et al. Model for seismic dynamic response of slope terrain site in layered half-space based on scaled boundary finite element method [J]. Chinese Journal of Geotechnical Engineering, 2023, 45(5): 1008-1016.
[29] 李艳朋，李志远，胡志强，等. 基于改进比例边界有限元法的成层场地中峡谷−地下孔洞体系的散射分析[J]. 岩土力学，2022，43(Z2): 553-562.
LI Yan-peng, LI Zhi-yuan, HU Zhi-qiang, et al. A modified scaled boundary finite element method for scattering analysis of canyon-underground cavity system in horizontally layered site [J]. Rock and Soil Mechanics, 2022, 43(Z2): 553-562.
[30] Liu Z X, Sun J, Meng S B, et al. 3D coupled FEM-IBIEM for non-linear seismic response of a sedimentary basin using equivalent-linearized soil model. Computers & Geosciences, 2022, 167: 105199.
[31] Liu Z X, Huang L, Liang JW, et al. A three-dimensional indirect boundary integral equation method for modeling elastic wave scattering in a layered half-space. International Journal of Solids Structures, 2019, 169:81-94.
[32] Zhang C, Liu Q J, Deng P. Surface motion of a half‐space with a semicylindrical canyon under P, SV, and Rayleigh waves [J]. Bulletin of the Seismological Society of America, 2017, 107(2): 809-820.
[33] Wong H L. Diffraction of plane P, SV and Rayleigh waves by surface topography [R]. Los Angeles: University