Saint-Venant torsion of cylindrical orthotropic bar

The Saint-Venant torsion of cylindrical orthotropic homogeneous linearly elastic bar is considered. The cross section of the bar is bounded by two straight lines and a curved arc. The geometry of the boundary curve of cross section depends on the shear moduli of the cylindrical orthotropic bar. An analytical method is presented to obtain the Prandtl’s stress functions, warping function, shearing stresses and the torsional rigidity.


Introduction
The Saint-Venant torsion problem for elastic bars has been addressed by engineers for the last 150 years. The uniform torsion was described by the rigorous theory developed by Saint-Venant [1,2]. According to the Saint-Venant's theory of uniform torsion only two shearing stresses appear nonzero, the axis of the bar remains straight and the cross-sectional warping does not depend on the position of the cross section. These conditions are valid for anisotropic bars if the one plane of elastic symmetry is orthogonal to the axis of bar [3][4][5][6][7][8].
The Saint-Venant torsion of anisotropic linearly elastic bars has been a subject of several studies. Books by Lekhnitskii [3,4], Sarkisyan [5,6], Rand and Rovenski [7], Milne-Thomson [9] give the detailed analyses of Saint-Venant torsion of anisotropic and of orthotropic bars. In these books, both the Prandtl' stress function and torsion function formulations are presented. The warping properties of twisted functionally graded Cartesian anisotropic linearly elastic bars are treated by Horgan [10]. Some important properties of torsion of Cartesian anisotropic bar for elliptical cross section have been studied in paper by Ecsedi and Baksa [11]. Paper [12] deals with the Saint-Venant's torsion of elastic cylindrical orthotropic solid elliptical cross section. By the use of principle of minimum of potential energy and principle of minimum of complementary energy, approximate analytical solutions are derived for the torsion function and for the Prandtl's stress function in paper by Ecsedi and Baksa [12].
Papers [13][14][15][16] present the formulation of uniform torsion and some analytical solutions for the cylindrical orthotropic elastic bars. The Saint-Venant torsion of the compound cylindrical orthotropic bar is studied in paper [17]. Analytical solution is given for the uniform torsion of the cylindrical orthotropic annular wedge-shaped bar whose curved boundary segments are strengthened by thin isotropic elastic shells (Ecsedi and Baksa [17]). The paper [18] deals with the Saint-Venant torsion of elastic cylindrical orthotropic composite bar whose cross section is a sector of a solid circle. An analytical method is presented to obtain the Prandtl's stress function and the torsion function of the composite cross section made of two different cylindrical orthotropic materials.
The present paper addresses the Saint-Venant torsion for cylindrical orthotropic homogeneous linearly elastic bars, whose cross section is bounded by two straight lines and a curved arc. The geometry of the boundary curve of solid cross section depends on the shear moduli of the cylindrical orthotropic bar. The cross section of the considered bar has an axis of symmetry (Fig. 1). The polar coordinate system Orϕz is positioned at the left and cross section of the bar as shown in Fig. 1. The unit vector in radial direction is e r , in circumferential direction the unit vector is denoted by e ϕ , and the unit vector in axial direction of the bar is e z . The cross section of the bar is bounded by two straight lines c 1 and c 2 and a curved arc c 3 .
The geometry of the cross section depends on the material properties of cylindrical orthotropic elastic material. The principle directions of cylindrical orthotropy at point P are given by the unit vectors e r , e ϕ and e z , and the shear moduli are [3][4][5][6] Let The positions of boundary line segments c 1 and c 2 are given by the angle α ( Fig. 1) and the equation of curved boundary curve c 3 is Since α < π/2, we have The shape of boundary curve segment c 3 as a function of β is shown in Fig. 2.

Governing equations
The torsion function is denoted by ω = ω(r, ϕ), the Prandtl's stress function of the considered triangle-like shape cross section is U = U (r, ϕ). The rate of twist is represented by ϑ, and the applied torque is T . The torsional rigidity S is defined as It is known that the Saint-Venant torsion of cylindrical orthotropic homogeneous linearly elastic bar leads to the following boundary value problem for the torsion function G r z n r ∂ω ∂r In Eqs. (7), (8), A denotes the cross section, ∂ A is the boundary curve of A; n r and n ϕ are the component of the boundary normal vector n = n r e r + n ϕ e ϕ [4][5][6][7]. The solution of Saint-Venant uniform torsion in terms of Prandtl's stress function obtained from the following Dirichlet type boundary value problem Knowing solutions of boundary value problems given by Eqs. (7), (8) and Eqs. (9), (10) the shearing stresses τ r z = τ r z (r, ϕ) and τ ϕz = τ ϕz (r, ϕ) can be computed from formulae (11) and (12) The torsional rigidity of the cross section can be obtained by the use of Eq. (13) or Eq. (14)

Solution of the Saint-Venant torsion by Prandtl's stress function
We shall try to solve the Dirichlet type boundary value problem formulated in Eqs. (9) and (10) for some computable constant C. The form of the function U = U (r, ϕ) is suggested by the homogeneous boundary condition (10). After some simple calculations from Eq. (9), it follows that that is, The shearing stresses can be computed by the application of formulae (11) and (12). We have From the formula (14), we obtain for the torsional rigidity Equation (20) for β = 1 (G r z = G ϕz ) gives the expression of torsional rigidity of isotropic equilateral triangle cross section [8]. The shearing stresses τ r z and τ ϕz are zero at point P (2a/3, 0), i.e., independently of the value of β. The torsional rigidity of the cross section in terms of cross-sectional area can be expressed as

Solution of the Saint-Venant torsion by torsion function
To obtain the torsion function ω = ω(r, ϕ), the following equations will be used The validity of these equations follows from the formulae (11) and (12). Solution of the system of partial differential equations (29) for ω = ω(r, ϕ) is as follows ω(r, ϕ) = 1 2aβ r 2 (4r cos 2 βϕ − 4a cos βϕ − r ) sin βϕ By substitution of the expression of torsion function ω = ω(r, ϕ) into Eqs. (7) and (8), it can be pointed out that the function given by Eq. (25) is a solution of the Neumann type boundary value problem formulated in Eqs. (7) and (8). In the checking computation it is used that, the components of the boundary curve normal vector are A rigid body like translational motion of the bar in direction of axis z is described by c. We note that the warping of the cross section depends on only one material parameter β.

Numerical example
The following data are used in the numerical example: The contour lines of the Prandtl's stress function are shown in Fig. 3. In Fig. 4, the contour lines of torsion function for c = 0 are given. The radial shearing stress on the boundary line segment c 1 is shown in Fig. 5 as a function of the radial coordinate r .
The circumferential shearing stress as a function of r at ϕ = 0 is illustrated in Fig. 6. In Fig. 7, the plots of torsion function are shown for five different values of radial coordinates (r = The contour lines of resultant shearing stress are shown in Fig. 8. The expressions of the τ (r, ϕ) in terms of T at r = a, ϕ = 0 and r = a √ 3 , ϕ = ± π 6β are as follows It is known that for β = 1 (32)

Conclusions
The object of this paper is the Saint-Venant torsion of cylindrical orthotropic homogeneous linearly elastic bar with solid cross section. The cross section of the considered bar is bounded by two straight lines and a curved arc. The geometry of the boundary curve of cross section depends on the shear moduli of the cylindrical orthotropic material. Closed-form solutions were derived for the Prandtl's stress function, torsion function, shearing stresses and torsional rigidity. An example illustrates the application of the obtained formulae.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Funding Open access funding provided by University of Miskolc.

Declarations
Conflict of interest On behalf of all authors, the corresponding author states that there is no conflict of interest.