Методы увеличения скорости снаряда
Поможем в ✍️ написании учебной работы
Поможем с курсовой, контрольной, дипломной, рефератом, отчетом по практике, научно-исследовательской и любой другой работой

Итак , имеем  или   (9.30)

Из полученного выражения (9.30) следует ,чтобы увеличить скорость снаряда необходимо :

1. увеличить среднее давление пороховых газов -Рср

2. увеличить длину ствола- lD

3. уменьшить поперечную нагрузку снаряда -

4. уменьшить .

Наличие пути снаряда –lD увеличивает использование энергии газов ( ) и поэтому являются одним из возможных методов повышения скорости снаряда , однако этот способ по двум причинам ограничен ,во-первых с увеличением lD убывает Рсри с некоторой величины lD дальнейшее увеличение будет мало влиять на произведение РсрlD ,а во-вторых ,орудийный ствол должен быть такой длины ,чтобы он не прогибался под собственным весом и удовлетворял другим техническим требованиям.

Увеличение среднего давления при заданном наибольшем давлении- Рm и заданном Wкн ,можно добиться разными путями .

1. применением прогрессивных форм пороха

2. применением флегматизированных порохов

3. применением цилиндро-конических каналов

4. применением камор с большим уширением

5.  многокамерных орудий

6. увеличение относительного веса заряда 

На рис. …… показаны различные кривые р(t) при одинаковых максимальных давлениях .

 Кривая 1-соответствует мгновенному сгоранию пороха .

 

Кривая 2- соответствует горению прогрессивного пороха

 Кривая 3-соответствует дегрессивного пороха

 Кривая 4-соответствует горению "толстого" пороха.

 Из этих кривых –кривая 2 имеет наибольшее значение среднего давления .

При заданном наибольшем давлении с возрастанием прогрессивности возрастает также плотность заряжания , достигая . В то время как для дегрессивных порохов .

Плотность заряжания.

Влияние на величину скорости незначительно, флегматизация способствует дополнительному увеличению прогрессивности ,т.е. будет увеличивать среднее давление .

Применение конического ствола увеличивает среднее значение –Рср ,что видно из расчетов проверенных И.П.Граве для винтовки .

Результаты представлены в таблице:

Таблица 19

R канал
2380 834 0,441 0,489 0,235 цилиндрический
2375 860 0,376 0,522 0,303 цилиндро-конический

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

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

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

 наибольшего давления будет возрастать и дульная скорость снаряда обычно R изменяется в пределах 0,4-0,6

 (9.31)

И, наконец , уменьшение коэффициента фиктивности - приводит к увеличению скорости снаряда . Насколько это возможно –рассмотрим подробнее.

9.4 Коэффициент учета второстепенных работ – коэффициент фиктивности - .Учет теплоотдачи и прорыва пороховых газов .

Рассмотрим виды работ совершаемые пороховыми газами :

 1. L1 - работа поступательного движение снаряда – главная работа .

 (9.32)

2. L2- работа, затрачиваемая на вращение снаряда   (9.33) ,

где

-радиус инерции снаряда , определяемый по формуле , где I- момент инерции относительно оси вращения

r=d/2- радиус сечения снаряда

-угол нарезки

- крутизна нарезов

 h- шаг нарезки

обычно  от 0,25% до 2,5%

К2=0,0025-0,025

 

 

L3- работа на преодоление трения между пояском снаряда и внутренней поверхности канала ствола , а также на преодоление трением между центрующими утолщениями снаряда и полями нарезов  (9.34) ,где 1 -коэффициент трения

4. L4- работа затрачиваемая на перемещение газов самого заряда и несгоревшего пороха . В предположении о постоянстве плотности  газопороховой смеси по всему за снарядному пространству .т.е. плотность  зависит только от времени и не зависит от координаты сечения и в самом сечении , энергия , затрачиваемая на перемещение заряда , будет выражаться зависимостью :   (9.35) , где , , где l0- приведенная длина каморы =

                                   L0- действительная длина каморы

                                   l- путь пройденный снарядом

Работа L4 вычисленная по формуле (9.35) будет больше действительной работы , т.к. плотность газопороховой смеси уменьшается в направлении к снаряду , что доказывают газодинамические расчеты .

5. Работа затрачиваемая на перемещение откатных частей.

(9.36)

, где  -масса откатных частей ; Q0 –вес откатных частей .

(9.37) ,где V –скорость откатных частей

 mг=  -масса заряда, т.к.   и , то

L6 –работа , расходуемая на врезание ведущего пояска в нарезы . Как правило не учитывается и может быть учтено косвенно в момент вылета снаряда из ствола .

L7 –работа, расходуемая на преодоление снарядом сопротивления воздуха , находящегося в канале орудия . Этой работой при малых скоростях пренебрегают . Однако при скоростях 2,5-3  противодавление уже составляет 150-300  и учет противодавления необходим .

- тепловая энергия , расходуемая во время выстрела на нагрев стенок ствола , гильзы и снаряда – потеря на теплоотдачу . Учитывается путем уменьшения силы пороха ,-f.

(9.38)

в момент времени t при t=0 ,

при t=tд

 ,где

; ; ; ; ;

; ; ; ; ; ;

 l –путь снаряда

-скорость снаряда

См% -определяется по кривой С% от tк при кг/дм2

d –калибр снаряда

 lкм –длина каморы

 Wкн –объём канала ,включая камору

 - вес заряда

q –вес снаряда

 Рсн – давление на снаряд

Таким образом , сила пороха по мере движения снаряда по каналу за счет теплоотдачи убывает ; она зависит от :

1) калибра-d

2) плотности заряжания -

3) коэффициента уширения каморы -

4) относительного пути снаряда -

5) отношения скоростей

 При выводе использовалось предположение , сделанное Мюреуром :

и соотношение  , сделанное профессором Вентцелем.

- энергия , теряемая газами , прорывающимися по зазорам между пояском снаряда и стенками канала орудия.

(9.39)

 где ; - относительная часть прорвавшихся газов

 

Если расход газа невелик , то , или с учетом теплоотдачи .

Значение  ,где a=1.09-1.03 и b= - для цилиндрических каналов и   и - для конических каналов.

9.5 Анализ изменения давления пороховых газов в канале ствола от условий заряжания .

Имея формулу для давления из основного уравнения пиродинамики (9.19)

,

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

(9.40)

,  или , , ,

Как видно из выражения (9.40) нарастание давления зависит от многих факторов .

В момент формообразования Р=Р0

, нарастание давления зависит от Р0,f и ,и обратно пропорционально  , при Р0=0 (для миномета ) ,тангенс угла наклона будет равен 0 . И далее тангенс угла наклона возрастает до точки перегиба и далее тангенс угла наклона убывает до 0 и далее становится отрицательным за счет значения .

На рис. ….. фиг. 91 показан характер этих кривых .При  получим соотношение , в момент ….. горения при ,  

 При переходе ко второму периоду выражение для давления имеет вид    , , ,

 Характер нарастания давления в функции от пути –l выразится общей формулой   (9.41)

Вначале движения ,когда тангенс угла наклона равен ∞ ,т.е. кривая Р(l) будет иметь касательную совпадающую с осью ординат рис.(фиг.92)

  

9.6 Влияние формы и размеров пороха на кривые давления газов и скорости снаряда.

Анализ формул (9.40) и (9.41) показывает ,что характер нарастания давления как во времени так и функции от пути снаряда зависит , главным образом от  при данной "силе" и природе пороха зависит от  

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

  форма зерна
1 трубка 1,003 0,0994 0,997
2 лента 1,06 0,89 0,943
3 пластинка 1,20 0,675 0,810
4 брусок 0 0 ~2,0 0 0
5 куб 3 0 0 3 0 0

 Нанеся на график изменение в зависимости от Z получим диаграмму изображенную на рис(фиг.93). На рис(фиг.94) приведены расчеты давлений газов при  и в функцию от пути . Диаграмма показывает ,что лента дает нормальное давление  и дульную скорость

