1、基于ANSYS的重力坝三维静动态结构分析基于ANSYS的重力坝三维静动态结构分析目 录1 引言 2 工程概况 3 基本资料 3.1 反应谱 3.2 材料参数 3.3 规范要求 4 分析简介 4.1 分析模型 4.2 边界条件 4.3 荷载工况 5 计算成果 5.1 工况一 5.2 工况二 5.3 工况三 5.4 工况四 5.5 工况五 5.6 工况六 5.7 结果总结及分析 6 结论及建议 7 分析命令流 1 引言重力坝是我国高坝中的主要坝型，在防洪、发电、灌溉、城镇供水、航运、养殖和旅游等方面发挥了巨大的作用，取得了显著的经济效益和社会效益。众所周知，重力坝主要依靠其自身的重力来维持稳定，其
2、坝体体积大，稳定性好。但由于各种原因，仍有可能失事。因此，重力坝的应力应变状态和坝基稳定性一直都是设计和施工十分重视的问题。此外，大坝多建于地震频发的地区，因而对重力坝进行地震荷载作用下的安全评估也十分必要。本次作业采用有限元方法，运用大型通用有限元分析软件ANSYS，对简化的三维重力坝的线弹性模型在静动力工况下进行有限元计算，并对结果加以分析，最后给出安全评价结论及建议。2 工程概况某水电站是以发电为主，兼有防洪，航运等综合效益的水电枢纽工程。该工程枢纽总体布置采用砼重力坝挡水，大坝基本坝剖面为上游坝坡铅直，下游坝坡为1:0.75。坝顶总长270m，坝高180m，坝顶宽18m，坝底宽139.
3、5m，正常蓄水位170m。重力坝坝低至坝高100m之间使用坝体混凝土，坝高100m至坝顶之间使用坝体混凝土。上游正常蓄水位为170m，下游无水。3 基本资料3.1 反应谱谱分析是一种将模态分析的结果与一个已知的谱联系起来计算结构的位移和应力的分析技术。在土木工程动力响应分析中，谱分析代替时间-历程分析，特别是抗震分析，主要用来确定结构对随机荷载或随时间变化荷载的动力响应。图3.1 大坝设计反应谱根据如图3.1所示的大坝设计的反应谱曲线图，可得大坝反应谱曲线方程：本次重力坝抗震性能分析中，取值为2，取值为0.3。本次谱分析采用的SV-FREQ曲线谱值点见表3.1。表3.1 SV-FREQ曲线谱值
4、点序号周期频率谱值101e5120.011001.130.110240.33.33333333250.42.51.6829327260.61.666666671.3195079170.81.251.110321528110.9711867591.20.833333330.87055056101.40.714285710.79364435111.60.6250.732539011220.50.64074431330.333333330.502377293.2 材料参数该重力坝坝体采用两种混凝土材料，具体材料参数见表3.2。其中混凝土动力分析弹性模量是静力分析弹性模量的1.5倍。计算时考虑为完全沉
5、降后作用，故不考虑基岩密度。表3.2 坝体混凝土材料属性参数混凝土材料编号密度（kg/m3）弹性模量（Pa）泊松比坝体混凝土124002.6E100.167坝体混凝土224002.85E100.167基岩302.9E100.33.3 规范要求如未特殊说明，本次分析参照的规范均指混凝土重力坝设计规范NB/T 35026-2014，以下为本次分析需用到的部分重要章节。1）正常使用极限状态按材料力学方法进行坝体上、下游面混凝土拉应力验算；必要时进行坝体及结构变形计算、复杂地基局部渗透稳定验算。正常使用极限状态作用效应采用下列设计表达式：对正常使用极限状态验算时，作用分项系数、材料性能分项系数都取1.
6、0 ，结构重要性系数不变。该规范提出了重力坝对正常使用极限状态的要求，规定坝踵及坝体上游面不产生垂直拉应力，施工期坝趾处垂直正应力可容许有不大于0.1MPa的拉应力，下游坝面主拉应力不大于0.2MPa2）承载能力极限状态承载能力极限状态，对坝体结构及坝基岩体进行强度和抗滑稳定计算，必要时进行抗浮、抗倾验算:抗震设防应满足DL5073 的有关规定。作用分项系数： 结构重要性系数： 设计状况系数： 0.85（地震）结构系数： 1.5（地震）材料性能分项系数： 3）抗滑稳定验算坝体混凝土与基岩接触面的抗滑稳定极限状态。作用效应函数为：抗滑稳定抗力函数为：式中： -坝基面上全部切向应力作用之和（KN）
7、； -坝基面抗剪断摩擦系数； - -坝基面抗剪断凝聚力(KPa)。本次分析系数取值： 4）混凝土强度本次分析的大坝坝体上部混凝土为C20，坝体下部为C25混凝土。用于混凝土重力坝承载能力极限状态计算的坝体混凝土强度标准值按表3.3取用。表3.3 大坝混凝土强度标准值4 分析简介4.1 分析模型坐标系定义为：X轴以水平向下游为正，Y轴以铅直向上为正，原点取在坝踵处，有限元计算范围：上游取1.5倍坝高，即上游长度为270m，下游取1.5倍坝高，即下游长度为坝轴线向下270m，铅直取2倍坝高，即向下取360m基岩深度。图4.1 重力坝分析模型剖面图图4.2 重力坝整体有限元模型图4.3 基岩有限元模
8、型 图4.4 坝体有限元模型图4.5 上部坝体有限元模型 图4.6 下部坝体有限元模型表4.1 有限元模型具体数据参数重力坝单元数节点数材料编号单元类型模型整体18560217141、2、3Solid45、Solid65坝基16800197343Solid45坝体176023461、2Solid65坝体上部88012241坝体下部880122424.2 边界条件基岩底部采用全约束，上下游坝基采用法向约束，即将X方向的两个侧面的X向位移进行约束，坝基两侧采用法向约束，即将Z方向的两个侧面的Z向位移进行约束，上部边界为自由边界，假设坝体与坝基固结。4.3 荷载工况(1) 自重：混凝土材料参数见3.
9、2节；(2) 静水压力：水容重取9.81KN/m3；(3) 扬压力： 扬压力计算简化坝踵处为无排水孔。(4)地震作用： 反应谱见3.1节。根据荷载的不同分为六种工况，工况一工况三为静力分析，工况四工况六为动力分析。表4.2 分析工况一览表荷载工况自重静水压力扬压力地震荷载备注工况一库空工况二正常使用极限状态工况三0.9511.1承载能力极限状态工况四模态分析工况五谱分析工况六0.9511.1偶然荷载注：表格中“”表示工况内含有该荷载，表中系数表示该荷载的作用分项系数；静水压力、扬压力按坝面分布荷载施加。5 计算成果5.1 工况一1）计算结果图5.1 整体模型总位移（m）图 图5.2 坝体UX（
10、m）图图5.3 坝体UY（m）图 图5.4 坝体S1（Pa）图2）本节小结表5.1 位移结果汇总表 结果位置 X方向位移(cm)Y方向位移(cm)总位移(cm)坝顶-1.54-2.022.54坝踵0.04-0.86坝趾0.04-0.28注：表中数据均指最大值，其中总位移为整体模型的总位移最大值。X方向为负值表明坝体倾向上游，Y方向位移为负值时表明坝体沉降，下同。表5.2 坝体应力结果汇总表 位置应力坝踵坝趾最大拉应力最大压应力S1（MPa）-0.670.203.27-0.67注：表中数据指该位置处的最大值，下同。库空工况下，坝体重心偏向上游，坝体整体向上游倾斜，由此产生的位移矢量总体趋势是倾向
11、上游，见图5.1。坝体水平变位从建基面至坝顶，随高程增加而加大，坝顶水平位移最大，其值为1.54cm(向上游)，见图5.2。对同一高程而言，坝体上游侧垂直变位大于下游侧垂直向变位，坝顶铅直向变位最大，其值为2.02cm(沉降)，坝踵铅直向变位0.86cm（沉降），坝趾铅直向变位0.28cm（沉降）。对于坝体而言，其最大主应力为拉应力，坝体最大主应力（图5.4）随越靠近两侧坝肩越大，同时，在坝肩附近有应力集中现象，最大拉应力为3.27MPa。5.2 工况二1）计算结果图5.5 整体模型总位移（m）图 图5.6 坝体UX（m）图图5.7 坝体UY（m）图 图5.8 坝体S1（Pa）图2）本节小结表
12、5.3 位移结果汇总表 结果位置 X方向位移(cm)Y方向位移(cm)总位移(cm)坝顶1.13-1.061.53坝踵0.36-0.28坝趾0.25-0.47表5.4 坝体应力结果汇总表 位置应力坝踵坝趾最大拉应力最大压应力S1（MPa）-0.230.172.57-1.03正常使用极限状态下，从图5.5中可以看出，坝体运行期间由于受库水压力+扬压力作用，产生的位移矢量总体趋势是偏向下游，坝体水平向变位倾向下游，垂直沉降变位量值在同一高程下游坝面大于上游坝面（图5.7）。正常蓄水位工况坝顶水平向位移为1.13cm（向下游），铅直向变位值为1.06cm（沉降）。由于坝基扬压力作用，沉降量值较库空工
13、况明显减小。坝体最大主应力为拉应力，就同一高程而言，上游侧最大主压应力大于下游侧（图5.8）。在坝肩靠近上游侧附近有应力集中现象，最大拉应力为2.57MPa。5.3 工况三1）计算结果图5.9 整体模型总位移（m）图 图5.10 坝体UX（m）图图5.11 坝体UY（m）图 图5.12 坝体S1（Pa）图2）本节小结表5.5 位移结果汇总表 结果位置 X方向位移(cm)Y方向位移(cm)总位移(cm)坝顶1.19-0.931.50坝踵0.37-0.13坝趾0.25-0.39表5.6 坝体应力结果汇总表 位置应力坝踵坝趾最大拉应力最大压应力S1（MPa）-0.230.162.54-1.02承载能
14、力极限状态下，坝体整体变位规律和应力分布规律与工况二的正常使用极限状态一致。但该状态重力、静水压力、扬压力分别乘以相应的作用分项系数后，坝体铅直向沉降量值较正常使用极限状态明显减小。5.4 工况四1）计算结果图5.13 坝体第一阶振型 图5.14 坝体第二阶振型图5.15 坝体第三阶振型 图5.16 坝体第四阶振型图5.17 坝体第五阶振型 图5.18 坝体第六阶振型2）本节小结表5.7 模态分析前10阶自振频率* INDEX OF DATA SETS ON RESULTS FILE *SETTIME/FREQLOAD STEPSUBSTEPCUMULATIVE13.364411125.583
15、812235.945113346.125214456.419015568.226916678.331417788.611418899.73991991011.1980110105.5 工况五1）计算结果图5.19 整体模型总位移（m）图 图5.20 坝体UX（m）图图5.21 坝体UY（m）图 图5.22 坝体S1（Pa）图2）本节小结表5.8 位移结果汇总表 结果位置 X方向位移(cm)Y方向位移(cm)总位移(cm)坝顶2.750.982.92坝踵0.040.14坝趾0.040.05表5.9 坝体应力结果汇总表 位置应力坝踵坝趾最大拉应力最大压应力S1（MPa）2.291.264.86从图
16、5.19中可以看出，坝体位移矢量总体趋势是偏向下，坝体水平向变位倾向下游，坝体水平变位从建基面至坝顶，随高程增加而加大，坝顶水平位移最大，其值为2.75cm(向上游)，见图5.20。对同一高程而言，坝体上游侧垂直变位大于下游侧垂直向变位，坝顶铅直向变位最大（图5.21）。对于坝体而言，最大主应力为拉应力，在坝肩靠近下游侧附近有应力集中现象最大拉应力为4.86MPa。5.6 工况六1）计算结果图5.23 整体模型总位移（m）图 图5.24 坝体UX（m）图图5.25 坝体UY（m）图 图5.26 坝体S1（Pa）图2）本节小结表5.10 位移结果汇总表 结果位置 X方向位移(cm)Y方向位移(c
17、m)总位移(cm)坝顶3.90-0.503.93坝踵0.20-0.05坝趾0.20-0.34表5.11 坝体应力结果汇总表 位置应力坝踵坝趾最大拉应力最大压应力S1（MPa）1.011.015.72-0.33自重+静水压力+扬压力+地震荷载条件下，坝体位移矢量总体趋势是倾向下游，见图5.23。坝体水平变位从建基面至坝顶，随高程增加而加大，坝顶水平位移最大，其值为3.90cm(向上游)，见图5.24。对同一高程而言，坝体上游侧垂直变位小于下游侧垂直向变位，坝体下游坡面中间位置铅直向变位最大。对于坝体而言，其最大主应力为拉应力，在坝肩附近有应力集中现象，最大拉应力为5.72MPa。5.7 结果总结
18、及分析表5.12 各工况位移结果一览表 结果工况 X方向位移(cm)Y方向位移(cm)总位移(cm)工况一坝顶-1.54-2.022.54坝踵0.04-0.86坝趾0.04-0.28工况二坝顶1.13-1.061.53坝踵0.36-0.28坝趾0.25-0.47工况三坝顶1.19-0.931.50坝踵0.37-0.13坝趾0.25-0.39工况五坝顶2.750.982.92坝踵0.040.14坝趾0.040.05工况六坝顶3.90-0.503.93坝踵0.20-0.05坝趾0.20-0.34表5.13 坝体第一主应力S1（MPa）结果一览表 位置工况坝踵坝趾最大拉应力最大压应力工况一-0.67
19、0.203.27-0.67工况二-0.230.172.57-1.03工况三-0.230.162.54-1.02工况五2.291.264.86工况六1.011.015.72-0.331）坝体变位分析坝体变位规律符合常规，库空条件下，坝体重心偏向上游，坝体整体向上游倾斜，由此产生的位移矢量总体趋势是倾向上游；正常使用极限状态下，从图5.5中可以看出，坝体运行期间由于受库水压力+扬压力作用，产生的位移矢量总体趋势是偏向下游，坝体水平向变位倾向下游，垂直沉降变位量值在同一高程下游坝面大于上游坝面。2）正常使用极限状态强度分析由图5.5和表5.13可知，坝体坝踵和坝趾处未产生拉应力，符合规范要求，但坝体
20、与两岸连接处坝肩的上游侧附近第一主拉应力达到2.57MPa1.78MPa，已超过混凝土的强度标准值，发生破坏，应采取加固措施。3）承载能力极限状态强度分析坝体材料C20，轴心抗压强度，轴心抗拉强度。 静力工况根据3.3节的公式和系数可得， 。即坝体压应力满足强度要求，拉应力不符合规范要求。 地震工况根据3.3节的公式和系数可得， 。即坝体压应力满足强度要求，拉应力不符合规范要求。4）抗滑稳定验算 使用ANSYS映射面功能（具体操作见命令流）算得：，（方向向下），根据3.3节的系数求出所以该重力坝满足稳定性要求。6 结论及建议本文使用有限元分析软件ANSYS对坝体进行了静力分析和动力分析，通过以
21、上的计算和分析，可以得出以下结论。(1) 重力坝变位规律符合常规，重力坝水平方向的最大位移发生在坝顶，地震荷载作用下，坝顶水平位移更大。对此，为了提高该重力坝的抗震性，应采取措施减小坝顶水平变位，如可以适当增大顶部刚度，减轻重量，提高混凝土标号，必要时可以在上下游面布设钢筋。（2）正常使用极限状态和承载能力极限状态下，在静力工况下，坝体的坝踵和坝趾位置处的应力值满足规范要求；在地震工况下，坝踵处出现1.01MPa的拉应力，可在此处适当增加配筋，防止地震之后裂缝发生和扩展；确保混凝土的施工质量，加强温控措施和养护，尽量减少表面裂缝。并且坝体的两岸坝肩处的拉应力较大，远超出混凝土的强度标准值，应采
22、取加固措施。(3) 通过稳定验算可知，该重力坝能够满足抗滑稳定性要求。为进一步提高该重力坝的抗滑稳定性，建议在靠近坝体上游面，采用深孔锚固高强度钢索，并施加预应力，既可增加坝体的抗滑稳定，又可消除坝踵处的拉应力。7 分析命令流本次分析以APDL语言为主，GUI操作为辅，主要的操作均编写成命令流形式（见下），但其运行时需注意其个别参数的调整（如作用分项系数），荷载工况组合主要以GUI操作为主，故此处未附该操作的命令流。FINISH$/CLEAR/FILNAME,Gravity dam$/PREP7!一.静力分析ET,1,SOLID65$ET,2,SOLID45$ET,3,MESH200$KEYO
23、PT,3,1,6 !静力分析材料属性MP,DENS,1,2400$MP,EX,1,2.6E10$MP,PRXY,1,0.167MP,DENS,2,2400$MP,EX,2,2.85E10$MP,PRXY,2,0.167MP,DENS,3,0$MP,EX,3,2.9E10$MP,PRXY,3,0.3!重力坝建模K,1$K,2,139.5$K,3,0,162$K,4,18,162$K,5,18,180$K,6,0,180$K,7,-270,-360K,8,0,-360$K,9,-270,$K,10,139.5,-360$K,11,409.5,-360$K,12,409.5$K,13,409.5,1
24、80K,14,-270,180A,3,4,5,6$A,1,2,4,3$A,1,3,6,14,9$A,2,12,13,5,4$A,1,9,7,8$A,1,2,10,8$A,2,12,11,10LSEL,S,2,4,2$CM,LN4,LINE$ALLSEL$LESIZE,LN4,4LSEL,S,1,5,2$LSEL,A,18$CM,LN5,LINE$ALLSEL$LESIZE,LN5,5LSEL,S,6,7$CM,LN18,LINE$ALLSEL$LESIZE,LN18,18LSEL,S,9,12,3$CM,LN22,LINE$ALLSEL$LESIZE,LN22,22LESIZE,8,10,2$
25、LESIZE,13,10,0.5$LESIZE,10,10,0.5$LESIZE,11,10,2LESIZE,14,10,2$LESIZE,16,10,0.5$LESIZE,17,10,2$LESIZE,19,10,2LESIZE,15,10,0.5$LESIZE,20,10,0.5LCCAT,2,6$LCCAT,4,7AATT,3$MSHKEY,1$AMESH,ALL$LDELE,21,22 $EPLOTTYPE,1$EXTOPT,ACLEAR,1$EXTOPT,ESIZE,8$EXTOPT,ATTR,0MAT,1$ESYS,0$VEXT,1,2,180*0.75TYPE,2$EXTOPT,
26、ACLEAR,1$EXTOPT,ESIZE,8$EXTOPT,ATTR,0MAT,3$ESYS,0$VEXT,5,7,1,180*0.75TYPE,2$EXTOPT,ACLEAR,1$EXTOPT,ESIZE,8$EXTOPT,ATTR,0MAT,3$ESYS,0$VEXT,1,7,1,-180WPOFFS,180*0.75$CSYS,4$VSYMM,Z,ALL$WPCSYS,-1$CSYS,0NUMMRG,ALL$NUMCMP,ALL!更改材料属性ESEL,S,TYPE,1$NSLE,R,ALL$NSEL,R,LOC,Y,0,100$ESLN,S,ALLEMODIF,ALL,MAT,2$AL
27、LSEL$EPLOTWPCSYS,-1$WPSTYLE,0$EPLOT!施加约束NSEL,S,LOC,Y,-360 !底部全约束的点CM,Y1,NODE$ALLSEL NSEL,S,LOC,X,409.5 !上下游X方向UX法向约束NSEL,A,LOC,X,-270$CM,X2,NODE$ALLSELNSEL,S,LOC,Z,450 !坝肩两侧Z方向UZ法向约束NSEL,A,LOC,Z,-180$CM,Z3,NODE$ALLSELD,Y1,ALL$D,X2,UX$D,Z3,UZ!施加荷载ACEL,9.81 !自重/PSF,PRES,NORM,2SFGRAD,PRES,Y,170,-9810 !施加静水压力ESEL,S,MAT,1,2$NSLE,S,ALL$NSEL,R,LOC,X,0$NSEL,R,LOC,Y,0,170CM,SHUI,NODE$SF,SHUI,PRES,0$ALLSELSFGRAD,PRES,X,139.5,-1.1*9.81*1000*170/139.5 !施加扬压力ESEL,S,MAT,1,2$NSLE,S,ALL$NS
