Structural Reliability Sensitivities under Nonstationary Random Vibrations

Response sensitivity evaluation is an important element in reliability evaluation and design optimization of structural systems. It has been widely studied under static and dynamic forcing conditions with deterministic input data. In this paper, structural response and reliability sensitivities are determined by means of the time domain covariance analysis in both classically and nonclassically damped linear structural systems. A time integration scheme is proposed for covariance sensitivity. A modulated, filtered, white noise input process is adopted to model the stochastic nonstationary loads. The method allows for the evaluation of sensitivity statistics of different quantities of dynamic responsewith respect to structural parameters. Finally, numerical examples are presented regarding a multistorey shear frame building.


Introduction
The problem of structural response to dynamic nondeterministic actions has to be posed when structures, such as rotating machinery components or constructions subject to sea/wind loads, are designed to operate in random vibration environments.In these cases, due to the intrinsic random nature of inputs, a random dynamic analysis is needed, as documented in random vibration theory (e.g., [1][2][3][4][5]).Random dynamic analysis offers valid information about structural response otherwise impossible to obtain through classic deterministic analysis.A suitable evaluation of structural integrity, fundamental to design, can be obtained through this approach, in which structural safety is evaluated by the probability that the analyzed structure will not collapse during its working or life time.More specifically, this probability regards one or more components of a given structure and the critical state or conditions of failure these components reach in a given time interval.In this work, failure is related to a first crossing threshold, so that the structural crisis is associated with a first overcoming of one or more structural response measures (from a safe domain).Typically, structural response measures are displacements, stresses, buckling loads or natural frequencies.Moreover, statistics of response and structural reliability including their sensitivities (i.e., evaluation of partial derivatives of a performance measure with respect to system parameters), allows us to evaluate structural performances when some design variables have parametric changes.In particular, sensitivity analysis is a suitable way to evaluate the response variation of structures and, due to its importance in structural design, it is actually an important research topic in the field of computational mechanics.The sensitivities of structural responses are a fundamental information for gradient-based optimization methods needed in structural reliability analysis, optimization, and identification.Even random structural analysis with uncertain system parameters, obtained by removing commonly used deterministic assumptions (often extremely unrealistic), needs structural sensitivities.The analysis and design sensitivity computations do not offer quite hard difficulties for static problems, and analytical formulations for design sensitivity analysis are being reported in the literature from early seventies [6][7][8].In dynamic problems, sensitivity analysis of structural response under transient dynamic loads has been studied especially in case of deterministic dynamic loads.Conversely not much research [9][10][11] seems to have been carried out for computation of design sensitivity under stochastic dynamic environment.Recently, an interesting analytical approach has been proposed by Chaudhuri and Chakraborty [12].It is based on response sensitivity evaluation of structures in seismic reliability evaluation.The formulation has been developed in double frequency domain to take nonstationary earthquake motion for obtaining the analytical sensitivity statistics of various dynamic response quantities with respect to structural parameters.A semianalytical approach to evaluate the sensitivity of the stochastic response of both classically and nonclassically damped structural systems under stationary or nonstationary stochastic Gaussian excitation has been proposed [13].The equations governing the evolution of the sensitivity of the first two statistical moments of the response for multi-degree-of-freedom (MDoF) structural systems are first derived in analytical form in the time domain, once the Kronecker algebra, the series, expansions, and a suitable modal expansion of the deterministic response are adopted.
Finally in [14] a new closed-form analytical approximation to the first-passage problem is presented in structural reliability by using the exact closed-form solutions for the spectral characteristics of nonstationary random processes.
The main objective of the present work is to propose an analytical method to evaluate in time domain response covariance and reliability sensitivity for general linear structures subjected to nonstationary stochastic loads.This analytical approach presents a serious reduction in computational cost, especially if compared with numerical procedures.The evolutive random input is modelled by a nonstationary modulated filtered white noise process, able to represent many real physic loads, as the earthquake loads here considered.In order to obtain a general nonstationary approach to be used in different contexts, the structural response is evaluated by a covariance approach.By knowledge of evolutionary covariance matrix in state-space, the reliability is derived with reference to first crossing failure events.Sensitivity of covariance matrix and reliability is then obtained.Therefore, the proposed formulation computes the sensitivity gradient of covariance matrix and structural reliability with respect to the structural design parameters.A time integration algorithm for solving this problem is proposed, and finally a multistory building idealized by a shear frame structure (Figure 2) is considered to study the sensitivity of dynamic responses.The rigid floors masses and interstorey lateral stiffness and damping are considered as design variables.

