A Vibration Reliability Analysis Method for the Uncertain Space Beam Structure

Considering that uncertainty is inherent and unavoidable in engineering practice and the available information about the uncertain parameters is always not sufficient, the paper tries to carry out the nonprobabilistic vibration reliability analysis so as to avoid resonance on uncertain structure with bounded parameters. The input uncertain-but-bounded parameters are treated as interval variables, and an interval model is adopted to describe bounded uncertainties. Then a theory of nonprobabilistic reliability is introduced, in which the dimensionless nonprobabilistic reliability index and system reliability index are defined. In order to investigate the resonance failure with reliability method, the resonance failure domains are stated according to the relationships between the natural frequencies and the excitation frequencies. Then the uncertain structure is modeled as a series system and a system reliability index is proposed to evaluate the safety of the structure.The paper also takes a frequency analysis on the uncertain space beam structure to get the resonance failure modes. A frequency analysis method based on the monotonicity discriminant of the frequency sensitivity is presented. Then an optimization algorithm is introduced to verify the validity of the former frequency analysis method. Two examples are provided to illustrate the effectiveness and feasibility of the presented method.


Introduction
Uncertainty is inherent and unavoidable in almost all engineering problems.In structural engineering, uncertainties may present in structural materials, geometrical properties, theoretical modeling, failure conditions, applied loads, and so forth [1].Reliability problems stem rightly from the existence of all these uncertainties.At present, in structural engineering, several mathematical theories are available for dealing with uncertainties, for example, probability theory, fuzzy set and evidence theory, interval analysis and convex-set models, and so forth.Then corresponding structural reliability method can be developed on the basis of these mathematical models.Considering that the method needs less sample data compared to the traditional probabilistic reliability method [2], the nonprobabilistic reliability method based on interval model has become one of the most widely used reliability methods in academic research and engineering in the past few years.Based on interval arithmetic, a nonprobabilistic measure and analysis methodology for structural reliability computation are presented by Shuxiang et al. [3].Elishakoff and Ren [4] modeled structural uncertainties by both the stochastic and interval methods to quantify the uncertainties in response quantities.An interval stress and interval strength interference model for structural interval reliability analysis was investigated by Yang and Sun [5].The investigation was based on a partition, which gave a full-scale description of comparison relations between interval strength and interval stress.A bilevel robust optimization model, combining interval extension of function algorithm and an order relation of interval number algorithm, was presented to investigate the effects of uncertainties on the objective and constraint functions for robust optimization by Sun et al. [6].A new method was developed for reliability analysis of the uncertain structures with both random and interval variables by Jiang et al. [7].And then a new nonprobabilistic structural reliability analysis method based on a multidimensional parallelepiped convex model was developed by Jiang et al. in the literature [8].The model was constructed using marginal intervals for the variables and correlation information between any two 2 Shock and Vibration variables.A nonprobabilistic robust reliability method for robust control of uncertain structures using static output feedback LQR approach was presented systematically by Guo and Li [9].And a nonprobabilistic reliability methodology was presented systematically by Guo and Lu in the literature [10] for analysis and reliability-based design optimization of structures with uncertain-but-bounded parameters.Undoubtedly, the interval nonprobabilistic reliability idea based on set models is creative and significant.Recently, some scholars utilize reliability method to avoid resonance of the structure.Zhang et al. [11] presented the probability of the resonance state of the structure which was regarded as a series system.Su et al. [12] carried out a vibration reliability analysis method for avoiding the resonance and studied the effects of random parameters on vibration reliability of the resonance structure.The reliability methods employed in the literature [11,12] are both traditional probabilistic reliability methods.Considering that the available information about the uncertain parameters may be not sufficient in engineering practice, the interval nonprobabilistic reliability method should be a good choice for resonance avoiding.
In order to define the resonance failure modes of the uncertain structure, a frequency analysis which aims to get the upper and lower bounds of the natural frequencies is needed.It is also known as the interval eigenvalue problem which has aroused the interest of many scholars.Rohn [13] studied the generalized interval eigenvalue problem and derived formulas for the interval eigenvalue of a symmetric interval matrix with an error matrix of rank one.Based on the invariance properties of the characteristic vector entries, Deif [14] presented a method to compute interval eigenvalues for the standard interval eigenvalue problem.Qiu et al. [15] extended Deif 's method to the generalized interval eigenvalue problem.By viewing the deviation amplitudes of the mass and stiffness matrices as perturbations around the nominal values of the interval matrix pair, a perturbation method for the solution of the generalized interval eigenvalues problem has been developed by Qiu et al. [16].Elishakoff [17] proposed a procedure for finding the range of eigenvalues due to uncertain elastic modulus and mass density by using the upper and lower stiffness and mass matrices.Gao [18] proposed the interval factor method to investigate the effects of geometrical and material interval uncertainties on the natural frequencies and mode shapes of truss structures.Modares et al. [19] proved that, in the presence of any physically allowable uncertainty in the structural stiffness, the solutions of two deterministic eigenvalue problems are sufficient to obtain the exact bounds of the system's fundamental frequencies without resort to any combinatorial solution procedures.Sofi et al. [20] presented an efficient procedure for the solution of the generalized interval eigenvalue problem arising from vibration analysis of linear undamped structures with uncertain-but-bounded parameters.
The aim of this paper is to propose an efficient method for the frequency analysis of the space beam structure and take an interval nonprobabilistic reliability analysis of avoiding resonance on the uncertain structure.The remainder of this paper is organized as follows.In Section 2, a series system vibration reliability based on the nonprobabilistic reliability, different from the traditional probabilistic reliability, is introduced.An analysis on the failure modes of resonance in uncertain structure is carried out in Section 3, and the system reliability index of the resonance structure is obtained.In Section 4, a frequency analysis on the space beam structure is carried out and an efficient frequency analysis method is presented.A space frame structure and a simplified GARTEUR plane model are used to demonstrate the presented methods in Section 5, and finally conclusions are given in Section 6.

