Metodo Romberg

Método de Romberg
Ir a la navegaciónIr a la búsqueda
En análisis numérico, el Método de Romberg genera una matriz triangular cuyos
elementos son estimaciones numéricas de la integral definida siguiente:
{\displaystyle \int _{a}^{b}f(x)\,dx}\int _{a}^{b}f(x)\,dx

usando la extrapolación de Richardson de forma reiterada en la regla del trapecio. El
método de Romberg evalúa el integrando en puntos equiespaciados del intervalo de
integración estudiado. Para que este método funcione, el integrando debe ser
suficientemente derivable en el intervalo, aunque se obtienen resultados bastante buenos
incluso para integrandos poco derivables. Aunque es posible evaluar el integrando en
puntos no equiespaciados, en ese caso otros métodos como la cuadratura gaussiana o la
cuadratura de Clenshaw–Curtis son más adecuados.
Índice
1 Método
2 Implementación en Python del método de Romberg
3 Ejemplo
4 Referencias
5 Enlaces externos
Método
El método se define de forma recursiva así:
{\displaystyle R(0,0)={\frac {1}{2}}(b-a)(f(a)+f(b))}{\displaystyle R(0,0)={\frac {1}

{2}}(b-a)(f(a)+f(b))}
{\displaystyle R(n,0)={\frac {1}{2}}R(n-1,0)+h_{n}\sum _{k=1}^{2^{n-1}}f(a+(2k-
1)h_{n})}{\displaystyle R(n,0)={\frac {1}{2}}R(n-1,0)+h_{n}\sum _{k=1}^{2^{n-
1}}f(a+(2k-1)h_{n})}
{\displaystyle R(n,m)=R(n,m-1)+{\frac {1}{4^{m}-1}}(R(n,m-1)-R(n-1,m-1))}
{\displaystyle R(n,m)=R(n,m-1)+{\frac {1}{4^{m}-1}}(R(n,m-1)-R(n-1,m-1))}
o
{\displaystyle R(n,m)={\frac {1}{4^{m}-1}}(4^{m}R(n,m-1)-R(n-1,m-1))}

{\displaystyle R(n,m)={\frac {1}{4^{m}-1}}(4^{m}R(n,m-1)-R(n-1,m-1))}
donde
{\displaystyle n\geq 1}{\displaystyle n\geq 1}

{\displaystyle m\geq 1}{\displaystyle m\geq 1}
{\displaystyle h_{n}={\frac {b-a}{2^{n}}}.}{\displaystyle h_{n}={\frac {b-a}
{2^{n}}}.}
La cota superior asintótica del error de R(n,m) es:
{\displaystyle O\left(h_{n}^{2^{m+1}}\right).}{\displaystyle
O\left(h_{n}^{2^{m+1}}\right).}
La extrapolación a orden cero {\displaystyle R(n,0)}{\displaystyle R(n,0)} es
equivalente a la Regla del trapecio con {\displaystyle n+2}n+2 puntos. a orden uno
{\displaystyle R(n,1)}{\displaystyle R(n,1)} es equivalente a la Regla de Simpson con
{\displaystyle n+2}n+2 puntos.
Cuando la evaluación del integrando es numéricamente costosa, es preferible

reemplazar la interpolación polinómica de Richardson por la interpolación racional
propuesta por Bulirsch & Stoer.
Implementación en Python del método de Romberg

La siguiente es una implementación del método de Romberg en Python:
def print_row(lst):
print ' '.join('%11.8f' % x for x in lst)
def romberg(f, a, b, eps = 1E-8):

"""Approximate the definite integral of f from a to b by Romberg's method.
eps is the desired accuracy."""
R = [[0.5 * (b - a) * (f(a) + f(b))]] # R[0][0]
print_row(R[0])
n=1
while True:
h = float(b-a)/2**n
R.append((n+1)*[None]) # Add an empty row.
R[n][0] = 0.5*R[n-1][0] + h*sum(f(a+(2*k-1)*h) for k in range(1, 2**(n-1)+1)) #
for proper limits
for m in range(1, n+1):
R[n][m] = R[n][m-1] + (R[n][m-1] - R[n-1][m-1]) / (4**m - 1)
print_row(R[n])
if abs(R[n][n-1] - R[n][n]) < eps:
return R[n][n]
n += 1
from math import *
# In this example, the error function erf(1) is evaluated.

print romberg(lambda t: 2/sqrt(pi)*exp(-t*t), 0, 1)
Ejemplo
Como ejemplo, se le integra la función gaussiana en el intervalo {\displaystyle [O,1]}
{\displaystyle [O,1]}, esto es la función error evaluada en 1, cuyo valor es {\displaystyle
{\rm {erf}}(1)\doteq 0.842700792949715}{\displaystyle {\rm {erf}}(1)\doteq
0.842700792949715}. Se calculan los elementos de la matriz triangular fila a fila,
terminando los cálculos cuando la diferencia entre las dos últimas filas es menor que
1E-8.
0.77174333
0.82526296 0.84310283
0.83836778 0.84273605 0.84271160
0.84161922 0.84270304 0.84270083 0.84270066
0.84243051 0.84270093 0.84270079 0.84270079 0.84270079
El resultado en la esquina inferior derecha de la matriz triangular es el resultado
correcto con la precisión deseada. Nótese que este resultado se deriva de
aproximaciones mucho peores obtenidas con la regla del trapecio mostradas aquí en la
primera columna de la matriz triangular.
Referencias
Richardson, L. F. (1911), «The Approximate Arithmetical Solution by Finite
Differences of Physical Problems Involving Differential Equations, with an Application
to the Stresses in a Masonry Dam», Philosophical Transactions of the Royal Society of
London. Series A 210: pp. 307-357.
Romberg, W. (1955), «Vereinfachte numerische Integration», Norske Videnskabers
Selskab Forhandlinger (Trondheim) 28 (7): pp. 30-36.
Thacher, Jr., Henry C. (julio de 1964), «Remark on Algorithm 60: Romberg
integration», Communications of the ACM 7 (7): 420-421.
Bauer, F.L.; Rutishauser, H.; Stiefel, E. (1963), «New aspects in numerical quadrature»,
en Metropolis, N. C., et al., ed., Experimental Arithmetic, high-speed computing and
mathematics, Proceedings of Symposia in Applied Mathematics (AMS) (15): pp. 199-
218.
Bulirsch, Roland; Stoer, Josef (1967), «Handbook Series Numerical Integration.
Numerical quadrature by extrapolation», Numerische Mathematik 9: 271-278.
Mysovskikh, I.P. (2002), «Romberg method», en Hazewinkel, Michiel, ed.,
Encyclopaedia of Mathematics, Springer-Verlag, ISBN 1-4020-0609-8.
Enlaces externos
ROMBINT -- código para MATLAB (autor: Martin Kacenak)
Module for Romberg Integration
Control de autoridades
Proyectos WikimediaWd Datos: Q1078449IdentificadoresMicrosoft Academic:
126976006
Categoría: Integración numérica
Menú de navegación
No has accedido
Discusión
Contribuciones
Crear una cuenta
Acceder
ArtículoDiscusión
LeerEditarVer historialBuscar
Buscar en Wikipedia
Portada
Portal de la comunidad
Actualidad
Cambios recientes
Páginas nuevas
Página aleatoria
Ayuda
Donaciones
Notificar un error
Herramientas
Lo que enlaza aquí
Cambios en enlazadas
Subir archivo
Deformación
(Redirigido desde «Deformacion plastica»)
Ir a la navegaciónIr a la búsqueda
La deformación es el cambio en el tamaño o forma de un cuerpo debido a tensiones
internas producidos por una o más fuerzas aplicadas sobre el mismo , la ocurrencia de
dilataciones térmicas o la aparición de efectos viscoelásticos.
Las deformaciones puede ser termodinámicamente reversibles o irreversibles. La

