simultaneously the best sites of both plants and warehouses and the best strategy for distributing the prod-
uct from the plants to the warehouses and from the warehouses to the customers.
A common objective in designing such a distribution network is to determine the least cost system design
such that the demands of all customers are satisfied without exceeding the capacities of the warehouses and
plants. This usually involves making tradeoffs inherent among the cost components of the system that in-
clude: (1) costs of opening and operating the plants and warehouses, and (2) the inbound and outbound
transportation costs.
Many researchers have extensively studied facility and demand allocation problems. Previous research
studies are well surveyed by Francis et al. [5], Aikens [1], Brandeau and Chiu [4], and Avella et al. [2]. More
recently, Jayaraman [8] studied the capacitated warehouse location problem that involves locating a given
number of warehouses to satisfy customer demands for different products. Pirkul and Jayaraman [9] extended
the previous problem by considering locating also a given number of plants. They formulated the problem as
a mixed integer model and developed a Lagrangean based heuristic solution procedure. The procedure was
tested using problem instances with up to 100 customers, 20 potential warehouses and 10 potential plants.
Tragantalerngsak et al. [11] considered a two-echelon facility location problem in which the facilities in
the first echelon are uncapacitated and the facilities in the second echelon are capacitated. The goal is to
determine the number and locations of facilities in both echelons in order to satisfy customer demand of
the product. They developed a Lagrangean relaxation based branch and bound algorithm to solve the prob-
lem and reported results of computational tests with up to 100 customers and 15 facilities. Gourdin et al. [7]
studied a particular type of the uncapacitated facility location problem where two customers allocated to
the same facility are matched. They developed several methods to solve the problem after deriving valid
inequalities, and optimality cuts for the problem.
One major drawback in most of past research studies like [7–9,11] is that they limit the number of capac-
ity levels available to each facility to just one. However, as it is the case in practice, there exist usually sev-
eral capacity levels to choose from for each facility. The use of different capacity levels makes the problem
more realistic and, at the same time, more complex to solve. Another major drawback in some previous
studies like [8,9] is that they limit the number of facilities to open to a pre-specified value. Moreover, these
studies fail to describe how this value can be determined in advance.
Our current study represents a significant improvement over past research by presenting a unified model
of the problem that includes the numbers, locations, and capacities of both warehouses and plants as var-
iables to be determined in the model and develops at the same time the best strategy for distributing the
product from the plants to the warehouses and from the warehouses to the customers. An efficient heuristic
solution procedure based on Lagrangean relaxation of the problem is developed and extensive computa-
tional tests with up to 500 customers, 30 potential warehouses, and 20 potential plants are reported.
The remainder of this paper is organized as follows. In Section 2, a mathematical formulation of the dis-
tribution design problem is presented. A Lagrangean relaxation of the problem is proposed in Section 3. An
efficient solution algorithm based on this relaxation is developed in Section 4. Computational results are
reported in Section 5. A summary of the work presented in this paper is given in Section 6.
2. Model formulation
The following notation is used in the formulation of the model.
Nindex set of customers/customer zones
Mindex set of potential warehouse sites
Lindex set of potential plant sites
Rindex set of capacity levels available to the potential warehouses
568 A. Amiri / European Journal of Operational Research 171 (2006) 567–576