# Order Allocation Research of Logistics Service Supply Chain with Mass Customization Logistics Service

Authors: Weihua Liu; Haitao Xu; Xinyu Sun; Yi Yang; Yuming Mo
Journal: Mathematical Problems in Engineering
 (2013)

---

## Abstract

This paper studies the order allocation between a logistics service integrator (LSI) and multiple functional logistics service providers (FLSPs) with MCLS. To maximize the satisfaction of FLSPs, minimize the total cost of LSI, and maximize the customized degree, this paper establishes a multiobjective order allocation model of LSSC that is constrained by meeting customer demand, customer order decoupling point, and order difference tolerance coefficient. Numerical analysis is performed with Lingo 12 software. This paper also discusses the influences of scale effect coefficient, order difference tolerance coefficient, and relationship cost coefficient on the comprehensive order allocation performance of the LSSC. Results show that LSI prefers FLSPs with better scale effect coefficients and does not need to set an extremely high order difference tolerance coefficient. Similarly, setting a high relationship cost coefficient does not necessarily correspond to better results. For FLSPs, the continuous improvement of large-scale operational capacity is required. When the comprehensive order allocation performance of the LSSC is high, the LSI should offer cost compensation to improve the satisfaction of the LSSC.


---

## Body

## 1. Introduction

Outsourcing service, such as information technology services, financial services, and logistics services, has become increasingly popular in recent years due to the fact it appears to be more profitable. Many service integrators and functional service providers develop a long-term service purchasing/supplier relationships to satisfy customers demand and obtain more profit with integrated services and then form a service supply chains (SSC). Specially, in a logistics service supply chain (LSSC), a logistics service integrator (LSI) provides customized logistics services by integrating the service capacities of multiple functional logistics service providers (FLSPs) [1]. For instance, Baogong Logistics Company, one of the largest LSI in China, integrates more than 500 storage companies, more than 1,200 highway transportation companies, and over 500 manual loading and unloading companies as its FLSPs and then accomplishes the integrated logistics services for its customers such as Procter & Gamble, and Philips.In recent years, customer requirements for specialized and customized logistics services have increased. Many logistics enterprises have also begun considering a change in logistics service mode. These logistics enterprises have attempted to provide mass customization logistics service (MCLS) instead of only mass logistics service [2]. MCLS is one of the most important trends in logistics service. The capability of offering MCLS to customers has become the key to improving market competitiveness. In the MCLS environment, many logistics enterprises spontaneously form a logistics service supply chain (LSSC) through unions and integrations to meet customized demand and offer large-scale services [1, 3]. As the core enterprise of LSSC, the logistics service integrator (LSI) integrates the advantages of other functional logistics service providers (FLSPs) in different processes and different functions and delivers the integrated service to customer. For instance, on 11 November, 2011, more than 2.67 million packages were delivered all over the country by the Yuantong Express Company in China; this number is more than four times the number of packages delivered on the same day in 2010. These 2.67 million packages were delivered not only within 3 days to 5 days to customers in 31 provinces in the country but also under customized logistics service. This process brought huge pressure to the Yuantong Express Company. To complete the courier service, Yuantong Express Company cooperated with other FLSPs (e.g., SF Express and ZTO Express) by outsourcing some customer orders.Under the MCLS environment, the LSI integrates different customer orders and allocates such orders to multiple FLSPs. In contrast to existing order allocation methods, the LSI needs to consider the influence of customer order decoupling points (CODPs). Each customer order has a unique CODP; therefore, when the LSI integrates multiple customer orders, the cost scale effect caused by large-scale services and the different CODPs of these orders should be considered. We also need to consider the minimum total cost of LSI and FLSP satisfaction in order allocation processes. Hence, a multiobjective decision-making problem exists. This problem involves addressing the diversity of customer orders, minimizing the total cost of the LSI, and considering the optimal satisfaction of FLSPs. This paper studies the order allocation problem (OAP) between LSI and multi-FLSPs, considers the characteristics of the MCLS by modeling and analysis, and explores the influence of related factors on the comprehensive performance of order allocation. Moreover, several management insights are provided.This paper is organized as follows. Section2 presents the literature review, where existing research on mass customization and order allocation are systematically summarized. Section 3 demonstrates the model building. An order allocation model of LSSC under the MCLS environment is established. Section 4 characterizes the model solution. A method for solving the multiobjective programming model is provided. Section 5 presents the numerical analysis. The influence of related factors on the comprehensive performance of order allocation is discussed in this section. Section 6 provides the main conclusions, management insights, and implications for LSI and FLSPs. Section 7 presents the limitations of the study and discusses suggestions for further research in this field.


## 2. Literature Review

This paper involves two main research areas: mass customization and order allocation. The literature review focuses on these two areas and elaborates their research progress and inadequacies.


### 2.1. Mass Customization

Studies on mass customization mainly include mass customization production and mass customization service. Mass customization production is primarily for manufacturing industry and mass customization service is mainly for service industry.


#### 2.1.1. Mass Customization Production

Most mass customization production studies have focused on the supply chain of the manufacturing industry [4, 5]. The main content of these studies is centered on the design of a delay system. Reference [6] proposed five kinds of postponement strategy: tags delay, packaging delay, assembly delay, manufacturing delay, and time delay. Thereafter, scholars studied the problem of choosing the appropriate postponement strategy in terms of the applicable conditions of time delay and cost delay in manufacturing [7]. Moreover, several scholars examined the influence of implementing postponement strategies on the cost of delay system [8, 9].Studies on the decision-making problem of CODP have increased in recent years. In enterprise manufacturing activities, CODP is the transition point from make-to-stock to make-to-order [10]. In our paper, CODP indicates the procedure to start customization service. The positioning model of CODP has progressed from single CODP to multi-CODP [11, 12] and from static CODP to dynamic CODP [13, 14]. Given the increasing number of studies on supply chains, several scholars have investigated that the CODP problem undersupply chain environments and achieved fruitful research results. For example, reference [15] studied the influence of CODP selection on inventory costs when the lead time and universal module average holding cost are the same regardless of CODP position, which considers inventory and other costs in the model. Reference [16] examined the best CODP selection model around the supply chain environment.


#### 2.1.2. Mass Customization Service

The increasing adoption of the service supply chain in recent years has attracted scholars in studying the CODP problem of the service industry under mass customization environments. These previous studies mainly examined how a single service enterprise positions the CODP in service processes while considering mass customization costs and ignoring inventory costs. Moreover, there is a special issue on service optimization and control in mathematical problems in engineering which includes a hot topic of mass customization service. There are also some scholars dedicating research perspective to specific industries, such as logistics services industry. References [17, 18] explored the scheduling model of LSSC in the context of MCLS.


### 2.2. Order Allocation on Supply Chain

To date, most OAPs are solved by using multiobjective linear programming or nonlinear programming methods, which spread around the produce supply chain. References [19, 20] studied the application of genetic algorithm for optimizing order allocation. Current studies on customer order allocation mainly consider the service level and minimizing procurement costs [13, 19, 21]. And most studies have focused on the manufacturing industry supply chain. Several scholars have investigated the integrated multiobjective decision-making method, which considers both supplier selection and OAPs [22]. In addition to the solving method, several scholars have focused on the order allocation model, such as research on the order allocation model of supply chains with FLSP load cases [17, 18, 23].The growth of LSSC has directed considerable research attention on its order allocation. Several studies have examined the OAP in two-echelon LSSC [12, 19]. Reference [1] investigated an emergency order allocation model based on multiproviders in two-echelon LSSC. Considering the multi-echelon nature of LSSC, references [17, 18] conducted a study on the order allocation of three-echelon LSSC.The literature review indicates that the MCLS environment has been left unexamined even though scholars have focused on the OAP of LSSC. The MCLS needs to consider the cost advantages of mass service and customized special requirements; therefore, the order allocation in this case is clearly distinct from general cases. Furthermore, cost control is not the most important decision-making objective when we consider the customization requirements. Thus, this paper focuses on the OAP around the MCLS environment by modeling its characteristics and exploring the influence of related factors on the comprehensive performance of order allocation. Furthermore, several management insights are also provided.


## 2.1. Mass Customization

Studies on mass customization mainly include mass customization production and mass customization service. Mass customization production is primarily for manufacturing industry and mass customization service is mainly for service industry.


### 2.1.1. Mass Customization Production

Most mass customization production studies have focused on the supply chain of the manufacturing industry [4, 5]. The main content of these studies is centered on the design of a delay system. Reference [6] proposed five kinds of postponement strategy: tags delay, packaging delay, assembly delay, manufacturing delay, and time delay. Thereafter, scholars studied the problem of choosing the appropriate postponement strategy in terms of the applicable conditions of time delay and cost delay in manufacturing [7]. Moreover, several scholars examined the influence of implementing postponement strategies on the cost of delay system [8, 9].Studies on the decision-making problem of CODP have increased in recent years. In enterprise manufacturing activities, CODP is the transition point from make-to-stock to make-to-order [10]. In our paper, CODP indicates the procedure to start customization service. The positioning model of CODP has progressed from single CODP to multi-CODP [11, 12] and from static CODP to dynamic CODP [13, 14]. Given the increasing number of studies on supply chains, several scholars have investigated that the CODP problem undersupply chain environments and achieved fruitful research results. For example, reference [15] studied the influence of CODP selection on inventory costs when the lead time and universal module average holding cost are the same regardless of CODP position, which considers inventory and other costs in the model. Reference [16] examined the best CODP selection model around the supply chain environment.


### 2.1.2. Mass Customization Service

The increasing adoption of the service supply chain in recent years has attracted scholars in studying the CODP problem of the service industry under mass customization environments. These previous studies mainly examined how a single service enterprise positions the CODP in service processes while considering mass customization costs and ignoring inventory costs. Moreover, there is a special issue on service optimization and control in mathematical problems in engineering which includes a hot topic of mass customization service. There are also some scholars dedicating research perspective to specific industries, such as logistics services industry. References [17, 18] explored the scheduling model of LSSC in the context of MCLS.


## 2.1.1. Mass Customization Production

