Compuestos de fibra de carbono: del laminado manual a la optimizacion FEA

El laminado no es chapa

Un compuesto de fibra de carbono no es un material isotropo como el aluminio. Es un apilado de laminas unidireccionales, cada una con rigidez extrema en la direccion de fibra (E1 ~ 120-180 GPa) y dos ordenes de magnitud menos en direccion transversal (E2 ~ 7-10 GPa). La secuencia de apilado determina todo: rigidez global, resistencia, modo de fallo y peso.

El disenador de composites dispone de demasiada libertad: elegir numero de capas, orientacion de cada una (0, 90, +45, -45, o cualquier angulo intermedio), material de cada lamina, y secuencia de apilado. La optimizacion manual es inviable mas alla de 10-15 capas. El FEA permite explorar este espacio de diseno de forma sistematica.

Teoria clasica de laminados (CLT)

La CLT asume que cada lamina esta en estado de tension plana y que las deformaciones son continuas a traves del espesor. La relacion constitutiva del laminado completo es:

{N}   [A  B] {epsilon_0}
{ } = [    ] {         }
{M}   [B  D] {kappa}

Donde:

• A: matriz de rigidez membranal (responde a cargas en el plano)
• B: matriz de acoplamiento membrana-flexion. Si B != 0, una carga de traccion produce curvatura no deseada
• D: matriz de rigidez a flexion

Un laminado simetrico ([0/90/45/-45]s) anula B, eliminando acoplamientos parasitos. Esta es la razon por la que los laminados simetricos son el estandar industrial.

Criterios de fallo

A diferencia de metales (Von Mises), los composites requieren criterios especificos que distinguen direccion de fibra, matriz y compresion/traccion:

• Tsai-Wu: criterio cuadratico interactivo. Simple pero no distingue modos de fallo
• Tsai-Hill: generalizacion del criterio de Von Mises a materiales anisotropos
• Hashin 3D: distingue 4 modos de fallo (fibra/traccion, fibra/compresion, matriz/traccion, matriz/compresion). Es el mas usado en industria aeroespacial
• Puck: basado en fisica de fallo interfacial fibra-matriz. Mas preciso, mas parametros

# Criterio de Hashin para fallo de fibra en traccion
def hashin_fiber_tension(sigma_11, sigma_12, sigma_13, Xt, S12, S13):
    return (sigma_11 / Xt)**2 + (sigma_12**2 + sigma_13**2) / S12**2

# donde Xt = resistencia a traccion en direccion fibra
#       S12, S13 = resistencias a cortante interlaminar

Delaminacion y fallo interlaminar

El talon de Aquiles de los composites laminados: la resina entre capas es debil en tension perpendicular (Z). Las tensiones interlaminares — cortante tau_xz, tau_yz y tension normal sigma_zz — pueden iniciar delaminacion incluso con cargas en el plano bajas.

Para capturar este fenomeno necesitamos elementos 3D (solid-shell o brick) con al menos un elemento por capa. Modelar un laminado de 20 capas con elementos solidos implica mallas de facilmente 500k-1M de elementos para una pieza simple.

Alternativa: elementos cohesivos (cohesive zone model, CZM) insertados entre capas, que permiten iniciacion y propagacion de grieta interlaminar usando curvas de traccion-separacion bilineales.

Optimizacion topologica + orientacion de fibra

El estado del arte combina dos optimizaciones anidadas:

1. Topologica: donde poner material y donde no (variables de densidad por elemento)
2. Orientacion de fibra: para cada punto con material, cual es el angulo optimo (variable continua por elemento)

Esto genera disenos biologicos — con fibras fluyendo siguiendo trayectorias de carga principales — que recuerdan a la madera, el hueso o los tendones. La biomimetica emerge espontaneamente del algoritmo matematico.

Verificacion experimental

Validamos el modelo FEA con ensayos de flexion en 4 puntos sobre probetas de 8 capas (CFRP T700/M21, curado en autoclave a 180C). La correlacion fuerza-desplazamiento entre FEA y ensayo tiene R2 = 0.97 hasta primera carga de fallo. La prediccion de carga de primera lamina fallada (FPF) esta dentro del +-5% del valor experimental.

Conclusiones

El FEA en composites no es un complemento opcional: es la unica forma de navegar el enorme espacio de diseno que ofrecen los materiales anisotropos. La CLT da una primera aproximacion rapida; el modelado 3D con criterios de Hashin y CZM captura los modos de fallo reales. La optimizacion conjunta topologia-orientacion promete disenos que ningun ingeniero dibujaria a mano.