Методика определения работоспособности энергонасыщенных материалов с использованием параметров вырожденных ударных волн

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

Рубрика Экономико-математическое моделирование
Вид статья
Язык русский
Дата добавления 31.08.2018
Размер файла 2,1 M

Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже

Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.

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

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

С.Ю. Ганигин, А.Ю. Мурзин, М.В. Ненашев, П.В. Письменный, И.Е. Хлыстова, А.А. Чеботаев Сергей Юрьевич Ганигин (к.т.н.), доцент, каф. технологии твердых химических веществ.

Андрей Юрьевич Мурзин (к.т.н.), доцент, каф. технологии твердых химических веществ.

Максим Владимирович Ненашев (д.т.н.), профессор, проректор по научной работе.

Павел Владимирович Письменный, старший преподаватель, каф. радиотехнических устройств.

Ирина Евгеньевна Хлыстова, аспирант, каф. технологии твердых химических веществ.

Александр Анатольевич Чеботаев, аспирант, каф. технологии твердых химических веществ.

Самарский государственный технический университет

443100, г. Самара, ул. Молодогвардейская, 244

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

Ключевые слова: работоспособность, энергонасыщенные материалы, акустический тротиловый эквивалент, ударные волны, звуковое давление.

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

Методика определения работоспособности энергонасыщенных материалов по параметрам вырожденных ударных волн

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

Как отмечено в [1], мерой относительной работоспособности ВВ могут быть также величины тротилового эквивалента, расчетные или полученные по экспериментальным измерениям параметров воздушных ударных волн (ВУВ). Такими параметрами могут быть импульс фазы сжатия волны или избыточное давление на ее фронте. В качестве импульсомеров используют маятниковые, поршневые и другие приборы. Избыточное давление измеряют напрямую или вычисляют по измеренной скорости ударных волн.

Сравнительная характеристика методов определения работоспособности энергонасыщенных материалов

Параметр

Метод

свинцовой бомбы

баллистической мортиры

баллистического маятника

по воронке выброса

Класс ВВ

ИВВ, БВВ

ИВВ, БВВ

ПВВ

ПВВ

Состояние испытуемых образцов

Порошкообразные, гранулированные вязко-текучие, жидкие, пластичные

Масса образцов

Не менее 10 г

Не менее 10 г

Не менее 200 г

Несколько кг

Состав технических средств

Свинцовая бомба Трауцля m=80 кг, средство инициирования (ОШ, ЭД), образец ВВ, бумажная гильза, кварцевый песок, средства измерения объема

Баллистическая мортира m=350 кг,

средство инициирования (ЭД и др.), образец ВВ, средства очистки и смазки рабочих поверхностей

мортиры

Баллистический маятник m=10 т, средство инициирования (ЭД и др.), образец ВВ, средства очистки и смазки рабочих поверхностей маятника

Средство инициирования (ЭД и др.),

образец ВВ

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

R* = R/G0,33, (1)

где R* - приведенное расстояние;

R - расстояние от эпицентра взрыва до заряда;

G - масса заряда.

На основе источников [2] следует, что многие эмпирические соотношения, связывающие давление во фронте ВУВ с приведенной массой или расстоянием, имеют вид

, (2)

где A, B, C, D - коэффициенты, найденные из опыта.

Подставляя выражение для приведенного расстояния, получаем:

, (3)

где МTNT - масса заряда тротила.

Приведем (3) к уравнению относительно массы заряда. Обозначим , тогда

. (4)

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

(5)

Таким образом, подставляя в (5) известное значение массы испытуемого заряда и найденное из (4) значение эквивалентного тротилового заряда, можно определить тротиловый эквивалент и принять его за меру работоспособности. Решение уравнения (4) может быть легко выполнено численными методами непосредственно в вычислительном устройстве системы измерения тротилового эквивалента.

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

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

Основными преимуществами предлагаемой методики определения работоспособности по сравнению со стандартными являются следующие:

1) определение работоспособности по предлагаемому методу осуществляется на образцах с массой в десятки и сотни раз меньшей, чем в стандартных методиках по ГОСТ 4546-81;

2) в процессе проведения испытаний не используются дорогостоящие материалы, как, например, цилиндрические блоки из свинца высокой чистоты массой 80 кг в методе свинцовой бомбы (повторно свинец не используется);

3) процесс инициирования и регистрации измерительной информации осуществляется автоматически и синхронно;

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

Качественно связь бризантного действия взрывчатых веществ с параметрами звуковых колебаний отражена в работе [3], где сопоставляли величины деформации оболочки зарядов образцов и мощность звуковых колебаний. Практический интерес представляет связь параметров акустических колебаний детонации в точках, расположенных на расстояниях, где скорость ВУВ становится равной скорости звука с работоспособностью бризантных ВВ. В экспериментах использовались вещества с разными значениями работоспособности, такие как ТЭН (работоспособность ДV = 500 см3), гексоген (ДV = 495 см3) и тетрил (ДV = 390 см3). Испытывались навески массой 50 мг. Инициирование осуществлялось сбросом груза массой 10 кг на копре К-44-2 в соответствии с методикой измерения чувствительности по ГОСТ 4545-88. Использовался электродинамический микрофон, установленный на расстоянии 3 м от образца. Зарегистрированные сигналы выравнивались по времени и определялись накопленные суммы квадратов отсчетов сигнала. Определялся акустический эквивалент работоспособности по параметрам взрыва опорного вещества (тротил) как

