4.- FUNCIONES DE INTERPOLACIONY CONVERGENCIA.4.1. FORMULACION INTEGRAL DEL PROBLEMA Y METODOS VARIACIONALES.

La idea general del mtodo de los elementos finitos es la divisin de un continuo en un conjunto de pequeos elementos interconectados por una serie de puntos llamados nodos. Las ecuaciones que rigen el comportamiento del continuo regirn tambin el del elemento. De esta forma se consigue pasar de un sistema continuo (infinitos grados de libertad), que es regido por una ecuacin diferencial o un sistema de ecuaciones diferenciales, a un sistema con un nmero de grados de libertad finito cuyo comportamiento se modela por un sistema de ecuaciones, lineales o no. En cualquier sistema a analizar podemos distinguir entre: Dominio. Espacio geomtrico donde se va ha analizar el sistema. Condiciones de contorno. Variables conocidas y que condicionan el cambio del sistema: cargas, desplazamientos, temperaturas, voltaje, focos de calor,... Incgnitas. Variables del sistema que deseamos conocer despus de que las condiciones de contorno han actuados sobre el sistema: desplazamientos, tensiones, temperaturas,...

El mtodo de los elementos finitos supone, para solucionar el problema, el dominio discretizado en subdominios denominados elementos. El dominio se divide mediante puntos (en el caso lineal), mediante lneas (en el caso bidimensional) o superficies ( en el tridimensional) imaginarias, de forma que el dominio total en estudio se aproxime mediante el conjunto de porciones (elementos) en que se subdivide. Los elementos se definen por un nmero discreto de puntos, llamados nodos, que conectan entre si los elementos. Sobre estos nodos se materializan las incgnitas fundamentales del problema. En el caso de elementos estructurales estas incgnitas son los desplazamientos nodales, ya que a partir de stos podemos calcular el resto de incgnitas que nos interesen: tensiones, deformaciones,... A estas incgnitas se les denomina grados de libertad de cada nodo del modelo. Los grados de libertad de un nodo son las variables que nos determinan el estado y/o posicin del nodo.

4.2. FORMA DE LOS ELEMENTOS, NODOS Y VARIABLES NODALES.

4.3. FUNCIONES DE INTERPOLACION. Discretizado el continuo, la idea es tomar un conjunto de funciones (funciones de interpolacin) que definan de manera nica el campo de desplazamientos dentro del elemento en funcin de los desplazamientos en los nodos del mismo. Es decir

Siendo {U} el vector con los desplazamientos nodales. Una vez conocidos los desplazamientos en todos los nodos se determinan las deformaciones

Donde [D] es el operador diferencial que depende del problema en estudio Sustituyendo el valor del desplazamiento tenemos que

Donde se obtiene el valor de las deformaciones en funcin de los desplazamientos nodales. Sustituyendo la ecuacin anterior en la expresin del PTV (1.4) tenemos

Reordenando esta ecuacin podemos llegar a un sistema de la forma

Donde se definen: Matriz de masa consistente

Matriz de rigidez

Matriz de cargas nodales consistentes

La expresin anterior es general y permite determinar las matrices elementales para cualquier tipo de discretizacin.

4.4. UTILIZACION DE POLINOMIOS.-

4.5. APLICACIN AL CASO DE TENSION PLANA.El estado de tensiones en la placa tiene cinco componentes: . Para facilitar el estudio posterior, estas tensiones se agrupan en dos vectores diferentes (figura 7.5): el vector contiene las tres tensiones contenidas en el plano de la placa, y el vector contiene las dos tensiones cortantes perpendiculares a ella.

Dado que z=0 y que slo existen las tres deformaciones unitarias situadas en el plano x,y, la relacin entre las tensiones y dichas deformaciones unitarias es:

Siendo D la matriz elstica y 0 el vector de deformaciones unitarias iniciales. Cuando est originado por una variacin de temperatura T este vector tiene el valor siguiente:

