Оптимизация траектории теннисного шарика Текст научной статьи по специальности «Физика»

ОПТИМИЗАЦИЯ ТРАЕКТОРИИ ТЕННИСНОГО ШАРИКА

Р.Н. Рудаков*, А.В. Каменских**, П.В. Шульгин**

* Россия, 614600, Пермь, Комсомольский проспект, 29а, Пермский государственный технический университет, кафедра теоретической механики

** Россия, 614600, Пермь, Комсомольский проспект, 45, Муниципальное образовательное учреждение средняя школа №9 им. А.С. Пушкина с углубленным изучением предметов физико-математического цикла

Аннотация: При игре в настольный теннис важную роль играет вращение шарика, которое создает подъемную силу (эффект Магнуса). Зависимость аэродинамических сил от скорости поступательного движения и угловой скорости вращения шара исследовалась ранее экспериментально [1,2]. В настоящей статье показано, что в настольном теннисе параметры движения шарика согласуются с диапазоном исследуемых параметров в [1,2], что позволило найти аэродинамические силы, действующие на теннисный шарик. Составлены дифференциальные уравнения движения центра масс теннисного шарика, которые решались численно методом пошагового интегрирования. Исследованы траектории движения шарика после его приема в ближней и дальней зонах при прямом центральном ударе по шару, а также приеме накатом (верхнее вращение). Найдены оптимальные способы приема шарика и оптимальные траектории, которые обеспечивают попадание шарика на сторону стола соперника с максимальной скоростью.

Ключевые слова: настольный теннис, аэродинамические силы, оптимальные траектории.

1. Система действующих на шарик сил

В данной работе рассматриваются плоские траектории движения шарика. На рис.1 приведена траектория движения мяча в вертикальной плоскости ОХУ. Однако приведенные ниже уравнения позволяют исследовать и пространственные движения.

Кинематические характеристики движения - вектор скорости поступательного движения V, направленный по касательной к траектории центра масс, и вектор угловой скорости ю, направленный по оси 2 в случае подрези.

При полете шарика на него действуют следующие силы: сила тяжести P=mg, где

т - масса шарика, g - ускорение свободного падения; сила лобового сопротивления Я и

подъемная сила Q

Я = 2 рС^у2 (1)

Q = 2рС\СТ'2 (2)

где р - плотность воздуха, £ - площадь наибольшего сечения шара, перпендикулярная V, Си и Сь - аэродинамические коэффициенты, соответственно, лобового сопротивления и подъемной силы, момент сил трения относительно центра масс не учитывался, ввиду его малости.

Сила тяжести Р направлена вертикально вниз. Сила лобового сопротивления Я. направлена противоположно вектору скорости V:

к = - V (3)

Сложнее обстоит дело с определением направления подъемной силы О, обусловленной эффектом Магнуса. Теоретического обоснования этого эффекта нет, однако основной механизм возникновения поперечной силы состоит в том, что при полете шарика встречный поток воздуха взаимодействует с потоком воздуха, приводимого в движение вращающимся шариком. В зоне смешанных потоков образуется область повышенного давления. Так, при полете шарика после подрезки область смешения находится вблизи нижней части шарика и поэтому возникает подъемная сила. Термин «подъемная сила» используется и при других направлениях поперечной аэродинамической силы. В эксперименте подъемная сила измеряется как сила направленная перпендикулярно векторам ю и V. Из опыта также известно, что она направлена в ту сторону, откуда кратчайший поворот от вектора ю к вектору V виден происходящим против часовой стрелки. Это позволяет применить правило векторного произведения двух векторов:

О = Q ^ • (4)

|Ш X V)

где прямые скобки означают модуль векторного произведения. Для случая подрезки вектор ю направлен параллельно оси 2 (рис.1). Для этого частного случая из (4) следует:

а•

Qy =^у • (5)

Q, = о.

Для наката знаки сил противоположны. При произвольном вращении могут быть отличны от нуля все три проекции подъемной силы.

2. Аэродинамические коэффициенты лобового сопротивления и подъемной силы

Определение аэродинамических коэффициентов - центральная часть проблемы исследования движения тел в сплошных средах. Для шара при отсутствии вращения коэффициент лобового сопротивления экспериментально найден в большом диапазоне

изменения скоростей [1], а при наличии вращения [2]-в более узкой области параметров.

2.1. Отсутствие вращения

Полагаем, что после приема мяча он летит поступательно или медленно вращается, когда эффект Магнуса невелик и можно принять, что коэффициент CL=0. Для поступательного движения шара известна зависимость аэродинамического сопротивления CD от числа Рейнольдса Re [1] (Re=Vd/v, d - диаметр шарика, V -кинематическая вязкость воздуха).

Эта зависимость приведена на рис. 2. На графике видно, что в области чисел Re от 1000 до 250 000 коэффициент лобового сопротивления CD практически постоянная величина, равная приблизительно 0.45. Для теннисного шарика, летящего в воздухе = 0.038 м, V = 0.15-10"4 м2/с) приведенные значения чисел Рейнольдса соответствуют скоростям от 0.4 м/с до 100 м/с, что с большим запасом охватывает диапазон скоростей шарика при игре в настольный теннис. Таким образом, при отсутствии вращения шарика можно принять для расчетов, что CD = 0.45, CL = 0.

2.2. Наличие вращения

Рассмотрим сложное движение шарика, состоящее из поступательного движения со скоростью V и вращательного движения вокруг оси, проходящей через центр масс с угловой скоростью ю. Важной характеристикой движения является отношение экваториальной скорости, обусловленной вращением, к скорости поступательного движения

ш d

n =

2V

log Re

Рис. 2. Зависимость аэродинамического коэффициента лобового сопротивления от числа Рейнольдса.

Рис. 3. Зависимость аэродинамических коэффициентов лобового сопротивления Св и подъемной силы Сь от отношения экваториальной скорости к скорости потока:

1 - Яв = 6.15-104, 2 - Яв = 7.74-104, 3 - Яв = 10.70-104

В работе [2] приведены зависимости значений коэффициентов Св и Сь от этого отношения для некоторых значений числа Рейнольдса. Эти зависимости показаны на рис. 3. На гранях видно, что при достаточно быстром вращении мяча аэродинамические коэффициенты слабо зависят от отношения п. Как показано в следующем параграфе, при ударе по неподвижному мячу достигается п = 1.5. При ударе по вращающемуся мячу это отношение может быть и больше. В соответствии с рис.3 примем, что Св=0.6, Сь=0.35.

Значения аэродинамических коэффициентов Св и Сь получены для больших скоростей мяча - от 25 до 40 м/с. Однако, поскольку они слабо зависят от Яв, то это дает основание принять такие же значения коэффициентов для меньших скоростей. Такое соотношение оправдывается также согласованием расчетных траекторий с реальными.

3. Оценка угловой скорости вращения теннисного шарика

Рассмотрим удар ракеткой по неподвижному шарику. Полагаем, что нет

проскальзывания шарика по ракетке, когда достигается наибольший эффект вращения

шарика.

Скорость ракетки и и составляет угол а с плоскостью ракетки. После удара центр масс ракетки получает скорость V, составляющие которой Ут и Уп, и шар вращается с угловой скоростью ю. Составляющая скорости ¥п находится через скорость и:

Уп = ки Бт а, (6)

где к - коэффициент восстановления при ударе (для мягких покрытий ракеток к « 0.5)

Для определения Ут и ю воспользуемся теоремами об изменении количества движения и кинетического момента шара при ударе. Так как до удара шар был неподвижен, то из теорем имеем:

шУх= 8Х, (7)

Зс та = £т г, Зс = -2 тг 2, (8)

где Зс - момент инерции шарика, полого внутри, относительно оси, проходящей через

центр масс; г = 2 - радиус шарика; £т - касательная компонента ударного импульса,

действующего на шарик со стороны ракетки. Система уравнений (7)-(8) дает решение ш г

—— = 1,5. Наиболее эффективное вращение происходит при малых углах а, когда ¥п

мало. Тогда допуская, что У = Ут окончательно получим отношение скоростей

шг

п= — =15-

используемое в предыдущем параграфе.

4. Уравнения движения центра масс шарика

т а с = Р + Я + О, (9)

где ас - ускорение центра масс шарика (ас = ), Р, Я, О - силы, описанные в

СII

параграфе 1. Проектируя левую и правую части (10) на оси координат и вводя обозначения Ух = и, Уу = V, получим систему уравнений:

т“ж _- к¥ * вг ■

&х___и

& , т| = ± ди,

І _ V •(>«)

R = 1 pCdsv 2, Q = 2 pClsv 2,

V = т/ и 2 + V2 ,

где х, у - координаты центра масс шарика, верхний знак в слагаемом с подъемной силой относится к подрезке, нижний - накату. Задаваемые параметры задачи:

3

р = 1.23кг/м , т = 0.0023кг , £ = ^ = 0.038м). В случае прямого удара

С в = 0.45, Сь = 0, при накате и подрезке Св = 0.6, Сь = 0.35.

Решается задача Коши с начальными условиями для функций х(1), у(1), и@), у(^). Решение находилось численно методом Эйлера.

5. Оптимизация приема мяча

Для отладки алгоритма первоначально было решено несколько тестовых задач -падение теннисного шарика без начальной скорости, движение шарика, брошенного вертикально вверх и под углом к горизонту. Полученные результаты совпадают с данными, приведенными в [3]. Однако в работе [3] не рассматривалось влияние

Рис. 4. Оптимальная траектория после приема мяча в ближней зоне прямым ударом. Параметры: а =10.5?, У0 = 11.2 м/с, V = 6.46 м/с, расстояние до стола - 0.3 м, высота - на уровне стола, удар - прямой.

У

Рис. 5. Оптимальная траектория после приема мяча в ближней зоне накатом. Параметры: а =11?, У0 = 28 м/с, Ук = 14.26 м/с, расстояние до стола - 0.3 м, высота - на

уровне стола, удар - накат.

Рис.7. Оптимальная траектория после приема мяча в ближней зоне подрезкой. Параметры: а =8?, У0 = 8.75 м/с, Ук = 5.57 м/с, расстояние до стола - 0.3 м, высота - на уровне стола, удар - подрезка.

вращения шарика на его траекторию. В данной работе исследовались траектории при подрезке и накате, и расчетные траектории сравнивались с опытными данными.

Оптимизация траекторий полета теннисного шарика проводилась в визуальном режиме. Задавалось начальное положение шарика, модуль и направление его начальной скорости. Расчетная траектория выводилась на экран компьютера, и подбирались такие начальные условия для скорости, чтобы шарик упал на сторону стола соперника с максимальной скоростью, не задев при этом сетки.

Исследовались приемы подачи прямым ударом, накатом и подрезкой в ближней зоне на высоте стола (рис. 4-6) и в дальней зоне также на высоте стола (рис 7-9). Прием высоко летящего шарика не рассматривался, так как в этом случае нет задачи оптимизации - конечная скорость шарика увеличивается неограниченно с ростом начальной скорости и сетка этому не мешает.

Выяснено, что при приеме мяча в ближней зоне наиболее эффективен накат, при котором скорость падения шарика на сторону соперника более чем в два раза превышает скорость при прямом ударе и приеме подрезкой (рис 7-9).

При приеме мяча в дальней зоне скорость падения шарика на сторону соперника существенно ниже и более эффективен прямой удар. Несколько меньшие скорости получается при приеме шарика накатом.

У

-Ч N. > N

Рис. 7. Оптимальная траектория после приема мяча в дальней зоне прямым ударом. Параметры: а =5.15?, У0 = 26.7 м/с, Ук = 9.2 м/с, расстояние до стола - 3.13 м, высота - на

уровне стола, удар - прямой.

А

У

х

Рис. 8. Оптимальная траектория после приема мяча в дальней зоне подрезкой. Параметры: а = -1?, У0 = 15.2 м/с, Ук = 5.88 м/с, расстояние до стола - 3.13 м, высота - на

уровне стола, удар - подрезка.

і \

У

N

N

Ч

\

Рис. 9. Оптимальная траектория после приема мяча в дальней зоне накатом. Параметры: а = 22.9?, У0 = 30 м/с, Ук = 8.39 м/с, расстояние до стола - 3.13 м, высота - на

уровне стола, удар - накат.

6. Заключение

Методами математического моделирования проведено исследование траекторий теннисного шарика после его приема прямым ударом, накатом и подрезкой. Учет силы аэродинамического сопротивления позволил получить результаты, согласующиеся с опытом. Выяснено, что наиболее эффективным является прием шарика накатом в ближней зоне. Если все-таки шарик принимается в дальней зоне, то расчеты отдают предпочтение прямому удару по шарику.

Литература

1. Современное состояние гидроаэродинамики вязкой жидкости. Под ред. С. Гольдштейна. Москва, Государственное издательство иностранной литературы, т. I, 1948.

2. Современное состояние гидроаэродинамики вязкой жидкости. Под ред. С. Гольдштейна. Москва, Государственное издательство иностранной литературы, т. II, 1948.

3. Р.Н. РУДАКОВ. Общие теоремы динамики и их приложение к решению задач биомеханики. Пермь, Пермский государственный технический университет, 1999.

OPTIMIZATION OF TRAJECTORY OF A TABLE TENNIS BALL

R.N. Rudakov, A.V. Kamenskhich, P.V. Shulgin (Perm, Russia)

The ball rotation plays an important role in the table tennis as it creates the lifting force (Magnus’ effect). The dependence of aerodynamical forces on translation velocity and angular velocity of the ball rotation has been studied earlier by an experiment. In this article it has been demonstrated that parameters of ball motion in the table tennis are in good agreement with the diapason of investigated parameters. It allows to determine the aerodynamical forces acting on the table tennis ball. The differential equations of motion of the center of mass of the table tennis ball were compiled which were solved numerically by the method of step integration. The trajectories of motion of the ball were investigated after its receiving in neighboring and distant zones in the case of direct central impact on the ball. The analogous solutions were obtained for upper rotation of the ball. The optimal methods of ball receiving and optimal trajectories were found which ensure ball hit on the table side of competitor with maximal velocity. Ref. 3.

Key words: table tennis, aerodynamical forces, optimal trajectories

Получено 10 май 2000