Brachistochrone
Introduction
To find the shape of the curve which the time is shortest possible. . .
We use WxMaxima to do the calculus part.
Theory
Variational Calculus and Euler--Lagrange Equation
The time that is needed for sliding from point to point is where is the Pythagorean distance measure and is determined from the the law of conservation of energy giving . Plugging these in, we get
,
where is the function subject to variational consideration.
According to the Euler--Lagrange differential equation the stationary value is to be found, if E-L equation is satisfied.
No Friction
Since does not depend on , we may use the simplified E--L formula Constant. The differentials are easy, and we have
So we have and multiplying this with the denominator and rearring, we have by redefining the constant. The standard solution to this differential equation is given by
and is the equation of a cycloid.
No Friction: Maxima
The details using WxMaxima:
energy : 1/2*m*v^2 = m*g*y;
v_sol : solve( energy, v);
v_sol : v_sol[2];
EL_f : rhs( sqrt(1+'diff(y,t)^2)/v_sol );
doof_dooyp : diff( EL_f, 'diff(y,t));
EL: EL_f - 'diff(y,t)*doof_dooyp = C;
radcan(%);
EL_dy : solve(EL, y);
ode2(EL_dy[1]^2,y,t);
but the ode2 solver cannot handle the nonlinear differential equation.
Rolling Ball: Angular momentum but no radius
The rotational energy is and by applying non-slipping condition we get . Note that actually the ball is rolling on a curve, and thus the given slipping condition is only an approximation.
For the simplified case, the calculation is similar to the previous one, and using Maxima, we get
energy : 1/2*m*v^2 + 1/2*I*v^2/r^2= m*g*y;
. . .
gives
and thus only the constant differs from the case with no angular momentum.
Friction
The normal force follows the path, and thus is given by , but The friction depends on the normal force of the path. The normal force is perpendicular to the previous, thus we have
The conservation of energy does not apply here, but we have Newton's Second Law, . We need the components along the curve . Thus we have
Clearly, for the left hand side of NII we have , and by including the differential part only, we have
and for the Euler--Lagrange equation is
Euler--Lagrange
NII : 1/2*v^2 = g*(y(x)-mu*x);
v_sol : solve( NII, v);
v_sol : v_sol[2];
EL_f : rhs( sqrt(1+'diff(y(x),x)^2)/v_sol );
df_dy : diff(EL_f, y(x));
df_dyp : diff(EL_f, 'diff(y(x),x));
d_dx : diff( df_dyp, x);
EL : df_dy - d_dx = 0;
Elrad : radcan( EL );
num( lhs(ELrad) )/sqrt(2)/sqrt(y(x)-mu*x)=0;
ratsimp(%);
Reduction
Remember that . Then, note that . Thus, we multiply EL equation by to obtain
The left hand side can be integrated:
The right hand side can be integrated using partial fraction decompisition
Together we have
that can be written as
and it gives finally
depends(y,x ); EL: 2*( y - mu*x )*diff( y,x,2) + (1 + diff(y,x)^2)*(1+mu*diff(y,x)) = 0; factor( ratsimp(solve(EL, diff(y,x,2))*(diff(y,x)-mu)*2/(1+diff(y,x)^2)/(1+mu*diff(y,x))) ); eq1 : integrate( rhs( EL_2[1]),x) + log(C); eq2 : integrate( partfrac( lhs( EL_2[1]), diff(y,x) ), x); exp(eq1)=exp(eq2);
Solution
The solution can be obtained by setting which implies and we have .
We solve for , and get
If which means no friction, we get , which is the result obtained earlier.
.
oo: (1+mu*cot(t/2))^2/(1+cot(t/2)^2);
oo1 : oo, t/2 = s;
trigrat(oo1);
trigexpand(%);
trigsimp(%);
oo2 : expand(%);
part(oo2,1) + part(oo2,2) + trigsimp( part(oo2, 3) + part(oo2,4) );
oo3 : factor(%);
oo3, s = t/2;
Maxima:
- trigreduce : product of sinuses and cosines as a Fourier sum (with terms containing only a single sin or cos).
- trigexpand : no multiple angles. Uses sum-of-angles formulas
- trigsimp : Pythagorean identity
- substitution, eg. eq1, 2*x = y;
- trigrat : does many things?
Rolling Ball with radius
The conservation of energy:
Beltrami Indentity
E-L states: , but , and now and by substituting the first result, we have
and thus Beltrami follows.
References
https://mathworld.wolfram.com/BrachistochroneProblem.html
https://physicscourses.colorado.edu/phys3210/phys3210_sp20/lecture/lec04-lagrangian-mechanics/
http://hades.mech.northwestern.edu/images/e/e6/Legeza-MechofSolids2010.pdf
https://www.tau.ac.il/~flaxer/edu/course/computerappl/exercise/Brachistochrone%20Curve.pdf
https://mate.uprh.edu/~urmaa/reports/brach.pdf The Nonlinear Brachistochrone Problem with Friction Pablo V. Negr´on–Marrero∗ and B´arbara L. Santiago–Figueroa
https://medium.com/cantors-paradise/the-famous-problem-of-the-brachistochrone-8b955d24bdf7
https://wiki.math.ntnu.no/_media/tma4180/2015v/calcvar.pdf BASICS OF CALCULUS OF VARIATIONS MARKUS GRASMAIR
http://www.doiserbia.nb.rs/img/doi/0354-5180/2012/0354-51801204697M.pdf
http://info.ifpan.edu.pl/firststep/aw-works/fsV/parnovsky/parnovsky.pdf Some Generalisations of Brachistochrone Problem. A.S. Parnovsky
[https://arxiv.org/pdf/1604.03021.pdf Tautochrone and Brachistochrone Shape Solutions for Rocking Rigid Bodies. Patrick Glaschke April 12, 2016]
https://issuu.com/nameou/docs/math_seminar_paper A complete detailed solution to the brachistochrone problem. N. H. Nguyen.
https://arxiv.org/pdf/1908.02224.pdf Brachistochrone on a velodrome. GP Benham, C Cohen, E Brunet and C Clanet
https://arxiv.org/pdf/1712.04647.pdf On the brachistochrone of a fluid-filled cylinder. Srikanth Sarma Gurram, Sharan Raja, Pallab Sinha Mahapatra and Mahesh V. Panchagnula.
https://arxiv.org/pdf/1001.2181.pdf A Detailed Analysis of the Brachistochrone Problem R.Coleman
https://www.jstor.org/stable/2974953?seq=1#metadata_info_tab_contents Exploring the Brachistochrone Problem. LaDawn Haws and Terry Kiser
https://math.stackexchange.com/questions/3077935/solving-the-euler-lagrange-equation-for-the-brachistochrone-problem-with-frictio.