基于有限元法的跌水闸闸室结构抗震安全复核研究

郭博文; 鲁立三; 王荆; 宋力; Bo-wen GUO; Li-san LU; Jing WANG; Li SONG

中国农村水利水电 ›› 2021 ›› (7) : 155-160.

农田水利

基于有限元法的跌水闸闸室结构抗震安全复核研究

作者信息 +

Research on the Seismic Safety Recheck of Drop Gate Chamber Structure Based on Finite Element Method

Author information +

稿件信息 +

摘要

抗震安全复核是水闸安全鉴定工作中的重要组成部分，为保证水闸结构安全运行，需要进行抗震安全复核研究。以某跌水闸为例，基于ADINA有限元分析软件，建立了跌水闸闸室结构三维有限元模型，采用振型分解反应谱法对闸室结构进行了地震动响应计算分析，并根据计算结果对跌水闸闸室结构进行了抗震安全评价。计算结果表明：①闸墩与闸底板相交处出现了2.00 MPa拉应力，根据结构力学方法复核，该处能承受的最大拉应力为4.41 MPa，满足安全需求；②跌水闸公路桥跨中底部位置和启闭机房底板部分区域出现了较大拉应力区，其数值远超混凝土动态轴心抗拉强度标准值，为安全起见，宜采取相应的抗震加固措施；③跌水闸闸室结构抗滑稳定安全系数大于《水闸设计规范》（SL265-2016）标准值，满足安全需求；④跌水闸闸室结构抗震安全满足《水工建筑物抗震设计标准》（GB51247-2018）要求，但存在不影响总体安全的缺陷，其抗震等级评定为B级。

Abstract

Seismic safety review is an important part of sluice safety appraisal. In order to ensure the safe operation of the sluice structure， seismic safety review research is required. Taking a drop sluice as an example， based on ADINA finite element analysis software， a three-dimensional finite element model of the drop gate chamber structure is established， and the seismic response of the chamber structure is calculated and analyzed by using the mode decomposition response spectrum method， and the seismic safety evaluation of the sluice is carried out according to the calculation results. The results show that： ①there is a tensile stress of 2.00 MPa at the intersection of the pier and the bottom plate of the gate. According to the structural mechanics method， the maximum tensile stress at this place is 4.41 MPa， which meets the safety requirements. ②There is a large tensile stress area at the bottom of the middle span of the highway bridge of the drop sluice and part of the bottom plate of the hoist room， which is far beyond the standard value of dynamic axial tensile strength of concrete. For the sake of safety， the corresponding anti-seismic reinforcement measures should be taken. ③The anti-sliding stability safety factor of the sluice chamber structure is greater than the standard value of sluice design code （SL 265-2016）， which meets safety requirements. ④The anti-seismic safety of the drop gate chamber structure meets the requirements of the standard for the seismic design of hydraulic structures （GB 51247-2018）， but there are seismic defects which do not affect the overall safety， and its seismic grade is evaluated as Grade B.

关键词

闸室结构 / 数值模拟 / 振型分解反应谱法 / 抗震安全复核

Key words

chamber structure / numerical simulation / mode decomposition response spectrum method / seismic safety review

基金

黄河水利科学研究院基本科研业务费专项项目(HKY-JBYW-2019-04)

引用本文

EndNote

Ris (Procite)

Bibtex

导出引用

郭博文 , 鲁立三 , 王荆 , 宋力. 基于有限元法的跌水闸闸室结构抗震安全复核研究[J].中国农村水利水电, 2021(7): 155-160

Bo-wen GUO , Li-san LU , Jing WANG , Li SONG. Research on the Seismic Safety Recheck of Drop Gate Chamber Structure Based on Finite Element Method[J].China Rural Water and Hydropower, 2021(7): 155-160

0 引言

水闸是调节水位、控制流量的低水头水工建筑物，具有挡水和泄水（引水）的双重功能，在防洪、治涝、灌溉、供水、航运、发电等方面应用十分广泛。水闸一般由闸室段、上游连接段和下游连接段三部分组成。其中，闸室段是水闸控制水流最主要的部分，通常包括闸墩、闸门、底板、胸墙、岸墙、交通桥和工作桥等几部分。目前，我国具有的各种类型的水闸数量巨大，其中流量大于等于5 m³/s的水闸数量就有105 283 座^［1］。《水闸安全鉴定管理办法》［水建管（2008）214号文］规定，水闸工程在竣工验收后5年内进行首次安全鉴定，以后应每隔10年进行一次全面安全鉴定。

