CC BY
65
МЕХАНИКА ДЕФОРМИРУЕМОГО ТВЕРДОГО ТЕЛА
УДК 623.4.012; 004.942
АНАЛИЗ ВЛИЯНИЯ ТЕХНОЛОГИЧЕСКИХ ПОГРЕШНОСТЕЙ МЕНИСКОВЫХ ОБЛИЦОВОК НА ДИНАМИКУ ВЗРЫВНОГО ФОРМИРОВАНИЯ ВЫСОКОСКОРОСТНЫХ СТЕРЖНЕВЫХ ЭЛЕМЕНТОВ
Н.А. Асмоловский, В.Д. Баскаков, О.В. Зарубина
МГТУ им. Н.Э. Баумана, Москва, Российская Федерация e-mail: Nikolay.Asmolovskiy@gmail.com; baskakov_vd@mail.ru; ov_zarubina@mail.ru
На основе численного трехмерного моделирования взрывного нагружения менисковых облицовок в системе координат Лагранжа исследовано влияние раз-нотолщинности облицовок в окружном направлении, представляемой в виде суммы элементарных тригонометрических гармоник, на кинематические характеристики и геометрические параметры формируемых высокоскоростных элементов. Реализован подход к моделированию малых по амплитуде гармонических составляющих разнотолщинности облицовки, основанный на принудительном незначительном искажении осесимметричной расчетной сетки. Для оценки формы высокоскоростных элементов, их поперечной и угловой скоростей движения разработана методика построения проекций контура высокоскоростного элемента на координатные плоскости и оценки положения его оси симметрии. С помощью разработанной методики установлена наиболее значимая кинематическая характеристика — угловая скорость, приобретаемая высокоскоростными элементами под влиянием разнотолщинности облицовок. Проведенные расчеты с учетом одновременного присутствия нескольких гармоник в структуре разнотолщинности облицовок позволили квалифицировать суперпозицию гармоник как одну из причин возможного нарушения симметрии высокоскоростных элементов в радиальном и осевом направлениях, а также их разрушения.
Ключевые слова: высокоскоростной элемент, периодические складки, искажение сеток, разнотолщинность облицовки.
RESEARCH INTO THE EFFECT OF TECHNOLOGICAL IMPERFECTIONS OF MENISCUS LINERS ON EXPLOSIVE FORMATION DYNAMICS OF HIGH-SPEED ROD ELEMENTS
N.A. Asmolovsky, V.D. Baskakov, O.V. Zarubina
Bauman Moscow State Technical University, Moscow, Russian Federation e-mail: Nikolay.Asmolovskiy@gmail.com; baskakov_vd@mail.ru; ov_zarubina@mail.ru
The paper presents the research into the effect of liners thickness variation in circumferential direction on kinematic characteristics and geometrical parameters of the generated high-speed elements. The liners thickness variation is presented as a sum of elementary trigonometric harmonics. The analysis is based on a threedimensional numerical simulation of explosive loading of the liners coatings in Lagrangian coordinates. The authors implement a method for simulation of small-amplitude harmonic components of the liners thickness variation based on the forced minor distortion of the axisymmetric mesh. The article discloses a technique
72 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2015. №5
projecting the contour of a high-speed element onto coordinate planes and estimating position of its axis of symmetry. The technique is used to estimate the shape of the high-speed elements as well as their transverse and angular velocities. With the help of this technique, the authors calculate the most important kinematic characteristic, that is an angular velocity, which the high-speed elements attain due to variation in the liners thickness. The calculations made with due account for simultaneous presence ofseveral harmonics in the structure of the liners thickness variation showed that the harmonics superposition can be considered as one of the reasons for the high-speed elements asymmetry in both radial and axial directions as well as for their destruction.
Keywords: high-speed element, periodical wrinkles, grid distortion, coating thickness variation.
Формируемые взрывом высокоскоростные элементы (ВСЭ) могут применяться для решения различных задач разрушения, например, для разделения на фрагменты крупных объектов, дробления горных пород и т.п. Применяются ВСЭ и в космических исследованиях. Например, стартовавшая 03.12.2014г. к астероиду 1999JU3 космическая научная миссия Хаябуса-2 с целью доставить на Землю образцы грунта планирует осуществить его дробление c помощью ВСЭ, который при взаимодействии с поверхностью сформирует большое число осколков различных фракций [1].
Высокоскоростные элементы образуются во время срабатывания взрывного устройства (ВУ) (1) при обжатии осесимметричной облицовки менисковой формы продуктами детонации заряда взрывчатого вещества (ВВ) в компактное или удлиненное тело с внутренней полостью и складчатой кормовой (хвостовой) частью (2), имеющее осевую скорость порядка 2000 м/с (рис. 1). Процесс формирования ВСЭ характеризуется высокой чувствительностью к технологическим погрешностям узлов и деталей (погрешностям формы и взаимного расположения поверхностей, неоднородностям свойств материалов и др.), составляющих ВУ. В результате возможна потеря устойчивости облицовки при нагружении продуктами детонации (аналогично потере устойчивости в классической задаче нагружения тонкостенного цилиндра внешним давлением [2]). Кроме того, материал облицовки находится в состоянии интенсивного пластического течения, поэтому неидеальные условия приложения нагрузки оказывают значительное влияние на форму и кинематику ВСЭ.
Несмотря на значительный объем исследований, выполненных для изучения воздействия технологических погрешностей ВУ на ВСЭ
Рис. 1. Устройство (7) для формирования ВСЭ (2); ТИ — точка инициирования
ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение” 2015. №5 73
(см., например, [3-5]), не в полной мере изучены закономерности формирования поперечной и угловой скоростей ВСЭ, особенности влияния на его форму различных гармонических составляющих погрешностей, а также их суперпозиции и ряд других важных вопросов.
Цель настоящей работы — детальное изучение влияния разностен-ности менисковой облицовки (ее разнотолщинности в окружном направлении) на формирование ВСЭ с помощью численных расчетов и с использованием пре- и постпроцессоров, адаптированных под особенности структуры технологических погрешностей ВУ и формы ВСЭ.
Методы решения задачи и принятые допущения. 1. Моделирование асимметрий технологической природы и обоснование метода численного моделирования. Технологические погрешности деталей ВУ принято представлять в виде гармоник тригонометрического ряда [3, 4]:
N
S ^ So + an cos (up + фп) , (1)
n=1
где S — изменение параметра детали ВУ (например, толщины облицовки) в плоскости его поперечного сечения; 80 — номинальное значение параметра конструкции ВУ в плоскости его поперечного сечения; n = 1, 2,... ,N — номер гармоники, N — максимальный номер значимой гармоники; ап — амплитуда n-й гармоники; p — угловая координата в плоскости поперечного сечения ВУ, изменяющаяся в диапазоне 0 < p < 2п; фп — начальная фаза n-й гармоники.
Ввиду неизбежного присутствия технологических погрешностей, их статистической природы, а также дороговизны экспериментов, целесообразно проводить исследования с помощью численного моделирования.
Формирование ВСЭ характеризуется следующими особенностями:
• высокими значениями деформаций и скоростей деформаций;
• критическим состоянием материала (пластическим течением и разрушением);
• высокой чувствительностью процесса к неравномерностям, в том числе и технологического происхождения;
• высокоскоростным взаимодействием деталей ВУ.
При моделировании задач такого рода общепринятым является подход Эйлера [5-7]. Однако его использование для решения трехмерных задач с учетом неравномерностей малой амплитуды не представляется возможным из-за чрезвычайно высоких требований к числу элементов дискретизации расчетной области.
Альтернативой является подход Лагранжа, при котором узлы сетки перемещаются совместно с материалом, обеспечивая более высокую точность расчета границы тела по сравнению с подходом Эйлера, что
74 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2015. №5
необходимо для решения задач о влиянии неравномерностей малой амплитуды. Однако ввиду высоких деформаций и скоростей деформаций необходимо применение дополнительных мер по контролю шага по времени, например, путем удаления сильно искаженных элементов с критически малым характерным размером. Кроме того, на границах раздела заряда ВВ и облицовки, заряда ВВ и корпуса были заданы граничные условия прилипания. Прочие особенности применения подхода Лагранжа к моделированию процесса взрывного формирования описаны в работах Джонсона [8, 9]. Так, например, было продемонстрировано, что качество и вид дискретизации оказывают значительное влияние на вид сформированного элемента [10]. Упрощенный подход Лагранжа также встречается в работах [11, 12].
Оценка чувствительности традиционно является одной из самых сложных задач численного моделирования. В технических системах, склонных к потере устойчивости, принято рассматривать состояние системы в зависимости от энергии пертурбаций (возмущений), причем особый интерес представляют возмущения малой амплитуды. Часто выделяют два типа возмущений: неравномерность геометрии и неравномерность нагрузки. Чувствительность оценивают путем внедрения в идеальную модель неравномерности различной природы и отслеживания параметров системы, варьируя значениями вносимой неравномерности. Такой подход использовался, например, для задачи оценки чувствительности оболочек [13], однако к исследованию поведения ВУ широко не применялся.
2. Особенности моделирования процесса функционирования ВУ и обработка результатов расчетов. Оценка влияния погрешностей облицовки осуществляется в несколько этапов.
1. Внедрение неравномерностей путем искажения расчетной сетки.
2. Решение задачи формирования стандартными решателями.
3. Оценка параметров кинематики и формы ВСЭ.
В настоящей работе рассматривается влияние двух типов разно-стенности облицовки, обусловленных несоосностью внешней и внутренней поверхностей облицовки (гармоника n = 1 в структуре зависимости (1)) и распределением толщины в периферийной части облицовки, вызванным одной из старших гармоник (n > 2) в соотношении (1). Погрешность n = 1 приводит к появлению радиальной и угловой скоростей ВСЭ, а погрешности n > 2 провоцируют формирование складок в его хвостовой части, улучшая аэродинамическую устойчивость [3, 4]. Технологические погрешности ВУ могут также привести к разрушению элемента во время полета [3].
В связи с повышенными требованиями к расчетной сетке в качестве исходной дискретизации ВУ используется сетка, состоящая из
ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение” 2015. № 5 75
элементов с 8 узлами. На первом этапе расчетные области облицовки и заряда имеют цилиндрическую форму [3]. Далее узлы сетки перемещаются в соответствии с заданными радиусами кривизны облицовки и гармониками погрешностей. Математическое описание алгоритма перемещения представлено ниже. Пояснительная расчетная схема изображена на рис. 2. Аналогичным образом происходит также преобразование сетки заряда ВВ.
3. Алгоритм генерирования расчетной сетки с погрешностями облицовки.
Входные данные:
• x,y,z — начальные координаты узла сетки облицовки, причем сетка облицовки имеет регулярную цилиндрическую конфигурацию,
• S0 — толщина облицовки,
• а1, ап — безразмерные амплитуды погрешностей гармоник n _ 1
и n > 1 (an an/SoX
• R1 — радиус кривизны облицовки,
• Rl — радиус внешней поверхности облицовки,
• OA — максимальная z-координата узла сетки облицовки,
• OB — минимальная z-координата узла сетки облицовки,
• n — число неравномерностей в окружном направлении. Выходные данные:
• X, Y, Z — обновленные координаты узла.
Алгоритм.
• Определение координат узла в цилиндрической системе координат — R _ у/x2 + у2, ^ = arctg (y/x).
• Определение перемещения узла, требуемого для задания радиуса кривизны облицовки, — h _ у/Ri2 — R2 — у/Rl2 — RL2.
^ Случай неравномерности с внешней стороны облицовки:
f _ So — (OA — z)
f So .
^ Случай неравномерности с внутренней стороны облицовки:
f _ OA — z f _ .
^ Случай неравномерности с обеих сторон облицовки: f _ 1.
Рис. 2. Пояснительная схема к алгоритму генерирования сетки
76 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2015. №5
Определение гармоники неравномерности толщины: 5 =
= (R/Rl )an5оcos (np + фп).
„ v _ z - OB
Введение погрешности несоосности: X = x + a1 ----.
5о
Введение погрешности разнотолщинности и кривизны облицовки: Z = z — (h + 50 f).
Изменения Y-координаты узлов не происходит: Y = у.
Результирующие сетки ВУ и облицовок показаны на рис. 3 (в целях визуализации погрешности облицовки на рис. 3, б и в увеличены).
Разностенность n > 2, образуемая выступами на наружной поверхности облицовки при гладкой внутренней поверхности, оказывает наибольшее влияние на форму ВСЭ [3]. Поэтому в настоящей работе, наряду с разностенностью, обусловленной несоосностью наружной и внутренней поверхностей облицовки (n = 1), рассматривался и этот вариант задания разностенности, в котором для определенности принималось n = 6. В производственных условиях гармоника n = 1 может формироваться на стадии штамповки мениска в инструментальном штампе вследствие несоосности рабочих поверхностей штампа. Причиной образования гармоники n = 6 может быть сила закрепления облицовки в 6-лепестковом цанговом патроне при механической обработке.
Расчет процесса взрывного формирования ВСЭ проводился с помощью решателя LS-Dyna. Подробное описание постановки задачи, включая параметры ВУ, начальные и граничные условия, а также принятые допущения, приведены в [3].
Для оценки кинематических параметров сформированного взрывом ВСЭ разработан специализированный алгоритм, учитывающий все особенности процесса формирования. Стоит отметить, что вследствие разрушения элемента и отсутствия симметрии в общем случае
Рис. 3. Результирующая сетка ВУ:
а — сетка заряда (половина); б — сетка облицовки при несоосности поверхностей; в — сетка облицовки с неравномерностью толщины
ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение” 2015. №5 77
применение стандартных методов оценки его кинематики может привести к неверным результатам.
Разработанный алгоритм состоит из следующих этапов. 1. Преобразование результатов расчетов LS-Dyna в матричный вид.
2. Выборка данных, описывающих ВСЭ.
3. Для каждого шага по времени: расчет координаты центра масс, скорости центра масс; определение геометрии внешнего контура; расчет угла поворота и угловой скорости.
Все промежуточные результаты вычислений сохраняются в целях увеличения скорости расчета при последующих обращениях. Рассмотрим подробнее некоторые из перечисленных этапов.
Выборка данных, описывающих ВСЭ. В процессе формирования облицовка претерпевает значительные деформации, в результате которых может происходить ее локальное разрушение. Например, на начальных этапах значительная доля периферийной части облицовки отделяется от более быстрого фрагмента, формирующего ВСЭ. Наибольший интерес представляет именно сам ВСЭ. Как показали расчеты, часть облицовки, из которой формируется головная часть ВСЭ, меньше всего подвержена разрушению. Базируясь на этом наблюдении, выделение элементов, представляющих ВСЭ, может быть проведено в несколько этапов.
1. Построение графа, в котором вершинами являются конечные элементы, а дугами соединены соседние элементы (т.е. элементы, имеющие общие узлы). В случае разрушения облицовки, полученный граф имеет области, которые не связаны друг с другом.
2. Все элементы, обнаруженные во время обхода графа, начиная с элемента, представляющего мыс головной части, принадлежат дискретизации ВСЭ. В настоящей работе применяется алгоритм поиска в ширину [10].
Расчет угла поворота ВСЭ непосредственно при формировании. За угол поворота принимался угол между осью ВУ и условной осью ВСЭ. Сформированный ВСЭ имеет лишь слабую симметрию, поэтому строгое определение положения его оси зачастую невозможно. За ось был принят вектор, соединяющий центр масс ВСЭ с его головной частью.
Следует отметить, что в процессе формирования происходит “выворачивание” облицовки. В таком случае существует состояние, при котором расстояние между центром масс и носовой частью ВСЭ мало, что приводит к сингулярности вычислений. Вследствие чего вычисление угла поворота ВСЭ в непосредственной близости его центра масс и носовой части по предлагаемому методу невозможно, поэтому угол определялся с момента времени t = 0,1 мс, т.е. после окончания формирования ВСЭ.
78 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2015. №5
Расчет угловой скорости. Вычисление угловой скорости уместно лишь после окончания процесса формирования ВСЭ. В целях упрощения расчета за результирующую угловую скорость целесообразно принять скорость изменения угла поворота ВСЭ.
Определение геометрии внешнего контура. Задача сводится к определению границы проекции ВСЭ на координатные плоскости. Многогранник, образованный узлами сетки ВСЭ, в общем случае не является выпуклым. Поэтому нахождение границы возможно при использовании сеточных алгоритмов, в соответствии с которыми: строятся проекции узлов сетки ВСЭ на координатную плоскость; совмещаются точки, находящиеся в малой окрестности друг от друга; строятся триангуляции Делоне [14] для множества спроецированных точек; для каждого ребра триангуляции определяется число элементов, содержащих данное ребро; ребра триангуляции, принадлежащие единственному элементу, — это фрагменты искомой границы.
Совмещение точек выполняется с помощью k-d деревьев [15]. В таком случае сложность алгоритма определения внешнего контура равна O(N log N), где N — число точек.
Результаты расчетов и их обсуждение. Расчеты проводились в несколько этапов. Вначале исследовалось влияние разностенности в варианте несоосности наружной и внутренней поверхностей облицовки (п = 1) на поведение ВСЭ. Далее с целью изучить влияние суперпозиции гармоник на ВСЭ, дополнительно к данной гармонике вводилась симметричная гармоника п = 6, не влияющая на радиальную и угловую скорости движения ВСЭ и вызывающая образование складок в его кормовой части [3].
Несоосность наружной и внутренней поверхностей облицовки. В качестве меры несоосности принимается безразмерный параметр а 1 = а,\/50 (см. зависимость (1)). Известно, что несоосность облицовки является одной из наиболее важных погрешностей, влияющих на кинематику ВСЭ [16].
В целях изучения влияния а1 на кинематические параметры ВСЭ проведены расчеты в диапазоне погрешностей, характерном для современных технологий изготовления облицовок. В первую очередь рассматривалось влияние несоосности на осевую V и радиальную Vp скорости. Расчетами установлено, что а1 не оказывает значительного влияния на Vo (коэффициент корреляции r = 0,471), в то время как между а1 и Vp имеется корреляция: r = 0,987. Однако радиальная скорость при значительной асимметрии а1 = 0,02 составляет лишь
1,5 м/с, что пренебрежимо мало по сравнению с осевой компонентой. Подобная задача была рассмотрена в работе [16], где Vp достигала 10 м/с. Приведенное различие может быть связано с расхождениями
ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение” 2015. № 5 79
Рис. 4. Изменение угла а поворота и угловой скорости ш во времени t:
•, х, о — численные расчеты;-кинематические параметры после завершения
процесса формирования ВСЭ
в алгоритмах определения радиальной скорости, а именно в определении области, для которой проводится расчет кинематических параметров.
Более существенно аг влияет на угол поворота а оси ВСЭ в плоскости полета, образованной взаимно перпендикулярными векторами о и р, а также на угловую скорость ш этого поворота. На рис. 4 показано изменение а и ш во времени t для различных значений асимметрии а\. На начальном этапе (от инициирования ВУ до t = 0,1 м/с) ввиду значительного деформирования облицовки угловые характеристики не применимы. На последующих стадиях формирования элемент приобретает компактную форму и продолжает деформироваться до тех пор, пока ВСЭ не будет окончательно сформирован. На протяжении всего процесса формирования угол ВСЭ постоянно увеличивается. Также заметны флуктуации значений угловой скорости, имеющие характер затухающих колебаний. Следует отметить наличие а и ш у ВСЭ при а± = 0, что обусловлено погрешностями численных расчетов. Коэффициент корреляции между а и а\ составляет r = 0,996, а между установившимися значениями ш и аг — r = 0,995.
Расчеты показали, что выбранный метод не привносит искусственных искажений в форму в осесимметричном случае (при отсутствии несоосности). Увеличение аг приводит к увеличению угла поворота ВСЭ и незначительному нарушению его осевой симметрии.
Суперпозиция гармоник погрешности. Влияние суперпозиции гармоник оценивалось аналогично влиянию несоосности, однако несоосность вводилась не в идеальную сетку, как ранее, а в сетку с периодической неравномерностью толщины облицовки в окружном направлении (n = 6). Гармоника, задающая неравномерность толщины, не влияет на симметрию движения ВСЭ, поэтому оценку влияния целесообразно проводить, варьируя а\ при фиксированной амплитуде а6 = а6/80 неравномерности толщины (1). Параметр а6 принимался равным 0,03 в целях гарантированого образования складок [3].
80 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2015. №5
Рис. 5. Влияние суперпозиции гармоник разностенности на а и ш:
I — только несоосность; II — Аф = п/3; III — Аф = 0; •, х, □ — численные расчеты; ------ регрессии
Очевидно также, что эффект от суперпозиции гармоник зависит от разности фаз Аф = |фх — ф61 между ними.
Расчеты не показали существенных различий во влиянии суперпозиции гармоник на осевую V0 и радиальную Vp скорости ВСЭ по сравнению с несоосностью наружной и внутренней поверхностей.
На рис. 5 (левые графики) приведено сравнение влияния несоосности на а и и в условиях наличия одной и двух гармоник с различной разностью фаз Аф. Характер зависимости для всех случаев имеет линейный вид (г > 0,99 для каждой отдельной конфигурации), но конфигурации с суперпозицией гармоник более чувствительны к уровню несоосности ах. Коэффициент корреляции для конфигураций с суперпозицией гармоник вне зависимости от разности фаз составляет г = 0,988 (см. рис. 5, правые графики).
Контуры элементов, сформированных из облицовок без погрешностей и с погрешностями различной конфигурации, приведены на рис. 6; ВСЭ, сформированный из облицовки идеальной формы (К1), имеет ярко выраженную осевую симметрию. Наличие неравномерности толщины облицовки (п = 6) приводит к образованию складок (К2) при этом изменения угла а поворота элемента не происходит. При суперпозиции гармоник п =1 и п = 6 (К3 и К4) с различной разностью фаз Аф наблюдается наклон элемента с искажением формы его кормовой части.
Полученные данные свидетельствуют о том, что наличие нескольких гармонических составляющих погрешностей повышает чувствительность ВСЭ к их влиянию.
ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение” 2015. №5 81
Рис. 6. Зависимость влияния погрешностей различных конфигураций на форму ВСЭ (слева направо):
K1 — идеальная осесимметричная облицовка; K2 — n = 6; K3 — n =1 и 6 при Аф = п/3; K4 — n = 1 и 6 при Аф = 0
Вследствие критического состояния материала в процессе формирования наиболее вероятным результатом наличия погрешностей значительной амплитуды является искажение формы и разрушение ВСЭ. Например, при а = 0,05 и а6 = 0,03 на промежуточных этапах формирования ВСЭ наблюдался его изгиб, приводивший в последствие к разрушению ВСЭ (рис. 7).
Выводы. 1. Предложен подход к численной оценке влияния погрешностей облицовки на кинематические параметры ВСЭ, основанный на внедрении неравномерности малой амплитуды в дискретизацию модели ВУ. Поскольку общепринятый метод расчета задач формирования ВСЭ с помощью подхода Эйлера не позволяет точно описать границу тел, то был применен подход Лагранжа. При этом параметры решателя были специально адаптированы для задач с высокой скоростью деформаций, в частности, применялся алгоритм удаления сильно искаженных элементов с критически малым характерным размером. В рамках исследований разработана математическая модель генерирования сетки облицовки с контролируемыми погрешностями несоосности и разнотолщинности. Разработан и реализован метод расчета параметров формы и кинематики ВСЭ.
X^—ir Рис. 7. Изгиб ВСЭ при наличии двух гармонических составляющих разностенно-сти облицовки
82 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2015. №5
2. Проведенные расчеты показали:
• влияние несоосности наружной и внутренней поверхностей облицовки (гармоники n = 1 разностенности) на осевую скорость ВСЭ незначительно;
• влияние гармоники n = 1 на радиальную скорость имеет линейный характер, однако оно в целом незначительно и им можно пренебречь;
• гармоника n =1 разностенности облицовки оказывает значительное влияние на угол поворота оси ВСЭ и его угловую скорость, например, при амплитуде гармоники n = 1 в 2 % толщины облицовки угловая скорость достигает 280 с-1;
• при суперпозиции гармоник n = 1 и n > 1 чувствительность ВСЭ к влиянию погрешностей ВУ возрастает. Например, при амплитуде гармоники n =1 в 2 % толщины облицовки и n = 6 в 3 % толщины облицовки угловая скорость ВСЭ достигает
420... 480 с-1 в зависимости от разности фаз между гармониками разностенности облицовки.
ЛИТЕРАТУРА
1. Takanao Saiki, Hirotaka Sawada, Chisato Okamoto, Hajime Yano, Yasuhiko Takagi, Yasuhiro Akahoshi, Makoto Yoshikawa. Small carry-on impactor of Hayabusa 2 mission, Acta Astronautica. March-April 2013. Vol. 84. P. 227-236.
2. Hutchinson J.W., Koiter W. T. Postbuckling theory // Applied Mechanics Reviews. 1970. Р 1353-1366.
3. Асмоловский Н.А., Баскаков В.Д., Тарасов В.А. Анализ влияния периодических возмущений на формирование высокоскоростных стержневых элементов // Известия высших учебных заведений. Машиностроение. 2013. № 8. С. 8-14.
4. Баскаков В.Д., Тарасов В.А., Колпаков В.И., Софьин А.С. Методика оценки влияния погрешностей конструкции снарядоформирующих зарядов на кучность попадания в цель и пробивное действие удлиненных поражающих элементов // Оборонная техника. 2010. № 1-2. С. 90-97.
5. Колпаков В.И., Баскаков В.Д., Шикунов Н.В. Математическое моделирование функционирования снарядоформирующих зарядов с учетом технологических асимметрий // Оборонная техника. 2010. № 1-2. С. 82-89.
6. Колпаков В.И. Математическое моделирование функционирования взрывных устройств // Электронный журнал “Наука и образование: электронное научнотехническое издание” МГТУ им. Н.Э. Баумана. 2012. № 02 (77-30569/334177).
7. Herve Couque and Remi Boulanger. EFP Simulations with Johnson-Cook Models. 23rd International Symposium on ballistics Tarragona. Spain 16-20 April 2007. Vol. I. P. 255-262.
8. Johnson G.R., Stryk R.A. Some considerations for 3D EFP computations // International Journal of Impact Engineering. October 2006. Vol. 32. Iss. 10. P. 16211634.
9. Johnson G.R., Stryk R.A. Symmetric contact and sliding interface algorithms for intense impulsive loading. Comput. Methods Appl. Mech. Eng. 2001. Vol. 190. Iss. 35-36. Р 4531-4549.
ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение” 2015. № 5 83
10. Beissel S.R., Johnson G.R. Large-deformation triangular and tetrahedral element formulations for unstructured meshes. Comput. Methods Appl. Mech. Eng. 2000. Vol. 18. Iss. 3-4. Р. 469-482.
11. Lou Jianfeng, Hong Tao, Liang Longhe, Han Bing. Numerical simulation of formation of EFP with charge of aluminized high explosive // 23rd International Symposium on ballistics Tarragona. Spain 16-20 April 2007. Vol. II. Р. 1265-1271.
12. Jian-qing Liu, Wen-bin Gu, Ming Lu, Hao-ming Xu, Shuang-zhang Wu. Formation of explosively formed penetrator with fins and its flight characteristics. Defence Technology. June 2014. Vol. 10. Iss. 2. P. 119-123.
13. Asmolovskiy N., TkachukA., BischoffM.Numerical approaches to stability analysis of cylindrical composite shells based on load imperfection // Engineering Computations. 2015. Vol. 32. Iss. 2.
14. Скворцов А.В. Триангуляция Делоне и ее применение. Томск: Изд-во Томского университета, 2002. 128 с.
15. Bentley J.L. Multidimensional binary search trees used for associative searching. Communications of the ACM. 1975. Vol. 18. Iss. 9. Р. 509-517.
16. Колпаков В.И., Баскаков В.Д., Кружков О.А., Шикунов Н.В. Оценка влияния технологических факторов на кинематические параметры удлиненного поражающего элемента кумулятивного заряда // Экстремальные состояния вещества. Детонация. Ударные волны: Труды междунар. конф. IX Харитоновские тематические научные чтения. 2007. С. 585-590.
REFERENCES
[1] Takanao Saiki, Hirotaka Sawada, Chisato Okamoto, Hajime Yano, Yasuhiko Takagi, Yasuhiro Akahoshi, Makoto Yoshikawa. Small carry-on impactor of Hayabusa 2 mission. Acta Astronautica, March-April 2013, vol. 84, pp. 227-236, IsSn 00945765.
[2] Hutchinson J.W., Koiter W.T. Postbuckling theory. Applied Mechanics Reviews, 1970, pp. 1353-1366.
[3] Asmolovskiy N.A., Baskakov V.D., Tarasov V.A. The Impact of Periodic Disturbanceson the Formation of High-Speed Rod Elements. Izv. Vyssh. Uchebn. Zaved., Mashinostr. [Proc. Univ., Mech. Eng.], 2013, no. 8, pp. 8-14 (in Russ.).
[4] Baskakov V.D., Tarasov V.A., Kolpakov V.I., Sof’in A.S. Assessment of the EFP Technological Imperfections on the Accuracy and Penetration of Elongated Projectiles. Oboronnaya tehnika [Defense Technology], 2010, no. 1-2, pp. 90-97 (in Russ.).
[5] Kolpakov V.I., Baskakov V.D., Shikunov N.V. Mathematical modeling of the functioning of shellmounting charges taking into account technological asymmetries. Oboronnaya tekhnika [Defense Technology], 2010, no. 1-2, pp. 82-89 (in Russ.).
[6] Kolpakov V.I. Mathematical simulation of the explosive devices’performance. Jelektr. Nauchno-Tehn. Izd “Nauka i obrazovanie” [El. Sc.-Tech. Publ. Science and Education], 2012, no. 2 (in Russ.). Availlable at: http://technomag.edu.ru/doc/334177.html
[7] Couque H., Boulanger R. EFP Simulations with Johnson-Cook Models, 23rd International Symposium on ballistics Tarragona. Spain, 16-20 April 2007, vol. I, pp. 255-262.
[8] Johnson G.R., Stryk R.A. Some considerations for 3D EFP computations, International J. of Impact Engineering, 2006, vol. 32, iss. 10, pp. 1621-1634, ISSN 0734-743X.
[9] Johnson G.R., Stryk R.A. Symmetric contact and sliding interface algorithms for intense impulsive loading. Comput. Methods Appl. Mech. Eng., 2001, vol. 190, iss. 35-36, pp. 4531-4549.
84 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2015. №5
[10] Beissel S.R., Johnson G.R. Large-deformation triangular and tetrahedral element formulations for unstructured meshes. Comput. Methods Appl. Mech. Eng., 2000, vol. 187, iss. 3-4, pp. 469-482.
[11] Jianfeng Lou, Tao Hong, Longhe Liang, Bing Han. Numerical simulation of formation of EFP with charge of aluminized high explosive. 23rd International Symposium on Ballistics Tarragona. Spain, 16-20 April, 2007, vol. II, pp. 12651271.
[12] Jian-qing Liu, Wen-bin Gu, Ming Lu, Hao-ming Xu, Shuang-zhang Wu. Formation of explosively formed penetrator with fins and its flight characteristics. Defence Technology, vol. 10, iss. 2, June 2014, pp. 119-123. ISSN 2214-9147.
[13] Asmolovskiy N., Tkachuk A., Bischoff M. Numerical approaches to stability analysis of cylin-drical composite shells based on load imperfections. Engineering Computations, 2014, vol. 32, iss. 2 (in press).
[14] Skvortsov A.V. Triangulyatsiya Delone i ee primenenie [Delaunay triangulation and its application]. Tomsk, Tomsk. Univer. Publ., 2002. 128 p.
[15] Bentley J.L. Multidimensional binary search trees used for associative searching. Communications of the ACM, 1975, vol. 18, iss. 9, pp. 509-517.
[16] Kolpakov V.I., Baskakov V.D., Kruzhkov O.A., Shikunov N.V. Assessing the impact of technological factors on the kinematic parameters of the elongated striking element of the shaped charge. Ekstremal'nye sostoyaniya veshchestva. Detonatsiya. Udarnye volny. Tr. Mezhdunar. Konf. 9 Kharitonovskie tematicheskie nauch. chteniya [Extreme states of matter. Detonation. Shock waves. Proc. of the International Conf. 9 Kharitonov thematic scientific reading], 2007, pp. 585-590 (in Russ.).
Статья поступила в редакцию 16.03.2015
Асмоловский Николай Александрович — аспирант кафедры “Технологии ракетнокосмического машиностроения” МГТУ им. Н.Э. Баумана.
МГТУ им. Н.Э. Баумана, Российская Федерация, 105005, Москва, 2-я Бауманская ул., д. 5.
Asmolovsky N.A. — Ph.D. student, Department of Technologies of Space-Rocket Mechanical Engineering, Bauman Moscow State Technical University.
Bauman Moscow State Technical University, 2-ya Baumanskaya ul. 5, Moscow, 105005 Russian Federation.
Баскаков Владимир Дмитриевич — д-р техн. наук, доцент, профессор кафедры “Технологии ракетно-космического машиностроения” МГТУ им. Н.Э. Баумана.
МГТУ им. Н.Э. Баумана, Российская Федерация, 105005, Москва, 2-я Бауманская ул., д. 5.
Baskakov V.D. — D.Sc. (Eng.), Professor, Department of Technologies of Space-Rocket Mechanical Engineering, Bauman Moscow State Technical University.
Bauman Moscow State Technical University, 2-ya Baumanskaya ul. 5, Moscow, 105005 Russian Federation.
Зарубина Ольга Васильевна — канд. техн. наук, доцент кафедры “Технологии ракетно-космического машиностроения” МГТУ им. Н.Э. Баумана.
МГТУ им. Н.Э. Баумана, Российская Федерация, 105005, Москва, 2-я Бауманская ул., д. 5.
Zarubina O.V. — Ph.D. (Eng.), Associate Professor, Department of Technologies of Space-Rocket Mechanical Engineering, Bauman Moscow State Technical University.
Bauman Moscow State Technical University, 2-ya Baumanskaya ul. 5, Moscow, 105005 Russian Federation.
ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение” 2015. № 5 85
Просьба ссылаться на эту статью следующим образом:
Асмоловский Н.А., Баскаков В.Д., Зарубина О.В. Анализ влияния технологических погрешностей менисковых облицовок на динамику взрывного формирования высокоскоростных стержневых элементов // Вестник МГТУ им. Н.Э. Баумана. Сер. Машиностроение. 2015. № 5. C. 72-86.
Please cite this article in English as:
Asmolovsky N.A., Baskakov V.D., Zarubina O.V. Research into the effect of technological errors of meniscus coatings on explosive formation dynamics of highspeed rod-shaped elements. Vestn. Mosk. Gos. Tekh. Univ. im. N.E. Baumana [Herald of the Bauman Moscow State Tech. Univ., Mech. Eng.], 2015, no. 5, pp. 72-86.
В Издательстве МГТУ им. Н.Э. Баумана вышла в свет книга
ФИЗИКО-МАТЕМАТИЧЕСКИЕ МОДЕЛИ И МЕТОДЫ РАСЧЕТА ВОЗДЕЙСТВИЯ МОЩНЫХ ЛАЗЕРНЫХ И ПЛАЗМЕННЫХ ИМПУЛЬСОВ НА КОНДЕНСИРОВАННЫЕ
И ГАЗОВЫЕ СРЕДЫ
Книга посвящена изучению физических процессов в веществе при взаимодействии с мощными источниками импульсного нагрева. Дан краткий исторический очерк, описана физика лазерного термоядерного синтеза. Основное внимание уделено построению математических моделей и численным исследованиям физических явлений в высокотемпературной плазме для инерционного и магнитноинерциального термоядерного синтеза. Приведены различные методы численных решений уравнений магнитно-радиационной плазмодинамики.
Для аспирантов физикотехнических специальностей университетов, а также научных сотрудников и инженеров, работающих в области астрофизики, систем управляемого термоядерного синтеза и физики газовых разрядов.
86 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2015. №5
