AN APPROXIMATE COST-INTEGRATED MODELLING APPROACH TO CROSS DOCK DESIGN

Cross docking is a form of operating break-bulk freight terminals in which the receiving products are not stocked but, instead, the incoming cargo is internally transferred to be immediately shipped in distribution trucks. With this approach, inventory costs and lead-times are significantly reduced. In this work, an approximate, continuous, non-stationary queueing model is established to estimate the number of cargo discharging doors in the receiving dock. Since the traditional first-come, first-served queueing discipline is not applicable to this specific problem, an alternative formulation, called least-distance, firstserved queue discipline, is adopted and analysed with a simulation approach. The problem is modelled to support a preliminary design of a cross-dock terminal. The objective function, to be minimized, is the sum of the rent cost of the terminal facility, the waiting cost of inbound trailers when being attended at the terminal, and the cost of handling and moving freight inside the terminal.


INTRODUCTION
Freight terminals, in the less-than-truckload (LTL) transport industry, are facilities in which shipments are unloaded, checked, sorted, consolidated, and ultimately loaded into outgoing trailers for delivery elsewhere in the supply chain.In addition to being costly, freight handling is important because the time that a shipment spends at the terminal is wasted in the sense that the shipment is not making progress toward its final destination (Gue, 1999).A special freight terminal operation is cross docking, a logistics strategy used by companies in different industries for transferring incoming product shipments directly to outgoing vehicles without significantly storing goods in between.Inbound trailers or trucks bring products from suppliers, which are unloaded, moved across the terminal and finally loaded into outbound trucks, which immediately leave the terminal to deliver the cargo to their ultimate destinations in the distribution network (Boysen & Fliedner, 2010).This practice has a number of objectives: the consolidation of shipments from diverse suppliers, the shortening of delivery lead-time, reducing inventory and order picking costs, etc. (Boysen & Fliedner, 2010;Van Belle et al., 2012;Buijs et al., 2014).In contrast to a traditional warehouse, the storage of goods is reduced to the greatest possible extent.As a result, the cross docking approach enables the consolidation of many less-than-truckload (LTL) shipments between multiple shippers and recipients, thus generating economies of transportation and inventory.Consequently, hub-and-spoke distribution networks, replacing traditional pointto-point deliveries, have been made possible (Stephan & Boysen, 2011).But handling freight in an LTL terminal is labour intensive and therefore expensive, since workers must quickly sort and move a variety of freight.An important way to reduce such work is to assign inbound trailers and outgoing trucks to rationally chosen doors of the terminal to take advantage of patterns of freight flows (Bartholdi & Gue, 2004).
As stated, the crossdocking terminal is a special type of warehouse.A traditional warehouse has four basic functions: cargo receiving, storage, order-picking and shipping.First, a bundle of freight is received at a warehouse, and then placed in a storage location.Next, order-picking occurs when an item is requested by a customer.Finally, shipping is activated when orders are consolidated and moved to trailers.Among all activities, storage and order picking together account for about 70% of warehouse operation expenses (Bartholdi & Hackman, 2014).The crossdocking system, like the just-in-time (JIT) system, is designed to reduce inventory and labour processing, and also aims to attain transportation economies within shorter order cycle times.Crossdocking attains these benefits by transferring shipments directly from inbound trailers to outbound trailers with minimum storage in between.Usually, cargo is shipped within 24 hours, though sometimes it could be less (Apte & Viswanathan, 2000).
Typically, a small or medium cross-docking terminal is composed of a rectangular open-area building with several dock doors located around its perimeter (Fig. 1).Other terminal layouts and formats, with X, L, T shapes, as well as some other forms, are also encountered in the real world and analyzed in the literature (Bartholdi & Gue, 2004, Van Belle et al., 2012), but the rectangular shape is predominant.Part of the dock doors is assigned to inbound trucks that bring freight to be unloaded.Another group of doors is reserved to outbound trucks, which will perform the immediate delivery of products to their destinations.Inbound vehicles are usually of a trailer type, whereas outbound vehicles are mostly trucks.But both terms are indifferently used for both situations in this text.The two dock door types are named strip doors, where trailers are attached to the dock to be unloaded, and stack doors, where outbound trucks are located to be subsequently loaded (Gue, 1999).A vehicle in a stack door receives freight for a same destination, or for a pre-defined group of destinations (for example, a group of stores in a supermarket chain), regardless of their origins or contents (Bartholdi & Gue, 2000).
When an inbound truck arrives at a cross dock terminal and all of the strip doors are occupied, it is assigned to a queue position on a parking place in order to wait to be served (Fig. 1).At due time, the vehicle is commanded to approach a strip door to start the cargo unloading process.
When the unloading process is under way, the cargo shipment is registered, checked, and sorted according to its destinations.Next, the cargo is moved inside the terminal with manual fourwheeled platform dollies, pallet jacks, forklift movers, or other material handling equipment, whose selection depends on the cargo characteristics and volume, up to a temporary staging area close to the respective shipping dock door.At this place, other services, such as labelling and addressing may be executed while the cargo waits to be loaded into the outbound trucks.As soon as an outbound vehicle is ready, it leaves the terminal to its final destinations (Stephan & Boysen, 2011).In this paper, an approximate, continuous, non-stationary queueing model (Newell, 1971) is established to estimate the number of cargo discharging doors in the receiving dock.Since the traditional first-come, first-served queueing discipline is not applicable to this specific kind of problem, an alternative formulation, called least-distance, first-served queue discipline, is adopted and analysed with a simulation approach.The objective function, to be minimized, is the sum of the rent cost of the terminal facility, the waiting cost of inbound trailers when being attended at the terminal, and the cost of handling and moving freight inside the terminal.The internal transport of cargo is alternatively performed with four-wheel platform dollies, or with electric pallet trucks.The problem is modelled to support a preliminary design of a cross-dock terminal.

LITERATURE REVIEW AND MOTIVATION
A great number of research work on cross docking is concerned with dock door assignment problems (Van Belle et al., 2012).It is because dock door assignment deeply affects the material flow and the freight carrying costs inside the facility.But, regarding truck scheduling problems, reported applications generally assume that either all inbound vehicles are present at the cross-dock facility at the beginning of a working period, or the time each truck arrives at the cross-dock facility is known in advance (Sadykov, 2009;Vahdani et al., 2010;Alpan et al. 2011; Konur & Golias, 2013;Van Belle et al., 2013).However, the assumption that truck arrival times are known beforehand with certainty is unrealistic, as truck arrivals are subject to uncertainties due to delays at the suppliers' premises, traffic congestions, weather conditions, or maintenance drawbacks (Boysen & Fliedner, 2010;Konur & Golias, 2013).As a consequence, current applications generally cover the problem only partially, reflecting realistic situations and maintaining the analysis restricted to some specific aspects, usually on the operational decision level (Agustina et al., 2010).Instead, the research work object of this paper assumes a tactical decision framework, with the objective of producing robust results, albeit approximate, to support a preliminary design of a cross dock terminal (Berghman & Leus, 2015).Such a preliminary design is normally intended to cover a gap between a general design conception and the ultimate detailed design, particularly when the level of initial conceptualization is not sufficient for a full evaluation of the project, mainly in economic terms.Approximate modelling is thus employed, covering the essential parts of the problem, as dock door assignment, truck scheduling and waiting time, as well as an integrated minimal cost formulation.In Brazil, although no specific works on crossdocking are available in the local OR literature, some articles cover related subjects, as integrated production and distribution problems (Reimann et al., 2014), and logistics cargo terminals (Vieira & Luna, 2016).The reviews are quite ample and provide a classification analysis of the most important problems encountered in the literature on cross docking modelling and applications.Diverse types of problems, covering strategic and tactical decision approaches are described and analysed, such as location of cross dock terminals, cross docking networks, facility layout design, as well as operational problems such as vehicle routing, dock door assignment, truck scheduling, and temporary storage (Van Belle et al., 2012).However, a large number of papers are not directly related to our work, and therefore will not be referred here.Readers interested in a more general review of the problem shall consult Agustina et al. (2010), Boysen & Fliedner (2010), and Van Belle et al. (2012).

