Published: December 22, 2025

Vibration mechanics of a semi-cylindrical cushion with inertial properties

Reyhan Akbarli1
Parviz Garayev2
Ruzigul Sayfutdinova3
Khadicha Urolova4
1, 2Department of Mechanics, Azerbaijan University of Architecture and Construction, Baku, Azerbaijan
3, 4Department of Civil Engineering, Tashkent State Transport University, Tashkent, Uzbekistan
Corresponding Author:
Reyhan Akbarli
Views 18
Reads 6
Downloads 37

Abstract

This study investigates the free and forced vibrations of an inertial semi-cylindrical cushion using the Lamé displacement–based formulation and the Hamilton-Ostrogradsky variational principle. Natural frequencies, vibration modes, and normal displacements were obtained for different boundary conditions. Special attention was given to the influence of inertial and elastic parameters, including density, Young’s modulus, Poisson’s ratio, and foundation inertia. Analytical relations describing the variation of vibration amplitudes and frequencies with respect to geometric and material characteristics were derived. A comparative analysis between the analytical solution and a 3D finite element model developed in ABAQUS demonstrated strong agreement, confirming the accuracy of the proposed formulation. The results provide practical guidelines for optimizing vibration-isolation elements and improve the design of structures exposed to dynamic loads. The study contributes to Sustainable Development Goal 9 by supporting the development of safer and more durable infrastructure components.

1. Introduction

Semi-cylindrical cushions are widely applied in engineering systems exposed to dynamic loading, where accurate prediction of vibration behavior is essential for ensuring structural safety and performance. Previous studies have focused primarily on circular and multilayer plates resting on elastic foundations. Leonenko [1-3] investigated free vibrations of three-layer circular plates on various foundations, providing analytical formulations for frequency prediction. Classical elasticity theory by Novatsky [4] and structural models developed by Amiro and Zarutsky [5] form the theoretical basis for analyzing complex curved bodies. Additional studies have examined shells interacting with elastic and liquid media [6, 7], highlighting the importance of including inertial effects in vibration analysis.

Recent research also emphasizes the role of material innovations and numerical modeling in dynamic systems. Studies on pump components [8], machine elements [9], water-related infrastructure [10, 11], and eco-concretes [12] demonstrate the growing need for reliable vibration models in engineering applications. Composite and fiber-based materials have also been analyzed for acoustic and vibrational properties [13], while surface-treatment studies [14] and deformation modeling in transportation engineering [15] further confirm the relevance of accurate structural vibration analysis across multiple domains.

However, despite these advances, the vibration behavior of inertial semi-cylindrical cushions has not been sufficiently studied. Unlike flat or circular plates, semi-cylindrical elements involve curvature-induced stiffness variation, complex boundary conditions, and significant sensitivity to foundation inertia. These factors necessitate a dedicated analytical approach capable of capturing both geometric and inertial characteristics.

Therefore, this study develops an analytical model for the free vibration of a semi-cylindrical cushion based on the Hamilton-Ostrogradsky variational principle and the Lamé displacement formulation. Material and geometric parameters are explicitly defined, and a numerical algorithm is constructed to compute natural frequencies and displacement fields. The analytical solution is validated through a 3D finite element model developed in ABAQUS to ensure accuracy and reliability.

The cushion material properties are: density ρs= 7850 kg/m3, Young’s modulus Es= 2.1×1011 Pa, and Poisson’s ratio νs= 0.30. The geometric parameters are radius R= 0.25 m, thickness hs= 0.012 m, and length L= 0.40-0.80 m. Foundation parameters include stiffness kφ= 1.2×106 N/m3 and inertia md= 150 kg·s2/m4.

2. Problem statement

This study investigates the vibrations of an inertial semi-cylindrical cushion positioned between a bridge span and its support pier (Fig. 1). The aim is to determine the cushion’s dynamic characteristics, assess the effects of its inertial and elastic properties, and establish how vibration amplitudes and frequencies depend on structural parameters.

Fig. 1Inertial semi-cylindrical cushion: p – load applied to the cushion from the bridge side, qr – load applied to the cushion from the support side

Inertial semi-cylindrical cushion: p – load applied to the cushion from the bridge side,  qr – load applied to the cushion from the support side

To solve the formulated problem, the variational Hamilton-Ostrogradsky principle [4, 5] is applied:

1
δW=0,

where, W=t'tJdt is the action functional, and, t'and t'' are arbitrarily chosen moments of time:

2
J=Vks-Vps-Ap-Aq,

where, Vps, Vks are the potential and kinetic energies of the semi-cylindrical cushion, respectively; Aq is the work of the force acting on the cushion from the bridge side during its vertical displacement, and Ap is the work of the force acting from the support side. These quantities are calculated as follows:

3
Vps=λs+2μs2e112+e222+e332+λs+2μse11e22+e22e33+e11e33,
4
Vks=ρs2Vsxt2+sθt2+srt2dxdydr,
5
Aq=-Sqrsrr=0dxdy,
6
Ap=-pxsrr=R,θ=π2dx.

The quantities sx, sθ, sr denote the displacements of the cushion points along the corresponding coordinate directions, and t represents time. The strain components e11, e22, e33 are expressed in terms of these displacements sx, sθ, sr as follows:

7
e11=srr,     e22=Rsθθ+sr,     e33=sxx.

It is assumed that the force qr, acting from the support side on the semi-cylindrical cushion during its vertical displacement sr, follows the Winkler-type law:

8
qr=kϑsr+md2srt2,

where kϑ is the stiffness coefficient, and md is the specific weight of the support material. The load from the bridge side acting on the cushion is defined by the function:

9
Px=p0 cosnπ2sinkxsinωt.

The displacements of the inertial semi-cylindrical cushion sx, sθ, sr are described by the Lamé equations in vector form [6, 7]:

10
al2graddiv u-at2rotrotu+ρs2ut2=0,

where at=λs+2μsρs, ae=μsρs are the longitudinal and transverse wave propagation velocities, ρs is the density, λs and μs are the Lamé coefficients.

3. Numerical method

The analytical expressions obtained from the assumed displacement functions were evaluated numerically using MATLAB. All integrals appearing in the kinetic and strain energy formulations were computed using standard Gaussian quadrature. Substituting these expressions into the Hamilton-Ostrogradsky functional leads to a system of algebraic equations that can be written in matrix form as AωC=B, where Aω depends on the vibration frequency and the material and geometric parameters of the cushion.

The natural frequencies were determined by solving the characteristic equation detAω= 0.

The nonlinear equation was solved using the Newton–Raphson method, ensuring rapid convergence. After obtaining each resonant frequency, displacement coefficients 𝐶 were calculated using Cramer’s rule and substituted into the assumed displacement field to determine mode shapes. Stability and convergence were verified by refining the integration scheme and testing multiple mesh densities, confirming the reliability of the algorithm.

4. Solution

The displacements of the inertial semi-cylindrical cushion sx, sθ, sr are given in the following form [6]:

11
sx=A~skInγer-C~sγt2μtInγtrcosnθcoskxsinωt,
sθ=-A~snrInγer-C~snkrμtInγtr-B~snInγtrrsinnθsinkxsinωt,
sr=A~sInγerr-C~skμtInγtrr+B~snrInγtrcosnθsinkxsinωt,

where A~s, C~s, B~s are unknown constant coefficients, and k, n, γe,γt – are the wave numbers and γe2=k2-μe2, γt2=k2-μt2.

By substituting Eq. (10) into Eqs. (3-6), the expressions for the potential and kinetic energies of the cushion are obtained:

12
Vps=πlλs+2μs16A~s20R[γl4In''2γlr+k4In2γlr-2k2γl2InγlrIn''γlr
      +β1R2+2γl2RIn''γlr-2k2RInγlr]dr+B~s20Rβ22+R2β3+β2β3dr
      +C~s20Rk2γt4μt2In''2γtr+k2γt4μt2In2γtr-k2γt4μt2InγtrIn''γtr+β4kγt2μt
      ×RInγtr-In''γtr]dr+A~sC~s0R-2kγl2γt2μtIn''γlrIn''γtr-2γt2k3μtInγlr
      ×Inγtr+kγl2γt2μtInγtrIn''γlr+γt2k3μtInγlrIn''γtr+2β4β5
      +β4γl2RIn''γlr-k2RInγlr+β5kγt2μtInγlr-RIn''γtr]dr
      +A~sB~s0R[2γl2β2In''γtr+2R2β1β3+γl2β3In''γlr+β1β2-k2β2Inγlr
      -k2Rβ3Inγlr]dr+B~sC~s0R-2kγt2β2μtIn''γtr+2R2β3β4-2kγt2Rβ3μt
      ×In''γtr+kγt2β2μtInγtr+kγt2Rβ3μtInγtr+Rβ2β4drsin2ωt,
