Fuzzy Robust Design of Dynamic Vibration Absorbers

This paper is dedicated to the development of robust optimization anddecisionmaking techniques taking into account the uncertain parameters of linear and nonlinear dynamic vibration absorbers. In this case, novel approaches are proposed regarding different fuzzy logic optimization strategies. The uncertain parameters of the considered mechanical systems are treated as fuzzy variables. Consequently, the associated optimization problem is described as a fuzzy function that maps the fuzzy inputs. The proposed techniques are applied to the design of dynamic vibration absorbers.This numerical study illustrates the versatility and convenience of the proposed fuzzy logic optimization strategies.


Introduction
The design process and operation limits of mechanical systems are currently obtained from the analysis of deterministic models, which should be able to represent the associated dynamic phenomena.However, parametric or nonparametric uncertainties are inherent to dynamic systems due to manufacturing errors and tolerances, damage, wear, and influence of environmental conditions that should be taken into account to obtain reliable mathematical models.These fluctuations are disregarded in ordinary deterministic approaches.In this sense, deterministic models can be considered inaccurate and should be replaced by models that take into account the uncertainties.
There are innumerous approaches for modeling uncertainties in mechanical systems.Among them, both the stochastic and fuzzy logic approaches are widely used.The fuzzy logic [1] is an intuitive approach that takes into account associated fuzzy sets and the possibility theory [2], while the stochastic methodology is based on the probability theory [3][4][5].Both methodologies have been successfully applied to various problems in engineering, leading to similar results.The stochastic approach shows a longer history of applications to mechanical systems [6][7][8][9][10] and, consequently, seems to be more established than the fuzzy logic.However, the fuzzy approach is receiving more attention in our days and the literature brings some successful applications [11][12][13][14].
The optimum design of a given mechanical system must be robust with respect to unavoidable uncertainties that affect its dynamic behavior.Regarding the optimization issues, robustness can be understood as a design that is frequently not the best solution obtained for a particular dynamic problem.However, the robust solution is a better solution for the cases in which the system is exposed to parameter fluctuations.Therefore, a robust design can be considered as the one wherein the goal is to obtain the best solution for a given uncertain scenario.The uncertainty analysis is associated with an optimization procedure in order to determine the robust design for the system (i.e., the robust optimization procedure).
Taguchi [15] is the pioneer in the area of robust optimization.He proposed the insertion of noise factors in the cost functions, thus defining a mean square deviation that has to be maximized in terms of the design variables.Commonly, the robust optimization approaches rely on stochastic techniques based on the statistical moments of the problem [16][17][18][19][20][21] and in the definition of a robust objective function.This approach leads to large combinatorial multiobjective problems, imposing a serious disadvantage regarding the associated computational cost.The fuzzy logic is a simpler mathematical logic for uncertainty analysis, as based on pessimistic and optimistic values (i.e., fuzzy sets).Therefore, the fuzzy logic seems to be well suited for robust optimization.
In this context, the present contribution aims at proposing different fuzzy robust optimization techniques.These methodologies were evaluated for the design of linear and nonlinear robust dynamic vibration absorbers (DVA and nDVA, resp.), considering parametric uncertainties affecting the mechanical system.The uncertainties are treated as fuzzy variables so that the optimization model becomes a fuzzy function that maps the fuzzy inputs.Therefore, a fuzzy evaluation of the optimization problem can be conceived through the evaluation of the pessimistic and optimistic values of the fuzzy sets (i.e., the uncertain information).This is the basic concept of the techniques conveyed in the present contribution.
This paper is organized as follows.First, the key concepts of fuzzy logic are discussed; then, the fuzzy robust optimization techniques are presented.In the sequence, two design problems are described and the robust optimization results are discussed.Finally, a few remarks are made and further research works are proposed.

Fuzzy Logic
Fuzzy logic is the term used for the intuitive formulation proposed by Zadeh [1] that generalizes the membership notion of an element in a mathematical set, creating a graded notion represented by a membership function.In a classical sense, a math set is the collection of elements that belongs to a given definition.Any element can belong or not to the mentioned set.Regarding the fuzzy approach, the Boolean notion of the membership function is replaced by a graded one.Thus, an element can belong, not belong, or partially belong to the so-called fuzzy set.
As presented in [22], the membership function of a fuzzy set can be basically understood through three different perspectives, as follows: the degree of similarity, in which the membership function quantifies how similar an element of the set is, as compared to a prototype element; the degree of preference, where the membership function quantifies the preference or the feasibility of an element of the set; the degree of uncertainty, in which the membership function is the degree of possibility that a determined parameter  assumes the value of a specific element  of the set.
Based on the concept of degree of uncertainty, a fuzzy set can be understood as the union of the pessimistic and optimistic values of an uncertain parameter considering each level of possibility (i.e., -levels of the membership function).It is important to point out that the membership function is a possibility distribution, not a probability distribution as used in the case of the stochastic theory.Additionally, possibility is the measure of whether an event can happen, while probability is a measure of whether an event will happen.Therefore, the possibility distribution of a given parameter  quantifies the possible values that this parameter can assume.Differently, a probability distribution quantifies the chances of the parameter  to assume a certain value .In the fuzzy set theory, the uncertainties are described as fuzzy inputs and the mathematical model is represented as a fuzzy function that maps those inputs.The extension principle [23] is a widely known method for mapping fuzzy inputs, in which the mapping process is made through the evaluation of an associated mathematical function.Another method used for mapping the fuzzy inputs is the so-called -level optimization [11].In this case, the fuzzy inputs are discretized by means of -cuts, thus creating a crisp subspace that is the search space of an associated optimization problem.The optimization process is carried out in order to find the minimum and maximum values of the model output for each -level (i.e., lower and upper limits of the correspondent -level, resp.; see Figure 1).The -level optimization is an eficient methodology for mapping fuzzy inputs.However, the effectiveness of the method is highlly related to the performance of the selected optimization technique.
The mapping of a fuzzy input generates an associated fuzzy output.Considering a decision making system, the fuzzy output should be defuzzified into a single value (i.e., a crisp value) that represents the uncertain information.Therefore, the defuzzification is the process of aggregation of the fuzzy information in a single value [24] that can be used to generate robust data.
Different defuzzification techniques exist in the literature to obtain robust data, as follows.The largest of the maximum method selects as crisp value the largest value among all realizations of the fuzzy number.The smallest of the maximum method chooses from the maximum membership function values the smallest one as crisp value.The mean of maximum method takes as crisp value the one obtained from the mean value between the smallest and largest maxima.However, the most commonly used defuzzification techniques are the so-called centroid and bisection methods [25].The centroid method determines the robust value from the area center (gravity center) of the membership function.Similarly, the bisection method calculates the bisector of the fuzzy output to obtain the crisp value.Figure 2 illustrates the presented defuzzification techniques.These aspects are the key points of fuzzy logic and fuzzy uncertainty analysis.

Fuzzy Robust Optimization
In general, robust optimization is performed relying on the stochastic theory, where a set of samples of the design variables (i.e., a population) is generated and robustness is evaluated by means of the statistical analysis of this set.The stochastic theory is an efficient robust optimization approach, especially when associated with heuristic optimization algorithms (genetic algorithms or differential evolution [26,27]) and sampling methods such as the Latin Hypercube [28].However, this procedure usually creates a large combinatorial problem imposing high computational costs.
The fuzzy uncertainty analysis can be used as an alternative approach to reduce the mentioned computational cost and, consequently, to increase the efficiency of the robust optimization process.Fuzzy uncertainty analysis through the -level optimization requires fewer evaluations of the cost function, presenting lower computational cost as compared with the stochastic approach.The associated minimization and maximization processes can be performed by using direct optimization procedures, such as the SQP (Sequential Quadratic Programing) [29], thus reducing even more the computational effort.Therefore, the development of fuzzy metrics to evaluate the robustness of a given design leads to an efficient tool for robust optimization.In this context, the next section is devoted to the discussion of novel robust optimization methodologies based on the fuzzy approach.of a model as fuzzy variables, a robust optimum should take into account the possible realizations of the uncertain values.In this context, a robust optimization procedure has to consider the pessimistic and optimistic values of the uncertain parameters during the optimization tasks.
A dual optimization problem can be conceived by treating the optimization problem as a fuzzy function and following a methodology similar to the -level optimization.The pessimistic and optimistic boundaries (i.e., the lower and upper bounds, resp.)associated with the uncertain parameters are created for each -level.A fuzzy optimum is then obtained considering two different optimization problems: (i) an optimization problem where the goal is to find an optimum by considering the pessimistic values of the uncertain information for the related -level and (ii) a second optimization problem dedicated to searching an optimum by taking into account the optimistic values for the same -level.This fuzzy optimum is a set of possible optima, that is, a fuzzy information that can be aggregated into a robust optimum through the defuzzification process.The dual optimization procedure is illustrated in Figure 3.
Defining x as the vector of deterministic design variables and p as the vector of uncertain values, a pseudo-objective function (sequential unconstrained minimization technique [29]) is expressed as where  is the vector of the penalty factors related with the penalty functions,  is the vector of the original cost function,  is the vector of constraints, and Φ is the resulting pseudoobjective function.
Based on the possibility theory, the dual optimization problem can be formulated as presented in (2) and (3).Consider the following: in which p and p represent the lower and upper bounds of the uncertain values, respectively,  and  are the lower and upper bounds of the cost functions vector,  and  are the bounderies of the constraints vector, and Φ and Φ are the pessimistic and optimistic pseudo-objective functions, respectively.
The optimum for the optimization problem expressed by (2) leads to a lower bound optimal value x * , while the optimization problem described by (3) results in an upper bound optimal value x * .Then, at each -level two sets of optimal results are generated, namely, the lower optima and the upper optima (i.e., pessimistic and optimistic values, resp.).The union of the two sets generates a fuzzy set or a fuzzy optimum, which possibly contain the optimal solutions of the design problem.This fuzzy information can be addressed to a robust solution through a defuzzification procedure, as previouslly mentioned.

Amplitude Minimization.
In this case, it is assumed that the robustness metric can be defined regarding the amplitude of the envelope determined by the cost functions given by ( 2) and (3).The higher the amplitude of the obtained envelope (i.e., the difference between the optimistic and pessimistic responses) the lower the robustness of the solution.Consequently, a more robust solution will present envelopes with lower amplitude (i.e., smaller difference between the optimistic and pessimistic responses).Thus, the robustness metric that measures the distance between the optimistic and pessimistic responses (measure of the response amplitude) is now defined.Equation (4) presents the proposed robustness metric for the amplitude minimization scheme: where  is the total number of points in which the pseudoobjective functions are evaluated.The distance between the optimistic and the pessimistic responses is then computed by using (4).
Therefore, a new objective function  (5) is formulated to find the solution associated with the highest robustness value that minimizes the proposed metric (4).Figure 4  (5) uncertainty analysis is performed generating a fuzzy output.As previously mentioned, this fuzzy output is the union of the pessimistic and optimistic responses, in which the level boundaries represent the lower and the upper bounds of all possible outcomes of the optimization problem.Consequently, it can be inferred that a robust solution will have the smallest difference between these bounds and therefore the smallest area for its fuzzy output.

