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

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

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

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

В настоящей работе представлена методика «Удар» для численного расчета напряженно-деформированного состояния железобетонных конструкций на заданный динамический импульс воздействия. Конструкция рассматривается в рамках трехмерного метода конечных элементов с явной схемой интегрирования по времени. Используются конечные элементы типа «гексаэдр» с восемью узловыми точками и восемью точками интегрирования по объему КЭ и типа «тетраэдр» с четырьмя узловыми точками и одной точкой интегрирования по объему КЭ.

Создание геометрии объекта, сетки армирования, генерация конечных элементов, задание нагрузки и граничных условий для численной модели производится при помощи программы FEMAP. Эта же программа используется для трехмерной визуализации результатов расчета после их выполнения в программе «Удар».

Математическая модель железобетона [1] учитывает неупругое деформирование бетона при сжатии и сдвиге, явное трещинообразование в зоне растяжения, дискретное расположение стержней армирования в бетоне.

Построение модели основано на принципе разделения механизмов динамического деформирования бетона.

Рассматриваются следующие механизмы деформирования [2]:

· механизм упругого деформирования бетона;

· механизм вязкопластического деформирования с накоплений пластических деформаций и повреждений бетона при сжатии и сдвиге;

· механизм неупругого уплотнения бетона от гидростатического сжатия;

· механизм накопления повреждений с образованием макротрещин при растяжении.

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

Рис. 1. Схемы исчерпания сопротивляемости бетонного элемента силовым воздействиям

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

Схема 2,3,4 (рис. 1б,в,г) - соответствует образованию от растягивающих напряжений одной (см. рис. 3.1а), двух (см. рис. 3.1б) или трех (см. рис. 3.1в) взаимно ортогональных трещин. После образования трещины бетонный элемент в перпендикулярных к площадкам направлениях не может воспринимать растягивающих напряжений. Однако возможна его работа на сжатия и сдвиг, с ограничением касательных напряжений.

В соответствии с принципом разделения процессов деформирования бетона полные деформации представляются в виде суммы следующих слагаемых:

 (1)

где - e - упругие деформации; p - пластические деформации;

k - компрессионные необратимые деформации, обусловленные закрытием микропор;

t - деформации трещинообразования, включающие в себя как необратимую часть (пластические) tp деформаций, связанных с процессами образования трещин, и обратимую часть te деформаций, обусловленных раскрытием, закрытием и относительным сдвигом берегов трещин.

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

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

Условие прочности в области одно-, двух-, трехосном сжатии описывается уравнением вида:

(2)

где - гидростатическое давления ;

Т - интенсивность касательных напряжений ;

J - второй инвариант дивиатора тензора напряжений

;

- угол вида напряженного состояния

;

J3 - третий инвариант дивиатора тензора напряжений

;

кривая поверхности прочности

.

Поверхность прочности в пространстве главных напряжений 1, 2, 3 (рис. 2) представляет собой форму усеченного конуса, сечение которой девиаторной плоскостью имеет вид треугольника со скругленными углами. Для наглядности ее удобно представить в координатах нормального касательного напряжений и угла вида напряженного состояния .

Рис. 2. К описанию поверхности прочности

Поверхность состоит из 6 одинаковых частей (лепестков) расположенных между =0 и =60 и симметрична относительно главных осей. В качестве направляющей выступает девиаторная кривая Fm (см. рис. 2б) которая является функцией Т от . Направляющими являются меридиональные кривые Т=F(0) . Каждому соответствует своя меридиональная кривая F. Линии F1(0) при =60 и F2(0) при =0 называются главными и являются базовыми кривыми для построения критерия прочности.

Кривая F1(0) проходит через реперную точку одноосносного сжатия () и аппроксимирует результаты трехосных испытаний бетона при , а кривая F2(0) проходит через точку двухосного равнокомпонентного сжатия () и аппроксимирует результаты трехосных испытаний бетона при . Девиаторной кривой обеспечивается согласование с экспериментальными данными для других видов напряженного состояния. В настоящее время исследователями предложены различные варианты для задания функций девиаторной кривой. Их анализ показывает что, наилучшего согласования с экспериментальными данными удается достичь при использовании критерия Вильяма-Варнке [4]:

, (3)

где параметр Надаи-Лоде связанный с углом вида напряженного состояния соотношением .

Выражения для меридиональных кривых прочности запишем виде:

, (4)

где R2 - относительная прочность бетона при двухосном равнокомпонентном сжатии; ai, bi, ci - безразмерные константы прочности, при этом задаются а1, с1, а2, с3, а остальные константы определяются исходя из непрерывности кривой F2 и ее первой производной, а также обязательного прохождения кривыми F1 и F2 через реперные точки (точку одноосного сжатия для - F1 , двухосном равнокомпонентного сжатия для F2).

Условие текучести с изотропным упрочнением и разупрочнением записывается в виде:

(5)