.Брусок –4 , имея большую начальную оголенность дает давление  и . Куб –5 в следствии втрое большей оголенности имеет давление  и . Если бы поставить задачу : сравнить дульные скорости при одинаковых Рm ,то ленточный порох показал наилучшие результаты .

Теперь рассмотрим влияние толщины свода при одинаковой форме зерна . Результаты расчетов сведены в таблицу :

2
1,5 1,414 1,256 3540 632
2,0 1,06 0,943 2040 575
2,5 0,848 0,744 1450 486

Влияние толщины пороха на кривые давления показаны на рис…(фиг.96)

 

10. Решение основной задачи внутренней баллистики (ОЗВБ).

 Установление закономерностей , связывающих разнообразные условия заряжания с зависящими от них величин , называемыми баллистическими элементами выстрела составляет общую задачу внутренней баллистики .

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

 К баллистическим элементам выстрела относятся : изменяющееся во времени путь снаряда –l, скорость снаряда , давление пороховых газов –Р, их температура – Т , а также количество газов , образовавшиеся к данному моменту ; а также относительная толщина горящего свода –Z.

При решении указанной выше ОЗВБ можно выделить две важнейших основных задач пиродинамики и ряд частных задач .

Первая основная задача пиродинамики состоит в определении расчетом изменения газов и скорости снаряда в канале ствола в функции от пути снаряда и от времени при заданных условиях заряжания . При этом наряду с кривыми Р(l),υ(l) или P(t),υ(t) и l(t) определяются две важнейшие баллистические характеристики орудия – наибольшее давление газов –Рm в канале ствола и дульная скорость снаряда -  ,т.е. скорость снарыда при вылете его из канала ствола . Эту задачу называют прямой задачей пиродинамики . При заданных условиях заряжания она иееет единственное решение. Изменяя условия заряжания можно провести анализ этих условий на изменение кривых давления газов и скорости снаряда ,т.е. решить ряд частных задач . Точность решения этой задачи зависит от выбранной математической модели выстрела и методов решения . Вторая основная задача пиродинамики – задача баллистического проектирования орудия состоит в определении конструктивных данных канала ствола и условий заряжания , при которых снаряд данного калибра-d и веса-q , получает при вылете определенную дульную скорость - .Эта скорость задается на основе тактико-технических требований ,предъявляемых к проектируемому орудию . При решении её обычно , задаются наибольшим давлением газов –Рm. Решение этой задачи многовариантно от целесообразности и рациональности выбранного варианта баллистического решения в значительной степени зависит дальнейшее проектирование всей артиллерийской системы в целом и боеприпасов к ней . По выбранным условиям заряжания производится расчет кривых давления и скорости . Полученная кривая Р(t) или P(l) используется конструкторами для расчета прочности стенок орудия и снаряда , лафета ,дистанционных трубок , взрывателей . Вместе с этим даются требуемая толщина и форма пороха , который должен быть изготовлен на заводе .

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

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

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

Основные допущения при решении ОЗВБ:

1. Горение пороха подчиняется геометрическому закону горения или физическому закону горения.

2. Порох горит при средних давлениях p , воспламенение мгновенное.

3. Состав продуктов горения не меняется (f и  - постоянные).

4. Скорость горения пороха пропорциональна давлению .

5. Учитываемые второстепенные работы пропорциональны главной работе поступательного движения снаряда и учитываются при помощи коэффициента .

6. Движение снаряда начинается , когда в каморе в результате сгорания части заряда разовьется давление форсирования –p0 , постепенность врезания в n не учитывается.

7. Работа врезания пояска отдельно не учитывается.

8. Растяжением стенок ствола при выстреле , прорывом газов через зазоры между ведущим пояском и стенками канала ствола и сопротивлением воздуха в канале ствола пренебрегаем.

9. Охлаждение газов в результате теплоотдачи стенкам ствола непосредственно не учитывается и может принято в расчет косвенно , снижением f и увеличением .

10. Движение снаряда рассматривается до момента прохождения его дна через дульный срез.

11. Величину  принимаем равной среднему значению для всего периода выстрела.

12. Плотность газопороховой смеси зависит только от времени и не зависит от координаты.

 Таким образом уравнения классической внутренней баллистики для усредненных значений давлений p, температуры T и относительного количества сгоревшего заряда  , при этом осреднение T и  получается как следствие осреднения давления . Иначе говоря , в классическом методе внутренней баллистики волновые процессы течения газа не учитываются , и применяется "термодинамический" закон расширения газов .

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

 Для нахождения элементов выстрела в классическом методе О.З.В.Б. имеем следующие зависимости :

1)  - основное уравнение пиродинамики , уравнение Резаля .

2)  или  - закон горения пороха .

3)  - двухчленный закон газообразования .

4)  - закон движения снаряда .

5)  - кинематическая связь между скоростью и путем .

Совокупность этих 5-ти уравнений позволяет найти 5-ть неизвестных p,z,υ,l,  как функции времени.

10.1 Система уравнений ОЗВБ для пороха простой дигрессивной формы.

 Для вычисления элементов выстрела по имеющейся специальной программе, по которой решается система обыкновенных дифференциальных уравнений 1-го порядка на ПЭВМ (macad) приведем наши уравнения к нормализованному виду (система 10,1)

(10,1)

где l  =  ;  .

Система уравнений ( 10,1 ) решается при следующих начальных данных :

при t = 0 , z =0 ,  = 0 , p = pв , υ = 0 , l = 0 , конец расчета l = lд .

 

Точность и правильность расчета проверяются по аналитическим зависимостям в момент p = p0 :

При горении пороха :

 

где pm , zm , Vm – значения давления , скорости и относительной толщины в момент максимального давления .

где υпр =

Правильность решения проверяется при различных условиях заряжания в постоянных переменных Н.Ф. Дроздова :

 

 

10.2 Системы уравнений для многоканального пороха прогрессивной формы .

 При решении ОЗВБ для многоканального пороха воспользуемся двухчленными формулами для 1-й фазы горения ( до распада зерна ) и 2-й фазы горения ( после распада зерна ) . Для 1-й фазы горения характеристики

 и  найдутся из условия , что при z = 1  и при z = 0,5 значения  будут совпадать как по трехчленной , так и по двухчленной формулам :

при z=0,5  

при z=1  

откуда имеем

 ,

 Для второй фазы горения имеем

 

 характеристики  и  найдутся из условия , что при z=zk  и при

z=zk поверхность горения должна обратиться в ноль . Откуда получим :

 

Решая находим эти уравнения:  

Подсчитав χ и λ по этим формулам для 7-ми канального стандартного пороха имеем :

=0,712 ; =0,225 ; = - 0,0237 ; =0,855 ; =1,375 ; ; =0,1873 ;

ошибка =0,004 .

Принимая для 7-ми канального пороха наружный радиус вписанного круга в наружную призмочку 0,532е1 получим ек1+0,532е1=1,532е1

 ; =0,855 ; = - 0,94 ;

Для зерна Уолша с 7-ю каналами  ; =0,95 ; = - 2,16 ;

=1,37 ; ; =0,218 .

Исходя из 2-х фаз горения : прогрессивного в первой фазе и дегрессивного во второй фазе система уравнений ( 10,2 ) будет иметь следующий вид

 

(10.2)

где    ;

Система ( 10,2 ) решается при следующих начальных данных : при t = 0 , z=0, =0 , p = pв , l = 0 , υ= 0 . Расчет заканчивается при l = lд . Точность и правильность расчета проверяется как при расчете системы ( 10,1 ) .

10.3 Системы уравнений для комбинированного заряда из дегрессивных порохов простой формы .

 Пусть комбинированный заряд состоит из n порохов дегрессивной простой формы . Характеристики i-того пороха обозначим с индексом i , так что

 , где i =1...n .

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

  ;   ;   ;   ;  

Расставив пороха по импульсу в конце горения , пороха по возрастающему значению J1<Ji<Jn и воспользовавшись системой уравнений ( 10,1 ) окончательно получим

 

( 10,3 )

где

