Decision-Making Framework for Implementing Safer Human-Robot Collaboration Workstations: System Dynamics Modeling

: Human-Robot Collaboration (HRC) systems are often implemented seeking for reducing risk of Work-related Musculoskeletal Disorders (WMSD) development and increasing productivity. The challenge is to successfully implement an industrial HRC to manage those factors, considering that non-linear behaviors of complex systems can produce counterintuitive effects. Therefore, the aim of this study was to design a decision-making framework considering the key ergonomic methods and using a computational model for simulations. It considered the main systemic inﬂuences when implementing a collaborative robot (cobot) into a production system and simulated scenarios of productivity and WMSD risk. In order to verify whether the computational model for simulating scenarios would be useful in the framework, a case study in a manual assembly workstation was conducted. The results show that both cycle time and WMSD risk depend on the Level of Collaboration (LoC). The proposed framework helps deciding which cobot to implement in a context of industrial assembly process. System dynamics were used to understand the actual behavior of all factors and to predict scenarios. Finally, the framework presented a clear roadmap for the future development of an industrial HRC system, drastically reducing risk management in decision-making.


Introduction
The recent advance in industrial technology is dealing with complex problems: to increase productivity without neglecting human factors; to look for specific improvements in spite of considering the whole system dynamics; and to solve today's problems without creating a new one to be addressed in the future [1,2]. The implementation of a Human-Robot Collaboration (HRC) system is often aimed at reducing Work-related Musculoskeletal Disorders (WMSD) at the same time that increases production, by combining the skills of both humans and robots [3,4]. Workers are exposed daily to physical and mental workplace hazards [5], especially in assembly workstations in industry [6]. Ergonomic risks lead to absenteeism due to WMSD [7,8]. In order to mitigate those issues, a possible ergonomic intervention is to implement an industrial HRC system, which combines both robot and human skills. Moreover, a cobot performs the task like a co-worker, with no physical barriers to help achieve a goal [9]. As a co-worker, it is also expected to increase productivity [10].
The implementation of an industrial HRC system is not simple and there are many organizational factors to consider [11,12]. Moreover, there are critical technical issues, ergonomic, safety, and economic aspects to consider when evaluating scenarios with a cobot [13,14]. By doing a correct assessment of the current workstation it is possible to run a simulation model [15], which brings more data for understanding the system and ultimately for a correct decision. Therefore, the aim of this study was to develop a decision-making framework for implementing HRC systems. The contribution highlights of this work follow: a framework to decide which collaborative robotic system to choose; a system dynamic model that considers Level of Collaboration (LoC), physical, and mental workloads with simulation of scenarios for absenteeism and productivity; and a case study applying a roadmap to reduce risk management in decision-making. The system dynamics model is limited to the boundaries chosen, which excludes other variables that might be useful depending on the characteristics screened for a specific workplace. The article is structured as follows: • An introduction to the fundamental concepts related to the problem being addressed, namely: HRC, ergonomics, and system dynamics. • A decision-making framework proposal, starting with the problem definition, followed by an assessment of the current workstation, a computational modeling of the system, the evaluation of possible solutions, and the final decision regarding an industrial HRC system. • A case study applying the aforementioned framework in a manual assembly workstation that intends to implement an industrial HRC.

Human-Robot Collaboration
The fourth industrial revolution is characterized by smart and autonomous systems [9]. Cobots allow workstations without barriers and their implementation is a potential solution when planning to improve workstation conditions. This solution combines typical robot's capabilities (accuracy, speed, and repeatability) with human skills (adaptability, dexterity, and perception) aimed at achieving a common goal [16]. According to [17], there are five LoC with increasing complexity, competences, and cost: Level 0 (cell), Level 1 (coexistence), Level 2 (synchrony), Level 3 (cooperation), and Level 4 (collaboration). The design of HRC workstations has to integrate some requirements in order to propose task allocations [18,19]: (Human) posture, strength, biomechanics and cognition; (Robot) joints position, velocity, acceleration and sensors; (Human and Robot interaction) force or feedback force; and (Process) manufacturing goals, such as cycle time and quality. Moreover, when designing an HRC system and considering ergonomic aspects, it is expected to reduce human workload, WMSD risk factors, and production costs. Therefore, many of the current workstations should consider the implementation of a cobot to work in a safe, ergonomic, and efficient way. Accordingly, a successful implementation of a HRC system depends on considering ergonomics [11].

