Research Article | Open Access

# Fuzzy Stochastic Vibrations of Double-Beam Complex System as Model Sandwich Beam with Uncertain Parameters

**Academic Editor:**S. He

#### Abstract

The dynamic behavior of a double Euler-Bernoulli beam system with uncertain parameters (fuzzy random variables) under a fuzzy stochastic excitation and axial compression is being considered. The beams are identical and parallel, one is above the other, and they are continuously coupled by a linear two-parameter (Pasternak subsoil) elastic element. This double Euler-Bernoulli beam system can be also treated as a theoretical model of a sandwich beam. The load process is fuzzy random both in space and time. The top beam carries a fuzzy stochastic load. The solution of the problem was found thanks to the fuzzy random dynamic influence function. The aim of the paper is to find the solution for the membership function of the probabilistic characteristics of the response of the structure.

#### 1. Introduction

Uncertainty modeling in computational mechanics has received significant attention in recent years. Dynamic analysis of structures often involves two kinds of uncertainty. One of them is referred to as aleatory uncertainty (randomness, stochastic variability) and the other as epistemic uncertainty which describes, among others, imprecision, vagueness, ambiguity, and lack of the data. The random variability is described by the use of the probability theory and the imprecision by the use of fuzzy sets. Very often sufficient statistical data is not available; in this case a fuzzy function (fuzzy process) or a fuzzy random variable (fuzzy stochastic process) is suitable for the modeling purposes. The basis of the probabilistic methods is primarily statistical data while in the theory of fuzzy sets we refer also to the human intuition and experience. Using fuzzy random variables or fuzzy stochastic processes we combine these two approaches. The classical problem of the dynamic response of structures under stochastic excitation has been presented in many monographs. In most cases it is assumed that the parameters of the structure are deterministic. On the other hand, the structural parameters like geometry characteristics and material and damping properties might be uncertain to some extent. Their uncertainty may have a strong influence on the reliability of a structure in the dynamic context. The idea of fuzzy set theory was initiated by Zadeh [1]. The concept of fuzzy random variables was introduced by Kwakernaak [2] and Puri and Ralescu [3] and combines both randomness and imprecision. The application of the uncertain forecasting in engineering and computational mechanics based on fuzzy stochastic processes is presented in the monograph [4]. Various types of structures and girders like beams, plates, shells, and frames have been considered. An important technological extension of a single string, beam, or plate is a double-string, double-beam, or string-beam system. Various aspects of the dynamic response of a double-string and double-beam system have been considered by Oniszczuk [5–9]. Free and forced vibrations of a double-beam system have been considered among others in [10, 11]. The problem of the vibration and buckling of a double-beam system under compressive axial loading is presented in [12]. Vibrations of a complex system under moving force have been studied in [13–15].

In this paper the dynamic behavior of a double Euler-Bernoulli beam system with uncertain parameters (fuzzy random variables) under a fuzzy stochastic excitation and axial compression is being considered. The beams are identical, parallel, one above the other and continuously coupled by a linear two-parameter (Pasternak subsoil) elastic element. This double Euler-Bernoulli beam system can be also treated as a theoretical model of a sandwich beam, in which the shear and compression of core are taken into account. The load process is fuzzy random both in space and time. The solution of the problem was found thanks to the fuzzy random dynamic influence function [16–18]. The aim of the paper is to find the solution for the membership function of the probabilistic characteristics of the response of the structure. The probabilistic characteristics of the response of the structure are sought in the form of the first two probabilistic moments, that is, the expected value and the correlation (covariance) function. The stochastic response of structure with random and uncertain parameters has been considered, among others, in [18–27]. Fuzzy stochastic finite element method based a spectral approach to analyze complex engineering structures under dynamic excitation has been presented in [28].

#### 2. Formulation of the Problem: General Solution

Let us consider the vibrations of a connected double Euler-Bernoulli beam continuous system connected by two-parameters elastic element under axial compression excited as is shown Figure 1.

The differential equations of motion of the beams system have the form where is the flexural rigidity of the beam, is Young’s modulus of elasticity, is the moment of inertia of the cross-section area of the beam, is the mass of the beam, is the damping coefficient, is the elastic stiffness (Winkler’s parameter), and is the shear stiffness of the elastic element.

