Integrated daylight and energy simulation in LBT

mostapha · January 7, 2021, 3:44pm

As I said above, this is how the DA is calculated. There is already a component for post-processing annual daylight simulation results. See the code and modify it for your purpose. For every hour that there is not enough daylight you should set the lighting schedule to 1. Here is the code:

github.com

ladybug-tools/honeybee-grasshopper-radiance/blob/7d0f2ae7aff468530d0098a396714dd5cb6a2a78/honeybee_grasshopper_radiance/src/HB Annual Daylight Metrics.py#L94-L113


def annual_metrics(ill_file, occ_pattern, total_occupied_hours,
                   threshold=300, min_t=100, max_t=2000):
    """Compute annual metrics for a given result file."""
    da = []
    udi = []
    with open(ill_file) as results:
        for pt_res in results:
            pda = 0
            pudi = 0
            for is_occ, hourly_res in zip(occ_pattern, pt_res.split()):
                if is_occ < 0.5:
                    continue
                value = float(hourly_res)
                if value > threshold:
                    pda += 1
                if min_t <= value <= max_t:
                    pudi += 1
            da.append(round(100.0 * pda / total_occupied_hours, 2))
            udi.append(round(100.0 * pudi / total_occupied_hours, 2))
    return da, udi