Оценка напряженного состояния защитной оболочки АЭС при воздействии высоких температур, вызванных горением топлива упавшего самолета

Размещено на http://www.allbest.ru/

Оценка напряженного состояния защитной оболочки АЭС при воздействии высоких температур, вызванных горением топлива упавшего самолета

В.Н. Медведев, А.Н. Ульянов, Александр С. Киселев, Алексей С. Киселев, А.С. Залесов, В.Ф. Стрижов

ИБРАЭ РАН, Москва, Россия

До настоящего времени такое техногенное событие, как падение самолета на здания гермооболочки АЭС рассматривалось только в виде динамической нагрузки. Однако практика показала, что последствия пожара, вызванного горением топлива упавшего самолета, могут привести к более разрушительным последствиям, чем непосредственно удар при падении самолета. Например, при столкновении аэробуса Boeing-737 со зданиями Всемирного торгового центра (WTC1 и WTC2) их обрушение произошло спустя 56 мин и 1 ч 43 мин, соответственно, от воздействия температурной нагрузки (пожара).

Продолжительность пожара и температура горения во многом зависит от типа самолета и количества топлива на борту. Следовательно, при падении самолета на защитную оболочку кроме механического воздействия твердого тела будет оказывать влияние и действие локальных температур. Вопрос влияния температурной нагрузки, вызванной горением топлива упавшего самолета на защитную оболочку, ранее не рассматривался.

В этой связи возникает необходимость оценить последствия горения топлива упавшего самолета на здания гермооболочки и определить, может ли это привести, например, к разгерметизации внутренней защитной оболочки, нарушению целостности конструкции или ее части, выходу из строя отдельных арматурных канатов системы преднапряжения.

В работе представлены некоторые результаты компьютерного моделирования процессов разрушения защитной оболочки АЭС с ВВЭР при локальном разогреве фрагмента ее наружной поверхности на цилиндре и в купольной части, который может иметь место при возникновении пожара вблизи стенки оболочки вследствие, например, падения самолета или терракта.

Нагрузка на конструкции и оборудование, возникающая от горения топлива упавшего самолета, должна определяться в зависимости от:

- количества топлива на борту самолета (при полной загрузке топливные баки аэробуса Boeing-767 вмещают 90 т керосина);

- объема помещений (высота, площадь), количества и размещения проемов, их площадь и высота;

- теплового потока в источнике.

В данном расчетном анализе температурная нагрузка принималась консервативно в соответствии с рекомендациями [1]. При возникновении стандартного пожара в зданиях общего назначения температура воздуха в зависимости от времени определяется по формуле:

где - температура воздуха, определяемая в С;

- время от начала аварии, час.

Результаты расчета приведены на рис. 1.

Особенности построения расчетной модели

Трехмерная конечноэлементная модель защитной оболочки АЭС с ВВЭР-1000, представленная на рис.2 адекватно отображает основные конструктивные особенности сооружения:

наличие отверстий диаметром 3 и 4 метра;

утолщения стенки оболочки в зоне этих отверстий;

закладные детали;

внутреннюю герметизирующую облицовку;

наличие внутренних и наружных слоев арматуры;

реальная трассировка арматурных канатов системы преднапряжения;

реальные значения усилий в арматурных канатах.

Рис. 1. Изменение температуры воздуха при стандартном пожаре в зависимости от времени.

Рассматривались зоны локального разогрева в купольной части в виде сектора, а на цилиндре - в виде пятна, которое на развертке цилиндра имеет форму прямоугольника.

Расчетная модель сечения стенки защитной оболочки (рис. 3) для цилиндрической и купольной части на удалении от зон разогрева составлена из 11 слоев конечных элементов, отличающихся друг от друга своими свойствами, включая внутреннюю металлическую облицовку (1 на рис. 2), слои бетона (2, 4, 6, 7, 9 и 11), слои меридиональной арматуры (3 и 10) и кольцевой арматуры (5 и 8). В зонах разогрева выполнено сгущение сетки конечных элементов для более точного моделирования градиентов температуры вблизи границ этих зон и по толщине (рис.4).

