拱坝的破坏分析及超载问题探讨

1、拱坝的破坏分析及超载问题讨论摘要：以拱坝和周围的岩体为研究对象，采用线弹性本构模型和不同的破坏追踪方式，对破坏过程中的不同的应力状态应用不同的强度准那么和功能函数，破坏过程的分析按三维随机有限元法分析可靠度指标场，从可靠度指示场中寻找并联络统中最易出现的破坏轨迹，给出其破坏形式和平安度。同时按逐步增加水的比重和进步水位两种方法给出拱坝的超载破坏过程，最后指出拱坝的超载是不现实的，也不能给出切合实际的拱坝平安性。关键词：破坏分析超载可靠度拱坝拱坝的超载才能很强，但对拱坝超载才能的研究始于20世纪60年代出现了电子计算机和有限元法之后，李新民提出用机动法分析拱坝的超载才能，他的做法是在直法线假定的

2、前提下求弹性厚壳的极限荷载，实际上，厚壳中还存在着薄膜应力，仅用铰链来描绘破坏点的受力状态，过于简单化。超载试验法是另一个拱坝极限荷载的研究方法，通过模型试验，对拱坝逐渐加载到超过正常水压几倍的荷载(或利用高密度的液体)，通过实际观测，可以得到拱坝裂缝开展的详细情况，同时，得到拱坝从初始开裂直至彻底破坏过程中任一破坏程度的超载系数，给出一个宏观的平安指标。但是，在拱坝的正常运行中，超载是不现实的，并且拱坝的真正危险不是荷载的成倍增加，而是材料强度的缺乏，因此，就不能用这个指标来衡量正常运行拱坝的平安度，并且也不能反映随机性。极限强度分析中，水压保持不变，减小材料的强度，直到拱坝出现破坏，极限强

3、度分析是从设计材料的强度剩余上观察拱坝的“超载才能，却看不到破坏的开展过程。通过数学模型和拱坝的弹塑性应力分析，可以建立一个拱坝从初始开裂直至彻底破坏全过程的数学模型。一般在上述分析的过程中把参与分析的参量视为定值，可以得到一个逐点破坏的平安系数序列和破坏的迹线，虽然在这个破坏的过程中每点的破坏并不是独立的，但是这个平安系数序列却是离散而独立的。因此，这种方法也不能给出一个破坏程度的指标。本文阐述了一种全新的关于拱坝破坏分析的思想：首先建立拱坝断裂的随机数学模型，然后进展数学模拟分析，在模拟分析的过程中把参与分析的主要参量取为随机变量，进展随机数学分析，那么可以记录诸点依次破坏的概率及条件概率

4、，得到拱坝从出现裂缝到每条裂缝开展稳定或不稳定溃坝的条件概率序列拱坝破坏轨迹的概率向量。这个由条件概率组成的概率向量中的任何一个分量(不妨定义为线失效概率)值都有明确的物理意义，即这个分量代表了一种破坏程度的概率。最后得到开裂稳定或到达溃坝的概率，即拱坝沿最大可能失效途径的失效概率和可靠度。但是它不同于构造的体系可靠度，体系可靠度是指构造所有主要失效形式的结合失效概率，但对于拱坝这种大型的体系来说计算其体系可靠度还不大可能，因为拱坝是一复杂的空间受力体系，每个点有不同的破坏形式，并且其节点数也非常多，因此其失效途径非常多。本文的目的在于研究拱坝一种较大可能的破坏过程，以研究溃坝的机理，并同工程

5、上常用的超载试验相比拟，分析它们之间的不同。1拱坝破坏轨迹的概念按连续体力学的概念，拱坝的破坏就是无限个不连续点的集合，然而用离散的有限元法分析连续体力学问题时，那么拱坝的破坏过程就应当是有限个不连续点的集合，按三维有限元法分析拱坝的效应尝功能函数场和可靠度指标场，记录可靠度指标的最小值及相应的失效概率，在指定可靠度指标最小单元破坏的前提下，重新计算与该单元有关的刚度矩阵，组装总刚，再次计算拱坝的效应尝功能函数场和可靠度指标场，再次记录可靠度指标的最小值及相应的失效概率和条件概率，诸点进展下去，直到破坏个单元后拱坝成为机构或开裂稳定，这个破坏单元的几何位置就构成了拱坝的破坏轨迹，这个条件概率的

