Construction Planning and Scheduling of a Renovation Project Using BIM-Based Multi-Objective Genetic Algorithm

Featured Application: Building information modeling with multi-objective genetic algorithm. Abstract: Renovation is known to be a complicated type of construction project and prone to errors compared to new constructions. The need to carry out renovation work while keeping normal business activities running, coupled with strict governmental building renovation regulations, presents an important challenge affecting construction performance. Given the current availability of robust hardware and software, building information modeling (BIM) and optimization tools have become essential tools in improving construction planning, scheduling, and resource management. This study explored opportunities to develop a multi-objective genetic algorithm (MOGA) on existing BIM. The data were retrieved from a renovation project over the 2018–2020 period. Direct and indirect project costs, actual schedule, and resource usage were tracked and retrieved to create a BIM-based MOGA model. After 500 generations, optimal results were provided as a Pareto front with 70 combinations among total cost, time usage, and resource allocation. The BIM-MOGA can be used as an efﬁcient tool for construction planning and scheduling using a combination of existing BIM along with MOGA into professional practices. This approach would help improve decision-making during the construction process based on the Pareto front data provided.


Introduction
Nowadays, construction management is facing several uncontrollable and controllable challenges. Unpredictable weather, local resource scarcity, and a range of strict regulations are common uncontrollable factors in the construction industry [1]. Traditional construction management also deals with controllable factors such as inconsistent design, lack of constructability, inaccurate documentation, and redundant work processes [2]. Building information modeling (BIM) has been widely used for a decade to improve the efficiency of a construction project, especially in construction planning and documentation. It provides data and visualization to project managers to make better decisions for improved quality and faster processes [3][4][5]. Through clash detection capability, BIM also helps resolve design conflict resulting from the complex interface between various building elements which commonly exist in large buildings and renovation work [6]. Construction complexity has a direct impact on project communication and performance [7]. Renovation projects are particularly complex because of physical constraints, such as obsolete conditions, limited access, and unknown conditions. Some renovation projects become more complex if there is a requirement for the business to remain open during the renovation [8]. The contractors are required to take extra care in planning in order to ensure that the construction activities do not interfere with on-going business.
The apparent benefit of BIM in dealing with the technical complexity of a construction project has seen a widespread adoption of BIM technology as a tool to improve efficiency in construction management. Recently, an artificial intelligence optimization approach has been suggested to further enhance BIM capability. Such an advancement represents a research opportunity to develop and apply a BIM-based optimization approach to existing traditional construction approaches [9].

Building Information Modeling
Building information modeling (BIM) was introduced in the early 2000s as an information model of building elements [10]. It has been widely used in the construction industry since then because it enables better project information flow throughout a building life cycle, from pre-construction to construction, post-construction, and operation [11]. It is used for the improvement of planning and design, clash detection, visualization, cost planning, and data management. A BIM model consists of (1) physical properties in three dimensions, such as material, density, weight, or location, and (2) embedded information, such as construction specifications or repair manuals [12]. Several researchers and practitioners have used BIM in a variety of aspects [13]. Most notably, it is used to improve communication between the project stakeholder over a traditional non-BIM approach [14] and to improve information synchronization of building coordination in a three-dimensional space. This is a significant improvement on traditional construction drawings and is very useful for complex projects, such as extra-large buildings that require synchronization between multi-disciplinary parties as well as multi-source documents [15].
BIM has been used for existing buildings in either renovation projects or for the improvement of existing facilities. Facility management [16] and sustainability improvement for existing buildings are the major benefits of BIM [17]. It also significantly improves construction processes when working on existing buildings; however, several challenges need to be addressed, including the difficulty of stakeholder involvement and the lack of automated data capturing [18]. In converting existing buildings to BIM, laser scanning surveyance and total station techniques are suggested [19]. This technique is also used for collecting historical building information [20]. It can be seen that using BIM as a tool for complex renovation projects is viable for making the project more successful.

