^{1}

^{2}

^{1}

The economic production quantity model is widely used as a decision making tool. However, climate change and energy supply instability have raised the transport costs of raw materials, thereby affecting production costs. Therefore, this paper presents a new model that considers the holding costs of raw materials under conditions of two-level trade credit and limited storage capacity. Four theorems are developed to characterize the optimal solutions according to a cost-minimization strategy. Finally, a sensitivity analysis is executed to investigate the effects of the various parameters on ordering policies and the annual total relevant costs.

The economic order quantity (EOQ) model [

[

[

Q order size;

P production rate;

D demand rate;

A ordering cost;

T cycle time.

ρ = 1 − D P > 0

L max storage maximum;

s unit selling price per item;

c unit purchasing price per item;

h m unit holding cost per item for raw materials in raw materials warehouse;

h o unit holding cost per item for product in the OW;

h r unit holding cost per item for product in the RW;

I p interest rate payable per $ unit time (year);

I e interest rate earned per $ unit time (year);

t s time in years at which production stops;

M the manufacturer’s trade credit offered by the supplier;

N the customer’s trade credit period offered by the manufacturer;

W storage capacity of the OW;

t w i time when the inventory level increases to W during the production period;

t w d time when the inventory level decreases to W during the production cease period is

T − W D ;

t w d − t w i time of rented warehouse is

{ D T ρ − W P − D + D T ρ − W D , if D T ρ > W 0 , if D T ρ ≤ W ;

T R C ( T ) total relevant cost per unit ti when T > 0 ;

T ∗ optimal solution of T R C ( T ) .

1) Demand rate D is known and constant. me of the model.

2) Production rate P is known and constant, P > D .

3) Shortages are not allowed.

4) A single item is considered.

5) Time period is infinite.

6) h r ≥ h o ≥ h m , M ≥ N and s ≥ c .

7) Storage capacity of raw materials warehouse is unlimited.

8) If the order quantity is larger than manufacturer s OW storage capacity, then the manufacturer will rent an RW in which storage capacity is unlimited. When demand occurs, it first is replenished from the RW which storages those exceeding items. RW takes first in last out (FILO).

9) During the period the account is not settled, generated sales revenue is deposited in and interest-bearing account.

a) When M ≤ T , the account is settled at T = M , the manufacturer pays off all units sold, keeps his or her profits, and starts paying for the higher interest payable on the items in stock with rate I p .

b) When T ≤ M , the account is settled at T = M and the manufacturer does not need to pay any interest payable.

10) If a customer buys an item from the manufacturer at time t ∈ [ 0, T ] , then the customer will have a trade credit period N − t and make the payment at time N.

11) The manufacturer can accumulate revenue and earn interest after his or her customer pays the amount of the purchasing cost to the manufacturer until the end of the trade credit period offered by the supplier. In other words, the manufacturer can accumulate revenue and earn interest during the period from N to M with rate I e under the condition of trade credit.

12) The manufacturer keeps the profit for the use in the other activities.

The model considers three stages of a supply chain system. This paper supposes that the supplier prepares raw materials for production, and the quantity of raw materials is expected to decrease with time (from time 0 to t s ). The quantity of products is expected to increase with time up to the maximum inventory level (from 0 to t s ); the products are sold on demand at the same time. After production stops (at time t s ), the products are sold only on demand until the quantity reaches zero (at time T), as shown in

The annual total relevant cost consists of the following element.

As shown in

( P − D ) t s − D ( T − t s ) = 0 ,

t s = D T P . (1)

0 − Q t s − 0 = − P , 0 ≤ t ≤ T .

Q = P t s = D T . (2)

L max = ( P − D ) × t s = ( P − D ) D T P = D T ρ . (3)

Annual ordering cost is

A T . (4)

Annual purchasing cost is

c Q × 1 T = c D . (5)

Annual holding cost

1) As shown in

h m × Q × t s 2 × 1 T = h m × D T × D T P 2 × 1 T = D 2 T h m 2 P . (6)

2) Two cases occur in annual holding costs of owned warehouse.

a) D T ρ ≤ W , as shown in

Annual holding cost in owned warehouse is

h o × T × L max 2 × 1 T = D T h o ρ 2 . (7)

b) D T ρ > W , as shown in

Annual holding cost in owned warehouse is

h o × [ ( t w d − t w i ) + T ] W 2 × 1 T = W h o − W 2 h o 2 D T ρ . (8)

3) Two cases occur in annual holding costs of rented warehouse.

a) D T ρ ≤ W , as shown in

Annual holding cost in rented warehouse is

0. (9)

b) D T ρ > W , as shown in

Annual holding cost in rented warehouse is

h r × ( t w d − t w i ) × ( L max − W ) 2 × 1 T = h r ( D T ρ − W ) 2 2 D T ρ . (10)

