溷凝土面板堆石坝的沉降分析.docx
- 文档编号:9779116
- 上传时间:2023-02-06
- 格式:DOCX
- 页数:17
- 大小:1.20MB
溷凝土面板堆石坝的沉降分析.docx
《溷凝土面板堆石坝的沉降分析.docx》由会员分享,可在线阅读,更多相关《溷凝土面板堆石坝的沉降分析.docx(17页珍藏版)》请在冰豆网上搜索。
溷凝土面板堆石坝的沉降分析
混凝土面板堆石坝的沉降分析
摘要:
土耳其的Kurtun坝是一座混凝土面板堆石坝,位于东黑海地区,人们对她的沉降过程做了调查和研究,在土耳其,这是第一座被监测的混凝土面板堆石坝。
人们利用计算机程序PLAXIS,通过二维平面应变有限元分析原理,对大坝施工和水库蓄水阶段进行了坝体应力和位移的监测,并将观测的数据与预测的应力和位移做了比较,比较的结果显示,在大坝的施工阶段,二者整体上十分吻合。
由于岸坡较陡,横跨山谷的拱对坝体的应力和位移有重要影响。
可以看出,在水库蓄水阶段,预测的位移比实际观测的要大,这就表明了,水库蓄水阶段相对于施工阶段,堆石坝的坝体强度表现的相对较高。
土的压缩特性曲线通常是非线弹性的,并且压力的大小和堆石坝的材料特性有关。
该模型的参数可以参考相近的工程实例从相关的规范中选取。
关键词:
混凝土面板坝位移有限元分析固结模型压力
1、引言
现在,混凝土面板堆石坝(CFRD)在世界范围内非常流行,特别是在多雨地区,和地基相对透水的地区。
在现行的技术阶段,混凝土面板堆石坝的设计主要是靠工程经验和实践(库克,1984).因为这些建筑物非常重要,所以他们在施工阶段和水库蓄水阶段的变化过程必须精确的估算出来,而有限元分析是预测建筑物变化过程的有效方法中的一个方法。
有限元分析的关键在于以三轴实验为基础,建立一个符合实际的堆石坝材料应力应变模型。
由于堆石坝材料的颗粒尺寸,直径可以达到1.2m,所以精确的实验数据很难得到,并且在工程实例中也很少。
虽然可用的实验数据很少,但是一些文献中的三轴实验数据表明堆石材料据有抗压强度高和应力应变曲线呈非线弹性的特征(Marsal,1967;Marachi等人,1972;Duncan等人,1980)。
通过近几年的研究,在众多适用的材料模型中,邓肯-张的双曲线函数模型最具有代表性,并且邓肯-张的双曲线函数模型在一些近似的实验中运用的非常成功(Khalid等人,1990;SaboyaandByrne,1993;Liu等人2002)。
在现有的研究中,常用“土的固结模型”来代替堆石材料的变化过程,并且土的固结模型实际上是对在PLAXIS程序中应用的双曲线函数模型的改(Schanz等人1999)。
2、大坝情况
2.1.一般情况
Kurtun坝修建在东黑海地区,Gümüshane省Torul镇西北27公里的Harsit河上,修建大坝的主要目的是用来发电,年发电量达198GWh。
大坝的坝体工程于1997年开始施工,1999年竣工。
在石料填充完毕后,施工过程中断了1.5年,直至堆石坝体大部分沉降完成以后,混凝土衬砌工程才开始施工。
水库蓄水开始于2002年2月8日,至2002年4月30日水库水面高程达到海拔630m。
Kurtun坝高出河床133m,上下游面的坡度分别为1.4:
1和1.5:
1(垂直投影长度与竖直投影长度相比),混凝土面板的坡度是1.406:
1,坝顶长300m。
大坝建在一个相对狭窄、岸坡较陡的峡谷里,山谷宽45m,左右两岸的平均坡度分别是61.32o和52.04o。
图一表示了坝的最大横剖面和不同分区。
图1Kurtun大坝最大横剖面
大坝坝址处的地质基本上是花岗闪长岩、辉绿细晶岩、安山岩以及石灰石,其中花岗闪长岩是大坝坝址处发现的最常见的一种岩石类型。
大坝建筑物和分区的详细资料见表1:
表一Kurtun坝的材料和建筑物的详细情况
2.2.大坝监测使用的仪器
因为kurtun坝是土耳其的第一座混凝土面板堆石坝,所以采用了很多仪器对大坝的变化进行监测。
用于监测大坝的四种类型的设备分别是:
-液压式沉降仪,位于堆石内(ZDO);
-液压式压力盒,在堆石坝体内(BO);
-应变片,在混凝土面板上(GO);
-表面式测缝计,位于混凝土面板内。
总共有33个液压沉降仪和22个土压力盒装在三个横断面上和四个不同的海拔高度,大坝最大横断面上的设备位置见图2,另外,在混凝土面板上还安装了16个测缝计和6个应变片。
图2堆石体的设备装置情况
3、分析和结果
3.1一般情况
正如前面的段落所阐述的,土的固结模型即邓肯-张双曲线函数模型的改进模型在研究中被广泛应用。
由于没有三轴实验的数据可以利用,土的固结模型参数一般从文献(SaboyaandByrne,1993)中选择。
表二例举了一些在初步分析中使用的一系列参数。
表二土的固结参数
在固结模型中,在首次加载的情况下,抗剪强度qf达到50%时对应的切割模量可以通过公式
求解。
式中E50ref是与侧限压力pref相对应的参考系数,E50ref取决于最小主应力σ3’即三轴实验中的侧限压力,若σ3’为负值则表示压力。
幂数m控制附加应力,类似于邓肯-张双曲线函数模型(1970)中的指数n。
在土的固结模型中,抗剪强度qf一般利用摩尔屈服理论估算。
图3双曲线函数应力应变曲线(负值表示压力)
2A分区和3A分区(见图1)不会对大坝坝体产生重大的影响,也没有用有限元进行分析。
在堆石坝的研究中,通常假定堆石材料间的凝聚力为0KPa,但是,为了提高计算精度,在这篇论文中,假定凝聚力为1Kpa。
把大坝的岩基视为刚体,不考虑大坝坝基的水平位移和竖直位移,同时对大坝施工阶段结束和水库蓄水完毕后进行分析。
3.2.有限元模型
利用基于有限元原理的计算机程序PLAXIS和二维平面应变原理对大坝进行分析,取15个三角形单元作为堆石坝坝体的代表并用土的固结模型作为材料模型。
这些单元有12个内部应力点,并且这些应力点的坐标和输出位移的单元节点不同。
混凝土面板取5个线弹性梁单元,这些单元有50cm厚和8个应力点。
混凝土的材料参数见表3。
初步研究时,对每一个施工阶段去一个10m厚的垫层作为研究对象,据观察,减小层厚可以增加模拟效果,同时延长运行时间,最终,垫层厚度选为5m。
由于防浪墙使坝的几何形状复杂化,所以坝顶最上面6m厚不在包含在网格内,尽管如此,在模型研究中,防浪墙的重量作为一个附加力仍被考虑在内,大小为100K
pa。
有限元网格如图4所示,网格包含了15417个计算节点和22608个应力点。
图4有限元网格
3.3.分析结果
3.3.1.沉降
3.3.1.1.施工完毕(EoC)。
完建情况下,大坝分析使用的一系列参数如表二所示,假定堆石坝体仅由3B区的材料组成。
若选取如下参数将会取得很好的效果:
表三混凝土面板的线弹性模型参数
当材料参数按表四中的数据选取时,监测数据和预测结果的比较见表五:
表四初步研究中使用的3B材料参数
表五施工完毕情况下的结果比较(仅用3B区材料)
从表五中的数据可知:
除了El555位置预测的位移比监测的数据要大,其他位置预测和监测的结果基本上相一致。
这种差异可能是由拱的效应使沉降增大引起的,由于本文采用的分析方法建立在二维平面应变假定的基础之上,所以并没有考虑拱的效应。
在以下的过程的中,3C区的堆石材料仍包含在模型中。
由于3C区材料的夯实层厚要比3B区材料的夯实厚度小,所以,开始时假定3C区材料的强度比3B区材料的强度低(HunterandFell,2003)。
分析过程中保持3B和3C区的幂数m(m=0.25)相同,并且3C区的E50ref值分别取19000KN/m2和23000KN/m2,以表示不同的强度。
观测表明:
当E50ref=19000KN/m2时,3C区F–F轴线上预测的沉降和观测的数值不是很符合。
同时发现当E50ref=23000KN/m2时,预测和观测的数据基本一致。
表六显示了三种不同条件下的分析结果:
表六施工完毕情况下3区采用不同材料参数时结果比较
3.3.1.2.水库蓄水完毕(RFC)。
首次蓄水对大坝来讲是个关键阶段,这是因为大部分后期沉降都发生在这个阶段,并且在这个阶段,水库的蓄水会引起堆石坝体的水平位移和竖直位移。
在蓄水期间,3B区的卸载变形和卸载参数建议取首次加载参数的合理倍数(NobariandDuncan,1972;Fitzpatrick等人,1985;SaboyaandByrne,1993)。
然而,由于使用的计算机程序PLAXIS7.1不要需要单独的卸载参数,所以分析时卸载过程并不考虑在内。
假定混凝土面板不透水、不开裂,水压力呈三角形分布,方向垂直指向混凝土面板。
分析的结果见表七:
表七蓄水完毕后预测和观测沉降比较
可以看出,沉降量最大的地方发生在上游面附近,并且向下游方向逐渐减小,和预期的一样。
同时,从表七中可以看出,预测的沉降量比观测的数值大一点,这是因为,首次蓄水时,库水被看作是引起堆石材料变形的附加荷载,所以电脑采用的是首次加载时的参数,事实上,卸载参数对沉降的影响远远大于库水对沉降的影响(SaboyaandByrne,1993)。
3.3.2.坝体压力
表八列举了施工完毕和蓄水完毕两种情况下坝体应力的比较结果,同时也列举了覆盖层的压力值:
表八施工完毕和蓄水完毕两种情况下坝体应力的比较
从表八中可以看出,在施工完毕的情况下,预测的数值和实际监测得到的数值,总的来说二者非常接近。
据悉,施工期间仪器BO-3出现了差错,所以监测的数据不见了。
同时位于混凝土面板下方的G-G轴线上的仪器BO-15也出现了偏差,可能是由于仪器的不当安装引起的。
在水库蓄水完毕情况下,除了与混凝土面板的BO-1轴线、BO-9轴线、BO-15轴线以及BO-20轴线平行的G-G轴线上的仪器所记录的数据之外,其他仪器实际监测的数值和预测的数值同样很接近。
这些偏差说明了一个事实,即这些仪器在混凝土面板施工之前也是在水库蓄水之前被重新安装过。
在水库蓄水期间,应力随着水库蓄水水位的升高而变大,靠近上游面的区域应力增加的快,而靠近下游面的区域应力增加则可以忽略不计。
同时,还可以看出,水库蓄水期间预测的数值和实际监测的数值之间的偏差比施工期间预测的数值和实际监测的数值之间偏差大一点儿,这也许可以通过水对堆石材料的软化效应来解释(NobariandDuncan,1972)。
3.3.3.应力和位移等值线
图5至图12表示了施工完毕(EoC)和水库蓄水完毕(RFC)两种情况下的应力分布。
水库蓄水对应力的影响可以从上游侧一半大坝的应力分布看出来。
水平方向的应力分布和竖直方向的应力分布很相似,特别是,水库的蓄水有使大坝向下游移动的趋势,并由此在上游一半坝体和下游一半坝体引起的正剪应力,这正如其他人所指出(Khalid等人,1990)。
图5施工完毕情况下的水平应力(单位Kpa,“-“表示压力)
图6施工完毕情况下的竖直应力(单位Kpa,“-”表示压力)
图7施工完毕情况下的剪应力(单位Kpa)
图8StressesinoutofplanedirectionforEoC
图9蓄水完毕情况下的水平应力(单位Kpa,“-“表示压力)
图10蓄水完毕情况下的竖直应力(单位Kpa,“-“表示压力)
图11水库蓄水完毕情况下的剪应力(单位Kpa,“-“表示压力)
图12StressesinoutofplanedirectionforRFC(单位Kpa,“-“表示压力)
图13至图16表示了竖直方向和水平方向的计算位移。
从图中可以看出来:
对于施工完毕的情况,大坝上游面的上部有向下游移动的趋势,而上游面的下部却有向上游移动的趋势,下游面的位移正好相反。
图13施工完毕情况下的水平位移(单位cm,指向下游为正)
图14施工完毕情况下的竖直位移(单位cm,“-”表示沉降)
图15蓄水完毕情况下的水平位移(单位cm,指向下游为正)
图16蓄水完毕情况下的竖直位移(单位cm,“-”表示沉降)
由于Kurtun坝水平方向的位移没有监测,所以没有机会将观测结果和计算的水平位移作比较。
对于上游一半坝体,计算的最大水平位移是19.5cm,发生在大坝轴线上游81.80m的El551.2位置处;对于下游一半坝体,计算的最大位移是17.3cm,发生在大坝轴线下游89.30m的El550.00位置处。
和竖直方向的位移比起来,水平方向的位移显得很小。
对于水库蓄水完毕的情况,水平方向的最大位移是36.09cm,发生在El560.63处。
水平方向的计算位移与先前的研究(Nobari和Duncan,1972)相比较,对于施工完毕和蓄水完毕两种情况下,计算的最大水平位移分别是199.76cm和52.79cm,分别发生在El580.00位置处和El565.63位置处。
4.结论
本文利用二维有限元原理对混凝土面板堆石坝Kurtun大坝做了分析,并将计算的位移和内部应力与相应位置的监测数据做了比较。
堆石材料特性用土的固结模型即改进的邓肯-张双曲线函数模型代替,材料的参数从相近的文献中选取合适的值,并对建筑物施工阶段和水库蓄水阶段进行了分析,发现利用土的固结模型和二维平面应变非线性有限元理论分析的结果和相应位置的实测数据非常吻合,计算机程序可以很好的反映堆石块体的变形特征。
由于u山谷比较狭窄,拱的效应对减小竖直位移有重要影响。
可以看出,当减小3C区的参数时,计算的沉降与实际观测的数值的偏差就会增大,这表明当3C区的材料处于受压状态时的变形特征和3B区材料的变性特征相似。
我相信对建筑物的侧向变形和应力的监测会对混凝土面板堆石坝的实际变形过程的研究会做出很大的贡献。
谢辞
衷心的感谢GeneralDirectorateofStateHydraulicWorks在提供Kurtun大坝监测数据方面的帮助。
参考文献
[1]Cooke,J.B.(1984)Progressinrock-filldams,(18thTerzaghiLecture),JournalofGeotechnicalEngineering,ASCE,110(10),1381–1414。
[2]Duncan,J.M.andChang,C.Y.(1970)Non-linearanalysisofstressandstraininsoil,JournalofSoilMechanicsandFoundationsDivision,ASCE,96(SM5),1629–1653。
[3]Duncan,J.M.,Byrne,P.,Wong,K.S.andBabry,P.(1980)Strength,
Stres-strainandBulkModulusParametersforFiniteElementAnalysesofStressesandMovementsinSoilMasses,inReportNo:
UCB/GT/80-01,UniversityofCaliforniaatBerkeley。
[4]Fitzpatrick,M.D.,Cole,B.A.,Kinstler,F.L.andKnoop,B.(1985)Designofconcrete-facedrock-filldams,Proc.oftheSymp.“ConcreteFaceRock-fillDams-Design,Construction,andPerformance”,JournalofGeotechnicalEngineering,ASCE,111,410–434。
[5]Hunter,G.andFell,R.(2003)Rock-fillmodulusandsettlementofconcretefacerock-filldams,JournalofGeotechnicalandGeoenvironmentalEngineering,ASCE,129(10),909–917。
[5]Khalid,S.,Singh,B.,Nayak,G.C.andJain,O.P.(1990)Non-linearanalysisofconcretefacerock-filldam,JournalofGeotechnicalEngineering,ASCE,116(5),822–837。
[7]Liu,X.,Wu,X.,Xin,J.andTian,H.(2002)ThreedimensionalstressanddisplacementanalysisofYutiaoconcretefacedrock-filldam,Proc.of2ndInt.Symp.OnFloodDefense,SciencePressNewYorkLtd,Beijing,Vol.1,pp.1494–1500。
[8]Marachi,N.D.,Chan,C.K.andSeed,H.B.(1972)Evaluationofpropertiesofrock-fillmaterials,JournalofSoilMechanicsandFoundationsDivision,ASCE,98(SM1),95-114。
[9]Marsal,R.J.(1967)Largescaletestingofrock-fillmaterials,JournalofSoilMechanicsandFoundationsDivision,ASCE,93(SM2),27–43。
[10]Nobari,E.S.andDuncan,J.M.(1972)Movementsindamsduetoreservoirfilling,Proc.ASCE,SpecialityConf.OnPerformanceofEarthandEarth-SupportedStructures,PurduUniv.,Lafayette,Indiana,pp.797–815。
[11]Saboya,F.Jr.andByrne,P.M.(1993)Parametersforstressanddeformationanalysisofrockfilldams,CanadianGeotechnicalJournal,30,690–701。
[12]Schanz,T.,Vermeer,P.A.andBonnier,P.G.(1999)TheHardeningSoilModel:
FormulationandVerification,Beyond2000inComputationalGeotechnics-10yearsofPlaxis,Balkema,Rotterdam。
(注:
本资料素材和资料部分来自网络,仅供参考。
请预览后才下载,期待您的好评与关注!
)
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 溷凝土 面板 堆石坝 沉降 分析