Application of fractional order theory of thermoelasticity to 3D time-dependent thermal shock problem for a half-space

A three-dimensional model of thermoelasticity with fractional order heat transfer is established. The resulting nondimensional coupled equations together with the Laplace and double Fourier transforms techniques are applied to a half space, which is assumed to be traction free and subjected to a thermal shock that is a function of time. The inverses of Fourier transforms and Laplace transforms are obtained numerically. Numerical results for the temperature, the stress, the strain, and the displacement distributions are represented graphically. The predictions of the theory are discussed and compared with those for the coupled and generalized theories of thermoelasticity.     

Dynamic response of a generalized piezoelectric-thermoelastic problem under fractional order theory of thermoelasticity

The dynamic response of a piezoelectric-thermoelastic rod made of piezoelectric ceramics (PZT-4) subjected to a moving heat source is dealt with in the context of the fractional order theory of thermoelasticity. The piezoelectric-thermoelastic governing equations for the problem are formulated and then solved by means of Laplace transform together with its numerical inversion. The distributions of the considered nondimensional temperature, displacement, stress, and electric potential are obtained and illustrated graphically. The effects of fractional order parameter and the velocity of heat source on the variations of the considered variables are investigated, and the results show that they have significant influence on the variations of the considered variables. The present investigation could be helpful for better understanding the multi-field coupling effect of mechanical, electric, and thermal fields in real piezoelectric ceramics structures, and provide some guidelines in the optimal design of actuators or sensors made of piezoelectric ceramics serving in a thermoelastic environment.     

Two-dimensional problem for a half-space with axi-symmetric distribution in the theory of generalized thermoelastic diffusion

In this work, we study a two-dimensional problem of axi-symmetric distribution of temperatures in a half-space with a permeating substance in contact with the bounding plane in the context of the theory of generalized thermoelastic diffusion with one relaxation time. The surface of the half-space is taken as traction free and subjected to axi-symmetric time-dependent thermal shock. The chemical potential is also assumed to be a known function of time on the bounding plane. The Laplace and Hankel transform techniques are used. The analytical solution in the transform domain is obtained by using a direct approach. The inverse of the double transform is obtained by using a numerical method based on Fourier expansion techniques. Numerical results for the temperature, displacement, stress, concentration, and chemical potential are carried out and represented graphically.     

Stress singularity analysis for orthotropic V‐notches in the generalised plane strain state

The singularity for the V‐notch under the generalised plane deformation is investigated by the combination of the asymptotic analysis with the interpolating matrix method developed by part of the authors before. The displacement asymptotic expansions at the vicinity of the V‐notch vertex are introduced into the equilibrium equations, which are transformed into a set of characteristic ordinary differential equations with respect to the notch singularity orders. The boundary conditions and interfacial compatibility conditions are also represented by the combination of the singularity orders and characteristic angular functions. The determination of the singularity orders and characteristic angular functions are transformed into solving the ordinary differential equations with variable coefficients, which are solved by the interpolating matrix method. The present method is suitable for the singularity analysis for isotropic and orthotropic V‐notches. It is versatile for analysing the stress singularity of single material V‐notches, bi‐material V‐notches, interface edges and cracks. The correctness of the results by the proposed method is ensured by the comparison with the published ones.     

Boundary element analysis for axisymmetric heat conduction and thermal stress in steady state

In this paper, axisymmetric heat conduction and thermal stress problems with three types of boundary conditions are analysed by the boundary element method. The temperature and thermal stress fields for the piston of a diesel engine are calculated using triangular finite elements and constant boundary elements, respectively, and the two results agree. However, BEM needs fewer data, less computer time (about one-sixth that of FEM) and storage volume. The advantages of BEM are sufficiently demonstrated.     

Asymptotic approach to transient thermal shock problem with variable material properties

