Next Article in Journal
Fuzzy Differential Subordination for Meromorphic Function Associated with the Hadamard Product
Next Article in Special Issue
A Game—Theoretic Model for a Stochastic Linear Quadratic Tracking Problem
Previous Article in Journal
The Existence, Uniqueness, and Multiplicity of Solutions for Two Fractional Nonlocal Equations
Previous Article in Special Issue
Bargaining-Based Profit Allocation Model for Fixed Return Investment Water-Saving Management Contract
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Joint Location–Allocation–Inventory Spare Part Optimization Model for Base-Level Support System with Uncertain Demands

1
School of Reliability and System Engineering, Beihang University, Beijing 100191, China
2
Science and Technology on Reliability and Environmental Engineering Laboratory, Beihang University, Beijing 100191, China
3
School of Aeronautic Science and Engineering, Beihang University, Beijing 100191, China
*
Author to whom correspondence should be addressed.
Axioms 2023, 12(1), 46; https://doi.org/10.3390/axioms12010046
Submission received: 23 November 2022 / Revised: 25 December 2022 / Accepted: 26 December 2022 / Published: 1 January 2023
(This article belongs to the Special Issue Advances in Uncertain Optimization and Applications)

Abstract

:
This paper copes with a joint Location-Allocation-Inventory problem in a three-echelon base-level spare part support system with epistemic uncertainty in uncertain demands of bases. The aim of the paper is to propose an optimization model under the uncertainty theory to minimize the total cost, which integrates crucial characterizations of the inventory control decisions and the location-allocation scheme arrangement under a periodic review order-up-to-S (T, S) policy. Uncertainty theory is introduced in this paper to characterize epistemic uncertainty, where demands are treated as uncertain variables and stockout loss is represented by value-at-risk in uncertain measurement. To solve the original uncertain optimization model, an equivalent deterministic model is derived and addressed by an improved bilevel genetic algorithm. Moreover, the proposed models and algorithm are encoded into numerical examples for supply chain programming. The results highlight the applicability of the model and the algorithm’s effectiveness in approaching the optimal solution compared with traditional genetic algorithm. Sensitivity analyses are further made for the impacts of review time and inventory capacity on different cost components.

1. Introduction

Base-level support systems provide spare parts storage, replacement, and other services, which play an essential role in the current maintenance strategy of sophisticated military and engineering systems [1,2]. Moreover, a collection of industries has witnessed that inventory control of spare parts and physical distribution of supply sites consume more than half of the total cost statistically [2]. In this sense, it is of great value for base-level support systems to explore how to appropriately manage the location of supply sites, inventory strategies, and allocation relationships.
The Location-Allocation-Inventory model (LAIM) is a sort of joint model aiming to simultaneously address the problem of location, allocation, and inventory control. Recent developments have revealed the significance of building such models. One of the earliest studies is proposed by Yao et al. [3]. Considering customer demands and safety storage stock, they established a model to shrink the supply chain costs by specifying the number of warehouses along with their corresponding locations, the allocation quantity for customers, and the control of inventory levels. Likewise, in 2012, Tsao et al. [4] presented an integrated facility-inventory allocation network to maximize overhead reduction without breaking the regulated area coverage restraints. Dai et al. [5] attempted to make a tradeoff between expenditure involved and stipulating stock capacity with carbon emissions for a perishable product location–inventory supply chain. In addition, this kind of joint model can also be utilized to improve the efficiency of municipal solid waste collection systems [6]. It can be seen that LAIMs have been applied to a variety of issues in industries and environmental fields with good enforcement performance. However, these models can’t be directly embedded in support systems due to the lack of consideration on the characteristics of support systems, such as supply availability, service level, and balanced dispatching regulations [7,8,9].
Uncertainty quantification of supply chain variables is also a non-negligible issue, which has a significant impact on the results of supply chain optimization. In the past decades, probability theory has been firstly employed to cope with such problems. For example, Asian and Nie [10] used normal distribution to describe demand’s uncertainty with consideration of market volatility. Islam et al. [11] delineated Poisson distributed demand and exponentially distributed supplier capacity in a three-layer inventory-allocation management system merging. Errico et al. [12] assumed service times to be mutually independent discrete triangular distributions in a vehicle routing problem. Nevertheless, probability theory is not always appropriate to disentangle indeterminate scenarios of spare part provision, especially for military resource supplements, because, for armament systems, the burstiness and variability of military missions often cause the historical knowledge unable to completely characterize the demand information. As for the abovementioned scene with sparse useful observational data, frequency probabilistic methods are limited by the law of the large numbers, namely the frequency tends to the probability when the number of independent repeated trials tends to infinity [13,14]. Moreover, it is epistemic uncertainty rather than aleatory uncertainty that becomes the principal uncertainty under such sparse or no data conditions. Therefore, two mathematical measures are introduced as alternatives to specify epistemic uncertainty.
One train of thought originating from fuzzy theory has been explored a lot for a few decades [15,16,17,18]. Nonetheless, the puzzle that fuzzy theory is not compliant with the duality axiom, has confused supply chain decision-makers for a long time. The other is the uncertainty theory built by Liu [19], which is regarded as a breakthrough in handling epistemic uncertainty problems that predicate normality, duality, subadditivity, and product axioms. In the wake of the theory’s development, it has been gradually discovered to be an excellent way to tackle the epistemic uncertainty in supply chains. For instance, Sheng et al. [20] referred to an infinite-horizon production–inventory optimal control problem with regulated production rates which is transformed into the value function illustrated by uncertainty theory owing to the scarcity of affordable experiment samples. Asim et al. [21] demonstrated an integrated multi-echelon multi-item production-transportation close-loop network, where cost, demand, and capacity are characteristic as uncertain variables to vanquish the inapplicability of probability or fuzzy set theory for the model with belief degree. Shen [22] pointed out that many environmental and social emergencies in supply chains incur the demand, cost, and capacity acting as uncertain variables. Hosseini and Pishvaee [23] put forward an α -maximum capacity path (UMCP) problem with inaccurate transportation reliability and operational capacity of links appearing as ’belief degrees’ caused by deficient information and computing simplification. In their work, links’ capacities are deemed as uncertain variables, and the stability analysis is launched in the framework of uncertainty theory. To this end, they utilized the uncertainty theory to establish an uncertain purveyor selection model concerning the elements of cost, environmental impact, and social benefits evaluated by uncertainty theory. Therefore, in this paper, uncertainty theory is expropriated to quantify the uncertainties of demand caused by the insufficiency of historical information in the programming stage and the shortage cost loss deriving from exigencies of long-distance delivery. Specifically, we suppose demands as normal-distributed uncertain variables mingled into the supportability chance constraints, while the penalty cost is weighted by the uncertain value-at-risk.
Since LAIMs need to concurrently determine the depot locations, the allocation relationships, and inventory parameters, the problem itself tends to be high non-linearity and NP-hardness with huge feasible space. Methodically, the kind of NP-hardness difficulty yields the rise of plenty of heuristic and metaheuristic algorithms which do well in scaling down the CPU time and model size [24]. Up to now, a handful of latest works have proponed some algorithms for addressing LAIMs. For example, Tirkolaee et al. [25] developed a Self-Learning Particle Swarm Optimization (SLPSO) algorithm for solving a multi-echelon capacitated LAIM, which tackled the PSO algorithm’s deficiency of lacking intelligent particles. Mousavi et al. [26] made a comparison of a Modified Genetic Algorithm (MGA) and a Particle Swarm Optimization (PSO) for gaining results of a random two-echelon LAIM. The gained results exemplified that the MGA with Taguchi parameter adjustment approach outperforms the PSO accordant with the fitness functions and computing time. Sahraei and Samouei [27] analyzed the distinct performance of genetic and electromagnetic meta-heuristic algorithms for a bi-level scenario-based LAIM and indicated the genetic algorithm is more effective. It can be seen that genetic algorithm and particle swarm optimization algorithm are both feasible approaches for solving the Location-Allocation-Inventory problems, and the introduction of parameter adjustment methods along with bilevel optimization methodology play positive roles in algorithm convergence and optimal solution acquisition. Consequently, we deploy an improved bilevel genetic algorithm to accomplish the model-solving, attempting to improve operation accuracy in contrast with the traditional genetic algorithm.
The highlights of this work include:
(i) Different from the previous literature, this work develops a joint Location-Allocation-Inventory model (LAIM) for a base-level spare part supply problem, which mingles with supportability indexes, namely, service level and supply availability.
(ii) Unlike the common way of demand characterization, the proposed model considers the vein with a lack of historical data about demand information, and thus the demands in the network were characterized by the uncertainty theory with value-at-risk introduced for stockout loss evaluation.
(iii) As for solving the uncertain optimization model, an equivalent deterministic model was derived with an improved bilevel genetic algorithm proposed for accelerating the optimal value acquirement.
The remainder of this paper is organized as follows: Section 2 presents the preliminaries on uncertainty theory. Section 3 is dedicated to describing the problem with an in-depth exploration of the provided three-echelon support system, and according to the hypothesis and analysis, a spare parts supply optimization model is developed with shortage events existing. Then, Section 4 emphasizes employing an improved bilevel genetic algorithm to help designate appropriate transportation and storage tactics. Section 5 shows an application of the proposed model in a specific case, a comparison between the proposed algorithm and the traditional GA, and a sensitivity analysis of various cost components. Concluding remarks are drawn along with a discussion for future research in Section 6.