6、数值就构成了描写这个破坏轨迹的概率向量Pf=pf1pf2pf(1)式中概率向量中的第k个分量为：(2)当k=1时为拱坝的初始开裂的概率；k=为开裂稳定(或溃坝)的概率；k=n(1n)为开裂程度为n的概率。假如将拱坝的开裂视为一个并联体系，那么这个轨迹是并联络统中最易出现的一条。上述算法没有考虑其他的失效途径，因此得不出构造的体系可靠度，但是点的可靠度计算仍然是有意义的。因为下一步的可靠度是在上一步最小可靠度单元破坏条件下的条件可靠度，在上一步向下一步转化的过程中，由于有单元失效，因此其效应场已发生了重分布。转贴于论文联盟.ll.2拱坝上任一点的破坏形式及准那么函数2.1拱坝上任一点的破坏形式根

7、据混凝土三轴受力情况，过镇海教授将混凝土的破坏形式归结为：拉断；片状劈裂；柱状破坏；斜剪破坏；挤压流动，如表1。当出现拉断和片状劈裂破坏时，该点降为二维应力状态，继续以二维应力状态参与拱坝的应力分析；当出现柱状破坏时，该点降为一维应力状态，继续以一维应力状态参与应力分析。在二维应力状态下还可以出现进一步的破坏，其破坏形式可以根据一点二维应力状态进展分析，当出现斜剪破坏和挤压流动时，该点将失去承载才能。表1典型三维破坏形式的划分破坏特征拉断柱状破坏片状劈裂主导应力132，3应力状态单轴1030二轴1201/30.51/30.052/30.22/30.2三轴12301/3，2/30.11/30.1

8、12301/3，2/30.11/3，2/30.151/30.11/3，2/30.15建议的分析方法按局标中单向弥散裂缝，物理正交异性单元处理按局标中双向弥散裂缝，物理正交异性单元处理按局标中单向弥散裂缝，物理正交异性单元处理注：表中T为拉应力，为压应力。2.2准那么函数在拱坝的破坏过程中，点的破坏准那么应当根据点的实际受力状况采用相应的破坏准那么。详细地说，对于未破坏而处于三维应力状态的单元，应当用混凝土材料的三轴破坏准那么，对于破坏后降为二维应力状态的单元，应当用双轴破坏准那么，对于破坏后降为一维应力状态的单元，应中选用单轴破坏准那么。对于三维应力状态本文选用俞茂宏双剪强度准那么。(3)(4

9、)(5)二轴一压一拉：(6)二轴受拉：g()1t-ft,g()2t-ft(7)g()-ft(8)为压应力时：g()-f(9)2.3弥散裂缝拱坝的诸点破坏过程中，每当出现一个开裂破坏点后，分析中用弥散裂缝来代替单独的裂缝。当出现拉断或片状劈裂时，弹性矩阵如下(10)式中：,，为抗剪残留系数(0.10.5)。当出现柱状破坏时，弹性矩阵如下：(11)式中：(i=4，5，6)。当出现斜剪破坏或挤压流动时，弹性矩阵为0。上述弹性矩阵都是在各自的部分坐标系中，对于不同的单元，由于主应力方向不同，部分坐标系也就不同，在进展构造矩阵的组装时，应进展坐标转换。转贴于论文联盟.ll.3随机分析开裂点的分析是基于点

