key: cord-0768964-5glpyvy4
authors: Vahedi-Nouri, Behdin; Tavakkoli-Moghaddam, Reza; Hanzálek, Zdeněk; Dolgui, Alexandre
title: Workforce planning and production scheduling in a reconfigurable manufacturing system facing the COVID-19 pandemic
date: 2022-04-28
journal: J Manuf Syst
DOI: 10.1016/j.jmsy.2022.04.018
sha: 4caf76136757596656adfc1c14979a8c09c665f7
doc_id: 768964
cord_uid: 5glpyvy4

Due to the outbreak of the COVID-19 pandemic, the manufacturing sector has been experiencing unprecedented issues, including severe fluctuation in demand, restrictions on the availability and utilization of the workforce, and governmental regulations. Adopting conventional manufacturing practices and planning approaches under such circumstances cannot be effective and may jeopardize workers’ health and satisfaction, as well as the continuity of businesses. Reconfigurable Manufacturing System (RMS) as a new manufacturing paradigm has demonstrated a promising performance when facing abrupt market or system changes. This paper investigates a joint workforce planning and production scheduling problem during the COVID-19 pandemic by leveraging the adaptability and flexibility of an RMS. In this regard, workers' COVID-19 health risk arising from their allocation, and workers' preferences for flexible working hours are incorporated into the problem. Accordingly, first, novel Mixed-Integer Linear Programming (MILP) and Constraint Programming (CP) models are developed to formulate the problem. Next, exploiting the problem’s intrinsic characteristics, two properties of an optimal solution are identified. By incorporating these properties, the initial MILP and CP models are considerably improved. Afterward, to benefit from the strengths of both improved models, a novel hybrid MILP-CP solution approach is devised. Finally, comprehensive computational experiments are conducted to evaluate the performance of the proposed models and extract useful managerial insights on the system flexibility.

Being a long-term crisis, the COVID-19 pandemic has imposed a massive worldwide disruption on the manufacturing sector and supply chains [1, 2] . Accordingly, global manufacturing recorded severe declines of 6%, 11.1%, and 1.1% in the output during the first three quarters of 2020 compared to the same periods in 2019 [3] . In addition to the demand volatilities, several workforce-related issues, such as workers' absence due to the infection and adhering to health regulations, enforce restrictions on workforce availability and utilization in manufacturing shop floors [4] . Having registered 1,438 epidemics between 2011 and 2018, the World Health Organization (WHO) pointed out that besides their death tolls, epidemics can trigger enormous threats to the world's economies [5] . Thus, it is critical to provide manufacturing systems with Based on these issues, a new Mixed-Integer Linear Programming (MILP) model is devised, in which multiple decisions are made simultaneously regarding the employment and allocation of workers to machines, reconfigurations of machines, and scheduling of jobs on machines. Due to the promising performance of the Constraint Programming (CP) approach, especially in solving scheduling problems [21] , it is also employed to model the problem. Afterward, two intrinsic properties of an optimal solution are identified and incorporated in the proposed MILP and CP models to improve their performances. Furthermore, to find a higher quality solution within a shorter run time, a hybrid model is devised to benefit from the advantages of both improved MILP and CP models. Finally, comprehensive computational experiments are conducted to evaluate the performance of the developed models and provide relevant managerial insights.

The rest of the paper is organized as follows. In the next section, some related studies are reviewed, and the contributions of this study are highlighted. In Section 3, the integrated workforce planning and production scheduling problem is elaborated, and the corresponding MILP model is presented. The CP model is introduced in Section 4. In Section 5, two properties of the problem are theoretically proved, and then, the MILP and CP models are accordingly improved, and finally, a hybrid model is developed. Section 6 concentrates on the computational experiments, and finally, conclusions and future research directions are given in Section 7.

Two streams of the literature are reviewed in this section. First, the current progress in integrated workforce planning and parallel machine scheduling problems is investigated. Second, the literature on reconfigurable machine scheduling problems is explored. Afterward, the research gaps are identified and the contributions of the papers are stated.

Activity scheduling and workforce planning are two critical and highly interrelated problems in many contexts including healthcare, project planning, and more specifically, manufacturing. Ignoring their connection can lead to poor or sub-optimal results [22] . Recently, some studies investigated these problems in an integrated manner. For example, Van Den Eeckhout et al. [23] integrated the project scheduling problem with the personnel staffing problem. Regarding the literature on parallel machine scheduling, the incorporation of the workforce can be classified into three categories that are described here.

