Puente Atirantado PDF

ESCUELA POLITÉCNICA
NACIONAL
FACULTAD DE INGENIERÍA CIVIL Y AMBIENTAL
ANÁLISIS MATRICIAL
ANÁLISIS MATRICIAL
EJERCICIO 7
PUENTE ATIRANTADO
Nombre:
Jhonny Vaca Pillajo
2020-02-17
PUENTE ATIRANTADO
Datos generales
%Puente Atirantado de Ejemplo 13-02-2020
%3 Vigas de L=10m Viga Central 30°
%2 Tirantes alineados con la viga Central y H=7m
% Viga T de 1.20x0.4x(.6+.2)
% Hormigon de f'c=280 kg/cm2
E=135000*280^0.5
miu=0.2
G=E/2/(1+miu)
A=.48
Ac=.32
Ix=0.013152
Iy=0.032
Iz=0.0272
Ea=2.04e7
Para el desarrollo del tanque atirantado necesitamos de una función denominada R_ME que nos
permitirá tener rigidices de los elementos en coordenadas reducidas de la estructura, para la
rigidez del tensor solo multiplicamos el área por el módulo de elasticidad y dividimos para la
longitud del cable.
%Rigideces de elementos en coord locales reducidas
%[KEME] = R_ME(A,Iz,Iy,Ix,L,E,G)
L=10
Lt=sqrt(10^2+5.7735^2)
[k]=R_ME(A,Iz,Iy,Ix,L,E,G)
Ra=(pi*.032^2/4*Ea/Lt)
Ingresamos las cargas en los elementos, para este caso es una carga distribuida solo en la sección
del puente y es de 2 ton/m2, hay que tomar en cuenta que en este caso el elemento está en
coordenadas locales.
%Acciones de empotramiento en coord locales completas
aep=[0 2*10/2 0 0 0 2*10^2/12 0 2*10/2 0 0 0 -2*10^2/12]'
Generamos nuestra matriz de transformación de acuerdo a las coordenadas establecidas en el

ejercicio, lo realizamos para el tensor y para la sección del puente.
%Matrices de transfrmacion MEspacial
%[TME,h,Taep]=MT_ME(p1,p2,p1p,p2p,Fi,bet1,bet2)
[T1,h,Ta1]=MT_ME([0 0 0],[10 0 0],[0 0 0],[10 0 0],0,0,0),T1, Ta1
[T2,h,Ta2]=MT_ME([10 0 0],[18.6602 5 0],[10 0 0],[18.6602 5 0],0,0,0),T2,
Ta2
[T3,h,Ta3]=MT_ME([18.6602 5 0],[28.6602 5 0],[18.6602 5 0],[28.6602 5
0],0,0,0),T3, Ta3
%Matrices de transfrmacion Arm Espacial

%[T]= MT_AE(p1,p2, beta1,beta2)
[T4]= MT_AE([0 -5.7735 7],[10 0 0],0,0)
[T5]= MT_AE([28.6602 10.7735 7],[18.6602 5 0],0,0)
Ya obtenido las matrices de transformación para cada elemento generado obtenemos la matriz
de rigidez de cada elemento en coordenadas de nudo, donde se genera mediante la siguiente
relación:
%Rigideces en coordenadas de nudo
K1=T1*k*T1'
K2=T2*k*T2'
K3=T3*k*T3'
K4=T4*Ra*T4'
K5=T5*Ra*T5'
Igual para este ejercicio necesitamos encontrar las acciones de empotramiento en coordenadas y
para ello lo encontramos con la matriz de transformación multiplicado por la carga generada
anteriormente en toda la sección.
%Acciones de empotramiento de MEsp en coord de nudos
AEP1=Ta1*aep
AEP2=Ta2*aep
AEP3=Ta3*aep
Ingresamos los resultados de los vectores de ensamble en la codificación principal y generamos

nuestro ensamble de la matriz de rigidez de toda la estructura, en donde vamos a introducir el
ensamble de la matriz de rigidez para cada elemento de la estructura, para ellos necesitamos los
grados de libertad de la estructura, los vectores de ensamble y las rigideces de los elementos
%Vectores de ensamble
lee1=[19:21 1:9];
lee2=[4:15];
lee3=[10:15 22:24 16:18];
lee4=[25:27 4:6]
lee5=[28:30 10:12]
leeME=[lee1;lee2;lee3]
leeAE=[lee4;lee5]
%Emsamble de la Matriz de Rigidez de la estructura
%[S,B]=Ens2(lee, K, ngdl, S, B)
S=zeros(18), B=zeros(30,1);
[S,B]=Ens2(lee1, K1, 18, S, B);

