Nonlinear transient heat conduction analysis of functionally graded materials in the presence of heat sources using an improved meshless radial point interpolation method

An improved meshless radial point interpolation method, for the analysis of nonlinear transient heat conduction problems is proposed. This method is implemented for the heat conduction analysis of functionally graded materials (FGMs) with non-homogenous and/or temperature dependent heat sources. The conventional meshless RPIM is an appropriate numerical technique for the analysis of engineering problems. One advantage of this method is that it is based on the global weak formulation, and also the associated shape functions possess the Kronecker delta function property. However, in the original form, the evaluation of the global domain integrals requires the use of a background mesh. The proposed method benefits from a meshless integration technique, which has the capability of evaluating domain integrals with a better accuracy and speed in comparison with the conventional integration methods, and therefore a truly meshless technique is attained. This integration technique is especially designed for the fast and accurate evaluation of several domain integrals, with different integrands, over a single domain. Some 2D and 3D examples are provided to assess the efficiency of the proposed method.     

Mathematical modelling of heat conduction in certain functionally graded composites

The aim of contribution is to formulate a certain extended version of the tolerance modelling technique for functionally graded composites. For the sake of simplicity, the considerations are restricted to the bidirectionally graded heat conductors. It is shown that the proposed approach enables to determine an entire class of mathematical models for which applications can be found in various specific problems. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Transient thermoelastic response in a cracked strip of functionally graded materials via generalized fractional heat conduction

This work is devoted to analyzing a thermal shock problem of an elastic strip made of functionally graded materials containing a crack parallel to the free surface based on a generalized fractional heat conduction theory. The embedded crack is assumed to be insulated. The Fourier transform and the Laplace transform are employed to solve a mixed initial-boundary value problem associated with a time-fractional partial differential equation. Temperature and thermal stresses in the Laplace transform domain are evaluated by solving a system of singular integral equations. Numerical results of the thermoelastic fields in the time domain are given by applying a numerical inversion of the Laplace transform. The temperature jump between the upper and lower crack faces and the thermal stress intensity factors at the crack tips are illustrated graphically, and phase lags of heat flux, fractional orders, and gradient index play different roles in controlling heat transfer process. A comparison of the temperature jump and thermal stress intensity factors between the non-Fourier model and the classical Fourier model is made. Numerical results show that wave-like behavior and memory effects are two significant features of the fractional Cattaneo heat conduction, which does not occur for the classical Fourier heat conduction.     

A model of the heat conduction for functionally graded laminate with uniformly distributed micro inclusions

The aim of the contribution is to formulate an asymptotic model of heat conduction for functionally graded two component laminate reinforced by periodically spaced micro inclusions. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A low-order meshless model for multidimensional heat conduction problems

In this paper, a perturbation method is used to solve a two-dimensional unsteady heat conduction problem. Low-order transfer functions are defined. Step responses are obtained and compared to the complete numerical solutions given by a meshless method. The analytical results are found to be in good agreement with numerical solutions which reveals the effectiveness and convenience of the used method.     

A CQM-based BEM for transient heat conduction problems in homogeneous materials and FGMs

A method for numerical solution of time-domain boundary integral formulations of transient problems governed by the heat equation is presented. The heat conduction problem is analyzed considering homogeneous and non-homogeneous media. In the case of the non-homogeneous media, the conductor material is assumed to be a functionally graded material, i.e., the material properties vary spatially according to known smooth functions. For some specific spatial variations of the material properties, the fundamental solution and the boundary integral equation of the problem are obtained thanks to a change of variables that transforms the original problem to the standard heat conduction problem for homogeneous materials. For the treatment of time-dependent terms, the convolution quadrature method is adopted to approximate numerically the integral equation of the time-domain boundary element method. In the case that the responses are required at a large number of interior points, the convolution performed to calculate them is very time consuming. It is shown that the discrete convolution of the proposed formulation can be computed by means of the fast Fourier transform technique, which considerably reduces the computational complexity. Results for some transient heat conduction examples are presented to validate the numerical techniques studied.     

Stochastic heat conduction analysis of a functionally graded annular disc with spatially random heat transfer coefficients

