ANALISIS ESTRUCTURAL DE UN RESERVORIO CIRCULAR APOYADO CON CUPULA DE CONCRETO ARMADO POR EL MTODO DE LOS ELEMENTOS FINITOS

Franz Condori Pacoricona 06-29883 ESIC-FAIN-UNJBG

PARTE I

OBJETIVOS

1. OBJETIVOS
1.1 General
Conocer y analizar estructuras formadas por slidos axisimtricos por el mtodo de los elementos finitos

1.2 Especficos
Utilizar el mtodo de los elementos finitos para el anlisis de un reservorio apoyado con cpula de concreto armado. Comprobacin de los resultados obtenidos con un software de anlisis estructural profesional (SAP2000). Programacin de una aplicacin en Python que permita estudiar los clculos de proceso que se realizan en un anlisis por el mtodo de los elementos finitos.

PARTE II

MARCO TERICO

2. ANLISIS DE ESTRUCTURAS CON ELEMENTOS FINITOS

2.1. Mtodo de los elementos finitos
2.1.1. Perspectiva histrica Segn los expertos no es posible determinar el momento exacto en que se inicia 2.1.2. Definicin Algunas definiciones dadas por algunos libros son:

O.C Zienkiewicz - R.L. Taylor(1993) es objeto de este texto presentar un panorama del mtodo de los elementos finitos como un procedimiento general de discretizacin de los problemas continuos planteadas por expresiones definidas matemticamente.

ANSYS-CivilFem(2003) el Mtodo de los Elementos Finitos, es uno de los procedimientos que existen para aproximar el comportamiento de una estructura con infinitos grados de libertad por el de otra, con aproximadamente las mismas propiedades fsicas y geomtricas, pero con un nmero finito de grados de libertad, cuyas ecuaciones de equilibrio pueden expresarse por un sistema algebraico de ecuaciones simultneas con un nmero limitado de incgnitas.

2.1.3. Esquema de presentacin En el presente trabajo se adopta el enfoque de la energa potencial para la formulacin del mtodo de los elementos finitos.

Principio de la energa potencial mnima.- Para sistemas conservativos, de todos los campos de desplazamiento cinemticamente admisibles, aquellos que corresponden a condiciones de equilibrio extremizan la energa potencial total. Si la condicin extrema es un mnimo, el estado de equilibrio es estable.

2.2. Modelamiento de elementos estructurales

Solo se comentarn brevemente los modelos para los elementos estructurales en un anlisis mediante el mtodo de los elementos finitos, debido a lo extenso del tema.

2.2.1. Elementos Adimensionales

Se consideran objetos adimensionales a los nudos, Mediante los nudos estn conectados los elementos en un discretizacin para el anlisis y representan los puntos, sobre la estructura, en las cuales se desea obtener los resultados requeridos, Figura 2.1 y 2.2.

2.2.2. Elementos unidimensionales El elemento unidimensional, la lnea, es un objeto que sirve para modelar el comportamiento de elementos estructurales como: armaduras, vigas, columnas, cables, etc. y todo elementos estructural en la que la longitud es mucho mayor que las dimensiones que definen su seccin transversal, Figura 2.1.

Figura 2.1 Para elementos estructurales unidimensionales, la formulacin del mtodo de los elementos finitos termina siendo el mtodo de la rigidez, adems de existir otros mtodos que pueden obtener resultados en menos tiempo y con gran exactitud.

2.2.3. Elementos Bidimensionales Se utiliza para modelar elementos estructurales de superficie, figura 2.2 y son: Cscaras (Shell).- Las cscaras son utilizadas para modelar estructuras bidimensionales o tridimensionales, como el comportamiento tipo cscaras de tanques cpulas o bvedas; como el comportamiento tipo placa de las losas o como el comportamiento tipo membrana de los muros de corte. El plano.- Para modelar slidos bidimensionales, en el anlisis del comportamiento de esfuerzos planos en estructuras planas y delgadas, o tambin para modelar el comportamiento de deformaciones planas en estructuras prismticas largas.

Slidos Axisimtricos.- Es usada para modelar slidos simtricos respecto a un eje, los cuales estn sujetos a cargas simtricas con respecto al mismo eje. ste es el tipo de estructuras que se analizar en el presente trabajo.

Figura 2.2

2.2.4. Elementos tridimensionales Para analizar estructuras utilizando modelos de slidos tridimensionales.

2.3. Slido Axisimtrico

Seguidamente se realiza el resumen del mtodo de los elementos finitos para estructuras que cumplen las condiciones para ser modelo como un slido de revolucin Figura 2.3. Se estudiar nicamente el anlisis de slidos con simetra axial sometidos a cargas de revolucin.

Figura 2.3 Eugenio Oate (1992).- Consideramos slidos en los que su geometra y propiedades mecnicas son independientes de la coordenada circunferencial Figura 2.3. Aunque

el comportamiento de dichos slidos es tridimensional, su estudio matemtico es generalmente bidimensional ya que en la mayora de los casos puede efectuarse utilizando variables que dependen de dos coordenadas cartesianas. Si las cargas exteriores son tambin de revolucin, el desplazamiento de un punto de una estructura considerada como un slido de revolucin tiene solo componentes en direccin radial ( ) y axial ( ).

2.3.1. Estructuras de simetra axial Las estructuras de revolucin representan un alto porcentaje en las estructuras de uso comn en ingeniera, entre ellos se puede mencionar: Los depsitos para gua, muros cilndricos, silos, techos de edificaciones especiales, etc.

2.3.2. Modelo como slido axisimtrico La Figura 2.3a muestra una estructura de revolucin general, debido a que la estructura tiene una completa simetra respecto al eje z del sistema de orientacin global, es posible modelar en el plano formado por su seccin caracterstica de revolucin Figura2.3b; en el modelo, para obtener la estructura entera, ser necesario girar una vuelta entera, 360 grados sexagesimales alrededor del eje z. Sobre esta seccin de revolucin se plantea el mtodo de los elementos finitos asignando todas las condiciones a la que est sometido la estructura completa Figura 2.3c

2.3.3. Campo de desplazamientos Para un slido que cumple las condiciones para ser modelado como un slido axisimtrico, el movimiento de un punto, drefinido en la seccin caracterstica de revolucin, queda perfectamente definido por las componentes de los desplazamiento radial ( ) y axial ( ) Figura 2.3.

(
2.3.4. Campo de deformaciones

(2.1)

El vector de deformaciones de un punto, tiene las cuatro componentes siguiente: r deformacin radial, z deformacin axial, deformacin circunferencial, deformacin tangencial.

( ) (
2.3.5. Campo de tensiones

(2.2)

Las tenciones no nulas se corresponden con las deformaciones no nulas, as, el vector de tenciones se escribe como:

(2.3)

Donde r, z, , rz,son respectivamente, las tensiones radial, axial, circunferencial y tangencial. La Figura 2.4 muestra el convenio de signos para las tenciones actuando sobre un elemento diferencial de un slido de revolucin.

Figura 2.4: Tensiones actuando sobre un elemento diferencial de un slido de revolucin con cargas de revolucin

2.3.6. Ecuacin constitutiva La relacin entre tensiones y deformaciones, se deduce de la elasticidad tridimensional, en presencia de tensiones y deformaciones iniciales se obtiene:

(2.4)

Para materiales istropos, la matriz de propiedades del elemento est dado por la ecuacin (2.5), en la cual E es el mdulo de elasticidad del material y v es el coeficiente de poisson.

(2.5)

2.3.7. Expresin del principio de la energa potencial mnima Considerando el volumen elemental mostrado en la figura (2.4) la energa potencial puede escribirse en la forma.

(2.6) Donde rdld es el rea de la superficie elemental, y la carga puntual Pi representa una carga lineal distribuida alrededor de un crculo, T es una carga distribuida en la superficie de la estructura y f es la fuerza del cuerpo; como se muestra en la Figura 2.3c. Todas las variables en las integrales son independientes de , entonces la ecuacin (2.6) puede escribirse mediante la ecuacin (2.7).

(2.7)

Se nota claramente de las ecuaciones (2.6) y (2.7), aunque el modelo analizado es un elemento plano Figura 2.3c, matemticamente se toma en cuenta la estructura completa en 3D.

2.4. Formulacin de elementos finitos

Se realiza una formulacin general del mtodo de los elementos finitos para el anlisis de estructuras, para iniciar es necesario dividir el dominio completo del modelo en anlisis, con elementos para los cuales son conocidos sus funciones de interpolacin,

este procedimiento es conocido en los libros especializados como tecnologa de elementos.

2.4.1. Funciones de forma Las funciones de forma son las que proporcionan las formas que puede adoptar localmente el campo incgnita. La eleccin de las funciones de forma est condicionado no slo por la forma que toman los dominios locales o elementos finitos, tambin por el tipo de problema de campo que intenta resolver, por ejemplo (elasticidad, plasticidad, transferencia de calor, electromagnetismo) y por la manera de

abordarlo(nmero de campos independientes).

2.4.1.1.

Elemento rectangular Con la finalidad de trabajar sobre un sub dominio comn en los elementos en que ser discretizado el dominio completo, se utiliza la formulacin izoparamtrica, en este dominio normalizado Figura 2.5, las funciones de interpolacin para un elemento cuadriltero de cuatro nodos est dado por la ecuacin (2.8).

(2.8)

Figura 2.5: Sistema de coordenadas cartesianas y normalizado

2.4.2. Coordenadas Las coordenadas en el interior de cada elemento finito estn dadas por la ecuacin (2.9)

(2.9)

2.4.3. Matriz Jacobiano El cambio desde el sistema de referencia global al sistema normalizado est expresado por la matriz jacobiana ecuacin (2.10)

(2.10)

2.4.4. Matriz de rigidez La matriz de rigidez se obtiene desde la ecuacin (2.7), desarrollando para un elemento de n nudos, el vector de desplazamientos para un elemento est ordenado convenientemente en la ecuacin (2.12) y el vector de funciones de forma ecuacin (2.11) es:

(2.11)

(2.12)

La matriz deformacin unitaria desplazamiento.

(2.13)

Las deformaciones estn dadas mediante

Bq, la matriz de rigidez en el

sistema normalizado para una estructura de simetra axial est dado por la ecuacin (2.14) para obtener las expresiones de las cargas se procede de manera similar.

(2.14)

Donde:

t: representa el ngulo de revolucin que gira cada elemento finito en

que fue discretizado el plano de la seccin caractersticas de revolucin.

B: es la matriz deformacin unitaria desplazamiento o simplemente

matriz de deformacin unitaria, ecuacin (2.13).

D: matriz de propiedades del elemento, ecuacin (2.5). r: coordenada radial para cada punto del elemento finito, ecuacin (2.9)

J: matriz jacobiana, ecuacin (2.10)

La mejor manera de obtener la matriz de rigidez, es resolver la ecuacin (2.14) numricamente.

3. Python Aplicado a la Ingeniera Civil

Python es un lenguaje de programacin interpretado creado por Guido van Rossumen el ao 1991. Entre las ventajas que nos da su empleo tenemos: Permite dividir el programa en mdulos reutilizables desde otros programas Python. Viene con una gran coleccin de mdulos estndar que se pueden utilizar como base de los programas. Tambin hay mdulos incluidos que proporcionan E/S de ficheros, llamadas al sistema, sockets y hasta interfaces a GUI (interfaz grfica con el usuario) como Tk,GTK, Qt entre otros. Se utiliza como lenguaje de programacin interpretado, lo que ahorra un tiempo considerable en el desarrollo del programa, pues no es necesario compilar ni enlazar. El intrprete se puede utilizar de modo interactivo, lo que facilita experimentar con caractersticas del lenguaje, escribir programas desechables o probar funciones durante el desarrollo del programa. Posee una licencia de cdigo abierto, denominada Python Software Foundation License, que es compatible con la licencia GPL a partir de la versin 2.1.1, e incompatible en ciertas versiones anteriores. Esta licencia no obliga a liberar el cdigo fuente al distribuir los archivos binarios. La posibilidad de ser un lenguaje libre, multiplataforma y con mdulos que pueden igualar al MATLAB que es comnmente usado en el campo de la ingeniera lo hace una excelente alternativa.

4. Normas Tcnicas empleadas en el diseo de reservorios circulares

En la norma peruana no se cuenta en detalle el diseo de reservorios, por lo tanto se har un principal uso de las normas del ACI: Code requirements for environmental engineering concrete structures (ACI 350) Requisitos de reglamento para concreto estructural (ACI 318S-08)

PARTE III

MARCO METODOLGICO

5. MARCO METODOLGICO
5.1. Condiciones de carga

5.2. Anlisis por el mtodo de los elementos finitos El mtodo de los elementos finitos es uno de los procedimientos que existen para aproximar el comportamiento de una estructura con infinitos grados de libertad por el de otra, con aproximadamente las mismas propiedades fsicas y geomtricas pero con un nmero finito de grados de libertad, cuyas ecuaciones de equilibrio pueden expresarse por un sistema algebraico de ecuaciones simultneas con un nmero limitado de incgnitas 5.2.1. Estructuras de simetra axial Las estructuras de revolucin representan un alto porcentaje en las estructuras de uso comn en ingeniera, como en nuestro caso. Figura 5.1. 5.2.1.1. Modelamiento Estructural

Es posible modelar la estructura con distintas formulaciones del mtodo de elementos finitos, una manera es considerando la formulacin como slidos de revolucin debido a que la estructura es simtrico respecto al eje z del sistema de referencia Figura 5.2, en la cual se elige una seccin caracterstica de revolucin, para esta seccin ser necesario girar 360 sexagesimales para formar la estructura entera, definida la seccin caracterstica de revolucin se le asignar todas las condiciones a las que se encuentra sometido la estructura completa y se analizar sobre esta seccin plana Figura 5.3.

Figura 5.1

Figura 5.2

Figura 5.3 Luego de definir el sistema de referencia global, se deber discretizar el dominio completo, se usarn elementos rectangulares de cuatro nodos y se ordenarn los argumentos de acuerdo a lo establecido para las coordenadas de los nudos, elementos, asignacin de propiedades, soportes y cargas.

Nudos.- Las coordenadas de los nudos deben ser establecidas en coordenadas globales Elementos.- Cada elemento debe ser identificado con cuatro nudos. Propiedades.- Se considerar que todos los elementos poseen las mismas propiedades, esto es requerido para los slidos de revolucin. E = 2.1*109; mdulo de elasticidad del material en Kg/m2 v = 0.3; coeficiente de poisson

t = 2 , espeso del slido de revolucin, expresado en radianes.

(4.1) Soportes.- Se deber de especificar las variables de contorno, definindose los grados de libertad. Cargas.- Se deber de asignar cargas equivalentes a los nudos de la estructura

5.2.1.2.

Proceso de anlisis de la estructura

o Funciones de forma.- para el sistema normalizado, para elementos cuadrilteros de cuatro nodos est dado por:

(4.2)

o Coordenadas.- Las coordenadas radial r y circunferencial z son expresados en funcin de las ecuaciones (4.2) y las coordenadas de los nudos en el sistema global.

(4.3)

o Desplazamiento.- Los desplazamiento en el interior del elemento son dependientes de las funciones de forma ecuaciones (4.2) y los desplazamiento de los nudos y est expresado mediante:

(4.4) o Matriz Jacobiana.- Representa el cambio de coordenadas ecuacin (4.5), depende de las derivadas parciales de las coordenadas globales expresadas en funcin del sistema normalizado.

(4.5) o Matriz Deformacin unitarias Unitaria-Desplazamiento.para elementos Las

deformaciones

estructurales

modelados como slidos de revolucin, est dado por la ecuacin (4.6), donde r es la deformacin normal en la direccin radial, z es la deformacin normal en la direccin axial, rz deformacin tangencial en el plano rz y es la deformacin normal en la direccin circunferencial.

(4.6) La relacin entre las deformaciones unitarias y los desplazamientos u est dado mediante:

(4.7) La matriz de deformacin unitaria depende de la inversa de la matriz jacobiana y las derivadas parciales de los desplazamientos expresados en el sistema normalizado para cada elemento en cada direccin.

(4.8)

La matriz B de la ecuacin (4.8), se divide en cuatro subamtrices, correspondientes a los cuatro nudos

Expresando todas las submatrices en una sola, se obtiene la matriz de deformacinunitaria desplazamiento ecuacin (4.8)

Donde, ecuacin (4.5).

es el jacobiano, determinante de la matriz jacobiana

o la matriz de rigidez del elemento.- Haciendo uso de las expresiones anteriores, la matriz de rigidez puede escribirse como:

(4.9) La solucin de la ecuacin (4.9) es demorosa inclusive para un programa de clculo, por lo tanto se evaluar numricamente mediante la cuadratura de Gauss, para lo cual se toma dos puntos por cada direccin sobre el dominio normalizado, los llamados puntos de Gauss, se evala la matriz k(m) para los cuatro puntos considerados, siendo la matriz del elemento la sumatoria de las matrices parciales Con la finalidad de evaluar numricamente la matriz de rigidez ecuacin (4.9), se define:

Los puntos de Gauss

en la direccin radial,

en la en la

direccin axial y los pesos asociados a cada punto, direccin radial y

en la direccin axial. Por comodidad se

expresan en vectores de cuatro elementos, bien no podra considerarse los pesos para cada direccin, ya que coincidentemente para un elemento cuadriltero de cuatro nodos el factor peso es la unidad.

Finalmente, la matriz de rigidez para un elemento m est dado por:

Evaluando la funcin para el elemento m=1 se obtiene el resultado siempre una matriz cuadrada simtrica de orden 8. Se procede de la misma manera para todos los elementos.

o Ensamblaje de la matriz de rigidez.- una vez obtenidas las matrices k(m) se procede al ensamble de la matriz de rigidez de la estructura. o Vector de fuerzas nodales equivalentes.- Se deber de construir un vector columna ordenada de acuerdo a los grados de libertad propuestos. o Desplazamiento en nudos.- Se procede al clculo de los desplazamientos de los nudos, por medio de la matriz K y el vector F (matriz de rigidez y de fuerzas respectivamente). o Reacciones en los apoyos.- Se procede al clculo de las reacciones en los apoyos, mediante R=K*Q-F, donde K es la matriz de rigidez de toda la estructura , Q los desplazamientos y F es el vector de cargas equivalentes, luego se deber de ordenar las reacciones convenientemente. o Tensiones.- Para obtener la tensin en un punto interior del elemento m se requiere los desplazamiento en los nudos del elemento.

PARTE IV

PROCEDIMIENTO DE CLCULO
6. PROCEDIMIENTO DE CLCULO

PARTE V

ANALISIS DE LOS RESULTADOS

7. CONCLUSIONES Y RECOMENDACIONES
Se ha considerado solo un anlisis esttico, por lo tanto se han considerado 3 cargas: la generada por peso propio, el peso hidrosttico y un peso distribuido sobre la cpula para muy posibles operaciones. Dichas cargas en conjunto y sus combinaciones han sido analizadas sin encontrarse valores excesivos de desplazamiento, an en la peor situacin. El diseo del reservorio se ha basado en un trabajo previo de concreto armado, comprobndose que prcticamente no existen deformaciones considerables como se puede observar en los resultados mostrados en el vector de desplazamientos. El modelamiento de slidos generados a partir de un perfil de revolucin presentan una facilidad de ser llevados a un plano, evitando as discretizar la estructura en slidos espaciales. Facilitndose el trabajo. El empleo de la cuadratura de gauss hizo posible el clculo de las integrales, que inclusive para un computador es bastante pesado. La programacin del mtodo de elementos finitos es bastante similar entre sus modelamientos, por lo tanto es posible una fcil programacin.

8. REFERENCIAS BIBLIOGRFICAS

Anlisis y diseo de un reservorio apoyado de concreto armado para una capacidad de 115 m3, www.arivte.com/comunidad

Gua para el diseo y construccin de reservorios apoyados, Organizacin Panamericana de la Salud CEPIS, Lima 2004

Anlisis aplicando elementos finitos, Larry J. Segerlind, Michigan State university

Anlisis del techo de un auditorio con FEM

Interfaz de programacin para aplicaciones sap2000 y visual C# fundamentado en mathcad, Edmundo Canchari Gutierrez, 2009

Anlisis de un reservorio por FEM, Edmundo Canchari Gutierrez, 2009