# Developing a Model with Ergonomic Aspects Using Endurance Time and Rest Allowance for Supporting the Optimization of Production Line Material Supply: A Case of Single-Operator Multi-Materials

Authors: Sezen Korkulu; Krisztián Bóna; Tamás Péter
Journal: Mathematical Problems in Engineering
 (2021)

---

## Abstract

Various studies on the issues regarding economic and environmental sustainability in inventory management have been investigated by many researchers in recent years. However, the integration of ergonomics as an aspect of social sustainability in an inventory model is still scarce. This paper presents an extension of ergonomic inventory modeling with a new mathematical model which integrates relaxation allowance and endurance time to the ergonomic inventory model. The rest allowance is determined by endurance time, contraction time, and relative force based on several items such as weight, and it helps to prevent ergonomic risks and fatigue of the back muscles. The model has been analyzed in a numerical work based on our specific simulator and parameteric analysis to present validness in different scenarios. Consequently, our results show that optimization with ergonomics in production line material supply provides notable advantages as it increases overall productivity.


---

## Body

## 1. Introduction

Inventory management and lot sizing are critical problems for many industries. The primary focus of the problems is economic sustainability. Especially, the supply of the production lines, the service time, and the availability of the materials are the most critical factors which may cause several problems. Fast service and continuous availability are the most common terms used in the production process to determine the aims of the main functions for optimization problems. In every case, the operators at the line must be able to reach the necessary materials for production. Otherwise, the production line may be stopped, and consequently, there may be a considerable cost and an impracticable problem. This problem can be handled by holding a lot of stocks at the lines, but this is not an economically sustainable solution as the specific cost of the reserved production area will be expensive. The other solution can be to supply the production lines with higher frequency and to hold a lower level of stocks at the lines. In this case, the total cost of the operations will increase. Many researchers have examined these problems extensively by beginning with the Economic Order Quantity (EOQ) model developed by Harris [1]. The EOQ and lot-sizing models have been studied and extended by many researchers in accordance with economic (time-varying demand, inflation, shortages, and backlogging), environmental (CO2 emission and environmental cost), and social (human rights, fair labour practices, living conditions, health, and safety) sustainability aspects [2]. Especially, the researchers have focused on the environmental and economic sustainability aspects. However, the social sustainability aspect has been neglected in inventory management, and there are very few studies focused on the social dimension of sustainability [3, 4]. Although there are many studies regarding the environmental and economic sustainability of inventory models, according to Kazemi et al. [5], there are few exceptions which investigate ergonomics as a social sustainability aspect [6–19]. Even though there are studies that integrate ergonomic aspects in inventory management, there is a lack of studies that investigate the maximum endurance time, rest allowance, and production line supply with multi-materials. It is necessary to integrate social sustainability in inventory management to optimize the decision-making process and performance. Ergonomics is an important component of the decision-making process in inventory management and lot sizing. Furthermore, ergonomics is necessary for the prevention of repetitive strain injuries and work-related musculoskeletal disorders, which may cause a long-term disability. Work-related musculoskeletal disorders (WRMDs) are the most common occupational disorders, such as sprains, strains, soreness, back disorders, pain, and carpal tunnel syndrome. Especially, work-related back disorders are associated with repetitive motion, sustained postures, excessive force, the weight of object lifted, high static muscle load, high force exertion on the hands and wrists, sudden application of force, inadequate rest or recovery periods, localized mechanical stresses to tissues, and poor spinal support [20]. We suppose that the causes are preventable with an ergonomic lot sizing in the material handling, which has an important impact on the inventories at the production lines. The official statistical reports show that there is a link between manual material handling and WRMDs. In accordance with the report issued by OSHA [21], 30% of European workers suffer from back pain, which is one of the most common work-related musculoskeletal disorders, and awkward postures, as well as lifting and handling loads, are the most common reasons. According to Driscoll et al. [22], low back pain as an occupational disability affects almost 26% of the world’s population. The report issued by the European Agency for Safety and Health at Work has emphasized that work-related musculoskeletal disorders are one of the most important reasons for work-related mortality in EU-28 (14.66%) and in the world (14.96%) [23]. It leads to the total cost of work-related musculoskeletal disorders up to 2% of GDP in Europe [24].In parallel with the purpose of this paper, the questions to be considered are as follows:(1)
How to integrate maximum endurance time and rest allowance in the calculation of optimal lot sizes for the production line supply with multi-materials to prevent ergonomic risks?(2)
What is the difference between optimal solution and nonoptimal solution after the integration of the ergonomic aspects?(3)
How do different parameters (height of placement, unit weight, and distance of movement) affect the optimal solution?This study aims to contribute a new model which extends the ergonomics studies in inventory management [10, 13] and integrates ergonomic aspects using endurance time and rest allowance for supporting the optimization of production line material supply of multi-materials. This paper integrates, for the first time, the maximum endurance time and rest allowance formulas [25, 26] in a lot-sizing model as an ergonomic measure. The rest allowance prevents disorders associated with manual material handling works in the production supply process that includes picking, storing, and pushing motions. The model determines the optimal lot sizes for multi-materials. The optimal lot size represents the load value which allows workers to work continuously without the need for rest allowance, and consequently, the total cost is minimized. The lot size values higher than the optimal lot size values mean that workers need to rest to prevent any ergonomic disorders and that there may be a rest time cost as they do not work during the rest. The total cost consists of rest time cost, holding cost, and material handling cost. The optimal solution of the model is based on the minimization of both economic and ergonomic costs. The main innovations considered in this study can be summarized as follows. A new modeling approach, which integrates endurance time and rest allowance to support the multi-objective optimization of the production line material supply process for single operator and multi-materials, has been developed, and a numerical study and parameteric analysis have been conducted to investigate the opportunities of the proposed model. Furthermore, the rest allowance [25] and maximum endurance time [26] formulations are linearized.This paper has been structured as follows. The following section will present a literature review on sustainable inventory management and ergonomic measures. In Section3, the problem will be examined in detail, and necessary assumptions will be explained. In Section 4, the mathematical formulation of the lot-sizing model developed will be shown with ergonomic aspects based on single operator and multi-materials. In Section 5, the results of a numerical example and parameteric analysis will be presented, and Section 6 will include the conclusion of the paper.


## 2. Literature Review

Sustainable inventory management and lot sizing aim to reduce social and environmental effects of an industry without affecting profitability [27]. There are many studies focused on the environmental and economic aspects of sustainability of inventories. Taft [28] has extended the EOQ model with continuous consumption and periodic production and has developed the Economic Production Quantity (EPQ) model. Wagner and Whitin [29] developed the Dynamic Economic Lot (DEL) model where the demand rate is supposed to be changing over a set number of periods. The EOQ model has also been extended for shortages, deterioration of items, payment delays, and growing items [30–33]. Especially, the environmental and economic sustainability aspects of inventory management have been studied. Inventory models have been developed with various economic aspects, such as imperfect production cost [34–36], incremental discounts [37, 38], and transportation cost [39–41]. The environmental sustainability aspect embodies an extensive investigation on CO2 emission and green emission [5, 27, 42–47]. For instance, an EPQ model with an environmental quality index, which considers economic and environmental aspects, was developed by Glock et al. [48]. Furthermore, an EOQ model, which considers CO2 emission as an environmental cost, was developed by Wahab et al. [49].In the literature, there are many studies on inventory models regarding economic and environmental sustainability aspects [34, 40, 42–45]. However, there are only a few contributions regarding human health and ergonomic aspects as social sustainability aspects [50].The ability of a muscle to maintain its function over a period of time with multiple contractions is endurance [51]. The maximum endurance time of the muscles is an important indicator for measuring muscular fatigue to prevent work-related musculoskeletal disorders. Many researchers have conducted several studies on endurance time. Rohmert [25] developed a general endurance time model for all muscle groups with a limit of 15% relative force, which means that, under this value, workers can continue to work without any fatigue. Monod and Scherrer [52] focused on the notions of maximum work and maximum time of work to define the critical power of muscular work. In addition, they developed an endurance time model for static, dynamic, and intermittent static work to define the maximum amount of work that can be performed in a specific time as well as the circumstances of work performed without fatigue. Sato et al. [53] investigated the relationship between the force level of isometric contractions and endurance time on the elbow flexors, shoulder abductors, and knee extensors. Rohmert [54] developed endurance time models for muscle groups, such as shoulder, elbow, trunk, and hand. Sjogaard [55] investigated endurance time and intramuscular changes during long-term contractions. The author has found that rest periods are needed for the recovery of muscles even after low-load muscular work. The relationship between endurance time and relative force for the static pull and static torque was investigated by Manenica [56]. Rose et al. [57] investigated endurance time in fully flexed postures. They found that the fully flexed postures may be assessed by more general prediction models for endurance. In addition, energy expenditure is one of the other ergonomic measurement methods used to estimate the required energy to complete a manual task from oxygen consumption [58]. The Ovako Working Posture Assessment System (OWAS), which analyzes working postures, has identified 72 different postures [59]. The NIOSH lifting equation, which assesses the ergonomic risks, has been used to evaluate lifting and lowering tasks [60]. The maximum endurance time has been defined as a key parameter to prevent muscle fatigue [61]. The maximum endurance time models have primarily investigated different muscle groups, such as elbow, hand, and shoulder or postures. As the aim of this study is to examine the load on the spine, various endurance time studies, which examine the muscle groups of the elbow, hand, and shoulder as well as different postures, are not suitable for our research. For this reason, the general endurance model developed by Rohmert [25] has been used in this study as it is suitable for all muscle groups and postures.Table1 shows the sustainable inventory models with ergonomic aspects in the literature. We have tried to summarize the main contributions based on the “number of objectives,” such as the “bi-objective” (A) and the “multi-objective” (B) aspects, and other important “ergonomic aspects,” such as the application of “human energy expenditure” (C), the “rest allowance” (D), the “energy expenditure rate based on Predetermined Motion Energy System” (E), the “lifting index” (F), the “vibration exposure” (G), the “OCRA (the occupational repetitive actions)” (H), the “maximum endurance time” (I), the “environmental ergonomics” (J), the “assembly system design with ergonomics” (K), and, last but not least, the “human availability” (L).Table 1
Sustainable inventory models with ergonomic aspects in the literature.
StudiesNumber of objectivesErgonomic aspectsABCDEFGHIJKLBattini et al. [8]●●Battini et al. [9]●●Battini et al. [10]●●●Battini et al. [11]●●Andriolo et al. [7]●●Battini et al. [12]●●●Battini et al. [13]●●●Botti et al. [14]●●●Zadjafar and Gholamian [18]●●Finco et al. [15]●●●Finco et al. [16]●●●Al-Araidah et al. [6]●●Finco et al. [17]●●●●Zhang et al. [19]●●This study●●●Battini et al. [8] analyzed the relationship between ergonomics and assembly system design techniques. They examined ergonomics in designing the assembly system and found that the design of the assembly system and ergonomics should complete each other for assessment of time efficiency in modern industries. In addition, they did not consider the maximum endurance time and rest allowance for assessment of the ergonomic risks. Battini et al. [10] developed a new measurement technique with consideration of energy expenditure equations [58] and rest allowance formulation [26]. Their method allows us to simplify the ergonomic assessment of each assembly task. Furthermore, they did not consider the determination of the optimal lot size and did not cover the multi-material supply. Battini et al. [9] developed various functions which are considered warehouse picking activities with the human availability and the rest allowance. They found that an improvement in ergonomic conditions has a positive impact on the total cost of the system. Their model did not cover motions other than picking associated with the WRMDs, and it did not consider the lot size. Battini et al. [11] developed a new multi-objective model for assembly line balancing, which includes energy expenditure rate based on the Predetermined Motion Energy System (PMES), and their model did not consider the lot size and multi-material supply. Battini et al. [12] developed a mixed-integer model which integrates assembly line balancing and parts feeding by incorporating ergonomic aspects based on rest allowance formulas of Garg et al. [58] and Price [62]. Their model did not consider the lot size and multi-material supply. Andriolo et al. [7] developed a lot-sizing model that considers multi-objective optimization of ergonomic aspects based on a lifting index (LI). Their model aims to minimize both cost and ergonomic risks, and it did not consider the rest allowance concept but instead used Pareto frontiers. In addition, their model is bi-objective, and they did not consider different motions which cause WRMDs and multi-material supply. Battini et al. [13] developed a mathematical model which investigates picking and storing motions with the energy expenditure rate for rest time assessment. Their ergonomic lot size model is based on the rest allowance formulation developed by Price [62]. Their model did not consider the multi-material supply. In addition, they have not considered motions other than picking and storing, which cause the WRMDs. Botti et al. [14] developed a mathematical model to design lean processes with ergonomics for hybrid assembly lines. They did not consider the lot size and manual material handling motions that cause WRMDs. Finco et al. [15] developed three heuristic methods that consider energy expenditure and rest allowance for assembly balancing problems. In addition, their model did not cover the determination of the lot size. Zadjafar and Gholamian [18] developed a mathematical model that considers environmental ergonomics and environmental pollution to minimize both cost and environmental pollution. Finco et al. [16] developed a bi-objective model that considers vibration exposure to eliminate ergonomic risks in assembly line design. Their model did not cover the multi-material supply and the determination of the optimal lot size. Al-Araidah et al. [6] developed a Monte Carlo simulation to estimate fatigue allowance [62] with the integration of energy expenditure formulas of Garg et al. [58] for female order pickers. Their study did not cover the manual material handling motions other than picking and did not consider the optimal lot size. Finco et al. [17] developed a mathematical model that integrates human energy expenditure as an ergonomic aspect in assembly line balancing problems with the smoothness index. However, they have did not consider the lot size and multi-material supply. Zhang et al. [19] developed a multi-objective mathematical model that integrates the OCRA method to eliminate the ergonomic risks for U-shaped assembly lines, and their model did not consider the determination of the optimal lot size.As can be seen in Table1, the research gap can be described as follows. Many kinds of research have focused on human energy expenditure as an ergonomic aspect of socially sustainable inventory management, but the maximum endurance time as an ergonomic measure has not been studied. Moreover, the issue of handling multi-materials in the production line material supply with ergonomic measures has not been investigated. Then again, pushing motion has been neglected in previous studies of sustainable inventory models which integrate ergonomic aspects. There are no previous studies on sustainable inventory models considering the maximum endurance time as an ergonomic measure and multi-material handling of production line material supply. There is a gap in the literature regarding studies that consider maximum endurance time and rest allowance, multi-material production material supply, and motions such as pushing or pulling as ergonomic aspects in the modeling. Furthermore, most of the inventory models apply special costs to the operations for refilling the inventory system with materials. Such refilling activities, especially for the supply of the production lines, are carried out in connection with the human workforce. However, the relevant activity carried out contains only the cost factors and does not contain the effects of ergonomic factors for the human workforce. Therefore, this paper presents the extension of the ergonomics studies in inventory management [10, 13] and the development of a new model for the determination of the lot sizes in material handling, especially in the case of the line supply with ergonomic aspects, to close this research gap. For the first time, the maximum endurance time has been integrated into a sustainable inventory model as an ergonomic measure. According to the literature, there is a lack of studies that consider pushing motions and multi-material production line supply. For this reason, the developed single-operator-multi-material inventory model considers rest allowance, maximum endurance time functions developed by Rohmert [25, 26], and picking, storing, and pushing motions for preventing work-related back disorders and ergonomic risks and supports economic and social sustainability.


