C-ALB (Collaborative Assembly Line Balancing): a new approach in cobot solutions

Collaborative robots can be a proper solution to improve the throughput of manual systems without reducing their flexibility. To effectively use cobots in productive systems, it is fundamental to develop a suitable task allocation model that considers collaboration. Hence, we present a model for collaborative assembly line balancing (C-ALB) which considers paralleling tasks and collaboration in the balancing resolution. Indexes that take into account both the product and process characteristics are defined to evaluate the quality of the proposed task allocation model and comparing it to others. The results confirm the influence of the product characteristics on the system performance, leading to the definition of a new paradigm for product design.

these systems is limited, with accuracy and repeatability needing further improvements. Moreover, the productivity of these systems is influenced by ergonomic problems [5].
On the other hand, automated systems [6] present higher throughput and better quality of the final product, reducing at the same time the cost of labor. Among the several improvements of flexible assembly systems (FAS) [7,8], collaborative assembly systems (CAS) [9] are one of the solutions that could be implemented to achieve the requested flexibility. As stated by Takata et al. [10], collaborative robots, or cobots, may help to implement a dynamic productive cell, capable of sustaining a multi-model production, with the greater potential to adapt to model and volume changes.
To achieve the maximum performance from these systems, and therefore develop a cost-effective solution, it is important to correctly assign the tasks between the resources. Indeed, the idle times of a balanced assembly system are as low as possible, improving the collaboration between cobots and human operators. The activity of assigning the tasks to the stations/resources while optimizing some criterion and avoiding the violation of a certain number of constraints is called assembly line balancing (ALB) problem [11]. ALB problems have been an active area of research for the last five decades. Hence, several techniques have been proposed for their solutions, and different elements that characterize the ALB have been considered, such as layout, objective type, number of models, stochastic or deterministic processing times, and buffer presence [40]. Among the different objective types, ALB problems could be grouped into different aims: to minimize the number of stations for a given cycle time, to minimize the cycle time for a given number of stations, and also to maximize the efficiency by shifting the tasks across stations. In the literature, the ALB problem is one of the most-studied topics for assembly systems, due to its practical implications, and many researchers presented effective solutions for exactly solving ALB methods [12].
However, traditional ALB models are not designed to take into account the possibility of two resources working in the same workstation. Despite previous ALB solutions have introduced a certain degree of parallelism [13], i.e., tasks distributed in parallel between the resources to improve the flexibility of the system, these methods aim to minimize the costs of the assembly line [14] since parallel stations require doubling the equipment and fixtures. Furthermore, parallel stations are used to achieve cycle times lower than the task duration, which is not the aim of CAS. Indeed, cobots are deployed for various reasons [15] and the resources are not necessarily executing the same tasks in parallel, but usually their tasks are different.
Hence, this work aims to develop a mathematical ALB model designed for collaborative systems, hereafter called C-ALB, which considers the typical constraints of HRC: -The task times may be different for each resource in the collaborative workcell; -Each resource may not be able to execute all tasks (technological constraints); -Precedence graph permitting, the resources can complete tasks in parallel.
The paper is organized as follows: Section 2 presents the state of the art for ALB problems applied to collaborative systems and Section 3 introduces the proposed model. Section 4 introduces the evaluation indexes used to estimate the quality of the obtained task allocation, whereas in Section 6 the model is compared with a traditional ALBP model. The results are presented in Section 5 and a practical application of the model is presented in Section 7. Lastly, Section 8 concludes the work.

