教育资源为主的文档平台

当前位置: 查字典文档网> 所有文档分类> 工程科技> 兵器/核科学> 向家坝水电站地下厂房围岩稳定的黏弹塑性有限元分析

向家坝水电站地下厂房围岩稳定的黏弹塑性有限元分析

上传者:沈虹
|
上传时间:2015-04-26
|
次下载

向家坝水电站地下厂房围岩稳定的黏弹塑性有限元分析

第3l卷第4期2010年4月

岩土力学

V01.31No.4

RockandSoilMechanics

Apr.2010

文章编号l

1000--7598(20LO)04一1276—07

向家坝水电站地下厂房围岩稳定的

黏弹塑性有限元分析

聂卫平,徐卫亚,周先齐,李良权,刘

(河海大学岩土工程研究所。南京210098)

摘要:考虑流变作用,建立三维有限元模型进行数值计算,以洞室变形和点抗滑安全系数为指标,针对向家坝水电站地下厂房雨岩的特殊性进行稳定性研究,结果表明,随着围压的增高,流变速率逐渐减小,初始应变逐渐减小;软弱夹层处流变速率较其他岩体减小缓慢,且开挖后流变达到稳态状态时软弱夹层最终流变位移较大;黏弹塑性下围岩位移分布及变化规律与弹塑性一致,但黏弹塑性下计算位移明显要比弹塑性大;流变效应对岩体变形和稳定性,以及对支护结构有重要影响;黏弹塑性情况下,洞室围岩特征点抗滑安全系数比弹塑性条件下小,软弱夹层出露处和拱顶点抗滑安全系数较低,点抗滑安全系数分析还表明,软弱夹层对其稳定性影响明显,验证了位移分析结果。关键词:弹颦性;黏弹塑性;流变;洞窜变形;点抗滑安全系数;流变位移中图分类号I

TV554;O242.21

文献标识码l

Viscoelasto—plasticfiniteelementanalysisofstabilityofsurroundingrock

forundergroundpowerhouseofXiangjiabahydro—powerstation

NIE

Wei ping,XUWei—ya,ZHOUXian—qi,LI

Liang—quan,LIU

Wei

(GeotechnicalResearchInstitute,HohaiUniversity,Nanjing210098,China)

elementmodelfornumericalcalculationisestablishedconsideringrheology.Accordingto

Abstract:Athree-dimensionalfinite

particularityofthesurroundingrockforstudiedbasedOIltheindexesofandinitialstraindecrease

cavern

underground

cavernof

Xiangjiabahydro-powerstation,stabilityofsurroundingrockis

deformationandpointsafetyfactoragainstsliding.Itisshownthattherheologicalspeedrate

graduallywitllconfiningpressureincreasing.Comparedwitllotherrocks.therheologicalspeedofsoft

interlayerislargerthanotherstateswhenrheologyreaches

rockof

interlayerdecreasesslowly.Finallyrheoiogicaldisplacementofsoftlowerstate.Thedisplacementdistribution

and

change

lawof

surrounding

viscoelastic—plasticityandelastoplasticityare

on

uniform;but

supportsliding

cavern

thecalculationdisplacement

isdone

ofviscoelastc。plasticityislarger.Animportantinfluence

rock

deformation,stabilityand

structure

byrheologicaleffect.Atthestateofviscoelastic-plasticity,characteristicpoints’safetyfactorsagainst

arelessthancharacteristicpoints’safetyfactorsofelastoplasticity,pointsafetyfactoragainstslidingofsoftinterlayerand

arch

crown.Itis

cavern

also

shown

thatthesoft

interlayerhasanimportantinfluence

on

stabilityof

surroundingrockfor

undergroundbystudyingpointsafetyfactoragainstsliding;accordingly,thedisplacementanalysisresultsaleverified.

Keywords:elastoplasticity;viscoelasto-plasticity;theology;deformationofcavern;pointdisplacement

safetyfactoragainstsliding;rheoiogical

引言

