Sub spr()

Call obieg(i_1)

End Sub


Sub obieg(i_1)

Worksheets("Arkusz1").Activate


'deklaracje entalpie
'i_1 As Double,

Dim i_2 As Double, i_2s As Double, i_3 As Double, i_4 As Double, i_5 As Double, _
 i_5s As Double, i_6 As Double, i_7 As Double, i_8 As Double, i_9 As Double, i_10 As Double, _
 i_11 As Double, i_12 As Double, i_13 As Double, i_14 As Double, i_15 As Double, i_16 As Double, _
 i_17 As Double, i_18 As Double, i_19 As Double, i_20 As Double, i_21 As Double, i_22 As Double, _
 i_23 As Double, i_24 As Double, i_25 As Double


' deklaracja entropie

Dim s_1 As Double, s_2 As Double, s_3 As Double, s_4 As Double, s_5 As Double, _
 s_10 As Double, s_11 As Double, s_12 As Double, s_14 As Double, s_15 As Double, s_16 As Double, _
 s_20 As Double, s_22 As Double, s_24 As Double


' deklaracja egzergii
Dim b_1 As Double, b_2 As Double, b_3 As Double, b_4 As Double, b_5 As Double, b_6 As Double, _
b_10 As Double, b_11 As Double, b_12 As Double, b_14 As Double, b_15 As Double, b_16 As Double, _
b_17 As Double, b_19 As Double, b_20 As Double, b_22 As Double, b_24 As Double



'Wczytanie danych z zadania
Q_g_0 = Range("C3")
eta_HT_CHP = Range("C4")
Q_t_0 = Range("C5")
p_1 = Range("C6")
t_1 = Range("C7")
p_2 = Range("C8")
p_5 = Range("C9")
eta_IHT = Range("C10")
eta_ILT = Range("C11")
N_P1 = Range("C12")
N_P2 = Range("C13")
t_15 = Range("C14")
t_14 = Range("C15")
t_10 = Range("C16")
t_11 = Range("C17")
t_12 = Range("C18")
t_13 = Range("C19")
t_17 = Range("C20")
t_18 = Range("C21")
t_19 = Range("C22")
t_20 = Range("C23")
t_24 = Range("C24")
t_16 = Range("C25")
t_22 = Range("C26")
T_sur = Range("C27")
t_n = Range("C28")
c_w = Range("c29")

'dane temperatury
T_sur = 268     'K
t_n = 273.15    'K
c_w = 4.19      'kJ/kgK

'***** OBLICZENIA *****

'Turbina pkt 1
 i_1 = h_pT(p_1, t_1)
 s_1 = s_pT(p_1, t_1)

' Turbina za częścią HT pkt 2s
i_2s = h_ps(p_2, s_1)
i_2 = i_1 - eta_IHT * (i_1 - i_2s)
s_2 = s_ph(p_2, i_2)
i_5s = h_ps(p_5, s_2)
i_5 = i_2 - eta_ILT * (i_2 - i_5s)
s_5 = s_ph(p_5, i_5)

i_3 = i_2
i_4 = i_2
s_3 = s_2
s_4 = s_2

i_6 = i_5
i_7 = i_5
s_6 = s_5
s_7 = s_5

'Entalpie dla wody w charakterystycznych punktach obiegu
i_11 = hL_T(t_11)
i_12 = hL_T(t_12)
i_17 = hL_T(t_17)
i_10 = hL_T(t_10)
i_16 = hL_T(t_16)
i_14 = hL_T(t_14)
i_15 = hL_T(t_15)
i_22 = hL_T(t_22)
i_19 = hL_T(t_19)
i_20 = hL_T(t_20)
i_24 = hL_T(t_24)



'Strumień ciepła dostarczany do USER G
Q_g = Q_g_0 / eta_HT_CHP


'Strumień ciepła dostarczany do USER T
Q_t = Q_t_0 / eta_HT_CHP


'Bilans energii dla USER G przekształcony w celu wyznaczenia strumienia m_11
'm_11=m_12
'uproszczenie te stosujemy do wyznaczenia strumienia m_11 z bilansu energii dla USER G
m_11 = Q_g / (i_11 - i_12)

'zachodzi równość
m_8 = m_11
m_9 = m_11
m_10 = m_11
m_12 = m_11

'w treści zadania
m_13 = 0.1 * m_11
m_15 = m_13


'Bilans substancji - Separator
m_14 = m_12 - m_13

'Bilans energii dla USER T przekształcony w celu wyznaczenia strumienia m_3
m_3 = Q_t / (i_3 - i_17)

'Bilans substancji - USER T
m_17 = m_3

