рефераты рефераты
Домой
Домой
рефераты
Поиск
рефераты
Войти
рефераты
Контакты
рефераты Добавить в избранное
рефераты Сделать стартовой
рефераты рефераты рефераты рефераты
рефераты
БОЛЬШАЯ ЛЕНИНГРАДСКАЯ БИБЛИОТЕКА
рефераты
 
МЕНЮ
рефераты Расчет системы водоснабжения рефераты

БОЛЬШАЯ ЛЕНИНГРАДСКАЯ БИБЛИОТЕКА - РЕФЕРАТЫ - Расчет системы водоснабжения

Расчет системы водоснабжения

Московский Государственный Университет Путей Сообщения

Институт Транспортной Техники и Организации Производства

(ИТТОП)

Кафедра: «Теплоэнергетика железнодорожного транспорта»

(ТЖТ)

Курсовая работа по дисциплине:

«Технологические энергоносители промышленных предприятий»

Выполнила:

Ст. гр. ТЭН-412

Проверил:

Москва 2009

Расчет системы водоснабжения

Необходимо:

- Построить график нагрузки сети и установить параметры режимов - максимального водопотребления и максимального транзита,

- Выбрать диаметры труб для участков сети, исходя из режима максимального транзита,

- Выполнить внутреннюю увязку сети для двух указанных режимов,

- Построить пьезометрические графики сети и вычислить характеристики водопитателей.

График нагрузки сети

Строится как сумма графиков потребления воды всеми потребителями. Режим максимального водопотребления соответствует часам максимальной нагрузки сети. В режиме максимального транзита нагрузки сети - минимально.

Потребитель

Расход, м3/с

Геодезическая высота

0-6

6-12

12-18

18-24

1

0,02

0,08

0,05

0,03

3,5

2

0,01

0,13

0,09

0,03

5,0

3

0,04

0,10

0,12

0,05

1,5

4

0,02

0,04

0,05

0,01

4,0

5

0,02

0,07

0,07

0,03

5,5

НС

0,5

Б

10,5

?

0,11

0,42

0,38

0,15

Qср - среднесуточная нагрузка системы определяется интегрированием графика нагрузки сети по времени

=(Q0-6*6+Q6-12*6+Q12-18*6+Q18-24*6)/24=(0,11+0,42+0,38+0,15)/4= 1,06/4=0,265 м3/с

Эта величина равна производительности насосной станции при условии, что в течении суток вода подается насосами равномерно.

Поток воды в водонапорную башню при минимальной нагрузке сети Qсmin =0,11 м3/с (режим максимального транзита)

QБ= QН - Qсmin=0,265 - 0,11 = 0,155 м3/с

В часы максимального водопотребления башня отдает воду в сеть в количестве:

QБ=Qcmax - QНС=0,42 - 0,265=0,155 м3/с,

где Qcmax =0,42 м3/с - максимальная нагрузка сети

Выбор диаметров труб для участков сети

Режим максимального транзита (Qс=min) является определяющим для выбора диаметров. Вначале задается начальное потокораспределение, которое удовлетворяет первому закону Кирхгофа. Для этого, рассматривая конфигурацию сети, как граф, рекомендуется построить дерево графа путем удаления части участков сети, число которых равно n

n=p-m+1 (1),

где p- число участков, m - число узлов, n - количество элементарных контуров

n=9-7+1=3

В результате сеть становится тупиковой (контуры в сети отсутствуют, n=0). Для нее просто задать начальное потокораспределение путем суммирования потребления воды, начиная с концевых узлов. На исключенных участках поток воды принимается равным нулю.

Для каждого i-того участка дерева диаметр труб выбирается из таблицы 1 по величине приведенного расхода qiп

,

где qi - поток на участке, м3/с

Э - действительный экономический фактор,

Эт - табличный экономический фактор,

х - доля единичного расхода на участке,

Единичный расход на участке:

Qс - нагрузка сети для рассматриваемого режима, м3/с (Qс=min):

1уч. НС-5 qп=0,265*(0,11*0,3/0,265)1/3= 0,198 м3/с=198 литр/с > D=450мм=0,45м

2уч. 5-2 qп=0,245*(0,11*1/0,245) 1/3=0,189 м3/с=189 литр/с > D=450мм=0,45м

3уч. 2-Б qп=0,155*(0,11*0,7/0,155) 1/3=0,123 м3/с=123 литр/с > D=350мм=0,35м

4уч. 2-1 qп=0,08*(0,11*0,3/0,08) 1/3=0,060 м3/с=60 литр/с > D=250мм=0,25м

5уч. 1-3 qп=0,06*(0,11*0,3/0,06) 1/3=0,040 м3/с=49 литр/с > D=250мм=0,25м

6уч. 3-4 qп=0,02*(0,11*0,3/0,02) 1/3=0,024 м3/с=24 литр/с > D=175мм=0,175м