The Role of Ergonomics in the Workstation Transformations
Ergonomics is the discipline that studies the relationship between man and work, seeking for adaptation of the workstation environment for the worker to carry out his activities [20]. New technologies and interventions in industry increases environment complexity, which means that ergonomics should be treated using a systematic view [21]. Organizations often do not get the most ergonomic benefits when they introduce new technologies [11]. It is confirmed by [22] that optimizing simultaneously task allocation while taking into account ergonomic aspects improves the efficiency and acceptance. Moreover, the concept design considers human mental and physical viewpoints at the same time by proposing a hand-guiding on the robot for the users to have control over the system, having the robot as a helper instead of a robot giving an object. A systemic view that predicts scenarios regarding ergonomics and productivity is needed, and not considering ergonomics in HRC systems may present undesired effects [23].

System Dynamics
System dynamics is a method to navigate the complex systems. It is grounded in the theory of nonlinear dynamics and feedback control developed in mathematics, physics, and engineering to understand a system's behavior [24]. Through this method, it is possible to design, to model, to simulate, to understand, and finally to make good decisions in complex systems taking both short-and long-term behavior into consideration. These systems exist in an environment within boundaries and are characterized by their structure, elements that interact, processes, goals, inputs, and outputs. As seen in [15], computer simulations are controlled experimental environments in which different possibilities can be tested without investing a great amount of resources before evaluating and deciding.
The interest in this topic has grown recently as stated in [25] that ergonomic study should always involve a systemic view, with any analysis or interventions considering the context as a whole. It was found in [26] that using system dynamics has potential to increase safety in complex systems that cannot be achieved with traditional approaches [27]. Similarly, ergonomics is also of crucial importance in the implementation of an industrial HRC. The three ergonomic dimensions have been considered in [28] when designing workstations for manual production and maintenance processes. In [29], employees suffering from injuries result in worse productivity and quality of products, which leads to losses and indirect costs. Indirect costs like absenteeism are often difficult to measure and are excluded from production costing models. A literature review found ergonomic frameworks with the following characteristics: focusing on the physical assessment [30]; an integrated framework addressing both physical and psychological levels [31]; an agent-based model included ergonomics in terms of assessing risk levels, considering how pain and fatigue affect performance [32]; an optimized hierarchical model to improve production time and ergonomics for decision-making [33]; and insights on the productivity and its ergonomic impact on workers [34]. It demonstrates that a framework is needed in order to avoid injury, poor quality, low production, and to decrease costs in complex systems.

Decision-Making Framework
A logical decision-making process was built analyzing the current situation, evaluating the available possibilities, and finally deciding what to do based on information about the system. This methodology was chosen to create a complete roadmap framework to be used together with system dynamics. It was divided in the following phases: (i) assessment of the current workstation by applying methods from the three ergonomic dimensions (physical, cognitive, and organizational) and collecting production data; (ii) computational modeling of the key factors of the system considering different LoC; (iii) simulation and prospection of scenarios regarding absenteeism and productivity when inserting a cobot into a production system; and (iv) technical, ergonomic, safety, and economic evaluation for a better decision-making.

Problem Definition
In the cases where this framework can be useful, usually the problem is related to both absenteeism due to WMSD and productivity. The reason is simple: by introducing an industrial HRC system, it is expected to reduce the human efforts regarding postures and repetitiveness. Additionally, by including automation into the production line, it is expected to increase productivity. Therefore, the framework is designed for those who need to choose what kind of cobot to install or what LoC fits better in the long term of this cooperation.

Ergonomic Assessment
Based on the characteristics of an HRC workstation, ergonomic methods from the three ergonomic domains have been used to assess the risks. The framework provides guidance, but does not suggest a strict set of tools, allowing more adaptability for the practitioner who must navigate within a particular context. Four physical risk levels and two mental Safety 2021, 7, 75 4 of 20 risk levels were defined by the simulation modelers in order to weight the overall human workload. It is suggested to use at least one from each of the following domains:

Organizational
Ergonomic methods to assess organizational factors are mandatory to allow building the causal loop diagram (CLD) [15], a diagram that represents the actual system under study. There are two main options to choose for a complete workstation description, according to preferences: • Ergonomic Analysis of Work [35]-an observational screening method to assess the workstation. It is divided into five stages: demand analysis, task analysis, activity analysis, diagnosis, and recommendations. By applying this method, it is possible to identify the main variables of a system and their interconnections in order to build the CLD of a dynamic system.