Literature review
As previously described, mass customization leads to the necessity for performing, and expensive, automatic systems, and hence the high practical implication of ALB problems. Since its first formalization by Salveson [18], many contributions are now available and several authors proposed different variations on the problem formulation, typically focusing on one or more specific features of the analyzed assembly system [19].
However, when considering the ALB problem applied to collaborative systems, the literature is limited. Table 1 reports a schematic representation of contributions on ALB applied to collaborative systems, focusing on works from the last decade, given the recent spread of this technology. The table groups previous researches on the basis of the objective of the task allocation, the input data required, and the proposed method. Regarding the latter, the majority of the considered works proposes a resolution by the development of frameworks for planning the collaborative process; despite this approach providing for easy-to-implement solutions, this advantage is counterbalanced by the qualitative nature of the framework. Even though the common objective of this work could be considered as improving the throughput of the system, our review identified further objectives, categorized in four main groups. In particular, Michalos et al. [20] proposed a framework with a multicriteria approach, aiming to solve the problem of task assignment considering also the need to design an optimal layout for human-robot collaboration, and comparing the different process planning solutions in terms of productivity, ergonomics, process quality, and layout efficiency.
Regarding the different objectives, the interaction between human operators and cobots has always raised safety concerns; therefore, it is reasonable that a high number of papers aims to improve operator safety during collaboration, while also considering the task ergonomics. These aspects are linked with the throughput in [26], where the authors proposed a planning method that regards both production time and physical strain. The authors introduced an objective function where a weight coefficient is defined, allowing the trade-off between these two parameters. This trade-off is also shown by Heydaryan et al. [34], who states that, despite a possible decrease in the throughput, collaboration can improve human ergonomics and reduce the risk of injury. Moreover, the framework proposed by Malik et al. [32] evaluates the complexity of assembly processes, also considering the assembly components characteristics, e.g., geometrical and physical properties, and feeding mechanisms alongside the safety issues of the tasks. The framework aims to help the designer to define the most appropriate task assignment in a human-robot collaborative system.
On the other hand, the number of mathematical models adopted is minimum and a majority of them are based on heuristic algorithms. Bogner et al. [24] imply that heuristics allows finding a reasonable solution in a short amount of time. The authors propose a matheuristic algorithm to obtain a trade-off between the goodness of the solution and the computational time. A similar approach is adopted by Weckenborg et al. [30], who deal with the problem of ergonomics and especially human fatigue issues by evaluating if the workload on the human operator is over a certain limit. However, despite presenting mathematical models based on integer linear programming (ILP) or mixed-integer linear programming (MILP), these works are focused on the comparison between the heuristic and mathematical approach, without presenting the influence of the product and process on the system performance, as seen in the considered input data, which in the majority of the studies is focused only on task data. Lastly, a fundamental aspect of HRC is indeed the collaboration between the resources. As identified by [25], cobots can share their tasks both in time and in space with the operators. Despite the lack of a quantitative parameter used to evaluate the collaboration between the resources, this capability has been considered in the literature. The authors in [27,31,35] considered different possibilities for task execution, i.e., human or cobot, or both at the same time. Similarly, [36] considered three levels of interaction: shared workspace without shared task, shared workspace and shared task without physical interaction, and shared workspace with shared task, i.e., the handing-over of parts. The authors identified that an increase in the interaction level leads to a decrease to the operator assembly time, and therefore the effects of the collaboration.
From the literature review, it is possible to observe how ALB problem resolution is one of the most important challenges in human-robot collaboration systems design. As we have tried to highlight, several approaches have been proposed, with different aims and methods, though a majority of the previous works focus on frameworks as resolution approaches and considering only the process characteristics, without considering the product influence. The lack of quantitative models, which address both product characteristics such as the precedence graph and process characteristics such as the different task times between the resources, leads to the aim of this study to develop a mathematical model based on ILP with the objective to increase not only the throughput but also the shared time between the resources.

