Selecting the Optimal LoA to Prevent the Expansion of COVID-19 in the Chemical Industry considering Sustainability Factors: A Fuzzy Mathematical Optimization Approach

Automation has attracted interest from the industry sector for its potential to improve energy eciency, cost eciency, and environmental performance. By elevating the LoA to the highest degree, associated costs will grow accordingly and its implementation will be far more complicated. is will also result in losing workers and decreasing environmental pollutants. On the other hand, increasing power consumption at high levels of automation leads to the production of greenhouse gases. is paper aims to increase the level of automation (LoA) considering the concept of sustainability. is study presents fuzzy multiobjective programming to determine the optimal LoA considering sustainability factors to achieve competitive advantages. To solve the model, the Zimmermann max-min approach was adopted and a cosmetics factory in Iran was chosen to optimize LoA according to this model. e results showed that it is possible to improve the LoA and also consider sustainability factors with the available resources without using the highest LoA. is study can help managers optimize the LoA in their organizations considering the current resources and sustainability issues, and control the company’s return on investment and cost of overhead. ey can run the model with every de nition of LoA proposed till now. is research can bene t the environment and the workers’ health in the production line by reducing environmental pollutants and prevent the dismissal of all personnel due to its negative social eects. It also reduces the risk of COVID-19 by minimizing the number of workers. So far, a mathematical model for selecting optimal LoA in the chemical industry considering sustainability has not been presented.


Introduction
e high competition among producers and a variety of products has caused organizations to optimize LoA which entails bene ts like an increase in quality as well as the production speed, more accurate and faster quality control (QC), a reduction in production waste, better interaction with business systems, an increase in the productivity of industrial units, an increase in the safety factor for manpower and the reduction of mental and physical stress among workers [1]. e implementation of high LoA is a very important issue due to the decrease in dangerous gases emitted from chemical interactions and the transportation hazards of chemicals that threaten the workers in the chemical industry [2]. However, implementing automation at high levels entails heavy expenditure which most organizations cannot a ord. Moreover, performing the jobs manually prolongs the production process that in turn increases the chance of in ltration of microbes into the product and can damage the environment and the consumers; therefore, chemical industries can reduce these hazards and the production time by increasing the LoA [3]. Although, at high levels of automation, the number of workers reduces which causes social concerns in the chemical industry. Due to the prevalence of COVID-19 as a pandemic, customer health has become more important than before, and to prevent the risk of this disease, manufacturers are forced to use high LOA in the process of production to prevent endangering the health of customers. It should be mentioned that COVID-19 spreads not only from person to person via close contact respiratory droplet transfer but also on contaminated surfaces. Coronavirus can persist on inanimate surfaces including metal, glass, or plastic for days. Another important reason to increase LoA is to decrease the danger of transferring COVID-19 and avoid transferring viruses to workers and customers. Also implementing optimal LoA can help organizations attain the goals of sustainability. Sustainability is a multifaceted concept. It forces organizations to set goals in three fully interrelated aspects as below: (i) Providing social responsibility in which the needs of all stakeholders are met; (ii) Effective protection of the environment and accurate use of natural resources; (iii) Providing economic growth and economic prosperity.
Few studies have focused on a specific aspect of the problem, such as environmental, social, and economic factors. is separates the model from the real world because these issues need to be considered together. erefore, the main objective of this research is to determine the optimal LoA for the chemical industry considering three dimensions of sustainability. So the main questions are: What is the optimal LoA in the chemical industry considering the different aspects of sustainability in the chemical industry? Which criteria and parameters should be considered for environmental and social objectives? To answer these questions we developed a bi-objective mathematical model that suggests the optimal LoA considering the concept of sustainability. Also, it keeps some workers in workstations and uses the current budget considering net present value (NAV). e model also reduces pollution in the factory.
(i) e main contribution of this work lies in developing a bi-objective mathematical model that suggests the optimal LoA considering current resources and the concept of sustainability. e remainder of the paper is as follows. e related literature is reviewed in Section 2 and in Section 3, the problem is discussed. e mathematical model is defined and the objective functions, constraints, and solution approach are described. In Section 4, the implementation of the model is presented, respectively, in an Iranian cosmetics factory. Section 5 provides managerial insights and practical implications. In conclusion, the findings and recommendations for further research are presented in Section 6.

Automation and Autonomy.
e term automation was coined in the early 1940s to determine the different mechanisms for the assignment of tasks requiring human monitoring, control, and intervention to automated machines and systems. e degree to which tasks are performed automatically in an industrial or service unit is called levels of automation [4]. e first definition of LoA presented by Sheridan and Verplank [5] defined ten levels of automation for industries. It is based on six activities that humans and systems do in the production process that involves getting, selecting, starting, requesting, approving, and telling. Many other definitions were given till now [6][7][8][9][10][11][12][13][14][15][16]. After that, Riley [17] offered a different definition. He considered a 2-D matrix, in which rows show the levels of automation while the columns correspond to the dimensions of automation. Vagia et al. [18] reviewed all research papers about LoA and developed a new taxonomy. In contrast to a conventional automated system designed to carry out a limited set of preprogrammed supervised tasks on behalf of the user, autonomy is the technology (either hardware or software) designed to carry out a user's goals and does not require supervision [1]. In this sense, successful autonomy is considered to be well designed and highly capable of automation and can adapt to a wider variety of conditions better [19].
is concept is in agreement with Hancock's idea [20] who described autonomy as a later evolution of automation which was historically more restricted in capability and scope than autonomy. Most authors do not differentiate between automation and autonomy but Fereidunian [21] and Parasuraman [14] made a difference between the levels of autonomy and automation. Determining the LoA has a wide range of applications in different industries for example in avionics, teleportation systems, remote control operations, and aircraft control.
In recent research on LOA, Mostafa et al. [22] reviewed 171 research papers and provided a fundamental understanding of adjustable autonomy and its application. Also, Malek [23] considered various definitions of the term automation. His paper aimed at the study, analysis, and discussion of several definitions of the LoA proposed in the literature. On the grounds of this analysis, a set of requirements to be gained by an accurate indicator was given out to propound a new definition of the LoA to be applied in the manufacturing domain to optimize workstations. At the end of this bibliographical study, a table summarizing the list of different definitions and the requirements that each of them meets is drawn. A proposal for a new definition of LoA as a time ratio is also presented. By studying the literature review, we developed a fuzzy multi-objective model that can improve LoA in the chemical industry by considering any one dimensions (1D) or two dimensions (2D) definition of LoA that has been defined before.

Appropriate
LoA. Despite the large benefits of automation, it should be mentioned that the highest LoA is not necessarily the best option and the most appropriate level for many industries; therefore, the optimal LoA for each industry should be determined by considering some issues [5]. e main research for optimizing LoA was done by Frohm [5] that proposed Dynamo's methodology to determine the levels of automation in organizations.
is methodology consists of the following stages: planning steps, a preliminary study for identifying activities, documenting the flow of products and activities, identifying the main activity of each production workstation, identifying the sub-activities of each production workstation, measuring the LoA, evaluating the LoA, and analyzing the results. Fasth et al. [9] followed Frohm's research [5] by presenting the Dynamo++ methodology in 12 steps that can measure the current LoA and provide suggestions for increasing the LoA and eventually help the executives to select the optimal LoA. e Dynamo and Dynamo ++ methodologies are time-consuming in terms of implementation and need to simulate the optimal scenario for their implementation. Many types of research are also done using Dynamo and Dynamo++ methodologies: Fasth and Stahre [9] examined six industrial groups that needed to change their automation level and used Dynamo methodology to evaluate the current LoA. In these industrial groups, the main parameter for changing the LoA was the flexibility or production time. ey concluded that it is not always necessary to change the LoA. Stahre and Fasth [24] calculated and analyzed the levels of automation in the assembly industry using Dynamo++ in 12 stages so that the best LoA could be determined considering efficiency and cost. e impact of automation on human resources, waste reduction, efficiency, and cost in different industries was investigated in a study by Wang [25], too.
Choe et al. [26] used the Dynamo++ methodology to calculate the LoA in a truck manufacturing factory and simulated the impact of automation on flexibility in the material transportation system. Before this study, the Dynamo methodology was used in assembling systems. Mehta and Subramanian [27] investigated and explored the barriers that a company would face while increasing the LoA in the preassembly production unit. To achieve the primary goal of investigating the barriers, their study took a threefold approach. ey measured first the current LoA for the preassembly workstations. is measurement was conducted by incorporating an existing methodology adapted from the Dynamo++ methodology. ey concluded that this methodology could be incorporated in measuring and analyzing the current LoA of the preassembly workstations. Hadi and Brillinger [28] declared that achieving high quality and high variety batch size production could be quite expensive. e focus of their research lies on high adaptive and cognitive aspects in the assembly along with qualitative aspects. ey presented a level of practical application matrix of all the possible adaptive technologies that were feasible to implement in the preassembly line using Frohm's [5] Dynamo methodology.
By studying the literature review, it can be said that there is no need to use the highest LoA in an organization, and according to the current conditions of organizations in terms of human resources, equipment, and financial resources, the optimal LoA can be selected to suit them. Most research studies show that the used Dynamo methodology is too complex and time-consuming. For this purpose, we have selected the optimal LoA in the desired factory for our case study using mathematical modeling that is user-friendly and can be run fast and considers sustainability. is model has not been used till now.
Recently, some researchers have been done in LoA such as: Krishnamoorthi [29] demonstrated a methodology for evaluating multiple construction processes and selecting an optimal solution. e methodology combined compositional modeling, case-based reasoning, and stochastic search. Potential solutions consisting of automated construction processes were explored by generating various combinations of process fragments. Malek [30] aimed to propose a new methodology for organizing and identifying the assembly operations that ought to be automated in an automotive assembly line. e different requirements of the methodology were defined, which led to a proposal for a method that respects all the requirements and allows not only the grouping of operations but also the analysis of the automation and the line balancing. Imset Marius [31] studied a drilling unit (MODU) from a subsea oil and gas well. He applied a framework for levels of automation to explore the critical decision process leading to an EQD. He also provided an overview of the benefits and drawbacks of existing automation and decision support systems vs. manual human decision-making. is paper summarizes the growth of Industry 4.0 during the last 5 years and provides a concise background overview of Industry 4.0-related works and its various application areas. Shastri [32] developed a model for Levels of automation for IT services based on Bloom's taxonomy taken as the reference and assessed the automation scope for all the processes of ITIL. Aryal [33] studied heating, ventilation, and air conditioning (HVAC) systems. In this article, he described the development and implementation of an Internet-of-ings (IoT)-based intelligent agent that learns individual occupant comfort requirements and controls the thermal environment using PCS (i.e., a local fan and a heater). e results showed that PCS used improved occupant satisfaction and including some LoA can improve occupant satisfaction further than what is possible with manually operated PCS. Among the levels of automation investigated, inquisitive automation, where the user approves/declines the control actions of the intelligent agent before execution, led to the highest occupant satisfaction with the thermal environment. References [34][35][36][37] studied LoA and social impacts on Self-Driving Vehicles.

Selecting Optimal LoA considering Sustainability Goals: A
Research Gap. Usually, Dynamo and Dynamo ++ methodologies have been adopted to select the optimal LoA. e implementation of these methodologies is extremely complex and time-consuming. To the best of our knowledge, no research has ever been conducted to determine the optimal LoA by using a mathematical modeling approach. Also, sustainable development factors (cost, social, and environmental factors ) have not yet been considered in selecting the optimal level and dimension of automation. In this study, a mathematical model will be presented to determine the optimal LoA based on both one dimension (1D) and two Discrete Dynamics in Nature and Society dimensions (2D) definitions of LoA considering the three dimensions of sustainability. In brief, the research gaps are as follows: Many studies have focused on a specific aspect of the problem, such as environmental, social, and economic aspects. is separates the model from the real world because these issues need to be considered together; A large number of studies have used examples generated with random numbers to validate their models, but it is better to validate models using realworld examples and to write models for real case studies; (ii) ere was no research conducted to determine the optimal LoA by using a fuzzy multi-objective model approach.
e goal of this research is to present a mathematical model for selecting the optimal LoA in the chemical industry based on sustainability considerations to reduce the danger of transmitting viruses like COVID-19. e research has been done in the chemical industry in Iran. In Table 1, a review of the LoA has been done and the research gap is shown.
In this study the researchers studied the previous research in LoA optimization, optimized LoA considering sustainable factors with a Fuzzy mathematical optimization approach, defuzzilized the model using LHS, used the Zimmermann max-min approach, conducted a sensitivity analysis (SA) to evaluate the importance of model inputs, and then examined the model in a case study. e differences between this study and the previous studies are as follows.
is study optimizes LoA with modeling, considers sustainable factors, and uses a case study.

Problem Statement
During chemical production, if the process of production is long and manually performed, which is usually done at low levels of automation, there is a high probability of infiltration of microbes and viruses like COVID-19 into the product and causing damage to workers and consumers. Consequently, to lessen the level of expected harm chemical manufacturing industries have to reduce the production time by increasing the LoA. It should be mentioned that using full automation has a high cost and most organizations cannot afford it [38]. Full automation  also leads to the downsizing of workers and has negative social effects, consumes more electricity, and produces greenhouse gases. Environmental pollution, in its place, is one of the main problems in society and governmental agencies impose heavy rules and fines for it. ere is even a risk of closure of the organization which cannot manage it as well. Taking these issues into account, it is essential to select the optimal LoA in chemical industries by considering the existing financial resources of organizations and maintaining a certain number of workers in production lines which reduces the social damage caused by downsizing the workforce and decreases the environmental pollutants. According to the mentioned issues, the major concern in this research is selecting the optimal LoA in the chemical industry considering sustainability as follows: Reducing production costs in different production workstations, including the costs of raw materials, labor costs, and overhead costs, maintaining a certain number of production line workers, and reducing environmental pollutants.
A picture of the problem statement is shown in Figure 1.
In this research, first, the literature review is done. e purpose of this section is to determine the research gap. In the second step, a mathematical model is designed. In this model, which was run for a long time, the Net Present Value (NPV) or the current value of cash flows with the desired return rate on the project will be positive relative to the initial capital. Also, the selection of the level and dimension of automation will be measurable with one of the taxonomies. e model applies to all taxonomies of LoA including one-dimensional taxonomy (1D) [39] or two-dimensional (2D) [17]. If the 2D taxonomy is used, the LoA in different dimensions is homogenous; for example, in one dimension it is not fully automatic and in the other dimension it is manual. In this study, the method developed by McKay [40]  Discrete Dynamics in Nature and Society  Discrete Dynamics in Nature and Society is employed to defuzzilize the model due to its effectiveness and computational efficiency. en, a case is stated and in the last step, to solve the mathematical model, the Zimmermann max-min approach is used and the results are analyzed. Finally, a sensitivity analysis (SA) is conducted to evaluate the importance of model inputs. e steps are shown in Figure 2.

Assumption of the Model
e amount of material consumption in each workstation is specified and predefined, and it is definite in each dimension and for every level of automation in a specific period; e number of specialists in each workstation is specific, definite, and predefined in each dimension and every LoA; e cost of materials is uncertain due to changes in exchange rates and world prices; e cost of specialists and workers in each workstation is certain in each dimension and for every level of automation in a specific period; (v) Overhead costs are uncertain because they vary in different geographical areas and different seasons; (vi) e chemical production process includes three workstations-manufacturing, filling, and packaging; (vii) e model runs for a long time.

Notations.
e sets, parameters, and variables used in the mathematical model are given as follows (see Table 2):

Objective Functions.
e first objective function (1) minimizes all costs associated with the implementation of the automation. Four different costs consist of fixed costs(line one), raw material costs(line two), overhead costs(line three), and labor costs(line four) formulated as follows: Line one represents the total fixed costs (purchasing machinery) at different levels and different dimensions in manufacturing workstations, filling workstations, and packaging workstations. Line two shows the total cost of raw material at different levels and different dimensions in manufacturing workstations, filling workstations, and packaging workstations.
Line three shows the total overhead costs that are equal to the multiplication of the overhead cost, activity cost driver, and the volume of product at different levels and different dimensions in manufacturing workstations, filling workstations, and packaging. Line four represents the total labor cost at different levels and different dimensions in manufacturing workstations, filling workstations, and packaging workstations. e equation is as follows: In the second objective function (2), the amounts of pollution in the manufacturing workstation, filling workstation, and packaging workstation during the period t are minimized, as can be seen in this equation:

Constraints
Discrete Dynamics in Nature and Society jm j�1 pv P ≥ DEMAND, Constraints (3)-(5) ensure that only one LoA can be selected in each dimension in the manufacturing workstation, filling workstation, and packaging workstation. Constrain (7) is introduced to consider social damage and states that the maximum fully automated level (without the operator) in the manufacturing workstation is ml LoA. Constrain (8) is introduced to consider social damage and states that the maximum fully automated level (without the operator) in the manufacturing workstation is fl LoA. Constrain (9) is    Discrete Dynamics in Nature and Society introduced to consider social damage and states that the maximum fully automated level (without the operator) in the manufacturing workstation is pl LoA. Constraints (9)- (11) monitor that fixed costs (purchase of machinery) do not exceed the intended budget of the organization (MAXCM), the consumption of raw materials does not exceed the intended budget of the organization (MAXCR), and overhead costs do not exceed the budget of the organization (MAXCO), respectively. Constraint (13) states that market demand should be met. Constraint (14) is the relationship between the number of products in the filling as well as the packaging workstation (i.e., the number of bottles in each carton packaging). Constraint (14) is the relationship between the number of products in the manufacturing and filling workstations (i.e., the amount of bulk used in each bottle that is filled in the filling workstation). Constraint (15) controls the company's return on investment during tm years. Constraints (16)- (18) are to keep the levels of automation in each workstation homogenous in different dimensions. Constraints (19)- (21) are to keep the levels of automation in all workstations homogenous. e difference in the LoA in the three workstations of manufacturing, filling, and packaging should not be great. For example, if the LoA in a manufacturing workstation is high but the LoA in the filling workstation is low, the filling workstation will not respond properly. Constraint (22) determines that the manual levels in workstations are not allowed. Constraint (23) determines the type of decision variables of the model.

Zimmermann Max-Min Approach.
To solve the problem, not all goals may be achieved simultaneously under system constraints. In this situation, some tolerance is defined in the model and for this reason, Zimmermann maxmin approach [41]is used to solve the bi-objective problem.
is method has the following steps: (1) Start (2) Solve first and second objective functions separately (3) Create a pay-off table for the first and second objective functions separately as follows Table 3: (4) Define the membership function for the first and second objective functions separately Table 4: (5) Figure 3 shows the linear membership function for minimizing the goal.   ). en, the membership functions are added to constraints for each objective function as follows: It should be noted that with the introduction of the Zimmermann maxmin and Lp-Metric methods, the introduced multiobjective model becomes a single-objective model. en, the model was run using the GAMS software on a system with a Corei7 processor.

Lp-Metric Method.
In this research, the Lp-metric method has been used to integrate objective functions. In this method, the deviation of the objective functions from their optimal value is considered. At first, the individual answers are calculated for the optimization of each objective function and then the objective function is minimized as follows: where w k is the degree of importance for the i-th objective function and 1 ≤ p < ∞ is the parameter defining the Lpmetric method. e value of p determines the degree of emphasis on the existing deviations so that the larger the value of p, the greater the emphasis on the largest deviation. Usually, the values of p � 1, p � 2, and p � ∞ are used in calculations; p � 1 indicates that the same importance is considered for all deviations and p � 2 indicates that each deviation has its weight so that the largest deviation takes the most weight. When p tends to infinity, the largest deviation indicates distance.
is method is converted to the Min-Max approach for the value p � ∞. e variable λ is defined as follows: erefore, the multi-objective model can be written as a single-objective model, as to the following relations:

Results
e case is a cosmetic factory that produces different types of hygienic-chemical products in Iran. e products manufactured at the factory include lipstick, eyeliner, mascara, creams, shampoos, varnishes, powder creams, and other products. Some workstations use a low LoA and employ many workers. However, in the Covid-19 epidemic, the managers have to reduce the number of workers to prevent transmission of the Coronavirus to other workers and customers. Also, due to the costly nature of chemical production and the competitive business environment in this industry concerning the quality level, increasing LoA is very important for reducing the final product's manufacturing costs and increasing the amount and quality of production. But using full automation has a very high cost and its implementation is too complex and time-consuming so the factory cannot perform it. Also, they have to fire many workers and pay a lot of penalties which can result in very negative social effects. On the other hand, it is important to reduce environmental pollutants due to the controls practiced by governmental agencies and the social responsibilities thereof in a way that if the rules are not observed, they will have to pay heavy penalties. erefore, the aim of implementing the proposed model in this factory is to determine the optimal LoA so that costs and environmental pollution are minimized and several workers are retained with a specific budget. Moreover, in the long run, the Net Present Value (NPV) or the current value of cash flows with the desired return rate on the project will be positive relative to the initial capital. To better understand the content, we examine the production processes as follows: e production line includes three workstations of production-manufacturing, filling, and packaging. e process flow of cosmetics manufacturing under study is shown in Figure 4.
Proud's definition [15] was used for the implementation of the model. is definition has four dimensions (observe, orient, decide, and act) and eight LoAs.
e LoA in the manufacturing workstation and filling workstation in dimensions observe, decide, and act is 3 and in the orient dimension is 2.
e LoA in the packaging workstation in the orient dimension is 2 and in the other dimensions is 3. e current LoA in this factory is shown in Figure 5. e raw material is added to the reactor in the manufacturing workstation. Unit costs of raw materials in manufacturing workstation, filling workstation, and packaging workstation are fuzzy variables. Also, overhead costs in workstations are fuzzy variables. We considered four overhead costs including the cost of electricity, the cost of water, depreciation cost, and repair cost as D1, D2, D3, and D4.
To deal with uncertainty, the LHS method was used to generate the scenario. Data from the last two years were used and the number of samples was 25 according to Morgan's table. Using the data, a cumulative distribution diagram was drawn and random scenarios were generated using coding in MATLAB as shown in Figures 6-8.
We used the data of three types of specialists who participated in the production processes at different automation levels including simple workers, excellent workers, and specialists. In manufacturing cosmetics, the mixing reactor is regarded as the heart of the process. In this process, the purified and deionized water is transferred into the mixer tank first. en, the materials are pumped or manually fed into the reactor to be mixed, being simultaneously heated. e amount and duration of heating depend on the process. e next phase is storing the materials in a large container followed by packaging which includes filling, capping, and labeling. e filling machine is a device with at least one tank and a nozzle for filling different bottles. After that, the final phases are packaging and putting the bottles in cartons, which can be performed manually or automatically. It should be noted that in low levels of automation, the workers are present in all production processes, especially packaging and for this reason the risk of transmission of viruses is high. During the manufacturing process, chemical reactions produce several gases including primary pollutants. e pollutants in the cosmetics industry include particulate matter (PM), oxides of carbon, sulfur, nitrogen, ozone, free radicals, and other airborne chemicals like pesticides, chemical sprays, and hydrocarbons. A secondary pollutant is not directly emitted as such but is formed when other  Table 6: Z and ƛ values were gathered using Zimmermann max-min approach in Scenario 1. pollutants (primary pollutants) react. Examples of secondary pollutants include ozone which is formed when hydrocarbons (HC) and nitrogen oxides (NOx) combine in the presence of sunlight, NO2 which is formed as NO combined with oxygen in the air, and acid rain which is formed when sulfur dioxide or nitrogen oxides react with rain. Another type of pollutant is the sewage made as the result of production processes, the most serious of which include surfactants, detergents, and phosphorus. Major pollutants comprise surfactants, detergents, and phosphorus which are required to be recycled according to environmental standards. In this study, we pay attention to CO2 emissions and greenhouse gases in the production process.

4.1.
Discussion. e developed model was coded in GAMS 24.9.2/CPLEX 12.7.1.0 as solver and all cases were run on a computer with a 3.00 GHz processor Intel ® core ™ i7 processor and 64 GB memory RAM 64 bit. For running the model, three scenarios were selected. e results of implementing various scenarios are presented here: Scenario 1. e planning horizon is 5 years (tm � 5). e number of workers in the manufacturing workstation is 1 (nm � 1), the number of workers in the filling workstation is 4 (nf � 4), and the number of workers in the packaging workstation is 4 (np � 4). e cost of buying machinery is $ 100,000 (b � 4; p � 6; d � 3). Scenario 2. e planning horizon is 4 years (tm � 4). e number of workers in the manufacturing workstation is 1 (nm � 1), the number of workers in the filling workstation is 1 (nf � 1), and the number of workers in the packaging workstation is 5 (np � 4). e cost of buying machinery is $ 200,000 (b � 4; p � 6; d � 3). Scenario 3. e planning horizon is 3 years (tm � 3). e number of workers in the manufacturing workstation is 0 (nm � 0), the number of workers in the filling workstation is 1 (nf � 1), and the number of workers in the packaging workstation is 1 (np � 1). e cost of buying machinery is $ 400,000 (b � 4; p � 6; d � 3).
It can be concluded from the scenario1 that if the organization is not able to pay the costs of purchasing machinery with a high LoA, it can compensate for it by increasing the production level over five years with a low LoA, but the number of production workers will still be really large. In scenario 2, with an increase in investment, the organization will be able to compensate for its costs for four years and some workers are maintained in packaging and filling workstations. On the contrary, in the manufacturing workstation which causes the greatest damage to workers and the environment, there is only one worker in the production lines and the number of pollutants in the environment is relatively low. By reducing the number of workers, we can reduce the danger of the microbial conditional of the product, which results in less damage to the environment. In scenario 3, with an increase in investment compared with scenario 2, the organization will be able to compensate for its costs over three years, and a few workers are maintained in packaging and filling workstations. On the contrary, in the manufacturing workstation that causes the greatest damage to workers and the environment, the completely automated level is used. At the full automated level, there is no worker in production lines and the amount of pollutants in the environment is almost zero. Also, the microbial potential of the product is almost zero and hence there will be no damage to the environment. It is suggested that scenario 3 be used in the factory. is does low damage to the environment and some production workers remain in production lines. It is proportionate to the rich organization's budget. If it is not possible, an increase in the investment amount of scenario 2 is proposed due to the short term of its return on investment. As it is seen, higher fixed costs bring about higher levels of automation with higher profitability and the period of return on investment will be shorter. But, the number of workers is reduced at a high LoA. To prevent social damage caused by downsizing the workers, medium levels of automation are proposed in the packaging workstation because there is no environmental pollution in these lines at lower levels of automation. However, we usually use high levels of automation in the manufacturing workstations because they have beneficial environmental impacts. It must be considered that LoA in different workstations must be homogenous. In running these scenarios, the difference between LoAs in workstations should be less or equal to 4. e results of implementing scenarios 1, 2, and 3 with the Zimmermann max-min approach in three workstations are given in Tables 5 and 6. In the third scenario, the implementation cost is high, but the amount of environmental pollutants is very low.
Results show that the best scenario is scenario 3. Figure 9 shows that in scenario 3 in the manufacturing workstation in all dimensions the LoA is eight and in the filling workstation and packaging workstation in four dimensions it is seven. In scenarios 1 and 2, the LoA is lower than in scenario 3 and the manufacturing workstation tries to use upper LoA because    Table 7 presents the results of the running of the proposed scenarios and the initial state of production. As can be seen, in the costly but profitable scenario, the production time, the amount of waste, and the number of workers in the production line are lower but the production quantity and profitability are significantly higher.
Other results of the implementation of the proposed scenarios are given in Table 8. Table 9 shows that in the lowcost scenario, the production time and the amount of product waste are slightly reduced. In the expert scenario, the amount of waste and the time of production are further reduced. In the costly but profitable scenario, the production time and the amount of waste are significantly reduced. In Figures 12(a,b,c, and d), the variations of production time, the amount of waste, production quantity, and profit in the initial state of production before the change (s0), and the three proposed scenarios (s1, s2, s3) are compared; it was found that the s3 scenario, i.e., costly but profitable, has the least waste, the least productive time, and the maximum profit as shown in Figure 13.
After increasing the number of experts, the results show that the amount of Z1 increases and the LoA decreases. e number of three types of experts is shown in three workstations in Table 9 and Figure 10; the amounts increased and the Z1 value increased, too.

Sensitivity Analysis.
e results show that the amount of Z1 increases after increasing the overhead costs. e results     Table 10. e results of sensitivity analysis show that increasing the cost of raw materials, the number of workers, and overhead costs leads to an increase in the value of Z1, but it is almost ineffective on the value of Z2 as shown in Table 11 and Figures 12, 14 and 15.
(i) Discussion: e Dynamo method that was previously used to determine the optimal LoA was very complicated and time-consuming, and simulation should be used to implement it, but the model in this research is very practical, simple, and low-cost and it also covers all aspects of sustainability.

Managerial Insights and Practical Implications
is model can help executives to select the appropriate level and dimension of automation proportionate to organization resources and requirements. Before optimizing LoA in  Figure 13: Changes in production parameters after implementing different scenarios s0(Before change), s1(Expert), s2(Low-cost), s3(Costly but profitable) (a) Changes in production time after implementing different scenarios (s0, s1, s2, s3) in the manufacturing, filling, and packaging workstations (b) Changes in production waste after implementing different scenarios (s0, s1, s2, s3) in the manufacturing, filling, packaging, and QC workstations. (c). Changes in profit after implementing different scenarios (s0, s1, s2, s3). (d). Changes in production quantity after implementing different scenarios (s0, s1, s2, s3) in the manufacturing, filling, and packaging workstations.
chemical industries, the level of environmental pollutants is high and may involve factory closure, heavy penalties, and environmental tax policies. is model helps managers to solve this problem as well. Due to the high cost of full implementation of automation, this model helps managers to create the best performance by choosing the optimal LoA and the dismissal of all workers is also prevented. By using the definition of LoA (2D), it is possible to improve the LoA in different dimensions according to the needs of the organization.

Conclusion
In this paper, a mathematical model is developed to determine the optimal LoA in the chemical industry considering triple aspects of sustainability, namely, economic, environmental, and social aspects. e chemical industry process contains three workstations manufacturing, filling, and packaging. Unit costs of raw materials and overhead costs in manufacturing workstations, filling workstations, and packaging workstations are fuzzy variables. e cost of specialists and workers in each workstation is certain.
is research develops a fuzzy multi-objective model that can improve the LoA by considering any one-dimensional (1D) or two-dimensional (2D) definition of LoA. e danger of transactions of COVID-19 is an important reason to increase LoA to avoid virus infection to workers and customers. To achieve these objectives, an organization was asked to consider the existing financial resources and try to maximize the net present value, minimize the environmental damage and maintain a certain number of production line workers to reduce the social damage of downsizing. In our model, NPV should be positive in the long run. By using a two-dimensional definition of LoA, automation levels can be defined in more precise dimensions, such as the definition proposed by Proud et al. [15]. Also, it should be noted that the proposed levels and dimensions of automation should be homogenous; for example, it cannot be completely manual in one dimension and completely automated in another dimension, and it should also be implementable. A bi-objective mathematical model is developed for this purpose.
e first objective of the model is to minimize the cost function including the cost of materials, labor cost, and overhead costs, and the second objective is to minimize the number of environmental pollutants. e model is solved using Zimmermann max-min approach. Finally, to validate the model, it is implemented in a real case in Iran and three scenarios are developed to cover the uncertainties that are present in some parameters.
(i) e results show that the optimal LoA can be improved with the available resources and the benefits of increasing the LoA can be used for reducing the social damage caused by downsizing the workforce and decreasing the environmental pollutants.
For further research, the authors can refer to optimizing LoA in other industries and optimizing LoA in the supply chain. ey can also consider recycling and optimizing LoA to improve profitability and reduce the related defects. Measuring manufacturing performance parameters can also be measured after improving LoA in organizations.

Data Availability
e data in the article are available in full in the tables.

Conflicts of Interest
e authors declare that they have no conflicts of interest. Discrete Dynamics in Nature and Society