УДК 621.7.044.2
ОПРЕДЕЛЕНИЕ ОПТИМАЛЬНОЙ МАССЫ ОБОЛОЧЕК ВЗРЫВНЫХ КАМЕР
© 2011 Ю.В. Скорняков, Ю.П. Мещеряков
Конструкторско-технологический филиал Института гидродинамики им. М.А. Лаврентьева СО РАН, г. Новосибирск
Поступила в редакцию 05.11.2011
Изучен вопрос оптимизации массы оболочки взрывной камеры. В работе рассматривались два расчетных случая - идеальная сферическая и цилиндрическая оболочки с постоянными толщинами стенок. Получены аналитические зависимости, позволяющие связать параметры заряда взрывчатого вещества (ВВ) с необходимыми прочностными свойствами материала оболочки и её геометрическими размерами. При некоторых допущениях определены максимальные и минимальные теоретически возможные значения массы оболочки для указанных случаев. Впервые получена аналитическая зависимость радиуса оболочки от радиуса заряда ВВ, которая позволяет на начальных этапах проектирования производить оценку необходимой по условию прочности геометрии камеры и определить её оптимальную массу.
Ключевые слова: взрывные камеры, сферическая оболочка, цилиндрическая оболочка, оптимальная масса
Металлические взрывные камеры являются уникальным оборудованием, позволяющим локализовать энергию взрывчатого вещества в замкнутом пространстве и использовать её в различных целях, в том числе технологических. Их использование позволяет нейтрализовать практически все виды воздействия взрыва на человека и окружающую среду - ударную волну, сейсмическое воздействие, разлет осколков и продуктов детонации ВВ. Оптимизация массы и габаритов оболочек необходима для создания более эргономичных взрывных камер и расширения области их применения, что является важной актуальной задачей в настоящее время. В специализированной литературе имеются публикации по вопросам оптимизации оболочек, но в них, как правило, применяется анализ на основе метода конечных элементов и громоздких численных расчетов и потребность в создании простых и эффективных методик оценочного расчета оболочек взрывных камер на начальных этапах проектирования остается. Некоторые рекомендации по данному вопросу приведены в работе [1], но они также нуждаются в уточнении и дополнении. В данной работе выведены простые аналитические зависимости, с помощью которых можно решать задачи оптимизации массы и размеров оболочек взрывных камер непосредственно в процессе проектирования.
Постановка задачи. Рассмотрим два расчетных случая - сферическая и цилиндрическая оболочки, внутри которых расположен заряд
Скорняков Юрий Викторович, аспирант. E-mail: scorni-akov@mail.ru
Мещеряков Юрий Павлович, кандидат физико-математических наук, ведущий научный сотрудник группы численного моделирования. E-mail: ura@kti-git.nsc.ru
взрывчатого вещества. Схемы случаев и необходимые параметры приведены на рис. 1 и 2. Сферическая оболочка (рис. 1) имеет радиус RC и толщину 5с, в её центре расположен сферический заряд взрывчатого вещества произвольной массы Mс, которое имеет удельную теплоту взрыва Q0. Цилиндрическая обечайка бесконечной длины (рис. 2) имеет радиус Rц, толщину 5ц и высоту H, на её оси симметрии расположен заряд ВВ произвольной массы Mц (выражается через г0ц и Н) который имеет удельную теплоту взрыва Q0. В расчетах будем считать Н=1 м, т.е. будем рассматривать погонные величины.
ВУ V
ВВ (HQ J
б
Рис. 1. Параметры сферической оболочки
Для обоих случаев определим оптимальную с точки зрения экономии массы толщину стенки оболочки при условии обеспечения заданной прочности по допускаемым напряжениям.
Будем считать, что в сферической оболочке реализуется двухосное напряженное состояние, напряжение ог примем равным нулю (оно действует вдоль радиуса). В цилиндрической оболочке напряженное состояние одноосное - напряжение ог и о2 примем равными нулю. В ходе расчета величины в формулах для сферической оболочки будут содержать индексы «с», а для цилиндрической - «ц».
Рис. 2. Параметры цилиндрической оболочки
Аналитическое решение. Итак, согласно [1] имеем:
2
Ро (Г0С )3 Ел/^
GM =
З-у/Г—Й р8с (Rc )2
Ро (оц ))у12О0Е р5ц ■ ^ а о
(1)
Здесь: р0 и г0 - соответственно плотность и радиус заряда ВВ. Е - модуль упругости. р - плотность материала оболочки. а0 - скорость звука в материале оболочки. Формулы (1) справедливы для тонкостенных оболочек (Я>105). Масса заряда ВВ
Mc = 4 пРо (гос )3
M Ц = пРо (г0Ц )3 H
(2)
=°ec =
i _ McEJQQ
2пл/1 - Й P8c (Rc )2 ао
=
M ц EfiQo
(3)
ттрдц ■ ^ ао Эквивалентные напряжения
J^p )2 + (°ec )2
-ачс°с =°С
Для обеспечения прочности оболочки, эквивалентные напряжения не должны превышать допустимых значений стэкв < [г]. За допустимые
значения примем предел ^упругости при циклическом нагружении (N>10 циклов). В результате получим следующие зависимости
1 М с Е^
м=
2^ 1 - Й Р8с (Rc )2
[ст]= M ц EjlQo пр8ц ■ ^ао
(4)
Введем некоторые коэффициенты К1 (сферическая оболочка) и К2 (цилиндрическая оболочка)
K1 =
1
__McEJQ
2^ 1 - Й РаоИ
M ц
K 2 =
пр■ а,
• аоИ
которые определяют физические свойства материала оболочки, а также свойства и параметры заряда ВВ. Формулы (4) примут следующий вид
¿c =
8 ц =
K1
R)
K 2
(5)
Rr,
Н - высота цилиндрической оболочки (а также, соответственно, и высота заряда). Будем считать Н=1 м, т.е. далее во всех формулах для цилиндра Мц будет являться величиной погонной массы заряда (размерность кг/м) и выражаться формулой Мц = пр0 (г0ц ) . С учетом (2), формулы (1) примут вид
Формулы (5) определяют необходимую толщину стенки, обеспечивающую заданную прочность в зависимости от радиуса оболочки. Выразим массу сферической и цилиндрической оболочки через известные размеры
mc =
4 пр[(с + 8с )3 -( )3
тц =
пр■ H [R +8ц )2-(ц )2
(6)
Высоту Н также примем равной единице. Тогда тц - погонная масса цилиндрической обечайки с
о
°рС= иС =
о
размерностью кг/м, тЦ = пр- [((Ц +3Ц )2 - ЯЦ 2 ]. Подставив в формулу (6) зависимость (5), полу-
чим
тС = —пр
тц = пр-
Re +-^1 -(Re
(Re
(7)
R„ +
K 2
Rrr
-(ц )2
Формулы (7) определяют необходимые для обеспечения заданной прочности массы оболочек при заданных радиусах и массах зарядов. Аналогичную зависимость вида т = m(S) можно получить, выразив из (5) R=R(S) и подставив в (6). Получим следующие выражения
mC = — пр
тц = пр -
K1 .
г i—л3
^ ' Kl ^
(8)
( K 2 Л '2
Vo Ц
( K 2 ^2
V Ц У
Формулы (8) определяют необходимые для обеспечения заданной прочности массы оболочек при заданных толщинах и массах зарядов.
Таблица 1. Механические и физические характеристики материалов оболочки камеры
Название Символ Величина
коэффициент Пуассона м 0,3
плотность р 7800 кг/м3
скорость звука в металле оболочки ао 5000 м/с
модуль упругости Е 2,1 ■1011Па
Анализ и выводы. Определим диапазон значений коэффициентов К1 и К2.
Kl = f (и, Р, ао, И E, M, Q0) K 2 = f (и, р, а0, И E, M, Q0)
Взрывные камеры изготавливаются из хорошо свариваемых, малоуглеродистых конструкционных сталей, а также нержавеющих сталей.
Величины, приведенные в таблице 1, меняются незначительно и могут считаться константами для всех марок сталей, используемых в производстве ВК. Принимая во внимание, что энергия высвобождающаяся при взрыве зарядов различных ВВ, определяется их тротиловым эквивалентом, в дальнейшем без ограничения общности можно считать удельную теплоту взрывчатого материала постоянной величиной равной Qo=4,2 106 Дж/кг, что соответствует
взрывчатому веществу тротилу. С учетом вышесказанного, получим
K1 = 2,099 -106 -
M_
И
6 M
K2 = 4,967 -106 --¡^
И
(9)
Анализ формул (9) показывает, что зависимости К1(М) и К2(М) при постоянном значении [о] линейны. Зависимости К1([о]) и К2([о]) при постоянном М- обратно пропорциональны.
Рассмотрим более детально зависимости т(М) и т([о]). Примем величину допускаемых напряжений постоянной и равной [о]=235 МПа (соответствует стали 09Г2С широко применяемой в производстве ВК). Теперь для различных значений массы заряда М определим коэффициенты К1 и К2 (табл. 2).
Таблица 2. Значения коэффициентов К1 и К2 для различных по массе М зарядов ВВ
М, Масса заряда ВВ, кг К1 К2
0,1 0,0007 0,0021
0,25 0,0019 0,0053
0,5 0,0037 0,0106
1 0,0075 0,0211
2 0,0150 0,0423
Подставив полученные значения в выражение (7), получим следующие графические зависимости для сферической оболочек (рис. 3). Графики функций т(Я) для цилиндрической оболочки будут подобными изображенным на рис. 3. и в данной работе не приводятся.
2000 1 750 ■ 1500 ■ 1 250 ■ j 1000 ■ с \ р=10б: 1 М=2
\
1
I
....... М=2
750 ■ 500 ■ 250 ■ •
/ М~ 0,5
М=0,25
\ ч.\ М=0,2
0,2 0,4 0,6 0,8 1 1,2 1,4 R, м
Рис. 3. Зависимость массы т сферической оболочки от радиуса Я при разных значениях массы М (в кг) заряда ВВ
Выводы (они правомерны и для сферической и для цилиндрической оболочек):
4
2
4
1. Для заданной массы заряда, геометрии оболочки и свойств материалов, из которых она сделана, существует минимально возможная по условию прочности масса отличная от нуля. Она достигается при Я—да (5—>0)
m„
= lim тС (R) = 4пр ■ K1
= lim тц( R) = 2пр- K 2
min Ц
R~ (10)
2. Линия R=10^ соответствует границе применимости формулы (1). В заштрихованной области оболочку следует считать толстостенной и вводить соответствующие поправки в расчетные формулы. Кроме того, в оболочке с геометрическими параметрами, лежащими вне области допустимых значений, возможно появление локальных пластических напряжений.
3. В обратном случае при R^-0 и ô^œ масса будет возрастать. Теоретическим пределом, ограничивающим сверху величину массы оболочки, будет случай R=r0. Т.е. заряд охватывается снаружи оболочкой, без зазора
mmax = т(Г0).
4. Очевидно, что при увеличении массы M заряда ВВ возрастает масса оболочки т. Как видно по рисунку, это отражается в смещении графика функции m(R) вверх и вправо при возрастании М.
5. При увеличении значения допустимого напряжения, график функции m(R) смещается вниз и влево. Подробнее об этом будет сказано в дальнейшем.
6. В [1] приведена оценка для массы оболочки в зависимости от массы заряда ВВ, из которой было определено, что на 1 кг массы заряда должно приходится 1000 кг металла оболочки. Это утверждение может быть уточнено. Во-первых, изменение размера оболочки приводит к изменению величины отношения m/M, что видно из графиков (рис. 3). Во-вторых, как уже отмечалось ранее, для каждого коэффициента К1 и К2 существует определенная минимальная масса mmin. Отношение mmin/M есть величина постоянная и её можно использовать для первичной оценки массы оболочки проектируемой камеры. Но для фиксированной величины радиуса Ri отношение m(Ri)/M будет меняться при изменении М.
7. Черным точкам на графиках функции соответствуют значения "¿¡Rm(R) =1. В них скорость
изменения массы совпадает со скоростью изменения размера оболочки. Левее, варьирование размера оболочки приводит к более резкому изменению её массы. Правее же, наоборот - масса меняется медленней, чем размер оболочки. Фактически этим точкам соответствует величина теоретического оптимума. Так как теоретически оптимальные соотношения не лежат в области
применимости исходных расчетных формул, следует считать, что оптимальными являются точки Я=10^, как наиболее близкие к теоретическим.
Для получения более удобных критериев оптимальности оболочек, выразим отношения масс оболочек к их минимальным возможным значениям в зависимости от соотношения Я к г0. С этой целью построим графики зависимости т(Я) в следующей системе координат. По оси абсцисс отложим величины, кратные г0 (т.е. Я будет выражен через радиус заряда ВВ - г0). Значениям на оси ординат будут соответствовать величины относительной массы вида
тг
(r ) = -
' V0C J
m
mC (r0C )
m
ОТНЦ
(r0 Ц )
min C (r0C ) mц (r0 Ц )
m
min Ц V 0Ц
(V )
Здесь r0 - радиус заряда ВВ, определяемый по формулам
r0C = 3
3Mr
4ПР
r =
' 0 Ц
M
Пр0
Тогда графики, изображенные на рис. 3, примут следующий вид - рис. 4, линия 2. Соответствующие графики функций для цилиндрической оболочки показаны на рис. 5.
1 1 . 1 Г г
\
\ \
ц
Rfr„
Рис. 4. Зависимость относительной массы сферической оболочки от её относительного радиуса:
1 - [с]=279 МПа, 2 - [с]=235 МПа, 3 - [с]=150 МПа Выводы:
1. Для указанных расчетных случаев (М=0,1-2 кг, [о]=235 МПа), все пять графиков функций т(Я) (рис. 3.) одинаковы в этих осях координат и соответствуют кривой 2 (рис. 4). Можно полагать, что получена универсальная графическая функция, связывающая массу оболочки с её размерами, которую можно применять для любых значений
массы заряда ВВ. Это утверждение также справедливо для цилиндрической оболочки (рис. 5.).
l : V V \ .Ш s;; \
——
О 10 20 30 <Ю 50 60 70 90 90 100 110 120 130
я/У,
Рис. 5. Зависимость относительной массы цилиндрической оболочки от её относительного радиуса: 1 - [с]=279 МПа, 2 - [с]=235 МПа, 3 - [о]=150 МПа
2. Квадратными маркерами на графиках обозначены границы применимости формул (1), выше их (графики показаны штриховыми линиями) оболочки нужно считать толстостенными и вносить соответствующие поправки в расчет.
3. Случаю теоретического оптимума соответствует величина равная примерно 1,1 для сферы. Это означает, что самая компактная оболочка (размером 8-10 г0) будет тяжелее самой легкой теоретически возможной (с массой mmin) лишь примерно на 10%. Это справедливо для любой массы заряда ВВ. В случае цилиндрической оболочки теоретический оптимум составит примерно 1,05 и максимальная разница при указанных условиях составит примерно 5%.
4. При увеличении значения допустимого напряжения общая масса оболочки уменьшается.
По графикам также видно, что использование более высококачественных марок сталей позволяет получать более оптимальные оболочки с массой, более близкой к минимальной (при заданном радиусе оболочки). Это позволяет получить двойной выигрыш при использовании качественных сталей.
5. Для величин Л>15г0 масса сферической оболочки меняется незначительно (не более чем на 1% при дальнейшем увеличении размеров). Следовательно, при проектировании сферических оболочек ВК целесообразно придерживаться величин 8-10 г0<К<15 г0, так как этим значениям соответствуют более компактные и оптимальные по массе сферические оболочки. В случае цилиндрических оболочек для величин R>60-70 г0 масса оболочки меняется незначительно (не более чем на 1% при дальнейшем увеличении размеров).
Следовательно, при проектировании цилиндрических оболочек ВК целесообразно придерживаться величин 30-40 г0<К<60-70 г0, так как этим значениям соответствуют более компактные и оптимальные по массе цилиндрические оболочки.
СПИСОК ЛИТЕРАТУРЫ:
1. Думчук, А.Ф. Металлические взрывные камеры: монография / А.Ф. Демчук, В.П. Исаков. - Красноярск: РИО КрасГУ, 2006. 297 с.
2. Феодосьев, В.И. Сопротивление материалов / В.И. Феодосьев. - М.: Наука, 1970. 544 с.
3. Работнов, Ю.Н. Механика деформируемого твердого тела / Ю.Н. Работнов. - М.: Наука, 1979. 744 с.
DEFINITION OF MASS OPTIMIZATION OF BLASTING CHAMBER SHELLS
© 2011 Yu.V. Skornyakov, Yu.P. Meshcheryakov
Design and Technological Branch of Institute of Hydrodynamics named after M.A. Lavrentyev SB RAS, Novosibirsk
Has arrived in edition 05.11.2011
The question of mass optimization of blasting chamber shell is studied. Two settlement cases - ideal spherical and cylindrical shells with constant thickness of walls were in-process observed. The analytical dependences are gained, allowing to connect explosive charge (EC) parameters with necessary of strength properties of shell material and its geometrical sizes. At some assumptions the maximum and minimum theoretically possible mass values of a shell for the specified cases are defined. For the first time analytical dependence of shell radius on EC radius which allows to make an estimation of chamber geometry concerning strength conditions at the designing initial stages and to define its optimum mass is gained.
Key words: blasting chambers, spherical shell, cylindrical shell, optimum mass
Yuriy Skornyakov, Post-graduate Student. E-mail: scorniakov@mail.ru Yuriy Meshcheryakov, Candidate of Physics and Mathematics, Leading Research Fellow at the Numeric Modeling Group. E-mail: ura@kti-git.nsc.ru