Abstract
Based on the Rayleigh beam theory, the forced transverse vibrations of an elastically connected simply supported doublebeam system with a Pasternak middle layer subjected to compressive axial load are investigated. It is assumed that the two beams of the system are continuously joined by a Pasternak layer. The dynamic responses of the system caused by arbitrarily distributed continuous loads are obtained. The effect of Pasternak layer on the forced vibrations of the Rayleigh doublebeam system are discussed for one particular case of excitation loading. The properties of the forced transverse vibrations of the system are found to be significantly dependent on the compressive axial load and shear foundation modulus of Pasternak layer. Vibrations caused by the harmonic exciting forces are discussed, and conditions of resonance and dynamic vibration absorption are formulated. The important result on which this paper puts emphasis is that the magnitudes of the steadystate vibration amplitudes become smaller when the shear Pasternak modulus increases and Pasternak layer can reduce the magnitudes of the steadystate vibration amplitudes more than a Winkler elastic layer. Thus the Rayleigh beamtype dynamic absorber with a Pasternak layer can be used to suppress the excessive vibrations of corresponding beam systems instead of those with a Winkler elastic layer
1. Introduction
A great number of mechanical systems are complex structures composed of two or more basic mechanical systems whose dynamic behaviors are conditioned by their interactions. The systems connected by an elastic layer constitute one group of such mechanical structures which are commonly encountered in mechanical systems, construction and aeronautical industry. Such mechanical systems are exposed to an increased possibility of creating resonance conditions which can cause breakage and damage. These models are important as they give the initial approximation of the solution and a general insight into a dynamic behavior of the system at slight motion. The problem concerning the vibrations of beams joined by a Winkler elastic layer has attracted the interest of a large group of scientists. The problem of two elastically connected beams jointed by the Winkler elastic layer has been investigated in order to determine the conditions for the behavior of the system acting as a dynamic absorber in technical practice.
A mathematical model was developed by Seelig and Hoppman [1]. They investigated the problem of an impulse load effect on a beam and generated a system of partial differential equations describing its vibration. The obtained theoretical and experimental results which confirmed a sound approximation of an analytical solution obtained for slender beams at small transverse motions using the EulerBernoulli theory.
Oniszczuk [2, 3] analyzed the problem of free and forced vibration of two elastically connected EulerBernoulli beams. He determined analytical solutions for amplitudes functions and vibration modes. He discussed the effect of the stiffness of the elastic interlayer on the frequencies and amplitudes of the system. He determined the conditions for the occurrence of resonance and the behavior of the system as a dynamic absorber.
The analysis of the system composed of two connected beams was carried on by Zhang et al. [4, 5]. They presented analytical solutions for natural frequencies of the system as a function of axial compression force impact and their effect on the vibration amplitude. They determined the codependency between the system’s critical force and the Euler critical load as a function of an axial force of the other beam.
The behavior of foundation materials in engineering practice cannot be represented by foundation model which consists of independent linear elastic springs. In order to find a physically close and mathematically simple foundation model, Pasternak proposed a socalled twoparameter foundation model with shear interactions. In an attempt to find a physically close and mathematically simple representation of an elastic foundation for these materials, Pasternak [6] proposed a foundation model consisting of a Winkler foundation with shear interactions. This may be accomplished by connecting the ends of the vertical springs to a beam consisting of incompressible vertical elements, which deforms only by transverse shear.
In Wang and Stephens [7] and De Rosa [8], the natural vibration of a Timoshenko beam on a Pasternaktype foundation is studied. Frequency equations are derived for beams with different end restraints. A specific example is given to show the effects of rotary inertia, shear deformation, and foundation constants on the natural frequencies of the beam.
Stojanovic et al. [9] analyzed free vibration and static stability of two elastically connected beams with the Winkler elastic layer inbetween considering the effects of rotary inertia and transverse shear, which leads to the rotation of the crosssection. They used an example to show the analytical solutions for natural frequencies and determined the critical force for the coupled beam system. Stojanovic et al. [10] discussed the case of forced vibration of two elastically connected beams with a Winkler elastic layer inbetween and the effect of an axial compression force on the amplitude ratio of system vibration for three types of external forcing (arbitrarily continuous harmonic excitation, uniformly continuous harmonic excitation and concentrated harmonic excitation). They determined general conditions of resonance and dynamic vibration absorption. Stojanovic et al. [11] discussed the analytical analysis of static stability of a system consisting of three elastically connected Timoshenko beams on an elastic foundation. They provided expressions for the critical force of the system under the influence of elastic Winkler layers. Stojanovic et al. [12] using the example of multiple elastically connected Timoshenko and ReddyBickford beams, determined the analytical forms of natural frequencies, their change under the effect of axial compression forces and the conditions for static stability for a different number of connected beams.
Li and Hua [13], using the spectral finite element method arrived at numerical solutions for natural frequencies of the two elastically connected Timoshenko beams with different supports. They determined vibration modes and amplitudefrequency dependence for forced vibration of the system. Also, Li et al. [14] addressed the problem of three elastically connected Timoshenko beams in their paper. They numerically determined natural frequencies for different types of support, mode shapes and the effect the stiffness of elastic interlayers had on vibration of the system.
Ariaei et al. [15] studied the problem of a movable body according to the system of multiple elastically connected Timoshenko beams on an elastic foundation using the Transfermatrix method. They concluded that the maximum deflection of the beam system decreases in case of a body moving on the beams closer to the surface.
Mao [16] determined mode shapes for the first ten system modes of multiple elastically connected Euler beams on an elastic foundation using the Adomian decomposition method.
Kozic et al. investigated the Lyapunov exponent and moment Lyapunov exponents of two degreesoffreedom linear systems subjected to white noise parametric excitation [17]. The method of regular perturbation used to determine the explicit asymptotic expressions for these exponents in the presence of small intensity noises. The Lyapunov exponent and moment Lyapunov exponents are important characteristics for determining the almostsure and moment stability of a stochastic dynamical system. As an example, they study the almostsure and moment stability of the doublebeam system under stochastic compressive axial loading. The validity of the approximate results for moment Lyapunov exponents checked by the numerical Monte Carlo simulation method for this stochastic system.
Stojanovic et al. studied the stochastic stability of the three elastically connected Euler beams on elastic foundation [18]. Stochastic stability conditions expressed by the Lyapunov exponent and moment Lyapunov exponents. The method of regular perturbation used to determine the asymptotic expressions for these exponents in the presence of small intensity noises. Analytical results presented for the almost sure and moment stability of a stochastic dynamical system. The results applied to study the moment stability of the complex structure with influence of the white noise excitation due to the axial compressive stochastic load. Also, Stojanovic et al. suggested a new $\rho $version finite element to deal with geometrically nonlinear vibrations of damaged Timoshenko beams [19]. The novelty of the $\rho $element comes from the use of new displacement shape functions, which are function of the damage location and, therefore, provide for more efficient models, where accuracy improved at lower computational cost. In numerical tests in the linear regime, coupling between cross sectional rotation and longitudinal vibrations discovered, with longitudinal displacements suddenly changing direction at the damage location and with a peculiar change in the cross section rotation at the same place. Geometrically nonlinear, forced vibrations then investigated in the time domain using Newmark’s method and further couplings between displacement components found.
Forced transverse vibration of a Rayleigh doublebeam system with effect of compressive axial load with a Pasternak layer inbetween are studied in the present paper.
2. Mathematical model
The following assumptions are considered: (a) the behavior of the beam material is linear elastic; (b) the crosssection is rigid and constant throughout the length of the beam and has one plane of symmetry; (c) shear deformations of the crosssection of the beam are considered while elastic axial deformations are ignored; (d) the equations are derived considering the geometric axial deformations; (e) the axial forces $F$ acting on the ends of the beam are not changed with time; (f) the two beams have the same effective material constants.
Fig. 1 shows doublebeam system with Pasternak layer inbetween with length of $l$ subjected to axial compressions ${F}_{1}$ and ${F}_{2}$ which are positive in compression, and arbitrarily distributed transverse continuous loads ${f}_{1}$ and ${f}_{2}$ that are positive when they act downward.
Fig. 1Doublebeam complex system with a Pasternak middle layer
An element of deflected differential layeredbeam of length $dx$ with Pasternak layer between two crosssections taken normal to the deflected axis to the beam is shown in Fig. 2.
Fig. 2Deflected differential layeredbeam element with a Pasternak middle layer
The equations of forced transverse vibration of a Rayleigh doublebeam system with the effect of compressive axial load with a Pasternak layer inbetween have the following form:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}+K\left({w}_{1}{w}_{2}\right){G}_{0}\left(\frac{{\partial}^{2}{w}_{1}}{\partial {x}^{2}}\frac{{\partial}^{2}{w}_{2}}{\partial {x}^{2}}\right)={f}_{1}\left(x,t\right),$
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}+K\left({w}_{2}{w}_{1}\right){G}_{0}\left(\frac{{\partial}^{2}{w}_{2}}{\partial {x}^{2}}\frac{{\partial}^{2}{w}_{1}}{\partial {x}^{2}}\right)={f}_{2}\left(x,t\right).$
The Eqs. (1) and (2) can be reduced to fourthorder partial differential equations for forced vibration of the Rayleigh doublebeam model:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}{C}_{r1}^{2}\frac{{\partial}^{4}{w}_{1}}{\partial {x}^{2}\partial {t}^{2}}+{H}_{1}\left({w}_{1}{w}_{2}\right)={m}_{1}{f}_{1}\left(x,t\right),$
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}{C}_{r2}^{2}\frac{{\partial}^{4}{w}_{2}}{\partial {x}^{2}\partial {t}^{2}}+{H}_{2}\left({w}_{2}{w}_{1}\right)={m}_{2}{f}_{2}\left(x,t\right),$
where:
The initial conditions in general form and boundary conditions for simply supported beams of the same length $l$ are assumed as follows:
${w}_{i}^{\mathrm{\text{'}}\mathrm{\text{'}}}\left(0,t\right)={w}_{i}\left(0,t\right)={w}_{i}\left(l,t\right)={w}_{i}^{\mathrm{\text{'}}\mathrm{\text{'}}}\left(l,t\right)=0,\mathrm{}\mathrm{}\mathrm{}i=1,\mathrm{}2.$
3. Solution of equations
The Eqs. (3) and (4) representing forced vibrations of a Rayleigh doublebeam system with Pasternak layer inbetween. The natural frequencies and the corresponding mode shapes of the system should be obtained by solving the undamped free vibration with appropriate boundary conditions. Assuming time harmonic motion and using separation of variables, the solutions to Eqs. (3) and (4) with the governing boundary conditions Eq. (5) can be written in the form:
where,
Substitution of Eq. (6) into Eqs. (3), (4) yields ordinary differential equations for the Rayleigh doublebeam system. Therefore:
where:
${N}_{1}={C}_{b1}^{2}{k}_{n}^{4}+{H}_{1},\mathrm{}\mathrm{}\mathrm{}\mathrm{}{N}_{2}={C}_{b2}^{2}{k}_{n}^{4}+{H}_{2}$
4. Forced vibration of a Rayleigh doublebeam system
First we should solve the undamped free vibration of the doublebeam system. Therefore:
The solutions of Eqs. (9) and (10) can be assumed to have the following forms:
Substituting Eq. (11) into Eqs. (9), (10) results in the following system of homogeneous algebraic equations for the unknown constants:
When the determinant of the coefficients in Eqs. (12), (13) vanishes, nontrivial solutions for the constants ${C}_{n}$ and ${D}_{n}$ can be obtained, which yields the following frequency (characteristic) equation:
$+({N}_{1}{F}_{1}{\eta}_{1})({N}_{2}{F}_{2}{\eta}_{2}){H}_{1}{H}_{2}{G}_{0}({H}_{1}{\eta}_{2}+{H}_{2}{\eta}_{1}{N}_{1}{\eta}_{2}$
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}{N}_{2}{\eta}_{1}+{F}_{1}{\eta}_{1}{\eta}_{2}+{F}_{2}{\eta}_{1}{\eta}_{2})=0.$
The roots of Eq. (14) are:
where:
$+2{J}_{2}^{2}{G}_{0}{\eta}_{1}({N}_{1}{F}_{1}{\eta}_{1})+2{J}_{1}^{2}{G}_{0}{\eta}_{2}({N}_{2}{F}_{2}{\eta}_{2})$
$+2{J}_{1}{J}_{2}{G}_{0}\left(2{H}_{1}{\eta}_{2}+2{H}_{2}{\eta}_{1}{N}_{1}{\eta}_{2}{N}_{2}{\eta}_{1}+{F}_{1}{\eta}_{1}{\eta}_{2}+{F}_{2}{\eta}_{1}{\eta}_{2}\right).$
For each of the natural frequencies, the associated amplitude ratio of vibration modes of the two beams is given by:
Applying the classical modal expansion method, the forced vibrations of beams subjected to arbitrarily distributed continuous loads can be determined. Particular solutions of nonhomogeneous differential Eqs. (3) and (4) representing forced vibrations of a Rayleigh double model can be assumed in the following:
Introduction of Eqs. (18), (19) into Eqs. (3), (4) results in:
By multiplying the Eqs. (20), (21) by the eigenfunction ${X}_{m}$, then integrating them with respect to $x$ from $0$ to $l$ and using orthogonality condition, we have:
where:
Now we have:
From Eqs. (17), (23), (24) after some algebraic manipulation we obtain:
From Eqs. (25), (26), we obtain:
where:
Also, ${M}_{1}={m}_{1}/{J}_{1}l$ and ${M}_{2}={m}_{2}/{J}_{2}l$. By combining Eqs. (18), (19) and (27), the forced vibrations of an elastically connected Rayleigh doublebeam system with Pasternak layer inbetween can be described by:
Now these general solutions Eqs. (30) and (31) are used to find the vibrations of the two coupled Rayleigh beams.
In the following, we conduct an analysis of forced vibrations for the case of uniformly distributed continuous harmonic load. For simplicity of further analysis, it is assumed that only one of the two beams is subjected to the exciting load. Without loss of generality, we suppose (see Fig. 3):
Fig. 3Doublebeam complex system with the Pasternak layer inbetween subjected to harmonic uniform distributed continuous load
Substituting Eq. (32) into Eqs. (28), (29), we obtain:
where:
Substituting Eqs. (33), (34) into Eqs. (30), (31) gives:
where:
${B}_{n\mathrm{I}}={M}_{Rn}\frac{{\alpha}_{n\mathrm{I}\mathrm{I}}}{{\alpha}_{n\mathrm{I}\mathrm{I}}{\alpha}_{n\mathrm{I}}}\left[\frac{\mathrm{\Omega}}{{\omega}_{n\mathrm{I}}({\mathrm{\Omega}}^{2}{\omega}_{n\mathrm{I}}^{2})}\right],\mathrm{}\mathrm{}\mathrm{}{B}_{n\mathrm{I}\mathrm{I}}={M}_{Rn}\frac{{\alpha}_{n\mathrm{I}}}{{\alpha}_{n\mathrm{I}\mathrm{I}}{\alpha}_{n\mathrm{I}}}\left[\frac{\mathrm{\Omega}}{{\omega}_{n\mathrm{I}\mathrm{I}}({\omega}_{n\mathrm{I}\mathrm{I}}^{2}{\mathrm{\Omega}}^{2})}\right].$
Neglecting the free response, the forced vibrations of the Rayleigh doublebeam system can be obtained by:
The following fundamental conditions of resonance and dynamic vibration absorption have practical significance:
1 – Resonance:
2 – Dynamic vibration absorption:
5. Numerical analysis
For simplicity, it is assumed that both beams are geometrically and physically identical:
Table 1The values of the parameters characterizing properties of the system
$l$ (m)  $E$ (N/m^{2})  $A$ (m^{2})  $K$ (N/m^{2})  $\rho $ (kg/m^{3})  $I$ (m^{4}) 
10  1×10^{10}  5×10^{2}  2×10^{5}  2×10^{3}  4×10^{4} 
If the axial compressions vanish, for the Rayleigh doublebeam we have:
where:
$+2{J}_{1}^{2}{G}_{0}{\eta}_{2}{N}_{2}+2{J}_{1}{J}_{2}{G}_{0}\left(2{H}_{1}{\eta}_{2}+2{H}_{2}{\eta}_{1}{N}_{1}{\eta}_{2}{N}_{2}{\eta}_{1}\right).$
Also:
To determine the effect of compressive axial load and shear foundation modulus of Pasternak layer on the steadystate vibration amplitudes ${A}_{n1}$ and ${A}_{n2}$ of the system, the results under compressive axial load and those without axial load are compared. Introducing the relation:
Using nondimensional ratio:
where, $P={P}_{ER}=EI{\pi}^{2}/{l}^{2}$ is known as the Euler load for Euler and Rayleigh beam, which is the smallest load at which the single beam ceases to be in stable equilibrium under axial compression. For the case of uniformly distributed harmonic load, the steadystate vibration amplitudes for the Rayleigh beams $({A}_{n1}$, ${A}_{n2}$, ${A}_{n1}^{0}$, ${A}_{n2}^{0})$ can be determined. We accounted the results for three cases of shear foundation modulus of Pasternak layer ${G}_{0}=\mathrm{}$0, 100, 1000. If we introduce ${G}_{0}=0$ in the whole equations of this paper, we can obtain the vibration equations of the doublebeam system on a Winkler elastic layer that the results are verified by comparing with those available in [12]. With the vibration mode number $n=\mathrm{}$3 and the exciting frequency $\mathrm{\Omega}=\text{0.6}{\omega}_{n\mathrm{I}\mathrm{I}}$, the effects of Pasternak layer on the steadystate vibration amplitudes ${A}_{n1}$ and ${A}_{n2}$ of the Rayleigh beam represented by the ratios ${\phi}_{1}$ and ${\phi}_{2}$, shown in Figs. 47, respectively.
Fig. 4Relationship between ratio φ1 and parameter s for different axial compression ratio ζ and G0= 100
Fig. 5Relationship between ratio φ2 and parameter s for different axial compression ratio ζ and G0= 100
As can be seen, the ratio ${\phi}_{1}$ decreases with the increase of the axial compression, which implies that the magnitude of the steadystate vibration amplitude ${A}_{n1}$ becomes smaller when the axial compression increases and the ratio ${\phi}_{2}$ increases with the increase of the axial compression, which implies that the magnitude of the steadystate vibration amplitude ${A}_{n2}$ becomes larger when the axial compression increases. It can be observed that the effect of compressive axial load on the magnitude of ${A}_{n1}$ is almost independent of the axial compression ratio $\zeta $ of the two beams whereas it is significantly dependent on the magnitude of ${A}_{n2}$.
Fig. 6Relationship between ratio φ1 and dimensionless parameter s for different shear foundation modulus of Pasternak layer and ζ= 0.5
Fig. 7Relationship between ratio φ2 and dimensionless parameter s for different shear foundation modulus of Pasternak layer and ζ= 0.5
Table 2Effects of shear foundation modulus of Pasternak layer on the steadystate vibration amplitudes ratios φ1 and φ2 for different compressive axial ratio ζ
$s$  
1  0.8  0.6  0.4  0.2  0  
${\phi}_{1}(\zeta =\text{0.1})$  
1.1999  1.1544  1.112  1.0723  1.0351  1  $({G}_{0}=\text{0})$ 
1.1998  1.1543  1.1119  1.0723  1.035  0.9999  $({G}_{0}=\text{100})$ 
1.1987  1.1533  1.111  1.0715  1.0343  0.9993  $({G}_{0}=\text{1000})$ 
${\phi}_{1}(\zeta =\text{0.5})$  
1.1965  1.1526  1.1112  1.0722  1.0351  1  $({G}_{0}=0)$ 
1.1963  1.1525  1.1111  1.0721  1.0351  0.9999  $({G}_{0}=\text{100})$ 
1.1952  1.1515  1.1102  1.0713  1.0343  0.9993  $({G}_{0}=\text{1000})$ 
${\phi}_{1}(\zeta =\text{0.9})$  
1.194  1.1513  1.1107  1.072  1.0352  1  $({G}_{0}=\text{0})$ 
1.1938  1.1512  1.1106  1.072  1.0351  0.9999  $({G}_{0}=\text{100})$ 
1.1927  1.1501  1.1097  1.0711  1.0344  0.9993  $({G}_{0}=\text{1000})$ 
${\phi}_{2}(\zeta =\text{0.1})$  
1.22  1.1699  1.1232  1.0796  1.0386  1  $({G}_{0}=\text{0})$ 
1.2108  1.1611  1.1147  1.0714  1.0308  0.9925  $({G}_{0}=\text{100})$ 
1.1287  1.0824  1.0393  0.999  0.9611  0.9255  $({G}_{0}=\text{1000})$ 
${\phi}_{2}(\zeta =\text{0.5})$  
1.302  1.2332  1.169  1.109  1.0528  1  $({G}_{0}=\text{0})$ 
1.2921  1.2238  1.1602  1.1006  1.0448  0.9925  $({G}_{0}=\text{100})$ 
1.2044  1.1408  1.0816  1.0261  0.9742  0.9255  $({G}_{0}=\text{1000})$ 
${\phi}_{2}(\zeta =\text{0.9})$  
1.3944  1.3028  1.2181  1.1398  1.0673  1  $({G}_{0}=\text{0})$ 
1.3839  1.2929  1.2089  1.1312  1.0592  0.9925  $({G}_{0}=\text{100})$ 
1.2897  1.2051  1.1269  1.0546  0.9876  0.9255  $({G}_{0}=\text{1000})$ 
As we said before in this paper the effects of Pasternak layer on forced transverse vibration of a Rayleigh doublebeam system with effect of compressive axial load is analyzed. Effects of Pasternak layer cannot be observed in Figs. 47. Numerical values of the ratios ${\phi}_{1}$ and ${\phi}_{2}$ for different shear foundation modulus of Pasternak layer with different axial compression ratio $\zeta $ of the two beams is shown in Table 2. Numerical results of the ratios ${\phi}_{1}$ and ${\phi}_{2}$ show that the difference with ratios ${\phi}_{1}$ and ${\phi}_{2}$ decrease with the increase of shear modulus of the Pasternak layer for different axial compression ratio $\zeta $, but It can be observed that the effect of shear foundation modulus of Pasternak layer on the magnitude of ${A}_{n1}$ is almost independent of the axial compression ratio $\zeta $ of the two beams whereas it is significantly dependent on the magnitude of ${A}_{n2}$.
Fig. 8Relationship between ratio φ1 and dimensionless parameter s for different mode number n and G0= 100
Fig. 9Relationship between ratio φ2 and dimensionless parameter s for different mode number n and G0= 100
Fig. 10Relationship between ratio φ1 and dimensionless parameter s for different shear modulus of Pasternak (n= 5)
Fig. 11Relationship between ratio φ2 and dimensionless parameter s for different shear modulus of Pasternak (n=5)
On the other hand, with the axial compression ratio $\zeta =\text{0.5}$ and the exciting frequency $\mathrm{\Omega}=\text{0.6}{\omega}_{n\mathrm{I}\mathrm{I}}$, the effects of Pasternak layer on the steadystate vibration amplitudes ${A}_{n1}$ and ${A}_{n2}$ are shown in Figs. 811 for different mode shape number $n$, respectively. It can be seen that, with the same axial compression, the ratios ${\phi}_{1}$ and ${\phi}_{2}$ diminish with the increasing vibration mode number $n$, which implies that the magnitudes of the steadystate vibration amplitudes ${A}_{n1}$ and ${A}_{n2}$ get smaller when the vibration mode number $n$ becomes larger.
Table 3 shows the effects of compressive axial load and shear modulus of Pasternak layer on difference between the steadystate vibration amplitudes ratios ${\phi}_{1}$ and ${\phi}_{2}$ of the Rayleigh beam for different mode number $n$. It can be observed that the differences between ratios ${\phi}_{1}$ and ${\phi}_{2}$ of the Rayleigh beam increase with increasing the dimensionless parameter $s$. Also we can see that the ratios ${\phi}_{1}$ and ${\phi}_{2}$ decrease with increasing of the shear foundation modulus of Pasternak layer for different mode shape numbers.
Table 3Effects of shear foundation modulus of Pasternak layer on the steadystate vibration amplitudes ratios φ1 and φ2 for different mode shape number n
$s$  
1  0.8  0.6  0.4  0.2  0  
${\phi}_{1}(n=\text{3})$  
1.1965  1.1526  1.1112  1.0722  1.0351  1  $({G}_{0}=\text{0})$ 
1.1963  1.1525  1.1111  1.0721  1.0351  0.9999  $({G}_{0}=\text{100})$ 
1.1952  1.1515  1.1102  1.0713  1.0343  0.9993  $({G}_{0}=\text{1000})$ 
${\phi}_{1}(n=\text{5})$  
1.0654  1.0518  1.0385  1.0254  1.0126  1  $({G}_{0}=\text{0})$ 
1.0654  1.0518  1.0384  1.0254  1.0125  1  $({G}_{0}=\text{100})$ 
1.065  1.0514  1.038  1.0250  1.0123  0.9998  $({G}_{0}=\text{1000})$ 
${\phi}_{1}(n=\text{7})$  
1.0327  1.026  1.0194  1.0129  1.0064  1  $({G}_{0}=\text{0})$ 
1.0327  1.026  1.0194  1.0128  1.0064  1  $({G}_{0}=\text{100})$ 
1.0324  1.0258  1.0192  1.0127  1.0063  0.9999  $({G}_{0}=\text{1000})$ 
${\phi}_{2}(n=\text{3})$  
1.302  1.2332  1.169  1.109  1.0528  1  $({G}_{0}=\text{0})$ 
1.2921  1.2238  1.1602  1.1006  1.0448  0.9925  $({G}_{0}=\text{100})$ 
1.2044  1.1408  1.0816  1.0261  0.9742  0.9255  $({G}_{0}=1000)$ 
${\phi}_{2}(n=\text{5})$  
1.0992  1.0784  1.0581  1.0382  1.0189  1  $({G}_{0}=\text{0})$ 
0.9331  0.9154  0.8982  0.8814  0.8649  0.8489  $({G}_{0}=\text{100})$ 
–0.5249  –0.515  –0.5053  –0.4959  –0.4867  –0.4777  $({G}_{0}=\text{1000})$ 
${\phi}_{2}(n=\text{7})$  
1.0493  1.0392  1.0292  1.0193  1.0096  1  $({G}_{0}=\text{0})$ 
–0.1343  –0.133  –0.1318  –0.1305  –0.1292  –0.1280  $({G}_{0}=\text{100})$ 
–10.2863  –10.1877  –10.0903  –9.9940  –9.8988  –9.8047  $({G}_{0}=\text{1000})$ 
6. Conclusions
In vibration analysis, a dynamic vibration absorber, such as Rayleigh double beam system, reduces or eliminates the vibration of a harmonically excited system. Rotating machines such as engines, motors, and pumps often incite vibration due to rotational imbalances. A dynamic absorber can be affixed to the rotating machine and tuned to oscillate in such a way that exactly counteracts the force from the rotating imbalance. This reduces the possibility that a resonance condition will occur, which can cause rapid catastrophic failure. Properly implemented, a dynamic absorber will neutralize the undesirable vibration, which would otherwise reduce service life or cause mechanical damage. Based on the Rayleigh beam theory, the forced transverse vibration of an elastically connected simply supported Rayleigh doublebeam, with a Pasternak interlayer under compressive axial loading for one case of particular excitation loading are studied. The dynamic response of the system caused by arbitrarily distributed continuous loads is obtained. The magnitudes of the steadystate vibration amplitudes of the beam are dependent on the axial compression and shear foundation modulus of Pasternak layer. Also we concluded that the ratios ${\phi}_{1}$ and ${\phi}_{2}$ decrease with increasing of the shear foundation modulus of Pasternak layer $\left({G}_{0}\right)$. Founded analytical forms can be used in the optimal design of a new type of a dynamic vibration absorber. Thus the Rayleigh beamtype dynamic absorber with Pasternak layer inbetween can be applied to suppress the excessive vibrations of corresponding beam systems and be more effective than the doublebeam system with a Winkler elastic layer. The beamtype dynamic damper is an accepted concept for a Continuous Dynamic Vibration Absorber (CDVA).
References

