книги / Численное моделирование нестационарных переходных процессов в активных и реактивных двигателях
..pdf112 |
|
Глава 2. Моделирование артиллерийского выстрела |
||||||
t, млс |
2,671 |
2,681 |
2,692 |
2,702 |
2,713 |
2,723 |
2,733 |
2,743 |
1 |
25 |
26 |
27 |
28 |
29 |
30 |
31 |
32 |
Т,К |
1074 |
1046 |
1064 |
1049 |
1081 |
1075 |
1074 |
1079 |
t, млс |
2,752 |
2,762 |
2,772 |
2,781 |
2,790 |
2,800 |
2,809 |
2,818 |
1 |
33 |
34 |
35 |
36 |
37 |
38 |
39 |
40 |
Т,К |
1048 |
1062 |
1082 |
1063 |
1070 |
1081 |
1073 |
1069 |
t, млс |
2,827 |
2,836 |
2,846 |
2,858 |
2,878 |
2,885 |
2,893 |
2,900 |
1 |
41 |
42 |
43 |
44 |
45 |
46 |
47 |
48 |
Т,К |
1069 |
1067 |
1090 |
1074 |
1061 |
1049 |
1075 |
1070 |
t, млс |
2,907 |
2,913 |
2,920 |
2,923 |
2,926 |
2,932 |
2,938 |
2,944 |
1 |
49 |
50 |
51 |
64 |
52 |
53 |
54 |
55 |
Т,К |
1067 |
1046 |
1070 |
1065 |
1059 |
1054 |
1053 |
1057 |
t, млс |
2,950 |
2,955 |
2,960 |
2,962 |
2,965 |
2,969 |
2,970 |
2,975 |
1 |
56 |
57 |
58 |
62 |
59 |
60 |
61 |
63 |
Т,К |
1067 |
1044 |
1032 |
1033 |
1031 |
1036 |
1036 |
1019 |
Зажигание первой расчётной точки /=36 на поверхности нижней сгораемой гильзы таблица 2.3 (конструктивно нижняя гильза состоит из несгораемой и сгораемой частей, сгораемая часть примерно на 1/3 охватывает нижний полузаряд со стороны верхнего полузаряда (см. рис.2.1)) происходит в момент времени /=2,166шс при температуре поверхности Тк=966К. Все остальные точки сгораемой гильзы вплоть до 1=55 зажигаются последовательно, постепенно увеличивая температуру поверх ности при вспышке до Тк =1043ЛГ. Время полного воспламенения нижней сгораемой гильзы - /=2,390жлс. Временной интервал от первой вспышки на поверхности гильзы (/=36) до последней (1=55) составляет Af = 0,224млс. Заметим, что темп воспла менения нижней сгораемой гильзы несколько выше, а температура вспышки несколько ниже, чем соответствующие параметры нижнего полузаряда. Это обусловлено иным, чем для нижнего полузаряда, химическим составом пороха.
Таблица 2.3. Время, координата и температура вспышки по длине нижней сгораемой гильзы (Т0= 293К ).
t, млс |
2,166 |
2,179 |
2,191 |
2,203 |
2,215 |
2,227 |
2,239 |
2,251 |
1 |
36 |
37 |
38 |
39 |
40 |
41 |
42 |
43 |
Т,К |
966 |
980 |
976 |
976 |
976 |
979 |
985 |
993 |
t, млс |
2,263 |
2,274 |
2,286 |
2,297 |
2,308 |
2,319 |
2,331 |
2,342 |
1 п |
44 |
45 |
46 |
47 |
48 |
49 |
50 |
51 |
2.5. Результаты численного моделирования |
|
|
|
113 |
||||
т,к |
1004 |
992 |
1009 |
1002 |
995 |
992 |
1018 |
1020 |
1, МЛС |
2,352 |
2,364 |
2,374 |
2,390 |
|
|
|
|
1 |
52 |
53 |
54 |
55 |
|
|
|
|
Т,К |
994 |
1023 |
1007 |
1043 |
|
|
|
|
|
Зажигание поверхности верхней сгораемой гильзы таблица |
2.4 происходит аналогично зажиганию поверхности верхнего полузаряда. Отличие в основном заключается в номиналах температуры вспышки при примерно равном темпе зажигания. Зажигание первой расчётной точки 1=1 на поверхности верхней сгораемой гильзы происходит в момент времени 1=2,4Юл/лс при
температуре поверхности |
Г0 = 1024#. Зажигание последней точки |
||
/=63 - в момент времени |
1=2,960.шс при температуре поверхности |
||
Т0 =968К . |
Временной |
интервал от первой вспышки на |
|
поверхности сгораемой гильзы до последней вспышки |
составляет |
||
At = 0,55млс. |
Температура вспышки изменяется в |
пределах |
|
Г0 = 948#... 1024#. |
|
|
Таблица 2.4. Время, координата и температура вспышки по длине
верхней сгораемой гильзы (Т0 = 293К ).
t, млс |
2,410 |
2,419 |
2,431 |
2,442 |
2,454 |
2,465 |
2,476 |
2,487 |
|
1 |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
|
Т ,К |
1024 |
998 |
1010 |
996 |
1012 |
1002 |
995 |
990 |
|
t, млс |
2,498 |
2,509 |
2,520 |
2,531 |
2,542 |
2,553 |
2,563 |
2,574 |
|
1 |
9 |
10 |
11 |
12 |
13 |
14 |
15 |
16 |
|
Т ,К |
987 |
986 |
987 |
990 |
995 |
998 |
986 |
995 |
|
t, млс |
2,584 |
2,595 |
2,604 |
2,613 |
2,617 |
2,628 |
2,640 |
2,651 |
|
1 |
17 |
18 |
19 |
20 |
21 |
22 |
23 |
24 |
|
Т ,К |
984 |
995 |
978 |
990 |
990 |
971 |
980 |
977 |
|
t, млс |
2,662 |
2,673 |
2,684 |
2,694 |
2,704 |
2,714 |
2,724 |
2,734 |
|
1 |
25 |
26 |
27 |
28 |
29 |
30 |
31 |
32 |
|
Т ,К |
979 |
986 |
990 |
984 |
978 |
975 |
974 |
974 |
|
ty млс |
2,744 |
2,754 |
2,764 |
2,773 |
2,782 |
2,792 |
2,801 |
2,810 |
|
1 |
33 |
34 |
35 |
36 |
37 |
38 |
39 |
40 |
|
Т, К |
979 |
985 |
991 |
984 |
975 |
989 |
984 |
981 |
|
ty млс |
2,819 |
2,828 |
2,838 |
2,851 |
2,866 |
2,875 |
2,882 |
2,889 |
|
1 |
41 |
42 |
43 |
44 |
45 |
46 |
47 |
48 |
|
Т ,К |
979 |
977 |
983 |
997 |
995 |
1000 |
995 |
994 |
|
t, млс |
2,896 |
2,899 |
2,903 |
2,909 |
2,915 |
2,921 |
2,927 |
2,933 |
|
|
---- —--- |
||||||||
__ 1 |
49 |
64 |
50 |
51 |
52 |
53 |
54 |
55 |
J |
114 |
|
Глава 2. Моделирование артиллерийского выстрела |
||||||
т,к |
993 |
984 |
1000 |
992 |
985 |
980 |
975 |
976 |
|
2,939 |
2,944 |
2,950 |
2,952 |
2,954 |
2,958 |
2,960 |
2,960 |
1 |
56 |
57 |
58 |
62 |
59 |
60 |
61 |
63 |
Т,К |
976 |
964 |
969 |
949 |
951 |
948 |
952 |
968 |
В таблицах 2.5-2.8 представлены данные по времени и температуре зажигания для каждой из расчётных точек по длине нижнего полузаряда, верхнего полузаряда, нижней сгораемой гильзы и верхней сгораемой гильзы при начальной температуре Т0= 223К . Развитие процесса воспламенения для каждого из компонентов порохового заряда артиллерийского выстрела в целом аналогичено развитию процесса при Т0 = 293К . Отличие заключается в большем времени задержки зажигания поверхности горения и в большем номинале температуры поверхности при вспышке. Процесс зажигания как бы тормозится во времени. Нестационарные эффекты при прогреве поверхности горения усиливаются.
Таблица 2.5. Время, координата и температура вспышки по длине нижнего полузаряда (Г0 = 223К ).
t, МПС |
0,948 |
1,082 |
1,202 |
1,295 |
1,373 |
1,440 |
1,499 |
1,553 |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
Т,К |
1070 |
1011 |
1037 |
1072 |
1057 |
1053 |
1033 |
1096 |
t, млс |
1,601 |
1,645 |
1,659 |
1,686 |
1,724 |
1,760 |
1,793 |
1,824 |
1 |
10 |
11 |
1 |
12 |
13 |
14 |
15 |
16 |
Т,К |
1059 |
1038 |
1236 |
1046 |
1052 |
1102 |
1078 |
1070 |
t, млс |
1,853 |
1,881 |
1,907 |
1,932 |
1,955 |
1,978 |
1,999 |
2,020 |
1 |
17 |
18 |
19 |
20 |
21 |
22 |
23 |
24 |
Т,К |
1053 |
1085 |
1072 |
1084 |
1062 |
1078 |
1052 |
1090 |
t, млс |
2,039 |
2,058 |
2,076 |
2,094 |
2,111 |
2,127 |
2,143 |
2,158 |
1 |
25 |
26 |
27 |
28 |
29 |
30 |
31 |
32 |
Т,К |
1049 |
1062 |
1061 |
1107 |
1119 |
1087 |
1099 |
1090 |
t, млс |
2,173 |
2,188 |
2,201 |
2,214 |
2,227 |
2,240 |
2,253 |
2,265 |
1 |
33 |
34 |
35 |
36 |
37 |
38 |
39 |
40 |
Т,К |
1071 |
1092 |
1069 |
1059 |
1063 |
1080 |
1109 |
1075 |
t, млс |
2,278 |
2,290 |
2,302 |
2,314 |
2,326 |
2,338 |
2,350 |
2,362 |
1 |
41 |
42 |
43 |
44 |
45 |
46 |
47 |
48 |
Т,К |
1122 |
1098 |
1085 |
1079 |
1081 |
1088 |
1101 |
1120 |
t, млс |
2,374 |
2,385 |
2,397 |
2,409 |
2,420 |
2,441 |
|
|
1 |
49 |
50 |
51 |
52 |
53 |
54 |
|
|
2.5. Результаты численного моделирования |
|
|
115 |
|||||
Т . К |
1144 |
1088 |
1121 |
1152 |
1141 |
1137 |
|
|
Таблица 2.6. Время, координата и температура вспышки по длине |
||||||||
t, млс |
2,466 |
верхнего полузаряда (Г0 = 2 2 3 К ). |
|
|
||||
2,478 |
2,490 |
2,503 |
2,515 |
2,527 |
2,539 |
2,551 |
||
1 |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
Т ,К |
1079 |
1120 |
1089 |
1138 |
1114 |
1100 |
1093 |
1094 |
t, млс |
2,563 |
2,575 |
2,586 |
2,598 |
2,609 |
2,621 |
2,632 |
2,643 |
1 |
9 |
10 |
11 |
12 |
13 |
14 |
15 |
16 |
Т ,К |
1101 |
1117 |
1071 |
1097 |
1068 |
1106 |
1086 |
1073 |
t, млс |
2,654 |
2,665 |
2,676 |
2,685 |
2,691 |
2,703 |
2,715 |
2,727 |
1 |
17 |
18 |
19 |
20 |
21 |
22 |
23 |
24 |
Т ,К |
1067 |
1086 |
1092 |
1065 |
1101 |
1083 |
1079 |
1095 |
t, млс |
2,738 |
2,750 |
2,761 |
2,772 |
2,783 |
2,793 |
2,804 |
2,815 |
1 |
25 |
26 |
27 |
28 |
29 |
30 |
31 |
32 |
Т ,К |
1066 |
1107 |
1097 |
1098 |
1111 |
1070 |
1100 |
1079 |
t, млс |
2,824 |
2,834 |
2,844 |
2,854 |
2,864 |
2,874 |
2,883 |
2,893 |
1 |
33 |
34 |
35 |
36 |
37 |
38 |
39 |
40 |
Т ,К |
1068 |
1065 |
1068 |
1077 |
1092 |
1115 |
1079 |
1111 |
t, млс |
2,902 |
2,912 |
2,921 |
2,934 |
2,955 |
2,962 |
2,970 |
2,977 |
1 |
41 |
42 |
43 |
44 |
45 |
46 |
47 |
48 |
Т ,К |
1085 |
1123 |
1071 |
1102 |
ИЗО |
1078 |
1097 |
1070 |
t, млс |
2,985 |
2,992 |
2,998 |
3,005 |
3,005 |
3,012 |
3,018 |
3,025 |
1 |
49 |
50 |
51 |
52 |
64 |
53 |
54 |
55 |
Т ,К |
1127 |
1112 |
1064 |
1079 |
1091 |
1102 |
1070 |
1099 |
t, млс |
3,031 |
3,037 |
3,043 |
3,049 |
3,049 |
3,054 |
3,056 |
3,068 |
1 |
56 |
57 |
58 |
59 |
62 |
60 |
61 |
63 |
Т ,К |
1074 |
1058 |
1050 |
1052 |
1082 |
1057 |
1053 |
1041 |
Таблица 2.7. Время, координата и температура вспышки по длине |
||||||||
t, млс |
|
нижней сгораемой гильзы (Т0 = 223К ). |
|
|
||||
2,212 |
2,224 |
2,236 |
2,249 |
2,261 |
2,273 |
2,285 |
2,297 |
|
1 |
36 |
37 |
38 |
39 |
40 |
41 |
42 |
43 |
Т ,К |
991 |
982 |
977 |
998 |
994 |
995 |
994 |
996 |
t, млс |
2,309 |
2,321 |
2,333 |
2,344 |
2,356 |
2,368 |
2,379 |
2,391 |
1 |
44 |
45 |
46 |
47 |
48 |
49 |
50 |
51 |
Т ,К |
1001 |
1008 |
1018 |
999 |
1010 |
1023 |
1004 |
1021 |
t,Mnc |
2,402 |
2,414 |
2,425 |
2,443 |
|
|
|
|
1 |
52 |
53 |
54 |
55 |
|
|
|
|
т ,к |
1010 |
1023 |
1031 |
1028 |
|
|
|
|
116 |
|
Гпава 2. Моделирование артиллерийского выстрела |
||||||
Таблица 2.8. Время, координата и температура вспышки по длине |
||||||||
|
|
верхней сгораемой гильзы (Т0 = 223К ). |
|
|
||||
t, млс |
2,466 |
2,476 |
2,488 |
2,500 |
2,513 |
2,525 |
2,537 |
2,549 |
1 |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
Т,К |
1079 |
1017 |
1004 |
995 |
1013 |
1007 |
1004 |
1004 |
t, млс |
2,561 |
2,573 |
2,584 |
2,596 |
2,607 |
2,619 |
2,630 |
2,641 |
1 |
9 |
10 |
11 |
12 |
13 |
14 |
15 |
16 |
Т,К |
1007 |
1008 |
994 |
1005 |
992 |
1004 |
1000 |
994 |
t, млс |
2,652 |
2,663 |
2,674 |
2,683 |
2,688 |
2,701 |
2,713 |
2,725 |
1 |
17 |
18 |
19 |
20 |
21 |
22 |
23 |
24 |
Т,К |
991 |
991 |
1001 |
989 |
979 |
995 |
994 |
997 |
t, млс |
2,736 |
2,747 |
2,759 |
2,770 |
2,780 |
2,791 |
2,802 |
2,812 |
1 |
25 |
26 |
27 |
28 |
29 |
30 |
31 |
32 |
Т,К |
980 |
981 |
998 |
998 |
982 |
990 |
999 |
991 |
t, млс |
2,822 |
2,832 |
2,842 |
2,852 |
2,862 |
2,872 |
2,881 |
2,891 |
1 |
33 |
34 |
35 |
36 |
37 |
38 |
39 |
40 |
Т,К |
989 |
987 |
989 |
993 |
999 |
1005 |
994 |
1004 |
Uмлс |
2,900 |
2,909 |
2,919 |
2,932 |
2,949 |
2,958 |
2,965 |
2,973 |
1 \ |
41 |
42 |
43 |
44 |
45 |
46 |
47 |
48 |
т,к |
996 |
987 |
990 |
1006 |
1024 |
1011 |
993 |
1006 |
t, млс |
2,980 |
2,987 |
2,991 |
2,994 |
3,001 |
3,007 |
3,014 |
3,030 |
1 |
49 |
50 |
64 |
51 |
52 |
53 |
54 |
55 |
Т,К |
995 |
995 |
987 |
998 |
1000 |
986 |
994 |
980 |
t, млс |
3,027 |
3,033 |
3,039 |
3,045 |
3,045 |
3,045 |
3,052 |
3,060 |
1 |
56 |
57 |
58 |
59 |
62 |
60 |
61 |
63 |
Т,К |
986 |
978 |
970 |
966 |
954 |
962 |
957 |
941 |
2. На рис.2.7-2.8 представлено изменение температуры в градусах Кельвина по своду порохового элемента (координата по своду - номер расчётной точки) в среднем сечении нижнего и верхнего полузаряда артиллерийского выстрела в различные моменты времени при начальной температуре Т0= 293К . Для других элементов порохового заряда изменения рассматриваемых параметров качественно аналогичны.
В момент времени /=2лхлс пороховой элемент в среднем сечении нижнего полузаряда (рис.2.7) прогрет на поверхности до температуры ~ Тк= 680АГ, но ещё не горит. В последующие моменты времени пороховой элемент в данном сечении находится в режиме горения. Максимальная температура по своду порохово-
2.5. Результаты численного моделирпяпм.а
117
Рис.2.7. Изменение температуры по своду порохового элемента нижнего полузаряда в различные моменты времени ( Т0 = 293К ).
Рис.2.8. Изменение температуры по своду порохового элемента верхнего
полузаряда в различные моменты времени (7^ —293К ).
118 |
Глава 2. Моделирование артиллерийского выстрела |
Рис.2.9. Изменение температуры по своду порохового элемента нижнего полузаряда в различные моменты времени (Г0 = 223К).
т ,к
— Змлс -*-4млс —А—5млс —к—бмлс ж 7млс ——8млс
Рис.2.10. Изменение температуры по своду порохового элемента верхнего полузаряда в различные моменты времени (TQ= 223К ).
2.5. Результаты численного моделирования |
119 |
го элемента - на поверхности горения. По мере роста давления в камере сгорания орудия растёт и температурное распределение по своду пороха. При t-Ьшс наблюдается пик давления в камере сгорания. Ему соответствует температура на поверхности горения
~Тк= 1680А:.
Вмомент времени t-Ъжс пороховой элемент в среднем сечении верхнего полузаряда (рис.2.8) прогрет на поверхности до
температуры |
~ Тк= 950К , |
а при t=6.же |
до температуры ~ |
Тк= 1480АГ. |
Максимальная |
температура по |
своду здесь также |
на поверхности горения. Изменение температуры по своду порохового элемента в силу особенностей воспламенительного периода (быстрое воспламенение) и меньшего перепада давления в данном сечении камеры сгорания по времени не столь значительно, чем для нижнего полузаряда.
На рис.2.9-2.10 представлено изменение температуры в градусах Кельвина по своду порохового элемента в среднем сечении нижнего полузаряда и верхнего полузаряда артил лерийского выстрела в различные моменты времени при
начальной температуре |
Т0= 223К . В целом характер изменения |
||||||
температуры |
по |
своду |
и |
по |
времени |
аналогичен |
случаю при |
Т0 = 293К . |
Но |
на тот |
же |
момент |
времени |
срабатывания |
артиллерийского выстрела номиналы температуры меньше. Так для нижнего полузаряда (рис.2.9) при t-2me температура на поверхности горения ~ Тк= 600 К (порох ещё не горит), а при
t=6mc |
~ Тк= 1550АГ. Для верхнего полузаряда (рис.2.10) при |
t-Зт с |
температура на поверхности горения ~ Тк= 940АГ, а при |
t=6mc ~ Тк= 1390АТ.
3. На рис.2.11 и рис.2Л2 показано изменение скорости горения (отношение текущей скорости горения к её стационарной величине £ = v*/v°) пороховых элементов заряда артил лерийского выстрела по времени (время дано в шс) при начальной температуре Т0 =293К . Изменение зафиксировано в среднем сечении нижнего и верхнего полузарядов (рис.2.11) и в среднем сечении верхней сгораемой гильзы и дополнительного заряда (рис.2.12). Данные графические зависимости хорошо иллюстри руют основные особенности изменения скорости горения заряда.
120 |
Глава 2. Моделирование артиллерийского выстрела |
Рис.2.11. Изменение по времени скорости горения нижнего полузаряда (крест) и верхнего полузаряда (звезда) при Т0= 293К .
Рис.2.12. Изменение по времени скорости горения верхней сгораемой гильзы (крест) и дополнительного заряда (звезда) при Г0 = 293К .
2.5. Результаты численного моделирования |
121 |
В первые моменты времени после зажигания рас сматриваемого сечения нижнего полузаряда (рис.2.11) (метка на графике - крест) скорость его горения существенно нестацио
нарная v* / v“ » 1 . При вспышке максимальное |
значение |
|||
отношения |
vk /v°k ~5. |
С течением |
времени нестационарность |
|
постепенно |
вырождается, и заряд горит в стационарном режиме |
|||
vk/ = 1 |
до полного |
сгорания. |
После зажигания |
рассмат |
риваемого сечения верхнего полузаряда (метка на графике - звезда) скорость его горения ещё более нестационарная. При вспышке максимальное значение отношения vk/v°k~ l. В
дальнейшем по времени нестационарность также вырождается, но полузаряд не горит в стационарном режиме - vk/v° > 1.
Вследствие обдува высокоскоростным спутным потоком продуктов сгорания поверхности горения скорость горения в рассматриваемом сечении верхнего полузаряда больше стационарной - полузаряд горит в турбулентном режиме. По мере уменьшения скорости обдува турбулентность вырождается и скорость горения верхнего полузаряда приближается к стационарной.
Влияние турбулентности на режим горения пороховых элементов заряда хорошо иллюстрируется зависимостью (рис.2.12). Верхняя сгораемая гильза (метка на графике - крест) и дополнительный заряд (метка на графике - звезда) изготовлены из одной марки пороха, но находятся по времени работы двигателя в различных условиях горения. Верхняя сгораемая гильза неподвижна, а дополнительный заряд смещается со снарядом от своего первоначального положения. После зажигания оба пороховых элемента в рассматриваемом сечении горят в нестационарном режиме - v* /v° » 1 . При вспышке максимальное
значение отношения - vk/vk = 6 . Со временем нестационарность вырождается. Поверхность горения верхней сгораемой гильзы интенсивно омывается спутным потоком продуктов сгорания и сгораемая гильза в рассматриваемом сечении горит в турбулентном режиме - vk/ v° > 1. Дополнительный заряд, имея скорость движения, близкую по величине к скорости спутного