Next Article in Journal
Enhancing Multi-Modal Perception and Interaction: An Augmented Reality Visualization System for Complex Decision Making
Next Article in Special Issue
AI Lifecycle Zero-Touch Orchestration within the Edge-to-Cloud Continuum for Industry 5.0
Previous Article in Journal
The Knowledge Spillover Effect of Multi-Scale Urban Innovation Networks on Industrial Development: Evidence from the Automobile Manufacturing Industry in China
Previous Article in Special Issue
Bibliometric Analysis of a Product–Service System’s Rebound Effect: Identification of a Potential Mitigation Strategy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling and Solving the Joint Replenishment Problem with Cross-Selling Effects Considering One Shared Minor Item

School of Economics and Management, Beihang University, Beijing 100191, China
*
Author to whom correspondence should be addressed.
Systems 2024, 12(1), 6; https://doi.org/10.3390/systems12010006
Submission received: 16 November 2023 / Revised: 16 December 2023 / Accepted: 18 December 2023 / Published: 22 December 2023
(This article belongs to the Special Issue Manufacturing and Service Systems for Industry 4.0/5.0)

Abstract

:
In this paper, we provide a model to handle multiple replenishment cycles and the cross-selling of multiple major items with one minor item, while allowing partial late delivery. The optimization analytic expression of the model is finally obtained by utilizing the convexity of cost function for F and using the first-order conditions in optimization theory. Numerical examples and sensitivity analysis demonstrate the effectiveness of the model and algorithm, which offers a competent solution for practical applications.

1. Introduction

Demand independence is one of the fundamental assumptions in classical economic order quantity (EOQ) models and joint replenishment problem (JRP) models [1,2,3]. However, in practice, managing products with independent demand often falls short of achieving the desired level of customer satisfaction [4]. Therefore, studies have long proposed that dependent demand should be considered in inventory management and production planning decisions. Stock-dependent demand [5], and demand for substitute products [6,7] are common non-independent demands in research. In recent years, there has been a significant increase in research on non-independent demands, leading to the development of a broad category of inventory and production planning models and decision-making methods for such demands [8,9].
Cross-selling implies that the demand for or sale of a major item will lead to an additional demand for its minor items, making it another significant factor in causing non-independent demand. Agrawal et al. [10] introduced the concept of “itemset” in the context of cross-selling, which generalized the description of cross-selling effects and expanded its application domain. Some scholars have addressed the issue of product selection in inventory by using association rules to reflect the impact of cross-selling effects [11,12,13]. Recently, the attention of scholars has focused on cross-selling in omnichannel strategies, cross-selling in recommender systems, and the usage of data mining methods in cross-selling. Liu et al. [14] have focused on omnichannel retailing with different order fulfillment and return options and found that the cross-selling benefit and the offline search cost have a significant impact on the retailer’s optimal omnichannel strategy. Yang and Ji [15] discussed the impact of cross-selling on managing consumer returns in omnichannel operations. Ghoshal et al. [16] have studied recommendations and cross-selling pricing strategies when personalizing firms’ cross-sell. Mokhtari [17] has addressed an economic order quantity (EOQ) model to determine the joint ordering policy for two products under completion and substitution conditions. Scholars have also used data mining methods to study multiproduct newsvendors with cross-selling and narrow-bracketing behavior [18]. Although these studies have given us many new insights, it is still order strategies that are most closely related to cross-selling.
In the field of inventory management, there have been some advancements in order policies that have considered cross-selling effects [19]. Zhang [20] provided a partial back-ordering EQO model to handle correlated demands caused by cross-selling in the context of the joint replenishment of multiple minor items. This led to the proposal of a series of JRP models that consider correlated demand caused by cross-selling, and a set of heuristics and exact solution algorithms have been generated around these models [21,22,23]. However, existing studies have only examined the management problem for one major item and one minor item [24] or the joint replenishment problem (JRP) for one major item and multiple minor item [21]. The many-to-one cross-selling JRP problem still needs to be studied further.
To address the gaps in existing research, this paper sets a more universal application scenario where multiple major products correspond to one minor product for the joint replenishment problem with cross-selling effects. Specifically, in the analysis of correlated demand caused by cross-selling, we divide the relationship of products into major items and minor items. One shared minor item can be applied to multiple major items. For example, a laptop computer may require multiple complementary products such as mice, keyboards, memory sticks, CPUs, and so on, while one of the minor product (e.g., mice, keyboards, memory sticks) can also be used in multiple electronic products such as laptops, desktop computers, and tablets. By proposing the JRP under this scenario and designing corresponding algorithms to solve the optimization problem, we believe that our research can solve the replenishment problem with the partial delayed delivery and related demand determinism caused by many-to-one cross-selling, which extends current relevant research.
The rest of this paper is organized as follows. Section 2 describes the model construction, that is, a joint replenishment problem with cross-selling effects considering one shared minor item. Section 3 addresses the analytical optimization algorithm for our proposed model, including objective function transformation, analytical process, and optimization strategy. Section 4 provides a numerical example for illustration. A detailed sensitivity analysis is performed in Section 5 around the numerical example to observe our model under various parameters. The paper concludes in Section 6.