我国是一个多地震灾害国家，抗震安全复核是水闸安全鉴定的重要组成部分；同时，闸室结构作为水闸工程的咽喉部位，对水利枢纽工程安全性的影响尤为突出。因此，有必要对水闸闸室结构抗震安全复核进行研究。目前，大部分学者采用拟静力法对水闸结构进行抗震复核计算，事实上，水闸属于三维空间薄壁结构，如按照常规平面方法简化计算，会忽略闸底板、闸墩、启闭机房以及交通桥等建筑物之间的联系作用，计算结果误差较大^［2］。近年来，随着计算机技术的快速发展，有限元数值模拟技术在水闸抗震分析中得到广泛应用^{［3］-［9］}。因此，有必要采用三维有限元数值模拟技术对闸室结构进行抗震复核研究。

1 研究方法

1.1 振型分解反应谱法

本次采用《水工建筑物抗震设计标准》（GB51247-2018）^［10］（以下简称抗震设计标准）中的标准设计反应谱对闸室结构进行抗震复核计算，具体如图1所示，其中阻尼比取7％，反应谱最大值的代表值β _max取为2.25。

图1 标准设计反应谱

Full size|PPT slide

根据抗震设计标准的规定，采用振型分解反应谱法计算地震作用效应时，可由各阶振型的地震作用效应按平方和方根组合。当两个振型的频率差的绝对值与其中一个较小的频率之比小于0.1时，地震作用效应宜采用完全二次型方根组合，具体计算公式如下：

S_{E} = \sqrt[]{\sum_{i}^{m} \sum_{j}^{m} ρ_{i j}} S_{i} S_{j}

（1）

ρ_{i j} = \frac{8 \sqrt[]{ζ_{i} ζ_{j}} (ζ_{i} + γ_{ω} ζ_{j}) γ_{ω}^{3 / 2}}{(1 - γ_{ω}^{2})^{2} + 4 ζ_{i} ζ_{j} γ_{ω} (1 + γ_{ω}^{2}) + 4 (ζ_{i}^{2} + ζ_{j}^{2}) γ_{ω}^{2}}

（2）

式中：S_E 为地震作用效应；

S_{i}

、

S_{j}

分别为第i阶、第j阶振型的地震作用效应；m为计算采用的振型数；ρ_ij 为第i阶、第j阶振型的振型相关系数；ζ_i 、ζ_j 分别为第i阶、第j阶振型的阻尼比；γ_ω 为圆频率比，γ_ω =ω_j /ω_i；ω_i 、ω_j 分别为第i阶、第j阶振型圆频率。

考虑到地震的往复性，由反应谱法得到的位移、应力等响应指标可正可负。因此，在进行结构响应的动静力叠加时，本次计算采用的叠加原则为：①静力计算结果直接加上反应谱计算结果（简称，静+动）；②静力计算结果直接减去反应谱计算结果（简称，静-动）。

1.2 拉应力复核

目前，如何对有限元计算结果中拉应力超过混凝土轴心抗拉强度标准值区域进行安全评价还研究较少，本文采用有限元数值模拟与结构力学计算相结合的方法对拉应力进行复核。总结开敞式水闸有限元抗震复核结果不难发现，受地震作用影响，闸墩与闸底板相交处易出现较大拉应力区，该区域拉应力一般会超过混凝土动态轴心抗拉强度标准值。地震作用下，闸墩为偏心受压构件，为安全考虑，考虑最不利情况，闸墩按纯弯构件考虑。

依据《水工混凝土结构设计规范》（SL191-2008）^［11］中正截面受弯承载力计算相关内容，矩形截面或翼缘位于受拉边的倒T形截面受弯构件，其正截面受弯承载力应符合下列规定：

