<p>Обработка исходных данных и зависимости для расчёта свойств газов</p>

0f6c11f3-bb6c-431b-a223-61323bb25515 7 4 5 false false 0 decimal

&[DATE] &[TIME] - &[FILENAME]

Расчет пропускной способности трубопровода с дросселем на выходе

Общие параметры

Параметры расчёта

ГАЗ

Абсолютное давление газа во входном сечении трубопровода

p.вх 16 :

gas

Диаметр дросселя

Температура газа во входном сечении трубопровода

d.др 5 :

T.тр 293 :

ЕДИНИЦЫ ДАВЛЕНИЯ

Абсолютное давление газа на выходе дросселя

Коэффициент расхода дросселя

un_p

p.вых 0.1 :

μ.тр 0.9 :

Диаметр сечения трубопровода

ЕДИНИЦЫ ТЕМПЕРАТУРЫ

d.тр 12 :

un_T

Гидравлическое сопротивление трубопровода, 10^8 м^(-4)

S.тр 187 :

ЕДИНИЦЫ РАДИУСОВ

Шероховатость стенок трубы, мм

Ed_D

Δ.тр 0.02 :

<p>Обработка исходных данных и зависимости для расчёта свойств газов</p>

1. Перевод заданных значений

p.вх p.вх un_p * : 1.6 107^* Pa *

p.вых p.вых un_p * : 1105^* Pa *

d.тр d.тр Ed_D * : 0.012 m *

d.др d.др Ed_D * : 0.005 m *

T.тр T.тр 273.15 un_T * + : 293

S.тр S.тр 108^* m 4 -^* : 1.87 1010^* 1 m 4^/ *

Δ.тр Δ.тр mm * : 2105 -^* m *

2. Определение параметров газа

R.gas gas 1 ≡ 287.1 J kg K * / * gas 2 ≡ 296.8 J kg K * / * gas 3 ≡ 2077.3 J kg K * / * gas 4 ≡ 259.835 J kg K * / * gas 5 ≡ 518.271 J kg K * / * gas 6 ≡ 4124.9 J kg K * / * 63.32 J kg K * / * if if if if if if 11 line :

k.gas gas 1 ≡ 1.4 gas 2 ≡ 1.4 gas 3 ≡ 1.66 gas 4 ≡ 1.4 gas 5 ≡ 1.31 gas 6 ≡ 1.41 1.7 if if if if if if 11 line :

3. Выбор матриц данных в соответствии с названием газа

ρ gas 1 ≡ ρ_vozduh_tab gas 2 ≡ ρ_azot_tab gas 3 ≡ ρ_geliy_tab gas 4 ≡ ρ_kislorod_tab gas 5 ≡ ρ_metan_tab gas 6 ≡ ρ_vod_norm_tab ρ_ksenon_tab if if if if if if 11 line :

4. Задание функции, реализующей двумерную кусочно-линейную интерполяцию по табличным данным о свойстве газа.

P T M Function_prop P P : T T : p M 2 M rows 11 submatrix : t M 112 M cols submatrix transpose : Func M 2 M rows 2 M cols submatrix : i 1 : i Func cols ≤ i i 1 + : vec\002A\ i el p Func i col P linterp : for d t vec\002A\ T linterp : d 81 line :

p T ρ.интерп p MPa / T ρ Function_prop kg m 3^/ * :

5. Определение абсолютного давления в выходном сечении требопровода

5.1 Расчёт коэффициента Дарси

λ.тр 0.11 Δ.тр d.тр / (0.25^* : 0.0222

5.2 Расчёт площади сечения трубы

f.тр π d.тр 2^* 4 / : 0.0001 m 2^*

5.3 Расчёт эквивалентной длинны трубопровода

L.экв S.тр f.тр 2^* d.тр λ.тр / * : 129.1436 m *

5.4 Расчёт критического отношения давлений на дросселе

σ.кр 2 k.gas 1 + / (k.gas k.gas 1 - /^: 0.5283

5.5 Расчёт абсолютного давления в выходном сечении трубы

p.др_вх f.тр p T.тр ρ.интерп p p.вх p.вых int p.др_вх T.тр ρ.интерп p.вх T.тр ρ.интерп / ln λ.тр 2 / L.экв d.тр / * - / sqrt * μ.др π d.др 2^* 4 / * 2 k.gas * k.gas 1 - / 1 R.gas / * sqrt * p.др_вх T.тр sqrt / * p.вых p.др_вх / σ.кр > p.вых p.др_вх / (2 k.gas /^p.вых p.др_вх / (k.gas 1 + k.gas /^- sqrt σ.кр 2 k.gas /^σ.кр k.gas 1 + k.gas /^- sqrt if 11 line * - p.др_вх p.вых p.вх solve err try : err

5.6 Расчёт массового расхода через трубу

G.тр f.тр p T.тр ρ.интерп p p.вх p.вых int p.др_вх T.тр ρ.интерп p.вх T.тр ρ.интерп / ln λ.тр 2 / L.экв d.тр / * - / sqrt * : #

5.7 Расчёт массового расхода через дроссель

G.др μ.др π d.др 2^* 4 / * 2 k.gas * k.gas 1 - / 1 R.gas / * sqrt * p.др_вх T.тр sqrt / * p.вых p.др_вх / σ.кр > p.вых p.др_вх / (2 k.gas /^p.вых p.др_вх / (k.gas 1 + k.gas /^- sqrt σ.кр 2 k.gas /^σ.кр k.gas 1 + k.gas /^- sqrt if 11 line * : #

5.7 Критерий адекватности расчёта - относительная разность расчётных значений расхода через трубу и через дроссель

Δ.расх 100 G.тр G.Др - G.тр / * : #

6 Формирование значений параметров для выдачи в таблицу

6.1 Плотность перед дросселем

ρ.др_вх p.др_вх T.тр ρ.интерп err try : err

6.2 Абсолютное давление перед дросселем

p.др_вх p.др_вх err try : err

6.4 Пропускная способность трубопровода

G.тр G.тр err try : kg s / err

6.5 Скорость потока газа перед дросселем

V.вых G.тр π d.тр 2^* 4 / ρ.др_вх * / err try : err

Обработка вывода данных

p.др_вх p.др_вх MPa / 3 round 10000 - try : 10000 -

G.тр G.тр g s / / 3 round 10000 - try : 10000 -

V.вых V.вых m s / / 3 round 10000 - try : 10000 -

Δ.расх Δ.расх 1 round 10000 - try : 10000 -

ρ.др_вх ρ.др_вх kg m 3^/ / 3 round 10000 - try : 10000 -