UNIVERSIDAD AUTÓNOMA DE COAHUILA

FACULTAD DE INGENIERÍA MECÁNICA Y ELÉCTRICA U. N.

NOTAS DE LA MATERIA DE DINÁMICA

PARA INGENIEROS MECÁNICOS ELECTRICISTAS

REALIZÓ:

Dr. en C. Sixtos Antonio Arreola Villa

e-mail:[email protected]

Monclova, Coahuila, México, Agosto 2022

 Índice general

1. Introducción al método de elementos finitos 13

1.1. Método numérico en la ingenierı́a mecánica . . . . . . . . . . . . . . . 15
1.1.1. Evolución del método de elemento finito . . . . . . . . . . . . 18
1.2. Conceptos básicos del MEF . . . . . . . . . . . . . . . . . . . . . . . . 20
1.2.1. Proceso de Discretización . . . . . . . . . . . . . . . . . . . . 20
1.2.2. Ventajas y limitaciones del MEF . . . . . . . . . . . . . . . . . 23
1.3. Procedimiento general del MEF . . . . . . . . . . . . . . . . . . . . . 25
1.4. Enfoque directo de análisis por MEF . . . . . . . . . . . . . . . . . . . 29
1.4.1. Elemento resorte lineal en 1D . . . . . . . . . . . . . . . . . . 29
1.4.2. Elemento resorte de torsión . . . . . . . . . . . . . . . . . . . 32
1.4.3. Elemento de conducción . . . . . . . . . . . . . . . . . . . . . 33
1.5. Enfoque energético de análisis por MEF . . . . . . . . . . . . . . . . . 36
1.5.1. Estructuras 1D de resortes lineales . . . . . . . . . . . . . . . . 37

2. Esfuerzos 40
2.1. Teorı́a básica y ecuaciones de elasticidad . . . . . . . . . . . . . . . . . 40
2.2. Modelado de armaduras por MEF . . . . . . . . . . . . . . . . . . . . 40
2.3. Modelado de barras por MEF en una y dos dimensiones . . . . . . . . . 40
2.4. Teorı́a básica para el análisis de vigas y marcos . . . . . . . . . . . . . 40
2.5. Modelado de vigas y marcos por MEF en una y dos dimensiones . . . . 40
2.6. Modelado en dos dimensiones de placas . . . . . . . . . . . . . . . . . 40
2.7. Simulación de elasticidad en tres dimensiones de placas utilizando soft-
ware basado en MEF . . . . . . . . . . . . . . . . . . . . . . . . . . . 40

3. Análisis térmico 41
3.1. Teorı́a básica y ecuaciones de transferencia de calor . . . . . . . . . . 41
3.2. Modelado de campos térmico por MEF . . . . . . . . . . . . . . . . . 41
3.3. Modelado de esfuerzos elásticos de origen térmico por MEF . . . . . . 41
3.4. Simulación de transferencia de calor y esfuerzos térmicos utilizando
software basado en MEF . . . . . . . . . . . . . . . . . . . . . . . . . 41

4. Análisis del campo fluidodinámico por MVF 42

4.1. Conceptos básicos del flujo de fluidos . . . . . . . . . . . . . . . . . . 42

II
ÍNDICE GENERAL III

4.2. Ecuaciones de Navier-Stokes . . . . . . . . . . . . . . . . . . . . . . . 42

4.3. Procedimiento general de análisis por MVF . . . . . . . . . . . . . . . 42
4.4. Modelado de problemas de flujo de fluidos por MVF . . . . . . . . . . 42
4.5. Simulación de flujo de fluidos utilizando software basado en MVF . . . 42
 Índice de figuras

1. Pensamiento de Isaac Newton, padre de la Mecánica Clásica. . . . . . . 12

1.1. Esquema conceptual sobre los métodos utilizados para la resolución de

ecuaciones diferenciales. . . . . . . . . . . . . . . . . . . . . . . . . . 16
1.2. Divisibilidad: a) espacio infinito y b) triángulo divisible infinitamente. . 20
1.3. Conceptos de continuidad: a) consecutivos, b)contiguos y c) continuo . 21
1.4. Concepto de convergencia . . . . . . . . . . . . . . . . . . . . . . . . 21
1.5. Convergencia para aproximar el área de un cı́rculo . . . . . . . . . . . . 22
1.6. Concepto de lı́mite . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
1.7. Discretización de un cı́rculo mediante un polı́gono inscrito . . . . . . . 27
1.8. Concepto de lı́mite . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
1.9. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
1.10. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
1.11. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
1.12. Elemento resorte de torsión, a) esquema del elemento, momento y ro-
taciones nodales, b)relación lineal momento-rotación . . . . . . . . . . 32
1.13. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
1.14. Concepto de lı́mite . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
1.15. Concepto de lı́mite . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
1.16. Concepto de lı́mite . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
1.17. Problema conceptual de transferencia de calor . . . . . . . . . . . . . 36
1.18. Conjunto estructural formado por dos resortes lineales. . . . . . . . . . 37
1.19. a)resorte lineal sujeto a cargas, b) barra empotrada en A sujeta a carga
en C . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39

IV
 Índice de cuadros

1.1. Valores numéricos de la solución exacta y aproximada. . . . . . . . . . 28

1
ÍNDICE DE CUADROS 2

Definición de Método de los elementos Finitos (MEF) se le llama al procedimiento

para encontrar una solución numérica de las ecuaciones que gobiernan los problemas
que se encuentran en la naturaleza. Usualmente el comportamiento de la naturaleza
puede describirse como un conjunto de ecuaciones expresadas en forma diferencial o
integral.
Por esta razón el MEF es entender temas matemáticos como técnica numérica para re-
solver ecuaciones diferenciales parciales o ecuaciones integrales. Generalmente el MEF
permite obtener la evolución en espacio tiempo de una o mas variables que representan
el comportamiento de un sistema fı́sico.
Cuando se hable del análisis de estructuras el MEF es un método potente para calcular
el desplazamiento, la deformación, y esfuerzo en estructuras sometidas a cargas. Siendo
precisamente uno de los objetivos del curso. Por otro lado, como objetivo quedará la
revisión del método de volumen finito (MVF), como técnica utilizada en la ingenierı́a
de fluidos.

Evaluación El presente curso será evaluado de acuerdo con los siguientes criterios:
Asistencia 10 %
Tres exámenes parciales con 25 % cada uno, dando un total de 75 %
Tareas y participaciones en clase 15 %

Reglamento El total de los porcentajes anteriores da el 100 % que valdrá la evalua-

ción de Método de los Elementos Finitos. Como responsabilidad principal del alumno
en esta materia, será asistir a todas las clases virtuales o presenciales a la hora que
se programa en el horario durante el transcurso del semestre, debido que el curso es-
colarizado, el profesor tomará lista de de asistentes durante la clase en el aula, y es
responsabilidad del alumno estar presente. En ese sentido, el total de clases impartidas
en el periodo de evaluación corresponderá al 10 % del valor de ponderación en el pro-
medio final. En caso de presentar faltas a clase durante los pases de lista ese porcentaje
será proporcional a las asistencias. Por otro lado, el alumno deberá entregar los traba-
jos completos en las fechas acordadas en hojas separas y realizadas a mano, donde el
alumno entregará evidencia de su trabajo o tarea por algún medio electrónico o fı́sico.
Un requisito primordial será que el alumno deberá asignar un cuadernillo exclusivo pa-
ra esta materia, debido que el profesor revisará que las notas de clase estén completas,
por lo tanto, le pedirá al alumno que lleve el cuaderno en los dı́as cercanos de evalua-
ción a la Facultad, por otro lado, el maestro estará atento para atender a los alumnos
por separado en caso de ser necesario y resolver dudas.
Los exámenes serán presenciales, por lo que se puede asignar alguna hora dentro o
fuera de clases para citar para examen en la Facultad, ya sea en un salón o laboratorio.
En caso de actividades escolares o extraescolares donde se involucre al grupo, que su-
cedan durante la hora de clase, será requisito solicitar el permiso al docente, de forma
verbal o escrita con anticipación. En ese sentido, una falta grupal, provocará una san-
ción el dı́a inmediato del curso, donde se evaluará el tema que no fue posible abordar
por esta situación.
Finalmente, el alumno se portará de forma respetuosa fuera y dentro del salón, con el
ÍNDICE DE CUADROS 3

docente y sus compañeros de clase. Por otro lado, no se permite generar polémica o
distracción en clase con temas ajenos a la clase. Cualquier tema particular que deba ser
tratado será entre el profesor y el alumno que presente dicha situación.
La responsabilidad el docente será completar e impartir el total del contenido de la ma-
teria, en este ciclo escolar, donde se estarán abordando y buscando diversas estrategias
para verificar cual funciona de la mejor manera, al alumno se le pide paciencia ya que
es la primera vez que se ofrece esta materia, siendo los temas relativamente nuevos.
Finalmente, se pide que el alumno sea responsable también de traer su equipo de
computo o conseguir un laboratorio cada vez que sea necesario.
 UNIVERSIDAD AUTÓNOMA DE COAHUILA
DIRECCIÓN DE ASUNTOS ACADÉMICOS
Subdirección de Superación Académica

Facultad de Ingeniería Mecánica y Eléctrica

Unidad Norte
Ingeniero Mecánico Electricista

Método de Los Elementos Finitos

Ingeniería Aplicada y Diseño de Ingeniería

1.- Datos Generales de la Asignatura

Nombre de la Asignatura: Método de los Elementos Finitos

Clave asignatura:

Nivel: Licenciatura

Carrera: Ingeniero Mecánico Electricista

