Authors:
Seyed Ali Hasheminejad School of Industrial Engineering, Iran University of Science and Technology, Tehran, Iran

Search for other papers by Seyed Ali Hasheminejad in
Current site
Google Scholar
PubMed
Close
https://orcid.org/0000-0001-5681-3818
,
Khadijeh Valipour School of Industrial Engineering, Iran University of Science and Technology, Tehran, Iran

Search for other papers by Khadijeh Valipour in
Current site
Google Scholar
PubMed
Close
, and
Hamid Khoshnood School of Industrial Engineering, Amirkabir University of Technology, Tehran, Iran

Search for other papers by Hamid Khoshnood in
Current site
Google Scholar
PubMed
Close
Open access

Abstract

Supply chain management intends to integrate supply chains' activities such as material flow, information flow and financial issues. Material flow management is the most significant issue since the inventory level in the whole supply chain could be optimized by an integrated plan. In other words, when one member of the supply chain plans to reduce its inventory level solely, despite reducing inventory in this node the inventory will be stocked in other partners' warehouses. Therefore, in this paper a new mathematical model has been developed to facilitate the process of finding the optimum solution in economic production, purchase and delivery lots and their schedules in a three-echelon supply chain environment; including raw material in suppliers, manufacturer and assembly facility as a customer. The manufacturer with a flow shop system provides its requirements from supplier, assemble multiple products, and delivers products to the customer (automotive OEM alike) on an optimum multiple delivery points. The delivery cycles would be identified through the production common cycle regarding the supply chain flexibility. Finally, a modified real-valued Genetic Algorithm (MRGA), and an Optimal Enumeration Method (OEM) are developed, and some numerical experiments have been done and compared as well.

Abstract

Supply chain management intends to integrate supply chains' activities such as material flow, information flow and financial issues. Material flow management is the most significant issue since the inventory level in the whole supply chain could be optimized by an integrated plan. In other words, when one member of the supply chain plans to reduce its inventory level solely, despite reducing inventory in this node the inventory will be stocked in other partners' warehouses. Therefore, in this paper a new mathematical model has been developed to facilitate the process of finding the optimum solution in economic production, purchase and delivery lots and their schedules in a three-echelon supply chain environment; including raw material in suppliers, manufacturer and assembly facility as a customer. The manufacturer with a flow shop system provides its requirements from supplier, assemble multiple products, and delivers products to the customer (automotive OEM alike) on an optimum multiple delivery points. The delivery cycles would be identified through the production common cycle regarding the supply chain flexibility. Finally, a modified real-valued Genetic Algorithm (MRGA), and an Optimal Enumeration Method (OEM) are developed, and some numerical experiments have been done and compared as well.

1 Introduction

Recently, it has been well perceived by individual businesses to work together as a supply chain in contrast with their conventional competitions of those as solely autonomous entities. Supply chain management has been developed to integrate the related activities such as material flow, information flow and financial issues in order to obtain reliable and long-term competitive advantages. Material flow management is the most significant issue in the supply chain concept. Whereas, the inventory level in the whole supply chain can be optimized by coordinated planning. When one member of the supply chain intends to reduce its inventory level solely, consequently the supply chain's inventory management level will be weakened. To put it another way, while the inventory of one node of a supply chain is reduced, this inventory will be stocked in other suppliers or distributor warehouses. Hence, an integrated plan is critical to schedule all production and deliveries in accordance with customer's demands.

Recently, the ELSP (Economic Lot Scheduling Problem) has been extended to the supply chain as ELDSP (Economic Lot and Delivery Problem). By optimizing production and delivery schedule, the supply chain's players will be harmonized and accordingly, inventory level will be minimized. Regarding the importance of coordination in the whole supply chain and high effect on competitive advantage, this problem has become popular in the academic area. Consequently, many researchers have paid attention to material flow management throughout the supply chain.

In this paper, it is desired to apply the ELSP model in a three-level supply chain which the manufacturer works in multi stages-multi products manner, a case that usually occurs in the real world. The key challenge is the changes, which should be applied in traditional ELSP to handle this environment and deliver a product with lower volume and high frequency. Concerning this issue, a new mathematical model has been developed to find optimum economic production, purchase and delivery lots and schedules in a three-echelon supply chain environment. Three-echelon supply chain consists of raw material supplier, manufacturer (could be automotive part supplier) and assembly facility (OEM) as a customer. The manufacturer provides its requirements from its supplier, transforms them into multiple products, and delivers products to customers on optimum delivery points.

2 Material and methods

Economic Lot and scheduling Problem was presented by Rogers in 1958 [1]. He assumed constant and continuous demand for each product, single stage production system and production capacity was considered finite. The ELSP is an important problem in scheduling which has been enormously considered by researchers during the last decades and due to its complexity, it is known as NP-hard problem demonstrated by Hsu in 1983 [2]. Santander-Mercado and Maria Jubiz-Diaz [3] surveyed ELSP and classified this problem in different aspects. They divided the problem into single facility and multi facility, based on the production system and regarding the scheduling approach, three different approaches were introduced: Common cycle approach, Basic/extended basic approach and Time-varying lot size approach. Huang, Yao [4], Ouenniche, and Boctor [5], each proposed two new heuristics approaches to solve the mathematical model for the ELSP with a flow shop system as a multi-product, multi-stage, economic lot-sizing problem. Recently, ELSP has been extended to sustainability challenges, Beck et al. [6] studied the effects of energy usage on ELSP with different approaches. Ferretti [7] considers ELSP with returns and sorting line regarding environmental legislations. Furthermore, Alaghebandha et al. [8] have explored ELSP to new branch and consider distributed flow shop in multi factories. As industrial competition increased drastically, it would be widely accepted that the competition is not just among companies and it is rival among supply chains. As ELSP entered in supply chain environment, it was named Economic Lot and Delivery Scheduling Problem. In the problem there is a simple supply chain consisting of one supplier and one customer, the supplier is in charge of producing some products and delivers them to the customer. This problem aims to optimize production and delivery lots and their schedule from supplier to the costumer.

Hahm and Yano [9, 10] have conducted some studies on Two-level supply chain that established the relationship between one supplier and one customer. Moreover, they developed two heuristic algorithms under common cycle policy. In the first algorithm, production and delivery cycles are the same but in the second one there were multi deliveries during a production cycle. Torabi et al. [11, 12] developed a new mathematical model in flexible flow shop and job shop manufacturing system with an assumption of constant demand rate and finite planning horizon. They presented two effective algorithms; Optimal Enumeration Method and algorithm genetic as a meta-heuristic method. Whereas, both production and delivery cycles are assumed equivalent. Furthermore, Torabi and Jenabi [13] extended the two-echelon model under two variants of the basic period approach and considered the parallel machines in each stage in the flow shop system. Finally, two efficient hybrid genetic algorithms have been applied. Lee et al. [14] developed a new model with infinite planning horizon. They considered a manufacturer who has a single-stage production system. The manufacturer purchases the required raw material from one supplier and after transformation to final product, delivers them to the customer. They upgraded a broader supply chain and considered the process of raw material purchase. Kim et al. [15] studied a three-level supply chain including a raw material supplier, manufacturer and multi customers. The manufacturer works with single stage–multi product system under a common cycle policy and all products have a common raw material. They introduced an efficient heuristic algorithm as a solution as well.

