SIMULADOR DE ESTRUCTURAS RETICULARES Autores Bohrquez vila C.A1 Bravo Alvarado c.bohor@gmeail.

com

RESUMEN El mtodo del elemento finito ha llegado a ser una herramienta poderosa en la solucin numrica de un amplio rango de problemas de ingeniera, las aplicaciones van desde el anlisis por deformacin y esfuerzo en aplicaciones como: automviles, aeronaves, edificios, y estructuras de puentes hasta el anlisis de campos de flujo de calor, de fluidos, filtraciones y otro tipo de problemas. Con los avances en la tecnologa de las computadoras y otros sistemas CAD, pueden modelarse problemas complejos con relativa facilidad. En una computadora pueden probarse varias configuraciones alternas antes de construir el primer prototipo. Estas caractersticas facilitan la generacin de programas en computador para solucionar problemas que normalmente podran demorarse horas o hasta das, en este caso se desarrolla un programa sobre la plataforma de visual Basic para el clculo de esfuerzos en estructuras en 2D y 3D.

codigo 1335

CELOCIAS En la Figura 1 se muestra una armadura plana tpica. Con sus grados de libertad por nodo.

Figura 1 Armadura en 2D tpica. Con sus grados de libertad Una armadura estructural slo consta de miembros sujetos a dos fuerzas, es decir, cada elemento de la armadura est en tensin o compresin directa. En una armadura se requiere que todas las cargas y reacciones estn aplicadas slo en los nodos y que todos los miembros estn conectados entre s en sus extremos por medio de articulaciones sin friccin, con el mtodo del elemento finito hay ventajas como: Se pueden hallar los desplazamientos en los nodos fcilmente. El mtodo se aplica tanto en armaduras estticamente determinadas como en indeterminadas. Anlisis en dos dimensiones de armaduras y en tres dimensiones con tan solo algunos cambios. Debido al planteamiento matricial del mtodo este es adaptable a la programacin en PC.

ANALISIS DE ARMADURAS EN DOS DIMENDIONES

Figura 2 Sistema de coordenadas locales y globales Los elementos con los cuales esta formada una armadura tienen varias orientaciones, para tomar en cuenta estas orientaciones, se manejan dos sistemas de coordenadas, uno denominado coordenadas globales y otro coordenadas locales como muestra la Figura 2. En esta figura se observa el esquema de numeracin local, en donde los dos nodos del elemento se numeran 1 y 2. el sistema local de coordenadas cosiste en el eje X que va en la misma direccin del elemento desde el nodo 1 al 2; en el sistema local se denotan todas las variables con ().

En el sistema global de coordenadas los ejes X y Y son fijos y no dependen de la ubicacin del elemento, hay que resaltar que el plano XY es perpenticular a la hoja. Para el sistema global cada nodo tiene dos grados de libertad (gdl); asi cada nodo quedar con los gdl 2j-1 y 2 y tendr asociado los desplazamientos del nodo j Q2j-1 y Q2j como se muestra en la figura 1. En el sistema de coordenadas local q1 y q2 son los deplazamientos de los nodos 1 y 2 repectivamente, el vector de desplezamiento en coordenadas locales sera: q =[q1 q2]T y en el sistema de coordenadas global sera : q =[q1 q2 q3 q4]T (2) De la figura 3 la relacion entre q y q se da asi : q1 es igual a la suma de las proyecciones de q1 y q2 sobre X q'1=q1 cos + q2 sen igualmente q'2=q3 cos + q4 sen (4) de lo anterior se puede introducir los cosenos directores l y m, los cosenos directores son los cosenos de los ngulos que el eje local X forma con los ejes globales X y Y las ecuaciones 3 y 4 se pueden escribir en forma matricial asi: q=Lq Y la matriz L de transformacin esta dada por: (5) (3) (1)

l m 0 0 L= 0 0 l m
CALCULO de l y m Para tal fin se debe tener en cuenta la figura 4 y las coordenadas nodos 1 y 2 .

(6)

de

los

l=

x x l
2 e

m=

y y l
2 e

(7)

Y la longitud le se obtiene con :

(x 2 x1) + (y 2 y1)
2

(8)

Figura 3 Cosenos Directores

Estos son los datos de cada elemento que se tienen en cuenta en el simulador para determinar la ubicacin de los elementos con respecto a los ejes X y Y. MATRIZ DE RIGIDEZ DE UN ELEMENTO Para un elemento considerado en el sistema de coordenadas locales se considera que es unidimencional, y la teoria del elemento finito da la rigidez del elemento unidimencional como sigue: K=

EA l
e e

