Mallado en FEA: calidad, densidad y plano medio

En el Análisis de Elementos Finitos (FEA), la fase de preprocesamiento —y en particular la generación de la malla— suele consumir hasta el 80% del tiempo total dedicado a un proyecto de simulación. Aunque las herramientas CAD y los generadores automáticos modernos permiten crear mallas con solo presionar un botón, el mallado sigue siendo un arte respaldado por la física y las matemáticas.

Una malla no es simplemente una red gráfica de colores; es el medio físico que discretiza el continuo en subdominios simples para evaluar el sistema de ecuaciones algebraicas del solucionador. Un diseño brillante con condiciones de contorno perfectas entregará resultados completamente erróneos si la malla es deficiente o numéricamente inestable.

En este artículo abordaremos las reglas prácticas para construir mallas de alta calidad, el concepto fundamental del plano medio (mid-plane) para estructuras delgadas y las métricas indispensables para evitar la distorsión numérica.

1. Densidad de malla y estrategia de refinamiento

La densidad de malla (el tamaño y cantidad de elementos) condiciona directamente la precisión de la solución frente a la realidad. Sin embargo, aumentar indiscriminadamente el número de nodos no está justificado: el consumo de memoria RAM escala con el cuadrado de los grados de libertad (N²) y el tiempo de cálculo lo hace aproximadamente con el cubo (N³).

Fig. 1 — Malla gruesa (izq.) frente a refinamiento local (dcha.): la densidad se concentra solo alrededor del agujero, donde está el gradiente de tensión, y se ahorran grados de libertad en el resto de la pieza.

Reglas de refinamiento local

Mallas gruesas (coarse): se emplean en fases preliminares del diseño, para la evaluación de deformaciones globales o en zonas alejadas de las áreas de interés donde los gradientes de tensión son bajos.
Mallas refinadas (fine): se aplican en regiones con altos gradientes de tensión o discontinuidades geométricas: agujeros, entallas, redondeos (fillets), esquinas reentrantes y zonas de aplicación de carga.
Transiciones suaves: la proporción de tamaño entre un elemento y su vecino no debe superar una relación de 1:3 para evitar saltos numéricos bruscos de rigidez que generen tensiones ficticias.

El estudio de convergencia de malla

Para garantizar que el tamaño de los elementos no esté falseando los resultados, el analista debe realizar un estudio de convergencia:

Se resuelve el modelo con una malla inicial promedio.
Se reduce progresivamente el tamaño de los elementos (método h) exclusivamente en las regiones críticas.
Se grafica la métrica de interés (generalmente la tensión de Von Mises o el desplazamiento) frente al número de nodos o tamaño de elemento.
La malla se considera convergida cuando la variación en los resultados entre dos iteraciones consecutivas es menor al 5%.

2. El concepto del plano medio (mid-plane) para chapas delgadas

Modelar piezas de pared delgada (como chapas estampadas, depósitos o paneles aeronáuticos) utilizando elementos sólidos 3D es uno de los errores más frecuentes. Para capturar la flexión en un sólido 3D sin sufrir bloqueo por cortante (shear locking), se requieren al menos de 3 a 4 elementos a lo largo del espesor, lo que infla la cantidad de ecuaciones a niveles inmanejables.

La solución idónea es el uso de elementos de cáscara 2D (Shells), los cuales requieren que la malla se construya sobre el plano medio (mid-plane) de la pieza.

Fig. 2 — La malla se construye sobre el plano medio y el software reparte ±t/2 a cada lado (izq.). Si se mallara sobre una cara externa, el espesor quedaría desplazado t/2 respecto a la pieza real (dcha.).

Concepto físico

Matemáticamente, el software distribuye la mitad del espesor asignado (+t/2) hacia la cara superior (Top) y la otra mitad (−t/2) hacia la cara inferior (Bottom) a partir de la superficie 2D mallada. Si la malla se genera sobre la cara externa o interna en lugar del plano medio, las dimensiones geométricas reales quedarán sobreestimadas o subestimadas, alterando la rigidez del modelo.

Las dos reglas de oro del plano medio

Compatibilidad en intersecciones: las superficies que se cruzan en el plano medio deben dividirse explícitamente en el CAD/preprocesador a lo largo de sus líneas de contacto para forzar que el algoritmo de mallado cree nodos compartidos en las uniones.
Espesores variables y offsets: en piezas escalonadas o mecanizadas con cambios de espesor, se debe definir el plano medio en el espesor mínimo y aplicar un parámetro de desvío (z-offset) para situar correctamente la fibra neutra del material.

