Wasfi Zakaria ha proporcionado una excelente descripción del enfoque que diferencia los elementos de primer orden de los de segundo orden.
Hay una complejidad sutil introducida en los elementos a medida que se vuelven de orden superior.
Veamos un triángulo en el espacio real.
La función de forma canónica en coordenadas reales para un elemento triangular lineal es:
P = a + bx + cy (3 parámetros y 3 nodos)
y
dP / dx = b o la deformación en la dirección x puede variar linealmente en y.
dP / dy = c o la deformación en la dirección y puede variar linealmente en x.
La función de forma canónica en coordenadas reales para un triángulo bilineal (segundo orden) es:
P = a + bx + cy + dx ^ 2 + ey ^ 2 + fxy (6 parámetros y 6 nodos)
y
dP / dx = b + dx + fv
dP / dy = c + ey + fx
Y nuevamente tenemos comportamientos de tensión simétrica.
Ahora veamos el elemento quad lineal:
P = a + bx + cy + dxy (cuatro parámetros, cuatro nodos)
y
dP / dx = b + dy
dP / dy = c + dx
Tenga en cuenta que existe una asimetría en los campos de deformación d / dx y d / dy.
Ahora veamos el elemento de casualidad biquadratica (ocho nodos):
P = a + bx + cy + dx ^ 2 + ey ^ 2 + fxy + gxy ^ 2 + hx ^ 2y (ocho parámetros, ocho nodos)
y los campos de deformación pueden ser determinados por
dP / dx = b + 2dx + fy + gy ^ 2 + 2hxy
dP / dy = c + 2ey + fx + 2gxy + hx ^ 2
y nuevamente los campos de deformación no son simétricos.
Por lo tanto, los elementos triangulares (y los elementos tetraédricos i 3D) tienen campos de deformación simétrica (y, por lo tanto, tensión), mientras que los elementos de serendipia cuádruple no.
¿Por qué eso importa?
Veamos un campo de desplazamiento constante puro (estiramiento constante). Todo el elemento solo exhibirá el término de deformación constante y todos se comportarán igualmente bien.
Veamos la deformación lineal a través de la sección (como se dice en flexión pura). el triángulo lineal es una deformación constante y, por lo tanto, coincide con la deformación real como un conjunto de funciones escalonadas y converge muy lentamente. Para ciertos problemas (plasticidad), estos elementos realmente se bloquean y se establecen correctamente, el comportamiento de convergencia es extraño. sin embargo, los elementos bilineales pueden representar explícitamente un campo de deformación que varía linealmente en x o y y los elementos convergen inmediatamente para un elemento.
Ahora echemos un vistazo a los campos de desplazamiento de orden superior, digamos un campo de desplazamiento cúbico que produce campos de deformación cuadráticos (flexión bajo carga final). El triángulo bilineal se ajustará al campo de desplazamiento con un conjunto de campos cuadráticos y la convergencia es relativamente rápida. Del mismo modo, la variación del campo de deformación se puede representar simétricamente a través del elemento y el campo de deformación se comporta bien. Veamos los elementos cuádruples. También asignarán el campo de desplazamiento como un conjunto de campos de desplazamiento cuadráticos y convergerán bastante rápido. Sin embargo, ahora hay componentes de deformación de segundo orden y estos pueden excitar los términos de segundo orden en la derivada de las funciones de forma. Y a medida que el campo de desplazamiento se vuelve más severo y complejo, estos campos de tensión de orden superior se excitan cada vez más. El resultado puede ser deformaciones oscilantes (y, por lo tanto, tensiones), ver más abajo.

tomado de:
Análisis estructural con el método de elementos finitos. Estática lineal
Esto se discute más en:
Alisamiento de deformación de mínimos cuadrados para el elemento de tensión del plano de serendipia de ocho nodos
y
Procedimientos de elementos finitos
y
Análisis estructural con el método de elementos finitos. Estática lineal
El suavizado de los mínimos cuadrados sobre el elemento (la línea recta en este caso) es una solución muy efectiva para este desafío.
Impacto:
1) los quads / rectángulos convergen más rápido que los triángulos / tetraedros
2) los elementos bilineales convergen mucho más rápido que los elementos lineales
3) los quads / rectángulos bilineales (o langrangianos o …) son susceptibles a las oscilaciones de estrés parasitario
4) el ajuste al cuadrado mínimo de los campos de esfuerzo / tensión sobre el elemento es muy efectivo para reducir esta oscilación