Sarker and Diponegoro [16] developed an optimal policy in three-level supply chain environment consisting of multi-suppliers, a manufacturer with single-facility single-product system and multi-buyers. The number of production cycles was given but their lengths were not necessarily the same and the manufacturer delivers the product to the buyers in fixed intervals and lots. As a solution method, an analytical method was derived to ensure higher potential to save cost in a supply chain. Jaber and Goyal [17] developed a mathematical model in a three-level supply chain with multiple suppliers, a manufacturer and multiple buyers. In their research, the manufacturer buys the raw material and delivers the finished product to the buyer in a cyclic manner. The manufacturer's production system was a single product-single machine in infinite planning horizon in which the product requires k items from m suppliers. Eventually, a centralized decision process as an effective solution has been proposed. Osman and Demirli [18] proposed a new three-level supply chain model, in which two of three stages consist of multiple suppliers and the third stage is a single assembly facility. They employed integer multipliers policy instead of the common cycle time. The inventory profile of the finished products in the assembly facility is not considered in the model and components at each supplier node are produced on a single-stage production system. In addition, all produced components are shipped in one lot at the end of the production cycle and the production setup cost and times at the assembly facility are neglected. In the end, a hybrid algorithm is applied to solve the developed nonlinear model. Goli and Alinaghian [19] have developed a new mathematical model where a manufacturer with an open-shop system procured raw materials from suppliers and after conversion into the final product transfer them to packaging companies. El-meehy et al. [20] developed a new mathematical model of incapacitated three-echelon supply chain and considered two raw material suppliers and two products. Finally, the model has been solved by Lingo and multiple samples by various capacities and has been experimented.

To the best of our knowledge, rare studies have been done to apply coordination in a three-level supply chain in spite of its high priority. Therefore, in this paper a new mathematical model is developed. There is a manufacturer with finite planning horizon and multi stage-multi product production system, which procures its optimum quantity of raw material from suppliers and, after assembling finished goods by flow shop system, delivers them to customers in accordance with optimum delivery cycles. The production system consists of production cost, work in process (WIP) holding cost, setup cost, purchase cost, raw material holding cost and finished goods delivery cost. As solution methodologies, an Optimal Enumeration Method (OEM) as heuristics method and a Modified Real-valued Genetic Algorithm (MRGA) as a metaheuristics method are proposed. In addition, the numerical examples are presented.

The remainder of this paper is organized as follows. In Section 3, a mathematical model is described. Section 4 proposes solution methods MRGA & OEM. Finally, in the last section numerical examples,conclusion and future study are provided.

In this paper, a new mathematical model is developed which shows a manufacturer with multi stage - multi product system and finite horizon time. The manufacturer procures its requirements from multiple suppliers and after transforming them to final product by flow shop system, delivers them to a customer. The objective function in this supply chain consists of eight costs while each of them would affect the model's output. The cost functions in the whole supply chain are:

  1. Raw material holding cost in suppliers' warehouse

  2. Raw material purchasing cost (consists of transportation cost except raw material price)

  3. Raw material holding cost in flow shop system

  4. Work in Process holding cost

  5. Setup cost

  6. Finished goods holding cost in flow shop system

  7. Finished goods delivery cost

  8. Finished goods holding cost in assembly facility warehouse

All costs are illustrated in Fig. 1.

Fig. 1.
Fig. 1.

Costs in considered supply chain

Citation: International Review of Applied Sciences and Engineering 13, 2; 10.1556/1848.2021.00334

The common cycle policy is applied in this model (i.e., the cycle times of all products are equal and there is a common cycle for all products). Therefore, an optimum common cycle time (T * ) should be calculated for the first step and, next step the production lot size of each product will be derivated (Q i = D i .T * ). Consequently, the purchase common cycle should be determined in a way that one lot of each product will be purchased in it. According to type of products in the automotive industry and associated logistics cost, they are usually provided more than the requirement of one production cycle. Hence, it may be assumed that the purchase common cycle is an integer multiple of the production common cycle (B = e.T).

It is supposed that each production common cycle (T) can consists of multiple delivery cycles (D), in order to decrease the finished goods holding cost at the end of flow shop system. The lots of products that have been finished before each delivery point are delivered to the customer. Otherwise, they should wait until the next delivery point. The relationship between production and delivery cycle is D * = T * /k (k is an integer).

The delivered product lot to a customer can cover its requirement for period T. So during the production common cycle the products are delivered to the customer based on production sequence. In other words, the type of product would be different in each delivery and this delivery sequence duplicates in the next cycles. In Fig. 2, the delivery points of products are depicted in a case that the production sequence is (2-1-3-4) with three delivery cycles in a production cycle (k = 3). Products 1 and 2 are delivered at the first delivery point while products 3 and 4 are delivered at the second and third points, respectively. It is noted that product 4 is delivered to the customer at the third delivery point because its lot was not prepared at the second point and this sequence will be repeated in the next steps.

Fig. 2.
Fig. 2.

Delivery cycle scheme

Citation: International Review of Applied Sciences and Engineering 13, 2; 10.1556/1848.2021.00334

There is the same policy to identify purchase lot size (B * ). It is desirable to find a common purchase cycle in which just one lot of each product is purchased. In accordance with this type of industry, such as automotive manufacturers, the raw materials are provided more than the requirement of one production cycle. Furthermore, to simplify the problem we assume that the purchase cycle is an integer multiple of production cycle time (B = e.T). In the next section, the mathematical model is developed. The assumption, parameters and variables are described below:

Problem assumptions:

  • Each component should be processed by one machine at each stage.

  • The lot size of each component is equal in all stages.

  • Machines are persistently available and each one can only process one component at a given time.

  • Backlogging is not allowed.

  • Setup times and costs are considered in flow shop system.

  • Setup times and costs for AF are negligible.

  • Delivery lead-time is assumed zero.

  • The sequence of Production for each machine at each stage is unique and is optimized by the model.

  • Both supplier and assembly facility have linear inventory holding costs on finished goods.

  • The supplier has linear inventory holding costs on WIP components.

  • Preemption is not allowed.

  • Lot splitting is not allowed.

  • There are unlimited buffers between successive stages and in-process inventories are allowed,

  • The capacity of the production system is sufficient to meet the demands.

  • Zero-switch rule is used. This means that the production of each product in each cycle begins when its inventory level reaches zero.

  • The lots of products are delivered at the end of each delivery cycle.

  • Raw material of each component is independent.

  • Raw material purchase cycle for all components is constant and is determined by the solution method.

  • The order cost for all materials is equal.

  • Finished goods demands are constant and continuous.

  • The production system is flow shop, semi-finished products are passing several sequential stages and there is a single machine at each stage.

Parameters:

n: number of components

m: number of stages

i, u: component indices

j: stage index

d i : Demand rate of component i

R: Order Cost of the component

p ij : Production rate of component i at stage j.

p im : Production rate at last stage.

s ij : Sequence-independent setup time of component i at stage j.

sc i : Total setup costs of component i over all stages.

h ij : time-dependent inventory holding cost per unit of component i between stages j and j + 1

h i : Inventory holding cost (for both supplier and assembler) per unit of final component i per unit time

