Instituto Tecnolgico y de Estudios Superiores de

Monterrey
Campus Ciudad de Mxico

PROYECTO SEGUNDO PARCIAL

PUENTE ATIRANTADO
ENERGIA

Mtodos Avanzados para la Resistencia de Materiales

Dr. Diego Crdenas Fuentes

EQUIPO:
Oscar Jess Torres Contreras A01152343

Objetivo: Analizar el puente del primer parcial, en este caso se analizara la

energa que puede absorber en un caso de emergencia.
Problemtica: Calcular la cantidad de energa que puede absorber el puente
que se diseo durante el primer parcial, la situacin a analizar ser en el
momento en el cual uno de los camiones suelta la carga que transporta desde
una altura establecida.
Requerimientos del diseo: En este caso solamente se buscara analizar y
obtener valores numricos, no importa que el puente pueda llegar a sufrir
rupturas o algo por el estilo, solamente se busca conocer que es lo que pasara
en esta situacin.

Imagen 1.1. Diagrama del planteamiento de la problemtica a analizar.

Imagen 1.2. Ejemplo de la problemtica planteada.

.5m

7m
Imagen 1.3. Diagrama de dimensiones del camin y carga.

500kg

1m
Imagen 1.4. Diagrama de dimensiones del la caja.

VISTA FRONTAL

VISTA DESDE ARRIBA

5m

20m

Imagen 1.5. Dimensionamiento del puente.

MATERIALES.

PTR (Travesaos):
PTR Rojo
Dimensiones: .1016m x .204m
Espesor: .00631m
Longitud: 5m
Momento de Inercia: 10x 106 m4
Modulo de elasticidad: 200 GPa.

VIGA I:
Kilos: 37.00 kg/m
Longitud: 20000 mm.
Seccin en mm. "h" Altura: 900 mm
Seccin en mm. "b" Ancho: 155 mm
Seccin en mm. "e": 10.6 mm
Seccin en mm. "e1": 12.20 mm
Calidad: S275JR
Modulo de elasticidad: 190-210 GPa.
Imagen 1.7. Diagrama y dimensionamiento Viga I.

CABLES DE ACERO:
Acero Inoxidable AISI316
Resistencia: 1570 N/ mm2
Modulo Elastico: 130000 N/ mm2
Acabado Pulido.
D nominal: 5.00cm
Peso: .005 kg/m
Seccin: .60 mm2
Carga Rotura: .82 kN

DESARROLLO 1.1

Para este desarrollo no se tomara en cuenta el peso del camin,

primero se buscara obtener la sumatoria de las energas, para
de esta manera obtener una relacin entre la cantidad de
energa que absorbe el puente y sus elementos y la w que
puede soportar.

Daremos valores predeterminados para el anlisis de energa,

se decidi utilizar un peso de 500kg que al multiplicarse por la
gravedad se convierte a Newtons, y as mismo se decidi
establecer el tamao de la caja y dividir el resultado de los

Newtons entre esa dimensin, esto con el fin de establecer una

carga distribuida.
Peso= 500kg
a= 1m
w 1=

5009.81 N
=4905 N / m
1m
Primero empezaremos analizando el Travesao, para despus
analizar la reaccin que este elemento genera en el Larguero y
de la misma forma en los Cables. Se realizaran diagramas de
cuerpo libre y sumatorias de fuerzas en combinacin con las
ecuaciones de deflexin, para terminar con las ecuaciones de la
energa.

TRAVESAO 1.1:
Diagrama de Cuerpo Libre Travesao.
a=1m
b=2m
w=4905N/m
W

Ra

5m

Sumatoria de fuerzas y momentos del Travesao:

Fy=0
Ra=

wa
2

Rb

Rb=

wa
2

M =0

( a2 )+ RbMb=0

wa b+

En este caso las reacciones salen de manera directa con las

sumatorias de fuerzas realizadas, pero para corroborar los
resultados se metieron los datos y se realizaron las siguientes
operaciones matriciales en Matlab:

clearall

E=200e9;
L=5;
h=0.204;
h1=1;
w=(500*9.81)/1;
a=1;
b=2;
I=10*10^6;
x=(0:.01:L);

A=[1100
0L11
L^2/20L0
L^3/60L^2/20];

B=[w;(w/2*(Lb)^2w/2*(L1b)^2);(w/6*(Lb)^3w/6*(L1b)^3);(w/24*(L
b)^4w/24*(L1b)^4)];

res=inv(A)*B;
Ra=res(1);
Rb=res(2);
Ma=res(3);
Mb=res(4);