3. Métricas de calidad de elementos: evitando la distorsión numérica

La precisión de los resultados en cada elemento depende de cuánto se desvíe su forma geométrica respecto al elemento isoparamétrico ideal (el cuadrado en 2D o el cubo/tetraedro regular en 3D). Elementos excesivamente distorsionados provocan un mal condicionamiento de la matriz de rigidez global [K], lo que se traduce en errores de cálculo o divergencia del solucionador.

Las principales métricas de calidad que se deben verificar en el preprocesador son:

Métrica de calidad	Definición y forma ideal	Límite aceptable recomendado	Impacto en el análisis
Relación de aspecto (Aspect ratio)	Cociente entre la dimensión más larga y la más corta del elemento. Ideal = 1.0 (cuadrado).	< 5.0 (hasta < 3.0 en zonas críticas)	Valores altos reducen drásticamente la precisión en gradientes de tensión y flexión.
Ángulo de sesgo (Skewness / Skew angle)	Desviación de los ángulos internos respecto a 90° (cuadriláteros) o 60° (triángulos). Ideal = 0°.	< 45°	Elementos muy agudos u obtusos alteran la interpolación de las funciones de forma.
Ángulo de alabeo (Warp angle)	Desviación fuera del plano para cuadriláteros 2D. Ideal = 0° (plano perfecto).	< 10°	Genera acoplamientos numéricos falsos entre membrana y flexión.
Jacobiano (Jacobian determinant)	Medida de la distorsión matemática al mapear el elemento real al sistema natural. Ideal = 1.0.	> 0.5 a 0.6 (valores ≤ 0 son inaceptables)	Un Jacobiano negativo indica un elemento invertido que detendrá la simulación.
Conicidad (Taper)	Desviación de la forma de cuadrilátero hacia una forma triangular. Ideal = 0.	< 0.5	Afecta la precisión del cálculo de deformaciones en los puntos de Gauss.

Fig. 3 — Elemento ideal (fila superior) frente a elemento distorsionado (inferior) para cada métrica, con el límite aceptable recomendado.

4. Los 4 errores de mallado 2D más comunes

Durante el proceso de discretización de superficies, los analistas suelen enfrentarse a cuatro errores recurrentes que deben ser detectados y corregidos:

Desviación geométrica (Geometrical deviation): ocurre cuando la malla no se adapta fielmente a la curvatura del CAD. Se debe controlar la tolerancia de cuerda (chord deviation) para evitar faceteados toscos en radios y cilindros.
Normales desalineadas (Inconsistent normals): todos los elementos Shell de una superficie continua deben tener sus vectores normales apuntando en la misma dirección. Normales invertidas causan que las presiones se apliquen en sentidos opuestos y que los resultados de tensión en las caras Top y Bottom se muestren intercambiados.
Bordes libres y nodos duplicados (Free edges / Duplicate nodes): ocurre cuando superficies adyacentes no comparten nodos en sus límites comunes. Esto crea una «grieta» o discontinuidad artificial en el modelo que destruye el camino de la carga.
Elementos duplicados (Duplicate elements): creación accidental de dos elementos superpuestos en la misma ubicación geométrica. Esto duplica artificialmente la rigidez y la masa de esa zona.

Fig. 4 — Dos errores típicos del mallado 2D: nodos duplicados que dejan un borde libre entre superficies adyacentes (izq.) y normales desalineadas dentro de una misma superficie (dcha.).

5. Reglas prácticas para crear mallas elegantes e inteligibles

Una malla bien estructurada no solo es visualmente limpia, sino que garantiza una distribución uniforme del error numérico.

Priorizar cuadriláteros (QUAD) sobre triángulos (TRIA): los elementos triangulares de primer orden (TRIA3) poseen deformación constante y sufren de excesiva rigidez. Se debe asegurar que al menos el 90% de los elementos de la malla 2D sean QUADs, reservando los triángulos únicamente para zonas de transición alejadas de las regiones de interés.
Modelado de agujeros:
- Usar siempre un número par de elementos alrededor del contorno (6, 8, 12, 16…) para mantener la simetría y las líneas de flujo de la malla.
- En regiones no críticas basta con 6 elementos; en regiones de concentración de tensiones se requieren al menos 12 elementos alrededor del agujero.
- El refinamiento radial debe extenderse a una distancia de 1.5 a 2 veces el diámetro del agujero, colocando al menos dos capas de elementos finitos en el perímetro.
Inspección sin promediar (Unsmoothed / Non-averaged plots): para evaluar la calidad de la malla tras el cálculo, se deben visualizar los mapas de tensiones sin promediar los nodos. Discontinuidades o saltos de color marcados entre elementos vecinos indican que la malla es demasiado gruesa en esa zona y requiere mayor refinamiento.

