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 7850 kg/m3, Young’s modulus 2.1×1011 Pa, and Poisson’s ratio 0.30. The geometric parameters are radius 0.25 m, thickness 0.012 m, and length 0.40-0.80 m. Foundation parameters include stiffness 1.2×106 N/m3 and inertia 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

To solve the formulated problem, the variational Hamilton-Ostrogradsky principle [4, 5] is applied:
where, is the action functional, and, and are arbitrarily chosen moments of time:
where, , are the potential and kinetic energies of the semi-cylindrical cushion, respectively; is the work of the force acting on the cushion from the bridge side during its vertical displacement, and is the work of the force acting from the support side. These quantities are calculated as follows:
The quantities , , denote the displacements of the cushion points along the corresponding coordinate directions, and represents time. The strain components , , are expressed in terms of these displacements , , as follows:
It is assumed that the force , acting from the support side on the semi-cylindrical cushion during its vertical displacement , follows the Winkler-type law:
where is the stiffness coefficient, and is the specific weight of the support material. The load from the bridge side acting on the cushion is defined by the function:
The displacements of the inertial semi-cylindrical cushion , , are described by the Lamé equations in vector form [6, 7]:
where , are the longitudinal and transverse wave propagation velocities, is the density, and 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 , where depends on the vibration frequency and the material and geometric parameters of the cushion.
The natural frequencies were determined by solving the characteristic equation 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 , , are given in the following form [6]:
where , , are unknown constant coefficients, and , , , – are the wave numbers and , .
By substituting Eq. (10) into Eqs. (3-6), the expressions for the potential and kinetic energies of the cushion are obtained:
By substituting all expressions from Eq. (11) into Eq. (2) for the and taking , , and applying the Hamilton-Ostrogradsky variational principle Eq. (1), the following expression is obtained:
From Eq. (12) for , it follows that the functional represents a quadratic form with respect to the constants , , . By varying the functional with respect to these unknowns , , , a system of non-homogeneous algebraic equations is obtained:
In expanded form, the system of equations is as follows:
where the coefficients (, 1, 2, 3) are obtained from the variational functional after variation as the coefficients of the constants , , . Using Cramer’s rule, these constants are determined as follows:
where, is the main determinant of the system, and ( 1, 2, 3) are the auxiliary determinants of the Eq. (14). By substituting the values of , , from Eq. (15) into Eq. (10) for , the final formula for the cushion displacement is obtained:
The resonance frequencies of the cushion are determined from the condition . The roots of the equation are found numerically. For the calculations, the following parameters characterizing the material of the semi-cylindrical cushion are adopted: 50 MPa/m, 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 1000 MPa/m, curve 2 to 1500 MPa/m, and curve 3 to 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

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

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
The authors have not disclosed any funding.
The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.
The authors declare that they have no conflict of interest.