Physical
The physical workload is directly related to high force loads, awkward postures, and repetitive movements, which lead to WMSD and absenteeism [38]. To understand the system, this framework suggests choosing at least one of the ergonomic methods in this domain to assess the risks related to the workstation. The method should be chosen according to preferences and the specific characteristics of an assembly task: • Rapid Upper Limb Assessment (RULA) [39]-an observational ergonomic tool that considers biomechanical and postural load requirements of job tasks. It is a good and widely used method to assess physical workload, except for not considering the duration of exposures.  Each one of the physical workload methods are based on their own scale to assess the physical risk of a task. These scales are usually named as acceptable risk, low risk, Safety 2021, 7, 75 5 of 20 medium risk, and high risk. As this framework intends to use the risk level in a simulation model, the different output values from the selected physical methods were converted into a numerical scale of four physical risk levels. According to the aforementioned references, Table 1 presents the four suggested physical risk levels that can be obtained for each ergonomic method chosen.

Cognitive
Mental workload is often ignored when assessing the risks related with a workstation. However, many authors emphasize the importance in considering this factor in a system [5,47,48], especially in the long term. There are three main methods to choose in this domain: -a widely applied questionnaire used to assess mental workload, including work systems with a high level of complexity. It evaluates mental demand, physical demand, temporal demand, effort, frustration, and performance. A numerical scale is used to assess the workload perceived by the worker for each of the six items. • Subjective Workload Assessment Technique (SWAT) [50]-it was originally designed to assess aircraft cockpit workload. It is divided in two phases: scale development and scale scoring. The three dimensions measured are: time, mental effort, and psychological stress.

•
Electrodermal Activity (EDA) [51]-a technique to identify changes in the skin's electricity by wearable sensors. It may be employed for assessing emotional states and to understand the worker's mental status. EDA is divided in electrodermal response (EDR) that reflects short-term stress, and electrodermal level (EDL) is more related to risk perception and a relevant indicator of long-term stress.
Similarly, to the physical workload, as input to the simulation model, the different output values from the selected mental methods were converted into a numerical scale. According to the aforementioned references, Table 2 presents the two suggested mental risk levels that can be obtained for each method chosen.

Production data
The current data of the production line are necessary to model the dynamics of a system. In this regard, based on the literature [15,24], the most important factors suggested by this framework are:

•
Production goals: number of pieces to be assembled in a period of time.

•
Takt time: assembly duration time needed to match the production goals. • Cycle time: the time it takes to complete one assembly. • Absenteeism due to WMSD: the sick leave rate due to musculoskeletal issues.
• Number of workers: the sum of workers in the production line.
It means that a correct assessment on production is sufficient with the above factors. However, for understanding the dynamics of an actual system other factors may be useful as well (e.g., gender, age, seniority, previous injuries, or illnesses).

System Dynamics
This framework suggests to model the system under study using both ergonomic and production data. Modeling is divided into conceptual definition, structural description, and simulation. It aims at representing the most important factors and their interconnections with a focus on solving a problem. In this case, the context also includes the insertion of a cobot considering different LoC into the system to evaluate scenarios. As a result, it is expected to maximize worker's well-being and system performance. An example of conceptual description, the CLD with the above-mentioned factors is represented in Figure 1. In the conceptual definition, four information feedback cycles are relevant: "production self-regulation", "learning by repetition", "limitation due to disease incidence", and "reinforcement of mental workload".
The current data of the production line are necessary to model the dynamics of a system. In this regard, based on the literature [15,24], the most important factors suggested by this framework are:

•
Production goals: number of pieces to be assembled in a period of time.

•
Takt time: assembly duration time needed to match the production goals. • Cycle time: the time it takes to complete one assembly. • Absenteeism due to WMSD: the sick leave rate due to musculoskeletal issues. • Number of workers: the sum of workers in the production line.
It means that a correct assessment on production is sufficient with the above factors. However, for understanding the dynamics of an actual system other factors may be useful as well (e.g., gender, age, seniority, previous injuries, or illnesses).

