基于振动孔隙效应的土体结构抗震稳定性分析.docxVIP

基于振动孔隙效应的土体结构抗震稳定性分析 1 基于地震滑移量的坝坝抗震稳定性评价方法 在土石坝抗灾设计的静态约束分析中，将以相同等效的静态约束法对整个坝的地震力进行积分，并使用土壤侵蚀稳定性的边界平衡方法计算系数。土石坝震害调查表明，拟静力分析方法并不能准确地评价坝坡的抗震安全性。实际上，一旦堤坝在地震作用过程中某一瞬间达到临界极限平衡状态，使得安全系数可能瞬间及其随后小于1，此时即称为超载，由于超载历时可能极短，堤坝边坡所产生的瞬时滑移或永久变形是有限度的，这并不能预示着坝体将完全失效，但是若整体滑移量或变形过大时，也将为正常使用所不容许。为此，Newmark于1965年建议，采用累计滑移量评价堤坝边坡的抗震稳定性，并提出了预测这种永久滑移量的刚塑性滑块模型。这种分析方法因其简便、实用而得到普遍接受，并且，在过去30多年中许多学者对这种简化的计算模型进行了改进和修正。尽管如此，这种以拟静力概念为基础的分析方法，既不能合理考虑坝体土料的动力非线性和滞回特性，也无法考虑震动中由于孔隙水压力的产生与上升对土的强度的衰减效应。 实际上，Newmark方法仅适用在地震荷载下动强度不会发生明显退化的土，如干砂、压实良好的粘土等。Sarma采用孔隙压力系数A和B近似地考虑孔隙水压力的影响；Wartman等采用十字板剪切试验测定破坏面上抗剪强度随剪切位移的变化关系，以此考虑振动中抗剪强度的退化效应，估算潜在滑动体的滑动位移。但是，这些改进方法尚未能够合理地考虑循环剪切荷载所引起的孔隙水压力上升而导致的土的强度的逐步退化效应。 为此，本文针对粘性土填筑坝，采用Yasuhara通过试验所得到的经验预测模型，考虑土体强度随振动孔隙压力上升所产生的渐进退化效应，将堤坝动力等价线性化有限元数值分析与潜在滑动体滑移量分析的Newmark模型相结合，提出了基于地震滑移量的堤坝抗震稳定性评价方法。在这种方法中，首先，基于动力等价线性化有限元分析确定坝坡潜在滑动土体的平均加速度时程；进而，基于拟静力概念采用极限平衡分析确定坝坡的安全系数随时间的变化历程，其中当安全系数瞬时小于1时，表明坝坡在地震中处于瞬时超载状态；采用Newmark刚体滑块模型估算瞬时超载所产生的滑动位移，将各个超载阶段的滑动位移叠加，求得设计地震动作用下堤坝边坡的累积滑移量。根据这种方法所进行的数值计算与分析表明,考虑强度循环退化效应后所得到的坝坡滑移量更为合理。 2 考虑到土壤强度的周期性衰减，水库地震滑动移量计算公式为 2.1 土的水动力响应分析 在地震动所产生的瞬时或循环荷载作用下，粘性土虽不会象砂土一样发生液化破坏，但由于在振动中将产生超静孔隙水压力，随着循环作用次数的增加，孔隙水压力不断上升，导致土的强度和刚度都会发生不同程度的渐进退化或弱化，从而，降低堤坝等土工建筑物与地基的抗震稳定性。 目前，关于计算中所采用的循环强度，一般直接选用静不排水强度的80%，或者根据Seed由土工循环试验所建立的强度预测关系曲线选取。并且认为，确定后的土的循环强度在地震作用过程中保持不变，从而，忽略了土的强度的动态变化效应，可能导致不合乎实际的分析结果。地震滑移量密切地依赖于土的循环强度及其可能的衰退。为合理地考虑这种强度退化效应，必须通过试验建立孔隙水压力增长模式与强度随孔隙水压力的退化模式。 由于地震动作用历时极短，因此，通常采用不排水分析方法。目前，一般基于不排水循环荷载引起的似超固结状态和由于应力释放引起的超固结等效的观点，通过土工试验建立粘性土动强度发展规律的预测模式[6,8～10]。Yasuhara在对现有研究工作进行系统分析与总结的基础上，给出了正常固结土在循环荷载作用下不排水强度的预测模式，以土的塑性指数为参数，通过三轴试验得到下列粘性土动强度与振动孔隙压力之间的依赖关系 式中su,cy,su,NC分别为循环不排水强度和正常固结条件下的初始静强度，(OCR)q为由孔隙水压力所引起的似超固结比， 式中 Ip为塑性指数；Cc与sC分别为压缩指数和回弹指数。 关于振动孔隙压力?u，这里直接采用Hyde所建议的下列经验模式 式中qr′为动偏应力的幅值；A,B,β为与固结比有关的系数，对正常固结土，A=-.1892,B=2.728,。 利用经验关系式（1）可以预测循环荷载作用后所导致的强度衰减，为了考虑地震过程中强度的渐进衰退，这里认为在循环荷载作用过程中的某一时刻，前期循环荷载及其作用次数所导致的循环强度衰减和前期累积孔隙压力也可简单地直接套用式（2），于是在动荷载作用过程中，某一时刻根据前期循环荷载的等效作用效应(包括循环应力与循环次数)就可以预测此时的不排水循环强度。将式（2）和式（3）代入式（1），并考虑荷载循环作用次数为N，则不排水循环强度比为 式中?u(N)与su,cy(N)分