Most mass customization production studies have focused on the supply chain of the manufacturing industry [4, 5]. The main content of these studies is centered on the design of a delay system. Reference [6] proposed five kinds of postponement strategy: tags delay, packaging delay, assembly delay, manufacturing delay, and time delay. Thereafter, scholars studied the problem of choosing the appropriate postponement strategy in terms of the applicable conditions of time delay and cost delay in manufacturing [7]. Moreover, several scholars examined the influence of implementing postponement strategies on the cost of delay system [8, 9].Studies on the decision-making problem of CODP have increased in recent years. In enterprise manufacturing activities, CODP is the transition point from make-to-stock to make-to-order [10]. In our paper, CODP indicates the procedure to start customization service. The positioning model of CODP has progressed from single CODP to multi-CODP [11, 12] and from static CODP to dynamic CODP [13, 14]. Given the increasing number of studies on supply chains, several scholars have investigated that the CODP problem undersupply chain environments and achieved fruitful research results. For example, reference [15] studied the influence of CODP selection on inventory costs when the lead time and universal module average holding cost are the same regardless of CODP position, which considers inventory and other costs in the model. Reference [16] examined the best CODP selection model around the supply chain environment.


## 2.1.2. Mass Customization Service

The increasing adoption of the service supply chain in recent years has attracted scholars in studying the CODP problem of the service industry under mass customization environments. These previous studies mainly examined how a single service enterprise positions the CODP in service processes while considering mass customization costs and ignoring inventory costs. Moreover, there is a special issue on service optimization and control in mathematical problems in engineering which includes a hot topic of mass customization service. There are also some scholars dedicating research perspective to specific industries, such as logistics services industry. References [17, 18] explored the scheduling model of LSSC in the context of MCLS.


## 2.2. Order Allocation on Supply Chain

To date, most OAPs are solved by using multiobjective linear programming or nonlinear programming methods, which spread around the produce supply chain. References [19, 20] studied the application of genetic algorithm for optimizing order allocation. Current studies on customer order allocation mainly consider the service level and minimizing procurement costs [13, 19, 21]. And most studies have focused on the manufacturing industry supply chain. Several scholars have investigated the integrated multiobjective decision-making method, which considers both supplier selection and OAPs [22]. In addition to the solving method, several scholars have focused on the order allocation model, such as research on the order allocation model of supply chains with FLSP load cases [17, 18, 23].The growth of LSSC has directed considerable research attention on its order allocation. Several studies have examined the OAP in two-echelon LSSC [12, 19]. Reference [1] investigated an emergency order allocation model based on multiproviders in two-echelon LSSC. Considering the multi-echelon nature of LSSC, references [17, 18] conducted a study on the order allocation of three-echelon LSSC.The literature review indicates that the MCLS environment has been left unexamined even though scholars have focused on the OAP of LSSC. The MCLS needs to consider the cost advantages of mass service and customized special requirements; therefore, the order allocation in this case is clearly distinct from general cases. Furthermore, cost control is not the most important decision-making objective when we consider the customization requirements. Thus, this paper focuses on the OAP around the MCLS environment by modeling its characteristics and exploring the influence of related factors on the comprehensive performance of order allocation. Furthermore, several management insights are also provided.


## 3. Model Building

Section3 introduces the model assumptions and variables. To maximize the satisfaction of FLSPs, minimize the total cost of the LSI, and maximize the customized degree of customers, we established a multiobjective order allocation model of LSSC that is constrained by meeting the demand, CODP, and order difference tolerance coefficient (Section 3.2).


### 3.1. Model Description

We considered a two-echelon LSSC, which consists of one LSI and many FLSPs. After receiving orders from multiple customers, the LSI analyzes the customized and universal service requirements of these customer orders and allocates them to more than one FLSP. Every logistic service for customer orders comprises multiple service processes. Each process needs one kind of service that can be provided by many cooperative FLSPs. The service capability of each FLSP is possibly different. We assumed that[aijk,bijk] denotes the service ability interval provided by FLSP ith for customer jth in the procedure kth; xijk stands for the jth customer’s order assigned to the ith FLSP in the kth procedure.Table1 shows the notations of the model. The other assumptions of the model are as follows.(1)
Each FLSPi  (i=1,2,…,n) can provide mass service and customization service for the LSI. The total service procedures of customer j are Lj  (j=1,2,…,M), and the order demand is Dj.(2)
Kj represents the procedure to start the customization service of the jth customer; this service is not necessarily the same for all customers. The LSI needs to consider the Kj of all customers and set a common CODP, which is signified by K′ for all customers. K′ is mass service, and K′ is customization service. K′ should be set before K0  (K0=min{K1,K2,…,KM}) because the CODP must be set before all Kj. Given the characteristics of the MCLS, Procedure 1 must be mass service; hence, 2≤K′≤K0).(3)
Given the importance of different customizations for each customer, the weight of thejth customer is βj (∑j=1Mβj=1).(4)
Customer orders are assigned to FLSPs through the LSI. FLSPs will accomplish the order assigned by the LSI. No game exists between LSI and FLSPs.Table 1
Model notations.
Notations
Description
[
a
i
j
k
,
b
i
j
k
]
The service ability interval of theith FLSP for jth customer in the kth procedure
[
a
i
j
k
1
,
b
i
j
k
1
]
The value of[aijk,bijk]  when the kth procedure is mass service
[
a
i
j
k
2
,
b
i
j
k
2
]
The value of[aijk,bijk] when the kth procedure is customization service
x
i
j
k
Thejth customer’s order assigned to the ith FLSP in the kth procedure
τ
Scale effect coefficient
D
j
The demand of thejth customer
L
j
The total service procedures of thejth customer
β
j
The weight of thejth customer
ρ
i
j
k
The satisfaction of theith FLSP for the order of the jth customer in the kth procedure
ρ
i
j
k
0
The initial satisfaction of theith FLSP in the kth procedure when the order is equal to the minimum service capacity of  ith FLSP
ρ
i
j
k
1
0
The value ofρijk0 when the kth procedure is mass service
ρ
i
j
k
2
0
The value ofρijk0 when the kth procedure is customization service
c
i
j
k
1
The unit operating cost of theith FLSP for the jth customer when the kth procedure is mass service
c
i
j
k
2
The unit operating cost of theith FLSP for the jth customer when the kth procedure is customization service
ε
i
1
The weight of the single satisfaction of theith FLSP
ε
i
2
The weight of the overall satisfaction of theith FLSP
λ
i
j
The preference of theith FLSP for the jth customer
k
Thekth service procedure
K
j
The procedure to start the customization service of thejth customer
K
0
The minimum of allKj
K
′
The CODP for all customers
θ
Order difference
ω
Order difference tolerance coefficient of LSI
δ
Relationship cost coefficient of LSI
f
j
k
The order quantity of thejth customer in the kth procedure
Z
1
The total satisfaction of FLSPs in LSSC
Z
2
The total cost of LSI in LSSC
Z
3
The total customized degree of LSSC
Z
2
*
The optimal solution of objective functionZ2 when not considering the objective function Z1 and Z3
F
1
*
The optimal solution of objective functionZ1 when not considering the objective function Z2 and Z3
F
3
′
The value ofZ3, when the optimal solution is F1*
F
3
*
The optimal solution of objective functionZ3 when not considering the objective function Z1 and Z2
F
1
′
The value ofZ1, when the optimal solution is F3*
α
1
Weight coefficients of the objective functionZ1 in weighted objective function
α
2
Weight coefficients of the objective functionZ3 in weighted objective function
Z
The objective function synthesized byZ1 and Z3 is also called the comprehensive performance of LSSC
Z
*
The optimal solution ofZ


### 3.2. Model Building

Previous order allocation models are mostly based on the costs and orders service level [13, 16, 24]. While conducting the order allocation of the LSSC under the MCLS environment, the LSI is required to consider multiobjectives, such as the total cost of the LSI, the total satisfaction of FLSPs, and the satisfaction of customized demand. To maximize the satisfaction of FLSPs, minimize the total cost of the LSI, and maximize the customization degree for customer demand, this paper establishes an order allocation model of LSSC that is constrained by meeting the demand, CODP, and order difference tolerance coefficient. The greatest difference between the previous model and ours is that we considered the impact of changing CODP on order allocation result, so that it can reflects the order allocation characteristics on MCLS better. The specific modeling procedure is described in the following sections. Sections 3.2.1 to 3.2.4 discuss the building processes of several aspects, namely, the objective function of FLSP satisfaction, LSI total cost, the objective function of customized degree for customer demand, and the three constraints, respectively.


#### 3.2.1. Maximizing the Satisfaction of FLSPs

Maximizing the satisfaction of FLSPs is considered the first process by the LSI in order allocation. Based on the satisfaction function of reference [1], the satisfaction of the ith FLSP for the jth customer in the kth procedure is as follows:
(1)ρijk={bijkxijk,xijk≻bijk,ρijk0+(1-ρijk0)*xijk-aijkbijk-aijk,aijk≤xijk≤bijk,xijkaijk*ρijk0,0≤x≤aijk,
where ρijk0 is the initial satisfaction of the ith FLSP in the kth procedure when the order value is equal to the minimum service capacity (1). Each FLSP has different service capacity intervals while conducting mass service and customization service, which can be expressed as k≤K′, [aijk,bijk]=[aijk1,bijk1]; k≻K′, [aijk,bijk]=[aijk2,bijk2]. In the same manner, each FLSP has different initial satisfaction on mass service and customization service. This initial satisfaction can be expressed as k≤K′, [aijk,bijk]=[aijk1,bijk1]; k≻K′, [aijk,bijk]=[aijk2,bijk2]. When 0≤x≤aijk, the satisfaction of the FLSP increases with increasing order value. Similarly, when aijk≤xijk≤bijk, the satisfaction of the FLSP increases with increasing order value (the satisfaction maximum is one). However, when xijk≥bijk, the satisfaction of the FLSP decreases with increasing order value because the order exceeds the upper limit of the service capacity.Given that each FLSP has different preferences for each customer, satisfaction comprises single satisfaction and overall satisfaction. The weight of the single satisfaction of theith FLSP is εi1, and the weight of the overall satisfaction of the ith FLSP is εi2 (εi1+εi2=1).(1)
Single satisfaction: the satisfaction of the FLSP for the order allocation result of one customer. Assuming that the preference of theith FLSP for the jth customer is λij, ∑j=1Mλij is equal to one, thus, the single satisfaction of the ith FLSP is εi1(∑j=1Mλij(1/Lj)∑k=1Ljρijk).(2)
Overall satisfaction: the satisfaction of the FLSP for the order allocation result of all customers. If the FLSP can service two customers, the FLSP expects that the satisfaction for both customers is high. Usually, we use the average satisfaction of multiple customers to represent the overall satisfaction. Thus, the overall satisfaction of theith FLSP is εi2((1/M)∑j=1M(1/Lj)∑k=1Ljρijk).The total satisfaction ofith FLSP is εi1(∑j=1Mλij(1/Lj)∑k=1Ljρijk)+εi2((1/M)∑j=1M(1/Lj)∑k=1Ljρijk).Therefore, the objective function of maximizing the satisfaction of FLSPs is expressed as follows:(2)maxZ1=1n∑i=1n(+εi2(1M∑j=1M1Lj∑k=1Ljρijk))εi1(∑j=1Mλij1Lj∑k=1Ljρijk)+εi2(1M∑j=1M1Lj∑k=1Ljρijk)).