h i0 : Inventory holding cost per unit of raw material of component i per unit time

A: transportation cost per delivery

M: A large real number Decision variables

T: Common production cycle length (time interval between setups)

b ij : Process beginning time of component i at stage j (after related setup operation)

Since processing each component at each stage, a value is supposed to be added for the given component, the set of h ij makes an increasing series, that is
h i > h i j ; h i j > h i , j 1 i = 1 , , n j = 2 , , m 1

Decision variables:

σ j : Production sequence vector at stage j

T: Production common cycle length (time interval between setups)

Q i : Production lot size of component i at different stages (Q i = d i . T)

F: The number of production cycles over the planning horizon

b ij : Process beginning time of component i at stage j (after related setup operation)

B: Raw material common purchasing cycle length that is considered greater than or equal to common production cycle length (B = e.T; e is an integer number)

α i : The delivery point number in which the component i is delivered to the costumer.

The problem (Problem P) can be formulated as a mixed integer nonlinear programming. As mentioned earlier, to formulate this problem, the common cycle policy has been applied. The objective of Problem P is to minimize the average of transportation, setup, work-in-process and raw material and end component inventory holding costs per unit time for this simple supply chain. It goes without saying that the average setup cost per unit time is i = 1 n s c i / T [12] and the average delivery cost per unit time is A.k/T in the objective function. The inventory holding costs are somewhat more complicated. It is obvious that the inventory holding costs are incurred at all three levels in the supply chain, the raw material suppliers, manufacturer and the assembly facility. Figure 3 shows the inventory level of finished component i at the assembly facility. According to Fig. 3, the total cost of inventory of component i per unit time at the assembly facility is i = 1 n d i . T . h i 2 [11].

Fig. 3.
Fig. 3.

Finished goods inventory level in AF

Citation: International Review of Applied Sciences and Engineering 13, 2; 10.1556/1848.2021.00334

Two types of inventory are considered on the manufacturer's side: work-in-process inventory and finished goods inventory. Figures 4 and 5 show the work-in-process inventory of component i between two successive stages j–1 and j, and the inventory level of finished component i, respectively. In Fig. 4, the average WIP inventory of component i between two successive stages j–1 and j per unit time is as below [11]:
T C w i p = i = 1 n j = 2 m h i , j 1 . d i . ( b i j + d i . T 2 p i j b i , j 1 d i . T 2 p i , j 1 )
Fig. 4.
Fig. 4.

Work-in-process inventory level

Citation: International Review of Applied Sciences and Engineering 13, 2; 10.1556/1848.2021.00334

Fig. 5.
Fig. 5.

Finished goods inventory level in manufacturer's site

Citation: International Review of Applied Sciences and Engineering 13, 2; 10.1556/1848.2021.00334

The finished goods inventory level in the manufacturer is depicted in Fig. 5. The producing process of semi-finished product i is begun at b im in stage m by production rate p im . As soon as passing t im = d iT /p im period, the inventory level reaches d i .T. Finally, inventory level vanishes as soon as production lot has been delivered to the customer at the end of the cycle.

Consequently, the total inventory holding cost of finished goods in manufacturer's site is depicted in Eq. (3) which “INT” refers to the integer part of a number in this equation.
TC F I = i = 1 n h i m . d i . { d i . T 2 p i m + ( T k ( 1 + I N T ( b i , m + d T i p i m T k ) ) ( b i m + d i T p i m ) ) } = i = 1 n h i m . d i . T . { 1 k ( d i 2. p i m ) } i = 1 n h i m . d i . b i m + i = 1 n h i m . d i . T k . I N T ( ( k . b i m T ) + ( k . d i p i m ) )

By having D, T and b ij of each item, the coefficient α i as a delivery point, will be obtained through “while cycle” as below:

For i = 1 to n

α = 1

While α D [ b i m + d i . T p i m ] < 0 do

  • α = α + 1

  • End while

  • α i = α

Next i

As it is depicted in Fig. 6, one lot of raw material of item i with value of e.d i .T (that is equal to the requirement of e production cycles) is entered into the manufacturer's raw material warehouse. At t = 0, the inventory level of item i is e.d i .T and at the end of each production cycle, it decreases by d i .T. This inventory is reducing in every cycle with velocity of p i1 (production rate of the first stage in flow shop system) and finally, after e cycle, the inventory level would be zero. The average inventory of an item is calculated as below. Whereas, the first part of the formulation is the area of trapezoid and the second one is for the area of rectangle in each cycle T.
I R I = 1 e . T { v = 1 e [ ( ( e + 1 v ) d i T + ( e v ) d i T ) . d i T 2 p i 1 ] + v = 1 e [ ( e v ) d i T . ( T d i T p i 1 ) ] } d i 2 T 2 p i 1 e v = 1 e ( 2 e 2 v + 1 ) + 1 e v 1 e [ ( e v ) d i T ( 1 d i p i 1 ) ]
Fig. 6.
Fig. 6.

Raw material inventory level in manufacturer's site

Citation: International Review of Applied Sciences and Engineering 13, 2; 10.1556/1848.2021.00334

Therefore, the total cost of raw material inventory holding cost on the manufacturer's side would be:
T C R I = i = 1 n h i 0 d i 2 . T 2 p i 1 e ( v = 1 e ( 2 e 2 v + 1 ) ) + 1 e i = 1 n h i 0 v = 1 e [ ( e v ) d i T ( 1 d i p i 1 ) ]
It is supposed that the order cost of all materials is given and equal. Therefore, the total cost would be calculated by the multiplication of the order cost and the number of order as follows:
T C R M I = R e T n
In raw material supplier's side, the supplier production rate is assumed equal to the demand rate of the manufacturer. That is because the model is the classic EOQ model in which the lots are delivered in each e cycle time (e.d i .T). The inventory chart is depicted in Fig. 7 and the total cost calculation is as follows:
T C R I 2 = i = 1 n h i 0 . e . d i . T 2
Fig. 7.
Fig. 7.

Raw material inventory level in supplier's site

Citation: International Review of Applied Sciences and Engineering 13, 2; 10.1556/1848.2021.00334

Regarding all previous calculations, the mathematical model of the proposed problem and its constraints are shown in problem 1 as follows:

Problem 1
( 10 ) M i n T C = i = 1 n h i m . d i . T . { 1 k ( d i 2. p i m ) } i = 1 n h i m . d i . b i m + i = 1 n h i m . d i . T k . I N T ( ( k . b i m T ) + ( k . d i p i m ) ) + i = 1 n j = 2 m h i , j 1 . d i . ( b i j + d i . T 2 p i j b i , j 1 d i . T 2 p i , j 1 ) + i = 1 n h i 0 . d i 2 . T 2 p i 1 e . ( v = 1 e ( 2 e 2 v + 1 ) ) + 1 e i = 1 n h i 0 . v = 1 e [ ( e v ) d i T . ( 1 d i p i 1 ) ] + i = 1 n d i . T . h i 2 + R e . T . n + i = 1 n s c i T + A . k T + i = 1 n h i 0 . e . d i . T 2
Subject to:
b i , j 1 + d i . T p i , j 1 b i j i = 1 , , n j = 2 , , m
b i j + d i . T p i j + s u j b u j M ( 2 z i , l , j z u , l + 1 , j ) ; i = 1 , , n u i , l < n
l = 1 n z i l j = 1 ; i = 1 , , n
i = 1 n z i l j = 1 ; l = 1 , , n
b i j s i j . z i 1 j ; i = 1 , , n
b i m + d i . T p i m T ; i = 1 , , n
F . T = H
e T = B
F , e 1 , integer
T , B 0 , b i j 0 , z i l j { 0,1 }