Seelig J. M., Hoppmann W. H. Impact on an elastically connected doublebeam system. Transactions of the American Society of Mechanical Engineers, Journal of Applied Mechanics, Vol. 31, 1964, p. 621626.

Oniszczuk Z. Free transverse vibrations of elastically connected simply supported doublebeam complex system. Journal of Sound and Vibration, Vol. 232, Issue 2, 2000, p. 387403.

Oniszczuk Z. Forced transverse vibrations of an elastically connected complex simply supported doublebeam system. Journal of Sound and Vibration, Vol. 264, Issue 2, 2003, p. 273286.

Zhang Y. Q., Lu Y., Wang S. L., Liu X. Vibration and buckling of a doublebeam system under compressive axial loading. Journal of Sound and Vibration, Vol. 318, Issue 1, 2008, p. 341352.

Zhang Y. Q., Lu Y., Ma G. W. Effect of compressive axial load on forced transverse vibrations of a doublebeam system. International Journal of Mechanical Sciences, Vol. 50, Issue 2, 2008, p. 299305.

Pasternak P. L. On a New Method of Analysis of an Elastic Foundation by Means of Two Foundation Constants. Moscow, 1954.

Wang T. W., Stephens J. E. Natural frequencies of Timoshenko beams on Pasternak foundations. Journal of Sound and Vibration, Vol. 51, Issue 2, 1977, p. 149155.

