Belgorod, Russian Federation
Russian Federation
CSCSTI 70.17
Russian Classification of Professions by Education 35.06.01
Russian Library and Bibliographic Classification 4
Russian Trade and Bibliographic Classification 5640
Russian Trade and Bibliographic Classification 5607
BISAC TEC003050 Agriculture / Irrigation
The relevance of the study is due to the need to ensure reliable operation of class IV hazard earth dams under conditions of changing groundwater levels, seismic and climatic loads, as well as compliance with the requirements of the Water Strategy of Russia until 2035. The object of the study is a zoned earth dam with a reinforced concrete core located on a water-saturated inclined foundation in the Kumensky district of the Kirov region. The model includes physical and mechanical properties of soils obtained from field surveys. For the analysis, numerical modeling by the finite element method in the MidasGTSNX software package and analytical calculations performed according to the methods of G.N. The results showed satisfactory agreement between the numerical and analytical depression curves (the difference is ≤0.06 m), while in the downstream pool there is a shift in the depression curve to the inclined foundation zone, indicating insufficient waterproofing of the core and the need for its strengthening. The analysis of saturation isofields revealed potential instability zones at the transitions of soil layers. The core of the hydraulic structure protects the downstream slope from excess water saturation, but its design should be improved (grooving, drainage prism) to increase the strength factor and stability of the slope. The obtained research results can be used in the design, reconstruction or major repairs of earth dams
inclined foundation, filtration, saturation, depression curve, earth dam, water saturation, hydraulic structure
Введение. Конструктивные исполнения грунтовых плотин IV класса опасности требуют тщательного проектирования и анализа при их эксплуатации, так как они подвержены различным внешним воздействиям, включая меняющийся уровень грунтовых вод, сейсмическую активность, атмосферные осадки и иные динамические нагрузки (от автомобильного или иного транспорта) [1]. Поддержание таких объектов гидротехнического назначения в рабочем состоянии является важнейшей задачей водной стратегии России на период до 2035 года [2].
В последние десятилетия численное моделирование стало основным инструментом для анализа поведения гидротехнических сооружений (ГТС) подобного рода с большим сроком эксплуатации (25 лет и более). Численные модели объектов гидротехнического назначения в мелиорации позволяет инженерам прогнозировать их работоспособность во времени, тем самым обеспечивать их дальнейшую безопасность на всём жизненном цикле [3].
Применение метода конечных элементов (МКЭ) [4] позволяет проводить детализированный анализ различных показателей устойчивости и прочности грунтовой плотины, что является критически важным для предотвращения аварий и разрушений гидротехнических сооружений [5].
В составе государственных и частных мелиоративных систем находятся различные по конструкции и назначению грунтовые плотины. Они обеспечивают безопасное использование водных ресурсов в водохранилищах для мелиорации и рыборазведения [6].
Актуальность настоящего исследования заключается в проверке адекватности результатов численного моделирования фильтрации эксплуатируемой более 25 лет грунтовой плотины с ядром на наклонном водопроницаемом основании для выявления возможных в дальнейшем прогрессирующих обрушений ее откосов, а также получения информации об их реальном состоянии на момент проведения работ по реконструкции или капитальному ремонту.
Цель исследования заключается в проведении сравнительного анализа результатов численного моделирования и аналитических расчетов фильтрации зонированной грунтовой плотины с ядром, расположенной на водопроницаемом (водонасыщенном) наклонном основании. Основные задачи исследования включают:
- оценку фильтрационной способности конструкции (по положению кривой депрессии) при различных уровнях воды в водохранилище, как численно, так и аналитически;
- численное определение сатурации в теле плотины при различном уровне воды в водохранилище.
Методы и объект исследования. Объектом исследования является грунтовая плотина с непроницаемым (железобетонным) ядром на водонасыщенном (водопроницаемом) основании, входящая в состав мелиоративной системы, расположенной в Куменском районе Кировской области. Расчетная модель гидротехнического сооружения в конечных элементах представлена на рис. 1. Для достижения поставленных целей использовались методы численного моделирования с применением программного обеспечения Midas GTS NX, которые в последующем верифицировались с аналитическими решениями задачи фильтрации по известным методиками ее расчета.
Рис. 1 - Расчетная модель гидротехнического сооружения без подтопления нижнего бьефа в конечных элементах
Для создания модели плотины были использованы данные о геометрии, свойствах грунта и условиях эксплуатации. Модель состоит из различных слоев грунта с различными физико-механическими свойствами (табл. 1), которые были взяты из геолого-геодезических изысканий перед началом реконструкции объекта исследования.
Табл. 1 - Физико-механические свойства грунтов исследуемой плотины
|
Параметры |
Величина |
|||
|
Наклонное основание |
Откосы |
Ядро |
Фильтр |
|
|
1 |
2 |
3 |
4 |
5 |
|
Удельный вес грунта при естественной влажности, кН/м3 |
17 |
16 |
22 |
18 |
|
Удельный вес грунта в водонасыщенном состоянии, кН/м3 |
21 |
20 |
25 |
22 |
|
Начальный коэффициент пористости |
0,5 |
0,5 |
0,5 |
0,5 |
|
Эффективный модуль Юнга, кН/м2 |
50·103 |
20·103 |
1,5·103 |
50·103 |
|
Коэффициент Пуассона |
0,3 |
0,33 |
0,35 |
0,3 |
|
Модуль сдвига, кН/м2 |
19,23·103 |
7519 |
555.600 |
19,23·103 |
|
Касательный одометрический модуль деформации, кН/м2 |
6,73·103 |
29,63·103 |
2407 |
6,73·103 |
|
Параметр референсного давления, кН/м2 |
1 |
5 |
5 |
1 |
|
Угол внутреннего трения, рад |
35 |
31 |
0 |
35 |
|
Угол дилатансии, рад |
5 |
1 |
0 |
5 |
|
Скорость прохождения поперечных волн, м/с |
105,3 |
67,9 |
18,46 |
106,2 |
|
Скорость прохождения продольных волн, м/с |
197,1 |
134,8 |
38,42 |
198,2 |
|
Коэффициент фильтрации по оси Х |
0,01 |
1 |
0,1·10–3 |
0,01 |
|
Коэффициент фильтрации по оси Y |
0,01 |
1 |
0,1·10–3 |
0,01 |
|
Глубинный параметр неравномерности распределения дилатансии, м |
10·103 |
10·103 |
10·103 |
10·103 |
|
Удельный параметр неравномерности распределения дилатансии, 1/м |
5,33·10–6 |
0,13·10–3 |
0,18·10–3 |
5,72·10–5 |
|
Реологический модуль |
1000·1012 |
1000·1012 |
1000·1012 |
1000·1012 |
В модели были заданы граничные условия, такие как уровень грунтовых вод, нагрузки от воды без учета сейсмичности района моделирования [7]. Эти параметры были выбраны на основе данных из аналитических расчетов по этому объекту и практики анализа аналогичных сооружений. Моделирование проводилось для нескольких сценариев при динамических нагрузках, возникающих в результате изменения уровня воды в водохранилище: при полном водохранилище; при уровне воды в водохранилище на отметке мертвого уровня и устоявшемся уровне воды на отм. 2,5 м (такая отметка водного уровня берется с учетом максимально возможных атмосферных осадков в период проведения работ по капитальному ремонту или реконструкции для данной местности). Временной промежуток изменения уровня воды в водохранилище относится к измерениям через 0 ч, 24 ч, 48 ч, 72 ч, 96 ч и 120 ч. При реконструкции необходимо спускать водохранилище до уровня мертвого объема (УМО, показано на рис. 1), соответственно, нормальный подпорный уровень воды (НПУ, показано на рис. 1) является началом отсчета и отметка 2,5 м – устоявшийся уровень воды (УМО).
После проведения расчетов было проанализировано положение кривой депрессии в теле плотины и сатурация её тела. Особое внимание уделялось зонам перехода между различными слоями грунта, так как именно здесь могут возникать в дальнейших численных и аналитических расчетах в реальном объекте потенциальные зоны неустойчивости.
Верификация положения кривой депрессии, полученной в результате численного моделирования, в исследуемой плотине проводилась в соответствии с методикой, изложенной в работе [8, C.376] при наклонном водоупоре по приближенной формуле Г. H. Каменского [9, С.265] и методу Н.Н. Павловского [9, C.266].
Основная часть. На рис. 2 показаны изополя степени насыщения «degree of saturation» для начальной стадии при установившемся режиме фильтрации и последующем снижении уровня воды. Максимальное значение, равное 1, наблюдается в «красной» зоне. Следовательно, в этой зоне устанавливается стационарная фильтрация для всех материалов, из которых выполнена плотина.
Рис. 2 - Изополя сатурации (водонасыщения) в зонированной грунтовой плотине:
а) – полное водохранилище (0 ч); б) – снижение уровня воды (24 ч); в) – снижение уровня воды (48 ч); г) – снижение уровня воды (72 ч); д) – снижение уровня воды (96 ч - УМО); е) – снижение уровня воды (120 ч - устоявшийся уровень воды на отм. 2,5 м)
Обратим внимание на ядро плотины. Часть ядра плотины находится в зоне сатурации с показателем равным 1. Это говорит о том, что конструкция ядра недостаточно изолирует нижний откос от избыточного влагонасыщения, следовательно, такая ситуация может привести к обрушению нижнего откоса, о чем свидетельствует «синяя» зона, расположенная на всей площади нижнего откоса.
Помимо прочего граница сатурации между максимально водонасыщенным наклонным основанием и нижним откосом устанавливает возможную зону обрушения откоса. Аналогичная ситуация может вызвать обрушение откоса в верхней части верхнего откоса при устоявшемся уровне воды на отм. 2,5 м (рис. 2 д, е), так как здесь образуются круглоцилиндрические поверхности, определяющие переход водонасыщенной зоны плотины (участок «красного» цвета) к неводонасыщенной (участок «сине-зеленого» цвета).
Кривые депрессии по телу плотины для различных стадий наполнения водохранилища показаны на рис. 3.
Рис. 3 - Кривая депрессии в исследуемой грунтовой плотине:
а) – полное водохранилище (0ч); б) – снижение уровня воды (24ч); в) – снижение уровня воды (48ч); г) – снижение уровня воды (72ч); д) – снижение уровня воды (96ч - УМО); е) – снижение уровня воды (120ч - устоявшийся уровень воды на отм. 2,5 м)
Характерной особенностью всех полученных линий кривых депрессий для различных уровней воды в водохранилище является то, что ниспадающая часть всех кривых находится в правой зоне ядра плотины (ближе к нижнему откосу), что также говорит о том, конструкция ядра (в том числе и материал ядра) плотины при капитальном ремонте или реконструкции объекта должна быть усовершенстована и больше гидроизолирована, в том числе за счет шпунтования основания ядра. Характер кривой депрессии, располагаемой под границей наклонного основания и нижнего откоса, совпадает с границей водонасыщения (сатурации), что говорит о необходимости изменения конструкции нижнего откоса и переустройства всего нижнего бьефа объекта. В нижней части нижнего откоса для данного исследуемого объекта необходимо предусмотреть дренажную призму, чтобы окончание кривой депрессии под нижним бьефом плотины находилось под границей наклонного основания и нижнего откоса [10].
Анализ численных данных с обязательной верификацией полученных результатов по параметрам фильтрации и сатурации объектов гидротехнического назначения в настоящее время является неотъемлемой частью реализации внедрения технологий BIM на всех треках строительной отрасли [11].
Результаты и их обсуждение. Проведем верификацию положения кривой депрессии по исследуемой плотине при максимальном уровне воды в водохранилище согласно методике изложенной в работах [8, C.376, 9, С.265, 9, C.266]. Получим графический результат в следующем виде – на рис. 5 показано сравнение положения кривой депрессии в исследуемой плотине при полном водохранилище с положением кривой депрессии, полученной в результате аналитического расчета. Максимальное расхождение положения кривой депрессии наблюдается в нижнем бьефе плотины и составляет по вертикальной оси 0,17 м (результаты аналитического расчета показывают, что кривая депрессии на этом участке находится в зоне наклонного основания, а не в зоне низового откоса, как это показано на результатах численного моделирования - рис. 5). На остальных участках сравнения кривых депрессии расхождение между ними по вертикальной оси не превышает 0,06 м. Кривые депрессии, полученные численными методами и аналитическим расчетом практически совпадают (расхождение их по вертикальной координате не превышает 0,01 м) на участке верхового откоса, ядра и нисходящего наклонного основания плотины. Считаем данные результаты по верификации положения кривой депрессии весьма удовлетворительными и достаточно точными (рис. 4).
Рис. 4 – Верификация положения кривой депрессии в теле исследуемой плотины: – (синяя линия) численный расчет положения кривой депрессии;
2 – (красная линия) аналитический расчет положения кривой депрессии
Подтвержденное положение кривой депрессии также говорит о том, что необходимо укреплять нижний бьеф плотины в месте ее сопряжения с наклонным основанием, как было об этом сказано при анализе рис. 3.
Таким образом, в процессе численного моделирования объекта гидротехнического назначения получены не только числовые значения рассматриваемых в исследовании параметров, но и сделаны выводы о необходимых проектно-изыскательских работах по совершенствованию конструкции плотины с целью повышения ее надежности и безопасности при дальнейшей эксплуатации.
Выводы. Проведена оценка фильтрационной способности конструкции грунтовой плотины с ядром на наклонном водопроницаемом основании при различных уровнях воды в водохранилище по полученным кривым депрессии в ходе численного моделирования объекта. Установлено, что ядро плотины защищает нижний откос от избыточного водонасыщения и обрушения, однако положение кривой депрессии целесообразно сместить в сторону верхового откоса для повышения коэффициента прочности и устойчивости низового откоса плотины.
Результаты верификации параметров фильтрации и сатурации исследуемой плотины считаем удовлетворительными. Полученные данные при численном моделировании объекта могут быть использованы при проектных работах в рамках реконструкции или капитального ремонта объекта. Полученные в результате численного моделирования параметры плотины позволяют сделать вывод о необходимых проектных решениях, повышающих надежность и безопасность гидротехнического сооружения с большим сроком эксплуатации.
1. Simonovich O. S., Snezhko V. L., Kozlov D. V. Periodichnost' remonta nizkonapornyh gruntovyh plotin dlya sohraneniya ih urovnya bezopasnosti // Prirodoobustroystvo. 2020. № 1. S. 59-64.
2. Vodnaya strategiya Rossiyskoy Federacii na period do 2035 goda [Elektronnyy resurs]. M., 2025. URL: https://raww.ru/assets/modckeditor/default/0/novaya-vodnaya-strategiya-do-2035.pdf
3. Kachaev A. E., Turapin S. S. Obosnovanie neobhodimosti razrabotki kompleksnyh raschetnyh modeley gruntovyh plotin meliorativnyh sistem // Nauka i mir. 2024. № 3. S. 1-5.
4. Gaydzhurov P. P., Al'-Hadzh M. A., Al'-Dzhabobi S. F. Metod konechnyh elementov v stroitel'stve // Colloquium-Journal. 2019. № 27-1(51). S. 23-27.
5. Kachaev A. E., Turapin S. S. Osobennosti rekonstrukcii zemlyanyh plotin meliorativnyh sistem // Nauka i mir. 2024. № 3. S. 6-10.
6. Shkura V. N., Shevchenko A. V. Kompleks rybovodnyh gidrotehnicheskih sooruzheniy dlya rybohozyaystvennoy melioracii evtrofnyh vodohranilisch // Vodnoe hozyaystvo Rossii: problemy, tehnologii, upravlenie. 2022. № 5. S. 95-106.
7. Ol'garenko G. V., Turapin S. S., Kashtanov V. V., Savushkin V.V., Gzhibovskiy S.A., Abramov V.V. Bezopasnaya ekspluataciya gidrotehnicheskih sooruzheniy meliorativnogo kompleksa. Kolomna: IP Lavrenov A.V., 2018. 228 s.
8. Rasskazov L.N., Orehov V.G., Aniskin N.A., Malahanov V.V., Bestuzheva A.S., Sainov M.P., Soldatov P.V., Tolstikov V.V. Gidrotehnicheskie sooruzheniya (rechnye): Uchebnik dlya vysshih uchebnyh zavedeniy. Chast' 1. 2-e izdanie, ispravlennoe i dopolnennoe. – Moskva: Izdatel'stvo Associacii stroitel'nyh vuzov, 2011. – 584 s.
9. Nesterov M. V. Gidrotehnicheskie sooruzheniya: uchebnik. — 2-e izd., ispr. i dop. — Moskva: INFRA-M, 2022. — 601 s.
10. Osnovaniya i fundamenty zdaniy i sooruzheniy. Fundamenty pri vibrodinamicheskih vozdeystviyah. Pravila proektirovaniya: TKP 45-5.01-264-2012. Vved. 28.05.12. – Minsk: RUP «Stroytehnorm», 2012. – 114 s.
11. Sidorenko D. A., Kachaev A. E. BIM-tehnologii v stroitel'stve: chto budet dal'she? // Novye tehnologii v uchebnom processe i proizvodstve : Materialy XXI Mezhdunarodnoy nauchno-tehnicheskoy konferencii, posvyaschennoy 35-letiyu poleta orbital'nogo korablya-raketoplana mnogorazovoy transportnoy kosmicheskoy sistemy "Buran", Ryazan', 12–14 aprelya 2023 goda / Pod redakciey A.N. Parshina. – Ryazan': Ryazanskiy institut (filial) federal'nogo gosudarstvennogo avtonomnogo obrazovatel'nogo uchrezhdeniya vysshego obrazovaniya "Moskovskiy politehnicheskiy universitet" , 2023. – S. 490-492.
