SciELO - Scientific Electronic Library Online

 
vol.90 número227Mercury and human health: probabilistic risk characterization in one of the oldest gold mining areas in EcuadorThe evolution and trends of hyperconvergence in the telecommunications sector: a competitive intelligence review índice de autoresíndice de assuntospesquisa de artigos
Home Pagelista alfabética de periódicos  

Serviços Personalizados

Journal

Artigo

Indicadores

Links relacionados

  • Em processo de indexaçãoCitado por Google
  • Não possue artigos similaresSimilares em SciELO
  • Em processo de indexaçãoSimilares em Google

Compartilhar


DYNA

versão impressa ISSN 0012-7353versão On-line ISSN 2346-2183

Dyna rev.fac.nac.minas vol.90 no.227 Medellín jul./set. 2023  Epub 19-Mar-2024

https://doi.org/10.15446/dyna.v90n227.107105 

Article

Occupational hazards and economic indicators in the scheduling of a make-to-order system

Peligros ocupacionales e indicadores económicos en la programación de un sistema por pedido

Germán Coca-Ortegóna 
http://orcid.org/0000-0002-6760-4111

aUniversidad EIA, Facultad de Ingeniería, Medellín, Colombia. german.coca@eia.edu.co


Abstract

This paper examines some specific occupational hazards and certain economic indicators of sustainability in a make-to-order manufacturing system. In this respect, two multi-objective protocols are structured. The first is a redesigned genetic algorithm; the second is based on normalized summations. The comparison establishes that the former outperforms the latter by 528%, concerning to the "C-Metric" indicator.

Keywords: multi-objective; sustainability; occupational hazards; economic indicators; C-Metric indicator

Resumen

Este artículo examina algunos peligros ocupacionales específicos y ciertos indicadores económicos de sostenibilidad en un sistema de manufactura por pedido. Al respecto, se estructuran dos métodos multiobjetivo. El primero es un algoritmo genético rediseñado; el segundo está basado en sumas normalizadas. La comparación establece que el primero supera al segundo en 528%, en lo concerniente al indicador “C-Metric”.

Palabras clave: multiobjetivo; sostenibilidad; peligros ocupacionales; indicadores económicos; indicador C-Metric

1. Introduction

In this contribution, the manufacturing operations of an organization engaged in producing "non-standardized products" in the metal-mechanic sector are examined. For this entity, its occupational hazards and manufacturing competitive advantages were analyzed. This way, the following factors were identified as priorities: "non-conforming" ratio generated by worker fatigue, sound intensity, lighting level penalty, processing time and direct personnel costs.

In particular, from the social dimension of sustainability, the evaluation of "occupational hazards" is a fundamental element that aims to maintain and improve people's physical and emotional health in the workplace [1-3]. In this sense, what has been previously explained is related to the specific content of the eighth sustainable development goal in the 17 goals designed in this area by the United Nations Organization and whose compliance goals are defined for the year 2030 [4]. This goal points out, among other aspects, the need to ensure people's health and the protection of healthcare systems in the workplace.

Based on the above, regarding the "nonconforming" proportion generated by worker fatigue, it is mentioned that the company is intensive in the use of labor, so the probability of obtaining defective products increases as the working day goes by [5-6]. Likewise, the company detected that fatigue in the workstations comes from exposure to handling heavy materials, vibrations and ambient temperature.

Likewise, concerning the sound intensity variable, it is stated that the company uses equipment with the capacity to emit sound intensity levels close to or higher than 90 decibels. Therefore, when the optimum process time is identified, this intensity level will increase because several pieces of equipment will operate simultaneously [7].

In turn, in the case of the lighting level penalty, the aim is to detect those work organization programs close to the optimum process time, whose implementation allows the tasks to be executed in the presence of natural lighting or, failing that, under normal operating conditions of the artificial lighting system.

Regarding the processing time, the company establishes it using dispatch rules (SOT, EDD, FCFS, LCFS, LOT). In this regard, implementing these rules causes a low level of productivity in the manufacturing systems. [8-9].

Regarding direct personnel costs, the entity's highly variable demand causes inconveniences in establishing the organizational form of the operations system associated with the lowest cost.

Globally, it is aimed to define the planning horizon of the production system, able to meet the following conditions: minimize labor risk (minimum exposure to vibration, weight handling and ambient temperature, which in turn results in obtaining a lower proportion of non-conforming production; the minimum level of sound intensity emitted; and lower level of lighting penalties), achieve the due delivery opportunity in the manufacturing process (identification of the minimum process time) and also, detect the minimum personnel cost.

By evaluating the above, it is stated that this study develops the process of multi-objective sustainability reduction of the following social and economic variables: "non-conforming" proportion generated by worker fatigue, sound intensity, lighting level penalties at the work site, processing time and direct personnel costs.

In addition, it is noted that the detection of the minimum value corresponding to the "non-compliant" fraction caused by worker fatigue will lead to the control of exposure conditions to weight handling, vibrations and ambient temperature.

Regarding the observed variables, it is necessary to establish the minimum processing time in the first step. Consequently, the parameter estimation in question, since there are "JT" jobs and "PC" processing centers, implies the analysis of JT!PC solutions [10-12]. Therefore, the difficulty arises to explore with the necessary timeliness JT!PC feasible points [13-6]. Thus, an algorithm structured from "multi-objective genetic algorithms" technique concepts [17-19] is developed in this paper. Thus, the designed algorithm would allow the simultaneous optimization of the selected factors once the processing time is established.