Inbound trailer waiting time
Queue modelling is a traditional branch of operational research, with reported application areas being healthcare, transport, and telecommunications, among others (Worthington, 2009).It involves analytical investigation of processes that lead to people, or other element, to wait before starting to be served.Its objective is to determine and analyse performance measures, which express the efficiency and effectiveness of these processes.Among those measures are: expected number of elements in the queue, average waiting time to be served, idle time of service providers, etc. (Fogliatti & Mattos, 2007).The average queue size and the time a truck spends in a queue waiting to be served is usually estimated with a mathematical queueing model conveyed in basic texts such as Prabhu, 1965;Page, 1972;Gross et al., 1985;Fogliatti & Mattos, 2007.A more recent review on queue modelling is presented by Worthington (2009).

Mathematical queue modelling is usually based on Markov Chains (Fogliatti & Mattos, 2007).
The queue, when is ergodic, converges to a stochastic equilibrium situation called a stationary state, in which state probabilities are independent of time (Fogliatti & Mattos, 2007).Such classical models, however, usually do not apply to the analysis of queues formed by incoming trucks to be unloaded at cross-docking facilities.In fact, in most cross-docking terminals the queueing process does not reach a stationary state.Three main reasons support this statement.First, cross dock terminals are usually operated within pre-defined time windows.A trailer arrives at a random time t respecting the time window t 0 ≤ t ≤ t E established by the terminal management, and within this time interval there is no guarantee that the stationary state will be attained.Second, the truck arrival rate may vary substantially over the time-window period, a situation quite different from the classic theoretical approach, in which the average arrival rate is constant.Finally, when the receiving strip doors are closed at time t E , although no additional incoming trailers are allowed, the unloading process may continue for a while in order to serve the vehicles still in the queue.
An interesting approximate approach to queueing investigation is due to Gordon F. Newell (1971), who developed a methodology to analyse non-stationary queues with simple, straightforward techniques, under a continuous approximation approach.Whereas most of the literature on queues deals with equilibrium distributions of queue lengths and waiting times, the main emphasis of the Newell's approach is on time-dependent behaviour, as in this paper, where the arrival rate of users changes along time.The techniques employed by Newell (1971) are mainly fluid and diffusion approximations.Much of the classical theory of queues deals with discrete models, but in simple continuous approximations such as Newell's, it is convenient to initially disregard the discrete nature of some elements, treat them as continuous or quasi-continuous variables, and return back to the discrete representation afterwards.
Another approach to queue modelling is simulation (Worthington, 2009).Today, due to advanced commercial simulation software, associated with technological evolution in computational techniques, many queue modelling applications make use of simulation models.A good number of papers using simulation techniques to solve cross dock operational problems are presented in the literature ( Shi et al., 2013).Such a technique will be partially used in this work because the traditional First-Come, First Served (FCFS) queue discipline assumed in Newell's (1971) approach is not adequate to represent the priority scheme of attending incoming trailers to be unloaded at strip doors (see Section 3.3).

Dock door assignment
Efficient operation of a cross-dock system requires an appropriate coordination of inbound and outbound flows.In fact, a substantial part of the cross-dock operating cost is due to the intense movement of workers and material handling equipment across the facility, when freight is transferred from incoming to outgoing trailers, and this operation is expensive (Gue, 1999).As a response to this, quite a number of quantitative modelling approaches have been developed and applied with the objective of defining adequate terminal layouts, in which one searches for the best arrangement of strip doors and stack doors (Gue, 1999; Miao et al., 2009Miao et al., , 2014;;Boysen & Fliedner, 2010;Vis & Roodbergen, 2011;Acar et al., 2012;Gelareh et al., 2016).
When an inbound trailer arrives at the cross dock terminal, it is necessary to decide to which strip door the vehicle should be assigned.On the other hand, outbound trucks must be assigned to stack doors.Thus, the dock door assignment problem aims to find the optimal positioning of strip and stack doors along the facility perimeter, with the objective of minimizing freight handling and displacement costs (Van Belle et al., 2012).A good door assignment can increase the productivity of the terminal and decrease the internal handling costs.
In some cases, as in this work, the assignment of outbound trucks to stack doors is executed on a mid-term horizon, with the result that each stack door will exclusively serve a specific outbound destination (or group of destinations) for a somewhat longer period of time (a month, or a season, for example).Such a fixed assignment eases the allocation of shipments to the outbound trucks, consequently yielding a learning effect on workers' movements.Such a fixed stack door assignment to outbound trailers is especially suited for steady commodity flows, with reliable distribution of cargo among inbound flows and outbound destinations (Boysen & Fliedner, 2010), as in the case of this application.Thus, on a typical working day, each outbound truck is assigned to a previously determined stack door, and this setting pattern is maintained during a somewhat long period of time (Van Belle et al., 2012).
Two recent reviews of dock door assignment are due to Van Belle (2012) and Acar et al. (2012).An early and extensively cited study about the assignment of dock doors was presented by Tsui & Chang (1992).They formulated the door assignment problem with the objective of minimizing the total distance travelled to move freight inside the facility.Bozer & Carlo (2008) formulated the door assignment problem as a mixed integer programming model with the objective of minimizing material handling cost.Simulated Annealing was also used to solve the proposed model, taking into consideration congestion and robustness.
Cohen & Keren (2009) also extended the approach by Tsui & Chang (1992).The mathematical model was adapted to allow that freight, for a certain destination, can be split and delivered to multiple doors assigned to that destination (the capacity of the outbound trucks is taken into account).The proposed formulation is a non-linear MIP model that is impractical for real size problems.So, the authors proposed a heuristic algorithm to solve it.Acar et al. (2012) proposed a robust door assignment method focussing operational uncertainties in LTL operations.Time variability associated with vehicle arrival and service time is treated in the model acting on the idle time subsequent to each truck unloading service.

Truck scheduling
In a truck scheduling problem, the goal is to determine the assignment of incoming vehicles to strip doors, and outbound ones to stack doors (Agustina et al., 2014) 2011), on the other hand, assume that truck arrival times are given.However, assuming truck arrival times are known with certainty is unrealistic as truck arrivals are subject to delays at the suppliers' premises, uncertainties due to traffic congestions, variable weather conditions, or truck failures (Boysen & Fliedner, 2010).These assumptions, therefore, may result in suboptimal scheduling strategies in practice (Boysen & Fliedner, 2010), and affect cross-docking efficiency as waiting and delivery times heavily depend on the scheduling strategies (Wang & Regan, 2008).In a more recent review of cross-docking studies, Van Belle et al. (2013) note that the assumption of deterministic arrival times is a shortcoming of the truck scheduling studies encountered in the literature.One important feature of this paper is to admit random arrivals of incoming trucks (Section 3.4), therefore overcoming this kind of restriction and making the analysis more realistic.
Van Belle et al. (2012) divide the research works about the truck scheduling problem in three categories.A first category considers a simplified cross dock pattern with a single strip and a single stack door.Truck scheduling reduces, in this case, to sequencing the inbound and outbound trucks.The works in the second category consider cross docks with multiple strip and stack doors, but deal only with the scheduling of inbound trucks.The last category is related to papers that consider both inbound and outbound truck scheduling at multiple dock doors.The truck scheduling problem investigated in this paper belongs to the second category.The objective is to minimize the sum of three costs, namely the real estate cost of the facility, the waiting cost of inbound trucks before being unloaded at the dock, and the cost of handling and moving freight inside the terminal (Section 4.2).

