Sneddon 1965
Sneddon 1965
Sneddon 1965
IAN N. SNEDDON
Abstract-A solution of the axisymmetric Boussinesq problem is &rived from which are deduced simple
formulae for the depth of penetration of the tip of a punch of arbitrary profile and for the total load which
must be applied to the punch to achieve this penetration. Simple expressions are also derived for the
distribution of pressure under the punch and for the shape of the deformed surface. The results are illustrated
by the evaluation of the expressions for several simple punch shapes.
1. INTRODUCTION
THE problem of determining (within the terms of the classical theory of elasticity) the
distribution of stress within an elastic half space when it is deformed by the normal pressure
against its boundary of a rigid punch seems to have been considered first by Boussinesq [l].
Using the methods of potential theory Boussinesq derived a solution of the problem corres-
ponding to the case of penetration by a solid of revolution whose axis was normal to the
original boundary of the half space, but the form of his solution did not lend itself to
practical computations and partial numerical results based on his solution were derived
only in the cases of a flat-ended cylindrical punch [2] and a conical punch [3].
After the publication of Boussinesqs solution several alternative solutions were
derived, an excellent account of which is given in Galins book [4]. From among these
we may mention (because they are particularly relevant to the present paper) the solutions
due to Harding and Sneddon [5] and Segedin [6].
The solution due to Harding and Sneddon uses the theory of Hankel transforms to
express the axisymmetric solution of the equations of elastic equilibrium in terms of Hankel
transforms of an arbitrary function and then to determine this arbitrary function by using
a solution due to Titchmarsh [7] of the dual integral equations to which the mixed boundary
problem can be reduced. The solution can be derived for the general case in which the
equation of the punch in cylindrical coordinates with origin at the tip of the punch is
z=w(p). The solutions appropriate to the cylindrical punch and the conical punch have
been discussed in full in [8] and [9] respectively. In each case expressions are derived for D,
the total depth of penetration of the tip of the punch, and for P, the total load on the punch,
in terms of a the radius of the circle of contact between the punch and the elastic solid, but
any attempt to express D and P in terms of an arbitrary function f is made difficult by the
complicated form of the Titchmarsh solution of the dual integral equations.
t This work was sponsored by the U.S. Joint Services Advisory Group under contract No. AF 49(638)
-1159 with the Applied Mathematics Research Group, North Carolina State, Raleigh.
47
48 IAN N. SNEDDON
In the present paper the solution of the Boussinesq problem is again derived in the form
of Hankel transforms of a function which is determined in terms of W(P) by the same pair
of dual integral equations, but, instead of using Titchmarshs solution, we make use of an
elementary solution [lo]. This enables us to derive the expressions
D= s
l f (x)dx
0J(1-x2)
P=4W--tt)-
s
1x f +)dx
OJ(1_x2)
where p and q are respectively the rigidity modulus and the Poisson ratio of the material
of the half space and the function f is defined by the relation w(p) =f(p/u).
Simple expressions
(equations (3.4) and (5.1) below) are also derived for the distribution of pressure under the
punch and for the shape of the deformed boundary of the half space.
In Section 6 the results corresponding to five special shapes of rigid punch are deduced
from the general formulae, and in Section 7 the results corresponding to the case in which
w(p)=Zc,,p are deduced. The expressions obtained for D and P in the latter case are in
agreement with those derived by Segedin [6].
The Boussinesq problem can be solved by means of a systematic use of Hankel trans-
forms and the theory of dual integral equations.
The boundary conditions are
in which the functionfis prescribed by the fact that, referred to the tip as origin, the punch
has equation z=f(p/u) so that f(O)=O;a is the radius of the circle of contact and D is a
parameter (as yet unspecified) whose physical significance is that it is the depth to which
the tip of the punch penetrates the elastic half space.
It is easily shown (see, e.g. [ll] p. 452 et seq.) that the field specified by the equations
u&, z) = 0 and
U,(P, 2) = - (2.3)
%(P,z)= - a %,[{2(1
2(1-I))
-~)+Tz}5-~(~~)e-5; <-*PI. (2.4)
is a possible displacement field. In these equations rl denotes the Poisson ratio of the
material of the half space, the function + is arbitrary and we have used the notation
*y UK 2) ; <+PI = ; MC WY 63) dt
s
to denote the Hankel transform of order v of the functionf(& z) with respect to the variable
Inserting this form for the displacement vector (u,, u+, u,) into the stress-strain relations
we obtain the equations
The relation between load and penetration in the axisymmettic Boussinesq problem 49
4=
~pz(P* (2.5)
for the z-components of the stress tensor; the remaining component a+=vanishes identically.
In these equations p represents the rigidity modulus (second Lame constant) which in terms
of Youngs modulus E is given by the equation p=.E/2( 1+q).
From equations (2.4), (2.5), (2.6) we obtain the boundary values
~,(P,O)=~,[s-l$(I); i-4
where x=p/u, so that the boundary conditions (2.1), (2.2) will be satisfied if we can find a
function +([) such that
Jl(C)=
~~t~~os(~~)dt (2.9)
s
we find that equation (2.8) is automatically satisfied and that equation (2.7) is equivalent to
the Abel integral equation
x xWdt
1
1
and since for any arbitrary function g(c) we have the relation
1 d
; dT;x~1C5-s(r); x]=&&(S); xl
we find that
1
1
If we substitute this expression into the right hand side of equation (3.2) we have the formula
~,,(P,0) = L
ax(l-q)
-
d
s
l
dx x&2-x2)
Now if x(t) is differentiable in the neighbourhood of t= 1 we see that this result can be
written as
x(l)=0 (3.5)
$x(t)=D-t
t
We can use this criterion to determine the constant D by writing equation (2.10) in the form
s f(x)dx
oJ(t2-x2
(3.6)
The relation between load and penetration in the axisymrnetric Boussinesq problem 51
D= fWdx (3.7)
f y/(1-x2)
0
The total load P on the punch required to produce the above penetration is given by the
equation
P= -2~ po,,(p,O)dp
s0
x(t)dt. (4.1)
Now if we combine equations (2.10) and (2.7) we find that x(t) is given by the equation
f(x)dx
(4.2)
o&-x2)
If we substitute this expression into equation (4.1) and carry out the integration we obtain
the formula
(4.3)
by means of which we can calculate the total load P necessary to produce the depth of pene-
tration (3.8) by a rigid punch whose profile is defined by the function f.
It is also of interest to know the form of u,(p, 0) when psa. We derive this from the
relation
1
1
t&O)=%', 5-l ~(~)~os(~~)d~; i-+x ,
t s 0
52 IAN N.SNEDWN
l xtt)dt
uz(ax,O)=
JoJ(X2_t2)> x>l.
If we substitute the expression (4.2) for x(t) into this equation and use the fact that
(5.1)
we find that
J YfWdy
1
u,(ax,O)=:sin-(l/x)+fJ(x2-1) (5.2)
0w - Y2>J(1 - Y2)
where D is the depth of penetration of the tip given by the equation (3.7).
In any practical problem it may be easier to calculate the normal component of the
surface displacement from the pair of formulae (5.1) and (4.2) than from the single formula
(5.2).
p=4CLaD (6.1)
1-q
Similarly from equation (3.4) we find that the distribution of pressure under the punch is
given b,y the equation
and from equation (5.1) we find that the shape of the deformed boundary is given by the
relation
These expressions are in agreement with known results; see, for example, [8] or (pp. 460461
of [ll]).
(b) Conical punch
For normal penetration by a cone of semivertical angle a we may take J(x) =6x where
E=a tan a.
From equation (3.7) we find that the vertex of the cone penetrates to a distance
D= $~E=~Iw tan a and from equation (4.3) that the total load necessary to effect this
penetration is given by the equation
7t/lll2
P= -tancc, (6.4)
1-V
which can be written in the alternative form
4pcotaD,
P=- . (6.5)
n(l -V)
From equation (2.10) we find that in this case
x(t)=+r).
Substituting this expression into equation (5.1) we find that the deformed shape of the free
surface is given by the relation
20
u,(p, 0) = G a sin- (a/P)-P+J(P2-a2) 9 p>a, (6.6)
{ I
and substituting it into equation (3.4) we find that the distribution of pressure under the
punch is given by the relation
2crD
b&O)= -(I _tl)na cash- (alp) 5 O<p<a. (6.7)
The equations (6.5), (6.6), (6.7) are in agreement with the results obtained by other methods,
see, for instance, [3], [9] or (pp. 462-466 of [l I]).
(c) Punch in the form of a paraboloid of revolution
If the punch is a paraboloid of revolution with equation p2 = 4kz we may take f (x) =&x2
where E= a2/4k. If we substitute this expression into equation (3.7) we find that the relation
between D and a is
a2=2kD, (6.8)
so that the depth of penetration of the paraboloid is twice the distance that the circle of
contact lies below the original boundary of the half space.
In a similar way equation (4.3) leads to the expression
p= 16w
3(1--r])
for the total load. In terms of the total penetration D this may be written in the form
so that we find from equation (3.3) that the distribution of stress under the punch is given
by the equation
(6.10)
and from equation (5.1) that the shape of the free surface is given by the equation
t _fWdx a2 xdx
s e J(t2 -x2> = s oJ~(t2-~2)(RZ-a2~2))
if we change the variable in the integral on the right from x to u= t-J(t -I-x2) we find that
t
s f(x)dx
,~(t2_~2)=~ut10gR_at
1 R+at
Putting t= 1 in this result and then substituting into equation (3.7) we find that the relation
(6.12)
R+a
D=ialog - . (6.13)
R-a
Substituting from equation (6.12) into equation (3.6) we find that in this case
2D at R+at
X(t)=~--pgRat. (6.14)
Using the fact that
1
R+atdt
tiog- 2
s 0 R-at 2
we see from equation (4.1) that the total load P on the sphere necessary to produce the
penetration (6.13) is given by the equation
-aR . (6.15)
The relation between load and penetration in the axisymmetric Boussinesq problem 55
~~q)B
p=(l [
(a2+f12)logp
a
-2aj?] . (6.17)
Segedin [6] has considered the case in which the profile of the punch has equation
z= 5 c,pn (7.1)
n=l
in which case we have
f(x) = zl cax .
If we substitute this expression into equation (3.7) we find that the appropriate formula for
the total penetration of the depth of the punch is
(7.2)
and if we substitute it into equation (4.3) we find that the total load P is given by the formula
pz2Jnpu m nT(+n+l)
-lx
1-q =I r($n++) CunY
(7.3)
(7.4)
Hence from equation (3.4) we find that the formula for the distribution of pressure under the
punch is
* 1)
-c2P
azz(p~o)=-a(l-tj),=l
nWn+
I-@++)
q&U4 , (x=pla, O<x<l), (7.5)
56 IAN N. SNEDDON
k(x)=
s l t-dr
xJ(tZ_x2) (n>l, O,<x<l),
and from equation (5.1) that the shape of the deformed surface is given by the relation
(7.6)
j,(x) =
s 1
tdt
0 J(x-- P)
The integrals i,(x) can be determined by the recursion relation
(n>l, x>l). (7.8)
REFERENCES
PI J. B~u=NE.w,Applications aks Potenticls d I&& de l&uiiibre et du Mouvement aks Soiidcs &stiques,
Gauthier-Villars, Paris, (1885); also see 1. TODHUNTERand K. PEARSON,History of the ?%eory of
Elasticity, 2,237 et seq. Cambridge University Press, (1893).
[2] A, E. H. LOVE,Phil. Tram. A. 22B, 377 (1929).
[3] A. E. H. LOVE,Qmt. J. Math. (Oxfod) 10, 161 (1939).
[4] L. A. GALIN,Ko~t~~e -&c&hi Teorii ~pr~osti. Gas. Izdat. Tech. -Tear. Liter., Moscow (1953).
[S] J. W. HARD~NQ and I. N. SNED~, Proc. Cam. Phi/+ Sot. 41, 12 (1945).
[6j C. M. SEOEDIN,Mat~~ti~ 4,156 (1957).
[7] E. C. TITCHMARSH, Introduction to the 77reory of Fourier Integrals, Clarendon Press, Oxford (1937).
[S] I. N. SNED~~N,P*oc. Cam. Phil. Sot. 44,29 (1946).
[9] I. N. SNEDDON,Proc. Cam. Phil. Sot. 44,419 (1948).
[IO] I. N. SNEDDON, Proc. Gfasgow Math. Assoc. 4, 108 (1960).
[ll] I. N. SNLDDON, Fourier lkmforms. McGraw-Hill, New York (1951).
R&m&---Lauteur a &abli une solution du prob&me axisym#rique de Boussirmq qui lti B penuis de
dtduire des form&s simples dormant la profon&ur de p&&ration duu p&&rateur de protil arb&raim
ainsi que la charge totale &eas&re pour assurer c&e caption. 11dome &galenxmtdes expm&ms
simples qui d&nis.snt la ~~~n de la pression sous le p&trateur aiusi que la forme de la svrfaoe
&form&e.
Lea r&Bats sont i&&r& par lapplication de tea expressions aux cas de pimieurs p&&trateursde
forme simple.
The relation between load and penetration in the axisymmetric Boussinesq problem 57
Sommari+Si deriva una soluzione de1 problema assisimmetrico de1 Boussinesq dalla quale si deducono
semplici formule per la profondita di penetrazione della punta di un punzone di protilo arbitrario e per il
carico totale the deve venire applicato al punzone per ottenere detta penetrazione. Si derivano inoltre
semplici espressioni per la distribuzione della pressione sotto il punzone e per il prolilo della super6cie
deformata. I risultati sono illustrati con la valutazione delle espressioni per vari profili semplici di punzone.