Problem 1 has the following set of constraints. Constraint (9) states that no component can be processed before it is completed at the previous stage. Constraint (10) shows that at each stage, no component can be processed before the completion of its predecessor in the related production sequence. Constraints (11) and (12) are assignment constraints at stages with only one machine; state that each component has a unique position in the sequence of these stages. Constraint (13) indicates that at each stage, processing the first component in the corresponding sequence cannot be started before setting up the corresponding machine. Constraint (14) assures that the resulting schedule is cyclic in a way that the process completion time for each component at the final stage is less than or equal to cycle time. Constraint (15) implies that the common cycle is such that the planning horizon H is a multiple of T. Constraints (17) show that F and e are some integers greater than or equal to one. Finally, Constraints (18) are the non-negativity variable constraints.

To solve the developed model, two separated heuristics and meta-heuristics methods are applied. The methods are explained in the next sections.

2.1 Optimal enumeration method

Regarding literature review, the developed model is NP-hard and it is impossible to obtain a global optimum solution. This model is a mixed integer non-linear programming. Solving this problem by means of exact method on small scale, not only the global optimum is not satisfied, but also the quality of its local optimum is drastically low. In order to generate better and more reasonable outputs, an optimal enumeration method (OEM) is presented. The OEM method has been widely applied in many problems in the supply chain; see Ounniche & Boctor [21] and Torabi et al. [11, 12].

The key principle of the OEM is that the economic lot, purchase and delivery scheduling problem consists of three variables F, e, h where production, purchase and delivery cycles can be derived by them. The variable F denotes the number of production common cycles and its length (T) in finite planning horizon. Since T is a real variable, in OEM we should use the integer variable F. Regarding the high complexity of the problem, LINGO software is not able to solve this problem obviously, even on a small scale. Therefore, LINGO should be applied in OEM. In this regard, to reduce the complexity of the problem, the values of F and e are given to LINGO and converted from variables to parameters. OEM consists of two interconnected cycles as follows:

Cycle 1

F = 1, go to the cycle 2 and obtain the value of Z F

set z * = z 1 , e * = e , T * = T 1 , k * = k 1 ;

Repeat until stop condition is satisfied

F=F+1 , and go to the cycle 2 and obtain the value of Z F

if Z F z * then

z * = z F , e * = e , T * = T F , k * = k F ;

else

  • stop

end if

end

Cycle 2

Consider the value of e variable equal to 1 and we will have Z F = Z 1,1

For e = 1 to F

  • Solve the resulting MILP

If z e z F then

  • ZF = z e , TF = T , kF = k;

End

End

In cycle 1 the value of F is given at first. Then Cycle 2 finds the most convenient (least total cost) value of e with respect to F. This procedure is iterated while the value of the objective function has a decreasing trend.

It is possible that the algorithm falls into the local optimum. This is because the algorithm is not a comprehensive search and it would be terminated in the first optimum situation, even if it was a local optimum. Due to its nonlinear nature, it is highly infeasible to get global optimum. Furthermore, it is proved that the mixed-integer linear problems have an exponential complexity, known as NP-complete problems [22]. So in a large-scale problem, this algorithm has low efficiency and it is recommended to apply Meta-heuristics methodology. The next section consists of more details about Genetic Algorithm principles and its application in the developed model.

2.2 Real valued genetic algorithm

In the previous section, a new mixed-integer nonlinear programming model was proposed. It has been proved that these problems have exponential complexity. So, with larger dimensions, the solution time is increasing exponentially. The meta-heuristics methods are new and effective in solving mixed integer programming problems where all of them use an intelligent random search process to find an exact or near optimum solution. There are sufficient studies in this area and researchers have developed a GA with real-valued Genes. Right [23] developed Real-valued GA in optimization problems and compared it with the traditional binary GA method. His research shows the high efficiency of RGA as an optimization method. Muhlenbein and Voosen [24] introduced GA with continuous numbers genes and various genetic operators have been applied, as well.

2.2.1 Modified genetic algorithm structure

In this section, the details of the designed GA are explained. The principle framework of this method is Simple Genetic Algorithm (SGA) depicted in Fig. 8. However, to modify the SGA, a convergence rate is considered which decreases the probability of using mutation generation by generation based on convergence strategy.

Fig. 8.
Fig. 8.

Pseudocode of SGA

Citation: International Review of Applied Sciences and Engineering 13, 2; 10.1556/1848.2021.00334

2.2.2 Chromosome structure

In this problem, the chromosome has two divisions or sub-chromosomes. The first sub-chromosome consists of m × n genes and illustrates a sequence of production begin time of n products in all m production stages. Hence, all of them are real numbers. The other sub-chromosome has three genes. These genes demonstrate three real valued genes (e, F, k) which are purchase, production and delivery cycles respectively. Thanks to finite horizon planning in this problem, all genes are finite and have upper and lower bounds.

In the second sub-chromosome, F illustrates the number of production common cycles in finite horizon time. Since usually the planning horizon time in an industrial environment is one year, therefore in this problem the finite horizon time has considered 52 weeks. In accordance with our assumption, the maximum value of F is 52. The e variable shows the number of production cycles that would be covered by each purchase. The value of e depends on F value, so it varies between 1 to F. The third variable in this sub-chromosome is k that shows the number of delivery cycles in the production common cycle. In this problem, since the literature review illustrated that the length of production cycles usually is between 1 and 4 weeks, the maximum value of delivery cycles in the production cycle is assumed to 10. So this upper bound for k would be reasonable and changes from 1 to 10. The first sub-chromosome is a set of production beginning times of all items in all stages where genes are real numbers. These real numbers have lower and upper bounds. The former is zero and the latter is set to be T (52/F). However, the beginning times never become zero and that is due to setup time before production. Note that this constraint is applied in the developed mathematical model. Figure 9 shows the representation of chromosome with n products and m stages.

Fig. 9.
Fig. 9.

Designed chromosome scheme

Citation: International Review of Applied Sciences and Engineering 13, 2; 10.1556/1848.2021.00334

In addition, the efficiency of GA depends on the quality of the first population. That is why in this paper some simple heuristics methods are applied to find high quality vector in the first population. The first population of proposed GA consists of three groups of chromosomes. To produce the first group, the heuristics method such as LPT (Longest Process Time), SPT (Shortest Process Time) and WSPT (Weighted Shortest Process Time) have been applied. In both LPT and SPT methods, m vectors are generated. The jth vector by SPT is obtained by sorting incremental order of pt ij in jth stage of production and decreasing order in LPT. The WSPT method is similar to SPT except that the sorting criterion is based on the value of pt ij /h j .

Consequently, we have produced 3m different sequential vectors. Whereas, the desired chromosome consists of begin times (b ij ). Therefore, to convert the sequential vector to begin times, a production common cycle (T) is required. By the given value of T, production lot sizes (d i .T) would be calculated. In addition, by means of these parameters the begin time and finish time of each item in all stages are derived and a high quality of population in GA would be generated.

