El esquema de Newton-Raphson para problemas unidimensionales
Dr. Alejo O. Sfriso Universidad de Buenos Aires SRK Consulting (Argentina) AOSA
materias.fi.uba.ar/6408 latam.srk.com www.aosa.com.ar
[email protected] [email protected] [email protected]
El esquema de Newton-Rapshon
El esquema de Newton-Raphson
Sea el problema no lineal Se parte de un estado convergido Se impone un incremento de def. El problema es hallar 𝑝'() que cumpla
𝑝 = 𝐾 𝑝 · 𝜖& 𝑝' = 𝐾 𝑝' · 𝜖&' 𝜖&'() = 𝜖&' + Δ𝜖& 𝑝'() = 𝐾 𝑝'() · 𝜖&'()
Se plantea la forma implícita
𝐹 𝑝 = 𝑝 − 𝐾 𝑝 · 𝜖&'() = 0
Por teorema de Taylor
𝐹 𝑝/() = 𝐹 𝑝/ +
De donde se despeja
Δ𝑝 =
41 25 01⁄02
01 3 02 /
Δ𝑝 = 0
5
Se actualiza 𝑝/() = 𝑝/ + Δ𝑝 y se itera hasta convergencia 2
1
El esquema de Newton-Rapshon
Cómo funciona
3
El esquema de Newton-Rapshon
Cómo se arregla (bisección)
• Se eligen topes ridículamente lejanos
pmin
pn
pmax
4
2
El esquema de Newton-Rapshon
Cómo se arregla (bisección)
• Se eligen topes ridículamente lejanos • Se resuelve N-R
pmin
p( ) 1
pn
pmax
5
El esquema de Newton-Rapshon
Cómo se arregla (bisección)
• Se eligen topes ridículamente lejanos • Se resuelve N-R • Si la presión cae fuera de los topes, se ignora el resultado y se adopta un promedio de los topes
pmin
pn
p( ) 1
pmax
p( ) 1
6
3
El esquema de Newton-Rapshon
Cómo se arregla (bisección)
• Se eligen topes ridículamente lejanos • Se resuelve N-R • Si la presión cae fuera de los topes, se ignora el resultado y se adopta un promedio de los topes • En función del valor de F se corrigen los topes pmin
p( ) = pmax 1
7
El esquema de Newton-Rapshon
Cómo se arregla (bisección)
El reducido radio de convergencia de NR originó el desarrollo de decenas de métodos de integración La combinación N-R con bisección es muy robusta, aunque no es la más eficiente en la mayoría de los casos Sin bisección no converge aún para un p0 razonable
Con bisección converge en pocos pasos aunque p0 sea extremo
8
4