De Rosa M. A. Free vibrations of Timoshenko beams on twoparameter elastic foundation. Computers and Structures, Vol. 57, Issue 1, 1995, p. 151156.

Stojanovic V., Kozic P., Pavlovic R., Janevski G. Effect of rotary inertia shear on vibration and buckling of a double beam system under compressive axial loading. Archive of Applied Mechanics, Vol. 81, Issue 12, 2011, p. 19932005.

Stojanovic V., Kozic P., Janevski G. Buckling instabilities of elastically connected Timoshenko beams on an elastic layer subjected to axial forces. Journal of Mechanics of Materials and Structures, Vol. 7, Issue 4, 2012, p. 363374.

Stojanovic V., Kozic P. Forced transverse vibration of Rayleigh and Timoshenko doublebeam system with effect of compressive axial load. International Journal of Mechanical Sciences, Vol. 60, Issue 1, 2012, p. 5971.

Stojanovic V., Kozic P., Janevski G. Exact closedform solutions for the natural frequencies and stability of elastically connected multiple beam system using Timoshenko and highorder shear deformation theory. Journal of Sound and Vibration, Vol. 332, Issue 3, 2013, p. 563576.

Li J., Hua H. Spectral finite element analysis of elastically connected doublebeam systems. Finite Elements in Analysis and Design, Vol. 43, Issue 15, 2007, p. 11551168.

Li J., Chen Y., Hua H. Exact dynamic stiffness matrix of a Timoshenko threebeam system. International Journal of Mechanical Sciences, Vol. 50, Issue 6, 2008, p. 10231034.

Ariaei A., ZiaeiRad S., Ghayour M. Transverse vibration of a multipleTimoshenko beam system with intermediate elastic connections due to a moving load. Archive of Applied Mechanics, Vol. 81, Issue 3, 2011, p. 263281.

Mao Q. Free vibration analysis of elastically connected multiplebeams by using the Adomian modified decomposition method. Journal of Sound and Vibration, Vol. 331, Issue 11, 2012, p. 25322542.

Kozic P., Janevski G., Pavlovic R. Moment Lyapunov exponents and stochastic stability of a doublebeam system under compressive axial loading. International Journal of Solids and Structures, Vol. 47, Issue 10, 2010, p. 14351442.

Stojanovic V., Petkovic M. Moment Lyapunov exponents and stochastic stability of a three dimensional system on elastic foundation using a perturbation approach. Journal of Applied Mechanics, Vol. 80, Issue 5, 2013, p. 051009.

Stojanovic V., Ribeiro P., Stoykov S. Nonlinear vibration of Timoshenko damaged beams by a new ρversion finite element method. Computers and Structurers, Vol. 120, 2013, p. 107119.