面板堆石坝有限元分析的现状与展望

1、    面板堆石坝有限元分析的现状与展望            作者：高辉时间：2007-11-25 12:24:00                     摘要：本文通过天生桥一级面板堆石坝有限元变形计算结果和实测结果的对比分析,阐述了

2、面板坝有限元应力变形分析的现状，并展望了其发展趋势。 关键词：面板坝 堆石体 有限元 本构关系  1. 前言     现代碾压面板堆石坝具有工程量小工期短造价低等优点，已成为坝型选择的主要类型之一，在我国已成大规模推广的趋势。国内已建成的有关门山、柯柯亚、成屏、西北口、天生桥面板堆石坝等十余座，正在规划设计的有10余座，广泛分布于全国各地。20世纪50年代以前设计的堆石坝一般是不计算坝体应力应变的，但随着计算技术的发展和坝体的日益增高，同时在一些坝体中发现裂缝，坝体应力应变计算才逐步为人们所重视。    从“七五”

3、科技攻关开始，国内专家用不同的计算模型，考虑面板受力特点，采用相应的参数，对高面板坝进行二维、三维有限元分析计算。但是到目前为止，国内外对这种坝主要凭经验设计。由于堆石体材料特性十分复杂，有限元的计算结果与实际观测值之间还存在一定的差距。可见，在这方面尚有许多研究工作有待开展。2堆石体材料的本构关系    堆石体的本构关系表达了堆石体应力应变之间的关系，它无疑对应力应变的计算结果起决定性作用。目前建立堆石体本构关系往往基于已有模型，再针对堆石的力学特性确定甚至调整本构关系中各种材料参数。面板坝有限元计算时常采用非线性弹性模型和弹塑性模型。2.1非线性弹性模型非线性

4、弹性模型包括邓肯模型、内勒 模型、修正邓肯模型和邓肯模型等。这些模型最初是针对土、砂等一类材料在常规三轴试验基础上提出的，仅适用于二维分析计算，其中计算结果和已建坝实测结果较为符合的是邓肯E-B模型和内勒 模型7。模型18康纳(Kondner)根据常规三轴试验结果发现，关系可用双曲线拟合。随后，邓肯(Duncan)引入摩尔-库仑(Mohr Coulomb)破坏准则，导出切线弹模为式中, 为凝聚力；为初始弹模；为破坏比,是破坏时主应力差和应力应变双曲线的渐近线的纵坐标的比值,即/；为内摩擦角,对于粗粒土假定是的函数，即,式中和为试验参数，大气压(下同)。   

5、60; 詹布(Janbu)据试验指出初始弹模和侧限压力有如下关系：式中 、为两个由试验确定的参数 。 相应于卸荷再加荷情况的切线弹模为，按下式计算式中、是由试验确定的两个参数， 一般较为大。     1980年，邓肯(Duncan)提出切线体积模量的计算式式中、为试验常数。    体模相当于修正邓肯模型中的,考虑到只能在之间变化，故有限元计算中应限定在范围内，否则得不到合理计算结果。    平面应变条件下邓肯 模型的增量型本构关系为式中、为法向应力增量；为剪应力增量；、为法向应变增量；

6、为剪应变增量。     邓肯模型的不足之处在于它不能反映压缩与剪切的交叉影响；不能反映各向异性；不能反映加荷卸荷对的变化；不能反映中主应力对强度指标的影响等。但总的来说，该模型反映了堆石体变形的主要规律。它反映了非线性；用于增量计算时，能反映应力路径对变形的影响；通过与的差别部分体现了加荷卸荷对变形的影响；通过假定是的函数，在一定程度上考虑了高固结压力的影响。再加上该模型简便直观，概念明确，并积累了相当多的应用经验，因此目前被工程设计人员广泛使用。模型7内勒非线性弹性模型是模型中较为简单的一种，模型参数容易由常规三轴试验确定。内勒(Naylor)认为土的

