En este trabajo, se presenta, implementa y prueba una novedosa formulación de alto orden para placas piezoeléctricas multicapa basada en la combinación de métodos de Galerkin discontinuo de penalización interior de orden variable y teorías generales de placas por capas. La característica clave de la formulación es la posibilidad de ajustar el orden de las funciones de base tanto en la aproximación en el plano como en la expansión a través del espesor de las variables primarias, es decir, desplazamientos y potencial eléctrico. Los resultados obtenidos de la aplicación a los casos de prueba considerados muestran precisión y robustez, confirmando así la técnica desarrollada como una herramienta computacional complementaria para el análisis y diseño de dispositivos laminados inteligentes.