System Dynamics
This framework suggests to model the system under study using both ergonomic and production data. Modeling is divided into conceptual definition, structural description, and simulation. It aims at representing the most important factors and their interconnections with a focus on solving a problem. In this case, the context also includes the insertion of a cobot considering different LoC into the system to evaluate scenarios. As a result, it is expected to maximize worker's well-being and system performance. An example of conceptual description, the CLD with the above-mentioned factors is represented in Figure 1. In the conceptual definition, four information feedback cycles are relevant: "production self-regulation", "learning by repetition", "limitation due to disease incidence", and "reinforcement of mental workload". As the problem is related to productivity and absenteeism due to WMSD, the hypothesis to be considered is at what extent the robot will collaborate. It is widely As the problem is related to productivity and absenteeism due to WMSD, the hypothesis to be considered is at what extent the robot will collaborate. It is widely demonstrated that the intervention of an industrial HRC system increases productivity and decreases ergonomic efforts [52][53][54]. Moreover, there are different competence levels for a cobot [55] to be chosen, which are very much related to the LoC. Which LoC is best for the current situation is what a simulation model can answer.
Therefore, by applying a structural and mathematical model, the Stock and Flow Map (SFM) is needed to understand what was established in the conceptual diagram, maintaining the cycles and the established concept itself. There are many examples of SFM in the literature that can be adapted to the current system; however, it must be tested and calibrated according to the situation. Table 3 presents the description of competences regarding safety and ergonomics, as well as the influence of LoC on the model.

Mental
Mental workload increases with the complexity of the task [47,52]. Physical Physical workload decreases when cobot assumes the tasks related with loads and repetitiveness [57,58]. Knowledge Knowledge of the task can assume different values depending on the specific workstation [15].
The system dynamics simulation phase is reserved to prospect scenarios for analysis of different possibilities. It aims at understanding the system to establish risk indicators, to propose organizational changes and ergonomic interventions, which finally support management evaluation of the system performance.

Management Evaluation and Decision
The evaluation process by the management is critical in deciding whether to include a robot in the system. After analyzing the prospected scenarios, there are other managerial aspects to consider. The minimal phases to be considered are shown in Figure 2.
Map (SFM) is needed to understand what was established in the conceptual diagram, maintaining the cycles and the established concept itself. There are many examples of SFM in the literature that can be adapted to the current system; however, it must be tested and calibrated according to the situation. Table 3 presents the description of competences regarding safety and ergonomics, as well as the influence of LoC on the model. Mental workload increases with the complexity of the task [47,52].
ysical Physical workload decreases when cobot assumes the tasks related with loads and repetitiveness [57,58]. wledge Knowledge of the task can assume different values depending on the specific workstation [15].
The system dynamics simulation phase is reserved to prospect scenarios for analysis of different possibilities. It aims at understanding the system to establish risk indicators, to propose organizational changes and ergonomic interventions, which finally support management evaluation of the system performance.

Management Evaluation and Decision
The evaluation process by the management is critical in deciding whether to include a robot in the system. After analyzing the prospected scenarios, there are other managerial aspects to consider. The minimal phases to be considered are shown in Figure 2.
Error! Reference source not found.  The technical evaluation considers the technologies available (robot arm, gripper, sensors), the new workstation configuration (shop floor space), computational systems to be integrated (interfaces, hardware, and software limitations), and even the experience and training of workers. Technical issues might be a critical variable when implementing a cobot. Ergonomic evaluations were already discussed during the first assessment.
Especially when deciding which tasks are to be allocated to the worker, the new workstation also needs special ergonomic considerations. Postures (workstation layout), frequency of repetitive actions, workloads, and time pressure as well as safety requirements. The new productivity goals are limited by ergonomic and safety issues. Therefore, when evaluating economics, every variable is important. The economic evaluation might be estimated considering the non-value added tasks that will be automated [14], as these tasks usually consume time and production resources. Therefore, it is possible to compare the savings with the investment. Another way to manage the economic evaluation is formulated in [55], where the cobot savings percentage of the manual takt time should be higher than the minimal required saving to justify implementing a cobot in the system.
There are other aspects that could be considered when evaluating the inclusion of a cobot for example, the product and process quality, which can be evaluated by the level of process variability and its standardization level [14]. A managerial decision is based on the available information. In order to avoid making decisions based on half information or mental models, this framework presented a logical roadmap to bring relevant data to be considered and help deciding at this last stage of the framework:

•
Summarizing the most important factors of the current situation; • Including the system dynamic model is key to understanding the whole system; • Considering the organizational recommendations from the simulation and prospection model; • Relying on the technical, ergonomic, safety, and economic evaluations; • Defining the new productivity and absenteeism goals.
The managerial team decides which cobot to implement. As the introduction of a new technology, it is important to have such a global approach. Figure 3 summarizes the proposed roadmap for a decision-making in HRC systems.