Overview of the Nonprobabilistic Reliability Method Based on Interval Model
Interval model is one of the most frequently used models in uncertainty analysis of structures.The model describes the uncertain-but-bounded parameters as interval variables X = ( where and   ∈ [−1,1] denotes the standard interval variables corresponding to   .  is the mean value of   , and   is the deviation of it.

Measure of Nonprobabilistic Reliability.
In structural reliability analysis, a limit-state function (LSF)  = (X) is usually defined by the performance or failure mode of the structure.(X) is usually a function of the input variables X = ( 1 ,  2 , . . .,   ), with (X) = 0 being the limit-state surface (LSS) which separates the variable space into two domains: the failure domain Ω  = {X | (X) < 0} and the safe domain In the situations that all the variables X = ( 1 ,  2 , . . .,   ) involved in the reliability problem under consideration can be bounded in certain intervals and represented by interval variables, as done above, all the interval variables are transformed into the standard interval variables  X = ( X1 ,  X2 , . . .,  X ) by means of transformation in (1).Then the LSF  = (X) can be transformed into the following normalized form [10]: and the reliability index can be defined as follows [3]: where ‖ ⋅ ‖ ∞ denotes the infinity norm of vector.It can be seen that the nonprobabilistic reliability index  in (4) represents the minimum distance which is measured by the infinity norm from the origin to the LSS (or to the failure domain) in the infinite topological space built by all the standard interval variables.The reliability index  is illustrated schematically by the two-dimensional case in Figure 1.When  > 1, the structure keeps safety; when  = 1, the structure is in the critical state; and when 0 <  < 1, part of the support domain is in the failure domain, the structure is considered to be unsafe for that it can't complete expected function convincingly.The case when  ≤ 0 is not included in the definition; it will be meaningless from the point of view of the reliability.

Nonprobabilistic Reliability Index of a Structure System.
In reliability analysis of a structure system, multiple LSFs which are defined by failure modes of the structure are usually involved.Assuming that there exist  potential significant failure modes which need to be taken into account, the corresponding LSFs are represented by   =   (X) ( = 1, 2, . . ., ), where   < 0 indicates the th failure mode.The structural system is usually regarded as a series system consisting of significant failure modes, because the occurrence of any one of the failure modes may lead to the failure of the structural system.Suppose that the input variables X = ( 1 ,  2 , . . .,   ) are independent interval variables.Transform all the interval variables into the standard ones  X = ( X1 ,  X2 , . . .,  X ) by means of (1) and then substitute that into the LSFs.Then the standardized LSFs   =   ( X ) ( = 1, 2, . . ., ) can be obtained.The nonprobabilistic reliability index of the structural system can be given as follows [3]: where   is the nonprobabilistic reliability index corresponding to the th LSF.The reliability index  sys of the structure Figure 2: The nonprobabilistic reliability of a structure system.system is illustrated schematically by the two-dimensional case in Figure 2. It can be seen from Figure 2 that  sys is the shortest distance from the origin to the failure surface. sys corresponds to the most likely occurred failure mode.So it is reasonable to take  sys in (5) as the measure of the structure system.

Nonprobabilistic Reliability of the Uncertain Resonance Structure
As we all know, resonance occurs when excitation frequency is close to the natural frequency of a vibrating system.The significant phenomenon of resonance is that the vibration amplitude and energy sharply increase, which may lead to the failure of the structure.The failure is called resonance failure here.

State Function of the Failure Resonance.
Resonance failure is different from the first passage problems [21]; it is also an important failure problem in vibration analysis.When resonance occurs, the vibration amplitude increases obviously, especially when the damp of the structure is low.The failure of resonance has two cases; one of them is that the amplitude is beyond the threshold, and the other is that the amplitude does not exceed the threshold due to the existence of damp [11].
In the first case, the structure is surely failure; in the second case, the propagation of fatigue crack is much accelerated by the resonance [22]; the structure is also regarded as failure here.
According to the reliability theory, the state function of failure resonance can be expressed as where   denotes the frequency of the th excitation,   denotes the th natural frequency of the structure, and  represents the critical spacing between natural frequency and excitation frequency.When the spacing between natural frequency and excitation frequency is less than , the resonance occurs.To the authors' knowledge, the value of  is related to the damp of the structure.The value of  should be reasonably defined in accordance with the effect of the amplitude on the structure and may be obtained based on previous experience in engineering practice.In this paper, a hypothetical value 1.5 Hz is given to .

Failure Mode of Resonance in Uncertain
Structure.For uncertain structures with interval parameters, the natural frequencies are implicit functions of the interval variables.So the natural frequencies   surely are interval numbers.Assuming that the excitation frequencies   are also interval numbers, the failure mode can be expressed as follows according to the state function in (6): where   and   can be transformed into standard ones as follows by means of (1): and   denote the mean value of   and   ;   and   denote the deviation of   and   , respectively;   and   are standard interval number equal to [−1, 1].Substituting ( 8) into (7), the expression of failure mode can be expressed as where  , = 0 denotes the two LSSs.
The failure modes are different due to the different relationships between   and   .Six different relationships are illustrated schematically in Figure 3.
The failure modes corresponding to different cases in Figure 3 are illustrated in Figure 4.
From Figure 4, it can be seen that the failure domains of resonance are infinite strip areas.In Figures 4(a) and 4(b), the origin of the standard space is in the failure domain.The structure in these two cases is considered to be failure for the absolute value of the difference between the mean value   and   is less than the critical spacing  and the reliability index  , is zero.In Figures 4(c) and 4(d), the failure domain passes through the support domain on either side of the origin, and the reliability index  , can be measured by the minimum distance (measured by the infinity norm) from the origin to the failure domain.In Figures 4(e) and 4(f), the midcourt line of the failure domain is out of the support domain, and the reliability index  , can also be measured by the minimum distance from the origin to the failure domain.
From Figure 4, it can also be known that the reliability index  , is the distance from the origin to upper boundary of the failure domain when the mean value   >   , or the distance from the origin to lower boundary of the failure domain when the mean value   <   .The reliability index  , can be expressed as follows according to (4): Failure domain Support domain Safe domain where  denotes the vector (  ,   ) and  , = 0 denote the LSSs.

Nonprobabilistic Reliability of the Uncertain Structure with
Resonance Failure.Assuming that the uncertain structure has  natural frequencies, with  excitation frequencies loading on it, the whole structure is considered to be in failure states because the resonance might occur when the excitation frequency is in the vicinity of any one of the natural frequencies.Thus, the structure system, to which excitation frequencies and natural frequencies are applied to analyze failure state, is a series system [11].The system reliability index is presented as follows according to (5): sys = min { , } ,  = 1, 2, . . ., ;  = 1, 2, . . ., . (11) System reliability index  sys is the minimum of  , , and it represents the most dangerous case for the structure when resonance occurs.

Frequency Analysis of Uncertain Space
Beam Structure

The Rayleigh Quotient Based on Finite Element Analysis.
For a space beam structure with  degrees of freedom, the natural frequency can be obtained by solving the following generalized Rayleigh quotient [23]: where   denotes the square of the th order natural frequency;  denotes the number of degrees of freedom; K denotes the global stiffness matrix of the structure, which is a  order real symmetric positive semidefinite matrix; M denotes the global mass matrix of the structure, which is a  order real symmetric positive definite matrix;   denotes the dimensionless th order standard mode.Taking Euler-Bernoulli beam structure into consideration,  space beam elements with  nodes are obtained by carrying a finite element analysis on the structure.And the global stiffness matrix and global mass matrix of the structure can be expressed as follows: where ] 6×6 , Shock and Vibration 7 Substituting K and M in ( 13) into ( 12), the square of the th order natural frequency is deduced as follows: where are the upper bounds of the interval parameters, respectively, then an investigation on the frequency sensitivity of material parameters is expanded.
Take the partial derivatives with respect to the interval variables   ,   , and   in (16).Considering that the error in the computed frequency by the Rayleigh quotient is proportional to the square of the deviation of the standard modes from their exact values [23], the effect of the modification of the standard mode on the natural frequency is a second-order small quantity.So that the effect can be ignored when taking the partial derivatives.Then the frequency sensitivities can be obtained: Considering the situation that the geometric parameters are certain parameters, then   ,   ,   ,   , and   are all certain positive values.According to ( 16) and ( 18), the global mass matrix M is positive definite, and the component matrices are given in ( 14) and (15).Based on the knowledge of linear algebra, a sufficient and necessary condition of positive semidefinite for the real symmetric matrices is that all of the principal minors are nonnegative.Checked by calculations, in ( 14), the nonzero principal minors of K ele (1) , K ele (2) , K ele (3) , and K ele ( 4) are all nonnegative numbers; in (15), the nonzero principal minors of M ele (1) , M ele (3) , and M ele ( 4) are all positive numbers and the nonzero principal minors of M ele (2) are all nonnegative numbers; therefore, the component matrices , and M ele (4) turn out to be positive semidefinite.
As the monotonicity of frequency sensitivity is obtained, the maximums and minimums of the frequencies can be expressed as follows by utilizing the generalized Rayleigh quotient in (12): where

Optimization Algorithm to Calculate the Frequency Interval.
For uncertain structure with bounded parameters, the square of the th order natural frequency   in ( 12) is in a certain interval range.Then the purpose of frequency analysis is to identify the interval ranges of the natural frequencies.Generally, this interval problem could be solved by an optimization algorithm.The optimization model can be established as follows: min where    denotes the lower bound of   and    denotes the upper bound of it.The global stiffness matrix K(X) and global mass matrix M(X) are functions of the interval variables vector X.By solving the optimization problems, the upper and lower bounds of the natural frequencies can be obtained.
Theoretically, an approximate value of the exact solution could be obtained via solving the former optimization model in (20).The optimization model is applicable to various uncertain structures, so it is utilized here to verify the correctness of the former presented monotonicity discriminant method.
Particle Swarm Optimization (PSO) is a metaheuristic evolutionary optimization technique which can be directly applied in a continuous global space environment and was first proposed by Kennedy and Eberhart [24].PSO is quite popular in the swarm intelligence community due to its simplest algorithmic structure, less parameter use, and being free from gradient use of an objective function.The PSO algorithm starts with randomly initialized population.It works on the social behavior of particle to get the best solution by adjusting each individual's position which can respect the global best position of the whole population.Each individual is adjusting by changing the velocity according to its own experience and by observing the position of the other particles by use of ( 21) [25].Consider where  is the inertia weight in [0.4,0.9];  1 and  2 are the acceleration factors in [0, 2] and  1 and  2 are random numbers in [0, 1];  is the population size and  is the stated iteration times; V  () denotes the velocity of the th particle after the th iteration and   () denotes the position of the th particle after the th iteration;   best is the best position   of the th particle and  best is the global best position in the population.
In this paper, a standard PSO algorithm is applied to calculate the maximums and minimums of the natural frequencies.The flowchart of the procedure for realization of the optimization process of ( 20) is depicted in Figure 5.After iteration calculation, the extremums of the natural frequencies can be worked out.It also can be seen that  ×  times of finite element analysis need to be accomplished, which may reduce the efficiency of calculation.

A Space Frame
Structure.The space frame structure (shown in Figure 6) is considered in this example.The frame is composed of 12 space beams.The sectional area of each beam is 0.01 × 0.01 m 2 .The elasticity modulus , Poisson's  the same results.Considering that the PSO has nice global searching capability, the result of the PSO is close to the exact solution.So the validity of the sensitivity monotonicity method presented in this paper is verified.
In addition, the sensitivity monotonicity method just experienced two times of FEM analysis to obtain the result, but the optimization algorithm in Section 4.4 took 20 × 5 times of FEM analysis away (20 represents the population; 5 represents the evolution times).So the sensitivity monotonicity method has much higher efficiency.
The changes of the first-order natural frequency along with the interval variables of the frame structure are shown in Figure 7.It can be seen that the frequency increases with the increasing of elasticity modulus and decreases with the increasing of Poisson's ratio and density, which shows the validity of the monotonicity of frequency sensitivity.
Already known are the minimal distance  = 1.5 Hz, excitation frequency  ∈ [58, 62] Hz, and the range of natural frequencies shown in Table 1.Transform the frequencies into standard ones.Then the system reliability index  sys can be calculated out according to ( 9), (10), and (11).The value of the system reliability index is 1.1633.The index is the distance from the origin to the lower boundary of the failure domain which is formed by  and  2 .The value of the index indicates Shock and Vibration  that the space frame structure is avoiding resonance while the excitation frequency  is loading on the structure.

