Journal of Materials

Volume 2013, Article ID 892050, 7 pages

http://dx.doi.org/10.1155/2013/892050

## Effects of ThermoMechanical Coupling in Tribological Surface Transformations: A One-Dimensional Modelling Including Irreversible Solid-Solid Phase Transformations and Classical Plasticity

Haute Ecole d'Ingénierie et de Gestion du Canton de Vaud, Institut COMATEC, Route de Cheseaux 1, 1401 Yverdon-les-Bains, Switzerland

Received 20 November 2012; Revised 1 February 2013; Accepted 1 February 2013

Academic Editor: Jacques Huot

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.

#### Abstract

Under compressive loads combined with friction, some materials undergo Tribological Surface Transformations (TSTs) on the surface of the loaded parts and in the immediately vicinity, which in the case of metals, are known as irreversible solid-solid phase transformations. During the solid-solid phase transformations occurring under mechanical loads, TRansformation Induced Plasticity (TRIP) processes are generated at much lower stress levels than those associated with the yield strength of the material in classical plasticity. In order to assess the effects of thermomechanical coupling in these TSTs, a one-dimensional modelling based on irreversible solid-solid phase transformations and classical plasticity is presented and discussed.

#### 1. Introduction

The TSTs which have been affecting some rails of the French railroad for the last twenty years consist of irreversible solid-solid phase transformations resulting in the development of a “White Etching Layer” (WELs) on the rail tread and which grows in depths from several nanometers to more than 100 *μ*m [1]. After the passing of some trains, the first “white” phase layers appear on the surface of the rails (which initially have a ferritic-pearlitic structure) before penetrating more deeply, and resulting in a “quasi-surfacic” martensitic transformation, that is, the material composing the rails is transformed into martensite [2]. The initiation and development of TSTs are induced by the combined thermomechanical loads in the wheel/rail contact area, where the normal and tangential stresses are combined with a temperature increase due to the friction occuring in the contact, that is, TSTs result from the thermomechanical coupling processes. During the steel phase transformations which occur under mechanical loads, an anomalous plastic flow known as TRansformation Induced Plasticity (TRIP) is generally superimposed on the classical plastic processes at stress levels which are much lower than the initial strength of the softest phase. Greenwood and Johnson [3] has explained TRIP process by suggesting that the microplasticity generated by the incompatibility between the volumes involved in the phase transformation may be oriented by the external stress applied. Many models have been presented to account for TRIP phenomena [4–8]. In order to assess the thermomechanical coupling associated with TSTs, a one-dimensional version of the model presented [9], which was developed in a thermodynamically consistent framework, is presented in Section 2. In this model, the irreversible phase transformations are taken to result from “strong” thermomechanical coupling, that is, from the combined effects of the change in temperature and the mechanical dissipation. A numerical example is illustrated in Section 3 in the case of a bar subjected to a thermomechanical loading, which is liable to undergo an irreversible solid-solid phase transformation, possibly with the classical plasticity.

#### 2. A 1D Phenomenological Modelling Approach Applied to Tribological Surface Transformations (TSTs): Thermomechanical Coupling

Based on the previous TRIP studies, which are extended here to the thermomechanical coupling associated with TSTs [9], the 1D phenomenological model presented in this Section is developed as follows:(i)the absolute temperature is explicitly considered and taken to be a state variable;(ii)the mass fraction of the daughter phase (i.e., the “white” martensite phase) is:
where is partial mass density of the daughter phase, is total mass density, and is partial mass density of the parent phase (i.e., the ferrite phase);(iii)the total infinitesimal strain can be written as follows:
where , , denote the elastic, thermal, classical (visco-) plasticity strains, respectively, is the transformation strain associated with the volume variation, is TRIP-like strain, is initial temperature, is the thermal expansion coefficient, and is a material parameter characterizing the change in the density occurring during phase transformation. Note that the solid-solid phase transformation under consideration here induces a volume change, or more specifically, a decrease in the volume, that is, , since the density mass of ferrite (the initial material of the rail before phase transformation) kg·m^{−3}, and that of martensite (material after phase transformation) kg·m^{−3}.

