Trnsys: Optimierung der Prozesse

Diese Seite verwendet Cookies. Durch die Nutzung unserer Seite erklären Sie sich damit einverstanden, dass wir Cookies setzen. Weitere Informationen

  • Trnsys: Optimierung der Prozesse

    Hallo,

    wir sind gerade dabei in der Universität mit Trnsys Simulationen zu gewissen realen Anwendungen zu erstellen. In meinem Fall (dient nur als Beispiel) brauche ich Hilfe bei einer, eventuell simplen, Frage: gibt es eine Möglichkeit Auswahlmöglichkeiten (also Schaltflächen, Konsoleneingaben) in der Simulation, also bei dem Graphen, zu erstellen, womit ich gewisse Typen, respektive Berechnungen, anhand von verschiedenen, von mir definierten, Fällen durchführen kann? In meinem Fall würde ich gerne je nach Benutzereingabe (in der Simulation) verschiedene Berechnungen für verschiedene Wärmepumpen durchführen können. Simpel gefragt: kann man in Trnsys selber Schaltflächen erstellen, also den Nutzer mit der Simulation interagieren lassen?


    Brainfuck-Quellcode

    1. Version 15.0
    2. *********************************************************************
    3. *
    4. * Brauchwasseranlage für IT-Projekt
    5. *
    6. * Version 1.1
    7. *
    8. *********************************************************************
    9. *
    10. * (basiert auf IEA High Performance House
    11. * Row_21_T_E_C_HP.DCK)
    12. *
    13. ASSIGN DHW_Listing.LST 6
    14. ASSIGN DHW_out.hour 93
    15. ASSIGN DHW_out.year 94
    16. ASSIGN DHW_out.month 95
    17. *
    18. * Brauchwasserprofil (Stundenweise)
    19. ASSIGN ..\Daten\brau1001.dat 22
    20. *
    21. * !!!!! Anpassen von TGROUND bei Wechsel der Wetterdaten!!!!!
    22. ASSIGN ..\Daten\MilaHour.dat 21
    23. *
    24. CONSTANCE 1
    25. * Mittlere Erdtemperatur
    26. TGROUND = 6.69
    27. *
    28. EQUATIONS 3
    29. * Simulationszeit: 1 Jahr (8760 h) --- ACHTUNG: START bei 1 (nicht 0)
    30. START = 1
    31. STOP = 8760
    32. STEP = 0.25
    33. *
    34. SIMULATION START STOP STEP
    35. WIDTH 72
    36. LIMITS 500 200
    37. *
    38. UNIT 9 TYPE 9 WETTERDATEN LESEN
    39. PARAMETERS 30
    40. * MODE SKIP NUMBER INTERVAL
    41. -1 0 8 1
    42. * Daten aus Datei -- Bei Minus keine Interpolation
    43. -1 1 0
    44. -2 1 0
    45. -3 1 0
    46. -4 3.6 0
    47. -5 3.6 0
    48. -6 1 0
    49. 7 1 0
    50. 8 1 0
    51. 21 -1
    52. *
    53. UNIT 10 TYPE 9 Einlesen Lastprofil Brauchwasser
    54. par 9
    55. -1 0 1 1
    56. -1 1 0
    57. 22 -1
    58. *
    59. EQUATIONS 1
    60. TA=[9,6]
    61. *
    62. *
    63. *********************************************************************
    64. * HEATING SYSTEM AND STORAGE TANK
    65. **********************************************************************
    66. *
    67. CONSTANTS 1
    68. * U-Wert Tank
    69. UTANK = 1.0
    70. *
    71. EQUATIONS 2
    72. * DHW load:
    73. * brau1001.dat is normalised to 100 l/d
    74. * 4 Persons x 40 l/d = 160 l/d at 40°C
    75. DHWload = 1.6*[10,1]*(40-TGROUND)/(MAX(30,[22,3])-TGROUND)
    76. * Power for DHW (in kJ/h)
    77. DHWpow = 3600*6
    78. *
    79. *---------------------------------------------------------------------
    80. *
    81. UNIT 22 TYPE 4 Storage Tank
    82. PARAMS 20
    83. * 1 2 3 4 5 6 7 8 9 10 11 12
    84. * mode Vol C2 rho Ut H1 AuxM heating1............
    85. 1 0.3 4.19 1000 UTank -1.6 1 1 1 51 2 DHWpow
    86. *
    87. * 13 14 15 16 17 18 19 20
    88. * heating2............... UAf Tflue Tboil
    89. 2 2 51 2 0 0 15 105
    90. * UAf=9.3 to model 60 W circulation loss
    92. INPUTS 7
    93. * T_hS M_hS T_load M_load T_env heat1 heat2
    94. 0,0 0,0 TGROUND DHWload 0,0 0,0 0,0
    95. 15 0 8.0 0 15 1 0
    96. DERIVATIVES 3
    97. * 1 Initial temp. of fluid: top tank segment ... bottom
    98. 50 30 10
    99. *
    100. *
    101. *
    102. *******************************************************
    103. * OUTPUT
    104. *******************************************************
    105. *
    106. EQUATIONS 1
    107. QINh = [22,8]/3600
    109. UNIT 40 TYPE 25 PRINTER hour
    110. PARAMETERS 4
    111. 1 START STOP 93
    112. INPUTS 6
    113. TA DHWload 22,3 22,13 22,1 QINh
    114. TOUT DHWload T1 T2 T3 QIN
    115. *
    116. *
    117. UNIT 41 TYPE 24 INTEGRATOR MONTH
    118. PARAMETERS 1
    119. -1
    120. INPUTS 3
    121. 22,5 22,6 22,8
    122. 0 0 0
    123. *
    124. EQUATIONS 4
    125. QLOSS = [41,1]/3600
    126. QLOAD = [41,2]/3600
    127. QIN = [41,3]/3600
    128. QOUT = QLOAD+QLOSS
    129. *
    130. UNIT 42 TYPE 25 PRINTER MONTH
    131. PARAMETERS 4
    132. -1 START STOP 95
    133. INPUTS 4
    134. QLOAD QLOSS QIN QOUT
    135. QLOAD QLOSS QIN QOUT
    136. *
    137. CONSTANCE 1
    138. *Primärenergiefaktor Strom
    139. PEF_EL = 2.35
    140. *
    141. EQUATIONS 1
    142. TEMP = -7 + (27/8760)*TIME
    144. UNIT 90 TYPE 113
    145. PARAMETERS 3
    146. 1 1 1
    147. INPUTS 2
    148. 0,0 TEMP
    149. 35 -7
    151. * 22,3 TA
    152. * 50 10
    154. EQUATIONS 4
    155. *
    156. * Hier kann der selbst errechnete COP eingebungen werden
    157. MY_COP=[90,1]
    158. * Energie für Brauchwasser
    159. DHWaux=[22,8]/3600
    160. * Primärenergie für Brauchwasser
    161. PRI_DHW=DHWaux*PEF_EL/MAX(1,MY_COP)
    162. * Gesamte Primärenergie mit Berücksichtigung von 5 kWh für Hilfsenergie
    163. PRIMARY=PRI_DHW+PEF_EL*5/8760
    165. *
    166. *
    167. UNIT 43 TYPE 24 INTEGRATOR YEAR
    168. PARAMETERS 1
    169. 8760
    170. INPUT 3
    171. DHWaux PRI_DHW PRIMARY
    172. 0 0 0
    173. *
    174. EQUATIONS 1
    175. * Integrierter Wert = Leistungsziffer
    176. COP_mean = PEF_EL*[43,1]/MAX(0.0001,[43,2])
    177. *
    178. UNIT 44 TYPE 25 PRINTER YEAR
    179. PARAMETERS 4
    180. 8760 START STOP 94
    181. INPUT 4
    182. 43,1 43,2 43,3 COP_mean
    183. DHWaux PRI_DHW PRIMARY COP_mean
    186. ***************************************************************
    187. *! ONLINE GRAPHICS
    188. ***************************************************************
    189. *
    190. UNIT 23 TYPE 65 ONLINE GRAPHICS
    191. PARAMETERS 10
    192. * 1 Number of variables displayed on left axis
    193. 2
    194. * 2 Number of variables displayed on right axis
    195. 1
    196. * 3 Min Y-axis value for left axis
    197. -20
    198. * 4 Max Y-axis value for left axis
    199. 80
    200. * 5 Min Y-axis value for right axis
    201. 2
    202. * 6 Max Y-axis value for right axis
    203. 5
    204. * 7 Number of plots per simulation
    205. 1
    206. * 8 X-axis gridpoints per plot
    207. 12
    208. * 9 shut off Online;
    209. * 0 or more: display Online ; -1 or less: do not display Online
    210. 0
    211. * 10 Output: 10 or more: logical unit number for output file
    212. * -1 or lower: do not generate output file
    213. -1
    214. *
    215. INPUTS 3
    216. TA 22,3 MY_COP
    217. T_luft_außen T_wasser_warm COP
    218. *
    219. LABELS 5
    220. C -
    221. Temperatur
    222. Leistungsziffer
    223. Online
    226. END
    Alles anzeigen



    Quellcode

    1. SUBROUTINE TYPE113 (TIME,XIN,OUT,T,DTDT,PAR,INFO,ICNTRL,*)
    3. C**************************************************************************
    4. C LUFT/WASSER-WÄRMEPUMPE
    5. C**************************************************************************
    7. DOUBLE PRECISION XIN,OUT
    8. INTEGER*4 INFO
    9. DIMENSION PAR(3),XIN(2),OUT(3),INFO(15)
    10. CHARACTER*3 YCHECK(3),OCHECK(3)
    12. C FIRST CALL - DO THE USUAL STUFF
    13. IF (INFO(7).GE.0) GO TO 100
    14. INFO(3)=2
    15. INFO(4)=3
    16. INFO(6)=3
    17. INFO(9)=1
    18. CALL TYPECK(1,INFO,2,3,0)
    20. DATA YCHECK /'TE1','MF1','TE1'/
    21. DATA OCHECK /'TE1','MF1','PW1'/
    22. CALL RCHECK(INFO,YCHECK,OCHECK)
    24. 100 CONTINUE
    26. tv = XIN(1)
    27. temp = XIN(2)
    29. IF (temp <= 2.) THEN
    30. COP_A=0.0555*temp+2.8889
    31. COP_ATop=0.0778*temp+3.544
    32. COP_ABot=0.088*temp+2.618
    33. ELSE IF ((temp > 2.) .AND. (temp <= 10.)) THEN
    34. COP_A=0.1625*temp+2.675
    35. COP_ATop=0.175*temp+3.35
    36. COP_ABot=0.088*temp+2.618
    37. ELSE IF ((temp > 10.) .AND. (temp<=15.)) THEN
    38. COP_A=0.075*temp+3.5
    39. COP_ATop=0.12*temp+3.9
    40. COP_ABot=0.04*temp+3.1
    41. ELSE IF (temp > 15.) THEN
    42. COP_A=0.075*temp+3.5
    43. COP_ATop=0.15*temp+3.45
    44. COP_ABot=0.04*temp+3.1
    45. END IF
    48. IF (COP_A<1.)THEN
    49. COP_A=1.
    50. END IF
    51. IF (COP_ABot<1.)THEN
    52. COP_ABot=1.
    53. END IF
    54. IF (COP_ATop<1.)THEN
    55. COP_ATop=1.
    56. END IF
    58. !50 Grad Celsius Vorlauftemperatur
    59. IF (temp <=2.)THEN
    60. COP_B=0.067*temp+2.371
    61. COP_BTop=0.0556*temp+2.689
    62. COP_BBot=0.033*temp+1.73
    63. ELSE IF ((temp >2.) .AND. (temp <=7.)) THEN
    64. COP_B=0.067*temp+2.371
    65. COP_BTop=0.0846*temp+2.631
    66. COP_BBot=0.14*temp+1.52
    67. ELSE IF ((temp >7.) .AND. (temp <=15.)) THEN
    68. COP_B=0.067*temp+2.371
    69. COP_BTop=0.0846*temp+2.631
    70. COP_BBot=0.044*temp+2.194
    71. ELSE IF(temp > 15.) THEN
    72. COP_B=0.049*temp+2.645
    73. COP_BTop=0.06*temp+3
    74. COP_BBot=0.03*temp+2.4
    75. END IF
    77. IF (COP_B<1.)THEN
    78. COP_B=1.
    79. END IF
    81. IF (COP_BTop<1.)THEN
    82. COP_BTop=1.
    83. END IF
    85. IF (COP_BBot<1.)THEN
    86. COP_BBot=1.
    87. END IF
    89. !COP durch Interpolation Berechnen
    90. COP=COP_A+((COP_B-COP_A)/(50.-35.))*(tv-35.)
    91. COP_Ideal=COP_ATop+((COP_BTop-COP_ATop)/(50.-35.))*(tv-35.)
    92. COP_Un=COP_ABot+((COP_BBot-COP_ABot)/(50.-35.))*(tv-35.)
    94. C SET THE OUTPUTS
    95. OUT(1) = COP
    96. OUT(2) = COP_Ideal
    97. OUT(3) = COP_Un
    98. RETURN 1
    99. END
    Alles anzeigen