[S,B]=Ens2(lee2, K2, 18, S, B);
[S,B]=Ens2(lee3, K3, 18, S, B);
[S,B]=Ens2(lee4, K4, 18, S, B);
[S,B]=Ens2(lee5, K5, 18, S, B);
Generamos un ensamble de la acumulación de las cargas con los vectores generados, para este
caso no se tiene cargas puntuales en B por lo que se hace el mismo procedimiento que el ejercicio
de la armadura.
%%Ensamble SAEP :[PR]=Acum2(lee, Ac, PR)
[B]=Acum2(lee1, AEP1, B)
[B]=Acum2(lee2, AEP2, B)
[B]=Acum2(lee3, AEP3, B);P=-B(1:18)
Resolvemos el sistema de ecuaciones esto con la finalidad de calcular los desplazamientos en los
grados de libertad existentes. Así verificamos la compatibilidad de los nudos para que exista la
continuidad.
%Solucion del sistema de ecuaciones
Delta=S\P; Delta(30)=0
Ya generada la matriz de desplazamientos Delta, se extrae los desplazamientos para cada

elemento en coordenadas de nudo en el extremo inicial y final. Esto se puede hacer mediante el
vector generado lee
%desplazamientos en cada elem en coord de nudos
D1=Delta(lee1);
D2=Delta(lee2);
D3=Delta(lee3);
DME=[D1 D2 D3]
D4=Delta(lee4);
D5=Delta(lee5);
DAE=[D4 D5]
Obtenemos el desplazamiento y los giros generados de un elemento en coordenadas reducidas,

esto si considerando que tienen cargas en los elementos de la dirección de gravedad. Entonces
utilizaos la matriz de transformación
% desplazamientos en cada elemento en coord locales reducidas
d1=T1'*D1;
d2=T2'*D2;
d3=T3'*D3;
dME=[d1 d2 d3]
d4=T4'*D4;
d5=T5'*D5;
dMA=[d4 d5]
A continuación, se genera los desplazamientos en acciones, mediante la siguiente multiplicación.

a1=k*d1;
a2=k*d2;
a3=k*d3;
aME=[a1 a2 a3]
a4=Ra*d4;
a5=Ra*d5;
aMA=[a4 a5]
Generamos las acciones de los elementos en las coordenadas locales completas multiplicando las
acciones locales reducidas por la altura y sumando las cargas establecidas.
%Acciones en coordenadas locales completas
ac1=h*a1+aep;
ac2=h*a2+aep;
ac3=h*a3+aep;
acME=[ac1 ac2 ac3]
las acciones en los elementos del sistema completo son calculadas con la matriz de equilibrio ya
que considera los efectos generados por cargas en los nudos iniciales y finales de cada elemento
y las acciones en coordenadas reducidas.
%Acciones de elem en coord de Nudos
A1n=Ta1*ac1;
A2n=Ta2*ac2;
A3n=Ta3*ac3;
AnME=[A1n A2n A3n]
A4n=T4*a4;
A5n=T5*a5;
AnMA=[A4n A5n]
%Obtener reaaciones y acciones en gdl
PR=zeros(30,1);
%[SAEP]=Ens_SAEP(lee,AEP,SAEP,ngdl)
[PR]=Ens_SAEP(lee1,A1n,PR,30);
[PR]=Ens_SAEP(lee5,A5n,PR,30);PR
Resultados Matlab
Matriz de ensamble S
Desplazamientos Cargas en grados de libertad y reacciones

Resultados SAP2000
Momento
Axial
Reacciones
El SAP2000 es una herramienta complementaria que puede ayudar a las personas a verificar sus
resultados provenientes de una programación realizada por él mismo en este caso, y corroborar
estos resultados, así como cuestionar la veracidad de los mismos, si se sabe utilizar de manera
correcta.
Los resultados obtenidos de la programación realizada en MATLAB, se diferencian en cierto grado

de los resultados obtenidos en el SAP2000, esto puede ser debido a que entre las propiedades de
estos materiales difieren sustancialmente, el ingreso de un dato o la ausencia del mismo, es lo
que podría generar esta diferencia.

Puente Atirantado PDF

Cargado por

Copyright:

Formatos disponibles

Puente Atirantado PDF

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Puente Atirantado PDF

Cargado por

Copyright:

Formatos disponibles

ESCUELA POLITÉCNICA

Generamos nuestra matriz de transformación de acuerdo a las coordenadas establecidas en el

%Matrices de transfrmacion Arm Espacial

Ingresamos los resultados de los vectores de ensamble en la codificación principal y generamos

[S,B]=Ens2(lee1, K1, 18, S, B);

Ya generada la matriz de desplazamientos Delta, se extrae los desplazamientos para cada

Obtenemos el desplazamiento y los giros generados de un elemento en coordenadas reducidas,

A continuación, se genera los desplazamientos en acciones, mediante la siguiente multiplicación.

Desplazamientos Cargas en grados de libertad y reacciones

Los resultados obtenidos de la programación realizada en MATLAB, se diferencian en cierto grado

También podría gustarte