Task allocation model for CAS
The objective of this work is to define an optimal task allocation between the resources in a CAS context. In greater detail, the addressed problem is: -To define a task allocation that minimizes the makespan ms, since the evaluation of the cycle time for a single station layout is meaningless; -To identify the effects of the product and process characteristics on the task allocation in a collaborative scenario.
For the product characteristics, i.e., the characteristics of the product that affect the performance and collaboration of the system, we mainly considered the precedence graph; indeed, it provides for information regarding the number of tasks and their precedence relations. Moreover, the task time depends on the characteristics of the product. Regarding the process characteristics, we identified the variables representing the effects due to the layout considered. Hence, we again considered the task time, since it also depends on the resource type (cobot or human operator), the task allocation considered, and therefore the achievable collaboration between the resources. As stated before, this work aims to solve the linear optimization problem by minimizing the objective function ms. To reach this goal, an optimization model based on linear programming is proposed considering the typical constraints of a cobot system intending to minimize not only the total makespan for the completion of all tasks, but also constraints related to the product and process characteristics of the system. The proposed C-ALB model is designed to optimize the collaboration between two resources, even if they are of the same type, i.e., human-human or robot-robot assembly systems. However, this work will focus only on human-robot collaborative assembly systems.

Hypotheses
The proposed model starts from the classic resolution of the ALB problem. Therefore, some of the hypotheses that characterize the model are retained in the proposed work, which are: -Mass production of one homogeneous product by performing J operations of a given product process (single-model line hypothesis); -Deterministic and integral operation time; -Each task is performed by only one resource.
As stated before, it is not optimal to apply the traditional model to a collaborative solution; therefore, other hypotheses should be defined to extend the resolution method. The hypotheses that characterize the model are: -The number of resources K is equal to 2; -The assembly line is composed by a single station consisting of one human operator and one cobot; -Collaborative resources that share workplace and task time.

C-ALB for collaborative assembly systems
We now present our model which aims to minimize the makespan ms in collaborative assembly station.
Subject to: where the optimization variable is represented by the binary variable x jkt : x jkt = 1 if task j is performed by the resource k at t 0 otherwise ∀j, k, t (9) and which by definition observes the integrality constraints in Eq. (8). Similarly to previous ALB problem resolution methods, x jkt is the task assignation variable; however, differently from previous formulations, our work introduces a third dimension (i.e., time t) to correctly evaluate ms (1). Moreover, the dimension k represents both the number of resource/station and the resource type, characterized by certain technological constraints and task times. Equation (2) represents the occurrence constraints, which, alongside the integrality constraints, ensure that in each temporal instant each task is performed by exactly one resource k. As for all the other constraints, the introduction of the temporal dimension requires modifying the classic constraint as seen in classical resolutions such as the Patterson and Albracht method. Similarly, Eq. (3) extends previous formulations of the precedence constraints. Indeed, given the definition of x jkt in Eq. (9) and the first set of constraints (2), the second member of Eq. (3) identifies the exact time instant when task j is performed by the resource k. Hence, the constraints ensure that task j starts after the execution of its predecessors i since its starting time should be greater than or equal to the corresponding ending of the execution of tasks i.
Equations (4-6) are a further set of precedence constraints newly introduced in this work and defined with a focus on collaborative systems. Indeed, Eq. (3) ensures that the precedence constraint is fulfilled only in the case of direct precedence between tasks i and j . Considering the precedence graph shown in Fig. 1, for a given task j the other J − 1 tasks can be divided into the following groups: -P j direct and transitive predecessors of task j , red in Fig. 1; -S j direct and transitive successors of task j , green in Fig. 1; -Tasks which are not direct or transitive predecessors/successors of task j .
Therefore, using only the constraints in Eqs. (2-3) does not consider the third group and tasks that belong to it could be performed at the same time by a resource k. To avoid this scenario, it is fundamental to ensure that the start time of task j is set after the end time of task i, that is: where start j represents the start time of tasks j , and end i represents the end time of task i and which corresponds to 2 sets of constraints, i.e., for j greater than i or j lesser than i. To correctly define these constraints, it is necessary to define the binary decision variable y ji such as: Hence, the constraints in Eq. (10) could be rewritten as: Hence, if task j is executed before task i, y ji equals to 1 and Eq. (12) is always true and Eq. (13) ensures that task i starts after the complete execution of task j . On the contrary, if task j is executed after the end of task i, Eq. (12) ensures that task j starts after the complete execution of task i, whereas Eq. (13) is always true. Lastly, Eq. (7) represents the technological constraints, i.e., it limits the possible solutions by considering that some Fig. 1 Identification of the P j predecessors and S j successors for a task j considering a precedence graph tasks j could not be carried out by the resource k. This is represented by the set U k , whose elements u jk are the tasks j that cannot be executed by resource k, differentiating between different types of resources. Indeed, by considering different task times and capabilities for the different types of resources, it implies that the model can work on layout characterized by not homogeneous resources, i.e., humanrobot collaboration. However, by considering comparable task times and capabilities, it is possible to apply the model to homogeneous scenarios.