群围岩稳定性开展的工作大多是对一些特定工程进行的研究,如Chryssanthakis等1111用有限元、DDA或二维和三维离散元方法对特定工程进行了稳定性分析:Sitharam等【4】对节理岩体中洞室围岩提出了较新的力学类型;Hibino和Motojmal5J用现场实测方法研究了日本十余个大型水电地下厂房的围岩变

我国是个多山国家,电站厂房群采用地下开挖的方法是一个发展趋势,因此,研究不同条件下地下洞室围岩稳定性就成为一个重要的研究课题。近年来,国内外学者或工程界人士发表的对地下洞室

收稿日期:2009-03.1I

基金项目:国家自然科学基金资助项1|(505391lO);国家十一五科技支撑项目(2006BAB04A06)。

第作者简介:聂卫平,男,1982年生,博士研究生,主要从事岩石力学与工程等方面的研究工作。E-mail:nwp_jzit@tom.com

万方数据 

第4期聂卫平等:向家坝水电站地下厂房围岩稳定的黏弹蠼性有限元分析

1277

形规律;安红刚16J综合应用遗传算法、人工智能、神经网络理论以及并行计算,提出了大型洞室群优化的并行进化神经网络有限元方法;朱维申等【7峭l对典型的洞室群型在不同条件和多因素影响下高边墙位移规律的预测进行了研究。而随着岩体工程规模的扩大及深部岩石工程的发展,某些地下洞室围岩变形的时间效应越来越突出,甚至成为影响洞室稳定性的关键因素。这种与时间相关的流变特性,单纯地运用弹塑性理论是无法模拟的,需要用新的

模型和方法对其进行研究。如Hsieh等p懈蠕变变

形(黏性变形)分成体积蠕变和剪切蠕变两部分,并分别采用Taylor体积蠕变规律和Singh.Mitchell偏态蠕变规律来求取;张玉军、刘谊平等IloJ建立了一种正交各向异性岩体的黏弹一黏塑性模型对层状岩体(横观各向同性介质)中的地下洞室围岩的流变动态进行了分析;邵国建、卓家寿等【lIj考虑了洞室开挖过程中的黏性效应,从理论和经验两方面讨论了迭代过程中时间步长的限制问题,对某洞室围岩稳定性所进行的分析;曹树刚、边金等¨2J通过对岩石的全应力.应变曲线和蠕变曲线的分析,提出了一种组合模型,能较好地反映岩石的非衰减蠕变特性。董志宏、丁秀丽等【l3】运用均匀设计.神经网络.遗传算法方法对构皮滩水电站地下厂房区的软岩流变参数进行了反演分析。虽然洞室流变研究较多,但向家坝地下洞室有其特殊性,本文基于前人研究的基础上,应用试验分析的结果,采用美国ANSYS公司的ANSYSl0.0软件建立三维有限元模型,并采用侧压力系数法对地下洞室区地应力进行施加,针对向家坝水电站地下洞室的特殊性开展流变分析,以研究流变作用下地下洞室围岩的稳定性,以满足工程施工、设计等的需要。文中选用洞室变形及点抗滑安全系数为考察指标,能很好地反映洞室围岩的稳定性变化规律。

2工程概况

向家坝水电站地下输水发电系统布置在坝址右岸,为一组大型地下洞室群,洞群结构复杂。主厂房高度为85.5m,跨度为31m,岩锚吊车梁以上最大跨度为33m;主变洞高度为23.74m,跨度为

26.3

m。洞室规模大,其中主厂房开挖跨度在国内

外已建、在建和拟建的水电站中居第一,其跨度和规模均居世界前列。主厂房第1层顶拱岩性主要为泥质粉砂岩、粉砂质泥岩透镜体及粗砂岩;洞室项拱以II类围岩为主,有部分围岩为IⅡ~Ⅳ类,厂房洞内出露的主要软弱夹层有JC2.2、JC2—3和JC2-4。

万 

方数据3

黏弹塑性流变本构模型及计算参数的确定

分析砂岩、粉砂质泥岩和软弱夹层试验曲线,

