A Generic WebLab Control Tuning Experience Using the Ball and Beam Process and Multiobjective Optimization Approach

: In control engineering education, the possibility of using a real control system in the learning process motivates professors to improve both students’ knowledge and skills, thus avoiding an approach only based on control theory. While considering that control engineering laboratories are expensive, mainly because educational plants should reproduce classical problems that are found in the industry, the use of virtual laboratories appears as an interesting strategy for reducing costs and improving the diversity of experiments. In this research, remote experimentation was assumed regarding the ball and beam process as an alternative didactic methodology. While assuming a nonlinear and unstable open-loop process, this study presents how students should proceed to control the plant focusing on the topic that is associated with multiobjective optimization. Proportional-Integral-Derivative (PID) controller was tuned considering the Non-dominated Sorting Genetic Algorithm (NSGA-II) to illustrate the WebLab learning procedures described in this research. The proposed strategy was compared to the Åström’s robust loop shaping method to emphasize the performance of the multiobjective optimization technique. Analyzing the feedback provided by the students, remote experimentation can be seen as an interesting approach for the future of engineering learning, once it can be directly associated with industry demand of connected machines and real-time information analysis.


Introduction
The control engineering is an area of study where scientific and technological development supports the quality improvement of process and products, thus guarantying repeatability and precision in industrial applications. In this context, automatic control can be considered to be essential in many fields of science and engineering [1], e.g., electronics, chemistry, robotics, manufacturing systems, etc. The two main difficulties that are presented in [2] are related to the educational process that is associated with automatic control, the first one consists in understanding how to draft control systems, and the second is how to associate mathematical results to practical problems.
Active learning methodologies are being used in many areas, and digital media associated with learning activities seem to increase both students' knowledge and skills [3]. The main idea is to provide 1. 1

.2. Related Works Associated with the Construction of the Ball and Beam Plant
The works that are presented in this subsection constitute the second part of the literature review.
Here, works focusing on the construction of a physical plant for the ball and beam process are described.
During simulations, disturbances do not affect the behavior of the system. On the other hand, when a real plant is considered, it is possible to evaluate problems that are associated with the network during communication, additional dynamics as transport delay, noisy signals, among other effects that are associated with the control of real systems. One example is the work that is presented in [19], where an apparatus was used during tests of both static and dynamic Sliding Mode Control (SMC) schemes. Additionally, four approaches assuming two modeling strategies, one simplified and other complete, were considered. The results showed that when the complete model was assumed, improved performance for controlling the ball position was obtained, mainly when the dynamic controller was considered.
An experimental approach was performed in [20], where the developed ball and beam apparatus was associated with a Linear Quadratic Regulator (LQR) controller. The main idea was to both stabilize and control the system by changing poles to optimal locations. While assuming a genetic algorithm to determine the optimal LQR configuration, two approaches were considered during the tests. The first one assumed a model-based configuration, while the second one considered a non-model-based controller. Both of the approaches used the PID controller and its variations, and the LQR strategy. The results showed that the LQR model-based approach provided better results when the steady-state error was evaluated.
Another ball and beam apparatus was presented in [21], where a Takashi-Sugeno fuzzy model-based adaptive dynamic surface controller was proposed. The proposed controller was compared to the conventional Dynamic Surface Control (DSC) strategy. The results showed that superior performance was obtained by the fuzzy approach in terms of computation robustness and feasibility. Other work considering a fuzzy controller was addressed in [22], but, in this case, real-time tuning was the objective of the study. Assuming the Differential Evolution (DE) algorithm [23], the possibility of real-time tuning while considering a classical PID structure was presented.
The association between the ball and beam plant and WebLabs can be found in [13]. This work was developed with the purpose of helping students of control engineering during their activities associated with advanced control theory. In this sense, students could compare the performance of their controllers with classical and well-established strategies. The WebLab system was developed, associating Easy Java Simulations (EJS) and Moodle learning management systems.
A state feedback LQR based control scheme, which consists of Linear Quadratic Estimator (LQE) and LQR, was presented in [24], to conclude this literature review section of this article. The ball position was considered the output, where all system states were calculated using both the ball position and its time derivatives. Additionally, the ball and beam apparatus was built while considering a servo motor used to modify the beam angle and a potentiometer to measure the ball position.
The literature review that is presented in the first part of this article showed that distinct control strategies can be associated with the ball and beam control problem. However, just a few works considered virtual experimentation as the one assumed in the present study. The research question associated with this work can be stated based on this statement. Is it possible to propose an interesting approach directed to advanced control learning, which is capable of reaching as many students as possible that are focused on learning through practical experiments, even considering a limited laboratory structure?
Previous works showed that the ball and beam process became a benchmark for testing different control strategies, showing to be an important process in control engineering learning. By associating virtual laboratories and control engineering, this article presents a WebLab approach for advanced control learning to make the most of the didactical processes available at the university, in order to reach as many students as possible. We considered the ball and beam process as a case study to illustrate all the procedures that should be performed by the students. The methodology consists of the development of advanced control strategies for the process using multiobjective optimization. While considering a wide range of algorithms available for this purpose, we have described all of the procedures using the Non-dominated Sorting Genetic Algorithm (NSGA-II). Finally, students that participate in the control WebLab experience could deal with distinct software, such as LabView and MATLAB, Python programming, among others, verifying the importance of system integration for remote data acquisition and control.
The next section of this article presents the ball and beam apparatus adopted in this study focusing on process modeling procedures. Section 3 presents the theoretical basis of PID control. Section 4 addresses the procedures to apply multiobjective optimization in control problems. When the classical PID control action is assumed, then the procedures to use evolutionary multiobjective optimization (EMO) algorithms to tune the controller are based on three steps: the definition of the multiobjective problem, the multiobjective optimization routine, and the multicriteria decision making. All of the previously mentioned steps are addressed in Section 4. Section 5 presents the experimental routine considering the ball and beam apparatus. Section 6 describes the Non-dominated Sorting Genetic Algorithm (NSGA-II), which is the algorithm that is assumed in this study. In the same section, the results that were obtained using NSGA-II were compared with Åström's robust loop shaping method to show the advantages of the multiobjective approach. Finally, Section 7 presents the conclusions and future works that are associated with this research.

