Thermodynamical Description of Running Discontinuities : Application to Friction and Wear

The friction and wear phenomena appear due to contact and relative motion between two solids. The evolution of contact conditions depends on loading conditions and mechanical behaviours. The wear phenomena are essentially characterized by a matter loss. Wear and friction are in interaction due to the fact that particles are detached from the solids. A complex medium appears as an interface having a strong effect on the friction condition. The purpose of this paper is to describe such phenomena taking account of different scales of modelization in order to derive some macroscopic laws. A thermodynamical approach is proposed and models of wear are analysed in this framework where the separation between the dissipation due to friction and that due to wear is made. Applications on different cases are presented.


Introduction
The wear phenomena are induced by contact and by relative motion between two solids.They depend on the loading conditions and on the materials properties; they are essentially characterized by a matter loss.
Particles are detached from the solids in contact and a complex medium takes place between the two solids forming a thin layer.The conditions of friction are then affected and therefore the conditions of wear evolve simultaneously.
The interface is a complex medium made of detached particles and eventually of a lubricant fluid.The damaged zones belonging to the solids can be also considered as a part of the interface.All these zones constitute a layer where the materials lose their own cohesion.
The evolution of this layer is complex especially in the transient phase of the interface formation.However, for particular geometries and for steady states, with a constant flux of matter, the wear states can be studied experimentally and conceptually.This framework is a first step of comprehensive study for understand the parameters governing friction and wear.
Following the terminology introduced in [1] the interface is called "third body" and must be considered as an aggregate thin layer of different particles with sometimes a lubricant fluid.This layer develops non-linear macroscopic rheology that must be characterized.
If macroscopic descriptions of such an interface are known in the literature [1,2], the connection with local mechanical quantities and discussion based on microscopic scale modelling are not currently developed, unless in some recent studies [3][4][5][6].
Some well known wear criteria, as Archard's law [7], are useful but such models cannot be predictive when the operating conditions are not sufficiently close to the experimental test conditions of their elaboration.
Experiments on wear and friction, as observed in Stribeck's curve, can provide relations between the friction coefficient µ and the lubricant coefficient L : µ = τ /p, L = ηU/p where τ is the shear stress, p the pressure of contact, η the fluid viscosity and U the relative velocity.This curve shows three particular regimes: the first regime corresponds to a Coulomb's friction law with constant coefficient, the second is an unstable regime with a strong decrease of the friction and the third is a linear regime corresponding to a state of mild wear.
These results suggest clearly a strong interaction between debris and the solids in relative motion.The modelization of this interaction can be explored through the evolution of specific internal state variables governing the third-body behaviour.
The simpler approach is made using of theory of mixtures, where the volume fraction of debris plays an important role.
The main purpose of this article is to propose a more general formulation using of thermodynamical considerations on wear phenomena studying the propagation of surfaces or layers inside sound bodies taking account of damage evolution and loss of sound matter.A thermodynamical approach of third body is then introduced and dissipation is analyzed.
For a macroscopic view point the behaviour of the interface is modelized by unit of surface of contact.The density of mass by surface unit is related to the volume fractions of detached particles and plays the role of an internal parameter.The driving force associated to wear is then deduced and criteria of wear are proposed.
Finally, particular situations are studied using some constitutive law of the interface.In each case the dissipation is analysed.The matter loss and the geometry change are determined according to specific criteria and associated laws.Analytical and semi-analytical solutions are presented.