Product and process characteristics indexes
To better consider the influence of the product and process characteristics on the system performance for a considered task allocation, four indexes have been defined. These values take into account both input and output variables and both product and process characteristics and were divided as such: -Parallelism index p % : input parameter based on the precedence graph (product characteristic); -Task time index t % : input parameter based on the difference of the resources task times (product and process characteristic); -Makespan index m % : output parameter which evaluates the quality of the ms achievable with the considered task allocation (process characteristic); -Collaboration parameter c % : output parameter that represents the shared activity between the resources [9] (process characteristic).
It should be noted that these parameters do not define the model; thus, they are not considered in the model description. This section aims to presents these four indexes, defining them and explaining briefly their influence on the system.

Parallelism index p %
The parallelism index (p % ) is used to evaluate the number of parallel branches in a precedence graph, thus characterizing the product influence on the CAS. This is evaluated as the arithmetic mean of the ratio between the number of arcs for each task, hence: where d j represents the number of arcs for task j and can be evaluated as: (a) p % = 0 (b) p % = 1 Fig. 2 Precedence graphs corresponding to the extreme cases for the parallelism index p % and corresponds to the number of tasks that cannot be executed in parallel with task j . The proposed formulation allows identifying the parallelism degree for a precedence graph between two boundaries: -If p % is null, the tasks must be carried out in a sequential order as seen in Fig. 2a; indeed: -If p % is 1, the tasks can be carried out independently as seen in Fig. 2b; thus:

Task time index t %
The task time index (t % ) is based on the model input data, i.e., the disparity between the resources task time, which depends on the process adopted to assembly the considered product. Since it is possible to observe an increase in the collaboration between the resources when their task times are comparable, the necessity to study this parameter is reasonable. t % is evaluated as the ratio between the minimum makespan, evaluated by assigning all tasks to the respective fastest resource, and the maximum one, evaluated by assigning all tasks to the respective slowest resource: This formulation takes into account a single parameter for both resources since it is independent of the k index.
Similarly to the p % index, the value of t % is limited between 0 and 1 where: -If t % is very small, albeit not null, the disparity between the resources task time is high. A null value cannot be reached unless the meaningless scenario where one of the resources has null task times; -If t % is equal to 1, the resources task times are equal on average; thus, this scenario is likely to increase the collaboration between the resources.
Indeed, since the number of tasks is the same for the resources, comparing the sum of the task times or the mean task time is equal. This is shown in Fig. 3, where considering a faster operator on average (dotted line in red, 5.42 s for the operator, 6.75 s for the cobot) does imply a smaller sum (65 s for the operator and 81 s for the robot) but not that all the task times are shorter, such as tasks 9 and 10.

Makespan index m %
To evaluate the effectiveness of the scheduling proposed by the model, it is important to compare the obtained ms with  Fig. 4 Task allocation for p % = 0, whose makespan is ms p % =0 a reference parameter. Hence, the makespan index m % is evaluated as the ratio between the estimated makespan ms and the achievable makespan for null p % : where ms p % =0 is the achievable makespan for a null parallelism index, i.e., the makespan evaluated for a precedence graph as seen in Fig. 2a. The minimization of this value can be solved by assigning each task to the fastest resource, as shown in Fig. 4. This formulation allows to compare the makespan obtained with different task allocation methods with a reference one; moreover, it is independent of the task times and J . Lastly, a suboptimal solution could lead to a value greater than 1, since this index is not limited like the previous ones.

