% UNSCH - FIMGC - EFPIC % PROGRAMA PARA EL ANALISIS ESTRUCTURAL DE PORTICOS USANDO LA MATRIZ A disp('==================================================================' disp('==================================================================' ) disp('==================================================================' disp('==================================================================' ) fpri!f ('" ('" ###################### PORTICOS $$$$$$$$$$$$$$$$$ "' ) fpri!f ('" ('" ######################## POR METODO DE RIGIDECES $$$$$$$$$$$$$$$$$$$ "') "' ) disp('==================================================================' disp('==================================================================' ) fpri!f ('" ('" ######################## ENTRADA DE DATOS $$$$$$$$$$$$$$$$$$$$$$ "' ) disp('==================================================================' disp('==================================================================' ) % ARREGLO CG& COORDENADAS GENERALES d=ip!('" d=ip!( '" DE NUDOS* ' ) r=ip!('" r=ip!( '" DE NUDOS RESTRINGIDOS*' ) +,= CG=.s(d/0) % ANALISIS DE RESTRICCIONES fr i=1*r fr i=1*r dr.s= ip! ( '" DEL NUDO RESTRINGIDO*' ) 21 = ip! ( ' " DESPLAZAMIENTO . 2 /si(s) ()* ' /'s' 's') ) if 21== if 21=='' '' CG(dr.s/1)= .,s./ .,s. /.d 31 = ip! ( '" DESPLAZAMIENTO . 3 /si(s) ()* ' /'s' 's') ) if 31== if 31=='' '' CG(dr.s/4)= .,s./ .,s. /.d R1 = ip! ( '" ROTACION /si(s) ()*' /'s' 's') ) if R1== if R1=='' '' CG(dr.s/0)= .,s./ .,s. /.d .d % GRADOS DE LI5ERTAD fr i=1*d fr i=1*d fr 6=1*0 fr 6=1*0 if CG(i/6)7= if CG(i/6)7= +,=+,81 CG(i/6)=+, .,s./ .,s. /.d .d .d CG % ARREGLOS 9ECTORIALES ii : fi (NUDOS INICIALES 3 FINALES) ;r ii ,.>r fi fr i=1*;
ii fi % ARREGLO 9C& 9ECTORES DE COLOCACION ,.>r 9C fr i=1*;r L ,.>r SENO ,.>r COSENO fr i=1*;r A fr i=1*;
.,s.if 9C(i/p)==CG(;/) p#0 ==4 i(i)=f(i)=SENO(i) i(i)=f(i)=COSENO(i) +i(i)=+f(i)= .,s.if 9C(i/p)==CG(;/) p#0 ==0 i(i)=f(i)= i(i)=f(i)= +i(i)=+f(i)=1 .,s./.d .d ,=CG(;/) fi,=(i-1)0 A(fi,81/,)=+i(i)-(f(i)-i(i))L(i) A(fi,84/,)=+f(i)-(f(i)-i(i))L(i) A(fi,80/,)=f(i)-i(i) .,s./.d
.d .d
.d fpri!f ('" MATRIZ A (MATRIZ DE COMPATI5ILIDAD)* "' ) A % CALCULO DE LA MATRIZ DE RIGIDEZ DE LA ESTRUCTURA J fpri!f ('" CARACTERISTICAS DE LOS MIEM5ROS*' ) E,>s= ip! ('" MODULO DE ELASTICIDAD T;4B*' ) G=&KE,>s 5.!>= ip! ('" FACTOR DE FORMA*' ) ,.>r J4I J=.rs(+,/+,) fr i=1*;(i) = 5(i)H(i) I.ri>(i) = 5(i)H(i)014 fi = (0E,>sI.ri>(i)5.!>)(GAr.>(i)L(i)4) EIK=((KE,>sI.ri>(i))L(i))((18fi)(18Kfi)) EI4=((4E,>sI.ri>(i))L(i))((1-4fi)(18Kfi)) EA=((E,>sAr.>(i))L(i)) % MATRIZ DE RIGIDEZ DE MIEM5RO J4 EN COORD& LOCALES J4(i/1/1)=EIKJ4(i/1/4)=EI4J4(i/1/0)= J4(i/4/1)=EI4J4(i/4/4)=EIKJ4(i/4/0)= J4(i/0/1)=J4(i/0/4)=J4(i/0/0)=EA % MATRIZ DE TRANSFORMACION T40(i/1/1)=COSENO(i)T40(i/1/4)=SENO(i)T40(i/1/0)=T40(i/1/K)=T4 0(i/1/)=T40(i/1/)= T40(i/4/1)=SENO(i)T40(i/4/4)=COSENO(i)T40(i/4/0)=T40(i/4/K)=T40(i/4/)=T 40(i/4/)= T40(i/0/1)=T40(i/0/4)=T40(i/0/0)=1T40(i/0/K)=T40(i/0/)=T4 0(i/0/)= T40(i/K/1)=T40(i/K/4)=T40(i/K/0)=T40(i/K/K)=COSENO(i)T40(i/K/ )=SENO(i)T40(i/K/)=
T40(i//1)=T40(i//4)=T40(i//0)=T40(i//K)=SENO(i)T40(i//)=COSENO(i)T40(i//)= T40(i//1)=T40(i//4)=T40(i//0)=T40(i//K)=T40(i//)=T4 0(i//)=1 % MATRIZ DE RIGIDEZ DEL MIEM5RO ? EN COORD& LOCALES fr ;=1*0 fr =1*0 J4I(;/)=J4(i/;/) .d .d fpri!f ('" MATRIZ DE RIGIDEZ ELEMENTO EN COORD& LOCALES*' ) ? = J4I fi, =(i-1)0 Ap=A(fi,81*fi,80/1*+,) Jp=Ap'?Ap J=J8Jp .d fpri!f ('" MATRIZ DE RIGIDEZ DE LA ESTRUCTURA*' ) J PORTICO4
Como en el caso del pórtico en 2d se proceden a ingresar los datos de la siguiente manera: INGRESO DE DATOS AL PROGRAMA
1) El número de nudos 2) El número de nudos con apoyo 2.1) Las restricciones en cada apoyo en el que n os piden para confrmar si hay desplazamiento en los respectios e!es la cual se de"e confrmar con un si#s) o un no#n) $) %úmero de miem"ros $.1) en cada miem"ro designamos la conectiidad de cada miem"ro designando el nudo inicial y el nudo fnal &) 'ngresamos las coordenadas de los nudos en metros teniendo como re(erencia la siguiente conención de signos
) Vj Y Uj y
Vi
x
Ui X
*) +espu,s ingresamos las caracter-sticas del material como elasticidad y rea .1) ingresamos la elasticidad del material en /0m 2 y las dimensiones de la sección. ) 'ngresamos las cargas en los miem"ros de acuerdo a como sean nuestras cargas. .1) i son cargas puntuales descomponemos en sus reacciones nodales empotrados per(ectamente e ingresamos esos datos. .2)i es una carga distri"uida ingresamos el alor de la carga con su respectio sigo de conención según este programa realizado.
SALIDA DE DATOS DEL PROGRAMA
1) Los nudos inicial y fnal ordenados 2) Las longitudes de cada elemento a partir de los nudos $) Cosenos directores de cada elemento &) 3atriz 4 #matriz de compati"ilidad ) ) 3atriz de rigidez elemento en coord. Locales *) 3atriz +e 5igidez +e La Estructura ) %os ensam"la las cargas nodales con la de las cargas en los apoyos 6) Luego hace la operación de cargas totales y con esto calcula 7) 8ector de desplazamientos 19) 8ector de de(ormaciones 11) Cargas internas 12) 4cciones en coordenadas locales #. complementario) 1$) 4cciones en coordenadas locales #. primario)
rosigamos con el e!ercicio propuesto