The mean and variance of the temperature are analytically obtained in a functionally graded annular disc with spatially random heat transfer coefficients (HTCs) on the upper and lower surfaces. This annular disc has arbitrary variations in the HTCs (i.e., arbitrary thermal interaction with the surroundings) and gradient material composition only along the radial direction and is subjected to deterministic axisymmetrical heating at the lateral surfaces. The stochastic temperature field is analysed by considering the annular disc to be multilayered with spatially constant material properties and spatially constant but random HTCs in each layer. A type of integral transform method and a perturbation method are employed in order to obtain the analytical solutions for the statistics. The correlation coefficients of the random HTCs are expressed in the form of a linear function with respect to the radial distance as a non-homogeneous random field of discrete space. Numerical calculations are performed for functionally graded annular discs composed of stainless steel and ceramic, which comprise two types of material composition distributions. The effects of the magnitude of the means of HTCs, volume fraction distributions of the constitutive materials and correlation strengths of the HTCs on the standard deviation of the temperature are discussed.     

Contact problems for functionally graded materials of complicated structure

An approximate analytical method allowing one to efficiently solve, to a preassigned accuracy, contact problems for materials with properties arbitrarily varying in depth is developed. Its possibilities are illustrated with the example of torsion of an elastic half-space, having a coating inhomogeneous across its thickness, by a circular stamp. All the results obtained are rigorously substantiated. For the approximate solutions constructed, their error is analyzed. The asymptotic properties of the solutions are investigated. The cases of a nonmonotonic change in the elastic properties are considered. In particular, the analytical solutions are examined in the case where the variation gradient of the elastic properties changes its sign many times. The results derived allow one to solve the inverse problems of elasticity theory of inhomogeneous media (e.g., the problem on controlling the variation in the elastic properties of a covering across its thickness).     

Simulation of brittle crack growth in functionally graded materials

We consider a thermodynamic consistent framework for crack propagation by applying a dissipation inequality to a time dependent migrating control volume. The direction of crack growth is obtained in terms of material forces as a result of the principle of maximum dissipation. In the numerical implementation a staggered algorithm – deformation update for fixed geometry followed by geometry update for fixed deformation – is employed within each time increment. The corresponding mesh is generated by combining Delaunay triangulation with local mesh refinement. A numerical example with inhomogeneous material properties illustrates the capability of the resulting algorithm. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Time-space boundary elements for transient heat conduction problems

This paper is concerned with a generalized time-space boundary element formulation for transient heat conduction problems in anisotropic media. A weighted residual form of the governing equation is used to obtain the boundary integral equation in terms of the fundamental solution. The resulting boundary integral equation is discretized by means of a wide variety of boundary elements from constant-elements to higher-order isoparametric elements located both in time and space.     

Some problems of heat conduction in transversally graded composites

The object of considerations is a two-component layer made of conductors non-periodically distributed in the form of laminas along the layer thickness. It is assumed that the distribution of the macroscopic properties of this laminate is approximated by continuous slowly-varying functions across laminas. Media of this kind can be treated as made of a functionally graded material. The aim of the paper is to apply the tolerance model, [8], to analyse one-directional, non-stationary heat conduction along the axis perpendicular to laminas. Moreover, received results are compared to the solutions obtained in the framework of the higher-order theory, [1]. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A boundary element approach for solving plane elastostatic equations of anisotropic functionally graded materials

A boundary element method is proposed for the numerical solution of an important class of boundary value problems governed by plane elastostatic equations of anisotropic functionally graded materials. The grading function of the material properties may be any general function that varies smoothly from point to point in the material. The proposed boundary element method is applied to solve some specific problems to check its validity and accuracy.     

Two collinear antiplane cracks in functionally graded magnetoelectroelastic composite materials

An analysis of two collinear antiplane cracks in inhomogeneous (functionally graded) anisotropic magnetoelectroelastic materials is presented. In designing components involving functionally graded materials, an important aspect of the problem is the fracture failure. The problem is formulated for transversely isotropic functionally graded magnetoelectroelastic materials. An integral transform is employed to reduce the problem to a singular integral equation that can be solved.     

