Thermal stresses in functionally graded materials caused by a laser thermal shock

 Mathematical simulation of a thermal shock method for reliability testing of functionally graded material (FGM) is performed with the end to determine operating parameters of the testing device (power of a laser, laser beam radius, duration of heating) and to investigate the effect of the composition of FGM on a magnitude of thermal stresses in a coating. An analytical method for solution of the thermal elasticity problem is developed whereby the approach of a multilayer plate is used for determining temperature and thermal stresses distributions in a coating. We considered the limiting case of the obtained solution when the thickness of a layer is infinitesimally small and the number of layers tends to infinity. This procedure allowed us to obtain the thermal stresses distribution in a FGM coating. The results for the FGM coating composed of WC (tungsten carbide) ceramics and HS-steel are presented. It is showed that variation of the volume content of ceramics strongly affects thermal stresses in a coating and they decrease significantly in the case of the uniform spatial distribution of ceramics. Received on 21 November 2000 / Published online: 29 November 2001     

Analytical and numerical analysis for the FGM thick sphere under combined pressure and temperature loading

Thermo-mechanical analysis of functionally graded hollow sphere subjected to mechanical loads and one-dimensional steady-state thermal stresses is carried out in this study. The material properties are assumed to vary non-linearly in the radial direction, and the Poisson’s ratio is assumed constant. The temperature distribution is assumed to be a function of radius, with general thermal and mechanical boundary conditions on the inside and outside surfaces of the sphere. In the analysis presented here, the effect of non-homogeneity in FGM thick sphere was implemented by choosing a dimensionless parameter, named β _i (i = 1, . . . , 3), which could be assigned an arbitrary value affecting the stresses in the sphere. It is observed that the solution process for β _i (i = 3) = −1 are different from those obtained for other values of β _i (i = 1, . . . , 3). Cases of pressure, temperature, and combined loadings were considered separately. It is concluded that by changing the value of β _i (i = 1 . . . 3), the properties of FGM can be so modified that the lowest stress levels are reached. The present results agree well with existing results. Using FEM simulations, the analytical findings for FGM spheres under the influence of internal pressure and temperature gradient were compared to finite element results.     

旋转中心刚体-FGM梁刚柔热耦合动力学特性研究

对旋转中心刚体-功能梯度材料(functionally graded material,FGM)梁刚柔热耦合动力学特性进行研究.FGM梁为物理性能参数沿厚度方向呈幂律分布的欧拉伯努利梁.考虑柔性梁的横向弯曲变形和轴向拉伸变形, 并计入横向弯曲变形引起的纵向缩短,即非线性耦合变形量.考虑变截面空心梁在外部高温、内冷通道冷却情况下的热力耦合对系统动力学特性的影响,求解得到FGM梁沿厚度方向分布的温度场, 进而在本构关系中计入热应变.采用假设模态法对柔性梁变形场进行离散,运用第二类拉格朗日方程推导得到系统的刚柔热耦合动力学方程,并编制动力学仿真软件, 然后通过仿真算例对系统的动力学问题进行研究.结果表明:不同截面梁动力学响应差异较大, 因此需对实际系统合理建模;大范围运动已知时, 考虑热冲击载荷的FGM梁将有效抑制横向弯曲变形,而大范围运动恒定时随热冲击的叠加会出现高频振荡; 大范围运动未知时,外力矩和热冲击载荷相互作用产生热力耦合效应, 导致系统呈现高频振荡,同时与中心刚体大范围旋转运动产生刚柔热耦合效应.      

Viscous dissipation effects on thermal entrance region heat transfer in pipes with uniform wall heat flux

The analytical solution to Graetz problem with uniform wall heat flux is extended by including the viscous dissipation effect in the analysis. The analytical solution obtained reduces to that of Siegel, Sparrow and Hallman neglecting viscous dissipation as a limiting case. The sample developing temperature profiles, wall and bulk temperature distributions and the local Nusselt number variations are presented to illustrate the viscous dissipation effects. It is found that the role of viscous dissipation on thermal entrance region heat transfer is completely different for heating and cooling at wall. In the case of cooling at wall, a critical value of Brinkman number, Br _c=−11/24, exists beyond which (−11/24<Br<0) the fluid bulk temperature will always be less than the uniform entrance temperature indicating the predominance of cooling effect over the viscous heating effect. On the other hand, with Br < Br _c the bulk temperature T _b will approach the wall temperature T _w at some downstream position and from there onward the bulk temperature T _b becomes less than the wall temperature T _w with T _w > B _b > T ₀ indicating overall heating effect for the fluid. The numerical results for the case of cooling at wall Br < 0 are believed to be of some interest in the design of the proposed artctic oil pipeline.     

Stochastic analysis of a thermoelastic problem in functionally graded plates with uncertain material properties

The statistics (i.e., mean and variance) of temperature and thermal stress are analytically obtained in functionally graded material (FGM) plates with uncertainties in the thermal conductivity and coefficient of linear thermal expansion. These FGM plates are assumed to have arbitrary nonhomogeneous thermal and mechanical properties through the entire thickness of plate and are subjected to deterministic convective heating. The stochastic temperature and thermal stress fields are analysed by assuming the FGM plate is multilayered with distinct, random thermal conductivity and coefficient of linear thermal expansion in each layer. Vodicka’s method, which is a type of integral transform method, and a perturbation method are employed to obtain the analytical solutions for the statistics. The autocorrelation coefficients of each random property and cross-correlation coefficients between different random properties are expressed in exponential function forms as a non-homogeneous Markov random field of discrete space. Numerical calculations are performed for FGM plates composed of partially stabilised zirconia (PSZ) and austenitic stainless steel (SUS304), which have the largest dispersion of the random properties at the place where the volume fractions of the two constituent materials are both 0.5. The effects of the spatial change in material composition, thermal boundary condition and correlation coefficients on the standard deviations of the temperature and thermal stress are discussed.     

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

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

Mach reflection of a plane shock wave passing a mountain of 45° inclination

The transition from regular reflection (RR) to Mach reflection (MR) as a plane shock wave diffracts around a triangular mountain of 45° inclination is analysed in this paper, both by optical measurement in a shock tube and by numerical simulation the numerical method developed by Li Yingfan^[1] is of the FLIC type with triangular mesh. The dependence of the critical transition point Lk ofRR→MR on shock Mach numberM _i is analyzed and the variations of the incidence angle ω _i of the impinging shock and the reflection angle ω _r with the distanceL ^* are investigated. Our experimental and numerical results agree well with the theoretical results of Iton and Italya.     

Natural convection in a vertical rectangular enclosure with localized heating and cooling zones

Experimental and numerical studies of natural convection in a single phase, closed thermosyphon were carried out using a vertical, rectangular enclosure model. Only one vertical plate plays the role of heat transfer surface having 100 mm height and 100 mm width, and others act as the adiabatic wall made of transparent plexi-glass. The heat transfer surface is separated into three horizontal zones with an equal height; top 1/3 and bottom 1/3 of the surface are cooling and heating zones, respectively and intermediate section is an adiabatic zone. Water is used as the working fluid. Variable parameters are distance D between the heat transfer surface and an adiabatic plate opposite to the heat transfer plate, and temperature difference ΔT between heating and cooling zones. By changing both D and ΔT, three regimes of the natural convection flow; quasi-two-dimensional steady, three-dimensional steady and unsteady flows are observed by means of thermo-sensitive liquid crystal powder and numerically simulated very well by solving a set of governing equations. Received on 17 January 2000     

Transient temperature fields in functionally graded materials with different shapes under convective boundary conditions

This paper presents analyses of the transient temperature fields in an infinite plate, an infinite solid cylinder and a solid sphere made of functionally graded materials (FGMs) under convective boundary conditions. The composition and the thermo-physical properties of the infinite FGM plate, the infinite FGM solid cylinder and the FGM solid sphere are of planar symmetric, axially symmetric and spherically symmetric distributions, respectively. The analytical formulae of the one-dimensional transient temperature fields for the three FGM solids are obtained respectively by using the separation-of-variables method and the variable substitution method. Numerical results reveal that the transient temperature fields of the FGM components exhibit similar shape effect to that of homogeneous components. The present work provides valuable basis for the investigation of the thermal shock resistance of FGMs with various shapes.     

Postbuckling of FGM plates with piezoelectric actuators under thermo-electro-mechanical loadings

A postbuckling analysis is presented for a simply supported, shear deformable functionally graded plate with piezoelectric actuators subjected to the combined action of mechanical, electrical and thermal loads. The temperature field considered is assumed to be of uniform distribution over the plate surface and through the plate thickness and the electric field considered only has non-zero-valued component E_Z. The material properties of functionally graded materials (FGMs) are assumed to be graded in the thickness direction according to a simple power law distribution in terms of the volume fractions of the constituents, and the material properties of both FGM and piezoelectric layers are assumed to be temperature-dependent. The governing equations are based on a higher order shear deformation plate theory that includes thermo-piezoelectric effects. The initial geometric imperfection of the plate is taken into account. Two cases of the in-plane boundary conditions are considered. A two step perturbation technique is employed to determine buckling loads and postbuckling equilibrium paths. The numerical illustrations concern the postbuckling behavior of perfect and imperfect, geometrically mid-plane symmetric FGM plates with fully covered or embedded piezoelectric actuators under different sets of thermal and electric loading conditions. The effects played by temperature rise, volume fraction distribution, applied voltage, the character of in-plane boundary conditions, as well as initial geometric imperfections are studied.     

Influence of anisotropy and boundary conditions on the thermoelastic equilibrium of nununiformly heated laminated plates

This paper deals with three-dimensional boundary-value problems of the heat transfer and thermoelasticity theory for laminated transversally isotropic rectangular plates affected by a nonuniform temperature field. An exact analytical solution is obtained for the case of both ideal and nonideal thermal and mechanical contacts of adjacent layers. Numerical calculations were performed for a two-layer rectangular plate of arbitrary thickness. The effect of material anisotropy as well as the boundary and contact conditions on the stress-strain state of a nonuniformly heated plate is revealed. S. P. Timoshenko Institute of Mechanics, National Academy of Sciences of Ukraine, Kiev Translated from Prikladnaya Mekhanika, Vol. 35, No. 11, pp. 96–105, November, 1999.     

Thermal shock resistance behavior of a functionally graded ceramic: Effects of finite cooling rate

This work presents a semi-analytical model to explore the effects of cooling rate on the thermal shock resistance behavior of a functionally graded ceramic (FGC) plate with a periodic array of edge cracks. The FGC is assumed to be a thermally heterogeneous material with constant elastic modulus and Poisson's ratio. The cooling rate applied at the FGC surface is modeled using a linear ramp function. An integral equation method and a closed form asymptotic temperature solution are employed to compute the thermal stress intensity factor (TSIF). The thermal shock residual strength and critical thermal shock of the FGC plate are obtained using the SIF criterion. Thermal shock simulations for an Al₂O₃/Si₃N₄ FGC indicate that a finite cooling rate leads to a significantly higher critical thermal shock than that under the sudden cooling condition. The residual strength, however, is relatively insensitive to the cooling rate.     

Nonlinear bending and post-buckling of a functionally graded circular plate under mechanical and thermal loadings

Based on the classical nonlinear von Karman plate theory, axisymmetric large deflection bending of a functionally graded circular plate is investigated under mechanical, thermal and combined thermal–mechanical loadings, respectively, and axisymmetric thermal post-buckling behavior of a functionally graded circular plate is also investigated. The mechanical and thermal properties of functionally graded material (FGM) are assumed to vary continuously through the thickness of the plate, and obey a simple power law of the volume fraction of the constituents. Governing equations for the problem are derived, and then a shooting method is employed to numerically solve the equations. Effects of material constant n and boundary conditions on the temperature distribution, nonlinear bending, critical buckling temperature and thermal post-buckling behavior of the FGM plate are discussed in details.     

Stability of a fluid-saturated porous medium contained in a vertical cylinder heated from below by forced convection

A linear stability analysis determining the critical Rayleigh number R _c for onset of convection in a bounded vertical cylinder containing a fluid-saturated porous medium is performed for insulated sidewalls, isothermal top surface, and bottom surface heated by forced convection. This Newtonian heating of the bottom surface involves a Biot number Bi that allows consideration of the continuum of boundary conditions ranging from constant heat flux, with global minimum R _min=27.096 found as Bi→0, to isothermal, with global minimum R _min=4π² found as Bi→ ∞. In both cases and for most cylinder aspect ratios, incipient convection sets in as an asymmetric mode, though islands of aspect ratio exist where the onset mode is symmetric. Sample three-dimensional renderings of disturbance temperature distributions showing preferred modes at onset of convection for fixed Bi are provided and an analytical fit to R _min as a function of Bi is given.     

Free vibration of functionally graded material beams with surface-bonded piezoelectric layers in thermal environment

Free vibration of statically thermal postbuckled functionally graded material （FGM） beams with surface-bonded piezoelectric layers subject to both temperature rise and voltage is studied. By accurately considering the axial extension and based on the Euler-Bernoulli beam theory, geometrically nonlinear dynamic governing equations for FGM beams with surface-bonded piezoelectric layers subject to thermo-electro- mechanical loadings are formulated. It is assumed that the material properties of the middle FGM layer vary continuously as a power law function of the thickness coordinate, and the piezoelectric layers are isotropic and homogenous. By assuming that the amplitude of the beam vibration is small and its response is harmonic, the above mentioned non-linear partial differential equations are reduced to two sets of coupled ordinary differential equations. One is for the postbuckling, and the other is for the linear vibration of the beam superimposed upon the postbuckled configuration. Using a shooting method to solve the two sets of ordinary differential equations with fixed-fixed boundary conditions numerically, the response of postbuckling and free vibration in the vicinity of the postbuckled configuration of the beam with fixed-fixed ends and subject to transversely nonuniform heating and uniform electric field is obtained. Thermo-electric postbuckling equilibrium paths and characteristic curves of the first three natural frequencies versus the temperature, the electricity, and the material gradient parameters are plotted. It is found that the three lowest frequencies of the prebuckled beam decrease with the increase of the temperature, but those of a buckled beam increase monotonically with the temperature rise. The results also show that the tensional force produced in the piezoelectric layers by the voltage can efficiently increase the critical buckling temperature and the natural frequency.     

Thermal buckling of functionally graded circular plates based on higher order shear deformation plate theory

In this research, thermal buckling of circular plates compose of functionally graded material (FGM) is considered. Equilibrium and stability equations of a FGM circular plate under thermal loads are derived, based on the higher order shear deformation plate theory (3rd order plate theory). Assuming that the material properties vary as a power form of the thickness coordinate variable z and using the variational method, the system of fundamental partial differential equations is established. A buckling analysis of a functionally graded circular plate (FGCP) under various types of thermal loads is carried out and the result are given in closed-form solutions. The results are compared with the critical buckling temperature obtained for FGCP based on first order (1st order plate theory) and classical plate theory (0 order plate theory) given in the literature. The study concludes that higher order shear deformation theory accurately predicts the behavior of FGCP, whereas the first order and classical plate theory overestimates buckling temperature.     

Effect of the Source Term on Steady Free Convection Boundary Layer Flows over an Vertical Plate in a Porous Medium. Part I

The Darcy free convection boundary layer flow over a vertical flat plate is considered in the presence of volumetric heat generation/absorption. In the present first part of the paper it is assumed that the heat generation/absorption takes place in a self-consistent way, the source term q ^′′′≡ S of the energy equation being an analytical function of the local temperature difference T − T _∞. In a forthcoming second part, the case of the externally controlled source terms S = S(x,y ) will be considered. It is shown that due to the presence of S, the physical equivalence of the up- and downflows gets in general broken, in the sense that the free convection flow over the upward projecting hot plate (“upflow”) and over its downward projecting cold counterpart (“downflow”) in general become physically distinct. The consequences of this circumstance are examined for different forms of S. Several analytical solutions are given. Some of them describe algebraically decaying boundary layers which can also be recovered as limiting cases of exponentially decayingones. This asymptotic phenomenon is discussed in some detail.     

Cut-off cooling velocity profiling inside a keyhole model using the Boubaker polynomials expansion scheme

The time dependent heating and cooling velocities are investigated in this paper. The temperature profile is found by using a keyhole approximation for the melted zone and solving the heat transfer equation. A polynomial expansion has been deployed to determine the cooling velocity during welding cut-off stage. The maximum cooling velocity has been estimated to be V _max ≈ 83°C s⁻¹.     

Temperature dependent viscosity and thermal conductivity effects on combined heat and mass transfer in MHD three-dimensional flow over a stretching surface with Ohmic heating

An analysis has been carried out to obtain the flow, heat and mass transfer characteristics of a viscous electrically conducting fluid having temperature dependent viscosity and thermal conductivity past a continuously stretching surface, taking into account the effect of Ohmic heating. The flow is subjected to a uniform transverse magnetic field normal to the plate. The resulting governing three-dimensional equations are transformed using suitable three-dimensional transformations and then solved numerically by using fifth order Runge–Kutta–Fehlberg scheme with a modified version of the Newton–Raphson shooting method. Favorable comparisons with previously published work are obtained. The effects of the various parameters such as magnetic parameter M, the viscosity/temperature parameter θ _r , the thermal conductivity parameter S and the Eckert number Ec on the velocity, temperature, and concentration profiles, as well as the local skin-friction coefficient, local Nusselt number, and the local Sherwood number are presented graphically and in tabulated form.     

Heat transfer during heat sterilization and cooling processes of canned products

In this paper, an analysis of transient heat transfer during heat sterilization and cooling processes of a cylindrical canned product is presented. In the analysis, most practical case including the boundary condition of third kind (i.e., convection boundary condition, leading to 0.1 ≤ Bi ≤ 100) was employed. A simple analytical model for determining effective heat transfer coefficients for such products is developed. For the heat sterilization process, heating coefficient is incorporated into heat transfer coefficient model. An experimental study was performed to measure the thermal center temperatures of the short-cylindrical canned products (i.e., Tuna fish) during heat sterilization at the retort medium temperatures of 115^∘C and 121^∘C, and during cooling process at 16^∘C. The effective heat transfer coefficient model used the experimental temperature data. Using these effective heat transfer coefficients the center temperature distributions were calculated and compared with the experimental temperature distributions. Agreement was found considerably high. The results of the present study indicate that the heat-transfer analysis technique and heat-transfer coefficient model are reliable, and can provide accurate results for such problems. Received on 12 November 1997