key: cord-1023052-z46oz1qm authors: Postolaki, Lesya; Tokovyy, Yuriy title: Axisymmetric residual stresses in a solid cylinder of finite length date: 2022-05-18 journal: J Eng Math DOI: 10.1007/s10665-022-10221-y sha: 4e9610ef8178fab8f4dd6af788fea2b394826b9c doc_id: 1023052 cord_uid: z46oz1qm A solution technique is presented for the determination of residual stresses in a finite-length solid cylinder subject to non-uniform axisymmetric distributions of incompatible residual strains. The problem is reduced to the sequential solving of three individual problems: a problem on the determination of residual stresses in an infinitely long cylinder (the basic state) and two auxiliary problems for evaluating the stresses induced by the end-face effects (the disturbed states). The variational method of homogeneous solutions is implemented in order to determine the disturbed states within the framework of two latter problems. The solution technique is verified numerically for typical distribution profiles of incompatible strains depending on both the radial and axial coordinates. The approach can be used to evaluate residual stresses in a solid cylinder of finite length due to the intense thermal treatment. The practical measurement of the residual strains or stresses in solids is typically based on either destructive or non-destructive testing. The destructive techniques imply the complete or partial destruction of structural members in order to perform the measurements (e.g., the hole drilling and ring coring techniques [5] [6] [7] , the slitting and contour methods [8] [9] [10] , etc.). On the other hand, the non-destructive methods (e.g., the X-ray and neutron diffraction [11] [12] [13] , the magnetic, ultrasonic, and optical methods [14] [15] [16] , etc.) allow for saving the structural integrity but have limited applicability. The most advanced way for predicting the actual levels of residual stresses corresponding to certain distribution profiles of incompatible strains combines both theoretical and experimental techniques [17, 18] . One of such approaches is a non-destructive computational-experimental method [19] based on solving the inverse elastoplastic problems [20] and fitting the results with the experimental evidence obtained by non-destructive techniques in order to evaluate the effective stress distribution parameters. The formulation and solution of the inverse problems is based on the so-called conventional-plastic strain hypothesis [21, 22] implying the base material of a solid to be elastic, while the expected zone of residual stress distribution exhibits specific elastic-plastic behavior so that the total strain can be represented in the following form: where ε is the total strain tensor, ε e is the tensor of elastic strains occurring within the entire solid, and ε 0 denotes the tensor of incompatible strain [23] distributed within the area affected by the impacts causing the residual stressstrain state. This approach has been efficiently employed for the analysis of residual stresses due to welding residual strains in infinite layers with rectilinear and circular welds [24, 25] , butt-welded cylindrical vessels [26, 27] , and rectangular plates [28] . Special attention in the latter paper was given to the end effects in the butt-welded rectangular plate by performing the exact analysis based on the direct integration method [29] . It was shown, in particular, that the presence of boundary that is transverse to the butt-weld axis effects significantly the aggregated stress state comparing with one for an infinitely long (with respect to the butt-weld length) plate. Thus, the end effects are very important for the accuracy in the evaluation of the residual stress-strain state of the bounded domains. In general, the exact analysis of the end effects in the stress state of finite solids with irregular points (lines) of the boundary (i.e., the corner points or the edge lines) presents a challenge for both analytical and numerical methods due to the fact that the resulting operator of the corresponding boundary value problem of the elasticity theory is not self-conjugated. This fact complicates significantly the variables separation in the governing equations and, consequently, exact satisfaction on the boundary conditions imposed on the dissimilar segments of the boundary separated with the corner points. An interested reader can find more detail on the earlier development history of the methods for two-dimensional problems of such kind presented by Meleshko [30] . In the case of general boundary conditions, the dominant analytical methods for solving the boundary value problems for finite domains with corner points imply the realization of one of the following strategies: A) the construction of general solutions to the governing equations which allow for approximate satisfaction of boundary conditions imposed on some or all sides of a finite domain, and B) the construction of a solution that exactly satisfies the full set of boundary conditions imposed, while the governing equations are satisfied within the required accuracy. Strategy B was realized, for example, with the method of direct integration [29, 31] , which has been recently used in [32] for solving a thermoelasticity problem in a cylinder of finite length. Strategy A is represented by at least two dominant methods, i.e., the method of cross-wise superposition [33, 34] and the method of homogeneous solutions [35] . The latter method presents a natural extension of the classical separation of variables implying the construction of eigenvalues, which belong to the complex plane C in the case of finite domains with corner points. The corresponding eigenfunctions are usually non-orthogonal which complicates the satisfaction of boundary conditions. In order to overcome this difficulty, the variational approaches with implication of the least square method are often employed. In [36] , a variational method of homogeneous solutions was used for the evaluation of residual stresses in the vicinity of an inplane weld joint connecting the dissimilar materials. This method was also used for the solution of problems on the non-destructive evaluation of the axisymmetric residual stresses near circular welds in long cylindrical shells on the basis of the available empirical data obtained through the magnetoelastic method [37] . 4 A variational method of homogeneous solutions was used in [38, 39] for the analysis of axisymmetric residual stresses in a finite-length cylinder. The residual stresses are modeled with the use of concept of incompatible residual strains (1) with the latter strains being regarded as the cause of the residual stresses appearance. The incompatible strains were presented by an isotropic tensor depending only on the radial coordinate. This method was also used for a variational formulation of the inverse problem for the residual stresses determination on the basis of empirical data obtained by a photoelasticity method. This paper is aimed to extend the technique on the evaluation of the axisymmetric residual stresses in a finite-length solid cylinder due to the incompatible strains depending on both the radial and axial coordinates. This provides a prominent tool for the analysis of the end effects in the cylinders of finite length. Consider an elastic cylinder of finite length 2H and radius R. In the dimensionless cylindrical coordinate system, the cylinder occupies domain Here, r = ρ/R, z = ζ /R, b = H/R, ρ and ζ are the dimensional radial and axial coordinates, and θ is the circumferential coordinate. Assume the lateral surface and the end faces of the cylinder to be free of the external force loadings, i.e., Here, σ rr (r, z), σ zz (r, z), and σ r z (r, z) are the normal and tangential stress-tensor components. Assume the interior of the cylinder to undergo a non-uniform axisymmetric (i.e., irrespective of coordinate θ ) distribution of incompatible [23] strains ε 0 rr (r, z), ε 0 θθ (r, z), and ε 0 zz (r, z) induced by thermal treatment and, hence, represented by the normal components only [40] . Then, according to the model (1), the components ε rr (r, z), ε θθ (r, z), ε zz (r, z), and ε r z (r, z) of the total strain are to be represented in the form as follows: Here, ε e rr (r, z), ε e zz (r, z), ε e θθ (r, z), and ε e r z (r, z) are the elastic strains. The total strains are expressed through the radial and axial displacements, u r (r, z) and u z (r, z), via the straindisplacement Cauchy equations [41] : Within the context of equations (3) and (4), the constitutive equations of Hooke's law [40, 41] take the following form: where ε = ∂u r /∂r + u r /r + ∂u z /∂z, ε 0 = ε 0 rr + ε 0 θθ + ε 0 zz , and λ and μ are the Lamé constants. Assuming the cylinder to be in the state of equilibrium, the stress-tensor components meet the following equilibrium equations [40, 41] 1 r Substituting equations (5) into the equilibrium equations (6) yields the following Lamé-Navier equations: where For the case of isotropic tensor of incompatible strains, i.e., when ε 0 rr = ε 0 θθ = ε 0 zz = ε 0 /3, Eqs. (7) take the following form: (5) and (8) present a mathematical model (a closed-form boundary value problem) for the determination of axisymmetric residual stresses and displacements in cylinder C due to non-uniform distributions of incompatible strains represented with the isotropic tensor induced by intensive thermal treatment depending on both the radial and axial coordinates. Below, we provide a solution technique for the formulated problem based on the variational method of homogeneous solutions. In our approach, we follow the solution strategy [39] by presenting the displacement field as a superposition of two states, i.e., Here,û r (r, z) andû z (r, z) are the displacements resulted from non-homogeneous equations (8) for an infinitely long cylinder (b 2 → ∞) of the unit radius (they reflect the so-called basic state of the cylinder), andũ r (r, z) and u z (r, z) are the displacements derived from the corresponding homogeneous equations formulated for cylinder C . Here,ε = ∂ũ r /∂r +ũ r /r + ∂ũ z /∂z. The latter displacements reflect the disturbance of the deformed state due to the end effects and thus are referred to as the "disturbed state." In order to construct a solution to the system of non-homogeneous equations (8) for the basic state, we implement the following representation [33] : where (r, z) is an unknown potential function that is analogous to the Papkovich -Goodier thermoelastic potential [42, 43] . Substitution of representations (10) into (8) yields A solution to Eq. (11) can be given in the form as follows [44] : where , which can be evaluated by making use of the following formula: Thus, the displacementsû r (r, z) andû z (r, z) can be computed by putting (12) into (10) . The corresponding stresstensor components,σ rr (r, z),σ θθ (r, z),σ zz (r, z), andσ r z (r, z), can be computed by making use of (5), (10) , and (11) in the following form: In such manner, the stresses and displacements are determined for the basic state of an infinitely long solid elastic cylinder of the unit radius under the isotropic incompatible strains distributed non-uniformly within the coordinates (r, z). These components induce the non-uniform boundary conditions on the lateral surface and the end faces of the finite cylinder C . In order to compensate these non-uniform boundary conditions, we introduce the "disturbed state," which can be determined by implementing the technique presented in the following section. The disturbed state is to be realized by solving the following two boundary value problems. The first problem implies solving system (9) for cylinder C with free lateral surfacẽ and the end faces loaded as follows: The second boundary value problem implies solving the system (9) accompanied with the boundary conditions of free end faces σ (2) zz (r, ±b) = 0,σ (2) r z (r, ±b) = 0 along with the conditions of loaded lateral surfacẽ The total components for the disturbed state can be computed by superposing the solutions of boundary value problems (9), (13) , and (14) , and (9), (15) , and (16) . A general solution to system (9) can be constructed by implementing the Love stress function χ(r, z) [45] satisfying the following biharmonic equation The displacements are represented via the Love function by the following expressions: and the stress-tensor components are expressed as follows: Here, ν is the Poisson ratio. In view of the linearity of the problem, we can construct the Love functions χ (1) (r, z) and χ (2) (r, z) corresponding to the boundary value problems (9), (13), (14) and (9), (15) , (16) , separately; then they can be superposed in order to determine the disturbed stress-strain state of the finite cylinder C . Assuming, for the simplicity sake, the symmetry of the problem with respect to the plane z = 0, we can represent function χ (1) (r, z) in the following form [46] : where B k are unknown complex constants, the overline denotes the complex conjugation, γ k are the complex roots of the transcendental characteristic equation Within the context of expressions (17) and (18), the stress field for the first disturbed problem (9), (13) , and (14) can be given in the form as follows: θθ,k (r, z) , where σ (1) rr,k (r, z), σ (1) θθ,k (r, z), σ (1) zz,k (r, z), and σ (1) r z,k (r, z) are given in Appendix A. Expressions (19) allow for exact satisfaction of the equations in system (9) and boundary conditions (13) under arbitrary constants B k , which are to be found by satisfying conditions (14) . The biharmonic function for the second disturbed problem χ (2) (r, z) can be constructed in the form as follows [47] : where , and λ k are the complex roots of the transcendental characteristic equation sinh (2λ k b) + 2λ k b = 0. Now, using (17) and (20), we can construct the stresses for the second disturbed problem (9), (15) , and (16) in the following form: rr,k (r, z) +Ā kσ (2) rr,k (r, z) + c, θθ,k (r, z) , zz,k (r, z) , Here, functions σ (2) rr,k (r, z), σ (2 θθ,k (r, z), σ (2) zz,k (r, z), and σ (2) r z,k (r, z) are given in Appendix B, and c is an arbitrary constant. Representations (19) of the stress-tensor components meet the boundary conditions (13) under arbitrary coefficients B k . Meanwhile, representations (21) meet the boundary conditions in (15) under arbitrary coefficients A k . Obviously, the sets of coefficients B k and A k , k = 1, 2, ..., can be determined by satisfying conditions (14) with expressions (19) and conditions (16) with expressions (21) , respectively. In order to do so, we implement the variational method of homogeneous solutions [39] by introducing the following quadratic functionals corresponding to the first and second disturbed problems: The functionals in (22) and (23) can be minimized under the following necessary conditions [48] : The last equation in (24) yields the following formula: rr,k (r, z) +Ā kσ (2) rr,k (r, z) dz. Within the context of the latter formula along with the representations (14), (16) , (19) , and (21) the remaining equations in (24) yield the following infinite systems of linear algebraic equations: where In view of the zero asymptotic of the coefficients in the systems given by equation (25) at k → ∞, the practical computation of the parameters A p k and B p k along with their complex conjugation can be performed by implementing the simple reduction algorithm [48] . The total values of the stress-tensor components of the first and second disturbed problems and the stress components of the initial state are the solution of the formulated direct problem of determining the residual stresses. Consider a numerical example of the isotropic tensor of residual strain given by formula Here, Fa(a, x) = (tanh(ax) + 1)/2, a = {a r , a z }, e l , r l , and b 1 are given constant parameters, l = 1, 2. Figure 1 presents the distribution of the residual strain (26) computed under the fixed parameters e 1 = 1, e 2 = −1, r 1 = 1/3, r 2 = 2/3, b = 1, a r = 18, a z = 10, and four different values b 1 = 1.00; 0.75; 0.50; 0.25. Note that the strain ε 0 given in the form (26) varies within both the radial and axial coordinates, r and z. It reflects extension on the cylinder axis, which is maximum at z = 0 and decreases when approaching the end faces z = ±b. On the lateral surface of the cylinder, the strain (26) verbalizes contraction which decreases when approaching the end faces. For smaller positive values of b 1 , the decrement occurs faster. The area r 1 < r < r 2 , whose width is defined by r 2 − r 1 , is not disturbed with the residual strain (26) . Thus, by controlling the parameters involved in (26) , it is possible Fig. 1 Full-field distributions of the residual strain (26) in a "cubic" cylinder (the radius equals the half-length) at e 1 = 1, e 2 = −1, r 1 = 1/3, r 2 = 2/3, b = 1, a r = 18, a z = 10, and b 1 = 1.00; 0.75; 0.50; 0.25 to control the distribution profile and the magnitude of the residual stress ε 0 . Note that decreasing the value of b 1 allows for reflecting rather "localized" effect of the residual strain (26) at a distance from the cylinder's end faces. So that for small values of b 1 , the cylinder can be regarded as an infinitely long. On the contrary, the greater values of b 1 allow for capturing the end effects in the distribution of residual stresses due to the strain (26) . Figures 2 -5 present the plots for the radial, axial, circumferential, and tangential stresses in a "cuboid" cylinder C (i.e., the one with ratio equal to one, b = 1) due to the residual strain (26) computed under the following parameters: e 1 = 1, e 2 = −1, r 1 = 1/3, r 2 = 2/3, b = 1, a r = 18, a z = 10. The Poisson ratio in all considered case studies is ν = 0.3. In the figures, curves 1 -5 correspond to the values r = 0; 1/6; 1/2; 5/6; 1, respectively, indicating the axis, cylindrical midsurface and lateral surface of the cylinder, as well as the middles of the disturbed zones. The solutions of systems in (25) were computed through the simple reduction algorithm by holding N = 15 terms in the sums by k. As we can observe in Fig. 2 , the computed radial stress exactly satisfies the boundary condition (2) on the lateral surface (curves 5) under all the considered distribution profiles of the residual strain (26) . For b 1 = 0.25, and smaller, the stress exhibits the trend typical for the distributions in infinitely long cylinder, i.e., tending to zero at cross-sections z = const, which are far enough from the plane z = 0. The magnitude of the residual radial stress at the cylinder axis decreases with decrement of b 1 . For the greater values of b 1 , curves 3 and 4 are remarkably close to one another. However, for the case studies of lesser end-effect, i.e., for b 1 = 0.25, the difference between the radial stress at these radii is more pronounced at z = 0. Quite similar behavior is observed in the circumferential Fig. 2 The radial stress versus axial coordinate z at the specific radii r = 0 (curve 1), r = 1/6 (curve 2), r = 1/2 (curve 3), r = 5/6 (curve 4), and r = 1 (curve 5) of a "cubic" cylinder C under the residual strains given by (26) , where e 1 = 1, e 2 = −1, r 1 = 1/3, r 2 = 2/3, b = 1, a r = 18, a z = 10, and b 1 = 1.00; 0.75; 0.50; 0.25 Fig. 3 The axial stress versus axial coordinate z at the specific radii r = 0 (curve 1), r = 1/6 (curve 2), r = 1/2 (curve 3), r = 5/6 (curve 4), and r = 1 (curve 5) of a "cubic" cylinder C under the residual strains given by (26) , where e 1 = 1, e 2 = −1, r 1 = 1/3, r 2 = 2/3, b = 1, a r = 18, a z = 10, and b 1 = 1.00; 0.75; 0.50; 0.25 Fig. 4 The circumferential stress versus axial coordinate z at the specific radii r = 0 (curve 1), r = 1/6 (curve 2), r = 1/2 (curve 3), r = 5/6 (curve 4), and r = 1 (curve 5) of a "cubic" cylinder C under the residual strains given by (26) , where e 1 = 1, e 2 = −1, r 1 = 1/3, r 2 = 2/3, b = 1, a r = 18, a z = 10, and b 1 = 1.00; 0.75; 0.50; 0.25 stress depicted in Fig. 4 . It is worth noting that the computed radial stress fits the necessary condition of the self-equilibration for the residual stress, i.e., The axial stress in Fig. 3 exactly satisfies the boundary condition (2) on the end faces of the cylinder for all values of the radial coordinate. Similarly to the radial stress, the disturbance due to the residual strain (26) occupies almost entire length of the cylinder for greater values of b 1 , i.e., b 1 = 1; 0.75. This effect is well pronounced, in particular, on the lateral surface of the cylinder (curves 5), while the magnitude of this stress on the axis of the cylinder (curves 1) remains almost the same. The tangential stress (Fig. 5) is an odd function of the longitudinal coordinate z satisfying the zero boundary conditions on both the lateral surface and end faces. For greater values of b 1 , this stress exhibits smaller maximum deviation from the zero distributions, which grows with narrowing the zone of the distributions of the residual strain (26) . Let us analyze the efficiency of the proposed computational algorithm with regard to solving the systems (25) by using the simple reduction method. The solution errors for the first and second disturbed problems depending on number of equations in the reduced systems (25) can be estimated with evaluating the corresponding functionalities calculated on the solution obtained for a selected value of N : Here, F N 1 and F N 2 are the values of the functionals (22) and (23) The tangential stress versus axial coordinate z at the specific radii r = 0 (curve 1), r = 1/6 (curve 2), r = 1/2 (curve 3), r = 5/6 (curve 4), and r = 1 (curve 5) of a "cubic" cylinder C under the residual strains given by (26) , where e 1 = 1, e 2 = −1, r 1 = 1/3, r 2 = 2/3, b = 1, a r = 18, a z = 10, and b 1 = 1.00; 0.75; 0.50; 0.25 Table 1 The errors in the minimization functionals (22) and (23) Table 1 . As we can see, the accuracy gradually grows for greater values of N . The convergence of the solutions depends on the type of boundary conditions. On the basis of conducted numerical experiments, we can conclude that accuracy sufficient for practical goals can be achieved at N 9. A technique for evaluating axisymmetric residual stresses in a finite-length solid cylinder is presented. The cylinder's surface is assumed to be free of force loadings. The residual stresses are induced by incompatible strains represented by an isotropic tensor, whose impact is encountered via the model of conventional-plastic strains. An analytical solution to this problem is constructed by superposing the solutions of three individual problems. The first problem implies determining the basic state for an infinitely long solid cylinder of the same radius, which undergoes the considered distribution profiles of the incompatible strains. In order to find an exact solution to this problem, an analytical technique was employed that is based upon the implementation of the Goodier-Papkovich potential constructed through the complete Legendre elliptic integral of the first kind. In order to compensate the stresses induced by the found stress field when "cutting" the original cylinder of finite length out of the infinitely long one, two auxiliary problems are considered. The solutions to these problems are constructed by implementing the Love biharmonic stress function and the general scheme of the method of homogeneous solutions, i.e., the ones satisfying the zero boundary conditions on the segments of the cylinder boundary, represented by different families of coordinate surfaces. These two problems are solved with the use of the variational method, in which the subordination of the solution to the boundary conditions is performed in the norm L 2 . As a result, the problems are reduced to infinite systems of linear algebraic equations regarding the unknown complex coefficients. Because of the zero asymptotic of the coefficients, the obtained systems can be solved within any given accuracy by implementing the simple reduction method. This makes the obtained system to be advantageous over the ones obtained through the method of cross-wise superposition [33] . The conducted numerical experiments for the considered case studies confirm high convergence of the suggested algorithm, as sufficient accuracy is reached at N 9. The constructed solution allows for analyzing the residual stresses in the finite-length cylinders with different ratios and various nature of the incompatible strains. The solution can be efficiently extended for the analysis of thermal stresses by implying the incompatible strains to be the thermal ones, i.e., ε 0 ll = α(T (r, z) − T 0 (r, z)), where l = {r, θ, z}, α is the linear thermal expansion coefficient, T 0 (r, z) is the reference temperature of the stress-free cylinder, and T (r, z) is the actual temperature. Finally, in view of the fact that the solution is constructed in an analytical form, it can be rather advantageous for further implementation with the development of non-destructive evaluation techniques for the residual stresses in elastic cylinders of finite length. In: Schajer GS (ed) Practical residual stress measurement methods Welding residual stresses and effects on fracture in pressure vessel and piping components: a millennium review and beyond Understanding residual stress and distortion in welds: an overview. In: Feng Z (ed) Processes and mechanisms of welding residual stress and distortion Lifetime estimation of welded joints In: Schajer GS (ed) Practical residual stress measurement methods Residual stress measurement combining blind-hole drilling and digital image correlation approach Deep hole drilling. In: Schajer GS (ed) Practical residual stress measurement methods Effect of employing buffer layer in repaired dissimilar welded joints on the residual stresses based on contour and slitting methods The slitting method. In: Schajer GS (ed) Practical residual stress measurement methods The contour method. In: Schajer GS (ed) Practical residual stress measurement methods Applied and residual stress determination using X-ray diffraction. In: Schajer GS (ed) Practical residual stress measurement methods Concepts for nondestructive and depth-resolved X-ray residual stress analysis in the near-surface region of nearly single crystalline materials with mosaic structure Neutron diffraction. In: Schajer GS (ed) Practical residual stress measurement methods Magnetic methods. In: Schajer GS (ed) Practical residual stress measurement methods Ultrasonics. In: Schajer GS (ed) Practical residual stress measurement methods Optical methods. In: Schajer GS (ed) Practical residual stress measurement methods Methods for determining residual welding stresses and their relief (a review) Welding residual stresses and effects on fracture in pressure vessel and piping components: a millennium review and beyond A nondestructive numerical-experimental method of determining the residual stresses in welded shells of revolution The inverse conditionally correct problem of determining the residual stresses in compound welded shells of revoltuion Thermoelastic distorsion problems Determination of the stress state of thin shells taking into account strains due to physicochemical phenomena Micromechanics of defects in solids Distribution of residual stresses induced by axially symmetric plastic strains in a layer Distribution of residual stresses in a layer containing a welded cylindrical disk Analysis of residual stresses in a long hollow cylinder Determination of the stressed state of main pipelines in the zones of circular welds Analytical determination of residual stresses in a butt-weld of two thin rectangular plates Investigation of the plane stressed state in a rectangular domain Selected topics in the history of the two-dimensional biharmonic problem The direct integration method for elastic analysis of nonhomogeneous solids Axisymmetric thermal stresses in an elastic hollow cylinder of finite length Axially symmetric temperature stresses in an elastic isotropic cylinder of finite length Equilibrium of an elastic finite cylinder under axisymmetric discontinuous normal loadings The biharmonic problem in the theory of elasticity Theoretical and experimental determination of residual stresses in plane joints Problem of nondestructive evaluation of residual stresses in a pipeline according to the data of magnetoelastic measurements Residual stresses in a finite cylinder. Direct and inverse problems and their solving using the variational method of homogeneous solutions Application of the variational method of homogeneous solutions for the determination of axisymmetric residual stresses in a finite cylinder Theory of elasticity and thermal stresses. Explanations, problems and solutions Elasticity. Theory, applications, and numerics General integral of thermal stresses (on the work 'Thermal stresses in the theory of elasticity On the integration of the thermo-elastic equations A dual-reciprocity boundary element approach for solving axisymmetric heat equation subject to specification of energy A treatise on the mathematical theory of elasticity Application of the least squares method in axisymmetric biharmonic problems Theory of elasticity New York Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations Acknowledgements This paper is devoted to the brilliant memory of Prof. Vasyl Chekurin (1951-2021) whose life was taken by COVID-19 in the middle of his research into this subject. Functions in expressions (19) for the stresses of the first disturbed problem: σ (1) rr,k (r, z) = 2μγ 2 k cosh(γ k z) κ k ((1 − 2ν)J 0 (γ k r ) − r γ k J 1 (γ k r )) + Functions in expressions (21) for the stresses of the second disturbed problem:θθ,k (r, z) = 2μrr, (z) = σ (2) rr, (1, z) , σ 2 rr, (z) =σ (2) rr, (1, z), σ 1 r z, (z) = σ (2) r z, (1, z) , σ 2 r z, (z) =σ