#### 3.2.2. Minimizing the Total Cost of the LSI

In the order allocation process, the LSI also considers minimizing the cost. The total cost of the LSI is related not only to the number of orders assigned but also to the position of the CODP. When CODP approaches the customer, the quantity of procedures and order of mass services increase; the order of mass service is also obvious. Therefore, the scale effect is more noticeable.The total cost of the LSI consists of two parts, namely, mass service cost and customization service cost. The unit operating cost of theith FLSP for the jth customer is cijk1 when the kth procedure is mass service (3). The unit operating cost of the ith FLSP for the jth customer is cijk2 when the kth procedure is customization service (4). Consider
(3)cijk1=Aik-Bikxijk,(4)cijk2=Ci,wherecijk1 displays a scale effect, which decreases with the increasing of xijk (3). A larger xijk corresponds to a smaller cijk1. However, because of the characteristics of customization, cijk2 does not change with the variation of xijk.When the procedure of mass service isK′ because of the existence of the scale effect, an increase in K′ corresponds to a smaller cijk1. Hence, we introduce the scale effect coefficient τ, which is related to the number of procedures. The cost of the ith FLSP for the jth customer in all procedures is expressed as follows:
(5)∑k=1K′cijk1*xijk-τK′∑k=1K′cijk1*xijk=(1-τK′)∑k=1K′cijk1*xijk.The cost of all FLSPs is summarized. Thereafter, the objective function of minimizing the total cost of the LSI is obtained:(6)minZ2=∑i=1n∑j=1M((1-τK′)∑k=1K′cijk1*xijk+∑k=K′+1Ljcijk2*xijk).


#### 3.2.3. Maximizing the Customized Degree for Customer Demand

Under the MCLS environment, each customer will intend to maximize their own customized degree. Therefore, the customized degree for customer demand is one of the most important objectives. However, given that the customized degree for customer demand will affect the total cost of the LSI, the former cannot be too high. Mass customization involves reducing total cost as much as possible and maximizing the customized demand (i.e., increase the cost of customization services). Thus, we multiply the ratio of the customization service cost and total cost by the weight to represent the customized degree of each customer. We then summarize the customized degree of each customer to obtain the customized degree of all customers:(7)maxZ3=∑j=1Mβj(+∑k=K′+1Ljcijk2*xijk((1-τK′)∑k=1K′cijk1*xijk))-1)(∑i=1n∑k=K′+1Ljcijk2*xijk)×(∑i=1n((1-τK′)∑k=1K′cijk1*xijk+∑k=K′+1Ljcijk2*xijk((1-τK′)∑k=1K′cijk1*xijk))-1).


#### 3.2.4. Order Allocation Constraints

In the order allocation process, the LSI also needs to meet several constraints. First, the total order quantity of all FLSPs should be equal to the customer demand in thekth procedure. Second, the K′ (CODP) should be between the second procedure and the K0 procedure (K0=min{K1,K2,…,KM}). Third, given that the service procedures of each customer and the procedures to start the customization service are different, the LSI needs to consider the order difference θ of customers while serving multiple customers. To conduct mass customization service, the order difference cannot be extremely high. That is, θ cannot exceed the upper limit ω. Finally, xijk must be nonnegative. The four constraints can be expressed as (8), (9), (10), and (11):
(8)∑i=1n∑j=1Mxijk=∑j=1Mfjk,(9)2≤K′≤K0,(10)θ=1M∑j=1MKj-K′Kj≤ω,(11)xijk≥0,i=1,2,…,n;j=1,2,…,M;k=1,2,…,Lj.As shown in (10), θ can be expressed as (12):
(12)θ=1M∑j=1MKj-K′Kj.Equation (12) shows that θ is related to Kj and K′. If Kj is equal to K′, θ is zero. A larger gap between Kj and K′ corresponds to a higher θ value.Combining the three objective functions given in Sections3.2.1 to 3.2.3, we establish a multiobjective order allocation model under the MCLS environment ((13) to (19)). Consider
(13)maxZ1=1n∑i=1n(εi1(∑j=1Mλij1Lj∑k=1Ljρijk)+εi2(1M∑j=1M1Lj∑k=1Ljρijk)),(14)minZ2=∑i=1n∑j=1M((1-τK′)∑k=1K′cijk1*xijk+∑k=K′+1Ljcijk2*xijk),(15)maxZ3=∑j=1Mβj(+∑k=K′+1Ljcijk2*xijk))-1)(∑i=1n∑k=K′+1Ljcijk2*xijk)×(∑i=1n((1-τK′)∑k=1K′cijk1*xijk+∑k=K′+1Ljcijk2*xijk))-1),(16)s.t∑i=1n∑j=1Mxijk=∑j=1Mfjk,(17)2≤K′≤K0,(18)θ=1M∑j=1MKj-K′Kj≤ω,(19)xijk≥0,i=1,2,…,n;j=1,2,…,M;k=1,2,…,Lj.In the objective functions, (13) maximizes the satisfaction of FLSPs; (14) minimizes the total cost of the LSI; (15) maximizes the customized degree for customer demand.For the constraints, (16) implies that the service capacity supply is equal to the order demand in each procedure; (17) means that K′ (CODP) should be between the second procedure and K0 procedure. Equation (18) indicates that the order difference should be less than the order difference tolerance coefficient. Equation (19) denotes that the order assigned to each FLSP must be nonnegative.


## 3.1. Model Description

We considered a two-echelon LSSC, which consists of one LSI and many FLSPs. After receiving orders from multiple customers, the LSI analyzes the customized and universal service requirements of these customer orders and allocates them to more than one FLSP. Every logistic service for customer orders comprises multiple service processes. Each process needs one kind of service that can be provided by many cooperative FLSPs. The service capability of each FLSP is possibly different. We assumed that[aijk,bijk] denotes the service ability interval provided by FLSP ith for customer jth in the procedure kth; xijk stands for the jth customer’s order assigned to the ith FLSP in the kth procedure.Table1 shows the notations of the model. The other assumptions of the model are as follows.(1)
Each FLSPi  (i=1,2,…,n) can provide mass service and customization service for the LSI. The total service procedures of customer j are Lj  (j=1,2,…,M), and the order demand is Dj.(2)
Kj represents the procedure to start the customization service of the jth customer; this service is not necessarily the same for all customers. The LSI needs to consider the Kj of all customers and set a common CODP, which is signified by K′ for all customers. K′ is mass service, and K′ is customization service. K′ should be set before K0  (K0=min{K1,K2,…,KM}) because the CODP must be set before all Kj. Given the characteristics of the MCLS, Procedure 1 must be mass service; hence, 2≤K′≤K0).(3)
Given the importance of different customizations for each customer, the weight of thejth customer is βj (∑j=1Mβj=1).(4)
Customer orders are assigned to FLSPs through the LSI. FLSPs will accomplish the order assigned by the LSI. No game exists between LSI and FLSPs.Table 1
Model notations.
Notations
Description
[
a
i
j
k
,
b
i
j
k
]
The service ability interval of theith FLSP for jth customer in the kth procedure
[
a
i
j
k
1
,
b
i
j
k
1
]
The value of[aijk,bijk]  when the kth procedure is mass service
[
a
i
j
k
2
,
b
i
j
k
2
]
The value of[aijk,bijk] when the kth procedure is customization service
x
i
j
k
Thejth customer’s order assigned to the ith FLSP in the kth procedure
τ
Scale effect coefficient
D
j
The demand of thejth customer
L
j
The total service procedures of thejth customer
β
j
The weight of thejth customer
ρ
i
j
k
The satisfaction of theith FLSP for the order of the jth customer in the kth procedure
ρ
i
j
k
0
The initial satisfaction of theith FLSP in the kth procedure when the order is equal to the minimum service capacity of  ith FLSP
ρ
i
j
k
1
0
The value ofρijk0 when the kth procedure is mass service
ρ
i
j
k
2
0
The value ofρijk0 when the kth procedure is customization service
c
i
j
k
1
The unit operating cost of theith FLSP for the jth customer when the kth procedure is mass service
c
i
j
k
2
The unit operating cost of theith FLSP for the jth customer when the kth procedure is customization service
ε
i
1
The weight of the single satisfaction of theith FLSP
ε
i
2
The weight of the overall satisfaction of theith FLSP
λ
i
j
The preference of theith FLSP for the jth customer
k
Thekth service procedure
K
j
The procedure to start the customization service of thejth customer
K
0
The minimum of allKj
K
′
The CODP for all customers
θ
Order difference
ω
Order difference tolerance coefficient of LSI
δ
Relationship cost coefficient of LSI
f
j
k
The order quantity of thejth customer in the kth procedure
Z
1
The total satisfaction of FLSPs in LSSC
Z
2
The total cost of LSI in LSSC
Z
3
The total customized degree of LSSC
Z
2
*
The optimal solution of objective functionZ2 when not considering the objective function Z1 and Z3
F
1
*
The optimal solution of objective functionZ1 when not considering the objective function Z2 and Z3
F
3
′
The value ofZ3, when the optimal solution is F1*
F
3
*
The optimal solution of objective functionZ3 when not considering the objective function Z1 and Z2
F
1
′
The value ofZ1, when the optimal solution is F3*
α
1
Weight coefficients of the objective functionZ1 in weighted objective function
α
2
Weight coefficients of the objective functionZ3 in weighted objective function
Z
The objective function synthesized byZ1 and Z3 is also called the comprehensive performance of LSSC
Z
*
The optimal solution ofZ


## 3.2. Model Building

Previous order allocation models are mostly based on the costs and orders service level [13, 16, 24]. While conducting the order allocation of the LSSC under the MCLS environment, the LSI is required to consider multiobjectives, such as the total cost of the LSI, the total satisfaction of FLSPs, and the satisfaction of customized demand. To maximize the satisfaction of FLSPs, minimize the total cost of the LSI, and maximize the customization degree for customer demand, this paper establishes an order allocation model of LSSC that is constrained by meeting the demand, CODP, and order difference tolerance coefficient. The greatest difference between the previous model and ours is that we considered the impact of changing CODP on order allocation result, so that it can reflects the order allocation characteristics on MCLS better. The specific modeling procedure is described in the following sections. Sections 3.2.1 to 3.2.4 discuss the building processes of several aspects, namely, the objective function of FLSP satisfaction, LSI total cost, the objective function of customized degree for customer demand, and the three constraints, respectively.


