西氣東輸課程設(shè)計(編輯修改稿)
2024-09-01 16:44 本頁面
【文章內(nèi)容簡介】 )。 %首站至倒數(shù)第二站的平均長度 LYP=Lb(Ldistance(NN1))。 %倒數(shù)第二站與末站間的距離 %判斷末站位置是否合適 if abs(LPJLYP)1e3 break。 end %改變壓比 if (LPJLYP0) XZYB=XZYB+。 else XZYB=。 end %首站壓比不改變 XZYB(1)=0。 %計算倒數(shù)第二站至末站的溫度分布 L=Ldistance(NN1)。 for j=1:1000 C=putlamda(Q(NN))*Z*T(NN1,j)*delta/(C0^2*d^5)。 P(NN1,j+1)=(P(NN1,j)^2C*1000*Q(NN)^2)^。 L=L+1000。 a=putK(L)*pi*D/(Q(NN)***Cp)。 T(NN1,j+1)=putT0(L)+(T(NN1,j)putT0(L))*exp(a*1000) ... Di/1e6/a/1000*(P(NN1,j)P(NN1,j+1))*(1exp(a*1000))。 if (L==Lb) break。 end end %更新末站進(jìn)站溫度和出站流量 T1(NN)=T(NN1,j)。 Q(NN+1)=Q(NN)putGU(Q(NN),T1(NN),putYB(NN)+XZYB(NN))。 Lb=round((ZCputLZ(Q(NN+1)))/1000)*1000。 end %末站位置position(NN)=Lb/Step。position(NN+1)=ZC/Step。distance(NN)=Lbsum(distance)+distance(NN)。distance(NN+1)=ZCLb。 %計算首站至門站的溫度壓力分布for i=1:(NN+1) P(i,1)=10e6。 %第i站出站壓力 T(i,1)=50+。 %第i站出站溫度 Q(i+1)=Q(i)putGU(Q(i),T1(i),putYB(i)+XZYB(i))。 for j=1:1000 C=putlamda(Q(i+1))*Z*T(i,j)*delta/(C0^2*d^5)。 P(i,j+1)=(P(i,j)^2C*1000*Q(i+1)^2)^。 L=position(i)*Stepdistance(i)+j*Step。 a=putK(L)*pi*D/(Q(i+1)***Cp)。 T(i,j+1)=putT0(L)+(T(i,j)putT0(L))*exp(a*1000) ... Di/1e6/a/1000*(P(i,j)P(i,j+1))*(1exp(a*1000))。 if (L==(position(i)*Step)) break。 end end T1(i+1)=T(i,j+1)。 P1(i+1)=P(i,j+1)。end %求壓力溫度分布aaaa=size(P)。M=0。MM=0。for i=1:aaaa(1) for j=1:aaaa(2) if (P(i,j)~=0) M=M+1。 PFB(M)=P(i,j)。 end if (T(i,j)~=0) MM=MM+1。 TFB(MM)=T(i,j)。 end endend %末端儲氣能力V=putMDCQNL(Q(NN+2))。 %輸出結(jié)果 for i=1:33 %壓氣站位置 if (i==1) BZ(i,1)=0。 else BZ(i,1)=position(i1)。 end %壓比 if (i==33) BZ(i,2)=0。 else BZ(i,2)=putYB(i)+XZYB(i)。 end %進(jìn)站壓力 BZ(i,3)=P1(i)。 %出站壓力 if (i==33) BZ(i,4)=PFB(max(size(PFB)))。 else BZ(i,4)=10e6。 end %進(jìn)站溫度 BZ(i,5)=T1(i)。 %出站溫度 if (i==33) BZ(i,6)=TFB(max(size(PFB)))。 else BZ(i,6)=T(i,1)。 end %耗氣量 if (i==33) BZ(i,7)=0。 else BZ(i,7)=Q(i)Q(i+1)。 end end 子函數(shù)：
醫(yī)療健康相關(guān)推薦
鄂ICP備17016276號-1