Caracterización cualitativa de algunos métodos de solución de sistemas de ecuaciones lineales. MROGINSKI, Javier L.
Download ReportTranscript Caracterización cualitativa de algunos métodos de solución de sistemas de ecuaciones lineales. MROGINSKI, Javier L.
Slide 1
Caracterización cualitativa de
algunos métodos de solución de
sistemas de ecuaciones lineales.
MROGINSKI, Javier L. – BENEYTO, Pablo A. –
DI RADO, H. Ariel – MANZOLILLO, Juan E.
Departamento de Mecánica Aplicada.
Facultad de Ingeniería. U.N.N.E.
Antecedentes
El Método de los Elementos Finitos (M.E.F.) es una
técnica de análisis numérico que permite la obtención de
soluciones aproximadas de ecuaciones diferenciales capaces
de representar una variedad de problemas de la Ingeniería y
Física.
El resultado de la aplicación del Método de los Elementos
Finitos a un dominio cualquier es un sistema de ecuaciones
lineales simultáneas de gran tamaño, correspondiendo una
ecuación por cada incógnita que se evalúe en cada nodo,
imposible de ser resuelto sin el auxilio de un computador.
El presente trabajo desarrolla algunos aspectos
significativos sobre los tres diferentes métodos de resolución
de sistemas de ecuaciones lineales simultáneas mas usados en
el Método de los Elementos Finitos que son el: Eliminación
de Gauss, Sistema Banda y Frontal entre otros.
Materiales y Métodos
En las diapositivas siguientes
expondrá una breve descripción de
métodos de solución de sistemas
ecuaciones lineales simultáneas que
estudian en el presente trabajo.
se
los
de
se
Materiales y Métodos
Eliminación de Gauss
Un sistema de ecuaciones lineales simultáneas
puede ser expresado en forma genérica de la
siguiente manera:
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
a 21 x 1 a 22 x 2 a 23 x 3 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... .
a n1 x 1 a n 2 x 2 a n 3 x 3 a n n x n C n
Materiales y Métodos
Eliminación de Gauss
Este método busca obtener un sistema de ecuaciones
triangular equivalente, mediante el uso de transformaciones
lineales. El cual se resuelve con facilidad aplicando la
denominada sustitución inversa.
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
´
´
´
´
a 22 x 2 a 23 x 2 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... ........
(n2)
(n2)
(n2)
a n 1 , n 1 x n 1 a n 1 , n x n C n 1
( n 1 )
( n 1 )
a n,n x n C n
Materiales y Métodos
Sistema banda de solución
Se dice que una matriz es matriz banda cuando posee
todos sus elementos no nulos concentrados en torno de
la diagonal principal. Esta es una propiedad muy
interesante que se presenta en las matrices de rigidez
locales y globales.
El método banda de resolución de sistemas de
ecuaciones lineales simultaneas se basa en esta
propiedad, permitiendo el ahorro de asignaciones de
memoria, ya que solo guarda los elementos no nulos.
Materiales y Métodos
Sistema banda de solución
El presente método centra su eficiencia en la forma en que
se realiza la numeración global de los nodos.
S em ian ch o d e b an d a
N ro. total de incognitas
X X X
0 0 0 0
X X X X 0 0 0
X X X X X 0 0
0 X X X X X 0
0 0 X X X X X
0 0 0 X X X X
0 0 0 0 X X X
0 0 0 0 0 X X
0 0 0 0 0 0 X
0 0 0 0 0 0 0
0 0 0 0 0 0 0
0 0 0 0 0 0 0
A n ch o d e b an d a
0
0
0
0
0
0
0
0
0
0
X 0
0
0
0
0
0
0
X X 0
0
0
0
0
0
0
0
X X X 0
0
0
0
0
0
0
0
0
X X X X 0
0 0 X X X
0 X X X X
A n ch o d e b an d a
D iag o n al
p rin cip al
X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X X X
X X X
0 X X X X
0 0 X X X
X X X X
0
X X X 0 0
X X
O rg an izació n co n v en cio n al d e la
m atriz d e rig id ez
0
X 0 0
O rg an izació n en b an d a p ara O rg an izació n en b an d a p ara
m atrices n o sim étricas
m atrices sim étricas
Materiales y Métodos
Método Frontal
Este método trata de evitar el ensamblaje de la matriz
de rigidez global. Consiste básicamente en ensamblar la
menor cantidad de elementos y eliminar las variables al
mismo tiempo, formando matrices de menor tamaño.
Esto es, inmediatamente después de que los coeficientes
de una ecuación se ensamblan a partir de la contribución
de todos los elementos relevantes, o que inciden sobre este
coeficiente, la variable de este puede eliminarse. Como
resultado, la matriz completa de toda la estructura jamás
se forma.
Materiales y Métodos
Método Frontal
A diferencia del método banda, la adecuada
enumeración de los elemento es de gran importancia, no
así la enumeración global de los nodos ya que estos son
tratados por su numeración local. Como primer medida
esta característica del Front lo hace mas versátil ya que
como la cantidad de elementos siempre es menor que la
cantidad de nodos, independientemente de la forma de la
estructura, siempre se tendrá mayos facilidad para
enumerar los elementos que los nodos.
Discusión de Resultados
Para analizar las cualidades de los distintos métodos
propuestos se planteó la resolución mediante el empleo del
FECCUND V 2.1 (software desarrollado por el Dto. De
Mecánica Aplicada de la Facultad de Ingeniería de la
UNNE) de dos ejemplos prácticos, representando los
resultados mediante el pre-post procesador GID 7.2.
Para dicho análisis se tomaró como variables el tiempo
de ejecución de cada método y su exactitud, la cual será
determinada por comparación con los resultados
obtenidos por el software denominado CALSEF provisto
por el Centro Internacional de Métodos Numéricos en
Ingeniería (CIMNE).
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Mallado con elementos
hexaédricos
Desplazamiento nodal
según eje y
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Tabla Nº 1
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
1,758
0,501953E-04
0,47619E-04
5,13
Banda
0,552
0,498874E-04
4,55
Frontal
0,441
0,479135E-04
0,61
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Mallado con elementos hexaédricos
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Desplazamiento nodal según eje y
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Tabla Nº 2
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
31,547
-0,00615923
-0.00544118
11,66
Banda
2,025
-0,00560842
2,98
Frontal
1,752
-0,00552264
1,47
Conclusiones
El método Frontal es superior al Gauss y al Banda
en cuanto a su exactitud y a la velocidad de
ejecución.
El efecto de la propagación de errores por redondeo
incrementa notablemente la deficiencia del método
de Gauss frente a problemas de mayor
envergadura.
El método de Banda puede ser recomendable en el
caso de contar con un software de pre-post
procesamiento ya que los mismos generan las
mallas enumerando los nodos en la forma mas
adecuada disminuyendo de ancho de banda.
Slide 2
Caracterización cualitativa de
algunos métodos de solución de
sistemas de ecuaciones lineales.
MROGINSKI, Javier L. – BENEYTO, Pablo A. –
DI RADO, H. Ariel – MANZOLILLO, Juan E.
Departamento de Mecánica Aplicada.
Facultad de Ingeniería. U.N.N.E.
Antecedentes
El Método de los Elementos Finitos (M.E.F.) es una
técnica de análisis numérico que permite la obtención de
soluciones aproximadas de ecuaciones diferenciales capaces
de representar una variedad de problemas de la Ingeniería y
Física.
El resultado de la aplicación del Método de los Elementos
Finitos a un dominio cualquier es un sistema de ecuaciones
lineales simultáneas de gran tamaño, correspondiendo una
ecuación por cada incógnita que se evalúe en cada nodo,
imposible de ser resuelto sin el auxilio de un computador.
El presente trabajo desarrolla algunos aspectos
significativos sobre los tres diferentes métodos de resolución
de sistemas de ecuaciones lineales simultáneas mas usados en
el Método de los Elementos Finitos que son el: Eliminación
de Gauss, Sistema Banda y Frontal entre otros.
Materiales y Métodos
En las diapositivas siguientes
expondrá una breve descripción de
métodos de solución de sistemas
ecuaciones lineales simultáneas que
estudian en el presente trabajo.
se
los
de
se
Materiales y Métodos
Eliminación de Gauss
Un sistema de ecuaciones lineales simultáneas
puede ser expresado en forma genérica de la
siguiente manera:
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
a 21 x 1 a 22 x 2 a 23 x 3 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... .
a n1 x 1 a n 2 x 2 a n 3 x 3 a n n x n C n
Materiales y Métodos
Eliminación de Gauss
Este método busca obtener un sistema de ecuaciones
triangular equivalente, mediante el uso de transformaciones
lineales. El cual se resuelve con facilidad aplicando la
denominada sustitución inversa.
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
´
´
´
´
a 22 x 2 a 23 x 2 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... ........
(n2)
(n2)
(n2)
a n 1 , n 1 x n 1 a n 1 , n x n C n 1
( n 1 )
( n 1 )
a n,n x n C n
Materiales y Métodos
Sistema banda de solución
Se dice que una matriz es matriz banda cuando posee
todos sus elementos no nulos concentrados en torno de
la diagonal principal. Esta es una propiedad muy
interesante que se presenta en las matrices de rigidez
locales y globales.
El método banda de resolución de sistemas de
ecuaciones lineales simultaneas se basa en esta
propiedad, permitiendo el ahorro de asignaciones de
memoria, ya que solo guarda los elementos no nulos.
Materiales y Métodos
Sistema banda de solución
El presente método centra su eficiencia en la forma en que
se realiza la numeración global de los nodos.
S em ian ch o d e b an d a
N ro. total de incognitas
X X X
0 0 0 0
X X X X 0 0 0
X X X X X 0 0
0 X X X X X 0
0 0 X X X X X
0 0 0 X X X X
0 0 0 0 X X X
0 0 0 0 0 X X
0 0 0 0 0 0 X
0 0 0 0 0 0 0
0 0 0 0 0 0 0
0 0 0 0 0 0 0
A n ch o d e b an d a
0
0
0
0
0
0
0
0
0
0
X 0
0
0
0
0
0
0
X X 0
0
0
0
0
0
0
0
X X X 0
0
0
0
0
0
0
0
0
X X X X 0
0 0 X X X
0 X X X X
A n ch o d e b an d a
D iag o n al
p rin cip al
X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X X X
X X X
0 X X X X
0 0 X X X
X X X X
0
X X X 0 0
X X
O rg an izació n co n v en cio n al d e la
m atriz d e rig id ez
0
X 0 0
O rg an izació n en b an d a p ara O rg an izació n en b an d a p ara
m atrices n o sim étricas
m atrices sim étricas
Materiales y Métodos
Método Frontal
Este método trata de evitar el ensamblaje de la matriz
de rigidez global. Consiste básicamente en ensamblar la
menor cantidad de elementos y eliminar las variables al
mismo tiempo, formando matrices de menor tamaño.
Esto es, inmediatamente después de que los coeficientes
de una ecuación se ensamblan a partir de la contribución
de todos los elementos relevantes, o que inciden sobre este
coeficiente, la variable de este puede eliminarse. Como
resultado, la matriz completa de toda la estructura jamás
se forma.
Materiales y Métodos
Método Frontal
A diferencia del método banda, la adecuada
enumeración de los elemento es de gran importancia, no
así la enumeración global de los nodos ya que estos son
tratados por su numeración local. Como primer medida
esta característica del Front lo hace mas versátil ya que
como la cantidad de elementos siempre es menor que la
cantidad de nodos, independientemente de la forma de la
estructura, siempre se tendrá mayos facilidad para
enumerar los elementos que los nodos.
Discusión de Resultados
Para analizar las cualidades de los distintos métodos
propuestos se planteó la resolución mediante el empleo del
FECCUND V 2.1 (software desarrollado por el Dto. De
Mecánica Aplicada de la Facultad de Ingeniería de la
UNNE) de dos ejemplos prácticos, representando los
resultados mediante el pre-post procesador GID 7.2.
Para dicho análisis se tomaró como variables el tiempo
de ejecución de cada método y su exactitud, la cual será
determinada por comparación con los resultados
obtenidos por el software denominado CALSEF provisto
por el Centro Internacional de Métodos Numéricos en
Ingeniería (CIMNE).
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Mallado con elementos
hexaédricos
Desplazamiento nodal
según eje y
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Tabla Nº 1
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
1,758
0,501953E-04
0,47619E-04
5,13
Banda
0,552
0,498874E-04
4,55
Frontal
0,441
0,479135E-04
0,61
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Mallado con elementos hexaédricos
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Desplazamiento nodal según eje y
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Tabla Nº 2
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
31,547
-0,00615923
-0.00544118
11,66
Banda
2,025
-0,00560842
2,98
Frontal
1,752
-0,00552264
1,47
Conclusiones
El método Frontal es superior al Gauss y al Banda
en cuanto a su exactitud y a la velocidad de
ejecución.
El efecto de la propagación de errores por redondeo
incrementa notablemente la deficiencia del método
de Gauss frente a problemas de mayor
envergadura.
El método de Banda puede ser recomendable en el
caso de contar con un software de pre-post
procesamiento ya que los mismos generan las
mallas enumerando los nodos en la forma mas
adecuada disminuyendo de ancho de banda.
Slide 3
Caracterización cualitativa de
algunos métodos de solución de
sistemas de ecuaciones lineales.
MROGINSKI, Javier L. – BENEYTO, Pablo A. –
DI RADO, H. Ariel – MANZOLILLO, Juan E.
Departamento de Mecánica Aplicada.
Facultad de Ingeniería. U.N.N.E.
Antecedentes
El Método de los Elementos Finitos (M.E.F.) es una
técnica de análisis numérico que permite la obtención de
soluciones aproximadas de ecuaciones diferenciales capaces
de representar una variedad de problemas de la Ingeniería y
Física.
El resultado de la aplicación del Método de los Elementos
Finitos a un dominio cualquier es un sistema de ecuaciones
lineales simultáneas de gran tamaño, correspondiendo una
ecuación por cada incógnita que se evalúe en cada nodo,
imposible de ser resuelto sin el auxilio de un computador.
El presente trabajo desarrolla algunos aspectos
significativos sobre los tres diferentes métodos de resolución
de sistemas de ecuaciones lineales simultáneas mas usados en
el Método de los Elementos Finitos que son el: Eliminación
de Gauss, Sistema Banda y Frontal entre otros.
Materiales y Métodos
En las diapositivas siguientes
expondrá una breve descripción de
métodos de solución de sistemas
ecuaciones lineales simultáneas que
estudian en el presente trabajo.
se
los
de
se
Materiales y Métodos
Eliminación de Gauss
Un sistema de ecuaciones lineales simultáneas
puede ser expresado en forma genérica de la
siguiente manera:
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
a 21 x 1 a 22 x 2 a 23 x 3 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... .
a n1 x 1 a n 2 x 2 a n 3 x 3 a n n x n C n
Materiales y Métodos
Eliminación de Gauss
Este método busca obtener un sistema de ecuaciones
triangular equivalente, mediante el uso de transformaciones
lineales. El cual se resuelve con facilidad aplicando la
denominada sustitución inversa.
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
´
´
´
´
a 22 x 2 a 23 x 2 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... ........
(n2)
(n2)
(n2)
a n 1 , n 1 x n 1 a n 1 , n x n C n 1
( n 1 )
( n 1 )
a n,n x n C n
Materiales y Métodos
Sistema banda de solución
Se dice que una matriz es matriz banda cuando posee
todos sus elementos no nulos concentrados en torno de
la diagonal principal. Esta es una propiedad muy
interesante que se presenta en las matrices de rigidez
locales y globales.
El método banda de resolución de sistemas de
ecuaciones lineales simultaneas se basa en esta
propiedad, permitiendo el ahorro de asignaciones de
memoria, ya que solo guarda los elementos no nulos.
Materiales y Métodos
Sistema banda de solución
El presente método centra su eficiencia en la forma en que
se realiza la numeración global de los nodos.
S em ian ch o d e b an d a
N ro. total de incognitas
X X X
0 0 0 0
X X X X 0 0 0
X X X X X 0 0
0 X X X X X 0
0 0 X X X X X
0 0 0 X X X X
0 0 0 0 X X X
0 0 0 0 0 X X
0 0 0 0 0 0 X
0 0 0 0 0 0 0
0 0 0 0 0 0 0
0 0 0 0 0 0 0
A n ch o d e b an d a
0
0
0
0
0
0
0
0
0
0
X 0
0
0
0
0
0
0
X X 0
0
0
0
0
0
0
0
X X X 0
0
0
0
0
0
0
0
0
X X X X 0
0 0 X X X
0 X X X X
A n ch o d e b an d a
D iag o n al
p rin cip al
X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X X X
X X X
0 X X X X
0 0 X X X
X X X X
0
X X X 0 0
X X
O rg an izació n co n v en cio n al d e la
m atriz d e rig id ez
0
X 0 0
O rg an izació n en b an d a p ara O rg an izació n en b an d a p ara
m atrices n o sim étricas
m atrices sim étricas
Materiales y Métodos
Método Frontal
Este método trata de evitar el ensamblaje de la matriz
de rigidez global. Consiste básicamente en ensamblar la
menor cantidad de elementos y eliminar las variables al
mismo tiempo, formando matrices de menor tamaño.
Esto es, inmediatamente después de que los coeficientes
de una ecuación se ensamblan a partir de la contribución
de todos los elementos relevantes, o que inciden sobre este
coeficiente, la variable de este puede eliminarse. Como
resultado, la matriz completa de toda la estructura jamás
se forma.
Materiales y Métodos
Método Frontal
A diferencia del método banda, la adecuada
enumeración de los elemento es de gran importancia, no
así la enumeración global de los nodos ya que estos son
tratados por su numeración local. Como primer medida
esta característica del Front lo hace mas versátil ya que
como la cantidad de elementos siempre es menor que la
cantidad de nodos, independientemente de la forma de la
estructura, siempre se tendrá mayos facilidad para
enumerar los elementos que los nodos.
Discusión de Resultados
Para analizar las cualidades de los distintos métodos
propuestos se planteó la resolución mediante el empleo del
FECCUND V 2.1 (software desarrollado por el Dto. De
Mecánica Aplicada de la Facultad de Ingeniería de la
UNNE) de dos ejemplos prácticos, representando los
resultados mediante el pre-post procesador GID 7.2.
Para dicho análisis se tomaró como variables el tiempo
de ejecución de cada método y su exactitud, la cual será
determinada por comparación con los resultados
obtenidos por el software denominado CALSEF provisto
por el Centro Internacional de Métodos Numéricos en
Ingeniería (CIMNE).
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Mallado con elementos
hexaédricos
Desplazamiento nodal
según eje y
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Tabla Nº 1
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
1,758
0,501953E-04
0,47619E-04
5,13
Banda
0,552
0,498874E-04
4,55
Frontal
0,441
0,479135E-04
0,61
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Mallado con elementos hexaédricos
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Desplazamiento nodal según eje y
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Tabla Nº 2
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
31,547
-0,00615923
-0.00544118
11,66
Banda
2,025
-0,00560842
2,98
Frontal
1,752
-0,00552264
1,47
Conclusiones
El método Frontal es superior al Gauss y al Banda
en cuanto a su exactitud y a la velocidad de
ejecución.
El efecto de la propagación de errores por redondeo
incrementa notablemente la deficiencia del método
de Gauss frente a problemas de mayor
envergadura.
El método de Banda puede ser recomendable en el
caso de contar con un software de pre-post
procesamiento ya que los mismos generan las
mallas enumerando los nodos en la forma mas
adecuada disminuyendo de ancho de banda.
Slide 4
Caracterización cualitativa de
algunos métodos de solución de
sistemas de ecuaciones lineales.
MROGINSKI, Javier L. – BENEYTO, Pablo A. –
DI RADO, H. Ariel – MANZOLILLO, Juan E.
Departamento de Mecánica Aplicada.
Facultad de Ingeniería. U.N.N.E.
Antecedentes
El Método de los Elementos Finitos (M.E.F.) es una
técnica de análisis numérico que permite la obtención de
soluciones aproximadas de ecuaciones diferenciales capaces
de representar una variedad de problemas de la Ingeniería y
Física.
El resultado de la aplicación del Método de los Elementos
Finitos a un dominio cualquier es un sistema de ecuaciones
lineales simultáneas de gran tamaño, correspondiendo una
ecuación por cada incógnita que se evalúe en cada nodo,
imposible de ser resuelto sin el auxilio de un computador.
El presente trabajo desarrolla algunos aspectos
significativos sobre los tres diferentes métodos de resolución
de sistemas de ecuaciones lineales simultáneas mas usados en
el Método de los Elementos Finitos que son el: Eliminación
de Gauss, Sistema Banda y Frontal entre otros.
Materiales y Métodos
En las diapositivas siguientes
expondrá una breve descripción de
métodos de solución de sistemas
ecuaciones lineales simultáneas que
estudian en el presente trabajo.
se
los
de
se
Materiales y Métodos
Eliminación de Gauss
Un sistema de ecuaciones lineales simultáneas
puede ser expresado en forma genérica de la
siguiente manera:
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
a 21 x 1 a 22 x 2 a 23 x 3 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... .
a n1 x 1 a n 2 x 2 a n 3 x 3 a n n x n C n
Materiales y Métodos
Eliminación de Gauss
Este método busca obtener un sistema de ecuaciones
triangular equivalente, mediante el uso de transformaciones
lineales. El cual se resuelve con facilidad aplicando la
denominada sustitución inversa.
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
´
´
´
´
a 22 x 2 a 23 x 2 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... ........
(n2)
(n2)
(n2)
a n 1 , n 1 x n 1 a n 1 , n x n C n 1
( n 1 )
( n 1 )
a n,n x n C n
Materiales y Métodos
Sistema banda de solución
Se dice que una matriz es matriz banda cuando posee
todos sus elementos no nulos concentrados en torno de
la diagonal principal. Esta es una propiedad muy
interesante que se presenta en las matrices de rigidez
locales y globales.
El método banda de resolución de sistemas de
ecuaciones lineales simultaneas se basa en esta
propiedad, permitiendo el ahorro de asignaciones de
memoria, ya que solo guarda los elementos no nulos.
Materiales y Métodos
Sistema banda de solución
El presente método centra su eficiencia en la forma en que
se realiza la numeración global de los nodos.
S em ian ch o d e b an d a
N ro. total de incognitas
X X X
0 0 0 0
X X X X 0 0 0
X X X X X 0 0
0 X X X X X 0
0 0 X X X X X
0 0 0 X X X X
0 0 0 0 X X X
0 0 0 0 0 X X
0 0 0 0 0 0 X
0 0 0 0 0 0 0
0 0 0 0 0 0 0
0 0 0 0 0 0 0
A n ch o d e b an d a
0
0
0
0
0
0
0
0
0
0
X 0
0
0
0
0
0
0
X X 0
0
0
0
0
0
0
0
X X X 0
0
0
0
0
0
0
0
0
X X X X 0
0 0 X X X
0 X X X X
A n ch o d e b an d a
D iag o n al
p rin cip al
X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X X X
X X X
0 X X X X
0 0 X X X
X X X X
0
X X X 0 0
X X
O rg an izació n co n v en cio n al d e la
m atriz d e rig id ez
0
X 0 0
O rg an izació n en b an d a p ara O rg an izació n en b an d a p ara
m atrices n o sim étricas
m atrices sim étricas
Materiales y Métodos
Método Frontal
Este método trata de evitar el ensamblaje de la matriz
de rigidez global. Consiste básicamente en ensamblar la
menor cantidad de elementos y eliminar las variables al
mismo tiempo, formando matrices de menor tamaño.
Esto es, inmediatamente después de que los coeficientes
de una ecuación se ensamblan a partir de la contribución
de todos los elementos relevantes, o que inciden sobre este
coeficiente, la variable de este puede eliminarse. Como
resultado, la matriz completa de toda la estructura jamás
se forma.
Materiales y Métodos
Método Frontal
A diferencia del método banda, la adecuada
enumeración de los elemento es de gran importancia, no
así la enumeración global de los nodos ya que estos son
tratados por su numeración local. Como primer medida
esta característica del Front lo hace mas versátil ya que
como la cantidad de elementos siempre es menor que la
cantidad de nodos, independientemente de la forma de la
estructura, siempre se tendrá mayos facilidad para
enumerar los elementos que los nodos.
Discusión de Resultados
Para analizar las cualidades de los distintos métodos
propuestos se planteó la resolución mediante el empleo del
FECCUND V 2.1 (software desarrollado por el Dto. De
Mecánica Aplicada de la Facultad de Ingeniería de la
UNNE) de dos ejemplos prácticos, representando los
resultados mediante el pre-post procesador GID 7.2.
Para dicho análisis se tomaró como variables el tiempo
de ejecución de cada método y su exactitud, la cual será
determinada por comparación con los resultados
obtenidos por el software denominado CALSEF provisto
por el Centro Internacional de Métodos Numéricos en
Ingeniería (CIMNE).
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Mallado con elementos
hexaédricos
Desplazamiento nodal
según eje y
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Tabla Nº 1
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
1,758
0,501953E-04
0,47619E-04
5,13
Banda
0,552
0,498874E-04
4,55
Frontal
0,441
0,479135E-04
0,61
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Mallado con elementos hexaédricos
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Desplazamiento nodal según eje y
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Tabla Nº 2
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
31,547
-0,00615923
-0.00544118
11,66
Banda
2,025
-0,00560842
2,98
Frontal
1,752
-0,00552264
1,47
Conclusiones
El método Frontal es superior al Gauss y al Banda
en cuanto a su exactitud y a la velocidad de
ejecución.
El efecto de la propagación de errores por redondeo
incrementa notablemente la deficiencia del método
de Gauss frente a problemas de mayor
envergadura.
El método de Banda puede ser recomendable en el
caso de contar con un software de pre-post
procesamiento ya que los mismos generan las
mallas enumerando los nodos en la forma mas
adecuada disminuyendo de ancho de banda.
Slide 5
Caracterización cualitativa de
algunos métodos de solución de
sistemas de ecuaciones lineales.
MROGINSKI, Javier L. – BENEYTO, Pablo A. –
DI RADO, H. Ariel – MANZOLILLO, Juan E.
Departamento de Mecánica Aplicada.
Facultad de Ingeniería. U.N.N.E.
Antecedentes
El Método de los Elementos Finitos (M.E.F.) es una
técnica de análisis numérico que permite la obtención de
soluciones aproximadas de ecuaciones diferenciales capaces
de representar una variedad de problemas de la Ingeniería y
Física.
El resultado de la aplicación del Método de los Elementos
Finitos a un dominio cualquier es un sistema de ecuaciones
lineales simultáneas de gran tamaño, correspondiendo una
ecuación por cada incógnita que se evalúe en cada nodo,
imposible de ser resuelto sin el auxilio de un computador.
El presente trabajo desarrolla algunos aspectos
significativos sobre los tres diferentes métodos de resolución
de sistemas de ecuaciones lineales simultáneas mas usados en
el Método de los Elementos Finitos que son el: Eliminación
de Gauss, Sistema Banda y Frontal entre otros.
Materiales y Métodos
En las diapositivas siguientes
expondrá una breve descripción de
métodos de solución de sistemas
ecuaciones lineales simultáneas que
estudian en el presente trabajo.
se
los
de
se
Materiales y Métodos
Eliminación de Gauss
Un sistema de ecuaciones lineales simultáneas
puede ser expresado en forma genérica de la
siguiente manera:
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
a 21 x 1 a 22 x 2 a 23 x 3 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... .
a n1 x 1 a n 2 x 2 a n 3 x 3 a n n x n C n
Materiales y Métodos
Eliminación de Gauss
Este método busca obtener un sistema de ecuaciones
triangular equivalente, mediante el uso de transformaciones
lineales. El cual se resuelve con facilidad aplicando la
denominada sustitución inversa.
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
´
´
´
´
a 22 x 2 a 23 x 2 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... ........
(n2)
(n2)
(n2)
a n 1 , n 1 x n 1 a n 1 , n x n C n 1
( n 1 )
( n 1 )
a n,n x n C n
Materiales y Métodos
Sistema banda de solución
Se dice que una matriz es matriz banda cuando posee
todos sus elementos no nulos concentrados en torno de
la diagonal principal. Esta es una propiedad muy
interesante que se presenta en las matrices de rigidez
locales y globales.
El método banda de resolución de sistemas de
ecuaciones lineales simultaneas se basa en esta
propiedad, permitiendo el ahorro de asignaciones de
memoria, ya que solo guarda los elementos no nulos.
Materiales y Métodos
Sistema banda de solución
El presente método centra su eficiencia en la forma en que
se realiza la numeración global de los nodos.
S em ian ch o d e b an d a
N ro. total de incognitas
X X X
0 0 0 0
X X X X 0 0 0
X X X X X 0 0
0 X X X X X 0
0 0 X X X X X
0 0 0 X X X X
0 0 0 0 X X X
0 0 0 0 0 X X
0 0 0 0 0 0 X
0 0 0 0 0 0 0
0 0 0 0 0 0 0
0 0 0 0 0 0 0
A n ch o d e b an d a
0
0
0
0
0
0
0
0
0
0
X 0
0
0
0
0
0
0
X X 0
0
0
0
0
0
0
0
X X X 0
0
0
0
0
0
0
0
0
X X X X 0
0 0 X X X
0 X X X X
A n ch o d e b an d a
D iag o n al
p rin cip al
X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X X X
X X X
0 X X X X
0 0 X X X
X X X X
0
X X X 0 0
X X
O rg an izació n co n v en cio n al d e la
m atriz d e rig id ez
0
X 0 0
O rg an izació n en b an d a p ara O rg an izació n en b an d a p ara
m atrices n o sim étricas
m atrices sim étricas
Materiales y Métodos
Método Frontal
Este método trata de evitar el ensamblaje de la matriz
de rigidez global. Consiste básicamente en ensamblar la
menor cantidad de elementos y eliminar las variables al
mismo tiempo, formando matrices de menor tamaño.
Esto es, inmediatamente después de que los coeficientes
de una ecuación se ensamblan a partir de la contribución
de todos los elementos relevantes, o que inciden sobre este
coeficiente, la variable de este puede eliminarse. Como
resultado, la matriz completa de toda la estructura jamás
se forma.
Materiales y Métodos
Método Frontal
A diferencia del método banda, la adecuada
enumeración de los elemento es de gran importancia, no
así la enumeración global de los nodos ya que estos son
tratados por su numeración local. Como primer medida
esta característica del Front lo hace mas versátil ya que
como la cantidad de elementos siempre es menor que la
cantidad de nodos, independientemente de la forma de la
estructura, siempre se tendrá mayos facilidad para
enumerar los elementos que los nodos.
Discusión de Resultados
Para analizar las cualidades de los distintos métodos
propuestos se planteó la resolución mediante el empleo del
FECCUND V 2.1 (software desarrollado por el Dto. De
Mecánica Aplicada de la Facultad de Ingeniería de la
UNNE) de dos ejemplos prácticos, representando los
resultados mediante el pre-post procesador GID 7.2.
Para dicho análisis se tomaró como variables el tiempo
de ejecución de cada método y su exactitud, la cual será
determinada por comparación con los resultados
obtenidos por el software denominado CALSEF provisto
por el Centro Internacional de Métodos Numéricos en
Ingeniería (CIMNE).
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Mallado con elementos
hexaédricos
Desplazamiento nodal
según eje y
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Tabla Nº 1
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
1,758
0,501953E-04
0,47619E-04
5,13
Banda
0,552
0,498874E-04
4,55
Frontal
0,441
0,479135E-04
0,61
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Mallado con elementos hexaédricos
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Desplazamiento nodal según eje y
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Tabla Nº 2
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
31,547
-0,00615923
-0.00544118
11,66
Banda
2,025
-0,00560842
2,98
Frontal
1,752
-0,00552264
1,47
Conclusiones
El método Frontal es superior al Gauss y al Banda
en cuanto a su exactitud y a la velocidad de
ejecución.
El efecto de la propagación de errores por redondeo
incrementa notablemente la deficiencia del método
de Gauss frente a problemas de mayor
envergadura.
El método de Banda puede ser recomendable en el
caso de contar con un software de pre-post
procesamiento ya que los mismos generan las
mallas enumerando los nodos en la forma mas
adecuada disminuyendo de ancho de banda.
Slide 6
Caracterización cualitativa de
algunos métodos de solución de
sistemas de ecuaciones lineales.
MROGINSKI, Javier L. – BENEYTO, Pablo A. –
DI RADO, H. Ariel – MANZOLILLO, Juan E.
Departamento de Mecánica Aplicada.
Facultad de Ingeniería. U.N.N.E.
Antecedentes
El Método de los Elementos Finitos (M.E.F.) es una
técnica de análisis numérico que permite la obtención de
soluciones aproximadas de ecuaciones diferenciales capaces
de representar una variedad de problemas de la Ingeniería y
Física.
El resultado de la aplicación del Método de los Elementos
Finitos a un dominio cualquier es un sistema de ecuaciones
lineales simultáneas de gran tamaño, correspondiendo una
ecuación por cada incógnita que se evalúe en cada nodo,
imposible de ser resuelto sin el auxilio de un computador.
El presente trabajo desarrolla algunos aspectos
significativos sobre los tres diferentes métodos de resolución
de sistemas de ecuaciones lineales simultáneas mas usados en
el Método de los Elementos Finitos que son el: Eliminación
de Gauss, Sistema Banda y Frontal entre otros.
Materiales y Métodos
En las diapositivas siguientes
expondrá una breve descripción de
métodos de solución de sistemas
ecuaciones lineales simultáneas que
estudian en el presente trabajo.
se
los
de
se
Materiales y Métodos
Eliminación de Gauss
Un sistema de ecuaciones lineales simultáneas
puede ser expresado en forma genérica de la
siguiente manera:
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
a 21 x 1 a 22 x 2 a 23 x 3 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... .
a n1 x 1 a n 2 x 2 a n 3 x 3 a n n x n C n
Materiales y Métodos
Eliminación de Gauss
Este método busca obtener un sistema de ecuaciones
triangular equivalente, mediante el uso de transformaciones
lineales. El cual se resuelve con facilidad aplicando la
denominada sustitución inversa.
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
´
´
´
´
a 22 x 2 a 23 x 2 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... ........
(n2)
(n2)
(n2)
a n 1 , n 1 x n 1 a n 1 , n x n C n 1
( n 1 )
( n 1 )
a n,n x n C n
Materiales y Métodos
Sistema banda de solución
Se dice que una matriz es matriz banda cuando posee
todos sus elementos no nulos concentrados en torno de
la diagonal principal. Esta es una propiedad muy
interesante que se presenta en las matrices de rigidez
locales y globales.
El método banda de resolución de sistemas de
ecuaciones lineales simultaneas se basa en esta
propiedad, permitiendo el ahorro de asignaciones de
memoria, ya que solo guarda los elementos no nulos.
Materiales y Métodos
Sistema banda de solución
El presente método centra su eficiencia en la forma en que
se realiza la numeración global de los nodos.
S em ian ch o d e b an d a
N ro. total de incognitas
X X X
0 0 0 0
X X X X 0 0 0
X X X X X 0 0
0 X X X X X 0
0 0 X X X X X
0 0 0 X X X X
0 0 0 0 X X X
0 0 0 0 0 X X
0 0 0 0 0 0 X
0 0 0 0 0 0 0
0 0 0 0 0 0 0
0 0 0 0 0 0 0
A n ch o d e b an d a
0
0
0
0
0
0
0
0
0
0
X 0
0
0
0
0
0
0
X X 0
0
0
0
0
0
0
0
X X X 0
0
0
0
0
0
0
0
0
X X X X 0
0 0 X X X
0 X X X X
A n ch o d e b an d a
D iag o n al
p rin cip al
X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X X X
X X X
0 X X X X
0 0 X X X
X X X X
0
X X X 0 0
X X
O rg an izació n co n v en cio n al d e la
m atriz d e rig id ez
0
X 0 0
O rg an izació n en b an d a p ara O rg an izació n en b an d a p ara
m atrices n o sim étricas
m atrices sim étricas
Materiales y Métodos
Método Frontal
Este método trata de evitar el ensamblaje de la matriz
de rigidez global. Consiste básicamente en ensamblar la
menor cantidad de elementos y eliminar las variables al
mismo tiempo, formando matrices de menor tamaño.
Esto es, inmediatamente después de que los coeficientes
de una ecuación se ensamblan a partir de la contribución
de todos los elementos relevantes, o que inciden sobre este
coeficiente, la variable de este puede eliminarse. Como
resultado, la matriz completa de toda la estructura jamás
se forma.
Materiales y Métodos
Método Frontal
A diferencia del método banda, la adecuada
enumeración de los elemento es de gran importancia, no
así la enumeración global de los nodos ya que estos son
tratados por su numeración local. Como primer medida
esta característica del Front lo hace mas versátil ya que
como la cantidad de elementos siempre es menor que la
cantidad de nodos, independientemente de la forma de la
estructura, siempre se tendrá mayos facilidad para
enumerar los elementos que los nodos.
Discusión de Resultados
Para analizar las cualidades de los distintos métodos
propuestos se planteó la resolución mediante el empleo del
FECCUND V 2.1 (software desarrollado por el Dto. De
Mecánica Aplicada de la Facultad de Ingeniería de la
UNNE) de dos ejemplos prácticos, representando los
resultados mediante el pre-post procesador GID 7.2.
Para dicho análisis se tomaró como variables el tiempo
de ejecución de cada método y su exactitud, la cual será
determinada por comparación con los resultados
obtenidos por el software denominado CALSEF provisto
por el Centro Internacional de Métodos Numéricos en
Ingeniería (CIMNE).
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Mallado con elementos
hexaédricos
Desplazamiento nodal
según eje y
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Tabla Nº 1
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
1,758
0,501953E-04
0,47619E-04
5,13
Banda
0,552
0,498874E-04
4,55
Frontal
0,441
0,479135E-04
0,61
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Mallado con elementos hexaédricos
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Desplazamiento nodal según eje y
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Tabla Nº 2
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
31,547
-0,00615923
-0.00544118
11,66
Banda
2,025
-0,00560842
2,98
Frontal
1,752
-0,00552264
1,47
Conclusiones
El método Frontal es superior al Gauss y al Banda
en cuanto a su exactitud y a la velocidad de
ejecución.
El efecto de la propagación de errores por redondeo
incrementa notablemente la deficiencia del método
de Gauss frente a problemas de mayor
envergadura.
El método de Banda puede ser recomendable en el
caso de contar con un software de pre-post
procesamiento ya que los mismos generan las
mallas enumerando los nodos en la forma mas
adecuada disminuyendo de ancho de banda.
Slide 7
Caracterización cualitativa de
algunos métodos de solución de
sistemas de ecuaciones lineales.
MROGINSKI, Javier L. – BENEYTO, Pablo A. –
DI RADO, H. Ariel – MANZOLILLO, Juan E.
Departamento de Mecánica Aplicada.
Facultad de Ingeniería. U.N.N.E.
Antecedentes
El Método de los Elementos Finitos (M.E.F.) es una
técnica de análisis numérico que permite la obtención de
soluciones aproximadas de ecuaciones diferenciales capaces
de representar una variedad de problemas de la Ingeniería y
Física.
El resultado de la aplicación del Método de los Elementos
Finitos a un dominio cualquier es un sistema de ecuaciones
lineales simultáneas de gran tamaño, correspondiendo una
ecuación por cada incógnita que se evalúe en cada nodo,
imposible de ser resuelto sin el auxilio de un computador.
El presente trabajo desarrolla algunos aspectos
significativos sobre los tres diferentes métodos de resolución
de sistemas de ecuaciones lineales simultáneas mas usados en
el Método de los Elementos Finitos que son el: Eliminación
de Gauss, Sistema Banda y Frontal entre otros.
Materiales y Métodos
En las diapositivas siguientes
expondrá una breve descripción de
métodos de solución de sistemas
ecuaciones lineales simultáneas que
estudian en el presente trabajo.
se
los
de
se
Materiales y Métodos
Eliminación de Gauss
Un sistema de ecuaciones lineales simultáneas
puede ser expresado en forma genérica de la
siguiente manera:
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
a 21 x 1 a 22 x 2 a 23 x 3 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... .
a n1 x 1 a n 2 x 2 a n 3 x 3 a n n x n C n
Materiales y Métodos
Eliminación de Gauss
Este método busca obtener un sistema de ecuaciones
triangular equivalente, mediante el uso de transformaciones
lineales. El cual se resuelve con facilidad aplicando la
denominada sustitución inversa.
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
´
´
´
´
a 22 x 2 a 23 x 2 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... ........
(n2)
(n2)
(n2)
a n 1 , n 1 x n 1 a n 1 , n x n C n 1
( n 1 )
( n 1 )
a n,n x n C n
Materiales y Métodos
Sistema banda de solución
Se dice que una matriz es matriz banda cuando posee
todos sus elementos no nulos concentrados en torno de
la diagonal principal. Esta es una propiedad muy
interesante que se presenta en las matrices de rigidez
locales y globales.
El método banda de resolución de sistemas de
ecuaciones lineales simultaneas se basa en esta
propiedad, permitiendo el ahorro de asignaciones de
memoria, ya que solo guarda los elementos no nulos.
Materiales y Métodos
Sistema banda de solución
El presente método centra su eficiencia en la forma en que
se realiza la numeración global de los nodos.
S em ian ch o d e b an d a
N ro. total de incognitas
X X X
0 0 0 0
X X X X 0 0 0
X X X X X 0 0
0 X X X X X 0
0 0 X X X X X
0 0 0 X X X X
0 0 0 0 X X X
0 0 0 0 0 X X
0 0 0 0 0 0 X
0 0 0 0 0 0 0
0 0 0 0 0 0 0
0 0 0 0 0 0 0
A n ch o d e b an d a
0
0
0
0
0
0
0
0
0
0
X 0
0
0
0
0
0
0
X X 0
0
0
0
0
0
0
0
X X X 0
0
0
0
0
0
0
0
0
X X X X 0
0 0 X X X
0 X X X X
A n ch o d e b an d a
D iag o n al
p rin cip al
X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X X X
X X X
0 X X X X
0 0 X X X
X X X X
0
X X X 0 0
X X
O rg an izació n co n v en cio n al d e la
m atriz d e rig id ez
0
X 0 0
O rg an izació n en b an d a p ara O rg an izació n en b an d a p ara
m atrices n o sim étricas
m atrices sim étricas
Materiales y Métodos
Método Frontal
Este método trata de evitar el ensamblaje de la matriz
de rigidez global. Consiste básicamente en ensamblar la
menor cantidad de elementos y eliminar las variables al
mismo tiempo, formando matrices de menor tamaño.
Esto es, inmediatamente después de que los coeficientes
de una ecuación se ensamblan a partir de la contribución
de todos los elementos relevantes, o que inciden sobre este
coeficiente, la variable de este puede eliminarse. Como
resultado, la matriz completa de toda la estructura jamás
se forma.
Materiales y Métodos
Método Frontal
A diferencia del método banda, la adecuada
enumeración de los elemento es de gran importancia, no
así la enumeración global de los nodos ya que estos son
tratados por su numeración local. Como primer medida
esta característica del Front lo hace mas versátil ya que
como la cantidad de elementos siempre es menor que la
cantidad de nodos, independientemente de la forma de la
estructura, siempre se tendrá mayos facilidad para
enumerar los elementos que los nodos.
Discusión de Resultados
Para analizar las cualidades de los distintos métodos
propuestos se planteó la resolución mediante el empleo del
FECCUND V 2.1 (software desarrollado por el Dto. De
Mecánica Aplicada de la Facultad de Ingeniería de la
UNNE) de dos ejemplos prácticos, representando los
resultados mediante el pre-post procesador GID 7.2.
Para dicho análisis se tomaró como variables el tiempo
de ejecución de cada método y su exactitud, la cual será
determinada por comparación con los resultados
obtenidos por el software denominado CALSEF provisto
por el Centro Internacional de Métodos Numéricos en
Ingeniería (CIMNE).
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Mallado con elementos
hexaédricos
Desplazamiento nodal
según eje y
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Tabla Nº 1
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
1,758
0,501953E-04
0,47619E-04
5,13
Banda
0,552
0,498874E-04
4,55
Frontal
0,441
0,479135E-04
0,61
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Mallado con elementos hexaédricos
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Desplazamiento nodal según eje y
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Tabla Nº 2
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
31,547
-0,00615923
-0.00544118
11,66
Banda
2,025
-0,00560842
2,98
Frontal
1,752
-0,00552264
1,47
Conclusiones
El método Frontal es superior al Gauss y al Banda
en cuanto a su exactitud y a la velocidad de
ejecución.
El efecto de la propagación de errores por redondeo
incrementa notablemente la deficiencia del método
de Gauss frente a problemas de mayor
envergadura.
El método de Banda puede ser recomendable en el
caso de contar con un software de pre-post
procesamiento ya que los mismos generan las
mallas enumerando los nodos en la forma mas
adecuada disminuyendo de ancho de banda.
Slide 8
Caracterización cualitativa de
algunos métodos de solución de
sistemas de ecuaciones lineales.
MROGINSKI, Javier L. – BENEYTO, Pablo A. –
DI RADO, H. Ariel – MANZOLILLO, Juan E.
Departamento de Mecánica Aplicada.
Facultad de Ingeniería. U.N.N.E.
Antecedentes
El Método de los Elementos Finitos (M.E.F.) es una
técnica de análisis numérico que permite la obtención de
soluciones aproximadas de ecuaciones diferenciales capaces
de representar una variedad de problemas de la Ingeniería y
Física.
El resultado de la aplicación del Método de los Elementos
Finitos a un dominio cualquier es un sistema de ecuaciones
lineales simultáneas de gran tamaño, correspondiendo una
ecuación por cada incógnita que se evalúe en cada nodo,
imposible de ser resuelto sin el auxilio de un computador.
El presente trabajo desarrolla algunos aspectos
significativos sobre los tres diferentes métodos de resolución
de sistemas de ecuaciones lineales simultáneas mas usados en
el Método de los Elementos Finitos que son el: Eliminación
de Gauss, Sistema Banda y Frontal entre otros.
Materiales y Métodos
En las diapositivas siguientes
expondrá una breve descripción de
métodos de solución de sistemas
ecuaciones lineales simultáneas que
estudian en el presente trabajo.
se
los
de
se
Materiales y Métodos
Eliminación de Gauss
Un sistema de ecuaciones lineales simultáneas
puede ser expresado en forma genérica de la
siguiente manera:
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
a 21 x 1 a 22 x 2 a 23 x 3 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... .
a n1 x 1 a n 2 x 2 a n 3 x 3 a n n x n C n
Materiales y Métodos
Eliminación de Gauss
Este método busca obtener un sistema de ecuaciones
triangular equivalente, mediante el uso de transformaciones
lineales. El cual se resuelve con facilidad aplicando la
denominada sustitución inversa.
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
´
´
´
´
a 22 x 2 a 23 x 2 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... ........
(n2)
(n2)
(n2)
a n 1 , n 1 x n 1 a n 1 , n x n C n 1
( n 1 )
( n 1 )
a n,n x n C n
Materiales y Métodos
Sistema banda de solución
Se dice que una matriz es matriz banda cuando posee
todos sus elementos no nulos concentrados en torno de
la diagonal principal. Esta es una propiedad muy
interesante que se presenta en las matrices de rigidez
locales y globales.
El método banda de resolución de sistemas de
ecuaciones lineales simultaneas se basa en esta
propiedad, permitiendo el ahorro de asignaciones de
memoria, ya que solo guarda los elementos no nulos.
Materiales y Métodos
Sistema banda de solución
El presente método centra su eficiencia en la forma en que
se realiza la numeración global de los nodos.
S em ian ch o d e b an d a
N ro. total de incognitas
X X X
0 0 0 0
X X X X 0 0 0
X X X X X 0 0
0 X X X X X 0
0 0 X X X X X
0 0 0 X X X X
0 0 0 0 X X X
0 0 0 0 0 X X
0 0 0 0 0 0 X
0 0 0 0 0 0 0
0 0 0 0 0 0 0
0 0 0 0 0 0 0
A n ch o d e b an d a
0
0
0
0
0
0
0
0
0
0
X 0
0
0
0
0
0
0
X X 0
0
0
0
0
0
0
0
X X X 0
0
0
0
0
0
0
0
0
X X X X 0
0 0 X X X
0 X X X X
A n ch o d e b an d a
D iag o n al
p rin cip al
X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X X X
X X X
0 X X X X
0 0 X X X
X X X X
0
X X X 0 0
X X
O rg an izació n co n v en cio n al d e la
m atriz d e rig id ez
0
X 0 0
O rg an izació n en b an d a p ara O rg an izació n en b an d a p ara
m atrices n o sim étricas
m atrices sim étricas
Materiales y Métodos
Método Frontal
Este método trata de evitar el ensamblaje de la matriz
de rigidez global. Consiste básicamente en ensamblar la
menor cantidad de elementos y eliminar las variables al
mismo tiempo, formando matrices de menor tamaño.
Esto es, inmediatamente después de que los coeficientes
de una ecuación se ensamblan a partir de la contribución
de todos los elementos relevantes, o que inciden sobre este
coeficiente, la variable de este puede eliminarse. Como
resultado, la matriz completa de toda la estructura jamás
se forma.
Materiales y Métodos
Método Frontal
A diferencia del método banda, la adecuada
enumeración de los elemento es de gran importancia, no
así la enumeración global de los nodos ya que estos son
tratados por su numeración local. Como primer medida
esta característica del Front lo hace mas versátil ya que
como la cantidad de elementos siempre es menor que la
cantidad de nodos, independientemente de la forma de la
estructura, siempre se tendrá mayos facilidad para
enumerar los elementos que los nodos.
Discusión de Resultados
Para analizar las cualidades de los distintos métodos
propuestos se planteó la resolución mediante el empleo del
FECCUND V 2.1 (software desarrollado por el Dto. De
Mecánica Aplicada de la Facultad de Ingeniería de la
UNNE) de dos ejemplos prácticos, representando los
resultados mediante el pre-post procesador GID 7.2.
Para dicho análisis se tomaró como variables el tiempo
de ejecución de cada método y su exactitud, la cual será
determinada por comparación con los resultados
obtenidos por el software denominado CALSEF provisto
por el Centro Internacional de Métodos Numéricos en
Ingeniería (CIMNE).
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Mallado con elementos
hexaédricos
Desplazamiento nodal
según eje y
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Tabla Nº 1
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
1,758
0,501953E-04
0,47619E-04
5,13
Banda
0,552
0,498874E-04
4,55
Frontal
0,441
0,479135E-04
0,61
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Mallado con elementos hexaédricos
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Desplazamiento nodal según eje y
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Tabla Nº 2
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
31,547
-0,00615923
-0.00544118
11,66
Banda
2,025
-0,00560842
2,98
Frontal
1,752
-0,00552264
1,47
Conclusiones
El método Frontal es superior al Gauss y al Banda
en cuanto a su exactitud y a la velocidad de
ejecución.
El efecto de la propagación de errores por redondeo
incrementa notablemente la deficiencia del método
de Gauss frente a problemas de mayor
envergadura.
El método de Banda puede ser recomendable en el
caso de contar con un software de pre-post
procesamiento ya que los mismos generan las
mallas enumerando los nodos en la forma mas
adecuada disminuyendo de ancho de banda.
Slide 9
Caracterización cualitativa de
algunos métodos de solución de
sistemas de ecuaciones lineales.
MROGINSKI, Javier L. – BENEYTO, Pablo A. –
DI RADO, H. Ariel – MANZOLILLO, Juan E.
Departamento de Mecánica Aplicada.
Facultad de Ingeniería. U.N.N.E.
Antecedentes
El Método de los Elementos Finitos (M.E.F.) es una
técnica de análisis numérico que permite la obtención de
soluciones aproximadas de ecuaciones diferenciales capaces
de representar una variedad de problemas de la Ingeniería y
Física.
El resultado de la aplicación del Método de los Elementos
Finitos a un dominio cualquier es un sistema de ecuaciones
lineales simultáneas de gran tamaño, correspondiendo una
ecuación por cada incógnita que se evalúe en cada nodo,
imposible de ser resuelto sin el auxilio de un computador.
El presente trabajo desarrolla algunos aspectos
significativos sobre los tres diferentes métodos de resolución
de sistemas de ecuaciones lineales simultáneas mas usados en
el Método de los Elementos Finitos que son el: Eliminación
de Gauss, Sistema Banda y Frontal entre otros.
Materiales y Métodos
En las diapositivas siguientes
expondrá una breve descripción de
métodos de solución de sistemas
ecuaciones lineales simultáneas que
estudian en el presente trabajo.
se
los
de
se
Materiales y Métodos
Eliminación de Gauss
Un sistema de ecuaciones lineales simultáneas
puede ser expresado en forma genérica de la
siguiente manera:
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
a 21 x 1 a 22 x 2 a 23 x 3 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... .
a n1 x 1 a n 2 x 2 a n 3 x 3 a n n x n C n
Materiales y Métodos
Eliminación de Gauss
Este método busca obtener un sistema de ecuaciones
triangular equivalente, mediante el uso de transformaciones
lineales. El cual se resuelve con facilidad aplicando la
denominada sustitución inversa.
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
´
´
´
´
a 22 x 2 a 23 x 2 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... ........
(n2)
(n2)
(n2)
a n 1 , n 1 x n 1 a n 1 , n x n C n 1
( n 1 )
( n 1 )
a n,n x n C n
Materiales y Métodos
Sistema banda de solución
Se dice que una matriz es matriz banda cuando posee
todos sus elementos no nulos concentrados en torno de
la diagonal principal. Esta es una propiedad muy
interesante que se presenta en las matrices de rigidez
locales y globales.
El método banda de resolución de sistemas de
ecuaciones lineales simultaneas se basa en esta
propiedad, permitiendo el ahorro de asignaciones de
memoria, ya que solo guarda los elementos no nulos.
Materiales y Métodos
Sistema banda de solución
El presente método centra su eficiencia en la forma en que
se realiza la numeración global de los nodos.
S em ian ch o d e b an d a
N ro. total de incognitas
X X X
0 0 0 0
X X X X 0 0 0
X X X X X 0 0
0 X X X X X 0
0 0 X X X X X
0 0 0 X X X X
0 0 0 0 X X X
0 0 0 0 0 X X
0 0 0 0 0 0 X
0 0 0 0 0 0 0
0 0 0 0 0 0 0
0 0 0 0 0 0 0
A n ch o d e b an d a
0
0
0
0
0
0
0
0
0
0
X 0
0
0
0
0
0
0
X X 0
0
0
0
0
0
0
0
X X X 0
0
0
0
0
0
0
0
0
X X X X 0
0 0 X X X
0 X X X X
A n ch o d e b an d a
D iag o n al
p rin cip al
X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X X X
X X X
0 X X X X
0 0 X X X
X X X X
0
X X X 0 0
X X
O rg an izació n co n v en cio n al d e la
m atriz d e rig id ez
0
X 0 0
O rg an izació n en b an d a p ara O rg an izació n en b an d a p ara
m atrices n o sim étricas
m atrices sim étricas
Materiales y Métodos
Método Frontal
Este método trata de evitar el ensamblaje de la matriz
de rigidez global. Consiste básicamente en ensamblar la
menor cantidad de elementos y eliminar las variables al
mismo tiempo, formando matrices de menor tamaño.
Esto es, inmediatamente después de que los coeficientes
de una ecuación se ensamblan a partir de la contribución
de todos los elementos relevantes, o que inciden sobre este
coeficiente, la variable de este puede eliminarse. Como
resultado, la matriz completa de toda la estructura jamás
se forma.
Materiales y Métodos
Método Frontal
A diferencia del método banda, la adecuada
enumeración de los elemento es de gran importancia, no
así la enumeración global de los nodos ya que estos son
tratados por su numeración local. Como primer medida
esta característica del Front lo hace mas versátil ya que
como la cantidad de elementos siempre es menor que la
cantidad de nodos, independientemente de la forma de la
estructura, siempre se tendrá mayos facilidad para
enumerar los elementos que los nodos.
Discusión de Resultados
Para analizar las cualidades de los distintos métodos
propuestos se planteó la resolución mediante el empleo del
FECCUND V 2.1 (software desarrollado por el Dto. De
Mecánica Aplicada de la Facultad de Ingeniería de la
UNNE) de dos ejemplos prácticos, representando los
resultados mediante el pre-post procesador GID 7.2.
Para dicho análisis se tomaró como variables el tiempo
de ejecución de cada método y su exactitud, la cual será
determinada por comparación con los resultados
obtenidos por el software denominado CALSEF provisto
por el Centro Internacional de Métodos Numéricos en
Ingeniería (CIMNE).
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Mallado con elementos
hexaédricos
Desplazamiento nodal
según eje y
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Tabla Nº 1
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
1,758
0,501953E-04
0,47619E-04
5,13
Banda
0,552
0,498874E-04
4,55
Frontal
0,441
0,479135E-04
0,61
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Mallado con elementos hexaédricos
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Desplazamiento nodal según eje y
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Tabla Nº 2
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
31,547
-0,00615923
-0.00544118
11,66
Banda
2,025
-0,00560842
2,98
Frontal
1,752
-0,00552264
1,47
Conclusiones
El método Frontal es superior al Gauss y al Banda
en cuanto a su exactitud y a la velocidad de
ejecución.
El efecto de la propagación de errores por redondeo
incrementa notablemente la deficiencia del método
de Gauss frente a problemas de mayor
envergadura.
El método de Banda puede ser recomendable en el
caso de contar con un software de pre-post
procesamiento ya que los mismos generan las
mallas enumerando los nodos en la forma mas
adecuada disminuyendo de ancho de banda.
Slide 10
Caracterización cualitativa de
algunos métodos de solución de
sistemas de ecuaciones lineales.
MROGINSKI, Javier L. – BENEYTO, Pablo A. –
DI RADO, H. Ariel – MANZOLILLO, Juan E.
Departamento de Mecánica Aplicada.
Facultad de Ingeniería. U.N.N.E.
Antecedentes
El Método de los Elementos Finitos (M.E.F.) es una
técnica de análisis numérico que permite la obtención de
soluciones aproximadas de ecuaciones diferenciales capaces
de representar una variedad de problemas de la Ingeniería y
Física.
El resultado de la aplicación del Método de los Elementos
Finitos a un dominio cualquier es un sistema de ecuaciones
lineales simultáneas de gran tamaño, correspondiendo una
ecuación por cada incógnita que se evalúe en cada nodo,
imposible de ser resuelto sin el auxilio de un computador.
El presente trabajo desarrolla algunos aspectos
significativos sobre los tres diferentes métodos de resolución
de sistemas de ecuaciones lineales simultáneas mas usados en
el Método de los Elementos Finitos que son el: Eliminación
de Gauss, Sistema Banda y Frontal entre otros.
Materiales y Métodos
En las diapositivas siguientes
expondrá una breve descripción de
métodos de solución de sistemas
ecuaciones lineales simultáneas que
estudian en el presente trabajo.
se
los
de
se
Materiales y Métodos
Eliminación de Gauss
Un sistema de ecuaciones lineales simultáneas
puede ser expresado en forma genérica de la
siguiente manera:
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
a 21 x 1 a 22 x 2 a 23 x 3 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... .
a n1 x 1 a n 2 x 2 a n 3 x 3 a n n x n C n
Materiales y Métodos
Eliminación de Gauss
Este método busca obtener un sistema de ecuaciones
triangular equivalente, mediante el uso de transformaciones
lineales. El cual se resuelve con facilidad aplicando la
denominada sustitución inversa.
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
´
´
´
´
a 22 x 2 a 23 x 2 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... ........
(n2)
(n2)
(n2)
a n 1 , n 1 x n 1 a n 1 , n x n C n 1
( n 1 )
( n 1 )
a n,n x n C n
Materiales y Métodos
Sistema banda de solución
Se dice que una matriz es matriz banda cuando posee
todos sus elementos no nulos concentrados en torno de
la diagonal principal. Esta es una propiedad muy
interesante que se presenta en las matrices de rigidez
locales y globales.
El método banda de resolución de sistemas de
ecuaciones lineales simultaneas se basa en esta
propiedad, permitiendo el ahorro de asignaciones de
memoria, ya que solo guarda los elementos no nulos.
Materiales y Métodos
Sistema banda de solución
El presente método centra su eficiencia en la forma en que
se realiza la numeración global de los nodos.
S em ian ch o d e b an d a
N ro. total de incognitas
X X X
0 0 0 0
X X X X 0 0 0
X X X X X 0 0
0 X X X X X 0
0 0 X X X X X
0 0 0 X X X X
0 0 0 0 X X X
0 0 0 0 0 X X
0 0 0 0 0 0 X
0 0 0 0 0 0 0
0 0 0 0 0 0 0
0 0 0 0 0 0 0
A n ch o d e b an d a
0
0
0
0
0
0
0
0
0
0
X 0
0
0
0
0
0
0
X X 0
0
0
0
0
0
0
0
X X X 0
0
0
0
0
0
0
0
0
X X X X 0
0 0 X X X
0 X X X X
A n ch o d e b an d a
D iag o n al
p rin cip al
X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X X X
X X X
0 X X X X
0 0 X X X
X X X X
0
X X X 0 0
X X
O rg an izació n co n v en cio n al d e la
m atriz d e rig id ez
0
X 0 0
O rg an izació n en b an d a p ara O rg an izació n en b an d a p ara
m atrices n o sim étricas
m atrices sim étricas
Materiales y Métodos
Método Frontal
Este método trata de evitar el ensamblaje de la matriz
de rigidez global. Consiste básicamente en ensamblar la
menor cantidad de elementos y eliminar las variables al
mismo tiempo, formando matrices de menor tamaño.
Esto es, inmediatamente después de que los coeficientes
de una ecuación se ensamblan a partir de la contribución
de todos los elementos relevantes, o que inciden sobre este
coeficiente, la variable de este puede eliminarse. Como
resultado, la matriz completa de toda la estructura jamás
se forma.
Materiales y Métodos
Método Frontal
A diferencia del método banda, la adecuada
enumeración de los elemento es de gran importancia, no
así la enumeración global de los nodos ya que estos son
tratados por su numeración local. Como primer medida
esta característica del Front lo hace mas versátil ya que
como la cantidad de elementos siempre es menor que la
cantidad de nodos, independientemente de la forma de la
estructura, siempre se tendrá mayos facilidad para
enumerar los elementos que los nodos.
Discusión de Resultados
Para analizar las cualidades de los distintos métodos
propuestos se planteó la resolución mediante el empleo del
FECCUND V 2.1 (software desarrollado por el Dto. De
Mecánica Aplicada de la Facultad de Ingeniería de la
UNNE) de dos ejemplos prácticos, representando los
resultados mediante el pre-post procesador GID 7.2.
Para dicho análisis se tomaró como variables el tiempo
de ejecución de cada método y su exactitud, la cual será
determinada por comparación con los resultados
obtenidos por el software denominado CALSEF provisto
por el Centro Internacional de Métodos Numéricos en
Ingeniería (CIMNE).
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Mallado con elementos
hexaédricos
Desplazamiento nodal
según eje y
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Tabla Nº 1
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
1,758
0,501953E-04
0,47619E-04
5,13
Banda
0,552
0,498874E-04
4,55
Frontal
0,441
0,479135E-04
0,61
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Mallado con elementos hexaédricos
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Desplazamiento nodal según eje y
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Tabla Nº 2
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
31,547
-0,00615923
-0.00544118
11,66
Banda
2,025
-0,00560842
2,98
Frontal
1,752
-0,00552264
1,47
Conclusiones
El método Frontal es superior al Gauss y al Banda
en cuanto a su exactitud y a la velocidad de
ejecución.
El efecto de la propagación de errores por redondeo
incrementa notablemente la deficiencia del método
de Gauss frente a problemas de mayor
envergadura.
El método de Banda puede ser recomendable en el
caso de contar con un software de pre-post
procesamiento ya que los mismos generan las
mallas enumerando los nodos en la forma mas
adecuada disminuyendo de ancho de banda.
Slide 11
Caracterización cualitativa de
algunos métodos de solución de
sistemas de ecuaciones lineales.
MROGINSKI, Javier L. – BENEYTO, Pablo A. –
DI RADO, H. Ariel – MANZOLILLO, Juan E.
Departamento de Mecánica Aplicada.
Facultad de Ingeniería. U.N.N.E.
Antecedentes
El Método de los Elementos Finitos (M.E.F.) es una
técnica de análisis numérico que permite la obtención de
soluciones aproximadas de ecuaciones diferenciales capaces
de representar una variedad de problemas de la Ingeniería y
Física.
El resultado de la aplicación del Método de los Elementos
Finitos a un dominio cualquier es un sistema de ecuaciones
lineales simultáneas de gran tamaño, correspondiendo una
ecuación por cada incógnita que se evalúe en cada nodo,
imposible de ser resuelto sin el auxilio de un computador.
El presente trabajo desarrolla algunos aspectos
significativos sobre los tres diferentes métodos de resolución
de sistemas de ecuaciones lineales simultáneas mas usados en
el Método de los Elementos Finitos que son el: Eliminación
de Gauss, Sistema Banda y Frontal entre otros.
Materiales y Métodos
En las diapositivas siguientes
expondrá una breve descripción de
métodos de solución de sistemas
ecuaciones lineales simultáneas que
estudian en el presente trabajo.
se
los
de
se
Materiales y Métodos
Eliminación de Gauss
Un sistema de ecuaciones lineales simultáneas
puede ser expresado en forma genérica de la
siguiente manera:
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
a 21 x 1 a 22 x 2 a 23 x 3 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... .
a n1 x 1 a n 2 x 2 a n 3 x 3 a n n x n C n
Materiales y Métodos
Eliminación de Gauss
Este método busca obtener un sistema de ecuaciones
triangular equivalente, mediante el uso de transformaciones
lineales. El cual se resuelve con facilidad aplicando la
denominada sustitución inversa.
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
´
´
´
´
a 22 x 2 a 23 x 2 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... ........
(n2)
(n2)
(n2)
a n 1 , n 1 x n 1 a n 1 , n x n C n 1
( n 1 )
( n 1 )
a n,n x n C n
Materiales y Métodos
Sistema banda de solución
Se dice que una matriz es matriz banda cuando posee
todos sus elementos no nulos concentrados en torno de
la diagonal principal. Esta es una propiedad muy
interesante que se presenta en las matrices de rigidez
locales y globales.
El método banda de resolución de sistemas de
ecuaciones lineales simultaneas se basa en esta
propiedad, permitiendo el ahorro de asignaciones de
memoria, ya que solo guarda los elementos no nulos.
Materiales y Métodos
Sistema banda de solución
El presente método centra su eficiencia en la forma en que
se realiza la numeración global de los nodos.
S em ian ch o d e b an d a
N ro. total de incognitas
X X X
0 0 0 0
X X X X 0 0 0
X X X X X 0 0
0 X X X X X 0
0 0 X X X X X
0 0 0 X X X X
0 0 0 0 X X X
0 0 0 0 0 X X
0 0 0 0 0 0 X
0 0 0 0 0 0 0
0 0 0 0 0 0 0
0 0 0 0 0 0 0
A n ch o d e b an d a
0
0
0
0
0
0
0
0
0
0
X 0
0
0
0
0
0
0
X X 0
0
0
0
0
0
0
0
X X X 0
0
0
0
0
0
0
0
0
X X X X 0
0 0 X X X
0 X X X X
A n ch o d e b an d a
D iag o n al
p rin cip al
X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X X X
X X X
0 X X X X
0 0 X X X
X X X X
0
X X X 0 0
X X
O rg an izació n co n v en cio n al d e la
m atriz d e rig id ez
0
X 0 0
O rg an izació n en b an d a p ara O rg an izació n en b an d a p ara
m atrices n o sim étricas
m atrices sim étricas
Materiales y Métodos
Método Frontal
Este método trata de evitar el ensamblaje de la matriz
de rigidez global. Consiste básicamente en ensamblar la
menor cantidad de elementos y eliminar las variables al
mismo tiempo, formando matrices de menor tamaño.
Esto es, inmediatamente después de que los coeficientes
de una ecuación se ensamblan a partir de la contribución
de todos los elementos relevantes, o que inciden sobre este
coeficiente, la variable de este puede eliminarse. Como
resultado, la matriz completa de toda la estructura jamás
se forma.
Materiales y Métodos
Método Frontal
A diferencia del método banda, la adecuada
enumeración de los elemento es de gran importancia, no
así la enumeración global de los nodos ya que estos son
tratados por su numeración local. Como primer medida
esta característica del Front lo hace mas versátil ya que
como la cantidad de elementos siempre es menor que la
cantidad de nodos, independientemente de la forma de la
estructura, siempre se tendrá mayos facilidad para
enumerar los elementos que los nodos.
Discusión de Resultados
Para analizar las cualidades de los distintos métodos
propuestos se planteó la resolución mediante el empleo del
FECCUND V 2.1 (software desarrollado por el Dto. De
Mecánica Aplicada de la Facultad de Ingeniería de la
UNNE) de dos ejemplos prácticos, representando los
resultados mediante el pre-post procesador GID 7.2.
Para dicho análisis se tomaró como variables el tiempo
de ejecución de cada método y su exactitud, la cual será
determinada por comparación con los resultados
obtenidos por el software denominado CALSEF provisto
por el Centro Internacional de Métodos Numéricos en
Ingeniería (CIMNE).
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Mallado con elementos
hexaédricos
Desplazamiento nodal
según eje y
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Tabla Nº 1
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
1,758
0,501953E-04
0,47619E-04
5,13
Banda
0,552
0,498874E-04
4,55
Frontal
0,441
0,479135E-04
0,61
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Mallado con elementos hexaédricos
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Desplazamiento nodal según eje y
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Tabla Nº 2
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
31,547
-0,00615923
-0.00544118
11,66
Banda
2,025
-0,00560842
2,98
Frontal
1,752
-0,00552264
1,47
Conclusiones
El método Frontal es superior al Gauss y al Banda
en cuanto a su exactitud y a la velocidad de
ejecución.
El efecto de la propagación de errores por redondeo
incrementa notablemente la deficiencia del método
de Gauss frente a problemas de mayor
envergadura.
El método de Banda puede ser recomendable en el
caso de contar con un software de pre-post
procesamiento ya que los mismos generan las
mallas enumerando los nodos en la forma mas
adecuada disminuyendo de ancho de banda.
Slide 12
Caracterización cualitativa de
algunos métodos de solución de
sistemas de ecuaciones lineales.
MROGINSKI, Javier L. – BENEYTO, Pablo A. –
DI RADO, H. Ariel – MANZOLILLO, Juan E.
Departamento de Mecánica Aplicada.
Facultad de Ingeniería. U.N.N.E.
Antecedentes
El Método de los Elementos Finitos (M.E.F.) es una
técnica de análisis numérico que permite la obtención de
soluciones aproximadas de ecuaciones diferenciales capaces
de representar una variedad de problemas de la Ingeniería y
Física.
El resultado de la aplicación del Método de los Elementos
Finitos a un dominio cualquier es un sistema de ecuaciones
lineales simultáneas de gran tamaño, correspondiendo una
ecuación por cada incógnita que se evalúe en cada nodo,
imposible de ser resuelto sin el auxilio de un computador.
El presente trabajo desarrolla algunos aspectos
significativos sobre los tres diferentes métodos de resolución
de sistemas de ecuaciones lineales simultáneas mas usados en
el Método de los Elementos Finitos que son el: Eliminación
de Gauss, Sistema Banda y Frontal entre otros.
Materiales y Métodos
En las diapositivas siguientes
expondrá una breve descripción de
métodos de solución de sistemas
ecuaciones lineales simultáneas que
estudian en el presente trabajo.
se
los
de
se
Materiales y Métodos
Eliminación de Gauss
Un sistema de ecuaciones lineales simultáneas
puede ser expresado en forma genérica de la
siguiente manera:
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
a 21 x 1 a 22 x 2 a 23 x 3 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... .
a n1 x 1 a n 2 x 2 a n 3 x 3 a n n x n C n
Materiales y Métodos
Eliminación de Gauss
Este método busca obtener un sistema de ecuaciones
triangular equivalente, mediante el uso de transformaciones
lineales. El cual se resuelve con facilidad aplicando la
denominada sustitución inversa.
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
´
´
´
´
a 22 x 2 a 23 x 2 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... ........
(n2)
(n2)
(n2)
a n 1 , n 1 x n 1 a n 1 , n x n C n 1
( n 1 )
( n 1 )
a n,n x n C n
Materiales y Métodos
Sistema banda de solución
Se dice que una matriz es matriz banda cuando posee
todos sus elementos no nulos concentrados en torno de
la diagonal principal. Esta es una propiedad muy
interesante que se presenta en las matrices de rigidez
locales y globales.
El método banda de resolución de sistemas de
ecuaciones lineales simultaneas se basa en esta
propiedad, permitiendo el ahorro de asignaciones de
memoria, ya que solo guarda los elementos no nulos.
Materiales y Métodos
Sistema banda de solución
El presente método centra su eficiencia en la forma en que
se realiza la numeración global de los nodos.
S em ian ch o d e b an d a
N ro. total de incognitas
X X X
0 0 0 0
X X X X 0 0 0
X X X X X 0 0
0 X X X X X 0
0 0 X X X X X
0 0 0 X X X X
0 0 0 0 X X X
0 0 0 0 0 X X
0 0 0 0 0 0 X
0 0 0 0 0 0 0
0 0 0 0 0 0 0
0 0 0 0 0 0 0
A n ch o d e b an d a
0
0
0
0
0
0
0
0
0
0
X 0
0
0
0
0
0
0
X X 0
0
0
0
0
0
0
0
X X X 0
0
0
0
0
0
0
0
0
X X X X 0
0 0 X X X
0 X X X X
A n ch o d e b an d a
D iag o n al
p rin cip al
X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X X X
X X X
0 X X X X
0 0 X X X
X X X X
0
X X X 0 0
X X
O rg an izació n co n v en cio n al d e la
m atriz d e rig id ez
0
X 0 0
O rg an izació n en b an d a p ara O rg an izació n en b an d a p ara
m atrices n o sim étricas
m atrices sim étricas
Materiales y Métodos
Método Frontal
Este método trata de evitar el ensamblaje de la matriz
de rigidez global. Consiste básicamente en ensamblar la
menor cantidad de elementos y eliminar las variables al
mismo tiempo, formando matrices de menor tamaño.
Esto es, inmediatamente después de que los coeficientes
de una ecuación se ensamblan a partir de la contribución
de todos los elementos relevantes, o que inciden sobre este
coeficiente, la variable de este puede eliminarse. Como
resultado, la matriz completa de toda la estructura jamás
se forma.
Materiales y Métodos
Método Frontal
A diferencia del método banda, la adecuada
enumeración de los elemento es de gran importancia, no
así la enumeración global de los nodos ya que estos son
tratados por su numeración local. Como primer medida
esta característica del Front lo hace mas versátil ya que
como la cantidad de elementos siempre es menor que la
cantidad de nodos, independientemente de la forma de la
estructura, siempre se tendrá mayos facilidad para
enumerar los elementos que los nodos.
Discusión de Resultados
Para analizar las cualidades de los distintos métodos
propuestos se planteó la resolución mediante el empleo del
FECCUND V 2.1 (software desarrollado por el Dto. De
Mecánica Aplicada de la Facultad de Ingeniería de la
UNNE) de dos ejemplos prácticos, representando los
resultados mediante el pre-post procesador GID 7.2.
Para dicho análisis se tomaró como variables el tiempo
de ejecución de cada método y su exactitud, la cual será
determinada por comparación con los resultados
obtenidos por el software denominado CALSEF provisto
por el Centro Internacional de Métodos Numéricos en
Ingeniería (CIMNE).
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Mallado con elementos
hexaédricos
Desplazamiento nodal
según eje y
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Tabla Nº 1
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
1,758
0,501953E-04
0,47619E-04
5,13
Banda
0,552
0,498874E-04
4,55
Frontal
0,441
0,479135E-04
0,61
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Mallado con elementos hexaédricos
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Desplazamiento nodal según eje y
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Tabla Nº 2
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
31,547
-0,00615923
-0.00544118
11,66
Banda
2,025
-0,00560842
2,98
Frontal
1,752
-0,00552264
1,47
Conclusiones
El método Frontal es superior al Gauss y al Banda
en cuanto a su exactitud y a la velocidad de
ejecución.
El efecto de la propagación de errores por redondeo
incrementa notablemente la deficiencia del método
de Gauss frente a problemas de mayor
envergadura.
El método de Banda puede ser recomendable en el
caso de contar con un software de pre-post
procesamiento ya que los mismos generan las
mallas enumerando los nodos en la forma mas
adecuada disminuyendo de ancho de banda.
Slide 13
Caracterización cualitativa de
algunos métodos de solución de
sistemas de ecuaciones lineales.
MROGINSKI, Javier L. – BENEYTO, Pablo A. –
DI RADO, H. Ariel – MANZOLILLO, Juan E.
Departamento de Mecánica Aplicada.
Facultad de Ingeniería. U.N.N.E.
Antecedentes
El Método de los Elementos Finitos (M.E.F.) es una
técnica de análisis numérico que permite la obtención de
soluciones aproximadas de ecuaciones diferenciales capaces
de representar una variedad de problemas de la Ingeniería y
Física.
El resultado de la aplicación del Método de los Elementos
Finitos a un dominio cualquier es un sistema de ecuaciones
lineales simultáneas de gran tamaño, correspondiendo una
ecuación por cada incógnita que se evalúe en cada nodo,
imposible de ser resuelto sin el auxilio de un computador.
El presente trabajo desarrolla algunos aspectos
significativos sobre los tres diferentes métodos de resolución
de sistemas de ecuaciones lineales simultáneas mas usados en
el Método de los Elementos Finitos que son el: Eliminación
de Gauss, Sistema Banda y Frontal entre otros.
Materiales y Métodos
En las diapositivas siguientes
expondrá una breve descripción de
métodos de solución de sistemas
ecuaciones lineales simultáneas que
estudian en el presente trabajo.
se
los
de
se
Materiales y Métodos
Eliminación de Gauss
Un sistema de ecuaciones lineales simultáneas
puede ser expresado en forma genérica de la
siguiente manera:
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
a 21 x 1 a 22 x 2 a 23 x 3 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... .
a n1 x 1 a n 2 x 2 a n 3 x 3 a n n x n C n
Materiales y Métodos
Eliminación de Gauss
Este método busca obtener un sistema de ecuaciones
triangular equivalente, mediante el uso de transformaciones
lineales. El cual se resuelve con facilidad aplicando la
denominada sustitución inversa.
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
´
´
´
´
a 22 x 2 a 23 x 2 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... ........
(n2)
(n2)
(n2)
a n 1 , n 1 x n 1 a n 1 , n x n C n 1
( n 1 )
( n 1 )
a n,n x n C n
Materiales y Métodos
Sistema banda de solución
Se dice que una matriz es matriz banda cuando posee
todos sus elementos no nulos concentrados en torno de
la diagonal principal. Esta es una propiedad muy
interesante que se presenta en las matrices de rigidez
locales y globales.
El método banda de resolución de sistemas de
ecuaciones lineales simultaneas se basa en esta
propiedad, permitiendo el ahorro de asignaciones de
memoria, ya que solo guarda los elementos no nulos.
Materiales y Métodos
Sistema banda de solución
El presente método centra su eficiencia en la forma en que
se realiza la numeración global de los nodos.
S em ian ch o d e b an d a
N ro. total de incognitas
X X X
0 0 0 0
X X X X 0 0 0
X X X X X 0 0
0 X X X X X 0
0 0 X X X X X
0 0 0 X X X X
0 0 0 0 X X X
0 0 0 0 0 X X
0 0 0 0 0 0 X
0 0 0 0 0 0 0
0 0 0 0 0 0 0
0 0 0 0 0 0 0
A n ch o d e b an d a
0
0
0
0
0
0
0
0
0
0
X 0
0
0
0
0
0
0
X X 0
0
0
0
0
0
0
0
X X X 0
0
0
0
0
0
0
0
0
X X X X 0
0 0 X X X
0 X X X X
A n ch o d e b an d a
D iag o n al
p rin cip al
X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X X X
X X X
0 X X X X
0 0 X X X
X X X X
0
X X X 0 0
X X
O rg an izació n co n v en cio n al d e la
m atriz d e rig id ez
0
X 0 0
O rg an izació n en b an d a p ara O rg an izació n en b an d a p ara
m atrices n o sim étricas
m atrices sim étricas
Materiales y Métodos
Método Frontal
Este método trata de evitar el ensamblaje de la matriz
de rigidez global. Consiste básicamente en ensamblar la
menor cantidad de elementos y eliminar las variables al
mismo tiempo, formando matrices de menor tamaño.
Esto es, inmediatamente después de que los coeficientes
de una ecuación se ensamblan a partir de la contribución
de todos los elementos relevantes, o que inciden sobre este
coeficiente, la variable de este puede eliminarse. Como
resultado, la matriz completa de toda la estructura jamás
se forma.
Materiales y Métodos
Método Frontal
A diferencia del método banda, la adecuada
enumeración de los elemento es de gran importancia, no
así la enumeración global de los nodos ya que estos son
tratados por su numeración local. Como primer medida
esta característica del Front lo hace mas versátil ya que
como la cantidad de elementos siempre es menor que la
cantidad de nodos, independientemente de la forma de la
estructura, siempre se tendrá mayos facilidad para
enumerar los elementos que los nodos.
Discusión de Resultados
Para analizar las cualidades de los distintos métodos
propuestos se planteó la resolución mediante el empleo del
FECCUND V 2.1 (software desarrollado por el Dto. De
Mecánica Aplicada de la Facultad de Ingeniería de la
UNNE) de dos ejemplos prácticos, representando los
resultados mediante el pre-post procesador GID 7.2.
Para dicho análisis se tomaró como variables el tiempo
de ejecución de cada método y su exactitud, la cual será
determinada por comparación con los resultados
obtenidos por el software denominado CALSEF provisto
por el Centro Internacional de Métodos Numéricos en
Ingeniería (CIMNE).
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Mallado con elementos
hexaédricos
Desplazamiento nodal
según eje y
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Tabla Nº 1
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
1,758
0,501953E-04
0,47619E-04
5,13
Banda
0,552
0,498874E-04
4,55
Frontal
0,441
0,479135E-04
0,61
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Mallado con elementos hexaédricos
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Desplazamiento nodal según eje y
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Tabla Nº 2
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
31,547
-0,00615923
-0.00544118
11,66
Banda
2,025
-0,00560842
2,98
Frontal
1,752
-0,00552264
1,47
Conclusiones
El método Frontal es superior al Gauss y al Banda
en cuanto a su exactitud y a la velocidad de
ejecución.
El efecto de la propagación de errores por redondeo
incrementa notablemente la deficiencia del método
de Gauss frente a problemas de mayor
envergadura.
El método de Banda puede ser recomendable en el
caso de contar con un software de pre-post
procesamiento ya que los mismos generan las
mallas enumerando los nodos en la forma mas
adecuada disminuyendo de ancho de banda.
Slide 14
Caracterización cualitativa de
algunos métodos de solución de
sistemas de ecuaciones lineales.
MROGINSKI, Javier L. – BENEYTO, Pablo A. –
DI RADO, H. Ariel – MANZOLILLO, Juan E.
Departamento de Mecánica Aplicada.
Facultad de Ingeniería. U.N.N.E.
Antecedentes
El Método de los Elementos Finitos (M.E.F.) es una
técnica de análisis numérico que permite la obtención de
soluciones aproximadas de ecuaciones diferenciales capaces
de representar una variedad de problemas de la Ingeniería y
Física.
El resultado de la aplicación del Método de los Elementos
Finitos a un dominio cualquier es un sistema de ecuaciones
lineales simultáneas de gran tamaño, correspondiendo una
ecuación por cada incógnita que se evalúe en cada nodo,
imposible de ser resuelto sin el auxilio de un computador.
El presente trabajo desarrolla algunos aspectos
significativos sobre los tres diferentes métodos de resolución
de sistemas de ecuaciones lineales simultáneas mas usados en
el Método de los Elementos Finitos que son el: Eliminación
de Gauss, Sistema Banda y Frontal entre otros.
Materiales y Métodos
En las diapositivas siguientes
expondrá una breve descripción de
métodos de solución de sistemas
ecuaciones lineales simultáneas que
estudian en el presente trabajo.
se
los
de
se
Materiales y Métodos
Eliminación de Gauss
Un sistema de ecuaciones lineales simultáneas
puede ser expresado en forma genérica de la
siguiente manera:
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
a 21 x 1 a 22 x 2 a 23 x 3 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... .
a n1 x 1 a n 2 x 2 a n 3 x 3 a n n x n C n
Materiales y Métodos
Eliminación de Gauss
Este método busca obtener un sistema de ecuaciones
triangular equivalente, mediante el uso de transformaciones
lineales. El cual se resuelve con facilidad aplicando la
denominada sustitución inversa.
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
´
´
´
´
a 22 x 2 a 23 x 2 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... ........
(n2)
(n2)
(n2)
a n 1 , n 1 x n 1 a n 1 , n x n C n 1
( n 1 )
( n 1 )
a n,n x n C n
Materiales y Métodos
Sistema banda de solución
Se dice que una matriz es matriz banda cuando posee
todos sus elementos no nulos concentrados en torno de
la diagonal principal. Esta es una propiedad muy
interesante que se presenta en las matrices de rigidez
locales y globales.
El método banda de resolución de sistemas de
ecuaciones lineales simultaneas se basa en esta
propiedad, permitiendo el ahorro de asignaciones de
memoria, ya que solo guarda los elementos no nulos.
Materiales y Métodos
Sistema banda de solución
El presente método centra su eficiencia en la forma en que
se realiza la numeración global de los nodos.
S em ian ch o d e b an d a
N ro. total de incognitas
X X X
0 0 0 0
X X X X 0 0 0
X X X X X 0 0
0 X X X X X 0
0 0 X X X X X
0 0 0 X X X X
0 0 0 0 X X X
0 0 0 0 0 X X
0 0 0 0 0 0 X
0 0 0 0 0 0 0
0 0 0 0 0 0 0
0 0 0 0 0 0 0
A n ch o d e b an d a
0
0
0
0
0
0
0
0
0
0
X 0
0
0
0
0
0
0
X X 0
0
0
0
0
0
0
0
X X X 0
0
0
0
0
0
0
0
0
X X X X 0
0 0 X X X
0 X X X X
A n ch o d e b an d a
D iag o n al
p rin cip al
X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X X X
X X X
0 X X X X
0 0 X X X
X X X X
0
X X X 0 0
X X
O rg an izació n co n v en cio n al d e la
m atriz d e rig id ez
0
X 0 0
O rg an izació n en b an d a p ara O rg an izació n en b an d a p ara
m atrices n o sim étricas
m atrices sim étricas
Materiales y Métodos
Método Frontal
Este método trata de evitar el ensamblaje de la matriz
de rigidez global. Consiste básicamente en ensamblar la
menor cantidad de elementos y eliminar las variables al
mismo tiempo, formando matrices de menor tamaño.
Esto es, inmediatamente después de que los coeficientes
de una ecuación se ensamblan a partir de la contribución
de todos los elementos relevantes, o que inciden sobre este
coeficiente, la variable de este puede eliminarse. Como
resultado, la matriz completa de toda la estructura jamás
se forma.
Materiales y Métodos
Método Frontal
A diferencia del método banda, la adecuada
enumeración de los elemento es de gran importancia, no
así la enumeración global de los nodos ya que estos son
tratados por su numeración local. Como primer medida
esta característica del Front lo hace mas versátil ya que
como la cantidad de elementos siempre es menor que la
cantidad de nodos, independientemente de la forma de la
estructura, siempre se tendrá mayos facilidad para
enumerar los elementos que los nodos.
Discusión de Resultados
Para analizar las cualidades de los distintos métodos
propuestos se planteó la resolución mediante el empleo del
FECCUND V 2.1 (software desarrollado por el Dto. De
Mecánica Aplicada de la Facultad de Ingeniería de la
UNNE) de dos ejemplos prácticos, representando los
resultados mediante el pre-post procesador GID 7.2.
Para dicho análisis se tomaró como variables el tiempo
de ejecución de cada método y su exactitud, la cual será
determinada por comparación con los resultados
obtenidos por el software denominado CALSEF provisto
por el Centro Internacional de Métodos Numéricos en
Ingeniería (CIMNE).
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Mallado con elementos
hexaédricos
Desplazamiento nodal
según eje y
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Tabla Nº 1
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
1,758
0,501953E-04
0,47619E-04
5,13
Banda
0,552
0,498874E-04
4,55
Frontal
0,441
0,479135E-04
0,61
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Mallado con elementos hexaédricos
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Desplazamiento nodal según eje y
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Tabla Nº 2
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
31,547
-0,00615923
-0.00544118
11,66
Banda
2,025
-0,00560842
2,98
Frontal
1,752
-0,00552264
1,47
Conclusiones
El método Frontal es superior al Gauss y al Banda
en cuanto a su exactitud y a la velocidad de
ejecución.
El efecto de la propagación de errores por redondeo
incrementa notablemente la deficiencia del método
de Gauss frente a problemas de mayor
envergadura.
El método de Banda puede ser recomendable en el
caso de contar con un software de pre-post
procesamiento ya que los mismos generan las
mallas enumerando los nodos en la forma mas
adecuada disminuyendo de ancho de banda.
Slide 15
Caracterización cualitativa de
algunos métodos de solución de
sistemas de ecuaciones lineales.
MROGINSKI, Javier L. – BENEYTO, Pablo A. –
DI RADO, H. Ariel – MANZOLILLO, Juan E.
Departamento de Mecánica Aplicada.
Facultad de Ingeniería. U.N.N.E.
Antecedentes
El Método de los Elementos Finitos (M.E.F.) es una
técnica de análisis numérico que permite la obtención de
soluciones aproximadas de ecuaciones diferenciales capaces
de representar una variedad de problemas de la Ingeniería y
Física.
El resultado de la aplicación del Método de los Elementos
Finitos a un dominio cualquier es un sistema de ecuaciones
lineales simultáneas de gran tamaño, correspondiendo una
ecuación por cada incógnita que se evalúe en cada nodo,
imposible de ser resuelto sin el auxilio de un computador.
El presente trabajo desarrolla algunos aspectos
significativos sobre los tres diferentes métodos de resolución
de sistemas de ecuaciones lineales simultáneas mas usados en
el Método de los Elementos Finitos que son el: Eliminación
de Gauss, Sistema Banda y Frontal entre otros.
Materiales y Métodos
En las diapositivas siguientes
expondrá una breve descripción de
métodos de solución de sistemas
ecuaciones lineales simultáneas que
estudian en el presente trabajo.
se
los
de
se
Materiales y Métodos
Eliminación de Gauss
Un sistema de ecuaciones lineales simultáneas
puede ser expresado en forma genérica de la
siguiente manera:
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
a 21 x 1 a 22 x 2 a 23 x 3 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... .
a n1 x 1 a n 2 x 2 a n 3 x 3 a n n x n C n
Materiales y Métodos
Eliminación de Gauss
Este método busca obtener un sistema de ecuaciones
triangular equivalente, mediante el uso de transformaciones
lineales. El cual se resuelve con facilidad aplicando la
denominada sustitución inversa.
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
´
´
´
´
a 22 x 2 a 23 x 2 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... ........
(n2)
(n2)
(n2)
a n 1 , n 1 x n 1 a n 1 , n x n C n 1
( n 1 )
( n 1 )
a n,n x n C n
Materiales y Métodos
Sistema banda de solución
Se dice que una matriz es matriz banda cuando posee
todos sus elementos no nulos concentrados en torno de
la diagonal principal. Esta es una propiedad muy
interesante que se presenta en las matrices de rigidez
locales y globales.
El método banda de resolución de sistemas de
ecuaciones lineales simultaneas se basa en esta
propiedad, permitiendo el ahorro de asignaciones de
memoria, ya que solo guarda los elementos no nulos.
Materiales y Métodos
Sistema banda de solución
El presente método centra su eficiencia en la forma en que
se realiza la numeración global de los nodos.
S em ian ch o d e b an d a
N ro. total de incognitas
X X X
0 0 0 0
X X X X 0 0 0
X X X X X 0 0
0 X X X X X 0
0 0 X X X X X
0 0 0 X X X X
0 0 0 0 X X X
0 0 0 0 0 X X
0 0 0 0 0 0 X
0 0 0 0 0 0 0
0 0 0 0 0 0 0
0 0 0 0 0 0 0
A n ch o d e b an d a
0
0
0
0
0
0
0
0
0
0
X 0
0
0
0
0
0
0
X X 0
0
0
0
0
0
0
0
X X X 0
0
0
0
0
0
0
0
0
X X X X 0
0 0 X X X
0 X X X X
A n ch o d e b an d a
D iag o n al
p rin cip al
X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X X X
X X X
0 X X X X
0 0 X X X
X X X X
0
X X X 0 0
X X
O rg an izació n co n v en cio n al d e la
m atriz d e rig id ez
0
X 0 0
O rg an izació n en b an d a p ara O rg an izació n en b an d a p ara
m atrices n o sim étricas
m atrices sim étricas
Materiales y Métodos
Método Frontal
Este método trata de evitar el ensamblaje de la matriz
de rigidez global. Consiste básicamente en ensamblar la
menor cantidad de elementos y eliminar las variables al
mismo tiempo, formando matrices de menor tamaño.
Esto es, inmediatamente después de que los coeficientes
de una ecuación se ensamblan a partir de la contribución
de todos los elementos relevantes, o que inciden sobre este
coeficiente, la variable de este puede eliminarse. Como
resultado, la matriz completa de toda la estructura jamás
se forma.
Materiales y Métodos
Método Frontal
A diferencia del método banda, la adecuada
enumeración de los elemento es de gran importancia, no
así la enumeración global de los nodos ya que estos son
tratados por su numeración local. Como primer medida
esta característica del Front lo hace mas versátil ya que
como la cantidad de elementos siempre es menor que la
cantidad de nodos, independientemente de la forma de la
estructura, siempre se tendrá mayos facilidad para
enumerar los elementos que los nodos.
Discusión de Resultados
Para analizar las cualidades de los distintos métodos
propuestos se planteó la resolución mediante el empleo del
FECCUND V 2.1 (software desarrollado por el Dto. De
Mecánica Aplicada de la Facultad de Ingeniería de la
UNNE) de dos ejemplos prácticos, representando los
resultados mediante el pre-post procesador GID 7.2.
Para dicho análisis se tomaró como variables el tiempo
de ejecución de cada método y su exactitud, la cual será
determinada por comparación con los resultados
obtenidos por el software denominado CALSEF provisto
por el Centro Internacional de Métodos Numéricos en
Ingeniería (CIMNE).
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Mallado con elementos
hexaédricos
Desplazamiento nodal
según eje y
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Tabla Nº 1
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
1,758
0,501953E-04
0,47619E-04
5,13
Banda
0,552
0,498874E-04
4,55
Frontal
0,441
0,479135E-04
0,61
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Mallado con elementos hexaédricos
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Desplazamiento nodal según eje y
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Tabla Nº 2
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
31,547
-0,00615923
-0.00544118
11,66
Banda
2,025
-0,00560842
2,98
Frontal
1,752
-0,00552264
1,47
Conclusiones
El método Frontal es superior al Gauss y al Banda
en cuanto a su exactitud y a la velocidad de
ejecución.
El efecto de la propagación de errores por redondeo
incrementa notablemente la deficiencia del método
de Gauss frente a problemas de mayor
envergadura.
El método de Banda puede ser recomendable en el
caso de contar con un software de pre-post
procesamiento ya que los mismos generan las
mallas enumerando los nodos en la forma mas
adecuada disminuyendo de ancho de banda.
Slide 16
Caracterización cualitativa de
algunos métodos de solución de
sistemas de ecuaciones lineales.
MROGINSKI, Javier L. – BENEYTO, Pablo A. –
DI RADO, H. Ariel – MANZOLILLO, Juan E.
Departamento de Mecánica Aplicada.
Facultad de Ingeniería. U.N.N.E.
Antecedentes
El Método de los Elementos Finitos (M.E.F.) es una
técnica de análisis numérico que permite la obtención de
soluciones aproximadas de ecuaciones diferenciales capaces
de representar una variedad de problemas de la Ingeniería y
Física.
El resultado de la aplicación del Método de los Elementos
Finitos a un dominio cualquier es un sistema de ecuaciones
lineales simultáneas de gran tamaño, correspondiendo una
ecuación por cada incógnita que se evalúe en cada nodo,
imposible de ser resuelto sin el auxilio de un computador.
El presente trabajo desarrolla algunos aspectos
significativos sobre los tres diferentes métodos de resolución
de sistemas de ecuaciones lineales simultáneas mas usados en
el Método de los Elementos Finitos que son el: Eliminación
de Gauss, Sistema Banda y Frontal entre otros.
Materiales y Métodos
En las diapositivas siguientes
expondrá una breve descripción de
métodos de solución de sistemas
ecuaciones lineales simultáneas que
estudian en el presente trabajo.
se
los
de
se
Materiales y Métodos
Eliminación de Gauss
Un sistema de ecuaciones lineales simultáneas
puede ser expresado en forma genérica de la
siguiente manera:
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
a 21 x 1 a 22 x 2 a 23 x 3 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... .
a n1 x 1 a n 2 x 2 a n 3 x 3 a n n x n C n
Materiales y Métodos
Eliminación de Gauss
Este método busca obtener un sistema de ecuaciones
triangular equivalente, mediante el uso de transformaciones
lineales. El cual se resuelve con facilidad aplicando la
denominada sustitución inversa.
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
´
´
´
´
a 22 x 2 a 23 x 2 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... ........
(n2)
(n2)
(n2)
a n 1 , n 1 x n 1 a n 1 , n x n C n 1
( n 1 )
( n 1 )
a n,n x n C n
Materiales y Métodos
Sistema banda de solución
Se dice que una matriz es matriz banda cuando posee
todos sus elementos no nulos concentrados en torno de
la diagonal principal. Esta es una propiedad muy
interesante que se presenta en las matrices de rigidez
locales y globales.
El método banda de resolución de sistemas de
ecuaciones lineales simultaneas se basa en esta
propiedad, permitiendo el ahorro de asignaciones de
memoria, ya que solo guarda los elementos no nulos.
Materiales y Métodos
Sistema banda de solución
El presente método centra su eficiencia en la forma en que
se realiza la numeración global de los nodos.
S em ian ch o d e b an d a
N ro. total de incognitas
X X X
0 0 0 0
X X X X 0 0 0
X X X X X 0 0
0 X X X X X 0
0 0 X X X X X
0 0 0 X X X X
0 0 0 0 X X X
0 0 0 0 0 X X
0 0 0 0 0 0 X
0 0 0 0 0 0 0
0 0 0 0 0 0 0
0 0 0 0 0 0 0
A n ch o d e b an d a
0
0
0
0
0
0
0
0
0
0
X 0
0
0
0
0
0
0
X X 0
0
0
0
0
0
0
0
X X X 0
0
0
0
0
0
0
0
0
X X X X 0
0 0 X X X
0 X X X X
A n ch o d e b an d a
D iag o n al
p rin cip al
X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X X X
X X X
0 X X X X
0 0 X X X
X X X X
0
X X X 0 0
X X
O rg an izació n co n v en cio n al d e la
m atriz d e rig id ez
0
X 0 0
O rg an izació n en b an d a p ara O rg an izació n en b an d a p ara
m atrices n o sim étricas
m atrices sim étricas
Materiales y Métodos
Método Frontal
Este método trata de evitar el ensamblaje de la matriz
de rigidez global. Consiste básicamente en ensamblar la
menor cantidad de elementos y eliminar las variables al
mismo tiempo, formando matrices de menor tamaño.
Esto es, inmediatamente después de que los coeficientes
de una ecuación se ensamblan a partir de la contribución
de todos los elementos relevantes, o que inciden sobre este
coeficiente, la variable de este puede eliminarse. Como
resultado, la matriz completa de toda la estructura jamás
se forma.
Materiales y Métodos
Método Frontal
A diferencia del método banda, la adecuada
enumeración de los elemento es de gran importancia, no
así la enumeración global de los nodos ya que estos son
tratados por su numeración local. Como primer medida
esta característica del Front lo hace mas versátil ya que
como la cantidad de elementos siempre es menor que la
cantidad de nodos, independientemente de la forma de la
estructura, siempre se tendrá mayos facilidad para
enumerar los elementos que los nodos.
Discusión de Resultados
Para analizar las cualidades de los distintos métodos
propuestos se planteó la resolución mediante el empleo del
FECCUND V 2.1 (software desarrollado por el Dto. De
Mecánica Aplicada de la Facultad de Ingeniería de la
UNNE) de dos ejemplos prácticos, representando los
resultados mediante el pre-post procesador GID 7.2.
Para dicho análisis se tomaró como variables el tiempo
de ejecución de cada método y su exactitud, la cual será
determinada por comparación con los resultados
obtenidos por el software denominado CALSEF provisto
por el Centro Internacional de Métodos Numéricos en
Ingeniería (CIMNE).
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Mallado con elementos
hexaédricos
Desplazamiento nodal
según eje y
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Tabla Nº 1
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
1,758
0,501953E-04
0,47619E-04
5,13
Banda
0,552
0,498874E-04
4,55
Frontal
0,441
0,479135E-04
0,61
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Mallado con elementos hexaédricos
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Desplazamiento nodal según eje y
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Tabla Nº 2
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
31,547
-0,00615923
-0.00544118
11,66
Banda
2,025
-0,00560842
2,98
Frontal
1,752
-0,00552264
1,47
Conclusiones
El método Frontal es superior al Gauss y al Banda
en cuanto a su exactitud y a la velocidad de
ejecución.
El efecto de la propagación de errores por redondeo
incrementa notablemente la deficiencia del método
de Gauss frente a problemas de mayor
envergadura.
El método de Banda puede ser recomendable en el
caso de contar con un software de pre-post
procesamiento ya que los mismos generan las
mallas enumerando los nodos en la forma mas
adecuada disminuyendo de ancho de banda.
Caracterización cualitativa de
algunos métodos de solución de
sistemas de ecuaciones lineales.
MROGINSKI, Javier L. – BENEYTO, Pablo A. –
DI RADO, H. Ariel – MANZOLILLO, Juan E.
Departamento de Mecánica Aplicada.
Facultad de Ingeniería. U.N.N.E.
Antecedentes
El Método de los Elementos Finitos (M.E.F.) es una
técnica de análisis numérico que permite la obtención de
soluciones aproximadas de ecuaciones diferenciales capaces
de representar una variedad de problemas de la Ingeniería y
Física.
El resultado de la aplicación del Método de los Elementos
Finitos a un dominio cualquier es un sistema de ecuaciones
lineales simultáneas de gran tamaño, correspondiendo una
ecuación por cada incógnita que se evalúe en cada nodo,
imposible de ser resuelto sin el auxilio de un computador.
El presente trabajo desarrolla algunos aspectos
significativos sobre los tres diferentes métodos de resolución
de sistemas de ecuaciones lineales simultáneas mas usados en
el Método de los Elementos Finitos que son el: Eliminación
de Gauss, Sistema Banda y Frontal entre otros.
Materiales y Métodos
En las diapositivas siguientes
expondrá una breve descripción de
métodos de solución de sistemas
ecuaciones lineales simultáneas que
estudian en el presente trabajo.
se
los
de
se
Materiales y Métodos
Eliminación de Gauss
Un sistema de ecuaciones lineales simultáneas
puede ser expresado en forma genérica de la
siguiente manera:
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
a 21 x 1 a 22 x 2 a 23 x 3 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... .
a n1 x 1 a n 2 x 2 a n 3 x 3 a n n x n C n
Materiales y Métodos
Eliminación de Gauss
Este método busca obtener un sistema de ecuaciones
triangular equivalente, mediante el uso de transformaciones
lineales. El cual se resuelve con facilidad aplicando la
denominada sustitución inversa.
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
´
´
´
´
a 22 x 2 a 23 x 2 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... ........
(n2)
(n2)
(n2)
a n 1 , n 1 x n 1 a n 1 , n x n C n 1
( n 1 )
( n 1 )
a n,n x n C n
Materiales y Métodos
Sistema banda de solución
Se dice que una matriz es matriz banda cuando posee
todos sus elementos no nulos concentrados en torno de
la diagonal principal. Esta es una propiedad muy
interesante que se presenta en las matrices de rigidez
locales y globales.
El método banda de resolución de sistemas de
ecuaciones lineales simultaneas se basa en esta
propiedad, permitiendo el ahorro de asignaciones de
memoria, ya que solo guarda los elementos no nulos.
Materiales y Métodos
Sistema banda de solución
El presente método centra su eficiencia en la forma en que
se realiza la numeración global de los nodos.
S em ian ch o d e b an d a
N ro. total de incognitas
X X X
0 0 0 0
X X X X 0 0 0
X X X X X 0 0
0 X X X X X 0
0 0 X X X X X
0 0 0 X X X X
0 0 0 0 X X X
0 0 0 0 0 X X
0 0 0 0 0 0 X
0 0 0 0 0 0 0
0 0 0 0 0 0 0
0 0 0 0 0 0 0
A n ch o d e b an d a
0
0
0
0
0
0
0
0
0
0
X 0
0
0
0
0
0
0
X X 0
0
0
0
0
0
0
0
X X X 0
0
0
0
0
0
0
0
0
X X X X 0
0 0 X X X
0 X X X X
A n ch o d e b an d a
D iag o n al
p rin cip al
X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X X X
X X X
0 X X X X
0 0 X X X
X X X X
0
X X X 0 0
X X
O rg an izació n co n v en cio n al d e la
m atriz d e rig id ez
0
X 0 0
O rg an izació n en b an d a p ara O rg an izació n en b an d a p ara
m atrices n o sim étricas
m atrices sim étricas
Materiales y Métodos
Método Frontal
Este método trata de evitar el ensamblaje de la matriz
de rigidez global. Consiste básicamente en ensamblar la
menor cantidad de elementos y eliminar las variables al
mismo tiempo, formando matrices de menor tamaño.
Esto es, inmediatamente después de que los coeficientes
de una ecuación se ensamblan a partir de la contribución
de todos los elementos relevantes, o que inciden sobre este
coeficiente, la variable de este puede eliminarse. Como
resultado, la matriz completa de toda la estructura jamás
se forma.
Materiales y Métodos
Método Frontal
A diferencia del método banda, la adecuada
enumeración de los elemento es de gran importancia, no
así la enumeración global de los nodos ya que estos son
tratados por su numeración local. Como primer medida
esta característica del Front lo hace mas versátil ya que
como la cantidad de elementos siempre es menor que la
cantidad de nodos, independientemente de la forma de la
estructura, siempre se tendrá mayos facilidad para
enumerar los elementos que los nodos.
Discusión de Resultados
Para analizar las cualidades de los distintos métodos
propuestos se planteó la resolución mediante el empleo del
FECCUND V 2.1 (software desarrollado por el Dto. De
Mecánica Aplicada de la Facultad de Ingeniería de la
UNNE) de dos ejemplos prácticos, representando los
resultados mediante el pre-post procesador GID 7.2.
Para dicho análisis se tomaró como variables el tiempo
de ejecución de cada método y su exactitud, la cual será
determinada por comparación con los resultados
obtenidos por el software denominado CALSEF provisto
por el Centro Internacional de Métodos Numéricos en
Ingeniería (CIMNE).
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Mallado con elementos
hexaédricos
Desplazamiento nodal
según eje y
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Tabla Nº 1
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
1,758
0,501953E-04
0,47619E-04
5,13
Banda
0,552
0,498874E-04
4,55
Frontal
0,441
0,479135E-04
0,61
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Mallado con elementos hexaédricos
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Desplazamiento nodal según eje y
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Tabla Nº 2
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
31,547
-0,00615923
-0.00544118
11,66
Banda
2,025
-0,00560842
2,98
Frontal
1,752
-0,00552264
1,47
Conclusiones
El método Frontal es superior al Gauss y al Banda
en cuanto a su exactitud y a la velocidad de
ejecución.
El efecto de la propagación de errores por redondeo
incrementa notablemente la deficiencia del método
de Gauss frente a problemas de mayor
envergadura.
El método de Banda puede ser recomendable en el
caso de contar con un software de pre-post
procesamiento ya que los mismos generan las
mallas enumerando los nodos en la forma mas
adecuada disminuyendo de ancho de banda.
Slide 2
Caracterización cualitativa de
algunos métodos de solución de
sistemas de ecuaciones lineales.
MROGINSKI, Javier L. – BENEYTO, Pablo A. –
DI RADO, H. Ariel – MANZOLILLO, Juan E.
Departamento de Mecánica Aplicada.
Facultad de Ingeniería. U.N.N.E.
Antecedentes
El Método de los Elementos Finitos (M.E.F.) es una
técnica de análisis numérico que permite la obtención de
soluciones aproximadas de ecuaciones diferenciales capaces
de representar una variedad de problemas de la Ingeniería y
Física.
El resultado de la aplicación del Método de los Elementos
Finitos a un dominio cualquier es un sistema de ecuaciones
lineales simultáneas de gran tamaño, correspondiendo una
ecuación por cada incógnita que se evalúe en cada nodo,
imposible de ser resuelto sin el auxilio de un computador.
El presente trabajo desarrolla algunos aspectos
significativos sobre los tres diferentes métodos de resolución
de sistemas de ecuaciones lineales simultáneas mas usados en
el Método de los Elementos Finitos que son el: Eliminación
de Gauss, Sistema Banda y Frontal entre otros.
Materiales y Métodos
En las diapositivas siguientes
expondrá una breve descripción de
métodos de solución de sistemas
ecuaciones lineales simultáneas que
estudian en el presente trabajo.
se
los
de
se
Materiales y Métodos
Eliminación de Gauss
Un sistema de ecuaciones lineales simultáneas
puede ser expresado en forma genérica de la
siguiente manera:
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
a 21 x 1 a 22 x 2 a 23 x 3 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... .
a n1 x 1 a n 2 x 2 a n 3 x 3 a n n x n C n
Materiales y Métodos
Eliminación de Gauss
Este método busca obtener un sistema de ecuaciones
triangular equivalente, mediante el uso de transformaciones
lineales. El cual se resuelve con facilidad aplicando la
denominada sustitución inversa.
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
´
´
´
´
a 22 x 2 a 23 x 2 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... ........
(n2)
(n2)
(n2)
a n 1 , n 1 x n 1 a n 1 , n x n C n 1
( n 1 )
( n 1 )
a n,n x n C n
Materiales y Métodos
Sistema banda de solución
Se dice que una matriz es matriz banda cuando posee
todos sus elementos no nulos concentrados en torno de
la diagonal principal. Esta es una propiedad muy
interesante que se presenta en las matrices de rigidez
locales y globales.
El método banda de resolución de sistemas de
ecuaciones lineales simultaneas se basa en esta
propiedad, permitiendo el ahorro de asignaciones de
memoria, ya que solo guarda los elementos no nulos.
Materiales y Métodos
Sistema banda de solución
El presente método centra su eficiencia en la forma en que
se realiza la numeración global de los nodos.
S em ian ch o d e b an d a
N ro. total de incognitas
X X X
0 0 0 0
X X X X 0 0 0
X X X X X 0 0
0 X X X X X 0
0 0 X X X X X
0 0 0 X X X X
0 0 0 0 X X X
0 0 0 0 0 X X
0 0 0 0 0 0 X
0 0 0 0 0 0 0
0 0 0 0 0 0 0
0 0 0 0 0 0 0
A n ch o d e b an d a
0
0
0
0
0
0
0
0
0
0
X 0
0
0
0
0
0
0
X X 0
0
0
0
0
0
0
0
X X X 0
0
0
0
0
0
0
0
0
X X X X 0
0 0 X X X
0 X X X X
A n ch o d e b an d a
D iag o n al
p rin cip al
X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X X X
X X X
0 X X X X
0 0 X X X
X X X X
0
X X X 0 0
X X
O rg an izació n co n v en cio n al d e la
m atriz d e rig id ez
0
X 0 0
O rg an izació n en b an d a p ara O rg an izació n en b an d a p ara
m atrices n o sim étricas
m atrices sim étricas
Materiales y Métodos
Método Frontal
Este método trata de evitar el ensamblaje de la matriz
de rigidez global. Consiste básicamente en ensamblar la
menor cantidad de elementos y eliminar las variables al
mismo tiempo, formando matrices de menor tamaño.
Esto es, inmediatamente después de que los coeficientes
de una ecuación se ensamblan a partir de la contribución
de todos los elementos relevantes, o que inciden sobre este
coeficiente, la variable de este puede eliminarse. Como
resultado, la matriz completa de toda la estructura jamás
se forma.
Materiales y Métodos
Método Frontal
A diferencia del método banda, la adecuada
enumeración de los elemento es de gran importancia, no
así la enumeración global de los nodos ya que estos son
tratados por su numeración local. Como primer medida
esta característica del Front lo hace mas versátil ya que
como la cantidad de elementos siempre es menor que la
cantidad de nodos, independientemente de la forma de la
estructura, siempre se tendrá mayos facilidad para
enumerar los elementos que los nodos.
Discusión de Resultados
Para analizar las cualidades de los distintos métodos
propuestos se planteó la resolución mediante el empleo del
FECCUND V 2.1 (software desarrollado por el Dto. De
Mecánica Aplicada de la Facultad de Ingeniería de la
UNNE) de dos ejemplos prácticos, representando los
resultados mediante el pre-post procesador GID 7.2.
Para dicho análisis se tomaró como variables el tiempo
de ejecución de cada método y su exactitud, la cual será
determinada por comparación con los resultados
obtenidos por el software denominado CALSEF provisto
por el Centro Internacional de Métodos Numéricos en
Ingeniería (CIMNE).
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Mallado con elementos
hexaédricos
Desplazamiento nodal
según eje y
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Tabla Nº 1
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
1,758
0,501953E-04
0,47619E-04
5,13
Banda
0,552
0,498874E-04
4,55
Frontal
0,441
0,479135E-04
0,61
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Mallado con elementos hexaédricos
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Desplazamiento nodal según eje y
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Tabla Nº 2
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
31,547
-0,00615923
-0.00544118
11,66
Banda
2,025
-0,00560842
2,98
Frontal
1,752
-0,00552264
1,47
Conclusiones
El método Frontal es superior al Gauss y al Banda
en cuanto a su exactitud y a la velocidad de
ejecución.
El efecto de la propagación de errores por redondeo
incrementa notablemente la deficiencia del método
de Gauss frente a problemas de mayor
envergadura.
El método de Banda puede ser recomendable en el
caso de contar con un software de pre-post
procesamiento ya que los mismos generan las
mallas enumerando los nodos en la forma mas
adecuada disminuyendo de ancho de banda.
Slide 3
Caracterización cualitativa de
algunos métodos de solución de
sistemas de ecuaciones lineales.
MROGINSKI, Javier L. – BENEYTO, Pablo A. –
DI RADO, H. Ariel – MANZOLILLO, Juan E.
Departamento de Mecánica Aplicada.
Facultad de Ingeniería. U.N.N.E.
Antecedentes
El Método de los Elementos Finitos (M.E.F.) es una
técnica de análisis numérico que permite la obtención de
soluciones aproximadas de ecuaciones diferenciales capaces
de representar una variedad de problemas de la Ingeniería y
Física.
El resultado de la aplicación del Método de los Elementos
Finitos a un dominio cualquier es un sistema de ecuaciones
lineales simultáneas de gran tamaño, correspondiendo una
ecuación por cada incógnita que se evalúe en cada nodo,
imposible de ser resuelto sin el auxilio de un computador.
El presente trabajo desarrolla algunos aspectos
significativos sobre los tres diferentes métodos de resolución
de sistemas de ecuaciones lineales simultáneas mas usados en
el Método de los Elementos Finitos que son el: Eliminación
de Gauss, Sistema Banda y Frontal entre otros.
Materiales y Métodos
En las diapositivas siguientes
expondrá una breve descripción de
métodos de solución de sistemas
ecuaciones lineales simultáneas que
estudian en el presente trabajo.
se
los
de
se
Materiales y Métodos
Eliminación de Gauss
Un sistema de ecuaciones lineales simultáneas
puede ser expresado en forma genérica de la
siguiente manera:
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
a 21 x 1 a 22 x 2 a 23 x 3 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... .
a n1 x 1 a n 2 x 2 a n 3 x 3 a n n x n C n
Materiales y Métodos
Eliminación de Gauss
Este método busca obtener un sistema de ecuaciones
triangular equivalente, mediante el uso de transformaciones
lineales. El cual se resuelve con facilidad aplicando la
denominada sustitución inversa.
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
´
´
´
´
a 22 x 2 a 23 x 2 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... ........
(n2)
(n2)
(n2)
a n 1 , n 1 x n 1 a n 1 , n x n C n 1
( n 1 )
( n 1 )
a n,n x n C n
Materiales y Métodos
Sistema banda de solución
Se dice que una matriz es matriz banda cuando posee
todos sus elementos no nulos concentrados en torno de
la diagonal principal. Esta es una propiedad muy
interesante que se presenta en las matrices de rigidez
locales y globales.
El método banda de resolución de sistemas de
ecuaciones lineales simultaneas se basa en esta
propiedad, permitiendo el ahorro de asignaciones de
memoria, ya que solo guarda los elementos no nulos.
Materiales y Métodos
Sistema banda de solución
El presente método centra su eficiencia en la forma en que
se realiza la numeración global de los nodos.
S em ian ch o d e b an d a
N ro. total de incognitas
X X X
0 0 0 0
X X X X 0 0 0
X X X X X 0 0
0 X X X X X 0
0 0 X X X X X
0 0 0 X X X X
0 0 0 0 X X X
0 0 0 0 0 X X
0 0 0 0 0 0 X
0 0 0 0 0 0 0
0 0 0 0 0 0 0
0 0 0 0 0 0 0
A n ch o d e b an d a
0
0
0
0
0
0
0
0
0
0
X 0
0
0
0
0
0
0
X X 0
0
0
0
0
0
0
0
X X X 0
0
0
0
0
0
0
0
0
X X X X 0
0 0 X X X
0 X X X X
A n ch o d e b an d a
D iag o n al
p rin cip al
X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X X X
X X X
0 X X X X
0 0 X X X
X X X X
0
X X X 0 0
X X
O rg an izació n co n v en cio n al d e la
m atriz d e rig id ez
0
X 0 0
O rg an izació n en b an d a p ara O rg an izació n en b an d a p ara
m atrices n o sim étricas
m atrices sim étricas
Materiales y Métodos
Método Frontal
Este método trata de evitar el ensamblaje de la matriz
de rigidez global. Consiste básicamente en ensamblar la
menor cantidad de elementos y eliminar las variables al
mismo tiempo, formando matrices de menor tamaño.
Esto es, inmediatamente después de que los coeficientes
de una ecuación se ensamblan a partir de la contribución
de todos los elementos relevantes, o que inciden sobre este
coeficiente, la variable de este puede eliminarse. Como
resultado, la matriz completa de toda la estructura jamás
se forma.
Materiales y Métodos
Método Frontal
A diferencia del método banda, la adecuada
enumeración de los elemento es de gran importancia, no
así la enumeración global de los nodos ya que estos son
tratados por su numeración local. Como primer medida
esta característica del Front lo hace mas versátil ya que
como la cantidad de elementos siempre es menor que la
cantidad de nodos, independientemente de la forma de la
estructura, siempre se tendrá mayos facilidad para
enumerar los elementos que los nodos.
Discusión de Resultados
Para analizar las cualidades de los distintos métodos
propuestos se planteó la resolución mediante el empleo del
FECCUND V 2.1 (software desarrollado por el Dto. De
Mecánica Aplicada de la Facultad de Ingeniería de la
UNNE) de dos ejemplos prácticos, representando los
resultados mediante el pre-post procesador GID 7.2.
Para dicho análisis se tomaró como variables el tiempo
de ejecución de cada método y su exactitud, la cual será
determinada por comparación con los resultados
obtenidos por el software denominado CALSEF provisto
por el Centro Internacional de Métodos Numéricos en
Ingeniería (CIMNE).
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Mallado con elementos
hexaédricos
Desplazamiento nodal
según eje y
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Tabla Nº 1
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
1,758
0,501953E-04
0,47619E-04
5,13
Banda
0,552
0,498874E-04
4,55
Frontal
0,441
0,479135E-04
0,61
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Mallado con elementos hexaédricos
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Desplazamiento nodal según eje y
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Tabla Nº 2
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
31,547
-0,00615923
-0.00544118
11,66
Banda
2,025
-0,00560842
2,98
Frontal
1,752
-0,00552264
1,47
Conclusiones
El método Frontal es superior al Gauss y al Banda
en cuanto a su exactitud y a la velocidad de
ejecución.
El efecto de la propagación de errores por redondeo
incrementa notablemente la deficiencia del método
de Gauss frente a problemas de mayor
envergadura.
El método de Banda puede ser recomendable en el
caso de contar con un software de pre-post
procesamiento ya que los mismos generan las
mallas enumerando los nodos en la forma mas
adecuada disminuyendo de ancho de banda.
Slide 4
Caracterización cualitativa de
algunos métodos de solución de
sistemas de ecuaciones lineales.
MROGINSKI, Javier L. – BENEYTO, Pablo A. –
DI RADO, H. Ariel – MANZOLILLO, Juan E.
Departamento de Mecánica Aplicada.
Facultad de Ingeniería. U.N.N.E.
Antecedentes
El Método de los Elementos Finitos (M.E.F.) es una
técnica de análisis numérico que permite la obtención de
soluciones aproximadas de ecuaciones diferenciales capaces
de representar una variedad de problemas de la Ingeniería y
Física.
El resultado de la aplicación del Método de los Elementos
Finitos a un dominio cualquier es un sistema de ecuaciones
lineales simultáneas de gran tamaño, correspondiendo una
ecuación por cada incógnita que se evalúe en cada nodo,
imposible de ser resuelto sin el auxilio de un computador.
El presente trabajo desarrolla algunos aspectos
significativos sobre los tres diferentes métodos de resolución
de sistemas de ecuaciones lineales simultáneas mas usados en
el Método de los Elementos Finitos que son el: Eliminación
de Gauss, Sistema Banda y Frontal entre otros.
Materiales y Métodos
En las diapositivas siguientes
expondrá una breve descripción de
métodos de solución de sistemas
ecuaciones lineales simultáneas que
estudian en el presente trabajo.
se
los
de
se
Materiales y Métodos
Eliminación de Gauss
Un sistema de ecuaciones lineales simultáneas
puede ser expresado en forma genérica de la
siguiente manera:
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
a 21 x 1 a 22 x 2 a 23 x 3 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... .
a n1 x 1 a n 2 x 2 a n 3 x 3 a n n x n C n
Materiales y Métodos
Eliminación de Gauss
Este método busca obtener un sistema de ecuaciones
triangular equivalente, mediante el uso de transformaciones
lineales. El cual se resuelve con facilidad aplicando la
denominada sustitución inversa.
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
´
´
´
´
a 22 x 2 a 23 x 2 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... ........
(n2)
(n2)
(n2)
a n 1 , n 1 x n 1 a n 1 , n x n C n 1
( n 1 )
( n 1 )
a n,n x n C n
Materiales y Métodos
Sistema banda de solución
Se dice que una matriz es matriz banda cuando posee
todos sus elementos no nulos concentrados en torno de
la diagonal principal. Esta es una propiedad muy
interesante que se presenta en las matrices de rigidez
locales y globales.
El método banda de resolución de sistemas de
ecuaciones lineales simultaneas se basa en esta
propiedad, permitiendo el ahorro de asignaciones de
memoria, ya que solo guarda los elementos no nulos.
Materiales y Métodos
Sistema banda de solución
El presente método centra su eficiencia en la forma en que
se realiza la numeración global de los nodos.
S em ian ch o d e b an d a
N ro. total de incognitas
X X X
0 0 0 0
X X X X 0 0 0
X X X X X 0 0
0 X X X X X 0
0 0 X X X X X
0 0 0 X X X X
0 0 0 0 X X X
0 0 0 0 0 X X
0 0 0 0 0 0 X
0 0 0 0 0 0 0
0 0 0 0 0 0 0
0 0 0 0 0 0 0
A n ch o d e b an d a
0
0
0
0
0
0
0
0
0
0
X 0
0
0
0
0
0
0
X X 0
0
0
0
0
0
0
0
X X X 0
0
0
0
0
0
0
0
0
X X X X 0
0 0 X X X
0 X X X X
A n ch o d e b an d a
D iag o n al
p rin cip al
X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X X X
X X X
0 X X X X
0 0 X X X
X X X X
0
X X X 0 0
X X
O rg an izació n co n v en cio n al d e la
m atriz d e rig id ez
0
X 0 0
O rg an izació n en b an d a p ara O rg an izació n en b an d a p ara
m atrices n o sim étricas
m atrices sim étricas
Materiales y Métodos
Método Frontal
Este método trata de evitar el ensamblaje de la matriz
de rigidez global. Consiste básicamente en ensamblar la
menor cantidad de elementos y eliminar las variables al
mismo tiempo, formando matrices de menor tamaño.
Esto es, inmediatamente después de que los coeficientes
de una ecuación se ensamblan a partir de la contribución
de todos los elementos relevantes, o que inciden sobre este
coeficiente, la variable de este puede eliminarse. Como
resultado, la matriz completa de toda la estructura jamás
se forma.
Materiales y Métodos
Método Frontal
A diferencia del método banda, la adecuada
enumeración de los elemento es de gran importancia, no
así la enumeración global de los nodos ya que estos son
tratados por su numeración local. Como primer medida
esta característica del Front lo hace mas versátil ya que
como la cantidad de elementos siempre es menor que la
cantidad de nodos, independientemente de la forma de la
estructura, siempre se tendrá mayos facilidad para
enumerar los elementos que los nodos.
Discusión de Resultados
Para analizar las cualidades de los distintos métodos
propuestos se planteó la resolución mediante el empleo del
FECCUND V 2.1 (software desarrollado por el Dto. De
Mecánica Aplicada de la Facultad de Ingeniería de la
UNNE) de dos ejemplos prácticos, representando los
resultados mediante el pre-post procesador GID 7.2.
Para dicho análisis se tomaró como variables el tiempo
de ejecución de cada método y su exactitud, la cual será
determinada por comparación con los resultados
obtenidos por el software denominado CALSEF provisto
por el Centro Internacional de Métodos Numéricos en
Ingeniería (CIMNE).
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Mallado con elementos
hexaédricos
Desplazamiento nodal
según eje y
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Tabla Nº 1
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
1,758
0,501953E-04
0,47619E-04
5,13
Banda
0,552
0,498874E-04
4,55
Frontal
0,441
0,479135E-04
0,61
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Mallado con elementos hexaédricos
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Desplazamiento nodal según eje y
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Tabla Nº 2
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
31,547
-0,00615923
-0.00544118
11,66
Banda
2,025
-0,00560842
2,98
Frontal
1,752
-0,00552264
1,47
Conclusiones
El método Frontal es superior al Gauss y al Banda
en cuanto a su exactitud y a la velocidad de
ejecución.
El efecto de la propagación de errores por redondeo
incrementa notablemente la deficiencia del método
de Gauss frente a problemas de mayor
envergadura.
El método de Banda puede ser recomendable en el
caso de contar con un software de pre-post
procesamiento ya que los mismos generan las
mallas enumerando los nodos en la forma mas
adecuada disminuyendo de ancho de banda.
Slide 5
Caracterización cualitativa de
algunos métodos de solución de
sistemas de ecuaciones lineales.
MROGINSKI, Javier L. – BENEYTO, Pablo A. –
DI RADO, H. Ariel – MANZOLILLO, Juan E.
Departamento de Mecánica Aplicada.
Facultad de Ingeniería. U.N.N.E.
Antecedentes
El Método de los Elementos Finitos (M.E.F.) es una
técnica de análisis numérico que permite la obtención de
soluciones aproximadas de ecuaciones diferenciales capaces
de representar una variedad de problemas de la Ingeniería y
Física.
El resultado de la aplicación del Método de los Elementos
Finitos a un dominio cualquier es un sistema de ecuaciones
lineales simultáneas de gran tamaño, correspondiendo una
ecuación por cada incógnita que se evalúe en cada nodo,
imposible de ser resuelto sin el auxilio de un computador.
El presente trabajo desarrolla algunos aspectos
significativos sobre los tres diferentes métodos de resolución
de sistemas de ecuaciones lineales simultáneas mas usados en
el Método de los Elementos Finitos que son el: Eliminación
de Gauss, Sistema Banda y Frontal entre otros.
Materiales y Métodos
En las diapositivas siguientes
expondrá una breve descripción de
métodos de solución de sistemas
ecuaciones lineales simultáneas que
estudian en el presente trabajo.
se
los
de
se
Materiales y Métodos
Eliminación de Gauss
Un sistema de ecuaciones lineales simultáneas
puede ser expresado en forma genérica de la
siguiente manera:
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
a 21 x 1 a 22 x 2 a 23 x 3 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... .
a n1 x 1 a n 2 x 2 a n 3 x 3 a n n x n C n
Materiales y Métodos
Eliminación de Gauss
Este método busca obtener un sistema de ecuaciones
triangular equivalente, mediante el uso de transformaciones
lineales. El cual se resuelve con facilidad aplicando la
denominada sustitución inversa.
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
´
´
´
´
a 22 x 2 a 23 x 2 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... ........
(n2)
(n2)
(n2)
a n 1 , n 1 x n 1 a n 1 , n x n C n 1
( n 1 )
( n 1 )
a n,n x n C n
Materiales y Métodos
Sistema banda de solución
Se dice que una matriz es matriz banda cuando posee
todos sus elementos no nulos concentrados en torno de
la diagonal principal. Esta es una propiedad muy
interesante que se presenta en las matrices de rigidez
locales y globales.
El método banda de resolución de sistemas de
ecuaciones lineales simultaneas se basa en esta
propiedad, permitiendo el ahorro de asignaciones de
memoria, ya que solo guarda los elementos no nulos.
Materiales y Métodos
Sistema banda de solución
El presente método centra su eficiencia en la forma en que
se realiza la numeración global de los nodos.
S em ian ch o d e b an d a
N ro. total de incognitas
X X X
0 0 0 0
X X X X 0 0 0
X X X X X 0 0
0 X X X X X 0
0 0 X X X X X
0 0 0 X X X X
0 0 0 0 X X X
0 0 0 0 0 X X
0 0 0 0 0 0 X
0 0 0 0 0 0 0
0 0 0 0 0 0 0
0 0 0 0 0 0 0
A n ch o d e b an d a
0
0
0
0
0
0
0
0
0
0
X 0
0
0
0
0
0
0
X X 0
0
0
0
0
0
0
0
X X X 0
0
0
0
0
0
0
0
0
X X X X 0
0 0 X X X
0 X X X X
A n ch o d e b an d a
D iag o n al
p rin cip al
X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X X X
X X X
0 X X X X
0 0 X X X
X X X X
0
X X X 0 0
X X
O rg an izació n co n v en cio n al d e la
m atriz d e rig id ez
0
X 0 0
O rg an izació n en b an d a p ara O rg an izació n en b an d a p ara
m atrices n o sim étricas
m atrices sim étricas
Materiales y Métodos
Método Frontal
Este método trata de evitar el ensamblaje de la matriz
de rigidez global. Consiste básicamente en ensamblar la
menor cantidad de elementos y eliminar las variables al
mismo tiempo, formando matrices de menor tamaño.
Esto es, inmediatamente después de que los coeficientes
de una ecuación se ensamblan a partir de la contribución
de todos los elementos relevantes, o que inciden sobre este
coeficiente, la variable de este puede eliminarse. Como
resultado, la matriz completa de toda la estructura jamás
se forma.
Materiales y Métodos
Método Frontal
A diferencia del método banda, la adecuada
enumeración de los elemento es de gran importancia, no
así la enumeración global de los nodos ya que estos son
tratados por su numeración local. Como primer medida
esta característica del Front lo hace mas versátil ya que
como la cantidad de elementos siempre es menor que la
cantidad de nodos, independientemente de la forma de la
estructura, siempre se tendrá mayos facilidad para
enumerar los elementos que los nodos.
Discusión de Resultados
Para analizar las cualidades de los distintos métodos
propuestos se planteó la resolución mediante el empleo del
FECCUND V 2.1 (software desarrollado por el Dto. De
Mecánica Aplicada de la Facultad de Ingeniería de la
UNNE) de dos ejemplos prácticos, representando los
resultados mediante el pre-post procesador GID 7.2.
Para dicho análisis se tomaró como variables el tiempo
de ejecución de cada método y su exactitud, la cual será
determinada por comparación con los resultados
obtenidos por el software denominado CALSEF provisto
por el Centro Internacional de Métodos Numéricos en
Ingeniería (CIMNE).
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Mallado con elementos
hexaédricos
Desplazamiento nodal
según eje y
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Tabla Nº 1
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
1,758
0,501953E-04
0,47619E-04
5,13
Banda
0,552
0,498874E-04
4,55
Frontal
0,441
0,479135E-04
0,61
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Mallado con elementos hexaédricos
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Desplazamiento nodal según eje y
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Tabla Nº 2
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
31,547
-0,00615923
-0.00544118
11,66
Banda
2,025
-0,00560842
2,98
Frontal
1,752
-0,00552264
1,47
Conclusiones
El método Frontal es superior al Gauss y al Banda
en cuanto a su exactitud y a la velocidad de
ejecución.
El efecto de la propagación de errores por redondeo
incrementa notablemente la deficiencia del método
de Gauss frente a problemas de mayor
envergadura.
El método de Banda puede ser recomendable en el
caso de contar con un software de pre-post
procesamiento ya que los mismos generan las
mallas enumerando los nodos en la forma mas
adecuada disminuyendo de ancho de banda.
Slide 6
Caracterización cualitativa de
algunos métodos de solución de
sistemas de ecuaciones lineales.
MROGINSKI, Javier L. – BENEYTO, Pablo A. –
DI RADO, H. Ariel – MANZOLILLO, Juan E.
Departamento de Mecánica Aplicada.
Facultad de Ingeniería. U.N.N.E.
Antecedentes
El Método de los Elementos Finitos (M.E.F.) es una
técnica de análisis numérico que permite la obtención de
soluciones aproximadas de ecuaciones diferenciales capaces
de representar una variedad de problemas de la Ingeniería y
Física.
El resultado de la aplicación del Método de los Elementos
Finitos a un dominio cualquier es un sistema de ecuaciones
lineales simultáneas de gran tamaño, correspondiendo una
ecuación por cada incógnita que se evalúe en cada nodo,
imposible de ser resuelto sin el auxilio de un computador.
El presente trabajo desarrolla algunos aspectos
significativos sobre los tres diferentes métodos de resolución
de sistemas de ecuaciones lineales simultáneas mas usados en
el Método de los Elementos Finitos que son el: Eliminación
de Gauss, Sistema Banda y Frontal entre otros.
Materiales y Métodos
En las diapositivas siguientes
expondrá una breve descripción de
métodos de solución de sistemas
ecuaciones lineales simultáneas que
estudian en el presente trabajo.
se
los
de
se
Materiales y Métodos
Eliminación de Gauss
Un sistema de ecuaciones lineales simultáneas
puede ser expresado en forma genérica de la
siguiente manera:
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
a 21 x 1 a 22 x 2 a 23 x 3 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... .
a n1 x 1 a n 2 x 2 a n 3 x 3 a n n x n C n
Materiales y Métodos
Eliminación de Gauss
Este método busca obtener un sistema de ecuaciones
triangular equivalente, mediante el uso de transformaciones
lineales. El cual se resuelve con facilidad aplicando la
denominada sustitución inversa.
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
´
´
´
´
a 22 x 2 a 23 x 2 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... ........
(n2)
(n2)
(n2)
a n 1 , n 1 x n 1 a n 1 , n x n C n 1
( n 1 )
( n 1 )
a n,n x n C n
Materiales y Métodos
Sistema banda de solución
Se dice que una matriz es matriz banda cuando posee
todos sus elementos no nulos concentrados en torno de
la diagonal principal. Esta es una propiedad muy
interesante que se presenta en las matrices de rigidez
locales y globales.
El método banda de resolución de sistemas de
ecuaciones lineales simultaneas se basa en esta
propiedad, permitiendo el ahorro de asignaciones de
memoria, ya que solo guarda los elementos no nulos.
Materiales y Métodos
Sistema banda de solución
El presente método centra su eficiencia en la forma en que
se realiza la numeración global de los nodos.
S em ian ch o d e b an d a
N ro. total de incognitas
X X X
0 0 0 0
X X X X 0 0 0
X X X X X 0 0
0 X X X X X 0
0 0 X X X X X
0 0 0 X X X X
0 0 0 0 X X X
0 0 0 0 0 X X
0 0 0 0 0 0 X
0 0 0 0 0 0 0
0 0 0 0 0 0 0
0 0 0 0 0 0 0
A n ch o d e b an d a
0
0
0
0
0
0
0
0
0
0
X 0
0
0
0
0
0
0
X X 0
0
0
0
0
0
0
0
X X X 0
0
0
0
0
0
0
0
0
X X X X 0
0 0 X X X
0 X X X X
A n ch o d e b an d a
D iag o n al
p rin cip al
X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X X X
X X X
0 X X X X
0 0 X X X
X X X X
0
X X X 0 0
X X
O rg an izació n co n v en cio n al d e la
m atriz d e rig id ez
0
X 0 0
O rg an izació n en b an d a p ara O rg an izació n en b an d a p ara
m atrices n o sim étricas
m atrices sim étricas
Materiales y Métodos
Método Frontal
Este método trata de evitar el ensamblaje de la matriz
de rigidez global. Consiste básicamente en ensamblar la
menor cantidad de elementos y eliminar las variables al
mismo tiempo, formando matrices de menor tamaño.
Esto es, inmediatamente después de que los coeficientes
de una ecuación se ensamblan a partir de la contribución
de todos los elementos relevantes, o que inciden sobre este
coeficiente, la variable de este puede eliminarse. Como
resultado, la matriz completa de toda la estructura jamás
se forma.
Materiales y Métodos
Método Frontal
A diferencia del método banda, la adecuada
enumeración de los elemento es de gran importancia, no
así la enumeración global de los nodos ya que estos son
tratados por su numeración local. Como primer medida
esta característica del Front lo hace mas versátil ya que
como la cantidad de elementos siempre es menor que la
cantidad de nodos, independientemente de la forma de la
estructura, siempre se tendrá mayos facilidad para
enumerar los elementos que los nodos.
Discusión de Resultados
Para analizar las cualidades de los distintos métodos
propuestos se planteó la resolución mediante el empleo del
FECCUND V 2.1 (software desarrollado por el Dto. De
Mecánica Aplicada de la Facultad de Ingeniería de la
UNNE) de dos ejemplos prácticos, representando los
resultados mediante el pre-post procesador GID 7.2.
Para dicho análisis se tomaró como variables el tiempo
de ejecución de cada método y su exactitud, la cual será
determinada por comparación con los resultados
obtenidos por el software denominado CALSEF provisto
por el Centro Internacional de Métodos Numéricos en
Ingeniería (CIMNE).
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Mallado con elementos
hexaédricos
Desplazamiento nodal
según eje y
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Tabla Nº 1
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
1,758
0,501953E-04
0,47619E-04
5,13
Banda
0,552
0,498874E-04
4,55
Frontal
0,441
0,479135E-04
0,61
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Mallado con elementos hexaédricos
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Desplazamiento nodal según eje y
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Tabla Nº 2
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
31,547
-0,00615923
-0.00544118
11,66
Banda
2,025
-0,00560842
2,98
Frontal
1,752
-0,00552264
1,47
Conclusiones
El método Frontal es superior al Gauss y al Banda
en cuanto a su exactitud y a la velocidad de
ejecución.
El efecto de la propagación de errores por redondeo
incrementa notablemente la deficiencia del método
de Gauss frente a problemas de mayor
envergadura.
El método de Banda puede ser recomendable en el
caso de contar con un software de pre-post
procesamiento ya que los mismos generan las
mallas enumerando los nodos en la forma mas
adecuada disminuyendo de ancho de banda.
Slide 7
Caracterización cualitativa de
algunos métodos de solución de
sistemas de ecuaciones lineales.
MROGINSKI, Javier L. – BENEYTO, Pablo A. –
DI RADO, H. Ariel – MANZOLILLO, Juan E.
Departamento de Mecánica Aplicada.
Facultad de Ingeniería. U.N.N.E.
Antecedentes
El Método de los Elementos Finitos (M.E.F.) es una
técnica de análisis numérico que permite la obtención de
soluciones aproximadas de ecuaciones diferenciales capaces
de representar una variedad de problemas de la Ingeniería y
Física.
El resultado de la aplicación del Método de los Elementos
Finitos a un dominio cualquier es un sistema de ecuaciones
lineales simultáneas de gran tamaño, correspondiendo una
ecuación por cada incógnita que se evalúe en cada nodo,
imposible de ser resuelto sin el auxilio de un computador.
El presente trabajo desarrolla algunos aspectos
significativos sobre los tres diferentes métodos de resolución
de sistemas de ecuaciones lineales simultáneas mas usados en
el Método de los Elementos Finitos que son el: Eliminación
de Gauss, Sistema Banda y Frontal entre otros.
Materiales y Métodos
En las diapositivas siguientes
expondrá una breve descripción de
métodos de solución de sistemas
ecuaciones lineales simultáneas que
estudian en el presente trabajo.
se
los
de
se
Materiales y Métodos
Eliminación de Gauss
Un sistema de ecuaciones lineales simultáneas
puede ser expresado en forma genérica de la
siguiente manera:
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
a 21 x 1 a 22 x 2 a 23 x 3 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... .
a n1 x 1 a n 2 x 2 a n 3 x 3 a n n x n C n
Materiales y Métodos
Eliminación de Gauss
Este método busca obtener un sistema de ecuaciones
triangular equivalente, mediante el uso de transformaciones
lineales. El cual se resuelve con facilidad aplicando la
denominada sustitución inversa.
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
´
´
´
´
a 22 x 2 a 23 x 2 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... ........
(n2)
(n2)
(n2)
a n 1 , n 1 x n 1 a n 1 , n x n C n 1
( n 1 )
( n 1 )
a n,n x n C n
Materiales y Métodos
Sistema banda de solución
Se dice que una matriz es matriz banda cuando posee
todos sus elementos no nulos concentrados en torno de
la diagonal principal. Esta es una propiedad muy
interesante que se presenta en las matrices de rigidez
locales y globales.
El método banda de resolución de sistemas de
ecuaciones lineales simultaneas se basa en esta
propiedad, permitiendo el ahorro de asignaciones de
memoria, ya que solo guarda los elementos no nulos.
Materiales y Métodos
Sistema banda de solución
El presente método centra su eficiencia en la forma en que
se realiza la numeración global de los nodos.
S em ian ch o d e b an d a
N ro. total de incognitas
X X X
0 0 0 0
X X X X 0 0 0
X X X X X 0 0
0 X X X X X 0
0 0 X X X X X
0 0 0 X X X X
0 0 0 0 X X X
0 0 0 0 0 X X
0 0 0 0 0 0 X
0 0 0 0 0 0 0
0 0 0 0 0 0 0
0 0 0 0 0 0 0
A n ch o d e b an d a
0
0
0
0
0
0
0
0
0
0
X 0
0
0
0
0
0
0
X X 0
0
0
0
0
0
0
0
X X X 0
0
0
0
0
0
0
0
0
X X X X 0
0 0 X X X
0 X X X X
A n ch o d e b an d a
D iag o n al
p rin cip al
X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X X X
X X X
0 X X X X
0 0 X X X
X X X X
0
X X X 0 0
X X
O rg an izació n co n v en cio n al d e la
m atriz d e rig id ez
0
X 0 0
O rg an izació n en b an d a p ara O rg an izació n en b an d a p ara
m atrices n o sim étricas
m atrices sim étricas
Materiales y Métodos
Método Frontal
Este método trata de evitar el ensamblaje de la matriz
de rigidez global. Consiste básicamente en ensamblar la
menor cantidad de elementos y eliminar las variables al
mismo tiempo, formando matrices de menor tamaño.
Esto es, inmediatamente después de que los coeficientes
de una ecuación se ensamblan a partir de la contribución
de todos los elementos relevantes, o que inciden sobre este
coeficiente, la variable de este puede eliminarse. Como
resultado, la matriz completa de toda la estructura jamás
se forma.
Materiales y Métodos
Método Frontal
A diferencia del método banda, la adecuada
enumeración de los elemento es de gran importancia, no
así la enumeración global de los nodos ya que estos son
tratados por su numeración local. Como primer medida
esta característica del Front lo hace mas versátil ya que
como la cantidad de elementos siempre es menor que la
cantidad de nodos, independientemente de la forma de la
estructura, siempre se tendrá mayos facilidad para
enumerar los elementos que los nodos.
Discusión de Resultados
Para analizar las cualidades de los distintos métodos
propuestos se planteó la resolución mediante el empleo del
FECCUND V 2.1 (software desarrollado por el Dto. De
Mecánica Aplicada de la Facultad de Ingeniería de la
UNNE) de dos ejemplos prácticos, representando los
resultados mediante el pre-post procesador GID 7.2.
Para dicho análisis se tomaró como variables el tiempo
de ejecución de cada método y su exactitud, la cual será
determinada por comparación con los resultados
obtenidos por el software denominado CALSEF provisto
por el Centro Internacional de Métodos Numéricos en
Ingeniería (CIMNE).
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Mallado con elementos
hexaédricos
Desplazamiento nodal
según eje y
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Tabla Nº 1
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
1,758
0,501953E-04
0,47619E-04
5,13
Banda
0,552
0,498874E-04
4,55
Frontal
0,441
0,479135E-04
0,61
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Mallado con elementos hexaédricos
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Desplazamiento nodal según eje y
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Tabla Nº 2
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
31,547
-0,00615923
-0.00544118
11,66
Banda
2,025
-0,00560842
2,98
Frontal
1,752
-0,00552264
1,47
Conclusiones
El método Frontal es superior al Gauss y al Banda
en cuanto a su exactitud y a la velocidad de
ejecución.
El efecto de la propagación de errores por redondeo
incrementa notablemente la deficiencia del método
de Gauss frente a problemas de mayor
envergadura.
El método de Banda puede ser recomendable en el
caso de contar con un software de pre-post
procesamiento ya que los mismos generan las
mallas enumerando los nodos en la forma mas
adecuada disminuyendo de ancho de banda.
Slide 8
Caracterización cualitativa de
algunos métodos de solución de
sistemas de ecuaciones lineales.
MROGINSKI, Javier L. – BENEYTO, Pablo A. –
DI RADO, H. Ariel – MANZOLILLO, Juan E.
Departamento de Mecánica Aplicada.
Facultad de Ingeniería. U.N.N.E.
Antecedentes
El Método de los Elementos Finitos (M.E.F.) es una
técnica de análisis numérico que permite la obtención de
soluciones aproximadas de ecuaciones diferenciales capaces
de representar una variedad de problemas de la Ingeniería y
Física.
El resultado de la aplicación del Método de los Elementos
Finitos a un dominio cualquier es un sistema de ecuaciones
lineales simultáneas de gran tamaño, correspondiendo una
ecuación por cada incógnita que se evalúe en cada nodo,
imposible de ser resuelto sin el auxilio de un computador.
El presente trabajo desarrolla algunos aspectos
significativos sobre los tres diferentes métodos de resolución
de sistemas de ecuaciones lineales simultáneas mas usados en
el Método de los Elementos Finitos que son el: Eliminación
de Gauss, Sistema Banda y Frontal entre otros.
Materiales y Métodos
En las diapositivas siguientes
expondrá una breve descripción de
métodos de solución de sistemas
ecuaciones lineales simultáneas que
estudian en el presente trabajo.
se
los
de
se
Materiales y Métodos
Eliminación de Gauss
Un sistema de ecuaciones lineales simultáneas
puede ser expresado en forma genérica de la
siguiente manera:
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
a 21 x 1 a 22 x 2 a 23 x 3 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... .
a n1 x 1 a n 2 x 2 a n 3 x 3 a n n x n C n
Materiales y Métodos
Eliminación de Gauss
Este método busca obtener un sistema de ecuaciones
triangular equivalente, mediante el uso de transformaciones
lineales. El cual se resuelve con facilidad aplicando la
denominada sustitución inversa.
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
´
´
´
´
a 22 x 2 a 23 x 2 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... ........
(n2)
(n2)
(n2)
a n 1 , n 1 x n 1 a n 1 , n x n C n 1
( n 1 )
( n 1 )
a n,n x n C n
Materiales y Métodos
Sistema banda de solución
Se dice que una matriz es matriz banda cuando posee
todos sus elementos no nulos concentrados en torno de
la diagonal principal. Esta es una propiedad muy
interesante que se presenta en las matrices de rigidez
locales y globales.
El método banda de resolución de sistemas de
ecuaciones lineales simultaneas se basa en esta
propiedad, permitiendo el ahorro de asignaciones de
memoria, ya que solo guarda los elementos no nulos.
Materiales y Métodos
Sistema banda de solución
El presente método centra su eficiencia en la forma en que
se realiza la numeración global de los nodos.
S em ian ch o d e b an d a
N ro. total de incognitas
X X X
0 0 0 0
X X X X 0 0 0
X X X X X 0 0
0 X X X X X 0
0 0 X X X X X
0 0 0 X X X X
0 0 0 0 X X X
0 0 0 0 0 X X
0 0 0 0 0 0 X
0 0 0 0 0 0 0
0 0 0 0 0 0 0
0 0 0 0 0 0 0
A n ch o d e b an d a
0
0
0
0
0
0
0
0
0
0
X 0
0
0
0
0
0
0
X X 0
0
0
0
0
0
0
0
X X X 0
0
0
0
0
0
0
0
0
X X X X 0
0 0 X X X
0 X X X X
A n ch o d e b an d a
D iag o n al
p rin cip al
X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X X X
X X X
0 X X X X
0 0 X X X
X X X X
0
X X X 0 0
X X
O rg an izació n co n v en cio n al d e la
m atriz d e rig id ez
0
X 0 0
O rg an izació n en b an d a p ara O rg an izació n en b an d a p ara
m atrices n o sim étricas
m atrices sim étricas
Materiales y Métodos
Método Frontal
Este método trata de evitar el ensamblaje de la matriz
de rigidez global. Consiste básicamente en ensamblar la
menor cantidad de elementos y eliminar las variables al
mismo tiempo, formando matrices de menor tamaño.
Esto es, inmediatamente después de que los coeficientes
de una ecuación se ensamblan a partir de la contribución
de todos los elementos relevantes, o que inciden sobre este
coeficiente, la variable de este puede eliminarse. Como
resultado, la matriz completa de toda la estructura jamás
se forma.
Materiales y Métodos
Método Frontal
A diferencia del método banda, la adecuada
enumeración de los elemento es de gran importancia, no
así la enumeración global de los nodos ya que estos son
tratados por su numeración local. Como primer medida
esta característica del Front lo hace mas versátil ya que
como la cantidad de elementos siempre es menor que la
cantidad de nodos, independientemente de la forma de la
estructura, siempre se tendrá mayos facilidad para
enumerar los elementos que los nodos.
Discusión de Resultados
Para analizar las cualidades de los distintos métodos
propuestos se planteó la resolución mediante el empleo del
FECCUND V 2.1 (software desarrollado por el Dto. De
Mecánica Aplicada de la Facultad de Ingeniería de la
UNNE) de dos ejemplos prácticos, representando los
resultados mediante el pre-post procesador GID 7.2.
Para dicho análisis se tomaró como variables el tiempo
de ejecución de cada método y su exactitud, la cual será
determinada por comparación con los resultados
obtenidos por el software denominado CALSEF provisto
por el Centro Internacional de Métodos Numéricos en
Ingeniería (CIMNE).
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Mallado con elementos
hexaédricos
Desplazamiento nodal
según eje y
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Tabla Nº 1
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
1,758
0,501953E-04
0,47619E-04
5,13
Banda
0,552
0,498874E-04
4,55
Frontal
0,441
0,479135E-04
0,61
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Mallado con elementos hexaédricos
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Desplazamiento nodal según eje y
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Tabla Nº 2
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
31,547
-0,00615923
-0.00544118
11,66
Banda
2,025
-0,00560842
2,98
Frontal
1,752
-0,00552264
1,47
Conclusiones
El método Frontal es superior al Gauss y al Banda
en cuanto a su exactitud y a la velocidad de
ejecución.
El efecto de la propagación de errores por redondeo
incrementa notablemente la deficiencia del método
de Gauss frente a problemas de mayor
envergadura.
El método de Banda puede ser recomendable en el
caso de contar con un software de pre-post
procesamiento ya que los mismos generan las
mallas enumerando los nodos en la forma mas
adecuada disminuyendo de ancho de banda.
Slide 9
Caracterización cualitativa de
algunos métodos de solución de
sistemas de ecuaciones lineales.
MROGINSKI, Javier L. – BENEYTO, Pablo A. –
DI RADO, H. Ariel – MANZOLILLO, Juan E.
Departamento de Mecánica Aplicada.
Facultad de Ingeniería. U.N.N.E.
Antecedentes
El Método de los Elementos Finitos (M.E.F.) es una
técnica de análisis numérico que permite la obtención de
soluciones aproximadas de ecuaciones diferenciales capaces
de representar una variedad de problemas de la Ingeniería y
Física.
El resultado de la aplicación del Método de los Elementos
Finitos a un dominio cualquier es un sistema de ecuaciones
lineales simultáneas de gran tamaño, correspondiendo una
ecuación por cada incógnita que se evalúe en cada nodo,
imposible de ser resuelto sin el auxilio de un computador.
El presente trabajo desarrolla algunos aspectos
significativos sobre los tres diferentes métodos de resolución
de sistemas de ecuaciones lineales simultáneas mas usados en
el Método de los Elementos Finitos que son el: Eliminación
de Gauss, Sistema Banda y Frontal entre otros.
Materiales y Métodos
En las diapositivas siguientes
expondrá una breve descripción de
métodos de solución de sistemas
ecuaciones lineales simultáneas que
estudian en el presente trabajo.
se
los
de
se
Materiales y Métodos
Eliminación de Gauss
Un sistema de ecuaciones lineales simultáneas
puede ser expresado en forma genérica de la
siguiente manera:
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
a 21 x 1 a 22 x 2 a 23 x 3 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... .
a n1 x 1 a n 2 x 2 a n 3 x 3 a n n x n C n
Materiales y Métodos
Eliminación de Gauss
Este método busca obtener un sistema de ecuaciones
triangular equivalente, mediante el uso de transformaciones
lineales. El cual se resuelve con facilidad aplicando la
denominada sustitución inversa.
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
´
´
´
´
a 22 x 2 a 23 x 2 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... ........
(n2)
(n2)
(n2)
a n 1 , n 1 x n 1 a n 1 , n x n C n 1
( n 1 )
( n 1 )
a n,n x n C n
Materiales y Métodos
Sistema banda de solución
Se dice que una matriz es matriz banda cuando posee
todos sus elementos no nulos concentrados en torno de
la diagonal principal. Esta es una propiedad muy
interesante que se presenta en las matrices de rigidez
locales y globales.
El método banda de resolución de sistemas de
ecuaciones lineales simultaneas se basa en esta
propiedad, permitiendo el ahorro de asignaciones de
memoria, ya que solo guarda los elementos no nulos.
Materiales y Métodos
Sistema banda de solución
El presente método centra su eficiencia en la forma en que
se realiza la numeración global de los nodos.
S em ian ch o d e b an d a
N ro. total de incognitas
X X X
0 0 0 0
X X X X 0 0 0
X X X X X 0 0
0 X X X X X 0
0 0 X X X X X
0 0 0 X X X X
0 0 0 0 X X X
0 0 0 0 0 X X
0 0 0 0 0 0 X
0 0 0 0 0 0 0
0 0 0 0 0 0 0
0 0 0 0 0 0 0
A n ch o d e b an d a
0
0
0
0
0
0
0
0
0
0
X 0
0
0
0
0
0
0
X X 0
0
0
0
0
0
0
0
X X X 0
0
0
0
0
0
0
0
0
X X X X 0
0 0 X X X
0 X X X X
A n ch o d e b an d a
D iag o n al
p rin cip al
X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X X X
X X X
0 X X X X
0 0 X X X
X X X X
0
X X X 0 0
X X
O rg an izació n co n v en cio n al d e la
m atriz d e rig id ez
0
X 0 0
O rg an izació n en b an d a p ara O rg an izació n en b an d a p ara
m atrices n o sim étricas
m atrices sim étricas
Materiales y Métodos
Método Frontal
Este método trata de evitar el ensamblaje de la matriz
de rigidez global. Consiste básicamente en ensamblar la
menor cantidad de elementos y eliminar las variables al
mismo tiempo, formando matrices de menor tamaño.
Esto es, inmediatamente después de que los coeficientes
de una ecuación se ensamblan a partir de la contribución
de todos los elementos relevantes, o que inciden sobre este
coeficiente, la variable de este puede eliminarse. Como
resultado, la matriz completa de toda la estructura jamás
se forma.
Materiales y Métodos
Método Frontal
A diferencia del método banda, la adecuada
enumeración de los elemento es de gran importancia, no
así la enumeración global de los nodos ya que estos son
tratados por su numeración local. Como primer medida
esta característica del Front lo hace mas versátil ya que
como la cantidad de elementos siempre es menor que la
cantidad de nodos, independientemente de la forma de la
estructura, siempre se tendrá mayos facilidad para
enumerar los elementos que los nodos.
Discusión de Resultados
Para analizar las cualidades de los distintos métodos
propuestos se planteó la resolución mediante el empleo del
FECCUND V 2.1 (software desarrollado por el Dto. De
Mecánica Aplicada de la Facultad de Ingeniería de la
UNNE) de dos ejemplos prácticos, representando los
resultados mediante el pre-post procesador GID 7.2.
Para dicho análisis se tomaró como variables el tiempo
de ejecución de cada método y su exactitud, la cual será
determinada por comparación con los resultados
obtenidos por el software denominado CALSEF provisto
por el Centro Internacional de Métodos Numéricos en
Ingeniería (CIMNE).
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Mallado con elementos
hexaédricos
Desplazamiento nodal
según eje y
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Tabla Nº 1
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
1,758
0,501953E-04
0,47619E-04
5,13
Banda
0,552
0,498874E-04
4,55
Frontal
0,441
0,479135E-04
0,61
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Mallado con elementos hexaédricos
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Desplazamiento nodal según eje y
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Tabla Nº 2
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
31,547
-0,00615923
-0.00544118
11,66
Banda
2,025
-0,00560842
2,98
Frontal
1,752
-0,00552264
1,47
Conclusiones
El método Frontal es superior al Gauss y al Banda
en cuanto a su exactitud y a la velocidad de
ejecución.
El efecto de la propagación de errores por redondeo
incrementa notablemente la deficiencia del método
de Gauss frente a problemas de mayor
envergadura.
El método de Banda puede ser recomendable en el
caso de contar con un software de pre-post
procesamiento ya que los mismos generan las
mallas enumerando los nodos en la forma mas
adecuada disminuyendo de ancho de banda.
Slide 10
Caracterización cualitativa de
algunos métodos de solución de
sistemas de ecuaciones lineales.
MROGINSKI, Javier L. – BENEYTO, Pablo A. –
DI RADO, H. Ariel – MANZOLILLO, Juan E.
Departamento de Mecánica Aplicada.
Facultad de Ingeniería. U.N.N.E.
Antecedentes
El Método de los Elementos Finitos (M.E.F.) es una
técnica de análisis numérico que permite la obtención de
soluciones aproximadas de ecuaciones diferenciales capaces
de representar una variedad de problemas de la Ingeniería y
Física.
El resultado de la aplicación del Método de los Elementos
Finitos a un dominio cualquier es un sistema de ecuaciones
lineales simultáneas de gran tamaño, correspondiendo una
ecuación por cada incógnita que se evalúe en cada nodo,
imposible de ser resuelto sin el auxilio de un computador.
El presente trabajo desarrolla algunos aspectos
significativos sobre los tres diferentes métodos de resolución
de sistemas de ecuaciones lineales simultáneas mas usados en
el Método de los Elementos Finitos que son el: Eliminación
de Gauss, Sistema Banda y Frontal entre otros.
Materiales y Métodos
En las diapositivas siguientes
expondrá una breve descripción de
métodos de solución de sistemas
ecuaciones lineales simultáneas que
estudian en el presente trabajo.
se
los
de
se
Materiales y Métodos
Eliminación de Gauss
Un sistema de ecuaciones lineales simultáneas
puede ser expresado en forma genérica de la
siguiente manera:
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
a 21 x 1 a 22 x 2 a 23 x 3 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... .
a n1 x 1 a n 2 x 2 a n 3 x 3 a n n x n C n
Materiales y Métodos
Eliminación de Gauss
Este método busca obtener un sistema de ecuaciones
triangular equivalente, mediante el uso de transformaciones
lineales. El cual se resuelve con facilidad aplicando la
denominada sustitución inversa.
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
´
´
´
´
a 22 x 2 a 23 x 2 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... ........
(n2)
(n2)
(n2)
a n 1 , n 1 x n 1 a n 1 , n x n C n 1
( n 1 )
( n 1 )
a n,n x n C n
Materiales y Métodos
Sistema banda de solución
Se dice que una matriz es matriz banda cuando posee
todos sus elementos no nulos concentrados en torno de
la diagonal principal. Esta es una propiedad muy
interesante que se presenta en las matrices de rigidez
locales y globales.
El método banda de resolución de sistemas de
ecuaciones lineales simultaneas se basa en esta
propiedad, permitiendo el ahorro de asignaciones de
memoria, ya que solo guarda los elementos no nulos.
Materiales y Métodos
Sistema banda de solución
El presente método centra su eficiencia en la forma en que
se realiza la numeración global de los nodos.
S em ian ch o d e b an d a
N ro. total de incognitas
X X X
0 0 0 0
X X X X 0 0 0
X X X X X 0 0
0 X X X X X 0
0 0 X X X X X
0 0 0 X X X X
0 0 0 0 X X X
0 0 0 0 0 X X
0 0 0 0 0 0 X
0 0 0 0 0 0 0
0 0 0 0 0 0 0
0 0 0 0 0 0 0
A n ch o d e b an d a
0
0
0
0
0
0
0
0
0
0
X 0
0
0
0
0
0
0
X X 0
0
0
0
0
0
0
0
X X X 0
0
0
0
0
0
0
0
0
X X X X 0
0 0 X X X
0 X X X X
A n ch o d e b an d a
D iag o n al
p rin cip al
X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X X X
X X X
0 X X X X
0 0 X X X
X X X X
0
X X X 0 0
X X
O rg an izació n co n v en cio n al d e la
m atriz d e rig id ez
0
X 0 0
O rg an izació n en b an d a p ara O rg an izació n en b an d a p ara
m atrices n o sim étricas
m atrices sim étricas
Materiales y Métodos
Método Frontal
Este método trata de evitar el ensamblaje de la matriz
de rigidez global. Consiste básicamente en ensamblar la
menor cantidad de elementos y eliminar las variables al
mismo tiempo, formando matrices de menor tamaño.
Esto es, inmediatamente después de que los coeficientes
de una ecuación se ensamblan a partir de la contribución
de todos los elementos relevantes, o que inciden sobre este
coeficiente, la variable de este puede eliminarse. Como
resultado, la matriz completa de toda la estructura jamás
se forma.
Materiales y Métodos
Método Frontal
A diferencia del método banda, la adecuada
enumeración de los elemento es de gran importancia, no
así la enumeración global de los nodos ya que estos son
tratados por su numeración local. Como primer medida
esta característica del Front lo hace mas versátil ya que
como la cantidad de elementos siempre es menor que la
cantidad de nodos, independientemente de la forma de la
estructura, siempre se tendrá mayos facilidad para
enumerar los elementos que los nodos.
Discusión de Resultados
Para analizar las cualidades de los distintos métodos
propuestos se planteó la resolución mediante el empleo del
FECCUND V 2.1 (software desarrollado por el Dto. De
Mecánica Aplicada de la Facultad de Ingeniería de la
UNNE) de dos ejemplos prácticos, representando los
resultados mediante el pre-post procesador GID 7.2.
Para dicho análisis se tomaró como variables el tiempo
de ejecución de cada método y su exactitud, la cual será
determinada por comparación con los resultados
obtenidos por el software denominado CALSEF provisto
por el Centro Internacional de Métodos Numéricos en
Ingeniería (CIMNE).
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Mallado con elementos
hexaédricos
Desplazamiento nodal
según eje y
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Tabla Nº 1
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
1,758
0,501953E-04
0,47619E-04
5,13
Banda
0,552
0,498874E-04
4,55
Frontal
0,441
0,479135E-04
0,61
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Mallado con elementos hexaédricos
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Desplazamiento nodal según eje y
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Tabla Nº 2
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
31,547
-0,00615923
-0.00544118
11,66
Banda
2,025
-0,00560842
2,98
Frontal
1,752
-0,00552264
1,47
Conclusiones
El método Frontal es superior al Gauss y al Banda
en cuanto a su exactitud y a la velocidad de
ejecución.
El efecto de la propagación de errores por redondeo
incrementa notablemente la deficiencia del método
de Gauss frente a problemas de mayor
envergadura.
El método de Banda puede ser recomendable en el
caso de contar con un software de pre-post
procesamiento ya que los mismos generan las
mallas enumerando los nodos en la forma mas
adecuada disminuyendo de ancho de banda.
Slide 11
Caracterización cualitativa de
algunos métodos de solución de
sistemas de ecuaciones lineales.
MROGINSKI, Javier L. – BENEYTO, Pablo A. –
DI RADO, H. Ariel – MANZOLILLO, Juan E.
Departamento de Mecánica Aplicada.
Facultad de Ingeniería. U.N.N.E.
Antecedentes
El Método de los Elementos Finitos (M.E.F.) es una
técnica de análisis numérico que permite la obtención de
soluciones aproximadas de ecuaciones diferenciales capaces
de representar una variedad de problemas de la Ingeniería y
Física.
El resultado de la aplicación del Método de los Elementos
Finitos a un dominio cualquier es un sistema de ecuaciones
lineales simultáneas de gran tamaño, correspondiendo una
ecuación por cada incógnita que se evalúe en cada nodo,
imposible de ser resuelto sin el auxilio de un computador.
El presente trabajo desarrolla algunos aspectos
significativos sobre los tres diferentes métodos de resolución
de sistemas de ecuaciones lineales simultáneas mas usados en
el Método de los Elementos Finitos que son el: Eliminación
de Gauss, Sistema Banda y Frontal entre otros.
Materiales y Métodos
En las diapositivas siguientes
expondrá una breve descripción de
métodos de solución de sistemas
ecuaciones lineales simultáneas que
estudian en el presente trabajo.
se
los
de
se
Materiales y Métodos
Eliminación de Gauss
Un sistema de ecuaciones lineales simultáneas
puede ser expresado en forma genérica de la
siguiente manera:
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
a 21 x 1 a 22 x 2 a 23 x 3 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... .
a n1 x 1 a n 2 x 2 a n 3 x 3 a n n x n C n
Materiales y Métodos
Eliminación de Gauss
Este método busca obtener un sistema de ecuaciones
triangular equivalente, mediante el uso de transformaciones
lineales. El cual se resuelve con facilidad aplicando la
denominada sustitución inversa.
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
´
´
´
´
a 22 x 2 a 23 x 2 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... ........
(n2)
(n2)
(n2)
a n 1 , n 1 x n 1 a n 1 , n x n C n 1
( n 1 )
( n 1 )
a n,n x n C n
Materiales y Métodos
Sistema banda de solución
Se dice que una matriz es matriz banda cuando posee
todos sus elementos no nulos concentrados en torno de
la diagonal principal. Esta es una propiedad muy
interesante que se presenta en las matrices de rigidez
locales y globales.
El método banda de resolución de sistemas de
ecuaciones lineales simultaneas se basa en esta
propiedad, permitiendo el ahorro de asignaciones de
memoria, ya que solo guarda los elementos no nulos.
Materiales y Métodos
Sistema banda de solución
El presente método centra su eficiencia en la forma en que
se realiza la numeración global de los nodos.
S em ian ch o d e b an d a
N ro. total de incognitas
X X X
0 0 0 0
X X X X 0 0 0
X X X X X 0 0
0 X X X X X 0
0 0 X X X X X
0 0 0 X X X X
0 0 0 0 X X X
0 0 0 0 0 X X
0 0 0 0 0 0 X
0 0 0 0 0 0 0
0 0 0 0 0 0 0
0 0 0 0 0 0 0
A n ch o d e b an d a
0
0
0
0
0
0
0
0
0
0
X 0
0
0
0
0
0
0
X X 0
0
0
0
0
0
0
0
X X X 0
0
0
0
0
0
0
0
0
X X X X 0
0 0 X X X
0 X X X X
A n ch o d e b an d a
D iag o n al
p rin cip al
X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X X X
X X X
0 X X X X
0 0 X X X
X X X X
0
X X X 0 0
X X
O rg an izació n co n v en cio n al d e la
m atriz d e rig id ez
0
X 0 0
O rg an izació n en b an d a p ara O rg an izació n en b an d a p ara
m atrices n o sim étricas
m atrices sim étricas
Materiales y Métodos
Método Frontal
Este método trata de evitar el ensamblaje de la matriz
de rigidez global. Consiste básicamente en ensamblar la
menor cantidad de elementos y eliminar las variables al
mismo tiempo, formando matrices de menor tamaño.
Esto es, inmediatamente después de que los coeficientes
de una ecuación se ensamblan a partir de la contribución
de todos los elementos relevantes, o que inciden sobre este
coeficiente, la variable de este puede eliminarse. Como
resultado, la matriz completa de toda la estructura jamás
se forma.
Materiales y Métodos
Método Frontal
A diferencia del método banda, la adecuada
enumeración de los elemento es de gran importancia, no
así la enumeración global de los nodos ya que estos son
tratados por su numeración local. Como primer medida
esta característica del Front lo hace mas versátil ya que
como la cantidad de elementos siempre es menor que la
cantidad de nodos, independientemente de la forma de la
estructura, siempre se tendrá mayos facilidad para
enumerar los elementos que los nodos.
Discusión de Resultados
Para analizar las cualidades de los distintos métodos
propuestos se planteó la resolución mediante el empleo del
FECCUND V 2.1 (software desarrollado por el Dto. De
Mecánica Aplicada de la Facultad de Ingeniería de la
UNNE) de dos ejemplos prácticos, representando los
resultados mediante el pre-post procesador GID 7.2.
Para dicho análisis se tomaró como variables el tiempo
de ejecución de cada método y su exactitud, la cual será
determinada por comparación con los resultados
obtenidos por el software denominado CALSEF provisto
por el Centro Internacional de Métodos Numéricos en
Ingeniería (CIMNE).
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Mallado con elementos
hexaédricos
Desplazamiento nodal
según eje y
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Tabla Nº 1
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
1,758
0,501953E-04
0,47619E-04
5,13
Banda
0,552
0,498874E-04
4,55
Frontal
0,441
0,479135E-04
0,61
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Mallado con elementos hexaédricos
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Desplazamiento nodal según eje y
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Tabla Nº 2
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
31,547
-0,00615923
-0.00544118
11,66
Banda
2,025
-0,00560842
2,98
Frontal
1,752
-0,00552264
1,47
Conclusiones
El método Frontal es superior al Gauss y al Banda
en cuanto a su exactitud y a la velocidad de
ejecución.
El efecto de la propagación de errores por redondeo
incrementa notablemente la deficiencia del método
de Gauss frente a problemas de mayor
envergadura.
El método de Banda puede ser recomendable en el
caso de contar con un software de pre-post
procesamiento ya que los mismos generan las
mallas enumerando los nodos en la forma mas
adecuada disminuyendo de ancho de banda.
Slide 12
Caracterización cualitativa de
algunos métodos de solución de
sistemas de ecuaciones lineales.
MROGINSKI, Javier L. – BENEYTO, Pablo A. –
DI RADO, H. Ariel – MANZOLILLO, Juan E.
Departamento de Mecánica Aplicada.
Facultad de Ingeniería. U.N.N.E.
Antecedentes
El Método de los Elementos Finitos (M.E.F.) es una
técnica de análisis numérico que permite la obtención de
soluciones aproximadas de ecuaciones diferenciales capaces
de representar una variedad de problemas de la Ingeniería y
Física.
El resultado de la aplicación del Método de los Elementos
Finitos a un dominio cualquier es un sistema de ecuaciones
lineales simultáneas de gran tamaño, correspondiendo una
ecuación por cada incógnita que se evalúe en cada nodo,
imposible de ser resuelto sin el auxilio de un computador.
El presente trabajo desarrolla algunos aspectos
significativos sobre los tres diferentes métodos de resolución
de sistemas de ecuaciones lineales simultáneas mas usados en
el Método de los Elementos Finitos que son el: Eliminación
de Gauss, Sistema Banda y Frontal entre otros.
Materiales y Métodos
En las diapositivas siguientes
expondrá una breve descripción de
métodos de solución de sistemas
ecuaciones lineales simultáneas que
estudian en el presente trabajo.
se
los
de
se
Materiales y Métodos
Eliminación de Gauss
Un sistema de ecuaciones lineales simultáneas
puede ser expresado en forma genérica de la
siguiente manera:
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
a 21 x 1 a 22 x 2 a 23 x 3 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... .
a n1 x 1 a n 2 x 2 a n 3 x 3 a n n x n C n
Materiales y Métodos
Eliminación de Gauss
Este método busca obtener un sistema de ecuaciones
triangular equivalente, mediante el uso de transformaciones
lineales. El cual se resuelve con facilidad aplicando la
denominada sustitución inversa.
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
´
´
´
´
a 22 x 2 a 23 x 2 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... ........
(n2)
(n2)
(n2)
a n 1 , n 1 x n 1 a n 1 , n x n C n 1
( n 1 )
( n 1 )
a n,n x n C n
Materiales y Métodos
Sistema banda de solución
Se dice que una matriz es matriz banda cuando posee
todos sus elementos no nulos concentrados en torno de
la diagonal principal. Esta es una propiedad muy
interesante que se presenta en las matrices de rigidez
locales y globales.
El método banda de resolución de sistemas de
ecuaciones lineales simultaneas se basa en esta
propiedad, permitiendo el ahorro de asignaciones de
memoria, ya que solo guarda los elementos no nulos.
Materiales y Métodos
Sistema banda de solución
El presente método centra su eficiencia en la forma en que
se realiza la numeración global de los nodos.
S em ian ch o d e b an d a
N ro. total de incognitas
X X X
0 0 0 0
X X X X 0 0 0
X X X X X 0 0
0 X X X X X 0
0 0 X X X X X
0 0 0 X X X X
0 0 0 0 X X X
0 0 0 0 0 X X
0 0 0 0 0 0 X
0 0 0 0 0 0 0
0 0 0 0 0 0 0
0 0 0 0 0 0 0
A n ch o d e b an d a
0
0
0
0
0
0
0
0
0
0
X 0
0
0
0
0
0
0
X X 0
0
0
0
0
0
0
0
X X X 0
0
0
0
0
0
0
0
0
X X X X 0
0 0 X X X
0 X X X X
A n ch o d e b an d a
D iag o n al
p rin cip al
X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X X X
X X X
0 X X X X
0 0 X X X
X X X X
0
X X X 0 0
X X
O rg an izació n co n v en cio n al d e la
m atriz d e rig id ez
0
X 0 0
O rg an izació n en b an d a p ara O rg an izació n en b an d a p ara
m atrices n o sim étricas
m atrices sim étricas
Materiales y Métodos
Método Frontal
Este método trata de evitar el ensamblaje de la matriz
de rigidez global. Consiste básicamente en ensamblar la
menor cantidad de elementos y eliminar las variables al
mismo tiempo, formando matrices de menor tamaño.
Esto es, inmediatamente después de que los coeficientes
de una ecuación se ensamblan a partir de la contribución
de todos los elementos relevantes, o que inciden sobre este
coeficiente, la variable de este puede eliminarse. Como
resultado, la matriz completa de toda la estructura jamás
se forma.
Materiales y Métodos
Método Frontal
A diferencia del método banda, la adecuada
enumeración de los elemento es de gran importancia, no
así la enumeración global de los nodos ya que estos son
tratados por su numeración local. Como primer medida
esta característica del Front lo hace mas versátil ya que
como la cantidad de elementos siempre es menor que la
cantidad de nodos, independientemente de la forma de la
estructura, siempre se tendrá mayos facilidad para
enumerar los elementos que los nodos.
Discusión de Resultados
Para analizar las cualidades de los distintos métodos
propuestos se planteó la resolución mediante el empleo del
FECCUND V 2.1 (software desarrollado por el Dto. De
Mecánica Aplicada de la Facultad de Ingeniería de la
UNNE) de dos ejemplos prácticos, representando los
resultados mediante el pre-post procesador GID 7.2.
Para dicho análisis se tomaró como variables el tiempo
de ejecución de cada método y su exactitud, la cual será
determinada por comparación con los resultados
obtenidos por el software denominado CALSEF provisto
por el Centro Internacional de Métodos Numéricos en
Ingeniería (CIMNE).
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Mallado con elementos
hexaédricos
Desplazamiento nodal
según eje y
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Tabla Nº 1
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
1,758
0,501953E-04
0,47619E-04
5,13
Banda
0,552
0,498874E-04
4,55
Frontal
0,441
0,479135E-04
0,61
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Mallado con elementos hexaédricos
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Desplazamiento nodal según eje y
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Tabla Nº 2
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
31,547
-0,00615923
-0.00544118
11,66
Banda
2,025
-0,00560842
2,98
Frontal
1,752
-0,00552264
1,47
Conclusiones
El método Frontal es superior al Gauss y al Banda
en cuanto a su exactitud y a la velocidad de
ejecución.
El efecto de la propagación de errores por redondeo
incrementa notablemente la deficiencia del método
de Gauss frente a problemas de mayor
envergadura.
El método de Banda puede ser recomendable en el
caso de contar con un software de pre-post
procesamiento ya que los mismos generan las
mallas enumerando los nodos en la forma mas
adecuada disminuyendo de ancho de banda.
Slide 13
Caracterización cualitativa de
algunos métodos de solución de
sistemas de ecuaciones lineales.
MROGINSKI, Javier L. – BENEYTO, Pablo A. –
DI RADO, H. Ariel – MANZOLILLO, Juan E.
Departamento de Mecánica Aplicada.
Facultad de Ingeniería. U.N.N.E.
Antecedentes
El Método de los Elementos Finitos (M.E.F.) es una
técnica de análisis numérico que permite la obtención de
soluciones aproximadas de ecuaciones diferenciales capaces
de representar una variedad de problemas de la Ingeniería y
Física.
El resultado de la aplicación del Método de los Elementos
Finitos a un dominio cualquier es un sistema de ecuaciones
lineales simultáneas de gran tamaño, correspondiendo una
ecuación por cada incógnita que se evalúe en cada nodo,
imposible de ser resuelto sin el auxilio de un computador.
El presente trabajo desarrolla algunos aspectos
significativos sobre los tres diferentes métodos de resolución
de sistemas de ecuaciones lineales simultáneas mas usados en
el Método de los Elementos Finitos que son el: Eliminación
de Gauss, Sistema Banda y Frontal entre otros.
Materiales y Métodos
En las diapositivas siguientes
expondrá una breve descripción de
métodos de solución de sistemas
ecuaciones lineales simultáneas que
estudian en el presente trabajo.
se
los
de
se
Materiales y Métodos
Eliminación de Gauss
Un sistema de ecuaciones lineales simultáneas
puede ser expresado en forma genérica de la
siguiente manera:
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
a 21 x 1 a 22 x 2 a 23 x 3 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... .
a n1 x 1 a n 2 x 2 a n 3 x 3 a n n x n C n
Materiales y Métodos
Eliminación de Gauss
Este método busca obtener un sistema de ecuaciones
triangular equivalente, mediante el uso de transformaciones
lineales. El cual se resuelve con facilidad aplicando la
denominada sustitución inversa.
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
´
´
´
´
a 22 x 2 a 23 x 2 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... ........
(n2)
(n2)
(n2)
a n 1 , n 1 x n 1 a n 1 , n x n C n 1
( n 1 )
( n 1 )
a n,n x n C n
Materiales y Métodos
Sistema banda de solución
Se dice que una matriz es matriz banda cuando posee
todos sus elementos no nulos concentrados en torno de
la diagonal principal. Esta es una propiedad muy
interesante que se presenta en las matrices de rigidez
locales y globales.
El método banda de resolución de sistemas de
ecuaciones lineales simultaneas se basa en esta
propiedad, permitiendo el ahorro de asignaciones de
memoria, ya que solo guarda los elementos no nulos.
Materiales y Métodos
Sistema banda de solución
El presente método centra su eficiencia en la forma en que
se realiza la numeración global de los nodos.
S em ian ch o d e b an d a
N ro. total de incognitas
X X X
0 0 0 0
X X X X 0 0 0
X X X X X 0 0
0 X X X X X 0
0 0 X X X X X
0 0 0 X X X X
0 0 0 0 X X X
0 0 0 0 0 X X
0 0 0 0 0 0 X
0 0 0 0 0 0 0
0 0 0 0 0 0 0
0 0 0 0 0 0 0
A n ch o d e b an d a
0
0
0
0
0
0
0
0
0
0
X 0
0
0
0
0
0
0
X X 0
0
0
0
0
0
0
0
X X X 0
0
0
0
0
0
0
0
0
X X X X 0
0 0 X X X
0 X X X X
A n ch o d e b an d a
D iag o n al
p rin cip al
X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X X X
X X X
0 X X X X
0 0 X X X
X X X X
0
X X X 0 0
X X
O rg an izació n co n v en cio n al d e la
m atriz d e rig id ez
0
X 0 0
O rg an izació n en b an d a p ara O rg an izació n en b an d a p ara
m atrices n o sim étricas
m atrices sim étricas
Materiales y Métodos
Método Frontal
Este método trata de evitar el ensamblaje de la matriz
de rigidez global. Consiste básicamente en ensamblar la
menor cantidad de elementos y eliminar las variables al
mismo tiempo, formando matrices de menor tamaño.
Esto es, inmediatamente después de que los coeficientes
de una ecuación se ensamblan a partir de la contribución
de todos los elementos relevantes, o que inciden sobre este
coeficiente, la variable de este puede eliminarse. Como
resultado, la matriz completa de toda la estructura jamás
se forma.
Materiales y Métodos
Método Frontal
A diferencia del método banda, la adecuada
enumeración de los elemento es de gran importancia, no
así la enumeración global de los nodos ya que estos son
tratados por su numeración local. Como primer medida
esta característica del Front lo hace mas versátil ya que
como la cantidad de elementos siempre es menor que la
cantidad de nodos, independientemente de la forma de la
estructura, siempre se tendrá mayos facilidad para
enumerar los elementos que los nodos.
Discusión de Resultados
Para analizar las cualidades de los distintos métodos
propuestos se planteó la resolución mediante el empleo del
FECCUND V 2.1 (software desarrollado por el Dto. De
Mecánica Aplicada de la Facultad de Ingeniería de la
UNNE) de dos ejemplos prácticos, representando los
resultados mediante el pre-post procesador GID 7.2.
Para dicho análisis se tomaró como variables el tiempo
de ejecución de cada método y su exactitud, la cual será
determinada por comparación con los resultados
obtenidos por el software denominado CALSEF provisto
por el Centro Internacional de Métodos Numéricos en
Ingeniería (CIMNE).
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Mallado con elementos
hexaédricos
Desplazamiento nodal
según eje y
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Tabla Nº 1
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
1,758
0,501953E-04
0,47619E-04
5,13
Banda
0,552
0,498874E-04
4,55
Frontal
0,441
0,479135E-04
0,61
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Mallado con elementos hexaédricos
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Desplazamiento nodal según eje y
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Tabla Nº 2
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
31,547
-0,00615923
-0.00544118
11,66
Banda
2,025
-0,00560842
2,98
Frontal
1,752
-0,00552264
1,47
Conclusiones
El método Frontal es superior al Gauss y al Banda
en cuanto a su exactitud y a la velocidad de
ejecución.
El efecto de la propagación de errores por redondeo
incrementa notablemente la deficiencia del método
de Gauss frente a problemas de mayor
envergadura.
El método de Banda puede ser recomendable en el
caso de contar con un software de pre-post
procesamiento ya que los mismos generan las
mallas enumerando los nodos en la forma mas
adecuada disminuyendo de ancho de banda.
Slide 14
Caracterización cualitativa de
algunos métodos de solución de
sistemas de ecuaciones lineales.
MROGINSKI, Javier L. – BENEYTO, Pablo A. –
DI RADO, H. Ariel – MANZOLILLO, Juan E.
Departamento de Mecánica Aplicada.
Facultad de Ingeniería. U.N.N.E.
Antecedentes
El Método de los Elementos Finitos (M.E.F.) es una
técnica de análisis numérico que permite la obtención de
soluciones aproximadas de ecuaciones diferenciales capaces
de representar una variedad de problemas de la Ingeniería y
Física.
El resultado de la aplicación del Método de los Elementos
Finitos a un dominio cualquier es un sistema de ecuaciones
lineales simultáneas de gran tamaño, correspondiendo una
ecuación por cada incógnita que se evalúe en cada nodo,
imposible de ser resuelto sin el auxilio de un computador.
El presente trabajo desarrolla algunos aspectos
significativos sobre los tres diferentes métodos de resolución
de sistemas de ecuaciones lineales simultáneas mas usados en
el Método de los Elementos Finitos que son el: Eliminación
de Gauss, Sistema Banda y Frontal entre otros.
Materiales y Métodos
En las diapositivas siguientes
expondrá una breve descripción de
métodos de solución de sistemas
ecuaciones lineales simultáneas que
estudian en el presente trabajo.
se
los
de
se
Materiales y Métodos
Eliminación de Gauss
Un sistema de ecuaciones lineales simultáneas
puede ser expresado en forma genérica de la
siguiente manera:
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
a 21 x 1 a 22 x 2 a 23 x 3 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... .
a n1 x 1 a n 2 x 2 a n 3 x 3 a n n x n C n
Materiales y Métodos
Eliminación de Gauss
Este método busca obtener un sistema de ecuaciones
triangular equivalente, mediante el uso de transformaciones
lineales. El cual se resuelve con facilidad aplicando la
denominada sustitución inversa.
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
´
´
´
´
a 22 x 2 a 23 x 2 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... ........
(n2)
(n2)
(n2)
a n 1 , n 1 x n 1 a n 1 , n x n C n 1
( n 1 )
( n 1 )
a n,n x n C n
Materiales y Métodos
Sistema banda de solución
Se dice que una matriz es matriz banda cuando posee
todos sus elementos no nulos concentrados en torno de
la diagonal principal. Esta es una propiedad muy
interesante que se presenta en las matrices de rigidez
locales y globales.
El método banda de resolución de sistemas de
ecuaciones lineales simultaneas se basa en esta
propiedad, permitiendo el ahorro de asignaciones de
memoria, ya que solo guarda los elementos no nulos.
Materiales y Métodos
Sistema banda de solución
El presente método centra su eficiencia en la forma en que
se realiza la numeración global de los nodos.
S em ian ch o d e b an d a
N ro. total de incognitas
X X X
0 0 0 0
X X X X 0 0 0
X X X X X 0 0
0 X X X X X 0
0 0 X X X X X
0 0 0 X X X X
0 0 0 0 X X X
0 0 0 0 0 X X
0 0 0 0 0 0 X
0 0 0 0 0 0 0
0 0 0 0 0 0 0
0 0 0 0 0 0 0
A n ch o d e b an d a
0
0
0
0
0
0
0
0
0
0
X 0
0
0
0
0
0
0
X X 0
0
0
0
0
0
0
0
X X X 0
0
0
0
0
0
0
0
0
X X X X 0
0 0 X X X
0 X X X X
A n ch o d e b an d a
D iag o n al
p rin cip al
X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X X X
X X X
0 X X X X
0 0 X X X
X X X X
0
X X X 0 0
X X
O rg an izació n co n v en cio n al d e la
m atriz d e rig id ez
0
X 0 0
O rg an izació n en b an d a p ara O rg an izació n en b an d a p ara
m atrices n o sim étricas
m atrices sim étricas
Materiales y Métodos
Método Frontal
Este método trata de evitar el ensamblaje de la matriz
de rigidez global. Consiste básicamente en ensamblar la
menor cantidad de elementos y eliminar las variables al
mismo tiempo, formando matrices de menor tamaño.
Esto es, inmediatamente después de que los coeficientes
de una ecuación se ensamblan a partir de la contribución
de todos los elementos relevantes, o que inciden sobre este
coeficiente, la variable de este puede eliminarse. Como
resultado, la matriz completa de toda la estructura jamás
se forma.
Materiales y Métodos
Método Frontal
A diferencia del método banda, la adecuada
enumeración de los elemento es de gran importancia, no
así la enumeración global de los nodos ya que estos son
tratados por su numeración local. Como primer medida
esta característica del Front lo hace mas versátil ya que
como la cantidad de elementos siempre es menor que la
cantidad de nodos, independientemente de la forma de la
estructura, siempre se tendrá mayos facilidad para
enumerar los elementos que los nodos.
Discusión de Resultados
Para analizar las cualidades de los distintos métodos
propuestos se planteó la resolución mediante el empleo del
FECCUND V 2.1 (software desarrollado por el Dto. De
Mecánica Aplicada de la Facultad de Ingeniería de la
UNNE) de dos ejemplos prácticos, representando los
resultados mediante el pre-post procesador GID 7.2.
Para dicho análisis se tomaró como variables el tiempo
de ejecución de cada método y su exactitud, la cual será
determinada por comparación con los resultados
obtenidos por el software denominado CALSEF provisto
por el Centro Internacional de Métodos Numéricos en
Ingeniería (CIMNE).
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Mallado con elementos
hexaédricos
Desplazamiento nodal
según eje y
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Tabla Nº 1
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
1,758
0,501953E-04
0,47619E-04
5,13
Banda
0,552
0,498874E-04
4,55
Frontal
0,441
0,479135E-04
0,61
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Mallado con elementos hexaédricos
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Desplazamiento nodal según eje y
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Tabla Nº 2
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
31,547
-0,00615923
-0.00544118
11,66
Banda
2,025
-0,00560842
2,98
Frontal
1,752
-0,00552264
1,47
Conclusiones
El método Frontal es superior al Gauss y al Banda
en cuanto a su exactitud y a la velocidad de
ejecución.
El efecto de la propagación de errores por redondeo
incrementa notablemente la deficiencia del método
de Gauss frente a problemas de mayor
envergadura.
El método de Banda puede ser recomendable en el
caso de contar con un software de pre-post
procesamiento ya que los mismos generan las
mallas enumerando los nodos en la forma mas
adecuada disminuyendo de ancho de banda.
Slide 15
Caracterización cualitativa de
algunos métodos de solución de
sistemas de ecuaciones lineales.
MROGINSKI, Javier L. – BENEYTO, Pablo A. –
DI RADO, H. Ariel – MANZOLILLO, Juan E.
Departamento de Mecánica Aplicada.
Facultad de Ingeniería. U.N.N.E.
Antecedentes
El Método de los Elementos Finitos (M.E.F.) es una
técnica de análisis numérico que permite la obtención de
soluciones aproximadas de ecuaciones diferenciales capaces
de representar una variedad de problemas de la Ingeniería y
Física.
El resultado de la aplicación del Método de los Elementos
Finitos a un dominio cualquier es un sistema de ecuaciones
lineales simultáneas de gran tamaño, correspondiendo una
ecuación por cada incógnita que se evalúe en cada nodo,
imposible de ser resuelto sin el auxilio de un computador.
El presente trabajo desarrolla algunos aspectos
significativos sobre los tres diferentes métodos de resolución
de sistemas de ecuaciones lineales simultáneas mas usados en
el Método de los Elementos Finitos que son el: Eliminación
de Gauss, Sistema Banda y Frontal entre otros.
Materiales y Métodos
En las diapositivas siguientes
expondrá una breve descripción de
métodos de solución de sistemas
ecuaciones lineales simultáneas que
estudian en el presente trabajo.
se
los
de
se
Materiales y Métodos
Eliminación de Gauss
Un sistema de ecuaciones lineales simultáneas
puede ser expresado en forma genérica de la
siguiente manera:
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
a 21 x 1 a 22 x 2 a 23 x 3 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... .
a n1 x 1 a n 2 x 2 a n 3 x 3 a n n x n C n
Materiales y Métodos
Eliminación de Gauss
Este método busca obtener un sistema de ecuaciones
triangular equivalente, mediante el uso de transformaciones
lineales. El cual se resuelve con facilidad aplicando la
denominada sustitución inversa.
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
´
´
´
´
a 22 x 2 a 23 x 2 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... ........
(n2)
(n2)
(n2)
a n 1 , n 1 x n 1 a n 1 , n x n C n 1
( n 1 )
( n 1 )
a n,n x n C n
Materiales y Métodos
Sistema banda de solución
Se dice que una matriz es matriz banda cuando posee
todos sus elementos no nulos concentrados en torno de
la diagonal principal. Esta es una propiedad muy
interesante que se presenta en las matrices de rigidez
locales y globales.
El método banda de resolución de sistemas de
ecuaciones lineales simultaneas se basa en esta
propiedad, permitiendo el ahorro de asignaciones de
memoria, ya que solo guarda los elementos no nulos.
Materiales y Métodos
Sistema banda de solución
El presente método centra su eficiencia en la forma en que
se realiza la numeración global de los nodos.
S em ian ch o d e b an d a
N ro. total de incognitas
X X X
0 0 0 0
X X X X 0 0 0
X X X X X 0 0
0 X X X X X 0
0 0 X X X X X
0 0 0 X X X X
0 0 0 0 X X X
0 0 0 0 0 X X
0 0 0 0 0 0 X
0 0 0 0 0 0 0
0 0 0 0 0 0 0
0 0 0 0 0 0 0
A n ch o d e b an d a
0
0
0
0
0
0
0
0
0
0
X 0
0
0
0
0
0
0
X X 0
0
0
0
0
0
0
0
X X X 0
0
0
0
0
0
0
0
0
X X X X 0
0 0 X X X
0 X X X X
A n ch o d e b an d a
D iag o n al
p rin cip al
X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X X X
X X X
0 X X X X
0 0 X X X
X X X X
0
X X X 0 0
X X
O rg an izació n co n v en cio n al d e la
m atriz d e rig id ez
0
X 0 0
O rg an izació n en b an d a p ara O rg an izació n en b an d a p ara
m atrices n o sim étricas
m atrices sim étricas
Materiales y Métodos
Método Frontal
Este método trata de evitar el ensamblaje de la matriz
de rigidez global. Consiste básicamente en ensamblar la
menor cantidad de elementos y eliminar las variables al
mismo tiempo, formando matrices de menor tamaño.
Esto es, inmediatamente después de que los coeficientes
de una ecuación se ensamblan a partir de la contribución
de todos los elementos relevantes, o que inciden sobre este
coeficiente, la variable de este puede eliminarse. Como
resultado, la matriz completa de toda la estructura jamás
se forma.
Materiales y Métodos
Método Frontal
A diferencia del método banda, la adecuada
enumeración de los elemento es de gran importancia, no
así la enumeración global de los nodos ya que estos son
tratados por su numeración local. Como primer medida
esta característica del Front lo hace mas versátil ya que
como la cantidad de elementos siempre es menor que la
cantidad de nodos, independientemente de la forma de la
estructura, siempre se tendrá mayos facilidad para
enumerar los elementos que los nodos.
Discusión de Resultados
Para analizar las cualidades de los distintos métodos
propuestos se planteó la resolución mediante el empleo del
FECCUND V 2.1 (software desarrollado por el Dto. De
Mecánica Aplicada de la Facultad de Ingeniería de la
UNNE) de dos ejemplos prácticos, representando los
resultados mediante el pre-post procesador GID 7.2.
Para dicho análisis se tomaró como variables el tiempo
de ejecución de cada método y su exactitud, la cual será
determinada por comparación con los resultados
obtenidos por el software denominado CALSEF provisto
por el Centro Internacional de Métodos Numéricos en
Ingeniería (CIMNE).
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Mallado con elementos
hexaédricos
Desplazamiento nodal
según eje y
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Tabla Nº 1
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
1,758
0,501953E-04
0,47619E-04
5,13
Banda
0,552
0,498874E-04
4,55
Frontal
0,441
0,479135E-04
0,61
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Mallado con elementos hexaédricos
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Desplazamiento nodal según eje y
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Tabla Nº 2
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
31,547
-0,00615923
-0.00544118
11,66
Banda
2,025
-0,00560842
2,98
Frontal
1,752
-0,00552264
1,47
Conclusiones
El método Frontal es superior al Gauss y al Banda
en cuanto a su exactitud y a la velocidad de
ejecución.
El efecto de la propagación de errores por redondeo
incrementa notablemente la deficiencia del método
de Gauss frente a problemas de mayor
envergadura.
El método de Banda puede ser recomendable en el
caso de contar con un software de pre-post
procesamiento ya que los mismos generan las
mallas enumerando los nodos en la forma mas
adecuada disminuyendo de ancho de banda.
Slide 16
Caracterización cualitativa de
algunos métodos de solución de
sistemas de ecuaciones lineales.
MROGINSKI, Javier L. – BENEYTO, Pablo A. –
DI RADO, H. Ariel – MANZOLILLO, Juan E.
Departamento de Mecánica Aplicada.
Facultad de Ingeniería. U.N.N.E.
Antecedentes
El Método de los Elementos Finitos (M.E.F.) es una
técnica de análisis numérico que permite la obtención de
soluciones aproximadas de ecuaciones diferenciales capaces
de representar una variedad de problemas de la Ingeniería y
Física.
El resultado de la aplicación del Método de los Elementos
Finitos a un dominio cualquier es un sistema de ecuaciones
lineales simultáneas de gran tamaño, correspondiendo una
ecuación por cada incógnita que se evalúe en cada nodo,
imposible de ser resuelto sin el auxilio de un computador.
El presente trabajo desarrolla algunos aspectos
significativos sobre los tres diferentes métodos de resolución
de sistemas de ecuaciones lineales simultáneas mas usados en
el Método de los Elementos Finitos que son el: Eliminación
de Gauss, Sistema Banda y Frontal entre otros.
Materiales y Métodos
En las diapositivas siguientes
expondrá una breve descripción de
métodos de solución de sistemas
ecuaciones lineales simultáneas que
estudian en el presente trabajo.
se
los
de
se
Materiales y Métodos
Eliminación de Gauss
Un sistema de ecuaciones lineales simultáneas
puede ser expresado en forma genérica de la
siguiente manera:
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
a 21 x 1 a 22 x 2 a 23 x 3 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... .
a n1 x 1 a n 2 x 2 a n 3 x 3 a n n x n C n
Materiales y Métodos
Eliminación de Gauss
Este método busca obtener un sistema de ecuaciones
triangular equivalente, mediante el uso de transformaciones
lineales. El cual se resuelve con facilidad aplicando la
denominada sustitución inversa.
a 11 x 1 a 12 x 2 a 13 x 3 a 1 n x n C 1
´
´
´
´
a 22 x 2 a 23 x 2 a 2 n x n C 2
.......... .......... .......... .......... .......... .......... ........
(n2)
(n2)
(n2)
a n 1 , n 1 x n 1 a n 1 , n x n C n 1
( n 1 )
( n 1 )
a n,n x n C n
Materiales y Métodos
Sistema banda de solución
Se dice que una matriz es matriz banda cuando posee
todos sus elementos no nulos concentrados en torno de
la diagonal principal. Esta es una propiedad muy
interesante que se presenta en las matrices de rigidez
locales y globales.
El método banda de resolución de sistemas de
ecuaciones lineales simultaneas se basa en esta
propiedad, permitiendo el ahorro de asignaciones de
memoria, ya que solo guarda los elementos no nulos.
Materiales y Métodos
Sistema banda de solución
El presente método centra su eficiencia en la forma en que
se realiza la numeración global de los nodos.
S em ian ch o d e b an d a
N ro. total de incognitas
X X X
0 0 0 0
X X X X 0 0 0
X X X X X 0 0
0 X X X X X 0
0 0 X X X X X
0 0 0 X X X X
0 0 0 0 X X X
0 0 0 0 0 X X
0 0 0 0 0 0 X
0 0 0 0 0 0 0
0 0 0 0 0 0 0
0 0 0 0 0 0 0
A n ch o d e b an d a
0
0
0
0
0
0
0
0
0
0
X 0
0
0
0
0
0
0
X X 0
0
0
0
0
0
0
0
X X X 0
0
0
0
0
0
0
0
0
X X X X 0
0 0 X X X
0 X X X X
A n ch o d e b an d a
D iag o n al
p rin cip al
X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X
X X X X X
X X X X X
X X X
0 X X X X
0 0 X X X
X X X X
0
X X X 0 0
X X
O rg an izació n co n v en cio n al d e la
m atriz d e rig id ez
0
X 0 0
O rg an izació n en b an d a p ara O rg an izació n en b an d a p ara
m atrices n o sim étricas
m atrices sim étricas
Materiales y Métodos
Método Frontal
Este método trata de evitar el ensamblaje de la matriz
de rigidez global. Consiste básicamente en ensamblar la
menor cantidad de elementos y eliminar las variables al
mismo tiempo, formando matrices de menor tamaño.
Esto es, inmediatamente después de que los coeficientes
de una ecuación se ensamblan a partir de la contribución
de todos los elementos relevantes, o que inciden sobre este
coeficiente, la variable de este puede eliminarse. Como
resultado, la matriz completa de toda la estructura jamás
se forma.
Materiales y Métodos
Método Frontal
A diferencia del método banda, la adecuada
enumeración de los elemento es de gran importancia, no
así la enumeración global de los nodos ya que estos son
tratados por su numeración local. Como primer medida
esta característica del Front lo hace mas versátil ya que
como la cantidad de elementos siempre es menor que la
cantidad de nodos, independientemente de la forma de la
estructura, siempre se tendrá mayos facilidad para
enumerar los elementos que los nodos.
Discusión de Resultados
Para analizar las cualidades de los distintos métodos
propuestos se planteó la resolución mediante el empleo del
FECCUND V 2.1 (software desarrollado por el Dto. De
Mecánica Aplicada de la Facultad de Ingeniería de la
UNNE) de dos ejemplos prácticos, representando los
resultados mediante el pre-post procesador GID 7.2.
Para dicho análisis se tomaró como variables el tiempo
de ejecución de cada método y su exactitud, la cual será
determinada por comparación con los resultados
obtenidos por el software denominado CALSEF provisto
por el Centro Internacional de Métodos Numéricos en
Ingeniería (CIMNE).
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Mallado con elementos
hexaédricos
Desplazamiento nodal
según eje y
Discusión de Resultados
Barra sometida a esfuerzos de tracción
Tabla Nº 1
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
1,758
0,501953E-04
0,47619E-04
5,13
Banda
0,552
0,498874E-04
4,55
Frontal
0,441
0,479135E-04
0,61
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Mallado con elementos hexaédricos
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Desplazamiento nodal según eje y
Discusión de Resultados
Viga sometida a esfuerzos de flexión
Tabla Nº 2
Tiempo de
ejecución (seg)
Desplazamient Desplazamien
Error
o en z s/
to en z s/
(%)
FECCUND (m) CALSEF (m)
Gauss
31,547
-0,00615923
-0.00544118
11,66
Banda
2,025
-0,00560842
2,98
Frontal
1,752
-0,00552264
1,47
Conclusiones
El método Frontal es superior al Gauss y al Banda
en cuanto a su exactitud y a la velocidad de
ejecución.
El efecto de la propagación de errores por redondeo
incrementa notablemente la deficiencia del método
de Gauss frente a problemas de mayor
envergadura.
El método de Banda puede ser recomendable en el
caso de contar con un software de pre-post
procesamiento ya que los mismos generan las
mallas enumerando los nodos en la forma mas
adecuada disminuyendo de ancho de banda.