American Journal of Operational Research
p-ISSN: 2324-6537 e-ISSN: 2324-6545
2015; 5(3): 47-56
doi:10.5923/j.ajor.20150503.01
Zohar Laslo, Gregory Gurevich
Faculty of Industrial Engineering and Management, SCE – Shamoon College of Engineering, Beer Sheva, Israel
Correspondence to: Gregory Gurevich, Faculty of Industrial Engineering and Management, SCE – Shamoon College of Engineering, Beer Sheva, Israel.
| Email: | ![]() |
Copyright © 2015 Scientific & Academic Publishing. All Rights Reserved.
The prevalent PERT/CPM, PERT/cost, and CPM time-cost trade-offs procedures are the backbone of current planning and controlling routines, but they require definition of critical paths. Such routines grossly mislead the project manager into thinking that the probability confidence level of attaining his time target is very good, when in reality it is very poor. Routines based on the expected lengths of the project paths also cause wasteful budget investments. We propose a stochastic approach using Monte Carlo simulations for planning and controlling projects with renewable resources. Simulation results show that the proposed planning and controlling routine may well enable the attainment of the project time target, and as much as possible, the cost target as well. This is particularly so in cases where prevalent planning and controlling routines based on so-called critical paths cannot produce similar results.
Keywords: Project Management, Critical Path, Uncertainty, Time-cost Trade-offs, Monte Carlo Simulations
Cite this paper: Zohar Laslo, Gregory Gurevich, Planning and Controlling Projects under Uncertainty, American Journal of Operational Research, Vol. 5 No. 3, 2015, pp. 47-56. doi: 10.5923/j.ajor.20150503.01.
denote the vector of budgets allocated to each of the
project's activities,
,
, where
and
correspond to the
th activity normal execution-mode,
, and crash execution-mode,
, respectively. The assessment of the expected project completion time related to
, by PDM, relies on single estimates, most being expected activity durations,
. PDM assumes that
equals the expected length of the critical path (CP),
, that is
, where
denotes the length of the path
. But where uncertainty is involved, as originally considered by PERT, network paths
with expected length
may be longer than
, and by that determine
. Wherefore, the assessment of
should take into consideration that there is a probability of additional project paths to determine
, thereby enabling participation in the determination of
. Ignoring the distributions of the activity durations and dealing with their expected durations as deterministic values is naïve. Such ignorance causes erroneous consideration that
, while it is obvious that![]() | (1) |
, is identical to the CP variance,
(see [9; 10; 21]). But it is obvious that this statement is groundless because
is not necessarily identical to
.Despite the fact that uncertainty is the crux in achieving the project time target, practitioners mostly implement PDM for assessing
and
. Present views throughout the project lifecycle that rely on definition of CPs cause errors of being certain that the time target is attainable. In point of fact, this is actually incorrect. PDM is thereby inadequate for assessing
because it may hinder exposing deviations from the planned trajectory while there is still time to take corrective actions.
when the initial planning or any present view throughout the project lifecycle indicates an unperformable time or cost target. As previously noted, these procedures are based on definitions of CPs that provide an overly optimistic assessment of
. Consequently, it is naïve to consider that shortening
by one unit of time assures that
will be shortened by one entire unit. Thus, inappropriate assurances that the project can be delivered on time, when this is not the case often leads to premature halt of the PERT/Cost procedure. The attainment of the time target (desired lead-time),
, may also be impeded earlier by a premature halt that derives from the procedure's dependency on the CP concept. Because the PERT/Cost procedure considers that
is determined by
, when all the durations of activities lying on one of the CPs have been shortened up to
investing additional budget for crashing durations of more activities seems to be disadvantageous; whereupon, at this point the crashing procedure is halted. But, crashing durations of more activities that are not lying on the alleged CP may reduce
, as demonstrated by Table 1.
|
is 14. The project is composed of two paths, each with two probable lengths. The expected normal and crash lengths of Path j are 16 and 15 respectively, and of Path k, 15 and 7. Following PERT/Cost, at the initial state
where Path j is identified as CP, which determines
as 16, while the realistic calculation shows that it is 18.5. Because Path j is recognized as CP,
was chosen to be shortened first by one unit of time. PERT/Cost assumes that
was crashed by one unit of time, but in fact it was crashed only by a half unit. The naïve assumption of the PERT/Cost procedure is that at this stage
is 15,
is satisfied, and thus, the coordination is accomplished. However, the realistic
is 18 and the crashing process should be continued. At this stage,
disables additional shortening of
related to Path j, which is considered as decisive CP, and this is what causes a premature halt of the PERT/Cost procedure. But, as we can see from Table 1, by continuing the shortening of
related to Path k, which is not defined as CP, the realistic
can be crashed further.It is clear that the naïve assessment of
and the premature halt of the PERT/Cost procedure disables full utilization of the coordination potential. Furthermore, as explained henceforth (Section 4.1.), the PERT/Cost procedure that omits uncertainty does not take into consideration that distributions of the activity durations may be affected by the crashing process, and by that, change the expected durations. This phenomenon may undermine the logic of the PERT/Cost procedure. These handicaps are troublesome, especially in complex networks where there are many parallel paths. Therefore, the deterministic techniques PERT/Cost and CPM time-cost trade-offs procedures seem to be inadequate when used as the basis of coordination routines.
. Many analytical approximations of the distribution of project completion time,
, have been developed (e.g., [4, 5, 7, 20, 24, 25, 29]), but presumably it is infeasible to compute exact results via analytic algorithms. When analytical algorithms are infeasible, Monte Carlo simulations (MC) tend to be used. MC enables receiving unbiased estimates of project completion time distribution, where other methods cannot.An upgraded version of a model for activity durations and costs under uncertainty [16] is introduced first. We next propose a simple method for reassessing duration and cost estimates for ongoing activities at inspections throughout the controlling process, and then, our proposal for a planning and controlling routine via MC is presented.
can be formulated as follows:• Each
with the target-effective execution duration,
, requires the allocation of budget
, where
and
are known values. • The random duration,
, related to
, is composed of the following random time components:ο The random effective execution duration,
, related to
and affected by ITU: ![]() | (2) |
versus
, are considered as pre-given continuous functions, with their estimated edge points corresponding to
and
respectively.-
is a random component that reflexes ITU; this variable is related to the relevant . It has zero expectation and a standard deviation proportional to
;
, where
is a known standard deviation of
related to
.ο The random wasted time caused by disturbances,
:![]() | (3) |
is considered as a known positive value.-
is a random component that reflexes ETU; the distribution of this variable does not depend on
. It has zero expectation and a known standard deviation
.Thus, the random duration related to
is defined as: ![]() | (4) |
, related to is composed of the following components:ο
, associated with
. ο A random variable that reflexes DDICU,
, is related to
. It has zero expectation and a standard deviation proportional to
,
, where
is a known standard deviation of DDICU related to
. ο A random component that reflexes DDDCU equal to the deviation from the target activity duration,
, multiplied by the activity execution cost-per-time unit related to
,
;
, where
is a known execution cost-per-time unit related to
. Thus, the random cost related to
is defined as:![]() | (5) |
,
and
, the generalized model for activity duration and cost enables obtaining the simulated activity time and cost performances related to each
.
and
related to known
and the evidence that this activity is still unaccomplished. The cumulative distribution functions of
and
(
and
respectively), directly follow from equations (4-5),
,
,
, and initial estimated distributions of
,
, and
. Actually, at inspection the additional information concerning the ongoing activity indicates that
is larger than the known value
, at the moment of the inspection. Thus, at inspection the distribution of the random duration of the ongoing activity should be modified using the information that
. That is, the adjusted ongoing activity duration is the conditional random variable
, with the cumulative distribution function
, where![]() | (6) |