Terminal layout
A strategic decision that has to be taken when starting a cross dock design is to choose its layout pattern.The layout is interpreted as the dimension and shape of the cross dock building, as well as the dimension and shape of the internal cross dock areas and their arrangement (Van Belle et al., 2012; Vis & Roodbergen, 2001).The layout problem of a cross dock terminal was investigated by Bartholdi & Gue (2004).In the analysis, the best shape of the terminal building (e.g., I, T, E, or X-shaped facilities) and the location of strip and stack doors are determined (Boysen & Fliedner, 2010).In practice, most of existing cross docks are long, narrow rectangles (I-shaped).Several experiments are reported in the literature in which the labour cost (estimated via the total internal travel distance) was measured for different shapes.The experiments suggested that an I-shape is the most efficient pattern for smaller cross-docks (fewer than about 150 dock doors).For docks of intermediate size, a T-shape is best, and for more than approximately 200 doors, an X-shape is best (Van Belle et al., 2012).
The layout of the terminal influences internal freight carrying costs in three ways.First, as the number of dock doors raises, the average distance between strip and stack doors also tends to increase, with parallel effect on freight handling costs.Second, workers are more productive when inbound trucks are nearer to the appropriate destination trucks, because this reduces travel distance and internal traffic conflicts.Thus, it is advised to assign incoming trailers to strip doors close to the destination vehicles for which they have the most freight.Third, the shape of the terminal (e.g.rectangular, I, T, X-shaped, etc.), as mentioned, also influences freight handling costs.
In addition to internal travel distances, Bartholdi & Gue (2000) introduced the congestion cost due to different types of traffic conflicts.Although of interest, most of the current published research on the subject takes into account just the internal travelling distance, taking this element as a proxy variable used to estimate the total cost of transferring freight from the strip doors to the shipping docks.In this paper, this approximation is also assumed.

PROBLEM FORMULATION
The problem studied in this paper belongs to the second category of truck scheduling analysis as indicated in Section 2.3, and considers a cross dock with multiple strip and stack doors, but dealing only with the scheduling of inbound trailers (Van Belle et al., 2012).In this case, the assignment of outbound trucks to stack doors is fixed for a somewhat long period of time (a month, or a season, for example), as indicated in Section 2.3.The objective function is the sum of three cost components: (a) the rent cost of the terminal building facility, which is related to its area and layout; (b) the waiting cost of inbound trucks when being attended at the terminal, and (c) the cost of handling and moving freight inside the terminal.These three costs are deeply interrelated, and must be analysed in an integrated way during the preliminary design phase of a cross docking facility (Rouwenhorst et al., 2000).The costs related to outbound trucks are considered fixed in our problem and therefore do not influence the search for the minimum terminal cost configuration.
The contributions of this paper are: (a) to develop a modelling framework of cross-docking basic operations in a simple and integrated way, involving different parts of the problem which have been treated separately in the literature; and (b) offering a simple, albeit robust, modelling approach that will allow researchers and practitioners to analyse prospective terminal configurations and compare alternative possibilities related to performance and costs.The basic assumptions of the truck scheduling problem are as follows (this listing is partially based on Van Belle et al., 2013): • An exclusive mode of service is considered, i.e. each dock door is either exclusively assigned to inbound or to outbound trucks; • A rectangular (I shaped) terminal is assumed (Bartholdi & Gue, 2004), with strip and stack doors assigned to either side of the facility (Fig. 1); • Arriving goods are unloaded from the inbound trailers and transferred to the appropriate outbound dock where they are ultimately loaded into outbound trucks.Other internal operations (like sorting and labelling) are not considered here; • Preemption of loading or unloading a truck is not allowed.Thus, a docked vehicle has to be completely processed before it leaves the dock; • Inbound trailers arrive randomly at the terminal, with arrival times following a known probability distribution (see Section 3.4); • Departure times of outbound trucks are established by terminal supervisors and are not included as decision variables in the model; • Freight is received in standardized cartons at the terminal; • Two alternative internal transport modes are considered: (I) freight in cartons, carried in four-wheeled platform dollies (Section 4.2.3-A);(II) cargo in pallets, with freight palletized upon reception at the terminal and carried by pallet trucks (Section 4.3.3-B); • The time needed (and the corresponding cost) to transfer freight from strip to stack doors is assumed to be directly proportional to the rectangular distance between the dock doors to which vehicles are assigned; • Intermediated storage inside the terminal is partially allowed.This means that goods can be unloaded from an inbound trailer before the appropriate outbound truck is available.The capacity of the storage is considered infinite; • The trailer changeover time is considered nil in the model; • Products are not interchangeable, i.e. any arriving product unit is assigned to a specific outbound truck; • The sequence in which goods are loaded or unloaded is not taken into account; • The problem consists of n V inbound trailers, m 1 strip doors, and m 2 freight destinations.
The main decision variable in the model is the number of strip doors.

Generating freight flows
In practice, for small and intermediate size cross-dock terminals, supervisors usually tend to assign inbound trailers to strip doors without explicit regard to the destinations of the freight, since usually they adopt the FCFS priority rule to unload the vehicles (Gue, 1999).Among other reasons to do so, they do not have a decision support system to help them in making better assignments, in which the average travelled distance within the facility is minimized, for instance.
Workers are more productive when incoming trailers are nearer the appropriate outbound trucks, because this situation decreases internal travel distances and consequently reduces labour costs.Thus, it is advised to assign incoming trailers to strip doors close to the destination trucks that will take the most freight (Gue, 1999).But, since inbound trailers usually carry different cargo mixes to different destinations, it is difficult for the supervisor to individually decide to which vacant strip door he will assign a loaded trailer without the support of a more elaborate evaluation methodology.Specially, when the average number of destinations per incoming trailer is low, and the dispersion of destinations varies significantly per incoming cargo load, the optimized assignment of strip doors based on cargo flow and internal distance can lead to significant savings in material handling costs within the terminal (Gue, 1999).
Real data on the mix of freight per inbound vehicle and per destination is hard to obtain in Brazil, due in part to proprietary reasons, but also because they are not collected and organized in a way as to allow their utilization in this kind of research.In addition, logistics operators in the LTL industry do not record the size of individual cargo segments within a shipment (Gue, 1999).Thus, and following Bartholdi & Gue (2004), empirical arrangements of freight flows carried by inbound trailers had to be specially generated through a Monte Carlo experiment, in order to supply quasi realistic data for modelling the problem.Freight is assumed to be transported in cartons.Specifically, in order to guarantee that the Monte Carlo generated configurations will show a relative low number of destinations per incoming trailer, as well as to have larger dispersion of destinations per vehicle as recommended by Bartholdi & Gue (2004), a specific Monte Carlo algorithm was developed to generate random loads of freight for the incoming trailers, which satisfy these requisites.
Since cargo is not palletized at the origin in this case, and usually is not loaded in carton units either, an intermediate situation was assumed, in which the trailers are loaded at the supplier's premises in carton packs of 20 units each.Since each destination receives different amounts of products, it was assumed that the relative quantity of freight sent to each destination would follow a negative exponential function, as in Bartholdi & Gue (2004) where f j is the fraction of the total flow of freight that goes to destination j , f max and f min are the maximum and minimum values of f j respectively, c f is a constant, and m 2 is the number of destinations.Destinations are classified into a decreasing order of cargo flow magnitude.
Selecting appropriate values for f max , f min , and c f and applying function (1) to the total daily flow of freight, one has the average distribution of cargo per destination as depicted in Figure 2.
meaning that all inbound trailers have a capacity C T and travel fully loaded.On the other hand, the total daily freight load ∅ j designated to destination j , and containing cargo from all origins, is where G is the total quantity of freight received at the terminal during the period {t 0 , t E }, in a typical working day.Since the Monte Carlo algorithm works sequentially, extracting random carton packs from a virtual inventory to be loaded into a trailer, the resulting load configurations are not necessarily independent in statistical terms.Therefore, and in order to avoid biased flow arrangements in terms of destinations and cargo quantities, the resulting load vectors V k, j are randomly sorted according to variable k(k = 1, 2, . . ., n v ) at the end of the Monte Carlo routine.
During a typical working day, each destination j receives a total of ∅ j cartons, according to the proportions shown in (1) and in Figure 2. The Monte Carlo algorithm (Fig. 3) is composed of two stages.The first one is a greedy routine that tries to insert a larger number of carton packs into the vehicle, per operation.But, closer to the end of the loading process, the possible combinations of carton packs remaining in the virtual stock and ready to be loaded into a trailer, might be reduced to a fewer combinations.At this point, the search for an adequate random number of packs to fit the remaining space in the trailer might take too long a time.Thus, a complementary step-by-step loading routine will insert one carton pack at a time into the trailer, by varying the destinations until all the cargo is completely loaded.Thus, cargo arrangements that eventually remained incomplete at the end of the first stage, will be fully loaded at the end of the one-by-one unit stage.Figure 4 shows two Monte Carlo generated trailer load configurations, each trailer having a capacity of C T = 1, 400 cartons.Configuration (a), in Figure 4, depicts freight assigned to 9 destinations, a number about 39% of the total number of destinations.One can see a significant variation of cargo volumes per destination in Figure 4a, although the dispersion among destinations is not very expressive in this case.Conversely, there are less destinations in Figure 4b (7 destinations, about 30% of the total), but more dispersed than in the other configuration.

