Torque: using modified balistic pendulum with spring

[aphophis]

Hi everyone i really need help with this, unfortunately my "signals and systems" professor dropped this on us without even covering it in class. From my earlier physics class i can solve for theta in a ballistic pendulum, but since this one has a string attached to it, i really have no clue how to write the torque equation for this.

Please any help or hints will be appreciated. Thanks

Quesion:
http://mywebpage.netscape.com/Drago%20Musilini/untitled1.JPG

A pendulum is initially at rest. It is fixed to a horizontal axis by a rigid rod length of l= 100cm (both masses of the axis and rod are negligible). A spring of torsional stiffness K=196n*m/rad and a torsional dashpot of damping coefficient B=88.54n*m*sec/rad have been connected to the axis. A bullet of mass m=100g traveling with a speed v is fired into the pendulum with a bob mass M = 19.9 kg at time t = 0 and remains lodged therein.

1.Find the differential equation relating the angular displacement θ to the input v. What is the order of the system? Determine the state and output equations in matrix form. Use Matlab to find and plot the unit impulse response and determine the total output response for the inputs of v(t): 200 m/s and 800 m/s
2.Discretize the system and repeat question 1.

 

[Andrew Mason]

What have you got so far?

What is the differential equation of motion without the spring and damping? Then add the spring and damping. Assume the bullet takes a very short time to stop in the block so that the initial angular velocity is:

v_0/R = \dot\theta_0 = \frac{mv_b}{R(m+M)}

AM

 

[aphophis]

so far i was able to get the final velocity of the system as

\frac{M+m}{2g} v_b^2 = (M+m) h where h = l(1-cos θ)

therefore v_b^2 = 2 g l (1-cos theta)

 

[Andrew Mason]

Try analysing the forces on the block tangential to the rod:

Gravity:

F_g = - (M+m)gsin\theta \approx - (M+m)g\theta

Spring:

F_s = - k\theta

Dashpot damper:

F_d = - B\dot\theta

So the torque on the system is:

F \times R = I\alpha = (F_g + F_s + F_d)R = (m+M)R^2\ddot\theta

So set up the differential equation and use the initial condition for v_0 at time t=0.

AM