TP M.E.F.
Esteban Manzolillo
Matricula: 44.695.922
Caso de análisis: 8; datos C)
1) Descripción del Problema y Resultados Esperados
El problema consiste en el análisis tensional de una planchuela metálica
con un orificio circular central, sometida a un momento flector puro en el
plano. El objetivo principal es determinar el factor de concentración de
tensiones alrededor del orificio, definido como la relación entre la tensión
máxima localizada cerca del agujero y la tensión nominal o de referencia
en zonas alejadas del mismo.
Se espera obtener como resultados principales:
Valores precisos de tensión máxima alrededor del orificio.
Valores de tensión nominal lejos de la zona del orificio.
Factor de concentración de tensiones para distintas configuraciones
geométricas.
Comparación y validación con soluciones analíticas disponibles en la
literatura.
Datos Geométricos y Materiales (Caso C)
Los datos geométricos específicos del caso C son:
Altura de la planchuela (D): 0.14 m.
Espesor de la planchuela (t): 0.10 m.
Radio del orificio (r): 0.025 m.
Momento flector aplicado (M): 1.4 kN·m.
Se eligió una relación de dimensiones longitud-altura (L/D): 2. Por lo
tanto, L = 2 × 0.14 m = 0.28 m.
Para la caracterización completa del problema, se utilizarán valores
estándar provenientes de bibliografía:
Módulo de elasticidad del material (E): 210 GPa (acero estructural).
Coeficiente de Poisson del material (ν): 0.3 (acero estructural).
Estos valores se justifican en base a propiedades estándar reportadas en
bibliografía técnica para aceros estructurales convencionales (como el
ASTM A36 o equivalente).
2) Descripción y Justificación del Modelo de Elementos Finitos
El modelo de elementos finitos propuesto aprovechará las simetrías
existentes debido al diseño geométrico y la carga aplicada. La planchuela
y su carga presentan doble simetría (vertical y horizontal), por lo cual se
discretizará únicamente un cuarto de la geometría original. Esto optimiza
considerablemente el esfuerzo computacional y mantiene la precisión del
análisis.
Dominio elegido:
Se discretiza un cuarto del dominio original aprovechando simetría
horizontal y vertical.
Condiciones de borde y simetría:
Borde inferior: restricción del desplazamiento vertical (Uy = 0).
==<< mal, es en x
Borde izquierdo: restricción del desplazamiento horizontal (Ux = 0).
Carga aplicada:
Considerando los datos específicos del caso C y la relación L/D = 2 elegida
anteriormente, el cálculo de la carga aplicada es el siguiente:
q = (4 × M) / (D × L)
Donde:
M = 1.4 kN·m = 1400 N·m
D = 0.14 m
L = 0.28 m (ya que L/D = 2)
Entonces:
q = (4 × 1400 N·m) / (0.14 m × 0.28 m)
q = 5600 N·m / 0.0392 m²
q = 142857.14 N/m = 142.86 kN/m
Esta carga distribuida horizontal se aplicará en el borde superior ajustado
a la simetría del cuarto del modelo.
Tipo de elemento utilizado:
Elementos bidimensionales triangulares de 6 nodos con
interpolación cuadrática para análisis tensional en régimen elástico
lineal, adecuados para modelar placas delgadas sometidas a cargas
en su plano.
Este modelo permitirá evaluar con precisión las distribuciones de tensión
alrededor del orificio, justificando y verificando posteriormente los valores
obtenidos frente a soluciones analíticas clásicas disponibles en la
literatura técnica.
3) Discretizaciones utilizadas:
Se utilizo en primera instancia una malla estructurada de 0.05m de
tamaño nominal, y luego de realizarse el calculo se obtuvo lo siguiente:
Se observó una concentración de tensiones en la esquina inferior derecha,
posiblemente debido a la restricción colocada. La mayor concentración se
encuentra en la parte superior del orificio, con una tensión mínima
(máxima en valor absoluto) de -1816.3 kPa. Para obtener resultados más
precisos, se refinó la malla en esa línea.
Comparación entre resultados numéricos y analíticos
Para validar los resultados obtenidos con el modelo de elementos finitos,
se comparan las tensiones máximas obtenidas numéricamente con la
solución analítica disponible en la literatura, según el manual de Roark y
Peterson.
Solución Analítica
Según la figura del manual consultado, el estado de tensiones en una
planchuela con orificio circular centrado sometida a flexión pura presenta
una tensión nominal definida como:
q = (4 × M) / (D × L)
Donde: M = 1.4 kN·m = 1400 N·m
D = 0.14 m
L = 0.28 m (ya que L/D = 2)
Entonces: q = (4 × 1400 N·m) / (0.14 m × 0.28 m)
q = 5600 N·m / 0.0392 m²
q = 142857.14 N/m = 142.86 kN/m
Resultado Numérico (GiD)
En la simulación de GiD con refinamiento local alrededor del orificio, la
tensión máxima obtenida (ver escala de colores en la segunda malla
refinada) es aproximadamente: Sigma max = 2483.6 kPa (aprox.)
Cálculo del factor de concentración de tensiones (Kt):
Kt = 3.61
Análisis de los resultados
Teóricamente, el manual indica que Kt=2 para flexión con orificio circular,
independientemente de la relación r/D. Sin embargo, esta solución se
refiere a una ubicación específica (puntos X) que no necesariamente
coincide con la que capta el software si no hay suficiente refinamiento o si
el criterio de evaluación no corresponde exactamente.
Numéricamente, el modelo con refinamiento más preciso muestra un Kt \
approx 3.6, lo cual puede explicarse por la presencia de singularidades
numéricas en bordes o vértices agudos, errores de interpolación en zonas
de alta curvatura, o diferencias entre la zona donde se evalúa sigma máx.
y la zona donde la solución analítica lo predice.
Conclusión:
El modelo numérico reproduce el comportamiento cualitativo de la
concentración de tensiones. La cuantificación del Kt está por encima del
valor teórico, lo cual puede deberse a discretización insuficiente o errores
localizados cerca del borde del agujero. Se recomienda realizar un análisis
de convergencia con sucesivas mallas refinadas en la zona del orificio
para verificar si Kt tiende hacia el valor analítico.
Los resultados muestran una concentración de tensiones en el orificio,
alcanzando los -1836.2 kPa, muy similar al análisis sin refinamiento de la
malla.