### 3.2.1. Maximizing the Satisfaction of FLSPs

Maximizing the satisfaction of FLSPs is considered the first process by the LSI in order allocation. Based on the satisfaction function of reference [1], the satisfaction of the ith FLSP for the jth customer in the kth procedure is as follows:
(1)ρijk={bijkxijk,xijk≻bijk,ρijk0+(1-ρijk0)*xijk-aijkbijk-aijk,aijk≤xijk≤bijk,xijkaijk*ρijk0,0≤x≤aijk,
where ρijk0 is the initial satisfaction of the ith FLSP in the kth procedure when the order value is equal to the minimum service capacity (1). Each FLSP has different service capacity intervals while conducting mass service and customization service, which can be expressed as k≤K′, [aijk,bijk]=[aijk1,bijk1]; k≻K′, [aijk,bijk]=[aijk2,bijk2]. In the same manner, each FLSP has different initial satisfaction on mass service and customization service. This initial satisfaction can be expressed as k≤K′, [aijk,bijk]=[aijk1,bijk1]; k≻K′, [aijk,bijk]=[aijk2,bijk2]. When 0≤x≤aijk, the satisfaction of the FLSP increases with increasing order value. Similarly, when aijk≤xijk≤bijk, the satisfaction of the FLSP increases with increasing order value (the satisfaction maximum is one). However, when xijk≥bijk, the satisfaction of the FLSP decreases with increasing order value because the order exceeds the upper limit of the service capacity.Given that each FLSP has different preferences for each customer, satisfaction comprises single satisfaction and overall satisfaction. The weight of the single satisfaction of theith FLSP is εi1, and the weight of the overall satisfaction of the ith FLSP is εi2 (εi1+εi2=1).(1)
Single satisfaction: the satisfaction of the FLSP for the order allocation result of one customer. Assuming that the preference of theith FLSP for the jth customer is λij, ∑j=1Mλij is equal to one, thus, the single satisfaction of the ith FLSP is εi1(∑j=1Mλij(1/Lj)∑k=1Ljρijk).(2)
Overall satisfaction: the satisfaction of the FLSP for the order allocation result of all customers. If the FLSP can service two customers, the FLSP expects that the satisfaction for both customers is high. Usually, we use the average satisfaction of multiple customers to represent the overall satisfaction. Thus, the overall satisfaction of theith FLSP is εi2((1/M)∑j=1M(1/Lj)∑k=1Ljρijk).The total satisfaction ofith FLSP is εi1(∑j=1Mλij(1/Lj)∑k=1Ljρijk)+εi2((1/M)∑j=1M(1/Lj)∑k=1Ljρijk).Therefore, the objective function of maximizing the satisfaction of FLSPs is expressed as follows:(2)maxZ1=1n∑i=1n(+εi2(1M∑j=1M1Lj∑k=1Ljρijk))εi1(∑j=1Mλij1Lj∑k=1Ljρijk)+εi2(1M∑j=1M1Lj∑k=1Ljρijk)).


### 3.2.2. Minimizing the Total Cost of the LSI

In the order allocation process, the LSI also considers minimizing the cost. The total cost of the LSI is related not only to the number of orders assigned but also to the position of the CODP. When CODP approaches the customer, the quantity of procedures and order of mass services increase; the order of mass service is also obvious. Therefore, the scale effect is more noticeable.The total cost of the LSI consists of two parts, namely, mass service cost and customization service cost. The unit operating cost of theith FLSP for the jth customer is cijk1 when the kth procedure is mass service (3). The unit operating cost of the ith FLSP for the jth customer is cijk2 when the kth procedure is customization service (4). Consider
(3)cijk1=Aik-Bikxijk,(4)cijk2=Ci,wherecijk1 displays a scale effect, which decreases with the increasing of xijk (3). A larger xijk corresponds to a smaller cijk1. However, because of the characteristics of customization, cijk2 does not change with the variation of xijk.When the procedure of mass service isK′ because of the existence of the scale effect, an increase in K′ corresponds to a smaller cijk1. Hence, we introduce the scale effect coefficient τ, which is related to the number of procedures. The cost of the ith FLSP for the jth customer in all procedures is expressed as follows:
(5)∑k=1K′cijk1*xijk-τK′∑k=1K′cijk1*xijk=(1-τK′)∑k=1K′cijk1*xijk.The cost of all FLSPs is summarized. Thereafter, the objective function of minimizing the total cost of the LSI is obtained:(6)minZ2=∑i=1n∑j=1M((1-τK′)∑k=1K′cijk1*xijk+∑k=K′+1Ljcijk2*xijk).


### 3.2.3. Maximizing the Customized Degree for Customer Demand

Under the MCLS environment, each customer will intend to maximize their own customized degree. Therefore, the customized degree for customer demand is one of the most important objectives. However, given that the customized degree for customer demand will affect the total cost of the LSI, the former cannot be too high. Mass customization involves reducing total cost as much as possible and maximizing the customized demand (i.e., increase the cost of customization services). Thus, we multiply the ratio of the customization service cost and total cost by the weight to represent the customized degree of each customer. We then summarize the customized degree of each customer to obtain the customized degree of all customers:(7)maxZ3=∑j=1Mβj(+∑k=K′+1Ljcijk2*xijk((1-τK′)∑k=1K′cijk1*xijk))-1)(∑i=1n∑k=K′+1Ljcijk2*xijk)×(∑i=1n((1-τK′)∑k=1K′cijk1*xijk+∑k=K′+1Ljcijk2*xijk((1-τK′)∑k=1K′cijk1*xijk))-1).


### 3.2.4. Order Allocation Constraints

In the order allocation process, the LSI also needs to meet several constraints. First, the total order quantity of all FLSPs should be equal to the customer demand in thekth procedure. Second, the K′ (CODP) should be between the second procedure and the K0 procedure (K0=min{K1,K2,…,KM}). Third, given that the service procedures of each customer and the procedures to start the customization service are different, the LSI needs to consider the order difference θ of customers while serving multiple customers. To conduct mass customization service, the order difference cannot be extremely high. That is, θ cannot exceed the upper limit ω. Finally, xijk must be nonnegative. The four constraints can be expressed as (8), (9), (10), and (11):
(8)∑i=1n∑j=1Mxijk=∑j=1Mfjk,(9)2≤K′≤K0,(10)θ=1M∑j=1MKj-K′Kj≤ω,(11)xijk≥0,i=1,2,…,n;j=1,2,…,M;k=1,2,…,Lj.As shown in (10), θ can be expressed as (12):
(12)θ=1M∑j=1MKj-K′Kj.Equation (12) shows that θ is related to Kj and K′. If Kj is equal to K′, θ is zero. A larger gap between Kj and K′ corresponds to a higher θ value.Combining the three objective functions given in Sections3.2.1 to 3.2.3, we establish a multiobjective order allocation model under the MCLS environment ((13) to (19)). Consider
(13)maxZ1=1n∑i=1n(εi1(∑j=1Mλij1Lj∑k=1Ljρijk)+εi2(1M∑j=1M1Lj∑k=1Ljρijk)),(14)minZ2=∑i=1n∑j=1M((1-τK′)∑k=1K′cijk1*xijk+∑k=K′+1Ljcijk2*xijk),(15)maxZ3=∑j=1Mβj(+∑k=K′+1Ljcijk2*xijk))-1)(∑i=1n∑k=K′+1Ljcijk2*xijk)×(∑i=1n((1-τK′)∑k=1K′cijk1*xijk+∑k=K′+1Ljcijk2*xijk))-1),(16)s.t∑i=1n∑j=1Mxijk=∑j=1Mfjk,(17)2≤K′≤K0,(18)θ=1M∑j=1MKj-K′Kj≤ω,(19)xijk≥0,i=1,2,…,n;j=1,2,…,M;k=1,2,…,Lj.In the objective functions, (13) maximizes the satisfaction of FLSPs; (14) minimizes the total cost of the LSI; (15) maximizes the customized degree for customer demand.For the constraints, (16) implies that the service capacity supply is equal to the order demand in each procedure; (17) means that K′ (CODP) should be between the second procedure and K0 procedure. Equation (18) indicates that the order difference should be less than the order difference tolerance coefficient. Equation (19) denotes that the order assigned to each FLSP must be nonnegative.


## 3.2.1. Maximizing the Satisfaction of FLSPs

Maximizing the satisfaction of FLSPs is considered the first process by the LSI in order allocation. Based on the satisfaction function of reference [1], the satisfaction of the ith FLSP for the jth customer in the kth procedure is as follows:
(1)ρijk={bijkxijk,xijk≻bijk,ρijk0+(1-ρijk0)*xijk-aijkbijk-aijk,aijk≤xijk≤bijk,xijkaijk*ρijk0,0≤x≤aijk,
where ρijk0 is the initial satisfaction of the ith FLSP in the kth procedure when the order value is equal to the minimum service capacity (1). Each FLSP has different service capacity intervals while conducting mass service and customization service, which can be expressed as k≤K′, [aijk,bijk]=[aijk1,bijk1]; k≻K′, [aijk,bijk]=[aijk2,bijk2]. In the same manner, each FLSP has different initial satisfaction on mass service and customization service. This initial satisfaction can be expressed as k≤K′, [aijk,bijk]=[aijk1,bijk1]; k≻K′, [aijk,bijk]=[aijk2,bijk2]. When 0≤x≤aijk, the satisfaction of the FLSP increases with increasing order value. Similarly, when aijk≤xijk≤bijk, the satisfaction of the FLSP increases with increasing order value (the satisfaction maximum is one). However, when xijk≥bijk, the satisfaction of the FLSP decreases with increasing order value because the order exceeds the upper limit of the service capacity.Given that each FLSP has different preferences for each customer, satisfaction comprises single satisfaction and overall satisfaction. The weight of the single satisfaction of theith FLSP is εi1, and the weight of the overall satisfaction of the ith FLSP is εi2 (εi1+εi2=1).(1)
Single satisfaction: the satisfaction of the FLSP for the order allocation result of one customer. Assuming that the preference of theith FLSP for the jth customer is λij, ∑j=1Mλij is equal to one, thus, the single satisfaction of the ith FLSP is εi1(∑j=1Mλij(1/Lj)∑k=1Ljρijk).(2)
Overall satisfaction: the satisfaction of the FLSP for the order allocation result of all customers. If the FLSP can service two customers, the FLSP expects that the satisfaction for both customers is high. Usually, we use the average satisfaction of multiple customers to represent the overall satisfaction. Thus, the overall satisfaction of theith FLSP is εi2((1/M)∑j=1M(1/Lj)∑k=1Ljρijk).The total satisfaction ofith FLSP is εi1(∑j=1Mλij(1/Lj)∑k=1Ljρijk)+εi2((1/M)∑j=1M(1/Lj)∑k=1Ljρijk).Therefore, the objective function of maximizing the satisfaction of FLSPs is expressed as follows:(2)maxZ1=1n∑i=1n(+εi2(1M∑j=1M1Lj∑k=1Ljρijk))εi1(∑j=1Mλij1Lj∑k=1Ljρijk)+εi2(1M∑j=1M1Lj∑k=1Ljρijk)).


