Condiciones de contorno: restricciones y cargas en FEA

En el desarrollo de un modelo de Análisis de Elementos Finitos (FEA), la definición de las condiciones de contorno (boundary conditions) representa el paso donde se cometen más errores conceptuales. Puedes contar con una geometría CAD perfecta, un material bien caracterizado y una malla refinada con métricas de calidad impecables, pero si las restricciones o las cargas no representan fielmente la física real, el modelo entregará resultados completamente incorrectos.

Las condiciones de contorno son la interfaz que conecta la simulación numérica con el mundo físico: indican al software dónde se apoya la estructura, cómo interactúa con los componentes vecinos y de qué manera se aplican las fuerzas externas.

1. El problema matemático: matriz singular y movimientos de cuerpo rígido

Para comprender por qué las restricciones son obligatorias, debemos recurrir al álgebra lineal subyacente. La ecuación estática global del MEF se expresa como:

Antes de aplicar restricciones, la matriz de rigidez global [K] de una estructura desvinculada es singular: su determinante es exactamente igual a cero (det([K]) = 0) y, por lo tanto, no posee inversa.

Físicamente, esto significa que el sistema no tiene un punto de referencia fijo en el espacio. Si se aplicara una fuerza sobre la estructura no restringida, esta experimentaría un movimiento de cuerpo rígido (rigid body motion): se trasladaría o rotaría en el espacio de forma infinita sin sufrir deformación elástica ni generar tensiones internas.

Los grados de libertad y los modos de cuerpo rígido

Un cuerpo sólido tridimensional en el espacio posee 6 modos de cuerpo rígido:

• 3 traslaciones puras en los ejes x, y, z (T_x, T_y, T_z).
• 3 rotaciones puras alrededor de los ejes x, y, z (R_x, R_y, R_z).

Para que la matriz [K] se vuelva invertible y el solucionador pueda calcular los desplazamientos, el analista debe aplicar suficientes restricciones cinemáticas para eliminar como mínimo estos 6 modos de cuerpo rígido.

2. Clasificación de las condiciones de contorno

Las condiciones de contorno se dividen formalmente en dos categorías matemáticas primarias:

1. Condiciones esenciales o de Dirichlet (primarias): imponen valores conocidos sobre los desplazamientos o rotaciones en la frontera del dominio (por ejemplo, u = 0 en un empotramiento).
2. Condiciones naturales o de Neumann (secundarias): imponen valores conocidos sobre las fuerzas, presiones o derivadas de la variable de campo (por ejemplo, una carga puntual F o una presión distribuida p).

3. El arte de no sobrerrestringir: apoyos isostáticos y la regla 3-2-1

Un error habitual en el modelado es la sobrerrestricción (over-stiffening). Ocurre cuando se aplican más restricciones de las necesarias o cuando se fijan nodos de forma rígida a tierra en regiones que en la realidad deberían poder deformarse libremente.

Por ejemplo, si se bloquean todas las traslaciones en la base de un bloque sometido a flexión o expansión térmica, se impedirá la contracción transversal por efecto Poisson (ε_y = −ν · ε_x) o la dilatación natural del material, introduciendo tensiones parásitas artificiales que no existen en la pieza real.

La solución isostática: la regla 3-2-1

Para evitar el sobrerrefuerzo y garantizar que la estructura se mantenga en equilibrio estático sin introducir rigidez adicional, se utilizan restricciones isostáticas (isostatic restraints) mediante la regla 3-2-1. Esta regla consiste en eliminar los 6 grados de libertad de cuerpo rígido (3 traslaciones y 3 rotaciones) utilizando exactamente tres puntos de apoyo sobre el plano de referencia:

1. En el nodo 1 (punto A): se restringen las 3 traslaciones [T_x, T_y, T_z]. Esto anula las traslaciones globales de la estructura, pero permite giros libres alrededor del punto A.
2. En el nodo 2 (punto B, no coincidente con A): se restringen 2 traslaciones [T_y, T_z]. En combinación con el nodo 1, esto elimina las rotaciones alrededor de los ejes perpendiculares a la línea A-B.
3. En el nodo 3 (punto C, fuera de la línea A-B): se restringe 1 traslación [T_z]. Esto bloquea el último modo de rotación restante, el giro alrededor de la línea A-B.

De este modo, se eliminan estrictamente los 6 grados de libertad de cuerpo rígido sin añadir ningún hiperestatismo al modelo, permitiendo que la pieza se contraiga, se dilate o se flecte naturalmente bajo carga.

4. Compatibilidad entre elementos y restricciones

Un aspecto crítico es verificar la compatibilidad entre el tipo de elemento y los grados de libertad restringidos:

• Elementos sólidos 3D (TET4, TET10, HEX8, HEX20): solo poseen 3 grados de libertad de traslación por nodo (u, v, w). No tienen rigidez rotacional en sus nodos. Si intentas aplicar un empotramiento restringiendo giros (θ_x, θ_y, θ_z) en una sola fila de nodos sólidos, el solucionador ignorará la restricción rotacional y el modelo actuará como una rótula cilíndrica, generando un mecanismo inestable. Más detalle en elementos sólidos 3D.
• Elementos de cáscara (shells) y vigas (beams): poseen de 5 a 6 grados de libertad por nodo (traslaciones y rotaciones). En estos elementos sí es físicamente posible y válido restringir momentos o rotaciones directamente en los nodos. Ver elementos Shell y barras y vigas.

5. Aplicación de cargas realistas: evitando singularidades

Del mismo modo que las restricciones deben ser realistas, la aplicación de cargas externas condiciona la validez de las tensiones locales.

A. El peligro de las cargas puntuales nodales

En la teoría elástica, aplicar una fuerza concentrada F sobre un único nodo (área A = 0) produce teóricamente una tensión infinita (σ = F/0 = ∞). En la malla, esto genera una singularidad de tensión que crecerá sin límite a medida que se refine la malla. Aprende a reconocerla en el estudio de convergencia.

B. Uso de elementos de acoplamiento R-Type (RBE2 y RBE3)

Para aplicar fuerzas o momentos remotos sin generar singularidades en un solo nodo ni alterar artificialmente la rigidez de la zona de carga, se utilizan conectores tipo R-Type:

• RBE2 (elemento rígido de N nodos): conecta un nodo maestro con un conjunto de nodos dependientes mediante una rigidez infinita. Se utiliza para simular uniones soldadas rígidas o placas de reparto extremadamente masivas. Atención: inmoviliza los nodos dependientes entre sí, impidiendo la deformación local.
• RBE3 (elemento de interpolación): distribuye la fuerza y el momento aplicados en un nodo de referencia hacia los nodos dependientes mediante una media ponderada, sin añadir rigidez al modelo. Permite que la estructura se deforme, se flecte o alabee libremente bajo la carga. Es la herramienta ideal para aplicar fuerzas resultantes o pares de torsión en orificios y caras.

6. Verificación y control de calidad de las condiciones de contorno

Antes de dar por válidos los resultados de tensión, el analista debe aplicar dos comprobaciones fundamentales de control de calidad:

1. Comprobación del equilibrio de reacciones: la suma de las fuerzas y momentos externos aplicados debe ser igual y opuesta a la suma de las reacciones calculadas en los nodos restringidos:
   ΣF_aplicadas = −ΣF_reacciones
   Si existe discrepancia, significa que el modelo contiene inconsistencias en las cargas o que la matriz no ha convergido adecuadamente.
2. Análisis modal libre-libre (free-free modal check): se ejecuta un análisis modal del modelo sin aplicar ninguna condición de contorno ni carga. Si el modelo está bien conectado internamente, los primeros 6 modos de vibración deben tener frecuencias numéricamente nulas (0 Hz). Si el análisis muestra más de 6 modos nulos, significa que existen elementos desconectados o mecanismos internos «sueltos» en la estructura que fallarán al ejecutar el cálculo estático.

Conclusión

Aplicar la regla 3-2-1 para evitar sobrerrestricciones, respetar la compatibilidad de grados de libertad de los elementos y utilizar elementos de interpolación RBE3 para distribuir las cargas son las prácticas que distinguen un modelo numérico estable y confiable de uno inútil.