## 3. Problem Description

### 3.1. A Framework for the Investigated Material Handling Process

The problem studied in this paper is a case of a production line supply process modeled as a so-called “single-operator-multi-material” system. The model covers transportation of “M” number of raw materials from a supermarket to a dedicated production line on the production area and transportation of “E” number of empties from this production line to the supermarket through manual material handling with a simple cart by a logistics operator. The investigated production environment has been given as a concrete layout that determines the place of the supermarket, the place of the dedicated production line, and the path between these two components of the layout. The physical process of the line supply will contain the picking, storing, and pushing activities for manual handling of the aforementioned raw materials and empties. Figure 1 shows the production layout investigated, and Figure 2 illustrates the “single-operator-multi-material” job cycle as a decomposition of the mentioned activities.Figure 1
The production layout of investigated cycle.Figure 2
Decomposition of the complex material handling activity investigated and the heights and distances applied in the model.The job cycle consists of ten steps; the first step is picking the materials from the supermarkethpi,s, the second step is storing the materials in the cart from the supermarket hst,c, the third step is pushing the cart from the supermarket to the production area ds,p, the fourth step is picking the materials from the cart in the production area hpi,c, the fifth step is storing the materials to the production line in the production area hst,p, the sixth step is picking the empties from the production line in the production area hpi,p, the seventh step is storing the empties in the cart in the production area hst,c, the eighth step is pushing the cart from the production area to the supermarket ds,p, the ninth step is picking the empties from the cart hpi,c, and the tenth step is storing the empties to the supermarket hst,s and resting. The double-sided arrows represent the height of placement and distance parameters applied by us to the AIM tables to determine the time necessities of the operations in the job cycle (see Section 4.1).The notations of the realized heights and distances during the material handling process in the model, which describe the physical topology, are given in Table2.Table 2
The notations of heights and distances of the material handling process.
SymbolsDefinitionhst,pHeight of placement of the production line for storing (m)dc,pDistance between the cart and the production line (m)ds,pDistance between the supermarket and the production line (m)hpi,cHeight of placement of the cart for picking (m)hpi,sHeight of placement of the supermarket for picking (m)dc,sDistance between the cart and the supermarket (m)hst,cHeight of placement of the cart for storing (m)hpi,pHeight of placement of the production line for picking (m)hst,sHeight of placement of the supermarket for storing (m)


### 3.2. Objectives, Assumptions, and Notations

The objectives of this study are as follows:(1)
To develop a new lot-sizing model for the production line supply operation in a “single-operator-multi-material” case based on rest allowance and endurance time formulations developed by Rohmert [25, 26] to integrate the ergonomic aspects in the calculation of the optimal lot size.(2)
To investigate the outcomes with different scenarios of endurance time and rest allowance on the lot sizes and the total cost of production line supply operation.The following assumptions are made for developing the new model:(1)
The model considers one logistics operator to handle the materials in the line supplying process.(2)
The model considers the handling of “M” types of raw materials between the supermarket and the production line. For this reason, based on this assumption and the previous one, we name this model variation as a multi-material-single-operator model (referred to as the MMSO).(3)
In all material handling activities, empties of the same kind (i.e., boxes) have been used to handle the raw materials.(4)
A simple cart is used for transportation of the raw materials and empties between the loading points with a pushing motion.(5)
The aforementioned line supplying process is a cyclical activity during the operation of the production line. The count of the necessary line supplying transactions will strongly depend on the material needs and the availability time of the raw materials on the production line. Therefore, all of the related costs will also strongly depend on the operational properties of the line supply.(6)
The developed model is a discrete-event-based model, and it investigates the actual values of the parameters in all of the “t” time slots. It means that if the start of a line supplying transaction is necessary for an investigated “t” time slot considering the actual availability rules of the raw materials in the production line, then based on the decomposition of the process (see Figure 2), all of the time necessities and the related cost components will be calculated.The following notations are used throughout the paper (Table3).Table 3
Notations used in the mathematical model.
SymbolsDefinitionOTtOperational time needed to be generated at a time slot “t” (min)OTmat,tOperational time needed for handling the material “mat” at a time slot “t” (min)OTem,tTotal operational time needed for handling the empty “em” at a time slot “t” (min)OTmat,tpi,sOperational time needed for picking the material “mat” from the supermarket at a time slot “t” (min)OTmat,tst,cOperational time needed for storing the material “mat” to the cart at a time slot “t” (min)OTmat,tpuOperational time needed for pushing the material “mat” to the production line at a time slot “t” (min)OTmat,tpi,cOperational time needed for picking the material “mat” from the cart at a time slot “t” (min)OTmat,tst,pOperational time needed for storing the material “mat” to the production line at a time slot “t” (min)OTem,tpi,pOperational time needed for picking the empty “em” from the production line at a time slot “t” (min)OTem,tst,cOperational time needed for storing the empty “em” to the cart at a time slot “t” (min)OTem,tpuOperational time needed for pushing the empty “em” to the supermarket at a time slot “t” (min)OTem,tpi,cOperational time needed for picking the empty “em” from the cart at a time slot “t” (min)OTem,tst,sOperational time needed for storing the empty “em” to the supermarket at a time slot “t” (min)wcaWeight of the cart (kg)wmatWeight of the material “mat” (kg/pcs)wemWeight of the empty “em” (kg/pcs)wmatuniWeight of the handled unit in case of material “mat” (kg)ttactTakt time of the production (min/product)rmatThe material requirement in the production line per product of material “mat” (pcs/product)xmat,tBinary variable for the material “mat” at a time slot “t”qmatLot size of the material “mat” (pcs)xem,tBinary variable for the empty “em” at a time slot “t”tmatavAvailability of the material “mat” on the production line (min)RmatThe total material request of the production line in the examined time period (pcs)RTtRest time needed to be generated at a time slot “t” (min)RTmat,tRest time needed for handling the material “mat” at a time slot “t” (min)RTem,tRest time needed for handling the empty “em” at a time slot “t” (min)RTmat,tpi,sRest time needed for picking the material “mat” from the supermarket at a time slot “t” (min)RTmat,tst,cRest time needed for storing the material “mat” to the cart at a time slot “t” (min)RTmat,tpuRest time needed for pushing the material “mat” to the production line at a time slot “t” (min)RTmat,tpi,cRest time needed for picking the material “mat” from the cart at a time slot “t” (min)RTmat,tst,pRest time needed for storing the material “mat” to the production line at a time slot “t” (min)RTem,tpi,pRest time needed for picking the empty unit “em” from the production line at a time slot “t” (min)RTem,tst,cRest time needed for storing the empty unit “em” to the cart at a time slot “t” (min)RTem,tpuRest time needed for pushing the empty unit “em” to the supermarket at a time slot “t” (min)RTem,tpi,cRest time needed for picking the empty unit “em” from the cart at a time slot “t” (min)RTem,tst,sRest time needed for storing the empty unit “em” to the supermarket at a time slot “t” (min)fMVCmatuniRelative force during the handling of one piece of unit load in case of material “mat” (%)fMVCemRelative force during the handling of one piece of empty load in case of empty “em” (%)Fmaxpi&stMaximum force load of the human spine for picking and storing (N)FmaxpuMaximum force load of the human spine for pushing (N)cwUnit worker wage ($/h)ch,matInventory holding cost ($/pcs/h)CtRTCost of the rest time at a time slot “t” ($)CtOPCost of the operation at a time slot “t” ($)COPThe total cost of the operation ($)CHOThe total cost of the holding ($)CtHOCost of the holding at a time slot “t” ($)CtOTCost of the operation time at a time slot “t” ($)Qmat,tClosing stock of the material “mat” on the line at a time slot “t” (pcs)CTotal cost ($)


## 3.1. A Framework for the Investigated Material Handling Process

The problem studied in this paper is a case of a production line supply process modeled as a so-called “single-operator-multi-material” system. The model covers transportation of “M” number of raw materials from a supermarket to a dedicated production line on the production area and transportation of “E” number of empties from this production line to the supermarket through manual material handling with a simple cart by a logistics operator. The investigated production environment has been given as a concrete layout that determines the place of the supermarket, the place of the dedicated production line, and the path between these two components of the layout. The physical process of the line supply will contain the picking, storing, and pushing activities for manual handling of the aforementioned raw materials and empties. Figure 1 shows the production layout investigated, and Figure 2 illustrates the “single-operator-multi-material” job cycle as a decomposition of the mentioned activities.Figure 1
The production layout of investigated cycle.Figure 2
Decomposition of the complex material handling activity investigated and the heights and distances applied in the model.The job cycle consists of ten steps; the first step is picking the materials from the supermarkethpi,s, the second step is storing the materials in the cart from the supermarket hst,c, the third step is pushing the cart from the supermarket to the production area ds,p, the fourth step is picking the materials from the cart in the production area hpi,c, the fifth step is storing the materials to the production line in the production area hst,p, the sixth step is picking the empties from the production line in the production area hpi,p, the seventh step is storing the empties in the cart in the production area hst,c, the eighth step is pushing the cart from the production area to the supermarket ds,p, the ninth step is picking the empties from the cart hpi,c, and the tenth step is storing the empties to the supermarket hst,s and resting. The double-sided arrows represent the height of placement and distance parameters applied by us to the AIM tables to determine the time necessities of the operations in the job cycle (see Section 4.1).The notations of the realized heights and distances during the material handling process in the model, which describe the physical topology, are given in Table2.Table 2
The notations of heights and distances of the material handling process.
SymbolsDefinitionhst,pHeight of placement of the production line for storing (m)dc,pDistance between the cart and the production line (m)ds,pDistance between the supermarket and the production line (m)hpi,cHeight of placement of the cart for picking (m)hpi,sHeight of placement of the supermarket for picking (m)dc,sDistance between the cart and the supermarket (m)hst,cHeight of placement of the cart for storing (m)hpi,pHeight of placement of the production line for picking (m)hst,sHeight of placement of the supermarket for storing (m)


## 3.2. Objectives, Assumptions, and Notations

The objectives of this study are as follows:(1)
To develop a new lot-sizing model for the production line supply operation in a “single-operator-multi-material” case based on rest allowance and endurance time formulations developed by Rohmert [25, 26] to integrate the ergonomic aspects in the calculation of the optimal lot size.(2)
To investigate the outcomes with different scenarios of endurance time and rest allowance on the lot sizes and the total cost of production line supply operation.The following assumptions are made for developing the new model:(1)
The model considers one logistics operator to handle the materials in the line supplying process.(2)
The model considers the handling of “M” types of raw materials between the supermarket and the production line. For this reason, based on this assumption and the previous one, we name this model variation as a multi-material-single-operator model (referred to as the MMSO).(3)
In all material handling activities, empties of the same kind (i.e., boxes) have been used to handle the raw materials.(4)
A simple cart is used for transportation of the raw materials and empties between the loading points with a pushing motion.(5)
The aforementioned line supplying process is a cyclical activity during the operation of the production line. The count of the necessary line supplying transactions will strongly depend on the material needs and the availability time of the raw materials on the production line. Therefore, all of the related costs will also strongly depend on the operational properties of the line supply.(6)
The developed model is a discrete-event-based model, and it investigates the actual values of the parameters in all of the “t” time slots. It means that if the start of a line supplying transaction is necessary for an investigated “t” time slot considering the actual availability rules of the raw materials in the production line, then based on the decomposition of the process (see Figure 2), all of the time necessities and the related cost components will be calculated.The following notations are used throughout the paper (Table3).Table 3
Notations used in the mathematical model.
SymbolsDefinitionOTtOperational time needed to be generated at a time slot “t” (min)OTmat,tOperational time needed for handling the material “mat” at a time slot “t” (min)OTem,tTotal operational time needed for handling the empty “em” at a time slot “t” (min)OTmat,tpi,sOperational time needed for picking the material “mat” from the supermarket at a time slot “t” (min)OTmat,tst,cOperational time needed for storing the material “mat” to the cart at a time slot “t” (min)OTmat,tpuOperational time needed for pushing the material “mat” to the production line at a time slot “t” (min)OTmat,tpi,cOperational time needed for picking the material “mat” from the cart at a time slot “t” (min)OTmat,tst,pOperational time needed for storing the material “mat” to the production line at a time slot “t” (min)OTem,tpi,pOperational time needed for picking the empty “em” from the production line at a time slot “t” (min)OTem,tst,cOperational time needed for storing the empty “em” to the cart at a time slot “t” (min)OTem,tpuOperational time needed for pushing the empty “em” to the supermarket at a time slot “t” (min)OTem,tpi,cOperational time needed for picking the empty “em” from the cart at a time slot “t” (min)OTem,tst,sOperational time needed for storing the empty “em” to the supermarket at a time slot “t” (min)wcaWeight of the cart (kg)wmatWeight of the material “mat” (kg/pcs)wemWeight of the empty “em” (kg/pcs)wmatuniWeight of the handled unit in case of material “mat” (kg)ttactTakt time of the production (min/product)rmatThe material requirement in the production line per product of material “mat” (pcs/product)xmat,tBinary variable for the material “mat” at a time slot “t”qmatLot size of the material “mat” (pcs)xem,tBinary variable for the empty “em” at a time slot “t”tmatavAvailability of the material “mat” on the production line (min)RmatThe total material request of the production line in the examined time period (pcs)RTtRest time needed to be generated at a time slot “t” (min)RTmat,tRest time needed for handling the material “mat” at a time slot “t” (min)RTem,tRest time needed for handling the empty “em” at a time slot “t” (min)RTmat,tpi,sRest time needed for picking the material “mat” from the supermarket at a time slot “t” (min)RTmat,tst,cRest time needed for storing the material “mat” to the cart at a time slot “t” (min)RTmat,tpuRest time needed for pushing the material “mat” to the production line at a time slot “t” (min)RTmat,tpi,cRest time needed for picking the material “mat” from the cart at a time slot “t” (min)RTmat,tst,pRest time needed for storing the material “mat” to the production line at a time slot “t” (min)RTem,tpi,pRest time needed for picking the empty unit “em” from the production line at a time slot “t” (min)RTem,tst,cRest time needed for storing the empty unit “em” to the cart at a time slot “t” (min)RTem,tpuRest time needed for pushing the empty unit “em” to the supermarket at a time slot “t” (min)RTem,tpi,cRest time needed for picking the empty unit “em” from the cart at a time slot “t” (min)RTem,tst,sRest time needed for storing the empty unit “em” to the supermarket at a time slot “t” (min)fMVCmatuniRelative force during the handling of one piece of unit load in case of material “mat” (%)fMVCemRelative force during the handling of one piece of empty load in case of empty “em” (%)Fmaxpi&stMaximum force load of the human spine for picking and storing (N)FmaxpuMaximum force load of the human spine for pushing (N)cwUnit worker wage ($/h)ch,matInventory holding cost ($/pcs/h)CtRTCost of the rest time at a time slot “t” ($)CtOPCost of the operation at a time slot “t” ($)COPThe total cost of the operation ($)CHOThe total cost of the holding ($)CtHOCost of the holding at a time slot “t” ($)CtOTCost of the operation time at a time slot “t” ($)Qmat,tClosing stock of the material “mat” on the line at a time slot “t” (pcs)CTotal cost ($)


## 4. The Methods and the Mathematical Model

### 4.1. The AIM Methodology

The AIM (“Anyagmozgatási Időszükséglet Meghatározás”; in English, “Determination of Material Handling Time Requirements”) tables were developed by Hungarian researchers at the Budapest University of Technology and Economics in the middle of the 1960s. The AIM tables are used to determine the time necessity of material handling systems operating periodically. The objectives of the AIM are to define the time necessity of a weight forwarding process, to determine how many machines, tools, and workers the system needs, to plan the schedule for tasks, and to analyze the working efficiency of existing systems [63]. The method allows for the planning of the material handling systems by hand without any tools and of the material handling systems by hand with tools and forklifts. Therefore, we apply the AIM method to estimate the material handling time necessities. The time estimation procedure used in the AIM method is very similar to the well-known MTM (Methods-Time Measurement) method; however, the AIM method has been specialized for logistics activities. The AIM method calculation parameters for the production line supply are shown in Figure 3 where the “h” parameters are the average levels of the picking and storing logistics units from and to the logistics equipment, and the “d” parameters refer to the estimated average length of the realized motions.Figure 3
The parameters of the AIM method for the job cycle.Additionally, we need to give the weight of the handled logistics units to be determined by the weight of the materials handled and the empty units applied. Furthermore, the weight of the logistics units handled may also depend on the lot sizes used indirectly. Based on these parameters, we can estimate the time necessities of the operations (see Appendix A).


### 4.2. Rest Allowance and Endurance Time

An excessive load on the spinal muscles affects neuromuscular spinal stability as it may cause muscle fatigue, musculoskeletal disorders, and low back pain [52]. The neuromuscular spinal stability is affected by vibration and repetitive loading [65, 66]. To reduce ergonomic risks, we calculate the rest time and maximum endurance time of the operator in accordance with the equation developed by Rohmert [25, 26]. The rest time and maximum endurance time equation developed by Rohmert [25, 26] can be defined as follows:(1)RT=1800⋅tTE1.4⋅fMVC−0.150.5,where “RT” is the rest time need if the contraction duration (working period of picking, storing, or pushing) is equal to “t” expressed in minutes, and the “TE” endurance time can be defined as follows:(2)TE=−1.5+2.1fMVC−0.6fMVC2+0.1fMVC3,where fMVC=Fload/Fmax is the relative force, Fload is the actual force load, and Fmax is defined as the maximum force load of the human spine based on the applied limitations in the case of the investigated activity.As given in Figure4, the force load on the spine is calculated based on the determination of the balance of the body. F1 represents the force load on the spine for picking and storing motions, and F3 represents the force load on the spine for pushing motion.Figure 4
The representation of the load on the spine for picking, storing, and pushing.We use the maximum endurance time as an ergonomic measurement method that is a valuable method for the assessment of specific muscle areas, body parts, and postures. As we aim to decrease the work-related musculoskeletal disorders of the back, the maximum endurance time is suitable for the ergonomic measurement of our model. Other ergonomic measurement methods have some limitations, such as the energy expenditure known as a common method, and have considered the whole-body energy expenditure for the assessment. Furthermore, it may not be applicable to assess the workload on the spine muscles as well as the disorders of the back. The OWAS, which analyzes working postures based on 72 different postures, has a limitation as it does not consider any aspects of repetition or duration of the sequential postures. The NIOSH Lifting Equation, which helps to assess the ergonomic risks of lifting and lowering tasks, has a limitation on the other motions, such as pulling, pushing, and carrying. The force load of picking and storing can be calculated as follows:(3)Floadpi&st=wpi&st⋅g⋅ba,where “wpi&st” is equal to the weight of the handled material during the picking and storing activity; “b” is the distance between the body axis and the load axis; “g” is the gravitational acceleration; and “a” is the length of the load on the spine. The force load of pushing can be calculated as follows:(4)Floadpu=wpu⋅g⋅μ⋅dc,where “wpu” is equal to the weight of the handled material during the pushing activity; “g” is the gravitational acceleration; “μ” is the coefficient of rolling friction; “d” is the distance between the shoulder and the waist; and “c” is the height of the load on the spine. Additionally, the Fmax parameter for the investigated activities can be calculated as follows:(5)Fmaxpi&st=wpi&stmax⋅g⋅ba,for picking and storing, and(6)Fmaxpu=wpumax⋅g⋅μ⋅dc,for pushing, where the maximum force load for picking and storing Fmaxpi&st is equal to 2452.5 N, and it is calculated according to ISO-11228 on Ergonomics-Manual Handling, Part 1: Lifting and Carrying [67], where the limit for two-handed lifting is 25 kg. The maximum force load for pushing Fmaxpu is calculated as 4326.2 N as suggested by Resnik and Chaffin [68] for the load limit of 225 kg for four-wheeled carts. According to the relaxation formulation, rest time is necessary if the fMVC is higher than 0.15. With this constraint, “RT” can be expressed by the linear function of the tangent plane approximation (see Appendix B):(7)RT=Are⋅t+Bre⋅fMVC−Cre,where Are is the coefficient of contraction time, Bre is the coefficient of relative force, and Cre is the constant of rest time. The linear function of rest time and endurance time can be defined as follows:(8)RT=3.6419⋅t+32.2251⋅fMVC−5.3225.


### 4.3. The Mathematical Model