Studies in this category aim to jointly allocate workers to machines and schedule jobs on the machines so that the workers will not be reallocated during the planning horizon. Thus, the number of workers is equal to or greater than the number of machines. In this regard, Hu [24] studied an identical parallel machine scheduling problem to minimize the total tardiness where the processing time of a job depends on the number of allocated workers to the associated machine. He developed a two-phase heuristic method to solve the problem, in which the first phase is related to job scheduling, and the second phase performs worker allocation. Later, Chaudhry and Drake [25] provided a Genetic Algorithm (GA) for a similar problem. Two metaheuristic algorithms, namely GA and particle swarm optimization (PSO), were developed by Lee et al. [26] to solve a uniform parallel machine problem with the makespan objective. In their problem, the learning effect was considered, in which the learning rates of workers are different. Zouba et al. [27] investigated a particular case for parallel machine scheduling, where the J o u r n a l P r e -p r o o f workers simultaneously supervise a few neighboring machines. They supposed that there are some predefined work modes for worker allocation, and the processing times of the jobs depend on the selected mode. A heuristic method was devised for the problem to determine the work mode in each period and the scheduling of jobs to minimize the makespan.

In this category, besides job scheduling on machines, a worker is determined for each job, and the number of workers is typically less than the number of machines. Thus, a worker can visit another machine after performing his/her current job. Depending on the automation level of the manufacturing system, a worker performs the whole operation of a job or a portion of it (e.g., setting up machines before performing operations, loading or unloading the jobs on machines, and cleaning and testing the jobs). Accordingly, Celano et al. [28] developed a simulationoptimization method for an unrelated parallel machine scheduling, in which a worker is needed to just set up a machine before performing a job. In their study, the workers are homogenous (i.e., their qualification and skill level are the same), and several objective functions were considered, including the total tardiness, makespan, and machine and worker workload. Lee and Jang [29] investigated a uniform parallel scheduling problem with job splitting on dedicated machines and a limited number of homogenous operators for setting up machines. They formulated the problem as an MILP model and provided a constructive heuristic for large-sized problems. Yepes-Borrero et al. [30] provided an MILP and three Greedy Randomized Adaptive Search Procedure (GRASP) algorithms for unrelated parallel machine scheduling to minimize the makespan with sequence and machine-dependent setup times and limited homogenous workers for setting up machines. Later, Yepes-Borrero et al. [31] extended this problem to a biobjective MILP model that minimizes the makespan and the number of workers and devised an iterated Pareto greedy algorithm to find Pareto solutions. Fang et al. [32] incorporated the limitation on the maximum consecutive working time, total working time, and minimum rest time of workers to an integrated scheduling and workforce allocation problem in an unrelated parallel machine setting. They developed an MILP model and a combinatorial Benders decomposition algorithm to minimize the machines' deprecation cost, workers-related costs, and the cost associated with the makespan. Recently, Zhang et al. [33] investigated unrelated parallel machine scheduling, in which workers are responsible for setting up machines. The setups are sequence and machine-dependent and the duration of a setup on a machine is influenced by the worker performing the setup and his/her learning procedure on the machine. They provided a combinatorial evolutionary algorithm to solve the problem.

Other studies have been considered the workforce just as an additional resource with limited capacity. Thus, there is no decision regarding workforce allocation, and only the workers' requirements for jobs are respected. In this line, Edis and Ozkarahan [34] probed an identical parallel machine scheduling problem with eligibility restriction, in which each job may require an entire worker, half of a worker, or no worker to be processed. They provided two solution methods based on a combination of two different Integer Programming (IP) models, and a combination of IP and CP models, in which the IP model determines the assignment of jobs to machines, then the scheduling of jobs on each machine is specified by either the IP or CP model. Afzalirad and Shafipour [35] also worked on this problem and presented a genetic algorithm (GA) to solve it. Lee et al. [36] and Lee and Kim [37] studied a special type of this problem J o u r n a l P r e -p r o o f where a job can be split into multiple parts and processed concurrently on different machines. In their works, workers are assumed as a limited resource to perform setups on machines. Accordingly, they provide MILP models, and heuristic algorithms to obtain decent solutions.