In this paper, an asymptotic approach is proposed to solve the transient thermal shock problem with variable material properties. The governing equations of isotropic elastic medium with temperature-dependent properties are derived in the context of the Lord–Shulman generalized theory of thermoelasticity, where the higher-order expansion with respected to increment temperature of the Helmholtz free energy is used to describe the relation of each material parameter with real temperature. Then, the layer method is used to deal with the variation of material properties with real temperature, and a system with discrete linear equations is obtained by ignoring some higher-order quantities. This system is then solved analytically by the integral transform method, where the Laplace transform technique and its limit theorem is employed to deal with these linear partial differential equations. This asymptotic approach is applied to solve the thermoelastic response of a thin plate with variable material properties. The asymptotic solutions of the displacement, temperature and stresses, induced by a sudden temperature rise at the boundary of thin plate, are obtained. The corresponding results for each physical field are discussed, as well as the comparison with the results obtained from the case with constant properties is also conducted to evaluate the effect of temperature dependency of material properties.     

基于分数阶热弹性理论的含有球型空腔无限大体的热冲击动态响应

基于新近提出的分数阶广义热弹性理论,研究了含有球型空腔的无限大体受热冲击作用时的动态响应。该文给出分数阶广义热弹性理论下的控制方程,通过拉普拉斯积分变换及其数值反变换对控制方程进行了求解,得到了带有球型空腔无限大体中的无量纲温度、位移、径向应力和环向应力等物理量的分布规律。计算中重点研究了分数阶参数对各物理量的影响效应。结果表明:含有球腔的无限大体内由于热冲击而出现了热弹耦合效应;分数阶参数显著地影响各物理量的分布规律。     

超急速传热过程中热弹性响应的解析分析

考虑超急速传热过程中诱发的热冲击效应,基于L-S广义热弹性理论,建立了温度突变加热条件下热弹性响应的控制方程组。借助于Laplace正逆变换,在适当简化的条件下推导了一维超急速传热问题热弹性响应的解析表达式。通过对温度场、位移场及应力场的解析求解,给出了超急速传热过程中热波和热弹性波在弹性体内的传递规律,并指出在超急速传热条件下,由于热波和热弹性波的相互叠加作用削弱了热作用产生的热冲击效应。     

Model of Fractional Heat Conduction in a Thermoelastic Thin Slim Strip under Thermal Shock and Temperature-Dependent Thermal Conductivity

The present paper paper, we estimate the theory of thermoelasticity a thin slim strip under the variable thermal conductivity in the fractional-order form is solved. Thermal stress theory considering the equation of heat conduction based on the time-fractional derivative of Caputo of order α is applied to obtain a solution. We assumed that the strip surface is to be free from traction and impacted by a thermal shock. The transform of Laplace (LT) and numerical inversion techniques of Laplace were considered for solving the governing basic equations. The inverse of the LT was applied in a numerical manner considering the Fourier expansion technique. The numerical results for the physical variables were calculated numerically and displayed via graphs. The parameter of fractional order effect and variation of thermal conductivity on the displacement, stress, and temperature were investigated and compared with the results of previous studies. The results indicated the strong effect of the external parameters, especially the time-fractional derivative parameter on a thermoelastic thin slim strip phenomenon.     

The method of fundamental solutions for heat conduction in layered materials

In this paper, we investigate the application of the Method of Fundamental Solutions (MFS) to two‐dimensional problems of steady‐state heat conduction in isotropic and anisotropic bimaterials. Two approaches are used: a domain decomposition technique and a single‐domain approach in which modified fundamental solutions are employed. The modified fundamental solutions satisfy the interface continuity conditions automatically for planar interfaces. The two approaches are tested and compared on several test problems and their relative merits and disadvantages discussed. Finally, we use the domain decomposition approach to investigate bimaterial problems where the interface is non‐planar and the modified fundamental solutions cannot be used. Copyright © 1999 John Wiley & Sons, Ltd.     

聚合物基复合材料的导热性能研究

