From 3bf674d107727b094046d094a1f07d761fd29866 Mon Sep 17 00:00:00 2001 From: guille Date: Wed, 8 Jun 2022 12:55:37 -0400 Subject: [PATCH] now schedules are re-calculated on change --- .../building_demand/thermal_zone.py | 412 +++++++++--------- 1 file changed, 206 insertions(+), 206 deletions(-) diff --git a/city_model_structure/building_demand/thermal_zone.py b/city_model_structure/building_demand/thermal_zone.py index 784bb0ab..eff21366 100644 --- a/city_model_structure/building_demand/thermal_zone.py +++ b/city_model_structure/building_demand/thermal_zone.py @@ -244,7 +244,7 @@ class ThermalZone: return None self._usage = '' for usage_zone in self._parent_internal_zone.usage_zones: - self._usage += str(round(usage_zone.percentage) * 100) + '-' + usage_zone.usage + '_' + self._usage += str(round(usage_zone.percentage*100)) + '-' + usage_zone.usage + '_' self._usage = self._usage[:-1] return self._usage @@ -308,39 +308,39 @@ class ThermalZone: """ if self.usage_zones is None: return None - if self._occupancy is None: - self._occupancy = Occupancy() - _occupancy_density = 0 - _convective_part = 0 - _radiative_part = 0 - _latent_part = 0 - for usage_zone in self.usage_zones: - if usage_zone.occupancy is None: - return None - _occupancy_density += usage_zone.percentage * usage_zone.occupancy.occupancy_density - if usage_zone.occupancy.sensible_convective_internal_gain is not None: - _convective_part += usage_zone.percentage * usage_zone.occupancy.sensible_convective_internal_gain - _radiative_part += usage_zone.percentage * usage_zone.occupancy.sensible_radiative_internal_gain - _latent_part += usage_zone.percentage * usage_zone.occupancy.latent_internal_gain - self._occupancy.occupancy_density = _occupancy_density - self._occupancy.sensible_convective_internal_gain = _convective_part - self._occupancy.sensible_radiative_internal_gain = _radiative_part - self._occupancy.latent_internal_gain = _latent_part - _occupancy_reference = self.usage_zones[0].occupancy - if _occupancy_reference.occupancy_schedules is not None: - _schedules = [] - for i_schedule in range(0, len(_occupancy_reference.occupancy_schedules)): - schedule = copy.deepcopy(_occupancy_reference.occupancy_schedules[i_schedule]) - new_values = [] - for i_value in range(0, len(_occupancy_reference.occupancy_schedules[i_schedule].values)): - _new_value = 0 - for usage_zone in self.usage_zones: - _new_value += usage_zone.percentage * usage_zone.occupancy.occupancy_schedules[i_schedule].values[i_value] - new_values.append(_new_value) - schedule.values = new_values - _schedules.append(schedule) - self._occupancy.occupancy_schedules = _schedules + self._occupancy = Occupancy() + _occupancy_density = 0 + _convective_part = 0 + _radiative_part = 0 + _latent_part = 0 + for usage_zone in self.usage_zones: + if usage_zone.occupancy is None: + return None + _occupancy_density += usage_zone.percentage * usage_zone.occupancy.occupancy_density + if usage_zone.occupancy.sensible_convective_internal_gain is not None: + _convective_part += usage_zone.percentage * usage_zone.occupancy.sensible_convective_internal_gain + _radiative_part += usage_zone.percentage * usage_zone.occupancy.sensible_radiative_internal_gain + _latent_part += usage_zone.percentage * usage_zone.occupancy.latent_internal_gain + self._occupancy.occupancy_density = _occupancy_density + self._occupancy.sensible_convective_internal_gain = _convective_part + self._occupancy.sensible_radiative_internal_gain = _radiative_part + self._occupancy.latent_internal_gain = _latent_part + + _occupancy_reference = self.usage_zones[0].occupancy + if _occupancy_reference.occupancy_schedules is not None: + _schedules = [] + for i_schedule in range(0, len(_occupancy_reference.occupancy_schedules)): + schedule = copy.deepcopy(_occupancy_reference.occupancy_schedules[i_schedule]) + new_values = [] + for i_value in range(0, len(_occupancy_reference.occupancy_schedules[i_schedule].values)): + _new_value = 0 + for usage_zone in self.usage_zones: + _new_value += usage_zone.percentage * usage_zone.occupancy.occupancy_schedules[i_schedule].values[i_value] + new_values.append(_new_value) + schedule.values = new_values + _schedules.append(schedule) + self._occupancy.occupancy_schedules = _schedules return self._occupancy @property @@ -351,47 +351,47 @@ class ThermalZone: """ if self.usage_zones is None: return None - if self._lighting is None: - self._lighting = Lighting() - _lighting_density = 0 - _convective_part = 0 - _radiative_part = 0 - _latent_part = 0 - for usage_zone in self.usage_zones: - if usage_zone.lighting is None: - return None - _lighting_density += usage_zone.percentage * usage_zone.lighting.density - if usage_zone.lighting.convective_fraction is not None: - _convective_part += usage_zone.percentage * usage_zone.lighting.density \ - * usage_zone.lighting.convective_fraction - _radiative_part += usage_zone.percentage * usage_zone.lighting.density \ - * usage_zone.lighting.radiative_fraction - _latent_part += usage_zone.percentage * usage_zone.lighting.density \ - * usage_zone.lighting.latent_fraction - self._lighting.density = _lighting_density - if _lighting_density > 0: - self._lighting.convective_fraction = _convective_part / _lighting_density - self._lighting.radiative_fraction = _radiative_part / _lighting_density - self._lighting.latent_fraction = _latent_part / _lighting_density - else: - self._lighting.convective_fraction = 0 - self._lighting.radiative_fraction = 0 - self._lighting.latent_fraction = 0 - _lighting_reference = self.usage_zones[0].lighting - if _lighting_reference.schedules is not None: - _schedules = [] - for i_schedule in range(0, len(_lighting_reference.schedules)): - schedule = copy.deepcopy(_lighting_reference.schedules[i_schedule]) - new_values = [] - for i_value in range(0, len(_lighting_reference.schedules[i_schedule].values)): - _new_value = 0 - for usage_zone in self.usage_zones: - _new_value += usage_zone.percentage * usage_zone.lighting.schedules[i_schedule].values[i_value] - new_values.append(_new_value) - schedule.values = new_values - _schedules.append(schedule) - self._lighting.schedules = _schedules + self._lighting = Lighting() + _lighting_density = 0 + _convective_part = 0 + _radiative_part = 0 + _latent_part = 0 + for usage_zone in self.usage_zones: + if usage_zone.lighting is None: + return None + _lighting_density += usage_zone.percentage * usage_zone.lighting.density + if usage_zone.lighting.convective_fraction is not None: + _convective_part += usage_zone.percentage * usage_zone.lighting.density \ + * usage_zone.lighting.convective_fraction + _radiative_part += usage_zone.percentage * usage_zone.lighting.density \ + * usage_zone.lighting.radiative_fraction + _latent_part += usage_zone.percentage * usage_zone.lighting.density \ + * usage_zone.lighting.latent_fraction + self._lighting.density = _lighting_density + if _lighting_density > 0: + self._lighting.convective_fraction = _convective_part / _lighting_density + self._lighting.radiative_fraction = _radiative_part / _lighting_density + self._lighting.latent_fraction = _latent_part / _lighting_density + else: + self._lighting.convective_fraction = 0 + self._lighting.radiative_fraction = 0 + self._lighting.latent_fraction = 0 + + _lighting_reference = self.usage_zones[0].lighting + if _lighting_reference.schedules is not None: + _schedules = [] + for i_schedule in range(0, len(_lighting_reference.schedules)): + schedule = copy.deepcopy(_lighting_reference.schedules[i_schedule]) + new_values = [] + for i_value in range(0, len(_lighting_reference.schedules[i_schedule].values)): + _new_value = 0 + for usage_zone in self.usage_zones: + _new_value += usage_zone.percentage * usage_zone.lighting.schedules[i_schedule].values[i_value] + new_values.append(_new_value) + schedule.values = new_values + _schedules.append(schedule) + self._lighting.schedules = _schedules return self._lighting @property @@ -402,47 +402,47 @@ class ThermalZone: """ if self.usage_zones is None: return None - if self._appliances is None: - self._appliances = Appliances() - _appliances_density = 0 - _convective_part = 0 - _radiative_part = 0 - _latent_part = 0 - for usage_zone in self.usage_zones: - if usage_zone.appliances is None: - return None - _appliances_density += usage_zone.percentage * usage_zone.appliances.density - if usage_zone.appliances.convective_fraction is not None: - _convective_part += usage_zone.percentage * usage_zone.appliances.density \ - * usage_zone.appliances.convective_fraction - _radiative_part += usage_zone.percentage * usage_zone.appliances.density \ - * usage_zone.appliances.radiative_fraction - _latent_part += usage_zone.percentage * usage_zone.appliances.density \ - * usage_zone.appliances.latent_fraction - self._appliances.density = _appliances_density - if _appliances_density > 0: - self._appliances.convective_fraction = _convective_part / _appliances_density - self._appliances.radiative_fraction = _radiative_part / _appliances_density - self._appliances.latent_fraction = _latent_part / _appliances_density - else: - self._appliances.convective_fraction = 0 - self._appliances.radiative_fraction = 0 - self._appliances.latent_fraction = 0 - _appliances_reference = self.usage_zones[0].appliances - if _appliances_reference.schedules is not None: - _schedules = [] - for i_schedule in range(0, len(_appliances_reference.schedules)): - schedule = copy.deepcopy(_appliances_reference.schedules[i_schedule]) - new_values = [] - for i_value in range(0, len(_appliances_reference.schedules[i_schedule].values)): - _new_value = 0 - for usage_zone in self.usage_zones: - _new_value += usage_zone.percentage * usage_zone.appliances.schedules[i_schedule].values[i_value] - new_values.append(_new_value) - schedule.values = new_values - _schedules.append(schedule) - self._appliances.schedules = _schedules + self._appliances = Appliances() + _appliances_density = 0 + _convective_part = 0 + _radiative_part = 0 + _latent_part = 0 + for usage_zone in self.usage_zones: + if usage_zone.appliances is None: + return None + _appliances_density += usage_zone.percentage * usage_zone.appliances.density + if usage_zone.appliances.convective_fraction is not None: + _convective_part += usage_zone.percentage * usage_zone.appliances.density \ + * usage_zone.appliances.convective_fraction + _radiative_part += usage_zone.percentage * usage_zone.appliances.density \ + * usage_zone.appliances.radiative_fraction + _latent_part += usage_zone.percentage * usage_zone.appliances.density \ + * usage_zone.appliances.latent_fraction + self._appliances.density = _appliances_density + if _appliances_density > 0: + self._appliances.convective_fraction = _convective_part / _appliances_density + self._appliances.radiative_fraction = _radiative_part / _appliances_density + self._appliances.latent_fraction = _latent_part / _appliances_density + else: + self._appliances.convective_fraction = 0 + self._appliances.radiative_fraction = 0 + self._appliances.latent_fraction = 0 + + _appliances_reference = self.usage_zones[0].appliances + if _appliances_reference.schedules is not None: + _schedules = [] + for i_schedule in range(0, len(_appliances_reference.schedules)): + schedule = copy.deepcopy(_appliances_reference.schedules[i_schedule]) + new_values = [] + for i_value in range(0, len(_appliances_reference.schedules[i_schedule].values)): + _new_value = 0 + for usage_zone in self.usage_zones: + _new_value += usage_zone.percentage * usage_zone.appliances.schedules[i_schedule].values[i_value] + new_values.append(_new_value) + schedule.values = new_values + _schedules.append(schedule) + self._appliances.schedules = _schedules return self._appliances @property @@ -453,54 +453,54 @@ class ThermalZone: """ if self.usage_zones is None: return None - if self._internal_gains is None: - _internal_gain = InternalGain() - _days = [cte.MONDAY, cte.TUESDAY, cte.WEDNESDAY, cte.THURSDAY, cte.FRIDAY, cte.SATURDAY, cte.SUNDAY, cte.HOLIDAY] - _average_internal_gain = 0 - _convective_fraction = 0 - _radiative_fraction = 0 - _latent_fraction = 0 - _schedules = None - _base_schedule = Schedule() - _base_schedule.type = cte.INTERNAL_GAINS - _base_schedule.time_range = cte.DAY - _base_schedule.time_step = cte.HOUR - _base_schedule.data_type = cte.ANY_NUMBER - _schedules_defined = True - values = numpy.zeros(24, 8) - for usage_zone in self.usage_zones: - for internal_gain in usage_zone.internal_gains: - _average_internal_gain += internal_gain.average_internal_gain * usage_zone.percentage - _convective_fraction += internal_gain.average_internal_gain * usage_zone.percentage \ - * internal_gain.convective_fraction - _radiative_fraction += internal_gain.average_internal_gain * usage_zone.percentage \ - * internal_gain.radiative_fraction - _latent_fraction += internal_gain.average_internal_gain * usage_zone.percentage \ - * internal_gain.latent_fraction - for usage_zone in self.usage_zones: - for internal_gain in usage_zone.internal_gains: - if len(internal_gain.schedules) == 0: - _schedules_defined = False - break - for day, _schedule in enumerate(internal_gain.schedules): - for v, value in enumerate(_schedule.values): - values[v, day] += value * usage_zone.percentage - if _schedules_defined: - _schedules = [] - for day in range(0, len(_days)): - _schedule = copy.deepcopy(_base_schedule) - _schedule.day_types = [_days[day]] - _schedule.values = values[:day] - _schedules.append(_schedule) + _internal_gain = InternalGain() + _days = [cte.MONDAY, cte.TUESDAY, cte.WEDNESDAY, cte.THURSDAY, cte.FRIDAY, cte.SATURDAY, cte.SUNDAY, cte.HOLIDAY] + _average_internal_gain = 0 + _convective_fraction = 0 + _radiative_fraction = 0 + _latent_fraction = 0 + _schedules = None + _base_schedule = Schedule() + _base_schedule.type = cte.INTERNAL_GAINS + _base_schedule.time_range = cte.DAY + _base_schedule.time_step = cte.HOUR + _base_schedule.data_type = cte.ANY_NUMBER + _schedules_defined = True + values = numpy.zeros(24, 8) + for usage_zone in self.usage_zones: + for internal_gain in usage_zone.internal_gains: + _average_internal_gain += internal_gain.average_internal_gain * usage_zone.percentage + _convective_fraction += internal_gain.average_internal_gain * usage_zone.percentage \ + * internal_gain.convective_fraction + _radiative_fraction += internal_gain.average_internal_gain * usage_zone.percentage \ + * internal_gain.radiative_fraction + _latent_fraction += internal_gain.average_internal_gain * usage_zone.percentage \ + * internal_gain.latent_fraction + for usage_zone in self.usage_zones: + for internal_gain in usage_zone.internal_gains: + if len(internal_gain.schedules) == 0: + _schedules_defined = False + break + for day, _schedule in enumerate(internal_gain.schedules): + for v, value in enumerate(_schedule.values): + values[v, day] += value * usage_zone.percentage - _internal_gain.convective_fraction = _convective_fraction / _average_internal_gain - _internal_gain.radiative_fraction = _radiative_fraction / _average_internal_gain - _internal_gain.latent_fraction = _latent_fraction / _average_internal_gain - _internal_gain.average_internal_gain = _average_internal_gain - _internal_gain.type = 'mean_value' - _internal_gain.schedules = _schedules - self._internal_gains = [_internal_gain] + if _schedules_defined: + _schedules = [] + for day in range(0, len(_days)): + _schedule = copy.deepcopy(_base_schedule) + _schedule.day_types = [_days[day]] + _schedule.values = values[:day] + _schedules.append(_schedule) + + _internal_gain.convective_fraction = _convective_fraction / _average_internal_gain + _internal_gain.radiative_fraction = _radiative_fraction / _average_internal_gain + _internal_gain.latent_fraction = _latent_fraction / _average_internal_gain + _internal_gain.average_internal_gain = _average_internal_gain + _internal_gain.type = 'mean_value' + _internal_gain.schedules = _schedules + self._internal_gains = [_internal_gain] return self._internal_gains @property @@ -511,55 +511,55 @@ class ThermalZone: """ if self.usage_zones is None: return None - if self._thermal_control is None: - self._thermal_control = ThermalControl() - _mean_heating_set_point = 0 - _heating_set_back = 0 - _mean_cooling_set_point = 0 - for usage_zone in self.usage_zones: - _mean_heating_set_point += usage_zone.percentage * usage_zone.thermal_control.mean_heating_set_point - _heating_set_back += usage_zone.percentage * usage_zone.thermal_control.heating_set_back - _mean_cooling_set_point += usage_zone.percentage * usage_zone.thermal_control.mean_cooling_set_point - self._thermal_control.mean_heating_set_point = _mean_heating_set_point - self._thermal_control.heating_set_back = _heating_set_back - self._thermal_control.mean_cooling_set_point = _mean_cooling_set_point - _thermal_control_reference = self.usage_zones[0].thermal_control - _types_reference = [] - if _thermal_control_reference.hvac_availability_schedules is not None: - _types_reference.append([cte.HVAC_AVAILABILITY, _thermal_control_reference.hvac_availability_schedules]) - if _thermal_control_reference.heating_set_point_schedules is not None: - _types_reference.append([cte.HEATING_SET_POINT, _thermal_control_reference.heating_set_point_schedules]) - if _thermal_control_reference.cooling_set_point_schedules is not None: - _types_reference.append([cte.COOLING_SET_POINT, _thermal_control_reference.cooling_set_point_schedules]) + self._thermal_control = ThermalControl() + _mean_heating_set_point = 0 + _heating_set_back = 0 + _mean_cooling_set_point = 0 + for usage_zone in self.usage_zones: + _mean_heating_set_point += usage_zone.percentage * usage_zone.thermal_control.mean_heating_set_point + _heating_set_back += usage_zone.percentage * usage_zone.thermal_control.heating_set_back + _mean_cooling_set_point += usage_zone.percentage * usage_zone.thermal_control.mean_cooling_set_point + self._thermal_control.mean_heating_set_point = _mean_heating_set_point + self._thermal_control.heating_set_back = _heating_set_back + self._thermal_control.mean_cooling_set_point = _mean_cooling_set_point - for i_type in range(0, len(_types_reference)): - _schedules = [] - _schedule_type = _types_reference[i_type][1] - for i_schedule in range(0, len(_schedule_type)): - schedule = copy.deepcopy(_schedule_type[i_schedule]) - new_values = [] - for i_value in range(0, len(_schedule_type[i_schedule].values)): - _new_value = 0 - for usage_zone in self.usage_zones: - if _types_reference[i_type][0] == cte.HVAC_AVAILABILITY: - _new_value += usage_zone.percentage * \ - usage_zone.thermal_control.hvac_availability_schedules[i_schedule].values[i_value] - elif _types_reference[i_type][0] == cte.HEATING_SET_POINT: - _new_value += usage_zone.percentage * \ - usage_zone.thermal_control.heating_set_point_schedules[i_schedule].values[i_value] - elif _types_reference[i_type][0] == cte.COOLING_SET_POINT: - _new_value += usage_zone.percentage * \ - usage_zone.thermal_control.cooling_set_point_schedules[i_schedule].values[i_value] - new_values.append(_new_value) - schedule.values = new_values - _schedules.append(schedule) - if i_type == 0: - self._thermal_control.hvac_availability_schedules = _schedules - elif i_type == 1: - self._thermal_control.heating_set_point_schedules = _schedules - elif i_type == 2: - self._thermal_control.cooling_set_point_schedules = _schedules + _thermal_control_reference = self.usage_zones[0].thermal_control + _types_reference = [] + if _thermal_control_reference.hvac_availability_schedules is not None: + _types_reference.append([cte.HVAC_AVAILABILITY, _thermal_control_reference.hvac_availability_schedules]) + if _thermal_control_reference.heating_set_point_schedules is not None: + _types_reference.append([cte.HEATING_SET_POINT, _thermal_control_reference.heating_set_point_schedules]) + if _thermal_control_reference.cooling_set_point_schedules is not None: + _types_reference.append([cte.COOLING_SET_POINT, _thermal_control_reference.cooling_set_point_schedules]) + + for i_type in range(0, len(_types_reference)): + _schedules = [] + _schedule_type = _types_reference[i_type][1] + for i_schedule in range(0, len(_schedule_type)): + schedule = copy.deepcopy(_schedule_type[i_schedule]) + new_values = [] + for i_value in range(0, len(_schedule_type[i_schedule].values)): + _new_value = 0 + for usage_zone in self.usage_zones: + if _types_reference[i_type][0] == cte.HVAC_AVAILABILITY: + _new_value += usage_zone.percentage * \ + usage_zone.thermal_control.hvac_availability_schedules[i_schedule].values[i_value] + elif _types_reference[i_type][0] == cte.HEATING_SET_POINT: + _new_value += usage_zone.percentage * \ + usage_zone.thermal_control.heating_set_point_schedules[i_schedule].values[i_value] + elif _types_reference[i_type][0] == cte.COOLING_SET_POINT: + _new_value += usage_zone.percentage * \ + usage_zone.thermal_control.cooling_set_point_schedules[i_schedule].values[i_value] + new_values.append(_new_value) + schedule.values = new_values + _schedules.append(schedule) + if i_type == 0: + self._thermal_control.hvac_availability_schedules = _schedules + elif i_type == 1: + self._thermal_control.heating_set_point_schedules = _schedules + elif i_type == 2: + self._thermal_control.cooling_set_point_schedules = _schedules return self._thermal_control