2. Modeling Analysis

2.1. Model Assumptions

Previous research on deterministic EOQ often limited one major item to be associated with one minor item [20,25]. In our work, this limitation is removed, making the model applicable to a wider range of problems. Assume that the demand for one shared minor item is related to the demand of N major items whose delivery can be partially delayed. Due to cross-selling, the sale of the major item may lead to additional demand for the shared minor item, which can be sold either separately or jointly with the major item. The inventory of the associated minor item will be reduced when the major item is shipped out.
For the other parts of the assumptions, we follow classic practices from the past around deterministic EOQ [24,25,26]. In other words, even though new work exists that improves the underlying model in various directions to bring it closer to reality, considering the complexity of the solving, we only extend the original classical models, keeping their assumptions intact. This paper mainly assumes that the major item has different cycles and the cycle of the minor item has no multiple relationship with the cycle of major item. At the same time, the following assumptions related to the model are given:
  • Replenishment is instantaneous;
  • Delivery time is zero;
  • The major item can be partially delayed, while the minor item should not be kept in stock;
  • The unit loss of the major item will reduce the demand for the minor item at a constant rate;
  • The demand and cost of each major item are independent of each other.
In real life, the replenishment cycle of various products cannot be exactly the same, especially for large manufacturers. For upstream manufacturers, similar to the retail sector, they need to adjust their production schedules for various products based on demand and source the required raw materials from downstream manufacturers. The same replenishment cycle will lead to excessive transportation pressure and cause problems such as capital turnover, so for different major products the replenishment cycle of each product is not the same. For example, laptop manufacturers produce computers with different functions, including game books, office books, and so on. If the replenishment cycle is the same, it is not only unrealistic in terms of customer demand but it is also unable to have a large inventory to store too many products for replenishment at one time.

2.2. Meaning of Variables in the Model

The meanings of the variables (exogenous variables and decision variables) to be used in this section are shown in Table 1 and Table 2 below.

2.3. Model Description

According to the basic deterministic EOQ model [25], which allows for stock-outs, the ordering cost per unit of time for major item i can be shown by Equation (1).
Γ i T i , F i = A i T i + C h i D i T i F i 2 2 + β i C b i D i ( 1 F i ) 2 2 + C o i D i 1 β i 1 F i
There has been some work introducing cross-selling into the inventory management problem [19], and we use the same approach to construct our model. Similarly, the cost per unit of time function for the minor item can be given by Equation (2).
Γ T 0 , F i = A T 0 + i = 1 N λ i C h 2 D i T i F i 2 + λ i C o D i 1 β i 1 F i + C h D λ i D i T i 2
The final objective function is to minimize the ordering cost of the items, i.e., the sum of the unit time cost of all major items plus the unit time cost of the minor item, and then the total ordering cost of the manufacturer per unit time can be obtained, as shown in Equation (3).
min Γ T , T 0 , F = i = 1 N Γ i T i , F i + Γ T 0 , F i

3. Model Solving

3.1. Objective Function Transformation

The analytical solution method is chosen in this section to derive the results of the model optimization problem established in Section 2. Compared with the computer simulation solution, the analytical solution does not require complex programming algorithms or simulation calculations, and the results obtained by mathematical derivation are more theoretical and more accurate than the numerical simulation results.
First of all, the objective function is written for the above optimization problem as shown in Equation (4), i.e., the expansion form of Equation (3).
Γ T , T 0 , F = i = 1 N A i T i + i N C h i + λ i C h 2 D i T i F i 2 + i N β i C b i D i ( 1 F i ) 2 2 + i N ( C o i + λ i C o ) D i ( 1 β i ) 1 F i + i N C h D λ i D i 2 T i + A T 0  
As assumed in the model in Section 2.1, the demand and cost of each major item are independent of each other, so the total optimization cost and the cost of each major item are equivalent to the cost generated by the major item and the optimization of the addition term can be converted into the optimization solution of each term. By collapsing the total cost with respect to the cost of each major item and the cost due to the major item as Γ ~ i T i , F i , Equation (4) can be written as:
Γ ~ T , T 0 , F = i N Γ ~ i T i , F i + A T 0
Γ ~ i T i , F i = A i T i + C h i + λ i C h 2 D i T i F i 2 + β i C b i D i 2 ( 1 F i ) 2 + ( C o i + λ i C o ) D i ( 1 β i ) ( 1 F i ) + C h D λ i D i 2 T i
For ease of reading, Γ ~ i T i , F i is abbreviated as Γ ~ i in the derivation and calculation section later in this section. In addition, in order to facilitate the solution and enhance the readable lines of the solution step, the non-decision variables in (6) are organized into the form of simplified coefficients, which are, respectively, G 1 i , G 2 i , G 3 i , G 4 i , as shown in Equations (7)–(10). According to the definition of each variable, it can be seen that they are positive real numbers greater than zero.
G 1 i = C h i + λ i C h 2 D i > 0
G 2 i = β i C b i D i 2 > 0
G 3 i = C 0 i + λ i C o D i 1 β i > 0
G 4 i = C h D λ i D i 2 > 0
Substituting G 1 i , G 2 i , G 3 i , G 4 i into Equation (6), a concise representation of Γ ~ i can be obtained, as shown in Equation (11).
Γ ~ i = A i T i + G 1 i T i F i 2 + G 2 i 1 F i 2 + G 3 i 1 F i + G 4 i T i
At this time, the total cost, Γ ~ T , T 0 , F , is a function only related to T, F, T 0 , and this function can roughly separate T and F, so that the two have only one multiplying term, which is a very good optimization function and can bring a significant amount of convenience in the derivation.