7、切线体积变形模量随侧限压力的增加而增加，土的切线剪切模量随剪应力的增加而减小，在破坏时变为零。因此切线体积变形模量为切线剪切模量为式中，为初始体积变形模量，由常规三轴试验确定；为初始切线剪切模量，由常规三轴试验确定；、为常规三轴试验确定的参数。   内勒非线性弹性模型在面板坝有限元分析中较少采用。2.2弹塑性模型4    弹塑性模型把总的变形分成弹性变形和塑性变形两部分，前者用虎克定律计算，后者用塑性理论求解。南京水科院1987年提出了用于堆石体的双屈服面弹塑性模型，一般称为南水模型。该模型可以反映应力引起的各向异性和堆石的剪缩特性，在理论方面

8、有其优越性。但限于试验设备，该模型中某些参数如、较难从试验中得到,因此尚未进入工程实用阶段。3.非线性有限元求解方法16      一般来说，求解方法有迭代法和增量法两种，较为常用的是增量法。用增量法计算时，荷载逐级递增，可以模拟施工过程，计算结果也符合观测结果。以邓肯 模型为例，假设现进行第级加荷计算，按中点增量法其计算步骤如下：1根据前一级应力全量确定弹性常数和，并形成刚度矩阵；2加本级荷载增量的一半于结构，用下式求位移增量：           &#

9、160;        并计算应力和应变增量，进而累计求得应力全量，据此计算；3施加全荷，求位移增量：                          相应地可求出应力应变增量，累计则可得位移和应力应变全量；   重复上述步骤，可得各级荷载增量下的

10、解答。   还有一个问题，即新填筑层各单元初始应力状态是，如果以此计算，则切线弹模，就无法进行计算。克拉夫(Clough)等人将新填土层视为重液体处理，即(为单元形心距土表的距离，为填土容重)。目前较为常用的处理方法式计算是引入侧压力系数，即。事实上，新填土层经过反复碾压，可视为超固结土，可取前期预固压力(视碾压轻重可取MPa)。4算例 天生桥一级面板堆石坝位于红水河上游南盘江干流，坝高178m，上游坝坡11.4, 下游坝坡11.29。按平面应变问题计算。4.1单元剖分和计算参数   面板顶厚0.3，底厚为0.9m，沿面板厚度方向剖分一层单元，沿坝高分2

11、0层，共剖分500个单元，1486个节点。其中堆石单元459个，面板单元20个，接触面单元20个，趾板单元1个。单元形态为等参元计算采用中点增量法，按20级加荷方式进行，每级荷载又分三次施加。这样相当于按60级加荷计算，但又比剖分60层的情况节省计算时间。   混凝土面板和趾板采用线弹性模型， =18000MPa，=0.167。    堆石各分区对应于模型的计算参数见表1。表1  应力应变计算参数材     料干密度g/cm3垫层料2.2010500.35422100.2040.7064760.23650

12、.67细堆石料2.109700.36120000.2100.7604400.19352.58灰岩主堆石料2.109400.35019800.2000.8493400.185413泥岩主堆石料2.155000.25010500.1200.72725004810灰岩堆石料2.057200.30315500.1600.798800-0.1854134.2计算成果从面板坝有限元的计算结果来看，依据不同模型所得的应力结果是比较一致的，不同之处在于位移的大小和分布。而实际工程中，设计人员也更关心位移的大小和分布。图2和图3分别为竣工时铅直和水平位移等值线图。最大铅直位移为2.35m，约位于1/2坝高处；最大水平位移：上游为0.58m，下游为1.03m，参见图2。 5.计算结果和实测值的比较分析除本文算例的计算结果外，还将引用文献3和文献4对天生桥一级面板坝的计算结果，来主要说明邓肯模型和弹塑性模型在计算结果上的差异。     天生桥一级面板坝坝体变形观测，沉降观测结果规律性强，十分可信；水平位移的观测结果，国内外专家尚有不同看法。因此，下文主要针对铅直位移进行讨论。表2 天生桥一级面板坝有限元