Точность решения системы уравнений проверяется по аналитическим зависимостям , по которым проверяется точность системы ( 10,1 ) . При этом необходимо брать характеристики эквивалентного фиктивного пороха f , Jk ,

 , Г ,  ,  , формулы которых приведены выше . Значения  и  такого пороха определяются по формулам ( 10,4 )

          ( 10,4 )

Уравнение газообразования примет вид :

( 10,5 )

10.4 Системы уравнений для комбинированного заряда , состоящего из 7-ми канального и трубчатого порохов .

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

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

 Пусть трубчаты порох имеет индекс – Т . 7-ми канальный имеет обозначения такие же как в системе ( 10,2 ) . ОЗВБ для такого заряда решается по системе уравнений ( 10,6 ) :

  (10,6 )

Фиктивный порох эквивалентный комбинированному заряду имеет следующие баллистические характеристики :

1)  при  где

2) или  при  где

 

Точность решения системы проверяется по аналитическим зависимостям , представленным выше , которые справедливы в 1-м случае – до распада зерна 7-ми канального пороха , во втором случае – до конца горения трубчатого пороха , который наступает раньше распада зерна .

10.5 Исходная система уравнений внутренней баллистики для миномётов, орудий и ракет.

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

За основу исследований при выводе системы уравнений принимаем частично уравновешенное орудие рис... , для которого справедливо соотношение :

( 2,13 )

где n – коэффициент уравновешенности ; S – площадь поперечного сечения канала ствола с учетом нарезов ; p – баллистическое давление ( среднее давление газов в за снарядном пространстве в данный момент времени ) ; Gp1- расход газов через сопловой блок орудия , имеющий размерность "кг/с" ;

Jr – удельный импульс , развиваемый пороховыми газами при истечении из сопла .

Величина коэффициента уравновешенности "n" ограничена пределами 0<=n<=1 . Для классического артиллерийского орудия n = 0 , для безоткатного n =1 .

 Будем полагать так же , что имеет место прорыв пороховых газов через ведущее устройство снаряда , количественно характеризующееся расходом

Gp2 , а так же имеется теплоотдача стенкам канала ствола , и стенка ствола расширяется упруго при выстреле .

Исходную систему уравнений запишем при следующих допущениях :

1. Горение пороха происходит параллельными слоями , т.е. справедливо уравнение газоприхода  и относительную поверхность горения :  где

 - характеристики формы пороха ; z =  - относительная толщина сгоревшего пороха ; e1 – половина толщины порохового зерна ; e – толщина слоя сгоревшего пороха ;  - сгоревшая часть порохового заряда ;  - вес заряда ; Sгор – горящая поверхность заряда ; Sгоро – начальная поверхность заряда .

2. Давление p , температура Т и плотность газопороховой смеси  в заснарядном пространстве для каждого момента времени t равны их среднему по объему значениям ( гипотеза квазистационарного процесса ) .

p , T и  связаны уравнением состояния :

где R – газовая постоянная ;  - коволюм газа .

3. Состав продуктов сгорания не меняется во время выстрела , а удельные теплоемкости Cp , C  равны их средним значениям для всего диапазона изменения температур .

 и = const .

4. Воспламенение порохового заряда происходит мгновенно .

5. Отсутствует выброс несгоревших частиц пороха .

6. Противодавлением воздуха в канале ствола пренебрегаем .

При выводе системы уравнений используем основные законы термодинамики :

закон сохранения энергии – первого закона термодинамики запишем в виде :

здесь  - скорость изменения тепла в газе , вес которого , к рассматриваемому моменту времени составит  , вследствие его взаимодействия с окружающей средой .

- скорость изменения внутренней энергии газа , где U= (2,15)

- мощность , развиваемая газом при его расширении или при сжатии .

Применительно к периоду движения снаряда при горящем заряде :

 ( 2,16 )

Здесь  - скорость подвода тепла вследствие сгорания порохового зерна , где Е = 4270  - механический эквивалент тепла;

 - калорийность пороха при воде жидкой , т.к.  .

- приход продуктов сгорания .

 - скорость оттока тепла из каморы орудия в атмосферу вследствие истечения газов через сопло ( Gp1 ) и в зазоры между ведущими устройствами снаряда и стенками канала ствола ( Gp2 ) .

- энтальпия одного кг газа .

 - скорость изменения тепла в следствие теплоотдачи между стенками ствола и газами ( символ  показывает , что не является полным дифференциалом ) .

Скорость отвода тепла из-за снарядного пространства в стенке канала ствола:

( 2,17 )

где  - коэффициент теплоотдачи от газа к стенкам ; Тсг – температура внутренней поверхности ствола ; F – поверхность , омываемая газами .

Точно уравнение ( 2,17 ) решается совместно с уравнением теплопроводности материала стенки при соответствующих краевых условиях. Величина  в общем случае выражается уравнением :

 ( 2,18 )

Здесь W – свободный объем за снарядного пространства ;

 - мощность , создаваемая газом , вследствие сгорания заряда с учетом истечения части газа из за снарядного пространства , где - удельный вес пороха ;

 - мощность , затрачиваемая на упругие деформации стенок , гильзы и ствола , где  - "упругое" приращение за снарядного объема . Обычно этой мощностью пренебрегают , хотя в пушках высоких давлений она может являться ощутимой ( ППН,ЛГУ ) .

( 2,19 )

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

Обычно в расчетах  принимается постоянным :

 ( 2,20 )

где к зависит от калибра и начальной скорости снаряда и типа снаряда .

С уменьшением калибра снаряда величина "к" – возрастает и уменьшается с ростом начальной скорости . Для снаряда с ведущим пояском "к" = 1,02-1,05.

Для пуль , имеющих калибр меньше 14,5 мм и не имеющих ведущего пояска

"к"=1,2-1,3 .

С учетом выражений ( 2,15 ) , ( 2,16 ) , ( 2,18 ) уравнение сохранения энергии применительно к рассматриваемой схеме ( рис... ) приводится к виду :

 

 ( 2,21 )

Соотношение ( 2,21 ) так же называют основным уравнением внутренней баллистики . Уравнение сохранения вещества может быть записано в виде

( 2,22 )

Приход газов в следствии сгорания пороха  определяется по формуле

 ( 2,23 )

Значения Gp2 , Gp1 зависят от значений полной температуры T00 и полного давления газов p00 в за снарядном пространстве , а так же площадей критических сечений сопла   и зазор Fз . При расчете Gp2 следует дополнительно учитывать скорости снаряда и формы зазоров во времени

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

или   ( 2,25 )

Уравнение состояния :

где

где

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

 

(2,26)

Накладывая определенные ограничения , с помощью полученной системы уравнений (2.26) можно описать процессы выстрела в следующих системах:

1. Безоткатные орудия (n=1).

2. Миномет (n=0 , Gp1=0).

3. Классическое орудие (n=0 , Gp1=0 , Gp2=0).

4. Ракетный двигатель (Gp2=0 , υ=0).

5. Бомба постоянного объема ( Gp1=0 , Gp2=0 , υ=0).

 Если параметры состояния газов определяются после окончательного горения заряда в системе (2.26) необходимо положить Sгор=0 и U=0.

После вылета снаряда из канала ствола расчет продолжается при Sгор=0 и U=0 и Fз=S и из системы (2.26) исключается уравнение движения снаряда.

  Внутренняя баллистика классического орудия.

Для закона сохранения энергии , когда Gp1=0 , Gp2=0   в момент горения пороха будем иметь:

 (2.27).

В предварительном периоде горения пороха идет при υ=0 и l=0 до момента p=p0 и  , где p0- давление форсирования снаряда;

 количество газа , образовавшегося в момент t=t0 – форсирования снаряда. Интегрируя уравнение (2.27) получим:

 ( 2,28 )

На 2-м периоде – периоде расширения пороховых газов ( Sгор = 0 , U = 0 ,

) , уравнение сохранения энергии примет следующий вид :

 

Для предварительного периода ( υ = 0 ) получим

 Для приближенного учета теплоотдачи  воспользуемся допущением , согласно которому процесс теплоотдачи можно считать квазистационарным,