3.2. Model Solving Procedure

Having simplified this equation, we can find the analytic solution of this function by taking the derivative of Γ ~ i . According to the following derivation result, Equations (12) and (13), it can be seen that this is a binary quadratic equation.
d Γ i d T i = A i T i 2 + G 1 i F i 2 + G 4 i
d Γ i d F i = 2 G 1 i T i F i 2 G 2 i + 2 G 2 i F i G 3 i
According to the first-order condition, the following optimization solution steps can be made. First, set Equation (12) to 0, i.e., A i T i 2 + G 1 i F i 2 + G 4 i = 0 , and the expression of T i can be calculated, as shown in Equation (14).
T i = A i G 1 i F i 2 + G 4 i
Similarly, if Equation (14) is set to 0, i.e., 2 G 1 i T i F i 2 G 2 i + 2 G 2 i F i G 3 i = 0 , the expression of optimal F i regarding T i can be calculated, as shown in Equation (15).
F i = 2 G 2 i + G 3 i 2 G 1 i T i + 2 G 2 i
Since T i is included in the expression of F i 2 , the expression of F i (the square of Equation (15)) can be substituted into Equation (14) to facilitate the solution of the expression of F i without T i . The square of Equation (15) is shown in Equation (16).
4 G 2 i 2 + 4 G 2 i G 3 i + G 3 i 2 4 G 1 i T i 2 + 8 G 1 i G 2 i T i + 4 G 2 i 2 = F i 2
Substituting the result of Equation (16) into Equation (14), an equation containing only one decision variable, T i , can be obtained, as shown in Equation (17).
G 1 i ( 2 G 2 i + G 3 i ) 2 ( 2 G 1 i T i + 2 G 2 i ) 2 = A i G 4 i T i 2 T i 2
It is easy to see that Equation (17) is a function only related to T, so it is inevitable that an analytical solution of T is obtained. The form of Equation (17) is rewritten into a polynomial form, as shown in Equation (18).
4 G 1 i G 2 i 2 T i 2 + 4 G 1 i G 2 i G 3 i T i 2 + G 3 i 2 T i 2 = 4 A i G 1 i T i 2 + 8 A i G 1 i G 2 i T i + 4 A i G 2 i 2 4 G 1 i G 4 i T i 4 8 G 1 i G 2 i G 4 i T i 3 4 G 2 i 2 G 4 i T i 2
After combining similar terms, the simplified form of Equation (18) can be obtained as Equation (19).
4 G 1 i G 4 i T i 4 + 8 G 1 i G 2 i G 4 i T i 3 + [ G 1 i ( 2 G 2 i + G 3 i ) 2 4 A i G 1 i + 4 G 2 i 2 G 4 i ] T i 2 8 A i G 1 i G 2 i T i + 4 A i G 2 i 2 = 0  
According as Equation (19), this is a quartic equation of T, and the optimal solution of T is the solution of this quartic equation.

3.3. The Optimization Strategy of the Model

