# Fuzzy Robust Design of Dynamic Vibration Absorbers

Authors: A. D. G. Silva; A. A. Cavalini Jr; V. Steffen Jr
Journal: Shock and Vibration
 (2016)

---

## Abstract

This paper is dedicated to the development of robust optimization and decision making 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.


---

## Body

## 1. 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–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–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–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–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.


## 2. 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 u assumes the value of a specific element x 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 u quantifies the possible values that this parameter can assume. Differently, a probability distribution quantifies the chances of the parameter u to assume a certain value x.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.Figure 1
Theα-level optimization.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., acrisp 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. Thelargest of the maximum method selects ascrisp value the largest value among all realizations of the fuzzy number. Thesmallest of the maximum method chooses from the maximum membership function values the smallest one ascrisp value. Themean of maximum method takes ascrisp value the one obtained from the mean value between the smallest and largest maxima. However, the most commonly used defuzzification techniques are the so-calledcentroid andbisection methods[25]. Thecentroid method determines the robust value from the area center (gravity center) of the membership function. Similarly, thebisection methodcalculates the bisector of the fuzzy output to obtain thecrisp value. Figure 2 illustrates the presented defuzzification techniques. These aspects are the key points of fuzzy logic and fuzzy uncertainty analysis.Figure 2
Defuzzification techniques.


## 3. 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.


### 3.1. Dual Optimization

As mentioned, the fuzzy set can be considered as the union of the pessimistic and optimistic values of a given fuzzy variable. Considering the uncertainties 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.Figure 3
Dual optimization procedure flowchart.Definingx 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(1)
Φ
x
,
p
~
,
r
p
=
F
x
,
p
~
+
r
p
∗
P
x
,
p
~
,where r
p is the vector of the penalty factors related with the penalty functions, F is the vector of the original cost function, P is the vector of constraints, and Φ is the resulting pseudo-objective function.Based on the possibility theory, the dual optimization problem can be formulated as presented in (2) and (3). Consider the following: (2)
Φ
_
x
,
p
_
,
r
p
=
F
_
x
,
p
_
+
r
p
∗
P
_
x
,
p
_
,
(3)
Φ
¯
x
,
p
¯
,
r
p
=
F
¯
x
,
p
¯
+
r
p
∗
P
¯
x
,
p
¯in which p
_ and p
¯ represent the lower and upper bounds of the uncertain values, respectively, F
_ and F
¯ are the lower and upper bounds of the cost functions vector, P
_ and P
¯ 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.


### 3.2. 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:(4)
r
o
i
=
Φ
¯
i
x
,
p
¯
,
r
p
-
Φ
_
i
x
,
p
_
,
r
p
2
i
=
1
,
…
,
n
,where n is the total number of points in which the pseudo-objective functions are evaluated. The distance between the optimistic and the pessimistic responses is then computed by using (4).Therefore, a new objective functionz (5) is formulated to find the solution associated with the highest robustness value that minimizes the proposed metric (4). Figure 4 summarizes the amplitude minimization procedure(5)
z
x
,
p
~
,
r
p
=
min
⁡
∑
i
=
1
n
r
o
i
.Figure 4
Flowchart of the amplitude minimization procedure.


### 3.3. Fuzzy Response Minimization

In the case of a robust optimization problem (considering uncertainties), for a given design (the configuration of the design variablesx) a fuzzy 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.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 Figure1). Defining Z as the function that represents the optimization problem, the uncertainty analysis can be defined by(6)
Z
x
,
p
~
=
Z
x
,
f
α
.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(7)
A
r
e
a
=
∫
0
1
Z
x
,
f
α
d
α
=
∫
0
1
Z
_
x
,
f
_
α
d
α
+
∫
0
1
Z
¯
x
,
f
¯
α
d
α
.Figure 5
Flowchart for 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, thecentroid method was used for defuzzification, once this method is the closest to what can be defined as the weighted average of a fuzzy set.


## 3.1. Dual Optimization

As mentioned, the fuzzy set can be considered as the union of the pessimistic and optimistic values of a given fuzzy variable. Considering the uncertainties 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.Figure 3
Dual optimization procedure flowchart.Definingx 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(1)
Φ
x
,
p
~
,
r
p
=
F
x
,
p
~
+
r
p
∗
P
x
,
p
~
,where r
p is the vector of the penalty factors related with the penalty functions, F is the vector of the original cost function, P is the vector of constraints, and Φ is the resulting pseudo-objective function.Based on the possibility theory, the dual optimization problem can be formulated as presented in (2) and (3). Consider the following: (2)
Φ
_
x
,
p
_
,
r
p
=
F
_
x
,
p
_
+
r
p
∗
P
_
x
,
p
_
,
(3)
Φ
¯
x
,
p
¯
,
r
p
=
F
¯
x
,
p
¯
+
r
p
∗
P
¯
x
,
p
¯in which p
_ and p
¯ represent the lower and upper bounds of the uncertain values, respectively, F
_ and F
¯ are the lower and upper bounds of the cost functions vector, P
_ and P
¯ 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.


## 3.2. 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:(4)
r
o
i
=
Φ
¯
i
x
,
p
¯
,
r
p
-
Φ
_
i
x
,
p
_
,
r
p
2
i
=
1
,
…
,
n
,where n is the total number of points in which the pseudo-objective functions are evaluated. The distance between the optimistic and the pessimistic responses is then computed by using (4).Therefore, a new objective functionz (5) is formulated to find the solution associated with the highest robustness value that minimizes the proposed metric (4). Figure 4 summarizes the amplitude minimization procedure(5)
z
x
,
p
~
,
r
p
=
min
⁡
∑
i
=
1
n
r
o
i
.Figure 4
Flowchart of the amplitude minimization procedure.


## 3.3. Fuzzy Response Minimization

In the case of a robust optimization problem (considering uncertainties), for a given design (the configuration of the design variablesx) a fuzzy 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.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 Figure1). Defining Z as the function that represents the optimization problem, the uncertainty analysis can be defined by(6)
Z
x
,
p
~
=
Z
x
,
f
α
.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(7)
A
r
e
a
=
∫
0
1
Z
x
,
f
α
d
α
=
∫
0
1
Z
_
x
,
f
_
α
d
α
+
∫
0
1
Z
¯
x
,
f
¯
α
d
α
.Figure 5
Flowchart for 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, thecentroid method was used for defuzzification, once this method is the closest to what can be defined as the weighted average of a fuzzy set.


## 4. Dynamic Vibration Absorbers: A Review

### 4.1. 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.Figure 6
Classical DVA scheme.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(8)
X
p
X
p
s
t
=
2
ζ
g
2
+
g
2
-
f
2
2
2
ζ
g
2
g
2
-
1
+
μ
g
2
2
+
μ
f
2
g
2
-
g
2
-
1
g
2
-
f
2
2
2
,where the following parameters are considered: μ
=
m
a
/
m
p is the mass ratio, f
=
ω
a
/
ω
p is the turning factor, g
=
ω
/
ω
p is the forcing frequency ratio, c
c
=
2
m
a
ω
p is the critical damping, ζ
=
c
a
/
c
c is the damping ratio, X
p
s
t
=
F
0
/
k
p is the static displacement of the primary mass, and X
p is the displacement of the primary mass. In this case, m
p, c
p, and k
p are the mass, damping factor, and stiffness coefficient of the primary system, respectively. For the DVA, m
a, c
a, and k
a are the mass, damping factor, and stiffness coefficient; ω is the frequency of the external force, ω
a
=
√
(
k
a
/
m
a
) is the natural frequency of the DVA, ω
p
=
√
(
k
p
/
m
p
) is the natural frequency of the primary system, and F
0 is the amplitude of the external force.Mathematically, the optimum design of a DVA consists in determining the parametersf 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).