瞬时变形性质上砂岩主要表现为弹性,粉砂质泥岩主要表现为弹塑性,软弱夹层主要表现为塑弹性;流变变形性质上,软弱夹层的黏性最突出,粉砂质泥岩次之,砂岩的黏性相对最小,在高应力水平下,三种岩石均具有黏弹性和黏塑性性质。因此,采用Burgers黏弹性流变模型和Mohr.Coulomb塑性屈服准则耦合的Burgers黏弹塑性流变本构模型进行向家坝地下厂房洞室群岩体的流变计算,模型示意图如图l。将地层按岩性归类,根据流变试验及流变参数辨识结果,确定地下厂房区各岩组Burgers黏弹翅性流变模型参数,如表1所示。弹塑性数值计算采用有限差分法,强度准则采用Mohr-Coulomb准则,各岩组力学参数建议值见表2。

☆^开尔文模型

马克斯韦尔模型

.广卜、

摩尔库伦模犁

圈1伯格斯黏弹塑性流变模型示意图

Fig.1

SketchofrheologicalmodelforBurgers

viscoelasto-plasticity

模型中的参数有:弹性体积模量膨开尔文剪

切模量GK;开尔文黏滞系数矿;马克斯韦尔剪切

模量G匕马克斯韦尔黏滞系数叮肼;黏聚力气;摩

擦系数£。

表1各岩组伯格斯黏弹塑性流变模型参数

Table1

Allpetrofabric

parametersofrheological

巴竺璺!!!竺!曼竺!g!!!!!!!竺!!苎!!竺:P!苎!望!!壁

T3’9.886.5017505231801.46lT3“lO.977.55

1750523

180

1.55

lT3”11.1l7.65

1750523180

1.55lT3“10.487.211750523180

1.52l

T3“

2.451.47

1604910

0.73

T3J

4.943.25

1604910

1.04OT3“3.2l2.02160

49100.87OT3“5.133.23160

49100.92O夹层

3.93

2.36

150

98

10

0.1I

注;c.、L为流变稳态时岩体黏聚力和摩擦系数

1278

岩土表2弹塑性数值计算岩体力学参数

Table2

MechanicaI

parametersofrock

massfor

elastoplastic

numericalcalculation

地层弹模

泊松比

重度

黏聚力c摩擦系数

参数

/GPaiz/(kN/m3)

/MPa

/l~4_PaT34

16.00

0.23

25.52

1.57

1.221.621"3’8.000.2325.34

1.23

1.05

1.23T32”18.800.22

25.491.701.281.76T1‘~18.01

0.2225.50

1.631.241.69Tl‘一

18.240.2225.50

1.651.251.7l

B‘+118.660.2225.491.671.26

1.73

T3。’3.670.2525.630.870.87

0.84T3“18.660.2225.491.671.26

1.73T3“5.000.2425.641.030.96

1.00I"3“17.600.2225.501.631.24

1.68夹层

2.00

O.25

20.00

0.10

0.25

0.00

注:C、厂为瞬时弹塑性岩体粘聚力和摩擦系数。

4模型的建立

三维模型研究范围为:X方向(主厂房横向,指向上游为正)0~350m;Y方向(主厂房纵向,指向山体内为正)0~550m;Z方向(竖直方向,向上为正)0.~550m,最大高度为550111。模型共划分单元172142个,节点79596个,采用实体单元模拟位移边界条件。地下厂房整体示意图如图2。根据设计院提供的资料,采用侧压力系数法对地下洞室区地应力进行施加,整个洞室开挖共分1l步,具体分步开挖步骤见图3和表3。

图2地下厂房洞室示意图

Fig.2

Diagrammaticsketchofunderground

powerhousecavern

图3厂房区三维仿真模型分层示意图

Fig.3

Laminationsketchofpowerhouseregion3-dfinite

elementmodel

万 

方数据力学2010笠

表3厂房区洞室开挖分步

Table3

Excavationstepsofpowerhouse

cavern

开挖步骤

洞事部位

第1步1号、3号施丁支洞和5号施r支洞、交通洞、通风洞

第2步

母线竖井

