Numerical Solution of a Problem of Thermal Stresses of a Magnetothermoelastic Cylinder with Rotation by Finite-Difference Method

The present article deals with the investigation thermal stress of a magnetothermoelastic cylinder subjected to rotation, open or closed circuit, thermal and mechanical boundary conditions. The outer and inner surfaces of the cylinder are subjected to both mechanical and thermal boundary conditions. A The transient coupled thermoelasticity in an infinite cylinder with its base abruptly exposed to a heat flux of a decaying exponential function of time is devised solve by the finite-difference method. The fundamental equations’ system is solved by utilizing an implicit finite-difference method. This current method is a second-order accurate in time and space; it is also unconditionally stable. To illustrate the present model’s efficiency, we consider a suitable material and acquire the numerical solution of temperature, displacement components, and the components of stresses with time t and through the radial of an infinite cylinder. The results indicate that the effect of coupled thermoelasticity, magnetic field, and rotation on the temperature, stresses, and displacement is quite pronounced. In order to illustrate and verify the analytical developments, the numerical solution of partial differential equations, stress components, displacement components and temperature is carried out and computer simulated results are presented graphically. This study is helpful in the development of piezoelectric devices.     

Transient thermal stress intensity factors for an internal longitudinal semi-elliptical crack in a thick-walled cylinder

In this paper, the stress intensity factors are derived for an internal semi-elliptical crack in a thick-walled cylinder subjected to transient thermal stresses. First, the problem of transient thermal stresses in a thick-walled cylinder is solved analytically. Thermal and mechanical boundary conditions are assumed to act on the inner and outer surfaces of the cylinder. The quasi-static solution of the thermoelasticity problem is derived analytically using the finite Hankel transform and then, the stress intensity factors are extracted for the deepest point and the surface points of the semi-elliptical crack using the weight function method. The results show to be in accordance with those cited in the literature in the special case of steady-state problem. Using the closed-form relations extracted for the transient thermal stress intensity factors, some conclusive results are drawn.     

Inverse problem of material distribution for desired fracture characteristics in a thick-walled functionally graded material cylinder with two diametrically-opposed edge cracks

This study concerns the inverse problem of evaluating the optimum material distribution for desired fracture characteristics in a thick-walled functionally graded material (FGM) cylinder containing two diametrically-opposed edge cracks emanating from the inner surface of the cylinder. The thermal eigenstrain developed in the cylinder material due to nonuniform coefficient of thermal expansion as a result of cooling from sintering temperature is taken into account. Based on a generalized method of evaluating stress intensity factors developed in a previous study, an inverse method is developed to optimize material distribution intending to realize prescribed apparent fracture toughness in the FGM cylinder. To present some numerical results, a TiC/Al₂O₃ FGM cylinder is considered and the inverse problems are solved to evaluate material distributions for two examples of prescribed apparent fracture toughness. The effect of cylinder wall thickness on the material distribution and comparison of material distributions corresponding to a single and two cracks are also discussed. The numerical results reveal that the apparent fracture toughness of FGM cylinders can be controlled by choosing the material distributions properly.     

On the singularity of temperature gradient near an inclined crack terminating at bimaterial interface

This paper deals with the singularity of temperature gradient near an inclined crack terminating at a bimaterial interface. The temperature field is solved by considering the continuity of temperature and heat flux at the interface and appropriate thermal boundary conditions on crack surfaces. The singularity of temperature gradient around the crack tip is then studied for the cases for which the temperature on crack surfaces is prescribed or crack surfaces are insulated. It is found that, unlike the oscillatory singularity of the stress field, no oscillatory character near the crack tip is observed for these problems. The dependence of the singularity of temperature gradient on the inclined angle of crack and thermal conductivity ratio of two dissimilar media is also shown.     

Transient thermal contact problem for axial crack in a hollow circular cylinder

The transient behavior of an axial-cracked hollow circular cylinder subjected to a sudden heating is investigated. It is shown that surface heating may induce compressive thermal stress near the inner surface of the cylinder which in turn may force the cracked surfaces to close together. Assuming that the existence of the crack does not alter the temperature distribution, this problem can be divided into two parts and solved by the principle of superposition. First, the temperature and transient thermal stress distributions along the axisymmetric surface of the imaginary cylinder without a crack are obtained by finite element/implicit time integration method. The calculated temperature and thermal stress distributions are in good agreement with the values predicted by the analytical method. Secondly, the opposite senses of the stress distributions along the cracked surfaces, which are obtained previously, are treated as the traction boundary conditions, and the contact length and contact pressure of the real cracked cylinder are obtained by a modified elimination finite element scheme. In this scheme, the concepts of contact-node-pairs' penetration, contact-double-forces and compliance matrix are introduced. The calculated results indicate that the contact length ratio becomes smaller when the crack length ratio increases, and becomes larger as the radius ratio increases. Finally, the normalized stress intensity factor for the crack tip of the cylinder is obtained. It is shown that the larger the crack length ratio the higher the stress intensity factor.     