Ball and Beam Process Description and Modelling
This section presents the ball and beam apparatus used during the development of this research, as well as two modeling approaches that were assumed to obtain the process model.

The Ball and Beam Apparatus
While considering that different hardware implementations may cause differences in the dynamic of the plant, which are mainly due to variations on both mechanical and electrical architectures, the coefficients of the transfer function that are associated with each plant could change [1]. In this way, this section describes the plant that is adopted in this research. The ball and beam apparatus corresponds to the CE 106 from Tecquipment [25] (see Figure 1), which was developed to be an educational plant in learning activities that is associated with control engineering.
The process that was adopted in this study is based on a beam that can rotate through an angle of ±10 • . The input signal represents the motor input voltage (0 to +10 VDC), and two outputs are available (0 to ±10 VDC): the ball position and the beam angle, respectively. When a signal is applied to the motor input, it varies the speed of the motor. A profiled cam, which is attached to the motor shaft, links the main beam via a secondary one. Additionally, a potentiometer allows for measuring the voltage that is associated with the beam angle. Two parallel wires are stretched along the top of the beam, where the ball may roll, to measure the ball position. When the ball rests between the wires, it allows for a fraction of the source voltage to be measured at the other wire while using a potentiometer.

Ball and Beam Process Description and Modelling
This section presents the ball and beam apparatus used during the development of this research, as well as two modeling approaches that were assumed to obtain the process model.

The Ball and Beam Apparatus
While considering that different hardware implementations may cause differences in the dynamic of the plant, which are mainly due to variations on both mechanical and electrical architectures, the coefficients of the transfer function that are associated with each plant could change [1]. In this way, this section describes the plant that is adopted in this research. The ball and beam apparatus corresponds to the CE 106 from Tecquipment [25] (see Figure 1), which was developed to be an educational plant in learning activities that is associated with control engineering.
The process that was adopted in this study is based on a beam that can rotate through an angle of 10°. The input signal represents the motor input voltage (0 to 10 VDC), and two outputs are available (0 to 10 VDC): the ball position and the beam angle, respectively. When a signal is applied to the motor input, it varies the speed of the motor. A profiled cam, which is attached to the motor shaft, links the main beam via a secondary one. Additionally, a potentiometer allows for measuring the voltage that is associated with the beam angle. Two parallel wires are stretched along the top of the beam, where the ball may roll, to measure the ball position. When the ball rests between the wires, it allows for a fraction of the source voltage to be measured at the other wire while using a potentiometer. There is no fixed point in the beam extremity, just in the center. It is possible to consider two modeling approaches, one that is based on a simplified formulation and a second one based on a complete formulation. Both of the models are presented in the sequence of this work.

Physical Modeling
This subsection describes the two modeling methods of the ball and beam process to obtain the simplified model or the full model, which can be used as parametric models for system simulation.

Simplified Model
The modeling process starts with the dynamic of the ball when rolling on the beam. While assuming the force, which accelerates the ball as it rolls, the component can be formulated based on the gravitational force, see Figure 2. There is no fixed point in the beam extremity, just in the center. It is possible to consider two modeling approaches, one that is based on a simplified formulation and a second one based on a complete formulation. Both of the models are presented in the sequence of this work.

Physical Modeling
This subsection describes the two modeling methods of the ball and beam process to obtain the simplified model or the full model, which can be used as parametric models for system simulation.

Simplified Model
The modeling process starts with the dynamic of the ball when rolling on the beam. While assuming the force, which accelerates the ball as it rolls, the component can be formulated based on the gravitational force, see Figure 2.  The simplified model considers no friction during the ball sliding over the beam. In this way, the mathematical equation that describes the movement is . (1) When small angles are considered, . Finally, the transfer function of the plant relates the beam angle, , and the ball position, , as presented in Equation (2). .