In our discrete-event-based model, we can simulate fort = 1 to T interval according to the actual values of the aforementioned parameters, which describe the effects of investigation in the production line supply. We introduce the following binary decision variables. xmat,t indicates whether the material “mat” is placed in the cart at the “t-th” time slot xmat,t=1 or not xmat,t=0. Likewise, xem,t indicates whether the empty unit “em” is placed in the cart at the “t-th” time slot xem,t=1 or not xem,t=0. Based on the problem description in the third section and the relevant cost components, the calculation of the different time components in our MMSO model is as follows.Operational time need (OT):(9)OTt=OTmat,t+OTem,t.The operational time need for the materials to be handled will be determined by the AIM method:(10)OTmat,t=∑mat=1MOTmat,tpi,s+OTmat,tst,c+OTmat,tpu+OTmat,tpi,c+OTmat,tst,p,where(11)OTmat,tpi,s=fAIMwmatuni,hpi,s,dc,s,OTmat,tst,c=fAIMwmatuni,hst,c,dc,s,OTmat,tpu=fAIM∑mat=1Mxmat,t⋅wmatuni,wca,ds,p,xmat,t=1,if the material″mat″is placed in the cart at the″t″−th time slot,0,otherwise,OTmat,tpi,c=fAIMwmatuni,hpi,c,dc,p,OTmat,tst,p=fAIMwmatuni,hst,p,dc,s,where(12)wmatuni=fwunitqmat,wmat,wem=qmat⋅wmat+wem,qmat=fqmattmatav,ttact,rmat=tmatavttact⋅rmat,where MINtmatav≥MAXOTt+RTt,where OTt is operational time need to be generated at a time slot “t.” OTmat,t is the operational time need for handling the material “mat” at a time slot “t.” OTmat,tpi,s is the operational time need for picking the material “mat” from the supermarket at a time slot “t,” and it is a function of the weight of the handled unit for material “mat” 0<wmatuni<25kg, the height of placement of the supermarket for picking 0≤hpi,s<2m, and the distance between the cart and the supermarket dc,s>0. OTmat,tst,c is the operational time need for storing the material “mat” to the cart at a time slot “t,” and it is a function of the weight of the handled unit for material “mat” wmatuni, the height of placement of the cart for storing 0≤hst,c<1,5m, and distance between the cart and the supermarket dc,s. OTmat,tpu is the operational time need for pushing the material “mat” to the production line at a time slot “t,” and it is a function of the total weight of the materials “mat” placed in the cart at the time slot “t” (calculated as ∑mat=1Mxmat,t⋅wmatuni), the weight of the cart 0<wca<50kg, and distance between the supermarket and the production line ds,p>0. OTmat,tpi,c is the operational time need for picking the material “mat” from the cart at a time slot “t,” and it is a function of the weight of the handled unit for material “mat” wmatuni, the height of placement of the cart for picking 0≤hpi,c<1,5m, and distance between the cart and the production line dc,p>0. OTmat,tst,p is the operational time need for storing the material “mat” to the production line at a time slot “t,” and it is a function of the weight of the handled unit for material “mat” wmatuni, the height of placement of the production line for storing 0≤hst,p<2m, and distance between the cart and the supermarket dc,s.fAIM means that the operational time needs are determined by the aforementioned AIM method, as a function of the weight, height, and distance parameters of the operation. fwunit means that the weight of the handled logistics unit is interpreted as a function of the lot size (qmat>0 and qmat∈ℤ+), and the weight of the material 0<wmat≤1kg, and the applied empty unit 0<wem≤1kg. fqmat means that the lot size of the material is interpreted as a function of the material requirement per product (rmat>0 and rmat∈ℤ+), the applied tact time of the production line ttact>0, and the previous availability time of the material on the production line tmatav>0. In the model developed, we aim to use this availability time as a decision variable for all of the materials investigated and, in this way, to determine the lot sizes.The operational time need for the empty units will be determined by the AIM method:(13)OTem,t=∑em=1EOTem,tpi,p+OTem,tst,c+OTem,tpu+OTem,tpi,c+OTem,tst,s,where(14)OTem,tpi,p=fAIMwem,hpi,p,dc,p,OTem,tst,c=fAIMwem,hst,c,dc,p,OTem,tpu=fAIM∑em=1Exem,t⋅wem,wca,ds,p,xem,t=1,if the empty“em”is placed in the cart at the“t”−th time slot,0,otherwise,OTem,tpi,c=fAIMwem,hpi,c,dc,s,OTem,tst,s=fAIMwem,hst,s,dc,s,where OTem,t is the total operational time need for handling the empty “em” at a time slot “t.” OTem,tpi,p is the operational time need for picking the empty “em” from the production line at a time slot “t,” and it is a function of the weight of the empty “em” wem, the height of placement of the production line for picking hpi,p, and distance between the cart and the production line dc,p. OTem,tst,c is the operational time need for storing the empty “em” to the cart at a time slot “t,” and it is a function of the weight of the empty “em” wem, the height of placement of the cart for storing hst,c, and distance between the cart and the production line dc,p. OTem,tpu is the operational time need for pushing the empty “em” to the supermarket at a time slot “t,” and it is a function of the total weight of the empties “em” placed in the cart at the time slot “t” (calculated as ∑em=1Exem,t⋅wem), the weight of the cart wca, and distance between the supermarket and the production line ds,p. OTem,tpi,c is the operational time need for picking the empty “em” from the cart at a time slot “t,” and it is a function of the weight of the empty “em” wem, the height of placement of the cart for picking hpi,c, and distance between the cart and the supermarket dc,s. OTem,tst,s is the operational time need for storing the empty “em” to the supermarket at a time slot “t,” and it is a function of the weight of the empty “em” wem, the height of placement of the supermarket for storing hst,s, and distance between the cart and the supermarket dc,s.Rest time need (RT):(15)RTt=RTmat,t+RTem,t.The rest time need for the materials to be handled:(16)RTmat,t=∑mat=1MRTmat,tpi,s+RTmat,tst,c+RTmat,tpu+RTmat,tpi,c+RTmat,tst,p,where(17)RTmat,tpi,s=Are⋅OTmat,tpi,s+Bre⋅fMVCmatuni−Cre,RTmat,tst,c=Are⋅OTmat,tst,c+Bre⋅fMVCmatuni−Cre,RTmat,tpu=Are⋅OTmat,tpu+Bre⋅fMVCmat,t−Cre,RTmat,tpi,c=Are⋅OTmat,tpi,c+Bre⋅fMVCmatuni−Cre,RTmat,tst,p=Are⋅OTmat,tst,p+Bre⋅fMVCmatuni−Cre,where(18)fMVCmatuni=wmatuni⋅g⋅b/aFmaxpi&st,fMVCmat,t=wca+∑mat=1Mxmat,t⋅wmatuni⋅g⋅μ⋅d/cFmaxpu,xmat,t=1,if the material“mat”is placed in the cart at the“t”−th time slot,0,otherwise,where RTt is the rest time need to be generated at a time slot “t.” RTmat,t is the rest time need for handling the material “mat” at a time slot “t.” RTmat,tpi,s is the rest time need for picking the material “mat” from the supermarket at a time slot “t,” and it is a function of the operational time need for picking the material “mat” from the supermarket at a time slot “t” OTmat,tpi,s, the relative force of the unit weight of the material “mat” for picking and storing motions fMVCmatuni, the coefficient of contraction time Are, the coefficient of relative force Bre, and the constant of rest time Cre. RTmat,tst,c is the rest time need for storing the material “mat” to the cart at a time slot “t,” and it is a function of the operational time need for storing the material “mat” to the cart at a time slot “t” OTmat,tst,c, the relative force of the unit weight of the material “mat” for picking and storing motions fMVCmatuni, the coefficient of contraction time Are, the coefficient of relative force Bre, and the constant of rest time Cre. RTmat,tpu is the rest time need for pushing the material “mat” to the production line at a time slot “t,” and it is a function of the operational time need for pushing the material “mat” to the production line at a time slot “t” OTmat,tpu, the relative force of the unit weight of the materials “mat” for pushing motion fMVCmat,t, the coefficient of contraction time Are, the coefficient of relative force Bre, and the constant of rest time Cre. RTmat,tpi,c is the rest time need for picking the material “mat” from the cart at a time slot “t,” and it is a function of the operational time need for picking the material “mat” from the cart at a time slot “t” OTmat,tpi,c, the relative force of the unit weight of the material “mat” for picking and storing motions fMVCmatuni, the coefficient of contraction time Are, the coefficient of relative force Bre, and the constant of rest time Cre. RTmat,tst,p is the rest time need for storing the material “mat” to the production line at a time slot “t,” and it is a function of the operational time need for storing the material “mat” to the production line at a time slot “t” OTmat,tst,p, the relative force of unit weight of the material “mat” for picking and storing motions fMVCmatuni, the coefficient of contraction time Are, the coefficient of relative force Bre, and the constant of rest time Cre.The relative force for picking and storing and pushing activities can be calculated using the previously described methodology in Section4.2, where in this case, wpi&st=wmatuni and wpu=wca+∑matMxmat,t⋅wmatuni. Based on this, it can be seen that the calculated values also depend on the weight of the units handled wmatuni and therefore the lot sizes calculated qmat. To calculate the rest time needed, a linearized rest time formula (equation (7) was used.The rest time need for the empty units to be handled:(19)RTem,t=∑em=1ERTem,tpi,p+RTem,tst,c+RTem,tpu+RTem,tpi,c+RTem,tst,s,where(20)RTem,tpi,p=Are⋅OTem,tpi,p+Bre⋅fMVCem−Cre,RTem,tst,c=Are⋅OTem,tst,c+Bre⋅fMVCem−Cre,RTem,tpu=Are⋅OTem,tpu+Bre⋅fMVCem,t−Cre,RTem,tpi,c=Are⋅OTem,tpi,c+Bre⋅fMVCem−Cre,RTem,tst,s=Are⋅OTem,tst,s+Bre⋅fMVCem−Cre,where(21)fMVCem=wem⋅g⋅b/aFmaxpi&st,fMVCem,t=wca+∑em=1Exem,t⋅wem⋅g⋅μ⋅d/cFmaxpu,xem,t=1,if the empty“em”is placed in the cart at the“t”−th time slot,0,otherwise,where RTem,t is the rest time need for handling the empty “em” at a time slot “t.” RTem,tpi,p is the rest time need for picking the empty unit “em” from the production line at a time slot “t,” and it is a function of the operational time need for picking the empty “em” from the production line at a time slot “t” OTem,tpi,p, the relative force for the empty “em” for picking and storing motions fMVCem, the coefficient of contraction time Are, the coefficient of relative force Bre, and the constant of rest time Cre. RTem,tst,c is the rest time need for storing the empty unit “em” to the cart at a time slot “t,” and it is a function of the operational time need for storing the empty “em” to the cart at a time slot “t” OTem,tst,c, the relative force for the empty “em” for picking and storing motions fMVCem, the coefficient of contraction time Are, the coefficient of relative force Bre, and the constant of rest time Cre. RTem,tpu is the rest time need for pushing the empty unit “em” to the supermarket at a time slot “t,” and it is a function of the operational time need for pushing the empty “em” to the supermarket at a time slot “t” OTem,tpu, the relative force for empties “em” for pushing motion fMVCem,t, the coefficient of contraction time Are, the coefficient of relative force Bre, and the constant of rest time Cre. RTem,tpi,c is the rest time need for picking the empty unit “em” from the cart at a time slot “t,” and it is a function of the operational time need for picking the empty “em” from the cart at a time slot “t” OTem,tpi,c, the relative force for the empty “em” for picking and storing motions fMVCem, the coefficient of contraction time Are, the coefficient of relative force Bre, and the constant of rest time Cre. RTem,tst,s is the rest time need for storing the empty unit “em” to the supermarket at a time slot “t,” and it is a function of the operational time need for storing the empty “em” to the supermarket at a time slot “t” OTem,tst,s, the relative force for the empty “em” for picking and storing motions fMVCem, the coefficient of contraction time Are, the coefficient of relative force Bre, and the constant of rest time Cre.The calculation of relative force and the rest time need can be realized in a similar way as in the previously mentioned case, but in this case,wpi&st=wem and wpu=wca+∑em=1Exem,t⋅wem equations must be considered.The cost of operation, the cost of rest time, and the cost of inventory holding can be described as follows:(22)COP=∑t=1TCtOP=∑t=1Tcw⋅OTt,CRT=∑t=1TCtRT=∑t=1Tcw⋅RTt,CHO=∑t=1TCtHO=∑t=1T∑mat=1Mch,mat⋅Qmat,t.The total cost can be calculated as follows:(23)C=COP+CHO+CRT.The objective function can be written as follows:(24)min∑t=1Tcw⋅OTt+RTt+∑mat=1Mch,mat⋅Qmat,t.


## 4.1. The AIM Methodology

The AIM (“Anyagmozgatási Időszükséglet Meghatározás”; in English, “Determination of Material Handling Time Requirements”) tables were developed by Hungarian researchers at the Budapest University of Technology and Economics in the middle of the 1960s. The AIM tables are used to determine the time necessity of material handling systems operating periodically. The objectives of the AIM are to define the time necessity of a weight forwarding process, to determine how many machines, tools, and workers the system needs, to plan the schedule for tasks, and to analyze the working efficiency of existing systems [63]. The method allows for the planning of the material handling systems by hand without any tools and of the material handling systems by hand with tools and forklifts. Therefore, we apply the AIM method to estimate the material handling time necessities. The time estimation procedure used in the AIM method is very similar to the well-known MTM (Methods-Time Measurement) method; however, the AIM method has been specialized for logistics activities. The AIM method calculation parameters for the production line supply are shown in Figure 3 where the “h” parameters are the average levels of the picking and storing logistics units from and to the logistics equipment, and the “d” parameters refer to the estimated average length of the realized motions.Figure 3
The parameters of the AIM method for the job cycle.Additionally, we need to give the weight of the handled logistics units to be determined by the weight of the materials handled and the empty units applied. Furthermore, the weight of the logistics units handled may also depend on the lot sizes used indirectly. Based on these parameters, we can estimate the time necessities of the operations (see Appendix A).


## 4.2. Rest Allowance and Endurance Time

An excessive load on the spinal muscles affects neuromuscular spinal stability as it may cause muscle fatigue, musculoskeletal disorders, and low back pain [52]. The neuromuscular spinal stability is affected by vibration and repetitive loading [65, 66]. To reduce ergonomic risks, we calculate the rest time and maximum endurance time of the operator in accordance with the equation developed by Rohmert [25, 26]. The rest time and maximum endurance time equation developed by Rohmert [25, 26] can be defined as follows:(1)RT=1800⋅tTE1.4⋅fMVC−0.150.5,where “RT” is the rest time need if the contraction duration (working period of picking, storing, or pushing) is equal to “t” expressed in minutes, and the “TE” endurance time can be defined as follows:(2)TE=−1.5+2.1fMVC−0.6fMVC2+0.1fMVC3,where fMVC=Fload/Fmax is the relative force, Fload is the actual force load, and Fmax is defined as the maximum force load of the human spine based on the applied limitations in the case of the investigated activity.As given in Figure4, the force load on the spine is calculated based on the determination of the balance of the body. F1 represents the force load on the spine for picking and storing motions, and F3 represents the force load on the spine for pushing motion.Figure 4
The representation of the load on the spine for picking, storing, and pushing.We use the maximum endurance time as an ergonomic measurement method that is a valuable method for the assessment of specific muscle areas, body parts, and postures. As we aim to decrease the work-related musculoskeletal disorders of the back, the maximum endurance time is suitable for the ergonomic measurement of our model. Other ergonomic measurement methods have some limitations, such as the energy expenditure known as a common method, and have considered the whole-body energy expenditure for the assessment. Furthermore, it may not be applicable to assess the workload on the spine muscles as well as the disorders of the back. The OWAS, which analyzes working postures based on 72 different postures, has a limitation as it does not consider any aspects of repetition or duration of the sequential postures. The NIOSH Lifting Equation, which helps to assess the ergonomic risks of lifting and lowering tasks, has a limitation on the other motions, such as pulling, pushing, and carrying. The force load of picking and storing can be calculated as follows:(3)Floadpi&st=wpi&st⋅g⋅ba,where “wpi&st” is equal to the weight of the handled material during the picking and storing activity; “b” is the distance between the body axis and the load axis; “g” is the gravitational acceleration; and “a” is the length of the load on the spine. The force load of pushing can be calculated as follows:(4)Floadpu=wpu⋅g⋅μ⋅dc,where “wpu” is equal to the weight of the handled material during the pushing activity; “g” is the gravitational acceleration; “μ” is the coefficient of rolling friction; “d” is the distance between the shoulder and the waist; and “c” is the height of the load on the spine. Additionally, the Fmax parameter for the investigated activities can be calculated as follows:(5)Fmaxpi&st=wpi&stmax⋅g⋅ba,for picking and storing, and(6)Fmaxpu=wpumax⋅g⋅μ⋅dc,for pushing, where the maximum force load for picking and storing Fmaxpi&st is equal to 2452.5 N, and it is calculated according to ISO-11228 on Ergonomics-Manual Handling, Part 1: Lifting and Carrying [67], where the limit for two-handed lifting is 25 kg. The maximum force load for pushing Fmaxpu is calculated as 4326.2 N as suggested by Resnik and Chaffin [68] for the load limit of 225 kg for four-wheeled carts. According to the relaxation formulation, rest time is necessary if the fMVC is higher than 0.15. With this constraint, “RT” can be expressed by the linear function of the tangent plane approximation (see Appendix B):(7)RT=Are⋅t+Bre⋅fMVC−Cre,where Are is the coefficient of contraction time, Bre is the coefficient of relative force, and Cre is the constant of rest time. The linear function of rest time and endurance time can be defined as follows:(8)RT=3.6419⋅t+32.2251⋅fMVC−5.3225.


## 4.3. The Mathematical Model

In our discrete-event-based model, we can simulate fort = 1 to T interval according to the actual values of the aforementioned parameters, which describe the effects of investigation in the production line supply. We introduce the following binary decision variables. xmat,t indicates whether the material “mat” is placed in the cart at the “t-th” time slot xmat,t=1 or not xmat,t=0. Likewise, xem,t indicates whether the empty unit “em” is placed in the cart at the “t-th” time slot xem,t=1 or not xem,t=0. Based on the problem description in the third section and the relevant cost components, the calculation of the different time components in our MMSO model is as follows.Operational time need (OT):(9)OTt=OTmat,t+OTem,t.The operational time need for the materials to be handled will be determined by the AIM method:(10)OTmat,t=∑mat=1MOTmat,tpi,s+OTmat,tst,c+OTmat,tpu+OTmat,tpi,c+OTmat,tst,p,where(11)OTmat,tpi,s=fAIMwmatuni,hpi,s,dc,s,OTmat,tst,c=fAIMwmatuni,hst,c,dc,s,OTmat,tpu=fAIM∑mat=1Mxmat,t⋅wmatuni,wca,ds,p,xmat,t=1,if the material″mat″is placed in the cart at the″t″−th time slot,0,otherwise,OTmat,tpi,c=fAIMwmatuni,hpi,c,dc,p,OTmat,tst,p=fAIMwmatuni,hst,p,dc,s,where(12)wmatuni=fwunitqmat,wmat,wem=qmat⋅wmat+wem,qmat=fqmattmatav,ttact,rmat=tmatavttact⋅rmat,where MINtmatav≥MAXOTt+RTt,where OTt is operational time need to be generated at a time slot “t.” OTmat,t is the operational time need for handling the material “mat” at a time slot “t.” OTmat,tpi,s is the operational time need for picking the material “mat” from the supermarket at a time slot “t,” and it is a function of the weight of the handled unit for material “mat” 0<wmatuni<25kg, the height of placement of the supermarket for picking 0≤hpi,s<2m, and the distance between the cart and the supermarket dc,s>0. OTmat,tst,c is the operational time need for storing the material “mat” to the cart at a time slot “t,” and it is a function of the weight of the handled unit for material “mat” wmatuni, the height of placement of the cart for storing 0≤hst,c<1,5m, and distance between the cart and the supermarket dc,s. OTmat,tpu is the operational time need for pushing the material “mat” to the production line at a time slot “t,” and it is a function of the total weight of the materials “mat” placed in the cart at the time slot “t” (calculated as ∑mat=1Mxmat,t⋅wmatuni), the weight of the cart 0<wca<50kg, and distance between the supermarket and the production line ds,p>0. OTmat,tpi,c is the operational time need for picking the material “mat” from the cart at a time slot “t,” and it is a function of the weight of the handled unit for material “mat” wmatuni, the height of placement of the cart for picking 0≤hpi,c<1,5m, and distance between the cart and the production line dc,p>0. OTmat,tst,p is the operational time need for storing the material “mat” to the production line at a time slot “t,” and it is a function of the weight of the handled unit for material “mat” wmatuni, the height of placement of the production line for storing 0≤hst,p<2m, and distance between the cart and the supermarket dc,s.fAIM means that the operational time needs are determined by the aforementioned AIM method, as a function of the weight, height, and distance parameters of the operation. fwunit means that the weight of the handled logistics unit is interpreted as a function of the lot size (qmat>0 and qmat∈ℤ+), and the weight of the material 0<wmat≤1kg, and the applied empty unit 0<wem≤1kg. fqmat means that the lot size of the material is interpreted as a function of the material requirement per product (rmat>0 and rmat∈ℤ+), the applied tact time of the production line ttact>0, and the previous availability time of the material on the production line tmatav>0. In the model developed, we aim to use this availability time as a decision variable for all of the materials investigated and, in this way, to determine the lot sizes.The operational time need for the empty units will be determined by the AIM method:(13)OTem,t=∑em=1EOTem,tpi,p+OTem,tst,c+OTem,tpu+OTem,tpi,c+OTem,tst,s,where(14)OTem,tpi,p=fAIMwem,hpi,p,dc,p,OTem,tst,c=fAIMwem,hst,c,dc,p,OTem,tpu=fAIM∑em=1Exem,t⋅wem,wca,ds,p,xem,t=1,if the empty“em”is placed in the cart at the“t”−th time slot,0,otherwise,OTem,tpi,c=fAIMwem,hpi,c,dc,s,OTem,tst,s=fAIMwem,hst,s,dc,s,where OTem,t is the total operational time need for handling the empty “em” at a time slot “t.” OTem,tpi,p is the operational time need for picking the empty “em” from the production line at a time slot “t,” and it is a function of the weight of the empty “em” wem, the height of placement of the production line for picking hpi,p, and distance between the cart and the production line dc,p. OTem,tst,c is the operational time need for storing the empty “em” to the cart at a time slot “t,” and it is a function of the weight of the empty “em” wem, the height of placement of the cart for storing hst,c, and distance between the cart and the production line dc,p. OTem,tpu is the operational time need for pushing the empty “em” to the supermarket at a time slot “t,” and it is a function of the total weight of the empties “em” placed in the cart at the time slot “t” (calculated as ∑em=1Exem,t⋅wem), the weight of the cart wca, and distance between the supermarket and the production line ds,p. OTem,tpi,c is the operational time need for picking the empty “em” from the cart at a time slot “t,” and it is a function of the weight of the empty “em” wem, the height of placement of the cart for picking hpi,c, and distance between the cart and the supermarket dc,s. OTem,tst,s is the operational time need for storing the empty “em” to the supermarket at a time slot “t,” and it is a function of the weight of the empty “em” wem, the height of placement of the supermarket for storing hst,s, and distance between the cart and the supermarket dc,s.Rest time need (RT):(15)RTt=RTmat,t+RTem,t.The rest time need for the materials to be handled:(16)RTmat,t=∑mat=1MRTmat,tpi,s+RTmat,tst,c+RTmat,tpu+RTmat,tpi,c+RTmat,tst,p,where(17)RTmat,tpi,s=Are⋅OTmat,tpi,s+Bre⋅fMVCmatuni−Cre,RTmat,tst,c=Are⋅OTmat,tst,c+Bre⋅fMVCmatuni−Cre,RTmat,tpu=Are⋅OTmat,tpu+Bre⋅fMVCmat,t−Cre,RTmat,tpi,c=Are⋅OTmat,tpi,c+Bre⋅fMVCmatuni−Cre,RTmat,tst,p=Are⋅OTmat,tst,p+Bre⋅fMVCmatuni−Cre,where(18)fMVCmatuni=wmatuni⋅g⋅b/aFmaxpi&st,fMVCmat,t=wca+∑mat=1Mxmat,t⋅wmatuni⋅g⋅μ⋅d/cFmaxpu,xmat,t=1,if the material“mat”is placed in the cart at the“t”−th time slot,0,otherwise,where RTt is the rest time need to be generated at a time slot “t.” RTmat,t is the rest time need for handling the material “mat” at a time slot “t.” RTmat,tpi,s is the rest time need for picking the material “mat” from the supermarket at a time slot “t,” and it is a function of the operational time need for picking the material “mat” from the supermarket at a time slot “t” OTmat,tpi,s, the relative force of the unit weight of the material “mat” for picking and storing motions fMVCmatuni, the coefficient of contraction time Are, the coefficient of relative force Bre, and the constant of rest time Cre. RTmat,tst,c is the rest time need for storing the material “mat” to the cart at a time slot “t,” and it is a function of the operational time need for storing the material “mat” to the cart at a time slot “t” OTmat,tst,c, the relative force of the unit weight of the material “mat” for picking and storing motions fMVCmatuni, the coefficient of contraction time Are, the coefficient of relative force Bre, and the constant of rest time Cre. RTmat,tpu is the rest time need for pushing the material “mat” to the production line at a time slot “t,” and it is a function of the operational time need for pushing the material “mat” to the production line at a time slot “t” OTmat,tpu, the relative force of the unit weight of the materials “mat” for pushing motion fMVCmat,t, the coefficient of contraction time Are, the coefficient of relative force Bre, and the constant of rest time Cre. RTmat,tpi,c is the rest time need for picking the material “mat” from the cart at a time slot “t,” and it is a function of the operational time need for picking the material “mat” from the cart at a time slot “t” OTmat,tpi,c, the relative force of the unit weight of the material “mat” for picking and storing motions fMVCmatuni, the coefficient of contraction time Are, the coefficient of relative force Bre, and the constant of rest time Cre. RTmat,tst,p is the rest time need for storing the material “mat” to the production line at a time slot “t,” and it is a function of the operational time need for storing the material “mat” to the production line at a time slot “t” OTmat,tst,p, the relative force of unit weight of the material “mat” for picking and storing motions fMVCmatuni, the coefficient of contraction time Are, the coefficient of relative force Bre, and the constant of rest time Cre.The relative force for picking and storing and pushing activities can be calculated using the previously described methodology in Section4.2, where in this case, wpi&st=wmatuni and wpu=wca+∑matMxmat,t⋅wmatuni. Based on this, it can be seen that the calculated values also depend on the weight of the units handled wmatuni and therefore the lot sizes calculated qmat. To calculate the rest time needed, a linearized rest time formula (equation (7) was used.The rest time need for the empty units to be handled:(19)RTem,t=∑em=1ERTem,tpi,p+RTem,tst,c+RTem,tpu+RTem,tpi,c+RTem,tst,s,where(20)RTem,tpi,p=Are⋅OTem,tpi,p+Bre⋅fMVCem−Cre,RTem,tst,c=Are⋅OTem,tst,c+Bre⋅fMVCem−Cre,RTem,tpu=Are⋅OTem,tpu+Bre⋅fMVCem,t−Cre,RTem,tpi,c=Are⋅OTem,tpi,c+Bre⋅fMVCem−Cre,RTem,tst,s=Are⋅OTem,tst,s+Bre⋅fMVCem−Cre,where(21)fMVCem=wem⋅g⋅b/aFmaxpi&st,fMVCem,t=wca+∑em=1Exem,t⋅wem⋅g⋅μ⋅d/cFmaxpu,xem,t=1,if the empty“em”is placed in the cart at the“t”−th time slot,0,otherwise,where RTem,t is the rest time need for handling the empty “em” at a time slot “t.” RTem,tpi,p is the rest time need for picking the empty unit “em” from the production line at a time slot “t,” and it is a function of the operational time need for picking the empty “em” from the production line at a time slot “t” OTem,tpi,p, the relative force for the empty “em” for picking and storing motions fMVCem, the coefficient of contraction time Are, the coefficient of relative force Bre, and the constant of rest time Cre. RTem,tst,c is the rest time need for storing the empty unit “em” to the cart at a time slot “t,” and it is a function of the operational time need for storing the empty “em” to the cart at a time slot “t” OTem,tst,c, the relative force for the empty “em” for picking and storing motions fMVCem, the coefficient of contraction time Are, the coefficient of relative force Bre, and the constant of rest time Cre. RTem,tpu is the rest time need for pushing the empty unit “em” to the supermarket at a time slot “t,” and it is a function of the operational time need for pushing the empty “em” to the supermarket at a time slot “t” OTem,tpu, the relative force for empties “em” for pushing motion fMVCem,t, the coefficient of contraction time Are, the coefficient of relative force Bre, and the constant of rest time Cre. RTem,tpi,c is the rest time need for picking the empty unit “em” from the cart at a time slot “t,” and it is a function of the operational time need for picking the empty “em” from the cart at a time slot “t” OTem,tpi,c, the relative force for the empty “em” for picking and storing motions fMVCem, the coefficient of contraction time Are, the coefficient of relative force Bre, and the constant of rest time Cre. RTem,tst,s is the rest time need for storing the empty unit “em” to the supermarket at a time slot “t,” and it is a function of the operational time need for storing the empty “em” to the supermarket at a time slot “t” OTem,tst,s, the relative force for the empty “em” for picking and storing motions fMVCem, the coefficient of contraction time Are, the coefficient of relative force Bre, and the constant of rest time Cre.The calculation of relative force and the rest time need can be realized in a similar way as in the previously mentioned case, but in this case,wpi&st=wem and wpu=wca+∑em=1Exem,t⋅wem equations must be considered.The cost of operation, the cost of rest time, and the cost of inventory holding can be described as follows:(22)COP=∑t=1TCtOP=∑t=1Tcw⋅OTt,CRT=∑t=1TCtRT=∑t=1Tcw⋅RTt,CHO=∑t=1TCtHO=∑t=1T∑mat=1Mch,mat⋅Qmat,t.The total cost can be calculated as follows:(23)C=COP+CHO+CRT.The objective function can be written as follows:(24)min∑t=1Tcw⋅OTt+RTt+∑mat=1Mch,mat⋅Qmat,t.


## 5. Results

### 5.1. Application of the New Model

We applied a numerical example to investigate the cost effects of rest time, where the total request of the production during the examined time period (480 min) related to the first materialRmat1 is equal to 1440 pcs, that related to the second material Rmat2 is equal to 960 pcs, that related to the third material Rmat3 is equal to 480 pcs, the maximum force load of the human spine Fmax for picking and storing is equal to 2452.5 N, and it has been calculated according to IS0-11228 on Ergonomics-Manual Handling, Part 1: Lifting and Carrying [67], where the limit for two-handed lifting is 25 kg. Fmax for pushing is 4326.2 N as suggested by Resnik and Chaffin [68] for the load limit of 225 kg for four-wheeled carts. The cart weight wca is 15 kg, the time for picking from the supermarket to the warehouse area is 0.11 min for three materials, the time for storing to the cart is 0.05 min for three materials, the time for pushing from the supermarket to the production area is 1.4 min, the time for picking the materials from the cart is 0.07 min, the time for storing the materials to the storage for production line supply is 0.05 min, the time for picking the empties from the storage near the production line is 0.07 min, the time for storing the empties to the cart is 0.05 min, the time for pushing the empties from the production area to the supermarket at the warehouse is 0.7 min, the time for picking the empties from the cart is 0.07 min, and the time for storing the empties to the supermarket is 0.09 min. Furthermore, the examined time period (T) is 480 min, the material need for the first material in the production line rmat1 is 3 pcs/product, the material need for the second material in the production line rmat2 is 2 pcs/product, the material need for the third material in the production line rmat3 is 1 pcs/product, the tact time ttact for the production is 1 min, the unit worker wage cw is 18 $/h, the cost of inventory holding for the first material ch,mat1 is 0.02 $/pcs/h, the cost of inventory holding for the second material ch,mat2 is 0.04 $/pcs/h, the cost of inventory holding for the third material ch,mat3 is 0.06 $/pcs/h, the weight of the first material wmat1 is 0.1 kg, the weight of the second material wmat2 is 0.2 kg, the weight of the third material wmat3 is 0.3 kg, the weight of the empties wem is 1 kg, “μ” is 0.28, and “g” is 9.81 m/s2. The height of placement of the production line for storing hst,p is 1.2 m, the height of placement of the cart for picking hpi,c is 0.8 m, the height of placement of the supermarket for picking hpi,s is 1.8 m, the height of placement of the cart for storing hst,c is 0.8 m, the height of placement of the production line for picking hpi,p is 1.2 m, the height of placement of the supermarket for storing hst,s is 1.8 m, the distance between the cart and the production line dc,p is 0.5 m, the distance between the supermarket and the production line in the production area ds,p is 100 m, and the distance between the cart and the supermarket dc,s is 0.5 m. Anthropometrical parameters have been calculated as follows: “a” and “c” are equal to 0.05 m for picking, storing, and pushing, “b” is equal to 0.5 m for picking and storing, and “d” is equal to 0.35 m for pushing according to the work-related musculoskeletal disorders mentioned in the journal published by the Ministry of Labour and Social Security of Turkey [69].We have used the following formal algorithm to calculate the necessary parameters:Fort = 1 to TIfxmat,t=1 thenFor mat = 1 toMCalculation ofOTmat,tIffMVCmatuni≥0,15 then calculation of RTmat,tpi,s,RTmat,tst,c,RTmat,tpi,c,RTmat,tst,pElseiffMVCmatuni<0,15 then RTmat,tpi,s=RTmat,tst,c=RTmat,tpi,c=RTmat,tst,p=0IffMVCmat,t≥0,15 then calculation of RTmat,tpuElseiffMVCmat,t<0,15 then RTmat,tpu=0Calculation ofRTmat,tIfxem,t=1 thenFor em = 1 toECalculation ofOTem,tIffMVCem≥0,15 then calculation of RTem,tpi,p,RTem,tst,c,RTem,tpi,c,RTem,tst,sElseiffMVCem≥0,15 then RTem,tpi,p=RTem,tst,c=RTem,tpi,c=RTem,tst,s=0IffMVCem,t≥0,15 then calculation of RTem,tpuElseiffMVCem,t≥0,15 then RTem,tpu=0Calculation ofRTem,tCalculation ofOTt and RTtIfMINtmatav≥MAXOTt+RTt then calculation of CtOP, CtRT, and CtHOElseifMINtmatav<MAXOTt+RTt then the situation is not valid, and no solution exists


#### 5.1.1. Calculation ofC for the Examination Period T

The applied formal algorithm above can be described withTn=4n2+4n2+n+n+1 function. Therefore, the time complexity of the algorithm based on the Big O principle is On2.We applied our simulation model (see Figure5) with the aforementioned algorithm to calculate the values of the parameters investigated and to find the ergonomically optimal lot sizes by using the solver application of Microsoft Excel for three materials where the total cost is minimum, MINtmatav≥MAXOTt+RTtandfMVC≥0,15. The minimum total cost equals to C=130,36$ and rest time cost equals toCRT=0where the ergonomically optimal lot size for the material one equals toqmat1 27 pcs, the material two equals toqmat2 14 pcs, and the material three equals toqmat3 9 pcs. The material availability of material one tmat1av is equal to 9 min, the material availability of material two tmat2av is equal to 7 min, the material availability of material three tmat3av is equal to 9 min, MINtmatav≥MAXOTt+RTt as MINtmatav=7min, and MAXOTt+RTt=4.39min where RTt=0. RTt=0 means that a logistics operator can work continuously without any need for rest time with these ergonomically optimal lot sizes. We have calculated the total cost of the maximum lifting and carrying limit (25 kg) according to ISO-11228 on Ergonomics-Manual Handling, Part 1: Lifting and Carrying [67], for two-handed lifting where each lot size is equal to 25 kg for the production supply process, and the lot size of material one qmat1 is equal to 240 pcs, the lot size of material two qmat2 is equal to 120 pcs, and the lot size of material three qmat3 is equal to 80 pcs. The total cost is equal to C=149,93$ and CRT=78,95$. CRT=78,95$ is the cost of the need for rest period of logistics operators to prevent any ergonomic risks. In our model, the saving achieved by using the ergonomic rest time and maximum endurance time in comparison to the total cost of the production supply process without ergonomics is equal to 13.05%.Figure 5
The result of a simulation experiment in case of the investigated numerical example in our specific simulator developed in MS Excel.


### 5.2. Parameteric Analysis

In this section, we have applied the mathematical model developed to different scenarios based on different parameter values derived from various industrial practices, and we have analyzed their impact on our new model. In addition, we have changed the weight of the three items from 0.01 to 3 kg in order to observe the effect of different materials. We have calculated the savings obtained by using our new model with differentwmatuni and different material needs in the production line for three different products. The savings have been calculated in comparison with the total cost of maximum lifting and carrying limit (25 kg) according to ISO-11228 on Ergonomics-Manual Handling, Part 1: Lifting and Carrying, for two-handed lifting where each lot size is equal to 25 kg and the optimal total cost of our model. Furthermore, we have calculated the savings where(25)rmat1=3pcs;rmat2=2pcs;rmat3=1pcs;rmat1=5pcs;rmat2=3pcs;rmat3=1pcs;rmat1=rmat2=rmat3=1pcs.Figure6 illustrates that savings increases or decreases by increase of the unit weight of items.. As we have calculated the savings regarding the maximum lifting and carrying limit determined by the ISO-11228 [67], the rest allowance stays constant as it increases each unit weight of materials “mat” but decreases the number of materials handled by causing an increase in the rest time cost and the operational cost.Figure 6
The savings for different weights of each material.We also conducted regression analysis for a better understanding of the relationship between variables. The regression analysis showed the relationship between the weight of handled materials and savings.rmat1=3pcs;rmat2=2pcs;rmat3=1pcs are statistically significant based on p value (0.007); rmat1=5pcs;rmat2=3pcs;rmat3=1pcs are statistically significant based on p value (0.025); and rmat1=rmat2=rmat3=1pcs is statistically significant based on p value (0.001).We have analyzed the material need of the first material in the production linermat1 as 1 pcs/product, the material need of the second material in the production line rmat2 as 1 pcs/product, and the material need of the third material in the production line rmat3 as 1 pcs/product, and, based on these different lot size values, we also analyzed how the total, operational, holding, and rest time costs change. Furthermore, we have applied the inventory holding cost of the first material ch,mat1 as 0.04 $/pcs/h, the inventory holding cost of the second material ch,mat2 as 0.08 $/pcs/h, the inventory holding cost of the third material ch,mat3 as 0.1 $/pcs/h, the weight of the first material wmat1 as 0.2 kg, the weight of the second material wmat2 as 0.3 kg, the weight of the third material wmat3 as 0.4 kg, the weight of the empties wem as 1 kg, “μ” as 0.28, and “g” as 9.81 m/s2, and other parameters were constant. The analysis results are given in Figure 7, where the lot size increases, especially the rest time, the cost of rest allowance, and the holding cost. As can be seen, the optimal solution of the model suggests that the lot size needs to be decreased for a decrease in the total cost, the holding cost, and the rest time cost.Figure 7
The total cost curve and additional cost curves for the holding cost, the operational cost, and the rest time cost of material one based on different lot sizes where the weights of the items are 0.2 kg, 0.3 kg, and 0.4 kg.To analyze the effect of the height of placement of the materials, we have applied different values:hst,p (1.7, 1.6, 1.5, 1.4, 1.3, 1.2, 1.1, 1, 0.9, 0.8, 0.7, 0.6, 0.5).hpi,c (1.4, 1.3, 1.2, 1.1, 1, 0.9, 0.8, 0.7, 0.6, 0.5, 0.4, 0.3, 0.2).hpi,s (2, 1.9, 1.8, 1.7, 1.6, 1.5, 1.4, 1.3, 1.2, 1.1, 1, 0.9, 0.8).hst,c (1.4, 1.3, 1.2, 1.1, 1, 0.9, 0.8, 0.7, 0.6, 0.5, 0.4, 0.3, 0.2).hpi,p (1.7, 1.6, 1.5, 1.4, 1.3, 1.2, 1.1, 1, 0.9, 0.8, 0.7, 0.6, 0.5).hst,s (2, 1.9, 1.8, 1.7, 1.6, 1.5, 1.4, 1.3, 1.2, 1.1, 1, 0.9, 0.8).As can be seen in Figure8, a higher height for placement and picking requires more rest time to prevent ergonomic risks and rest time cost because it may cause a more serious contraction of the muscles when the operator reaches a higher height. Therefore, it increases the need for rest time and the rest time cost as the operator needs to have rest to recover from a longer contraction. We also conducted regression analysis for a better understanding of the relationship between the height of placement for picking and storing and rest time. The regression analysis shows that the relationship between the height of placement for picking of materials hpi,c, hpi,s, hpi,p and rest time is statistically significant based on p value (0.001) and the relationship between the height of placement for storing of materials hst,c, hst,s, hst,p and rest time is statistically significant based on a p value (0.001).Figure 8
The rest time cost for different heights of placement and picking.Furthermore, we have analyzed the effect of the unit item weight on the lot size of materials. As shown in Figure9, the increase in the unit weight of the items reduces the optimal lot size values. The optimal lot size values represent the lot sizes where a logistics operator can work continuously and the need for rest time is equal to zero. This means that a lot size higher than the optimal lot size increases the need for rest time and that the logistics operator must have necessary rest time after handling the lot to prevent ergonomic risks due to higher load values. Higher unit item weight increases the contraction of the muscles that leads to increased fatigue and injury risks. Therefore, the model suggests that the lot size needs to be decreased when the item weight increases.Figure 9
Optimal lot sizesqmat for different weights of each item.In this analysis, the distance value is considered in different cases, such as the distance between the cart and the production line area is longer or shorter for the analysis of picking and storing (lifting and lowering) motions. The distance values have been changed fordc,p between 2.6 and 0.2 m, and for dc,s, they are between 1.4 and 0.2 m. As shown in Figure 10, an increase in the distance between the cart and production area and the distance between the supermarket and the production line causes an increase in the need for rest time and the rest time cost. Furthermore, an increase in the item weight causes an increase in the rest time cost. We also conducted regression analysis for a better understanding of the relationship between the distance between the cart and the supermarket dc,s, the distance between the cart and the production line dc,p, and rest time. The regression analysis showed the relationships between the distance between the cart and the supermarket dc,s and rest time.Figure 10
The rest time cost for different distances of transportationdc,s and dc,p.wmat1=0.1kg; wmat2=0.2kg; wmat3=0.3kg are statistically significant based on p value (0.001); wmat1=0.2kg; wmat2=0.3kg; wmat3=0.4kg are statistically significant based on p value (0.001); and wmat1=0.3kg; wmat2=0.4kg; wmat3=0.5kg are statistically significant based on p value (0.001).The regression analysis showed the relationships between the distance between the cart and the production linedc,p and rest time.wmat1=0.1kg;0.1wmat2=0.2kg; wmat3=0.3kg are statistically significant based on p value (0.001); wmat1=0.2kg; wmat2=0.3kg; wmat3=0.4kg are statistically significant based on p value (0.001); and wmat1=0.3kg; wmat2=0.4kg; wmat3=0.5kg are statistically significant based on p value (0.001).We applied different distance values to analyze the pushing motion between the supermarket area and the production line area. The distance value ofds,p changes between 300 and 5 m. The unit item weight of the materials has been found as wmat1=0.1kg, wmat2=0.2kg, and wmat3=0.3kg, wmat1=0.3kg, wmat2=0.4kg, and wmat3=0.5kg, and wmat1=0.4kg, wmat2=0.5kg, and wmat3=0.6kg.Figure11 illustrates that the increase in the distance of ds,p will increase the rest time and the rest time cost. Furthermore, the increase in the distance ds,p will increase the rest time cost more for heavier items. We also conducted regression analysis for a better understanding of the relationship between the distance between the supermarket and the production line ds,p and rest time. The regression analysis showed the relationships between the distance between the supermarket and the production line ds,p and rest time.Figure 11
The rest time cost for different distances of material handling.wmat1=0.1kg; wmat2=0.2kg; wmat3=0.3kg are statistically significant based on p value (0.001); wmat1=0.3kg; wmat2=0.4kg; wmat3=0.5kg are statistically significant based on p value (0.001); and wmat1=0.4kg; wmat2=0.5kg; wmat3=0.6kg are statistically significant based on p value (0.001).


### 5.3. Discussion

Application of the model shows that using the ergonomic rest time and maximum endurance time in this production supply process decreases the total cost of the process with a 13.05% saving compared to the total cost of the production supply process without the ergonomic rest time and maximum endurance time considerations. The reason is that the ergonomic production supply process model decreases the holding cost and the ergonomic cost of the process. Figure7 shows the total cost curve and additional cost curves for the holding cost, the operational cost, and the rest time cost of material one based on different lot sizes and weight of the items. The savings and results have been calculated in comparison with the total cost of maximum lifting and carrying limit (25 kg) according to ISO-11228 on Ergonomics-Manual Handling, Part 1: Lifting and Carrying, for two-handed lifting where each lot size is equal to 25 kg. The optimal total cost of our model shows that the savings increases or decreases by increase of unit item weights, especially for materials where the weight is under 0,1 kg in the production line (see Figure 6). We investigated the different material requirements in the production line (see Table 4), and the results obtained by application of the new model were consistent with the results from saving analysis in case of different materials where the ergonomic production supply process model decreases the holding cost and the ergonomic cost of the process. It has proved that the model developed is decreasing the cost of the production process and it has provided an opportunity to increase the productivity of the workforce by decreasing the ergonomic risks associated with work-related back disorders.Table 4
The input parameters of different analyses.
ParameterTesting of the modelSaving analysisThe cost analysisHeight of placement analysisOptimal lot size analysisDistance of transportation analysisDistance of material handling analysiswmat1kg0.10.01–30.01–30.01–30.20.40.01–10.10.20.30.10.30.4wmat2kg0.20.01–30.01–30.01–30.30.50.01–10.20.30.40.20.40.5wmat3kg0.30.01–30.01–30.01–30.40.60.01–10.30.40.50.30.50.6wcakg15151515151515151515151515wemkg1111111111111rmat1pcs/prod.3531131111111rmat2pcs/prod.2321121111111rmat3pcs/prod.1111111111111cws/h18181818181818181818181818ch,mat1s/pcs/h0.020.020.020.020.040.020.020.020.020.020.020.020.02ch,mat2s/pcs/h0.040.040.040.040.080.040.040.040.040.040.040.040.04ch,mat3s/pcs/h0.060.060.060.060.10.060.060.060.060.060.060.060.06hst,pm1.21.51.51.51.21.7–0.51.21.21.21.21.21.21.2hpi,cm0.80.50.50.50.81.4–0.20.80.80.80.80.80.80.8hpi,sm1.82221.82–0.81.81.81.81.81.81.81.8hst,cm0.80.80.80.80.81.4–0.20.80.80.80.80.80.80.8hpi,pm1.21.21.21.21.21.7–0.51.21.21.21.21.21.21.2hst,sm1.81.81.81.81.82–0.81.81.81.81.81.81.81.8dc,pm0.50.50.50.50.50.50.52.6–0.22.6–0.22.6–0.20.50.50.5ds,pm100100100100100100100100100100300–5300–5300–5dc,sm0.50.50.50.50.50.50.51.4–0.21.4–0.21.4–0.20.50.50.5As the optimal lot size analysis shows, the increase in the unit weight of the items reduces the optimal lot size values, where the optimal lot size values represent the value that a logistics operator can work continuously without the need for rest time (see Figure9). There are two possible solutions during the realization of logistics operations. The first one is the logistics operator can work without rest time needed with the optimal lot size, and the second one is if the lot size is larger than the optimal lot size; therefore, proper rest time must be applied, which the model suggests. This means that a lot size higher than the optimal lot size increases the need for rest time and that the logistics operator must have necessary rest time after handling the lot to prevent ergonomic risks due to higher load values. An increase in unit item weight increases the contraction of the muscles that leads to increased fatigue and injury risks. Therefore, the model suggests that the lot size needs to be decreased when the item weight increases to eliminate the ergonomic risks due to load on the spine associated with the work of the logistics operator, which will also prevent the medical costs associated with work-related disorders and injuries.As the distance of transportation analysis shows, an increase in the distance between the cart and production line (dc,p) and the distance between the supermarket and the cart (dc,s) causes an increase in the need for rest time and the rest time cost (see Figure 10). Furthermore, an increase in the item weight causes an increase in the rest time cost and rest time need of an operator. This is because the increase in item weight increases the contraction of the muscles that leads to a higher injury risk. Simultaneously, the analysis related to the distance of material handling shows that the increase in the distance between the supermarket and the production line ds,p will increase the rest time and the rest time cost (see Figure 11). Furthermore, the increase in this distance will increase the rest time cost more for heavier items. These analyses could help managers to evaluate different lot sizes according to the distance of transportation and distance of material handling to prevent the ergonomic risks due to long contraction time. It is important to highlight that this method can help in the efficient design of the processes, which decreases the ergonomic risk associated with the distance of transportation and material handling.According to the height of placement analysis, a higher height for placement—wherehst,p represents the height of placement of the production line for storing, hst,c represents the height of placement of the cart for storing, hst,s represents the height of placement of the supermarket for storing, where hpi,c represents the height of placement of the cart for picking, hpi,s represents the height of placement of the supermarket for picking, and hpi,p represents the height of placement of the production line for picking—requires more rest time to prevent ergonomic risks because it may cause a more serious contraction of the muscles when the operator reaches a higher height. Furthermore, the increase of height increases the need for rest time, and it leads to more rest time cost because the operator needs to have rest to recover from a longer contraction. As a result, it can be said that the consideration of this effect can help in the process design to prevent the ergonomic risks associated with the height of placement of materials for improving the worker’s health and productivity.The results of the analyses show that the model, despite its limitations, constitutes a useful guide for the search for the optimum lot sizes for the production line supply process of multi-materials. It improves both social sustainability as it prevents ergonomic risks and economic sustainability as it increases the productivity of the operator, decreases the medical cost associated with work-related back disorders, and decreases the total cost of the process.


### 5.4. Limitations of the Study

In this study, there is a limitation that can be addressed in future research. The study focuses on an average person as a worker for the calculation of rest allowance and maximum endurance time. However, several other factors associated with human performance and health can also have an impact on the rest time needs and maximum endurance time, and they should be addressed in future studies. Further research can investigate human diversity, such as age, gender, and health condition of the worker for the calculation of rest time and maximum endurance time. Besides, the study focused mainly on motions directly associated with the execution of the production line supply tasks by the workers that are picking, storing, and pushing. However, the other motions such as pulling and carrying can also have an impact on human performance and health, and they should be addressed in future studies.


## 5.1. Application of the New Model

We applied a numerical example to investigate the cost effects of rest time, where the total request of the production during the examined time period (480 min) related to the first materialRmat1 is equal to 1440 pcs, that related to the second material Rmat2 is equal to 960 pcs, that related to the third material Rmat3 is equal to 480 pcs, the maximum force load of the human spine Fmax for picking and storing is equal to 2452.5 N, and it has been calculated according to IS0-11228 on Ergonomics-Manual Handling, Part 1: Lifting and Carrying [67], where the limit for two-handed lifting is 25 kg. Fmax for pushing is 4326.2 N as suggested by Resnik and Chaffin [68] for the load limit of 225 kg for four-wheeled carts. The cart weight wca is 15 kg, the time for picking from the supermarket to the warehouse area is 0.11 min for three materials, the time for storing to the cart is 0.05 min for three materials, the time for pushing from the supermarket to the production area is 1.4 min, the time for picking the materials from the cart is 0.07 min, the time for storing the materials to the storage for production line supply is 0.05 min, the time for picking the empties from the storage near the production line is 0.07 min, the time for storing the empties to the cart is 0.05 min, the time for pushing the empties from the production area to the supermarket at the warehouse is 0.7 min, the time for picking the empties from the cart is 0.07 min, and the time for storing the empties to the supermarket is 0.09 min. Furthermore, the examined time period (T) is 480 min, the material need for the first material in the production line rmat1 is 3 pcs/product, the material need for the second material in the production line rmat2 is 2 pcs/product, the material need for the third material in the production line rmat3 is 1 pcs/product, the tact time ttact for the production is 1 min, the unit worker wage cw is 18 $/h, the cost of inventory holding for the first material ch,mat1 is 0.02 $/pcs/h, the cost of inventory holding for the second material ch,mat2 is 0.04 $/pcs/h, the cost of inventory holding for the third material ch,mat3 is 0.06 $/pcs/h, the weight of the first material wmat1 is 0.1 kg, the weight of the second material wmat2 is 0.2 kg, the weight of the third material wmat3 is 0.3 kg, the weight of the empties wem is 1 kg, “μ” is 0.28, and “g” is 9.81 m/s2. The height of placement of the production line for storing hst,p is 1.2 m, the height of placement of the cart for picking hpi,c is 0.8 m, the height of placement of the supermarket for picking hpi,s is 1.8 m, the height of placement of the cart for storing hst,c is 0.8 m, the height of placement of the production line for picking hpi,p is 1.2 m, the height of placement of the supermarket for storing hst,s is 1.8 m, the distance between the cart and the production line dc,p is 0.5 m, the distance between the supermarket and the production line in the production area ds,p is 100 m, and the distance between the cart and the supermarket dc,s is 0.5 m. Anthropometrical parameters have been calculated as follows: “a” and “c” are equal to 0.05 m for picking, storing, and pushing, “b” is equal to 0.5 m for picking and storing, and “d” is equal to 0.35 m for pushing according to the work-related musculoskeletal disorders mentioned in the journal published by the Ministry of Labour and Social Security of Turkey [69].We have used the following formal algorithm to calculate the necessary parameters:Fort = 1 to TIfxmat,t=1 thenFor mat = 1 toMCalculation ofOTmat,tIffMVCmatuni≥0,15 then calculation of RTmat,tpi,s,RTmat,tst,c,RTmat,tpi,c,RTmat,tst,pElseiffMVCmatuni<0,15 then RTmat,tpi,s=RTmat,tst,c=RTmat,tpi,c=RTmat,tst,p=0IffMVCmat,t≥0,15 then calculation of RTmat,tpuElseiffMVCmat,t<0,15 then RTmat,tpu=0Calculation ofRTmat,tIfxem,t=1 thenFor em = 1 toECalculation ofOTem,tIffMVCem≥0,15 then calculation of RTem,tpi,p,RTem,tst,c,RTem,tpi,c,RTem,tst,sElseiffMVCem≥0,15 then RTem,tpi,p=RTem,tst,c=RTem,tpi,c=RTem,tst,s=0IffMVCem,t≥0,15 then calculation of RTem,tpuElseiffMVCem,t≥0,15 then RTem,tpu=0Calculation ofRTem,tCalculation ofOTt and RTtIfMINtmatav≥MAXOTt+RTt then calculation of CtOP, CtRT, and CtHOElseifMINtmatav<MAXOTt+RTt then the situation is not valid, and no solution exists


### 5.1.1. Calculation ofC for the Examination Period T

The applied formal algorithm above can be described withTn=4n2+4n2+n+n+1 function. Therefore, the time complexity of the algorithm based on the Big O principle is On2.We applied our simulation model (see Figure5) with the aforementioned algorithm to calculate the values of the parameters investigated and to find the ergonomically optimal lot sizes by using the solver application of Microsoft Excel for three materials where the total cost is minimum, MINtmatav≥MAXOTt+RTtandfMVC≥0,15. The minimum total cost equals to C=130,36$ and rest time cost equals toCRT=0where the ergonomically optimal lot size for the material one equals toqmat1 27 pcs, the material two equals toqmat2 14 pcs, and the material three equals toqmat3 9 pcs. The material availability of material one tmat1av is equal to 9 min, the material availability of material two tmat2av is equal to 7 min, the material availability of material three tmat3av is equal to 9 min, MINtmatav≥MAXOTt+RTt as MINtmatav=7min, and MAXOTt+RTt=4.39min where RTt=0. RTt=0 means that a logistics operator can work continuously without any need for rest time with these ergonomically optimal lot sizes. We have calculated the total cost of the maximum lifting and carrying limit (25 kg) according to ISO-11228 on Ergonomics-Manual Handling, Part 1: Lifting and Carrying [67], for two-handed lifting where each lot size is equal to 25 kg for the production supply process, and the lot size of material one qmat1 is equal to 240 pcs, the lot size of material two qmat2 is equal to 120 pcs, and the lot size of material three qmat3 is equal to 80 pcs. The total cost is equal to C=149,93$ and CRT=78,95$. CRT=78,95$ is the cost of the need for rest period of logistics operators to prevent any ergonomic risks. In our model, the saving achieved by using the ergonomic rest time and maximum endurance time in comparison to the total cost of the production supply process without ergonomics is equal to 13.05%.Figure 5
The result of a simulation experiment in case of the investigated numerical example in our specific simulator developed in MS Excel.


## 5.1.1. Calculation ofC for the Examination Period T

The applied formal algorithm above can be described withTn=4n2+4n2+n+n+1 function. Therefore, the time complexity of the algorithm based on the Big O principle is On2.We applied our simulation model (see Figure5) with the aforementioned algorithm to calculate the values of the parameters investigated and to find the ergonomically optimal lot sizes by using the solver application of Microsoft Excel for three materials where the total cost is minimum, MINtmatav≥MAXOTt+RTtandfMVC≥0,15. The minimum total cost equals to C=130,36$ and rest time cost equals toCRT=0where the ergonomically optimal lot size for the material one equals toqmat1 27 pcs, the material two equals toqmat2 14 pcs, and the material three equals toqmat3 9 pcs. The material availability of material one tmat1av is equal to 9 min, the material availability of material two tmat2av is equal to 7 min, the material availability of material three tmat3av is equal to 9 min, MINtmatav≥MAXOTt+RTt as MINtmatav=7min, and MAXOTt+RTt=4.39min where RTt=0. RTt=0 means that a logistics operator can work continuously without any need for rest time with these ergonomically optimal lot sizes. We have calculated the total cost of the maximum lifting and carrying limit (25 kg) according to ISO-11228 on Ergonomics-Manual Handling, Part 1: Lifting and Carrying [67], for two-handed lifting where each lot size is equal to 25 kg for the production supply process, and the lot size of material one qmat1 is equal to 240 pcs, the lot size of material two qmat2 is equal to 120 pcs, and the lot size of material three qmat3 is equal to 80 pcs. The total cost is equal to C=149,93$ and CRT=78,95$. CRT=78,95$ is the cost of the need for rest period of logistics operators to prevent any ergonomic risks. In our model, the saving achieved by using the ergonomic rest time and maximum endurance time in comparison to the total cost of the production supply process without ergonomics is equal to 13.05%.Figure 5
The result of a simulation experiment in case of the investigated numerical example in our specific simulator developed in MS Excel.


## 5.2. Parameteric Analysis

In this section, we have applied the mathematical model developed to different scenarios based on different parameter values derived from various industrial practices, and we have analyzed their impact on our new model. In addition, we have changed the weight of the three items from 0.01 to 3 kg in order to observe the effect of different materials. We have calculated the savings obtained by using our new model with differentwmatuni and different material needs in the production line for three different products. The savings have been calculated in comparison with the total cost of maximum lifting and carrying limit (25 kg) according to ISO-11228 on Ergonomics-Manual Handling, Part 1: Lifting and Carrying, for two-handed lifting where each lot size is equal to 25 kg and the optimal total cost of our model. Furthermore, we have calculated the savings where(25)rmat1=3pcs;rmat2=2pcs;rmat3=1pcs;rmat1=5pcs;rmat2=3pcs;rmat3=1pcs;rmat1=rmat2=rmat3=1pcs.Figure6 illustrates that savings increases or decreases by increase of the unit weight of items.. As we have calculated the savings regarding the maximum lifting and carrying limit determined by the ISO-11228 [67], the rest allowance stays constant as it increases each unit weight of materials “mat” but decreases the number of materials handled by causing an increase in the rest time cost and the operational cost.Figure 6
The savings for different weights of each material.We also conducted regression analysis for a better understanding of the relationship between variables. The regression analysis showed the relationship between the weight of handled materials and savings.rmat1=3pcs;rmat2=2pcs;rmat3=1pcs are statistically significant based on p value (0.007); rmat1=5pcs;rmat2=3pcs;rmat3=1pcs are statistically significant based on p value (0.025); and rmat1=rmat2=rmat3=1pcs is statistically significant based on p value (0.001).We have analyzed the material need of the first material in the production linermat1 as 1 pcs/product, the material need of the second material in the production line rmat2 as 1 pcs/product, and the material need of the third material in the production line rmat3 as 1 pcs/product, and, based on these different lot size values, we also analyzed how the total, operational, holding, and rest time costs change. Furthermore, we have applied the inventory holding cost of the first material ch,mat1 as 0.04 $/pcs/h, the inventory holding cost of the second material ch,mat2 as 0.08 $/pcs/h, the inventory holding cost of the third material ch,mat3 as 0.1 $/pcs/h, the weight of the first material wmat1 as 0.2 kg, the weight of the second material wmat2 as 0.3 kg, the weight of the third material wmat3 as 0.4 kg, the weight of the empties wem as 1 kg, “μ” as 0.28, and “g” as 9.81 m/s2, and other parameters were constant. The analysis results are given in Figure 7, where the lot size increases, especially the rest time, the cost of rest allowance, and the holding cost. As can be seen, the optimal solution of the model suggests that the lot size needs to be decreased for a decrease in the total cost, the holding cost, and the rest time cost.Figure 7
The total cost curve and additional cost curves for the holding cost, the operational cost, and the rest time cost of material one based on different lot sizes where the weights of the items are 0.2 kg, 0.3 kg, and 0.4 kg.To analyze the effect of the height of placement of the materials, we have applied different values:hst,p (1.7, 1.6, 1.5, 1.4, 1.3, 1.2, 1.1, 1, 0.9, 0.8, 0.7, 0.6, 0.5).hpi,c (1.4, 1.3, 1.2, 1.1, 1, 0.9, 0.8, 0.7, 0.6, 0.5, 0.4, 0.3, 0.2).hpi,s (2, 1.9, 1.8, 1.7, 1.6, 1.5, 1.4, 1.3, 1.2, 1.1, 1, 0.9, 0.8).hst,c (1.4, 1.3, 1.2, 1.1, 1, 0.9, 0.8, 0.7, 0.6, 0.5, 0.4, 0.3, 0.2).hpi,p (1.7, 1.6, 1.5, 1.4, 1.3, 1.2, 1.1, 1, 0.9, 0.8, 0.7, 0.6, 0.5).hst,s (2, 1.9, 1.8, 1.7, 1.6, 1.5, 1.4, 1.3, 1.2, 1.1, 1, 0.9, 0.8).As can be seen in Figure8, a higher height for placement and picking requires more rest time to prevent ergonomic risks and rest time cost because it may cause a more serious contraction of the muscles when the operator reaches a higher height. Therefore, it increases the need for rest time and the rest time cost as the operator needs to have rest to recover from a longer contraction. We also conducted regression analysis for a better understanding of the relationship between the height of placement for picking and storing and rest time. The regression analysis shows that the relationship between the height of placement for picking of materials hpi,c, hpi,s, hpi,p and rest time is statistically significant based on p value (0.001) and the relationship between the height of placement for storing of materials hst,c, hst,s, hst,p and rest time is statistically significant based on a p value (0.001).Figure 8
The rest time cost for different heights of placement and picking.Furthermore, we have analyzed the effect of the unit item weight on the lot size of materials. As shown in Figure9, the increase in the unit weight of the items reduces the optimal lot size values. The optimal lot size values represent the lot sizes where a logistics operator can work continuously and the need for rest time is equal to zero. This means that a lot size higher than the optimal lot size increases the need for rest time and that the logistics operator must have necessary rest time after handling the lot to prevent ergonomic risks due to higher load values. Higher unit item weight increases the contraction of the muscles that leads to increased fatigue and injury risks. Therefore, the model suggests that the lot size needs to be decreased when the item weight increases.Figure 9
Optimal lot sizesqmat for different weights of each item.In this analysis, the distance value is considered in different cases, such as the distance between the cart and the production line area is longer or shorter for the analysis of picking and storing (lifting and lowering) motions. The distance values have been changed fordc,p between 2.6 and 0.2 m, and for dc,s, they are between 1.4 and 0.2 m. As shown in Figure 10, an increase in the distance between the cart and production area and the distance between the supermarket and the production line causes an increase in the need for rest time and the rest time cost. Furthermore, an increase in the item weight causes an increase in the rest time cost. We also conducted regression analysis for a better understanding of the relationship between the distance between the cart and the supermarket dc,s, the distance between the cart and the production line dc,p, and rest time. The regression analysis showed the relationships between the distance between the cart and the supermarket dc,s and rest time.Figure 10
The rest time cost for different distances of transportationdc,s and dc,p.wmat1=0.1kg; wmat2=0.2kg; wmat3=0.3kg are statistically significant based on p value (0.001); wmat1=0.2kg; wmat2=0.3kg; wmat3=0.4kg are statistically significant based on p value (0.001); and wmat1=0.3kg; wmat2=0.4kg; wmat3=0.5kg are statistically significant based on p value (0.001).The regression analysis showed the relationships between the distance between the cart and the production linedc,p and rest time.wmat1=0.1kg;0.1wmat2=0.2kg; wmat3=0.3kg are statistically significant based on p value (0.001); wmat1=0.2kg; wmat2=0.3kg; wmat3=0.4kg are statistically significant based on p value (0.001); and wmat1=0.3kg; wmat2=0.4kg; wmat3=0.5kg are statistically significant based on p value (0.001).We applied different distance values to analyze the pushing motion between the supermarket area and the production line area. The distance value ofds,p changes between 300 and 5 m. The unit item weight of the materials has been found as wmat1=0.1kg, wmat2=0.2kg, and wmat3=0.3kg, wmat1=0.3kg, wmat2=0.4kg, and wmat3=0.5kg, and wmat1=0.4kg, wmat2=0.5kg, and wmat3=0.6kg.Figure11 illustrates that the increase in the distance of ds,p will increase the rest time and the rest time cost. Furthermore, the increase in the distance ds,p will increase the rest time cost more for heavier items. We also conducted regression analysis for a better understanding of the relationship between the distance between the supermarket and the production line ds,p and rest time. The regression analysis showed the relationships between the distance between the supermarket and the production line ds,p and rest time.Figure 11
The rest time cost for different distances of material handling.wmat1=0.1kg; wmat2=0.2kg; wmat3=0.3kg are statistically significant based on p value (0.001); wmat1=0.3kg; wmat2=0.4kg; wmat3=0.5kg are statistically significant based on p value (0.001); and wmat1=0.4kg; wmat2=0.5kg; wmat3=0.6kg are statistically significant based on p value (0.001).


## 5.3. Discussion

Application of the model shows that using the ergonomic rest time and maximum endurance time in this production supply process decreases the total cost of the process with a 13.05% saving compared to the total cost of the production supply process without the ergonomic rest time and maximum endurance time considerations. The reason is that the ergonomic production supply process model decreases the holding cost and the ergonomic cost of the process. Figure7 shows the total cost curve and additional cost curves for the holding cost, the operational cost, and the rest time cost of material one based on different lot sizes and weight of the items. The savings and results have been calculated in comparison with the total cost of maximum lifting and carrying limit (25 kg) according to ISO-11228 on Ergonomics-Manual Handling, Part 1: Lifting and Carrying, for two-handed lifting where each lot size is equal to 25 kg. The optimal total cost of our model shows that the savings increases or decreases by increase of unit item weights, especially for materials where the weight is under 0,1 kg in the production line (see Figure 6). We investigated the different material requirements in the production line (see Table 4), and the results obtained by application of the new model were consistent with the results from saving analysis in case of different materials where the ergonomic production supply process model decreases the holding cost and the ergonomic cost of the process. It has proved that the model developed is decreasing the cost of the production process and it has provided an opportunity to increase the productivity of the workforce by decreasing the ergonomic risks associated with work-related back disorders.Table 4
The input parameters of different analyses.
ParameterTesting of the modelSaving analysisThe cost analysisHeight of placement analysisOptimal lot size analysisDistance of transportation analysisDistance of material handling analysiswmat1kg0.10.01–30.01–30.01–30.20.40.01–10.10.20.30.10.30.4wmat2kg0.20.01–30.01–30.01–30.30.50.01–10.20.30.40.20.40.5wmat3kg0.30.01–30.01–30.01–30.40.60.01–10.30.40.50.30.50.6wcakg15151515151515151515151515wemkg1111111111111rmat1pcs/prod.3531131111111rmat2pcs/prod.2321121111111rmat3pcs/prod.1111111111111cws/h18181818181818181818181818ch,mat1s/pcs/h0.020.020.020.020.040.020.020.020.020.020.020.020.02ch,mat2s/pcs/h0.040.040.040.040.080.040.040.040.040.040.040.040.04ch,mat3s/pcs/h0.060.060.060.060.10.060.060.060.060.060.060.060.06hst,pm1.21.51.51.51.21.7–0.51.21.21.21.21.21.21.2hpi,cm0.80.50.50.50.81.4–0.20.80.80.80.80.80.80.8hpi,sm1.82221.82–0.81.81.81.81.81.81.81.8hst,cm0.80.80.80.80.81.4–0.20.80.80.80.80.80.80.8hpi,pm1.21.21.21.21.21.7–0.51.21.21.21.21.21.21.2hst,sm1.81.81.81.81.82–0.81.81.81.81.81.81.81.8dc,pm0.50.50.50.50.50.50.52.6–0.22.6–0.22.6–0.20.50.50.5ds,pm100100100100100100100100100100300–5300–5300–5dc,sm0.50.50.50.50.50.50.51.4–0.21.4–0.21.4–0.20.50.50.5As the optimal lot size analysis shows, the increase in the unit weight of the items reduces the optimal lot size values, where the optimal lot size values represent the value that a logistics operator can work continuously without the need for rest time (see Figure9). There are two possible solutions during the realization of logistics operations. The first one is the logistics operator can work without rest time needed with the optimal lot size, and the second one is if the lot size is larger than the optimal lot size; therefore, proper rest time must be applied, which the model suggests. This means that a lot size higher than the optimal lot size increases the need for rest time and that the logistics operator must have necessary rest time after handling the lot to prevent ergonomic risks due to higher load values. An increase in unit item weight increases the contraction of the muscles that leads to increased fatigue and injury risks. Therefore, the model suggests that the lot size needs to be decreased when the item weight increases to eliminate the ergonomic risks due to load on the spine associated with the work of the logistics operator, which will also prevent the medical costs associated with work-related disorders and injuries.As the distance of transportation analysis shows, an increase in the distance between the cart and production line (dc,p) and the distance between the supermarket and the cart (dc,s) causes an increase in the need for rest time and the rest time cost (see Figure 10). Furthermore, an increase in the item weight causes an increase in the rest time cost and rest time need of an operator. This is because the increase in item weight increases the contraction of the muscles that leads to a higher injury risk. Simultaneously, the analysis related to the distance of material handling shows that the increase in the distance between the supermarket and the production line ds,p will increase the rest time and the rest time cost (see Figure 11). Furthermore, the increase in this distance will increase the rest time cost more for heavier items. These analyses could help managers to evaluate different lot sizes according to the distance of transportation and distance of material handling to prevent the ergonomic risks due to long contraction time. It is important to highlight that this method can help in the efficient design of the processes, which decreases the ergonomic risk associated with the distance of transportation and material handling.According to the height of placement analysis, a higher height for placement—wherehst,p represents the height of placement of the production line for storing, hst,c represents the height of placement of the cart for storing, hst,s represents the height of placement of the supermarket for storing, where hpi,c represents the height of placement of the cart for picking, hpi,s represents the height of placement of the supermarket for picking, and hpi,p represents the height of placement of the production line for picking—requires more rest time to prevent ergonomic risks because it may cause a more serious contraction of the muscles when the operator reaches a higher height. Furthermore, the increase of height increases the need for rest time, and it leads to more rest time cost because the operator needs to have rest to recover from a longer contraction. As a result, it can be said that the consideration of this effect can help in the process design to prevent the ergonomic risks associated with the height of placement of materials for improving the worker’s health and productivity.The results of the analyses show that the model, despite its limitations, constitutes a useful guide for the search for the optimum lot sizes for the production line supply process of multi-materials. It improves both social sustainability as it prevents ergonomic risks and economic sustainability as it increases the productivity of the operator, decreases the medical cost associated with work-related back disorders, and decreases the total cost of the process.


## 5.4. Limitations of the Study

In this study, there is a limitation that can be addressed in future research. The study focuses on an average person as a worker for the calculation of rest allowance and maximum endurance time. However, several other factors associated with human performance and health can also have an impact on the rest time needs and maximum endurance time, and they should be addressed in future studies. Further research can investigate human diversity, such as age, gender, and health condition of the worker for the calculation of rest time and maximum endurance time. Besides, the study focused mainly on motions directly associated with the execution of the production line supply tasks by the workers that are picking, storing, and pushing. However, the other motions such as pulling and carrying can also have an impact on human performance and health, and they should be addressed in future studies.


## 6. Conclusion

Inventory management and lot sizing play a critical role in business success. Various kinds of research have been done to find the optimum plan or lot size to reach sustainable success. Especially, economic and environmental sustainability aspects have been investigated extensively and integrated by many researchers. However, there are only a few contributions regarding human health and social sustainability. This paper demonstrates an extended economic and ergonomic multi-material inventory model to investigate this research gap. The model developed integrates the rest allowance to reduce ergonomic risks and work-related disorders of the back. The rest allowance is determined by endurance time, contraction time, and relative force based on the weight and number of the items. The rest time is the time required to prevent the risks of fatigue associated with a load value higher than the optimal lot size. The new model aims to minimize both the total cost of the production line supply process and ergonomic risks. The case study based on our model presented in this paper investigates the supply of three different raw materials with one operator to the production line with a simple cart as a single-operator-multi-material model. The efficiency of the model has been investigated with the simulation and parameteric analysis through different datasets. The analyses have proved that the model developed is suitable for different practices, with the constituted rest time provided with an opportunity to increase the productivity of the workforce and to reduce the cost of a company. Furthermore, the analyses have shown that a higher distance of movement and height of placement for manual material handling will cause a higher need for rest time. The relevant approach can help managers to evaluate lot sizes with the distance of movement (production line and supermarket) and the height of placement and the effects of picking to prevent ergonomic risks. In addition, it will decrease high costs for work-related musculoskeletal disorders, work-related disorders of the back, absenteeism, productivity loss, increased health care, disabilities, and workers’ compensation costs. Our model also can help managers to maintain the decision-making process based on a schedule of work-rest time. The results of this study have demonstrated the viability of the inventory model optimization with multiple objectives and complex constraints. With regard to the maximum endurance time, the pushing motion and the multi-material production line supply in an inventory model have been considered for the first time. Most notably, this study helps to interpret the optimal ergonomic lot size by using special aspects for multi-materials that accommodate social and economic sustainability in the production line supply.Furthermore, future research can investigate a scenario handling many materials by multi-operator (MMMO) and single material by multi-operator (SMMO) for a better understanding of the relationship between lot sizing and ergonomics in the production line supply and can investigate how to improve the ergonomic and economic conditions in industrial practices. In addition, future research can investigate a lot-sizing model with ergonomics for different manual material handling tasks, and it would be valuable to extend this study with the integration of other motions such as carrying and pulling.



---

*Source: 9957299-2021-08-29.xml*