土坝浸水变形分析论文.docx

土坝浸水变形分析论文摘要摘要：基于土坝浸水变形分析的增量有限元法并结合邓肯E－B模型编制了相应的土坝应力及变形有限元计算程序。应用实例的计算结果表明所采用的增量有限元分析方法是合理的所得计算结果也是接近实际的。摘要：浸水变形本构模型增量有限元法1引言众所周知土坝对水的功能十分敏感尤其在水库蓄水初期土坝会发生某些非凡的变形。比如水库初次蓄水期土坝上游坝壳表面往往产生向上游的位移这和水压力功能的方向正好相反同时还伴随一定的下沉而不是想象中的随浮托力的增大而上抬。大量的探究成果表明产生这种现象的根本原因就在于坝体土料的浸水变形或称湿化变形［12］。浸水变形往往使坝顶发生横向的伸长变位从而引起坝体发生不均匀沉降甚至导致坝体产生纵向裂缝［1］。因此如何对初次蓄水期土坝的浸水变形进行准确而合理地分析在此基础上把握其发生和发展的规律对确保土坝的正常运行具有重要意义。2浸水变形的分析方法目前有关浸水变形增量分析方法的探究还不够深入。但据笔者分析文献［2］所建议的一种增量分析方法即增量有限元法相对而言还是较为合理的。其基本原理如下摘要：设单元在浸水前的应力状态为｛σd｝假定它是由n级应力增量按比例增加达到的则每级应力增量为式中刚度矩阵［Dd］是和当前的应力状态有关的干态非线性弹性或弹塑性矩阵。将各级增量下的｛Δε｝累加即得浸水前的总应变｛ε｝。假定浸水前后应变相同则浸水后每级的应力增量可按下式计算摘要：｛Δσw｝＝［Dw］｛Δε｝（3）式中［Dw］为浸水饱和状态的刚度矩阵。将各级增量下的｛Δσ｝累加即得浸水后的总应力｛σ｝。按假想约束的思路可确定由浸水变形产生的“初应力”为｛Δσ｝＝｛σd｝－｛σw｝（4）然后将此假想的“初应力”约束释放转化为等效结点荷载即式中［B］为单元几何矩阵；负号表示由于浸水变形转化的等效结点荷载实际上不存在故需从单元结点荷载中予以扣除。由等效结点荷载｛F｝即可求得土体由于浸水变形所引起的附加位移及附加应变。为简化计算可将｛F｝和水压力或渗透压力、浮托力等所转化的结点荷载叠加在一起进行计算。3应力—应变本构模型目前非线性弹性的邓肯E—B模型是土石坝应力变形计算中最常用的模型［3］。该模型的基本原理如下摘要：材料的切线模量4应用实例某水库大坝坝型为砂砾石坝壳粘土心墙坝最大坝高127．5m坝顶长440m。上游坝坡1∶2．2下游坝坡1∶1．8。心墙两侧设反滤层上游侧宽3m下游侧宽5m。大坝心墙土料自1999年10月开始填筑2001年底大坝竣工、水库开始蓄水。基于上述的邓肯E－B本构模型及浸水变形的增量有限元方法笔者编制了相应的二维有限元计算程序用此程序对该水库初次蓄水时的大坝变形进行了有限元计算。计算取用的断面为大坝设计桩号为0＋226m的最大断面。计算中模拟逐级加载过程共采用17个荷载级其中施工阶段坝体填筑分11级水库初次蓄水从现有的几个观测水位直到正常高水位594．0m共分6级。计算所用的坝体材料模型参数见表1所列。心墙在竣工期及和各蓄水位相应的最大竖向变位的计算结果见表2此表中还同时列出相应的部分竖向变位观测成果以资对照。从图1、图2及表2可以看出随着水库蓄水水位的逐渐升高土体颗粒受到水的润滑进而产生越来越大的浸水变形从而使心墙产生越来越大的竖向变位；同时可以看出计算得到的竖向变位和实测结果较为接近且二者的变化趋向基本一致。根据位移计算结果计算得到的竖向位移均是铅直向下的这说明浸水变形所引起的下沉功能总大于浮托力所引起的上抬功能。从图3、图4及其它蓄水位时的水平位移计算结果可以看出上游坝壳土体的水平位移受浸水变形的功能较为明显在竣工期及较低蓄水位时由于无水平推力或水平推力尚小致使上游坝壳部分单元产生向上游的位移但随着蓄水位升高几乎全部上游坝壳的单元水平位移又变为向下游的位移。5结语本文基于土坝浸水变形分析的增量有限元法并结合邓肯E—B模型编制了相应的土坝应力及变形有限元计算程序。应用实例的计算结果表明所采用的增量有限元分析方法是合理的所得计算结果也是接近实际的。显然本文有关土坝浸水变形的探究思路和分析方法对类似工程新问题的探究具有重要的参考价值。参考文献［1］朱百里沈珠江．计算土力学［M］．上海摘要：上海科学技术出版社1990．［2］钱家欢殷宗泽．土工原理和计算［M］．北京摘要：中国水利水电出版社2000．［3］蒋国澄等．混凝土面板坝工程［M］．武汉摘要：湖北科学技术出版社19