Research Article  Open Access
Guofeng Ma, Keke Hao, Yu Xiao, Tiancheng Zhu, "Critical Chain Design Structure Matrix Method for Construction Project Scheduling under Rework Scenarios", Mathematical Problems in Engineering, vol. 2019, Article ID 1595628, 14 pages, 2019. https://doi.org/10.1155/2019/1595628
Critical Chain Design Structure Matrix Method for Construction Project Scheduling under Rework Scenarios
Abstract
Rework risks have been a major challenge in the construction industry that constantly affects project schedules and threatens ontime project completion. Traditional project scheduling methods are not capable of modeling rework relationships between activities and mitigating the impact of resulting uncertainties during the development of project schedules. To address this challenge, a critical chain design structure matrix (CCDSM) method is proposed in this paper. The CCDSM method aims to develop construction project schedules that are adaptive to rework scenarios and robust against rework risks. The CCDSM method models and displays largescale rework relationships among activities and introduces a new rework buffer to quantitatively represent the impact of rework instances in project schedules. A maxplus algorithm is adopted in CCDSM to transform complex logic relationships into simple matrix operations, reducing computational load of schedule generation. A case study was conducted to demonstrate the implementation of the CCDSM method and assess its effectiveness in addressing rework risks. The results showed that the CCDSM is a promising tool to generate schedules, which could improve ontime project completion rate and reduce impacts of varying rework scenarios on project execution.
1. Introduction
Rework has been regarded as one of the major challenges that can adversely affect project performance in the construction industry [1]. Hwang et al. [2] surveyed about ClientRelated Rework (CRR) in 381 projects conducted by 51 construction companies and concluded that more than 80% companies and 59% projects experienced CRR, which increased project costs by 7.1% and caused 3.3 weeksâ€™ delay on average. Simpeh et al. [3] investigated 78 construction professionals and found that if the mean rework cost of a project reached 5.12%, the likelihood of the project exceeding its budget would be as high as 76%.
Extensive efforts have been made to manage reworkcaused project delay and minimize the rework impact in construction projects. When used to predict the project duration as well as its ontime completion probability, traditional scheduling techniques, such as CPM, PERT, and Gantt Charts, cannot quantitatively measure activity rework time and its impact on project completion time [4] and hence are not able to directly take into account the rework risks during project planning. Although buffers in the CCPM can aggregate uncertainties in project execution, which are reflected as blocks of resource or time redundancy in project schedules, it is still difficult to model largescale rework relationships between activities using the CCPM [5]. To further address this challenge, a design structure matrix (DSM) method has been introduced, which is designed to represent the dependency or information flow between activities, providing an effective representation of rework instances [5]. However, one limitation of the DSM method exists that the matrix format is too concise to fully display certain attributes of activities, which leads the DSM method usually to be adapted as a project duration calculator rather than a process control tool [6]. Overall, most of the current tools for managing rework are reactive, can barely compute and diagram rework instances with satisfactory accuracy, and focus on measuring the rework risk of activity duration rather than that of overall project completion time. Due to these limitations, although in many cases project managers expect rework to occur, they lack reliable tools to fully assess and address these risks at the planning stage of the projects and to mitigate potential impacts of rework with proactive measures.
Accordingly, a new method, termed critical chain design structure matrix (CCDSM), is proposed in this paper. Built on the strength of the CCPM and DSM in project scheduling, the CCDSM aims to integrate these two methods to overcome their respective limitations in developing construction project schedules that are adaptive to rework instances. The CCDSM is designed in such a way that it not only models largescale rework relationships among activities, improving the accuracy and effectiveness of project management, but also displays various attributes of activities, enabling analysis and visualization of the project schedule under rework scenarios. Moreover, the CCDSM adopts maxplus algorithm, a straightforward linear algebra algorithm, to transform complex logic relationships into simple matrix operations, which effectively reduces the computational load of generating schedules. The detailed design of the CCDSM and a case study that demonstrates the effectiveness of the CCDSM in addressing rework risks in construction projects are presented in the remainder of this paper.
2. Review of Relevant Works
2.1. Buffer Management and Critical Chain Scheduling in the CCPM
The CCPM, which has proved to be a powerful tool for project scheduling in many fields [7, 8], has been increasingly used in project scheduling and subject to a number of refinements and extensions [9]. There are two main components in the implementation of the CCPM, including buffer management and critical chain scheduling.
There are three major types of buffers, including resource buffer (RB), feeding buffer (FB), and project buffer (PB) in CCPM, which are distinguished by their positions and functions in the schedules [10]. RBs are set to protect the critical chain from the tightness of critical resources. They are used as warnings, consuming no time. FBs are set to protect the critical chain from the variation of tasks not on the critical chain and are placed where noncritical and critical activities converge. PBs are placed at the end of the critical chain to protect against exceeding project delivery dates. The concentration principle of buffer determination in CCPM is shown as Figure 1, where represents the estimated duration of activity i and denotes the safety time of . In fact, another two major time buffering methods existed for developing reliable schedules and protect them from effects of uncertainties. One is optimizationbased time buffering method such as starting time criticality and Tabu search, which insert buffers in front of activities dispersedly to ensure that each activity can be conducted on time in accordance with the project schedules. However, it is unknown whether such scattered buffers can ensure ontime completion of the entire project [11]. The other is simulationbased time buffering method. However, a major drawback of simulationbased methods exists that simulation is highly computationally demanding, which could become a significant problem in practice when large projects are involved or when project managers want to perform whatif analysis for a wide variety of scenarios [12]. The CCPMbased time buffering used in this study has several advantages. It calculates time buffers efficiently, provides explicit protection against stochastic variation and presents a set of reasonable control guidelines to ensure ontime completion of the entire project [13].
Determining the size of buffers depends on various factors, such as managerial experience and preferences of the project team, project circumstances, personnel and equipment capabilities, and so on [14]. A number of buffer sizing approaches have been developed in CCPM, among which the cut and paste method (C&PM) and the root square error method (RSEM) are widely used. The C&PM reduces the safety time of each activity by 50% and estimates half of the sum of safety time to be a buffer. However, the C&PM adopts a linear procedure, and the size of the calculated buffer increases linearly with the length of the critical chain. The RSEM calculates the square root of the total squares of the difference between two estimates, a safe estimate, and an average estimate, for each activity in the critical chain as the buffer size [15]. Tukel et al. [16] argued that the RSEM could be less affected by the length of the critical chain than the C&PM. In addition, Icmeli and Erenguc [17] proposed a resource utilization factor (RF) to calculate buffers. Roghanian et al. [18] took into account variation of fuzzy numbers for required resources when buffer sizing. Ghaari and Emsley [9] employed Monte Carlo simulation method to validate the efficiency of C&PM in buffer sizing and put forward buffer sizing schemes in multiproject environments considering the level of resource capacity. A major limitation of these existing buffer sizing approaches, however, is that they operate based on an assumption that project activity durations are independent of each other (Li et al., 2012). These approaches, therefore, are not appropriate to cope with rework, which involves extensively the mutual relationships between activities. In a recent study, Zhang et al. [19] proposed a CCPMbased scheduling method that, by integrating the DSM, took into account rework risks in the calculation of resource tightness and information constraints. However, this method failed to recognize that the rework risks could also impact the critical chain, and it identified the critical chain in a traditional manner.
With respect to critical chain scheduling, Francisco and Alfonso [20] extended traditional critical chain scheduling methods by using critical set and critical cloud to avoid ambiguity in the identification of critical tasks. Ma et al. [21] proposed a critical chain scheduling method with a scenariobased proactive robustness optimization. The method was able to yield higher probabilities of ontime project completion than traditional CCPMbased scheduling methods. Goto [22] developed a scheduling framework for projects under limited resources by representing the CCPM method in maxplus algebra to achieve a short lead time. Salama et al. [23] took into account the impacts of resource continuity and variability on activity durations and integrated linear scheduling and CCPM methods to identify multiple critical sequences. Despite active research in CCPMbased linear process scheduling in recent years, the challenge of rework in CCPMbased schedules has barely been addressed [24], where further research is needed.
2.2. DSM as a Process Management Technique
Since Steward [25] first introduced the DSM for the squarematrixbased models of processes, process architecture DSM models (process DSMs) have received the most attention among different application areas of DSM [6]. The DSM uses matrix format to uniformly represent three different logical relationships and dependencies, including sequential, parallel, and coupled, between activities, as shown in Figure 2. An offdiagonal element of the DSM matrix represents with a mark or a number the information flow from activity to activity or the dependency relationship between these two activities. If no dependency exists between activities and , it will be null or zero. Diagonal element is occasionally used to represents certain characteristics, such as duration, of the activity , or simply nothing in most cases.
(a) Sequential
(b) Parallel
(c) Coupled
In this matrixbased display, elements above the diagonal are usually regarded as information feedback or iterative activities relationship [5]. Such display provides a useful tool for analyzing rework and has given rise to a number of relevant studies. The first DSMbased, discrete event Monte Carlo simulation model was proposed by Browning [26], and later extended by Browning and Eppinger [27], to estimate project duration and cost and their variances. They argued that processes with the fewest feedback marks in the DSM may not be necessarily optimal and could be sped up with appropriate increase in overlapping and iterations [6]. They introduced rework probability (RP) matrix, rework impact (RI) matrix, and learning curve (LC) to represent and calculate rework duration [27]. For the element in row and column in , if , which means the element belongs to upper triangular, it describes the probability that the completion of activity causing rework of activity j; if , which means the element belongs to lower triangular, it describes the probability that after activity completes rework, the activity will be influenced. represents the possibility of activity to be reworked when rework is caused by activity for . represents the ratio of activity to be reworked because the participants may benefit from learning and adaption when rerunning the activity. Based on their definition, the values of the elements in these two matrices are between 0 and 1. The values of elements in the two matrices can be collected from historical data or estimated based on risk preferences of the project team. Then the expected rework time of activity caused by activity can be calculated by , where represents the duration of activity i.
A number of other DSM simulations have adopted the above framework, with certain extensions to account for additional constraints in process scheduling. For instance, Cho and Eppinger [28] proposed a heuristic for solving stochastic, resourceconstrained project scheduling problems in an iterative project network. LĂ©vĂˇrdy and Browning [29] accounted for technical performance characteristics besides duration and cost, by setting up a superset of general classes of activities, each with modes that vary in terms of inputs, duration, cost, and expected benefits. Meanwhile, a number of studies have looked into transferring process DSM to a more effective process management tool that is closer to traditional process management tools and expressions. A time factor was therefore introduced [30]. The time factor is defined based on the overlapping relationship between the activities. As shown in Figure 3, and are defined as the duration of activities and j, respectively, as the duration between the start of activity and the time when information transition from activity to activity is completed, and as the duration between the start of activity j and the time the information transition is completed. A predecessor time factor is defined as and a successor time factor as , where . Two matrices, P and S, are defined by these two factors. The project duration is then calculated according to (1), and a nomenclature is provided in Table 1.