Enriched finite element spaces for transient conduction heat transfer

This paper presents an alternative approach via finite elements to treat numerically the thermal shocks in heat transfer finite element analysis. The method consists in using the standard enriched finite element approaches with time-interpolation. It will be applied here to the transient conduction heat equation where the classical Galerkin method is shown to be unstable. The proposed method consists in adding and eliminating bubbles to the finite element space and then to interpolate the solution to the real time step. This modification is equivalent to the addition of a stabilizing term tuned by a local time-dependent stability parameter, which ensures an oscillating-free solution. To validate this approach, the numerical results obtained in classical 2D and 3D benchmark problems are compared with the Galerkin and the analytical solutions.     

Periodically varying heat source response in a functionally graded microelongated medium

The present paper deals with the response due to periodically varying heat sources in the neighborhood of the origin of a functionally graded isotropic unbounded microelongated medium, in the context of generalized thermoelastic theory. The expressions for displacement, microelongation and temperature fields have been obtained in Laplace-Fourier transformed domain. After computing the inverse Fourier transforms by contour integration technique, the inversion of Laplace transforms has been obtained numerically. The changes of displacement, microelongation, and normal strain have been shown graphically for different types of heat source.     

End effects for a generalized biharmonic equation with applications to functionally graded materials

In this paper we consider a generalized biharmonic equation modeling two-dimensional inhomogeneous elastic state in the curvilinear rectangle a<r<b, 0<θ<α, where (r,θ) denote plane polar coordinates. Such an arch-like region is maintained in equilibrium under self-equilibrated traction applied on one of the edges, while the other three edges are traction free. Our aim is to derive some explicit spatial estimates describing how some appropriate measures concerning the specific Airy stress function evolve with respect to the distance to the loaded edge. Two types of smoothly varying inhomogeneity are considered: (i) the elastic moduli vary smoothly with the polar distance, (ii) they vary smoothly with the polar angle. Such types of smoothly varying inhomogeneous elastic materials provide a model for technological important functionally graded materials. The results of the present paper prove how the spatial decay rate varies with the constitutive profile.     

A BEM for transient thermoelastic analysis of a functionally graded layer on a homogeneous substrate under thermal shock

Transient thermoelastic analysis of isotropic and linear thermoelastic bimaterials, which are constituted by a functionally graded (FG) layer attached to a homogeneous substrate, subjected to thermal shock is presented in this paper. For this purpose, a boundary element method for transient linear coupled thermoelasticity is developed. The material properties of the FG layer are assumed to be continuous functions of the spatial coordinates. The boundary-domain integral equations are derived by using the fundamental solutions of linear coupled thermoelasticity for the corresponding isotropic, homogeneous and linear thermoelastic solids in the Laplace-transformed domain. For the numerical solution, a collocation method with piecewise quadratic approximation is implemented. Numerical results for the dynamic stress intensity factors are presented and discussed. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Application of a hybrid meshless technique for natural frequencies analysis in functionally graded thick hollow cylinder subjected to suddenly thermal loading

A hybrid meshless technique based on composition of meshless local Petrov–Galerkin method (for spatial variables) and Newmark finite difference method (for time domain) is developed for natural frequencies analysis of thick cylinder made of functionally graded materials (FGMs). The FG cylinder is assumed to be under suddenly thermal loading, axisymmetric and plane strain conditions. The dynamic behaviors and time history of displacements are obtained in time domain using Green–Naghdi (GN) theory of coupled thermo-elasticity (without energy dissipation). Using fast Fourier transform (FFT) technique, the displacements are transferred to frequency domain and all natural frequencies are illustrated for various grading patterns of FGMs. The variations of mechanical properties in FG thick hollow cylinder are considered to be in nonlinear volume fraction law through radial direction. The presented hybrid meshless technique furnishes a ground to analyze the effects of various grading patterns of FGMs on natural frequencies, which are obtained employing GN coupled thermo-elasticity governing equations. Also, the frequency history and natural frequencies are illustrated for various grading patterns at several points across thickness of cylinder.