EX Pablo Gumiel

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 9

Mecánica del medio continuo

Pablo Andrés Gumiel Rivera


EXAMEN DE MECÁNICA DEL MEDIO CONTINUO

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

Deformación Deformación Esfuerzo


Carga (kg)
(cm) Unitaria (kg/cm2)
0 0 0 0
1000 0,01 0,000666667 3333,333333
3000 0,03 0,002 10000
5000 0,05 0,003333333 16666,66667
7000 0,07 0,004666667 23333,33333
9300 0,1 0,006666667 31000
10900 0,15 0,01 36333,33333
11000 0,2 0,013333333 36666,66667
10000 0,3 0,02 33333,33333

Á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

10000 0,002; 10000

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

10000 0,002; 10000

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

10000 0,002; 10000

0,000666667; 3333,333333
0 0; 0
0 0,005 0,01 0,015 0,02 0,025
DEFORMACION UNITARIA

Limite de Fluencia (0,2%) = 35000 𝐾𝑔/𝑐𝑚2

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:

a) hallar las esfuerzos y direcciones principales

 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)

b) calcular el esfuerzo en un plano inclinado cuyo vector normal es (1,1,1)

b) calcular el esfuerzo en un plano inclinado cuyo vector normal es (1, 1, 1)


1  1 
1  
  3
1  1   1 1 1 
u= =   ut =  
12 + 12 + 12  3   3 3 3
 1 
 3
 
Esfuerzo longitudinal Unitaria
 n = ut * * u
 1 
 
3
 1 0 −2  
 1 1 1     1 
n =  * 0 3 1 *
 3 3 3 3 
 −2 1 −2  
 1 
 3 

 n = 5,5511*10−3
Mecánica del medio continuo
Pablo Andrés Gumiel Rivera
c) las direcciones que sólo experimentan esfuerzo transversal.

Ordenamos los esfuerzos principales

 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

Código Matlab para resolver


a)
clear
clc
% Definir la matriz de tensiones
q = [1 0 -2;0 3 1;-2 1 -2]
% a) Hallar las tensiones y direcciones principales
% Usar la función eig para obtener los valores y vectores propios
[t,v] = eig(q)
[V, D] = eig(q);
% Mostrar resultados
disp('Autovalores:');
disp(D);
disp('Autovectores:');
disp(V);

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

Para una de las mitades se pide: a) Las condiciones de Borde

 x = f1 ( x, y ), xy = f 2 ( x, y ), y = f 3 ( x, y ) Hallar funciones desconocidas


(No hay simetria horizontal)

 H H
Arriba y Abajo  y ( x, y ) = 0 (1, xy ( x, y ) = 0 (2

 L L 1 L 1 L 1

H H H
M
CB  Derecha  x ( x , y ) = 0 Px =   x ( x , y ) dz dy = 0 (3, Py =   xy ( x , y ) dz dy = − (4, M z =   x ( x , y ) y dz dy = 0 (5
 −H −H L −H

 H 0 1 H 0 1 H 0 1

I zquierda Px =   x ( x , y ) dz dy = 0 (6, Py =   xy ( x , y ) dz dy = 0 (7, M z =   x ( x , y ) y dz dy = M (8


 −H −H −H

SOLO TOMAREMOS DERECHA

a) La función polinómica de Airy

 ( 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 )
 xy
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 = − =− = −
 xy 2H L2 2
2H L 2H L

M y2 M x y2 M x y
 ( x, y ) = − +
4H2 4H2 L 2H L

b) La expresión de los esfuerzos

  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 = − =− = −
 xy 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

Código Matlab para simplificar solución

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)

disp('esfuerzos SIN SIMPLIFICAR ')


qxx=(diff(q, y, 2))
qyy=(diff(q, x, 2))
txy=-(diff(diff(q, y, 1), x, 1))
disp('esfuerzos')
qxx=simplify(diff(q, y, 2))
qyy=simplify(diff(q, x, 2))
txy=-simplify(diff(diff(q, y, 1), x, 1))
Mecánica del medio continuo
Pablo Andrés Gumiel Rivera
disp('izquierda esfuerzo qxx X=f')
disp('ver el valor de f')
% Evalúa qxx en x=g
L = sym('L');
qxx_remplaza = simplify(subs(qxx, x, (g)))
qxx_remplaza=(qxx_remplaza==f)
disp('arriba qyy Y=a')
% Evalúa qyy en y=a
H = sym('H');
disp('REMPLAZANDO a')
qyy_remplaza = simplify(subs(qyy, y, (a)))
qyy_remplaza=(qyy_remplaza==b)
i=qyy_remplaza
disp('aBAJO qyy Y=-a')
disp('REMPLAZANDO -a')
H = sym('H');
qyy_remplaza = simplify(subs(qyy, y, (-a)))
qyy_remplaza=(qyy_remplaza==c)
ii=qyy_remplaza
disp('suma qyy si es igual a cero i + ii')
i+ii
simplify(i+ii)
disp('resta qyy si es igual a cero i - ii')
i-ii
simplify(i-ii)
disp('CORTANTE ARRIBA txy Y=a')
% Evalúa arriba qxy en y=a
H = sym('H');
disp('txy REMPLAZANDO H')
txy_remplaza = simplify(subs(txy, y, (a)))
txy_remplaza=(txy_remplaza==d)
iii=txy_remplaza

% Calcular la integral definida con respecto a y izquierda


disp('txy REMPLAZANDO izquierda PY X=j')
L = sym('L');
txy_remplaza = subs(txy, x, (j))
Py = int(txy_remplaza, y, (li), (ls))
Py =(int(txy_remplaza, y, (li), (ls))==k)
L = sym('L');
disp('qxx MOMENTO REMPLAZANDO Mz X=n')
qxx_remplaza = y*subs(qxx, x, (n))
Mz=int(qxx_remplaza, y, (li), (ls))
Mz = (int(qxx_remplaza, y, (li), (ls))==r)
disp('qxx REMPLAZANDO PX X=g')
L = sym('L');
qxx_remplaza = subs(qxx, x, (g))
Px = int(qxx_remplaza, y, (li), (ls))
Px = (int(qxx_remplaza, y, (li), (ls))==f)

También podría gustarte