Since its introduction, lots of studies have been performed on different aspects of RMSs, including conceptual design, system design, performance evaluation, and production and process planning. In this regard, comprehensive reviews were provided by Bortolini et al. [8] , Yelles-Chaouche et al. [9] , and Koren et al. [38] . However, a limited number of studies explored scheduling problems in RMSs. Due to the importance of scheduling on operational decisions and consequently on the performance and profitability of RMSs, reconfigurable machine scheduling problems have recently received a lot of attention. Accordingly, Bensmaine et al. [39] presented a heuristic for an integrated process planning and scheduling problem in a shop floor consisting of RMTs. In this problem, each operation of a job can be processed with different machineconfiguration combinations and there are precedence constraints among the operations of a job. Dou et al. [40] devised a multi-objective PSO algorithm for a cyclic process planning and scheduling problem in a reconfigurable flow line with the objectives of the total capital and reconfigurable cost, and the total tardiness. In this problem, the configuration of the system (i.e., the number and type of the machines for each stage of the flow line) and RMTs are determined and fixed for the entire planning period. Mahmoodjanloo et al. [41] explored a flexible job shop environment with multiple RMTs in which a reconfiguration on a machine can be implemented between every two consecutive operations. They developed two MILP models and a Self-Adaptive Differential Evolution Algorithm (SADEA) to solve the problem considering the makespan as the objective function. Later, Mahmoodjanloo et al. [42] devised a self-adaptive hybrid Equilibrium Optimizer (EO) algorithm for a dynamic and multi-factory extension of this problem with the objective of the total weighted lateness. Fan et al. [43] also studied the scheduling problem of a reconfigurable flexible job shop with the objective of total weighted tardiness. They developed an MILP model and an improved GA for the problem assuming that the number of auxiliary modules to provide different configurations on RMTs are limited.

Although integrated workforce planning and production scheduling has widely been investigated, there is a lack of research to assimilate restrictions and health concerns arising from an infectious disease like the COVID-19 into the problem. Thus, the existing approaches in the literature are not applicable during such crises and may jeopardize workers' health and satisfaction as well as the continuity of businesses. Moreover, despite the substantial impact of the workforce on the performance of all manufacturing systems, to the best of the authors' knowledge, no study has integrated workforce-related decisions and human factors in RMS scheduling problems.

To address these gaps, this study explores an integrated workforce planning and production scheduling during the COVID-19 pandemic leveraging the adaptability of an RMS. The main contributions of this study can be summarized below:

 Investigating an integrated workforce planning and production scheduling problem in an RMS consisting of several non-identical parallel RMTs addressing workers' COVID-19 health risk arising from their allocation and workers' preferences for flexible working hours.

 Formulating novel MILP and CP models for the problem, which concurrently makes decisions regarding workforce allocation, reconfigurations of machines, and job scheduling.  Identifying two properties of an optimal solution for the problem and incorporating them in the proposed MILP and CP models to enhance their performance.  Developing a hybrid MILP-CP model to take advantage of both models' strengths.  Conducting comprehensive computational experiments to evaluate the performance of the proposed solution methods and provide some managerial insights as well.

In the following subsections, first, the investigated problem is stated, and its characteristics and assumptions are described. Afterward, the notations and corresponding MILP model are provided.

In this paper, the integrated production scheduling and workforce planning in an RMS is investigated. This system consists of a set of non-identical parallel RMTs. Moreover, there is a set of heterogeneous workers (i.e., workers with different capabilities in operating a set of specific machines) to be allocated to the machines. The planning horizon is assumed to be one day, in which a set of jobs must be performed. The main assumptions of the problem are listed below.

 Each worker is capable of working on a specific subset of machines.  At most, one worker is allocated to a machine during the day, and likewise, a worker is allocated to at most one machine.  The RMTs are non-identical, and their functionality and speed can be changed through a reconfiguration. The reconfiguration time on a machine depends on the current and new configurations, and the reconfiguration is performed by the allocated worker. Moreover, the triangle inequality holds for the reconfiguration times.  Each job can be performed by at least one machine-configuration combination. In other words, some configurations of some machines may be eligible to perform a specific job. On the other hand, each configuration of a machine may be eligible for performing more than one job.  The processing time of a job depends on the assigned machine-configuration combination and the related allocated worker.  Each machine cannot perform more than one job at any moment, and the processing of a job cannot be preempted.  Each machine has a COVID-19-related safety risk level that can be determined by an occupational safety and health expert based on the proximity of a machine to other machines and corridors, the presence of barriers, and the ventilation condition.  There is a COVID-19-related vulnerability level associated with each worker. This can be evaluated by the existing measure (like ALAMA), which is a comprehensive web-based tool provided by the Association of Local Authority Medical Advisors [17] .  Regarding the occupational regulation, a maximum allowed working time is considered for each worker.

 Complying with the curfew regulation, a maximum allowed end time of the day is respected.  Workers' priorities to enter into the system at different arrival times are regarded (a lower value of worker priority is more favorable).  Regarding governmental regulations, a limit on the number of workers entering into the system at the same arrival time is established.  Each worker needs to sanitize his/her allocated machine, table, and equipment at the beginning of his/her work.

This problem aims to provide a plan to fulfill the demand, by considering the issues arising from the COVID-19 pandemic. Accordingly, the following decisions are simultaneously made in the problem.

 Allocation of workers to the machines and determination of their arrival times. It should be noted that considering the reduction of the demand during the COVID-19 pandemic and the vulnerability level of workers, the presence of all workers may not be necessary. Thus, the decision on the presence of each worker is also made in the problem.  Assignment of jobs to machines  Determination of a configuration to perform each job on its assigned machine.  Scheduling of jobs on their assigned machines.

