Abstract

The vast majority of the research efforts in project risk management tend to assess cost risk and schedule risk independently. However, project cost and time are related in reality and the relationship between them should be analyzed directly. We propose an integrated cost and schedule risk assessment model for complex product systems R&D projects. Graphical evaluation review technique (GERT), Monte Carlo simulation, and probability distribution theory are utilized to establish the model. In addition, statistical analysis and regression analysis techniques are employed to analyze simulation outputs. Finally, a complex product systems R&D project as an example is modeled by the proposed approach and the simulation outputs are analyzed to illustrate the effectiveness of the risk assessment model. It seems that integrating cost and schedule risk assessment can provide more reliable risk estimation results.

1. Introduction

Complex product systems can be defined as high cost, engineering-intensive products, systems, networks, and constructs [1]. Over the past two decades, the importance of complex product systems has been recognized and investigated by researchers (e.g., [2, 3]). Risk assessment of cost and schedule, which plays an important role in complex product systems R&D project management, is concerned and studied closely [4].

The traditional approach to estimating risk of cost and schedule has been to estimate them independently (e.g., [5, 6]). However, a number of limitations of independent estimation are pointed out by researchers. For example, Xu et al. [7] set up the estimation models of cost risk, schedule risk, and integrated risk, which show that the integrated risk probability estimation is different from independent estimation. Additionally, Isidore and Back [8] verified that choosing a schedule value having a high confidence level does not guarantee that the associated cost estimate, corresponding to the selected schedule value, will also have a comparably high level of confidence, and vice versa. Hence, it is valuable to study the integrated risk.

In recent years, many approaches have been proposed to carry out relative research into integrated analysis problems of stochastic cost and schedule. Carr [9] presented equations which described parallel hierarchical cost and schedule control systems that calculate cost, schedule, and time variances at different levels of detail. Lee and Yi [10] utilized matrix analysis method which integrated the data of cost and duration in the form of matrix to study project integrated plan. Sen et al. [11] proposed an integrated probabilistic risk analysis for systems with multiple state variables. Robert and Randal [12] discussed a technique for integrating the uncertainties associated with cost and schedule when using the Monte Carlo simulation tool. A method to study integrated optimization plan of cost and schedule with a higher confidence level by using activity based costing simulation and multiple-simulation analysis technique was presented by Isidore and Back [13, 14]. Xu et al. [15] proposed an integrated confidence regression model of cost and schedule with Monte Carlo multiple simulation and regression analysis method. Poh and Tah [16] modeled risk impacts with integrated cost and schedule influence network for the construction project. Seyedhoseini et al. [17] introduced an integrated methodology with regard to time, quality, and cost which performed the project risk analysis in an effective manner. Furthermore, many other researchers discussed the techniques for analysis of risk to schedule and cost (e.g., [1821]).

Xu et al. [7] first introduced the concept of integrated cost and schedule risk as the risk of failing to achieve the projects under a certain cost objective () and schedule plan (). More exactly, it is expressed as an integrated bivariate probability distribution function of cost and schedule: They also take advantage of frequency statistical analysis method to model integrated risk probability estimation with Monte Carlo simulation.

In complex product systems R&D project, rework of an activity is common, which is usually caused by probabilistic failure to meet the planned design objective. GERT considered such failure probability of activities and could be used to model and analyze the process of complex product systems with Monte Carlo simulation [22, 23]. Taylor and Moore [24] studied R&D project planning with Q-GERT network modeling and simulation. In their research, a failure probability is given to each activity, and this activity can be fed back to its previous stage to be reworked.

In the paper, motivated by the successful use of the GERT, we employ this technique with multifeedback branches to describe the process of complex systems R&D project. On the basis of the theory of probability and the relationship of marginal probability distribution function, conditional probability distribution function, and integrated probability distribution function, we construct an integrated cost and schedule risk estimation model. Finally, an example of a real-life project is analyzed with the integrated risk estimation model in detail. Compared with the prediction estimation using approximate curve surface, our method can estimate risk probability more rapidly and accurately by processing Monte Carlo simulation results.

The remainder of this paper is organized as follows. Section 2 describes the independent cost and schedule risk estimation model. Section 3 presents the integrated cost and schedule risk estimation model. An application of the model to a real-life project example is given in Section 4. Finally, Section 5 concludes the paper.

2. Independent Cost and Schedule Risk Estimation

Suppose that the simulation runs times, , and the simulation output results for total cost and total duration are and , respectively. Utilizing frequency statistical analysis method, the output results are divided into the following intervals as and . Then, we can count the frequency for each interval.

Let be the frequency of the cost interval and let be that of the duration interval . Thus, the marginal probability density functions of cost and schedule are defined as As a result, the marginal cumulative probability distribution functions are Accordingly, the independent risk probability distribution functions are

3. Integrated Cost and Schedule Risk Estimation

3.1. Integrated Risk Probability Estimation

The integrated cost and schedule risk is defined as the probability of failing to complete the project under a specified objective with respect to cost and schedule. The integrated cost and schedule frequency statistics is shown in Figure 1.

Let be the frequency of the shadow region in Figure 1. Then the integrated cost and schedule probability density function is

The integrated probability distribution function in region is where is the cumulative frequency of simulation outputs which fall on the shadow region .

Accordingly, the integrated cost and schedule risk probability distribution function is

3.2. Conditional Risk Probability Estimation

and are components of two-dimensional random vector of cost and schedule, and on the condition that one ( or ) is assigned to a certain value, the probability density function of the other ( or ) is defined as conditional probability density function.

The conditional probability density functions of cost and schedule are The conditional probability distribution functions of cost and schedule are

4. Application Research

4.1. GERT-Based Model

Our integrated cost and schedule risk estimation model is demonstrated with example data from an obstacle clearance armored vehicle (OCAV) R&D project. We utilize GERT to model the OCAV R&D process, which is usually divided into five phases: project feasibility study, project evaluation, subsystems development, system assembly, and integrated testing. The OCAV is composed of six subsystems. To focus on the simulation-based method presented in this paper, we name the six subsystems as subsystems 1 to 6 instead of their actual names.

Figure 2 shows a stochastic multifeedback GERT model of the OCAV R&D project. In this model, we assume that the parameters of duration and cost are random variables and both of them follow the triangular distribution: (as shown in Table 1), where is the best-case estimate, is the most likely estimate, and is the worst-case estimate. The parameters are estimated by project managers in a subjective manner or assessed through historical data.

4.2. Statistical Analysis of Simulation Output

Simulation runs times. From the simulation output, we can get 1000 pairs of data of total cost and duration , . Using frequency statistical method, the simulation outputs are grouped as and , respectively. According to statistics, we obtain the marginal frequency distribution of total cost , the marginal frequency distribution of total duration , and the integrated bivariate frequency distribution , and the results are shown in Table 2.

We can also obtain the following results. (1)According to the simulation results, the mean and standard deviations of the total cost and duration are calculated as follows:  (2)Frequency histograms of total cost and duration can be drawn and are shown in Figures 3 and 4. (3)We perform the goodness of fit test and deduce that the probability distribution of the total cost and duration is normal distribution. More specifically, the theoretical probability distribution functions of the total cost and duration are subject to and , respectively. Besides, the distribution curves are shown in Figures 3 and 4. (4)From the result in (3), we believe that the total cost and duration are all subject to normal distribution. Hence, interval estimation of parameter can be carried out.The confidence interval of 95% for the total cost is The confidence interval of 95% for the total duration is

4.3. Independent Estimation of Risk

We can get the marginal probability density functions of the total cost and total duration and , the marginal probability distribution functions and , and the marginal risk probability distribution functions and by substituting data into formulas (1)–(3); all the results are shown in Table 3.

The marginal probability distribution and the marginal risk probability distribution of the total cost and the total duration are shown in Figures 5 and 6, respectively.

4.4. Integrated Estimation of Risk

According to formulas (6) and (7), we can get the integrated probability distribution of the total cost and duration in Table 4 and the integrated risk probability distribution of total cost and duration in Table 5.

According to the data in Table 5, we draw the curved surface of integrated risk probability distribution of total cost and duration in Figure 7.

4.5. Conditional Estimation of Risk

According to formula (9), we can get the conditional probability density distributions of the total cost and total duration and in Tables 6 and 7.

According to the above results, the smooth curves of conditional probability distribution of the total cost and the total duration are drawn in Figures 8 and 9.

Using goodness of fit test method, we infer that the conditional probability distributions of the total cost and the total duration follow normal distributions, and the goodness of fit is good. Meanwhile, the mean and the standard deviation of the conditional probability distribution are obtained. The results are shown in Table 8; the computational process is similar to the following example:

4.6. Prediction Estimation of Integrated Risk Probability
4.6.1. Prediction Estimation Based on Curved Surface of Integrated Risk Probability

Through the data in Table 5, we establish the integrated risk space surface of two-dimensional continuous random vector of the cost and duration. Also, we obtain the integrated risk probability approximation under cost and duration constraints by surface interpolation method. For example, is the risk probability value under constraints of the total cost thousand RMB and the total duration months. In Figure 6, when the number of simulations is large enough, the small curved surface which consisted of , , , and can be considered as a plane. In order to avoid underestimating the probability value, we choose , , and which have higher risk probability to establish equation:

After that, the risk probability value is acquired under constraints of total cost and total duration . We can substitute , , into (14) to gain .

4.6.2. Prediction Estimation Based on Theory of Integrated Risk Probability Distribution

Following the definition of the conditional probability distribution function and the marginal probability distribution function, the integrated risk probability of the total cost and duration is calculated as follows:

For example, we also get the risk probability value under the integrated constraints of the total cost ten thousand RMB and the total duration months. The specific calculation steps are as follows.(1)With the data from Table 8, the mean estimation and the standard variance estimation of the conditional distribution functions and are determined by interpolation: (2)Calculate the conditional probability distribution values by normal distribution: (3)According to formula (3), the marginal probability distributions and are obtained as (4)Through formula (15), we calculate the values of the integrated probability distribution and and the integrated risk probability values and ; that is,

4.6.3. Error Analysis of Prediction

By analyzing the output data of 1000 times simulation, we can calculate the frequency of occurrence that the total cost ten thousand RMB and the total duration months are 162. By formula (15), we get

For purposes of comparison, we also report results based on the curved surface method of integrated risk probability. The prediction result is , and the relative error of prediction is . Meanwhile, based on the conditions of risk probability distribution functions proposed in this paper, we have (a) , where the relative error is and (b) , where the relative error is .

The above comparison results of relative error reveal that both of the methods are practical. However, it is too complicated to build the equation when using the approximate curve surface method. And the method proposed in this paper can estimate risk probability more rapidly and accurately by processing one set of Monte Carlo simulation results.

5. Conclusions

Based on GERT multifeedback simulation and the theory of probability distribution, we have presented an integrated risk model of cost and schedule. Using statistical analysis and Monte Carlo techniques, we get the marginal probability distribution functions and the conditional probability distribution functions. Finally, the integrated probability distribution functions of cost and schedule and risk probability distribution function are obtained. This method is proven to be more accurate. However, we know that there are many other factors in complex product systems except cost and schedule, such as environment and project resources. In the future, we will study the risk assessment about more elements.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This research is supported by the Innovation Foundation of BUAA for Ph.D. Graduates and the National Science Foundation of China under Grant 71271019.