с коэффициентом теплоотдачи – линейно зависящим от удалённого веса газа в за снарядном пространстве

 ( 2,29 )

где  - постоянный коэффициент теплоотдачи . Т.к.  , то

 =

 

Введем обозначение  , тогда

 =  ( 2,30 )

 Значение  находятся из экспериментов . Для автоматических пушек допустимо принимать   , . С учетом этого допущения система ( 2,26 ) примет вид :

( 2,31 )

где Jk =  

Учитывая , что  

где f – сила пороха , - доля твердых остатков в продуктах сгорания ( для дымных порохов  = 0,5 ; для бездымных порохов  = 0 ).

 

 

 

 

( 2,32 )

где  

т.к.  и

 то получим

( 2,33 )

Подставляя выражение ( 2,33 ) в первое уравнение системы ( 2,32 ) окончательно получим :

( 2,34 )

 

где F = S0 +

 ;

Если пренебречь растяжением стенок ствола (  ) окончательно получим систему уравнений с учетом теплоотдач :

 

 

( 2,35 )

где  ;

 .

Рассмотрим учет теплоотдачи при выстреле предложенный Мюрауром для бомбы , принимая во внимание постепенное возрастание охлаждающей поверхности стенки :

( 2,36 )

где Cm – экспериментально найденный коэффициент в бомбе при сжигании дымного пороха по времени сгорания его при  = 0,2 кг/дм3 .

Имеется кривая Cm=Cm(tk) или таблица , по которой Cm можно определить .

С другой стороны имеем :

 ( 2,37 )

Согласно формуле профессора Мамонтова М.А.

( 2,38 )

где  - коэффициент теплопередачи для газов при U = 0 .

   - плотность пороховых газов

      

 

  U – скорость течения газов у стенки

   n – показатель степени : n = 0,5-1

   - скоростной коэффициент , определяемый из опытов .

Сравнивая ( 2,36 ) и ( 2,37 ) получим :

( 2,39 )

Коэффициент Cm/ учитывает скорость течения пороховых газов в орудии .

В мfнометрической бомбе U=0 и Cm/ = Cm ,тогда

 ( 2,41 )

где - потеря температуры пороховых газов за счет теплоотдачи ;

   T1 – температура горения пороха ;

  7,774 – переходной коэффициент от бомбы Мюраура к нашим условиям .

Поверхность теплоотдачи F =  ( 2,42 )

где F0 =  - поверхность каморы орудия ;

  lкам – длина каморы ;

  Д – диаметр каморы ;

  d – диаметр канала ствола ;

  l – пройденный путь снарядом в канале ствола ;

 - коэффициент , учитывающий поверхность граней нарезов ;

tn – глубина нарезов ;

  n – число нарезов .

Объем каморы

где

- приведенная длина каморы ;

S- площадь сечения ствола ;

 - уширение каморы .

Учитывая , что сила пороха f = RT1 , то введя значение f1 =  (2,43)

мы придем к системе уравнений имеющей вид ( 11,1 ) принимает вид системы ( 2,44 )

 

 

 

( 2,44 )

 где

;

    

W0 =

;   ;

Начальные условия :

при t = 0 , p = pв , V = 0 , l = 0 ,  = 0 , z = 0 .

Расчет ведется до l = lд . Шаг интегрирования не более

 и уточняется в процессе расчета .

Входными данными являются :

Параметры ствола – Д , lкам , d , tn , a , n , lд .

Параметры снаряда – q , p0 .

Параметры порохового заряда -  .

Марка пороха : 2е1 ,  или Jк ,  .

Коэффициент фиктивности  или k :

10.6 Баллистическое проектирование артиллерийских стволов .

Задача обычно расчленяется на две части :

1. Устанавливается калибр орудия , тип снаряда и его начальная скорость , обеспечивающая решение поставленной боевой задачи .

2. Определяются размеры канала ствола и характеристики заряда .

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

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

10.5.1. Порядок баллистического проектирования артиллерийского ствола .

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

2. Определяются исходные данные для баллистического решения , с этой целью вычисляется дульная энергия снаряда – Ед :

Коэффициент могущества – CE :

CE  =

причем величина дульной скорости принимается равной начальной скорости Vд = V0 , т.е. с небольшим запасом на импульс после действия пороховых газов . Величина CЕ определяет в первом приближении исходные данные системы . Профессор В.Е. Слухоцкий , проанализировав большое число артиллерийских систем установил приближенные зависимости между коэффициентом могущества CЕ и основными характеристиками орудия . Результаты этого исследования приведены в таблице .

 

CЕ ,  , pm ,  , lд , клб. CЕ ,  , pm ,  , lд , клб
100 124 1700 0,5 1,02 14 900 106 3250 0,69 1,85 71
200 120 1950 0,55 1,09 23 1000 105 3350 0,69 1,98 78
300 117 2200 0,59 1,18 31 1100 104 3450 0,70 2,11 85
400 114 2400 0,62 1,28 38 1200 104 3550 0,71 2,25 91
500 112 2600 0,64 1,39 44 1300 103 3650 0,71 2,40 98
600 110 2800 0,66 1,50 51 1400 103 3750 0,72 2,57 105
700 108 2950 0,67 1,61 57 1500 102 3900 0,73 2,75 112
800 107 3100 0,68 1,73 64 1600 102 4000 0,74 2,95 119

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

В связи с повышением качества ствольной стали установление с помощью таблицы ( ) значение pmкр следует увеличить примерно на 10-15% .

3. Установление изменения исходных данных .

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

Рекомендуемые значения плотностей заряжания  М.Е. Серебряковым и М.С. Гороховым приведены в таблице .

пороха

pm кг/см2

 ,

кг/дм3

  2000 2400 2800 3200 3600 4000
лента 0,53 0,6 0,66 0,71 0,76 0,8 1,1
трубка 0,54 0,62 0,68 0,73 0,77 0,82 1,0
7-ми кан. 0,66 0,74 0,78 0,84 0,88 0,93 1,35

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

CЕ , тм/дм3 100-1000 1200 1400 1600
86 85 84 83

Орудие наименьшего объема оказывается реально приемлемым лишь при

V0 > 1300 м/с . При меньших начальных скоростях следует отойти от него , увеличивая  , т.е. увеличивая длину ствола . Таким образом  берется в пределах , указанных в таблицах .

Зная  нетрудно получить все исходные параметры по формулам :

,   ,   ,   ,   ,

,  .

Профессор М.Е. Серебряков для установления вариантов решения предлагает использовать специальный график  называемый директивной диаграммой рис... Центром графика является точка с координатами   и  , соответствующего орудия наименьшего объема .

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

При проектировании систем малого калибра , следует обратить внимание на то , что эти системы имеют большую плотность заряжания  = 0,8-0,9 кг/дм3 , при которых сокращаются габариты и вес патрона , что облегчает проектирование механизмов автоматики . Такие плотности заряжения обеспечиваются применением мелких трубчатых или 7-ми канальных или сферических порохов .

4. Расчет вариантов .

В каждом варианте необходимо определить длину нарезной части и полную длину ствола с затвором , а так же примерное число выстрелов за время баллистической жизни ствола . Первая часть легко решается с помощью ТБР по величинам pm таб и Vд таб при данном  ( или решается система уравнений 10,1 или 10,2 ) .

lств=lд+lкам+(1-2)d

Для оценки живучести можно использовать упрощенную формулу профессора В.Е. Слухоцкого :

, где

Nусл – живучесть ствола ; A – постоянный коэффициент для всех вариантов ;

где tдтаб – время выстрела ;

5. Сравнение вариантов и выбор лучшего .

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

где M = const , кратное десяти .

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

6. Уточненный расчет выбранного варианта и построение кривых

p(l) , V(l) , t(l) .

7. Рассчет ствола на прочность .

Для этого найти кривую наибольших давлений на систему ствола , построенную при 3-х температурах заряда t0max заряда , t0нор.зар. и t0min.зар. .

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

и

 

11. Газодинамический метод решения ОЗВБ.

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

