Analysis of residual stresses in bars and tubes of cylindrically orthotropic materials

One of the most frequently used techniques of determining residual stresses in cylindrical bodies of elastically isotropic materials is the Mesnager-Sachs boring-out method. This method is adequate for many engineering materials; however, certain materials possess elastic properties which vary with direction. Materials having one fairly common kind of directional variation are known as cylindrically orthotropic materials. Heretofore, no means has been available for determining the residual stresses from strain data taken after successive borings of members made of such materials. This paper gives the derivation of equations for this purpose along with graphs depicting the error to be expected if the data had been reduced by the Sachs equations.     

Measurement of Residual Stresses in Nuclear-grade Zircaloy-4(R) Tubes—Effect of Heat Treatment

Nuclear-grade Zircaloy-4(R) tubes are produced by a unique manufacturing process known as pilgering, which leaves the material in a work-hardened state containing a pattern of residual stresses. Moreover, such tubes exhibit elastic anisotropy as a result of the pilgering process. Therefore, standard equations originally proposed by Sachs (Z Met Kd, 19: 352–357, 1927; Sachs, Espey, Iron Age, 148: 63–71, 1941). for isotropic materials do not apply in this situation. Voyiadjis et al. (Exp Mech, 25: 145–147, 1985) proposed a set of equations for treating elastically anisotropic materials, but we have determined that there are discrepancies in their equations. In this paper, we present the derivation for a set of new equations for treating elastically anisotropic materials, and the application of these equations to residual stress measurements in Zr-4(R) tubes. To this end, through thickness distribution of residual stress components in as-received and heat treated (500°C) Zr-4(R) tubes was measured employing the Sachs’ boring-out technique in conjunction with electrochemical machining as the means of material removal, and our new equations. For both as-received and the heat treated materials, the axial and tangential residual stresses were significantly higher than the radial and shear residual stresses. The largest residual stress was the tangential stress component in the as-received material, showing a tensile value at the outer surface and a compressive value at the inner surface. At high values of von Mises equivalent stress, the principal directions of residual stress coincided with the principal axes of the tube for the as-received material, as well as for the material heat treated at 500°C.     

Shear of rubber tube springs

Rubber tube springs consist basically of cylindrical rubber tubes bonded on their inner and outer curved surfaces to rigid cylindrical tubes. They are widely used as flexible linkages, for example in vehicle suspensions. Rotation of one rigid tube with respect to the other about their common axis subjects the rubber tube to azimuthal shear. Displacement of one rigid tube with respect to the other along their common axis puts the rubber tube into axial shear. Using FEA, we have calculated the stresses set up in both cases, for a long rubber tube of a non-linearly elastic (neo-Hookean) material. The results are compared for the two modes of deformation, and with analytical predictions where available. For a long tube the shear stresses are substantially independent of the end conditions, but the normal stresses are strongly affected, as found previously for sheared rectangular blocks [A.N. Gent, J.B. Suh, S.G. Kelly III, Mechanics of rubber shear springs, Int. J. Nonlinear Mech. 42 (2007) 241-249]. If the end surfaces are stress-free, unexpectedly large normal stresses are generated, even in azimuthal shear. These high tensile stresses are attributed to restraints at the inner and outer cylindrical boundaries that compensate for the absence of stresses on the end surfaces that would be needed to maintain a simple shear deformation. Thus, the boundary conditions affect the stresses everywhere (in contrast to an “end effect” that would diminish away from the ends). Small departures from complete incompressibility are found to lower the internal stresses markedly, and even cause the sign of the stresses to be reversed.     

Analysis of residual stresses in cylindrically anisotropic materials

Metal-forming operations leave residual stresses in formed parts due to nonuniform deformation occurring during the process. An exact method of determining the longitudinal, radial and circumferential (tangential) residual stresses in axisymmetric specimens was proposed by Mesnager¹ and further developed by Sachs². The boring-out technique can be complemented by a similar procedure in which strains are measured on the inner surface of the tube when material is removed from the outer surface.The work proposed in this paper extends previous analyses of residual stresses to the case where the material exhibits cylindrical elastic anisotropy, i.e., the principal axes of anisotropy correspond to the longitudinal, radial and circum-ferential directions of the tube. In addition, the present analysis considers the case in which a residual-shear stress, developed by twisting the tube about its axis, exists in the tube. When such shearing stresses are present, the principal axes of the residual-stress distribution are not parallel to the principal axes of the tube.     