Collaboration parameter c %
The collaboration parameter (c % ) represents the interference between the human and the robot during the assembly process. Previous studies presented different formulation for c % and this work adopts the one presented in [39]: where T coll is the shared time between the resources, as seen in Fig. 5.

Influence of the product and process characteristics on the CAS system
The following paragraph presents the influence of the product and process characteristics, i.e., p % and t % , on the quality of the task allocation for the CAS system, whose results will prove is not negligible. As previously described, p % depends on the product design; thus, it is a product characteristic, whereas t % represents both the product and process characteristics. As an example, we consider a screwing process: the difference in the task times depends on both product characteristics such as screw size, part positioning, and process characteristics, e.g., presence of fixtures, force sensors for the cobots, tools. Moreover, the influence of t % and especially of p % is a peculiarity of CAS, which is due to its definition of collaborative resources. Lastly, the technological constraints in Eq. (7) were not considered to avoid adding further parameters in the influence analysis.

Influence of p % on m % and c %
The first tests aim to identify the influence of the precedence graph, represented by p % , on the system performance, i.e., m % and c % . The tests were carried out with t % constant and varying J ; for each J , 10 values of p % were identified as seen in Table 2.
Regarding the other input value, t % , the used task times were considered between 2 and 10 s, and with a constant t % equal to 97%. The obtained results show that the influence of J on the system performance is less significant than the influence of p % ; therefore, it is not considered in the following graph. The results were fitted with different models, such as hyperbolic, power, and through the coefficient of determination R 2 , a 3rd degree polynomial form was considered the most suitable one, with an R 2 value of at least 0.95. The obtained non-linear equations suggest that the precedence graph could present a modest value of p % to greatly improve the performance of the system. Figure 6 represents the effect of p % on the makespan, represented by the evaluation index m % . As defined, for low values of p % , the makespan is comparable with ms p % =0 ;  thus, m % tends towards 1. Moreover, for a higher value of p % , m % presents an asymptotic behavior, identified by the red dashed line, with a limit equal to 1 K , which corresponds to a value of 0.5 for the considered input values. When m % reaches the asymptotic value, the performance of the system is at their utmost limit and the collaboration parameter c % is equal to 1. Figure 7 shows the effect of p % on the collaboration parameter c % . A null collaboration value is correctly expected for a null value of p % , since a similar precedence graph does not allow shared tasks. On the other end, for a high value of p % , c % tends towards the maximum value.
The obtained results prove that a higher value of p % improves the performance of a CAS since it is possible to execute more tasks in parallel, thus increasing the shared time. Moreover, it proves the need to study the influence of the precedence graph when solving the task allocation problem for a collaborative solution to improve the performance of the system, as seen in Fig. 8, where the relation between c % and m % is represented. Indeed, the obtained curve proves that an increase in collaboration leads to a lower makespan, therefore proving the necessity of an ALB methodology focused on collaborative application to improve the collaboration, with the improvements shown in

Influence of t % on m % and c %
Differently from p % , it is not possible to distinguish the effect of J and t % ; therefore, the obtained graphs are presented for different values of J . Indeed, for a certain value of t % , a lower value of J means that the difference between the tasks times is less distributed; thus, the possibility that the resources are not synchronized is higher, leading to a higher makespan. Therefore, not considering J would lead to an unreliable analysis. Furthermore, the simulation tests showed the influence of p % on the effects of t % , even if the overall behavior does not change. Hence, the tests were carried out for different values of p % , represented in Table 3, where for each value of J , 3 values of p % were considered to represent low, medium, and high degrees of parallelism. For each J , 21 values for t % defined between 0.2 and 1 were considered; these are omitted for shortness.
Similar to the previous tests, the results were fitted adopting a 3rd degree polynomial form, since it resulted in a higher value of R 2 . Figure 9 represents the influence of t % for low values of p % , i.e., lower than 0.27. As previously seen with the p % -tests, the non-linear behavior shows that modest changes in t % greatly improve the system performance; however, in this case for moderate task time  difference, i.e., t % greater than 0.6, m % is constant, and thus t % does not influence the makespan. The effects of J depend on the adopted value of p % , as seen in Fig. 10: for p % about 0.49, the effect of J is limited with respect to the behavior shown for p % about 0.27. This shows that for low values of parallelism a higher number of tasks J improves the system performance, represented by lower values of m % . Indeed, even if p % is low, a higher J allows carrying out a higher number of tasks, thus reducing the makespan, further confirming the importance of considering parallel tasks when defining the task allocation for collaborative workcells. This is further proved by Fig. 11, where the influence of J is minimal for a moderate task time difference (t % greater than 0.6). Moreover, an increased p % leads to a lower value of m % in the plateau region, further proving the impact of p % on the system performance. Lastly, the results did not show a great influence between c % and t % .