### 4.2. 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-of-freedom 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 theVan der Pol Transformation was used. Consider(9)
1
+
μ
-
ω
p
2
u
p
+
μ
u
a
-
2
ζ
p
ω
p
v
p
-
3
4
ζ
p
u
p
u
p
2
+
v
p
2
+
F
=
0
,
μ
u
p
+
μ
-
μ
ω
a
2
u
a
-
μ
2
ζ
a
ω
a
v
a
-
3
4
μ
ζ
a
u
a
u
a
2
+
v
a
2
=
0
,
ω
p
2
-
1
-
μ
v
p
+
μ
v
a
-
2
ζ
p
ω
p
u
p
+
3
4
ζ
p
v
p
u
p
2
+
v
p
2
=
0
,
μ
v
p
+
μ
-
μ
ω
a
2
v
a
-
μ
2
ζ
a
ω
a
u
a
-
3
4
μ
ζ
a
v
a
u
a
2
+
v
a
2
=
0
,where the displacements u and velocities v of the primary system (p) and nDVA (a) (vibration responses u
p, u
a, v
p, and v
a, resp.) are assumed to be slow functions of a normalized time. In this case, the formulation of new parameters is also required, as follows: ζ
i
=
k
i
n
l
x
c
2
/
(
m
i
ω
2
), where x
c
=
F
0
/
k
p and k
i
n
l is the nonlinear stiffness coefficient of the spring (i
=
p, a), and F
=
F
0
/
(
m
p
x
c
ω
2
). The solution of this algebraic system is obtained by defining the parameters β
=
F
/
ω
p
2, Ω
=
ω
/
ω
p, so that in the algebraic system one has ω
p
=
1
/
Ω and ω
a
=
ρ
/
Ω, where ρ
=
ω
a
/
ω
p.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. Figure7 illustrates the two objective functions of the considered optimization process.Figure 7
Objective functions of the nDVA design problem (adapted from [30]).


## 4.1. 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.Figure 6
Classical DVA scheme.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(8)
X
p
X
p
s
t
=
2
ζ
g
2
+
g
2
-
f
2
2
2
ζ
g
2
g
2
-
1
+
μ
g
2
2
+
μ
f
2
g
2
-
g
2
-
1
g
2
-
f
2
2
2
,where the following parameters are considered: μ
=
m
a
/
m
p is the mass ratio, f
=
ω
a
/
ω
p is the turning factor, g
=
ω
/
ω
p is the forcing frequency ratio, c
c
=
2
m
a
ω
p is the critical damping, ζ
=
c
a
/
c
c is the damping ratio, X
p
s
t
=
F
0
/
k
p is the static displacement of the primary mass, and X
p is the displacement of the primary mass. In this case, m
p, c
p, and k
p are the mass, damping factor, and stiffness coefficient of the primary system, respectively. For the DVA, m
a, c
a, and k
a are the mass, damping factor, and stiffness coefficient; ω is the frequency of the external force, ω
a
=
√
(
k
a
/
m
a
) is the natural frequency of the DVA, ω
p
=
√
(
k
p
/
m
p
) is the natural frequency of the primary system, and F
0 is the amplitude of the external force.Mathematically, the optimum design of a DVA consists in determining the parametersf 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).


## 4.2. 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-of-freedom 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 theVan der Pol Transformation was used. Consider(9)
1
+
μ
-
ω
p
2
u
p
+
μ
u
a
-
2
ζ
p
ω
p
v
p
-
3
4
ζ
p
u
p
u
p
2
+
v
p
2
+
F
=
0
,
μ
u
p
+
μ
-
μ
ω
a
2
u
a
-
μ
2
ζ
a
ω
a
v
a
-
3
4
μ
ζ
a
u
a
u
a
2
+
v
a
2
=
0
,
ω
p
2
-
1
-
μ
v
p
+
μ
v
a
-
2
ζ
p
ω
p
u
p
+
3
4
ζ
p
v
p
u
p
2
+
v
p
2
=
0
,
μ
v
p
+
μ
-
μ
ω
a
2
v
a
-
μ
2
ζ
a
ω
a
u
a
-
3
4
μ
ζ
a
v
a
u
a
2
+
v
a
2
=
0
,where the displacements u and velocities v of the primary system (p) and nDVA (a) (vibration responses u
p, u
a, v
p, and v
a, resp.) are assumed to be slow functions of a normalized time. In this case, the formulation of new parameters is also required, as follows: ζ
i
=
k
i
n
l
x
c
2
/
(
m
i
ω
2
), where x
c
=
F
0
/
k
p and k
i
n
l is the nonlinear stiffness coefficient of the spring (i
=
p, a), and F
=
F
0
/
(
m
p
x
c
ω
2
). The solution of this algebraic system is obtained by defining the parameters β
=
F
/
ω
p
2, Ω
=
ω
/
ω
p, so that in the algebraic system one has ω
p
=
1
/
Ω and ω
a
=
ρ
/
Ω, where ρ
=
ω
a
/
ω
p.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. Figure7 illustrates the two objective functions of the considered optimization process.Figure 7
Objective functions of the nDVA design problem (adapted from [30]).