The material in question is initially nontransformed (, and ) and nonplasticized (, and where is the isotropic hardening variable associated with classical plasticity) at the temperature , and the Helmholtz free energy potential per unit mass is written as: where is the specific heat capacity, is Young’s modulus, is the solid-solid phase transformation temperature, is a material parameter characterizing the linear isotropic hardening associated with classical plasticity, is a material parameter associated with the latent heat of the phase transformation, is the indicator function ( when and otherwise) and and are the initial density and the initial Helmholtz free energy of the material per unit mass, respectively. In a first approximation, the thermoelastic parameters, , , and , are not assumed to depend on the temperature and are taken to be identical in the both phases. Equation (3) is strongly inspired by SMAs model with a 1-D-Helmholtz free energy potential developed by [10] where the thermomechanical coupling process considered here is extended to include an irreversible solid-solid phase transformation. Note that in (3), the first two terms correspond to the thermal leaks (irreversible part), the third one corresponds to the free energy part at constant temperature which may be immediately recoverable by an elastic unloading (reversible part), and the next two terms correspond to the hardening associated with the classical plasticity (isotropic hardening) and the latent heat of the phase transformation, respectively.

Neglecting the viscoelastic effects and noting that the specific entropy is a state function, the state equations based on the above representation of the Helmholtz energy (3) can be written as: where is the Cauchy stress tensor, and denote the “thermodynamic forces” associated with the variables .

Using (3)-(4) and assuming that , the Cauchy stress tensor is: The local expression for the First Principle of Thermodynamics can be written (see e.g., [11]) as: where denotes the internal energy potential per unit mass, is the heat flux vector; the remote heat production rate has been omitted here.

The Second Principle of Thermodynamics must be checked, whatever the local thermomechanical state of the material, whatever the evolution locally undergone by the material, and whatever the temperature gradient acting locally on the material (see [10]) as follows: Adopting Fourier’s law (, where denotes the thermal conductivity), and combining (3)-(4) and (7), the Second Principle of Thermodynamics must satisfy: The laws of evolution of the internal state variables are.

(i) In the irreversible solid-solid phase transformation (TST) process: with and where denotes the yield strength of the softest phase (parent phase) in the TST when , and are a material parameter associated with TST’s yield strength and the hardening process, respectively, is the absolute value function ( when and when ), denotes the Macaulay brackets ( when and when ), is the characteristic time of the viscous effects associated with the TST process, and denotes the Heaviside step function ( when and when ).

It must be stressed that: (i) TSTs (or WELs), which correspond to the ferrite-pearlite martensite transformations, do not fit the framework of conventional metallurgical transformations. According to these metallurgical theories, the white layer (martensite) may result from the martensitic transformation (high velocity cooling) after a sufficiently large temperature increase to produce an austenitic transformation. Dang Van et al. [12] have established that in the case of a steadily moving train, the temperature will increase at the surface of the rails by or so at most. Therefore, this temperature increase alone, which is caused by friction involved in the contact, is obviously too small to generate solid-solid phase transformations. The assumption adopted here is therefore that the thermomechanical coupling, where the thermal stresses superimposed on the purely mechanical (isothermal) stresses, can be taken to generate TSTs. (ii) Equation (11) means that the irreversible solid-solid phase transformation is activated by the stress which acts as “driving force,” and accelerated by temperature (see [9] for more details). TST occurs at stress levels which trigger of dislocation motions due to the strain incompatibilities between the ferrite-pearlite and martensite phases. The triggering of these dislocation motions is only facilitated by the increase in the temperature. In agreement with (i), (10) reflects the fact that the validity field of the law of evolution governing the TST process requires the presence of both the stress and the temperature.

(ii) In the classical plasticity process: with and where denotes the classical yield strength of the material when , is the yield strength of the hardest phase (daughter phase) and is the characteristic time of the viscous effects associated with the classical plasticity process. Note that when (resp., ), the yield strength is only composed to one phase, that is, the softest phase (resp., the hardest phase) and it consists of two phases when .

Combining (3), (4), and (6) and Fourier’s law, the heat equation is given by: where denotes the Laplacian operator. The quantities on the right-hand side of (15) are the “source terms” corresponding to the thermoelasticity, the classical plasticity, and the irreversible solid-solid phase transformation, respectively.

It is worth noting the thermodynamic consistency of this model [9] is checked since the laws of evolution of the internal state variables verify systematically (8).

#### 3. Numerical Example

*Preliminary Remark*. The case of the structure studied in this section is not representative of the complex processes that could be expected between wheel and rail, such as nonlinear contact and friction. This one-dimensional problem was simply intended to show that the thermomechanical coupling can generate the initiation and development of TSTs in the immediate vicinity where a sufficient thermomechanical load level is applied on the structure. Therefore, it leads on naturally to first the step to the setting up of a TSTs benchmark.

