A novel theoretical model is proposed to investigate the torsional response of a pile in fractional-order viscoelastic unsaturated transversely isotropic soil with imperfect contact. This model employs Biot's framework for three-phase porous media along with the theory of fractional derivatives. Unlike previous models that assume continuous displacement at the pile-soil interface, this study uses the Kelvin model to simulate relative slippage between pile-soil contact surfaces (imperfect contact). Incorporating fractional-order viscoelastic and transversely isotropic models to describe the stress-strain relationship, comprehensive dynamic governing equations are derived. Using the separation of variables method, inverse Fourier transform, and convolution theory, analytical solutions for the frequency domain response and semi-analytical solutions for the time domain response of the pile head under semi-sine pulse excitation are obtained. Using numerical examples, the effects of model parameters in the fractional-order viscoelastic constitutive model, pile-soil relative slip and continuity model, and soil anisotropy on the torsional complex impedance, twist angle, and torque are presented.

Geo-materials naturally display a certain degree of anisotropy due to various effects such as deposition. Besides, they are often two-phase materials with a solid skeleton and voids filled with water, and commonly known as poroelastic materials. In the past, despite numerous studies investigating the vibrations of strip foundations, dynamic impedance functions for multiple strip footings bonded to the surface of a multi-layered transversely isotropic poroelastic half-plane have never been reported in the literature. They are first presented in this paper. All strip foundations are assumed to be rigid, fully permeable, and subjected to three types of time-harmonic loadings. The dynamic interaction problem is investigated by using an exact stiffness matrix method and a discretization technique. The flexibility equations are established by enforcing the appropriate rigid body displacement boundary conditions at each footing-layered soil interface. Numerical results for dynamic impedance functions of two-strip system are presented to illustrate the influence of various effects on dynamic responses of multiple rigid strip foundations.