Abstract
The problem of repeated imposition of large tension-compression and torsion deformations is solved for a sequentially formed multilayered cylinder, each layer of which is made of a viscoelastic incompressible material. After attaching each successive layer, tension-compression and torsion of the composite cylinder is performed. Based on the description of the formation and deformation of the cylinder, a mathematical formulation of the problem is given for the quasi-static case and the main steps of its exact analytical solution are described. The results of numerical calculations are presented.
Keywords
1. Introduction
In nonlinear mechanics of a deformable solid, great importance is attached to the study of effects that arise during torsion of cylindrical bodies. One of these effects is the Poynting effect [1,2], which consists in changing the length of the cylinder during torsion in the absence of axial forces applied to the bases of the cylinder. Another manifestation of this effect is the occurrence of stresses directed along its axis during torsion of the cylinder, when the bases of the cylinder cannot move along this axis.
When solving the problem of torsion of a cylinder made of a nonlinear elastic material, exact analytical solutions are of great importance. For incompressible isotropic nonlinear elastic materials, an exact solution to the problem of torsion and tension-compression of a circular cylinder is known, which is one of the universal solutions [2–9]. Torsion of cylinders of non-circular cross-section under large deformations is studied by Zubov and Bogachkova [10], in which, in particular, approximate analytical solutions for cylinders of elliptical and rectangular cross-sections were obtained. The influence of surface effects on the stress-strain state of rods during torsion and tension-compression was studied by Sigaeva and Czekanski [11].
For compressible nonlinear elastic materials, solving the problem of torsion of rods is more complicated. Exact analytical solutions are known for the case of pure torsion for some special classes of materials [12]. Approximate analytical solutions using series were proposed by De Pascalis [13]. Dell’Isola et al. [14,15] applied the Signorini method to the approximate calculation of the stress-strain state of pre-deformed rods. Hajhashemkhani and Hematiyan [16] formulated the problem of inflation, stretching and torsion of a functionally gradient compressible hyperelastic tube as a one-dimensional problem and solved by the finite difference method.
Effects similar to the Poynting effect have been revealed at large deformations for inelastic materials. In plasticity theory, this effect is known as the “Swift effect” [17]. An experimental and numerical study of this effect for polymers in a glassy state was carried out by Wu and van der Giessen [18]. Bruhns et al. [19] obtained an approximate analytical solution to the problem of torsion of a circular cylinder made of an elastoplastic material within the framework of a plasticity model with isotropic hardening. Cazacu et al. [20] studied within the framework of the elastic-plasticity model, the Swift effect during cyclic deformation of a sample. Colak and Krempl [21] simulated the Swift effect in a hollow cylinder made of viscoplastic material, including the case of cyclic deformation. Khan et al. [22] studied the torsion of a polymer material within the framework of nonlinear viscoelasticity models. The combined action of high pressure and torsion is used to create metals with very small grain sizes [23]. To model this phenomenon, both numerical methods and approximate analytical solutions are used [24,25]. An analytical solution to the problem of large elastoplastic deformations of a cylinder under the combined action of pressure and torsion was obtained by Sevastyanov [26].
When studying joint torsion and tension-compression, problems of loading of bodies in several stages are of interest, when at each stage a new part of the body is attached to the loaded body. These problems belong to the class of problems on the surface growth of deformable solids. A detailed literature review of the problems on surface growth is given by Sozio and Yavari [27]. The models of surface growth (accretion) include the models of continuous growth and the models of discrete growth [28–31]. If the deformations are sufficiently large, then the theory of multiple superposition of large deformations is used when these problems are formulated and solved [32–34]. Zingerman and colleagues [35,36] considered the problems of tension-compression and torsion of composite cylinders made of nonlinear elastic materials. These problems were solved for the case of discrete growth when one or more cylindrical layers were attached to the cylinder. Similar problems have been solved for hypoelastic materials [37].
The continuous surface growth of nonlinear-elastic circular cylindrical bars subjected to finite torsion is analyzed by Yavari and Pradhan [38]. The residual stresses are computed, and nonlinear effects are estimated. The surface growth of cylindrical bars, growing continuously in radial direction and subjected to axial stretch or the action of axial force, is investigated by Yavari et al. [39]. In these papers, it is assumed that strains are finite. For viscoelastic bodies the models of continuous accretion under finite and small strains are developed by Drozdov [40]. In particular, the continuous growth of viscoelastic bars in this book is analyzed under small torsion strains.
The models of surface accretion can be used for the description of growth of biological tissues [41,42] and for the strength analysis of additive manufacturing products [43–45].
In the current work, the discrete accretion of a composite viscoelastic cylinder subjected to tension/compression and torsion is analyzed under finite strains. The cylinder is considered as a body subjected to the imposition of large deformations. The mechanical properties of the cylinder material are described by constitutive relations of integral type, which, in the absence of viscosity, are reduced to the constitutive relations for a neo-Hookean material [46,47].
2. Notation
In this work, the notation of the theory of repeated superposition of large deformations is used [32–37]. The states (configurations) of the body are numbered from 0 (initial state) to n (final state).
Relationships between cylindrical and Cartesian coordinates:
3. The mechanical problem statement
The problem is formulated within the framework of the theory of repeated superposition of large strains [32–35]. There is a cylindrical bar made of incompressible isotropic viscoelastic material. Initially, there are no stresses and deformations in this bar. Then, at a given time interval, the bar undergoes an initial torsional deformation together with tension or compression in the axial direction under the action of a torque and axial force. This deformation can change over time according to a given rule. The lateral surface of the cylinder is unloaded. Due to the incompressibility of the material, tension or compression of the rod along its axis is accompanied, respectively, by compression or tension of the rod in the radial direction. It is assumed that at the end of the given time interval the bar is in the first intermediate state. In this state, a hollow undeformed viscoelastic cylinder (layer) joins it. The inner radius of this hollow cylinder coincides with the radius of the bar after the initial deformation, and the length coincides with the length of the bar after the initial deformation. The cylinders are glued along their common boundary surface (for the initial cylinder this is the outer and only boundary, for the attached cylinder—the inner one), and in what follows it is assumed that the conditions of ideal contact are satisfied on this surface.
Furthermore (during the next time interval), there is a joint additional deformation of the cylinders. It is also a deformation of torsion together with tension or compression along the axis, which occurs as a result of the action of a torque together with an axial force. At the end of the second time interval, the compound cylinder passes to the second intermediate state. Then the next hollow viscoelastic cylinder (layer) is attached to it. This layer is not initially deformed.
This procedure of sequential connection of the cylinders with their subsequent deformation of torsion and tension-compression can be repeated several times.
The solution is obtained for the quasi-static processes. It is assumed that each layer is added instantly, in sense of a study by Christensen [48], and the dynamical effects are not taken into account.
4. Calculation of strains and stresses in a composite cylinder
We consider a multilayered circular cylinder deformed over several stages. The mth stage (
At the beginning of the
The coordinates of particles of the cylinder during the mth stage of deformation are related to the coordinates of these particles in the previous state as follows [36,37]:
where
Considering relations (1) for several successive stages, we obtain:
where
To characterize the deformation of the lth layer during the mth stage, the deformation gradient is used:
The Green deformation tensor during transition from state l to m is computed as follows:
It is assumed that the material of the multilayered cylinder is viscoelastic. The generalized stress tensor satisfies the constitutive relations:
where μ, η, β are the mechanical parameters of the material,
The true stress tensor is related to the generalized stress tensor as follows:
When writing formula (5), the incompressibility of the material is taken into account.
Using formulas (3)–(5), we can write down the stresses
We write the equilibrium equations for the mth state in the form:
Differential equation (6) is supplemented with boundary conditions. For the outer layer
For the internal layers the boundary conditions are written as follows:
By solving equation (6) with boundary condition (7) or (8) for each
The axial force and torque in the mth state are computed using the relations [29–31]
For the verification of the obtained results, the computations were performed for the case of pure elastic material [
5. Numerical studies
Calculations were performed for the case of
The formulas for the deformation parameters have the form
where
The distribution of true stresses in the cylinder at the final moment of time

Distribution of

Distribution of

Distribution of

Distribution of

Distribution of

Distribution of

Distribution of

Distribution of
The dependences of hoop stresses at the outer boundary of the cylinder on the deformation parameters K and L at various times are shown in Figures 9 and 10. The values of time, referred to T, are given in the graph legends.

Dependence of

Dependence of
The dependences of axial force and torque on the deformation parameters K and L at various times are shown in Figures 11–14. The values of time, referred to T, are given in the graph legends.

Dependence of

Dependence of

Dependence of

Dependence of
One can see from Figures 9–13 that significant nonlinear effects take place for the given values of parameters. The dependence of torque on the deformation parameter K is linear (Figure 14). Figures 9–14 show that the hoop stress, axial force, and torque are decreased in time by absolute value. This effect is related with the relaxation processes.
6. Conclusion
An exact analytical solution is constructed for the problem of tension-compression and torsion of composite viscoelastic cylinder undergoing large strains. It is assumed that the cylinder is composed of some layers and these layers are added successively. The problem is formulated using the theory of multiple imposition of large strains. The results of numerical studies are obtained. Significant nonlinear effects are found. The effects of viscosity are analyzed. The results of solution can be used for the verification and validation of software packages for engineering analysis. For example, this solution can be used to verify calculations of the stress–strain state of structural elements manufactured using additive technologies [45,51]. It is of interest to take into account surface effects [52]. It is possible to apply the proposed approach to modeling the surface growth of deformable bodies [28,29,31] and to modeling phase transitions [53].
Footnotes
Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was completed at Moscow State University and supported by the Russian Science Foundation (project no. 22-11-00110).
