[发明专利]一种砂砾石土在管涌侵蚀中的本构关系描述方法

说明书

技术领域

本发明属于水利水电工程堤坝的设计、研究及控制的技术领域，具体地涉及一种砂砾石土在管涌侵蚀中的本构关系描述方法，主要用于堤坝在渗流作用下管涌侵蚀动态模拟，及其对堤坝渗流和应力变形及安全性的影响分析。

背景技术

部分江河大堤和一些水库大坝，都不可避免地坐落于渗流作用下内部不稳定的砂砾石覆盖层上，这些堤坝在高水头渗流作用下可能发生内部管涌侵蚀，进而威胁堤坝的安全。砂砾石土的管涌侵蚀，会改变土体的渗透性、强度和应力变形关系，导致坝基及坝体发生变形和应力调整。因此，亟需能够定量地分析管涌侵蚀对坝基及坝体应力变形和安全性的影响程度，而这就需要一个可以描述土体在管涌侵蚀过程中及侵蚀完成后应力应变关系的方法，用于模拟侵蚀造成的坝基和坝体应力变形。

发明内容

本发明的技术解决问题是：克服现有技术的缺陷，提供一种砂砾石土在管涌侵蚀中的本构关系描述方法，其能够定量地描述土体在管涌侵蚀过程中及侵蚀完成后的应力-应变关系，用于堤坝在渗流作用下管涌侵蚀动态模拟，及其对堤坝渗流和应力变形及安全性的影响分析。

本发明的技术解决方案是：这种砂砾石土在管涌侵蚀中的本构关系描述方法，该方法以颗粒流失量为自变量，以本构模型的强度和变形参数为因变量，通过建立模型参数与这个自变量的关系来模拟对砂砾石土不同大小的颗粒流失量时的应力变形。

由于在渗流作用下，内部不稳定的砂砾石土，由粗颗粒构成土体的骨架，因此，在渗流作用下部分细颗粒流失后的应力应变特性与未侵蚀时基本相同，可以用相同的模型来模拟，侵蚀改变的是模型参数，这些参数的改变可以用颗粒流失量作为自变量，通过建立起模型参数与这个自变量的关系实现对砂砾石土不同大小的颗粒流失量时的应力变形关系，从而能够定量地描述土体在管涌侵蚀过程中及侵蚀完成后应力应变关系，用于堤坝在渗流作用下管涌侵蚀动态模拟，及其对堤坝渗流和应力变形及安全性的影响分析。

附图说明

图1a示出了覆盖层邓肯E-B模型的体变参数K_b与侵蚀率的关系，图1b示出了体变参数m与侵蚀率的关系。

图2a示出了覆盖层邓肯E-B模型的剪切变形参数K与侵蚀率的关系，图2b示出了剪切变形参数n与侵蚀率的关系。

图3示出了覆盖层邓肯E-B模型的应力水平对剪切模量的影响参数R_f与侵蚀率的关系。

具体实施方式

这种砂砾石土在管涌侵蚀中的本构关系描述方法，该方法以颗粒流失量为自变量，以本构模型的强度和变形参数为因变量，通过建立模型参数与这个自变量的关系来模拟对砂砾石土不同大小的颗粒流失量时的应力变形。

由于在渗流作用下，内部不稳定的砂砾石土，由粗颗粒构成土体的骨架，因此，在渗流作用下部分细颗粒流失后的应力应变特性与未侵蚀时基本相同，可以用相同的模型来模拟，侵蚀改变的是模型参数，这些参数的改变可以用颗粒流失量作为自变量，通过建立起模型参数与这个自变量的关系实现对砂砾石土不同大小的颗粒流失量时的应力变形关系，从而能够定量地描述土体在管涌侵蚀过程中及侵蚀完成后应力应变关系，用于堤坝在渗流作用下管涌侵蚀动态模拟，及其对堤坝渗流和应力变形及安全性的影响分析。

另外，所述模型参数通过侵蚀前土样和侵蚀后细颗粒减少的土样的常规三轴试验和侧限压缩试验来获得。

另外，当采用邓肯E-B模型时，参数与自变量的关系的描述方式为：

C＝a(1-β)^b

其中a、b为参数，β为砂砾石土颗粒侵蚀率，C是模型中的参数的统称，不同的模型参数，以上公式中的a、b参数不相同。

另外，管涌侵蚀流失后的土样，在未侵蚀试样配料表的基础上，扣减掉细颗粒的流失重量制样获得。

下面以某砂砾石土管涌过程中的应力应变关系的情况，用邓肯E-B模型为例来举例说明。

邓肯模型的变形参数有5个，包括体变参数k_b和m，初始剪切模量参数k、n，应力水平影响参数R_f。

如图1a、1b所示，邓肯E-B模型的两个体变参数k_b和m用如下公式描述：