Elementos Shell: los híbridos del FEA que todo analista debe dominar

En el Análisis de Elementos Finitos (FEA), la eficiencia computacional es tan importante como la precisión de los resultados. Para modelar estructuras tridimensionales de gran tamaño sin saturar los recursos de hardware, los ingenieros recurren frecuentemente a una de las herramientas más potentes y elegantes de la biblioteca de elementos: los elementos Shell (cáscara o lámina).

Estos elementos actúan como verdaderos híbridos entre el mundo 2D y el 3D. Geométricamente, se definen en un espacio tridimensional, lo que les permite orientarse en cualquier dirección y soportar cargas multidireccionales, pero se formulan matemáticamente sobre una superficie bidimensional. Esto significa que, en lugar de modelar físicamente el espesor de una estructura con miles de pequeños ladrillos sólidos, el espesor se trata simplemente como un parámetro numérico asociado al elemento.

¿Cuándo usar elementos Shell? (la regla de la esbeltez)

La decisión de cuándo reducir un modelo de sólido 3D a elemento Shell se rige por la relación de esbeltez de la geometría. Cuando una estructura tiene dos de sus dimensiones (longitud y ancho) que son mucho mayores que la tercera (el espesor), se clasifica como una estructura de pared delgada. La regla general en la ingeniería práctica es aplicar elementos Shell cuando la relación entre el ancho (o luz de la pieza) y el espesor es mayor que 20.

Ejemplos clásicos de estas estructuras abundan en diversas industrias:

Fuselajes y revestimientos de alas en la industria aeroespacial.
Depósitos, silos, tuberías y tanques de almacenamiento en la industria de procesos.
Chapas, paneles de carrocería y perfiles estampados en la automoción.

b / t > 20 ⟹ elemento Shell · b / t ≤ 20 ⟹ sólido 3D

¿Por qué no modelar estas piezas con sólidos 3D? Intentar modelar una chapa delgada con elementos sólidos suele ser un grave error metodológico. Para capturar correctamente los efectos de flexión sin sufrir de bloqueo por cortante (shear locking), necesitarías colocar al menos tres o cuatro elementos sólidos a lo largo del espesor, lo que dispararía de forma prohibitiva el número de grados de libertad, el tiempo de cálculo y los recursos de memoria.

Los dos mundos del elemento Shell: membrana y flexión

La verdadera potencia del elemento Shell reside en su formulación matemática, la cual combina de manera única dos comportamientos en un solo elemento de 5 o 6 grados de libertad por nodo: el comportamiento de membrana y el de placa (flexión).

Estado de membrana: ocurre cuando las cargas actúan de forma paralela al plano medio del elemento. En este estado, el elemento transmite tracción, compresión y cortante en su propio plano. Las deformaciones y las tensiones resultantes son completamente uniformes a lo largo de todo el espesor de la chapa.
Estado de flexión: ocurre cuando las cargas actúan de forma perpendicular al plano de la chapa. En este caso, el elemento se comporta como una placa sometida a flexión y torsión fuera de su plano. A diferencia del estado de membrana, las tensiones varían de forma lineal a través del espesor, alcanzando su valor máximo de tracción en una de las superficies externas, compresión máxima en la opuesta, y un valor nulo en la fibra neutra (el plano medio).

Fig. 1 — En estado de membrana (izq.) la tensión es constante en todo el espesor; en flexión (dcha.) varía linealmente, con tracción máxima en una cara, compresión máxima en la opuesta y valor nulo en la fibra neutra.

Al combinar estos dos estados, el elemento Shell es capaz de modelar con precisión tanto el comportamiento de un globo inflado (membrana pura) como el de una placa de acero que se dobla bajo el peso de una carga (flexión pura).

El arte del modelado: la superficie media (Middle Surface)

Para implementar con éxito un modelo Shell, el analista debe comprender tres conceptos geométricos críticos:

La superficie media (Middle Surface): debido a que la malla no tiene un volumen físico, se genera sobre el plano medio de la chapa real. Esto requiere que, al importar un modelo 3D de CAD, el ingeniero realice una preparación de la geometría para extraer las superficies medias (mid-planes). En zonas con cambios de espesor, es fundamental definir correctamente estos planos y, si es necesario, aplicar desfases geométricos (offsets) para alinear las superficies.
Las caras top (superior) y bottom (inferior): dado que el elemento se representa como una superficie 2D sin espesor visual en la malla, el solucionador FEA debe designar matemáticamente qué lado de la superficie representa la cara superior (+t/2 o z₁) y cuál la inferior (−t/2 o z₂). La orientación de estas caras viene determinada por la dirección del vector normal del elemento, el cual se calcula mediante la regla de la mano derecha basada en el orden de numeración de los nodos del elemento.