2. Preliminaries on Uncertainty Theory

For the sake of appropriately forecasting the demand quantity based on subjective experience, we adopt the uncertainty theory which is a novel axiomatic mathematical theory proposed by Liu for resolving epistemic uncertainty [19]. Below are some definitions and theorems related to uncertain assessment.
Definition 1 
(Uncertain variable). An uncertain variable ξ is a function from an uncertainty space ( Γ , L , M ) to the set of real numbers such that ξ B is an event for any Borel set B of real numbers.
Definition 2 
(Uncertain distribution). The uncertainty distribution Φ of an uncertain variable ξ is defined by
Φ ( α ) = M { ξ x }
for any real number x
Definition 3 
(Normal uncertainty distribution). An uncertain variable ξ is called a normal variable if it has a normal uncertainty distribution
Φ ( x ) = 1 + e x p π ( e x ) 3 σ 1
denoted by N ( e , σ ) ,where e and σ are real numbers with σ > 0 .
Definition 4 
(Inverse uncertain distribution). Let ξ be an uncertain variable with regular uncertainty distribution Φ ( x ) . Then, the inverse function Φ 1 ( α ) is called the inverse uncertainty distribution of ξ. A function Φ 1 is an inverse uncertainty distribution of an uncertain variable ξ if and only if
M { ξ Φ 1 ( α ) } = α
for all α [ 0 , 1 ] .
Definition 5 
(Inverse normal uncertain distribution). The inverse uncertainty distribution of normal uncertain variable N ( e , σ ) is
Φ 1 ( α ) = e + σ 3 π ln α 1 α .
Definition 6 
(Inverse normal uncertain distribution). Let ξ be an uncertain variables. Then, the expected value of ξ is defined by
E [ ξ ] = 0 + M { ξ x } d x 0 M { ξ x } d x
provided that at least one of the two integrals is finite.
Definition 7 
(Inverse normal uncertain distribution). Assume that a system contains uncertain factors ξ 1 , ξ 2 , , ξ n and has a loss function f. Then, the value-at-risk is defined as
VaR ( α ) = s u p { x | M { f ( ξ 1 , ξ 2 , , ξ n ) x } α } .
Note that VaR (α) represents the maximum possible loss when α percent of the right tail distribution is ignored.
Theorem 1 
(Inverse uncertain distribution operational law). Let ξ 1 , ξ 2 , , ξ n be independent uncertain variables with regular uncertainty distributions Φ 1 , Φ 2 , , Φ n , respectively. If f ( x 1 , x 2 , , x n ) is continuous, strictly increasing with respect to x 1 , x 2 , , x m and strictly decreasing with respect to x m + 1 , x m + 2 , , x n , then ξ = f ( ξ 1 , ξ 2 , , ξ n ) has an inverse uncertainty distribution
Ψ ( α ) = f ( Φ 1 1 ( α ) , Φ 2 1 ( α ) , , Φ m 1 ( α ) , Φ m + 1 1 ( 1 α ) , Φ m + 2 1 ( 1 α ) , , Φ n 1 ( 1 α ) ) .
Theorem 2 
(Normal uncertain distribution operational law). Let ξ 1 and ξ 2 be independent normal uncertain variables N ( e 1 , σ 1 ) and N ( e 2 , σ 2 ) , respectively. Then, the sum ξ 1 + ξ 2 is also a normal uncertain variable N ( e 1 + e 2 , σ 1 + σ 2 ) , i.e.,
N ( e 1 , σ 1 ) + N ( e 2 , σ 2 ) = N ( e 1 + e 2 , σ 1 + σ 2 )