Overall results
The overall behavior is presented in Fig. 12, where the influence of both p % and t % on m % is presented, considering a product with J = 10 as an example. The overall behavior is represented by adopting a linear interpolation between the scenarios; although it is not the most faithful fit, it immediately shows the effect of the product and process characteristics on the system performance, which was represented more faithfully in the previous figures.
The obtained results show that it is necessary to take into account the difference of the task times and the precedence graph parallelism to improve the system performance. Moreover, it further proves that, although the graph shows a steeper slope along the direction of t % than that of p % , which hints to a greater influence of t % , a higher value of p % allows achieving better improvements, especially in the case of a limited number of tasks, as represented in the proposed scenario.

ALB problem model comparison: traditional versus C-ALB
In this section, the proposed model is compared with a variant of the ALB's resolution method studied by Patterson and Albracht [40] by analyzing their resulting task allocation and assembly system performance. Indeed, the proposed task allocation model is needed since the traditional ALB resolution, which minimizes the cycle time, does not correctly represent a collaborative application, since the modeled layout is a simple assembly line and the constraints do not consider the parallel branches in the precedence graph. Therefore, it is understandable to presume that a task allocation obtained with this model is not optimal to both productivity and collaboration.
To observe the differences between the two models, several precedence graphs were considered to obtain the values presented in Table 4 for the input variables J , p % , and t % . Therefore, the comparison tests were conducted by combining the input data; these values were chosen to point out how the presence of parallel branches and different task times between the resources influence the need for a proper collaborative model.
The task scheduling resulting from the simulations confirms the assumption made before. For example in Table 5, the results obtained with J = 15, p % = 0.27, t % = 0.4, are presented: both collaboration time and makespan are smaller when the model used is the one studied for collaborative systems. Moreover, the goodness of the allocation can also be evaluated starting from the makespan index, which is always less than 1 with the model presented in Section 3, while this is not true for the traditional formulation.  The obtained results show that the traditional approach, i.e., Patterson-Albracht, that we considered as an example due to its spread, is not suitable to solve the task allocation problem when considering an HRC application. Indeed, the comparison shows a reduction of about 31% when adopting the C-ALB proposed model, and an increase of the shared task time of 520%. This last result shows that adopting traditional approaches that do not consider the collaboration between the resources, i.e., parallel task execution, cannot be used for CAS, further proving the necessity of an ALB model designed for CAS. To reach this goal, as previously described, it is fundamental to develop a model that takes into account product characteristics such as p % , which are unique to collaborative systems.