Four cases to occur in costs of annual interest payable for the items kept in stock.

1) 0 < T ≤ N .

Annual interest payable is

0. (11)

2) N ≤ T ≤ M .

Annual interest payable is

0. (12)

3) M ≤ T ≤ P M D , as shown in

Annual interest payable is

c I p × ( ( T − M ) × D ( T − M ) 2 ) × 1 T = c I p D ( T − M ) 2 2 T . (13)

4) M ≤ P M D ≤ T , as shown in

Annual interest payable is

c I p × ( T × D T ρ 2 − M × ( P − D ) M 2 ) × 1 T = c I p ρ ( D T 2 − P M 2 ) 2 T . (14)

Three cases to occur in annual interest earned.

1) 0 < T ≤ N , as shown in

Annual interest earned is

s I e × D T ( M − N ) × 1 T = s I e D ( M − N ) . (15)

2) N ≤ T ≤ M , as shown in

Annual interest earned is

s I e × ( ( D N + D T ) ( T − N ) 2 + D T ( M − T ) ) × 1 T = s I e D ( 2 M T − N 2 − T 2 ) 2 T . (16)

3) N ≤ M ≤ T , as shown in

Annual interest earned is

s I e × ( ( D N + D M ) ( M − N ) 2 ) × 1 T = s I e D ( M 2 − N 2 ) 2 T . (17)

From the above arguments, the annual total relevant cost for the manufacturer can be expressed as T R C ( T ) = annual ordering cost + annual purchasing cost + annual holding cost + annual interest payable + annual interest earned.

Because storage capacity W = D T ρ , there are four cases arise:

1) W D ρ < N ,

2) N ≤ W D ρ < M ,

3) M ≤ W D ρ < P M D ,

4) P M D ≤ W D ρ .

Case 1. W D ρ < N .

According to Equations (1)-(17), the annual total relevant cost T R C ( T ) can be expressed by

