土石坝拟静力抗震稳定分析的强度折减有限元法.doc





《土石坝拟静力抗震稳定分析的强度折减有限元法.doc》由会员分享,可在线阅读,更多相关《土石坝拟静力抗震稳定分析的强度折减有限元法.doc(25页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、【精品文档】如有侵权,请联系网站删除,仅供学习与交流土石坝拟静力抗震稳定分析的强度折减有限元法.精品文档.土石坝拟静力抗震稳定分析的强度折减有限元法 第31卷第5期 岩 土 力 学 Vol.31 No.5 2010年5月 Rock and Soil Mechanics May 2010 文章编号:10007598 (2010) 05150307 李 湛1, 2,栾茂田2 (1.中国建筑科学研究院地基所,北京 100013;2.大连理工大学 土木水利学院岩土工程研究所,大连 116024) 摘 要:基于拟静力抗震设计概念,提出利用强度折减有限单元法分析土石坝的抗震稳定性,给出了两种确定地震惯性力
2、的方法:(1)依据水工建筑物抗震设计规范1,并结合有关土石坝动态分布系数计算了沿坝高分布的地震惯性力;(2)直接利用土石坝有限元地震动力反应分析得到的单元节点加速度反应,依据建议的方法确定坝体各单元节点的地震惯性力。将上述计算确定的地震惯性力与其他形式的外荷载共同作用到土石坝上,采用强度折减有限元法确定土石坝坝体的拟静力抗震安全系数。对于稳定渗流期,水位降落期等不同工况,或需要考虑振动孔隙水压力作用的饱和无黏性土填筑坝等不同计算条件,给出了使用折减强度有限元法分析坝体抗震稳定性的实现途径和方法。研究表明,有限元法对边界条件、复杂断面条件和材料分区及荷载组合均具有较强的适应能力,因此,使用有限元
3、法分析土石坝抗震稳定性具有显著的优越性。 关 键 词:土石坝;拟静力抗震稳定性;地震惯性力;折减强度有限元法;弹塑性有限元法 中图分类号:TV 641 文献标识码:A Pseudo-static stability calculation of earth and rockfill dams using shear strength reduction finite element method LI Zhan1, 2 , LUAN Mao-tian2 (1. China Academy of Building Research, Beijing 100013, China; 2. Insti
4、tute of Geotechnical Engineering, School of Civil and Hydraulic Engineering, Dalian University of Technology, Dalian 116024, China) Abstract: An analysis procedure based on shear strength reduction finite element method for pseudo-static seismic stability calculation of earth and rockfill dam (ERFD)
5、 is presented. In the proposed procedure, two methods are proposed to determine the seism-induced inertia force: (1) the inertia forces are computed with the same method adopted in the pseudo-static limit equilibrium slices method; (2) with the absolute acceleration response of element nodes obtaine
6、d from the finite element dynamic response analysis of the dam, the inertia forces corresponding to the node acceleration are computed with the proposed empirical method. With the inertia forces achieved above, the self-gravity, and/or the forces corresponding to the seepage force or dynamic pore wa
7、ter pressure, the elastoplastic finite element analysis is conducted to achieve the pseudo-static factor of safety of the dam by reducing the effective shear strength parameters of soil until the limit equilibrium is reached. Finally, numerical computations are conducted to illustrate the proposed a
8、nalysis method; and the results are compared with those obtained from the pseudo-static limit equilibrium method. It is concluded that the utilization of the finite element method for aseismic stability analysis of ERFD is feasible. Key words: earth and rockfill dams; pseudo-static seismic stability
9、 analysis; strength reduction finite element method; elastoplastic finite element analysis 1 引 言 目前已建立的土石坝抗震稳定分析方法主要有 静力法、拟静力极限平衡法、剪切条法、有限元动 力分析方法等,其中拟静力极限平衡法是我国土石 坝抗震设计规范中推荐的主要设计方法1,在多年的工程实践与理论研究中其理论、方法及评价体系等方面都已经较为成熟与完善,且积累了大量的工程与理论经验。拟静力极限平衡法一般采用某种动力分析方法,如剪切条法确定地震加速度分布系数沿坝高的分布,将其与土体的重量相乘,得到地震惯性力的
10、大小,通常还要在上述地震系数的基础上收稿日期:2008-10-22 基金项目:国家自然科学基金资助项目(No. 50179006);教育部跨世纪优秀人才培养计划研究基金资助项目和中国科学院武汉岩土力学研究所前沿领域基础研究基金资助项目(No. Q110305)。 第一作者简介:李湛,男,1975年生,博士,主要从事岩土地震工程、地基与基础方面的科研和生产工作。E-mail: lz-xj 1504 岩 土 力 学 2010年 乘以一个修正系数,以弥补实际结构同理想结构之间动力反应差异所带来的影响。从上述分析可以看出,拟静力极限平衡法在确定地震力进行拟抗震稳定分析中存在如下的问题:(1)只给出了地
11、震力沿坝高的分布,不能考虑同一坝高不同位置处地震惯性力的不同,因此不能全面地反映地震惯性力在坝内的非均匀分布;(2)没有区分坝体几何形状、坝体类型、土体材料特性等的不同,以及坝体-地基相互作用的影响等,对所有土石坝都采用了相同的地震加速度动态分布系数的分布模式;(3)忽视了沿坝高不同位置处地震惯性力方向的不同,假定地震惯性力全部指向上游或下游;(4)不能考虑土体应力-应变特性,土体的强度-变形-稳定之间的联系等。 近年来,先后发展了有限元、刚体弹簧元、拉格朗日有限差分等理论上更为严密的数值方法,这些方法也被用来进行边坡的稳定计算与分析。在各种现代数值分析方法中,有限元法应用最为广泛和成熟。在边
12、坡稳定性分析中,有限元方法主要应用在以下两方面:(1)首先通过有限元数值计算确定给定荷载作用下边坡中的实际应力场,由此可以确定沿某一潜在滑裂面上的正向有效应力与剪应力,并进一步地确定沿该潜在滑裂面失稳时边坡的某种意义上的整体安全系数2,进而按照极限平衡法原理,通过搜索或试算确定临界滑裂面及其相应的安全系数。(2)将Bishop的安全系数定义3与弹塑性有限元法相结合,按照某一系数逐步将土的实际抗剪强度参数进行折减,基于折减后的抗剪强度进行一系列的弹塑性有限元分析,直至所选定的抗剪强度折减系数使得边坡达到临界失稳状态,与之相应的抗剪强度折减系数与Bishop安全系数定义中的整体安全系数是一致的,这
13、类方法通常称为强度折减弹塑性有限元方法4 6 地震惯性力沿坝高的分布;(2)利用非线性有限元动力反应分析确定的单元节点的绝对加速度反应,按建议的经验方法直接计算单元各节点的地震惯性力。将上述计算得到的地震惯性力与土的自重、渗透作用力及振动孔隙水压力等荷载一齐作用到坝体上,按照折减强度有限元方法计算坝体的拟静力抗震安全系数。针对具体算例进行了计算与分析,并将所得到结果与拟静力极限平衡法进行了比较,给出了一些有益的结论及建议。 2 土石坝稳定分析的强度折减弹塑 性有限单元法 将土的抗剪强度参数c和同时折减某一个相同的系数F,得到一组新的强度参数: =cm ctan? =arctan? , m? (
14、1)FF? 和m为基以如此折减后的抗剪强度参数cm 础,根据理想弹塑性理论进行弹塑性有限元数值分析,并逐步判断边坡是否发生失稳或者达到极限平衡状态。一般先从某一个较小的假定折减系数F开始,逐渐增大F,进行弹塑性有限元计算,直到边坡达到极限平衡状态,此时所对应的强度参数折减系数F就是边坡的安全系数Fs。 在本文的弹塑性有限元数值计算与分析中,采用下列的Mohr-Coulomb准则作为屈服准则: f( m,q,)= ?sinsinq? 3msin?3ccos=0 (2) 式中:m、q分别为平均应力和广义剪应力;为应力Lode角。同时采用非关联流动法则,所选用的塑性势函数与屈服函数在形式上一致,而差
15、别在于将式(2)中摩擦角由剪胀角代替。当=时,即相当于关联流动法则。 在这种弹塑性有限元方法中,目前已经提出了各种用于临界失稳评判的标准7 8 ,这种方法无需人为事先 假定破裂面,也无须对于临界滑裂面进行随机搜索或多次试算等过程。用有限元方法来进行边坡稳定性计算与分析,具有如下的优点:(1)能够考虑土体的应力应变特性,例如材料的硬化、软化和渐进破坏等;(2)适用于具有马道等结构的复杂形状、具有复杂材料分区或边界条件的土工建筑物稳定分析;(3)能够了解强度发生恶化时土工结构的渐进破坏过程等。 基于上述分析,本文将基于折减强度的弹塑性有限元分析方法应用于土石坝拟静力抗震稳定分析。地震惯性力的确定使
16、用了两种方法:(1)采用水工建筑物抗震设计规范1建议的方法,计算 ,本文采用坡顶 水平位移发生突变时所对应的状态作为临界失稳状态,与之相对应的抗剪强度折减系数作为边坡的整体安全系数。 3 基于规范法的土石坝拟静力抗震稳 定分析 3.1 地震惯性力 土石坝抗震稳定计算最早采用静力法,即对坝体全断面都采用相同的地震系数,将地震系数与土 体重量相乘后作为相应的地震惯性力,然后按照极限平衡条分法计算土石坝的抗震安全系数。我国水工建筑物抗震设计规范规定采用拟静力法进行土石坝抗震稳定分析,与静力法不同,拟静力法吸取了剪切条等动力分析方法关于坝体动力响应计算的成果,考虑了土石坝地震加速度反应沿坝高呈上大下小
17、分布这一实际情况,并用坝体动态分布系数来反映这一调整的概念。水工建筑物抗震设计规范对于土石坝,规定了质点i的动态分布系数i沿坝高的分布形式,如图1所示,其中坝顶的最大动态分布系数m在设计烈度度时分别取3.0、2.5和2.0。 3.2 土石坝拟静力有限元抗震稳定分析 由密实的黏性土、堆石等土料填筑的土石坝,在地震作用过程中不产生或仅产生很小的振动孔隙水压力,对于这种情况,在土石坝抗震稳定计算时,可不计算振动孔隙水压力的影响,只需计算稳定渗流期或水位正常降落情况下的静孔隙水压力。此时,可以根据有限元渗流计算得出的静孔隙水压力或渗透力,依据有限元理论计算与其相应的等价节点力9。 用弹塑性有限元数值分
18、析计算土石坝的拟静力抗震安全系数时,所需计算的荷载包括土体自重、地震惯性力荷载、与渗流作用力相应的节点力荷载以及其他静力计算需要考虑的荷载形式。有限 m1)/3 元分析中采用土体的静力抗剪强度参数。 3.3 土石坝拟静力总应力法及有效应力法有限元 抗震稳定分析 对于由饱和无黏性土、饱和砂卵石等土体建造的土石坝,在地震作用下土体中产生振动孔隙水压力并逐渐累积,造成土体有效抗剪强度的降低,并 降低坝体的稳定性,此时应采用适当的方法估算坝体各部位的振动孔隙水压力,并在坝体抗震稳定分析中考虑孔隙水压力的影响,可以通过总应力法或有效应力法测量土体的抗剪强度指标,相应地,可以采用基于地震总应力抗剪强度指标
19、或地震有效应力抗剪强度指标的弹塑性有限元分析确定坝体的抗震安全系数10。 采用地震总应力法分析坝体的抗震稳定性,需要采用总应力法测量的动力抗剪强度指标,可以根据震级先确定出场地的等效振动次数,然后,通过振动三轴试验测量单元主应力比与动总应力内摩擦角之间的关系曲线。在使用折减强度弹塑性有限元法分析之前,先对坝体进行有限元静力计算,确定出坝体单元的主应力比,然后由上述主应力比与动总应力内摩擦角的关系曲线确定坝体各单元的地震总应力内摩擦角,根据折减强度弹塑性有限元法确定坝体的拟静力抗震安全系数。 采用有效应力法进行抗震稳定分析,需要先由振动三轴试验测量动孔隙水压力与振动加速度及固结主应力比之间的关系
20、曲线。根据图1的坝体动态分布系数i和设计水平加速度代表值h计算坝体不同高程的绝对加速度,结合土石坝有限元静力分析得到的各单元的固结主应力比,由上述动孔隙水压力-振动加速度-固结主应力比关系曲线确定坝体各单元的振动孔隙水压力,将由此确定的孔隙水压 力转化为等价的节点力荷载,并与其他形式的荷载 (a) 坝高H40 m (b) 坝高H>40 m 图1 土石坝动态分布系数i Fig.1 Dynamic distribution coefficient i for earth and rockfill dam 采用拟静力法计算地震作用效应时,沿土石坝高度作用于质点i的水平向地震惯性力为 Fi=ah
21、GEiaig (3) 式中:Fi为作用于质点i的水平向地震惯性力;ah为水平向设计地震加速度代表值,设计烈度为度时分别取0.1g、0.2g与0.4g,g为重力加速度; 为地震作用效应的折减系数,一般情况下取为 0.25;GEi为集中在质点i的重力作用标准值。 采用式(3)确定的地震惯性力具有如下特点:在同一坝高的水平截面上,地震惯性力是均匀分布的;从坝基到坝顶地震惯性的作用方向都相同。由此确定的地震惯性力为体积力,根据有限元法的原理,可以将各单元的地震惯性力转化为作用在单元节点的等价节点力荷载,将与地震惯性力等效的节点力荷载与其他形式的荷载一起施加到坝体上,采用折减强度弹塑性有限元数值计算即可


- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 土石 静力 抗震 稳定 分析 强度 有限元

限制150内