2020
Математика и механика
№ 68
УДК 531.332
Б01 10.17223/19988621/68/5
С.О. Гладков, С.Б. Богданова
К ТЕОРИИ ПРОСТРАНСТВЕННОЙ БРАХИСТОХРОНЫ
Приводится общий вид функционала для задачи о брахистохроне в случае п-мерного евклидового пространства. С помощью основных принципов вариационного исчисления получена система дифференциальных уравнений, позволяющая найти формальное решение поставленной задачи. Найдено параметрическое решение системы уравнений в трехмерном и п-мерном случаях. Аналитически и численно доказано, что в пространственном случае, так же, как и в двухмерном, линией наибыстрейшего скатывания будет плоская кривая, что подтверждается графической иллюстрацией полученных решений.
Ключевые слова: трехмерная брахистохрона, п-мерный случай, функционал, экстремаль.
Во всех известных нам публикациях, касающихся анализа задач, так или иначе связанных с брахистохроной, речь всегда ведется о двухмерной задаче (см., к примеру, работы [1-7] и монографию [8]) за исключением работы [9], в которой был рассмотрен трехмерный случай вращающейся брахистохроны.
В этом коротком сообщении мы проведем подробное исследование в том случае, когда брахистохрона многомерная, но вначале подробно остановимся на трехмерном случае при условии, что координаты у, 2 являются независимыми и
могут зависеть лишь от третьей координаты х, как от параметра.
При этом совершенно понятно, что геометрия желоба, по которому скатывается тело массой т, имеет абсолютно произвольную ориентацию в пространстве.
В рамках поставленной задачи общий функционал в результате будет иметь вполне очевидный вид
где 5/ - интервал времени в условиях скатывании тела, g - ускорение силы тяжести, И - максимальная высота тела над Землей.
С целью решения поставленной задачи удобно воспользоваться уравнениями Эйлера - Лагранжа [10], которые имеют вид
Решение и анализ задачи в 3-мерном случае
(1)
Vi+y^+Z'2
где функция F = ——. -. Заметим, что в силу линейности уравнений (2)
vh - z
множитель 1_ не существен, и мы его учитывать не будем. V2g
Из верхнего уравнения системы (2) следует, что
Fy' = C1 = const, (3)
а нижнее, в силу независимости подынтегральной функции от x, дает
F - ZFz, = C2, (4)
где C2 - константа интегрирования.
После простых вычислений приходим к следующей системе уравнений:
[ y' = 1 + У'2 + z'2,
[i + y2 = C2yfh—zyl 1 + y'2 + z'2.
(5)
C
Поделив их друг на друга и вводя новую константу A = —1, приходим к соот-
C2
ношению
= a , (6)
У C2
где A = const .
Отсюда сразу же следует решение
y = ax + b , (7)
A IA2
где а = — - у —— 1 , а Ь - константа интегрирования.
Заменяя в верхнем уравнении системы (5) производную у' на а, получаем
а = С^Н - 2^1 + а2 + г'2 . (8)
Уравнение (8) тривиально решается с помощью подстановки
г ' = 41 + а2 гяР, (9)
где угловая переменная р представляет собой параметр.
В результате простого интегрирования немедленно находим решение уравнения (8) в параметрическом виде
1 + а2
2 = Н--2—СОБ2 Р . (10)
С1
Поскольку согласно (9), dz = 41 + a2 tg р dx,
то отсюда следует, что
x = C3 +
-=!= J ctg P ded P,
где С3 - константа интегрирования. В результате подстановки сюда (10) приходим к следующему решению:
х = С3 + + ^Ш2Р|.
(11)
Таким образом, собирая решения (7), (10) и (11) в единую систему, получаем решение задачи о пространственной брахистохроне в виде
^ у/1 + а2 (а 1 . „„ х = С3 +--—I В+-Б1п2р
3 С12 I 2
у = ах + Ь,
г = И -
1 + а2
2С2
(1 + соб2Р).
(12)
Выбирая константу С12 в виде
С 2 = 1 + а1
С1 = 7
(13)
вместо (12) находим
х = С3 +
VI
И ^Р+ ■1Б1п2Р
+ а
у = ах + Ь, И
г = - (1 - соб2Р).
(14)
Считая, наконец, что брахистохрона «привязана» к началу координат, получим Ь = С3 = 0, и решение становится вполне компактным:
1
И Гр + ^1п2р
У =
+ а аИ
1
+ а
2 ^р+ -281п2р
(15)
2 = - (1 - СОБ 2Р).
Полагая в предельном случае а = 0, приходим к параметрическим уравнениям обычной двухмерной брахистохроны
х = И |р + 2б1п2Р И
2 = -(1 - СОБ2Р).
Полагая в уравнениях (15), например, а = л/3 , получаем
h f 1 x = 2 lp+2sin2P
W3 f i . y=—lP+2sm2P
h
z = - (1 - cos2p).
(17)
Из решений (17) видно, что на плоскостях х - 2 и у - 2 будут наблюдаться типичные двухмерные брахистохроны, только на плоскости у - 2 амплитуда будет в л/3 раз больше (см. рис. 1).
0.2-
0.4
0.6
0.8:
0
1
2
3
Рис. 1. Решения системы (17): линия 1 построена в предположении, что h = 1, a = л/3 , линия 2 построена при h = 1, a = -Л5 . Параметр t принадлежит сегменту [0,3п] Fig. 1. Solutions to system (17): line 1 corresponds to a case of h = 1, a = V3 and line 2, to a case of h = 1, a = л/Г5. The parameter t belongs to the interval of [0,3n]
2
4
1
Рис. 2. Проекции на плоскость xoy трех
пространственных брахистохрон: линия 1 -проекция брахистохроны с параметрами h = 1, a = л/3 ; линия 2 - h = 1, a = л/15 ; линия 3 - h = 1, a = ^18Ô для t е [0; 3п]
Fig. 2. Projections of three nonpolar brachistochrones on the plane xoy : line 1
corresponds to a brachistochrone at
h = 1, a = л/3, line 2, at h = 1, a = л/15, and
line 3, at h = 1, a = V8Ô for t е [0;3n]
и-Мерный случай
Рассмотрим теперь случай «-мерного евклидового пространства. В этом случае функционал (1) должен быть заменен на
8' w=v2F pjh-
2 + x3 + ...xn
dx1
(18)
где вектор x = ( x2, x3,...xn ).
0
0
4
z
6
8
1
3
0
X
X
В соответствии с уравнениями Эйлера - Лагранжа получаем в результате следующую систему уравнений:
: с14н-7п41+х22 + х32 + ...х'2,
х2 = <
х3 = с2^1 И - х„41 + х22 + х32 + ...х;2,
1 + х22 + х32 + ... + х;-1 = сп-1л/И - х„41 + х22 + х32 + ...х'п2. Деля почленно верхние п - 2 уравнения друг на друга, находим
= а,
(19)
х2 с:
1 с2
х2 с!
х4 сз
х2 с1
х5 с4
(20)
с,
г' С
т-2 т-3
т-4*
где т > 4, а коэффициент а0 = 1.
Таким образом, имеется следующее рекуррентное соотношение:
х2
т-2 ~ '
ат-4
где число т подчиняется неравенству
4 < т < п + 2 .
(21)
(22)
Подставляя теперь (21) в нижнее уравнение системы (19), немедленно получаем такое уравнение:
х0 х0
2 2 а а.
Или
где
1+х22 = Сп1+х22 +-у +%.. +4-+х
2 ап-1 V а3 а4 ап-1
1 + ёп-1 х22 = Сп-ЛИ - х„41 + ёп- х22 + х'п2 ,
ёп-1 = 1 + ^ + ^ > 1.
(23)
п-1
Учитывая верхнее уравнение системы (19), приходим к замкнутой системе из двух уравнений, а именно
С1Л/И - хп^1 + ёп-1
х2 + хм ,
I1 + ёп-1 х22 = сп-г4И - х^1 + ёп-1 х22 + х'п2.
Решение этой системы точно такое же, как и системы (5). Поэтому, деля почленно уравнения друг на друга, будем иметь
x2 = C1 = 1 1 + gn-1 x22 Cn-1 pn-1
Откуда x2 = An-1 x1 + B, (25)
где B - константа интегрирования, а новый коэффициент
An1 = pn-1 - Vp2-i ~4gn-L (26)
2gn-1
Поэтому из нижнего уравнения системы (24) получим
1 + gn-1 A„2-1 = Cl4h-7n41 + gn-1 A„2-! + x" . (27)
Совершая здесь подстановку
x'„=41 + gn-1 A- tg0 , (28)
находим
(1 + gn-1 A„2-!) cos2 0 xn = h--. (29)
C2 1
Поэтому согласно (28), получаем
V1 + gn-1 A2-:
1
1 С12 ^ 2
где С - константа интегрирования.
И, значит, согласно (25) из (21) следует, что
+ -sin20 1 + C , (30)
An-1 -xl B +—,
a1 a1
A- -xl B + —
a2 a2
An-1 -xl B +-
a3 a3
=A-1
xn-1 = x1 +
(31)
ап-3 ап-3
Решения (29) - (31) отвечают на вопрос о форме брахистохроны в многомерном случае.
Заключение
В заключение отметим:
1. Дано решение задачи о брахистохроне в п-мерном случае и, как пример, приводится подробное решение для трехмерного случая.
2. Отмечено принципиальное отличие многомерной брахистохроны от классической двухмерной.
ЛИТЕРАТУРА
1. Гладков С.О., Богданова С.Б. Геометрический фазовый переход в задаче о брахистохроне // Ученые записки физического факультета МГУ. 2016. № 1. С. 161101-1-5.
2. Гладков С. О. О траектории движения тела, входящего в жидкость под произвольным углом // Ученые записки физического факультета МГУ. 2016. № 4. С. 164002-1-5.
3. Гладков С.О., Богданова С.Б. К теории движения тел с переменной массой // Вестник Томского государственного университета. Математика и механика. 2020. № 65. С. 83-91. DOI: 10.17223/19988621/65/6.
4. Denman Harry H. Remarks on brachistochrone - tautochrone problems // Amer. J. Phys. 1985. V. 53. P. 224-227.
5. Scarpello G.M. and Ritelli D. Planar brachistochrone of a particle attracted in vacuo by an infinite rod // New Zealand J. Math. 2007. V. 36. P. 241-252.
6. Goldstein H. and Bender C. Relativistic brachistochrone // J. Math. Phys. 1986. V. 27. N 2. P. 507-511.
7. Scarpello G.M. and Ritelli D. Relativistic brachistochrones under electric or gravitational uniform fields // Z. Angew. Math. Mech. 2006. V. 86. No. 9. P. 736-743.
8. Эльсгольц Л.Э. Дифференциальные уравнения и вариационное исчисление. М.: Наука, 1969. 424 с.
9. Гладков С.О., Богданова С.Б. К теории движения шарика по вращающейся брахистохроне с учетом сил трения // Ученые записки физического факультета МГУ. 2017. С. 172101-1-6.
10. Ландау Л.Д., Лифшиц Е.М. Механика . М.: Физматлит, 2004. 210 с.
Статья поступила 14.12.2019
Gladkov S.O., Bogdanova S.B. (2020) TO THE THEORY OF n-DIMENSIONAL BRACHISTOCHRONE. Vestnik Tomskogo gosudarstvennogo universiteta. Matematika i mekhanika [Tomsk State University Journal of Mathematics and Mechanics]. 68. pp. 53-60
DOI 10.17223/19988621/68/5
Keywords: three-dimensional brachistochrone, n-dimensional case, functional, extremal.
In this paper, a solution to the problem of the motion of a brachistochrone in the n-dimensional Euclidean space is firstly presented. The very first formulation of the problem in a two-dimensional case was proposed by J. Bernoulli in 1696. It represented an analytical description of the trajectory for the fastest rolling down under gravitational force only. Thereafter, a number of problems devoted to a brachistochrone were considered with account for gravitational forces, dry and viscous drag forces, and a possible variation in the mass of a moving body.
Analytical solution to the formulated problem is presented in details by an example of the body moving along a brachistochrone in three-dimensional Cartesian coordinates. The obtained parametric solution is confirmed by a graphical interpretation of the calculated result. The formulated problem is solved for an ideal case when drag forces are neglected. If dry and viscous friction forces are taken into account, the plane shape of the brachistochrone remains the same, while the analysis of the solution becomes more complicated. When, for example, a side air flow is taken into account, the plane curve is replaced by a three-dimensional brachistochrone.
Sergey O. GLADKOV (Doctor of Physics and Mathematics, Professor, Moscow Aviation Institute, Moscow, Russian Federation). E-mail: sglad51@mail.ru
Sofiya B. BOGDANOVA (Candidate of Physics and Mathematics, Associate Professor, Moscow Aviation Institute, Moscow, Russian Federation). E-mail: sonjaf@list.ru
REFERENCES
1. Gladkov S.O., Bogdanova S.B. (2016) Geometricheskiy fazovyy perekhod v zadache o brakhistokhrone [Geometrical phase transition in a problem on brachistochrone]. Uchenye zapiski fizicheskogo fakul'teta Moskovskogo universiteta - Memoirs of the Faculty of Physics, Lomonosov Moscow State University. 1. pp. 161101-1-5.
2. Gladkov S.O. (2016) O traektorii dvizheniya tela, vkhodyashchego v zhidkost' pod proizvol'nym uglom [On the trajectory of the body entering the fluid at an arbitrary angle]. Uchenye zapiski fizicheskogo fakul'teta Moskovskogo universiteta - Memoirs of the Faculty of Physics, Lomonosov Moscow State University. 4. pp. 164002-1-5.
3. Gladkov S.O., Bogdanova S.B. (2020) K teorii dvizheniya tel s peremennoy massoy [To the theory of motion of bodies with variable mass]. Vestnik Tomskogo gosudarstvennogo universiteta. Matematika i mekhanika - Tomsk State University Journal of Mathematics and Mechanics. 65. pp. 83-91. DOI: 10.17223/19988621/65/6.
4. Denman H.H. (1985) Remarks on brachistochrone - tautochrone problems. American Journal of Physics. 53. pp. 224-227.
5. Scarpello G.M., Ritelli D. (2007) Planar brachistochrone of a particle attracted in vacuo by an infinite rod. New Zealand Journal of Mathematics. 36. pp. 241-252.
6. Goldstein H., Bender C. (1986) Relativistic brachistochrone. Journal of Mathematical Physics. 27(2). pp. 507-511.
7. Scarpello G.M., Ritelli D. (2006) Relativistic brachistochrones under electric or gravitational uniform fields. Journal of Applied Mathematics and Mechanics. 86(9). pp. 736-743.
8. Elsholz L.E. (1969) Differentsial'nye uravneniya i variatsionnoe ischislenie [Differential equations and calculus of variations]. Moscow: Nauka.
9. Gladkov S.O., Bogdanova S.B. (2017) K teorii dvizheniya sharika po vrashchayushcheysya brakhistokhrone s uchetom sil treniya [To the question of analysis equations for moving due to the rotating of the plane brahistohron]. Uchenye zapiski fizicheskogo fakul'teta Moskovskogo universiteta - Memoirs of the Faculty of Physics, Lomonosov Moscow State University. 2. pp. 172101-1-6.
10. Landau L., Lifshitz E. (1976)Mechanics. Oxford: Butterworth-Heinemann.
Received: December 14, 2019