β1r=-n2rInγlr+γlIn'γlr,       β2r=nγtrIn'γtr-nr2Inγtr,
β3r=-γtIn'γtr+nrInγtr,        β4r=-kμtn2rInγtr-γtIn'γtr,
β5r=-n2rInγlr+γlIn'γlr,
Vks=πlω2ρs8A~s20Rk2In2γlr+n2r2 In2γlr+γl2In'2γlrdr
      +B~s20Rγt2n2In'2γlr+n2r2 In2γtrdr+C~s20Rγt4μt2Inγtr+n2k2r2μt2In2γtr
      +k2γt2μt2In'2γtrdr+A~sC~s0R-2kγt2InγlrInγtr+2n2kr2μtInγlrInγtr
      -2kγlγtμtIn'γlrIn'γtrdr+A~sB~s0R2γtrInγlrIn'γtr
      +2nγlrInγtrIn'γlrdr+B~sC~s0R2kγtμtInγtrIn'γtr
      -2nkγtrμtInγtrIn'γtrdrsin2ωt,
Aq=-πlR4[As20Rγl2In'2γlrdr+C~s2k2γt2μt20RIn'2γtrdr+B~s2n20RIn2γtrr2dr
      -A~sC~s2kγlγtμt0RIn'γlrIn'γtrdr+A~sB~s2γln0RInγtrIn'γlrrdr
      -B~sC~s2nkγtμt0RInγtrIn'γtrrdrkϑ-ω2mdsin2ωt,
Ap=-p0l2A~sInγeRr-C~skμtInγtRr+B~snRInγtRsin2ωt.

By substituting all expressions from Eq. (11) into Eq. (2) for the J and taking t'=0, t''=π/ω, and applying the Hamilton-Ostrogradsky variational principle Eq. (1), the following expression is obtained:

13
W=πlω2ρs8A~s20Rk2In2γlr+n2r2 In2γlr+γl2In'2γlrdr
      +B~s20Rγt2n2In'2γlr+n2r2 In2γtrdr+C~s20Rγt4μt2Inγtr+n2k2r2μt2In2γtr
      +k2γt2μt2In'2γtr]dr+A~sC~s0R-2kγt2InγlrInγtr+2n2kr2μtInγlrInγtr
      -2kγlγtμtIn'γlrIn'γtr]dr+A~sB~s0R2γtrInγlrIn'γtr
      +2nγlrInγtrIn'γlr]dr+B~sC~s0R2kγtμtInγtrIn'γtr
      -2nkγtrμtInγtrIn'γtr]dr-πlλs+2μs16A~s20R[γl4In''2γlr+k4In2γlr
      -2k2γl2InγlrIn''γlr+β1R2+2γl2RIn''γlr-2k2RInγlr]dr
      +B~s20Rβ22+R2β3+β2β3dr+C~s20Rk2γt4μt2In''2γtr+k2γt4μt2In2γtr
      -k2γt4μt2InγtrIn''γtr+β4kγt2μt×RInγtr-In''γtrdr
      +A~sC~s0R-2kγl2γt2μtIn''γlrIn''γtr-2γt2k3μtInγlr×Inγtr
      +kγl2γt2μtInγtrIn''γlr+γt2k3μtInγlrIn''γtr+2β4β5
      +β4γl2RIn''γlr-k2RInγlr+β5kγt2μtInγlr-RIn''γtrdr
      +A~sB~s0R[2γl2β2In''γtr+2R2β1β3+γl2β3In''γlr+β1β2-k2β2Inγlr
      -k2Rβ3Inγlr]dr+B~sC~s0R-2kγt2β2μtIn''γtr+2R2β3β4-2kγt2Rβ3μt
      ×In''γtr+kγt2β2μtInγtr+kγt2Rβ3μtInγtr+Rβ2β4]dr
      +πlR4As20Rγl2In'2γlrdr+C~s2k2γt2μt20RIn'2γtrdr+B~s2n20RIn2γtrr2dr
      -A~sC~s2kγlγtμt0RIn'γlrIn'γtrdr+A~sB~s2γln0RInγtrIn'γlrrdr
      -B~sC~s2nkγtμt0RInγtrIn'γtrrdrkϑ-ω2md
      +p0l2A~sInγeRr-C~skμtInγtRr+B~snRInγtRπ2ω.

