Job rotation and human–robot collaboration for enhancing ergonomics in assembly lines by a genetic algorithm

Currently, the largest percentage of the employed workforce in the manufacturing industry is involved in the assembly process, making ergonomics a key factor when dealing with assembly-related problems. During these processes, repetitive tasks and heavy component handling are frequent for workers, who may result overloaded from an energetic point of view, thus affecting several aspects not only relating to the human factor but also to potentially reduced productivity. Different organizational strategies and technological solutions could be adopted to overcome these drawbacks. For these purposes, the present paper proposes a genetic algorithm for solving the typical problem of assembly line balancing, taking into account job rotation and human–robot collaboration for enhancing ergonomics of workers. The objectives of the problem are related to both economic aspects and human factor: (i) the cost for implementing the assembly line is minimized, evaluated on the basis of the number of workers and differentiated by skill levels and on equipment installed on workstations, including collaborative robots, and (ii) the energy load variance among workers is also minimized, so as to smooth their energy expenditure in performing the assigned assembly operations, calculated according to their movements, physiological characteristics, job rotations and degree of collaboration with robots. The paper finally presents and discusses the application of the developed tool to an industrial assembly case.


Introduction
The human factor is becoming ever more central in industrial production systems. According to the 2030 Agenda for Sustainable Development [1], a specific goal is to promote safe and secure working environments for all workers, thus making ergonomics a key factor nowadays in the design of workplaces to ensure workers' well-being [2][3][4]. Motivations rely not only on safeguarding health and safety of personnel but also on economic aspects related to medical costs caused by physical and psychological damage and potential absenteeism, for which workers need to be replaced by unexperienced substitutes who may reduce productivity [5][6][7].
Among different processes, assembly currently involves the largest percentage of the employed workforce in the manufacturing industry, as most of the operations are manual. During assembly, repetitive movements and heavy components handling are frequent for workers, who may result physically stressed and overloaded from an energetic point of view [8,9].
Based on these arguments, it can be stated that the design of assembly lines should undergo ergonomic considerations, on which the process efficiency highly depends [8,10]. The simple assembly line balancing problem of type 1 (SALBP-1), whose main objective is to assign the operations necessary to assemble a mass product to the minimum number of workstations of a production line for a given cycle time, is one of the major problems to be solved in the field of assembly line design. Reasons rely on the mathematical combinatorial complexity of balancing, which is a NP-hard optimization problem, and on the great impact on economic factors its resolution has [11,12].
To foster an ergonomic planning of assembly lines, different organizational strategies and technological solutions 1 3 could be implemented, also taking into consideration personal features and physiological characteristics, such as gender and age, which differentiate each worker's physical capabilities [5].

