close

Вход

Забыли?

вход по аккаунту

?

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

код для вставкиСкачать
ФИО соискателя: Реут Игорь Игоревич Шифр научной специальности: 01.04.17 - химическая физика, горение и взрыв, физика экстремальных состояний вещества Шифр диссертационного совета: Д 212.217.01 Название организации: Самарский государственный техниче
На правах рукописи
Реут Игорь Игоревич
ВЛИЯНИЕ ПАРАМЕТРОВ ДЕТОНАЦИИ КОНДЕНСИРОВАННЫХ
ВЗРЫВЧАТЫХ ВЕЩЕСТВ И ДИНАМИЧЕСКИХ СВОЙСТВ МЕТАЛЛОВ НА
ЭФФЕКТИВНОСТЬ МЕТАТЕЛЬНОГО ДЕЙСТВИЯ
01.04.17 – Химическая физика, горение и взрыв,
физика экстремальных состояний вещества
Автореферат
диссертации на соискание ученой степени
кандидата технических наук
Самара – 2012
Работа выполнена в Федеральном государственном бюджетном
образовательном учреждении высшего профессионального образования
«Самарский государственный технический университет»
Научный руководитель: д.т.н., доцент Кривченко Александр Львович.
Официальные оппоненты:
Епифанов Владимир Борисович д.т.н., доцент, ФГБОУ ВПО «Самарский государственный технический университет», профессор кафедры Химия и технология полимерных и композиционных материалов;
Смирнов Александр Сергеевич д.т.н., ОАО «Государственный научноисследовательский институт машиностроения имени Бахирева», начальник отдела экспериментальной физики.
Ведущая организация Учреждение Российской академии наук «Институт структурной макрокинетики и проблем материаловедения РАН» (г. Черноголовка).
Защита диссертации состоится «16» мая 2012 г. в 14 часов 00 мин. на заседании диссертационного совета Д 212.217.01 ФГБОУ ВПО «Самарский государственный технический университет» по адресу: 443100 г. Самара, ул. Молодогвардейская, 244, Главный
корпус, аудитория 500.
Отзывы по данной работе в двух экземплярах, заверенные печатью, просим направлять
по адресу: Россия, 443100 г. Самара, ул. Молодогвардейская, 244, Главный корпус на
имя ученого секретаря диссертационного совета Д 212.217.01; факс: (846) 242-28-89.
С диссертацией можно ознакомиться в библиотеке ФГБОУ ВПО «Самарский государственный технический университет» (ул. Первомайская, 18).
Автореферат разослан «____» апреля 2012 г.
Ученый секретарь
диссертационного совета
Самборук Анатолий Романович
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Актуальность работы.
Сложные динамические взаимодействия представляют самостоятельный интерес
для различных областей науки и техники, в том числе и для разработки новых систем
динамического оружия, у которого поражающие элементы баллистической формы
должны иметь достаточно высокие характеристики по форме и массе, которые определяются соотношением массы и квадрата скорости. D-металлы и их сплавы являются
наиболее перспективными в этом плане, так как, во-первых, неполно изучены их ударно-волновые свойства, а во-вторых, в данное время в термодинамике рассмотрены только полуэмпирические модели уравнения состояния металлов, что требует дополнительного проведения дорогостоящих экспериментов по определению динамических характеристик металлов и их сплавов. Однако именно в этой концепции возможно решение
проблем разработки эффективных поражающих элементов кинетического оружия, зарядов взрывного бурения и кумулятивных перфораторов.
Для использования на практике эффектов динамического нагружения необходимо
знать такие характеристики как скорость звука во взрывчатом веществе и в метаемом
материале, детонационные и массовые скорости и ряд других характеристик. Задачи
разработки методов теоретического прогнозирования скоростей детонации и метания
требуют комплексной диагностики ударно-волновых процессов при динамическом
нагружении. Прогнозирование таких процессов является актуальной научной задачей.
Исследование поведения металлических материалов при ударном сжатии имеют большое значение для решения задач по прогнозированию взрывных воздействий на различные материалы.
В работе также уделено внимание решению актуальной задачи по теоретическому
исследованию параметров детонации и метательной способности некоторых перспективных взрывчатых веществ на основе 1,3,5-триазинов.
Цель работы.
1. Комплексное исследование взаимодействия продуктов детонации (ПД) с метаемой
пластиной;
2. Разработка методов экспериментального и теоретического прогнозирования метательной способности;
3. Определение влияния материала метаемого элемента на эффективность отбора энергии им у продуктов детонации, а также выбор критерия его оценки.
Задачи исследования.
1. Определение параметров взаимодействия продуктов детонации с метаемой пластиной и влияния динамических характеристик металла на параметры его нагружения;
2. Исследование аддитивных свойств d-элементов таблицы Д.И. Менделеева с целью
определения их динамических адиабат;
3. Разработка теоретических методов определения скорости торцевого метания и расширения цилиндрической оболочки;
4. Определение детонационных и метательных свойств производных 1,3,5-триазина.
Методы исследования.
В работе использован аппарат ударно-волнового взаимодействия ПД и нагружаемых материалов, теоретический расчѐт параметров детонации, детонационной оптики.
Ударно-волновые импульсы меди и латуни, в условиях ударного нагружения,
изучались с помощью электромагнитного метода регистрации скоростей метания исследуемых образцов в процессе ударного сжатия. Расчѐт параметров исследуемых взрывчатых веществ и металлов осуществлялся в рамках полуэмпирического подхода.
1
Объекты исследования – прессованные заряды взрывчатых веществ: A-IX-1 и ТГ50/50;
металлические пластины: медь М1 и латунь Л75.
Научная новизна работы.
– Предложены методы прогнозирования параметров метания гидродинамическим методом и методом эквивалентных масс;
– Предложен способ по усреднению результатов, который позволяет уменьшить абсолютную ошибку по отношению к экспериментальным данным;
– Разработана предварительная экспресс-методика определения смесевых динамических адиабат сплавов, которая даѐт возможность определить наиболее оптимальный для
метания сплав;
– Впервые определены детонационные параметры взрывчатых веществ ряда 1,3,5триазинов и их возможная метательная способность.
Практическая значимость работы определяется предложениями по уменьшению количества экспериментальных исследований за счѐт научно-обоснованного подхода к
расчѐтно-экспериментальной оценке характеристик эффективности.
– Разработан метод теоретического прогнозирования определения скоростей метания;
– Разработана новая полуэмпирическая методика, позволяющая определить ударные
адиабаты металлов и сплавов для поражающих элементов;
– Разработанная экспресс-методика определения смесевых динамических адиабат
сплавов позволяет расширить круг различных металлов и сплавов и выделить из них
наиболее оптимальные для необходимых целей;
– Предложен перспективный экспериментальный метод регистрации ударноволновых импульсов взаимодействия продуктов детонации с метаемым элементом при
динамическом нагружении;
– Разработанные в работе методы приведут к снижению затрат на разработку боеприпасов с использованием перспективных взрывчатых веществ и оптимальных материалов
облицовок для снарядоформирующих зарядов.
На защиту выносятся следующие положения:
– Разработка электромагнитного метода фиксации скорости метаемого элемента и характер взаимодействия продуктов детонации с материалом пластины;
– Разработка принципа расчѐта скорости метания гидродинамическим методом;
– Разработка способа расчѐта коэффициента λ ударной адиабаты D=α+λU на основе
аддитивных свойств элементов периодической системы Д.И. Менделеева для определения динамических адиабат d-металлов и их сплавов;
– Способы расчѐта параметров торцевого метания гидродинамическим методом и методом детонационной оптики, и скорости расширения цилиндрической оболочки с использованием единой энергетической характеристики – массовой скорости;
– Результаты теоретического расчѐта параметров детонации и скоростей метания для
перспективных взрывчатых веществ производных 1,3,5-триазинов.
Достоверность научных результатов работы подтверждается совпадением результатов эксперимента с результатами теоретического прогнозирования по разным методикам и экспериментальным данным, а также выполненными оценками погрешностей измерений.
Апробация работы.
Научные результаты работы апробированы на Международных, Всероссийских и
Межвузовских научно-технических конференциях и симпозиумах: VII Всероссийской
научно-технической конференции «Наука. Промышленность. Оборона» (2006 г.), г. Новосибирск; XXXII Самарской областной студенческой научно-технической конферен-
2
ции памяти профессора Х.С. Хазанова. «Общественные, естественные и техн. науки»
(2006 г.), г. Самара; XVII Менделеевской конференции молодых учѐных (2007 г.), г. Самара; Международной конференции «Ударные волны в конденсированных средах»
(2008 г.), г. Санкт-Петербург; XI, XIII Харитоновских тематических научных чтениях –
Международная конференция «Экстремальные состояния вещества. Детонация. Ударные волны» (2009 г., 2011 г.), г. Саров; I Международной научно-практической конференции «Наука и современность» (2010 г.), г. Новосибирск; Международной конференции «Ударные волны в конденсированных средах» (2010 г.), г. Великий Новгород; Московский семинар по физике взрыва (2011 г.), г. Дзержинск.
Публикации.
По результатам диссертации опубликовано 12 работ, докладов и тезисов конференций, в том числе три статьи в изданиях, рекомендованных Высшей аттестационной
комиссией РФ.
Личный вклад автора.
В выборе направления исследований; в постановке задач диссертации; личном
проведении экспериментов; автор принимал участие в количественной обработке полученных экспериментальных данных; формулировке выводов и положений, выносимых
на защиту, а также в подготовке публикаций в печать.
Структура и объем работы.
Диссертационная работа состоит из введения, четырѐх глав, выводов, библиографического списка. Общий объем работы составляет 117 страниц, включая 16 рисунков,
25 таблиц и библиографического списка включающего 101 наименование.
ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ
Во введении дана общая характеристика диссертационной работы, обоснована
актуальность изучаемой проблемы, сформулированы цели работы, перечислены полученные в диссертации новые результаты, их практическая ценность, изложено основное
содержание работы и изложены основные положения, выносимые на защиту.
В первой главе анализируется современное состояние проблемы. Представлен
анализ способов определения метательной способности заряда взрывчатого вещества.
Рассмотрены практические и теоретические методы определения скорости метаемых
элементов. Рассмотрен способ определения оптимальной толщины пластины для торцевого метания. Показано, что все приведѐнные методы основаны на едином принципе, в
котором основным параметром является скорость звука в веществе.
На основании проведенного обзора сформулированы цели и задачи исследований.
Во второй главе разрабатываются способы теоретического определения метательной способности взрывчатых веществ (ВВ). Показано, что для гидродинамического
метода определения скорости метания необходимым является знание динамических
адиабат продуктов детонации (ПД) заряда взрывчатого вещества и ударных адиабат метаемых пластин. Для прогнозирования динамических адиабат ПД требуется определение параметров детонации ВВ. По результатам анализа известных экспресс-методик
сделан выбор направлений разработки для новых методов расчѐта скорости метания, а
также сформулированы научные принципы разработки методов расчѐтной оценки основных свойств ВВ и металлов. Анализ 10 экспресс-методик показал, что для расчета
параметров детонации ВВ оптимально использовать выражения:
для скорости детонации D:
1,65
(1)
D 1,2C0 o 2,55 k
oQm (1 Д ) , м/с
5,5
для массовой скорости продуктов детонации U:
3
U
k 1,65
5,5
oQm (1 Д ) , м/с
(2)
где C0 – объѐмная скорость звука во взрывчатом веществе, м/с; ρ – плотность заряда,
г/см3; ρ0 – предельная плотность, г/см3; αk – кислородный коэффициент; Qm – максимальная теплота взрыва, Дж/кг; βД - массовая доля добавки.
Для вычисления объѐмной скорости звука во ВВ, значение которой входит в выражение (1), проведѐн анализ методов расчѐта объѐмной скорости звука в органических
веществах. В качестве оптимального способа расчѐта выбран метод Рао, который применяется как для твѐрдых, так и для жидких органических веществ и имеет наиболее
близкое совпадение данных расчѐта и эксперимента
3
Ra С0 ,
Mм (3)
где Ra – инкремент химической связи данного вида; ρ – плотность, г/см3; Mм – молекулярная масса вещества, по которому распространяется звуковая волна.
Эксперименты показывают, что для различных материалов зависимость между
скоростью фронта и скоростью вещества за фронтом является линейной и форма задания ударной адиабаты в D-U-координатах, в которых кривая может быть приближѐнно
аппроксимирована линейной функцией:
D=α+λU
(4)
где D – скорость ударной волны; U – массовая скорость; α и λ – эмпирические коэффициенты, причѐм константа α численно равна объѐмной скорости звука C0.
Следует также отметить, что объемная и продольная скорости звука являются
чувствительным индикатором для фазового состояния ударно сжатых веществ. Знание
скоростей звука в материалах требуется также для правильной постановки динамических экспериментов в условиях, которые позволяют исключать влияние волн разгрузки.
Для определения значений эмпирических коэффициентов α и λ уравнения ударной адиабаты D=α+λU нами разработано «Правило симметричного окружения», применяемое для периодического закона Д.И. Менделеева. Расчѐт коэффициента для элемента
проводится по выражению:
Y
1 8
Yi ,
8 i 1
(5)
где Y – эмпирический коэффициент ударной адиабаты (4).
На рис. 1 представлен пример схемы определения коэффициента λ в уравнении
ударной адиабаты по «правилу ближайшего окружения».
24
Cr
1,47
42
Mo
1,24
74
W
1,29
Rh 25
Mn
26
Fe
43
Tc
44
Ru
75
Re
76
Os
27
Co
1,33
45
Rh
1,64
77
Ir
1,94
28
Ni
29
Cu
46
Pd
47
Ag
78
Pt
79
Au
Cr Co Zn Cd Hg Ir W Mo
30
Zn
1,56
48
Cd
1,67
80
Hg
2,20
1,59
8
Рис. 1. Схема расчета коэффициента λ по «правилу симметричного окружения» на примере элемента №45 Rh.
4
Разработанный способ определения скорости звука в металлах, а также обоснованно отобранные методы вычисления параметров детонации и скорости звука для ВВ,
легли в основу определения метательной способности ВВ, а также нахождению ударных
адиабат для металлов и сплавов.
При разработке теоретического метода оценки метательной способности индивидуальных ВВ по методике Т-20 предполагалось, что кинетическая энергия ускоряемого тела пропорциональна объѐмному энергосодержанию ВВ. Известно, что скорость разлѐта оболочки должна зависеть от динамических характеристик ВВ и массовой
скорости ПД, т.е. скорости звука в том или ином виде, в соотношениях, определяющих
свойства метаемой оболочки α/λ, а также коэффициента нагрузки β в соотношении
β=mВВ/mоб (где mВВ – масса ВВ; mоб – масса оболочки, приходящаяся на 1 см3 заряда ВВ).
Для разработки метода оценки скорости расширения оболочки был использован
массив данных, состоящий из 25 ВВ. Нами предлагается выражение по определению
скорости расширения цилиндрической оболочки:
C U , м/с
(6)
Wcyl ВВ
2,6 где CВВ – скорость звука во ВВ, определяемая по методу Рао (3); U – массовая скорость
ПД, определяемая по выражению (2); β – коэффициент нагрузки; значение 2,6 численно
связано с коэффициентами α и λ уравнения ударной адиабаты (4) для Cu (ρ=8,93 г/см3),
т.е. α/λ = 3,915/1,495 = 2,6.
В табл. 1 представлены результаты вычисления Wcyl и экспериментальные данные, описанные в литературе, а также приведены средние значения для двух методик.
Для уменьшения процента абсолютной ошибки был использован метод усреднения расчѐтных значений полученных по:
1) интерполяционному уравнению А.С. Смирнова:
WТ-20 = 20,50ρз0,418Qк0,458Nг0,245βВВ0,310Δr0,1
(7)
3
где WТ-20 – скорость движения оболочки, м/с; ρз – плотность заряда ВВ, г/см ; Qк – теплота взрыва (калориметрическая), ккал/кг; Nг – число молей газообразных продуктов
взрыва, моль/кг; βВВ – коэффициент нагрузки; Δr = R-R0 , мм (R, R0– текущий и начальный радиусы расширяющейся трубки). Коэффициент множественной корреляции
R=0,986, среднеквадратическая ошибка Δ=±52 м/с (±3,29%);
2) разработанному приближѐнному методу определения скорости расширения цилиндрической оболочки Wcyl.
Таблица 1
Экспериментальные и расчѐтные данные скорости расширения медной цилиндрической оболочки по методике Т-20
№
п/п
Скорость, м/с
расч.
расч.
WТ-20
Wcyl
5
6
ρз,
г/см3
эксп.
2
3
4
1.
Октоген
1,73
1590
1598
1610
1604
+0,9
2.
Октоген
1,81
1650
1722
1661
1691
+2,5
3.
Октоген
1,91
1783
1812
1716
1764
-1,1
4.
ТНТ
1,63
1180
1270
1182
1226
+3,9
5.
ТНТ
1,66
1230
1287
1231
1259
+2,3
6.
CL-20
2,04
2005
2120
2040
2080
+3,7
1
ВВ
5
Усреднение
7
%
ош.
8
1
2
3
4
5
6
7
8
7.
Бензотрифуроксан
1,80
1580
1673
1527
1600
+1,3
8.
1,80
1573
1611
1536
1574
+0,1
1,78
1510
1545
1454
1499
-1,7
10.
Гексоген
Тринитроэтиловый эфир тринитромасляной кислоты
Бис(тринитроэтил)нитрамин
1,90
1570
1623
1492
1557
-0,8
11.
ТАТБ
1,80
1290
1330
1346
1338
+3,7
12.
Тетрил
1,70
1305
1394
1344
1369
+4,9
13.
Гексанитростильбен
1,65
1240
1320
1227
1273
+2,7
14.
Гексанитробензол
1,78
1703
1610
1627
1618
-4,9
15.
Пентаэритриттетранитрат
1,76
1560
1585
1492
1538
-1,4
16.
1,77
1567
1593
1504
1548
-1,2
1,70
1590
1682
1511
1596
+0,4
18.
Пентаэритриттетранитрат
N,N-бис(тринитроэтил)этилендинитроамин
Диэтанолнитраминдинитрат
1,67
1418
1470
1481
1475
+4,0
19.
Бис(нитраминометил)нитрамин
1,69
1485
1530
1509
1519
+2,2
20.
Бис(нитраминометил)нитрамин
1,85
1675
1717
1610
1663
-0,7
21.
Бис(нитраминометил)нитрамин
1,89
1783
1760
1699
1729
-3,1
22.
N,N'-динитропиперазин
1,63
1292
1280
1373
1326
+2,7
23.
5-оксо-3-нитро-1,2,4-триазол
1,85
1330
1401
1262
1331
+0,1
24.
Нитрометан
1,14
1092
1081
1016
1048
-4,0
25.
Метилнитрат
1,21
1133
1217
1074
1145
+1,1
9.
17.
Среднеквадратическая ошибка для 25 ВВ по уравнению (6) составляет 3,56%.
Способ усреднения приводит к уменьшению процента ошибки, при этом среднеквадратическая ошибка составила 2,55%. Таким образом, усреднение результатов расчѐта по
двум методикам увеличивает достоверность прогнозирования скорости метания.
При разработке методики расчѐта относительной метательной способности
взрывчатых веществ при торцевом метании по методике М-60 использовалась модель «активной массы» ВВ. Нами предложен метод эквивалентных масс, в котором энергетика ВВ, участвующего
в разгоне пластины, сопоставляется с метаемым элементом в соотношении 1:1 (см. Рис. 2). В качестве критерия
принимается эквивалентное взаимодействие ВВ с метаемым элементом.
Известно, что при постоянной высоте ВВ h еѐ энер- Рис. 2. Схема методики эквивалентных масс
гия E, участвующая в метании, определяется соотношением:
(8)
E mU 2 ,
где U – массовая скорость продуктов детонации; m – активная масса ВВ.
С учѐтом бокового и осевого разлѐта ВВ активная масса, участвующая в разгоне
пластины, будет занимать объѐм конуса с радиусом основания r и высотой h:
6
2
m r 30 ,
3
(9)
2
E r 30U 2 .
3
(10)
где ρ0 – плотность заряда ВВ.
Подставив (9) в (8) получим
При этом метаемая пластина радиуса r с толщиной hпл, получившая энергию
Eпл mплW 2
,
2
(11)
имеет массу:
mпл r 2hпл пл ,
где ρпл – плотность метаемой пластины.
Тогда выражение (11) принимает следующий вид:
r 2hпл плW 2
.
Eпл 2
(12)
(13)
Приравнивая E и Eпл и учитывая, что m=mпл, получим что скорость метания при эквивалентных массах заряда и пластины определяется из выражения:
4 r 0
Wэкв U
.
(14)
3 hпл пл
Проведѐнные расчѐты по уравнению (14) показывают, что в случае эквивалентности масс пластины и заряда ВВ при торцевом метании справедливо следующее выражение:
Wэкв 1,41U .
(15)
В случае неэквивалентного метания, то есть когда масса пластины будет больше
или меньше массы заряда, скорость метания соответственно изменится. Данный метод
может быть предложен как наиболее удобный и обобщѐнный.
По методу эквивалентных масс рассчитана относительная метательная способность некоторых ВВ и сопоставлена с экспериментальными данными торцевого метания. В табл. 2 представлены расчѐтные и экспериментальные данные относительной
скорости метания медной пластины по методике М-60 относительно октогена. Расчѐтные данные скорости метания пластины были получены с использованием следующих
полуэмпирических моделей:
1) интерполяционное уравнение для определения скорости метания пластины, полученное Смирновым А.С. и др. при помощи регрессионного анализа массива экспериментальных данных:
ηМ-60 = 1,23ρз0,871Qк0,432Nг0,230 , %
(16)
где ηМ-60 – эффективность метательного действия по методике М-60 относительно октогена; ρз – плотность заряда, г/см3; Qк – теплота взрыва (калориметрическая), ккал/кг; Nг –
число молей газообразных продуктов взрыва, моль/кг. Коэффициент множественной
корреляции R=0,990, среднеквадратическая ошибка Δ=±0,5 (±2,65%).
2) гидродинамический метод обобщѐнных ударных адиабат для определения
скорости метания пластины WГД, основанный на объѐмной скорости звука;
3) метод эквивалентных масс Wэкв, разработанный по модели детонационной оптики.
7
Таблица 2
Расчетные и экспериментальные данные относительной скорости метания по методике М-60
№
п/п
ВВ
ηМ-60,
%
WГД,
%
Wэкв,
%
эксп.,
%
1
Октоген
100
100
100
100
Усреднение
%
100
2
Гексоген
94,3
96,2
97,4
97,0
95,9
+1,2
3
ТНТ
74,8
74,3
79,9
76,0
76,3
+0,4
4
Бис-(фтординитроэтил)нитрамин
94,7
95,1
101,2
93,7
97,0
+3,5
5
Бис-(хлординитроэтил)нитрамин
91,3
92,4
96,8
90,3
93,5
+3,5
6
78,3
75,9
77,7
77,2
77,3
+0,1
95,3
95,4
102,5
95,3
97,7
+2,5
8
1,4-динитропиперазин
2',2',2'-тринитроэтил-4,4,4тринитробутират
Бис-(тринитроэтил)формаль
89,5
89,1
95,0
87,7
91,2
+3,9
9
ТЭН
95,2
92,0
101,4
93,7
96,2
+2,6
10
ГНБ
105,8
110,3
118,7
109,7
111,6
+1,7
11
ТНБ
79,5
80,0
85,4
79,4
81,6
+2,7
12
ТАТБ
83,4
80,8
79,2
82,5
81,1
-1,7
13
Тетрил
85,8
88,2
90,2
85,7
88,0
+2,7
14
Гексанитростильбен
76,0
80,5
84,5
79,0
80,3
+1,6
15
3,3'-динитро-4,4'-азоксифуразан
102,6
103,7
105,4
101,7
103,9
+2,2
16
Бензотрифуроксан
3,4-бис(4-нитрофуразанил-3)фуроксан
Нитроглицерин
102,7
101,1
104,1
101,5
102,6
+1,1
92,9
105,4
102,1
100,8
100,1
+5,3
100,0
82,1
83,1
76,0
102,5
+1,2
7-амино-4,6-динитробензофуроксан
N,N-бис(тринитроэтил)этилендинитроамин
Бис-(2,2,2-тринитроэтил)нитрамин
88,8
85,9
89,6
88,2
88,1
-0,2
96,6
103,4
105,7
97,8
101,9
+4,2
95,5
99,5
103,4
95,4
99,5
+4,3
7
17
18
19
20
21
%
ош.
-
Среднеквадратическая ошибка по гидродинамическому методу: 2,75%; среднеквадратическая ошибка по методу эквивалентных масс: 3,33%.
Для уменьшения процента ошибки был использован метод усреднения расчѐтных
значений по значениям трѐх методик, что позволило уменьшить среднеквадратическую
ошибку до 1,69%. Таким образом, усреднение результатов расчѐта по трѐм методикам
увеличивает достоверность прогнозирования скорости метания пластины.
В третьей главе исследуется определение влияния материала пластины на параметры торцевого метания. В этом плане параметры соударения метаемого элемента с
преградой определяется ударной адиабатой соударяющихся тел. Одним из вопросов
разработки систем кинетического оружия является разработка ударника и определении
8
его свойств. Это также актуально для зарядов взрывного бурения, ударных возбудителей сейсмосигнала и кумулятивных перфораторов. В настоящее время наиболее изученным в качестве ударников являются медь (ударные ядра), сталь и вольфрам (бронебойные и подкалиберные снаряды). Характеристики же сплавов, в качестве ударников, изучены неполно, хотя многие, особенно на основе d-металлов, являются наиболее перспективными в этом плане.
Для определения скорости торцевого метания пластины по динамическим адиабатам метаемых элементов и дальнейшего сравнения результатов расчѐта с экспериментальными данными нами проведѐн и обоснован выбор метаемого материала. Наиболее
часто в качестве облицовки для СФЗ используется медь, а ближайшими к меди аналогом
по физико-механическим свойствам является семейство сплавов латуни. Благодаря высоким технологическим и механическим свойствам латуни являются самыми распространѐнными и дешѐвыми из медных сплавов. Проведѐнный нами анализ свойств двойных латуней показывает, что для исследования следует рассматривать сплавы с содержанием цинка менее 25%. Для определения ударных адиабат гидродинамическим методом в качестве объекта исследования выбрана латунь Л75 и еѐ компоненты медь и цинк.
В табл. 3 приведены расчѐтные значения латуни Л75 и еѐ компонентов.
Нами предложен алгоритм расчѐта скорости метания пластины гидродинамическим методом, который осуществляется в следующей последовательности:
1. Расчѐт параметров детонации по экспресс-методике (1), (2);
2. Построение динамической адиабаты с использованием политропы ЛандауСтанюковича P=Aρn для продуктов взрыва;
3. Построение динамической адиабаты метаемого элемента;
4. Нахождение параметров входящей ударной волны;
5. Нахождение скорости метаемого элемента W=2U.
Таблица 3
Металл
Zn
Cu
Cu+Zn
Плотность
, г/см3
7,14
8,90
8,63
Характеристики металлов
Коэффициент ударОбъемная скорость
звука С0, м/с
ной адиабаты 1,45
3300
1,49
3980
1,51
3725
Массовая
доля 0,235
0,765
-
При определении скорости торцевого метания пластины по динамическим адиабатам метаемых элементов нами показано влияние материала пластины на параметры
метания. Расчѐтная и экспериментальная ударные адиабаты для латуни Л75 и еѐ компонентов в D-U координатах приведены на рис. 2.
Как видно из рис. 3 имеет место близость хода экспериментальной и полученной
ударных адиабат латуни, с использованием расчѐтных значений объѐмной скорости звука и коэффициента λ, полученных методом «правила симметричного окружения» и
практически полного их совпадения в диапазоне массовых скоростей от 700 м/с до 1600
м/с. Именно в этом диапазоне лежат массовые скорости материала пластины, которые
для меди при еѐ нагружении зарядом A-IX-1 составляет 955 м/с, а для латуни 987 м/с, а
при нагружении данных облицовок зарядом ТГ50/50 – 730 м/с и 763 м/с соответственно.
Соответственно, кинетический модуль для единицы массы на ≈4,0% выше при метании
облицовки из латуни в случае А-IX-1, и на ≈8,0% в случае метания зарядом ТГ50/50.
9
Рис. 3. Ударные адиабаты. 1 – экспериментальная ударная адиабата меди; 2 – экспериментальная ударная адиабата цинка; 3 – экспериментальная адиабата латуни Л75; 4
– расчѐтная адиабата латуни Л75.
В третьей главе также описывается методика и проводится экспериментальное
определение способов регистрации взаимодействий между ВВ и метаемой пластиной
электромагнитным методом при торцевом метании. Усовершенствованный метод отличается от экспериментального измерения массовой скорости продуктов детонации,
предложенным Шведовым К.К., тем, что в данном случае на торец метаемой пластины
устанавливаются стальные иглы с базовым сечением 5 мм. Схема опыта приведена на
рис. 4.
В опытах использовались заряды ВВ диаметром 20 мм и длиной 100
мм, которые устанавливались на монтажной плате. Детонация инициировалась плосковолновым генератором.
Заряд компоновался из двух таблеток
прессованного порошка ВВ. В качеРис. 4. Схема экспериментальной сборки для
стве зарядов ВВ применялись: прессоопределения скорости торцевого метания:
ванные заряды А-IX-1 и состав
1 – электродетонатор; 2 – плосковолновой геТГ50/50. Шашки ВВ изготавливались
нератор; 3 – заряд ВВ; 4 – метаемая пластина;
методом двухстороннего прессования
5 – игольчатые датчики.
для
равномерного
распределения
плотности. В табл. 4 приведены характеристики используемых зарядов. В эксперименте
использовались медные и латунные пластины, которые крепились на торец заряда ВВ.
Характеристики метаемых пластин приведены в табл. 5.
Таблица 4
Характеристики зарядов
Характеристика
Заряд
А-IX-1
ТГ50/50
Скорость
Плотность,
детонации,
г/см3
м/с
1,67
8380
1,65
7560
10
Теплота взрыва,
кДж/кг
5518
4765
Давление в точке
Чепмена-Жуге,
ГПа
32,4
25,1
Монтажная плата размещалась в однородном поле постоянного магнита с напряженностью 337-342 эрстед, которое создавалось электромагнитами поперечного сечения
320х320 мм. Напряжѐнность контролировалась в каждой серии опытов.
Таблица 5
Характеристики метаемых пластин
Характеристика
Металл
Медь М1
Латунь Л75
Плотность,
г/см3
8,90
8,63
Масса,
г
9,79
9,48
Толщина,
мм
3,5
3,5
Диаметр,
мм
20
20
Из таблицы видно, что толщина метаемых пластин составляла 3,5 мм, что обеспечивает затухание химпика в материале пластины.
При взрыве баратоловой линзы диаметром 20 мм в основном заряде, состоящем
из двух таблеток ВВ толщиной 40 мм, инициировалась детонационная волна с плоским
фронтом. К метаемому элементу вплотную устанавливались два игольчатых датчика
диаметром 1 мм и длиной 50 мм. Расстояние между датчиками составляло 5 мм. При
движении метаемая пластина с игольчатыми датчиками пересекает магнитные силовые
линии, при этом генерируется ЭДС, которая фиксировалась на цифровом осциллографе.
База, на которой проводилось измерение, равна длине игольчатых датчиков и составляет
50 мм.
Во всех проведѐнных экспериментах чѐтко фиксировались три ударно-волновых
взаимодействия ПД с пластиной. На рис. 5 представлены примеры осциллограмм экспериментов.
Рис. 5. Осциллограммы скорости торцевого метания продуктами взрыва заряда A-IX-1
при метании: слева) медной пластины М1; справа) латунной пластины Л75.
На рис. 6 приведена схема типичной осциллограммы многократного разгона металлической
пластины в зависимости скорости метания W от
времени t. Скорость метаемых пластин рассчитывалась из соотношения
W
lH
10 8 , см/с,
(17)
где ε – величина наведѐнной ЭДС, В; H – напряжѐнность магнитного поля, Э; l – длина датчика,
см.
11
Рис. 6. Схема типичной осциллограммы трѐхкратного взаимодействия ПД
с метаемой пластиной
Экспериментальные значения скорости метания пластин продуктами детонации
приведены в табл. 6.
Таблица 6
Экспериментальные значения скорости метания пластин
продуктами детонации
№
Материал пла- 1-й импульс
ВВ
эксп.
стины
W1, м/с
1.
Медь (М1)
1930
2.
Медь (М1)
1901
3.
Медь (М1)
1912
4.
Медь (М1)
1895
A-IX-1
5.
Латунь (Л75)
2041
6.
Латунь (Л75)
2064
7.
Латунь (Л75)
2070
8.
Латунь (Л75)
2053
9.
Медь (М1)
1444
10.
Медь (М1)
1450
11.
Медь (М1)
1439
12.
Медь
(М1)
1450
ТГ
50/50 Латунь (Л75)
13.
1544
14.
Латунь (Л75)
1532
15.
Латунь (Л75)
1526
16.
Латунь (Л75)
1538
2-ой импульс
W2, м/с
2398
2357
2368
2345
2532
2550
2561
2544
1795
1795
1784
1801
1918
1912
1889
1901
3-й импульс
W3, м/с
2567
2544
2561
2526
2737
2737
2749
2743
1912
1930
1918
1942
2058
2035
2018
2041
W2/W1 W3/W1
1,242
1,239
1,238
1,237
1,240
1,235
1,237
1,239
1,243
1,237
1,239
1,242
1,242
1,248
1,238
1,236
1,330
1,338
1,339
1,333
1,341
1,326
1,328
1,336
1,324
1,331
1,333
1,339
1,333
1,328
1,322
1,327
Как видно из табл. 6 во всех случаях имеет место доразгон метаемых пластин при
их трѐхкратном нагружении. В случае увеличения диаметра заряда ВВ число взаимодействий, вероятно, будет больше, однако вклад в энергетику пластины будет минимальным.
В табл. 7 приведено сравнение полученных нами экспериментальных данных
скорости метания зарядом ТГ50/50 меди М1 с экспериментальными данными Lefrancois
A. (Lefrancois A., Baudin G., Bouinot P.), полученными при помощи лазерного интерферометра VISAR.
Таблица 7
Сравнение экспериментальных значений скоростей метания медных пластин зарядом ТГ50/50
ВВ
Экспериментальные данные
ТГ50/50
Экспериментальные данные
ТГ50/50 Lefrancois A.
м/с
Ошибка
%
1-ый импульс
W1, м/с
2-ой импульс
W2, м/с
3-й импульс
W3, м/с
1446
1794
1925
1500
1860
2010
-54
3,6
-66
3,6
-85
4,2
Как видно из табл. 7 значения импульсов зафиксированных при помощи осциллографа имеют удовлетворительную сходимость с импульсами, зафиксированными лазерным интерферометром.
При многократном ударно-волновом нагружении медные и латунные пластины
доразгоняются во втором импульсе примерно на 24%, в третьем – на 8%. При этом
средний доразгон пластин составляет ≈33% от скорости метания в первом импульсе.
12
Эти данные, наряду с измерением скорости ударной волны, позволяют определить всю картину нагружения среды. Данный метод представляется целесообразным
использовать для оценочных практических расчѐтов скорости торцевого метания пластины, при этом можно использовать значение равные приблизительно 2,66U от массовой скорости выходящей на свободную поверхность ударной волны. Метод не накладывает ограничений на размеры заряда ВВ, используемого для создания в материале ударной волны. Погрешность измерений определяется точностью измерения геометрических
размеров метаемой пластины и ВВ, а также взаимного расположения источников постоянного магнитного поля и игольчатых датчиков.
Оптимальными методами теоретического прогнозирования характеристик метания следует считать, как было отмечено выше, регрессионный анализ массива экспериментальных данных, метод детонационной оптики и гидродинамический метод оценки,
последний является наиболее важным, но с физической точки зрения наиболее адекватным, хотя и требует знания целого ряда параметров ВВ и метаемых элементов.
Интенсивные поисковые работы последних десятилетий привели к синтезу и исследованию свойств большого числа новых классов энергетических соединений, обладающих чрезвычайно разнообразным и весьма различающимся строением в очень широком диапазоне энергетических, взрывчатых и физико-химических свойств. Одним из
таких классов соединений являются производные 1,3,5-триазинов. Поэтому представляется целесообразным применить разработанные нами методы определения скорости метания элемента (15) и разлѐта осесимметричной оболочки (6) на классе производных
1,3,5-триазинов.
В четвѐртой главе проведено определение влияния ВВ ряда 1,3,5-триазинов на
параметры торцевого метания. Осуществлена расчѐтная оценка свойств производных
1,3,5-триазинов. Оценка перспективности ВВ заключается в проведении комплекса исследований, по результатам которых возможно прогнозирование поведения взрывного
изделия при различных видах воздействий.
Оценка взрывчатых характеристик производных 1,3,5-триазинов показала, что эти
вещества имеют перспективы для использования во взрывчатых составах. Структурные
формулы производных 1,3,5-триазинов приведены на рис. 6.
OCH2C(NO2)3
OCH2C(NO2)3
N
(NO2)3C
N
N
N
OCH2C(NO2)3
I
2,4-ди(тринитроэтокси)-6тринитрометил-1,3,5-триазина
(NO2)3CH2CO
H3CO
N
OCH2C(NO2)3
VI
2-метокси-4,6-бис(тринитроэтокси)-1,3,5-триазин
N
N
OCH2C(NO2)3
N3
H2N
OCH2C(NO2)3
N
OCH2C(NO2)3
V
2-амин-4,6-бис(тринитроэтокси)-1,3,5-триазин
OCH2C(NO2)3
III
2-азид-4,6-бис(тринитроэтокси)-1,3,5-триазин
N
N
N
N
OCH2C(NO2)3
N
N
N
II
2,4,6-трис-(тринитроэтокси)1,3,5-триазин
OCH2C(NO2)3
N
OCH2C(NO2)3
F
N
N
OCH2C(NO2)3
VI
2-фтор-4,6-бис(тринитроэтокси)-1,3,5-триазин
Рис. 6. Структурные формулы производных 1,3,5-триазинов
Экспериментальные и расчѐтные характеристики производных 1,3,5-триазинов
приведены в табл. 8.
13
Таблица 8
Экспериментальные и расчѐтные характеристики
производных 1,3,5-триазинов
Вещество
Характеристика
k
0, г/cм
3
Тпл, K
С0, м/с
Qm, кДж/кг
D1,8, м/с
U, м/с
Qдет.вз, кДж/кг
Р, ГПа
N
I
II
III
IV
V
VI
1,111
1,86
442 (с разл.)
1685
4376
7840
2065
4264
30,7
2,53
1,000
1,84
414 (с разл.)
1703
5625
8410
2278
5189
32,9
2,45
0,875
1,80
389
1596
5637
8340
2201
4844
29,8
2,42
0,769
1,74
332
1664
5537
8200
2099
4406
26,8
2,50
0,824
1,82
389
1694
5247
8070
2114
4469
28,6
2,51
0,903
1,80
385
1573
4741
7890
2030
4121
25,8
2,48
Как видно из табл. 8 вещества I-VI все рассмотренные ВВ имеют плотность в
пределах 1,74-1,86 г/см3 и различный кислородный коэффициент αk от 0,77 до 1,11 и
имеют различные при этом температуры плавления Тпл. Вещество IV имеет низкую
Тпл=59oC и в принципе перспективен как плавкий взрывчатый компонент, а образец I
перспективен в плане взрывчатого окислителя.
Скорость звука в I-VI близки между собой и лежат в пределах от 1573 до 1703
м/с. По этому показателю они не относятся к классу динамиков. Максимальная теплота
взрывчатого превращения Qm лежит в пределах 4376-5637 кДж/кг, а детонационная теплота взрыва Qдет.вз, определяющаяся квадратом массовой скорости, имеет пределы 41215189 кДж/кг.
Производные 1,3,5-триазинов представляют из себя интересные ВВ, которые
имеют достаточно высокую массовую скорость и низкую скорость звука, поэтому представляется интересным определить скорости метания данных веществ по разработанным в диссертации методикам и статистическому методу Смирнова А.С.
В табл. 9 приведены результаты расчѐтов:
– относительной скорости метания пластины по методу эквивалентных масс Wэкв
с использованием уравнения (15);
– относительной скорости торцевого метания ηМ-60 по интерполяционному уравнению Смирнова А.С. (16);
– скорости расширения цилиндрической оболочки Wcyl по формуле (6);
– скорости расширения цилиндрической оболочки по интерполяционному уравнению Смирнова А.С. WТ-20 по уравнению (7).
Таблица 9
Расчѐтные значения скоростей метания меди по методикам М-60 и Т-20 производными 1,3,5-триазина
Производное
1,3,5-триазина
Wэкв, % по отношению к октогену
ηМ-60, % по отношению к октогену
Wcyl, м/с
WТ-20, м/с
I
II
III
IV
V
VI
93,6
91,5
1410
1496
102,5
98,5
1497
1606
97,9
93,1
1427
1534
91,8
89,4
1389
1478
94,5
94,9
1432
1564
90,3
88,5
1355
1455
Как видно из табл. 9 наиболее перспективным является вещество II и, вероятно,
III, а в плане окислителя вещество I. Вероятно, комбинация взрывчатой смеси из III и V
14
превзойдѐт по метательной способности II, так как скорость звука в системе будет
меньше, чем во II. Таким образом, данные вещества могут найти применение, например,
в качестве облицовок для снарядоформирующих зарядов.
В табл. 10 приведены параметры относительной метательной способности данных
ВВ и скорости расширения цилиндрической оболочки из латуни Л75.
Таблица 10
Расчѐтные значения скоростей метания латуни Л75
по методикам М-60 и Т-20 производными 1,3,5-триазина
Производное 1,3,5-триазина
Wэкв, % по отношению к октогену
Wcyl, м/с
I
93,6
1524
II
102,5
1618
III
97,9
1543
IV
91,8
1500
V
94,5
1548
VI
90,3
1465
Как видно из табл. 10 при нагружении латуни Л75 зарядами производных 1,3,5триазинов наблюдается увеличение расчѐтной скорости расширения цилиндрической
облицовки на 8%.
Таким образом, разработанный экспресс-метод оценки эффективности действия
вполне применим для металлов, которые могут быть использованы в качестве метаемых
элементов типа ударное ядро, и даѐт возможность определить наиболее оптимальный
материал, хотя и требует экспериментальной проверки.
ВЫВОДЫ:
1. С целью разработки практических и теоретических способов определения метательной способности проведѐн критический анализ существующих методов определения параметров метания пластин продуктами детонации. Показано, что существующие
практические методы определения скорости торцевого метания обладают существенным
недостатком – методы не регистрируют динамику разгона пластины, а из имеющихся
теоретических методов определения относительной метательной способности наиболее
применимыми являются регрессионный анализ, метод детонационной оптики и гидродинамический метод. Последние два метода не разработаны досконально.
2. Разработан принцип расчѐта скорости метания гидродинамическим методом,
включая оценку параметров детонации ВВ на основе прогнозирования скорости звука в
зарядах ВВ и в материале метаемой пластины. Показано, что наиболее приемлемой является методика расчѐта параметров детонации с использованием скорости звука ВВ, так
как позволяет рассчитать массовую скорость с максимально возможной точностью.
3. Зарегистрировано трѐхкратное ударно-волновое взаимодействие продуктов
детонации с метаемым элементом и определены процентные вклады каждого взаимодействия продуктов детонации с пластиной, которые составляют около 24% во втором импульсе и ≈10% в третьем. На основании разработанной методики определения доразгона
пластины при торцевом метании экспериментально определены скорости медной и латунной пластин при многократном нагружении их продуктами детонации. Показано, что
кроме параметров ВВ на эффективность действия влияет материал ударника. В частности, использование ударника из латуни Л75 позволяет увеличить импульс на 4-8% по
сравнению с медью М1.
4. Усовершенствован способ прогнозирования физических и динамических характеристик металлических элементов на основе периодического закона Д.И. Менделеева. Разработана экспресс-методика определения смесевых динамических адиабат сплавов, которая даѐт возможность определить наиболее оптимальный состав компонентов
для реализации максимальной эффективности метательного действия;
15
5. Предложены методики расчѐта относительной скорости торцевого метания
гидродинамическим методом и методом детонационной оптики, а так же приближѐнный
метод определения скорости расширения цилиндрической облицовки Т-20. Разработаны
теоретические методы прогнозирования скоростей метания по методам Т-20 и торцевого
метания М-60;
6. Проведено сравнение предложенных методик расчѐта скоростей метания с
известными интерполяционными уравнениями Смирнова А.С.. Показано, что предложенные методы расчѐта метательного действия имеют равную точность с методом
Смирнова А.С., а усреднение результатов теоретического прогнозирования скорости метания, полученных по разным методикам, снижает среднеквадратическую ошибку в
среднем на 1,5%.
7. На основе разработанных методик произведѐн расчѐт параметров детонации и
скоростей метания по методикам Т-20 и М-60 для ряда перспективных взрывчатых веществ производных 1,3,5-триазинов. Показано, что вещества обладают хорошей метательной способностью, а 2,4,6-трис-(тринитроэтокси)-1,3,5-триазин имеет скорость метания порядка октогена.
СПИСОК РАБОТ, ОПУБЛИКОВАННЫХ ПО ТЕМЕ ДИССЕРТАЦИИ
Публикации в изданиях, рекомендованных ВАК РФ:
1. Реут И.И., Гидаспов А.А., Кривченко А.Л., Кожевников Е.А. Расчет детонационных характеристик и параметров метательной способности производных 1,3,5триазинов // Вест. Сам. гос. техн. ун-та. Сер.: Технические науки. 2011. №4 (32), С. 110114.
2. Реут И.И., Кривченко А.Л. Расчѐт метательной способности взрывчатых веществ при цилиндрическом и торцевом метании металла // Вест. Сам. гос. техн. ун-та.
Сер.: Физ.-мат. науки. 2011. №4 (25), С. 173-177.
3. Реут И.И., Кривченко А.Л. О методе оценки латуни как материала для облицовки боеприпасов типа «ударное ядро» // Вест. Сам. гос. техн. ун-та. Сер.: Физ.-мат.
науки. 2011. №4 (25), С. 188-200.
Публикации в других изданиях и материалах научно-технических конференций:
4. Гидаспов А.А., Ермаков С.Ю. Реут И.И. Замещение спиртами и фенолами
тринитрометильной группы в тринитрометил-1,3,5-триазинах без применения в реакции
оснований // Тр. VII Всероссийской н.-т. конф. «Наука. Промышленность. Оборона». Новосибирск: НГТУ. 2006. С. 155-156.
5. Гидаспов А.А., Каемов С.А., Ермаков С.Ю., Реут И.И. Синтез 2-изопропокси4-тринитроэтокси-6-тринитрометил-1,3,5-триазина // Тез. докл. XXXII Сам. обл. студ. н.т. конф. памяти проф. Х.С. Хазанова. «Общественные, естественные и техн. науки». Самара: 2006. ч.1. С. 126.
6. Гидаспов А.А., Ермаков С.Ю., Реут И.И. Катализ солями тринитрометана замещения тринитрометильных групп спиртами и фенолами в тринитрометил-1,3,5триазинах // Сб. науч. тр. «Успехи в химии и химической технологии». М.: РХТУ им.
Д.И. Менделеева. 2007. т. XXI, №5. С. 9-11.
7. Кривченко А.Л., Кривченко Д.А., Реут И.И., Чуркин О.Ю. Расчѐт ударных
адиабат d-металлов и их сплавов с использованием периодического закона Д.И. Менделеева // Сб. тез. междунар. конф. «Ударные волны в конденсированных средах». СПб.
2008. С. 253-255.
16
8. Кривченко А.Л., Кривченко Д.А., Реут И.И., Чуркин О.Ю. О возможности
расчѐта динамических характеристик d-металлов и их сплавов // Сб. тез. докл. междунар.
конф. «XI Харитоновские тематические научные чтения». Саров. 2009. С. 229-230.
9. Кривченко А.Л. Гидаспов А.А. Реут И.И., Заломленков В.А.О способе прогнозирования характеристик производных 1,3,5-триазинов // Сб.тез.межд.конф.«Shock
waves in condensed matter».Novgorod.2010.С.55-57
10. Реут И.И., Бертяев Б.И. Об одном феноменологическом подходе к расчѐту
кинетической энергии ядра атома в кристаллах с ОЦК и ГЦК решѐтками // Сб. тез. междунар. конф. «Shock waves in condensed matter». Novgorod. 2010. С. 308-310.
11. Кривченко А.Л., Кривченко Д.А., Реут И.И., Чуркин О.Ю. Определение параметров детонации в гомогенных и гетерогенных системах // Сб. тез. докл. междунар.
конф. «XIII Харитоновские тематические научные чтения». Саров. 2011. С. 70-72.
12. Кривченко А.Л., Гидаспов А.А., Реут И.И., Заломленков В.А. Метод расчѐта
параметров торцевого метания пластин из различных материалов // Сб. тез. докл. междунар. конф. «XIII Харитоновские тематические научные чтения». Саров. 2011. С. 255257.
17
Автореферат отпечатан с разрешения диссертационного совета Д 212.217.01
ФГБОУ ВПО «Самарский государственный технический университет»
(протокол №5 от «20» марта 2012 г.)
Заказ №314. Тираж 100 экз.
Отпечатано на ризографе.
ФГБОУ ВПО «Самарский государственный технический университет»
Отдел типографии и оперативной печати
443100 г. Самара ул. Молодогвардейская, 244
Документ
Категория
Технические науки
Просмотров
65
Размер файла
699 Кб
Теги
кандидатская
1/--страниц
Пожаловаться на содержимое документа