Door assignment
Most small cross docks are I-shaped because this type of layout offers the chance to move freight directly across the dock from a strip door to a stack door.But, in general, the assignment of strip and stack doors along the terminal perimeter is mostly based on getting the least average distance among all the doors on the dock (Bartholdi & Gue, 2004).It is easy to see by this measure that doors in the centre of an I-shaped dock are the most convenient, and this measure of convenience decreases from the centre doors to the doors at each end of the dock (Bartholdi & Gue, 2004).
In practice, an approximate evaluation criterion is to consider the door with the smallest average distance to all other doors as the "best door", because it is the one to/from which a freight lot can be expected to move internally with the least travel distance.According to Bartholdi & Gue (2004), the objective is to minimize the internal average travelled distance weighted by the freight volume The average distance d is assumed to be an approximation of the total internal travel cost across the dock, and therefore a proxy estimate of the variable labour cost to move freight through the facility (Bartholdi & Gue, 2000, 2004).The optimization problem (4) is a combinatorial problem with the following formulation: (a) from the total of m 1 + m 2 available dock doors, select m 1 doors to be used as strip doors according to some appropriate criterion, and m 2 doors to be stack doors; (b) assign each one of the m 2 stack doors to a specific destination.Strip doors, on the other hand, are assigned to inbound trailers along the working day according to the queue discipline and other factors, not having a specific identification related to freight origin, or other.
In the literature, this type of problem, with a number of variations, has been solved with a mathematical integer programming model, a genetic algorithm formulation, a simulated annealing approach, or other (Agustina et al., 2010; Van Belle et al., 2012).Depending on the total number of dock doors, the problem can be somewhat difficult to solve in mathematical terms.Bartholdi & Gue (2004) propose an approximate computational heuristic easier to be solved, and based on the measurement of average internal travel distance.In the Bartholdi & Gue (2004) formulation, one starts assigning strip doors to the "best doors", i.e. the doors having the smallest average distance to all other doors.Then, one assigns the destinations which receive the most freight to the next-best doors, and so successively.
The approach to treat our specific problem is similar to Bartholdi & Gue (2004) heuristic.We start the solution process by assuming a standard rectangular layout as in Figure 1, and assign an even number to the total number of dock doors m 1 + m 2 , thus allowing an equal number of dock doors at each side of the facility (Fig. 1).Next, one assigns strip doors to the "best doors", i.e. the doors having the smallest average distance to all other doors, as in Bartholdi & Gue (2004).
Considering the random nature of trailer arrivals, the randomness of the unloading times, and the varying occurrence of trailers waiting in the queue, as well as other probabilistic factors, there is no information beforehand on which strip door an arriving trailer will be assigned.Since all strip doors are already assigned at the beginning of the computation model, one can substitute equation ( 4) by taking the average distance d ( j ) from all the strip doors i = 1, 2, . . .m 1 to stack door j Now, we sort destinations j ( j = 1, 2, . . ., m 2 ) according to the decreasing order of the freight flow ∅ j ( j = 1, 2, . . ., m 2 ).We also sort stack doors according to the increasing order of the average distance d ( j ) ( j = 1, 2, . . ., m 2 ), given by ( 5).The destination with the largest value of ∅ j is assigned to the stack door with the minimum value of d ( j ) .Next, the second largest value of ∅ j in the list is assigned to the second value of d ( j ) in the corresponding list, and so successively.
With this procedure we assign all the m 2 destinations j to the corresponding best stack doors, thus minimizing (5).This process leads to the minimum value of d, as required in (4).This procedure has shown to be a simple, but a robust approximate heuristic for the assignment of dock doors.
A door assignment example with m 1 = 15 strip doors is shown in Figure 5. First, one assigns strip doors to the "best doors", i.e. the doors having the smallest average distance to all other doors, as in Bartholdi & Gue (2004).The number of stack doors in the example is m 2 = 23, but the volumes of freight to destinations 1 and 2 are significantly higher, requiring one additional door each.Each additional door will be located next to the basic one, forming pairs of twin stack doors.Thus, there will be m = m 1 + m 2 + 2 = 15 + 23 + 2 = 40 doors.Internal rectangular distances are measured on a system of coordinates with its origin at the southwestern edge of the facility (Fig. 5).Furthermore, the 40 doors are numbered in two layers along the south and the north edge of the terminal, as shown in Figure 5.As long as the doors positions are defined, a matrix with the distances between doors is easily computed.Compute now the average distances d ( j ) , ( j = 1, 2, . . ., m 2 ) according to (5), and sort the stack doors according to the increasing order of magnitude of d ( j ) , as shown in Table 1.Then, assign the door with the least value of d ( j ) , (door #9, with d (9) = 29.52m), to the destination j with the greatest value of ∅ j , i.e. ∅ 1 = 7, 740, and so successively, as shown in Table 1, columns (d) and (h) respectively.Notice that door 10, adjacent to door 9, was also assigned to destination 1, as explained, forming a pair of contiguous doors.The same occurs with stack doors 11 and 12, assigned jointly to destination #2.Table 1 shows the assignment of doors for this example, a routine repeated for the other alternative layouts with varying values of m 1 .

Trailer unload scheduling
Although Newell (1971) mentions some simple examples of queueing problems with disciplines different from the FCFS configuration, it is implicit in his method the adoption of the latter discipline.Besides, the adoption of a different queueing discipline in this work, to be described next, involves more complex mathematical formulation as compared to the FCFS alternative.Thus, a supplementary simulation model was developed in order to bypass such difficulty.For this, the Arena simulation software, Version 14.0, was used.In the simulation model, when an inbound trailer k(1, 2, . . ., n V ) arrives at the terminal it is assigned to a module identified as Reception.Then, a virtual badge is assigned to the trailer k, and it remains stuck on it until the vehicle is finally assigned to a strip door in order to be unloaded.Let us assume that the vehicle k has been ultimately assigned to strip door i.This vehicle brings a mix of freight represented by vector V k, j (Section 3.1).Each destination j is associated to a specific stack door, as seen in Section 3.2.Thus, the average weighted distance to transport freight V k, j from strip door i to all stack doors j ( j = 1, 2, . . ., m 2 ) is noticing that the simplification in the second hand side of ( 6) is due to (2).Therefore, the virtual badge assigned to the incoming vehicle k contains a vector π k,i that indicates the average weighted distance from every strip door i to all stack doors contemplated in vector V k, j .Thus, the new queue discipline is to put the arriving trailer into the waiting line, and when a generic strip door i becomes available, select the vehicle k in the queue (including an eventual arriving vehicle) that shows the least average distance π k,i in its virtual badge.Accordingly, the vehicle k will be assigned to the vacant strip door i.In Arena's terms, any incoming trailer is assigned to a sub-model identified as Queue, where it is assigned with the virtual badge π k,i , and waits its turn to be served.On the other hand, if more than one strip door is vacant, the first selected vehicle from the queue and the corresponding strip door will be such that the value of π k,i is minimum among all the candidates among k and i.In another words, the arrival order of the trailer will not influence the servicing order, the selecting criteria being only the expected average distance given by ( 6).This queueing discipline is named LDFS -Least-Distance, First-Served.Gue (1999) argues that this discipline could cause high waiting cost, with trailers stuck in the queue for a long period of time.The effect of this consideration will depend on the freight flow pattern, and on the size of the terminal, the latter reflected by the number of dock doors.Of course, such possible effect can be evaluated through the analysis of the overall cost to be developed in Section 4.3, together with a comparison with the traditional FCFS discipline.A mixed queue discipline, imposing an upper waiting time limit, could also be developed.In the Arena software, special queue discipline algorithms can be introduced into the simulation model as Visual Basic sub-programs.Such VB addition, adopting the LDFS discipline, has been performed in this work, in the trailer strip door scheduling module.

Inbound trailer arrival process
Let t 0 < τ 1 < τ 2 • • • < t E be the arrival times of the incoming trailers, where λ(τ ) is a continuous function representing the trailer arrival rate.Assume the vehicles are numbered according to the order of their arrival.Figure 6 depicts a typical diagram representing the truck arrival process at a cross-docking terminal, with arrival times restricted to the time window {t 0 , t E }.Such a diagram is based on Odoni & de Neufville (1992), who dealt with passenger arrivals at airport check-ins.Although it refers to the arrival of passengers at an airport counter, this scheme represents quite well the human behaviour in other operational cases, when clear servicing limitations, in the form of time windows, are observed.First, there is a minimum arrival time restriction in both cases.In the airport situation, the check-in process starts at an established time period before airplane take-off, and analogously, the truck unloading process starts at a pre-defined time t 0 (Fig. 6).On the other hand, passenger check-in terminates at a given instant governed by the airplane departure time.A similar situation occurs at a cross-docking terminal, with a definite deadline time t E for incoming trailers.Otherwise, if an arrival deadline were not imposed, it might not have enough time to unload the incoming trailers, process the cargo, and load the shipping vehicles in due time.According to Odoni & de Neufville (1992), the users (truck operators in this case) adopt, in general terms, the following behaviour pattern.Each incoming vehicle driver tries to choose an arrival time t such that, according to his/her perceptions, it guarantees a 95% probability (or greater) of respecting the terminal deadline.Thus, one is not considering the absolute value of t , but, instead, its perceived or anticipated value, which is strongly influenced by the driver's worse previous experiences, including road traffic conditions, adverse weather situations, etc., as well as the punitive measures imposed by the terminal administration to faults, whenever users do not observe the deadline restriction.
When the strip doors are initially opened to unload incoming trailers, the arrival rate is low, and λ(t ) increases until time t * , which is the point of maximum value of λ(t ) (Fig. 6).It is seen in Figure 6 that t * is not located at the middle point of the time window {t 0 , t E }.In fact, the curve of λ(t ) is skewed to the right.This situation can be explained by looking at the driver's behaviour, by which, when t > t * , his/her preoccupation in not missing the unloading operation increases substantially.Of course, the λ(t ) pattern depends on the application, but the curve depicted in Figure 6 represents reasonably well the truck driver's behaviour when approaching the cross docking terminal.
A third degree polynomial function was fitted to the curve depicted in Figure 6 λ(t ) = ω(at + bt 2 + ct 3 ), (7) where ω = 0, and a, b and c are constants.A third degree equation can be expressed as a function of its roots where τ A , τ B and τ C are the roots.Letting t 0 = 0 and, in accordance to Figure 6, λ(t 0 ) = 0, one makes τ A = 0 in (8).Another root of ( 8) is τ B = t E − t 0 = L, since λ(L) = 0 too.It can be shown that equation ( 8) has a double root at point Expanding (9), it results In order to determine the value of ω one imposes that the integral of λ(t ) in the interval {t 0 , t E }, must be equal to G, which is the total quantity of freight received at the terminal during the period {t 0 , t E }, in a typical working day leading to Substituting ( 12) in ( 10), one has One observes in Figure 6 that the curve λ(t ) shows a point of maximum in the interval {t 0 , t E } approximately at t * ∼ = 2 3 L.In order that t * be a maximum, the first order derivative of λ(t ) must be zero at that point Pesquisa Operacional, Vol.37(1), 2017 However, t * = 2 3 L does not satisfy (14), meaning that such point is not a maximum of ( 13).However, considering the point t * * = 1  3 L, which is symmetric to the point t * = 2 3 L when one takes the middle point of the {t 0 , t E } segment as reference, satisfies (14).But, in order that t * * be a point of maximum of ( 13), the second order derivative of λ(t ) must be negative at that point which is indeed satisfied by t * * = 1 3 L. Therefore, we make a scale transformation (L − t ) ← t , leading to Now, the point of maximum of λ(t ) in the interval 0 ≤ t ≤ L is located at t = 2 3 L, as expected.An important function in our analysis is the integral of ( 13), which will be used further in this paper to compute the average queue length Furthermore, the cumulative probability distribution of λ is As mentioned in Section 3.3, the Arena software allows for the introduction of input programs written in Visual Basic, whenever they are not available in the original package.This was necessary with regard to the trailer arrival process represented by the λ(t ), equation ( 16).Accordingly, a Visual Basic "Routine.Lambda" sub program was developed and introduced into the Arena software, whose pseudo-code is exhibited in Figure 7.

Inbound trailer queueing model
In the first part of this section the average queue length is determined using Newell's (1971) methodology, which uses continuous approximate queueing models with a FCFS discipline.The application of Newell's (1971) queueing approach to the FCFS configuration presupposes that the underlying system is subject to heavy congestion conditions, meaning in practice that the queue is not empty (Newell, 1971).Optimization is not the objective of the model but, instead, it is a hybrid discrete, continuous, approximate, simulation model.In addition, and using Little's Law (Little, 1961;Little & Graves, 2008), the corresponding average waiting times are also estimated.

Queue length
Since part of the incoming trailers do not find a vacant strip door upon their arrival at the facility, they may have to wait in order to be served.As mentioned in Section 2.1, classical queueing models generally do not apply to the analysis of queues formed by incoming trucks to be unloaded at cross-docking facilities.In fact, in most cross-docking terminals the queueing process does not reach a stationary state and, furthermore, the trailer arrival rate may vary substantially over the time window period.Thus, a different queue modelling method had to be used instead.
In our research, the formulation developed Newell (1971) has been employed.Newell's approach is intended to analyse non-stationary queues with approximate, straight-forward techniques.Whereas most of the literature on queues deals with equilibrium distributions of queue length and waiting time, the main emphasis of the Newell (1971) approach is on time-dependent behaviour as in this paper, where the arrival rate of users changes along time.The techniques employed by Newell (1971) are mainly fluid and diffusion approximations.Much of the classical theory of queues deals with discrete models, but in simple approximations such as the one adopted here, it is convenient to disregard the discrete nature of some elements and treat them as continuous, or quasi-continuous variables.The fluid approximation adopted in this queuing analysis assumes that the discrete, random arrivals of incoming trailers can be approximated by a deterministic quasi continuum, as it were a fluid flowing into a reservoir (Newell, 1971), and expressed here in cartons.The conversion between trailers to cartons, and vice-versa, is simply done in this application by considering that an incoming trailer always travels fully loaded with C T cartons.
Figure 8 shows a typical cumulative arrival rate A(t ), as represented in (17).Assuming a constant service rate μ per strip door, and m 1 strip doors in operation, the overall service rate is μm 1 .
Since the service rate is assumed constant, the cumulative service rate D(t ) is a straight line function shown in (Fig. 8).
where a is a constant, whose value will be analysed further in this section.Assuming no queue exists at time t 0 , i.e.Q(0) = 0, the arrival rate λ(t ) is less than μm 1 in the interval t 0 ≤ t ≤ t 1 , when the queue remains essentially nil (Fig. 8).At time t 1 the slope of A(t ) matches μm 1 , i.e. λ(t 1 ) = μm 1 , and the straight line D(t ) is tangent to A(t ) at that point (Newell, 1971).
Beyond t 1 , one has λ(t ) > μm 1 , and a queue starts to form, with the service rate varying linearly according to (19).The queue length, represented by Q(t ), is given by the difference (Newell, 1971) Figure 8 -Queue formation when t B > t E .
The process ends when D(t ) = G (Fig. 8), where G is the total flow of freight received at the terminal during the time interval {t 0 , t E }.The value of a in ( 19) is determined by means of a local search algorithm.Time t * is the point of maximum λ(t ) value in the interval t 0 ≤ t ≤ t E (Figs. 6 and 8), and thus t * is the inflexion point of curve A(t ).Therefore, the search for a is limited to the interval 0 ≤ t 1 ≤ t * , since otherwise the straight line D(t ) would not cross the curve A(t ), and consequently no queue would be formed.To numerically determine the value of t 1 , one solves the following search problem where λ(t ) is given by ( 16).The function Abs [λ(t ) − μm 1 ] is unimodal in the interval 0 ≤ t 1 ≤ t * , and problem ( 21) is solved by means of a golden section search technique (Gerald & Wheatley, 2004).Constant a is determined by putting t = t 1 in ( 19) Substituting ( 23) into ( 19), one has An important element in this queue analysis is the point t B where D(t B ) = G, i.e. when the last trailer in the queue starts its unloading process, as shown in Figure 8. Varying the number m 1 of strip doors, the D(t ) pattern changes accordingly.In Figure 8, t B > t E , but other situations may occur, as will be discussed further in this section.Putting t = t B and D(t B ) = G in (24), and solving for t B , one has Three queue cases are analysed next.

Case I: t B ≥ t E
As shown in Figure 8, depending on the value of m 1 , and for given values of the other parameters of the problem, a queue may still occur for t > t E .No more trailers are accepted to be unloaded when t > t E , but a group of remaining trailers are still in the queue, and must be served until the waiting line reaches zero.In this case, the queue is analysed in the time interval t 0 ≤ t ≤ t B , separated in three parts (Fig. 8) The average queue length, expressed now in number of vehicles, is (Fig. 8) Case II: t B < t E and t 1 < t * If t 1 < t * , the tangent D(t ) to the curve A(t ), at time t 1 (point A, Fig. 9), crosses curve A(t ) at point C, with abscissa t Z .Thus, the phase with queue formation comprises the time interval {t 1 , t Z }, as shown in Figure 9. Accordingly, the queue is thus defined Pesquisa Operacional, Vol.37(1), 2017 To determine t Z , a golden section search technique was employed in the interval t * ≤ t ≤ t E (Gerald & Wheatley, 2004), as follows f) The average queue length is (Fig. 9) Case III: t B < t E and t 1 < t * As long as the number of strip doors m 1 increases, it is observed that lim And, since t 1 converges to the inflexion point t * , the tangent D(t ) to curve A(t ) tends to be above D(t ) in the interval t * ≤ t ≤ t E , and therefore, in the limit, D(t ) will not cross A(t ).Consequently, in this specific case, Q(t ) → 0 as t 1 → t * .Nevertheless, Newell (1971) points out that the fluid analogy requires that the queue will remain large enough so that the cumulative number of departures might approximately be represented by a non-random (quasi) continuous variable.Otherwise, when Case III occurs with λ(t ) < μ m 1 , the fluid analogy should not be applied (Newell, 1971).In practical terms, when searching for a number of strip doors with minimum total cost, the optimal solution will normally be under Case I or II.