By having sequence through heuristics algorithm, production common cycle and random begin time of the first item in the first stage, the begin times of all items in all stages could be extracted. The rule to get all begin times is that the begin time of one lot of item i at stage j is equal to the summation of setup time of item i at stage j (s ij ) with the maximization between finish time of one lot of item i at stage j–1 and the finish time of one lot of the predecessor item (u) at stage j. This rule is as follows:
b i j = s i j + m a x { b i , j 1 + d i T p i , j 1 , b u j + d u T p u , j }

The second group of chromosomes in the first population is similar to the first group's structure. The sequence is generated randomly in the first stage and the next stages are the same. The third group might be generated to produce a more scattered population and cover more solution space. In this group, the second sub-chromosome genes are selected randomly between their ranges just like the previous groups. In addition, in the first sub-chromosome the genes are selected randomly between zero and T. Since the genes in the first sub-chromosome are produced randomly, as a result the probability of feasibility is low. To prevent transferring of these infeasible chromosomes to the next generations, their fitness value should be low. Regarding this issue, a penalty method has been developed as a part of the fitness function.

To survive the most qualified chromosomes in the next generations and get rid of low quality chromosomes in simulation of a natural process, a fitness function in GA is developed. This function evaluates the quality of each chromosome of the population. In the proposed problem, the solution with lower total cost is considered a better solution. Newton & Sarker [25] have implemented a genetic algorithm in ELSP. They utilized from penalty method to satisfy the constraints of the problem. The proposed penalty methods include Static penalty method, Dynamic penalty method and Adaptive penalty method. They executed the problem with all three methods and the experiments showed that the Static penalty method is the most efficient method in Economic Lot Scheduling problems.

Subsequently, the fitness function consists of two parts in this paper. The first part is objective function that is made of eight types of cost in the supply chain area. The second part is consist of the constraints of the problem associated with the penalty. Certainly, a chromosome that satisfies all constraints is qualified as a feasible solution. To prevent the transmission of infeasible chromosome to the next generation, static penalty method has been implemented. In the static penalty method, a penalty coefficient is assigned to each constraint. If one of the constraints is not satisfied, the absolute value of its deviation multiplied by its penalty coefficient will add to the objective function in fitness value of MRGA. The coefficients are depicted in Table 1.

Table 1.

Penalty coefficient for each constraint

constraint Penalty Coefficient constraint Penalty Coefficient
97 9 × 10 9 101 3 × 10 6
98 9 × 10 7 102 3 × 10 6

In this study, the GA algorithm is developed in MATLAB [26]. To find the most suitable crossover, several crossover operators have been examined. As depicted in Fig. 10, a simple two-point crossover is applied.

Fig. 10.
Fig. 10.

Two point crossover operator scheme

Citation: International Review of Applied Sciences and Engineering 13, 2; 10.1556/1848.2021.00334

Finally, because of the real number nature of chromosomes, the Gaussian mutation is applied. In the next step, to identify effective crossover probability (p crossover ), the GA is evaluated by various probabilities and the GA is executed 10 times for each value. As depicted in Fig. 11, the best crossover probability is equal to 0.8.

Fig. 11.
Fig. 11.

Fitness function changes based on crossover probability value

Citation: International Review of Applied Sciences and Engineering 13, 2; 10.1556/1848.2021.00334

In Gaussian mutation, a random number from Gaussian distribution is added to each chromosome to produce a new child. The Gaussian mutation in real valued chromosome is shown in Fig. 12.

Fig. 12.
Fig. 12.

Gaussian mutation

Citation: International Review of Applied Sciences and Engineering 13, 2; 10.1556/1848.2021.00334

Finally, a pre-defined stopping criterion identifies the algorithm stopping time. In the proposed Genetic Algorithm, the optimum iteration is determined equal to 200 by evaluating various iterations.

3 Result and discussion

To evaluate the efficiency of two proposed solution methods (optimal enumeration method (OEM) and Modified Real valued genetic algorithm), seven different examples with various dimensions have been analyzed. Every parameter value is extracted from uniform probability distribution function as depicted in Table 2.

Table 2.

The distributions of parameters

Parameters Description Distribution
d i Demand rate of component i U(100,1000)
p ij Production rate of component i at stage j U(1000,15000)
s ij Sequence-independent setup time of component i at stage j U(0.01,0.25)
h ij time dependent inventory holding cost per unit of component i between stages j and j + 1 U(1,20)
A transportation cost per delivery U(5000,15000)
R Order Cost of the component U(3000,6000)

All examples are run on a Celeron 2.4 GHz system. Each sample has been solved by OEM and Modified Real-valued GA 10 times. By increasing the problem dimensions, solving the problem by OEM becomes too much more time consuming. Consequently, the best obtained solutions for each problem after maximum 24-h runtime by OEM are in Table 3. It may be concluded that the OEM is not efficient specially for medium and large dimension problems. Because of the nonlinear nature of this problem, the obtained solutions by OEM method and LINGO software is considered local optimum and in general the greater dimensions of problems resulted from the lower quality of the solutions. As Fig. 13 shows, the efficiency of MRGA would be clarified by scaling up the problem.

Table 3.

Calculation result

No. Problem Dimensions (stages*items) Number of times (up to 10) that GA result is better than OEM Mean cost reduction in GA rather than OEM GA calculation time (per second)
Mean time Max. time Scatter index
1 3 * 3 7 911.155 26.76 27.225 0.155
2 2 * 6 8 2,409.895 29.212 31.452 0.746
3 3 * 5 8 2,913.495 29.663 32.198 0.845
4 3 *6 9 3,496.265 31.956 35.149 1.06
5 2 * 11 10 3,968.755 34.264 38.568 1.43
6 5 * 5 10 4,381.846 36.869 43.247 2.12
7 3 * 10 10 4,966.825 39.548 46.894 2.44
Fig. 13.
Fig. 13.

The performance of MRGA in various samples

Citation: International Review of Applied Sciences and Engineering 13, 2; 10.1556/1848.2021.00334

Regarding the empirical result, it would be deduced that:

  1. Thanks to cost reduction and mean time indexes, it is obvious that the MRGA as a metaheuristic method is more efficient than OEM as heuristic method especially on a larger scale.

  2. In case of larger scale of samples, by considering mean time and max time, the scatter index of solutions grows up exponentially. It would be concluded that by increasing the scale of the problem, the calculation time of the proposed methodology is raised and the efficiency of MRGA might be weakened.

4 Conclusions

Due to the recent development in supply chain concepts and the importance of suitable decision models in order to obtain high-level coordination between supply chain partners and cost reduction in the supply chain, the Lot-sizing and Scheduling Problems are widely prioritized in this area.

The literature review shows that there are three branches of using Lot sizing and Scheduling Models. The first branch is Economic Lot and Scheduling Problem (ELSP). In this problem, a supplier with a cyclic production system produces several parts at a constant rate and delivers them to customers continuously. In this category, optimization of supplier's production system is the main goal and the whole supply chain is neglected. The second branch is related to Economic Lot, Delivery and Scheduling Problem (ELDSP). In these problems, a supplier produces various parts in a cyclic production system and delivers the production lot to the customer at the end of each cycle. In this model, the optimization of supplier and customer is considered simultaneously and cost reduction in two-level supply chain is desirable.