Frecuencia (h/semana) Teóricas: 2 Prácticas: 2 Total: 4

Requisitos Obligatorios:

Laboratorio: Diseño Avanzado y Simulación matemática

Total horas en el período escolar: 64

Créditos SATCA: 4

Docentes Participantes: Dr. Sixtos Antonio Arreola Villa

Dr. Ismael Calderón Ramos

MC Mario Alberto Barrera Moreno

Fecha de Elaboración: Mayo 2019

2.-Presentación

La asignatura Método de los elementos Finitos (MEF) se presenta, como una herramienta para el diseño
de elementos mecánicos. Abordando el análisis estructural, térmico y dinámica de fluidos. Esta materia se
recomienda en los cursos finales, siendo necesario que el estudiante tenga acreditadas las materias de
Mecánica de Materiales, Transferencia de Calor y Mecánica de Fluidos. También, se presenta de forma
introductoria el Método de Volumen Finito (MVF), para el estudio de flujo de fluidos.

Departamento de Desarrollo Curricular

 UNIVERSIDAD AUTÓNOMA DE COAHUILA
DIRECCIÓN DE ASUNTOS ACADÉMICOS
Subdirección de Superación Académica

Intención didáctica.

El temario se encuentra organizado en cuatro unidades. En la primera unidad se aborda, de forma

introductoria a los métodos numéricos utilizándolos en la ingeniería. Posteriormente se presentan los
conceptos y procedimientos de análisis por MEF, mostrando el enfoque directo y energético.

En la segunda unidad se estudian diversos tipos de elementos mecánicos en una, dos y tres dimensiones,
enfocando el análisis a esfuerzos y deformaciones elásticas producidas por fuerzas externas, utilizando el
MEF. Para solucionar las matrices de rigidez de la solución de problemas de pizarrón se utilizará la
programación en software. Finalmente, en la solución de sistemas mecánicos en dos y tres dimensiones se
utilizará un software basado en MEF.

La tercera unidad se plantea de la misma manera que la unidad anterior, pero con el enfoque a estudiar
los esfuerzos y deformaciones elásticas producidas por gradientes térmicos.

En la cuarta unidad se tratan los conceptos básicos del MVF, para el modelado de fluidodinámica.
Primeramente, se abordan los problemas de difusión pura, después los problemas difusivos-convectivos y
finalmente se realizan modelos numéricos en tres dimensiones utilizando software basado en MVF.

Propósito del Curso: El alumno aprenderá a utilizar métodos y herramientas numéricas de cálculo y
modelado para sistemas mecánicos.

3.-Aportación de la asignatura al Perfil del Egresado:

Conocimientos: Conocimiento y desarrollo de Métodos numéricos de modelado en partes y equipos

mecánicos, sometidos a cargas.

Habilidades: Manejo del modelado y simulación numérica como herramienta en la construcción de

equipos mecánicos, con y sin manejo de software especializado.

Actitudes: Carácter para decidir sobre posibles modificaciones en el diseño de equipos mecánicos,
utilizando la teoría elástica.
Materias de prerrequisitos Obligatorio Sugerido

Departamento de Desarrollo Curricular

 UNIVERSIDAD AUTÓNOMA DE COAHUILA
DIRECCIÓN DE ASUNTOS ACADÉMICOS
Subdirección de Superación Académica

4.-Objetivo General del Curso:

Presentar el MEF como una herramienta para el diseño de sistemas mecánicos utilizando análisis
estructural y térmico. Además de introducir el MVF para el modelado de campos fluidodinámicos,
dejando bases para la utilización de software de simulación.

5.- Objetivos específicos del curso

5.1 Objetivos Conceptuales 5.2 Objetivos Procedimentales 5.3 Objetivos Actitudinales y

(aprender a hacer) Valorativos (aprender a convivir y a
(aprender a aprender)
ser)
El estudiante aplicará los
conceptos y criterios de diseño, El estudiante fortalecerá las
El estudiante aprenderá a utilizando MEF y el MVF, para el actitudes hacia los demás (saber
conocer y comprender los cálculo de elementos de convivir) y sus valores (saber ser),
conceptos del MEF, para máquinas y sistemas mecánicos mediante la aplicación y ejecución
utilizarlos junto con los reales, generando de actividades propias de la materia
conocimientos adquiridos en conocimientos y habilidades adquirirá habilidades y
cursos anteriores, en las para aumentar su capacidad y conocimientos que fortalecerán su
materias de mecánica de visión en la toma de decisiones. criterio como futuro ingeniero y su
materiales, transferencia de confianza en la toma de decisiones.
calor y mecánica de fluido, en la Estas actividades promoverán el
solución de sistemas mecánico trabajo colaborativo y la resolución
de forma numérica. Por otro conjunta de problemas, buscando
lado, el estudiante conocerá de que el estudiante tenga la
forma introductoria el MVF autoestima, la capacidad de juicio
para el análisis de flujo de moral y la autonomía en la toma de
fluidos. decisiones al trabajar en equipo.

6.- Contenidos Temáticos

7.- Estrategias de Aprendizaje-Enseñanza

6.1 Nombre 6.2 6.3 Actividades Actividades Medios

de Unidades Objetivos Descripción realizadas por el realizadas por didácticos y
Temáticas específicos detallada del alumno el profesor recursos
de las contenido de utilizados

Departamento de Desarrollo Curricular

 UNIVERSIDAD AUTÓNOMA DE COAHUILA
DIRECCIÓN DE ASUNTOS ACADÉMICOS
Subdirección de Superación Académica

Unidades las Unidades

Temáticas Temáticas

Unidad 1 - Estudiar los 1.1 Métodos Resolver Exponer ante el Pintarrón,

conceptos numéricos en problemas grupo de forma marcadores y
fundamental la ingeniería conceptuales de general los cañón.
“INTRODUCCI es de la mecánica. elementos distintos
ÓN AL formulación mecánicos métodos
1.2 Conceptos
MÉTODO DE directa en el utilizando la numéricos
básicos del
ELEMENTOS cálculo de formulación utilizados en la
MEF.
FINITOS” elementos básica del MEF. ingeniería.
mecánicos 1.3
Exponer los
elásticos. Procedimiento
antecedentes
general del
- Estudiar los de MEF.
MEF.
conceptos
primarios del 1.4 Enfoque
MEF. directo de
análisis por
MEF.

1.5 Enfoque
energético de
análisis por
MEF.

Unidad 2 -Desarrollar 2.1. Teoría -Resolución de -Exponer ante Pintarrón,

el MEF para básica y elementos grupo la marcadores y
diseñar ecuaciones de mecánicos formulación del cañón.
“ESFUERZOS” elementos elasticidad. enfocados alMEF, en la
mecánicos, estudio de resolución de
2.2 Modelado
enfocados al esfuerzos yelementos
de armaduras
estudio de deformaciones mecánicos
por MEF.
esfuerzos y elásticas, enfocados al
deformacion 2.3 Modelado utilizando el MEF. análisis de
es elásticas de barras por esfuerzos y
-Resolución de
en una, dos y MEF en una y deformaciones
matrices de
elásticas, en
rigidez
Departamento de Desarrollo Curricular
 UNIVERSIDAD AUTÓNOMA DE COAHUILA
DIRECCIÓN DE ASUNTOS ACADÉMICOS
Subdirección de Superación Académica

tres dos ensambladas una, dos y tres

dimensiones. dimensiones. utilizando dimensiones, de
software. forma
2.4 Teoría
numérica,
básica para el -Resolución de
mostrando el
análisis de sistemas
procedimiento
vigas y mecánicos
para discretizar
marcos. enfocados al
el dominio de
estudio de
2.5 Modelado cálculo y el
esfuerzos y
de vigas y proceso de
deformaciones
marcos por ensamble-
elásticas,
MEF en una y solución de
utilizando
dos matrices.
software
dimensiones.
comercial basado -En el caso de
2.6 Modelado en MEF. tres
en dos dimensiones
dimensiones desarrollar
de placas. modelos de
simulación
2.7 Simulación numérica en
de elasticidad Software
en tres basado en MEF.
dimensiones
utilizando
Software
basado en
MEF.

Unidad 3 -Desarrollar 3.1 Teoría -Resolución de -Exponer ante Pintarrón,

el MEF para básica y elementos grupo la marcadores y
“ANÁLISIS
diseñar Ecuaciones de mecánicos formulación del cañón.
TÉRMICO”
elementos transferencia enfocados al MEF, en la
mecánicos, de calor. estudio de resolución de
enfocados al esfuerzos y sistemas
3.2 Modelado
estudio de deformaciones mecánicos
de campos
esfuerzos y elásticas de enfocados al
térmico por
deformacion origen térmico, análisis de
MEF.
es elásticas utilizando el MEF. esfuerzos y

Departamento de Desarrollo Curricular

 UNIVERSIDAD AUTÓNOMA DE COAHUILA
DIRECCIÓN DE ASUNTOS ACADÉMICOS
Subdirección de Superación Académica

de origen 3.3 Modelado -Resolución de deformaciones

térmico en de esfuerzos sistemas elásticas, en
una, dos y elásticos de mecánicos una, dos y tres
tres origen térmico enfocados al dimensiones, de
dimensiones. por MEF. estudio de forma
esfuerzos y numérica,
3.4 Simulación
deformaciones
de -En el caso de
elásticas de
transferencia tres
origen térmico,
de calor y dimensiones
utilizando
esfuerzos desarrollar
software
térmicos modelos de
comercial basado
utilizando simulación
en MEF.
Software numérica en
basado en Software
MEF. basado en MEF.

