При подъеме УВБ от отм. 500 до отм. 540 м (мм)

Расчет Епл=30000 МПа	Натурные измерения 1993 г.	Натурные измерения 1994 г.
18 секция	33 секция	45 секция	18 секция	33 секция	45 секция	18 секция	33 секция	45 секция
			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 циклов натурных измерений некоторого диагностического показателя (например, приращения радиальных перемещений dХ) в I точках и введем следующую традиционную меру близости D расчетных диагностических показателей (расчетных приращений перемещений) к измеренным:

D = S( - )²,

где — измеренные значения приращений перемещений в 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