讨论了聚合物基高导热高绝缘纳米复合材料的导热机理与常用的导热理论模型。考虑到填充率、温度等的影响，用不同的理论模型计算了氧化铝纳米颗粒填充环氧树脂的热导率，并结合相关研究实验对不同的导热理论模型进行分析比较。     

Multi-material topology optimization for the transient heat conduction problem using a sequential quadratic programming algorithm

Transient heat conduction analysis involves extensive computational cost. It becomes more serious for multi-material topology optimization, in which many design variables are involved and hundreds of iterations are usually required for convergence. This article aims to provide an efficient quadratic approximation for multi-material topology optimization of transient heat conduction problems. Reciprocal-type variables, instead of relative densities, are introduced as design variables. The sequential quadratic programming approach with explicit Hessians can be utilized as the optimizer for the computationally demanding optimization problem, by setting up a sequence of quadratic programs, in which the thermal compliance and weight can be explicitly approximated by the first and second order Taylor series expansion in terms of design variables. Numerical examples show clearly that the present approach can achieve better performance in terms of computational efficiency and iteration number than the solid isotropic material with penalization method solved by the commonly used method of moving asymptotes. In addition, a more lightweight design can be achieved by using multi-phase materials for the transient heat conductive problem, which demonstrates the necessity for multi-material topology optimization.     

Boundary element method analyses of transient heat conduction in an unbounded solid layer containing inclusions

The boundary element method (BEM) is used to compute the three-dimensional transient heat conduction through an unbounded solid layer that may contain heterogeneities, when a pointwise heat source placed at some point in the media is excited. Analytical solutions for the steady-state response of this solid layer when subjected to a spatially sinusoidal harmonic heat line source are presented when the solid layer has no inclusions. These solutions are incorporated into a BEM formulation as Greens functions to avoid the discretization of flat media interfaces. The solution is obtained in the frequency domain, and time responses are computed by applying inverse (Fast) Fourier Transforms. Complex frequencies are used to prevent the aliasing phenomena. The results provided by the proposed Greens functions and BEM formulation are implemented and compared with those computed by a BEM code that uses the Greens functions for an unbounded media which requires the discretization of all solid interfaces with boundary elements. The proposed BEM model is then used to evaluate the temperature field evolution through an unbounded solid layer that contains cylindrical inclusions with different thermal properties, when illuminated by a plane heat source. In this model zero initial conditions are assumed. Different simulation analyses using this model are then performed to evaluate the importance of the thermal properties of the inclusions on transient heat conduction through the solid layer.     

Simulation of transient heat conduction using one‐dimensional mapped infinite elements

Many engineering problems exist in physical domains that can be said to be infinitely large. A common problem in the simulation of these unbounded domains is that a balance must be met between a practically sized mesh and the accuracy of the solution. In transient applications, developing an appropriate mesh size becomes increasingly difficult as time marches forward. The concept of the infinite element was introduced and implemented for elliptic and for parabolic problems using exponential decay functions. This paper presents a different methodology for modeling transient heat conduction using a simplified mesh consisting of only two‐node, one‐dimensional infinite elements for diffusion into an unbounded domain and is shown to be applicable for multi‐dimensional problems. A brief review of infinite elements applied to static and transient problems is presented. A transient infinite element is presented in which the element length is time‐dependent such that it provides the optimal solution at each time step. The element is validated against the exact solution for constant surface heat flux into an infinite half‐space and then applied to the problem of heat loss in thermal reservoirs. The methodology presented accurately models these phenomena and presents an alternative methodology for modeling heat loss in thermal reservoirs. Copyright © 2010 John Wiley & Sons, Ltd.     

Investigation on a thermo-piezoelectric problem with temperature-dependent properties under fractional order theory of thermoelasticity

