[发明专利]一种模拟砂砾石坝基管涌侵蚀对大坝应力变形影响的方法

说明书

技术领域

本发明属于水利水电工程堤坝的设计、研究及控制的技术领域，具体地涉及一种模拟砂砾石坝基管涌侵蚀对大坝应力变形影响的方法。

背景技术

部分江河大堤和一些水库大坝都不可避免地坐落于在渗流作用下内部不稳定的砂砾石覆盖层上，这些堤坝在高水头渗流作用下可能发生内部管涌侵蚀，严重时威胁堤坝的安全。砂砾石土的管涌侵蚀，会改变土体的渗透性、强度和应力变形关系，导致坝基及坝体渗流场的变化、导致大坝发生变形和应力调整，危害坝基的渗透稳定性和坝基的抗滑稳定性，危害上部及坝基内部结构的正常工作，最终危害整个堤坝的安全运行。

砂砾石土的管涌侵蚀的类型包括：管道型(流土型)管涌和管涌型管涌。管涌型管涌是指粗颗粒构成骨架而细颗粒在粗颗粒之间的孔隙中发生移动侵蚀。管涌型管涌可以在堤坝基础内部首先发生，并溯源发展。管涌型管涌侵蚀发生对应力变形有两个方面的影响：(1)侵蚀改变土体的渗透性，从而改变坝基的渗流场，而导致坝基及其上部结构产生应力变形；(2)侵蚀改变土体的强度和变形模量，土体在应力状态不变时将产生应变，从而引起坝基的应力动态调整和迁移。为了保证堤坝的安全运行，亟需一种模拟砂砾石坝基管涌侵蚀对大坝应力变形影响的方法，以便定量描述砂砾石坝基管涌侵蚀对大坝应力变形造成的影响。

发明内容

本发明的技术解决问题是：克服现有技术的不足，提供一种模拟砂砾石坝基管涌侵蚀对大坝应力变形影响的方法，其能够定量描述砂砾石坝基管涌侵蚀对大坝应力变形造成的影响，用于分析和控制管涌型管涌对大坝应力变形造成的影响，从而保证整个堤坝安全运行。

本发明的技术解决方案是：这种模拟砂砾石坝基管涌侵蚀对大坝应力变形影响的方法，该方法包括以下步骤：

(1)仿真计算坝体填筑和蓄水过程中坝体的渗流和应力变形，当坝基中砂砾石土中的渗透坡降超过临界管涌坡降时，判断管涌侵蚀开始发生；

(2)进行管涌侵蚀计算和大坝应力变形计算，并对二者进行耦合，获得管涌侵蚀对大坝应力变形影响的结果。

本发明通过仿真计算坝体填筑和蓄水过程中坝体的渗流和应力变形，进行管涌侵蚀计算和大坝应力变形计算，并对二者进行耦合，获得管涌侵蚀对大坝应力变形影响的结果，所以能够定量描述砂砾石坝基管涌侵蚀对大坝应力变形造成的影响，用于分析和控制管涌型管涌对大坝应力变形造成的影响，从而保证整个堤坝安全运行。

附图说明

图1a示出了根据本发明一个优选实施例的典型剖面在正常高水位管涌侵蚀100小时后由于管涌渗流场变化引起的水平位移增量；图1b示出了根据本发明一个优选实施例的典型剖面在正常高水位管涌侵蚀100小时后由于管涌渗流场变化引起的垂直位移增量。

图2a示出了图1a的管涌渗流场变化后大坝的水平正应力增量；图2b示出了图1b的管涌渗流场变化后大坝的垂直正应力增量。

图3a示出了图1a的管涌土体模量衰化引起的水平位移增量；图3b示出了图1b的管涌土体模量衰化引起的垂直位移增量。

具体实施方式

这种模拟砂砾石坝基管涌侵蚀对大坝应力变形影响的方法，该方法包括以下步骤：

(1)仿真计算坝体填筑和蓄水过程中坝体的渗流和应力变形，当坝基中砂砾石土中的渗透坡降超过临界管涌坡降时，判断管涌侵蚀开始发生；

(2)进行管涌侵蚀计算和大坝应力变形计算，并对二者进行耦合，获得管涌侵蚀对大坝应力变形影响的结果。

本发明通过仿真计算坝体填筑和蓄水过程中坝体的应力变形，进行管涌侵蚀计算和大坝应力变形计算，并对二者进行耦合，获得管涌侵蚀对大坝应力变形影响的结果，所以能够定量描述砂砾石坝基管涌侵蚀对大坝应力变形造成的影响，用于分析和控制管涌型管涌对大坝应力变形造成的影响，从而保证整个堤坝安全运行。

另外，所述步骤(2)在每一个计算时步中耦合管涌侵蚀计算和大坝应力变形计算。

或者，所述步骤(2)包括以下分步骤：

(2.1)以较小的时步进行管涌侵蚀计算，以较大的时步进行应力变形计算，管涌侵蚀计算不耦合应力变形计算结果，而应力变形计算耦合管涌侵蚀计算结果；

(2.2)以较大的时步时的渗流场、管涌侵蚀计算结果为基础，引入描述土体管涌侵蚀过程中应力应变关系的本构模型，分两部分计算坝体的应力变形增量：(a)采用未侵蚀时土体的变形参数计算渗流场变化引起的应力和位移增量；(b)计算由于侵蚀引起变形参数变化导致的应力调整和迁移产生的应力和位移增量。

另外，所述部分(a)中渗流场变化引起的应力和位移增量，基于管涌后的渗流场结果，采用渗流应力耦合算法。