第3步主厂房、主变洞、尾水洞2段、尾水洞3段、引水洞第1层第4步主厂房、主变洞、尾水洞2段、尾水洞3段、引水洞第2层第5步主厂房、主变洞、尾水洞2段、尾水洞3段第3层、

第6步

主厂房、尾水洞3段第4层、

第7步

主厂房第五层第8步主厂房第六层第9步主厂房第七层第10步

尾水洞1段第1层第11步

尾水洞1段第2层

结果分析

5.1特征点位置

为分析洞室应力位移变化趋势,选取了l撑、2撑、

3群、斜4个典型断面进行计算分析。l撑断面和2撑断面为主厂房断面,位置分别为Y=278.44

m、

Y=162.43

m;3撑断面和4j|}断面为主变洞断面,位

置分别为Y=278.44m、y=162.43m。每个断面都布置一些特征点进行位移及点抗滑安全系数信息的记录,断面及特征点位置见图4。5.2特征点流变位移分析

由数值计算呵得各特征点流变位移曲线,选列主厂房l{|}断面流变合位移图,见图5。

由流变合位移图可以发现:流变位移随时间的增长逐渐增加,但增大的速度逐渐减小,即随着围压的增高,流变速率逐渐减小,初始应变逐渐减小;软弱夹层流变速率减小的缓慢,且随各步开挖后流变达到稳定状态,软弱夹层流变位移较其他岩体流变位移大。

5.3各断面流变位移分析

根据上述计算模型、流变参数、支护参数,通过有限元数值分析,计算结果最终整理如图6---7所示。综合考虑洞室围岩轴向、切向与径向三个方向的变形,据此分析如下:

(1)弹塑性条件下,洞室围岩在开挖完毕后不同部位向临空面变形大小不一致,地下厂房区l撑断面拱顶最大变形为35.24mm,底板最大变形为

24.04

mm,上游侧边墙最大变形为51.07mm,位于JC2-2上游侧边墙出露处,下游侧边墙最大变形为

38.39

mm,位于JC2.2下游侧边墙出露处;2群断面拱项最大变形为31.87mm,底板最大变形为

26.28

mm,上游侧边墙最大变形为49.55mm,位于JC2.3上游侧边墙出露处,下游侧边墙最大变形

50.23

内容需要下载文档才能查看 内容需要下载文档才能查看

mm,位于JC2.3下游侧边墙出露处。主变洞

第4期聂卫平等:向家坝水电站地下厂房围岩稳定的黏弹塑性有限元分析

1279

区3稃断面拱顶最大变形为37.35mm,底板最大变形为22.15inm,上游侧边墙最大变形为22.84

mm,

位于上游侧边墙中上部,下游侧边墙最大变形

32.19

mm,位于下游侧边墙中上部:4群断面拱顶最

图4断面及特征点位置

Fig.4

G。aph0f3ectiousand。ha。acteristic

pointspositions

万 

方数据BIB3

量静趔制爝

0200400600

800

l000

i200

流变时间/d

(a)非软弱夹层特征点流变位移

出露P2出露P4

出露P6

姜静趟锹

∞如∞加∞如柏如加m

内容需要下载文档才能查看

200

400

600

80010001200

流变时间,d

(”软弱夹层特征点流变位移

图5各点流变合位移图

Fi昏5

Rheologicalresultantdisplacementgraphsofcharacteristicpoints

mm,底板最大变形25.56mm,上游mm,位于上游侧边墙中上mm,位于JC2 2下(2)黏弹塑性条件下:软弱夹层对洞室围岩变mm,位于2撑断面下游mm,位于(3)黏弹塑性洞周位移分布及变化规律与弹塑

大变形27.04侧边墙最大变形24.72部,下游侧边墙最大变形30.11游侧边墙出露处。围岩主要向临空面变形,顶板下沉,边墙向洞内变形,底板向上,形成“底鼓”现象。

形影响比较明显,主厂房边墙围岩在软弱夹层处变形比较大,最大达到68.24侧边墙JC2.2出露处。主厂房位移都表现为边墙最大,顶拱和底板次之;主变洞底板围岩在软弱夹层处变形比较大,支护前最大达到60.313#断面拱顶,3群断面由于不受软弱夹层影响,位移主要表现为顶拱最大,边墙次之,底板最小,4撑断面由于受到软弱夹层影响,位移表现为边墙最大,项拱次之,底板最小。由此可知,软弱夹层对洞室围岩变形的影响较大。

性一致,洞室围岩拱顶位移向下,底板位移向上,

边墙向洞内变形;对于主厂房区域1}≠断面上游边墙