3. Problem Description and Model Formulation

In this section, we first collate the problem of a three-level supply chain network for the base-level support system. With the introduction of the assumptions and notations, an uncertain optimization model is formulated for the problem.

3.1. Problem Statement

This paper specializes in providing a three-echelon Location-Allocation-Inventory network for supporting a kind of base-level spare part. As Figure 1 demonstrated, the network is made up of one supplier and several bases scattered within a relatively concentrated region. Supplier is responsible for producing and offering spare parts required by bases, which is usually unique and relatively far from the bases in the support system. The bases are responsible for replacing spare parts for the machine [28]. Therefore, to realize the timely support of spare parts, it is essential for the bases to have a certain amount of storage capacity. Each base in the network has a depot, which is of course the ideal situation. However, due to the constraints of resources, only a certain number of depots can be allowed to be built. In this paper, we identify the number of depots as a given number n with the consideration of depot-covered areas and affordable support resources. As a result, the bases can be divided into two categories: bases without depot, and bases with depot. The latter should undertake the responsibility for launching periodical orders from the supplier and responding to real-time demand from bases including itself. Note that each base accepts spare parts from a single depot, each depot can serve more than one base.
In respect of demand estimation, we confront the situation that historical demand data is scarce or untrustworthy with epistemic uncertainties domaining the parameter characterization. Because of this, we introduce uncertainty theory to mathematically describe the demand quantity and cost constituents.
In the articulated supply chain for base-level spare parts provision, the coordinates of bases are predefined, while depot locations and the depot-base supply relationships need to be determined. In addition, we suppose the selected centralized depots are amenable to the (T, S) policy defined as the periodic review order-up-to-S [29]. As Figure 2 demonstrated, the floating demand of bases (Figure 2a) transfers to upper-layer centralized depots, which brings about the descent of inventory level, whereas the periodical ordering and the followed inventory replenishment provoke the frequent rebounds of inventory level. It is worth noting that there is a certain length of time (the so-called lead time) between the ordering and inventory replenishment stages, which emanates from the manufacturing process and long-distance delivery [30]. For simplification, we set it to be a constant value. Figure 2 also underlines one scene where the net inventory level is inferior to zero, and the scene can be explained as the trigger of a risk shortage. From the abovementioned management of Location-Allocation-Inventory activities, we intend to minimize the total supply cost, which covers the amortized expenditure arising from depot setup/renovation along with safeguard, delivery expense, inventory holding cost, stockout risk loss, and order outlay.
Compared with the provided services of common merchandise, spare parts support systems are bound to meet some supportability standards. In this article, we pick up service level and supply availability as supportability evaluation indexes. Meanwhile, reasonable allotment of safeguard assets is amalgamated to the equilibrium supply relationship of depots.

3.2. Assumptions

In accordance with the problem background, the underlying assumptions are summarized as follows:
  • Lead time is regarded as a constant/
  • Bases have no storage capacity initially but have the probability to construct a depot and each of them should be supplied by only one depot including itself.
  • The difference value between the number of service bases among centralized depots cannot exceed one for equilibrating provision;
  • (T, S) inventory policy is adopted for each depot with a review period regulated within a certain range.
  • Base demand is independent and normally distributed uncertain variable.
  • Unmet demand in a period is reckoned as backlogging loss which needs to be replenishment in the next period.
  • The total number of depots in the network is fixed.
  • Depots are considered as distributing centers that place orders to the supplier and satisfy the real-time demand from depots.

3.3. Abbreviations and Notations

Abbreviations
LAIM: Location–Allocation–Inventory Model;
ULAIM: Uncertain Location–Allocation–Inventory Model;
EBO: Expected backorder (stockout order);
GA: Genetic Algorithm;
Parameters
I: set of bases which are also viewed as candidate depots;
b i : fix the section of unit time safeguards expenditure of candidate depot i, i I ;
c i 1 : inventory-related section of unit time unit inventory capacity safeguards expenditure of candidate depot i, i I ;
c i 2 : unit distance unit spare part allocation expense of the candidate depot i, i I ;
c i 3 : unit inventory order outlay of the candidate depot i, i I ;
h i : inventory cost per spare part of the candidate depot i, i I ;
g i : risk shortage loss coefficient of the candidate depot i, i I ;
k i : single inventory review outgoing of the candidate depot i, i I ;
( x i , y i ) : coordinate of base i, i I ;
n: number of depots, a constant;
L: lead time of spare part centralized depots, a constant;
A: lower limit of inventory supply availability for all bases;
α : belief degree for satisfying the criteria of depot service level for each depot;
β : belief degree for guaranteeing the supply available for each base;
γ : stockout risk level for the proposed system;
Z: unit time installation number of spare parts for each equipment;
N i : equipment quantity of base i, i I ;
C: unit time total cost;
C i : unit time total cost of depot i;
C B : unit time overall maintenance expenditure of depots;
C T : unit time overall allocation expense;
C H : unit time overall holding cost of depots;
C R : unit time overall stockout risk loss of depots;
C O : unit time overall order outlay of depots;
ξ i : unit time uncertain demands of the base i, i I ;
f i : shortage value-at-risk quantity of the candidate depot i, i I ;
Φ i x : uncertainty distribution of ξ i , i I ;
Φ i 1 x : inverse uncertainty distribution of ξ i , i I ;
E x : expected value of the uncertain variable x;
Decision variable
S i : designed inventory level of the selected base i with a depot, i I ;
T i : designed review period of the selected base i with a depot, i I ;
X i = 1 , i f  b a s e i h a s  c e n t r a l i z e d  d e p o t 0 , i f  b a s e i c a n t  s t o r e  s p a r e  p a r t s , i I ;
Y i j = 1 , i f  b a s e i i s  s u p p l i e d  b y  s e l e c t e d  b a s e j 0 , i f  b a s e i  i s  n o t  s u p p l i e d  b y  s e l e c t e d  b a s e j , i , j I ;

3.4. Model

According to the problem background, we formulate an uncertain optimization model to realize the total cost minimality without violating the appointed constraints. The mission of the model is to determine the centralized depot locations, inventory parameters, and allocation relationships. The following subsection presents the objective and constraints of the model.

3.4.1. Objective Function

