Next Article in Journal
The Impact of Institutional Dimensions on Entrepreneurial Intentions of Students—International Evidence
Next Article in Special Issue
A Discount Technique-Based Inventory Management on Electronics Products Supply Chain
Previous Article in Journal
Multi-Factorized Semi-Covariance of Stock Markets and Gold Price
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Sustainable Economic Recycle Quantity Model for Imperfect Production System with Shortages

1
Department of Mechanical and Industrial Engineering, Prince Sattam Bin Abdulaziz University, Al-Kharj 16278, Saudi Arabia
2
Department of Mathematics, Jahangirnagar University, Saver, Dhaka 1342, Bangladesh
3
Department of Mathematics, Hajee Mohammad Danesh Science and Technology University, Dinajpur 5200, Bangladesh
4
Industrial and Systems Engineering Department, Chung Yuan Christian University, No. 200, Chung Pei Road, Chungli, Taoyuan City 32023, Taiwan
5
Kishore Bharati Bhagini Nivedita College, Ramkrishna Sarani, Behala, Kolkata 700060, India
6
Department of Mathematics and Statistics, University of Victoria, Victoria, BC V8W 3R4, Canada
7
Department of Medical Research, China Medical University Hospital, China Medical University, Taichung 40402, Taiwan
8
Department of Mathematics and Informatics, Azerbaijan University, 71 Jeyhun Hajibeyli Street, AZ1007 Baku, Azerbaijan
9
Section of Mathematics, International Telematic University Uninettuno, I-00186 Rome, Italy
*
Author to whom correspondence should be addressed.
J. Risk Financial Manag. 2021, 14(4), 173; https://doi.org/10.3390/jrfm14040173
Submission received: 6 March 2021 / Revised: 4 April 2021 / Accepted: 6 April 2021 / Published: 10 April 2021
(This article belongs to the Special Issue Sustainable Mathematical Modelling in Business Analysis)

Abstract

:
Recycling of products has a great impact on contemporary sustainable business strategies. In this study, a sustainable recycling process in a production-inventory model for an imperfect production system with a fixed ratio of recyclable defective products is introduced. The piecewise constant demand rates of the non-defective items are considered under production run-time, production off-time with positive stock, and production off-time with shortages under varying conditions. Based on the production process, two cases are studied using this model. The first case does not consider recycling processes, while the second case picks up all defective items before sending these items to recycling during the production off-time; the recycled items are added to the main inventory. The aim of this study is to minimize the total cost and identify the optimal order quantity. The manufacturing process with the recycling process provides a better result compared to without recycling in the first case. Some theoretical derivations are developed to enunciate the objective function using the classical optimization technique. To validate the proposed study, sensitivity analysis is performed, and numerical examples are given. Finally, some managerial insights and the scope of future research are provided.

1. Introduction and Literature Review

Inventory management control monitors and determines the optimal inventory order to minimize the total inventory cost. In every business, retailers often implement numerous strategies to intensify the profitability and sustainability of the system (Mashud et al. 2021c; Popescu and Popescu 2019). The familiar economic order quantity (EOQ) model was first initiated by Harris (1915). Later, Alharkan et al. (2020) worked on a two-echelon supply chain model with lot-size-dependent lead-time, and Taft (1918) extended the EOQ model to consider incremental orders during the production process. Inventory of goods is essential to the production processes to meet customer demand. The standard economic production quantity (EPQ) models typically assume perfect product quality and perfect production processes. In reality, however, ideal production processes are rare, and defective items will occur. Many industries recycle these faulty items. Some examples of recyclable materials are jute, paper, glass, and cotton and other fabrics.
During the last two decades, many researchers have studied inventory models with imperfect production systems. Defective items are mainly generated due to faulty production process (Rosenblatt and Lee 1986; Hsu and Hsu 2014). Salameh and Jaber (2000) attempted to modify the traditional EPQ model to consider poor-quality items, while Chiu et al. (2014) investigated lot-sizing problems with arbitrary defective rates in finite manufacture systems under rework processes and shortages. Defective items also deteriorate in some cases and emit carbon into the environment (Daryanto and Wee 2019; Mashud et al. 2020a). Despite deterioration, some products are wasted during the production process in the whole supply chain. Oh and Hwang (2006) anticipated an inventory model considering recycled raw materials without shortage. They assumed unequal holding costs for the raw materials and serviceable objects. Benkherouf and Omar (2017) investigated optimal production batch sizing with rework procedures, while Mokhtari (2018) outlined an order lot size combination between producers’ and suppliers’ defective manufacturing systems considering production rework processes. Al-Salamah (2019) designed an EPQ model for defective manufacturing procedures using synchronous and asynchronous elastic rework rates. A model for a long-run single-stage manufacturing system with imperfect items was developed by Kang et al. (2018) who incorporated rework operation, the inspection process, and planned backordering. In that model, the unfulfilled customer demands during the production phase due to process imperfections were satisfied either at the end of the inspection process or after reworking the imperfect products. Hasan et al. (2020) provided a model for agricultural products where to prevent product deterioration, they used product separation and offered different discount policies, while Nobil et al. (2020) presented a model for imperfect production process with reorder point in an EOQ model. An imperfect multi-stage production system considering defective proportions in the production process and uncertain product demand, and incorporating lean philosophy, was introduced by Tayyab et al. (2019b) who tried to reduce system costs. Onwude et al. (2020) analyzed a food wastage supply chain for fresh agricultural products. However, none of these models reworked defective items at each production stage. This model assumes the reworking of defective items at each production stage to bridge the aforementioned gap, and consequently the defective items (scraped) are perfect after reworking. Very few studies have been carried out with consideration of the recycling process with shortages, so the mentioned model also tries to incorporate shortages in the proposed ERQ model.
Demand is a fundamental key factor in inventory models because it determines the perplexity of optimal results. Sometimes demand becomes uncertain (Mashud et al. 2020b) and sometimes it depends on the price of the products (Mashud et al. 2021a, 2021c). Fuzzy demand is also very popular among the practitioners of inventory research (Tayyab et al. 2019a). Bai and Varanasi (1996) formulated an optimal control problem with indefinite piecewise constant demand and developed an algorithm to determine the optimal production strategy using Pontryagin’s minimum principle. Hsu and Hsu (2016) considered a model for defective items with stock-reliant demand with inflation for a Weibull distribution deterioration. Palanivel and Uthayakumar (2016) modified their model by adding payment delays. Viji and Karthikeyan (2016) developed a model for a three-level production process for Weibull distribution deterioration. Recently, Manna et al. (2017) offered an EPQ model for a production rate-reliant defective rate of items under advertisement with demand dependent on a defective production system.
In the collected works of imperfect production systems, most models were developed considering reworks, repair, remanufacture, and different production and supply policies. Some researchers considered two-level and three-level production rates, but production inventory models of three-level constant demand patterns are not available in the literature survey. The main contributions of the proposed model are:
(i)
The proposed model is developed focusing on the recycling of defective items collected from regular production after the proper screening with 100% recovery of raw materials used for production in the corresponding item category.
(ii)
An inventory model is developed for three-level piecewise constant demand, which varies under three different production time parts, from production run-time to production off-time with positive stock and production off-time with shortages.
The remainder of the paper is organized as follows: Notation and assumptions of the model are described in Section 2. Section 3 presents the mathematical formulation of the models. Theoretical results along with optimal solution uniqueness are shown in Section 4. Section 5 comprises numerical examples and convexity graphs. Sensitivity with cost–benefit is analyzed in Section 6. Managerial implementations are delivered in Section 7. Lastly, the conclusion and future research are discussed in Section 8.

