Transient thermal stresses in two-dimensional functionally graded thick hollow cylinder with finite length

In this paper transient thermal stresses in a thick hollow cylinder with finite length made of two-dimensional functionally graded material (2D-FGM) based on classical theory of thermoelasticity are considered. The volume fraction distribution of materials, geometry and thermal load are assumed to be axisymmetric but not uniform along the axial direction. The finite element method with graded material properties within each element is used to model the structure. Temperature, displacements and stress distributions through the cylinder at different times are investigated. Also the effects of variation of material distribution in two radial and axial directions on the thermal stress distribution and time responses are studied. The achieved results show that using 2D-FGM leads to a more flexible design so that time responses of structure, maximum amplitude of stresses and uniformity of stress distributions can be modified to a required manner by selecting suitable material distribution profiles in two directions.     

Natural frequency analysis of 2D-FGM thick hollow cylinder based on three-dimensional elasticity equations

In this paper, natural frequencies characteristics of a thick hollow cylinder with finite length made of two-dimensional functionally graded material (2D-FGM) based on three-dimensional equations of elasticity is considered. The axisymmetric conditions are assumed for the 2D-FGM cylinder. The material properties of the cylinder are varied in the radial and axial directions with power law functions. Effects of volume fraction distribution and FGM configuration on the natural frequencies of a simply supported cylinder are analyzed. Also, the effects of length and thickness of the cylinder are considered for different material distribution profiles. Three-dimensional equations of motion are used and the eigen value problem is developed based on direct variational method. Finite element method with graded material characteristics within each element of the structure is used for solution. The study shows that the 2D-FGM cylinder exhibit interesting frequency characteristics when the constituent volume fractions and its configuration are varied.     

Temperature and stress analysis of an elastic circular cylinder in contact with heated rigid stamps

Summary  The present paper discusses a plane strain problem of transient thermoelasticity in a circular cylinder which is in partial contact with two heated rigid stamps, in the case where the coefficient of relative heat transfer on the contact surface of the cylinder is different from that on the traction-free surface. A finite difference method with respect to the time variable and Airy's thermal stress function is employed to analyze the temperature and thermoelastic fields. The problem is formulated in terms of two dual-series equations derived not only from the thermal boundary conditions but also from the mechanical boundary conditions. Since the radial, hoop and axial stresses have singularities at the end of the contact surface of the cylinder, the stress singularity coefficients are defined and then the relationship among these three coefficients is also obtained. Finally, numerical results are illustrated graphically. Received 3 March 2000; accepted for publication 12 July 2000     

Thermo-Mechanical Model and Thermal Analysis of Hollow Cylinder Planetary Roller Screw Mechanism

Planetary roller screw mechanism (PRSM) is widely used for rapid and precise motion translation from rotary into linear motion due to its high stiffness and high position accuracy. However, a high speed PRSM drive system naturally generates significant amount of frictional heat at the contact interfaces, which causes thermal deformation and thermal error and reduces motion accuracy. Preload is usually applied to remove the axial backlash of the PRSM for achieving high accuracy and great stiffness. However, more frictional heat is produced by such preload. On the other hand, larger numbers of angular contact bearings are needed to support the heavy axial load. The friction heat generated in support bearings has also to be investigated. In order to estimate the thermal distribution and thermal error of the hollow cylinder PRSM, a thermo-mechanical model based on finite element method (FEM) is developed, where heat generation from the two main sources of the PRSM, the parameters calculation of heat transfer coefficient and other thermal boundary conditions were studied. The presented model is proven capable of investigating temperature distribution, thermal error, and cooling performances of coolants of the PRSM system.     

Transient hyperbolic heat conduction in thick-walled FGM cylinders and spheres with exponentially-varying properties

