COUPLED THERMO-MECHANICAL ANALYSIS OF FUNCTIONALLY GRADIENT WEAK/MICRO-DISCONTINUOUS INTERFACE WITH GRADED FINITE ELEMENT METHOD

Coupled thermo-mechanical analysis of two bonded functionally graded materials subjected to thermal loads is conducted in this study with the graded finite element method. The thermal-mechanical properties of the bi-material interfaces are classified based on discontinuity degrees of their material properties and their derivatives at the interfaces. Numerical results indicate that discontinuity exerts remarkable effect on the temperature profile and stress value at the interface of two bonded functionally-graded materials. Under the thermal flux loading conditions, the stronger the interface discontinuity is, the smaller the heat flux is.     

THERMAL POST-BUCKLING OF FUNCTIONALLY GRADED MATERIAL TIMOSHENKO BEAMS

Analysis of thermal post-buckling of FGM (Functionally Graded Material) Timoshenko beams subjected to transversely non-uniform temperature rise is presented. By accurately considering the axial extension and transverse shear deformation in the sense of theory of Timoshenko beam, geometrical nonlinear governing equations including seven basic unknown functions for functionally graded beams subjected to mechanical and thermal loads were formulated. In the analysis, it was assumed that the material properties of the beam vary continuously as a power function of the thickness coordinate. By using a shooting method, the obtained nonlinear boundary value problem was numerically solved and thermal buckling and post-buckling response of transversely non-uniformly heated FGM Timoshenko beams with fixed-fixed edges were obtained. Characteristic curves of the buckling deformation of the beam varying with thermal load and the power law index are plotted. The effects of material gradient property on the buckling deformation and critical temperature of beam were discussed in details. The results show that there exists the tension-bend coupling deformation in the uniformly heated beam because of the transversely non-uniform characteristic of materials.     

THERMAL FRACTURE OF FUNCTIONALLY GRADED PLATE WITH PARALLEL SURFACE CRACKS

This work examines the fracture behavior of a functionally graded material （FGM） plate containing parallel surface cracks with alternating lengths subjected to a thermal shock. The thermal stress intensity factors （TSIFs） at the tips of long and short cracks are calculated using a singular integral equation technique. The critical thermal shock △Tc that causes crack initiation is calculated using a stress intensity factor criterion. Numerical examples of TSIFs and △Tc for an Al2O3/Si3N4 FGM plate are presented to illustrate the effects of thermal property gradation, crack spacing and crack length ratio on the TSIFs and △Tc. It is found that for a given crack length ratio, the TSIFs at the tips of both long and short cracks can be reduced significantly and △Tc can be enhanced by introducing appropriate material gradation. The TSIFs also decrease dramatically with a decrease in crack spacing. The TSIF at the tips of short cracks may be higher than that for the long cracks under certain crack geometry conditions. Hence, the short cracks instead of long cracks may first start to grow under the thermal shock loading.     

Elasticity solutions for functionally graded plates in cylindrical bending

The plate theory of functionally graded materials suggested by Mian and Spencer is extended to analyze the cylindrical bending problem of a functionally graded rectangular plate subject to uniform load. The expansion formula for displacements is adopted. While keeping the assumption that the material parameters can vary along the thickness direction in an arbitrary fashion, this paper considers orthotropic materials rather than isotropic materials. In addition, the traction-free condition on the top surface is replaced with the condition of uniform load applied on the top surface. The plate theory for the particular case Of cylindrical bending is presented by considering an infinite extent in the y-direction. Effects of boundary conditions and material inhomogeneity on the static response of functionally graded plates are investigated through a numerical example.     

A MULTI-COUPLED FINITE ELEMENT ANALYSIS OF RESISTANCE SPOT WELDING PROCESS

