Analysis of Slope Terrain Effect on the Properties of Ground Motion
-
摘要: 采用有限元有限差分方法,结合人工透射边界理论,研究局部坡地地形对地震动特性的影响,分析坡高、坡角对地形放大效应的影响。研究结果表明:坡地地形斜坡段各点反应谱谱比最大值沿坡高逐渐增大,坡脚点对地震动反应谱谱比呈缩小效应;坡底段各点反应谱谱比接近1;坡起平台段空间点受地形效应的影响较大;坡高和坡角对地形效应的影响较明显,当坡角不变时,不同坡高反应谱谱比曲线形状基本相同,峰值点对应的特征周期随着坡高的增加呈增长趋势,峰值点对应的反应谱谱比随着坡高的增加呈增大趋势,斜坡顶点阿里亚斯烈度比值在地形坡高超过一定数值后随着坡高的增加呈增大趋势;当坡高不变时,大部分频段内反应谱谱比最大值随着坡角的增加逐渐增大,斜坡顶点阿里亚斯烈度比值在坡角达到一定数值后随着坡角的增加呈增大趋势。Abstract: Based on the finite element finite difference method and the theory of artificial transmission boundary, the influence of local slope terrain on ground motion characteristics (including response spectrum, peak acceleration, peak velocity, peak displacement and arias intensity) is studied, and the influence of slope height and slope angle on ground shape amplification effect is analyzed. The results show that: the maximum response spectrum ratio of each point on the slope section of slope terrain gradually increases along the slope height, and the maximum value is not more than 1.3; the slope toe point has a narrowing effect on the response spectrum ratio of ground motion; the spectrum ratio of each point on the bottom section of slope is close to 1; the spatial point on the platform section is greatly affected by the terrain effect, and gradually decreases with the distance from the top of slope. The height and slope angle of slope terrain have obvious influence on terrain effect. When the slope angle is constant, the shape of the spectral ratio curve of different slope heights is basically the same, the characteristic period corresponding to the peak point increases with the increase of height, and the spectral ratio corresponding to the peak point also increases with the increase of height. When the terrain height exceeds a certain value, the ratio of arias intensity at the top of slope increases with the increase of height; when the slope height is constant, the maximum value of response spectrum ratio in most frequency bands also increases with the increase of slope angle. And the ratio of arias intensity at the top of the slope increases with the increase of the slope after reaching a certain slope angle.
-
Key words:
- Terrain effect /
- Slope terrain /
- Response spectrum ratio /
- Terrainheight /
- Terrain slope /
- Arias intensity
-
引言
重力坝作为重要的基础设施,一旦发生破坏,可能会导致潜在的生命危险和巨大的经济损失,对强震作用下的混凝土重力坝进行抗震安全研究极其重要。
现有研究表明,影响重力坝抗震性能的主要因素有(Chopra,2012;Løkke等,2018):(1)坝体-库水动力相互作用,包含库水可压缩性和库底沉积物引起的库底吸收系数(Hall等,1982;Fenves等,1985);(2)坝体-基岩相互作用,包括岩体的惯性作用(Fenves等,1985;Tan等,1996);(3)无限地基辐射阻尼效应(Tan等,1996;Zhang等,2009);(4)地震动的空间分布和不确定性(Chopra,2010);(5)坝体和基岩的材料非线性(El-Aidi等,1989;Fenves等,1992;Bhattacharjee等,1993;Cervera等,1995;Pan等,2011)。众多研究者通过考虑以上部分或全部因素开展了混凝土重力坝抗震性能评估等相关研究。Alembagheri(2016)采用静力弹塑性分析方法,从线性地震分析结果出发,提出一种系统且合理的损伤程度评估方法;将大坝混凝土的拉伸开裂作为主要破坏模式,以3座现有混凝土重力坝为例,对该方法进行说明,并讨论了可能的非线性响应和破坏机理。Hariri-Ardebili等(2013)首次以专业的方式提出基于应变的混凝土坝结构性能判别准则,并探讨了其在拱坝地震破坏评估中的适用性,利用需求能力比DCR、累积非弹性持续时间和超应力/超应变区等指标对大坝的抗震性能进行研究,结果表明,采用基于应力的准则对拱向作用的评估偏向保守,而采用基于应变准则的抗震安全性评估对梁向的作用评估偏向保守。郑晓东等(2016)基于混凝土塑性损伤模型考虑大坝混凝土材料非线性,针对强震持时对混凝土重力坝损伤破坏累积进行了研究,结果表明,混凝土大坝损伤累积随着地震动持时的增加逐渐增大,局部损伤指标可用于确定大坝抗震薄弱部位,整体损伤指标能够用于大坝地震整体损伤破坏评价。潘坚文等(2010)针对强震输入方式对重力坝的地震响应进行了讨论,分别采用无质量地基模型和弹簧-阻尼边界模型对不同地震荷载、不同基岩和混凝土弹性模量比值下重力坝的地震响应进行对比分析,并提出等效结构阻尼理念和近似方法。殷琳等(2019)开展了水平分层土层系统的等效阻尼比近似计算方法研究,并建议采用基于三角分布的加权函数计算等效阻尼比。杜修力等(2017)针对软、中、硬3种土层场地,选取100条实测地震动记录调幅至0.1 g、0.2 g和0.3 g,并基于一维等效线性方法开展场地随机地震反应研究。
综上可知,上述研究通过考虑影响大坝抗震性能的主要因素,从不同角度对重力坝抗震性能进行研究,获得了丰硕的研究成果,但上述研究均是在重力坝坝基为地质条件较好的基岩前提下开展的。实际工程中,重力坝可能处于软弱覆盖层基础、砂砾石基础等工程地质条件较差的场址,覆盖土层地基情况下的地震动输入及重力坝地震响应与基岩场地存在较大的差异。鉴于此,本文结合国外某强震区深厚覆盖层场地重力坝工程,采用成层状地基地震动输入计算方法、粘弹性边界模型和接触非线性模型,开展超强地震作用下覆盖层场地重力坝的非线性动力分析,结合DCR评价指标,对重力坝抗震安全性进行评估,为重力坝工程设计提供支持。
1. 计算理论
1.1 粘弹性边界与地震动输入
1.1.1 粘弹性边界
众多研究者根据不同假设条件从不同角度开展了考虑坝体和地基相互作用的研究。目前应用较广泛的是粘弹性人工边界(Deeks等,1994;刘晶波,1998),包括在两侧和底部边界每个节点增加弹簧和阻尼器。在有限元计算中,通过在两侧和底部考虑边界弹簧刚度和阻尼系数实现粘弹性人工边界的施加,垂直于边界方向的弹簧系数KN和阻尼系数CN分别为
$ \dfrac{E}{{2{r_{\rm{b}}}}}A $ 、$ \rho {c_{\rm{p}} }A $ ,平行于边界方向的弹簧系数KT和阻尼系数CT分别为$ \dfrac{G}{{2{r_{\rm{b}}}}}A $ 、$ \rho {c_{\rm{s}}}A $ ,其中,E为弹性模量,G为剪切模量,$ \rho $ 为密度,A为人工边界节点影响面积,$ {r_{\rm{b}}} $ 表示从边界底部到顶部的距离,$ {c_{\rm{p}} } $ 和$ {c_{\rm{s}}} $ 分别为有限元模型外侧介质的压缩波波速和剪切波波速。1.1.2 地震动输入
对于成层状地基,其地震动输入采用自由场模型输入,详细求解方法见文献(Idriss等,1992)。本文采用的一维波动系统如图1所示,该系统由N个在水平方向上可延伸至无穷远的水平层组成,每一层均匀且各向同性,材料特性包括厚度h、密度ρ、剪切模量G和阻尼系数β。图1所示剪切波竖向传播会产生水平向位移:
$$ u = u\left( {x,t} \right) $$ (1) 频率为
$ \omega $ 的剪切波水平向位移为:$$ u\left( {x,t} \right){\text{ = }}U\left( x \right){{\rm{e}}^{i\omega t}} $$ (2) 位移
$ u\left( {x,t} \right) $ 须满足波动方程:$$ \rho \frac{{{\partial ^2}u}}{{\partial {t^2}}} = G\frac{{{\partial ^2}u}}{{\partial {x^2}}} + \eta \frac{{{\partial ^2}u}}{{\partial x\partial t}} $$ (3) 由式(2)、(3)可得:
$$ \left( {G + i\omega \eta } \right)\frac{{{{\rm{d}}^2}U}}{{{\rm{d}}{x^2}}}{\text{ = }}\rho {\omega ^2}U $$ (4) 其一般解为:
$$ U\left( x \right) = E{{\rm{e}}^{ikx}} + F{{\rm{e}}^{ - ikx}} $$ (5) 其中,
$ {k^2} = \dfrac{{\rho {\omega ^2}}}{{G + i\omega \eta }} = \dfrac{{\rho {\omega ^2}}}{{{G^ * }}} $ ,为复波数;$ {G^ * } $ 为复剪切模量。1.2 接触非线性
对于接触非线性问题, ABAQUS有限元分析软件中通过2种模型对接触压力进行定义。首先是基于拉格朗日乘子法的硬接触模型(ABAQUS,2010),该方法对接触压力
$ p $ 的定义如下:(1) 当
$ h < 0 $ 时,$ p = 0 $ ,表示张开;(2)当
$ h{\text{ = }}0 $ 时,$ p > 0 $ ,表示闭合。其次是基于指数关系的软接触模型,模型中接触面由主面和从面组成,接触面的接触压力-过盈曲线遵循指数关系,如图2所示,表达式如下所示:
$$ \left\{ {\begin{array}{*{20}{l}} {p = 0\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\qquad h \leqslant - c} \\ {p = \dfrac{{{p_0}}}{{e - 1}}\left[ {\left( {\dfrac{h}{c} + 1} \right)\left( {\exp \left( {\dfrac{h}{c} + 1} \right) - 1} \right)} \right]\;\;\;\;\;h > - c} \end{array}} \right. $$ (6) 式中,
$ h $ 表示张开度;$ {p_0} $ 为零过盈时的典型压力值;$ c $ 是初始接触距离。2. 算例与分析
2.1 工程背景
尼泊尔某水电站位于加德满都东北方向约75 km处,距中国边境陆路樟木口岸约5 km,于1997年开工建设,2001年1月建成发电。该水电站为低坝长引水式电站,挡水坝为溢流式混凝土重力坝,电站装机2台,单机容量22.5 MW。
该水电站经历了2次较严重的自然灾害事件:(1)在2015年5月尼泊尔里氏7.8级地震及余震中,重力坝完好无损,水库淤积在右岸1#、2#重力坝段上游区域,基本接近坝顶高程。(2)2016年7月5日因冰川湖溃决引发泥石流,冲击坝址首部区,造成重力坝严重破坏。泥石流严重损毁了1#坝段的坝肩和地基,浆砌石护岸和通往大坝的道路被破坏,中尼高速公路被切断,靠近岸坡的1#坝段地基被掏空。在泥石流的冲刷作用下,右坝肩处形成一条绕过重力坝的新河道。2#坝段下游沿坝顶有裂纹及剥落现象。除1#坝段右侧因河道侵蚀而露出外,坝体上游大部分被冲积层淤积掩盖,需修复或重建以恢复重力坝挡水发电功能。鉴于工程处于高地震烈度区,需对修复结构进行抗震性能评估。
2.2 有限元模型
修复后的重力坝几何模型如图3所示,坝顶高程为1440.5 m,建基面最低高程为1413.0 m,最大坝高为27.5 m,覆盖土层以下为基岩,厚度为50 m。依据几何模型构建重力坝坝体-地基有限元计算模型,如图4所示。基于所构建的有限元模型分别开展考虑覆盖土层地震动输入的线弹性和非线性动力有限元计算分析,对其抗震安全性进行全面论证评估。基岩刚度远大于覆盖层刚度,覆盖层底部的地震动几乎不受上部土层影响,可以考虑为基岩露头处的地震动输入,同时通过以上成层状地基地震动输入计算方法获取两侧边界自由运动,采用粘弹性边界考虑地基辐射阻尼效应。本文所考虑地震动输入模型如图5所示。
2.3 静动力荷载与材料参数
(1)静、动力荷载
静态荷载包括坝体自重、上游水和淤沙荷载、下游水荷载以及扬压力。正常运行上游水位为1434 m,淤沙高程为1425 m,下游水位为1425 m。扬压力从上游坝踵到下游坝趾沿坝基交界面线性分布。地震加速度时程如图6所示,其中运行基准地震OBE和最大设计地震MDE水平向地震动峰值加速度分别为0.65 g、1.2 g,相应的竖向峰值加速度分别为0.54 g、0.99 g。
(2)坝体和地基材料参数
动力计算时,混凝土及覆盖层的阻尼比分别取5%、7%,混凝土、地基材料及各类接触面参数如表1、2所示。
表 1 混凝土及地基材料参数Table 1. Material parameters of concrete and foundation材料 容重γ/kN·m−3 剪切
模量G/MPa泊松比$ \mu $ 弹性模量E/MPa 摩擦角φ/° 容许承载力/kPa 抗压/抗拉强度/MPa 砼C20 24.0 12000 0.167 28000 - 12500 20/2.40 毛石砼
C1224.0 9500 0.167 22000 - 7500 12/1.71 覆盖层 19.0 - 0.200 25 32.5 400 - 岩石 26.5 1800 0.275 4500 - - - 表 2 各类接触面参数Table 2. Parameters of contact surfaces接触面 粘聚力c/kPa 摩擦角φ/° 砼-砼 0 45.0 砼-毛石砼 0 40.0 砼-覆盖层 0 28.8 毛石砼-覆盖层 0 28.8 2.4 结果分析
本工程按照《Gravity dam design》(EM 1110-2-2200)(US Army Corps of Engineers, 1995)、《Time-history dynamic analysis of concrete hydraulic structures》(EM 1110-2-6051)(US Army Corps of Engineers, 2003)、《Stability analysis of concrete structures》(EM 1110-2-2100)(US Army Corps of Engineers, 2005)、《Earthquake design and evaluation of concrete hydraulic structures》(EM 1110-2-6053)(US Army Corps of Engineers, 2007)进行设计,将混凝土应力性能评估利用需求能力比DCR作为关键绩效指标。Alembagheri(2016)、Hariri-Ardebili等(2013)基于DCR方法对混凝土坝的抗震性能进行了评估。通过开展OBE作用坝体-地基线弹性动力时程分析得到坝体应力、坝顶位移及坝基交界面滑动安全系数,如图7~图10所示。由图可知,坝体最大主应力为1.18 MPa,小于允许值2.4 MPa(DCR=1);最小主应力为3.38 MPa,亦小于允许值11.5 MPa;坝顶最大位移为1.32 m;滑动安全系数为0.19。
根据线弹性动力时程分析可得,在OBE作用下,重力坝坝体应力均在允许范围内,但其抗震稳定安全系数仅0.19,难以满足抗震稳定性。为全面评估重力坝的抗震性能,需进一步开展考虑坝体和地基接触非线性的动力时程分析。根据非线性动力时程分析可得OBE和MDE工况下坝体最终滑移值分别为1.92 m、16.06 m,如图11、图12所示。线性和非线性动力分析结果表明,在OBE和MDE作用下,重力坝均不能保持稳定。
为增强重力坝抗震稳定性,在重力坝坝后回填土,同时将地基范围延伸扩展至基岩,回填土材料参数与覆盖土层一致(表1),坝体-地基几何模型及预设接触面如图13所示,相应的坝体-地基非线性有限元计算模型如图14所示。图15、图16分别给出了OBE和MDE作用下坝基接触面滑移时程,由图可知,在OBE作用下,坝基交界面最大滑移量为0.138 m(沿上游方向),残余滑移量为0.059 m(沿上游方向);在MDE作用下,坝基交界面最大滑移量为0.41 m(沿下游方向),残余滑移量为0.025 m(沿下游方向)。综上所述,重力坝坝体下游坝后回填土能够有效增强其抗震稳定性。
3. 结论
本文依据成层状地基地震动输入计算方法得到覆盖层边界的自由场运动,采用粘弹性边界考虑地基辐射阻尼效应,通过线弹性和非线性动力有限元分析,详细论证、评估了超强地震作用下国外某覆盖土层重力坝的抗震安全性,本研究可为强震区覆盖土层重力坝抗震分析提供参考,主要结论如下:
(1)线弹性OBE工况下,坝体最大主应力为1.18 MPa,最小主应力为3.38 MPa,均小于允许值;坝顶最大位移为1.32 m,滑动安全系数为0.19,难以满足抗震稳定性;
(2)通过非线性动力时程分析得到OBE和MDE工况下,坝体最终滑移值分别为1.92 m、16.06 m,重力坝-地基体系无法保持稳定;
(3)为加强重力坝抗震稳定性,在坝体下游坝后回填土,通过建立新的分析模型得到,在OBE作用下,坝基交界面最大滑移量为0.138 m(沿上游方向),残余滑移量为0.059 m(沿上游方向);在MDE作用下,坝基交界面最大滑移量为0.41 m(沿下游方向),残余滑移量为0.025 m(沿下游方向),重力坝的抗震稳定性得到了有效加强。
-
-
[1] 巴振宁, 陈昊维, 梁建文, 2018. 任意多个沉积谷地对平面SV波的散射. 振动与冲击, 37(5): 98—107Ba Z. N., Chen H. W., Liang J. W., 2018. Scattering of arbitrary alluvial valleys to incident plane SV waves. Journal of Vibration and Shock, 37(5): 98—107. (in Chinese) [2] 郝明辉, 张郁山, 2014. 凸起地形对地震动特性的影响. 地震学报, 36(5): 883—894Hao M. H., Zhang Y. S., 2014. Analysis of terrain effect on the properties of ground motion. Acta Seismologica Sinica, 36(5): 883—894. (in Chinese) [3] 郝明辉, 张郁山, 2019. 基于DEM数据的地形效应经验预测模型研究. 土木工程学报, 52(2): 86—96Hao M. H., Zhang Y. S., 2019. Research on empirical prediction model of terrain effects based on DEM. China Civil Engineering Journal, 52(2): 86—96. (in Chinese) [4] 蒋涵, 周红, 高孟潭, 2015. 三维地形中地震动的频域特征——以芦山地区为例. 震灾防御技术, 10(1): 59—67 doi: 10.11899/zzfy20150106Jiang H., Zhou H., Gao M. T., 2015. The characteristics of frequency domain of ground motion in 3-D topography-A case study of Lushan area. Technology for Earthquake Disaster Prevention, 10(1): 59—67. (in Chinese) doi: 10.11899/zzfy20150106 [5] 李家祥, 温瑞智, 刘启方等, 2021. SH波垂直入射下的斜坡地形放大因子分析. 世界地震工程, 37(1): 176—188Li J. X., Wen R. Z., Liu Q. F., et al., 2021. Amplification factor analysis of slope topography under normal incidence of SH wave. World Earthquake Engineering, 37(1): 176—188. (in Chinese) [6] 李小军, 廖振鹏, 杜修力, 1992. 有阻尼体系动力问题的一种显式差分解法. 地震工程与工程振动, 12(4): 74—80Li X. J., Liao Z. P., Du X. L., 1992. An explicit finite difference method for viscoelastic dynamic problem. Earthquake Engineering and Engineering Vibration, 12(4): 74—80. (in Chinese) [7] 李英民, 王丽萍, 赵耀, 2010. 岩质坡地建筑设计水平地震动放大系数的确定. 地震工程与工程振动, 30(4): 159—165Li Y. M., Wang L. P., Zhao Y., 2010. Amplification factors of design horizontal ground motions for structures on rock slopes. Journal of Earthquake Engineering and Engineering Vibration, 30(4): 159—165. (in Chinese) [8] 李英民, 王肖巍, 宋维举, 2019. 动态子结构法在凸起地形地震动响应中的应用. 哈尔滨工业大学学报, 51(6): 156—161, 170 doi: 10.11918/j.issn.0367-6234.201805166Li Y. M., Wang X. W., Song W. J., 2019. Application of dynamic substructure method to protruding topography on characteristics of ground motion. Journal of Harbin Institute of Technology, 51(6): 156—161, 170. (in Chinese) doi: 10.11918/j.issn.0367-6234.201805166 [9] 梁建文, 张郁山, 顾晓鲁等, 2003. 圆弧形层状凹陷地形对平面SH波的散射. 振动工程学报, 16(2): 158—165 doi: 10.3969/j.issn.1004-4523.2003.02.006Liang J. W., Zhang Y. S., Gu X. L., et al., 2003. Scattering of plane SH waves by a circular-Arc layered canyon. Journal of Vibration Engineering, 16(2): 158—165. (in Chinese) doi: 10.3969/j.issn.1004-4523.2003.02.006 [10] 刘甲美, 高孟潭, 陈鲲, 2015. 地形效应影响下地震动参数与斜坡稳定性的相关性研究. 地震学报, 37(5): 865—874 doi: 10.11939/jass.2015.05.014Liu J. M., Gao M. T., Chen K., 2015. On the correlation of ground motion parameters with slope stability incorporating topographic effects. Acta Seismologica Sinica, 37(5): 865—874. (in Chinese) doi: 10.11939/jass.2015.05.014 [11] 刘中宪, 张征, 王少杰等, 2019. 基于快速多极子间接边界元法的城市区域沉积盆地三维地震响应宽频模拟. 岩土力学, 40(10): 4101—4110Liu Z. X., Zhang, Z., Wang, S. J., et al., 2019. 3D seismic response broadband-simulation of the alluvial basin in urban region based on the FMM-IBEM. Rock and Soil Mechanics, 40(10): 4101—4110. (in Chinese) [12] 孙纬宇, 汪精河, 严松宏等, 2019. SV波斜入射下河谷地形地震动分布特征分析. 振动与冲击, 38(20): 237—243, 265Sun W. Y., Wang J. H., Yan S. H., et al., 2019. Characteristic analysis of ground motions of a canyon topography under obliquely incident SV waves. Journal of Vibration and Shock, 38(20): 237—243, 265. (in Chinese) [13] 殷跃平, 2008. 汶川八级地震地质灾害研究. 工程地质学报, 16(4): 433—444 doi: 10.3969/j.issn.1004-9665.2008.04.001Yin Y. P., 2008. Researches on the geo-hazards triggered by Wenchuan Earthquake, Sichuan. Journal of Engineering Geology, 16(4): 433—444. (in Chinese) doi: 10.3969/j.issn.1004-9665.2008.04.001 [14] 袁晓铭, 廖振鹏, 1996. 任意圆弧形凸起地形对平面SH波的散射. 地震工程与工程振动, 16(2): 1—13Yuan X M, Liao Z P, 1996. Scattering of plane SH waves by a cylindrical hill of circular-arc cross-section. Earthquake Engineering and Engineering Vibration, 16(2): 1—13. (in Chinese) [15] 中华人民共和国住房和城乡建设部, 中华人民共和国国家质量监督检验检疫总局. 2010. GB 50011—2010 建筑抗震设计规范. 北京: 中国建筑工业出版社.Ministry of Housing and Urban-Rural Development of the People’s Republic of China (MOHURD), General Administration of Quality Supervision, Inspection and Quarantine of the People's Republic of China. 2010. GB 50011—2010 Code for Seismic Design of Buildings. Beijing: China Architecture & Building Press. (in Chinese) [16] 周国良, 李小军, 侯春林等, 2012. SV波入射下河谷地形地震动分布特征分析. 岩土力学, 33(4): 1161—1166 doi: 10.3969/j.issn.1000-7598.2012.04.029Zhou G. L., Li X. J., Hou C. L., et al., 2012. Characteristic analysis of ground motions of canyon topography under incident SV seismic waves. Rock and Soil Mechanics, 33(4): 1161—1166. (in Chinese) doi: 10.3969/j.issn.1000-7598.2012.04.029 [17] Ashford S. A., Sitar N., Lysmer J., et al., 1997. Topographic effects on the seismic response of steep slopes. Bulletin of the Seismological Society of America, 87(3): 701—709. [18] Boore D. M., 1972. A note on the effect of simple topography on seismic SH waves. Bulletin of the Seismological Society of America, 62(1): 275—284. [19] Borcherdt R. D., 1970. Effects of local geology on ground motion near San Francisco Bay. Bulletin of the Seismological Society of America, 60(1): 29—61. 期刊类型引用(5)
1. 于淳蛟. 基于坝体安全评价模型的复合地基坝体质量安全分析. 水利科技与经济. 2024(02): 161-166 . 百度学术
2. 徐乐意,黄海斌,郁怀光. 主余震对重力坝损伤特性影响. 华北地震科学. 2024(02): 9-13 . 百度学术
3. 汪海洋,郭涛,刘海龙,吴文禄,马鸿泽,虎雪洁. 强震作用下重力坝非线性动力响应分析. 中国水运(下半月). 2024(07): 111-113 . 百度学术
4. 汪海洋,郭涛,刘海龙,吴文禄,马鸿泽,虎雪洁. 强震作用下重力坝非线性动力响应分析. 中国水运. 2024(14): 111-113 . 百度学术
5. 王伟,范建朋. 尼泊尔上博迪克西水电站重力坝修复. 西北水电. 2022(04): 58-62 . 百度学术
其他类型引用(0)
-