Three objective functions are followed in the problem, which gets normalized and integrated into a single objective. The first objective minimizes the workers' COVID-19 health risk (briefly COVID-19 risk), which is the total product of the vulnerability level of each worker and the safety risk level of his/her allocated machine. The second objective respects workers' preferences, which minimizes the total workers' priorities based on their allocated arrival times. Moreover, the third one is the minimization of the maximum completion time (i.e., makespan). with the highest vulnerability level is not allocated, and machine 2 with the highest safety risk level is not utilized. Worker , who arrived at 8 AM, is allocated to machine . After sanitization, jobs , and are processed on machine with configuration , then the machine is reconfigured to configuration (with the reconfiguration time equal to ), and then jobs and are performed. In Subsections 3.2 and 3.3, the notations and the corresponding MILP model are presented.

Sets:

Set of jobs Set of machines Set of workers Set of configurations for machine Set of times

Index of job Index of machine Index of worker Index of configuration Index of time

Processing time of job performed by worker on machine with configuration Covid-19 vulnerability level associated with worker Covid-19 safety risk level associated with machine Reconfiguration time on machine from configuration to 1 if machine with configuration is eligible to process job ; 0, otherwise 1 if worker is capable of operating machine ; 0, otherwise Arrival time (i.e., the moment that a relevant worker enters into the system) Sanitization time Maximum portion of the number of workers that are allowed to enter into the system at each arrival time Maximum allowed end time of the day Worker 's priority to enter into the system at arrival time ; a lower value is more favorable for a worker Maximum 

Subject to:

J o u r n a l P r e -p r o o f

Equation (1) minimizes the integrated objective function. The first term is related to the COVID-19 risk, the second is related to the workers' preferences, and the third one is the makespan. The terms are normalized by dividing them to their maximum possible values, and then, they are integrated considering the weight of each term. Constraints (2) and (3) are related to the assignment of jobs to machine-configuration combinations respecting the eligibility restriction. Constraint (4) indicates that a worker can be allocated to at most one machine-arrival time combination. It also means that the model may decide a worker be absent. The restriction for workers' capability is regarded by Constraint (5). Constraint (6) relates that at most, one worker can be allocated to a machine during the day. It also allows the idleness of a machine. Constraint (7) expresses that jobs can be assigned to a machine only if a worker is allocated to the respected machine. The completion time of a job is calculated by Constraints (8)- (10) according to the determined machine-configuration combination, the previous configuration on the machine, the related allocated worker, his/her arrival time, and sanitization time. Constraint (11) specifies the makespan, while Constraint (12) restricts it to the maximum allowed end time of the day. The working time of a worker is determined by Constraint (13), while the maximum allowed working time of a worker is enforced by Constraint (14) . Constraint (15) restricts the number of workers that can arrive at each arrival time. Finally, Constraint (16) defines the domains of the decision variables.

CP has been recognized as a powerful solution approach to tackle combinatorial optimization problems. Although it was initially developed to find high-quality feasible solutions within a reasonable time, the current commercial CP solvers can find optimal solutions and prove the optimality [21] . This approach has successfully been applied in various optimization problems, more noticeably in production scheduling and sequencing (e.g., Safarzadeh and Kianfar [44] , Mokhtarzadeh et al. [45] , Ham [46] ). In this paper, IBM ILOG CP Optimizer is used to model and solve the problem.

Two fundamental variables in the CP Optimizer are interval and sequence variables, which are very useful to formulate scheduling problems. An interval represents an activity, task, or operation possessing some characteristics, including startOf, endOf, sizeOf, lengthOf, and presenceOf. The sizeOf an interval is the actual processing time, while the lengthOf indicates the duration between the startOf and endOf an interval. An interval can be optional, which can represent the possible allocations of activity to resources. When an optional interval is present, its presenceOf is equal to true (one). A sequence variable indicates the permutation for a set of preset intervals. In the CP Optimizer, several specific constraints exist, which will be described J o u r n a l P r e -p r o o f later in the CP model. The sets, indices, and parameters of the developed CP model are the same as the described MILP model. The decision variables of the CP model are as follows.

Interval variable representing job Optional interval variable representing the assignment of job to machine with configuration Optional interval variable representing worker Optional interval variable representing the allocation of worker to machine who enters into the system at arrival time (

): Sequence variable representing the sequence of jobs assigned to machine . Moreover, the configurations of jobs processed by machine in each position of the sequence are tracked by . This will be useful for imposing reconfiguration time between every two consecutive jobs on a machine.

Subject to:

J o u r n a l P r e -p r o o f

Equation (17) indicates the objective function of the problem, which is similar to Equation (1); it consists of the three normalized terms, including the COVID-19 risk, workers' preferences, and makespan, respectively. Through alternative function, Constraint (18) determines the assignment of machine-configuration combination for each job respecting the eligibility restriction, and the start time and end time of the selected interval ( ) are equal to its associated job interval ( ). Likewise, in Constraint (19), a worker's allocated machine and arrival time are determined (if he/she is present) respecting the capability restriction. Constraint (20) expresses that at most, one worker is allocated to a machine. Constraint (21) ensures that a job can only be processed by an eligible machine-configuration combination if a capable worker is allocated to it. Constraints (22)-(24) are equivalent to Constraints (8)-(10) that determine the start time, size, and sequence of jobs on machines. The function in Constraint (24) guarantees that the intervals (jobs) assigned to a machine do not overlap, and the corresponding reconfiguration time between every two consecutive intervals is imposed (if a reconfiguration is needed). Constraint (25) calculates the makespan while Constraint (26) restricts it to the maximum allowed end time of the day. Constraints (27) and (28) specify the start time and end time of a worker, while the maximum allowed working duration of a worker is respected by Constraint (29) . Finally, Constraint (30) limits the maximum number of workers who can arrive at each arrival time.

In this section, first, two properties of an optimal solution for the investigated problem are presented. Then, the improved MILP and CP models are developed accordingly. Finally, a hybrid model is devised based on these improved MILP and CP models.

Proof. To prove this property, two cases should be considered:

Machines with completion times equal to : proceeding with a proof by contradiction, let indicate an optimal schedule on a machine with the completion time of , in which Property 1 is not held. As a result, at least one job is performed separately from other jobs with the same configuration on the corresponding machine. Fig. 2(a) 

Since there is no precedence between jobs, inserting , -right after , -, a feasible solution, namely , which is shown in Fig. 2(b) , can be achieved. Accordingly, the makespan of the new solution is calculated by: 

Thus,