The boundary conditions of the simply supported beam are The above double-beam complex system can be treated also as a theoretical model of a sandwich beam which consists of two bendable outer layers and a sheared middle layer, which is also compressed.

The following system parameters are assumed to be uncertain: the Young’s modulus , the mass of the beam , the damping coefficient , and the parameters of the elastic element , . The symbol denotes an elementary event, and for simplicity we will skip it in the paper. The symbol denotes the fuzziness of the parameters. The uncertain parameters of the system are assumed to be fuzzy random variables. The load process is assumed to be a fuzzy stochastic process. The responses of the system and are also fuzzy stochastic processes as well. The aim is to find the solution for the membership function of the probabilistic characteristics of the response, and of the double-beam complex system. The probabilistic characteristics of the response of the structure are sought in the form of the first two probabilistic moments, that is, the expected value and the correlation (covariance) function.

In order to decouple (1) and (2), let us introduce two new functions

From (1) and (2), we obtain two new differential equations for functions and :

Equation (5) describes vibrations of a single beam, while (6) describes vibrations of a single beam resting on an elastic Pasternak support with parameters and . From (4), it follows that

In the case when the parameters of (5) and (6) are random, we can solve the problem if the right-hand side is deterministic. To overcome these difficulties, we introduce the space-time fuzzy random dynamic influence functions and (FRDIF) which satisfy the following equations with a deterministic right-hand side:

After solving (8) one obtains

If the FRDIF and are known, the response of the system and can be presented in the following form: where , and if , then one considers transition vibrations, and for , one considers steady-state vibration case.

Thus, in order to determine the probabilistic characteristics of the displacement of the beams, one can apply the expectancy operator to (10) and consequently obtain the expected values as follows: Using -level optimization procedure [4] for arbitrary or the max-min operator in the extension principle [1], the smallest and the largest expected values at an established point and time can be found.

Due to relationship (11), one obtains

Taking into account (10), the covariance functions of the beams displacement on the fuzziness level have the form where

The lower and upper endpoints of the covariance could be defined using (13) as

In order to find the expected values and , we can use the perturbation method or Monte Carlo method. In this paper the perturbation method has been used for a particular solution of the expected values and variances and is presented in Section 3 for the response of the double-beam system.

#### 3. Particular Solutions

The solutions of (8) for boundary conditions (3) are assumed to be in the form of the sine series: where .

By substituting expression (17) into (8) and using orthogonalization method, one obtains the following set of uncoupled ordinary differential equations: where, for , and for , and .

The dots denote differentiation with respect to the time. These functions fulfill the initial conditions The solution of (18) taking into account the initial conditions (19) has the form where .

Taking into account (9), (17), and (20) we have

Let us consider the steady-state vibration , assuming that the excitation load is a fuzzy weakly stationary “white noise” stochastic process both in time and space. In this case we have . and . We assume that the Young’s modulus is a fuzzy random variable. The other system parameters are assumed to be deterministic. In order to find the probabilistic characteristics the function of the random variables has been expanded into Taylor series around the mean value and restricted to three items of the expansion. The expected value for steady-state case is equal to where is the variance of the Young modulus.

Formula (13) gives the variance of the system displacement for the steady-state vibrations in the form where The numerical calculations have been done under the following assumptions.(i)The mean values and are fuzzy numbers with triangular membership function (Figures 2 and 3),(ii)The variation coefficients are constant and do not depend on the level of fuzziness.(iii)Other quantities describing the structure are deterministic.

**(a)**

**(b)**

**(a)**

**(b)**

The following values of the parameters are used in the following numerical calculations: , , , , , , , and , , ., , and .

In (23) it has been assumed that .

In Figures 2, 3, 4, 5, 6, and 7 the membership of the expected value and variance of the beams response are presented.

**(a)**

**(b)**

**(c)**

**(d)**

**(e)**

**(a)**

**(b)**

**(a)**

**(b)**

**(a)**

**(b)**

**(c)**

**(d)**

**(e)**

#### 4. Fuzzy Stochastic Moving Load