The last branch that is considered by researchers recently, is Economic lot and scheduling problem in three-echelon supply chain which consists of customer (Assembly Facility), manufacturer and raw material suppliers. The manufacturer works with single-stage system in an infinite horizon and the production lot is delivered to the customer at the end of each cycle. In all the described models, the manufacturer's production system was considered single-stage. However, in this paper this restriction is relaxed and the system is set to be multi stage (flow shop system). In the proposed Supply Chain; Assembly Facility purchases the raw material from their suppliers in optimized lot, changes them to finished goods and finally delivers them to customers in multi delivery points during the production cycle in an optimum quantity.

In this paper, we developed a new mixed integer nonlinear programming model in accordance with the described scenario. After that, two solution methods were proposed, OEM and MRGA. Some empirical experiments were done and two proposed methods were compared together and finally, the calculation showed that MRGA is more effective than OEM in time-consumption and solution quality.

Regarding future research, there are various fields to improve decision-making in economic lot and scheduling in the supply chain. As far as the structure of the supply chain is concerned, the multi tiers of suppliers could be considered and optimization of production lots and scheduling in a comprehensive supply chain (from tier n suppliers to end users) would be desired to be more compatible with the real world. Another area for improvement is to develop a mathematical model and consider more limitations such as dependent setup times, backlogs, physical inventory limits, etc. It would be worthwhile to include costs related to the value of stocks in the model such as depreciation losses on inventories, insurance costs, etc. Improvement more effective optimizing methods is also considerable.

