Accurate dynamic analysis of beams supported on soil is essential for the safe and economical design of vibration-resistant and earthquake-resistant structures. Earlier contributions, such as the semi-analytical framework developed by Guenfoud, Bosakov, and Laefer (2009), provided reliable predictions for beams on elastic half-spaces with inertial properties. However, these formulations are mathematically intensive and often difficult to apply directly in engineering practice. Conventional soil models, including Winkler’s subgrade reaction method and continuum-based elastic models, either oversimplify soil behaviour or involve complex derivations, further limiting their practical utility. The present study introduces a modified version of the semi-analytical method for beams resting on an elastic half-space with inertial properties, reformulated to reduce mathematical complexity without compromising accuracy. The approach incorporates soil inertia effects while enabling more straightforward evaluation of displacements under vertical loading. To illustrate its effectiveness, the dynamic response of a beam measuring 10 m in length, 1 m in width, and 0.5 m in depth is investigated. The corresponding Eigenfrequencies and natural shapes of the beam are determined, and the results show close agreement with those obtained from established analytical methods. This refinement not only simplifies computational effort but also makes the analysis more accessible for engineering applications. The method thus provides a practical and reliable tool for studying beam–soil interaction dynamics, with direct relevance for the design of safer and more cost-effective structural systems.
Introduction
A fundamental solution (or Green’s function) describes how a solid reacts to a static or dynamic point load at another location.
These solutions are crucial for understanding soil–structure interaction (SSI), especially in the context of elastic and dynamic behavior of soil.
?? 2. Approaches to Modeling Beams on Elastic Soil
Two Main Methods:
Subgrade Reaction Model (Winkler model):
Soil is modeled as independent springs.
Deflection is local — no interaction between adjacent points.
Simplified but lacks realism.
Elastic Continuum Model:
Uses Boussinesq’s solution for a half-space.
Considers soil interaction across regions.
More accurate, but computationally complex.
???? 3. Categories of Dynamic Analysis Methods
Simplified Procedures
Use charts and empirical approximations.
Easy for practicing engineers.
Semi-analytical Methods
Discretize foundation surface.
Use Fast Fourier Transforms (FFT) and impedance functions.
Analytical Methods
Based on wave equations with boundary assumptions.
Dynamic Finite Element Models (FEM)
Include absorbing boundaries.
Suitable for layered or inhomogeneous soils.
Boundary Element Methods (BEM)
Reduce the problem to algebraic equations.
Handle wave propagation and far-field effects.
Experimental Methods
Field/lab tests.
Validate theoretical/numerical models.
???? 4. Method Used in the Study
A semi-analytical method is adopted for modeling beam-soil interaction, focusing on:
Dynamic response
Natural frequencies
Mode shapes
Based on the approach by Zhemochkin and Sinitsyn, and extended by Guenfoud et al.
The semi-analytical method provides a balance between accuracy and computational efficiency for modeling dynamic soil–beam systems.
Green’s functions are central to capturing how soil reacts to forces.
The system is modeled as discrete elements rather than a continuous contact surface.
The approach enables evaluation of vibration modes, natural frequencies, and responses to loads, incorporating soil inertia but neglecting damping and friction.
Conclusion
This paper presents a semi-analytical approach for the dynamic analysis of beams resting on an elastic half-space with inertial properties. The method is developed to evaluate the Eigenfrequencies, mode shapes, and dynamic response of the beam subjected to external excitation.
The analysis considers the interaction forces within the contact zone, which are essential for determining other physical parameters. The formulation is based on Green’s functions to represent the displacements in the contact region. Incorporating the inertial effects of the half-space, also known as Lamb’s problem, introduces mathematical challenges, particularly due to instabilities associated with hypergeometric functions. The present work overcomes these limitations and provides stable, reliable expressions for the solution. Results show that the computed mode shapes of the beam are irregular, and the proposed semi-analytical method is both computationally efficient and practical for engineering applications. Furthermore, the approach can serve as an approximation function in numerical models, enabling its application to more complex problems of wave propagation and dynamic soil–structure interaction.
References
[1] Lamb H. On the propagation of tremors over the surface of an elastic solid. Philosophical Transactions of the Royal Society of London 1904;A 203:1–42.
[2] Graff Karl F. Wave motion in elastic solids. Courier Dover Publications;1991.
[3] Guenfoud S, Amrane MN, Bosakov SV, Ouelaa N. Semi-analytical evaluation of integral forms associated with Lamb’s problem. Soil Dynamics and Earth- quake Engineering 2008;29;438-43.
[4] Guenfoud S, Bosakov SV, Laefer DF. Dynamic analysis of a beam resting on an elastic half-space with inertial properties. Soil Dynamics and Earth- quake Engineering 2009; 29; 1198-207
[5] Reissner E. Stationare, axialsymmetrische, druch ein schuttelnde Masse erregte Schwingungeneines homogenen elastische Halbraumes. Ingenieur Archiv 1936; 7(6):381–396.
[6] Shekhter O.Y. The taking into account of the inertial properties of the ground in the calculation of a foundation unit on forced vibration. The work’s collection of N.I.I.O., vibration of constructions and foundations 1948;(2):p.72–89;
[7] Bycroft G.N. Forced vibrations of a rigid circular plate on a semi infinite elastic space and on an elastic stratum. Philos Trans R Soc Lond 1956; 248(A.948):327–68.
[8] Pekeris CL. The seismic buried pulse. Proceedings of the National Academy of Sci USA 1955;41:469–80.
[9] Guan F., Novak M. Transient response of an elastic homogeneous halfspace to suddenly applied rectangular loading. Journal of Applied Mechanics: ASME 1994;61:256–263.
[10] Shinozuka M., Deodatis G., Zhang R., Papageorgiou A.S. Modeling, synthetics and engineering applications of strong earthquake wave motion. Soil Dynamics and Earthquake Engineering 1999;18:209–228.
[11] Jin B., Liu H. Exact solution for horizontal displacement at center of the surface of an elastic half space under horizontal impulsive punch loading. Soil Dynamics and Earthquake Engineering 1999;18(7): p.495–498.
[12] Jin B. Elastic halfspace under impulsive, distributed, vertical loading at surface: exact solution at center for a punch-like distribution. Soil Dynamics and Earthquake Engineering 1998;17(5):311–316.
[13] Jin B. Exact solution for horizontal displacement at center of the surface of an elastic halfspace under horizontal impulsive punch loading. Soil Dynamics and Earthquake Engineering 1999;18:495–508.
[14] Zhou X.L., Wang J.H., Lu J.F. Transient foundation solution of saturated soil to impulsive concentrated loading. Soil Dynamics and Earthquake Engineering 2002;22:273–281.
[15] Pradhan P.K., Baidya D.K., Ghosh D.P. Dynamic response of foundations resting on layered soil by cone model. Soil Dynamics and Earthquake Engineering 2004;24:425–434.
[16] Gazetas G. Formulas and charts for impedances of surface and embedded foundations. Journal of Geotechnical Engneering ASCE1991;117(9):1363–81.
[17] Zhemochkin BN, Sinitsyn AP. Practical methods of the calculation of beams and plates resting on an elastic foundation. Moscow: Stroyizdat publishing company; 1962 .
[18] Johnson K. Contact interaction’s mechanics. Moscow: Mir publishing company; 1989 (Russian Edition).
[19] Georgiadis HG, Vamvatsikos D, Vardoulakis I. Numerical implementation of the integral-transform solution to Lamb’s point-load problem. Comput Mech 1999; 24:90–9.
[20] Zeng S, Liang R. A fundamental solution of a multilayered half-space due to an impulsive ring source. Soil Dyn Earthquake Eng 2002; 22(7):541–50
[21] Zimmermann, H., 1888. The calculation of the railway superstructure. Ernst & Korn, Berlin.
[22] Winkler, E., 1867. The theory of elasticity and strength with special reference to their application in the technology. Dominicus, Prag.
[23] Boussinesq (1885). Potential application of the study of equilibrium and motion of elastic solids.
[24] Nam-Gyu Park, Sang-Youn Jeon, Shin-Ho Lee, Kyeong-Lak Jeon, Joon-Ro Lee, Kyu-Tae Kim Mark. E. Dye. Dynamic Analysis of a Beam with Discretely Spaced Elastic Supports. Journal of Mechanical Science and Technology ; Volume 18 Number 5 (2004), 733-741
[25] J. C. O. Nielsen. Eigenfrequencies and Eigenmodes of a beam structure on an elastic foundation. International Journal for Numerical Methods in Engineering. 1992; 33: 1843-1859