From the cost point of view, this paper summarizes five types of cost, namely, depot maintenance expenditure cost, allocation expense, inventory holding cost, risk shortage loss, and order outlay cost. As for the long-term service of the supply chain, we aspire to minimize the unit time total cost equal to the aggregation of the above five cost compositions, which can be expressed as:
C = C B + C T + C H + C R + C O .
Explicitly, the detailed interpretations of different unit time cost components are listed below.
  • The unit time overall maintenance expenditure of depots
    If base i is appointed as depot location equivalent to X i = 1 , it will undergo two types of cost-consuming maintenance activities. One contains the initial setup and multitudinous renovations which are believed to spend a certain expenditure that is in line with the depot capacity S i . The linear correlation cost coefficient is recorded as c i 1 . Additionally, each support depot on guard should pay some safeguard cost for maintaining its normal operation, which is simplified to be a fixed value b j . Therefore, the unit time overall maintenance expenditure of depots can be evaluated by the sum of the two cost sections of all depots, i.e.,
    C B = i I X i ( b i + c i 1 · S i ) .
  • The unit time overall allocation expense
    On the condition that base i obtains service from depot j ( Y i j = 1 ) , there will be an amount of allocation expense which indicates the cost spent on spare parts transportation between depot j and base i. The expense is associated with the average demands E ξ i , the two sites transportation distance, and the unit distance unit spare part allocation expense of depot j, (denoted by c j 2 ). Hence, the unit time overall allocation expense satisfies the expression:
    C T = j I i I c j 2 · ( x i x j ) 2 + ( y i y j ) 2 · E ξ i · Y i j .
  • The unit time overall holding cost of depots
    The influences of the time-varying inventory level of depot j on cost consumption can be integrated into two aspects under the assumption that X j = 1 .
    On the one hand, when the depot j has a positive inventory level, it needs to pay for managing the inventory and the corresponding expense depending on the stock quantity. In particular, the unit time holding cost of depot j can be computed by multiplying its expected net inventory quantity with a given cost coefficient h j . The expected net inventory is computed by the line integral. As Figure 2 displayed, the inventory level is replenished to S j at the beginning of review periods while the expected inventory level at end of the periods is S j i I E ( ξ i ) · Y i j · T j . Therefore, the expected inventory level is S j 1 2 i I E ( ξ i ) · Y i j · T j . In addition, the spare parts at the lead time will not incur the inventory cost, the corresponding quantifies of which should be deducted from the expected inventory level. Thus, the unit time overall holding cost C H can be assessed by:
    C H = j I X j · h j · S j 1 2 i I E ( ξ i ) · Y i j · T j i I E ( ξ i ) · L · Y i j +
    where the function ( x ) + represents the superior value of x and 0.
  • The unit time overall stockout risk loss of depots
    On the other hand, the negative inventory level attributes to the shortage loss. According to uncertainty theory, the shortage value-at-risk quantity f j as a substitute for the risk index indicates the maximum possible loss number of spare parts in a review period. Through the multiplication of value f j with an allotted coefficient g j , we gain the whole review period shortage loss of depot j. Then, the unit time overall stockout risk loss of depots is demonstrated below:
    C R = j I X j · g j · 1 T j · sup f j | M { i I ξ i · Y i j · T j S j f j } γ +
    where the function ( x ) + represents the superior value of x and 0.
  • The unit time overall order outlay of depots
    Order outlay hails from the rapid inventory counting and the followed order requesting procedures at the end of each period. For depot j, every inventory counting brings about a fixed payment k j , plus an order batch related spend paid for launching periodical orders, where an order batch is the demand totality of the bases set that accords with Y i j = 1 , i I . Accordingly, the unit time order outlay can be written as:
    C O = j I ( i I c j 3 · E ( ξ i ) · Y i j + k j · X j / T j ) .
    In the formula, c j 3 are already known.

3.4.2. Constraints

Subsequently, we will explicit constraints in the proposed model:
  • Supportability constraints
    As a matter of fact, the support system is bound to vouch for the effectiveness of supply. Here, we judge the relative system feature by two supportability criteria. The first is the service level, which can be interpreted as the chance of inventory shortage meeting a given belief degree α :
    M { S j i I ξ i · Y i j · T j 0 } α j I .
    Secondarily, the support demand-side often contracts with the service side that the machines requesting spare parts are in action to some degree, which implies the extent of task completion. In accordance with the definition equation proposed by Sherbrooke 2006, the ‘expected availability’ can be interpreted: if excluding cannibalization and failure coupling, the supply availability A can be indicated as [31]:
    A = 1 E B O N Z Z ,
    where E B O indicates the expected stockout number of spare parts, N is the total number of machines, and Z denotes the number of spare parts per machine.
    Through embedding ‘supply availability’ into our supportability measure architecture, we obtain the constraint as below:
    M i I ξ i · Y i j · T j S j inf { ( 1 A 1 / Z ) · N i · Z · T j | Y i j = 1 } β j I ,
    where β is a predetermined confidence level.
  • Centralized depot number constraint
    In the preliminary scheme planning phase, we assume the number of depots n has been arranged, that is,
    i I X i = n .
  • Supply equilibrium constraint
    The balanced utilization can avoid overusing some of the depots while desolating others. Therefore, the difference in the number of bases served by each depot is to be less than 1. To this end, the supply equilibrium constraint can be written as:
    X j · X k · i I Y i j i I Y i k 1 j , k I , j k .
  • Demand serving constraint
    According to assumption NO. 2, one base should be provided by only one depot, which can be expressed as:
    j I Y i j = 1 i I .

3.4.3. Uncertain Location–Allocation–Inventory Model (ULAIM)

In summary, with the unit time multicomponent cost as the objective function, and the before-mentioned requirements as constraints, we establish the supply relationship selection and inventory parameters optimization model as follows:
min C = i I X i · b + c i ( 1 ) · S i + j I ( i I c j 3 · E ( ξ i ) · Y i j + k j · X j / T j ) + j I i I c j 2 · ( x i x j ) 2 + ( y i y j ) 2 · E ξ i · Y i j + j I X j · h j · S j 1 2 i I E ( ξ i ) · Y i j · T j i I E ( ξ i ) · L · Y i j + + j I X j · g j · 1 T j · sup f j | M i I ξ i · Y i j · T j S j f j γ + s . t . M { S j i I ξ i · Y i j · T j 0 } α j I M i I ξ i · Y i j · T j S j inf { ( 1 A 1 / Z ) · N i · Z · T j | Y i j = 1 } β j I i I X i = n X j · X k · i I Y i j i I Y i k 1 j , k I , j k j I Y i j = 1 i I Y i j { 0 , 1 } , X i { 0 , 1 } i , j I