Fig. 2 — La superficie media no tiene volumen: el software distingue la cara top (+t/2) de la bottom (−t/2) según el vector normal, calculado con la regla de la mano derecha a partir del orden de numeración de los nodos (1→2→3→4).

La importancia crítica de definir correctamente las normales

Definir y comprobar la consistencia de las normales de tus elementos Shell antes de correr el análisis es un paso crítico por dos razones fundamentales:

La dirección de las cargas de presión: las cargas de presión en elementos Shell se aplican en función del vector normal. Una presión positiva empujará en la dirección de la normal (+z), mientras que una presión negativa lo hará en sentido opuesto (−z). Si las normales de tu depósito o fuselaje están desalineadas (algunos elementos apuntan hacia afuera y otros hacia adentro), la presión actuará de forma inconsistente, destruyendo la coherencia física del modelo.
La interpretación de resultados: cuando el software calcula las tensiones, las reporta por separado para la cara top (z₁) y bottom (z₂). Si un elemento de tu malla está invertido respecto a sus vecinos, el post-procesador mostrará las tensiones de tracción en la cara donde debería haber compresión, creando una discontinuidad artificial en tus mapas de colores que puede enmascarar fallos estructurales graves y falsear tu validación de diseño.

La normal no es un detalle del renderizado: es el eje que decide hacia dónde empuja la presión y qué cara es tracción cuando el mapa de tensiones empieza a mentir.

# Reglas prácticas — elementos Shell
shell.dof_nodo        = "5-6"      # 3 traslación + 2-3 rotación por nodo
shell.esbeltez        = ancho/espesor > 20  # umbral frente a sólido 3D

membrana.tension      = "constante"   # uniforme en todo el espesor
flexion.tension       = "lineal"      # max en top/bottom, nula en z=0

normal.regla          = "mano derecha" # orden de nodos → dirección +z (top)
presion.signo         = "+p empuja según +n, -p según -n"

Membrana y flexión, lado a lado

Estado	Tipo de carga	Distribución de tensión	Ejemplo típico
Membrana	En el plano medio de la chapa	Constante en todo el espesor	Depósito a presión interna, globo
Flexión	Perpendicular al plano de la chapa	Lineal: máxima en top/bottom, nula en la fibra neutra	Placa apoyada bajo carga puntual

En resumen, los elementos Shell son una obra de arte de la simplificación matemática que permiten resolver problemas masivos de manera óptima. Sin embargo, esta abstracción traslada la responsabilidad directamente al analista. Comprobar la consistencia de las normales, verificar la alineación del plano medio y entender la diferencia entre las tensiones de membrana y flexión son las claves para que tus simulaciones de pared delgada sean verdaderamente confiables. Recuerda que el software siempre te dará un resultado, pero es tu criterio de ingeniería el único capaz de validar si ese resultado es físicamente coherente. ¿Te gustaría que en el próximo artículo exploremos cómo conectar correctamente estos elementos Shell con sólidos 3D para evitar problemas de incompatibilidad de grados de libertad?

FIN

In Finite Element Analysis (FEA), computational efficiency matters just as much as the accuracy of the results. To model large three-dimensional structures without saturating hardware resources, engineers frequently reach for one of the most powerful and elegant tools in the element library: shell elements.

These elements act as true hybrids between the 2D and 3D worlds. Geometrically, they're defined in three-dimensional space, letting them orient in any direction and carry multidirectional loads, but they're formulated mathematically over a two-dimensional surface. That means that instead of physically modelling a structure's thickness with thousands of small solid bricks, the thickness is simply treated as a numerical parameter attached to the element.

When to use shell elements (the slenderness rule)

The decision of when to reduce a 3D solid model to a shell element is governed by the geometry's slenderness ratio. When a structure has two of its dimensions (length and width) far larger than the third (thickness), it's classified as a thin-walled structure. The general rule of thumb in practical engineering is to apply shell elements when the ratio between the part's width (or span) and its thickness is greater than 20.

Classic examples of these structures abound across industries:

Fuselages and wing skins in the aerospace industry.
Tanks, silos, pipes and storage vessels in the process industry.
Sheet metal, body panels and stamped profiles in the automotive industry.

b / t > 20 ⟹ shell element · b / t ≤ 20 ⟹ 3D solid

Why not model these parts with 3D solids? Trying to model a thin sheet with solid elements is usually a serious methodological mistake. To correctly capture bending effects without suffering from shear locking, you'd need to place at least three or four solid elements through the thickness, which drives up the number of degrees of freedom, computation time and memory usage prohibitively.

The shell element's two worlds: membrane and bending

The real power of the shell element lies in its mathematical formulation, which uniquely combines two behaviours in a single element with 5 or 6 degrees of freedom per node: membrane behaviour and plate (bending) behaviour.

Membrane state: occurs when loads act parallel to the element's mid-plane. In this state, the element carries tension, compression and shear within its own plane. The resulting strains and stresses are completely uniform across the entire thickness of the sheet.
Bending state: occurs when loads act perpendicular to the plane of the sheet. Here, the element behaves like a plate subjected to out-of-plane bending and twisting. Unlike the membrane state, stresses vary linearly through the thickness, reaching their maximum tension value on one outer surface, maximum compression on the opposite one, and zero at the neutral fibre (the mid-plane).

Fig. 1 — In the membrane state (left) stress is constant through the whole thickness; in bending (right) it varies linearly, with maximum tension on one face, maximum compression on the other, and zero at the neutral fibre.

By combining these two states, the shell element can accurately model both the behaviour of an inflated balloon (pure membrane) and that of a steel plate bending under a load (pure bending).

The art of modelling: the middle surface

To successfully build a shell model, the analyst must understand three critical geometric concepts:

The middle surface: because the mesh has no physical volume, it's generated on the mid-plane of the real sheet. This means that when importing a 3D CAD model, the engineer must prepare the geometry to extract the mid-planes. In zones with thickness changes, it's essential to define these planes correctly and, if needed, apply geometric offsets to align the surfaces.
The top and bottom faces: since the element is represented as a 2D surface with no visible thickness in the mesh, the FEA solver must mathematically designate which side of the surface is the top face (+t/2 or z₁) and which is the bottom (−t/2 or z₂). The orientation of these faces is determined by the direction of the element's normal vector, computed with the right-hand rule based on the element's node numbering order.

Fig. 2 — The middle surface has no volume: the solver tells the top face (+t/2) from the bottom (−t/2) by the normal vector, computed with the right-hand rule from the node numbering order (1→2→3→4).

Why getting your normals right is critical

Defining and checking the consistency of your shell elements' normals before running the analysis is a critical step, for two fundamental reasons:

The direction of pressure loads: pressure loads on shell elements are applied relative to the normal vector. A positive pressure will push in the direction of the normal (+z), while a negative pressure will push the opposite way (−z). If the normals on your tank or fuselage are misaligned — some elements pointing outward, others inward — the pressure will act inconsistently, destroying the model's physical coherence.
Interpreting results: when the software computes stresses, it reports them separately for the top face (z₁) and the bottom face (z₂). If one element in your mesh is flipped relative to its neighbours, the post-processor will show tensile stress on the face where there should be compression, creating an artificial discontinuity in your colour maps that can mask serious structural failures and invalidate your design check.

The normal isn't a rendering detail: it's the axis that decides which way the pressure pushes and which face reads as tension when the stress map starts lying to you.

# Practical rules — shell elements
shell.dof_node        = "5-6"      # 3 translation + 2-3 rotation per node
shell.slenderness     = width/thickness > 20  # threshold vs. 3D solid

membrane.stress       = "constant"   # uniform through thickness
bending.stress        = "linear"      # max at top/bottom, zero at z=0

normal.rule           = "right-hand"  # node order → +z direction (top)
pressure.sign         = "+p pushes along +n, -p along -n"

Membrane and bending, side by side

State	Load type	Stress distribution	Typical example
Membrane	Within the sheet's mid-plane	Constant through the thickness	Pressurised tank, balloon
Bending	Perpendicular to the sheet's plane	Linear: max at top/bottom, zero at the neutral fibre	Plate under a point load

In short, shell elements are a work of art in mathematical simplification, letting you solve massive problems in an optimal way. But that abstraction shifts the responsibility straight onto the analyst. Checking the consistency of your normals, verifying the alignment of the mid-plane, and understanding the difference between membrane and bending stresses are the keys to making your thin-wall simulations genuinely trustworthy. Remember: the software will always hand you a result, but your engineering judgement is the only thing that can validate whether that result is physically coherent. Should the next article cover how to properly connect these shell elements to 3D solids, avoiding degree-of-freedom incompatibility issues?

END