Deep Learning Approach of Energy Estimation Model of Remote Laser Welding
Deep Learning Approach of Energy Estimation Model of Remote Laser Welding
Deep Learning Approach of Energy Estimation Model of Remote Laser Welding
net/publication/333018026
CITATION READS
1 136
3 authors, including:
Jumyung Um
Technologie-initiative SmartFactory KL e.V.
31 PUBLICATIONS 303 CITATIONS
SEE PROFILE
Some of the authors of this publication are also working on these related projects:
MAYA: Methodologies and Tools for Design, Management and Simulation of CPS-based Factories View project
All content following this page was uploaded by Jumyung Um on 11 May 2019.
Abstract: Due to concerns about energy use in production systems, energy-efficient processes have
received much interest from the automotive industry recently. Remote laser welding is an innovative
assembly process, but has a critical issue with the energy consumption. Robot companies provide only
the average energy use in the technical specification, but process parameters such as robot movement,
laser use, and welding path also affect the energy use. Existing literature focuses on measuring
energy in standardized conditions in which the welding process is most frequently operated or on
modularizing unified blocks in which energy can be estimated using simple calculations. In this
paper, the authors propose an integrated approach considering both process variation and machine
specification and multiple methods’ comparison. A deep learning approach is used for building
the neural network integrated with the effects of process parameters and machine specification.
The training dataset used is experimental data measured from a remote laser welding robot producing
a car back door assembly. The proposed estimation model is compared with a linear regression
approach and shows higher accuracy than other methods.
Keywords: remote laser welding; energy-efficient process; machine learning; welding process;
neural network
1. Introduction
Manufacturing activities account for one-third of the total energy consumption and CO2 emissions
worldwide. Thus, energy-efficient production has been the focus of research over the last decade
and is considered one of the major decision-making attributes of manufacturing today. This paper
takes into account both the energy efficiency of the manufacturing process and the system. In mass
production industries energy-efficient processes have received much interest, for example from the
automotive industry.
The energy use depends on the process technology of the machine conducting each process.
This paper concerns energy efficiency in welding, specifically with RLW (Remote Laser Welding) as
compared to spot welding. Figure 1 illustrates one aspect of movement differences between spot
welding and remote welding. A spot welding robot visits all positions of stitches (Figure 1 (Left))
and consumes energy for movement to welding positions as well as for the welding itself whereas in
remote welding the robot performs welding at a distance and hence can use simpler movement paths
(Figure 1 (Right)).
Bearing this in mind, it can be seen that the process technology influences the energy consumption
of other components in a machine. To reduce the total energy use, measurement and estimation
of the whole process are required. As pointed out in [1], when new processes are introduced into
factory production lines, energy consumption is a significant factor in machine selection in addition to
productivity needs. In this respect, RLW is interesting because it dramatically reduces both lead time
and energy use in the assembly phase of production. Because of these benefits, RLW is of interest in
mass production, notably in automotive assembly lines. However, it is necessary to be able to evaluate
accurately potential benefits in order to see when the high initial investment is defrayed with savings.
The current approach to new machine selection is to compare the technical specifications of energy use,
which is a simple static number. The specification usually seldom shows the variation of the energy
profile occurring in the process of specific products.
Figure 1. (Left) Short-distance welding path; (Right) remote welding path. The differences between
short-distance and remote welding.
It is necessary for approaches to estimate the energy profile occurring when a welding robot
attaches a specific set of stitches, as well as the energy efficiency of the overall processes. To do this,
an energy estimation model is required before a new welding robot is introduced into the assembly
line. The literature on laser welding focuses on measuring the energy profile of welding a specific
assembly. Most research about welding energy consumption concerns estimation of the robot kinetic
energy. Currently, there are many energy models of robot arms in the literature. The major energy
consumers of RLW are the robot arm, the laser machine because laser generation energy is higher
than conventional lasers, and cooling equipment. The ratio of robot motion among major consumers
is relatively smaller than for a spot welding robot, which is measured with only kinetic energy use.
For this reason, the energy use of each assembly process is different from others because the laser
equipment is an additional variable causing fluctuations in the energy profile compared with the spot
welding process.
The mechanism of the laser machine consists of the series of laser generation, concentration,
delivery, and focusing, which are relatively more complicated than robot motion. The energy
consumption model also reflects the mechanism and main consumption by using weighting factors.
Using rule-based approaches, it is difficult to develop common models for different processes and
machines. Thus, an energy estimation model for laser machines is proposed in this paper using a
machine learning approach. In particular, a deep learning model is utilized for training and correlation
with process parameters and extracted feature data.
This paper proposes a comprehensive energy estimation model for RLW and its energy efficiency
model. In Section 2, the basic concept and existing models for energy estimation are introduced.
Section 3 describes the deep learning approach to simulate estimated energy use, followed by Section 4
explaining the results of experiments with the welding process for an automotive rear door assembly.
Then, in Section 5, the difference between the proposed approach and conventional models is described.
Conclusions are given in Section 6.
2. Literature Survey
In this section, previous models of energy efficiency of manufacturing resources, which are
specifically for machine tools, are discussed. Then, some methods and strategies for estimating energy
Energies 2019, 12, 1799 3 of 20
use, which is a key input for calculation of energy efficiency, are introduced from aspects of energy
efficiency and energy consumption of functional analysis and process analysis. Finally, the use of
machine learning to estimate energy profiles is evaluated in terms of pros and cons.
• How do we count the amount of work? Typically, the work of metal cutting is measured by the
removal volume. This factor is directly correlated to the feed rate, depth of cut, and cutting speed
of milling and turning.
• How do we measure the productivity of each shape of a stitch? Each stitch has a shape, i.e., a line,
S-line, C-curve, etc., with different depths and widths. The energy use depends on all aspects.
• The velocity of the robot arm is another factor, but with the same effect on energy consumption as
for metal cutting.
Energies 2019, 12, 1799 5 of 20
• Even though the same process parameters of laser welding are applied, the absorption rate of
different materials is a factor that should be considered.
Figure 2. Overall process of the eco-efficiency of the Remote Laser Welding (RLW) process.
Figure 3. Compensation between the functional approach and the process approach
• Function level: This approach classifies functions into predefined general parts that consist of
machine operation, workpiece manipulation, process conditioning, tool handling and die change,
machine cooling and heating, and recyclables and waste handling. In the case of RLW, the
machine operation function is laser generation.
• Component level: RLW consists of a 24-V power supply, CNC air-conditioning, drive power
supply, two cooling fans, a generator, and a cooling pump, which are described in the e-scheme
of RLW.
Each component was assigned to predefined functions according to ISO 14955-1 as shown in
Figures 4 and 5. Representing the connected load of the component, the e-scheme was used for
selecting components from the whole machine. The estimation results showed which components
were significant energy consumers. The proposed functional analysis was that machining and cooling
conditioning were major consumers. Physically, machining functions were conducted by a robot arm
and laser generator. For cooling conditioning, the chiller was the major part. The sub-components
of major consumers were joint motors, laser diodes, laser disks, and the water compressor. The laser
scanner of the robot arm was a relatively minor consumer compared with joint motors.
Energies 2019, 12, 1799 7 of 20
Figure 4. Cluster matrix for functional analysis according to ISO 14955 [19].
The structure of the laser source influences the energy produced by laser generation. For the case
study, a Trumpf TruDisk 4002 was used. This is a high-power laser source concentrating 176 low-power
diode laser sources through four disks [20]. A single stack consisted of 11 diode bars generating 100
W of laser power. Sixteen stacks transferred the laser into 4fourserial disks chilled by the coolants.
The overall structure is illustrated in Figure 6. Each disk was connected with the optical structure to
collect small diode laser sources. Collected laser power reached 4 KW and a 1-Mfocal length. This
performance enabled welding robots to cover a wide area with only laser scanner movement [21].
Energies 2019, 12, 1799 8 of 20
With this performance specification, the proposed energy estimation model was designed with
the following assumptions.
• Each laser diode was activated when a laser machine was ready to launch. Only optical control,
such as mirror blocking, was conducted without deactivating the didoes if the laser was not used.
• Otherwise, a laser machine was equipped with a chiller to cool disks by coolant. The chiller
consumed more energy if laser disks concentrated more laser power than before. This means that
the chiller energy increased if laser output power rose.
monitoring program. The power cell was connected with the power supply of the robot controller and
the laser machine to measure the current and voltage used. The calculated power was transferred to
the data acquisition device, which converted the power signal into digital values sampled at a rate of
5000 Hz. The energy monitoring program stored the resulting energy profile with a time series dataset.
Feature Engineering
Feature engineering is a means for utilizing the process parameters so as to estimate the detail
of the energy profile by overcoming the limitation of localization. In the case of RLW, the given
parameters were smaller than for the other cases of big data problems. It is necessary to focus on
the specific data and find the important relationships. As mentioned in the analysis of the electric
connection, cooling and generation of the laser were the major components. Laser generation showed
rapid jumps in energy consumption. Cooling followed the motor behaviors to reduce the temperature
of disks and diodes. The key features of laser energy consumption were the current status of the laser
machine and laser control data. In order to build a common model of laser cooling, the details of the
cooling system were included in the formulation. For this reason, the model followed the assumption
that the electrical energy of the motor had a positive correlation with the current status of the laser
machine. The given data had the current power as the current status.
As discussed in the paper by Zhou et al. [3], the localization of the deep learning model is a
limitation even though its accuracy is high. The derivation of key features is needed to overcome the
localized learning. The authors of this paper proposed that the following features be used as the means
to train the neural network with the characteristics of laser welding: (1) mode is the feature to separate
each operation of the laser controller or status of the laser machine; (2) starting time and (3) starting
power: starting time and power are the factors influencing the current power by each initial status of
laser launch.
Generally, a time-series dataset is applied to neural networks because of the correlation between
current status and previous status. The energy profile of laser welding has fluctuations at the
microscopic level from 2500 Hz–5000 Hz. This feature easily causes an over-fitting problem of the
neural network. The collected data included the laser controller program. By using the controller
program, the fluctuation was reduced. The procedure of feature engineering consisted of smoothing,
peak detection, and synchronizing with the controller program.
The Savitzky–Golay filter was used as the smoothing filter. The filter windows were estimated by
the welding time derived from the controller program. For example, the window was set as 500 Hz
if each welding process time was 0.1 s. This number was the half of each peak profile. The welding
energy profile followed the fourth polynomial. This was determined by the analysis of experimental
data. Figure 9 shows raw data and the smoothing result.
Energies 2019, 12, 1799 11 of 20
Figure 9. Raw data (left) and the result (right) after the smoothing filter.
The filtered data were applied to the peak detection algorithm, which utilized continuous wavelet
transform-based pattern matching to find the highest point in each welding process. The algorithm
was described by Pan et al. [22] with detail. Low peaks were collected in {tli |i = 1...n}, while high
peaks were stored in {thi |i = 1...n}. The measured power data were represented in { p(t)|t = 1...n}.
Mode 1 used process time ∆t(= t − ti ) and starting power P(tli ) of each welding profile as the input of
neural network. Mode 2 divided each welding profile into welding status and cooling status. Welding
mode and cooling mode are defined in the following equation.
0 if t ≤ thi
mi ( t ) = (1)
1 if t > thi
The next section describes the proposed estimation methods of energy consumption and the
calculation of the energy efficiency of RLW.
Linear Regression
The first estimation model was linear regression using the modes of laser launching and cooling
defined in the feature engineering section. The energy profile consisted of a series of curves. Therefore,
the regression model was applied to find coefficients in order to generalize the experimental data.
After trying a linear model, exponential model, and power model, the highest accuracy model was
utilized as the estimation model. Each piece of the estimated profile was merged into a single profile
to show the overall energy consumption.
The learning rate is one of the significant factors in neural network optimization. In a pre-test,
1000 cycles, a 0.1 learning rate, and 0.2 momentum were effective for the convergence of the network
learning. Each model was compared with linear regression and compared with two different modes.
The root mean squared and R squared were calculated on each model and were the criteria for selection
of the optimal model to be used for energy efficiency.
In the efficiency calculation, comparison with alternative processes was needed. The welding
speed of each stitch was already fixed because of the effect on welding quality. The rapid movement
between stitches was the means to vary the process characteristics. The case study included variation
of rapid movement.
4. Case Study
The proposed methodology described in the previous section was applied to the RLW of a car rear
door assembly. This case study introduced the method for application of the integrated approach to a
real product and showed how the proposed approach can improve the accuracy. The estimation model
was implemented using RapidMiner 9.2. The process approach was applied to robot programming
written in COMAUPDL2. Assembly of the car rear door parts involved a combination of welding
positions in a set of 10 continuous stitches and 19 dimples made by a single laser welding robot with a
fixed work-in-process on a single fixture (Figure 11). The COMAU SmartLaser consisted of a robot arm
and laser machine. The robot arm had two parts. The first part had four revolute joints and the second
part another three laser manipulators. Figure 12 shows the experimental setting of the SmartLaser and
the rear door assembly on the fixture. Further details of the robot arm can be found in Erdős et al. [23].
Energies 2019, 12, 1799 13 of 20
This section analyzes the results of all experiments and compares these with each other.
The estimation results for energy consumption and the errors in the estimated energy profile were
investigated. First of all, the estimated profile was generated with a linear formulation by using a
regression model for fluctuation. A first order polynomial was used for drawing the profile. Each
cooling and launching cycle was divided into decrease direction and increase direction, respectively.
Linear regression was applied to find the coefficients of the exponential function on laser launch and
of the log function on the cooling operation.
The second approach is the mean to utilise the neural network model to estimate energy profile
with hyperbolic tangential activation function. The overall structure of single network is shown in
Figure 13. the variety of the neural network which are generated iteratively from one layer to three
layers with increasing of each layer from 4 nodes to 12 nodes. As mentioned in the feature engineering
section, laser mode, starting time and starting power are added into the training data by derivation
from the time-series power data of the training set. To compare the effectiveness of each aspect, the
following four conditions are investigated.
Energies 2019, 12, 1799 14 of 20
• (1) Linear formulation generated by the chain of linear regressions cooling and launching
• (2) Profile generated by neural network trained only using mode and time
• (3) Profile generated by neural network trained with additional data, which are Starting power
and starting time (mode 1)
• (4) Profile generated by neural network trained with the additional conditions on cooling and
launching (mode 2)
Condition (2) did not lead to convergence of the loss function. The comparison between (1) and (3)
is shown in Figures 14 and 15. The features of exponential and log functions caused the high peaks of
each stitch. In the comparison between Mode 1 and Mode 2, the 5 × 10 two-layered network showed
the highest accuracy in Mode 1, while the 10 × 15 × 5 network was better in Mode 2, as shown in
Figures 16 and 17 and Table 1. In the table, MAPE (Mean Absolute Percentage Error) showed the gap
between the estimated value and the measured value with the following equation.
N
1 | Emeasure (i ) − Estandard (i )|
MAPE( MeanAbsolutePercentageError ) =
N ∑ Estandard (i )
(2)
i =1
The detail of this function was described by Wang et al. [24] Both networks showed high
performance in terms of drawing an accurate profile; three layers showed a fine graph, then the
one of high RMS and R squared.
Category Mode 1 (5, 10) Mode 2 (5, 10) Mode 1 (10, 15, 5) Mode 2 (10, 15, 5) Measured
Energy (KWh) 5.355 5.393 5.384 5.393 5.407
Accuracy (%) 87.94 99.27 97.98 99.26 100.0
MAPE (%) 1.072 0.226 0.509 0.254 0.000
Energies 2019, 12, 1799 15 of 20
In terms of the energy efficiency, welding speed cannot be changed because it is optimized for the
material characteristics and shape of the parts being welded as a precondition. Rapid movement is
the only option for changing the process parameters. Figure 18 shows the estimated energy profile
generated by three rapid movements using the linear regression model. Scan speed was ten-times
faster than joint motion. If rapid movements were only done by scanning, then the total efficiency
could be increased. Rapid movement had a tremendous influence with low speeds, but for high
speeds, the variation was small. The optimal speed of rapid movement was related to cooling time to
get steady state energy. A longer time than this for movement was not necessary because it caused
Energies 2019, 12, 1799 16 of 20
meaningless energy loss. If the cooling time was long, it caused long steady state and increased energy,
while shorter time was better because cooling dropped exponentially.
Figure 16. Comparison with Mode 1 (left) and Mode 2 (right) in two-layer networks.
Figure 17. Comparison with Mode 1 (left) and Mode 2 (right) in three-layer networks.
Energies 2019, 12, 1799 17 of 20
Figure 18. (a) Relationship between welding speed and energy use; (b) energy profiles for varying
speeds. Comparison of the energy profile in different operating conditions.
5. Discussion
This section summarizes the major aspects of the proposed estimation model for RLW.
The proposed energy estimation model involved component selection by functional analysis and
machine learning compensation by using process parameters and experimental data. To conclude this
paper, we highlight those aspects of the proposed energy estimation model that can produce energy
profiles for a given welding program and a small amount of experimental data.
5.1. Provision of the Correlation with Laser Power and Process Time
This paper has formalized the effects of functional specification and process parameters in
an energy consumption model. This is based on the fact that the laser source determines energy
efficiency in the RLW processes. By providing a functional analysis and process analysis, this paper
shows how energy profiles were measured and the determination of key factors that influence total
energy consumption.
As shown, laser power and process time affected energy consumption in the RLW process.
Both analyses showed that the laser source and the chiller used more electrical energy than the robot
motors. Laser power is a significant criterion for estimating the energy profile. It is noted that low
power laser welding leads to low energy efficiency. However, door materials determine the laser
power and welding speed. Hence, rapid movements between welding stitches should be reduced.
The integrated approach showed the correlation between laser power and process time in a single
dimension. For this reason, the proposed approach is beneficial for optimization of path planning in
terms of process time and energy use.
The integration of both analyses is an interesting illustration. In some areas, for example in
“energy use distribution”, the laser machine had a higher rate of energy consumption than the rate
of the industrial robot. This was the opposite to the rate of metal cutting. The steady state energy
consumption was relatively lower than for machine tools. In the “reduction of process time”, the
shortest rapid movement time was not the best solution to minimize the total energy use. Very long
continuous welding caused an unstoppable increase in the temperature of the laser disks and, hence,
an increase in cooling.
Energies 2019, 12, 1799 18 of 20
6. Conclusions
This paper proposed an integrated approach merging a functional model and a process model with
machine learning-based compensation. The research was done in the domain of remote laser welding,
an assembly technique that has benefits for assembly in terms of energy savings. The characteristics
of remote laser welding were applied to the features of the machine learning model and showed
improved learning performance. This lead to optimization of performance for the process.
As described in the Discussion Section, the correlation between laser power and process power
had been formalized and the use of machine learning described. The formalization led to a better
understanding of the parameter settings in planning, as well as evaluation of the benefits of new
equipment. This is important in order to allow manufacturers to appreciate the benefits of new
equipment and processes, as well as quantifying energy savings. The machine learning aspects
showed an example of the application of artificial intelligence technology to overcome the lack of
formal understanding of the problem domain. The Discussion Section described the benefits and
shortcomings encountered in this application domain. The paper also described the differences
between the RLW application methods and those applied to metal cutting processes.
Future work is intended to adapt more features that have two cases: domain knowledge-based
feature (model driven) and data science-based feature (data driven). Treatment of data was done from
two perspectives. The first was the classical approach for data, which is called the data-driven method.
The second method used deep learning and is powerful also for use in other fields.
References
1. Um, J.; Stroud, I.A. Total energy estimation model for remote laser welding process In Proceedings of the
46th CIRP Conference on Manufacturing Systems, Sesimbra, Portugal, 29–31 May 2013.
2. Sebastian, T. Energy Efficiency in Manufacturing Systems; Springer: Berlin, Germany, 2012.
Energies 2019, 12, 1799 19 of 20
3. Zhou, L.; Li, J.; Li, F.; Meng, Q.; Li, J.; Xu, X. Energy consumption model and energy efficiency of machine
tools: A comprehensive literature review. J. Clean. Prod. 2016, 112, 3721–3734.
4. Pagone, E.; Salonitis, K.; Jolly, M. Energy and material efficiency metrics in foundries. In Proceedings of the
15th Global Conference on Sustainable Manufacturing, Haifa, Israel, 25–27 September 2017.
5. Gutowski, T.; Dahmus, J.; Thiriez, A. Electrical energy requirements for manufacturing processes.
In Proceedings of the 13th CIRP International Conference on Life Cycle Engineering, Leuven, Belgium,
31 May–2 June 2006.
6. Gutowski, T. The Carbon and Energy Intensity of Manufacturing. In Proceedings of the 40th CIRP
International Manufacturing Systems Seminar, Liverpool, UK, 30 May–2 June 2007.
7. Gutowski, T.; Branham, M.; Dahmus, J.; Jones, A.; Thiriez, A. Thermodynamic analysis of resources used in
manufacturing processes. Environ. Sci. Technol. 2009, 43, 1584–1590.
8. Gontarz, A.; Weiss, L.; Wegener, K. Evaluation approach with function-oriented modeling of machine tools.
In Proceedings of the 1st International Conference on Sustainable Intelligent Manufacturing, Leiria, Portugal,
29 June–1 July 2011.
9. Winter, M.; Herrmann, C. Eco-Efficiency of alternative and conventional cutting fluids in external cylindrical
grinding. In Proceedings of the 21st CIRP Conference on Life Cycle Engineering, Trondheim, Norway,
18–20 June 2014.
10. Winter, M.; Thiede, S.; Herrmann, C. Influence of the cutting fluid on process energy demand and surface
roughness in grinding—A technological, environmental and economic examination. Int. J. Adv. Manuf.
Technol. 2015, 77, 2005–2017.
11. Ma, F.; Zhang, H.; Hon, K.K.B.; Gong, Q. An optimization approach of selective laser sintering considering
energy consumption and material cost. J. Clean. Prod. 2018, 199, 529–537, doi:10.1016/j.jclepro.2018.07.185.
12. Zhu, Y.; Peng, T.; Jia, G.; Zhang, H.; Xu, S.; Yang, H. Electrical energy consumption and mechanical properties
of selective-laser-melting-produced 316L stainless steel samples using various processing parameters.
J. Clean. Prod. 2019, 208, 77–85, doi:10.1016/j.jclepro.2018.10.109.
13. Gregory, J.; Olivares, A.; Staffetti, E. Energy-optimal trajectory planning for robot manipulators with
holonomic constraints. Syst. Control Lett. 2012, 61, 279–291.
14. Avram, O.I.; Xirouchakis, P. Evaluating the use phase energy requirements of a machine tool system.
J. Clean. Prod. 2011, 19, 699–711, doi:10.1016/j.jclepro.2010.10.010.
15. Duflou, J.R.; Sutherland, J.W.; Dornfeld, D.; Herrmann, C.; Jeswiet, J.; Kara, S.; Hauschild, M.; Kellens, K.
Towards energy and resource efficient manufacturing: A processes and systems approach. CIRP Ann. 2012,
61, 587–609, doi:10.1016/j.cirp.2012.05.002.
16. Newman, S.T.; Nassehi, A.; Imani-Asrai, R.; Dhokia, V. Energy efficient process planning for CNC machining.
CIRP J. Manuf. Sci. Technol. 2012, 5, 127–136, doi:10.1016/j.cirpj.2012.03.007
17. Marvuglia, A.; Messineo, A. Using recurrent artificial neural networks to forecast household electricity
consumption. In Proceedings of the 2011 2nd International Conference on Advances in Energy Engineering,
Bangkok, Thailand, 27–28 December 2011.
18. Biswas, M.R.; Robinson, M.D.; Fumo, N. Prediction of residential building energy consumption: A neural
network approach. Energy 2016, 117, 84–92, doi:10.1016/j.energy.2016.10.066.
19. Um, J., Gontarz, A., and Stround, I., Developing energy estimation model based on Sustainability KPI of
machine tools. Procedia CIRP 2015, 26, 217–222, doi:10.1016/j.procir.2015.03.002.
20. Giesen, A.; Speiser, J. Fifteen years of work on thin-disk lasers: Results and scaling laws. IEEE J. Quantum
Electron. 2007, 13, 598–609, doi:10.1109/JSTQE.2007.897180.
21. Reinhart, G.; Munzert, U.; Vogl, W. A programming system for robot-based remote-laser-welding with
conventional optics. CIRP Ann. Manuf. Technol. 2008, 57, 37–40, doi:10.1109/JSTQE.2007.897180.
22. Du, P.K.; Warren, A.; Lin, S.M. Improved peak detection in mass spectrum by incorporating continuous
wavelet transform-based pattern matching. Bioinformatics 2006, 22, 2059–2065, doi:10.1093/bioinformatics/
btl355.
23. Erdős, G.; Kardos, C.; Kemény, Z.; Kovács, A.; Váncza, J. Process planning and offline programming for
robotic remote laser welding systems. Int. J. Comput. Integr. Manuf. 2016, 29, 1287–1306, doi:10.1080/
0951192X.2015.1033753.
24. Wang, S.; Liang, Y.C.; Li, W.D.; Cai, X.T. Big Data enabled Intelligent Immune System for energy efficient
manufacturing management. J. Clean. Prod. 2018, 195, 507–520.
Energies 2019, 12, 1799 20 of 20
c 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access
article distributed under the terms and conditions of the Creative Commons Attribution
(CC BY) license (http://creativecommons.org/licenses/by/4.0/).