Диаметры исключенных участков принимаются равными средним значениям диаметров тех участков, которые они соединяют.

7уч. Б-1 D=(D2-Б+D2-1)/2=(0,35+0,25)/2=0,3м

8уч. 3-2 D=(D1-2+D1-3)/2=(0,25+0,25)/2=0,25м

9уч. 5-4 D=(D2-5+D3-4)/2=(0,45+0,175)/2=0,3125?0,3м

Режим минимального транзита (Qс=max) :

1уч. НС-5 qп=0,265*(0,42*1/0,265)1/3= 0,309 м3/с=309 литр/с

2уч. 5-2 qп=0,245*(0,42*1/0,245) 1/3=0,293 м3/с=293 литр/с

3уч. 2-Б qп=0,155*(0,42*0,7/0,155) 1/3=0,192 м3/с=192 литр/с

4уч. 2-1 qп=0,08*(0,42*0,3/0,08) 1/3=0,093 м3/с=93 литр/с

5уч. 1-3 qп=0,06*(0,42*0,3/0,06) 1/3=0,077 м3/с=77 литр/с

6уч. 3-4 qп=0,02*(0,42*0,3/0,02) 1/3=0,037 м3/с=37 литр/с

Внутренняя увязка сети:

Увязка выполняется по методу Лобачева-Кросса. Для каждого участка сети вычисляется гидравлическая характеристика si

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

Падение давления на участке связано с расходом квадратической зависимостью

pi=siqi2

где pi - падения давления, м вод ст; si - гидравлическая характеристика; qi - расход воды на участке, м3/с

По соотношению Эйлера (1) определяется количество элементарных контуров сети. Увязка выполняется следующих образом:

1. Для каждого участка сети задается положительное направление потока и его величина, причем все потоки на участках должны удовлетворять первому закону Кирхгофа; для этого удобно рассматривать его как нулевое приближение,

2. В каждом контуре вычисляется невязка ?hn, определяющая меру несоответствия потоков воды в сети второму закону Кирхгофа,

?hn=?siqiabs(qi)

Суммирование выполняется по участкам сети, составляющим контур, с учетом знаков и направления обхода, принятых при формулировке второго закона Кирхгофа,

3. Если модуль максимальной невязки меньше допустимой величины, увязка заканчивается; в противном случае процесс итераций продолжается,

4. Для каждого контура вычисляется поправка ?qn, положительное направление которой противоположно обходу контура по второму закону Кирхгофа

Суммирование в знаменателе выполняется по участкам контура,

5. Контурные поправки используются для коррекции значений потоков воды на участках

qik+1=qik+??qn

Суммируются поправки только тех контуров, в состав которых входит рассматриваемый участок; суммирование выполняется с учетом направлений, принятых для потоков и контурных поправок; далее вычисления повторяются, начиная с пункта 2.

По программе получилось, что выбранные направления на участках 5,6,8,9 надо изменить на противоположные. Это будет выглядеть так:

Пьезометрический график

Номер узла

Геодез. Высота Z, м

Свободный Напор (предвар) Н,м

Пьез. Высота (предвар) П, м

Н, м

П, м

НС

0,5

32,636

33,136

63,8

64,3

1

3,5

1,504

5,004

21,5

25 (дикт)

2

5,0

19

24

38,996

43,996

3

1,5

21,728

23,228

41,724

43,224

4

4,0

20

24 (дикт)

51,164

55,164

5

5,5

18,453

23,953

49,617

55,117

Б

10,5

-5,507

4,993

14,489

24,989

Выбираю диктующего потребителя: 4. Задаю П4=24 м, H4=П4-z4=24-4=20 м

П5=П4-?p4-5=24-0,047=23,953

П3=П4-?p3-4=24-0,772=23,228

П1=П3-?p1-3=23,228-18,224=5,004

П2=П3+?p2-3=23,228+0,772=24

ПНС=П5+?pНС-5=23,953+9,183=33,136

ПБ=П1-?pБ-1=5,004-0,001=5,003

ПБ=П2-?pБ-2=24-19,017=4,983 ПБ?4,993

ННС=ПНС-zНС=33,136-0,5=32,636

Н1=П1-z1=5,004-3,5=1,504

Н2=П2-z2=24-5=19

Н3=П3-z3=23,228-1,5=21,728

Н5=П5-z5=23,953-5,5=18,453

НБ=ПБ-zБ=4,993-10,5= -5,507

Минимальный свободный напор получился для потребителя 1

Принимаю П1=25м

П3=П1+?p1-3=25+18,244=43,224

П2=П3+?p2-3=43,224+0,772=43,996

П4=П3+?p3-4=43,224+11,940=55,164

П5=П4-?p4-5=55,164-0,047=55,117

