动孔压试验曲线法在土坝地震反应中的应用.docxVIP

动孔压试验曲线法在土坝地震反应中的应用 0 土体液化判断方法 当沉积物土和沉积物土（包括地基）受到地震的影响时，内部孔水压（以下简称“孔压”）随着振动时间的推移而逐渐增加，有效荷载相应减少，动切割模型和循环切割强度相应降低。当动孔压继续增长到某个程度时，堤坝可能出现液化或滑坡现象，其中液化是地震作用引起的最主要的震害之一。研究表明，堤坝是否发生液化与土体类型、初始应力状态、动荷载及排水条件等因素关系密切。目前判断土质堤坝是否发生液化的常用方法有剪切波速法、标贯法、临界孔隙比法及剪应力对比法等，其中剪应力对比法有较为系统的理论基础，实际应用较多。但剪应力对比法中的液化抗剪强度并不一定是真正液化时的抗剪强度，而可能是采用动三轴试验中轴向应变达到某一特定值时的剪应力强度。另外，剪应力对比法是一种总应力判断方法，没有考虑孔隙水压力对剪切模量等的影响。因此，剪应力对比法用于液化判断有明显的理论缺陷。笔者认为，比较理想的液化计算方法应该不仅在有效应力框架内考虑动孔压变化的影响而且还能简便合理地计算和判定是否发生液化。 此外，滑坡也是土质堤坝地震时的主要震害之一。地震时，堤坝的边坡是否能保持稳定与边坡土体中动孔压有密切关系，即与土体的液化度相关性很大。一般而言，在液化度未达到100%之前，堤坝边坡就有可能发生大面积塌滑。因此，在研究土质堤坝地震反应时，不能遗漏边坡动力稳定性的分析。本文采用考虑动孔压影响的液化判定方式和边坡动力稳定计算模式，对某土坝进行了非线性有效应力动力有限元法计算，对大坝可能液化范围和坝坡动力稳定进行了预测，并对所采用的抗震加固措施进行了抗震效果的计算评价。 1 动孔压试验曲线法 用有效应力法进行土质堤坝液化分析的关键是动孔隙水压力的计算，但孔隙水压力的增大、扩散和耗散是一个非常复杂的过程，不仅受土体的物理特性、颗粒级配、填筑密实程度、透水性能、初始应力等内部因素影响，而且还与地震级别、震中距离、振动持续时间等外部条件有关。虽然已有不少学者提出了许多孔压计算模型，如应力模型、应变模型、能量模型、内时模型等，但这些计算模型都做了一定的假设和简化，往往只能考虑部分主要影响因素，而且参数选取是否合理对计算结果影响很大。事实上，土体的液化过程是一个灰色过程，很难采用一个计算式对上述主要影响因素进行综合考虑。如果采用合适的试验方法，合理模拟土体液化过程，并对试验成果采用适当的方式进行整理表达，则可以较为全面地考虑上述主要的液化影响因素，尤其是那些不可“量化”的复杂因素。在此基础上，建立孔压计算方法，则能较为全面地反映液化前孔压动态变化的综合过程。顾淦臣提出的动孔压试验曲线法就是基于上述思想的一种较为理想的动孔压计算方法。利用该方法进行动孔压计算的基本原理和主要步骤为： (1）对符合设计标准的土样进行不同初始应力状态下的动三轴试验，整理出不同振动次数（N）下的动孔压比与动剪应力比的关系曲线（见图1）。图1中，横坐标为破坏面上动剪应力τd与平均有效应力σm′之比，即动剪应力比αd；纵坐标是动孔压力ud与围压σ3之比，即振动孔隙水压力比，以下简称“动孔压比”；Kc为固结应力σ1与σ3之比。 (2）根据坝址设计震级，确定地震历时t和等效振动次数Neq。将地震历时t分成若干时段△Ti，相应第i时段的振动次数为△ni=Neq△Ti/t，则到第m时段时累计振动次数为nm= (3）根据地震反应计算得到的至第m时段为止各单元破坏面上动剪应力时程曲线峰值τmax的65%作为平均动剪应力τav，除以该时段的平均有效应力σm′后，得到动剪应力比αd=τav/σm′。根据该时段累计振动次数，从相应振动次数的图1曲线查取各单元动孔压力比（或从不同振动次数曲线间插值），然后乘以σ3即可得该时段孔压力ud。 (4）利用求得的孔压ud，对第m时段各单元的有效应力σm′进行修正，重新计算该时段的剪切模量和阻尼比，并判断是否满足收敛准则，如不满足，迭代至满足，然后进入第m+1时段动剪模量及固结应力比等的计算。重复上述过程直至地震结束，即可得到整个地震历时过程中大坝各处的孔压时程曲线和大坝的孔压等值线分布情况。 值得一提的是，图1中动孔压比与动剪应力比试验曲线是综合考虑土体振动过程中各种主要内外因素得到的，概念清晰，计算结果与实际情况较为符合。同时，该方法又避开了如取轴向应变5%时的动剪应力作为液化剪应力的任意性。通过试验曲线查取动孔压并进行液化判断分析的方法直观简捷，容易被广大工程人员接受。 2 等效线性粘弹性模型 当前较为合理地反映土体动应力动应变特性的动力本构模型主要有双曲线模型、粘弹性模型和理想弹塑性模型。其中，等效线性粘弹性模型因概念明确，参数确定方便，应用较为普遍。此类模型中，最常用的当属哈定-德涅维契模型，其主要计算公式如下： 式中，γr为参考剪应变，由