Intensity of singular stress fields causing interfacial debonding at the end of a fiber under pullout force and transverse tension

In this study, singular stress fields at the ends of fibers are discussed by the use of models of rectangular and cylindrical inclusions in a semi-infinite body under pullout force. Those singular stresses have not been discussed yet in the previous studies for pullout problems although they are important for causing interfacial initial debonding. The body force method is used to formulate those problems as a system of singular integral equations where unknowns are densities of the body forces distributed in a semi-infinite body having the same elastic constants as those of the matrix and inclusions. In order to compare the results with the previous solutions, tension problems of a fiber in a semi-infinite body are also considered. Then, generalized stress intensity factors at the corner of rectangular and cylindrical inclusions are systematically calculated for various geometrical conditions with varying the elastic ratio, length, and spacing of the location from edge to inner of the body. The effects of elastic modulus ratio and aspect ratio of inclusion upon the stress intensity factors are discussed for pullout problems.     

Analysis of laminated piezoelectric cylindrical shells

A new method is developed for three-dimensional stress analysis of laminated piezoelectric cylindrical shell with simple support. The shell can be subjected to various applied loadings, including distributed body force, inner and outer surface traction and potential. Each layer of the shell can be piezoelectric or elastic/dielectric, with perfect bonding assumed between each interface. The governing equations are solved by the state-space technique. Numerical results are presented to show the sensing and actuating effects of three-layered piezoelectric cylindrical shell. The project supported by the National Natural Science Foundation of China (19572027)     

A mixed problem of plane elasticity for a domain with partially unknown boundary

The paper addresses a problem of plane elasticity for a doubly connected body with outer and inner boundaries in the form of regular polygons with a common center and parallel sides. The neighborhoods of the vertices of the inner boundary are unknown equal full-strength smooth arcs symmetric about the rays coming from the vertices to the center. It is assumed that this elastic body is inserted into a hole of a rigid body, with the hole boundary coinciding with the outer boundary of the elastic body. Absolutely smooth rigid punches with rectilinear bases are pressed into all the rectilinear sections of the inner polygonal boundary of the elastic body. There is no friction between the elastic and rigid bodies. The unknown arcs are free from external stresses. Complex variable theory is used to determine the unknown arcs and the stress state of the elastic body __________ Translated from Prikladnaya Mekhanika, Vol. 42, No. 3, pp. 110–118, March 2006.     

Time-dependent thermoelastic creep analysis of rotating disk made of Al–SiC composite

Time-dependent creep stress redistribution analysis of rotating disk made of Al–SiC composite is investigated using Mendelson’s method of successive elastic solution. All mechanical and thermal properties except Poisson’s ratio are radial dependent based on volume fraction percent of SiC reinforcement. The material creep behavior is described by Sherby’s constitutive model using Pandey’s experimental results on Al–SiC composite. Loading is an inertia body force due to rotation and a distributed temperature field due to steady-state heat conduction from inner to outer surface of the disk. Using equations of equilibrium, stress strain, and strain displacement, a differential equation, containing creep strains, for displacement is obtained. History of stresses and deformations are calculated using method of successive elastic solution. It is concluded that the uniform distribution of SiC reinforcement does not considerably influence on stresses. However, the minimum and most uniform distribution of circumferential and effective thermoelastic stresses belongs to composite disk of aluminum with 0% SiC at inner surface and 40% SiC at outer surface. It has also been found that the stresses, displacement, and creep strains are changing with time at a decreasing rate so that after almost 50 years the solution approaches the steady-state condition.     

夹层FGM圆柱壳在扭转载荷作用下的弹性稳定性

采用半解析方法研究了两端简支的功能梯度夹层圆柱壳在端部扭转载荷作用下的弹性稳定性.考虑圆柱壳的里外表层为均匀材料,中间层为材料性质沿厚度方向连续变化的功能梯度材料,并且在界面处的材料性质保持连续. 基于Flügge薄壳理论,建立了位移形式的结构静态屈曲控制方程.根据边界条件将位移表示为三角级数形式,获得包含柱壳端部扭转载荷参数的近似线性代数特征值问题,并通过数值方法求得了表征结构失稳特征的临界载荷. 数值结果表明,临界载荷随着半径与厚度比的增加而减小,随着功能梯度中间层的弹性模量的平均值的增加而增加.      

Nonlinear dynamic responses of sandwich functionally graded porous cylindrical shells embedded in elastic media under 1:1 internal resonance