deformación elástica es siempre reversible, mientras que la deformación plástica es
termodinámicamente irreversible.
Índice
1 Medidas de la deformación
1.1 Deformación unidimensional
1.2 Deformación de un cuerpo
2 Deformaciones elástica y plástica
3 Conceptos relacionados
3.1 Desplazamientos
3.2 Energía de deformación
4 Véase también
5 Referencias
5.1 Bibliografía
5.2 Enlaces externos
Medidas de la deformación
Deformación unidimensional
La magnitud más simple para medir la deformación es lo que en ingeniería se llama
deformación axial o deformación unitaria. Se define como el cambio de longitud por
unidad de longitud:
(*)de la misma magnitud {\displaystyle e\ ={\frac {\Delta s}{s}}={\frac {s'-s}{s}}}

{\displaystyle e\ ={\frac {\Delta s}{s}}={\frac {s'-s}{s}}}
Donde {\displaystyle s}s es la longitud inicial de la zona en estudio y {\displaystyle s'}

{\displaystyle s'} la longitud final o deformada. Es útil para expresar los cambios de
longitud de un cable o un prisma mecánico. La deformación calculada de acuerdo a (*)
se llama deformación ingenieril. En la práctica se pueden usar otras medidas
relacionadas con estas como el estiramiento:
{\displaystyle \lambda ={\frac {s'}{s}}=1+e}{\displaystyle \lambda ={\frac {s'}

{s}}=1+e}
La deformación axial logarítmica o deformación de Hencky que se define como:
{\displaystyle \varepsilon _{H}=\ln(\lambda )=\ln(1+e)}{\displaystyle \varepsilon

_{H}=\ln(\lambda )=\ln(1+e)}
La deformación de Green-Lagrange viene dada por:

{\displaystyle \varepsilon _{G}={\frac {1}{2}}(\lambda ^{2}-1)={\frac {(1+e)^{2}-1}
{2}}=e+{\frac {e^{2}}{2}}={\frac {\exp(2\varepsilon _{H})-1}{2}}}{\displaystyle
\varepsilon _{G}={\frac {1}{2}}(\lambda ^{2}-1)={\frac {(1+e)^{2}-1}{2}}=e+{\frac
{e^{2}}{2}}={\frac {\exp(2\varepsilon _{H})-1}{2}}}
La deformación de Euler-Almansi viene dada por:
{\displaystyle \varepsilon _{E}={\frac {1}{2}}\left(1-{\frac {1}{\lambda

^{2}}}\right)={\frac {1+e/2}{1+e}}e={\frac {1-\exp(-2\varepsilon _{H})}{2}}}
{\displaystyle \varepsilon _{E}={\frac {1}{2}}\left(1-{\frac {1}{\lambda
^{2}}}\right)={\frac {1+e/2}{1+e}}e={\frac {1-\exp(-2\varepsilon _{H})}{2}}}
Deformación de un cuerpo
En la mecánica de sólidos deformables la deformación puede tener lugar según diversos
modos y en diversas direcciones, y puede además provocar distorsiones en la forma del
cuerpo, en esas condiciones la deformación de un cuerpo se puede caracterizar por un
tensor (más exactamente un campo tensorial) de la forma:
{\displaystyle [D]={\begin{pmatrix}\varepsilon _{11}&\varepsilon _{12}&\varepsilon

_{13}\\\varepsilon _{21}&\varepsilon _{22}&\varepsilon _{23}\\\varepsilon
_{31}&\varepsilon _{32}&\varepsilon _{33}\end{pmatrix}}}{\displaystyle
[D]={\begin{pmatrix}\varepsilon _{11}&\varepsilon _{12}&\varepsilon
_{13}\\\varepsilon _{21}&\varepsilon _{22}&\varepsilon _{23}\\\varepsilon
_{31}&\varepsilon _{32}&\varepsilon _{33}\end{pmatrix}}}
Donde cada una de las componentes de la matriz anterior, llamada tensor deformación
representa una función definida sobre las coordenadas del cuerpo que se obtiene como
combinación de derivadas del campo de desplazamientos de los puntos del cuerpo.
Deformaciones elástica y plástica

Tanto para la deformación unitaria como para el tensor deformación se puede
descomponer el valor de la deformación en:
Deformación plástica, irreversible o permanente. Modo de deformación en que el

material no regresa a su forma original después de retirar la carga aplicada. Esto sucede
porque, en la deformación plástica, el material experimenta cambios termodinámicos
irreversibles al adquirir mayor energía potencial elástica. La deformación plástica es lo
contrario a la deformación reversible.
Deformación elástica, reversible o no permanente, el cuerpo recupera su forma original
al retirar la fuerza que le provoca la deformación. En este tipo de deformación, el sólido,
al variar su estado tensional y aumentar su energía interna en forma de energía potencial
elástica, solo pasa por cambios termodinámicos reversibles.
Comúnmente se entiende por materiales elásticos, aquellos que sufren grandes
elongaciones cuando se les aplica una fuerza, como la goma elástica que puede estirarse
sin dificultad recuperando su longitud original una vez que desaparece la carga. Este
comportamiento, sin embargo, no es exclusivo de estos materiales, de modo que los
metales y aleaciones de aplicación técnica, piedras, hormigones y maderas empleados
en construcción y, en general, cualquier material, presenta este comportamiento hasta un
cierto valor de la fuerza aplicada; si bien en los casos apuntados las deformaciones son
pequeñas, al retirar la carga desaparecen.
Al valor máximo de la fuerza aplicada sobre un objeto para que su deformación sea
elástica se le denomina límite elástico y es de gran importancia en el diseño mecánico,
ya que en la mayoría de aplicaciones es este y no el de la rotura, el que se adopta como
variable de diseño (particularmente en mecanismos). Una vez superado el límite elástico
aparecen deformaciones plásticas (que son permanentes tras retirar la carga)
comprometiendo la funcionalidad de ciertos elementos mecánicos.
Conceptos relacionados
Desplazamientos
Cuando un medio continuo se deforma, la posición de sus partículas materiales cambia
de ubicación en el espacio. Este cambio de posición se representa por el llamado vector
desplazamiento, u = (ux, uy, uz). No debe confundirse desplazamiento con
deformación, porque son conceptos diferentes aunque guardan una relación matemática
entre ellos:
{\displaystyle \varepsilon _{ij}={1 \over 2}\left({\partial u_{i} \over \partial x_{j}}+

{\partial u_{j} \over \partial x_{i}}+\sum _{k}{\partial u_{k} \over \partial x_{i}}
{\partial u_{k} \over \partial x_{j}}\right)\approx {1 \over 2}\left({\partial u_{i}
\over \partial x_{j}}+{\partial u_{j} \over \partial x_{i}}\right)}{\displaystyle
\varepsilon _{ij}={1 \over 2}\left({\partial u_{i} \over \partial x_{j}}+{\partial u_{j}
\over \partial x_{i}}+\sum _{k}{\partial u_{k} \over \partial x_{i}}{\partial u_{k}
\over \partial x_{j}}\right)\approx {1 \over 2}\left({\partial u_{i} \over \partial x_{j}}+
{\partial u_{j} \over \partial x_{i}}\right)}
donde {\displaystyle \varepsilon _{ij}}{\displaystyle \varepsilon _{ij}} son las

componentes del tensor deformación de Green-Lagrange. La diferencia entre
deformación y desplazamiento puede ilustrarse considerando un voladizo o ménsula
empotrada en un extremo y libre en el otro, las deformaciones son máximas en el
extremo empotrado y cero en el extremo libre, mientras que los desplazamientos son
cero en el extremo empotrado y máximos en el extremo libre.
Energía de deformación
Artículo principal: Energía de deformación
La deformación es un proceso termodinámico en el que la energía interna del cuerpo
acumula energía potencial elástica. A partir de unos ciertos valores de la deformación se
pueden producir transformaciones del material y parte de la energía se disipa en forma
de plastificado, endurecimiento, fractura o fatiga del material.
Véase también
Esfuerzo interno
Elasticidad (mecánica de sólidos)
Ley de Hooke
Ensayo de tracción
Galga extensiométrica
Referencias
Bibliografía
Luis Ortiz Berrocal (2007). Resistencia de materiales, Madrid: Ed. McGraw-Hill. ISBN
9788448156336.
Enlaces externos
Wikcionario tiene definiciones y otra información sobre deformación.
Página donde se enumeran aspectos relevantes en el trabajo de metales por deformación
plástica.
Control de autoridades
Proyectos WikimediaWd Datos: Q193514Commonscat Multimedia: Deformation
IdentificadoresLCCN: sh85036465Microsoft Academic: 46585869OntologíasNúmero
IEV: 113-03-57
Categorías: Magnitudes físicasMecánica de sólidos deformables
Menú de navegación
No has accedido
Discusión
Contribuciones
Crear una cuenta
Acceder
ArtículoDiscusión
LeerEditarVer historialBuscar
Buscar en Wikipedia
Portada
Portal de la comunidad
Actualidad
Cambios recientes
Páginas nuevas
Página aleatoria
Ayuda
Donaciones
Notificar un error
Herramientas
Lo que enlaza aquí
Cambios en enlazadas
Subir archivo
Páginas especiales
Enlace permanente
Información de la página
Citar esta página
Elemento de Wikidata
Imprimir/exportar
Crear un libro
Descargar como PDF
Versión para imprimir
En otros proyectos
Wikimedia Commons
En otros idiomas
‫العربية‬
Deutsch
English
Français
हिन्दी
Bahasa Indonesia
Português
Русский
中文
42 más
Editar enlaces
Esta página se editó por última vez el 19 may 2021 a las 16:00.
El texto está disponible bajo la Licencia Creative Commons Atribución Compartir Igual
3.0; pueden aplicarse cláusulas adicionales. Al usar este sitio, usted acepta nuestros
términos de uso y nuestra política de privacidad.
Wikipedia® es una marca registrada de la Fundación Wikimedia, Inc., una organización
sin ánimo de lucro.
Política de privacidadAcerca de WikipediaLimitación de responsabilidadVersión para
móvilesDesarrolladoresEstadísticasDeclaración de cookiesWikimedia
FoundationPowered by MediaWiki

Metodo Romberg

Cargado por

Copyright:

Formatos disponibles

Metodo Romberg

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Metodo Romberg

Cargado por

Copyright:

Formatos disponibles

Método de Romberg

{\displaystyle \int _{a}^{b}f(x)\,dx}\int _{a}^{b}f(x)\,dx

{\displaystyle R(0,0)={\frac {1}{2}}(b-a)(f(a)+f(b))}{\displaystyle R(0,0)={\frac {1}

{\displaystyle R(n,m)={\frac {1}{4^{m}-1}}(4^{m}R(n,m-1)-R(n-1,m-1))}

{\displaystyle n\geq 1}{\displaystyle n\geq 1}

Cuando la evaluación del integrando es numéricamente costosa, es preferible

Implementación en Python del método de Romberg

def romberg(f, a, b, eps = 1E-8):

from math import *

# In this example, the error function erf(1) is evaluated.

Las deformaciones puede ser termodinámicamente reversibles o irreversibles. La

(*)de la misma magnitud {\displaystyle e\ ={\frac {\Delta s}{s}}={\frac {s'-s}{s}}}

Donde {\displaystyle s}s es la longitud inicial de la zona en estudio y {\displaystyle s'}

{\displaystyle \lambda ={\frac {s'}{s}}=1+e}{\displaystyle \lambda ={\frac {s'}

La deformación axial logarítmica o deformación de Hencky que se define como:

{\displaystyle \varepsilon _{H}=\ln(\lambda )=\ln(1+e)}{\displaystyle \varepsilon

La deformación de Green-Lagrange viene dada por:

La deformación de Euler-Almansi viene dada por:

{\displaystyle \varepsilon _{E}={\frac {1}{2}}\left(1-{\frac {1}{\lambda

{\displaystyle [D]={\begin{pmatrix}\varepsilon _{11}&\varepsilon _{12}&\varepsilon

Deformaciones elástica y plástica

Deformación plástica, irreversible o permanente. Modo de deformación en que el

{\displaystyle \varepsilon _{ij}={1 \over 2}\left({\partial u_{i} \over \partial x_{j}}+

donde {\displaystyle \varepsilon _{ij}}{\displaystyle \varepsilon _{ij}} son las

También podría gustarte

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.