(6)

где sx - сигнал, полученный при «пробе» вещества с искомой работоспособностью;

s0 - сигнал, полученный при «пробе» вещества с эталонной работоспособностью (тротил).

Процедура обработки данных содержит механизм временного выравнивания реализаций сигналов. При этом за нулевой отсчет принимается момент времени, начиная с которого дисперсия реализации сигнала в окне из m отсчетов, соответствующем выбранному интервалу времени анализа, превышает предыдущее значение на заданную величину е1. Сигнал B(t) при одинаковых статистических характеристиках шума асимптотически стремится к некоторому фиксированному значению. Работоспособность (последний отсчет) определяется в точке, где дисперсия сигнала не превышает заданной величины е2, соответствующей дисперсии шума:

; . (7)

На рис. 1 представлен результат эксперимента по определению относительной работоспособности, нормированной по тротилу.

Р и с. 1. Результат обработки эксперимента по определению относительной работоспособности, нормированной по тротилу

Алгоритм расчета (6) - (8) позволяет использовать только нормированные значения работоспособности и требует необходимости предварительных испытаний образцов эталонных веществ, например тротила. Это является обычной практикой и для других методов: поскольку, например, величина расширения объема свинцовой бомбы и отклонения баллистического маятника, а также высота подъема мортиры не находятся в линейной зависимости от истинной работоспособности [4], то получаемые по всем методам результаты характеризуют собой лишь относительную работоспособность. Их обычно сопоставляют с работоспособностью кристаллического тротила, принятого за эталон.

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

. (8)

Приравнивая производную этого соотношения к нулю, находим точки экстремума:

, (9)

т. к. величина, обратная времени, не обращает это выражение в ноль нигде, кроме как на бесконечности, можно написать:

; . (10)

Окончательно значение работоспособности определяем как максимальное значение нормированной по времени удельной энергии (9)

, (11)

где i - номера точек экстремума;

N - общее количество точек экстремума на усредненной реализации сигналов, полученных в эксперименте.

Для реализации алгоритма (8-11) были проведены испытания в соответствии с разработанной методикой определения работоспособности ВВ с использованием копра К-44-2. Некоторые результаты испытаний приведены на рис. 2. Для примера на рис. 2 приведен ансамбль усредненных реализаций по 3-м опытам сигналов звукового давления, полученных при испытаниях образцов ТЭНа, тротила, гексогена, октогена и ТГ-50 массой 50 мг.

Р и с. 2. Ансамбль усредненных реализаций сигналов звукового давления, полученный при испытаниях образцов ТЭНа, гексогена, ТГ-50, октогена массой 50 мг

акустический тротиловый эквивалент

Использовался датчик звукового давления BF40 GRAS, установленный на расстоянии 1 м и направленный на образец. Ввиду особенностей работы емкостного микрофона реализации сигналов инвертированы. Максимальное зарегистрированное значение звукового давления не превышает 3000 Па или 165 дБ. Длительность первой положительной полуволны колебания не превышает 150 мкс. Длительность первой отрицательной полуволны не превышает 800 мкс. Реализации сигналов звукового давления приводились к нулевому уровню вычитанием постоянной составляющей. Расчет работоспособности осуществлялся на интервале 14 мс.

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

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

Временные диаграммы нормированной по времени удельной энергии сигнала на интервале 90 мс и 2 мс приведены на рис. 3, а, б.

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

Как видно из рис. 3, б, недостатком алгоритма расчета (8) - (11) является неоднозначность определения максимального значения этой характеристики, т. к. максимум может приходиться на первую, вторую или третью полуволну колебаний среды. Получаемые значения абсолютного эквивалента работоспособности имеют размерность [Па2]. Точность и повторяемость результатов измерений можно повысить, если работоспособность определять как отношение полученного значения к дополнительно измеренному акустическому импедансу среды. Вводимая таким образом величина получает размерность плотности потока мощности.

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

Р и с. 4. а - иллюстрация к определению значения работоспособности по звуковому эффекту взрыва образца ВВ (используются накопленные суммы удельной энергии); б - соотношение акустического эквивалента работоспособности ВВ и работоспособности, полученной методом свинцовой бомбы

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

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

ЗАКЛЮЧЕНИЕ

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

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1. Физика взрыва / Под ред. Л.П. Орленко. - Т.1. - М.: ФИЗМАТЛИТ, 2002. - 832 с.

2. Гельфанд Б.Е., Сильников М.В. Фугасные эффекты взрывов. - СПб.: Полигон, 2002. - 272 с.

