INGENIERÍA ELÉCTRICA

MÉTODOS NUMÉRICOS

CATEDRÁTICO:

MOLINA FALCÓN JOSÉ FRANCISCO

INVESTIGACIÓN

UNIDAD 4 “INTEGRACIÓN NUMÉRICA”

AUTORES:

ESTRELLA PEREZ SAÚL

HERNÁNDEZ ROJAS JUVENCIO ANTONIO

LLINAS ACEVEDO JOSÉ ÁNGEL DE JESÚS

GRUPO:

12:00-13:00 pm
FECHA DE ENTREGA:

H. VERACRUZ, VER. A 30 DE MAYO DEL 2018

Contenido
4 INTEGRACIÓN NUMÉRICA ................................................................................ 3
4.1 FORMULAS DE INTEGRACIÓN DE NEWTON COTES................................ 3
4.1.1 REGLA TRAPECIAL. ............................................................................... 5
4.2 APLICACIONES DE LA INTEGRACIÓN NUMÉRICA. ................................... 9
4.3 USO DE HERRAMIENTAS COMPUTACIONALES ..................................... 12
4 INTEGRACIÓN NUMÉRICA

INTRODUCCIÓN

La integral resuelve el problema de calcular el ´área bajo la gráfica de una función

positiva definida sobre un intervalo cerrado. El cálculo elemental de funciones de
una variable real proporciona un método elegante de calcular la integral de una
función. El teorema fundamental del cálculo nos dice que el problema de calcular la
integral de una función continua se reduce al de buscar una segunda función cuya
derivada sea la función dada, es decir una primitiva de ella. Sin embargo, el
problema de hallar una primitiva de una función dada puede resultar muy difícil si
no imposible. De hecho, sabemos que existen funciones elementales —es decir,
combinaciones algebraicas de funciones trigonométricas y logarítmicas y sus
inversas— cuyas primitivas no son expresables de esta forma (p. ej. e −x 2). Por
esta razón es por la que estudiamos métodos numéricos que aproximen el valor de
la integral buscada. Ya la definición de integral de Riemann proporciona un método
de aproximación numérica: las sumas de Riemann. Sin embargo, su convergencia
es muy lenta y no resultan ´útiles para obtener resultados prácticos. Los métodos
numéricos que vamos a estudiar consisten en sustituir la función dada por una
aproximación suya y tomar como valor de la integral de la función el valor de la
integral de su aproximada.

4.1 FORMULAS DE INTEGRACIÓN DE NEWTON COTES.

Las fórmulas de Newton - Cotes son los tipos de integración numérica más
comunes. Se basan en la estrategia de reemplazar una función complicada o datos
tabulados por un polinomio de aproximación que es fácil de integrar:

donde fn(x) es un polinomio de la forma

donde n es el grado del polinomio. Por ejemplo, en la Figura 1 se utiliza un polinomio
de primer grado como una aproximación, mientras que en la Figura 2, se emplea
una parábola con el mismo propósito.

La integral también se puede aproximar usando un conjunto de polinomios

aplicados por pedazos a la función o datos, sobre segmentos de longitud constante.
Así, en la Figura 3, se usan tres segmentos de línea recta para aproximar la integral.
 Existen formas cerradas y abiertas de las fórmulas de Newton - Cotes. En
esta sección sólo se analizarán las formas cerradas. En ellas, se conocen los datos
al inicio y al final de los límites de integración.

Grado de precisión de una fórmula de integración numérica

El grado de precisión de una fórmula de integración numérica es el número

natural n que verifica que el error de truncamiento E[Pi]=0 para todos los polinomios
Pi(x) de grado i ≤ n, y existe un polinomio Pn+1(x) de grado n+1 tal que E[Pn+1]≠0.

4.1.1 REGLA TRAPECIAL.

La regla del trapecio es la primera de las fórmulas cerradas de Newton - Cotes.

Corresponde al caso donde el polinomio de aproximación es de primer grado:

Teniendo en cuenta que la ecuación de la recta que pasa por los puntos (a;f(a))
y (b;f(b)) es:

el área bajo esta línea recta en el intervalo [a;b] está dada por:
 Esta integral constituye una aproximación de la integral de
f(x) en dicho intervalo. El resultado de la integral anterior es:

que se denomina regla del trapecio.

Geométricamente, la regla del trapecio consiste en

aproximar el área debajo de la curva definida por f(x), por el
área bajo la recta que une los puntos (a, f(a)) y (b, f(b)).
Recuerde que la fórmula para calcular el área de un trapecio
es la altura por el promedio de las bases.

ERROR DE LA REGLA DEL TRAPECIO

Cuando se emplea la integral bajo un segmento de línea

recta para aproximar la integral bajo una curva, obviamente se
tiene un error que puede ser importante. Una estimación del
error de truncamiento E de la regla del trapecio es:

donde ξ Є (a;b). La expresión anterior indica que, si la función

a integrar es lineal, la regla del trapecio será exacta. Es decir,
la regla del trapecio tiene grado de precisión n = 1. Además,
sólo es posible aplicar la regla del trapecio si f(x) es de clase
C2[a;b].
LA REGLA COMPUESTA DEL TRAPECIO

Una forma de mejorar la precisión de la regla del trapecio consiste en dividir el

intervalo de integración [a;b] en varios segmentos, y aplicar el método a cada uno
de ellos. Las áreas de los segmentos se suman después para obtener una
aproximación de la integral en todo el intervalo. Las expresiones resultantes se
llaman fórmulas de integración, de aplicación múltiple o compuestas.

La siguiente figura muestra el formato general y la nomenclatura que se usará

para obtener integrales de aplicación múltiple.

Si hay n+1 puntos igualmente espaciados, existen n segmentos del mismo

ancho:

Si a y b se designan como x0 y xn, respectivamente, la integral completa se

representará como:
 Sustituyendo la regla del trapecio en cada integral se obtiene:

o, agrupando términos:

El error que se comete al aplicar la regla compuesta del trapecio está dado
por:

Esto significa que el error es de orden O(h2). Además, cuando las derivadas de f(x)
se conocen, es posible estimar el número de subintervalos necesarios para alcanzar
la precisión deseada.
4.2 APLICACIONES DE LA INTEGRACIÓN NUMÉRICA.
En la actualidad, gracias a la gran evolución que han tenido los métodos numéricos
y su implementación en potentes computadoras, es posible, por ejemplo, modelar
el choque de un vehículo o hacer el análisis aerodinámico estructural de un avión,
resolviendo en cada caso sistemas algebraicos de ecuaciones con varios cientos
de miles (a veces de millones) de incógnitas. Se presentan a continuación algunas
aplicaciones de los métodos numéricos a diversos problemas de ingeniería.

Desde tiempos ancestrales el papel del ingeniero ha sido básicamente el mismo,

tratar de conocer e interpretar los mecanismos de la naturaleza para así poder
modificarla al servicio del hombre. Para ello ha utilizado sus conocimientos,
intuición, experiencia y los medios naturales a los que en cada momento ha tenido
disponibles. Con el gran poder de cómputo que se tiene en estos días, el ingeniero
dispone de grandes ventajas para poder llevar a cabo su misión y abordar cada día
retos más ambiciosos en la solución de nuevos problemas, cuyos aspectos políticos,
económicos, científicos o tecnológicos pueden tener un mayor impacto en la mejora
de la calidad de vida del hombre. Encontramos así aplicaciones de los métodos
numéricos en los ámbitos más diversos desde sectores tecnológicos tan clásicos
como la ingeniería estructural o la aerodinámica de aviones, hasta aplicaciones más
sofisticadas como ingeniería de alimentos, ingeniería médica, diseño de fármacos,
biología, etc.

Mecánica de Sólidos

Existen hoy en día, un gran número de estructuras en ingeniería, que son

modelados desde su concepción utilizando técnicas de elementos finitos. Ejemplos
de ellas puede ser el edificio de la Unidad de Ciencias Económico Administrativas
de la Universidad de Guanajuato (UCEA), el cual se muestra en la Figura 1. Para
modelarlo se utilizaron más de 10,000 elementos de lámina plana, los cuales
representan los materiales que forman la estructura (concreto y acero), en Figura 2
pueden verse los resultados de iso-contornos correspondientes a los momentos
flectores máximos provocados por cargas de servicio estáticas en la estructura. En
la Figura 3 pueden verse los primeros modos de vibración de la estructura, los
cuales se pueden excitar bajo un movimiento sísmico. Dichos modos de vibración
son obtenidos resolviendo un problema de valores y vectores propios generalizados
[*].

Figura1.- Edificio de la Unidad Ciencias Económico-Administrativas de la

Universidad de Guanajuato (UCEA).

Figura 2.- Malla de Elementos Finitos para el modelado del Edificio de la Unidad
Ciencias Económico-Administrativas de la Universidad de Guanajuato (UCEA).
Figura 3.- Diagrama de momentos del edificio, resultado de la modelación numérica
utilizando elementos finitos del Edificio de la Unidad Ciencias Económico-
Administrativas de la Universidad de Guanajuato (UCEA).

En el análisis anterior se puede observar un claro ejemplo de la aplicación de los

métodos numéricos aplicados, utilizando software inteligente capaz de interactuar y
permitir simular procesos complejos, sin necesidad de realizarlos empíricamente,
con lo cual, se inicia el siguiente subtema.
4.3 USO DE HERRAMIENTAS COMPUTACIONALES

Gracias al avance tecnológico de las últimas décadas, en la actualidad podemos

realizar complejos y tardados análisis, en breves periodos y de forma muy práctica,
utilizando software especialmente diseñado para realizar estos complejos trabajos,
tales como se presentan algunos a continuación.

• NASTRAN es un programa de cálculo de elementos finitos.

Los elementos finitos son pequeñas triangulaciones o cuadraturas en las que un

programa como Catia V5 crea las piezas. Cuando la pieza recibe una carga, una
fuerza, un esfuerzo, una presión… la pieza se deforma. Y gracias a esos pequeños
triangulitos podemos saber cómo es esa deformación. Dependiendo del material,
no solo sabremos su deformidad, sino que también sabremos si la pieza se rompe
o no.

Una práctica normal en la

que hacemos cálculos de
este tipo es en las roturas por
fatiga. Una pieza se rompe
cuando repetimos y
repetimos varias veces un
mismo movimiento.

¿Quién no ha roto alguna vez un alambre con la mano?, doblamos y doblamos hasta
que al final el alambre se parte. Bien, eso es rotura por fatiga. Nastran es capaz de
saber cuántas veces debemos doblar el alambre antes de que se parta en dos. Catia
V5 tiene un módulo de cálculo de elementos finitos al cual llamamos MEF o FEM
(en inglés), pero la verdad, no es muy bueno. Por eso las grandes empresas como
Airbus, utilizan Nastran y Patran.
El hermano de Nastran, de la misma compañía, es Patran. Son programas
independientes. Patran sirve para simular todos los cálculos de Patran, simplemente
eso. Una aplicación práctica de Nastran-Patran que nos sirva de ejemplo puede ser
con las piezas diseñadas con Catia V5 para Airbus de materiales compuestos.

Los materiales compuestos son los materiales hechos de fibras (fibra de vidrio, fibra
de carbono, fibra de algodón…).

Normalmente son piezas diseñadas con el módulo Generative Shape Design,

piezas de surperficies. El calculista mediante Nastran-Patran efectúa los cálculos
para saber las zonas conflictivas y se las hace llegar de nuevo al diseñador, que,
mediante el módulo de Composites Design asigna a cada zona de la pieza, los datos
del calculista. A partir de aquí Catia realiza sus cálculos y nos dice cuántas capas
de fibra hay que poner, en que disposición debemos ponerlas y que orientación
deben tener las fibras.

• COMSOL MULTIPHYSICS

Comsol es un software CAE para modelado, análisis y simulación de fenómenos

físicos 3D en ingeniería, como problemas con fluidos, estructurales, térmicos,
electromecánicos entre otros, que permite definir la geometría 3D especificando el
mesh o malleado, cargas y la visualización previa del análisis para luego ejecutar el
postproceso y ver reportes finales.

Tiene interfases con diferentes programas CAD, tiene además una serie de
módulos por aplicación y solvers específicos que se pueden agregar y que a
continuación describimos:

COMSOL AC / DC MODULE

Módulo de ingeniería electrónica que contempla el modelado de capacitores,

inductores, motores y microsensores. Estos elementos también son afectados por
fenómenos térmicos y electromecánicos, las capacidades de AC/DC pueden
ejecutar análisis de estos componentes bajo los fenómenos descritos.

COMSOL MEMS MODULE

Módulo de ingeniería electrónica de COMSOL con capacidades multifísicas para el

diseño de sistemas micro electromecánicos como actuadores, acelerómetros,
sensores, capacitores, etc. (requiere COMSOL Multiphysics).
COMSOL RF MODULE

La ingeniería óptica y las microondas requieren resolver la intensidad de transmitir

de un dispositivo mientras captura efectos de gran magnitud, el módulo RF ofrece
las herramientas para atender estos fenómenos, como resultado puedes modelar
3D antenas, guías de ondas, microondas y componentes ópticos. El módulo provee
un post procesador para la solución de las ondas electromagnéticas.

COMSOL OPTIMIZATION LAB

Es un solver para diseños mecánicos con un amplio rango de aplicaciones de

multifisica, aquí puedes aplicar métodos de optimización a tus diseños incluyendo
optimización de topología y modelado inverso. El programa está basado en el
código SNOPT desarrollado en la Universidad de Stanford y de UC San Diego.

COMSOL ACOUSTICS MODULE

Aplicación de ingeniería mecánica para solucionar los análisis relacionados con

acústica, provee herramientas para modelar la propagación de las ondas en aire,
agua, otros fluidos y sólidos. Su aplicación se incluye en micrófonos bocinas,
ayudas acústicas para humanos, control de ruidos y diseño de edificios.

Octave o GNU Octave es