Рис. 2. Трехмерная конечноэлементная модель защитной оболочки АЭС с ВВЭР-1000, которая включает 113006 объемных восьмиузловых конечных элементов и 122316 узлов.

Рис. 3. Расчетная модель сечения цилиндрической стенки защитной оболочки. R - расстояние от оси цилиндра до границы каждого расчетного слоя стенки оболочки, мм.

Конечноэлементная модель защитной оболочки включает 113006 объемных восьмиузловых конечных элементов и 122316 узлов. Расчеты выполнялись с применением программного комплекса CONT, реализующего суперэлементный алгоритм и позволяющего проводить анализ прочности защитной оболочки с учетом нелинейного поведения бетона при образовании в нем трещин вследствие воздействия локальных полей температуры.

Рис. 4. Сгущение конечноэлементной сетки вблизи зоны локального разогрева

Определение изменения температур по сечению стенки

При анализе поведения оболочки в условиях локального разогрева первоначально рассчитывались нестационарные поля температуры, связанные с заданием условий конвективного теплообмена на участке поверхности модели, соответствующей зоне локального разогрева. Расчеты полей температуры выполнены в нелинейной постановке с учетом зависимости теплофизических свойств материалов от температуры [2].

На рис. 5 представлена зависимость коэффициента теплоотдачи от времени, которые задаются в качестве параметров конвективного теплообмена на фрагменте поверхности цилиндрической части или купола в зависимости от температуры среды (см. рис. 1).

Рис. 5. Зависимость коэффициента теплоотдачи от времени.

Характерная картина распределения температуры при разогреве локальной зоны купола представлена на рис. 6 (а), а на рис. 6(б) аналогичная информация для локального разогрева цилиндрической части. На рис. 6(б) часть модели удалена для того, чтобы показать градиент температуры по толщине стенки оболочки. На рис. 7 и рис. 8 представлено распределение температуры по толщине купола и цилиндра в центре зоны локального разогрева в различные моменты времени после начала разогрева.

Рис. 6. Распределение температуры в зоне локального разогрева купольной части на момент времени 8.6 часа (а), распределение температуры в зоне локального разогрева цилиндрической части на моменты времени 7.4 часа (б).

Рис. 7 Распределение температуры по толщине в купольной части (номера шагов соответствуют моментам времени в часах: 1- 0.03, 2- 0.08, 3 - 0.14, 4 - 0.27, 5 - 0.46, 6 - 0.55, 7 - 1.30, 8 - 3.05, 9 - 4.90, 10 - 6.75, 11 - 8.6).

Рис. 8 Распределение температуры по толщине в цилиндрической части (номера шагов соответствуют моментам времени в часах: 1- 0.35, 2- 1.20, 3 - 1.82, 4 - 2.72, 5 - 3.65, 6 - 5.55, 7 - 7.40).

Определение напряженно-деформированного состояния оболочки

Полученная информация о распределении температуры использовалась в качестве нагружающего фактора на следующем этапе расчета НДС с учетом образования трещин в бетоне. Кроме температурного воздействия в качестве нагрузки задавался вектор узловых усилий от системы преднапряжения, обеспечивающей обжатие оболочки. Нагрузка от преднапряжения приложена к срединной поверхности оболочки. При расчете значений усилий учитывается реальное значение усилия натяжения каждого каната, его траектория и потери усилия по длине каната от трения и анкеровки. Средняя величина усилия натяжения каждого из канатов составляет ~800 тонн. На нижнем торце модели задаются условия закрепления типа «заделка» - запрещены перемещения узлов в направлении осей XYZ.

В расчетах учитывались изменения физико-механических характеристик бетона и арматуры.