In order to facilitate the representation and calculation in the process of model optimization and enhance the simplicity of the operation process, Equation (6) was re-simplified and written in the following form as Equation (20),
Γ i ~ T i , F i = A i T i + T i u F i + v F i ,
where u F i and v F i in Equation (20) are two functions of F i , as shown in Equations (21) and (22).
u F i = G 1 i F i 2 + G 4 i
v F i = G 2 i 1 F i 2 + G 3 i 1 F i
Therefore, Equation (14) can be written as Equation (23) in the following form:
T i F i = A i u F i
Substituting Equation (23) into Equation (20), the form of Equation (24) can be obtained as follows:
Γ i ¯ F i = A i A i u F i + A i u F i u F i + v F i = A i [ u ( F i ) ] 1 2 A i 1 2 + A i 1 2 [ u ( F i ) ] 1 2 u F i + v F i = A i 1 2 [ u ( F i ) ] 1 2 + A i 1 2 [ u F i ] 1 2 + v F i = 2 A i u ( F i ) + v ( F i )
It is easy to show that the function Γ i ¯ ( F i ) is continuous on an interval of F [ 0 ,   1 ] , so the optimal solution can be obtained by taking the derivative.
Firstly, the first and second derivatives of functions u ( F i ) and v ( F i ) are given, as shown in Equations (25)–(28).
u F i = 2 G 1 i F i
u F i = 2 G 1 i
v F i = 2 G 2 i F i 2 G 2 i G 3 i
v F i = 2 G 2 i
Then, the first and second derivatives of Γ i ¯ ( F i ) can be calculated, as shown in Equations (29) and (30).
Γ i ¯ F i = A i u F i u F i + v F i = A i 2 G 1 i F i u F i + 2 G 2 i F i 2 G 2 i G 3 i
2 Γ i ¯ 2 F i = A i [ 2 u ( F i ) u ( F i ) ( u ( F i ) ) 2 ] 2 [ u ( F i ) ] 3 2 + v ( F i ) = A i 8 G 1 i 2 F i 4 G 1 i 2 F i 2 2 [ u ( F i ) ] 3 2 + 2 G 2 i
According to Equations (7)–(9) and (21) and the definition of satisfying rate F i , the following conditions can be easily obtained:
G 1 i > 0 ,   0 F i 1 ,   2 G 2 i > 0 ,   u ( F i ) > 0  
Therefore, it is easy to draw the inequality, as in Equation (31).
8 G 1 i 2 F i 4 G 1 i 2 F i 2 > 0
As the inequality relation of Equation (31), it is easy to deduce 2 Γ i ¯ 2 F i > 0 . According to the second-order condition of convex function, it is obvious that Γ ¯ i ( F i ) is convex on the domain of F i .
Since F i [ 0 ,   1 ] , three cases are discussed below.
Case 1: F i = 0
As shown in Figure 1, when F i = 0 , Γ i ¯ F i = 2 G 2 i G 3 i 0 can be calculated. It follows from the first-order property that Γ i ¯ ( F i ) is monotonically decreasing near F i = 0 , which means that the optimal solution can never be taken at F = 0.
Therefore, the following conclusion can be drawn: the objective function cannot obtain the optimal value when F i = 0 .
Case 2: F i = 1
If F i = 1 , Equation (32) can be calculated.
Γ i ¯ F i   F i = 1 = A i 2 G 1 i G 1 i + G 4 i + 2 G 2 i 2 G 2 i G 3 i
Since Equation (32) depends on parameter values, its positive and negative values are difficult to determine, so it needs to be discussed by case.
As shown in Figure 2, similar to the case of a major item with one minor item, the optimal solution F i 0 , 1 when Equation (32) > 0, i.e., the optimal solution, T i , is the positive real root of the quartic equation of one variable, Equation (19), and then the optimal solution, F i , is obtained by Equation (15).
In this case, we can take the derivative of Γ i ¯ , Equation (33).
Γ i ¯ F i F 1 = A i G 1 i + G 4 i G 3 i > 0
By simplifying Equations (5) and (33), we can obtain:
A i G 1 i + G 4 i G 3 i > 0
Substituting the original value of G 1 i , G 3 i , G 4 i into Equation (34), a critical value of the out-of-stock rate can be obtained as follows:
A i G 1 i + G 4 i ( C o i + λ i C o ) D i + ( C o i + λ i C o ) D i β i > 0 ( C o i + λ i C o ) D i β i > ( C o i + λ i C o ) D i A i G 1 i + G 4 i β i > 1 A i ( C o i + λ i C o ) D i G 1 i + G 4 i β i > 1 1 C o i + λ i C o D i A i G 1 i + G 4 i
Since the optimal solution to T is T = A i G 1 i F i 2 + G 4 i , Equation (34) can be rewritten as Equation (35), as follows:
β i > β C i = 1 1 C o i + λ i C o T ,
where a β C i is determined by Equation (35); it is the lower threshold of the out-of-stock rate.
Therefore, the process of the optimal strategy is as follows:
  • For each major item, a β C i is determined to judge the relationship between the shortage rate and the critical value.
  • When the shortage rate is far less than the critical value, the optimal solution of F i is 0 or 1 and the corresponding value of T i is calculated according to Equation (14), then the value of the objective function is calculated. A group of F i and T i that makes the objective function smaller is selected as the optimal solution.
  • When the shortage rate is greater than the critical value, the optimal T i is a positive real root of the quadratic equation in Equation (19) and the corresponding F i is calculated according to Equation (15).
  • The strategy of one major item and one minor item is used to determine the T i , F i of each major item.
  • T 0 is determined by assuming that T 0 is the largest period of all T i .
  • All the cost items are totaled to find the final total cost.

4. Numerical Computations

A numerical example is used to validate the algorithm in this section, which includes different cases where the shortage rate is greater than or less than the critical value. To ensure that this numerical example is appropriate for the context of our study, i.e., the demand for one minor item is related to the demand for N major items, we introduce a scenario that is close to reality: Consider the minor item as one type of mouse of a certain brand, and the major items as three different models of laptop computers. Product 1 represents the laptop with a long life cycle and a general market demand. Product 2 represents the emerging laptop with a short life cycle and high market interest. Product 3 represents the laptop with a long life cycle and relatively low market demand for a specific group of people. According to practical experience, one kind of mouse is generally applicable to many models of laptops, and the demand for the mouse tends to increase after the sale of laptops. Different models of laptops are applicable to different groups of people, and there is usually a 10-fold or more difference in the demands. Specific parameters are shown in Table 3, Table 4 and Table 5:
Major items:
Table 3. Major item Dataset.
Table 3. Major item Dataset.
Product Requirement   D (Pieces/Year) Fixed   Ordering   Cost   A ($/Time) Cost   of   Sales   Losses   C o ($/Piece) Holding   Cost   C h ($/Piece/Year) Delivery   Delaying   Cost   C b ($/per/Year) Allowable   Delivery   Delaying   Rate   β
i = 1 40020052100.7
i = 2 800600153080.2
i = 3 104002150.96
Minor item:
Table 4. Minor item Dataset.
Table 4. Minor item Dataset.
Product Requirement   D
(Pieces/Year)
Fixed   Ordering   Cost   A ($/Time) Cost   of   Sales   Losses   C o ($/Piece) Holding   Cost   C h ($/Piece/Year)
Mouse40020051
The proportion of related product sales loss caused by the shortage of a major item:
Table 5. The proportion of related product sales loss caused by the shortage of major item i .
Table 5. The proportion of related product sales loss caused by the shortage of major item i .
i = 1 i = 2 i = 3
λ i 0.10.50.3
The critical value of β of each major item can be calculated as Table 6:
According to the algorithm, the optimal T and F of each major item can be calculated as Table 7:
The higher the cycle T 0 of the minor item is, the lower the total cost is. Set T 0 = max T i = 2.4759 . Substituting all the T i , F j , T 0 obtained into the original objective function, the total cost of the optimal solution of this example can be obtained as m i n Γ i T i , F j , T 0 = 1370.83 .

5. Model Sensitivity Analysis

5.1. Sensitivity Analysis of Model Optimal Cost to Ordering Cost

The numerical example in the numerical calculation experiment in the previous section will continued to be used in this section. The ordering cost A of the minor item and the ordering cost A i of the major item i are respectively generated, and 10 different examples are calculated according to the range [100, 10,000] to obtain the optimal solution by keeping other parameters in the numerical example unchanged, as shown in Figure 3 and Figure 4.
As can be seen from Figure 3 and Figure 4, the optimal result calculated by the model is greatly influenced by it (reflected in the high rate of change in the figure) when the ordering cost of minor item A is within the interval of [100, 2000]. At the same time, the optimal result calculated by the model is influenced very little by ordering cost, A i , of the major item i , which is close to linear influence.

5.2. Sensitivity Analysis of Model Optimal Cost to Demand Rate

By keeping other parameters in the numerical example unchanged, the demand rate, D , of the minor item and the ordering cost, D i , of major item i are respectively calculated for 10 different examples in accordance with the ranges [400, 2200] and [200, 2000] to obtain the optimal solution, as shown in Figure 5 and Figure 6.
As can be seen from the two figures above, when the demand rate, D , of the minor item is within the range of [400, 800], the optimal result obtained by the model calculation is greatly affected by it. At the same time, the optimal result calculated by the model is influenced very little by the ordering cost, D i , of major item i , which is close to linear influence.

5.3. Sensitivity Analysis of Model Optimal Cost to Order Backlog Cost

By keeping other parameters in the numerical example unchanged, 10 different examples of order backlog cost, C b i , of major item i are generated in accordance with the range [1, 30] and the optimal solution is obtained, as shown in Figure 7.
As can be seen from the figure above, the optimal cost calculated by the model proposed in this section will not be affected by the order backlog cost, C b i , of major item i .

5.4. Sensitivity Analysis of Model Optimal Cost to Inventory Carrying Cost

By keeping other parameters in the numerical examples unchanged, 10 different examples are respectively generated for the inventory holding cost, C h , of the minor item and the inventory holding cost, C h i , of major item i in accordance with the range [1, 10] to obtain the optimal solution, as shown in Figure 8 and Figure 9.
As can be seen from the two curves above, the optimal result obtained by this model is influenced by the inventory holding cost, C h , and the inventory holding cost, C h i , of major item i , which is very small and close to linear.