Construction Planning and Scheduling Optimization
Construction planning and scheduling are two of the most important activities in all phases of construction. A well-planned and optimized construction schedule contributes significantly to a successful project. Several researchers adopted optimization techniques to solve construction planning and scheduling problems [21][22][23]. In the past, linear programming was one of the techniques most widely used to solve this issue. Multi-objective linear programming is suggested for optimization when faced with constrained resource problems, including construction costs, project duration, resource idle time, and project delivery time [24]. A particle swarm optimization is also suggested for construction scheduling problems in order to solve resource constraint issues [25]. An ant colony optimization, which focuses on minimizing project duration while varying activity sequences, is another approach for resource-constrained projects [26]. Genetic algorithms (GA) have been used in several construction scheduling research problems, such as resource allocation and leveling optimization [27], concrete precast production optimization [28], nonproductive resource determination [29], and building performance assessment for long time usages [8]. The basic concept and overview of the applications of optimization algorithm in construction Appl. Sci. 2021, 11, 4716 3 of 18 processes can be further found in References [9,30,31], where advantages and limitations were also discussed.
In real construction projects, there are always multiple objectives forming a set of solutions as optimization data for project managers to make decisions. Consequently, the field of multi-objective optimization has been studied intensively [32,33]. Multi-objective genetic algorithms (MOGA) were first introduced by Tadahiko Murata in 1995 to provide decision-makers with Pareto optimization instead of constant weights [34]. It has been discussed in several construction optimization studies. For example, it is used for construction time-cost trade-off optimization [35] by defining optimal time and cost with variable crane locations and crews by minimizing crane stoppage time [36]. It is also used to improve safety by optimizing safety costs on different project layouts [37]. Even construction quality is used as an objective function for MOGA along with construction project time and cost, with the results showing their weights and Pareto solutions for the repetitive scheduling method [38]. The optimization results can be used for the decision-maker in other applications, such as planning, scheduling, or even unpredictable situations such as construction logistics network during a hazard [39].
A combination of BIM and optimization techniques has been discussed in the past [30], such as genetic algorithms (GA), particle swarm optimization (PSO), and fuzzy theory working with Autodesk Revit as a BIM tool [9] used to improve information flow and visualization. An optimization technique such as a genetic algorithm was used with BIM for possibilities usage of information from BIM directly and parsing the data to GA directly [40]. More than the fuzzy-logic dataset was possible usages in GA for construction scheduling [41]. Other applications of BIM-based optimization were used for the construction assembly line [42], building energy performance evaluation [43] and planning of construction site layout [44].
Regarding such characteristics of renovation projects, more flexible planning models unique to the extra-large-scale renovation projects are needed in order to improve the planning efficiency in terms of actual construction time, total cost, and worker utilization fluctuation [27,[45][46][47], especially in real construction projects where there are always multiple objectives forming a set of solutions for project managers to make decisions. Therefore, with the performance of BIM, this paper explored opportunities to develop a systematic BIM-based multi-objective approach for construction planning and scheduling of renovation projects.
In this paper, the proposed BIM-based MOGA model is discussed. A two-year building renovation project is used as a case study. The project construction planning and scheduling data are optimized and building information modeling data are combined and optimized. To achieve this aim, the paper is organized as follows. In the next section, the related works and research methodology are discussed. The BIM-MOGA process is described along with an experimental design based on a case study. This is followed by discussion of the possibility of applying the proposed model to complex construction projects. The project started on 17 October 2018, with a construction duration of 720 days. The proposed schedule was divided into four phases: (1) the 4th floor renovation from day 1 to 180, (2) the 3rd floor renovation from day 181 to 360, (3) the 2nd floor renovation from day 361 to 540, and (4) the 1st floor renovation and landscape from day 541 to 720.

Project Background
The project was finally completed on 2 November 2020, after an extension period of 36 days. The project scopes were to improve the current look and atmosphere of the building, update MEP (mechanical, electrical, and plumbing) and ICT systems to modern standards, and to revamp the landscape and surrounding area with a new concrete patio.
Renovations are often carried out in an existing building that remains in operation. Unique challenges for the renovations are, for instance, an additional level of uncertainty (unexpected problems that suddenly arise when the demolition phase has started, the unavailability or inaccuracy of existing as-built drawings) and particular management issues (e.g., difficulty in managing the interaction between the construction team and users due to limited working space) [8,48]. The project site of Chiang Mai University Library during the renovation is displayed in Figure 1. The project started on 17 October 2018, with a construction duration of 720 days. The proposed schedule was divided into four phases: (1) the 4th floor renovation from day 1 to 180, (2) the 3rd floor renovation from day 181 to 360, (3) the 2nd floor renovation from day 361 to 540, and (4) the 1st floor renovation and landscape from day 541 to 720.
The project was finally completed on 2 November 2020, after an extension period of 36 days. The project scopes were to improve the current look and atmosphere of the building, update MEP (mechanical, electrical, and plumbing) and ICT systems to modern standards, and to revamp the landscape and surrounding area with a new concrete patio.
Renovations are often carried out in an existing building that remains in operation. Unique challenges for the renovations are, for instance, an additional level of uncertainty (unexpected problems that suddenly arise when the demolition phase has started, the unavailability or inaccuracy of existing as-built drawings) and particular management issues (e.g., difficulty in managing the interaction between the construction team and users due to limited working space) [8,48]. The project site of Chiang Mai University Library during the renovation is displayed in Figure 1.

Construction Scheduling
Before the project commenced, the project schedule was initially set according to the project contract. Information was extracted from the contract, which consisted of a project delivery system, bill of quantities (BOQ), resources, phasing, and payment. Then, during the construction, the actual schedule and resources usage were retrieved on-site from daily site progress reports. The first phase focused on the demolition of existing walls and removing existing ceilings and floors. Then, the project closed each floor, starting with the 4th floor, while the 1st, 2nd, and 3rd floors remained operational as usual.
The original scheduling data were retrieved from the main contractor in the form of Microsoft Project and Microsoft Excel files. Over two years, on-site data were collected, including the number of workers, crew sizes, actual durations, and project progress. Then, all data from different sources were used to create the BIM.

Construction Scheduling
Before the project commenced, the project schedule was initially set according to the project contract. Information was extracted from the contract, which consisted of a project delivery system, bill of quantities (BOQ), resources, phasing, and payment. Then, during the construction, the actual schedule and resources usage were retrieved on-site from daily site progress reports. The first phase focused on the demolition of existing walls and removing existing ceilings and floors. Then, the project closed each floor, starting with the 4th floor, while the 1st, 2nd, and 3rd floors remained operational as usual.
The original scheduling data were retrieved from the main contractor in the form of Microsoft Project and Microsoft Excel files. Over two years, on-site data were collected, including the number of workers, crew sizes, actual durations, and project progress. Then, all data from different sources were used to create the BIM.

BIM Model Development
Over 50 years, there had been several minor renovations and modifications which were not recorded in any hardcopy format. A BIM was developed based on the as-is condition from the original as-built drawings, terrestrial laser scanning survey data, and on-site survey data.
The BIM model of the renovation project was initially created based on the as-built drawings; however, several details were missing either due to faded blueprints or nonrecorded renovation activities. Moreover, information from a traditional two-dimensional blueprint was obscure and prone to error due to several unconnected figures, tables, and remarks. Then, a traditional on-site survey along with a terrestrial laser scan was carried out in order to collect the as-is condition of the building and surrounding area. Three-dimensional point cloud data were retrieved from the survey. Then, the point cloud data were laid on top of the existing as-built drawings to complete the building information. When data from different sources were combined, discrepancy data appeared, as illustrated in Figure 2, which shows an example of an opening position error when comparing the asbuilt drawings and as-is condition. Today, this approach is a typical initial BIM procedure for gathering data in renovation projects.

BIM Model Development
Over 50 years, there had been several minor renovations and modifications which were not recorded in any hardcopy format. A BIM was developed based on the as-is condition from the original as-built drawings, terrestrial laser scanning survey data, and onsite survey data.
The BIM model of the renovation project was initially created based on the as-built drawings; however, several details were missing either due to faded blueprints or nonrecorded renovation activities. Moreover, information from a traditional two-dimensional blueprint was obscure and prone to error due to several unconnected figures, tables, and remarks. Then, a traditional on-site survey along with a terrestrial laser scan was carried out in order to collect the as-is condition of the building and surrounding area. Threedimensional point cloud data were retrieved from the survey. Then, the point cloud data were laid on top of the existing as-built drawings to complete the building information. When data from different sources were combined, discrepancy data appeared, as illustrated in Figure 2, which shows an example of an opening position error when comparing the as-built drawings and as-is condition. Today, this approach is a typical initial BIM procedure for gathering data in renovation projects. Following the above procedure, the BIM model was created for visualization of the project (3D BIM) and the project schedule (4D BIM). Figure 3 displays the 3D BIM model for visualization. It displayed the final appearance when the renovation was completed. This was used for communication between project stakeholders. The scheduling and cost data model are displayed in the BIM model from a combination of the 3D BIM and a project Gantt chart, as illustrated in Figure 4. The appearance of building elements was changed and highlighted according to a specific date. This allowed the project stakeholders to discuss and monitor the current project schedule in an efficient manner. Moreover, the 4D BIM renovation was helpful since several demolition tasks were not displayed in the model but required significant project resources. This allowed the project stakeholders to plan and schedule for non-element construction tasks much more easily than in a traditional approach. In this 4D BIM, a demolition model and a new renovation model were overlapped for planning and scheduling in different project phases.
During construction, schedule and resource information was collected on site and put back into the model for control and monitoring. Finally, total cost data were then used for project management throughout the project life cycle. Following the above procedure, the BIM model was created for visualization of the project (3D BIM) and the project schedule (4D BIM). Figure 3 displays the 3D BIM model for visualization. It displayed the final appearance when the renovation was completed. This was used for communication between project stakeholders. The scheduling and cost data model are displayed in the BIM model from a combination of the 3D BIM and a project Gantt chart, as illustrated in Figure 4. The appearance of building elements was changed and highlighted according to a specific date. This allowed the project stakeholders to discuss and monitor the current project schedule in an efficient manner. Moreover, the 4D BIM renovation was helpful since several demolition tasks were not displayed in the model but required significant project resources. This allowed the project stakeholders to plan and schedule for non-element construction tasks much more easily than in a traditional approach. In this 4D BIM, a demolition model and a new renovation model were overlapped for planning and scheduling in different project phases.
During construction, schedule and resource information was collected on site and put back into the model for control and monitoring. Finally, total cost data were then used for project management throughout the project life cycle. Sci. 2021, 11, x FOR PEER REVIEW 6 of 19

Optimization Definitions
The optimization problem for the multi-objective optimization was defined with the goal of minimizing total cost, construction time, and variance in the number of workers per day. The relevant parameters and associated formula are presented in Equations (1)-(10).

Problem Definitions
Optimizing project cost, duration, and resource usage while maintaining the BIM data collection is crucial for this research. Typically data in the BIM model itself are interconnected to the 3D model and schedule to enable 4D modeling with an automatic update of scheduling information. However, when dealing with a complex solution, BIM requires external optimization tools and needs to leave the BIM data pipeline. This causes redundant manual tasks which are prone to error. The goal was to optimize those objectives and exchange data back and forth via BIM. To satisfy these requirements, the BIM model can be used to embed information into a three-dimensional representation [40,49] and then connect to the multi-objective genetic algorithm (MOGA) to solve multiobjective problems. At the end, optimal Pareto front is provided for decision-makers instead of static objective weights for cost, duration, and resources.

Optimization Definitions
The optimization problem for the multi-objective optimization was defined with the goal of minimizing total cost, construction time, and variance in the number of workers per day. The relevant parameters and associated formula are presented in Equations (1)-(10).

Problem Definitions
Optimizing project cost, duration, and resource usage while maintaining the BIM data collection is crucial for this research. Typically data in the BIM model itself are interconnected to the 3D model and schedule to enable 4D modeling with an automatic update of scheduling information. However, when dealing with a complex solution, BIM requires external optimization tools and needs to leave the BIM data pipeline. This causes redundant manual tasks which are prone to error. The goal was to optimize those objectives and exchange data back and forth via BIM. To satisfy these requirements, the BIM model can be used to embed information into a three-dimensional representation [40,49] and then connect to the multi-objective genetic algorithm (MOGA) to solve multiobjective problems. At the end, optimal Pareto front is provided for decision-makers instead of static objective weights for cost, duration, and resources.

Optimization Definitions
The optimization problem for the multi-objective optimization was defined with the goal of minimizing total cost, construction time, and variance in the number of workers per day. The relevant parameters and associated formula are presented in Equations (1)-(10).

Problem Definitions
Optimizing project cost, duration, and resource usage while maintaining the BIM data collection is crucial for this research. Typically data in the BIM model itself are interconnected to the 3D model and schedule to enable 4D modeling with an automatic update of scheduling information. However, when dealing with a complex solution, BIM requires external optimization tools and needs to leave the BIM data pipeline. This causes redundant manual tasks which are prone to error. The goal was to optimize those objectives and exchange data back and forth via BIM. To satisfy these requirements, the BIM model can be used to embed information into a three-dimensional representation [40,49] and then connect to the multi-objective genetic algorithm (MOGA) to solve multi-objective problems. At the end, optimal Pareto front is provided for decision-makers instead of static objective weights for cost, duration, and resources.

BIM-MOGA Model
In this section, the development of BIM-MOGA is described. After retrieving project information and constraints from the BIM model, which is prepared in Excel spreadsheet format, the spreadsheet is created by formulating actual construction time, resource utilization fluctuation, and total cost calculations, as shown in the following subsections. It is then transferred to the MOGA module. These two modules are connected by Activity ID. A Python-based modeling tool was written to automate the MOGA process. The MOGA process solves multi-objective optimization problems using GA. In doing so, a Pareto front must firstly be defined. Once the Pareto front has been defined, the best solution will then be chosen by the decision-makers. The renovation scheduling was generated by Microsoft Project. After that, the researchers exported the file to Autodesk Navisworks and exported the BIM model from Autodesk Revit to Autodesk Navisworks in order to create a 4D BIM for the renovation project. This allowed the concerned parties to monitor the process and reduce workloads and conflicts while working at the site. It also helped improve collaboration between the concerned parties. Each related part is explained in detail in the following sections.
The proposed model is a multi-objective optimization for minimizing project duration, cost, and resources simultaneously. The steps to generate an optimized schedule are illustrated in Figure 5.

BIM-MOGA Model
In this section, the development of BIM-MOGA is described. After retrieving project information and constraints from the BIM model, which is prepared in Excel spreadsheet format, the spreadsheet is created by formulating actual construction time, resource utilization fluctuation, and total cost calculations, as shown in the following subsections. It is then transferred to the MOGA module. These two modules are connected by Activity ID. A Python-based modeling tool was written to automate the MOGA process. The MOGA process solves multi-objective optimization problems using GA. In doing so, a Pareto front must firstly be defined. Once the Pareto front has been defined, the best solution will then be chosen by the decision-makers. The renovation scheduling was generated by Microsoft Project. After that, the researchers exported the file to Autodesk Navisworks and exported the BIM model from Autodesk Revit to Autodesk Navisworks in order to create a 4D BIM for the renovation project. This allowed the concerned parties to monitor the process and reduce workloads and conflicts while working at the site. It also helped improve collaboration between the concerned parties. Each related part is explained in detail in the following sections.
The proposed model is a multi-objective optimization for minimizing project duration, cost, and resources simultaneously. The steps to generate an optimized schedule are illustrated in Figure 5.

Optimization Model
The optimization model was developed and organized into three main modules: (1) Initialization Module: Construction data are calculated, i.e., an initial project schedule and the number of total float days, (2) BIM Module, and (3) MOGA Module. These three main modules are described in more detail, as follows.

Optimization Model
The optimization model was developed and organized into three main modules: (1) Initialization Module: Construction data are calculated, i.e., an initial project schedule and the number of total float days, (2) BIM Module, and (3) MOGA Module. These three main modules are described in more detail, as follows.

Initialization Module
There are three main calculations for this module, as listed in the following subsections.

Actual Construction Time (Ta) Calculations
The calculation used construction data from each activity and project calendar constraint, i.e., activity duration (D i ), activity sequence, and exception date (e i ) from the baseline schedule. The precedence networks were created using the Precedence Diagram Method (PDM). In the PDM, there are four precedence relationship constraints: Finish to Start (FS), Start to Start (SS), Finish to Finish (FF), and Start to Finish (SF), with overlapping time (L i ). The PDM calcuation procedure used sequential forward and backward calculations via the network to calculate the early start times (ST i ) to early finish times (FT i ), and late start times (LS i ) to late finish times (LF i ) for each activity. From these values, actual construction time (Ta) and total float time (TF i ) were determined using Equations (1) and (2):

Resource Utilization Fluctuation (Mx) Calculations
The resource utilization fluctuation, also known as the static moment, Mx, was calculated by summing the total numbers of daily fluctuations in Equations (3) and (4): [Number o f workers per day t ] 2 (3)

Number o f workers per day t
where n is the total number of activities in day t, Q i is the quantity of activity i, Di is the duration of activity i, PDRi is the productivity rate of the worker for activity i, and T is the total number of project working days.

Total Cost (Ct) Calculations
Total cost (Ct) was calculated using Equation (5), which is referred to as direct cost (DC), indirect cost (IC), and late penalty fee (LPF). DC was calculated using Equation (6), which includes material cost and labor cost, and related costs including freight and transportation. Material cost was a function of material quantity from the BIM model, as shown in Equation (9), and labor cost was the function of daily payrate of workers, as shown in Equation (10). Indirect costs were majorly calculated as overhead cost based on working time with a project delay causing surplus indirect cost. IC was calculated using Equation (7) and LPF was calculated using Equation (8).
LPF = Daily penalty fee × (Ta − Tc) (8) where Tc is the contract construction time.
[Number o f workers per day t x Pay rate t ]

BIM Module
Once the construction data were calculated, they were transferred to the BIM module. BIM enabled the concerned parties to see the 3D visualization of the project in its current stage. A collision check was conducted before the start of the construction to identify any visible clash between building elements. The schedule planning was correlated with the model. The BIM-based schedule was integrated to the 4D model and used as a visualization tool. The Activity ID generated by the 4D BIM was used to link schedules with 3D objects. To create the actual construction stimulation, actual work start time and end time were fed into the model. Resources and cost were also monitored. The "BIM Model" was used during the process as a dynamic management process simulation of project progress. A baseline schedule providing required BIM data was subsequently imported to the optimization module.

MOGA Module
The MOGA module starts by defining optimization objectives, constraints, and project parameters. The parameters which were retrieved from construction data were project activities, activity predecessors, successors, resource availability, cost data, and project calendars. Then, a genetic algorithm optimizer was employed. Once the MOGA process is complete, several optimal solutions are obtained as a Pareto front. The data is then used in BIM for scheduling, budgeting, and planning for different optimal alternatives. This allows a decision-maker to select an optimum solution from a set of solutions based on construction project constraints.

Optimization Model
Considering the multi-objective optimization problem with different scales and units, the nondominated sorting genetic algorithm 2 (NSGA-II) [50] was adopted in order to minimize the three objectives listed in Equations (11)- (13). Furthermore, because of deficiencies with respect to depicting the optimal Pareto front, the weighted-sum genetic algorithm approach [51] was not an option for the multi-objective optimization problem. Several optimal solutions for multi-objective problems were obtained with the following objectives and constraints.
Objectives: Minimize f1(s, x): T a (11) Minimize f2(s, x): Ct (12) Minimize f3(s, x): Mx (13) where Ta is the actual construction time, Ct is the total cost, Mx is the worker utilization fluctuation moment, s is the shifting time set, and x is the predecessor option set. Constraints: where r xt is the worker demand type-x on day t, MWA x is the maximum worker availability of resource type x, and LPF is the late penalty fee. This study considers four types of resource, namely general workers, skilled carpenters, sanitary plumbers, and electricians. LPF was calculated from the construction contract, and was not allowed to be greater than 10% of the total contract cost, otherwise the contract will be terminated.

Decision Variables
Each solution comprises an N number of blocks, as displayed in Figure 6. Each block represents individual problem decision variables providing two different values. The first value is the number of shifting times from the original plan, while the second value is a predecessor activity option. 0 ≤ S i ≤ TF i (17) 0 ≤ x ihg ≤ 1 (18) where S i is the shifting time of activity i, TF i is the total float of activity i, and X ihg is the predecessor option between h or g of activity i.
contract, and was not allowed to be greater than 10% of the total contract cost, otherwise the contract will be terminated.

Decision Variables
Each solution comprises an N number of blocks, as displayed in Figure 6. Each block represents individual problem decision variables providing two different values. The first value is the number of shifting times from the original plan, while the second value is a predecessor activity option.
where Si is the shifting time of activity i, TFi is the total float of activity i, and Xihg is the predecessor option between h or g of activity i. The predecessor option is defined (Xihg) and optimal scheduling is now processed by considering both Set1 and Set2 options. The Set2 option is proposed based on the renovation project characteristics, including additional levels of uncertainties and particular management issues. If the interactive planner's viewpoint is obtained, selecting the Set2 option can provide more efficient optimal scheduling. The shifting time and predecessor option are two sets of decision variables with 250 total number of activities for all decision variables. The population size is chosen as 500, as suggested in previous research for double or quadruple design variables [52]. The predecessor option is defined (X ihg ) and optimal scheduling is now processed by considering both Set1 and Set2 options. The Set2 option is proposed based on the renovation project characteristics, including additional levels of uncertainties and particular management issues. If the interactive planner's viewpoint is obtained, selecting the Set2 option can provide more efficient optimal scheduling. The shifting time and predecessor option are two sets of decision variables with 250 total number of activities for all decision variables. The population size is chosen as 500, as suggested in previous research for double or quadruple design variables [52].

Genetic Operations
By performing GA operator, a next generation population is created based on fitness value calculations using GA, including a selection, crossover, and mutation. The selection technique depends on a uniform random mechanism. The crossover is performed using a one-point crossover routine. The mutation is performed by uniform randomization around old variable values. The number of mutations proportionally changes based on the setting of the mutation rate between the value of 0 and 1 [53]. The GA steps are listed below, with the pseudocode displayed in Figure A1 in the Appendix A.

1.
Population creation: in each solution, the shifting time set values (Set1) and predecessor option set values (Set2) are randomized. Set1 is uniformly randomized from 0 to a sigma parameter value, while Set2 is uniformly randomized from 0 to 1. From continuous adjustment, the sigma parameter value in this research was set at 10, which is within the upper and lower bounds.

2.
Parent selection: the uniform random technique is used to select two parents with equal selection probability. After the parents are selected, they breed with the onepoint crossover method.

3.
Crossover: one crossover point is selected randomly with the one-point crossover method. The initial solution values are copied from the first parent, while the other solution values are copied from the second parent. 4.
N-point mutation: this paper also proposed an n-point mutation technique to reduce computation complexity. The technique initially selects n number of decision variables in a solution by uniform randomization. Then, they are mutated by uniform randomization around old values with upper and lower bounds. The upper bound is the old value that adds the sigma value, while the lower bound is the old value that subtracts the sigma value. If the mutated values are out of bounds, then the value is uniformly randomized in possible value bounds. The parameters set, in this research, were n as 1 and possible values were between 0 and TF value. However, if the shifting time was out of bounds, the sigma value was used instead. 5.
Fitness calculation: the NSGA-II algorithm is used to calculate fitness scores, which are used to find non-dominated solutions, called a Pareto front. A challenge of the fitness score calculation in this research was the requirement for the PDM network to re-create everytime after mutation, resulting in complexity of the computation. 6.
Pareto front selection: the Pareto front selection uses the crowding distance technique to reduce the number of Pareto front which is over the upper bound. This is based on a tournament of crowding distances. 7.
Predecessor option (X ihg ): In renovation projects, construction phases are normally considered based on work spaces. This is also defined as Set2 decision variables. Figure 7 illustrates the initial predecessor (X h ), where construction sequences start with the 4th floor, followed by the 3rd, 2nd, and 1st floors, and then the basement. A suggestion of renovation sequences is displayed as an alternative (X ihg ). The renovation starts from the 4th floor, followed by the the 3rd, 2nd, and 1st floors, while simultaeneously working on the basement.

Results
A Pareto front displays optimum points of the proposed BIM-MOGA model. The data were analyzed and coded in Python with the library of Pandas, Numpy, and Math [54]. The Pareto front was calculated based on the relationship between time, cost, and resource utilization fluctuation moment (Mx), with parameters and values listed in Table   Figure 7. Initial and optional construction sequences.

Results
A Pareto front displays optimum points of the proposed BIM-MOGA model. The data were analyzed and coded in Python with the library of Pandas, Numpy, and Math [54]. The Pareto front was calculated based on the relationship between time, cost, and resource utilization fluctuation moment (Mx), with parameters and values listed in Table 1. Table 1. MOGA simulation parameters and values.

Parameter Name Values
The size of start population 500 The size of minimum population 400 The size of maximum population 500 The maximum generation 500 The crossover rate (Pc) 0.5 The mutation rate (Pm) 0.004 n 1 sigma 10 Maximum days (from construction contract) 780 days Figure 8 presents data in the Pareto front for the optimal solutions with the axes of total cost, time usage, and resource allocation. The charts in Figure 9 are separated into three different objective pairs.    There are 70 total numbers of optimal solutions, as displayed in Table 2. Due to the construction contract termination, optimal solutions for a duration over 780 days are removed. The ability to extend the duration without contract termination would result in better outcomes for both cost and resources fluctuations. Figure 10 displays optimal resource fluctuation moment over GA generations. There are 70 total numbers of optimal solutions, as displayed in Table 2. Due to the construction contract termination, optimal solutions for a duration over 780 days are removed. The ability to extend the duration without contract termination would result in better outcomes for both cost and resources fluctuations. Figure 10 displays optimal resource fluctuation moment over GA generations.  Figure 10. Relation graph between generation and optimal resource fluctuation. Figure 10. Relation graph between generation and optimal resource fluctuation.
Pareto front gives numbers of optimal plans which could be appropriate options for the project manager to make a decision. Cost, time, and resource utilization fluctuation were considered as the objective functions. These parameters need to be correlatively considered to generate the optimal scheduling. Moreover, they were also considered as constraints in the optimization process.

Discussion
This study shows the possibility of applying BIM-MOGA to complex construction projects. It provides optimization data for a project manager to make decisions based on existing data underneath the project during the project life cycle. NSGA-II is a simple MOGA algorithm which is easy to understand and implement.

Cost
Cost is one of the analyzed factors displayed in an axis of the Pareto front chart. It relates directly to project time and resource utilization fluctuation. This relationship is depicted as a renovation process where cost data is varying according to time, visualized by 4D BIM. Although indirect costs are changed based on the project time, fluctuation of resources causes both direct and indirect costs. These relationships are clearly organized by using the BIM-MOGA technique to produce optimum solutions based on minimizing total cost. This is useful for the project with limited cost, especially the govermental project under a lumpsum contract.

Time
Time generally has an inverse relationship with cost in construction projects, while resource utilization fluctuation can have a positive or negative relationship with the time.
The Pareto front displays optimized solution between these relationships. Minimum time shown in Table 2 can be reduced by either changing time constraint (i.e., predecessor option and activity relationship) or resource constraint (i.e., resource availability). This is useful for the project with limited time preventing delay liquidated damages or penalties.

Resource Utilization Fluctuation (Mx)
Resource utilization fluctuation (Mx) is an effective improvement in the process. It can have either positive or negative impacts on both cost and time factors. Pareto front displays optimum results retrieved from BIM-MOGA. This is useful for further analysis, such as human resource risk management. Table 3 displays an example of a comparison between the original plan and an optimal plan with a set of selected data from the Pareto front. The optimal plan is selected based on the original cost with the maximum time without delay in order to avoid blacklisting for a following governmental project. When considering the cost and time of both plans, which are equal, it shows that the optimized Mx objective is clearly better than the non-optimized Mx objective. This means we can reduce worker density on the construction site, leading to reduced risk-related worker factors. Furthermore, there are many optimal plans for utilizing that are better than one original plan.

Conclusions
This paper provided a systematic BIM-based multi-objective genetic algorithm approach for construction planning and scheduling of renovation projects. The results showed a potential application of a multi-objective optimization problem using BIM with MOGA. The outputs displayed several alternatives for different scenarios for further analysis. It creates opportunities for construction experts to develop more robust and efficient tools for construction planning and scheduling using a combination of existing BIM along with MOGA in professional practice. The construction process will benefit from useful information for better decision-making depending on the strategies based on the Pareto front data provided.
It was simulated for a governmental project of an extra-large building renovation between 2018 and 2020. Direct and indirect project costs, actual scheduling, and resource usage were tracked and retrieved by the researchers. After the information was modeled as a BIM-based MOGA, the optimization results were displayed as a Pareto front with 70 combinations among total cost, time usage, and resource allocation. The increase in generation provided better results, as shown in Figure 10, which this research used at 500 generations. BIM and MOGA worked together through standard spreadsheet files, such as .xls, .xlsx, and .csv extensions. The standard file formats are easy for reading, editing, and sharing, which is beneficial when the model's results need to be sent and the receiver only has standard spreadsheet programs such as Microsoft Excel, Google Sheets, Smartsheet, or Numbers.