基于ANSYS的重力坝三维静动态结构分析
目录
1 引言 (1)
2 工程概况 (1)
3 基本资料 (1)
3.1 反应谱 (1)
3.2 材料参数 (2)
3.3 规范要求 (2)
4 分析简介 (4)
4.1 分析模型 (4)
4.2 边界条件 (6)
4.3 荷载工况 (6)
5 计算成果 (7)
5.1 工况一 (7)
5.2 工况二 (8)
5.3 工况三 (10)
5.4 工况四 (11)
5.5 工况五 (12)
5.6 工况六 (14)
5.7 结果总结及分析 (15)
6 结论及建议 (17)
7 分析命令流 (17)
1 引言
重力坝是我国高坝中的主要坝型，在防洪、发电、灌溉、城镇供水、航运、养殖和旅游等方面发挥了巨大的作用，取得了显著的经济效益和社会效益。众所周知，重力坝主要依靠其自身的重力来维持稳定，其坝体体积大，稳定性好。但由于各种原因，仍有可能失事。因此，重力坝的应力应变状态和坝基稳定性一直都是设计和施工十分重视的问题。此外，大坝多建于地震频发的地区，因而对重力坝进行地震荷载作用下的安全评估也十分必要。本次作业采用有限元方法，运用大型通用有限元分析软件ANSYS，对简化的三维重力坝的线弹性模型在静动力工况下进行有限元计算，并对结果加以分析，最后给出安全评价结论及建议。
2 工程概况
某水电站是以发电为主，兼有防洪，航运等综合效益的水电枢纽工程。该工程枢纽总体布置采用砼重力坝挡水，大坝基本坝剖面为上游坝坡铅直，下游坝坡为1:0.75。坝顶总长270m，坝高180m，坝顶宽18m，坝底宽139.5m，正常蓄水位170m。重力坝坝低至坝高100m之间使用坝体混凝土Ⅱ，坝高100m至坝顶之间使用坝体混凝土Ⅰ。上游正常蓄水位为170m，下游无水。
3 基本资料
3.1 反应谱
谱分析是一种将模态分析的结果与一个已知的谱联系起来计算结构的位移和应力的分析技术。在土木工程动力响应分析中，谱分析代替时间-历程分析，特别是抗震分析，主要用来确定结构对随机荷载或随时间变化荷载的动力响应。
图3.1 大坝设计反应谱
根据如图3.1所示的大坝设计的反应谱曲线图，可得大坝反应谱曲线方程：
本次重力坝抗震性能分析中，max β取值为2，g T 取值为0.3。本次谱分析采用的SV-FREQ 曲线谱值点见表3.1。
表3.1 SV-FREQ 曲线谱值点
3.2 材料参数
该重力坝坝体采用两种混凝土材料，具体材料参数见表3.2。其中混凝土动力分析弹性模量是静力分析弹性模量的1.5倍。计算时考虑为完全沉降后作用，故不考虑基岩密度。
表3.2 坝体混凝土材料属性
3.3 规范要求
如未特殊说明，本次分析参照的规范均指《混凝土重力坝设计规范》 NB/T 35026-2014，以下为本次分析需用到的部分重要章节。 1）正常使用极限状态
按材料力学方法进行坝体上、下游面混凝土拉应力验算；必要时进行坝体及结构变形计算、复杂地基局部渗透稳定验算。
正常使用极限状态作用效应采用下列设计表达式：
(),,,0K K K k S G Q f a C
γ≤
对正常使用极限状态验算时，作用分项系数、材料性能分项系数都取1.0 ，结构重要性系数不变。
该规范提出了重力坝对正常使用极限状态的要求，规定坝踵及坝体上游面不产生垂直拉应力，施工期坝趾处垂直正应力可容许有不大于0.1MPa 的拉应力，下游坝面主拉应力不大于0.2MPa 2）承载能力极限状态
承载能力极限状态，对坝体结构及坝基岩体进行强度和抗滑稳定计算，必要时进行抗浮、抗倾验算:抗震设防应满足DL5073 的有关规定。
(),,,01
K G K Q K K K d m f S G Q a R a γψγγγγ⎛⎫
≤
⎪⎝⎭ 作用分项系数：0.95G γ= 1Q γ= 1.1u γ= 结构重要性系数：01γ=
设计状况系数：1ψ= 0.85（地震） 结构系数： 1.8d γ= 1.5（地震） 材料性能分项系数： 1.5m γ= 3）抗滑稳定验算
坝体混凝土与基岩接触面的抗滑稳定极限状态。 作用效应函数为：
()R S P •=∑
抗滑稳定抗力函数为：
()R
R R R R f W C A ''•=∑+ 式中：R P ∑----坝基面上全部切向应力作用之和（KN ）； R f '-------坝基面抗剪断摩擦系数；
R
C '---- --坝基面抗剪断凝聚力(KPa)。 本次分析系数取值：
Rk R mf
f f γ'
'=
1.7mf γ= 1.1Rk f '=
Rk R mc
C C γ''=
2.0mc γ= 1.05Rk C MPa '= 4）混凝土强度
本次分析的大坝坝体上部混凝土为C20，坝体下部为C25混凝土。用于混凝土重力坝承载能力极限状态计算的坝体混凝土强度标准值按表3.3取用。
表3.3 大坝混凝土强度标准值
4 分析简介 4.1 分析模型
坐标系定义为：X 轴以水平向下游为正，Y 轴以铅直向上为正，原点取在坝踵处，有限元计算范围：上游取1.5倍坝高，即上游长度为270m ，下游取1.5倍坝高，即下游长度为坝轴线向下270m ，铅直取2倍坝高，即向下取360m 基岩深度。
图4.1 重力坝分析模型剖面图
图4.2 重力坝整体有限元模型