Los resultados de las operaciones pasadas dieron los siguientes

resultados:

Ra= 2.4525e3 N
Rb= 2.4525e3 N
Ma= 3.0248e3 N
Mb= 3.0248e3 N

Con los resultados obtenidos podemos pasar a plantear la

ecuacin de momentos internos de nuestro travesao:

M ( x )=Ma+ Rax

w
2 w
2

x b + xab
2
2

Teniendo la ecuacin de momentos, integramos una ves para

obtener la ecuacin del ngulo de deflexin.

( x ) =Max +

3 w
3
Rax 2 w
xb + x ab
2
6
6

Ahora volvemos a integrar la ecuacin para obtener la

educacin de deflexin.
2

4
4
Max Rax
w
w
+
xb + xab
2
6
24
24
Las mismas ecuaciones las planteamos en el programa de
Matlab y nos dan los siguientes resultados:

V (x )=

Mmax= 1.5328e4 Nm

Angulo Deflexin max= 9.32e-4

Deflexin Max= .0015m

Continuamos con la ultima parte que es sacar la energa total

que el travesao absorbe:
L

U=
0

M (x)
dx
2 EI

El programa de Matlab es relativamente sencillo y es lo

siguiente:

M2=Mx.^2;
intM2=trapz(x,M2);
U=1/(2*E*I)*intM2;

U=199.5847 Joules
LARGUERO 1.1:
Diagrama de Cuerpo Libre Larguero.
L=20m
Lc1=6m
Lc2=9m
Lc3=12m
Ra= 2.4525e3 N
T1

T2

60

M2
M2

Lc1
R2

T3

60

Lc2

60

Ra
20m

M3
Lc3

R2
Debido a la complejidad del planteamiento de este elementoR3
se
requiri meter todos los datos de dimensiones y materiales al
programa Matlab y de esa manera poder tener lo necesario
para resolver la problemtica. Primero se sacan las reacciones
por medio de matrices:

clearall
Ac=.001963;
E=210e9;
Ec=130e9;
L=20;
l=10;
x1=4;
x2=8;
x3=10;
x4=14;
lc1=6;
lc2=9;
lc3=12;
th3=60;
th2=60;
th1=60;
h=0.3048;
I=0.000010338239;
Ra=2.4525e3;
x=(0:.01:L);

A=[01sind(th1)sind(th2)sind(th3)10
10lc1*sind(th1)lc2*sind(th2)lc3*sind(th3)L1
LL^2/2(((Llc1)^2)/2)*sind(th1)(((Llc2)^2)/2)*sind(th2)(((L
lc3)^2)/2)*sind(th3)00
L^2/2L^3/6(((Llc1)^3)/6)*sind(th1)(((Llc2)^3)/6)*sind(th2)
(((Llc3)^3)/6)*sind(th3)00
lc1^2/2lc1^3/6(sind(th1)*lc1)/(Ec*Ac)0000
lc2^2/2lc2^3/6(((lc2lc1)^3)/6)*sind(th1)(sind(th2)*lc2)/
(Ec*Ac)000
lc3^2/2lc3^3/6(((lc3lc1)^3)/6)*sind(th1)(((lc3
lc2)^3)/6)*sind(th2)(sind(th3)*lc3)/(Ec*Ac)00];

B=[Ra;(l*Ra);(Ra/2*(Ll)^2);(Ra/6*(Ll)^3);(Ra/6*(lc1l)^3);
(Ra/6*(lc2l)^3);(Ra/6*(lc3l)^3)];

res=inv(A)*B;

M2=res(1);
R2=res(2);
T1=res(3);
T2=res(4);
T3=res(5);
R3=res(6);
M3=res(7);

Los resultados de las operaciones pasadas dieron los siguientes

resultados:

R2= 5.99e3 N
M2= 8.107e3 Nm
T1= 3.0248e3 N
T2= 1.039e4 N
T3= 2.18e3 N
R3= 140.62 N
M3= 375 Nm

Con los resultados obtenidos podemos pasar a plantear la

ecuacin de momentos internos de nuestro travesao:

M ( x )=Ma+ Rax+T 1 sen 60 x lc 1 +T 2 sen 60 xlc2 +T 3 sen 60 xlc 3 Ra x10

Teniendo la ecuacin de momentos, integramos una ves para

obtener la ecuacin del ngulo de deflexin.

Rax 2 T 1 sen 60 xlc 1 T 2 sen 60 xlc 2 T 3 sen 60 x lc 3 Ra x10