The GARTEUR Plane
Model.GARTEUR plane model was developed by the French Aerospace Research Institute in 1995.The model has the same characteristic with the actual aircraft such as high flexibilities and modal density [26].As shown in Figure 8(a), the model has a total length of 1.5 m and a wingspan of 2.0 m, the material is aluminum, and there is a damping layer attached to the wing.The density, elastic modulus, and Poisson ratio of the wing changed due to the existence of damping layer [27].The whole model is simplified to a space beam model, and there are a total of 67 nodes and 66 beam elements; the restricted model is shown in Figure 8(b).The elastic modulus of the connection element was set as interval variable   ∈ [100, 1000] GPa after examining the simplification of the connection between the wing and the fuselage.Considering the existence of the damping layer, the elastic modulus of the wing was set as interval variable   ∈ [67.5, 68.5] GPa, the density of the wing as interval variable   ∈ [2700, 2900] Kg/m 3 , and Poisson's ratio of the wing as interval variable   ∈ [0.29, 0.31].The gust frequency loaded on the wing is assumed as  ∈ [10,15] Hz.By use of the sensitivity monotonicity method and the optimization algorithm, the frequency analysis was, respectively, carried out on the GARTEUR plane model.The frequency intervals were, respectively, obtained by calculating.The first 6 order intervals are shown in Table 2. From Table 2 we can find that the results of the two methods are nearly the same.Considering that the result of the PSO is close to the exact solution, the validity of the sensitivity monotonicity method is verified.
In addition, the sensitivity monotonicity method only experienced two times of FEM analysis to obtain the results, but the optimization algorithm in Section 4.4 took 20 × 10 times of FEM analysis away (20 represents the population; 10 represents the evolution times).So the sensitivity monotonicity method is much more efficient.
The changes of the first-order natural frequency along with the interval variables of the GARTEUR plane model are shown in Figure 9.It can be seen that the frequency increases with the increasing of elasticity modulus of the wing and decreases with the increasing of density of the wing.The elasticity modulus of the connection and Poisson's ratio of the wing have little impact on the first-order natural frequency.It is caused by the special structure of the connection and the wing.The analysis result of Figure 9 also shows the validity of the monotonicity of frequency sensitivity.
Already known are the minimal distance  = 1.5 Hz, the gust frequency  ∈ [8,14] Hz, and the range of natural frequencies shown in Table 2. Transform the frequencies into standard ones.Then the system reliability index  sys can be calculated out according to (9), (10), and (11).The value of the system reliability index is 1.4123.The index is the distance from the origin to the upper boundary of the failure domain which is formed by  and  1 .This value of the index indicates that the GARTEUR plane model is avoiding resonance while the gust frequency is loading on the wing.

Conclusion
Considering the facts that available information about the uncertain parameters is often limited in practical engineering and it is relatively easy to construct a conservative interval model, a series system vibration reliability based on the interval nonprobabilistic reliability is introduced to avoid resonance of the structures with uncertain-but-bounded parameters in this paper.According to the interval theory, the failure domains similar to infinite strips are presented and the LSSs are obtained corresponding to the different relationship between the natural frequencies and the excitation 6.75 6.76 6.77 6.78 6.79 6.8 6.81 6.82 6.83 6.84 6.85 4.25 frequencies.Therefore, the series system vibration reliability index is obtained.Although this series system has the same failure mechanism with the traditional stochastic series system, it has a different system reliability index due to the different variable types.The series system vibration reliability index represents the shortest distance from the origin to the LSSs which are the boundaries of the failure domain.It also can be regarded as a robust reliability index because the structure is much more robust when the index is larger.In order to get the failure modes of the uncertain structure, a frequency analysis is needed.In this paper, the space beam structure with material parameters as the interval variables is considered as the object.A frequency analysis method based on the monotonicity discriminant of the frequency sensitivity is presented.This method only needs two times of FEM analysis to calculate the upper and lower bounds of the frequencies, so it is much more efficient.Then the global optimization algorithm PSO is introduced to verify the validity of the method.At last, a space frame structure and a simplified GARTEUR plane model are provided to illustrate the validity and feasibility of the presented method.

Figure 1 :
Figure 1: The nonprobabilistic reliability based on two-dimensional interval model.
velocities and positions of the particles: Yes No

Figure 5 :
Figure 5: Flowchart of the PSO procedure to calculate ranges of the natural frequencies.

FirstFigure 7 :
Figure 7: The change of the first-order natural frequency along with the variables of the frame.(a) The change along with elasticity modulus.(b) The change along with Poisson's ratio.(c) The change along with density.

FirstFigure 9 :
Figure 9: The change of the first-order natural frequency along with the variables of the GARTEUR model.(a) The change along with elasticity modulus of the connection.(b) The change along with elasticity modulus of the wing.(c) The change along with Poisson's ratio of the wing.(d) The change along with density of the wing.
where    is the lower bound of   and    is the upper bound of it.Then the interval variables can be transformed into standard ones by the following normalized transformation: 1 ,  2 , . . .,   ) satisfying   ∈ [   ,    ] ( = 1, 2, . . ., ), denotes the 6 × 6 global stiffness matrix; M denotes the 6×6 global mass matrix;   is the beam crosssection area of the th element;   is the  axial moment of inertia;   is the  axial moment of inertia;   is the polar moment of inertia; and   is the beam element length.For isotropic material,   =   /(2(1+  )) holds.  is the elasticity modulus of the th element;   is Poisson's ratio and   is the density; K ele  is the expanded 6×6 element stiffness matrix of the th element; M ele  is the expanded 6×6 element mass matrix of the th element; K ele gets the minimum value   min when   ,   get the upper bounds  up  ,  up  and simultaneously   gets the lower bound  low  .
max and   min are the th order standard modes when   gets the maximum and minimum values, respectively.  gets the maximum value   max when   ,   get the lower bounds  low  ,  low  and simultaneously   gets the upper bound  up  ;

Table 1 :
Natural intervals of the space frame structure.

Table 1 .
According to the sensitivity monotonicity method in Section 4, the natural frequencies get the maximum values when the elasticity modulus  is on the upper bound meanwhile the Poisson ratio  and density  are on the lower bounds, and the natural frequencies get the minimum values when the elasticity modulus  is on the lower bound meanwhile the Poisson ratio  and density  are on the upper bounds.From Table1, we can find that the sensitivity monotonicity method and the optimization algorithm have

Table 2 :
Natural frequency intervals of GARTEUR model.