1280

弹塑性位移

鼯蕊

6I.

c49.36

40/

Ⅶ叫薹彤

吼¨“妇土

溉鬈m

24‘3

l例婴鲫薹垦

12.

(a)l#断面

弹塑性位移哩乃沁。

踟汕

题甜9§泳

彤如m

"弘

煅68.24醢

内容需要下载文档才能查看

24.75

蚋刚

30.47

69

l:Z!i望;22:5,.39n5

(b)2#断面

一瓣挚%黧兰一缝辫蒜瓣

(c)3#断面

一404”.di薪严23弋遗

12

熬杂穗霾

(d)锥断面

图6弹塑性和黏弹塑性总位移对比

Fig.6

Resultantdisplacement

comparisonof

elastoplasticityandviscoelasto-plasticity

万 

方数据力

2010生

圈7弹塑性和黏弹塑性开挖后围岩位移矢量

(长度代表位移量值)

Fig.7

Rock

displacement

vectorsafterexcavatIon

completion

ofelastoplasticityandviseoelasto-plasticity

(1ength

for

displacements)

部位的位移高于下游边墙,2撑断面下游边墙部位的位移高于上游边墙,下游边墙和上游边墙位移中部较大,两侧偏小,这主要是由于边墙中部有软弱夹层穿过所引起;对于主变洞区域3}|}和4撑断面,由于距离地下厂房较近,上游侧边墙和底部易受厂房的挠动影响,上游边墙部位的位移高于下游边墙;黏弹塑性计算所得位移明显要比弹塑性计算所得位移大,对于主厂房位移,黏弹塑性位移顶拱平均增幅为55.56%,底板平均增幅为37.83%,上游侧边墙平均增幅为34.62%,下游侧边墙平均增幅为32.05%;对于主变洞位移,黏弹塑性位移项拱平均增幅为69.11%,底板平均增幅为29.75%,上游侧边墙平均增幅为45.63%,下游侧边墙平均增幅为78.11%;受结构面切割影响,洞室拐角和分岔处都有比较大变形。

5.4特征点抗滑安全系数分析

通过有限元数值分析,可以得到断面特征点抗滑安全系数对比,对比结果见图8。据此分析如下:

(1)特征点抗滑安全系数均大于l;黏弹塑性条件下,各断面特征点抗滑安全系数较弹塑性条件

下小,说明考虑时间的影响,由于产生滞后位移,使得洞室围岩稳定性降低。

(2)弹塑性条件下:l群断面点抗滑安全系数最小为1.33,位于JC2.2下游侧边墙出露Pl点:2稃断面点抗滑安全系数最小为1.0,位于JC2.2两侧边墙出露P1、P2点。主变洞区3撑断面边墙无软弱夹层出露,拱顶点抗滑安全系数较边墙小,拱顶G2点最小1.22:4撑断面点抗滑安全系数最小为1.O,位于JC2—2上游侧边墙出露P2点。黏弹塑性条件下:1撑断面点抗滑安全系数最小1.04,位于JC2.4底板出露P5点;2撑断面点抗滑安全系数最小为1.0,位

版权声明:此文档由查字典文档网用户提供,如用于商业用途请与作者联系,查字典文档网保持最终解释权!

下载文档

热门试卷