5.5. Sensitivity Analysis of Model Optimal Cost to Opportunity Cost

By keeping other parameters in the numerical examples unchanged, 10 different examples are respectively generated for the inventory holding cost, C h , of related products and the inventory holding cost, C h i , of major item i in accordance with the range [1, 10] to obtain the optimal solution, as shown in Figure 10 and Figure 11.
As can be seen from the two figures above, the optimal result calculated by this model is influenced very little by the opportunity cost of the minor item, C o , and the inventory carrying cost, C o i , of major item i and only fluctuates within a very small range.

5.6. Sensitivity Analysis of Model Optimal Cost to Order Backlog Rate

By keeping other parameters in the numerical examples unchanged, 10 different examples are respectively generated for the order backlog rate β i of major item i in accordance with the range [0.1, 0.95] to obtain the optimal solution, as shown in Figure 12.
As can be seen from the two figures above, the optimal result calculated by this model is influenced very little by the opportunity cost of the minor item, C o , and the inventory carrying cost, C o i , of major item i and only fluctuates within a very small range. The optimal cost calculated by the model presented in this section will not be affected by the order backlog rate, β i , of major item i . Through the analysis of the solution steps, it can be seen that once the β i exceeds the critical value, it will not affect the solution process, and then will not affect the optimal solution obtained by the model calculation.

5.7. Sensitivity Analysis of Model Optimal Cost to Sales Loss Rate

By keeping other parameters in the numerical examples unchanged, 10 different examples are respectively generated for the loss ratio,   λ i , of major item i to minor item sales in accordance with the range [0.1, 0.95] to obtain the optimal solution, as shown in Figure 13.
As can be seen from the figure above, the optimal cost calculated by the model proposed in this section has a great influence on the optimal cost in the interval λ i [ 0.8 ,   0.95 ] , and it is close to linear influence in other intervals.

6. Conclusions

In this paper, we consider solving the many-to-one cross-selling JRP problem, which leads to the certainty of partial late delivery and related demand. We extend the proposed similar model to make it capable of handling multiple replenishment cycles and the cross-selling of multiple major items with one minor item, allowing partial late delivery. In this paper, various possibilities are discussed under the condition of the demand satisfaction rate of major items’ F. The optimization analytic expression of the model is finally obtained by utilizing the convexity of cost function for F and using the first-order conditions in optimization theory.
The proposed method has been validated in numerical examples of the model, and the acquisition of accurate solutions provides strong evidence for the effectiveness of the proposed approach. The model and approach offer a competent solution for practical applications.
Additionally, if the replenishment cycle of major items is assumed to be consistent in real life, the decision space will often be reduced. Therefore, the relaxation of the conditions of an inconsistent replenishment cycle can also effectively increase the decision space and give the model proposed in this chapter a certain application prospect.
We also considered a widespread scenario in reality, that is, multiple major items versus one shared minor item under the cross-selling effect. From the perspective of practice, our paper may be able to provide some value to enterprises under appropriate conditions. Based on the numerical experiments in Section 5, we can see that some metrics of the minor items can significantly affect total ordering cost when compared to the major items. Specifically, for the ordering cost, demand rate, and inventory carrying cost, under certain conditions (being in a certain interval, controlling for other metrics) total ordering cost is strongly influenced by the minor items and weakly influenced by major items. The insight from this is that if an enterprise wants to reduce its total costs in a more efficient way, under similar conditions, it should focus more on adopting a means of controlling the indicators for minor items. To avoid sudden losses, the enterprise should prepare in advance to prevent large fluctuations caused by minor item suppliers.

Author Contributions