Case study
The following section presents an industrial application of the proposed model to demonstrate its practical implication; Fig. 13 Assembly case of a box Place the support 4 Place the support 5 Place the support 6 Place the outline border 7 Place the outline border 8 Place the outline border 9 Place the outline border 10 Place the cover board    1  1  3  3  2  3  6  5  3  3  6  5  4  3  6  5  5  3  6  5  6  7  14  10  7  7  14  10  8  7  14  10  9  7  14  10  10  10  22  20 the considered assembly is a simplified adaptation of a real industrial process. This simplification allowed us to consider two different values of p % without changing the tasks. A description of the task activities needed to assembly the product in Fig. 13 is presented in Table 6.
The considered product is defined by a number of tasks J equal to 10 and with the following input characteristics: p % equal to 0.13 (Fig. 14a) and 0.48 (Fig. 14b) -t % equal to 0.13 and 0.91 Regarding t % , in this work, the operator task times were taken from experimental evaluations and thus considered constant, whereas the robot ones were evaluated on the basis of the chosen value of t % and assigned to the manipulator. This is due since it was difficult to change the operator's speed, differently from the robot. The adopted task times are presented in Table 7.  16 Task allocation for the case study: the red cross indicates the task assigned to the operator, the blue ones are assigned to the cobot A collaborative workstation has been developed in the Robotics and Automation Laboratory at the Department of Management and Engineering of the University of Padova, as represented in Fig. 15. The tests were carried out using a KUKA LBR iiwa 14R820 cobot to perform the tasks assigned offline by the proposed model.
As a result, Fig. 16 presents the spatial task allocation for the case study for one of the scenarios: the points marked with a red cross are the ones assigned to the human operator, the others are carried out by the cobot, whereas Fig. 17 presents the task allocation between the human operator (in red) and the cobot (green); moreover, the blue bar presents the achieved collaboration. The case study confirms the need for certain characteristics, i.e., t % and p % , to obtain the maximum advantages from the adoption of a cobot. Indeed, for the considered scenario, i.e., t % = 0.6 and p % = 0.48, the system achieves a makespan reduction, with m % equal to 0.76, corresponding to a reduction of the makespan from 57 to 39 s, and an increase in collaboration with c % = 0.61.

Conclusion
Collaborative assembly systems are increasingly common in the industry since collaborative robots are considered an optimal solution to increase both the flexibility and throughput of the system. However, using traditional ALB models is not optimal since a CAS workstation is composed of multiple resources working in parallel in the same workstation. Hence, this work presents a model for assembly line balancing, called C-ALB, developed for collaborative systems. From the proposed model, the following results could be obtained: -A proper model for collaborative systems should be considered; indeed, when comparing with traditional ALB models, our C-ALB model achieves a greater throughput and collaboration. -A difference in the task time between the resources has a minor influence on the system performance unless it is considerable (t % less than 0.4). Moreover, the number of tasks influences the performance for a lower value of p % . -Regarding p % , precedence graphs with more parallel branches are more suitable for CAS since the proposed model shows that p % has a great influence on m % .
As previously observed, this last result shows the importance of the product design when considering a collaborative system; hence, a new paradigm called design for collaborative assembly (DFCA) must be investigated and developed in future works, but this work proves its importance. Moreover, the proposed model has been defined under the hypothesis of mass production of one homogeneous product; future works will extend the model to Mixed-Model Assembly Line Balancing (MALB) problems. This extension will allow taking advantage of the cobot ability to be flexible in processing a range of product variations. Further future works will consider CAS composed of multiple resources and workstations in the assembly line. We will especially investigate the case of big size products, where their dimensions make the use of multiple resources within the same workstation feasible [41].

Fig. 17
Task allocation for the case study (t % = 0.6, p % =0.48): the red bar indicates the task assigned to the operator, the green ones are assigned to the cobot and in blue the collaboration Author contribution The corresponding author M. Faccio has been responsible for planning and coordinating the steps of the research. R. Minto, M. Milanese, and M. Faccio have been responsible for writing and reviewing this paper, together with the quantitative model proposal and the analysis. G. Boschetti has been responsible for revising the paper, including the suitable structure and contents of this paper.
Funding Open access funding provided by Università degli Studi di Padova within the CRUI-CARE Agreement. Open Access funding provided by Università degli Studi di Padova.

Declarations
Ethics approval and consent to participate This paper is new. Neither the entire paper nor any part of its content has been published or has been accepted elsewhere. It is not being submitted to any other journal as well.

Conflict of interest
The authors declare no competing interests.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons. org/licenses/by/4.0/.