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

УДК 531/534:57+612.7

ОПТИМИЗАЦИЯ СЛАЛОМНОЙ ТРАЕКТОРИИ НА КРИВОЛИНЕЙНОМ СКЛОНЕ

Р.Н. Рудаков*, А.Ф. Лисовский**, А.Р. Гайсина*, В.В. Хитрюк***

*Пермский государственный технический университет, кафедра теретической механики, 614600, Россия, г. Пермь, Комсомольский проспект, д. 29а, e-mail: rrn@theormech.pstu.ac.ru **Чайковский институт физической культуры, 617762, Россия, г. Чайковский, ул. Ленина, д. 67 ***Средняя школа № 9 г. Перми, 614600, Россия, г. Пермь, Комсомольский проспект, д. 45

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

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

Введение

Решается задача математического моделирования прохождения горнолыжником слаломной трассы. Ранее в работах [1, 2] из сопоставления опытных и расчетных данных были найдены аэродинамический параметр системы лыжник-лыжи при спуске горнолыжника по снежному склону и усредненное значение коэффициента сопротивления скольжению лыж, обусловленного врезанием кантов лыж в снег при выполнении поворотов.

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

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

© Р.Н. Рудаков, А.Ф. Лисовский, А.Р. Гайсина, В.В. Хитрюк, 2003

Уравнения движения точки по криволинейной поверхности в естественных

координатах

Рассмотрим движение лыжника по цилиндрической поверхности, образующая которой горизонтальна, а касательная т1 к линии наибольшего ската I составляет угол а с горизонтальной плоскостью (рис. 1). Профиль склона в этом случае можно задать уравнением

а = а0 I рад, (1)

яс

где I - криволинейная координата, отсчитываемая от I = 10 по линии наибольшего ската, а0 - значение угла при I = 10, а Яс - переменный радиус кривизны склона.

Рассмотрим движение лыжника как движение материальной точки, которая движется по траектории

г=г() > (2)

где г - ее отклонение от осевой линии наибольшего ската I. В произвольной точке М касательная к траектории т составляет угол р с касательной т1 к линии наибольшего ската в этой же точке. Значение угла ф находится по формуле

р = аг^ г'. (3)

Здесь и далее штрих означает производную по I.

Рассматривая движение лыжника как движение материальной точки,

совпадающей с центром масс, запишем II закон И. Ньютона

та = Р, (4)

где т - масса системы лыжник-лыжи, а - ускорение ее центра масс, Р - главный вектор всех внешних сил, действующих на эту систему.

Левую и правую части уравнения (4) спроектируем на оси, связанные как с траекторией, так и с поверхностью горы [3]. Выберем в качестве естественных осей

следующие оси: т - касательную к траектории, п1 - нормаль к поверхности горы,

направленную к центру кривизны линии наибольшего ската, п2 - нормаль к траектории, лежащую в плоскости, касательной к поверхности горы в данной точке траектории. Главная нормаль п к траектории отходит от поверхности из-за ее изгиба и составляет угол у с осью п2 (рис. 1).

Введем орты касательной ф0 и главной нормали n0 к траектории. Тогда ускорение точки может быть записано в виде

а = атф 0 + ап n o, (5)

а уравнение (4) с учетом формул касательного и нормального ускорений примет вид:

(dV V2 'I

m-Н>0 + — П0 = F, (6)

^dt р )

где V - скорость точки, р - радиус кривизны траектории. Проектируя левую и правую части (6) на оси естественного трехгранника Мтп1п2, получим уравнения:

dV V2 . V2 ds

m----= FT, m—sin/ = Fn , m — cos/ = Fn , V =—. (7)

dt р 1 р 2 dt

Здесь s - криволинейная координата, отсчитываемая от начального значения s0

по траектории.

Проведем через касательную т плоскость, перпендикулярную к поверхности в точке M. Получим нормальное сечение поверхности, радиус кривизны которого в точке М обозначим через RN. Кроме того, рассечем поверхность соприкасающейся

плоскостью Мтп и получим наклонное сечение, которое составляет угол 90° - / с нормальным сечением. По теореме Менье радиус кривизны наклонного сечения, совпадающий с радиусом кривизны траектории р , находится по формуле:

р = rn sin /. (8)

Радиус кривизны RN может быть найден через главные кривизны K1 и K2 в данной точке поверхности по формуле Эйлера:

= K1 cos2 р + K2sin2 р, (9)

rn

где р - угол между произвольно ориентированным нормальным сечением и сечением с кривизной K1 .

Подставляя формулу (8) в уравнения (7), получим естественные уравнения движения по траектории, лежащей на любой криволинейной поверхности:

dV г, V2 V2 ds

m = FT, m-------= Fn , m----------= Fn , V =—. (10)

dt Rn ni Rn tg/ n2 dt K '

Для цилиндрического склона горы эти уравнения существенно упрощаются.

Главные кривизны склона K1 = 1/ R'c, K2 = 0 . Тогда из формулы (9) получим

Rn =-R2-. (11)

cos р

Из геометрических построений легко находится угол /. Он удовлетворяет уравнению

rc

rn

где гс - радиус кривизны траектории (2), лежащей на плоском склоне, который может быть определен по формуле:

Гс (13)

Подставим (11) и (12) в уравнения (10) и получим естественные уравнения движения точки по цилиндрической поверхности:

dV р V2 2 г? mV2 ds

m-= ft, m cos p= Fn , ---------------= Fn , V = —. (14)

dt T Rc * n1 rc n2 dt V '

В частном случае плоского склона достаточно в уравнениях (14) положить

Rc = го, тогда ось n2 совпадет с главной нормалью к траектории n .

Силы, действующие на лыжника при движении по слаломной траектории

На рис. 1 изображены внешние силы, действующие на лыжника при его движении по слаломной траектории. На него действует сила тяжести Р, сила лобового сопротивления Я, подъемная сила Q, сила трения .тр и реакция снежного покрова

склона с составляющими Ы1 и Ы2 :

Р = т%, Я = 1- рУСсГ2, Q = | рС^'2, .^р = /д/лт.2 + N2,

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

Р = Р + Я + О + ¥тР + N1 + N2. (15)

Его проекции на естественные оси координат в соответствии с рис. 1 запишутся

в виде

. = Р БтасоБ^-Я - .тр, (16)

( а' ^

=(- mg cosa + Q + N1)

v a 1 у

(17)

.п =-mg втавт^т^г + Ы2. (18)

2 1л

В формуле (17) в рассматриваемом случае вогнутой горы а' < 0. В формуле (18) знак первого слагаемого в правой части зависит от угла между осями г1 и п2. Если угол острый, то знак плюс, в противном случае - минус.

Если подставить формулы (16)-(18) в уравнения (14), то получатся естественные уравнения движения лыжника по цилиндрическому склону горы с учетом действующих на него сил. Эти уравнения приводятся к двум дифференциальным уравнениям первого порядка следующего вида:

^ ^ ,, /Т72 772 ds

m— = mgsinacosp-R-/Jn12 + N2, — = V. (19)

dt dt

Здесь

тГ2 соб2 р|а"

Ы1 =-------------!—[ + mgсова - Q > 0,

Яса

тГ2 у"

Ы2 =----------------------------------------ъ mg б1п аъ\пр^-., (20)

гс У1

1 2 1 2

Я = 2р£Сиу2, Q = -ръсг2.

Искомые функции удовлетворяют начальным условиям:

^ = 0: ^ = 50, V = Г0. (21)

Задача Коши (19)—(21) решается численно методами пошагового

интегрирования. На каждом шаге по времени, кроме V и s, определяются криволинейная координата I, углы а и р .

Определение коэффициента сопротивления скольжению лыж в повороте

При врезании кантов лыж в снег на поворотах приходит значительное торможение лыжника, существенно возрастает сила трения лыж о снег. В работе [2] выполнена оценка коэффициента сопротивления / по сравнению времени спуска горнолыжника по модельной трассе и времени, найденного из решения задачи (19)-(21). Было найдено среднее значение коэффициента, хотя очевидно, что его величина зависит от глубины врезания кантов лыж в снег. В данной работе мы предлагаем новую методику оценки коэффициента сопротивления скольжению лыж в поворотах.

Во время практических экспериментов, проведенных на горнолыжной базе «Снежинка» (г. Чайковский Пермской области), был промерен профиль склона на участке горнолыжной трассы. Обработка результатов измерений показала, что вогнутый профиль склона близок к поверхности цилиндра и имеет радиус Яс , равный

Я

520 метрам [1]. Начальный угол склона а0 составил 32,8 ------------- радиан. При

180

произвольном значении криволинейной координаты I угол склона определяется по формуле (1).

На промеренном участке были проведены два эксперимента. В первом эксперименте измерялось время прямого спуска горнолыжника в средней стойке без поворотов на трассе длиной 180 м. Для системы лыжник - лыжи общей массой 78,08 кг среднее время спуска составило 10,32±0,03 с. Из сопоставления этого времени с расчетным временем на трассе был найден аэродинамический параметр системы ц, который характеризует воздействие всех аэродинамических сил:

Ц = £ с - с ). (22)

Коэффициент ц оказался равным 0,46 [2].

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

у = а I), (23)