Similarly, some developments of the technique called "genetic algorithms" are illustrated below. For this purpose, the following "classes" are defined. Operator redefinition in the production of "non-standardized goods," determination of algorithms in the production of "non-standardized goods," sustainability analysis with a primary focus on energy consumption evaluation, and, finally, multi-objective evaluation in the production of "non-standardized goods." These "classes" are shown below:

Operators in the production of "non-standardized goods" are redefined. This category includes three documents. The papers seek to ensure the highest dominance of solutions in the new generations. The above is achieved by restructuring the crossover and mutation operations [20-21-22].

Determination of algorithms in the production of "non-standardized goods." Four contributions are recorded for this purpose, whose perspective seeks parallel minimization of processing time and direct personnel costs [23-26].

Sustainability analysis with a primary focus on the evaluation of energy consumption. Eleven articles are identified in this case that considers the exclusive study of energy consumption as a criterion for sustainability analysis [13,27-36]. The above documents do not examine other environmental aspects or other relevant factors of a social nature, as is the particular case of occupational health and safety.

Multi-objective evaluation in the production of "non-standardized goods." In this class, ten contributions are identified. These contributions can be categorized into exact techniques or heuristic techniques. Concerning the perspective of exact techniques, two articles oriented to examine multi-objective sequencing in "Job Shop" were detected. Thus, the first of these articles incorporates the following methods in its analysis: Monte Carlo simulation, Markov chains, and waiting lines [14], while the second article uses the Branch and Bound technique [37].

Exact technique implementation evolved towards the design of heuristic techniques due to the existing limitations in treating complex problems. Accordingly, the other eight contributions correspond to the field of heuristics. In particular, three papers deal with the structuring of multi-objective algorithms. In turn, four other developments propose the solution of specific problems through the design of hybrid algorithms. Finally, the last contribution exposes the trends of multi-objective sequencing in obtaining "non-standardized" goods.

The review of the three articles based on the definition of multi-objective algorithms shows that they focus on analyzing certain parameters not reported in the literature in relation to some relevant research situations.

Thus, in the first application, an algorithm is implemented in a dynamic environment (random occurrence of adverse events) to minimize the manufacturing system's processing time and instability [38]. Regarding the second development, an algorithm that improves the protocol's performance defined in NGSA-II is structured [39]. Regarding the third one, an algorithm that minimizes both the processing time and the deterioration of the machines according to their utilization level in real time is defined [40].

The following four developments correspond to a series of hybrid protocols that aim to detect the minimum values of the variables discussed below: processing time, total working weight, and maximum working weight [41-44].

In addition, the eighth development synthesizes, among others, the evolution of the production-sequencing problem of "non-standardized" goods and obtains the projection of its trends [45].

The literature review shows that the simultaneous analysis of the variables examined in this study has been carried out superficially in the contributions evaluated. Similarly, the issue of occupational health and safety has been scarcely studied in depth in the sequencing of manufacturing processes at a multi-objective level [46-48].

Additionally, it is stated that the organization under review will also achieve other benefits, such as a reduction of the level of noise emitted to the community of influence. In Addition, attenuation of the proportion of people affected by occupational hearing loss in the regional and national health systems since the sound intensity coming from the operation of machines and equipment will be controlled. Reduction of metal waste and carbon dioxide emissions generated during the reprocessing of "nonconforming" products in the manufacturing system; and reduction of the number of people affected by musculoskeletal diseases since the level of vibrations, ambient temperature, and weight handling will be reviewed.

Finally, this paper will set out the materials and methods; later, the results will be recorded, followed by their discussion, and, finally, the conclusions will be structured. Next, the previously mentioned stages will be presented.

2. Materials and methods

In this section, steps 1 to 6 describe the operation mode of the redesigned genetic algorithm (first method). In contrast, in step 7, the operation mode of the algorithm corresponding to the execution of weighted sums (second method) is shown. Similarly, step 8 compares the results obtained by executing both methods.

Step 1 Obtaining start-up population

Based on the characteristics of the evaluated company, n orders to be examined are observed. These orders by exchanging positions constitute n! solutions. Each solution includes six subsolutions. The segments of Table 1 represent a subsolution: "Solution," "work rotation," "working hours," and "works 6 or 7 days a week".

Thus, in the "solution" segment, the individual is illustrated as an example: "6-3-8-9-9-1-4-7-5-2". This solution comes from one of 9! possible permutations (nine working orders in total: from the first order in the example, which corresponds to order 6, to the last order in the same example, i.e., order 2). In the individual named above, it is pointed out that "order 6" would be equivalent to the first order to be transformed; thus, "order 2" would correspond to the order processed in ninth place. Six "ways" of structuring the work shifts are identified for each solution. These "forms" can be seen in the segments of Table 1: "Work rotations" and "working hours." According to the data in these segments, it is defined that a sollows: one 8-hour shift, two 8-hour shifts or three 8-hour shiflution can be scheduled in three types of shifts each day, as fots (first shift: 6 a.m. to 2 p.m.; second shift: 2 p.m. to 10 p.m.; third shift: 10 p.m. to 6 a.m.). Similarly, in the segment: "Work 6 or 7 days a week", it is indicated that the solution can be analyzed in two work schedules: up to Saturday (6 days) or with a Sunday holiday (7 days). Table 1 is recorded below.

