基于Wolffersdor...面板堆石坝应力应变数值计算_黄小华.pdf

设计施工水利规划与设计 年第 期：.基于 模型面板堆石坝应力应变数值计算黄小华，吴乐海（福建水利电力职业技术学院，福建 永安）摘要：作为面板堆石坝受力主体的堆石料具有非线性、剪胀剪缩性等特点，选择一个合理的堆石料本构模型对面板堆石坝数值计算至关重要。模型理论以连续介质力学为基础，在模型中加入体应变项，很好的反应了堆石料的剪胀剪缩特性。文章分别采用 模型和沈珠江模型对福建省仙游抽水面板堆石坝下库进行数值计算，两者模型在反应坝体以及面板的应力应变规律上相似，数值偏差不大，证明了 模型在面板堆石坝实际工程中的应用是可行的。关键词：模型；面板堆石坝；有限元中图分类号：.文献标识码：文章编号：（）收稿日期：基金项目：福建省教育厅教科研课题（）作者简介：黄小华（年），男，讲师。：.通讯作者：吴乐海（年），男，讲师。：.在全球能源紧缺，碳中和政策的背景下，新能源的开发已变得越来越受关注。利用流量和落差发电的水电能源能实现零碳排放，符合当下的能源要求。作为当前坝工界主要坝型之一的面板堆石坝因其具有施工方便，投资相对较低，在高度提升情况下占地较小等优点，在水利工程中得到广泛应用。由于面板堆石坝防渗的主体是混凝土面板，受力主体为堆石料，若在设计阶段对堆石料变形预测不当，后期可能会导致面板脱空、拉裂等情况的发生，从而造成严重渗漏，还有可能造成溃坝的风险，因此对堆石料应力应变的研究显得尤为重要。由于堆石料具有剪胀剪缩、非线性等性质，因此在堆石料的应力应变分析中选择一个合理的本构模型至关重要。目前，弹塑性本构模型和非线性弹性模型是国内外坝工界对混凝土面板堆石坝力学分析的 大主要本构模型。但是，这 类本构模型都不能很好地反映堆石料剪胀剪缩、非线性等力学特征，在考虑颗粒破损方面也较为欠缺，且忽略了屈服面与应力路径的关系。近年来，国内外有学者将亚塑性理论应用到无黏性土的本构模型研究中，提成各种亚塑性本构模型，该类模型能弥补传统本构模型的不足，能更好地反映堆石料特殊的力学特征。本文简单介绍 亚塑性本构模型及其有限元计算方程。并将该亚塑性本构模型应用于福建仙游抽水蓄能电站下库混凝土面板堆石坝力学数值计算中，同时将该模型计算的数值成果与沈珠江双屈服面弹塑性本构模型计算的数值成果进行对比。旨在验证该本构模型应用于面板堆石坝应力应变分析的可行性。亚塑性本构模型亚塑性本构模型理论是在寻求颗粒骨架之间的级配性质的基础上，以连续介质力学理论，通过建立应力率与应变率的张量函数，在函数中加入其他本构理论，如非线性理论、考虑颗粒破损、湿化效应、剪胀性等情况，对亚塑性本构模型进行不断完善，从而得到适用于堆石料的本构模型。对砂砾料进行研究，在前人研究压缩性模型的基础上引入 临界准则，加入了临界状态函数“”，使得模型在模拟颗粒破坏状态下的应力应变关系更加真实。模型表达式见式（）所示。年第 期水利规划与设计设计施工?（）?（?）（）?（）式中，?应力率张量；规范化的应力率张量 （）；?变形率张量；?欧式范数?（?）；临 界 状 态 函 数，轴对称状态下为；向密因子，；标 量 因 子，与颗粒骨架的平均应力 和 孔 隙 比 有 关；临 界 参 数，（），由临界摩擦角 确定。通过推导，可得 亚塑性模型参数为，若直接将模型用于堆石料应力应变曲线拟合，主堆石料和次堆石料的拟合曲线无法很好反应堆石料本身剪胀的特性。刘国明在 模型的基础上加入体应变项（?）（），使改进后的 模型能够更好反应堆石料的减胀特性。因此，为改善 亚塑性模型反应堆石料剪胀性不佳的情况，加入体应变项（?）（），改进后的 亚塑性模型见式（）：?（）?（?）（?）（）（）?（）模型中参数 和 的求取，吴长彬在常规三轴压缩试验下得到与平均压力 有关，为简化其关系，引入参数、，如式（）（）所示：（）（）（）（）简化后的参数共 个：、，使得参数获取更加容易，根据福建仙游抽水蓄能电站的筑坝材料试验报告数据计算得到主堆石料和次堆石料的模型参数见表：亚塑性材料矩证模量的确定有限元的理论是将整体进行分割成单元，将单元利用虚功原理外力功与应力功平衡，列出平衡方程，再对单元进行叠加求得整体平衡方程。单元的平衡方程如下：（）（）式中，节点虚位移；虚应变列阵，可表示为 ；单元应力列阵 。将式（）进行转换可得单元劲度矩阵，再叠加单元劲度矩阵得到整体劲度矩阵，得到整体平衡方程式（）：（）（）式中，整体节点位移列阵；整体结点荷载列阵；整体劲度矩阵；材料矩阵模量。因此上式有限元平衡方程的求解取决于，由率型本构方程：（）：（）（）根据 张 量 计 算 方 法，：，带入上式（）简化得到：（）（）（）因此，可得到矩阵的表达式为：（）（）（）对式（）进行展开：（）表 模型参数表材料名称 （）主堆石料.次堆石料.设计施工水利规划与设计 年第 期 提出式（）中的，结合式（），可得到矩阵模量分量：，（）式中，下标 ，；，从而可得整体的劲度矩阵求解出有限元平衡方程（）。亚塑性本构模型在工程中的应用.工程概况及计算方法为了验证改进的 亚塑性本构模型应用于面板堆石坝应力应变分析的可行性，本次选用的工程案例为福建仙游抽水蓄能电站下库混凝土面板 堆 石 坝，该 大 坝 的 基 本 参 数 如 下：坝 高.，坝顶宽度.，坝长，前坡坡比：.，后坡坡比：.，混凝土面板的材料特性参数取 .、.、.。大坝的断面图如图 所示，对大坝有限元网格剖分如图 所示。对式（）式非线性有限元方程的求解，采用中点增量法，对仙游面板堆石坝的下库大坝进行分 级加载，施加本级荷载的一半力于结构，得到本级的平均值，使得结构的精度得到提高，中点增量法示意图如图 所示。.计算结果分析.坝体应力应变结果分析在对福建仙游抽水蓄能电站下库面板堆石坝的三维有限元分析中，堆石料的本构模型分别采用改进的 亚塑性本构模型和沈珠江弹塑性模型，计算工况为竣工期和正常蓄水期。种本构模型计算的坝体应力及变形结果见表，各种工况下的水平位移、垂直位移及大小主应力等值线图如图 所示。其中，水平方向的位移以向下游为正，垂直方向的位移以向下为正，应图 大坝断面图（单位：）图 大坝有限元网格图图 中点增量法 年第 期水利规划与设计设计施工表 不同模型坝体应力应变结果对比坝体应力变形沈珠江模型 模型竣工期 蓄水期竣工期蓄水期变形最大水平位移.最小水平位移.最大垂直位移.应力第一主应力最大值.第三主应力最大值.力以压为正。由表 可以看出，竣工期 模型计算出的最大水平位移为.，最大垂直位移为.，沈珠江模型最大水平位移.，最大垂直位移.。竣工期 模型计算出的最大水平位移比沈珠江模型更大，增幅达到，种本构模型计算出的最大垂直位移较为接近，约占坝高的.；蓄水期模型最大水图 大坝河床断面水平位移等值线分布图图 大坝河床断面垂直位移等值线分布图图 大坝河床断面第一主应力等值线分布图设计施工水利规划与设计 年第 期图 大坝河床断面第三主应力等值线分布图平位移.，最大垂直位移.，沈珠江模型最大水平位移.，最大垂直位移.。蓄水期 模型计算出的最大水平位移比沈珠江模型增大，最大垂直位移的比较中，模型也比沈珠江模型更大.。从计算结果上看，种模型位移的计算结果偏差不大，在水压力作用下，下游向的水平位移比上游向的水平位移来的更大。从位移等值线分布云图分析，两者模型分布特征基本相似，在堆石的自重作用下，坝体垂直位移最大值出现在 坝高处，并向四周逐渐减小。水平位移大致成对称性，呈下凸上凹形状，略有不同在于蓄水期 模型坝体上游侧中上部位指向下游侧，沈珠江模型指向上游侧。由表 可知，竣工期 模型第一主应力最大值为.，第三主应力最大值为.，沈 珠 江 模 型 相 应 的 应 力 为.，.，从计算结果可以看出，在竣工期，种本构模型第一主应力计算数值非常接近，第三主应力二者差值为.，从应力的等值线分布图来看，两者模型的应力分布规律相似，大小主应力均出现在大坝的底部。蓄水期 模型第一主应力最大值为.，第三主应力最大值为.，对比沈珠江模型为.、.，种模型计算结果非常接近。从图 还可以看出，种模型的第三主应力最大值较竣工期均有往上游移动的趋势。.面板应力应变结果分析蓄水期 种本构模型计算的面板应力变形计算成果见表、如图 所示，规定应力以压为正，变形以下游向为正。表 不同模型蓄水期面板应力应变结果对比应力变形沈珠江模型 模型变形最大挠度.面板最大坝轴线向位移.面板最小坝轴线向位移.应力面板最大顺坡向向应力.面板最小顺坡向应力.面板最大坝轴线向应力.面板最小坝轴线向应力.如图 蓄水期面板位移云图所示，模型和沈珠江模型挠度位移和轴线向位移分布规律相同，最大挠度部分均出现在坝体 区域，并从坝体两侧向中间逐渐变大。从云图可以看出，混凝土面板最大扰度的位置与坝体堆石料最大垂直位置非常接近。从表 的变形计算结果可知，沈珠江模型最大挠度为.，轴线向位移最大为.，最小为.。模型最大挠度为.，轴线向位移最大为.，最小为.。模型计算的最大扰度比沈珠江模型增加，相比下 模型的计算结果偏于安全。如图 蓄水期面板应力云图所示，种模型的顺坡向和轴线向的应力分布规律相似，都能够体现顺坡向应力中部受压，顶部和底部受拉的状态，轴线向应力中部受压，两边受拉状态，应力状态分布合理。沈珠江模型最大顺坡向应力为.，最小为 .，轴线向最大应力为.，最小为 .。模型最大顺坡向应力为.，最小为 .，轴线向最大应力为.，最小为.。从计算结果可以 年第 期水利规划与设计设计施工图 蓄水期面板位移云图图 蓄水期面板应力云图看出，模型计算的面板应力分布比沈珠江模型更为合理。沈珠江模型计算的混凝土面板底部顺坡向拉应力存在应力过大的问题，与实际情况不符。结论（）模型的建立是以连续介质力学为基础，沈珠江模型是以弹塑性模型为基础，两者理论不同，但两模型通过对仙游面板堆石坝的坝体及面板应力应变计算对比分析，两者结果相差不大，从实际验证了改进后的 模型可以应用于实际工程的应力应变分析。（）在混凝土面板应力应变计算方面，相比于沈珠江弹塑性模型，模型计算参数简单，容易确定，使用简便，且能够反映堆石料剪胀性，应力应变非线性等力学特征。但是该模型也存设计施工水利规划与设计 年第 期在一些不足，如不能考虑堆石料流变特征等，有待进一步研究。参考文献 黄小华 在流变效应影响下面板堆石坝应力变形分析 浙江水利水电学院学报，（）：钮新强 高面板堆石坝安全与思考 水力发电学报，（）：袁磊，向佐年，杨杰，等 阿尔塔什面板堆石坝混凝土面板设计关键技术 水利规划与设计，（）：杨智，陈立宝 深厚覆盖层基础建面板堆石坝的应力变形分析 水利技术监督，（）：王华俊，朱智勇，杨兰强 堆石料基本性质及其本构模型研究进展 科技通报，（）：刘京茂 堆石料和接触面弹塑性本构模型及其在面板堆石坝中的应用研究 大连理工大学，王洪波，张学增，熊保林 亚塑性理论简介 岩土力学，（）：明华军，孙开畅，徐小峰，等 考虑颗粒破碎对特征孔隙比影响 的堆石体亚塑性本构模型 岩土力学，（）：岑威钧，考虑湿化效应的堆石料 亚塑性模型应用研究 岩土力学，（）：刘国明，陈泽钦，吴乐海 堆石料 亚塑性本构模型改进及参数确定方法 岩土力学，（）：，（）：吴乐海 亚塑性模型及其在面板堆石坝应力变形分析中的应用研究 福州大学，刘国明，陈泽钦，吴乐海 堆石料 亚塑性本构模型改进及参数确定方法 岩土力学，（）：吴长彬，张岩，许小东 应用于堆石坝的弹塑性模型与亚塑性模型的理论对比分析 灾害与防治工程，（）：南京水利科学研究院 福建仙游抽水蓄能电站筑坝材料试验报告 王智，吴洋，朱先发，等 基于亚塑性本构模型的砂土剪切特性分析 隧道与地下工程灾害防治，（）：（上接第 页）东 东厢分洪路径方案的运用成本更低，分洪造成的经济损失为.亿.亿元。因此，意东分洪推荐意东 东厢分洪路径方案。东厢分洪方案最大分洪流量 ，分洪造成的经济损失为.亿.亿元。总体上，东厢分洪方案的运用效果最好，且运用成本最低，比选为推荐方案。（）当韩江遭遇 年一遇洪水，南北堤一旦溃决，造成的经济损失约为 亿 亿元（按 年地区经济发展水平计算）。东厢分洪方案可成为防御韩江流域超标准洪