Výpočet výkonu čerpadla

Ptop35={ 10.95,13.21,15.73,16.56,17.66,18.50,19.05,19.81,20.40};
Ptop50={11.50,14.40,16.00,16.45,16.77,17.40,17.97,18.85,19.40};
Ptop60={11.95,15.05,16.15,16.35,16.45,16.90,17.45,18.20,18.75};

COP35={2.71,3.14,3.62,3.80,3.99,4.40,4.60,4.75,4.90};
COP50={2.20,2.50,2.80,2.92,3.02,3.20,3.35,3.45,3.55};
COP60={1.90,2.10,2.32,2.40,2.45,2.58,2.70,2.80,2.85};

n1=Length[Ptop35];
n2=Length[Ptop50];
n3=Length[Ptop60];

teplota={-15,-7,2,5,7,10,15,20,25};

radekTop35[i_]:={teplota[[i]],Ptop35[[i]]};
radekTop50[i_]:={teplota[[i]],Ptop50[[i]]};
radekTop60[i_]:={teplota[[i]],Ptop60[[i]]};

radekCOP35[i_]:={teplota[[i]],COP35[[i]]};
radekCOP50[i_]:={teplota[[i]],COP50[[i]]};
radekCOP60[i_]:={teplota[[i]],COP60[[i]]};

DataPtop35=Table[radekTop35[i],{i,n1}];
DataPtop50=Table[radekTop50[i],{i,n2}];
DataPtop60=Table[radekTop60[i],{i,n3}];

DataCOP35=Table[radekCOP35[i],{i,n1}];
DataCOP50=Table[radekCOP50[i],{i,n2}];
DataCOP60=Table[radekCOP60[i],{i,n3}];

Ptopna35[Tvzduch_]=Interpolation[DataPtop35] [Tvzduch];
Ptopna50[Tvzduch_]=Interpolation[DataPtop50] [Tvzduch];
Ptopna60[Tvzduch_]=Interpolation[DataPtop60] [Tvzduch];

TopnyCOP35[Tvzduch_]=Interpolation[DataCOP35] [Tvzduch];
TopnyCOP50[Tvzduch_]=Interpolation[DataCOP50] [Tvzduch];
TopnyCOP60[Tvzduch_]=Interpolation[DataCOP60] [Tvzduch];

Pel35[Tvzduch_]=((Ptopna35[Tvzduch])/(TopnyCOP35[Tvzduch]));
Pel50[Tvzduch_]=((Ptopna50[Tvzduch])/(TopnyCOP50[Tvzduch]));
Pel60[Tvzduch_]=((Ptopna60[Tvzduch])/(TopnyCOP60[Tvzduch]));


ts=0.006;

grafTopny=Plot[{Ptopna35[Tvzduch],Ptopna50[Tvzduch],Ptopna60[Tvzduch]},{Tvzduch,-15,25},AxesLabel->{"(Subscript[θ, e]) (°C)","Topný výkon (kW)"},PlotStyle->{{Darker[Blue], Thickness[ts]},{Darker[Red], Thickness[ts]},{Darker[Green], Thickness[ts]}},AxesOrigin->{0,10},PlotLegends-> {"Teplota vody 35 °C","Teplota vody 50 °C","Teplota vody 60 °C"}, GridLines->Automatic]


grafCOP=Plot[{TopnyCOP35[Tvzduch],TopnyCOP50[Tvzduch],TopnyCOP60[Tvzduch]},{Tvzduch,-15,25},AxesLabel->{"(Subscript[θ, e]) (°C)","COP (-)"},PlotStyle->{{Darker[Blue], Thickness[ts]},{Darker[Red],Thickness[ts]},{Darker[Green],Thickness[ts]}},AxesOrigin->{0,0},PlotLegends-> {"Teplota vody 35 °C","Teplota vody 50 °C","Teplota vody 60 °C"}, GridLines->Automatic]

grafPel=Plot[{Pel35[Tvzduch],Pel50[Tvzduch],Pel60[Tvzduch]},{Tvzduch,-15,25},AxesLabel->{"(Subscript[θ, e]) (°C)","Dodaný výkon (el.) (kW)"},PlotStyle->{{Darker[Blue], Thickness[ts]},{Darker[Red],Thickness[ts]},{Darker[Green],Thickness[ts]}},AxesOrigin->{0,0},PlotLegends-> {"Teplota vody 35 °C","Teplota vody 50 °C","Teplota vody 60 °C"}, GridLines->Automatic]