Fig. 5 — Número par de elementos alrededor del agujero: 6 bastan en zonas no críticas (izq.); en concentración de tensiones se usan al menos 12 y dos capas radiales (dcha.).

Fig. 6 — Tensiones sin promediar entre elementos vecinos: un salto de color marcado (izq., 150 → 220 MPa) delata una malla tosca; una transición suave (dcha., 195 → 200 MPa) indica una malla adecuada y convergida.

Conclusión

El mallado no debe considerarse un mero trámite antes de presionar el botón de cálculo (Solve), sino el paso donde el ingeniero aplica su verdadero criterio técnico. Ajustar las métricas de aspecto y sesgo, definir correctamente el plano medio en chapas delgadas y realizar estudios de convergencia son las garantías para obtener resultados representativos y confiables.

En el modelado estructural se cumple la máxima de la simulación:

Si la malla parece limpia y ordenada, la solución numérica será robusta (If it looks good, it is good).

FIN

In Finite Element Analysis (FEA), the preprocessing phase —and meshing in particular— can consume up to 80% of the total time spent on a simulation project. Although modern CAD tools and automatic mesh generators can create meshes at the push of a button, meshing is still an art backed by physics and mathematics.

A mesh is not just a coloured graphic grid; it is the physical medium that discretises the continuum into simple subdomains so that the solver's system of algebraic equations can be evaluated. A brilliant design with perfect boundary conditions will deliver completely wrong results if the mesh is poor or numerically unstable.

In this article we cover the practical rules for building high-quality meshes, the fundamental concept of the mid-plane for thin structures, and the essential metrics for avoiding numerical distortion.

1. Mesh density and refinement strategy

Mesh density (the size and number of elements) directly determines how close the solution is to reality. Increasing the number of nodes indiscriminately is not justified, however: RAM consumption scales with the square of the degrees of freedom (N²) and computation time roughly with the cube (N³).

Fig. 1 — Coarse mesh (left) versus local refinement (right): density is concentrated only around the hole, where the stress gradient is, saving degrees of freedom in the rest of the part.

Local refinement rules

Coarse meshes: used in early design stages, for evaluating global deformations, or in regions far from the areas of interest where stress gradients are low.
Fine meshes: applied in regions with high stress gradients or geometric discontinuities: holes, notches, fillets, re-entrant corners and load application zones.
Smooth transitions: the size ratio between an element and its neighbour should not exceed 1:3, to avoid abrupt numerical jumps in stiffness that create spurious stresses.

The mesh convergence study

To make sure the element size is not distorting the results, the analyst should carry out a convergence study:

Solve the model with an average initial mesh.
Progressively reduce the element size (h-method) exclusively in the critical regions.
Plot the metric of interest (usually the von Mises stress or the displacement) against the number of nodes or the element size.
The mesh is considered converged when the variation in the results between two consecutive iterations is less than 5%.

2. The mid-plane concept for thin sheets

Modelling thin-walled parts (such as stamped sheets, tanks or aircraft panels) with 3D solid elements is one of the most frequent mistakes. To capture bending in a 3D solid without suffering shear locking, at least 3 to 4 elements through the thickness are required, which inflates the number of equations to unmanageable levels.

The right solution is to use 2D shell elements, which require the mesh to be built on the part's mid-plane.

Fig. 2 — The mesh is built on the mid-plane and the software distributes ±t/2 on each side (left). If it were meshed on an outer face, the thickness would end up shifted by t/2 relative to the real part (right).

Physical concept

Mathematically, the software distributes half of the assigned thickness (+t/2) towards the top face and the other half (−t/2) towards the bottom face, starting from the meshed 2D surface. If the mesh is generated on the outer or inner face instead of the mid-plane, the real geometric dimensions will be over- or underestimated, altering the stiffness of the model.

The two golden rules of the mid-plane

Compatibility at intersections: surfaces that cross at the mid-plane must be explicitly split in the CAD/preprocessor along their contact lines, to force the meshing algorithm to create shared nodes at the joints.
Variable thicknesses and offsets: in stepped or machined parts with thickness changes, the mid-plane must be defined at the minimum thickness and an offset parameter (z-offset) applied to place the neutral fibre of the material correctly.

