2.2. The Case of Rebounds
The four cases analyzed in the first paper assumed that the coin is not rebounding after impact. This is obviously the case by definition for the first two cases with inelastic surfaces, the velocity component normal to the landing surface becomes nil, and the coin does not rebound.
For the last two cases, where elastic surfaces are involved, one has to be more cautious. A simple condition of no-rebound can be found by writing that the normal impulse of the landing surface at impact must be smaller than the integral over time of the coin weight considered as an impulsive force during the short period of impact, i.e.,
For this condition to be respected in the two cases of elastic surfaces, the duration of impact
must be such as:
where
is the coefficient of restitution. Relation (6) yields, respectively, for the cases of elastic and perfectly smooth bodies and of elastic and partially rough bodies (with
μ as the coefficient of friction)
where:
For the numerical values given in (1) and assuming, for a 50 Euro cent coin, coefficients of restitution = 0.5 and of friction = 0.05, (7) and (8) yield impact duration longer than several seconds, which is practically impossible for usual surfaces. This shows that the condition (5) cannot be practically fulfilled, and the coin will rebound.
In this case, one can still calculate the probability that the coin will stay on its edge at the second impact by considering that the first impact will be such as to deliver the favorable initial conditions (Event 1) for the second impact to result in the coin staying on its edge (Event 2). The probability of the first event
is obviously independent of the probability of the second event
. The opposite is, of course, not true and the probability of the second event
depends on the first event
having occurred. Therefore, the probability that the coin will stay on its edge at the second impact, knowing that the first impact has delivered the favorable conditions for the second event to occur, is:
where indexes 1 and 2 refer, respectively, to the first and second impacts.
With the following hypotheses:
- -
between the first and second impacts, the coin attains a maximum height:
- -
the rotation velocities and moments of inertia of the coin, respectively, before the second impact and after the first impact are the same, yielding:
- -
the vertical components of the coin velocity before and after the second impact are, respectively:
- -
the coin and landing surface are elastic and partially rough bodies, such as the normal and tangential impulses read:
and the impact equation reduces to:
yielding:
where:
The coin will stay on its edge if the angular velocity after the second impact
is nil while
. Solving the Equation (19)
= 0 for
yields:
under the condition that the denominator of (19) is different from zero, i.e., for:
Note that this condition can be extended to the other configurations of the position of the instantaneous axis of rotation after impacts, as was discussed for (66) in [
1]. This would introduce additional numerical coefficients in the expressions of the moments of inertia, but would not fundamentally alter the value of the angle
.
For
to be real, the condition on the root in the numerator of (21) reads:
For the condition
to hold, i.e.,
the argument of the arctan function in (21) must be positive and smaller than
, i.e.,
For the positive condition, the numerator and the denominator must be of the same sign. Considering the negative sign in front of the numerator root in (25):
- -
both denominator and numerator are positive if
that combined with (23) yields:
which is true only for a positive sign in front of
, i.e., a coin counterclockwise rotation;
- -
both denominator and numerator are negative if:
which includes condition (23), and is true either for a positive sign in front of
, i.e., a coin counterclockwise rotation after the first impact, as long as values of
fulfil the condition (28), or for all of the negative values of
, i.e., a coin clockwise rotation.
For the positive sign in front of the numerator root of (21), the denominator must be positive, which is true if condition (28) holds, meaning that the positive sign in front of must be chosen, i.e., a coin counterclockwise rotation. However, the solution, in this case, yields a too large value of .
The other part of the condition, i.e., the argument of the arctan in (21) smaller than
yields, as in relations (72) to (74) in [
1]:
These three conditions include condition (23), as the right part of (23) is always greater than the right parts of (29) to (31).
Summarizing the conditions for this case:
- -
if , then (29) includes (27) and must be positive, i.e., a coin counterclockwise rotation after first impact;
- -
if , then (30) is the limiting case of (26) and (27) and is positive, i.e., a coin counterclockwise rotation after first impact;
- -
if , then (31) includes (27) and is positive, i.e., a coin counterclockwise rotation after first impact;
- -
if , then (31) includes (28) and can be either positive or negative, i.e., a coin counterclockwise rotation or clockwise rotation after first impact.
The three conditions (29) to (31) constrain the favorable values of
, which can be translated into conditions on
through (63) of [
1] (where the index
has been added to
,
and
).
Assuming that
, the left and right parts of (29) yield, respectively:
with
and:
The analytical solutions of (32) and (33) for
require solving analytically fourth-degree equations in
, which is not an easy task. However, the numerical solutions can be found to define the range
of the allowed values of
at first impact to deliver favorable conditions. Relations (32) and (33) yield numerical values of, respectively,
and
, involving functions
and
of
,
,
,
,
and
, respectively, through
and
and through
and
:
The probability that the first impact delivers favorable conditions for the second impact to bring the coin on its edge reads, then:
For this case
, the condition (29) constrains the favorable coin rotation velocity
after the first impact between the minimum and maximum values
and
corresponding to
and
through (63) of [
1]:
Replacing these two values
and
in (21) yield the two limiting values
and
:
with:
and a similar expression for
with
replacing
.
The probability that the coin stays on its edge at the second impact, knowing that the first impact has delivered the appropriate favorable conditions is then:
where the formula
was used.
The total probability that the coin stays on its edge at the second impact, after the first rebound has delivered the appropriate conditions, is given by (9).
Assuming coefficients of restitution
= 0.5 and of friction
= 0.05 for a 50 Euro cent coin yields that
. One considers further the two extreme cases of the position of the instantaneous rotation axis after the first impact, i.e., respectively, along the horizontal axis passing first through the coin center of mass, and second, through the coin contact point with the landing surface (see discussion leading to (66) in [
1]), yielding values of
varying between 2.955 × 10
−3 and 1.475 × 10
−2, while
= 2.955 × 10
−3 under the initial hypothesis.
The condition (29) imposes a quite stringent range of values for the ratio
, namely:
respectively, for
= 2.955 × 10
−3 m and
= 1.475 × 10
−2 m, yielding, respectively, for
and for
m:
The allowed range for the coin angular velocities after the first impact thus increases for the increasing initial height .
Within the value ranges (1) of initial parameters, the largest value of the probability is obviously attained for the largest values of and . The largest value of is obtained for the largest value of and the smallest value of . is the largest for and , i.e., for = 0 and = 0, while is the smallest for and , i.e., for = 10π rad/s, = 5 m/s, and .
For the first case of = 2.955 × 10−3 m, a series of coin tosses having initial conditions in the ranges (1), one finds, successively:
= 0.0101, = −0.0588 (rad)
= 0.0950, = 0.0261 (rad)
= 190.311, = 8.616 (rad/s)
= −0.0849, = −0.0150 (rad)
= 2.193 × 10−2, = 2.225 × 10−2, = 4.880 × 10−4
or approximately a throw every 2050 throws.
For the second case of = 1.475 × 10−2 m, one has similarly:
= 0.0099, = −0.0588 (rad)
= 0.0948, = 0.0261 (rad)
= 38.135, = 1.727 (rad/s)
= −0.0849, = −0.0150 (rad)
= 2.188 × 10−2, = 2.225 × 10−2, = 4.868 × 10−4
or approximately a throw every 2054 throws.
For more common values, a series of throws from an initial height = 1.5 m with a velocity = 1 m/s under an angle of and an initial rotation velocity varying between 0.5 and 5 turns/s yields, successively, first for = 2.955 × 10−3 m:
= 0.0829, = 0.0585 (rad)
= 164.814, = 68.840 (rad/s)
= −0.0849, = −0.0150 (rad)
= 7.774 × 10−3, = 2.224 × 10−2 and = 1.729 × 10−4
or approximately a throw every 5783 throws, and second for = 1.475 × 10−2 m:
= 0.0829, = 0.0586 (rad)
= 33.026, = 13.794 (rad/s),
= −0.0849, = −0.0150 (rad),
= 7.759 × 10−3, = 2.224 × 10−2 and = 1.726 × 10−4
or approximately a throw every 5794 throws.
For the sake of completeness, for all of the above values, the term under the root sign in condition (22) is negative, yielding imaginary values for the argument of the arctan function, which shows that condition (22) is fulfilled for all of the practical values of .
The cases of successive rebounds can be treated similarly, albeit with more and more complexity in the various relations.