Práctica 01
Asignatura: Dinámica
Fecha trabajo: 23/10/2024
Fecha de entrega: 23/10/2024
Tiempo planificado en el sílabo: 1 hora
Tiempo de práctica por grupo de estudiantes: No aplica
Número de estudiantes por grupo: Individual
Responsables: Prof. José Cuenca
Resultados de aprendizaje: Calcula y analiza los parámetros cinemáticos en partículas, sistemas de
partículas y sólidos rígidos con movimiento plano general y relativo, organizando racionalmente los procesos
de cálculo y justificando la validez de los resultados, para su aplicación en equipos e instalaciones con
pertinencia social.
Tema
Visualización y análisis de problemas de cinemática de partículas empleando Matlab.
Objetivos
• Conocer la metodología de resolución de ejercicios de movimiento de partículas.
• Conocer el lenguaje de programación de Matlab (y/o Python) para graficar la posición, velocidad y
aceleración de una partícula de un movimiento rectilíneo.
Materiales y reactivos
• Presentaciones en power point sobre el tema de la práctica.
Equipos y herramientas
• Matlab estudiantil
• Google Colab
Preparatorio y marco teórico
1. Movimiento rectilineo uniforme.
2. Movimiento rectilineo uniformemente variable.
3. Gráficas del movimiento: , ,
4. Sistemas de partículas: Análisis de la dinámica de sistemas compuestos por múltiples partículas o
cuerpos.
1
5. Conocimientos básicos de Matlab y Python, incluyendo la sintaxis, manejo de variables, uso de
funciones y comandos para graficar.
Actividades guiadas
A continuación se presenta el desarrollo del Problema 11.61 Pag. 635.
Una partícula se mueve en linea recta con la aceleración que se muestra en la figura. Si se sabe que la
partícula inicia desde el origen con ,
a) construya las curvas y para ,
b) determine la posición y la velocidad de la partícula y la distancia total recorrida cuando .
figure
t1 = linspace(0, 4);
t2 = linspace(4,10);
t3 = linspace(10,20);
a1 = t1*0+3;
a2 = t2*0+6;
a3 = t3*0-5;
plot(t1, a1, 'Color', 'b');
hold on
plot(t2, a2, 'Color', 'b');
plot(t3, a3, 'Color', 'b');
yline(0,'--r','Tiempo (s)')
xline(0,'--m','t=0')
xline(4,'--m','t=4')
xline(10,'--m','t=10')
xline(20,'--m','t=20')
xlabel('Tiempo (s)');
ylabel('Aceleración (ft/s^2)');
title('$a=f\left(t\right)$', 'Interpreter', 'latex', 'FontSize', 15);
2
xlim([-0.5 21.5])
ylim([-7.0 6.7])
grid on
hold off
figure
v1 = 3*t1-18;
v2 = 6*t2-30;
v3 = -5*t3+80;
plot(t1, v1, 'Color', 'b');
hold on
plot(t2, v2, 'Color', 'b');
plot(t3, v3, 'Color', 'b');
yline(0,'--r','Tiempo (s)')
xline(0,'--m','t=0')
xline(4,'--m','t=4')
xline(10,'--m','t=10')
xline(20,'--m','t=20')
xlabel('Tiempo (s)');
ylabel('Velocidad (ft/s)');
title('$v=f\left(t\right)$', 'Interpreter', 'latex', 'FontSize', 15);
xlim([-0.5 21.5])
ylim([-25 35])
3
grid on
hold off
figure
x1 = (3/2).*t1.^2-18*t1+0;
x2 = 3*t2.^2-30*t2+24;
x3 = (-5/2)*t3.^2+80*t3-526;
plot(t1, x1, 'Color', 'b');
hold on
plot(t2, x2, 'Color', 'b');
plot(t3, x3, 'Color', 'b');
yline(0,'--r','Tiempo (s)')
xline(0,'--m','t=0')
xline(4,'--m','t=4')
xline(10,'--m','t=10')
xline(20,'--m','t=20')
xlabel('Tiempo (s)');
ylabel('Posición(ft)');
title('$x=f\left(t\right)$', 'Interpreter', 'latex', 'FontSize', 15);
xlim([-0.5 21.5])
ylim([-60 120])
grid on
4
hold off
Si desea trabajar en Python, el código sugerido lo puede encontrar en el siguiente archivo
de colab_google: https://colab.research.google.com/drive/1ei3Rszz1sIIYQ6MBxHnet17uIy-0e1cO?usp=sharing
(Tiene que ingresar con la cuenta de la universidad)
A continuación se presenta el desarrollo del Ejercicio 11.81
El registro de aceleración que se muestra en la figura se obtuvo durante las pruebas de rapidez de un
automóvil deportivo. Si se sabe que el automóvil inicia desde el reposo, determine de manera aproximada
a) la velocidad del automóvil en t = 8 s,
b) la distancia recorrida por el automóvil en t=20s.
5
Usar https://apps.automeris.io/wpd/ para leer los datos de la gráfica. También se puede emplear https://
plotdigitizer.com/app. Cuando se emplea este último software online, se debe ejecutar el siguiente código.
Borrar la primera fila o encabezado del archivo csv.
filename = 'd1181_2.csv';
delimiter = ',';
fileID = fopen(filename, 'r');
data = textscan(fileID, '%s%s', 'Delimiter', delimiter);
fclose(fileID);
columna1 = data{1};
columna2 = data{2};
t = str2double(columna1);
a = str2double(columna2);
figure
plot(t,a)
hold on
scatter(t,a)
f= fit(t,a,'poly4')
f =
Linear model Poly4:
f(x) = p1*x^4 + p2*x^3 + p3*x^2 + p4*x + p5
Coefficients (with 95% confidence bounds):
p1 = 3.918e-05 (3.357e-05, 4.48e-05)
p2 = -0.0026 (-0.002863, -0.002336)
p3 = 0.06965 (0.06568, 0.07361)
p4 = -0.9546 (-0.9754, -0.9337)
p5 = 6.614 (6.584, 6.644)
plot(f,t,a)
6
ecuacion = formula(f);
disp("Ecuación del ajuste:");
Ecuación del ajuste:
disp(ecuacion);
p1*x^4 + p2*x^3 + p3*x^2 + p4*x + p5
v = zeros(size(t));
x = zeros(size(t));
for i = 2:length(t)
dt = t(i) - t(i-1);
v(i) = v(i-1) + a(i-1) * dt;
x(i) = x(i-1) + v(i-1) * dt + 0.5 * a(i-1) * dt^2;
end
figure;
plot(t, v, 'b-');
xlabel('Tiempo');
ylabel('Velocidad');
grid on;
7
plot(t, x, 'b-');
xlabel('Tiempo');
ylabel('Posición');
grid on;
8
Actividades a desarrollar por los estudiantes
• Desarrollar manualmente los ejercicios 11.63 y 11.65 de la página 635 del libro Beer, F. P.,
& Johnston, E. R. (2021). Mecánica Vectorial para Ingenieros: Dinámica (9ª ed.), con base a la
metodología empleada para resolver el ejercicio 11.61 adjunto.
A continuación se presenta el desarrollo del Ejercicio 11.63
Una partícula se mueve en línea recta con la velocidad que se muestra en la figura. Si se sabe que
en , a) construya las curvas y para , y determine b) la distancia total recorrida por la
partícula cuando , c) los dos tiempos en los cuales
9
a) construya las curvas y para
figure;
t1 = linspace(0, 10);
t2 = linspace(10, 26);
t3 = linspace(26, 41);
t4 = linspace(41, 46);
t5 = linspace(46, 50);
v1 = 60 + 0 * t1;
v2 = linspace(60, -20, numel(t2));
v3 = -20 + 0 * t3;
v4 = linspace(-20, -5, numel(t4));
v5 = -5 + 0 * t5;
plot(t1, v1, 'Color', 'b');
hold on
plot(t2, v2, 'Color', 'b');
plot(t3, v3, 'Color', 'b');
plot(t4, v4, 'Color', 'b');
plot(t5, v5, 'Color', 'b');
yline(0, '--r', 'Tiempo (s)')
xline(0, '--m', 't=0')
xline(10, '--m', 't=10 s')
xline(26, '--m', 't=26')
xline(41, '--m', 't=41 s')
xline(46, '--m', 't=46 s')
xline(50, '--m', 't=50 s')
xlabel('Tiempo (s)');
ylabel('Velocidad (ft/s)');
title('$v=f\left(t\right)$', 'Interpreter', 'latex', 'FontSize', 15);
xlim([-1 51])
ylim([-25 65])
grid on
hold off
10
figure;
a1 = gradient(v1, t1);
a2 = gradient(v2, t2);
a3 = gradient(v3, t3);
a4 = gradient(v4, t4);
a5 = gradient(v5, t5);
plot(t1, a1, 'Color', 'b');
hold on
plot(t2, a2, 'Color', 'b');
plot(t3, a3, 'Color', 'b');
plot(t4, a4, 'Color', 'b');
plot(t5, a5, 'Color', 'b');
yline(0, '--r', 'Tiempo (s)')
xline(0, '--m', 't=0')
xline(10, '--m', 't=10 s')
xline(26, '--m', 't=26 s')
xline(41, '--m', 't=41 s')
xline(46, '--m', 't=46 s')
xline(50, '--m', 't=50 s')
xlabel('Tiempo (s)');
ylabel('Aceleración (ft/s^2)');
11
title('$a=f\left(t\right)$', 'Interpreter', 'latex', 'FontSize', 15);
xlim([-1 51])
ylim([-7 7])
grid on
hold off
figure;
x0 = -540;
x1 = cumtrapz(t1, v1) + x0;
x2 = cumtrapz(t2, v2) + x1(end);
x3 = cumtrapz(t3, v3) + x2(end);
x4 = cumtrapz(t4, v4) + x3(end);
x5 = cumtrapz(t5, v5) + x4(end);
plot([t1, t2(2:end), t3(2:end), t4(2:end), t5(2:end)], [x1, x2(2:end), x3(2:end),
x4(2:end), x5(2:end)], 'Color', 'b');
yline(0, '--r', 'Tiempo (s)')
xline(0, '--m', 't=0')
xline(10, '--m', 't=10 s')
xline(26, '--m', 't=26 s')
xline(41, '--m', 't=41 s')
xline(46, '--m', 't=46 s')
xline(50, '--m', 't=50 s')
12
xlabel('Tiempo (s)');
ylabel('Posición (ft)');
title('$x=f\left(t\right)$', 'Interpreter', 'latex', 'FontSize', 15);
xlim([-1 51])
grid on
13
14
b) la distancia total recorrida por la partícula cuando
La partícula ha recorrido una distancia total de 1382.5 ft al transcurrir los 50 segundos.
c) los dos tiempos en los cuales
Los dos tiempos en los cuales son y .
A continuación se presenta el desarrollo del Ejercicio 11.65
Un paracaidista cae libremente a razón de 200 km/h cuando abre su paracaídas a una altura de 600 m.
Luego de una rápida y constante desaceleración, desciende a una razón constante de 50 km/h desde 586 m
hasta 30 m, donde maniobra el paracaídas en el viento para frenar aún más su descenso. Si se sabe que el
15
paracaidista aterriza con una velocidad descendente despreciable, determine a) el tiempo que requiere para
aterrizar después de abrir su paracaídas, b) la desaceleración inicial.
figure
t1 = linspace(0, 0.403);
t2 = linspace(0.403, 40.43);
t3 = linspace(40.43, 44.8);
x1 = (103.39/2)*t1.^2-5.55*t1+600;
x2 = -13.89*t2+591.598;
x3 = (3.21/2)*t3.^2-143.6703*t3+3215.133;
plot(t1, x1, 'Color', 'b');
hold on
plot(t2, x2, 'Color', 'b');
plot(t3, x3, 'Color', 'b');
xline(0,'--m','t=0')
xline(0.403,'--m','t=0.403')
xline(40.43,'--m','t=40.43')
xlabel('Tiempo (s)');
ylabel('Posición(m)');
title('$x=f\left(t\right)$', 'Interpreter', 'latex', 'FontSize', 15);
xlim([-1 55])
ylim([-1 650])
16
grid on
hold off
figure
v1 = 103.39*t1-55.55;
v2 = -13.89+0*t2;
v3 = 3.21*t3-143.67;
plot(t1, v1, 'Color', 'g');
hold on
plot(t2, v2, 'Color', 'b');
plot(t3, v3, 'Color', 'b');
yline(0,'--r','Tiempo (s)')
xline(0,'--m','t=0')
xline(0.403,'--m','t=0.403')
xline(40.43,'--m','t=40.43')
xlabel('Tiempo (s)');
ylabel('Velocidad (m/s)');
title('$v=f\left(t\right)$', 'Interpreter', 'latex', 'FontSize', 15);
xlim([-1 55])
ylim([-60 3])
grid on
hold off
17
figure
a1 = 103.39+0*t1;
a2 = 0*t2;
a3 = 3.21+0*t3;
plot(t1, a1, 'Color', 'b');
hold on
plot(t2, a2, 'Color', 'b');
plot(t3, a3, 'Color', 'b');
xline(0,'--m','t=0')
xline(0.403,'--m','t=0.403')
xline(40.43,'--m','t=40.43')
xlabel('Tiempo (s)');
ylabel('Aceleración (m/s^2)');
title('$a=f\left(t\right)$', 'Interpreter', 'latex', 'FontSize', 15);
xlim([-1 50])
ylim([-1 10])
grid on
hold off
18
Luego, desarrollar las gráficas solicitadas empleando el software de Matlab, Python u otro software
similar.
• Desarrollar los ejercicios 11.82 y 11.87 de la página 640 del libro Beer, F. P., & Johnston, E. R.
(2021). Mecánica Vectorial para Ingenieros: Dinámica (9ª ed.), con base a la metodología empleada
para resolver el ejercicio 11.81 adjunto.
A continuación se presenta el desarrollo del Ejercicio 11.82
Se requieren dos segundos para dejar en reposo la varilla de un pistón neumático; el registro de aceleración
de la varilla del pistón durante los 2 s es como se indica en la figura. Determine de manera aproximada a)
la velocidad inicial de la varilla del pistón, b) la distancia recorrida por la varilla del pistón antes de quedar en
reposo.
19
filename = '1182.csv';
delimiter = ';';
fileID = fopen(filename, 'r');
data = textscan(fileID, '%s%s', 'Delimiter', delimiter);
fclose(fileID);
columna1 = data{1};
columna2 = data{2};
columna1 = strrep(columna1, ',', '.');
columna2 = strrep(columna2, ',', '.');
t = str2double(columna1);
a = str2double(columna2);
figure
plot(t,a)
hold on
scatter(t,a)
f= fit(t,a,'poly4')
f =
Linear model Poly4:
f(x) = p1*x^4 + p2*x^3 + p3*x^2 + p4*x + p5
Coefficients (with 95% confidence bounds):
p1 = 0.4636 (0.402, 0.5252)
p2 = -2.871 (-3.108, -2.633)
p3 = 6.876 (6.577, 7.175)
p4 = -7.963 (-8.098, -7.827)
p5 = 3.981 (3.964, 3.998)
plot(f,t,a)
20
ecuacion = formula(f);
disp("Ecuación del ajuste:");
Ecuación del ajuste:
disp(ecuacion);
p1*x^4 + p2*x^3 + p3*x^2 + p4*x + p5
v = zeros(size(t));
x = zeros(size(t));
for i = 2:length(t)
dt = t(i) - t(i-1);
v(i) = v(i-1) + a(i-1) * dt;
x(i) = x(i-1) + v(i-1) * dt + 0.5 * a(i-1) * dt^2;
end
figure;
plot(t, v, 'b-');
xlabel('Tiempo');
ylabel('Velocidad');
grid on;
21
plot(t, x, 'b-');
xlabel('Tiempo');
ylabel('Posición');
grid on;
22
A continuación se presenta el desarrollo del Problema 11.87
Durante las pruebas realizadas a una nueva lancha salvavidas, un acelerómetro adherido a la lancha
proporciona el registro que se muestra en la figura. Si la lancha tiene una velocidad de 7.5 ft/s en t 0 y
llega al reposo en el tiempo t1, utilice el método de la sección 11.8 para determinar a) el tiempo t1, b) la
distancia que recorre la lancha antes de quedar en reposo.
23
filename = '1187.csv';
delimiter = ';';
fileID = fopen(filename, 'r');
data = textscan(fileID, '%s%s', 'Delimiter', delimiter);
fclose(fileID);
columna1 = data{1};
columna2 = data{2};
columna1 = strrep(columna1, ',', '.');
columna2 = strrep(columna2, ',', '.');
t = str2double(columna1);
a = str2double(columna2);
figure
plot(t,a)
hold on
yline(0,'--r','Tiempo (s)')
xline(0.75,'--m','t=0.75')
scatter(t,a)
f= fit(t,a,'poly4')
f =
Linear model Poly4:
f(x) = p1*x^4 + p2*x^3 + p3*x^2 + p4*x + p5
Coefficients (with 95% confidence bounds):
p1 = 0.975 (-0.9168, 2.867)
p2 = -15.11 (-26.53, -3.689)
p3 = 72.81 (50.56, 95.07)
p4 = -134.3 (-149.8, -118.7)
p5 = 65.39 (62.28, 68.51)
plot(f,t,a)
24
ecuacion = formula(f);
disp("Ecuación del ajuste:");
Ecuación del ajuste:
disp(ecuacion);
p1*x^4 + p2*x^3 + p3*x^2 + p4*x + p5
v = zeros(size(t));
x = zeros(size(t));
for i = 2:length(t)
dt = t(i) - t(i-1);
v(i) = v(i-1) + a(i-1) * dt;
x(i) = x(i-1) + v(i-1) * dt + 0.5 * a(i-1) * dt^2;
end
figure;
plot(t, v, 'b-');
xlabel('Tiempo');
ylabel('Velocidad');
grid on;
yline(0,'--r','Tiempo (s)')
xline(0.75,'--m','t=0.75')
25
plot(t, x, 'b-');
xlabel('Tiempo');
ylabel('Posición');
yline(0,'--r','Tiempo (s)')
xline(0.75,'--m','t=0.75')
grid on;
26
Resultados Obtenidos
Los resultados serán los siguientes:
• Ejercicios resueltos empleando Matlab como herramienta de análisis y graficación, y
• Gráficas del movimiento: , , ejecutadas en Matlab que facilitan la comprensión de los
distintos movimientos rectilíneos: MR y MRUV.
Discusión
27
En esta práctica, el uso de Matlab facilitó el análisis de los parámetros cinemáticos de manera precisa. Las
herramientas de Matlab permitieron representar gráficamente la posición, velocidad y aceleración de una
partícula, lo que contribuyó a una mejor comprensión del movimiento rectilíneo uniforme y uniformemente
acelerado.Las fortalezas del uso de Matlab radican en su facilidad para realizar cálculos numéricos y generar
gráficos, lo que agilizó el análisis.
Conclusiones
• A través del desarrollo de esta práctica, se logró aplicar una metodología clara y precisa para resolver
problemas de cinemática de partículas.
• El uso de Matlab en esta práctica ha permitido aprovechar su fortaleza en la manipulación de datos y
generación de gráficos, facilitando el análisis dinámico del sistema.
Recomendaciones
No aplica
Preguntas de control
• No aplica.
Bibliografía
• Beer, F. P., & Johnston, E. R. (2021). Mecánica Vectorial para Ingenieros: Dinámica (9ª ed.). Editorial
McGraw-Hill.
• https://colab.research.google.com/drive/1ei3Rszz1sIIYQ6MBxHnet17uIy-0e1cO?usp=sharing
• https://drive.google.com/file/d/1ptRhufJQMAKaN7xM-2k9rSUTKdhqtX_u/view?usp=sharing
28
