EX Pablo Gumiel
EX Pablo Gumiel
EX Pablo Gumiel
1.- En la tabla se muestran los esfuerzos y deformaciones de un ensayo a tracción de una probeta. Hallar a)
Limite de Proporcionalidad, b) Modulo de Elasticidad c) Limite de Fluencia (,002%) y d) Limite de Rotura
Área 0,3
Longitud Inicial 15
ESFUERZO-DEFORMACION
40000
0,01; 36333,33333
0,013333333; 36666,66667
35000
0,02; 33333,33333
0,006666667; 31000
30000
25000
0,004666667; 23333,33333
ESFUERZO
20000
0,003333333; 16666,66667
15000
5000
0,000666667; 3333,333333
0 0; 0
0 0,005 0,01 0,015 0,02 0,025
DEFORMACION UNITARIA
a) Limite de Proporcionalidad
Mecánica del medio continuo
Pablo Andrés Gumiel Rivera
LIMITE DE PROPORCIONALIDAD
25000 y = 5E+06x
R² = 1 0,004666667;
23333,33333
20000
0,003333333;
16666,66667
TÍTULO DEL EJE
15000
5000
0,000666667;
3333,333333
0 0; 0
0 0,0005 0,001 0,0015 0,002 0,0025 0,003 0,0035 0,004 0,0045 0,005
TÍTULO DEL EJE
El punto donde la gráfica deja un comportamiento lineal se puede analizar en base a una regresión
lineal.
𝜎𝑝 = 23333,33 𝐾𝑔/𝑐𝑚2
b) Modulo de Elasticidad
23333,33 − 0
𝐸= = 5.000.000,00𝐾𝑔/𝑐𝑚2
0,0046 − 0
c) Limite de Fluencia (0,2%)
ESFUERZO-DEFORMACION
60000
50000
40000
0,01; 36333,33333
0,013333333; 36666,66667
ESFUERZO
0,02; 33333,33333
30000 0,006666667; 31000
0,004666667; 23333,33333
20000
0,003333333; 16666,66667
0,000666667; 3333,333333
0 0; 0
0 0,005 0,01 0,015 0,02 0,025
DEFORMACION UNITARIA
d) Limite de Rotura
= 33333,33𝐾𝑔/𝑐𝑚2
Mecánica del medio continuo
Pablo Andrés Gumiel Rivera
2.- La matriz de esfuerzos de un material en un punto, referida a cierto sistema de referencia OXYZ, es:
1 0 −2
= 0 3 1
−2 1 −2
1 − 0 −2
0= 0 3− 1
−2 1 −2 −
hallando autovalores
esfuerzos
1 = -3.1313
2 = 1.8481
3 = 3.2832
hallando autovectores
0, 4311 -0.8719 0.2322
-0.1452 -0.3210 -0.9359
0.8905 0.3697 -0.2650
direcciones principales
u 1 (0, 4311;-0.1452;0.8905)
u 2 (-0.8719;-0.3210;0.3697)
u 3 (0.2322;-0.9359;-0.2650)
3 = 3.2832
2 = 1.8481
1 = -3.1313
Usando la ecuación de las cuádricas indicatrices
3 x2 + 2 y 2 + 1z 2 = 1
De la condición que dos esfuerzos principales son positivos y otra negativa
3 0; 2 0;1 0;
Ambos hiperboloides tienen el mismo cono asintótico
Las direcciones que solo experimentan esfuerzos transversal son las que coinciden con las generatrices del
cono asintótico.
( , , ) son las componentes del vector que definen la dirección que solo presenta esfuerzos trasversal
3 x2 + 2 y 2 + 1z 2 = 0
3.2832x 2 + 1.8481y 2 -3.1313z 2 = 0
b)
% b) Calcular el esfuerzo en un plano inclinado cuyo vector normal es (1,1,1)
% Usar la fórmula del esfuerzo normal: sigma_n = n' * q * n
n = [1;1;1] % Vector normal del plano
n = n / norm(n)% Normalizar el vector
sigma_n = n' * q * n % Esfuerzo normal en el plano
Mecánica del medio continuo
Pablo Andrés Gumiel Rivera
3.- En la figura se muestra una viga simplemente apoyada con un momento en su parte central
( x, y ) = a1 x 2 + a2 xy + a3 y 2 + a4 x 3 + a5 x 2 y + a6 xy 2 + a7 y 3
GRADO 3
2 3
( x, y ) = a1 x + a2 xy + a3 y + a4 x + a5 x y + a6 xy 2 + a7 y 3
2 2 3 2
( x, y ) = a1 x 2 + a2 xy + a3 y 2 + a4 x 3 + a5 x 2 y + a6 xy 2 + a7 y 3
4 4 4
= 4 + 2 2 2 + 4 = 0 → 0 + 2 (0) + 0 = 0
4
x x y y
2
x = = 2a3 + 2a6 x + 6a7 y
y 2
2
y = 2 = 2a1 + 6a4 x + 2a5 y
x
2
xy = − = −(a2 + 2a5 x + 2a6 y )
xy
Mecánica del medio continuo
Pablo Andrés Gumiel Rivera
DERECHA
L
xx ( x , y ) = 0
2a3 + 2a6 ( L) + 6a7 y = 0
a7 = 0
a3 = −a6 ( L)
H
Arriba yy ( x, y ) = 0
2
y = = 2 a1 + 6 a4 x + 2 a5 y = 0 = 2a1 + 6a4 x + 2a5 H
x 2
0 = 2a1 + 6a4 x + 2a5 H (i )
−H
Abajo yy ( x, y ) = 0
0 = 2a1 + 6a4 x − 2a5 H (ii )
sumamosi + ii
a1 + 3 a4 x = 0
a1 = 0 ; a4 = 0
restamosi − ii
4 H a5 = 0
a5 = 0
H
xy ( x, y ) = 0 (2
−(a2 + 2a6 H ) = 0
a2 + 2a6 H = 0
a2 = −2a6 H
H L 1
M
Py =
−H
xy ( x , y ) dz dy = −
L
H
M
Py =
−H
2 a ( H − y )dy = − L
6
H
M
(2 a Hdy − 2 a ydy) = − L
−H
6 6
y2 M
2 a6 Hy − 2 a6 H
−H
H
−H =−
2 L
M
4 a6 H 2 = −
L
M
a6 = −
4 LH 2
REMPLAZANDO
Mecánica del medio continuo
Pablo Andrés Gumiel Rivera
M
a3 = −a6 ( L) =
4H2
M
a2 = −2a6 H =
2H L
2 M ( L − x ) M Mx
x = = = −
y 2 2
2H L 2H 2
2H2 L
2
y = 2 = 0
x
2 M (H − y) My M
xy = − =− = −
xy 2H L2 2
2H L 2H L
M y2 M x y2 M x y
( x, y ) = − +
4H2 4H2 L 2H L
2 M ( L − x ) M Mx
x = = = −
y 2 2
2H L 2H 2
2H2 L
2
y = =0
x 2
2 M (H − y) My M
xy = − =− = −
xy 2H L2 2
2H L 2H L
M y2 M x y2 M x y
( x, y ) = − +
4H2 4H2 L 2H L
bh3 h3
La inercia es I= y en vista que b=1 → I= → 12I=H 3
12 12
2 MH MH x MH x
x = 2 = − = (1 − )
y 24I 24I L 24I L
2
y = 2 = 0
x
2 M Hy MH 2 M H 2 y
xy = − = − = ( − 1)
x y 24 I L 24I L 24I L H
MH y 2 M Hx y 2 MH 2 x y
( x, y ) = − +
48I 48I L 24I L
clc
clear
syms y x H L M P a1 a2 a3 a4 a5 a6 a7 a8 a9 a10 a11 a12 a13 a14 a15 a16 a17 a18
%ARRIBA H(y) qyy
a=H
%sigma yy valor de arriba
b=0
%sigma yy valor de abajo
Mecánica del medio continuo
Pablo Andrés Gumiel Rivera
c=0
%arriba H(Y) Txy
%tau xy valor de arriba
d=0
%tau xy VALOR abajo
e=0
%cambia en funcion de derecha-izquerda
%Px-qxx valor de x
g=L
%px=
f=0
%Py-tau xy valor de x
j=L
%py=
k=-M/L
%momento valor de x
n=L
%Mz=
r=0
%limites de la intergral
ls=H
li=-H
%tercer grado
a8=0;a9=0;a10=0;a11=0;a12=0;a13=0;a14=0;a15=0;a16=0;a17=0;a18=0;
a1=0
a2=-2*a6*H
a3=-a6*L
a4=0
a5=0
a6=-M/(4*L*H^2)
a7=0
a2=-2*a6*H
a3=-a6*L
q2=a1*x^2+a2*x*y+a3*y^2;
q3=a4*x^3+a5*x^2*y+a6*x*y^2+a7*y^3;
q4=a8*x^4+a9*x^3*y+a10*x^2*y^2+a11*x*y^3+a12*y^4;
q5=a13*x^5+a14*x^4*y+a15*x^3*y^2+a16*x^2*y^3+a17*x*y^4+a18*y^5;
q=q2+q3+q4+q5
%funcion air
d4f_dx4 = diff(q, x, 4)
d4f_dy4 = diff(q, y, 4)
d4f_dx2dy2 = diff(diff(q, y, 2), x, 2)
air=simplify(d4f_dx4+2*d4f_dx2dy2+d4f_dy4==0)