Linear Elastic MDoF Subject to Nonstationary Random Vibration
Many real cases of structural problem deal with structural configurations that could be represented by linear viscoelastic lumped masses system, subject to stationary or nonstationary actions.In this work, both of these input conditions are considered to evaluate structural response statistics by means of the covariance approach.This offers serious advantages both in stationary and more in nonstationary conditions if inputs are modelled as white noises, eventually filtered to increase their accordance with real dynamic loads.
The structural response of a deterministic second order   -degree of freedom linear mechanical system subject to a probabilistic dynamic input is determined by solving the dynamic equilibrium system equations: where   , Ẋ , and Ẍ are the structural displacement, velocity, and acceleration process vectors, and M, C, and K are the mass, viscous and stiffness matrices.
The forcing vector collects  stochastic excitations applied to the structure, and G s is a    matrix coupling the excitation components of forcing vector to the structural degrees of freedom.If system excitations vector elements are white noises, the first and second order statistical moments are time invariant, so that where ⟨⋅⟩ is for the mathematical expectation and (⋅) is Dirac's delta function.
If    () are Gaussian excitations, the responses   () and their time derivatives Ẋ () constitute a Markov vector in the 2  dimension phase state.
The matrix J, related to a vector satisfying the shot noise properties (usually denoted as a shot noise vector), has diagonal elements [J]  = 2  0 equal to the auto covariance intensity of each force and extra diagonal elements equal to , representing the level of correlation between two generic different forces, so that [J]  can vary from 2 0 , if   and   are completely correlated, to zero, if   and   are completely uncorrelated.
Then, in case of complete uncorrelated forcing loads, the matrix J is replaced by the more simple diagonal matrix of components [D]  = 2  0 , where the elements   0 are the input power spectral density of each entry.
A widely adopted approach consists in representing a nonstationary input by means of an intensity modulation of a stationary process, also denoted as uniform modulation, in which it is assumed that process intensity changes during time as a deterministic function (), without any variation in spectral contents.Thus, in case of time modulation, a stationary forcing process vector is replaced by the following nonstationary vector: with the stochastic characterization and, in case of uncorrelated excitations, the covariance matrix 2.1.Prefilters Technique.The prefilter approach is adopted when the autocorrelation function of the input is not known.In this approach, input processes are modelled by solving filter differential equations, whose input is a white noise process.So doing, it is possible to give a representation of physical phenomena whose frequency characterization is subject to seriously variations in time, potentially causing structural resonant consequences, while maintaining the advantages offered by a shot noise input.
In particular, in the prefilter approach, the filter response is described by the 2  filter state-space vector   = (  , Ẋ )  , solution of the 2  set of differential equations that generally could have a time-dependent form, when not only the frequency, but also the amplitude of loads, has an intrinsic evolutive nature.In (8),   () is a vector of   white noise processes (stationary or nonstationary), H 1  () and H 2  () are two     matrices, and G  is a 2    matrix that couples the excitation components of the forcing vector to the filter with a degree of freedom.Then, adopting the prefilter technique, the following system of motion differential equations is written in the state-space: which can be written also as ) and summarized as where is a new global 2(  +   ) state-space vector (structure plus filter).
In a design sensitivity problem, both the structural matrix and the response vector depend on a design parameter vector , whose elements may be structural stiffnesses, dampings, or masses, or any other mechanical parameters, as crosssections, Young modulus, boundary conditions, and so forth.Generally, also filter parameters and input intensity can be taken into account as design parameters.Thus, the matrix equation ( 11) must be explicitly re-written as: In case of zero initial condition, the solution of (13) has the general expression: where the matrix Φ(,  1 ,  2 ) (e.g., see [1]) is usually called a transition matrix.The mean space-state vector   () can be determined by the differential vectorial equation: and the covariance matrix R  (,  1 ,  2 ) of the state-space (, ) is

Mathematical Problems in Engineering
This second-order symmetric statistical moments matrix R  (,  1 ,  2 ) has (2 2 + ) independent elements and can be evaluated by the well-known Lyapunov differential matrix equation where and where P can be written as The matrix B(, ) can be simplified when the forcing vector is a white noise process as defined in (4) or in (7).In these cases, due to the properties of the Dirac function and of the transition matrix, the matrix P() is equal to Finally, B(, ) can be written as where the     submatrix L() is a diagonal matrix whose generic elements are For some applications, covariance information about structural acceleration is needed so that the matrix must be determined.This matrix is easily obtainable by the relation where

Measure of Structural Reliability
As stated previously, a wide class of structural dynamic problems deal with reliability evaluation.Due to the essential random nature of actions, a random dynamic analysis is necessary and the natural way of testing the structural integrity is to evaluate the probability whether a structure may not have a failure during its lifetime.Mechanical safety or reliability  at a fixed time  is then defined as a failure survival probability, where the failure is a partial or total static damage in a given time [0, ].It is then clear that the definition of failure plays a central role in the evaluation of reliability, given the many possible interpretations.Normally, only two different kinds of mechanical failure are considered: fatigue failure, due to cumulative damage, and threshold crossing failure, determined by the first time crossing of a structural response parameter  (able to describe failure conditions) through a given threshold value .The present work deals with threshold crossing failure.Regarding a generic mechanical system subject to a stochastic process, the reliability () is defined as the survival probability for the exceeding of the given threshold value  and, under the assumption that the initial survival probability is equal to one ((0) = 1), is given by The hazard function ℎ(, ) is defined as the probability of having a threshold crossing, in a unit time, in absence of previous threshold crossings.Its exact formulation is still an open question.In the Rice original formulation [15] for a single degree of freedom system, the hazard function has the following general form: in which   Ṡ (, ṡ | ()) is the joint probability density of () and Ṡ () processes, and () is the condition of absence of excursions of the fixed barrier  before the time .Difficulties related to the determination of this joint probability often impose the use of approximate solutions and one of which commonly used is the replacing of the conditional probability of failure with the unconditional probability   Ṡ (, ṡ ).By means of this assumption, the hazard function ℎ(, ) is replaced by the threshold crossing rate V +  (, ) which, under the further assumption that the above stochastic processes are Gaussian with null mean values, can be written as (1) (, )  (2) (, ) ×  (3) (, , )  [  (, , )] , (28) where  (1) (, ) =  Ṡ (, )   (, ) , ) In (30),   (, ) and  Ṡ (, ) are the standard deviations of (, ) and Ṡ (, ), and   Ṡ (, ) is their correlation factor.For a structural system, the threshold crossing rate V + s of a global safety domain  ⊂   can be determined by means of the following extension of the Rice formula [16], regarding a generic multidimensional stochastic process, and here referred to as    1 internal force vector s, whose elements are suitable structural response parameters (internal forces or stresses evaluated in a finite number of critical sections or points): where Γ  is a safety domain boundary and  is the outer normal to the boundary.Though an exact evaluation of reliability based on (31) is conceptually possible, it presents some difficulties.Many different simplified approaches have been proposed in last twenty years, when outcrossing rates of a vector random process have been studied in the context of problems in load combinations and in structural reliability [17][18][19][20][21][22][23].
In this work, an approximate and conservative evaluation of V + s is preferred as it can be easily adopted when the safe domain  is a hyper rectangle in   (Figure 1) and V + s is evaluated as the sum of the mean crossing rate on each face of the hyper rectangle .In this simplified approach, it is assumed that V + s is the sum of the crossing rates through each hyper plane on which the faces of  lie.For instance, in the two-dimensional space  2 , the safe domain  is the rectangle  1 ∩  2 , as shown in Figure 1, and ( 32) is an integral evaluated on the rectangular boundary Γ  .The safe domain boundary Γ  is replaced by Γ  1 ∪ Γ  2 , and V + s is evaluated as the sum of the crossing rates V + 1 and V + 2 out of the regions  1 and  2 , respectively.Generally, the safe domain can be obtained as where Since the integrand function in (31) is always nonnegative and given we have that where and in this particular case with symmetric barriers, we have In the following, V + s will always be evaluated as ∑ V +  .It must be noted that ( 35) is strictly satisfied in case of rare events.
The final structural reliability is evaluated as the product of all reliability evaluated separately for each structtural response parameter.Using this approach, it is possible to define a reliability vector () as the collection of independent reliability  ℎ () related to the first time crossing of a structural response parameter  ℎ through given symmetric threshold values ± ℎ that induce a structural failure.
It is then necessarily required to evaluate the probability  ,ℎ that the absolute value | ℎ | exceeds  ℎ , at least once during all input duration.The generic reliability vector element  ℎ () is defined as where, under the Poisson hypothesis, we have: where , and   Ṡ ℎ are, respectively, the ℎ and the  + ℎ diagonal elements of R     (), where   = (  Ṡ )  and   ℎ =  ℎ /  ℎ .With reference to a vector of response functions of interest () in reliability evaluation, related to the state space vector () by the linear transformation () = T 1 ()(), Ṡ () = T 2 ()(), the general form to express it become where and the covariance matrix in structural response parameters space-state   is The final global structural reliability is then The equation represents the probability that the ℎth interfloor displacement will cross the maximum acceptable value  ℎ during time interval [0, ].All the barriers can be collected in the barrier vector .
In case of failure due to a double symmetric threshold crossing, the threshold crossing rate V +  in the generic th critical node is given by where V +  (  , ) is determined by (37).