A two-dimensional axisymmetric finite element model is developed to analyze the transient thermal and mechanical behaviors of the Resistance Spot Welding (RSW) process using commercial software ANSYS. Firstly a direct-coupled electrical-thermal Finite Element Analysis (FEA) is performed to analyze the transient thermal characteristics of the RSW process. Then based on the thermal results a sequential coupled thermo-elastic-plastic analysis is conducted to determine the mechanical features of the RSW process. The thermal history of the whole process and the temperature distribution of the weldment are obtained through the analysis. The mechanical features, including the distributions of the contact pressure at both the faying surface and the electrode-workpiece interface, the stress and strain distributions in the weldment and their changes during the RSW process, the deformation of the weldment and the electrode displacement are also calculated.     

Transient thermal cracking associated with non-classical heat conduction in cylindrical coordinate system

This paper studies the fracture behavior of a thermoelastic cylinder subjected to a sudden temperature change on its outer surface within the framework of non-classical heat conduction.The heat conduction equation is solved by separation of variable technique.Closed form solution for the temperature field and the associated thermal stress are established.The critical parameter governing the level of the transient thermal stress is identified.Exact expression for the transient stress intensity factor is obtained for a crack in the cylinder.The difference between the non-classical solutions and the classical solution are discussed.It is found that the traditional classical heat conduction considerably underestimates the transient thermal stress and thermal stress intensity factor.     

Effect of Stress-Dependent Thermal Conductivity on Thermo-Mechanical Coupling Behavior in GaN-Based Nanofilm Under Pulse Heat Source

The thermal properties of a nanostructured semiconductor are affected by multi-physical fields,such as stress and electromagnetic fields,causing changes in temperature and strain distributions.In this work,the influence of stress-dependent thermal conductivity on the heat transfer behavior of a GaN-based nanofilm is investigated.The finite element method is adopted to simulate the temperature distribution in a prestressed nanofilm under heat pulses.Numerical results demonstrate the effect of stress field on the thermal conductivity of GaN-based nanofilm,namely,the prestress and the thermal stress lead to a change in the heat transfer behavior in the nanofilm.Under the same heat source,the peak temperature of the film with stress-dependent thermal conductivity is significantly lower than that of the film with a constant thermal conductivity and the maximum temperature difference can reach 8.2 K.These results could be useful for designing GaN-based semiconductor devices with higher reliability under multi-physical fields.     

Interaction between heat dipole and circular interfacial crack

The heat dipole consists of a heat source and a heat sink. The problem of an interracial crack of a composite containing a circular inclusion under a heat dipole is investigated by using the analytical extension technique, the generalized Liouville theorem, and the Muskhelishvili boundary value theory. Temperature and stress fields are formulated. The effects of the temperature field and the inhomogeneity on the interracial fracture axe analyzed. As a numerical illustration, the thermal stress intensity factors of the interfacial crack are presented for various material combinations and different positions of the heat dipole. The characteristics of the interfacial crack depend on the elasticity, the thermal property of the composite, and the condition of the dipole.     

Conduction-radiation effect on transient natural convection with thermophoresis

This paper discusses the effect of thermophoretic particle deposition on the transient natural convection laminar flow along a vertical flat surface,which is immersed in an optically dense gray fluid in the presence of thermal radiation.In the analysis,the radiative heat flux term is expressed by adopting the Rosseland diffusion approximation.The governing equations are reduced to a set of parabolic partial differential equations.Then,these equations are solved numerically with a finite-difference scheme in the entire time regime.The asymptotic solutions are also obtained for sufficiently small and large time.The obtained asymptotic solutions are then compared with the numerical solutions,and they are found in excellent agreement.Moreover,the effects of different physical parameters,i.e.,the thermal radiation parameter,the surface temperature parameter,and the thermophoretic parameter,on the transient surface shear stress,the rate of surface heat transfer,and the rate of species concentration,as well as the transient velocity,temperature,and concentration profiles are shown graphically for a fluid(i.e.,air) with the Prandtl number of 0.7 at 20 C and 1.013 × 10 5 Pa.     

Nanoparticle aggregation effects on radiative heat transport of nanoliquid over a vertical cylinder with sensitivity analysis