Т.к. волны , равносильно малым возмущениям , распространяется как вправо , так и влево с местной скоростью звука – c по отношению подвижного наблюдателя , двигающегося со скоростью движения газа – U (вправо) . Закон движения волн , идущих направо , относительно неподвижного наблюдателя ( например стенки ствола ) будет :

,

а скорость движения волн , идущих налево будет :

.

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

 1-й случай .Скорость снаряда стала на столько велика и ствол очень короткий , что отраженная волна от дна каморы не успевает догнать снаряд и отраженная волна , дойдя до дульного среза констатирует , что снаряд вылетел из ствола .

 2-й случай . Скорость снаряда стала велика , но ствол длинный , то отраженная волна успевает догнать снаряд и отразившись от дна не может двигаться влево по отношению неподвижного наблюдателя ( U>c ) и тем самым не сможет достигнуть дна каморы второй раз .

 3-й случай . Скорость снаряда мала , а ствол достаточно длинен , что отраженная от дна каморы волна догоняет снаряд , отразившись от дна снаряда успевает достигнуть дна каморы , снова отразившись от дна каморы может второй раз догнать снаряд .Такое многократное отражение волны , когда скорость снаряда U<c приводит к нивелированию характеристик течения газа , когда можно будет считать , что средняя плотность газа в за снарядном пространстве будет падать обратно пропорционально объему этого пространства , а давление с плотностью будет связано законом изэнторы ( адиабаты ) , т.е. установится "термодинамический" режим расширения продуктов сгорания . Этот случай и описывают уравнения систем ( 10,1 ; 10,2 ; 10,3 ; 10,6 ) при решении ОЗВБ при допущении что параметры газа ( p , T ,  ) в за снарядном пространстве не зависят от координаты сечения , т.е. являются функциями только времени и одинаковы по всему за снарядном пространстве .

В отличии от "термодинамического" режима течения газа в газодинамическом методе решения ОЗВБ это допущение снимается , т.е. параметра газа ( p , T ,  ) и скорость течения – U и скорость звука в газе – с зависят не только от времени , но и от координаты за снарядного пространства . Строгая постановка этой задачи приводит к осемметричному движению газопороховой смеси . Однако , учитывая , что объем переходного конуса по сравнению с объемом цилиндрической части каморы и цилиндро–конического канала ствола и толщина пограничного слоя пороховых газов , прилегающих к стенкам каморы и канала ствола , мала по сравнению с диаметрами . Решение можно искать в одномерной постановке, т.е. предполагая , что параметры газа ( p , T , , U , c ) одинаковы в каком-то сечении и равны параметрам на оси канала ствола и каморы . Решение задачи Лагранжа , выполненное кандидатом физико-математических наук Зинченко Ю.К. в осемметричной постановке показало , что при  < 15 0 (угол конусности канала ствола) приводит к разнице результатов не более 10% в случаи решения этой задачи в одномерной постановке .

11.1 Основные уравнения одномерного движения газа и их решение.

Выделим в газе , находящемся в трубе , некоторый объем v . Сила , действующая на него равна - . В случае одномерного движения

- = -

где p – давление ; f – площадь поверхности , ограничивающей объем ; dx – элемент длины ; - элемент объема .

Уравнение движения будет иметь вид :

где  - плотность газа , U – скорость газа ,  - ускорение газа .

Производная  определяет не изменение скорости в данной неподвижной точке пространства , а изменение скорости данной частицы газа передвигающейся в пространстве , т.е. :

где  , t – время , x – координата .

Таким образом приходим к уравнению движения в форме Эйлера :

 

Рассмотрим закон сохранения массы газа , в нашем объеме она равна

. Полная масса газа , вытекающего ( или втекающего ) через поверхность f за единицу времени будет равна:

С другой стороны , уменьшение ( или увеличение ) массы газа в объеме за единицу времени будет - . Приравнивая оба выражения имеем :

Это равенство справедливо для любого объема . Поэтому

( 11,2 )

Уравнение ( 11,2 ) называется уравнением неразрывности . Для изоэнтропического ( адиабатического ) течения газа уравнение энергии записывается в виде идеального газа :

 ( 11,3 )

Для реального газа , для которого справедливо уравнение состояния в виде

      (пороховой газ).

уравнение энергии будет иметь вид :

 ( 11,4 )

где - коволюм единицы массы газа .

 Для изоэнтропического течения преобразуем уравнение движения и уравнение неразрывности , выражая p и  через теплосодержание i (энтальпию ) или скорость звука в газе с . Для этого воспользуемся следующими термодинамическими соотношениями :

     

А случае изотропы :

     

Для пороховых газов :

        ( 11,5 )

 Отметим: скорость звука реального газа выше скорости звука идеального газа ( ).

Кроме с - скорости звука С.А. Бетехтеным был введен аналог скорости звука

     

 или  (11.6). Из соотношения ( 11,6 ) следует, что аналог скорости звука совпадает со скоростью звука в идеальный газе . Разрешая соотношения (11,5) , (11,4) , (11,6) относительно p и  получим :

                  (11.7)

        

Если при помощи соотношений ( 11,6 ) и ( 11,7 ) исключить p и , то уравнения ( 11,1 ) и ( 11,2 ) будут иметь вид :

Для пороховых газов :

 

 ( 11,8) - уравнение движения .

 ( 11,9)- уравнение непрерывности .

Для идеального газа ( ) .

 (11,10) и

 ( 11,11 )

Складывая или вычитая почленно уравнение (11,8) и (11,9) или (11.10) и (11.11)получим окончательно :

 ( 11,12 ) для реального газа .

 (11,13) для идеального газа .

В таком виде уравнение ( 11,12 ) или ( 11,13 ) удобны для исследования . Очевидно , что заданные значения величины  или  распространяются вправо со скоростью U+c , а заданное значение

 или  распространяется влево со скоростью U-c .

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

Система уравнений ( 11,12 ) или ( 11,13 ) есть система уравнений в частных производных гиперболического тела . Методов решения их достаточно много. Мы остановимся на методе характеристик , который более наглядно дает физическую картину течения газа .

Решение задачи Лагранжа .

 Задача Лагранжа в классическом ее понимании заключается в следующем.

В цилиндре , ограниченном слева неподвижной стенкой и безгранично продолженном вправо на расстоянии l0 от дна цилиндра находится поршень весом q . Площадь поперечного сечения цилиндра равна S0 . В пространстве между дном цилиндра и поршнем находится  кг газа под давлением p0 с плотностью 0 . Газ однороден и неподвижен . В момент t = 0 поршень получает возможность двигаться без сопротивления под действием давления газа . Движение газа является одномерным , влияние теплоотдачи не учитывается . Требуется определить возможное в этих условиях движение газа и движение поршня. Таким образом задача Лагранжа сводится к основной задачи внутренней баллистики в предположении о мгновенном сгорании заряда и отсутствии теплоотдачи сопротивления движению и других второстепенных работ , за исключением движения газа . Задача поставленная в 1790г. неоднократно привлекала к себе внимание многих ученых : Пиддек ( 1921г. ) , Госсо и Лиувиль(1922г. ) , Фок( 1935г. ) , Платрие ( 1936г. )-вот далеко не полный перечень ученых , посветивших свои работы решению этой задачи . Наиболее полно она решена С.А. Бетехтиным (1948г.) с учетом коволюма газа для каморы с уширением , с учетом теплоотдачи . Позднее она решалась Л.Л. Поповым , Зайченко Ю.И. ( в оссиметричной постановке ) , Никулиным О.А. в новых относительных переменных , Ушаковым В.М. , Комаровским А.В. и другими .

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

С газодинамической точки зрения задача Лагранжа является задачей об одномерном неустановившемся движении газа при соответствующих начальных и граничных условиях . Если поместить начало координат у дна цилиндра и направить ось X в сторону движения поршня , то начальными условиями будут следующие : при t =0 0=<x<l0 , Р=Р0 , , U=0 , T=T0 .

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

X=0 , U=0

X=l+l0 ,    ( 11,14 )

где l - путь , пройденный поршнем к рассматриваемому моменту времени .