Расчет НДС защитной оболочки, обжатой системой преднапряжения, выполнялся за несколько шагов (7 шагов решения в случае локального разогрева цилиндрической части и 11 шагов решения в случае разогрева купольной части) с переменным полем температуры, вызванным локальным разогревом. При этом для моделирования образования трещин в бетоне на каждом шаге решения был реализован специальный итерационный алгоритм. После расчета НДС на очередном шаге по нагрузке для каждого «бетонного» конечного элемента проверяется условие его растрескивания от напряжений сжатия (предел прочности 18 МПа) или растяжения (предел прочности 1.5 МПа) в меридиональном и окружном направлениях. Если растрескивание произошло в результате растяжения, то модуль упругости «бетонного» элемента при растяжении составляет 0.001 от исходного модуля в данном направлении - окружном или меридиональном. В случае превышения предела прочности на растяжение в обоих направлениях модуль упругости элемента задается близким к нулю в обоих направлениях. Однако растрескавшийся в результате растяжения элемент работает как обычный «бетонный» при смене знака напряжений, т.е. при сжатии этого элемента, что отражает процесс закрытия трещин при смене знака напряжений. В случае превышения сжимающими напряжениями предела прочности на сжатие происходит его разрушение и модуль упругости снижается до 0.2 от исходного значения. При этом, если элемент разрушился от сжатия то на растяжение он не работает. Т.к. модуль упругости элемента снижается не до нуля, то на последующих шагах нагружения при превышении окружными или осевыми напряжениями предела прочности на сжатие в этих элементах текущее значения модуля упругости в соответствующем направлении вновь умножается на 0.2, что соответствует умножению на 0.04 исходного значения модуля. Аналогичная процедура выполняется на всех последующих шагах расчета.

Если хотя бы в одном из элементов упругие свойства из-за образования тещин изменились, то выполняется следующая итерация: снова рассчитывается матрица жесткости и решается система конечноэлементных уравнений, рассчитываются значения напряжений в элементах. Если свойства не изменялись, то переходят к следующему шагу расчета. Поскольку расчеты выполнялись с применением суперэлементного алгоритма, то конечноэлементная модель представляется в виде совокупности блоков - суперэлементов. Это дает существенные вычислительные преимущества в случае локального растрескивания бетона, поскольку проводить вновь расчет матриц суперэлементов нужно не для всех, а только для отдельных суперэлементов.

На рисунках 9, 10 и 11 показана кинетика роста зон разрушенного бетона при воздействии пожара. Номера материалов, отличающиеся цветом, соответствуют следующим типам трещин:

- материал номер 11 - трещины в бетоне от окружных напряжений растяжения;

- материал номер 12 - трещины в бетоне от меридиональных напряжений растяжения;

- материал номер 13 - трещины в бетоне от окружных и меридиональных напряжений растяжения;

- материал номер 14 - трещины в бетоне от окружных напряжений сжатия;

- материал номер 15 - трещины в бетоне от меридиональных напряжений сжатия;

- материал номер 16 - трещины в бетоне от окружных и меридиональных напряжений сжатия;

- материал номер 17 - трещины в бетоне от окружных напряжений растяжения и меридиональных напряжений сжатия;

- материал номер 18 - трещины в бетоне от меридиональных напряжений растяжения и окружных напряжений сжатия.

На рис. 9 показана последовательность разрушения бетона в зоне локального разогрева купольной части на моменты времени 0.55, 1.30, 3.05, 4.90, 6.75, 8.60 часа. Справа на рисунках приводится таблица номеров материалов и соответствующих им цветов.

Представленная информация показывает, что в начальные моменты времени происходит образование трещин в бетоне от сжимающих окружных и меридиональных напряжений в поверхностном слое непосредственно в зоне разогрева, т.к. здесь из-за температурного расширения бетона и значительного градиента температуры реализуется высокий уровень сжимающих напряжений.

По мере прогрева стенки оболочки и смещении границы разогретого пятна вглубь наблюдается появление вблизи внутренней поверхности слоев бетона разрушившихся от воздействия растягивающих окружных и меридиональных компонентов напряжений, которые уравновешивают сжимающие напряжения, возникающие вблизи наружной поверхности непосредственно под зоной разогрева. Вблизи границ разогретого пятна реализуются зоны разрушения от отдельных сжимающих меридиональных или окружных компонентов напряжений, а также зоны смешанного разрушения - от сжимающих меридиональных и растягивающих окружных напряжений и наоборот.