'z treści zadania
m_18 = 0.75 * m_3
m_20 = m_18


'Bilans dla wymiennika ciepła - HE 1
'm_4=m_16
'uproszczenie te stosujemy do wyznaczenia strumienia m_4 z bilansu energii dla tego wymiennika
m_4 = (m_11 * i_11 - m_10 * i_10) / (i_4 - i_16)
'Bilans substancji
m_16 = m_4
m_2 = m_3 + m_4

'Bilans substancji - MIXER 1
m_8 = m_15 + m_14
'Bilans energii dla MIXER 1 przekształcony w celu wyznaczenia entalpii i_8
i_8 = (m_15 * i_15 + m_14 * i_14) / m_8


'Bilans substancji dla pompy - PUMP 1
'm_8 = m_9

'Bilans energii dla pompy przekształcony w celu wyznaczenia entalpii i_9
i_9 = (m_8 * i_8 + N_P1) / m_9


'Bilans dla wymiennika ciepła - HE2
'Przekształcony bilans energii z zastosowaniem podstawienia m_22=m_6+m_16
m_6 = (m_16 * i_22 + m_10 * i_10 - m_9 * i_9 - m_16 * i_16) / (i_6 - i_22)

m_22 = m_6 + m_16

'Bilans oddzielacza SEPARATOR 2
'Bilans substancji
m_18 = m_20
m_3 = m_17
m_19 = m_17 - m_18

'Bilans mieszalnika MIXER 2
'Bilans substancji
m_21 = m_3
'Bilans energii mieszalnika przekształcony w celu wyznaczenia entalpii i_21
i_21 = (m_20 * i_20 + m_19 * i_19) / m_21


'Bilans mieszalnika MIXER 3
m_23 = m_21 + m_22
'Bilans energii mieszalnika przekształcony w celu wyznaczenia entalpii i_23
i_23 = (m_22 * i_22 + m_21 * i_21) / m_23


'Bilans odgazowywacza - DEAERATOR
'm_7=m_24-m_23
'Bilans energii odgazowywacza przekształcony w celu wyznaczenia strumienia m_24
m_24 = (m_23 * (i_23 - i_7)) / (i_24 - i_7)

'Z bilansu substancji
m_7 = m_24 - m_23


'Bilans substancji dla pompy - PUMP2
m_1 = m_24
m_25 = m_24
'Bilans energii pompy przekształcony w celu wyznaczenia entalpii i_25
i_25 = (m_24 * i_24 + N_P2) / m_25

'Bilans substancji turbiny
m_1 = m_25
m_1 = m_24
m_5 = m_6 + m_7
m_2 = m_1 - m_5
'Bilans energii dla turbiny
N_IHT = m_1 * (i_1 - i_2) 'kW
N_ILT = m_5 * (i_2 - i_5) 'kW
N_IT = N_IHT + N_ILT      'kW





'BILANS EGZERGII

'TURBINA
b_1 = i_1 - T_sur * s_1     'kJ/kg
b_2 = i_2 - T_sur * s_2     'kJ/kg
b_5 = i_5 - T_sur * s_5     'kJ/kg

delta_B_HT = m_1 * (b_1 - b_2) - N_IHT  'kW
delta_B_LT = m_5 * (b_2 - b_5) - N_ILT  'kW
delta_B_Turbine = delta_B_HT + delta_B_LT     'kW


'USER G
s_11 = c_w * (Log((t_11 + 273.15) / t_n))
s_12 = c_w * (Log((t_12 + 273.15) / t_n))

b_11 = i_11 - T_sur * s_11
b_12 = i_12 - T_sur * s_12
delta_B_G = m_11 * (b_11 - b_12)        'kW


'USER T
s_3 = s_2
s_17 = c_w * (Log((t_17 + 273.15) / t_n))
b_3 = i_3 - T_sur * s_3
b_2 = b_3
b_17 = i_17 - T_sur * s_17
delta_B_T = m_3 * (b_3 - b_17)          'kW


'HEAT EXCHANGER 1
b_4 = b_2
s_10 = c_w * (Log((t_10 + 273.15) / t_n))
s_11 = c_w * (Log((t_11 + 273.15) / t_n))
s_16 = c_w * (Log((t_16 + 273.15) / t_n))

b_10 = i_10 - T_sur * s_10
b_11 = i_11 - T_sur * s_11
b_16 = i_16 - T_sur * s_16

delta_B_HE_1 = m_4 * (b_4 - b_16) + m_10 * (b_10 - b_11)    'kW


'HEAT EXCHANGER 2
b_6 = b_5
s_14 = c_w * (Log((t_14 + 273.15) / t_n))
s_15 = c_w * (Log((t_15 + 273.15) / t_n))
s_22 = c_w * (Log((t_22 + 273.15) / t_n))