with the cumulative distribution function
, which can be straightforwardly calculated using the form of the cumulative distribution function of the adjusted ongoing activity duration,
. In particular, assuming normal distributions for
,
, and
(i.e., normal distributions for duration and cost of the project activities), the cumulative distribution function of an activity duration has the form![]() | (7) |
is the cumulative distribution function of the standard normal distribution. Then, the cumulative distribution function of an adjusted active activity duration (given
) is defined as![]() | (8) |
defined by equation (8), we note that
, where
is the inverse function of the cumulative distribution function
, and
is a random variable uniformly distributed on the interval
(
). Therefore, by equation (8),![]() | (9) |
![]() | (10) |
is the inverse function of the cumulative distribution function of the standard normal distribution
. Thus, to generate an observation from
one can generate an observation from the uniform distribution
and use equation (10).
. The proposed planning and controlling routine is based on a previously developed algorithm [16], in which the activity is prioritized in receiving an additional unit of budget if it leads to maximal shortening of the
simulated fractile of the project completion time,
. This algorithm was found superior to an alternative algorithm where the priority for receiving additional budget was defined by the activity criticality index (ACI), i.e., the approximated by MC probability of the activity to lie on CP [8]. Assuming some distributions for
,
,
and the generalized model (4-5) at the initial planning and at each inspection, the proposed procedure is composed of the following steps:
Regarding the number of MC repetitions, we clarify that
is obtained as a mean of 10,000 generated
,
is obtained as such value that 9,500 generated
are less than or equal to it and 500 generated
are greater than or equal to it. Note that the obtained estimator for
based on 10,000 repetitions is stable (the approximate confidence interval for
is
and for actual projects it is reasonable to assume that
. The obtained estimator for
is more problematic because the approximate confidence interval for
is
,
, and
is an unknown density function of
(see [26]). Thus, for
this estimator might not be stable enough, in which case one must consider a more essential number of repetitions. In our study we consider several
,
, and use 10,000 repetitions.
, or any cost fractile with probability confidence level,
,
, subject to a pre-determined
or each pre-determined
. A what if? analysis on a large project (30 activities; 160 paths) with renewable resources was performed with the aim of evaluating the possible consequence of implementing alternative planning and controlling routines for delivering projects on time, and as much as possible, within budget. Assuming the generalized model (4-5) for durations and costs of the project activities as well as normal distributions for
,
,
, project execution was simulated and accompanied by control routines based on PDM and on the proposed routine, with several pre-determined
. The simulation results for the project, as shown in Table 2, derive from a scanty budget that disables accomplishing the project on time within budget. But the proposed control routine for planning and controlling, even with modest demand for
, enables the project accomplishment on time with cost overflows, while routines based on PERT/Cost cannot consume the available budget for on time project delivery. Naïve project completion assessments of PDM cause a late response to deviation from the planned trajectory, which requires drastic crashing of project completion. Drastic completion crashing, while the number of the remaining not-yet-started activities is small, mostly invokes a premature halt of the crashing procedure. Due to the fact that
is quite ample, time re-establishment obstructions caused by crash execution-modes do not occur at the beginning of project execution, and as
is lower, such obstructions occur later, if ever. Later, when the present view indicates attainable
, the focus is turned to the cost objective, i.e., accomplishing the project on time with minimal cost. At this stage,
is re-crashed in order to cut costs, but this re-crashing is obstructed by the situation of
in all not-yet-started activities..
|
, but this conclusion is valid only for this specific scenario. We can also conclude that more generous budgeting will enable the accomplishment of the project on time within budget with the proposed routine, while routines based on so-called critical paths cannot accomplish the desired result on time, with any budgeting.