Допустим , что газ идеальный . Рассмотрим область движения в плоскости X,t . Слева она ограничена осью t , а справа - неизвестной еще нам кривой , изображающей закон движения поршня X=f(t) . Начальными условиями определены все искомые функции на отрезке 0 - l0 оси X . Известно , что f'(0)=0 , f''(0)>0 , f''(t)>0 , т.е. поршень начинает двигаться с нулевой скоростью , но с конечным ускорением , причем ускорение поршня остается положительным во время его движения . Поскольку слой газа , примыкающий к дну поршня не может ни оторваться , ни обогнать поршень в своем движении , будет всегда выполняться условие :

Ux=f(t)=f'(t)

где Ux=f(t) - скорость слоев газов , прилегающих к дну поршня . Таким образом скорость , а следовательно и остальные параметры переднего слоя будут непрерывно изменяться в результате движения поршня .

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

(11.15) , где С –скорость звука в газе (уравнение 11.15 нельзя записать в конечном виде т.к. "U" и "c" могут быть функциями "x" и"t" ).

 Допустим ,что нам удалось найти частное решение системы (11.13), т.е. отыскать такие две функции :  и  , которые удовлетворяют уравнениям (11.13) и граничному условию , имея это решение , мы можем построить в плоскости (x,t) систему линий

, которые называются характеристиками (рис… фиг.2).

В условиях задачи Лагранжа , линии будут прямыми , идущими расходящимися пучками , т.е. на одно элементарное возмущение не будет в процессе своего перемещения обгонять предыдущее ,т.к. скорость звука падает при расширении . А так как каждое последующее элементарное возмущение распостраняется по газу , все более и более разряженному , воздействием предыдущих возмущений , то отсюда и следует , что одно последующее возмущение не сможет догнать предыдущее . Система значений U=0 и с=с0 нетрудно убедится , также будет одним из частных решений системы U=U1(x,t) и c=c1(x,t) будут справедливы лишь выше линии ОА . Таким образом вдоль линии ОА происходит переход от одного частного решения системы к другому .

Как было показано выше вдоль характеристик выполняются соотношения , аналогично вдоль характеристики , идущей вправо  будет выполнятся соотношение или окончательно можно записать :  (11.16)