2. Assumptions and Notation

In this section, some assumptions and notation necessary to the development of the paper are presented.

2.1. Assumptions

The following assumptions are worthy of mention in formulating the models discussed.
(i)
The production rate is finite and fixed.
(ii)
The defective rate is known, and constant and defective items are produced randomly alongside the perfect product.
(iii)
The demand rate of the perfect product is a piecewise constant function (motivated from Bai and Varanasi (1996)):
The demand during production-run and production off-times and stock-out period is as follows:
{ D ;                                 during   production run   time x D ;     during   production off   time   with   positive   stock y D ;           during   production off   time   with   shortages
(iv)
The sum of the demand rate and defective rate is less than the rate of production.
(v)
Lead-time is negligible, and the number of shortages is acceptable while the model is considered for a single items production system.
(vi)
Defective items are recyclable with 100% recovery of raw material usable for the production of same product in the next cycle time, and the value of the recycled materials is higher than that of purchased materials.
(vii)
The holding cost per defective item is same as that of fresh items, while the recycle cost per defective item is equal irrespective of its number.

2.2. Notation

The following notations (Table 1) are used in the development of the model:

3. Model Formulation

During the production run-time of cycle time T, a lot size of Q units is produced at a production rate P, and with it, W units of defective items are produced at a defective rate d. The production process runs during [0, K], and during the production-off time, the stock depletes due to customer demand. At time C, the stock becomes zero, and shortages start until G. Current demand is expressed at the demand rate “D” during production run-time, which is calculated again using a ratio coefficient to obtain a rate “ x D ” during production off-time while the stock is positive. Shortages also occur at the demand rate “ y D ” during the time t 3 while the stock is negative. Thus, several inventory flow situations exist for different values of the ratios “x” and “y”, for which two examples are depicted in the following Figure 1 and Figure 2.
The stock is positive during the replenishment period, which increases by “(P − D − d)” units of non-defective items over time through the line OA of Figure 2 during the time period t1. Replenishment ends at point A with a maximum number of Q d (non-defective) units. Therefore,
t 1 =   Q d P D d
Production stops at the end of the time period t1 at point K, and the inventory decreases at another demand rate of xD units per unit time, expressed in Figure 2 as line AC during t2. At the end of time period t2, the inventory reaches point C with zero on-hand stock. In this case, if the demand rate remains the same as “D”, then the inventory decreases as indicated by the dashed line AB. Therefore,
t 2 = Q d x D
Therefore,
t 1 + t 2 = { P ( 1 x ) D d } Q d x D ( P D d )
At the end of time t2, thestock is 0 and then decreases at the demand rate yD units per unit time. During time t3, shortages reach a total of Q s units.
t 3 = Q s y D
As soon as the shortage level reaches Qs at point E, production starts. The shortage amount Qs is served at the rate of (P − D − d) units with the current demand, which is fulfilled at the rate of D during timet4. After this period, the stock becomes zero again at the point F. Therefore,
t 4 =   Q s P D d
Hence,
t 3 + t 4 = { P ( 1 y ) D d } Q s y D ( P D d )
A total of Q units is produced in t 1 + t 4 time where W units are defective. Therefore,
t 1 + t 4 = W d = Q P
The following relation is obtained using Equation (7) in addition to Equations (1) and (5).
Q d +   Q s = ( P D d ) W d
If there are no defective items (i.e., d = 0) and no demand (i.e., D = 0), then the inventory is replenished through the line OH.
Given that the ratios of demand rates x and y modify the total demand D during production off-time with positive and negative stock respectively. The average demand rate in a cycle becomes ( 2 + x + y 4 m D )   and requires a total of (Q − W) units of fresh items for the cycle time T. Therefore,
T = 4 m ( Q W ) ( 2 + x + y ) D   =   4 m ( P d ) W ( 2 + x + y ) d D
Thus
t 1 + t 2 T = ( 2 + x + y ) d { P ( 1 x ) D d } Q d 4 m x ( P D d ) ( P d ) W
and
t 3 + t 4 T = ( 2 + x + y ) d { P ( 1 y ) D d } Q s 4 m y ( P D d ) ( P d ) W
m” is an arbitrary constant that has a single value for every set of system parameter values. The theoretical value “m” is explained in Section 4.3.
The foremost target of this work is to minimize the probable total cost of inventory where defective items are formed with perfect goods, but defective items are recyclable, and recycled raw materials are reusable for the production of the same category of new product. Firstly, the manufacturer must define all costs, production features, all abilities of the production procedure, and their position on recycling. The inventory holding cost of the inventory depends on average inventory. Average inventory also increases if defective items are stored properly during production run-time before they are sent to recycling. Recycling processes also incur a cost per defective item.
The following two cases are discussed regarding recycling of defective items.
Case 1:
EPQ model for defective items with three levels of piecewise constant demand under shortages
Case 2:
ERQ model for defective items with three levels of piecewise constant demand under shortages

3.1. Case 1 (EPQ Model for Defective Items with Three Levels of Piecewise Constant Demand under Shortages)

The total cost function of the inventory has the following components.
(1)
Average setup cost or fixed cost (FC):
FC   =   1 T C o = ( 2 + x + y ) d D 4 m ( P d ) W × C o
(2)
Average production cost (PC):
PC   =   1 T Q C p = ( 2 + x + y ) P D 4 m ( P d ) C P
(3)
Average raw material cost (RMC):
Each production lot size is a quantity for which raw materials are purchased for exact production in each time cycle. Therefore,
RMC   =   1 T Q C R = ( 2 + x + y ) P D 4 m ( P d ) C R
(4)
Average holding cost (HC):
Defective items have no holding cost in this case. Therefore,
Total inventory = 1 2 Q d ( t 1 + t 2 )
And average inventory = 1 2 Q d ( t 1 + t 2 T ) = d ( 2 + x + y ) { P ( 1 x ) D d } 4 m x ( P D d ) ( P d ) · Q d 2 2 W
Therefore,
HC   =   d ( 2 + x + y ) { P ( 1 x ) D d } 8 x m ( P d ) ( P D d ) W { ( P D d ) W d   Q s } 2 C h
(5)
Average shortage cost (SC):
SC   =   Q s 2 × C s × t 3 + t 4 T   =   ( 2 + x + y ) d { P ( 1 y ) D d } 4 m y ( P D d ) ( P d ) ( Q s 2 2 W ) × C s
Hereafter, the total cost function is
TCD = FC + PC + RMC + HC + SC
T C D = ( 2 + x + y ) 4 m ( P d ) [ d D W C o + P D C P + P D C R + d { P ( 1 y ) D d } 2 y ( P D d ) W Q s 2 C s +   d { P ( 1 x ) D d } 2 x ( P D d ) W { ( P D d ) W d   Q s } 2 C h ]
The objective is to discover the optimal number of defective items and the shortage level in order to minimize total cost of the inventory per unit time. The corresponding optimal production lot size (EPQ) is also consequently determined.

3.2. Case 2 (ERQ Model for Defective Items with Three Levels of Piecewise Constant Demand under Shortages)

The producer receives a fixed portion of recyclable defective items during regular production, and defective items are sent to be recycled during production off-time. Recyclable materials become raw materials again for the new production cycle, yet to meet the demand, supplementary raw materials must be purchased from the suppliers. The target of the study is to control raw material procurement and the recycle size and production procedure to fulfill demand while minimizing the total inventory cost. The recycling system is shown in Figure 3.
In this case, the production cycle starts with a shortage of amount Q s in the inventory, which accumulates at the rate ( P D d ) over the time t 4 . During this time, the current demand is also fulfilled at the demand rate D . The inventory accumulates at the rate ( P D d ) over time t 1 , and production stops at the end of t 1 . After that, the inventory level starts to decrease due to the demand rate xD over the time t 2 and due to the third level of demand rate yD, while shortages reach Q s overtime t 3 . In the production-runtime ( t 1 + t 4 ) , W units of defective items are produced and are recycled during production off-time ( t 2 + t 3 ) and added to the raw material of the next production cycle. Consequently, the inventory of the ERQ model procures raw materials for ( Q W ) units of items in each cycle time, and the process repeats.
The total cost function of the inventory system comprises the following cost components:
(1)
Average production cost (PC) is:
PC   =   1 T Q C p = ( 2 + x + y ) P D 4 m ( P d ) C P
(2)
Average raw material cost (RMC): Each production lot size is Q , but raw materials for ( Q W ) units of items in each time cycle are purchased. Therefore,
RMC   =   1 T ( Q W ) C R = ( 2 + x + y ) D 4 m C R
(3)
Average holding cost (HC):
Holding cost of defective items is included during time ( t 1 + t 4 ) before the items are sent for recycling. Therefore,
Total inventory
= 1 2 Q d × t 1 + 1 2 Q d × t 2 + 1 2 d ( t 1 + t 4 ) × ( t 1 + t 4 )
= 1 2 Q d ( t 1 + t 2 ) + 1 2 d ( t 1 + t 4 ) 2
Average inventory
= 1 2 Q d ( t 1 + t 2 T ) + 1 2 d ( t 1 + t 4 ) 2 T
= ( 2 + x + y ) 4 m ( P d ) · [ d { P ( 1 x ) D d } 2 x ( P D d ) W { ( P D d ) W d Q s } 2 + D W 2 ]
Hence,
HC   =   ( 2 + x + y ) 4 m ( P d ) [ d { P ( 1 x ) D d } 2 x ( P D d ) W { ( P D d ) W d Q s } 2 + D W 2 ] C h
(4)
Average shortage cost (SC):
SC   =   Q s 2 × C s × t 3 + t 4 T   =   ( 2 + x + y ) d { P ( 1 y ) D d } 4 m y ( P D d ) ( P d ) Q s 2 2 W × C s
(5)
Average recycle cost (RC):
RC   =   1 T W C r = ( 2 + x + y ) d D 4 m ( P d ) C r
Hence, the total cost function is
TCR = FC + PC + RMC + HC + SC + RC
TCR   =   ( 2 + x + y ) 4 m ( P d ) [ d D W C o + P D C P + ( P d ) D C R   + d { P ( 1 x ) D d } 2 x ( P D d ) W { ( P D d ) W d Q s } 2 C h + D W 2 C h + d { P ( 1 y ) D d } y ( P D d ) Q s 2 2 W × C s +   d D C r ]
The objective is to find out the optimal number of defective items (ERQ) to be recycled and the optimal shortage level so as to minimize the total cost of the inventory per unit time. The corresponding optimal production lot size (EPQ) is also consequently determined.
In the next section, the uniqueness of optimal solution of the total cost functions TCD and TCR are explained and justified.

4. Theoretical Derivations

Convexity of total cost functions TCD and TCR relating to the decision variables was established using the Hessian matrix, allowing for the unique optimal solution to be derived. Special cases of the optimal values are also deliberated.

4.1. Case 1 (EPQ Model for Defective Items with Three Levels of Piecewise Constant Demand under Shortages)

The convexity of TCD and uniqueness of optimal solution regarding the decision variables are established to prove the following theorem.
Theorem 1.
The total cost function TCD gives minimum value regarding W and Q s simultaneously when ( 1 x ) D + d > 0 , P ( 1 y ) D d > 0 , and hence TCD emits a unique optimal solution W* and Q s * .
Proof. 
TCD is a function of W and Qs, and therefore the first order partial derivatives of TCD with regard to W and Q s are
( T C D ) W = ( 2 + x + y ) 4 m ( P d ) [ d D C o W 2   ( P D d ) { P ( 1 x ) D d } C h 2 x d +   d { P ( 1 x ) D d } C h Q s 2 2 x ( P D d ) W 2 d { P ( 1 y ) D d } C s Q s 2 2 y ( P D d ) W 2 ]
And
( T C D ) Q s = ( 2 + x + y ) 4 m ( P d ) [   d { P ( 1 x ) D d } x ( P D d ) W { ( P D d ) W d   Q s } C h + d { P ( 1 y ) D d } C s y ( P D d ) W Q s ]
The second order partial derivatives of TCD are
2 ( T C D ) W 2 = d D ( 2 + x + y ) 2 m ( P d ) W 3 C o + d ( 2 + x + y ) [ y { P ( 1 x ) D d } C h + x { P ( 1 y ) D d } C s ] Q s 2 4 m x y ( P d ) ( P D d ) W 3
2 ( T C D ) Q s 2 = d ( 2 + x + y ) [ y { P ( 1 x ) D d } C h + x { P ( 1 y ) D d } C s ] 4 m x y ( P d ) ( P D d ) W
Also
2 ( T C D ) W Q s = d ( 2 + x + y ) [ y { P ( 1 x ) D d } C h + x { P ( 1 y ) D d } C s ] Q s 4 m x y ( P d ) ( P D d ) W 2
And
2 ( T C D ) Q s W = d ( 2 + x + y ) [ y { P ( 1 x ) D d } C h + x { P ( 1 y ) D d } C s ] Q s 4 m x y ( P d ) ( P D d ) W 2
The Hessian matrix of TCD is
H i j = [ 2 ( T C D ) W 2 2 ( T C D ) W Q s 2 ( T C D ) Q s W 2 ( T C D ) Q s 2 ]
Hence, the first principal minor
| H 11 | = 2 ( T C D ) W 2   =   ( 2 + x + y ) 4 m ( P d ) [ 2 d D W 3 × C o +   d { P ( 1 x ) D d } C h x ( P D d ) W 3 Q s 2 + d { P ( 1 y ) D d } C s y ( P D d ) W 3 Q s 2 ]
It is assumed that, P D d > 0 . Therefore,
P ( 1 x ) D d > 0 ,   P ( 1 y ) D d > 0 and   P d > 0 .   Consequently ,   | H 11 | > 0 .
Once more the second principal minor
| H 22 |   =   2 ( T C D ) W 2 × 2 ( T C D ) Q s 2 2 ( T C D ) Q s W × 2 ( T C D ) W Q s =   d 2 D ( 2 + x + y ) 2 [ y { P ( 1 x ) D d } C h + x { P ( 1 y ) D d } C s ] C o 8 m 2 x y ( P d ) 2 ( P D d ) W 4 >   0 .
Since the first and second principal minors of the Hessian matrix for TCD are positive, the Hessian matrix is positive definite. Hence TCD is a nonnegative, differentiable, and strictly convex function regarding the decision variables W and Qs simultaneously, and TCD is minimum at the unique optimal values W * and Q s * . Solving the necessary conditions ( T C D ) W = 0 and ( T C D ) Q s = 0 from Equations (14) and (15) the following results are obtained.
Q s W =   y ( P D d ) { P ( 1 x ) D d } C h d [ x { P ( 1 y ) D d } C s + y { P ( 1 x ) D d } C h ]
And
W =   d 2 D C O C h x { P ( 1 y ) D d } C s + y { P ( 1 x ) D d } C h ( P D d ) { P ( 1 x ) D d } { P ( 1 y ) D d } C s
Therefore, the optimal solution is given by
W *   =   d 2 D C O C h x { P ( 1 y ) D d } C s + y { P ( 1 x ) D d } C h ( P D d ) { P ( 1 x ) D d } { P ( 1 y ) D d } C s
And
Q s * = y 2 D C O ( P D d ) { P ( 1 x ) D d } { P ( 1 y ) D d } C h [ x { P ( 1 y ) D d } C s + y { P ( 1 x ) D d } C h ] C s
Then the corresponding optimal production lot size is
Q *   =   P 2 D C O 1 ( P D d ) { P ( 1 x ) D d } { P ( 1 y ) D d } x { P ( 1 y ) D d } C s + y { P ( 1 x ) D d } C h C h C s
And the optimal cycle time is
T * =   2 C O D . { P ( 1 x ) D d } { P ( 1 y ) D d } ( P D d ) [ x { P ( 1 y ) D d } C s + y { P ( 1 x ) D d } C h ] · C s + C h C h C s
Also the maximum on-hand stock is
Q d * = x 2 D C O ( P D d ) { P ( 1 y ) D d } { P ( 1 x ) D d } C s [ x { P ( 1 y ) D d } C s + y { P ( 1 x ) D d } C h ] C h
 □

Special Cases

(i)
If d = 0, then W*= 0 .
Therefore, no defective items will be produced in the inventory if the defective rate is zero.
(ii)
If x = y = 1 and d = 0, then we get
Q *   = 2 D C O C h P ( P D ) · C s + C h C s
This is the lot size in the standard EPQ model. Therefore, the optimal solution of the total cost function of this model conforms to the standard EPQ model.
(iii)
If x = y, then,
W *   =   d 2 D C O C h x ( C s + C h ) ( P D d ) { P ( 1 x ) D d } C s
and Q* = P 2 D C O C h x ( C s + C h ) ( P D d ) { P ( 1 x ) D d } C s
These are the results of the inventory model for two levels of piecewise constant demand.
(iv)
If x = y = 1, then,
W *   =   d 2 D C O C h ( C s + C h ) ( P D d ) ( P d ) C s
And Q* = P 2 D C O C h ( C s + C h ) ( P D d ) ( P d ) C s
These are the results of the inventory model of defective items for constant demand.

4.2. Case 2 (ERQ Model for Defective Items with Three Levels of Piecewise Constant Demand under Shortages)

The convexity of TCR and hence the uniqueness of optimal values of the decision variables are determined to prove the following theorem.
Theorem 2.
The total cost function TCR gives a minimum value regarding the decision variables W and Q s and hence a unique solution exists for W* and Q s * .
Proof. 
TCR is a function of W and Qs, and the relation between TCR and TCD is obtained as
TCR = TCD ( 2 + x + y ) D 4 m ( P d ) [ d ( C R C r ) W 2 C h ]
From the first order and second order partial derivatives, the following results are obtained from Equation (27).
2 ( T C R ) W 2 = 2 ( T C D ) W 2 ;   2 ( T C R ) Q s 2 = 2 ( T C D ) Q s 2 ; 2 ( T C R ) W Q s = 2 ( T C D ) W Q s   and   2 ( T C R ) Q s W = 2 ( T C D ) Q s W
The Hessian matrix of TCR is
H i j = [ 2 ( T C R ) W 2 2 ( T C R ) W Q s 2 ( T C R ) Q s W 2 ( T C R ) Q s 2 ] = [ 2 ( T C D ) W 2 2 ( T C D ) W Q s 2 ( T C D ) Q s W 2 ( T C D ) Q s 2 ]
Therefore, from Theorem 1, it can be concluded that the first principal minor
| H 11 | = 2 ( T C R ) W 2 = 2 ( T C D ) W 2 > 0
And the second principal minor
| H 22 |   =   2 ( T C R ) W 2 × 2 ( T C R ) Q s 2 2 ( T C R ) Q s W × 2 ( T C R ) W Q s   =   2 ( T C D ) W 2 × 2 ( T C D ) Q s 2 2 ( T C D ) Q s W × 2 ( T C D ) W Q s >   0 .
Since the first and second principal minor of Hessian matrix for TCR are positive, the Hessian matrix is positive definite. Hence TCR is a nonnegative, differentiable, and strictly convex function regarding W and Qs concurrently and then TCR is minimum at the unique optimal solution W * and Q s * .
The necessary conditions of minimization of TCR are the first order partial derivatives of TCR with respect to W and Q s are equivalent to zero, and hence the following results are obtained.
Q s W = y ( P D d ) { P ( 1 x ) D d } C h d [ x { P ( 1 y ) D d } C s + y { P ( 1 x ) D d } C h ]
And
W =   d 2 D C O C h x { P ( 1 y ) D d } C s + y { P ( 1 x ) D d } C h ( P D d ) { P ( 1 x ) D d } { P ( 1 y ) D d } C s + d D [ x { P ( 1 y ) D d } C s + y { P ( 1 x ) D d } C h ]
Therefore, the required economic recycle quantity (ERQ) is
W *   =   d 2 D C O C h x { P ( 1 y ) D d } C s + y { P ( 1 x ) D d } C h ( P D d ) { P ( 1 x ) D d } { P ( 1 y ) D d } C s + d D [ x { P ( 1 y ) D d } C s + y { P ( 1 x ) D d } C h ]
And shortage level is
Q s * =   y ( P D d ) { P ( 1 x ) D d } C h [ x { P ( 1 y ) D d } C s + y { P ( 1 x ) D d } C h ] C h 2 D C O ( P D d ) { P ( 1 x ) D d } { P ( 1 y ) D d } C s + d D [ x { P ( 1 y ) D d } C s + y { P ( 1 x ) D d } C h ]
The corresponding production lot size of the inventory is
Q *   =   P 2 D C O C h x { P ( 1 y ) D d } C s + y { P ( 1 x ) D d } C h ( P D d ) { P ( 1 x ) D d } { P ( 1 y ) D d } C s + d D [ x { P ( 1 y ) D d } C s + y { P ( 1 x ) D d } C h ]
And the optimal cycle time is
T *   =   2 C O D . { P ( 1 x ) D d } { P ( 1 y ) D d } ( P D d ) { P ( 1 x ) D d } { P ( 1 y ) D d } C s + d D [ x { P ( 1 y ) D d } C s + y { P ( 1 x ) D d } C h ] × C h + C s [ x { P ( 1 y ) D d } C s + y { P ( 1 x ) D d } C h ] C h
Also, the maximum on-hand inventory is
Q d * = x ( P D d ) { P ( 1 y ) D d } C s [ x { P ( 1 y ) D d } C s + y { P ( 1 x ) D d } C h ] C h 2 D C O ( P D d ) { P ( 1 x ) D d } { P ( 1 y ) D d } C s + d D [ x { P ( 1 y ) D d } C s + y { P ( 1 x ) D d } C h ]
 □

Special Cases

(i)
If d = 0, then W* = 0 .
Therefore, no defective items will be produced in the inventory if the defective rate is zero.
(ii)
If x = y = 1 and d = 0, then,
Q *   = 2 D C O C h P ( P D ) C s + C h C s
This is the lot size in the standard EPQ model.
Therefore, the optimal outcomes are conformable with the standard EPQ model.
(iii)
If x = y then,
W *   =   d 2 D C O C h x ( C s + C h ) ( P D d ) { P ( 1 x ) D d } C s + x d D ( C s + C h )
Q *   =   P 2 D C O C h x ( C s + C h ) ( P D d ) { P ( 1 x ) D d } C s + x d D ( C s + C h )
These are the results of the ERQ model for two levels of piecewise constant demand.
(iv)
If x = y = 1, then,
W *   =   d 2 D C O C h ( C s + C h ) ( P D d ) ( P d ) C s + d D ( C s + C h )
Q *   =   P 2 D C O C h ( C s + C h ) ( P D d ) ( P d ) C s + d D ( C s + C h )
These are the results of the ERQ model of defective items for constant demand.

4.3. Value of “m”

It is assumed that the arbitrary constant “m” determines the accurate average demand as well as total inventory cost but has no significance on the optimal solution. It is introduced to avoid the complexity of the theoretical results. The assumption is justified by the obtained optimal solutions of TCD and TCR and their special cases. The value of “m” is calculated in the following manner.
Addition of Equations (3) and (6) gives the cycle time
T =   x { P ( 1 y ) D d } Q s +   y { P ( 1 x ) D d } Q d x y D ( P D d )
Then, Equation (9) implies
m = d ( 2 + x + y ) [ x { P ( 1 y ) D d } Q s +   y { P ( 1 x ) D d } Q d ] 4 x y ( P d ) ( P D d ) W
Putting the optimal values of Q s and Q d of both of the total cost function TCD and TCR, the same value of “m” is obtained as follows:
m =   ( 2 + x + y ) { P ( 1 x ) D d } { P ( 1 y ) D d } ( C s + C h ) 4 ( P d ) [ x { P ( 1 y ) D d } C s + y { P ( 1 x ) D d } C h ]

5. Numerical Illustrations and Case Study

The inventory model is analyzed considering a piecewise constant demand rate that may vary from one level to another. Certain products, for instance, such as bricks, mustard oil, and shopping mall products, have different demands during their production-run periods than during production-off shortage periods. The practical implications of the proposed model can be observed in a brick production company. For the production time of brick raw materials, some imperfect bricks are made before firing, for instance, during the setting and drying periods. These defective bricks are collected and stored properly in stock before they are sent for recycling during production off-time and added to the inventory for the next production cycle. Furthermore, the demand for bricks during the production off-time and shortage time may increase or decrease from the demand during their production run-time.
In order to study the applicability of the model, numerical illustrations are provided below considering the following static input values of parameters taken from MJB brick production company in Bangladesh where: P = 5000, D = 4500, d = 100, Co = 1000, Cp = 50, Ch= 10, CR = 50, Cr = 5, Cs = 3, x = 0.75, and y = 0.5.

5.1. For Case 1

Optimal decision variable values:
W* = 136, Q* = 6822, Qs* = 414, Qd* = 131, and T * = 1.5877
Optimal cycle time level values:
t 1 * = 0.3275, t 2 * = 0.0388, t 3 * = 0.1843, t 4 * = 1.0370
Optimal costs and the total inventory cost:
FC = 630, PC = 214,859, SC = 478, RMC = 214,859, HC = 151, and TCD = 430,978
Cycle time verification for TCD:
t 1 * +   t 2 * +   t 3 * +   t 4 * =   0.3275 + 0.0388+ 0.1843+ 1.0370 = 1.5876 = T *
Hence, the model is consistent.
For the values of independent variables in the range 120 < W < 150 and 400 < Qs < 430, Figure 4 depicts the convexity of the total cost function TCD with respect to W and Qs.
For Qs = 414 and 120 < W < 150, Figure 5 depicts the convexity of the total cost function TCD vs. W.
Fixing W = 136 and 250 < Qs < 350, the following Figure 6 shows that TCD is convex with respect to Qs, and the total inventory cost reaches minimum at “Qs = 414”.

5.2. For Case 2

Optimal decision variable values: W* = 98, Q* = 4910, Qs* = 298, Qd* = 94, T * = 1.1426
Optimal cycle time level values: t 1 * = 0.2357, t 2 * = 0.0279, t 3 * = 0.1326, t 4 * = 0.7462
Optimal costs and total inventory cost: FC = 875, PC = 214,859, SC = 344, RC = 430, RMC = 210,562, HC = 530 and TCR = 427,602
Cycle time verification: t 1 * +   t 2 * +   t 3 * +   t 4 * =   0.2357 + 0.0279 + 0.1326 + 0.7462 = 1.1425 = T * . Hence, the model is consistent.
For the values of independent variables in the range 80 < W < 120 and 250 < Qs < 350, the following Figure 7 depicts the convexity of the total cost function TCR for W and Qs.
For Qs = 300 and 80 < W < 120, the following figure (Figure 8) depicts the convexity of the total cost function with respect to “W”.
For W = 98 and 250 < Qs < 350, Figure 9 shows that TCR is convex with respect to Qs and the total cost of the inventory is minimum at “Qs = 300”.

6. Sensitivity Analysis

In this section, a sensitivity analysis of the anticipated model of various system parameters is examined. When the sensitivity of a parameter is analyzed, only changes to a parameter are considered keeping other parameters as fixed. The sensitivity analyses of the key parameters are given in Table 2, Table 3, Table 4, Table 5, Table 6, Table 7, Table 8 and Table 9.
From Table 2, it is observed that total cost increases with increasing values of “x”, and as a result, total order quantity upsurges. In contrast, the total cycle length “T” decreases. The significant increases in all the associated costs, e.g., PC, HC, RC are also noted. However, the maximum shortage Qs is notably reduced.
Explanations: From Table 3, it is clear that with the rise of “y”, the total cost is amplified, while the total order quantity is also increased. In contrast, the total cycle length “T” declines. A noteworthy intensification in all the related costs, e.g., PC, HC, RC, can also be observed. Moreover, the maximum shortage Qs notably rises.

Cost–Benefitof ERQ Model

Cost–benefit is a tool used to gain insights into the decision-making process. In the cost–benefit analysis, we considered the comparative total cost of two inventory systems: the system without recycling and system with recycling of defective items. Thus,
Average cost–benefit ACB is
ACB   =   TCD     TCR   =   ( 2 + x + y ) D 4 m ( P d ) [ d ( C R   C r ) W 2 C h ]
Percentage of cost–benefit (PCB) is
PCB   =   ( A C B T C D × 100 ) %
Therefore, from Equation (38), we see that the cost–benefit increases when the cost of raw materials increases and the cost–benefit decreases when holding cost, recycling cost, or both costs increase.
While the effect of a parameter is analyzed, only changes to one parameter is considered while other parameters remain static.

7. Managerial Insights

The proposed model helps to control inventories and to determine the economic recycle quantity (ERQ), optimal production lot size, cycle time, total inventory cost, maximum on-hand stock, and maximum backlogs. From the numerical analysis, sensitivity analysis (Table 2, Table 3, Table 4, Table 5, Table 6, Table 7, Table 8 and Table 9), and cost–benefit analysis, the following findings provide managerial insights to inventory managers.
Recycling larger or smaller quantities of defective items than the ERQ will increase the total inventory cost. Therefore, inventory managers must evaluate the ERQ of the item before recycling.
When the company decides to recycle, production lot sizes should be set up in favor of the ERQ, and this strategy will provide greater facility to the manager to decline the total cost.
If raw material costs increase significantly, recycling will benefit the company because it will diminish the load on raw-materials and use the wastage items.
The total cost is more sensitive to demand increases during shortage time than it is to increases during production time with positive stock.
If the defective rate is high, companies will receive better cost–benefit from recycling because it produces more defective items, which are usually treated as rejected, but the recycling process converts them to useable, which provides huge revenue for the managers.

8. Conclusions

In this study, we developed an imperfect production model with and without the recycling of defective items. The model is a generalization of models considering recycled raw materials. It is worth mentioning that recycled materials can be an alternative source of raw materials in the event of any uncertainties. Moreover, recycling defective production reduces waste and hence is also environmentally friendly. The study can be applied to bricks, textile, glass, paper, and jute production. It is demonstrated that the total cost function is strictly convex with a unique optimal solution. While this model considers a single item production system, multi-item inventory systems of the same recyclable raw material can be developed for further research. This paper outlines when to stop the production process and when to start it. Proper management of the recycling process is also illustrated in the model with some insights for managers.
Furthermore, the model can be extended to consider time-dependent defective rates and selling prices (Mashud et al. 2020a). Considering the environmental emissions (Mashud et al. 2021b; Mishra et al. 2021) would be another interesting extension of this study. As the recycling process is shown in the model, a proper waste management (Tsai et al. 2021) could be another important discussion that is missing in this study. Some marketing strategies, for instance, discount on imperfect items (Mashud et al. 2020a), trade credit financing (Liao et al. 2018, 2020; Srivastava et al. 2018; Mashud et al. 2021a, 2021d) can be implemented in the proposed model to make it more lucrative to the practitioners. Other imperfect process in both the EOQ and EPQ models (Lin and Srivastava 2015; Srivastava et al. 2021) can also be consider as an interesting extension.

Author Contributions

Conceptualization, A.A., M.M.M., M.S.U., A.H.M.M. and H.-M.W.; methodology, A.A., M.M.M., M.S.U., A.H.M.M., H.-M.W., S.S.S. and H.M.S.; software, A.A., M.M.M., M.S.U. and A.H.M.M.; validation, A.A., M.M.M., M.S.U., A.H.M.M. and H.-M.W.; formal analysis, A.A., M.M.M., M.S.U., A.H.M.M. and H.-M.W.; investigation, A.A., M.M.M., M.S.U., A.H.M.M. and H.-M.W.; data curation, A.A., M.M.M., M.S.U., A.H.M.M. and H.-M.W.; writing—original draft preparation, A.A., M.M.M., M.S.U., A.H.M.M. and H.-M.W.; writing—review and editing, S.S.S. and H.M.S.; supervision, S.S.S. and H.M.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was completed at the authors’ own cost. Funding from a third party was not used to carry out this research work.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data used to justify the proposed model are given in the manuscript.

Conflicts of Interest

We do hereby declare that we do not have any conflicts of interest with other works.

References

  1. Alharkan, Ibrahim, Mustafa Saleh, Mageed Ghaleb, Abdulsalam Farhan, and Ahmed Badwelan. 2020. Simulation-Based Optimization of a Two-Echelon Continuous Review Inventory Model with Lot Size-Dependent Lead Time. Processes 8: 1014. [Google Scholar] [CrossRef]
  2. Al-Salamah, Muhammad. 2019. Economic Production Quantity in an Imperfect Manufacturing Process with Synchronous and Asynchronous Flexible Rework Rates. Operations Research Perspectives 6: 100103. [Google Scholar] [CrossRef]
  3. Bai, S. X., and S. Varanasi. 1996. An optimal production flow control problem with piecewise constant demand. Mathematical and Computer Modelling 24: 87–107. [Google Scholar] [CrossRef]
  4. Benkherouf, Lakdere, and Mohamed Omar. 2017. Optimal manufacturing batch size with rework for a finite-horizon and time-varying demand rates inventory model. RAIRO-Operations Research 51: 173–87. [Google Scholar] [CrossRef]
  5. Chiu, S. W., C. T. Tseng, M. F. Wu, and P. C. Sung. 2014. Multi-item EPQ model with scrap, rework and multi-delivery using common cycle policy. Journal of Applied ResearchandTechnology 12: 615–22. [Google Scholar] [CrossRef]
  6. Daryanto, Yosef, and Hui Ming Wee. 2019. Low carbon economic production quantity model for imperfect quality deteriorating items. Journal of Industrial Engineering and Management 1: 1–8. [Google Scholar] [CrossRef]
  7. Harris, Ford Whitman. 1915. What quantity to make at once. The Library of Factory Management 5: 47–52. [Google Scholar]
  8. Hasan, Rakibul, Abu Hashan Md Mashud, Yosef Daryanto, and Hui Ming Wee. 2020. A non-instantaneous inventory model of agricultural products considering deteriorating impacts and pricing policies. Kybernetes. ahead-of-print. [Google Scholar] [CrossRef]
  9. Hsu, Jia Tzer, and Lie Fern Hsu. 2014. A supplement to an EOQ model with imperfect quality items, inspection errors, shortage backordering, and sales returns. International Journal of Industrial Engineering Computation 5: 199–210. [Google Scholar] [CrossRef] [Green Version]
  10. Hsu, Lie Fern, and Jia Tzer Hsu. 2016. Economic production quantity (EPQ) models under an imperfect production process with shortages backordered. International Journal of Systems Science 47: 852–67. [Google Scholar]
  11. Kang, Chang Wook, Misbah Ullah, and Biswajit Sarkar. 2018. Optimum ordering policy for an imperfect single-stage manufacturing system with safety stock and planned backorder. The International Journal of Advanced Manufacturing Technology 95: 109–20. [Google Scholar] [CrossRef]
  12. Liao, Jui Jung, Kuo Nan Huang, Kun Jen Chung, Shy Der Lin, Sheng Tu Chuang, and Hari Mohan Srivastava. 2020. Optimal ordering policy in an economic order quantity (EOQ) model for non-instantaneous deteriorating items with defective quality and permissible delay in payments. Revista de la Real Academia de Ciencias Exactas, Físicas y Naturales. Serie A Matemáticas 114: 1–26. [Google Scholar] [CrossRef]
  13. Liao, Jui Jung, Kuo Nan Huang, Kun Jen Chung, Shy Der Lin, Pin Shou Ting, and Hari Mohan Srivastava. 2018. Retailer’s optimal ordering policy in the EOQ model with imperfect-quality items under limited storage capacity and permissible delay. Mathematical Methods in the Applied Sciences 41: 7624–40. [Google Scholar] [CrossRef]
  14. Lin, Tien Yu, and Hari Mohan Srivastava. 2015. A two-warehouse inventory model with quantity discounts and maintenance actions under imperfect production processes. Applied Mathematics and Information Sciences 9: 2493–505. [Google Scholar]
  15. Manna, Amalesh Kumar, Joyanta Kumar Dey, and Shyamal Kumar Mondal. 2017. Imperfect production inventory model with production rate dependent defective rate and advertisement dependent demand. Computers & Industrial Engineering 104: 9–22. [Google Scholar]
  16. Mashud, Abu Hashan Md, Dipa Roy, Yosef Daryanto, and Hui Ming Wee. 2021a. Joint pricing deteriorating inventory model considering the product life cycle and advance payment. RAIRO-Operations Research 55: S1069–S1088. [Google Scholar] [CrossRef]
  17. Mashud, Abu Hashan Md, Dipa Roy, Yosef Daryanto, and Mohd Helmi Ali. 2020a. A Sustainable Inventory Model with Imperfect Products, Deterioration, and Controllable Emissions. Mathematics 8: 2049. [Google Scholar] [CrossRef]
  18. Mashud, Abu Hashan Md, Dipa Roy, Yosef Daryanto, Ripon Kumar Chakraborty, and Ming Lang Tseng. 2021b. A controllable carbon emission and deterioration inventory model with advance payments scheme. Journal of Cleaner Production 296: 126608. [Google Scholar] [CrossRef]
  19. Mashud, Abu Hashan Md, Hui Ming Wee, and Chiao Ven Huang. 2021c. Preservation technology investment, trade-credit and partial backordering model for a non-instantaneous deteriorating inventory. RAIRO-Operations Research 55: 551–77. [Google Scholar] [CrossRef] [Green Version]
  20. Mashud, Abu Hashan Md, Hui Ming Wee, Chiao Ven Huang, and J. Z. Wu. 2020b. Optimal Replenishment Policy for Deteriorating Products in a Newsboy Problem with Multiple Just-in-Time Deliveries. Mathematics 8: 1981. [Google Scholar] [CrossRef]
  21. Mashud, Abu Hashan Md, Md Rakibul Hasan, Yosef Daryanto, and Hui Ming Wee. 2021d. A resilient hybrid payment supply chain inventory model for post Covid-19 recovery. Computers & Industrial Engineering, 107249. [Google Scholar] [CrossRef]
  22. Mishra, Umakanta, Abu Hashan Md Mashud, Ming Lang Tseng, and Jei Zang Wu. 2021. Optimizing a Sustainable Supply Chain Inventory Model for Controllable Deterioration and Emission Rates in a Greenhouse Farm. Mathematics 9: 495. [Google Scholar] [CrossRef]
  23. Mokhtari, Hadi. 2018. A joint internal production and external supplier order lot size optimization under defective manufacturing and rework. The International Journal of Advanced Manufacturing Technology 95: 1039–58. [Google Scholar]
  24. Nobil, Amir Hossein, Amir Hossein Afsar Sedigh, and Leopoldo Eduardo Cardenas Barron. 2020. Reorder point for the EOQ inventory model with imperfect quality items. Ain Shams Engineering Journal 11: 1339–43. [Google Scholar] [CrossRef]
  25. Oh, Yong Hui, and Hark Hwang. 2006. Deterministic inventory model for recycling system. Journal of Intelligent Manufacturing 17: 423–28. [Google Scholar] [CrossRef]
  26. Onwude, Daniel, Guangnan Chen, Nnanna Eke-Emezie, Abraham Kabutey, Alfadhl Yahya Khaled, and Barbara Sturm. 2020. Recent Advances in Reducing Food Losses in the Supply Chain of Fresh Agricultural Produce. Processes 8: 1431. [Google Scholar] [CrossRef]
  27. Palanivel, M., and Ramasamy Uthayakumar. 2016. An inventory model with imperfect items, stock dependent demand and permissible delay in payments under inflation. RAIRO-Operations Research 50: 473–89. [Google Scholar] [CrossRef]
  28. Popescu, Cristina Raluca Gh, and Georghe N. Popescu. 2019. An Exploratory Study Based on a Questionnaire Concerning Green and Sustainable Finance, Corporate Social Responsibility, and Performance: Evidence from the Romanian Business Environment. Journal of Risk and Financial Management 12: 162. [Google Scholar] [CrossRef] [Green Version]
  29. Rosenblatt, M. J., and H. L. Lee. 1986. Economic production cycles with imperfect production processes. IIE Transaction 18: 48–55. [Google Scholar]
  30. Salameh, M. K., and M. Y. Jaber. 2000. Economic production quantity model for items with imperfect quality. International Journal of Production Economics 64: 59–64. [Google Scholar] [CrossRef]
  31. Srivastava, Hari Mohan, Ghi Feng Yen, An Kuo Lee, Yi Xiu Wu, and Shy Der Lin. 2018. The optimal retailer’s economic production quantity (EPQ) policies with two-level trade credit under alternate due date of payment and limited storage capacity. Applied Mathematics and Information Sciences 12: 1073–89. [Google Scholar] [CrossRef]
  32. Srivastava, Hari Mohan, Kun Jen Chung, Jui Jung Liao, Shy Der Lin, and Shi Fang Lee. 2021. An accurate and reliable mathematical analytic solution procedure for the EOQ model with non-instantaneous receipt under supplier credits. Revista de la Real Academia de Ciencias Exactas, Físicas y Naturales. Serie A Matemáticas 115: 1–22. [Google Scholar] [CrossRef]
  33. Tayyab, Muhammad, Biswajit Sarkar, and Bernardo Nugroho Yahya. 2019a. Imperfect Multi-Stage Lean Manufacturing System with Rework under Fuzzy Demand. Mathematics 7: 13. [Google Scholar] [CrossRef] [Green Version]
  34. Tayyab, Muhammad, Jihed Jemai, Han Lim, and Biswajit Sarkar. 2019b. A sustainable development framework for a cleaner multi-item multi-stage textile production system with a process improvement initiative. Journal of Cleaner Production 246: 119055. [Google Scholar] [CrossRef]
  35. Taft, E. W. 1918. The most economical production lot. The Iron Age 101: 1410–12. [Google Scholar]
  36. Tsai, Feng Ming, Tat Dat Bui, Ming Lang Tseng, Ming K. Lim, Kuo Jui Wu, and Abu Hashan Md Mashud. 2021. Assessing hierarchical sustainable solid waste management structure in qualitative information: Policy and regulations drive social impacts and stakeholder participation. Resources, Conservation and Recycling 168: 105285. [Google Scholar] [CrossRef]
  37. Viji, G., and K. Karthikeyan. 2016. An economic production quantity model for three levels of production with Weibull distribution deterioration and shortage. Ain Shams Engineering Journal 9: 1481–87. [Google Scholar] [CrossRef] [Green Version]
Figure 1. On-hand stock flow for x > 1 and y < x.
Figure 1. On-hand stock flow for x > 1 and y < x.
Jrfm 14 00173 g001
Figure 2. On-hand stock flow for x < 1 and y > x.
Figure 2. On-hand stock flow for x < 1 and y > x.
Jrfm 14 00173 g002
Figure 3. Framework of recycling.
Figure 3. Framework of recycling.
Jrfm 14 00173 g003
Figure 4. Convexity of TCD with respect to W and Qs.
Figure 4. Convexity of TCD with respect to W and Qs.
Jrfm 14 00173 g004
Figure 5. TCD vs. W graph.
Figure 5. TCD vs. W graph.
Jrfm 14 00173 g005
Figure 6. TCD vs. Qs graph.
Figure 6. TCD vs. Qs graph.
Jrfm 14 00173 g006
Figure 7. Convexity of TCR relating to W and Qs.
Figure 7. Convexity of TCR relating to W and Qs.
Jrfm 14 00173 g007
Figure 8. TC vs. W graph.
Figure 8. TC vs. W graph.
Jrfm 14 00173 g008
Figure 9. TCR vs. Qs graph.
Figure 9. TCR vs. Qs graph.
Jrfm 14 00173 g009
Table 1. Notation used in the model.
Table 1. Notation used in the model.
NotationUnitDescription
PConstantProduction rate
DConstantRate of demand of non-defective items during production run-time
dConstantProduction rate defective items
Co$/CycleSetup cost
Ch$/UnitsInventory holding cost
Cp$/UnitsPer unit production cost
Cs$/UnitsShortage cost
CR$/UnitsRaw material cost per unit item
Cr$/UnitsRecycle cost per unit item of defective items
TCD$/Unit timeTotal cost of inventory (for case 1)
TCR$/Unit timeTotal cost of inventory (for case 2)
t1Time unitsProduction run-time with positive stock
t2Time unitsProduction off-time with positive stock
t3Time unitsProduction off-time with negative stock
t4Time unitsProduction run-time with negative stock
xConstantRatio of demand rates of production off-time with positive stock and production run-time
yConstantRatio of demand rates of production off-time with negative stock and production run-time
mConstantShape parameter of the accurate average demand
Decision variables
WUnitsQuantity of defective items produced per cycle time
QUnitsVolume of lot
QsUnitsMaximum shortage level
QdUnitsMaximum on-hand stock of non-defective items
TTime unitsProduction cycle time. Therefore, T = t1 + t2 + t3 + t4
Table 2. Sensitivity of “x” to different decision variables and costs. The static input values of the parameters are: P = 5000, D = 4500, d = 100, Co = 1000, Cp = 50, Ch = 10, CR = 50, Cr = 5, Cs = 3, y = 1.
Table 2. Sensitivity of “x” to different decision variables and costs. The static input values of the parameters are: P = 5000, D = 4500, d = 100, Co = 1000, Cp = 50, Ch = 10, CR = 50, Cr = 5, Cs = 3, y = 1.
x0.50.751.755Observations
W*98989999Increase
Q*4894490949294940Increase
Qs*301299295293Decrease
Qd*909499102Increase
t 1 * 0.2250.2350.2480.254Increase
t 2 * 0.04100.0270.0120.004Decrease
t 3 * 0.1330.1320.1310.130Decrease
t 4 * 0.7520.7460.7370.733Decrease
T*1.1521.1421.1291.122Decrease
FC867875885890Increase
RMC208,019210,562213,820215,565Increase
PC212,264214,859218,184219,965Increase
HC519530545552Increase
SC347344340338Increase
RC424430436440Increase
TCR422,442427,602434,210437,755Increase
Table 3. Sensitivity of “y” to different decision variables and costs. The static input values of the parameters are: P = 5000, D = 4500, d = 100, Co = 1000, Cp = 50, Ch = 10, CR = 50, Cr = 5, Cs = 3, x = 1.
Table 3. Sensitivity of “y” to different decision variables and costs. The static input values of the parameters are: P = 5000, D = 4500, d = 100, Co = 1000, Cp = 50, Ch = 10, CR = 50, Cr = 5, Cs = 3, x = 1.
y0.50.81.5025Observations
W*98.3699.36100.4100.84101.2Increase
Q*49184968502050425058Increase
Qs*297305309313316Increase
Qd*9691.7290.7289.588Decrease
t1*0.240980.22930.226820.223830.22157Decrease
t2*0.021420.025470.020160.019890.01969Decrease
t3*0.132020.0849920.046050.027890.01404Decrease
t4*0.742650.76430.777210.784550.79016Increase
T*1.137091.104061.070251.056191.04548Decrease
FC879905934947956Increase
RMC211,938220,500229,841233,914237,093Increase
PC216,263225,000234,532238,688241,931Increase
HC536552575584591Increase
SC342352359362364Increase
RC432450469477483Increase
TCR430,392447,762466,711474,972481,421Increase
Table 4. Sensitivity of defective rate to cost–benefit. The static input values of the parameters are: P = 5000, D = 4500, Co = 1000, Cp = 50, Ch = 10, CR = 50, Cr = 5, Cs = 3, x = 1.5, y = 1.5.
Table 4. Sensitivity of defective rate to cost–benefit. The static input values of the parameters are: P = 5000, D = 4500, Co = 1000, Cp = 50, Ch = 10, CR = 50, Cr = 5, Cs = 3, x = 1.5, y = 1.5.
d100110120130140
PCB0.797%0.878%0.96%1.042%1.125%
RemarksCost–benefit increases with the increase in “d”.
Table 5. Sensitivity of “recycle cost” to cost–benefit. The static input values of the parameters are: P = 5000, D = 4500, d = 100, Co = 1000, Cp = 50, Ch = 10, CR = 50, Cs = 3, x = 1.5, y = 1.5.
Table 5. Sensitivity of “recycle cost” to cost–benefit. The static input values of the parameters are: P = 5000, D = 4500, d = 100, Co = 1000, Cp = 50, Ch = 10, CR = 50, Cs = 3, x = 1.5, y = 1.5.
Cr510152025
PCB0.797%0.697%0.59%0.497%0.398%
RemarksCost–benefit decreases with increase in recycle cost.
Table 6. Sensitivity of “holding cost” to cost–benefit. The static input values of the parameters are: P = 5000, D = 4500, d = 100, Co = 1000, Cp = 50, CR = 50, Cr = 5, Cs = 3, x = 1.5, y = 1.5.
Table 6. Sensitivity of “holding cost” to cost–benefit. The static input values of the parameters are: P = 5000, D = 4500, d = 100, Co = 1000, Cp = 50, CR = 50, Cr = 5, Cs = 3, x = 1.5, y = 1.5.
Ch1020304050
PCB0.797%0.737%0.689%0.648%0.613%
RemarksCost–benefit decreases with increase in holding cost.
Table 7. Sensitivity of “x” to cost–benefit. The static input values of the parameters are: P = 5000, D = 4500, d = 100, Co = 1000, Cp = 50, Ch = 10, CR = 50, Cr = 5, Cs = 3, y = 1.5.
Table 7. Sensitivity of “x” to cost–benefit. The static input values of the parameters are: P = 5000, D = 4500, d = 100, Co = 1000, Cp = 50, Ch = 10, CR = 50, Cr = 5, Cs = 3, y = 1.5.
x0.511.522.5
PCB0.7976%0.7972%0.7970%0.7969%0.7969%
Remarksx” is not sensitive to cost–benefit.
Table 8. Sensitivity of “y” to cost–benefit. The static input values of the parameters are: P = 5000, D = 4500, d = 100, Co = 1000, Cp = 50, Ch = 10, CR = 50, Cr = 5, Cs = 3, x = 1.5.
Table 8. Sensitivity of “y” to cost–benefit. The static input values of the parameters are: P = 5000, D = 4500, d = 100, Co = 1000, Cp = 50, Ch = 10, CR = 50, Cr = 5, Cs = 3, x = 1.5.
y0.511.522.5
Percentage of cost–benefit0.7990%0.7975%0.7970%0.7968%0.7966%
Remarksy” is not sensitive to cost–benefit.
Table 9. Sensitivity of raw material cost to cost–benefit. The static input values of the parameters are: P = 5000, D = 4500, d = 100, Co = 1000, Cp = 50, Ch= 10, Cr = 5, Cs = 3, x = 1.5, y = 1.5.
Table 9. Sensitivity of raw material cost to cost–benefit. The static input values of the parameters are: P = 5000, D = 4500, d = 100, Co = 1000, Cp = 50, Ch= 10, Cr = 5, Cs = 3, x = 1.5, y = 1.5.
C R 5055606570
PCB0.797%0.854%0.906%0.954%0.997%
RemarksCost–benefit increases with the increase in raw material cost.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

AlArjani, A.; Miah, M.M.; Uddin, M.S.; Mashud, A.H.M.; Wee, H.-M.; Sana, S.S.; Srivastava, H.M. A Sustainable Economic Recycle Quantity Model for Imperfect Production System with Shortages. J. Risk Financial Manag. 2021, 14, 173. https://doi.org/10.3390/jrfm14040173

AMA Style

AlArjani A, Miah MM, Uddin MS, Mashud AHM, Wee H-M, Sana SS, Srivastava HM. A Sustainable Economic Recycle Quantity Model for Imperfect Production System with Shortages. Journal of Risk and Financial Management. 2021; 14(4):173. https://doi.org/10.3390/jrfm14040173

Chicago/Turabian Style

AlArjani, Ali, Md. Maniruzzaman Miah, Md. Sharif Uddin, Abu Hashan Md. Mashud, Hui-Ming Wee, Shib Sankar Sana, and Hari Mohan Srivastava. 2021. "A Sustainable Economic Recycle Quantity Model for Imperfect Production System with Shortages" Journal of Risk and Financial Management 14, no. 4: 173. https://doi.org/10.3390/jrfm14040173

Article Metrics

Back to TopTop