Minimal polynomial (linear algebra): Difference between revisions

From formulasearchengine
Jump to navigation Jump to search
m remove unreferenced tag
plain html to {{math| ...}} and consistency of notation ''F'' and '''F''' also there was an error with the syntax of one of the formulas
Line 1: Line 1:
{{Technical|date=March 2011}}
Andrew Berryhill is what his spouse loves to call him and he completely digs that title. Her family members lives in Alaska but her husband wants them to move. She is really fond of caving but she doesn't have the time lately. Office supervising is my profession.<br><br>My homepage :: telephone psychic ([http://fashionlinked.com/index.php?do=/profile-13453/info/ fashionlinked.com])
 
{{continuum mechanics|cTopic=[[Solid mechanics]]}}
[[File:Kontakt Spannungsoptik.JPG|290px|thumb|Stresses in a contact area loaded simultaneously with a normal and a tangential force. Stresses were made visible using [[photoelasticity]].]]
 
'''Contact mechanics''' is the study of the [[Deformation (mechanics)|deformation]] of [[solids]] that touch each other at one or more points.<ref name=Johnson>Johnson, K. L, 1985, '''Contact mechanics''', Cambridge University Press.</ref><ref name=Popov>Popov, Valentin L., 2010,  ''Contact Mechanics and Friction. Physical Principles and Applications'', Springer-Verlag, 362 p., ISBN 978-3-642-10802-0.</ref> The physical and mathematical formulation of the subject is built upon the [[mechanics of materials]] and [[continuum mechanics]] and focuses on computations involving [[Elasticity (physics)|elastic]], [[viscoelasticity|viscoelastic]], and [[Plastic Deformation|plastic]] bodies in [[Statics|static]] or [[Dynamics (physics)|dynamic]] contact. Central aspects in contact mechanics are the [[pressure]]s and [[adhesion]] acting perpendicular to the contacting bodies' surfaces, the [[normal vector|normal direction]], and the [[friction]]al [[Stress (mechanics)|stresses]] acting [[Tangential and normal components|tangentially]] between the surfaces. This page focuses mainly on the normal direction, i.e. on frictionless contact mechanics. [[Frictional contact mechanics]] is discussed separately.
 
Contact mechanics is foundational to the field of mechanical [[engineering]]; it provides necessary information for the safe and energy efficient design of technical systems and for the study of [[tribology]] and [[indentation hardness]]. Principles of contacts mechanics can be applied in areas such as locomotive wheel-rail contact, [[coupling]] devices, [[Brake|braking]] systems, [[tire]]s, [[bearing (mechanical)|bearings]], [[Internal combustion engine|combustion engines]], mechanical [[Linkage (mechanical)|linkage]]s, [[gasket]] seals, [[metalworking]], metal forming, [[ultrasonic welding]], [[Electrical connector|electrical contacts]], and many others. Current challenges faced in the field may include [[stress analysis]] of contact and coupling members and the influence of [[lubrication]] and material [[design]] on [[friction]] and [[wear]]. Applications of contact mechanics further extend into the [[Microtechnology|micro]]- and [[Nanotechnology|nanotechnological]] realm.
 