Response Reliability Sensitivity.
After covariance matrix R  has been determinated, its first order sensitivity, whose determination is necessary in determining the structural reliability sensitivity, is evaluated (in the following, the notations ⋅ ,  = {⋅} ,  =  ⋅ /  will be adopted).Consider In case of nonstationary conditions, the design sensitivity R  ,   of the covariance matrix R  is determined by solving the following matrix differential equation: The 2(  +   )2(  +   ) matrix A ,  has the following expressions which, without loss of generalities, can been simplified as: in which it has been assumed that input parameters do not depend on the design vector .In (47), we have To determine reliability sensitivity in structural response parameter space , it is necessary to obtain and thus the generic th reliability sensitivity is where the mean rate crossing failure sensitivity is (see details in Appendix B) In view of (42), the first order sensitivity of the global reliability is then (, )) .(52)

Numerical Example
In the present application of the proposed method, a multistorey building subject to an earthquake excitation is taken into consideration.For simplicity, and without loss of the generalities in the proposed technique, the building model adopted is a generic shear type plane frame structure, as shown in Figure 2. When floor slabs in a building generally have a very large degree of stiffness in-plane, then they can be regarded as rigid diaphragms.This increases analysis efficiency greatly without any significant loss of accuracy in the assessment of the response resulting from ground excitations.In addition to the rigid floor diaphragm assumption, the efficiency of computation can be further increased by applying the matrix condensation technique.
The main building mechanical model is assumed to be linear.This is an acceptable hypothesis within the horizontal displacement limitations required to guarantee an operative service level.In other words, this simplified hypothesis is perfectly functional when the assigned maximum interstorey drift is quite close to the elastic limit of the structural displacement (e.g., in case of full operatives demand).The design vector is assumed to be in order composed by each floor mass, stiffness, and damping: where The nonstationary character of the problem is introduced by the deterministic temporal modulation function (), which controls the intensity variation, without changing earthquake frequency contents.In this case, the modulation function proposed by Jennings [25] is adopted as follows: (57) The system of motion equations for the complete structural system then are It can be written in matrix form as where  = [1, 1, 1, . . ., 1, 1]  is the   1 drag vector, and M(), C(), and K() are, respectively, the mass, viscosity, and stiffness     structural matrices, whose general expression is The three vectors Ẍ (, ), Ẋ (, ), and   (, ) are acceleration, velocity, and displacement   1 vectors, relative to the ground.In the analyzed structure, the mass matrix is diagonal and the two viscous and stiffens matrices are tridiagonal.The matrix differential equation ( 59) is also written as

Ż (𝑏, 𝑡) = A (𝑏) 𝑍 (𝑏, 𝑡) + 𝐹 (𝑡) .
(61)  Thus, assuming the Poissonian crossing hypothesis and the final global structural reliability is then determined as: Although exact analytical solutions for  global (, , ) are generally definable from the theoretical point of view, it is known that ( 64) is an approximate upper-bound of global reliability and provides a conservative estimation of  global (), as previous stated.Moreover, it can be efficiently used for design and predesign purposes in practical engineering problems.
In order to evaluate the reliability vector related to the interfloor relative displacement threshold crossing, the interstorey drift vector is introduced as and the storey displacement vector (, ), where the transform matrix T is bidiagonal and is independent from design vector ( ) .
The covariance matrix R     (, ) is then related to R  (, ) through the following relations: The reliability vector   previously defined can be evaluated as the collection of where