Case Study
A workstation for the manual assembly of electrical components in Portugal was planning to incorporate an industrial HRC system. According to the company's information, the workstation deals with recurrent absenteeism due to WMSD. Therefore, this framework was applied in order to help decision makers which system to choose. Figure   Figure 3. Phases of the developed framework.

Case Study
A workstation for the manual assembly of electrical components in Portugal was planning to incorporate an industrial HRC system. According to the company's information, the workstation deals with recurrent absenteeism due to WMSD. Therefore, this framework was applied in order to help decision makers which system to choose. Figure 4 presents the steps towards the objectives of this case study.

Objectives
Absenteeism is often the starting point of ergonomic interventions in a production line, as it directly impacts in human being as well as productivity [59]. According to the company, the main goal is to reduce absenteeism due to WMSD by 25% and simultaneously increase productivity by 20%. In order to prospect scenarios with HRC systems, a computational model was used in this case study. Therefore, a workstation assessment regarding ergonomics and productivity was carried out to design and to feed the model.

Ergonomic Work Analysis
Direct observation was performed on the workstation under analysis in order to characterize the environment, the process, and the main variables as well as their interconnections. The most critical WMSD risk factors were identified, such as awkward postures and repetitive tasks. Therefore, a replica-workstation for manual assembly of electronic components was built considering the same dimensions in order to replicate body ranges and postures.
The workstation layout is presented in Figure 5. and consists of reaching pieces on both the left and right sides at different heights, and to assemble them in the center (in front of the worker). Thus, the characteristic movements are vertical reach movements at shoulder level, sometimes with lateral flexion of the trunk and displacement in steps to the side.

Objectives
Absenteeism is often the starting point of ergonomic interventions in a production line, as it directly impacts in human being as well as productivity [59]. According to the company, the main goal is to reduce absenteeism due to WMSD by 25% and simultaneously increase productivity by 20%. In order to prospect scenarios with HRC systems, a computational model was used in this case study. Therefore, a workstation assessment regarding ergonomics and productivity was carried out to design and to feed the model.