Each sub-solution recorded in Table 1 will be treated from this point of the document as a specific solution. Consequently, the scheme concerning the organization of a solution is recorded in Table 2. In the scheme, 11 segments are observed. In new of them, the work orders are arranged. For the case, "segment 1" is associated with "order 6", representing that such "order" would be scheduled first. When the segment number is increased to "segment 9", it is established that "order 2" is assigned to this segment; therefore, the order above will enter the production plant in ninth place. Likewise, "segment 10" shows the number of hours worked each day (in the example, 8 hours: starts at 6 a.m. and ends at 2 p.m.). Finally, "segment 11" refers to the number of days scheduled per week (7 for the example, since it includes the term "with Sunday").

According to Table 2, "S" solutions are randomly generated.

Table 1 Ways to structure work shifts. 

Sub-solution number Solution Job rotations Workdayhours Work 6 or 7 days a week
1 6-3-8-9-1-4-7-5-2 6 a 14 8 Up to Saturday
2 6-3-8-9-1-4-7-5-2 6 a 14 - 14 a 22 16 Up to Saturday
3 6-3-8-9-1-4-7-5-2 6 a 14 - 14 a 22 - 22 a 6 24 Up to Saturday
4 6-3-8-9-1-4-7-5-2 6 a 14 8 Inc Sunday
5 6-3-8-9-1-4-7-5-2 6 a 14 - 14 a 22 16 Inc Sunday
6 6-3-8-9-1-4-7-5-2 6 a 14 - 14 a 22 - 22 a 6 24 Incl Sunday

Source: prepared by the authors.

Table 2 Outline and implementation of a solution. 

Source: prepared by the authors.

Step 2 Group definition

Size A population is divided into five groups. Each group will allow the evaluation of one objective and will be composed of A/5 solutions. The solutions per group are analyzed and sorted in ascending order. Subsequently, 30% of the solutions with a higher value are removed per group.

Then the processing time of each group is established. The variable in question is estimated from the application of eq. (1):

In eq. (1): "p" = the number of orders; "q" = the number of equipment; "r": index varying from 1 to total orders; "t": index fluctuating from 1 to total equipment; "Ost" = operation time of order "p" on equipment "q"; and FitnessMakespan = processing time of the evaluated orders.

The groups are then examined with estimation techniques as follows:

Estimation technique for the "non-conforming" proportion. This technique is based on the data processing corresponding to the "non-conforming" fraction established by the organization examined. Thus, the "non-conforming" proportion allocated over the scheduling horizon, given that there are orders in process at the hourly level, is totaled and multiplied by 100.

Sound intensity estimation technique. The first step to implementing the method under consideration is identifying the sound intensity in each working hour. Then, its percentage participation to 90 decibels is determined. This participation is recorded in each of the hours that make up the production horizon. Then, the sum of the participation previously mentioned is added up. Finally, the previous sum is related to the projected number of shifts.

Lighting level estimation technique. The hourly Lighting level in a particular work area is first established in this case. Later, the degree of hourly compliance is estimated as far as lighting level is concerned. This degree of compliance is determined by finding the mathematical equivalence between the actual and required lighting levels in the evaluated area. For example, if the lighting needs of such an area amount to 2000 lumens and the direct measurement yields a result of 500 lumens, a degree of compliance of 25% would be achieved. Subsequently, the penalty level per hour is obtained, which is given by calculating the difference between the ideal goal of 100% and the degree of compliance defined for the area under analysis. Finally, the penalties identified for each of the hours that make up the projected production plan are consolidated. Thus, the total score to be examined is quantified.

Technique for estimating direct personnel costs. This estimation is obtained from the following work schedules: 6 a.m. to 10 p.m. (16 hours for the first and second shifts combined) or 10 p.m. to 6 a.m. (8 hours for the third shift) on working days or Sundays. Each production schedule's hours, depending on the day (working or Sunday), is assigned a specific score. This score comes from the percentage overpayment concerning the standard working day wage base. Consequently, the sum of the score defined for each hour and work center observed in the projected production plan is added.

Step 3 Establish the first population structure

The first population structure is integrated: each group of 0.7A/5 solutions is divided into five segments. From each of the segments, 0.02A solutions are randomly removed. Together the groups provide 0.5A solutions to the first population. These 0.5A solutions are affected with a crossover ratio of 97%. The task is run until an additional amount corresponding to 42% of the population is generated. Such crossover is performed by exchanging those genes that include the work orders between solutions. In turn, 7% of the solutions are determined by applying the mutation operation with a 12% probability to 90% of the population of size A already created. In addition, the remaining 1% of the solutions are extracted from the total number of possible solutions to strengthen diversity. Additionally, it is indicated that to execute the operations of the genetic algorithms, crossover, and mutation, the segments (genes) related to the hours or days projected respectively per day or week are not considered.

Step 4 Determine the Pareto structure in the defined population

The solutions of the first population are evaluated in the various objective functions; then, the estimated values for each solution are compared with those determined for the other solutions in the population in question. Consequently, the "dominant" solutions are segregated, i.e., the solutions whose values are strictly lower in the various objective functions or, failing that, which, being equal in some of them, have lower values in the other functions. Concerning the above, the first group of solutions obtained constitutes the first Pareto front. Then, the previously explained process is applied to the dominant solutions of the population until a new level of best individuals is identified each time.

Step 5 Establish the second population structure

The second population structure is made up as follows: 0.5A solutions come from the first Pareto fronts established in step 4; for this purpose, the level of dominance corresponding to each front is used as the order of priority assignment. In addition, 0.3A solutions originate from the random crossing of the individuals that make up the first three fronts; if there were only two or a single front, the crossing operation would be carried out between the solutions of the existing fronts. In addition, 0.19A solutions come from the random mutation of the 0.8A individuals previously obtained. To close 0.01A solutions, they are taken randomly from the complete set of solutions to ensure population divergence.