Let us notice that the general solution presented in Section 2 can also be used, after some modifications, if the double-beam system is loaded by a fuzzy stochastic load moving with velocity (Figure 8).

In this case the right-hand side in (1) is equal to .

The FRDIF () in solutions (10)–(16) should be replaced by fuzzy random dynamic moving influence function (FRDMIF) which satisfies, for , the equations

The fuzzy random dynamic moving influence functions have the forms where , , , and .

The general solutions for moving fuzzy stochastic load after modification of (10)–(16) have the forms where , , where

The lower and upper endpoints of the covariance could be defined using (32) as where and .

For steady-state solutions and it should be assumed that .

#### 5. Conclusion

In the paper, the dynamic behavior of a double Euler-Bernoulli beam complex system with uncertain parameters (fuzzy random variables) under a fuzzy stochastic excitation and axial compression has been studied. The load process is fuzzy random both in space and time. In order to find the solution for the membership function of the probabilistic characteristics of the response of the structure, the idea of the fuzzy random dynamic influence function has been used. The probabilistic characteristics of the response of the structure are sought in the form of the first two probabilistic moments, that is, the expected value and the correlation (covariance) function. This double Euler-Bernoulli beam system can be also treated as a theoretical model of a sandwich beam. The algorithm (similar to the one presented in the paper for the double Euler-Bernoulli beam complex system) can be used in stochastic dynamic analysis of other complex systems like strings, beams, plates, and so on, with uncertain parameters.

#### Acknowledgment

The second author gratefully acknowledges the support of the Polish National Science Centre under Grant number ID: 147573.

#### References

