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.     

Generalized thermoelastic functionally graded solid with a periodically varying heat source

This paper deals with the problem of thermoelastic interactions in a functionally graded isotropic unbounded medium due to the presence of periodically varying heat sources in the context of the linear theory of generalized thermoelasticity without energy dissipation (TEWOED). The governing equations of generalized thermoelasticity without energy dissipation (GN model type II) for a functionally graded materials (FGM) (i.e. material with spatially varying material properties)are established. The governing equations are expressed in Laplace–Fourier double transform domain and solved in that domain. Now, the inversion of the Fourier transform is carried out by using residual calculus, where poles of the integrand is obtained numerically in complex domain by using Laguerre’s method and the inversion of Laplace transform is done numerically using a method based on Fourier series expansion technique. The numerical estimates of the displacement, temperature, stress and strain are obtained for a hypothetical material. The solution to the analogous problem for homogeneous isotropic material is obtained by taking nonhomogeneity parameter suitably. Finally the results obtained are presented graphically to show the effect of nonhomogeneity on displacement, temperature, stress and strain.     

Transient thermopiezoelectric response of a one-dimensional functionally graded piezoelectric medium to a moving heat source

The thermopiezoelectricity problem of a one-dimensional (1-D), finite length, functionally graded medium excited by a moving heat source is investigated in this paper. The Lord and Shulman theory of generalized coupled thermoelasticity is employed to account for both the finite speed of thermal waves and coupling of temperature field with displacement and electric fields. Except thermal relaxation time and specific heat, which are taken to be constant for simplicity, all other properties are assumed to vary exponentially along the length through an arbitrary non-homogeneity index. Laplace transform has been used to eliminate the time effect, and three coupled fields, namely, displacement, temperature, and electric fields are obtained analytically in the Laplace domain. The solutions are then inverted to time domain using a numerical Laplace inversion method. Numerical examples are displayed to illustrate the effects of non-homogeneity index, length and thermal relaxation time on the results. When the medium is homogeneous, the results of the current paper are reduced to exactly the same results available in the literature.     

Dynamic analysis of functionally graded plates using the hybrid Fourier-Laplace transform under thermomechanical loading

In this study, the analytical solution is presented for dynamic response of a simply supported functionally graded rectangular plate subjected to a lateral thermomechanical loading. The first-order and third-order shear deformation theories and the hybrid Fourier-Laplace transform method are used. The material properties of the plate, except Poisson’s ratio, are assumed to be graded in the thickness direction according to a power-law distribution in terms of the volume fractions of the constituents. The plate is subjected to a heat flux on the bottom surface and convection on the upper surface. A third-order polynomial temperature profile is considered across the plate thickness with unknown constants. The constants are obtained by substituting the profile into the energy equation and applying the Galerkin method. The obtained temperature profile is considered along with the equations of motion. The governing partial differential equations are solved using the finite Fourier transformation method. Using the Laplace transform, the unknown variables are obtained in the Laplace domain. Applying the analytical Laplace inverse method, the solution in the time domain is derived. The computed results for static, free vibration, and dynamic problems are presented for different power law indices for a plate with simply supported boundary conditions. The results are validated with the known data reported in the literature. Furthermore, the results calculated by the analytical Laplace inversion method are compared with those obtained by the numerical Newmark method.     

Hybrid graded element model for transient heat conduction in functionally graded materials

This paper presents a hybrid graded element model for the transient heat conduction problem in functionally graded materials (FGMs). First, a Laplace transform approach is used to handle the time variable. Then, a fundamental solution in Laplace space for FGMs is constructed. Next, a hybrid graded element is formulated based on the obtained fundamental solution and a frame field. As a result, the graded properties of FGMs are naturally reflected by using the fundamental solution to interpolate the intra-element field. Further, Stefest’s algorithm is employed to convert the results in Laplace space back into the time-space domain. Finally, the performance of the proposed method is assessed by several benchmark examples. The results demonstrate well the efficiency and accuracy of the proposed method.     

黏弹性功能梯度材料裂纹问题的有限元方法

