Jump to ContentJump to Main Navigation
Show Summary Details
More options …

Open Physics

formerly Central European Journal of Physics

Editor-in-Chief: Seidel, Sally

Managing Editor: Lesna-Szreter, Paulina


IMPACT FACTOR 2018: 1.005

CiteScore 2018: 1.01

SCImago Journal Rank (SJR) 2018: 0.237
Source Normalized Impact per Paper (SNIP) 2018: 0.541

ICV 2017: 162.45

Open Access
Online
ISSN
2391-5471
See all formats and pricing
More options …
Volume 15, Issue 1

Issues

Volume 13 (2015)

Chain on a cone

Oleg Zubelevich
  • Corresponding author
  • Dept. of Theoretical mechanics, Mechanics and Mathematics Faculty, M. V. Lomonosov Moscow State University, 119899, Moscow, Vorob’evy gory, MGU Russia
  • Email
  • Other articles by this author:
  • De Gruyter OnlineGoogle Scholar
Published Online: 2017-11-10 | DOI: https://doi.org/10.1515/phys-2017-0078

Abstract

We consider a loop of a chain thrown like a lasso on a fixed right circular cone. The system is in the standard homogeneous gravity field. The axis of the cone is vertical. It is shown that under certain vertex angles chain’s loop has an oblique equilibrium.

Keywords: mechanics of chains

PACS: 02.40.-k

1 Statement of the Problem and the main Theorem

In inheritance from the classical epoch we have obtained the two most famous variational problems.

The first one is the catenary problem, formulated and solved by Johann Bernoulli in 1690.

The brachistochrone curve problem was independently considered and solved also approximately about 1690 by Johann Bernoulli, Christiaan Huygens and Gottfried Wilhelm Leibniz.

These both variational problems are simultaneously beautiful non trivial and can be solved by pure analytical means.

This is absolutely exclusive situation: there are very few variational problems that are not artificially composed but shows up from mechanics and can be solved by hands.

In this short note we propose one of such problems.

Ends of thin homogeneous chain are connected to obtain a loop. The mass of the chain is m and its length is l. This loop is putted on a right cone with vertex angle 2 α, α ∈ (0,π/2); see Figure 1. The surface of the cone is smooth. The axis of the cone is vertical. The system is in the standard mg–gravity field.

Chain on a cone
Figure 1

Chain on a cone

Find all equilibriums of chain’s loop on the cone.

This problem has trivial solution: the loop forms a circle lying in horizontal plane. But do another equilibriums exist?

Theorem 1

If π/6 < α < π/4, then there is a unique (up to rotations about cone’s axis) oblique equilibrium (see Figure 1 overflow).

For all another α the chain has only trivial equilibrium.

Remark 1

Numerical simulations show that at this oblique equilibrium the chain is not contained in a plane.

2 Proof

Introduce a Cartesian frame Oxyz such that the axis Z is directed down along the cone axis and the origin O is cone’s vertex. Introduce also a cylindrical coordinate frame such that (z,r,ψ),x=rcosψ,y=rsinψ.

The cone is given by the equation z=ar,a=cotα>0.

Let r=r(ψ),r(ψ+2π)=r(ψ)

be the equation of the curve described by the chain. Infinitesimal arclength element of such a curve is as follows ds=dx2+dy2+dz2=(1+a2)(r)2+r2dψ,r=ddψr(ψ).

Then Z–coordinate of the centre of mass of the chain is expressed by the formula. Z[r()]=al02πr(ψ)ds.

So we are looking for conditional extremals of the functional r(⋅)↦ Z[r(⋅)] in class of 2π–periodic functions r(ψ) under the condition 02πds=l.(2.1)

The corresponding Lagrange function is L(r,r)=(alr+λ)(1+a2)(r)2+r2, with a Lagrange multiplier λ [1].

Since the Lagrangian L does not depend on ψ we have the following integral of the Lagrange equations h=rLrL,

or (alr+λ)r2=h(1+a2)(r)2+r2.(2.2)

So that the constant h and the expression (alr+λ) must have the same sign for all ψ ∈ ℝ i.e. sgnh=sgn(alr+λ).(2.3)

Under this assumption (we justify this assumption in the sequel) take square from both sides of (2.2) (1+a2)(r)2+r21h2(alr+λ)2r4=0,

and perform a change of variables rρ by the formula r=λlaρ,ρ=ρ(ψ)(2.4)

to obtain (1+a2)(ρ)2+ρ2u2(ρ+1)2ρ4=0,(2.5)

where u2=λ4l2h2a2.(2.6)

The condition (2.1) takes the form 02π(1+a2)(ρ)2+ρ2dψ=a|λ|.(2.7)

Our plan is as follows. Choosing a constant u2 we will find 2π–solution ρ(ψ) of equation (2.5) and then take |λ| to satisfy (2.7). By reason that will be clear in the sequel we assume that λ<0.(2.8)

Then already known u2 and |λ| we will substitute to (2.6) and find a constant h2.

The solution ρ must also be such that the expression (alr+λ)=λ(ρ+1)

is sign-definite for all ψ. If the solution ρ(ψ) provides this condition then (2.3) is satisfied by choosing the sign of h.

Introducing new variable t=ψ/a2+1, rewrite equation (2.5) as follows 1(ρ+1)2ρ4ρ˙2+1ρ2(ρ+1)2=u2,ρ˙=dρdt.(2.9)

Now we are looking for the solution ρ(t) of equation (2.9) that has period 2πa2+1.

Equation (2.9) has the form of energy integral of a classical mechanical system with kinetic energy E=1(ρ+1)2ρ4ρ˙2

and a potential energy V(ρ)=1ρ2(ρ+1)2. Use this observation to analyze system (2.9).

Having a graph of the function V we see that the point C = (−1/2,0) on the phase plane (ρ,ρ˙) is an equilibrium of the type ”centre”. This equilibrium is surrounded with closed curves that are squeezed between vertical lines ρ = −1 and ρ = 0. These trajectories are marked with parameter u2. There are no periodic orbits in other domains of the phase space.

Therefore all the periodic solutions satisfy the condition ρ(t)(1,0),tR

consequently one has ρ(t)+1 > 0. By virtue of formulas (2.8) and (2.4) it follows that r > 0.

The equilibrium C corresponds to the value u2 = 16 and gives the trivial equilibrium of the chain.

Separating variables in (2.9), find period of the trajectory ρ(t) by the formula: T(u)=2ρρ+dρρu2(ρ+1)2ρ21,ρ±=1±14/u2.

The function T(u) is defined for u > 4. It is not hard to show that the function T is continuous and decreased.

Thus each root u of the equation T(u)=2π1+a2

corresponds to the solution of (2.9) with desired period. This solution gives the oblique equilibrium.

The proof of the theorem is concluded by a lemma.

Lemma 1

The following formulas hold

limu4+T(u)=2π,limuT(u)=π.(2.10)

Indeed, the oblique equilibrium shows up iff π<2π1+a2<2π,

or π/6<α<π/4.

The theorem is proved.

2.1 Proof of the Lemma

The first one of the formulas (2.10) is the most simple; we prove it by using a little bit of informal mechanical argument. Rigorous proof of this formula follows by the same manner as we employ below to prove the second formula of (2.10).

Up to third order terms in the neighbourhood of the point C equation (2.9) has the form ξ˙2+2ξ2=const,ρ=12+ξ.

So that the period of small oscillations is equal to 2π/2 or limu4+T(u)=2π.

Let us check the second formula. Observe that T(u)=2uρρ+dρρ(ρρ)(ρρ+)(ρρ~)(ρρ~+),

where ρ~±=1±1+4/u2.

Introducing a small parameter ε = 1/u → 0 as u → ∞, we get ρ+=ϵ+O(ϵ2),ρ=1+ϵ+O(ϵ2),ρ~+=ϵ+O(ϵ2),ρ~=1ϵ+O(ϵ2).

For small ε the following inequalities hold ρ+<ϵ,ρ>1+ϵ.

It is easy to see that T(u)=2ϵ12ρ+dρρ(ρρ)(ρρ+)(ρρ~)(ρρ~+)+O(ϵ).

Transform this integral as follows: 12ρ+dρρ(ρρ)(ρρ+)(ρρ~)(ρρ~+)=12ρ+(1+O(ϵ))dρρ(ρ+1)(ρρ+)(ρ+ρ+ρ+ρ~+).

Sinse ρ+ρ~+=O(ϵ2) the last integral equals 12ρ+(1+O(ϵ))dρρ(ρ+1)ρ2ρ+2.

The integral 12ρ+dρρ(ρ+1)ρ2ρ+2

is computed explicitly, nevertheless the formula is very large and we do not bring it; write down the asymptotic 12ρ+dρρ(ρ+1)ρ2ρ+2=π2ϵ+O(ln1ϵ).

Gathering all these formulas we yield T(u)=π+O(ϵ).

The Lemma is proved.

Acknowledgement

Partially supported by grants RFBR 12-01-00441., Science Sch.-2964.2014.1.

References

  • [1]

    van Brunt B., The Calculus of Variations, Springer, New York, 2004 Google Scholar

About the article

Received: 2017-05-08

Accepted: 2017-09-26

Published Online: 2017-11-10


Citation Information: Open Physics, Volume 15, Issue 1, Pages 671–673, ISSN (Online) 2391-5471, DOI: https://doi.org/10.1515/phys-2017-0078.

Export Citation

© 2017 Oleg Zubelevich. This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License. BY-NC-ND 4.0

Comments (0)

Please log in or register to comment.
Log in