Выбор математических моделей грунтов в статических и сейсмических расчетах грунтовых плотин
- Авторы: Ляпичев Ю.П.1,2
-
Учреждения:
- АО «Институт “Гидропроект”»
- Международная комиссия по большим плотинам (ICOLD)
- Выпуск: Том 16, № 4 (2020)
- Страницы: 261-270
- Раздел: Численные методы расчета конструкций
- URL: https://journals.rudn.ru/structural-mechanics/article/view/24468
- DOI: https://doi.org/10.22363/1815-5235-2020-16-4-261-270
Цитировать
Полный текст
Аннотация
Актуальность. Численные статические и сейсмические (динамические) расчеты напряженно-деформированного состояния (НДС) грунтовых плотин обычно выполняются с применением различных вычислительных программ, в которых используются разные, часто сложные математические модели грунтов. Однако правильный выбор этих моделей остается за пользователем программы, обычно не имеющим достаточного опыта, поэтому результаты расчетов часто бывают непонятными и ошибочными. Таким образом, цели настоящего исследования - разработка рекомендаций по выбору достоверных математических моделей грунтов в численных расчетах грунтовых плотин и их сравнение с наиболее распространенными моделями в современной практике применения. Методы. Глубокий сравнительный анализ многих моделей грунтов проводился по результатам их применения в расчетах грунтовых плотинах во время работы в Комитете СИГБ (ICOLD) по численным расчетам и проектированию плотин. Результаты. На основе оценки достоверности моделей грунтов получены рекомендации по выбору и применению моделей в численных расчетах НДС грунтовых плотин при статических и сейсмических воздействиях; определены взаимодействия между результатами расчетов и данными мониторинга поведения плотин.
Полный текст
Введение[5] Будучи членом комитета СИГБ (ICOLD) в 1993-2013 гг. по численным расчетам и проектированию плотин, нам довелось отвечал за подготовку больших разделов двух Бюллетеней СИГБ, посвященных грунтовым плотинам и математическим моделям грунтов [1-2]. В настоящей статье рассматриваются математические модели грунтов, применяемые в статических и сейсмических (динамических) численных расчетах грунтовых плотин. Динамический расчет означает, что сейсмическое воздействие используется в виде записанных или синтетических акселерограмм землетрясений. Любая численная модель сооружения представляет собой упрощенное представление реальной физической системы, отклонениями от реальности которой являются: 1) упрощение геометрии (конечные границы, упущенные детали); 2) упрощение свойств материалов в принятых определяющих законах; 3) упрощение пространства изменений поведения материалов: например, в модели изменения в пространстве перемещений необходимо следовать функциям формы (со степенью 1 или 2), которые по сути лишь приближения, часто слишком грубые; 4) схематизация распределения нагрузок в пространстве и во времени. Таким образом, результаты численных расчетов, которые обеспечивают применение численной модели, могут быть достаточно далекими от того, что происходит в реальности, в зависимости от принятых упрощений в расчетах. Сама численная модель не дает оценки собственного уровня неточности по отношению к реальной моделируемой системе. 1. Сравнение численных расчетов НДС различных грунтовых плотин при сейсмических и статических воздействиях Выбор методов численных расчетов при сейсмических (динамических) воздействиях При сейсмическом воздействии для определения реакции грунтовых плотин используются более сложные методы, чем при статическом. Реакция плотины зависит от нагрузок, которые она получает от основания и водохранилища (верхнего бьефа - ВБ), и изменения механических характеристик материалов, вызванных сейсмическим воздействием. Простые методы, такие как линейный эквивалентный метод Сида (ЛЭМС), модифицируют модуль сдвига и коэффициент затухания колебаний как функцию искривления (деформации сдвига). ЛЭМС прямо не оценивает сейсмически возникающее давление в поровой воде грунта, хотя оно играет очень важную роль, когда грунты не плотные и могут уплотняться при землетрясении. Современные динамические методы учитывают поведение поровой воды при и после землетрясения. Они используют реологические законы, описывающие дилатансию (разуплотнение при сдвиге) грунта, эффекты его размягчения, деформационное упрочнение и др. В простых динамических задачах необратимые перемещения грунтовой плотины рассчитываются по методу Ньюмарка [3]. Это включает определение предельного горизонтального ускорения колебаний, которое соответствует коэффициенту запаса устойчивости плотины 1,0 для заданной круговой или многоугольной поверхности скольжения. Согласно СП39.13330.2012 коэффициент запаса устойчивости грунтовой плотины превышает 1,0 с учетом класса плотины, коэффициента условий работы и сочетания нагрузок. Затем, используя акселерограмму землетрясения в центре тяжести неустойчивой зоны плотины, дважды интегрируется предельное горизонтальное ускорение для получения горизонтального необратимого перемещения. В сложных динамических задачах метод Ньюмарка может использоваться для обеспечения неразрывности этих задач со стандартной инженерной практикой. Методы расчета устойчивости грунтовых плотин являются простым способом использования исходных механических параметров грунта для определения коэффициента запаса устойчивости плотины. Однако опытный инженер по механике грунтов получит максимальную отдачу от расчетов, используя сложные численные модели, которые будут наиболее точными инструментами. Выбор типов моделей грунтов и методов их испытаний Численные расчеты, проведенные с помощью линейной определяющей модели грунта, являются неполными, если не учтены сейсмические нагрузки низкой интенсивности. Поэтому следует применять связанные или несвязанные расчеты с нелинейными определяющими моделями с использованием пошагового интегрирования или по крайней мере эквивалентного линейного метода, последний является адекватным для оценки максимального ускорения, НДС, но не способен оценить остаточное НДС после сейсмического воздействия. В динамических расчетах, как несвязанных, так и связанных с нелинейными определяющими моделями грунтов, основная трудность заключается в использовании этих моделей, способных достоверно описывать при динамических циклических нагрузках следующие аспекты: а) реакцию НДС на соответствующее обратное преобразование модели нагрузки; многие модели усиливают гистерезис колебаний по сравнению с данными экспериментов; б) рассеивание колебаний в системе «плотина - основание», вызванное пластическими деформациями и вязкими явлениями; в) явления циклической деградации (например, изменчивость НДС с ростом количества циклов); г) фазовые эффекты взаимосвязи (достоверный прогноз избыточного порового давления воды имеет определяющее значение для прогноза максимального ускорения и оценки постсейсмических эффектов). Поэтому нельзя считать, что связанные расчеты с использованием упруго-пластичных моделей грунтов являются более надежными, чем несвязанные, учитывая, что последние имеют критерии прогноза порового давления. Граница возможностей этого моделирования представлена трехфазной средой с возможностью описания поведения частично водонасыщенных грунтов. Если плотина находится в высокосейсмичном районе, свойства грунтов обычно определяются с помощью специальных экспериментов, таких как полевые динамические (на плотине и основании) и лабораторные (грунтов плотины и основания) испытания, обеспечивающих получение параметров деформации с требуемой точностью. Менее надежным является прогноз параметров деформации, используемых в статическом расчете. Эти параметры сложно получить из полевых геотехнических испытаний, и они имеют низкую надежность при получении из лабораторных испытаний (например, трехосных). Относительно затухания колебаний, учитываемого в линейном эквивалентном и нелинейном расчетах, то общие оценки этого параметра вызывают высокий разброс результатов расчета. Численные расчеты земляных плотин Уровень сложности динамического расчета для прогноза поведения земляных плотин при сейсмических воздействиях выше относительно статических расчетов. Снижение надежности прогнозов объясняется повышенной сложностью определения условий нагружения в расчете и неопределенностями ряда факторов, среди которых наиболее важными являются: а) трудность определения и моделирования отражающих/поглощающих границ расчетной области; б) взаимодействие плотины с грунтом основания: в случае его сжимаемости это сильно меняет сейсмическую реакцию плотины по сравнению с жестким грунтом, необходимо обнаружить в аллювиальном основании прослойки с низкой плотностью и механическими параметрами, но трудно получить эти параметры; в) зависимость механических свойств грунта, определяющих сейсмическую реакцию плотины, от ее статического НДС. При оценке реакции земляной плотины при и после сейсмического воздействия расчет ее НДС составляет основу любого прогноза. Общая устойчивость плотины, возможность разжижения грунта, поведение плотины после сейсмического воздействия, сосредоточенная фильтрация, гидравлический разрыв глинистого ядра и т. п. могут определяться и в псевдодинамических расчетах. С другой стороны, даже если выполняются псевдо-статические расчеты, они не могут достоверно учитывать такие важные явления, как эффект усиления амплитуды колебаний, остаточное избыточное поровое давление воды после сейсмического воздействия, разжижение грунта, перераспределение напряжений. Однако взаимодействие земляной плотины с ВБ существенно не изменяет частотную реакцию плотины. Численные расчеты каменно-земляных плотин с глинистым ядром Для этих плотин с водонепроницаемым глинистым ядром только ядро и дренажи могут быть смоделированы с помощью реологического закона. Механическое поведение дренажей не является важным аспектом, поскольку они узкие, но их гидравлическое поведение следует учесть, так как это определяет граничные условия поведения воды в ядре. Продвинутые модели обеспечивают помощь в анализе риска гидравлического разрыва ядра путем расчета главного минимального напряжения и порового давления воды. Арочные эффекты в ядре (в его верховой или низовой грани), когда боковые зоны слабо деформируемые, или от берега к берегу в каньонах, или сочетание того и другого могут быть проанализированы с помощью этих методов. В настоящее время не существует адекватного закона поведения ядра, описывающего влияние наполнения ВБ на верховую зону плотины. Эта зона оседает при наполнении ВБ, хотя она насыщена водой и поэтому в ней низкие межблочные напряжения. Численные расчеты каменно-насыпных плотин с непроницаемым экраном В каменно-насыпной плотине с непроницаемым экраном вода ВБ обычно не проникает в плотину. Если вода, просачивающаяся из боковых примыканий плотины, и вода ВБ, поступающая в обход экрана, и/или грунтовые воды, поступающие с берегов, удовлетворительно контролируются цементационными и дренажными завесами, плотина может считаться неводонасыщенной. Проблема, которую следует решить, заключается в численных расчетах поведения экрана при статических и динамических нагрузках. Следовательно, принципиально, чтобы численная модель могла использовать реологический закон с деформационным упрочнением, который дает адекватное моделирование поведения каменной насыпи как функции ее напряженного состояния. Модель должна моделировать этап строительства плотины, так как закон деформационного упрочнения описывает историю нагружения материала плотины. 2. Классификация и выбор моделей грунтов в численных расчетах Нелинейно-упругие модели Этот тип моделей грунтов является наиболее простым и широко применяемым. Упругие модули (модуль Юнга, E, и коэффициент Пуассона, µ, или модули объемной деформации, B, и сдвига, G) принимаются нелинейно зависимыми от напряжений. Эти модели включают до девяти параметров, легко определяемых из стандартных трехосных испытаний грунтов. Самая известная среди них - модель Дункана - Чанга [4], которая часто используется в статических расчетах грунтовых плотин, особенно каменно-насыпных с глинистыми ядрами, диафрагмами из асфальтобетона и экранами из железобетона. Главным недостатком этих моделей является то, что все деформации грунта принимаются упругими и направление приращения напряжений зависит от направления приращения деформаций. Однако это допущение несправедливо для грунтовых плотин, в которых во время наполнения и опорожнения ВБ происходит поворот осей главных напряжений и изменение направления нагрузки. Таким образом, эти модели не учитывают влияние пути нагружения и не могут моделировать важный эффект дилатансии (разуплотнения при сдвиге) уплотненных грунтов и их нелинейное поведение при разгрузке и повторной нагрузке. Поэтому эти модели лучше всего использовать для расчетов однородных земляных плотин из несвязных грунтов. Упругие идеально пластичные модели Среди них наиболее упрощенные модели - упругие идеально пластичные модели с критериями прочности Mора - Кулона, Мизеса и Друкера - Прагера. Известно, что эти приближения дают плохие результаты и не имеют реальных преимуществ в численных расчетах грунтовых плотин. Дилатансия несвязных грунтов контролирует предельные нагрузки грунтовых плотин. Тем не менее прогноз, выполненный с критериями прочности Мора - Кулона с ассоциированным правилом течения, нереально высокий. Упруго-пластичные модели с упрочнением (модели УПУ) Изотропные, кинематические и изотропно-кинематические правила упрочнения вследствие пластических объемных и сдвиговых деформаций приводят к моделям разной сложности. Модель критического состояния (КС или CS), разработанная Роско, Шофилдом и Берландом [5-6], модель бесконечного числа поверхностей (БЧП или INS), разработанная Мрозом, Норрисом и Зенкевичем [7], модели Лейда, Превоста и Зарецкого [8-11] - ряд примеров в этой группе моделей. На данном этапе разработки моделей только модели УПУ могут достоверно описывать поведение грунтов. Список этих моделей с диапазоном их применимости, типов нагрузок и формулировки моделей приведены в нашей статье [12]. Хорошо разработанная модель КС (модифицированная модель Кэм-Клей) является основой около 30 моделей, предложенных для монотонных и циклических нагрузок. В формулировках моделей указаны правила течения (связанные или несвязанные) и правила упрочнения (изотропные или кинематические), так как это оказывает большое влияние на компьютерную реализацию моделей. Модель КС, разработанная для нормально уплотненных (Н-У или N-C) и слабо переуплотненных (коэффициент переуплотнения не более 2,0) глин, представляет собой простую модель с 4-5 параметрами и способна качественно предсказать основные аспекты поведения грунтов. Далее рассмотрены наиболее распространенные в практике проектирования грунтовых плотин математические модели Кэм-Клей, применяемые в статических расчетах этих плотин. Модифицированная модель Кэм-Клей (МКК) Модель МКК была первоначально разработана для трехосных условий нагрузки. Экспериментальные измерения на мягких глинах послужили основой для развития конститутивной модели, выражающей изменение коэффициента пористости е (объемная деформация εv) в зависимости от логарифма эффективной длины напряжения , как видно двух графиков на рис. 1. Оба графика связаны следующим образом: λ* = λ / (1 + e); k* = k / (1 + e), где k - наклон набухания; λ - наклон нормали уплотнения; e - текущий коэффициент пористости. Первый график состоит из нормали уплотнения (NCL) и ряда кривых набухания. При первичной загрузке грунт деформируется вниз NCL. Далее он уплотняется до определенного напряжения (давления предварительного уплотнения pc) и разгружается до текущего набухания. Затем после перезагрузки грунта на первых этапах он деформируется к кривой набухания и достигает напряженного состояния (параметра pc), который был до разгрузки. В этот момент грунт начинает деформироваться до NLC (первичная нагрузка - кривая сжатия). Параметры κ и λ определяются по формулам κ = 1,3Сс (1 - υ) / (1+ υ); λ = Сs / 2,3; где Сс - коэффициент одномерной сжимаемости; Сs - коэффициент одномерного расширения; κ - угловой коэффициент линии критического состояния; λ - угловой коэффициент нормали уплотнения; υ - коэффициент Пуассона. Линия сжатия первичной загрузки [Compression line of primary loading] Кривая набухания зазгрузка/перегрузка [Curve of swelling unloading/overloading] NCL - нормаль уплотнения [Consolidation normal] Рис. 1. Реакция грунта при изотропном уплотнении (основополагающая модель) [Figure 1. Soil response for the isotropic consolidation (fundamental model)] Сверхкритическое расширение [Supercritical expansion] Окончание изотропного уплотнения [End of isotropic consolidation] Линия критического состояния [Critical state line] Докритическое сжатие [Subcritical compression] Уплотнение [Consolidation] Размягчение [Softening] Рис. 2. Отображение функции текучести в меридиональных и девиаторных плоскостях [Figure 2. Representation of the yield function in meridional and deviatoric planes] Поверхность текучести является гладкой без возможности развитии растягивающих напряжений. Модель МКК позволяет, в отличие от других моделей, прямое моделирование деформационного упрочнения или размягчения для нормально уплотненных или переуплотненных грунтов, учитывает нелинейную зависимость объемной деформации грунта от эффективного среднего напряжения и предельные условия идеальной пластичности. При использовании модели МКК грунт нагружается и при сдвиге в нем могут возникнуть пластические деформации (на рис. 2 точки 1 и 2 для не затвердевания и размягчения соответственно) до достижения критического состояния (точки 3 и 2 для затвердевания и размягчения соответственно). Далее грунт деформируется при сдвиге при условии идеальной пластичности без изменения e и . После разгрузки предполагается линейная (упругая) деформация грунта. Развитие поверхности текучести (упрочнение/размягчение) определяется текущим значением давления предварительного уплотнения рс: = Pci exp [-∆ / (λ* - k*)], где Pci+1 - текущее давление предварительного уплотнения; ∆ - прирост объемной пластической деформации. Помимо параметров κ и λ, собственного веса грунта и коэффициента Пуассона, модель МКК использует следующие три параметра: Mcs - наклон линии критического состояний; OCR - коэффициент переуплотнения; e0 - начальный коэффициент пористости. Наклон линии критического состояния Mcs определяется из выражений Mcs+30°(φcv) = (2√3 sin φcv) / / (3 - sin φcv); Mcs-30°(φcv) = (2√3 sin φcv) / (3 - sin φcv), где первая формула относится к трехосному сжатию, вторая - к трехосному расширению; φcv - угол внутреннего трения грунта, соответствующий его критическому состоянию. Обобщенная модель Кэм-Клей (ОКК) Эта модель является усовершенствованием модели MКК при моделировании грунта в сверхкритической области (рис. 3), где условие разрушения следует моделям Мора - Кулона, Друкера - Прагера и модифицированной модели Мора - Кулона (MMК). В отличие от модели МКК (штриховая линия на рис. 3), проекция ОКК на девиаторную плоскость имеет тот же вид, что и MMК. Проекция модели МКК и модели Друкера - Прагера имеет вид окружности. В докритической области поведение моделей МКК и ОКК идентично. При разгрузке предполагается линейная реакция грунта. Параметры, которые нужно задавать, применяя ОКК, аналогичны параметрам МКК. Сверхкритическая область [Supercritical domain] Докритическая область [Subcritical domain] Рис. 3. Проекция области пластичности МКК и ОКК на меридианную и девиаторическую плоскости: ---- CSL - линия критического состояния; --- Y.F. GCC - поверхность текучести в обобщенной модели Кэм-Клей (ОКК); …… Y.F. MCC - поверхность текучести в модифицированной модели Кэм-Клей (МКК) [Figure 3. Plasticity projection of MCC and GCC on meridian and deviatographic planes: ---- CSL - critical state line; --- Y.F. GCC - yield surface in generalized Cam Clay Model (GCC); …… Y.F. MCC - yield surface in modified Cam Clay Model (MCC)] Достоверность численных расчетов грунтовых плотин Обобщим проблемы грунтовых плотин, которые пока не могут быть решены путем математического моделирования (индексы достоверности R.I. 3 и 4 подчеркнуты). 1. Внутренняя суффозия при строительстве, первом наполнении водохранилища и эксплуатации (статическое и постсейсмическое поведение, R.I. = 4). 2. Гидравлический разрыв глинистого ядра при первом наполнении водохранилища и эксплуатации (статическое и постсейсмическое поведение, R.I. = 3). 3. Оценка НДС при строительстве, первом наполнении водохранилища и эксплуатации (динамическое поведение, R.I. = 3). 4. Оценка общей устойчивости примыканий и откосов при строительстве, первом наполнении водохранилища и эксплуатации (динамическое поведение, R.I. = 3). 5. Разжижение грунта основания при строительстве, первом наполнении водохранилища и эксплуатации (динамическое поведение, R.I. = 3). 6. Постсейсмическая фильтрация при строительстве, первом наполнении водохранилища и эксплуатации (динамическое поведение, R.I. = 3). Таблица 1 Достоверность численных расчетов грунтовых плотин: статическое поведение Проблемы безопасности плотины Элементы системы «плотина - основание - водохранилище» Индексы достоверности, RI Строительство Наполнение ВБ Эксплуатация Тип плотины: a, b, c, d, e Тип плотины: a, b, c, d, e Тип плотины: a, b, c, d, e НДС (смещения, деформации, напряжения) Тело плотины; основание RI = 2 RI = 2 a, b, c, d: RI = 2; e: RI = 1 RI = 2 RI = 1 RI = 2 Общая устойчивость плотины Тело плотины; примыкания, откосы RI = 1 RI = 2 RI = 1 RI = 2 RI = 1 RI = 2 Окончание табл. 1 Проблемы безопасности плотины Элементы системы «плотина - основание - водохранилище» Индексы достоверности, RI Строительство Наполнение ВБ Эксплуатация Тип плотины: a, b, c, d, e Тип плотины: a, b, c, d, e Тип плотины: a, b, c, d, e Фильтрация в плотине и основании Тело плотины; основание RI = 2 RI = 2 RI = 2 RI = 1 RI =2 Внутренняя суффозия грунта Тело плотины; основание RI = 4 RI = 4 RI = 4 RI = 4 RI = 4 Гидравлический разрыв глинистых элементов Тело плотины; основание RI = 3 a, b, d: RI = 3 RI = 3 a, b, d: RI = 3 RI =3 Статическое разжижение грунта Основание плотины RI = 2 RI = 2 RI = 2 Растворение грунта при фильтрации Основание плотины RI = 4 RI = 4 RI = 4 Примечания. Индексы RI (1, 2, 3, 4) даны в матрице [i, j], i = 1, …, 12 - число рядов, где рассмотрена проблема; j = 1, …, 3 - число столбцов, где даны периоды службы плотины. Типы грунтовых плотин: a - однородная плотина; b - зонированная плотина; c - плотина с негрунтовым экраном или диафрагмой; d - зонированная каменно-земляная плотина; e - каменно-насыпная плотина с негрунтовым экраном или диафрагмой. Table 1 Reliability of numerical analyses of embankment dams: static behavior Problems of dam safety Elements of system “dam - foundation - reservoir “ Indices of reliability, RI Construction Filling of reservoir Operation Type of dam: a, b, c, d, e Type of dam: a, b, c, d, e Type of dam: a, b, c, d, e Stress-strain state (displacements, deformations, stresses) Dam body; foundation RI = 2 RI = 2 a, b, c, d: RI = 2; e: RI = 1 RI = 2 RI = 1 RI = 2 General stability of dam Dam body; abutments, slopes RI = 1 RI = 2 RI = 1 RI = 2 RI = 1 RI = 2 Seepage in dam and foundation Dam body; foundation RI = 2 RI = 2 RI = 2 RI = 1 RI = 2 Inner piping of soil Dam body; foundation RI = 4 RI = 4 RI = 4 RI = 4 RI = 4 Hydraulic fracturing of clay elements Dam body; foundation RI = 3 a, b, d: RI = 3 RI = 3 a, b, d: RI = 3 RI = 3 Static liquefaction of soil Dam foundation RI = 2 RI = 2 RI = 2 Dilution of soil during seepage Dam foundation RI = 4 RI = 4 RI = 4 Notes. Indices RI (1, 2, 3, 4) are given in matrix [i, j], i = 1, …, 12 - number of rows, where problem is considered; j = 1, …, 3 - number of columns, where periods of dam service are given. Types of embankment dams: a - homogeneous dam; b - zoned dam; c - dam with non-soil face or diaphragm; d - rockfill dam with clay core; e - rockfill dam with non-soil face or diaphragm. Таблица 2 Достоверность численных расчетов грунтовых плотин: динамическое поведение Проблемы безопасности плотины Элементы системы «плотина - основание - водохранилище» Индексы достоверности, RI Строительство Наполнение ВБ Эксплуатация Тип плотины: a, b, c, d, e Тип плотины: a, b, c, d, e Тип плотины: a, b, c, d, e НДС (смещения, деформации, напряжения) Тело плотины; основание RI = 2 RI = 2 a, b, c, d: RI = 2; e: RI = 1 RI = 2 RI = 1 RI = 2 Общая устойчивость плотины Тело плотины; примыкания, откосы RI = 1 RI = 2 RI = 1 RI = 2 RI = 1 RI = 2 Фильтрация в плотине и основании Тело плотины; основание RI = 2 RI = 2 RI = 2 RI = 1 RI = 2 Внутренняя суффозия грунта Тело плотины; основание RI = 4 RI = 4 RI = 4 RI = 4 RI = 4 Гидравлический разрыв глинистых элементов Тело плотины; основание RI = 3 a, b, d: RI = 3 RI = 3 a, b, d: RI = 3 RI = 3 Статическое разжижение грунта Основание плотины RI = 2 RI = 2 RI = 2 Растворение грунта при фильтрации Основание плотины RI = 4 RI = 4 RI = 4 Примечания к табл. 2 аналогичны примечаниям к табл. 1. Table 2 Reliability of numerical analyses of embankment dams: dynamic behavior Problems of dam safety Elements of system “dam - foundation - reservoir” Indices of reliability, RI Construction Filling of reservoir Operation Type of dam: a, b, c, d, e Type of dam: a, b, c, d, e Type of dam: a, b, c, d, e Stress-strain state (displacements, deformations, stresses) Dam body; foundation RI = 2 RI = 2 a, b, c, d: RI = 2; e: RI = 1 RI = 2 RI = 1 RI = 2 General stability of dam Dam body; abutments, slopes RI = 1 RI = 2 RI = 1 RI = 2 RI = 1 RI = 2 Seepage in dam and foundation Dam body; foundation RI = 2 RI = 2 RI = 2 RI = 1 RI = 2 Inner piping of soil Dam body; foundation RI = 4 RI = 4 RI = 4 RI = 4 RI = 4 Hydraulic fracturing of clay elements Dam body; foundation RI = 3 a, b, d: RI = 3 RI = 3 a, b, d: RI = 3 RI = 3 Static liquefaction of soil Dam foundation RI = 2 RI = 2 RI = 2 Dilution of soil during seepage Dam foundation RI = 4 RI = 4 RI = 4 Notes for Table 2 are similar to notes for Table 1. Заключение Модифицированная математическая модель Кэм-Клей (МКК-модель) до сих пор является самой проверенной и применяемой из всех упруго-пластических моделей грунтов в численных расчетах грунтовых плотин, массивов и оснований. Поэтому в численных расчетах грунтовых плотин и оснований рекомендуется использовать модифицированную или обобщенную модель Кэм-Клей и одну-две более сложные модели (Зарецкого, Превоста или Лейда). Подготовка плотины к пуску в эксплуатацию [Preparing dam for commissioning] Корректировка некоторых расчетных параметров изменения проекта [Adjustment of some estimated parameters of project change] Различия < Допуски [Differences < Tolerances] Нет [No] Сравнение расчетных и замеренных эквивалентных параметров поведения плотины [Comparison of calculated and measured equivalent parameters of dam behavior] Данные мониторинга поведения плотины [Monitoring data of dam behavior] Численные расчеты поведения плотины [Numerical analyses of dam behavior] Расчетные предположения [Estimated assumptions] Да [Yes] Рис. 4. Блок-схема взаимодействия между данными мониторинга и численного моделирования поведения плотины при первом наполнении водохранилища (согласно СИГБ) [Figure 4. Flow chart of the interaction between monitoring data and numerical modelling of dam behavior during the first filling of reservoir (as per ICOLD)] Ограничения эксплуатации водохранилища. Обеспечение безопасности плотины [Restrictions on reservoir operation. Ensuring dam safety] При опасном поведении плотины запуск системы оповещения в НБ [In case of dangerous dam behavior launching of alert system in downsteam region] Продолжение нормальной эксплуатации плотины [Continuation of normal
Об авторах
Юрий Петрович Ляпичев
АО «Институт “Гидропроект”»; Международная комиссия по большим плотинам (ICOLD)
Автор, ответственный за переписку.
Email: lyapichev@mail.ru
SPIN-код: 3096-6362
эксперт АО «Институт “Гидропроект”» по зарубежным проектам, член Комитета СИГБ по компьютерным аспектам расчета и проектирования плотин, доктор технических наук, профессор
Французская Республика, 75016, Париж, Avenue Kleber, 61Список литературы
- ICOLD Bulletins. No. 122. Computational procedures for dams. ICOLD Edition, Paris; 2001.
- ICOLD Bulletins. № 155. Guidelines for use of numerical models in dams. ICOLD Edition, Paris; 2013.
- Newmark N.M., Rosenblueth E. Fundamentals of Earthquake Engineering. Prentice-Hall, Inc., Englewood Cliffs, N. J.; 1971.
- Duncan J.M., Chang Y.Y. Nonlinear analysis of stress and strain in soils. Journal of the Soil Mechanics and Foundations Division. 1970;96(5):1629-1653.
- Mroz Z., Norris V., Zienkiewicz O. Anisotropic hardening model for soils and its application to cyclic loading. Int. J. Num. & Anal. Methods in Geomechanics. 1978;2:203-221.
- Roscoe K.H., Schofield A.N. Mechanical behaviour of an idealized ‘wet clay’. Proc. of 2nd European Conf. on Soil Mechanics (Wiesbaden). 1963;1:47-54.
- Roscoe K.H., Burland J.B. On the generalized stress-strain behaviour of ‘wet clay’. In: Heyman J., Leckie F.A. (eds.) Engineering Plasticity. Cambridge University Press, Cambridge; 1968. p. 535-609.
- Lade P.V. Elastoplastic stress-strain theory for cohesionless soil with curved yield surfaces. Int. Journal of Solids and Structures. 1977;13:1019-1035.
- Prevost J.H. Plasticity theory for soil stress-strain behavior. Journal of the Engineering Mechanics Division. 1978;104(5):1177-1194.
- Prevost J.H. Anisotropic undrained stress-strain behavior of clays. Journal of the Geotechnical Engineering Division. 1978;104(8):1075-1090.
- Zaretsky Yu. Soil viscoplasticity and design of structures. Balkema, Holland; 1996.
- Lyapichev Yu.P. Ocenka dostovernosti matematicheskih modelej gruntov dlya chislennyh raschetov povedeniya gruntovyh plotin [Estimation of reliability of mathematical models of soils for numerical calculations of the behavior of soil dams]. RUDN Journal of Engineering Researches. 2000;(3):110-115. (In Russ.)