Состояние вопроса и задачи исследований. Цель исследования: разработать методику расчета долговечности железобетонных конструкций для климатических условий прибрежной зоны морей Дальнего Востока от комплексного воздействия карбонизации и хлоридной агрессии. Задачи исследования:Проанализировать результаты исследований комплексного воздействия карбонизации и хлоридной агрессии на морской бетон;Усовершенствовать теоретические модели карбонизации и хлоридной агрессии с учетом зарождения и развития трещины в защитном слое бетона;Разработать вероятностный метод прогнозирования срока службы железобетонных портовых конструкций с учетом комплексного воздействия карбонизации и хлоридной агрессии;Экспериментально исследовать техническое состояние железобетонных элементов эксплуатируемых портовых сооружений, и выявить особенности их деградации в морской среде, в том числе, с учетом технологических факторов;Верификация детерминистической модели расчета совместного действия карбонизации и хлоридной агрессии на морской бетон на основе 2-го закона Фика.Принята модель решения дифференциального уравнения диффузии, выведенная J. Crank, с учетом влияния карбонизации на перенос хлор-ионов в бетоне. Основное уравнение имеет вид:∂Cf∂t=DCl*1+(1we)∙(αL(1-d∙ac)(1+βL∙Cfb)2)  ∂2Cf∂x2               (1)где αL  и βL  - эмпирические константы; ac  – степень карбонизации бетона; d  – коэффициент понижения связывающей способности хлорида за счет карбонизации; Ccl  – общее содержание хлорида в бетоне; Cf – содержание свободного хлорида в бетоне; t – время эксплуатации; b  – масса вяжущего; x  – глубина защитного слоя. Степень карбонизации бетона определяется из пропорции %Xc - 100% ac %КС - x%ac где Xc – предельная величина карбонизации,КС - карбонатная составляющая  КС=mCaCO3mн100%                   (2)Здесь mCaCO3  – масса карбонатной составляющей; mн  – масса навески пробы, определяется экспериментально.Предельная величина карбонизации Xc  определяется по формуле Xc(t)=2∙DCO2ac∙1tfTt∙fWt∙CCO2tdt∙(t0t)0.12  (3)     где fTt, fWt, CCO2t  – функции влияния температуры, влажности и концентрации CO2  на коэффициент диффузии; ac  – коэффициент, определяющий количество CO2 , необходимое для превращения всех способных карбонизироваться продуктов гидратации; DCO2  – начальный коэффициент диффузии углекислого газа в бетоне.Для учета влияния СО2 выполнена оценка содержания углекислого газа в воздухе с учетом срока службы железобетонных конструкций, согласно данным кривой Килинга. Поскольку концентрация хлоридов в морской среде изменяется в зависимости от погодных условий, модель воздействия морской воды на сооружения модифицирована вводом зависимости от расстояния между сооружением и берегом [1–4].Для верификации модели совместного действия карбонизации и хлоридной агрессии взяты железобетонные конструкции классов XC4 и XS3 по условиям эксплуатации со средними значениями параметров бетонной смеси согласно EN 206:2013 и минимальными толщинами защитного слоя бетона по СП 28.13330.2012. Согласно предложенной методике, проведены расчеты для данных табл. 1.Для моделирования использована и модифицирована программа в Mathcad (авторы Д. Шестовицкий и Э. Карапетов), выполненная на основе конечно-разностного подхода (рис. 1, табл. 2).Результаты показали, что карбонизация привела к высвобождению хлоридных ионов в поровый раствор, в следствии чего, уменьшается срок службы конструкции.  Методика вероятностного расчёта совместного действия карбонизации и хлоридной агресси.  представлена методика вероятностного расчета совместного воздействия карбонизации и хлоридной агрессии на бетон. В основе лежит уравнение вероятности отказа:Pf=P(R-S≤0)≤P                    (4)где Pf  – вероятность отказа; P  – допустимая вероятность наступления предельного состояния;S – функция нагружения; R – функция сопротивления конструкции.Для хлоридной коррозии, в вероятностной постановке, S  представляет собой значение концентрации хлоридов C(x, t) . Ccrit  – параметр критической (пороговой) концентрации хлорида на уровне залегания арматуры, превышение которого приводит к инициированию коррозии. В этом случае вероятность безотказной работы:Pf=P(Ccrit-C(x, t)≤0)≤P            (5)C(x, t)  определяется с использованием модели, которая базируется на решении 2-го закона диффузии А. Фика с помощью функции ошибки Гаусса (C. Andrade) C1(x,t)=C0n=0∞anerfc2∙n∙e+x2D1(t)∙t-a∙erfc(2n+2∙e-x2D1(t)∙t)                                (6)C2(x,t)=2∙k∙C0k+1n=0∞anerfc(2n+1∙e+k∙x2D1∙t)                                           (7)a=1-k1+k                                                                            (8)k=D1D2                                                                       (9) Модель базируется на различии коэффициентов диффузии в одном срезе («скин-эффект»), что происходит либо в результате ремонтного восстановления защитного слоя бетона, либо при действии множества агрессивных факторов внешней среды на конструкцию. В случаях совместного действия карбонизации и хлоридной агрессии, ур. (6) и (7) преобразуются к виду: Cclcb(x,t)=Csn=0∞anerfc2∙n∙Xc+x2Dcl,cb(t)∙t-a∙erfc(2n+2∙Xc-x2Dcl,cb(t)∙t)                         (10)Cclucb(x,t)=2∙k∙Csk+1n=0∞anerfc(2n+1∙x+k∙(x-Xc)2Dcl,cb∙t)                                          (11) где Cs  – поверхностная концентрация хлоридов, %; x  – толщина защитного слоя бетона, мм; erfc (x) – обратная функция ошибок  Гаусса; Dcl,cb  – коэффициент диффузии хлоридов карбонизированного бетона; Xc =Xc(t) -глубина карбонизации бетона, мм; t – время, годы;  a и k – коэффициенты из ур. (8) и (9); Dcl,ucb  – коэффициент диффузии хлоридов некарбонизированного бетона.  Таблица 1Исходные данные конечно-разностной моделиПараметр              Ед. изм.Участок о. СахалинСеверный Центральный Южный Tmax °C18,320,517,7Tmin °C-7,3-6,2-2,4W max%868185W min%747671w/b 0,40,40,4bкг/м3 350350350Карбонизацияge,-2,52,52,5fe-555ЕкДж/моль 404040RкДж/К 8,314·10-38,314·10-38,314·10-3СSкг/м33,890·10-43,890·10-43,890·10-4DCO2см2/с3,399·10-43,399·10-43,399·10-4nm-0.120.120.12Хлоридная агрессияЕкДж/моль 41,841,841,8RкДж/К 8,314·10-38,314·10-38,314·10-3αL –0,11850,11850,1185βL -0,090,090,09Wref%656565Cenv(L) кг/м36,26,26,2m-0,40,40,4t0дней (лет)28 (0,0767)28(0,0767)28(0,0767)tлет505050 Таблица 2Результаты моделированияПараметрЕд. изм.Место эксплуатацииСеверный Центральный Южный Фронт карбонизации (t=50 лет)мм30,829,629,4Степень карбонизации-0,610,60,6Концентрация хлоридов на глубине залегания арматуры без учета карбонизации (при t=50 лет)%0,550,650,44Тоже с учетом карбонизации (при t=50 лет)%0,650,60,54Время инициации хлоридной коррозии без учета карбонизации годы504043Тоже с учетом карбонизации годы453530   – концентрация ионов хлорида на глубине защитного слоя X см без учета карбонизация в зависимости от времени t, кг-м3.  - концентрация ионов хлорида на глубине защитного слоя X см с учетом карбонизации, в зависимости от времени t, кг-м3. Критическая концентрация хлорида принята 0,4 % или 1,4 кг/м3 по массе вяжущего Рис.1. График изменения концентрации хлоридов в защитном слое бетона с учетом и без учета карбонизации (южный участок)  В исследовании также учитывается сопротивление между слоями, которое возникает в результате разницы коэффициентов диффузии в одном срезе:Cclucbx,t= =2∙k∙Cs∙Rk+1n=0∞anerfc(2n+1∙x+k∙(x-Xc)2Dcl,cb∙t)     (12)где R  – сопротивление между слоями; C(x, t)  с учетом действия карбонизации рассчитывается, как система из ур. (10) и (12). В расчетную модель комбинированного действия карбонизации и хлоридной агрессии входит ряд базисных переменных. Предложены рекомендации по средним значениям этих переменных и их типам распределения. Для южной части о. Сахалин в табл. 3 даны их средние значения, стандартное отклонение и тип распределения. Таблица 3Данные для вероятностного моделирования концентрации хлоридов Южная часть о. Сахалин Тип распределенияСреднее значениеСтандартное отклонение ПараметрЕд. изм. Cs %Const2,5– x ммConstвектор из множества 0-50 – Dcl,cb0 м2/сNormal11,689·10-121,2·10-12 Dcl,ucb0 –Normal 2,387·10-121,2·10-12 ke –Normal0,670,05  КNormal-–  КConst273– kt –Normal0,800,05 kc –Normal10,125 t0 годConst0,0767– tгодConstвектор из множества t0-50 – ncl –Beta0,3a=0; b=1 Ccrit %Normal0,40,063  Для аналитического решения прямой задачи определения вероятности ресурсного отказа и обратной задачи определения процентного ресурса конструкций использовано имитационное моделирование с вычислением необходимых функционалов, например, содержания хлоридов на заданной глубине, срока службы и т.п. Для компьютерной реализации в программе Matlab написан код для расчета по модели проникновения хлоридов с учетом эффекта карбонизации. Результат расчета программы - вероятности ресурсного отказа конструкции и индексы надежности в течение срока службы для различных значений толщин защитного слоя бетона. На первом этапе программа определяет глубины карбонизации и изменение концентрации хлоридов по глубине (рис. 2–4). Далее, после построения графиков рис. 2–4 программа вычисляет вероятность отказа железобетонной конструкции и ее индекс надежности (рис. 5 и 6, табл. 4). Для 50 лет эксплуатации на южном участке о. Сахалин вероятность отказа получена pf = 98 %. Рис. 2. Рост глубины карбонизации с течением времени Рис. 3. Изменение концентрации  хлоридов в приарматурной зоне для всего срока эксплуатаии  Рис. 4. Профиль концентрации  хлоридов в приарматурной зонев последний год эксплуатации (50 лет)   Рис. 5. Вероятность отказа конструкции   Рис. 6. Индекс надежности конструкции   Таблица 4Вероятность отказа и индекс надежности железобетонных конструкций в зависимости от срока эксплуатации для южного участка о. СахалинСрок эксплуатацииВероятность отказа, pfИндекс надежности, β100,00013,688200,0051,102300,19-0,173400,849-1,02500,981-1,66Результаты натурных исследований Работы включали: визуальный осмотр, определение критических элементов и областей, определение прочности и толщины защитного слоя бетона, выбор тестовых зон по результатам измерений.В тестовых зонах проведены: визуальный осмотр для выбора мест тестирования и отбора проб, проверка глубины защитного слоя бетона, выбор мест отбора проб в «наихудших местах» конструкций, определение глубины карбонизации способом фенолфталеиновой пробы в 6 и более местах (на свежеобломанных и просверленных участках); отбор пластин для хлоридных профилей: (минимум шесть пластин в каждой тестовой зоне минимального размера 70×70 мм и минимальной глубины 50 мм). В лаборатории определены: с помощью ионоселективного электрода значение концентрации хлоридов по глубине образцов, способом фенолфталеиновой пробы – глубина карбонизации.Результаты моделирования совместного действия карбонизации и хлоридной агрессии и их сравнение с экспериментальными данными (табл. 5). Таблица 5Сравнение расчетных и экспериментальных результатовТип Место отбора пробыВозраст конструкцииГлубина защитного слоя, ммИзмеренная концентрация Cl, % Расчетная концентрация Cl, % Концентрация Cl, % по вероятностной моделиХолмский морской торговый портХ4Железобетонная колонна  пешеходного моста33102,242,802,20201,972,201,90301,641,811,62401,101,201,03500,510,580,49Х5Железобетонная балка пешеходного моста33102,252,282,21201,982,211,94301,681,821,62401,101,201,03500,500,580,49Корсаковский морской торговый портК4Фундамент под знак СНО44102,322,762,13201,812,191,77301,461,771,44401,101,171,07500,550,630,53К5Железобетонная опорная под трубы.46102,322,792,15201,811,221,80301,481,811,50401,101,231,10500,550,650,55 Получена хорошая сходимость с вероятностной моделью и удовлетворительная с конечно-разностной (последняя не учитывает «скин-эффект»). Коэффициент диффузии в конечно-разностной модели постоянен на всей глубине защитного слоя, т.е. нет разбивки на 2 слоя с разной диффузией и нет сопротивления между этими слоями, поэтому концентрации хлоридов в приарматурной зоне завышена. Однако, эта модель может быть применима для расчета совместного действия карбонизации и хлоридной агрессии при малой глубине карбонизации, когда «скин-эффект» не оказывает значительного влияния на концентрацию хлоридов в приарматурной зоне.Таким образом, при незначительной глубине карбонизации (до 8 мм) расчет можно вести по конечно-разностной модели, при значительной – значения завышены и расчет следует вести по вероятностной модели. На практике при глубине карбонизации до 8 мм эффект совместного действия карбонизации и хлоридной агрессии не учитывается, и используется модель DuraCrete. Таким образом, наиболее точная модель – вероятностная.Вероятностная модель для определения параметров для ремонта защитного слоя бетона. В этом случае предполагается, что необратимые последствия, приводящие к хлоридной коррозии арматуры в карбонизированном бетоне, могут начинаться уже при концентрации ионов хлорида 0,2 %. Используя это значение как критическое, при котором необходимо проводить ремонт защитного слоя бетона, а также используя модель для двойной среды ур. (10) и (12), разработана программа расчета среднего времени и глубины ремонта поврежденного защитного слоя, которая также позволяет прогнозировать срок службы конструкции, но уже с учетом ремонта (табл. 4 и рис. 7). В качестве материала для ремонта выбран раствор, аналогичный исходному составу бетона.После 50 лет эксплуатации в наиболее неблагоприятном районе о. Сахалин по условиям эксплуатации, вероятность отказа составилаpf =58 %. Таким образом, ремонт конструкции способом замены карбонизированного слоя новым с аналогичными характеристиками увеличивает долговечность. Например, в конструкции, в которой инициировалась коррозия после 29 лет эксплуатации, новое время инициации составляет 45 лет. Таблица 4Срок службы железобетонной конструкции с учетом ремонта защитного слоя бетонаПараметрЕд. изм.Место эксплуатацииЮжная часть о. СахалинВремя инициации хлоридной коррозии без учета замены карбонизированного слоягод29Срок службы конструкции без учета замены карбонизированного слоягод33Среднее время замены карбонизированного слоягод16Средняя глубина замены карбонизированного слоямм24,5Время инициации хлоридной коррозии с учетом замены карбонизированного слоягод45Срок службы конструкции с учетом замены карбонизированного слоягод49    Рис. 7. Изменение концентрации  хлоридовв приарматурной зоне для всего срока эксплуатаии с учетом ремонта   Рис. 8. Вероятность отказа отремонтированной конструкции  Рис. 9. Индекс надежности отремонтированной конструкции Таблица 5Вероятность отказа и индекс надежности отремонтированной железобетонной конструкции в зависимости от срока эксплуатации для наиболее неблагоприятного района о. Сахалин по условиям эксплуатацииСрок эксплуатацииВероятность отказа, pfИндекс надежности, β100,00014,27200,00083,29300,0122,00400,1090,44500,582-0,57 Заключение.На основе анализа моделей совместного действия карбонизации и хлоридной агрессии защитного слоя бетона и верификации с экспериментальными данными определена модель для оценки долговечности морских железобетонных конструкций, учитывающая следующие факторы: толщину защитного слоя бетона; коэффициенты диффузии хлоридов в карбонизированном и некарбонизированном бетоне; критическое содержание и поверхностное содержание хлоридов, поверхностное содержание CO2, их время воздействия; морские условия; фронт карбонизации и др. Часть значений являются случайными величинами с заданными законами распределения. Разработана методика определения срока ремонта конструкции и глубины ремонта защитного слоя бетона сооружений.Разработана методика прогнозирования долговечности железобетонных конструкций при воздействии агрессивной морской среды с учетом ремонта конструкции и с использованием вероятностной модели расчета.Выполнена верификация результатов расчетов вероятности отказа железобетонных элементов по предложенной вероятностной модели. Рекомендации по практическому использованию результатовРазработанные модели позволяют рассчитывать глубину карбонизации, концентрацию ионов хлорида на заданной глубине, срок эксплуатации конструкции, время ее возможного ремонта и глубину возможного восстановления защитного слоя бетона для прибрежной и шельфовой зоны Дальнего Востока.Разработанную методику прогнозирования долговечности железобетонных конструкций при совместном воздействии карбонизации и агрессивной хлоридсодержащей среды с использованием конечно-разностной модели расчета предлагается использовать при глубине карбонизации до 8 мм. Разработанную методику прогнозирования долговечности железобетонных конструкций при совместном воздействии карбонизации и агрессивной хлоридсодержащей среды с использованием вероятностной модели расчета можно использовать:- при оценке эксплуатационной пригодности (безопасности) при обследовании железобетонных конструкций прибрежных и шельфовых сооружений;- при прогнозировании срока службы вновь проектируемых железобетонных конструкций;- при расчете необходимой толщины защитного слоя бетона проектируемых железобетонных конструкций при заданном сроке службы и условиях эксплуатации;- при расчете срока службы бетона в конкретных условиях эксплуатации.- при прогнозировании срока ремонта эксплуатируемых сооружений;Полученные результаты могут использоваться в проектировании новых сооружений и/или ремонте (реконструкции) существующих сооружений, эксплуатируемых в агрессивных условиях морской среды, а также в учебном процессе.