紫坪铺面板堆石坝三维有限元弹塑性分析.pdfVIP

现代水利水电工程抗震防灾研究与进展 紫坪铺面板堆石坝三维有限元弹塑性分析 孔宪京1，2邹德高1，2徐斌1，2周扬2刘京茂2 (1大连理工大学海岸和近海工程国家重点试验室，辽宁大连116024； 2大连理工大学水利工程学院，辽宁大连116024) 摘要：发展能够合理评价面板堆石坝在施工期、蓄水期与地震期应力与变形的数值 分析方法对于面板堆石坝的设计和安全运行具有重要的意义。本文开发了三维弹塑性 有限元分析程序，根据紫坪铺筑坝堆石料单调和循环荷载大型三轴试验确定堆石料的 广义塑性模型参数，对紫坪铺面板堆石坝的施工、蓄水过程及汶川地震震害进行了数 值模拟，并与实际施工监测和震害调查结果进行了对比。结果表明：数值计算得到的 大坝施工期、地震期的变形，地震后面板应力分布及施工缝错台均与实测结果规律一 致。表明采用考虑压力相关性的堆石料广义塑性模型进行面板堆石坝的施工期、蓄水 期与地震期的全过程模拟是可行的。 关键词：水工结构；弹塑性分析；三维有限元分析；紫坪铺面板堆石坝；全过程模拟 0引言 合理评价面板堆石坝在施工期、蓄水期以及地震过程中的应力和变形对于面板堆石坝的设 计和安全运行具有重要的意义。目前。针对面板堆石坝施工填筑及蓄水过程的应力和变形分析主 要采用邓肯E-B模型【ll，动力反应分析主要采用等效线性模型【2】，并在参数确定和计算程序开 发方面积累了丰富的经验。邓肯模型参数物理意义明确且易于通过室内三轴试验确定；等效线 性模型计算得到的地震时大坝加速度反应及应力分布较为合理【3】。但邓肯模型不能反应堆石体 的剪胀和剪缩特性，也不能合理考虑复杂的应力路径：等效线性模型不能直接得到大坝的地震 永久变形，需要借助滑动体分析方法【4J或基于应变势【5】的整体变形分析方法来计算大坝地震永久 变形。 模型)，P．Z模型具有许多优点，包括：不需要定义塑性势面函数直接确定塑性流动方向；不需 要定义加载面函数直接确定加载方向；不需要依据相容性条件直接确定塑性模量；可以考虑剪 胀和剪缩以及循环累计残余变形。此外，P-Z模型框架清晰，便于在有限元程序中实现，用一 套参数即可完成土构筑物的静、动力分析过程。 自P-Z模型提出以来，在粘土和砂土例方面得到了广泛的应用。一些学者根据原有模型的 思路，对原始模型做了进一步的改进。如考虑土的各向异性【lo】、主应力轴旋榭111、应力水平及 循环硬化对土的变形特性的影响【12J及将临界状态引入到模型中以反映较大围压和较大密实度范 围的土的变形特性【l¨41，颗粒破碎对堆石料变形特性的影响【15】以及在接触面上的应用推广【16】 等。 目前，P-Z模型在地下管线、地铁、加筋挡土墙和心墙堆石坝等方面均有所应用D7-23l。 P．Z模型不能适应较大的平均主应力变化范围，邹德高等【24】对P．Z模型进行了改进，在弹性 模量、加载模量和卸载模量方面考虑了筑坝材料的应力相关性。本文在GEODYNA[25】软件平 台上开发了改进广义塑性模型的土石坝三维静、动力计算模块，进行了筑坝堆石料单调和循 环加载大型三轴试验并确定了模型参数，对紫坪铺面板堆石坝施工、蓄水及地震过程进行了 三维有限元数值分析，对改进广义塑性模型用于高面板堆石坝施工期、蓄水期及地震期全过 程应力与变形分析的适应性进行了探讨。 现代水利水电工程抗震防灾研究与进展 1 广义塑性模型改进及参数确定 1．1广义塑性模型 根据塑性理论，总应变增量可以表示弹性应变增量和塑性应变增量之和 dE=d毛+dq (1) 其中，d乞为弹性应变增量；d‘为塑性应变增量。 应力应变关系可以表示为 dcr’D印：de (2) 弹塑性刚度阵可以表示为 Dep=De一等焉杀 ㈤ 其中，do。为有效应力增量，％为弹塑性刚度距阵，见为弹性刚度矩阵，ne,t．／u为加载或卸 载塑性流动方向，以为相当于屈服面法线方向，Hx：o为加载或卸载模量，其中“L”