Анализ прочности метаемых тел при импульсном нагружении

Страницы работы

5 страниц (Word-файл)

Содержание работы

 АНАЛИЗ ПРОЧНОСТИ МЕТАЕМЫХ ТЕЛ ПРИ ИМПУЛЬСНОМ НАГРУЖЕНИИ

Д.П. Александров*, В.А. Санников** (*СПб государственный технический университет, **Балтийский гос. технический университет, СПб, Россия)

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

В основу положена схема численного интегрирования одномерных дифференциальных уравнений внутренней баллистики, описывающих динамику пороховых газов в частных производных при постоянном поперечном сечении [2]:

Уравнение движения газа записывается в лагранжевой форме:

            (1)

где скорость, tвремя, т—масса на единицу площади, Р—давление и Q—искусственная вязкость для «размазыва­ния» ударных волн [I]. Уравнение сохранения энергии имеет вид

                    (2)

где ЕВвнутренняя энергия, отнесенная к единице массы газа, r-1=1/r—удельный объем (величина, обратная плот­ности) и VE скорость поступления энергии в газ за счет сгорания заряда. Движение снаряда в стволе описывается уравнением

              qд×S  .                         (3)

Здесь МТмасса снаряда, —его скорость, qддавление газа на основание снаряда и S — площадь поперечного сечения ствола.

Уравнение (3) при решении задачи определения напряженно-деформированного состояния (НДС) приобретает вид

[M(t)]d{(t)}/dt + [K(t)]{U(t)} = {F(t)}  ,         (4)

где: [М(t)], [K(t)] - матрицы массы, жесткости метаемого тела; {F(t)}vвектор столбец нагрузок вызванных внутренними массовыми силами или {F(t)}т тепловыми изменениями в объеме V;{q(t)}S1 – распределенными поверхностными нагрузками (например qд) на поверхности S1; {U}S2 –перемещениями на S2. В целом для общего вектора {F(t)}= {F(t)}v + {F(t)}т + {q(t)}S1 .

Решение уравнений (1) и (3), с учетом уравнения состояния (3), описывающих функции переменные в пространстве и вре­мени, выполняется методом конечных элементов, а (4) - методом Бубнова, индифферентного к граничным условиям, что соответствует характеру рассматриваемых процессов

 (4,а)

Векторы перемещений Ui-1искоростей  в начальный момент времени, при t=tн=0, равны значениям U(0) и (0); {F(t)}- текущий вектор нагрузки; i- шаг итерации прямого решения, a a11, a12, b11, b12, g11  и т. д. - коэффициенты рекуррентного со­отношения, учитывающие координату времени на исследуемом отрезке D(ti); [K]i, [M]i-1, Fi-1, Fi - матрицы и векторы параметров исследуемого объекта, размерностью n´n входящие в уравнение (1) на текущем и предыдущем шагах.

Для обоснования применимости МКЭ к расчету внутренней баллистики используются экспериментальные и теоретические данные, приведенные в работе [2]: длина камеры орудия Lо=100 см, площадь его поперечного сечения S=130 см2, масса снаряда МТ = 3,8×105 г, полное расстояние, проходимое снарядом, L = 600 см. Порох (90% нитроцеллюлозы и 10% добавок) и состоит из цилиндрических зерен диаметром ØD=0,69 см и плотностью твердого вещества ro =1,5 г/см3, полный вес заряда W= 9,0×103 г.

Термическое уравнение состояния газа представляется в виде

Р×(V-b)=пRТ        ,                   (5)

где b»0,95 см3/г—объемная поправка; получающееся при сгорании 1 г нитроцеллюлозы количество газа п»0,04 грам-моль/г; R газовая постоянная, равная ~83,15 бар×см3/ (граммоль×°К).

Процесс горения цилиндрических зерен пороха диаметром fD, где 0 £f£1 при сгорании пороха в окружающем газе, описывается эмпирическим законом

     ,         (6)

где b—постоянная скорость горения [4]. Полагая процесс выгорания топлива как j = 1 — f2, с учетом                       

                 ,           (7)

и получаемой скорости поступления энергии в газ равной  В=E0 зависимость плотности газа от j  будет описывается формулой

ro =Wj / [S(L+x)-(1-j)W/ro],            (9)

где х=х(j)расстояние, которое проходит снаряд под давлением.

Похожие материалы

Информация о работе