Статическая поверхность текучести подобна поверхности прочности и описывается соотношениями (2) и (3). При этом уравнения меридианных кривых и определяется следующим образом.

, (6)

где ai, bi, ci - константы прочности бетона, q - параметр упрочнения (разупрочнения) характеризующий уровень достигнутых напряжений по отношению к максимальной прочности, для одноосного сжатия q = /Rb.

Пластическое деформирование, имеющее место при f>0, сопровождается либо увеличением параметра q от qу (предел упругости) до qm = 1, либо снижением его значения от qm = 1 до qe = const (где qe - предел остаточной прочности).

Значения предела упругости qу, задающего начальное положение поверхности текучести, принято qу = 0,40,5, что физически приблизительно соответствует нижней границе микротрещинообразования в бетоне. При q = qm = 1 статическая поверхность текучести в области одно-, двух- и трехосном сжатии совпадает с поверхностью прочности. Параметр qе - определяет нижнее положение поверхности текучести на участке запредельного деформирования, значение. Отметим, что экспериментальные данные по запредельному деформированию носят противоречивый характер. В большинстве опубликованных работ на диаграммах деформирования - приводится нисходящая ветвь деформирования или показывается ее начальный участок. Хотя и существует понятие «остаточной прочности», физически соответствующее полному разрушению бетона и превращению его в несвязанную раздробленную среду, однако, в нормах стран ЕС и рекомендациях НИИЖБ принято за критерий разрушения бетона считать достижения деформациями бетона значений соответствующих уровню напряжений на нисходящей ветви диаграммы в диапазоне 0,30,75 от максимальной прочности.

Параметр упрочения (разупрочнения) принят как функция от некоторой скалярной величины , характеризующей уровень накопления пластических деформаций сдвига и принимающей значения (рис. 3):

0<c<1 - на участке упрочения;

1<c<2 - на запредельном участке деформирования.

Рис. 3. К описанию пластических деформаций и заданию параметра упрочнения

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

 (7)

где p - интенсивность скорости пластических деформаций сдвига:

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

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

Выбор функции q(c) основан на анализе экспериментальных диаграмм - при пропорциональном нагружении бетона в области одно-двухосного сжатия. Помимо экспериментальных данных использовались также рекомендованные европейскими нормами эталонные кривые для одноосного сжатия. Судя по рассматриваемым данным, функция q(c) на участке упрочнения существенно не линейна, а на участке запредельного деформирования близка к линейной зависимости. При этом непосредственный переход от восходящей ветви к запредельному деформированию должен быть плавным. Для получения гладких кривых () необходимо и достаточно, чтобы сама функция и ее первая производная были непрерывны, а на границах C = 0 и C = 1 производная принимала следующие значения:

при С=0 ; при С=1 . (8)

Этим требованиям отвечает функция вида (рис. 4):

(9)

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

Динамическое условие текучести описывается уравнением вида:

T - F(qdn,,)=0 (10)

Скорости пластических деформаций определяются следующим образом:

(11)

где - скорости интенсивности деформаций сдвига.

(12)

f =T - F(q,,)

Здесь отношение qdnq/q численно равно так называемому коэффициенту динамичности.

Ф - функция вязкозкопластичности, являющаяся обратной функцией коэффициента динамического упрочнения.

При численной реализации параметр qdn определяется следующим образом (рис. 4):

(13)

Рис. 4. К описанию динамического деформирования бетона

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

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

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

Описания данного процесса базируется на теории пластического течения и подходе Пэжины [6] к материалам чувствительным к скорости деформирования. В пространстве напряжений вводится плоскость текучести перпендикулярная гидростатической оси (рис. 6). Условие текучести записывается в виде

(14)

где Rk - константа модели, характеризующая начало неупругого уплотнения материала Rk=(0.5 -0.7)Rb.

Рис. 6. К описанию механизма неупругого объемного деформирования бетона

Функция упрочнения задает кривую напряжения - деформаций при гидростатическом сжатии (рис. 6)

Скорости необратимых компрессионных деформации вычисляются следующим образом:

(15)

Остановимся на описании механизма образования макротрещин от растяжения и деформирования бетона с трещинами. До начала нагружения в бетоне трещины отсутствуют, а деформации =0.

При действии растягивающих или смешанных напряжений типа «сжатие-растяжение» в каждом элементарном объеме (трехмерном КЭ, слое бетона пластины или оболочки) может образовываться семейство макротрещин трех взаимно перпендикулярных, но произвольно ориентированных в пространстве плоскостях.

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

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

Условие текучести (процесса образования трещин) при растяжении определяется неравенством:

, (16)

где Rbt - прочность бетона при одноосном растяжении;

- функция, задающая вид статической диаграммы деформирования при одноосном растяжении:

(17)

;

btm - деформации бетона, соответствующие концу участка текучести;

bte - предельное значение диаграммы, соответствующее образованию трещин;