Trailer waiting time
Little's Law (Little, 1961;Little & Graves, 2008) states that, under steady state condition, the average number Q of items in a queue line equals the average arrival rate λ, multiplied by the average time W t that an item spends in the queue or in the system, depending on the application and on the analytic framework (Little & Graves, 2008).The law is This relationship is remarkably simple and general.Theoretically, Little's Law requires stationarity conditions about the underlying stochastic processes, but it is quite surprising that it can be employed more generally (Little & Graves, 2008).Little's Law provides a fundamental relationship among three key parameters in a queuing system.In fact, the system can be very general.An important feature of Little's Law is that, by knowing two of the three parameters, the third one can be calculated directly.This is a useful property since measurements of all three parameters may be difficult in certain applications.Little's Law presents great robustness and generality.
It is exact even in cases when arrival service processes are not stationary.The essential condition is to have a finite time window of observations that starts and stops when the system is empty (Little & Graves, 2008).This is the case of the queuing process studied in this paper.For this, let us consider Case I first (Section 3.5.1).The total number of units that arrive during the interval t 0 ≤ t ≤ t B is G (Fig. 9), letting λ = G/(t B − t 0 ).On the other hand, the average number of freight units Q in the queue is given by expression (29).Conversely, for Case II one has λ = G/(t E − t 0 ), with Q given by expression (33).