На рис. 10 показана последовательность разрушения бетона в зоне локального разогрева цилиндрической части на моменты времени 1.20, 1.80, 2.70, 3.65, 5.50, 7.40 часа. На рис.11 представлена аналогичная информация на виде изнутри оболочки (облицовка удалена, чтобы показать «бетонные элементы) для моментов времени 1.80, 2.70, 3.65, 5.50, 7.40 часа. Кинетика роста разрушенных зон аналогична той, которая наблюдалась при разогреве купольной части. Различие заключается в том, что в этом случае растрескивание внутренних слоев бетона от напряжений растяжения на поздних стадиях разогрева связано появлением только кольцевых трещин от действия меридиональной компоненты напряжений. Это объясняется более высоким уровнем сжимающих окружных напряжений от преднапряжения, которые практически в два раза превышают меридиональную компоненту. Такое соотношение меридиональной и окружной компоненты напряжений от преднапряжений обусловлено необходимостью компенсации повышения давления под оболочкой при аварии. В то же время при локальном разогреве значения компонентов окружных и меридиональных напряжений возрастают в одинаковой степени. При этом растягивающие окружные компоненты компенсируются системой преднапряжения, а меридиональные - нет.

Таким образом, разработанные методики и программное обеспечение позволяют проводить анализ прочности защитной оболочки с учетом нелинейного поведения бетона с учетом образования трещин вследствие воздействия локальных полей температуры связанных с разогревом фрагментов внешней поверхности стенки.

Результаты исследований показывают, что для разработки рекомендаций по назначению мероприятий, снижающих уровень опасных последствий от воздействия повышенных и высоких температур, необходимо выполнить детальный анализ НДС защитной оболочки, а нормативный документ НП-010-98 «Правила устройства и эксплуатации локализующих систем безопасности атомных станций» не в полной мере отражает современные требования в части внешних воздействий и требует пересмотра.

В тоже время нормы МАГАТЕ по безопасности № NS-G-3.1 «Внешние события техногенного происхождения в оценке площадки для атомных электростанций» указывают на необходимость учета воздействий, вызванных горением топлива самолета (п. 5.16 - 5.19). Аналогичные требования указаны и в нормах МАГАТЭ № NS-G-1.10 «Design of Reactor Containment Systems for Nuclear Power Plants» и № NS-G-1.5 «External Events Excluding Earthquakes in the Design of Nuclear Power Plants» в случае, если расположение атомной станции требует рассмотрения воздействия авиационной аварии.

Пожары в результате таких сценариев, как авиационная авария, должны рассматриваться в той же комбинации нагрузки и с теми же проектными предположениями. Постулированные характеристики пожаров для моделирования могут быть описаны в терминах излучаемой энергии, площади возгорания и формы пламени, угловых факторов, скорости распространения и длительности воздействия.

Выводы

1. Разработанные методики и программное обеспечение позволяют проводить анализ прочности защитной оболочки с учетом нелинейного поведения бетона при образовании трещин вследствие воздействия локальных полей температуры.

2. Результаты расчета показали, что при воздействии высоких температур от горения топлива упавшего самолета защитная оболочка может потерять несущую способность.

3. В нормах МАГАТЕ по безопасности указывается на необходимость учета воздействий, вызванных горением топлива самолета.

4. Для учета влияния температурной нагрузки, вызванной горением топлива упавшего самолета на защитную оболочку АЭС, в Российские нормативные документы по безопасности необходимо внести соответствующие изменения.

атомный электростанция оболочка деформация

Список литературы

1. СНиП 2.03.04-84. Бетонные железобетонные конструкции, предназначенные для работы в условиях воздействия повышенных и высоких температур. НИИЖБ Госстороя СССР, 1985 г.

2. Рекомендации по расчету пределов огнестойкости бетонных и железобетонных конструкций. НИИЖБ Госстороя СССР, 1986 г.

Размещено на Allbest.ru

Рис. 9. Последовательность растрескивания бетона в купольной части при локальном разогреве на моменты времени: 0.55, 1.30, 3.05, 4.90, 6.75, 8.6 часа.

Рис. 10 Последовательность растрескивания бетона в цилиндрической части при локальном разогреве на моменты времени: 1.20, 1.80, 2.70, 3.65, 5.50, 7.40 часа.

Рис. 11. Последовательность растрескивания бетона в цилиндрической части при локальном разогреве на моменты времени: 1.20, 1.80, 2.70, 3.65, 7.40 часа (вид со стороны внутренней поверхности под облицовкой).