Thermomechanical Analysis of Elastoplastic Bodies in a Sliding Spherical Contact and the Effects of Sliding Speed, Heat Partition, and Thermal Softening

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.

Author(s):  
W. Wayne Chen ◽  
Q. Jane Wang

Frictional heating leads to the temperature rise, thermal expansion, and the thermo-elasto-plastic deformation, which may be responsible for the failure of components under contact and relative motion. This paper reports the development of a three-dimensional thermo-elasto-plastic contact model of counterformal bodies, which account for the transient heat transfer, temperature-dependent strain hardening behavior of materials, and realistic heat partition between surfaces. An extensive study on the contact of a sliding half-space over a stationary ball is conducted using this model.


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.


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.


2000 ◽  
Vol 123 (3) ◽  
pp. 595-602 ◽  
Author(s):  
Geng Liu ◽  
Qian Wang ◽  
Shuangbiao Liu

The rough surface contact in a tribological process involves frictional heating and thermoelastic deformations. A three-dimensional thermal-mechanical asperity contact model has been developed, which takes into account steady-state heat transfer, asperity distortion due to thermal and elastic deformations, and material yield. The finite-element method (FEM), fast Fourier transform (FFT), and conjugate gradient method (CGM) are employed as the solution methods. The model is used to analyze the thermal-mechanical contact of typical rough surfaces and investigate the importance of thermal effects on the contact performance of surface asperities.


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.


1971 ◽  
Vol 38 (4) ◽  
pp. 899-905 ◽  
Author(s):  
L. B. Freund

Three-dimensional wave propagation in an elastic half space is considered. The half space is traction free on half its boundary, while the remaining part of the boundary is free of shear traction and is constrained against normal displacement by a smooth, rigid barrier. A time-harmonic surface wave, traveling on the traction free part of the surface, is obliquely incident on the edge of the barrier. The amplitude and the phase of the resulting reflected surface wave are determined by means of Laplace transform methods and the Wiener-Hopf technique. Wave propagation in an elastic half space in contact with two rigid, smooth barriers is then considered. The barriers are arranged so that a strip on the surface of uniform width is traction free, which forms a wave guide for surface waves. Results of the surface wave reflection problem are then used to geometrically construct dispersion relations for the propagation of unattenuated guided surface waves in the guiding structure. The rate of decay of body wave disturbances, localized near the edges of the guide, is discussed.


Author(s):  
L. Salles ◽  
M. Vahdati

The aim of this paper is to study the effects of mistuning on fan flutter and to compare the prediction of two numerical models of different fidelity. The high fidelity model used here is a three-dimensional, whole assembly, time-accurate, viscous, finite-volume compressible flow solver. The Code used for this purpose is AU3D, written in Imperial College and validated for flutter computations over many years. To the best knowledge of authors, this is the first time such computations have been attempted. This is due to the fact that, such non-linear aeroelastic computations with mistuning require large amount of CPU time and cannot be performed routinely and consequently, faster (low fidelity) models are required for this task. Therefore, the second model used here is the aeroelastic fundamental mistuning model (FMM) and it based on an eigenvalue analysis of the linearized modal aeroelastic system with the aerodynamic matrix calculated from the aerodynamic influence coefficients. The influence coefficients required for this algorithm are obtained from the time domain non-linear Code by shaking one blade in the datum (tuned) frequency and mode. Once the influence coefficients have been obtained, the computations of aero damping require minimal amount of CPU time and many different mistuning patterns can be studied. The objectives of this work are to: 1. Compare the results between the two models and establish the capabilities/limitations of aeroelastic FMM, 2. Check if the introduction of mistuning would bring the experimental and computed flutter boundaries closer, 3. Establish a relationship between mistuning and damping. A rig wide-chord fan blade, typical of modern civil designs, was used as the benchmark geometry for this study. All the flutter analyses carried out in this paper are with frequency mistuning, but the possible consequences of mistuned mode shapes are briefly discussed at the end of this paper. Only the first family of modes (1F, first flap) is considered in this work. For the frequency mistuning analysis, the 1F frequency is varied around the annulus but the 1F mode shapes remain the same for all the blades. For the mode shape mistuning computations, an FE analysis of the whole assembly different mass blades is performed. The results of this work clearly show the importance of mistuning on flutter. It also demonstrates that when using rig test data for aeroelastic validation of CFD codes, the amount mistuning present must be known. Finally, it should be noted that the aim of this paper is the study of mistuning and not steady/unsteady validation of a CFD code and therefore minimal aerodynamic data are presented.


2017 ◽  
Vol 84 (11) ◽  
Author(s):  
Yilan Huang ◽  
Guozhan Xia ◽  
Weiqiu Chen ◽  
Xiangyu Li

Exact solutions to the three-dimensional (3D) contact problem of a rigid flat-ended circular cylindrical indenter punching onto a transversely isotropic thermoporoelastic half-space are presented. The couplings among the elastic, hydrostatic, and thermal fields are considered, and two different sets of boundary conditions are formulated for two different cases. We use a concise general solution to represent all the field variables in terms of potential functions and transform the original problem to the one that is mathematically expressed by integral (or integro-differential) equations. The potential theory method is extended and applied to exactly solve these integral equations. As a consequence, all the physical quantities of the coupling fields are derived analytically. To validate the analytical solutions, we also simulate the contact behavior by using the finite element method (FEM). An excellent agreement between the analytical predictions and the numerical simulations is obtained. Further attention is also paid to the discussion on the obtained results. The present solutions can be used as a theoretical reference when practically applying microscale image formation techniques such as thermal scanning probe microscopy (SPM) and electrochemical strain microscopy (ESM).


2001 ◽  
Vol 27 (8) ◽  
pp. 513-520
Author(s):  
Ugur Tanriver ◽  
Aravinda Kar

This note is concerned with the three-dimensional quasi-steady-state heat conduction equation subject to certain boundary conditions in the wholex′y′-plane and finite inz′-direction. This type of boundary value problem arises in laser welding process. The solution to this problem can be represented by an integral using Fourier analysis. This integral is approximated to obtain a simple analytic expression for the temperature distribution.


Sign in / Sign up

Export Citation Format

Share Document