1-FEAP Guion PDF
1-FEAP Guion PDF
1-FEAP Guion PDF
DEL
FEAP: APLICACION
METODO
DE LOS
ELEMENTOS FINITOS
Autores:
Enero de 2012
Departamento de Mecanica de Estructuras e Ingeniera Hidraulica
Universidad de Granada
Introducci
on
Tras estudiar la teora del Metodo de los Elementos Finitos, >sabra denirlo?
El Metodo de los Elementos Finitos es un metodo numerico que permite resolver de forma aproximada numerosos problemas de ingeniera. Cosiste en la discretizacion de un solido
en una serie de elementos interconectados, los elementos nitos. Dentro de cada elemento
se denen una serie de nodos que tienen asignados una serie de incognitas, deniminadas
grados de libertad. Las relaciones existentes entre los grados de libertad forman un sistema de ecuaciones. En el analisis de estructucturas mediante el Metodo de los Elementos
Finitos, las incognitas del sistema son los desplazamientos. Tras la obtencion del campo de
desplazamientos, se obtienen el campo de deformaciones y tensiones, al aplicar relaciones
cinematicas y ley de Hooke, respectivamente.
Una de las propiedades mas importantes del metodo es la convergencia. Al considerar
discretizaciones mas nas, la solucion numerica calculada converge rapidamente hacia la
solucion exacta del sistema de ecuaciones. Si bien, hay que tener en cuenta que los tiempos
de calculo aumentan con el numero de incognitas del problema por lo que habra que llegar
a un equilibrio entre ambos.
1.1.
Objetivos generales
Los principales objetivos que se pretenden alcanzar con la ejecucion de esta practica se
enumeran a continuacion:
1. Aprender a discretizar una estructura
2. Distinguir los tipos de estructuras
3. Analizar los tipos de materiales y su comportamiento
4. Resolver un problema mediante el Metodo de los Elementos Finitos
1.2.
Contenido
El contenido de esta gua se divide en dos partes principales. En la pirmera se hara una
aproximacion al calculo de estructuras de barras, incidendiendo en los contenidos del problema y enfatizaondo en los aspecto especcos de este tipo de estructural. La segunda parte
consiste en la aproximacion de calculo de elasticidad bidimensional, enfatizando con mayor
interes los nuevos aspectos a considerar en el calculo que dieran de los calculos anteriores.
Ambos apartados concluyen con unos ejemplos de autoevaluacion para el alumno.
Actividades
En este punto nos centraremos en el programa de elementos nitos que vamos a emplear
para la resolucion de problemas. FEAPpv (Finite Element Analisys Program personal version ) desarrollado por el profesor R. L. Taylor de la Universidad de Berkeley:
http://www.ce.berkeley.edu/projects/feap/feappv/ es un programa academico que
permite al alumno acercarse al empleo de diferentes tipos de elementos y modelos de resolucion de problemas mediante el metodo de los elementos nitos.
2.1.
Objetivos
1. Introducir los conceptos basicos para la formulacion del Metodo de los Elementos
Finitos para estructuras de barras, distinguiendo las estructuras de barras articuladas
y nudo rgidos.
2. Aplicar estos conceptos en la implementacion en FEAP.
3. Aproximacion al interfaz graco de FEAP.
4. Entender la organizacion del programa FEAP.
Definici
on del problema
Para el analisis de una estructura mediante el metodo de los elementos nitos, el primer
paso es proceder a la discretizacion de la region en elementos y nodos. Dentro de este proceso
hay que considerar tambien los aspectos siguientes: (1) el tipo de elementos a utilizar, (2)
donde localizar los nodos, (3) como aplicar las condiciones de contorno y cargas, (4) el
2
modelo del material y que propiedades son necesarias para su denicion y (5) aquellos
aspectos particulares que requiera el problema que estemos analizando.
Aproximaci
on a FEAP
FEAPpv es un programa de elementos nitos programado en codigo Fortran. Los datos
del problema se deenen mediante el empleo de comandos que se iran redactando en un
archivo de texto: Abra el archivo Iej0 con el bloc de notas y observe la estructura general
de un problema de elementos nitos desarrollado para su calculo con FEAPpv.
1. Cada lnea puede tener como maximo 255 caracteres en formato libre.
2. El tama~no maximo para un valor es de 15 caracteres.
1.234567890123*1.e-5
3. Cada conjunto de datos comienza con la denicion de un comando, del que el programa
reconoce los cuatro primeros caracteres y termina con una lnea en blanco.
4. Se pueden incluir comentarios mediante la inclusion del caracter !
5. Se pueden denir parametros como un caracter simple, dos caracteres o un caracter y
un numero. El comando para su denicion es PARAmeter.
PARAmeter
a = 3.567
nu = 0.3
e1 = 200.0e9
! fin de lnea
6. Se pueden emplear operaciones matematicas (^, , /, +, ) con comandos o directamente con valores numericos y es posible el uso de las principales funciones trigonometricas (abs, exp, int, log, sqrt, sin, cos, tan, asin, acos, atan). La denicion de
los angulos de las funciones trigonometricas se realiza en radianes, si bien se permite
la denicion en grados sexagesimales mediante el empleo de las funciones sind, cosd,
tand, atand, asind y acosd.
Cada funcion solo tiene un argumento entre parentesis, no es posible poner funciones
o mas paretesis en el argumento.
Comando de inicio
Para preparar el problema a resolver, se debe denir un archivo en formato texto sin
extensi
on1 . El nombre dado al archivo debe comenzar con la letra I, por ejemplo, vamos a
Abrimos el Bloc de Notas y comenzamos con la denicion del problema mediante dos
lneas:
1. La primera lnea debe ser el comando FEAPpv. Indistintamente se puede escribir en
mayusculas o minusculas. Las columnas desde la 5 a la 80 se reservan para denir el
ttulo del problema.
2. La lnea segunda contiene informacion sobre el tama~no del problema2 :
b)
c)
d)
e)
f)
NUMEL: N
umero total de elementos;
NUMMAT: N
umero total de materiales;
NDM: Dimension espacial del problema;
NDF: N
umero maximo de incognitas por nodo;
NEN: N
umero maximo de nodos por elemento.
FEAPpv **titulo**
NUMNP, NUMEL, NUMMAT, NDM, NDF, NEN
! fin de lnea
Definici
on de Nodos y Elementos
El comando ELEMent es empleado para denir las conexiones entre los nodos. Para
elementos con numero de nodos inferior a 13, la denicion del comando es la que sigue:
ELEMent
1, 1G, MA, (NDi, i=1,NEN)
2, 2G, MA, (NDi, i=1,NEN)
Las condiciones de contorno nodales pueden denirse mediante los comandos BOUndary,
FORCe, DISPlacement y ANGLe. Para estos comandos es necesario especicar el n
umero del
nodo. Aunque tambien es posible denir las condiciones de contorno a partir de coordenadas
o bordes (vea el Manual, pp. 32)
La estructura para estos comandos es la que sigue:
BOUNdary
nodo1, ngen1,
nodo2, ngen2,
! fin de lnea
FORCe
nodo1, ngen1,
nodo2, ngen2,
! fin de lnea
DISPlacement
nodo1, ngen1,
nodo2, ngen2,
! fin de lnea
ANGLe
nodo1, ngen1,
nodo2, ngen2,
! fin de lnea
(id(i,nodo1),i=1,ndf)
(id(i,nodo2),i=1,ndf)
(f(i,nodo1),i=1,ndf)
(f(i,nodo2),i=1,ndf)
(d(i,nodo1),i=1,ndf)
(d(i,nodo2),i=1,ndf)
angl(nodo1)
angl(nodo2)
El comando BOUNdary se usa para especicar las condiciones de apoyo. Para cada nodo
i, se especica el n
umero de nodo nodoi, el incremento de generacion al siguiente nodo ngeni,
si es necesario, en caso contrario se dene con 0 y el codigo para la restriccion del grado de
libertad desde la direccion 1 hasta la direccion ndf. El codigo de la restriccion puede ser 0
si se impide el desplazamiento y 1 si la fuerza aplicada es nula.
El comando FORCe permite especicar los valores de las fuerzas aplicadas en los nodos.
Para cada nodo i, se dene el numero de nodo nodoi, el incremento de generacion al siguiente
nodo ngeni, si es necesario, en caso contrario se dene con 0 y el valor de la fuerza para la
direccion 1 hasta la direccion ndf.
El comando DISPlacement indica el valor de los desplazamientos impuestos en nodos.
Para un nodo i, se dene el numero de nodo nodoi, el incremento de generacion al siguiente
nodo ngeni, si es necesario, en caso contrario se dene con 0 y el valor del desplazamiento
en la direccion 1 hasta la direccion ndf.
El comando ANGLe se emplea para especicar el angulo en grados que gira un nodo.
Para un nodo i, se dene el numero de nodo nodoi, el incremento de generacion al siguiente
nodo ngeni, si es necesario, en caso contrario se dene con 0 y el angulo que gira el nodo
angl(nodoi).
Tipo de material y sus propiedades
FEAPpv contiene una librera con elementos estandar y modelos de material que pueden
ser empleados en multimple problemas de mecanica de estructuras y solidos. El tipo de
elemento a emplear en el analisis se especica dentro del comando MATErial. Los tipo de
elementos incluidos en esta version del programa son los recogidos en la pagina 40 del manual.
As mismo, los modelos de material que pueden ser utilizados se pueden ver descritos con
detalle en la pagina 47 del manual.
MATErial, ma, etiqueta
tipo, iel, id,(idf(i),i=1,ndf)
! fin de lnea
La estructura que estamos analizando esta compuesta por barras que trabajan a axil
y todas tienen las mismas propiedades elasticas, por lo que unicamente sera necesaria la
denicion de un tipo de material.
MATErial, 1
TRUSs
Ya hemos denido la malla, condiciones de contorno y material del problema. Es, en este
punto, cuando se debe emplear el comando END, que termina la entrada de datos y le da el
control de ejecucion al programa.
Comandos de lenguaje del programa
Para proceder a solucionar el problema, el archivo de entrada hay que proceder a emplear
el comando BATCh. Los comandos necesarios para resolver un problema estatico lineal son
los siguientes:
BATCh ! pasa al modo computacional
TANG ! forma la matriz del sistema
FORM ! forma el vector de terminos independientes del sistema
SOLVe ! resuelve el sistema de ecuaciones
END ! finaliza el modo computacional
Los resultados pueden obtenerse de forma numerica o graca. Los datos numericos que
podemos escribir son los desplazamientos, tensiones y reacciones mediante los comandos
DISP,ALL, STRE,ALL y REAC,ALL, respectivamente. Los resultados se escribiran en pantalla
y, a la vez, en el archivo de salida que, normalmente, tendra el nombre del archivo de entrada
pero comenzando por O.
Para la creacion de gracos es necesario pasar al modulo de dibujo, para ello escribimos
el comando PLOT
7
A continuacion, complete la tabla la descripcion del comando ejecutado dentro del modulo PLOT.
MESH
NODE
ELEMents
LOAD
REACtions
BOUNdary
WIPE
DEFO,100,1
MESH
DISPlacements
STRess,1
>Como puede comprobar que el resultado es valido?
8
2.2. CONCLUSIONES
Autoevaluaci
on
Tiene que tener superados los siguientes aspectos:
1. Tiene la base del concepto de FEAP del empleo de comandos.
2. Sabe seguir el analisis tpico de un problema de elementos nitos con: (1) fase de
preproceso: donde especica la geometra, material y tipo de elemento al programa,
junto a condiciones de contorno y acciones, (2) fase de proceso mediante la resolucion
y (3) fase de postproceso: con el analisis de resultados.
3. Sabe trabajar con una estructura de nudos articulados. Sabe que propiedades denen
el material y cuantos grados de libertad tiene cada nodo.
2.2.
Conclusiones
En este punto es muy importante tener claro que se pretende con este Trabajo:
1. Tiene las bases para poder empezar a resolver una estructura mediante el Metodo de
los Elementos Finitos con cualquiera de los programas que podra encontrar a lo largo
de su vida profesional.
2. Sabe cuales son los pasos a seguir para resolver un problema mediante el metodo de
los elementos nitos.
b ) Analisis.
1) Tareas previas al calculo: numeracon global de grados de libertad y comprobaciones respecto al calculo (memoria suciente, etc.)
2) Calculo de la matriz de rigidez global K y vector de cargas f.
3) Resolucion del sistema Ku = f.
2.2. CONCLUSIONES
10
Bibliografa
Ramon Arguelles Alvarez. Fundamentos de Elasticidad y su Programacion por Elementos
Finitos. Bellisco, Madrid. 1992.
O.C. Zienkiewicz, R.L. Taylor. El Metodo de los Elementos Finitos. Vols 1 y 2. CIMNEMc Graw Hill, 1994.
E. O~nate. Calculo de Estructuras por el Metodo de los Elementos Finitos. CIMNE,
Barcelona. 1995
11