radekCelyCOP35[i_]:={teplota[[i]],35.,COP35[[i]]};
radekCelyCOP50[i_]:={teplota[[i]],50.,COP50[[i]]};
radekCelyCOP60[i_]:={teplota[[i]],60.,COP60[[i]]};

radekCelyTop35[i_]:={teplota[[i]],35.,Ptop35[[i]]};
radekCelyTop50[i_]:={teplota[[i]],50.,Ptop50[[i]]};
radekCelyTop60[i_]:={teplota[[i]],60.,Ptop60[[i]]};

n=Length[teplota];

dataCOPcompl=Flatten[Table[{radekCelyCOP35[i],radekCelyCOP50[i],radekCelyCOP60[i]},{i,n}],1];
dataTopcompl=Flatten[Table[{radekCelyTop35[i],radekCelyTop50[i],radekCelyTop60[i]},{i,n}],1];

COPcompl[Tvzduch_,Tvoda_]=Interpolation[dataCOPcompl,InterpolationOrder->{2,2}][Tvzduch,Tvoda];
Ptopcompl[Tvzduch_,Tvoda_]=1000*Interpolation[dataTopcompl,InterpolationOrder->{2,2}][Tvzduch,Tvoda];

Pelcompl[Tvzduch_,Tvoda_]=((Ptopcompl[Tvzduch,Tvoda])/(COPcompl[Tvzduch,Tvoda]));

Stanovení bodu bivalence - 1.část

MernaZtrataProstupem=4311.78;
MernaZtrataVetranim=2388.13;
StaleTepelneZisky=35300;
Tin=20;
MinSolar=11000;
CelkovaZtrata=(MernaZtrataProstupem*(Tin-Tout)+MernaZtrataVetranim*(Tin-Tout)-StaleTepelneZisky)/1000;
FunkceZtrata[TeplotaOut_]:=(1/1000)(MernaZtrataProstupem*(Tin-TeplotaOut)+MernaZtrataVetranim*(Tin-TeplotaOut)-StaleTepelneZisky-MinSolar);
PouzeZtrata[TeplotaOut_]:=(MernaZtrataProstupem*(Tin-TeplotaOut)+MernaZtrataVetranim*(Tin-TeplotaOut));
GrafZtraty=Plot[FunkceZtrata[Tout],{Tout,-15,20}, PlotStyle->{Red,Thickness[ts]}]

Průběh venkovní teploty a délka topného období

PrubehTeploty={-1.0,1.0,4.0,9.0,14.6,17.0,18.2,18.8,13.8,9.4,4.0,-0.5,-1.0};
PrubehTeploty2={-1.0,1.0,4.0,9.0,14.6,17.0,18.2,18.8,13.8,9.4,4.1,-0.5};
Mesice={1,2,3,4,5,6,7,8,9,10,11,12,13};
DataTeploty=Riffle[Mesice,PrubehTeploty];
Temperature=Partition[DataTeploty,2];

Interpolace[mesic_]=Interpolation[Temperature][mesic];
Pomocna[m_]:=Interpolace[m]/;1≤m<13;
Pomocna[m_]:=Pomocna[m-12]/;m≥13;
Pomocna[m_]:=Pomocna[m+12]/;m<1;

JedenMesic=30.42*24*3600;
rok=365*24*60*60;

Mesic[t_]:=Pomocna[(t/JedenMesic)];
strhodn=(1/rok)NIntegrate[Mesic[t],{t,0,rok}];
SkutecnaStrhodn=9.07;
Korekce[t_]:=Mesic[t]-strhodn+SkutecnaStrhodn;
TopnaTeplota=13;
Plot[{Korekce[t],TopnaTeplota},{t,0,3*rok},AxesLabel->{"t (s)","(Subscript[θ, e])(°C)"},PlotStyle->{{Thickness[ts],Darker[Blue]},{Thickness[ts],Darker[Red]}},GridLines->Automatic]


ZacatekNetopne=t/.FindRoot[Korekce[t]==TopnaTeplota,{t,1*(10^7),2*(10^7)}];
ZacatekTopne=t/.FindRoot[Korekce[t]==TopnaTeplota,{t,2*(10^7),3*(10^7)}];
KonecTopne=t/.FindRoot[Korekce[t]==TopnaTeplota,{t,4*(10^7),5*(10^7)}];
ZacatekTopne2=t/.FindRoot[Korekce[t]==TopnaTeplota,{t,5*(10^7),6*(10^7)}];
DelkaTopne=KonecTopne-ZacatekTopne;
DelkaTopneDny=Round[(KonecTopne-ZacatekTopne)/(3600*24)];
Style["Délka topného období je " <>ToString[DelkaTopneDny]<> " dní.",FontSize->18,FontColor->Red, FontFamily->"Times"]

Off[InterpolatingFunction::"dmval",FindRoot::"cvmit",InverseFunction::"ifun",NIntegrate::"ncvb",NIntegrate::izero,NIntegrate::slwcon,ReplaceAll::reps];
Pocasi=WeatherData[{"Most", "Ustecky", "CzechRepublic"},"MeanTemperature", {{2012, 1, 1}, {2014,12,30}, "Day"}];
HodnotyTeplot=Pocasi["Values"];
HodnotybezJednotky=QuantityMagnitude[HodnotyTeplot,"Celsius"];
Dny=Table[i,{i,1,Length[HodnotyTeplot],1}];
Delka=Length[Dny];
TeplotaDny[i_]:={Dny[[i]],HodnotybezJednotky[[i]]};
DataTeplotaRok=Table[TeplotaDny[i],{i,Delka}];
InterpolaceTeplota[Den_]=Interpolation[DataTeplotaRok] [Den];
JedenDen=24*3600;
Den[t_]:=InterpolaceTeplota[(t/JedenDen)];
grafTeplota=Plot[{Den[t],TopnaTeplota},{t,0,3*rok},AxesLabel->{"t (s)","(Subscript[θ, e]) (°C)"},PlotStyle->{{Darker[Blue], Thickness[ts]},{Darker[Red], Thickness[ts]}}, GridLines->Automatic]

Potřeba tepla během topného období - norma

UcinnostSdileni=0.9;
UcinnostDistribuce=0.89;
UcinnostNadoba=0.98;
VyuzitiZisku={999,996,968,752,333,179,111,78,483,912,998,1000,999};
DataZisky=Riffle[Mesice,VyuzitiZisku];
Gains=Partition[DataZisky,2];

Interpolace2[mesic_]=Interpolation[Gains][mesic];
Pomocna2[m_]:=Interpolace2[m]/;1≤m<13;
Pomocna2[m_]:=Pomocna2[m-12]/;m≥13;
Pomocna2[m_]:=Pomocna2[m+12]/;m<1;

JedenMesic=30.42*24*3600;
rok=365*24*60*60;

Mesic2[t_]:=Pomocna2[(t/JedenMesic)];
strhodn2=(1/rok)NIntegrate[Mesic2[t],{t,0,rok}];
SkutecnaStrhodn2=648.6;
Korekce2[t_]:=Mesic2[t]-strhodn2+SkutecnaStrhodn2;
Plot[Korekce2[t]/1000,{t,0,2*rok},AxesLabel->{"t (s)"},PlotStyle->{Thickness[ts],Darker[Blue]}]

Zisky={(40227/(31*24*3600)),(43273/(28*24*3600)),(58838/(31*24*3600)),(67481/(30*24*3600)),(80701/(31*24*3600)),(80544/(30*24*3600)),(80282/(31*24*3600)),(75991/(31*24*3600)),(61728/(30*24*3600)),(48505/(31*24*3600)),(35363/(30*24*3600)),(34935/(31*24*3600)),(40227/(31*24*3600))}*3600000;
DataAbsolutni=Riffle[Mesice,Zisky];
GainsAll=Partition[DataAbsolutni,2];

Interpolace3[mesic_]=Interpolation[GainsAll][mesic];
Pomocna3[m_]:=Interpolace3[m]/;1≤m<13;
Pomocna3[m_]:=Pomocna3[m-12]/;m≥13;
Pomocna3[m_]:=Pomocna3[m+12]/;m<1;

JedenMesic=30.42*24*3600;
rok=365*24*60*60;

Mesic3[t_]:=Pomocna3[(t/JedenMesic)];
strhodn3=(1/rok)NIntegrate[Mesic3[t],{t,0,rok}];
SkutecnaStrhodn3=80807;
Korekce3[t_]:=Mesic3[t]-strhodn3+SkutecnaStrhodn3;
Plot[Korekce3[t],{t,0,2*rok},AxesLabel->{"t (s)"},PlotStyle->{Thickness[ts],Darker[Blue]}]

PotrebaTepla=NIntegrate[PouzeZtrata[Korekce[t]]-((Korekce2[t])/1000)*Korekce3[t],{t,ZacatekTopne,KonecTopne}]/((10^6)*3.6*UcinnostSdileni*UcinnostDistribuce*UcinnostNadoba);




Style["Potřeba tepla během jednoho topného období i po započtení účinností je: " <>ToString[PotrebaTepla]<> " kWh.",FontSize->18,FontColor->Red, FontFamily->"Times"]

(*Hodnota se mírně liší od vypočtené hodnoty v rámci technických systémů hodnocené budovy, protože se použila jiná účinnost zdroje - akumulační nádoby - dříve výměník. V rámci výpočotvého kroku jeden měsíc se také nedostanou úplně stejné hodnoty jako v rámci výpočtu za jedno topné období*)

Matematický popis topného okruhu s tepelným čerpadlem + stanovení bodu bivalence II. část

Off[InterpolatingFunction::"dmval",FindRoot::"cvmit",InverseFunction::"ifun",NIntegrate::"ncvb",NIntegrate::izero,NIntegrate::slwcon,ReplaceAll::reps];
Pk=8800*12;
c=4186.;
TeplotaOhrevTUV=50;
LitruDenne=8400;
PotrebaTeplaTUVdenne=TeplotaOhrevTUV*LitruDenne*c;
StalyVykonTUV=(PotrebaTeplaTUVdenne/(24*3600));
DelkaPotrubiTUV=650; (*menší než u původního provedení - nevede již do výměníku*)
MernaZtrataRozvod=185*3600;
ZtrataRozvod=(1/(24*3600))MernaZtrataRozvod*DelkaPotrubiTUV;
SpotrebaTeplaTUV=StalyVykonTUV+ZtrataRozvod;

FunkceZtrataUpravena[TeplotaOut_]:=(MernaZtrataProstupem*(Tin-TeplotaOut)+MernaZtrataVetranim*(Tin-TeplotaOut)-StaleTepelneZisky)+SpotrebaTeplaTUV;


ReseniTUV[Tvzduch_]:=Solve[SpotrebaTeplaTUV==StridaCerpadlo2*12000*Ptopna60[Tvzduch],StridaCerpadlo2][[1]];



rce1=m*c*(Tvoda-T2)==(StridaCerpadlo-StridaCerpadlo2/.ReseniTUV[Tvzduch])*(12*Ptopcompl[Tvzduch,Tvoda]);
rce2=m*c*(T1-Tvoda)==Pk*StridaKotel;
rce3=m*c(T1-T2)==1000*FunkceZtrata[Tvzduch];
rce4=kTop*(((T1+T2)/2)-Tint)==1000*FunkceZtrata[Tvzduch];

ParametrySoustavy={rce3,rce4};
ZnamyStav={Tvzduch->-15,Tint->20.,T1->60,T2->50};

StavSoustavy=Solve[ParametrySoustavy/.ZnamyStav][[1]];

VsechnyRovnice[Tvzduch_]={rce1,rce2,rce3,rce4}/.StavSoustavy;
Bivalence=FindRoot[VsechnyRovnice[Tvzduch]/.{StridaKotel->0,StridaCerpadlo->1,Tint->20},{{Tvzduch,-15,0},{T1,35.,60.},{T2,30,55.},{Tvoda,35.,60.}},MaxIterations->(10^2)];
BodBivalence =Tvzduch/.Bivalence;
Style["Teplota bivalence je: " <>ToString[Round[BodBivalence,0.01]]<> " °C",FontSize->18,FontColor->Red, FontFamily->"Times"]
GrafBivalence=Plot[{(1000*FunkceZtrata[Vzduch]+SpotrebaTeplaTUV),12*Ptopcompl[Vzduch,Tvoda]/.Bivalence},{Vzduch,-15,25},PlotStyle->{{Thickness[ts],Darker[Red]},{Thickness[ts],Darker[Green]}},PlotLegends->Placed[{"Potřebný tepelný výkon","Výkon tepelných čerpadel"},Below],AxesLabel->{"(Subscript[θ, e]) (°C)","P (W)"}]

(*Vypocet trva dele*)

Reseni[Tvzduch_]:=Union[FindRoot[VsechnyRovnice[Tvzduch]/.{StridaKotel->0,Tint->20},{{StridaCerpadlo,0,1},{T1,32.,58.},{T2,45.,55.},{Tvoda,32.,70.}},MaxIterations->(10^2)],{StridaKotel->0}]/;Tvzduch≥BodBivalence;
Reseni[Tvzduch_]:=Union[FindRoot[VsechnyRovnice[Tvzduch]/.{StridaCerpadlo->1,Tint->20},{{StridaKotel,0,1},{T1,32.,58.},{T2,45.,55.},{Tvoda,32.,55.}},MaxIterations->(10^2)],{StridaCerpadlo->1}]/;Tvzduch<BodBivalence;


(*Použití denních průměrných teplot*)
VykonTUV[t_]:=(StridaCerpadlo2/.ReseniTUV[Den[t]])*12000*Ptopna60[Den[t]];
CelkovyVykonTUV=NIntegrate[VykonTUV[t],{t,rok,2*rok}]/(1000000*3.6);

PrikonTUV[t_]:=(StridaCerpadlo2/.ReseniTUV[Den[t]])*12000*Pel60[Den[t]];
CelkovyPrikonTUV=NIntegrate[PrikonTUV[t],{t,rok,2*rok}]/(1000000*3.6);

TopnyFaktorTUV=CelkovyVykonTUV/CelkovyPrikonTUV;

TepelnyVykonTC[t_]:=(StridaCerpadlo-StridaCerpadlo2/.ReseniTUV[Den[t]])*12*Ptopcompl[Den[t],Tvoda]/.Reseni[Den[t]];
EnergieTC=NIntegrate[TepelnyVykonTC[t],{t, rok,KonecTopne}]/(1000000*3.6)+NIntegrate[TepelnyVykonTC[t],{t, ZacatekTopne2,2*rok}]/(1000000*3.6);

ElPrikonTC[t_]:=(StridaCerpadlo-StridaCerpadlo2/.ReseniTUV[Den[t]])*12*Pelcompl[Den[t],Tvoda]/.Reseni[Den[t]];
DodanaEnergieTC=NIntegrate[ElPrikonTC[t],{t,rok,KonecTopne}]/(1000000*3.6)+NIntegrate[ElPrikonTC[t],{t,ZacatekTopne2,2*rok}]/(1000000*3.6);

VykonElektroK[t_]:=StridaKotel*Pk/.Reseni[Den[t]];
PrikonElektroK[t_]:=StridaKotel*Pk/.Reseni[Den[t]];
DodanaEnergieEK=NIntegrate[VykonElektroK[t],{t,rok,KonecTopne}]/(1000000*3.6)+NIntegrate[VykonElektroK[t],{t,ZacatekTopne2,2*rok}]/(1000000*3.6);
CelkovyPrikonEK=NIntegrate[PrikonElektroK[t],{t,rok,KonecTopne}]/(1000000*3.6*0.95)+NIntegrate[PrikonElektroK[t],{t,ZacatekTopne2,2*rok}]/(1000000*3.6*0.95);

CelkovyPrikon=(DodanaEnergieTC+CelkovyPrikonEK);
TopnyFaktorTC=EnergieTC/DodanaEnergieTC;
CelkovyTepelnyVykonTC=EnergieTC+CelkovyVykonTUV;
CelkovyPrikonTC=DodanaEnergieTC+CelkovyPrikonTUV;
CelkovyCOP=CelkovyTepelnyVykonTC/CelkovyPrikonTC;

CelkovyVykonEK=DodanaEnergieEK;
SpotrebaZasobnikTUV=((2.3*365*10)/0.95);
VykonZasobnikTUV=2.3*365*10;
CelkovaSpotrebaElEnergie=SpotrebaZasobnikTUV+CelkovyPrikonTC+CelkovyPrikonEK;
BivalentniZdroje=SpotrebaZasobnikTUV+CelkovyPrikonEK;
CelkovaDodanaEnergie=CelkovyTepelnyVykonTC+CelkovyVykonEK+VykonZasobnikTUV;

(*Použití měsíčních průměrných teplot*)

VykonTUV2[t_]:=(StridaCerpadlo2/.ReseniTUV[Korekce[t]])*12000*Ptopna60[Korekce[t]];
CelkovyVykonTUV2=NIntegrate[VykonTUV2[t],{t,rok,2*rok}]/(1000000*3.6);

PrikonTUV2[t_]:=(StridaCerpadlo2/.ReseniTUV[Korekce[t]])*12000*Pel60[Korekce[t]];
CelkovyPrikonTUV2=NIntegrate[PrikonTUV2[t],{t,rok,2*rok}]/(1000000*3.6);

TopnyFaktorTUV2=CelkovyVykonTUV2/CelkovyPrikonTUV2;

TepelnyVykonTC2[t_]:=(StridaCerpadlo-StridaCerpadlo2/.ReseniTUV[Korekce[t]])*12*Ptopcompl[Korekce[t],Tvoda]/.Reseni[Korekce[t]];
EnergieTC2=NIntegrate[TepelnyVykonTC2[t],{t, rok,KonecTopne}]/(1000000*3.6)+NIntegrate[TepelnyVykonTC2[t],{t, ZacatekTopne2,2*rok}]/(1000000*3.6);

ElPrikonTC2[t_]:=(StridaCerpadlo-StridaCerpadlo2/.ReseniTUV[Korekce[t]])*12*Pelcompl[Korekce[t],Tvoda]/.Reseni[Korekce[t]];
DodanaEnergieTC2=NIntegrate[ElPrikonTC2[t],{t,rok,KonecTopne}]/(1000000*3.6)+NIntegrate[ElPrikonTC2[t],{t,ZacatekTopne2,2*rok}]/(1000000*3.6);

VykonElektroK2[t_]:=StridaKotel*Pk/.Reseni[Korekce[t]];
PrikonElektroK2[t_]:=StridaKotel*Pk/.Reseni[Korekce[t]];
DodanaEnergieEK2=NIntegrate[VykonElektroK2[t],{t,rok,KonecTopne}]/(1000000*3.6)+NIntegrate[VykonElektroK2[t],{t,ZacatekTopne2,2*rok}]/(1000000*3.6);
CelkovyPrikonEK2=NIntegrate[PrikonElektroK2[t],{t,rok,KonecTopne}]/(1000000*3.6*0.95)+NIntegrate[PrikonElektroK2[t],{t,ZacatekTopne2,2*rok}]/(1000000*3.6*0.95);

CelkovyPrikon2=(DodanaEnergieTC2+CelkovyPrikonEK2);
TopnyFaktorTC2=EnergieTC2/DodanaEnergieTC2;
CelkovyTepelnyVykonTC2=EnergieTC2+CelkovyVykonTUV2;
CelkovyPrikonTC2=DodanaEnergieTC2+CelkovyPrikonTUV2;
CelkovyCOP2=CelkovyTepelnyVykonTC2/CelkovyPrikonTC2;

CelkovyVykonEK2=DodanaEnergieEK2;
SpotrebaZasobnikTUV2=((2.3*365*10)/0.95);
VykonZasobnikTUV2=2.3*365*10;
CelkovaSpotrebaElEnergie2=SpotrebaZasobnikTUV+CelkovyPrikonTC2+CelkovyPrikonEK2;
BivalentniZdroje2=SpotrebaZasobnikTUV+CelkovyPrikonEK2;
CelkovaDodanaEnergie2=CelkovyTepelnyVykonTC2+CelkovyVykonEK2+VykonZasobnikTUV2;