3.4.4. Deterministic Equivalence of the Model

In order to convenient model solving, the uncertain model is transformed into the equal deterministic form, which is expressed as:
min C = i I X i · b + c i ( 1 ) · S i + j I ( i I c j 3 · E ( ξ i ) · Y i j + k j · X j / T j ) + j I i I c j 2 · ( x i x j ) 2 + ( y i y j ) 2 · E ξ i · Y i j + j I X j · h j · S j 1 2 i I E ( ξ i ) · Y i j · T j i I E ( ξ i ) · L · Y i j + + j I X j · g j · 1 T j · i I Φ i 1 ( 1 γ ) · Y i j · T j S j + s . t . i I Φ i 1 ( α ) · Y i j · T j S j j I i I Φ i 1 ( β ) · Y i j · T j S j + inf { ( 1 A 1 / Z ) · N i · Z · T j | Y i j = 1 } j I i I X i = n X j · X k · i I Y i j i I Y i k 1 j , k I , j k j I Y i j = 1 i I Y i j { 0 , 1 } , X i { 0 , 1 } i , j I

4. Algorithm

So far, diverse modified genetic algorithms that import parameter optimization methods have been widely employed to puzzle out the NP-hard problems in supply chains [32]. In an effort to ease model resolving and optimal in the research, we propose an improved heuristic algorithm that adopted the population renewal idea in the GA and commingled the retrieval thought springing from neighborhood search. The algorithm intends to exhaust the possible search domain more directionally than the traditional GA. During simulation, we have dug out the reality that the supply–demand relationship makes a more remarkable difference to the inventory parameters than the opposite effect. It is the reason why the problem can be disassembled into nested bilevel population sequencing. In the outer layer, our optimization concentrates on determining depot location and supply–demand relationship. In contrast, the inter-layer optimization highlighted in the blue frame aims at programming the inventory parameters.
As shown in Figure 3, the detailed simulation steps are exhibited as follows:
1.
Iterate over all X j , Y i j combinations to generate outer layer populations, and the number of the entire combining samples is treated as the population volume.
2.
For all population samples, compute unit time depot supply quantity and the floor level of S j / T j that satisfies the supportability requirements.
3.
According to the range and required precision, choose the appropriate interval T of the review period.
4.
Generate new populations owing distinctive T j that displays as an arithmetic sequence with T as the difference, like the left section in Figure 4.
5.
For all depots in the generated new populations, initialize the inventory level S j as the minimum supportable inventory level S ( l b , j ) , compute the order outlay, stockout risk loss, and inventory cost of each depot, and then proceed to acquire the addictive total cost of each depot.
6.
Record the current population, these cost components and the additive total costs.
7.
Within the specified number of epochs, mutate S j and solve the correspondent cost components. Compare them with the previous total cost of each depot. If any depot’s total cost is lower, replace the recorded inventory level and the corresponding cost values.
8.
Pick up the population samples with the lowest total cost for each depot and reform the population set in Step 1 (Figure 4 right side shown).
9.
If T does not satisfy the T j ’s precision, make T j T and T j + T to be the lower and upper bound of the review period, and then repeat Step 3~Step 9. If the precision regulation is met, exert the next step.
10.
Sequence population samples and filtrate the best-centralized depot location and provision network architecture with optimized inventory parameters.

5. Numerical Analyses

In the present section, numerical examples of a base-level support system is discussed and analyzed to certify the validity of our proposed model and algorithm.

5.1. Settings

The sizes of the studied system are listed in Table 1 with known coefficients and parameters enumerated in Table 2. Noticeably, the bases’ demands are assumed to be normal uncertain variables. For the optimization of the decision variable, the inventory capacity of each depot should be an integer variable, while, for simplicity, the review period is quoted in two decimal places. Moreover, the proposed algorithm does not utilize crossover process (crossover possibility = 0) and mutation of the inventory level is set to happen each epoch (mutation possibility = 1). This is because out-layer populations are gained by traversing all the possible location-allocation combinations, so that the crossover process is no longer needed. The other configuration parameters are listed in Table 3. The inner populations are generated through taking different values of parameter T at equal intervals period, the intervals T are changed from 1 for roughly search to 0.1 for more accurate optimization and then to 0.01 which conforms to the regulated precision. For each inner population, we alter the inventory level value S for 30 epochs (epoch = 30). If the total cost is lower than the recorded value in any epoch, the corresponding S will replace the recorded value.

5.2. Results

According to the above settings, the procedure is coded in the MATLAB R2019a and tested on a workstation with 256 GB RAM and 2.9 GHz MAD EPYC Processor. It consumes 173.26928 s for solving the case and the results show that the minimum total cost is 416.08 with base 1, base 2, and base 4 chosen as the centralized depot sites. The according second supply–demand relationships are portrayed in Figure 5, while the relative parameters and detailed costs are presented in Table 4 and Table 5.
Table 4 manifests the optimal parameters of the selected depot and their distinctive cost components, while Table 5 summarizes the outcomes of the cost distribution in the entire supply chain.

5.3. Algorithm Comparison

In an effort to prominent the efficacy of the proposed algorithm, we itemize outcomes of the proposed heuristic algorithm and the typical genetic algorithm (populations: 75,600; Epochs: 30; crossover probability: 0.4; mutation probability: 0.95) in Table 6. The conclusion shows that the proposed algorithm is better at approaching the optimal value. Besides, it costs 1103.632415 s to complete computation which is roughly 7 times as long as the proposed algorithm. The results also verify the efficiency of the method.

5.4. Extensive Case Studies

To better explore the general applicability of the model, we also try to change the supportability indexes and the results are exhibited in Table 7. Note that the bases’ positions and demands remain unchanged and the varied parameters contain α , β , γ .
From Table 7, we can summarize that the increase of belief degree α or β may lead to the adjustment of inventory parameters and the increase of total cost and vice versa.The higher belief degree α means the higher service level requirement, and then the corresponding cost will increase. Similarly, the rise of belief degree β equals to the stricter supply availability standard with the sacrifice of cost. In practice, the balanced level of cost and supportability depends on the selection of α and β . The decline of stockout risk level γ does not change the inventory parameter but it represents the descending penalty for shortage risk, which leads to the drop of stockout loss and causes a decrease of the total cost.

5.5. Sensitivity Analysis of Inventory Management

