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.     

热对流作用下筒壁涂层的边裂行为

边裂(边缘开裂)是涂层热致损伤的主要模式之一. 边缘裂纹穿透涂层后,常导致界面脱粘从而驱使涂层与基体剥离,最终丧失对基体的保护作用. 本文以热应力强度因子表征边缘裂纹的扩展驱动力,研究筒壁涂层在热对流作用下的边裂行为. 首先,利用拉普拉斯变换法,得到了瞬态温度场及热应力场的封闭解. 其次,运用Fett等的三参数法确定了筒壁涂层边缘裂纹的权函数. 最后,基于叠加原理和权函数方法计算了边缘裂纹的热应力强度因子. 探讨了无量纲时间、边缘裂纹深度、基体/涂层厚度比、热对流强度等参数对热应力强度因子的影响规律. 结果表明:热应力强度因子的峰值既非发生在热载荷初始时刻,也非发生在热稳态时刻,而出现在时间历程的中间时刻;增大热对流强度不仅可提高热应力强度因子的峰值,而且使峰值提前出现;其他条件相同时,热应力强度因子随着边缘裂纹长度的增大而降低;增大涂层厚度或减小基体厚度可增强涂层抵抗瞬态热载荷的能力.      

Transient thermal elastic fracture of a piezoelectric cylinder specimen

A finite piezoelectric cylinder with an embedded penny-shaped crack is investigated for a thermal shock load on the outer surface of the cylinder. The theory of linear electro-elasticity is applied to solve the transient temperature field and the associated thermal stresses and electrical displacements without crack. These thermal stresses and electrical displacements are added to the surfaces of the crack to form an electromechanical coupling and mixed mode boundary-value problem. The electrically permeable crack face boundary condition assumption is used, and the thermal stress intensity factor and electrical displacement intensity factor at the crack border are evaluated. The thermal shock resistance of the piezoelectric cylinder is evaluated for the analysis of piezoelectric material failure in practical engineering applications.     

Periodic heat conduction with relaxation time in cylindrical geometry

Steady-periodic heat conduction with relaxation time in an infinitely long hollow cylinder is considered. Four boundary value problems, with boundary conditions of the first and of the second kind, are solved analytically. The solution for a solid cylinder with a sinusoidally varying surface temperature is obtained as a special case of a solution found for the hollow cylinder. The effects of the relaxation time on the steady-periodic temperature field are analysed, in details, for a solid cylinder with a sinusoidally varying surface temperature and for a hollow cylinder with a sinusoidally varying heat flux at the inner surface and with a constant temperature at the outer surface. The results show that thermal resonances may occur and suggest that accurate measurements of the relaxation time could be obtained by means of experiments on steady-periodic heat conduction in cylindrical geometry. Received on 15 April 1997     

Properties and appropriate conditions of stress reduction factor and thermal shock resistance parameters for ceramics

Through introducing the analytical solution of the transient heat conduction problem of the plate with convection into the thermal stress field model of the elastic plate, the stress reduction factor is presented explicitly in its dimensionless form. A new stress reduction factor is introduced for the purpose of comparison. The properties and appropriate conditions of the stress reduction factor, the first and second thermal shock resistance (TSR) parameters for the high and low Biot numbers, respectively, and the approximation formulas for the intermediate Biot number-interval are discussed. To investigate the TSR of ceramics more accurately, it is recommended to combine the heat transfer theory with the theory of thermoelasticity or fracture mechanics or use a numerical method. The critical rupture temperature difference and the critical rup- ture dimensionless time can be used to characterize the TSR of ceramics intuitively and legibly.     

Transient temperature and thermally induced stress distributions in a partly-circumferentially heated cylindrical workpiece

This paper presents the numerical solutions of the transient temperature and thermally induced stress distributions in a partly-circumferentially heated cylindrical hollow workpiece (steel) with conjugate heat transfer. Outer surface of the workpiece is heated partly-circumferentially heat flux as its remainder outer surface is circumferentially cooled with fluid (water). Three phenomena have been considered as; (1) conduction inside the cylinder, (2) convection from the cylinder surface to the surrounding fluid, and (3) thermal stress produced by high temperature gradient inside the cylinder. The governing flow and energy equations have been solved numerically by using a control volume approach. The PHOENICS 3.2 and HEATING7 computer codes have been used for the numerical evaluation. The transient calculations have been performed individually for four fluid inlet velocities, u_i = 0.005, 0.01, 0.015 and 0.020 m/s, until the system attains steady-state. The results of this study clearly demonstrate that the temperature contours in the low inlet velocity cases are more near to a symmetric case with respect to the y = 0 plane than that in the high inlet velocity cases, and the increment of the inlet velocity exponentially reduces the temperatures and thermally induced stresses in the workpiece. The effective thermal stress differences occurring in the workpiece can be significantly reduced by the high fluid inlet velocity.     

Heat wave phenomena in solids subjected to time dependent surface heat flux

Hyperbolic heat conduction in a plane slab, infinitely long solid cylinder and solid sphere with a time dependent boundary heat flux is analytically studied. The solution is based on the separation of variables method and Duhamel’s principle. The temperature distribution, the propagation and reflection of the temperature wave and the effect of geometry on the shape of the wave front are studied for the case of a rectangular pulsed boundary heat flux. Comparisons with the solution obtained for Fourier heat conduction are performed by considering the limit of a vanishing thermal relaxation time.     

Modeling and simulation of first and second sound in solids

We investigate the propagation of mechanical and thermal waves in solids at cryogenic temperatures. The latter are known as the second sound phenomenon. It occurs, e.g., in dielectric solids and differs greatly from the classical case in which the heat transport proceeds by diffusion. Since Fourier’s law of heat conduction fails for modeling second sound, we apply a non-classical one. During the last two decades, the non-classical thermoelastodynamic theory of Green and Naghdi enjoys steadily growing research activities.     

Transient heat conduction in two-dimensional functionally graded hollow cylinder with finite length

In this paper a thick hollow cylinder with finite length made of two dimensional functionally graded material (2D-FGM) subjected to transient thermal boundary conditions is considered. The volume fraction distribution of materials, geometry and thermal boundary conditions 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 and the Crank–Nicolson finite difference method is implemented to solve time dependent equations of the heat transfer problem. Two-dimensional heat conduction in the cylinder is considered and variation of temperature with time as well as temperature distribution through the cylinder are investigated. Effects of variation of material distribution in two radial and axial directions on the temperature distribution and time response are studied. The achieved results show that using two-dimensional FGM leads to a more flexible design so that transient temperature, maximum amplitude and uniformity of temperature distributions can be modified to achieve required specifications by selecting a suitable material distribution profile in two directions.     

Heat conduction and thermal stress induced by an electric current in an infinite thin plate containing an elliptical hole with an edge crack

A uniform electric current at infinity was applied to a thin infinite conductor containing an elliptical hole with an edge crack. The electric current gives rise to two states, i.e., uniform and uneven Joule heat. These two states must be considered to analyze the heat conduction problem. The uneven Joule heat gives rise to uneven temperature and thus to heat flux, and to thermal stress.Using a rational mapping function, problems of the electric current, the Joule heat, the temperature, the heat flux, the thermal stress are analyzed, and each of their solutions is obtained as a closed form. The distributions of the electric current, the Joule heat, the temperature, the heat flux and the stress are shown in figures.The heat conduction problem is solved as a temperature boundary value problem. Solving the thermal stress problem, dislocation and rotation terms appear, which complicates this problem. The solutions of the Joule heat, the temperature, the heat flux and the thermal stress are nonlinear in the direction of the electric current. The crack problems are also analyzed, and the singular intensities at the crack tip of each problem are obtained. Mode II (sliding mode) stress intensity factor (SIF) is produced as well as Mode I (opening mode) SIF, for any direction of the electric current. The relations between the electric current density and the melting temperature and between the electric current density and SIF are investigated for some crack lengths in an aluminum plate.     

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.     

短脉冲激光加热分数阶导热及其热应力研究

短脉冲激光加热引起材料内部复杂的传热过程及热变形,现有的以Fourier定律或Cattaneo-Vernotte松弛方程结合弹性理论为框架建立起来热应力理论在刻画其热物理过程存在严重缺陷.本文基于分数阶微积分理论,以半空间为研究对象,建立了分数阶Cattaneo热传导方程和相应的热应力方程,给出了问题的初始条件和边界条件,采用拉普拉斯变换方法,给出了非高斯时间分布激光热源辐射下温度场和热应力场的解析解,研究了短脉冲激光加热的温度场及热应力场的热物理行为.数值计算中,首先对理论解进行数值验证,然后取分数阶变量p=0.5研究温度场和热应力场的变化特点及激光参数对温度和热应力的影响,最后数值计算分数阶参数对温度和热应力场的影响.计算结果表明,分数阶Cattaneo传热方程和热应力方程描述的温度和热应力任然具有波动特性,与经典的Fourier传热模型和标准的Cattaneo传热模型相比,分数阶阶次越大,热波波速越小,热波波动性越明显;反之,则热波波速越大,热扩散性越强.激光加热和冷却的速度越快,温度上升和下降的速度越快,压应力和拉应力交替变化越快,温度变化幅值越小,热应力幅值影响不明显.     

Stochastic transient analysis of thermal stresses in solids by explicit time-domain method

Stochastic heat conduction and thermal stress analysis of structures has received considerable attention in recent years. The propagation of uncertain thermal environments will lead to stochastic variations in temperature fields and thermal stresses. Therefore, it is reasonable to consider the variability of thermal environments while conducting thermal analysis. However, for ambient thermal excitations, only stationary random processes have been investigated thus far. In this study, the highly efficient explicit time-domain method(ETDM) is proposed for the analysis of non-stationary stochastic transient heat conduction and thermal stress problems. The explicit time-domain expressions of thermal responses are first constructed for a thermoelastic body. Then the statistical moments of thermal displacements and stresses can be directly obtained based on the explicit expressions of thermal responses. A numerical example involving non-stationary stochastic internal heat generation rate is investigated. The accuracy and efficiency of the proposed method are validated by comparison with the Monte-Carlo simulation.     

瞬态热冲击下层合材料板界面的热弹性行为

基于带有两个热松弛时间的G-L广义热弹性理论, 利用有限元方法研究了零阻抗理想界面层合板在瞬态热冲击诱导的位移、应力和温度等通过界面时的热弹性行为. 通过比较不同层中材料的比热容、热导系数、热松弛时间和密度等对界面处的位移、应力和温度的影响, 研究了不同材料参数对复合材料热力学行为影响, 发现不同材料参数将导致热穿过界面时界面处温度、位移和应力发生突变, 研究结果可以为由热引起的层合板挠曲变形提供理论依据.      

Similarity between unsteady-state conduction in a planar slab for short times and steady-state conduction in a uniform,straight fin

The principal goal of this study is two-fold. First, to elucidate an analogy between unsteady-state conduction in a planar slab for short times and steady-state conduction in a straight fin of uniform cross section. Second, to present approximate analytical solutions of the transient heat conduction equation for short times in a plane having a uniform initial temperature and subjected to a uniform surface temperature (Dirichlet boundary condition). Use was made of a hybrid computational method, theTransversal Method Of Lines (TMOL) to bypass the classical solution techniques for partial differential equations and exploit the physical analogy with the steady-state, heat conduction in a straight fin. The resulting quasi-steady, approximate analytical solution is very easy to employ and is suitable for obtaining quality short-time temperature distributions in the slab.     

Thermal response of a partially insulated interface crack in a graded coating-substrate structure under thermo-mechanical disturbance: Energy release and density

Consider the thermal fracture problem of a functionally graded coating-substrate structure of finite thickness with a partially insulated interface crack subjected to thermal-mechanical supply. A new model is proposed that the heat conduction through the crack region occurs and the temperature drop across the crack surfaces is the result of the thermal resistance. For the first time, real fundamental solutions are derived for the fracture analysis of functionally graded materials. The complicated mixed boundary problems of equations of heat conduction and elasticity are converted analytically into singular integral equations, which are solved numerically. The asymptotic expressions with higher order terms for the singular integral kernels are considered to improve the accuracy and efficiency of the numerical integration. Explicit expressions of various failure modes including stress intensity factors, energy release rate and strain energy density, are provided. Numerical results are presented to illustrate the effects of non-homogeneity parameters and the dimensionless thermal resistance on the temperature distribution along the crack surfaces and extended crack line, the thermal stress intensity factors and minimum strain energy density.     

无限长圆柱体中广义电磁热弹耦合问题研究

基于L-S广义热弹性理论,研究了处于磁场中无限长理想圆柱导体在边界受热冲击作用时的电磁热弹耦合问题.建立了广义电磁热弹耦合的有限元方程,为避免积分变换方法求解带来的精度丟失.采用将有限元方程直接在时间域求解的方法,得到了圆柱体中的温度、位移、应力、感应磁场和感应电场的分布规律,反映了热的波动性及电磁热弹的耦合效应.结果表明,将有限元方程直接在时间域求解,可以获得各物理量的准确分布.得到温度在热波波前处的阶跃,准确地反应热波的波动效应.     

Practical Green's function for the thermal stress field induced by a heat source in plane thermoelasticity

The classical Green’s functions used in the literature for a heat source in a homogeneous elastic medium cannot lead to finite remote thermal stresses in the medium, so that they may not work well in practical thermal stress analyses. In this paper, we develop a practical Green’s function for a heat source disposed eccentrically into an elastic disk/cylinder subject to plane deformation. The edge of the disk/cylinder is assumed to be thermally permeable and traction-free. The full thermal stress field induced by the heat source in the disk/cylinder is determined exactly and explicitly via the Cauchy integral techniques. In particular, a very simple formula is obtained to describe the hoop thermal stress on the edge of the disk/cylinder, which may be conveniently useful for analyzing the thermal stresses in microelectronic components.     

Time-fractional radial heat conduction in a cylinder and associated thermal stresses

The theory of thermal stresses based on the heat conduction equation with the Caputo time-fractional derivative of order 0 < α ≤ 2 is used to investigate axisymmetic thermal stresses in a cylinder. The solution is obtained applying the Laplace and finite Hankel integral transforms. The Dirichlet and two types of Neumann problems with the prescribed boundary value of the temperature, the normal derivative of the temperature, and the heat flux are considered. Numerical results are illustrated graphically.     

Stochastic thermoelastic problem of a functionally graded plate under random temperature load

This study attempts to derive the statistics of temperature and thermal stress in functionally graded material (FGM) plates exposed to random external temperatures. The thermomechanical properties of the FGM plates are assumed to vary arbitrarily only in the plate thickness direction. The external temperatures are expressed as random functions with respect to time. The transient temperature field in the FGM plate is determined by solving a nonhomogeneous heat conduction problem for a multilayered plate with linear nonhomogeneous thermal conductivity and different homogeneous heat capacity in each layer. The autocorrelations and power spectrum densities (PSDs) of temperature and thermal stress are derived analytically. These statistics for FGM plates composed of partially stabilised zirconia (PSZ) and austenitic stainless steel (SUS304) are computed under the condition that the fluctuation in the external temperature can be considered as white noise or a stationary Markov process.