Wedge, Block, and Spring

A triangular wedge has height h h and mass M M , and it makes an angle θ \theta with the horizontal. The wedge sits on a smooth horizontal surface, and can slide freely on the surface. The right side of the wedge is connected to a horizontal spring of force constant k k . The ambient gravitational acceleration g g is downward.

There is a block of mass m m on the slanted part of the wedge, where it can slide freely. Initially, both the block and the wedge are at rest, the block is at the top of the wedge, and the spring is at its natural length.

When the block gets to the bottom of the wedge, how far is the wedge from its initial position?

Details and Assumptions:
1) h = 10 h = 10
2) M = 1 M = 1
3) θ = π / 6 \theta = \pi/6
4) m = 1 m = 1
5) g = 10 g = 10
6) k = 4 k = 4
7) Assume that the block is confined to the surface of the wedge


The answer is 0.806.

This section requires Javascript.
You are seeing this because something didn't load right. We suggest you, (a) try refreshing the page, (b) enabling javascript if it is disabled on your browser and, finally, (c) loading the non-javascript version of this page . We're sorry about the hassle.

2 solutions

Steven Chase
Jun 28, 2020

Let the horizontal coordinate of the bottom left vertex of the wedge be x 0 x_0 . Let the distance of the block up the ramp be α \alpha . The coordinates of the block are:

x = x 0 + α cos θ y = α sin θ x = x_0 + \alpha \cos \theta \\ y = \alpha \sin \theta

Let x 0 = 0 x_0 = 0 be the initial location of the wedge. The kinetic energy, potential energy, and lagrangian are:

T = 1 2 m ( x ˙ 0 2 + 2 cos θ x ˙ 0 α ˙ + α ˙ 2 ) + 1 2 M x ˙ 0 2 V = m g α sin θ + 1 2 k x 0 2 L = T V = 1 2 m ( x ˙ 0 2 + 2 cos θ x ˙ 0 α ˙ + α ˙ 2 ) + 1 2 M x ˙ 0 2 m g α sin θ 1 2 k x 0 2 T = \frac{1}{2} m (\dot{x}_0^2 + 2 \cos \theta \, \dot{x}_0 \dot{\alpha} + \dot{\alpha}^2) + \frac{1}{2} M \dot{x}_0^2 \\ V = m g \alpha \sin \theta + \frac{1}{2} k x_0^2 \\ L = T - V = \frac{1}{2} m (\dot{x}_0^2 + 2 \cos \theta \, \dot{x}_0 \dot{\alpha} + \dot{\alpha}^2) + \frac{1}{2} M \dot{x}_0^2 - m g \alpha \sin \theta - \frac{1}{2} k x_0^2

Equations of motion:

d d t L x ˙ 0 = L x 0 d d t L α ˙ = L α \frac{d}{dt} \frac{\partial{L}}{\partial{\dot{x}_0}} = \frac{\partial{L}}{\partial{x}_0} \\ \frac{d}{dt} \frac{\partial{L}}{\partial{\dot{\alpha}}} = \frac{\partial{L}}{\partial{\alpha}}

Evaluating the equations of motion results in:

( m + M ) x ¨ 0 + m cos θ α ¨ = k x 0 m cos θ x ¨ 0 + m α ¨ = m g sin θ (m+M) \ddot{x}_0 + m \cos \theta \, \ddot{\alpha} = - k x_0 \\ m \cos \theta \, \ddot{x}_0 + m \ddot{\alpha} = - m g \sin \theta

Solve for both second derivatives on each time step of a numerical simulation routine. The initial conditions are:

x 0 = 0 x ˙ 0 = 0 α = h sin θ α ˙ = 0 x_0 = 0 \\ \dot{x}_0 = 0 \\ \alpha = \frac{h}{\sin \theta} \\ \dot{\alpha } = 0

Plots of x 0 x_0 and α \alpha over time are given below. The simulation is run until α = 0 \alpha = 0 .

Simulation code:

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
import math
import numpy as np

h = 10.0
M = 1.0
theta = math.pi/6.0

m = 1.0
g = 10.0
k = 4.0

dt = 10.0**(-5.0)

#############################

x0 = 0.0

# alpha*math.sin(theta) = h
alpha = h/math.sin(theta)

U0 = m*g*h

x0d = 0.0
alphad = 0.0

x0dd = 0.0
alphadd = 0.0

t = 0.0
count = 0

#############################

print "t x0 alpha"

while alpha >= 0.0:

    x0 = x0 + x0d*dt
    alpha = alpha + alphad*dt

    x0d = x0d + x0dd*dt
    alphad = alphad + alphadd*dt

    J11 = m + M
    J12 = m*math.cos(theta)

    J21 = m*math.cos(theta)
    J22 = m

    J = np.array([[J11,J12],[J21,J22]])
    vec = np.array([-k*x0,-m*g*math.sin(theta)])

    Sol = np.linalg.solve(J, vec)

    x0dd = Sol[0]
    alphadd = Sol[1]

    t = t + dt
    count = count + 1

    if count % 100 == 0:
        print t,x0,alpha

#############################

print ""
print ""

vx = x0d + alphad*math.cos(theta)
vy = alphad*math.sin(theta)

print dt
print t
print ""
print x0
print x0d
print ""
print vx
print vy
print ""

E = 0.5*m*(vx**2.0 + vy**2.0) + 0.5*M*(x0d**2.0) + 0.5*k*(x0**2.0)
print (E/U0)
print ""
print ""

#############################

@Steven Chase Sir I am planning to post a good question, but I am not able to solve, will you help me???
Please reply.

A Former Brilliant Member - 11 months, 2 weeks ago

Log in to reply

Perhaps. What is it?

Steven Chase - 11 months, 2 weeks ago

@Steven Chase sir I have uploaded a new problem.

A Former Brilliant Member - 11 months, 2 weeks ago
Karan Chatrath
Jun 28, 2020

I considered the X-axis to be placed horizontally to the left and the Y-axis vertically upwards. The origin is at the initial position of the wedge.

The coordinates of the wedge at any instant of time are ( x 1 , y 1 ) = ( x , 0 ) (x_1,y_1)=(x,0)

Let the distance moved downward by the block along the hypotenuse be s s , then the coordinates of the block are:

( x 2 , y 2 ) = ( h s sin θ , x + s cos θ ) (x_2,y_2)=(h - s \sin{\theta}, x + s \cos{\theta})

The kinetic and potential energies of the system are:

T = 1 2 M ( x ˙ 1 2 + y ˙ 1 2 ) + 1 2 m ( x ˙ 2 2 + y ˙ 2 2 ) T = \frac{1}{2}M\left(\dot{x}_1^2 + \dot{y}_1^2\right) + \frac{1}{2}m\left(\dot{x}_2^2 + \dot{y}_2^2\right) V = m g y 2 + k 2 x 2 V = mgy_2 + \frac{k}{2}x^2

The Lagrangé equations give two linear differential equations as such:

x ¨ = 16 x 5 2 3 \ddot{x} = -\frac{16x}{5}-2\sqrt{3} s ¨ = 8 3 x 5 + 8 \ddot{s} = \frac{8\sqrt{3}x}{5} +8

Now:

x ( 0 ) = s ( 0 ) = x ˙ ( 0 ) = s ˙ ( 0 ) = 0 x(0)=s(0)=\dot{x}(0)=\dot{s}(0)=0

Applying and solving:

x = 5 3 8 ( cos ( 4 t 5 ) 1 ) \boxed{x = \frac{5\sqrt{3}}{8} \left(\cos \left(\frac{4t}{\sqrt{5}}\right) -1\right)}

The solution for s s can be obtained by substituting x x into the equation for s ¨ \ddot{s} and double integrating twice. This gives:

s = 15 cos ( 4 t 5 ) 40 t 2 15 16 \boxed{s = -\frac{15 \cos \left(\frac{4t}{\sqrt{5}}\right) -40t^2-15}{16}}

A plot of s s vs. x x as plotted by @Steven Chase will give the required answer. However, the plot will be scaled by -1 in this case as I chose the axes to be differently oriented.

Solution edited and plot has been added as well.

Karan Chatrath - 11 months, 2 weeks ago

0 pending reports

×

Problem Loading...

Note Loading...

Set Loading...