From Eq. (12) for W, it follows that the functional represents a quadratic form with respect to the constants A~s, B~s, C~s. By varying the functional with respect to these unknowns A~s, B~s, C~s, a system of non-homogeneous algebraic equations is obtained:

WA~s=φ1,     WB~s=φ2,     WC~s=φ3.

In expanded form, the system of equations is as follows:

14
φ11As+φ12Bs+φ13Cs=-lp02,     φ21As+φ22Bs+φ23Cs=-lp02,
φ31As+φ32Bs+φ33Cs=-lp02nInγtRR,

where the coefficients φij (i, j= 1, 2, 3) are obtained from the variational functional W after variation as the coefficients of the constants A~s, B~s, C~s. Using Cramer’s rule, these constants are determined as follows:

15
A~s=1,     B~s=2,     C~s=3,

where, is the main determinant of the system, and i (i= 1, 2, 3) are the auxiliary determinants of the Eq. (14). By substituting the values of A~s, B~s, C~s from Eq. (15) into Eq. (10) for sr, the final formula for the cushion displacement is obtained:

sr=A~sInγeRr-C~skμtInγtRr+B~snRInγtRcosnπ2sinkxsinωt.

The resonance frequencies of the cushion are determined from the condition  =0. The roots of the equation  =0 are found numerically. For the calculations, the following parameters characterizing the material of the semi-cylindrical cushion are adopted: kϑ= 50 MPa/m, md= 1000 MPa/m.

The dependence of the natural frequencies of the semi-cylindrical cushion on its length is shown in Fig. 2. The normal displacements of surface points of the cushion as functions of the coordinates for different values of the foundation’s specific weight are shown in Fig. 3. In all graphs: curve 1 corresponds to md= 1000 MPa/m, curve 2 to md= 1500 MPa/m, and curve 3 to md= 2000 MPa/m.

A decrease in frequency is observed with increasing length due to reduced stiffness-to-mass ratio.

Fig. 2Fundamental natural frequency of the semi-cylindrical cushion as a function of its length

Fundamental natural frequency of the semi-cylindrical cushion as a function of its length

Fig. 3Vertical deflection profile of the semi-cylindrical cushion for different values of foundation inertia md. Curve 1 – low inertia; Curve 2 – medium; Curve 3 – high inertia

Vertical deflection profile of the semi-cylindrical cushion for different values of foundation inertia md. Curve 1 – low inertia; Curve 2 – medium; Curve 3 – high inertia

5. Model validation

A 3D FEM model of the semi-cylindrical cushion was created in ABAQUS to verify the analytical model, replicating all boundary conditions and material properties. Using 8-node C3D8 elements and mesh convergence, natural frequencies were computed via the Lanczos method. Comparison with analytical results showed differences below 2.8 % for the first three modes, confirming the validity of the displacement functions and variational equations. This demonstrates the model’s reliability for parameter studies, optimization, and vibration isolation design.

6. Discussion

Analytical and FEM results reveal the dynamic behavior of the semi-cylindrical cushion, influenced by geometry, material properties, and foundation inertia. Longer cushions lower natural frequencies, while heavier foundations increase vertical displacements. Higher stiffness raises frequencies, whereas higher density lowers them. Close agreement between analytical and FEM models (< 3 % deviation) confirms the formulation’s validity, emphasizing the need to optimize geometry, material, and foundation parameters for effective bridge vibration isolation.

7. Conclusions

1) Increasing the length of the semi-cylindrical cushion leads to a decrease in its natural vibration frequencies.

2) An increase in the specific weight of the support results in higher natural vibration frequencies.

3) An increase in the specific weight of the foundation causes greater deflection of the cushion.