This paper focuses on non-Fourier hyperbolic heat conduction analysis for heterogeneous hollow cylinders and spheres made of functionally graded material (FGM). All the material properties vary exponentially across the thickness, except for the thermal relaxation parameter which is taken to be constant. The cylinder and sphere are considered to be cylindrically and spherically symmetric, respectively, leading to one-dimensional heat conduction problems. The problems are solved analytically in the Laplace domain, and the results obtained are transformed to the real-time space using the modified Durbin’s numerical inversion method. The transient responses of temperature and heat flux are investigated for different inhomogeneity parameters and relative temperature change values. The comparisons of temperature distribution and heat flux between various time and material properties are presented in the form of graphs.     

Heat conduction and heat wave propagation in functionally graded thick hollow cylinder base on coupled thermoelasticity without energy dissipation

In this paper, heat wave propagation and coupled thermoelasticity without energy dissipation in functionally graded thick hollow cylinder is presented based on Green–Naghdi theory. The material properties are supposed to vary as a power function of radius across the thickness of cylinder. The cylinder is considered in axisymmetry and plane strain conditions and it is divided to many sub-cylinders (layers) across the thickness. Each sub-cylinder is considered to be made of isotropic material and functionally graded property can be created by suitable arrangement of layers. The Galerkin finite element method and Newmark finite difference method are employed to solve the problem. The time history of second sounds and displacement wave propagation are obtained for various values of power function. Computed results agree well with the published data.     

Non-fourier temperature field in a solid homogeneous finite hollow cylinder

Analytical solution of the non-Fourier Axisymmetric temperature field within a finite hollow cylinder is investigated considering the Cattaneo-Vernotte constitutive heat flux relation. The solution is found for the most general linear time-independent boundary conditions. The material is assumed to be homogeneous and isotropic with temperature-independent thermal properties. The standard method of separation of variables is used. The present solution can be reduced to special problems of interest by choosing appropriate boundary condition parameters. The solution is applied for two special cases including constant heat flux and the Gaussian distribution heating of a cylinder, and their respective non-Fourier thermal behavior is studied.     

高温流体与变物性固体耦合系统温度场及热变形的数值计算

在腔体内部有高温流体,其外壳等固体材料的物性温度而为化的瞬态非线性问题中,其边界条件的确定十分困难。本文采用流一固“混合流”模式进行温度场的计算,由于方程计及了瞬态、复杂形状、浮动外边界条件等因系,采用变步长的控制容积法及时域有限元进行计算。对某个高温燃气配气阀的温度场及热变形进行了数值计算,得到与实测相符的结果。     

Multiple axial cracks in a coated hollow cylinder due to thermal shock

The transient thermal stress problem of an inner-surface-coated hollow cylinder with multiple pre-existing surface cracks contained in the coating is considered. The transient temperature, induced thermal stress, and the crack tip stress intensity factor (SIF) are calculated for the cylinder via finite element method (FEM), which is exposed to convective cooling from the inner surface. As an example, the material pair of a chromium coating and an underlying steel substrate 30CrNi2MoVA is particularly evaluated. Numerical results are obtained for the stress intensity factors as a function of normalized quantities such as time, crack length, convection severity, material constants and crack spacing.     

四边形等参单元温度场有限元分析中的混合边界条件

在温度场有限元分析中，边界条件的合理确定是一个非常重要的问题。本文以四边形等参单元为基本单元，采用基于微分方程等效积分原理的Galerkin加权余量法，建立了热传导问题的有限元方程，推导出混合边界条件的有限元计算公式。最后，根据求得的温度场可得截面温度应力分布。     

Postbuckling analysis of axially loaded functionally graded cylindrical panels in thermal environments