References

  • [1]

    J. Rogers , “A computational approach to the economic lot scheduling problem,” Manage. Sci., vol. 4, no. 1, pp. 264291, 1958.

  • [2]

    W. Hsu , “On the general feasibility test of scheduling lot sizes for several products on one machine,” Manage. Sci., vol. 29, no. 1, pp. 93105, 1983.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [3]

    A. Santander-Mercado and M. Jubiz-Diaz , “The economic lot scheduling problem: a survey,” Int. J. Prod. Res., vol. 54, no. 16, pp. 49734992, 2016.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [4]

    J. Y. Huang and M. J. Yao , “A new algorithm for determining production schedules when solving the multi-product economic lot sizing problem in flow shops,” Int. J. Adv. Manuf. Technol., vol. 33, pp. 560578, 2007.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [5]

    J. Ouenniche and FF Boctor , “The two-group heuristic, to solve the multi-product, economic lot-sizing and scheduling problem in flow shops,” Eur. J. Oper. Res., vol. 129, no. 3, pp. 539554, 2001.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [6]

    F. G. Beck , K. Biel , and C. H. Glock , “Integration of energy aspects into the economic lot scheduling problem,” Int. J. Prod. Econ., vol. 209, pp. 399410, 2019.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [7]

    I. Ferretti , “Multiproduct economic lot scheduling problem with returns and sorting line,” Systems, vol. 8, no. 16, pp. 120, 2020.

  • [8]

    M. Alaghebandha , B. Naderi , and M. Mohammadi , “Economic lot sizing and scheduling in distributed permutation flow shops,” J. Optim. Ind. Eng., vol. 12, no. 1, pp. 103117, 2018.

    • Search Google Scholar
    • Export Citation
  • [9]

    Hahm and Yano , “The economic lot and delivery scheduling problem: models for nested schedules,” IIE Trans., vol. 27, no. 2, pp. 126139, 2007.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [10]

    Hahm and Yano , “The economic lot and delivery scheduling problem: the single item case,” Int. J. Prod. Econ., vol. 28, no. 2, pp. 235252, 1992.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [11]

    S. A. Torabi and S. M. T. Fatemi Ghomi , “A hybrid genetic algorithm for the infinite horizon economic lot and delivery scheduling in supply chains,” Eur. J. Oper. Res., vol. 173, no. 1, pp. 173189, 2006.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [12]

    S. A. Torabi and S. M. T. Fatemi Ghomi , “The common cycle economic lot scheduling in flexible job shops: The finite horizon case,” Int. J. Prod. Econ., vol. 97, no. 1, pp. 5265, 2005.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [13]

    S. A. Torabi and M. Jenabi , “Multiple cycle economic lot and delivery-scheduling problem in a two-echelon supply chain,” Int. J. Adv. Manuf. Technol., vol. 43, no. 7, pp. 785798, 2009.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [14]

    W. Lee , “A joint economic lot size model for raw material ordering, manufacturing, and finished goods delivering,” Omega Int. J. Manag. Sci., vol. 3, no. 2, pp. 163174, 2005.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [15]

    T. Kim , Y. Hong , and S.Y. Chang , “Joint economic procurement-production-delivery policy for multipule items in a single-manufacturer, multiple retailer system,” Int. J. Prod. Econ., vol. 103, no. 1, pp. 199208, 2006.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [16]

    B. R. Sarker and A. Diponegoro , “Optimal production plans and shipment schedules in a supply-chain system with multiple suppliers and multiple buyers,” Eur. J. Oper. Res., vol. 194, no. 3, pp. 753773, 2009.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [17]

    M. Y. Jaber and S. K. Goyal , “Coordinating a three-level supply chain with multiple suppliers, a vendor and multiple buyers,” Int. J. Prod. Econ., vol. 116, no. 1, pp. 95103, 2008.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [18]

    H. Osman and K. Demirli , “Economic lot and delivery scheduling problem for multi-stage supply chains,” Int. J. Prod. Econ., vol. 36, no. 15, pp. 275286, 2012.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [19]

    A. Goli and M. Alinaghian , “A new mathematical model for production and delivery scheduling problem with common cycle in a supply chain with open-shop system,” Int. J. Manuf. Technol. Manage., vol. 34, no. 2, pp. 174187, 2020.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [20]

    A. O. El-meehy , A. K. El-Kharbotly , and M. M. El-Beheiry , “Joint lot sizing and scheduling of a multi-period supply chain,” Indep. J. Manage. Prod., vol. 10, no. 5, pp. 15161535, 2019.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [21]

    J. Ouenniche and F. F. Boctor , “The multi-product , economic lot-sizing problem in flow shops : the powers- of- two heuristic,” Comput. Operations Res., vol. 28, no. 12, pp. 11651182, 2001.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [22]

    S. Vavasis , Nonlinear Optimization: Complexity Issues. New York, USA: Oxford University Press, 1991.

  • [23]

    A. H. Wright , “Genetic algorithms for real parameters optimization,” Found. Genet. Algorithms, vol. 1, no. 1, pp. 205218, 1991.

    • Search Google Scholar
    • Export Citation
  • [24]

    H. Muhlenbein and D. S. Voosen , “Predictive models for the breeder genetic algorithm,” Evol. Comput., vol. 1, no. 1, pp. 2549, 1993.

  • [25]

    R. Sarker , “A genetic algorithm for solving economic lot size scheduling problem,” Comput. Ind. Eng., vol. 42, no. 2, pp. 189198, 2002.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [26]

    A. Chipperfield , P. Fleming , and H. Pohlheim , “Genetic algorithm TOOLBOX for use with MATLAB, user’s guide,” in International Conference on System Engineering, Coventry, UK, Sep. 6–8, 1994.

    • Search Google Scholar
    • Export Citation
  • [1]

    J. Rogers , “A computational approach to the economic lot scheduling problem,” Manage. Sci., vol. 4, no. 1, pp. 264291, 1958.

  • [2]

    W. Hsu , “On the general feasibility test of scheduling lot sizes for several products on one machine,” Manage. Sci., vol. 29, no. 1, pp. 93105, 1983.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [3]

    A. Santander-Mercado and M. Jubiz-Diaz , “The economic lot scheduling problem: a survey,” Int. J. Prod. Res., vol. 54, no. 16, pp. 49734992, 2016.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [4]

    J. Y. Huang and M. J. Yao , “A new algorithm for determining production schedules when solving the multi-product economic lot sizing problem in flow shops,” Int. J. Adv. Manuf. Technol., vol. 33, pp. 560578, 2007.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [5]

    J. Ouenniche and FF Boctor , “The two-group heuristic, to solve the multi-product, economic lot-sizing and scheduling problem in flow shops,” Eur. J. Oper. Res., vol. 129, no. 3, pp. 539554, 2001.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [6]

    F. G. Beck , K. Biel , and C. H. Glock , “Integration of energy aspects into the economic lot scheduling problem,” Int. J. Prod. Econ., vol. 209, pp. 399410, 2019.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [7]

    I. Ferretti , “Multiproduct economic lot scheduling problem with returns and sorting line,” Systems, vol. 8, no. 16, pp. 120, 2020.

  • [8]

    M. Alaghebandha , B. Naderi , and M. Mohammadi , “Economic lot sizing and scheduling in distributed permutation flow shops,” J. Optim. Ind. Eng., vol. 12, no. 1, pp. 103117, 2018.

    • Search Google Scholar
    • Export Citation
  • [9]

    Hahm and Yano , “The economic lot and delivery scheduling problem: models for nested schedules,” IIE Trans., vol. 27, no. 2, pp. 126139, 2007.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [10]

    Hahm and Yano , “The economic lot and delivery scheduling problem: the single item case,” Int. J. Prod. Econ., vol. 28, no. 2, pp. 235252, 1992.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [11]

    S. A. Torabi and S. M. T. Fatemi Ghomi , “A hybrid genetic algorithm for the infinite horizon economic lot and delivery scheduling in supply chains,” Eur. J. Oper. Res., vol. 173, no. 1, pp. 173189, 2006.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [12]

    S. A. Torabi and S. M. T. Fatemi Ghomi , “The common cycle economic lot scheduling in flexible job shops: The finite horizon case,” Int. J. Prod. Econ., vol. 97, no. 1, pp. 5265, 2005.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [13]

    S. A. Torabi and M. Jenabi , “Multiple cycle economic lot and delivery-scheduling problem in a two-echelon supply chain,” Int. J. Adv. Manuf. Technol., vol. 43, no. 7, pp. 785798, 2009.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [14]

    W. Lee , “A joint economic lot size model for raw material ordering, manufacturing, and finished goods delivering,” Omega Int. J. Manag. Sci., vol. 3, no. 2, pp. 163174, 2005.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [15]

    T. Kim , Y. Hong , and S.Y. Chang , “Joint economic procurement-production-delivery policy for multipule items in a single-manufacturer, multiple retailer system,” Int. J. Prod. Econ., vol. 103, no. 1, pp. 199208, 2006.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [16]

    B. R. Sarker and A. Diponegoro , “Optimal production plans and shipment schedules in a supply-chain system with multiple suppliers and multiple buyers,” Eur. J. Oper. Res., vol. 194, no. 3, pp. 753773, 2009.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [17]

    M. Y. Jaber and S. K. Goyal , “Coordinating a three-level supply chain with multiple suppliers, a vendor and multiple buyers,” Int. J. Prod. Econ., vol. 116, no. 1, pp. 95103, 2008.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [18]

    H. Osman and K. Demirli , “Economic lot and delivery scheduling problem for multi-stage supply chains,” Int. J. Prod. Econ., vol. 36, no. 15, pp. 275286, 2012.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [19]

    A. Goli and M. Alinaghian , “A new mathematical model for production and delivery scheduling problem with common cycle in a supply chain with open-shop system,” Int. J. Manuf. Technol. Manage., vol. 34, no. 2, pp. 174187, 2020.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [20]

    A. O. El-meehy , A. K. El-Kharbotly , and M. M. El-Beheiry , “Joint lot sizing and scheduling of a multi-period supply chain,” Indep. J. Manage. Prod., vol. 10, no. 5, pp. 15161535, 2019.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [21]

    J. Ouenniche and F. F. Boctor , “The multi-product , economic lot-sizing problem in flow shops : the powers- of- two heuristic,” Comput. Operations Res., vol. 28, no. 12, pp. 11651182, 2001.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [22]

    S. Vavasis , Nonlinear Optimization: Complexity Issues. New York, USA: Oxford University Press, 1991.

  • [23]

    A. H. Wright , “Genetic algorithms for real parameters optimization,” Found. Genet. Algorithms, vol. 1, no. 1, pp. 205218, 1991.

    • Search Google Scholar
    • Export Citation
  • [24]

    H. Muhlenbein and D. S. Voosen , “Predictive models for the breeder genetic algorithm,” Evol. Comput., vol. 1, no. 1, pp. 2549, 1993.

  • [25]

    R. Sarker , “A genetic algorithm for solving economic lot size scheduling problem,” Comput. Ind. Eng., vol. 42, no. 2, pp. 189198, 2002.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • [26]

    A. Chipperfield , P. Fleming , and H. Pohlheim , “Genetic algorithm TOOLBOX for use with MATLAB, user’s guide,” in International Conference on System Engineering, Coventry, UK, Sep. 6–8, 1994.

    • Search Google Scholar
    • Export Citation
  • Collapse
  • Expand

Senior editors

Editor-in-Chief: Ákos, LakatosUniversity of Debrecen, Hungary

Founder, former Editor-in-Chief (2011-2020): Ferenc Kalmár, University of Debrecen, Hungary

Founding Editor: György Csomós, University of Debrecen, Hungary

Associate Editor: Derek Clements Croome, University of Reading, UK

Associate Editor: Dezső Beke, University of Debrecen, Hungary