General Features on Moving Surfaces and Moving Layers
The wear phenomena due to contact and to relative motion between two solids V 1 , V 2 depend on loading conditions and on material mechanical constitutive laws of the solids in contact.Wear is essentially characterized by a matter loss.Each solid V i is decomposed into a undamaged part Ω i and a damaged zone Ω i3 .The boundary between Ω i and Ω i3 is Particles are moved from the solids Ω i when some criteria are satisfied at the boundaries Γ i between Ω i and a complex medium Ω 3 is generated called the interface or the third body.
The third body Ω 3 is made of detached particles, damaged zone Ω i3 of the two solids as depicted in Figure 2.
The first step of modelization is to characterize the motion of the surface Γ i .The third body Ω 3 is then Ω 13 ∪ Ω 23 ∪ Ω 33 , where Ω 33 is a complex mixture of detached particles and eventually of a lubricant fluid.
Along the surface Γ i3 the displacement is continuous and the stress vector too.The local behaviour changes from those of the undamaged material Ω i to those of a damaged solid Ω i3 .The surface Γ i3 moves with the normal velocity c i = ϕ i ν i in the reference state of the solid, ν i is the outward unit normal vector to Ω i3 and then ϕ i is positive.The value of ϕ i is given by a constitutive law, this is a part of the wear description.
This transition can be brutal or diffuse.In the first situation, the transition is governed by the motion of a surface Γ i3 along which the material characteristics endure change by an irreversible process.Due to damage, the elastic moduli in domain Ω i3 are lower then those in Ω i , then during this brutal transformation the elastic moduli are discontinuous along Γ i3 .These discontinuities contribute to the dissipation by a surface term.
In the second situation, the variation of the quantities is smooth and the quantities remain continuous.The damaged zone Ω i3 constitutes a thin layer where the damage and the internal variables vary continuously.The boundary Γ i3 between the undamaged material and the damaged zone is moving but all the mechanical quantities are continuous, the surface term in the dissipation disappears.
Each situation implies different type of dissipation and for each case a driving force associated to the wear process is proposed.
The conditions of loading are chosen such that the inertia effects are small.Extension to dynamical problems uses dynamical concept for the driving force applied for defects motion [8], interface or surface motion [9][10][11][12].
The local constitutive law The state of each body is characterized by the displacement u, from which the strain field ε is derived.The other parameters are the temperature θ and a set of internal parameters α.The behaviour of V i is defined by the free energy density ψ as a function of strain ε, the temperature θ and the set of internal parameters α.The mass density of each phase (undamaged and damaged) is the same ρ.The state equations of each phase are The boundary Γ i3 , between Ω i and Ω i3 , is a perfect interface.The external boundary ∂Ω i /Γ i3 is decomposed into ∂V u and ∂V T on which the displacement u d and the loading T d are prescribed respectively.
The internal state of stresses satisfies the conservation of the momentum This state σ is decomposed in the reversible part σ r and an irreversible part σ ir .The irreversible stresses σ ir are essentially due to viscosity.In non linear mechanics, the internal state is generally associated with irreversibility.The evolution of internal state must satisfy the second law of thermodynamics.Such requirement is fulfilled by the existence of a potential of dissipation.
Analysis of the dissipation and potential of dissipation The fundamental inequality of thermodynamics implies that the internal production of entropy must be non negative.The equations of state do not provide all the constitutive equations; some complementary laws are necessary to describe the irreversibility.In the total dissipation D T the contribution of the conduction and those of internal forces must be distinguished.
The two parts are assumed to be separately non-negative.The mechanical part has the form: and using of equations of state, the thermodynamical forces associated to irreversibility are obtained and we recognize the role of the irreversible stresses σ ir .
To determine the evolution of the internal state, relation between velocity and driving forces must be given.These relations must be compatible with the positivity of the internal production of entropy D m .
Let us assume that the local behaviour belongs to the class of the so-called "generalized standard" materials [13].For this class, a potential of dissipation d( ε, α) defines the local irreversibility.The potential d is a convex function of these arguments with a minimum value at the origin.The evolution of the internal state is given by the normality rule: this means that the subdifferential ∂d( ε, α) of d is the set of state (σ ir , A) such that: for all admissible fields (ε * , α * ).The existence of such a potential for the dissipation ensures the positivity of the entropy production:

The conservation laws
Classical laws of conservation are written, taking discontinuities into account: Energy balance where e is the internal energy e = w + sθ, s is the entropy, q the heat flux

