METODO DEL ELEMENTO INFINITO
El metodo del elemento finito ofrece una alternativa que es mas adecuada para
tales sistemas. A diferencia de las tecnicas por diferencias finitas, la tecnica del
elemento finito divide el dominio de la solucion en regiones con formas sencillas o
[Link] puede desarrollar una solucion aproximada de la EDP para cada
uno de estos elementos. La solucion total se genera uniendo, o ensamblando,
las soluciones individuales, teniendo cuidado de asegurar la continuidad de las
fronteras entre los elementos. De este modo, la EDP se satisface por secciones.
El uso de elementos, en lugar de una malla rectangular, proporciona una mejor
aproximacion para sistemas con forma irregular. Ademas, se pueden generar
continuamente valores de las incognitas a traves de todo el dominio de la solucion
en lugar de puntos aislados.
La implementacion del metodo del elemento finito usualmente sigue un
procedimiento estandar paso a paso. A continuacion se presenta un panorama
general de cada uno de estos pasos.
31.1.1 Discretizacion
Este paso consiste en dividir el dominio de la solucion en elementos finitos. En la
figura 31.2 se muestran ejemplos de los elementos empleados en una, dos y tres
dimensiones. Los puntos de interseccion de las lineas que forman los lados de los
elementos se conocen como nodos, y los mismos lados se denominan lineas o
planos nodales.
31.1.2 Ecuaciones de
los elementos
El
siguiente
paso
consiste en desarrollar
ecuaciones
para
aproximar la solucion de cada elemento y consta de dos pasos. Primero, se debe
elegir una funcion apropiada con coeficientes desconocidos que aproximara la
solucion. Segundo, se evaluan los coeficientes de modo que la funcion aproxime
la solucion de manera optima.
Eleccion de las funciones de aproximacion. Debido a que son faciles de
manipular matematicamente, a menudo se utilizan polinomios para este proposito.
En el caso uni- dimensional, la alternativa mas sencilla es un polinomio de primer
grado o linea recta.
donde u(x) = la variable dependiente, a 0 y a1= constantes y x = la variable
independiente. Esta funcion debe pasar a traves de los valores de u(x) en los
puntos extremos del elemento x1 y x2. Por lo tanto,
donde u1 = u(x1) y u2 = u(x2). De estas ecaciones, usando la regla de Cramer, se
obtiene
Estos resultados se sustituyen en la ecuacion (31.1) la cual, despues de reagrupar
terminos, se escribe como
donde,
La ecuacion (31.2) se conoce como una funcion de aproximacion, o de forma, N 1
y N2 se denomina funciones de interpolacion. Una inspeccion cuidadosa revela
que la escion(31.2) es, en realidad, el polinomio de interpolacion de primer grado
de Lagrange. Esta ecuacion ofrece un medio para predecir calores intermedios (es
decir, para interpolar) entre valores dados u 1 y u2 en los [Link] figura 31.3
muestra la funcion de forma junto con las funciones de interpolacion
correspondientes. Observe que la suma de las funciones de interpolacion es igual
a uno. Ademas, el hecho de que estemos tratando con ecuaciones lineales facilita
las ope- raciones como la diferenciacion y la integracion. Tales manipulaciones
seran importantes en secciones posteriores. La derivada de la ecuacion (31.2) es
(31.5)
De acuerdo con las ecuaciones (31.3) y (31.4), las derivadas de las N se
calculan como sigue:
(31.6)
y, por lo tanto, la derivada de u es
(31.7)
En otras palabras, es una diferencia dividida que representa la pendiente de la
linea recta que une los nodos.
La integral se expresa como:
Cada uno de los terminos del lado derecho es simplemente la integral de un
triangulo rectangulo con base x x y altura u. Es decir,
Asi, la integral completa es
( 31.8)
En otras palabras, esto es simplemente la regla del trapecio.
Obtencion de un ajuste optimo de la funcion a la solucion.
Una vez que se ha elegido la funcion de interpolacion, se debe desarrollar la
ecuacion que rige el compor- tamiento del elemento. Esta ecuacion representa un
ajuste de la funcion a la solucion de la ecuacion diferencial de que se trate. Existen
varios metodos para este proposito; entre los mas comunes estan el metodo
directo, el metodo de los residuos ponderados y el metodo variacional. Los
resultados de todos estos metodos son analogos al ajuste de curvas. Sin
embargo, en lugar de ajustar funciones a datos, estos metodos especifican
relaciones entre las incognitas de la ecuacion (31.2) que satisfacen de manera
optima la EDP.
Matematicamente, las ecuaciones del elemento resultante a menudo consisten en
un sistema de ecuaciones algebraicas lineales que puede expresarse en forma
matricial,
[k]{u} = {F}
(31.9)
donde [k] = una propiedad del elemento o matriz de rigidez, {u} = vector columna
de las incognitas en los nodos y {F} = vector columna determinado por el efecto de
cualquier influencia externa aplicada a los nodos. Observe que, en algunos casos,
las ecuaciones pueden ser no lineales. Sin embargo, en los ejemplos elementales
descritos aqui, asi como en muchos problemas practicos, los sistemas son
lineales.
31.1.3 Ensamble
Una vez obtenidas las ecuaciones de elementos individuales, estas deben unirse o
ensam- blarse para caracterizar el comportamiento de todo el sistema. El proceso
de ensamble esta regido por el concepto de continuidad. Es decir, las soluciones
de elementos contiguos se acoplan, de manera que los valores de las incognitas
(y algunas veces las derivadas) en sus nodos comunes sean equivalentes. Asi, la
solucion total sera continua. Cuando finalmente todas las versiones individuales de
la ecuacion (31.9) estan ensambladas, el sistema completo se expresa en forma
matricial como
[K]{u} = {F}
(31.10)
donde [K] = la matriz de propiedades de ensamble y {u} y {F} = vectores columna
de las incognitas y de las fuerzas externas, marcadas con apostrofos para denotar
que son ensamble de los vectores {u} y {F} de los elementos individuales.
31.1.4 Condiciones de frontera
Antes de resolver la ecuacion (31.10) debe modificarse para considerar las
condiciones de frontera del sistema. Dichos ajustes dan como resultado
[K]{u} = {F}
(31.11)
donde la barra significa que las condiciones de frontera se han incorporado.
31.1.5 Solucion
Las soluciones de la ecuacion (31.11) se obtienen con las tecnicas que se
describieron en la parte tres, tal como la descomposicion LU. En muchos casos,
los elementos pueden configurarse de manera que las ecuaciones resultantes
sean bandeadas. Asi, es posible utilizar los esquemas de solucion altamente
eficientes para estos sistemas.
31.1.6 Procesamiento posterior
Una vez obtenida la solucion, esta se despliega en forma tabular o de manera grafica. Ademas,
pueden determinarse las variables secundarias y tambien mostrarse.
Aunque los pasos anteriores son muy generales, son comunes a la mayoria de las
implementaciones del metodo del elemento finito. En la siguiente seccion ilustraremos como se
aplican para obtener resultados numericos de un sistema fisico simple (una barra calentada).
31.2 APLICACIN DEL ELEMENTO FINITO EN UNA DIMENSIN
En la figura 31.4 se muestra un sistema que puede modelarse mediante la forma unidi- mensional
de la ecuacion de Poisson
(31.12)
donde f(x) = una funcion que define una fuente de calor a lo largo de la barra, y
donde los extremos de la barra se mantienen a temperaturas fijas,
T(0, t) = T1
Y
T(L,t) = T2
Observe que esta no es una ecuacion diferencial parcial, sino una EDO con valor
en la frontera. Se usa este modelo sencillo porque nos permitira introducir el
metodo del elemen- to finito sin algunas de las complicaciones de una EDP,
bidimensional por ejemplo.
31.2.1 Discretizacion
Una configuracion simple para modelar el sistema consiste en una serie de
elementos de igual longitud. Asi, el sistema se trata como cuatro elementos de
igual longitud y cinco nodos.
31.2.2 Ecuaciones de los elementos
El metodo directo. En el caso donde f(x) = 0, se utiliza un metodo directo para
gene- rar las ecuaciones de los elementos. La relacion entre el flujo de calor y el
gradiente de temperatura puede representarse mediante la ley de Fourier :
donde q1 es el flujo de calor en el nodo 1. De manera similar, para el nodo 2,
Estas dos ecuaciones expresan la relacion de la distribucion de la temperatura
interna de los elementos (determinada por las temperaturas nodales) con el flujo
de calor en sus extremos. En consecuencia representan nuestras ecuaciones de
los elementos deseadas. Se simplifican aun mas reconociendo que la ley de
Fourier se puede utilizar para expre- sar los flujos de los extremos en terminos de
los gradientes de temperatura en las fron- teras. Es decir,
que se sustituyen en las ecuaciones de los elementos para dar
(31.4)
Observe que la ecuacion (31.14) se presento en el formato de la ecuacion (31.9). Asi, logramos
generar una ecuacion matricial que describa el comportamiento de un elemen- to tipico de nuestro
sistema.
El metodo directo resulta muy intuitivo. Se utiliza en areas como la mecanica, para resolver
problemas importantes. Sin embargo, en otros contextos a menudo es dificil, si no es que
imposible, obtener directamente las ecuaciones del elemento finito. En con- secuencia, como se
describe a continuacion, se cuenta con tecnicas matematicas mas generales.
El metodo de los residuos ponderados. La ecuacion diferencial (31.12) se reexpresa como
La solucion aproximada [ecuacion (31.13)] se sustituye en esta ecuacion. Como la
ecuacion (31.13) no es la solucion exacta, el lado izquierdo de la ecuacion
resultante no sera cero, sino que sera igual a un residuo,
(31.15)
El metodo de los residuos ponderados (MRP) consiste en encontrar un minimo
para el residuo, de acuerdo con la formula general
(31.16)
donde D = de la solucion y funciones de ponderacion linealmente independientes. Aqui, se tienen. Aqui multiples opciones para las funciones de
ponderacion (cuadro 31.1). El procedimiento mas comun para el metodo del
elemento finito consiste en emplear las funciones de interpolacion Ni como
lasfunciones de ponderacion. Cuando estas se sus- tituyen en la ecuacion (31.16),
el resultado se conoce como el metodo de Galerkin,
En nuestra barra udimensional, la ecuacion (31.15) se sustituye en esta
formulacion para dar
que se pueden reexpresar como sigue:
(31.17)
Ahora, se aplicaran varias manipulaciones matematicas para simplificar y evaluar
la ecuacion (31.17). Una de las mas importantes es la simplificacion del lado
izquierdo usando la integracion por partes. Del calculo, recuerde que esta
operacion se expresa como
Se puede elegir entre varias funciones de ponderacion para la ecuacion (31.16).
Cada una representa un procedimiento alter- nativo para el MRP.
En el metodo de la colocacion, elegimos tantas posiciones como coeficientes
desconocidos existan. Despues, se ajustan los coefi- cientes hasta que los
residuos desaparezcan en cada una de estas posiciones. En consecuencia, la
funcion de aproximacion dara resultados perfectos en las posiciones elegidas,
pero en las posi- ciones restantes tendremos un residuo diferente de cero. Asi,
este metodo es parecido a los de interpolacion del capitulo 18. Observe que la
colocacion corresponde a usar la funcion de ponderacion
donde n= el numero de coeficientes desconocidos y (x x i) = la funcion delta de
Dirac, ques igual a cero en todas partes excepto en x = x i, donde es igual 1.
En el metodo del subdominio, el intervalo se divide en tantos segmentos, o
subdominios, como coeficientes desconocidos existan. Despues, se ajustan los
coeficientes hasta que el valor promedio del residuo sea cero en cada subdominio.
Asi, en cada subdominio, la funcion de ponderacion sera igual a 1, y la ecua- cion
(31.16) se convierte en
donde xi-1 y xi son las fronteras del subdominio.
En el caso de minimos cuadrados, los coeficientes se ajustan hasta minimizar la
integral del cuadrado del residuo. De manera que las funciones de ponderacion
son
al sustituir las Wi en la escuacion (31.16), se obtiene
La comparacion de esta formulacion con la del capitulo 17 muestra que esta es la
forma continua de la regresion.
El metodo de Galerkin emplea las funciones de interpolacion Ni como funciones
de ponderacion. Recuerde que estas funciones siempre suman 1 en cualquier
posicion en un elemento. En muchos problemas, el metodo de Galerkin da los
mismos resultados que los que se obtienen con los metodos variacionales. En
con- secuencia, esta es la version del MRP que se emplea con mas frecuencia en
el analisis del elemento finito.
Si u y v se eligen adecuadamente, la nueva integral en el lado derecho sera mas
facil de evaluar que la integral original del lado izquierdo. Esto se puede hacer
para el termino del lado izquierdo de la ecuacion (31.17), escogiendo N i (x) como
u, y (d2T/dx2) dx como dv se obtiene,
(3
1.18)
Asi, hemos dado el importante paso de bajar el orden en la formulacion: de una
segunda a una primera derivada.
A continuacion, se evalua cada uno de los terminos que hemos creado en la
ecuacion (31.18). Para i = 1, el primer termino del lado derecho de la ecuacion
(31.18) se evalua como sigue
Sin embargo, de la figura 31.3 recuerde que N i(x2) = 0 y N1(x1) = 1 y, por lo tanto,
(31.19)
Demanerasimilar,parai=2,
(31.20)
si, el primer termino en el lado derecho de la ecuacion (31.18) representa las
condicio- nes de frontera naturales en los extremos de los elementos.
Ahora, antes de continuar, reagrupemos sustituyendo en la ecuacion original los
terminos correspondientes por nuestros resultados. Empleamos las ecuaciones
(31.18) a (31.20) para hacer las sustituciones correspondientes en la ecuacion
(31.17); para i = 1,
(31.21)
yparai=2,
(31.22)
Observe que la integracion por partes nos llevo a dos importantes resultados. Primero, ha incorporado las condiciones de frontera directamente dentro de las
ecuaciones del elemento. Segundo, ha bajado la evaluacion de orden superior, de
una segunda a una primera derivada. Este ultimo resultado tiene como
consecuencia significativa que las funciones de aproximacion necesitan preservar
continuidad de valor, pero no pendiente en los nodos.
Observe tambien que ahora podemos comenzar a darles significado fisico a cada
uno de los terminos que obtuvimos. En el lado derecho de cada ecuacion, el
primer termino representa una de las condiciones de frontera del elemento; y el
segundo es el efecto de la funcion de fuerza del sistema, en este caso, la fuente
de calor f(x). Como ahora sera evidente, el lado izquierdo representa los
mecanismos internos que rigen la
distribucion de la temperatura del elemento. Es decir, en terminos del metodo del
ele- mento finito, el lado izquierdo sera la matriz de propiedad del elemento.
Para ver esto nos concentramos en los terminos del lado izquierdo. Para i = 1, el
termino es
(31.23)
Recordemos de la seccion 31.1.2 que la naturaleza lineal de la funcion hace que la
dife- renciacion y la integracion sean sencillas. Si empleamos las ecuaciones
(31.6) y (31.7) para hacer las sustituciones correspondientes en la ecuacion
(31.23), obtenemos
(31.24)
Demanerasimilarparai=2[ecuacion(31.22)],
(31.25)
Una comparacion con la ecuacion (31.14) nos muestra que estas son similares a
las relaciones obtenidas con el metodo directo usando la ley de Fourier, lo cual se
aclara mas al expresar las ecuaciones (31.24) y (31.25) en forma matricial como
sigue:
Si este resultado se sustituye en las ecuaciones (31.21) y (31.22), y despues se
expresa en forma matricial, obtenemos la version final de las ecuaciones de los
elementos.
(31.26)
Observe que las ecuaciones del elemento pueden obtenerse no solo mediante los
metodos directo y de los residuos ponderados, sino tambien usando el calculo de
variaciones (por ejemplo, vease Allaire, 1985). En el caso presente, este metodo
proporciona ecuaciones identicas a las deducidas arriba.
