Research Article  Open Access
Muhammad AlSalamah, "Temporary Workforce Planning with Firm Contracts: A Model and a Simulated Annealing Heuristic", Mathematical Problems in Engineering, vol. 2011, Article ID 209693, 18 pages, 2011. https://doi.org/10.1155/2011/209693
Temporary Workforce Planning with Firm Contracts: A Model and a Simulated Annealing Heuristic
Abstract
The aim of this paper is to introduce a model for temporary staffing when temporary employment is managed by firm contracts and to propose a simulated annealingbased method to solve the model. Temporary employment is a policy frequently used to adjust the working hour capacity to fluctuating demand. Temporary workforce planning models have been unnecessarily simplified to account for only periodic hiring and laying off; a company can review its workforce requirement every period and make hirefire decisions accordingly, usually with a layoff cost. We present a more realistic temporary workforce planning model that assumes a firm contract between the worker and the company, which can extend to several periods. The model assumes the traditional constraints, such as inventory balance constraints, worker availability, and labor hour mix. The costs are the inventory holding cost, training cost of the temporary workers, and the backorder cost. The mixed integer model developed for this case has been found to be difficult to solve even for small problem sizes; therefore, a simulated annealing algorithm is proposed to solve the mixed integer model. The performance of the SA algorithm is compared with the CPLEX solution.
1. Introduction
Temporary employment has been advocated as a means of flexible employment by which companies can effectively adjust the workforce level when faced by fluctuating demand [1, 2]. Most manufacturing companies, such as, and especially, in hitech industry, face the frequently occurring problem of manpower planning. Production planners need to plan the production by determining the delivery dates and invested capital, and one of the most important and most serious planning activities is the planning of the workforce size needed to meet the production demands. Manpower planning is not an easy task given the high costs involved. Most of these companies revert to temporary employment to absorb and accommodate shortterm requirements for manpower. Essentially, manpower planning, among many things, tries to obtain the right mix of companyowned workers who have low setup costs and high variable costs and contracted workers who have high setup costs and low variable costs. The significance of the problem comes from the fact that this problem confronts the hardwareassembly industry, such as the PC assembly companies in Saudi Arabia.
These manufacturing companies have fixed capital investments in production equipment; hence their fixed costs are primarily a function of the workforce size. Shortterm employment entails two expenses a company should be willing to pay: the training cost and nonproductive pay. At the beginning of the contract, contracted workers will go through a training period during which they are not carrying out any productive work. The worker contracts are considered to be firm and the worker has to serve and will be paid for the duration of the contract. The contract covers the training period, which typically lasts for a small fraction of the duration of the contract.
These PC assembly companies hire workers with technical college degrees, who are skilled in the technical matters of assembly of electronic components and wiring. When those workers are hired on temporary contracts, they receive training on the assembly procedures in the plant, and, because of their technical background, the workers pick up the work very quickly and become as skillful as the regular workers. Other tasks, such as packaging and delivery, do not require a technical background, and the workers hired for these tasks are trained on the production procedures of the company, and they can master these tasks in a very short period of time.
A typical PC assembly plant has the design shown in Figure 1. The PC assembly process consists of three main stages. The first stage is the assembly stage, where the components of the PC are assembled and added to the PC case. The finished cases are delivered to the next stage, which is the software installation stage. In this stage, the hardware is tested before the software can be installed. Then, the cases are taken to the final stage, which is packaging. From packaging, the PC units are delivered to storage. The two stages that require technical education are the assembly and software installation stages, which can only be fulfilled by workers with a technical college education.
The forms and types of temporary employment vary across countries and industries according to employment laws and manpower needs and strategies [3]. In the major European countries, the duration of the contract can range from twelve months to an unlimited period, with provisions for repeated contracts [4].
Thelabor laws in the Kingdom have a provision for temporary employment. The laws permit firms and companies to employ workers on a nonpermanent basis for a specific period or for the duration of a particular job. The Labor Law (LABOR LAW, Royal Decree no. M/51, 27 September 2005, https://mol.gov.sa/en/Pages/LaborLaw.aspx, accessed on November 5, 2010) created by the Royal Decree no. M/51 defines temporary work as follows:
Work considered by its nature to be part of the employer’s activities, the completion of which requires a specific period or relates to a specific job and ends with its completion. (Labor Law, 2005, page 7)
The law hence permits companies to engage in temporary work that would require temporary employment. Under this law companies are not limited to a specific number of temporary workers, and the law does not limit the applicability of temporary employment contracts to specific types of jobs. As stated in the Labor Law, the provisions and laws apply to any contract whereby a person commits himself to work for an employer and under his management or supervision for a wage. However, the law under Article 26 of Chapter 1, Part II, puts restrictions on the manner in which workers are recruited, and, in fact, firms are required to have plans for the nationalization of jobs. The article states that firms, regardless of their number of workers, should exert every effort to attract and employ nationals, provide conditions to keep them in the job, and avail them of an adequate opportunity to prove their suitability for the job by guiding, training, and qualifying them for their assigned jobs. The article also requires firms to let the share of nationals in the total employed be no less than seventyfive per cent.
Article 6 of Chapter 2 of Part I of the Labor Law identifies the provisions that are related to temporary workers. The article states that temporary workers are subject to the provisions on duties and disciplinary rules, the maximum working hours, daily and weekly rest intervals, overtime work, official holidays, safety rules, occupational health, work injuries, and compensation as well as other provisions that are decided by the Ministry of Labor. Under this article, temporary workers have the same privileges and responsibilities as regular workers. In terms of workload, a worker may not work for more than eight hours a day on a daily basis or more than fortyeight hours a week on a weekly basis. During the month of Ramadan, the actual working hours for Muslims is reduced to a maximum of six hours a day or thirtysix hours a week (Article 98, Chapter 2, Part VI). In special cases, the number of working hours can be raised to nine hours a day or reduced to seven hours (Article 99, Chapter 2, Part VI).
The law specifies the rest times for workers as a matter of duration and timing. By law, workers should have a break at every five hours at a minimum, and the break should not last for less than thirty minutes. In terms of work duration, the daily workload should not exceed eleven hours. Workers should be compensated for overtime working hours by an additional amount equal to the hourly wage plus fifty per cent of the basic wage (Articles 101 and 107, Chapter 3, Part VI). The law specifies the duration of the temporary employment. Article 2 of Chapter 1, Part 1, declares that temporary employment contracts should not exceed ninety days. Implicit in this law is the possibility of contract renewal.
Traditional periodic workforce review models permit planners to adjust the workforce level every period, with the flexibility to hire temporary workers without longterm commitment; in most cases, a layoff cost is incurred. This flexibility may not always be available to a company; laws, as explained above, or the company’s own policies may impose fixedterm contracts. We offer an alternative formulation of the traditional temporary manpower planning model with fixed contract duration that can extend to several time periods.
2. Literature on Periodic Workforce Review Models
The literature on temporary employment management can be classified as in the following: work commitment (i.e., [4–7]), worker attitude (i.e., [8, 9]), economics of temporary employment (i.e., [10–13]), and worker accident (i.e., [14, 15]). Traditionally, human management decisions have been solved using optimization mathematical models. Recently, many mathematical models have been proposed that would evaluate in economic terms the benefits of temporary employment and find the optimal mix of the regular and temporary workforce. Some notable periodic review models on temporary employment include Berman and Larson [16], Herer and Harel [17], Pinker and Larson [18], and Techawiboonwong et al. [19]. An extensive review of research related to temporary workers is reported by Connelly and Gallagher [7]. Generally, mathematical models are constructed to optimize a measure of benefit or load, subject to workload demand and manpower balance constraints.
Lee and Vairaktarakis [20] have studied a class of serial assembly system lines which are formed by arranging several production cells or stations in series where all stations have the same production cycle. The developed model will find a sequence of job, which minimizes the maximum workforce requirements over all production cycles.
Foote and Folta [10] studied the use of temporary workers as an option for a workforce expansion decision. They concluded that there must be a particular reason for a tradeoff between these two types of workers. Kandel and Pearson [6] studied the optimal balance between permanent employment and temporary contracts between a firm and its workers: the tradeoff between higher productivity on the one hand and flexibility on the other. Through an examination of the optimal policy using an irreversible investment model, they found that higher expected growth in demand implies increases in the proportion of permanent workers, while an increase in the variability of demand lowers their proportion. An increase in the absolute cost advantage of the permanent workers naturally increases their optimal proportion in the labor force.
Wirojanagud et al. [21] developed a mixed integer programming model to determine the amount of hiring, firing, and crosstraining in a staffing decision problem, in which workers are inherently different, to minimize the total costs, which include training, salary, firing, and missed production costs over multiple time periods.
Song and Huang [22] considered the planning problem of transferring, hiring, or firing employees among different departments or branches of an organization under an environment of uncertain workforce demands and turnover, with the objective of minimizing the expected cost over a finite planning horizon. The organization could transfer the employees from one department to another, hire, or fire employees to meet the workforce demands, under the conditions that there may be constraints on the number of employees that could be transferred, hired, or fired in each period.
Techawiboonwong et al. [19] studied temporary employment in a laborintensive industry. To aid the decision makers, an MIP model was proposed to determine the type (skilled or unskilled) and the number of temporary workers to hire or lay off in each period. The skill levels are mandated by the nature of the manufacturing process. In this treatment of the problem, temporary workers can be laid off in any period but there is a cost associated with the layoff. It is assumed that the pool of regular workers is crosstrained, and therefore they are not restricted to being assigned only to specific workstations.
3. A Temporary Employment Model with Extended Contract Durations
We model the temporary employment planning problem with fixedterm temporary employment contracts that can extend to several periods as an MIP model. We argue that firms can lower their workforcerelated costs if they adopt a fixedterm contract policy. The assumptions of the model are illustrated next. Training is required for temporary workers; but we assume, from an observed case, that temporary workers, after training, become as skillful as regular workers and that their performance is comparable to those of the regular workers when considering important manufacturing variables such as quality, scrap rate, and setups. This assumption holds true especially in the case presented earlier of the PC assembly process, in which temporary workers have a technical degree education and, when they are hired, receive training on the internal procedures of the company and the assembly process, which can be mastered very quickly. Therefore, we assume that the training period is negligibly short. For the sake of a wider coverage of the problem described here and all similar problems, we consider a multiproduct, multistage manufacturing/assembly process.
The manufacturing/assembly process under consideration is illustrated in Figure 2. This general process consists of a number of production lines: each line is for a specific product. Each line has a number of manufacturing/assembly stages. Every product line is equipped with an individual inventory. Regular workers are crosstrained and they can be freely assigned to stages and product lines. Each stage of any line can be operated by a pool of regular and temporary workers. While regular workers are crosstrained, temporary workers are hired for specific tasks and stages. This assumption is close to the case of the PC assembly process, in which temporary workers are hired to perform specific jobs. Following the assumption of Techawiboonwong et al. [19], work in process is not considered.
Unlike the model by Techawiboonwong et al. [19], our MIP model does not have the flexibility of overtime allowance. In our manufacturing/assembly process, there are product lines; for each line, there are stages or workstations. Let and be the indices for product and stage, respectively. The length of the planning horizon is , and is the time index. Define the following parameters:: unit processing time for product in stage ,: production time in period ,: labor productivity rating of temporary and regular workers, respectively,: total regular workers,: wage of temporary worker assigned to stage of product ,: training cost of temporary worker assigned to stage of product ,: inventory holding cost of product per period,: backorder cost of product per unit,: inventory level limit for product ,: duration of the contract,: demand forecast for product in time .
The decision variables are as follows:: number of temporary workers available for stage of product in time ,: number of regular workers assigned to stage of product in time ,: quantity of production from product in time ,: number of temporary workers hired in time for stage of product ,: inventory level of product in time ,: backorder quantity of product in time ,: assigned hours of regular workers to product in stage ,: assigned hours of temporary workers to product in stage .
The MIP model for temporary workforce planning with fixedterm contracts, referred to henceforth as TWPSFTC, is as follows.
TWPSFTC: The objective function consists of the temporary employment wages, training costs of temporary workers, inventory holding costs, and backorder costs. Constraint set (3.2) ensures that the number of available temporary workers in time is equal to the number of temporary workers hired in the current period and those whose contracts end in the current period, and all the temporary workers in between. As initial values, for all and . Since hiring a temporary worker whose contract duration is going to extend beyond the end of the planning horizon is not permissible, we require for for all and . Constraint set (3.3) is the set of inventory balance constraints. We impose the restriction that there should not be a backorder in the last period; hence, for all . The inventory capacity constraint is managed by constraint set (3.4). Constraint set (3.5) restricts the production hours to the available temporary and regular worker hours. Constraint sets (3.6) and (3.7) ensure the availabilities of temporary and regular workers to meet the planned production hours, respectively. The level of the regular workforce is imposed by constraint (3.8).
4. Simulated Annealing Solution Approach
The model TWPSFTC has integer variables, where is the cardinality of the set. Even for relatively small problems, MIP models are hard to solve to optimality. The simulated annealing heuristic has been successfully applied to manpower planning, such as that mentioned by Brusco and Jacobs [23, 24], Thompson [25], Ahire et al. [26], Lim et al. [27], Li et al. [28], Seçkiner and Kurt [29], and others. The simulated annealing heuristic that is proposed here to solve TWPSFTC generates an all feasible solution in every run of the heuristic. The SA heuristic has only two free variables to change: assignment set of regular workers ( variables) and the production set ( variables). Therefore, the dependent variables , , , , , and are calculated from the determined sets.
The simulated annealing algorithm proposed to solve TWPSFTC requires four rules. The first rule specifies the generation of the initial assignment set of regular workers and the production set. This rule is required to start the simulated annealing algorithm. The second rule specifies the way the neighborhood solution is generated from the current solution. The third rule specifies the generation of a feasible solution. The fourth rule determines the acceptance criterion of a solution and the termination criterion of the algorithm.
In the following subsections, we will explain how the initial sets are generated randomly and we show how the values of , , , , , and are generated so that the generated solution becomes feasible. The main steps of the SA algorithm can be outlined as follows.(1)Generate a random initial set.(2)Generate a feasible solution.(3)Loop until the stopping criterion is satisfied:(a)generate a neighborhood solution,(b)generate a feasible solution,(c)evaluate the fitness function,(d)compare and accept or reject the new solution.
The fact that the production assignment determines inventory and backorder values [30] will make it possible to generate a feasible solution, as will be explained in the following sections.
4.1. Generation of the Initial Assignment Set of Regular Workers
In each period of the planning horizon, the company has regular workers to assign to the different product lines and stages. The initial assignment set of the regular workers can best be generated randomly. The following steps illustrate the mechanism for the random generation of the assignment set.(1)Start with the first period in the planning horizon.(2)Repeat until all regular workers are assigned:(a)randomly choose a product line and a stage,(b)assign to the chosen product line and stage one regular worker. If regular workers have been assigned to this particular line stage before, add this worker to the existing workers.(3)Advance to the next period, and repeat step (2) until the last period of the planning horizon is covered.
Let and be the randomly determined values of the product line and stage, respectively. Then, and , where is the floor function and is a uniform number generator between and .
4.2. Generation of a Neighborhood Solution from the Assignment Set of Regular Workers
The neighborhood search structure (NSS) generates a new solution from the current candidate solution by slightly modifying its formation. The structure we adopt in this algorithm is of reassignment, or a move operation as in Jin et al. [31]. A neighborhood solution from the assignment set of regular workers is generated from a current set by randomly reassigning a regular worker to a different stage in the same product line or to a different line and stage. The steps involved in generating a neighborhood solution are as follows.(1)Repeat the following steps until a stage is found with regular workers assigned:(a)randomly choose a period in the planning horizon,(b)randomly choose a product line,(c)randomly choose a stage.(2)Reduce the regular workers assigned by one.(3)Randomly choose a product line and stage in the same period.(4)Increase the assigned number of regular workers by one.
These steps are illustrated in Figure 3. The assignment of regular workers can be organized in a matrix; the rows are for the product lines and the columns are for the stages. In this representation, there are matrices equal to the number of periods in the planning horizon. The choice of the period is made randomly by computing . Within this period, random product line and stage are selected by computing and , respectively. If no regular workers are assigned to this stage, a different period is again selected randomly and the process is repeated until a stage with regular workers is found. Then, one regular worker is reassigned to a different product line and stage in the same period.
4.3. Generation of the Initial Set of the Production
The only requirement for a production plan to be feasible is that the sum of the productions in all periods is equal to the sum of demands in all periods for a specific product; mathematically, we have to have for product . As an initial feasible solution, the sum of the demands can be distributed evenly throughout the periods.
4.4. Generation of a Neighborhood Solution from the Production Set
In generating a neighborhood solution from the production set, we apply a load transfer operation. For a random period and random product, a random production quantity is transferred and loaded into another period’s production. Assume that and are generated randomly; generate . Make , generate another ; and make .
4.5. Generation of Inventory and Backorder Sets
The production set determines the inventory and backorder sets. We apply a lookahead mechanism for computing the inventory and backorder. Any backorder quantity represents an additional demand that is put into a future period, and inventory serves as extended production capacity that will be available in a future period.
We apply the following steps to compute the inventory and backorder for a given production set.(1)Initialization: set , and set inventory and backorder to zero for all products and periods.(2)While , do(i)Initialization: (ii)While , do(1)Compute (2)Check(a)If , set (b)Else if , set (c)Else, do not change and (3)Decrement by 1(iii)End while(b) Increment by 1.(3)End while.
4.6. Generation of Regular Labor Hours and Temporary Labor Hours
The regular labor hours are determined from the regular labor assignment set, which is determined by the generation steps in regular worker assignment set. The regular labor hours are computed by reversing constraint (3.7), so that we can write that the regular labor hours in time , product line , and stage are . Constraint (3.5) defines the relationship between the demand hours, regular worker hours, and temporary labor hours, so we can write , for all , , and .
4.7. Generation of Temporary Worker Utilization Plan Set and Temporary Worker Hiring Set
From the temporary worker hours, we can calculate the number of temporary workers required to meet the production demands. The number of temporary workers needed in period , product line , and stage is .
When the temporary worker utilization plan set is determined, the temporary worker hiring set can be determined. In the longterm contract agreement, additional temporary workers are hired only when previously hired temporary workers who are still within the contract validity are not adequate to meet the production demands. Hence, the number of temporary workers that have to be hired in the current period is merely the difference between the required temporary workers in period and the total hired workers in the interval . We can write that for all , , and . It is assumed for convenience that .
4.8. The Fitness Function
In simulated annealing algorithms, the fitness function is usually defined by the objective function of the underlying model, as in the studies by Chen and Shahandashti [32], He et al. [33], and Sayarshad and Ghoseiri [34]. In our proposed simulated annealing algorithm, the fitness function is defined by the objective function in model TWPSFTC.
4.9. Acceptance of a Solution and Termination of the Algorithm
In simulated annealing, a neighborhood solution is accepted in two conditions: the neighborhood solution is better than the best solution in the current iteration or else it is accepted with some probability. If the neighborhood solution is not accepted, it is discarded and a new neighborhood solution is generated. In our simulated annealing algorithm, there are two independent decision variables: the regular worker assignment and production . This situation causes a difficulty because generating a neighborhood solution that is a combination of two generated neighborhood solutions is intractable; the improvement in the fitness function cannot be attributed to a single neighborhood solution; and the improving neighborhood solution cannot be known. Therefore, we propose a twostage procedure that will systematically generate a tractable neighborhood solution.
The twostage procedure calls for the arbitrary selection of either regular worker assignment or production to be in the first stage in the procedure. The first stage will generate a neighborhood solution for the first variable and pass it to the second stage where several neighborhood solutions are generated for the second variable before the execution returns to the first stage where a new neighborhood solution for the first variable is generated. How many neighborhood solutions have to be generated in the second stage is a decision variable and only through experimentation can the best value be determined.
Simulated annealing was formally introduced in 1983 by Kirkpatrick et al. [35] to solve combinatorial optimization problems. It is motivated by an analogy between combinatorial optimization and physical annealing. Metals are first heated to a temperature that permits many molecules to move freely with respect to each other, then they are cooled slowly until the material freezes into a completely ordered form. In combinatorial optimization, simulated annealing techniques use an analogous cooling operation for transforming a poor, unordered solution into an ordered, desirable solution, so as to optimize the objective function.
The outer loop of the simulated annealing algorithm is terminated when a specific number of iterations are consumed. In each iteration of the algorithm, a random neighborhood solution is generated; if its fitness is better than the fitness of the current solution, it is accepted; otherwise it is accepted with a probability that is a function of the annealing temperature.
The literature lists many temperature cooling schemes [36]. In our simulated annealing algorithm, we apply the cooling scheme , where is called the cooling rate, which controls the speed at which the temperature falls [31]. The value given to should be chosen carefully so that the temperature does not drop to zero prematurely; values that are close to zero will make the drop very slow. Again, the choice of can be determined by experimentation.
4.10. The Simulated Annealing Algorithm
Our simulated annealing algorithm starts by randomly generating initial sets of regular workers and production. At the beginning, both the local and the global best solutions are set to the initial solution. A neighborhood solution, based on the outlined procedures in Sections 4.2, 4.4, and 4.9, is generated; this solution is temporarily accepted when it is better then the best local solution, or else it is accepted with some probability, which is influenced by the cooling temperature. The mainstream in simulated annealing research talks about updating the temperature in every run of the algorithm whether there has been a successful move. However, we propose that the temperature is updated only when the local solution has been changed; this mechanism will prevent the temperature from being unnecessarily reduced and prematurely too small. The algorithm is terminated when the loop counter exceeds a prespecified limit.
Define to be the assignment set of regular workers in iteration , to be the production set in step of the twostep procedure as explained in Section 4.9 in iteration , to be the fitness in iteration , to be the fitness of the local solution in iteration , to be the fitness of the global solution, to be the local assignment set in iteration , to be the global assignment set, to be the local production set in iteration , and to be the global production set. The counter limit is set to . The steps of the algorithm are explained in the following procedure.
Simulated Annealing Algorithm
Step 1. Initialization step: set the outer loop counter and cooling temperature .Step 2. Generate initial sets and as outlined in Sections 4.1 and 4.3. Generate the inventory and backorder sets as in Section 4.5, the regular labor hours and temporary labor hours as in Section 4.6, and the temporary worker utilization plan set and temporary worker hiring set as in Section 4.7. Evaluate the fitness of the initial solution ; let ; set ; set .Step 3. Do the outer loop.(a)Set the inner loop counter .(b)Generate a neighborhood assignment set from as in Section 4.2.(c)While , do the inner loop.(i)Set .(ii)Generate a neighborhood production set from as in Section 4.4.(iii)Generate the inventory and backorder sets, the regular labor hours and temporary labor hours, and the temporary worker utilization plan set and temporary worker hiring set. Evaluate the fitness of the solution . Set .(iv)Check If , then accept the solution and set , , and . If , then set , , and , reduce the cooling temperature .(d)End while.Step 4. Set .Step 5. Check(a)If , then set and go back to Step 3.(b)Else terminate the algorithm.
In Section 4.9, we outlined a twolevel procedure to overcome the difficulty of the existence of two independent decision variables in the simulated annealing formulation. In the algorithm, we choose the assignment of regular workers in the first level of the procedure and the production in the second level.
5. Case Study and Computational Analysis
In this section, a temporary staffing plan for a local PC assembly company will be determined using our modeling approach; while the second purpose of this case study is to illustrate the potential of our simulated annealing heuristic in reducing the solution time with an acceptable solution quality. To protect the company’s confidential information, data have been altered. The PC assembly employs a threestage process, consisting of hardware assembly, software installation, and testing. It is assumed three models of PC are to be produced, which differ in the hardware and the software specifications as requested by customers. The monthly demands for the three PC models are assumed as shown in Table 1.

The labor times in minutes for the three stages are shown in Table 2.

Workers operate 7 hours a day for 5 days a week. There are 8,400 total available labor minutes in a month. The worker productivity values for temporary workers are are shown in Table 3.

The regular worker productivity is assumed to be 1. Other parameters have the following values are shown in Table 4.

The wages of the temporary workers in Saudi riyals (1 euro = 5.2 riyals) are estimated to be as shown in Table 5.

The training cost of a temporary worker represents half of the wage of the worker, which accounts for the nonproductive work of the worker during training. The inventory and backorder costs are assumed to be as shown in Table 6.

The parameters of the simulated annealing algorithm have been chosen as shown in Table 7.

Assessment of the Simulated Annealing Algorithm and Parameter Sensitivity Analysis
The simulated annealing algorithm has been coded in the GAMS environment version 22.6. The model was run under the Windows XP operating system with a 3 GHz processor and 1 GH RAM. For the sole purpose of assessing the quality of the solution that is provided by the simulated annealing algorithm, the model TWPSFTC with the data given above is solved using the CPLEX solver. The relative and absolute optimality criteria are set to zero and the solution time is set to 5 hours. The node iteration limit is set to 5,000. The solver terminates prematurely with the best cost equal to SR 107,709.70. The lower bound to the solution, which is given by the relaxed problem, is equal to SR 34,033.53. The solution given by the simulated annealing algorithm has a cost value of SR 110,788, which is within 3% of the solution obtained by the prematurely terminated CPLEX. In terms of solution time, the simulated annealing algorithm terminated in exactly 29 minutes. Since the effort in a single iteration of the algorithm is invariable, the simulated annealing algorithm time cannot be influenced by the choice of the seed of the random number generator.One the most influencing parameters of simulated annealing is the choice of the initial temperature . High initial temperatures will allow the algorithm to accept every randomly generated neighborhood solution; hence, the ability of the algorithm to distinguish between a good and a bad solution is dismantled, which will permit the exploration of a large set of solutions. Therefore, we set the to a relatively high value (100,000); other values have been tested as shown in Figure 4.
The effect of has been tested by varying the values given to it from 1 to 20. When is set to 1, we will allow a simultaneous change in the assignment of regular workers as well as production. The value of gives the best cost value, as shown in Figure 5.
6. Conclusion
This paper proposes a more realistic temporary workforce planning model that assumes an extended contract duration. The state of the art in temporary workforce planning proposes a periodic review of the temporary workforce; which is based on the flexibility to hire and lay off temporary workers in any period of the planning horizon when the company faces fluctuations in the demand for the workforce. This research provides a temporary workforce planning tool for the analysis of the workforce requirement when the company can hire workers on contracts that can extend to several months, and when Saudi Arabia’s labor laws prohibit companies from opting out of the signed contracts. The developed tool is an MIP model that has an objective to minimize the training, wage, inventory, and backorder costs. The model has seven sets of constraints. The model has been inspired by a local PC assembly company. It has been observed that the model is not solvable in a reasonable time to arrive at the optimal solution. Hence, a simulated annealing procedure has been developed, which has been found to be relatively as good as the prematurely terminated CPLEX solver.
Acknowledgment
The author would gratefully acknowledge support received from King Fahd University of Petroleum & Minerals, under Research Grant Project no. IN080399.
References
 D. Francas, N. Löhndorf, and S. Minner, “Machine and labor flexibility in manufacturing networks,” International Journal of Production Economics. In press. View at: Publisher Site  Google Scholar
 J. G. Polavieja, “Flexibility or polarization? Temporary employment and job tasks in Spain,” SocioEconomic Review, vol. 3, no. 2, pp. 233–258, 2005. View at: Publisher Site  Google Scholar
 J. Burgess and J. Connell, “Temporary work and human resources management: issues, challenges and responses,” Personnel Review, vol. 35, no. 2, pp. 129–140, 2006. View at: Publisher Site  Google Scholar
 P. Cahuc and F. PostelVinay, “Temporary jobs, employment protection and labor market performance,” Labour Economics, vol. 9, no. 1, pp. 63–91, 2002. View at: Publisher Site  Google Scholar
 D. G. Gallagher and J. M. Parks, “I pledge thee my troth ... contingently: commitment and the contingent work relationship,” Human Resource Management Review, vol. 11, no. 3, pp. 181–208, 2001. View at: Publisher Site  Google Scholar
 E. Kandel and N. D. Pearson, “Flexibility versus commitment in personnel management,” Journal of the Japanese and International Economies, vol. 15, no. 4, pp. 515–556, 2001. View at: Publisher Site  Google Scholar
 C. E. Connelly and D. G. Gallagher, “Emerging trends in contingent work research,” Journal of Management, vol. 30, no. 6, pp. 959–983, 2004. View at: Publisher Site  Google Scholar
 A. Engellandt and R. T. Riphahn, “Temporary contracts and employee effort,” Labour Economics, vol. 12, no. 3, pp. 281–299, 2005. View at: Publisher Site  Google Scholar
 D. E. Guest, P. Oakley, M. Clinton, and A. Budjanovcanin, “Free or precarious? A comparison of the attitudes of workers in flexible and traditional employment contracts,” Human Resource Management Review, vol. 16, no. 2, pp. 107–124, 2006. View at: Publisher Site  Google Scholar
 D. A. Foote and T. B. Folta, “Temporary workers as real options,” Human Resource Management Review, vol. 12, no. 4, pp. 579–597, 2002. View at: Publisher Site  Google Scholar
 A. L. Kalleberg, J. Reynolds, and P. V. Marsden, “Externalizing employment: flexible staffing arrangements in US organizations,” Social Science Research, vol. 32, no. 4, pp. 525–552, 2003. View at: Publisher Site  Google Scholar
 J. K. Stratman, A. V. Roth, and W. G. Gilland, “The deployment of temporary production workers in assembly operations: a case study of the hidden costs of learning and forgetting,” Journal of Operations Management, vol. 21, no. 6, pp. 689–707, 2004. View at: Publisher Site  Google Scholar
 J. Friesen, “Statutory firing costs and layoffs in Canada,” Labour Economics, vol. 12, no. 2, pp. 147–168, 2005. View at: Publisher Site  Google Scholar
 M. Guadalupe, “The hidden costs of fixed term contracts: the impact on work accidents,” Labour Economics, vol. 10, no. 3, pp. 339–357, 2003. View at: Publisher Site  Google Scholar
 B. Fabiano, F. Currò, A. P. Reverberi, and R. Pastorino, “A statistical study on temporary work and occupational accidents: specific risk factors and risk management strategies,” Safety Science, vol. 46, no. 3, pp. 535–544, 2008. View at: Publisher Site  Google Scholar
 O. Berman and R. C. Larson, “Determining optimal pool size of a temporary callin work force,” European Journal of Operational Research, vol. 73, no. 1, pp. 55–64, 1994. View at: Google Scholar
 Y. T. Herer and G. H. Harel, “Determining the size of the temporary workforce: an inventory modeling approach,” Human Resource Planning, vol. 21, pp. 20–32, 1998. View at: Google Scholar
 E. J. Pinker and R. C. Larson, “Optimizing the use of contingent labor when demand is uncertain,” European Journal of Operational Research, vol. 144, no. 1, pp. 39–55, 2003. View at: Publisher Site  Google Scholar
 A. Techawiboonwong, P. Yenradee, and S. K. Das, “A master scheduling model with skilled and unskilled temporary workers,” International Journal of Production Economics, vol. 103, no. 2, pp. 798–809, 2006. View at: Publisher Site  Google Scholar
 C. Y. Lee and G. L. Vairaktarakis, “Workforce planning in mixed model assembly systems,” Operations Research, vol. 45, no. 4, pp. 553–567, 1997. View at: Google Scholar
 P. Wirojanagud, E. S. Gel, J. W. Fowler, and R. Cardy, “Modelling inherent worker differences for workforce planning,” International Journal of Production Research, vol. 45, no. 3, pp. 525–553, 2007. View at: Publisher Site  Google Scholar
 H. Song and H. C. Huang, “A successive convex approximation method for multistage workforce capacity planning problem with turnover,” European Journal of Operational Research, vol. 188, no. 1, pp. 29–48, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 M. J. Brusco and L. W. Jacobs, “A simulated annealing approach to the cyclic staffscheduling problem,” Naval Research Logistics, vol. 40, no. 1, pp. 69–84, 1993. View at: Google Scholar
 M. J. Brusco and L. W. Jacobs, “A simulated annealing approach to the solution of flexible labour scheduling problems,” Journal of the Operational Research Society, vol. 44, no. 12, pp. 1191–1200, 1993. View at: Google Scholar
 G. M. Thompson, “A simulatedannealing heuristic for shift scheduling using noncontinuously available employees,” Computers and Operations Research, vol. 23, no. 3, pp. 275–288, 1996. View at: Publisher Site  Google Scholar
 S. Ahire, G. Greenwood, A. Gupta, and M. Terwilliger, “Workforceconstrained preventive maintenance scheduling using evolution strategies,” Decision Sciences, vol. 31, no. 4, pp. 833–859, 2000. View at: Google Scholar
 A. Lim, B. Rodrigues, and L. Song, “Manpower scheduling with time windows,” in Proceedings of the ACM Symposium on Applied Computing, pp. 741–746, Melbourne, Fla, USA, March 2003. View at: Google Scholar
 Y. Li, A. Lim, and B. Rodrigues, “Manpower allocation with time windows and jobteaming constraints,” Naval Research Logistics, vol. 52, no. 4, pp. 302–311, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 S. U. Seçkiner and M. Kurt, “A simulated annealing approach to the solution of job rotation scheduling problems,” Applied Mathematics and Computation, vol. 188, no. 1, pp. 31–45, 2007. View at: Publisher Site  Google Scholar
 E. V. Denardo, Dynamic Programming: Models and Applications, Dover Publications, Mineola, NY, USA, 2003.
 F. Jin, S. Song, and C. Wu, “A simulated annealing algorithm for single machine scheduling problems with family setups,” Computers and Operations Research, vol. 36, no. 7, pp. 2133–2138, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 P. H. Chen and S. M. Shahandashti, “Hybrid of genetic algorithm and simulated annealing for multiple project scheduling with multiple resource constraints,” Automation in Construction, vol. 18, no. 4, pp. 434–443, 2009. View at: Publisher Site  Google Scholar
 Z. He, N. Wang, T. Jia, and Y. Xu, “Simulated annealing and tabu search for multimode project payment scheduling,” European Journal of Operational Research, vol. 198, no. 3, pp. 688–696, 2009. View at: Publisher Site  Google Scholar
 H. R. Sayarshad and K. Ghoseiri, “A simulated annealing approach for the multiperiodic railcar fleet sizing problem,” Computers and Operations Research, vol. 36, no. 6, pp. 1789–1799, 2009. View at: Publisher Site  Google Scholar
 S. Kirkpatrick, C. D. Gelatt Jr., and M. P. Vecchi, “Optimization by simulated annealing,” Science, vol. 220, no. 4598, pp. 671–680, 1983. View at: Google Scholar
 B. Naderi, M. Zandieh, A. K. Balagh, and V. Roshanaei, “An improved simulated annealing for hybrid flowshops with sequencedependent setup and transportation times to minimize total completion time and total tardiness,” Expert Systems with Applications, vol. 36, no. 6, pp. 9625–9633, 2009. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2011 Muhammad AlSalamah. 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.