T R C ( T ) = { T R C 1 ( T ) , if 0 < T < W D ρ ( 18 a ) T R C 2 ( T ) , if W D ρ ≤ T < N ( 18 b ) T R C 3 ( T ) , if N ≤ T < M ( 18 c ) T R C 4 ( T ) , if M ≤ T < P M D ( 18 d ) T R C 5 ( T ) , if P M D ≤ T (18e)

where

T R C 1 ( T ) = A T + c D + D 2 T h m 2 P + D T h o ρ 2 − s I e D ( M − N ) , (19)

T R C 2 ( T ) = A T + c D + D 2 T h m 2 P + W h o − W 2 h o 2 D T ρ + h r ( D T ρ − W ) 2 2 D T ρ − s I e D ( M − N ) , (20)

T R C 3 ( T ) = A T + c D + D 2 T h m 2 P + W h o − W 2 h o 2 D T ρ + h r ( D T ρ − W ) 2 2 D T ρ − s I e D ( 2 M T − N 2 − T 2 ) 2 T , (21)

T R C 4 ( T ) = A T + c D + D 2 T h m 2 P + W h o − W 2 h o 2 D T ρ + h r ( D T ρ − W ) 2 2 D T ρ + c I p D ( T − M ) 2 2 T − s I e D ( M 2 − N 2 ) 2 T , (22)

T R C 5 ( T ) = A T + c D + D 2 T h m 2 P + W h o − W 2 h o 2 D T ρ + h r ( D T ρ − W ) 2 2 D T ρ + c I p ρ ( D T 2 − P M 2 ) 2 T − s I e D ( M 2 − N 2 ) 2 T . (23)

T R C ( T ) is continuous at T, T ∈ [ 0, ∞ ) because of T R C 1 ( W D ρ ) = T R C 2 ( W D ρ ) , T R C 2 ( N ) = T R C 3 ( N ) , T R C 3 ( M ) = T R C 4 ( M ) , T R C 4 ( P M D ) = T R C 5 ( P M D ) .

Case 2. N ≤ W D ρ < M .

According to Equations (1)-(17), the annual total relevant cost T R C ( T ) can be expressed by

T R C ( T ) = { T R C 1 ( T ) , if 0 < T < N ( 24 a ) T R C 6 ( T ) , if N ≤ T < W D ρ ( 24 b ) T R C 3 ( T ) , if W D ρ ≤ T < M ( 24 c ) T R C 4 ( T ) , if M ≤ T < P M D ( 24 d ) T R C 5 ( T ) , if P M D ≤ T (24e)

where

T R C 6 ( T ) = A T + c D + D 2 T h m 2 P + D T h o ρ 2 − s I e D ( 2 M T − N 2 − T 2 ) 2 T . (25)

T R C ( T ) is continuous at T, T ∈ [ 0, ∞ ) because of T R C 1 ( N ) = T R C 6 ( N ) , T R C 6 ( W D ρ ) = T R C 3 ( W D ρ ) , T R C 3 ( M ) = T R C 4 ( M ) , T R C 4 ( P M D ) = T R C 5 ( P M D ) .

Case 3. M ≤ W D ρ < P M D .

According to Equations (1)-(17), the annual total relevant cost T R C ( T ) can be expressed by

T R C ( T ) = { T R C 1 ( T ) , if 0 < T < N ( 26 a ) T R C 6 ( T ) , if N ≤ T < M ( 26 b ) T R C 7 ( T ) , if M ≤ T < W D ρ ( 26 c ) T R C 4 ( T ) , if W D ρ ≤ T < P M D ( 26 d ) T R C 5 ( T ) , if P M D ≤ T (26e)

where

T R C 7 ( T ) = A T + c D + D 2 T h m 2 P + D T h o ρ 2 + c I p D ( T − M ) 2 2 T − s I e D ( M 2 − N 2 ) 2 T . (27)

T R C ( T ) is continuous at T, T ∈ [ 0, ∞ ) because of T R C 1 ( N ) = T R C 6 ( N ) , T R C 6 ( M ) = T R C 7 ( M ) , T R C 7 ( W D ρ ) = T R C 4 ( W D ρ ) , T R C 4 ( P M D ) = T R C 5 ( P M D ) .

Case 4. P M D ≤ W D ρ .

According to Equations (1)-(17), the annual total relevant cost T R C ( T ) can be expressed by

T R C ( T ) = { T R C 1 ( T ) , if 0 < T < N ( 28 a ) T R C 6 ( T ) , if N ≤ T < M ( 28 b ) T R C 7 ( T ) , if M ≤ T < P M D ( 28 c ) T R C 8 ( T ) , if P M D ≤ T < W D ρ ( 28 d ) T R C 5 ( T ) , if W D ρ ≤ T (28e)

where

T R C 8 ( T ) = A T + c D + D 2 T h m 2 P + D T h o ρ 2 + c I p ρ ( D T 2 − P M 2 ) 2 T − s I e D ( M 2 − N 2 ) 2 T . (29)

T R C ( T ) is continuous at T, T ∈ [ 0, ∞ ) because of T R C 1 ( N ) = T R C 6 ( N ) , T R C 6 ( M ) = T R C 7 ( M ) , T R C 7 ( P M D ) = T R C 8 ( P M D ) , T R C 8 ( W D ρ ) = T R C 5 ( W D ρ ) .

For convenience, all T R C i ( T ) ( i = 1 ∼ 8 ) are defined on T > 0 .

Equations (19)-(23), (25), (27) and (29) yield the first order and second-order derivatives as follows.

T R C ′ 1 ( T ) = − A T 2 + D 2 ( D h m P + ρ h o ) , (30)

T R C ″ 1 ( T ) = 2 A T 3 > 0 , (31)

T R C ′ 2 ( T ) = − 2 A + W 2 ( h o − h r ) D ρ 2 T 2 + D 2 ( D h m P + ρ h r ) , (32)

T R C ″ 2 ( T ) = 2 A + W 2 ( h r − h o ) D ρ T 3 > 0 , (33)

T R C ′ 3 ( T ) = − 2 A + W 2 ( h o − h r ) D ρ − s I e D N 2 2 T 2 + D 2 ( D h m P + ρ h r + s I e ) , (34)

T R C ″ 3 ( T ) = 2 A + W 2 ( h r − h o ) D ρ + s I e D N 2 T 3 > 0 , (35)

T R C ′ 4 ( T ) = − 2 A + W 2 ( h o − h r ) D ρ + D M 2 ( s I e − c I p ) − s I e D N 2 2 T 2 + D 2 ( D h m P + ρ h r + c I p ) , (36)

T R C ″ 4 ( T ) = 2 A + W 2 ( h r − h o ) D ρ + D M 2 ( c I p − s I e ) + s I e D N 2 T 3 , (37)

T R C ′ 5 ( T ) = − 2 A + W 2 ( h o − h r ) D ρ + D M 2 ( s I e − c I p ) − s I e D N 2 + c I p P M 2 2 T 2 + D 2 ( D h m P + ρ ( h r + c I p ) ) , (38)

T R C ″ 5 ( T ) = 2 A + W 2 ( h r − h o ) D ρ + D M 2 ( c I p − s I e ) + s I e D N 2 − c I p P M 2 T 3 , (39)

T R C ′ 6 ( T ) = − 2 A − s I e D N 2 2 T 2 + D 2 ( D h m P + ρ h o + s I e ) , (40)

T R C ″ 6 ( T ) = 2 A + s I e D N 2 T 3 > 0 , (41)

T R C ′ 7 ( T ) = − 2 A + D M 2 ( s I e − c I p ) − s I e D N 2 2 T 2 + D 2 ( D h m P + ρ h o + c I p ) , (42)

T R C ″ 7 ( T ) = 2 A + D M 2 ( c I p − s I e ) + s I e D N 2 T 3 , (43)

T R C ′ 8 ( T ) = − 2 A + D M 2 ( s I e − c I p ) − s I e D N 2 + c I p P M 2 2 T 2 + D 2 ( D h m P + ρ ( h o + c I p ) ) , (44)

and

T R C ″ 8 ( T ) = 2 A + D M 2 ( c I p − s I e ) + s I e D N 2 − c I p P M 2 T 3 . (45)

Let

G 4 = 2 A + W 2 ( h r − h o ) D ρ + D M 2 ( c I p − s I e ) + s I e D N 2 , (46)

G 5 = 2 A + W 2 ( h r − h o ) D ρ + D M 2 ( c I p − s I e ) + s I e D N 2 − c I p P M 2 , (47)

G 7 = 2 A + D M 2 ( c I p − s I e ) + s I e D N 2 , (48)

and

G 8 = 2 A + D M 2 ( c I p − s I e ) + s I e D N 2 − c I p P M 2 . (49)

Equations (46)-(49) imply

G 8 < G 5 < G 4 , (50)

and

G 8 < G 7 < G 4 . (51)

Equations (30)-(45) reveal the following results.

Lemma 1. 1) T R C i ( T ) is convex on T > 0 if i = 1 , 2 , 3 , 6 .

2) T R C i ( T ) is convex on T > 0 if G i > 0 . Otherwise T R C ′ i ( T ) is increasing on T > 0 for all i = 4 , 5 , 7 , 8 . Solving

T R C ′ i ( T ) = 0, i = 1 ∼ 8 (52)

then

T 1 * = 2 A D ( D h m P + ρ h o ) , (53)

T 2 * = 2 A + W 2 ( h r − h o ) D ρ D ( D h m P + ρ h r ) , (54)

T 3 * = 2 A + W 2 ( h r − h o ) D ρ + s I e D N 2 D ( D h m P + ρ h r + s I e ) , (55)

T 4 * = 2 A + W 2 ( h r − h o ) D ρ + D M 2 ( c I p − s I e ) + s I e D N 2 D ( D h m P + ρ h r + c I p ) if G 4 > 0 , (56)

T 5 * = 2 A + W 2 ( h r − h o ) D ρ + D M 2 ( c I p − s I e ) + s I e D N 2 − c I p M 2 D ( D h m P + ρ ( h r + c I p ) ) if G 5 > 0 , (57)

T 6 * = 2 A + s I e D N 2 D ( D h m P + ρ h o + s I e ) , (58)

T 7 * = 2 A + D M 2 ( c I p − s I e ) + s I e D N 2 D ( D h m P + ρ h o + c I p ) if G 7 > 0 , (59)

and

T 8 * = 2 A + D M 2 ( c I p − s I e ) + s I e D N 2 − c I p M 2 D ( D h m P + ρ ( h o + c I p ) ) if G 8 > 0 , (60)

are the respective solutions of Equation (52). Furthermore, if T i * exists, then T R C i ( T ) is convex on T > 0 and

T R C ′ i ( T ) = { < 0 if 0 < T < T i * ( 61 a ) = 0 if T = T i * ( 61 b ) > 0 if T i * < T (61c)

Equations (61a)-(61c) imply that T R C i ( T ) is decreasing on ( 0, T i * ] and increasing on [ T i * , ∞ ) for all i = 1 ∼ 8 .

Case 1. W D ρ < N .

Equations (30), (32), (34), (36) and (38) yield

T R C ′ 1 ( W D ρ ) = T R C ′ 2 ( W D ρ ) = Δ 12 2 ( W D ρ ) 2 , (62)

T R C ′ 2 ( N ) = T R C ′ 3 ( N ) = Δ 23 2 N 2 , (63)

T R C ′ 3 ( M ) = T R C ′ 4 ( M ) = Δ 34 2 M 2 , (64)

T R C ′ 4 ( P M D ) = T R C ′ 5 ( P M D ) = Δ 45 2 ( P M D ) 2 , (65)

where

Δ 12 = − 2 A + D ( W D ρ ) 2 ( D h m P + ρ h o ) , (66)

Δ 23 = − 2 A + W 2 ( h o − h r ) D ρ + D N 2 ( D h m P + ρ h r ) , (67)

Δ 34 = − 2 A + W 2 ( h o − h r ) D ρ − s I e D N 2 + D M 2 ( D h m P + ρ h r + s I e ) , (68)

Δ 45 = − 2 A + W 2 ( h o − h r ) D ρ + D M 2 ( s I e − c I p ) − s I e D N 2 + D ( P M D ) 2 ( D h m P + ρ h r + c I p ) . (69)

Equations (66)-(69) imply

Δ 12 < Δ 23 < Δ 34 < Δ 45 . (70)

Case 2. N ≤ W D ρ < M .

Equations (30), (34), (36), (38) and (40) yield

T R C ′ 1 ( N ) = T R C ′ 6 ( N ) = Δ 16 2 N 2 , (71)

T R C ′ 6 ( W D ρ ) = T R C ′ 3 ( W D ρ ) = Δ 63 2 ( W D ρ ) 2 , (72)

T R C ′ 3 ( M ) = T R C ′ 4 ( M ) = Δ 34 2 M 2 , (73)

T R C ′ 4 ( P M D ) = T R C ′ 5 ( P M D ) = Δ 45 2 ( P M D ) 2 , (74)

where

Δ 16 = − 2 A + D N 2 ( D h m P + ρ h o ) , (75)

Δ 63 = − 2 A − s I e D N 2 + D ( W D ρ ) 2 ( D h m P + ρ h o + s I e ) . (76)

Equations (68), (69), (75) and (76) imply

Δ 16 ≤ Δ 63 < Δ 34 < Δ 45 . (77)

Case 3. M ≤ W D ρ < P M D .

Equations (30), (36), (38), (40) and (42) yield

T R C ′ 1 ( N ) = T R C ′ 6 ( N ) = Δ 16 2 N 2 , (78)

T R C ′ 6 ( M ) = T R C ′ 7 ( M ) = Δ 67 2 M 2 , (79)

T R C ′ 7 ( W D ρ ) = T R C ′ 4 ( W D ρ ) = Δ 74 2 ( W D ρ ) 2 , (80)

T R C ′ 4 ( P M D ) = T R C ′ 5 ( P M D ) = Δ 45 2 ( P M D ) 2 , (81)

where

Δ 67 = − 2 A − s I e D N 2 + D M 2 ( D h m P + ρ h o + s I e ) , (82)

Δ 74 = − 2 A + D M 2 ( s I e − c I p ) − s I e D N 2 + D ( W D ρ ) 2 ( D h m P + ρ h o + c I p ) . (83)

Equations (64), (70), (77) and (78) imply

Δ 16 ≤ Δ 67 ≤ Δ 74 < Δ 45 . (84)

Case 4. P M D ≤ W D ρ .

Equations (30), (38), (40), (42) and (44) yield

T R C ′ 1 ( N ) = T R C ′ 6 ( N ) = Δ 16 2 N 2 , (85)

T R C ′ 6 ( M ) = T R C ′ 7 ( M ) = Δ 67 2 M 2 , (86)

T R C ′ 7 ( P M D ) = T R C ′ 8 ( P M D ) = Δ 78 2 ( P M D ) 2 , (87)

T R C ′ 8 ( W D ρ ) = T R C ′ 5 ( W D ρ ) = Δ 45 2 ( W D ρ ) 2 , (88)

where

Δ 78 = − 2 A + D M 2 ( s I e − c I p ) − s I e D N 2 + D ( P M D ) 2 ( D h m P + ρ h o + c I p ) , (89)

Δ 85 = − 2 A + D M 2 ( s I e − c I p ) − s I e D N 2 + c I p P M 2 + D ( W D ρ ) 2 ( D h m P + ρ ( h o + c I p ) ) . (90)

Equations (69), (75), (89) and (90) imply

Δ 16 ≤ Δ 67 ≤ Δ 78 ≤ Δ 85 . (91)

Based on the above arguments, the following results holds.

Lemma 2. (A) If Δ_{34} ≤ 0, then

(a1) G 4 > 0 ,

(b2) T 4 * exists,

(c3) T R C 4 ( T ) is convex on T > 0 .

(B) If Δ 45 ≤ 0 , then

(b1) G 4 > 0 and G 5 > 0 ,

(b2) T 4 * and T 5 * exist,

(b3) T R C 4 ( T ) and T R C 5 ( T ) are convex on T > 0 .

(C) If Δ 67 ≤ 0 , then

(c1) G 7 > 0 ,

(c2) T 7 * exists,

(c3) T R C 7 ( T ) is convex on T > 0 .

(D) If Δ 74 ≤ 0 , then

(d1) G 4 > 0 and G 7 > 0 ,

(d2) T 4 * and T 7 * exist,

(d3) T R C 4 ( T ) and T R C 7 ( T ) are convex on T > 0 .

(E) If Δ 78 ≤ 0 , then

(e1) G 7 > 0 and G 8 > 0 ,

(e2) T 7 * and T 8 * exist,

(e3) T R C 7 ( T ) and T R C 8 ( T ) are convex on T > 0 .

(F) If Δ 85 ≤ 0 , then

(f1) G 5 > 0 and G 8 > 0 ,

(f2) T 5 * and T 8 * exist,

(f3) T R C 5 ( T ) and T R C 8 ( T ) are convex on T > 0 .

Proof

(A) (a1) If Δ 34 ≤ 0 , then

2 A ≥ W 2 ( h o − h r ) D ρ − s I e D N 2 + D M 2 ( D h m P + ρ h r + s I e ) . (92)

Equation (92) implies

G 4 ≥ D M 2 ( D h m P + ρ h r + c I p ) > 0. (93)

(a2) Equation (56) and lemma 1 imply that T 4 * exists.

(a3) Equation (37) and lemma 1 imply that T R C 4 ( T ) is convex on T > 0 .

(B) (b1) If Δ 45 ≤ 0 , then

2 A ≥ W 2 ( h o − h r ) D ρ + D M 2 ( s I e − c I p ) − s I e D N 2 + D ( P M D ) 2 ( D h m P + ρ h r + c I p ) . (94)

Equation (94) implies

G 5 ≥ D ( P M D ) 2 ( D h m P + ρ ( h r + c I p ) ) > 0. (95)

Equations (50) and (95) demonstrate G 4 > G 5 > 0 .

(b2) Equations (56), (57) and lemma 1 imply that T 4 * and T 5 * exist.

(b3) Equations (37), (39) and lemma 1 imply that T R C 4 ( T ) and T R C 5 ( T ) are convex on T > 0 .

(C) (c1) If Δ 67 ≤ 0 , then

2 A ≥ s I e D N 2 + D M 2 ( D h m P + ρ h o + s I e ) , (96)

Equation (96) implies

G 7 ≥ D M 2 ( D h m P + ρ h o + c I p ) > 0. (97)

(c2) Equation (59) and lemma 1 imply that T 7 * exists.

(c3) Equation (43) and lemma 1 imply that T R C 7 ( T ) is convex on T > 0 .

(D) (d1) If Δ 74 ≤ 0 , then

2 A ≥ D M 2 ( s I e − c I p ) − s I e D N 2 + D ( W D ρ ) 2 ( D h m P + ρ h o + c I p ) . (98)

Equation (98) implies

G 7 ≥ D ( W D ρ ) 2 ( D h m P + ρ h o + c I p ) > 0. (99)

Equations (51) and (99) demonstrate G 4 > G 7 > 0 .

(d2) Equations (56), (59) and lemma 1 imply that T 4 * and T 7 * exist.

(d3) Equations (37), (43) and lemma 1 imply that T R C 4 ( T ) and T R C 7 ( T ) are convex on T > 0 .

(E) (e1) If Δ 78 ≤ 0 , then

2 A ≥ D M 2 ( s I e − c I p ) − s I e D N 2 + D ( P M D ) 2 ( D h m P + ρ h o + c I p ) . (100)

Equation (100) implies

G 8 ≥ D ( P M D ) 2 ( D h m P + ρ ( h o + c I p ) ) > 0. (101)

Equations (51) and (101) demonstrate G 7 > G 8 > 0 .

(e2) Equations (59), (60) and lemma 1 imply that T 7 * and T 8 * exist.

(e3) Equations (43), (45) and lemma 1 imply that T R C 7 ( T ) and T R C 8 ( T ) are convex on T > 0 .

(F) (f1) If Δ 85 ≤ 0 , then

2 A ≥ D M 2 ( s I e − c I p ) − s I e D N 2 + c I p P M 2 + D ( W D ρ ) 2 ( D h m P + ρ ( h o + c I p ) ) . (102)

Equation (102) implies

G 8 ≥ D ( W D ρ ) 2 ( D h m P + ρ ( h o + c I p ) ) > 0. (103)

Equations (51) and (103) demonstrate G 5 > G 8 > 0 .

(f2) Equations (57), (60) and lemma 1 imply that T 5 * and T 8 * exist.

(f3) Equations (39), (45) and lemma 1 imply that T R C 5 ( T ) and T R C 8 ( T ) are convex on T > 0 .

Incorporate the above arguments, we have completed the proof of lemma 2.

Theorem 1. Suppose W D ρ < N . Hence,

(A) if 0 < Δ 12 , then T R C ( T * ) = T R C 1 ( T 1 * ) and T * = T 1 * .

(B) if Δ 12 ≤ 0 < Δ 23 , then T R C ( T * ) = T R C 2 ( T 2 * ) and T * = T 2 * .

(C) if Δ 23 ≤ 0 < Δ 34 , then T R C ( T * ) = T R C 3 ( T 3 * ) and T * = T 3 * .

(D) if Δ 34 ≤ 0 < Δ 45 , then T R C ( T * ) = T R C 4 ( T 4 * ) and T * = T 4 * .

(E) if Δ 45 ≤ 0 , then T R C ( T * ) = T R C 5 ( T 5 * ) and T * = T 5 * .

Proof (A) If 0 < Δ 12 , then 0 < Δ 12 < Δ 23 < Δ 34 < Δ 45 . So, Equations (61a)-(61c), lemma 1 and 2 imply

(a1) T R C 1 ( T ) is decreasing on ( 0, T 1 * ] and increasing on [ T 1 * , W D ρ ] .

(a2) T R C 2 ( T ) is increasing on [ W D ρ , N ] .

(a3) T R C 3 ( T ) is increasing on [ N , M ] .

(a4) T R C 4 ( T ) is increasing on [ M , P M D ] .

(a5) T R C 5 ( T ) is increasing on [ P M D , ∞ ) .

Since T R C ( T ) is continuous on T > 0 , Equations (18a)-(18e) and (a1)-(a5) reveal that T R C ( T ) is decreasing on ( 0, T 1 * ] and increasing on [ T 1 * , ∞ ) . Hence, T * = T 1 * and T R C ( T * ) = T R C 1 ( T 1 * ) .

(B) If Δ 12 ≤ 0 < Δ 23 , then Δ 12 ≤ 0 < Δ 23 < Δ 34 < Δ 45 . So, Equations (61a)-(61c), lemma 1 and 2 imply.

(b1) T R C 1 ( T ) is decreasing on [ 0, W D ρ ] .

(b2) T R C 2 ( T ) is decreasing on [ W D ρ , T 2 * ] and increasing on [ T 2 * , N ] .

(b3) T R C 3 ( T ) is increasing on [ N , M ] .

(b4) T R C 4 ( T ) is increasing on [ M , P M D ] .

(b5) T R C 5 ( T ) is increasing on [ P M D , ∞ ) .

Since T R C ( T ) is continuous on T > 0 , Equations (18a)-(18e) and (a1)-(a5) reveal that T R C ( T ) is decreasing on ( 0, T 2 * ] and increasing on [ T 2 * , ∞ ) . Hence, T * = T 2 * and T R C ( T * ) = T R C 2 ( T 2 * ) .

(C) If Δ 23 ≤ 0 < Δ 34 , then Δ 12 < Δ 23 ≤ 0 < Δ 34 < Δ 45 . So, Equations (61a)-(61c), lemma 1 and 2 imply.

(c1) T R C 1 ( T ) is decreasing on [ 0, W D ρ ] .

(c2) T R C 2 ( T ) is decreasing on [ W D ρ , N ] .

(c3) T R C 3 ( T ) is decreasing on [ N , T 3 * ] and increasing on [ T 3 * , M ] .

(c4) T R C 4 ( T ) is increasing on [ M , P M D ] .

(c5) T R C 5 ( T ) is increasing on [ P M D , ∞ ) .

Since T R C ( T ) is continuous on T > 0 , Equations (18a)-(18e) and (c1)-(c5) reveal that T R C ( T ) is decreasing on ( 0, T 3 * ] and increasing on [ T 3 * , ∞ ) . Hence, T * = T 3 * and T R C ( T * ) = T R C 3 ( T 3 * ) .

(D) If Δ 34 ≤ 0 < Δ 45 , then Δ 12 < Δ 23 < Δ 34 ≤ 0 < Δ 45 . So, Equations (61a)-(61c), lemma 1 and 2 imply

(d1) T R C 1 ( T ) is decreasing on [ 0, W D ρ ] .

(d2) T R C 2 ( T ) is decreasing on [ W D ρ , N ] .

(d3) T R C 3 ( T ) is decreasing on [ N , M ] .

(d4) T R C 4 ( T ) is decreasing on [ M , T 4 * ] and increasing on [ T 4 * , P M D ] .

(d5) T R C 5 ( T ) is increasing on [ P M D , ∞ ) .

Since T R C ( T ) is continuous on T > 0 , Equations (18a)-(18e) and (d1)-(d5) reveal that T R C ( T ) is decreasing on ( 0, T 4 * ] and increasing on [ T 4 * , ∞ ) . Hence, T * = T 4 * and T R C ( T * ) = T R C 4 ( T 4 * ) .

(E) If Δ 45 ≤ 0 , then Δ 12 < Δ 23 < Δ 34 < Δ 45 ≤ 0 . So, Equations (61a)-(61c), lemma 1 and 2 imply.

(e1) T R C 1 ( T ) is decreasing on [ 0, W D ρ ] .

(e2) T R C 2 ( T ) is decreasing on [ W D ρ , N ] .

(e3) T R C 3 ( T ) is decreasing on [ N , M ] .

(e4) T R C 4 ( T ) is decreasing on [ M , T 4 * ] .

(e5) T R C 5 ( T ) is decreasing on [ P M D , T 5 * ] and increasing on [ T 5 * , ∞ ) .

Since T R C ( T ) is continuous on T > 0 , Equations (18a)-(18e) and (e1)-(e5) reveal that T R C ( T ) is decreasing on ( 0, T 5 * ] and increasing on [ T 5 * , ∞ ) . Hence, T * = T 5 * and T R C ( T * ) = T R C 5 ( T 5 * ) .

Incorporating all argument above arguments, we have completed the proof of Theorem 1.

Theorem 2. Suppose N ≤ W D ρ < M . Hence,

(A) if 0 < Δ 16 , then T R C ( T * ) = T R C 1 ( T 1 * ) and T * = T 1 * .

(B) if Δ 16 ≤ 0 < Δ 63 , then T R C ( T * ) = T R C 6 ( T 6 * ) and T * = T 6 * .

(C) if Δ 63 ≤ 0 < Δ 34 , then T R C ( T * ) = T R C 3 ( T 3 * ) and T * = T 3 * .

(D) if Δ 34 ≤ 0 < Δ 45 , then T R C ( T * ) = T R C 4 ( T 4 * ) and T * = T 4 * .

(E) if Δ 45 ≤ 0 , then T R C ( T * ) = T R C 5 ( T 5 * ) and T * = T 5 * .

Applying lemmas 1, 2 and Equations (24a)-(24e), the following results hold.

Theorem 3. Suppose M ≤ W D ρ < P M D . Hence,

(A) if 0 < Δ 16 , then T R C ( T * ) = T R C 1 ( T 1 * ) and T * = T 1 * .

(B) if Δ 16 ≤ 0 < Δ 67 , then T R C ( T * ) = T R C 6 ( T 6 * ) and T * = T 6 * .

(C) if Δ 67 ≤ 0 < Δ 74 , then T R C ( T * ) = T R C 7 ( T 7 * ) and T * = T 7 * .

(D) if Δ 74 ≤ 0 < Δ 45 , then T R C ( T * ) = T R C 4 ( T 4 * ) and T * = T 4 * .

(E) if Δ 45 ≤ 0 , then T R C ( T * ) = T R C 5 ( T 5 * ) and T * = T 5 * .

Applying lemmas 1, 2 and Equations (26a)-(26e), the following results hold.

Theorem 4. Suppose P M D ≤ W D ρ . Hence,

(A) if 0 < Δ 16 , then T R C ( T * ) = T R C 1 ( T 1 * ) and T * = T 1 * .

(B) if Δ 16 ≤ 0 < Δ 67 , then T R C ( T * ) = T R C 6 ( T 6 * ) and T * = T 6 * .

(C) if Δ 67 ≤ 0 < Δ 78 , then T R C ( T * ) = T R C 7 ( T 7 * ) and T * = T 7 * .

(D) if Δ 78 ≤ 0 < Δ 85 , then T R C ( T * ) = T R C 8 ( T 8 * ) and T * = T 8 * .

(E) if Δ 85 ≤ 0 , then T R C ( T * ) = T R C 5 ( T 5 * ) and T * = T 5 * .

Applying lemmas 1, 2 and Equations (28a)-(28e), the following results hold.

We present the sensitivity analysis by Maple 18.00 for determining the unique solution T i * when T R C ′ i ( T ) = 0, i = 1 ∼ 8 , given P = 4500 units / year , D = 3000 units / year , A = $ 1000 / order , s = $ 30 / unit , c = $ 10 / unit , h m = $ 0.5 / unit / year , h o = $ 1.5 / unit / year , h r = $ 5 / unit / year , I p = $ 0.3 / year , I e = $ 0.08 / year , M = 90 days = 90 / 365 year , N = 45 days = 45 / 365 year , W = 350 units . We increase/decrease 25% and 50% of the parameters simultaneously to execute the sensitivity analysis.

Although traditional models facilitate calculation, they ignore the importance of raw materials. However, the related costs of raw materials directly affect the total relevant costs, thereby generating significant errors in the overall consideration.

Therefore, this paper presents a new inventory model that applies raw materials with two-level of trade credit―the finite replenishment rate and limited storage capacity―by considering the trade credit condition proposed by [

Impact | T^{*} | TRC(T) | ||
---|---|---|---|---|

this paper | [ | this paper | [ | |

Positive & Major | A | A | c | A |

Positive & Minor | h_{r}, I_{e} | h_{r}, I_{e} | A, s, h_{m}, h_{o}, I_{p} | s, c, h_{o}, I_{p} |

Negative & Minor | s, c, h_{m}, h_{o} | s, c, h_{o} | h_{r}, I_{e} | h_{r}, I_{e} |

Negative & Major | I_{p} | I_{p} | - | - |

The authors declare no conflicts of interest regarding the publication of this paper.

Yen, G.-F., Lin, S.-D. and Lee, A.-K. (2019) EPQ Policies Considering the Holding Cost of Raw Materials under Conditions of Two-Level Trade Credit and Limited Storage Capacity. Open Access Library Journal, 6: e5140. https://doi.org/10.4236/oalib.1105140