Conceptualization, R.Z.; Methodology, M.Y. and R.Z.; Writing—original draft, M.Y.; Writing—review & editing, M.Y.; Supervision, R.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data is contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Harris, F.W. How many parts to make at once. Oper. Res. 1990, 6, 947–950. [Google Scholar] [CrossRef]
  2. Lippman, S.A.; McCardle, K.F. The competitive newsboy. Oper. Res. 1997, 1, 54–65. [Google Scholar] [CrossRef]
  3. Khouja, M.; Goyal, S. A review of the joint replenishment problem literature: 1989–2005. Eur. J. Oper. Res. 2008, 1, 1–16. [Google Scholar] [CrossRef]
  4. Song, J.S. On the order fill rate in a multi-item inventory system. Oper. Res. 1998, 6, 831–845. [Google Scholar] [CrossRef]
  5. Whitin, T.M. The Theory of Inventory Management; Princeton University Press: Princeton, NJ, USA, 1957. [Google Scholar]
  6. Veinott, A.F. Optimal policy for a multi-product dynamic non-stationary inventory problem. Manag. Sci. 1965, 3, 206–222. [Google Scholar] [CrossRef]
  7. Huang, D.; Zhou, H.; Zhao, Q. A competitive multiple-product newsboy problem with partial product substitution. Omega 2010, 3, 302–312. [Google Scholar] [CrossRef]
  8. Pan, Q.; He, X.; Skouri, K.; Chen, S.-C.; Teng, J.-T. An inventory replenishment system with two inventory-based substitutable products. Int. J. Prod. Econ. 2018, 204, 135–147. [Google Scholar] [CrossRef]
  9. Pando, V.; San-José, L.A.; Sicilia, J. A new approach to maximize the profit/cost ratio in a stock-dependent demand inventory model. Comput. Oper. Res. 2020, 120, 104940. [Google Scholar] [CrossRef]
  10. Agrawal, R.; Imieliński, T.; Swami, A. Mining Association Rules between Sets of Items in Large Database. ACM SIGMOD Record. 1993, 2, 207–216. [Google Scholar] [CrossRef]
  11. Brijs, T.; Swinnen, G.; Vanhoof, K.; Wets, G. Using Association Rules for Product Assortment Decisions: A Case Study. In Proceedings of the KDD, San Diego, CA, USA, 15–18 August 1999. [Google Scholar]
  12. Wang, K.; Su, M.T. Item Selection by “Hub-Authority” Profit Ranking. In Proceedings of the SIGKDD’02, Edmonton, AB, Canada, 23–26 July 2002. [Google Scholar]
  13. Wong, R.C.; Fu, A.W.; Wang, K. MPIS: Maximal-Profit Item Selection with Cross-Selling Consideration. In Proceedings of the Third IEEE International Conference on Data Mining (ICDM’03), Melbourne, FL, USA, 22 November 2003. [Google Scholar]
  14. Liu, Y.; Xiao, Y.; Dai, Y. Omnichannel retailing with different order fulfillment and return options. Int. J. Prod. Res. 2023, 15, 5053–5074. [Google Scholar] [CrossRef]
  15. Yang, G.; Ji, G. The impact of cross-selling on managing consumer returns in omnichannel operations. Omega 2022, 111, 102665. [Google Scholar] [CrossRef]
  16. Ghoshal, A.; Mookerjee, V.S.; Sarkar, S. Recommendations and Cross-selling: Pricing Strategies when Personalizing Firms Cross-sell. J. Manag. Inf. Syst. 2021, 2, 430–456. [Google Scholar] [CrossRef]
  17. Mokhtari, H. Economic order quantity for joint complementary and substitutable items. Math. Comput. Simul. 2018, 154, 34–47. [Google Scholar] [CrossRef]
  18. Chen, J.; Feng, X.; Kou, G.; Mu, M. Multiproduct newsvendor with cross-selling and narrow-bracketing behavior using data mining methods. Transp. Res. Part E Logist. Transp. Rev. 2023, 169, 102985. [Google Scholar] [CrossRef]
  19. Netessine, S.; Savin, S.; Xiao, W.Q. Revenue management through dynamic cross-selling in e-commerce retailing. Oper. Res. 2006, 5, 893–913. [Google Scholar] [CrossRef]
  20. Zhang, R.-Q. An extension of partial backordering EOQ with correlated demand caused by cross-selling considering multiple minor items. Eur. J. Oper. Res. 2012, 3, 876–881. [Google Scholar] [CrossRef]
  21. Zhang, R.-Q.; Kaku, I.; Xiao, Y.-Y. Model and heuristic algorithm of the joint replenishment problem with complete backordering and correlated demand. Int. J. Prod. Econ. 2012, 1, 33–41. [Google Scholar] [CrossRef]
  22. Zhang, R.-Q.; Yi, M.; Wang, Q.-Q.; Xiang, C. Polynomial algorithm of the inventory model with complete backordering and correlated demand caused by cross-selling. Int. J. Prod. Econ. 2018, 199, 193–198. [Google Scholar] [CrossRef]
  23. Liu, R.; Zhou, Z.; Qin, Q.; Fan, B. Centralized drug procurement operation scheduling with a capacitated joint replenishment and delivery strategy: Evidence from China. Comput. Ind. Eng. 2022, 172, 108584. [Google Scholar] [CrossRef]
  24. Zhang, R.-Q.; Kaku, I.; Xiao, Y.-Y. Deterministic EOQ with partial backordering and correlated demand caused by cross-selling. Eur. J. Oper. Res. 2011, 3, 537–551. [Google Scholar] [CrossRef]
  25. Axsäter, S. Using the deterministic EOQ formula in stochastic inventory control. Manag. Sci. 1996, 42, 830–834. [Google Scholar] [CrossRef]
  26. Pentico, D.W.; Drake, M.J. The deterministic EOQ with partial backordering: A new approach. Eur. J. Oper. Res. 2009, 194, 102–113. [Google Scholar] [CrossRef]
