Enabling Knowledge Discovery in Multi-Objective Optimizations of Worker Well-Being and Productivity

Usually, optimizing productivity and optimizing worker well-being are separate tasks performed by engineers with different roles and goals using different tools. This results in a silo effect which can lead to a slow development process and suboptimal solutions, with one of the objectives, either productivity or worker well-being, being given precedence. Moreover, studies often focus on finding the best solutions for a particular use case, and once solutions have been identified and one has been implemented, the engineers move on to analyzing the next use case. However, the knowledge obtained from previous use cases could be used to find rules of thumb for similar use cases without needing to perform new optimizations. In this study, we employed the use of data mining methods to obtain knowledge from a real-world optimization dataset of multi-objective optimizations of worker well-being and productivity with the aim to identify actionable insights for the current and future optimization cases. Using different analysis and data mining methods on the database revealed rules, as well as the relative importance of the design variables of a workstation. The generated rules have been used to identify measures to improve the welding gun workstation design.


Introduction
Simulation is widely used in industries such as the automotive industry because it can efficiently create, test, and optimize the design of products and production systems in the virtual world without any need to create, test, and optimize prototypes in the physical world [1][2][3]. Using simulation saves time and money and allows a more thorough investigation of the solution space. These improvements in productivity have led to simulation being used in the design of workstations [4,5]. Simulation tools using digital human modeling (DHM) are also used in designing workstations to assess workers' wellbeing [6].
However, simulations for optimizing productivity and simulations for assessing worker well-being are usually performed by people with different roles (production engineers solve productivity problems, while ergonomics specialists consider worker wellbeing). They have different focuses or objectives and use different tools. This can create silo effects, leading to slow development processes and suboptimal solutions.
Productivity and worker well-being often go hand in hand, because improving working conditions usually increases productivity [7][8][9][10]. However, sometimes the goals of productivity and well-being are at odds. Companies need to find and implement solutions in their production facilities to maintain profitability, output, and quality, as well as the well-being of workers. Previous research has identified the core elements of DHM tools and suggested a structured process for applying DHM tools in the design and development the study, based on key anthropometric variables aimed to represent the anthropometric diversity of the workers at the factory, i.e., stature and elbow height ( Figure 1). The virtual model was modeled in the DHM tool IPS IMMA [21] by representing the welding gun workstation by imported CAD geometries. A family of 14 manikins (7 female, 7 male) was created by using an anthropometric module [22] in the DHM tool used in the study, based on key anthropometric variables aimed to represent the anthropometric diversity of the workers at the factory, i.e., stature and elbow height ( Figure 1). The welding process for each welding gun at each welding spot was simulated for the entire manikin family. Table 1 shows the values for the stature and elbow height of the manikins simulated for the use case.

Well-Being and Productivity Evaluation
We evaluated the productivity of the welding workstation by measuring the cycle time of the welding sequence from the simulation. The cycle time was calculated by considering the actions of the workers performing the sequence. These actions included both value-adding operations (e.g., the time to weld a welding spot) and non-value-adding The welding process for each welding gun at each welding spot was simulated for the entire manikin family. Table 1 shows the values for the stature and elbow height of the manikins simulated for the use case.

. Well-Being and Productivity Evaluation
We evaluated the productivity of the welding workstation by measuring the cycle time of the welding sequence from the simulation. The cycle time was calculated by considering the actions of the workers performing the sequence. These actions included both value-adding operations (e.g., the time to weld a welding spot) and non-value-adding operations (e.g., changing the welding gun, changing the welding side, and moving between welding spots).
We used the Rapid Upper Limb Assessment (RULA) method [23] to evaluate the risk of WMSDs. We considered RULA appropriate since the work stresses mainly involve postural stresses on the upper limbs, as the weight of the welding guns is supported by a lifting device. For each posture assessed, RULA gives a risk score from 1 to 7, which results in four action levels. A score of 1 to 2 is acceptable, a score of 3 to 4 suggests changes may be required, 5 to 6 indicates changes will soon be required, and 7 indicates that changes are required immediately.

Mathematical Modeling of Optimization
The optimization model considers both worker well-being and productivity objectives, that is, it is a multi-objective optimization model. The indices, parameters, variables, and objectives of the optimization model are shown in Table 2. We used the Rapid Upper Limb Assessment (RULA) method [23] to evaluate the risk of WMSDs. We considered RULA appropriate since the work stresses mainly involve postural stresses on the upper limbs, as the weight of the welding guns is supported by a lifting device. For each posture assessed, RULA gives a risk score from 1 to 7, which results in four action levels. A score of 1 to 2 is acceptable, a score of 3 to 4 suggests changes may be required, 5 to 6 indicates changes will soon be required, and 7 indicates that changes are required immediately.

Mathematical Modeling of Optimization
The optimization model considers both worker well-being and productivity objectives, that is, it is a multi-objective optimization model. The indices, parameters, variables, and objectives of the optimization model are shown in Table 2.

