On the basis of the assumptions explained above, this section presents the results of fire and evacuation co-simulation under each of the construction scenarios. The evacuation times associated with generated snapshots as well as the duration and cost of each construction schedule are obtained and discussed in this section. In addition, the safety factor estimation and schedule evaluation are explained. for the sake of comparing the impacts of a fire incident on the construction plans, the RSETs of both scenarios (with and without fire) were utilized as the standard evaluation metrics. Further investigations can include the ASET metric as well for the scenarios with fire.
5.1. Fire and Evacuation Co-Simulations
After running the fire and evacuation co-simulation, the evacuation time of the entire building and the 8th floor, as the scope of the investigation, were obtained for 90 snapshots with fire plus 18 snapshots without fire. Results showed that in eight snapshots with fire, the simulation did not converge, due to the blockage of the only available exit door (when the fire location is right in front of that exit and occupants are stuck). For a better illustration of the data associated with the 82 remaining with-fire cases, a ∆T parameter was defined as the ratio between evacuation times of each with-fire (RSET
With Fire) and without fire (RSET
Without Fire) scenario in the same snapshot. Based on this definition, ∆T can be in one of these three conditions: ∆T > 1, i.e., fires negatively impacting the evacuation time hence longer evacuation at the with-fire snapshots compared to without-fire ones; ∆T = 1, i.e., fire has no impact on the evacuation time; and ∆T < 1, i.e., the fire incident positively impacting the evacuation time, hence making the evacuation shorter.
Figure 5 shows the distribution of the number of snapshots in each ∆T category in the entire building, as well as the 8th floor.
The ∆T > 1 category would be originally the expected case to approves the hypothesis of a longer time being required for the evacuation under the fire emergency. It is expected that, by implementing the fire, the agents’ visibility reduces and their speed would decline [
30]. As a result, agents exposed to the smoke will have a slower evacuation and the whole evacuation time should be increased accordingly. Since the ‘without fire’ scenario’s evacuation time is considered as the safe baseline for each snapshot, the closer the evacuation times to the baseline (∆T to 1) the safer the scenario will be from the viewpoint of the fire evacuation. Therefore, ∆T can be considered as an indicator of the unsafe condition. The variation of ∆T values for both 8th floor and the whole building are illustrated in
Figure 6.
The two remaining categories, i.e., the same and faster evacuation time under the fire compared to the normal condition, appear contradictory to the hypothesis of longer evacuation in the presence of fire. Although the value of ∆T for most cases of ∆T < 1 is very close to 1.0 (between 0.87 and 0.998, with an average at 0.96 for the entire building and between 0.94 and 0.99, with an average of 0.98 for the 8th floor), and could have been attributed to the stochastic speed settings for the occupants; further investigations were made in the associated cases for both fire zone and the entire building, to identify the root cause. It was observed that three triggers control the evacuation time, which are critical egress, i.e., width of the door or staircase; the renovation operation location; and the fire location. For the cases where ∆T = 1 (i.e., the fire did not affect the evacuation time), it was concluded that the width of the door (Exit Door-2) and the staircase (Staircase-02) were the main influencers behind that. The exit door and staircase were narrower than the others, and thus occupants were more congested and evacuated more slowly through them. Accordingly, in these cases, the fire location was at the spiral stairs or the Exit Door-4 (please refer to
Figure 3), the evacuation of these narrower locations was not impacted and thus the time was not changed. This phenomenon occurred in 22% and 8% of the cases related to the evacuation of the eighth floor and the entire building, respectively. To help support this analysis,
Figure 7 shows the visibility range in different areas of the floor during the evacuation for different time steps, when the fire is in front of Exit Door-04 (
Figure 7a) or the spiral Stair (
Figure 7b). In all captured heatmaps, the red color means the horizontal visibility of 3.0 m, i.e., the full visibility range of an agent, and dark blue means the visibility of zero. In both figures, the area of congestion during the evacuation is marked. For all snapshots in the category of ∆T = 1, according to the fire slices results (
Figure 7), propagation of smoke shows that in both mentioned fires, during the entire evacuation process, smoke does not reach the congestion area, hence will not impact the occupants’ speeds. All the evacuations are completed before 184 sec, which is when the visibility decreases in the congestion area. That is why the evacuation duration remains the same as in the case of not having any fire. Of course, the panic effect was not included in this study, and that can be considered one limitation of the simulation.
For the cases of ∆T < 1 (i.e., faster evacuation than without fire), the analysis showed different causes for the 8th floor and entire building. The agents are modeled to opt for the closest exit (shortest path). In the case of whole-building analysis, when narrower staircases are blocked due to construction or fire, the occupants select the other (wider) exits, which speeds up the overall evacuation. This indicated that the width of stairs was a more influential factor than the fire, when the fire was placed in the E-block of the building. In the case of 8th floor, by the time the occupants get into the congestion, there is yet no smoke there, because the source of fire is far from the critical exit door. Hence, the with and without fire evacuation scenarios are almost the same, with a minor difference in the number of people in the congestion (∆T in most of these cases is slightly less than 1). Another interesting observation was when the smoke completely reaches the congestion, at t = 57 sec as shown in
Figure 8. Checking the frame-by-frame evacuation behavior, in this case, suggests that the speed reduction due to the smoke has helped the congestion by moving the crowd more smoothly and letting the occupants (accumulated from the 8th floor and floors above) evacuate in less density and a shorter time. However, this behavior of the agent-based model may not be realistic and requires further investigation to be validated and justified. It is evident for the fire and evacuation community that due to the lack of collected behavioral datasets, the validation process which should be a primary step in evacuation modeling becomes very challenging. There is a general absence of uniformity for evaluating evacuation models, which is mostly subjective depending on the users’ acceptance. However, the International Standards Organization (ISO) published initially the ISO 20414:2020 standard named “Verification and validation protocol for building fire evacuation models” through the committee of FSE (ISO/TC 92/SC 4) which should help in the future with this issue [
48,
49].
5.2. Construction Planning Results (Time, Cost and Safety)
Three criteria were considered for the schedule evaluation, i.e., cost; time; and safety, each of which has a separate measure and evaluation process. While several different metrics can be used to measure each criterion, and also there would be various ways to combine the metrics, the multicriteria decision-making aspect is not within the scope of this paper, and future studies will be needed in this regard. Nevertheless, we test some possibilities, mainly to show how the proposed method of this paper can be used in quantitative decision analysis.
Table 7 compares the three schedule types in terms of the metrics used for cost and time. Starting with the time criterion, the final 32 schedules had a total project duration of 65 days for schedule type 1 (1 crew–1 crew) and 36 days for both type 2 and 3 alternatives (i.e., 2 crews–1 crew and 2 crews–2 crews). Given the high demand for access to the labs, it was assumed, as an owner requirement, that the desired project duration is 60 days or less. Moving on to the cost metric, considering that the only cost component modeled in the project is the direct cost (leading to the independence of project budget from duration), the authors used the ‘mode cost’, i.e., cash flow analysis, as a metric for cost evaluation. The mode cost was considered as the ‘highest most frequent daily cost’ throughout a schedule. Schedule type 1 with
$60, schedule type 3 with
$120, and schedule type 2 with
$160, had the lowest to highest mode costs, respectively. Following the regular logic of construction projects, it was assumed that any contractor would prefer the lowest mode cost. Accordingly, the assumed metric for the cost was having a mode value less than
$160.
For the third criterion, i.e., safety, two parameters were considered: fatalities and evacuation time. The occupants’ lives were considered paramount; therefore any schedule that had a snapshot with fatality was eliminated. Twelve schedules consisted of eight snapshots with fatalities, all of which were considered unacceptable. On the other hand, for the evacuation times (the remaining 82 snapshots), those with ∆T < 1 were assumed to be safe (until further studies will make better clarifications regarding the anomalies explained earlier). After that, two different methods were applied to provide a metric for the safety of each schedule.
Method A, Average total evacuation extension—In the first method, five ∆Ts per snapshot (one ∆T for every fire incident modeled), as expressed by Equation (1), are averaged to obtain the total evacuation time extension due to the fire, as
in Equation (2).
where
is the required egress time in fire scenario and
is required egress time in the baseline scenario (without fire).
where
is calculated from Equation (1) and
i is associated with the fire scenarios (from 1 to 5 in our case study with considering five different fire incidents). The
for each snapshot is then multiplied by the duration of that snapshot to calculate the fire risk of the snapshot (
) as shown by Equation (3).
where
is the summation of all
related to one specific snapshot with five different fire scenarios and
is the time of each snapshot occurrence in the schedule. It must be noticed that we are assuming similar likelihoods for the incident of all five different fires.
In each schedule, all values of
are summed up, and the
is eventually calculated as the factor of safety for the alternative schedule, as shown by Equation (4).
where
is derived from Equation (3) and
i is the number of different snapshots existing/ considered in one schedule.
Method B, Maximum evacuation extension—In the second method, the maximum of five different ∆Ts is selected as
, as suggested by Equation (5). In each schedule, the maximum
for different snapshots of the schedule is considered as the safety factor (
) as shown by Equation (6).
where
is the difference between evacuation time in with and without fire scenarios and
i is associated with the fire scenario (from 1 to 5 for our case study).
where
comes from Equation (5) and
i is associated with the fire scenarios (from 1 to 5 in our case study).
Both methods A and B were applied to both fire zone and the entire building evacuation times separately, for our case study. The results were combined with the other two criteria, i.e., time and cost, to compare the construction schedule alternatives, as will be explained in the next section.
5.3. Construction Planning Evaluation
It is worth noting again that developing a multi-criteria schedule evaluation has been beyond the scope of this paper. The main objective of this study was to introduce co-simulation as a tool for quantitative analysis of safety in construction renovation projects. However, to show how the result can be mixed with other traditionally used metrics to evaluate schedules, i.e., time and cost, one simple combination option, i.e., weighted summation, is tested here. Further studies are necessary to determine the optimal weights, as well as better options for the metrics and multi-criteria function.
After calculating the three performance metrics for the assessment of the schedule, we define ‘TCS’ as an index to evaluate the goodness of construction scenarios. TCS in this study was defined as the weighted sum of the range-normalized values for these three metrics, for each alternative schedule, as shown by Equation (7). The aim is to identify the schedule with the lowest TCS. We used both versions of FR and FR’ and calculated the TCS twice for the case study. Detailed results are shown in
Appendix A.
Table 8 shows the results of Method (A) for the eighth floor.
where
is the normalized total duration of the schedule completion,
is the normalized most repetitive daily cost (mode cost) in the associated schedule,
is the normalized safety metric explained earlier, and the weights w
1 through w
3 are user-defined (for the sake of this paper, we applied equal weights to each metric).
As shown in
Table 9, any schedule that had a snapshot with fatality was eliminated from the selection process, leaving schedules 16, 18, 20, 22, and 24 to decide between. Logically, schedule 16 provided the lowest TCS for both the 8th floor fire zone and the whole building; however, it violated the core owner requirement of a duration less than 60 days and hence was also excluded from the comparison. Consequently, this table shows the TCS values for each evaluation method for the 17 schedules and the corresponding best schedule.