e31a8284-d8c5-442d-bb85-a4b05e1f753a 253 6 5 auto

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

Дополнительно см. РМ 4-163-77 Расчет и применение регулирующих органов в системах автоматизации технологических процессов

T 0 °C * 10 °C * 20 °C * 30 °C * 40 °C * 50 °C * 60 °C * 70 °C * 80 °C * 90 °C * 100 °C * 110 °C * 120 °C * 130 °C * 140 °C * 150 °C * 160 °C * 170 °C * 118 mat transpose :

Абсолютное давление насыщения водяного пара, [Па] p.sw 611.2 1228.2 2339.2 4246.7 7384.4 1235119946312014741570182101420143380198670270260361500476100618140792050118 mat transpose :

t p.sw T p.sw t K / ainterp eval Pa * :

t p.sw_изб t p.sw 1 atm * - :

145 °C * p.sw_изб bar 3.1432

48 °C * p.sw bar 0.1121

Удельный объём воды при p = 1 МПа, [×10⁻³ м³/кг] v.w 0.9997 0.9999 1.0014 1.0040 1.0074 1.0117 1.0167 1.0223 1.0286 1.0355 1.0430 1.0511 1.0599 1.0693 1.0794 1.0902 1.1017 1.1141 118 mat transpose :

Удельный объём воды при p = 1 МПа t v.w T v.w t K / linterp 103 -^* eval m 3^kg / * :

Плотность воды при p = 1 МПа t ρ.w 1 t v.w / :

ρ.w_max 1000 kg m 3^/ * :

Коэффициент температурного расширения t β.w ρ.w_max t ρ.w / :

Динамическая вязкость воды при p = 1 МПа, [×10^-6 Па с] μ.w 179013061002 797.5 653.2 547.1 466.7 404.2 354.6 314.7 282 254.9 232.3 213.1 196.7 182.6 170.3 159.6 118 mat transpose :

t μ.w T μ.w t K / linterp eval 106 -^* Pa * s * :

Кинематическая вязкость, [м²/с] μ ρ η μ ρ / :

Объёмный расход через
массовый расход, [м³/c] G ρ q G ρ / :

Площадь сечения
круглой трубы, [м²] d F π 4 / d 2^* :

Скорость потока расхода q
через трубу сечением F, [м/с] q F v.F q F / :

Скорость потока расхода q через
круглую трубу диаметром d, [м/с] q d v q d F v.F :

ммHg 101325760 / Pa * :

мH2O 0.1 kgf cm 2^/ * :

Число Рейнольдса v d η Re v d * η / :

