Assessment of Thermomechanical Couplings in Tribological Surface Transformations: Application to the Irreversible Near-Surface Solid-Solid Phase Transformations

Grégory Antoni Haute Ecole d’ingénierie et de Gestion du Canton de Vaud, Institut COMATEC, Route de Cheseaux 1, 1401 Yverdon-les-Bains, Switzerland Correspondence should be addressed to Grégory Antoni; antoni.gregory@yahoo.fr Received 27 November 2012; Accepted 27 January 2013 Academic Editors: V. E. Buchanan, Y. Desplanques, Z. Wen, H. Xing, and Y. Zhang Copyright © 2013 Grégory Antoni. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Tribological Surface Transformations (TSTs), which are irreversible near-surface solid-solid phase transformations, tend to occur on railroads frequented by heavy freight trains. The present study is proposed to assess the contribution of thermomechanical coupling processes to the emergence and development of TSTs near the surface of the rails.


Introduction
The irreversible quasi-surface solid-solid phase transformations known as Tribological Surface Transformations (TSTs), which occur near the surface of rails [1][2][3], are induced by contact fatigue processes [4,5], occurring after the passing of several trains on the rail tread. In order to account for the emergence and development of TSTs, it was suggested in a previous study [6] that TSTs result from repeated normal and stress loads, combined with highly localized thermal effects caused by friction between the wheels and rails. It is proposed here to test this assumption, namely that the thermomechanical coupling processes may generate TSTs near the surface of the rails.

A Thermodynamic Framework for Modelling Irreversible Near-Surface Solid-Solid Phase Transformations
In this section, the constitutive equations involved in the thermomechanical model for irreversible solid-solid phase transformations [6] are recalled: =̇t e +̇p c +̇p z , where is the absolute temperature, = 1/3Tr( )G + e is the small-strain tensor, is the elastic strain tensor, 2 ISRN Tribology = 1/3Tr( )G + e (where = (te, pc, pz)) are the thermoelastic, classical (visco-)plastic and TRIP-like strain tensors, respectively, G denotes the metric tensor, denotes the isotropic hardening variable associated with classical (visco-)plasticity, > 0 and > −(2/3) are the Lamé constants, > 0 is the thermal expansion coefficient, = / ∈ [0, 1] is the mass fraction of the daughter phase (where = martensite is partial mass density of the daughter phase, = + is total mass density, and = ferrite is partial mass density of the parent phase (i.e., the ferrite phase)), > 0 is a material parameter characterizing the change in the density occurring during phase transformations, is the solid-solid phase transformation temperature when the pressure is zero, > 0 is a material parameter characterizing a "sensitive pressure" level, ℎ > 0 is a material parameter characterizing the linear isotropic hardening associated with classical plasticity, ≥ 0 is a material parameter associated with the latent heat of the phase transformation (see [7]), is the Cauchy stress tensor, − G = 1/3Tr( )G and s are the spherical and deviatoric parts of the Cauchy stress tensor, respectively, eq = (3/2 s : s) 1/2 is the Von Mises equivalent stress, is the classical yield strength, and are the characteristic time of the viscous effects associated with TRIP-like and the classical plasticity processes, respectively, ⟨⋅⟩ denotes the Macaulay brackets (⟨ ⟩ = when ≥ 0 and ⟨ ⟩ = 0 when < 0), and (⋅) denotes the Heaviside step function ( ( ) = 1 when ≥ 0 and ( ) = 0 when < 0).
Based on the assumption that TSTs are induced by "strong" thermomechanical coupling-the combined effects of the temperature and mechanical processes-the heat equation can be written as follows (see [7]), takinġto denote the increase in the temperature: where > 0 is the specific heat capacity, is the initial temperature, is the thermal conductivity, and Δ denotes the Laplacian. The mass density is roughly an equal value to the initial mass density, ≃ . The parameters, , , , , , and are not assumed to depend on the temperature and are taken to be identical in both phases.
Note that the thermodynamic consistency of this thermomechanical model was checked in our previous study [6].

