По параметрам пара перед турбиной определим значение энтальпии острого пара:

Определим теоретический располагаемый теплоперепад в первом отсеке:

Определим значение действительного теплоперепада в первом отсеке:

Определим энтальпию пара в первом регенеративном отборе:

По аналогии найдем значения теплоперепада в остальных отсеках и значения энтальпий в остальных отборах:

тепловая паровая турбина электростанция

По найденным значениям строем процесс расширения пара в турбине на h-s диаграмме (Рисунок 2).

Рисунок 2 - Процесс расширения пара в проточной части турбине

2. Определение значения расхода острого пара на турбину

Подставим известные параметры

,

где - расход пара в -ом отсеке.

D₁=D0;

D₂=D0-D1;

D₃=D₂-D2;

D₄=D₃-D3;

D₅=D₄-D4;

По найденной выполним оценку расчетной погрешности приближения

6. Текст программы

#include <vcl. h>

#pragma hdrstop

#include "hsdgr2. cpp"

#include "gaussjor. cpp"

#include "Unit1. h"

void __fastcall TForm1:: Button1Click (TObject *Sender)

{ float We,P0,t0,P1,P2,P3,P4,Pn,Pd,Dn,dDn,tvk,Pdv,txov,tcl,n0i,nem,h0,h1,H1,h2,H2,h3,H3,h4,H4,H5,hn,D0,Kreg,Gpv,Gxov,Hpr,Dpr1,Gpr1,Tcl,Txov1,Txov,Top,Tcp,Dp1,h1p1,H1pr,H11pr

,Td,We1,d,T1cp,T1kn,dTkn,Tkn,Tod,T1d,D1p2,Tpv,Tdv,dTdn,T1dn,dTpv,Dp3,Dp4,D11p2,Dp,T1pv,tpn,T11d,Tpn,T1pn,dTpn,tpv,Tdn,H31;

h0=eheatHPT (P0,t0);

H1=eheatH0 (P0,h0,P1) *n0i;

h1=h0-H1;

H2=eheatH0 (P1,h1,P2) *n0i;

h2=h1-H2;

H3=eheatH0 (P2,h2,P3) *n0i;

h3=h2-H3;

H4=eheatH0 (P3,h3,P4) *n0i;

h4=h3-H4;

H5=eheatH0 (P4,h4,Pn) *n0i;

hn=h4-H5;

Kreg=1.225;

D0=Kreg*We*1000/ (nem* (H1+H2+H3+H4+H5));

Gpv=D0+0.03*D0+0.015*D0+0.035*D0;

Hpr=eheatH1 (eheatTs (1.1*P0));

H1pr=eheatH1 (eheatTs (7));

H11pr=eheatH11 (eheatTs (7));

Tcl=eheatHPT (7,tcl);

Txov=4.19*txov;

Tcp=eheatHPT (P1,eheatTs (P1) - 3);

h1p1=eheatHPT (P1,eheatTs (P1) +12.5);

T1cp=eheatHPT (P2,eheatTs (P2) - 3);

Td=eheatH1 (eheatTs (P1));

T1d=eheatH1 (eheatTs (P2));

T1pv=eheatH1 (eheatTs (Pd));

dTpn= ( (1.3*P0-Pd) *eheatV1 (eheatTs (Pd)) *100) /0.7;

Tpv=T1pv+dTpn;

tpv=eheatT (1.3*P0,Tpv);

T11d=eheatHPT (P2,tpv+8);

T1dn=eheatHPT (P3,eheatTs (P3) - 5);

H31=eheatH1 (eheatTs (P3));

Tdv=eheatH1 (eheatTs (Pdv));

dTpn= ( (12-Pdv) *eheatV1 (Pdv) *100) /0.7;

Tdn=Tdv+dTdn;

tpn=eheatT (12,Tdn);

Tkn=eheatHPT (Pn,tvk);

dTkn= ( (5-Pn) *eheatV (Pn,eheatHPT (Pn,tcl)) *100) /0.7;

T1kn=Tkn+dTkn;

float A [11] [11] ={{1,1,0,0,0,0,0,0,0,0,0},

{H11pr,H1pr,0,0,0,0,0,0,0,0,0},

{0,0,0,Gpv, (h1p1-h1) *0.985,0,0,0,0,0,0},

{0,0,0,0, (h1p1-Td) *0.985,0,0,0,0,0,0},

{0,0,0,0, (Td-T1d) *0.985,Gpv, (h2-T1d) *0.985,0,0,0,0},

{0,0,0,0, (T1d-T11d) *0.985,-Gpv, (T1d-T11d) *0.985,0,0,0,0},

{1,0,1,0,1,0,1,1,1,1,dDn},

{H11pr*0.985,0,T1dn,0,T11d*0.985,0,T11d*0.985,h2*0.985,T1dn,T1dn,dDn*T1dn},

{0,0,Tdn-T1dn,0,0,0,0,0, (h3-H31) *0.985+ (Tdn-T1dn),Tdn-T1dn,dDn* (Tdn-T1dn) },

{0,H1pr-Tcl,Txov-Tdv,0,0,0,0,0,H31*0.985-Tdv,h4*0.985-Tdv,dDn* (T1kn-Tdv) },

{0,0,1,0,0,0,0,0,0,0, (dDn-1) }};

float B [11] ={0.035*D0,0.035*D0*Hpr,Gpv*Tcp,Gpv* (Tcp-T1cp),Gpv*T1cp,-Gpv*Tpv,Gpv,Gpv*T1pv,0,0,0.08*D0};

float X [11];

GaussJordan (A [0],B,11,X);

We1=D0*H1*nem+ (D0-X [4]) *H2*nem+ (D0-X [4] - X [6] - X [7]) *H3*nem+ (D0-X [4] - X [6] - X [7] - X [8]) *H4*nem+X [10] *H5*nem;

d= (We-We1) *100/We

7. Вывод результатов

Полный расход топлива:

,

где - расход острого пара, кг/с;

- энтальпия питательной воды перед котельным агрегатом;

- КПД котельного агрегата ;

- калорийная способность условного топлива, .

КПД брутто по выработке электроэнергии на ТЭЦ:

КПД брутто по выработке тепла на ТЭЦ:

Удельный расход на выработку электрической энергии:

Список используемых источников