""" TestOccupancyFactory test and validate the city model structure schedules parameters SPDX - License - Identifier: LGPL - 3.0 - or -later Copyright © 2020 Project Author Soroush Samareh Abolhassani - soroush.samarehabolhassani@mail.concordia.ca """ from geomeppy import IDF import os import esoreader from pathlib import Path import helpers.constants as cte class Idf: _THERMOSTAT = 'HVACTEMPLATE:THERMOSTAT' _IDEAL_LOAD_AIR_SYSTEM = 'HVACTEMPLATE:ZONE:IDEALLOADSAIRSYSTEM' _SURFACE = 'BUILDINGSURFACE:DETAILED' _CONSTRUCTION = 'CONSTRUCTION' _MATERIAL = 'MATERIAL' _MATERIAL_NOMASS = 'MATERIAL:NOMASS' idf_surfaces = { 'Wall': 'wall', 'Ground': 'floor', 'Roof': 'roof' } idf_usage = { 'residential': 'residential_building' } def __init__(self, idf_file_path, idd_file_path, epw_file_path): self._idd_file_path = str(idd_file_path) self._idf_file_path = str(idf_file_path) self._epw_file_path = str(epw_file_path) IDF.setiddname(self._idd_file_path) self._idf = IDF(self._idf_file_path, self._epw_file_path) self._idf.epw = self._epw_file_path def _add_material(self, layer): for material in self._idf.idfobjects[self._MATERIAL]: if material.Name == layer.material.name: return for material in self._idf.idfobjects[self._MATERIAL_NOMASS]: if material.Name == layer.material.name: return if layer.material.no_mass: self._idf.newidfobject(self._MATERIAL_NOMASS, Name=layer.material.name, Roughness=cte.ROUGHNESS, Thermal_Resistance=layer.material.thermal_resistance, Thermal_Absorptance=layer.material.thermal_absorptance, Solar_Absorptance=layer.material.solar_absorptance, Visible_Absorptance=layer.material.visible_absorptance ) else: self._idf.newidfobject(self._MATERIAL, Name=layer.material.name, Roughness=cte.ROUGHNESS, Thickness=layer.thickness, Conductivity=layer.material.conductivity, Density=layer.material.density, Specific_Heat=layer.material.specific_heat, Thermal_Absorptance=layer.material.thermal_absorptance, Solar_Absorptance=layer.material.solar_absorptance, Visible_Absorptance=layer.material.visible_absorptance ) def add_schedule(self, building): schedule_occupancy = self._idf.newidfobject("SCHEDULE:DAY:HOURLY".upper()) schedule_occupancy.Name = 'occupant schedules' schedule_occupancy.Hour_1 = building.usage_zones[0].occupants.occupant_schedule[0] schedule_occupancy.Hour_2 = building.usage_zones[0].occupants.occupant_schedule[1] schedule_occupancy.Hour_3 = building.usage_zones[0].occupants.occupant_schedule[2] schedule_occupancy.Hour_4 = building.usage_zones[0].occupants.occupant_schedule[3] schedule_occupancy.Hour_5 = building.usage_zones[0].occupants.occupant_schedule[4] schedule_occupancy.Hour_6 = building.usage_zones[0].occupants.occupant_schedule[5] schedule_occupancy.Hour_7 = building.usage_zones[0].occupants.occupant_schedule[6] schedule_occupancy.Hour_8 = building.usage_zones[0].occupants.occupant_schedule[7] schedule_occupancy.Hour_9 = building.usage_zones[0].occupants.occupant_schedule[8] schedule_occupancy.Hour_10 = building.usage_zones[0].occupants.occupant_schedule[9] schedule_occupancy.Hour_11 = building.usage_zones[0].occupants.occupant_schedule[10] schedule_occupancy.Hour_12 = building.usage_zones[0].occupants.occupant_schedule[11] schedule_occupancy.Hour_13 = building.usage_zones[0].occupants.occupant_schedule[12] schedule_occupancy.Hour_14 = building.usage_zones[0].occupants.occupant_schedule[13] schedule_occupancy.Hour_15 = building.usage_zones[0].occupants.occupant_schedule[14] schedule_occupancy.Hour_16 = building.usage_zones[0].occupants.occupant_schedule[15] schedule_occupancy.Hour_17 = building.usage_zones[0].occupants.occupant_schedule[16] schedule_occupancy.Hour_18 = building.usage_zones[0].occupants.occupant_schedule[17] schedule_occupancy.Hour_19 = building.usage_zones[0].occupants.occupant_schedule[18] schedule_occupancy.Hour_20 = building.usage_zones[0].occupants.occupant_schedule[19] schedule_occupancy.Hour_21 = building.usage_zones[0].occupants.occupant_schedule[20] schedule_occupancy.Hour_22 = building.usage_zones[0].occupants.occupant_schedule[21] schedule_occupancy.Hour_23 = building.usage_zones[0].occupants.occupant_schedule[22] schedule_occupancy.Hour_24 = building.usage_zones[0].occupants.occupant_schedule[23] def _add_construction(self, thermal_boundary): for construction in self._idf.idfobjects[self._CONSTRUCTION]: if construction.Name == thermal_boundary.construction_name: return for layer in thermal_boundary.layers: self._add_material(layer) layers = thermal_boundary.layers if len(layers) == 2: self._idf.newidfobject(self._CONSTRUCTION, Name=thermal_boundary.construction_name, Outside_Layer=layers[0].material.name, Layer_2=layers[1].material.name) elif len(thermal_boundary.layers) == 3: self._idf.newidfobject(self._CONSTRUCTION, Name=thermal_boundary.construction_name, Outside_Layer=layers[0].material.name, Layer_2=layers[1].material.name, Layer_3=layers[2].material.name) elif len(thermal_boundary.layers) == 4: self._idf.newidfobject(self._CONSTRUCTION, Name=thermal_boundary.construction_name, Outside_Layer=layers[0].material.name, Layer_2=layers[1].material.name, Layer_3=layers[2].material.name, Layer_4=layers[3].material.name) else: raise Exception("This libs version cannot handle more than 4 layers") def _add_heating_system(self, building): for usage_zone in building.usage_zones: thermostat_name = f'Thermostat {building.name}' # todo: this will fail for more than one usage zone static_thermostat = self._idf.newidfobject(self._THERMOSTAT, Name=thermostat_name, Constant_Heating_Setpoint=usage_zone.heating_setpoint, Constant_Cooling_Setpoint=usage_zone.cooling_setpoint, ) for zone in self._idf.idfobjects['ZONE']: if zone.Name.find(building.name) != -1: self._idf.newidfobject(self._IDEAL_LOAD_AIR_SYSTEM, Zone_Name=zone.Name, Template_Thermostat_Name=static_thermostat.Name, ) @staticmethod def _matrix_to_list(points): points_list = [] for point in points: point_tuple = (point[0], point[1], point[2]) points_list.append(point_tuple) return points_list @staticmethod def _matrix_to_2d_list(points): points_list = [] for point in points: point_tuple = (point[0], point[1]) points_list.append(point_tuple) return points_list def add_block(self, building): _points = IdfHelper._matrix_to_2d_list(building.foot_print.points) self._idf.add_block(name=building.name, coordinates=_points, height=building.max_height, num_stories=int(building.storeys_above_ground)) self._add_heating_system(building) self._idf.intersect_match() def add_surfaces(self, building): index = 0 for zone in building.thermal_zones: zone_name = f'Building {building.name} usage zone {index}' self._idf.newidfobject('ZONE', Name=zone_name) for boundary in zone.bounded: self._add_construction(boundary) idf_surface = self.idf_surfaces[boundary.surface.type] wall = self._idf.newidfobject(self._SURFACE, Name=f'{building.name}-{boundary.surface.name}', Surface_Type=idf_surface, Zone_Name=zone_name, Construction_Name=boundary.construction_name) coordinates = IdfHelper._matrix_to_list(boundary.surface.points) wall.setcoords(coordinates) index += 1 self._add_heating_system(building) self._idf.intersect_match() def add_lighting(self): for zone in self._idf.idfobjects["ZONE"]: self._idf.newidfobject("LIGHTS", Name=zone.Name + "_" + "lighting", Zone_or_ZoneList_Name=zone.Name, Schedule_Name="Lighting_Weekday", Design_Level_Calculation_Method='Watts/Area', Watts_per_Zone_Floor_Area=15, Return_Air_Fraction=0, Fraction_Radiant=0.7, Fraction_Visible=0.2, Fraction_Replaceable=1, # Lights.End-Use_Subcategory='LightsWired', Return_Air_Fraction_Calculated_from_Plenum_Temperature='No' ) def add_occupancy(self): for zone in self._idf.idfobjects["ZONE"]: self._idf.newidfobject("PEOPLE", Name=zone.Name + "_" + "schedules", Zone_or_ZoneList_Name=zone.Name, Number_of_People_Schedule_Name='occupant schedules', Number_of_People_Calculation_Method="People", Number_of_People=500, # Zone_Floor_Area_per_Person=31.41, # People_per_Zone_Floor_Area=17*0.05, # Zone_Floor_Area_per_Person=1.65880764E+01, Fraction_Radiant=0.3, Activity_Level_Schedule_Name='occupant schedules' ) def add_equipment(self): for zone in self._idf.idfobjects["ZONE"]: self._idf.newidfobject("ELECTRICEQUIPMENT", Name=zone.Name + "_" + 'electricload', Zone_or_ZoneList_Name=zone.Name, Schedule_Name='ElectricalEquipment', Design_Level_Calculation_Method='EquipmentLevel', Design_Level=566000 ) def add_infiltration(self): for zone in self._idf.idfobjects["ZONE"]: self._idf.newidfobject("ZONEINFILTRATION:DESIGNFLOWRATE", Name=zone.Name + "_" + "infiltration", Zone_or_ZoneList_Name=zone.Name, Schedule_Name='Infiltration_schedule', Design_Flow_Rate_Calculation_Method='AirChanges/Hour', Air_Changes_per_Hour=0.35, Constant_Term_Coefficient=0.606, Temperature_Term_Coefficient=3.6359996E-02, Velocity_Term_Coefficient=0.1177165, Velocity_Squared_Term_Coefficient=0.0000000E+00 ) def run(self, output_directory, window_ratio=0.35, display_render=False, output_prefix=None, keep_file=None): self._idf.set_default_constructions() self._idf.set_wwr(window_ratio, construction="Project External Window") self._idf.translate_to_origin() if display_render: self._idf.view_model() # Run self._idf.newidfobject("OUTPUT:METER", Key_Name="Heating:DistrictHeating", Reporting_Frequency="hourly") self._idf.newidfobject("OUTPUT:METER", Key_Name="Cooling:DistrictCooling", Reporting_Frequency="hourly") idf_path = None if keep_file is not None: idf_path = (keep_file / 'in.idf').resolve() self._idf.saveas(str(idf_path)) if idf_path is None: idf_path = (Path(__file__).parent / 'in.idf').resolve() # There is a bug in the IDF class, when called, it return an error, as a work around we call call energy+ directly run_command = f"C:\EnergyPlusV9-4-0\energyplus.exe --weather {self._epw_file_path} --output-directory {output_directory} --idd " \ f"{self._idd_file_path} --expandobjects --output-prefix {output_prefix} {idf_path}" os.system(run_command) if keep_file is None: os.remove(idf_path) return @staticmethod def read_eso(eso_file_path): dd, data = esoreader.read(eso_file_path) list_values = [v for v in data.values()] heating = [(float(x)) / 3600000.0 for x in list_values[0]] cooling = [(float(x)) / 3600000.0 for x in list_values[1]] return heating, cooling