Quiero integrar una expresión polinómica sobre un elemento de 4 nodos en 3D. Varios libros sobre FEA cubren el caso en el que la integración se realiza sobre un elemento plano arbitrario de 4 nidos. El procedimiento habitual en este caso es encontrar la matriz de Jacobi y usar su determinante para cambiar la base de integración a la normalizada en la que tengo los límites de integración más simples [-1; 1] y la técnica de cuadratura de Gauss-Legendre se usa fácilmente.
En otras palabras, se reduce a la forma de
Pero en el caso 2D, cambio el elemento arbitrario plano al plano plano pero cuadrado bien formado 2 por 2.
El elemento 3D de 4 nodos no es plano en general, pero supongo que todavía se puede mapear con un sistema de coordenadas 2D que de alguna manera está relacionado con el sistema de coordenadas cartesianas. No puedo entender cómo expresar {x, y, z} en términos de {e, n} y cuál sería el tamaño de la matriz de Jacobi en este caso (se supone que es cuadrado).