In-line kitting for part feeding of assembly lines: workload balancing and storage assignment to reduce the workers’ walking effort

An efficient part feeding is among the top challenges of many mass producers applying mixed-model assembly lines, for instance, in the automotive industry. This paper introduces a novel part feeding policy applied by a large German assembly plant for car engines: In-line kitting. Under this policy, the first stations of the line do not execute assembly operations, but are reserved for picking parts while passing containers of stock-keeping units (SKUs) arranged along the line. In this way, the parts are collected in traveling kits moving along with each workpiece on the conveyor, so that later assembly stations have the required parts directly available and do not lose precious labor time for unproductive parts handling. A major operational challenge when applying this part feeding policy is the walking effort for the human pickers while putting the SKUs of their respective stations into the traveling kits of the passing workpieces. Due to a high product variety, a large number of comparatively bulky SKU containers have to fit into each station, so that the walking distance to be covered by a worker during a work shift exceeds multiple kilometers. We show that this physical burden can be reduced significantly by balancing the workload among stations and optimizing the storage assignment of SKU containers within each in-line kitting station. We formulate the resulting optimization problem and provide suited solution procedures. Our computational study shows that the walking distance of pickers can be reduced significantly without producing any additional costs.


Introduction
The mid-sized assembly plant of German car producer BMW in Dingolfing has to timely supply their assembly stations with parts arriving in more than 13,000 containers delivered by about 600 suppliers on more than 400 trucks on each average workday (Battini et al. 2013). These figures elucidate that part logistics is among the greatest challenges of today's mass-producers, especially in the automotive industry (Boysen et al. 2015). Therefore, it is anything but surprising that selecting the right part feeding policy has received much attention both among practitioners and scientists in the recent years. Part feeding is the in-house logistics process which moves parts from receipt after their arrival at the assembly plant up to their respective workstations where they are finally assembled into the workpieces. Note that a similar definition is provided by Choi and Lee (2002).
A part feeding policy specifies this in-house logistics process for a subset of parts. To do so, a specific policy has to choose among the alternative elements of the four basic process elements defined in Fig. 1. Note that previous papers (i.e., Battini et al. 2013;Boysen et al. 2015;Kilic and Durmusoglu 2015;Schmid and Limère 2019) apply similar, yet slightly different elements of part feeding: • Initially, the part feeding process is triggered either by previous consumption or by prospective demands, so that parts are either pulled (e.g., via a Kanban mechanism or by part inventory falling below a predefined reorder level) or pushed (i.e., according to the demands defined by the workpieces produced in the next production cycles) toward their respective assembly stations. • Once the demand for parts is defined, the respective parts need to be retrieved from their storage positions. This can either be a central (receiving) warehouse, a decentralized logistics area directly on the shop floor (often called supermarket, see , or local storage space directly next to the assembly line. • Then, parts need to be moved by a suitable transportation device, which may be a traditional forklift, a tow train (i.e., a manned or automated tugger vehicle towing a handful of wagons (see , or a conveyor belt. • Finally, the parts need to be positioned next to the line, which is also called their line-side presentation (see also Limère et al. 2012; Schmid and Limère 2019). In-line kitting for part feeding of assembly lines: workload… Line stocking places containers (e.g., large unit loads or smaller bins) each filled with homogeneous pieces of the same stock-keeping unit (SKU) directly in the respective workstations (e.g., multiple bins each containing a specific exterior mirror variant), so that the assembly workers have to identify and fetch the respective parts required in each production cycle. To avoid unproductive part sorting of (highly paid) assembly workers stationary kits are placed next to the line and contain pre-sorted pieces of different SKUs ordered according prospective demand (e.g., the exterior mirror variants sequenced for the next production cycles). Traveling kits, instead, contain pre-sorted pieces of multiple SKUs for a single workpiece (e.g., the exterior and interior mirrors and other parts for a specific car) and accompany their workpiece on the conveyor.
The most widespread part feeding policies resulting from choosing a specific element for each process step are traditional bulk feeding and the supermarket concept (see Fig. 1). Triggered by previous part consumption, the former applies forklifts to deliver large homogeneous unit loads from central storage to the line and leaves the sorting of parts to the assembly workers (Boysen and Bock 2011). Decentralized storage in supermarkets allows for a more reactive part supply in smaller bins of pre-sorted products and applies tow trains delivering multiple workstations along fixed routes . The supermarket concept can be applied in a push and pull environment and can service all kinds of line-side presentation. Identifying a suitable part feeding policy for each part is an important practical decision task, and the pros and cons of each policy are vividly discussed in the scientific literature (e.g., Hanson and Brolin 2013; Sali et al. 2015). Moreover, quantitative decision models for selecting the right feeding policy (e.g., Limère et al. 2012;Caputo et al. 2015Caputo et al. , 2018 and its interdependency with assembly line balancing (Sternatz 2015;Battini et al. 2017;Nourmohammadi et al. 2019) have been investigated. The manifold literature in this area is summarized by the survey papers of Boysen et al. (2015), Kilic and Durmusoglu (2015), and Schmid and Limère (2019). This paper treats a part feeding policy, which (to the best of the authors' knowledge) has been overlooked by previous research. We elaborate this novel policy, which combines the elements marked by the dotted arrow in Fig. 1, in the following section.

In-line kitting
We saw in-line kitting at work in the assembly plant of MDC Power in Kölleda, which is one of Germany's largest assembly plants for car engines and belongs to the Daimler group. The engines are produced on multiple mixed-model assembly lines in lot-size one. Each engine moves from station to station on a hanging workpiece carrier transported by a trolley conveyor. The first stations of the line are reserved for in-line kitting, which means that the specific parts dedicated to each engine are directly added to the cases of their carriers. The carriers, thus, also serve as traveling kits. In this way, the parts required in later assembly stations are directly available on each workpiece carrier, so that unproductive walking of assembly workers for fetching parts from racks and containers next to the line is eliminated.
The in-line kitting segment of the assembly line, whose schematic layout is depicted in Fig. 2, is separated into distinct areas called stations each operated by a single dedicated worker. To better distinguish these workers solely concerned with part logistics from assembly workers of later stations, we call them the pickers. Each station consists of a given set of successively arranged skeleton containers each containing homogeneous pieces of a specific SKU. Pick-to-light displays indicate which SKUs (and how many of each) have to be added to the cases of the engine carrier currently entering a station. The picker accompanies the engine, which is continuously moved forward by the conveyor, adds the indicated parts whenever they pass the respective SKU, and accredits the picks (by pushing a button on the pick-to-light display) until the end of the station is reached. Naturally, there are no precedence constraints that restrict the sequence in which the SKUs are packed on the carriers, so that the pickers are advised to follow the sequence in which the SKU containers are arranged along the line. Then, the picker rushes backward toward the successive engine of the next cycle. Since engines are produced in lot-size one, each engine has its different part demands. Given the specific part demands of the next cycle, the picker repeats the kitting process for the next workpiece. Note that the kitting line segment itself is replenished via a bulk feeding process (see Fig. 1), where forklifts replace (almost) empty with filled skeleton containers delivered from central storage once a specific reorder level is reached.
The biggest advantage of in-line kitting is that assembly workers, who are specially trained and earn higher wages, are relieved from part fetching. The parts are directly available on the cases of each engine carrier, so that unproductive time for walking to part containers next to the line as well as identifying and retrieving the required parts from there is eliminated. These work contents are executed by less costly pickers in the kitting area, who are typically assigned to another pay scale group. Another advantage is that there is no need to store part containers in assembly stations, which is especially valuable for small sized workpieces such as car engines where stations are much smaller compared to automobile assembly. These advantages, however, could also be realized by moving traveling kits assembled in In-line kitting for part feeding of assembly lines: workload… some remote logistics area onto the conveyor. In engine assembly, however, there is no belt conveyor, so that a traveling kit would have to be attached to the carrier, which requires a suited (and stable) technical solution. Furthermore, in-line kitting moves part logistics under the strict regime of a paced flow process, so that part feeding is under tight control and picking errors can easily be reduced, e.g., by adding an extra weighing mechanism for the added parts to the workpiece carrier.
On the negative side, in-line kitting requires an extended conveyor and blocks space directly on the shop floor where room is notoriously scarce. Car engines and the part variety they require, however, are smaller compared to typical assembly line products such as cars, so that the prolongation of the line due to the additional kitting segment is moderate. Furthermore, trolley conveyors are less costly than large ground conveyors for cars, so that in the specific case of engine assembly the additional investment for in-line kitting seems acceptable. These characteristics question the suitability of in-line kitting for car assembly, but make it a valid part feeding alternative for small product assembly, e.g., in the electronics industry.

How to reduce the pickers' walking effort
For our engine producer, one of the greatest operational challenges when applying in-line kitting is the high physical stress for the picker workforce. Since the parts in engine assembly are small and light-weighed, ergonomic stress caused by lifting, which is a major cause for long-term back injuries in many industries (e.g., Grosse et al. 2015;Otto et al. 2017), is less of an issue. Instead, accompanying the workpieces while putting the indicated SKUs into the engine carriers and walking back to the start of the station again and again during a work shift accumulates to a considerable total walking distance for each picker. The station each picker services is relatively large due to the following two reasons: • The standardized skeleton containers for the SKUs applied by our engine producer are comparatively large (i.e., they have the same outline as a europallet with a width of 800 mm. Due to the small size of the parts, these containers have the advantage that plenty pieces of a SKU fit into each container, so that the logistics effort for replenishing the kitting stations is reduced. • Each assembly line produces not only a single type of engine, but multiple variants in a mixed-model setting (see, e.g., Boysen et al. 2009). Therefore, each station contains many SKUs, but only some of them are required by any type of engine. Recall that the part variety is small compared to car assembly but still considerable.
This results in station sizes of dozens of meters. At a cycle time between 60 and 90 seconds, which are typical values in the automotive industry (see Boysen et al. 2015), the walking distance quickly accumulates to several kilometers per working shift. If a worker has to service 50 SKUs whose containers have an identical width of 800 mm, the resulting area width is 40 m. Under the assumption that the picker always has to traverse the complete area at a cycle time of 60 seconds, the worst-case walking distance of a 7-hour working shift amounts to 17 km. We admit that this is rather a ballpark calculation to highlight the potential worst-case walking effort of workers. In the real-world, these extremes are seldom reached due to more generous break regulations, job rotation, or longer picking times that enforce smaller stations. However, the managers of our engine producer told us that more than 10-km walking per shift regularly occur, which they suspect as a major source for the comparatively high level of absenteeism and labor turnover in the kitting area. Furthermore, the kitting area employs many impaired workers not suitable for the even more demanding workplaces in the assembly stations, so that the managers (and the representatives of the trade union) aim to reduce the pickers' walking distances.
Obvious levers for reducing the walking distance in the kitting segment are a larger workforce or an extended cycle time. The former, however, increases wage costs and the latter reduces the output of the assembly line. Hence, the management rather prioritize less costly countermeasures, such as the following. Pickers need not traverse their complete stations in every cycle. Once pickers recognize that the remaining pick-to-light displays indicate no further demand for the current engine, they can let the workpiece pass the remainder of their station unattended and can prematurely return back toward the next workpiece. Such a reduction in the walking effort can be influenced by the following two decisions tasks: • Workload balancing: Naturally, only SKUs required by the current workpiece and located within the respective station have to be serviced by a picker. Thus, distributing the SKUs among the kitting stations in a fair manner enables a workload balancing for pickers. In this way, exceptionally large walking distances of single pickers can be avoided. Thus, the partitioning of SKUs among kitting stations influences the walking effort of the pickers, and this decision task is tackled by the optimization procedures developed in this paper. Workload balancing alone (without the subsequent storage assignment within each station, see below) resembles the famous assembly line balancing problem (see Boysen et al. 2022Boysen et al. , 2008, which decides on the division of labor among subsequent assembly stations. In our case, we face a special line balancing problem without precedence constraints (i.e., the SKUs can be processed in any sequence) and including the model mix (see Merengo et al. 1999), because each engine type requires a different subset of SKUs. The biggest difference to existing research, however, is the interdependence with the following decision task. • Storage assignment: Another lever, which reduces picker walking without deteriorating picking performance, is the storage assignment of SKUs within each kitting station. It decides on the sequence, in which the set of SKUs of a station are arranged along the conveyor. In our setting, the picking orders, i.e., the subset of SKUs demanded by each engine, are known with certainty. If we find a storage assignment, such that the spread of each picking order, i.e., its distance from the first to the last container required by the respective order, is reduced, then pickers need not traverse their complete stations, but only a smaller subsegment. For this purpose, we integrate the storage assignment of SKUs to reduce the pickers' walking effort into our optimization procedures. Storage assignment of SKUs is 1 3 In-line kitting for part feeding of assembly lines: workload… one of the classical decision problems of warehousing research (see De Koster et al. 2007). This also holds true for a special order picking system called trolley line picking (see Boysen et al. 2021). Here, hanging trolleys each assigned to a specific customer order pass by SKUs, and it is the duty of human pickers to add requested SKUs to the trolleys. Thus, trolley line picking directly transfers the basic setup of in-line kitting to the warehouse domain. Here, however, customer orders are typically not known with certainty, when having to decide on the storage assignment, we have no continuously moving line, but trolleys automatically stop in front of requested SKUs, and pickers need not strictly pick one order after the other but can switch among waiting trolleys (see . Furthermore, the main lever to improve these systems is the sequence, in which the orders are processed (see ). This lever, however, is not available to us, as is argued in the following.
The production sequence, in which the engines are sent down the line, can also relieve pickers. If the storage position of the last SKU within a picker's station demanded by the previous engine is close to the storage position of the first SKU of the subsequent engine, the walk-back distance between successive engines can be reduced. In-line kitting, however, is typically dominated by the subsequent assembly stations, which have a complex demand for specific production sequences obeying different technical restrictions and leveling constraints (see Boysen et al. 2009).
Thus, determining the model sequences according to the walking effort of kitting workers is, typically, not an option. Note that other industrial settings (e.g., where heavier or bulkier parts are to be handled) suggest a more holistic view on the physical strain of workers that also includes the ergonomic burden caused by lifting and unfavorable body postures. Blueprints on appropriate measures are provided by the order picking literature (see, e.g., Calzavara et al. 2017;Glock et al. 2019). An interesting future research issue in this direction would, for instance, be field tests that investigate for which types of products the total walking distance is a sufficient proxy for the total physical burden and for which more detailed models should be selected.

Contribution and paper structure
After the paper's first contribution, namely introducing and discussing the in-line kitting part feeding policy, the remainder of the paper is concerned with reducing the walking effort of pickers putting together the traveling kits under the in-line kitting policy. For this purpose, we formulate the joint workload balancing and storage assignment problem, which combines both levers for reducing the pickers' total walking distances in a holistic problem setting. Unfortunately, our computational results will show that a mixed-integer model for the holistic problem cannot be solved by an off-the-shelf solver for instances of real-world size. For solving these instances, we rather introduce different decomposition approaches addressing both decision parts in an iterative manner. Once a suitable heuristic solution procedure is available (and proven successful in a comprehensive performance test), we can apply this algorithm to explore to what extent the total walking distance of pickers can be reduced in our case study. Our results indicate that the walking distances of pickers can be reduced significantly without deteriorating picking performance. The remainder of the paper is structured as follows. Section 2 elaborates our optimization task and provides a mixed-integer model for the holistic problem. Two competing heuristic decomposition approaches, namely, balance first, sequence second and (vice versa) sequence first, balance second, are presented in Sect. 3. The computational performance of our solution methods is investigated in Sect. 4. This section also contains further tests on managerial issues, where we explore the possible reduction of the pickers' walking effort. Finally, Sect. 5 concludes the paper.

The joint workload balancing and storage assignment problem
This section treats the joint workload balancing and storage assignment problem, which we dub WBSAP. Section 2.1 defines the problem, and Sect. 2.2 discusses our basic assumptions and investigates its complexity status. Finally, Sect. 2.3 presents a mixed-integer programming (MIP) model for WBSAP.

Problem definition
Consider a mixed-model assembly line producing n different types of car engines. Each engine i ∈ I = {1, 2, … , n} is produced with a given frequency f i within the planning horizon and requires a specific subset of parts for assembly. We define the total SKU set of parts by S = {1, 2, … , |S|} and the subset of SKUs required by engine i ∈ I constitutes its associated picking order O i ⊂ S with O = {O 1 , O 2 , … , O n } . Each SKU is available in a separate skeleton container. All containers are standardized and have identical size; they only differ in the dedicated SKU they contain. Note that we discuss this and all further assumptions in more detail below. The storage positions P = {1, 2, … , |P| = |S|} along the line are sorted in ascending order of their index. The part demands have to be fulfilled by a kitting line segment operated by a picking workforce of given size. Each picker operates a separate station, so that we have a given set K = {1, … , |K|} of kitting stations.
Our decision task is to partition the total SKU set among the stations and to decide on the storage assignment within each station. These two decisions can be encoded by (Φ, ) , i.e., • a family of sets Φ = {Φ 1 , Φ 2 , … , Φ |K| } with Φ k defining the subset of SKUs assigned to kitting station k and • a SKU sequence , i.e., a permutation of S defining the storage assignment of each SKU along the line.
We call a solution (Φ, ) feasible, if the following three conditions are met:

3
In-line kitting for part feeding of assembly lines: workload… that is, a SKU cannot be assigned to more than one station at a time, • ⋃ k∈K Φ k = S , that is, each SKU is assigned to at least one station, and • ∃k ∈ K, p, p � ∈ P ∶ p < p � ∧ (p), (p � ) ∈ Φ k → (p + 1) ∈ Φ k , that is, only consecutive SKUs can be assigned to the same station.
Our aim is to reduce the pickers' walking effort. Among all feasible solutions, we, thus, seek those which allow to assemble each picking order in a small subarea of each station. In this way, each picker does not have to traverse her complete kitting station in each cycle. We formalize this aim by the order spread i,k (Φ, ) , which defines the number of SKU containers between the first and the last occurrence of a SKU required by order O i within the storage assignment of station k. Specifically, where function −1 (s) returns the storage position p ∈ P = S of SKU s within storage assignment .
Given the standardized container sizes we presuppose, the order spread directly defines the distance each picker has to accompany an order while picking. Unfortunately, at the point in time storage assignments are planned, the production sequence of engines is, typically, not available, so that we cannot exactly anticipate the total walking distance. To do so, the storage assignment, which defines the walking distances during order picking, and the exact order sequence, which defines the walkback distance between two successive orders for walking from the last position of the predecessor to the first of the successor, both need to be available. As the latter information is not at hand, we have to do without and aim to minimize the total weighted order spread In our case, it is well-known how many of each engine type are to be produced within the next working shifts (but not their exact production sequence), so that it seems advisable to weigh the order spreads ik (Φ, ) with picking frequencies f i . We dub the problem that seeks the minimum total weighted order spread Z(Φ, ) among all feasible storage assignments , the joint workload balancing and storage assignment problem (WBSAP).
Example: Consider the example data given in Fig. 3a. Four engines i = 1, … , 4 each produced once during the planning horizon have different parts demands for six SKUs S = {A, B, … , F} . These picking orders have to be fulfilled by two workers in kitting stations k = 1, 2 . Solution one depicted in Fig. 3b shows an objective of Z(Φ, ) = 10 , constituted by the total walking distance of Station 2 exceeding that of Station 1. Solution two (Fig. 3c) reduces the maximum total walking distance to Z(Φ, ) = 8 by altering the storage assignment within both stations, whereas the division of labor among both stations remains unaltered. Balancing the workload by reassigning SKUs among both stations leads to Solution three and a further reduction to Z(Φ, ) = 7 depicted in Fig. 3d.

Assumptions and complexity
After having defined our optimization problem, we summarize the simplifying assumptions (explicitly and implicitly) contained in our problem setting: • We aim at a reduction in the physical effort of pickers induced by their total walking distance. We consider this general aim by minimizing the maximum walking distance among all stations. This is considered to be more fair compared to a min-sum objective where the total walking distance accumulated over all stations may be smaller, but some pickers may receive exceptionally large walking distances. • Furthermore, we want to relieve pickers without causing additional costs, so that we assume a given picker workforce. In line with the setting of our engine producer, each picker operates a dedicated kitting station defined by a fixed subset of SKUs. Dynamic workload sharing within changing station borders, for instance enabled by the bucket brigade protocol (Bartholdi and Eisenstein 1996), or the relieve promised by job rotation (Otto and Scholl 2013) are not considered. • Assigning storage positions is a planning task executed over a mid-term planning horizon (De Koster et al. 2007). Thus, the composition O i of picking orders and their picking frequencies f i may be available, but not the specific order sequences. In the automotive industry, production sequences frequently need to be altered briefly before production starts, e.g., due to missing parts not timely delivered (Boysen et al. 2009). Thus, we cannot exactly quantify the picker's total walking distance, which also depends on the walk-back distances between successive jobs. Our order spreads, therefore, only measure one part of the total walking distance, namely the picker's company of an engine from its first to the last SKU container. In Sect. 4.3, we simulate the total walking distances by randomly drawing production sequences in order to check whether our proxy is indeed a good surrogate for the pickers' total walking distances. • As is the case at our engine manufacturer, we presuppose standardized SKU containers each having identical size. This allows us to calculate the objective value in number of passed containers. Many manufacturers aim to apply standardized containers, because this reduces their logistics and handling costs (Boysen et al. 2015). If, nonetheless, differently sized containers are applied, an order spread no longer solely depends on the first and last storage position, but also on the size of the containers placed in between. We rather aim at the most basic problem setting, which is also relevant for our real-world case. But adapting our models and solution approaches by this aspect seems truly straightforward. However, smaller containers are certainly another potential lever to reduce the pickers' walking effort, which comes at the price of more frequent replenishments. Evaluating this trade-off could be an interesting issue for future research. • We restrict our view on ground storage of containers along the conveyor, because this is the most basic setting. At our engine manufacturer this setting is applied to allow for quick container exchanges without double handling for repacking items or maneuvering containers stockpiled on top of each other. The alterations required to adapt our setting to the two-dimensional case are straightforward, so that we abstain from a detailed description. • We assume that the total in-line kitting area has enough space to house containers for all SKUs. If this is not given, depending on the parts demanded by subsequent workpieces, SKU containers must be added to and removed from stations on short notice. Given the short cycle times of the automotive industry these dynamic SKU swaps seem rather unrealistic, but they could be an interesting field for future research in other industries. • Missing parts in the assembly process produce excessive costs, in the worst case, the line needs to be stopped with hundreds of workers being idle (Boysen et al. 2009). Therefore, our engine manufacturer takes great care that there are always enough parts available in each SKU container. Consequently, we can abstract from the replenishment process and assume that there are always enough pieces in a container to satisfy demand. Now, we consider the computational complexity of WBSAP, which has the following complexity status.
Theorem 1 WBSAP is NP-hard in the strong sense, even if we have just a single kitting station |K| = 1.
The proof is by transformation from the linear arrangement problem (LAP), which is well-known to be strongly NP-hard (Garey and Johnson 1979) and stated as follows: LAP: Given a graph G = (V, E) and a positive integer M, is there a one-to-one- Proof The transformation scheme for generating an instance of WBSAP from an LAP instance is as follows. Since we have only a single station |K| = 1 , the complete workload is to be processed by a single worker, and the problem reduces to the storage assignment part seeking a sequence of SKUs within the station. For each node of LAP, we introduce a SKU, i.e., |S| = |V| , and for each edge, we generate an engine and its associated picking order with unit weight f i = 1 exclusively demanding the two SKUs represented by the adjacent nodes. Thus, we have a direct mapping between SKUs and nodes, picking orders and edges as well as storage positions and node numbers, so that a one-to-one mapping between both problems is readily available. ◻ Table 1 Notation for WBSAP-MIP In-line kitting for part feeding of assembly lines: workload…

Mixed-integer model for WBSAP
Given the notation summarized in Table 1, we are able to formulate WBSAP as a MIP with objective function (3) and constraints (4) to (14).

WBSAP-MIP:
subject to (3) minimizes the maximum weighted order spread over all kitting stations. Constraints (4) to (7) enable proper assignments of SKUs and storage positions: (4) assigns each SKU to exactly one storage position, (5) assigns exactly one SKU to each storage position, (6) assigns each storage position to exactly one kitting station, and (7) assigns at least one storage position to each kitting station. In inequalities (8), our z s,k -variables are aligned with x s,p and y p,k . Note that in case SKU s is not assigned to kitting station k, the z s,k -variables can attain either value 0 or 1, if this does not negatively affect solutions. (9) ensures connected kitting stations along the line. Specifically, if two positions p and p ′ are assigned to the same station, all other positions between p and p ′ must be assigned to this station as well. Constraints (10) define the order spread for each order and for each kitting station, and (11) defines the maximum order spread over all kitting stations, which is minimized in the objective function. Finally, variable domains are set in (12) to (14).
To further strengthen this formulation of WBSAP, we introduce several valid inequalities that eliminate symmetric solutions (so-called symmetry breakers): Constraints (15) sort the kitting stations along the line in increasing index order. In this way, solutions that only differ in the numbering of stations are eliminated. Note that we have also tried out the alternative formulation of these constraints of Ritt and Costa (2018), which can be borrowed from their model for assembly line balancing, on a small sample of instances. Since we were not able to detect any remarkable performance differences, we abstain from reporting these tests in detail. Inequalities (16) and (17) specify the assignment of SKU sets to stations: The indexes of the first SKUs in each station are sorted in increasing order. Hence, solutions that only differ in the assignments of SKU sets to stations are eliminated. Finally, (18) to (20) eliminate solutions that only differ in the direction of the SKU sequence within a station, since this direction is not relevant for the order spread. The first SKU in each station must have a smaller index than the second SKU in that station. Note that constraints (15) to (20) are only necessary, if we have more than one kitting station, i.e., |K| ≥ 2 . In case of |K| = 1 , (21) can be applied to define the direction of the SKU sequence: Preliminary computational tests have shown that all of these additional constraints, individually and in combination, reduce the solution time of a standard solver. However, the combination of all symmetry breakers performs best. Therefore, all further tests are conducted with symmetry breaking constraints (15) to (20) resp. constraints (21). Note that we also tried out another MIP based on variables x s,p,k , which directly encode the assignment SKU s to storage position p in kitting station k. This MIP, however, was outperformed by the one elaborated above, so that we abstain from reporting further details here and refer to Appendix A instead.

Decomposition approaches for the holistic problem
In order to solve a complex optimization problem such as WBSAP, decomposition is often a promising approach. We will follow this general idea in this section by decomposing WBSAP into its two core subproblems, storage assignment and workload balancing. If we have a solution for one of the subproblems, we will find it much easier to solve the other one. However, we still have to decide on the order in which we solve the problems, as both approaches go hand in hand with advantages and disadvantages: • Balance first, sequence second (BFSS): In this approach, we first decide on the assignment of SKUs to kitting stations. Given this assignment, we can determine the final SKU sequence for each station individually. The advantage of this approach is that the remaining sequencing problems can be small enough to solve them in reasonable time. However, this strongly depends on the amount and size of kitting stations since the remaining sequencing problem for a single station is still NP-hard. In Sect. 3.1, we introduce a meta-heuristic procedure for the workload balancing on the upper level and a beam search approach for the remaining storage assignment problem. • Sequence first, balance second (SFBS): If we have a sequence of SKUs along the storage area given, the remaining workload balancing problem reduces to the decision on sizing and positioning of kitting stations, a problem solvable in polynomial time. The downside of this approach, however, is the bigger size of (21) x s,1 + x s � ,|P| ≤ 1 ∀ s, s � ∈ S ∶ s > s � the SKU sequence to be determined upfront. In Sect. 3.2, we provide a metaheuristic procedure for the storage assignment problem on the upper level and a dynamic programming procedure for the remaining workload balancing problem.
Each of these approaches is dedicated a separate section in the following.

Balance first, sequence second
We explain this decomposition approach in reverted sequence. We start with the subsequent problem and elaborate how we determine a storage assignment problem for a given station workload. Then, we turn to the upper level and introduce a metaheuristic for determining the station workloads.

Storage assignment for a given station workload
This section defines the storage assignment problem (SAP) for a single kitting station once the workload for this station is already given. Thus, we have a single kitting station k operated by a dedicated picker and the set S k of SKUs to be placed within the station. The SKU containers have to be placed along the conveyor system, which has |S k | storage positions, so that each SKU receives exactly one position. Such a storage assignment can, thus, be represented by a SKU sequence k , i.e., a permutation of SKU set S k = {1, 2, … , |S k |} . The order spread SA i,k ( k ) of the picking order O i,k ⊂ S k associated with the SKUs of engine i ∈ I stored in station k is defined by: where function −1 k (s) returns the storage position of SKU s within storage assignment k . By aggregating all order spreads, we receive the total weighted order spread within the station: In-line kitting for part feeding of assembly lines: workload… Example (cont.): Consider the given workload of Station 1 within Solution three of Fig. 3d in Sect. 2.1. From the overall workload of six SKUs Station 1 only has to process SKU set S = {A, C, E, F} and the adapted order set depicted in Fig. 4a. Solution one and two depicted in Fig. 4b and c lead to a total order spread of 8 and 6, respectively.
Naturally, problem SAP can also be formulated as a MIP, which we present in Appendix B. In preliminary computational tests that are-for a matter of conciseness-not reported in this paper, however, the following approach proved much more successful. Recall that this problem has been shown to be strongly NP-hard within Theorem 1. Therefore, we introduce a dynamic programming (DP) scheme in the following, which is subsequently applied within a beam search heuristic. The following DP, which is an adaption of the basic sequencing DP of Held and Karp (1962), solves SAP to optimality.
Our DP consists of P k + 1 stages, each corresponding to a storage position in our line setup (plus a virtual starting stage). Each stage p includes states S k ⊆ S k , each defining a subset of SKUs stored in the first p storage containers along the line. The initial state is represented by an empty set S k,0 = � . The partial objective value of a state S k is denoted by z SA (S k ) and corresponds to the cumulative weighted order spread over all picking orders with respect to the first |S k | storage positions. Furthermore, we have a transition S k →S ′ k from state S k to state S ′ k if S k ⊂S ′ k and S � k ⧵S k = {s} with s ∈ S k , that is, the successive state contains the same SKUs as its predecessor plus an additional SKU. The additional weighted order spread for such a transition amounts to With the transition's contributions to the objective value on hand the Bellman recursion is defined by with z(�) = 0 . After a stage-wise forward recursion, the last stage P k + 1 contains only a single state S k = S k including all SKUs with objective value Z SA k = z(S k ) . Finally, we can extract the optimal storage assignment by a simple backward recursion.
Regarding the computational effort of our DP, we have O(2 |S k | ) states and O(|S k | ⋅ 2 |S k | ) transitions. The implied exponential runtime of O(|S k | ⋅ 2 |S k | ) is in line with our complexity result for SAP.
Due to the exponential runtime (in the number of SKUs), DP struggles with larger instances of our storage assignment problem. To accelerate the procedure, we modify DP and introduce a heuristic beam search (BS) approach. BS requires less runtime and memory as it only branches the most promising (with respect to the partial objective value) states of each stage. is called the beam width. With this modification, runtime and memory are polynomial bounded. However, optimality is not guaranteed hereby, since states leading to an optimal solution can be discarded during the forward recursion due to less promising partial objective values. Preliminary tests have shown that a beamwidth of = 100 performs well regarding solution quality and runtime for smaller and larger instances as well.
Example (cont.): Recall the example from above (see Fig. 4), where SKUs A,C,E, and F have to be sequenced within Station 1. The resulting beam search graph with a beamwidth of = 2 is depicted in Fig. 5. The procedure determines the four equally good-and in this case optimal-SKU sequences ⟨C,  In-line kitting for part feeding of assembly lines: workload…

Workload balancing on the upper level
With the approach for solving the sequencing problem for a given kitting station k at hand, we can now tackle the balancing problem on the upper level. For this purpose, we implemented a straightforward simulated annealing (SA) approach to determine the assignment of SKUs to kitting stations. SA is a stochastic metaheuristic inspired by thermal processes for obtaining low-energy states in heat baths. Based on the probabilistic acceptance of neighboring solutions, SA is able overcome local optima (e.g., Kirkpatrick et al. 1983;Aarts et al. 1997).
Our SA operates on an array of sets (S 1 , … , S |K| ) , whereby set k contains all SKUs dedicated to station k. Hence, each SKU must be element of one of the sets, i.e., ⋃ k∈K S k = S , each set must contain at least one SKU, i.e., |S k | >= 1, ∀k ∈ K , and the sets must be disjunct, i.e., S k ∩ S k � = �, ∀k ≠ k � ∈ K . The initial solution is determined by randomly assigning SKUs to sets, such that the above conditions are fulfilled. For obtaining a neighboring solution, we randomly perform either a swap move or a switch move: • Swap: Randomly choose two SKUs from different stations and swap their assignment to sets. • Switch: Move a randomly chosen SKU from one set to a different, also randomly chosen one.
Given the new (neighboring) assignment of SKUs to stations, we can obtain the weighted order spread Z k for each station k by applying the BS approach introduced in the last section. The objective value of our (balancing) solution Z(S 1 , … , S k ) can be determined by calculating the maximum weighted order spread over all stations Z(S 1 , … , S k ) = max k∈K Z k . The decision whether or not a neighboring solution (S � 1 , … , S � k ) should be accepted is made according to the following traditional probability scheme (see Aarts et al. 1997): If accepted, our current solution (S 1 , … , S k ) is replaced by neighbor (S � 1 , … , S � k ) as the new starting point for further iterations.
We applied a simple static cooling schedule (see Kirkpatrick et al. 1983) for steering our SA. The initial value T init for our control parameter T, the temperature, is given by T init = Z((S 1 , … , S k ) init ) , where (S 1 , … , S k ) init denotes our randomly generated initial solution. For each value of T, we perform three iterations of constructing and comparing neighboring solutions. Afterward, T is decreased by multiplying it with the factor 0.999. We continue this procedure until the T <= 0.001 ⋅ T init . If the given time limit for the procedure has not been reached, we restart SA with a new random solution and the initial temperature up to a maximum of five times.
Note that in our computational study, we invariably used control parameter values as described above as preliminary studies indicated that this parameter constellation outperforms other settings and delivers a reasonable compromise between solution quality and time.

Sequence first, balance second
The next decomposition approach solves the two subproblems in opposite order, which we again describe in reverted order. First, we elaborate the workload balancing problem for a given SKU sequence, and only afterward, we present a heuristic procedure to determine such a sequence.

Workload balancing for a given storage assignment
If the sequence of SKUs along the line is given, the remaining workload balancing problem (WBP) simply decides on the range of storage positions covered by each kitting station. Hence, we look for a solution WB = (l 1 , l 2 , … , l |K| ) , i.e., a sequence of storage positions with l k defining the last storage position of station k ∈ K . Such a solution is feasible, if l k+1 ≥ l k + 1 with l 0 = 0 , that is each station contains at least on storage position. With the given SKU sequence , we are now able to determine the order spread for each order in kitting station k: with l 0 = 0 . The final objective value, i.e., the maximum of weighted order spreads across the kitting stations, can then be determined by In-line kitting for part feeding of assembly lines: workload… Naturally, WBP can also be modeled as a MIP. However, since this problem can be solved to optimality in polynomial time, we present such a MIP only in Appendix C. Instead, we elaborate a more efficient DP approach in the following.
This DP consists of |K| + 1 stages, each corresponding to a kitting station (plus a virtual starting stage). Each stage k includes states (k, p) that each define the ending position p of kitting station k. The initial state is represented by (0, 0). The partial objective value of a state (k, p) is denoted by z WB (k, p) and corresponds to the maximum weighted order spread over the first k kitting stations along the line. Furthermore, we have a transition (k, p) → (k + 1, p � ) � from state (k, p) to state (k + 1, p � ) , if p < p ′ , that is, the ending position of a station is smaller than the ending position of the subsequent kitting station, and p � ≤ |P| − |K| + k + 1 , that is, the remaining number of positions is not smaller than the remaining number of kitting stations. The order spread for the new station added by such a transition is given by w p+1,p � . With the transition's contributions to the objective value on hand the Bellman recursion is defined by with z WB (0, 0) = 0 . After a stage-wise forward recursion, the last stage |K| contains only a single state (|K|, |P|) with objective value Z = z WB (|K|, |P|) . Finally, we can extract the optimal workload balancing by a simple backward recursion.
Regarding the computational effort of our DP, we have O(|K|) stages each containing O(|P|) states. Since the number of transitions leaving a state is in O(|P|) and |P| = |S| , the resulting runtime is in O(|K| ⋅ |S| 2 ) , thus polynomial. However, in order to apply the DP procedure, the values w p,p ′ have to be determined beforehand, which requires a runtime in O(|I| ⋅ |S| 3 ) . Since |K| ≤ |S| , the total runtime is in O(|I| ⋅ |S| 3 ) , thus polynomial. This allows us to solve even larger sized instances in a short runtime.
Example: Recall the example above depicted in Fig. 6 with the given SKU sequence = ⟨E, A, C, F, D, B] and |K| = 3 kitting stations. The aggregated weighted order spreads are given by Applying the introduced DP procedure, we obtain the graph depicted in Fig. 7. The optimal solution with an objective value of 5 is highlighted by the bold transitions. Here, the first kitting station includes the storage positions 1, 2, and 3 (state (1, 3)) with a workload of w 1,3 = 5 , the second station includes positions 4 and 5 (state (2, 5)) with a workload of w 4,5 = 5 and the last (third) station includes only position In-line kitting for part feeding of assembly lines: workload… 6 (state (3, 6)) with a workload of w 6,6 = 3 . The maximum workload over all three stations is max{5, 5, 3} = 5.

Storage assignment on the upper level
In the first stage of our sequence first, balance second approach, we address the sequencing problem, i.e., the ordering the SKUs along the kitting area. Again, we implemented a straightforward simulated annealing heuristic, which operates on a permutation (s 1 , s 2 , … , s |P| of SKU set |S|, with s p defining the SKU stored at storage position p. To obtain an initial solution, we simply assign SKUs randomly. A neighboring solution (s � 1 , … , s � |P| ) is then determined by randomly performing one of the following moves: With the new (neighboring) sequence on hand, we can determine the weighted order spread Z(s � 1 , … , s � |P| ) of our (sequencing) solution by applying the DP approach introduced in the previous section. The probability for accepting the neighboring solution-and therefore replacing the former solution-is given by

Computational study
First, we evaluate the computational performance of our different solution methods in Sect. 4.2. Then, we answer the question whether our surrogate objective (i.e., the order spread) is indeed a suitable proxy for the workers' actual total walking distances (see Sect. 4.3). Furthermore, we address several managerial aspects and analyze the benefit of our optimization approach in Sect. 4.4. Beforehand, however, we have to elaborate on our data instances in Sect. 4.1.

Instance generation
Since there is no established testbed for our WBSAP and confidentiality issues prohibit the application of real-world data, we had to generate our own instance generator. To do so, Table 2 lists the main parameters that are handed over to our generator as its own input data. Table 2 lists the parameter values for four different sets of instances employed during the subsequent sections. For a given combination of parameter values, each single instance is generated as follows: For given values of |S|, |I|, and |K|, we construct a set of different orders that refer to different engine types to be produced. For being able to generate order sets that are either more homogeneous or rather heterogeneous, we first derive a reference order with max(1, rd(|S| ⋅ r)) randomly chosen SKUs, whereby rd(x) yields the closest integer to x. Subsequently, we construct |I| orders O i , whereby each SKU (not) contained in the reference order is also demanded by an order O i with a probability of 1 − 2 ( 2 ). Smaller (larger) values of , with ∈ (0, 1] , thus lead to rather homogeneous (heterogeneous) orders. As each engine type demands an individual set of parts, we make sure that all orders are different. To ensure this, |I| must be smaller or equal to 2 |S| − |S| − 2 . If the constructed order is already contained in our order set, we dispose it and restart order generation. We also make sure that each order contains at least one SKU and that each SKU is used at least once. We, then, determine the |I| order frequencies randomly. Hereby, we presuppose a 7-hour working shift and a cycle time of 90 s, which equals the situation at our engine producer. Therefore, during a shift, about 280 engines are produced. We determine order frequencies f i ∈ (1 − , 1 + ) and normalize them, so that ∑ i∈I f i ≈ 280 . Hence, we set

3
In-line kitting for part feeding of assembly lines: workload… Analogously to our order set, smaller (larger) values of , with ∈ [0, 1) , lead to rather homogeneous (heterogeneous) frequencies. Note that due to rounding the total number of engines may differ slightly from 280. For the first three instance sets (Small, Large, and Surrogate), we repeat the construction process five times for each parameter combination, resulting in sets of 750, 225, and 625 instances, respectively. For the fourth set (Managerial), we only vary one parameter at a time and fix the others to their default values listed in parentheses within Table 2. We repeat the generation process 25 times for each parameter combination and obtain 1425 instances in this data set.

Computational performance
To evaluate our solution methods for WBSAP, we benchmark their computational performance for differently sized sets of instances: Small instances are still manageable by a commercial solver (to be solved to proven optimality in acceptable time) and large instances represent real-world problem sizes (see Table 2, Columns 'Small' and 'Large').
We solved the 750 small and 225 large instances each with both the balance first, sequence second (BFSS) algorithm and the sequence first, balance second (SFBS) procedure. Our industry partner currently determines a random SKU assignment to storage positions and assigns all stations the same amount of SKUs. We dub this approach RND. All procedures were implemented in Visual Basic, and tests were performed on a 64-bit system on an Intel(R) Core(TM)2 Quad CPU with 2.83 gigahertz and 8 gigabytes memory. Moreover, we apply optimization software 'Gurobi' (GRB, version 9, see Gurobi Optimization, LLC 2021) with a time limit of 300 s (900 s) for small (large) instances. The criteria with which we evaluate computational performance is defined at the bottom of Tables 3 and 4. The results summarized in these tables suggest the following findings.
• Gurobi: Among the 750 small instances, Gurobi solves 443 (59%) to proven optimality within the given time frame of 300 s. However, even for the small instances, Gurobi regularly hits the timeout when more stations |K| and SKUs |S| are involved. On average over all small instances, Gurobi is clearly outperformed by our decomposition approaches both with regard to solution quality and runtime. This finding becomes even more pronounced for the large instances. Here, Gurobi does not terminate prior to the timeout of 900 s in any instance, and for 109 out of 225 instances (48.4%), Gurobi is not even able to find a feasible solution during this time. Hence, we conclude that Gurobi is not suitable for instances of real-world size, and heuristic solutions seem preferable instead. • BFSS vs. SFBS: We can observe no clear dominance of either of our decomposition approaches. In some cases, BFSS delivers better results than SFBS; in others, it is the other way round. This holds true for both runtime and solution quality. On average, however, SFBS leads to better results regarding both test sets, small and large, and both performance criteria, runtime and solution quality. Thus, SFBS seems the better option. In-line kitting for part feeding of assembly lines: workload… #opt: Number of instances (out of 25) solved to optimality, gapOpt: Optimality gap, #best: Number of instances (out of 25) where the method performed best compared to others, gapBest: Objective gap to best performing procedure  Table 4 Performance results of Gurobi (GRB), random approach Rseq, balance first, sequence second (BFSS), and sequence first, balance second (SFBS) on large instances #noSol: Number of instances (out of 15) without feasible solution within time limit, #best:

Number of instances (out of 15)
where method performed best compared to others, gapBest: Objective gap to best performing procedure In-line kitting for part feeding of assembly lines: workload… • Real-world method: Finally, we can observe that the status quo method applied by our engine producer seems not advisable. Although its runtime requirement is barely measurable, it produces gaps to the best solutions of up to 47.81%. How this translates into walking effort for the pickers is investigated in more detail in Sect. 4.4. Note that a single random SKU sequence is a rather weak competitor from an optimization perspective. In Appendix D, we explore the performance impact if more than a single random solution are drawn. Our results show that especially for large instances such a simple approach cannot compete with our more sophisticated optimization approaches.
To conclude, our results suggest that both our decomposition approaches BFSS and SFBS deliver good solutions in a short amount of time, so that they seem wellsuited even if large instances of real-world size are to be solved. SFBS is a bit faster and delivers better average objective values, so that for all further tests SFBS is our method of choice.

On the appropriateness of the surrogate objective
In this section, we aim to answer the question whether our surrogate objective, namely, minimizing the maximum weighted order spread, is indeed a good proxy for the actual objective, which is to minimize the workers' actual maximum walking distance. Recall that when planning the storage assignment, which is done on a midterm basis (e.g., twice a year), we do not know the actual daily production sequences of engines. Instead, we only have reliable forecasts on the frequencies, in which each engine type is produced. Note that our engine producer has long-term contracts with car manufactures and rather long lead times, so that aggregate demands for the next months are indeed well predictable in this case. This allows us to quantify the distances each worker has to accompany the engines over the planning horizon with the help of our surrogate objective. However, we cannot add the walk-back distances between subsequent engines, because they depend on the production sequences. In this section, we investigate if our proxy, in spite of this inaccuracy, still supports our 'true' objective.
To do so, we evaluate 625 instances, including 25 different settings of station numbers |K| ∈ {1, 2, … , 5} and SKU numbers |S| ∈ {10, 20, … , 50} , each generated with the data generator defined in Sect. 4.1, applying the parameter settings stated in Column 'Surrogate' of Table 2. For each of these instances, 100 random WBSAP solutions, each consisting of a random storage assignment of SKUs along the line and random station borders, are obtained. Each solution is evaluated with our surrogate objective to determine the resulting maximum order spread over all kitting stations. These values are compared to the actual maximum walking distances that include the walk-back distances. To obtain them, we draw 250 production sequences per instance, where each engine type occurs in this sequence according to the given production frequencies. For a given production sequence, the actual total walking distance per station and thus the maximum over all stations can easily be obtained. Finally, we average the maximum actual walking distances over all 250 production sequences. For the resulting 2 ⋅ 100 values, namely of our surrogates and the actual objectives, we calculate Pearson's correlation coefficient. The minimum, average and maximum correlation coefficients, rounded to six digits, are listed in Table 5 for all 25 |S|-|K|-parameter settings. To get an idea on the absolute values of our surrogate objective and the (simulated) total walking distance (in m), we also present the respective average values in the table as well. Note that the actual forward distance (in m) can be determined by 0.8 ⋅ (objective value) as SKUs have an assumed width of 80 cm.
For every parameter combination, we observe correlation coefficients close to 1. This result verifies that we have a strong positive correlation between our surrogate 1 3 In-line kitting for part feeding of assembly lines: workload… 'maximum weighted order spread' and actual objective 'maximum actual walking distance.' We can thus conclude that our proxy seems well-suited and, if applied to guide an optimization task, can successfully support the reduction in the workers' actual walking distances.

Managerial issues: reduction in walking effort
In this section, we evaluate whether our optimization task, i.e., optimizing the storage assignment and division of labor among kitting stations, actually has the potential to considerably reduce the workers' walking effort. To explore this, our computational experiment is designed as follows. We start with the 1425 instances generated with our instance generator (see Sect. 4.1) and the parameter settings from Column 'Managerial' in Table 2, and derive WBSAP solutions with the status quo method RND of our engine producer (i.e., random storage assignment with equally sized stations) and our best-performing decomposition approach SFBS for each instance. To obtain the actual maximum walking distances for each of these solutions that also include the walk-back distances, we again derive 250 random production sequences of 7-hour shifts with engine occurrences according to the given frequencies. The resulting actual maximum walking distances are finally averaged over all 250 production sequences. First, we take a look on the impact of parameters |S| (number of SKUs), |K| (number of kitting stations), and |I| (number of engine types to be produced) when benchmarking our two competitors. In Fig. 8, we report the actual maximum walking distances of our two solution methods depending on the above three parameters. These results suggest the following findings: • Impact of SKUs: First, we can confirm an expectable result in Fig. 8 (left). If the engines require more parts and |S| increases, then more SKU containers have to be arranged along the line. Obviously, this leads to a linear increase in the walking distances. • Impact of stations: Figure 8 (middle) displays the impact of another lever to reduce the walking effort of the workers: additional kitting stations. This, however, also increases the kitting workforce and thus wage costs. We can observe diminishing returns of additional workers. When having a small workforce every additional worker leads to a considerable relieve. If however the workforce is already substantial, the extra relieve of yet another worker is small. • Impact of engine types: Finally, Fig. 8 (right) shows that the number |I| of different engine types that are produced on the line has only moderate impact. It is not the sheer number of engines but rather the level of heterogeneity of their part demands that impacts the walking effort, which can be seen in our further analysis.
Our previous results of Fig. 8 have shown a consistent improvement of our optimization approach over status quo method RND. We analyze these gains in more detail with the help of Fig. 9. Here, we display the percental reduction of the actual maximum walking distances of our optimization approach SFBS compared to the baseline of status-quo method RND of our engine producer. The results of Fig. 9 suggest the following findings: • Impact of reference order: During the discussions with our practice partner, we had to learn that in engine production there are quite a few basic parts that are required by most engine types. The reference order represents these parts. A larger r-value (with default = 0.5 ) mainly results in larger order sets per engine (see Fig. 9 (left)). In this case, large parts of each kitting station must be passed in each cycle anyway, so that optimization cannot gain much. If only a few SKUs are needed (i.e., at a smaller r-value), optimization has much more flexibility to cut the accompanying walks for specific engines short. • Impact order heterogeneity: Parameter controls the probability with which a SKU of the reference set is contained in a specific order (see Fig. 9 (middle)).
Thus, a small (large) leads to more (fewer) SKUs of the reference set being contained in each order. Thus, smaller -values lead to rather homogeneous SKU demands of all engine types. Hence, we benefit substantially by optimiz- ing the arrangement of the most common SKUs. Optimization can build clusters of SKUs jointly demanded in similar orders and the workers must only traverse subsections of their stations for each engine. In this case, the walking effort of pickers can considerably be reduced by up to 30%. For more heterogeneous orders, it is harder to find SKU arrangements that lead to small order spread for each order, and therefore, the benefit of optimization decreases with increasing values of . • Impact of frequency: In Fig. 9 (right), we can observe that more or less variation in the frequencies, in which the engine types are produced (i.e., regulated by parameter ), has only negligible impact. The improvements of optimization over the status quo are consistently over 20 %. This result indicates that when applying the default values for all other parameters, this obviously leads to rather homogeneous demands for parts among different engine types. Hence, varying their frequencies still leads to similar part demands.
To conclude, our study reveals the following main take-home message. On average over our instances, sophisticated optimization can reduce the maximum walking distances of pickers by more than 20 %, which is good news for both workers and their employers. For workers, it makes a considerable difference, especially on their long-term well-being, if their daily walking distances reduce by several kilometers. Healthy and more satisfied workers are certainly also important for their employers, especially in the aging societies of many industrialized countries. Furthermore, these gains can be realized without more wage costs for additional workers or a reduction in the output; merely a low-cost rearrangement of SKU containers along the line and altered station borders are required.

Conclusions and outlook
In this paper, we introduce a novel part feeding policy to the scientific literature: in-line kitting. The first part of an assembly line is reserved to part logistics. Here, pickers have to add the requested parts that are arranged in containers along the line to traveling part kits that move with their workpieces on the conveyor. To reduce the ergonomic effort for these pickers when having to walk along their kitting stations, we formulate a novel optimization problem that arranges the SKU containers along the line and determines the borders of kitting stations. We introduce different solution methods, and our computational performance tests identify decomposition approach sequence first, balance second (SFBS) as an appropriate approach even if large instances of real-world size must be solved. In further tests on managerial issues, we show that applying this method can reduce the walking distances of pickers by about 20 % without producing any additional costs or reducing the output.
Future research could challenge our solution methods and provide even more powerful (especially exact) optimization approaches. Furthermore, it could be beneficial to provide the same SKU not only at a single kitting station, but also at multiple ones. This adds an operational selection problem, from which station a specific SKU demand should be satisfied, and increases the space demand for kitting. On the other hand, duplicate SKUs could-if properly placed-further reduce the walking effort of the workers. Furthermore, reducing the walking effort could not only lower the physical effort of the existing workforce. These savings could also be used to improve productivity and reduce the kitting workforce. The resulting balancing problem that minimizes the workforce for a given cycle time is an interesting challenge for future research. Finally, future research should evaluate whether in-line kitting could also be a suitable alternative for other production lines (beyond engine production). This suggests an extensive benchmark study of all alternative part feeding policies in different production environments.

Appendix A: An alternative MIP for WBSAP
This appendix presents an alternative MIP model, which due to its three-dimensional assignment variables seems better accessible but turned out non-competitive. Applying the notation summarized in Table 6, this MIP (dubbed WBSP-MIP2) consists of objective function (32) subject to constraints (33) to (46)  In-line kitting for part feeding of assembly lines: workload…

SAP-MIP:
subject to Objective function (47) minimizes the weighted order spread in kitting station k summarized over all picking orders. Constraints (48) and (49) ensure that each storage position receives exactly one SKU and vice versa. The order spreads for all picking orders are dimensioned by (50). Here, each pair of SKUs per order is evaluated, and the maximum distance among any such pair defines the order spread. Finally, constraints (51) set the domain of the binary variables. We have carefully tested the application of this MIP within our decomposition approach. However, especially for larger instances solving this model took default solver Gurobi too much time to evaluate a large number of storage assignments. Our beam search approach delivers only heuristic results but is much faster, so that it turned our as the better option for this stage. For a matter of conciseness, we abstain from a detailed elaboration of these computational results.
) x s,p ∈ {0, 1} ∀s ∈ S k ; p ∈ P k We applied these procedures to the small and large set of instances from Subsect. 4.2 for #sol ∈ {1, 10, 100, 1000, 10 4 , 10 5 , 10 6 } . The results are presented in Table 9. Note that to improve readability, we only list the average values for the entire data set rather than the results for each parameter combination. gapBest denotes the average objective gap to the best performing procedure.
The results of Table 9 clearly indicate that even drawing a significant amount of random sequences, so that the solution times are about as high as those of our more sophisticated optimization approach SFBS, is not competitive. This supports our choice of applying method SFBS (see Sect. 3.2) for solving WBSAP.
Funding Open Access funding enabled and organized by Projekt DEAL.

Data availability
The datasets generated and analyzed during the current study are available from the corresponding author upon reasonable request.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.