CC BY
58
Russian Journal of Biomechanics
www.biomech.ac.ru
ON INFLUENCE OF THE RESISTANCE CRISIS ON THE FOOTBALL
FLIGHT DYNAMICS
R.N. Rudakov, V.A. Lochov, A.E. Fedorov
Department of Theoretical Mechanics, Perm State Technical University, 29a, Komsomolskii Prospect, 614600, Perm, Russia, e-mail: rrn@theormech.pstu.ru
Abstract. It is described the mathematical model of football movement when employing standard positions (kick from the center of the football ground and corner kick). An important feature of the model is an account of dependence of aerodynamic coefficients of drag and lift forces on the ball velocity. It is found the optimal conditions which provide a football hit a rival goal for the minimal time. Also the problem of the impact theory is considered which allow to describe the football boot kick on the ball. It is proposed the method of determination of aerodynamic coefficients in three-dimensional movement of a ball.
Key words: biomechanics in sports, model of football movement, aerodynamic coefficients, optimal solution
1. Introduction
The major role in football belongs to the trajectory of a ball for a correct pass or to score a goal. Often realization of so called standard positions leads to success. The trajectory of a ball represents a three-dimensional curve which depends on the initial conditions of ball's movement. The aim of this work is the analysis of optimal initial conditions at which a ball hits the goal for the minimal time. In literature on sports biomechanics there are papers on mathematical modeling of movement of volley-ball and tennis-ball [1] - [4]. In this work a theoretical research of a football's flight phase is made and several problems are solved which are connected with execution of standard positions: kick from the center of the football ground and corner cock. An important feature of the model is an account of dependence of aerodynamic coefficients of drag and lift forces on the ball velocity. This account allows to receive results which correlate with experiments. Also the problem of the impact theory is considered which allows to describe the football boot kick on the ball. It is proposed the method of determination of aerodynamic coefficients in three-dimensional movement of a ball. These results have not only scientific but also an applied importance. They can be useful for coaches and football-players in training of standard positions. For clear understanding of obtained results the scheme of the football ground is presented (Fig. 1).
2. System of forces acting on a football. Aerodynamic coefficients of drag and lift forces
In this work we consider flat and spatial trajectories of the ball movement. In Fig. 2 a flat trajectory of a ball movement in vertical plane Oxy at the initial velocity V0, lying in
plane Oxy is represented. The rotation of the ball takes place around the axis parallel to axis z .
The important kinematical parameters of a movement are: the velocity vector of translation V directed at a tangent of the center of mass trajectory and the angular velocity vector « directed along the axis of rotation to that side from where the rotation is seen counter-clockwise. The rotation can be arbitrary depending on the kick on the ball.
During the ball flight there are forces acting on the ball: P = mg - the force of gravity, where m - the mass of the ball, g - the acceleration of the free fall; R - the drag and Q - the lift force. The forces R and Q depend on velocity squared and are defined as follows:
Fig. 2. The system of forces acting on a football at the reverse rotation.
R = 2 p CDSV2, (2.1)
Q = ± p CLSV 2, (2.2)
where p - the air density, S -the area of the largest section perpendicular to the velocity vector V (middle section), CD and CL - the aerodynamical coefficients: drag coefficient and lift coefficient. On account of the small magnitude the friction force moment with respect to the center of mass was ignored. In this case the angular velocity vector is constant in time during the flight. The force of gravity P is directed vertically down. The drag force R is opposed to the velocity vector V :
R = - R V. (2.3)
VI
It is more difficult to define the direction of the lift force Q due to Magnus effect. There is no theoretical substantiation of this effect but the main mechanism of the lift force occurrence lies in the fact that in the flight of the ball the contrary air flow interacts with the air flow due to rotating ball. In mixture zone the domain of elevated pressure forms. For example, in the reverse rotation of a ball the velocity of the lower point of the ball coincides with the direction of V . Therefore the mixture zone is located in the vicinity of the lower part of the ball where the domain of elevated pressure arises. The term "lift force" is used also under other directions of the lateral aerodynamic force. In experiment the lift force is measured as a force directed perpendicular to vectors œ and V. Also from experiment it is known that force is directed to that side from where the shortest turn of vector œ to vector V is seen counter-clockwise. This allows to use the rule of cross product of two vectors [4]:
Q = Q . (2.4)
œ x V
The central part of problem of investigation of solid movement in continuum is determination of aerodynamic coefficients. The drag coefficient for the sphere without rotation was found experimentally in large range of velocities. In paper [5] the drag coefficient CD dependence of Reynolds number Re (Re=Vd/v, d - the diameter of a ball, v -the kinematic air viscosity) was investigated. This dependence is shown in Fig. 3. In this diagram it can be seen that at Re values from 1000 to 140000 (it corresponds to football velocity 0.1 ^14m/s) drag coefficient CD is almost constant and it is equal to 0.45. In football the velocity changes in the large range, thereby the account must be taken of coefficient CD dependence of the great values of Reynolds number. As can be seen from Fig. 3 the sharp decrease of coefficient CD begins at value log Re > 5.4 - "resistance crisis". This effect allows to increase significantly the ball flight distance at increase of the initial velocity.
In Fig. 4 CD and CL dependence of the angular velocity of ball rotation for different Reynolds numbers (up to Re = 105 ) is represented [6]. It is clear that at n = > 1 (ratio of
the rotational velocity to the velocity of translation) coefficients CD and CL depend weakly on angular velocity and Reynolds number. According to Fig. 4 it is possible to take CD = 0.6 and CL = 0.35 at the fast ball rotation.
log CD
2.75 2.25 1.75 1.25 0.75 0.25 -0.25 -0.75 -1.25 -1.75
-1.5
-0.5
0.5
1.5 2.5
log Re
3.5
4.5
5.5
6.5
Fig. 3. Drag coefficient dependence of Reynolds number.
0.6
0.4
0.2
3 ""1
2 /2
3^ Cl "1
n
Fig. 4. Drag coefficient CD and lift coefficient tt dependence of the ratio of the rotational velocity to the velocity of translation: 1 - Re = 6.15-104, 2 - Re = 7.74-104, 3 - Re = 10.70-104.
3. Differential equations of the ball center of mass movement. Cauchy problem
By the theorem of the center of mass movement
m ac = P + R + Q, (3.1)
where ac - the ball center of mass acceleration (ac = dV), P, Q, R - forces described in
c v c dt'
section 2. By projecting the left and right hand sides (3.1) on system axes we can receive with the help of (2.1) - (2.4) the system of equations of center of mass movement (for convenience the projections of center of mass velocity V on axis x, y, z (Fig. 2) and denoted as u, v, w):
mdu = -Ru + Q X V)x , dt V ^ x v|
mf = - Rv + Q ^ - mg,
dt V ^ xV|
mdw = -rw + Q (^XV>L,
dt V ^ xV|
u = dx v = , w = dz, (3 .2)
dt dt dt y J
V = -\l u2 + v2 + w2 ,
where R, Q are defined by (2.1) and (2.2), x, y, z - the ball center of mass coordinates. The
2 3 d2 given parameters of the problem are: g = 9.8 m/s , p = 1.23 kg/m , m = 0.4kg, S = k—
(d = 0.22 m). In the case of the central impact CD = CD (Re), CL = 0, at the eccentric impact CD = 0.6, CL = 0.35. For the initial-value problem it is necessary to preset coordinates and velocity of the center of mass at the initial time moment:
x — x 0, u — u 0,
t = 0: y = y 0, v = ^ . (3.3)
z = ^ w = ^
Cauchy problem (3.2) - (3.3) was solved numerically using the method of step-by-step integration (Euler's method). The integration step was chosen so that calculation result was independent of the step parameter (At = 0.001 s).
4. Ball trajectories calculation at standard positions
For a finding optimal ball trajectories it is necessary to solve a problem of optimization. At kick on a ball we may vary three parameters: the module of the ball initial velocity V0 , the angle between the plane of a ground (Oxz) and the vector of the initial velocity a and the angle between the plane of a goal (Oyz) and the vector of the initial velocity P (Fig. 1). The aim of this optimization is reception of ball trajectories, moving on which the ball will hit the goal under the goal crossbeam for the minimal time. Then the problem is written as:
t (V0, a, p) ^ min, (4.1)
on condition that
H < y(z = L) < H, (4.2)
where H - is the goal height minus the ball radius, H - is the minimal allowable height of a ball in the goal. Also it is necessary to impose restriction on the maximal speed of the kick:
V < Vmax, (4.3)
were Vmax - is the maximal allowable initial ball velocity.
4.1. Center kick
It happened several times that a football-player kicked a ball from the center of a ground and scored a goal. Basically it was goalkeeper's fault in most occasions, but forward's skill is the main thing. All football funs know that kick of Pellet, when he scored a goal to Bulgaria right from the center of the football ground.
We calculated the possible trajectories of a ball and founded an amazing fact. If drag coefficient CD is constant and equal to 0.45 (as usually accepted), the ball initial velocity should be more than 50 m/s to reach the goal. Even the strongest players can receive hardly such velocity.
Test of CD coefficient dependence of Reynolds's number (Fig. 3) shows that so-called drop of resistance starts at Re* = 3.5 -105 i.e. value of aerodynamic drag coefficient slumps. The physical reason of this event is that air flows around a ball moving at a high speed separates from the ball surface closer to the rear part. At the same time the eddy diameter decreases and improves the ball streamlining. For football it starts at the velocity:
Re*v
K =
d
(4.4)
the ball diameter d = 0.22 m, so
Since air the kinematic viscosity v = 0.15 • 10 4 critical speed V = 14m/s. We consider CD dependence of the ball velocity. This dependence is shown in Table 1.
Table 1.
V0 m/c 14 15 16 17 18 19 20 25 30 35 40 45
0.45 0.41 0.38 0.35 0.28 0.2 0.14 0.12 0.09 0.09 0.1 0.095
If the velocity is under 14m/s, then CD = 0.45. This is a result of data processing in
Fig. 3.
The optimal ball trajectory is shown in Fig. 5. We used Table 1 data to calculate this trajectory taking into account CD dependence of the ball velocity. The minimal initial velocity of our ball should be 30 m/s at kick angle a = 29O to score a goal. The ball flight time is 2.7 s. By the way it took 4 seconds to score the "quickest" goal from the center of the ground. According to our calculations the ball flight time from the center to the goal is t < 3 s. Real and estimated time comparison gives a certain assurance in the exactness of our calculations.
Curve 1 in Fig. 5 represents the ball trajectory, line 2 is drawn at the goal' height. The origin of the coordinates is in the ground corner (Fig. 1). It is clear that the trajectory is a ballistic type of curve: it is gentle at the beginning and deep at the end. A goalkeeper can catch a ball staying basically in the goal range. Our calculations let us estimate a safe distance for a goalkeeper to come back to his goal. We tried to determine the height of optimal trajectory and the time of ball flight dependence of its initial velocity (Figs. 6 and 7).
y,m
2 %
a 7s
10
20
30
40
50
x, m
Fig. 5. Optimal ball trajectory at center kick K0 = 30m/s, a = 29
0
o
H, m
30
35
40
45
50 Vq, m/s
t, s
Fig. 6. Dependence of maximum ball trajectory height of its initial velocity.
30
35
40
45
50
x, m
Fig. 7. Ball flight time in dependence of its initial velocity.
4
3
2
1
y, m
0 B 1 %
2 r-
0 10 20 30 40 50 X, m
Fig. 8. Ball trajectory at center kick Vq = 45m/s, a = 11o.
A dependence of maximum ball trajectory height of its initial velocity is shown in Fig. 6. Figure 7 represents a ball flight time in dependence of its initial velocity. Our measurements have shown that at height H > 10 meters the ball flight time is approximately 3 seconds and these results correspond to our calculations. If we increase the initial kick velocity up to V0 = 45 m/s (as follows from Fig. 7), the ball flight time decreases more than twice (see the trajectory in Fig. 8). In this case the trajectory is more gentle, its height is 3.5 m and the flight time is 1.4 sec.
It is obviously that the initial velocity increase from 28 m/s to 35 m/s reduces the flight time in 1.8 times, so a possibility to score a goal increases also as compared to computational results. We used records of TV football match to measure the flight time and estimate the trajectory height of a ball kicked from the goal. According to our analysis the maximum velocity can reach 35 m/s. Although it is possible that the strongest football-players could impart to a ball the higher velocity.
4.2. Corner kick
In this problem we calculated the possible ball trajectories for corner kick. In our calculations (as follows from Fig. 4) we assumed that the drag coefficient CD = 0.6 and the lift coefficient CL = 0.45 .
The ball trajectories in vertical and horizontal projections are shown in Fig. 9 and Fig. 10. The minimal initial velocity should be 30 m/s to score a goal. So big value of velocity is explained by the rotation of a ball, thereby drag coefficient CD is increased. Also the trajectory is extended due to a deviation of a ball from the vertical plane. According to Fig.9 the maximum flight height is about 6 m and according to Fig.10 the maximum ball deviation from the goal line is about 2 m.
y, m 10
5
V0
0 5 10 15 20 25 30 z, m
Fig. 9. Ball trajectory at corner kick in goal plane projection (Oyz), Vq = 30m/s, a = 30o, /3 = 19o
z, m 5
x, m
Fig. 10. Ball trajectory at corner kick in horizontal plane projection (Oxz), VQ = 30m/s, a = 30o, / = 19o
(top view).
5. Eccentric kick on a ball
Using an eccentric kick a football-player can generate rotary and translational motion. It is necessary to examine an impact interaction of player's foot and a ball to determine the initial velocity. Let us assume that the velocity of football boot does not change at impact (player's foot mass exceeds ball's mass many times over) and there is no slipping of football boot on a ball.
A ball is at rest before impact. The ball's mass center after impact acquires velocity U and the ball also starts the rotation around the mass center with the angular velocity a . The angle p that determines the impact point and the boot velocity V are known.
We have to find impact angle a, the ball's mass center velocity after impact U and the angular velocity a (Fig. 11).
Let us use the impact theory to solve the task. We write down the theorems of momentum and of angular momentum for the impact:
mUx = Sx,
Ic ©= Sx r, (5.1)
where ST is the tangent component of impulse of the impact force, m is the ball mass, r is the 2
ball radius, IC = — mr2 is the moment of ball's inertia with respect to the central axis.
0
The coefficient of restitution is equal to the ratio of relative velocities of a ball after and before the impact:
U - V
k = Un--Vn, (5.2)
Vn
where Vn is the boot velocity V projection on axis n:
Vn = V cos(a + 9). (5.3)
The tangent component of velocity V :
Vx = V sin(a + 9). (5.4)
Considering no slipping of football boot on a ball the theorem of addition of velocities
yields:
Vx = Ut+q r . (5.5)
We can find all sought quantities from equations (5.1) - (5.4):
tg (a+ 9) = 2.5(1 + k) tgp, (5.6)
U = (1 + k) V cos(a + 9)
cos 9 ' '
ro= 5 r (1 + k )V cos(a + 9) tg9. (5.8)
Equation (5.5) defines angle a, (5.7) and (5.8) - translational and rotary velocities of a ball. For the central impact we have a = 0, 9 = 0. Using formulas (5.6) and (5.8) we can receive the following values for the ball's mass center velocity U and angular velocity a :
U = (1 + k)V , © = 0. (5.9)
We have got the experimental data of coefficient of restitution k for volley-ball rebound from solid surface k = 0.65 [2].
We use this value to estimate the football take-off velocity. With (possible for football player) the kick velocity V = 20 m/s we can calculate the ball's take-off velocity after impact U = 33 m/s, it is sufficiently to score a goal at center and corner kick.
6. Method of aerodynamic coefficients determination
There are no experimental data for coefficients value CD and CL at high velocity with ball rotation. According to Fig. 4 they are known up to V = 6.8m/s i.e. up to usual ball velocity. It is actually to conduct an experiment to receive those coefficients.
The velocities higher than V = 6.8m/s were used for examination of ball rotation. Reynolds's number basically does not influence on drag and lifting force coefficients at high rotation rate (see Fig. 4). Nevertheless there is a need for new experiments to determine coefficients CD and CL. We offer using methods suggested in paper [1] to determine CD for volley-ball flight and developed in paper [2] to determine CD and CL for ski-jumper flight. This model is based on object movement record analysis that determines its mass center velocities and accelerations. Differential equations of the mass center movement are used to determine drag and lift force coefficients. Only plain trajectories are considered in papers [1] and [4].
We consider the three-dimensional movement of the football. The differential equations of the center of mass movement are:
mxc = Rx + Qx,
myc = Ry + Qy, (6.1)
mzc =-mg + Rz + Qz. Substituting results of numerical differentiation of ball coordinates in (6.1) yields:
2
md^ = CDSV 2 u + £ CLSV2 ,
dt2 2 D V 2 L xV| '
m^C = -£ CDSV2 £ + £ CLSV2 -^^f, (6.2)
dt2 2 D V 2 L x V ' V '
m^-Zt = -mg-H 2 W + £ QSF 2(" X V) z
dt2 0 ^ 2^"' I« x V
dxc dy„ dzc u = —c, v = , w = ■ c
dt ' dt ' dt
V = V x? + y? + z^ .
We suggest to video-record the ball flight from two points at the right angles. This record will determine the coordinates of the center of mass as functions of time:
xc = xc (t X
yc = yc (t), (6.3)
zc = zc (t).
Thus only two unknown values CD and CL are contained in three equations (6.3). Those values could be determined from first two equations (6.3) and the third equation could be used to check the method accuracy. It should be fully identical while substituting the coefficients determined.
This paper shows only the determination method of CD and CL.
7. Conclusions
The mathematical model of the football mass center motion under the action of gravity, drag and lift forces is developed. Several standard ball positions are examined. The optimal trajectory is calculated for ball kicked from the center of the ground to score a goal. The drop of resistance is the main point to consider. Calculations without considering drop of resistance cause the results that do not correspond to practice. According to the calculations the initial velocity increase causes the essential reduction of the ball flight time. In this case ball moves basically along whole trajectory with the drop of resistance and this cause the lesser resistance. Optimal corner kick trajectory with ball rotation is calculated to aim a ball into the goal. The football kick is examined. The ball initial velocity after kick is estimated. А method of drag and lift coefficients determination is suggested. This method is based on ball flight video-records processing. The results could be used by trainers and football-players to improve performance in the standard positions.
8. References
1. DEPRA P., BRENZIKOFER R., GOES M., BARROS R. Fluid mechanics analysis in volleyball services. Proceedings of XVI International Symposium on Biomechanics in Sports. Book 2. Riehle H.J., Vieten M.M. (Editors), UVK-Universitatsverlag Konstanz GmbH, 161-164, 1998.
2. РУДАКОВ Р.Н., ШИРИНКИН Т.А., АРТЮШКОВ Н.Ю., БРАТЧИКОВ А.И., МИХЕЕВА С.А. Математическое моделирование подачи волейбольного мяча. Российский журнал биомеханики, 4(2): 50-55, 2000 (in Russian).
3. РУДАКОВ Р.Н., КАМЕНСКИХ А.В., ШУЛЬГИН П.В. Оптимизация траектории теннисного шарика. Российский журнал биомеханики, 4(2): 42-49, 2000 (in Russian).
4. RUDAKOV R.N., NYASHIN Y.I., PODGAETS A.R., LISOVSKI A.F., MIHEEVA S.A. The influence of aerodynasmical forces on the movement of sportsmen and sport balls. Russian Journal of Biomechanics, 5(2): 83-94, 2001.
5. Современное состояние гидроаэродинамики вязкой жидкости. Под ред. С. Гольдштейна. Москва, Государственное издательство иностранной литературы, т. I, 1948 (in Russian).
6. Современное состояние гидроаэродинамики вязкой жидкости. Под ред. С. Гольдштейна. Москва, Государственное издательство иностранной литературы, т. II, 1948 (in Russian).
ВЛИЯНИЕ КРИЗИСА СОПРОТИВЛЕНИЯ НА ДИНАМИКУ ПОЛЕТА
ФУТБОЛЬНОГО МЯЧА
Р.Н. Рудаков, В.А. Лохов, А.Е. Федоров (Пермь, Россия)
При игре в футбол важнейшую роль играет траектория мяча, обеспечивающая правильный пас или поражение ворот соперника. Часто к успеху приводит реализация так называемых стандартных положений. Траектория мяча представляет собой пространственную кривую, вид которой зависит от начальных условий движения мяча. Цель работы - исследование оптимальных начальных условий, при которых мяч попадет в ворота соперника за минимальное время. В литературе по биомеханике спорта известны работы по математическому моделированию движения волейбольного мяча и теннисного шарика [1] - [4]. В настоящей работе проведено теоретическое исследование фазы полета футбольного мяча и решены задачи исполнения стандартных положений - удар с центра поля по воротам и углового удара. Особенностью полученных решений является то, что учтена зависимость аэродинамического коэффициента сопротивления от скорости мяча. Это позволило получить результаты, согласующиеся с экспериментальными данными. Рассмотрен также удар бутсой по мячу и представлена методика определения аэродинамических коэффициентов лобового сопротивления и подъемной силы для общего случая пространственной траектории мяча. Полученные результаты имеют не только научное, но и прикладное значение. Они могут быть использованы тренерами и футболистами при отработке техники исполнения стандартных положений. Библ. 6.
Ключевые слова: биомеханика спорта, модель движения футбольного мяча, аэродинамические коэффициенты, оптимальное решение
Received 01 May 2002