In short, prior research has explored the possibility of using the DSM to predict and control rework in project scheduling and planning and has achieved noticeable progress. However, DSM lacks the capability of functioning as a standalone project management technique, for that DSMâ€™s matrixbased expression is difficult to be effectively converted to network diagrambased expression, sometimes leading to ambiguity and confusion [5]. Therefore, it is challenging to apply DSM to the process management of construction projects, and DSM needs to be extended and improved with further research efforts.
3. The Critical Chain Design Structure Matrix Method
In this section, we are going to demonstrate how to use the CCDSM method to generate the project schedule. We improve the maxplus algorithm to reduce computational load in generating project schedules and introduce rework buffers to evaluate and address rework risks.
3.1. Construction Project Scheduling with the MaxPlus Algorithm
A discrete event system is statediscrete and eventdriven. Its state evolution depends on asynchronous discrete events occurring at discrete points over time [31]. Simple discrete event systems are usually linear systems in which state and output variables for all possible input variables and the initial state satisfy superposition principle of the system. A construction project can generally be seen as a linear system that consists of a number of discrete events.
The maxplus algorithm provides an alternative way to describe the discrete event system. Based on the maxplus algorithm, the status of all discrete events, denoted as vector , in a system can be used to describe the status of the system, and adjacent statuses are interchangeable with a limited number of linear changes, as follows:Based on the above equation, if the linear changes at every stage, , are known, then the final status of the system, , can be derived from the initial status of the system . Using this as a starting point, the maxplus algorithm introduces four types of operations to describe the discrete event systems. For two matrices and denoted as , whose elements are nonnegative numbers or negative infinity, these operations are defined as follows:where could be further simplified as , and .
To apply the maxplus algorithm to schedule management, three new matrices need to be developed to describe relevant project information. These matrices include ) activity duration matrix, denoted as D, which is developed by filling the diagonal elements with each activityâ€™s duration and the nondiagonal elements with negative infinity; activity deviation matrix, denoted as F, which is developed based on the rule that represents the deviation of the start time of activity and activity i; and activity relationship matrix, denoted as R, which is developed based on the rule that if there is information delivered from activity to activity i, then the value of is 0; otherwise the value is negative infinity.
According to the maxplus algorithm [32], we have the following.
The early finish time for each activity is calculated as where . The th element of matrix means the largest deviation of the start time of activity and activity i, if activities and are on one or more paths at the same time; otherwise, it is . For simplicity, the elements of vector are set as e, which means the beginning of activities in the project will not be affected by other construction projects.
The length of critical path, denoted as the maximum of early finish time for all activities, is calculated as
The early start time is calculated as the difference between early finish time for each activity and its duration according to
The late start time of activity is calculated as the difference between length of critical path and the sum of activity â€™s duration as well as subsequent critical activitiesâ€™ durations according to
The total float is calculated as
The critical chain is then determined by the set of activities that satisfy .
It needs to be noted that when determining floats, the maxplus algorithm is run in simple matrix form which significantly improves the computational efficiency and is more applicable to largescale projects [33]. The definitions of parameters in the above equations are summarized in the nomenclature in Table 1.
3.2. Determining Project Rework Buffer
It is widely believed that information uncertainty is the main cause of rework [34]. Such uncertain information may transmit from upstream activities to downstream activities or feedback in the opposite direction, creating information flow. Information flow interactions lead to rework risks in construction projects and consequently bring about more rework time. There are two main situations that rework exists.
First rework: as depicted in Figure 4(a), i is a precedent activity of j, and they perform sequentially according to relationships. After it is completed, the performance information of will be generated and transmitted to i, which may result in rework of . The first rework time, shown as the shadow area in Figure 4(a), can be calculated by [27]Then the total first rework time of caused by all its downstream activities can be obtained by where represents the number of activities in the schedule.
(a) First rework
(b) Second rework caused by information transition from its precedent activity
(c) Second rework caused by information feedback from its successive activity
Second rework: as depicted in Figures 4(b) and 4(c), after completing the rework of caused by , transmits or feeds back some revised information to , which may cause rework of . Such two forms are defined as the second rework in the paper. The second rework times, shown as the shadow areas in Figures 4(b) and 4(c), can be calculated byThen the total second rework time of caused by all its downstream activities can be obtained by
To proactively account for the uncertainty caused by interactivity relationships and the associated rework risks, a rework buffer () is introduced in the proposed CCDSM. Conceptually, assesses and compresses the duration uncertainty caused by rework and is placed after the completion of rework activities in project schedules. is designed to warn project participants from how much workload may increase for each rework activity, so as to prepare the project participants for sufficient time and resources.
The rework buffer of can be obtained by
For the last critical activity in the schedule, considering all reworks it causes other activities to generate will take place after its completion and have an impact on project completion time; therefore, its rework buffer is defined as the sum of the total rework time it causes other activities to generate and its total rework time caused by other activities.
The total rework time it causes all other activities to generate can be calculated by The rework buffer of then can be obtained by
3.3. Determining Project Buffers and Feeding Buffers in MaxPlus Representation
3.3.1. Determining Project Buffers
Previous studies on the CCPMmaxplus representation adopted the C&PM method for calculating project buffers, shown as (16), and feeding buffers [33]. where th element of matrix means the cumulative time of duration of activity , which is on the critical chain and range from to j; th element of matrix satisfies that if activity is on the critical chain, then th element is duration of activity I; otherwise it is .
However, the C&PM method lacks sound mathematical foundation and overestimates project durations, resulting in a waste of time and resources. Alternatively, the C&PM is replaced by the RSEM, which is based on the large number law and central limit theorem, in the CCDSM to calculate the PB: where is the variance of duration of activity on the critical path.
To calculate the PB in maxplus representation, a matrix denoted as is introduced. For the element in row i and column in , if activity is critical activity, it will be filled with the variance of duration of activity i; else, it will be zero. For other elements in , it is filled with zero. The PB then can be represented in maxplus algebra aswhere .
3.3.2. Determining Feeding Buffers
In order to insert feeding buffers in place to protect critical path from delays in noncritical paths, two operators, denoted as and #, are first introduced. For a vector , a new matrix can be obtained by the operator as , and a new vector can be obtained by the operator # as .
In preparation of the location and size of the feeding buffer, another two vectors are introduced, denoted as and . Matrix is developed based on the rule that for a diagonal element; if the corresponding activity is on the critical path, then the element value is 0; otherwise, the element value is . Matrix is developed based on the rule that, for a diagonal element, if the corresponding activity is on the noncritical path, then the value is 0; otherwise, the value is . Moreover, an adjacency matrix is transformed from matrix by . Vector is then introduced to locate the feeding buffers as follows:where is the adjacency matrix, which represents transitions from noncritical activities to critical ones and can be obtained by . represents the noncritical activity, one of whose successors is critical activity. Therefore, a feeding buffer should be inserted behind the activity .
The feeding buffer can then be determined by square rooting of the sum of , where is the set of activities on a certain noncritical chain. The formula is expressed in maxplus algebra as
Additionally, the calculation of feeding buffer should be adjusted using the following equation to meet the constraint that the size of feeding buffer cannot exceed total float:
3.4. Implementation of the CCDSM with MaxPlus Linear Expression
The maxplusbased implementation of the CCDSM is summarized in the following seven steps. It needs to be noted that, as the maxplus algorithm cannot solve nonlinear conversions such as multiplication, the input data for the project rework problem including matrices A, P, S, RP, and RI needs to be preprocessed:
Step 1 (preprocess ). Extract all diagonal elements of matrix and form a diagonal matrix D.
Step 2 (preprocess and ). Combine matrices and to generate matrix according to the following equation:Extract all nonnull elements of matrix and use them to form matrix . Then update R, by replacing all nonnull elements with and all null elements with .
Step 3 (preprocess RP and RI). Based on (9)(15), merge matrices RP and RI to generate the vector , whose th element is the corresponding rework buffer of activity i.
Step 4 (calculate the length of critical path). Based on (4)(8), the length of critical path TT and other parameters, including the total float time TF and critical path are determined.
Step 5 (calculate matrices and ). According to the critical path determined in Step 4, matrices and can be generated. Meanwhile, the activity duration variance matrix SD can be generated based on data collected from the project to calculate the buffers size;
Step 6 (calculate buffer sizes). Based on (17)(21), project buffer PB and feeding buffers FBs are calculated with RSEM and represented in maxplus algebra.
Step 7 (generate the schedule). As the principle of CCPM, the feeding buffers should be placed on the noncritical chains prior to the joints of the critical chain and noncritical chains, and the project buffer should be placed at the end of the schedule to protect the whole project process. Then, the estimated project time PT can be calculated based on
4. Case Study
In this section, the proposed CCDSM method is implemented in a case project to demonstrate its feasibility and effectiveness in addressing rework risks in project schedule management. The settings of the case project are first presented, followed by descriptions of implementation of the CCDSM in the case project. The schedule generated with the CCDSM method is assessed in detail and compared with the schedules generated with traditional CCPM and DESbased methods.
4.1. Case Project Settings
The case used in this study was derived from a modular real estate development project firstly introduced in [35] and further described in [36]. The matrix of the case is shown in Figure 5(a). This case consisted of 19 major activities and 183 interactivity relationships, including 65 rework relationships represented by elements above the diagonal of matrix . It is assumed that rework can propagate up to twice to avoid infinite loop in the computation. The parameters for each activity and the interactivity relationships were derived based on literature and empirical evidence, and are shown in Figures 5(b), 5(c), 5(d), and 5(e) and summarized in Table 2. Largescale rework relationships and complex connections between activities were observed in the case project, which was representative of typical construction projects in reality.