Model flow chart
The flow chart of the underlying model is presented in Figure 10.After the introduction of the input data (Section 4.2), the model starts defining the elements concerning the inbound trailer arrival process, analysed in Section 3.4 (block 1, Fig. 10).The trailer servicing process follows a log-normal distribution of the unloading time, with parameters defined in Table 2. Next, the inbound freight flows are generated through a Monte Carlo process, as described in Section 3.1 (block 2, Fig. 10).The basic decision variable of the model is the number m 1 of strip doors (block 3).The variation of m 1 was set by empirically evaluating the congestion level of the dock, leading to a minimum of 9, and a maximum of 19 strip doors.Next, in block 4, the queue discipline is contemplated.
For the FCFS discipline, the Newell approach was employed (blocks 5 to 10).The average queue length Q is estimated according to three alternative formulations, represented by Cases I, II, and III, and analysed in Section 3.5.1.On the other hand, the average waiting time W was estimated using Little's law, Section 3.5.3.
Costs are computed in block 15 (Fig. 10) following the methodology to be described in Section 4.3.Finally, the results of the work (block 16) are analysed in Section 4.4.Table 2 shows the main input data used in the analysis.

Incoming trailer queue simulation with LDFS discipline
The structure of the queuing model with LDFS discipline was described in section 3.3.The simulation was developed in Arena, version 14.0.In order to define the LDFS structure in detail, it was necessary to first schematize the layout of the terminal, setting its dimensions and door positions, and defining the distance matrices between doors for all the strip door alternatives (block 11, Fig. 10).Next, the strip doors are first assigned (block 12, Fig. 10), followed by the stack doors, as discussed in Section 3.2.Computations in blocks 11 and 12 (Fig. 10) were developed in Visual Basic and inserted in the Arena software.Blocks 13 and 14 correspond to Arena proper computations.Block 13 deals with the trailer scheduling process analysed in Section 3.3.Finally, the average queue length Q and the average waiting time W are computed by the Arena software (block 14, Fig. 10).
Since the number of stack doors is fixed and equal to 25 (Table 2), we adopted a varying odd number of strip doors in such that the total number of operating doors (strip plus stack) would be an even number, thus letting a fully symmetric door configuration on both sides of the terminal.
Accordingly, the number m 1 of strip doors covers the range from 9 to 19, in odd increments, as shown in Table 2.
The central variable of interest here is the mean number of incoming trailers in the queue.The number of replications in the simulation was set with the help of Arena's Output Analyzer, assuming a α = 95% confidence interval and adopting an intermediate configuration of m 1 = 15 strip doors.Tentatively assuming a number n where is the critical value of the t statistics.Approximating the t variable by the corresponding normal critical value z 1−α/2 , one has For α = 95%, z 1−α/2 = 1.96, and putting s = 0.896 and h = 0, 255 in (37), one gets n (2)  r = 47.Thus, the round number of replications n r = 50 attends well the 95% confidence interval requirement.
The hardware to process the simulation has an Intel i5 3570 processor, with an 8 GB RAM, and SSD 500gb Samsung 850 Evo.The average time to run 50 replications of each configuration, represented by the number m 1 of strip doors, was 119.3 seconds, ranging from 91 to 135 seconds.Average queue length, queue length standard deviation, and trailer waiting time, both for the FCFS and the FDFS disciplines, and obtained via simulation, are presented in Table 3 and depicted in Figure 11.