3. Калашников В.В., Ненашев М.В., Керов А.В., Ибатуллин И.Д., Ганигин С.Ю., Мурзин А.Ю., Письменный П.В., Чеботаев А.А., Кобякина О.А., Хлыстова И.Е. Использование акустических колебаний для качественного и количественного анализа ударно-волновых процессов, порождаемых детонацией взрывчатых веществ, газовых смесей и аэрозолей // Известия Самарского научного центра РАН. Том 13, № 1. - С. 355-360.

4. Авакян Г.А., Дубнов Л.В., Мельников А.А., Ким Ю.М. Об экспериментальных и расчетных методах оценки работоспособности (фугасности) взрывчатых веществ // Взрывное дело. - М.: Недра, 1978. - № 80/37. - С. 22-29.

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

...

Подобные документы

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

    дипломная работа [238,1 K], добавлен 10.06.2011

  • Особенности расчета параметров уравнений линейной, степенной, полулогарифмической, обратной, гиперболической парной и экспоненциальной регрессии. Методика определения значимости уравнений регрессии. Идентификация и оценка параметров системы уравнений.

    контрольная работа [200,1 K], добавлен 21.08.2010

  • Основные параметры уравнения регрессии, оценка их параметров и значимость. Интервальная оценка для коэффициента корреляции. Анализ точности определения оценок коэффициентов регрессии. Показатели качества уравнения регрессии, прогнозирование данных.

    контрольная работа [222,5 K], добавлен 08.05.2014

  • Особенности группировки экономических данных. Методика определения средних показателей, мод, медиан, средней арифметической, индексов товарооборота, цен и объема реализации, абсолютных приростов, темпов роста и прироста. Анализ цен реализации товара.

    контрольная работа [51,1 K], добавлен 03.05.2010

  • Методика определения параметров линейной регрессии, составления экономической интерпретации коэффициентов регрессии. Проверка выполнения предпосылок МНК. Графическое представление физических и модельных значений. Нахождение коэффициентов детерминации.

    контрольная работа [218,0 K], добавлен 25.05.2009

  • Оценка параметров шестимерного нормального закона распределения с помощью векторов средних арифметических и среднеквадратического отклонений и матрицы парных коэффициентов корреляции (по программе Statistica). Методика определения Z-преобразования Фишера.

    контрольная работа [33,6 K], добавлен 13.09.2010

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

    лабораторная работа [35,7 K], добавлен 12.11.2014

  • Определение параметров линейной регрессии и корреляции с использованием формул и табличного процессора MS Excel. Методика расчета показателей парной нелинейной регрессии и корреляции. Вычисление значений линейных коэффициентов множественной детерминации.

    контрольная работа [110,4 K], добавлен 28.07.2012

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

    лабораторная работа [2,1 M], добавлен 01.12.2009

  • Особенности определения пространственных и временных данных на примере диаграмм структуры использования денежных доходов. Понятие парной регрессии, порядок ее решения. Методика составления матрицы переходных вероятностей для средних годовых изменений.

    контрольная работа [62,9 K], добавлен 10.09.2010

  • Методика определения максимального числа подзадач без пределов и максимального числа целочисленных решений без пределов. Метод "затраты-эффект" и особенности его применения в управлении проектами. Задача привлечения дополнительных финансовых ресурсов.

    контрольная работа [37,9 K], добавлен 09.07.2014

  • Общая характеристика и порядок определения коэффициента корреляции, методика и этапы его оценки. Описание автокорреляционных функций. Сущность критерия Дарбина-Уотсона. Примеры практических расчетов с помощью макроса Excel "Автокорреляционная функция".

    курсовая работа [1,7 M], добавлен 03.07.2010

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

    презентация [197,0 K], добавлен 02.05.2014

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

    презентация [107,3 K], добавлен 22.08.2015

  • Использование метода оценки параметров в стандартных масштабах для определения неизвестных параметров линейной модели множественной регрессии. Специфика изучения взаимосвязей по временным рядам. Моделирование взаимосвязей и тенденций в финансовой сфере.

    контрольная работа [326,7 K], добавлен 22.04.2016

  • Построение модели управления запасами в условиях детерминированного спроса. Методы и приемы определения оптимальных партий поставки для однопродуктовых и многопродуктовых моделей. Определение оптимальных параметров системы управления движением запасов.

    реферат [64,5 K], добавлен 11.02.2011

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

    дипломная работа [372,5 K], добавлен 29.11.2014

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

    дипломная работа [1,7 M], добавлен 19.03.2014

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

    курсовая работа [230,6 K], добавлен 18.04.2014

  • Решение задач линейного программирования с применением алгоритма графического определения показателей и значений, с использованием симплекс-метода. Использование аппарата теории двойственности для экономико-математического анализа оптимального плана ЗЛП.

    контрольная работа [94,6 K], добавлен 23.04.2013

Работы в архивах красиво оформлены согласно требованиям ВУЗов и содержат рисунки, диаграммы, формулы и т.д.
PPT, PPTX и PDF-файлы представлены только в архивах.
Рекомендуем скачать работу.