un programa y lenguaje de
programación para realizar cálculos
numéricos. Como su nombre indica, Octave
es parte del proyecto GNU. Es considerado
el equivalente libre de MATLAB. Entre
varias características que comparten, se
puede destacar que ambos ofrecen
un intérprete, permitiendo ejecutar órdenes
en modo interactivo. Nótese que Octave no
es un sistema de álgebra computacional,
como lo es máxima, sino que está orientado
al análisis numérico.
El proyecto fue creado alrededor del año 1988, pero con una finalidad diferente: ser
utilizado en un curso de diseño de reactores químicos. Posteriormente, en el
año 1992, se decidió extenderlo, y comenzó su desarrollo a cargo de John W. Eaton.
La primera versión alpha fue lanzada el 4 de enero de 1993. Un año más tarde,
el 17 de febrero de 1994, apareció la versión 1.0.

El nombre surge de Octave Levenspiel, profesor de uno de los autores y conocido

por sus buenas aproximaciones, por medio de cálculos elementales, a problemas
numéricos en ingeniería química.

SINTAXIS

o La sintaxis es casi idéntica a la utilizada en MATLAB.

o Es un lenguaje interpretado.

o No permite pasar argumentos por referencia. Siempre son pasados por valor.

o No permite punteros.

o Se pueden generar scripts.

o Soporta gran parte de las funciones de la biblioteca estándar de C.

o Puede ser extendido para ofrecer compatibilidad con las llamadas al

sistema UNIX.

o El lenguaje está pensado para trabajar con matrices, y provee mucha

funcionalidad para trabajar con estas.

o Soporta estructuras similares a los "struct"s de C.

• MATLAB

MATLAB es un entorno de computación y desarrollo de aplicaciones totalmente

integrado orientado para llevar a cabo proyectos en donde se encuentren
implicados elevados cálculos matemáticos y la visualización gráfica de los mismos.
MATLAB integra análisis numérico, cálculo matricial, proceso de señal y
visualización gráfica en un entorno completo donde los problemas y
sus soluciones son expresados del mismo modo en que se escribirían
tradicionalmente, sin necesidad de hacer uso de la programación tradicional.

MATLAB dispone también en la actualidad de un amplio abanico de programas de

apoyo especializados, denominados Toolboxes, que extienden significativamente el
número de funciones incorporadas en el programa principal. Estos Toolboxes
cubren en la actualidad prácticamente casi todas las áreas principales en el mundo
de la ingeniería y la simulación, destacando entre ellos el 'toolbox' de proceso
de imágenes, señal, control robusto, estadística, análisis
financiero, matemáticas simbólicas, redes neurales, lógica difusa, identificación
de sistemas, simulación de sistemas dinámicos, etc. es un entorno de cálculo
técnico, que se ha convertido en estándar de la industria, con capacidades no
superadas en computación y visualización numérica.

De forma coherente y sin ningún tipo de fisuras, integra los requisitos claves de
un sistema de computación técnico: cálculo numérico, gráficos, herramientas para
aplicaciones específicas y capacidad de ejecución en múltiples plataformas.
Esta familia de productos proporciona al estudiante un medio de carácter único,
para resolver los problemas más complejos y difíciles.

Manipulación y Reducción de Datos

MATLAB tiene un rango completo de funciones para pre-procesar datos para
análisis, incluyendo:
• Decimando
• Secciones de datos
• Promediando
• Procesando umbrales
• Filtrando
Numerosas operaciones para manipular arreglos multidimensionales, incluyendo
reticulación e interpolación de datos, están también disponibles.
 • SALOME
Es una plataforma genérica para
pre y post procesado
para simulación numérica. Tiene
una arquitectura abierta y flexible
realizada a partir de
componentes reutilizables. Es
una aplicación disponible
bajo código abierto, licenciada
bajo la LGPL, cuyos ejecutables
y código fuente pueden ser
descargados de su página oficial.

Soporta
o interoperabilidad
entre software de modelado CAD
y de computación software) (enlace CAD-CAE).

o Simplifica la integración de nuevos componentes en sistemas heterogéneos

de análisis numérico.

o Proporciona una interfaz genérica, amigable y eficiente, que permite reducir

los costes y retrasos of carrying out the studies.

o Reduce el tiempo de aprendizaje al tiempo específico necesario para

aprender el software basado en esta plataforma.

o Todas las funcionalidades están accesibles a través de la

consola Python integrada.

Actualmente, el software está disponible para las distribuciones Debian

(Etch/Lenny/Squeeze), Mandriva (2008/2010, 32/64 bits), Red Hat Enterprise 4 (64
bits), Scientific Linux 5.1 (64 bits) y CentOS 5.5 (64 bits). Sin embargo, el código
fuente está disponible para cualquier programador que desee desarrollarlo y ponerlo
a disposición de otras plataformas, entre ellas Microsoft Windows.
BIBLIOGRAFÍA

[Link]

[Link]

[Link]

[Link]