Ergonomic Work Analysis
Direct observation was performed on the workstation under analysis in order to characterize the environment, the process, and the main variables as well as their interconnections. The most critical WMSD risk factors were identified, such as awkward postures and repetitive tasks. Therefore, a replica-workstation for manual assembly of electronic components was built considering the same dimensions in order to replicate body ranges and postures.
The workstation layout is presented in Figure 5. and consists of reaching pieces on both the left and right sides at different heights, and to assemble them in the center (in front of the worker). Thus, the characteristic movements are vertical reach movements at shoulder level, sometimes with lateral flexion of the trunk and displacement in steps to the side. In order to obtain data from all different dimensions of the Portuguese population, the worker's sample was divided in three classes of percentiles. In total, six workers were chosen (one female and one male from each of the percentile classes presented in Table 4). The anthropometric dimension "stature" was the dimension used as reference [60]. RULA was the ergonomic method chosen for physical workload assessment. Workers were equipped with wearable sensors and RULA score was obtained automatically with an algorithm developed by Xsens ® . The algorithm estimates the percentage of time that the subject was in each of the four risk levels defined by RULA during the task. The The piece to be assembled consists of seven components which are presented in Figure 6. From P1 to P3, the components are on the left side, from P4 to P6, on the right side, and P7 is placed in front of the worker (close to the assembly table). In order to obtain data from all different dimensions of the Portuguese population, the worker's sample was divided in three classes of percentiles. In total, six workers were chosen (one female and one male from each of the percentile classes presented in Table 4). The anthropometric dimension "stature" was the dimension used as reference [60]. RULA was the ergonomic method chosen for physical workload assessment. Workers were equipped with wearable sensors and RULA score was obtained automatically with an algorithm developed by Xsens ® . The algorithm estimates the percentage of time that the subject was in each of the four risk levels defined by RULA during the task. The In order to obtain data from all different dimensions of the Portuguese population, the worker's sample was divided in three classes of percentiles. In total, six workers were chosen (one female and one male from each of the percentile classes presented in Table 4). The anthropometric dimension "stature" was the dimension used as reference [60]. RULA was the ergonomic method chosen for physical workload assessment. Workers were equipped with wearable sensors and RULA score was obtained automatically with an algorithm developed by Xsens ® . The algorithm estimates the percentage of time that the subject was in each of the four risk levels defined by RULA during the task. The procedure involved the following steps: sociodemographic data (age, height, weight) of the six workers were collected; body dimensions were collected as input to the software; biomechanical data were collected in real time of work activity [61] at a frequency of 120 Hz [62]; data were transmitted wirelessly to a computer loaded with a software that allows the movements to be observed, recorded, and analyzed; the speed of the movements was defined by each participant, according to their individual abilities; workers performed the task continuously during six cycle times. Figure 7 shows the sensors placement that was followed by calibration and recording using Xsens ® MVN software. Figure 8 shows workers performing the assembly tasks with the sensors capturing movements.
Safety 2021, 7, x FOR PEER REVIEW 12 of 21 procedure involved the following steps: sociodemographic data (age, height, weight) of the six workers were collected; body dimensions were collected as input to the software; biomechanical data were collected in real time of work activity [61] at a frequency of 120 Hz [62]; data were transmitted wirelessly to a computer loaded with a software that allows the movements to be observed, recorded, and analyzed; the speed of the movements was defined by each participant, according to their individual abilities; workers performed the task continuously during six cycle times. Figure 7 shows the sensors placement that was followed by calibration and recording using Xsens ® MVN software. Figure 8 shows workers performing the assembly tasks with the sensors capturing movements. The average percentages of time exposure obtained for each risk level were the following: 0.2% acceptable, 77.5% low risk, 19.3% medium risk, and 3.0% high risk, which was considered RULA = 4. Therefore, according to Table 1, the model used "Physical Risk Level II" for simulation purposes. procedure involved the following steps: sociodemographic data (age, height, weight) of the six workers were collected; body dimensions were collected as input to the software; biomechanical data were collected in real time of work activity [61] at a frequency of 120 Hz [62]; data were transmitted wirelessly to a computer loaded with a software that allows the movements to be observed, recorded, and analyzed; the speed of the movements was defined by each participant, according to their individual abilities; workers performed the task continuously during six cycle times. Figure 7 shows the sensors placement that was followed by calibration and recording using Xsens ® MVN software. Figure 8 shows workers performing the assembly tasks with the sensors capturing movements. The average percentages of time exposure obtained for each risk level were the following: 0.2% acceptable, 77.5% low risk, 19.3% medium risk, and 3.0% high risk, which was considered RULA = 4. Therefore, according to Table 1, the model used "Physical Risk Level II" for simulation purposes. The average percentages of time exposure obtained for each risk level were the following: 0.2% acceptable, 77.5% low risk, 19.3% medium risk, and 3.0% high risk, which was considered RULA = 4. Therefore, according to Table 1, the model used "Physical Risk Level II" for simulation purposes.

Mental Workload Assessment
The ergonomic method chosen for mental workload assessment was NASA-TLX [49]. The questionnaire was applied to the six volunteers after a work shift. The procedure was: to explain the definition of the six items (mental demand, physical demand, temporal demand, performance, effort, and frustration); to ask the workers to choose in a scale their own perception on the task; and finally, to ask each worker to choose the item they consider most relevant for each of the fifteen pairs of items presented. The average results are shown in Figure 9b. The overall workload is 44, which means "moderate", and no item is especially critical. It was classified as low mental risk level according to Table 2. Therefore, the model considers "Mental Risk Level I" for simulation purposes.

Mental Workload Assessment
The ergonomic method chosen for mental workload assessment was NASA-TLX [49]. The questionnaire was applied to the six volunteers after a work shift. The procedure was: to explain the definition of the six items (mental demand, physical demand, temporal demand, performance, effort, and frustration); to ask the workers to choose in a scale their own perception on the task; and finally, to ask each worker to choose the item they consider most relevant for each of the fifteen pairs of items presented. The average results are shown in Figure 9b. The overall workload is 44, which means "moderate", and no item is especially critical. It was classified as low mental risk level according to Table 2. Therefore, the model considers "Mental Risk Level I" for simulation purposes.

Cycle Time
Time study is used to determine the time required by a qualified and well-trained person working at a normal pace to do a specified task [63]. A digital stopwatch was used to record 20 minutes of task executions. The average cycle time in the replica-workstation was 23.10s.
The results obtained in the assessment phase are summarized in Table 5. to be used in the modeling and simulation phases.