K M_{s} \leq f_{c} b x (h_{0} - \frac{x}{2}) + f_{y}^{'} A_{S}^{'} (h_{0} - a_{s}^{'})

（3）

f_{c} b x = f_{y} A_{s} - f_{y} A_{S}^{'}

（4）

（5）

式中：K为承载力安全系数；M_s 为弯矩设计值，N·m；f_c 为混凝土轴心抗压强度设计值，Pa；A_s 为纵向受拉钢筋的截面面积，m²；

A_{S}^{'}

为纵向受压钢筋的截面面积，m²；f_y 为钢筋抗拉强度设计值，Pa；

f_{y}^{'}

为钢筋抗压强度设计值，Pa；h ₀为截面有效高度，m；b为巨型截面的宽度或T形截面的腹板宽度，m；x为受压区计算高度，m；

a_{S}^{'}

为受压钢筋合力点至受压区边缘的距离，m。

同时，对于纯弯构件，其截面上任一点处正应力的计算公式为：

σ = \frac{M y}{I_{z}}

（6）

式中：M为横截面上的弯矩，为得到截面上所能承受的最大正应力，此处M可取为M_s，N·m；I_z 为横截面对中性轴z的惯性矩，m⁴；y为所求应力的点到中性轴z的距离，m。

将根据上述公式计算得到的应力结果与有限元计算结果进行对比，即可对拉应力超过混凝土动态轴心抗拉强度标准值区域进行复核。

1.3 抗滑稳定计算

根据《水闸设计规范》（SL265-2016）^［12］，土基上沿闸室基底面的抗滑稳定安全系数按以下公式计算：

K_{c} = \frac{f \sum G}{\sum H} = \frac{f \sum_{i = 1}^{n} g_{i}}{\sum_{i = 1}^{n} h_{i}}

（7）

式中：K_c 为抗滑稳定安全系数；f为基础底面与地基土之间的摩擦系数；∑G为作用在闸室上的全部竖向荷载，N；∑H为作用在闸室上的全部水平荷载，N；n为闸室底板节点个数；g_i 为第

i

个闸室底板节点上所受的竖向荷载，N；h_i 为第i个闸室底板节点上所受的水平向荷载，N。

2 有限元计算模型

2.1 有限元模型

某跌水闸工程共4孔，为2孔一联钢筋混凝土开敞式水闸，每联之间设置有分缝。闸室段总长度为12.50 m，闸室总宽度为47.45 m，边墩厚度为1.40 m，中墩厚度为1.40 m，缝墩厚度为1.85 m，每孔净宽为10.00 m。考虑到边墩两侧所受土压力和水压力的不同，本次计算主要对边孔一联跌水闸闸室结构进行计算分析。

根据跌水闸闸室结构实际尺寸，建立包括闸底板、闸墩、钢闸门、横梁和启闭机排架结构的边孔三维有限元模型，具体如图2所示，模型共计48 802个节点，34 825个单元，全部采用八结点六面体单元进行空间离散。采用笛卡尔坐标系，横河向为X方向，顺水向方向为Y方向，铅直方向为Z方向，计算时闸墩和闸底板底部施加三向固定约束。需要指出的是，由于该枢纽工程已运行50余年，根据监测资料显示，地基沉降已基本稳定，同时，根据考虑地基时的初步试算结果，静、动力作用下闸室结构闸底板三向位移均较小，最大位移数值不足0.5 mm。因此，为便于分析计算，本次忽略地基的影响。

图2 某跌水闸边孔有限元模型

Full size|PPT slide

另外，由于跌水闸结构每孔闸墩上的公路桥和启闭机房是独立存在的，在相邻两节公路桥和启闭机房之间设置薄层单元，具体如图3所示，在有限元计算中使薄层单元不参与运算，即可实现两节启闭机房之间以及公路桥之间的相对独立性。

图3 跌水闸启闭机房之间以及公路桥之间的薄层单元

Full size|PPT slide

2.2 计算参数

本次计算中没有考虑钢筋单元，所有的混凝土单元代表的是素混凝土，为了体现钢筋的存在对混凝土弹性模量的影响，采用等效弹性模量来模拟钢筋混凝土的弹性模量。其中，在线弹性阶段，钢筋和混凝土是协调变形的，因此等效的原则如下：

E_{d} = E_{c} + E_{s} \frac{A_{s}}{A}

（8）

式中：E_d 为钢筋混凝土材料等效的弹性模量，GPa；E_c 为素混凝土弹性模量，GPa；E_s 为钢筋弹性模量，GPa；A_s 为钢筋截面面积，m²；A为混凝土截面面积，m²。

本次计算所采用的材料参数见表1。

表1 混凝土材料参数

材料号	材料名	密度/（kg·m^-3）	弹模/GPa	泊松比	动态轴心抗压强度标准值/MPa	动态轴心抗拉强度标准值/MPa
1	闸底板	2 548.00	33.23	0.167	19.66	1.97
2	闸墩	2 548.00	31.20	0.167	16.22	1.62
3	公路桥	2 548.00	31.52	0.167	16.72	1.67
4	启闭机房	2 548.00	31.52	0.167	16.72	1.67

注：①表中给出的是各材料参数均为等效后的材料参数；同时，本次计算中各闸室结构不同材料的弹性模量均是采用钻芯取样测得的参数，并根据《水工混凝土结构设计规范》（SL191-2008）^[11]查表计算得到；②根据《水工建筑物抗震设计标准》（GB51247-2018）^[10]中内容规定，进行动力计算时，混凝土材料动态弹性模量可较其静态弹性模量提高50%，其动态轴心抗拉强度标准值可取为动态轴心抗压强度标准值的10%。

2.3 计算工况

本次计算工况为正常蓄水位运行期加Ⅷ度地震作用。其中，正常蓄水位运行时跌水闸闸室结构闸前水深4.03 m，闸后无水；Ⅷ度地震作用时地震动峰值加速度为0.20 g，地基反应谱特征周期T_g 取0.35 s。

2.4 计算荷载

本次计算中主要考虑模型自重、水荷载、淤沙压力、浪压力、土压力、扬压力、风荷载、启闭机自重以及地震荷载等荷载作用。其中，扬压力由浮托力和渗透压力组成，渗透压力可采用改进阻力系数法计算。

3 有限元计算结果与分析

3.1 自振特性计算结果与分析

运用结构自振特性分析方法，考虑闸前水体对闸室结构的影响，对结构进行自振分析，得到跌水闸闸室结构前五阶频率、振型等自振特性参数，其中闸前水体对结构的影响采用Westergard的附加质量法进行模拟。正常蓄水位工况下跌水闸闸室结构前五阶自振频率和周期如表2所示。从表2中可以看出，在正常蓄水工况下跌水闸闸室结构自振的基频为4.038 Hz，且第二阶自振频率与基频数值相近，这主要是因为闸室结构每孔闸墩上的启闭机房是相互独立的，而闸室结构第一阶和第二阶振型分别为两节相互独立的启闭机房沿顺河向方向的振动。

表2 不同工况下跌水闸闸室结构自振频率周期表

阶数	正常蓄水位工况
阶数	频率/Hz	周期/s
1	4.038	0.248
2	4.178	0.239
3	7.708	0.130
4	13.013	0.077
5	13.820	0.072

针对正常蓄水位工况，图4给出了边孔闸室结构前五阶振型图。可以看出，跌水闸闸室结构前五阶振型主要体现为启闭机房和公路桥的振动，闸墩参与较少。

图4 跌水闸边孔闸室结构各阶振型图

Full size|PPT slide

3.2 位移计算结果与分析

图5至图6给出了动静叠加下边孔闸室结构在X向、Y向和Z向位移等值线图。由图5、6可知，正常蓄水工况地震作用下跌水闸闸室结构位移在不同的叠加方式下最大数值分别为12.0和10.0 mm，位置出现在启闭机房顶部位置，主要为顺河向位移。

图5 边孔闸室结构动静叠加（静+动）下位移等值线图（单位：m）

Full size|PPT slide

图6 边孔闸室结构动静叠加（静-动）下位移等值线图（单位：m）

Full size|PPT slide

3.3 应力计算结果与分析

图7至图9分别给出了动静叠加下跌水闸闸室结构各部位应力计算结果等值线图。可以看出，在正常蓄水位地震作用下，受地震作用影响，跌水闸闸墩与闸底板相交处出现了2.00 MPa左右的拉应力，超过了跌水闸闸墩和闸底板混凝土动态轴心抗拉强度标准值（具体数值参照表1）。图10给出了跌水闸边墩单位长度截面配筋图，基于2.2节中所介绍的拉应力复核方法，根据式（1）至式（3）计算可得单位长度边墩横截面上最大承载弯矩为2 879.6 kN·m，由式（4）计算得到与跌水闸边墩与闸底板相交处能承受的最大拉应力为4.41 MPa，大于2.00 MPa，满足安全需求。

受自身重量、公路桥荷载以及地震作用影响，跌水闸公路桥跨中底部位置处出现了较大拉应力区，其数值达到了5.00 MPa左右，远超跌水闸公路桥混凝土动态轴心抗拉强度标准值，考虑到跌水闸公路桥运行时间较长，且桥梁底部钢筋混凝土也出现脱落、露筋和胀裂等情况，具体如图11所示，应对跌水闸公路桥进行相应的抗震处理措施。另外，受地震作用影响，启闭机底板部分区域出现了8.00 MPa拉应力区，为安全起见，宜采取相应的抗震加固措施。

图11 跌水闸桥梁底部混凝土脱落、漏筋和胀裂现象

Full size|PPT slide

同时，机公路桥和启闭机房与闸墩联结位置处有较大压应力区，最大压应力数值达到了9.00 MPa左右，未超过跌水闸闸墩、公路桥和启闭机房混凝土动态轴心抗压强度标准值，满足安全要求。

3.4 稳定计算结果与分析

表3给出了动静叠加下跌水闸闸室结构抗滑稳定计算结果。由表3可知，根据有限元动力计算结果，作用在边孔闸室上的水平向地震荷载为2 863.72 kN。由于地震作用是随机往复的，当该水平向地震惯性力朝下游时，结合静力计算结果，此时跌水闸闸室结构抗滑稳定安全系数K_c 为2.30，满足安全需求；当该水平向地震惯性力朝上游时，结合静力计算结果，此时跌水闸闸室结构抗滑稳定安全系数K_c 为4.64，满足安全需求。

表3 动静叠加下跌水闸闸室结构抗滑稳定计算分析表

工况	静力工况下竖向荷载/kN	静力工况下水平向荷载/kN	水平向地震惯性力/kN	摩擦系数	抗滑稳定系数	规范值
水平向地震惯性力朝上游	25 197.52	964.93	-2 863.72	0.35	4.64	1.10
水平向地震惯性力朝下游	25 197.52	964.93	2 863.72	0.35	2.30	1.10

3.5 抗震安全复核

综合上述计算结果，根据《水闸安全评价导则》（SL214-2015）^［13］，跌水闸抗震安全满足《水工建筑物抗震设计标准》（GB51247-2018）要求，但存在不影响总体安全的缺陷，其抗震等级评定为B级。

4 结论

基于有限元数值模拟技术，对某跌水闸闸室进行了抗震复核计算，根据计算结果，得到如下结论。

（1）正常蓄水位工况下，跌水闸闸室结构基频为3.57 Hz，其第一阶振型表现为启闭机房沿顺河向振动；

（2）跌水闸闸墩与闸底板相交处出现了2.00 MPa拉应力，结合此处的配筋量，该处能承受的最大拉应力为4.41 MPa，满足安全需求；同时，跌水闸公路桥跨中底部位置出现了较大拉应力区，其数值远超混凝土动态轴心抗拉强度标准值，且桥梁底部钢筋混凝土出现了脱落、露筋和胀裂等情况，为安全起见，宜采取相应的抗震加固措施；

（3）地震作用最不利工况下，跌水闸闸室结构抗滑稳定安全系数为2.30，满足安全需求。根据《水闸安全评价导则》（SL214-2015），跌水闸抗震安全满足《水工建筑物抗震设计标准》（GB51247-2018）要求，但存在不影响总体安全的缺陷，其抗震等级评定为B级。

参考文献

原文顺序 | 文献年度倒序 | 文中引用次数倒序

1	宋力，汪自力. 水闸安全评价及加固修复技术指南［M］. 郑州：黄河水利出版社，2018．本文引用 [1]

2	马利嘉，蔡荨，历丹丹. 基于有限元法的闸室抗震安全评价［J］. 水利水电技术，2017，48（3）：30-33，70. 本文引用 [1]

3	曾迪，邸庆霜. 水闸震害调查与抗震研究进展［J］. 水利水电技术， 2014，45（10）：42-44. 本文引用 [1]

4	卢羽平. 水闸参数化三维有限元分析软件的研究［D］. 南京：河海大学， 2005.

5	孙瑞祥. 高层建筑结构竖向地震作用下振型分解反应谱法分析研究［D］. 兰州：兰州理工大学， 2011.

6	胡晓. 水闸震害调查与分析［J］. 水力发电， 2009，35（5）：18-20.

7	王雪奎. 黄河中下游水闸抗震能力分析［D］. 郑州：华北水利水电大学， 2016.

8	许新勇，康迎宾，罗全胜. 洹河景观水闸工程抗震安全评价研究［J］. 人民黄河，2012，34（ 8）：147-148.

9	朱庆华，顾美娟. 水闸闸室抗震动力分析及措施［J］. 水电能源科学，2012，30（1）：114-116，208. 本文引用 [1]

10	中华人民共和国住房和城乡建设部等. 水工建筑物抗震设计标准：GB51247-2018 ［S］. 北京：中国计划出版社，2018. 本文引用 [2]

11	中华人民共和国水利部. 水工混凝土结构设计规范：SL191-2008 ［S］. 北京：中国水利水电出版社，2008. 本文引用 [2]

12	中华人民共和国水利部. 水闸设计规范：SL265-2016 ［S］. 北京：中国水利水电出版社，2016. 本文引用 [1]

13	中华人民共和国水利部. 水闸安全评价导则：SL214-2015 ［S］. 北京：中国水利水电出版社，2015. 本文引用 [1]

PDF(2124 KB)

访问

引用

详细情况

段落导航

摘要
Abstract
关键词
Key words
基金
引用本文
0 引言
1 研究方法
1.1 振型分解反应谱法
图1 标准设计反应谱
1.2 拉应力复核
1.3 抗滑稳定计算
2 有限元计算模型
2.1 有限元模型
图2 某跌水闸边孔有限元模型
图3 跌水闸启闭机房之间以及公路桥之间的薄层单元
2.2 计算参数
表1 混凝土材料参数
2.3 计算工况
2.4 计算荷载
3 有限元计算结果与分析
3.1 自振特性计算结果与分析
表2 不同工况下跌水闸闸室结构自振频率周期表
图4 跌水闸边孔闸室结构各阶振型图
3.2 位移计算结果与分析
图5 边孔闸室结构动静叠加（静+动）下位移等值线图（单位：m）
图6 边孔闸室结构动静叠加（静-动）下位移等值线图（单位：m）
3.3 应力计算结果与分析
图7 跌水闸闸墩和底板结构主应力等值线图（单位：Pa）
图8 横梁动静叠加下主应力等值线图（单位：Pa）
图9 机架桥动静叠加下主应力等值线图（单位：Pa）
图10 跌水闸边墩单位长度截面配筋图
图11 跌水闸桥梁底部混凝土脱落、漏筋和胀裂现象
3.4 稳定计算结果与分析
表3 动静叠加下跌水闸闸室结构抗滑稳定计算分析表
3.5 抗震安全复核
4 结论
参考文献

收稿日期	出版日期
2020-09-22	2021-07-15
发布日期
2021-08-22

选择文件类型/文献管理软件名称

选择包含的内容

摘要

Abstract

关键词

Key words

基金

引用本文

0 引 言

1 研究方法

1.1 振型分解反应谱法

图1 标准设计反应谱

1.2 拉应力复核

1.3 抗滑稳定计算

2 有限元计算模型

2.1 有限元模型

图2 某跌水闸边孔有限元模型

图3 跌水闸启闭机房之间以及公路桥之间的薄层单元

2.2 计算参数

表1 混凝土材料参数

2.3 计算工况

2.4 计算荷载

3 有限元计算结果与分析

3.1 自振特性计算结果与分析

表2 不同工况下跌水闸闸室结构自振频率周期表

图4 跌水闸边孔闸室结构各阶振型图

3.2 位移计算结果与分析

图5 边孔闸室结构动静叠加（静+动）下位移等值线图（单位：m）

图6 边孔闸室结构动静叠加（静-动）下位移等值线图（单位：m）

3.3 应力计算结果与分析

图7 跌水闸闸墩和底板结构主应力等值线图（单位：Pa）

图8 横梁动静叠加下主应力等值线图（单位：Pa）

图9 机架桥动静叠加下主应力等值线图（单位：Pa）

图10 跌水闸边墩单位长度截面配筋图

图11 跌水闸桥梁底部混凝土脱落、漏筋和胀裂现象

3.4 稳定计算结果与分析

表3 动静叠加下跌水闸闸室结构抗滑稳定计算分析表

3.5 抗震安全复核

4 结 论

{{custom_sec.title}}

{{custom_sec.title}}

参考文献

0 引言

4 结论