压水堆核电厂二回路热力系统课程设计

运行结果如下：

Ne=10^6; %发电功率（KW）

Ed=0.0105; %排污率

n1=0.994; %一回路能量利用系数

nli=0.89; %低压缸内效率

nm=0.98; %汽轮机组机械效率

nge=0.98; %发电机效率

nh=0.98; %加热器效率

Hfwp=6.7; %给水泵扬程

pw=918.19; %给水密度

nfwpp=0.58; %给水泵效率

nfwpti=0.8; %给水泵汽轮机内效率

nfwptm=0.9; %给水泵汽轮机机械效率

nfwptg=0.98; %给水泵汽轮机减速器效率

Hfh=2774.9846; %蒸汽发生器出口新蒸汽比焓（kj/kg）

Hg=842.86; %蒸汽发生器给水比焓（kj/kg）

Hps=1241.17; %蒸汽发生器运行压力下的饱和水焓（kj/kg）

xhz=0.8615; %高压缸排气干度

xz1=0.995; %第一级再热器的干度

Hho=2480.12; %高压缸出口焓

Hhi=2774.9782; %高压缸进口蒸汽焓

Hli=2988.3463; %低压缸进口焓

Hlo=2296.4958; %低压缸出口焓

Hsrh1=2614.86; %汽水分离再热器第一级抽气焓

Hzs1=1018.9756; %汽水分离再热器第一级疏水焓

Hsrh2=2783.348; %汽水分离再热器第二级抽气焓

Hzs2=1219.8773; %汽水分离再热器第二级疏水焓

Hfw=707.227; %分离器分离水比焓

h1=99.46; %高低压加热器中的平均焓升

h2=117.6; %再热器中的平均焓升

Hfwi=[146.6448 246.10 345.56 445.02 544.48 643.94 743.40]; %各级加热器进口焓（kj/kg）

Hfwo=[246.10 345.56 445.02 544.48 643.94 743.40 842.86]; %各级加热器出口焓(kj/kg) Hc=[2377.40 2525.49 2666.18 2803.096 2480.12 2448.787 2531.89]; %每一级的给水加热器的加热蒸汽比焓

Hw=[254.06 353.61 453.235 552.88 0 753.57 854.18]; %每一级给水加热器疏水比焓 H=100;

for Gslp=100:3000

if abs(H-Hfwo(5))/Hfwo(5)<0.005

H=(X*Hfwo(4)+(Ghes(1)+Ghes(2)+Gzc1+Gzc2)*Hw(6)+Gsdea*Hho+Hfw*Gslp*(xz1-xhz)/xhz)/Gfw; break

else

end

for Gfw=Gslp:3000;

Nfwpp=1000*Gfw*Hfwp/pw; %给水泵有效功率

Nfwpt=Nfwpp/(nfwpp*nfwptm*nfwptg); %给水泵实际功率

Gsfwp=Nfwpt/(Hfh-Hho); %给水泵耗气量

Ds=Gfw/(1+Ed); Gd=Gfw-Ds; %Ds为新蒸汽气量，Gd为排污量

X=Gslp+Gd+Gsfwp; %假设给水量

Gles(4)=h1*X/(Hc(4)-Hw(4))/nh;

Gles(3)=(h1*X-nh*Gles(4)*(Hw(4)-Hw(3)))/(Hc(3)-Hw(3))/nh;

Gles(2)=(h1*X-nh*(Gles(3)+Gles(4))*(Hw(3)-Hw(2)))/(Hc(2)-Hw(2))/nh;

Gles(1)=(h1*X-nh*(Gles(2)+Gles(3)+Gles(4))*(Hw(2)-Hw(1)))/(Hc(1)-Hw(1 ))/nh;

w1=(x-Gles(1)-Gles(2)-Gles(3)-Gles(4))*(Hli-Hlo)+Gles(1)*(Hli-Hc(1))+Gles(2)*(Hli-Hc(2))+Gles(3)*(Hli-Hc(3))+Gles(4)*(Hli-Hc(4));

W1=w1*nm*nge;

W2=10^6-W1; %上述程序Gles指的是低压给水再热器抽气，括号内数字为级数。W1和W2为高低压缸发电量

Gsrh1=Gslp*h2/(nh*(Hsrh1-Hzs1)); %第一级再热抽汽量

Gsrh2=Gslp*h2/(nh*(Hsrh2-Hzs2)); %第二级再热抽汽量

Ghes(2)=(Gfw*h1-nh*Gsrh2*(Hzs2-Hw(7)))/(nh*(Hc(7)-Hw(7)));

Ghes(1)=(Gfw*h1-nh*((Ghes(2)+Gsrh2)*(Hw(7)-Hw(6))+Gsrh1*(Hzs1-Hw(6))))/(nh*(Hc(6)-Hw(6)));

Gsdea=Gfw-X-Ghes(1)-Ghes(2)-Gsrh1-Gsrh2-Gslp*(xz1-xhz)/xhz;

Gt=Gsdea+Gslp*xz1/xhz; %上述Ghes为高压给水再热器抽气量，括号内数字为级数。Gsdea为除氧器耗气量

w2=Gt*(Hfh-Hho)+Ghes(1)*(Hfh-Hc(6))+Ghes(2)*(Hfh-Hc(7))+Gsrh1*(Hfh-Hzc1);

U2=w2*nm*nge; %高压缸实际发电量

if abs(W2-U2)/W2<0.001 %根据任务书要求其误差小于0.1%时迭代终止，输出结果 H=(X*Hfwo(4)+(Ghes(1)+Ghes(2)+Gsrh1+Gsrh2)*Hw(6)+Gsdea*Hho+Hfw*Gslp*(xz1-xhz)/xhz)/Gfw;

break

else

end

end

end

Qr=(Ds*(Hfh-Hg)+Gd*(Hps-Hg))/n1;

Nepp=Ne/Qr;

Gshp=Gt+Ghes(1)+Ghes(2)+Gsrh1;

disp('核电厂效率')

Nepp

disp('反应堆热功率 MW')

Qr

disp('蒸汽发生器总蒸汽产量 kg/s')

Ds

disp('汽轮机高压缸耗气量 kg/s')

Gshp

disp('汽轮机低压缸耗气量 kg/s')

Gslp

disp('第一级再热器耗气量 kg/s')

Gsrh1

disp('第二级再热器耗气量 kg/s')

Gsrh2

disp('除氧器耗气量 kg/s')

Gsdea

disp('给水泵汽轮机耗气量 kg/s')

Gsfwp

disp('给水泵给水量 kg/s')

Gfw

disp('给水泵扬程 MPa')

Hfwp

disp('第1级高压给水加热器抽汽量kg/s')

Ghes(1)

disp('第2级高压给水加热器抽汽量kg/s')

Ghes(2)

disp('第1级低压给水加热器抽汽量kg/s')

Gles(1)

disp('第2级高压给水加热器抽汽量kg/s')

Gles(2)

disp('第3级高压给水加热器抽汽量kg/s')

Gles(3)

disp('第4级高压给水加热器抽汽量kg/s') Gles(4)