This example will be addressed in order to assess the thermomechanical coupling occurring in the TSTs using the approach presented in Section 2, subject to the assumptions stated therein. Let us consider a solid occupying a one-dimensional space , measuring of length, that is, a solid bar where the material considered here is made of steel. On the left end of the bar, a thermomechanical loading , which progresses linearly in conjunction with time , is applied whereas both zero displacement () and a set temperature () are imposed constant with time on the right end of the bar, that is, , and , (see Figure 1). The mechanical component of the loading is a compression which decreases from zero to in and its thermal component, and it increases from to in the same time. Bearing in mind that the one-dimensional problem studied cannot take into account both the non-uniformity of normal and shear stresses in contact area including the temperature increase caused by friction, the maximum values chosen to describe the thermomechanical loading would seem to be realistic enough at a contact point under moderate load for pressure and temperature in increase. Furthermore, in a first approximation, the lateral surfaces of the bar is assumed to be adiabatic thereby both the temperature gradient and heat flux are one-dimensional (here, -direction), that is, , and , (where is a unit vector).

Assuming that the convective terms of particular derivative are zero, (15) can be written in the following form: where and denote the first and second partial derivatives with respect to and (temporal and spatial quantities).

Neglecting the effects of both the inertia and gravity, the local expression of quasi-static equilibrium reduces to:
In line with (17), the stress field is uniform, that is, . The geometry involved and the boundary conditions are summarized in Figure 1. The constitutive thermomechanical parameters used are: ; ; kg·mm^{−3}; ; ; ; ; W·mm^{−1}·K^{−1}; ; ; ; ; ; ; bmbox?>J·mm^{−3}·K^{−1}. Note that , , , and are classical thermoelastic parameters associated with a ferritic or martensitic steel, which are considered here constants in temperature range studied. As pointed out in Section 2, the material parameter characterizing the change in the density occurring during phase transformation (where density mass of ferrite kg·m^{−3} and martensite kg·m^{−3}). Therefore, based on these experimental data, the transformation ferrite martensite is accompanied by a volume contraction (for a constant mass given) to be taken into account, that is,
On the other hand, the equation of conservation of mass (under small perturbations assumption) can be written as:
Combining (5), (9), and (12) and considering both a return to thermal and mechanical equilibrium (i.e., and ), we can write that:
Assuming both that and and using (9) (here for a mechanical evolution and then a return to the mechanical state ), (20) can be rewritten:
In the case of full phase transformation (), (21) can be reduced to:
Combining (19) and (22), the expression of (18) can be put in the following form:

From (18), we can deduce that . Note also that the characteristic times of the viscous effects associated with the TRIP-like process and the classical plasticity are rather low which means that the viscous effects are also small (when and , these are nonexistent). The reader is invited to refer to [9] the physical interpretations associated with other parameters including their chosen set value.

The results of the numerical simulations are presented in Figures 2, 3, 4, and 5. They show that due to both the nonstationary heat equation with source terms (16) and quasi-static equilibrium (17) combined with the thermomechanical loading (see Figure 1), a uniform -field (see Figure 2(a)) and a non-uniform temperature distribution strongly localized (near where the load is applied, see Figure 2(b)) are obtained throughout of the bar. These source terms present in the heat equation will provide a maximum temperature increase of a few degrees (around 15 K) which occurs near the left end of the bar (at time ), whereas no such increase occurred elsewhere. The increase of about 10 K, which can be observed in the close vicinity of the point where the thermomechanical loading is applied (in the same time), is due only to the thermoelastic coupling. These results indicate that a “new” phase () develops in the immediate vicinity of the point on the bar where the thermomechanical loading is applied: from to about at time , at time , and in depth at time (see Figure 3(a)). Figure 3(b) brings out the distribution of the TRIP-like strain in the depth of the bar at time . Note that at time , , . A full phase transformation () develops up to about at time (see Figure 4(a)), where there is only one phase, the daughter phase (martensite). At depths of , the partial transformation () persists, that is, the material is a mixture between the parent phase (ferrite-pearlite) and the daughter phase (martensite). Distribution of transformation strain in total bar length is also depicted in the same time. Figure 4(b) shows the distribution of the TRIP-like and classical plasticity strain in the depth of the bar at time . Near the left end of the bar, the TRIP-like strain predominates, reaching with a extremum value of at time , whereas the strain associated with classical plasticity amounted to at the same time. It can be observed that at the point (Figure 5): (i) on the loading path (where ), the response of the material is purely thermoelastic (i.e., ). The thermal effects gives a non-linear stress-strain curve; (ii) on the loading path (where ), thermo-elasto-viscoplastic behavior occurs due to the solid-solid phase transformation and classical plasticity (i.e., , , , ); (iii) on the loading path (where ), a return to the thermoelastic response of the material (i.e., ) occurs.

