_ ПРОБЛЕМЫ ЯДЕРНОЙ, РАДИАЦИОННОЙ _
И ЭКОЛОГИЧЕСКОЙ БЕЗОПАСНОСТИ -
УДК 539.2/.6:623.45:662.22
ОПРЕДЕЛЕНИЕ УДАРНЫХ УСКОРЕНИЙ НА ЭЛЕМЕНТАХ
КОНСТРУКЦИИ МАЛОГАБАРИТНОЙ ЯДЕРНОЙ ЭНЕРГЕТИЧЕСКОЙ УСТАНОВКИ ПРИ ЕЕ ПАДЕНИИ НА
ПОВЕРХНОСТЬ
© 2020 О.А. Губеладзе, А.Р. Губеладзе
Донской государственный технический университет (ДГТУ), Ростов-на-Дону, Ростовская обл., Россия
В статье проводится исследование одного из этапов развития аварийной ситуации с межконтинентальной баллистической ракетой, сопровождающейся взрывом носителя, выбросом головной части из шахты и последующим ее падением. Выполнено моделирование соударения головной части со слабодеформируемой преградой при падении, проведен анализ результатов расчета и эксперимента.
Ключевые слова: межконтинентальная баллистическая ракета, шахтная пусковая установка, головная часть, слабодеформируемая преграда, ударное ускорение.
Поступила в редакцию 29.11.2019 После доработки 03.02.2020 Принята к публикации 10.02.2020
Основной целью обеспечения безопасности на этапе эксплуатации ядерно- и радиационно опасных объектов (ЯРОО) является принятие эффективных мер, направленных на предотвращение тяжелых аварийных ситуаций (АС). Определение слабых мест этого объекта и формулировка основных принципиальных решений по повышению безопасности осуществляется на основе результатов вероятностного анализа для ЯРОО. Перечень АС, составленный с учетом видов опасностей ядерного оружия (ЯО) [1], значения характеристик внешних воздействующих факторов в моделях аварийных ситуаций определяются с учетом опыта эксплуатации ЯРОО, а также выполненных исследований (экспериментальных и теоретических работ) их безопасности [2]. К настоящему времени широко известно о многочисленных авариях и катастрофах, связанных с ядерным оружием, среди которых особый интерес представляет связаная с взрывом компонентов топлива ракеты с ядерной головной частью (ГЧ). Рост температуры и давления в стволе шахты [3, 4] могут привести к выбросу ГЧ с последующим ее падением.
В 1980 году в шахте взорвалась межконтинентальная баллистическая ракета (МБР) «Титан-2» (США) [5]. Утром 19 сентября 1980 года произошёл взрыв компонентов топлива ускорителя первой ступени МБР, ударная волна вырвала закрытую крышу массой 740 тонн железобетонного защитного устройства шахты (рис. 1). Ускоритель второй ступени вместе со штатной ГЧ Мк.6 с ядерным боеприпасом (ЯБП) W53 были выброшены из шахты. В результате взрыва второй ступени ГЧ (масса - 3690 кг) была подброшена в воздух на 200 метров и упала в 30 метрах от строительного сооружения [6]. Следует отметить, что в СССР также происходили взрывы жидкостных ракет в шахтных пусковых установках, но во всех случаях носители были без ГЧ, которые на время проведения регламентных работ (технического обслуживания) отстыковывались [7, 8]. В подавляющем большинстве
© Национальный исследовательский ядерный университет «МИФИ», 2020
АС, связанных с взрывом носителей ЯО в шахтных пусковых установках, инициирующим событием являлись ошибочные действия персонала.
а) б)
Рисунок 1 - Стартовая позиция межконтинентальной ракеты «Титан-2» (США): а) - до аварии; б) - после взрыва МБР [The starting position of the Intercontinental missile «Titan-2» (USA) a) - before the
accident; b) - after the explosion of ICBM]
Для исследования соударения ЯРОО со слабодеформируемой преградой при падении из перечня возможных этапов развития АС [9] выбраны представленные в таблице 1.
Таблица 1 - Этапы развития аварийной ситуации [The stages of emergency situation development]
№ этапа Этап Перечень параметров вторичных воздействий
Э-6 Взрыв ракетного топлива температура, давление; продолжительность воздействия; скорость выброса ГЧ из шахты ПУ
Э-2 Удар ГЧ о поверхность (преграду) скорость удара; характеристики преград
Э-13 Взрыв, сгорание
Проведем моделирование взаимодействия (соударения) объекта со слабодеформируемой преградой (поверхностью) при его падении. Используем математическую модель, приведенную в монографии [10], проведя ее уточнение. Рассмотрим свободное падение объекта по нормали к поверхности. При решении контактной динамической задачи по определению сил, действующих в зоне контакта объекта и преграды, исходим из допущения, что кинетическая энергия падающего объекта в основном перейдет в тепловую энергию удара и потенциальную энергию деформации преграды (условие не разрушения конструкции при падении). Потенциальной энергией упругой деформации конструкции пренебрегаем. В этом случае перегрузка в центре масс (координата Х) изделия получится несколько завышенной. На рисунке 2б также представлены: I - общая длина объекта и Я - радиус притупления носка. Также примем допущения: форма той части поверхности объекта, которая встречается с преградой, близка к сферической; ось симметрии объекта при соударении нормальна к поверхности преграды; радиус притупления равен Х; движение в области преграды происходит в вертикальной плоскости; скорость встречи объекта с преградой существенно меньше скорости распространения упругих волн в
области преграды. Расчетная схема, выбранная с учетом принятых допущений, представлена на рисунке 2в.
л
J
Рисунок 2 - Объект исследования: а) - общий вид; б) - основные геометрические параметры; в) - расчетная схема [Object of research a) - general view; b) - basic geometric parameters;
c) - calculation scheme]
При составлении уравнения движения объекта в области преграды для случая, когда вектор скорости направлен перпендикулярно к плоскости поверхности падения, предположим, что вертикальная составляющая реактивного давления со стороны преграды на поверхность контакта постоянна и равна q, причем интенсивность этого давления зависит от скорости движения изделия в области преграды. Уравнение равновесия сил в направлении вертикальной оси изделия будет иметь вид:
т ■ — + О = 0, (1)
йХ
где т, G - масса и сила тяжести объекта;
V - скорость движения в области преграды; N (у) = жг2 д, - равнодействующая сила реакции преграды; г - радиус лунки.
Если глубина лунки к по сравнению с радиусом поверхности контакта объекта небольшая, то:
г2 = 2ВИ . (2)
Удельное давление q будет пропорционально второй степени скорости движения V:
д(у )= дс
( У
V
V у0 J
(3)
Здесь q0 - наибольшее удельное давление со стороны преграды при вдавливании в него круглого штампа; Vo - начальная скорость входа изделия в преграду. Тогда
т--ъ 2лКд0
йХ
к - О = 0.
V у0 J
(4)
йк
Из соотношения — = V йХ
к = | уйХ + С
(5)
Подставляя (5) в (4) (с учетом, что в начальный момент времени к = 0, то и С = 0), получим формулу движения объекта в области преграды:
т--ъ 2лкq0
йХ
(
V V J 0
| vdХ - О = 0.
(6)
Уравнение (3) решается приближенно [11,12] с учетом
2 Л
v = у0 соб
2Хп
(7)
где Х0 - время движения объекта в области преграды до остановки. Тогда ускорение:
йу лТуп . ЛХ
а = — = --
-Бт-
йХ 2Х0 2Х0
(8)
Интеграл в выражении (6):
[уйХ = v гсоб2 ^—йХ = — V. 1 01 2Х0 2 '
( Х0 . ЛХЛ Х Ъ—°81П —
V Л Х 0 J
(9)
Выражения (7), (8) и (9) подставим в (6) и решим, используя способ ортогонализации:
лти, . Л -Б1П--ълKq0v0
(
2Х„
Хп . 711 Х Ъ—°Б1П — Л Х
\
0 J
4 Л
соб--О
2Х„
2 Л
соб -йХ = 0.
2Х„
Выражение для определения времени Х0 будет иметь вид:
0
2
V
0
0
о,71эа
'о = -
Яд о Уо
1 +
1 +
4Яд о Уо2 ^
(11)
где g - ускорение свободного падения.
Максимальное значение перегрузки в центре масс объекта
а„
п =
шах
жЯд 0 V2
g
1,43Gg
1 +
1 + 4Яг^
1,43Gg
2 Л'
(12)
Если пренебречь в уравнении (6) влиянием силы тяжести и принять G = 0, то оно примет вид:
ёу 2жЯд о 2, ^
т— +-у к = о.
Л у2
(13)
Проведем замену переменных, тогда у--1-
ёк тУ2
ёу 2жЯд„ 2, ё
■ У2 к = о или —
Л,2 Л
ёк
V 2 у
+ ау2 к = о,
где а =
2лЯд0
ту,
2 '
Обозначив 2 = —, запишем 2
ё ёк
/„2 Л
V 2 у
0 ёг + ау2 к = — + 2агк = о. ёк
(14)
Разделим переменные и, интегрируя, получим:
г = С ■ ехр (- ак 2 ), (15)
где С - постоянная интегрирования. Если к=0, то г = у2, тогда г = у^ ■ ехр (— ак2).
ёу
Значение перегрузки в центре масс объекта будет п =а = , но
£ £
ëv ëv ё1
— = V--= —, тогда:
ё1 ёк ёк
ёу
п = а = А. = — к ехр (— ак2) . £ £ £
(16)
ёп
Наибольшее значение перегрузки находим из условия — = о. Глубина лунки,
ёк
соответствующая ёп ау!
максимальному
значению
перегрузки
ёк £
ехр (— ак2 )2ак2 — 1)= о,
0
2
к = (17)
л/2а
Тогда:
п_ =-1,075^ ^. (18)
£ V т
Из формулы (18) видно, что величина перегрузки зависит от радиуса оболочки изделия, его массы и механических свойств преграды. Интенсивность удельного давления со стороны преграды на проникающую лобовую часть изделия и величина радиуса лунки при наибольшей перегрузке, будут равны:
д = д0
с у
V
V ^ J
= д0 exp(-ак2)=—(19)
= 42Як =
exp (0,5)'
тЯ
Л?0
(20)
Для приближенных прочностных расчетов характеристики различных преград можно получить следующим образом, так как поверхность контакта объекта имеет форму, близкую к сферической, то осадка поверхности преграды при проникновении в нее круглого штампа определяется по формулам:
- в центре лунки к = 1,13(1 - ¡2 (21)
- на границе к = 0,72(1 - ¡л2 (22)
Здесь Е - площадь поперечного сечения штампа;
Л - коэффициент Пуассона материала преграды; Е - модуль упругости преграды.
Из формул (21) и (22) можно получить усредненное значение осадки преграды:
К = 1 (к + к )= 0,925(1 -л2 УЁ^. (23)
2 Е
Формула (23) определяет взаимосвязь между физико-механическими характеристиками грунта и глубиной лунки. С учетом (17) получим выражение для определения д :
0,643Е 4л[а
20 =706-^=. (24)
(1 -л2}1ЛЁа '
Подставляя (24) в (18) и используя v0=gХ, для определения максимального значения перегрузки получим следующее выражение:
= -0,904
(1 -*2)
(25)
m
g
где Н - высота падения объекта.
С целью проверки адекватности расчетной модели была проведена серия экспериментов с использованием стенда [13]. На рисунке 3 представлен массово-габаритный макет (МГМ), масса которого могла изменяться (т = 247^980 кг), со сменным носком (Я = 0,1; 0,3; 0,6 м). На рисунке 4 показан фрагмент экспериментального стенда (вертикальный копр с тросом подвески и замком сброса макета), позволяющий обеспечить высоту сброса МГМ в диапазоне 0,28 ^ 4,30 м на поверхность преграды (бетон, мелкий песок и твердая глина). В испытаниях применялся измерительный комплекс, в который входили: система измерения и регистрации ударных ускорений, состоящая из датчиков ускорения, устанавливаемых на элементах конструкции МГМ, и регистратора; скоростная киносъемка сброса, приземления и отскока МГМ.
Рисунок 3 - Массово-габаритный макет: а) - компоновочная схема; б) - сменный носок; в) - днище [Mass-dimensional layout а) - Layout diagram; b) - replaceable toe; c) - bottom]
Рисунок 4 - Фрагмент экспериментального стенда [Fragment of the experimental stand]
Результаты экспериментальных исследований и расчетов по формуле (25) для объекта т = 400 кг и Я = 0,1; 0,3 м представлены на рисунке 5.
Hiri* 1000
750
500
БпФп м- 0,17: Е- 4,:-10J
\ ; : ' 1 1 » * : : i--1---1-
■ J* 1 * jr- , 2
rji Л * ■ / * w * p 4
1 * . * I 4 :
а)
Иных
500 400
300 200 100
ПзИХ 100
75
50
25
А Н, м
Глина И-0,41: Е-0,76-10*
JltiOK II- 0.2 P: Еы 0.0M0'
f k « it ! i * V
i-vi2
i £ * * * #
* * 4
; ; . : : x 1 * ' У» * 1 / т : !/ i
A J ™ /V ' / X i / У i / / 1 i / / 1 1
V ' ' / k 1 о
н;м
л н,м
^ 0 эксперимент
_ 1 -ш= 400 кг, Е.= 0,3 ы
_ 2 - ш= 400 кг, R. = 0,1 ы
в)
Рисунок 5 - Результаты моделирования падения объекта: а) - на бетон; б) - на мелкий песок; в) - на твердую глину [The results of falling object simulation: a) - concrete; b) - sand; с) - clay solid]
Анализ полученных результатов позволил сделать следующие выводы:
1. Расхождение расчетных значений с результатами эксперимента составило: для бетона - до 17%; для песчаной преграды - не более 10%. Большее расхождение для бетонной поверхности обусловлено принятым допущением об условии не разрушения конструкции при падении.
2. Для одной и той же преграды при наименьшей глубине проникновения отмечается наибольшее значение перегрузки в центре масс объекта.
3. Увеличение радиуса носка притупления (при m = const) приводит к росту значений перегрузки.
4. При достижении ударных ускорений на элементах МГМ 160 ед. в зонах крепления грузов отмечались значительные механические повреждения.
5. Результаты расчетов для высот падения 120^200 метров позволяют прогнозировать полное разрушение конструкции объекта с возможным радиоактивным заражением местности осколками заряда. Однако аварийный взрыв ЯБП маловероятен [5, 10].
СПИСОК ЛИТЕРАТУРЫ
1. Губеладзе, О. А. Экспресс-оценка результатов нерегламентированных деструктивных воздействий на ядерно- и радиационноопасный объект / О. А. Губеладзе, А. Р. Губеладзе // Глобальная ядерная безопасность. - 2018. - № 4. - С. 24-30.
2. Губеладзе, О. А. Оценка последствий воздействия кинетических ударников на взрывоопасный объект с малогабаритной ядерной энергетической установкой / О. А. Губеладзе, А. Р. Губеладзе // Глобальная ядерная безопасность. - 2019. - № 3. -С. 33-40.
3. Энергетические конденсированные системы. Краткий энциклопедический словарь / Под редакцией Б. П. Жукова. - Москва : Янус К, 2000. - 483 с.
4. Орленко, Л. П. Физика взрыва и удара / Л. П. Орленко. - Москва : ФИЗМАТЛИТ, 2008. - 304 с.
5. Кириллов, В. М. Физические основы радиационной и ядерной безопасности / В. М. Кириллов. - Москва : РВСН, 1992. - 212 с.
6. Cochran, Thomas B.; Arkin, William M.; Hoenig, Milton M. Nuclear Weapons Databook, Vol. II: U.S. Nuclear Warhead Production. Cambridge, Mass.: Ballinger Publishing Company, 1987. 240 p.
7. Гуляев, А. А. Герои и жертвы «холодной войны» / А. А. Гуляев // Военно-исторический журнал. - 2014. - № 12. - С.47-54.
8. Малиновский, Г. Н. Записки ракетчика / Г. Н. Малиновский. - Москва : ЦИПК РВСН, 1999. -207 с.
9. Михайлов, В. Н. Безопасность ядерного оружия России / Под редакцией В. Н. Михайлова. -Москва : Министерство по атомной энергии, 1998. - 148 с.
10. Денисов, О. В. Комплексная безопасность населения и территорий в чрезвычайных ситуациях. Проблемы и решения : монография / О. В. Денисов, О. А. Губеладзе, Б. Ч. Месхи, Ю. И. Булыгин; под общей редакцией Ю. И. Булыгина. - Ростов-на-Дону : Издательский центр ДГТУ, 2016. - 278 с.
11. Вайнберг, М. М. Интегро-дифференциальные уравнения / М. М. Вайнберг // Итоги науки. Сер. Мат. анал. Теор. вероятн. Регулир. - 1962, ВИНИТИ, Москва, 1964. - С. 5-37.
12. Busa J., Hnati M., Honkonen J., Lucivjansky T. Numerical Solution of a Nonlinear Integro-Differential Equation. Mathematical Modeling and Computational Physics. Published at: EPJ Web of Conferences 108, 02017 (2016).
13. ГОСТ Р ИСО 8568-2010 Национальный стандарт РФ. «Стенды ударные» (введен в действие Приказом Федерального агентства по техническому регулированию и метрологии № 441-ст от 12.11. 2010 г.). - URL : https://internet-law.ru/gosts/gost/50555/ (дата обращения: 17.07.2019).
REFERENCES
[1] Gubeladze O.A. Express-otsenka rezul'tatov nereglamentirovannykh destruktivnykh vozdeystviy na yaderno- i radiatsionnoopasnyy ob'ekt [Express Assessment of Results of Independent Destructive Impacts on Nuclear and Radiation-Hazardous Object] Globalnaya yadernaya bezopasnost [Global Nuclear Safety]. 2018. №4 (29). Р. 24-30 (in Russian).
[2] Gubeladze O.A. Otsenka posledstvij vozdeystviya kineticheskih udarnikov na vzryvoopasnyj ob'ekt s malogabaritnoj yadernoj energeticheskoj ustanovkoj [Assessment of the Influence Consequence of Kinetic Impactors at Explosive Object with Small-Scale Nuclear Power Facility] GlobaFnaya yadernaya bezopasnost [Global Nuclear Safety]. 2019. №3 (32). Р. 33-40 (in Russian).
[3] Energeticheskie kondensirovannye sistemy. Kratkiy entsiklopedicheskiy slovar' [The Power Condensed Systems. Short Encyclopedic Dictionary]. Moskva [Moscow]. Yanus K. 2000. 483 p. (in Russian).
[4] Orlenko L.P. Fizika vzryva i udara [Physics of Explosion and Blow]. Moskva. FIZMATLIT [Moscow. PHIZMATLIT]. 2008. 304 p. (in Russian).
[5] Kirillov V.M. Fizicheskie osnovy radiatsionnoy i yadernoy bezopasnosti [Physical Bases of Radiation and Nuclear Safety]. Moskva [Moscow]. RVSN. 1992. 212 p. (in Russian).
[6] Cochran, Thomas B.; Arkin, William M.; Hoenig, Milton M. Nuclear Weapons Databook, Vol. II: U.S. Nuclear Warhead Production. Cambridge, Mass.: Ballinger Publishing Company, 1987. 240 p.
[7] Gulyaev A.A. Geroi i gertvy «holodnoj vojny» [Heroes and Victims of «the Cold War»] Voyenno-istorichesky zhurnal [Military Historical Journal]. 2014. № 12. P.47-54 (in Russian).
[8] Malinovskiy G.N. Zapiski raketchika [Notes of a Rocketeer]. Moskva [Moscow]. CIPK RVSN. 1999. 207 p. (in Russian).
[9] Mikhaylov V.N. Bezopasnost' yadernogo oruzhiya Rossii [Safety of Nuclear Weapon of Russia]. Moskva. Ministerstvo po atomnoy energii [Moscow. Ministry of Nuclear Energy]. 1998. 148 p. (in Russian).
[10] Denisov O.V., Gubeladze O.A., Meskhi B.Ch., Bulygin Yu.I. Kompleksnaya bezopasnost' naseleniya i territory v chrezvy'chajny'x situaciyax. Problemy' i resheniya: monografiya [Complex Safety of the Population and Territories in Emergency Situations. Problems and Solutions]. Rostov-na-Donu. Izdatel'skij centr DGTU [Rostov-on-Don. Publishing center Don State Technical University]. 2016. 278 p. (in Russian).
[11] Vainberg M.M. Integro-differencial'ny'e uravneniya [Integro-Differential Equations]. Itogi Nauki. Ser. Mat. Anal. Teor. Ver. Regulir [The Results of Science. Series Mathematical Analysis Probability Theory Regulation]. 1962. Moskva [Moscow]. VINITI, 1964, pp. 5-37 (in Russian).
[12] Busa J., Hnati M., Honkonen J., Lucivjansky T. Numerical Solution of a Nonlinear Integro-Differential Equation. Mathematical Modeling and Computational Physics. Published at: EPJ Web of Conferences 108, 02017 (2016).
[13] GOST R ISO 8568-2010 Nazionalnyj standart RF. «Stendy udarnye» [National Standard of the Russian Federation. Shock Testing Machines. Declaration and Verification of Characteristics]. URL : https://internet-law.ru/gosts/gost/50555/ (in Russian).
The Determination of Shock Accelerations on a Small Nuclear Power Plant Structural Elements When Falling to the Surface
12 O.A. Gubeladze , A.R. Gubeladze
Don State Technical University, Gagarin square 1, Rostov-on-Don, Russia, 344000 'ORCID iD: 0000-0001-6018-4989 WoS Researcher ID: F-6921-20' e-mail: buba2602'966@yandex.ru 2ORCID iD: 0000-0002-6966-639' WoS Researcher ID: F-72'5-20'7 e-mail: buba2602'966@yandex.ru
Abstract - The article conducts research one of the expansion phase of emergency with an Intercontinental Ballistic Missile, accompanied by an carrier explosion, an ejection of the reentry vehicle from a silo and its subsequent fall. The simulation of the collision of the reentry vehicle with a weakly deformable barrier in the event of fall is performed, the result calculation valuation and experiment are carried out.
Keywords: intercontinental ballistic missile, silo launcher, head part, weakly deformable barrier, shock acceleration.