ex-5-9-meriam.sce

theta = 60
beta = 18.02

// Omega en radianes
w = 1500*2*%pi/60


// Analisis de velocidades
//
// V_A = V_B + V_A/B
//
// Omega vector (solo tiene componente en -z por regla de la mano derecha)
W = [0 0 -w]

// Vector r_BO en ft
r_BO = 5/12*[-cosd(theta) sind(theta) 0]

// Vector r_AB en ft
r_AB = 14/12*[-cosd(beta) -sind(beta) 0]

//Vector velocidad V_B en ft/seg
V_B = cross(W,r_BO)

// Definicion de variables simbolicas
// velocidad de A como variable (desconocida)
va = poly(0,'va')
// velocidad angular de AB
wab = poly(0,'wab')
// wAB como vector (solo tiene componente en z)
W_AB = [0 0 1]*wab

//Vector V_A
// V_A solo tiene componente en la direccion x (i)
V_A = [1 0 0 ]*va

//Velocidad relativa de A con respecto a B
V_AB = cross(W_AB,r_AB)

//Igualando en y, se tiene que wab = 29.5 rad/seg
// Igualando en x, se tiene que va = 67.3 ft/sec
disp(V_B+V_AB,'=',V_A)
//disp('es igual a:)
//disp(V_B+V_AB)

wab = 32.72/1.109
va = 56.68+0.36*(29.5)

//Se vuelven a calcular los vectores W_AB, V_A y V_AB (ya no son incognitas)
V_A = [1 0 0 ]*va
W_AB = [0 0 1]*wab
V_AB = cross(W_AB,r_AB)


// Analisis de aceleraciones
//
// A_A = A_B + A_A/B
//
// Definicion de variables simbolicas
// aceleracion de A como variable (desconocida)
aA = poly(0,'aA')
// aceleracion angular de AB
alfab = poly(0,'alfab')
// alfa_AB como vector (solo tiene componente en z)
Alfa_AB = [0 0 1]*alfab
// aceleracion de A como vector (solo tiene componente en x)
A_A = [1 0 0]*aA

//aceleracion de B, (solo tiene componente normal)
A_B = cross(W,V_B)

//aceleracion de A/B, tiene comp normal y tangencial
A_AB_n = cross(W_AB,V_AB)
A_AB_t = cross(Alfa_AB,r_AB)

disp(A_B+A_AB_n+A_AB_t,'=',A_A)

//
//Igualando en y, se tiene que alfab = 7742 rad/seg^2
// Igualando en x, se tiene que aA = 3312 ft/sec^2
alfab = 8589.301/(-1.10944)
aA = 6106.1748 + 0.3609071*(alfab)