基于ANSYS的重力坝三维静动态结构分析Word文档下载推荐.docx
- 文档编号:6176328
- 上传时间:2023-05-06
- 格式:DOCX
- 页数:29
- 大小:1.71MB
基于ANSYS的重力坝三维静动态结构分析Word文档下载推荐.docx
《基于ANSYS的重力坝三维静动态结构分析Word文档下载推荐.docx》由会员分享,可在线阅读,更多相关《基于ANSYS的重力坝三维静动态结构分析Word文档下载推荐.docx(29页珍藏版)》请在冰点文库上搜索。
取值为2,
取值为0.3。
本次谱分析采用的SV-FREQ曲线谱值点见表3.1。
表3.1SV-FREQ曲线谱值点
序号
周期
频率
谱值
1
1e5
2
0.01
100
1.1
3
0.1
10
4
0.3
3.33333333
5
0.4
2.5
1.68293272
6
0.6
1.66666667
1.31950791
7
0.8
1.25
1.11032152
8
0.97118675
9
1.2
0.83333333
0.87055056
1.4
0.71428571
0.79364435
11
1.6
0.625
0.73253901
12
0.5
0.6407443
13
0.33333333
0.50237729
3.2材料参数
该重力坝坝体采用两种混凝土材料,具体材料参数见表3.2。
其中混凝土动力分析弹性模量是静力分析弹性模量的1.5倍。
计算时考虑为完全沉降后作用,故不考虑基岩密度。
表3.2坝体混凝土材料属性
参数
混凝土
材料编号
密度
(kg/m3)
弹性模量
(Pa)
泊松比
坝体混凝土Ⅰ
2400
2.6E10
0.167
坝体混凝土Ⅱ
2.85E10
基岩
2.9E10
3.3规范要求
如未特殊说明,本次分析参照的规范均指《混凝土重力坝设计规范》
NB/T35026-2014,以下为本次分析需用到的部分重要章节。
1)正常使用极限状态
按材料力学方法进行坝体上、下游面混凝土拉应力验算;
必要时进行坝体及结构变形计算、复杂地基局部渗透稳定验算。
正常使用极限状态作用效应采用下列设计表达式:
对正常使用极限状态验算时,作用分项系数、材料性能分项系数都取1.0,结构重要性系数不变。
该规范提出了重力坝对正常使用极限状态的要求,规定坝踵及坝体上游面不产生垂直拉应力,施工期坝趾处垂直正应力可容许有不大于0.1MPa的拉应力,下游坝面主拉应力不大于0.2MPa
2)承载能力极限状态
承载能力极限状态,对坝体结构及坝基岩体进行强度和抗滑稳定计算,必要时进行抗浮、抗倾验算:
抗震设防应满足DL5073的有关规定。
作用分项系数:
结构重要性系数:
设计状况系数:
0.85(地震)
结构系数:
1.5(地震)
材料性能分项系数:
3)抗滑稳定验算
坝体混凝土与基岩接触面的抗滑稳定极限状态。
作用效应函数为:
抗滑稳定抗力函数为:
式中:
----坝基面上全部切向应力作用之和(KN);
-------坝基面抗剪断摩擦系数;
------坝基面抗剪断凝聚力(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基岩有限元模型图4.4坝体有限元模型
图4.5上部坝体有限元模型图4.6下部坝体有限元模型
表4.1有限元模型具体数据
重力坝
单元数
节点数
单元类型
模型整体
18560
21714
1、2、3
Solid45、Solid65
坝基
16800
19734
Solid45
坝体
1760
2346
1、2
Solid65
坝体上部
880
1224
坝体下部
4.2边界条件
基岩底部采用全约束,上下游坝基采用法向约束,即将X方向的两个侧面的X向位移进行约束,坝基两侧采用法向约束,即将Z方向的两个侧面的Z向位移进行约束,上部边界为自由边界,假设坝体与坝基固结。
4.3荷载工况
(1)自重:
混凝土材料参数见3.2节;
(2)静水压力:
水容重取9.81KN/m3;
(3)扬压力:
扬压力计算简化坝踵处为无排水孔。
(4)地震作用:
反应谱见3.1节。
根据荷载的不同分为六种工况,工况一~工况三为静力分析,工况四~工况六为动力分析。
表4.2分析工况一览表
荷载
工况
自重
静水压力
扬压力
地震荷载
备注
工况一
√
库空
工况二
正常使用极限状态
工况三
0.95
承载能力极限状态
工况四
模态分析
工况五
谱分析
工况六
偶然荷载
注:
表格中“√”表示工况内含有该荷载,表中系数表示该荷载的作用分项系数;
静水压力、扬压力按坝面分布荷载施加。
5计算成果
5.1工况一
1)计算结果
图5.1整体模型总位移(m)图图5.2坝体UX(m)图
图5.3坝体UY(m)图图5.4坝体S1(Pa)图
2)本节小结
表5.1位移结果汇总表
结果
位置
X方向位移
(cm)
Y方向位移
总位移
坝顶
-1.54
-2.02
2.54
坝踵
0.04
-0.86
坝趾
-0.28
表中数据均指最大值,其中总位移为整体模型的总位移最大值。
X方向为负值表明坝体倾向上游,Y方向位移为负值时表明坝体沉降,下同。
表5.2坝体应力结果汇总表
位置
应力
最大拉应力
最大压应力
S1(MPa)
-0.67
0.20
3.27
表中数据指该位置处的最大值,下同。
库空工况下,坝体重心偏向上游,坝体整体向上游倾斜,由此产生的位移矢量总体趋势是倾向上游,见图5.1。
坝体水平变位从建基面至坝顶,随高程增加而加大,坝顶水平位移最大,其值为1.54cm(向上游),见图5.2。
对同一高程而言,坝体上游侧垂直变位大于下游侧垂直向变位,坝顶铅直向变位最大,其值为2.02cm(沉降),坝踵铅直向变位0.86cm(沉降),坝趾铅直向变位0.28cm(沉降)。
对于坝体而言,其最大主应力为拉应力,坝体最大主应力(图5.4)随越靠近两侧坝肩越大,同时,在坝肩附近有应力集中现象,最大拉应力为3.27MPa。
5.2工况二
图5.5整体模型总位移(m)图图5.6坝体UX(m)图
图5.7坝体UY(m)图图5.8坝体S1(Pa)图
表5.3位移结果汇总表
1.13
-1.06
1.53
0.36
0.25
-0.47
表5.4坝体应力结果汇总表
-0.23
0.17
2.57
-1.03
正常使用极限状态下,从图5.5中可以看出,坝体运行期间由于受库水压力+扬压力作用,产生的位移矢量总体趋势是偏向下游,坝体水平向变位倾向下游,垂直沉降变位量值在同一高程下游坝面大于上游坝面(图5.7)。
正常蓄水位工况坝顶水平向位移为1.13cm(向下游),铅直向变位值为1.06cm(沉降)。
由于坝基扬压力作用,沉降量值较库空工况明显减小。
坝体最大主应力为拉应力,就同一高程而言,上游侧最大主压应力大于下游侧(图5.8)。
在坝肩靠近上游侧附近有应力集中现象,最大拉应力为2.57MPa。
5.3工况三
图5.9整体模型总位移(m)图图5.10坝体UX(m)图
图5.11坝体UY(m)图图5.12坝体S1(Pa)图
表5.5位移结果汇总表
1.19
-0.93
1.50
0.37
-0.13
-0.39
表5.6坝体应力结果汇总表
0.16
-1.02
承载能力极限状态下,坝体整体变位规律和应力分布规律与工况二的正常使用极限状态一致。
但该状态重力、静水压力、扬压力分别乘以相应的作用分项系数后,坝体铅直向沉降量值较正常使用极限状态明显减小。
5.4工况四
图5.13坝体第一阶振型图5.14坝体第二阶振型
图5.15坝体第三阶振型图5.16坝体第四阶振型
图5.17坝体第五阶振型图5.18坝体第六阶振型
表5.7模态分析前10阶自振频率
*****INDEXOFDATASETSONRESULTSFILE*****
SET
TIME/FREQ
LOADSTEP
SUBSTEP
CUMULATIVE
3.3644
5.5838
5.9451
6.1252
6.4190
8.2269
8.3314
8.6114
9.7399
11.1980
5.5工况五
图5.19整体模型总位移(m)图图5.20坝体UX(m)图
图5.21坝体UY(m)图图5.22坝体S1(Pa)图
表5.8位移结果汇总表
2.75
0.98
2.92
0.14
0.05
表5.9坝体应力结果汇总表
2.29
1.26
4.86
从图5.19中可以看出,坝体位移矢量总体趋势是偏向下,坝体水平向变位倾向下游,坝体水平变位从建基面至坝顶,随高程增加而加大,坝顶水平位移最大,其值为2.75cm(向上游),见图5.20。
对同一高程而言,坝体上游侧垂直变位大于下游侧垂直向变位,坝顶铅直向变位最大(图5.21)。
对于坝体而言,最大主应力为拉应力,在坝肩靠近下游侧附近有应力集中现象最大拉应力为4.86MPa。
5.6工况六
图5.23整体模型总位移(m)图图5.24坝体UX(m)图
图5.25坝体UY(m)图图5.26坝体S1(Pa)图
表5.10位移结果汇总表
3.90
-0.50
3.93
-0.05
-0.34
表5.11坝体应力结果汇总表
1.01
5.72
-0.33
自重+静水压力+扬压力+地震荷载条件下,坝体位移矢量总体趋势是倾向下游,见图5.23。
坝体水平变位从建基面至坝顶,随高程增加而加大,坝顶水平位移最大,其值为3.90cm(向上游),见图5.24。
对同一高程而言,坝体上游侧垂直变位小于下游侧垂直向变位,坝体下游坡面中间位置铅直向变位最大。
对于坝体而言,其最大主应力为拉应力,在坝肩附近有应力集中现象,最大拉应力为5.72MPa。
5.7结果总结及分析
表5.12各工况位移结果一览表
工况
表5.13坝体第一主应力S1(MPa)结果一览表
1)坝体变位分析
坝体变位规律符合常规,库空条件下,坝体重心偏向上游,坝体整体向上游倾斜,由此产生的位移矢量总体趋势是倾向上游;
正常使用极限状态下,从图5.5中可以看出,坝体运行期间由于受库水压力+扬压力作用,产生的位移矢量总体趋势是偏向下游,坝体水平向变位倾向下游,垂直沉降变位量值在同一高程下游坝面大于上游坝面。
2)正常使用极限状态强度分析
由图5.5和表5.13可知,坝体坝踵和坝趾处未产生拉应力,符合规范要求,但坝体与两岸连接处坝肩的上游侧附近第一主拉应力达到2.57MPa>
1.78MPa,已超过混凝土的强度标准值,发生破坏,应采取加固措施。
3)承载能力极限状态强度分析
坝体材料C20,轴心抗压强度
,轴心抗拉强度
。
①静力工况
根据3.3节的公式和系数可得
,
即坝体压应力满足强度要求,拉应力不符合规范要求。
②地震工况
4)抗滑稳定验算
使用ANSYS映射面功能(具体操作见命令流)算得:
(方向向下),根据3.3节的系数求出
所以该重力坝满足稳定性要求。
6结论及建议
本文使用有限元分析软件ANSYS对坝体进行了静力分析和动力分析,通过以上的计算和分析,可以得出以下结论。
(1)重力坝变位规律符合常规,重力坝水平方向的最大位移发生在坝顶,地震荷载作用下,坝顶水平位移更大。
对此,为了提高该重力坝的抗震性,应采取措施减小坝顶水平变位,如可以适当增大顶部刚度,减轻重量,提高混凝土标号,必要时可以在上下游面布设钢筋。
(2)正常使用极限状态和承载能力极限状态下,在静力工况下,坝体的坝踵和坝趾位置处的应力值满足规范要求;
在地震工况下,坝踵处出现1.01MPa的拉应力,可在此处适当增加配筋,防止地震之后裂缝发生和扩展;
确保混凝土的施工质量,加强温控措施和养护,尽量减少表面裂缝。
并且坝体的两岸坝肩处的拉应力较大,远超出混凝土的强度标准值,应采取加固措施。
(3)通过稳定验算可知,该重力坝能够满足抗滑稳定性要求。
为进一步提高该重力坝的抗滑稳定性,建议在靠近坝体上游面,采用深孔锚固高强度钢索,并施加预应力,既可增加坝体的抗滑稳定,又可消除坝踵处的拉应力。
7分析命令流
本次分析以APDL语言为主,GUI操作为辅,主要的操作均编写成命令流形式(见下),但其运行时需注意其个别参数的调整(如作用分项系数),荷载工况组合主要以GUI操作为主,故此处未附该操作的命令流。
FINISH$/CLEAR
/FILNAME,Gravitydam$/PREP7
!
一.静力分析
ET,1,SOLID65$ET,2,SOLID45$ET,3,MESH200$KEYOPT,3,1,6
静力分析材料属性
MP,DENS,1,2400$MP,EX,1,2.6E10$MP,PRXY,1,0.167
MP,DENS,2,2400$MP,EX,2,2.85E10$MP,PRXY,2,0.167
MP,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,-360
K,8,0,-360$K,9,-270,$K,10,139.5,-360$K,11,409.5,-360$K,12,409.5$K,13,409.5,180
K,14,-270,180
A,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,10
LSEL,S,,,2,4,2$CM,LN4,LINE$ALLSEL$LESIZE,LN4,,,4
LSEL,S,,,1,5,2$LSEL,A,,,18$CM,LN5,LINE$ALLSEL$LESIZE,LN5,,,5
LSEL,S,,,6,7$CM,LN18,LINE$ALLSEL$LESIZE,LN18,,,18
LSEL,S,,,9,12,3$CM,LN22,LINE$ALLSEL$LESIZE,LN22,,,22
LESIZE,8,,,10,2$LESIZE,13,,,10,0.5$LESIZE,10,,,10,0.5$LESIZE,11,,,10,2
LESIZE,14,,,10,2$LESIZE,16,,,10,0.5$LESIZE,17,,,10,2$LESIZE,19,,,10,2
LESIZE,15,,,10,0.5$LESIZE,20,,,10,0.5
LCCAT,2,6$LCCAT,4,7
AATT,,,3$MSHKEY,1$AMESH,ALL$LDELE,21,22$EPLOT
TYPE,1$EXTOPT,ACLEAR,1$EXTOPT,ESIZE,8$EXTOPT,ATTR,0
MAT,1$ESYS,0$VEXT,1,2,,,,180*0.75
TYPE,2$EXTOPT,ACLEAR,1$EXTOPT,ESIZE,8$EXTOPT,ATTR,0
MAT,3$ESYS,0$VEXT,5,7,1,,,180*0.75
MAT,3$ESYS,0$VEXT,1,7,1,,,-180
WPOFFS,,,180*0.75$CSYS,4$VSYMM,Z,ALL$WPCSYS,-1$CSYS,0
NUMMRG,ALL$NUMCMP,ALL
更改材料属性
ESEL,S,TYPE,,1$NSLE,R,ALL$NSEL,R,LOC,Y,0,100$ESLN,S,ALL
EMODIF,ALL,MAT,2$ALLSEL$EPLOT
WPCSYS,-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$ALLSEL
NSEL,S,LOC,Z,450!
坝肩两侧Z方向UZ法向约束
NSEL,A,LOC,Z,-180$CM,Z3,NODE$ALLSEL
D,Y1,ALL$D,X2,UX$D,Z3,UZ
施加荷载
ACEL,,9.81!
/PSF,PRES,NORM,2
SFGRAD,PRES,,Y,170,-9810!
施加静水压力
ESEL,S,MAT,,1,2$NSLE,S,ALL$NSEL,R,LOC,X,0$NSEL,R,LOC,Y,0,170
CM,SHUI,NODE$SF,SHUI,PRES,0$ALLSEL
SFGRAD,PRES,,X,139.5,-1.1*9.81*1000*170/139.5!
施加扬压力
ESEL,S,MAT,,1,2$NSLE,S,ALL$NS
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 ANSYS 重力坝 三维 动态 结构 分析