contact mechanics and elements of tribology · contact mechanics and elements of tribology lecture...
TRANSCRIPT
![Page 1: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/1.jpg)
Contact Mechanics and Elements of Tribology
Lecture 4.Surface roughness
Vladislav A. Yastrebov
MINES ParisTech, PSL Research University, Centre des Matériaux, CNRS UMR 7633,Evry, France
@ Centre des MatériauxFebruary 20, 2018
![Page 2: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/2.jpg)
Outline •
Introduction
Measurement techniques
Classifications
Main characteristics
PDF and PSD
Random process model of roughness
Computational roughness models
Reading
V.A. Yastrebov Lecture 4 2/100
![Page 3: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/3.jpg)
Roughness •Natural and industrial
surfaces are rough:
processingpolishingcoatingmicrostructuresurface energydeformationagingenvironment
Fig. Examples of rough surfaces
V.A. Yastrebov Lecture 4 3/100
![Page 4: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/4.jpg)
Roughness •Natural and industrial
surfaces are rough:
processingpolishingcoatingmicrostructuresurface energydeformationagingenvironment
Fig. Persistent slip marks [1]
Fig. Rumpling (thermal cycling induced roughness in air)[2]
Fig. Epitaxial surface growth [3,4]
[1] J.Polák, J. Man & K. Orbtlík, Int J Fatigue 25 (2003)[2] V.K. Tolpygo, D.R. Clarke, Acta Mat 52 (2004)[3] M. Einax, W. Dieterich, P. Maass, Rev Mod Phys 85 (2013)[4] J.R. Arthur, Surf Sci 500 (2002)
V.A. Yastrebov Lecture 4 4/100
![Page 5: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/5.jpg)
Roughness •Natural and industrial
surfaces are rough:
processingpolishingcoatingmicrostructuresurface energydeformationagingenvironment
Roughness affects:
stress-strain statedry frictionwearadhesionfluid flowsealingenergy transfer
Fig. True contact area and stress fluctuations
V.A. Yastrebov Lecture 4 5/100
![Page 6: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/6.jpg)
Roughness •Natural and industrial
surfaces are rough:
processingpolishingcoatingmicrostructuresurface energydeformationagingenvironment
Roughness affects:
stress-strain statedry frictionwearadhesionfluid flowsealingenergy transfer
Fig. True contact area and stress fluctuations
V.A. Yastrebov Lecture 4 6/100
![Page 7: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/7.jpg)
Roughness •Natural and industrial
surfaces are rough:
processingpolishingcoatingmicrostructuresurface energydeformationagingenvironment
Roughness affects:
stress-strain statedry frictionwearadhesionfluid flowsealingenergy transfer
Fig. Numerical simulation of airflow around a (dimpled) golf ball[5]
[5] C.E. Smith, PhD thesis (2011)
V.A. Yastrebov Lecture 4 7/100
![Page 8: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/8.jpg)
Roughness •Natural and industrial
surfaces are rough:
processingpolishingcoatingmicrostructuresurface energydeformationagingenvironment
Roughness affects:
stress-strain statedry frictionwearadhesionfluid flowsealingenergy transfer
Fig. Fluid passage through free volume between rough surfaces
Fig. Heat transfer between rough surfaces(asperity-based model)
V.A. Yastrebov Lecture 4 8/100
![Page 9: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/9.jpg)
Surface metrology techniques •
Stylus measurements
• Mechanica contact of a tip with surface• Force ≥ 3µg, tip radius ≥ 50 nm• Mainly for profile measurements y(x)
Optical measurements
• Confocal (laser scanning) microscopy- highest lateral resolution
• Interferometry (WLI):- highest vertical resolution- 10 to 100 times faster than CM
• Scanning Electronic Microscopy(SEM):- in secondary electron emission- electrons penetrate in the matter→roughness smoothing- conducting materials
Nano-contact measurements
• Atomic Force Microscopy (AFM)roughness + adhesive and elastic properties
• Scanning Tunneling Microscope (STM)
Stylus profilometer
Valley measurement error
V.A. Yastrebov Lecture 4 9/100
![Page 10: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/10.jpg)
Surface metrology techniques •
Stylus measurements
• Mechanica contact of a tip with surface• Force ≥ 3µg, tip radius ≥ 50 nm• Mainly for profile measurements y(x)
Optical measurements
• Confocal (laser scanning) microscopy- highest lateral resolution
• Interferometry (WLI):- highest vertical resolution- 10 to 100 times faster than CM
• Scanning Electronic Microscopy(SEM):- in secondary electron emission- electrons penetrate in the matter→roughness smoothing- conducting materials
Nano-contact measurements
• Atomic Force Microscopy (AFM)roughness + adhesive and elastic properties
• Scanning Tunneling Microscope (STM)
Stylus profilometer
Modern stylus profilometerwww.bruker.com
Roughness measurements (∆z ≈ 30 µm)www.icryst.com
V.A. Yastrebov Lecture 4 10/100
![Page 11: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/11.jpg)
Surface metrology techniques •
Stylus measurements
• Mechanica contact of a tip with surface• Force ≥ 3µg, tip radius ≥ 50 nm• Mainly for profile measurements y(x)
Optical measurements
• Confocal (laser scanning) microscopy- highest lateral resolution
• Interferometry (WLI):- highest vertical resolution- 10 to 100 times faster than CM
• Scanning Electronic Microscopy(SEM):- in secondary electron emission- electrons penetrate in the matter→roughness smoothing- conducting materials
Nano-contact measurements
• Atomic Force Microscopy (AFM)roughness + adhesive and elastic properties
• Scanning Tunneling Microscope (STM)
Confocal microscopy
Principle of confocal microscopyadapted from www.wikipedia.org
V.A. Yastrebov Lecture 4 11/100
![Page 12: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/12.jpg)
Surface metrology techniques •
Stylus measurements
• Mechanica contact of a tip with surface• Force ≥ 3µg, tip radius ≥ 50 nm• Mainly for profile measurements y(x)
Optical measurements
• Confocal (laser scanning) microscopy- highest lateral resolution
• Interferometry (WLI):- highest vertical resolution- 10 to 100 times faster than CM
• Scanning Electronic Microscopy(SEM):- in secondary electron emission- electrons penetrate in the matter→roughness smoothing- conducting materials
Nano-contact measurements
• Atomic Force Microscopy (AFM)roughness + adhesive and elastic properties
• Scanning Tunneling Microscope (STM)
Confocal microscopy
1 euro surface www.wikipedia.org
Stainless steel machined with micro-electricdischarge www.laserfocusworld.org
V.A. Yastrebov Lecture 4 12/100
![Page 13: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/13.jpg)
Surface metrology techniques •
Stylus measurements
• Mechanica contact of a tip with surface• Force ≥ 3µg, tip radius ≥ 50 nm• Mainly for profile measurements y(x)
Optical measurements
• Confocal (laser scanning) microscopy- highest lateral resolution
• Interferometry (WLI):- highest vertical resolution- 10 to 100 times faster than CM
• Scanning Electronic Microscopy(SEM):- in secondary electron emission- electrons penetrate in the matter→roughness smoothing- conducting materials
Nano-contact measurements
• Atomic Force Microscopy (AFM)roughness + adhesive and elastic properties
• Scanning Tunneling Microscope (STM)
White Light
Interferometry
Diamond-turned optics www.zygo.com
US quarter surface www.zygo.com
V.A. Yastrebov Lecture 4 13/100
![Page 14: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/14.jpg)
Surface metrology techniques •
Stylus measurements
• Mechanica contact of a tip with surface• Force ≥ 3µg, tip radius ≥ 50 nm• Mainly for profile measurements y(x)
Optical measurements
• Confocal (laser scanning) microscopy- highest lateral resolution
• Interferometry (WLI):- highest vertical resolution- 10 to 100 times faster than CM
• Scanning Electronic Microscopy(SEM):- in secondary electron emission- electrons penetrate in the matter→roughness smoothing- conducting materials
Nano-contact measurements
• Atomic Force Microscopy (AFM)roughness + adhesive and elastic properties
• Scanning Tunneling Microscope (STM)
AFM
Modern AFMwww.bruker.com
Roughness and elastic moduli (color) ofpolymer blend
www.bruker.com
V.A. Yastrebov Lecture 4 14/100
![Page 15: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/15.jpg)
Surface metrology techniques •
Stylus measurements
• Mechanica contact of a tip with surface• Force ≥ 3µg, tip radius ≥ 50 nm• Mainly for profile measurements y(x)
Optical measurements
• Confocal (laser scanning) microscopy- highest lateral resolution
• Interferometry (WLI):- highest vertical resolution- 10 to 100 times faster than CM
• Scanning Electronic Microscopy(SEM):- in secondary electron emission- electrons penetrate in the matter→roughness smoothing- conducting materials
Nano-contact measurements
• Atomic Force Microscopy (AFM)roughness + adhesive and elastic properties
• Scanning Tunneling Microscope (STM)
STM
Fig. Center for NanoScience logo imprintedat atomic scalewww.cens.de
Atomic steps on platinum surface(500×500 nm)
www.icryst.com
V.A. Yastrebov Lecture 4 15/100
![Page 16: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/16.jpg)
Roughness: classification •
Fig. Roughness classification according to Nayak[1]
[1] Nayak, J. Lub. Tech. (ASME) 93:398 (1971)
V.A. Yastrebov Lecture 4 16/100
![Page 17: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/17.jpg)
Roughness: classification •
Fig. Roughness classification according to Nayak[1]
[1] Nayak, J. Lub. Tech. (ASME) 93:398 (1971)
V.A. Yastrebov Lecture 4 17/100
![Page 18: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/18.jpg)
Roughness and geometry/form •Roughness vs geometry of surfaces
Sometimes macroscopic geometry is substracted (filtered out):form→error of form→wavyness→roughness
Non-trivial to remove macroscopic shape
Most roughness measurement tools enable shape removal
Fig. Circular metalic seal with turned copper surface[1]
[1] F.P. Rafòls, Licentiate Thesis, LTU 2016.
V.A. Yastrebov Lecture 4 18/100
![Page 19: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/19.jpg)
Roughness and geometry/form •Roughness vs geometry of surfaces
Sometimes macroscopic geometry is substracted (filtered out):form→error of form→wavyness→roughness
Non-trivial to remove macroscopic shape
Most roughness measurement tools enable shape removal
Fig. (left) impact crater, (right) shape is filtered out
V.A. Yastrebov Lecture 4 19/100
![Page 20: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/20.jpg)
Roughness and geometry/form •Roughness vs geometry of surfaces
Sometimes macroscopic geometry is substracted (filtered out):form→error of form→wavyness→roughness
Non-trivial to remove macroscopic shape
Most roughness measurement tools enable shape removal
Fig. (left) spherical indenter, (right) z = a(x2 + y2) shape is substracted
V.A. Yastrebov Lecture 4 20/100
![Page 21: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/21.jpg)
Main characteristics •Integral quantities
Average of absolute values [l.u.] (profile - Ra, surface - Sa)
Sa =1
A
∫
A
|z(x, y) − z| dA, Sa =1
N2
N∑
i=1
N∑
j=1
|zij − z|
Standard deviation of height [l.u.] (σ or Rq for profile, Sq for surface)
σ =
√
√
1
A
∫
A
(z(x, y) − z)2 dA, σ =1
N
√
√
√
N∑
i=1
N∑
j=1
(zij − z)2
Maximal valley depth Rv,Sv, maximal peak height Rp,Sp [l.u.]very sensitive to sample area
Skewness [adim] (γ1 or Rsk,Ssk)
γ1 =1
Aσ3
∫
A
(z(x, y) − z)3 dA, γ1 =1
N2σ3
N∑
i=1
N∑
j=1
(zij − z)3
Kurtosis [adim] (κ or Rku,Sku)
κ =1
Aσ4
∫
A
(z(x, y) − z)4 dA, κ =1
N2σ4
N∑
i=1
N∑
j=1
(zij − z)4
V.A. Yastrebov Lecture 4 21/100
![Page 22: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/22.jpg)
Main characteristics II •Integral quantities II
Average of absolute value of gradient (slope) [adim] (profile - Rdq,surface - Sdq)
Sda = 〈∣
∣
∣∇z(x, y) − ∇z∣
∣
∣〉 = 1
A
∫
A
∣
∣
∣∇z(x, y) − ∇z∣
∣
∣ dA
Sda =1
N2
N∑
i=1
N∑
j=1
∣
∣
∣
∣
∣
∣
∣
zi+1,j − zi,j − ∆zx
∆x
∣
∣
∣
∣
∣
∣
∣
+
∣
∣
∣
∣
∣
∣
∣
zi,j+1 − zi,j − ∆zy
∆y
∣
∣
∣
∣
∣
∣
∣
Standard deviation of gradient (slope) [adim] (profile - Rdq, surface - Sdq)
Sdq =
√
〈∣
∣
∣∇z(x, y) − ∇z∣
∣
∣
2〉 =√
√
1
A
∫
A
∣
∣
∣∇z(x, y) − ∇z∣
∣
∣
2dA
Sdq =
√
√
√
√
1
N2
N∑
i=1
N∑
j=1
∣
∣
∣
∣
∣
∣
∣
zi+1,j − zi,j − ∆zx
∆x
∣
∣
∣
∣
∣
∣
∣
2
+
∣
∣
∣
∣
∣
∣
∣
zi,j+1 − zi,j − ∆zy
∆y
∣
∣
∣
∣
∣
∣
∣
2
Often in integrated slope measurements a smoothing filter is used, forexample, according to ASME B46.1 standard
∂z
∂x≈ 1
60∆x(zi+3,j − 9zi+2,j + 45zi+1,j − 45zi−1,j + 9zi−2,j − zi−3,j)
V.A. Yastrebov Lecture 4 22/100
![Page 23: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/23.jpg)
Main characteristics: probability density •
Probability density of heights P(z)
Properties and moments
1 =
∞∫
−∞
P(z) dz, z =
∞∫
−∞
zP(z) dz, σ =
√
√
√
√
√
∞∫
−∞
(z − z)2P(z) dz
µq =
∞∫
−∞
zqP(z) dz then µ0 = 1, µ1 = z, µ2 = σ2 + z2
Link to skewness
V.A. Yastrebov Lecture 4 23/100
![Page 24: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/24.jpg)
Main characteristics: probability density II •Distribution examples
Normal (Gaussian): P(x) =1
σ√
2πexp
[
− (x − µ1)2
2σ2
]
, x ∈ R
adapted from wikipedia.org
V.A. Yastrebov Lecture 4 24/100
![Page 25: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/25.jpg)
Main characteristics: probability density II •Distribution examples
Lognormal: P(x) =1
xσ√
2πexp
[
− (log(x) − µ1)2
2σ2
]
x ∈ R+
adapted from wikipedia.org
V.A. Yastrebov Lecture 4 25/100
![Page 26: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/26.jpg)
Main characteristics: probability density II •Distribution examples
Weibull: P(x) =k
λ
(
x
λ
)k−1
exp−(x/λ)k x ∈ R+
adapted from wikipedia.org
V.A. Yastrebov Lecture 4 26/100
![Page 27: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/27.jpg)
Main characteristics: comments on PDF •Real rough surfaces
Turning, scratching, shaping changes macroscopic distribution Pmacro
but might keep microscopic distribution intact Pmicro
Wear, polishing, flattening results in removal of the right distributiontail P(z > z0)→ 0: negative curtosis
Macro- and microscopic roughness (left: diamond-turned surface, right: cross-hatched surface)Images from www.zygo.com used
V.A. Yastrebov Lecture 4 27/100
![Page 28: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/28.jpg)
Main characteristics: comments on PDF •Real rough surfaces
Turning, scratching, shaping changes macroscopic distribution Pmacro
but might keep microscopic distribution intact Pmicro
Wear, polishing, flattening results in removal of the right distributiontail P(z > z0)→ 0: negative curtosis
Turned surface topography and a sketch of height PDFPérez-Ràfols, Larsson, Almqvist, Tribol Int 94 (2016)
V.A. Yastrebov Lecture 4 28/100
![Page 29: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/29.jpg)
Main characteristics: comments on PDF •Real rough surfaces
Turning, scratching, shaping changes macroscopic distribution Pmacro
but might keep microscopic distribution intact Pmicro
Wear, polishing, flattening results in removal of the right distributiontail P(z > z0)→ 0: negative curtosis
Initial surface with Gaussian PDF and surface after polishing
V.A. Yastrebov Lecture 4 29/100
![Page 30: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/30.jpg)
Main characteristics: comments on PDF •Real rough surfaces
Turning, scratching, shaping changes macroscopic distribution Pmacro
but might keep microscopic distribution intact Pmicro
Wear, polishing, flattening results in removal of the right distributiontail P(z > z0)→ 0: negative curtosis
Very fresh asphalt concrete Normal asphalt concrete
V.A. Yastrebov Lecture 4 30/100
![Page 31: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/31.jpg)
Main characteristics: comments on PDF •Real rough surfaces
Turning, scratching, shaping changes macroscopic distribution Pmacro
but might keep microscopic distribution intact Pmicro
Wear, polishing, flattening results in removal of the right distributiontail P(z > z0)→ 0: negative curtosis
Very fresh asphalt concrete Normal asphalt concrete . . . with a bolt ,
V.A. Yastrebov Lecture 4 31/100
![Page 32: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/32.jpg)
Main characteristics: comments on PDF •Real rough surfaces
Turning, scratching, shaping changes macroscopic distribution Pmacro
but might keep microscopic distribution intact Pmicro
Wear, polishing, flattening results in removal of the right distributiontail P(z > z0)→ 0: negative curtosis
Old asphalt concrete with worn out bitumen
V.A. Yastrebov Lecture 4 32/100
![Page 33: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/33.jpg)
Main characteristics: comments on PDF •Real rough surfaces
Turning, scratching, shaping changes macroscopic distribution Pmacro
but might keep microscopic distribution intact Pmicro
Wear, polishing, flattening results in removal of the right distributiontail P(z > z0)→ 0: negative curtosis
Atlas of machined surfaces (with height distributions)
V.A. Yastrebov Lecture 4 33/100
![Page 34: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/34.jpg)
Autocorrelation function •
Continuous autocorrelation function
R(∆x,∆y) = limL→∞
1
L2
L∫
0
L∫
0
z(x + ∆x, y + ∆y)z(x, y) dxdy
Discrete autocorrelation function for a surface N ×N
R(∆x,∆y) =1
N2
N−1∑
i=0
N−1∑
j=0
z(x + ∆x, y + ∆y)z(x, y)
V.A. Yastrebov Lecture 4 34/100
![Page 35: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/35.jpg)
Power spectral density (PSD) •
Recall: Fourier Transform: f (k) =
∞∫
−∞
f (x) exp(−2πikx) dx
Recall: Discrete Fourier Transform: fk =
N−1∑
n=0
xn exp(−2πikn/N)
where x is the spatial coordinate, k = 2π/λ is the wavenumber and λ isthe wavelength.
PSD is the Fourier Transform of R
Φ(kx, ky) ≡ R(kx, ky) = FFT[
z(x + ∆x, y + ∆y) ∗ z(x, y)]
Using convolution theorem
Φ(kx, ky) = z(kx, ky) z∗(kx, ky) = z2(kx, ky)
Interpretation: energy distribution by frequencies
Usage: signal analysis, seismology, microstructure characterization,roughness.
V.A. Yastrebov Lecture 4 35/100
![Page 36: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/36.jpg)
Spectral moments •Spectral moment mpq, p, q ∈N:
mpq =
∞"
−∞
kpxk
qyΦ(kx, ky) dkxdky
mpq =
[
2π
L
]p+q N−1∑
i=0
N−1∑
j=0
ipjq Φ(2πi/L, 2πj/L)
Generalized spectral moment mpq, p, q ∈ R+
For isotropic surface: m2 = m20 = m02, m4 = 3m22 = m40 = m04
Averaging:
m2 =m20 +m02
2, m4 =
m40 + 3m22 +m04
3
Physical meaning:
Height variance1: m0 =⟨
(z − 〈z〉)2⟩
Gradient variance: 2m2 =⟨
(∇z − 〈∇z〉)2⟩
Curvature variance: m4 =⟨
(∇ · ∇z − 〈∇ · ∇z〉)2⟩
1Variance is a squared standard deviation
V.A. Yastrebov Lecture 4 36/100
![Page 37: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/37.jpg)
Nayak parameter •
Recall:
Height variance2: m0 =⟨
(z − 〈z〉)2⟩
Gradient variance: 2m2 =⟨
(∇z − 〈∇z〉)2⟩
Curvature variance: m4 =⟨
(∇ · ∇z − 〈∇ · ∇z〉)2⟩
Nayak parameter (interpretation: breadth of the spectrum):
α =m0m4
m22
, α ∈ (1.5,∞]
α < ∞: the spectrum is bounded in terms of the mean curvature.α→∞: all harmonics affect the curvature, the mean curvature tends toinfinity.
2Variance is a squared standard deviation
V.A. Yastrebov Lecture 4 37/100
![Page 38: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/38.jpg)
Comment on spectral density description •
PSD cannot distinguish deepvalleys and high peaks[1], i.e.Φ(z) = Φ(−z)
From the mechanical point ofview two surfaces with z(x, y)and −z(x, y) can be different
Wear, polishing, flatteningresults in removal of highasperities (results in negativecurtosis)
[1] Borodich, Pepelyshev, Savencu. Tribol Int 103 (2016).
[2] Yastrebov, Durand, Proudhon, Cailletaud. CR Mech
339 (2011).
V.A. Yastrebov Lecture 4 38/100
![Page 39: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/39.jpg)
Comment on spectral density description •
PSD cannot distinguish deepvalleys and high peaks[1], i.e.Φ(z) = Φ(−z)
From the mechanical point ofview two surfaces with z(x, y)and −z(x, y) can be different
Wear, polishing, flatteningresults in removal of highasperities (results in negativecurtosis)
However, if at small scales theheight probability isGaussian[1,2], the PSD descriptioncan be used at this scale
[1] Borodich, Pepelyshev, Savencu. Tribol Int 103 (2016).
[2] Yastrebov, Durand, Proudhon, Cailletaud. CR Mech
339 (2011). Asphalt is non-Gaussian at macro-scale,but Gaussian at small (aggregates’ tip) scale.
V.A. Yastrebov Lecture 4 39/100
![Page 40: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/40.jpg)
Summary •
Fractal (self-affine) roughness
Power spectral density (PSD)
Φ(k) ∼ k−2(H+1)
k is a wavenumber,H is the Hurst exponent.
Isotropic/anisotropic surfaces
Gaussian/non-Gaussian heightdistribution P(h)
Fig. 3D and radial power spectral densities
V.A. Yastrebov Lecture 4 40/100
![Page 41: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/41.jpg)
Summary •
Fractal (self-affine) roughness
Power spectral density (PSD)
Φ(k) ∼ k−2(H+1)
k is a wavenumber,H is the Hurst exponent.
Isotropic/anisotropic surfaces
Gaussian/non-Gaussian heightdistribution P(h)
Fig. 3D and radial power spectral densities
V.A. Yastrebov Lecture 4 41/100
![Page 42: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/42.jpg)
Summary •
Fractal (self-affine) roughness
Power spectral density (PSD)
Φ(k) ∼ k−2(H+1)
k is a wavenumber,H is the Hurst exponent.
Isotropic/anisotropic surfaces
Gaussian/non-Gaussian heightdistribution P(h)
Fig. Power spectral density, measurements
[1] Majumdar, Tien, Wear 136 (1990)[2] Schmittbuhl, Jørgen Måløy, Phys. Rev. Lett. 78 (1997)[3] Vallet, Lasseux, Sainsot, Zahouani, Tribol. Int. 42 (2009)
V.A. Yastrebov Lecture 4 42/100
![Page 43: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/43.jpg)
Summary •
Fractal (self-affine) roughness
Power spectral density (PSD)
Φ(k) ∼ k−2(H+1)
k is a wavenumber,H is the Hurst exponent.
Isotropic/anisotropic surfaces
Gaussian/non-Gaussian heightdistribution P(h)
Fig. Power spectral density, geological scales
Adapted from[4] Renard, Candela, Bouchaud, Geophys. Res. Lett. 40 (2013)
V.A. Yastrebov Lecture 4 43/100
![Page 44: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/44.jpg)
Summary •
Fractal (self-affine) roughness
Power spectral density (PSD)
Φ(k) ∼ k−2(H+1)
k is a wavenumber,H is the Hurst exponent.
Isotropic/anisotropic surfaces
Gaussian/non-Gaussian heightdistribution P(h)
Fig. Height distribution P(h)
Fig. Height distribution of a polished metalsurface
V.A. Yastrebov Lecture 4 44/100
![Page 45: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/45.jpg)
Summary •
Fractal (self-affine) roughness
Power spectral density (PSD)
Φ(k) ∼ k−2(H+1)
k is a wavenumber,H is the Hurst exponent.
Isotropic/anisotropic surfaces
Gaussian/non-Gaussian heightdistribution P(h)
Characteristics:
•√
〈z2〉 - rms heights
•√
〈|∇z|2〉 - rms slope (surface gradient)
• α = m00m40/m220
- breadthof the spectrum (Nayak’s parameter[B]),
spectral moments mpq =
∞!−∞
kpxk
qyΦ(kx , ky) dkxdky
Random process theory
[A] Longuet-Higgins, Philos. Trans. R. Soc. A 250:157 (1957)[B] Nayak, J. Lub. Tech. (ASME) 93:398 (1973)[C] Greenwood, Wear 261: 191 (2006)
Fig. Height distribution P(h)
Fig. Height distribution of a polished metalsurface
V.A. Yastrebov Lecture 4 45/100
![Page 46: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/46.jpg)
Summary •
Fractal (self-affine) roughness
Power spectral density (PSD)
Φ(k) ∼ k−2(H+1)
k is a wavenumber,H is the Hurst exponent.
Isotropic/anisotropic surfaces
Gaussian/non-Gaussian heightdistribution P(h)
Characteristics:
•√
〈z2〉 - rms heights
•√
〈|∇z|2〉 - rms slope (surface gradient)
• α = m00m40/m220
- breadthof the spectrum (Nayak’s parameter[B]),
spectral moments mpq =
∞!−∞
kpxk
qyΦ(kx , ky) dkxdky
Random process theory
[A] Longuet-Higgins, Philos. Trans. R. Soc. A 250:157 (1957)[B] Nayak, J. Lub. Tech. (ASME) 93:398 (1973)[C] Greenwood, Wear 261: 191 (2006) Distribution of asperity curvatures
V.A. Yastrebov Lecture 4 46/100
![Page 47: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/47.jpg)
Summary •
Fig. Example of a rough surface for H = 0.3 Fig. Example of a rough surface for H = 0.8
Recall: the Hurst exponent H and the fractal dimension D in 2D space are interconnected via D = 3 −H
V.A. Yastrebov Lecture 4 47/100
![Page 48: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/48.jpg)
Fractals •
Flight over a rough surface
V.A. Yastrebov Lecture 4 48/100
![Page 49: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/49.jpg)
Fractals •
Romanesco broccoli www.fourmilab.ch
V.A. Yastrebov Lecture 4 49/100
![Page 50: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/50.jpg)
Fractals •
Mandelbrot set (not a fractal)
Recursive functionzi+1 = z2
i + z, z ∈ C
V.A. Yastrebov Lecture 4 50/100
![Page 51: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/51.jpg)
Fractals •
Mandelbrot set (not a fractal)
Recursive functionzi+1 = z2
i + z, z ∈ C
V.A. Yastrebov Lecture 4 51/100
![Page 52: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/52.jpg)
Fractals •
Mandelbrot set (not a fractal)
Recursive functionzi+1 = z2
i + z, z ∈ C
V.A. Yastrebov Lecture 4 52/100
![Page 53: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/53.jpg)
Fractals •
Mandelbrot set (not a fractal)
Recursive functionzi+1 = z2
i + z, z ∈ C
V.A. Yastrebov Lecture 4 53/100
![Page 54: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/54.jpg)
Fractals •
Mandelbrot set (not a fractal)
Recursive functionzi+1 = z2
i + z, z ∈ C
V.A. Yastrebov Lecture 4 54/100
![Page 55: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/55.jpg)
Fractals •
Mandelbrot set (not a fractal)
Recursive functionzi+1 = z2
i + z, z ∈ C
V.A. Yastrebov Lecture 4 55/100
![Page 56: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/56.jpg)
Fractals •
Mandelbrot set (not a fractal)
Recursive functionzi+1 = z2
i + z, z ∈ C
V.A. Yastrebov Lecture 4 56/100
![Page 57: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/57.jpg)
Fractals •
Mandelbrot set (Wikipedia)
V.A. Yastrebov Lecture 4 57/100
![Page 58: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/58.jpg)
Fractals •
Element of Mandelbrot set (Wikipedia)
V.A. Yastrebov Lecture 4 58/100
![Page 59: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/59.jpg)
Fractals •
Element of Mandelbrot set (Wikipedia)
V.A. Yastrebov Lecture 4 59/100
![Page 60: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/60.jpg)
Fractals •
Element of Mandelbrot set (Wikipedia)
V.A. Yastrebov Lecture 4 60/100
![Page 61: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/61.jpg)
Fractals •
Element of Mandelbrot set (Wikipedia)
V.A. Yastrebov Lecture 4 61/100
![Page 62: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/62.jpg)
Fractals •
Element of Mandelbrot set (Wikipedia)
V.A. Yastrebov Lecture 4 62/100
![Page 63: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/63.jpg)
Fractals •
Animation
V.A. Yastrebov Lecture 4 63/100
![Page 64: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/64.jpg)
Synthesized rough surfaces: in pictures •
[1] Y. Z. Hu and K. Tonder, Int. J. Machine Tools Manuf. 32, 83 (1992)
V.A. Yastrebov Lecture 4 64/100
![Page 65: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/65.jpg)
Synthesized rough surfaces: in pictures •
[1] Y. Z. Hu and K. Tonder, Int. J. Machine Tools Manuf. 32, 83 (1992)
V.A. Yastrebov Lecture 4 65/100
![Page 66: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/66.jpg)
Synthesized rough surfaces: in pictures •
[1] Y. Z. Hu and K. Tonder, Int. J. Machine Tools Manuf. 32, 83 (1992)
V.A. Yastrebov Lecture 4 66/100
![Page 67: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/67.jpg)
Synthesized rough surfaces: in pictures •
[1] Y. Z. Hu and K. Tonder, Int. J. Machine Tools Manuf. 32, 83 (1992)
V.A. Yastrebov Lecture 4 67/100
![Page 68: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/68.jpg)
Synthesized rough surfaces: in pictures •
[1] Y. Z. Hu and K. Tonder, Int. J. Machine Tools Manuf. 32, 83 (1992)
V.A. Yastrebov Lecture 4 68/100
![Page 69: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/69.jpg)
Synthesized rough surfaces: in equations •White noise:w(xi, yj), 〈w〉 = 0, 〈w2〉 = Φ0
Transform in Fourier space:
wij = w(kx, ky) =
N−1∑
i=0
N−1∑
j=0
w(xi, yi) exp[−i(kxxi + kyyj)], 〈ww∗〉 = 〈w2〉 = Φ0
Create a filter
fij = f (kx, ky) =
[
K2x+K2
y
k2l
]−(1+H)/2
, for 1 ≤√
K2x+K2
y
kl≤ ζ
0, elsewhere,,
where Kx = (s+ 1)π/L− skx, Ky = (t+ 1)π/L− tky for s, t ∈ {−1, 1}, ζ = ks/kl
Filter white noise:
zij = z(kx, ky) =ℜ(fij)[
ℜ(wij) + iℑ(wij)]
Back to real space:
z(xi, yj) =N−1∑
l=0
N−1∑
m=0zlm exp[i2π(lxi +myj)/L]
V.A. Yastrebov Lecture 4 69/100
![Page 70: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/70.jpg)
Synthesized rough surfaces: in equations II •
Power spectral density:
Φ(kx, ky) = z(kx, ky)z∗(kx, ky) = f 2(kx, ky)w2(kx, ky)
Averaging over multiple samples:
〈Φ(kx, ky)〉 = 〈w2(kx, ky)〉f 2(kx, ky) =
Φ0
[√
K2x+K2
y
kl
]−2(1+H)
, for 1 ≤√
K2x+K2
y
kl≤ ζ
0, elsewhere,
For isotropic surface:
〈Φ(K)〉 =
Φ0(K/kl)−2(1+H), if 1 ≤ K/kl ≤ ζ
0, otherwise.
V.A. Yastrebov Lecture 4 70/100
![Page 71: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/71.jpg)
Effect of parameters: illustration •
• Effect of the high frequency cutoff ks
Fig. Power spectral density (Fourier space)and corresponding rough surface (real space) for
kl = 4, ks = 8
V.A. Yastrebov Lecture 4 71/100
![Page 72: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/72.jpg)
Effect of parameters: illustration •
• Effect of the high frequency cutoff ks
Fig. Power spectral density (Fourier space)and corresponding rough surface (real space) for
kl = 4, ks = 16
V.A. Yastrebov Lecture 4 72/100
![Page 73: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/73.jpg)
Effect of parameters: illustration •
• Effect of the high frequency cutoff ks
Fig. Power spectral density (Fourier space)and corresponding rough surface (real space) for
kl = 4, ks = 32
V.A. Yastrebov Lecture 4 73/100
![Page 74: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/74.jpg)
Effect of parameters: illustration •
• Effect of the high frequency cutoff ks
Fig. Power spectral density (Fourier space)and corresponding rough surface (real space) for
kl = 4, ks = 64
V.A. Yastrebov Lecture 4 74/100
![Page 75: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/75.jpg)
Effect of parameters: illustration •
• Effect of the high frequency cutoff ks
Fig. Power spectral density (Fourier space)and corresponding rough surface (real space) for
kl = 4, ks = 128
V.A. Yastrebov Lecture 4 75/100
![Page 76: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/76.jpg)
Effect of parameters: illustration •
• Effect of the lower frequency cutoff kl for ks/kl = const
Fig. Power spectral density (Fourier space)and corresponding rough surface (real space) for
kl = 1, ks = 43
V.A. Yastrebov Lecture 4 76/100
![Page 77: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/77.jpg)
Effect of parameters: illustration •
• Effect of the lower frequency cutoff kl for ks/kl = const
Fig. Power spectral density (Fourier space)and corresponding rough surface (real space) for
kl = 4, ks = 171
V.A. Yastrebov Lecture 4 77/100
![Page 78: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/78.jpg)
Effect of parameters: illustration •
• Effect of the lower frequency cutoff kl for ks/kl = const
Fig. Power spectral density (Fourier space)and corresponding rough surface (real space) for
kl = 12, ks = 512
V.A. Yastrebov Lecture 4 78/100
![Page 79: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/79.jpg)
Effect of parameters: illustration •
• Effect of the ratio of the higher cutoff to the discretization ks/N
Fig. Power spectral densities (Fourier space)and corresponding rough surfaces (real space) for
kl = 12, ks/N = 1/8 VS kl = 12, ks/N = 1/2
V.A. Yastrebov Lecture 4 79/100
![Page 80: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/80.jpg)
Effect of parameters: illustration •
• Effect of the ratio of the higher cutoff to the discretization ks/N
Fig. Power spectral densities (Fourier space)and corresponding rough surfaces (real space) for
kl = 12, ks/N = 1/8 (fine) VS kl = 12, ks/N = 1/2 (too coarse)for mechanical simulations
V.A. Yastrebov Lecture 4 80/100
![Page 81: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/81.jpg)
Effect of parameters: illustration •
• Effect of the discretisation (single asperity)
Fig. Effect of the mesh on mechanical response
V.A. Yastrebov Lecture 4 81/100
![Page 82: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/82.jpg)
Effect of parameters: illustration •
• Data interpolation (Shanon, bi-cubic Bézier surfaces)
Fig. Bi-cubic Bézier interpolation of an experimental rough surface
[1]Hyun, Robbins, Tribol. Int. (2007)[2] Yastrebov, Durand, Proudhon, Cailletaud, C.R. Mécan. (2011)
V.A. Yastrebov Lecture 4 82/100
![Page 83: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/83.jpg)
Effect of parameters •Effect of parameters:
kl low frequency cutoff- representativity/normality[1,2,3]
ks high frequency cutoff- smoothness and densityof asperities
ζ = ks/kl ratio[3]
- breadth of the spectrum
α ∼ ζ2H
Nayak’s parameter α is the centralcharacteristic of roughness in asperitybased mechanical models.
[1] Vallet, Lasseux, Sainsot, Zahouani, Tribol. Int. (2009)
[2] Yastrebov, Durand, Proudhon, Cailletaud, C.R. Mécan. (2011)
[3] Yastrebov, Anciaux, Molinari, Phys. Rev. E (2012)
[4] Yastrebov, Anciaux, Molinari, Int. J. Solids Struct. (2015)
Fig. 3D and radial power spectral densities
V.A. Yastrebov Lecture 4 83/100
![Page 84: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/84.jpg)
Effect of parameters •Effect of parameters:
kl low frequency cutoff- representativity/normality[1,2,3]
ks high frequency cutoff- smoothness and densityof asperities
ζ = ks/kl ratio[3]
- breadth of the spectrum
α ∼ ζ2H
Nayak’s parameter α is the centralcharacteristic of roughness in asperitybased mechanical models.
[1] Vallet, Lasseux, Sainsot, Zahouani, Tribol. Int. (2009)
[2] Yastrebov, Durand, Proudhon, Cailletaud, C.R. Mécan. (2011)
[3] Yastrebov, Anciaux, Molinari, Phys. Rev. E (2012)
[4] Yastrebov, Anciaux, Molinari, Int. J. Solids Struct. (2015)
V.A. Yastrebov Lecture 4 84/100
![Page 85: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/85.jpg)
Effect of parameters •Effect of parameters:
kl low frequency cutoff- representativity/normality[1,2,3]
ks high frequency cutoff- smoothness and densityof asperities
ζ = ks/kl ratio[3]
- breadth of the spectrum
α ∼ ζ2H
Nayak’s parameter α is the centralcharacteristic of roughness in asperitybased mechanical models.
[1] Vallet, Lasseux, Sainsot, Zahouani, Tribol. Int. (2009)
[2] Yastrebov, Durand, Proudhon, Cailletaud, C.R. Mécan. (2011)
[3] Yastrebov, Anciaux, Molinari, Phys. Rev. E (2012)
[4] Yastrebov, Anciaux, Molinari, Int. J. Solids Struct. (2015)
V.A. Yastrebov Lecture 4 85/100
![Page 86: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/86.jpg)
Interconnection of parameters •
Spectral moment and kl, ks,H:
m0p ≈ mp0 ≈ Φ0
ks∫
kl
2π∫
0
[
k cos(ϕ)]p
(k/kl)−2(1+H) kdkdϕ = Φ0k
p+2
l
ζp−2H − 1
p − 2HT(p)
with T(p) =
2π∫
0
cosp(ϕ)dϕ =
2π, if p = 0;
π, if p = 2;
3π/4, if p = 4.
Nayak’s parameter
α(H, ζ) =3
2
(1 −H)2
H(H − 2)
(ζ−2H − 1)(ζ4−2H − 1)
(ζ2−2H − 1)2
Asperity density
D =
√3
18π
m4
m2=
√3
24π
1 −H
2 −H
ζ4−2H − 1
ζ2−2H − 1k2
l
V.A. Yastrebov Lecture 4 86/100
![Page 87: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/87.jpg)
Interconnection of parameters •
Numerical verification on 100 000 generated rough surfaces with 2048×2048 points
V.A. Yastrebov Lecture 4 87/100
![Page 88: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/88.jpg)
Interconnection of parameters •
12
48
1616
32
64
128
256
5121
510
50100
5001000
5000
Nayak's parameter
k l
ks
H=0.8
H=0.3
32
Nayak’s parameter test
V.A. Yastrebov Lecture 4 88/100
![Page 89: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/89.jpg)
Interconnection of parameters •
12
48
1616
32
64
128
256
5121
510
50100
5001000
5000
32
Asperity density
k l
ks
H=0.8
H=0.3
Asperity density test
V.A. Yastrebov Lecture 4 89/100
![Page 90: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/90.jpg)
Asperity analysis •
Detect summits (zij higher than neighbouring points) and evaluatesecond derivatives
∂2z
∂x2=
zi+1j + zi−1j − 2zij
2∆x2;∂2z
∂y2=
zi+1j + zi−1j − 2zij
2∆x2
∂2z
∂x∂y=
zi+1j+1 + zi+1j+1 − zi+1j−1 − zi−1j+1
4∆x2
Principal curvatures κ1,2:
κ1,2 =12
(
∂2z∂x2 +
∂2z∂y2
)
±√
(
∂2z∂x∂y
)2+ 1
4
(
∂2z∂x2 − ∂2z
∂y2
)2
Saddle point κ1κ2 < 0, extrema κ1κ2 > 0
Mean curvature which can be safely used in Hertz theory: κ =√κ1κ2
V.A. Yastrebov Lecture 4 90/100
![Page 91: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/91.jpg)
Asperity analysis •
Rough surfaces and associated asperities
V.A. Yastrebov Lecture 4 91/100
![Page 92: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/92.jpg)
Asperity analysis •
Rough surfaces and associated asperities
V.A. Yastrebov Lecture 4 92/100
![Page 93: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/93.jpg)
Asperity analysis •
Rough surfaces and associated asperities
V.A. Yastrebov Lecture 4 93/100
![Page 94: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/94.jpg)
Examples •
[1] Yastrebov et al, Three-level multi-scale modeling of electrical contacts sensitivity study and experimentalvalidation, Proceedings of Holm Conference, 2015.
V.A. Yastrebov Lecture 4 94/100
![Page 95: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/95.jpg)
Examples •
Asperity curvatures
V.A. Yastrebov Lecture 4 95/100
![Page 96: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/96.jpg)
Examples •
[1] Yastrebov et al, Three-level multi-scale modeling of electrical contacts sensitivity study and experimentalvalidation, Proceedings of Holm Conference, 2015.
V.A. Yastrebov Lecture 4 96/100
![Page 97: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/97.jpg)
Examples •
Fatigue & creep fracture surfaces (Ti-alloy)in collaboration with A. Marchenko
V.A. Yastrebov Lecture 4 97/100
![Page 98: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/98.jpg)
Examples •
Fatigue fracture surface (Co-alloy), particles WCCourtesy of V. Esin
V.A. Yastrebov Lecture 4 98/100
![Page 99: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/99.jpg)
Examples •
Crater topography and PSDIn collaboration with D. Tkalich (NTNU, Sintef)
V.A. Yastrebov Lecture 4 99/100
![Page 100: Contact Mechanics and Elements of Tribology · Contact Mechanics and Elements of Tribology Lecture 4. Surface roughness Vladislav A. Yastrebov MINES ParisTech, PSL Research University,](https://reader033.vdocument.in/reader033/viewer/2022050606/5fada8da950fba58d4195c04/html5/thumbnails/100.jpg)
, Thank you for your attention!