## 3.2.2. Minimizing the Total Cost of the LSI

In the order allocation process, the LSI also considers minimizing the cost. The total cost of the LSI is related not only to the number of orders assigned but also to the position of the CODP. When CODP approaches the customer, the quantity of procedures and order of mass services increase; the order of mass service is also obvious. Therefore, the scale effect is more noticeable.The total cost of the LSI consists of two parts, namely, mass service cost and customization service cost. The unit operating cost of theith FLSP for the jth customer is cijk1 when the kth procedure is mass service (3). The unit operating cost of the ith FLSP for the jth customer is cijk2 when the kth procedure is customization service (4). Consider
(3)cijk1=Aik-Bikxijk,(4)cijk2=Ci,wherecijk1 displays a scale effect, which decreases with the increasing of xijk (3). A larger xijk corresponds to a smaller cijk1. However, because of the characteristics of customization, cijk2 does not change with the variation of xijk.When the procedure of mass service isK′ because of the existence of the scale effect, an increase in K′ corresponds to a smaller cijk1. Hence, we introduce the scale effect coefficient τ, which is related to the number of procedures. The cost of the ith FLSP for the jth customer in all procedures is expressed as follows:
(5)∑k=1K′cijk1*xijk-τK′∑k=1K′cijk1*xijk=(1-τK′)∑k=1K′cijk1*xijk.The cost of all FLSPs is summarized. Thereafter, the objective function of minimizing the total cost of the LSI is obtained:(6)minZ2=∑i=1n∑j=1M((1-τK′)∑k=1K′cijk1*xijk+∑k=K′+1Ljcijk2*xijk).


## 3.2.3. Maximizing the Customized Degree for Customer Demand

Under the MCLS environment, each customer will intend to maximize their own customized degree. Therefore, the customized degree for customer demand is one of the most important objectives. However, given that the customized degree for customer demand will affect the total cost of the LSI, the former cannot be too high. Mass customization involves reducing total cost as much as possible and maximizing the customized demand (i.e., increase the cost of customization services). Thus, we multiply the ratio of the customization service cost and total cost by the weight to represent the customized degree of each customer. We then summarize the customized degree of each customer to obtain the customized degree of all customers:(7)maxZ3=∑j=1Mβj(+∑k=K′+1Ljcijk2*xijk((1-τK′)∑k=1K′cijk1*xijk))-1)(∑i=1n∑k=K′+1Ljcijk2*xijk)×(∑i=1n((1-τK′)∑k=1K′cijk1*xijk+∑k=K′+1Ljcijk2*xijk((1-τK′)∑k=1K′cijk1*xijk))-1).


## 3.2.4. Order Allocation Constraints

In the order allocation process, the LSI also needs to meet several constraints. First, the total order quantity of all FLSPs should be equal to the customer demand in thekth procedure. Second, the K′ (CODP) should be between the second procedure and the K0 procedure (K0=min{K1,K2,…,KM}). Third, given that the service procedures of each customer and the procedures to start the customization service are different, the LSI needs to consider the order difference θ of customers while serving multiple customers. To conduct mass customization service, the order difference cannot be extremely high. That is, θ cannot exceed the upper limit ω. Finally, xijk must be nonnegative. The four constraints can be expressed as (8), (9), (10), and (11):
(8)∑i=1n∑j=1Mxijk=∑j=1Mfjk,(9)2≤K′≤K0,(10)θ=1M∑j=1MKj-K′Kj≤ω,(11)xijk≥0,i=1,2,…,n;j=1,2,…,M;k=1,2,…,Lj.As shown in (10), θ can be expressed as (12):
(12)θ=1M∑j=1MKj-K′Kj.Equation (12) shows that θ is related to Kj and K′. If Kj is equal to K′, θ is zero. A larger gap between Kj and K′ corresponds to a higher θ value.Combining the three objective functions given in Sections3.2.1 to 3.2.3, we establish a multiobjective order allocation model under the MCLS environment ((13) to (19)). Consider
(13)maxZ1=1n∑i=1n(εi1(∑j=1Mλij1Lj∑k=1Ljρijk)+εi2(1M∑j=1M1Lj∑k=1Ljρijk)),(14)minZ2=∑i=1n∑j=1M((1-τK′)∑k=1K′cijk1*xijk+∑k=K′+1Ljcijk2*xijk),(15)maxZ3=∑j=1Mβj(+∑k=K′+1Ljcijk2*xijk))-1)(∑i=1n∑k=K′+1Ljcijk2*xijk)×(∑i=1n((1-τK′)∑k=1K′cijk1*xijk+∑k=K′+1Ljcijk2*xijk))-1),(16)s.t∑i=1n∑j=1Mxijk=∑j=1Mfjk,(17)2≤K′≤K0,(18)θ=1M∑j=1MKj-K′Kj≤ω,(19)xijk≥0,i=1,2,…,n;j=1,2,…,M;k=1,2,…,Lj.In the objective functions, (13) maximizes the satisfaction of FLSPs; (14) minimizes the total cost of the LSI; (15) maximizes the customized degree for customer demand.For the constraints, (16) implies that the service capacity supply is equal to the order demand in each procedure; (17) means that K′ (CODP) should be between the second procedure and K0 procedure. Equation (18) indicates that the order difference should be less than the order difference tolerance coefficient. Equation (19) denotes that the order assigned to each FLSP must be nonnegative.


## 4. Model Solution

The model is a multiobjective programming problem with three objectives and three constraints. Multiobjective programming problems have numerous mature solutions, such as the evaluation function method including linear weighting method, reference target method, maximin method [9], goal programming method [25], delaminating sequence method [26], and subordinate function method [27]. For the specific issues, it is necessary to choose an appropriate solution method to solve the practical problems. Given that conflicts and incommensurability between goals exists, finding an absolute optimal solution is difficult. Therefore, a compromise between the various objectives is commonly performed to obtain the Pareto-optimal solution for decision makers. Specifically, two approaches can be used to solve multiobjective programming problem. One is to balance between objectives, assign weight to each of the objective functions and switch the multiobjective functions to a single objective function. The other method is to reduce the number of objective functions and transform the reduced objective functions into constraints. However, on MCLS, the LSI should not only consider the total cost target but also the satisfaction target and the customized degree target. Generally, supply chain management emphasizes the strategic partnership. In order to maintain long-term cooperation with FLSPs and customers, in order allocation, the LSI should keep the total cost within a specific range instead of merely seeking for the minimum. Maximizing the satisfaction and the customized degree should be the target. As a result, we transform the total cost objective into a constraint to reduce the number of objective functions and transform the target of satisfaction customized degree into a single objective, after which the model is solved.Based on the previous solving methods, in order allocation, considering these actual situations, we introduce a parameter called the relationship cost coefficientδ into our model. This parameter is used to represent the increasing cost limit of the LSI to enhance the cooperation and satisfaction of FLSPs and customers. Hence, Z2 can become a constraint. Thereafter, we choose a linear weighting method, which is the most typical method, to transform the multiobjective Z1, Z3 into a single objective Z and solve the model. The specific steps in solving the model are as follows.Step 1.
The consideration of the objective functionsZ1, Z3 is set aside to consider the case of objective function Z2. Thereafter, the relationship cost coefficient δ is introduced. The objective can then be translated into a constraint:
(20)Z2≤(1+δ)Z2*.
This approach denotes that the LSI can minimize and maintain the total cost within a certain range.δ=0.2 means that the LSI can exist with a total cost less than 1.2 times of Z2* to increase the satisfaction and customized degree. The LSI considers this input to be conducive in the completion of orders.Step 2.
The model hasZ1 and Z3. The linear weighing method is used to transform the objective Z1 and Z3 into a new objective function maxZ=α1Z1+α2Z3 (α1+α2=1). Determining the weight (α1,α2) is the main problem of this approach. α1, α2 can be determined by the following equation:
(21)α1=F3*-F3′F1*-F1′+F3*-F3′,α2=F1*-F1′F1*-F1′+F3*-F3′.Step 3.
After determiningα1 and α2, the model is transformed into a single-objective planning model:
(22)maxZ=α1Z1+α2Z3,s.tZ2≤(1+δ)Z2*,∑i=1n∑j=1Mxijk=∑j=1Mfjk,2≤K′≤K0,θ=1M∑j=1MKj-K′Kj≤ω,xijk≥0,i=1,2,…,n;j=1,2,…,M;k=1,2,…,Lj.
In solving the objective functionZ, the xijk obtained is the optimal solution Z*. Z reflects the comprehensive performance of the LSSC. A bigger Z corresponds to better comprehensive performance.


## 5. Numerical Analysis

This section illustrates the validity of the model by numerical analysis, explores the influence of related factors on time scheduling results, and provides several management insights. First, we use the Lingo 12 software to implement the numerical analysis. The results are shown in Section5.1. From the perspective of LSI order management decision, we choose and study the influence of following factors (τ,ω,δ) on order allocation decision (Section 5.2 to Section 5.5). τ is an important factor to reflect the scale effect on MCLS. ω is also an important factor to reflect the influence of order difference on LSI’s decision making. δ indicates the cost volatility with which the LSI can live. We will study the influence of the three (τ,ω,δ) on order allocation separately.


### 5.1. Model Solution