Editorial Board

  • Mohammad Nazir AHMAD, Institute of Visual Informatics, Universiti Kebangsaan Malaysia, Malaysia

    Murat BAKIROV, Center for Materials and Lifetime Management Ltd., Moscow, Russia

    Nicolae BALC, Technical University of Cluj-Napoca, Cluj-Napoca, Romania

    Umberto BERARDI, Toronto Metropolitan University, Toronto, Canada

    Ildikó BODNÁR, University of Debrecen, Debrecen, Hungary

    Sándor BODZÁS, University of Debrecen, Debrecen, Hungary

    Fatih Mehmet BOTSALI, Selçuk University, Konya, Turkey

    Samuel BRUNNER, Empa Swiss Federal Laboratories for Materials Science and Technology, Dübendorf, Switzerland

    István BUDAI, University of Debrecen, Debrecen, Hungary

    Constantin BUNGAU, University of Oradea, Oradea, Romania

    Shanshan CAI, Huazhong University of Science and Technology, Wuhan, China

    Michele De CARLI, University of Padua, Padua, Italy

    Robert CERNY, Czech Technical University in Prague, Prague, Czech Republic

    Erdem CUCE, Recep Tayyip Erdogan University, Rize, Turkey

    György CSOMÓS, University of Debrecen, Debrecen, Hungary

    Tamás CSOKNYAI, Budapest University of Technology and Economics, Budapest, Hungary

    Anna FORMICA, IASI National Research Council, Rome, Italy

    Alexandru GACSADI, University of Oradea, Oradea, Romania

    Eugen Ioan GERGELY, University of Oradea, Oradea, Romania

    Janez GRUM, University of Ljubljana, Ljubljana, Slovenia

    Géza HUSI, University of Debrecen, Debrecen, Hungary

    Ghaleb A. HUSSEINI, American University of Sharjah, Sharjah, United Arab Emirates

    Nikolay IVANOV, Peter the Great St. Petersburg Polytechnic University, St. Petersburg, Russia

    Antal JÁRAI, Eötvös Loránd University, Budapest, Hungary

    Gudni JÓHANNESSON, The National Energy Authority of Iceland, Reykjavik, Iceland

    László KAJTÁR, Budapest University of Technology and Economics, Budapest, Hungary

    Ferenc KALMÁR, University of Debrecen, Debrecen, Hungary

    Tünde KALMÁR, University of Debrecen, Debrecen, Hungary

    Milos KALOUSEK, Brno University of Technology, Brno, Czech Republik

    Jan KOCI, Czech Technical University in Prague, Prague, Czech Republic

    Vaclav KOCI, Czech Technical University in Prague, Prague, Czech Republic

    Imre KOCSIS, University of Debrecen, Debrecen, Hungary

    Imre KOVÁCS, University of Debrecen, Debrecen, Hungary

    Angela Daniela LA ROSA, Norwegian University of Science and Technology, Trondheim, Norway

    Éva LOVRA, Univeqrsity of Debrecen, Debrecen, Hungary

    Elena LUCCHI, Eurac Research, Institute for Renewable Energy, Bolzano, Italy

    Tamás MANKOVITS, University of Debrecen, Debrecen, Hungary

    Igor MEDVED, Slovak Technical University in Bratislava, Bratislava, Slovakia

    Ligia MOGA, Technical University of Cluj-Napoca, Cluj-Napoca, Romania

    Marco MOLINARI, Royal Institute of Technology, Stockholm, Sweden

    Henrieta MORAVCIKOVA, Slovak Academy of Sciences, Bratislava, Slovakia

    Phalguni MUKHOPHADYAYA, University of Victoria, Victoria, Canada

    Balázs NAGY, Budapest University of Technology and Economics, Budapest, Hungary

    Husam S. NAJM, Rutgers University, New Brunswick, USA

    Jozsef NYERS, Subotica Tech College of Applied Sciences, Subotica, Serbia

    Bjarne W. OLESEN, Technical University of Denmark, Lyngby, Denmark

    Stefan ONIGA, North University of Baia Mare, Baia Mare, Romania

    Joaquim Norberto PIRES, Universidade de Coimbra, Coimbra, Portugal

    László POKORÁDI, Óbuda University, Budapest, Hungary

    Roman RABENSEIFER, Slovak University of Technology in Bratislava, Bratislava, Slovak Republik

    Mohammad H. A. SALAH, Hashemite University, Zarqua, Jordan

    Dietrich SCHMIDT, Fraunhofer Institute for Wind Energy and Energy System Technology IWES, Kassel, Germany

    Lorand SZABÓ, Technical University of Cluj-Napoca, Cluj-Napoca, Romania

    Csaba SZÁSZ, Technical University of Cluj-Napoca, Cluj-Napoca, Romania

    Ioan SZÁVA, Transylvania University of Brasov, Brasov, Romania

    Péter SZEMES, University of Debrecen, Debrecen, Hungary

    Edit SZŰCS, University of Debrecen, Debrecen, Hungary

    Radu TARCA, University of Oradea, Oradea, Romania

    Zsolt TIBA, University of Debrecen, Debrecen, Hungary

    László TÓTH, University of Debrecen, Debrecen, Hungary

    László TÖRÖK, University of Debrecen, Debrecen, Hungary

    Anton TRNIK, Constantine the Philosopher University in Nitra, Nitra, Slovakia

    Ibrahim UZMAY, Erciyes University, Kayseri, Turkey

    Andrea VALLATI, Sapienza University, Rome, Italy

    Tibor VESSELÉNYI, University of Oradea, Oradea, Romania

    Nalinaksh S. VYAS, Indian Institute of Technology, Kanpur, India

    Deborah WHITE, The University of Adelaide, Adelaide, Australia

International Review of Applied Sciences and Engineering
Address of the institute: Faculty of Engineering, University of Debrecen
H-4028 Debrecen, Ótemető u. 2-4. Hungary
Email: irase@eng.unideb.hu

Indexing and Abstracting Services:

  • DOAJ
  • ERIH PLUS
  • Google Scholar
  • ProQuest
  • SCOPUS
  • Ulrich's Periodicals Directory

 

2023  
Scimago  
Scimago
H-index
11
Scimago
Journal Rank
0.249
Scimago Quartile Score Architecture (Q2)
Engineering (miscellaneous) (Q3)
Environmental Engineering (Q3)
Information Systems (Q4)
Management Science and Operations Research (Q4)
Materials Science (miscellaneous) (Q3)
Scopus  
Scopus
Cite Score
2.3
Scopus
CIte Score Rank
Architecture (Q1)
General Engineering (Q2)
Materials Science (miscellaneous) (Q3)
Environmental Engineering (Q3)
Management Science and Operations Research (Q3)
Information Systems (Q3)
 
Scopus
SNIP
0.751


International Review of Applied Sciences and Engineering
Publication Model Gold Open Access
Online only
Submission Fee none
Article Processing Charge 1100 EUR/article
Regional discounts on country of the funding agency World Bank Lower-middle-income economies: 50%
World Bank Low-income economies: 100%
Further Discounts Limited number of full waivers available. Editorial Board / Advisory Board members: 50%
Corresponding authors, affiliated to an EISZ member institution subscribing to the journal package of Akadémiai Kiadó: 100%
Subscription Information Gold Open Access

International Review of Applied Sciences and Engineering
Language English
Size A4
Year of
Foundation
2010
Volumes
per Year
1
Issues
per Year
3
Founder Debreceni Egyetem
Founder's
Address
H-4032 Debrecen, Hungary Egyetem tér 1
Publisher Akadémiai Kiadó
Publisher's
Address
H-1117 Budapest, Hungary 1516 Budapest, PO Box 245.
Responsible
Publisher
Chief Executive Officer, Akadémiai Kiadó
ISSN 2062-0810 (Print)
ISSN 2063-4269 (Online)

Monthly Content Usage

Abstract Views Full Text Views PDF Downloads
Aug 2024 0 31 6
Sep 2024 0 25 5
Oct 2024 0 50 6
Nov 2024 0 8 2
Dec 2024 0 11 10
Jan 2025 0 11 1
Feb 2025 0 5 2