In this article, the nonlinear dynamic responses of sandwich functionally graded(FG) porous cylindrical shell embedded in elastic media are investigated. The shell studied here consists of three layers, of which the outer and inner skins are made of solid metal, while the core is FG porous metal foam. Partial differential equations are derived by utilizing the improved Donnell's nonlinear shell theory and Hamilton's principle. Afterwards, the Galerkin method is used to transform the governing equations into nonlinear ordinary differential equations, and an approximate analytical solution is obtained by using the multiple scales method. The effects of various system parameters,specifically, the radial load, core thickness, foam type, foam coefficient, structure damping,and Winkler-Pasternak foundation parameters on nonlinear internal resonance of the sandwich FG porous thin shells are evaluated.     

应用增量理论考虑应变硬化效应包辛格效应和温度对机械性质影响的自增强厚壁圆筒分析

本文提出的采用增量理论的分析模型和方法可以定量分析应变硬化、包辛格效应和材料机械性质随温度变化等因素对自增强厚壁圆筒中残余应力和操作应力分布的影响.厚壁圆筒看成由N个同轴薄壁圆筒套在一起组成的构件;采用包括弹性、塑性和温度应变的增量型本构关系和相容条件导出了自增强圆筒的基本方程;设计了计算机程序并给出了分析实例。分析结果表明,应变硬化会减小塑性区并降低残余应力;包辛格效应使反向屈服容易出现;温度升高将使残余应力和热应力松驰.     

Stress intensity factors around a cylindrical crack in an interfacial zone in composite materials

Axisymmetric stresses around a cylindrical crack in an interfacial cylindrical layer between an infinite elastic medium with a cylindrical cavity and a circular elastic cylinder made of another material have been determined. The material constants of the layer vary continuously from those of the infinite medium to those of the cylinder. Tension surrounding the cylinder and perpendicular to the axis of the cylinder is applied to the composite materials. To solve this problem, the interfacial layer is divided into several layers with different material properties. The boundary conditions are reduced to dual integral equations. The differences in the crack faces are expanded in a series so as to satisfy the conditions outside the crack. The unknown coefficients in the series are solved using the conditions inside the crack. Numerical calculations are performed for several thicknesses of the interfacial layer. Using these numerical results, the stress intensity factors are evaluated for infinitesimal thickness of the layer.     

On Foundations of the Ultrasonic Nondestructive Method for Determination of Stresses in Near-Surface Layers of Solid Bodies

The ultrasonic nondestructive method for determining stresses in near-surface layers of solid bodies is based on the laws governing the propagation of elastic surface waves in bodies with initial (residual) stresses. These laws are established within the framework of the linearized three-dimensional theory of waves in bodies with initial (residual) stresses. The dispersion equations in associated problems are solved by computational methods. The nondestructive method and measuring instruments and devices are described. Some examples of nondestructive determination of welding-induced residual stresses and operating stresses in near-surface layers of materials are presented Presented at the International Conference on Computational and Experimental Engineering and Sciences (ICCES'04) (Madeira, Portugal, July 26–29, 2004) and published in the journal Computer Modeling in Engineering and Sciences (CMES). __________ Published in Prikladnaya Mekhanika, Vol. 41, No. 8, pp. 130–144, August 2005.     

Development of a computer program for the prediction of flow in a rotating cavity

A computer program has been developed to predict laminar source-sink flow in a rotating cylindrical cavity. Although the program is based on a standard finite difference technique for recirculating flow, it incorporates two novel features. Step changes in grid size are employed to obtain sufficient resolution in the boundary layers and special treatment is given to the solution of the pressure correction equations, in the ‘SIMPLE’ algorithm, in order to improve the convergence properties of the method. Results are presented both for the flow in an infinite rotating cylindrical annulus and a finite rotating cylindrical cavity, with the inner cylindrical surface acting as a uniform source and the outer cylinder as a sink. These show good agreement with existing analytical solutions and illustrate some of the problems associated with the computation of rapidly rotating flows.     

Dispersion of radial vibrations of an orthotropic cylindrical tube placed in a magnetic field

The present paper studies the dispersion relation of the radial vibrations of an orthotropic cylindrical tube. The effects of the magnetoelastic interaction on the problem are investigated. The problem is represented by the equations of elasticity taking into account the effect of the magnetic field as given by Maxwell's equations in the quasi-static approximation. The stress free conditions on the inner and outer surfaces of the hollow cylindrical cube are satisfied to form a dispersion relation in terms of the wavelength, the cylinder radii and the material constants. This study shows that waves in a solid body propagating under the influence of a superimposed magnetic field can differ significantly from those propagating in the absence of a magnetic field. The results have been verified numerically and represented graphically.     

Continuum Dynamics on Manifolds: Application to Elasticity of Residually-Stressed Bodies

This paper is concerned with the dynamics of continua on differentiable manifolds. We present a covariant derivation of the equations of motion, viewing motion as a curve in an infinite-dimensional Banach manifold of embeddings of a body manifold in a space manifold. Our main application is the motion of residually-stressed elastic bodies, where the residual stresses result from a geometric incompatibility between body and space manifolds. We then study a particular example of elastic vibrations of a two-dimensional curved annulus embedded in a sphere.     

WAVE PROPAGATION IN A SYSTEM OF COAXIAL TUBES FILLED WITH INCOMPRESSIBLE MEDIA: A MODEL OF PULSE TRANSMISSION IN THE INTRACRANIAL ARTERIES

The propagation of harmonic waves through a system formed of coaxial tubes filled with incompressible continua is considered as a model of arterial pulse propagation in the craniospinal cavity. The inner tube represents a blood vessel and is modelled as a thin-walled membrane shell. The outer tube is assumed to be rigid to account for the constraint imposed on the vessels by the skull and the vertebrae. We consider two models: in the first model the annulus between the tubes is filled with fluid; in the second model the annulus is filled with a viscoelastic solid separated from the tubes by thin layers of fluid. In both models, the elastic tube is filled with fluid. The motion of the fluid is described by the linearized form of the Navier–Stokes equations, and the motion of the solid by classical elasticity theory. The results show that the wave speed in the system is lower than that for a fluid-filled elastic tube free of any constraint. This is due to the stresses generated to satisfy the condition that the volume in the system has to be conserved. However, the effect of the constraint weakens as the radius of the outer tube is increased, and it should be insignificant for the typical physiological parameter range.     

Plasticity Effects in the Hole-Drilling Residual Stress Measurement in Peened Surfaces

The incremental hole-drilling technique (IHD) is a widely established and accepted technique to determine residual stresses in peened surfaces. However, high residual stresses can lead to local yielding, due to the stress concentration around the drilled hole, affecting the standard residual stress evaluation, which is based on linear elastic equations. This so-called plasticity effect can be quantified by means of a plasticity factor, which measures the residual stress magnitude with respect to the approximate onset of plasticity. The observed resultant overestimation of IHD residual stresses depends on various factors, such as the residual stress state, the stress gradients and the material’s strain hardening. In peened surfaces, equibiaxial stresses are often found. For this case, the combined effect of the local yielding and stress gradients is numerically and experimentally analyzed in detail in this work. In addition, a new plasticity factor is proposed for the evaluation of the onset of yielding around drilled holes in peened surface layers. This new factor is able to explain the agreement and disagreement found between the IHD residual stresses and those determined by X-ray diffraction in shot-peened steel surfaces.     

Bounds for the torsional rigidity of inhomogeneous cylindrical bars

The object of this paper is the uniform torsion of inhomogeneous, isotropic, linearly elastic cylindrical bar. The aim is to give lower and upper bounds for the torsional rigidity of the bar with doubly connected cross section. The outer and inner boundary curves of cross section are similar curves. The level lines of the function which gives the change of the shear rigidity on the cross section are also similar curves to the boundary curves. The application of derived bounding formulae is illustrated by examples. An approximated formula to determine the shear stresses is also presented.     

An Exact Derivation from Three-Dimensions of the Linear Equations for the Bending of an Elastically Monoclinic Plate and Associated Three-Dimensional Solutions

The exact linear three-dimensional equations for a elastically monoclinic (13 constant) plate of constant thickness are reduced without approximation to a single 4th order differential equation for a thickness-weighted normal displacement plus two auxiliary equations for weighted thickness integrals of a stress function and the normal strain. The 4th order equation is of the same form as in classical (Kirchhoff) theory except the unknown is not the midsurface normal displacement. Assuming a solution of these plate equations, we construct so-called modified Saint-Venant solutions—“modified” because they involve non-zero body and surface loads. That is, solutions of the exact three-dimensional elasticity equations that exhibit no boundary layers and that are subject to a special set of body and surface loads that leave the analogous plate loads arbitrary.