Completely Analytical Tools for the next Generation of Surface and Coating Optimization

Usually, some severe efforts are required to obtain tribological parameters like Archard's wear depth parameter k d. Complex tribological experiments have to be performed and analyzed. The paper features an approach where such parameters are extracted from effective interaction potentials in combination with more physical-oriented measurements, such as Nanoindentation and physical scratch. Thereby, the effective potentials are built up and fed from such tests. By using effective material potentials one can derive critical loading situations leading to failure (decomposition strength) for any contact situation. A subsequent connection of these decomposition or failure states with the corresponding stress or strain distributions allows the development of rather comprehensive tribological parameter models, applicable in wear and fatigue simulations, as demonstrated in this work. From this, a new relatively general wear model has been developed on the basis of the effective indenter concept by using the extended Hertzian approach for a great variety of loading situations. The models do not only allow to analyze certain tribological experiments, such as the well known pin-on disk test or the more recently developed nano-fretting test, but also to forward simulate such tests and even give hints for structured optimization or result in better component lifetime prediction. The work will show how the procedure has to be applied in general and a small selection of practical examples will be presented.


Introduction
In order to achieve the goal set in the headline, namely, to analytically optimize surface and coating structures, which includes the extraction of generic tribological or wear parameters, it is necessary to combine quite a few fields and/or concepts of material science.Therefore, this introduction needs to cover the following issues:  first principle based effective interatomic potential description of mechanical material behavior [1] and its extension into the time-domain [2];  the effective indenter concept [3][4][5][6], made time dependent;  the extension of the Oliver and Pharr method to analyze nanoindentation data to layered materials and time dependent mechanical behavior [7][8][9][10];  the physical scratch and/or tribological test and its analysis [11][12][13][14].

Simple First Principle Based Interatomic Potential Description of Mechanical Material Behavior
It was shown in [15] that on the basis of an effective potential function like the Morse potential given as A contact problem can be evaluated using the mechanical parameters derived from such a potential.Here, p, , σ are material parameters and r 0 usually denotes the equilibrium bond length.In such a case, the potential would define the pair interaction.Here, however, as in [1,15], we will apply the potential as an effective one with r 0 denoting the lattice constant (see also [16]).With respect to molecular dynamic simulation such an effective potential could be the basis for the extraction of the necessary pair and higher order interactions as demonstrated in [16].For our study, however, we will not need this, because we are only interested in the mechanical constants, especially the Young's modulus and the decomposition strength.Having this, one can apply the method described in [15] to simulate a mechanical contact problem, thereby even taking nonlinear effects like the pressure, shear and temperature dependency of the Young's modulus, and other mechanical constants into account [2].We will show that especially within tribological applications temperature fields can play an important role in influencing the mechanical properties of the material.

Brief Story about the "Effective Indenter Concept" and Its Extension to Layered Materials
"Surface" usually is not "bulk" and, so, subsequently, the behavior of surfaces, especially their mechanical behavior can be dramatically different from what one might expect by just applying bulk concepts to surface problems, such as contact situations with very surface-located, surface-dominated stress and strain fields.In tribo and wear problems, this very often compromises our ability to simulate and understand the physical processes taking place in certain tests (e.g., [17][18][19][20][21]).In order to improve this situation we, therefore, resort to a stringent application of a layered material model considering and modeling the mechanical surface always as-at least potentially-having a property profile starting from the mechanical bulk values in depth usually well known for a certain material to rather often completely different properties on the top-most surface layer.In the case of coated materials this model extension, of course, can easily be justified by the explicit coating structure.However, we should stress the point of also carefully and critically considering apparent "homogenous surfaces" as-rather often-being of gradient or somehow layered character.In addition tribological applications often come with temperature fields due to friction, be it the internal friction caused by shear stresses or the external one coming from the multiple tribological contact situations.As shown in [2], such temperature fields can couple back into the mechanical properties influencing the tribological contact conditions.
The "effective indenter concept" itself can mathematically be described or understood as some kind of quasi conform coordinate transformation transforming the difficult problem of a curved surface contacted by a well defined indenter (like a cone) into a flat surface loaded with a complexly formed indenter.Already, in 1994, Bolshakov, Oliver, and Pharr [3] introduced this "Concept of the Effective Indenter" and refined it in a series of wonderful publications until, in 2002, the paper about "Understanding of nanoindentation unloading curves" [4] was published by Pharr and Bolshakov.The extension of this concept is simply performed by substituting the homogenous half space model describing the loaded sample body by a layered half space model [22].

About the Extension of the Oliver and Pharr Method to Analyze Nanoindentation Data of Layered Materials and Time Dependent Mechanical Behavior
As it is a well established fact, that the classical Oliver and Pharr method [7], as an approach based upon the homogenous half space model, cannot directly be applied to layered materials and small structures, the author here refers to the literature (e.g., [23]).Soon after the publication of the Oliver and Pharr method it became clear that there is a physical concept on which this method can be based.It was called the "Concept of the Effective Indenter" [3,4].During a small conference in Italy in 1999, the author learned of that concept and decided to work out a theory solving not only the problem for the mechanical contact of an indenter with general shape of symmetry of revolution, but also to extend this solution to layered structures [24].Later on, he also presented a variety of applications together with Pharr, Chudoba, Richter and others [5,6,8,9,[25][26][27].Soon, it became clear, however, that the "effective indenter theory", even though powerful, was not something one could easily give to the engineer or an indenter experimentalist and expect her or him to use it as a tool for the analysis of indentation data.The reason for this lays in the complexity of the formulae building the solution.Thus, the whole approach was brought into a software package named FilmDoctor ® [28].

Making the Classical Oliver and Pharr Method Fit for Time Dependent Mechanical Behavior
In an indentation test we always find complex three-dimensional stress states with usually all stress components being non-zero.Thus, as Fischer-Cripps put it, "the nature of the loading is a complex mixture of hydrostatic compression, tension, and shear" [29].In the case of viscous behavior we also have to understand that we even have complex mixtures of stress and strain rates.These however, usually influence the time dependent mechanical parameters [30].This automatically means that at different positions we find different stress and strain states and rates subsequently leading to different mechanical parameters (time dependent) at different positions within the material.This automatically makes the system of usually linear partial differential equations of Elasticity non-linear in the moment viscosity (even elastic viscosity) comes into play.In order to keep things simple, however, we will not go here for a non-linear basic solution, but try to find effective, phenomenological descriptions for visco-elastic, visco-plastic contact problems, based on the concept of a time dependent effective indenter.As an additional difficulty coming into play especially in connection with tribological applications we again have to consider temperature fields possibly influencing the mechanical parameters.It was shown in [2] that even smallest temperature gradients can have a dramatic influence and as tribological contacts usually and mainly inevitably produce heat this effect cannot be ignored when trying to physically simulate tribo experiments.

Introduction into the Physical Scratch and/or Tribological Test and Its Analysis
The standard scratch test is a widely used method to test the mechanical stability of coatings on different types of substrates and has become a sensitive technique to control the reliability of the manufacturing process.It is based on various standards [31,32].
Unfortunately, these standards do not allow to design and/or perform the tests in a truly physical way, meaning to extract generic physical parameters being later applicable in life-time prediction and optimization.For this, a more sophisticated approach is required, which assures physical parameter identification at each and every step towards physical tribological tests.The procedure is outlined in Figure 1.
Following the flowchart in Figure 1 we have started our surface characterization process with nanoindentation.The section above covered the calculation of true mechanical (potentially layered) surface properties out of these nanoindentation tests.Thus, now all mechanical values are known which are not only required to properly dimension (fine-tune) a scratch or tribology test for specific surface structures, but also are of need for a physical analysis of these tests.Again the flow chart (Figure 1) might give an illustrative understanding into such a test procedure.
The reader will find illustrative examples and a much more comprehensive elaboration of the method elsewhere [11][12][13][14].

Figure 1.
A flow chart of the procedure of mechanical characterization and optimization of arbitrary structured surfaces with respect to wear and fretting.See references about features and tools [28,[33][34][35][36] been given in the flow chart.

First Principle Based Interatomic Potential Description of Mechanical Material Behavior
As shown in [1], the pressure P and bulk modulus B can be derived from an effective potential.Thereby it is convenient to express P and B in units of B 0 (B at P = 0).By substituting the lattice distance r by r = c × r 0 , P and B result in the relations: ; (2) As estimates for the critical r leading to decomposition (cf.[1]) we take the following expressions which have to be solved numerically (r 00 = p × r 0 ) to extract the critical c-value c m for maximum P(c): ( As a purely mathematically based measure for the critical bond length or in our case of an effective potential the lattice distance, the inflexion point for c ifp > c m could be used.This can be numerically obtained for the Morse potential via: (4) As shown in [2], the presence of deviatoric or shear stress components reduces the theoretical strength towards tension by up to 10%.This is of special importance for multiaxial contact situations where the stress tensor always is of mixed form and fully filled.
One also needs to check for temperature gradients caused by external and internal friction due to the tribological shear.It was shown in [2] that, caused by the inhomogeneity of the temperature T(t,r) and subsequently the relaxation τ(t,r) field, an originally simple approach for the bond interaction, e.g., the Standard Linear Solid (SLS) becomes rather complicated regarding its bulk modulus behavior.Here, t and r are denoting time and space coordinates, respectively.Instead of a phenomenological SLS with the well known bulk modulus of the kind: (5) (N gives the number of particles and c defines a parameter connected with the volume of an individual particle) and a subsequent Young's modulus time dependency like we have a resulting phenomenological Young's modulus structure of the form: (6) With a known principle structure for the temperature distribution function deduced from the shear stress distribution (adiabatic approach) such an approach seems to be manageable by the means of iterations.With a typical dependency for the temperature dependent relaxation time like [37,38]: (7) (with a constant W A and the T = T 0viscosity η 0 and the Boltzmann-constant k B ) one can evaluate the mechanical contact situations applying the techniques elaborated in [2] (theory part II) under the assumption of completely adiabatic temperature fields resulting from shear and completely following the second stress invariant distribution.Please note, the positive sign in the exponential function for the relaxation time assuring decreasing viscosity η with increasing temperature.

The Effective Indenter Concept
In order to have a sufficiently great variability for the definition of differently shaped "effective indenters", we apply the extended Hertzian approach as shown in [39] for instance.With this approach normal and even tangential load distributions of the form [24] with n = 0, 2, 4, 6 and arbitrary constants c can be solved completely (by following the instructions of the mathematical procedures for obtaining the complete potential functions as given in [24] even arbitrarily high, but only even N).
Together with lateral loads (occurring in all scratch-, tribotesters, or the next generation of nanoindenters and their applications, see, e.g., [39][40][41][42]) one often faces tilting moments leading to a normal surface stress distribution of the form: These stresses can for example occur when the indenter shaft is dragged over the surface.As the shaft itself is elastic and, thus, would be bent during the lateral loading, an unavoidable tilting moment results and acts on the contacted surface.In addition, curved surfaces (e.g., due to roughness) can lead to such tilting moments.
In order to simulate the internal complex material structure of porous or composite materials, certain defect fields must be developed and combined with the external loads.Such fields can be extracted from the mathematical apparatus originally developed for the effective indenter contact by a simple generalization procedure.By developing such a defect model, one also obtains a very comprehensive tool for the construction of relatively general intrinsic stress distributions caused by internal inhomogeneities.Circular, disc-like inclusions could for example be simulated by the use of plane defects within the layered half space.Thus, introducing circular defects of radii a i of the loading type: (with x i , y i , z i denoting the center of the defect and n = 0, 2, 4, 6) directly allows us the application of the extended Hertzian approach [24] that provides a complete solution of the elastic field of the defect loading given above.By superposing a multitude of such "defect dots", one could model (simulate) a very great variety of material inhomogeneities and intrinsic stress distributions.The evaluation of the complete elastic field is straight forward.It only requires the evaluation of certain derivatives of the potential functions, given in [24].Finally, we need to take into account the curvature of the surfaces in order to its effect on the resulting contact pressure distribution [40].As the theoretical approach would not find enough space in this short note, the author will publish the necessary details elsewhere.However, the interested reader may derive the results presented here by comparing the solutions of the Laplace equation in Cartesian and Paraboloidal coordinates.

The Extension of the Oliver and Pharr Method to Analyze Nanoindentation Data to Layered Materials and Time Dependent Mechanical Behavior
Oliver and Pharr [7] have shown that the force removal curve from elasto-plastic deformations with a Berkovich indenter can be described by a power function: (13) Applying now the concept of the effective indenter, as shown in Figure 2, one can easily deduce that such an unloading curve can be connected with an Indenter of the shape.(14) with: (15) The indices i and s are standing for the indenter and sample, respectively.The evaluation above is based on the well known results of Sneddon [43], which require a linear elastic model.For this paper, we need to reformulate the basic contact equation incorporating the plastically deformed surface of the sample with the locally defined shape function (Figure 3): Together with the following approach for the indenter part being in contact with the sample surface of the kind: we obtain: as the new governing contact equation with the parameters B = B I + B S and n defining the effective indenter shape Z(r), as given above.Similarly, also more general effective indenter concepts can be derived, such as this one: ( ) ( ) based on a locally paraboloid surface with r-terms r 0 , r 2 , r 4 , r 6 , and r 8 .The reader might easily recognize the extended Hertzian character [6] of the basic contact equation given above.The simpler Hertzian contact would read: Some of the following examples and discussions will be based on these more general approaches.However, due to the wide use of the power law fit given above we will explicitly concentrate on the power law approach in connection with our practical examples.
In [6] and [24], it is shown how these general contact approaches can be applied and how the complete elastic fields have to be evaluated while the extension to the case of layered materials is been elaborated in [22].Extension to tilting and lateral loads is given in [6] and [24].
The next extension of need now is a time dependent analysis method for ordinary quasi-static nanoindentation tests (e.g., [2,44] with application [45]).Our approach will be of the following kind [29]: with t denoting the time.The next step is the introduction of a time dependent material model.Here, for the reason of simplicity, we at first resort to the well known three parameter approach given by a Young's modulus of the following kind Now we substitute the time dependent Young's modulus into the function C(t) of Equation ( 21) and obtain: Please note that in cases of substantial preloading situations integration with respect to the time t is required.
This makes the classical fit of the three constants h 0 , m and C a time dependent 6-parameter fit of h 0 , n, B and the material constants E 0 , E 1 and τ.It is clear that so many parameters automatically lead to numerical difficulties and instabilities if applied on data not providing sufficient linear independence.Thus, even in the pure visco-plastic case it is strongly suggested to use a variety of unloading curves (3 might be good number) obtained with different unloading speeds or at different maximum loads with similar unloading times.
In cases of time dependent inelastic behavior (like visco-plasticity) or more complex constitutive laws also n, B and h 0 have to be taken as n(t), B(t), h 0 (t) making the fit even more complex.The same holds in the visco-elastic case for just h 0 when there are great differences between unloading time and τ (comparable strain rates).The parameters n and B on the other hand are then only geometrical parameters and do not explicitly depend on time.A more detailed derivation and discussion of this extension can be found in [2].
In recent years more and more experimental concepts have been introduced trying to generalize the classical normal indentation process into a multiaxial, combined tilted or twisted contact test [46][47][48][49][50][51].In such cases also asymmetric and tilted contact situations need to be manageable.We also need to address the question of how to handle gradient structures, because especially with respect to tribo-protective polymer coatings gradient layers are of use.However, as closed form solutions of this kind are unwieldy and difficult to handle numerically (cf.[52,53] regarding elliptical contact situation or [54][55][56] in the case of a certain gradient structure) or simply do not exist in sufficiently general manner, we prefer the use of load dots as introduced in [57] or [58].With this approach, a solution of the problem given above can be constructed by combining the superposition of various elastic or visco-elastic spaces as given in [57] with the introduction of interconnected load dots [57,58].In our case, the new elastic spaces and loads to be superposed would be differently loaded elastic spaces at suitably chosen time frames with adapted viscose properties along the virtual time-axis.While each frame satisfies a homogeneous elastic problem, the superposition does not necessarily, at least not if the superposition is done in a non-linear way.For some the concept of superposed spaces might be a bit strange, but we remind the reader that in other fields of physics even stranger objects are superposed, such as living and dead cats in quantum mechanics, for instance.This way an approximated solution even for the nonlinear problem given above could be constructed.The method should be demonstrated by the means of interconnected Hertzian load dots (see Appendix) with different contact radii a i [58] being made time dependent as follows: Here we have used the complex presentation of the lateral displacements with u c = u + iv (with , u displacement in x-and v in y-direction, while w gives the displacement in z-direction).The coordinates r i have to be understood as with x i and y i defining the positions of the various load dots.The material functions, respectively, constants H i and , are defined through: the parameters  i have to be fitted to the integral contact equations.They are considered to be only time-dependent, meaning  I =  i (t).The reader should note that this approach is just a superposition of Hertzian loads with different contact radii a i at various positions possibly acting on viscose-elastic spaces being characterized by H i = H i (t).From these displacements the strain tensor can be evaluated in the usual manner by using: The stress components in the case of linear elasticity can be found using the following identities: with , , , ;; In the visco-elastic case applying a 3-parameter Standard Linear Solid model given by a Young's modulus of the following kind E(t) = E 0 + E 1 exp(−t/τ) and assuming a strain driven loading situation we have instead of Equation ( 27): with the tensor of relaxation G jkl being a function of time in the case of linear viscosity.Now, we know, however (e.g., [1,2], and from what was said above), that the assumption of linearity with respect to real life problems is a rather flawed one.In addition, it is almost hopeless to make proper assumption or estimates for the complete tensor function G jklm .Especially in tribological contact situations we will find stress, strain and temperature fields dramatically influencing the tensor function.However, with our inter-atomic based approach from the beginning of the theory section we are in principle able to give a more correct description by the means of our field dependent Young's modulus for the sample: (29) This results in q q jk ll jk s t s t jj q sq sq qq q jk t st j q jk ll jk sq V q q st j In some cases it is more convenient to use a continuous relaxation function, which results in: With the assumption of adiabatic temperature fields following the shear field, such an approach appears to be manageable, but its solution strongly depends on the actual shear distribution (cf.[2]).For the reason of simplicity, we will continue the evaluation with the phenomenological SLS-approach, which should be applicable in situ ations without significant temperature fields [2].
We start as follows: We see that by performing the following substitution underneath the time integrals: s t s t jj q sq sq qq q tt ts st j q sq V q q st j s sq q we are always able to switch back to the more general form.Further on, in order to have maximum flexibility for the construction of a great variety of contact situations combined with sufficient simplicity to keep the model analytical, we introduce the method of time dependent Hertzian load dots as outlined in the Appendix.
The approach has the advantage, that the evaluation of the complete elastic field for the whole body would only be a question of summing up a series of potentially time dependent Hertzian fields.In addition, the solution can easily be extended to layered materials [28].In order to fit the model even better to the real test or application one might also make the constants c i time-dependent.This, of course, requires reevaluation of the Integral in Equation (33) but caused by the distribution structure (properties of the Dirac delta function) in the integrand any such additional time-dependency would be very easy to handle.In addition one could go for extended Hertzian load dots instead of Hertzian ones (e.g., [6]), incorporate tilting, twisting and lateral loads, defects, and intrinsic stresses as done in [12,39,40].Then the method is also applicable to asymmetric scratch-or tribo-like contact situations.

Theory for the Physical Scratch and/or Tribological Test
The governing contact equation and stress distributions for more general loading conditions as occurring during tests like scratch and pin on disc have already been given in the Section "the effective indenter concept".The evaluation of the complex stress and strain fields for these tests is elaborated in the references given there.However, as these evaluations are rather cumbersome and lengthy, the reader is also referred to a software package performing such calculations in an automated and quick manner [28].
It should explicitly be pointed out that knowledge of the complete stress and strain field is essential for a proper failure characterization (scratch) or wear mechanism analysis (tribo).An example is presented in the figures below.There, Figures 4 and 5 present the evolution of the stress field for the von Mises and the normal stress, respectively in the case of a typical scratch test performed on hard coatings on steel substrates (cf.[10][11][12] regarding material and test parameters).Figure 6 shows the typical failure mechanism as being observed in the tests shown in the tests presented in Figures 4 and 5.For more information, the reader is referred to [12].Hence, only such a physical analysis of mechanical contact measurements like instrumented indentations, scratch, and tribo-tests enable one to find out why a surface structure fails under certain loading conditions.These results provide indications on how the investigated coating or surface structure can be improved (structurally be optimized of adapted with respect to the material selection).Now we need to establish the relationship between quasi-static characteristics, like hardness, yield strength, Young's modulus, etc., and highly dynamic characteristics like wear, fretting and other tribological processes and effects.
First of all, it must be understood that within the concept of this approach the process of wear, fretting or a general tribological process is a to be considered as a multi-physical, multi body (asperity) ensemble of contact situations (called load dots, cf.[57,58]) with each having its own geometry, load conditions (including tilting, normal, lateral, twisting loads, plus temperature fields caused by internal or external friction) and-in complex cases with debris-also multiple time scales, meaning various parts of the global tribological contact are running in their own speed.Of course, we will then also have to take temperature effects into account.These are often coupling back into the mechanical properties and time dependent material behavior (cf.Section "The extension of the Oliver and Pharr method to analyze nanoindentation data to layered materials and time dependent mechanical behavior").However, the mathematics of partial differential equations does provide an interesting short cut here.It is due to the similarity of the governing differential equations that by having found the solution for the contact problem also the temperature field or any other "diffusion-like" problem can be considered solved.Thus, the layered solution for any diffusion problem can be constructed from the elastic solution by letting the Poisson's ratio go to infinity, which leads to perfectly non-singular field solution for the diffusion problem in question.
In addition to that, the connection between non-physical parameters (like hardness) towards wear is not the intention of the paper.In fact, there is no such connection-strictly speaking.Something not generic, like hardness, simply cannot, not generally, be extended or applied to a dynamic process like wear.Hardness, after all, is a mixed parameter.There are so many things contributing to hardness that, in general, it is not clear what effects affecting the hardness are also influencing wear (and in what manner), for instance.Thus, yield strength, critical fracture stresses, fracture toughness (energy loss caused by fracture propagation), in short, parameters, which could clearly be connected to certain stress fields and stress components, are much better suited to also interpret and discuss tribological effects.
We are now concentrating on some questions that often pop up in connection with tribological tests and observations.
Where is the role of debris particles?They are just adding up to the complex jumble of contact situations mentioned above.How are strain-rate or even more general time-dependent effects taken into account?This is usually been done by applying time-functional dependencies for the mechanical parameters.However, as shown in [2], in the presence of inhomogenous fields possibly influencing the mechanical parameters the resulting governing equations are becoming highly nonlinear.This can either be treated by the means of piecewise X-parameter models or by solving the nonlinear governing system of equations with time-stress-and strain-dependent material parameters (cf.Section "The extension of the Oliver and Pharr method to analyze nanoindentation data to layered materials and time dependent mechanical behavior").First principle derived basics are given in [2].
How are intrinsic stresses been taken into account?Intrinsic stresses are just adding to the external stresses (and strain) fields and lead to a subsequent shift of onset of inelastic behavior, be it fracture, plastic flow, or phase transition.Therefore it is very important to know and consider intrinsic stresses as accurately as possible.Here, "accurately" explicitly means that intrinsic stresses are usually not homogeneously distributed but also present a complex field adding up to the deformation field coming from external loading situations.
Due to their importance we here refer to a more comprehensive elaboration about the modeling of these stresses especially in connection with external loads (e.g., [41,57,59]).
How can the many faces of wear and tribo effects be incorporated into one theoretical apparatus?
The straight forward answer would be: By decomposition limits extracted from first principle approaches.Hereby, the most general way would be to extract decomposition limits from the first principle approaches as described in [2].These limits have to be compared with the deformation fields obtained in the multiple complex contact model describing our tribo, fretting or wear experiments, where nonlinear effects, as described in [1] and [2], should not be neglected.Even though we will demonstrate how this, in principle, can be done within the procedure-plan below, it is eminent that such a procedure is not necessarily fit for daily or even industrial use.Therefore, we also need a phenomenological description of the physics behind tribo-effects in an as general as possible manner.
At first, we explicitly point out, that wear or any other tribological effect cannot be connected by a simple k d -value or any other scalar value (e.g., Archard's law) to the deformation or stress field.This is by far too simple and not general enough.Instead, in the simplest (linear) case, the tribo-effect is a tensor, coupling with wear-moduli to every deformation field component in a fully covariant manner.Consequently, any tribo-process can be generalized as: It should be made clear, that the sensitivity of a given tribological pairing towards certain components of the deformation field strongly depends on the initial failure mechanisms driving, respectively dominating the process.Thus, we will probably have a higher dependency on the deviatoric stress field (given by the second invariant or the von Mises stress) where the tribology is mainly determined by plastic flow.The examples given in the paper are apparently governed by these effects.In those cases however, where we expect fracture to dominate the tribological process, wear components coupling to the normal stresses or directed shear stresses will increase over the other components of the tribo-or wear tensor.Here, we can easily deduce that mode I fracture would lead to higher components connected with the normal stresses, while mode II and mode III fracture driven tribology requires dominated wear tensors.In some cases also pressure driven phase transitions could be responsible for the observed tribology.Then the tribo-or wear tensor will have dominant terms coupling to the first invariant of the stress field, because this gives the pressure field multiplied by three.
Here we will concentrate mainly on linear dependencies [first line in Equation ( 34)], where we used the following denotations: k xx ijkl -tensors are tensors coupling to the various field values or tensors like the stress σ kl , strain kl , displacement-vector u i or scalar values S n , like free or distortion energy strain work etc.The symbol δ ij is the Kroenecker symbol.In some cases, like simple wear tests for instance, it should be sufficient to consider only the stresses: Wherefrom the scalar wear-depth h w has to be evaluated via: with n i denoting the surface normal unit vector.
We have to point out explicitly, that all the stresses, strains, vectors, energies etc. as being used in Equation (34) are field values and not just numbers.Thus, they are to be understood as functions of the coordinates.However, in the end and in order to have an apparatus one can compare with real experiments one prefers to obtain numbers for the tribo-effect or wear instead of complicated functions.Thus, for many applications the components of the k-tensors must even have operational character.For example, if considering the maximum of a certain stress it might be of great importance where exactly that maximum lays within the material.Let us assume a blunt indenter contact situation and a tribo-effect being dominated by fatigue caused by the total shear.In such a case, we will always find the von Mises stress maximum underneath the surface.It seems reasonable in this situation to assume that the amount of destruction respectively the mass of worn off material is not just connected with the absolute value of the von Mises stress maximum σ vM , but also its distance z to the surface respectively the contact centre r .One could formulate this as: with λ k being yet another parameter in addition to the linear k dvM characterizing the wear performance in dependency of the von Mises stress distribution.We point out that, in using the local load acting on a certain contacted surface point, we have the following local wear depth: Using the quantum mechanical marking to note operators we could use the following generalization for an operational form of (34): For the wear-example given above we can formulate this as: As we can see, the Archard's law given with a scalar wear coefficient k d by the simple relation 33 is nothing but a rather dramatic simplification of Equation ( 35) being possible wherever either stress is dominated by its normal component in the surface normal direction (here, we named it σ 33 ) or where the coefficient tensor k ijkl is zero except for those components coupling to the normal surface stress σ 33 , which would then read: We can deduce now that for complex contact conditions, where the stress tensor is fully set and no component is dominant against all the others, one should be rather careful with the assumption of having the wear-tensor being of the most simple, Archard's-law-like kind One should also take the other stresses into account and investigate their possible influence regarding the resulting global wear, which, in this case, has to be taken as the sum over all stress-components in connection with the wear-tensor:   Here, we have made use of the symmetry of the stress tensor, also requiring a symmetric wear-tensor.We also point out, that in the general law as given above in Equation (35) the hydrostatic (sphere) and deviatoric stress parts are distinguishable.Such a simplified law might read: _ 11 33 or even simpler, and further scalarized: with σ vM , σ H , denoting the von Mises and the hydrostatic stress, respectively.

Experimental Section
Within a variety of examples the procedure or parts of it have already been applied [7][8][9][10][11][12][13][14].Here now, we want to deepen the analysis and discussion by the means of these previous examples and try to improve the simulation by taking nonlinear and temperature effects into account.An illustrative elaboration of the procedure is been given in [60].
It was shown in [13] and [14] that the use of an operator-like wear or tribo-law as introduced in the section above brought some improvement regarding the comparison with experimental data.It was also shown and elaborated that the analytical method as outlined here has the big advantage of delivering quasi "in-sight-views" about the deformation fields, strains, stresses, and so on, of the surfaces subjected to such tribology tests.Due to the completely analytical character of the tools used the calculation usually is extremely fast and robust.This way weak spots within the material composition can be found more easily and efficiently.Subsequently, also optimizing surface systems towards a better tribo-performance could be done quicker and still more accurately.
We have demonstrated this on the wear track example presented in [13] and [14], where an oscillating wear test with a spherical indenter had been performed.The indenter was sliding forth and back along a certain track thereby wearing off material.In [60] (Figures 11-19) it was shown how the theoretical models mirrors the gradual "digging" in of the indenter with increasing cycle number and the subsequent evolution of the stress field for a variety of stress components.The evaluation was performed with the software FilmDoctor [28].
Still, there had been deviations from the experiment which could not been explained in a satisfactory manner.However, it was already mentioned (cf.[13,14] and also see [2]) that nonlinear or temperature effects could play a role.Especially with respect to temperature one would expect to see some influence as the tests performed were reciprocating wear tests where a spherical indenter is performing an oscillatory sliding movement over the, in some cases coated, sample surface.Here, the friction will automatically lead to a temperature increase.With such an increase of temperature, we expect a change of the mechanical parameters, especially Young's modulus and yield strength.
Simply assuming such an influence and incorporating it into the wear simulation as described in [13] and [14] by the means of a Young's modulus function of the kind given in Equation ( 29) thereby assuming a temperature dependency of the relaxation time or viscosity like: We hope to result into a better fit of our wear or tribo model to the experimental data.Using a similar dependency for the yield strength Y as we did for the Young's modulus and adapting our wear-law (40) such that the wear should increase with decreasing Y, we set: ˆe () Here, we need to pay attention to the fact, that with having a temperature dependent Young's modulus by the means of Equations ( 29) and (45) we also result in a temperature dependency for the von Mises stress maximum and its position.This makes the evaluation of Equation ( 46) relatively complicated and we therefore restrict ourselves here to the task of finding a temperature function T(n) with n denoting the number of cycles or oscillations of the indenter which would fit the real wear data significantly better than an assumed constant temperature.For both examples, we had at hand [13,14] and we received significantly better fit results in this way.Interestingly, the best function T(n) to use was of the kind T(n) = c 0 × tanh(c 1 × n) which simply coincides with a temperature rise at the beginning of the test and reaching a plateau value after only a view cycles.This surprisingly simple result does not only have some beauty with respect to its low number of free parameters, but also can easily be understood by the interpretation of the indenter producing friction and leading to a temperature increase until a state of equilibrium is reached where the temperature in the vicinity of the contact does not increase any further.The results are shown in Figures 7 and 8.The data have been reused from [14].
It should be pointed out that, especially in scratch and wear tests, often, effects like work hardening, grain refinement, intrinsic stress, and general defect accumulation lead to significant material parameter changes during the test.In principle the model allows to account for such effects simply by making the material parameters which are affected by such effects dependent on the evolution of the test.Such a dependency could either be steered via time, wear cycle or any other evolutional parameter counting the test progress.Within the examples given in Figures 7 and 8 we only found enough experimental significance for the temperature coupling into the material properties leading to a parameter dependency for the Young's modulus and the yield strength.As elaborated above, this dependency was accounted for by using the wear cycle number n as evolutional parameter.

Conclusions
By combining the global incremental wear model with the concept of the effective indenter, the extended Hertzian approach and a layered half space solution for contact problems with rather arbitrary combinations of normal, lateral and tilting loads a general, quick and powerful wear model has been created.
Inversion of the model makes it fit for parameter identification and optimization problems.Further, adding first principle models, like effective interatomic interaction potentials allows a deeper understanding of those failure mechanisms being responsible for wear results observable in practical tests.
Temperature effects had been taken into account by allowing the temperature to couple into the mechanical properties Young's modulus and yield strength.
Possible subsequent time dependent material behavior leads to a relatively complex inhomogeneous Young's modulus distribution which must be taken into account if one aims for a better simulation of wear experiments.
The analytical method as outlined here has the advantage of delivering quasi "in-sight-views" about the deformation fields, strains, stresses, and so on, of the surfaces subjected to such tribology tests.Due to the completely analytical character of the tools used the calculation usually is extremely fast and robust.This way weak spots within the possibly layered material composition can be found more easily and efficiently.Subsequently, also optimizing surface systems towards a better tribo-performance could be done quicker and still more accurately.

Arbitrary Hertzian Load Dots
By introducing Hertzian load-dots as: we obtain the rather simple expression for the stress tensor with index i: Now we derive the situation for a stress driven deformation state.Again applying the 3-parameter model together with our simple load-dot approach we have:   ( )e d 1, 2,... ; , 1, 2,3 where the stresses  are now to be understood as Generalized with respect to field dependent mechanical properties and sticking to a 3-parameter potential interaction the Equation (A7) would read: Extension to discrete or continuous retardation gives:   1 e d () For reasons of simplicity however, we will here proceed with the gradient free simple SLS-model.Together with the normalization condition and the total normal load p = p(t) one obtains Both cases are relatively easy whenever there will be no change of the contact area during the evolution of the deformation time being considered.In most practical applications however, this is not the case.Concentrating on contact experiments we can often assume that they are load controlled.However, due to the complex mixed stress states underneath the contact zone, there are also strain effects involved and it is not clear yet how to take these into account.Therefore Equations (A6)-(A11) and (A16) are applicable but need to be extended with respect to strain driving effects (e.g., [2], Appendix).There are several options how to make the contact area evolvable during the time slot being considered.One could either change the number of load dots being taken into account respectively contributing to the whole contact or make the contact radii a i time dependent.As the latter way immediately leads to much more complicated integrals we are going for the first option not by the apparent simple way of changing the total load-dot number N but by introducing Heaviside distributions to the γ i (t) (or respectively the  i (t) in the strain drive case) as follows: and δ giving the Dirac delta distribution.With this approach the integration in Equation (A7) can be performed easily as Equation (A19).For many cases in contact applications we are facing mixed stress-strain driven conditions requiring a combination of Equations (A4) and (A7).Choosing a sufficiently high number of load dots N one can approximate any displacement (time dependent effective indenter shape) of symmetry of revolution.By placing the load dots in a non concentric manner also asymmetric contact situation can be described. :

Figure 2 .
Figure2.The effective indenter concept (see:[4]) transferring the theoretically difficult problem of a well defined sharp indenter on (a) an elasto-plastically deformed surface with complex shape to (b) an effective indenter on a flat surface.

Figure 3 .
Figure 3. Formulating the contact equation h = w + z + w' + z' with respect to partial plastic deformation z of the sample surface.

Figure 4 .Figure 5 .
Figure 4.The evolution of von-Mises stress during the scratch test shown at three measurement positions: (a) at the beginning of the scratch test, (c) in the moment of L C failure, and (b) in between.The black cross hairs indicate the location of the maximum.

Figure 6 .
Figure 6.Illustrative scheme of the failure mechanism (a) and an optical graph of the post-scratch surface (b) in which the corresponding L C position is marked by the red dashed line [12].

Figure 7 .
Figure 7.Comparison of various wear laws with the experimental data (cf.text in Section 3).The data have been reused from[13].

Figure 8 .
Figure 8.Comparison of various wear laws with the experimental data (cf.text in Section 3).The data have been reused from[14].
driven case, while one derives in the stress driven situation with 