Figures1 and 2 show that each customer (A, B, and C) has an order and needs the LSI Q to provide transportation services. Q assigns the three orders to the FLSPs (a, b, c, d, and e). The demand of each customer is 60, 100, and 80. The number of service procedures of each customer is 8, 7, and 7. K1=6, K2=5, K3=5, K0=min(K1,K2,K3)=5, τ=0.05, ω=0.4, and δ=20%. The values of other parameters are shown in Tables 2 and 3.Table 2
Parameter data (1).
[
a
i
j
k
,
b
i
j
k
]
[
a
i
j
k
1
,
b
i
j
k
1
]
[
a
i
j
k
2
,
b
i
j
k
2
]
i
=
1
[
40
,
80
]
[
20
,
40
]
i
=
2
[
30
,
70
]
[
20
,
35
]
i
=
3
[
35
,
90
]
[
25
,
40
]
i
=
4
[
30
,
60
]
[
25
,
50
]
i
=
5
[
30
,
80
]
[
10
,
40
]
c
i
j
k
c
i
j
k
1
c
i
j
k
2
i
=
1
10
-
0.02
x
i
j
k
15
i
=
2
11
-
0.02
x
i
j
k
18
i
=
3
9
-
0.02
x
i
j
k
16
i
=
4
12
-
0.02
x
i
j
k
20
i
=
5
14
-
0.02
x
i
j
k
22
ρ
i
j
k
ρ
i
j
k
1
ρ
i
j
k
2
i
=
1
0.2
0.3
i
=
2
0.35
0.35
i
=
3
0.15
0.3
i
=
4
0.3
0.2
i
=
5
0.25
0.3
ε
i
ε
i
1
ε
i
2
i
=
1
0.6
0.4
i
=
2
0.5
0.5
i
=
3
0.7
0.3
i
=
4
0.4
0.6
i
=
5
0.35
0.65Table 3
Parameter data (3).
λ
i
j
j
=
1
j
=
2
j
=
3
i
=
1
0.3
0.3
0.4
i
=
2
0.2
0.4
0.4
i
=
3
0.35
0.35
0.3
i
=
4
0.25
0.35
0.4
i
=
5
0.3
0.4
0.3
β
j
0.3
0.4
0.3
D
j
60
100
80
K
j
6
5
5
L
j
8
7
7Figure 1
Structure of the LSSC.Figure 2
Logistics service procedure of customer A, B, and C. Note: The red circle denotes the procedure that starts the customization service of each customer.By using the solution method in Section4, Tables 4 and 5 show the result and performance of the order allocation, respectively.Table 4
Order allocation results.
Procedure
1–5
6–7
8
CustomerOrder valueFLSP
A
B
C
A
B
C
A
a
0
4.31
0
31.59
18.20
37.79
31.59
b
0
6.33
0
6.14
37.86
0.53
6.14
c
60
4.99
80
4.63
40
37.82
4.63
d
0
4.93
0
6.41
0.35
0
6.41
e
0
79.44
0
11.23
3.60
3.86
11.23Table 5
Order allocation performance.
Total cost of LSSC
17472.020
CODP
5
Customized degree
0.522
Z
0.391


### 5.2. Influence of Scale Effect Coefficientτ on Order Allocation

To analyze the influence ofτ on the related factors of order allocation (Z2,CODP,Z3,Z), we make τ equal to 0, 0.025, 0.05, 0.075, 0.1, 0.125, and 0.15. The data is shown in Table 6. Figures 3 to 6 are plotted with the data from Table 6.Table 6
Influence ofτ on order allocation.
τ
Z
2
Z
3
CODP
Z
0
20280
0.596
4
0.403
0.025
18960.000
0.524
5
0.388
0.05
17472.020
0.522
5
0.391
0.075
16320.000
0.580
5
0.417
0.1
15000.000
0.638
5
0.469
0.125
13680.000
0.646
5
0.472
0.15
12360.000
0.783
5
0.582Figure 3
Influence ofτ on Z2 (total cost of LSI).Figure3 shows that τ has significant influence on the total cost of LSI Z2. A greater τ corresponds to a smaller Z2. The scale effect of the cost is more obvious with increasing τ; thus, the total cost decreases.Figure4 shows that τ has a certain influence on CODP. If τ=0, the scale effect is nonexistent and the CODP is in the fourth procedure. With τ≥0.025, the CODP changes to the fifth procedure because the scale effect of cost is more obvious with increasing τ. Therefore, the LSI will delay the CODP to acquire a considerable scale effect.Figure 4
Influence ofτ on CODP.Ifτ<0.05, the customized degree Z3 decreases with increasing τ, which is contrary to τ≥0.05 (Figure 5). Z3 initially decreases and then increases. This observation differs from the observation of previous studies that an increasing τ corresponds to a more obvious scale effect and smaller cost of mass service. Furthermore, previous studies show that the customized degree is equal to the ratio of the customization service cost and total cost (mass service cost plus customization service cost). Therefore, the customized degree should always increase.Figure 5
Influence ofτ on Z3 (customized degree).Figure 6
Influence ofτ on Z (comprehensive performance of LSSC).Our results differ from those of previous studies that did not include the CODP. The LSI will delay the CODP to acquire the scale effect, which leads to a significant increase in mass service cost and a decrease in customized degree (Figure4). Furthermore, the CODP becomes stable when τ increases to a certain value. A greater τ corresponds to a more obvious scale effect. The increasing extent of the mass service cost becomes smaller and the customized degree becomes larger. The mass service cost decreases and the CODP moves backward with increasing τ. The customized degree initially decreases and then increases.Ifτ≺0.025, Z decreases with increasing τ, which is contrary to τ≥0.025 (Figure 6). The comprehensive performance of LSSC initially decreases and then increases.Further analysis shows that whenZ is greater than 0.403 (τ=0.0615), Z increases with increasing τ; thus, the LSI prefers the FLSP with a bigger τ. When Z is less than 0.403 (τ=0.0615), two τ’s equal one performance. The LSI still prefers the FLSP with a bigger τ because a bigger τ leads to smaller total costs.


### 5.3. Influence of Order Difference Tolerance Coefficientω on Order Allocation

To analyze the influence ofω on the related factors of order allocation (Z2,CODP,Z3,Z), we equate ω to 0.4, 0.45, 0.5, 0.55, 0.6, 0.65, 0.7, 0.8, 0.9, and 1. The data is shown in Table 7. Figures 7 to 10 are plotted by using the data from Table 7.Table 7
Influence ofω on order allocation.
ω
Z
2
Z
3
CODP
Z
0.4
17472.020
0.522
5
0.391
0.45
17472.020
0.522
5
0.391
0.5
17499.530
0.613
5
0.471
0.55
17499.530
0.613
5
0.471
0.6
17499.530
0.613
5
0.471
0.65
17640.000
0.681
4
0.561
0.7
17640.000
0.681
4
0.561
0.8
17640.000
0.681
4
0.561
0.9
17640.000
0.681
4
0.561
1
17640.000
0.681
4
0.561Figure 7
Influence ofω on Z2 (total cost of LSI).The order difference tolerance coefficientω has a certain influence on the total cost of the LSI (Figure 7). The total cost of the LSI increases with increasing ω. The order difference also increases with increasing ω. The CODP should be set before K0=min{K1,K2,…,KM}. The CODP may move forward, thus increasing the procedure of customization service and total cost. Therefore, the LSI should choose the customer with low order difference and maintain the CODP in the same point to reduce total cost.Whenω is greater than 0.65, the total cost is stable at 17,640 (Figure 7). Increasing ω with no limitations is unnecessary because of the upper limit (ω=0.65 in numerical analysis). Consequently, the LSI should set a reasonable upper limit for ω in the MCLS.Whenω is less than 0.65, the CODP is in the fifth procedure; when ω is larger than 0.65, the CODP is in the fourth procedure (Figure 8). These findings show that ω has a certain influence on CODP. A larger ω corresponds to a larger CODP range. The CODP moves forward with increasing ω to a certain value.Figure 8
Influence ofω on CODP.Figure9 demonstrates that ω affects Z3. Z3 increases with increasing ω. The CODP moves to the second procedure, thus increasing the procedure of customization service. Consequently, the customized degree will also increase. When ω is greater than 0.65, increasing ω with no limitation is unnecessary because of the upper limit. Thus, the LSI should set a reasonable upper limit for ω in the MCLS.Figure 9
Influence ofω on Z3 (customized degree).Figure 10
Influence ofω on Z (comprehensive performance of LSSC).Figure10 indicates that Z increases steadily with increasing ω and becomes stable at 0.561. The comprehensive performance of LSSC shows a step growth instead of a linear growth with increasing ω. This observation is attributed to the forward movement of the CODP to the second procedure. Variables ω and Z have ranges; however, the CODP cannot move forward indefinitely because it must move after the second procedure, which determines the maximum value of ω; ω is not a constraint if it is greater than the maximum value.


### 5.4. Influence of Relationship Cost Coefficientδ on Order Allocation

To analyze the influence ofδ on the related factors of order allocation (Z2,CODP,Z3,Z), we equated δ to 0.15, 0.2, 0.25, 0.3, 0.35, 0.4, and 0.45. The data are shown in Table 8. Figures 11 to 14 are plotted by using the data from Table 8.Table 8
Influence ofδ on order allocation.
δ
Z
2
Z
3
CODP
Z
0.15
15400.000
0.575
5
0.421
0.2
17472.020
0.522
5
0.391
0.25
17990.84
0.570
5
0.413
0.3
18887.720
0.652
4
0.497
0.35
19845.000
0.668
4
0.465
0.4
20580.000
0.663
4
0.477
0.45
21002.100
0.685
4
0.494Figure 11
Influence ofδ on Z2 (total cost of LSI).Figure11 shows that δ has a significant influence on the total cost of the LSI. A larger δ corresponds to a larger Z2. This conclusion is consistent with the conclusion of Liu et al. [1]. A larger δ corresponds to a greater Z2 range, thus increasing the total cost.Figure12 indicates that δ has a certain influence on the CODP. The CODP moves forward to the second procedure with increasing δ. A larger δ corresponds to a greater Z2 range, which increases the procedure of the customization service and decreases the procedure of the mass service. Therefore, the CODP will move forward.Figure 12
Influence ofδ on CODP.Whenδ<0.2, Z3 decreases along with increasing δ, which is contrary to δ≥0.2 (Figure 13). The relationship cost coefficient δ influences the customized degree positively when δ is greater than a certain value.Figure 13
Influence ofδ on Z3 (customized degree).Figure 14
Influence ofδ on Z (comprehensive performance of LSSC).Further analysis shows that if the customized degree is greater than 0.575,Z3 increases with increasing δ (Figure 13). Thus, the LSI prefers to increase the value of δ. However, when the customized degree is less than 0.575, two δs will be present. For this case, the LSI will prefer the smaller δ. For example, when the customized degree is 0.55, two δ will be present: C1 (0.174,0.55) and C2 (0.22,0.55) (Figure 13). In this case, the LSI will choose C1.Figure14 shows that Z decreases before increasing with increasing δ. The breaking point is when δ=0.2. This observation is different from previous studies that indicate that a greater δ corresponds to a looser constraint. Hence, Z will increase.This observation is attributed to the two parts ofZ: customized degree and satisfaction. Figure 13 shows that if δ is less than 0.2, the customized degree decreases with increasing δ. Consequently, Z decreases. When δ is more than 0.2, Z increases with increasing customized degree (Figure 13).Figure14 shows that Z only exhibits a positive correlation with δ when Z is greater than 0.421. When Z is less than 0.421, Z initially decreases and then increases. When Z is less than 0.421, two δ’s are present. For this case, the LSI will prefer the smaller δ. For example, when Z equals 0.4, two δ’s exist: D1 (0.18, 0.4) and D2 (0.22, 0.4). In this case, the LSI will choose D1.


### 5.5. Relationship of Order Allocation Results and Related Factors

In this section, the relationship of order allocation results (total costZ2, customized degree Z3, CODP, and comprehensive performance Z) and related factors (scale effect coefficient τ, order difference tolerance coefficient ω, and relationship cost coefficient δ) is summarized.


## 5.1. Model Solution

Figures1 and 2 show that each customer (A, B, and C) has an order and needs the LSI Q to provide transportation services. Q assigns the three orders to the FLSPs (a, b, c, d, and e). The demand of each customer is 60, 100, and 80. The number of service procedures of each customer is 8, 7, and 7. K1=6, K2=5, K3=5, K0=min(K1,K2,K3)=5, τ=0.05, ω=0.4, and δ=20%. The values of other parameters are shown in Tables 2 and 3.Table 2
Parameter data (1).
[
a
i
j
k
,
b
i
j
k
]
[
a
i
j
k
1
,
b
i
j
k
1
]
[
a
i
j
k
2
,
b
i
j
k
2
]
i
=
1
[
40
,
80
]
[
20
,
40
]
i
=
2
[
30
,
70
]
[
20
,
35
]
i
=
3
[
35
,
90
]
[
25
,
40
]
i
=
4
[
30
,
60
]
[
25
,
50
]
i
=
5
[
30
,
80
]
[
10
,
40
]
c
i
j
k
c
i
j
k
1
c
i
j
k
2
i
=
1
10
-
0.02
x
i
j
k
15
i
=
2
11
-
0.02
x
i
j
k
18
i
=
3
9
-
0.02
x
i
j
k
16
i
=
4
12
-
0.02
x
i
j
k
20
i
=
5
14
-
0.02
x
i
j
k
22
ρ
i
j
k
ρ
i
j
k
1
ρ
i
j
k
2
i
=
1
0.2
0.3
i
=
2
0.35
0.35
i
=
3
0.15
0.3
i
=
4
0.3
0.2
i
=
5
0.25
0.3
ε
i
ε
i
1
ε
i
2
i
=
1
0.6
0.4
i
=
2
0.5
0.5
i
=
3
0.7
0.3
i
=
4
0.4
0.6
i
=
5
0.35
0.65Table 3
Parameter data (3).
λ
i
j
j
=
1
j
=
2
j
=
3
i
=
1
0.3
0.3
0.4
i
=
2
0.2
0.4
0.4
i
=
3
0.35
0.35
0.3
i
=
4
0.25
0.35
0.4
i
=
5
0.3
0.4
0.3
β
j
0.3
0.4
0.3
D
j
60
100
80
K
j
6
5
5
L
j
8
7
7Figure 1
Structure of the LSSC.Figure 2
Logistics service procedure of customer A, B, and C. Note: The red circle denotes the procedure that starts the customization service of each customer.By using the solution method in Section4, Tables 4 and 5 show the result and performance of the order allocation, respectively.Table 4
Order allocation results.
Procedure
1–5
6–7
8
CustomerOrder valueFLSP
A
B
C
A
B
C
A
a
0
4.31
0
31.59
18.20
37.79
31.59
b
0
6.33
0
6.14
37.86
0.53
6.14
c
60
4.99
80
4.63
40
37.82
4.63
d
0
4.93
0
6.41
0.35
0
6.41
e
0
79.44
0
11.23
3.60
3.86
11.23Table 5
Order allocation performance.
Total cost of LSSC
17472.020
CODP
5
Customized degree
0.522
Z
0.391


## 5.2. Influence of Scale Effect Coefficientτ on Order Allocation

To analyze the influence ofτ on the related factors of order allocation (Z2,CODP,Z3,Z), we make τ equal to 0, 0.025, 0.05, 0.075, 0.1, 0.125, and 0.15. The data is shown in Table 6. Figures 3 to 6 are plotted with the data from Table 6.Table 6
Influence ofτ on order allocation.
τ
Z
2
Z
3
CODP
Z
0
20280
0.596
4
0.403
0.025
18960.000
0.524
5
0.388
0.05
17472.020
0.522
5
0.391
0.075
16320.000
0.580
5
0.417
0.1
15000.000
0.638
5
0.469
0.125
13680.000
0.646
5
0.472
0.15
12360.000
0.783
5
0.582Figure 3
Influence ofτ on Z2 (total cost of LSI).Figure3 shows that τ has significant influence on the total cost of LSI Z2. A greater τ corresponds to a smaller Z2. The scale effect of the cost is more obvious with increasing τ; thus, the total cost decreases.Figure4 shows that τ has a certain influence on CODP. If τ=0, the scale effect is nonexistent and the CODP is in the fourth procedure. With τ≥0.025, the CODP changes to the fifth procedure because the scale effect of cost is more obvious with increasing τ. Therefore, the LSI will delay the CODP to acquire a considerable scale effect.Figure 4
Influence ofτ on CODP.Ifτ<0.05, the customized degree Z3 decreases with increasing τ, which is contrary to τ≥0.05 (Figure 5). Z3 initially decreases and then increases. This observation differs from the observation of previous studies that an increasing τ corresponds to a more obvious scale effect and smaller cost of mass service. Furthermore, previous studies show that the customized degree is equal to the ratio of the customization service cost and total cost (mass service cost plus customization service cost). Therefore, the customized degree should always increase.Figure 5
Influence ofτ on Z3 (customized degree).Figure 6
Influence ofτ on Z (comprehensive performance of LSSC).Our results differ from those of previous studies that did not include the CODP. The LSI will delay the CODP to acquire the scale effect, which leads to a significant increase in mass service cost and a decrease in customized degree (Figure4). Furthermore, the CODP becomes stable when τ increases to a certain value. A greater τ corresponds to a more obvious scale effect. The increasing extent of the mass service cost becomes smaller and the customized degree becomes larger. The mass service cost decreases and the CODP moves backward with increasing τ. The customized degree initially decreases and then increases.Ifτ≺0.025, Z decreases with increasing τ, which is contrary to τ≥0.025 (Figure 6). The comprehensive performance of LSSC initially decreases and then increases.Further analysis shows that whenZ is greater than 0.403 (τ=0.0615), Z increases with increasing τ; thus, the LSI prefers the FLSP with a bigger τ. When Z is less than 0.403 (τ=0.0615), two τ’s equal one performance. The LSI still prefers the FLSP with a bigger τ because a bigger τ leads to smaller total costs.


## 5.3. Influence of Order Difference Tolerance Coefficientω on Order Allocation

To analyze the influence ofω on the related factors of order allocation (Z2,CODP,Z3,Z), we equate ω to 0.4, 0.45, 0.5, 0.55, 0.6, 0.65, 0.7, 0.8, 0.9, and 1. The data is shown in Table 7. Figures 7 to 10 are plotted by using the data from Table 7.Table 7
Influence ofω on order allocation.
ω
Z
2
Z
3
CODP
Z
0.4
17472.020
0.522
5
0.391
0.45
17472.020
0.522
5
0.391
0.5
17499.530
0.613
5
0.471
0.55
17499.530
0.613
5
0.471
0.6
17499.530
0.613
5
0.471
0.65
17640.000
0.681
4
0.561
0.7
17640.000
0.681
4
0.561
0.8
17640.000
0.681
4
0.561
0.9
17640.000
0.681
4
0.561
1
17640.000
0.681
4
0.561Figure 7
Influence ofω on Z2 (total cost of LSI).The order difference tolerance coefficientω has a certain influence on the total cost of the LSI (Figure 7). The total cost of the LSI increases with increasing ω. The order difference also increases with increasing ω. The CODP should be set before K0=min{K1,K2,…,KM}. The CODP may move forward, thus increasing the procedure of customization service and total cost. Therefore, the LSI should choose the customer with low order difference and maintain the CODP in the same point to reduce total cost.Whenω is greater than 0.65, the total cost is stable at 17,640 (Figure 7). Increasing ω with no limitations is unnecessary because of the upper limit (ω=0.65 in numerical analysis). Consequently, the LSI should set a reasonable upper limit for ω in the MCLS.Whenω is less than 0.65, the CODP is in the fifth procedure; when ω is larger than 0.65, the CODP is in the fourth procedure (Figure 8). These findings show that ω has a certain influence on CODP. A larger ω corresponds to a larger CODP range. The CODP moves forward with increasing ω to a certain value.Figure 8
Influence ofω on CODP.Figure9 demonstrates that ω affects Z3. Z3 increases with increasing ω. The CODP moves to the second procedure, thus increasing the procedure of customization service. Consequently, the customized degree will also increase. When ω is greater than 0.65, increasing ω with no limitation is unnecessary because of the upper limit. Thus, the LSI should set a reasonable upper limit for ω in the MCLS.Figure 9
Influence ofω on Z3 (customized degree).Figure 10
Influence ofω on Z (comprehensive performance of LSSC).Figure10 indicates that Z increases steadily with increasing ω and becomes stable at 0.561. The comprehensive performance of LSSC shows a step growth instead of a linear growth with increasing ω. This observation is attributed to the forward movement of the CODP to the second procedure. Variables ω and Z have ranges; however, the CODP cannot move forward indefinitely because it must move after the second procedure, which determines the maximum value of ω; ω is not a constraint if it is greater than the maximum value.


## 5.4. Influence of Relationship Cost Coefficientδ on Order Allocation

To analyze the influence ofδ on the related factors of order allocation (Z2,CODP,Z3,Z), we equated δ to 0.15, 0.2, 0.25, 0.3, 0.35, 0.4, and 0.45. The data are shown in Table 8. Figures 11 to 14 are plotted by using the data from Table 8.Table 8
Influence ofδ on order allocation.
δ
Z
2
Z
3
CODP
Z
0.15
15400.000
0.575
5
0.421
0.2
17472.020
0.522
5
0.391
0.25
17990.84
0.570
5
0.413
0.3
18887.720
0.652
4
0.497
0.35
19845.000
0.668
4
0.465
0.4
20580.000
0.663
4
0.477
0.45
21002.100
0.685
4
0.494Figure 11
Influence ofδ on Z2 (total cost of LSI).Figure11 shows that δ has a significant influence on the total cost of the LSI. A larger δ corresponds to a larger Z2. This conclusion is consistent with the conclusion of Liu et al. [1]. A larger δ corresponds to a greater Z2 range, thus increasing the total cost.Figure12 indicates that δ has a certain influence on the CODP. The CODP moves forward to the second procedure with increasing δ. A larger δ corresponds to a greater Z2 range, which increases the procedure of the customization service and decreases the procedure of the mass service. Therefore, the CODP will move forward.Figure 12
Influence ofδ on CODP.Whenδ<0.2, Z3 decreases along with increasing δ, which is contrary to δ≥0.2 (Figure 13). The relationship cost coefficient δ influences the customized degree positively when δ is greater than a certain value.Figure 13
Influence ofδ on Z3 (customized degree).Figure 14
Influence ofδ on Z (comprehensive performance of LSSC).Further analysis shows that if the customized degree is greater than 0.575,Z3 increases with increasing δ (Figure 13). Thus, the LSI prefers to increase the value of δ. However, when the customized degree is less than 0.575, two δs will be present. For this case, the LSI will prefer the smaller δ. For example, when the customized degree is 0.55, two δ will be present: C1 (0.174,0.55) and C2 (0.22,0.55) (Figure 13). In this case, the LSI will choose C1.Figure14 shows that Z decreases before increasing with increasing δ. The breaking point is when δ=0.2. This observation is different from previous studies that indicate that a greater δ corresponds to a looser constraint. Hence, Z will increase.This observation is attributed to the two parts ofZ: customized degree and satisfaction. Figure 13 shows that if δ is less than 0.2, the customized degree decreases with increasing δ. Consequently, Z decreases. When δ is more than 0.2, Z increases with increasing customized degree (Figure 13).Figure14 shows that Z only exhibits a positive correlation with δ when Z is greater than 0.421. When Z is less than 0.421, Z initially decreases and then increases. When Z is less than 0.421, two δ’s are present. For this case, the LSI will prefer the smaller δ. For example, when Z equals 0.4, two δ’s exist: D1 (0.18, 0.4) and D2 (0.22, 0.4). In this case, the LSI will choose D1.


## 5.5. Relationship of Order Allocation Results and Related Factors

In this section, the relationship of order allocation results (total costZ2, customized degree Z3, CODP, and comprehensive performance Z) and related factors (scale effect coefficient τ, order difference tolerance coefficient ω, and relationship cost coefficient δ) is summarized.


## 6. Main Conclusions and Management Insights

In this section, we elaborate the main conclusions and implications for researchers. We also explore the management insights for LSI and FLSP and provide some suggestions.


### 6.1. Main Conclusions

On the basis of the analysis on Section5, we obtained the following conclusions concerning order allocation under an MCLS environment.(1)
The total cost of LSIZ2 decreases with increasing scale effect coefficient τ and increases with increasing order difference tolerance coefficient ω and relationship cost coefficient δ. However, the influence of ω on Z2 is finite. Z2 becomes stable when ω increases to a certain value.(2)
The CODP moves back to the customers with increasing scale effect coefficientτ and moves forward to the second procedure with increasing order difference tolerance coefficient ω and relationship cost coefficient δ. LSI can move the CODP forward and increase the customized degree by increasing the value of δ. When τ is greater, the LSI prefers moving the CODP backward and acquiring the scale effect.(3)
The customized degreeZ3 initially decreases before increasing with increasing scale effect coefficient τ and relationship cost coefficient δ. When Z3 is greater than a certain value, τ and δ have a positive correlation with Z3. ω always has a positive influence on Z3. When ω increases to a certain value, Z3 tends to be stable.(4)
The comprehensive performance of order allocationZ initially decreases before increasing with increasing scale effect coefficient τ and relationship cost coefficient δ. When Z is greater than a certain value, τ and δ have a positive correlation with Z; ω always has a positive effect on Z. When ω increases to a certain value, Z tends to be a stable.


### 6.2. Implications for Researchers

In line with academic research, this paper elaborates on the management of order allocation under the MCLS environment and the influences caused by allocation results. This paper provides a theoretical basis for the study of control methods about the performance of order allocation in MCLS. For instance, the type of order allocation method needs to be determined in choosing the FLSP. How will the comprehensive performance of the LSSC improve with minimum cost? How will the contradiction between cost and customization be solved? What factors have more influence on the comprehensive performance of the LSSC? This study provides the necessary theoretical basis for further research in the application and scheduling of the theoretical and empirical studies of LSSC under the MCLS environment.


### 6.3. Implications for Managers

LSI and FLSP managers should understand the factors that need focus and establish long-term strategic partnership for better cooperation.The LSI should be careful in order allocation and choose the FLSP with a greater scale effect. The order difference tolerance factor should not be set extremely high. Excessive order difference tolerance results in a higher comprehensive performance of order allocation; however, excessive order difference tolerance can also lead to excessive total LSI cost. The order difference tolerance has a reasonable upper limit. Beyond this upper limit, the comprehensive performance of order allocation and the degree of customized experience initially decrease before increasing with increasing relationship cost coefficient. Therefore, a bigger relationship cost coefficient of the LSI does not indicate a better setup; instead, we should consider both the customized degree of the target and comprehensive performance.For FLSP, the scale effect coefficient is important to enterprises for order allocation. Thus, enterprises need to improve large-scale operational capabilities continuously. When the comprehensive performance is high, FLSPs should request for cost compensation from the LSI by increasing the relationship cost coefficient, which improves FLSP satisfaction.


## 6.1. Main Conclusions

On the basis of the analysis on Section5, we obtained the following conclusions concerning order allocation under an MCLS environment.(1)
The total cost of LSIZ2 decreases with increasing scale effect coefficient τ and increases with increasing order difference tolerance coefficient ω and relationship cost coefficient δ. However, the influence of ω on Z2 is finite. Z2 becomes stable when ω increases to a certain value.(2)
The CODP moves back to the customers with increasing scale effect coefficientτ and moves forward to the second procedure with increasing order difference tolerance coefficient ω and relationship cost coefficient δ. LSI can move the CODP forward and increase the customized degree by increasing the value of δ. When τ is greater, the LSI prefers moving the CODP backward and acquiring the scale effect.(3)
The customized degreeZ3 initially decreases before increasing with increasing scale effect coefficient τ and relationship cost coefficient δ. When Z3 is greater than a certain value, τ and δ have a positive correlation with Z3. ω always has a positive influence on Z3. When ω increases to a certain value, Z3 tends to be stable.(4)
The comprehensive performance of order allocationZ initially decreases before increasing with increasing scale effect coefficient τ and relationship cost coefficient δ. When Z is greater than a certain value, τ and δ have a positive correlation with Z; ω always has a positive effect on Z. When ω increases to a certain value, Z tends to be a stable.


## 6.2. Implications for Researchers

In line with academic research, this paper elaborates on the management of order allocation under the MCLS environment and the influences caused by allocation results. This paper provides a theoretical basis for the study of control methods about the performance of order allocation in MCLS. For instance, the type of order allocation method needs to be determined in choosing the FLSP. How will the comprehensive performance of the LSSC improve with minimum cost? How will the contradiction between cost and customization be solved? What factors have more influence on the comprehensive performance of the LSSC? This study provides the necessary theoretical basis for further research in the application and scheduling of the theoretical and empirical studies of LSSC under the MCLS environment.


## 6.3. Implications for Managers

LSI and FLSP managers should understand the factors that need focus and establish long-term strategic partnership for better cooperation.The LSI should be careful in order allocation and choose the FLSP with a greater scale effect. The order difference tolerance factor should not be set extremely high. Excessive order difference tolerance results in a higher comprehensive performance of order allocation; however, excessive order difference tolerance can also lead to excessive total LSI cost. The order difference tolerance has a reasonable upper limit. Beyond this upper limit, the comprehensive performance of order allocation and the degree of customized experience initially decrease before increasing with increasing relationship cost coefficient. Therefore, a bigger relationship cost coefficient of the LSI does not indicate a better setup; instead, we should consider both the customized degree of the target and comprehensive performance.For FLSP, the scale effect coefficient is important to enterprises for order allocation. Thus, enterprises need to improve large-scale operational capabilities continuously. When the comprehensive performance is high, FLSPs should request for cost compensation from the LSI by increasing the relationship cost coefficient, which improves FLSP satisfaction.


## 7. Research Limitations and Future Work

Under increasing customer demand for specialized and customized logistics services, the competitiveness of the LSSC depends on its ability to meet the need for customized requirements with the operating expense of mass service. The satisfaction of FLSPs will be maximized by minimizing the total cost of the LSI and maximizing the customized degree for customer demand. This paper establishes an order allocation model for the LSSC that is constrained by meeting the demand, CODP, and order difference tolerance coefficient. Numerical analysis is conducted with Lingo 12 software. The influences of the scale effect coefficient, order difference tolerance coefficient, and relationship cost coefficient on the order allocation and comprehensive performance of the LSSC are discussed. Management insights are also proposed.However, the order allocation model also has some disadvantages. For instance, we assumed that each service process in order allocation requires only one service capability, such as transport service capability. In many cases, logistics service is the combination of different services including transportation, storage, loading, and unloading; each process may also be supported by a variety of capacities. The equivalent problem of each capacity needs to be considered in detail in the future. We also assumed that customer demand and service capacity are stable; however, customer demand may be random in actuality. Order allocation with random customer demand is an important direction for future studies. Furthermore, we assumed that the members of the LSSC have long-term cooperation, transparency, and mutual trust. In reality, we need to consider the game behavior among members because of its influence on order allocation. These problems need to be explored in detail in future research.



---

*Source: 957260-2013-11-12.xml*