Cost analysis
The objective of the comprehensive model is to minimize the sum of three costs: (a) the rent cost of the terminal building facility, which is a function of its area and layout; (b) the waiting cost of inbound trailers when being attended at the terminal, and (c) the cost of handling and moving freight inside the terminal, which depends upon the material handling equipment in use and upon the facility layout.The costs related to outbound trucks are fixed in our problem, since the assignment of those trucks to stack doors is done externally to the model, as indicated in Section 2.2, and therefore they do not influence the search for the minimum cost.

Economic rent of the facility
Let m be the total number of dock doors located on both sides of the rectangular facility.Usually, the number of stack doors is greater than the number of destinations m 2 , because some destina-  tions may show a larger freight flow level, requiring more than one stack door for shipping the cargo (Bartholdi & Gue, 2004).In this application, destinations 1 and 2 have shown larger volumes (see Fig. 2), thus requiring one additional stack door each, as to handle excess freight to these destinations.So, the total number of dock doors is m = m 1 + m 2 + 2 in this application.Since m was assumed in the model to be an even number, there will be m/2 doors on each side of the building.An additional area of 400 m 2 was assigned for administrative activities.Therefore, the facility area will be in general The operating cost of the terminal was based on a similar facility located in Tucson, AZ, USA (Boyd Co., 2013), with an approximate annual cost of US$ 266.70/year/sq.meter,which con-verted to the BRL Brazilian currency (at an exchange rate of 3.28 BRL$ to the US dollar) lead to BRL$ 874.73/year/sq.meter.On the other hand, inbound vehicles require a parking space when waiting to be unloaded, circa 120 sq.meter per trailer.A margin of 20% has been added to cover circulation, control space, and eventualities, leading to an estimated area of 144 sq.meter of parking space per trailer.In its turn, the land to be reserved as a parking lot for the inbound trailers was assumed to be worth BRL$ 400.00 per m 2 , to be amortized in 15 years, with a yield annual rate of 8.5%.
Table 4 shows the total annual rent cost of the terminal.Since the total freight flow through the terminal is G = 84, 000 cartons per working day, the last column in Table 4 shows the average terminal rent cost per carton.It was assumed 300 working days per year.

Inbound trailer waiting cost
The estimated trailer waiting times for queue disciplines FCFS and LDFS are exhibited in Table 3 and Figure 11.Analysing empirical data related to Brazilian freight transport, an estimated average rental cost of C trl = BRL$ 40.00 per hour, per trailer was defined.Thus, the total average waiting cost of incoming trailers is where W t is the average inbound trailer waiting cost per carton, for queue discipline δ, n V is the number of inbound trailers per day; C trl is the average trailer rental cost per hour; W (∂) is the average trailer waiting time for queue discipline δ, and G = 84, 000 cartons is the total daily freight flow through the terminal.Table 5 presents the values of C (δ) W t as a function of the number of strip doors and queue discipline.

Cost for handling freight inside the terminal
Since the cross dock terminal studied in this paper is a small one by international standards, and since in developing countries there is quite a number of terminals still manually operated, we decided to examine two situations.The first one, using manually operated four-wheel platform dollies, and the second, employing electric pallet trucks.

A -Freight handling with four-wheel platform dollies
The four-wheel platform truck is 0.81 m long and 0.73 m wide, carrying loads up to 1,800 kg and taking about 20 cartons per voyage.A worker takes circa 6 minutes to arrange 20 cartons on a dolly, with the opposite operation taking approximately the same time.Thus, the effort to load and unload the dollies, for a total flow of 84,000 cartons per day, will require a total of 840 manhours/day.Assuming a speed of 3 meter/sec for a worker moving a dolly across the facility, with an average travelled distance of d meters between the strip and the stack door, the estimated dolly cycle time is t Cycle = 2d/3 seconds.Since 20 cartons are transported in each cycle, the number of cycles per day will be 84,000/20 = 4,200 cycles.Thus, the daily labour effort to displace the freight across the facility is 0.778 man-hours/day/meter.
Introducing two correcting factors, the total daily labour effort of handling and moving freight across the terminal, with four-wheel dollies, is where H d is the total labour effort in man-hours per working day, d is the average internal distance between strip and stack doors, in meters as estimated in Section 3.3, α = 0.3 is a margin to cover unexpected demand variations, and = 0.7 is the estimated labour efficiency.
Personnel assigned to this kind of job receives an estimate of 1.9 minimum legal wages, or BRL$

B -Freight handling with electric pallet trucks
The electric pallet truck is a Hyster Spectrum P20S, with a 112.5 cm long fork and about 2,000 kg load capacity.General operating data are as follows • Travelling speed: 220 ft/min = 1.12 m/sec; • Pallet length 1.20 m, width 1.00 m; • Cartons per pallet: 50; • Time to reach a pallet and lift it: 3 sec; • Time to lower a pallet and liberate it: 3 sec; • Time margin related to delays, traffic idle time, etc.: 30%; • Overall operating efficiency: 70%.
Thus, the pallet truck cycle time is where t trk is the pallet truck cycle time in minutes, α = 0.30 is the time margin, d is the average rectangular distance from the strip door to the stack door in meters, and = 0.7 is the overall operating efficiency.Since a pallet truck operates T = 8 hours = 480 minutes per day, the number of cycles it will perform per day is The total expected number of pallet truck cycles per day n T is equal to the total daily freight flow G = 84, 000 cartons divided by the number of cartons per pallet, or 84, 000 ÷ 50 = 1, 680 cycles/day.Since the pallet truck operates T = 8 hours = 480 minutes per day, the required number of trucks is which includes an excess of 10% to take care for out-of-order trucks, periodic maintenance, etc.Since freight does not arrive palletized, the cargo must be palletized just upon the unloading operation at the terminal.It requires 5 man-hours to palletize 1,000 cartons.Since a total of 84,000 cartons are handled per day, an effort of 420 of man-hours/day is necessary to palletize all the freight.For a total operation period of T = 8 hours per day, it will require 420/8 ∼ = 52, 2 men.Adding a 10% margin, it leads to a total of 58 workers to perform the palletization task.The labour cost related to this kind of service is the same adopted in Section 4.2.3A,i.e.BRL$ 122.70 per day, per man, or BRL$ 7,116.60 per working day.
Table 7 presents the costs for handling and moving freight inside the terminal for the electric pallet truck configuration.