b_14 = i_14 - T_sur * s_14
b_15 = i_15 - T_sur * s_15
b_22 = i_22 - T_sur * s_22

delta_B_HE_2 = m_6 * b_6 + m_16 * b_16 + m_15 * b_15 + m_14 * b_14 - m_10 * b_10 - m_22 * b_22


'DEAERATOR
b_7 = b_5
s_19 = c_w * (Log((t_19 + 273.15) / t_n))
s_20 = c_w * (Log((t_20 + 273.15) / t_n))
s_24 = c_w * (Log((t_24 + 273.15) / t_n))

b_19 = i_19 - T_sur * s_19
b_20 = i_20 - T_sur * s_20
b_24 = i_24 - T_sur * s_24

delta_B_D = m_7 * b_7 + m_20 * b_20 + m_19 * b_19 + m_22 * b_22 - m_24 * b_24


'Zczytywanie wartości do arkusza
'zczytywanie entalpii
Range("G4") = i_1
Range("G5") = i_2
Range("G6") = i_2s
Range("G7") = i_3
Range("G8") = i_4
Range("G9") = i_5
Range("G10") = i_5s
Range("G11") = i_6
Range("G12") = i_7
Range("G13") = i_8
Range("G14") = i_9
Range("G15") = i_10
Range("G16") = i_11
Range("G17") = i_12
Range("G18") = " "
Range("G19") = i_14
Range("G20") = i_15
Range("G21") = i_16
Range("G22") = i_17
Range("G23") = " "
Range("G24") = i_19
Range("G25") = i_20
Range("G26") = i_21
Range("G27") = i_22
Range("G28") = i_23
Range("G29") = i_24
Range("G30") = i_25


'zczytywanie entropii
Range("H4") = s_1
Range("H5") = s_2
Range("H6") = s_2s
Range("H7") = s_3
Range("H8") = s_4
Range("H9") = s_5
Range("H10") = s_5s
Range("H11") = s_6
Range("H12") = s_7
Range("H13") = s_8
Range("H14") = s_9
Range("H15") = s_10
Range("H16") = s_11
Range("H17") = s_12
Range("H18") = s_13
Range("H19") = s_14
Range("H20") = s_15
Range("H21") = s_16
Range("H22") = s_17
Range("H23") = s_18
Range("H24") = s_19
Range("H25") = s_20
Range("H26") = s_21
Range("H27") = s_22
Range("H28") = s_23
Range("H29") = s_24
Range("H30") = s_25


'zczytywanie strumieni
Range("I4") = m_1
Range("I5") = m_2
Range("I6") = " "
Range("I7") = m_3
Range("I8") = m_4
Range("I9") = m_5
Range("I10") = " "
Range("I11") = m_6
Range("I12") = m_7
Range("I13") = m_8
Range("I14") = m_9
Range("I15") = m_10
Range("I16") = m_11
Range("I17") = m_12
Range("I18") = m_13
Range("I19") = m_14
Range("I20") = m_15
Range("I21") = m_16
Range("I22") = m_17
Range("I23") = m_18
Range("I24") = m_19
Range("I25") = m_20
Range("I26") = m_21
Range("I27") = m_22
Range("I28") = m_23
Range("I29") = m_24
Range("I30") = m_25

'zczytywanie egzergii
Range("J4") = b_1
Range("J5") = b_2
Range("J6") = " "
Range("J7") = b_3
Range("J8") = b_4
Range("J9") = b_5
Range("J10") = " "
Range("J11") = b_6
Range("J12") = b_7
Range("J13") = b_8
Range("J14") = b_9
Range("J15") = b_10
Range("J16") = b_11
Range("J17") = b_12
Range("J18") = b_13
Range("J19") = b_14
Range("J20") = b_15
Range("J21") = b_16
Range("J22") = b_17
Range("J23") = b_18
Range("J24") = b_19
Range("J25") = b_20
Range("J26") = b_21
Range("J27") = b_22
Range("J28") = b_23
Range("J29") = b_24
Range("J30") = b_25

'pozostałe
Range("N4") = Q_g
Range("N5") = Q_t
Range("N6") = N_IHT
Range("N7") = N_ILT
Range("N8") = N_IT
Range("N9") = delta_B_HT
Range("N10") = delta_B_LT
Range("N11") = delta_B_Turbine
Range("N12") = delta_B_G
Range("N13") = delta_B_T
Range("N14") = delta_B_HE_1
Range("N15") = delta_B_HE_2
Range("N16") = delta_B_D

MsgBox "Yes."

End Sub