- L. A. Zadeh, “Fuzzy sets,”
*Information and Computation*, vol. 8, pp. 338–353, 1965. View at: Google Scholar | Zentralblatt MATH | MathSciNet - H. Kwakernaak, “Fuzzy random variables. I. Definitions and theorems,”
*Information Sciences*, vol. 15, no. 1, pp. 1–29, 1978. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - M. L. Puri and D. A. Ralescu, “Fuzzy random variables,”
*Journal of Mathematical Analysis and Applications*, vol. 114, no. 2, pp. 409–422, 1986. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - B. Möller and M. Beer,
*Fuzzy Randomness. Uncertainty in civil engineering and computational mechanics*, Springer, Berlin, Germany, 2004. View at: MathSciNet - Z. Oniszczuk, “Transverse vibrations of elastically connected double-string complex system, Part I: free vibrations,”
*Journal of Sound and Vibration*, vol. 232, no. 2, pp. 355–366, 2000. View at: Publisher Site | Google Scholar - Z. Oniszczuk, “Transverse vibrations of elastically connected double-string complex system, Part II: forced vibrations,”
*Journal of Sound and Vibration*, vol. 232, no. 2, pp. 367–386, 2000. View at: Publisher Site | Google Scholar - Z. Oniszczuk, “Free transverse vibrations of elastically connected simply supported double-beam complex system,”
*Journal of Sound and Vibration*, vol. 232, no. 2, pp. 387–403, 2000. View at: Publisher Site | Google Scholar - Z. Oniszczuk, “Damped vibration analysis of an elastically connected complex double-string system,”
*Journal of Sound and Vibration*, vol. 264, no. 2, pp. 253–271, 2003. View at: Publisher Site | Google Scholar - Z. Oniszczuk, “Forced transverse vibrations of an elastically connected complex simply supported double-beam system,”
*Journal of Sound and Vibration*, vol. 264, no. 2, pp. 273–286, 2003. View at: Publisher Site | Google Scholar - S. Kukla, “Free vibration of the system of two beams connected by many translational springs,”
*Journal of Sound and Vibration*, vol. 172, no. 1, pp. 130–135, 1994. View at: Publisher Site | Google Scholar - H. V. Vu, A. M. Ordóñez, and B. H. Karnopp, “Vibration of a double-beam system,”
*Journal of Sound and Vibration*, vol. 229, no. 4, pp. 807–822, 2000. View at: Publisher Site | Google Scholar - Y. Q. Zhang, Y. Lu, S. L. Wang, and X. Liu, “Vibration and buckling of a double-beam system under compressive axial loading,”
*Journal of Sound and Vibration*, vol. 318, no. 1-2, pp. 341–352, 2008. View at: Publisher Site | Google Scholar - W. Szcześniak, “Vibration of elastic sandwich and elastically connected double-beam system under moving loads,”
*Scientific Works of Warsaw University of Technology, Civil Engineering*, vol. 132, pp. 111–151, 1998 (Polish). View at: Google Scholar - M. Abu-Hilal, “Dynamic response of a double Euler-Bernoulli beam due to a moving constant load,”
*Journal of Sound and Vibration*, vol. 297, no. 3–5, pp. 477–491, 2006. View at: Publisher Site | Google Scholar - J. Rusin, P. Śniady, and P. Śniady, “Vibrations of double-string complex system under moving forces. Closed solutions,”
*Journal of Sound and Vibration*, vol. 330, no. 3, pp. 404–415, 2011. View at: Publisher Site | Google Scholar - G. Adomian,
*Stochastic systems*, vol. 169, Academic Press, Orlando, Fla, USA, 1983. View at: MathSciNet - P. Śniady, R. Adamowski, G. Kogut, and W. Zielichowski-Haber, “Spectral stochastic analysis of structures with uncertain parameters,”
*Probabilistic Engineering Mechanics*, vol. 23, no. 1, pp. 76–83, 2008. View at: Publisher Site | Google Scholar - K. Mazur-Śniady, P. Śniady, and W. Zielichowski-Haber, “Dynamic response of micro-periodic composite rods with uncertain parameters under moving random load,”
*Journal of Sound and Vibration*, vol. 320, no. 1-2, pp. 273–288, 2009. View at: Publisher Site | Google Scholar - R. G. Ghanem and P. D. Spanos,
*Stochastic Finite Elements: A Spectral Approach*, Springer, New York, NY, USA, 1991. View at: Publisher Site | MathSciNet - W. K. Liu, A. Mani, and T. Belytschko, “Finite element methods in probabilistic mechanics,”
*Probabilistic Engineering Mechanics*, vol. 2, no. 4, pp. 201–213, 1987. View at: Google Scholar - J. Li and S. Liao, “Response analysis of stochastic parameter structures under non-stationary random excitation,”
*Computational Mechanics*, vol. 27, no. 1, pp. 61–68, 2001. View at: Publisher Site | Google Scholar - W. Mironowicz and P. Śniady, “Dynamics of machine foundations with random parameters,”
*Journal of Sound and Vibration*, vol. 112, no. 1, pp. 23–30, 1987. View at: Google Scholar - H. Jensen and W. D. Iwan, “Response of systems with uncertain parameters to stochastic excitation,”
*Journal of Engineering Mechanics*, vol. 118, no. 5, pp. 1012–1025, 1992. View at: Google Scholar - I. D. Gupta and R. G. Joshi, “Response spectrum superposition for structures with uncertain properties,”
*Journal of Engineering Mechanics*, vol. 127, no. 3, pp. 233–241, 2001. View at: Publisher Site | Google Scholar - F. Yamazaki, M. Shinozuka, and G. Dasgupta, “Neumann expansion for stochastic finite element analysis,”
*Journal of Engineering Mechanics*, vol. 114, no. 8, pp. 1335–1354, 1988. View at: Google Scholar - C. Soize, “Vibration damping in low-frequency range due to structural complexity: a model based on the theory of fuzzy structures and model parameters estimation,”
*Computers & Structures*, vol. 58, no. 5, pp. 901–915, 1996. View at: Publisher Site | Google Scholar | MathSciNet - M. Gładysz and P. Śniady, “Spectral density of the bridge beam response with uncertain parameters under a random train of moving forces,”
*Archives of Civil and Mechanical Engineering*, vol. 9, no. 3, pp. 31–47, 2009. View at: Google Scholar - P. Sniady, K. Mazur-Sniady, R. Sieniawska, and S. Zukowski, “Fuzzy stochastic elements method: spectral approach,”
*Mechanical Systems and Signal Processing*, vol. 37, pp. 152–162, 2013. View at: Google Scholar

#### Copyright

Copyright © 2013 Krystyna Mazur-Śniady et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.