( x ) =Max+
+
+
+

2
2
2
2
2
2

Ahora volvemos a integrar la ecuacin para obtener la

educacin de deflexin.
Ma x 2 Rax 3 T 1 sen 60 xlc 1 T 2 sen 60 xlc2 T 3 sen 60 xlc 3 Ra x 10
+
+
+
+

2
6
6
6
6
6
3

V (x )=

Las mismas ecuaciones las planteamos en el programa de

Matlab y nos dan los siguientes resultados:

Mmax= 1.44e3 Nm

Angulo Deflexin max= .0048

Deflexin Max= .0128m

Continuamos con la ultima parte que es sacar la energa total

que el larguero y los cables absorben:
L

M (x)
Ular=
dx
0 2 EI
Ucab=

T 2Lcable
2EA

El programa de Matlab es relativamente sencillo y es lo

siguiente:

M2=Mx.^2;
intM2=trapz(x,M2);
Ular=1/(2*E*I)*intM2;
Uc1=(((T1)^2)*12)/(2*Ec*Ac);
Uc2=(((T2)^2)*18)/(2*Ec*Ac);
Uc3=(((T3)^2)*24)/(2*Ec*Ac);

Ular=40.2285 Joules
Uc1= 15.4527 Joules
Uc2= 38.0892 Joules
Uc3= 2.2411 Joules

Por ultimo, para esta primer parte, tenemos que sumar toda la
energa de cada elemento y de esta manera igualar a la energa
potencial, para despus sacar la W.
Ut=Utrave +2Ularg+ Ucables
Ut=391.608 J
Teniendo la Energa Total que puede absorber el sistema,
pasamos a calcular nuestra W, la cual se calcula dividiendo el
peso inicial que utilizamos, entre la cantidad de energa que el
sistema completo absorbi. El resultado que nos de, lo
multiplicamos por la Energa Potencial de nuestra caja, y
obtenemos nuestra nueva W.
4905
(5009.81.5)=30718.3 N
391.608

Ya que obtenemos el valor de nuestra W, procedemos a realizar

el mismo anlisis que se elaboro anteriormente, pero ahora
tomando el valor de W para nuestro planteamiento.

Para este paso solamente se pegara el cdigo y se explicaran

los resultados obtenidos.

Travesao 1.2.
clearall

E=210e9;
L=5;
h=0.204;
h1=.5;
W=30718.3;
a=1;
b=2;
I=10*10^6;

x=(0:.01:L);

Ra=(W*a)/2;
Rb=(W*a)/2;
Ma=(W*a)*(2.5)+Rb;

Mx=Ma+(Ra*x)(W/2)*(xb).^2.*heaviside(xb)+(W/2)*(xa
b).^2.*heaviside(xab);
T=(Ma*x+Ra*x.^2/2W/6*(xb).^3.*heaviside(xb)+W/6*(xa
b).^3.*heaviside(xab))/(E*I);
V=(Ma*x.^2/2+Ra*x.^3/6W/24*(xb).^4.*heaviside(xb)+W/24*(xa
b).^4.*heaviside(xab))/(E*I);

Mmax=max(abs(Mx));
Tmax=max(abs(T));
Vmax=max(abs(V));

sigmax=((Mmax*h)/2)/I;
Fs=(250e6)/sigmax;

M2=Mx.^2;
intM2=trapz(x,M2);
U=1/(2*E*I)*intM2;

Lo que nos interese en este caso es el factor de seguridad, ya

que si esta por debajo de 1 absoluto, esto significa que nuestro
travesao se rompera.

FStrave= .2553
LO CUAL NOS INDICA QUE NUESTRO TRAVESAO SE ROMPE,
CON LA FUERZA DE LA W.

De igual forma obtenemos la energa absorbida por el

travesao:

M2=Mx.^2;
intM2=trapz(x,M2);
U=1/(2*E*I)*intM2;

Utrave= 15.6558e3 J

Larguero 1.2.
clearall

Ac=.0001963;
E=210e9;
Ec=130e9;
L=20;
l=10;
x1=4;
x2=8;
x3=10;
x4=14;

lc1=6;
lc2=9;
lc3=12;
th3=60;
th2=60;
th1=60;
h=0.3048;
I=0.000010338239;
Ra=1.5359e4;
x=(0:.01:L);