Figure 1. The situation of Γ ¯ i ( F i ) when F i = 0 .
Figure 1. The situation of Γ ¯ i ( F i ) when F i = 0 .
Systems 12 00006 g001
Figure 2. The situation of Γ ¯ i ( F i ) when F i = 1 .
Figure 2. The situation of Γ ¯ i ( F i ) when F i = 1 .
Systems 12 00006 g002
Figure 3. Optimal total cost curve varying with ordering cost of minor items.
Figure 3. Optimal total cost curve varying with ordering cost of minor items.
Systems 12 00006 g003
Figure 4. Optimal total cost curve varying with ordering cost of major items.
Figure 4. Optimal total cost curve varying with ordering cost of major items.
Systems 12 00006 g004
Figure 5. Optimal total cost curve varying with the demand rate of minor item.
Figure 5. Optimal total cost curve varying with the demand rate of minor item.
Systems 12 00006 g005
Figure 6. Optimal total cost curve varying with the ordering cost of major item.
Figure 6. Optimal total cost curve varying with the ordering cost of major item.
Systems 12 00006 g006
Figure 7. Optimal total cost curve varying with the ordering cost of major item.
Figure 7. Optimal total cost curve varying with the ordering cost of major item.
Systems 12 00006 g007
Figure 8. Optimal total cost curve varying with the ordering cost of major item.
Figure 8. Optimal total cost curve varying with the ordering cost of major item.
Systems 12 00006 g008
Figure 9. Optimal total cost curve changing with inventory carrying cost of the major item.
Figure 9. Optimal total cost curve changing with inventory carrying cost of the major item.
Systems 12 00006 g009
Figure 10. Optimal total cost curve changing with the opportunity cost of minor item.
Figure 10. Optimal total cost curve changing with the opportunity cost of minor item.
Systems 12 00006 g010
Figure 11. Optimal total cost curve changing with inventory carrying cost of major item.
Figure 11. Optimal total cost curve changing with inventory carrying cost of major item.
Systems 12 00006 g011
Figure 12. Optimal total cost curve with the change of order backlog rate β i of major item i .
Figure 12. Optimal total cost curve with the change of order backlog rate β i of major item i .
Systems 12 00006 g012
Figure 13. Optimal total cost curve changing with the loss rate of sales of major item to minor item.
Figure 13. Optimal total cost curve changing with the loss rate of sales of major item to minor item.
Systems 12 00006 g013
Table 1. Variable Definition.
Table 1. Variable Definition.
DiDemand per unit of time for major item i per unit of time, in pieces/year
AiFixed order cost of major item i, in $/order
COiOpportunity cost of lost sales of major item i, in $/unit
ChiHolding cost per unit of time per unit of major item i, in $/piece/year
CbiCost per backorder of major item i per unit of time, in $/unit/year
βiBackorder rate of major item i
DjDemand per unit time of associated minor item j per unit time, in piece/year
AFixed order cost of minor item, in $/order
CoOpportunity cost of lost sales per unit of minor item per order cycle time, in $/piece
ChHolding cost per unit of minor item per unit of time, in $/piece/year
λiProportion of sales loss of minor items caused by sales loss of major item i
TOrder cycle
QiOrder volume
FiInventory/actual demand
Table 2. Table of interpretation of decision variables.
Table 2. Table of interpretation of decision variables.
T 0 Order cycle time of related products
T Order cycle time of major items, i.e., basic order cycle, in years
F Demand satisfaction rate of major items
Table 6. The critical value of β of each major item.
Table 6. The critical value of β of each major item.
i = 1 i = 2 i = 3
β i 0.99950.99980.9595
Table 7. The optimal T and F of each major item.
Table 7. The optimal T and F of each major item.
Optimal SolutionMajor Item
i = 1
Major Item
i = 2
Major Item
i = 3
T i 1.05412.47591.4173
F i 110.0541
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

Yi, M.; Zhang, R. Modeling and Solving the Joint Replenishment Problem with Cross-Selling Effects Considering One Shared Minor Item. Systems 2024, 12, 6. https://doi.org/10.3390/systems12010006

AMA Style

Yi M, Zhang R. Modeling and Solving the Joint Replenishment Problem with Cross-Selling Effects Considering One Shared Minor Item. Systems. 2024; 12(1):6. https://doi.org/10.3390/systems12010006

Chicago/Turabian Style

Yi, Meng, and Renqian Zhang. 2024. "Modeling and Solving the Joint Replenishment Problem with Cross-Selling Effects Considering One Shared Minor Item" Systems 12, no. 1: 6. https://doi.org/10.3390/systems12010006

APA Style

Yi, M., & Zhang, R. (2024). Modeling and Solving the Joint Replenishment Problem with Cross-Selling Effects Considering One Shared Minor Item. Systems, 12(1), 6. https://doi.org/10.3390/systems12010006

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