针对组分材料体积含量任意分布的黏弹性功能梯度材料裂纹问题建立有限元分析途径. 通过Laplace变换,将黏弹性问题转化到象空间中求解,基于反映材料非均匀的梯度单元和裂纹尖端奇异特性的奇异单元计算象空间中的位移、应力和应变场,应用虚拟裂纹闭合方法得到应变能释放率,分别由应力和应变能释放率确定应力强度因子. 给出这些断裂参量在物理空间和象空间之间的对应关系,由数值逆变换求出其在物理空间的相应值. 文中分析两端均匀受拉的黏弹性边裂纹板条,首先针对松弛模量表示为空间函数和时间函数乘积的特殊梯度材料进行计算,结合对应原理验证方法的有效性. 然后分析组分材料体积含量具有任意梯度分布的情形,由Mori-Tanaka方法预测象空间中的等效松弛模量. 计算结果表明,蠕变加载条件下,应变能释放率随时间增加,其增大程度与黏弹性组分材料体积含量相关. 由于梯度材料的非均匀黏弹性性质,产生应力重新分布,导致应力强度因子随时间变化,其变化范围与组分材料的体积含量分布方式有关.     

Coupled thermoelasticity of functionally graded cylindrical shells

The coupled thermoelastic response of a functionally graded circular cylindrical shell is studied. The coupled thermoelastic and the energy equations are simultaneously solved for a functionally graded axisymmetric cylindrical shell subjected to thermal shock load. A second-order shear deformation shell theory that accounts for the transverse shear strains and rotations is considered. Including the thermo-mechanical coupling and rotary inertia, a Galerkin finite element formulation in space domain and the Laplace transform in time domain are used to formulate the problem. The inverse Laplace transform is obtained using a numerical algorithm. The shell is graded through the thickness assuming a volume fraction of metal and ceramic, using a power law distribution. The results are validated with the known data in the literature.     

The Anti-plane Dynamic Fracture Analysis of Functionally Graded Materials with Arbitrary Spatial Variations of Material Properties

Anti-plane dynamic fracture analysis is presented for functionally graded materials (FGM) with arbitrary spatial variations of material properties. The FGM with the material properties varying continuously in an arbitrary manner is modeled as a multi-layered medium with the elastic modulus and mass density varying linearly in each sub-layer and continuous at the interfaces between two adjacent sub-layers. With this linearly inhomogeneous multi-layered model, the problem of a crack in a graded interfacial zone bonded to two homogeneous half-spaces or in a coating bonded to a homogeneous half-space subjected to the anti-plane shear impact load is investigated. Laplace and Fourier transforms and transfer matrix are applied to reduce the associated mixed boundary value problem to a Cauchy singular integral equation which is solved numerically in the Laplace transformed domain. The dynamic stress intensity factors (DSIF) are obtained by using the numerical technique of Laplace inversion.     

Axisymmetric non-fourier temperatures in cylindrically bounded domains

Heat conduction solutions are presented for the case where the material obeys a non-Fourier conduction law. In contrast to the Fourier law which predicts an infinite speed of heat propagation, the non-Fourier theory implies that the speed of thermal signals are finite. Axisymmetric problems for regions interior and exterior to a circular cylinder are investigated by using methods of Laplace transformation and asymptotic analysis. Comparisons of the temperature profiles are made with Fourier theory for the case of step function temperature boundary conditions.     

AN ANALYSIS OF TRANSIENT HEAT CONDUCTION FOR FUNCTIONALLY GRADED SPHERES BASED ON STATE SPACE THEORY1)

基于状态空间理论研究功能梯度圆球的球对称瞬态热传导问题。根据热传导方程和热流密度的定义,取温度场和热流密度为系统的状态向量,通过将圆球分层和在时域内应用差分格式对控制方程进行离散,建立了系统的状态方程,给出了功能梯度圆球瞬态热传导问题的半解析解。算例分析表明：本文解不但结果正确、计算效率高,而且适用于材料参数沿径向任意梯度变化的圆球瞬态热传导分析。     

