0% found this document useful (0 votes)
29 views5 pages

3.4.5 Soln

The document discusses the dynamics of a bead on a wire, including the forces acting on the bead, the equation of motion, equilibrium solutions, stability analysis, and a bifurcation diagram. It provides equations for the restoring force, component of force along the wire, Newton's second law of motion, equilibrium solutions, stability of fixed points, and a dimensionless form of the equation of motion. Graphs and analysis are presented to determine the number of equilibrium points as a parameter is varied.

Uploaded by

Vivek
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
29 views5 pages

3.4.5 Soln

The document discusses the dynamics of a bead on a wire, including the forces acting on the bead, the equation of motion, equilibrium solutions, stability analysis, and a bifurcation diagram. It provides equations for the restoring force, component of force along the wire, Newton's second law of motion, equilibrium solutions, stability of fixed points, and a dimensionless form of the equation of motion. Graphs and analysis are presented to determine the number of equilibrium points as a parameter is varied.

Uploaded by

Vivek
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 5

TFY4305 solutions exercise set 5 2014

Problem 3.5.4
a) Looking at Fig. 1 in the textbook, we see that the restoring force is given

F = ( h2 + x2 − L0 )k . (1)

The force is in the direction of the spring. We need the component of the force along the
wire, which is given by

Fwire = F sin θ
x
= F√ 2 , (2)
h + x2
Newton’s second law of motion then becomes
!
L0
mẍ + bẋ + kx 1 − √ 2 = 0. (3)
h + x2

b) Equilibrium solutions are given by ẍ = ẋ = 0. This yields


!
L0
kx 1 − √ 2 = 0. (4)
h + x2
q
Thus x = 0 and 1− √hL2 +x
0
2 = 0. This gives x = 0 and x = ± L20 − h2 . The nonzero solutions
exist only if L0 ≥ h.

c) If m = 0, the equation of motion reads


!
kx L
ẋ = √ 0 −1 . (5)
b h + x2
2

The function that determines the stability of the fixed points is


!
kx L
f (x) = √ 0 −1 . (6)
b h + x2
2

1
2

The derivative is given by


L 0 x2
!
k L
0
f (x) = √ 0 − 1 − 3 . (7)
b h2 + x2 (k 2 + x2 ) 2
This yields
k L0
 
0
f (0) = −1 (8)
b h
The origin is stable for L0 < h and unstable for L0 > h. Moreover
L2
!
0
q k
f (± L20 − h2 ) = 1 − 20 . (9)
b h
q
Thus the fixed points x = ± L20 − h2 are stable since L0 ≥ h. The stability of the fixed
points depends on the ratio y ≡ L0 /h. We therefore plot the bifurcation diagram as a func-
tion of this dimensionles variable.
x

L0 h
-1 0 2 3

-1

-2

Figure 1: Bifurcation diagram for the bead on a wire. The origin changes stability in x = 0.

d) We can rewrite the equation of motion as


m d2 x b dx
!
L0
+ + x 1 − √ = 0. (10)
k dt2 k dt h2 + x2
Defining a new dimensionless time variable τ via
b
τ. t = (11)
k
In terms of the new time variable, Newton’s equation reads
m k 2 d2 x dx
!
L0
2 2
+ +x 1− √ 2 = 0. (12)
k b dτ dτ h + x2
One can therefore ignore the second-order term if
mk
1. (13)
b2
3

Problem 3.6.5
The forces acting on the bead are gravity and the force from the spring (see Fig. 3.6.7). The
component of gravity along the wire is F = mg sin θ. The force from the spring is
√ 
Fspring = k x2 + a2 − L0 . (14)
The component of this force along the wire is
x
Fwire = − √ Fspring . (15)
x 2 + a2
The sum of F and Fwire is zero when the system is in equilibrium. This yields
!
L0
mg sin θ − kx 1 − √ 2 = 0. (16)
x + a2
b) Dividing Eq. (16) by kx, we can write
 
mg L0 1
sin θ = 1 − q  . (17)
kx a 1 + ( x )2
a

Rearranging terms, we find


mg a L0 1
1− sin θ = q . (18)
ka x a 1 + ( x )2
a

This is on the form


h R
1− = √ , (19)
u 1 + u2
mg
if we identify u = x/a, R = L0 /a, and h = ka
sin θ.

c) Note that the variable u can be positive and negative. Without loss of generality we
can restrict h to positive values (negative values correspond to tilting the wire the other
way) The function √ g(u) = 1 − h/u approaches zero as u → ±∞. The maximum of the
function h(u) = R/ 1 + u2 is R. For R < 1, it is then clear that the functions intersect at
a single point and hence there is one fixed point. This is shown in Fig. 2
For R > 1, the function g(u) still intersects the function h(u) for positive values of u.
However, there may be one or two intersections for negative values of u depending on h. In
fact, as Figs. 3 and 4 suggest, a saddle-node bifurcation is taking place at a critical value of
h. Since the derivatives of the functions g(u) and h(u) are equal at the bifurcation, h and u
satisfy the two equations
h R
1− = √ , (20)
u 1 + u2
h Ru
2
= − 3 . (21)
u (1 + u2 ) 2
4

2.0

1.5

1.0

0.5

u
-4 -2 0.0 2 4

1
Figure 2: Graphical analysis. R = 2
and h = 12 .

2.0

1.5

1.0

0.5
u
-4 -2 2 4

0.0

Figure 3: Graphical analysis. R = 2 and h = 0.44.

d) The right-hand-side of Eq. (19) can be expanded in powers of u. This yields


h 1
1− ≈ (r + 1)(1 − u2 ) , (22)
u 2
where r = R − 1. Ignoring the term ∼ ru2 which is much smaller than the term ∼ u2 , we
obtain
h 1
− ≈ r − u2 , (23)
u 2
or
1 3
u − ru − h = 0 . (24)
2

e) We first define the functions f (u) = 12 u3 − ru and g(u) = h(u). The bifurcation takes
place when f (u) = g(u) and f 0 (u) = g 0 (u). The latter equation gives
3 2
u −r = 0. (25)
2
This yields R(u) = 1 + 32 u2 . The equation g(u) = h(u) gives h(u) = 12 u3 − ru = −u3
5

2.0

1.5

1.0

0.5
u
-4 -2 2 4

0.0

Figure 4: Graphical analysis. R = 2 andh = 0.4.

f) The exact equations for the bifurcation are given by Eqs. (20) and (21). Solving the
latter with respect to uh and inserting it into the former, one finds

Ru2 R
1+ 3 = √ . (26)
(1 + u2 ) 2 1 + u2

Solving this equation with respect to R gives


3
R(u) = (1 + u2 ) 2 . (27)

Using Eq. (20), we find

Ru3
h(u) = − 3 = −u3 (28)
(1 + u2 ) 2

Note that the exact Eqs. (27) and (28) reduce to the approximate solutions upon expanding
R(u) to first nontrivial order in u for small u.

g) and h) were not a part of the exercise.

You might also like

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy