Integrated Production-Inventory Supply Chain Model

Abstract

In this paper, an integrated production-inventory model is presented for a supplier, manufacturer, and retailer supply chain under conditionally permissible delay in payments in uncertain environments. The supplier produces the item at a certain rate, which is a decision variable, and purchases the item to the manufacturer. The manufacturer has also purchased and produced the item in a finite rate. The manufacturer sells the product to the retailer and also gives the delay in payment to the retailer. The retailer purchases the item from the manufacture to sell it to the customers. Ideal costs of supplier, manufacturer, and retailer have been taken into account. An integrated model has been developed and solved analytically in crisp and uncertain environments, and finally, corresponding individual profits are calculated numerically and graphically.

Source: Dipak Kumar Jana, Kalipada Maity, and Tapan Kumar Roy, https://juaa-journal.springeropen.com/articles/10.1186/2195-5468-1-6
This work is licensed under a Creative Commons Attribution 2.0 License.

Introduction

Supply chain management has taken a very important and critical role for any company in the increasing globalization and competition in the market. A supply chain model (SCM) is a network of suppliers, producers, distributors, and customers which synchronizes a series of interrelated business process in order to have (1) optimal procurement of raw materials from nature, (2) transportation of raw materials into a warehouse, (3) production of goods in the production center, and (4) distribution of these finished goods to retailers for sale to the customers. With a recent paradigm shift to the supply chain (SC), the ultimate success of a firm may depend on its ability to link supply chain members seamlessly.

One of the earliest efforts to create an integrated SCM has been developed by Oliver and Webber, Cohen and Baghanan, and Cachon and Zipkin. They developed a production, distribution, and inventory (PDI) planning system that integrated three supply chain segments comprising supply, storage/location, and customer demand planning. The core of the PDI system was a network model and diagram that increased the decision maker's insights into supply chain connectivity. The model however was confined to a single-period and single-objective problem. Viswanathan and Piplani were concerned an integrated inventory model through common replenishment in the SC.  Khouja was the first to consider a three-stage supply chain with one or more firms at each stage. Agarwal et al. have developed a dynamic balancing of inventory model in supply chain management. Rau et al. developed an integrated SCM of a deteriorating item with shortages. Lee added a new dimension to the single vendor-single buyer problem by setting the number of raw material shipments received by the vendor per cycle to be a decision variable. Ben-Daya et al. have developed an integrated production-inventory model with raw material replenishment considerations in a three-layer supply chain. Sana has integrated a production-inventory model of imperfect quality products in a three-layer supply chain. Recently, Pal et al. have developed a three-layer supply chain model with production-inventory model for reworkable items. All of the abovementioned SCMs are considered with constant, known demand and production rates in a crisp environment.

Different types of uncertainty such as fuzziness, randomness, and roughness are common factors in SCM. In many cases, it is found that some inventory parameters involve fuzzy uncertainty. For example, inventory-related costs such as holding cost and setup cost, demand, and selling price depend on several factors such as bank interest, stock amount, and market situation which are uncertain in a fuzzy sense. To be more specific, inventory holding cost is sometimes represented by a fuzzy number, and it depends on the storage amount which may be imprecise and range within an interval due to several factors such as scarcity of storage space, market fluctuation, human estimation, and/or thought process. The following papers have been developed in these environments.

Wang and Shu developed a fuzzy decision methodology that provides an alternative framework to handle SC uncertainties and to determine SC inventory strategies, while there is a lack of certainty in data or even a lack of available historical data. Fuzzy set theory is used to model SC uncertainty. A fuzzy SC model based on possibility theory is developed to evaluate SC performances. Based on the proposed fuzzy SC model, a genetic algorithm approach is developed to determine the order-up-to levels of stock-keeping units in the SC to minimize SC inventory cost subject to the restriction of fulfilling the target fill rate of the finished product. The proposed model allows decision makers to express their risk attitudes and to analyze the trade-off between customer service level and inventory investment in the SC, so that better SC inventory strategies can be made.

