книги / Последействие газов на ствол. Расчет и моделирование дульных тормозов
.pdfростях течения, близких к звуковым, обе модели дают сравнительно малоотличающиеся результаты.
1.3. Аналитические формулы моделей последействия газов
Последействие с политропным и адиабатическим истечением газа
Очевидно, что процесс истечения газов из канала носит политропный характер с неизвестным заранее показателем п. Для получения аналитических выражений для P t , t , T t ,
Uд t ,Gд t , Rд t в периоде последействия используются из-
вестные законы термодинамического состояния газа и изменения состояний в политропном процессе:
P T – уравнение состояния идеального газа (1.5);
P |
|
|
P |
|
n |
|
n |
|
|||
n |
const |
или |
P |
|
|
|
|
, |
P |
|
, |
|
|
|
д |
|
|
д |
|
|
|
|
где п – показатель политропного процесса.
Соотношения между относительными параметрами газа в политропном процессе имеют вид
|
n |
; |
|
|
n |
|
|
|
1 |
|
|
P |
|
|
n 1 |
; |
|
n 1 |
; |
||||
|
|
|
P T |
T |
|
||||||
|
1 |
|
|
|
n 1 |
|
|
|
n 1 |
(1.9) |
|
P n ; |
T P n ; |
|
. |
||||||||
T |
|
Дифференцируя уравнение состояния (1.5) по времени, получаем
P T T или |
P |
|
T |
|
T |
|
|
|
|
|
|
. |
|
|
|
|
||||
|
P |
|
P |
|
P |
Скорость изменения плотности имеет вид (1.3) обык-
новенной дифференциальной модели и связана с расходом G через дульное отверстие:
11
Gд G.
Учитывая это и используя уравнения состояния, можно записать
P Gд G T .
P T
Соотношение TT найдем, используя выражение (1.9),
|
|
n 1 |
P |
1 |
P |
|
|
|
|
|||
T |
|
n |
n 1 |
P |
|
|||||||
|
|
|
||||||||||
|
n |
|
|
|
|
. |
||||||
|
|
|
n 1 |
|
|
|||||||
T |
|
|
|
|
|
|
|
|
|
n P |
||
|
|
P 2n |
|
|
||||||||
|
|
|
|
|
|
|
|
Далее, используя (1.8), найдем соотношение меду давлением и расходом газа:
G |
P |
|
P |
|
или G |
P |
n 1 |
|
||
|
|
|
|
|
|
. |
|
|||
|
|
|
|
2n |
(1.10) |
|||||
T |
|
n i |
|
|||||||
|
|
2n |
|
|
|
|
|
|
||
|
|
|
P |
|
|
|
|
|
|
|
Подставляя найденные соотношения в дифференциальную форму уравнения состояния, получим
|
|
|
G |
|
|
n 1 |
|
n |
|
|
|
|
|
|
|
|
G |
|
|
|
2n |
|
|
|
|
|
|
||||||||||
P |
|
|
д P |
|
|
|
1P |
|
или |
P |
n |
д |
. |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
P |
|
|
|
1 |
|
n P |
P |
3n 1 |
|
|||||||||
|
|
Pn |
|
|
|
2n |
|
|
Предлагается самостоятельно доказать, что полученные уравнения есть закон сохранения энергии (1.4) при n k и при отсутствии учета боковой теплоотдачи.
Интегрируя последнее уравнение, получаем
P |
dP |
t |
|
Gд |
|
|
|
|
|
|
||||
|
|
|
|
n |
|
|
dt, |
|
|
|
(1.11) |
|||
|
3n 1 |
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||
1 |
P 2n |
0 |
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
2n |
Gд n 1 |
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
||||
чтоприводит к решению P |
1 Bt |
n 1 , где B |
|
2 |
|
пока- |
||||||||
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
затель интенсивностиистечения газов.
12
Нетрудно найти и текущее время процесса истечения:
tn 1 Pn2n1
B
1 . (1.12)
Учитывая, что теоретическая продолжительность периода последействия соответствует падению давления до величины 2Ратм, её расчетное значение получим, используя (1.12):
|
1 2P |
|
2n |
||
|
|
|
|
|
n 1 |
tn B |
P |
|
|
||
|
|
|
атм |
|
|
|
|
|
д |
|
|
|
|
|
|
|
|
1 . (1.13)
Далее приводим выражения для плотности, температуры и расхода потока:
2 |
|
|
|||
1 Bt |
|
|
, |
(1.14) |
|
n 1 |
|||||
T 1 Bt 2 , |
(1.15) |
||||
G 1 Bt |
n 1 |
(1.16) |
|||
|
. |
||||
n 1 |
Последействие с изотермическим истечением газа
Профессор Е.Л. Бравин описывал последействие газов из предположения изотермичности процесса истечения. Зависимость, выведенная Е.Л. Бравиным,
P e bt . |
(1.17) |
Действительно, если положить n 1 (изотерма), то выражение (1.11) для политропного истечения примет экспоненциальный вид, близкий к (1.17),
P e |
Gд |
|
Gд |
|
(1.18) |
||
|
, т.е. |
в |
. |
||||
|
|||||||
|
|
|
|
|
|
Однако Е.Л. Бравин предлагает несколько отличающийся подход к определению показателя интенсивности:
13
в |
SPд |
, |
(1.19) |
|
0,5 Vд |
||||
|
|
|
где β – коэффициент полного действия газов на ствол, определяется как отношение средней эффективной скорости истечения газа из ствола за период последействия к скорости истечения
в начале периода, т. е. β = Ucp .
Uд
Для адиабатического процесса истечения
|
2 |
|
2 1 |
|
. |
(1.20) |
||
|
|
|
|
|
|
|||
k |
|
k 1 M |
|
|||||
|
|
д |
|
Однако на практике часто пользуются эмпирической формулой, так называемой французской формулой:
|
1300 приV 1000м/с; |
|
|
д |
|
|
Vд |
|
|
1400 приVд 1000м/с. |
(1.21) |
|
Vд |
|
Продолжительность периода последействия при изотермическом истечении
|
1 |
2Pатм |
|
|
tn |
в ln |
|
. |
(1.22) |
P |
||||
|
|
д |
|
Сравнение аналитических моделей последействия по интенсивности истечения газов показано на рис. 1.2. Расчеты проведены для 130-миллиметрового орудия с 15,5кг,
Vд 1030м/с, |
Pд 123МПа. При этом Gд 2590 кг/с в соот- |
ветствии с (1.6).
14
Рис. 1.2. Изменение давления P t в канале ствола в период
последействия при различных процессах: 1 – политропный процесс, n 1,45, tn 0,0167 с;
2 – адиабатический процесс, n k 1,25, tn 0,0267 с; 3 – изотермический процесс, n 1, tn 0,0384 с;
4 – формула E.Л. Бравина, tn 0,0534 с
Графики показывают, что теплоотдача от порохового газа стенкам ствола ускоряет интенсивность падения давления. Например, отличия в давлениях, посчитанных по политропной зависимости (1.11) и формуле E.Л. Бравина, достигают 60–70 %.
1.4.Измерение количества движения газов в стволе
иреакция истечения в дульном отверстии
Ранее была получена связь между силой отдачи, реакцией газов в дульном отверстии и изменением их количества движения в периоде последействия (1.2). Получим аналитическое выражение для текущего количества движения L и скорости изме-
нения dLdt . Введем основные допущения:
1.Плотность газа по всей длине канала ствола в каждый момент времени есть величина постоянная, т.е. ddx 0.
2.Скорость фиксированного слоя газа в момент начала истечения Ux (tд) сохраняет своё значение при движении этого
слоя к дульному отверстию. Следовательно, скорость истечения данного слоя в текущий момент Ue t Ux tд (рис. 1.3).
15
Рис. 1.3. Расчетная схема, поясняющая допущения аналитической модели
3. Распределение скоростей течения газа по длине канала всегда подчиняется линейному закону, согласно гипотезе Пиобера (см. рис. 1.3).
Ux t Ue t |
x |
Ux t Ue t |
x |
|
e |
e |
|||
|
|
|||
|
с |
|
с |
В результате получаем соотношение
Ux t Ue t exс
Очевидно, что масса газа гх , находящегося левее выде-
ленного слоя, в соответствии с принятыми допущениями также линейно зависит от координаты х и остается неизменной.
При движении к дульному отверстию
г t гх tд , |
гх t г t |
x |
, |
гх t г tд |
x |
, |
|
e |
e |
||||||
|
|
|
|
|
|||
|
|
с |
|
|
с |
|
г t г tд exс exс .
Определим начальное количество движения газов в стволе в момент вылета снаряда:
16
Lд |
1 |
U tд или |
Lд |
1 |
Uд. |
(1.23) |
|
2 |
|
|
2 |
|
|
В текущий момент времени
L t 12 г t U t
или с учетом ранее полученных соотношений
L |
U |
|
|
x 2 |
|
U |
2 |
|
|
U |
|
P |
. |
||||
|
д |
|
|
|
д |
|
|
|
|
|
|
д |
|||||
|
|
|
|
|
|
|
|
|
г |
|
|
|
|
|
|
2 |
|
|
2 |
|
|
eс |
|
|
2 |
|
|
|
|
2 |
|
|
|
||
Учитывая (1.23), (1.9), имеем |
|
|
|
|
|
|
|
|
|
||||||||
|
L Lд 2, |
|
или |
|
L Lд |
2 |
. |
|
|
|
(1.24) |
||||||
|
|
|
n |
|
|
|
Выражения (1.23), (1.24) определяют общий характер измерения количества движения газов в стволе в зависимости от текущей плотности или давления.
Политропное и адиабатическое истечение
Используя зависимость (1.11), получим
|
|
|
|
|
|
|
4 |
|
|
|
|
|
|
|
L Lд 1 Bt |
|
. |
|
|
(1.25) |
|||||
|
|
n 1 |
||||||||||
|
dl |
|
Lд |
|
4 |
|
1 Bt |
n 3 |
|
|||
|
|
|
n 1 |
|
|
|||||||
|
dt |
B n 1 |
(1.26) |
|||||||||
|
|
|
|
|
|
|
||||||
|
dL |
GдUд 1 |
|
n 3 |
||||||||
|
|
|
||||||||||
или |
Bt n 1 . |
|
||||||||||
|
dt |
|
|
|
|
|
|
|
|
|
|
|
Тогда, учитывая уравнения (1.2), (1.11), можно записать выражения реакции газов в дульном отверстии:
|
R P |
dL |
|
|
|
|
||
|
отд |
dt |
|
|
|
(1.27) |
||
|
|
|
|
|
|
|||
|
|
2n |
|
|
n 3 |
|||
или |
R SPд 1 Bt |
GдUд 1 |
Bt |
, |
||||
n 1 |
n 1 |
|
17
где
U |
|
|
k |
2 |
|
|
|
|
1 |
k 1M 2RT |
. |
(1.28) |
|||
д |
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
2 |
д |
д |
|
|||||
|
|
|
|
k 1 |
|
|
|
|
|
|
|
||||
При t 0 |
|
|
R Rд |
SPд GдUд. |
|
|
|
|
Изотермическое истечение
В соответствии с (1.17) или (1.18) и принимая n 1,
|
или L L e 2 |
Gд |
|||
L L e 2вt |
|
|
t . |
||
|
|
||||
д |
д |
|
|
|
|
Обозначая через a 2В, илиa 2 |
Gд |
|
, |
|
|
|
|
|
|||
|
|
|
|
|
изменения количества движения:
(1.29)
получим скорость
dL aL e at |
или |
dL G U |
e at . |
(1.30) |
|||||
dt |
д |
|
|
dt |
д |
д |
|
|
|
|
|
|
|
|
|
|
|||
Реакция газов в дульном отверстии |
|
|
|
||||||
|
R SPдe вt |
GдUдe at . |
|
|
(1.31) |
||||
Все полученные выше зависимости отличаются показа- |
|||||||||
тельным характером функции |
Pотд t , |
dL t |
, R t . |
На рис 1.4 |
|||||
dt |
|||||||||
|
|
|
|
|
|
|
|
показаны эти зависимости для периода последействия 130-мил- лиметрового орудия с представленными в параграфе 1.4 исходными баллистическими данными.
Из графика видно, что нестационарность истечения, характеризуемая скоростью изменения количества движения газов
в стволе dLdt , значительно (на 40–60 %) уменьшает силу отдачи
орудия Pотд в сравнении с реакцией газов в дульном отверстии Rд. В этом заключается преимущественное отличие про-
18
цесса отдачи артиллерийского ствола от силы тяги реактивного двигателя, в котором dLdt 0 и Pотд R.
Рис. 1.4. Зависимость силовых характеристик периода последействия от времени
1.5. Параметры отдачи оружия при выстреле
Процессы отдачи оружия при выстреле характеризуются текущей силой Pотд t и интегральной величиной, называемой
импульсом отдачи, Jотд t t Pотд t dt. При этом импульс отда-
0
чи Jотд есть величина, аккумулируемая стволом как механиче-
ской системой. При откате ствола этот импульс преобразуется в количество движения.
Зависимость силы отдачи от времени определяется решением задач внутренней и промежуточной баллистики:
Pотд t SP t или |
|
SPд 1 Bt |
2n |
|
||
Pотд |
|
; |
|
|||
n 1 |
(1.32) |
|||||
P |
SP e вt . |
|||||
|
||||||
отд |
д |
|
|
|
|
19
Импульс отдачи есть площадь под кривой Pотд t в каж-
дый текущий момент времени.
В период движения снаряда по каналу ствола импульс отдачи зависит от скоростей и масс снаряда и газопороховой смеси.
Учитывая уравнение движения снаряда
dV |
|
SP t |
|
V |
|
|
|
, |
имеем Jотд t qdV qV. |
(1.33) |
|||
dt |
q |
|||||
|
|
0 |
|
Этот же импульс можно получить и на основе закона сохранения количества движения масс внутри канала ствола с учетом гипотезы Пиобера,
Jотд q 0,5 V . |
(1.34) |
Показанные выше зависимости дают достаточно близкие результаты,однаконапрактикенаиболеечастоиспользуется(1.34).
Вмоментвылетаснарядаимпульсотдачивосновномпериоде
Jотд.л |
q 0,5 Vд. |
|
|
|
|
|
(1.35) |
|||||||||
Для периода последствия, учитывая (1.11) для политроп- |
||||||||||||||||
ного истечения, имеем |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t |
|
2n |
|
|
n 1SPд |
|
|
n 1 |
|
|
|
|||||
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
||||||||||
Jотд SPд 1 Bt |
|
n 1 dt |
|
|
|
|
|
1 1 |
Bt |
|
n 1 |
. |
(1.36) |
|||
|
b 1 B |
|
||||||||||||||
0 |
|
|
|
|
|
|
|
|
|
|
|
|
||||
При t Jотд |
|
– импульс отдачи периода последейст- |
||||||||||||||
вия, примет вид |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
J |
отд |
|
n 1SPд |
. |
|
|
|
|
|
(1.37) |
||||
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
n 1 B |
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Для изотропного истечения |
|
|
|
|
|
|
|
|
||||||||
Jотд SPд t e вt |
SPд |
|
1 e вt |
. |
|
|
|
|
(1.38) |
|||||||
|
|
|
|
|
||||||||||||
|
|
|
|
0 |
|
|
B |
|
|
|
|
|
|
|
|
20