Objectives
RULA score for a manikin on a side with a welding gun at a welding spot Cycle time of welding process (s) Average RULA score per manikin in the welding process Average RULA score of all manikins in the welding process A multi-objective optimization of the average of the RULA scores for all the manikins and the cycle time was performed. The mean RULA score of the 14 manikins was calculated as a mean value for all the manikins using all welding spots. The risk of WMSDs is therefore calculated by a single objective: The cycle time was calculated as the sum of the welding time at each welding spot and the time to change welding gun, welding spot, and welding side: Average RULA score of all manikins in the welding process A multi-objective optimization of the average of the RULA scores for all the manikins and the cycle time was performed. The mean RULA score of the 14 manikins was calculated as a mean value for all the manikins using all welding spots. The risk of WMSDs is therefore calculated by a single objective:

MI N
Sustainability 2022, 14, x FOR PEER REVIEW operations (e.g., changing the welding gun, changing the welding tween welding spots).
We used the Rapid Upper Limb Assessment (RULA) method [2 of WMSDs. We considered RULA appropriate since the work stresse tural stresses on the upper limbs, as the weight of the welding gu lifting device. For each posture assessed, RULA gives a risk score from in four action levels. A score of 1 to 2 is acceptable, a score of 3 to 4 s be required, 5 to 6 indicates changes will soon be required, and 7 i are required immediately.

Mathematical Modeling of Optimization
The optimization model considers both worker well-being an tives, that is, it is a multi-objective optimization model. The indices, p and objectives of the optimization model are shown in Table 2. A multi-objective optimization of the average of the RULA scor and the cycle time was performed. The mean RULA score of the 14 lated as a mean value for all the manikins using all welding spots. T therefore calculated by a single objective: The cycle time was calculated as the sum of the welding time and the time to change welding gun, welding spot, and welding sid The cycle time was calculated as the sum of the welding time at each welding spot and the time to change welding gun, welding spot, and welding side:

Optimization Method
Many real-world optimization problems involve multiple conflicting objectives and therefore lead to multiple Pareto-optimal solutions, that is, the optimal solution for one optimization objective may not be optimal for another objective [24]. Therefore, solving a MOOP is about balancing the trade-offs among conflicting objectives. This complexity also affects the process of obtaining optimal solutions while performing optimization. Due to their population-based nature, multi-objective evolutionary algorithms are often used for solving MOOPs. The evolutionary algorithm NSGA-II was used to optimize this application study because of its efficiency in multi-objective optimizations [25]. The parameters used for the optimization are presented in Table 3.

Knowledge Discovery
We sought to extract general knowledge from the solutions to the MOOP. A MOOP solution can be seen as inhabiting two different spaces: the decision space, where the decision variables exist, and the objective space where the objective values exist. While optimization algorithms drive the search to find higher performing solutions in the objective space, data mining methods can be used to find knowledge about preferred solutions in terms of the decision space. To generate knowledge, we defined the preference by analyzing the solutions, and then applied a data mining method to generate decision rules to describe those preferences. The process of knowledge discovery is done in several steps: (1) data filtering, (2) data clustering, (3) data visualization, (4) rule extraction, and (5) knowledge interpretation.
The first step, data filtering, was done by removing the duplicates in the database to remove any bias in the resulting knowledge ( Figure 2). In the second step, data clustering, we identified four principal objectives: lowest cycle time, lowest average RULA score for all manikins, balance between average RULA score for all manikins and cycle time, and worker diversity inclusion rules. In the next step, data visualization, we analyzed the database by using parallel coordinates, boxplots, and 2D and 3D scatter plots. Later, knowledge discovery was performed using influence score by rank (InfS-R) and flexible pattern mining (FPM) methods [26]. In the final step, knowledge interpretation, we analyzed the generated rules and related them to the specific implications in the application study to generate reproducible knowledge for future workstation designs.    [27] is a recently developed method for obtaining decision rules in the decision space, based on selections made by the decision maker in the objective space. FPM is an extension of sequential pattern mining [28] and can find rules of the forms {x i < c 1 }, {x i > c 2 }, and {x i = c 3 } for a decision variable x i and constant values c 1 , c 2 , and c 3 . The FPM procedure finds decision rules about the decision space that describes selections made in the objective space. To run the process, the decision maker must supply a set of selected solutions and a set of unselected solutions. The resulting FPM rules can be seen as a tuple containing the following elements: label, referring to the name of the decision variable; sign, referring to one of the signs <, >, or = ; value, referring to the constant value of the rule; sig., referring to the significance of the rule, that is, the percentage of the solutions in the selected set that the rule covers; and unsig., referring to the unselected significance (unsignificance), that is, the percentage of solutions in the unselected set that the rule covers. Note that since the significance and unsignificance can be found simply by counting the solutions covered by a single rule, the same can be done for interactions of two or more rules. In Section 3 we show three levels of FPM rule interactions. An interesting and descriptive FPM rule would have high significance and low unselected significance. The ratio of significance over unselected significance serves as a general metric for a good rule.

InfS-R
Factor screening to find which variable has the most influence in a model can be a useful tool for understanding the model and how the variables interact [29]. Recently, an approach for finding influential variables in the resulting solution set from a multi-objective optimization problem was presented in [26]. The approach defines two methods. One finds an influence score for the decision variables in terms of how they affect the solutions on the Pareto-optimal front; this method is called influence score by Pareto front (InfS-P). The other method finds an influence score for decision variables in terms of how they affect the convergence on the Pareto front by considering the rank of the solutions after non-dominated sorting, called InfS-R. Both methods divide the solutions into different selections and find FPM rules to describe the differences between them. The frequency of and significance of the different variables is then used to determine how influential they are. We use InfS-R to find the relative influence of the contribution of certain variables to generating Pareto-optimal solutions.

Results
The solution space from the optimization is presented in a 2D scatter plot ( Figure 3) of the colliding objectives of average RULA for all manikins and cycle time. are. We use InfS-R to find the relative influence of the contribution of certain variables to generating Pareto-optimal solutions.

Results
The solution space from the optimization is presented in a 2D scatter plot ( Figure 3) of the colliding objectives of average RULA for all manikins and cycle time. The solutions that correspond to the non-dominated solutions represent the Pareto front for this optimization (marked by "+" in Figure 3) and are the best solutions of the optimization. These solutions are shown in Table 4.  We used the Rapid Upper Limb Assessm of WMSDs. We considered RULA appropriate tural stresses on the upper limbs, as the weig lifting device. For each posture assessed, RULA in four action levels. A score of 1 to 2 is accept be required, 5 to 6 indicates changes will soo are required immediately.

Mathematical Modeling of Optimization
The optimization model considers both tives, that is, it is a multi-objective optimizatio and objectives of the optimization model are s The solutions that correspond to the non-dominated solutions represent the Pareto front for this optimization (marked by "+" in Figure 3) and are the best solutions of the optimization. These solutions are shown in Table 4. Table 4. Results from optimization.

Result Selected CT (s)
Welding gun used at each welding spot Welding side at each welding spot Objectives R w Cycle time of welding process (s) Average RULA score per manikin in the welding process Average RULA score of all manikins in the welding process A multi-objective optimization of the average of t and the cycle time was performed. The mean RULA lated as a mean value for all the manikins using all w therefore calculated by a single objective: The cycle time was calculated as the sum of the and the time to change welding gun, welding spot, an Sequence Lowest CT 47 3.09 Compromise between CT and

Objectives
RULA score for a man with a welding gun at Cycle time of welding process (s) Average RULA score per manikin in the welding process Average RULA score of all manikins in the welding process A multi-objective optimization of the average of the RULA scores for a and the cycle time was performed. The mean RULA score of the 14 mani lated as a mean value for all the manikins using all welding spots. The ris therefore calculated by a single objective: The cycle time was calculated as the sum of the welding time at eac and the time to change welding gun, welding spot, and welding side:

Objectives
RULA score for a manikin on a s with a welding gun at a welding Cycle time of welding process (s) Average RULA score per manikin in the welding process Average RULA score of all manikins in the welding process A multi-objective optimization of the average of the RULA scores for all the man and the cycle time was performed. The mean RULA score of the 14 manikins was c lated as a mean value for all the manikins using all welding spots. The risk of WMS therefore calculated by a single objective: The cycle time was calculated as the sum of the welding time at each welding and the time to change welding gun, welding spot, and welding side: Spot sequence: 4-5-6-7-1-2-3 Gun sequence: 4-5-6-7-1-2-3 Side sequence: 2-2-2-1-1-1-1

Data Filtering
The database of the application study was filtered by removing duplicates. Also, in order to help the data mining process, related variables were merged. In this case, the welding guns and the welding sides were converted into a single variable "gun and side" before proceeding with the data mining. To convert gun and side to the same variable, the variables were converted into a single integer per spot, GS w . Due to the availability of the guns and sides in different welding spots, the range was different for every welding spot (Table 5). Table 5. Corresponding value for every GS w depending on the gun and side availability at each spot.

Data Clustering
The selected clusters that need to be analyzed represent the (1)  operations (e.g., changing t tween welding spots). We used the Rapid Up of WMSDs. We considered R tural stresses on the upper lifting device. For each postu in four action levels. A scor be required, 5 to 6 indicates are required immediately.

Mathematical Modelin
The optimization mod tives, that is, it is a multi-obj and objectives of the optimi  operations (e.g., changing the welding gun, changing the welding side, and moving between welding spots). We used the Rapid Upper Limb Assessment (RULA) method [23] to evaluate the risk of WMSDs. We considered RULA appropriate since the work stresses mainly involve postural stresses on the upper limbs, as the weight of the welding guns is supported by a lifting device. For each posture assessed, RULA gives a risk score from 1 to 7, which results in four action levels. A score of 1 to 2 is acceptable, a score of 3 to 4 suggests changes may be required, 5 to 6 indicates changes will soon be required, and 7 indicates that changes are required immediately.

Mathematical Modeling of Optimization
The optimization model considers both worker well-being and productivity objectives, that is, it is a multi-objective optimization model. The indices, parameters, variables, and objectives of the optimization model are shown in Table 2.

Objectives
RULA score for a manikin on a side with a welding gun at a welding spot Cycle time of welding process (s) Average RULA score per manikin in the welding process Average RULA score of all manikins in the welding process A multi-objective optimization of the average of the RULA scores for all the manikins and the cycle time was performed. The mean RULA score of the 14 manikins was calculated as a mean value for all the manikins using all welding spots. The risk of WMSDs is therefore calculated by a single objective: The cycle time was calculated as the sum of the welding time at each welding spot and the time to change welding gun, welding spot, and welding side: were defined by selecting the solutions with the lowest scores. The cluster of (3) (balance between CT and Sustainability 2022, 14, x FOR PEER REVIEW operations (e.g., changing the welding gun, changing the welding side tween welding spots).
We used the Rapid Upper Limb Assessment (RULA) method [23] t of WMSDs. We considered RULA appropriate since the work stresses m tural stresses on the upper limbs, as the weight of the welding guns lifting device. For each posture assessed, RULA gives a risk score from 1 in four action levels. A score of 1 to 2 is acceptable, a score of 3 to 4 sugg be required, 5 to 6 indicates changes will soon be required, and 7 indic are required immediately.

Mathematical Modeling of Optimization
The optimization model considers both worker well-being and p tives, that is, it is a multi-objective optimization model. The indices, para and objectives of the optimization model are shown in Table 2. A multi-objective optimization of the average of the RULA scores fo and the cycle time was performed. The mean RULA score of the 14 ma lated as a mean value for all the manikins using all welding spots. The therefore calculated by a single objective: The cycle time was calculated as the sum of the welding time at e ) was defined as all solutions that had a lower CT than cluster (2), and lower Sustainability 2022, 14, x FOR PEER REVIEW operations (e.g., changing the welding gun, changing the welding side, and mo tween welding spots).
We used the Rapid Upper Limb Assessment (RULA) method [23] to evaluate of WMSDs. We considered RULA appropriate since the work stresses mainly invo tural stresses on the upper limbs, as the weight of the welding guns is support lifting device. For each posture assessed, RULA gives a risk score from 1 to 7, which in four action levels. A score of 1 to 2 is acceptable, a score of 3 to 4 suggests chan be required, 5 to 6 indicates changes will soon be required, and 7 indicates that are required immediately.

Mathematical Modeling of Optimization
The optimization model considers both worker well-being and productivit tives, that is, it is a multi-objective optimization model. The indices, parameters, v and objectives of the optimization model are shown in Table 2. Objectives RULA score for a manikin on with a welding gun at a weldi Cycle time of welding process (s) Average RULA score per manikin in the welding process Average RULA score of all manikins in the welding process A multi-objective optimization of the average of the RULA scores for all the m and the cycle time was performed. The mean RULA score of the 14 manikins wa lated as a mean value for all the manikins using all welding spots. The risk of W therefore calculated by a single objective: = ∑ ∑ than cluster (1) (Figure 4).
The selected clusters that need to be analyzed represent the (1) lowest cycle time ( ) (marked by "+"), (2) lowest average RULA score for all manikins ( ) (marked by "-"), (3) the balance between average RULA score for all manikins and cycle time (marked by "x"), and (4) worker diversity inclusion rules. The clusters of (1) lowest and (2) lowest were defined by selecting the solutions with the lowest scores. The cluster of (3) (balance between and ) was defined as all solutions that had a lower than cluster (2), and lower than cluster (1) (Figure 4). An initial analysis of the manikins was necessary to select the cluster for (4) (diversity inclusion rules) since the critical manikins were still not identified.

Data Visualization
A boxplot analysis of the average RULA score of each manikin ( ) was made to evaluate the critical manikins and consider worker anthropometric diversity ( Figure 5).    An initial analysis of the manikins was necessary to select the cluster for (4) (diversity inclusion rules) since the critical manikins were still not identified.

Data Visualization
A boxplot analysis of the average RULA score of each manikin (ER m ) was made to evaluate the critical manikins and consider worker anthropometric diversity ( Figure 5).  Manikins 13 and 14 have higher average RULA scores than other manikins for all solutions ( Figure 5). These are the manikins with the lowest stature in the population considered ( Table 1). The cluster for (4) (diversity inclusion rules) was therefore formed by the best solutions for the of manikins 13 and 14, and (marked by "x" in Figure 6). Manikins 13 and 14 have higher average RULA scores than other manikins for all solutions ( Figure 5). These are the manikins with the lowest stature in the population considered ( Table 1). The cluster for (4) (diversity inclusion rules) was therefore formed by the best solutions for the ER m of manikins 13 and 14, ER 13 and ER 14 (marked by "x" in Figure 6). Manikins 13 and 14 have higher average RULA scores than other manikins for all solutions ( Figure 5). These are the manikins with the lowest stature in the population considered ( Table 1). The cluster for (4) (diversity inclusion rules) was therefore formed by the best solutions for the of manikins 13 and 14, and (marked by "x" in Figure 6).

Knowledge Discovery
To obtain rules from the results of the application study, InfS-R [26] was used to determine the relative importance of the variables in the optimization. The first InfS-R ranked the non-dominated sorting for the objectives and . The variables studied were the selection of welding gun and side for every welding spot ( or SpotW G/S in figures) and the welding sequence ( ) (Figure 7).

Knowledge Discovery
To obtain rules from the results of the application study, InfS-R [26] was used to determine the relative importance of the variables in the optimization. The first InfS-R ranked the non-dominated sorting for the objectives CT and are required immediately.

Mathematical Modeling of Optimiz
The optimization model considers tives, that is, it is a multi-objective optimi and objectives of the optimization mode

Objectives
Cycle time of welding process (s) Average RULA score per manikin in the welding process Average RULA score of all manikins in the welding process A multi-objective optimization of th and the cycle time was performed. The m lated as a mean value for all the manikin therefore calculated by a single objective = The cycle time was calculated as th and the time to change welding gun, we . The variables studied were the selection of welding gun and side for every welding spot (GS w or SpotW G/S in figures) and the welding sequence (X w ) (Figure 7).     After the InfS-R analysis, an FPM analysis was run for the four clusters (1) low cycle time (marked by "+" in Figure 4), (2) lowest average RULA score for all manik (marked by "-" in Figure 4), (3) balance between average RULA score for all manikins a cycle time (marked by "x" in Figure 4), and (4) worker diversity inclusion rules (mark by "x" in Figure 6). The FPM analysis was run for three levels of rule-interactions and minimum significance and minimum support. After that, the rules were filtered to higher ratio of significance/unsignificance to obtain the most relevant rules for the f clusters. The significance, unsignificance, and ratio of the obtained rules in the four cl ters are presented in Table 6.  After the InfS-R analysis, an FPM analysis was run for the four clusters (1) lowest cycle time (marked by "+" in Figure 4), (2) lowest average RULA score for all manikins (marked by "-" in Figure 4), (3) balance between average RULA score for all manikins and cycle time (marked by "x" in Figure 4), and (4) worker diversity inclusion rules (marked by "x" in Figure 6). The FPM analysis was run for three levels of rule-interactions and 0.5 minimum significance and minimum support. After that, the rules were filtered to the higher ratio of significance/unsignificance to obtain the most relevant rules for the four clusters. The significance, unsignificance, and ratio of the obtained rules in the four clusters are presented in Table 6. operations (e.g., changing the welding gun, changing the welding side, and moving between welding spots). We used the Rapid Upper Limb Assessment (RULA) method [23] to evaluate the risk of WMSDs. We considered RULA appropriate since the work stresses mainly involve postural stresses on the upper limbs, as the weight of the welding guns is supported by a lifting device. For each posture assessed, RULA gives a risk score from 1 to 7, which results in four action levels. A score of 1 to 2 is acceptable, a score of 3 to 4 suggests changes may be required, 5 to 6 indicates changes will soon be required, and 7 indicates that changes are required immediately.

Mathematical Modeling of Optimization
The optimization model considers both worker well-being and productivity objectives, that is, it is a multi-objective optimization model. The indices, parameters, variables, and objectives of the optimization model are shown in Table 2.

Objectives
RULA score for a manikin on a side with a welding gun at a welding spot Cycle time of welding process (s) Average RULA score per manikin in the welding process Average RULA score of all manikins in the welding process A multi-objective optimization of the average of the RULA scores for all the manikins and the cycle time was performed. The mean RULA score of the 14 manikins was calculated as a mean value for all the manikins using all welding spots. The risk of WMSDs is therefore calculated by a single objective: The cycle time was calculated as the sum of the welding time at each welding spot and the time to change welding gun, welding spot, and welding side: operations (e.g., changing the welding gun, changing the welding side, and moving between welding spots). We used the Rapid Upper Limb Assessment (RULA) method [23] to evaluate the risk of WMSDs. We considered RULA appropriate since the work stresses mainly involve postural stresses on the upper limbs, as the weight of the welding guns is supported by a lifting device. For each posture assessed, RULA gives a risk score from 1 to 7, which results in four action levels. A score of 1 to 2 is acceptable, a score of 3 to 4 suggests changes may be required, 5 to 6 indicates changes will soon be required, and 7 indicates that changes are required immediately.

Mathematical Modeling of Optimization
The optimization model considers both worker well-being and productivity objectives, that is, it is a multi-objective optimization model. The indices, parameters, variables, and objectives of the optimization model are shown in Table 2.

Objectives
RULA score for a manikin on a side with a welding gun at a welding spot Cycle time of welding process (s) Average RULA score per manikin in the welding process Average RULA score of all manikins in the welding process A multi-objective optimization of the average of the RULA scores for all the manikins and the cycle time was performed. The mean RULA score of the 14 manikins was calculated as a mean value for all the manikins using all welding spots. The risk of WMSDs is therefore calculated by a single objective: The cycle time was calculated as the sum of the welding time at each welding spot and the time to change welding gun, welding spot, and welding side:

Knowledge Interpretation
The boxplot analysis showed that manikins 13 and 14 gave higher values for ER m than the other manikins ( Figure 5). It can also be seen that the minimum values for manikins 13 and 14 were higher than the maximum values for the rest of the manikins. The values for ER 13 and ER 14 were high due to the stature of these manikins (Table 1). This meant that the cluster for (4) (diversity inclusion rules) was defined by the solutions that have low values for manikins 13 and 14. Once the clusters were defined, two InfS-R analyses were run to identify the relative importance of the selection of guns, sides, and welding sequence. The first analysis was run for the objectives of CT and = A multi-objective optimization of the average of the RULA scores and the cycle time was performed. The mean RULA score of the 14 m lated as a mean value for all the manikins using all welding spots. Th therefore calculated by a single objective: The cycle time was calculated as the sum of the welding time at and the time to change welding gun, welding spot, and welding side: . The variables with the highest relative importance were the selected gun and side in spots 1, 2, 3, and 7 (GS 1 , GS 2 , GS 3 , and GS 7 ) and the welding sequence. The second analysis was run for the objectives of ER 13 and ER 14 . The variables with the highest relative importance were also the selected gun and side in spots 1, 2, 3, and 7 (GS 1 , GS 2 , GS 3 , and GS 7 ). However, for the objectives of ER 13 and ER 14 , the welding sequence did not have high relative importance since the average RULA values of the manikins are only affected by the gun and side used in the welding.
After the InfS-R analyses, four FPM analyses were run for the clusters (1) lowest cycle time (CT), (2) lowest average RULA score for all manikins ( and objectives of the optimization mod

Objectives
Cycle time of welding pr cess (s) Average RULA score pe manikin in the welding process Average RULA score of manikins in the weldin process A multi-objective optimization of and the cycle time was performed. Th lated as a mean value for all the manik therefore calculated by a single objectiv The cycle time was calculated as and the time to change welding gun, w ), (3) balance between average RULA score for all manikins and cycle time, and (4) worker diversity inclusion rules. The rules were produced by running the FPM analyses for 3 levels of rule-interactions. The rules were filtered by selecting the ones with the highest ratio of significance/unsignificance. The resulting rules apply in clusters (1), (3) and (4) to GS 1 , GS 2 , and GS 3 , as expected from the first InfS-R study. Also, a rule was found for the welding sequence X 5 in the cluster (2) for the compromise between CT and be required, 5 to 6 indicates changes will soon be required, a are required immediately.

Mathematical Modeling of Optimization
The optimization model considers both worker well-be tives, that is, it is a multi-objective optimization model. The ind and objectives of the optimization model are shown in Table 2

Objectives
RULA sco with a wel Cycle time of welding process (s) Average RULA score per manikin in the welding process Average RULA score of all manikins in the welding process A multi-objective optimization of the average of the RULA and the cycle time was performed. The mean RULA score of lated as a mean value for all the manikins using all welding s therefore calculated by a single objective: The cycle time was calculated as the sum of the welding and the time to change welding gun, welding spot, and weldi . However, it has a high unsignificance and therefore is not very interesting since it does not distinguish the selection.
The cluster with the highest ratio in the created rules is (4), worker diversity inclusion. The rules define the welding gun and side used in spots 1, 2, and 3, as expected from the second InfS-R analysis (Figure 8). The rules GS 1 == 2, GS 2 == 2 and GS 3 == 1 have a ratio of 2.46, 3.48, and 3.62 and a significance of 93.84%, 93.84%, and 82.35%, respectively, and a combined ratio of 237.34 and significance of 75.85%. This means that these rules strongly define cluster (4) against the rest of the solutions, with an unsignificance of 0.32%. If the decision maker was to select a solution of (4) (worker diversity inclusion), then welding gun 2 should be used in spots 1 and 2, and welding gun 1 should be used in spot 3. This would mean changing the welding gun at least twice in the welding sequence, with a consequent increase in CT.
The same rules were applied in cluster (2) to find the lowest average RULA score for all manikins ( operations (e.g., changing the welding gun, changing the welding side, and moving between welding spots). We used the Rapid Upper Limb Assessment (RULA) method [23] to evaluate the risk of WMSDs. We considered RULA appropriate since the work stresses mainly involve postural stresses on the upper limbs, as the weight of the welding guns is supported by a lifting device. For each posture assessed, RULA gives a risk score from 1 to 7, which results in four action levels. A score of 1 to 2 is acceptable, a score of 3 to 4 suggests changes may be required, 5 to 6 indicates changes will soon be required, and 7 indicates that changes are required immediately.

Mathematical Modeling of Optimization
The optimization model considers both worker well-being and productivity objectives, that is, it is a multi-objective optimization model. The indices, parameters, variables, and objectives of the optimization model are shown in Table 2. Table 2. Indices, parameters, variables, and objectives of the optimization model.

Welding guns
Time to change welding gun (s) = 1 …

Welding sides
Time to change welding side (s) = 1 …

Manikins
Time to move to a far position (s) = 1 …

Welding sequence
Time to move to a near position (s) Variables

Objectives
RULA score for a manikin on a side with a welding gun at a welding spot Cycle time of welding process (s) Average RULA score per manikin in the welding process Average RULA score of all manikins in the welding process A multi-objective optimization of the average of the RULA scores for all the manikins and the cycle time was performed. The mean RULA score of the 14 manikins was calculated as a mean value for all the manikins using all welding spots. The risk of WMSDs is therefore calculated by a single objective: The cycle time was calculated as the sum of the welding time at each welding spot and the time to change welding gun, welding spot, and welding side: ). While in cluster (2) the ratio for the rules combined was lower than for cluster (4), in this case the significance of each of the three rules was 100%. This means that all the solutions that have the lowest average RULA score for all manikins require using welding gun 2 in spots 1 and 2 and welding gun 1 in spot 3. This is also reflected in the lowest Sustainability 2022, 14, x FOR PEER REVIEW 4 of 15 operations (e.g., changing the welding gun, changing the welding side, and moving between welding spots). We used the Rapid Upper Limb Assessment (RULA) method [23] to evaluate the risk of WMSDs. We considered RULA appropriate since the work stresses mainly involve postural stresses on the upper limbs, as the weight of the welding guns is supported by a lifting device. For each posture assessed, RULA gives a risk score from 1 to 7, which results in four action levels. A score of 1 to 2 is acceptable, a score of 3 to 4 suggests changes may be required, 5 to 6 indicates changes will soon be required, and 7 indicates that changes are required immediately.

Mathematical Modeling of Optimization
The optimization model considers both worker well-being and productivity objectives, that is, it is a multi-objective optimization model. The indices, parameters, variables, and objectives of the optimization model are shown in Table 2.

Objectives
RULA score for a manikin on a side with a welding gun at a welding spot Cycle time of welding process (s) Average RULA score per manikin in the welding process Average RULA score of all manikins in the welding process A multi-objective optimization of the average of the RULA scores for all the manikins and the cycle time was performed. The mean RULA score of the 14 manikins was calculated as a mean value for all the manikins using all welding spots. The risk of WMSDs is therefore calculated by a single objective: The cycle time was calculated as the sum of the welding time at each welding spot and the time to change welding gun, welding spot, and welding side: solution of the Pareto front (Table 4). For the cluster (1) (lowest cycle time-CT), the rules generated are GS 1 == 4, GS 2 > 2, and GS 3 == 3, which have a significance of 91.53%, 100%, and 79.1% and a ratio of 3.38, 2.45, and 2.74, respectively. With a significance of 100%, rule GS 2 > 2 indicates that welding gun 3 (on any side) must be used to obtain low cycle time scores at welding spot 2. This collides with the rules obtained for clusters (2) and (4), which show that for a low average RULA score for all manikins ( Sustainability 2022, 14, x FOR PEER REVIEW operations (e.g., changing the welding gun, changing the welding side, and m tween welding spots).
We used the Rapid Upper Limb Assessment (RULA) method [23] to evalua of WMSDs. We considered RULA appropriate since the work stresses mainly in tural stresses on the upper limbs, as the weight of the welding guns is suppo lifting device. For each posture assessed, RULA gives a risk score from 1 to 7, wh in four action levels. A score of 1 to 2 is acceptable, a score of 3 to 4 suggests cha be required, 5 to 6 indicates changes will soon be required, and 7 indicates tha are required immediately.

Mathematical Modeling of Optimization
The optimization model considers both worker well-being and productiv tives, that is, it is a multi-objective optimization model. The indices, parameters, and objectives of the optimization model are shown in Table 2.

Objectives
RULA score for a manikin with a welding gun at a wel Cycle time of welding process (s) Average RULA score per manikin in the welding process Average RULA score of all manikins in the welding process A multi-objective optimization of the average of the RULA scores for all the and the cycle time was performed. The mean RULA score of the 14 manikins w lated as a mean value for all the manikins using all welding spots. The risk of W therefore calculated by a single objective: The cycle time was calculated as the sum of the welding time at each we ), it is necessary to use welding gun 1. Also, rules GS 1 == 4 and GS 3 == 3 for a low cycle time indicates that welding gun 3 should also be used in spots 1 and 3. The lowest CT solution in the Pareto front (Table 4) also shows that using welding gun 3 in all welding spots gives the lowest cycle time. This is due to the time saved by not changing the welding gun during the process.
Of the four clusters, the rules for cluster (3) have the lowest ratio. These rules have a high unsignificance due to the difficulty in finding rules that describe this cluster. Rule GS 1 < 4 only discards using welding gun 3 on side 2 for spot 1, and rule GS 3 > 2 only discards using welding gun 1 in spot 3. In the case of rule X 5 > 2, the rule only defines that spot 5 should be welded after two other spots. These three rules have a ratio of 1.03, 1.41, and 1.18 respectively, and a ratio of 3.02 when combined. With these low ratios, the rules do not describe the cluster (3) and distinguish it from the other solutions. This could be because defining the balance between CT and Welding gun used at each welding spot Welding side at each welding spot Objectives R w Cycle time of welding process (s) Average RULA score per manikin in the welding process Average RULA score of all manikins in the welding process A multi-objective optimization of the average of and the cycle time was performed. The mean RULA lated as a mean value for all the manikins using all w therefore calculated by a single objective: The cycle time was calculated as the sum of the and the time to change welding gun, welding spot, a cannot be done by straightforward descriptions but requires a mix of the rules of clusters (1) and (2).

Discussion
In this study we aimed to use data mining in multi-objective optimizations of worker well-being and productivity to discover knowledge that could be applied in future workstation designs. The initial optimization showed that the consideration of cycle time together with RULA scores allowed analysis of the impact of different configurations of the welding sequence. Also, optimization allowed consideration of the anthropometric diversity of the workers, helping workstation designers accommodate the diversity of the workforce. The use of different analysis and data mining methods on the database generated from the multi-objective optimization of worker well-being and productivity allowed rules to be discovered, and showed the relative importance of the design variables of the workstation.
The boxplot analysis of the individual manikins average RULA scores (ER m ) showed that manikins 13 and 14, the manikins representing the lowest stature percentiles of female and male populations, had the highest RULA scores (ER 13 and ER 14 ). The workstation design is clearly more adapted to high stature percentile populations than to low stature percentile populations. Redesigning the workstation to be at a lower height could benefit shorter workers, improving their posture while welding. Due to this result for manikins 13 and 14, we defined a cluster for the best solutions for manikins 13 and 14.
We performed two InfS-R analyses, one for the CT and in four action levels. A score of 1 to 2 is acc be required, 5 to 6 indicates changes will are required immediately.

Mathematical Modeling of Optimiza
The optimization model considers b tives, that is, it is a multi-objective optimiz and objectives of the optimization model a

Objectives
Cycle time of welding process (s) Average RULA score per manikin in the welding process Average RULA score of all manikins in the welding process A multi-objective optimization of the and the cycle time was performed. The m lated as a mean value for all the manikins therefore calculated by a single objective:

=
The cycle time was calculated as the and the time to change welding gun, weld objectives, and a second one for ER 13 and ER 14 to consider the critical manikins. In both cases, GS 1 , GS 2 , GS 3 and GS 7 had the highest relative importance. We realized that the lack of rules that included GS 4 , GS 5 , and GS 6 implied that those spots did not allow the use of welding guns 1 and 2 (Table 5). Therefore, it was impossible to improve any objective by changing only the side where welding gun 3 was used, leading to a low diversity of solutions. Removing some of the constraints that impede using welding guns 1 and 2 in spots 4, 5, and 6 would increase the diversity in the solutions, allowing new solutions to be found that could benefit both worker well-being and productivity.
After using InfS-R, FPM was used in the four clusters to find rules. The rules found apply to the actual design of the workstation, that is, the workstation that is constrained in spots 4, 5 and 6. The rules, therefore, apply mostly in the selection of gun and side at spots 1, 2, and 3. There are no rules that include the selection of gun and side for welding spot 7 since welding gun 3 obtained both better CT and Sustainability 2022, 14, x FOR PEER REVIEW operations (e.g., changing the welding gun, chang tween welding spots).
We used the Rapid Upper Limb Assessment (R of WMSDs. We considered RULA appropriate since tural stresses on the upper limbs, as the weight o lifting device. For each posture assessed, RULA giv in four action levels. A score of 1 to 2 is acceptable, be required, 5 to 6 indicates changes will soon be are required immediately.

Mathematical Modeling of Optimization
The optimization model considers both work tives, that is, it is a multi-objective optimization mo and objectives of the optimization model are show

Objectives
Cycle time of welding process (s) Average RULA score per manikin in the welding process Average RULA score of all manikins in the welding process A multi-objective optimization of the average and the cycle time was performed. The mean RUL lated as a mean value for all the manikins using al therefore calculated by a single objective: The cycle time was calculated as the sum of t and the time to change welding gun, welding spot results at that spot. The rules for clusters (3) and (4) contradict the rules of cluster (1). This means that when it comes to the actual design of the workstation, there is a clash between productive solutions (solutions with low cycle time) and solutions relating to worker well-being and inclusivity (average RULA score of the manikins). This clash does not imply that one solution should be chosen without regard for the other (equally valid) objective. Instead, the design of the workstation should be modified. In this case, welding spots 1, 2, and 3 appeared in the rules of all clusters creating the conflict between objectives. Therefore, modifying the design so that other design solutions are generated could benefit both objectives at the same time.
The knowledge discovered in this study could be applied in future designs of workstations, where workstation designers should try to not constrain welding positions for different guns. In order to keep the cycle time to a minimum, the welding gun should be changed as few times as possible. To keep low RULA average scores, the workers should be able to use the welding gun in different positions depending on their anthropometric measures. Therefore, it is critical to design the workstation for specific welding guns taking into account the available welding positions. This would allow optimizations with higher diversity in the solution space, increasing the number of design solutions of the workstation, which would lead to improved knowledge discovery. The data mining method used in this article has previously been used in other fields for knowledge discovery [18,19]. With this study we aim to include knowledge discovery methods in context of workstation design, and by that supporting engineers to find successful workstation designs and getting better understanding of what makes certain design alternatives better than others, especially in regards to worker well-being and consideration of workforce diversity.

Conclusions
Using optimization algorithms to find optimized workstation designs allows the solution space to be explored by a strategic search through feasible solutions without manually processing each of all possible configurations. The use of the InfS-R and FPM methods in the database of all solutions generated in the optimization provides a deeper understanding of the behavior of the workstation design, allowing engineers to identify critical factors in the workstations and later improve them. By using these methods in the welding gun use case, we were able to identify the critical design improvements necessary to improve both workers' well-being and productivity. We discovered that the workstation design constrained welding guns 1 and 2, removing the constraints could provide better solutions for welding spots 4, 5, and 6. This could lead to better solutions for CT and and objectives of The cycle tim and the time to ch , and also considers anthropometric diversity in the workstation. The knowledge discovered in this study could be applied in the design of future workstations, so that engineers can avoid constrained positions of welding guns and generate better workstation design solutions to improve productivity and worker well-being in factories.
The use of knowledge discovery for multi-objective optimizations of worker wellbeing and productivity can be used in relation to different workstation designs. Such knowledge can help to engineers find good design solutions for all future workstations, including other types of workstations such as assembly-line workstations.
Knowledge discovery requires engineers to have some expertise in performing multiobjective optimizations and extracting the knowledge from the databases of the optimizations. In order to further support engineers, the optimization setup and knowledge discovery process should be implemented in a digital tool that considers users' expertise and guides the users through the entire process.