The dynamic response of a one-dimensional thermo-piezoelectric problem with variable material properties is investigated in the context of fractional order theory of thermoelasticity. The system of governing equations for the thermo-piezoelectric problem is formulated and then solved by means of Laplace transform together with its numerical inversion. The distributions of the considered non-dimensional displacement, temperature, and stress are obtained and illustrated graphically. The effects of fractional order parameter and temperature-dependent properties on the considered quantities are investigated and the results show that they have significant influence on the variations of the considered quantities.     

瞬态热传导温度场不确定性区间数值分析

针对具有区间参数的瞬态热传导问题,建立了一种瞬态热传导温度场不确定性分析的求解模式。在空间上,采用八节点等参元技术进行离散,并且考虑了非均质和参数分布的影响。利用基于单元分析的区间有限元法和矩阵摄动理论,将不确定性系统进行分离,分别建立确定性部分和不确定性部分的有限元模型。时间域上,引入时域精细算法进行离散求解,进而得到瞬态温度场响应的区间范围。数值算例取得了满意的结果,结果表明所提求解模式的可行性和有效性。     

A 2D problem of time-fractional heat order for two-temperature thermoelasticity under hydrostatic initial stress

The present study solves the problem of thermoelastic interactions in a half-space medium under hydrostatic initial stress in the context of a fractional order heat conduction model with two-temperature theory. The analytical solutions of the field variables are obtained by using the normal mode analysis. The obtained solutions are then applied to a specific problem for a thermally insulated surface which is acted upon by a load. The distributions of the two temperatures, displacements, and the stress components inside the half-space are studied. The graphical results depict that the fractional parameter has significant effects on all the studied field variables. Comparisons are made within the theory in the presence and absence of the hydrostatic initial stress. Thus, we can conclude that the fractional order generalized thermoelasticity model may be an improvement on studying elastic materials.     

The transient response of a functionally graded piezoelectric rod subjected to a moving heat source under fractional order theory of thermoelasticity

The transient thermo-piezoelectric response of a functionally graded piezoelectric rod subjected to a moving heat source is investigated in the context of fractional order theory of thermoelasticity proposed by Sherief. The material properties of the functionally graded piezoelectric rod are assumed to vary exponentially along the length, except for the thermal relaxation time and the specific heat, which are taken to be constant. To solve the governing equations of the problem, Laplace transform is applied, eliminating the time effect; the analytical solutions of the displacement, stress, temperature, and electric field in Laplace domain are obtained. Subsequently, the solutions of the considered variables in time domain are obtained by numerical Laplace inversion and illustrated graphically. In calculation, the effect of the fractional order parameter on the variations of the considered variables is presented.     

Analytical solution for the axisymmetric plane strain electroelastic dynamics of a special non-homogeneous piezoelectric hollow cylinder

By virtue of the introduction of a dependent variable and the separation of variables technique, the axisymmetric plane strain electroelastic dynamic problem of a special non-homogeneous piezoelectric hollow cylinder is transformed to a Volterra integral equation of the second kind about a function with respect to time, which can be solved successfully by means of the interpolation method. Then the solutions of displacements, stresses, electric displacements and electric potential are obtained. The present method is suitable for a piezoelectric hollow cylinder with an arbitrary thickness subjected to arbitrary mechanical and electrical loads. Numerical results are finally presented.     

A meshless analysis of three-dimensional transient heat conduction problems

In this paper, we consider a numerical modeling of a three-dimensional transient heat conduction problem. The modeling is carried out using a meshless reproducing kernel particle (RKPM) method. In the mathematical formulation, a variational method is employed to derive the discrete equations. The essential boundary conditions of the formulated problems are enforced by the penalty method. Compared with numerical methods based on meshes, the RKPM needs only scattered nodes, rather than having to mesh the domain of the problem. An error analysis of the RKPM for three-dimensional transient heat conduction problem is also presented in this paper. In order to demonstrate the applicability of the proposed solution procedures, numerical experiments are carried out for a few selected three-dimensional transient heat conduction problems.