Register to reply

Help with equations of motion

by boeing_737
Tags: equations, motion
Share this thread:
boeing_737
#1
Jul7-13, 03:43 AM
P: 12
Hi,

I am having a little bit of conceptual trouble with this problem and would appreciate your help. The problem setup is given in the figure. Let's say we have a slender uniform rigid arm(mass m, length l) in space, with a coordinate system [itex]B[/itex] attached to the left end of the arm as shown. C is the center of mass of the arm. We have a moment [itex]M_{z_b}[/itex] acting about the [itex]\hat{z}_{b}[/itex] axis.



Let [itex](u,v,w)[/itex] and [itex](p,q,r)[/itex] be the inertial velocity and inertial angular velocity vectors expressed in [itex]B[/itex]. I get the scalar equations of motion as (assuming that the angular velocity is only along [itex]\hat{z}_b[/itex]):

[itex]m \dot{u} - \frac{ml}{2} r^2 = F_{x_b}[/itex]
[itex]m \dot{v} + \frac{ml}{2} \dot{r} = F_{y_b}[/itex]
[itex]m \dot{w} = F_{z_b}[/itex]

[itex]0 = M_{x_b}[/itex]
[itex]-\frac{ml}{2} \dot{w} = M_{y_b}[/itex]
[itex]\frac{ml^2}{3} \dot{r} + \frac{ml}{2} \dot{v} = M_{z_b}[/itex]

The applied moment is given as : [itex]M_{z_b}(t) = 160 \left(1 - \cos \left(\frac{2 \pi t}{15} \right) \right)[/itex]. For [itex]t > 15, M_{z_b} = 0[/itex]. See figure below :


Integrating these equations using MATLAB's ode45, I get the following plot :


From the above figure :
1) There is only one component of angular velocity (yaw rate) which is as expected. But is the magnitude correct (ie should it reach 24 rad/s)?

2) I am not able to figure out what's going on with u. Why is it increasing so rapidly?

Any help would be really appreciated.

yogesh
Phys.Org News Partner Physics news on Phys.org
Step lightly: All-optical transistor triggered by single photon promises advances in quantum applications
The unifying framework of symmetry reveals properties of a broad range of physical systems
What time is it in the universe?
Simon Bridge
#2
Jul7-13, 10:36 PM
Homework
Sci Advisor
HW Helper
Thanks
Simon Bridge's Avatar
P: 12,961
1. it's your math.
The total change in momentum is the area of the force-time graph ... so you can check.
I don't know why you are not doing this in polar coordinates.

2. I imagine because the moment is quite high.
Have you got any reason to expect u to increase less rapidly?
What sort of value were you expecting and why?


Register to reply

Related Discussions
Need help with motion: equations by rote and applying equations General Physics 2
Equations of motion Introductory Physics Homework 4
6 equations motion Classical Physics 2
Motion Equations Introductory Physics Homework 6
Equations of Motion (Deriving equations) Introductory Physics Homework 1