A postbuckling analysis is presented for a functionally graded cylindrical panel of finite length subjected to axial compression in thermal environments. Material properties are assumed to be temperature dependent, and graded in the thickness direction according to a simple power law distribution in terms of the volume fractions of the constituents. The governing equations of a functionally graded cylindrical panel are based on Reddy’s higher order shear deformation shell theory with a von Kármán–Donnell-type of kinematic nonlinearity and including thermal effects. Two cases of the in-plane boundary conditions are considered. The nonlinear prebuckling deformations and initial geometric imperfections of the panel are both taken into account. A boundary layer theory of shell buckling, which includes the effects of nonlinear prebuckling deformations, large deflections in the postbuckling range, and initial geometric imperfections of the shell, is extended to the case of functionally graded cylindrical panels under axial compression. A singular perturbation technique is employed to determine the buckling loads and postbuckling equilibrium paths. The numerical illustrations concern the postbuckling behavior of axially loaded, perfect and imperfect, functional graded cylindrical panels with two constituent materials and under different sets of thermal environments. The influences played by temperature rise, volume fraction distributions, the character of in-plane boundary conditions, transverse shear deformation, panel geometric parameters, as well as initial geometric imperfections are studied.     

Non-Fourier heat conduction in a finite hollow cylinder with periodic surface heat flux

Analytical solution of the non-Fourier axisymmetric temperature field within a finite hollow cylinder exposed to a periodic boundary heat flux is investigated. The problem studied considering the Cattaneo–Vernotte (CV) constitutive heat flux relation. The material is assumed to be homogeneous and isotropic with temperature-independent thermal properties. The standard method of separation of variables is used for solving the problem with time-independent boundary conditions, and the Duhamel integral is used for applying the time dependency. The solution is applied for the special cases of harmonic uniform heat flux and an exponentially pulsed heat flux with Gaussian distribution in outer surface for modeling a laser pulse, and their respective non-Fourier thermal behavior is studied.     

瞬态耦合热弹塑性接触有限元分析方法

在弹性接触问题有限元混合法的基础上,把材料非线性和表面非线性两种迭代过程耦合,在瞬态温度场分析中将伽辽金法和向后差分法结合,并用混合法进行热接触迭代,把瞬态温度场分析和弹塑性接触分析耦合。提出了一种瞬态耦合热弹塑性接触有限元分析方法,并已成功地用于核容器的密封分析。     

Thermal stresses factors of orthotropic graded layers with a finite line bond

In the paper, the problem of a finite line bond between two orthotropic functionally graded strips under thermal loading is considered. Considering some new boundary conditions, it is assumed that the temperature drop across the finite line bond is the result of the thermal conductivity index controlling heat conduction through the bond region. Using Fourier transforms technique, the therm-elastic mixed boundary value problems are reduced to a system of singular integral equations which can be solved approximately by applying the Chebyshev polynomials. The numerical results for the temperature and displacement field as well as the thermal stress intensity factors (TSIFs) are presented. The influence of the thickness of the layers and the thermo-elastic nonhomogeneity parameters on the temperature distribution and the TSIFs is discussed in detail. These results can be expected to be used for the purpose of gaining better understanding of the thermo-mechanical behavior of layered structures.     

Interface crack problem in layered orthotropic materials under thermo-mechanical loading

In the present paper, the behavior of an interface crack for a homogeneous orthotropic strip sandwiched between two different functionally graded orthotropic materials subjected to thermal and mechanical loading is considered. It is assumed that interface crack is partly insulated, and the temperature drop across the crack surfaces is the result of the thermal resistance due to the heat conduction through the crack region. The elastic properties of the material are assumed to vary continuously along the thickness direction. The principal directions of orthotropy are parallel and perpendicular to the crack orientation. The complicated mixed boundary problems of equations of heat conduction and elasticity are converted analytically into singular integral equations, which are solved numerically. The main objective of the paper is to study the effects of material nonhomogeneity parameters and the dimensionless thermal resistance on the thermal stress intensity factors for the purpose of gaining better understanding of the thermal behavior of graded layer.     

Thermoelastic buckling behavior of thick functionally graded rectangular plates