Unidad 4 Desarrollar 4.1 Conceptos -Modelado de -Exponer ante Pintarrón,

las bases de básicos del sistemas de grupo la marcadores y
“ANÁLISIS DE
la flujo de fluidos difusión pura, en formulación del cañón.
CAMPOS
metodología una y dos MVF para las
FLUIDODINÁ 4.2 Ecuaciones
numérica dimensiones por ecuaciones de
MICOS POR de Navier-
empleada en MVF. conservación-
MVF” Stokes.
el estudio
-Modelado de -Exponer ante
computacion 4.3
sistemas grupo la
al de la Procedimiento
difusivos- solución de
mecánica de generar de
convectivos, en problemas de
fluidos, análisis por el
una y dos difusión pura
basándose MVF.
dimensiones por utilizando el
en el MVF
4.4 Modelado MVF. MVF.
de problemas
-Simulación -Exponer la
de flujo de
numérica de solución de los
fluidos por
dinámica de problemas de
MVF.
fluidos utilizando difusión-
4.5 Simulación software convección, por
de flujo de comercial basado el MVF.
fluidos en MVF.
utilizando
Departamento de Desarrollo Curricular
 UNIVERSIDAD AUTÓNOMA DE COAHUILA
DIRECCIÓN DE ASUNTOS ACADÉMICOS
Subdirección de Superación Académica

Software - Desarrollo de
basado en modelos de
MVF. simulación
numérica en
Software
basado en MVF.

8.- Procedimiento de Evaluación:

Porcentaje
Tipo de evaluación Desarrollo del 90
Conocimiento
50
Exámenes parciales (promedio global, con parciales)
Examen Final 5
Tareas 30
Proyectos 5
Participación en el aula

Desarrollo de Habilidades 5
Trabajo en equipo
Comunicación oral y escrita
Planteamiento y solución de problemas

Desarrollo de Actitudes 5
Responsabilidad
Colaboración
Compromiso

TOTAL 100

BIBLIOGRAFÍA:

1. Finite Element Modeling and Simulation with Ansys Workbench. Xiaolin C. & Yijun L.. Ed. CRC Press. 2015.
2. Engineering Analysis with Ansys Software. Stolarski T., Nakasone Y. & Yoshimoto S. L. Ed. ELseviern Butterworth-
Heinemann. 2006.
3. Finite Element of Heat Transfer. Bergheau J. M & Fortunier R. Wiley 2008.

Departamento de Desarrollo Curricular

 UNIVERSIDAD AUTÓNOMA DE COAHUILA
DIRECCIÓN DE ASUNTOS ACADÉMICOS
Subdirección de Superación Académica

4. An Introduction to Computational Fluid Dynamics, The Finite Volume Method. Versteeg H. K. & Malalasekera W.
Segunda edición. Pearson 2007.
5. Técnicas Numéricas en Ingeniería de Fluidos. Introducción a la Dinámica de Fluidos Computacionales (CFD) por el
Método de Volúmenes Finitos. Fernández Oro J. M. Ed. Reverte. 2012.

Departamento de Desarrollo Curricular

ÍNDICE DE CUADROS 12

Figura 1: Pensamiento de Isaac Newton, padre de la Mecánica Clásica.

 C AP ÍTULO P RIMERO

I NTRODUCCI ÓN AL M ÉTODO DE

ELEMENTOS FINITOS

Objetivo Estudiar los conceptos fundamentales de la formulación directa en el cálcu-

lo de elementos mecánicos elásticos. Estudiar los conceptos primarios del MEF

Introducción Como introducción al tema, se empezará con la definición de un mode-

lo como una idealización de la realidad, una abstracción que permite si representación
a través de expresiones matemáticas. Haciendo una analogı́a neurolingüı́stica, un mo-
delo podrı́a considerarse como un mapa que representa el territorio que forma parte de
la realidad. A mayor calidad en los detalles ofrecidos por el mapa, una mejor repre-
sentación de la realidad, y, en consecuencia, una más fácil orientación al recorrer eo
territorio.
En términos generales, la complejidad de un modelo matemático esta directamente
relacionada con el nivel de fidelidad con el cual se intenta llevar a cabo la descripción de
la realidad. Existen, sin embargo, una vriedad de modelos simples pero suficientemente
representativos de los fenómenos naturales, que se utilizan con frecuencia para efectuar
cálculos y simulaciones en el ámbito de las ciencias e ingenierı́as. De hecho, muchos
modelos matemáticos sencillos han representado la realidad con gran fidelidad a lo
largo del tiempo, demostrando un alto grado de robustez y repetitividad en su rango de
aplicación. Algunos ejemplos de tales modelos son:
1. La ley de Hooke generalizada (elasticidad lineal) para materiales sólidos.
2. Las leyes de Fourier y de Newton para describir los flujos de calor por conducción
y convección por el interior y el contorno de cuerpos sólidos.
3. La ley constitutiva de Newton para fluidos viscosos newtonianos
4. La ley de Darcy de flujo en medios porosos
Los resultados obtenidos a partir de estos modelos le han dado forma en gran medida,
al mundo en el que vivimos, y han servido como puntos de referencia para el desarrollo
de nuevos modelos más complejos.
La existencia de una modelo está asociada unı́vocamente a la presencia del error. Es-
te puede deberse a múltiples razones, que van desde la compresión incompleta del

13
 14

fenómeno fı́sico, al despreciar y/o descartar variables que son esenciales, hasta caer en
un exceso de simplicidad a la hora de describir matemáticamente la realidad.
Otras fuentes de error se encuentra intrı́nsecamente contenida en el método estrategia
utilizados para llevar a cabo la resolución del modelo matemático. Este punto en par-
ticular será objeto de análisis y discusión a lo largo del curso, a la hora de abordar
las problemáticas y limitaciones inherentes a las metodologı́as analı́ticas o numéricas
disponibles para resolver modelos matemáticos.
El modelo de una determinado fenómeno fisico puede llevarse a cabo desde diferentes
enfoques. La elección de uno u otro dependerá de los objetivos perseguidos, del grado
de detalle requerido , del tiempo disponible para efectuar el análisis, etc. Un caso con-
creto, donde puede apreciarse la existencia de diferentes estrategias de modelado, es el
análisis de un mecanismo o de un sistema estructural, siendo este considerado como
sistema discreto o como sistema continuo.
Los sistemas discretos son aquellos en los que es posible reconocer a simple visita a
cada uno de los miembros o elementos que lo componen. Algunos ejmplos de este tipo
de sistemas son:
1. Circuitos eléctricos con componentes tales como resistencias, capacitores, ba-
terı́as, etc.
2. Redes de conducción de fluidos, donde se observan tramos y accesorios de cañerı́a
de diferentes caracterı́sticas, depósitos de reserva, ,fuentes y sumideros de flui-
dos, etc.
3. Sistemas estructurales reticulados formados por barras, vigas conectores elásti-
cos etc.
Por otra parte, los sistemas continuos son aquellos que representan porciones de materia
en estado sólido, lı́quido o gaseoso. Estos sistemas, a escala macroscópica, presentan
una apariencia continua y regular, ocultando ası́ las caracterı́sticas discontinuas de la
escala molecular de la materia a nivel microscópico.
La descripción matemática de los sistemas continuos se encuadra dentro de la mecáni-
ca de medios continuos. Esta disciplina se caracteriza por el estudio de las variables
macroscópicas de un medio (variables de campo). las cuales describen el comporta-
miento fenomenológico de la materia por medio de funciones continuas. Tales mode-
los matemáticos están representados fundamentalmente por un conjunto de ecuaciones
diferenciales ordinarias (ODE) o por ecauciones en derivadas parciales (PDE), depen-
diendo esto del número de dimensiones del dominio espacial bajo análisis.
Los modelos matemáticos descritos por ecuaciones diferenciales quedan completamen-
te definidos al establecer las respectivas condiciones de contorno e iniciales en caso de
tratarse de un problema transitorio o de evolución temporal. En consecuencia, las ecua-
ciones diferenciales tienen validez en el interior del dominio de análisis, mientras que
las condiciones de borde hacen lo propio en el contorno que rodea al dominio.
En caso de tratarse de un problema de evolución, las condiciones de borde deben com-
pletarse con las condiciones iniciales del problema, que se definen en la totalidad del
dominio para instante de tiempo inicial. Dichas condiciones alcanzan a todas aquellas
variables o magnitudes fundamentales que describen al problema.
1.1. MÉTODO NUMÉRICO EN LA INGENIERÍA MECÁNICA 15

Como consecuencia de ello, el análisis de sistemas continuos, también llamados de

parámetros distribuidos en la dimensión del espacio, requiere de un proceso primario
de discretización del dominio espacial en una cantidad finita de elementos o celdas;
sobre las cuales, a posteriori, se aplican las ecuaciones de equilibrio que surgen del
tratamiento matemático realizado sobre las ecuaciones diferenciales que describen el
fenómenos fı́sico. Estas pueden describir fenómenos de tipo térmico, mecánico, elec-
tromagnético y quı́mico, entre otros.
En caso de que el problema analizado sea transitorio, tanto para sistemas discretos
como continuos, es preciso efectuar adicionalmente la discretización del dominio tem-
poral, en el cual se va a desarrollar la evolución del proceso.
Frecuentemente, los modelos discretos se utilizan para llevar a cabo un análisis glo-
bal de sistemas extensos, o para utilizar sus resultados como condición de borde para
la resolución de problema más detallado, basado en modelos continuos. Un ejemplo
práctico de estos último es la evolución a la fatiga de la soldadura circunferencial de
una cañerı́a, en cuyo caso se realiza, en primer lugar, un análisis de tensiones a novel
global de trazado completo de la cañerı́a, a fin de obtener las condiciones de carga y/o
de desplazamientos y rotaciones en las inmediaciones de la soldadura; para más tarde
aplicárselas a un modelo continuo de un tramo de cañerı́a destinado a realizar el estudio
detallado de la fatiga en la soldadura.

1.1. Método numérico en la ingenierı́a mecánica

Las ecuaciones diferenciales pueden resolverse por medio de diferentes metodologı́as,
que se van desde las puramente analı́ticas hasta las completamente numéricas. Por lo
tanto, ante la variedad de metodologı́as y estrategias de resolución existentes, cabe lle-
var a cabo una breve clasificación y descripción de cada una de ellas, para generar crite-
rios de selección en función de los objetivos perseguidos y de los recursos de cálculo y
tiempo disponible. En tal sentido, el esquema conceptual de la figura 1.1 brinda una de
las posibles clasificaciones de los principios métodos que se utilizan para la resolución
de ecuaciones diferenciales.
Los principales grupos de métodos surgen de la clasificación de las soluciones resul-
tantes en exactas o aproximadas, y dentro de estas últimas, en las de tipo analı́tico
o numérico. Cabe destacar que en la actualidad, y como consecuencia de su elevada
versatilidad, robustez y grado de automatización, las soluciones aproximadas de tipo
numérico han ganado terreno por sobre las soluciones exactas y analı́ticas. Claros ejem-
plos de ello son las ampliamente utilizados métodos de diferencias finitas, volúmenes
finitos y elementos finitos para la resolución de problemas de interés tecnológico.
Las metodologı́as numéricas tienen en común el hecho de que las ecuaciones a resolver
deben ser discretizadas en el dominio espacial y/o temporal, según sea el caso. Sin em-
bargo, las expresiones a discretizar y resolver presentan ciertas particularidades propias
de cada metodologı́a.
En el caso del método de diferencias finitas (FDM), la ecuación diferencial que go-
bierna al problema se discretiza de manera directa a través de aproximaciones de las
derivadas en forma de diferencias. Por lo tanto, el proceso de discretización se redu-
ce a reemplazar las ecuaciones en diferencias en la expresión diferencial de gobierno
1.1. MÉTODO NUMÉRICO EN LA INGENIERÍA MECÁNICA 16

del problema. Mientras tanto, su posterior resolución se basa en la aplicación de las

expresiones discretas en todos los puntos interiores del dominio discretizado: De este
proceso surge un sistema de ecuaciones, en el que, luego de imponer las condiciones
de contorno e iniciales, es posible obtener la solución del problema.
Las principales desventajas del método de diferencias finitas son, por un lado, las difi-
cultades asociadas a la imposición de las condiciones de contorno naturales del proble-
ma, y, por otro lado, la reducida flexibilidad de las tı́picas discretizaciones estructurales
(cuadrı́culas) para adaptarse a dominios de geometrı́a compleja e irregular.

Figura 1.1: Esquema conceptual sobre los métodos utilizados para la resolución de ecuaciones
diferenciales.

El método de diferencias finitas se utiliza con mayor frecuencia para la resolución de

problemas de transferencia de calor y de análisis estructural en sólidos, de flujo poten-
cial y de difusión-reacción en fluidos, de electromagnetismo y de todos aquellos casos
donde sea necesario llevar a cabo la rápida discretización de dominios geométricamente
simples, sin importar el grado de complejidad de la ecuación diferencial a resolver.
Por su parte, el método de volúmenes finitos (FVM) representa el concepto de balance
o conservación en celdas que pueden ser irregulares. Sobre dichas celdas se plantea la
1.1. MÉTODO NUMÉRICO EN LA INGENIERÍA MECÁNICA 17

variación volumétrica (volumen de control) de una determinada variable, en términos

de los flujos que se producen a través de sus contornos (superficies de control). Por sus
caracterı́sticas intrı́nsecas, esta metodologı́a se aplica fundamentalmente a problema de
fluidos gobernados por flujos convectivos o problemas que involucran flujos de calor.
La discretización se efectúa de manera directa sobre la ecuación diferencial de gobierno
del problema, para la cual se consideran puntos del interior de la celda para considerar
los términos volumétricos y puntos sobre las caras de la misma, para estudiar los térmi-
nos convectivos. A través de la aplicación este proceso a un número finito de celdas, y
de la imposición de las condiciones de contorno e iniciales, se alcanza la solución del
problema.
El método de los elementos finitos (FEM) es la metodologı́a numérica actualmente
más utilizada para resolver una gran variedad de problemas perteneciente a la ciencias
y a la fı́sica. A diferencia de los métodos enumerados anteriormente, antes de proceder
con la discretización, el FEM requiere de un proceso de transformación de la ecuación
diferencial de gobierno del problema en una forma integral equivalente.
En el lenguaje utilizado normalmente en el ámbito de los métodos numéricos, un pro-
blema planteado en términos de su ecuación diferencial se denomina forma fuerte del
problema; mientras tanto, si el mismo problema está expresado de manera integral, se
denomina forma débil o variacional del problema.
La forma débil implica un enfoque energético o variacional del problema, por medio
del cual se relajan las condiciones impuestas por la ecuación diferencial en cada punto
del dominio hacia una solución generalizada, con la que se tiene a minimizar el error
ponderado por medio total en el dominio.
De acuerdo con las diferentes maneras de escribir el modelo matemático asociado a
un determinado problema, y en función de los explicado previamente, se observa que
los métodos de diferencias finitas y volúmenes finitos se aplican sobre la forma fuerte
de un determinado problema, mientras que el método de los elementos finitos lo hace
sobre la forma débil de este.
El FEM tiene como ventaja ser extremadamente flexible para ser aplicado en dominios
geométricamente irregulares y permitir una sencilla imposición de las condiciones de
contorno naturales. La discretización del dominio se lleva a cabo con particiones que,
en general, tomán de triángulos o cuadriláteros en 2D y de tetraedros o hexaedros 3D.
Una caracterı́stica y propiedad que tienen en común los métodos de diferencias fini-
tas, volúmenes finitos y elementos finitos es que, a medida que el dominio espacial se
discretiza en una cantidad mayor de particiones, celdas o elementos, el resultado se
aproxima más a la solución exacta del problema.
Esta propiedad se denomina comúnmente convergencia de la malla, y es objeto de
análisis en algunos de los problemas de aplicación mecánica.
El concepto de convergencia de la discretización espacial se extiende también a la dis-
cretización temporal, para la cual un mayor grado de subdivisión implica una mayor
precisión y calidad del resultado en el dominio temporal; esto se denomina comúnmen-
te convergencia del incremento de tiempo.
1.1. MÉTODO NUMÉRICO EN LA INGENIERÍA MECÁNICA 18

1.1.1. Evolución del método de elemento finito

Los elementos finitos se utilizaron por primera vez en la década de los años cincuenta
en el diseño de un aeroplano: Los principios fundamentales del MEF se conocen desde
el siglo XIX, pero la matemática que lleva implı́cito el método y la necesidad de la
computación conllevó a que se comenzará a aplicar en forma global en la década del
70 en las distintas oficinas de diseño. El método tuvo en sus inicios un mayor impacto
en la ingenierı́a civil, la industria automovilı́stica, la industria aeroespacial y el diseño
de buques. La generación del método, el cual no está limitado sólo al análisis tensional
y deformacional, ha permitido solucionar un número grande de problemas de potencial
eléctrico, flujo de fluidos, conducción de calor, procesos de conformación y problemas
no lineales en general.
La aplicación que tiene el MEF con sus paquetes de programas profesionales en las
industrias grandes y pequeñas, oficinas de diseño de los paı́ses productores de equipos
y tecnologı́as, implica la necesidad cada vez mayor de ingenieros conocedores de esta
técnica por lo que aparece cada vez con mayor frecuencia dentro de los planes de es-
tudio y programas de la enseñanza universitaria, lo que inidiscutiblemente adecuado a
las necesidades y exigencias de cada paı́s y universidad.
En los problemas de ingenierı́a mecánica relacionados con el diseño mecánico estruc-
tural, el problema consiste en estudiar las tensiones y deformaciones que ocurren en el
interior de un elemento estructural (biela, pieza bajo sistema de cargas) o sucesión de
cuerpos o elementos estructurales (bloque de un motor, vagón de ferrocarril) originado
como consecuencia de acciones externas: dinámicas o estáticas o las combinaciones de
ambas.
El MEF ha adquirido una gran importancia en la solución de problemas ingenieriles,
fı́sicos, etc., ya que permite resolver casos que hasta hace poco tiempo eran práctica-
mente imposibles de resolver por métodos matemáticos tradicionales.
Esta circunstancia obligada a realizar prototipos, ensayarlos e ir realizando mejoras de
forma iterativa, lo que traı́a consigo un elevado costo tanto económico como un mayor
tiempo de desarrollo del producto.
La concepción del MEF permite realizar un modelo matemático de cálculo del sistema
real, más fácil y económico de modificar que un prototipo. Sin embargo, no deja de ser
un método aproximado de cálculo debido a la hipótesis básicas del método. Los pro-
totipos, por tanto, siguen siendo necesarios, pero en menor número, ya que el primero
puede acercarse bastante al diseño óptimo.
La estructura básica del MEF en los últimos años ha sufrido un gran desarrollo debido
a los avances computacionales. Han sido precisamente estos adelantos tecnológicos lo
que han puesto a disposición de los clientes gran cantidad de programas que permiten
realizar cálculos con elementos finitos. El manejo correcto de este tipo de programas
exige un profundo conocimiento no solo de la temática que se trabaja, si no también de
los principios del MEF. De esta manera los ingenieros o profesionales especializados
en general están en condiciones de garantizar que los resultados obtenidos se ajusten a
la realidad.
1.1. MÉTODO NUMÉRICO EN LA INGENIERÍA MECÁNICA 19

Breve historia
Este método ha evolucionado mucho en las décadas desde su aparición, durante las
décadas de los 40. En el siguiente párrafo se recogen las contribuciones de ingenieros,
matemáticos y fı́sicos que han hecho posible el desarrollo y el estudio de dicha técni-
ca. El origen del método del elemento finito se remonta hacia los primeros años del
siglo XX, cuando algunos investigadores aproximaron y modelaron un medio elástico
continuo discretizando barras elásticas equivalentes. Courant (1943) ha sido reconoci-
do como la primera persona que desarrollo el método de elemento finito, por medio
de una publicación, en la cual utiliza iterpolaciones polinomiales sobre subregiones
triangulares para investigar problemas de torsión. En los años 50, la empresa Boeing
ya empleaba elementos triangulares de esfuerzos para simular las alas de los aviones.
Aunque, el método del elemento finito comenzó a ganar popularidad hacia los años
60, fue hasta el año 1967 cuando se publicó el primer libro especializado en el tema
por Zienkiewicz y Cheung. El método de elemento finito es un procedimiento numérico
que puede emplearse para solucionar una gran cantidad de problemas de ingenierı́a, por
ejemplo, en transferencia de calor, análisis de esfuerzos, electromagnetismo, mecánica
de fluidos, acústica, etc. En transferencia de calor por conducción, el método se emplea
para determinar la distribución de temperatura a través de un cuerpo con condicio-
nes lı́mites conocidas y las propiedades termo fı́sicas del material en estado estable o
transiente (señal en forma de onda que empieza con amplitud cero). En mecánica de
fluidos se puede solucionar problemas en estado estable o transiente de fluidos com-
presibles o incompresibles. Los análisis estructurales permiten encontrar la solución
de esfuerzos y desplazamientos de estructuras o piezas mecánicas sometidas a diver-
sas cargas. En general existen muchas áreas de aplicación y lo más importante es que
se pueden encontrar soluciones a problemas que combinan condiciones de operación
diferentes mediante un solo análisis.
Ahora el MEF ha alcanzado ya un grado de madurez que no se presta a progresos
espectaculares como los vividos en las décadas anteriores. Las lı́neas de investigación
a lo largo de las cuales se está desarrollando el MEF en la actualidad son:
1. Adaptación de algoritmos de cálculo a las nuevas arquitecturas de las compu-
tadoras, con objeto de aumentar la velocidad de cálculo y, por tanto, el tamaño
máximo de los problemas que se van a abordar
2. Desarrollo de medidas de error, mallados autoadaptivos y elementos de altas pres-
taciones, con objeto de aumentar la precisión y fiabilidad de los resultados obte-
nidos por clientes inexpertos en entorno de cálculo integrado con el CAD.
3. Modelos numéricos de leyes de comportamiento de materiales, sobre todo para
la predicción de fallo y para la representación del comportamiento de nuevos
materiales.
Actualmente el método se encuentra en una fase de gran expansión: es ampliamente
utilizado en la industria y continúan apareciendo cientos de trabajos de investigación
en este campo. Las computadoras has aportado el medio eficaz de resolver la multitud
de ecuaciones que se plantean en el MEF, cuyo desarrollo práctico ha ido caminando
parejo de las innovaciones obtenidas en el campo de la arquitectura de las computado-
ras. Entre éstas, además de permitir la descentralización de los programas de elementos
finitos (EF), ha contribuido a favorecer su uso a través de sofisticados paquetes gráficos
1.2. CONCEPTOS BÁSICOS DEL MEF 20

que facilitan el modelado y la sı́ntesis de resultados.

1.2. Conceptos básicos del MEF

Sin el empleo de las técnicas numéricas, probablemente muchos de los problemas de
ingenierı́a serı́an muy complicados de resolver con un buen grado de aproximación. La
mayorı́a de las técnicas numéricas en mecánica del medio continuo se basan en el prin-
cipio de que es posible derivar ecuaciones y relaciones que describen aproximadamente
el comportamiento de partes diferentes del cuerpo. Por medio de la división total del
cuerpo en proporciones pequeñas y usando algunas otras relaciones para ligar o ensam-
blarlas se puede obtener una predicción aproximada de los valores para las variables,
tales como esfuerzos y deformaciones sufridas en el cuerpo.

1.2.1. Proceso de Discretización

La discretización implica la aproximación de la realidad a un medio continuo. Se Em-
plea un número de términos para esquematizar la discretización, tales como; subdivi-
sión, continuidad, compatibilidad, convergencia, lı́mite superior e inferior, error, poten-
cial estacionario y mı́nimo residual.

Subdivisión
El espacio es finito e infinitamente divisible, además se sabe que para cada cosa existe
una magnitud. La figura 1.2 muestra el concepto de subdivisión en un espacio finito. Por
ejemplo, la tierra se encuentra contenida en el espacio inmenso, sin embargo, esta ocupa
un espacio finito dentro de ese espacio que puede subdividirse en un gran númemro de
partes. Otros ejemplos, es la subdivisión de un triángulo en un número de componentes
infinitos.

Figura 1.2: Divisibilidad: a) espacio infinito y b) triángulo divisible infinitamente.

1.2. CONCEPTOS BÁSICOS DEL MEF 21

Continuidad
Una cantidad esta formada de elementos divisibles. Por ejemplo, entre los dos puntos
de una recta, existen otros puntos, ası́ como, otros momentos existen entre dos mo-
mentos de un periodo de tiempo. Por lo tanto, el espacio y el tiempo son continuos e
infinitamente divisibles, y las cosas son consecutivas, contiguas y continuas. Figura 1.3

Figura 1.3: Conceptos de continuidad: a) consecutivos, b)contiguos y c) continuo

Esta ideas de espacios finitos, divisibilidad y continuidad permiten dividir las cosas
continuas dentro de pequeños componentes, unidades o elementos,

Convergencia
Para determinar el valor aproximado de π, o el área de un cı́rculo, se puede dibujar un
polı́gono dentro y fuera de cı́rculo, tal como se observa en las figuras siguientes.

(a) Convergencia y lı́mite para aproximar(b) Convergencia y lı́mite para aproximar

el área de un cı́rculo con polı́gonos inscri-el área de un cı́rculo con polı́gonos cir-
tos cunscritos

Figura 1.4: Concepto de convergencia

1.2. CONCEPTOS BÁSICOS DEL MEF 22

Si se comienza a incrementar el número de lados de los polı́gonos, estos serán cada

vez más pequeños, lo cual tiende a aproximarse a la circunferencia o área del cı́rculo.
Este proceso de moverse sucesivamente hacia la solución exacta más aproximada, se le
llama convergencia.
En el caso del cı́rculo de las figuras anteriores 1.4, la convergencia implica que como
el polı́gono interior y exterior incrementa su número de lados, entonces se aproxima o
converge al área del cı́rculo.
La figura 1.6, representa el gráfico de la aproximación sucesiva para los valores del área
del cı́rculo para polı́gonos internos y externos.

Figura 1.5: Convergencia para aproximar el área de un cı́rculo

Lı́mites o fronteras
Dependiendo el curso de acción que se tome desde el interior o exterior, se puede apro-
ximar a la solución del área del circulo. Sin embargo, el resultado para cada método
será diferente. La figura 1.5 muestra que la convergencia a partir de los polı́gonos ex-
ternos dan un valor más aproximado que cuando se aproxima con polı́gonos internos.
Lo cual indica que las condiciones lı́mite juegan un papel importante en el proceso de
aproximación, tal como se observa en la figura 1.6

Figura 1.6: Concepto de lı́mite

1.2. CONCEPTOS BÁSICOS DEL MEF 23

Error
Es evidente que la discretización involucra una aproximación. Consecuentemente, el
resultado obtenido no es una solución exacta, pero representa una aproximación a la
solución. A esta cantidad en la cual se difiere se le llama error.
Por ejemplo, las áreas (o perı́metros) de los polı́gonos inscritos en los cı́rculos de la
figura son siempre menores que el área (o perı́metro) del cı́rculo, y las áreas (o perı́me-
tros) de los polı́gonos circunscritos son siempre mayores que el área (o perı́metro) del
cı́rculo.

Principios y leyes
Para describir el comportamiento de un fenómeno o comportamiento fı́sico, es nece-
sario establecer leyes basadas en principios. Una ley puede establecerse a través de
un enunciado o por medio de una expresión matemática. Los principios por lo gene-
ral son propuestos a través de la intuición, los cuales son comprobados a través de la
experimentación. Ası́ por ejemplo, la segunda ley de Newton, establece:

d2 u
F = ma = m (1.1)
dt2
Donde F es la fuerza, m es la masa y a es la aceleración o la segunda derivada de
desplazamiento u con respecto de t.
El principio es que para el tiempo dado, el cuerpo esta en equilibrio dinámico y una
medida de energı́a contenida en el cuerpo supone un valor estacionario.

1.2.2. Ventajas y limitaciones del MEF

Como todo método aproximado de solución de problemas de ingenierı́a, el método de
los elementos finitos posee ventajas y limitaciones que deben ser conocidas por los
profesionales que pretenden aplicar el mismo

Ventajas
El MEF puede manejar elementos de geometrı́a y restricciones complejas.
El MEF puede trabajar con solicitudes complejas de carga desde cargas pun-
tuales, uniformemente distribuidas y dinámicas. El MEF puede resolver varios
problemas de ingenierı́a desde mecánica de sólidos, dinámica, transferencia de
calor, campos electrostáticos y problemas no lineales.

Limitaciones MEF
Las soluciones por el MEF a menudo son aproximadas. No obstante, a más refi-
namiento de la malla, mayor exactitud en las soluciones.
Las soluciones del MEF pueden contener errores inherentes a la computación
como resultado de la acumulación de errores durante la computación numérica.
Las soluciones del MEF pueden contener errores fatales como resultado de una
incorrecta modelación de estructuras, cargas ó de las condiciones lı́mites.
1.2. CONCEPTOS BÁSICOS DEL MEF 24

Integración de disciplinas para MEF

La Mecánica Computacional representa la integración de varias disciplinas, como se
muestra en cuadro anterior. Se pueden distinguir varias ramas de la mecánica compu-
tacional según la escala fı́sica
 del foco de atención:


 Nanomecánica
 Micromecánica


 
Mecánica computacional  Sólidos

 Mecánica continua Fluidos
 
Multifisı́ca




Sistemas
La nanomecánica se ocupa de fenómenos a nivel molecular y atómico. Como tal, está
estrechamente relacionado con la fı́sica y la quı́mica de partı́culas. En la escala atómica
se pasa a la mecánica cuántica. La micromecánica analiza principalmente los niveles
cristalográficos y granulares de la materia. Su principal aplicación tecnológica es el
diseño y fabricación de materiales y microdispositivos.
La mecánica continua estudia los cuerpos a nivel macroscópico, utilizando modelos
continuos en los que la microestructura se homogeneiza mediante promedios fenome-
nológicos. Los dos campos de aplicación tradicionales son la mecánica de sólidos y de
fluidos. La mecánica estructural es una rama conjunta de la mecánica de sólidos, ya
que las estructuras, por razones obvias, se fabrican con sólidos. La mecánica compu-
tacional de sólidos favorece un enfoque de ciencias aplicadas, mientras que la mecánica
estructural computacional enfatiza las aplicaciones tecnológicas al análisis y diseño de
estructuras.
La mecánica de fluidos computacional se ocupa de problemas que involucran el equili-
brio y el movimiento de lı́quidos y gases. Las subáreas relacionadas bien desarrolladas
son la hidrodinámica, la aerodinámica, la fı́sica atmosférica, la propulsión y la combus-
tión.
La multifı́sica es un recién llegado. Esta área pretende incluir sistemas mecánicos que
trasciendan los lı́mites clásicos de la mecánica de sólidos y fluidos. Un ejemplo clave es
la interacción entre fluidos y estructuras, que tiene importantes subáreas de aplicación
como la aeroelasticidad y la hidroelasticidad. Los problemas de cambio de fase, como
el derretimiento del hielo y la solidificación de metales, encajan en esta categorı́a, al
igual que la interacción de los sistemas de control, mecánicos y electromagnéticos.
Finalmente, el sistema identifica objetos mecánicos, ya sean naturales o artificiales,
que realizan una función distinguible. Ejemplos de sistemas creados por el hombre son
aviones, edificios, puentes, motores, automóviles, microchips, radiotelescopios, robots,
patines y aspersores de jardı́n. Los sistemas biológicos, como la ballena, la ameba, el
virus o el pino, se incluyen si se estudian desde el punto de vista de la biomecánica.
Las entidades ecológicas, astronómicas y cosmológicas también forman sistemas.

Estático vs Dinámico
Los problemas de mecánica continua se pueden subdividir según se tengan en cuenta o
no los efectos inerciales:
1.3. PROCEDIMIENTO GENERAL DEL MEF 25
 (
Tiempo invariante


Mecánica continua Estático
 Cuasi estático

Dinámico
En estática las fuerzas de inercia se ignoran o desprecian. Estos problemas pueden
subclasificarse en invariantes en el tiempo y cuasiestáticos. Para el primer tiempo no
es necesario considerarlo explı́citamente; cualquier parámetro de ordenación de res-
puestas similar al tiempo (si es necesario) servirá. En problemas cuasiestáticos como
asentamientos de cimientos, flujo de fluencia, plasticidad dependiente de la velocidad
o ciclos de fatiga, se requiere una estimación más realista del tiempo, pero las fuerzas
de inercia se ignoran siempre que los movimientos sigan siendo lentos.
En dinámica, la dependencia del tiempo se considera explı́citamente porque el cálculo
de las fuerzas de inercia (y/o de amortiguación) requiere que se tomen derivadas con
respecto al tiempo real.

Lineal vs no lineal
Una clasificación
( de problemas estáticos puede ser

Estático Lineal
No lineal
El análisis estático lineal se ocupa de problemas estáticos en los que la respuesta es
lineal en el sentido de causa y efecto. Por ejemplo: si se duplican las fuerzas aplicadas,
también se duplican los desplazamientos y las tensiones internas. Los problemas fuera
de este dominio se clasifican como no lineales.

1.3. Procedimiento general del MEF

El MEF es una herramienta de análisis numérico que se utiliza para resolver en forma
aproximada problemas de borde de sistemas discretos, tales como las estructuras, y
medios continuos, como son los campos escalares y vectoriales. Dicho método es una
técnica más eficiente y precisa que el método de diferencia finita dado que este se adapta
de mejor manera a los dominios irregulares. Por otra parte, en el MEF se discretiza el
dominio del problema mientras que en el método de diferencia finita se discretiza la
ecuación diferencial. Actualmente, se considera al método de las diferencias finitas
como una subclase del MEF y de hecho se puede demostrar que el MEF se reduce al
método de diferencia finita cuando las mallas son regulares.
La premisa básica del MEF, es que una región (dominio) donde se resuelve una ecua-
ción diferencial se modela analı́ticamente reemplazándola con un arreglo de elementos
discretos. Esto permite reducir un número infinito de incógnitas del problema a uno
con un número finito de incógnitas. Esto se hace dividiendo el dominio en elementos
y expresando las variables de campo incógnitas en términos de funciones aproximadas
dentro de cada elemento.
Las funciones aproximadas se expresan en términos de valores de variables de campo
para ciertos puntos, llamados nodos o puntos nodales. El conjunto de nodos y elementos
configura una malla o rejilla de solución para el problema. Esta malla puede o no seguir
la configuración fı́sica del campo. Por ejemplo, se puede trasladar el problema al campo
1.3. PROCEDIMIENTO GENERAL DEL MEF 26

de solución matemático, cuyas fronteras pueden coincidir con las orillas del cuerpo
fı́sico.
Existen básicamente cuatro maneras de formular las ecuaciones del sistema, que son:
1. Aproximación directa, las ecuaciones del sistema se ensamblan directamente de
las ecuaciones que gobiernan el problema. La desventaja de está formulación es
que solo se puede analizar elementos de formas o geometrı́as simples.
2. Aproximación variacional, en esta alternativa, se requiere usar funciones obte-
nidas del cálculo varriacional, es decir, encontrando los valores crı́ticos de una
funcional (por ejemplo la Energı́a potencial). Esta formulación tiene la ventaja
de que se puede usar en formas simples y complejas.
3. Aproximación de residuos ponderados, en esta solución se traslada el problema
del campo de solución fı́sico al campo de solución puramente matemático. Tie-
ne como ventajas que puede aplicarse donde no se encuentra con una funcional
adecuado.
4. Aproximación de balance de energı́a, descansa en el balance de energı́a térmica
o mecánica del sistema. Su ventaja es similar a la formulación anterior.
Los pasos a seguir, de manera general, para resolver un problema del medio continuo
usando el MEF son los siguientes:
1. Discretizar el dominio mediante elementos y nodos, es decir construir una malla
sobre el dominio.
2. Seleccionar y plantear las ecuaciones de los elementos (funciones de interpola-
ción apropiada)
3. Ensamblar las propiedades de los elementos para formar las ecuaciones del sis-
tema (trasladar las ecuaciones del esquema local, elemento por elemento, al es-
quema global del sistema).
4. Modificar las ecuaciones del sistema para tomar en cuenta las condiciones de
frontera.
5. Resolver el sistema de ecuaciones.
6. Estudio de convergencia y estimación de error.
7. Presentación gráfica de los resultados.

Cálculo de área de un cı́rculo usando el concepto de elementos finitos Para expli-

car de forma sencilla se calculará el área de un cı́rculo utilizando elementos triangulares
que forman un polı́gono inscrito en el cı́rculo de radio R. Los pasos a seguir son los
siguientes:
1.3. PROCEDIMIENTO GENERAL DEL MEF 27

(a) (b)

Figura 1.7: Discretización de un cı́rculo mediante un polı́gono inscrito

1. Discretización del dominio, la figura 1.7a muestra una discretización de la cir-

cunferencia, mediante un polı́gono inscrito formado por elementos triangulares.
Cada triángulo esta formado por tres nodos ubicados en los vértices.
Se debe destacar que la discretización se puede realizar también con elementos
rectangulares u otra forma geométrica.
2. Planteamiento de las ecuaciones de los elementos
La ecuación del elemento para calcular el área del triángulo se obtiene a partir de
la figura 1.7b donde el área del elemento (Ae ) es:

bh
Ae = (1.2)
2
como:

θ
b = 2R sin (1.3)
2

θ
h = R cos (1.4)
2
Se obtiene al sustituir la ecuación 1.3 y 1.4 en 1.3 que:

R2
Ae = sin θ (1.5)
2

3. Ensamble, modificación y resolución de las ecuaciones de los elementos

Para calcular el área total del cı́rculo se suman el área de cada elemento, es decir:

N
X
AN = Ae (1.6)
e=1

Como se tienen N elementos de área iguales, resulta que:

R2
AN = N A e = N sin θ (1.7)
2
1.3. PROCEDIMIENTO GENERAL DEL MEF 28

donde

2π
θ= (1.8)
N

N R2
 
2π
AN = sin (1.9)
2 N

4. Estudio de convergencia y estimación de error, el error absoluto se define como

la diferencia en valor absoluto de la solución exacta y solución aproximada por
elementos finitos esto es:

A − AN
kekN, % = × 100 (1.10)
A

donde kekN, % es el error cometido con una aproximación de N elementos, A es

la solución exacta y AN es la solución aproximada del área de la circunferencia.
Se debe notar que la solución exacta se alcanza cuando

A = lı́m AN ⇒ lı́m eN = 0, (1.11)

N →∞ N →∞

Para verificar la ecuación 1.11 se calcula el área exacta de la circunferencia A =

πR2

N A AN eN, %
3
6
9
12
15
18
20
30
40
50
100
150
200
300

Tabla 1.1: Valores numéricos de la solución exacta y aproximada.

Al graficar los valores de la tabla se observarı́a que cuando aumenta el número de

elementos al infinito el área converge al valor exacto de la solución (A = π). También
note que el error absoluto disminuye cuando el numero de elementos aumenta este se
aproxima a cero.
1.4. ENFOQUE DIRECTO DE ANÁLISIS POR MEF 29

1.4. Enfoque directo de análisis por MEF

El método directo de la rigidez brinda una herramienta inestimable a la hora de lle-
var acabo análisis simplificados del comportamiento global de estructuras o elementos
complejos.
En el presente tema se desarrollan las formulaciones para el análisis estático con ele-
mentos discretos de tipo resorte, barra y vigas; con ejemplos ilustrativos de la aplicación
del método directo en casos estructurales tı́picos, de torsión y transferencia de calor.

1.4.1. Elemento resorte lineal en 1D

El resorte lineal puede ser considerado como el elemento finito más simple, sin em-
bargo, detrás de su aparente simplicidad, se ocultan los conceptos fundamentales del
FEM. El resorte lineal presenta dos nodos vinculados entre si por el cuerpo resisten-
te del elemento,
 N que
 se representa matemáticamente por medio de la rigidez de valor
constante kr mm 2 ver la figura 1.8c). Su comportamiento mecánico se caracteriza por
la relación proporcional entre la deflexión de resorte δ = (u2 − u1 )[mm], que mide
su acortamiento o alargamiento, y la fuerza F [N ] necesaria para deformarlo ver figura
1.8b)

Figura 1.8: Concepto de lı́mite

El resorte lineal presenta dos grados de libertad dados por los desplazamientos nodales
u1 y u2 en la dirección de los elementos
Para deducir las ecuaciones que representan el comportamiento mecánico de l resor-
te lineal, es posible plantear el equilibrio de fuerzas usando la 2da ley de newton o
conservación del momento lineal), para el caso estático se obtiene

f1 + f2 = 0 (1.12)
f 1 = −f 2 (1.13)

donde f1 y f2 son las fuerzas correspondientes a cada uno de los nodos y grados de
libertad del resorte.
Por otra parte, de la definición de las fuerzas elásticas del resorte (ley constitutiva de
Hooke), se tiene que

f 1 = kr δ = kr (u1 − u2 ) (1.14)
f 2 = kr δ = kr (u2 − u1 ) (1.15)
1.4. ENFOQUE DIRECTO DE ANÁLISIS POR MEF 30

De lo que resulta el sistema de ecuaciones asociado a la condición de equilibrio estático

del resorte lineal en 1D
    
kr −kr u1 f
= 1 (1.16)
−kr kr u2 f2
   
kr −kr u1
Donde, la primera es la matriz de rigidez , ur = es el vector de
  −kr kr u2
f
desplazamiento nodal y fr = 1 es el vector de fuerzas nodales del elemento resorte
f2
lineal.

Problema 1 Considere una barra con sección transversal variable, que soporta una
carga P, como se muestra en la figura. La barra esta fija a uno de los extremos y en el
otro extremo se aplica la carga. Designamos el espesor de la barra como w1 en la parte
superior y como w2 en la parte inferior, con un espesor t homogéneo y una longitud L.
El material de la barra tiene modulo de elasticidad E. Aproximar la deformación de la
barra en varios puntos a lo largo de su longitud si supone una carga considerablemente
mayor al peso de la barra se puede despreciar su peso. Asuma que E = 10.46 inlb2 ,
w1 = 2in, w2 = 1in, t = 0.125in, L = 10in y P = 1000lb

(a)

Figura 1.9

Problema 2 Para la estructura formada por tres resortes y tres masas que se muestra
en la figura, determinar los desplazamientos de cada una de las masas y el valor de la
reacción de vı́nculo en el empotramiento superior.
1.4. ENFOQUE DIRECTO DE ANÁLISIS POR MEF 31

(a)

Figura 1.10

Problema 3 Una placa de acero esta sujeta a una carga axial, como se muestra en la
figura. Aproximar la deformación y el esfuerzo en la placa. La placa tiene 1/16in de
espesor y un modulo de elasticidad E = 296 inlb2 .

(a)

Figura 1.11
1.4. ENFOQUE DIRECTO DE ANÁLISIS POR MEF 32

1.4.2. Elemento resorte de torsión

El elemento resorte de torsión tiene un comportamiento fı́sico similar al del resorte
lineal, pero, a diferencia de este, su mecanismo de funcionamiento se basa en la apli-
cación de momentos antágicos en sus extremos, a fin de modificar el ángulo existente
entre ellos. Por lo tanto, el resorte de torsión presenta dos nodos, caracterizados cada
uno de ellos por el grado de libertad de rotación θ que se produce alrededor del eje Z,
que es perpendicular al plano XY como se muestra en la figura 1.12

Figura 1.12: Elemento resorte de torsión, a) esquema del elemento, momento y rotaciones no-
dales, b)relación lineal momento-rotación

Repitiendo el mismo procedimiento seguido anteriormente para el elemento resorte

lineal, pero considerando en este caso el grado de libertad de rotación, e imponiendo el
equilibrio de momentos nodales como m1 + m2 = 0, se tiene que

m1 = kt δ = kt (θ1 − θ2 ) (1.17)
m2 = kt δ = kt (θ2 − θ1 ) (1.18)

De lo que resulta el sistema de ecuaciones asociados al resorte de torsión:

    
kt −kt θ1 m1
= (1.19)
−kt kt θ2 m2
   
kt −kt θ
Donde, la primera es la matriz de rigidez de torsión , θt = 1 es el vector
  −k t kt θ2
m1
de rotaciones nodales y mt = es el vector de momentos nodales del elemento
m2
resorte de torsión y kt = JG
l

Problema 4 Considere un árbol circular sometido a torsión, que se muestra en la

figura. Recuerde de su anterior estudio de la mecánica de materiales que el ángulo de
torsión para un eje con área de sección transversal constante con momento polar de
inercia J y longitud L, hecha de material homogéneo con un módulo de rigidez G ,
sujeto a torsión T
1.4. ENFOQUE DIRECTO DE ANÁLISIS POR MEF 33

(a)

Figura 1.13

1.4.3. Elemento de conducción

En los análisis de transferencia de calor, una formulación de elementos finitos para
problemas de transferencia de calor puede desarrollarse. Las ecuaciones principales
que se involucran en la conducción es la ley de Fourier y la ley de enfriamiento de
Newton, considerando que la transferencia de calor en un sólido sucede desde una
medio o hasta un medio convectivo. Suponiendo un problema donde el dominio de
elementos empieza su discretizado, desde un ambiente de alguna atmósfera de fluido
móvil a determinada temperatura, donde el calor pasa a través de uno o varios tipos
sólido con una conductividad térmica dada.
Las ecuaciones se pueden escribir para la conducción cuando existe un gradiente térmi-
∂T
co en el medio, partiendo de la ley de Fourier en una dimensión como qx = −kA ∂X ,
donde qx es la componente X de la velocidad de transferencia de calor, k es la conduc-
∂T
tividad térmica del medio, A es el área normal al flujo de calor y ∂X es el gradiente
térmico.

Figura 1.14: Concepto de lı́mite

El signo menos de la ecuación significa que la temperatura disminuye en dirección del

flujo de calor. La ecuación puede reescribirse como una diferencia entre temperatu-
ra nodal entre dos nodos separados por una longitud correspondiente a una longitud
elemental, de acuerdo a la ecuación,
1.4. ENFOQUE DIRECTO DE ANÁLISIS POR MEF 34

kA(Ti+1 − Ti )
q= (1.20)
l

En el campo de la transferencia de calor, es común escribir la ecuación 1.20 en términos

de coeficiente global de transferencia de calor U donde U = kt , donde U representa la
transferencia de calor por unidad de área, en el sistema inglés U tiene unidades de
Btu/hr · f t2 · F , de esta manera

q = U A(Ti+1 − Ti ) (1.21)

En estado estable los elementos convectivos pueden ser modelados con la ley de en-
friamiento de Newton q = hA(Ts − Tf ) donde h es un coeficiente convectivo, Ts es la
temperatura en la superficie, y Tf representa la temperatura del fluido móvil. La ley de
Newton puede reescribirse en términos de U como

q = U A(Ts − Tf ) (1.22)

Donde U = h, y representa el inverso de la resistencia térmica debido a la condición de

frontera de la convección. En estado estable de conducción, la aplicación del balance
de energı́a a una superficie, con una transferencia de calor convectiva, requiere que la
energı́a transferida a la superficie por conducción sea igual a la energı́a transferida por
convección, esto significa que,

∂T
−kA = hA [Ts − Tf ] (1.23)
∂X

Figura 1.15: Concepto de lı́mite

Ahora podemos entender los dos modos de transferencia de calor envueltos en el pro-
blema, y aplicar el balance de energı́a entre las diferentes superficies, empezando desde
1.4. ENFOQUE DIRECTO DE ANÁLISIS POR MEF 35

la superficie exterior las pérdidas de calor por conducción que será igual a la pérdida
de calor producida por la corriente de aire, considerando tres nodos, el primero en la
parte fluido con Tf , el segundo en la interfase sólido-lı́quido con Ti = Ts y el tercero
con T = Ti+1 nuestro balance será

U2 A(T3 − T2 ) = U1 A(T2 − T1 ) (1.24)

En general, en los problemas de conducción, la velocidad de transferencia de calor qi

y qi+1 y las temperaturas nodales serán Ti y Ti+1 para cada elemento de acuerdo a la
ecuación,

kA
qi = (Ti − Ti+1 ) (1.25)
l
kA
qi+1 = (Ti+1 − Ti ) (1.26)
l

Figura 1.16: Concepto de lı́mite

Estas ecuaciones pueden representarse como temperaturas nodales,

    
qi kA 1 −1 Ti
=
qi+1 l −1 1 Ti+1

Ahora la matriz puede llamarse de conductancia y podrá escribirse como,

 
 (e) kA 1 −1
K =
l −1 1

Esta matriz puede escribirse como

 
 (e) 1 −1
K = UA
−1 1
1.5. ENFOQUE ENERGÉTICO DE ANÁLISIS POR MEF 36

Similarmente, bajo condiciones de estado estacionario, la aplicación de la conservación

de la energı́a en los nodos convectivos será:
    
qi Ti 1 −1
= hA
qi+1 Ti+1 −1 1

La matriz de conductancia puede escribirse,

 
 (e) Ti
K = hA
Ti+1

La ecuación anterior, puede escribirse en términos de U (U = h):

 
 (e) 1 −1
K = UA
−1 1

Problema 5 Una pared de una casa esta construida de los materiales mostrados en la
tabla de la imagen, si se asume una temperatura interior de 70°F y exterior de 20°F, con
una área expuesta de 150f t2 . Determine la distribución de temperaturas en el espesor.

Figura 1.17: Problema conceptual de transferencia de calor

1.5. Enfoque energético de análisis por MEF

Una manera alternativa de deducir el sistema de ecuaciones que representa el equilibrio
estático de un resorte lineal es haciendo uso del principio de mı́nima energı́a potencial.
Para ello, se parte de la definición de la energı́a potencial elástica de un resorte lineal:

1 1
er = kr δ 2 = kr (u2 − u1 )2 (1.27)
2 2
Y de la definición del trabajo de las fuerzas nodales se tiene que:

wr = u1 f1 + u2 f2 (1.28)

Luego, se define el funcional de trabajo y energı́a como:

1.5. ENFOQUE ENERGÉTICO DE ANÁLISIS POR MEF 37

πr = er − wr (1.29)

Finalmente, para encontrar la condición extrema, se impone el requisito de que la pri-

mera variación del funcional sea nula, es decir que dπr = 0, de lo que surge su mini-
mización:

∂πr
=0
∂u1
(1.30)
∂πr
=0
∂u2

En consecuencia, reemplazando las expresiones de la energı́a de potencia elástica y del

trabajo de las fuerzas nodales, y calculando posteriormente sus derivadas, se tiene que:

∂πr
= kr (u1 − u2 ) − f1 = 0
∂u1
(1.31)
∂πr
= kr (u2 − u2 ) − f2 = 0
∂u2

De este proceso, resulta un sistema de ecuaciones que es coincidente con el obtenido

por la aplicación del equilibrio de fuerzas en la Ec. 1.16

1.5.1. Estructuras 1D de resortes lineales

En los casos en los que se tienen conjuntos estructurales formados por más de un resor-
te, como el que se muestra en la Fig. 1.18, es preciso obtener el sistema de ecuaciones
que describe el equilibrio de fuerzas global de la estructura. Para ello, una posibilidad
es aplicar el principio de mı́nima energı́a potencial descrito en el apartado ??, el cual
para el sistema de la Fig. 1.18, resulta en el conjunto de expresiones de la Ec. 1.32, la
Ec. 1.33 y la Ec. 1.34

Figura 1.18: Conjunto estructural formado por dos resortes lineales.

π =e−w (1.32)

ne
X 1 1 1
e= ki δi2 = k1 (u2 − u1 )2 + k2 (u3 − u2 )2
i=1
2 2 2
1.5. ENFOQUE ENERGÉTICO DE ANÁLISIS POR MEF 38

donde ne = 2

nn
X
w= ui fi = u1 f1 + u2 f2 + u3 f3
i=1

donde nn = 3

1 1
π = k1 (u2 − u1 )2 + k2 (u3 − u2 )2 − (u1 f1 + u2 f2 + u3 f3 )
2 2
Luego, para determinar el comportamiento fı́sico del sistema, se encuentra la condición
estacionaria de mı́nima energı́a potencial al calcular la primera variación del funcional
de trabajo y energı́a dπ = 0,

∂πr
dπ = 0 → = 0 i = 1, . . . , nn
∂ui
∂π
= k1 (u1 − u2 ) − f1 = 0
∂u1
(1.33)
∂π
= k1 (u2 − u1 ) + k2 (u2 − u3 ) − f2 = 0
∂u2
∂π
= k2 (u3 − u2 ) − f3 = 0
∂u3

Finalmente, el sistema de ecuaciones resultante que representa el comportamiento mecáni-

ca del sistema de dos resortes de la Fig. 1.18 se puede expresar de forma matricial
como:
Si bien esta metodologı́a es posible aplicarla a cualquier sistema mecánico siguiendo
el procedimiento descrito, esto puede resultar en la práctica algo inabordable, o hasta
inviable, cuando la cantidad de elementos resortes presentes en el sistema es elevada.
Para solventar esta situación, dentro de la metodologı́a FEM se desarrolló un formato
automatizado del procedimiento denominado proceso de ensamblaje de las matrices de
rigidez y de los vectores de cargas elementales. Este proceso de ensamblaje permite
combinar los sistemas de ecuaciones elementales como los de la Ec. 1.16 o la Ec. 1.19,
obtenidos de manera independiente, en un solo sistema de ecuaciones global acoplado
como la de la Ec. 1.34
    
k1 −k1 0 u1  f1 
−k1 k1 + k2 −k2  u2 = f2 (1.34)
0 −k2 k2 u3 f3
   

Problema 6 Considerar las siguientes situaciones: a) Donde se aplica una fuerza F a

un resorte lineal como se muestra en la Fig. 1.19a. Dependiendo del valor de la rigidez
del resorte, el resorte se estira un cierto valor x. El equilibrio estático requiere que las
fuerzas aplicadas F sea igual a las fuerzas internas del resorte kx,

F
F = kx ó x = (1.35)
x
1.5. ENFOQUE ENERGÉTICO DE ANÁLISIS POR MEF 39

b) La barra mostrada en la figura pesa 8N y es soportada por un resorte con rigidez

k = 20N/cm cuando una fuerza P = 12N es aplicada en el punto C como se muestra
en la figura, calcular la deformación del resorte.

(a) (b)

Figura 1.19: a)resorte lineal sujeto a cargas, b) barra empotrada en A sujeta a carga en C
 C AP ÍTULO S EGUNDO

E SFUERZOS

Objetivo Desarrollar el MEF para diseñar elementos mecánicos, enfocados al estudio

de esfuerzos y deformación elástica en una, dos y tres dimensiones.

Introducción

2.1. Teorı́a básica y ecuaciones de elasticidad

2.2. Modelado de armaduras por MEF
2.3. Modelado de barras por MEF en una y dos dimen-
siones
2.4. Teorı́a básica para el análisis de vigas y marcos
2.5. Modelado de vigas y marcos por MEF en una y dos
dimensiones
2.6. Modelado en dos dimensiones de placas
2.7. Simulación de elasticidad en tres dimensiones de
placas utilizando software basado en MEF

40
 C AP ÍTULO T ERCERO

A N ÁLISIS T ÉRMICO

objetivo Desarrollar el MEF para diseñar elementos mecánicos, enfocados al estudio

de esfuerzos y deformación elásticas de origen térmico en una, dos y tres dimensiones.

Introducción

3.1. Teorı́a básica y ecuaciones de transferencia de ca-

lor
3.2. Modelado de campos térmico por MEF
3.3. Modelado de esfuerzos elásticos de origen térmico
por MEF
3.4. Simulación de transferencia de calor y esfuerzos
térmicos utilizando software basado en MEF

41
 C AP ÍTULO C UARTO

A N ÁLISIS DEL CAMPO

FLUIDODIN ÁMICO POR MVF

objetivo Desarrollar las bases de la metodologı́a numérica empleada en el estudio

computacional de la mecánica de fluidos, basándose en el MVF.

Introducción

4.1. Conceptos básicos del flujo de fluidos

4.2. Ecuaciones de Navier-Stokes
4.3. Procedimiento general de análisis por MVF
4.4. Modelado de problemas de flujo de fluidos por MVF
4.5. Simulación de flujo de fluidos utilizando software
basado en MVF

42
Bibliografı́a

43