10、的失效概率分析。点失效概率分析的一次二阶矩方法如下：设可靠度分析的根本变量X=(x1,x2,xn)为一组互相独立的变量，假设某个变量为非正态变量，那么当量正态化为正态变量。令Yi=(xi-xi)/xi，那么Y(y1,y2,yn)为互相独立的标准正态变量。记Y=T(X)。设功能函数为g(X)=g(R(X),S(X)，将其转化为标准正态空间为G(X)=g(R(T-1(Y),S(T-1(Y)，按迭代格式(12)某双曲混凝土拱坝，最大坝高292，坝顶长度805，坝顶厚度12，坝底厚度72.9。坝体采用8结点六面体单元(360个)，沿岸坡为6结点五面体单元(80个)充填，沿坝厚分为5层，地基的单元全部是

11、8结点六面体单元，沿建基面向周围扩散为4层，坝和地基的总结点数为1434个，有限元模型图如图1。计算工况为：水库正常蓄水位(1240.00)及相应的尾水位(1010.00)和设计正常温降，自重，泥沙压力(淤沙高程1120.30)。图1有限元模型考虑混凝土拉、压强度的随机性，取混凝土的抗拉强度均值为2.95Pa，方差为0.59Pa，抗压强度均值为31.0Pa，方差为6.2Pa。4.1分析思路拱坝的破坏分析中采用3种方法，第一种是按可靠度最小进展破坏追踪，即假定可靠度最小的单元破坏，计算效应场，然后再假定可靠度最小的单元破坏，直到拱坝成为机构为止；第二种方法是逐渐增大水的比重，直至拱坝破坏；第三种

12、是逐步进步水位直到拱坝破坏。4.2结果分析4.2.1破坏过程分析拱坝的破坏过程如图2、3、4所示(第一层为上游，第五层为下游，图2中的数字为破坏的单元序号；图3中的数字为水密度每增大0.5倍后破坏的单元，即1为正常水密度，2为1.5倍水密度，3为2倍水密度，以此类推；图4中的数字为水位每增高50后的破坏单元，即1为正常设计水位，2为正常设计水位进步50，以此类推)。按可靠度最小进展破坏追踪分析拱坝的破坏过程是拱坝的上游面坝踵区首先开裂，裂缝向深处和沿周边开展到一定的程度，拱坝下游坝址区破坏。当坝底裂缝贯穿以后，裂缝沿上游面周边继续向上开展直至坝顶，这样上游面形成周边缝。在此过程中下游坝址区的裂

13、缝沿周边向上也有所开展，但范围不大。然后是下游左岸大面积的开裂直至破坏。按水比重增加进展溃坝分析的破坏顺序也是首先上游坝踵区先裂，随着裂缝向深处和沿周周边向上开展，下游坝址区破坏，随着荷载的增加裂缝按上述顺序诸层向上开展，直至破坏。破坏时上游图2按最小可靠度指标追踪溃坝分析图3按水比重增加进展溃坝分析图4按水位进步进展溃坝分析面大部分单元破坏，其他层破坏的单元数也很多，上下游都形成了周边缝，但上游面的形成先于下游面。按水位进步进展溃坝分析也是首先上游坝踵区先裂，随着裂缝向深处和沿周边向上游开展，下游坝趾区破坏，但是裂缝沿周边开展的速度要比前两种要快，随着水位的加高，周边缝上下游贯穿，同时上游坝

14、面中部大部分范围破坏。4.2.2破坏过程的应力场分析图5给出了拱坝破坏过程中的坝体最大及最小主应力，按可靠度最小进展破坏追踪分析时随着开裂单元的增加，开场坝域中的应力状态一步步地得到改善，最大主拉应力降到3Pa以下，最小主压应力大幅度减小，当破坏到27个单元时最小主压应力降到-33Pa，这恰恰迎合了混凝土材料适于抗压、不适于抗拉的材料性质，而这正是拱坝构造的优越所在它的开裂使应力重分布的趋势向着拱的受力性能转化，似乎拱坝超载才能的本质就在这里体形开场的约束并非全是合理的，随着某些单元的破坏，应力的调整，逐渐解除了那些不合理的约束，拉应力的数值和区域都减小，压应力的绝对值和区域都增大。然后随着破

15、坏的开展最小主压应力上升，坝体内的应力情况正常。按其他两种方法分析时坝体最大主应力是随着单元的开裂应力逐步增加，特别是图5三种破坏方式的坝体最大主应力比照图6拱冠梁顺河向位移图7拱坝破坏过程的破坏概率、破坏可靠度指标及指定破坏单元的可靠度指标从以上应力场及位移场的分析可以看出三种方法分析的结果是极不一样的，虽然后两种方法能给出工程上易于承受的平安概念，由于在拱坝的实际运行中，无论何种方式的超载都是极不现实的，拱坝的破坏在除去根底沉陷、坝肩失稳等因素的影响外，主要是材料强度的缺乏，特别是随着拱坝的运行某些区域的混凝土力学性质要逐步下降，第一种方法正是按照这一思路逐步指定可靠度最小单元破坏而实现破坏追踪。4.3可靠度指标及平安评价从以上应力、位移效应场看出按可靠度最小进展破坏追踪分析拱坝的破坏是现实可行的，图7给