To clarify the impacts of inventory management on diverse cost components, we implement some sensitivity analysis of inventory parameters and then draw the results in Figure 6 and Figure 7. Since the change rules among disparate costs and inventory parameters are similar for each depot, so the following analysis is applicable to every designated depot in the paper.
Firstly, we explore the impact of inventory capacity on disparate expense elements portrayed in Figure 6. According to the supportability constraints and the expression of stockout risk loss, the whole optimal process can be divided into three regions, namely, the unreliable region, the multifactor region, and the none stockout region. The detailed description is as follows:
  • Region I: Unreliable region means that the centralized depot cannot meet the given supportability constraints, which makes the supply system unreliable. In this scene, owing to the substandard inventory level, comparatively frequent shortages and relatively long downtime cause the whole region undesirable.
  • Region II: Multifactor region indicates that the five cost components are in coexistence and exert their synergy effects on the total cost. Precisely, the increment of the total cost is identical to the escalation of inventory plus order expenditure minus the drop of the stockout risk loss. Because of the incompatible offset of the decline in contrast with the upward trend, the total cost exhibits continuing inflation.
  • Region III: None stockout region signifies that shortage can be neglected. Although it is a fact that stockouts are likely to occur in any circumstance, we believe that if the inventory level reaches a certain high level, shortages are too rare to be ignored. Therefore, the transition of the total cost in the third region has a steeper amplification in contrast with region two.
To sum up, along with the accretion of inventory capacity, the cost spent on depot maintenance C B j (j = 1∼3) and parts allocation C T j (j = 1∼3) stays steady, while the expense for inventory holding C H j (j = 1∼3) and ordering C O j (j = 1∼3) keeps growing. The diversity of regions predominantly manifests in shortage cost C O j (j = 1∼3) and the total cost C j (j = 1∼3).
From Figure 7, the red line manifested that the optimal inventory capacity rises linearly with the review period. Simultaneously, with the growth of the review period, the shortage loss (blue line) and the allocation expense (pink line) maintain immutable, whereas the expenditure for inventory (purple line) and maintenance (green line) remain boosting. Nonetheless, since the average per unit time inventory counting cost is proportional to counting frequency 1/ T j , there exists an inversely proportional downward trend to order outlay (yellow line). The aforementioned trends add up to a dip preceding a persistent rise in total cost (the black curve shown). The minimum extreme value is regarded as the final choice of the parameter values marked on the horizontal axis.
Conclusively, the above discussions on the inventory parameter values help to evaluate the behavior of distinctive costs and find the optimal inventory policy that guarantees the system supportability requirements.

6. Conclusions

In this study, we came up with an optimization Location-Allocation-Inventory model for addressing a three-layer base-level spare part supply chain. The goal of the model is to minimize the total cost via appropriate management of centralized depot location, allocation relationships, and inventory policy. On the premise of practical conditions, depots were nominated to be controlled by (T, S) policy and balanced one-to-multiple allocation mode. Regarding supportability judgment, we supposed that the system was restricted by regulated service level and supply supportability. In view of the epistemic uncertainty in the model, uncertainty theory was introduced to quantify demand with stockout loss assessed by risk-at-value. To overcome the obstacle of solution acquirement, we put forward an improved bilevel genetic algorithm whose two stages target to accomplish supply–demand matching and inventory policy updating separately. Finally, numerical analyses are carried out to verify the efficiency of the proposed model and algorithm. The main conclusions can be summarized as follows:
(1) The joint location-inventory-allocation model was solvable with the gained decision variables containing the selected depot locations, allocation relations and inventory parameters.The results also demonstrated the distribution of different cost components.
(2) The suggested algorithm was validated to perform efficiently in optimal solution approaching compared to the traditional genetic algorithm.
(3) Extensive and sensitive analyses explored the general practicality of the model, and illustrated that the inventory parameters T, S had certain impacts on different cost components and the selection of suitable parameters could reduce the total expenditure.
With respect to future explorations, the following research can concentrate on the limitations of this paper. Firstly, other approaches for uncertainty measurement can be wielded for uncertain factor evaluation and compared with the current approach. Secondly, other supportability indicators, like logistics delay time and supply reliability [33], can also be mingled into the spare part provision models. Finally, production and maintenance affairs can also be conducted on the supply chain establishment for support systems.

Author Contributions