Thermoelastic buckling behavior of thick rectangular plate made of functionally graded materials is investigated in this article. The material properties of the plate are assumed to vary continuously through the thickness of the plate according to a power-law distribution. Three types of thermal loading as uniform temperature raise, nonlinear and linear temperature distribution through the thickness of plate are considered. The coupled governing stability equations are derived based on the Reddy’s higher-order shear deformation plate theory using the energy method. The resulted stability equations are decoupled and solved analytically for the functionally graded rectangular plates with two opposite edges simply supported subjected to different types of thermal loading. A comparison of the present results with those available in the literature is carried out to establish the accuracy of the presented analytical method. The influences of power of functionally graded material, plate thickness, aspect ratio, thermal loading conditions and boundary conditions on the critical buckling temperature of aluminum/alumina functionally graded rectangular plates are investigated and discussed in detail. The critical buckling temperatures of thick functionally graded rectangular plates with various boundary conditions are reported for the first time and can be served as benchmark results for researchers to validate their numerical and analytical methods in the future.     

COUPLED THERMO-MECHANICAL ANALYSIS OF FUNCTIONALLY GRADIENT WEAK/MICRO-DISCONTINUOUS INTERFACE WITH GRADED FINITE ELEMENT METHOD

Coupled thermo-mechanical analysis of two bonded functionally graded materials subjected to thermal loads is conducted in this study with the graded finite element method. The thermal-mechanical properties of the bi-material interfaces are classified based on discontinuity degrees of their material properties and their derivatives at the interfaces. Numerical results indicate that discontinuity exerts remarkable effect on the temperature profile and stress value at the interface of two bonded functionally-graded materials. Under the thermal flux loading conditions, the stronger the interface discontinuity is, the smaller the heat flux is.     

Nonlinear transient transfinite element thermal analysis of thick-walled FGM cylinders with temperature-dependent material properties

An algorithm for investigation of nonlinear systems by the transfinite element method is presented. Basically, the transformation techniques have been developed for linear systems. Nonlinear transient heat transfer of a thick FGM cylinder with temperature-dependent material properties is investigated in the present paper to clarify the proposed algorithm. Two main novelties of the present research are: (1) incorporating the temperature-dependency of the material properties in the thermal analysis which lead to highly non-linear governing equations and (2) proposing an updating numerical transfinite element procedure to solve the resulted highly nonlinear governing equations. To reduce the effect of the artificial local heat source generation at the mutual boundaries of the elements, second order elements are used. Influences of various boundary conditions, geometric parameters, and volume fraction indices on the temperature distribution are investigated. Results of the proposed transfinite element technique show a good agreement with those obtained using the iterative time integration or analytical method. Furthermore, results reveal the significant effect of the temperature-dependency of the material properties. The present solution algorithm prevents numerical oscillations and damping, and accumulated time integration errors. The present technique may be used to obtain relatively accurate and stable results in a less computational time.     

Numerical simulation of natural convective heat transfer and fluid flow around a heated cylinder inside an enclosure

 A collocated, non-orthogonal grid based finite volume technique has been applied for investigating the two dimensional natural convective flow and heat transfer around a heated cylinder kept in a square enclosure. The effects of different enclosure wall thermal boundary conditions, fluid Prandtl number and the ratio between enclosure and cylinder dimensions (aspect ratio) upon the flow and thermal features, have been systematically studied. It is observed that the patterns of recirculatory flow and thermal stratification in the fluid are significantly modified, if any of these parameters is varied. The overall heat transfer rates are also affected due to the changes in the flow and temperature patterns. The study presents useful observations regarding the variation of local Nusselt number along each wall, for the different cases considered. Received on 2 August 2000 / Published online: 29 November 2001     

The transient piezothermoelastic problem of a thick functionally graded thermopiezoelectric strip due to nonuniform heat supply

This paper is concerned with the theoretical treatment of the transient piezothermoelastic problem involving a thick functionally graded thermopiezoelectric strip due to nonuniform heat supply in the width direction. The thermal, thermoelastic and piezoelectric constants of the strip are assumed to vary exponentially in the thickness direction. The transient two-dimensional temperature is analyzed by the methods of Laplace and finite sine transformations. We obtain the exact solution for a simply supported strip under the state of plane strain. Some numerical results for the temperature change, the displacement, the stress and electric potential distributions are presented in figures and table. Furthermore, the influence of the nonhomogeneity of the material and that of the electric boundary conditions are investigated.