(a) Matrix A
(b) Matrix P
(c) Matrix S
(d) Matrix RP
(e) Matrix RI
4.2. Implementation of the CCDSM Method
The proposed CCDSM method was implemented in the case project following the steps explained in the last section and the results were shown below. Noticeably, for better understanding of scheduling process, the implementation process of the CCDSM method is reorganized as below.
Step 1. Matrices and were derived from matrix A, as shown in Figure 6. The duration of each activity was presented in the diagonal of the matrix D.
(a) Matrix D
(b) Matrix R (element is omitted)
Step 2. Matrix F (see Figure 7) was derived from matrices and according to (22). This step calculated the activity deviation matrix based on the predecessor time matrix and the successor time matrix.
Step 3. Optimize activity sequence and update matrices D, R, RP, RI, and according to new activity sequence. To reduce project rework and obtain nearoptimal project completion time, the genetic algorithm (GA) was applied to optimize activity sequence. The GA is a metaheuristic method that searches for optimal solutions using processes similar to those in natural selections and genetics [37]. In the paper, minimization of total length of rework path first proposed by Gebala and Eppinger [38] was used as the objective function to calculate the optimal sequence of activities:where represents rework probability of activity caused by activity .
The parameter settings were selected as follows: population size set as 50, number of generations set as 150, crossover probability set as 0.95, and mutation probability set as 0.08. The GA process reported an optimal schedule as . The GA convergence process is shown in Figure 8.
Step 4. was generated from matrices RP and RI based on (9)(15). The goal of this step was to factor in and calculate the rework time of the case project. was .
Step 5. To generate reliable project schedules, the most likely duration of each activity, which has been widely used and accepted by project teams in prior research [39], was selected to calculate the estimated project duration. The critical chain was reported as .
Step 6. The activity duration variance matrix SD was generated with 1,000 Monte Carlo simulations of the activitiesâ€™ duration that follows a beta distribution.
Step 7. Buffers were calculated based on (17)(21). The project buffer was 36 days, and the feeding buffer was .
Step 8. The estimated project duration was calculated as 102 days, based on (23). The project buffer was placed at the end of the schedule, and the feeding buffer was placed on the noncritical chains prior to the joints of the critical chain and noncritical chains, following the principle of CCPM.
It was worth mentioning that overlapping was observed among various activities, including critical activities, in the above schedule, which illustrated realworld highfrequency interactions among activities in the project.
To assess the performance of the CCDSM method and the quality of the project schedule it generated, the execution of case project was simulated 1000 times using Monte Carlo approach, where the sampling of actual activity durations followed beta distribution. The simulated actual project durations were depicted in Figure 9. It can be seen in the figure that under rework assumption, the process schedule of 102 days had a completion probability of 98.1%, indicating that the proposed CCDSM method was able to provide satisfactory assurance of ontime project completion under varying scheduling scenarios.
4.3. Comparison of the CCDSM Method with Traditional CCPM and DES Methods
This subsection further compares the performance of the CCDSM method with traditional CCPM and DES methods to assess its effectiveness in addressing rework risks.
A CCPMbased schedule was generated with traditional buffers including project buffer and feeding buffer and without rework buffer. The results indicated that the estimated project durations of the CCPMbased schedule were 93 days, including project buffers of 36 days. When the estimated project duration of CCPMbased schedule was superimposed in the Monte Carlosimulated project durations, as depicted in Figure 10, it can be seen that the ontime completion probability was 49.3%. The results suggested that the CCDSMbased schedule, which took into consideration rework relationships and had a duration of 102 days, significantly outperformed the CCPMbased schedule in ensuring ontime completion of the case project that was faced with typical rework risks.
Furthermore, a DESbased schedule was generated by implementing DES algorithm, in which activity durations following a beta distribution were generated at random to compute the project duration. When the DES algorithm was used to generate project schedules under rework scenarios, a probabilistic judgment is performed firstly on whether or not each activity causes rework, and then an iterative simulation process began until all activities were finished. Specifically, the DES algorithm proposed by Browning and Eppinger [27] was adopted in this study, with the following minor modifications to make it more applicable to the case project and comparable to the CCDSM method: the optimized activity sequence in Section 4.2 was employed in the DES process; activity did not started until STS logic tie with its each predecessor activity ; that is , was satisfied; rework probability of activity would decrease by 50% each time it reworked [28].
To compare the CCDSMbased schedule with the DESbased schedule, the execution of the case project using the DESbased schedule was simulated 1,000 times using Monte Carlo approach, and the simulated actual project durations are depicted in Figure 11. The results showed that the estimated project duration of the DESbased schedule was 115 days, which was 13 days longer than that of the CCDSMbased schedule, with the same ontime completion probability of 98.1%. In other words, the CCDSM method outperformed the DES method by ensuring the same ontime completion probability with a shorter project schedule. Moreover, the CCDSM method was also advantageous over the DESbased method in that it enabled effective control of project schedules during project implementation, which could be done by monitoring consumption of time buffers, leading to reduction in project overrun risks.
5. Conclusions
Construction projects are constantly challenged by rework risks, which have largely remained unaddressed by existing project scheduling. This paper proposed the CCDSM method for developing project schedules that are adaptive to rework instances during project execution. This method models largescale rework relationships among activities with the introduction of a new rework buffer in traditional CCPM method. The method allows analysis and visualization of the schedules and utilizes maxplus algorithm to transform complex logic relationships into simple matrix operations, reducing the computational load of generating process schedules. A case study was conducted to demonstrate the implementation of the CCDSM method and assess its effectiveness in addressing rework risks. The results showed that CCDSMbased schedule outperformed the CCPMbased schedule in ensuring ontime completion of the case project that was faced with typical rework risks. In addition, in comparison to the DES method, the CCDSM performed better by providing a smaller project duration with the same probability of ontime project completion.
The CCDSM method contributes to construction management in the following ways. First, the CCDSM method absorbs risks of rework and various other uncertainties to generate the reliable project schedule, enabling managers to predict project durations more accurately before projects start. Based on accurate prediction of project durations, the negotiation of contracts and supply of resources can be conducted in a more reasonable manner, avoiding costs overruns, schedule delays, and even project failures. Second, based on generated project schedule with time buffers, effective monitoring and control of construction progress can be realized. Overall speaking, the proposed CCDSM method provides a promising solution to mitigate rework risks in construction projects and protect projects from undesirable delays caused by possible rework instances.
Meanwhile, the CCDSM bears several limitations that would benefit from future improvement. First, only the STS logic tie between activities is tested in this study. Techniques for integrating multiple logic ties including STS, STF, FTS, and FTF deserve further attention. Second, resource conflicts and levering are considered beyond the scope of this study, but it is a critical challenge that should be addressed in the context of CCDSMbased scheduling in future research. Third, the parameters required to implement the CCDSM method are hard to obtain for a specific project. To employ CCPM or maxplus algorithm for project scheduling, it is necessary to obtain the estimated activity durations under various completion rates and the logic tie between activities. Similarly, the rework relative parameters including RP, RI, and LC are required for implementation of DSMbased method. As an integration method of CCPM, maxplus algorithm, and DSM, the implementation of CCDSM requires all the above parameters to be obtained. Lastly, methods for dynamic updating of buffer sizes during project execution are needed in order to better prepare project teams for changing project environment and dynamic project progress.
Data Availability
All data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (NSFC) [Grant no. #71671128].
References
 Y. Li and T. R. Taylor, â€śModeling the impact of design rework on transportation infrastructure construction project performance,â€ť Journal of Construction Engineering and Management, vol. 140, no. 9, p. 04014044, 2014. View at: Publisher Site  Google Scholar
 B. Hwang, X. Zhao, and K. J. Goh, â€śInvestigating the clientrelated rework in building projects: The case of Singapore,â€ť International Journal of Project Management, vol. 32, no. 4, pp. 698â€“708, 2014. View at: Publisher Site  Google Scholar
 E. K. Simpeh, R. Ndihokubwayo, P. E. Love, and W. D. Thwala, â€śA rework probability model: a quantitative assessment of rework occurrence in construction projects,â€ť International Journal of Construction Management, vol. 15, no. 2, pp. 109â€“116, 2015. View at: Publisher Site  Google Scholar
 P. W. Morris, J. K. Pinto, and S. Jonas, The Oxford Handbook of Project Management, Oxford University Press, Oxford, UK, 2012.
 S. D. Eppinger and T. R. Browning, Design Structure Matrix Methods and Applications, MIT Press Books, 2012. View at: Publisher Site
 T. R. Browning, â€śDesign structure matrix extensions and innovations: a survey and new opportunities,â€ť IEEE Transactions on Engineering Management, vol. 63, no. 1, pp. 27â€“52, 2016. View at: Publisher Site  Google Scholar
 A. Kulkarni, D. K. Yadav, and H. Nikraz, â€śAircraft maintenance checks using critical chain project path,â€ť Aircraft Engineering Aerospace Technology, vol. 89, no. 6, pp. 879â€“892, 2017. View at: Publisher Site  Google Scholar
 J. Trojanowska and E. Dostatni, â€śApplication of the theory of constraints for project management,â€ť Management and Production Engineering Review, vol. 8, no. 3, pp. 87â€“95, 2017. View at: Publisher Site  Google Scholar
 M. Ghaari and M. Emsley, â€śBufer sizing in CCPM portfolios with dierent resource capacities,â€ť International Journal of Information Technology Project Management, vol. 8, no. 3, pp. 40â€“51, 2017. View at: Google Scholar
 E. M. Goldratt, Critical Chain: A Business Novel, MA: North River Press, Great Barrington, United States, 1997.
 S. Van de Vonder, E. Demeulemeester, and W. Herroelen, â€śProactive heuristic procedures for robust project scheduling: An experimental analysis,â€ť European Journal of Operational Research, vol. 189, no. 3, pp. 723â€“733, 2008. View at: Publisher Site  Google Scholar
 O. Lambrechts, E. Demeulemeester, and W. Herroelen, â€śTime slackbased techniques for robust project scheduling subject to resource uncertainty,â€ť Annals of Operations Research, vol. 186, no. 1, pp. 443â€“464, 2011. View at: Publisher Site  Google Scholar
 D. Trietsch and K. R. Baker, â€śPERT 21: Fitting PERT/CPM for use in the 21st century,â€ť International Journal of Project Management, vol. 30, no. 4, pp. 490â€“502, 2012. View at: Publisher Site  Google Scholar
 G. Ma, A. Wang, N. Li, L. Gu, and Q. Ai, â€śImproved critical chain project management framework for scheduling construction projects,â€ť Journal of Construction Engineering and Management, vol. 140, no. 12, p. 04014055, 2014. View at: Publisher Site  Google Scholar
 R. C. Newbold, Project Management in The Fast Lane: Applying The Theory of Constraints, CRC Press, Boca Raton, Florida, London, 1998. View at: Publisher Site
 O. I. Tukel, W. O. Rom, and S. D. Eksioglu, â€śAn investigation of buffer sizing techniques in critical chain scheduling,â€ť European Journal of Operational Research, vol. 172, no. 2, pp. 401â€“416, 2006. View at: Publisher Site  Google Scholar
 O. Icmeli and S. S. Erenguc, â€śA branch and bound procedure for the resource constrained project scheduling problem with discounted cash flows,â€ť Management Science, vol. 42, no. 10, pp. 1395â€“1408, 1996. View at: Publisher Site  Google Scholar
 E. Roghanian, M. Alipour, and M. Rezaei, â€śAn improved fuzzy critical chain approach in order to face uncertainty in project scheduling,â€ť International Journal of Construction Management, vol. 18, no. 1, pp. 1â€“13, 2017. View at: Publisher Site  Google Scholar
 J. Zhang, X. Song, and E. DĂaz, â€śProject buffer sizing of a critical chain based on comprehensive resource tightness,â€ť European Journal of Operational Research, vol. 248, no. 1, pp. 174â€“182, 2016. View at: Publisher Site  Google Scholar
 F. A. Rivera and A. Duran, â€śCritical clouds and critical sets in resourceconstrained projects,â€ť International Journal of Project Management, vol. 22, no. 6, pp. 489â€“497, 2004. View at: Publisher Site  Google Scholar
 G. Ma, L. Gu, and N. Li, â€śScenariobased proactive robust optimization for criticalchain project scheduling,â€ť Journal of Construction Engineering and Management, vol. 141, no. 10, p. 04015030, 2015. View at: Publisher Site  Google Scholar
 H. Goto, â€śForwardcompatible framework with criticalchain project management using a maxplus linear representation,â€ť OPSEARCH, vol. 54, no. 1, pp. 201â€“216, 2017. View at: Publisher Site  Google Scholar
 T. Salama, A. Salah, and O. Moselhi, â€śIntegration of linear scheduling method and the critical chain project management,â€ť Canadian Journal of Civil Engineering, vol. 45, no. 1, pp. 30â€“40, 2018. View at: Publisher Site  Google Scholar
 C. W. Hu, X. D. Chen, and L. H. Wu, â€śThe setting method of project buffer in critical chain management of mould manufacturing project based reworking,â€ť Advanced Materials Research, vol. 317319, pp. 418â€“422, 2011. View at: Publisher Site  Google Scholar
 D. V. Steward, Systems Analysis and Management: Structure, Strategy and Design, Irwin/McGrawHill, Boston, 1981.
 T. R. Browning, Modeling and analyzing cost, schedule, and performance in complex system product development [Ph.D. thesis], Massachusetts Institute of Technology, Cambridge, Massachusetts, USA, 1998.
 T. R. Browning and S. D. Eppinger, â€śModeling impacts of process architecture on cost and schedule risk in product development,â€ť IEEE Transactions on Engineering Management, vol. 49, no. 4, pp. 428â€“442, 2002. View at: Publisher Site  Google Scholar
 S. H. Cho and S. D. Eppinger, â€śA simulationbased process model for managing complex design projects,â€ť IEEE Transactions on Engineering Management, vol. 52, no. 3, pp. 316â€“328, 2005. View at: Publisher Site  Google Scholar
 V. LĂ©vĂˇrdy and T. R. Browning, â€śAn adaptive process model to support product development project management,â€ť IEEE Transactions on Engineering Management, vol. 56, no. 4, pp. 600â€“620, 2009. View at: Publisher Site  Google Scholar
 J. U. Maheswari and K. Varghese, â€śProject scheduling using dependency structure matrix,â€ť International Journal of Project Management, vol. 23, no. 3, pp. 223â€“230, 2005. View at: Publisher Site  Google Scholar
 C. G. Cassandras and S. Lafortune, Introduction to Discrete Event Systems, Springer, New York, NY, USA, 2nd edition, 2008. View at: Publisher Site  MathSciNet
 S. Yoshida, H. Takahashi, and H. Goto, â€śModified maxplus linear representation for inserting time buffers,â€ť in Proceedings of the IEEE International Conference on Industrial Engineering & Engineering Management, IEEE, 2010. View at: Google Scholar
 S. Yoshida, H. Takahashi, and H. Goto, â€śResolution of time and worker conflicts for a single project in a maxplus linear representation,â€ť Industrial Engineering & Management Systems, vol. 10, no. 4, pp. 279â€“287, 2011. View at: Google Scholar
 S. Muralidharan, â€śFactors affecting rework in construction project,â€ť International Journal of Engineering Sciences & Research Technology, vol. 1, no. 5, pp. 578â€“584, 2016. View at: Publisher Site  Google Scholar
 B. B. E. Bulloch and J. Sullivan, Application of the Design Structure Matrix (DSM) to the real estate development process. [Master, thesis], Massachusetts Institute of Technology, Cambridge, Massachusetts, USA, 2009.
 S. D. Eppinger, S. Bonelli, and A. M. Gonzalez, â€śManaging iterations in the modular real estate development process,â€ť Reducing Risk in Innovation: Proceedings of the 15th International DSM Conference Melbourne, Australia, 2930 August 2013, pp. 37â€“44, 2013. View at: Google Scholar
 D. T. Pham and D. Karaboga, Intelligent Optimisation Techniques: Genetic Algorithms, Tabu Search, Simulated Annealing and Neural Networks, Springer Science & Business Media, 2000. View at: MathSciNet
 D. A. Gebala and S. D. Eppinger, â€śMethods for analyzing design procedures,â€ť in Proceedings of the AMSE Third International Conference On Design Theory and Methodology, vol. 31, pp. 227â€“233, Miami, Florida, 1991. View at: Google Scholar
 Ă–. HazÄ±r, M. Haouari, and E. Erel, â€śRobust scheduling and robustness measures for the discrete time/cost tradeoff problem,â€ť European Journal of Operational Research, vol. 207, no. 2, pp. 633â€“643, 2010. View at: Publisher Site  Google Scholar  MathSciNet
Copyright
Copyright © 2019 Guofeng Ma et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.