On moving interface
To study this problem, the analysis is made on one body V i .The solid is decomposed into Ω i and Ω i3 .Along the boundary Γ i3 , ν i is the outward unit normal vector to Ω i3 , c is the normal velocity and ϕ i is the magnitude of c.
When the surface Γ i3 is moving, a mechanical quantity f can suffer a discontinuity The volume V i is decomposed in two domains, separated by a moving boundary Γ i3 .The volume average of f has a rate defined by Entropy production Using the definition of the entropy of the system S = ∫ V i ρs dΩ the second law of thermodynamic is written as and after integration by part, the value of the volume dissipation D T is recovered and a surface contribution appears Taking the balance equation of energy and the continuity of temperature into account, the dissipation by unit of surface is obtained This dissipation is a linear function of the normal velocity ϕ i of Γ i3 .The quantity G i is the driving force associated to the movement of Γ i3 .The displacement is continuous along Γ i3 then the discontinuity [∇u] Γ i satisfies the Hadamard's relation and the release rate of energy G i becomes The condition m i = 0, which explains that there is no loss of matter, induces that ϕ i = 0 then there is no discontinuity [v] Γ i = 0, so there is no dissipation along the surface Γ i3 , because this surface does not move.
Interface propagation law To control the matter loss, a criterion based on G i can be formulated.For example, we can considered a Griffith's type law for the propagation of the moving interface When the interface moves hence G i = G c , this quantity must be conserved during the motion.
Convected differentiation To study the evolution of the mechanical state around the moving surface Γ i3 a convected derivative of any mechanical quantity f is needed.A point X Γ is on Γ i3 if its coordinates satisfy the scalar equation: It is well known that the gradient ∇S(X Γ , t) For any function In particular, the evolution of the middle surface satisfies and the convected derivatives of X Γ and ν i are The convected derivative is useful for write the Hadamard's equation of compatibility: for any continuous function For the displacement, we have The convected derivative of u(X + Γ , t) = u + and of u(X − Γ , t) = u − are given by the classical relation on discontinuity is recovered as Consistency condition The law of propagation shows that the condition ϕ i ̸ = 0 implies G i = G c , that relation is conserved during the propagation.The consistency condition associated with the Griffith's propagation law is then written at point s ∈ Γ i3 where G i = G c in the form This condition determines the velocity ϕ i (X Γ , t) ; the surface equation is then implicitly defined by the condition This description is based on the fact that the transition between the initial solid and the damaged material is brutal, that transition induces discontinuities.

On moving layer
The situation of diffuse damage is now considered.For example, a damage constitutive law governed by a continuous scalar function d is assumed and a free energy is taken as the damage parameter varies from 0 (no damage) to 1 (totally damaged).Over Ω i , d = 0 and d varies continuously from 0 to 1 over Ω i3 .
In this description along the surface Γ i3 , between sound material and damaged zone, d = 0.The stress vector ([σ] Γ i .νi = 0), the displacement and the elastic moduli are continuous quantities.These properties induce that ∇u is continuous and so the energy release rates G i is zero along the surface Γ i3 .The previous results on dissipation must be revisited.The loss of matter is governed by the motion of the surface d(x, t) = 1.
The volume V i is decomposed as previously, but now the displacement and the internal parameters are continuous along the surface Γ i3 .Inside the domain Ω i3 , the evolution of a mechanical quantity f (ε, α) function of the strain ε and of the set of internal parameters α must be characterized.The volume average of f is decomposed as and has a time-derivative expressed in term of a classical time derivative in Ω i and a convected derivative on Ω i3 .Convective derivative is useful to describe hypothesis of self-similarity process in the damaged zone during the wear.The domain Ω i3 evolves simultaneously with two surfaces: one is Γ i3 associated with the boundary d(x, t) = 0 and the other is Γ o along which d(x, t) = 1.
A point x in Ω i3 is locally defined by it's normal projection X Γ on the surface Γ i3 and the normal coordinate z along the unit normal vector ν i .The thickness of the layer is H(S, t) and 0 ≤ z ≤ H.
When the surface Γ i3 moves with the normal velocity ϕ i , the local frame τ i , ν i moves simultaneously.For example a 2D motion of a curve satisfies the geometric relations where γ i is the curvature of the surface Γ i3 at point S, ν i the normal unit vector et τ i the tangent unit vector to the curve.It can be noticed that the elementary length is function of the curvature in the reference coordinates (X Γ , z).Similar relations exist for a surface motion.
The point x has a convected derivative During the motion, the thickness evolves and the velocity of the point M = x(H) is The evolution of each term of F can now be given.
In these expressions, The function f is continuous along Γ, then there is no contribution of discontinuities on the time derivative of F , then The only surface contribution is due to the flux of matter along the surface Γ o .Along Γ o , the loss matter occurs simultaneously with the full damage d = 1.
Over Ω i3 , the convected derivative of f is used and we have ∫ then the rate of F is finally The global dissipation is expressed using these expressions.