A sensitivity analysis is performed to analyze the effects of the nanoparticle(NP)aggregation and thermal radiation on heat transport of the nanoliquids(titania based on ethylene glycol)over a vertical cylinder.The optimization of heat transfer rate and friction factor is performed for NP volume fraction(1%≤φ≤3%),radiation parameter(1≤R_t≤3),and mixed convection parameter(1.5≤λ≤2.5)via the facecentered central composite design(CCD)and the response surface methodology(RSM).The modified Krieger and Dougherty model(MKDM)for dynamic viscosity and the Bruggeman model(BM)for thermal conductivity are utilized to simulate nanoliquids with the NP aggregation aspect.The complicated nonlinear problem is treated numerically.It is found that the temperature of nanoliquid is enhanced due to the aggregation of NPs.The friction factor is more sensitive to the volume fraction of NPs than the thermal radiation and the mixed convection parameter.Furthermore,the heat transport rate is more sensitive to the effect of radiative heat compared with the NP volume fraction and mixed convection parameter.     

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

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.     

SOLUTIONS FOR CYLINDRICAL CAVITY IN SATURATED THERMOPOROELASTIC MEDIUM

Based on the thermodynamics of irreversible processes, the mass conservation equation and heat energy balance equation are established. The governing equations of thermal consolidation for homogeneous isotropic materials are presented, accounting for the coupling effects of the temperature, stress and displacement fields. The case of a saturated medium with a long cylindrical cavity subjected to a variable thermal loading and a variable hydrostatic pressure （or a variable radial water flux） with time is considered. The analytical solutions are derived in the Laplace transform space. Then, the time domain solutions are obtained by a numerical inversion scheme. The results of a typical example indicate that thermodynamically coupled effects have considerable influences on thermal responses.     

Thermal stresses in infinite circular cylinder subjected to rotation

The present investigation is concerned with the effect of rotation on an infinite circular cylinder subjected to certain boundary conditions.An analytical procedure for evaluation of thermal stresses,displacements,and temperature in rotating cylinder subjected to thermal load along the radius is presented.The dynamic thermal stresses in an infinite elastic cylinder of radius a due to a constant temperature applied to a variable portion of the curved surface while the rest of surface is maintained at zero temperature are discussed.Such situation can arise due to melting of insulating material deposited on the surface cylinder.A solution and numerical results are obtained for the stress components,displacement components,and temperature.The results obtained from the present semi-analytical method are in good agreement with those obtained by using the previously developed methods.     

MD simulation of a copper rod under thermal shock

In this paper, thermoelastic problem of onedimensional copper rod under thermal shock is simulated using molecular dynamics method by adopting embedded atom method potential. The rod is on axis x, the left outermost surface of which is traction free and the right outermost surface is fixed. Free boundary condition is imposed on the outermost surfaces in direction y and z. The left and right ends of the rod are subjected to hot and cold baths, respectively. Temperature, displacement and stress distributions are obtained along the rod at different moments, which are shown to be limited in the mobile region, indicating that the heat propagation speed is limited rather than infinite. This is consistent with the prediction given by generalized thermoelastic theory. From simulation results we find that the speed of heat conduction is the same as the speed of thermal stress wave. In the present paper, the simulations are conducted using the large-scale atomic/molecular massively parallel simulator and completed visualization software.     

TRANSIENT THERMAL RESPONSE IN THICK ORTHOTROPIC HOLLOW CYLINDERS WITH FINITE LENGTH： HIGH ORDER SHELL THEORY

The transient thermal response of a thick orthotropic hollow cylinder with finite length is studied by a high order shell theory. The radial and axial displacements are assumed to have quadratic and cubic variations through the thickness, respectively. It is important that the radial stress is approximated by a cubic expansion satisfying the boundary conditions at the inner and outer surfaces, and the corresponding strain should be least-squares compatible with the strain derived from the strain-displacement relation. The equations of motion are derived from the integration of the equilibrium equations of stresses, which are solved by precise integration method （PIM）. Numerical results are.obtained, and compared with FE simulations and dynamic thermo-elasticity solutions, which indicates that the high order shell theory is capable of predicting the transient thermal response of an orthotropic （or isotropic） thick hollow cylinder efficiently, and for the detonation tube of actual pulse detonation engines （PDE） heated continuously, the thermal stresses will become too large to be neglected, which are not like those in the one time experiments with very short time.     

The axisymmetric torsional contact problem of a functionally graded piezoelectric coated half-space

In this article, we study the axisymmetric tor-sional contact problem of a half-space coated with func-tionally graded piezoelectric material (FGPM) and subjected to a rigid circular punch. It is found that, along the thick-ness direction, the electromechanical properties of FGPMs change exponentially. We apply the Hankel integral trans-form technique and reduce the problem to a singular integral equation, and then numerically determine the unknown con-tact stress and electric displacement at the contact surface. The results show that the surface contact stress, surface azimuthal displacement, surface electric displacement, and inner electromechanical field are obviously dependent on the gradient index of the FGPM coating. It is found that we can adjust the gradient index of the FGPM coating to modify the distributions of the electric displacement and contact stress.     

TRANSIENT TORSIONAL WAVE IN FINITE HOLLOW CYLINDER WITH INITIAL AXIAL STRESS

An analytical solution is obtained for transient torsional vibration of a finite hollow cylinder with initial axial stress. The cylinder is subjected to dynamic shearing stress at the internal surface and is fixed at the external surface. The basic equations are presented and the solution is obtained by means of Fourier series expansion technique and the separation of variables method. The effects of the initial stress on the natural frequencies and transient torsional responses are presented and discussed.     

THERMOMECHANICAL RESPONSE OF SMA FIBER COMPOSITES WITH NON-UNIFORM TEMPERATURE DISTRIBUTION

A micromechanics method based on the High-Order-Theory developed by Aboudi et al. is used to predict the thermomechanical response of composites reinforced by shape memory alloy (SMA) fibers, and the non-uniform thermal distribution in composite arising from the process of heating or cooling is considered. The numerical development based on this model was coded to predict the thermomechanical response of shape memory alloy fiber/elastomer matrix composite subjected to thermal cycle loading. When the composite is heated, two heating ways, thermal gradients and heat source by passing an electric current through the SMA fibers are imposed on the composite respectively. Upon cooling, the first thermal boundary condition and the second thermal boundary condition are subjected to the composite respectively. A series of stress distributions and temperature distributions for different instants are calculated to reveal the interaction between the SMA material and matrix. It is useful to analyze and design the SMA actuator driven by heat source or the surface temperature.     

Thermoelastic analysis of functionally graded annulus with arbitrary gradient

A thermoelastic problem of a circular annulus made of functionally graded materials with an arbitrary gradient is investigated. Different from previous works, our analysis neither requires a special form of the gradient of material properties nor demands partitioning the entire structure into a multilayered homogeneous structure. Instead, we propose a new method for solving the thermoelastic problem of a functionally graded circular annulus by transforming it to a Fredholm integral equation. The distribution of thermal stresses and radial displacement can be obtained by solving the resulting equation. Illustrative examples are given to show the effects of varying gradients on the thermal stresses and radial displacement for given temperature changes at the inner and outer surfaces. The results indicate that the thermal stresses can be relaxed for specified gradients, which is beneficial to design an inhomogeneous annulus to maintain structural integrity.     

THERMAL CONSOLIDATION OF LAYERED POROUS HALF-SPACE TO VARIABLE THERMAL LOADING

An analytical method was derived for the thermal consolidation of layered, saturated porous half-space to variable thermal loading with time. In the coupled governing equations of linear thermoelastic media, the influences of thermo-osmosis effect and thermal filtration effect were introduced. Solutions in Laplace transform space were first obtained and then numerically inverted. The responses of a double-layered porous space subjected to exponential decaying thermal loading were studied. The influences of the differences between the properties of the two layers (e.g., the coefficient of thermal consolidation, elastic modulus) on thermal consolidation were discussed. The studies show that the coupling effects of displacement and stress fields on temperature field can be completely neglected, however, the thermo-osmosis effect has an obvious influence on thermal responses.