где: амплитуда а = 4 м, а период синуса Ь = 20 м. Общая длина спуска по линии наибольшего ската составляет 180 м. Среднее время спусков для того же лыжника, найденное в опыте, оказалось равным 39,60±0,04 с. [2].

Для определения расчетного времени спуска было принято, что коэффициент сопротивления движению лыж в поворотах линейно зависит от тангенса угла наклона лыжника к нормали к поверхности. В одноточечной модели лыжника (материальная точка на конце невесомого стержня) получается, что тангенс этого угла равен отношению составляющих реакции лыжни Ы2/ N1 (20). Тогда

/ = 0,05 + К^. (24)

N

При этом вблизи точек перегиба траектории горнолыжника при переходах от одного поворота вокруг шестов разметки трассы к другому составляющая Ы2 проходит через нулевое значение и коэффициент сопротивления / равен коэффициенту трения 0,05, как это принималось при прямом спуске без поворотов.

При решении задачи Коши (19)-(21) было принято, что коэффициент сь = 0, а ранее найденный аэродинамический параметр ц = 8св составил величину 0,46. Численное решение задачи Коши дало время спуска I = 39,6 с при К = 0,25. На каждом периоде синусоиды коэффициент / изменяется от 0,05 до максимального значения. Как и следовало ожидать, наибольшее значение коэффициент сопротивления движению принимает при обходе горнолыжником шеста разметки трассы, в этом случае угол наклона лыжника к нормали по отношению к поверхности доходит до 66 градусов, а значение /тах до 0,6.

Модель траектории (23) при заданных а и Ь не имеет варьируемых параметров. Поэтому была предпринята попытка построения периодической траектории движения горнолыжника из полиномов достаточно высокого порядка таким образом, чтобы некоторые параметры траектории можно было изменять.

Варьируемая модель слаломной траектории

Сначала для заданных амплитуды а и периода Ь построена одна четверть периода Ь траектории у1 = у1(/1). Считалось, что на осевой линии (11 = 0 ) задается угол входа в поворот рс и эта точка кривой является точкой перегиба, а при обходе флажка (/1 = Ь /4 ) радиус кривизны траектории составляет величину Ь:

/1 = 0: у1 = а у" = Х^рс, У " = 0

/ Ь , 0 . 1 (25)

11 =-: У1 =^ У" =0, У =-т.

4 Ь

Построим траекторию в виде полинома

у1 = с1 + с2/1 + с3/1 + с

Выполнение условий (25) определяет коэффициенты с

У1 = с1 + С211 + С3/12 + С413 + С51 + С0/15. (26)

60

А

р р, град

40

,.и с

20

- 10 х Ь, м

0

Рис. 2. Зависимость оптимального угла входа в поворот р Р (сплошная линия), радиуса кривизны траектории Ь при обходе шеста (штриховая линия) и времени спуска t (штрих-пунктир) от расстояния ворот а до осевой линии при периоде траектории 1=20 м

Полученная часть траектории достраивается до полного периода. Пусть переменная / изменяется от 0 до Ь. При построении траектории у(/) учтена следующая симметрия кривой:

Строится требуемое число периодов кривой и получается зависимость у (/) на всей модельной трассе.

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

уравнений (19) при начальных условиях (21), соответствующих условиям проведения эксперимента.

Расчеты показали, что существует оптимальный угол входа в поворот рР = 52 ° при Ь = 1,37 м, дающий минимальное время спуска t = 34 с. Эти результаты получены

С1 = °, с 2 = tgрс, с 3 = 0,

(27)

0 < / < Ь/4: /1 = /, у(1) = У1 (/1);

Ь/4 < / < Ь/2: /1 = Ь/2 -/, у(/) = у1(/1); Ь/2 < / < 3Ь/4: /1 = / - Ь/2, у(/) = - у х (/х); 3Ь/4 < / < Ь : /х = Ь-/, у(/) = -у1(/1).

Минимизация времени спуска по модельной слаломной трассе

при следующей расстановке ворот: a = 4 м, L = 20 м. На рис. 2 показана зависимость оптимального угла входа в поворот (pP и соответствующего радиуса кривизны траектории b при обходе шеста, а также времени спуска t от отклонения ворот a от осевой линии при периоде L = 20 м.

В проведенном исследовании получен оптимальный угол входа в поворот (pP, осредненный по всем периодам трассы. Учет различия углов (pP на разных периодах траектории существенно усложняет задачу оптимизации. Эти различия будут тем меньше, чем меньше различаются скорости входа в поворот на разных периодах траектории движения центра масс системы лыжник-лыжи.

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

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

Заключение

В результате проведенной работы рассмотрено движение горнолыжника по модельной горнолыжной трассе на криволинейном аналоге реального склона. Составлены естественные дифференциальные уравнения движения горнолыжника, учитывающие кривизну склона горы. Из сопоставления расчетного времени спуска по горнолыжной трассе со временем, полученном в эксперименте, был найден коэффициент сопротивления скольжению лыж на снегу в предположении его зависимости от угла наклона лыжника к нормали к склону горы.

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

Таким образом, в работе впервые дано теоретическое обоснование и решена задача минимизации времени прохождения горнолыжником размеченной шестами горнолыжной трассы в сделанных авторами допущениях.

Полученные результаты представляют интерес для теории и методики горнолыжного спорта, для тренеров и спортсменов-горнолыжников и могут быть использованы для постановки новых экспериментов.

Список литературы

1. Rudakov R.N., Nyashin Yu. I., Podgayets A.R., Lisovski A.F., Miheeva S.A. The influence of aerodynamic forces on the movement of sportsmen and sport balls // Russian Journal of Biomechanics. 2001. V. 5. № 2. P. 83-103.

2. Рудаков Р.Н., Гайсина А.Р., Лисовский А.Ф., Подгаец А.Р. Математическое моделирование специального слалома // Сб. мат. Всероссийской научно-практ. конф. «Перспективные технологии и методики в спорте, физической культуре и туризме». Чайковский, 2002. С. 194-198.

3. Бухгольц Н.Н. Основной курс теоретической механики. М.: Наука, 1967.

THE SLALOM PATH OPTIMIZATION ON THE CURVILINEAR SLOPE

R.N. Rudakov, A.F. Lisovski, A.R. Gaysina, V.V. Hitryuk (Perm, Russia; Tchaikovsky, Russia)

The problem of optimization of an alpine skier trajectory along the model skiing track on a curvilinear analogue of a real hill surface is under consideration in this paper. Comparison of calculated and experimental data produced hydrodynamic forces and forces of resistance to movement of skis in snow at turns. A piecewise-polynomial trajectory with two parameters was constructed. The following parameters of optimal trajectory with minimal time required for finishing the track have been found: the angle of skier’s entry into the turn and radius of curvature of trajectory while turning around the pole. The dependence of optimal trajectory parameters on the deviation of the gates from the axis line was investigated.

Key words: slalom path, optimization, aerodynamic forces, frictional forces.

Получено 1 мая 2003