Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Численное моделирование нестационарных переходных процессов в активных и реактивных двигателях

..pdf
Скачиваний:
4
Добавлен:
20.11.2023
Размер:
21.12 Mб
Скачать

112

 

Глава 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. Дополнительный заряд, имея скорость движения, близкую по величине к скорости спутного

Соседние файлы в папке книги