qt(c) - функция, учитывающая снижение сопротивления бетона растягивающим усилиям вследствие напряжения повреждений (микротрещин) от сжимающих напряжений:

 (18)

Для определения скоростей необратимых деформаций воспользуемся вязкопластическим подходом [6].

(19)

Ниже, на рис. 7 показано сопоставление полученной в результате расчета кривой динамического упрочения при растяжении с зависимостью, предложенной в [5] Ю.В. Баженовым, аппроксимирующей экспериментальные данные.

Численные расчеты проводились при нагружении с постоянной скоростью деформирования .

Рис. 7. Сравнения результатов расчета коэффициента динамичности при растяжении с экспериментальной кривой

Из сопоставлений видно, что в диапазоне скоростей результаты хорошо соответствуют друг другу. Аналогичным образом при и осуществляется процесс образования трещин и накопления необратимых деформаций по двум другим направлениям и . Условия сдвига (проскальзывания) берегов макротрещины относительно друг друга определяется неравенством вида:

(20)

где - касательные напряжения в плоскости трещины;cy - прочность контактов берегов трещин.

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

Вклад усилий в узле арматуры i в значения узловых сил конечного элемента вычисляется, например для глобальной компоненты х, по формуле:

, (21)

где Fxi - усилие в i-ом узле арматурного стержня;

Nj - функции формы конечного элемента;

- локальные координаты узла i стержня арматуры в конечном элементе бетона;

k - количество узловых точек в конечном элементе бетона.

В качестве иллюстраций применения численной методики «Удар» ниже приведены ряд примеров расчета защитных конструкций АЭС при падении самолета.

На рис. 8 показана защитная оболочка АЭС, радиус которой составляет 25м, а толщина - 0,5м. Расчет производится на падение самолета легкого класса - «Lear Jet». На рис. 9 приведен график импульса воздействия. Пятно воздействия площадью 12м2 находится в зоне сопряжения сферической и цилиндрической частей защитного купола реакторного отделения.

Рис. 8. Расчет защитного купола АЭС на падение самолета Lear Jet

Рис. 9. График нагрузки в пятне воздействия при падении самолета типа «Lear Jet»

На рис. 10 представлен график изменения во времени прогиба оболочки в центре пятна нагружения. Максимальный прогиб был достигнут в момент времени t = 85мс - 48мм. Рис. 11 показывает изополя меридиональных деформаций оболочки при t = 85 мс. Видны зоны трещинообразования на внутренней части конструкции. Максимальные деформации растяжения составили 0,29%, а сжатия - 0,12 %.

Рис. 10. Изменение во времени прогиба оболочки при падении самолета Lear Jet

Рис. 11. Деформации в защитной оболочке при падении самолета Lear Jet (t = 85 мс)

На рис. 12 показан фрагмент большепролетной железобетонной плиты толщиной 1,5м с проемом. По оси y проходит плоскость симметрии, а на остальных внешних границах приняты условия защемления. Расчет производится на воздействие при падении самолета типа Phantom массой 20т (рис.13). Пятно воздействия площадью 7м2 находится вблизи проема.

Рис. 12. Фрагмент защитной конструкции с проемом

Рис. 13. Изменение во времени нагрузки при падении военного самолета Phantom RF-4E

Рис. 14 демонстрирует распределение по объему конструкции параметра накопления ущерба в бетоне С (см. формулу (7)). Превышение данного параметра величины 1 указывает на выход на ниспадающую ветвь деформирования, а превышение величины 2 свидетельствует о наличии зон разрушения бетона. В рассматриваемом случае мы наблюдем зоны сдвигового разрушения бетона в угловых зонах возле проема - областей концентрации напряжений в конструкции.

Рис. 14. Изополе параметра поврежденности бетона конструкции в момент максимального прогиба

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

1. Голяков В.И., Обизюк В.А. «Математическое моделирование деформирования бетона при интенсивных динамических нагрузках». Сборник научных трудов 26 ЦНИИ МО РФ, выпуск II/2, - М.: 26 ЦНИИ, 2006, с.56-77.

2. Голда Ю.Л. «Математическое моделирование процессов динамического деформирования железобетонных конструкций при действии интенсивных динамических нагрузок».//НТС «Некоторые проблемы механики инженерных сооружений и конструкций», - М.:26 ЦНИИ, 1998, с.126-186.

3. Берг О.Я. «Физические основы теории прочности бетона и железобетона». - М.: Стройиздат, 1961, 95с.

4. Willam K. J. and Warnke E. P. «Constitutive model for the triaxial behavior of concrete», Italy, Int. Assoc. Bridge Struct. Eng. Proc. Vol.19, pp.1-30, 1975.

5. Баженов Ю.М. «Бетон при динамическом нагружении». - М.: Изд. литературы по строительству, 1968, 271с.

6. Пэжина П. «Основные вопросы вязкопластичности». - М.: Мир, 1968, 176с.

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