Functionally graded penny-shaped cracks under dynamic loading

Dynamic load is applied to a functionally graded material with penny-shaped cracks. The materials are also transversely isotropic depending only on the axial coordinate z. The elastic region may be regarded to consist of many thin layers such that properties are constants within each layer, but they may vary from layer to layer. Laplace and Hankel transform are used in conjunction with the stiffness matrix approach. The Dual integral equations are then obtained by application of appropriate boundary and interface conditions. Stress intensity factors are then determined in the Laplace transform domain. Inversion yields the results in the time domain. Numerical examples show that multiple crack configurations in functionally graded materials can be treated where the continuously varied material properties can be divided into a finite number of layers with different properties.     

固体非傅立叶温度场的时域间断Galerkin有限元法

运用时域间断Galerkin有限元法[1],对高频非傅立叶热波动问题[2-3]进行分析。其主要特点是:取温度及温度的时间导数为基本未知量,对其分别进行3次Hermite插值和线性插值。在保证节点温度自动保持连续的基础上,温度的时间导数在离散时域存在间断。数值结果表明所提出的方法能够滤掉虚假的数值震荡,能够良好地模拟固体中的非傅立叶热波动行为。     

功能梯度材料热传导问题的近场动力学模型

当材料中存在不连续性缺陷的时候,传统的基于连续介质理论的方法在研究热传导问题时存在诸多不便．我们基于近场动力学方法(Peridynamics)建立了功能梯度材料的热传导模型,并且研究了在温度荷载作用下功能梯度材料的温度场的变化表现．该文概述了PD方法应用于热传导问题时的详细理论基础,描述了其建模思路以及计算体系,给出了使用PD方法模拟结构承受温度荷载时的计算格式,讨论了不同梯度形式对功能梯度材料内热传导的影响．算例结果表明：通过PD模型所得到的温度场与解析解吻合较好,证明了PD方法在分析功能梯度材料热传导行为等问题时的可行性．     

基于时域间断有限元方法的生物组织非傅立叶热行为数值分析

在热传导分析中,当热流与温度梯度存在时间延迟时,需采用非傅立叶热传导模型进行分析。生物组织具有较强的热松弛时间系数,承受激光、微波及烧烫等作用时,其呈现出较强的非傅立叶行为。本文对脉冲热源作用下生物组织的非傅立叶热传导进行研究,针对强脉冲引起的温度场在空间域的高梯度变化、波阵面的间断行为以及通用传统时域数值方法会带来虚假数值振荡的特点,提出采用所发展的时域间断Galerkin有限元法（DG-FEM ）进行求解计算。对多种脉冲热源作用下的非傅立叶热传导过程进行数值模拟,通过考量强脉冲作用下温度场分布和热致生物组织损伤行为的影响,表明了本文所发展的DGFEM 能够有效、准确地描述温度场空间分布和热传导过程以及非傅立叶行为下的生物热损伤更为明显,在生物组织热行为分析中应该受到重视。     

Dynamic investigation of a functionally graded layered structure with a crack crossing the interface

The dynamic response of a functionally graded layered structure with a crack crossing the interface is analyzed. The in-plane impact loading condition is considered. By using the Laplace and Fourier integral transforms, singular integral equation method and residue theory, the present problem is reduced to a singular integral equation in the Laplace transform domain. The influences of Young’s modulus ratio, thickness ratio, and crack length and location on the dynamic stress intensity factors (DSIFs) are investigated. Particularly, the DSIFs corresponding to different crack locations are shown in the case when the crack center moves from one layer to another layer through the interface. The peak and static values and overshoot characteristics of the DSIFs are analyzed. It is found that these values typically exhibit kinking behavior when the crack tips arrive at the interface. This study is different from previous other investigations in the following respects: (1) the dynamic response of a crack crossing the interface of a functionally graded structure is studied analytically, which has hardly been done in the past and (2) the present model can be reduced to some important problems, such as a functionally graded coating-substrate structure with a crack in the graded coating or homogeneous substrate or one intersecting the interface.     