Interfacial Effects on Stress Transfer in Fiber-Reinforced Composites

A three-dimensional solution for a fiber-reinforced circular cylinder is presented for axisymmetric force and displacement boundary conditions. The solution, which can satisfy all the boundary conditions prescribed on the curved and end surfaces of the cylinder, can be used directly in the micromechanical analysis of fiber-reinforced composites to investigate the typical representative volume element (RVE). The element consists of a combined circular cylinder composed of a solid inner circular cylinder of transversely isotropic fiber, a concentric outer circular cylinder of isotropic matrix material, and an interface layer between the fiber and the matrix. The radial and tangential flexural compliances of the interfacial material are considered, and their effects on the stress transfer at the interface are studied parametrically. The numerical results presented show that the material properties of the interfacial layer have significant influences on the stress distribution within the RVE, particularly at the cross sections near the ends of the cylinder, where external loads are applied.     

Wave propagation modeling in human long bones

Summary The dynamic behavior of a dry long bone that has been modeled as a piezoelectric hollow cylinder of crystal class 6 is investigated. The solution for the wave propagation problem is expressed in terms of a potential function which satisfies an eighth-order partial differential equation, whose solutions lead to the derivation of the explicit solution of the wave equation. The mechanical boundary conditions correspond to those of stress free lateral surfaces, while the electrical boundary conditions correspond to those of short circuit. The satisfaction of the boundary conditions leads to the dispersion relation which is solved numerically. The frequencies obtained are presented as functions of various parameters and they compare well with other researchers' theoretical results.     

Meshless local Petrov–Galerkin method for coupled thermoelasticity analysis of a functionally graded thick hollow cylinder

In this article, coupled thermoelasticity (without energy dissipation) based on Green–Naghdi model is applied to functionally graded (FG) thick hollow cylinder. The meshless local Petrov–Galerkin method is developed to solve the boundary value problem. The Newmark finite difference method is used to treat the time dependence of the variables for transient problems. The FG cylinder is considered to be under axisymmetric and plane strain conditions and bounding surfaces of cylinder to be under thermal shock loading. The mechanical properties of FG cylinder are assumed to vary across thickness of cylinder in terms of volume fraction as nonlinear function. A weak formulation for the set of governing equations is transformed into local integral equations on local subdomains by using a Heaviside test function. Nodal points are regularly distributed along the radius of the cylinder and each node is surrounded by a uni-directional subdomain to which a local integral equation is applied. The Green–Naghdi coupled thermoelasticity equations are valid in each isotropic subdomain. The temperature and radial displacement distributions are obtained for some grading patterns of FGM at various time instants. The propagation of thermal and elastic waves is discussed in details. The presented method shows high capability and efficiency for coupled thermoelasticity problems.     

Weight function for an external circumferential semielliptical crack in a cylinder

In this paper, the weight function was extracted at the deepest point of a semielliptical circumferential crack. The crack is assumed to exist on the outer surface of the cylinder. For this purpose, the three‐dimensional finite element method was accomplished to specify two reference loads, which are indispensable for determining the weight function. The verification study confirms the accuracy of the derived weight function under prescribed mechanical loading on the crack surfaces. There is consistency among the solution results compared with those in the literature. The second part describes the application of the weight function for the thermal boundary conditions. Steady‐state thermal stress intensity factors are demonstrated using the weight function and presented as a closed‐form solution. The results were compared with the finite element data on the special case of thermal loading, and good agreement is obtained.     

Static deformations of functionally graded polar-orthotropic cylinders with elliptical inner and circular outer surfaces

Functionally graded materials (FGMs) enable one to tailor the spatial variation of material properties so as to fully use the material everywhere. For example, in a hollow circular cylinder one can vary, in the radial direction, the material moduli to make the hoop stress constant. Whereas the problem for a hollow cylinder with the inner and the outer surfaces circular has been studied, that of a cylinder with a circular outer surface and a non-circular inner surface or vice versa has not been investigated. We study here such a plane-strain problem when the cylinder material is polar-orthotropic, material properties vary exponentially in the radial direction, and deformations are independent of the axial coordinate. The problem is challenging since the cylinder thickness varies with the angular position of a point, and the cylinder material is inhomogeneous. Equilibrium equations are solved by expanding the radial and the circumferential displacements in Fourier series in the angular coordinate. The method of Frobenius series is used to solve ordinary differential equations for coefficients of the Fourier series, and boundary conditions are satisfied in the sense of Fourier series. A parametric study has been conducted that delineates effects on stresses of the eccentricity of the ellipse, the material property gradation index and loads applied on boundaries of the cylinder. The analytical solutions presented here will serve as benchmarks for comparing solutions derived by numerical methods.     