ПНС=П5+?pНС-5=55,117+9,183=64,3

ПБ=П1-?pБ-1=25-0,001=24,999

ПБ=П2-?pБ-2=43,996-19,017=24,979 ПБ?24,989

ННС=ПНС-zНС=64,3-0,5=63,8

Н1=П1-z1=25-3,5=21,5

Н2=П2-z2=43,996-5=38,996

Н3=П3-z3=43,224-1,5=41,724

Н4=П4-z4=55,164-4=51,164

Н5=П5-z5=55,117-5,5=49,617

НБ=ПБ-zБ=24,989-10,5=14,489

Характеристики водопитателей

Рабочий объем водонапорной башни Vp находится с помощью графика нагрузки сети. Необходимо вычислить интеграл

Vp=?(Qc(ф)-Qc)d ф=6*[(0,42-0,265)+(0,38-0,265)]=1,62 м3

Если диаметр бака башни и максимальный уровень его заполнения относятся как 2:1, то с учетом 10% запаса воды диаметр бака равен

DБ=1,46Vp1/3=1,46*1,61/3=1,715 м,

h0=,5*DБ=0,8575 м

где h0 - максимальный уровень заполнения бака водой,

DБ - диаметр бака

Vp - рабочий объем

H'НС=Нтр+(zA-zНС)+??p'i=20+(3,5-,5)+89,434=112,434

H”НС=HБ+(zБ-zНС)+h0+??p''i=14,489+(10,5-0,5)+0,8575+36,907=62,2535

Где ?p' падения давления при Qc=max, ?p'' падения давления при Qc=min, а

HБ, H'НС - max ; H”НС - min

Итоговая таблица

Уч

D, м

L, м

q', м3/с

q'',м3/с

?p', м вод ст

?p'', м вод ст

1 (НС-5)

0,45

700

0,309

0,198

9,183

3,771

2 (5-2)

0,45

900

0,301

0,194

11,215

4,671

3 (2-Б)

0,35

1000

0,191

0,123

19,017

7,817

4 (2-1)

0,25

800

0,088

0,056

19,015

7,817

5 (1-3)

0,25

850

0,083

0,053

18,244

7,472

6 (3-4)

0,175

700

0,029

0,019

11,940

4,995

7 (Б-1)

0,3

650

0,001

0,000

0,001

0,000

8 (3-2)

0,25

1200

0,014

0,010

0,772

0,344

9 (5-4)

0,3

600

0,008

0,005

0,047

0,020

?

89,434

36,907

Где ?p' и q'- падения давления и расход при Qc=max, а ?p'' и q” - падения давления и расход при Qc=min

Газопровод среднего давления

потребитель

Расход по часам суток, куб.м /час

 

0-6

6-12

12-18

18-24

1

150

550

400

200

2

200

950

750

250

3

150

400

800

100

4

100

500

550

200

5

250

850

850

150

6

100

150

150

100

7

100

200

250

100

8

150

350

250

200

9

200

100

200

100

сумма

1400

4050

4200

1400

Газопровод низкого давления

потребитель

Расход по часам суток, куб.м /час

 

0-6

6-12

12-18

18-24

6

100

150

150

100

7

100

200

250

100

8

150

350

250

200

9

200

100

200

100

сумма

550

800

850

500

Диаметры труб участков сети выбираются с помощью номограммы, исходя из найденного оптимального pi , а также Li , и расхода газа Vi

Участок

Длина, м

Расход, м3/ч

Удельные потери, Па/м

Диаметр стандартный, м

1 - ГРП1

900

1 - 2

900

2 - 3

1000

2 - 4

550

2 - 5

750

5 - ГРП2

500

6 - ГРП2

250

6 - 7

400

6 - 8

150

9 - ГРП 2

250

Испр_max

N=9

DIM flux(N), diam(N), length(N),s(N), press(N)

READ dHmin, a

DATA 0.001, 0.15

for i=1 to N

READ b

flux(i)=b

READ b

length(i)=b

READ b

diam(i)=b

NEXT i

DATA 0.309, 700, 0.450

DATA 0.293,900, 0.450

DATA 0.192,1000, 0.350

DATA 0.093, 800,0.250

DATA 0.077, 850,0.250

DATA 0.037, 700, 0.175

DATA 0, 650, 0.300

DATA 0, 1200, 0.250

DATA 0, 600, 0.300

for i=1 to N

s(i)=0.001735*(1+a)*length(i)/diam(i)^5.3

next i

[A]

dH1=s(3)*flux(3)*Abs(flux(3))-s(7)*flux(7)*Abs(flux(7))-s(4)*flux(4)*Abs(flux(4))

dH2=s(4)*flux(4)*Abs(flux(4))-s(5)*flux(5)*Abs(flux(5))-s(8)*flux(8)*Abs(flux(8))