Das et al. have presented a joint performance of an SC with two warehouse facilities in a fuzzy environment. A realistic two-warehouse and multi-collection-production-inventory model with constant/stock-dependent demand, defective production system, and fuzzy budget constraint has been formulated and solved in an SC context. Later Chen et al. developed a multi-criteria fuzzy optimization for locating warehouses and distribution centers in a supply chain network.

Peidro et al. developed a fuzzy linear programming model for tactical supply chain planning in a multi-echelon, multi-product, multi-level, multi-period supply chain network in a fuzzy environment. In this approach, the demand, process, and supply uncertainties are jointly considered. The aim is to centralize multi-node decisions simultaneously to achieve the best use of the available resources along the time horizon so that customer demands are met at a minimum cost. This proposal is tested using data from a real automobile SC. The fuzzy model provides the decision maker with alternative decision plans with different degrees of satisfaction.

Chu developed the supply chain flexibility that has become increasingly important. This study thus builds a group decision-making structure model of flexibility in supply chain management development. Recently, Jana et al. have developed a fuzzy simulation via contractive mapping genetic algorithm approach to an imprecise production-inventory model under volume flexibility. This study presents a framework to evaluate the supply chain flexibility that comprises two parts: (1) an evaluation hierarchy with flexibility dimensions and related metrics and (2) an evaluation scheme that uses a three-stage process to evaluate the supply chain flexibility. This study then proposes an algorithm to determine the degree of supply chain flexibility using a fuzzy linguistic approach. Evaluations of the degree of supply chain flexibility can identify the need to improve supply chain flexibility and identify specific dimensions of supply chain flexibility as the best directions for improvement. The results of this study are more objective and unbiased for two reasons. First, the results are generated by group decision-making with interactive consensus analysis. Second, the fuzzy linguistic approach used in this study has more advantages to preserve no loss of information than other methods. Additionally, this study presents an example using a case study to illustrate the availability of the proposed methods and compare it with other methods.

Kristianto et al. developed an adaptive fuzzy control application to support a vendor-managed inventory (VMI). This paper also guides the management in allocating inventory by coordinating with suppliers and buyers to ensure minimum inventory levels across a supply chain. Adaptive fuzzy VMI control is the main contribution of this paper.

However, the uncertainty theory was developed by Liu, and it can be used to handle subjective imprecise quantity. Much research work has been done on the development of the uncertainty theory and related theoretical work. You proved some convergence theorems of uncertain sequences. Liu has defined uncertain process and Liu has discussed uncertain theory. In this paper, we developed for the first time a three-layer supply chain model under delay in payment in an uncertain environment.

In the traditional economic order quantity (EOQ) model, it often assumed that the retailer must pay off as soon as the items are received. In fact, the supplier offers the retailer a delay period, known as trade credit period, in paying for the purchasing cost, which is a very common business practice. Suppliers often offer trade credit as a marketing strategy to increase sales and reduce on-hand stock level. Once a trade credit has been offered, the amount of the tied up retailer's capital in stock is reduced, and that leads to a reduction in the retailer's holding cost of finance. In addition, during the trade credit period, the retailer can accumulate revenues by selling items and by earning interests. As a matter of fact, retailers, especially of small businesses which tend to have a limited number of financing opportunities, rely on trade credit as a source of short-term funds. In this research field, Goyal was the first to establish an EOQ model with a constant demand rate under the condition of permissible delay in payments. Khanra, Ghosh, and Chaudhuri have developed an EOQ model for a deteriorating item with time-dependent quadratic demand under permissible delay in payment. Also, Maihami and Abadi have established joint control of inventory and its pricing for non-instantaneously deteriorating items under permissible delay in payments and partial backlogging.

The proposed model considers a three-layer supply chain involving the supplier, manufacturer, and retailer who are responsible in transforming the raw materials into finished product and making them available to satisfy the customer's demand time. Inventory and production decisions are made at the supplier, manufacturer, and retailer levels in uncertain environments. The problem is to coordinate production and inventory decisions across the supply chain so that the total profit of the chain is maximized.

Necessary knowledge about uncertain variables

Assumptions and notations

Model description and diagrammatic representation

Mathematical formulation of the model

Formulation of suppliers' individual average profit