## 5. Numerical Results

In this paper, the proposed fuzzy robust optimization techniques were applied to the design of the systems depicted in Figure8 (DVA, Figure 8(a), and nDVA, Figure 8(b)). Regarding the DVA design, uncertainties affecting its stiffness k
a and damping c
a coefficients were taken into account. For the design of the nDVA, the parameters ξ
a, μ, β, and ρ were considered uncertain.Figure 8
Linear (DVA) and nonlinear (nDVA) dynamic vibration absorbers.
(a)
(b)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. Table1 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 k
p also affects the parameters f, ζ, and the forcing frequency factor g (see (8)). In the following, the masses adopted for the DVA were the same as those used for the nDVA.Table 1
Parameters used for the design of the linear DVA.
m
p [kg]
m
a [kg]
k
p [N/m]
F
0 [N]
Lower∗
Nominal
Upper∗
Lower∗
Nominal
Upper∗
100
2
4500
5000
5500
0.9
1
1.1
∗
Lower stands for lower limit. Upper stands for upper limit.Table 2
Parameters used for the design of the nDVA.
Parameters
ξ
p
ξ
a
ζ
p
ζ
a
µ
B
ρ
Nominal
0.0001
0.01
0.01
0.01
0.05
0.1
1.0
Lower∗
—
0.007
—
—
0.035
0.07
0.7
Upper∗
—
0.013
—
—
0.065
0.13
1.3
∗
Lower stands for lower limit. Upper stands for upper limit.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.3 Hz), comprising 300 frequency steps. Regarding the solution of the algebraic system defined by (9), the initial conditions of the nDVA were considered as follows: u
p = 5 m, u
a = 5 m, v
p = 2 m/s, and v
a = 4 m/s. The DVA was considered at rest for t
=
0.Table3 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 3
Optimization results obtained by the DVA design problem.
Approach
k
a [N/m]
c
a [Ns/m]
Cost function
Robustness
Deterministic
96.2744
1.9550
0.0000
24.2396
Dual optm.
96.2371
2.0474
0.0020
23.4711
Amplitude min
94.6159
2.1958
1.7131
22.0145
Fuzzy response
95.4366
2.1794
0.8235
22.2918Figure 9
Robustness envelopes determined in the design of the DVA.Table4 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.Table 4
Optimization results obtained by the nDVA design problem.
Approach
ξ
a
µ
β
ρ
F
1
F
2
Robustness
Deterministic
0.0092
0.0550
0.090
1.1000
8.2487
0.1043
245.90
Dual optm.
0.0078
0.0402
0.0780
1.1982
3.7370
0.1365
172.87
Amplitude min
0.0070
0.0350
0.0700
1.2993
3.4338
0.1485
142.42
Fuzzy response
0.0071
0.0350
0.0700
1.2996
3.4345
0.1485
143.51As 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 ofF
1 and maximization of F
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.Figure 10
Robustness envelopes determined in the design of the DVA.


## 6. 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 absorbers. The originality regarding the robust optimization methodology is related to the formulation and use of fuzzy metrics to evaluate the design robustness, handling the robust optimization problem as a monoobjective optimization instead of imposing a new cost function to the design problem and thus leading to a multiobjective optimization problem.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.



---

*Source: 2081518-2016-01-05.xml*