Modelling results
In this section we present the results of a numerical experiment carried out to evaluate the performance of a cross-docking terminal (small by international standards), employed by a supermarket organization in the distribution of freight to local stores.Since the main objective of the study has been to investigate the effects of different terminal settings and formulations, we have opted to generate some empirical data to test the model, as explained in Section 3.1.
The application of the methodology described in this research work has led to the following approximate characteristics of the analysed cross docking terminal: • A total of 40 dock doors is proposed, distributed along both sides of the I-shaped facility (see Fig. 5); • A total of 15 strip doors; • A total of 23 destinations, with 25 stack doors, in which destinations 1 and 2 have two stack doors each; • Stack doors layout and assignment as in indicated in Figure 5; • Inbound vehicle unload scheduling is performed according to a LDFS queue discipline; • Inbound trailers wait an average of 51 minutes to be unloaded (Table 3); • The expected queue size of inbound trailers is 1.2 vehicles (Table 3); • An internal terminal area of about 2,992 square meters is necessary; • Use electric pallet trucks to move freight inside the terminal; • Total expected cost of BRL$ 0.253 per carton.
In Table 8 and Figure 12 it is shown the total cost per carton as a function of the number of strip doors, and considering the two alternative queue disciplines FCFS and LDFS.The two alternative freight handling alternatives are also analysed.In Figure 11 one can see that the queue length for the LDFS discipline is smaller than in the FCFS alternative, because the trailer scheduling process is better organized in the latter discipline, but the waiting time comparatively increases, as expected when compared to the classic discipline FCFS.
For the four-wheel dollies, the best solution is to assign 15 strip doors associated with the FCFS queueing discipline.When adopting the electric pallet trucks, the best solution is also 15 strip doors, but with the LDFS discipline.However, although the analysed cross-docking terminal is seen as a small facility by international terms, the unit cost of the pallet truck version is much lower than the manual solution.In fact, the four-wheel dolly best alternative costs about BRL$ 0.41 per carton, while the pallet truck version costs about BRL$ 0.25 per carton, a value about 39% smaller.
Table 8 and Figure 12 show that the employment of pallet trucks to carry freight in the terminal shows a much smaller cost than the four-wheeled platform dolly alternative.But, with regard to the use of either the FCFS or the LDFS queuing disciplines, Table 8 and Figure 12 show that both disciplines present approximately the same cost in this application.Apart from the factors discussed in Section 3.1, the advantage of employing the LDFS discipline increases with terminal size.In fact, the configuration analysed here is related to a small facility.As long as the number of dock doors increases to adjust to the growing demand, the average distance between strip doors and stack doors increases accordingly, and the cost of internally moving freight tends to grow significantly.Thus, if the terminal freight movement tends to grow in the future, the facility will probably be enlarged at some point of time, and demanding for a LDFS orientation.So, an earlier adoption of such discipline can be tried out by the workers and management, letting it to be adapted to the company needs over a somewhat longer period of time, with positive benefits.So, the LDFS discipline seems to be better than the FCFS one.

CONCLUSIONS AND RESEARCH PROSPECTS
The methodological structure developed in this paper involves a number of analytic approaches intended to produce elements to help in the design of a cross docking terminal.Of course, the methodology is not universal, since different terminal problem settings will require adaptations and further developments in some specific parts of the model.But one important and positive result of this work is that some parts of the problem, normally considered cumbersome in terms of mathematical treatment as, for example, queueing analysis, door assignment, and truck scheduling, can be treated in a robust and simple modelling way, with reasonable results, as shown in the present paper.
Some aspects of the problem will receive special attention in our future research efforts.First, the time window of 8 hours set up to receive freight at the terminal can be split into shorter time-windows, establishing, for example, four periods of two hours each.Such a measure would help in improving arrival reliability, and also leading to less strip doors and smaller facility areas, with reduced costs, but it would require more integration among the terminal organization, the suppliers, and the transport operators.
Regarding the use of simulation to solve this type of problem, the available commercial software, as Arena and others, provide simple evaluations such as in this application, in which the number of split doors constitutes the basic decision variable.But simple evaluations may not provide enough details for optimal decision making.One alternative is to combine simulation with optimization techniques to get better performance evaluations.For example, the Arena optimization module OptQuest allows to model and experiment with several alternative scenarios, letting the researcher to select the one that best meets the problem objectives.Thus, one research development to follow further will be to enlarge and detail the simulation framework of the problem, employing OptQuest-likely optimization search technologies in order to exploit better solutions.
Another point to investigate further is to focus on the benefits of anticipating the labelling and addressing of freight units to vendor premises, which will reduce freight handling and facility space at the terminal, with positive impacts in cost reduction.This, of course, will depend on specific arrangements with suppliers.Another research possibility is to introduce the congestion effects in terminal performance and costs, due to different types of internal traffic conflicts, as reported by Bartholdi & Gue (2000).

Figure 1 -
Figure 1 -A typical rectangular cross dock terminal.
In general, three extensively cited papers in the literature present a good review of cross docking research: Agustina et al. (2010), Boysen & Fliedner (2010), and Van Belle et al. (2012).

Figure 2 -
Figure 2 -Distribution of expected freight flow per destination.

Figure 3 -
Figure 3 -Pseudocode of the Monte Carlo routine to generate random freight loads for inbound trailers.

Following
Barthold & Gue (2004), let n V = number of inbound trailers per working day; m 1 = number of strip doors; m 2 = number of stack doors (equal to the number of destinations); i = a strip door, i = 1, 2, . . ., m 1 ; j = the destinations to which the freight loads V k, j , with k = 1, 2, . . ., n V and j = 1, 2, . . ., m 2 , are assigned to; d i, j = the rectangular distance from strip door i to stack door j ; C T = cargo capacity of one inbound trailer (in cartons); ∅ j = total quantity of freight sent to destination j , in cartons per working day (3).

Figure 7 -
Figure 7 -Pseudo-code to generate a sequence of inbound trailer arrival times.

( 1 )
r = 50 of replications, the Arena Output Analyzer yielded the half-width confidence interval h = 0, 255, and the standard deviation of the sample s = 0.896.A recursive approximate formula (Kelton et al., 2014; Rubisntein & Kroese, 2017) was used to further estimate the value of n r , based on the first guess

Figure 11 -
Figure 11 -Inbound trailer queue length and waiting time.

Figure 5 (
Figure 5 (Section 3.2) shows the internal arrangement of the dock doors as well as the basic terminal dimensions for an example with m 1 = 15, m 2 = 23 and m = 25.The door width is a = 3.60 m, and the internal width of the facility is b = 36 m.Since m was assumed in the model to be an even number, there will be m/2 doors on each side of the building.An additional area of 400 m 2 was assigned for administrative activities.Therefore, the facility area will be in general

Figure 12 -
Figure 12 -Cost per carton as a function of the number of strip doors.

Table 2 -
Main input data to the model.

Table 4 -
Facility rent cost as a function of the number of strip doors.

Table 5 -
Inbound trailer waiting cost as a function of the number of strip doors.
1,680.00 per month, plus legal obligations of 83.3%, leading to a gross wage of BRL$ 3,062.80 per month.Considering a total of 25 working days per month, a worker will get BRL$ 122.50 per working day, or BRL$ 15.30 per hour of work.A dolly acquisition cost is BRL$ 600.00, with an estimated life of 3 years.Considering 300 working days per year, one has a dolly cost of BRL$ 0.73/day.About 200 dollies were estimated to attend the internal transport service, leading to a subtotal cost of BRL$ 146.00/day.The cost of handling freight per day for this version is depicted in Table6.

Table 6 -
Cost for moving freight inside the terminal with four-wheel platform dollies.
The annual rent of a pallet truck is BRL$ 84,820.60,or BRL$ 282.73 per working day.The basic wage of a pallet truck operator is BRL$ 2,361.70 per month.Adding legal obligations of 83.3%, one has a total wage value of BRL$ 4,329.00 per worker and per month, or BRL$ 173.16 per working day.

Table 7 -
Cost for moving freight inside the terminal with electric pallet trucks.

Table 8 -
Expected cost per carton.