Коэффициент гидравлического трения (по Альтшулю) q d Δ.э η λ 0.11 Δ.э d / 68 q d v d η Re / + (0.25^* :

Кинематическая вязкость
воды, [м²/с] t η.w t μ.w t ρ.w η :

usgal 231 in 3^* :

Динамическое давление, [Па] ρ v p.d ρ v 2^* 2 / :

Удельное динамическое давление, [Па/(кг/с)²]
(динамическое давление ед. расхода) ρ d A ρ 1 ρ / d v p.d :

Динамическое давление через удельное динамическое давление A G p.dA A G 2^* :

Потери давления на участке

Формула Вейсбаха ζ p.d Δp ζ p.d * :

Пересчёт расхода на другие параметры q.1 Δp.1 ρ.1 Δp.2 ρ.2 q.2 q.1 Δp.2 Δp.1 / ρ.1 ρ.2 / * sqrt * :

Формула Вейсбаха ζ q d ρ Δp ζ ρ q d v p.d Δp :

ζ d Δp ρ q d F 2 Δp * ζ ρ * / sqrt * :

Потеря напора через потерю давления Δp ρ Δh.ρ Δp ρ g.e * / :

Потеря давления через потерю напора Δh ρ Δp.ρ ρ g.e * Δh * :

Пересчёт потерь на другие параметры q.1 Δp.1 ρ.1 q.2 ρ.2 Δp.2 Δp.1 q.2 q.1 / (2^* ρ.2 ρ.1 / * :

Через потерю напора Δh v ζ.Δh Δh 2 g.e * v 2^/ * :

Через потерю давления Δp ρ v ζ.Δp Δp ρ v p.d / :

Участка трубы (Формула Дарси) l d λ ζ.l λ l d / * :

Участка трубы (Формула Дарси) q d l Δ.э η ζ.l l d q d Δ.э η λ ζ.l :

ζ d λ l.э ζ d λ / * :

Формула Дарси-Вейсбаха q d l Δ.э t Δp q d l Δ.э t η.w ζ.l q d t ρ.w Δp :

Характеристика сопротивления
(потери давления на  ед. расхода)
[Па/(кг/с)²] ζ ρ d S ρ d A ζ * :

Через характеристику сопротивления S S G Δp.S S G 2^* :

Через характеристику сопротивления S S ρ d ζ.S S ρ d A / :

Коэффициент расхода
ε - коэффициент сжатия струи
a - коэффициент Кориолиса
http://gidravl.narod.ru/istechenie.html ζ ε a μ ε a ζ + sqrt / :

http://www.nex7.narod.ru/gidra/8.htm μ F Δp ρ q.μ μ F * 2 Δp * ρ / sqrt * :

μ F ΔH q.μ μ F * 2 g.e * ΔH * sqrt * :

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

Av [m²]

ζ F Av F 2 ζ / sqrt * :

Av F ζ.Av 2 F Av / (2^* :

Av Δp ρ q.Av Av Δp ρ / sqrt * :

Пропускная способность Kv (Flow factor) (по ГОСТ Р 52720-2007) [m³/hr]

ГОСТ Р 52720-2007 ρ.Kv 1000 kg m 3^/ * :

ГОСТ Р 52720-2007 Δp.Kv 0.1 MPa * :

Kv (flow factor) из ζ ζ d Kv ζ d F Av Δp.Kv ρ.Kv q.Av :

Kv (flow factor) из Δp q Δp ρ Kv q Δp ρ Δp.Kv ρ.Kv q.2 :

Коэффициент местного сопротивления из Kvs Kv d ζ Δp.Kv ρ.Kv Kv d v ζ.Δp :

Kv Δp ρ q Kv Δp.Kv ρ.Kv Δp ρ q.2 :

q Kv ρ Δp Kv Δp.Kv ρ.Kv q ρ Δp.2 :

Cv (Flow coefficient) [usgal/min]

Cv->Kv (безразмерные) Cv Kv Cv usgal min / psi sqrt / * m 3^hr / bar sqrt / (1 -^* :

Kv->Cv (безразмерные) Kv Cv Kv m 3^hr / bar sqrt / * usgal min / psi sqrt / (1 -^* :

Cv->Kv (gal/min->m³/hr) Cv.u Kv.u Cv.u bar psi / sqrt * :

Kv->Cv (m³/hr->gal/min) Kv.u Cv.u Kv.u psi bar / sqrt * :

Kv регулятора

Регулирующая функция 
h   - положение регулятора (0..1)
ep - вектор, первый элемент которого -
            h начала равнопроцентной хар-ки,
            второй - его пропускная способность h ep KvTrim h 0 < (h 1 > (| ∞ h ep 1 el ≤ h ep LFunc h ep EPFunc if if :

h ep LFunc h ep 1 el / ep 2 el * :

h ep EPFunc ep 2 el (1 h - 1 ep 1 el - /^:

h KvEPHelper 0.0183 4 h * exp * :

h KvEPHelper 150 / (1 h -^:

EPline 1121 mat :

Kv при относительном положении регулятора h Kvs h ep Kv.h Kvs h ep KvTrim * :

Последовательное соединение гидравлических сопротивлений Kv.vec Kv.посл Kv.vec k l el (2 -^(l 1 Kv.vec cols sum k 1 Kv.vec rows sum 2 - nthroot :

Последовательное соединение гидравлических сопротивлений Kv.vec Kv.посл Kv.vec 2 -^vectorize eval sum 2 - nthroot :

Предельно допустимый перепад давлений на клапане (по кавитации).
Z - коэффициент начала кавитации (характеристика клапана по VDMA 24 422)
P1 - избыточное давление теплоносителя перед клапаном
t - температура теплоносителя Z P.1 t Δp.кл_max Z P.1 t p.sw_изб - (* :

ζ.задв 0.5 :

ζ.вент_кос 0.5 :

ζ.вент_верт 6.0 :

ζ.кл_обр 7.0 :

ζ.кл_обр_захл 3.0 :

ζ.кран_пр 2 :

ζ.комп_сальн 0.3 :

ζ.комп_П 2.8 :

Значения коэффициентов сопротивления взяты по данным справочной системы ZuluThermo См. также Идельчик И.Е. Справочник по гидравличесим сопротивлениям

ζ.отв_90_скл_3d 0.8 :

ζ.отв_90_скл_4d 0.5 :

ζ.отв_90_гл_1d 1.0 :

ζ.отв_90_гл_3d 0.5 :

ζ.отв_90_гл_4d 0.3 :

ζ.отв_30_1ш 0.2 :

ζ.отв_45_1ш 0.3 :

ζ.отв_60_1ш 0.7 :

ζ.отв_90_2ш 0.6 :

ζ.отв_90_3ш 0.5 :

ζ.тр_сл_пр 1.2 :

ζ.тр_сл_отв 1.8 :

ζ.тр_разв_пр 1.0 :

ζ.тр_разв_отв 1.5 :

ζ.тр_встр 3.0 :

ζ.расш 1.0 :

ζ.суж 0.5 :

ζ.гряз 10 :

ζ.кр_ш_нпп 1.8 :

ζ.кр_ш_стд 0.8 :

ζ.кр_ш_пп 0.5 :

1. Характеристика регулируемого участка Здесь необходимо определить общую пропускную способность участка с учётом местных и линейных сопротивлений. Результат нужно сохранить в ΣKv

Располагаемый напор на врезке ΔP 43 мH2O * :

Температура подачи, расход
и температура обратки; число ТО Режим 115 °C * 1688 m 3^hr / * 70 °C * 4115 °C * 1088 m 3^hr / * 70 °C * 469 °C * 2367 m 3^hr / * 48 °C * 469 °C * 1083 m 3^hr / * 48 °C * 469 °C * 1284 m 3^hr / * 48 °C * 454 mat :

t Режим 1 col :

Диаметр трубы d 250 mm * :

Длина участка l 15 m * :

ΔP.доЦТП # :

Q Режим 2 col :

D 700 mm * :

t.обр Режим 3 col :

Эквивалентная шероховатость Δ.э 0.5 mm * :

L 470 m * 2 1.6 * (* :

n.ТО Режим 4 col :

ΔP.доЦТП k 1 n range M0 k el Q k el D L Δ.э t k el Δp : for M0 21 line :

ΔP.доЦТП Q D L Δ.э t Δp vectorize :

Сетчатый фильтр Kvs.ф 1600 m 3^hr / * :

n Режим rows :

Q 0.468889 m 3^s / * 0.302222 m 3^s / * 0.6575 m 3^s / * 0.300833 m 3^s / * 0.356667 m 3^s / * 51 mat

q Q n.ТО / vectorize eval :

ΔP.ЦТП ΔP ΔP.доЦТП - eval :

ΔP.ЦТП мH2O 40.2098 41.8364 37.3258 41.8009 41.319 51 mat

Теплообменник Kvs.то 592 m 3^hr / * 4.31 мH2O * 69 °C * ρ.w Kv :

q m 3^hr / #

n.ТО 4444451 mat

Отводы и арматура Σζ.м ζ.отв_90_гл_3d 6 * ζ.отв_90_гл_3d 2 / 2 * + 0.5 + ζ.кр_ш_пп 3 * + :

Искомый авторитет клапана A 0.5 :

271 m 3^hr / * Kvs.то 69 °C * ρ.w Δp мH2O 0.9032

Трубы q t Δp.тр q d l Δ.э t Δp :

Отводы и арматура q t Δp.м Σζ.м q d t ρ.w Δp :

Сетчатый фильтр q t Δp.ф q Kvs.ф t ρ.w Δp :

Теплообменник q t Δp.то q Kvs.то t ρ.w Δp :

Весь участок без регулирующего клапана q t ΣΔp q t Δp.тр q t Δp.м + q t Δp.ф + q t Δp.то + :

Коэффициент сопротивления труб ζ.тр q k el d l Δ.э t k el η.w ζ.l k 1 n sum n / eval :

Коэффициент сопротивления фильтра ζ.ф Kvs.ф d ζ eval :

Коэффициент сопротивления теплообменника ζ.то Kvs.то d ζ eval :

ζ.ф 2.4397

Σζ ζ.тр ζ.ф + ζ.то + Σζ.м + eval :

ΣKv ζ.тр d Kv Σζ.м d Kv Kvs.то Kvs.ф 14 mat Kv.посл eval :

ΣKv m 3^hr / #

2. Подбор Kvs и располагаемого напора

k.max kk 1 : k 2 n range q k el t k el ΣΔp q kk el t kk el ΣΔp > kk k : kk kk : if for kk 31 line eval :

q.max q k.max el :

t.max t k.max el :

Перепад напора на регулируемом сегменте
(кроме клапана) в максимальном режиме ΣΔp.max q.max t.max ΣΔp eval :

Kvs Q.helper Σζ Kvs d ζ + d q.max Kvs 1.2 / t.max Δp ΣΔp.max + t.max q :

ΣΔp.max мH2O #

Δp.кл_рек ΣΔp.max A 1 A - / * :

Kvs.кл_рек Kvs m 3^hr / * Q.helper Kvs m 3^hr / * t.max Δp Σζ Kvs m 3^hr / * Q.helper d t.max Δp / A 1 A - / ≡ Kvs 01000 solve m 3^hr / * :

Δp.кл_рек мH2O #

Kvs.кл_рек q.max Δp.кл_рек t.max ρ.w Kv :

Kvs.кл_рек m 3^hr / #

Рекомендованный общий перепад в регулируемом сегменте Δp.c_рек q.max Kvs.кл_рек 1.2 / t.max ρ.w Δp ΣΔp.max + :

Δp.c_рек мH2O #

q.max Kvs.кл_рек 1.2 / t.max ρ.w Δp ΣΔp.max / #

q.max 400 m 3^hr / * t.max ρ.w Δp мH2O #

Δp.бал_рек ΔP.ЦТП k.max el Δp.c_рек - :

Δp.бал_рек мH2O #

Принятый перепад на сегменте Δp.c 0.16 MPa * :

Перепад на регуляторе перепада Δp.бал ΔP.ЦТП Δp.c - :

Коэффициент запаса K.з 1.2 :

Δp.бал мH2O 23.8943 25.5209 21.0103 25.4855 25.0036 51 mat

Kv.кл k 1 n range M 1 k el q k el Δp.c q k el t k el ΣΔp - t k el ρ.w Kv : M 2 k el Q k el Δp.бал k el t k el ρ.w Kv : 21 line for M 21 line :

Kv.кл m 3^hr / #

Регулирующий клапан Kvs.1 Kv.кл 1 row max K.з * :

Регулятор перепада давлений Kvs.2 Kv.кл 2 row max K.з * :

Kvs.1 m 3^hr / #

Kvs.2 m 3^hr / 1958

dd.2 50 mm * :

Принятый для клапана Kvs 1000 m 3^hr / * :

Точка EP EP 0.4 0.15 21 mat :

Расчётная температура tt 115 °C * :

Расход через всю систему при заданном Kv клапана Kv t Q ΣKv Kv 12 mat Kv.посл Δp.c t ρ.w q :

Определить положение регулятора для достижения желаемого расхода Q1 ep solveQ Lh Kvs h ep LFunc * tt Q Q1 ≡ h 01 ep 2 el / ep 1 el * solve : Lh ep 1 el > Kvs h ep EPFunc * tt Q Q1 ≡ h ep 1 el 1 solve Lh if 21 line :

При открытом клапане Q.max Kvs tt Q eval :

Q.max m 3^hr / #

Относительный расход системы h ep q' Kvs h ep Kv.h tt Q Q.max / :

h ep dp.кл Kvs h ep Kv.h tt Q Kvs h ep Kv.h tt ρ.w Δp :

Относительное падение напора на клапане h ep dp'.кл h ep dp.кл Δp.c / :

В остальной системе h ep dp'.с 1 h ep dp'.кл - :

k 1 n range h.ep k el q k el EP solveQ : h.l k el q k el EPline solveQ : 21 line for

h.ep1 h.ep min :

h.l1 h.l min :

h.ep2 h.ep max :

h.l2 h.l max :

Равнопроцентная характеристика Ось x - положение регулятора зелёная кривая - падение напора на клапане синяя кривая - собственная хар-ка клапана красная кривая - действит. хар-ка клапана

x EP KvTrim x EP q' x EP dp'.кл h.ep1 h.ep1 EP q' o 10 red h.ep1 h.ep1 EP dp'.кл o 10 green 2 5 mat h.ep2 h.ep2 EP q' x 10 red h.ep2 h.ep2 EP dp'.кл x 10 green 2 5 mat 0 h.ep1 EP dp'.кл h.ep1 h.ep1 EP dp'.кл 2 2 mat 0 h.ep1 EP q' h.ep1 h.ep1 EP q' 2 2 mat 0 h.ep2 EP dp'.кл h.ep2 h.ep2 EP dp'.кл 2 2 mat 0 h.ep2 EP q' h.ep2 h.ep2 EP q' 2 2 mat h.ep 1 el 0 h.ep 1 el 1 2 2 mat h.ep 1 el 0 1 1 3 mat h.ep 2 el 0 h.ep 2 el 1 2 2 mat h.ep 2 el 0 2 1 3 mat h.ep 3 el 0 h.ep 3 el 1 2 2 mat h.ep 3 el 0 3 1 3 mat h.ep 4 el 0 h.ep 4 el 1 2 2 mat h.ep 4 el 0 4 1 3 mat h.ep 5 el 0 h.ep 5 el 1 2 2 mat h.ep 5 el 0 5 1 3 mat 1 0 0 0 0 1 3 2 mat 0.05 - 0 0 10 red 1 0 1 10 red 0.06 - 1 1 10 red 1 0.06 h 10 red 0.075 - 1.06 Δp;q 10 red 5 5 mat 21 1 sys

Линейная характеристика Ось x - положение регулятора зелёная кривая - падение напора на клапане синяя кривая - собственная хар-ка клапана красная кривая - действит. хар-ка клапана

x EPline KvTrim x EPline q' x EPline dp'.кл h.l1 h.l1 EPline q' o 10 red h.l1 h.l1 EPline dp'.кл o 10 green 2 5 mat h.l2 h.l2 EPline q' x 10 red h.l2 h.l2 EPline dp'.кл x 10 green 2 5 mat 0 h.l1 EPline dp'.кл h.l1 h.l1 EPline dp'.кл 2 2 mat 0 h.l1 EPline q' h.l1 h.l1 EPline q' 2 2 mat 0 h.l2 EPline dp'.кл h.l2 h.l2 EPline dp'.кл 2 2 mat 0 h.l2 EPline q' h.l2 h.l2 EPline q' 2 2 mat h.l 1 el 0 h.l 1 el 1 2 2 mat h.l 1 el 0 1 1 3 mat h.l 2 el 0 h.l 2 el 1 2 2 mat h.l 2 el 0 2 1 3 mat h.l 3 el 0 h.l 3 el 1 2 2 mat h.l 3 el 0 3 1 3 mat h.l 4 el 0 h.l 4 el 1 2 2 mat h.l 4 el 0 4 1 3 mat h.l 5 el 0 h.l 5 el 1 2 2 mat h.l 5 el 0 5 1 3 mat 1 0 0 0 0 1 3 2 mat 0.05 - 0 0 10 red 1 0 1 10 red 0.06 - 1 1 10 red 1 0.06 h 10 red 0.075 - 1.06 Δp;q 10 red 5 5 mat 21 1 sys

Действительный авторитет a 1 EPline dp.кл Δp.c / :

a #

h.ep #

Δp.c q 1 el t 1 el ΣΔp - мH2O #

Δp.c q 2 el t 2 el ΣΔp - мH2O #

Δp.c q 3 el t 3 el ΣΔp - мH2O #

h.l #

Δp.c q 4 el t 4 el ΣΔp - мH2O #

Δp.c q 5 el t 5 el ΣΔp - мH2O #

h.ep min #

h.ep max #

h.l min #

h.l max #

Проверка на кавитацию

1. Клапан после ТО

n P.изб 130 мH2O * 64 мH2O * - ΔP.доЦТП n el 2 / - q n el t n el ΣΔp - :

k 1 n range Кав_зап k 1 el 0.2 k P.изб t.обр k el Δp.кл_max : Кав_зап k 2 el Δp.c q k el t k el ΣΔp - : Кав_зап k 3 el Кав_зап k 1 el Кав_зап k 2 el - : 31 line for

Кав_зап мH2O #

1 P.изб мH2O #

q 1 el t 1 el ΣΔp мH2O #

2 P.изб мH2O #

q 2 el t 2 el ΣΔp мH2O #

3 P.изб мH2O #

q 3 el t 3 el ΣΔp мH2O #

4 P.изб мH2O #

q 4 el t 4 el ΣΔp мH2O #

5 P.изб мH2O #

q 5 el t 5 el ΣΔp мH2O #

2. Клапан перед ТО

n P.изб 130 мH2O * 64 мH2O * - ΔP.доЦТП n el 2 / - :

k 1 n range Кав_зап k 1 el 0.2 k P.изб t k el Δp.кл_max : Кав_зап k 3 el Кав_зап k 1 el Кав_зап k 2 el - : 21 line for

Кав_зап мH2O #

1 P.изб мH2O 64.6049

2 P.изб мH2O 65.4182

3 P.изб мH2O 63.1629

4 P.изб мH2O 65.4005

5 P.изб мH2O 65.1595

Расчёт спускника

Диаметр водовода D 720 mm * 29 * mm * - :

Диаметр спускника D.сп 219 mm * 28 * mm * - :

l.сп 1 m * :

l 58.5 m * :

H.1 64.55 m * :

H.2 58.35 m * :

H H.1 H.2 - 720 mm * + l.сп + :

D.сп 37 mm * :

H 7.92 m *

Kv.тр 150 L s / * D l Δ.э 40 °C * η.w ζ.l D Kv :

Kv.тр 4.3583 m 3^s / *

ζ.тр_сп 150 L s / * D.сп l.сп Δ.э 40 °C * η.w ζ.l :

http://armtek-msk.ru/zadvizhka-30s41nzh ζ.30с41нж 0.8 :

Идельчик диаг. 11-1 - скоростной напор
(эквивалентно коэф. Кориолиса) ζ.изл 1 :

Спускник: труба + задвижка + излив в атмосферу Kv.сп ζ.тр_сп ζ.30с41нж + ζ.изл + D.сп Kv :

Kv.сп 0.3302 m 3^s / *

150 L s / * D.сп Δ.э 40 °C * η.w λ 0.0246

ζ.вн_суж 0.5 :

Kv.суж ζ.вн_суж D.сп Kv :

Kv.суж 0.6473 m 3^s / *

Осн. тр. + сужение + спускник Kv.Σ Kv.суж Kv.сп 12 mat Kv.посл :

Kv.Σ 0.2942 m 3^s / *

ζ.вн_суж ζ.тр_сп + ζ.30с41нж + 11 μ 0.6427

ρ.воз 1.24 kg m 3^/ * :

D.воз 45 mm * 23 * mm * - :

l.воз 1 m * :

η.воз 15.11 106 -^* m 2^s / * :

http://www.pnzprom.ru/valves/kzap/15s22njbpa21004.html ζ.15с22нж 4.5 :

Идельчик диаг. 6-1 (180°) ζ.изг 1.4 0.17 * 1 * :

ζ.тр_воз 150 L s / * D.воз l.воз Δ.э η.воз ζ.l :

ζ.воз ζ.15с22нж ζ.изг + ζ.тр_воз + ζ.изл + :

ζ.воз 6.6909

ζ.воз 0 :

Располагаемый напор (разность отметок) расходуется на преодоление гидравлических сопротивлений водой, воздухом и создание скоростного напора D.1 D.2 D.воз Kv.Σ ζ.воз ρ q Δp.спуск q Kv.Σ ρ Δp ρ 2 / q D.2 v (2^q D.1 v (2^- (* + ζ.воз q D.воз ρ.воз Δp + :

D.1 D.2 D.воз H.1 H.2 Kv.Σ ζ.воз q.спуск ρ.40 40 °C * ρ.w eval : Δp H.1 H.2 - ρ.40 Δp eval : q.max Kv.Σ Δp ρ.40 q eval s m 3^/ * : D.1 D.2 D.воз Kv.Σ ζ.воз ρ.40 q m 3^s / * Δp.спуск Δp ≡ q 0 q.max solve m 3^s / * 41 line :

D.воз H.1 H.2 Kv.Σ ζ.воз ρ.вод ρ.воз q.спуск H.1 H.2 - ρ.вод Δp.ρ 1 bar * ρ.вод * Kv.Σ 2^ρ.w_max * / 8 ζ.воз * ρ.воз * π 2^D.воз 4^* / + / sqrt :

H.1 D + H.2 l.сп - (- H.2 H.2 l.сп - (- / 7.902

H1 H.1 D + :

H2 H.2 l.сп - :

ρ.w40 40 °C * ρ.w :

q.нач D.воз H1 H2 Kv.Σ ζ.воз ρ.w40 ρ.воз q.спуск :

q.нач 0.137744735730105 m 3^s / *

q.нач D.сп v 4.2559 m s / *

q.нач D v 0.3559 m s / *

Kv.Σ Kv.суж Kv.сп 12 mat Kv.посл :

q.кон D.воз H.2 H2 Kv.Σ ζ.воз ρ.w40 ρ.воз q.спуск :

q.кон 0.049001175381122 m 3^s / *

q.кон D.сп v 1.514 m s / *

q.нач q.кон / 2.811

V l π 4 / * D 2^* :

V q.нач q.кон + (0.5 * / min 4.0415

V 22.6423 m 3^*

ζ.воз q.нач D.воз ρ.воз Δp 55155.7568 Pa *

q.нач Kv.Σ ρ.w40 Δp 21767.1619 Pa *

ζ.воз q.нач D.воз ρ.воз Δp мH2O 5.6243

V 350 m 3^* :

V q.нач q.кон + (0.5 * / hr 1.0412

q.нач Kv.Σ ρ.w40 Δp мH2O 2.2196

ρ.w40 2 / q.нач D.сп v (2^q.нач D v (2^- (* 8927.0233 Pa *

H.1 D + H.2 l.сп - (- ρ.w40 Δp.ρ 76922.9187 Pa *

ρ.w40 2 / q.нач D.сп v (2^q.нач D v (2^- (* мH2O 0.9103

H.1 D + H.2 l.сп - (- ρ.w40 Δp.ρ мH2O 7.844

МУ 34-747-76 Указания по проектированию стальных трубопроводов гидротехнических сооружений, прил. XIII

Внутренний диаметр основного трубопровода, м D.0 D :

Внутренний диаметр трубопровода опорожнения, м d.0 D.сп :

Общая длина опоражниваемых участков, считая по оси трубопровода, м L.1 l :

L.2 0.00001 m * :

L.3 0.00001 m * :

Длина трубопровода опорожнения l.0 l.сп :

Коэффициент трения для определения потери напора в трубопроводе опорожнения λ.0 0.03 :

H.1 H1 H2 - :

H.1 7.902 m *

H.2 0.1 m * :

a D.0 d.0 / (4^1 λ.0 l.0 d.0 / * + (* 1.3 D.0 d.0 / (3.4^* + 2 g.e * / :

k L.2 L.1 + L.3 / :

k 5.85 106^*

c D.0 d.0 / (4^k 2^* 1 λ.0 l.0 d.0 / * + (* 2.6 k 2^* D.0 d.0 / (3.4^* + 2 g.e * / :

t.1 2 L.1 * a sqrt * H.1 / H.1 sqrt H.2 sqrt - (* 2 L.2 * H.2 sqrt / c sqrt * + :

t.1 1674.7079 sec *

t.1 min 27.9118

v.нач1 H.1 a / sqrt :

v.нач1 0.7835 m s / *

t.2 2 L.1 H.1 sqrt / * c sqrt * :

v.нач2 H.1 c / sqrt :

v.нач2 1.1528 107 -^* m s / *

t.3 2 L.1 H.1 sqrt / * a sqrt * :

t.3 149.3247 sec *

t.3 min 2.4887

T 0102030405060708090100110120130140150160170118 mat transpose :

Абсолютное давление насыщения водяного пара, [Па] P 611.2 1228.2 2339.2 4246.7 7384.4 1235119946312014741570182101420143380198670270260361500476100618140792050118 mat transpose :

t p T 0.001 P * t ainterp eval :

25 p 3.175261

x p

Create clean vector/matrices PT x 01701 range : i 1 x rows range pt i el T 0.001 P * x i el ainterp eval : for x pt augment 31 line :

PT 0 0.6112 1 0.65067 2 0.69508 3 0.74443 4 0.79872 5 0.85795 6 0.92212 7 0.99123 8 1.06528 9 1.14427 10 1.2282 11 1.317239 12 1.411481 13 1.510814 14 1.615124 15 1.724299 16 1.838226 17 1.956793 18 2.079886 19 2.207392 20 2.3392 21 2.480559 22 2.636015 23 2.804398 24 2.984536 25 3.175261 26 3.3754 27 3.583784 28 3.799242 29 4.020604 30 4.2467 31 4.485807 32 4.745888 33 5.025299 34 5.322398 35 5.635541 36 5.963085 37 6.303386 38 6.654801 39 7.015687 40 7.3844 41 7.772267 42 8.190231 43 8.636072 44 9.107568 45 9.602501 46 10.11865 47 10.653794 48 11.205714 49 11.772189 50 12.351 51 12.957127 52 13.605093 53 14.291988 54 15.014906 55 15.770937 56 16.557173 57 17.370706 58 18.208627 59 19.068028 60 19.946 61 20.861744 62 21.833941 63 22.858899 64 23.93293 65 25.052343 66 26.213449 67 27.412557 68 28.645978 69 29.910023 70 31.201 71 32.542892 72 33.959092 73 35.445027 74 36.996123 75 38.607807 76 40.275507 77 41.99465 78 43.760661 79 45.568969 80 47.415 81 49.327981 82 51.336508 83 53.435062 84 55.618121 85 57.880168 86 60.21568 87 62.61914 88 65.085026 89 67.60782 90 70.182 91 72.842434 92 75.623303 93 78.518057 94 81.520147 95 84.623025 96 87.82014 97 91.104945 98 94.470889 99 97.911424 100 101.42 101 105.037373 102 108.803596 103 112.71107 104 116.752195 105 120.91937 106 125.204996 107 129.601471 108 134.101197 109 138.696574 110 143.38 111 148.1984 112 153.197953 113 158.369943 114 163.705654 115 169.196372 116 174.83338 117 180.607962 118 186.511404 119 192.534988 120 198.67 121 204.96966 122 211.486458 123 218.210577 124 225.132203 125 232.241523 126 239.528721 127 246.983982 128 254.597492 129 262.359436 130 270.26 131 278.358777 132 286.714601 133 295.316514 134 304.153558 135 313.214777 136 322.489215 137 331.965912 138 341.633914 139 351.482262 140 361.5 141 371.752958 142 382.306247 143 393.14783 144 404.265669 145 415.647727 146 427.281968 147 439.156352 148 451.258844 149 463.577406 150 476.1 151 488.87803 152 501.971326 153 515.380494 154 529.106141 155 543.148873 156 557.509294 157 572.188012 158 587.185631 159 602.502759 160 618.14 161 634.09685 162 650.3724 163 666.96665 164 683.8796 165 701.11125 166 718.6616 167 736.53065 168 754.7184 169 773.22485 170 792.05 1712 mat

this is not a function definition

x D x p (x diff :

diff((p(x)),x) is an analytical operator
It does not evaluate  nor plot  because it
it not a "numerical scalar" operator and
none of: linterp, cinterp, ainterp are
scalar functions in Smath 25 D #

x p ( x diff