Step 6. Define the Pareto structure of the population identified in step 5

The Pareto structure of the population identified in step 5 is determined as recorded in step 4.

Step 7. Performance evaluation of a linear function by solution (second method)

The activities described below are carried out in each program cycle. Randomly generate A solutions; estimate the objective functions for each solution generated; establish the performance of the linear function per solution; and, finally, choose the solution that corresponds to the highest value of such linear function. Then, consolidate the selected solutions into a single population; re-estimate the performance of the linear function per solution for the chosen solutions; and, finally, identify the best-performing solution, i.e. the solution with the highest performance in terms of the linear function applied.

Eq. (2) shows the form of estimation relative to the weighted valuation of each solution in each function analyzed.

In eq. (2): Fd= Each function analyzed; Fdmxm = upper value in the set of response solutions per function examined; Fdmnm = lower value in the set of response solutions per function evaluated; Fd(Zn) = response value of each Zn solution per function; Zn = number of revised solution in each function; n = number of solutions in each response population (n varies from 1 to the total number of response solutions in each population); WV(Fd) = weighted value of each function (its determination comes from the decision of the work team in the company, being able to fluctuate between 0 and 100%); WT(Fd) = total weighted factors of the functions (the consolidated is 1); TF = total functions to be evaluated; d = index varying from 1 to the total functions to be examined (TF); FWd(Zn) = weighted valuation of each solution in each function analyzed.

After estimating FWd(Zn), the total weighted valuation of each solution for the set of evaluated solutions is established in eq. (3).

In eq. (3): FWT (Zn) = weighted total rating of each solution for the set of evaluated functions.

Step 8 Check the designed methods

The analysis of the populations identified according to the proposed methods is carried out using the value corresponding to the C-Metric indicator (dominance indicator when comparing the results of the multi-objective methods) [49-50-51] as a comparison parameter [49-51].

3. Results

The populations determined through the application of the methods designed and the methods in question are then analyzed.

Population established according to redesigned genetic algorithm implementation (first method)

As described in the section "materials and methods", Table 3 shows the results of two solutions concerning the 220 specimens that make up the response population structure. The values corresponding to the average and variation coefficient for each variable examined are also shown for this population. The estimates are recorded in points in Tables 3, 4, 5 and 6.

According to the contents of Table 3, the following characterization in the number of solutions in the Pareto structure is determined: 71 on the first front, 63 on the second front, 36 on the third front, 29 on the fourth front and 21 on the fifth front.

The population is determined according to method implementation based on a linear function by weighted sums (second method)

Table 4 reports the results from step 7 of the "materials and methods" section. Two solutions of the 220 solutions detected in the final response population are recorded in this regard. Likewise, the values concerning the average and the variation coefficient per variable are recorded for the latter population.

Table 3 Results summary (first method). 

Position solution in the population Makespan (hours) Makespan (weeks) Estimate costs Estimated ratio Estimate sound Estimated lighting
8 210 2,18 0 1720 842 1522
137 218 1,51 2015 1985 881 1738
Average 248,70 2,27 1408,21 1832,26 786,26 1758,23
Coefficient Variation 0,12 0,63 1,04 0,23 0,47 0,42

Source: prepared by the authors.

Table 4 Results summary (Second method). 

Position solution in the population Makespan (hours) Makespan (weeks) Estimate costs Estimated ratio Estimate sound Estimated Lighting
6 228 2,38 0 1524 773 1396
205 195 2,03 0 1650 846 1462
Average 203,44 3,08 0 1628,37 828,25 1478,23
Coefficient Variation 0,042 0,32 0/0 0,15 0,18 0,17

Source: prepared by the authors.

Likewise, the dominance analysis of the second response population allows us to conclude that the solutions as a whole constitute a single Pareto front.

Comparison of proposed methods

When evaluating the information recorded in Tables 5 and 6, it is determined that the values examined (coefficient of variation and maximum value) are higher when comparing "method one" with "method two". Due to the fact described above, it is pointed out that the modified genetic algorithm shows better ability concerning the detection of the solutions dispersion that integrates the whole set of possible answers. Consequently, the C-Metric indicator corresponding to "method 1" amounts to 88.57 points, i.e. 88.57% of its "non-dominated" solutions (62 out of 71 solutions on the first Pareto front) retain this characteristic after comparison with the "dominant" solutions found using "method 2" (220 solutions on the optimal front). Thus, it is noted that, for the indicator in question, the "modified genetic algorithm" outperforms the "weighted summation" method by 528%.

Based on the above, the examined company decided to implement the "first method" in its production system's scheduling. Thus, for the solutions belonging to the first Pareto front, a specific grouping of data by process time (in weeks) is established in consensus with the company's management team. This grouping is shown in Table 7.

The company must determine the market demand trend for the time of the year under analysis by studying the data recorded in Table 7. This way, it will proceed to segregate that individual whose processing time corresponds to the shortest time within the respective interval. Subsequently, the company will define and apply decisions related to the other variables of interest: direct personnel costs, "non-conforming" ratio, sound intensity, and lighting level. Table 8 shows some of the general actions that the company could take depending on the individual chosen for each specific interval.

Table 5 Dispersion and C-Metric (processing time and direct personnel costs). 

Source: prepared by the authors

Table 6 Dispersion ("non-compliant", sound and lighting). 

Source: prepared by the authors.

Table 7 Systematization of first front information (method one). 