Conceptualization, M.W. and R.K.; methodology, P.L. and M.W.; software, P.L.; validation, P.L., M.W. and T.Z.; formal analysis, P.L., M.W. and T.Z.; investigation, P.L.; resources, M.W. and R.K..; data curation, P.L.; writing—original draft preparation, P.L.; writing—review and editing, P.L. and T.Z.; visualization, P.L. and T.Z.; supervision, T.Z., M.W. and R.K.; project administration, M.W. and R.K.; funding acquisition, T.Z., M.W. and R.K. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Stable Supporting project of Science & Technology on Reliability & Environmental Engineering Laboratory (WDZC20220102), Foundation of Science & Technology on Reliability & Environmental Engineering Laboratory (6142004210101), and National Natural Science Foundation of China (No. 62073009, 71671009).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wu, J.; Liu, H.; Zuo, H.; Yang, Y.; Ma, Y.; Kong, L. The demand supply steady-state process-based multi-level spare parts optimization. Sensors 2021, 21, 8324. [Google Scholar] [CrossRef] [PubMed]
  2. VanHoutum, G.J.; Kranenburg, B. Spare Parts Inventory Control under System Availability Constraints; Springer: Berlin/Heidelberg, Germany, 2015; Volume 227. [Google Scholar]
  3. Yao, Z.; Lee, L.H.; Jaruphongsa, W.; Tan, V.; Hui, C.F. Multi-source facility location–allocation and inventory problem. Eur. J. Oper. 2010, 207, 750–762. [Google Scholar] [CrossRef]
  4. Tsao, Y.C.; Mangotra, D.; Lu, J.C.; Dong, M. A continuous approximation approach for the integrated facility-inventory allocation problem. Eur. J. Oper. Res. 2012, 222, 216–228. [Google Scholar] [CrossRef]
  5. Dai, Z.; Aqlan, F.; Zheng, X.; Gao, K. A location-inventory supply chain network model using two heuristic algorithms for perishable products with fuzzy constraints. Comput. Ind. Eng. 2018, 119, 338–352. [Google Scholar] [CrossRef]
  6. Tirkolaee, E.B.; Mahdavi, I.; Esfahani, M.M.S.; Weber, G.W. A robust green Location-Allocation-Inventory problem to design an urban waste management system under uncertainty. Waste Manag. 2020, 102, 340–350. [Google Scholar] [CrossRef]
  7. Basten, R.J.I.; Van Houtum, G.J. System-oriented inventory models for spare parts. Surv. Oper. Res. Manag. Sci. 2014, 19, 34–55. [Google Scholar] [CrossRef] [Green Version]
  8. Gerchak, Y.; Magazine, M.J.; Gamble, A.B. Component commonality with service level requirements. Manag. Sci. 1988, 34, 753–760. [Google Scholar] [CrossRef]
  9. Tompkins, J.A. The Distribution Management Handbook; McGraw-Hill: New York, NY, USA, 1994. [Google Scholar]
  10. Asian, S.; Nie, X. Coordination in supply chains with uncertain demand and disruption risks: Existence, analysis, and insights. IEEE Trans. Syst. Man Cybern. Syst. 2014, 44, 1139–1154. [Google Scholar] [CrossRef]
  11. Islam, M.; Azeem, A.; Jabir, M.; Paul, A.; Paul, S.K. An inventory model for a three-stage supply chain with random capacities considering disruptions and supplier reliability. Ann. Oper. Res. 2022, 315, 1703–1728. [Google Scholar] [CrossRef]
  12. Errico, F.; Desaulniers, G.; Gendreau, M.; Rei, W. The vehicle routing problem with hard time windows and stochastic service times. EURO J. Transp. Logist. 2018, 7, 223–251. [Google Scholar] [CrossRef]
  13. Zhang, B.; Li, H.; Li, S.; Peng, J. Sustainable multi-depot emergency facilities location-routing problem with uncertain information. Appl. Math. Comput. 2018, 333, 506–520. [Google Scholar] [CrossRef]
  14. Loeve, M. Probability Theory; Courier Dover Publications: New York, NY, USA, 2017. [Google Scholar]
  15. Gholamian, N.; Mahdavi, I.; Mahdavi-Amiri, N.; Tavakkoli-Moghaddam, R. Hybridization of an interactive fuzzy methodology with a lexicographic min-max approach for optimizing a multi-period multi-product multi-echelon sustainable closed-loop supply chain network. Comput. Ind. Eng. 2021, 158, 107282. [Google Scholar] [CrossRef]
  16. Gai, L.; Liu, H.; Wang, Y.; Xing, Y. Green supplier selection and order allocation using linguistic z-numbers multimoora method and bi-objective non-linear programming. Fuzzy Optim. Decis. Mak. 2022. [Google Scholar] [CrossRef]
  17. De, A.; Singh, S.P. Analysis of fuzzy applications in the agri-supply chain: A literature review. J. Clean. Prod. 2021, 283, 124577. [Google Scholar] [CrossRef]
  18. Díaz-Curbelo, A.; Andrade, R.A.E.; Municio, Á.M.G. The role of fuzzy logic to dealing with epistemic uncertainty in supply chain risk assessment: Review standpoints. Int. J. Fuzzy Syst. 2020, 22, 2769–2791. [Google Scholar] [CrossRef]
  19. Liu, B. Uncertainty Theory, 2nd ed.; Springer: Berlin/Heidelberg, Germany, 2007. [Google Scholar]
  20. Sheng, L.; Zhu, Y.; Wang, K. Uncertain dynamical system-based decision making with application to production-inventory problems. Appl. Math. 2018, 56, 275–288. [Google Scholar] [CrossRef]
  21. Asim, Z.; Jalil, S.A.; Javaid, S. An uncertain model for integrated production-transportation closed-loop supply chain network with cost reliability. Sustain. Prod. Consum. 2019, 17, 298–310. [Google Scholar] [CrossRef]
  22. Shen, J. An uncertain sustainable supply chain network. Appl. Math. Comput. 2020, 378, 125213. [Google Scholar] [CrossRef]
  23. Hosseini, A.; Pishvaee, M.S. Capacity reliability under uncertainty in transportation networks: An optimization framework and stability assessment methodology. Fuzzy Optim. Decis. Mak. 2022, 21, 479–512. [Google Scholar] [CrossRef]
  24. Wu, W.; Zhou, W.; Cheng, Z. A hybrid metaheuristic algorithm for location inventory routing problem with time windows and fuel consumption. Expert Syst. Appl. 2021, 166, 114034. [Google Scholar] [CrossRef]
  25. Tirkolaee, E.B.; Mahmoodkhani, J.; Bourani, M.R.; Moghaddam, R.T. A self-learning particle swarm optimization for robust multi-echelon capacitated location–allocation–inventory problem. J. Adv. Manuf. 2019, 18, 677–694. [Google Scholar] [CrossRef]
  26. Mousavi, S.M.; Pardalos, P.M.; Niaki, S.T.A.; Fügenschuhd, A.; Fathi, M. Solving a continuous periodic review inventory-location allocation problem in vendor-buyer supply chain under uncertainty. Comput. Ind. Eng. 2019, 128, 541–552. [Google Scholar] [CrossRef]
  27. Sahraei, Z.; Samouei, P. Bi-level scenario-based Location-Allocation-Inventory models for post-crisis conditions and solving with electromagnetic and genetic algorithms. Int. J. Eng. 2022, 35, 1803–1819. [Google Scholar] [CrossRef]
  28. Auweraer, S.V.; Boute, R.N.; Syntetos, A.A. Forecasting spare part demand with installed base information: A review. Int. J. Forecast. 2019, 35, 181–196. [Google Scholar] [CrossRef] [Green Version]
  29. Aisyati, A.; Jauhari, W.A.; Rosyidi, C.N. Periodic review model for determining inventory policy for aircraft consumable spare parts. Int. J. Bus. Res. Manag. 2014, 5, 41–45. [Google Scholar]
  30. Cuc, V.T.K.; Thi, N.T.; LeThuy, N.T. A simulation modelling approach for selection of inventory policy in a supply chain. In Proceedings of the 2019 International Conference on System Science and Engineering (ICSSE), Dong Hoi, Vietnam, 20–21 July 2019; pp. 349–354. [Google Scholar]
  31. Sherbrooke, C.C. Optimal Inventory Modeling of Systems: Multi-Echelon Techniques; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2006. [Google Scholar]
  32. Maiti, M.K.; Maiti, M. Two storage inventory model in a mixed environment. Fuzzy Optim. Decis. Mak. 2007, 6, 391–426. [Google Scholar] [CrossRef]
  33. Jones, J.V. Supportability Engineering Handbook; McGraw-Hill Education: New York, NY, USA, 2007. [Google Scholar]