Fuzzy Response
Thus, in the fuzzy response minimization scheme, the robust optimization procedure is formulated aiming at finding the robust optimum, which is represented by the samallest area of the fuzzy output (see Figure 1).Defining  as the function that represents the optimization problem, the uncertainty analysis can be defined by  (x, p) =  (x,  ()) . ( The area of the fuzzy output can be viewed as the robustness metric wherein its minimization leads to a robust optimum.Therefore, the metric is calculed from (7). Figure 5 illustrates the main concepts of the fuzzy response minimization procedure All the robust optimization procedures conceived in this contribution are formulated as single-objective optimization problems in which the objective function is directly related to the robustness of the problem.However, all these procedures can be easily extended to multiobjective optimization problems by considering a cost function related to the original optimization problem and another one related to the robustness evaluation (one of the proposed criteria).
In order to model an uncertain quantity, the perspective of the degree of uncertainty is very interesting.However, the degree of preference seems to be more appropriate for optimization purposes.Nevertheless, in the fuzzy robust optimization criteria presented in this paper, the interpretation of the degree of uncertainty is used both to model the uncertain parameters (fuzzy variables) and to define the optimization task.
In addition, the centroid method was used for defuzzification, once this method is the closest to what can be defined as the weighted average of a fuzzy set.

Linear Vibration Absorber (DVA).
A linear dynamic vibration absorber (DVA) is a mechanical system composed by inertia, stiffness, and damping elements, which is attached to a primary system aiming at attenuating the vibration response [31].The classical one-degree-of-freedom DVA coupled to the primary system is presented by Figure 6.
The vibration amplitude of the primary system vanishes when the natural frequency of the DVA matches the excitation frequency of the force acting on the primary system [31].The DVA is tuned by properly choosing its mechanical properties, so that the vibration energy of the primary system can be mitigated.Equation ( 8) presents the dimensionless vibration amplitude of the mechanical system presented in Figure 6.Consider where the following parameters are considered:  =   /  is the mass ratio,  =   /  is the turning factor,  = /  is the forcing frequency ratio,   = 2    is the critical damping,  =   /  is the damping ratio, (  ) st =  0 /  is the static displacement of the primary mass, and   is the displacement of the primary mass.In this case,   ,   , and   are the mass, damping factor, and stiffness coefficient of the primary system, respectively.For the DVA,   ,   , and   are the mass, damping factor, and stiffness coefficient;  is the frequency of the external force,   = √(  /  ) is the natural frequency of the DVA,   = √(  /  ) is the natural frequency of the primary system, and  0 is the amplitude of the external force.Mathematically, the optimum design of a DVA consists in determining the parameters  and  that minimize the maximum amplitude of the frequency response function of the complete system described by (8) (i.e., the primary system plus the DVA; see Figure 6).

Nonlinear Dynamic Vibration Absorber (nDVA).
The primary system, the DVA, or both systems may contain nonlinear elements (particularly the springs) that exhibit nonlinear behavior.In [30] the authors analyzed the limitations and capabilities of a nonlinear dynamic vibration absorber (nDVA).In this work, the robust design of a one-degree-offreedom nDVA was carried out by taking into account the response of the system as described in the frequency domain.
Equations ( 9) present a system of algebraic equations representing the dynamic behavior of an nDVA coupled to a nonlinear primary system (stiffness nonlinearity).This mathematical formulation was adopted in [30], in which the Van der Pol Transformation was used.Consider where the displacements  and velocities V of the primary system () and nDVA () (vibration responses   ,   , V  , and V  , resp.) are assumed to be slow functions of a normalized time.In this case, the formulation of new parameters is also required, as follows: , where   =  0 /  and    is the nonlinear stiffness coefficient of the spring ( = , ), and  =  0 /(     2 ).The solution of this algebraic system is obtained by defining the parameters  = / 2  , Ω = /  , so that in the algebraic system one has   = 1/Ω and   = /Ω, where  =   /  .
For the optimum design of the nDVA, a multiobjective optimization problem is formulated considering a vector  containing two objective functions.The first is devoted to the minimization of the vibration amplitude associated with the first vibration mode of the nonlinear damped system.The second objective function of the problem is related to the maximization of the suppression bandwidth, as achieved by the nDVA.Figure 7 illustrates the two objective functions of the considered optimization process.

Numerical Results
In this paper, the proposed fuzzy robust optimization techniques were applied to the design of the systems depicted in Figure 8 (DVA, Figure 8(a), and nDVA, Figure 8(b)).Regarding the DVA design, uncertainties affecting its stiffness   and damping   coefficients were taken into account.For the design of the nDVA, the parameters   , , , and  were considered uncertain.
Considering the robust design of the DVA, the uncertain parameters were modeled as fuzzy triangular numbers with an uncertainty of ±10% with respect to the nominal value of the parameter.Table 1 shows the uncertainty scenario considered in the present contribution.In the design of the nDVA, the uncertain parameters were also treated as fuzzy triangular number.However, an uncertainty range of ±30% was adopted with respect to the nominal values.Table 2 shows the uncertainty scenario considered for the nDVA design.Additionally, the masses for the primary and DVA subsystems are also given.It is worth mentioning that for the design of the DVA, the uncertainty associated with the stiffness of the primary system   also affects the parameters , , and the forcing frequency factor  (see (8)).In the following, the masses adopted for the DVA were the same as those used for the nDVA.
The simulation of the DVA was performed considering a frequency bandwidth of 5 Hz (i.e.,  = 5-10 Hz), with a sample frequency of 0.001 Hz.The design of the nDVA was carried out by taking into account a frequency bandwidth of 0.6 Hz (i.e., Ω = 0.7-1.3Hz), comprising 300 frequency steps.Regarding the solution of the algebraic system defined by ( 9), the initial conditions of the nDVA were considered as follows:   = 5 m,   = 5 m, V  = 2 m/s, and V  = 4 m/s.The DVA was considered at rest for  = 0. Table 3 presents the obtained results applying both the deterministic and robust optimization procedures related to the design of the DVA.The robustness metrics, presented in (4), were used to evaluate the performance of the proposed fuzzy techniques.Note that using the amplitude minimization procedure the smallest robust parameter resulted, which is similar to the one determined for the case in which the fuzzy response minimization was applied.Additionally, it can be observed that the dual optimization procedure led to the parameters with the lowest associated robustness.Regarding the cost function (i.e., the difference between the invariant points), the amplitude minimization procedure led to the worst results.The envelopes obtained for each robust procedure proposed in this work are presented in Figure 9.
Table 4 presents the results (taken from [30]) obtained considering the deterministic and robust optimization problem related to the design of the nDVA.The design of the nDVA is a multiobjective optimization problem that was solved by using the compromise programming technique [29] considering equal weight cost functions (i.e., the same weighting factor for both cost functions).The targets of the optimization problem were the results obtained by analyzing each cost function separately.The worst known value for each objective function was adopted as the value of each cost function evaluated at the initial point.Therefore, the multiobjective problem was handled as a representative monoobjective optimization problem, resulting in a single optimum instead of a Pareto frontier.
As observed in the DVA problem, the amplitude minimization technique presented the most robust optimum value, followed by the fuzzy response minimization procedure.In this case, the differences between the robustness metrics were most likely caused by numerical approximations.The dual optimization approach led to the worst result.In contrast with the design of the linear system, the robustness of the nDVA increases significantly by applying the proposed fuzzy techniques as compared with the deterministic solution.Concerning the cost functions (minimization of  1 and maximization of  2 ; see Table 4), all the proposed fuzzy methodologies presented better solutions than the ones (taken from [30]) obtained from the deterministic approach.Figure 10 presents the envelopes obtained for each proposed fuzzy approach for the design of the nDVA.

Concluding Remarks
In this paper, three robust optimization methodologies based on the fuzzy logic approach were proposed and evaluated for the robust design of linear and nonlinear dynamic vibration Amplitude (m) Amplitude (m) Amplitude (m) The results demonstrated that, for all the considered design problems, the proposed methodologies found solutions that are more robust than the ones obtained from   the deterministic approach.Additionally, the so-called amplitude minimization procedure obtained the most robust solution for all considered design problems.Finally, the methodologies conveyed demonstrate that fuzzy optimization techniques are viable for robust optimization problems and may even be more suitable for this kind of design problem as compared with other competing alternatives, due to its mathematical simplicity and less number of evaluations of the objective function.

Figure 4 :
Figure 4: Flowchart of the amplitude minimization procedure.

Figure 9 :
Figure 9: Robustness envelopes determined in the design of the DVA.

Figure 10 :
Figure 10: Robustness envelopes determined in the design of the DVA.

Table 1 :
Parameters used for the design of the linear DVA.Lower stands for lower limit.Upper stands for upper limit. *

Table 2 :
Parameters used for the design of the nDVA.
* Lower stands for lower limit.Upper stands for upper limit.

Table 3 :
Optimization results obtained by the DVA design problem.

Table 4 :
Optimization results obtained by the nDVA design problem.