Number Segment Lower elevation (Time weeks) Upper elevation (Time weeks) Quantity solutions Percentage (%) Demand trends
1 0,90 1,95 14 19,72 In constant growth
2 1,95 3,15 37 52,11 Constant with periods of growth
3 3,15 3,85 08 11,27 Constant
4 3,85 4,65 07 09,86 Constant with decreasing periods
5 4,65 5,85 05 07,04 Decreasing with normalization periods.

Source: prepared by the authors.

Table 8 Description of general actions to take. 

Number Segment Performance overall process time Overall description of performance of the other variables examined General actions to be taken
1 Very small Very high level of risks and direct personnel costs. Coordinate the use of sound isolation methods at the source. Evaluate the possibility of using another solution that, in the same interval, makes better use of natural lighting. Adjust the number of lamps according to the specific lighting needs of the workstations. Review the operating conditions of the workstations to reduce exposure to high ambient temperature, handled weight and vibrations.
2 Small High intensity level, in terms of the behavior of the variables recorded in segment 1. Actions to implement would be equivalent to the actions recorded in interval 1. In this regard, it is clarified that the levels of investment required to implement the solution chosen in interval two would be lower than the levels of investment required to apply the individual chosen in interval one since, in this second interval, less exposure to the occupational hazards evaluated is generated.
3 Middle Middle intensity level, around the performance of the variables examined. The actions to be executed are similar to those recorded for intervals 1 and 2. Likewise, the levels of investment required for this case are lower than the levels required for applying the solutions chosen in the intervals mentioned above.
4 High Low level of intensity of the variables evaluated. In this interval, the exposure to occupational hazards would be significantly reduced since the processing time is high. However, given the decrease in demand, the company, in order to seek a balance in terms of sustainability, should project its products to new markets.
5 Very high Very low level of intensity of the variables analyzed. Exposure to occupational hazards would decrease dramatically in the present interval since the longer process time would be implemented. In turn, the company should develop new products to target either current or new markets.

Source: prepared by the authors.

4. Results discussion

Independently of the fact that in the literature reviewed [9,10,13,28,41-43,45,46,52] and its theoretical support, research developments in manufacturing processes concerning the production of "non-standardized" goods using techniques based on metaheuristic methods have been widely extended since the late 1990s. It is noted that the variables associated with certain occupational hazards (high sound level, poor lighting and factors affecting the behavior of the proportion of "non-conforming" production -exposure to vibrations, weight handled and ambient temperature-) and, also, with certain economic indicators (processing time and direct personnel costs), have not been assessed simultaneously, nor have they been analyzed in the required depth: vibrations, weight handled and ambient temperature). Moreover, likewise, certain indicators of an economic nature (processing time and direct personnel costs) have not been simultaneously assessed, nor have they been analyzed with the required depth. Similarly, the problem area is of vital importance since it allows the application of sustainability in the operations of manufacturing organizations, based on two complementary theoretical references on the one hand, the management of occupational hazards and, on the other hand, the manufacturing missions (flexibility and costs). Consequently, it seeks to ensure the projection of the organization in a global market that requires entities to act socially responsibly.

When considering two critical evaluation indicators, it is found that the first method outperforms the weighted sum method (second method). Thus, the first of the indicators mentioned corresponds to the C-Metric indicator. In this case, it is found that the value of this parameter is higher by 528% when the respective comparison of the "redesigned genetic algorithm" is made about the "weighted sums" method. Consequently, it is stated that, between the two protocols examined, the first method is more skillful in detecting dominant solutions. In turn, the second indicator refers to the higher number of fronts established through the application of the first method, thus evidencing its greater ability to represent diverse solutions in a Pareto optimal structure.

Likewise, it is also detected that the minimum values in the several variables examined are lower when method one is compared with the weighted summation method, as follows:

By 10.5% as regards processing time (expressed in hours), by 94.44% as regards processing time (expressed in weeks), by 4.21% as relates to the proportion of "non-compliant" production, by 13% as regards the sound intensity variable and, finally, by 11.13% as regards the lighting level penalty variable.

It is also expressed that the company should evaluate market trends, selecting from the interval of interest the individual with the shortest manufacturing time. However, it is necessary to analyze whether implementing the previous individual generates occupational risk situations in the production system, whose impact could be potentially harmful to people's health. Based on the above, it is indicated that the organization must take appropriate measures, given the identification of the risk above situations. Thus, the above actions should focus on noise, the isolation of the source, and lighting, the installation of efficient systems, under which natural lighting is prioritized as far as possible.

5. Conclusions

Based on what has been defined in the previous sections, it is concluded that: 1) the structured method shows relevant performance characteristics such as a significant number of dominant solutions, a representative amount of Pareto fronts, and detection of minimum values in the different variables of interest. 2) The detailed exposition of the steps related to the designed methods facilitates their replication by the academic community interested in this topic. 3) Given that a certain value of process time is selected, to adapt the organization to the predominant market trends, it is necessary to evaluate the possible effects on workers due to the potential increase in exposure to certain occupational hazards. Thus, for those cases required, the corresponding preventive measures must be defined and implemented. 4) The decision to give greater importance to a certain value of process time is related to the behavior of the economic system; therefore, it is essential to complement this decision by incorporating relevant criteria of a social nature.

In addition, the need to integrate relevant concepts of uncertainty and some relevant trends of dynamic programming at the multi-objective level into the present study is also established. In this way, a robust model capable of holistically explaining the behavior of the operations subsystem, both in the context of the organization and the context of the market, would be structured under a sustainability perspective.

References

[1] Hale, J., Legun, J., Campbell, H., and Carolan, M., Social sustainability indicators as performance. Geoforum, 103, pp. 47-55, 2019. DOI: https://doi.org/10.1016/j.geoforum.2019.03.008Links ]

[2] Hutchins, M., Richter, J., Henry, M. and Sutherland, J., Development of indicators for the social dimension of sustainability in a U.S. business context. Journal of Cleaner Production, 212, pp. 687-697, 2019. DOI: https://doi.org/10.1016/j.jclepro.2018.11.199Links ]

[3] Khalid, A., Khan, Z., Idrees, M., Kirisci, P., Ghrairi, Z., Thoben, K. and Pannek, J., Understanding vulnerabilities in cyber physical production systems, International Journal of Computer Integrated Manufacturing, 35(6), pp. 569-582, 2021. DOI: https://doi.org/10.1080/0951192X.2021.1992656Links ]

[4] Naciones Unidas, La Agenda 2030 y los objetivos de desarrollo sostenible: una oportunidad para América Latina y el Caribe LC/G. 2681-P/Rev, [online]. 2018. Available at: https://www.un.org/sustainabledevelopment/es/objetivos-de-desarrollo-sostenible/Links ]

[5] Berti, N., Artigues, C., Battaïa, O., Guillaume, R. and Battini, D., Heuristic approaches for scheduling manufacturing tasks while talking into account accumulated human fatigue. IFAC-PapersOnLine, 52(13), pp. 963-968, 2019. DOI: https://doi.org/10.1016/j.ifacol.2019.11.319Links ]

[6] Yung, M., Kolus, A., Wells, R., and Neumann, W.P., Examining the fatigue-quality relationship in manufacturing. Applied Ergonomics, 82, art. 102919, 2020. DOI: https://doi.org/10.1016/j.apergo.2019.102919Links ]

[7] Kempen, E., Casas, M., Pershagen, G. and Foraster, M., Who environmental noise guidelines for the European Region: a systematic review on environmental noise and cardiovascular and metabolic effects. A summary. International Journal of Environmental Research and Public Health, 15(2), art. 5020379, 2018. DOI: https://doi.org/10.3390/ijerph15020379Links ]

[8] Amiri, F., Shirazi, B. and Tajdin, A., Multi-objective simulation optimization for uncertain resource assignment and job sequence in automated flexible job shop. Applied Soft Computing Journal, 75, pp. 190-202, 2019. DOI: https://doi.org/10.1016/j.asoc.2018.11.015Links ]

[9] Gu, X., Huang, M. and Liang, X., A discrete particle swarm optimization algorithm with adaptive inertia weight for solving multiobjective flexible job-shop scheduling problem. IEEE Access, 8, pp. 33125-33136, 2020. DOI: https://doi.org/10.1109/ACCESS.2020.2974014Links ]

[10] Bissoli, D., Zufferey, N. and Amaral, A., Lexicographic optimization-based clusteringsearch metaheuristic for the multiobjective flexible job shop scheduling problem. International Transactions in Operational Research, 28(3), pp. 1-26, 2019. DOI: https://doi.org/10.1111/itor.12745Links ]

[11] Rahmati, S., Zandieh, M., and Yazdani, M., Developing two multi-objective evolutionary algorithms for the multi-objective flexible job shop scheduling problem. The International Journal of Advanced Manufacturing Technology, 64(5-8), pp. 915-932, 2012. DOI: http://doi.org/10.1007/s00170-012-4051-1.Links ]

[12] Zhang, F., Bail, J., Yang D., and Wang, O., Digital twin data-driven proactive job-shop scheduling strategy towards asymmetric manufacturing execution decision. Scientific Reports, 12, art. 1546, 2022. DOI: https://doi.org/10.1038/s41598-022-05304-wLinks ]

[13] Guo, Y., Huang, M., Wang, Q. and León, V., Single-machine rework rescheduling to minimize total waiting time with fixed sequence of jobs and release times, in: IEEE Access, 9, pp. 1205-1218, 2021. DOI: https://doi.org/10.1109/ACCESS.2019.2957132Links ]

[14] Nicoara, E., Multi-objective flexible job sho scheduling optimization models. Economic Inshights - Trends and Challenges, 49(2), pp. 79-86, 2015. Available at: https://upg-bulletin-se.ro/old_site/archive/2012-2/7.%20Nicoara.pdfLinks ]

[15] Ozturk, G., Bahadir, O., and Teymourifar, A., Extracting priority rules for dynamic multiobjective flexible job shop scheduling problems using gene expression programming. International Journal of Production Research, 57(19), pp. 3121-3137, 2019. DOI: https://doi.org/10.1080/00207543.2018.1543964Links ]

[16] Valenzuela, V., Cosío, M., and Romero, A., A cooperative coevolutionary algorithm approach to the no-wait job shop scheduling problem. Expert Systems With Applications, 194(15), art. 116498, 2022. DOI: https://doi.org/10.1016/j.eswa.2022.116498Links ]

[17] Deb, K., and Jain, H., An evolutionary many-objective optimization algorithm using reference point-based nondominated sorting approach, Part I: Solving problems with box constraints. IEEE Transactions on Evolutionary Computation, 18(4), pp. 577-601, 2014. DOI: https://doi.org/10.1109/TEVC.2013.2281535Links ]

[18] Deb, K., Pratap, A., Agarwal, S., and Meyarivan, T., A fast and elitist multiobjective genetic algorithm: NSGA-II. IEEE Transactions on Evolutionary Computation, 6(2), pp. 182-197, 2002. Available at: https://ieeexplore.ieee.org/document/996017Links ]