2016年四川省内江市中考化学试卷
广西钦州市高新区2017届高三11月月考政治试卷
浙江省湖州市2016-2017学年高一上学期期中考试政治试卷
浙江省湖州市2016-2017学年高二上学期期中考试政治试卷
辽宁省铁岭市协作体2017届高三上学期第三次联考政治试卷
广西钦州市钦州港区2016-2017学年高二11月月考政治试卷
广西钦州市钦州港区2017届高三11月月考政治试卷
广西钦州市钦州港区2016-2017学年高一11月月考政治试卷
广西钦州市高新区2016-2017学年高二11月月考政治试卷
广西钦州市高新区2016-2017学年高一11月月考政治试卷
山东省滨州市三校2017届第一学期阶段测试初三英语试题
四川省成都七中2017届高三一诊模拟考试文科综合试卷
2017届普通高等学校招生全国统一考试模拟试题(附答案)
重庆市永川中学高2017级上期12月月考语文试题
江西宜春三中2017届高三第一学期第二次月考文科综合试题
内蒙古赤峰二中2017届高三上学期第三次月考英语试题
2017年六年级(上)数学期末考试卷
2017人教版小学英语三年级上期末笔试题
江苏省常州西藏民族中学2016-2017学年九年级思想品德第一学期第二次阶段测试试卷
重庆市九龙坡区七校2016-2017学年上期八年级素质测查(二)语文学科试题卷
江苏省无锡市钱桥中学2016年12月八年级语文阶段性测试卷
江苏省无锡市钱桥中学2016-2017学年七年级英语12月阶段检测试卷
山东省邹城市第八中学2016-2017学年八年级12月物理第4章试题(无答案)
【人教版】河北省2015-2016学年度九年级上期末语文试题卷(附答案)
四川省简阳市阳安中学2016年12月高二月考英语试卷
四川省成都龙泉中学高三上学期2016年12月月考试题文科综合能力测试
安徽省滁州中学2016—2017学年度第一学期12月月考​高三英语试卷
山东省武城县第二中学2016.12高一年级上学期第二次月考历史试题(必修一第四、五单元)
福建省四地六校联考2016-2017学年上学期第三次月考高三化学试卷
甘肃省武威第二十三中学2016—2017学年度八年级第一学期12月月考生物试卷

网友关注

网友关注视频

沪教版八年级下册数学练习册21.3(3)分式方程P17
冀教版小学英语五年级下册lesson2教学视频(2)
七年级英语下册 上海牛津版 Unit5
【部编】人教版语文七年级下册《逢入京使》优质课教学视频+PPT课件+教案,安徽省
北师大版数学四年级下册第三单元第四节街心广场
3月2日小学二年级数学下册(数一数)
二年级下册数学第二课
【部编】人教版语文七年级下册《老山界》优质课教学视频+PPT课件+教案,安徽省
外研版英语七年级下册module3 unit2第二课时
沪教版八年级下次数学练习册21.4(2)无理方程P19
冀教版小学英语四年级下册Lesson2授课视频
沪教版八年级下册数学练习册20.4(2)一次函数的应用2P8
沪教版牛津小学英语(深圳用) 四年级下册 Unit 3
七年级英语下册 上海牛津版 Unit9
每天日常投篮练习第一天森哥打卡上脚 Nike PG 2 如何调整运球跳投手感?
苏科版八年级数学下册7.2《统计图的选用》
北师大版数学四年级下册3.4包装
【部编】人教版语文七年级下册《过松源晨炊漆公店(其五)》优质课教学视频+PPT课件+教案,江苏省
二次函数求实际问题中的最值_第一课时(特等奖)(冀教版九年级下册)_T144339
二年级下册数学第一课
冀教版小学数学二年级下册第二周第2课时《我们的测量》宝丰街小学庞志荣
沪教版牛津小学英语(深圳用) 五年级下册 Unit 12
外研版八年级英语下学期 Module3
冀教版小学数学二年级下册第二单元《余数和除数的关系》
【部编】人教版语文七年级下册《过松源晨炊漆公店(其五)》优质课教学视频+PPT课件+教案,辽宁省
六年级英语下册上海牛津版教材讲解 U1单词
沪教版八年级下册数学练习册一次函数复习题B组(P11)
沪教版牛津小学英语(深圳用) 四年级下册 Unit 2
【部编】人教版语文七年级下册《逢入京使》优质课教学视频+PPT课件+教案,辽宁省
苏科版数学八年级下册9.2《中心对称和中心对称图形》