Aufgrund ihrer mechanischen Eigenschaften ersetzen Verbundwerkstoffe zunehmend klassische metallische Materialien. Aus diesem Grund stellt die Simulation von faserverstärkten Kunststoffen durch numerische Verfahren wie die Finite-Elemente-Methode ein stetig wachsendes Forschungsfeld dar. Diese Arbeit liefert einen Beitrag zur Modellierung dieser Materialien. Es wird ein Schalenelement zur Berechnung von geschichteten Strukturen vorgestellt, das auf einem Mehrfeldfunktional basiert. Das Verschiebungsfeld beinhaltet neben den Anteilen aus der Reissner-Mindlin-Theorie zusätzliche Verwölbungen und Verschiebungen senkrecht zur Schalenfläche. Durch diese Größen ist es möglich, interlaminare Schub- und Normalspannungen zu bestimmen. Die Verwölbungen und die Verschiebungen in Dickenrichtung der Schale werden in den lokalen Gleichgewichtsbeziehungen berücksichtigt. Somit ist die Anwendung auf physikalisch und geometrisch nichtlinearer Probleme möglich. Eine Nebenbedingung erzwingt den korrekten Verlauf des erweiterten Verschiebungsfeldes. Aus statischer Kondensation resultiert ein gemischtes 4-Knoten-Schalenelement mit den \textcolor{black}{Standard-Verschiebungsfreiheitsgraden}. Dies ermöglicht die Anwendung der Schalenformulierung auf komplexe Geometrien. Innerhalb des Laminats kann an beliebiger Stelle eine Zwischenschicht zur Abbildung von Delamination platziert werden. Ein isotropes Schädigungsmodell für kleine Verzerrungen wird in dieser Zwischenschicht für \textcolor{black}{die Bestimmung von Delaminationen} verwendet. Anhand mehrerer numerischer Beispiele wird unter anderem der Vorteil gegenüber der Verwendung von Volumenschalenelementen in Bezug auf die Rechenzeit verdeutlicht. Daneben werden die Berechnung interlaminarer Spannungen verschiedener Laminate, ein 3-Punkt-Biegeversuch unter Berücksichtigung von Delamination und eine Reihe weiterer Beispiele gezeigt.