Differential equation for the supplier in Figure 2 in [0,T] is given by

$\dfrac{dq_s}{dt} =\left\{ \begin{array}{lll} p_s−p_m, & {0 ≤ t ≤ t_s} \\ −p_m, & {t_s ≤ t ≤ T_s} \\ 0, & {T_s ≤ t ≤ T} \end{array} \right.$

with boundary conditions $q_s (t) = 0$ and $t = 0,T_s$. Solving the differential equation with the boundary condition, we have

${q_s}(t) =\left\{ \begin{array}{lll} (p_s−p_m)t, & 0 ≤ t ≤ t_s \\ p_m (T_s - t), & t_{s} ≤ t ≤ T_s \\ 0, & T_s ≤ t ≤ T \end{array} \right.$

(3)

$H_s = \text{Holding cost of supplier}$

$= h_s [∫^{ts}_0 (p_s−p_m) tdt + ∫^{Ts}_{t_s} p_m (T_s−t) dt ]$

$=h_s [\dfrac{p_s t^{2}_{s}} {p_m} −p_s t^{2}_{s}]$.

Figure 2 Inventory level of supplier.

The total idle cost = $id_s [TR+P_s t_s(\dfrac{1}{D_c}−\dfrac{1}{p_m})]$, purchase cost = $c_s p_m T_s$, selling price = $c_m p_m T_s$, and ordering cost = $A_s$.

$APS = \text{Average profit for supplier}$

$=\dfrac{1}{T} [\text{revenue from sale-purchase cost-holding cost-idle cost-ordering cost}]$

$=\dfrac{1}{T} [(c_m−c_s) p_s t_s − h_s (\dfrac{p^{2}_{s} t^{2}_{s}} {p_m} − p_s t^{2}_{s}) − id_s (T_R + p_s t_s (\dfrac{1}{D_c} −\dfrac{1}{p_m}))−A_s]$

(4)

Formulation of manufacturer individual average profit

Inventory level of manufacturer in Figure 3 in [0,T] is given by

${q_m}(t) =\left\{ \begin{array}{lllllll} p_mt, & 0 ≤ t ≤ T_R \\ p_mT - iD_R, & iT_R < t ≤ (i+1) \quad i = 1, 2, ..., (r-1)\\ p_mT - rD_R, & rT_R < t ≤ T_s \\ p_m T_s - rD_r, & T_{s} < t ≤ (r+1) T_R \\ p_mT_s - iD_R, & iT_R < t ≤ (i+1) T_R \quad i = r+1, r+2, ..., (n-1)\\ p_m T_s - nD_R, &nT_R < t ≤ (n+1) T_R\\ 0, & (n+1)T_R < t ≤ T \end{array} \right.$

(5)

with boundary conditions $q_m (0) = 0$ and $q_m (i T _R  + 0) = q_m (i T_R ) − D_R$.

Figure 3 Inventory level of manufacturer.

$H_m=\text{Holding cost for manufacturer}$.

$=hm[∫^{TR}_{0}p_mtdt+∑^{1}_{r−1}∫^{(i+1)T_R}_{iT_R}(p_mt−iD_R)dt+∫^{(T_s}_{rT_R}(p_mt−rD_R)dt \\ + ∫^{(r+1)T_R}_{T_s}(p_mT_s−rD_R)dt+∑_{r+1}^{n−1}∫^{(i+1)T_R}_{iT_R} (p_mT_s−iD_R)dt+∫^{(n+1)T_R}_{nT_R} (p_mT_s−nD_R)dt] \\ =h_m[np_mT_sT_R−\dfrac{n^2+n−2r−2}{2}T_RD_R−\dfrac{p^{2}_{s}t^{2}_{s}}{2p_m}]$.

The total idle cost = $id_m[\dfrac{pm_{Tm}−nD_R}{D_c}]$, purchase cost =$c_m p_m T_s$, selling price = $c_r p_m T_s$, and ordering cost = $A_m$.

Case 1

(When $M ≤ T′ ≤ T R$)

$I_{em}=I_{pr}=\text{Amount of interest earned by the manufacturer in [0,T] from retailer}$.

$= \text{Amount of interest paid by the retailer to the manufacturer in [0,T]}$.

$=c_rI_p[n∫^{T_R}_{M}(D_R−D_ct)dt+∫^{T′}_{M}(p_mT_s−nD_R−D_ct)dt]$

$=\dfrac{nc_rI_p}{2}[T_RD_R+D_cM_2−2MDR]+c_rI_p[(\dfrac{(p_mT_s−nD_R)^2}{2D_c} + (p_mT_s−nD_R)M+\dfrac{D_cM^2}{2}]$

$APM_1=\text{Average profit of manufacturer}$.

$=\dfrac{1}{T} \text{[revenue from sale-purchase cost-holding cost-idle cost+earned interest-ordering cost.]}$

$=\dfrac{1}{T}[(c_r−c_m)p_mT_h_m(np_mT_sT_R−\dfrac{n^2+n−2r−2}{2}T_RD_R−\dfrac{p^{2}_{s}t^{2}_{s}}{2pm}) \\ −id)m(\dfrac{p_mT_m−nD_R}{D_c})+\dfrac{nc_rI_p}{2}[T_RD_R+D_cM^2−2MD_R] \\ +c_rI_p(\dfrac{(p_mT_s−nD_R)^2}{2D_c}+(p_mT_s−nD_R)M+\dfrac{D_cM^2}{2})−Am]$
(6)

Case 2

(When $T′ ≤ M ≤ T R$)
$I_{em}=I_{pr}=\text{Amount of interest earned by the manufacturer in [0,T]  from retailer.}$

$=\text{Amount of interest paid by the retailer to the manufacturer in [0,T.]$</div><div style="text-align: left;">$=c_rI_p[n∫^{T_R}_{M}(D_R−D_ct)dt+∫^{T′}_{M}(p_mT_s−nD_R−D_ct)dt]$

$=\dfrac{nc_rI_p}{2}[T_RD_R+D_cM^2−2MD_R]+c_rI_p[(\dfrac{(p_mT_s−nD_R)^2}{2D_c} \\+(p_mT_s−nD_R)M+\dfrac{D_cM^2}{2}]$

$APR_2=\text{Average profit of retailer.}$

$=\dfrac{1}{T} \text{[revenue from sale-purchase cost-holding cost+earned interest-payable interest-idle cost-ordering cost].}$

$=\dfrac{1}{T}[(cr_1−c_r)p_mT_s−\dfrac{h_r}{2}(\dfrac{p^{2}_{m}T^{2}_{s}}{D_c}−2np_mT_sT_R−(2n+1)T_RD_R) \\ +\dfrac{nc_{r+1}i_eD_cM^2}{2}+ \dfrac{c_{r_1}{i_e}}{2} (p_mT_s−nD_R)(2M−T′) \\ −\dfrac{nc_rI_p}{2}[T_RD_R+D_cM_2−2MD_R]−id_rT_R−A_r]$
(7)

Formulation of retailer individual average profit

Inventory level of retailer in Figure 4 in [0,T] is given by

$q_r(t)= \left\{ \begin{array}{ll} D_ct, & iT_R ≤ t ≤ (i+1)T_R  \\ p_mT_s nD_r - d_ct, & (n+1)T_R ≤ t ≤ T \end{array} \right.$

(8)
with boundary conditions $q _r ((n + 1)T_R ) = 0$ and $q_r (T) = 0$.

Figure 4 Inventory level of retailer.

$Hr = \text{Holding cost of retailer.}$

$=nh_r[∫^{T_R}_{0}(D_R−D_ct)dt+∫^{T′}_{0}(p_mT_s−nD_R−D_ct)dt]$

$=\dfrac{2}[\dfrac{p^{2}_{m}T^{2}_{s}}{D_c}−2np_mT_sT_R−(2n+1)T_RD_R]$.

The total idle cost = $id_r T_R$, purchase cost = $c_r p_m T_s$, selling price = $c_{r1} p_m T_s$, and ordering cost = $A_r$.

Case 1

(When $M ≤ T′ ≤ T R$)

Interest earned by the retailers for ($n + 1$) cycle is given by

$I_{er}=\text {Amount of interest earned by the retailer from the bank in} (n+1) \text {cycle}$.

$=(n+1)c_{r1}i_e[∫^{M}_{0}(M−t)D_cdt]$

$=\dfrac{(n+1)c_{r1}i_eD_cM^2}{2}$,

$I_{pr}=\text{Amount of interest paid by the retailer to the manufacturer in [0,T]}$.

$=c_rI_p[n∫^{T_R}{M}(D_R−D_ct)dt+∫^{T′}_{M}(p_mT_s−nD_R−D_ct)dt]$

$=nc_rI_p[\dfrac{T_RD_R+D_cM^2−2MD_R}{2}]+c_rI_p[\dfrac{(p_mT_s−nD_R)^2}{2D_c} \\ +(p_mT_s−nD_R)M+\dfrac{D_cM^2}{2}]$,

$APR_1=\text{Average profit for retailer}$.

$=\dfrac{1}{T} [\text {revenue from sale-purchase cost-holding cost +earned interest-payable interest-idle cost-ordering cost.]}$

$\dfrac{1}{T}[c_{r1}p_mT_s−c_rp_mT_s−\dfrac{h_r}{2}(\dfrac{p^{2}_{m}T^{2}_{s}}{D_c}−2np_mT_sT_R−(2n+1)T_RD_R) \\ +\dfrac{(n+1)c_{r1}i_eD_cM^2}{2}−\dfrac{n_{cr}I_p}{2}[T_RD_R+D_cM^2−2MD_R]−id_rT_R−A_r \\ +(p_mT_s−nD_R)M+\dfrac{D_cM^2}{2})−c_rI_p(\dfrac{(p_mT_s−nDR)^2}{2D_c}]$
(9)

Case 2

(When $T′ ≤ M ≤ T R$)

Interest earned by the retailers for ($n + 1$) cycle is given by

$I_{re}=\text{Retailers' earned interest.}$

$\begin{aligned}&=c_{r_{1}} i_{e}\left[n \int_{0}^{M}(M-t) D_{c} d t+\int_{0}^{T^{\prime}}\left(T^{\prime}-t\right) D_{c} d t+\left(M-T^{\prime}\right)\left(p_{m} T_{s}-n D_{R}\right)\right] \\&=\frac{n c_{r_{1}} i_{e} D_{c} M^{2}}{2}+\frac{c_{r_{1}} i_{e}}{2}\left(p_{m} T_{s}-n D_{R}\right)\left(2 M-T^{\prime}\right)\end{aligned}$

Interest payable by the retailers for the first n cycle is given by
$I_{rp}=\text{Retailers' payable interest.}$

$\begin{aligned}&=c_{r} I_{p}\left[n \int_{M}^{T_{R}}\left(D_{R}-D_{c} t\right) d t\right] \\&=\frac{n c_{r} I_{p}}{2}\left[T_{R} D_{R}+D_{c} M^{2}-2 M D_{R}\right]\end{aligned}$

(10)

$APR_2=\text{Average profit for retailer.}$

$$=\frac{1}{T}[$ revenue from sale-purchase cost-holding cost+earned interest-payable interest-idle cost-ordering cost]\(\begin{aligned}=& \frac{1}{T}\left[\left(c_{r_{1}}-c_{r}\right) p_{m} T_{s}-\frac{h_{r}}{2}\left(\frac{p_{m}^{2} T_{s}^{2}}{D_{c}}-2 n p_{m} T_{s} T_{R}-(2 n+1) T_{R} D_{R}\right)\right.\\&+\frac{n c_{r_{1}} i_{e} D_{c} M^{2}}{2}+\frac{c_{r_{1}} i_{e}}{2}\left(p_{m} T_{s}-n D_{R}\right)\left(2 M-T^{\prime}\right) \\&\left.-\frac{n c_{r} I_{p}}{2}\left[T_{R} D_{R}+D_{c} M^{2}-2 M D_{R}\right]-\mathrm{id}_{r} T_{R}-A_{r}\right]\end{aligned}$\)

(11)

Crisp environment

Case 1

($M ≤ T′ ≤ T R$)
$ATP_1=\text{Total average profit for integrated model}$

$=AP_S+APM_1+APR_1$

$\begin{aligned}=& \frac{D_{c}}{p_{m} T_{s}}+D_{R}\left[\left(c_{m}-c_{s}\right) p_{m} T_{s}-h_{s}\left(\frac{p_{s} t_{s}^{2}}{p_{m}}-p_{s} t_{s}^{2}\right)-\mathrm{id}_{s}\left(T_{R}+p_{s} t_{s} \frac{1}{D_{c}}-\frac{1}{p_{m}}\right)\right.\\&+\left(c_{r}-c_{m}\right) p_{m} T_{s}-h_{m}\left(n p_{m} T_{s} T_{R}-\frac{n^{2}+n-2 r-2}{2} T_{R} D_{R}-\frac{p_{s}^{2} t_{s}^{2}}{2 p_{m}}\right) \\&-\mathrm{id}_{m}\left(\frac{p_{m} T_{m}-n D_{R}}{D_{c}}\right) \\&-A_{m}+\left(c_{r_{1}}-c_{r}\right) p_{m} T_{s}-\frac{h_{r}}{2}\left(\frac{p_{m}^{2} T_{s}^{2}}{D_{c}}-2 n p_{m} T_{s} T_{R}-(2 n+1) T_{R} D_{R}\right) \\&\left.+\frac{(n+1) c_{r_{1}} I_{e} D_{c} M^{2}}{2}-\mathrm{id}_{r} T_{R}-A_{r}\right] \\=& \frac{D_{c}}{p_{m} T_{s}+D_{R}}\left[A p_{s}^{2} t_{s}^{2}-B p_{s} t_{s}+E\right]\end{aligned}$

(12)

where $P_m T_s  = P_s t_s$ , $A=\dfrac{h_m−h_s}{2p_m}−\dfrac{h_r}{2D_c}$ (13)

$B=(c_{r1}−c_s)+n(h_r−h_m)T_R−\dfrac{id_m+id_s}{D_c}$ (14)

and

$E=[h_m \dfrac{n^2+n−2r−2}{2}+h_r \dfrac{2n+1}{2}]T_RD_R+\dfrac{(n+1)c_{r1}I_eD_cM^2}{2} \\ +(nid_m−id_s−id_r)T_R+id_sT_s−(A_s+A_m+A_r$) (15)

 $\begin{aligned}\frac{d}{d p_{s}}\left(\mathrm{ATP}_{1}\right) &=0 \\\Rightarrow p_{s} &=\frac{2 D_{R}-B \pm \sqrt{\left(2 t_{s}-A\right)^{2}-4\left(2 t_{s}-A\right)\left(E-B D_{R}\right)}}{2\left(2 t_{s}-A\right)}\end{aligned}$ (16)

$\frac{d^{2}}{d p_{s}^{2}}\left(\mathrm{ATP}_{1}\right)$ (17)

Therefore, $ATP_1$ is concave if $2Ap_s+B < 4p_st_s+Bt_s+2D_R$. (18)

Case 2

($T′ ≤ M ≤ T R$)
$ATP_2= \text{Total average profit for integrated model.}$

$\begin{aligned}=& \frac{D_{c}}{p_{m} T_{s}+D_{R}}\left[\left(c_{m}-c_{s}\right) p_{m} T_{s}-h_{s}\left(\frac{p_{s} t_{s}^{2}}{p_{m}}-p_{s} t_{s}^{2}\right)-\mathrm{id}_{s}\left(T_{R}+p_{s} t_{s} \frac{1}{D_{c}}-\frac{1}{p_{m}}\right)-A_{s}\right.\\&+\left(c_{r}-c_{m}\right) p_{m} T_{s}-h_{m}\left(n p_{m} T_{s} T_{R}-\frac{n^{2}+n-2 r-2}{2} T_{R} D_{R}-\frac{p_{s}^{2} t_{s}^{2}}{2 p_{m}}\right) \\&-\operatorname{id}_{m}\left(\frac{p_{m} T_{m}-n D_{R}}{D_{c}}\right) \\&-A_{m}+\left(c_{r 1}-c_{r}\right) p_{m} T_{s}-\frac{h_{r}}{2}\left(\frac{p_{m}^{2} T_{s}^{2}}{D_{c}}-2 n p_{m} T_{s} T_{R}-(2 n+1) T_{R} D_{R}\right) \\&\left.+\frac{n c_{r_{1}} I_{e} D_{c} M^{2}}{2}+\frac{c_{r_{1}} I_{e}}{2}\left(p_{m} T_{s}-n D_{R}\right)(2 M-T)-\mathrm{id}_{r} T_{R}-A_{r}\right] \\=& \frac{D_{c}}{p_{m} T_{s}+D_{R}}\left[A p_{s}^{2} t_{s}^{2}-B p_{s} t_{s}+F\right]\end{aligned}$

(19)

where A and B are given in (16) and (17), respectively and

F=[hmn2+n−2r−22+hr2n+12]TRDR+(n+1)cr1IeDcM22

+(nidm−ids−idr)TR+idsTs−(As+Am+Ar) (20)

$\begin{aligned}\frac{d}{d p_{s}}\left(\mathrm{ATP}_{2}\right) &=0 \\\Rightarrow p_{s} &=\frac{2 D_{R}-B \pm \sqrt{\left(2 t_{s}-A\right)^{2}-4\left(2 t_{s}-A\right)\left(F-B D_{R}\right)}}{2\left(2 t_{s}-A\right)}\end{aligned}$ (21)

$\frac{d^{2}}{d p_{s}^{2}}\left(\mathrm{ATP}_{1}\right)$ (22)

Therefore, $ATP_1$ is concave if $2Ap_s+B < 4p_st_s+Bt_s+2D_R$ (23)

Proposed inventory model in uncertain environment

Let us consider $\widetilde {id}_s$, $\widetilde{id}_r$, $\widetilde{id}_m$, and $\widetilde{I}_{re}$ as zigzag uncertain variables where $\widetilde{id}_s=L(id_{s1},id_{s2},id_{s3})$, $\widetilde{id}_r =L(id_{r1},id_{m2},id_{r3})$, $\widetilde{id}_m= L(id_{m1},id_{m2},id_{m3})$, and $\widetilde{I}_{re}=L(I_{re1},I_{re2},I_{re3})$. Then, the objective is reduce to the following:

•  For Case 1 ($M ≤ T′ ≤ T R$ )

$A\widetilde{T}P_1 =\dfrac{D_c}{p_mT_s+D_R}[0(c_m−c_s)p_mT_s−h_s(\dfrac{p_st^{2}_{s}}{p_m}−p_st^{2}_{s})−\widetilde{id}_s(T_R+p_st_s(\dfrac{1}{D_c}−\dfrac{1}{p_m})) \\ \begin{aligned} &-A_{s}+\left(c_{r}-c_{m}\right) p_{m} T_{s}-h_{m}\left(n p_{m} T_{s} T_{R}-\frac{n^{2}+n-2 r-2}{2} T_{R} D_{R}-\frac{p_{s}^{2} t_{s}^{2}}{2 p_{m}}\right) \\ &-A_{m}+\left(c_{r_{1}}-c_{r}\right) p_{m} T_{s}-\frac{h_{r}}{2}\left(\frac{p_{m}^{2} T_{s}^{2}}{D_{c}}-2 n p_{m} T_{s} T_{R}-(2 n+1) T_{R} D_{R}\right) \\ &\left.-\widetilde{i d}_{m}\left(\frac{p_{m} T_{m}-n D_{R}}{D_{c}}\right)+\frac{(n+1) c_{r_{1}} \widetilde{I}_{r e} D_{c} M^{2}}{2}-\widetilde{i d}_{r} T_{R}-A_{r}\right] \end{aligned}$

(24)

•  For Case 2 ($T′ ≤ M ≤ T R$)

$A\widetlide{T}P_2 = \frac{D_{c}}{p_{m} T_{s}+D_{R}}\left[\left(c_{m}-c_{s}\right) p_{m} T_{s}-h_{s}\left(\frac{p_{s} t_{s}^{2}}{p_{m}}-p_{s} t_{s}^{2}\right)-\tilde{i d}_{s}\left(T_{R}+p_{s} t_{s} \frac{1}{D_{c}}-\frac{1}{p_{m}}\right)−A_s \right. \\\begin{aligned}&+\left(c_{r}-c_{m}\right) p_{m} T_{s}-h_{m}\left(n p_{m} T_{s} T_{R}-\frac{n^{2}+n-2 r-2}{2} T_{R} D_{R}-\frac{p_{s}^{2} t_{s}^{2}}{2 p_{m}}\right) \\&-\widetilde{i d}_{m}\left(\frac{p_{m} T_{m}-n D_{R}}{D_{c}}\right) \\&-A_{m}+\left(c_{r 1}-c_{r}\right) p_{m} T_{s}-\frac{h_{r}}{2}\left(\frac{p_{m}^{2} T_{s}^{2}}{D_{c}}-2 n p_{m} T_{s} T_{R}-(2 n+1) T_{R} D_{R}\right) \\&\left.+\left\{\frac{n c_{r_{1}} D_{c} M^{2}}{2}+\frac{c_{r_{1}}}{2}\left(p_{m} T_{s}-n D_{R}\right)\left(2 M-T^{\prime}\right)\right\} \widetilde{I}_{r e}-\widetilde{i d}_{r} T_{R}-A_{r}\right]\end{aligned}$

(25)

The equivalent crisp model

Using Lemma 1 and applying Theorem 2, the expected total average profit is given by the following:

•  For Case 1 $(M ≤ T′ ≤ T R )$
  $E A [\widetilde {T} P_1] =\dfrac{D_c}{pm^{T_s+D_R}} [0(c_m−c_s)p_mT_s−h_s(\dfrac{p_{s}t^{2}_{s}} {p_m} −p_st^{2}_{s}) \\ − E[\widetilde{id}_s](T_R+p_st_s(\dfrac{1}{D_c}−\dfrac{1}{p_m})) \\ − A_s+(c_r−c_m)p_mT_s−h_m(np_mT_sT_R−\dfrac{n_2+n−2r−2}{2}T_RD_R−\dfrac{p^{2}_{s}t^{2}_{s}}{2p_m}) \\ − A_m+(cr_1−c_r)p_mT_s−\dfrac{h_r}{2}(\dfrac{p^{2}_{m}T^{2}_{s}}{D_c}−2np_mT_sT_R−(2n+1)T_RD_R)\\ − E[\widetilde{id}_m](\dfrac{p_mT_m−nD_R}{D_c})+\dfrac{(n+1)cr_1E[\widetilde{Ire}]D_cM^2}{2}−E[\widetilde{id}_r]T_R−A_r]$
  

(26)

•  For Case 2 (T′ ≤ M ≤ T R )

$E A [\widetilde {T} P_2] =\dfrac{D_c}{pm^{T_s+D_R}} [(c_m−c_s)p_mT_s−h_s(\dfrac{p_{s}t^{2}_{s}} {p_m} −p_st^{2}_{s}) −E[\widetilde {id}_s](T_R+p_st_s \dfrac{1}{D_c} −\dfrac{1}{p_m})−A_s \\ + (c_r−c_m)p_mT_s−h_m(np_mT_sT_R−\dfrac{n^2+n−2r−2}{2}T_RD_R−\dfrac{p^{2}_{s}t^{2}_{s}} {2p_m}) \\ −E[\widetilde {id}_m](\dfrac{p_mT_m−nD_R}{D_c}) −A_m+(cr_1−c_r)p_mT_s \\ −\dfrac{h_r}{2} (\dfrac{p^{2}_{m}T^{2}_s}{D_c} −2np_mT_sT_R−(2n+1)T_RD_R) \\ +{\dfrac{nc_{r1}D_cM^2}{2}+\dfrac{c_{r1}}{2}(p_mT_s−nD_R)(2M−T′)}E[\widetilde {Ire}]−E[\widetilde {id}r]T_R−A_r]$

(27)