Figure 1. The framework of a three-echelon supply chain with product and information flow.
Figure 1. The framework of a three-echelon supply chain with product and information flow.
Axioms 12 00046 g001
Figure 2. The time-variant characteristic sketch maps: (a) the distribution of demand from base ξ and (b) inventory level S.
Figure 2. The time-variant characteristic sketch maps: (a) the distribution of demand from base ξ and (b) inventory level S.
Axioms 12 00046 g002
Figure 3. Key flows of the proposed heuristic algorithm.
Figure 3. Key flows of the proposed heuristic algorithm.
Axioms 12 00046 g003
Figure 4. Internal structure of the bilevel population samples.
Figure 4. Internal structure of the bilevel population samples.
Axioms 12 00046 g004
Figure 5. Allocation relationships gained from simulation.
Figure 5. Allocation relationships gained from simulation.
Axioms 12 00046 g005
Figure 6. Sensitivity analysis of inventory capacity.
Figure 6. Sensitivity analysis of inventory capacity.
Axioms 12 00046 g006
Figure 7. Sensitivity analysis of inventory period.
Figure 7. Sensitivity analysis of inventory period.
Axioms 12 00046 g007
Table 1. The size of the problem.
Table 1. The size of the problem.
SupplierDepotsBase Zone
1310
Table 2. Specified parameters for the supply chain.
Table 2. Specified parameters for the supply chain.
i( x i , y i ) ξ i h i g i k i N i
e i s i
1(44, 98)83160.230.187435
2(16, 21)78130.240.171406
3(40, 53)85100.260.138449
4(96, 14)84150.220.199448
5(74, 67)77130.240.189457
6(24, 42)84120.250.185423
7(5, 2)78170.270.152414
8(54, 83)83120.250.172468
9(64, 48)70130.240.144463
10(98, 99)85160.260.184477
ξ i follows a normal uncertainty distribution, namely, ξ i N ( e i , s i ). Other parameters for the supply chain: n = 3, L = 0.01, α = 0.9, β = 0.9, γ = 0.01, A = 0.85, b = 5. Cost set: ∀iI, c i ( 1 ) = 0.01, c i ( 2 ) = 0.002, c i ( 3 ) = 0.5. Limited by the ability about supplying and inventory checking, the review time T i is in the range from 0.5 to 5.
Table 3. The configuration of the proposed algorithm.
Table 3. The configuration of the proposed algorithm.
Outer-Layer Populations T Inter-Layer PopulationsEpoches of Mutating S j
75,600(1, 0.1, 0.01)(5, 20, 20)30
Table 4. The parameters and various cost values of centralized depots.
Table 4. The parameters and various cost values of centralized depots.
ParameterDepot 1Depot 2Depot 3
Site selection1(44, 98)2(16, 21)4(96, 14)
Supply relationship1(1, 3, 8, 10)2(2, 6, 7)4(4, 5, 9)
Review period T j 0.860.951.06
Stocked capacity S j 346277298
Maintenance expenditure8.467.777.98
Transportation expense9.933.607.69
Inventory cost45.5838.5438.12
Stockout risk loss13.599.6510.92
Order outlay83.6066.1164.61
Total cost161.15129.51129.31
Table 5. The different cost components.
Table 5. The different cost components.
C B C T C H C R C O Total Cost C
24.2021.21122.0134.34214.31416.08
Table 6. Configurations and outcomes of the two algorithms.
Table 6. Configurations and outcomes of the two algorithms.
OutcomesThe Proposed AlgorithmThe Typical Genetic Algorithm
Site selection1, 2, 41, 2, 4
Supply relationship1(1, 3, 8, 10), 2(2, 6, 7), 4(4, 5, 9)1(1, 3, 8, 10), 2(2,6,7), 4(4,5,9)
Review time of depot1(0.86), 2(0.95), 4(1.06)1(0.92), 2(0.81), 4(0.98)
Stocked capacity1(346), 2(582), 4(562)1(371), 2(582), 4(562)
Maintenance expenditure1(8.46), 2(7.77), 4(7.98)1(8.71), 2(7.43), 4(8.13)
Transportation expense1(9.93), 2(3.60), 4(7.69)1(9.93), 2(3.60), 4(7.69)
Inventory cost1(45.58), 2(38.54), 4(38.12)1(49.01), 2(34.42), 4(43.45)
Stockout risk loss1(13.59), 2(9.65), 4(10.92)1(13.42), 2(8.21), 4(3.31)
Order outlay1(83.60), 2(66.11), 4(64.61)1(80.34), 2(73.38), 4(68.00)
Total cost416.08419.01
Table 7. Extensive case studies.
Table 7. Extensive case studies.
α β γ Inventory ParameterTotal CostRunning Time
0.90.90.011(0.86, 346); 2(0.95, 277); 4(1.06, 298)416.083173.269280 s
0.850.90.011(0.86, 345); 2(0.95, 276); 4(1.06, 297)415.952171.329017 s
0.90.950.011(0.85, 364); 2(0.95, 293); 4(1.05, 315)418.338172.475892 s
0.90.90.051(0.86, 346); 2(0.95, 277); 4(1.06, 298)393.071171.735646 s
The location-allocation relations keep unchanged in the programs.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Li, P.; Wen, M.; Zu, T.; Kang, R. A Joint Location–Allocation–Inventory Spare Part Optimization Model for Base-Level Support System with Uncertain Demands. Axioms 2023, 12, 46. https://doi.org/10.3390/axioms12010046

AMA Style

Li P, Wen M, Zu T, Kang R. A Joint Location–Allocation–Inventory Spare Part Optimization Model for Base-Level Support System with Uncertain Demands. Axioms. 2023; 12(1):46. https://doi.org/10.3390/axioms12010046

Chicago/Turabian Style

Li, Peixuan, Meilin Wen, Tianpei Zu, and Rui Kang. 2023. "A Joint Location–Allocation–Inventory Spare Part Optimization Model for Base-Level Support System with Uncertain Demands" Axioms 12, no. 1: 46. https://doi.org/10.3390/axioms12010046

APA Style

Li, P., Wen, M., Zu, T., & Kang, R. (2023). A Joint Location–Allocation–Inventory Spare Part Optimization Model for Base-Level Support System with Uncertain Demands. Axioms, 12(1), 46. https://doi.org/10.3390/axioms12010046

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop