CHPTER STOCHSTIC SCHEDULING In some situations, estimating activity duration becomes a difficult task due to ambiguity inherited in and the risks associated with some work. In such cases, the duration of an activity is estimated as a range of time values rather than being a single value. This chapter deals with the scheduling of the project under uncertain activity duration. The program evaluation and review technique is presented as one of the methods that can be used to schedule projects with uncertain activities durations. lso, an introduction to Monte Carlo simulation technique is presented..1 Scheduling with Uncertain Durations Some scheduling procedures explicitly consider the uncertainty in activity duration estimates by using the probabilistic distribution of activity durations. That is, the duration of a particular activity is assumed to be a random variable that is distributed in a particular fashion. For example, an activity duration might be assumed to be distributed as a normal or a beta distributed random variable as illustrated in Figure.1. This figure shows the probability or chance of experiencing a particular activity duration based on a probabilistic distribution. The beta distribution is often used to characterize activity durations, since it can have an absolute minimum and an absolute maximum of possible duration times. The normal distribution is a good approximation to the beta distribution in the center of the distribution and is easy to work with, so it is often used as an approximation. Construction Management 0 Dr. Emad Elbeltagi
Figure.1: Beta and normally distributed activity durations If a standard random variable is used to characterize the distribution of activity durations, then only a few parameters are required to calculate the probability of any particular duration. Still, the estimation problem is increased considerably since more than one parameter is required to characterize most of the probabilistic distribution used to represent activity durations. For the beta distribution, three or four parameters are required depending on its generality, whereas the normal distribution requires two parameters. The most common formal approach to incorporate uncertainty in the scheduling process is to apply the critical path scheduling process and then analyze the results from a probabilistic perspective. This process is usually referred to as the Program Evaluation and Review Technique (PERT) method. s noted earlier, the duration of the critical path represents the minimum time required to complete the project. Using expected activity durations and critical path scheduling, a critical path of activities can be identified. This critical path is then used to analyze the duration of the project incorporating the uncertainty of the activity durations along the critical path. The expected project duration is equal to the sum of the expected durations of the activities along the critical path. Construction Management 1 Dr. Emad Elbeltagi
ssuming that activity durations are independent random variables, the variance or variation in the duration of this critical path is calculated as the sum of the variances along the critical path. With the mean and variance of the identified critical path known, the distribution of activity durations can also be computed..1.1 Program Evaluation and Review Technique Both CPM and PERT were introduced at approximately the same time and, despite their separate origins, they were very similar. The PERT method shares many similarities with CPM. Both require that a project be broken down into activities that could be presented in the form of a network diagram showing their sequential relationships to one another. Both require time estimates for each activity, which are used in routine calculations to determine project duration and scheduling data for each activity. CPM requires a reasonably accurate knowledge of time and cost for each activity. In many situations, however, the duration of an activity can not be accurately forecasted, and a degree of uncertainty exists. Contrary to CPM, PERT introduces uncertainty into the estimates for activity and project durations. It is well suited for those situations where there is either insufficient background information to specify accurately time and cost or where project activities require research and development. In the original development of PERT approach, O notations are used. However, ON diagramming can be easily used alternatively. The method is based on the well-known central limit theorem. The theorem states that: Where a series of sequential independent activities lie on the critical path of a network, the sum of the individual activity durations will be distributed in approximately normal fashion, regardless of the distribution of the individual activities themselves. The mean of the distribution of the sum of the activity durations will be the sum of the means of the individual activities and its variance will be the sum of the activities variances. The primary assumptions of PERT can be summarized as follows: Construction Management Dr. Emad Elbeltagi
1. ny PERT path must have enough activities to make central limit theorem valid.. ny PERT path must have enough activities to make central limit theorem valid. 3. The mean of the distribution of the path with the greatest duration, from the initial node to a given node, is given by the maximum mean of the duration distribution of the paths entering the node.. PERT critical path is longer enough than any other path in the network. PERT, unlike CPM, uses three time estimates for each activity. These estimates of the activity duration enable the expected mean time, as well as the standard deviation and variance, to be derived mathematically. These duration estimates are: - Optimistic duration (a); an estimate of the minimum time required for an activity if exceptionally good luck is experienced. - Most likely or modal time (m); the time required if the activity is repeated a number of times under essentially the same conditions. - Pessimistic duration (b); an estimate of the maximum time required if unusually bad luck is experienced. These three time estimates become the framework on which the probability distribution curve for the activity is erected. Many authors argue that beta distribution is mostly fit construction activities. The use of these optimistic, most likely, and pessimistic estimates stems from the fact that these are thought to be easier for managers to estimate subjectively. The formulas for calculating the mean and variance are derived by assuming that the activity durations follow a probabilistic beta distribution under a restrictive condition. The probability density function of beta distributions for a random variable x is given by: α β f(x) = k(x - a) (x - b), a x b, α, β > 1 (.1) Construction Management 3 Dr. Emad Elbeltagi
where k is a constant which can be expressed in terms of α and β. Several beta distributions for different sets of values of α and β are shown in Figure.. Figure.: Illustration of several Beta distributions Using beta distribution, simple approximations are made for the activities mean time and its standard deviation. Using the three times estimates, the expected mean time (t e ) is derived using Eq... Then, t e is used as the best available time approximation for the activity in question. The standard deviation is given by Eq..3, and hence the variance (ν) can be determined as ν = σ. a + m + b t e = (.) b a σ = (.3) By adopting activity expected mean time, the critical path calculations proceed as CPM. ssociated with each duration in PERT, however, is its standard deviation or its Construction Management Dr. Emad Elbeltagi
variance. The project duration is determined by summing up the activity expected mean time along the critical path and thus will be an expected mean duration. Since the activities on the critical path are independent of each other, central limit theory gives the variance of the project duration as the sum of the individual variances of these critical path activities. Once the expected mean time for project duration (T X ) and its standard deviation (σ X ) are determined, it is possible to calculate the chance of meeting specific project duration (T S ). Then normal probability tables are used to determine such chance using Equation.. Z T T S X = (.) σ X T = T + Z * σ is an equivalent form of Equation., which enables the scheduled S X X time for an event to be determined based on a given risk level. The procedure for hand probability computations using PERT can be summarized in the following steps: 1. Make the usual forward and backward pass computations based on a single estimate (mean) for each activity.. Obtain estimates for a, m, and b for only critical activities. If necessary, adjust the length of the critical path as dictated by the new t e values based on a, m, and b. 3. Compute the variance for event x (ν X ) by summing the variances for the critical activities leading to event x.. Compute Z using Equation. and find the corresponding normal probability. Consider the nine activity example project shown in Table.1. Suppose that the project have very uncertain activity time durations. s a result, project scheduling considering this uncertainty is desired. Construction Management Dr. Emad Elbeltagi
Table.1: Precedence relations and durations for a 9-activity project example ctivity Description Predecessors Duration B C D E F G H I Site clearing Removal of trees General excavation Grading general area Excavation for trenches Placing formwork and RFT for concrete Installing sewer lines Installing other utilities Pouring concrete --- --- B, C B, C D, E D, E F, G 3 9 1 Table. shows the estimated optimistic, most likely and pessimistic durations for the nine activities. From these estimates, the mean, variance and standard deviation are calculated. In Figure.3, PERT calculations are performed very similar to that of CPM, considering the mean duration of each activity. Table.: ctivity duration estimates for the 9-activity project ctivity Duration Standard Mean a m b Deviation B C D E F G H I 3 3 9 1 1 1.0 3..0. 9.3 1.0.3.3.0 0.33 0.0 0. 0.0 1.33 0. 0.33 0. 0. Construction Management Dr. Emad Elbeltagi
.0 D. 1.9. 1. 1.3 H...3 30.0 0.0.0.0 C 1 1.0 E 1.3 1.3 G 3.3.0 I 30.0 0.0.0.0.0.0 1.0 1. 9.3 1. 1..3.0.0.0 30.0 0.0 B 3. 1.0 F.0. 3. 1.0 1.0 1.0 Figure.3: PERT calculations for 9-activity example The critical path for this project ignoring uncertainty in activity durations consists of activities, C, F and I. pplying the PERT analysis procedure suggests that the duration of the project would be approximately normally distributed. The sum of the means for the critical activities is.0 +.0 + 1.0 +.0 = 30.0 days, and the sum of the variances is (0.33) + (0.) + (0.) + (0.) = 1. leading to a standard deviation of 1. days. With normally distributed project duration, the probability of meeting a project deadline can be computed using Equation (.). For example, the probability of project completion within 3 days is: 3 30 Z = =.1 1. where z is the standard normal distribution tabulated value of the cumulative standard distribution, which can be determined form standard tables of normal distribution. From Table.3, the probability of completing the project in 3 days is 0%. Construction Management Dr. Emad Elbeltagi
-3σ -σ -σ 0 σ σ 3σ Figure.: Normal distribution curve Table.3: rea under the normal curve measured from the center SD rea % from the center SD rea % from the center 0.1σ 0. 0.3 0. 0. 0. 0. 0. 0.9 1.0 1.1 1. 1.3 1. 1..0.9 11. 1. 19.... 31. 3.1 3. 3. 0.3 1.9 3.3 1. 1. 1. 1.9.0.1..3......9 3.0....1....9 9. 9. 9. 9. 9. 9.9 9.99 Example.1 Suppose that a network has been developed for a particular project with nondeterministic durations for the activities and the completion time for that network is 30 days and the sum of the standard deviation for the activities on the critical path is 130. Find the probability that the project will be completed in 300 days. Construction Management Dr. Emad Elbeltagi
Solution First, convert the normal random variable to the standard normal random variable. T S T Z = X = (300 30) /. = - 0.3 σ X From Table 3, the corresponding probability = 1.% Then, the probability to complete the project in 300 days equals = 0 1. = 33.%. Example. Given the information from the previous example, what is the duration that you can give with 90 percent assurance? Solution From tables find the value of z corresponding to probability of 0%, thud yields z = 1. then, apply z into equation : 1. = (t 30) /. or t =. x 1. + 30 = 30 days..1. Criticism to Program Evaluation and Review Technique While the PERT method has been made widely available, it suffers from three major problems. First, the procedure focuses upon a single critical path, when many paths might become critical due to random fluctuations. For example, suppose that the critical path with longest expected time happened to be completed early. Unfortunately, this does not necessarily mean that the project is completed early since another path or sequence of activities might take longer. Similarly, a longer than expected duration for an activity not on the critical path might result in that activity suddenly becoming Construction Management 9 Dr. Emad Elbeltagi
critical. s a result of the focus on only a single path, the PERT method typically underestimates the actual project duration. s a second problem with the PERT procedure, it is incorrect to assume that most construction activity durations are independent random variables. In practice, durations are correlated with one another. For example, if problems are encountered in the delivery of concrete for a project, this problem is likely to influence the expected duration of numerous activities involving concrete pours on a project. Positive correlations of this type between activity durations imply that the PERT method underestimates the variance of the critical path and thereby produces over-optimistic expectations of the probability of meeting a particular project completion deadline. Finally, the PERT method requires three duration estimates for each activity rather than the single estimate developed for critical path scheduling. Thus, the difficulty and labor of estimating activity characteristics is multiplied threefold.. Monte Carlo Simulation..1 Monte Carlo Simulation Characteristics - Replaces analytic solution with raw computing power. - voids need to simplify to get analytic solution - No need to assume functional form of activity/project distributions. - llows determining the criticality index of an activity (Proportion of runs in which the activity was in the critical path). - Hundreds to thousands of simulations needed... Monte Carlo Simulation Process - Set the duration distribution for each activity. No functional form of distribution assumed. Could be joint distribution for multiple activities. - Iterate: for each trial ( realization ) Sample random duration from each distributions Construction Management 1 Dr. Emad Elbeltagi
Find critical path & durations with standard CPM; Record these results - Report recorded results Report recorded results. Duration distribution Per--node criticality index (% runs where critical) Example.3 Table.: Data for Example.3 ctivity Optimistic time, a Most likely time, m Pessimistic time, b Expected value, d Standard deviation, s 1 B 1 3 3 0. C 9 0.33 D 1 E 0.33 F 0. G 0.33 Construction Management 111 Dr. Emad Elbeltagi
Table.: Summary of simulation runs for Example.3 Run ctivity Duration Number B C D E F G 1 3 9.3.1..3..1....0. 1..9.3. 0... 1.1 3.....9...9.9..3..0.0 9.....0.9.3......1.3..9.1...0..3.. 3.0.9 3.1...9.....0.9 3.1 Critical Path -C-F-G -D-F-G -C-F-G -D-F-G -C-F-G -C-F-G -C-F-G -C-F-G -C-F-G -C-F-G Completion Time. 1...0 3.0.3..1 1.9 19. Figure.: Project duration distribution Construction Management 11 Dr. Emad Elbeltagi
Then the probability that a project ends in a specific time (t) equals number of times the project finished in less than or Equal to t divided by the total number of replications. For example, the probability that the project ends in 0 weeks or less is: P( 0 ) = 13 / 0 = %...3 Criticality Index Criticality index is defined as the proportion of runs in which the activity was in the critical path. PERT assumes binary (either 0% or 0%)..3 Exercises 1. Recently, you were assigned to manage a project for your company. You have constructed a network diagram depicting various activities in the project. In addition, you asked various managers and subordinates to estimate the amount of time they would expect each activity to take. Their responses (in days) were as follows: ctivity Optimistic (a) Duration (days) Most likely (m) Pessimistic (b) B 1 C 3 D E 9 C D 0 B E 3 Construction Management 113 Dr. Emad Elbeltagi
a. Compute the mean and variance in time for each activity. b. Determine the critical path and the expected length of the critical path. c. ssume that the time required to complete a path is normally distributed. What is the probability of completing the critical path in less than 1 days? d. If you wanted to be at least 9 percent sure of completing the project on time, what schedule durations would you quote?. Consider the project given in the next table. Find the probability that the project will be completed within 0 and 0 days. What is the project expected duration corresponding to 0% assurance. ctivity Optimistic (a) Duration (days) Most likely (m) Pessimistic (b) Dependencies 1 - B 3 - C 3 0 D 1 1 1 B E 0 3 C, D F 13 1 19, B 3. You have just decided to open your own geotechnical consulting firm. To secure financing from the bank, you must present a plan for when you expect some income from your firm. You sit down with your banker and develop this plan of things that must be accomplished prior to starting your fist job. Estimated duration and precedence relationships for these activities are shown in the table below. Durations are given in days. a. Draw a PERT diagram and determine the expected duration and the critical path. Construction Management 11 Dr. Emad Elbeltagi
b. ssume that the time required to complete a path is normally distributed. Determine the probability of being able to start your fist job within days. c. To provide an allowance for unforeseen problems, you want to present your banker with a time before you start operations based upon a 9 percent confidence level. How many days from receipt of the loan will you tell banker to expect you to start your first job? ctivity Description Optimistic (a) B C D E F G H I Purchase land Hire staff Obtain permits Obtain business license Site preparation Construct office Paving & landscaping Stock soil testing equipment Test equipment 30 9 0 3 1 9 11 Duration (days) Most likely (m) 0 1 30 1 Pessimistic (b) 90 3 1 11 1 1 1 1 Immediate predecessor - C, D E F B, G H Construction Management 11 Dr. Emad Elbeltagi