Modeling
The use of system dynamics allows the systemic contextualization of factors assessed in the workstation, a mathematical description of their interrelationship, and the prospection of scenarios in non-linear relationships. A system dynamics was designed in the Vensim software [64] to model the whole environment (See Appendix A- Figure A1 and Table A1). The production system was based on the model described by Sterman [24] (e.g., "work in process", "production start rate" and "production rate"). Furthermore, the construction of the model was adapted to the characteristics of the context under study. In this sense, the "work for process" is a complement factor inserted along with its input rate and influence variables. The structure of employees' dynamics and its reflection on the production system follows an analogous context [15], in which the same phenomenon is analyzed. The model calibration was performed aiming to drive the system's behavior to

Cycle Time
Time study is used to determine the time required by a qualified and well-trained person working at a normal pace to do a specified task [63]. A digital stopwatch was used to record 20 min of task executions. The average cycle time in the replica-workstation was 23.10 s.
The results obtained in the assessment phase are summarized in Table 5. to be used in the modeling and simulation phases.

Modeling
The use of system dynamics allows the systemic contextualization of factors assessed in the workstation, a mathematical description of their interrelationship, and the prospection of scenarios in non-linear relationships. A system dynamics was designed in the Vensim software [64] to model the whole environment (See Appendix A- Figure A1 and Table A1). The production system was based on the model described by Sterman [24] (e.g., "work in process", "production start rate" and "production rate"). Furthermore, the construction of the model was adapted to the characteristics of the context under study. In this sense, the "work for process" is a complement factor inserted along with its input rate and influence variables. The structure of employees' dynamics and its reflection on the production system follows an analogous context [15], in which the same phenomenon is analyzed. The model calibration was performed aiming to drive the system's behavior to values close to those observed and recorded for the real system. Therefore, certain equations receive an adjustment constant (e.g., equation Time to gain Knowledge). The "Time to gain Knowledge (TK)" is inversely proportional to the "Cycle Frequency (CF)" since frequent cycles (with greater repetition) result in less time to gain knowledge (faster learning through repetition), being mathematically formulated as: TK = constant/CF. The constant value calibrated in this study was 480, resulting in TK = 480/CF. The variables "Mental Overload" (MO), "Physical Overload" (PO), and "Knowledge Required" (KR) were entered as a function of the LoC ( Table 6). The variable "Postural Requirement" (PR) is an intermediate variable between the LoC and PO for the model to support the insertion of values obtained from any ergonomic method chosen to assess physical effort. The integration increment was set to 1/5 of the value of the shortest time constant. Since the analysis is focused on operational dynamics whose effects occur in a relatively short time, the simulation horizon was defined in a work shift (8 hours). The total human workload value was subjectively weighted by the modelers at 33% for mental demands (Risk Level I) and 67% for physical demands (Risk Level II), which means that physical factors were considered twice more critical in the current study. In addition, a very important aspect for this analysis is the insertion of the LoC influences in the simulation model. In this case, the insertion was done by the definition of values in ordinal scales of five levels, in which level 0 of collaboration means a workstation without robot, and level 4 means maximum collaboration. As LoC influences on mental workload, postural demand, and knowledge necessary to perform the task, the values defined by the modelers to configure the model were assessed and are shown in Table 6. The relationship among LoC, cycle time, and sick leaves over the simulation horizon is presented in Table 7 which turns out to be very useful for the interpretation of the system. The simulation shows that the most favorable scenario to avoid the incidence of occupational diseases, which leads to employees on leave, corresponds to the highest LoC. Although this level means a greater frequency of cycles (due to the reduced cycle time) and the performance of mental workload, the reduction in postural load seems to compensate for these factors. The results demonstrate that inserting an industrial HRC system is complex. Both physical and mental workload may affect leave rate and productivity, which significantly change the prospective scenarios.

Technical Evaluation
Some critical checks were performed in order to verify if an industrial HRC system was technically practicable in the workstation with the mentioned characteristics. The main possibilities for a cobot to be useful are delivering, handling, and assembling. In this regard, the main complexities are those related to the work pieces: geometry, dimensions, loads, and materials; and related to the organization: assembly location, layout, ranges, and sequences [65]. The critical technical issues follow the list presented in [14].

•
Delivering: without perceived restrictions; • Handling: three small components would be difficult for a cobot to manipulate; • Assembly: some components are resistant to insertion, the assembly table is overconstrained, the assembly process demands reorientation of previous assembled components, and components must be compressed during assembly.
In summary, the small pieces must be manipulated by the worker and most of the assembly tasks require human skills to guarantee quality in the production line. Other tasks are considered feasible in terms of robot execution. For example, it is technically capable of handling and delivering the four larger pieces.

