Глубина раскрытия швов (м)

Вид шва	УВБ (м)
	510	517	524	531	536	540
Контактный	2	9	13	17	33	27
Строительный на отм. 355 м	-	-	-	-	6	22

Выполненная корректировка расчетной модели позволила существенно сблизить результаты расчетов и данные натурных измерений. В частности, расчеты показали, что в системе плотина-основание со стороны верховой грани на отм. 308 (контакт плотины с основанием) и на отм. 350 м образуются две магистральных трещины. Характер трещинообразования в расчетах и в натуре оказался весьма схожим.

Таким образом, основное отличие линейно упругой проектной модели для плотины Саяно-Шушенской ГЭС от выявленной на стадии начальной эксплуатации реальной работы плотины заключалась в том, что не было учтено образование и развитие горизонтальных трещин со стороны верховой грани вблизи контакта плотины с основанием. В дальнейшем (начиная с 1996 года) все поверочные расчеты плотины выполнялись в неупругой постановке с учетом образования и развития трещин в теле плотины и ее основании.

3.2. Калибровка расчетных моделей

В п. 4.8 «Методики» описана общая схема калибровки (уточнения физико-механических характеристик материалов плотины и основания, заложенных в детерминистические расчетные модели). Рассмотрим возможные процедуры калибровки на примере плотины Саяно-Шушенской ГЭС.

Первые сопоставления данных натурных измерений с результатами расчетов, выполненных на стадии технического проекта и рабочих чертежей, показали, что расчетные приращения перемещений при росте УВБ систематически превосходят измеренные примерно на 30% (см. таблицу П.П.6).

Таблица П.П.6

Приращения радиальных перемещений гребня 18-ой, 33-ей, 45-ой секций плотины при подъеме увб от отм. 500 до отм. 540 м (мм)

Расчет Епл=30000 МПа			Натурные измерения 1993 г.			Натурные измерения 1994 г.
18 секция	33 секция	45 секция	18 секция	33 секция	45 секция	18 секция	33 секция	45 секция
93	115	97	66,1	90,9	59,8	61,7	83,1	51,8

Из таблицы видно, что для ключевой секции 33 расчетные перемещения больше измеренных в 1993 г. — в 115/90,9 =1,26 раза, в 1994 г. — в 115/83,1 = 1,38 раза. На основании подобных сравнений экспертами (специалистами Ленгидропроекта и СШ ГЭС) был сделан вывод о том, что реальная жесткость (модуль деформации) бетона плотины выше проектной. Проектный модуль деформации бетона плотины был принят равным 25000 МПа в соответствии с рекомендациями действовавших норм СНиП II-54-77 «Плотины бетонные и железобетонные». На основании сделанного вывода, начиная с 1994 г., во всех расчетах конструкционный модуль деформации бетона плотины был принят равным 32000 МПа. Однопараметрическое увеличение модулей деформации плотины и основания в 1,3 раза в линейно упругой задаче приводит к уменьшению расчетных перемещений во столько же раз. Такая простейшая процедура калибровки позволила существенно сблизить данные расчетов и натурных измерений. В этом легко убедиться, разделив расчетные перемещения таблицы П.П.6 на 1,3.

Рассмотрим более строгую и формализованную процедуру калибровки.

Примечание. Проектные расчеты выполнялись в трехмерной постановке. В рассматриваемом примере для простоты далее использована двухмерная расчетная модель теории оболочек средней толщины (см. Л.А. Гордон, И.К. Соколовский, Л.Х. Цовикян «Прогноз перемещений арочной плотины на основе идентифицированной расчетной модели», — Известия ВНИИГ, т. 214, 1989 г. В этой модели были приняты следующие модули деформации материалов плотины и основания: Е_пл = 30000 МПа, Е_осн = 20000 МПа.

Выберем J циклов натурных измерений некоторого диагностического показателя (например, приращения радиальных перемещений Х) в I точках и введем следующую традиционную меру близости D расчетных диагностических показателей (расчетных приращений перемещений) к измеренным:

D = ( - )²,

где — измеренные значения приращений перемещений в i-ой точке на j-ом цикле измерений; — расчетные значения перемещений в i-ой точке, соответствующие УВБ и температурам j-го цикла.

Если измеренные и расчетные показатели совпадают, то D = 0. Чем сильнее различаются измеренные и расчетные показатели, тем больше D. Например, для данных измерений, приведенных в таблице П.П.6, где два цикла измерений (мах УВБ 1993 и 1994 гг.) и три точки (гребень 18-ой, 33-ой,45-ой секций), функционал D будет:

D = {[F(18)-66,l]²+[F(33)-90,9]²+[F(45)-59,8]²+[F(18)-61,7]²+[F(33)-83,1]²+[F(45)-51,8]²},

где величины F(18), F(33), F(45) — расчетные приращения перемещений гребня 18-ой, 33-ей, 45-ой секций при подъеме УВБ от отм. 500 м до отм. 540 м. Для расчетной модели, данные которой приведены в таблице П.П.6 (Е_пл = 30000 МПа, Е_осн = 20000 МПа), значения этих параметров были: F(18) = 93 мм, F(33) = 115 мм, F(45) = 97 мм.

Вычислим меру погрешности D₀, подставив в D расчетные значения F(18), F(33), F(45):

D₀ = (93-66,1)² + (115-90,9)²+(97-59,8)²+(93-61,7)²+(115-83,1)²+(97-51,8)² = 723,6 + 580,8 + 1383,8 + 979,7 + 1017,6 + 2043,0 = 6728,5

Откалибруем расчетную модель (подберем такие значения Е_пл и Е_осн, чтобы мера погрешности D уменьшилась), приняв за нулевое приближение приведенные проектные значения приращений перемещений.

В соответствии с процедурой, описанной в п. 4.8 «Методики», проварьируем Е_пл, Е_осн: дадим Е_пл, Е_осн приращения, увеличив их однопараметрически, например, в 1,2 раза, т.е. примем:

=30000  1,2 = 36000 МПа; = 20000  1,2 = 24000 МПа

Поскольку все расчетные модули деформации в линейно деформируемой системе увеличились в 1,2 раза, то все расчетные перемещения уменьшились в 1,2 раза. Поэтому, не пересчитывая систему плотина-основание, можно получить расчетные перемещения при модулях , делением на 1,2 данных расчета таблицы П.П.6 (см. таблицу П.П.7):

Таблица П.П.7