土石坝抗震极限能力有限元分析的上限法.docVIP

  土石坝抗震极限能力有限元分析的上限法# 迟世春，杨昕光* （大连理工大学水利工程学院） 5 10 15 20 25 30 摘要：有限元极限分析一般将 Mohr–Coulomb 强度线性化，不适合具有明显的非线性特征 的土石坝极限分析。本文将极限分析的标准二次锥规划数学模型转化为对偶问题的“静力形 式”，可求解坝体内部应力，从而用迭代的方法确定材料在非线性强度准则下的等效剪切强 度，应用于土石坝极限抗震能力分析，得到该问题的上限解。 关键词：土石坝；极限分析；上限定理；非线性破坏准则 中图分类号：TV312；TV641.4 Upper bound method for aseismatic ultimate capacity analysis of rock-fill dams using finite elements CHI Shichun, Yang Xinguang (Shool of hydraulic engineering, Dalian university of technology, Dalian 116024) Abstract: The Mohr-Coulomb criterion is taken as a linearized yield surface in limit analysis using finite elements, but earth-rock material possess nonlinear strength properties which can’t be used directly by upper bound method. A static form, which is corresponding dual second-order cone programming (SOCP) problem of the upper bound limit analysis, is formulated. By solving the dual problem, the distribution of stress of slopes can be obtained. Then, the equivalent shear strength parameters of materials with nonlinear criterion can be determined by means of iterative solutions. So the method is applied to rock-fill dams with nonlinear criterion, the aseismatic ultimate capacity of dam can be determined. Key words: earth rock dam; limit analysis; upper bound theorem; nonlinear failure criterion 0 引言 汶川地震发生后，国家加强了水电工程抗震防震工作。要求对已建、在建及待建水坝进 行抗震复核，除在设计地震工况下满足“可修复”的要求外，增加了在校核地震工况下满足 “不溃坝”的要求，对特别重要的挡水建筑物，还应研究极限抗震能力和地震破坏模式。但 如何进行土石坝抗震极限能力分析，目前国内外还没有一致的看法。 塑形极限分析具有严密的理论基础，它不考虑结构受力过程中的塑性演变，直接求解结 构的极限承载能力[1 、2] ，已经结构分析的有力工具。在塑形极限分析中，为了计算简便，材 料的剪切强度一般由黏聚力及内摩擦角两个参数描述。这种破坏准则由于与应力无关而被称 35 40 为线性破坏准则。但根据大量的实验表明，筑坝土石料的强度准则远非线性，而是具有应力 相关性。由于基于有限元的极限分析在数学上是线性规划问题，一切方程均需描述成线性函 数，这给土石坝极限抗震能力的分析带来了困难。 我们在 Makrodimopoulos 和 Martin 等人[1]的研究基础之上，对坝坡抗震极限能力问题进 行上限有限元分析，通过对偶问题的“静力形式”迭代求解材料非线性破坏准则下的等效剪切 强度参数，得到问题的上限解。 基金项目： 高等学校博士学科点专项科研基金 20090041110016 作者简介：迟世春（196408），男，教授，主要研究方向：土石坝地震工程. E-mail: schchi@ - 1 -   1 极限分析上限理论 设存在一个完全刚塑性结构 V ，其边界为 S ，则当结构达到极限状态时，必定存在一允 许速度场 u ，使得内能耗散不大于外力做功，即  V T (1) 45 其中，