Ergonomic Evaluation
An ergonomic intervention in the workstation by implementing an industrial HRC system will definitively change organizational, physical, and cognitive domains. In the new layout, work pieces must be delivered to the worker in a location closer to him, thus avoiding inappropriate postures on his part. By decreasing RULA scores, less WMSD is expected in the production line [66]. On the other hand, mental workload also increases with repeatability of movements as psychological well-being depends on production flexibility [67].
The simulations showed the behavior of the system regarding the occurrence of sick leave and productivity for different LoC. Task subdivision increased worker's skills to accomplish the tasks and directly increased productivity. However, the influence of physical overload on output showed that productivity is also associated with stereotyped movement, which is often consequence of repeatability or excessive task subdivision. Therefore, it is not recommended to overpush. Moreover, physiological recovery is needed in order to decrease the influence of physical overload on leave rate [15].
In general, the inclusion of an industrial HRC system is beneficial for ergonomics. The main attention points are related to the new pace of work and the division of tasks between the worker and the cobot. These factors are better discussed in the economic evaluation.

Safety Evaluation
Safety requirements: Due to the risk of collisions, a robot working without fences represents a risk in the HRC system [68]. Safety management, sensors, and motion planning and control are the most important to avoid collisions [27]. According to ISO/TS 15,066 [56], the safety methods may vary for different LoC: safety-rated monitored stop, hand guiding, speed and separation monitoring, and power and force limiting.
Risk assessment: Similar assembly application with HRC evaluated as low risk for the human [9]. However, an evaluation is important due to individual scenarios of a collaboration [69]. In this work, it has been considered that this evaluation does not change as the light-weighted pieces are delivered always at the same place by the robot, and safety requirements may follow robot's capabilities in different LoC, according to Table 3.

Economic Evaluation
Considering the above mentioned regarding technical, ergonomic, and safety evaluations, an economic evaluation aims at recognizing the tasks for which a cobot provides value to the production process. As stated in [14], it is important to integrate organizational and economic factors for a better collaborative workstation. Table 8 presents the division of the workstation tasks between human and cobot. In this case study, the assembly tasks must be performed by the worker, which means that all the value-added tasks were allocated to the worker. However, the cobot is capable of handling and delivering P1, P4, P5, and P6 (larger pieces), which means time saving in the production line. Moreover, according to the prospections of the model, less absenteeism is expected due to WMSD, which means maintaining human resources and knowledge.

Decision
Considering the objectives, assessments, simulations, evaluations, and investments to be made, the workstation condition justifies the implementation of an industrial HRC system. It has been taken into account that by thinking on the whole system, higher levels of automation do not necessarily guarantee higher levels of productivity, safety, or operator well-being [70]. Higher LoC may induce other problems, such as confusion, complacency, or loss of certain skills. According to the classification suggested in [55], in order to achieve management expectations in short-and long-terms, a meaningful decision would be upon a cobot with capabilities in the Level 3 presented in Table 3.
Here are some final considerations. The new workstation must: ensure the cobot arm to reach the work pieces; consider workers anthropometric dimensions in the layout design; choose adequate technologies for gripping, safety, and recognition and awareness sensors.

Conclusions
Inserting a cobot into an assembly line is a complex decision. The dynamics of the elements that interact with each other can result in counter-intuitive effects. In this context, decision-making needs to be carefully thought out and developed. The current framework addresses the challenge of implementing a cobot. Ergonomics, safety, and productivity aspects were discussed based on system dynamics in a context of industrial assembly process. The system behavior and its prospective scenarios with HRC were key to managerial evaluation. Finally, a decision was made based on reliable data instead of a mental model, drastically reducing the risk of failure when deciding upon a HRC system. Therefore, this framework was presented to improve managerial decision assertiveness. Future work is intended to be developed in order to simplify the model and maintaining its usefulness. We also recommend the practical application of this framework in other assembly lines that meet the demand for the insertion of a cobot.

Data Availability Statement:
The Appendix A is the source of data generated during the study. Researchers can replicate the results by following the Stock and Flow Map together with equations.

Conflicts of Interest:
The authors declare no conflict of interest. The sponsors had no role in the design, execution, interpretation, or writing of the study. Table A1. Model equations.

Label Equation Unit
Work for process (WFP) WFP =