Gegenstand der vorliegenden Arbeit ist die numerische Simulation piezoelektrischen Materialverhaltens, sowie spezieller Probleme aus der Mechanik (inkompressibles Materialverhalten) unter Anwendung der Methode der finiten Elemente. Hierbei wird die Strategie der adaptiven Netzsteuerung angewendet, welche mit Hilfe einer lokalisierten a-posteriori Fehlerschätzung erlaubt, den lokalen Feinheitsgrad der Diskretisierung den Besonderheiten der Aufgabenstellung anzupassen. Beide betrachteten Problemklassen führen nach der Diskretisierung und FEM auf Gleichungssysteme in spezieller Blockstruktur, die insgesamt symmetrisch, aber nicht positiv definit ist. Als Löser kann nicht der gewöhnliche CG verwendet werden, stattdessen wird eine Variante des Bramble-Pasciak-CGs benutzt, welcher als Speziallöser die Matrizenstruktur ausnutzt. Für diesen Löser wird eine Strategie zur Parameterwahl vorgeschlagen sowie die Wirksamkeit einer Vorkonditionierung im piezoelektrischen Fall theoretisch nachgewiesen. Weiterhin wird die FEM einschließlich Adaptivität für Piezomaterialien auf rotationssymmetrische Probleme erweitert, so daß diese spezielle Problemklasse zweidimensional gerechnet werden kann. Numerische Vergleiche mit echter 3D-Rechnung illustrieren enorme Vorteile in Genauigkeit und Rechenaufwand. Im letzten Kapitel werden in piezoelektrische Materialien hineinwachsende Risse betrachtet und entsprechende Anpassungen vorgenommen. Mit Wahl geeigneter Datenstrukturen und einer passenden Vorkonditionierung ist es möglich, eine Simulationssoftware bereitzustellen welche als Grundlage zum Test von Bruchkriterien verwendet werden kann. Die beschriebenen numerischen Methoden wurden in ein bestehendes adaptives 2D-FEM-Programm implementiert, und an ausgewählten Beispielen ein Vergleich mit einer analytischen Lösung durchgeführt sowie die Effektivität der Rechnung getestet.

Gegenstand der vorliegenden Arbeit ist die numerische Simulation piezoelektrischen Materialverhaltens, sowie spezieller Probleme aus der Mechanik (inkompressibles Materialverhalten) unter Anwendung der Methode der finiten Elemente. Hierbei wird die Strategie der adaptiven Netzsteuerung angewendet, welche mit Hilfe einer lokalisierten a-posteriori Fehlerschätzung erlaubt, den lokalen Feinheitsgrad der Diskretisierung den Besonderheiten der Aufgabenstellung anzupassen. Beide betrachteten Problemklassen führen nach der Diskretisierung und FEM auf Gleichungssysteme in spezieller Blockstruktur, die insgesamt symmetrisch, aber nicht positiv definit ist. Als Löser kann nicht der gewöhnliche CG verwendet werden, stattdessen wird eine Variante des Bramble-Pasciak-CGs benutzt, welcher als Speziallöser die Matrizenstruktur ausnutzt. Für diesen Löser wird eine Strategie zur Parameterwahl vorgeschlagen sowie die Wirksamkeit einer Vorkonditionierung im piezoelektrischen Fall theoretisch nachgewiesen. Weiterhin wird die FEM einschließlich Adaptivität für Piezomaterialien auf rotationssymmetrische Probleme erweitert, so daß diese spezielle Problemklasse zweidimensional gerechnet werden kann. Numerische Vergleiche mit echter 3D-Rechnung illustrieren enorme Vorteile in Genauigkeit und Rechenaufwand. Im letzten Kapitel werden in piezoelektrische Materialien hineinwachsende Risse betrachtet und entsprechende Anpassungen vorgenommen. Mit Wahl geeigneter Datenstrukturen und einer passenden Vorkonditionierung ist es möglich, eine Simulationssoftware bereitzustellen welche als Grundlage zum Test von Bruchkriterien verwendet werden kann. Die beschriebenen numerischen Methoden wurden in ein bestehendes adaptives 2D-FEM-Programm implementiert, und an ausgewählten Beispielen ein Vergleich mit einer analytischen Lösung durchgeführt sowie die Effektivität der Rechnung getestet.

Optimal control problems for finite-strain elasticity are considered. An inner pressure or an inner fiber tension is acting as a driving force. Such internal forces are typical, for instance, for the motion of heliotropic plants, and for muscle tissue. Non-standard objective functions relevant for elasticity problems are introduced. Optimality conditions are derived on a formal basis, and a limited-memory quasi-Newton algorithm for their solution is formulated in function space. Numerical experiments confirm the expected mesh-independent performance.

Optimal control problems for finite-strain elasticity are considered. An inner pressure or an inner fiber tension is acting as a driving force. Such internal forces are typical, for instance, for the motion of heliotropic plants, and for muscle tissue. Non-standard objective functions relevant for elasticity problems are introduced. Optimality conditions are derived on a formal basis, and a limited-memory quasi-Newton algorithm for their solution is formulated in function space. Numerical experiments confirm the expected mesh-independent performance.