dH3=s(2)*flux(2)*Abs(flux(2))+s(8)*flux(8)*Abs(flux(8))-s(6)*flux(6)*Abs(flux(6))-s(9)*flux(9)*Abs(flux(9))

if Abs(dH1)<dHmin and Abs(dH2)<dHmin and Abs(dH3)<dHmin goto [B]

dq1=0.5*dH1/(s(4)*Abs(flux(4))+s(7)*Abs(flux(7))+s(3)*Abs(flux(3)))

dq2=0.5*dH2/(s(4)*Abs(flux(4))+s(5)*Abs(flux(5))+s(8)*Abs(flux(8)))

dq3=0.5*dH3/(s(8)*Abs(flux(8))+s(2)*Abs(flux(2))-s(9)*Abs(flux(9))+s(6)*Abs(flux(6)))

flux(1)=flux(1)

flux(2)=flux(2)-dq3

flux(3)=flux(3)-dq1

flux(4)=flux(4)+dq1-dq2

flux(5)=flux(5)+dq2

flux(6)=flux(6)+dq3

flux(7)=flux(7)+dq1

flux(8)=flux(8)+dq2-dq3

flux(9)=flux(9)-dq3

goto [A]

[B]

for i=1 to N

press(i)=s(i)*flux(i)*Abs(flux(i))

next i

for i=1 to N

PRINT "i="; i;

PRINT ";lenght="; length(i);

PRINT ";diam="; diam(i);

PRINT ";flux="; using ("#.###", flux(i));

PRINT ";press="; using ("##.###", press(i))

next i

PRINT "dH1="; using ("#.#####", dH1) ;";dH2="; using ("##.#####",dH2); ";dH3="; using("#.#####", dH3)

end

Испр_min

N=9

DIM flux(N), diam(N), length(N),s(N), press(N)

READ dHmin, a

DATA 0.001, 0.15

for i=1 to N

READ b

flux(i)=b

READ b

length(i)=b

READ b

diam(i)=b

NEXT i

DATA 0.198, 700, 0.450

DATA 0.189,900, 0.450

DATA 0.123,1000, 0.350

DATA 0.060, 800,0.250

DATA 0.049, 850,0.250

DATA 0.024, 700, 0.175

DATA 0, 650, 0.300

DATA 0, 1200, 0.250

DATA 0, 600, 0.300

for i=1 to N

s(i)=0.001735*(1+a)*length(i)/diam(i)^5.3

next i

[A]

dH1=s(3)*flux(3)*Abs(flux(3))-s(7)*flux(7)*Abs(flux(7))-s(4)*flux(4)*Abs(flux(4))

dH2=s(4)*flux(4)*Abs(flux(4))-s(5)*flux(5)*Abs(flux(5))-s(8)*flux(8)*Abs(flux(8))

dH3=s(2)*flux(2)*Abs(flux(2))+s(8)*flux(8)*Abs(flux(8))-s(6)*flux(6)*Abs(flux(6))-s(9)*flux(9)*Abs(flux(9))

if Abs(dH1)<dHmin and Abs(dH2)<dHmin and Abs(dH3)<dHmin goto [B]

dq1=0.5*dH1/(s(4)*Abs(flux(4))+s(7)*Abs(flux(7))+s(3)*Abs(flux(3)))

dq2=0.5*dH2/(s(4)*Abs(flux(4))+s(5)*Abs(flux(5))+s(8)*Abs(flux(8)))

dq3=0.5*dH3/(s(8)*Abs(flux(8))+s(2)*Abs(flux(2))-s(9)*Abs(flux(9))+s(6)*Abs(flux(6)))

flux(1)=flux(1)

flux(2)=flux(2)-dq3

flux(3)=flux(3)-dq1

flux(4)=flux(4)+dq1-dq2

flux(5)=flux(5)+dq2

flux(6)=flux(6)+dq3

flux(7)=flux(7)+dq1

flux(8)=flux(8)+dq2-dq3

flux(9)=flux(9)-dq3

goto [A]

[B]

for i=1 to N

press(i)=s(i)*flux(i)*Abs(flux(i))

next i

for i=1 to N

PRINT "i="; i;

PRINT ";lenght="; length(i);

PRINT ";diam="; diam(i);

PRINT ";flux="; using ("#.###", flux(i));

PRINT ";press="; using ("##.###", press(i))

next i

PRINT "dH1="; using ("##.#####", dH1) ;";dH2="; using ("##.#####",dH2); ";dH3="; using("##.#####", dH3)

end

РЕКЛАМА

рефераты НОВОСТИ рефераты
Изменения
Прошла модернизация движка, изменение дизайна и переезд на новый более качественный сервер


рефераты СЧЕТЧИК рефераты

БОЛЬШАЯ ЛЕНИНГРАДСКАЯ БИБЛИОТЕКА
рефераты © 2010 рефераты