. In addition, since this movement has no impact on the allocation of workers to machines and their arrival times, the other two terms of the objective function (the COVID-19 risk, and workers' preferences) are unchanged. As a result, is better than , which is in contradiction to the assumption. It is worth mentioning that if more than one job does not comply with Property 1, the difference between the makespan of and can be even widened.

Machines with completion times less than : suppose that Fig. 2(a) represents an optimal schedule on a machine with a completion time less than . By performing , -right after , -, the completion time of the machine is improved. However, since it does not worsen , and the other terms of the objective function are unchanged, Property 1 is held, and it completes the proof.  

Proof. Following Property 1, since there is no reconfiguration time between the jobs with the same configuration on a machine, changing the order of jobs within a block does not change the completion time of the respected machine. Also, the two other terms of the objective function (The COVID-19 risk, and workers' preferences) remain unchanged. Since instead of sequencing jobs on machines, just sequencing of blocks is required, it can be helpful to reduce the number of relevant binary variables in the mathematical model. Fig. 2(c) presents a schedule ( ) where job , -in the optimal solution related to ( Fig. 2(b) ) is swapped with , -. As can be seen, the makespan does not change. 

Following Properties 1 and 2, besides allocating workers to machine-arrival time and assigning jobs to machine-configuration combination, just sequencing the blocks on the relevant machines is required. This considerably improves the model by reducing the number of related constraints and variables. In this regard, the required decision variables and the Improved MILP (I-MILP) model are provided as follows.

1 if block on machine is formed; 0, otherwise 1 if on machine , block is scheduled before block ; 0, otherwise Completion time of block on machine

Equation (1) Subject to: Constraints (2)-(6), (12) , (14) , and (15)

J o u r n a l P r e -p r o o f

Constraints (33) and (34) ensure that a job can be assigned to a machine-configuration combination if a worker is allocated to the machine and the relevant block is formed. Similar to Constraints (8)-(10), Constraints (35) and (36) calculate the completion time of each block on each machine based on its allocated worker, assigned jobs, and previous configuration. The makespan is determined by Constraint (37) . Constraint (38) calculates the working time of each worker. Lastly, the domains of the variables are defined in Constraint (39).

Similarly, new variables are defined for Improve CP (I-CP), and the modified formulation is presented as follows.

: Optional interval variable indicating configuration is utilized on machine . ( ): Sequence variable representing the sequence of jobs assigned to machine with configuration .

( ): Sequence variable representing the sequence of blocks utilized in machine . Also, the configuration in each position of the sequence is tracked by to impose reconfiguration time between every two consecutive configurations on a machine.

Equation (17) Subject to:

Constraint (40) guarantees that there is no overlap among the jobs performed on a machine with the same configuration. Constraint (41) expresses that if is present, it spans all jobs with configuration on machine so that it starts with the first job and finishes together with the last job. Constraint (42) indicates that on each machine, there is no overlap among and makes sure the reconfiguration times between the blocks are respected.

Providing an initial solution (starting point) to a CP model as an input may considerably improve its performance to obtain a better solution in a shorter computational time. The starting point is not necessary to be a complete or feasible solution. The CP Optimizer first visits this solution. If it is infeasible or incomplete, the information included in the starting point is utilized as a guideline for the search [21] . This technique, known as warm starting, has been successfully implemented in several studies (e.g., Ham A warm start solution can be generated by a heuristic method, an MILP, or even a CP model. To take advantage of the strength of both MILP (in the assignment) and CP (in sequencing) models, in this paper, a combination of the I-MILP and I-CP models is considered. Hereafter, this hybrid model is called I-CP with warm starting (I-CP-WS). The pseudo-code of the I-CP-WS is presented in Algorithm 1. First, the I-MILP solves the problem within a specified time limit (line 1). Then, the obtained assignments of jobs to machine-configuration combination, and allocations of workers to machine-arrival times (i.e., the values of and ) are used to set the present values of the related interval variables in the I-CP (i.e., and ) (Lines 2 to 10). Finally, the CP Optimizer employs this partial solution as a starting point for the I-CP model and then proceeds its procedure to solve the problem within the time limit (Lines 11 and 12). It should be noted that according to our preliminary experiments, ignoring the third term of the objective function (related to ) in the I-MILP leads to more satisfactory results in the I-CP-WS. 

This section provides comprehensive computational experiments to assess the performance of the developed MILP models, CP models, and the hybrid model as well. Furthermore, some relevant managerial insights are extracted through a sensitivity analysis conducted on a few influential parameters of the investigated problem. All the models are implemented in IBM ILOG CPLEX Optimization Studio Version: 20.1.0.0, and executed on a personal computer with Processor Intel (R) Core (TM) i5-7300HQ 2.5GHz and 12 GB of RAM.

Several instances are generated to conduct the computational experiments. Different numbers of jobs (15, 20, 30, 40, 50, 60, 70, 80, 90 , 100, and 110), machines/workers (3, 5, 6, 8, 9, 12, 15, and 16) , and configurations per machine (3, 5, and 7) are considered. A three-field name is employed to represent each instance. For example, 30j-5m-3c indicates an instance with 30 jobs, 5 machines/workers, and 3 configurations per machine. The instances are categorized into three groups based on their sizes. If the value of is less than 1200, the instance is considered small-size; medium-size for the value between 1200 and 4000; and large-size for the value greater than 4000. The following values are utilized to randomly generate the instances: ; ;

, which equals to 6:00 PM; ;

. Moreover, 20% of machine-configuration combinations are regarded eligible to process a job, and ⌊ ⌋ workers are considered being capable of working on a machine.

To examine the performance of the developed models, 33 instances are solved. Based on the size of an instance, the computational time limits of 100s, 300s, and 900s are applied for small-, medium-, and large-size instances, respectively. The Relative Percentage Deviation (RPD) is employed as an evaluation metric according to Equation (43) , in which is the objective found by a model, and is the best objective found among all models.

( ) Table 1 reports the computational time limits, the best objectives found among all models, the RPD values of each model, and the times when each model reached its best solutions. Based on the results, it can be observed that the MILP cannot find any feasible solution within the time limits for a large portion of the instances, especially for the large-size instances due to the memory outage (indicated with NS in Table 1 ). In contrast, implementing the proposed properties in I-MILP leads to a considerable improvement as I-MILP can find proper solutions for all the instances within reasonable computational times. The average RPD for the instances solved by the MILP is 62.51, while the I-MILP reaches the average RPD of 0.59 for these instances. Likewise, comparing the CP and I-CP, incorporating these properties results in about a 21% reduction in the average RPD (from 41.56 to 32.84). Moreover, applying the warm start technique in the I-CP-WS bring about the best results in the majority of instances. Fig. 3 depicts the computational times the I-MILP and I-CP-WS models reach their best solutions. Accordingly, not only the I-CP-WS obtained the best average RPD (0.91) compared to the I-MILP (4.44), but also it managed to reach its solutions in notably shorter computational times. As a result, the I-CP-WS model demonstrated itself as an efficient solution method, especially when a high-quality solution is required within a limited time.

To statistically evaluate the performance of the developed models, the individual value plot, and Games-Howell simultaneous tests for differences of means with 95% Confidence Intervals (CIs) are employed (using Minitab 18) based on the values of the RPD reported in Table 1 . Fig. 4 indicates the distribution of the obtained RPD for each model. As it is evident, the I-CP-WS is the most reliable model with the lowest variation and average RPDs followed by I-MILP. Furthermore, each pair of the models is compared in Fig. 5 . According to the Games-Howell simultaneous test, if an interval does not contain zero, there is a significant difference between the relevant pair of methods, and the method with a better average is superior. Accordingly, the I-CP-WS significantly outperformed the other models. Moreover, the superiority of the I-MILP over the CP and I-CP is evident while there is no statistical difference between the CP and I-CP. 

This subsection focuses on analyzing the impact of three influential and controllable factors on the objectives of the problem, including i) worker flexibility to work with different machines, ii) machine flexibility to perform various jobs based on different configurations, and iii) reconfiguration times. In this regard, instance 40j-6m-3c is considered to conduct the analysis.