The original work in contact mechanics dates back to 1882 with the publication of the paper "On the contact of elastic solids"<ref>H. Hertz, Über die berührung fester elastischer Körper (On the contact of rigid elastic solids). In: ''[http://www.archive.org/details/cu31924012500306 Miscellaneous Papers]''. Jones and Schott, Editors, ''J. reine und angewandte Mathematik'' 92, Macmillan, London (1896), p. 156 English translation: Hertz, H.</ref> ([http://gdz.sub.uni-goettingen.de/no_cache/dms/load/img/?IDDOC=251917 "Ueber die Berührung fester elastischer Körper"]) by [[Heinrich Hertz]]. Hertz was attempting to understand how the optical properties of multiple, stacked [[lens (optics)|lenses]] might change with the [[force]] holding them together. Hertzian contact stress refers to the localized stresses that develop as two curved surfaces come in contact and deform slightly under the imposed loads. This amount of deformation is dependent on the [[modulus of elasticity]] of the material in contact. It gives the contact stress as a function of the normal contact force, the radii of curvature of both bodies and the modulus of elasticity of both bodies. Hertzian contact stress forms the foundation for the equations for load bearing capabilities and [[fatigue (material)|fatigue]] life in bearings, gears, and any other bodies where two surfaces are in contact.
 
== History ==
[[File:Hertz contact animated.gif|290px|thumb|widthpx|When a sphere is pressed against an elastic material, the contact area increases.]]
Classical contact mechanics is most notably associated with Heinrich Hertz.<ref name=Hertz>Hertz, H. R., 1882, ''Ueber die Beruehrung elastischer Koerper (On Contact Between Elastic Bodies),'' in '''Gesammelte Werke''' (Collected Works), Vol. 1, Leipzig, Germany, 1895.</ref> In 1882, Hertz solved the contact problem of two elastic bodies with curved surfaces. This still-relevant classical solution provides a foundation for modern problems in contact mechanics. For example, in [[mechanical engineering]] and [[tribology]], ''Hertzian contact stress'' is a description of the stress within mating parts. The Hertzian contact stress usually refers to the stress close to the area of contact between two spheres of different radii.
 
It was not until nearly one hundred years later that [[Kenneth L. Johnson|Johnson]], Kendall, and Roberts found a similar solution for the case of [[Adhesion|adhesive]] contact.<ref name=jkr>K. L. Johnson and K. Kendall and A. D. Roberts, Surface energy and the contact of elastic solids, Proc. R. Soc. London A 324 (1971) 301-313</ref>  This theory was rejected by [[Boris Derjaguin]] and co-workers<ref name=Maugis>D. Maugis, Contact, Adhesion and Rupture of Elastic Solids, Springer-Verlag, Solid-State Sciences, Berlin 2000, ISBN 3-540-66113-1</ref> who proposed a different theory of adhesion<ref name=dmt>B. V. Derjaguin and V. M. Muller and Y. P. Toporov, Effect of contact deformations on the adhesion of particles, J. Colloid Interface Sci. 53 (1975) 314--325</ref> in the 1970s. The Derjaguin model came to be known as the DMT (after Derjaguin, Muller and Toporov) model,<ref name=dmt/> and the Johnson et al. model came to be known as the JKR (after Johnson, Kendall and Roberts) model for adhesive elastic contact. This rejection proved to be instrumental in the development of the Tabor<ref name=Tabor>D. Tabor, The hardness of solids, J. Colloid Interface Sci. 58 (1977) 145-179</ref> and later Maugis<ref name=Maugis/><ref name=Maugis1>D. Maugis, Adhesion of spheres: The JKR-DMT transition using a Dugdale model, J. Colloid Interface Sci. 150 (1992) 243--269</ref> parameters that quantify which contact model (of the JKR and DMT models) represent adhesive contact better for specific materials.
 
Further advancement in the field of contact mechanics in the mid-twentieth century may be attributed to names such as Bowden and [[David Tabor|Tabor]]. Bowden and Tabor were the first to emphasize the importance of surface roughness for bodies in contact.<ref name=bowden1939>, Bowden, FP and Tabor, D., 1939, ''The area of contact between stationary and between moving surfaces'', Proceedings of the Royal Society of London. Series A, Mathematical and Physical Sciences, 169(938), pp. 391--413.</ref><ref name=bowden2001>Bowden, F.P. and Tabor, D., 2001, '''The friction and lubrication of solids''', Oxford University Press.</ref> Through investigation of the surface roughness, the true contact area between friction partners is found to be less than the apparent contact area. Such understanding also drastically changed the direction of undertakings in tribology. The works of Bowden and Tabor yielded several theories in contact mechanics of rough surfaces.
 
The contributions of Archard (1957)<ref name=archard1957>Archard, JF, 1957, ''Elastic deformation and the laws of friction'',  Proceedings of the Royal Society of London. Series A, Mathematical and Physical Sciences, 243(1233), pp.190--205.</ref> must also be mentioned in discussion of pioneering works in this field. Archard concluded that, even for rough elastic surfaces, the contact area is approximately proportional to the [[normal force]]. Further important insights along these lines were provided by Greenwood and Williamson (1966),<ref>Greenwood, JA and Williamson, JBP., 1966, ''Contact of nominally flat surfaces'', Proceedings of the Royal Society of London. Series A, Mathematical and Physical Sciences, pp. 300-319.</ref> Bush (1975),<ref>Bush, AW and Gibson, RD and Thomas, TR., 1975, ''The elastic contact of a rough surface'', Wear, 35(1), pp. 87-111.</ref> and Persson (2002).<ref>Persson, BNJ and Bucher, F. and Chiaia, B., 2002, ''Elastic contact between randomly rough surfaces: Comparison of theory with numerical results'', Physical Review B, 65(18), p. 184106.</ref> The main findings of these works were that the true contact surface in rough materials is generally proportional to the normal force, while the parameters of individual micro-contacts (i.e., pressure, size of the micro-contact) are only weakly dependent upon the load.
 
== Classical solutions for non-adhesive elastic contact ==
The theory of contact between elastic bodies can be used to find contact areas and indentation depths for simple geometries.  Some commonly used solutions are listed below.  The theory used to compute these solutions is discussed later in the article.
 
=== Contact between a sphere and an elastic half-space ===
[[File:Kontakt Kugel Ebene.jpg|thumb|Contact between a sphere and an elastic half-space]]
An elastic [[sphere]] of [[radius]] <math>R</math> indents an elastic [[Half-space (geometry)|half-space]] to depth <math>d</math>, and thus creates a contact area of radius
:<math>a=\sqrt{2Rd-d^2}</math>
:<math>a=\sqrt{2Rd}</math> , for d<<R
The applied force <math>F</math> is related to the displacement <math>d</math> by
:<math>
  F=\tfrac{4}{3} E^*R^{1/2}d^{3/2}
</math>
where
:<math>
  \frac{1}{E^*}=\frac{1-\nu^2_1}{E_1}+\frac{1-\nu^2_2}{E_2}
</math>
and  <math>E_1</math>,<math>E_2</math> are the [[Elastic modulus|elastic moduli]] and <math>\nu_1</math>,<math>\nu_2</math> the [[Poisson's ratio]]s associated with each body.
 
=== Contact between two spheres ===
[[File:Kontakt Kugel Kugel.jpg|120px|thumb|Contact between two spheres]]
[[File:Kontakt gekreuzter Zylinder.jpg|thumb|Contact between two crossed cylinders of equal radius]]
For contact between two spheres of radii <math>R_1</math> and <math>R_2</math>, the area of contact is a circle of radius <math>a</math>.  The distribution of normal traction in the contact area as a function of distance from the center of the circle is<ref name="Johnson"/>
:<math>
  p(r) = p_0\left(1-\frac{r^2}{a^2}\right)^{1/2}
</math>
where <math>p_0</math> is the maximum contact pressure given by
:<math>
  p_0 = \cfrac{3F}{2\pi a^2} = \cfrac{1}{\pi}\left(\cfrac{6F{E^*}^2}{R^2}\right)^{1/3}
</math>
where the effective radius <math>R</math> is defined as
:<math>
  \frac{1}{R} = \frac{1}{R_1}+\frac{1}{R_2}
</math>
The radius of the circle is related to the applied load <math>F</math> by the equation
:<math>
  a^3 = \cfrac{3 F R}{4 E^*}
</math>
The depth of indentation <math>d</math> is related to the maximum contact pressure by
:<math>
  d = \cfrac{a^2}{R} = \left(\cfrac{9F^2}{16R{E^*}^2}\right)^{1/3}
</math>
The maximum shear stress occurs in the interior at <math>z\approx 0.49a</math> for <math>\nu = 0.33</math>.
 
===Contact between two crossed cylinders of equal radius <math>R</math>===
This is equivalent to contact between a sphere of radius <math>R</math> and a [[Plane (geometry)|plane]] (see above).
 
=== Contact between a rigid cylinder and an elastic half-space ===
[[File:Kontakt Zylindrischer Indenter Ebene.jpg|thumb|Contact between a rigid cylindrical indenter and an elastic half-space]]
 
If a rigid [[Cylinder (geometry)|cylinder]] is pressed into an elastic half-space, it creates a pressure distribution described by<ref name=Sneddon>Sneddon, I. N., 1965, ''The Relation between Load and Penetration in the Axisymmetric Boussinesq Problem for a Punch of Arbitrary Profile.'' Int. J. Eng. Sci. v. 3, pp. 47–57.</ref>
:<math>
  p(r)=p_0\left(1-\frac{r^2}{a^2}\right)^{-1/2}
  </math>
where <math>a</math> is the radius of the cylinder and
:<math>
  p_0=\frac{1}{\pi}E^*\frac{d}{a}
</math>
The relationship between the indentation depth and the normal force is given by
:<math>
  F=2aE^*d\,
</math>
 
=== Contact between a rigid conical indenter and an elastic half-space===
[[File:Kontakt Kegel Ebene.jpg|thumb|Contact between a rigid conical indenter and an elastic half-space]]
In the case of [[indentation]] of an elastic half-space of Young's modulus <math>E</math> using a rigid [[Cone (geometry)|conical]] indenter, the depth of the contact region <math>\epsilon</math> and contact radius <math>a</math> are related by<ref name=Sneddon/>
:<math>
  \epsilon=a\tan\theta
    </math>
with <math>\theta</math> defined as the angle between the plane and the side surface of the cone.  The total indentation depth <math>d</math> is given by:<math>
d=\frac{\pi}{2}\epsilon</math>
 
The total force is
:<math>
  F =\frac{\pi E}{2 \left(1-\nu^2\right)} a^2 \tan \theta=\frac{2E}{\pi\left(1-\nu^2\right)}\frac{d^2}{\tan \theta}
</math>
 
 
The pressure distribution is given by
:<math>
  p {\left( r \right)}=\frac{Ed}{\pi a\left(1-\nu^2\right)}\ln\left(\frac{a}{r}+\sqrt{\left(\frac{a}{r}\right)^2-1}\right)
= \frac{Ed}{\pi a\left(1-\nu^2\right)} \cosh^{-1}\left(\frac{a}{r}\right)
</math>
The stress has a [[logarithm]]ic [[Mathematical singularity|singularity]] at the tip of the cone.
 
===Contact between two cylinders with parallel axes===
[[File:Kontakt paralleler Zylinder.jpg|thumb|Contact between two cylinders with parallel axes]]
In contact between two cylinders with parallel axes, the force is linearly proportional to the indentation depth:
:<math>
  F=\frac{\pi}{4}E^*Ld
</math>
The radii of curvature are entirely absent from this relationship. The contact radius is described through the usual relationship
:<math>a=\sqrt{Rd}
</math>
with
:<math>\frac{1}{R}=\frac{1}{R_1}+\frac{1}{R_2}</math>
as in contact between two spheres. The maximum pressure is equal to
:<math>
  p_0=\left(\frac{E^*F}{\pi LR}\right)^{1/2}
</math>
 
=== The Method of Dimensionality Reduction ===
[[File:WikipediabilderKap 4.jpg|thumb|Contact between a sphere and an elastic half-space and one-dimensional replaced model]]
Many contact problems can be solved easily with the Method of Dimensionality Reduction.
In this method, the initial three-dimensional system is replaced with a contact of a body with a linear elastic or viscoelastic foundation (see Fig). The properties of one-dimensional systems coincide hereby exactly with those of the original three-dimensional system, if the form of the bodies is modified and the elements of the foundation are defined according to the rules of the RMD.
<ref name=Popov2>Popov, V.L., ''Method of reduction of dimensionality in contact and friction mechanics: A linkage between micro and macro scales,'' Friction, 2013, v.1, N. 1, pp.41–62. </ref>
<ref name=Popov3>Popov, V.L. and Heß, M., '''Methode der Dimensionsreduktion in Kontaktmechanik und Reibung''', Springer, 2013. </ref>
 
== Hertzian theory of non-adhesive elastic contact ==
The classical theory of contact focused primarily on non-adhesive contact where no tension force is allowed to occur within the contact area, i.e., contacting bodies can be separated without adhesion forces.  Several analytical and numerical approaches have been used to solve contact problems that satisfy the no-adhesion condition.  Complex forces and [[moment (physics)|moments]] are transmitted between the bodies where they touch, so problems in contact mechanics can become quite sophisticated.  In addition, the contact stresses are usually a nonlinear function of the deformation.  To simplify the solution procedure, a [[frame of reference]] is usually defined in which the objects (possibly in motion relative to one another) are static. They interact through surface tractions (or pressures/stresses) at their interface.
 
As an example, consider two objects which meet at some surface <math>S</math> in the (<math>x</math>,<math>y</math>)-plane with the <math>z</math>-axis assumed normal to the surface. One of the bodies will experience a normally-directed [[pressure]] distribution <math>p_z=p(x,y)=q_z(x,y)</math> and in-plane [[surface traction]] distributions <math>q_x=q_x(x,y)</math> and <math>q_y=q_y(x,y)</math> over the region <math>S</math>. In terms of a [[Newtonian mechanics|Newtonian]] force balance, the forces:
:<math>
  P_z = \int_S p(x,y)~ \mathrm{d}A ~;~~ Q_x = \int_S q_x(x,y)~ \mathrm{d}A ~;~~ Q_y = \int_S q_y(x,y)~ \mathrm{d}A
</math>
must be equal and opposite to the forces established in the other body. The moments corresponding to these forces:
:<math>
  M_x = \int_S y~p(x,y)~ \mathrm{d}A ~;~~ M_y = \int_S x~p(x,y)~ \mathrm{d}A ~;~~ M_z = \int_S [x~q_y(x,y) - y~q_x(x,y)]~ \mathrm{d}A
  </math>
are also required to cancel between bodies so that they are [[kinematics|kinematically]] immobile.
 
=== Assumptions in Hertzian theory===
The following assumptions are made in determining the solutions of '''Hertzian''' contact problems:
* the strains are small and within the elastic limit,
* each body can be considered an elastic half-space, i.e., the area of contact is much smaller than the characteristic radius of the body,
* the surfaces are continuous and non-conforming, and
* the surfaces are frictionless.
Additional complications arise when some or all these assumptions are violated and such contact problems are usually called '''non-Hertzian'''.
 
=== Analytical solution techniques ===
[[File:Hertz.svg|right|thumb|290px|Contact between two spheres.]]
Analytical solution methods for non-adhesive contact problem can be classified into two types based on the geometry of the area of contact.<ref name=Shigley>Shigley, J.E., Mischke, C.R., 1989, '''Mechanical Engineering Design''', Fifth Edition, Chapter 2, McGraw-Hill, Inc, 1989, ISBN 0-07-056899-5.</ref>  A '''conforming contact''' is one in which the two bodies touch at multiple points before any deformation takes place (i.e., they just "fit together").  A '''non-conforming contact''' is one in which the shapes of the bodies are dissimilar enough that, under zero load, they only touch at a point (or possibly along a line). In the non-conforming case, the contact area is small compared to the sizes of the objects and the [[stress (physics)|stresses]] are highly concentrated in this area. Such a contact is called ''concentrated'', otherwise it is called ''diversified''.
 
A common approach in [[linear elasticity]] is to [[Superposition principle|superpose]] a number of solutions each of which corresponds to a point load acting over the area of contact.  For example, in the case of loading of a [[half-plane]], the [[Flamant solution]] is often used as a starting point and then generalized to various shapes of the area of contact.  The force and moment balances between the two bodies in contact act as additional constraints to the solution.
 
==== Point contact on a (2D) half-plane ====
{{Main|Flamant solution}}
[[File:Point_contact_fig.svg|thumb|right|290px|Schematic of the loading on a plane by force P at a point (0,0).]]
A starting point for solving contact problems is to understand the effect of a "point-load" applied to an isotropic, homogeneous, and linear elastic half-plane, shown in the figure to the right.  The problem may be either be [[plane stress]] or [[plane strain]].  This is a [[boundary value problem]] of linear elasticity subject to the traction [[boundary condition]]s:
:<math>\sigma_{xz} (x,0)=0 ~;~~ \sigma_z (x,z)=-P\delta (x,z)</math>
where <math>\delta(x,z)</math> is the [[Dirac delta function]].  The boundary conditions state that there are no shear stresses on the surface and a singular normal force P is applied at (0,0). Applying these conditions to the governing equations of elasticity produces the result
:<math>
  \begin{align}
  \sigma_{xx} & =-\frac{2P}{\pi}\frac{x^2z}{(x^2+z^2)^2} ~;~~
  \sigma_{zz} =-\frac{2P}{\pi}\frac{z^3}{(x^2+z^2)^2} \\
  \sigma_{xz} & =-\frac{2P}{\pi}\frac{xz^2}{(x^2+z^2)^2}
  \end{align}
</math>
for some point, <math>(x,y)</math>, in the half-plane. The circle shown in the figure indicates a surface on which the maximum shear stress is constant. From this stress field, the [[strain (materials science)|strain]] components and thus the [[displacement field (mechanics)|displacements]] of all material points may be determined.
 
==== Line contact on a (2D) half-plane ====
 
===== Normal loading over a region <math>(a,b)</math>=====
Suppose, rather than a point load <math>P</math>, a distributed load <math>p(x)</math> is applied to the surface instead, over the range <math>a<x<b</math>. The principle of linear superposition can be applied to determine the resulting stress field as the solution to the [[integral]] equations:
:<math>
  \begin{align}
    \sigma_{xx} & =-\frac{2z}{\pi}\int_a^b\frac{p(x')(x-x')^2\, dx'}{[(x-x')^2+z^2]^2} ~;~~
    \sigma_{zz} =-\frac{2z^3}{\pi}\int_a^b\frac{p(x')\, dx'}{[(x-x')^2+z^2]^2} \\
    \sigma_{xz} & =-\frac{2z^2}{\pi}\int_a^b\frac{p(x')(x-x')\, dx'}{[(x-x')^2+z^2]^2}
  \end{align}
  </math>
 
===== Shear loading over a region <math>(a,b)</math>=====
The same principle applies for loading on the surface in the plane of the surface. These kinds of tractions would tend to arise as a result of friction. The solution is similar the above (for both singular loads <math>Q</math> and distributed loads <math>q(x)</math>) but altered slightly:
:<math>
  \begin{align}
    \sigma_{xx} & =-\frac{2}{\pi}\int_a^b\frac{q(x')(x-x')^3\, dx'}{[(x-x')^2+z^2]^2} ~;~~
    \sigma_{zz} =-\frac{2z^2}{\pi}\int_a^b\frac{q(x')(x-x')\, dx'}{[(x-x')^2+z^2]^2} \\
    \sigma_{xz} & =-\frac{2z}{\pi}\int_a^b\frac{q(x')(x-x')^2\, dx'}{[(x-x')^2+z^2]^2}
  \end{align}
  </math>
These results may themselves be superposed onto those given above for normal loading to deal with more complex loads.
 
====  Point contact on a (3D) half-space ====
Analogously to the Flamant solution for the 2D half-plane, fundamental solutions are known for the linearly elastic 3D half-space as well. These were found by [[Joseph Valentin Boussinesq|Boussinesq]] for a concentrated normal load and by Cerutti for a tangential load. See the section on this in [[Linear elasticity#Solutions for elastostatic cases|Linear elasticity]].
 
=== Numerical solution techniques ===
Distinctions between conforming and non-conforming contact do not have to be made when numerical solution schemes are employed to solve contact problems. These methods do not rely on further assumptions within the solution process since they base solely on the general formulation of the underlying equations
<ref name=Kalker>Kalker, J.J. 1990, ''Three-Dimensional Elastic Bodies in Rolling Contact''. (Kluwer Academic Publishers: Dordrecht).</ref>
<ref name=Wriggers>Wriggers, P. 2006, ''Computational Contact Mechanics''. 2nd ed. (Springer Verlag: Heidelberg).</ref>
<ref name=Laursen>Laursen, T. A., 2002, ''Computational Contact and Impact Mechanics: Fundamentals of Modeling Interfacial Phenomena in Nonlinear Finite Element Analysis'', (Springer Verlag: New York).</ref>
<ref name=Acary>Acary V. and Brogliato B., 2008,''Numerical Methods for Nonsmooth Dynamical Systems. Applications in Mechanics and Electronics''. Springer Verlag, LNACM 35, Heidelberg.</ref>
.<ref name=Popov1>Popov, Valentin L., 2009, ''Kontaktmechanik und Reibung. Ein Lehr- und Anwendungsbuch von der Nanotribologie bis zur numerischen Simulation'', Springer-Verlag, 328 S., ISBN 978-3-540-88836-9.</ref> Besides the standard equations describing the deformation and motion of bodies two additional inequalities can be formulated. The first simply restricts the motion and deformation of the bodies by the assumption that no penetration can occur. Hence the gap <math>g_N</math> between two bodies can only be positive or zero
:<math>g_N \ge 0</math>
where <math>g_N = 0</math> denotes contact. The second assumption in contact mechanics is related to the fact, that no tension force is allowed to occur within the contact area (contacting bodies can be lifted up without adhesion forces). This leads to an inequality which the stresses  have to obey at the contact interface. It is formulated for the contact pressure <math>p_N = \mathbf{t} \cdot \mathbf{n}</math>
:<math>p_N \le 0 \,.</math>
Since for contact, <math>g_N = 0 </math>, the contact pressure is always negative, <math>p_N < 0 </math>, and further for non contact the gap is open, <math>g_N > 0 </math>,  and the contact pressure is zero, <math>p_N = 0 </math>, the so-called [[Karush–Kuhn–Tucker conditions|Kuhn–Tucker]] form of the contact constraints can be written as
:<math>g_N \ge 0\,, \quad p_N \le 0\,, \quad p_N\,g_N = 0\,.</math>
These conditions are valid in a general way. The mathematical formulation of the gap depends upon the kinematics of the underlying theory of the solid (e.g., linear or nonlinear solid in two- or three dimensions, [[beam theory|beam]] or [[plate theory|shell]] model).
 
== Non-adhesive contact between rough surfaces ==
When two bodies with rough surfaces are pressed into each other, the true contact area <math>A</math> is much smaller than the apparent contact area <math>A_0</math>. In contact between a "random rough" surface and an elastic half-space, the true contact area is related to the normal force <math>F</math> by<ref name=Johnson/><ref name=Greenwood/><ref name=Mikic/><ref name=Hyun>Hyun, S., and M.O. Robbins, 2007, ''Elastic contact between rough surfaces: Effect of roughness at large and small wavelengths.'' Tribology International, v.40, pp. 1413-1422.</ref>
:<math>
  A=\frac{\kappa}{E^*h'}F
</math>
with  <math>h'</math>  equal to the root mean square (also known as the quadratic mean) of the surface slope and <math>\kappa \approx2</math> . The median pressure in the true contact surface
:<math>
  p_{\mathrm{av}} =\frac{F}{A}\approx\frac{1}{2}E^*h'
</math>
can be reasonably estimated as half of the effective elastic modulus <math>E^*</math> multiplied with the root mean square of the surface slope <math>h'</math> .
 
For the situation where the asperities on the two surfaces have a Gaussian height distribution and the peaks can be assumed to be spherical,<ref name=Greenwood>Greenwood, J. A. and Williamson, J. B. P., (1966), ''Contact of nominally flat surfaces,'' Proceedings of the Royal Society of London. Series A, Mathematical and Physical Sciences, vol. 295, pp. 300--319.</ref> the average contact pressure is sufficient to cause yield when <math> p_{\mathrm{av}} = 1.1\sigma_y \approx 0.39 \sigma_0</math> where <math>\sigma_y</math> is the uniaxial [[yield stress]] and <math>\sigma_0</math> is the indentation hardness.<ref name=Johnson/>  Greenwood and Williamson<ref name=Greenwood/> defined a dimensionless parameter <math>\Psi</math> called the '''plasticity index''' that could be used to determine whether contact would be elastic or plastic.
 
The Greenwood-Williamson model requires knowledge of two statistically dependent quantities; the standard deviation of the surface roughness and the curvature of the asperity peaks.  An alternative definition of the plasticity index has been given by Mikic.<ref name=Mikic>Mikic, B. B., (1974), ''Thermal contact conductance; theoretical considerations'', International Journal of Heat and Mass Transfer, 17(2), pp. 205-214.</ref>  Yield occurs when the pressure is greater than the uniaxial yield stress. Since the yield stress is proportional to the indentation hardness <math>\sigma _0</math>, Micic defined the plasticity index for elastic-plastic contact to be
:<math>\Psi = \frac{E^*h'}{\sigma _0}> \tfrac{2}{3}~.</math>
In this definition <math>\Psi</math> represents the micro-roughness in a state of complete plasticity and only one statistical quantity, the rms slope, is needed which can be calculated from surface measurements. For <math>\Psi <\tfrac{2}{3}</math>, the surface behaves elastically during contact.
 
In both the Greenwood-Williamson and Mikic models the load is assumed to be proportional to the deformed area.  Hence, whether the system behaves plastically or elastically is independent of the applied normal force.<ref name=Johnson/>
 
== Adhesive contact between elastic bodies==
When two solid surfaces are brought into close proximity, they experience attractive [[van der Waals force]]s.  Bradley's van der Waals model<ref name=Bradley>Bradley, RS., 1932, ''The cohesive force between solid surfaces and the surface energy of solids'', Philosophical Magazine Series 7, 13(86), pp. 853--862.</ref> provides a means of calculating the tensile force between two rigid spheres with perfectly smooth surfaces.  The Hertzian model of contact does not consider adhesion possible.  However, in the late 1960s, several contradictions were observed when the Hertz theory was compared with experiments involving contact between rubber and glass spheres.
 
It was observed<ref name=jkr/> that, though Hertz theory applied at large loads, at low loads
* the area of contact was larger than that predicted by Hertz theory,
* the area of contact had a non-zero value even when the load was removed, and
* there was strong adhesion if the contacting surfaces were clean and dry.
This indicated that adhesive forces were at work.  The Johnson-Kendall-Roberts (JKR) model and the Derjaguin-Muller-Toporov (DMT) models were the first to incorporate adhesion into Hertzian contact.
 
=== Bradley model of rigid contact ===
It is commonly assumed that the surface force between two atomic planes at a distance <math>z</math> from each other can be derived from the [[Lennard-Jones potential]]. With this assumption
:<math>
  F(z) = \cfrac{16\gamma}{3 z_0}\left[\left(\cfrac{z}{z_0}\right)^{-9} - \left(\cfrac{z}{z_0}\right)^{-3}\right]
</math>
where <math>F</math> is the force (positive in compression), <math>2\gamma</math> is the total surface energy of ''both'' surfaces per unit area, and <math>z_0</math> is the equilibrium separation of the two atomic planes.
 
The Bradley model applied the Lennard-Jones potential to find the force of adhesion between two rigid spheres. The total force between the spheres is found to be
:<math>
  F_a(z) = \cfrac{16\gamma\pi R}{3}\left[\cfrac{1}{4}\left(\cfrac{z}{z_0}\right)^{-8} - \left(\cfrac{z}{z_0}\right)^{-2}\right] ~;~~  \frac{1}{R} = \frac{1}{R_1} + \frac{1}{R_2}
</math>
where <math>R_1,R_2</math> are the radii of the two spheres.
 
The two spheres separate completely when the ''pull-off force'' is achieved at <math>z = z_0</math> at which point
:<math>
  F_a = F_c = -4\gamma\pi R .
</math>
 
=== Johnson-Kendall-Roberts (JKR) model of elastic contact===
[[File:JKRModel.svg|400px|thumb|Schematic of contact area for the JKR model.]]
[[File:JKR-full-cycle.gif|thumb|right|JKR test with a rigid bead on a deformable planar material: complete cycle]]
To incorporate the effect of adhesion in Hertzian contact, Johnson, Kendall, and Roberts<ref name=jkr/> formulated the JKR theory of adhesive contact using a balance between the stored [[elastic energy]] and the loss in [[surface energy]].  The JKR model considers the effect of contact pressure and adhesion only inside the area of contact. The general solution for the pressure distribution in the contact area in the JKR model is
:<math>
  p(r) = p_0\left(1 - \cfrac{r^2}{a^2}\right)^{1/2} + p_0'\left(1 - \cfrac{r^2}{a^2}\right)^{-1/2}
</math>
Note that in the original Hertz theory, the term containing <math>p_0'</math> was neglected on the ground that tension could not be sustained in the contact zone.  For contact between two spheres
:<math>
  p_0 = \cfrac{2 a E^*}{\pi R} ~;~~
  p_0' = -\left(\cfrac{4\gamma E^*}{\pi a}\right)^{1/2}
</math>
where <math>a\,</math> is the radius of the area of contact, <math>F</math> is the applied force, <math>2\gamma</math> is the total surface energy of '''both''' surfaces per unit contact area,
<math>R_i, E_i, \nu_i,~~i=1,2</math> are the radii, Young's moduli, and Poisson's ratios of the two spheres, and
:<math>
  \frac{1}{R} = \frac{1}{R_1} + \frac{1}{R_2} ~;~~ \frac{1}{E^*} = \frac{1-\nu_1^2}{E_1} + \frac{1-\nu_2^2}{E_2}
</math>
The approach distance between the two spheres is given by
:<math>
  d = \cfrac{\pi a}{2 E^*}(p_0 + 2p_0') = \cfrac{a^2}{R}
</math>
The Hertz equation for the area of contact between two spheres, modified to take into account the surface energy, has the form
:<math>
  a^3 = \cfrac{3R}{4E^*}\left(F + 6\gamma\pi R + \sqrt{12\gamma\pi R F + (6\gamma\pi R)^2}\right)
</math>
When the surface energy is zero, <math>\gamma = 0</math>, the Hertz equation for contact between two spheres is recovered.  When the applied load is zero, the contact radius is
:<math>
  a^3 = \cfrac{9R^2\gamma\pi}{E^*}
</math>
The tensile load at which the spheres are separated, i.e., <math>a = 0</math>, is predicted to be
:<math>
  F_c = -3\gamma\pi R\,
</math>
This force is also called the '''pull-off force'''.
Note that this force is independent of the moduli of the two spheres.  However, there is another possible solution for the value of <math>a</math> at this load.  This is the critical contact area <math>a_c</math>, given by
:<math>
  a_c^3 = \cfrac{9R^2\gamma\pi}{4E^*}
</math>
 
If we define the '''work of adhesion''' as
:<math>
  \Delta\gamma = \gamma_1 + \gamma_2 - \gamma_{12}
</math>
where <math>\gamma_1,\gamma_2</math> are the adhesive energies of the two surfaces and <math>\gamma_{12}</math> is an interaction term, we can write the JKR contact radius as
:<math>
  a^3 = \cfrac{3R}{4E^*}\left(F + 3\Delta\gamma\pi R + \sqrt{6\Delta\gamma\pi R F + (3\Delta\gamma\pi R)^2}\right)
</math>
The tensile load at separation is
:<math>
  F = -\cfrac{3}{2}\Delta\gamma\pi R\,
</math>
and the critical contact radius is given by
:<math>
  a_c^3 = \cfrac{9R^2\Delta\gamma\pi}{4E^*}
</math>
The critical depth of penetration is
:<math>
  d_c = \cfrac{a_c^2}{R} = \left(\cfrac{9}{4}\right)^{\tfrac{2}{3}}(\Delta\gamma)^{\tfrac{2}{3}}\left(\cfrac{\pi^{\tfrac{2}{3}}~R^{\tfrac{1}{3}}}{{E^*}^{\tfrac{2}{3}}}\right)
</math>
 
=== Derjaguin-Muller-Toporov (DMT) model of elastic contact ===
The Derjaguin-Muller-Toporov (DMT) model<ref name=DMT>Derjaguin, BV and Muller, VM and Toporov, Y.P., 1975, ''Effect of contact deformations on the adhesion of particles'', Journal of Colloid and Interface Science, 53(2), pp. 314-326.</ref><ref>Muller, VM and Derjaguin, BV and Toporov, Y.P., 1983, ''On two methods of calculation of the force of sticking of an elastic sphere to a rigid plane'', Colloids and Surfaces, 7(3), pp. 251-259.</ref> is an alternative model for adhesive contact which assumes that the contact profile remains the same as in Hertzian contact but with additional attractive interactions outside the area of contact.
 
The area of contact between two spheres from DMT theory is
:<math>
  a^3 = \cfrac{3R}{4E^*}\left(F + 4\gamma\pi R\right)
</math>
and the pull-off force is
:<math>
  F_c = -4\gamma\pi R\,
</math>
When the pull-off force is achieved the contact area becomes zero and there is no singularity in the contact stresses at the edge of the contact area.
 
In terms of the work of adhesion <math>\Delta\gamma</math>
:<math>
  a^3 = \cfrac{3R}{4E^*}\left(F + 2\Delta\gamma\pi R\right)
</math>
and
:<math>
  F_c = -2\Delta\gamma\pi R\,
</math>
 
=== Tabor coefficient ===
In 1977, Tabor<ref name=Tabor1>Tabor, D., 1977, ''Surface forces and surface interactions'', Journal of Colloid and Interface Science, 58(1), pp. 2-13.</ref> showed that the apparent contradiction between the JKR and DMT theories could be resolved by noting that the two theories were the extreme limits of a single theory parametrized by the '''Tabor coefficient''' (<math>\mu</math>) defined as
:<math>
  \mu := \cfrac{d_c}{z_0} \approx \left[\cfrac{R(\Delta\gamma)^2}{{E^*}^2 z_0^3}\right]^{1/3}
</math>
where <math>z_0</math> is the equilibrium separation between the two surfaces in contact.  The JKR theory applies to large, compliant spheres for which <math>\mu</math> is large.  The DMT theory applies for small, stiff spheres with small values of <math>\mu</math>.
 
=== Maugis-Dugdale model of elastic contact===
[[File:MaugisDugdale.svg|400px|thumb|Schematic of contact area for the Maugis-Dugdale model.]]
Further improvement to the Tabor idea was provided by Maugis<ref name=Maugis1>Maugis, D., 1992, ''Adhesion of spheres: the JKR-DMT transition using a Dugdale model'', Journal of Colloid and Interface Science, 150(1), pp. 243-269.</ref> who represented the surface force in terms of a Dugdale [[cohesive zone]] approximation such that the work of adhesion is given by
:<math>
  \Delta\gamma = \sigma_0~h_0
</math>
where <math>\sigma_0</math> is the maximum force predicted by the Lennard-Jones potential and <math>h_0</math> is the maximum separation obtained by matching the areas under the Dugdale and Lennard-Jones curves (see adjacent figure).  This means that the attractive force is constant for <math>z_0 \le z \le z_0 + h_0</math>.  There is not further penetration in compression.  Perfect contact occurs in an area of radius <math>a</math> and adhesive forces of magnitude <math>\sigma_0</math> extend to an area of radius <math>c > a</math>.  In the region <math>a < r < c</math>, the two surfaces are separated by a distance <math>h(r)</math> with <math>h(a) = 0</math> and <math>h(c) = h_0</math>.  The ratio <math>m</math> is defined as
:<math> m := \cfrac{c}{a}</math>.
 
In the Maugis-Dugdale theory,<ref>Johnson, KL and Greenwood, JA, 1997, ''An adhesion map for the contact of elastic spheres'', Journal of Colloid and Interface Science, 192(2), pp. 326-333.</ref> the surface traction distribution is divided into two parts - one due to the Hertz contact pressure and the other from the Dugdale adhesive stress.  Hertz contact is assumed in the region <math>-a < r < a</math>.  The contribution to the surface traction from the Hertz pressure is given by
:<math>
  p^H(r) = \left(\cfrac{3F^H}{2\pi a^2}\right)\left(1 - \cfrac{r^2}{a^2}\right)^{1/2}
</math>
where the Hertz contact force <math>F^H</math> is given by
:<math>
  F^H = \cfrac{4E^*a^3}{3R}
</math>
The penetration due to elastic compression is
:<math>
  d^H = \cfrac{a^2}{R}
</math>
The vertical displacement at <math>r = c</math> is
:<math>
  u^H(c) = \cfrac{1}{\pi R} \left[a^2(2 - m^2)\sin^{-1}\left(\cfrac{1}{m}\right) + a^2\sqrt{m^2-1}\right]
</math>
and the separation between the two surfaces at <math>r = c</math> is
:<math>
  h^H(c) = \cfrac{c^2}{2R} - d^H + u^H(c)
</math>
The surface traction distribution due to the adhesive Dugdale stress is
:<math>
  p^D(r) = \begin{cases}
              -\cfrac{\sigma_0}{\pi}\cos^{-1}\left[\cfrac{2-m^2-\cfrac{r^2}{a^2}}{m^2\left(1-\cfrac{r^2}{m^2a^2}\right)}\right] & \quad \mathrm{for} \quad r \le a\\
              -\sigma_0 & \quad \mathrm{for} \quad a \le r \le c
            \end{cases}
</math>
The total adhesive force is then given by
:<math>
  F^D = -2\sigma_0 m^2a^2\left[\cos^{-1}\left(\cfrac{1}{m}\right) + \frac{1}{m^2}\sqrt{m^2 - 1}\right]
</math>
The compression due to Dugdale adhesion is
:<math>
  d^D = -\left(\cfrac{2\sigma_0 a}{E^*}\right)\sqrt{m^2-1}
</math>
and the gap at <math>r = c</math> is
:<math>
  h^D(c) = \left(\cfrac{4\sigma_0 a}{\pi E^*}\right)\left[\sqrt{m^2-1}\cos^{-1}\left(\cfrac{1}{m}\right) + 1-m\right]
</math>
The net traction on the contact area is then given by <math>p(r) = p^H(r) + p^D(r)</math> and the net contact force is <math>F = F^H + F^D</math>.  When <math>h(c) = h^H(c) + h^D(c) = h_0</math> the adhesive traction drops to zero.
 
Non-dimensionalized values of <math>a, c, F, d</math> are introduced at this stage that are defied as
:<math>
  \bar{a} = \alpha a ~;~~ \bar{c} := \alpha c ~;~~ \bar{d} := \alpha^2 R d ~;~~ \alpha := \left(\cfrac{4E^*}{3\pi\Delta\gamma R^2}\right)^{1/3} ~;~~ \bar{A} := \pi c^2 ~;~~ \bar{F} = \cfrac{F}{\pi\Delta\gamma R}
</math>
 
In addition, Maugis proposed a parameter <math>\lambda</math> which is equivalent to the Tabor coefficient.  This parameter is defined as
:<math>
  \lambda := \sigma_0\left(\cfrac{9R}{2\pi\Delta\gamma{E^*}^2}\right)^{1/3} = 1.16\mu
</math>
where the step cohesive stress <math>\sigma_0</math> equals to the theoretical stress of the Lennard-Jones potential
:<math>
  \sigma_{th} = \cfrac{16\Delta\gamma}{9\sqrt{3}z_0}
</math>
Zheng and Yu <ref name=Zheng>Zheng, Z.J. and Yu, J.L., 2007, Using the Dugdale approximation to match a specific interaction in the adhesive contact of elastic objects'', Journal of Colloid and Interface Science, 310(1), pp. 27-34.</ref> suggested another value for the step cohesive stress
:<math>
  \sigma_{0} = \exp\left(-\cfrac{223}{420}\right)\cdot\cfrac{\Delta\gamma}{z_0} \approx 0.588\cfrac{\Delta\gamma}{z_0}
</math>
to match the Lennard-Jones potential, which leads to
:<math>
  \lambda \approx 0.663\mu
</math>
Then the net contact force may be expressed as
:<math>
  \bar{F} = \bar{a}^3 - \cfrac{4}{3}~\lambda \bar{a}^2\left[\sqrt{m^2 -1} + m^2 \sec^{-1} m\right]
</math>
and the elastic compression as
:<math>
  \bar{d} = \bar{a}^2 - \cfrac{4}{3}~\lambda \bar{a}\sqrt{m^2-1}
</math>
The equation for the cohesive gap between the two bodies takes the form
:<math>
  \cfrac{\lambda \bar{a}^2}{2}\left[(m^2-2)\sec^{-1} m + \sqrt{m^2-1}\right] + \cfrac{4\lambda\bar{a}}{3}\left[\sqrt{m^2-1}\sec^{-1} m - m + 1\right] = 1
</math>
This equation can be solved to obtain values of <math>c</math> for various values of <math>a</math> and <math>\lambda</math>.  For large values of <math>\lambda</math>, <math>m \rightarrow 1</math> and the JKR model is obtained.  For small values of <math>\lambda</math> the DMT model is retrieved.
 
=== Carpick-Ogletree-Salmeron (COS) model ===
The Maugis-Dugdale model can only be solved iteratively if the value of <math>\lambda</math> is not known a-priori.  The Carpick-Ogletree-Salmeron approximate solution <ref name=COS>Carpick, R.W. and Ogletree, D.F. and Salmeron, M., 1999, ''A general equation for fitting contact area and friction vs load measurements'', Journal of colloid and interface science, 211(2), pp. 395-400.</ref> simplifies the process by using the following relation to determine the contact radius <math>a</math>:
:<math>
  a = a_0(\beta) \left(\cfrac{\beta + \sqrt{1 - F/F_c(\beta)}}{1 + \beta}\right)^{2/3}
</math>
where <math>a_0</math> is the contact area at zero load, and <math>\beta</math> is a transition parameter that is related to <math>\lambda</math> by
:<math>
  \lambda = -0.924 \ln(1-1.02\beta)
</math>
The case <math>\beta = 1</math> corresponds exactly to JKR theory while <math>\beta = 0</math> corresponds to DMT theory.
For intermediate cases <math>0 < \beta < 1</math> the COS model corresponds closely to the Maugis-Dugdale solution for <math>0.1 < \lambda < 5</math>.
 
== See also ==
{{colbegin|4}}
* [[Adhesive]]
* [[Adhesive bonding]]
* [[Adhesive dermatitis]]
* [[Adhesive surface forces]]
* [[Bearing capacity]]
* [[Bioadhesives]]
* [[Contact dynamics]]
* [[Dispersive adhesion]]
* [[Electrostatic generator]]
* [[Energetically modified cement]]
* [[Frictional contact mechanics]]
* [[Galling]]
* [[Goniometer]]
* [[Non-smooth mechanics]]
* [[Plastic wrap]]
* [[Shock (mechanics)]]
* [[Signorini problem]]
* [[Surface tension]]
* [[Synthetic setae]]
* [[Unilateral contact]]
* [[Wetting]]
{{colend}}
 
==References==
{{reflist|2}}
 
== External links ==
* [http://ntrs.nasa.gov/archive/nasa/casi.ntrs.nasa.gov/19970025228_1997043322.pdf]: More about contact stresses and the evolution of bearing stress equations can be found in this publication by NASA Glenn Research Center head the NASA Bearing, Gearing and Transmission Section, Erwin Zaretsky.
 
* [http://www.amesweb.info/HertzianContact/HertzianContact.aspx Online engineering calculator to calculate contact stress for spherical and cylindrical contacts.]
 
{{DEFAULTSORT:Contact Mechanics}}
[[Category:Bearings]]
[[Category:Mechanical engineering]]
[[Category:Solid mechanics| ]]
 
[[bg:Механика на контакта]]
[[de:Kontaktmechanik]]
[[es:Mecánica de contacto]]
[[fr:Mécanique des contacts]]
[[hr:Mehanika kontakta]]
[[ja:ヘルツの接触応力]]
[[ru:Механика контактного взаимодействия]]
[[sr:Механика контакта]]
[[sv:Kontaktmekanik]]
[[zh:接触力学]]

Revision as of 07:10, 24 February 2014

Andrew Berryhill is what his spouse loves to call him and he completely digs that title. Her family members lives in Alaska but her husband wants them to move. She is really fond of caving but she doesn't have the time lately. Office supervising is my profession.

My homepage :: telephone psychic (fashionlinked.com)