1 1 1 1

(9)

Donde Ae es el area de seccion transversal del elemento y Ee es el modulo de young. Ahora para obtener una expresin para la matriz de rigidez del elemento en el sistema coordenado global se tiene en cuenta la energa de deformacin unitaria en el elemento; la energa de deformacin unitaria del elemento en coordenadas locales es:

Ue =

1 T q kq 2

(10)

Sutituyendo q=Lq en (10) se obtiene

Ue =

1 T T q [L kL]q 2
1 T q kq 2

(11)

La energa de deformacin unitaria en coordenadas globales es: Ue = (12)

donde k es la matriz de rigidez en coordenadas globales, de lo anterior obtenemos dicha matriz como k=LTkL sustituyendo L de (6) y k (9) se obtiene
l2 lm lm m2 K = E e Ae 2 l e l lm 2 lm m l 2 lm lm m 2 l2 lm lm m2

(13)

(14)

CALCULO DE ESFUERZOS Para los esfuerzos en el elemento pueden obtenerse expresiones considerando que un elemento armadura en coordenadas locales es un miembro simple con dos fuerzas figura 1, as el esfuerzo en un elemento armadura esta dado por

= Ee
Como la deformacin unitaria es el cambio de logitud por unidad de longitud original,

(15)

= Ee

q' q ' E = [ 1 1] q l l
2 1

(16)

la ecuacion anterior puede escribirse en terminos de los desplazamientos globales q usando la transformacion q=Lq

=
La sustitucin de L de la ecuacin (6) da

E [ 11]L q l
e e

(17)

E [ l l
e e

m l

m] q

(18)

Una vez determinados los desplazamientos al resolver las ecuaciones de los elementos finitos, los esfuerzos pueden calcularse de la ecuacin (18) para cada elemento. Hay que notar que un esfuerzo positivo implica que el elemento esta en tensin y un esfuerzo negativo implica compresin. ARMADURAS TRIDIMENCIONALES El elemento armadura tridimencional puede manejarse como una generalizacin directa del elemento armadura bidimencional que se meciono antes, los sistemas de coordenadas local y global para un elemento armadura tridimencional se muestran en la figura 4, hay que notar que el sistema de coordenadas local es de

nuevo el eje X y esta alineado a lo largo del elemento, ya que un elemento armadura es sismplemente un miembro con dos fuerzas. En consecuencia el vector de desplazamiento nodal en coordenadas locales es q =[q1 q2]T el vector de desplazmineto nodal en coordenadas globales es ahora q =[q1 q2 q3 q4 q5 q6]T con referencia a la figura 4 encontrmos que la transformacin entre las coordenadas local y global es q=Lq donde la matriz de transformacion L esta dada por (21) (20) (19)

l m n 0 0 0 L= 0 0 0 l m n

(22)

Donde l, m y n son los cosenos derectores del eje X local con respecto a los ejes x, y , z globales, respectivamente. La matriz de rigidez del elemento en coordenadas globales esta dada por la ecuacin (13) y resulta
l2 l 2 lm nl lm nl lm m2 mn lm m 2 mn 2 2 nl mn n mn n EA nl K= L l 2 lm nl l2 lm nl lm m 2 mn lm m2 mn 2 nl mn n2 nl mn n

(23)

Las

l=

formulas x2 x1

calcular l m y n y y1 n = Z 2 Z 1 (24) m= 2

para

son

La longitud le es

(x 2 x1) + (y 2 y1) + (Z 2 Z 1)
2 2

(25)

Figura 4 elemento armadura tridimensional en coordenadas local y global VIGAS Las vigas son miembros esbeltos que se usan para soportar cargas transversales. Los miembros horizontales largos usados en edificios y puentes, y las flechas apoyadas en cojinetes son algunos ejemplos de vigas. Las vigas con secciones transversales que son simtricas respecto al plano de carga sern consideradas en este artculo. En la figura 5 se muestra una viga horizontal comn, muestra tambin la seccin transversal y la distribucin del esfuerzo por flexin. Para deflexiones pequeas, recordamos de a teora elemental de vigas la siguiente ecuacin:

= =
2

E
2

M y I

(26)

(27) (28)

d v=M dx EI

Figura 5 Viga comn con cargas Donde es el esfuerzo normal, es la deformacin unitaria normal, M es el momento flexionante en al seccin, v es la deflexin del eje centroidal en la posicin X e I es el momento de inercia de la seccin respecto al eje neutro. FORMULACIN DEL ELEMENTO FINITO Las vigas se dividen en elementos, como se muestra en la figura 6, Cada nodo tiene dos grados de libertad. Tipicamente los grados de libertad del nodo i son Q2i-1 y Q2i . el grado de libertad Q2i-1 es un desplazamiento transversal y Q2i es una pendiente o una rotacin. El vector:

Q=[Q1, Q2,..,Q10]T

(32)

Representa el vector de desplazamientos globales. Para un solo elemento, los grados de libertad locales estan representados por: q=[q1, q2, q3, q4]T (33) la correspondencia local-global es facil de ver en la tabla de la figura 6. q es lo mismo que [v1,v1,v2,v2]T las funciones de forma para interpolar v en un elemento estan definidas en terminos de de -1 a , Las funciones de forma para elementos de viga difieren de los analizados en armaduras, como estan implicados valores nodales y pendientes nodales, se definen las funciones de forma de Hermite, que satisface valores nodales y requisitos de continuidad de pendiente. Cada una de las funciones de forma es de orden cbico representadas por: Hi=ai+bi+ci2+di3 i=1,2,3,4 (34) Figura 6 funcin de forma de hermite H2 1 0 H3 0 1 H3 0 0 H4 0 0 H4 0 1 H=[H1, l e H2,H3, l e H4]
2

Se debe satisfacer las condiciones que se dan en la siguiente tabla. H1 =-1 =1 1 0 H1 0 0 H2 0 0

De donde se halla el vector de carga H: (35)

Para hallar la matriz de rigidez del elemento ke es

(36)

FUERZA CORTANTE Y MOMENTO FLEXIONANTE Usando las ecuaciones de fuerza cortante y momento flexionante
M = EI

dv dx
2

V =

dM dx

v = Hq

Se obtiene el momento flexionante y la fuerza cortante del elemento: M=

EI [6q1+(3-1)leq2-6q3+(3+1)leq4] l e2
6 EI (2q1+ leq2-2q3+ leq4) l e3

(37)

V =

(38)

Los valores anteriores del momento flexionante y la fuerza cortante son para las cargas modeladas usando cargas puntuales equivalentes. Asumiendo las cargas del equilibrio en los extremos del elemento como R1, R2, R3, R4, se deduce que:
R1 R 2 = EI R3 l e3 R 41 12 6l e 12 6 le 6 le 4l 2 e 6 le 2l 2 e 12 6 le 12 6 le pl e 2 6 l e q1 pl 2 e 2 l 2 e q 2 12 + 6 l e q 3 pl e 2 4 l e q 4 122 pl e 12 e

(39)

Se puede observar que el primer termino de la derecha es k q. hay que tener en cuenta que el segundo termino de la derecha es el vector de carga y se debe tener en cuenta solo en elementos con carga distribuida; en libros de anlisis matricial estructural, las ecuaciones anteriores se describen directamente del equilibrio del elemento, tambien el ltimo vector al lado derecho de (39) consta de terminos que se llaman reaccioones de extremo fijo. Las fuerzas cortantes en los extremos del elemento son: V1=R1 y V2=-R3, los momentos flexionantes en los extremos son M1=-R2 y M2=R4. MARCOS PLANOS Aqu se consideran las estructuras planas con miembros conectados rgidamente. Esos miembros sern similares a las vigas excepto que se tendrn presentes cargas axiales y deformaciones axiales. Los elementos tambin tendrn diferentes orientaciones. La figura 7 muestra el elemento de un marco. En cada nodo se tiene dos desplazamientos y una deformacin rotacional. El vector de desplazamiento nodal esta dado por:

q= [q1,q2,q3,q4,q5,q6]T (40) Tambien se define un sistema coordenado locar o del cuerpo x, y tal que x este orientada a lo largo de 1-2, con cosenos directores l, m (donde l=cos , m=sen). Estos se evaluan usando las relaciones dadas para el elemento armadura, mostrado en la figura 3. el vector de desplazamiento nodales en el sistema local es:

q=[q1,q2,q3,q4,q5,q6]T

(41)

FFigura 7 elemento de marco reconociendo que q3= q3 y q6= q6, que son rotaciones con respecto al cuerpo, obtenemos la transformacin local-global: q=Lq donde
m 0 0 0 0 l m l 0 0 0 0 0 0 1 0 0 0 L= 0 0 l m 0 0 0 0 0 m l 0 0 0 0 0 1 0

(42)

(43)

Se observa que q2, q3, q5, y q6] son los grados de libertad de las vigas, mientras que q1 y q4] son similares a los dezplazamientos de una barra,como se vio en paginas anteriores. Combinando las dos rigideces y situndolas en las posiciones apropiadas, obtenemos la rigidez del elemento para un elemento de marco como:
EA L e 0 0 = EA Le 0 0 0 12 EI 3 Le 6 EI 3 Le 0 12 EI 3 Le 6 EI 2 Le 0 6 EI 2 Le 4 EI Le 0 6 EI 2 Le 2 EI Le EA Le 0 0 EA Le 0 0 0 12 EI 3 Le 6 EI 2 Le m 12 EI 3 Le 6 EI 2 Le 0 6 EI 2 Le 2 EI Le 0 6 EI 2 Le 4 EI Le

(44)

K re

PORTICOS TRIDIMENSIONALES Cada nodo de prtico tridimensional se compone de seis grados de libertad, en lugar de tres que componen los nodos de los prticos en dos dimensiones un ejemplo de prtico en tres dimensiones se puede ver en la figura 11. La numeracin de los grados de libertad se realiza de la siguiente manera; para el nodo j hay tres desplazamientos y tres rotaciones, para los desplazamientos en X, Y, Z los grados de libertad son 6j5, 6j4 y 6j3, para las rotaciones 6j 2, 6j 1 y 6jrespectivamente. Los elementos se ubican en coordenadas locales y con la ayuda de la matriz de transformacin (45) se pasan a coordenadas globales, para poder realizar el ensamble de la matriz. Global del modelo con las matrices de rigidez (49) en Coordenadas globales. La matriz de transformacin es De 12 x 12 y se compone como Figura 11 Prtico Tridimensional sigue:
0 L= 0

(45)

Es la matriz de csenos directores y la matriz general es de 12 x 12

CARACTERISTICAS SIMULADOR DE ESTRUCTURAS El simulador sirve para disear los tipos de elementos que se aplican en un estructura; en esta etapa el simulador disea elementos tipo Celosa 2D, celosa 3D y Viga 2D por el mtodo de esfuerzos admisibles, pero en la siguiente etapa se disearn por el mtodo de factores de carga y resistencia. Para estar acorde con la normatividad al respecto se disean de Acuerdo a la Norma Tcnica Colombiana NSR-98. Las Normas NSR-98 estipulan que los algunos materiales pueden disearse de acuerdo con esta metodologa. Se aplican tambin los conceptos del diseo estructural basados en factores de Carga y seguridad Visual Basic fue el programa que se uso para compilar el cdigo que genera la aplicacin, es de fcil manejo y ofrece gran variedad de controles para el desarrollo de aplicaciones, la aplicacin no tiene vnculos con algn otro programa y es compatible con cualquier equipo que tenga Windows como sistema operativo. A continuacin se muestra el diagrama de flujo y algunas impresiones de las pantallas que aparecen en la ejecucin del programa

Bibliografa [1] Introduction to Finite Elements in Engineering. Tirupathi R. Chandrupatla, 2002. [2] Mecnica de materiales.,Beer y Johnston, McGRAW-HILL, Edicin 2,1993 . [3] Anlisis de estructuras.,Jairo Uribe Escamilla, Escuela Colombiana de ingeniera, Edicin 2, 2004. [4] Tentative Provisions for the Development of Seismic Regulations for Building., Palo Alto,Appied technologi council,1978. [5] Estudio general de amenaza sismica de Colombia.,Comit AIS 300,Asociacin Colombiana de ingeniera,1996. [6] Building code Requirements for reinforced concrete.,ACI 318-95 American concrete institute,1995. [7] MSDN Library ., Microsoft Corp,2000. [8] Madsen, H.O, Krenk, S. Lind,N.C. ., Pretince-Hall,1986. [9] Some Thoughts on Maximum Probable Loads and Structural Safety Insurance., Cornell C.a.,Pretince-Hall,1967.

[10] Normas Colombianas de Diseo y Construccin Sismo Resistente., Asociacin Colombiana de ingeniera ssmica (AIS),1998. [11] Safety and Limit States for Reinforced Concrete.,MacGregor J.G. ,Canadian Journal of Civil Engineering,1976. [12] Control de calidad en estructuras de acero.,Quiroga. ,Tesis para obtener el titulo de Magster en ingeniera civil, universidad de los Andes,1990. [13] Aseguramiento de la calidad en estructuras metlicas.,Aragn ,Proyecto de grado para obtener el titulo de ingeniero Civil, universidad de los Andes,1995. [14] Diseo de estructuras de acero utilizando coeficientes de carga y resistencia., Uribe J.,Sociedad Colombiana de ingenieros,1979. [15] Load and Resistance Factor Design for steel., Ravindra y Galambos,Journal of the Structural Division,1979.