To obtain more solid results, 10 different sets of parameters are randomly generated for this problem (based on Section 6.1), and in each case, the average changes of the 10 problems are reported. The results are provided in the following subsections, which can provide helpful insights for managers to shape their decisions to be prepared during the COVID-19 pandemic or even similar situations. Fig. 6 illustrates the percentages of changes in the three objectives and the integrated objective based on the number of machines each worker is capable of working with. As depicted, by increasing the number from 1 to 2 (3), the integrated objective (Obj) and the COVID-19 risk are significantly dropped by 22% (27%) and 25% (30%), respectively. However, the further increase does not noticeably impact them. Thus, worker flexibility can be a helpful remedy in such situations. As a result, managers should pay more attention to this fact by implementing some training strategies, like cross-training, to enhance the skillset and flexibility of their workforce up to a reasonable level, respecting their budget limits.

According to Fig. 7 , the positive impact of machine flexibility is evident. Respectively, by increasing the number of eligible machine-configuration combinations for each job from 2 to 4, the integrated objective, COVID-19 risk, and workers' preferences will be improved by 26%, 48.5%, and 26%. In this regard, this aspect can be appreciated by a reasonable investment in RMTs with higher levels of flexibility and providing more multi-function modules and tools for them. Undoubtedly, making a trade-off between budget and machine flexibility is crucial.

Despite the substantial advantages of RMSs, the drawbacks regarding reconfiguration, like reconfiguration times, can be noticeable. For this reason, a sensitivity analysis based on different ranges of reconfiguration times is conducted, and the results are demonstrated in Fig. 8 . Accordingly, the negative impact of an increase in reconfiguration times is evident. By doubling the reconfiguration times, the integrated objective, COVID-19 risk, and workers' preferences will be worsened by19.5%, 59.5%, and 22.5%, respectively. For this reason, the proper design or purchasing of modules and tools for RMTs, as well as adequate training of workers to quickly perform reconfigurations are of great importance. 

This paper delved into a novel integrated workforce planning and production scheduling problem, considering issues arising from the COVID-19 pandemic. In this regard, due to the promising capability of RMSs to adapt to sudden changes, a shop floor with multiple nonidentical parallel RMTs was studied. Accordingly, MILP and CP models were developed to formulate the problem while taking workers' COVID-19 health risk and workers' preferences for flexible working hours into account. Then, two properties of an optimal solution were excavated and incorporated in the developed models, namely I-MILP and I-CP models. Next, to enjoy the advantages of both models, a hybrid model (I-CP-WS) was introduced by warm starting the I-CP model with a partial solution provided by the I-MILP model. According to the computational experiments, the I-CP-WS model significantly outperformed the other models regarding both RPD and computational time criteria. Moreover, regarding the RPD, the superiority of the I-MILP model over the MILP, CP, and I-CP models is evident. Based on the sensitivity analysis, the worker and machine flexibility, and reconfiguration times have sizeable impacts on the objective of the problem, and particularly on the COVID-19 risk. The following recommendations can be applied to strengthen the health risk mitigation and workers' preferences against situations similar to the COVID-19 pandemic: a) implementing worker training strategies, like cross-training, to enhance the flexibility of workers, b) investing in RMTs with proper levels of flexibility, and modules and tools with wider functionality, respecting the budget limit, and c) careful design and selection of RMTs' modules as well as adequate training of workers to reduce reconfiguration times. Although the proposed models were formulated based on the specific conditions during the COVID-19 pandemic, they can be adjusted to other similar situations (e.g., an outbreak of a disease), in which a health risk can influence the allocation of workers. Applying collaborative robots (cobots) can also be fruitful to mitigate the related health risks. As a result, extending the problem to incorporate cobots and providing relevant models can be considered for future research. Moreover, exploring the problem in other manufacturing settings like cellular manufacturing and assembly lines, where the density of workers is high, is valuable. Furthermore, expanding the planning horizon to multi-day planning can lead to more practical results. Therefore, it can be regarded as a promising research direction.

The authors, whose names are listed immediately below, certify that they have no affiliations with or involvement in any organization or entity with any financial interest (such as honoraria; educational grants; participation in speakers' bureaus; membership, employment, consultancies, stock ownership, or other equity interest; and expert testimony or patent-licensing arrangements), or non-financial interest (such as personal or professional relationships, affiliations, knowledge or beliefs) in the subject matter or materials discussed in this manuscript.

 Studying an integrated workforce planning and production scheduling in an RMS facing COVID-19 crisis.  Devising novel MILP and CP models to formulate the problem.  Employing two properties of an optimal solution of the problem to enhance the proposed models.  Developing a hybrid MILP-CP model to take advantage of both models' strengths.  Conducting computational experiments to evaluate the performance of the proposed solution methods.

OR-methods for coping with the ripple effect in supply chains during COVID-19 pandemic: Managerial insights and research implications

Past, present, and future barriers to digital transformation in manufacturing: A review

United Nations Industrial Development Organization

Occupational Safety and Health Administration

Global Preparedness Monitoring Board. A world at risk: annual report on global preparedness for health emergencies

Design of reconfigurable manufacturing systems

Manufacturing system architecture for cost-effective mass-individualization. Manufacturing letters

Reconfigurable manufacturing systems: Literature review and research trend

Reconfigurable manufacturing systems from an optimisation perspective: a focused review of literature

Flexible and reconfigurable manufacturing systems paradigms. International journal of flexible manufacturing systems

A novel methodology to measure the responsiveness of RMTs in reconfigurable manufacturing system

Synchronization-oriented reconfiguration of FPAI under graduation intelligent manufacturing system in the COVID-19 pandemic and beyond

Reconfiguring and ramping-up ventilator production in the face of COVID-19: Can robots help?

Susceptibility to coronavirus (COVID-19) in occupational settings: The complex interplay between individual and workplace factors

Coronavirus disease (COVID-19): Health and safety in the workplace

Protecting Workers: Guidance on Mitigating and Preventing the Spread of COVID-19 in the Workplace

The Association of Local Authority Medical Advisors. Covid-19 Medical Risk Assessment

European Centre for Disease Prevention and Control. Heating, ventilation and air-conditioning systems in the context of COVID-19

Ventilation and air conditioning during the coronavirus (COVID-19) pandemic

Centers for Disease Control and Prevention. guidance for businesses and employers to plan and respond to coronavirus disease 2019

IBM ILOG CP optimizer for scheduling

A resource type analysis of the integrated project scheduling and personnel staffing problem

A decomposed branch-and-price procedure for integrating demand planning in personnel staffing problems

Minimising total tardiness for the worker assignment scheduling problem in identical parallelmachine models

Minimizing total tardiness for the machine scheduling and worker assignment problems in identical parallel machines using genetic algorithms

Uniform parallel-machine scheduling to minimize makespan with position-based learning curves

Scheduling identical parallel machines and operators within a period based changing mode

Scheduling of unrelated parallel manufacturing cells with limited human resources

Uniform Parallel Machine Scheduling with Dedicated Machines

GRASP algorithm for the unrelated parallel machine scheduling problem with setup times and additional resources

Bi-objective parallel machine scheduling with additional resources during setups

A combinatorial Benders decomposition algorithm for parallel machine scheduling with working-time restrictions

A combinatorial evolutionary algorithm for unrelated parallel machine scheduling problem with sequence and machine-dependent setup times, limited worker resources and learning effect

Solution approaches for a real-life resource-constrained parallel machine scheduling problem

Design of an efficient genetic algorithm for resource-constrained unrelated parallel machine scheduling problem with machine eligibility restrictions

Iterative job splitting algorithms for parallel machine scheduling with job splitting and setup resource constraints

A heuristic algorithm for identical parallel machine scheduling: splitting jobs, sequence-dependent setup times, and limited setup operators

Reconfigurable manufacturing systems: Principles, design, and future trends

A new heuristic for integrated process planning and scheduling in reconfigurable manufacturing systems

A multi-objective particle swarm optimisation for integrated configuration design and scheduling in reconfigurable manufacturing system

Flexible job shop scheduling problem with reconfigurable machine tools: An improved differential evolution algorithm

Distributed job-shop rescheduling problem considering reconfigurability of machines: a self-adaptive hybrid equilibrium optimiser

An improved genetic algorithm for flexible job shop scheduling problem considering reconfigurable machine tools with limited auxiliary modules

Job shop scheduling with the option of jobs outsourcing

Scheduling of human-robot collaboration in assembly of printed circuit boards: a constraint programming approach

Transfer-robot task scheduling in job shop

Acknowledgment. This work was supported by the European Regional Development Fund under the project AI&Reasoning (reg. no. CZ.02.1.01/0.0/0.0/15_003/0000466).

☒ The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.☐The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: