Problemas de física resueltos usando Matlab
Ejercicio1: cajas y fricción
clc;
clear all;
close all;
% Definición de parámetros
m =100; % kg
g=9.81; % m/s^2
L = 0.5; % m
H =1.5; % m
F_max = 100; % Fuerza máxima, N
mu=0.2; % Coeficiente de fricción sin unidades
N =m *g; % Fuerza normal, N
% Cálculo inicial para F = 100 N
x= F_max * H / (m *g) +L / 2; % m
if F_max <= mu *N
if x >L
str="vuelca";
else
str= "no se vuelca";
end
else
if x> L
str ="vuelca";
else
str ="arrastramos la caja sin volcar";
end
end
% Mostrar resultado para el caso inicial
disp(['si aplicamos esta f: ', num2str(F_max), ' resultado: ', str]);
% vectores de fuerza
F = linspace(0, F_max, 101);
x = F *H / (m *g) + L/2; % Vector de posiciones x para cada F
% Gráfica
figure;
plot(F, x,'b', 'LineWidth', 2);
hold on;
yline(L, 'r--','LineWidth', 2); % Línea de límite de vuelco
title('Análisis de estabilidad de la caja');
EXPLICACIÓN:
En este ejercicio calculamos si la caja
rectangular puede volcar, para
comprobar eso cambiamos la cantidad
de fuerza que le ejercemos.
De este modo nos damos cuenta que
cuanta mas fuerza ejerzamos la
posibilidad de que la caja vuelque es mayor.
EJ 2
% Datos del problema
m = 50; % masa del candelabro en kg
g = 9.81; % gravedad en m/s^2
W = m * g; % peso del candelabro en Newtons
% Ángulos en grados
angulo_BD = 45;
angulo_CD = 30;
angulo_AB = 30; % Ángulo en el cable AB
angulo_BC = 45; % Ángulo en el cable BC
% Calcular cosenos y senos de los ángulos
cos_30 = cosd(30);
sin_30 = sind(30);
cos_45 = cosd(45);
sin_45 = sind(45);
% 1. Resolver las tensiones en BD y CD (equilibrio en el punto D)
% Ecuaciones:
% En x: T_BD * cos(45) = T_CD * cos(30)
% En y: T_BD * sin(45) + T_CD * sin(30) = W
% Resolviendo para T_CD y T_BD
T_CD = W / (sin_45 + sin_30 * cos_45 / cos_30); % Despejamos T_CD
T_BD = T_CD * cos_30 / cos_45; % Despejamos T_BD
% Mostrar los resultados de T_BD y T_CD
fprintf('Tensión en el cable BD: %.2f N\n', T_BD);
fprintf('Tensión en el cable CD: %.2f N\n', T_CD);
% 2. Resolver las tensiones en AB y BC (equilibrio en el punto B)
% Ecuaciones:
% En x: T_AB * cos(30) = T_BC * cos(45) + T_BD * cos(45)
% En y: T_AB * sin(30) + T_BC * sin(45) = T_BD * sin(45)
% Despejar T_BC usando la ecuación en y
T_BC = (T_BD * sin_45 - T_AB * sin_30) / sin_45;
% Despejar T_AB usando la ecuación en x
T_AB = (T_BC * cos_45 + T_BD * cos_45) / cos_30;
% Mostrar los resultados de T_AB y T_BC
fprintf('Tensión en el cable AB: %.2f N\n', T_AB);
fprintf('Tensión en el cable BC: %.2f N\n', T_BC);
Explicación del Código:
1. Declaración de Variables Simbólicas: Utilizamos syms para definir las tensiones TAB,
TBC, TBD, y TCD como variables simbólicas.
2. Definición de Constantes: Calculamos los valores de las funciones trigonométricas
para 30º y 45º, y definimos el peso W.
3. Equilibrio en el Punto D: Establecemos dos ecuaciones de equilibrio en x e y en el
punto D y resolvemos para TBD y TCD.
4. Equilibrio en el Punto B: Usamos los valores obtenidos de TBD y TCD para
establecer las ecuaciones de equilibrio en el punto B y resolver para TAB y TBC.
5. Mostrar Resultados: Imprimimos los valores finales de todas las tensiones en los
cables.
EJ3 momentos
% Datos del problema
F = 6; % Fuerza en kN
a = 3; % Longitud en m
% Descomposición de la fuerza en componentes horizontales y verticales
FH = F * cosd(θ); % Componente horizontal
FV = F * sind(θ); % Componente vertical
% Momentos respecto al punto A
M_FA = FH * a; % Momento debido a la componente horizontal
M_FV = FV * a; % Momento debido a la componente vertical
% Ecuación del momento total
syms θ
eqn = -36 * sind(θ) + 18 * cosd(θ) == 0;
% Resolviendo para el ángulo θ
sol = solve(eqn, θ);
sol = double(sol);
% Mostrando la solución que está en el rango correcto (0 a 360 grados)
sol_valid = sol(sol >= 0 & sol <= 360);
disp('El ángulo θ que satisface la condición de equilibrio es:')
disp(sol_valid)
Explicación del Código
1. Se define una fuerza F de 6 kN.
2. La fuerza F se descompone en sus componentes horizontal (FH) y vertical (FV)
usando trigonometría.
3. Se usa cosd(45) y sind(45) porque el ángulo está en grados. cosd y sind
calculan el coseno y el seno respectivamente en grados.
4. Aquí se calcula el momento que cada componente de la fuerza (FH y FV) genera
respecto a un punto A.
5. Se supone que las distancias desde A hasta las líneas de acción de FH y FV.
EJ 4:RESORTES CON PESOS
clear all
clc
close all
% Definición de constantes
k = 500; % N/m
x = 0.3; % m
d = 2; % m
g = 9.81; % m/s^2
% Cálculo de la fuerza en el resorte
Fs = k * x;
fprintf('Fuerza en el resorte, Fs = %.2f N\n', Fs);
% Cálculo del ángulo theta en grados
theta = asind(0.7 / 1.4);
fprintf('Ángulo θ = %.2f grados\n', theta);
% Cálculo de la tensión en el cable
T = Fs / (2 * sind(theta));
fprintf('Tensión en el cable, T = %.2f N\n', T);
% Cálculo de la masa de cada bloque
m = T * cosd(theta) / g;
fprintf('Masa de cada bloque, m = %.2f kg\n', m);
Explicación del Código
1. Definición de constantes: Se define la rigidez del resorte k, la elongación x, y la
gravedad g.
2. Cálculo de la fuerza en el resorte (Fs):
Fs=k×xFs = k \times xFs=k×x
Esto determina la fuerza que ejerce el resorte debido a su elongación.
3. Cálculo del ángulo theta: Usando trigonometría, se obtiene el ángulo theta en
grados que corresponde a una relación de lados dada.
4. Cálculo de la tensión en el cable (T):
T=Fs2×sin(θ)T = \frac{Fs}{2 \times \sin(\theta)}T=2×sin(θ)Fs
Este valor representa la fuerza que debe soportar cada cable.
5. Cálculo de la masa de cada bloque (m):
m=T×cos(θ)gm = \frac{T \times \cos(\theta)}{g}m=gT×cos(θ)
Calcula la masa que cada bloque debería tener para equilibrar la tensión en el
sistema.