A three-dimensional treatment of transient thermal stresses in a transversely isotropic semi-infinite circular cylinder subjected to an asymmetric temperature on the cylindrical surface

Summary This paper deals with a three-dimensional transient thermal stress problem in a transversely isotropic semi-infinite circular cylinder subjected to an asymmetrical temperature on the cylindrical surface and a convective heat loss on the surfaces. In analyzing the problem, the generalized Fourier transform to a temperature and the modified transversely isotropic potential functions method to stresses are used. For the numerical example, the temperature and the thermal stresses are calculated for a grafite which belongs under a transverse isotropy, and compared with the values under isotropic conditions. The effects of the various anisotropies of the material properties on the temperature and the thermal stresses are studied.With 15 Figures     

Shear induced loss of saturation in a fluid infused swollen hyperelastic cylinder

We discuss a continuum model for the absorption and redistribution of fluid in swollen elastomers due to mechanical loading; and in particular distinguish between systems that are saturated with liquid and those that are not. To this end we consider a boundary value problem of radial displacement combined with azimuthal shear for an annular cylinder consisting of a fluid infused hyperelastic media. In the absence of load the elastomer deforms by free swelling, giving a homogeneous expansion in which the imbibed fluid is uniformly distributed. This free swelling is described by the theory of Flory and Rehner. We then consider the effect of various boundary displacements and tractions so as to study how this alters the uniform fluid distribution. This problem was previously considered by Wineman and Rajagopal for the case in which the lateral surfaces maintained the radius associated with free swelling. Here we consider certain generalizations in which the lateral surfaces may undergo not only a prescribed relative twist, but also radial displacement. This permits fluid to either enter or exit the cylinder. A numerical method is invoked to solve these boundary value problems using representative material parameters. Certain boundary tractions generate an overall volume increase after the free swelling. If the amount of available fluid is limited, this gives rise to the possibility of complete fluid absorption, whereupon the system is no longer saturated. It is found that the overall mechanical response after loss of saturation is stiffer than it would be if the system had remained saturated.     

Coupled transient thermoelastic contact problems for axial cracks in hollow cylinders

A coupled transient thermoelastic behaviour of an axial-cracked hollow circular cylinder subjected to a sudden heating is investigated in this study. It is shown that surface heating may induce the compressive thermal stress near the inner surface of the cylinder which in turn may force the cracked surfaces to close together. Assuming that the existence of the crack does not alter the temperature distribution, we can divide this problem into two parts and solve it by the principle of superposition. First, the temperature and transient thermal stress distributions along the axisymmetric surface of the imaginary cylinder without crack are obtained by finite element implicit time integration method Secondly, the opposite sense of the stress distributions along the cracked surfaces, which is obtained previously, is treated as the traction boundary conditions; the contact length and contact pressure of the real cracked cylinder are obtained by modified elimination finite element scheme. Finally, we also obtained the normalized stress intensity factor for the crack tip of the cylinder. It is concluded that the effect due to thermoelastic coupling term on stress intensity factor becomes more important for higher coupling coefficient, and this coupling term also results in a small time lag in temperature, thermal stress and stress intensity factor.     

Generalized Magneto-thermoelastic Interaction in a Fiber-Reinforced Anisotropic Hollow Cylinder

In this article, an estimation is made to investigate the transient phenomena in the magneto-thermoelastic model in the context of the Lord and Shulman theory in a perfectly conducting medium. A finite element method is proposed to analyze the problem and obtain numerical solutions for the displacement, temperature, and radial and hoop stresses. The boundary conditions for the mechanical and Maxwell’s stresses at the internal and outer surfaces are considered. An application of a hollow cylinder is investigated where the inner surface is traction free and subjected to thermal shock, while the outer surface is traction free and thermally isolated. The displacement, incremental temperature, and the stress components are obtained and then presented graphically. Finally, the effects of the presence and absence of reinforcement on the temperature, stress, and displacement are studied.     

Thermal shock analysis and thermo-elastic stress waves in functionally graded thick hollow cylinders using analytical method

