""" 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 class Idf: _THERMOSTAT = 'HVACTEMPLATE:THERMOSTAT' _IDEAL_LOAD_AIR_SYSTEM = 'HVACTEMPLATE:ZONE:IDEALLOADSAIRSYSTEM' _SURFACE = 'BUILDINGSURFACE:DETAILED' _CONSTRUCTION = 'CONSTRUCTION' _MATERIAL = 'MATERIAL' _MATERIAL_NOMASS = 'MATERIAL:NOMASS' _ROUGHNESS = "MediumRough" _HOURLY_SCHEDULE = "SCHEDULE:DAY:HOURLY" _ZONE = "ZONE" _LIGHTS = "LIGHTS" _PEOPLE = "PEOPLE" _ELECTRIC_EQUIPMEN = "ELECTRICEQUIPMENT" _INFILTRATION = "ZONEINFILTRATION:DESIGNFLOWRATE" _BUILDING_SURFACE = "BuildingSurfaceDetailed" idf_surfaces = { # todo: make an enum for all the surface types 'Wall': 'wall', 'Ground': 'floor', 'Roof': 'roof' } idf_usage = { # todo: make an enum for all the usage types 'residential': 'residential_building' } def __init__(self, city, idf_file_path, idd_file_path, epw_file_path): self._city = city 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) @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_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=self._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=self._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, usage_zone, schedule_type): for schedule in self._idf.idfobjects[schedule_type]: if schedule.Name == f'{schedule_type} schedules {usage_zone.usage}': return schedule = self._idf.newidfobject(self._HOURLY_SCHEDULE) schedule.Name = f'{schedule_type} schedules {usage_zone.usage}' schedule.Schedule_Type_Limits_Name = 'Any Number' schedule.Hour_1 = usage_zone.schedules[schedule_type]["WD"][0] schedule.Hour_2 = usage_zone.schedules[schedule_type]["WD"][1] schedule.Hour_3 = usage_zone.schedules[schedule_type]["WD"][2] schedule.Hour_4 = usage_zone.schedules[schedule_type]["WD"][3] schedule.Hour_5 = usage_zone.schedules[schedule_type]["WD"][4] schedule.Hour_6 = usage_zone.schedules[schedule_type]["WD"][5] schedule.Hour_7 = usage_zone.schedules[schedule_type]["WD"][6] schedule.Hour_8 = usage_zone.schedules[schedule_type]["WD"][7] schedule.Hour_9 = usage_zone.schedules[schedule_type]["WD"][8] schedule.Hour_10 = usage_zone.schedules[schedule_type]["WD"][9] schedule.Hour_11 = usage_zone.schedules[schedule_type]["WD"][10] schedule.Hour_12 = usage_zone.schedules[schedule_type]["WD"][11] schedule.Hour_13 = usage_zone.schedules[schedule_type]["WD"][12] schedule.Hour_14 = usage_zone.schedules[schedule_type]["WD"][13] schedule.Hour_15 = usage_zone.schedules[schedule_type]["WD"][14] schedule.Hour_16 = usage_zone.schedules[schedule_type]["WD"][15] schedule.Hour_17 = usage_zone.schedules[schedule_type]["WD"][16] schedule.Hour_18 = usage_zone.schedules[schedule_type]["WD"][17] schedule.Hour_19 = usage_zone.schedules[schedule_type]["WD"][18] schedule.Hour_20 = usage_zone.schedules[schedule_type]["WD"][19] schedule.Hour_21 = usage_zone.schedules[schedule_type]["WD"][20] schedule.Hour_22 = usage_zone.schedules[schedule_type]["WD"][21] schedule.Hour_23 = usage_zone.schedules[schedule_type]["WD"][22] schedule.Hour_24 = usage_zone.schedules[schedule_type]["WD"][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 # The constructions should have at least one layer _kwargs = {"Name": thermal_boundary.construction_name, "Outside_Layer": layers[0].material.name} for i in range(1, len(layers) - 1): _kwargs[f'Layer_{i + 1}'] = layers[1].material.name self._idf.newidfobject(self._CONSTRUCTION, **_kwargs) def _add_zone(self, usage_zone): for zone in self._idf.idfobjects['ZONE']: if zone.Name == usage_zone.id: return # todo: what does we need to define a zone in energy plus? self._idf.newidfobject(self._ZONE, Name=usage_zone.id) self._add_heating_system(usage_zone) def _add_thermostat(self, usage_zone): thermostat_name = f'Thermostat {usage_zone.usage}' for thermostat in self._idf.idfobjects[self._THERMOSTAT]: if thermostat == thermostat_name: return thermostat return self._idf.newidfobject(self._THERMOSTAT, Name=thermostat_name, Constant_Heating_Setpoint=usage_zone.heating_setpoint, Constant_Cooling_Setpoint=usage_zone.cooling_setpoint) def _add_heating_system(self, usage_zone): thermostat = self._add_thermostat(usage_zone) # todo: doesn't the air system have name? self._idf.newidfobject(self._IDEAL_LOAD_AIR_SYSTEM, Zone_Name=usage_zone.id, Template_Thermostat_Name=thermostat.Name) def _add_occupancy(self, usage_zone): self._idf.newidfobject(self._PEOPLE, Name=f'{usage_zone.id}_occupancy', Zone_or_ZoneList_Name=usage_zone.id, Number_of_People_Schedule_Name=f'Occupancy schedules {usage_zone.usage}', Number_of_People_Calculation_Method="People", Number_of_People=500, # todo: get people from where? Fraction_Radiant=0.3, # todo: howto get this from InternalGains Activity_Level_Schedule_Name='occupant schedules' ) def _add_equipment(self, usage_zone): self._idf.newidfobject(self._ELECTRIC_EQUIPMENT, Name=f'{usage_zone.id}_electricload', Zone_or_ZoneList_Name=usage_zone.id, Schedule_Name=f'Electrical schedules {usage_zone.usage}', # todo: add electrical schedules Design_Level_Calculation_Method='EquipmentLevel', Design_Level=566000 # todo: change it from usage catalog ) def _add_infiltration(self, usage_zone): for zone in self._idf.idfobjects["ZONE"]: self._idf.newidfobject(self._INFILTRATION, Name=f'{usage_zone.id}_infiltration', Zone_or_ZoneList_Name=usage_zone.id, Schedule_Name=f'Infiltration schedules {usage_zone.usage}', Design_Flow_Rate_Calculation_Method='AirChanges/Hour', Air_Changes_per_Hour=0.35, # todo: change it from usage catalog Constant_Term_Coefficient=0.606, # todo: change it from usage catalog Temperature_Term_Coefficient=3.6359996E-02, # todo: change it from usage catalog Velocity_Term_Coefficient=0.1177165, # todo: change it from usage catalog Velocity_Squared_Term_Coefficient=0.0000000E+00 # todo: change it from usage catalog ) def export(self, city, export_type="Surfaces"): """ Export the idf file into the given path export type = "Surfaces|Block" """ for building in city.buildings: for usage_zone in building.usage_zones: self._add_schedule(usage_zone, "Infiltration") self._add_schedule(usage_zone, "Lights") self._add_schedule(usage_zone, "Occupancy") self._add_zone(usage_zone) self._add_heating_system(usage_zone) self._add_equipment(usage_zone) self._add_construction(usage_zone) if export_type == "Surfaces": self._add_surfaces(building) else: self._add_block(building) def _add_block(self, building): _points = self._matrix_to_2d_list(building.foot_print.coordinates) self._idf.add_block(name=building.name, coordinates=_points, height=building.max_height, num_stories=int(building.storeys_above_ground)) for surface in self._idf.idfobjects[self._SURFACE]: for thermal_zone in building.thermal_zones: for boundary in thermal_zone.bounded: if surface.Type == self.idf_surfaces[boundary.surface.type]: surface.Construction_Name = boundary.construction_name break for usage_zone in thermal_zone.usage_zones: surface.Zone_Name = usage_zone.id break break self._idf.intersect_match() def _add_surfaces(self, building): for thermal_zone in building.thermal_zones: for boundary in thermal_zone.bounded: idf_surface = self.idf_surfaces[boundary.surface.type] for usage_zone in thermal_zone.usage_zones: surface = self._idf.newidfobject(self._SURFACE, Name=f'{boundary.surface.name}', Surface_Type=idf_surface, Zone_Name=usage_zone.id, Construction_Name=boundary.construction_name) coordinates = self._matrix_to_list(boundary.surface.coordinates) surface.setcoords(coordinates) self._idf.intersect_match()