State of the art and motivations of the work
In the scientific literature, different studies investigate ergonomics in assembly systems, mainly with the objective of reducing the ergonomic risk in terms of worker posture [13].
As an example, in [14], the design of assembly systems is aimed at improving ergonomics by avoiding work-related musculoskeletal disorders, potentially generated from certain tasks entailing frequent uncomfortable movements and positions. The outcome of the study highlights how productivity improvements can be reached with ergonomic improvements on assembly workstations. Other works develop digital tools for simulating the assembly process and for the design of system layouts to deal with ergonomics of workers [10,[15][16][17][18]. Among the solutions proposed for enhancing ergonomics, the implementation of job rotation strategies in assembly systems proved to effectively balance workload distribution among operators and reduce assembly errors [13,[19][20][21][22][23].
In the last years, several studies have been conducted to introduce the ergonomic risk also in the assembly line balancing problem [13]. In [24], different clusters of SALBP-1, named ErgoSALBP, are modelled to include the ergonomic risk into the assembly line balancing. The posture of the worker is evaluated in the work proposed in [25], through postural risk assessment techniques, with the objective of maximizing the physical workload balance and minimizing the accumulated risk of postures (ARPs), which evaluates the load on back, arms and legs. Studies presented in [26] propose a new objective function to minimize the maximum ergonomic risk on the line. The work developed in [27] deals with the optimization of traditional objectives for line balancing, i.e. the number of workstations and the workload variance, and of the ergonomic risk at each station, based on rapid upper limb assessment technique, which evaluates the strain of upper limbs, neck and torso. Also in other works [28,29], the ergonomic risk is minimized, together with assembly process times. The study conducted in [30] considers a line balancing problem with ergonomic constraints for the minimization of the number of workers. Mokhtarzadeh et al. [31] propose a two-stage framework to first classify the risk of each task using various standards, such as the National Institute for Occupational Safety and Health (NIOSH), Occupational Repetitive Actions (OCRA) and Ergonomic Assessment Worksheet (EAWS), to then smooth the tasks of each class assigned to each workstation. In [32], the problem objectives are balancing the cycle time and the physical workload of operations among stations, evaluated through rapid entire body assessment (REBA) method. Ergonomic risk factors are also considered in the balancing problem faced in other works employing OCRA assessment technique [33,34].
A different approach for considering ergonomics in the line balancing problem, mainly based on the assignment of assembly operations to workstations, is modelled in [35], where a time smoothness index and an energy smoothness index are optimized for SALBP-2 (in which the cycle time is minimized for a given number of workstations). In this case, ergonomic evaluations are based on energy expenditure of workers in carrying out assembly operations. The energy smoothness index is used to distribute the energetic workload among the different stations of the assembly line, disregarding however personal features of assigned workers, which can influence their physical capabilities. The ergonomic approach based on energy expenditure of workers has then been treated in other works concerning line balancing, such as [36].
In the academic research relating to assembly, different studies have also been carried out on human-robot collaboration (HRC) [37], the implementation of which could lead to advantages for workers' ergonomics. As a matter of fact, collaborative solutions for assembly lines could enhance not only process efficiency and flexibility, by combining human dexterity and cognitive capabilities with robot repetitiveness, but also reduce workers' physical effort, thanks to the robot support in heavy component handling [38][39][40][41][42]. Some assembly line balancing studies developed in the field propose the optimization of cost-related [43,44] and time-related aspects [45,46], while few others include ergonomic considerations [47,48]. Table 1 summarizes the contributions characterizing cited references, showing, among others, if the work deals with assembly line balancing (ALB) or with a different approach concerning assembly, such as analysis or system design, job scheduling or process planning.
As resulting from Table 1, in the scientific literature relating the complex formulation of assembly line balancing, job rotation and HRC are never treated simultaneously, although both currently contribute to a sustainable development of assembly processes. The two aspects are of fundamental importance to allow a balanced distribution of the workload on operators, aimed at safeguarding them, protecting their health, welfare and safety.
For these reasons, the present work aims at facing the previous gap by proposing a software tool based on a novel genetic algorithm (GA) to solve a multi-objective version of the SALBP-1 for collaborative assembly lines including job rotation. The objectives of the problem are related to both economic aspects and human factor. In particular, the cost for implementing the assembly line is minimized, evaluated on the basis of the number of workers and differentiated by 1 3 technical skill levels, and of equipment installed on workstations, including collaborative robots. The other objective to be minimized is the energy load variance among workers, so as to smooth their energy expenditure in performing the assigned assembly operations and not to exceed their physical capabilities. Energy expenditure parameter is calculated according to the movements in executing assembly tasks, to personal characteristics that differentiate each worker and to the degree of collaboration with robots.

Proposed approach
In the present study, assembly line workers are subject to a job rotation, an organizational strategy to periodically move workers among line stations, thus enabling different advantages for them, such as multi-tasking skills, reduction of monotony and improvement of ergonomics. Among various strategies analysed from both academic and industrial applications, a Ban rotation has been implemented, meaning that all workers can be assigned to all workstations [49]. To contrast the repetitiveness of movements, at each rotation, every worker is placed on a different workstation from the one of the previous rotation, given that a maximum of one worker can be assigned to each station. As regards the duration, the rotation of workers every 2 h on an 8-h shift is employed, as it is appointed as the most efficient strategy from various applications for assembly lines [50,51]. Furthermore, workers are characterized by different technical skills, which are intended in the present approach to range from a low to an intermediate until a high skill level. Workers can thus be assigned to workstations in which operations they are able to perform are allocated, so requiring a lower or equal skill level.
HRC has also been implemented for reducing workers physical effort by employing collaborative robots for burdensome operations. A workplace and time sharing system is considered, where human worker and collaborative robot are working in the same workplace and collaborating for assembly tasks [52]. Operations are characterized by the possibility to be performed manually and/or automatically and/or jointly by the human worker and the collaborative robot [53]. Given this, line stations can be (i) manual, if the assigned assembly tasks are performed by a worker; (ii) automatic, if tasks are performed by a robot; and (iii) collaborative, if tasks require workers and robots, working jointly or independently. Once assigned, robots are fixed and do not undergo any modification or rotation along the line. Moreover, automatic workstations, in which tasks are entirely executed by robots, are not considered for assigning workers during job rotation.
As an ergonomic factor, the present paper deals with the energy expenditure parameter, which is a measure of the worker's physical fatigue in performing manual operations. The energy expenditure evaluation is based on different aspects. In particular, it depends on personal characteristics and physiological parameters of workers, such as gender, age and weight, on body movements for completing manual tasks and on the eventual support of the robot in reducing the worker's physical load. The energy expenditure of workers is then compared to their energetic limits, again calculated according to personal characteristics, so as to ensure that no worker is physically overloaded or operating beyond capabilities.
Lastly, among the various resolution techniques for the problem proposed in the years, the software tool uses a GA to solve the SALBP-1, for the proven efficiency in reaching optimal solutions for complex and multiple constraint optimization problems [54]. Design ✓ 20 Scheduling Planning

Problem modelling
This section presents the problem formulation for the proposed multi-objective version of the SALBP-1. The notations considered in the mathematical model for indices are defined as follows: The parameters of the problem are: Parameters associated with each movement are present when applicable, based on the type of movement (e.g. lifting or lowering an object, walking).
As far as assembly operations are concerned, precedence constraints are contained in a binary square matrix, in which Equipment include tools, automatic screwdrivers, pressing machines, vibratory feeders and collaborative robots necessary to accomplish assembly operations. These data are coded by an equipment matrix, where each element EQ i,q can be As regards workers, gender is binary coded; in particular Technical skill level S z are coded as follows: The energetic limit, evaluated to verify that no worker is operating beyond physical capabilities, is calculated as a function of the maximum power of a worker maxW z (W), according to the calculations for professional applications proposed in [55]. Specifically, the energetic limit EL z (kcal) of worker z is evaluated for a work shift of 8 h proportional to the 35% of maxW z : where T (min) is the cycle time and 0.06 is a conversion factor of units of measure. maxW z depends on maximum oxygen consumption O2 z (l/min), one of the most used physical condition index. The energy equivalent of oxygen at maximum exercise is 21.25 kJ/l. Thus, properly converting this value into consistent units, maxW z , can be calculated as [55]: The maximum oxygen consumption of workers depends on their peculiar characteristics and is calculated as proposed in [56,57]: Variables of the problem, whose value depends on the assignment resulting by the solution, are evaluated as follows.
Execution times can assume from 1 to 3 values, according to the possibility of performing movements and operations manually and/or automatically and/or collaboratively. For manual execution, time T i (min) of operation i is given by the sum of times t k (min) of the various movements k that the worker performs during operation i. For automatic execution, T i is equal to the time for the robot to perform operation i. For collaborative execution, T i (min) is given by the sum of times t k (min) of the various manual 1 if worker z has a low skill level; 2 if worker z has an intermediate skill level; 3 if worker z has a high skill level.
if worker z is a woman movements k that the worker performs during operation i, plus the time for performing collaborative movements. Skill level can assume different values as well, varying based on the manual and/or collaborative execution. In particular, SL i identifies the level of technical ability the worker must have at least to carry out operation i, codified as follows: The energy expenditure of workers is evaluated according to the calculations proposed by Garg [58]. The Garg model relates energy expenditure to the operator's oxygen consumption in executing the various movements that compose each manual operation. Once energy expenditure of movements is calculated for each assembly operation, the total energy expenditure E z (kcal) of worker z can be evaluated. E z (kcal) is given as the sum of the energy expenditure E z , w (kcal) of the various movements manually performed that compose assembly operations allocated to workstation w where worker z is placed, plus the energy at rest for maintaining the position till the end of cycle time: with where e k,i,z is the energy expenditure for movement k of operation i allocated to workstation w where worker z is placed.
Ė 0 z (kcal/min) is the energy expenditure rate for maintaining the position over time, calculated as a function of the worker's weight W z (kg) and of the working position, as viewable in [58], while IT w is the idle time of workstation w (min).
As regards the objective function F, the SALBP-1 version the present work deals with is to assign assembly operations to workstations in order to minimize the cost for implementing the assembly line and the energy load variance among workers. To this purpose, the following expression has been formulated: Cand̃ are the values of the cost of the line C and of the energy load variance σ properly normalized between their lower and upper bounds, evaluated through simulations and empirical considerations; and (1-α) are the relative weight of objectives; p is a penalty value related to the satisfaction of problem constraints, equal to The first objective to be minimized is the cost C (€/h) for implementing the assembly line, given by being nw is the number of workstations established by allocating the operations into workstations such that the sum of execution times does not exceed the cycle time; C z,w (€/h) is the hourly cost of worker z placed in workstation w that range from 0 if no worker is present and then growing with technical skill level growing; nq is the number of equipment units allocated to workstation w; and C q,w (€/h) is the hourly cost of the equipment q at workstation w. This formulation of the cost objective allows not only avoiding duplication of resources, but also minimizing the number of workstations, consistently with SALBP-1 formulation.
The other objective to be minimized is the energy load variance σ among workers, given by where nz is the number of workers and ΔE is the mean energy gap among workers, being ΔE z the normalized value of energy gap between the physical limit of worker z and the energy expenditure for performing the assembly operations. Given the implemented job rotation, each worker in an 8-h shift is assigned to 4 workstations (one every 2 h) where different assembly operations are allocated, thus requiring different energy expenditure at each rotation. Therefore, ΔE z is obtained as follows: with ΔE z,r is normalized value of energy gap between the physical limit of worker z and the energy expenditure during rotation r, calculated for each rotation as The problem is subject to various constraints that, besides traditional ones on cycle time and precedence relationships to be respected, have been specifically introduced to deal with this specific version of SALBP-1. In particular, for each worker z at each job rotation r: 0.50 if more than 5 constraints are not respected 0.65 if 3 or 4 constraints are not respected 0.80 if 1 or 2 constraints are not respected; 1 otherwise • The energy expenditure must not exceed the energetic limit: • The maximum power maxW z , must not be exceeded: being W i,z (W) the power required to perform operation i for worker z. • The technical skill level S z , must be equal or higher than the maximum one maxSL w required to carry out assembly operations at workstation w where worker z is placed: • The workstation w z,r where worker z is placed during rotation r must be different from the one occupied at the previous rotation w z,r-1 :

Genetic algorithm architecture
The software tool based on the use of a GA for the resolution of the modelled problem has been developed on MATLAB®. The chromosome structure for the proposed GA is based on a tailored encoding method for the particular problem that uses the task-oriented representation for the assembly sequence. With each task-oriented chromosome, generated by a sequence planner that ensures the respect of precedence constraints, two other elements are randomly created and associated: a string named equipment-oriented chromosome with a number of elements equal to the number of workstations (nw) containing the necessary equipment for each station, including collaborative robots, and a matrix named job rotation-oriented chromosome of dimensions 4xnw with information on job rotation.
An example of chromosome, structured as proposed, is shown in Fig. 1. The task-oriented chromosome is composed of 10 genes encoding 10 assembly operations placed in order of execution. The same figure also shows the equipment-oriented chromosome with equipment (S for automatic screwdriver, P for pressing machine, V1 for vibratory feeder for screws, V2 for vibratory feeder for nuts) and robots (R) assigned to the 5 workstations of the line (nw = 5), distinct by different colours. In the underlying job rotation-oriented chromosome, for each of the 4 rotations, the sequence of workers, among the ones available in a company, is reported. As examples, in the first workstation, in yellow, assembly operations 3 and 2 are allocated and need to be manually performed with the use of a screwdriver (S) and vibratory feeder for screws (V1) by worker A at the first rotation, C at the second, D at the third and so on. In the second workstation, only operation 4 is automatically performed. Conversely, the third workstation is collaborative, and robot (R) and worker (first C, then A and so on) working together helped by pressing machine (P) and vibratory feeder for nuts (V2).
To evaluate each chromosome, the objective function F described by Eq. 11 is used as fitness function. The lower the value of F is, the higher the probability is for the chromosome to survive to following generations of the GA.
Roulette wheel selection is applied at each iteration of the GA to select chromosomes for producing offspring, according to their fitness value. The order-based crossover is then employed to ensure feasible assembly sequences in terms of precedence constraints.
As shown in Fig. 2, this genetic operator is applied to the task-oriented chromosome by selecting a couple of parents, each of them randomly divided into initial, intermediate and final part; the two parents produce two children, each of them made by initial and final part of one parent and by missing genes ordered as they are in the other parent. Then, the swap mutation is used for exchanging two randomly selected genes in task-oriented chromosome. This genetic operator is exemplified in Fig. 3 and has been implemented to enable proper randomness in the algorithm, which could lead to improved solutions through generations. Equipment and job-oriented parts of the chromosome remain unchanged after genetic operations, so that the fitness function is eventually penalized if some constraints are violated. Strong elitism is then applied to only preserve chromosomes with outperforming fitness values, for enforcing a steady enhancement of solutions.

Validation
Different numerical examples and real case studies have been analysed for testing the software tool. In this paper, results for the assembly process of a product taken from the automotive industry are presented. In particular, the assembly of a vehicle front-end has been selected because they Task-oriented   chromosome  3  2  4  5  1  6  9  7  8 10 Fig. 1 Example of chromosome structure, divided in task-oriented chromosome with assembly sequence; equipment-oriented chromosome with equipment and robots in workstations; job rotationoriented chromosome with workers rotating among manual and collaborative stations are representative in terms of both market diffusions, so as to justify the investment for a collaborative production line, and necessary equipment, especially robots for supporting workers in handling large components.

Dataset
The exploded view of the vehicle front-end is shown in Fig. 4. The assembly sequence for the 36 depicted components is formed by 29 operations. The assumptions for the problem include a production rate of 12 products/hour, a line efficiency of 96% and thus a cycle time of 4.80 min/product. The experiments have been conducted with a population size of 100 and a number of iterations of 150. Crossover and mutation probability have been set to 0.98 and 0.1, respectively, while 0.5 has been used as value (Eq. 11) to prioritize both objectives during the optimization through the same weight.
Initial data for the problem are given in Table 2 as regards labour and equipment costs while in Table 3 as regards the assembly process. In particular, for each task, the following input are reported: precedence constraints; manual execution time; automatic execution time; collaborative execution time; required skill level for manual execution; required skill level for collaborative execution; and needed equipment (robots are optional). To calculate energy expenditure, movements have been identified for every assembly operation, reaching a total of 177 movements for the 29 operations. As an example, movements have been detailed for the manual execution of task n.10 in Table 4, where related energy expenditures are also calculated for two workers (I and R) selected from the group of available operators simulated for the validation, whose peculiar characteristics are shown in Table 5.
As stated before, energy expenditures are calculated according to the formulas reported by Garg in [58]. To provide a better understanding, let consider, as examples, movements 1 and 3 in Table 4, walking and lifting, respectively. The values of e k,,i,z reported in the table are calculated as follows [58]: For walking: For lifting:

Results and discussion
The optimal solution obtained with the developed algorithm using the dataset depicted in Sect. 4.1 is represented in Fig. 5, where the relative chromosome structure is shown. The solution resulted to be the one with the best fitness e k,,i,z = 10 −2 (51 + 2.54W z sw 2 Fig. 2 Example of order-based crossover applied to task-oriented chromosome Fig. 3 Example of swap mutation applied to task-oriented chromosome of the first child of Fig. 2 value, equal to 0.83, very close to the maximum value of 1, meaning that both objectives of the problem have been achieved. As viewable from Fig. 5, the chromosome returned by the proposed software tool provides an assembly line made of 6 workstations, differentiated by colours, and requiring 6 workers, selected from the group of Table 5. The computed solution is also represented in Fig. 6 in terms of the energy load histogram for the 6 workers. In this graph, the height of bars corresponds to ΔE z , the mean energy expenditure of workers in carrying out the assigned operations during the 4 rotations, explained by Eq. 15. The same figure shows, below the graph, the values of the energy gap between the physical limit and the energy expenditure of every worker for each of the 4 rotations ( ΔE z,r (Eq. 16)). Energy expenditure is reported in normalized value between 0 and 1, to be comparable among workers.
To better clarify, in the fifth workstation of the line (the green one in Fig. 5), operations 23, 29, 24 and 26, presented in Table 3, are performed through the use of a robot (R) and a screwdriver (S) by the following workers selected from the group of Table 5: W at the first rotation, O at the second, G at the third and U at the fourth. These workers are subject to different energy consumptions. As an example, the energy expenditure of worker U (the first column of Fig. 6) is 0.27  at the first rotation, 0.41 at the second, 0.27 at the third and 0.36 at the fourth, producing a mean value of 0.33 that corresponds to the height of the bar. As regards the first objective relating costs, the software shows to be able to assign to each workstation a worker having the exact skill level required by the allocated tasks, so as to minimize the number of high-skilled workers on the line and thus reduce costs. As an example, workers J an L are the only 2, among the 6 workers selected, to have a high skill level (S J = S L = 3) and are alternatively placed on workstations 2 and 3, as viewable from Fig. 5, where the operations requiring a high skill level are performed, such as 1, 2, 17, 18 and 9, 15 (SL 1 = SL 2 = SL 9 = SL 15 = S L 17 = SL 18 = 3). Similar considerations apply to equipment resources; the only operations requiring a pressing machine (i.e. 5, 6 and 10) are all executed in the first workstation, where a pressing machine will be placed; operations 3, 4, 9 and 15 to be performed with both a vibratory feeder for  Table 4 Example of movement decomposition referred to manual task n.10 of Table 3  Relating to the other objective of energy load variance, in the proposed solution, the last workstation is completely manual, while the others entail HRC, as viewable from Fig. 5. As a matter of fact, 47% of operations in the proposed configuration require the support of a robot, highlighting the importance of collaborative work in manual assembly processes. Indeed, the strong presence of robots is reflected in data related to workers' energy expenditure. As shown in Fig. 6, the physical limit, given by the maximum energy expenditure value ( ΔE = 1 ), is always respected; furthermore, the use of collaborative robots allows a reduction in energy expenditure, which is largely below half of the energetic limit for all workers ( ΔE z ≤ 0.4∀z) that are thus working beyond their physical capabilities. Also, Fig. 6 reports a well-balanced workload for operators, which are subject to similar energy expenditure during rotations, thus demonstrating to reduce the energy load variance.
The correct functioning of the software tool can also be deduced from Fig. 7, where a comparison on the energy load histogram of Fig. 6, obtained at the final iteration number 150 (in orange), with the one of solution found with the algorithm running, precisely at iteration 25 (in blue), is shown. Workers reported in Fig. 7 refer to the blue solution. The fitness value is in this case 0.48, nearly half of the fitness of the best solution. As shown, the energetic workload is not well distributed among workers, which are almost all subject to higher energy expenditure. These types of assignments are discarded during the algorithm execution to prioritize and evolve towards improved solutions, as the one given in Fig. 6.
To better appreciate these results, the basic scenario described in Sect. 4.1 has been also used for some comparisons within the same proposed approach, given the novelty of the objectives introduced in this work. In particular, the solution obtained in Fig. 5 and represented in Fig. 6 has been compared with two different configurations, obtained maintaining the remaining conditions and re-executing the GA: in the first, workers are not subject to a job rotation but remain fixed in the assigned workstations; in the second, collaborative robots are not employed, so the entire assembly process is carried out manually. The configuration without job rotation is shown in Fig. 8 and can be noticed to produce a worst distribution of the energy load, if compared to Fig. 6, where job rotation is implemented. In fact, the system is still able not to violate energetic limits of workers, but the absence of job rotation clearly worsens workload smoothness among them. In the case of employing job rotation but not HRC, results are shown in Fig. 9. As viewable, if collaborative robots are not present, one worker and one workstation more are necessary on the line because some assembly operations (e.g. 1,9,20) may require longer times if manually performed, compared to a collaborative or automatic execution

Fig. 5
Chromosome structure for the best solution returned by the algorithm for product of Fig. 2 with robots. These increased times may result in the cycle time to be exceeded in certain points of the assembly line, so as to justify the need for more workstations. Moreover, workers have a higher mean energy expenditure (Fig. 9) if compared with the line configuration involving HRC (Fig. 6). This is due to the support collaborative robots can give to relieve workers from an energetic overload. Although Fig. 8 and 9 represent the best solution in terms of fitness function for the two further tested configurations, the energy load histograms may appear not so clearly balanced, mainly due to the different constraints of the problem. Nonetheless, the proposed approach is more effective in the basic scenario and demonstrates, through the computed results, to simultaneously minimize costs and energy load variance, achieving optimal compromises on the proposed objectives, thanks to the adoption of job rotation and HRC.
Therefore, the computed results show how the software tool simultaneously minimizes costs and energy load variance, achieving optimal compromises on the proposed objectives, thanks to the adoption of job rotation and HRC.
In Table 6, the characteristics of other examples and case studies tested for the validation are presented, together with the best value of the fitness reached. The results achieved highlight the versatility of the system and show its capability to deal with various problems and address different issues.

Conclusions
This paper presents a genetic algorithm-based software tool aimed at solving a novel version of the SALBP-1. The system optimizes not only production parameters, but also   ergonomics of workers, still little discussed in the scientific literature in terms of balancing the assembly line taking into account distinguishing features of workers for determining their physical capabilities. Assembly workers overloaded from an energetic point of view could undergo physical damage, also leading to medical costs and decreased process productivity. The innovative system has thus been developed to reduce the costs of implementing the assembly line and the energy load variance among workers, through to the use of job rotation and collaborative robots.
The computed results demonstrate the capability of the software tool to simultaneously minimize both objectives. Cost-efficient configurations are generated thanks to an optimal allocation of operations to workstations, whose number is minimized, according to required resources, namely, skilled workers, collaborative robots and equipment. Assembly operations needing a particular technical skill are grouped in the same workstation where the algorithm assigns a worker, from the ones at disposal for a company, having the exact skill level needed. Robots and other equipment, such as machines and tools, are allocated in a few workstations as well, where tasks requiring a particular resource are performed. Grouping operations requiring the same skill level and the same equipment allows cost-effective configurations of the assembly line. An ergonomic distribution of operations is also achieved, as demonstrated by the obtained solutions. As found in the literature, the implementation of job rotation confirms to improve ergonomics of workers. For the proposed system, experimental results show a better and smoothed distribution of energy expenditure among workers, calculated according to their personal characteristics, compared to the case of absence of job rotation. Even collaborative robots prove to enhance ergonomics, by relieving humans from heavy components handling and burdensome activities. In all tested scenarios, the tool is able to assign operations to workers not exceeding their physical capabilities and simultaneously respecting the other constraints of such a combinatorial and complex problem. This confirms its conceptual efficacy and demonstrates the context-based capability of providing the optimal solution.
Future research could be conducted with regard to other aspects improving ergonomics for workers of assembly lines, such as the reduction of perceived noise during the process. In addition, further developments of the tool could relate the extension to multi-product assembly lines and the introduction of various types of collaborative robots, differentiated by the level of dexterity.
Acknowledgements The authors wish to thank the MSc student S. Corsinovi from the University of Pisa (Italy) for his support in the experimental tests.
Author contribution Both authors MDM, GD have been responsible for planning and coordinating the steps of the research. The corresponding author MDM has been responsible for writing and reviewing the paper. GD has been responsible for revising the entire article, including structure and contents.
Funding Open access funding provided by Università di Pisa within the CRUI-CARE Agreement.

Data availability Not applicable.
Code availability Not applicable.

Conflict of interest
The authors declare no competing interests.

Ethical approval
The research activity envisaged in this work has been conducted applying fundamental ethical principles. We confirm that all the authors involved in the writing of this article are aware of this work and approve all its contents. In addition, the paper has not been published or submitted to any other journal.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.