Quiet@Remove["Global`*"];
$HistoryLength=2;
SetDirectory[NotebookDirectory[]]

(***********************************************************************************)
(* Výpočet tepelných ztrát objektu*)
(***********************************************************************************)
(*Definovani vypoctovych teplot*)
Tout1=-12;(* venkovní výpočtová hodnota*)
Tout2=-18;(*venkovni výpočtová hodnota 2*)
T2=5.2;(* prumerna rocni teplota*)
Tin=21;(*vnitrni navrhova teplota*)
Tz=5;(*teplota země*)
deltaT=Tin-Tout1;(*rozdíl vnitřní a venkovní teploty*)
deltaT2=Tin-Tout2;(*rozdíl vnitřní a venkovní teploty*)
k=1000;
mm=1/1000;


(*Zakladni informace o objektu
 celkova podlahova plocha vyhrivaneho prostoru *)
plocha=8.480*(1.710+0.900+0.620+0.900+1.710)+3.380*(1.710+9.00+6.20+9.00+1.710);
kontrola = 6.20+22.10+1.30+5.40+1+1.10+16.2;
(*celkova plocha
izolace 
podlaha - v podlaze 8 cm mineralni vlny, pod zakladovou desku polozena tepelna izolace z podlahoveho polystyrenu o tl 100 mm, ve stropnim kci podkrovi - 300 mm mineralni vaty,
v obvodovych stenach - 250 mm ¨mineralni vaty.
celkova plocha obalky budovy 316.9 m2
*)
(***********************************************************************************)
(*součet odpor při při přestupu tepla na vnitřní a vnější straně stěn - Rsi + Rse*)
Rsi=0.04;
RseDown=0.17;
RseUp=0.10;
RseHorizont=0.13;

(***********************************************************************************)
(*zadani jednotlivych ploch vnejsi konstrukce - steny, strop a podlaha*)
plochaPodlaha=52.6;
plochaPodlahaZasklepena=16.6;
plochaSvisleSteny=147.1651;
strechaSikma=43.8;
strechaVodorovna=38.5;

(***********************************************************************************)
(*Materialy sten*)
(*fermacel tl. 12.5 mm*)
deltaF=12.5mm;
lambdaF=0.32;(*W/m*K*)
tkapacitaF=1.1k;(*J/kg*K*)
hustF=1150; (*+-50 kg/m3, jedna se o objemovou hmotnost*)
(*OSB deska tl. 18 mm*)
deltaO=18mm;
lambdaO=0.13;(* W/m*K *)
hustO=600;(*kg/m3, objemova hustota*)
tkapacitaO;
(*OSB deska 22mm*)
deltaOt=22mm;
(*deskový podlahový polystyren 100mm, Polystyren EPS 100 S Stabil*)
deltaP=100mm;
lambdaP=0.038;
(*kontaktni zateplovaci system (vnejsi strana)*)
lambdaK=0.04;(*FKL, Knaufinsulation*)
deltaK=30mm;
(*mineralni vlna / knauf classic 039*)
(*součinitel tepelné vodivosti izolace je 0.039, součinitel upraven pro respektování podílu dřeva ve stěnách (10%)*)
lambdaV=0.0571;
deltaVstr=300mm;(*ve stropu*)
deltaVp=100mm;(*v podlaze*)
deltaVste=250mm;(*ve svislych stenach*)

(***********************************************************************************)

(*Součinitel prostupu tepla svislé stěny*)
Rsvisla=deltaF/lambdaF+deltaO/lambdaO+deltaK/lambdaK+deltaVste/lambdaV;
Usvisla=1/(Rsi+    RseHorizont+Rsvisla);
(*Součinitel prostupu tepla vodorovný strop*)
RstropHorizontal=deltaF/lambdaF+deltaOt/lambdaO+deltaVstr/lambdaV;
UstropHorizontal=1/(Rsi+RseUp+RstropHorizontal);
(*Součinitel prostupu tepla šikmý strop*)
RstropSikmy=deltaF/lambdaF+deltaO/lambdaO+deltaVstr/lambdaV;
UstropSikmy=1/(Rsi+RseUp+RstropSikmy);
(*Součinitel prostupu tepla podlaha*)
Rpodlaha=deltaOt/lambdaO+deltaVp/lambdaV+deltaP/lambdaP;
Upodlaha=1/(Rsi+RseDown+Rpodlaha);

(***********************************************************************************)

(*okna - plastova okna zasklena izolacnim trojsklem*)
Uw=0.8; (* W/m2*K - trojsklo*)
Ur=1;     (* W/m2*K - ram *)
Uwmale=1.1; (* W/m2*K*)
Uwstresni=1.1;(* W/m2*K - Střešní okno Velux GLL 1061 *)
pocetOken=11;(*pocet beznych oken 900x1200*)
pocetOs=1;(*pocet oken stresnich*)
pocetOm=1;(*pocet oken malych*)
plochaOken=1.2*0.9;(*plocha beznych oken*)
plochaOs=0.87*0.98;(*plocha stresniho okna*)
plochaOm=0.9*0.5;(*plocha maleho okna*)
plochaBalkonoveDvere=3.4;
UBalkonoveDvere=1.10;
UvchodoveDveře=1.2;
plochaVchodoveDvere=2.1;

(***********************************************************************************)
(*Výpočet tepelných ztrát prostupem*)
Qp1=Usvisla*plochaSvisleSteny*deltaT+UstropHorizontal*strechaVodorovna*deltaT+UstropSikmy*strechaSikma*deltaT+Upodlaha*plochaPodlaha*(Tin-Tz)+(1/(Rsi+RseHorizont+1/Uw))*deltaT*pocetOken*plochaOken*Ur+plochaOs*(1/(Rsi+RseHorizont+1/Uwstresni))*pocetOs*Ur*deltaT+plochaOm*(1/(Rsi+RseHorizont+1/Uwmale))*deltaT+plochaBalkonoveDvere*(1/(Rsi+RseHorizont+1/UBalkonoveDvere))*deltaT+plochaVchodoveDvere*(1/(1/UvchodoveDveře+Rsi+RseHorizont))*deltaT;

Qp2=Usvisla*plochaSvisleSteny*deltaT2+UstropHorizontal*strechaVodorovna*deltaT2+UstropSikmy*strechaSikma*deltaT2+Upodlaha*plochaPodlaha*(Tin-Tz)+(1/(Rsi+RseHorizont+1/Uw))*deltaT2*pocetOken*plochaOken*Ur+plochaOs*(1/(Rsi+RseHorizont+1/Uwstresni))*pocetOs*Ur*deltaT2+plochaOm*(1/(Rsi+RseHorizont+1/Uwmale))*deltaT2+plochaBalkonoveDvere*(1/(Rsi+RseHorizont+1/UBalkonoveDvere))*deltaT2+plochaVchodoveDvere*UvchodoveDveře*deltaT2;
(***********************************************************************************)
(*výpočet ztrát větráním*)
objemBudovy=381;
ztratyVetranim1=1300*0.0277778*deltaT/10(*W*)(*0.0277778*)

ztratyVetranim2=1300*0.0277778*deltaT2/10(*W*);
0.0125;
0.015278;

(***********************************************************************************)

CelkoveTepelneZtraty1=ztratyVetranim1+Qp1(*W*)
CelkoveTepelneZtraty2=ztratyVetranim2+Qp2(*W*);

(***********************************************************************************)
teplozaden=CelkoveTepelneZtraty2*60*60*24/(10^6)/3.6;

(***********************************************************************************)
(*Stanovení tepelných zisků zasklení ze slunečního záření*)
(*Ezm=Egm*plochaPrusvitneCasti*celkovaPropustnost*cinitelVyuziti*cinitel dopadu*)
plochaPrusvitneCastiJih=0.88*1.180*5;
plochaPrusvitneCastiVychod=0.88*1.180*4+0.77*0.97+1.8+0.6*0.78;(*do plochy okne orientovanych na vychod byla zapocita i plocha oken orientovanych na zapad*)
celkovaPropustnost=0.73*0.42*0.9;(*jedná se o trojskla, se světlou reflexní fólií, 0.9 repektuje znečištění skla*)
AVcinitelVyuziti1=0.91;(*pro východ a západ*)
AVcinitelVyuziti2=0.80; (*pro jih*)
ciniteldopadu=0.9;
EgmAVv= 211.23; (*průměrná hodnota globálního sluneční záření za topné období pro východ a západ, topné období je od října do dubna*)
EgmAVj= 416.99; (*průměrná hodnota globálního sluneční záření za topné období pro jih*)
EzmJih=EgmAVj*plochaPrusvitneCastiJih*celkovaPropustnost*AVcinitelVyuziti2*ciniteldopadu;(*kWh/topné období*)
EzmVychod=EgmAVv*plochaPrusvitneCastiVychod*celkovaPropustnost*AVcinitelVyuziti2*ciniteldopadu;(*výpočet pro západ i východ*)(*kWh/topné období*);
SlunecniZiskCelk=EzmJih+EzmVychod;
SlunecniZiskAV=SlunecniZiskCelk/211(*Průměrný denní sluneční zisk - kWh*);

Ptopny=Fit[{{-12,CelkoveTepelneZtraty1},{-18,CelkoveTepelneZtraty2}},{1,T},T]

Ptopny>>"Ptopny"