The transient fracture behavior for a functionally graded layered structure subjected to an in-plane impact load

The transient fracture behavior of a functionally graded layered structure subjected to an in-plane impact load is investigated. The studied structure is composed of two homogeneous layers and a functionally graded interlayer with a crack perpendicular to the boundaries. The impact load is applied on the face of the crack. Fourier transform and Laplace transform methods are used to formulate the present problem in terms of a singular integral equation in Laplace transform domain. Considering variations of parameters such as the nonhomogeneity constant, the thickness ratio and the crack length, the dynamic stress intensity factors (DSIFs) in time domain are studied and some meaningful conclusions are obtained.The project supported by the National Science Foundation for Excellent Young Investigators (10325208), the National Natural Science Foundation of China (10432030) and the China Postdoctoral Science Foundation (2004036018)The English text was polished by Ron Marshall.     

Transient heat conduction in functionally graded thick hollow cylinders by analytical method

In this article, transient heat conduction in a cylindrical shell of functionally graded material is studied by using analytical method. The shell is assumed to be in axisymmetry conditions. The material properties are considered to be nonlinear with a power law distribution through the thickness. The temperature distribution is derived analytically by using the Bessel functions. To verify the proposed method the obtained numerical results are compared with the published results. The comparisons of temperature distribution between various time and material properties are presented.     

Hamilton体系下功能梯度梁的热冲击动力屈曲分析

在Hamilton体系下,基于Euler梁理论研究了功能梯度材料梁受热冲击载荷作用时的动力屈曲问题;将非均匀功能梯度复合材料的物性参数假设为厚度坐标的幂函数形式,采用Laplace变换法和幂级数法解析求得热冲击下功能梯度梁内的动态温度场:首先将功能梯度梁的屈曲问题归结为辛空间中系统的零本征值问题,梁的屈曲载荷与屈曲模态分别对应于Hamilton体系下的辛本征值和本征解问题,由分叉条件求得屈曲模态和屈曲热轴力,根据屈曲热轴力求解临界屈曲升温载荷。给出了热冲击载荷作用下一类非均匀梯度材料梁屈曲特性的辛方法研究过程,讨论了材料的梯度特性、结构几何参数和热冲击载荷参数对临界温度的影响。     

Analytical solution of time fractional Cattaneo heat equation for finite slab under pulse heat flux

A fractional Cattaneo model is derived for studying the heat transfer in a finite slab irradiated by a short pulse laser. The analytical solutions for the fractional Cattaneo model, the classical Cattaneo-Vernotte model, and the Fourier model are obtained with finite Fourier and Laplace transforms. The effects of the fractional order parameter and the relaxation time on the temperature fields in the finite slab are investigated. The results show that the larger the fractional order parameter, the slower the thermal wave. Moreover, the higher the relaxation time, the slower the heat flux propagates. By comparing the fractional order Cattaneo model with the classical Cattaneo-Vernotte and Fourier models, it can be found that the heat flux predicted using the fractional Cattaneo model always transports from the high temperature to the low one, which is in accord with the second law of thermodynamics. However, the classical Cattaneo-Vernotte model shows that the unphysical heat flux sometimes transports from the low temperature to the high one.     

热冲击载荷作用下的含球形空腔的广义热弹性功能梯度球形各向同性体

This paper is concerned with the determination of thermoelastic displacement, stress and temperature in a functionally graded spherically isotropic infinite elastic medium having a spherical cavity, in the context of the linear theory of generalized thermoelasticity with two relaxation time parameters （Green and Lindsay theory）. The surface of cavity is stress-free and is subjected to a time-dependent thermal shock. The basic equations have been written in the form of a vector-matrix differential equation in the Laplace transform domain, which is then solved by an eigenvalue approach. Numerical inversion of the transforms is carried out using the Bellman method. Displacement, stress and temperature are computed and presented graphically. It is found that variation in the thermo-physical properties of a material strongly influences the response to loading. A comparative study with a corresponding homogeneous material is also made.