( характеристики I семейства

(11.17)

( характеристики II семейства

 S и R называются инвариантами Римана .

Уравнения (11.16) и (11.17) имеют вполне определенный физический смысл. Уравнения характеристик в плоскости x, t как следует выше описывают законы перемещения элементарных возмущений . Уравнения же характеристик в плоскости искомых функций определяют связь между изменениями скорости течения "U" и местной скорости звука –"с" , которая должна выполнятся при перемещении данного элементарного возмущения .

Значения "U" и "с" в любой точке области движения можно выразить через S и R в этой точке , т.е. через значения постоянных интегрирующих вдоль характеристик , противоположных семейств , проходящих через данную точку ,   (11.18)

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

1. В некоторой части плоскости (x, t) S имеет на всех характеристиках I семейства одно и то же значение S=Idcm, а R=Idm. одно и то же значение на всех характеристиках II семейства в этой же части плоскости ,т.е. во всех точках плоскости U=const и c=const , т.е. поток газа будет в этой части плоскости однородным . Если при этом S=R,то U=0 ,а с=с0 рассматриваемая область будет областью покоя ( рис…)

2. Постоянные интегрирования одного из семейств ( предположим ,что I семейства ) сохраняет постоянное значение в некоторой области (S=Idm.) постоянные же интегрирования другого семейства (II семейства ) R меняется от характеристики к характеристике. Поэтому в пределах рассматриваемой области вдоль каждой характеристики II семейства в плоскости x, t скорость течения "U" и местная скорость звука будут сохранять постоянные значения (S=Idm. R=const вдоль этой характеристики ) I семейства , то сами характеристики II семейства будут прямыми линиями. Течения такого рода носят названия волн одного направления (рис..) .Так как вдоль линии ОА все характеристики I семейства имеют одно и то же значение S=Idm, то в области волны одного направления характеристики II семейства будут прямыми линиями , а характеристики I семейства, где в каждой точке только S=const ,а R –меняется ,эти характеристики I семейства будут кривыми линиями . Характеристики I семейства –кривая АВ будет отделять область волны одного направления от области , где S, R меняться от характеристики к характеристики ( область общего случая неустановившегося движения в газе ).

3.  Третий случай , когда в пределах рассматриваемой области оказываются переменными как S так и R это является наиболее общим случаем неустановившегося движения газа . Для этого случая нельзя установить каких-либо иных закономерностей , кроме тех, которые заложены в самих уравнениях характеристик.

В области волны одного направления параметры газа и скорость течения , а также закон движения поршня находятся по аналитическим зависимостям . В общим случае , задачи решаются численно . Если ОА разбить на 256 отрезков , то значения  находятся с точностью до 0,1% в области волны одного направления , эту точность можно перенести и в общий случай . Как правило , чтобы уменьшить количество вариантов задача Лагранжа решается в относительных переменных . Обычно за отнсительные переменные принимаются :

, , , , , , (11.19)

Такие относительные переменные неудобны при оценке влияния коволюма газа , ускорения поршня длины ствола , уширения каморы .

В частности ,в этих относительных переменных решена задача Лагранжа ,С.А. Бетехтиным ,Н.Н. Поповым и другими . Из результатов решения трудно выявить волновую картину процесса расширения газа .

О.А. Никулиным предложены новые относительные переменные , которые свободны от перечисленных выше недостатков .

За единицу измерения параметров принимаются : , где - аналог скорости звука в газе , равный

(11.20)

- начальное ускорение поршня (11.21)

S0 –площадь сечения канала ствола

q –вес снаряда

Р0 –начальное (максимальное давление газа)

Т0 –начальная температура(max)

R –газовая постоянная

- ковалюм единицы массы газа

-начальная плотность газа (max)

K –показатель адиабаты

 Тогда , , , , , , , (11.22)

В области волны одного направления относительные переменные определяются по следующим формулам ( в зависимости от )

(11.23)

                  (11.24)

(11.25)

            (11.26)

             (11.27)

                (11.28)        (11.29)

 Уравнение справедливо до момента т. В ,когда отраженная от дна канала ствола волна догонит снаряд

              (11.30)

 Откуда оптимальный относительный вес газа

(11.31) , т.е. когда отраженная от дна каморы волна догонит снаряд у дульного среза .

Из выражения (11.20-11.31) следует ,что скорость снаряда  и параметры газа , а также оптимальное значение  при одних и тех же значениях  и длине ствола S0 не зависят от ковалюма газов влияет только на объём каморы , или длину каморы при S0=const увеличивая её на величину  . Численные расчеты показывают справедливость такого вывода и в более общем случае .

газ водород

Таблица 25

 

1 0 0 1 0 0
0,7881 0,1671 0,00157 0,7882 0,1671 0,0016
0,6351 0,3137 0,00650 0,6352 0,3137 0,0065
0,5242 0,4440 0,0150 0,5214 0,4440 0,0152
0,4344 0,5614 0,0276 0,4346 0,5614 0,0275
0,3668 0,6673 0,0437 0,3669 0,6669 0,0437
0,3134 0,7638 0,0641 0,3135 0,7634 0,0640
0,2704 0,8519 0,0885 0,3705 0,8519 0,0885
0,2553 0,9336 0,1172 0,2354 0,9332 0,1172
0,0881 1,1314 0,2606 0,0883 1,1311 0,2602
0,0361 1,262 0,5385 0,0361 1,2613 0,5373
0,01558 1,354 1,0472 0,0155 1,3533 1,0439
0,00689* 1,421* 1,9520* 0,0069* 1,4205* 1,9456*
0,003116 1,472 3,548 0,0031 1,4718 3,5328
0,001427 1,512 6,3531 0,0014 1,5113 6,3222

 

В таблице отмечены изменения давления на снаряд , относительная скорость снаряда и относительный путь снаряда для идеального водорода ( ) и реального водорода ( ). Из таблицы видно ,что коволюм не влияет на эти характеристики и так как расчеты проводились в абсолютных значениях величин , а обработка в новых относительных переменных , то видно ,что точность расчета выше чем 0,1% (при сравнении с аналитическими формулами (11.20)-(11.31)).

На рис…показано изменение относительного давления от скорости снаряда , а на рис….. изменение относительной скорости от относительного пути снаряда . Из рис. 1 и 2 видно ,что огибающие кривых соответствуют случаю когда  ∞ , точками 0.5;1.0;1.5;2.0;2.5;3.0;3.5;4.0;5.0;6.0;7.0 соответствует значению  ,   и т.д. и они показывают :

1. при заданной относительной длине ствола , соответствует максимальной относительной скорости снаряда и оптимальное значение веса газа , т.е.  и

2. Если мы увеличим количество газа , то скорость снаряда ни найоту не увеличивается ,т.е. дополнительное количество газа просто не участвует в передачи энергии снаряду т.к. дополнительного слоя волна разряжения идущая от дна снаряда не успела дойти , а снаряд вылетел из ствола .

 Если мы уменьшим количество газа т.е. , то получим резкое падение скорости снаряда при заданной длине ствола

Например, при  , , ,  , ,

при  и , , , то есть ни дульная скорость снаряда, ни дульное давление не увеличивалось.

Увеличилась только длина каморы в 1,5 раза . Этот вывод не был бы сделан если бы не делали пересчет при использовании прежних относительных переменных , при  при  , которое не зависит от (это важно ,при обычных относительных переменных  увеличилась бы в 2 раза ) дульная скорость , .

На рис. 3 приведены оптимальные значения для разных газов , принимая

 показатель адиабаты для гелия  , для водорода и пороховых газов . Уравнения (11.23),(11.26),(11.31)примут вид :

для гелия

для водорода

для пороховых газов

 Результаты расчетов представлены в таблице 26

  5/3 7/5 11/9 5/3 7/5 11/9 5/3 7/5 11/9
0,5 0,2316 0,206 0,191 0,550 0,434 0,364 0,4996 0,4997 0,4998
1,0 2,016 1,469 1,225 1,56 1,11 0,902 0,9949 0,9956 0,9987
1,5 12,75 6,50 4,615 3,75 2,24 1,64 1,484 1,476 1,484
2,0 102 25,8 14,5 10 4,24 2,76 1,884 1,93 1,95
2,5 2269 107,5 42,34 43,8 8,17 4,5 2,19 2,333 2,386
3,0 - 529,8 122,0 - 17,0 7,25 - 2,67 2,78
3,5     358,9     11,8     3,075
4,0     1108     19,7     3,42

В таблице приведены результаты скорости снаряда  в случае использования классического метода расчета –" термодинамического " расширения газа .

,где  или окончательно

(11.32)

при одинаковых значениях  и .

Из таблицы в частности видно ,что при дульных скоростях снаряда для пороховых газов  .

Термодинамическое расширение даёт дульную скорость на 500-600  ниже , а при скорости 3,5-3,8  реально получены на установках ППН для легких снарядов ( ) .Из соотношения (11.32) ,найдем оптимальное значение , которое дает максимальное значение скорости при заданной длине ствола .

Дифференцируя и приравнивая , получим соотношение :

(11.33)

 Разрешая это уравнение , находим .

 Например . Для пороховых газов  получим вместо  значение  при этом , т.е. 2.5% больше , что не может быть т.к. в передаче энергии снаряду участвует лишь  или пусть  получим вместо  значение  при этом  вместо , однако количество газа участвующая в передаче энергии снаряду равно . Остальное количество газа лишнее .

11.3 Влияние уширения каморы.

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

Результаты расчетов при больших значениях и резком переходе из газовой каморы l ствол приведены на рис.3 .

Приближенно оценить увеличение скорости снаряда за счет уширения каморы можно по выведенным О.А. Никулиным формулам

(11.34)

 увеличение инвариантно S

где "U1", и "U" скорости газа на входе и выходе из соединительного конуса между каморой стволом .

 Формула (11.34) выведена при предположении ,что в переходном конусе .

(11.35)

 Подставляя (11.35) в уравнении (11.34) окончательно получим

(11.36)

где  

Для бесконечного уширения каморы x=∞ получим

(11.37)

Результаты полученные для водорода и гелия при  и разных значениях уширения каморы приведены в таблице 27

 

Таблица 27

водород

гелий

приближенное точное приближенное точное
4 0,388 0,302 0,347 0,288
16 0,587 0,402 0,570 0,388
0,605   0,645  

Точное решение полученное при решении задачи Лагранжа с уширением каморы и 16 при нахождении точного решения радиуса переходного конуса имеем вид   (11.38)

 где -высота переходного конуса R и r радиусы каморы и ствола , соответственно . В таблице 28 приведены результаты расчетов для реального водорода ( ) с уширением каморы x=16 и разных значениях

Таблица 28

1 0,1921 0,9239 0,1358 0,0563 0,019
2 0,3757 0,8232 0,238 0,1370 0,081
3 0,5469 0,7180 0,260 0,2871 0,187
4 0,7046 0,6212 0,3190 0,3854 0,346
10 1,1032 0,3854 0,4684 0,6663 1,206
20 1,6860 0,1668 0,5574 1,1286 5,03
32 2,1028 0,0800 0,5889 1,5139 13,684
40 2,207 0,0544 0,595 1,702 21,80
64 73,065 0,0074 0,5458 2,519 144,17

Как видно из таблицы в случае реального водорода (с учетом ковалюма ) значение  приближается к 0,595 т.е. к вычисленной по приближенной формуле 0,587 для идеального водорода . Таим образом реально получить приращение 500-600м/с при одинаковых и максимальных давлениях используя уширение каморы x=16 или .

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

11.4 Условия постоянного давления на снаряд по всей длине ствола.

При горении пороха рассмотрим характеристику I семейства в области волны одного направления  откуда

 т.к. имеется зависимость  получим

 или

(11.39)

 где используется соотношение  при малых значениях  .В случае пороховых газов  ,

, т.е. ошибка равна  т.е. ошибка составляет 13% в сторону занижения .( , где относительная ошибка ). Рассмотрим пороховую камору большого объёма , когда количеством газа можно пренебречь , тогда

 (11.40)

при  получим  т.к. , то окончательно   (11.41), т.е. зная максимальное давление в каморе и давление на снаряд, которое постоянно по стволу и , легко получить марку пороха и все остальные характеристики .

 ,принимая

и К=11/9 получим  ;  при .

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

 

12.Пороховые и легко газовые установки.

 Пороховые и легко газовые установки нашли широкое применение при гиперзвуковых исследованиях .

Пороховые установки рис. 55 как и орудие имеет ствол , пороховую камору и затвор . В пороховых установках ППН , созданных под руководством О.А. Никулина : на ствол, навинчивался пороховой стакан , представляющий собой толстостенную болванку со сплошным дном . В болванке имеется камора с уширением x=3,6 переходящая в цилиндр диаметром равным диаметру канала ствола . В ствол помещалось метаемое тело в поддоне. Иногда между стволом и пороховым стаканом помещалась мембрана со специальной насечкой для увеличения давления форсирования снаряда . В дне болванки имелось отверстие , через которое вставлялся инициирующее устройство в камору . Ствол со стаканом размещался на легкой раме . При выстреле установка отскакивала назад . Калибры ствола 23;34;50. Длина ствола 70-80клб. Давление в каморе 10000-15000кг/см2. Скорости ,полученные на установках со снарядом Сq=1,0-2,0кг/дм3 были 3000-2600 м/с соответственно.

На установках использовались штатные пороха. Были проведены специальные стрельбы по выяснению влияния различных фактов влияющих на баллистику выстрела , в частности применение различных добавок и основному заряду , использование комбинированных зарядов и желеобразных топлив и т.д. Проектирование установок базировалось на результатах решения задачи Лагранжа в новых относительных переменных ,предложенных О.А. Никулиным . В дальнейшем установки ППН были модернизированы с увеличением объёма каморы , и её уширения до x=8-9 при этом достигнуты скорости  на тех же калибрах . При решении ОЗВБ использовался газодинамический метод , разработанный профессором В.М. Ушаковом.

На установках ППН проведены несколько десятков тысяч опытов . Типичный характер изменения давления от пройденного пути снарядом представлены на рис 56 легкогазовые установки.

Как известно , максимальная возможная скорость неустановившегося течения : , которая прямо пропорционально скорости звука в газе . При одинаковой температуре газа наибольшей энергией будут обладать атомарный водород , гелий , так называемые легкие газы . Скорость звука молекулярного водорода при комнатной температуре примерно 1300м/с , выше скорости звука пороховых газов ( 1000м/с) .Если нагреть легкие газы до температуры 2000-2500К , то скорость звука 3000-4000м/с, а следовательно и максимальная скорость снаряда будет значительно выше , чем она достигнута на пороховых установках . Нагреть водород или гелий можно различными способами : электрическим разрядом , ударной волной, стехиометрической смесью кислорода и водорода (геливопаровая пушка) и движущимся поршнем , который сжимает легкий газ и в процессе сжатия нагревает его до нужной температуры . Поршневые легко газовые установки оказались наиболее перспективные и на них достигнуты скорости до 12 км/с ,скорости 7-8км/с являются рабочими скоростями ЛГУ , где используется деформируемый поршень . На рис.57 показано схематически ЛГУ с легким поршнем .

ЛГУ состоит из пороховой каморы поз.1 , в которой размещен пороховой заряд и газовой каморы поз.4, в которую накачивается легкий газ ( водород или гелий )при комнатной температуре до давления 5-30кг/см2 . Пороховая камора отделена от газовой поршнем , выполненного из легко деформируемого материала (полиэтилена ). Газовая камера соединяется со стволом (поз.8) с помощью конического переходника (поз.5) , выдерживающей высокие давления газа (8000-15000кг/см2) . В ствол вставляется снаряд (поз.9) , состоящей из метаемого тела , которое помещается в поддоне из легкого прочного материала , так что относительный вес снаряда Сq=1-3кг/дм3. Метаемое тело отделено от легкого газа диафрагмой (поз.10) , которая раскрывается при определенном заданном давлении ( давления форсирования –Рф). Установка работает следующим образом : при подаче электрического импульса на инициирующее устройства, устройство срабатывает и воспламеняет пороховой заряд . При достижении давления форсирования , которое больше или равно начальному давлению водорода , поршень начинает двигаться по газовой камере . По мере сгорания порохового заряда поршень двигаясь с большой скоростью сжимает легкий газ и нагревает его : создаются определенные условия по давлению и температуре . При достижении заданного давления легкого газа –Рф диафрагма разрывается и снаряд – метаемая сборка под действием легкого газа разгоняется по вакуумированому стволу до заданной скорости . Такая легко газовая установка называется двух ступенчатой ЛГУ с "легким " поршнем . В первую ступень входит пороховая камера и поршень. Выстрел из первой ступени подобен выстрелу из обычного орудия с противодавлению снаряду . Решение задачи внутренней баллистики первой ступени и методы идентичны ОЗВБ орудия. Здесь могут быть использованы как газодинамический , так " термодинамический " методы решения . Разделение поршней на "тяжелый " и "легкий" связано с достижением поршнем скорости больше или меньше скорости звука . При малых скоростях поршня ("тяжелый" поршень ) ударной волны не образуется в легком газе и метод характеристик , изложенный выше как правило используется при решении ОЗВБ второй ступени . В случае использования "легкого " поршня газодинамический подход к решению ОЗВБ второй ступени обязательные методы "сквозного " счета , где ударная волна "размазывается " за счет введения "искусственной " вязкости газа .

Типичная картина изменения давления в пороховой и газовой камерах представлены на рис…..

Исследование внутренней баллистики ЛГУ усложняется не только за счет увеличения количества параметров , от которых зависит процесс выстрела , но и конструктивным оформлением поршня, диафрагмы , метаемой сборки , наличие конического переходника , в котором происходит торможение поршня в целом при одновременном ускорении его переднего торца за счет "гидроэффекта " . Приборы и аппаратура используемая при внутри баллистических исследованиях базируется на последних достижениях науки и техники в области быстропротекающих процессов и главное –результаты исследований , достижения в ЛГУ вполне могут быть перенесены на баллистику орудия . В этом смысле исследование в ЛГУ –это "форпост" в баллистике орудия. В доказательство этого вывода являются докторские диссертации Л.В. Комаровского , Ю.П. Хоменко , В.М. Ушакова , В.В. Жаровцева и др.

                               Список используемых источников

1. Благонравов А.А. Основы проектирования автоматического оружия

2. "Баллистические установки и их применение в экспериментальных исследованиях" под. ред. Златина Н.А. и Мишина О.А. издательство "Наука" , Москва 1974г.

3. Бетехтин С.А. , Виницкий М.А. , Горохов М.С. и др. "Газодинамические основы внутренней баллистики" . Издательство оборонной промышленности , Москва – 1957г.-384с.

4. "Бог войны" . Сборник статей . Автор-составитель полковник Латухин А.И. Издательство "Молодая гвардия" , Москва , 1879. , 256с.

5. Профессор Горохов М.С. "Внутренняя баллистика" . Отчет СФТИ при ТГУ им. В.В. Куйбышева , г. Томск , 1950г., 273с.

6. Корнер

7. Иванников Л.С. и Никулин О.А.

"Эксперементально-теоретическое определение силы сопротивления для 37ми мм. снаряда при врезании в нарезы канала ствола орудия" . Дипломная работа , ТГУ , г. Томск , 1957г. –96с.

8. Никулин О.А. "Приложение" к диссертации К.Т.Н. по спецтеме ( инв. №4034 ) , г. Бийск , 1968-40с.

9. Профессор Орлов Б.В. и Мазинг Г.Ю.

"Термодинамические и баллистические основы проектирования ракетных двигателей на твердом топливе" . Второе издание , переработанное и дополненное . Издательство "Машиностроение" , Москва ,1968г. –536с.

10. "Проектирование ракетных и ствольных систем" под редакцией д.т.н. проф. Орлова Б.В. Издательство "Машиностроение" , Москва , 1974г.-828с.

11. Академик Зельдович Я.Б. , Лейпунский О.И. , Либрович В.Б.

"Теория нестационарного горения пороха" , М. "Наука" 1975г.-131с.

12. Христенко Ю.Ф. и Жалнин Е.В.

"Монометрическая бомба" . Заявка на изобретение №2000118367

приоритет от 10.7.2000г.

13. Христенко Ю.Ф. "Экспериментальные методы исследования нестационарных эффектов при горении заряда зерненного пороха // исследование по баллистики и смежным вопросам механики" Вып.3 , Томск ,

издательство ТГУ , 1999г. , стр. 38-39 .

14. Христенко Ю.Ф. "Экспериментальные методы исследования закономерностей горения зерненных порохов в широком диапазоне изменения плотностей заряжания" . Доклады региональных конференций

Волж.Р.Ц. РАРАН "Современные методы проектирования и обработки ракетно-артиллерийского вооружения" , г. Саратов , издательство

РФЯ.Ц. - ... 2000г. с341-354 .

15. Проф. Серебряков М.Е. "Внутренняя баллистика" – М. , Оборонгиз , 1949г. , 458с.

16. Проф. Серебряков М.Е. "Внутренняя баллистика ствольных систем и пороховых ракет" 3-е издание , дополненное и переработанное ГНТИ оборонгиз , Москва 1962г. – 704с.

17. Тюлина И.А. "Жозеф Луи Лагранж" издательство "Наука" , Москва 1977г. 224с.

18. ЭРР "Артиллерия в прошлом , настоящем и будущем" . Военное издательство НК Обороны Союза ССР , Москва 1941г. – 348с.

19. Фундаментальные и прикладные проблемы современной механики (доклады 2-ой всероссийской научной конференции) . г. Томск, 6-8 июня 2000 г.

Дата: 2019-07-24, просмотров: 317.