Uℎ
,   ℎ Uℎ , and   ℎ =  ℎ /  ℎ .The first two are, respectively, the ℎ and the  + ℎ main diagonal elements of R     (), the third is the ℎ,  + ℎ element of the same matrix, and finally  ℎ is the ℎth barrier.Equation (68) represents the probability that the ℎth interfloor displacement will cross the maximum acceptable value  ℎ during time interval [0, ].All the barriers are collected in the vector .In the example here assumed its elements are all equal to 3.0 cm for each floor, which corresponds to a lateral drift of 1.0% in case of 3 m interstorey.
In Figure 3 the structural interstorey covariances are shown, in terms of displacement (a) and velocity (b).Continuous lines are for   = 10 sec., dashed lines are for   = 20 sec., dash-dot lines are for   = 30 sec., and finally dotted lines are for   = 40 sec.Blue lines are for first level covariance, magenta lines are for second level, and red lines are for the 3 dh level.It must be noticed that structural responses reach a stationary level if   is longer than 10 sec.Moreover, for this specific structural configuration, covariance responses of the first level are the greatest; the second is only slightly smaller than the first and the third is about the half.
In Figure 4, reliabilities are shown and structural safeties are evaluated for each lateral interstorey drift threshold crossing, and global reliability is evaluated as approximate upper-bound global reliability.Also in this figure, line type and line colours are used as in Figure 3.Moreover, here black lines are for global system reliability.It should be noticed that the first interstorey drift failure is the more probable, the second one has just a smaller failure probability, and finally the third one has a quite negligible failure probability (i.e.,  3 = 1).
In Figure 5, first order sensitivities of lateral interstorey drifts evaluated by the proposed approach are shown, for different time durations and storey levels.Results are plotted by using a nondimensional measure, represented by the parameter Sensitivities are evaluated with regards to three sets of mechanical parameters: the three floor masses (in graphs a, d, and g), the three floor stiffness (in graphs b, e, and h), and finally the three floor dampings (in graphs c, f, and i).
With reference to the first set of mechanical parameters, it must be noticed that the increase of floor masses (blue for the first one, magenta for the second, and red for the third one) has the effect of increasing lateral drift covariance responses.This effect takes place in any analysed case.Moreover, it must be also noticed that the first interstorey drift sensitivity has about the same values with respect to each three structural masses, as it is extremely evident in Figure 5(a).On the contrary, the second interstorey drift Figure 5(b) is more sensible to  2 (magenta) and  3 (red), then to  1 (blue).Finally, the third interstorey drift Figure 5(c) is more sensible with respect to  3 than to the other two masses.This is due to the fact that in general, a change of masses in a position lower     Mathematical Problems in Engineering than that of the analyzed interstorey drift produces a limited response variation.
A different sensitivity of interstorey drifts could be noticed in analysing the effect of stiffness variations, as shown in Figures 6(a), 6(b), and 6(c).It is evident that at each stiffness variation, a reduction of the interiorly variance    of the   -level occurs only when the stiffness   is increased (blue for  1 , magenta for  2 , and red for  3 ).On the contrary, for each -level, the increase of a rigidity   ( ̸ = ) in levels different from those investigated has an opposite effect, by increasing the lateral displacement variance    .For example, in Figure 6(b), it should be noticed that   2 decreases when  2 increase (magenta line) and meanwhile grows up (with slow velocity) when  1 (blue) and  3 (red) are increased.
A further effect on    has the variation of interstorey dampings.It is evident in Figures 7(a), 7(b), and 7(c) that increasing each of them (blue for  1 , magenta for  2 , and red for  3 ) has the effect of inducing a reduction in all floors lateral displacement.This is a direct effect of the energy dissipation that is induced by increasing structural damping.Then, the growing up of dissipated energy induces in any case the reduction of elastic and kinetic one, and, therefore, of elastic deformations.With regards to this specific structural system, maximum sensitivities in (   ) AD ,  have about the same magnitude (0.4-1.0) in case of masses and stiffness variations, and a smaller one (0.20-0.25) for damping variations.
First order sensitivities of reliability associated with interstorey drift exceeding over a given barrier, evaluated by (64), are shown in Figures 8, 9 Keeping in mind that positive sensitivities indicate an increase of survival probability, it is evident that the above sensitivities have the same trend of lateral drift covariance shown in Figure 3.Besides this, it has to be observed that the reliability sensitivities of the first two floors have comparable values (Figures 8(a c), 9(c) and 10(c) ) is smaller.This is related to the circumstance that  3 practically has a unitary value, and so its variation due to changing in structural parameters is practically negligible.
Finally in Figure 11, the first order sensitivities of global reliability are reported.These are evaluated with respect to variations in structural masses in Figure 11

