SeminarcenterThermostat example¶
In this example, a rule based controller is used to control the model “SeminarcenterThermostat-v0”. At first we import the controller:
[1]:
from energym.examples.Controller import SimpleController
Next, we import Energym and create the simulation environment by specifying the model, a weather file and the number of simulation days.
[2]:
import energym
weather = "DNK_MJ_Horsens1"
env = energym.make("SeminarcenterThermostat-v0", weather=weather, simulation_days=30)
c:\users\psh\documents\gitlab\best\energym\schedules\EProductionSchedule.py:125: SettingWithCopyWarning:
A value is trying to be set on a copy of a slice from a DataFrame.
Try using .loc[row_indexer,col_indexer] = value instead
See the caveats in the documentation: http://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy
- oldmin * random_max_fac * oldmax / delta
[OK] fmi2Instantiate: The Resource location of FMU with instance name %s is %s.
[WARNING] fmi2Instantiate: Argument loggingOn is set to %d
. This is not supported. loggingOn will default to '0'.
[OK] The current working directory is %s
[OK] fmi2Instantiate: Path to fmuUnzipLocation is not %s.
[OK] fmi2Instantiate: Path to fmuUnzipLocation is not %s.
[OK] fmi2Instantiate: Path to fmuUnzipLocation is not %s.
[OK] fmi2Instantiate: Path to fmuUnzipLocation %s
[OK] fmi2Instantiate: Path to fmuResourceLocation %s
[OK] Command executes to copy content of resources folder: %s
[OK] fmi2Instantiate: Path to model description file is %s.
[OK] fmi2Instantiate: The FMU modelIdentifier is %s.
[OK] fmi2Instantiate: The FMU modelGUID is %s.
[OK] fmi2Instantiate: Slave %s is instantiated.
[OK] fmi2Instantiate: Instantiation of %s succeded.
[OK] fmi2EnterInitializationMode: The sockfd is %d.
[OK] fmi2EnterInitializationMode: The port number is %d.
[OK] fmi2EnterInitializationMode: This hostname is %s.
[OK] fmi2EnterInitializationMode: TCPServer Server waiting for clients on port: %d.
[OK] fmi2EnterInitializationMode: The number of input variables is %d.
[OK] fmi2EnterInitializationMode: The number of output variables is %d.
[OK] Get input file from resource folder %s.
[OK] Searching for following pattern %s
[OK] Read directory and search for *.idf, *.epw, or *.idd file.
[OK] Read directory and search for *.idf, *.epw, or *.idd file.
[OK] Read directory and search for *.idf, *.epw, or *.idd file.
[OK] Read directory and search for *.idf, *.epw, or *.idd file.
[OK] Read directory and search for *.idf, *.epw, or *.idd file.
[OK] Read directory and search for *.idf, *.epw, or *.idd file.
[OK] Found matching file %s.
[OK] done searching pattern %s
[OK] Get input file from resource folder %s.
[OK] Searching for following pattern %s
[OK] Read directory and search for *.idf, *.epw, or *.idd file.
[OK] Read directory and search for *.idf, *.epw, or *.idd file.
[OK] Read directory and search for *.idf, *.epw, or *.idd file.
[OK] Found matching file %s.
[OK] done searching pattern %s
[OK] Get input file from resource folder %s.
[OK] Searching for following pattern %s
[OK] Read directory and search for *.idf, *.epw, or *.idd file.
[OK] Read directory and search for *.idf, *.epw, or *.idd file.
[OK] Read directory and search for *.idf, *.epw, or *.idd file.
[OK] Read directory and search for *.idf, *.epw, or *.idd file.
[OK] Found matching file %s.
[OK] done searching pattern %s
[OK] This version uses the **energyplus** command line interface to call the EnergyPlus executable. **RunEPlus.bat** and **runenergyplus** , which were used in earlier versions, were deprecated as of August 2015.
[OK] fmi2EnterInitializationMode: The connection has been accepted.
[OK] fmi2EnterInitializationMode: Slave %s is initialized.
The control inputs can be inspected using the get_inputs_names() method and to construct a controller, we pass the list of inputs and further parameters. This controller determines inputs to get close to the temperature setpoints and uses fixed setpoints during the night.
[3]:
inputs = env.get_inputs_names()[1:] # only thermostat setpoints
print(inputs)
controller = SimpleController(control_list=inputs, lower_tol=0.3, upper_tol=0.8, nighttime_setback=True, nighttime_start=22, nighttime_end=9, nighttime_temp=17)
['Z02_T_Thermostat_sp', 'Z03_T_Thermostat_sp', 'Z04_T_Thermostat_sp', 'Z05_T_Thermostat_sp', 'Z06_T_Thermostat_sp', 'Z08_T_Thermostat_sp', 'Z09_T_Thermostat_sp', 'Z10_T_Thermostat_sp', 'Z11_T_Thermostat_sp', 'Z13_T_Thermostat_sp', 'Z14_T_Thermostat_sp', 'Z15_T_Thermostat_sp', 'Z18_T_Thermostat_sp', 'Z19_T_Thermostat_sp', 'Z20_T_Thermostat_sp', 'Z21_T_Thermostat_sp', 'Z22_T_Thermostat_sp']
To run the simulation, a number of steps is specified (here 144 steps per day for 10 days) and the obtained control inputs are passed to the simulation model with the step() method. To generate some plots later on, we save all the inputs and outputs in lists.
[4]:
steps = 144*10
out_list = []
outputs = env.step(env.sample_random_action())
hour = 0
controls = []
for _ in range(steps):
control = controller.get_control(outputs, 22, hour)
controls +=[ {p:control[p][0] for p in control} ]
outputs = env.step(control)
_,hour,_,_ = env.get_date()
out_list.append(outputs)
Since the inputs and outputs are given as dictionaries and are collected in lists, we can simply load them as a pandas.DataFrame.
[5]:
import pandas as pd
out_df = pd.DataFrame(out_list)
cmd_df = pd.DataFrame(controls)
To generate plots, we can directly get the data from the DataFrames, by using the key names. Displayed are the zone temperatures and the setpoints determined by the controller for zone 1, the external temperature and the total power demand.
[6]:
import matplotlib.pyplot as plt
%matplotlib notebook
f, (ax1,ax2,ax3) = plt.subplots(3,figsize=(10,15))#
ax1.plot(out_df['Z02_T'], 'r')
ax1.plot(out_df['Z02_T_Thermostat_sp_out'], 'b--')
ax1.set_ylabel('Temp')
ax1.set_xlabel('Steps')
ax2.plot(out_df['Ext_T'], 'r')
ax2.set_ylabel('Temp')
ax2.set_xlabel('Steps')
ax3.plot(out_df['Fa_Pw_All'], 'g')
ax3.set_ylabel('Power')
ax3.set_xlabel('Steps')
plt.subplots_adjust(hspace=0.4)
plt.show()
[7]:
env.close()
[OK] fmi2Terminate: fmiFreeInstanceSlave must be called to free the FMU instance.
[OK] fmi2FreeInstance: The function fmi2FreeInstance of instance %s is executed.
[OK] freeInstanceResources: %s will be freed.
[ ]: