Optimization of the Automated Production Process Using Software Simulation Tools

: The purpose of this article is to point out the need to use software simulation tools in industrial practice to optimize the production process and assess the economic effectiveness of investment, including risk. The goal of the research is to ﬁnd an optimal investment variant to ensure an increase in the production volume of at least 50% and to achieve the maximum economic efﬁciency of the investment, even considering the risk. The article presents a comprehensive approach that enables the achievement of the set research goal. The selection of the optimal version of the investment is carried out in three steps. Firstly, the versions of the investment variants are assessed from the production point of view using the program Tecnomatix Plant Simulation. Subsequently, the versions of the investment variants are assessed from an economic point of view and from a risk point of view. Economic efﬁciency is assessed using the ﬁnancial criteria net present value (NPV), proﬁtability index (PI), and discounted payback period (DPP), and risk analysis is carried out using Monte Carlo simulations. Finally, the accepted outputs are evaluated overall using a multi-criteria method, namely the method of partial order.


Introduction
The investment activity of enterprises, especially investments in production technologies, is a necessary condition for the competitive ability of a production company.Decision-making about investments belongs to strategic decisions, and therefore should be supported by a thorough analysis of production possibilities and production capacity utilization, as well as an analysis of the economic efficiency of the investment, including risk.Several research works deal with the issue of investment decision-making in various aspects.A research work by Freiberg and Scholz [1] presented a comprehensive analysis of the benefits generated by actual investment in modern manufacturing equipment covering all relevant effects related to the investment.A frequent subject of research works on this topic is production optimization aimed at balancing production lines.For example, in a paper by Aguilar et al. [2], parallel assembly lines balancing problem (PALBP) studies were surveyed.Likewise, Van Hop [3] addressed the mixed-model line balancing problem with fuzzy processing time.They developed a fuzzy heuristic to solve the problem based on aggregating fuzzy numbers and combined precedence constraints.In the study by [4], the factors that affect the implementation of intelligent systems in motor production lines were analyzed.The implementation of intelligent systems is heading toward Industry 4.0.In view of Industry 4.0, Zhao et al. [5] also utilized the semi-physical simulation technology of digital twins to verify a proposed design scheme to achieve the balance optimization of a production line.Daneshjo and Malega presented in [6] the usage of automation in an improvement project for the Cinematic production line.A study by [7] Processes 2023, 11, 509 2 of 17 proposed an approach to simultaneously designing products and processes by integrating assembly line balancing, design, and sensitivity analysis of assembly line design.Lu and Yang [8] proposed an approach for solving integrated assembly sequence planning and assembly line balancing with an ant colony algorithm.The research conducted by Liu et al. [9] utilized uncertainty theory to model uncertain demand.The uncertain demand was described using scenario probability and triangular fuzzy numbers.Finally, a new optimization model for mixed-model assembly line balancing was established.In a further publication, Kalayci et al. [10] presented a fuzzy extension of the disassembly line balancing problem with fuzzy task processing times.In a paper by [11], a heuristic assembly line balancing procedure was used to solve the limitations of maximum working team numbers in assembly line stations.The workload smoothness of assembly lines was the subject of further research.The authors proposed an algorithm to solve a two-sided assembly line balancing problem [12] and developed linear mathematical models to assign tasks to the stations of an assembly line [13].Fathi et al. [14] investigated the efficiency of performance measures for minimizing the number of workstations in approaches addressing assembly line balancing problems for both straight and U-shaped lines.Similarly, using the branch and bound algorithm, a simple assembly line balancing problem was solved by Azizo glu and ˙Imat [15].
In decision-making processes about investments in production equipment, software tools designed for production process modelling and simulations are a great help.Through simulations, it is possible to review production capacities in a virtual environment, assess their utilization, and thus improve the quality of investment decisions.In this area, e.g., Chai et al. [16] published their experience with software tools for production process simulation.In their study, a production line 3D visualization monitoring system based on OpenGL modelling was established on a VC++6.0 platform.Production line 3D visualization provided production capacity simulation to optimize the parameter settings.Research on an assembly line using the Arena Simulation program to address the balancing problem was performed by Kayar and Akalin [17].Likewise, to address the balancing issue, Wang et al. [18] proposed a method for optimization of the assembly line based on the genetic algorithm and system simulation, and Bongomin et al. [19] conducted simulationbased optimization under stochastic task times using the Arena simulation software.The paper [20] illustrated how the performance of a balanced manufacturing line could be improved.The authors made the analysis with a discrete-time simulation program, which applied the next-event-time advance mechanism.Based on simulation results, Yang [21] also solved the balancing problem for product processing in a mixed flow assembly line.Song et al. [22] analyzed the efficiency of the assembly line, where the simulation model of the logistics distribution system for the assembly line was built based on Witness software.The simulation and optimization of production lines allow efficiently optimized production parameters.The paper by [23] presented an application of simulation to the virtual design and optimization of an industrial production line.Broad possibilities for analysis and optimization of production processes are offered by the Tecnomatix Plant Simulation software, which is used in practice as well as in research in this field.The Tecnomatix Plant Simulation software tool was used in [24] to detect bottlenecks and deficiencies in the company's production, in [25] to optimize the production system, and, in [26], to design a simulation model for intermediate stock with material flow representation.A comparison analysis of production line simulation results obtained in the two simulation tools FlexSim and Tecnomatix Plant Simulation was performed by [27].
In the context of investment decision-making in production technologies, it is also necessary to mention the financial aspect, which, in practice, is usually considered first.Investments supposed to ensure the viability of the enterprise for many years must be decided after considering all related economic aspects [28] and with the acceptance of a tolerable level of risk.In practice, the effectiveness and riskiness of investments are often assessed by mathematical-economic methods, e.g., using a method based on a genetic neural network [29], based on failure mode and effect analysis (FMEA), or based on the interval-valued intuitionistic fuzzy analytic hierarchy process (AHP) [30].For economic indicators, the application of programs based on Monte Carlo simulations appears to be effective for risk assessment.Khalfi and Ourbih-Tari [31] are among the authors presenting investment risk analysis using Monte Carlo simulations.They proposed a decision-making model to evaluate investment risk using @RISK software.The authors of [32] suggested a method that uses Monte Carlo simulation to estimate the period of value at risk of an investment, and Tobisova et al. [33] analysed investment risk via the net present value parameter.Abba et al. [34] reviewed renewable energy investment risks.The methods used for risk assessment were systems dynamics methods and quantitative methods such as agent-based modelling and Monte Carlo simulation.A nuclear power investment model employing real options theory with the Monte Carlo method was established by Zhu in [35].
This paper aims to introduce a comprehensive, three-stage methodology consisting of production, economic, and risk assessments to support decision-making on investment in production technologies.The assessment of production capacities is based on simulations in Tecnomatix Plant Simulation.Specifically, the workplaces' occupation is analyzed, bottlenecks are identified, and line productivity is evaluated.The economic and risk assessments are based on Monte Carlo simulations.The final investment decision is made using the multi-criteria evaluation method.

Current State of Production
The object of investigation is the production of fire dampers, which form one component of building ventilation.Throttle valves are produced with chimney diameters of 100, 125, 160, 200, 250, and 315 mm.The production description is presented through the product-a throttle valve (Figure 1) with a chimney diameter of 160 mm (TV160)-which has the largest share of the total production volume of throttle valves.
Processes 2023, 11, x FOR PEER REVIEW 3 of 18 neural network [29], based on failure mode and effect analysis (FMEA), or based on the interval-valued intuitionistic fuzzy analytic hierarchy process (AHP) [30].For economic indicators, the application of programs based on Monte Carlo simulations appears to be effective for risk assessment.Khalfi and Ourbih-Tari [31] are among the authors presenting investment risk analysis using Monte Carlo simulations.They proposed a decisionmaking model to evaluate investment risk using @RISK software.The authors of [32] suggested a method that uses Monte Carlo simulation to estimate the period of value at risk of an investment, and Tobisova et al. [33] analysed investment risk via the net present value parameter.Abba et al. [34] reviewed renewable energy investment risks.The methods used for risk assessment were systems dynamics methods and quantitative methods such as agent-based modelling and Monte Carlo simulation.A nuclear power investment model employing real options theory with the Monte Carlo method was established by Zhu in [35].This paper aims to introduce a comprehensive, three-stage methodology consisting of production, economic, and risk assessments to support decision-making on investment in production technologies.The assessment of production capacities is based on simulations in Tecnomatix Plant Simulation.Specifically, the workplaces' occupation is analyzed, bottlenecks are identified, and line productivity is evaluated.The economic and risk assessments are based on Monte Carlo simulations.The final investment decision is made using the multi-criteria evaluation method.

Current State of Production
The object of investigation is the production of fire dampers, which form one component of building ventilation.Throttle valves are produced with chimney diameters of 100, 125, 160, 200, 250, and 315 mm.The production description is presented through the product-a throttle valve (Figure 1) with a chimney diameter of 160 mm (TV160)-which has the largest share of the total production volume of throttle valves.Currently, the production of TV160 throttle valves takes place on a non-automated production line with Welding, Assembly1, and Assembly2 workstations.One worker works at each workstation of the production line, i.e., there are three in total.The sequence of production operations and their durations according to individual workplaces on the production line are recorded in Table 1.Between the workplaces, there are Buffer1 and Buffer2 interoperation warehouses with a capacity of 10 units, which serve to dampen fluctuations in the line's activity.At the end of the line there is Buffer3, where the completed production is collected.Currently, the production of TV160 throttle valves takes place on a non-automated production line with Welding, Assembly1, and Assembly2 workstations.One worker works at each workstation of the production line, i.e., there are three in total.The sequence of production operations and their durations according to individual workplaces on the production line are recorded in Table 1.Between the workplaces, there are Buffer1 and Buffer2 interoperation warehouses with a capacity of 10 units, which serve to dampen fluctuations in the line's activity.At the end of the line there is Buffer3, where the completed production is collected.
The TV160 production process was recorded using the Tecnomatix Plant Simulation program, as shown in Figure 2, including a graph of the utilization of individual workplaces.The simulation time is eight hours, i.e., one work shift.
Based on the facts mentioned above, it can be concluded that 208 TV160 products are produced on the given production line in one eight-hour work shift.This means the hourly production of the line is 26 pieces.It can be concluded that the line is relatively well-balanced, and the workplaces of the line are utilized 80-100%.The bottleneck is the Welding workplace, where the total operating time is 138 s.Interoperation warehouses sized for ten pieces are practically unused, and no stocks larger than one piece are created.Based on the facts mentioned above, it can be concluded that 208 TV160 products are produced on the given production line in one eight-hour work shift.This means the hourly production of the line is 26 pieces.It can be concluded that the line is relatively well-balanced, and the workplaces of the line are utilized 80-100%.The bottleneck is the Welding workplace, where the total operating time is 138 s.Interoperation warehouses sized for ten pieces are practically unused, and no stocks larger than one piece are created.

Methodical Procedure of Solution
The problem in the production of throttle valves is the long-term failure to achieve the planned production, i.e., the inability of the company to meet demand in full.Eliminating this problem is possible by investing in the automation of selected production operations.The introduction of automation into production stabilizes and increases product quality, improves the organization and management of production, humanizes human work, and saves resources.Thus, automation enables the saving of materials, reducing of rejects, reducing of the consumption of technological substances, etc.In addition to the benefits resulting from production cost savings, automation also has secondary effects.These include labor reduction in production preparation and in operative production planning and shortening of the production time of products.

Methodical Procedure of Solution
The problem in the production of throttle valves is the long-term failure to achieve the planned production, i.e., the inability of the company to meet demand in full.Eliminating this problem is possible by investing in the automation of selected production operations.The introduction of automation into production stabilizes and increases product quality, improves the organization and management of production, humanizes human work, and saves resources.Thus, automation enables the saving of materials, reducing of rejects, reducing of the consumption of technological substances, etc.In addition to the benefits resulting from production cost savings, automation also has secondary effects.These include labor reduction in production preparation and in operative production planning and shortening of the production time of products.
The procedure for solving the problem is documented in Figure 3 and defined in more detail in the following text.
The main goal of the research is to find the optimal solution related to investment in the automation of selected production operations of the line.The intention is to ensure an increase in production volume by at least 50% (from 208 pcs/8 h to at least 312 pcs/8 h) and, at the same time, ensure that the investment is economically efficient, even taking into account the risk.The company is not considering the transition from one-shift to two-shift operations.
Assessment of the production capacities of investment variants.A workplace simulation model was created for each investment variant (V1, V2) to assess production capacities.Subsequently, the production process was simulated.The processing time of individual workplaces was entered based on the sum of the times of operations performed at each workplace (1).
where Tp is cumulative processing time, t i is the operating time of the i-th operation, and n is the number of operations at the workplace.
Processes 2023, 11, x FOR PEER REVIEW 5 of 18 The procedure for solving the problem is documented in Figure 3 and defined in more detail in the following text.
The main goal of the research is to find the optimal solution related to investment in the automation of selected production operations of the line.The intention is to ensure an increase in production volume by at least 50% (from 208 pcs/8 h to at least 312 pcs/8 h) and, at the same time, ensure that the investment is economically efficient, even taking into account the risk.The company is not considering the transition from one-shift to twoshift operations.

Assessment of the production capacities of investment variants.
A workplace simulation model was created for each investment variant (V1, V2) to assess production capacities.Subsequently, the production process was simulated.The processing time of individual workplaces was entered based on the sum of the times of operations performed at each workplace (1).
where Tp is cumulative processing time, ti is the operating time of the i-th operation, and n is the number of operations at the workplace.Individual operation times and the cumulative processing time across every workplace are presented in Table 2.The simulation output is a chart showing the utilization of workplaces in % (work time, waiting, and blocking time).Another output is the production volume produced during the simulation, i.e., in an 8 h shift, obtained from workplace statistics.Individual operation times and the cumulative processing time across every workplace are presented in Table 2.The simulation output is a chart showing the utilization of workplaces in % (work time, waiting, and blocking time).Another output is the production volume produced during the simulation, i.e., in an 8 h shift, obtained from workplace statistics.
Assessment of the economic efficiency of investment variants, including risk.In the first step, economic efficiency evaluation of the proposed investment variants is carried out based on the incremental method.This means that the expected effects and resources expended due to the automation of selected operations on the production line are expressed as the difference between the expected values when implementing the investment option and the values achieved with current non-automated production.The calculation is carried out in MS Excel using a deterministic financial model.Economic efficiency is assessed using three dynamic financial criteria, namely net present value (NPV), profitability index (PI) and discounted payback period (DPP).The mentioned financial criteria support the same decision when evaluating investment variants-to accept or not to accept an investment variant.Their numerical results are different because each financial criterion assesses the investment variant from a different point of view.NPV expresses the profitability of an investment option in an absolute amount for the entire life of the project.PI monitors the return on one euro of capital expenditure.DPP monitors the liquidity of the investment option, i.e., the time for which it will be repaid.The NPV calculation method is given by relation (2), PI relation (3), and DPP relation (4).
The values of the financial criteria are determined by the increase in annual cash flow (CF) from the investment variants converted to the present value at the discount rate (d r ), and the amount of one-time investment costs (IC).N is the period of the economic life of the investment, n are the individual years of the economic life of the investment.The relationship for calculating the annual CF from the operating activity is expressed as follows (5): where EBITDA is the annual increase in the economic result before interest, taxation, and depreciation, ODP is the annual depreciation, and d is the coefficient of the income tax rate.EBITDA is influenced by production and sales activities.The calculation method for individual investment variants takes into account the number of robots included in the production line, the number of production workers involved in the production process, the annual production volume of the TV160, and more.The calculation of the EBITDA indicator, specifically its annual increase, is expressed by relation (6): where c is the price per unit of production, Q is the increase in annual production volume, N FIX is the annual increase in fixed costs, and N VAR is the annual increase in variable costs.The annual increase in fixed costs includes the increase in annual depreciation and the increase in annual costs for voluntary insurance of production equipment.The annual increment of variable costs is the sum of the annual increments of direct material consumption, auxiliary material consumption, personnel costs, robot spare parts costs, robot repair and maintenance costs, and other variable costs.The second step of economic efficiency evaluation considers risk in the proposed investment variants assessment.This is a stochastic approach.Risk consideration will be done using the Monte Carlo method in the Crystal Ball software product.In this step, it is necessary to transform the deterministic financial model into a stochastic one.It is necessary to determine input risk factors and their probability distributions, establish statistical dependencies between the selected risk factors, and define the probability distribution of the output variable (for the selected financial criterion).The simulation process itself is defined by the number of simulation steps, and it is true that the higher the number of repetitions, the higher the accuracy of the results.The simulation results refer to the output variable, which will be the mean value of the financial criterion NPV.The simulation results are generated in numerical form, in the form of statistical characteristics, and in graphical form as a histogram, trend chart, sensitivity chart, overlay chart, and others.
Complex evaluation of the investment variants.It is advisable to choose the optimal version of the investment variant using a multi-criteria evaluation using the partial ranking method.The method is based on determining the sum of the partial order of the production and financial criteria for the compared versions of the investment variants.Partial orders depend on the desired tendency of the criterion value.If the value is descending (for example, DPP), then the criterion with the lowest value will come first.If, on the other hand, the value of the criterion is to increase (for example, NPV), the highest value of the criterion is best-placed.The most advantageous variant is the one whose sum of partial orders is the lowest.

Proposals for Investment Variants of Production Line Automation
Achieving the research goal, i.e., increasing the production volume by at least 50% in the production of TV160, is possible only in the production operations of welding and sealing.For this reason, modification of the current TV160 production process is necessary.The Welding workplace will be abolished, and two new workplaces will be created, namely the Pre-assembly workplace and the Robot workplace.At the Pre-assembly workplace, the operations of placing the chimney and placing it in the fixture will take place.The sidewall riveting operation, which took place at the Welding workplace, has been cancelled.The Robot workplace will replace the welding performed at the workplace Welding and the sealing of the chimney performed at the Assembly1 workplace.
Two investment variants for the automated production line (for one and two robots) are proposed in two versions according to the number of workers (three and four).The automation of the production process is ensured at the Robot workplace.Workplaces Assembly1 and Assembly2 do not change.That means the automated production line consists of four workplaces: Pre-assembly, Robot, Assembly1, and Assembly2.The sequence of operations performed at the listed workplaces, including their operating times, is presented in Table 2.
The determination of times for welding and sealing operations performed at the Robot workplace is based on the input data listed in Table 3.The time of the welding operation is the product of the number of welds and the welding time of one weld.The distance between the two welds is set at 77 mm.The time of the sealing operation is determined by the length of the sealant used and the time required to apply it per unit length.Individual investment variants are simulated in the Tecnomatix Plant Simulation software.The simulation output is a chart showing in percentages the time utilization of workplaces, in terms of work time, waiting time, and blocking time.Another output is the volume produced during the simulation, i.e., in 8 h, obtained from workplace statistics.The description of individual investment variants is given in the following text.
Variant 1a (V1a).In investment variant 1 (Figure 4), the Robot workplace consists of one robot performing both welding and sealing operations.One worker works at the line workplaces Pre-assembly, Assembly1, and Assembly2, so there are three in total.The usability of the robotic workplace is 72%, and the workplace is blocked (fully occupied but not processing) for 28% of the available working time.The bottleneck is the Assembly2 workplace.Low utilization is found at the Pre-assembly workplace, which is only 34% utilized.The production of the line is 292 pcs/8 h, which does not meet the requirement for a 50% increase in the line production capacity.Variant 1b (V1b).In investment variant V1b, the welding and sealing operation is performed by a single robot in the robotic workplace.There are four workers on the line.By adding a fourth worker who works at the Assembly1 and Assembly2 workplaces (in the simulation, a 50%/50% work distribution is considered), the bottleneck of the line will be moved to the robotic workplace (see Figure 5).The utilization of the Assembly workplaces is more than 80%; the Pre-assembly workplace is still only 47% loaded.With the addition of the fourth worker, the line productivity is increased to 434 pcs/8 h, thus meeting the requirement to increase production by at least 50%.In this variant, production will increase by approximately 109%.Variant 1b (V1b).In investment variant V1b, the welding and sealing operation is performed by a single robot in the robotic workplace.There are four workers on the line.By adding a fourth worker who works at the Assembly1 and Assembly2 workplaces (in the simulation, a 50%/50% work distribution is considered), the bottleneck of the line will be moved to the robotic workplace (see Figure 5).The utilization of the Assembly workplaces is more than 80%; the Pre-assembly workplace is still only 47% loaded.With the addition of the fourth worker, the line productivity is increased to 434 pcs/8 h, thus meeting the requirement to increase production by at least 50%.In this variant, production will increase by approximately 109%.
By adding a fourth worker who works at the Assembly1 and Assembly2 workplaces (in the simulation, a 50%/50% work distribution is considered), the bottleneck of the line will be moved to the robotic workplace (see Figure 5).The utilization of the Assembly workplaces is more than 80%; the Pre-assembly workplace is still only 47% loaded.With the addition of the fourth worker, the line productivity is increased to 434 pcs/8 h, thus meeting the requirement to increase production by at least 50%.In this variant, production will increase by approximately 109%.Variant 2a (V2a).In investment variant V2a, there are two robots in a parallel arrangement for welding and sealing operations.Robots perform both operations.There are three workers on the line.The bottleneck is again the Assembly2 workplace (see Figure 6).The line, arranged in this way, can produce 291 pcs/8 h.The utilization of the Assem-bly1 workplace is more than 80%.The Pre-assembly workplace is only 35% loaded.The parallel robotic workstation is only 36% utilized, 49% waiting, and 15% of the time blocked.In the V2a variant, production will increase by approx.109% (see Table 4).Variant 2a (V2a).In investment variant V2a, there are two robots in a parallel arrangement for welding and sealing operations.Robots perform both operations.There are three workers on the line.The bottleneck is again the Assembly2 workplace (see Figure 6).The line, arranged in this way, can produce 291 pcs/8 h.The utilization of the Assembly1 workplace is more than 80%.The Pre-assembly workplace is only 35% loaded.The parallel robotic workstation is only 36% utilized, 49% waiting, and 15% of the time blocked.In the V2a variant, production will increase by approx.109% (see Table 4).Variant 2b (V2b).Investment variant V2b has two robots at the Robot workplace connected in parallel (see Figure 7), while both robots perform welding and sealing operations.There are four workers on the production line (one Pre-assembly, one at Assembly1, one at Assembly2, and the fourth worker helps at Assembly1) and two workplaces in a ratio of 50%/50%.The number of robots and the distribution of workers will make it possible to increase the performance of the robotic workplace by 100%, which increases line productivity to 440 pcs/8h.Workplace utilization is more than 80% in the case of Assem-bly1.The line bottleneck is Assembly2, whose capacity is fully utilized.The robotic workplace and the Pre-assembly workplace are about 50% utilized.Compared to the current state, production increases by approximately 112%, thus fulfilling the requirement to increase production by at least 50%.Variant 2b (V2b).Investment variant V2b has two robots at the Robot workplace connected in parallel (see Figure 7), while both robots perform welding and sealing operations.There are four workers on the production line (one Pre-assembly, one at Assembly1, one at Assembly2, and the fourth worker helps at Assembly1) and two workplaces in a ratio of 50%/50%.The number of robots and the distribution of workers will make it possible to increase the performance of the robotic workplace by 100%, which increases line productivity to 440 pcs/8 h.Workplace utilization is more than 80% in the case of Assembly1.The line bottleneck is Assembly2, whose capacity is fully utilized.The robotic workplace and the Pre-assembly workplace are about 50% utilized.Compared to the current state, production increases by approximately 112%, thus fulfilling the requirement to increase production by at least 50%.Variant 2b (V2b).Investment variant V2b has two robots at the Robot workplace connected in parallel (see Figure 7), while both robots perform welding and sealing operations.There are four workers on the production line (one Pre-assembly, one at Assembly1, one at Assembly2, and the fourth worker helps at Assembly1) and two workplaces in a ratio of 50%/50%.The number of robots and the distribution of workers will make it possible to increase the performance of the robotic workplace by 100%, which increases line productivity to 440 pcs/8h.Workplace utilization is more than 80% in the case of Assem-bly1.The line bottleneck is Assembly2, whose capacity is fully utilized.The robotic workplace and the Pre-assembly workplace are about 50% utilized.Compared to the current state, production increases by approximately 112%, thus fulfilling the requirement to increase production by at least 50%.From the production point of view, it can be concluded that the requirement to increase the production volume by at least 50% was met by investment variants 1b and 2b (see Table 4).Of the listed investment variants, Variant 2b enables the greatest increase in production volume.The production volume of TV160 increases from 207 pcs/8 h (current state) to 440 pcs/8 h, which means an increase of 112.6%.From the production point of view, it can be concluded that the requirement to increase the production volume by at least 50% was met by investment variants 1b and 2b (see Table 4).Of the listed investment variants, Variant 2b enables the greatest increase in production volume.The production volume of TV160 increases from 207 pcs/8 h (current state) to 440 pcs/8 h, which means an increase of 112.6%.

Economic Evaluation of Investment Variants of the Production Line
Deterministic approach to evaluating the economic efficiency of investment variants.The company will incur investment costs by procuring robots.Their amount is influenced by the purchase price of the robot/robots and the costs related to their procurement.Selected input data necessary for defining a deterministic financial model and subsequent assessment of the economic efficiency of investment variants are presented in Table 5.
When calculating the financial criteria, which are used to assess the economic efficiency of investment proposals, the annual time fund of the robot in the amount of 2014 h is considered.Its utilization is expected to be 85%.Robot, according to Act no.595/2003 Coll. on income tax, belongs to the second depreciation group with a depreciation period of six years, which is also the lifetime of the investment.Economic calculations consider the price of TV160 at the amount of EUR 9.50 per piece, the hourly labor costs per worker at the amount of EUR 10.50, and the material costs of TV160 at the amount of EUR 10.00 per two pieces.The resulting values of the financial criteria, according to the proposed investment variants V1, V2 and their versions, are recorded in Table 6.It can be concluded that, from an economic point of view, except for the investment variant V2a, the other proposed versions of the investment variants are acceptable.The best economic results were achieved with version V1b.Its implementation will achieve profitability in the form of monetary income for a period of six years in the amount of EUR 459,578, profitability per one euro of invested investment costs in the amount of EUR 4.31, and a payback period of 1.22 years.The second in line is the investment variant V2b.Marginal economic results are achieved by the investment variant V1a.Stochastic approach to evaluating the economic efficiency of investment variants.Taking risk into account in evaluating the economic efficiency of investment variants V1, V2, and their versions is implemented using the Monte Carlo simulation method.Individual input risk factors defined using statistical characteristics and probability distributions are listed in Table 7.The selection of the most advantageous investment variant, taking into account the risk, is based on the mean value of the financial criterion NPV (mean NPV).The probability distributions of versions of investment variants V1 and V2 are determined by a mix of input risk factors, and their courses are shown in Figure 8.They provide information about the expected value of NPV and the probability of not achieving it or exceeding it.

Product price
EUR/pc Likeliest 9.50; Min.9.00; Max.10.50 Triangular The selection of the most advantageous investment variant, taking into account the risk, is based on the mean value of the financial criterion NPV (mean NPV).The probability distributions of versions of investment variants V1 and V2 are determined by a mix of input risk factors, and their courses are shown in Figure 8.They provide information about the expected value of NPV and the probability of not achieving it or exceeding it.Based on the data presented in Table 8, it can be concluded that the symmetry of probability distributions of NPV is approximately the same for all versions of the investment variants.The NPV probability distributions for all versions are positive, that is, right-skewed, which means that higher values prevail over lower values.Due to the fact Based on the data presented in Table 8, it can be concluded that the symmetry of probability distributions of NPV is approximately the same for all versions of the investment variants.The NPV probability distributions for all versions are positive, that is, rightskewed, which means that higher values prevail over lower values.Due to the fact that the mean NPVs of the individual versions of the investment variants differ significantly, the rule of mean value and the coefficient of variation are used as decision criteria for determining their advantage.According to the mentioned rule, the version of the investment option that has a higher mean value for the financial criterion and a lower variability is better.The use of this rule confirms that the investment variant version V1b is dominant over the other versions of the investment variants, because it is better in both the considered characteristics compared to them.The NPV_V1b histogram is shown in Figure 9.This documents that, with a probability of 39.54%, higher values of the financial criterion NPV_V1b than in its base case will be achieved.From the sensitivity graph (see Figure 10), it is possible to determine the key risk factors both in terms of the overall forecast uncertainty and their impact on the final value of the financial criterion NPV_V1b.It is clear from the graph that the uncertainty of the NPV_V1b forecast is positively influenced by the price of the TV160 product to the extent of 53.8% and the amount of the nominal time fund to the extent of 0.5%.The uncertainty of the NPV_V1b forecast is negatively affected in the amount of 40.1% by material consumption for the TV160 product, in the amount of 5.1% by the investment costs associated with the procurement of the robot, and in the amount of 0.4% by hourly labor costs.From the sensitivity graph (see Figure 10), it is possible to determine the key risk factors both in terms of the overall forecast uncertainty and their impact on the final value of the financial criterion NPV_V1b.It is clear from the graph that the uncertainty of the NPV_V1b forecast is positively influenced by the price of the TV160 product to the extent of 53.8% and the amount of the nominal time fund to the extent of 0.5%.The uncertainty of the NPV_V1b forecast is negatively affected in the amount of 40.1% by material consumption for the TV160 product, in the amount of 5.1% by the investment costs associated with the procurement of the robot, and in the amount of 0.4% by hourly labor costs.The tornado chart (see Figure 11) shows the strengths of the individual influences of the input variables on the change in the financial criterion NPV_V1b.It shows the variables according to their impact on NPV_V1b for a considered change of ±10%.In the chart, the variables are arranged according to the size of their impact on the financial criterion NPV_V1b.Among the most serious factors causing the NPV_V1b deviation are, in order, product price, material consumption, nominal time fund, investment costs, and labor costs.Less important factors include electricity consumption, electricity price and planned downtime for repair and maintenance of robots.The tornado chart (see Figure 11) shows the strengths of the individual influences of the input variables on the change in the financial criterion NPV_V1b.It shows the variables according to their impact on NPV_V1b for a considered change of ±10%.In the chart, the variables are arranged according to the size of their impact on the financial criterion NPV_V1b.Among the most serious factors causing the NPV_V1b deviation are, in order, product price, material consumption, nominal time fund, investment costs, and labor costs.Less important factors include electricity consumption, electricity price and planned downtime for repair and maintenance of robots.
Due to the effect of taking into account the risk in the most economically advantageous investment variant V1b (see Table 9), there was a decrease in mean NPV_V1b compared to the expected value of NPV_V1b by EUR 23,352, which represents a decrease of 5.1%.
the input variables on the change in the financial criterion NPV_V1b.It shows the variables according to their impact on NPV_V1b for a considered change of ±10%.In the chart, the variables are arranged according to the size of their impact on the financial criterion NPV_V1b.Among the most serious factors causing the NPV_V1b deviation are, in order, product price, material consumption, nominal time fund, investment costs, and labor costs.Less important factors include electricity consumption, electricity price and planned downtime for repair and maintenance of robots.Due to the effect of taking into account the risk in the most economically advantageous investment variant V1b (see Table 9), there was a decrease in mean NPV_V1b compared to the expected value of NPV_V1b by EUR 23,352, which represents a decrease of 5.1%.

Complex Evaluation of Investment Variants
Selection of the optimal version of the investment variant is carried out using a multicriteria evaluation, using the partial order method.Only those versions of investment variants that meet three conditions are included in the multi-criteria evaluation: The production condition: i.e., they will ensure an increase in production volume by at least 50%.Only the V1b and V2b versions of the investment variants met this condition.
The condition of economic efficiency, determined by acceptable values of the financial criteria NPV, PI, and DPP.In this case, the condition was met by the versions of investment variants V1a, V1b, and V2b.Version V1a achieves only marginally acceptable values for the financial criteria and at the same time does not meet the production condition; therefore, it is not included in the multi-criteria evaluation.
The condition of economic efficiency, with consideration of risk using the rule of mean NPV and the coefficient of variation.The best values are achieved by the versions of variants V1b and V2b.
The V1b version was determined to be the optimal investment variant based on a multi-criteria evaluation.It achieves the best economic parameters through dynamic financial criteria.In terms of production parameters, it allows one robot to achieve almost identical production volume values as with the V2b version with two robots.
Case studies that have presented the optimization of production processes, either using simulation tools or not, were mentioned in the introduction.This was mainly one critical assessment focused on the production or economic aspects of the project, which either did not include uncertainty in the assessment or only in connection with the economics assessment.In the presented study, the production, finances, and risk aspects are in a complex multi-criteria evaluation (Table 10).However, other authors of scientific papers also point to the need for risk integration in multi-criteria evaluation, albeit in different contexts.For example, the authors of [36] presented a fuzzy multi-criteria decision approach in the selection of a sustainable risk management strategy.Similarly, [37] used a multicriteria assessment system when choosing the method of operational risk management, and [38] performed a multi-criterion analysis based on integrated process-risk optimization.The sum of the partial orders 9 14 Overall order 1. 2.

Conclusions
Automation of production cannot be done without investment in new equipment or technologies.Although these will expand production possibilities and improve the quality of production operations, on the other hand, they require significant financial resources and imply new risks.A decision on the optimal investment variant is based on an assessment of the investment from several sides.Computer simulation programs are also a supporting tool in decision-making, the potential of which is very large at the moment.
In this article, a comprehensive approach to investment decision-making using simulations was presented through a case study of automation of the production line workplace.Two variants of the investment were proposed according to the equipment of the workplace, with one or two robots, and for each variant two versions of the line service (with three or four workers) were evaluated.Tecnomatix Plant Simulation was used to assess the line in terms of production capacity.The capacity of the line, the bottleneck, and the workload of the workplaces were checked on the simulation models for the individual investment variants.The Crystal Ball program was used for the economic assessment of the investment, including risk factors, where the economic efficiency of the investment was evaluated and risk factors were analyzed using Monte Carlo simulations.The investment variants were assessed first from the point of view of production, where the goal was to achieve production 50% higher than current production.Subsequently, they were evaluated from an economic point of view through the financial indicators of NPV, PI, and DPP.From the point of view of risk, the statistical indicators mean value, standard deviation, and coefficient of variation were selected.The outputs from the individual phases of the investment assessment were ultimately evaluated as a whole by a multi-criteria assessment of the investment variants using the method of partial order.
The primary contribution of the paper is the application of a complex approach to the design of production systems using simulations.However, simulation programs, as a progressive tool in decision-making, were applied only to a limited extent, which is considered a limit of the present research.Their possibilities are much broader, including, e.g., analysis and optimization of all production logistics processes, including interoperation storage and transport.In the case of Monte Carlo simulations, which are part of the Crystal Ball program, it is also appropriate to use the OptQuest module to optimize the process while taking risk into account.A broader application of simulation tools is the intention of future research in this issue.

Figure 1 .
Figure 1.3D view of a throttle valve.

Figure 1 .
Figure 1.3D view of a throttle valve.
"Labeling" transferred to an off-line worker (operator)The TV160 production process was recorded using the Tecnomatix Plant Simulation program, as shown in Figure2, including a graph of the utilization of individual workplaces.The simulation time is eight hours, i.e., one work shift.

Figure 3 .
Figure 3. Procedure for defined problem solving.

Figure 3 .
Figure 3. Procedure for defined problem solving.

Figure 8 .
Figure 8. Probability distributions of NPV versions of investment variants.

Figure 8 .
Figure 8. Probability distributions of NPV versions of investment variants.

Figure 10 .
Figure 10.Graph of sensitivity of risk factors to NPV_V1b.

Table 1 .
The current state of production of TV160-workplaces, production operations, and operating times.

Table 1 .
The current state of production of TV160-workplaces, production operations, and operating times.

Table 2 .
Proposed state of production-workplaces, production operations, and total times at workplaces in seconds.

Workplace Assembly1 Production operations Stickers incl. submissions Chimney sealing Insert the seal Flap installation Sealing the cover Closing the lid
Note: Production operations marked in red are removed compared to the current state.

Table 3 .
Input data for the calculation of times for welding and sealing operations at the Robot workplace.

Table 4 .
Comparison of the production parameters of the proposed investment variants with the current state.

Table 5 .
Selected input data for assessing the economic efficiency of investment variants.

Table 6 .
Financial criteria of economic efficiency-deterministic approach.

Table 7 .
Probability distributions and statistical characteristics of risk factors.

Table 8 .
Selected statistical characteristics of NPV versions of investment variants.

Table 9 .
Comparison of NPV_V1b and mean NPV_V1b values.

Table 10 .
Selection of the optimal version of the investment variant by the method of partial orders.