Style["Dodané teplo tepelnými čerpadly během topné sezony na vytápění: " <>ToString[EnergieTC]<> " kWh.",FontSize->18,FontColor->Red, FontFamily->"Times"]
Style["Dodané teplo tepelnými čerpadly během topné sezony na vytápění: " <>ToString[EnergieTC2]<> " kWh.",FontSize->18,FontColor->Green, FontFamily->"Times"]
Style["Spotřeba elektrické energie během topné sezony tepelnými čerpadly: " <>ToString[DodanaEnergieTC]<> " kWh.",FontSize->18,FontColor->Red, FontFamily->"Times"]
Style["Spotřeba elektrické energie během topné sezony tepelnými čerpadly: " <>ToString[DodanaEnergieTC2]<> " kWh.",FontSize->18,FontColor->Green, FontFamily->"Times"]
Style["Dodané teplo tepelnými čerpadly na ohřev TUV během jednoho roku: " <>ToString[CelkovyVykonTUV]<> " kWh.",FontSize->18,FontColor->Red, FontFamily->"Times"]
Style["Dodané teplo tepelnými čerpadly na ohřev TUV během jednoho roku: " <>ToString[CelkovyVykonTUV2]<> " kWh.",FontSize->18,FontColor->Green, FontFamily->"Times"]
Style["Spotřeba elektrické energie na ohřev TUV teplenými čerpadly za rok: " <>ToString[CelkovyPrikonTUV]<> " kWh.",FontSize->18,FontColor->Red, FontFamily->"Times"]
Style["Spotřeba elektrické energie na ohřev TUV teplenými čerpadly za rok: " <>ToString[CelkovyPrikonTUV2]<> " kWh.",FontSize->18,FontColor->Green, FontFamily->"Times"]
Style["Spotřeba elektrické energie zásobníky TUV: " <>ToString[SpotrebaZasobnikTUV]<> " kWh.",FontSize->18,FontColor->Red, FontFamily->"Times"]
Style["Spotřeba elektrické energie zásobníky TUV: " <>ToString[SpotrebaZasobnikTUV2]<> " kWh.",FontSize->18,FontColor->Green, FontFamily->"Times"]
Style["Celkové dodané teplo tepelnými čerpadly za rok: " <>ToString[CelkovyTepelnyVykonTC]<> " kWh.",FontSize->18,FontColor->Red, FontFamily->"Times"]
Style["Celkové dodané teplo tepelnými čerpadly za rok: " <>ToString[CelkovyTepelnyVykonTC2]<> " kWh.",FontSize->18,FontColor->Green, FontFamily->"Times"]
Style["Celková spotřeba elektrické energie tepelnými čerpadly za rok: " <>ToString[CelkovyPrikonTC]<> " kWh.",FontSize->18,FontColor->Red, FontFamily->"Times"]
Style["Celková spotřeba elektrické energie tepelnými čerpadly za rok: " <>ToString[CelkovyPrikonTC2]<> " kWh.",FontSize->18,FontColor->Green, FontFamily->"Times"]
Style["Celková spotřeba elektrické energie systémem za rok: " <>ToString[CelkovaSpotrebaElEnergie]<> " kWh.",FontSize->18,FontColor->Red, FontFamily->"Times"]
Style["Celková spotřeba elektrické energie systémem za rok: " <>ToString[CelkovaSpotrebaElEnergie2]<> " kWh.",FontSize->18,FontColor->Green, FontFamily->"Times"]
Style["Celková spotřeba elektrické energie bivalentních zdrojů za rok: " <>ToString[BivalentniZdroje]<> " kWh.",FontSize->18,FontColor->Red, FontFamily->"Times"]
Style["Celková spotřeba elektrické energie bivalentních zdrojů za rok: " <>ToString[BivalentniZdroje2]<> " kWh.",FontSize->18,FontColor->Green, FontFamily->"Times"]
Style["Průměrný topný faktor za rok: " <>ToString[CelkovyCOP]<> ".",FontSize->18,FontColor->Red, FontFamily->"Times"]
Style["Průměrný topný faktor za rok: " <>ToString[CelkovyCOP2]<> ".",FontSize->18,FontColor->Green, FontFamily->"Times"]
Style["Celkové teplo dodané elektrokotli: " <>ToString[CelkovyVykonEK]<> " kWh.",FontSize->18,FontColor->Red, FontFamily->"Times"]
Style["Celkové teplo dodané elektrokotli: " <>ToString[CelkovyVykonEK2]<> " kWh.",FontSize->18,FontColor->Green, FontFamily->"Times"]
Style["Celková spotřeba elektrické energie elektrokotli: " <>ToString[CelkovyPrikonEK]<> " kWh.",FontSize->18,FontColor->Red, FontFamily->"Times"]
Style["Celková spotřeba elektrické energie elektrokotli: " <>ToString[CelkovyPrikonEK2]<> " kWh.",FontSize->18,FontColor->Green, FontFamily->"Times"]
Style["Celkové dodané teplo celou soustavou: " <>ToString[CelkovaDodanaEnergie]<> " kWh.",FontSize->18,FontColor->Red, FontFamily->"Times"]
Style["Celkové dodané teplo celou soustavou: " <>ToString[CelkovaDodanaEnergie2]<> " kWh.",FontSize->18,FontColor->Green, FontFamily->"Times"]