A Fast Approximate Method for Heat Conduction in an Inhomogeneous Half-Space Subjected to Frictional Heating

2018 ◽  
Vol 140 (4) ◽  
Author(s):  
Xiujiang Shi ◽  
Liqin Wang ◽  
Qinghua Zhou ◽  
Qian Wang

This paper reports a new three-dimensional model for heat conduction in a half-space containing inhomogeneities, applicable to frictional heat transfer, together with a novel combined algorithm of the equivalent inclusion method (EIM) and the imaging inclusion approach for building this model. The influence coefficients (ICs) for temperature and heat flux are obtained via converting the frequency response function (FRF) and integrating Green's function. The model solution is based on the discrete convolution and fast Fourier transform (DC-FFT) algorithm using the ICs, convenient for solving problems involving multiple elliptical inhomogeneities with arbitrary orientations. A group of parametric studies are conducted for understanding the thermal fields in the inhomogeneous half-space due to surface frictional heating, influenced by the properties of the inhomogeneity, its depth, and orientation.

2013 ◽  
Vol 135 (4) ◽  
Author(s):  
Zhanjiang Wang ◽  
Xiaoqing Jin ◽  
Leon M. Keer ◽  
Qian Wang

Contacts involving partial slip are commonly found at the interfaces formed by mechanical components. However, most theoretical investigations of partial slip are limited to homogeneous materials. This work proposes a novel and fast method for partial-slip contact involving a material with an inhomogeneity based on the equivalent inclusion method, where the inhomogeneity is replaced by an inclusion with properly chosen eigenstrains. The stress and displacement fields due to eigenstrains are formulated based on the half-space inclusion solutions recently derived by the authors and solved with a three-dimensional fast Fourier transform algorithm. The effectiveness and accuracy of the proposed method is demonstrated by comparing its solutions with those from the finite element method. The partial slip contact between an elastic ball and an elastic half space containing a cuboidal inhomogeneity is further investigated. A number of in-depth parametric studies are performed for the cuboidal inhomogeneity with different sizes and at different locations. The results reveal that the contact behavior of the inhomogeneous material is more strongly influenced by the inhomogeneity when it is closer to the contact center and when its size is larger.


Author(s):  
Zhanjiang Wang ◽  
Xiaoqing Jin ◽  
Shuangbiao Liu ◽  
Leon M. Keer ◽  
Jian Cao ◽  
...  

This paper presents a new method of contact plasticity analysis based on Galerkin vectors to solve the eigenstresses due to eigenstrain. The influence coefficients relating eigenstrains to eigenstresses thus can be divided into four terms the one due to the eigenstrains in the full space, others due to the mirrored eigenstrains in the mirror half space. Each term can be solved fast and efficient by using the three-dimensional discrete convolution and fast Fourier transform (DC-FFT) or the three-dimensional discrete correlation and fast Fourier transform (DCR-FFT). The new method is used to analyze the contact plastic residual stresses in half space.


2002 ◽  
Vol 125 (1) ◽  
pp. 33-43 ◽  
Author(s):  
Shuangbiao Liu ◽  
Qian Wang

Computing the thermoelastic stress field of a material subjected to frictional heating is essential for component failure prevention and life prediction. However, the analysis for three-dimensional thermoelastic stress field for tribological problems is not well developed. Furthermore, the pressure distribution due to rough surface contact is irregular; hence the frictional heating can hardly be described by an analytical expression. This paper presents a novel set of frequency-domain expressions (frequency response functions) of the thermoelastic stress field of a uniformly moving three-dimensional elastic half-space subjected to arbitrary transient frictional heating, where the velocity of the half-space, its magnitude and direction, can be an arbitrary function of time. General formulas are expressed in the form of time integrals, and important expressions for constant velocities are given for the transient-instantaneous, transient-continuous, and steady-state cases. The thermoelastic stress field inside a translating half-space with constant velocities are illustrated and discussed by using the discrete convolution and fast Fourier transform method when a parabolic type or an irregularly distributed heat source is applied.


2008 ◽  
Vol 130 (4) ◽  
Author(s):  
W. Wayne Chen ◽  
Q. Jane Wang

A thermomechanical analysis of elasto-plastic bodies is a necessary step toward the understanding of tribological behaviors of machine components subjected to both mechanical loading and frictional heating. A three-dimensional thermoelastoplastic contact model for counterformal bodies has been developed, which takes into account steady state heat flux, temperature-dependent strain hardening behavior, and interaction of mechanical and thermal loads. The fast Fourier transform and conjugate gradient method are the underlying numerical algorithms used in this model. Sliding of a half-space over a stationary sphere is simulated with this model. The friction-induced heat is partitioned into two bodies based on surface temperature distributions. In the simulation, the sphere is considered to be fully thermoelastoplastic, while the half-space is treated to be thermoelastic. Simulation results include surface pressure, temperature rise, and subsurface stress and plastic strain fields. The paper also studies the influences of sliding speed and thermal softening on contact behaviors for sliding speed ranging three orders of magnitude.


2016 ◽  
Vol 08 (02) ◽  
pp. 1650014 ◽  
Author(s):  
Kun Zhou ◽  
Qingbing Dong

This paper develops a three-dimensional (3D) model for a heterogeneous half-space with inclusions distributed periodically beneath its surface subject to elastohydrodynamic lubrication (EHL) line-contact applied by a cylindrical loading body. The model takes into account the interactions between the loading body, the fluid lubricant and the heterogeneous half-space. In the absence of subsurface inclusions, the surface contact pressure distribution, the half-space surface deformation and the lubricant film thickness profile are obtained through solving a unified Reynolds equation system. The inclusions are homogenized according to Eshelby’s equivalent inclusion method (EIM) with unknown eigenstrains to be determined. The disturbed half-space surface deformations induced by the subsurface inclusions or eigenstrains are iteratively introduced into the lubricant film thickness until the surface deformation finally converges. Both time-independent smooth surface contact and time-dependent rough surface contact are considered for the lubricated contact problem.


2000 ◽  
Vol 123 (1) ◽  
pp. 17-26 ◽  
Author(s):  
Shuangbiao Liu ◽  
Qian Wang

A necessary step in understanding failure problems of tribological elements is to investigate the contact performance of rough surfaces subjected to frictional heating. It is essential that the interfacial variables are obtained through solving the interactive thermomechanical contact problem. This paper studies the three dimensional thermomechanical contact of non-conforming rough surfaces, the model of which includes the normal surface displacements caused by the contact pressure, frictional shear, and frictional heating. Influence coefficients and frequency response functions for elastic and thermoelastic displacements, as well as those for temperature rises, are investigated for model construction. In order to develop an accurate and efficient solver, the numerical algorithms with the discrete convolution and fast Fourier transform techniques and the single-loop conjugated gradient method are used. The model modules are numerically verified and the thermomechanical performance of the rough surfaces in a point contact is studied.


Author(s):  
B. S. Yilbas ◽  
M Sami

Lasers are widely used as a machine tool in the metal industry. One of the important areas of laser application is surface treatment of engineering metals. To improve the process parameters in the laser heating process, an exploration of the heating mechanism is fruitful. The present study is carried out to develop a three-dimensional model for a laser pulsed heating process using the electron kinetic theory approach. The heating model introduced relies on successive electronphonon collisions; therefore, it is this process that describes the heat conduction mechanism. This study is limited to heat conduction only. Consequently, the phase change process is not taken into account. To validate the theoretical predictions, an experiment is conducted to measure the surface temperature using an optical method. Moreover, a one-dimensional heating model developed previously is also considered and the predictions of three- and one-dimensional heating models as well as experimental results are compared. It is found that the three-dimensional model gives lower surface temperatures compared with the one-dimensional model considered. However, experimental results agree well with the results obtained from the three-dimensional model. In addition, an equilibrium time is introduced. In that case, energy gain of electrons via incident beam absorption balances the energy losses due to conduction through successive electron-phonon collisions.


2000 ◽  
Vol 123 (2) ◽  
pp. 330-342 ◽  
Author(s):  
Wei Peng ◽  
Bharat Bhushan

A new numerical model for the three-dimensional contact analysis of a layered elastic–perfectly plastic half space with another rough surface is presented. The model is based on a variational principle in which the real area of contact and contact pressure distribution are those which minimize the total complementary potential energy. A quasi-Newton method is used to find the minimum. The influence coefficients matrix is determined using the Papkovich–Neuber potentials with fast Fourier transformation. The model is extended to elastic–perfectly plastic contacts in dry and wet conditions. Contact analyses have been performed to predict contact statistics of layered elastic/plastic solids with rough surfaces using this model. The effects of the stiffness of the layer and the substrate, layer thickness, as well as normal load are studied. Optimum layer parameters are identified to provide low friction/stiction and wear.


Sign in / Sign up

Export Citation Format

Share Document