Transient stress field and thermo-elastic stress wave propagation are studied in functionally graded thick hollow cylinder under arbitrary thermo-mechanical shock loading, in this article. Thermo-mechanical properties of functionally graded (FG) cylinder are assumed to be temperature independent and vary continuously and smoothly in the radial direction. The governing dynamic equations are analytically solved in temperature and elastic fields. To solve the problem, Laplace transform is used respect to time in all constitutive equations and boundary conditions. At first, temperature field equation analytically solved using Laplace transform and series method. The dynamic behaviors of thermo-elastic stresses are illustrated and discussed for various grading patterns of thermo-mechanical properties in several points across the thickness of FG cylinder. Time history of temperature field and thermal stresses are obtained using the residual theorem and the fast Laplace inverse transform method (FLIT), respectively. Also, the effects of the cylinder thickness and convection heat transfer coefficient on dynamic response of FG cylinder are revealed and discussed. The presented analytical method provides a ground to study the time histories of radial and hoop stresses in FG cylinders with different thickness and various volume fraction exponents. The advantage of this method is its mathematical ability to support simple and complicated mathematical function for the thermo-mechanical boundary conditions. A reasonable agreement can be seen in comparison of obtained results based on the presented analytical method with published data.     

Thermoelastic state of a two-component hollow cylinder with edge radial cracks

Stationary two-dimensional axisymmetric problems of thermal conductivity and thermoelasticity for a hollow two-component cylinder with cracks are studied by the method of singular integral equations. The cross section of the cylinder has the form of a circular concentric ring with a layer of another material that also has the form of a concentric ring and contains edge radial cracks. The surfaces of the cylinder are free of stresses. Thermal processes on these surfaces are characterized by temperature conditions of the third kind. Conditions of ideal thermal and mechanical contact are satisfied on the interface of the two media. A numerical solution is obtained for the case where the inner and outer cylindrical surfaces are kept at different constant temperatures. Stress intensity factors near the tip of one or two edge cracks were found for various values of thermal and mechanical characteristics of the cylinder.Karpenko Physicomechanical Institute, Ukrainian Academy of Sciences, L'viv. Translated from Fiziko-Khimicheskaya Mekhanika Materialov, Vol. 30, No. 4, pp. 76–80, July – August, 1994.     

Deformation due to mechanical and thermal sources in generalised orthorhombic thermoelastic material

A dynamical two-dimensional problem of thermoelasticity has been considered to investigate the disturbance due to mechanical (horizontal or vertical) and thermal source in a homogeneous, thermally conducting orthorhombic material. Laplace-Fourier transforms are applied to basic equations to form a vector matrix differential equation, which is then solved by eigenvalue approach. The displacements, stresses and temperature distribution so obtained in the physical domain are computed numerically and illustrated graphically. The numerical results of these quantities for zinc crystal-like material are illustrated to compare the results for different theories of generalised thermoelasticity for an insulated boundary and a temperature gradient boundary.     

Thermal stress intensity factors of an infinite orthotropic layer with a crack

Stress intensity factors are determined for a crack in an infinite orthotropic layer. The crack is situated parallel to the plane surfaces of the layer. Stresses are solved for two kinds of the boundary conditions with respect to temperature field. In the first problem, the upper surface of the layer is heated to maintain a constant temperature T ₀, while the lower surface is cooled to maintain a constant temperature –T ₀. In the other problem, uniform heat flows perpendicular to the crack. The surfaces of the crack are assumed to be insulated. The boundary conditions are reduced to dual integral equations using the Fourier transform technique. To satisfy the boundary conditions outside the crack, the difference in temperature at the crack surfaces and differences in displacements are expanded in a series of functions that vanish outside the crack. The unknown coefficients in each series are evaluated using the Schmidt method. Stress intensity factors are then calculated numerically for a steel layer that behaves as an isotropic material and for a tyrannohex layer that behaves as an orthotropic material.     

Boundary element analysis of thermal stress intensity factors for interface griffith and cusp cracks

The thermal stress intensity factors for interface cracks of Griffith and symmetric lip cusp types under vertical uniform heat flow in a finite body are calculated by the boundary element method. The boundary conditions on the crack surfaces are insulated or fixed to constant temperature. The relationship between the stress intensity factors and the displacements on the nodal point of a crack-tip element is derived. The numerical values of the thermal stress intensity factors for an interface Griffith crack in an infinite body are compared with the previous solutions. The thermal stress intensity factors for a symmetric lip cusp interface crack in a finite body are calculated with respect to various effective crack lengths, configuration parameters, material property ratios and the thermal boundary conditions on the crack surfaces. Under the same outer boundary conditions, there are no appreciable differences in the distribution of thermal stress intensity factors with respect to each material property. However, the effect of crack surface thermal boundary conditions on the thermal stress intensity factors is considerable.     

A note on the steady state distribution of thermal stress in circular cylinders and thick walled tubes

A simple approach is given to the problem of calculating the steady state thermal stresses in a heated circular tube or cylinder. The displacement vector is expressed in terms of four harmonic functions, one of which is directly related to the temperature field. The remaining harmonics are then completely determined by the mechanical boundary conditions. The note concludes with an example illustrating the procedure.