Inhalt der Arbeit ist die numerische Simulation von Systemen mit stochastischen Parametern, die durch stochastische partielle Differentialgleichungen (SPDGLn) beschrieben werden. Es werden die Theorie linearer und nichtlinearer elliptischer SPDGLn sowie Diskretisierungsverfahren beschrieben. Für die räumliche Diskretisierung wird eine existierende Simulationssoftware verwendet, während die stochastische Diskretisierung durch die direkte numerische Integration von Statistiken unter Verwendung von Monte Carlo- und Smolyak-Quadraturverfahren oder durch Reihenentwicklungen in Tensorprodukten finiter Elemente und stochastischer Ansatzfunktionen erfolgt. Die Reihenentwicklung wird dabei durch orthogonale Projektionen oder durch Galerkinverfahren gewonnen. Bei der Anwendung stochastischer Galerkinvervahren entstehen große Systeme gekoppelter Blockgleichungssysteme, welche hier durch iterative Verfahren gelöst werden. Zur Lösung linearer SPDGln werden effiziente Darstellungen der Gleichungssysteme und iterative Löser entwickelt. Aufgrund der Größe der entstehenden Gleichungssysteme wird ein paralleler Löser bereitgestellt. Die Lösung nichtlinearer SPDGLn geschieht durch approximative und Quasi-Newtonverfahren. Ein duales Verfahren ermöglicht die adaptive Verfeinerung der Lösung. Diese Verfahren werden in einer Allzwecksoftware für stochastische finite Elemente implementiert, die es erlaubt, existierende Simulationscodes um stochastische Unsicherheiten zu erweitern.