4) The inertial properties of the semi-cylindrical cushion lead to a reduction in its natural frequencies.

References

  • L. D. Leonenko, “Free vibrations of a three-layer circular plate on an inertial foundation under thermal effects,” Izvestiya of F. Skorina Gomel State University, Natural Sciences, No. 6, pp. 174–179, 2019.
  • D. V. Leonenko, “Free vibrations of three-layer circular plates on an elastic foundation,” Ecological Bulletin of Scientific Centers of Black Sea Economic Cooperation, No. 3, pp. 42–47, 2008.
  • D. V. Leonenko, “Free vibrations of three-layer circular plates on a Pasternak elastic foundation,” Ecological Bulletin of Scientific Centers of Black Sea Economic Cooperation, No. 1, pp. 59–63, 2014.
  • V. Novatsky, Theory of Elasticity. Moscow: Mir Publishers, 1975.
  • I. Y. Amiro and V. A. Zarutsky, Theory of Ribbed Shells: Methods for Shell Calculation. Kyiv: Naukova Dumka, 1980.
  • F. S. Latifov, Vibrations of Shells with Elastic and Liquid Medium. Baku, 1999.
  • Y. A. Amanzade, Theory of Elasticity. Moscow: Higher School Publishing, 1976.
  • O. Glovatskii et al., “Modernization of submersible pump designs for sustainable irrigation: A bibliometric and experimental contribution to Sustainable Development Goals (SDGs),” Indonesian Journal of Science and Technology, Vol. 10, No. 3, pp. 427–438, Jun. 2025, https://doi.org/10.17509/ijost.v10i3.85476
  • A. Mamahonov, I. Khikmatillaev, and D. Kholbaev, “Automatic determination and control of working body loads of technological machines,” in International Conference on Physical Research and Engineering Technology Problems, Vol. 3304, p. 030083, Jan. 2025, https://doi.org/10.1063/5.0269315
  • D. Hammoumi, H. S. Al-Aizari, Z. Alkhawlani, S. Chakiri, and Z. Bejjaji, “water quality assessment using the water quality index, and geographic information systems in Nador Canal, Morocco,” Journal of Environmental and Earth Sciences, Vol. 6, No. 2, pp. 1–16, May 2024, https://doi.org/10.30564/jees.v6i2.6310
  • C. Yerli and U. Sahin, “Reducing gray, blue, green, and total water footprint in wastewater irrigated silage maize in zero tillage,” Journal of Water and Climate Change, Vol. 15, No. 12, pp. 5778–5792, Dec. 2024, https://doi.org/10.2166/wcc.2024.278
  • H. Bentegri et al., “Assessment of compressive strength of eco-concrete reinforced using machine learning tools,” Scientific Reports, Vol. 15, No. 1, Feb. 2025, https://doi.org/10.1038/s41598-025-89530-y
  • V. Sekar, S. Palaniyappan, S. Y. Eh Noum, A. Putra, S. Sivanesan, and D. D. Chin Vui Sheng, “Acoustic absorbers made of wood fiber composites developed by compression molding and additive manufacturing,” Wood Research, Vol. 68, No. 1, pp. 68–82, Feb. 2023, https://doi.org/10.37763/wr.1336-4561/68.1.6882
  • Z. Aringozhina, B. Rakhadilov, D. Yeskermessov, Y. Tabiyeva, and W. Ahmed, “Study of the Influence of Electrolytic-Plasma Treatment on the Tribological Properties of High-Speed Steel,” Advances in Chemical and Materials Engineering, pp. 22–82, Sep. 2023, https://doi.org/10.4018/978-1-6684-6830-2.ch002
  • Y. Y. Rustamov and G. N. Akhundov, “Study of deformations caused by passenger and freight trains in the soil bed during movement,” Scientific Works, 2023.

About this article

Received
October 27, 2025
Accepted
November 18, 2025
Published
December 22, 2025
SUBJECTS
Mechanical vibrations and applications
Keywords
dynamic behavior
inertial properties
natural frequencies
semi-cylindrical cushion
vibrations
Acknowledgements

The authors have not disclosed any funding.

Data Availability

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Conflict of interest

The authors declare that they have no conflict of interest.