Analysis of dissipation The dissipation has the classical expression
as previously we have The dissipation contains to terms, one in the volume and the last one defined on the surface.In the transient phase this terms is zero, because the material in Ω i3 is only partially damaged d < 1, then D Γ M.ν H = 0.When rupture occurs, this term contributes to the dissipation, the value of energy at point d = 1 is lost.Some others expressions for dissipation can be formulated.This tensor satisfies the relations

Others expressions of dissipation
And using the convected frame with Ω i3 with the velocity D Γ x, the dissipation is rewritten as This expression has a great interest when local conditions of the stationarity of the steady state is realized.If the damage d is a function of the distance z to the curve Γ i3 the local stationarity is expressed as and for local stationarity, D Γ α = 0.

Comments on local stationarity
If Y is the driving force associated to d, then the contribution of damage to dissipation is Therefore, if the distribution of ∇d is known inside the layer, the average value Ḡi is the driving force associated to the velocity ϕ i , The quantity j(x) takes the shape of the layer into account and it defines the measure of area at x with respect to the measure of area at X Γ so that dΩ(x) = ds dz = j(x) dS dz (58) Assuming now that the profile d(z) is a given function, the propagation law is generalized as Combining the hypothesis of stationarity and an imposed profile for d, the thickness H is given, the Γ o interface is moving with the normal velocity ϕ i of Γ i3 , when − ∫ H ∇d.ν i j(x)dz = 1.In this case, the iso-d curves are orthogonal to the normal ν i and each iso-d curve has the same velocity ϕ i .

The Macroscopic Interface Study
From a macroscopic point of view, it is not necessary to model the microscopic behaviour which occurs in Ω 3 .Micromechanical considerations can be useful in order to describe in a realistic manner the behaviour of the interface.This type of analysis can provide a relevant interface law at the mesoscopic level and at macroscopic scale too.The multi scale approach of the interface is based on fundamental characteristics at three different scales: • At the microscopic scale, the contact between asperities govern the wear mechanisms.Some studies with plastic strains and with micro cracks propagation have been tempted [5,14].These are the fundamental ideas of the decomposition of Ω 3 in Ω 13 ∪ Ω 33 ∪ Ω 23 .
• At the mesoscopic scale, this is the description of the third body.This model was proposed in [1] developed by [15] and [16].The local physics is that Ω i3 is a porous medium.In Ω 33 the solid particles are in suspension forming a shear layer as inside a viscous fluid flow.This point is developed in [3].In other situations, a very large plastic shear deformation is present, so the profile of local deformation in this thin layer depends strongly on the loading conditions.
• The macroscopic modelling is based on models of friction law which are depending on parameters to account of the evolution of the interface.These models can be inferred from the smaller scale by some averaging technics as those proposed under local stationarity hypothesis.
Some other studies are founded on cyclic loading taking account of cyclic asymptotic behaviour of elastoplastic materials generalised to the theory of damage in mechanics [17].
The system of two bodies in relative motion separated by the third body Ω 3 is considered.As previously in each domain the local behaviour is described by a free energy ψ and a potential d of dissipation.The dissipation is considered only for domain Ω 3 .So the behaviour in Ω 1 and Ω 2 is reversible.Therefore dissipation is concentrate into the interface.
The proposed description is depicted in Figure 3.
Sharp transition At the mesoscopic level, the dissipation is given by where , and d 3 is the volume dissipation due to the irreversibility processes inside Ω 3 : Comments If matter loss occurs, m 1 or m 2 is positive, then the corresponding discontinuity [v] Γ i exists too and the local quantity There is a dissipation due to the loss of material: that is a characterization of wear.The volume dissipation d 3 contains two contributions, one is due to heat conduction, the other is associated to mechanical irreversibility.This contribution can be associated with the viscosity of the fluid which carry the detached particle.When the shear reaches a critical value τ c the velocity field inside the interface is reduced to a gliding motion v = U/H ze x and the dissipation in the mesoscopic scale is evaluated as The second term is due to the viscoplasticity governed by the yield stress τ c , it plays the role of friction [18,19].The control of the wear dissipation obeys to the normality rule This simple description is based on the brutal transition on the mechanical characteristic from initial material to a damaged one.The interface Γ i3 moves according to this normality rule.
Diffuse damage Consider now that Ω i3 is a zone of diffuse damage.The description of damage can be made with a smooth transition governed by a damage parameter d in Ω i3 , a set of internal parameters α, which contains the plastic strain ε p .For example a free energy of the form is considered.In this case, the damage parameter governs continuously the change of moduli of elasticity and it describes the degradation of the stiffness.
As previously, the state equations are given by and a potential of dissipation d i ( ε, α, ḋ) is given so that The thickness of the interface At this stage possible choice for the local constitutive behaviour can be used.And literature provide many papers on this subject with computational results.But the main difficulty is to determine the physical characteristics of the behaviour at this scale.
A macro level approach based on this local description can be tempted to understand the specific structure of the thin layer and to estimate some necessary parameters to a relevant description of the mechanisms inside the interface.