A=[01sind(th1)sind(th2)sind(th3)10
10lc1*sind(th1)lc2*sind(th2)lc3*sind(th3)L1
LL^2/2(((Llc1)^2)/2)*sind(th1)(((Llc2)^2)/2)*sind(th2)(((L
lc3)^2)/2)*sind(th3)00
L^2/2L^3/6(((Llc1)^3)/6)*sind(th1)(((Llc2)^3)/6)*sind(th2)
(((Llc3)^3)/6)*sind(th3)00
lc1^2/2lc1^3/6(sind(th1)*lc1)/(Ec*Ac)0000
lc2^2/2lc2^3/6(((lc2lc1)^3)/6)*sind(th1)(sind(th2)*lc2)/
(Ec*Ac)000
lc3^2/2lc3^3/6(((lc3lc1)^3)/6)*sind(th1)(((lc3
lc2)^3)/6)*sind(th2)(sind(th3)*lc3)/(Ec*Ac)00];

B=[Ra;(l*Ra);(Ra/2*(Ll)^2);(Ra/6*(Ll)^3);(Ra/6*(lc1l)^3);
(Ra/6*(lc2l)^3);(Ra/6*(lc3l)^3)];

res=inv(A)*B;
M2=res(1);
R2=res(2);
T1=res(3);
T2=res(4);
T3=res(5);
R3=res(6);
M3=res(7);
M=M2+R2*x+T1*sind(th1)*(xlc1).*heaviside(xlc1)+T2*sind(th2)*(x
lc2).*heaviside(xlc2)+T3*sind(th3)*(xlc3).*heaviside(xlc3)Ra*(x
l).*heaviside(xl);
A=(M2*x+R2*x.^2/2+T1/2*sind(th1)*(xlc1).^2.*heaviside(x
lc1)+T2/2*sind(th2)*(xlc2).^2.*heaviside(xlc2)+T3/2*sind(th3)*(x
lc3).^2.*heaviside(xlc3)Ra/2*(xl).^2.*heaviside(xl))/(E*I);
V=(M2*x.^2/2+R2*x.^3/6+T1/6*sind(th1)*(xlc1).^3.*heaviside(x
lc1)+T2/6*sind(th2)*(xlc2).^3.*heaviside(xlc2)+T3/6*sind(th3)*(x
lc3).^3.*heaviside(xlc3)Ra/6*(xl).^3.*heaviside(xl))/(E*I);

Mmax=max(abs(M));
Amax=max(abs(A));
Vmax=max(abs(V));
sigmax=((Mmax*h)/2)/I;
Fs=(250e6)/sigmax;

Fsc=(29.62e3)/abs(T1);
Fsc=(29.62e3)/abs(T2);
Fsc=(29.62e3)/abs(T3);

M2=M.^2;
intM2=trapz(x,M2);

U=1/(2*E*I)*intM2;

Uc1=(((T1)^2)*12)/(2*Ec*Ac);
Uc2=(((T2)^2)*18)/(2*Ec*Ac);
Uc3=(((T3)^2)*24)/(2*Ec*Ac);

El factor de seguridad del larguero y cables da lo siguiente:

Fs=(250e6)/sigmax;
Fsc=(29.62e3)/abs(T1);
Fsc=(29.62e3)/abs(T2);
Fsc=(29.62e3)/abs(T3);

FSlarg= .3840
FScables= 2.1165
LO CUAL NOS INDICA QUE EL LARGUERO SE ROMPE Y LOS
CABLES NO SE ROMOPEN.

La energa absorbida por el larguero y los cables fue la

siguiente:

M2=M.^2;
intM2=trapz(x,M2);
U=1/(2*E*I)*intM2;
Uc1=(((T1)^2)*12)/(2*Ec*Ac);
Uc2=(((T2)^2)*18)/(2*Ec*Ac);
Uc3=(((T3)^2)*24)/(2*Ec*Ac);

Ularg= 3155.6 J
Ucable1= 1212.11 J
Ucable2= 2987.8 J
Ucable3= 175.794 J
Con todos los valores de energa, volvemos a realizar la
sumatoria correspondiente de acuerdo a la cantidad de
elementos que se tiene. Esta sumatoria nos tiene que dar un
valor de 30718.3N.
15.6558e3+2 ( 3155.6 )+ 2 ( 1212.11 )+ 2 ( 2987.8 ) +2 ( 175.794 ) =30718.4 J

Despus del anlisis realizado con la W que sacamos, damos

por terminado nuestro desarrollo, ya que el puente no soporta
la W y con esta fuerza ejercida ya se sobre pasa el esfuerzo
mximo permitido. Por lo cual no tendremos que tomar en
cuenta el peso y las reacciones generadas por el camin.