Al ser la variacin de las deformaciones unitarias lineal en el espesor, tambin las tensiones varan de forma lineal en l. Cada punto del material est sometido a un estado plano de tensiones situado en el plano x,y, que vara de forma lineal a lo largo del espesor de la placa. Las tensiones cortantes verticales son proporcionales a las deformaciones unitarias correspondientes , que segn la hiptesis de deformacin efectuada son ambas nulas:

Pero en realidad deben existir tensiones cortantes verticales , para que pueda existir esfuerzo cortante vertical. Esta incoherencia se resuelve considerando que el mdulo de elasticidad de cortadura G es infinito, es decir que la placa no se deforma ante la cortadura vertical, aunque existan tensiones en dicha direccin. Ms adelante se obtendr una expresin de las tensiones cortantes verticales, basndose en criterios de equilibrio en lugar de en la ecuacin constitutiva. 4.6. UTILIZACION DE FUNCIONES TRIGONOMETRICAS.Dado que los distintos mtodos de discretizacin desarrollados en la presente tesis permiten la generacin de mallas triangulares, cuadrilteras (con ciertas restricciones) o mixtas, es interesante describir el mtodo de descomposicin de un elemento cuadriltero en elementos triangulares antes de presentar los distintos procedimientos de generacin de malla. Esta descomposicin consiste en dividir el cuadriltero en dos tringulos por una de sus diagonales, es decir, crear una nueva arista y dos elementos triangulares. Para elegir la diagonal por la que se divide el cuadriltero hay que tener presente el criterio de mxima regularidad posible de los tringulos generados. Como se puede observar en la figura 2.1, este criterio se cumple si se divide por la diagonal ms corta. No obstante, este criterio es mejorable porque, como se puede observar en la figura 2.2, en cuadrilteros con grandes diferencias entre sus ngulos, la eleccin de la diagonal ms corta genera en ocasiones peores resultados. Se puede observar que, en el primer caso, se crea un tringulo bastante regular y otro muy deficiente en ese aspecto. Si se opta por la diagonal larga, los dos tringulos generados no son demasiado regulares, pero s aceptables. Parece claro, por tanto, que la diagonal elegida debe partir del mayor ngulo del cuadriltero. Este criterio es ms general que el de la diagonal corta, coincidiendo con l en la mayora de los casos.

4.7. CONVERGENCIA.Usando elementos finitos obtenemos soluciones aproximadas, idealmente si incrementamos el nmero de elementos empleados en un anlisis, la solucin debera converger a la solucin exacta del problema. Hasta ahora solo se ha requerido que de las funciones de aproximacin sean integrables, pero esto no es suficiente para garantizar la convergencia. Notemos que los integrandos de las integrales de los mtodos de residuos ponderados y de las integrales que provienen de principios variacionales deben ser constantes en el lmite cuando el tamao de los elementos tiende a cero. Esto es, a medida que el tamao de los elementos disminuye estos integrandos deben tender a un valor constante sobre cada elemento. Luego una condicin necesaria para la convergencia de las funciones de aproximacin es que para un tamao finito los elementos finitos deben ser capaces de reproducir cualquier valor constante de los integrandos que pueda ocurrir en el limite del tamao de los elementos tendiendo a cero. As, por ejemplo, para un anlisis de tensiones, la condicin de convergencia implica que un elemento finito debe ser capaz de reproducir cualquier estado de deformacin constante. Ya que los integrandos, en este caso, son funciones de las derivadas primeras de los desplazamientos, la condicin de convergencia implica que las funciones de aproximacin deben ser capaces de reproducir cualquier variacin lineal de desplazamientos sobre el elemento. La construccin de funciones de prueba con continuidad mayor que C0 puede ser bastante complicada, en particular para dos dimensiones. Esto ha hecho que en algunos casos se utilicen elementos con menor continuidad que la mnima requerida, a estos elementos se los denomina incompatibles. La utilizacin de elementos incompatibles requiere que estos elementos cumplan requisitos adicionales para poder ser utilizados con seguridad. Uno de estos requisitos adicionales es el llamado test de la parcela[2] (patch test), que consiste en que el elemento debe poder reproducir a escala finita cualquier comportamiento bsico infinitesimal para asegurar la convergencia.

En general, este test consiste en testear con la computadora estados bsicos a escala finita, y no sirven solamente para verificar la convergencia sino tambin el correcto funcionamiento del programa. 4.8. IDEA DE APROXIMACION.El primer paso para obtener una aproximacin por elementos finitos es realizar una discretizacin del dominio. Esto es, debemos generar una malla de elementos finitos que cubra todo el dominio. Adems debemos numerar los nodos de la malla, que son aquellos puntos que tienen asociadas coordenadas generalizadas. Para el caso particular de anlisis de tensiones las coordenadas generalizadas son los desplazamientos nodales. As para el nodo i sus desplazamientos nodales sern ui y vi. Consideremos una aproximacin por elementos finitos de los desplazamientos u, como

donde ui, vi son los desplazamientos nodales. Cada funcin de prueba Ni se compone de las funciones de forma asociadas al nodo i de todos los elementos que contienen ese nodo, esto es

donde nel es el nmero de elementos de la malla. Adems, los desplazamientos ue, ve en cada elemento se pueden expresar como

donde nnod es el nmero de nodos del elemento. Ntese que en este caso el ndice j se refiere a la numeracin local del nodo en el elemento. As, por ejemplo, para un tringulo de 3 nodos los desplazamientos sobre el elemento son

esta ecuacin puede escribirse matricialmente como

siendo Ne la matriz de funciones de forma del elemento

y de es el vector de desplazamientos nodales del elemento.

En forma particionada la matriz de funciones de forma se puede escribir como

donde las submatrices Nie que estn asociadas a cada nodo del elemento son

El vector de de desplazamientos nodales del elemento se puede tambin expresar en forma particionada como

donde los vectores die que estn asociados a los desplazamientos de cada nodo i del elemento son

siendo uie, vie los desplazamientos del nodo i del elemento

4.9. ANALISIS DE ERRORES EN EL M.E.F. ERRORES DE APROXIMACIN Los errores introducidos al interior de la solucin de una ecuacin diferencial dada pueden atribuirse a tres fuentes bsicas. 1. 2. 3. Error de aproximacin del dominio, debido a la aproximacin del dominio. Errores de cuadratura y aritmtica finita, debidos a la evaluacin numrica de integrales y al clculo numrico en una computadora. Error de aproximacin, debido a la aproximacin de la solucin.

dondeU, representa el valor de z/en el nodo global /, global asociada con el nodo global.

^representa la funcin interpolacin

En los problemas de una dimensin discutidos, los dominios considerados han sido lneas rectas. Por lo tanto la aproximacin del dominio no ha sido necesaria. En problemas de dos dimensiones que involucran dominios no rectangulares, errores de aproximacin del dominio son introducidos en las soluciones del elemento finito. En general pueden interpretarse como errores en la especificacin de los datos del problema debido a que resolvemos ahora la ecuacin diferencial dada sobre un dominio modificado. Como retinemos la malla, el dominio es representado ms exactamente, y, por lo tanto, ios errores de aproximacin frontera se espera se aproximen a cero. Cuando los clculos del elemento finito se desarrollan en una computadora, errores redondeados en el clculo de nmeros y errores debido a la evaluacin numrica de integrales son introducidos en la solucin. En a mayora de los problemas lineales con nmero pequeo razonablemente de grados de libertad totales en el sistema, esos errores se espera sean pequeos (o cero cuando solamente se desea exactitud en un punto decimal). El error introducido en la solucin del elemento finito U debido a la aproximacin de la e variable dependiente u en un elemento Cl es inherente a cualquier problema:
e

Donde w es la solucin de la solucin del elemento finito sobre el dominio (uh -f/ euQ \Nes el nmero de elementos en la malla M es el nmero total de nodos globales, y n es el nmero de nodos en un elemento. Deseamos conocer como el error E = u-uh medido de manera significativa se comporta conforme se incremente el nmero de elementos en la malla. Se puede demostrar que el error aproximacin es cero para ecuaciones simples de segundo orden y de cuarto orden con coeficientes constantes del elemento discreto.