A Macrolevel Approach
The preceding sections emphasise the fact that in a more general global approach, we can tempt to characterize the thin layer behaviour by considering it is in a state of equilibrium under given conditions applied to Γ i3 .Then the behaviour of the thin layer can be considered as affected to a surface Γ, by considering that the thickness is an infinitesimal quantity compared to a characteristic length of the bodies in interaction.Time derivatives of integral of varying surface or volume domains are studied in [20].Some such results are used here to express the dissipation.
The interface Ω 3 is a thin layer composed by Ω 13 ∪ Ω 33 ∪ Ω 23 .The geometry of the layer is described by its middle surface Γ and the thickness H = 2h at each point of this surface.
The free energy ψ s per unit of area of the middle surface Γ of the layer is defined by x = X Γ + zν, X Γ is the position of the middle surface, z is the normal coordinate −h ≤ z ≤ h, the middle surface has a curvature κ, and j(x) is the variation of area surface due to curvature with respect to the middle surface coordinates.The local volume is then dΩ = j(x) dz dS.The mass ρ s by unit of area is defined by The mass of the overall system is not conserved.The mass loss ρ i ϕ i from Ω i is added to Ω 3 hence we obtain d dt But in the same time, a flux of matter is present in the tangent plane of Γ with the velocity C s Then the rate of the surface density is given by This contributes locally to the matter flux.Assuming that the conservation is locally true, we have The local problem The layer is in equilibrium with external loading and contact conditions.The displacement is continuous on Γ i , the surface energy ψ s depends upon the given displacements along Γ i .The other parameters are a set of internal parameters α(X Γ , z), the thickness H(X Γ ) and the temperature.The local strain ε derives from the displacement u 3 .The displacement is continuous along the interfaces Γ i then Then, it is obvious that the surface free energy is a function of u i (x i , t) and of the thickness H. Time derivatives of integral of varying surface or volume domains are studied in [20].Some such results are used here to express the dissipation.
The strain is rate is ε = 1 2 (∇v + ∇ T v), where v = u.The dissipative function D s associated to the local potential d 3 = d( ε, α) is defined by the average We consider that the behaviour inside Ω i is reversible.The steady state equations over Ω i are then For the interface, we have where T r i is the reversible tension along the interface.The balance momentum equations using of the local behaviour is given by and the complementary laws are The normal ν i along the interface Γ i3 is not the normal to the middle surface because h depends on X Γ .
The macroscopic description as depicted on 4 ignores the details at the microscopic level.
The local thickness H(X Γ ) = 2h is obviously the relevant parameter geometrically associated with the interface, while the internal parameters govern the physical properties of the layer.For example, the set of volume fraction of debris from materials 1 and 2 are such parameters.

Analysis of dissipation
The global dissipation is derived from the macroscopic description.It is useful to introduce the global free energy P of the tribologic system The rate of the free energy is then which becomes after rearrangements where γ is the mean curvature of the middle surface and ϕ the normal velocity of the mean surface.
Taking account of equilibrium and conservation of mass, we have The displacement is continuous along the interfaces Γ i3 then where The total dissipation D m is hence where A term due to the variation of the thickness appears But ϕ i , H are dependent.The velocities of propagation are linked And the dissipation is equivalent to the second term is due to D Γ ν and it is very small because H is small, and Expansion of displacement with respect to z This macroscopic point of view suggests to develop the internal state over Ω 3 as an asymptotic expansion of the coordinate z The relative displacement w = u 1 − u 2 is then a function of the gradient of u on the middle surface at the first order in z.
The free energy is expanded with to z and the expression is given by where ε = ε o + zε 1 at first order.In the same way, the potential of dissipation is Under this approximation so the energy depends on the value of w and on the gradient of the displacement along Γ.
In this case, the global state of equilibrium is revisited to take into account of the dependance of the energy with gradient of displacement ρ s ψ s = f (u i , ∇u i ).The variations of ∫ Γ ρ s ψ s dS with respect to u i are then by integration by parts we obtain Hence the equilibrium of the interface is now given by if there is no viscosity, if not the right part contains new terms This approximation shows that the conditions of continuity of the displacement given by the asymptotic expansion and the displacement of the undamaged bodies implies that the imposed shear U/H e x is related to the discontinuity v 1 − v 2 of the two bodies in contact.These expressions show that global approach based on relative displacement or relative rate of displacement can be justified by micromechanical considerations, but they emphasise the role plays by the mass density by unit of area.This quantity is an important internal parameter which is governed by the mass conservation and the wear criterion.
An important task is to develop models of the layer based on micromechanical hypotheses.Using homogenization theory for thin layer under global shear loading, taking account of the mass flux, must be investigated.This process will provide different models of interface behaviour depending on the dissipative mechanisms evolving inside the third body.These models must be derived taking account of conditions of sliding contact, in particular in the presence of viscous fluid, the flow between the bodies must be characterized.This shows the emergence of specific time scales according to each models of interface depending on loading conditions.The mechanisms of degradation can also be modified during the loading history.This approach should be useful for determination of the domains described by the Stribeck's curve or also to study the transition from a regular state of low rate of wear to a state of abrasive wear.

Examples and Applications
A linear elastic half-plane in plane strain is considered.The purpose of the study is to analyse the contact wear under a rigid punch.The studies are made in two cases of loading : • the sliding wear under steady relative motion [3], • the sliding wear under cyclic loading [17].
In order to obtain displacements, stresses and strains at the surface of the half-space, which is covered in the contact area by the interface model, it is useful to consider integral equations.
The displacement (u x , u y ) of the upper boundary of an linear elastic half-space with elasticity characteristics Young modulus E, Poisson's ratio ν, satisfies the integral equations [24], for which the contact area is (−a, a): where the constants c i are and the principal value P v is defined as

Sliding contact in steady relative motion
The rigid punch has a vertical displacement u p y = δ + x 2 /2R and we assume that wear occurs only in the half space.Ahead the punch there is no debris, then the volume fraction of debris f = 0.The thickness of the interface is H o , which corresponds to the sum of rugosities of the solids and contains the thickness of the incompressible fluid.Due to wear, the thickness evolves.The mass conservation and the fluid incompressibility give the relations between the wear rate ϕ(x), the fraction of debris f (x) and the thickness of the thin layer H(x).
All the equations of conservation are written in the moving frame with the punch at the velocity U e by integration over the thickness H we have and because the layer is submitted to a local shear v = U/Hze x we obtain The mass of fluid is also conserved then The constitutive law The free energy of the mixture is given by A potential of dissipation is given to determine the irreversible contribution, essentially due to viscosity This constitutive law generalizes the law use in [3] in which η a = k t = 0.For (η a = k t = 0) we have an interfacial behaviour given by k(f ) and η(f ) are chosen from typical homogenized value of the phases.
that the homogenized Reuss's model for the stiffness, and the Einstein's law for the viscosity.Introducing these equations in the equilibrium equation for a given profile ϕ(x) the answer of the half-space will be determined.The wear rate ϕ must satisfy a complementary law, as proposed before.For the sake of simplicity we take ϕ = k p σ yy (109) This relation determines the velocity ϕ.The solution is obtained analytically by an asymptotic expansion in series of the volume fraction f of particles.
• At zero order, the Hertz's solution is recovered • At first order, a dependance with f is obtained.Wear occurs, and the profile of the pressure σ yy (x) evolves.The presence of viscous fluid induces a displacement of the maximum of pressure like under the dry contact with friction [25].
This analytical solution is studied in paper [3] and discussed in [26].

Cyclic loading
The matter loss is determined with respect a criterion of wear in the case of cyclic loading.For this the results of cyclic plasticity are generalized for elastic brittle materials.
The vertical displacement δ of the punch is prescribed and a periodic horizontal displacement of the punch is imposed.
The wear rate is related to the Griffith's law: During the loading, wear occurs, the surface evolves and reaches an asymptotic shape Γ.
The loss of matter is determined by the surface shape y = Ψ(x) of the half space such that the volume of loss matter is the smaller volume compatible with the criterion of wear: during the motion of the punch the criterion (G(y, t) − G c ≤ 0, ∀y = Ψ(x)) is fulfilled all time.
Assume that ϕ = Ψ(x) is small.During the motion, the applied loading is cyclic and is convected by the motion to the surface.Using of derivative D Γ , on the initial configuration the loading is given by σ.ν = F (x, t) and is now convected along the moving interface with the velocity ϕ, then The rate solution is determined on the initial geometry, with components of stress and displacement rate given by these relations.
Since the variation of the geometry is small, the kernel N o (x, z) of the integral equation must be developed at least at first order in ϕ [17].
For an isotropic material, the solution for displacement is found to be: This problem is solved using for the driving force G = 1 2 ε : C : ε−n(C : ε).∇u.nwhere n is the normal vector to the surface (x, ϕ(x)).In the considered case, the indenter is rigid so there is no contribution of the indenter to the driving force.
The stabilized answer is obtained numerically, using an iterative algorithm [17].The conditions of contact include Signorini's condition with friction and comparisons are made with experiments [27].

Conclusions
In this article, a review of different descriptions of wear is presented.The approach emphasise the role of micromechanics on the elaboration of a relevant modelisation of wear by introducing the most important parameters in the description of the third body.The modelization shows that the mass density by unit of surface plays an crucial role, especially in transient phase.
The thermodynamical analysis suggests also to consider local criterion for wear based on driving forces associated with the local phenomena involving the final rupture.At different scale the relations between local driving forces and global driving forces to describe the behaviour of the thin layer must be established solving a complex boundary value problem.This problem is a problem of equilibrium of a thin layer under macroscopic shear deformation with non-linear behaviour including damage and viscoplasticity.According to non-linear behaviour and to presence of fluid, different time scales must appear determining domains of validity of each model of interface.
This multi-scale approach has been shortly discussed.According to different choice of modelisation, different expressions for dissipation are obtained.The driving force associated to wear process have been expressed and wear laws have been proposed.
Simpler formulations have been achieved, taking account of the most important features of the wear process and semi-explicit answers are obtained.These results show the ability of such an approach.Some of the ideas and techniques used here (like integral equations) could be of some interest for solving other problems on wear.

Figure 2 .
Figure 2. The microscopic description of wear process.

Figure 3 .
Figure 3.The transition from microscopic to mesoscopic description of wear process.

Figure 4 .
Figure 4.The transition from mesoscopic description to macroscopic description of wear process.