[19] Zitzler, E., Laumanns, M., and Thiele, L., SPEA2: improving the strength pareto evolutionary algorithm. ETH Zentrum. Research Collection, 103, pp. 1-22, 2001. DOI: https://doi.org/10.3929/ethz-a-004284029Links ]

[20] Huang, W., Zhao. Y., and Ma, X., An improved genetic algorithm for job-shop scheduling problem with process sequence flexibility. International. Journal of Simulation Modeling, 13(4), pp. 510-522, 2014. Available at: http://m.ijsimm.com/Full_Papers/Fulltext2014/text13-4_510-522.pdfLinks ]

[21] Wang, Y., Stein, N., Bäck, T., and Emmerich, M., Improving NSGA-III for flexible job shop scheduling using automatic configuration, smart initialization and local search. GECCO '20: proceedings of the 2020 genetic and evolutionary computation conference companion, [online]. 2020. pp. 181-182. Available at: https://scholarlypublications.universiteitleiden.nl/handle/1887/3618523Links ]

[22] Xie, J., Li, X., Gao, L. and Gui, L., A hybrid algorithm with a new neighborhood structure for job shop scheduling problems. Computers and Industrial Engineering, 169, art. 108205, 2022. DOI: https://doi.org/10.1016/j.cie.2022.108205Links ]

[23] Aschauer, A., Roetzer, F., Steinboeck, A., and Kugi, A., Efficient scheduling of a stochastic no-wait job shop with controllable processing times. Expert Systems with Applications, 162, art. 113879, 2020. DOI: https://doi.org/10.1016/j.eswa.2020.113879Links ]

[24] Geng, Z., Yuan, J., and Yuan, J., Scheduling with or without precedence relations on serial-batch machine to minimize makespan and maximum cost. Applied Mathematics and Computation, 332, pp. 1-18, 2018. DOI: https://doi.org/10.1016/j.amc.2018.03.001Links ]

[25] Safarzadeh, H., and Kianfar, F., Job shop scheduling with the option of jobs outsourcing. International Journal of Production Research, 57(10), pp. 3255-3272, 2019. DOI: https://doi.org/10.1080/00207543.2019.1579934Links ]

[26] Zheng, F., and Sui, Y., Bi-objective optimization of multiple-route job shop scheduling with route cost. IFAC-PapersOnLine, 52(13), pp. 881-886, 2019. DOI: https://doi.org/10.1016/j.ifacol.2019.11.241Links ]

[27] Golpîra, H., Rehman, K., and Zhang, Y., Robust Smart energy efficient production planning for a general job-shop manufacturing system under combined demand and supply uncertainty in the presence of grid-connected microgrid. Journal of Cleaner Production, 202, pp. 649-665, 2018. DOI: https://doi.org/10.1016/j.jclepro.2018.08.151Links ]

[28] Huo, D., Xiao, X., and Pan, Y., Multi-objective energy-saving job-shop scheduling based on improved NSGA-II. International Journal of Simulation Modelling, [online]. 19(3), pp. 494-504, 2020. Available at: http://www.ijsimm.com/Full_Papers/Fulltext2020/text19-3_CO12.pdfLinks ]

[29] Masmoudi, O., Delorme, X., and Gianessi, P., Job-shop scheduling problem with energy consideration. International Journal of Production Economics, 216, pp. 12-22, 2019. DOI: https://doi.org/10.1016/j.ijpe.2019.03.021Links ]

[30] Ren, J., Ye, C., and Li, Y., A two-stage optimization algorithm for multi-objective job-shop scheduling problem considering job transport. Journal Europeen des Systemes Automatises, 53(6), pp. 915-924, 2020. DOI: https://doi.org/10.1007/s10845-009-0294-6Links ]

[31] Renna, P. and Materi, S., Switch off policies in job-shop manufacturing systems including workload evaluation. International Journal of Management Science and Engineering Management, 16(4), pp. 254-263, 2021. DOI: https://doi.org/10.1080/17509653.2021.1941369Links ]

[32] Salido, M., Escamilla, J., Barber, F., and Giret, A. rescheduling in job-shop problems for sustainable manufacturing systems. Journal of Cleaner Production, 162(20), pp. 121-132, 2017. DOI: https://doi.org/10.1016/j.jclepro.2016.11.002. [ Links ]

[33] Wen, X., Wang, K., Li, H., Sun, H., Wang, H., and Jin, L., Two-stage solution method based on NSGA-II for Green Multi-Objective integrated process planning and scheduling in a battery packaging machinery workshop. Swarm and Evolutionary Computation, 61, art. 100820, 2021. DOI: https://doi.org/10.1016/j.swevo.2020.100820Links ]

[34] Zhao, J., Peng, S., and Li, T., Energy-aware fuzzy job-shop scheduling for engine remanufacturing at the multi -machine level. Front. Mech. Eng. 14, pp. 474-488, 2019. DOI: https://doi/org.ezproxy.unal.edu.co/10.1007/s11465-019-0560-zLinks ]

[35] Duan, J., and Wang, J., Robust scheduling for flexible machining job shop subject to machine breakdowns and new job arrivals considering system reusability and task recurrence. Systems with Applications, 203, art. 117489, 2022. DOI: https://doi.org/10.1016/j.eswa.2022.117489Links ]

[36] Para, J., Del Ser, J., and Nebro, A., Energy-aware multi-objective job shop scheduling optimization with metaheuristics in manufacturing industries: a critical survey, results, and perspectives. Applied Sciences. 12(3), art. 1491, 2022. DOI: https://doi.org/10.3390/app12031491Links ]

[37] Daneshamooz, F., Fattahi, P., and Hosseini, S., Mathematical modeling and two efficient branch and bound algorithms for job shop scheduling problem followed by an assembly stage. Kybernetes. ahead-of-print. 2021. DOI: https://doi.org/10.1108/K-08-2020-0521Links ]

[38] Shen, X., Zhang, M., and Fu, J., Multi-objective dynamic job shop scheduling a survey and prospects. International Journal of Innovative, [online]. 10(6), pp. 2113-2126, 2014, Available at: http://www.ijicic.org/ijicic-14-01022.pdfLinks ]

[39] Mu, H., Disruption management of flexible job shop scheduling considering behavior perception and machine fault based on improved NSGA-II algorithm. Journal Européen des Systèmes Automatisés, 52(2), pp. 149-156, 2019, DOI: https://doi.org/10.18280/jesa.520206Links ]

[40] Santos, V., Carvalho, F., Assis, L., Weiss-Cohen, M., and Guimarães, F., Multi-objective iterated local search based on decomposition for job scheduling problems with machine deterioration effect. Engineering Applications of Artificial Intelligence, 112, art. 104826, 2022. DOI: https://doi.org/10.1016/j.engappai.2022.104826Links ]

[41] Shokouhi, E., Integrated multi-objective process planning and flexible job shop scheduling considering precedence constraints. Production & Manufacturing Research, 6(1), pp. 61-89, 2018. DOI: https://doi.org/10.1080/21693277.2017.1415173Links ]

[42] Zhang, S., Li, X., Zhang, B., and Wang, S., Multi-objective optimization in flexible assembly job shop scheduling using a distributed ant colony system. European Journal of Operational Research, 283(2), pp. 441-460, 2020. DOI: https://doi.org/10.1016/j.ejor.2019.11.016Links ]

[43] Zhu, H., He, B., and Li, H., Modified bat algorithm for the multi-objective flexible job shop scheduling problem. International Lournal of Performability Endineering, [online]. 13(7), pp. 999-1012, 2017. Available at: http://www.ijpe-online.com/EN/10.23940/ijpe.17.07.p1.9991012Links ]

[44] Gayathri, D., Mishra, R., and Madan, A., A dynamic adaptive firefly algorithm for flexible job shop scheduling. Intelligent Automation & Soft Computing. 31(1), pp. 429-448, 2022. DOI: https://doi.org/10.32604/iasc.2022.019330Links ]

[45] Renke, L., Piplani, R., and Toro, C., A review of dynamic scheduling: context, techniques and prospects. In: Toro, C., Wang, W., and Akhtar, H., Eds. Implementing Industry 4.0. Intelligent Systems Reference Library. 202. Springer, Cham, [online]. 2021. DOI: https://doi/org.ezproxy.unal.edu.co/10.1007/978-3-030-67270-6_9Links ]

[46] Dabbagh, R., and Yousefi, S., Hybrid decision-making approach based on FCM and MOORA for occupational health and safety risk analysis. Journal of Safety Research, 71, pp. 111-123, 2019. DOI: https://doi.org/10.1016/j.jsr.2019.09.021Links ]

[47] Villicaña, E., and Ponce, J., Sustainable strategic planning for a national natural gas energy system accounting for unconventional sources. Energy Conversion and Management, [online]. 181, pp. 382-397, 2019. Available at: https://www.scopus.com/inward/record.uri?eid=2-s2.0-Links ]

[48] García, F., Medina, S., Gonzales, R., Huertas, I., Ferrari, A., and Settembre, D., Industry 4.0-based dynamic Social Organizational Life Cycle Assessment to target the social circular economy in manufacturing. Journal of Cleaner Production. 327, art. 129439, 2021. DOI: https://doi.org/10.1016/j.jclepro.2021.129439Links ]

[49] Autuori, J., Hnaien, F., and Yalaoui, F., Three metaheuristics improved by a mapping method. IFAC - International Federation of Automatic Control, 49(12), pp. 1472-1477, 2016. DOI: https://doi.org/10.1016/j.ifacol.2016.07.779Links ]

[50] Zhang, R., and Chiong, R., Solving the energy-efficient job shop scheduling problem: a multi-objective genetic algorithm with enhanced local search for minimizing the total weighted tardiness and total energy consumption. Journal of Cleaner Production, 112(4), pp. 3361-3375, 2016. DOI: https://doi.org/10.1016/j.jclepro.2015.09.097Links ]

[51] Rashno, A., Shafipour, M., and Fadaei S., Particle ranking: an efficient method for multi-objective particle swarm optimization feature selection. Knowledge-Based Systems, 245, art. 108640, 2022. DOI: https://doi.org/10.1016/j.knosys.2022.108640Links ]

[52] Zhang, F. , Bail, J. , Yang, D., andWang, O. , Digital twin data-driven proactive job-shop scheduling strategy towards asymmetric manufacturing execution decision. Scientific Reports, [online]. 12, art. 1546, 2022, Available at: https://www.nature.com/articles/s41598-022-05304-wLinks ]

How to cite: Coca-Ortegón, G., Occupational hazards and economic indicators in the scheduling of a make-to-order system. DYNA, 90(227), pp. 117-125, July - September, 2023.

G. Coca-Ortegón Is a professor in the Universidad EIA, Colombia. He is a PhD in Engineering, MSc. in Engineering - Industrial Engineering from the Universidad Nacional de Colombia. He is consultant in quality, productivity and sustainability. He worked for several manufacturing and financial companies. ORCID: 0000-0002-6760-4111

Received: February 02, 2023; Revised: June 29, 2023; Accepted: July 21, 2023

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License