Heat Equation:
Source Terms. In this Section, it is proposed to assess the effects of the "source" terms on the evolution of the temperature on the emergence and development of TSTs. In our previous study [6], the source terms corresponding to the thermoelasticity, the classical (visco-)plasticity, and the irreversible solid-solid phase transformations were omitted in the first approximation, and in the second approximation, we dealt only with the steady case of the heat equation, that is, Δ .
However, the orders of magnitude of these terms can be assessed, a posteriori, and compared. In the incremental form, the full heat equation (2) can written as follows: where (⋅) +1 (resp. (⋅) ) denote the quantities at time +1 (resp. ) and Δ = +1 − is the time step.
Although these source terms (Φ te , Φ pc and Φ pz ) were obtained from a temperature field (see [6]) resulting from the resolution of the heat equation in the steady case, the results obtained in the previous study [6] can be used to assess, in a first approach, the heat sources at any time, as defined by the term on the right-hand side of (3). Note that the physical interpretations associated with each parameter of (3), including their chosen set value, are discussed in [7]. Figure 1 shows the heat field sources at time = 10 s under the thermomechanical loading conditions studied in [6]. Note that the heat sources thus obtained are much greater than the power of the internal stresses ( Figure 2). The heat source caused by the thermoelastic processes alone is plotted in Figure 3. In fact, in the volume measuring 25 mm 3 (10 mm × 2.5 mm × 1 mm) in the immediate vicinity of the surface where the thermomechanical loading is applied, the power dissipation of the heat sources caused by the thermoelastic and classical (visco-)plastic processes and the irreversible solid-solid phase transformation coupling is approximately 425.10 −3 W, whereas the power of the internal stresses is only of 145.10 −3 W. On the other hand, the power dissipation of the heat source due the thermoelastic effects alone is approximately 402.10 −3 W. It can be seen here that under the boundary conditions defined in [6], the thermoelastic effects are nonnegligible. It is now proposed to assess these thermoelastic coupling effects in the case of a unidimensionnal problem (Section 3.2).

Thermoelastic Coupling: 1-D Problem.
In order to estimate the thermoelastic coupling effects, it is proposed to take a 1-D thermoelastic problem, that of a bar in steel (occupying a domain Ω∈[0, ]⊂R) exposed to thermomechanical loads (Figure 4), assuming the presence of small perturbations.

ISRN Tribology
The Helmholtz free energy potential per unit mass of the material considered (see [8]) can be written as follows: where and denote the absolute temperature and the smallstrain, respectively, is the initial mass density, is the Young's modulus, > 0 is the specific heat capacity, is the initial temperature, is the thermal expansion coefficient, and is the initial Helmholtz free energy potential per unit mass.
H1: the mass density in the actual ( ) and initial ( ) configurations is approximately equal, that is, ≃ .
Assuming H1, the state functions associated with both the temperature and the strain variable are: where and denote the specific entropy per unit mass and the Cauchy stress. The local expression of First Principle of Thermodynamics (FPT) is: (7) where = + TS denotes the internal energy potential per unit mass and q is the heat flux vector. Note that the remote heat production rate will be omitted from now on, that is, = 0.
H2: Fourier's law is adopted: where q denotes the heat flux vector. Combining H1, H2, and (5), the FPT can be written as follows (also called "Heat Equation"): H3: the absolute temperature in the actual ( ) and initial ( ) configurations are fairly similar, that is, / ≅ 1.
Using (6) and based on H2 and H3, (9) can be written as follows:̇− H4: the term " " is very large in comparison with " 2 " (because ⋘ 1), that is, ⋙ 2 . The Heat Equation (10) obtained under the assumption H4, neglecting the convective terms of the particular derivative and taking the lateral thermal leaks in the bar into account, can be rewritten in the following form: where denotes the characteristic time of the lateral thermal leaks in the bar. Assuming that both inertial and gravity effects are nonpresent, the local expression of quasi-static equilibrium can be reduced to: Based on (12), the stress field ( , ) is uniform throughout the bar, that is, for all ∈ Ω, ( , ) = 0 ( ). The boundary conditions are summarized in Figures 4 and 5, and the constitutive thermomechanical parameters used are: = 78 × 10 −7 kg ⋅ mm −3 ; = 12 × 10 −6 K −1 ; The results obtained in this section show that the source term due to the thermoelasticity is responsible for a variably marked increase in the temperature increase near the end of the bar and throughout the bar, but that the thermomechanical loading conditions considered here (see Figure 5) have greater effects in the immediate vicinity of the end of the bar.
In line with what was observed in our previous study [6], these terms can therefore be reasonably neglected in the first approximation.

Conclusion
This paper deals with the contribution of thermomechanical coupling processes to the emergence and development of irreversible near-surface solid-solid phase transformations. Based on the results obtained here, the thermoelastic effects can be said to predominate over the effects of the other thermomechanical heat source resulting from classical (visco-)plasticity and irreversible solid-solid transformations. Although these source terms cause the temperature to  increase at the surface and in the depth of the structure, it emerges clearly from this study that the thermomechanical loads applied constitute the main driving force responsible for generating TSTs.