Conclusions
A new method for evaluating the gradient with respect to structural parameters of stochastic response and structural reliability of a multi-degree-of-freedom system subject to nonstationary random loads has been proposed.The random input is modelled by a nonstationary modulated filtered white noise process, able to represent many real physic loads.In order to obtain a general nonstationary approach to be used in different contexts, the structural response is evaluated by a covariance approach, and the time domain differential matrix equations governing the sensitivity of the stochastic response are derived.In particular, the Lyapunov equation is differentiated with respect to the assigned design variables, which could be assumed as structural or input parameters.Based on this result, the structural reliability is also determined.This is defined as the first passage failure probability, expressed by the standard nonstationary expression based on Poissonian first crossing events.Sensitivity of covariance and reliability are then obtained.The proposed formulation computes the sensitivity gradient of covariance matrix and structural reliability with respect to the structural design parameters.A time integration algorithm to solve this problem is proposed.
Finally, a multistorey building idealized by a shear frame structure is considered to study the sensitivity of dynamic responses to a systemic action.The rigid floors masses and interstorey lateral stiffness and damping are considered as design variables.The structural behavior after a modification of these significant parameters has been evaluated by means of non-dimensional sensitivity parameters of covariance response and structural reliability.

Figure 2 :
Figure 2: Mechanical scheme of analysed plane shear type.

Figure 3 :
Figure 3: Interstorey displacement ((a),   ) and velocities ((b),  V ) response of a 3-DoF (different colors for different levels) system subject to modulated filtered white noise with dissimilar time durations (different line-types).

Figure 4 :
Figure 4: Structural reliability of a 3-DoF system, evaluated as the probability of maximum interstorey drift exceeds over a given threshold of 3 cm.

Figure 5 :
Figure 5: Covariance displacement sensitivity of a 3-DoF subject to different seismic duration (different line-types, continuous for   = 10, lines for   = 20, dot-line for   = 30, and dotted for   = 40) with respect to structural masses (different line colours for each storey mass) evaluated at each interstorey drift ((a) for the first, (b) for the second, and (c) for the third).

Figure 6 :
Figure 6: Covariance displacement sensitivity of a 3-DoF subject to different seismic duration (different line-types, continuous for   = 10, lines for   = 20, dot-line for   = 30, and dotted for   = 40) with respect to structural stiffness (different line colours for each storey   ) evaluated at each interstorey drift ((a) for the first, (b) for the second, and (c) for the third).

Figure 7 :Figure 8 :
Figure 7: Covariance displacement sensitivity of a 3-DoF subject to different seismic duration (different line-types, continuous for   = 10, lines for   = 20, dot-line for   = 30, and dotted for   = 40) with respect to structural damping (different line colours for each storey   ) evaluated at each interstorey drift ((a) for the first, (b) for the second, and (c) for the third).

Figure 9 :Figure 10 :
Figure 9: Level reliability   sensitivity of a 3-DoF subject to different seismic duration (different line-types, continuous for   = 10, lines for   = 20, dot-line for   = 30, and dotted for   = 40) with respect to structural stiffness (different line colours for each storey   ) evaluated at each interstorey drift ((a) for the first, (b) for the second, and (c) for the third).

Figure 11 :
Figure 11: Global reliability sensitivity of a 3-DoF system subject to different seismic duration (different line-types, continuous for   = 10, lines for   = 20, dot-line for   = 30, and dotted for   = 40) with respect to structural masses (a), structural stiffness (b), and structural damping (c).
, and 10.Also in this case, they are evaluated for different time durations and storey levels, with the same legend structure of previous figures.Results are plotted by using the non-dimensional measure represented by(   ()) AD ,  =    ()        () .(70)As in Figures5, 6, and 7, also in this case, the sensitivities are evaluated with regards to the three floor masses, stiffnesses, and dampings.Variations are evaluated assuming that   =   in Figures8(a), 8(b), and 8(c),   =   in Figures 9(a), 9(b), and 9(c), and finally   =   in Figures 10(a), 10(b), and 10(c).Also in these figures Blue lines are for  = 1, magenta lines are for  = 2, and red lines are for  = 3.
), 9(a), and 10(a) and Figures 8(b), 9(b), and 10(b) ), while the sensitivity of  3 , regarding the top floor, (Figures 8( (a), stiffiness in Figure11(b), and dampings in Figure11(c).It should be noticed that the global reliability has almost the same behaviour of  1 in Figures8, 9, and 10, being it the smallest one compared with the other two.