Full Model
The full model can be obtained while considering the Lagrangian method, where the Lagrangian, , for the process can be stated as * , ( In Equation (3), * represents the kinetic energy and is the potential energy in the system. When an open-loop equation is considered, the potential energy becomes zero, since there are no springs in the system. The potential energy that is provided by the gravity is considered and the Lagrangian is given by the kinetic energy, which represents both the translational and the rotational kinetic energy of the ball. * 1 2 where is the ball mass and is the ball moment of inertia. The translational velocity of the ball, , and the angular velocity of the ball, , have to be re-written in terms of the generalized coordinates and to obtain the transfer function. While considering an angular distance , which was rolled by the ball, the angular velocity can be obtained, see Figure 3. In this case, represents the rolling radius of the ball.
.   The simplified model considers no friction during the ball sliding over the beam. In this way, the mathematical equation that describes the movement is x. (1) When small angles are considered, sin θ = θ. Finally, the transfer function of the plant relates the beam angle, θ, and the ball position, x, as presented in Equation (2).

Full Model
The full model can be obtained while considering the Lagrangian method, where the Lagrangian, L, for the process can be stated as In Equation (3), U * represents the kinetic energy and T is the potential energy in the system. When an open-loop equation is considered, the potential energy becomes zero, since there are no springs in the system. The potential energy that is provided by the gravity is considered and the Lagrangian is given by the kinetic energy, which represents both the translational and the rotational kinetic energy of the ball.
where m is the ball mass and I is the ball moment of inertia. The translational velocity of the ball, v, and the angular velocity of the ball, ω, have to be re-written in terms of the generalized coordinates x and θ to obtain the transfer function. While considering an angular distance x, which was rolled by the ball, the angular velocity ω can be obtained, see Figure 3. In this case, r represents the rolling radius of the ball.
Information 2020, 11, 132 6 of 22 The simplified model considers no friction during the ball sliding over the beam. In this way, the mathematical equation that describes the movement is . (1) When small angles are considered, . Finally, the transfer function of the plant relates the beam angle, , and the ball position, , as presented in Equation (2). .

Full Model
The full model can be obtained while considering the Lagrangian method, where the Lagrangian, , for the process can be stated as * , ( In Equation (3), * represents the kinetic energy and is the potential energy in the system. When an open-loop equation is considered, the potential energy becomes zero, since there are no springs in the system. The potential energy that is provided by the gravity is considered and the Lagrangian is given by the kinetic energy, which represents both the translational and the rotational kinetic energy of the ball. * 1 2 where is the ball mass and is the ball moment of inertia. The translational velocity of the ball, , and the angular velocity of the ball, , have to be re-written in terms of the generalized coordinates and to obtain the transfer function. While considering an angular distance , which was rolled by the ball, the angular velocity can be obtained, see Figure 3. In this case, represents the rolling radius of the ball.
.     While considering the angular velocity that is presented in Equation (5), and the resultant velocity of Equation (6), the Lagrangian of Equation (4) can be re-written, as presented in Equation (7), while the moment of inertia of the ball is presented in Equation (8).
where represents the ball radius. While considering that the ball rest between two wires, the radius is always larger than the rolling radius , see Figure 5. .
The complete model of the ball and beam process, which associates the radius , the rolling radius, , the beam angle, , the ball position, , and the gravitational acceleration, , can be expressed by: .
(10) Equation (10) represents a nonlinear model, as it involves a nonlinear term that associates both the beam velocity and its angle, and the ball position. In this case, the model must be linearized while considering small-angle variations to obtain an approximated linear system. Moreover, the velocity of the beam and its acceleration are also both considered to be small. In this way, the transfer function of the ball and beam process can be obtained while considering the full formulation associating the ball position (output of the system), and the beam angle (input of the system), see Equation (11). While considering the angular velocity that is presented in Equation (5), and the resultant velocity of Equation (6), the Lagrangian of Equation (4) can be re-written, as presented in Equation (7), while the moment of inertia of the ball is presented in Equation (8).
where R represents the ball radius. While considering that the ball rest between two wires, the radius R is always larger than the rolling radius r, see Figure 5. While considering the angular velocity that is presented in Equation (5), and the resultant velocity of Equation (6), the Lagrangian of Equation (4) can be re-written, as presented in Equation (7), while the moment of inertia of the ball is presented in Equation (8).
where represents the ball radius. While considering that the ball rest between two wires, the radius is always larger than the rolling radius , see Figure 5. .
The complete model of the ball and beam process, which associates the radius , the rolling radius, , the beam angle, , the ball position, , and the gravitational acceleration, , can be expressed by: . (10) Equation (10) represents a nonlinear model, as it involves a nonlinear term that associates both the beam velocity and its angle, and the ball position. In this case, the model must be linearized while considering small-angle variations to obtain an approximated linear system. Moreover, the velocity of the beam and its acceleration are also both considered to be small. In this way, the transfer function of the ball and beam process can be obtained while considering the full formulation associating the ball position (output of the system), and the beam angle (input of the system), see Equation (11). Finally, the model of the system can be obtained based on the Lagrangian equation.
The complete model of the ball and beam process, which associates the radius R, the rolling radius, r, the beam angle, θ, the ball position, x, and the gravitational acceleration, g, can be expressed by: Equation (10) represents a nonlinear model, as it involves a nonlinear term that associates both the beam velocity and its angle, and the ball position. In this case, the model must be linearized while considering small-angle variations to obtain an approximated linear system. Moreover, the velocity of the beam and its acceleration are also both considered to be small. In this way, the transfer function of Information 2020, 11, 132 8 of 23 the ball and beam process can be obtained while considering the full formulation associating the ball position (output of the system), and the beam angle (input of the system), see Equation (11).

PID Control and Performance Measures
Proportional-Integral-Derivative (PID) controllers are used in 90% of control applications [26], as reported by the Japan Electric Measuring Instrument Manufacturers' Association. The main motivation to be considered is that other approaches often need mathematical models as a state-space representation or transfer functions to be formalized, short implementation time, and when the implementation effort is evaluated, other techniques frequently present small advantages over the PID controllers [27].
The classical PID controller transfer function, C(s), can be defined based on three gains: the proportional gain, K p , the derivative gain, K d , which can be associated with the derivative time, T d , and the integral gain, K i , which can be associated with the integral time, T i .
Frequently, one of the most common objectives in control systems is to reduce the errors between the reference and the system output in a shorter time [18,24]. In the case of the ball and beam system, the objective is to control the ball position that is based on an error signal, e(t), which represents the difference between the desired ball position and the real ball position. A classical PID work on the error signal, which is processed to generate a control signal. These procedures form an iterative process until the stabilization of the ball on the desired position.
In control systems, the analysis of the response of the system over time consists of two sections: the transient response and the steady-state response. The first one is related to a change in the reference of the system, the response goes from the initial to a final state. Performance metrics can be assumed to evaluate the transient response of the system, such as the rise time (Rt); the time spent by the system to reduce the output error from 90% to 10%; the peak amplitude or overshoot (Pa), which is the maximum value reached by the system response caused by a reference change; and, the settling time (St), representing the time that is taken for the error to be reduced to a value smaller than 2% or 5%. The second one, or the steady-state response, refers to the output system behavior is stabilized, i.e., when the system do not present an oscillatory behavior or the time approaches to infinity, in this situation the stationary error (Se), which indicates the percentage error after output stabilization, can be evaluated [1]. Figure 6 illustrates the previously mentioned performance coefficients for a second-order system step response, where a change in the reference is performed and the system presents an oscillatory behavior until the stabilization, more information regarding the previously mentioned performance measures can be found in [1].
i.e., when the system do not present an oscillatory behavior or the time approaches to infinity, in this situation the stationary error ( ), which indicates the percentage error after output stabilization, can be evaluated [1]. Figure 6 illustrates the previously mentioned performance coefficients for a secondorder system step response, where a change in the reference is performed and the system presents an oscillatory behavior until the stabilization, more information regarding the previously mentioned performance measures can be found in [1].  The PID control system has block diagram variations that can be found in the literature [28]. Some of them are the PI-D ( Figure 7a) and Integral-Proportional Derivative (I-PD) (Figure 7b). The PI-D variation has the derivative portion removed from the forward branch and placed on the feedback branch. The forward branch corresponds to the direct way between the input and the output of the system. The feedback branch corresponds to a measure of the output signal that is used to calculate the error between the reference and the real system's output. This modification results in a control strategy where the derivative term acts on the plant output signal y(t), not in the difference x(t) − y(t). In this way, no differentiation is applied over x(t) step input, because step differentiation would generate a high amplitude signal at the actuator, thus leading to sorter lifetime or damages. In other cases, PI-D is used to avoid the amplification of high-frequency noise that might be present in the reference or in the error signals.
Information 2020, 11, 132 9 of 22 PI-D variation has the derivative portion removed from the forward branch and placed on the feedback branch. The forward branch corresponds to the direct way between the input and the output of the system. The feedback branch corresponds to a measure of the output signal that is used to calculate the error between the reference and the real system's output. This modification results in a control strategy where the derivative term acts on the plant output signal , not in the difference . In this way, no differentiation is applied over step input, because step differentiation would generate a high amplitude signal at the actuator, thus leading to sorter lifetime or damages. In other cases, PI-D is used to avoid the amplification of high-frequency noise that might be present in the reference or in the error signals.   (16), where represents the control signal and is the current experiment time in seconds [29]. Some metrics that are commonly used in experiments that are related to the error are the Integral Absolute Error (IAE) (13) and the Integral of Time multiplied by Absolute of the Error (ITAE) (14) to evaluate the controller's performance. Related to the control effort, it can be cited the Integral of Time multiplied by the Squared Error (ITSE) (15) and the Integral of the Absolute Variation of Control Signal (IAVU) (16), where u represents the control signal and t is the current experiment time in seconds [29].

Multiobjective Optimization Applied to Control Engineering
When trying to solve an engineering problem, it is common to search for solutions that provide the best performance when compared to other possible solutions. This is the principle of an optimization problem, where a solution based on predetermined assessment conditions, which improves the system performance, is required [30].
Most of the engineering problems can be stated as optimization problems. When multiple objectives are considered, conflicting objectives often appear during the problem formulation. An example of this case was addressed in [31], where the authors focused on the improvement of the system performance, trying to avoid the increase of energy consumption.
Assuming the condition of conflicting objectives, some strategies can be stated in order to obtain a reasonable solution for the problem. Most of these strategies consist of aggregating multiple objectives in a single function, thus reducing the complexity of analyzing a set of solutions. Aggregating strategies have an articulation of preferences prior to the optimization process, e.g., by defining weighting factors among all of the objectives, where the solution will reflect the preferences of the decision-maker. Examples of the previously mentioned case are the weighted global criterion method, the weighted sum method, the lexicographic method, the weighted min-max method, and others (for a detailed description see [32]).
Some recent techniques that can be used to solve this kind of problems are the Jaya algorithm [33], the Proactive Particles in Swarm Optimization (PPSO) [34], the Hybrid Population-Based Mean-Variance Mapping Optimization (MVMO-PH) [35], the Successful-History based Adaptative Differential Evolution with Linear size reduction and four adaptation strategies with Improved ε constraint-handling method (LSHADE44-IEpsilon) [36], the Hunting Optimization (HO) [37], and the Modified Adaptative Bats Sonar Algorithm with Doppler Effect (MABSA-DE) [38].
Another approach for solving multiobjective optimization problems is though the simultaneous optimization of different objective functions. In this case, each objective function is treated as a distinct problem, and a set of solutions is obtained at the end of the optimization process. In this case, it is not necessary to normalize the cost functions, as the analysis and selection of a solution will be made a posteriori of the optimization process. Some of the newer techniques that can be used are the Bayesian Multi-Objective Optimization (BMOO) [39], the Adaptative Gradient Multiobjective Particle Swarm Optimization (AGMOPSO) [40], the Multiobjective Simulated Kalman Filter algorithm (MOSKF) [41], the Hybrid Multi-Objective Cuckoo Search (HMOCS) [42], among others.
The main difference between the two presented approaches is that the second one needs a multicriteria decision making (MCDM) strategy. In classical decision-making researches, it is common to find some formulations to aggregate multiple cost functions into a single objective function, where the focus is to select an optimal solution. The MCDM method evaluates the trade-offs of the candidate solutions to assess the contribution of each one for the problem, ranking them [43].
While considering MOPs in the control engineering area, the process of improving performance measures, those that are mentioned at the end of the previous section, are directly related to the selected control action. When the classical PID control action is assumed, the procedures to use evolutionary multiobjective optimization (EMO) to tune this controller were provided as a general framework in [44]. These procedures consist of three main stages that are illustrated in Figure 8. [43].
While considering MOPs in the control engineering area, the process of improving performance measures, those that are mentioned at the end of the previous section, are directly related to the selected control action. When the classical PID control action is assumed, the procedures to use evolutionary multiobjective optimization (EMO) to tune this controller were provided as a general framework in [44]. These procedures consist of three main stages that are illustrated in Figure 8. In a two objective optimization problem, the Pareto Front can be represented in a two-dimensional (2D) cartesian plot, where each cost function represents a specific axis. In this case, the Pareto Front, as presented in Figure 2, is represented by the solutions that are plotted in green color. The first stage is the problem statement, where both the parametric model and the cost functions are defined. The second stage consists of the multiobjective optimization, where an algorithm (optimization strategy) is assumed to obtain the Pareto Front. The Pareto Front consists of a set of non-dominated solutions when the cost functions are analyzed according to the decision variables. In a two objective optimization problem, the Pareto Front can be represented in a two-dimensional (2D) cartesian plot, where each cost function represents a specific axis. In this case, the Pareto Front, as presented in Figure 2, is represented by the solutions that are plotted in green color.
The relation between dominated and non-dominated solutions consists of a comparison of the cost functions for each candidate solution. If there is no solution that is better in all the evaluated cost functions, that candidate solution is non-dominated, otherwise, it is a dominated solution [45]. The set of all the non-dominated solutions will form the Pareto Front; Figure 9 presents an example of dominated and non-dominated solutions obtained during the optimization process. The relation between dominated and non-dominated solutions consists of a comparison of the cost functions for each candidate solution. If there is no solution that is better in all the evaluated cost functions, that candidate solution is non-dominated, otherwise, it is a dominated solution [45]. The set of all the non-dominated solutions will form the Pareto Front; Figure 9 presents an example of dominated and non-dominated solutions obtained during the optimization process. Finally, the last stage assumes a multicriteria decision-making method to select a feasible solution for the problem.
While considering the wide range of algorithms that are available to solve multiobjective optimization problems [46], the EMO methods are interesting for consideration. The use of EMO algorithms to obtain the Pareto Front is motivated by the idea that these techniques allow for obtaining multiple Pareto optimal sets in a single execution, thus reducing the complexity of the problems or when compared with direct methods [47]. These methods are based on the Pareto Front approximation, where the Pareto optimal solutions represent the result of the optimization process Finally, the last stage assumes a multicriteria decision-making method to select a feasible solution for the problem.
While considering the wide range of algorithms that are available to solve multiobjective optimization problems [46], the EMO methods are interesting for consideration. The use of EMO algorithms to obtain the Pareto Front is motivated by the idea that these techniques allow for obtaining multiple Pareto optimal sets in a single execution, thus reducing the complexity of the problems or when compared with direct methods [47]. These methods are based on the Pareto Front approximation, where the Pareto optimal solutions represent the result of the optimization process [48]. Through the Pareto Front approximation, it is possible to analyze and select a suitable solution that is based on distinct criteria. The procedure to select the candidate solution can be performed through an empirical analysis or considering well-founded decision-making strategies [49], some of these strategies are discussed in Section 4.3 of this article.
When applied to control problems, the use of MCDM strategies, those that can relate multidimensional data, is indicated [48]. In this case, it is possible to emphasize the trade-offs analysis through data visualization. The sequence of this section describes the multiobjective optimization procedures (MOOP), showing, in detail, the three stages, as presented in Figure 8.

Multiobjective Problem Statement
The Equation (17) shows a general MOP formulation [32], where the goal is to simultaneously minimize n objectives, J i (ϕ), i = [1, 2, . . . , n], of a variable vector ϕ (1 × m) with l inequality constraints and q equality constraints, see Equation (18). In a control problem, the coefficients of the minimization function [x 1 , x 2 , . . . , x m ] represent the gains of the controller. where The first step is associated with the definition of a parametric model that will be used to perform computational simulations to calculate the cost functions. In control applications, a transfer function is assumed as a parametric model. The transfer function can be obtained through a system identification procedure [50], or through physical modeling. In Section 2, two models were described and assumed for the ball and beam apparatus while considering the physical modeling approach.
The second step is the cost functions definition. In this case, at least one error criterion, this being related to the performance of the system, and one robustness criterion, which is associated with controllers' confidence, are assumed. Frequently, the performance criterion is related to the minimization of the error signal, improving the system capability to follow the reference and reject disturbances in the system, as presented in Equations (13) and (14). On the other hand, the robustness criterion considers the behavior of the control signal that is subject to noises and unexpected variations in the working cycles, see Equations (15) and (16).

Multiobjective Optimization Process
The multiobjective optimization stage is dedicated to obtaining a set of possible solutions for the control problem, which is represented by the parametric model and associated with the minimization of the cost functions. In the case of non-conflicting objectives, a single solution can be considered to be an optimal solution, and the use of additional methods to analyze the set of solutions is not necessary. However, when conflicting objectives are considered, the final result of the optimization process is a set of solutions that are defined as non-dominated solutions, named the Pareto Front [45]. In this case, a ranking criterion can be assumed to define the optimal solution.
The multiobjective optimization techniques are strategies that can be used to approximate the Pareto Front. These techniques are classified into two basic types: the direct methods, and multiobjective evolutionary algorithms [51].
The direct methods are a class of algorithms that are oriented by the iterations of previous results according to the descending information, e.g., gradient information, presenting fast convergence and searching efficiency. Nevertheless, this kind of algorithm has some problems when the gradient is hard to, or cannot be, calculated, frequently getting trapped in local optimal solutions [52]. Some strategies for approximating the Pareto Front through direct methods are the Normal-Boundary Intersection (NBI) [53], the Physical Programming [54], the Normalized Normal Constraint method [55], among others.
In another way, the multiobjective evolutionary algorithms are stochastic optimization techniques, where, at each iteration, the cost functions are evaluated for each individual, obtaining the dominance relations among the solutions. This information is used to select potentially better solutions and produce the next generation [56]. Some metaheuristics that can be used for this purpose are the Nondominated Sorting Genetic Algorithm II (NSGA-II) [57], the Multi-Objective Artificial Bee Colony Algorithm with Regulation Operators (RMOABC) [58], and the Adaptive Multiswarm Comprehensive Learning Particle Swarm Optimization (AMCLPSO) [59], among others.

Multicriteria Decision Making
The multicriteria decision-making stage is the last step of the multiobjective optimization design procedure, where a method is applied to select a reasonable solution for the problem that is based on the conflicting criteria of the objective functions.
There are two main approaches for applying MCDM techniques, one considering a single method and others assuming multiple multicriteria methods. In the first approach, the decision-maker should have knowledge regarding the MOP to decide which objective function will be emphasized. In the second approach, while considering multiple MCDM methods, the main idea is to adequately visualize the decision scene defining a mathematical model to reflect the decision-maker preferences. In this way, it is possible to reduce the set of solutions, and to then solve the mathematical model while using another MCDM technique, i.e., the set of solutions is successively reduced [49].
Some interesting examples of MCDMs that can be directly applied to select a candidate solution in the Pareto Front are the cooperative solutions of the bargaining games [60]. These can be divided into the Nash solution (Ns), the Kalai-Smorodinsky solution (K-Ss), and the Egalitarian solution (Es). Both cooperative solutions need the specification of a utopian solution and a disagreement point to be founded, and these points consist of the best (even impossible) situation and the combination of the worst value for each cost function for all of the candidate solutions, respectively. The Ns is a solution selected from the set of candidates where the biggest volume (n-dimensional, depending on the number of cost functions) is obtained, starting from the disagreement point until the candidate solution evaluated. The K-Ss is the candidate solution in the point of intersection between a line connecting the points of disagreement and the utopian solution and the Pareto Front, and the Es is obtained through the intersection of a line, which started in the disagreement point and had an angle of 45 degrees, and the Pareto Front.

Remote Experiment Description
This section describes the remote experiment structure, followed by procedures that have to be performed by the students along an advanced control course.

Experiment Structure
The control system of the remote ball and beam experiment was hosted on a Raspberry Pi 3 microcomputer, which was connected to the university network. This microcomputer, called a lab server, is intended to receive commands from a user while considering an interface hosted on a cloud server. This cloud server communicates with the lab server through an MQTT protocol-based architecture.
The user must connect to a service broker to access the remote experiment. After authorized access, the date and time scheduling routine can be configured. After that, the user can perform a remote experiment. Once authorized to run the experiment, the user is redirected to the interface webpage of the experiment, where he will perform the actions and collect datasets of his results that will be stored on a storage server for later access if needed. The architecture that is described here is called ISA-iLab Shared Architecture [61].
Additionally, the operating system that is used for both the lab server and the interfaces located on the cloud server consists of Linux distributions. For programming, the control systems of the remote experiment, as well as communication via MQTT, the Python programming language was used. Additionally, for the programming of the controller of the remote experiment, the C language was considered. Figure 10 illustrates all of the architecture described above, which is composed of the cloud server, the lab server, the communication protocols, and the intermediation of this communication through ISA. This architecture was called ELSA-SP.
First, the user accesses the ISA server through a university-provided internet address by entering login and password data. In the sequence, the user selects the desired date and time to access the experiment. Once the date and time were configured, the user requests the access and the execution of the experiment. ISA will send an access coupon containing data, such as user ID, server ID, and passkey to the lab server. The lab server checks on a new ISA connection whether the coupon is correct and, if it is, ISA provides a ticket to the lab server containing access data and the duration of the experiment. The lab server authorizes the execution of the experiment by informing ISA of the ticket data. ISA, which has ticket data, finally authorizes the user to access and operate the experiment interface.
The main contribution that is expected through experimentation procedures in control engineering courses is to assist in the critical review in studies of advanced control. In this way, the ball and beam system can be used in the comparison between the configuration and analyses of the controller's parameters. Additionally, the controller's performance obtained while considering two methods can also be verified. In this case, the first approach usually assumes a classical controller tuning rules, and the second one is based on multiobjective optimization with evolutionary algorithms. methods can also be verified. In this case, the first approach usually assumes a classical controller tuning rules, and the second one is based on multiobjective optimization with evolutionary algorithms.

Experimental Procedures
This section presents the experimental procedures that are associated with the control of the ball and beam system. A sampling time of 100 ms was considered, and the additional parameters are the ball radius, 2.0 cm, its mass, 50.0 g; the rolling radius was measured as 1.65 cm, and the beam length of 1 m. For the remote experimental procedures, we suggest two distinct control structures that are based on PD control action. The first one considers a cascade controller with an integrator as a highfrequency filter. Additionally, the beam angle control was performed, as shown in Figure 11. The second one with the derivative part in the feedback loop, as shown in Figure 12. The second structure

Experimental Procedures
This section presents the experimental procedures that are associated with the control of the ball and beam system. A sampling time of 100 ms was considered, and the additional parameters are the ball radius, R = 2.0 cm, its mass, m = 50.0 g; the rolling radius was measured as r = 1.65 cm, and the beam length of 1 m.
For the remote experimental procedures, we suggest two distinct control structures that are based on PD control action. The first one considers a cascade controller with an integrator as a high-frequency filter. Additionally, the beam angle control was performed, as shown in Figure 11. The second one with the derivative part in the feedback loop, as shown in Figure 12. The second structure with the derivative term in the feedback loop is an interesting approach, as students can analyze and understand how to reduce the noise amplification.

Experimental Procedures
This section presents the experimental procedures that are associated with the control of the ball and beam system. A sampling time of 100 ms was considered, and the additional parameters are the ball radius, 2.0 cm, its mass, 50.0 g; the rolling radius was measured as 1.65 cm, and the beam length of 1 m. For the remote experimental procedures, we suggest two distinct control structures that are based on PD control action. The first one considers a cascade controller with an integrator as a highfrequency filter. Additionally, the beam angle control was performed, as shown in Figure 11. The second one with the derivative part in the feedback loop, as shown in Figure 12. The second structure with the derivative term in the feedback loop is an interesting approach, as students can analyze and understand how to reduce the noise amplification.   Figure 13 presents the sequence for the experimental approach. The basis of the study is structured in four main stages to be followed: (i) to simulate the models presented in Section 2 of this article, the simplified model (3), and the full model (13); (ii) to use a classical controller tuning rule to define the controllers' parameters. Some interesting methods for this purpose can be found in [63]; (iii) to study MOOP designed for control engineering problems, as described in Section 4. In the sequence, students must define the cost functions, those that are adopted to minimize the performance measures described in Section 3. In this case, MOEA and MCDM methods are both considered; and, (iv) to perform tests of the controllers tuned through distinct methods for the ball and beam WebLab system, when comparing the simulated solutions with the real system performance. Figure 13. Experimental procedures to be developed by the students.

Results and Discussion
This section is dived into two subsections in order to illustrate all of the procedures described in Section 4 associated with multiobjective optimization, and results that were obtained while applying the NSGA-II algorithm. The first one shows how the three steps, those that were addressed in Section 4 for multiobjective optimization, were applied to the ball and beam control problem. In the sequence, the results were presented and discussed, showing promising results in assuming EMO strategies for control tuning.

Multiobjective Optimization Procedures
This subsection describes the procedures applied to the ball and beam process to optimize the controller gains that are based on MOOP.

Problem Definition
The authors assumed the parametric model presented in Figure 11, once it was possible to compare the multiobjective approach to a reference controller already described in [25], in order to define the multiobjective problem associated with the control of the ball and beam apparatus. In this  Figure 13 presents the sequence for the experimental approach. The basis of the study is structured in four main stages to be followed: (i) to simulate the models presented in Section 2 of this article, the simplified model (3), and the full model (13); (ii) to use a classical controller tuning rule to define the controllers' parameters. Some interesting methods for this purpose can be found in [63]; (iii) to study MOOP designed for control engineering problems, as described in Section 4. In the sequence, students must define the cost functions, those that are adopted to minimize the performance measures described in Section 3. In this case, MOEA and MCDM methods are both considered; and, (iv) to perform tests of the controllers tuned through distinct methods for the ball and beam WebLab system, when comparing the simulated solutions with the real system performance.  Figure 13 presents the sequence for the experimental approach. The basis of the study is structured in four main stages to be followed: (i) to simulate the models presented in Section 2 of this article, the simplified model (3), and the full model (13); (ii) to use a classical controller tuning rule to define the controllers' parameters. Some interesting methods for this purpose can be found in [63]; (iii) to study MOOP designed for control engineering problems, as described in Section 4. In the sequence, students must define the cost functions, those that are adopted to minimize the performance measures described in Section 3. In this case, MOEA and MCDM methods are both considered; and, (iv) to perform tests of the controllers tuned through distinct methods for the ball and beam WebLab system, when comparing the simulated solutions with the real system performance. Figure 13. Experimental procedures to be developed by the students.

Results and Discussion
This section is dived into two subsections in order to illustrate all of the procedures described in Section 4 associated with multiobjective optimization, and results that were obtained while applying the NSGA-II algorithm. The first one shows how the three steps, those that were addressed in Section 4 for multiobjective optimization, were applied to the ball and beam control problem. In the sequence, the results were presented and discussed, showing promising results in assuming EMO strategies for control tuning.

Multiobjective Optimization Procedures
This subsection describes the procedures applied to the ball and beam process to optimize the controller gains that are based on MOOP.

Problem Definition
The authors assumed the parametric model presented in Figure 11, once it was possible to compare the multiobjective approach to a reference controller already described in [25], in order to define the multiobjective problem associated with the control of the ball and beam apparatus. In this Figure 13. Experimental procedures to be developed by the students.

Results and Discussion
This section is dived into two subsections in order to illustrate all of the procedures described in Section 4 associated with multiobjective optimization, and results that were obtained while applying the NSGA-II algorithm. The first one shows how the three steps, those that were addressed in Section 4 for multiobjective optimization, were applied to the ball and beam control problem. In the sequence, the results were presented and discussed, showing promising results in assuming EMO strategies for control tuning.

Multiobjective Optimization Procedures
This subsection describes the procedures applied to the ball and beam process to optimize the controller gains that are based on MOOP.

Problem Definition
The authors assumed the parametric model presented in Figure 11, once it was possible to compare the multiobjective approach to a reference controller already described in [25], in order to define the multiobjective problem associated with the control of the ball and beam apparatus. In this way, a minimization problem was assumed to optimize the controller gains.
where ϕ = K p , K d , K h , K p1 are the controller's gain vector, being K p and K d the proportional and derivative gains, K h , the high-frequency filter, and K p1 the proportional gain that controls the beam angle. All of the gains are restricted to the [0, 50] interval, since the critical gain of the system transfer function tends to infinity. The selected interval was used to simplify the example regarding the search space.
The two cost functions that were presented in Equation (19) were considered. As described in Section 3, the first one represents an error criterion, the Integral of Absolute Error (IAE), and the second one is the control robustness criteria, the Integral of the Absolute Variation of Control Signal (IAVU). While assuming what is presented in Section 4.1, to select both the error and the robustness criteria in the optimization statement improves the confidence in the obtained controller.

Multiobjective Optimization through NSGA-II
Due to interesting results that were presented by the NSGA-II in works that are associated with control applications [64][65][66], this algorithm was assumed in this study. We used the classical multiobjective algorithm NSGA-II, which was described in Algorithm 1 [57], to tune the PD controller applied to the ball and beam apparatus. The PD controller was selected, since the transfer function of the full model described in Equation (13) has two integrators. A population size of 200 was assumed, and the stopping criterion was based on the number of generations, in this case, being equal to 100. The probability of genetic operations was applied, as proposed in [57], for crossover 90% and, for mutation, the inverse of the number of decision variables, i.e., 25%.

Multicriteria Decision Making Strategy
A set of non-dominated solutions was obtained as the considered optimization strategy is an EMO algorithm. The normalized approximation result was obtained through a linear transformation and is shown in Figure 14 at every 20 iterations, showing the method convergence. The axes of Figure 14 represent the selected cost functions IAE and IAUV.
While considering the set of solutions in the Pareto Front, the Nash solution, as described in Section 4.3, was selected and compared to the reference controller that was provided by the manufacturer [25]. Furthermore, for comparison purposes, the Åström's robust loop shaping method [29] was applied, and the controller's gains of the method were obtained while considering a joint sensitivity of 1.4. This value was found on the previously cited reference for a double integer system (see Table 1). We considered, for the time delay, t m , a small approximation of 10 ms, since the method needs a time delay to obtain the controller's gains. This approximation is one-tenth of the sampling time, and it was just assumed for calculation purpose, as it does not interfere in the system dynamics. For the gain K m , the gravitational force in the city of Curitiba (Brazil) equal to 9.78 m/s 2 was considered.   Table 2 presents the controllers' gain and the cost function values for all of the controllers. As can be verified in Table 2, the optimized controller performed better than the other controllers for both of the cost functions, with a reduction of 44.87% in the reference tracking error and 20.17% in the control effort when compared to the reference controller. Additionally, a reduction of 45.46% and 17.82% for the error and the control effort were obtained, respectively, when compared to the Åström robust loop shaping method.
In Figure 15, comparisons were performed while considering changes in the desired position of the ball. The output signal of the ball and beam process, which represents the real ball position, was   Table 2 presents the controllers' gain and the cost function values for all of the controllers. As can be verified in Table 2, the optimized controller performed better than the other controllers for both of the cost functions, with a reduction of 44.87% in the reference tracking error and 20.17% in the control effort when compared to the reference controller. Additionally, a reduction of 45.46% and 17.82% for the error and the control effort were obtained, respectively, when compared to the Åström robust loop shaping method.
In Figure 15, comparisons were performed while considering changes in the desired position of the ball. The output signal of the ball and beam process, which represents the real ball position, was presented while considering the three controllers. Moreover, Figure 16 represents the control signal that is associated with the results that were obtained in Figure 15. In this case, it is important to emphasize the control effort of each strategy. It is also possible to verify that the optimized controller presents an overshoot of 2.6% versus 10.92% for the Åström Controller. Even observing that the Reference Controller did not present an overshoot, its settling time, reduction of the error to a value lesser than 2% from reference, is 2.37 s. When the optimized controller is analyzed in terms of settling time, just 0.83 s was obtained, while the Åström controller provided a 5.30 s for this criterion in positive reference changes.  As previously presented, the use of multiobjective optimization procedures in control engineering can be used to improve the performance of the system. The combination of these procedures can be assumed to be a generalist tool in control tuning problems. A more specific analysis of the ball and beam control case presented in this study showed that the application of the multiobjective optimization procedures that are described in Section 4 increases the possibility of considering distinct control objectives. These objectives could be the classical ones presented in Figure 6 (Section 3) or additional ones that can be defined according to optimization requirements that are associated with the problem.
When remote experimentation is considered, multiobjective optimization can be used to improve students' knowledge and the quality of the courses that are associated with industrial  As previously presented, the use of multiobjective optimization procedures in control engineering can be used to improve the performance of the system. The combination of these procedures can be assumed to be a generalist tool in control tuning problems. A more specific analysis of the ball and beam control case presented in this study showed that the application of the multiobjective optimization procedures that are described in Section 4 increases the possibility of considering distinct control objectives. These objectives could be the classical ones presented in Figure 6 (Section 3) or additional ones that can be defined according to optimization requirements that are associated with the problem.
When remote experimentation is considered, multiobjective optimization can be used to improve students' knowledge and the quality of the courses that are associated with industrial As previously presented, the use of multiobjective optimization procedures in control engineering can be used to improve the performance of the system. The combination of these procedures can be assumed to be a generalist tool in control tuning problems. A more specific analysis of the ball and beam control case presented in this study showed that the application of the multiobjective optimization procedures that are described in Section 4 increases the possibility of considering distinct control objectives. These objectives could be the classical ones presented in Figure 6 (Section 3) or additional ones that can be defined according to optimization requirements that are associated with the problem.
When remote experimentation is considered, multiobjective optimization can be used to improve students' knowledge and the quality of the courses that are associated with industrial processes simulations and its challenges, especially when a new and specific process has to be inserted into the academic environment. Moreover, remote experimentation provides a more complex environment in terms of technology integration, as students will not only work with control theory, but with data acquisition, data communication, and virtual instrumentation, technologies that are being assumed in the recent context of industry 4.0.

Conclusions
This research started presenting an overview of the essential learning methods that are used in the control engineering learning, as well as the importance that is associated with the approximation of the academic to the industrial area. As the next step, a brief description of the ball and beam process, which represents an example of a complex system that is frequently found in industry, was addressed. The third section presented the PID control theory, showing two variations that can be found in the literature to control the ball and beam system. Additionally, metrics that were commonly used as performance measures related to both the error and the control effort were described. Section 4 presents the multiobjective optimization approach that was related to control engineering problems, while Section 5 showed the experimental structure for the WebLab and experimental procedures assumed in this study. Finally, Section 6 reported how the presented method is effective in comparison with a classical control tuning approach, emphasizing the advantages of applying multiobjective optimization in the educational context that is associated with control engineering.
The use of WebLabs in education is an interesting tool for improving the quality of control engineering courses, since it is possible to improve the diversity of experiments and approximate the challenges of industrial areas to the academic space. The present research proposed a sequential method for the use of WebLabs that is associated with PID tuning problems, where the use of computational intelligence through multiobjective optimization design was applied.
The control of the ball and beam plant, which is a nonlinear and unstable open-loop process, is a common problem that students will find in industry, and the comparison between the classical direct tuning method with the multiobjective optimization strategy provides a better understanding of the advantages of assuming EMO techniques for similar problems that are associated to control systems.
WebLabs can be used to improve the availability of the laboratory structures for experimentation processes while considering a collaboration network. The implementation costs for this solution is around US$3000.00. The ball and beam plant cost is US$2500.00, the lab server, US$60.00, the cloud server, US$200.00/year, and other devices can be estimated at US$240.00. These values are related to the implementation of the structure in a Brazilian university.
For future works, we will expand on the use of intelligent control systems through multiobjective optimization for other processes, while also including the use of fuzzy systems for specialist control, and machine learning methods that can be used in the identification procedures while focusing on the obtention of parametric models, these procedures that are always associated with virtual laboratories. Additionally, a formal analysis focusing on the feedback provided by the students will also be performed. Previous results that are associated with the students' feedback indicated that this type of methodology is understood as the future of engineering learning, once the industry demand nowadays is directly associated with connected machines, integrated systems, and real-time information analysis.