Trabajo Práctico N° 2
Carrera: Ingeniería Industrial.
Segundo Cuatrimestre 2024.
ALUMNOS:
1 - B-5847/5 Borgatello, Juan Cruz jborgatello98@[Link]
2 - G-5223/9 Ghella, Matías Nicolás [Link]@[Link]
3 - L- 3002/3 Lares, Brian brianlares10@[Link]
4 - R-4386/9 Raies Mendi, Tadeo tadeoraiesmendi@[Link]
INTRODUCCIÓN:
En este trabajo nos centraremos en estudiar el endurecimiento de los materiales,
poniendo en práctica distintos métodos para encontrar raíces y así sacar
conclusiones sobre los resultados.
El endurecimiento es el fenómeno mediante el cual un material aumenta su
resistencia a medida que se deforma plásticamente. En este trabajo se propone la
resolución de una ecuación no lineal derivada de un modelo de endurecimiento
plástico.
1) Conversión y Derivación
𝐸*ϵ𝑝 𝑁
Sea la ecuación de endurecimiento: σ𝑓 = σ𝑦(1 + σ𝑦
)
Donde: Tensión de fluencia, σ𝑓= 200 MPa
Límite elástico, σ𝑦= 150 MPa
Módulo de young, E= 210 GPa
exponente de endurecimiento, N=0.15
𝐸*ϵ𝑝 𝑁
Podemos definirla como una función: 𝑓(ϵ𝑝) = σ − σ𝑦(1 + σ𝑦
)
𝑓
Sobre la cual nos interesa conocer sus raíces y el comportamiento de f’ y f’’ en las
proximidades de la misma.
Figura 1: Gráficos de f,f' y f’’.
f, f’ y f’’ se graficaron en el intervalo (0;0.005) y como primera aproximación
podemos decir que la raíz de f se encuentra muy próxima a ϵ𝑝= 0.004.
𝑓(ϵ𝑝+ℎ)−𝑓(ϵ𝑝−ℎ)
Derivada primera numérica: 𝑓'(ϵ𝑝) ≈ 2ℎ
𝐸𝑁 𝐸*ϵ𝑝 𝑁−1
Derivada primera analitica: 𝑓'(ϵ𝑝) =− σ𝑦
(1 + σ𝑦
)
𝑓(ϵ𝑝+ℎ)−2𝑓(ϵ𝑝)+𝑓(ϵ𝑝−ℎ)
Derivada segunda numérica: 𝑓''(ϵ𝑝) ≈ 2
ℎ
2
−𝐸 𝑁(𝑁−1) 𝐸*ϵ𝑝 𝑁−2
Derivada segunda analitica: 𝑓''(ϵ𝑝) = 2 (1 + σ𝑦
)
σ𝑦
El paso de derivación propuesto, h=0.0001. Con valores mayores la curva no se
llegaba a apreciar, parecía más un polígono.
2) Método de Newton
En este apartado se pide obtener el valor de la deformación plástica acumulada que
satisfaga la ecuación de endurecimiento con los valores dados utilizando el método
de Newton.
Para comenzar la resolución debemos verificar la hipótesis adicional del método
Newton-Raphson, la cual dice “La derivada segunda de la función no se anula en un
intervalo [a,b] que contenga a la raíz”.
Para realizar dicha verificación usaremos el gráfico de la derivada segunda de la
función obtenida en el apartado anterior:
Figura 2: Gráfico de f’’
Como podemos observar se cumple la hipótesis debido a que nunca toma el valor
cero. Por lo tanto continuamos con el análisis.
Para calcular el punto inicial determinamos que si f(a) y f´´(a) tienen igual signo el
punto de arranque adecuado es “a” y en caso contrario, es “b”. En nuestro caso
obtuvimos que el punto inicial es 0.
Analizaremos la siguiente ecuación:
𝐸*ϵ𝑝 𝑁
𝑓(ϵ𝑝) = σ − σ𝑦(1 + σ𝑦
)
𝑓
210000 𝑀𝑝𝑎*ϵ𝑝 0.15
𝑓(ϵ𝑝) = 200 𝑀𝑝𝑎 − 150 𝑀𝑝𝑎(1 + 150 𝑀𝑝𝑎
)
Una vez aplicado el método a la función, obtuvimos una raíz en Ep=0.004148 en 5
iteraciones. Para calcular la constante Asintótica A necesitamos el valor real de la
raíz, el cual se obtuvo haciendo más precisa la tolerancia y aumentando el número
de iteraciones, obteniendo Epreal=0.0041475.
Como se puede observar, el valor de A va disminuyendo con cada iteración, algo
que es coherente debido a que va haciendo “pasos” más cortos a medida que se va
acercando más a la raíz.
3) Método de Punto Fijo
Para encontrar la raíz con el método de punto fijo primero definimos la siguiente
1
σ𝑦 σ𝑓 𝑁
1
generatriz: 𝑔(ϵ𝑝) = ( 𝐸
(( σ ) − 1) + ϵ𝑝) 2
𝑦
Que cumple con: 𝑔(ϵ𝑝) = ϵ𝑝
Aplicando el método a g en el intervalo (0;0.005) con 0 como el punto de arranque y
−4
una tolerancia de 10 la raíz estimada es ϵ𝑝= 0.004083.
En cuanto a la cantidad máxima de iteraciones despejamos n de la siguiente
𝑛
|
𝐾 𝑃1−𝑃0 |
|
desigualdad: 𝑃 − 𝑃𝑛 ≤ | 1−𝐾
−4
10
𝑙𝑜𝑔((1−𝐾)( −3 ))
Quedando: 𝑛 ≥ 𝑙𝑜𝑔(𝐾)
2*10
Nos da una cantidad de iteraciones necesarias de
n=6.
Figura 3: G(x) vs Y=x
4) Método de Bisección
En este punto obtendremos el valor del límite elástico σ𝑦 que satisfaga la ecuación
de endurecimiento con el siguiente despeje, se utilizó el Ep encontrado en el
apartado anterior:
𝐸*ϵ𝑝 𝑁
𝑓(σ𝑦) = σ − σ𝑦(1 + σ𝑦
)
𝑓
210000 𝑀𝑝𝑎*0.004083 0.2
𝑓(σ𝑦) = 200 𝑀𝑝𝑎 − σ𝑦(1 + σ𝑦
)
Para determinar los puntos a y b probamos la función con distintos valores de σ𝑦 ,
los cuales debían dar como resultado valores de distinto signo para saber que entre
ellos se encuentra la raíz. Establecimos a=130, b=140 y una tolerancia de 10-3
Como resultado se obtuvo una valor de σ𝑦=134.035645, en 12 iteraciones:
Quiere decir que con dicho valor de σ𝑦 Se cumple la ecuación de endurecimiento.
5) Análisis de Sensibilidad
Para concluir el trabajo se realizará un análisis de sensibilidad graficando la curva
del factor σ𝑦 en función de N utilizando el método de falsa posición. Para ello
analizamos el siguiente despeje en función de σ𝑦 variando N de 0 a 1 y σ𝑦 de 2 a
200:
𝐸*ϵ𝑝 𝑁
𝑓(σ𝑦) = σ − σ𝑦(1 + σ𝑦
)
𝑓
210000 𝑀𝑝𝑎*0.02 𝑁
𝑓(σ𝑦) = 200 𝑀𝑝𝑎 − σ𝑦(1 + σ𝑦
)
Buscamos los valores de σ𝑦 tal que se cumpla :
210000 𝑀𝑝𝑎*0.02 𝑁
200 𝑀𝑝𝑎 − σ𝑦(1 + σ𝑦
) = 0 , para cada valor del vector definido como
N=0:1/100:1.
Figura 4:Gráfico de Ty(n)
Viendo el comportamiento de la función Ty(n) podemos afirmar que el método es
estable entre 0 ≤ n ≤ 0.6 ya que pequeños cambios en n producen pequeños
cambios en Ty, es decir, la curva tiene un comportamiento suave.
Para valores mayores a N=0.6 se obtienen resultados imaginarios o negativos algo
que carece de sentido teniendo en cuenta que estamos analizando el límite elástico.
Concluimos que a medida que aumenta el coeficiente de endurecimiento N
disminuye el límite elástico, por lo que con un N mayor, será más sencillo deformar
plásticamente el material.