Modelos Matemáticos y Simulación
Ciclo 02-2017
Guía de Laboratorio 5. GiD como software para pre y post proceso.
Una ecuación diferencial en derivadas parciales (EDP), por su semejanza con las ODE, es
una ecuación donde una cierta función incógnita u viene definida por una relación entre sus
derivadas parciales con respecto a las variables independientes.
Si u = u(x,y,z), una ecuación diferencial en derivadas parciales sería
El orden de una EDP es el de la derivada parcial de mayor orden que aparece en la ecuación.
Se dice que una ecuación diferencial parcial es lineal, si es lineal en la función
desconocida y en todas sus derivadas, con coeficientes que dependen sólo de las variables
independientes.
Debido a su amplia aplicación en ingeniería, nuestro estudio de las EDP se concentrará en
las ecuaciones diferenciales lineales de segundo orden. Para dos variables independientes,
tales ecuaciones se pueden expresar de la forma general siguiente, donde A, B y C son
funciones de x y y, y D es una función de x, y, u, ∂u/∂x y ∂u/∂y.
Dependiendo de los valores de los coeficientes de los términos de la segunda derivada
(A, B y C):
Estas ecuaciones se pueden clasificar en
B2-4AC Categoría Ejemplo
<0 Elíptica Ecuación de Laplace (estado estacionario con dos dimensiones
espaciales)
=0 Parabólica Ecuación de conducción del calor (variable de tiempo y una
dimensión espacial)
>0 Hiperbólica Ecuación de onda (variable de tiempo y una dimensión espacial)
Modelos Matemáticos y Simulación
Ciclo 02-2017
Cada una de las categorías de ecuaciones diferenciales parciales corresponde a una clase
específica de problemas en ingeniería.
Comúnmente, las ecuaciones elípticas se utilizan para caracterizar sistemas en estado
estacionario. Como en la ecuación de Laplace, esto se indica por la ausencia de una derivada
con respecto al tiempo. Así, estas ecuaciones se emplean para determinar la distribución en
estado estacionario de una incógnita en dos dimensiones espaciales
Las ecuaciones parabólicas determinan cómo una incógnita varía tanto en el espacio como
en el tiempo, lo cual se manifiesta por la presencia de las derivadas espacial y temporal, como
la ecuación de conducción de calor. Tales casos se conocen como problemas de propagación,
puesto que la solución se “propaga”, o cambia, con el tiempo.
La categoría hiperbólica, también tiene que ver con problemas de propagación. Sin embargo,
una importante diferencia manifestada por la ecuación de onda, es que la incógnita se
caracteriza por una segunda derivada con respecto al tiempo. En consecuencia, la solución
oscila
La técnica general para hallar soluciones numéricas de sistemas de ecuaciones diferenciales
e integrales es el Método de Elementos Finitos (FEM). Esta técnica divide el dominio de la
solución en regiones con formas sencillas o “elementos”.
En esta técnica el dominio de la solución se divide en una malla con puntos discretos o
nodos(discretización). Entonces, se aplica la EDP en cada nodo, donde las derivadas parciales
se reemplazan por diferencias finitas divididas. Aunque tal aproximación por puntos es
conceptualmente fácil de entender, tiene varias desventajas. En particular, es difícil de aplicar
a sistemas con una geometría irregular, con condiciones de frontera no usuales o de
composición heterogénea. El método del elemento finito ofrece una alternativa que es más
adecuada para tales sistemas. A diferencia de las técnicas por diferencias finitas, la técnica
del elemento finito divide el dominio de la solución en regiones con formas sencillas o
“elementos”. Se puede desarrollar una solución aproximada de la EDP para cada uno de estos
elementos. La solución total se genera uniendo, o “ensamblando”, las soluciones
individuales, teniendo cuidado de asegurar la continuidad de las fronteras entre los
elementos. De este modo, la EDP se satisface por secciones.
Para resolución de ecuaciones diferencias parciales se empleará la ayuda de un software para
resolver directamente: GiD.
GiD
Es un programa gráfico diseñado para simulación numérica de pre a post procesamiento:
modelado geométrico, definición y preparación de los datos de análisis, así como a la
visualización de resultados.
Modelos Matemáticos y Simulación
Ciclo 02-2017
GiD fue desarrollado por International Center for Numerical Methods in Engineering
(CIMNE), es un centro de investigación con sede en Barcelona, España. CIMNE fue creado
en 1987, vinculado a la prestigiosa Universidad Politécnica de Cataluña. CIMNE se
especializa en el desarrollo aplicaciones de métodos numéricos y software para la búsqueda
de soluciones a una variedad de problemas en ingeniería y ciencias aplicadas. Toda
información correspondiente a GiD pueden acceder a: [Link]
GiD como pre-procesador es un sistema CAD que cuenta ampliamente con superficies
NURBS, sirve para:
▪ Diseñar geometrías/dominios de problemas de elemento finito en 1D, 2D y 3D.
▪ Asignar materiales.
▪ Definir condiciones iniciales y de frontera.
▪ Generar mallas para el dominio con diferentes tipos de elementos: líneas, triángulos,
cuadriláteros, tetraedros, hexaedros y prismas.
▪ Crear interfaces de usuario para declarar materiales, condiciones, parámetros, etc.
▪ Generar archivos de datos con la información anterior.
▪ Puede llamar programas externos (“solvers”).
GiD como post-procesador sirve para:
▪ Puede leer datos de resultados generados por los “solvers”.
▪ Visualizar resultados: desplazamientos, campos escalares, campos vectoriales, flujos.
Generar animaciones.
▪ Efectos de visualización: cortes, texturas en superficies, etc.
Entorno Grafico
Modelos Matemáticos y Simulación
Ciclo 02-2017
Donde:
1. Top menú: ofrece varios tipos de comandos. Las opciones que se muestran dependen
en el ambiente de usuario del análisis que se está ejecutando ya sea: pre-proceso o
post-proceso. Para cambiar la configuración actual de análisis de GiD se realiza el
cambio en:
Top menú: Archivos → Postproceso
o
Top menú: Archivos → Preproceso
2. Barra de Herramientas: entre las que integran el entorno de GiD al ejecutarse son:
Operaciones de Geometría y Vista, y barra de herramientas estándar. Para configurar
su posición en la ventana o mostrarlas o no mostrarlas, se realiza el cambio en:
Top menú: Utilidades → Herramientas → Barra de Herramientas
3. Línea de Comandos: es la opción al usuario que permite la ejecución de todos los
comandos de GiD, sin tener que acceder a través de los menús. No se hace distinción entre
el uso de letras mayúsculas y minúsculas.
4. Botones Derechos: es un menú rápido que provee al usuario de botones de opciones
de acuerdo a la herramienta que se usando. La función de la tecla Escape es de importancia
en el uso del programa que permite al usuario en casi todas las acciones realizadas declarar
que la acción ha sido terminada y salir de ella.
Preproceso
A través de un ejemplo se introducirá las herramientas básicas para la creación de entidades
geométricas y la generación de mallas. Se diseña de forma jerárquica: punto, línea, superficie,
volúmenes. En GiD no tiene definida sus unidades.
Modelos Matemáticos y Simulación
Ciclo 02-2017
Para crear un punto:
■ Desde Top menu: Geometría ^ Crear ^ Punto.
■ Se introduce por medio del mouse en el área de trabajo cada punto; O se
introduce en la línea Comandos las coordenadas del punto con formato: “x,y” o
“x,y,z” y [Enter]
■ [Esc] para terminar.
■ Es de color negro.
Para crear una linea:
■ Desde Top menú: Geometría ^ Crear ^ Línea recta; O hacer clic en el
botón de Crear Línea ubicado en la barra de herramientas de Geometría.
■ Se introduce por medio del mouse en el área de trabajo punto a punto; O se
introducen los puntos que definen la línea a través de la Línea Comandos
digitando las coordenadas del punto con formato: “x,y” o “x,y,z” y [Enter]
■ [Esc] para terminar.
■ Es de color azul.
Para crear una superficie:
■ Desde Top menu: Geometría ^ Crear ^ Superficie NURBS; O hacer clic
en el boton de Crear Superficie NURBS ubicado en la barra de herramientas
de Geometría.
■ Se seleccionan las lineas que encierran una superficie.
■ [Esc] para terminar.
■ Es de color magneta
Para crear un volumen:
■ Desde Top menu: Geometría ^ Crear ^ Volumen; O hacer clic en el
boton de Crear Volumen ubicado en la barra de herramientas de Geometría.
■ Se seleccionan las superficies que encierran un volumen.
■ [Esc] para terminar.
■ Es de color celeste.
Para crear un mallado:
■ Desde Top menu: Malla ^ Generar Malla.
■ Se ingresa el tamaño de los elementos a generar
■ En la etapa de modelado de la malla la version gratuita de GiD permite 1010 nodos.
■ [Ver malla] para terminar.
■ El boton de Conmutar vista geometria-malla permite ocultar o visualizar la malla.
■ Es de color gris.
Modelos Matemáticos y Simulación
Ciclo 02-2017
Ejemplo 1:
Disena la geometria en 2D: un cuadrado. Y generar la malla para el dominio con
diferentes tipos de elementos: triangulos, cuadrilateros y circulos.
Opcion 1 para crear el cuadrado:
■ Clic en el boton de Crear Linea.
■ Introducir coordenadas primera línea en la Linea Comandos: “0,0,0” ^
[Enter] y “10,0,0” ^ [Enter] ^ [Esc]
■ Clic en el boton de Crear Linea.
■ Introducir coordenadas segunda línea en la Línea Comandos: “0,10,0” ^
[Enter] y “10,10,0” ^ [Enter] ^ [Esc]
■ Clic en el boton de Zum Todo.
■ Para unir los puntos para crear el cuadrado:
Clic en boton de Crear Linea ^ boton derecho del mouse ^
Contextual ^ Juntar
Opcion 2 para crear el cuadrado:
■ Clic en el boton de Crear Linea ubicado en la barra de herramientas de Geometría.
■ Para crear el cuadrado introducir uno a uno en la Linea Comandos las
coordenadas: “0,0,0” ^ [Enter] ^ “10,0,0” ^ [Enter] “10,10,0” ^ [Enter] ^
“0,10,0” ^ “0,0,0” ^ [Enter] ^ Juntar ^ [Esc]
Para generar na superficie y la malla:
■ Clic en boton derecho del mouse ^ Etiquetas ^ Todo
■ Clic en boton de Crear Superficie
■ Se seleccionan una a una las lineas o las cuatro lineas a la vez ^ [Esc]
■ Clic en Malla ^ Generar malla ^ tamaño del elemento: 1 ^ Aceptar ^ Ver malla
Para generar una malla más uniforme:
■ Clic en Malla ^ Eliminar malla ^ Aceptar
■ Clic en Malla ^ Estructurada ^ Superficies ^ Asignar número de divisiones
■ Seleccionar superficie ^ [Esc]
■ Numero de divisiones: 10 ^ Asignar ^ seleccionar las lmeas ^ [Esc] ^ Cerrar
■ Clic en Malla ^ Generar malla ^ tamano del elemento: 1 ^ Aceptar ^ Ver malla
■ Clic en Malla ^ Tipo de elemento ^ Triangulo/Cuadrilatero/Circulo
■ Seleccionar superficie y^ [Esc]
■ Clic en Malla ^ Generar malla ^ tamano del elemento: 1 ^ Aceptar ^ Ver malla Para
concentrar en el mallado:
■ Clic en Malla ^ Estructurada ^ Lineas ^ Concentrar elementos
■ Seleccionar lineas ^ [Esc]
■ Peso inicio: 0.5
■ Aceptar ^ [Esc]
■ Clic en Malla ^ Generar malla ^ tamano del elemento: 1 ^ Aceptar ^ Ver malla
Guarde el proyecto con el nombre [Link].
■ Clic en Archivo -> Guardar como -> Nombre del archivo: Cuadrado -> Salvar
Modelos Matemáticos y Simulación
Ciclo 02-2017
Ejemplo 2:
Disena la geometria en 3D: un cubo. Y generar la malla para el dominio con
diferentes tipos de elementos: tetraedros y hexaedros.
■ Superficie Ejemplo 1
■ Clic en Malla ^ Eliminar malla ^ Aceptar
■ Clic en Malla ^ Eliminar datos mallado ^ Aceptar
■ Clic Utilidades -> Copiar
o Tipo de Entidades:
Superficies o Transformation:
Traslacion o Primer punto:
“0,0,0” o Segundo punto:
“0,0,10” o Extruir: Superficies
o Seleccionar: superficie del
ejemplo 1. o Terminar o
Cancelar
■ Clic en Vista ^ Rotar ^ Isometrico
■ Clic en Crear Volumen
■ Seleccionar todas las superficies ^ [Esc]
■ Clic en Malla ^ Generar malla ^ tamano del elemento: 2 ^ Aceptar ^ Ver
malla Para generar una malla más uniforme:
■ Clic en Malla ^ Estructurada ^ Volumenes ^ Asignar numero de divisiones
■ Seleccionar volumen ^ [Esc]
■ Numero de divisiones: 7 ^ Asignar ^ seleccionar las lineas ^ [Esc] ^ Cerrar
■ Clic en Malla ^ Generar malla ^ tamaño del elemento: 2 ^ Aceptar ^ Ver malla
■ Clic en Malla ^ Tipo de elemento ^ Tetraedros/ Hexaedros
■ Seleccionar volumen ^ [Esc]
■ Clic en Malla ^ Generar malla ^ tamano del elemento: 2 ^ Aceptar ^ Ver malla
Guarde el proyecto con el nombre [Link].
■ Clic en Archivo -> Guardar como -> Nombre del archivo: Cubo -> Salvar
Calculo
Un Tipo de Problema es un conjunto de archivos configurados por un desarrollador
para que el programa pueda preparar los datos a analizar. Permiten capturar datos
para problemas específicos de elemento finito para después exportarlos en archivo
de texto.
Pueden definir:
• Condiciones iniciales.
• Condiciones de frontera.
• Parametros para los solvers. o Materiales o Unidades de medida
Modelos Matemáticos y Simulación
Ciclo 02-2017
Tipo de Problema se define mediante archivos de texto.
La configuracion de GiD se realiza a través de archivos con formato de texto.
Se requieren los siguientes archivos:
• .prb: configuración del parámetro general (no asociado a entidades). o .mat:
configuración de materiales y sus propiedades. o .cnd: configuración de las
condiciones impuestas en el cálculo. o .bas: (archivo de plantilla) el archivo para
configurar el formato del archivo de intercambio que media entre GiD.
• .bat: el archivo que se puede ejecutar llamado desde GiD. Este archivo inicia el
módulo de cálculo.
Ejemplo de Problem type:
Centromasa3D.
Realiza el cálculo del centro de masa de una figura compleja en 3D. GiD en el pre-proceso
realiza una discretización del objeto en estudio y genera una malla de elementos. A cada uno
de los elementos se les asigna un material y algunas condiciones. Esta información de pre
procesamiento en GiD permite al módulo de cálculo generar resultados.
Módulos Mat-fem
Mat-fem es un programa informático diseñado para introducir a los usuarios en el uso
práctico del Método de Elementos Finitos (FEM) para el análisis de estructuras y problemas
de campo gobernados por las ecuaciones de Poisson y Laplace (transferencia de calor,
acústica, infiltración, electromagnetismo, etc.). Toda información correspondiente puede
acceder a: [Link]
Mat-fem está escrito en Matlab y tiene una interfaz gráfica de usuario amigable con el GiD.
Entre las opciones de módulos de Métodos de Elementos Finitos desarrollados para usar en
GiD se encuentran:
o 2D&3D: Análisis estructural para deformación plana en 2D y 3D. o Beams: Análisis
de vigas. o Plates: Análisis de placas.
o Shells: Análisis de estructuras como puentes, cubiertas, tanques, fuselaje de
aviones, carrocerías de vehículos, etc.
o Axi-Shells: Análisis de estructuras de cascaras asimétricas.
o Heat: Análisis mediante la implementación de la ecuación de Poisson tal como el
mecanismo de conducción de calor a través de un cuerpo en un medio permeable. o
Modelos Matemáticos y Simulación
Ciclo 02-2017
Sound: Análisis mediante la implementación de la ecuación bidimensional de
Helmholtz.
Ejemplo 3:
Calcular el centro de masa del Cubo del Ejemplo 2 por medio del Tipo de problema llamado
[Link] para resolver.
■ Cargue el archivo del proyecto realizado en el Ejemplo 2.
■ Clic en Malla ^ Eliminar malla ^ Aceptar
■ Clic en Malla ^ Eliminar datos mallado ^ Aceptar
■ Clic en Datos ^ Tipo de Problema ^ Cargar
Saldrá una advertencia: “Toda la información adicional se perderá”, Clic en Aceptar.
Buscar en el directorio la carpeta [Link] que contiene la definición del
Tipo de problema y abrir.
■ Clic en Datos ^ Materiales
Estos son los datos generales del tipo de problema. Para asignar un determinado
material y asociar a los elementos, se abre la ventana que posee la información de las
densidades de distintos materiales. Seleccione uno.
■ Clic en Asignar ^ Volúmenes ^ seleccionar Volumen ^-Terminar ^ Cerrar
■ Clic en Malla ^ Generar malla ^ tamaño del elemento: 2 ^ Aceptar ^ Ver malla
■ Clic en Calcular ^ Calcular
Se procede al cálculo del centro de masa del Cubo asi como la distancia de cada uno
de los nodos generados por la malla al centro de masa. Automáticamente GiD correrá
el programa hecho en MatLab para calcular el centro de masa.
■ Saldra la ventana: El proceso ha finalizado ^ Aceptar.
■ Los resultados se han guardado para visualizarse el postproceso.
Iniciación a GiD – Postproceso
El solver crea los archivos: .[Link] y .[Link] para visualizar los resultados en GiD.
Para cargar el modelo con los resultados de la información de la simulación se debe cambiar
la configuración actual de análisis de GiD:
■ Clic en Archivos ^ Postproceso
■ Clic Archivos^- Abrir ^ se seleccionan el archivo .[Link] o .[Link]
Los resultados están contenidos en sets. Los sets son las capas que contienen volúmenes,
superficies o cortes de los resultados.
Con la ventana de Estilos de visualización se pueden especificar diversas opciones para los
volúmenes, superficies y cortes de la figura en análisis. Entre estas opciones se pueden activar
y desactivar volúmenes, superficies y cortes; las propiedades de sus colores pueden ser
cambiadas, y su transparencia también.
Desde esta ventana los sets de volúmenes, superficies o cortes pueden ser eliminados o
modificar sus nombres. Para acceder a esta ventana en el análisis de postproceso:
■ Clic en Ventana ^ Estilo de visualización
Modelos Matemáticos y Simulación
Ciclo 02-2017
Ahora se podrá visualizar todos los resultados del Ejemplo 3 abra los archivos
correspondientes y siga los siguientes pasos:
■ Clic Archivos^- Abrir ^ carpeta [Link] ^ se seleccionan el archivo .[Link] o
.[Link]
■ Clic en Ver resultados ^ Áreas coloreadas ^ Distancias CM
