Abstract
This study presents an analytical solution to examine the mechanical behavior of an incompressible, functionally graded hyperelastic cylinder under combined extension and torsion. The exp-exp strain energy density function characterizes the hyperelastic material, with parameters varying exponentially along the radial direction. To validate the solution, finite element simulations using a custom UHYPER in ABAQUS are performed. The analytical and numerical results show strong agreement across different stretch and twist levels. The stress distribution and maximum stress are significantly influenced by the exponential parameter governing material gradients. Unlike axial stretch, torsion induces a more intricate longitudinal stress distribution, with large twisting producing two extrema that shift toward the cylinder's center and outer surface. Longitudinal stress primarily governs von Mises stress and strain energy density variations across the radial direction. A critical axial stretch is identified, below which torsion-induced axial force transitions to compression, elongating the cylinder during twisting. Beyond this stretch, the axial force shifts from tensile to compressive with increasing twist, causing initial shortening before further elongation.