The results presented in Section 3 are in line with the experimental available data on TSTs [1], since a full layer transformed is strongly localized and achieved a realistic TST thickness, that is, at a maximum depth of around 100 *μ*m (see Section 1). The thermomechanical model proposed here can describe both the initiation and development of TSTs at the surface and in the immediate vicinity where the thermomechanical loading was applied. Therefore, this model will be examined on the wheel/rail contact where both the nonlinear contact and friction can take place.

#### 4. Conclusion

The thermodynamically consistent model presented in this study can be used to simulate the initiation and development of Tribological Surface Transformations (TSTs). The results obtained here show that the thermomechanical coupling can generate a thin layer of fully irreversible solid-solid phase transformation in the immediate vicinity of the point on the bar where the thermomechanical loading is applied. These results also highlight that the layer thickness transformed is in line with the experimental available data on TSTs since a depth of around one hundred micrometers has been reached. In the future studies, it is proposed to test the present thermomechanical model in more realistic numerical simulations involving issues relatied to wheel/rail contact problems.

#### Acknowledgment

The author is indebted to Professor Thierry Désoyer and Professor Frédéric Lebon for the valuable discussions.

#### References

- Österle, Rooch, Pyzalla, and Wang, “Investigation of white etching layers on rails by optical microscopy, electron microscopy, X-ray and synchroton X-ray diffraction,”
*Materials Science and Engineering A*, vol. 303, no. 1-2, pp. 150–157, 2001. View at Publisher · View at Google Scholar · View at Scopus - G. Baumann, H. J. Fecht, and S. Liebelt, “Formation of white-etching layers on rail treads,”
*Wear*, vol. 191, no. 1-2, pp. 133–140, 1996. View at Google Scholar · View at Scopus - G. W. Greenwood and R. H. Johnson, “The deformation of metals under small stresses during phase transformations,”
*Proceedings of the Royal Society of London A*, vol. 293, pp. 403–422, 1965. View at Google Scholar - Leblond, Devaux, and Devaux, “Mathematical modelling of transformation plasticity in steels I: case of ideal-plastic phases,”
*International Journal of Plasticity*, vol. 5, no. 6, pp. 551–572, 1989. View at Google Scholar · View at Scopus - J. B. Leblond, J. Devaux, and J. C. Devaux, “Mathematical modelling of transformation plasticity in steels II: coupling with strain hardening phenomena,”
*International Journal of Plasticity*, vol. 5, no. 6, pp. 573–591, 1989. View at Google Scholar · View at Scopus - J.-C. Videau, G. Cailletaud, and A. Pineau, “Experimental study of the transformation induced plasticity in a Cr-Ni-Mo-Al-Ti steel,”
*Journal de Physique*, vol. 6, Supplément au Journal de Physique III, no. 1, pp. 465–474, 1996. View at Google Scholar - L. Taleb and F. Sidoroff, “A micromechanical modeling of the Greenwood-Johnson mechanism in transformation induced plasticity,”
*International Journal of Plasticity*, vol. 19, no. 10, pp. 1821–1842, 2003. View at Publisher · View at Google Scholar · View at Scopus - L. Taleb and S. Petit, “New investigations on transformation induced plasticity and its interaction with classical plasticity,”
*International Journal of Plasticity*, vol. 22, no. 1, pp. 110–130, 2006. View at Publisher · View at Google Scholar · View at Scopus - G. Antoni,
*Transformations Tribologiques de Surface: Une Approche Thermo-Mécanique [Ph.D. thesis]*, Université de Provence, Marseille, France, 2010. - A. Chrysochoos, C. Licht, and R. Peyroux, “A one-dimensional thermomechanical modeling of phase change front propagation in a SMA monocrystal,”
*Comptes Rendus Mecanique*, vol. 331, no. 18, pp. 25–32, 2003. View at Google Scholar - J. Garrigues,
*Fondements de la Mécanique des Milieux Continus*, Hermes Science Publications, Paris, France, 2007. - K. Dang Van, M. H. Maitournam, and B. Prasil, “Elastoplastic analysis of repeated moving contact application to railways damage phenomena,”
*Wear*, vol. 196, no. 1-2, pp. 77–81, 1996. View at Google Scholar · View at Scopus