3. Element quality metrics: avoiding numerical distortion

The accuracy of the results in each element depends on how far its geometric shape deviates from the ideal isoparametric element (the square in 2D or the regular cube/tetrahedron in 3D). Excessively distorted elements cause ill-conditioning of the global stiffness matrix [K], which results in calculation errors or solver divergence.

The main quality metrics to check in the preprocessor are:

Quality metric	Definition and ideal shape	Recommended acceptable limit	Impact on the analysis
Aspect ratio	Ratio between the longest and the shortest dimension of the element. Ideal = 1.0 (square).	< 5.0 (down to < 3.0 in critical zones)	High values drastically reduce accuracy in stress and bending gradients.
Skewness / Skew angle	Deviation of the internal angles from 90° (quadrilaterals) or 60° (triangles). Ideal = 0°.	< 45°	Very sharp or obtuse elements distort the interpolation of the shape functions.
Warp angle	Out-of-plane deviation for 2D quadrilaterals. Ideal = 0° (perfectly flat).	< 10°	Creates spurious numerical coupling between membrane and bending behaviour.
Jacobian determinant	Measure of the mathematical distortion when mapping the real element to the natural coordinate system. Ideal = 1.0.	> 0.5 to 0.6 (values ≤ 0 are unacceptable)	A negative Jacobian indicates an inverted element that will stop the simulation.
Taper	Deviation of a quadrilateral towards a triangular shape. Ideal = 0.	< 0.5	Affects the accuracy of the strain calculation at the Gauss points.

Fig. 3 — Ideal element (top row) versus distorted element (bottom row) for each metric, with the recommended acceptable limit.

4. The 4 most common 2D meshing errors

During surface discretisation, analysts regularly run into four recurring errors that must be detected and corrected:

Geometrical deviation: occurs when the mesh does not follow the curvature of the CAD faithfully. The chord deviation tolerance must be controlled to avoid coarse faceting on radii and cylinders.
Inconsistent normals: all shell elements on a continuous surface must have their normal vectors pointing in the same direction. Reversed normals cause pressures to be applied in opposite directions and the stress results on the Top and Bottom faces to appear swapped.
Free edges / duplicate nodes: occurs when adjacent surfaces do not share nodes along their common boundaries. This creates an artificial “crack” or discontinuity in the model that destroys the load path.
Duplicate elements: accidental creation of two overlapping elements at the same geometric location. This artificially doubles the stiffness and mass of that zone.

Fig. 4 — Two typical 2D meshing errors: duplicate nodes that leave a free edge between adjacent surfaces (left) and misaligned normals within the same surface (right).

5. Practical rules for clean, readable meshes

A well-structured mesh is not only visually clean; it also ensures a uniform distribution of numerical error.

Prefer quadrilaterals (QUAD) over triangles (TRIA): first-order triangular elements (TRIA3) have constant strain and are excessively stiff. Make sure that at least 90% of the 2D mesh elements are QUADs, reserving triangles only for transition zones away from the regions of interest.
Modelling holes:
- Always use an even number of elements around the contour (6, 8, 12, 16…) to preserve the symmetry and the flow lines of the mesh.
- In non-critical regions 6 elements are enough; in stress concentration regions at least 12 elements around the hole are required.
- Radial refinement should extend to a distance of 1.5 to 2 times the diameter of the hole, placing at least two layers of finite elements on the perimeter.
Unsmoothed (non-averaged) plots: to assess mesh quality after the calculation, display the stress maps without nodal averaging. Marked discontinuities or colour jumps between neighbouring elements indicate that the mesh is too coarse in that zone and needs further refinement.

Fig. 5 — Even number of elements around the hole: 6 are enough in non-critical zones (left); at stress concentrations use at least 12 and two radial layers (right).

Fig. 6 — Unsmoothed stresses between neighbouring elements: a marked colour jump (left, 150 → 220 MPa) betrays a coarse mesh; a smooth transition (right, 195 → 200 MPa) indicates an adequate, converged mesh.

Conclusion

Meshing should not be treated as a mere formality before pressing the Solve button, but as the step where the engineer applies real technical judgement. Tuning the aspect and skew metrics, correctly defining the mid-plane in thin sheets and running convergence studies are the guarantees of representative and reliable results.

In structural modelling, the maxim of simulation holds true:

If the mesh looks clean and orderly, the numerical solution will be robust (If it looks good, it is good).

END