Convergent Numerical Schemes for the Compressible Hyperelastic Rod Wave Equation

We propose a fully discretised numerical scheme for the hyperelastic rod wave equation on the line. The convergence of the method is established. Moreover, the scheme can handle the blow-up of the derivative which naturally occurs for this equation. By using a time splitting integrator which preserves the invariants of the problem, we can also show that the scheme preserves the positivity of the energy density

Convergent numerical schemes for the compressible hyperelastic rod wave equation

We propose a fully discretised numerical scheme for the hyperelastic rod wave equation on the line. The convergence of the method is established. Moreover, the scheme can handle the blow-up of the derivative which naturally occurs for this equation. By using a time splitting integrator which preserves the invariants of the problem, we can also show that the scheme preserves the positivity of the energy densit

Theoretical and numerical comparison of hyperelastic and hypoelastic formulations for Eulerian non-linear elastoplasticity

The aim of this paper is to compare a hyperelastic with a hypoelastic model describing the Eulerian dynamics of solids in the context of non-linear elastoplastic deformations. Specifically, we consider the well-known hypoelastic Wilkins model, which is compared against a hyperelastic model based on the work of Godunov and Romenski. First, we discuss some general conceptual differences between the two approaches. Second, a detailed study of both models is proposed, where differences are made evident at the aid of deriving a hypoelastic-type model corresponding to the hyperelastic model and a particular equation of state used in this paper. Third, using the same high order ADER Finite Volume and Discontinuous Galerkin methods on fixed and moving unstructured meshes for both models, a wide range of numerical benchmark test problems has been solved. The numerical solutions obtained for the two different models are directly compared with each other. For small elastic deformations, the two models produce very similar solutions that are close to each other. However, if large elastic or elastoplastic deformations occur, the solutions present larger differences.Comment: 14 figure

A numerical study of variational discretizations of the Camassa-Holm equation

We present two semidiscretizations of the Camassa-Holm equation in periodic domains based on variational formulations and energy conservation. The first is a periodic version of an existing conservative multipeakon method on the real line, for which we propose efficient computation algorithms inspired by works of Camassa and collaborators. The second method, and of primary interest, is the periodic counterpart of a novel discretization of a two-component Camassa-Holm system based on variational principles in Lagrangian variables. Applying explicit ODE solvers to integrate in time, we compare the variational discretizations to existing methods over several numerical examples.Comment: 45 pages, 14 figure

Novel Discretization Schemes for the Numerical Simulation of Membrane Dynamics

Motivated by the demands of simulating flapping wings of Micro Air Vehicles, novel numerical methods were developed and evaluated for the dynamic simulation of membranes. For linear membranes, a mixed-form time-continuous Galerkin method was employed using trilinear space-time elements, and the entire space-time domain was discretized and solved simultaneously. For geometrically nonlinear membranes, the model incorporated two new schemes that were independently developed and evaluated. Time marching was performed using quintic Hermite polynomials uniquely determined by end-point jerk constraints. The single-step, implicit scheme was significantly more accurate than the most common Newmark schemes. For a simple harmonic oscillator, the scheme was found to be symplectic, frequency-preserving, and conditionally stable. Time step size was limited by accuracy requirements rather than stability. The spatial discretization scheme employed a staggered grid, grouping of nonlinear terms, and polygon shape functions in a strong-form point collocation formulation. Validation against existing experimental data showed the method to be accurate until hyperelastic effects dominate

Inverse dynamics of underactuated flexible mechanical systems governed by quasi-linear hyperbolic partial differential equations

Diese Arbeit befasst sich mit der inversen Dynamik unteraktuierter, flexibler, mechanischer Systeme, welche durch quasi-lineare hyperbolische partielle Differentialgleichungen beschrieben werden können. Diese Gleichungnen, sind zeitlich veränderlichen Dirichlet-Randbedingungen unterworfen, welche durch unbekannte, räumlich disjunkte, also nicht kollokierte Neumann-Randbedingungen erzwungen werden. Die zugrundeliegenden Gleichungen werden zunächst abstrakt hergeleitet, bevor verschiedene mechanische Systeme vorgestellt werden können, die mit der eingangs postulierten Formulierung übereinstimmen. Hierzu werden geometrisch exakte Theorien hergeleitet, welche in der Lage sind große Bewegungen schlanker Strukturen wie Seile und Balken, aber auch ganz allgemein, dreidimensionaler Festkörper zu beschreiben. In der Regel werden Anfangs-Randwertprobleme, die in der nichtlinearen Strukturdynamik auftreten, durch Anwendung einer sequentiellen Diskretisierung in Raum und Zeit gelöst. Diese Verfahren basieren für gewöhnlich auf einer räumlichen Diskretisierung mit finiten Elementen, gefolgt von einer geeigneten zeitlichen Diskretisierung, welche meist auf finiten Differenzen beruht. Ein kurzer Überblick über derartige sequentielle Integrationsverfahren für das vorliegende Anfangs-Randwertproblem wird zunächst anhand der direkten Formulierung des Problems gegeben werden. D.h. es wird zunächst das reine Neumann-Randproblem betrachtet, bevor anschließend ganz allgemein, verschiedene Möglichkeiten zur Einbindung etwaiger Dirichlet-Randbedingungen diskutiert werden. Darauf aufbauend wird das Problem der inversen Dynamik im Kontext räumlich diskreter mechanischer Systeme, welche rheonom-holonomen Servo-Bindungen unterliegen, eingeführt. Eine ausführliche Untersuchung dieser Art von gebundenen Systemen soll die grundlegenden Unterschiede zwischen Servo-Bindungen und klassischen Kontakt-Bindungen herausarbeiten. Die daraus resultierenden Folgen für die Entwicklung geeigneter numerisch stabiler Integrationsverfahren können dabei ebenfalls angesprochen werden, bevor zahlreich ausgewählte Beispiele vorgestellt werden können. Aufgrund der sehr eingeschränkten Anwendbarkeit der sequentiellen Lösung der inversen Dynamik in Raum und Zeit, wird eine eingehende Analyse des vorliegenden Anfangs-Randwertproblems unternommen. Vor allem durch die Freilegung der hyperbolischen Struktur der zugrundeliegenden partiellen Differentialgleichungen werden sich weitere Einblicke in das vorliegende Problem erhofft. Die Erforschung der daraus resultierenden Mechanismen der Wellenausbreitung in kontinuierlichen Strukturen öffnet die Tür zur Entwicklung numerisch stabiler Integrationsverfahren für die inverse Dynamik. So kann unter anderem eine Methode vorgestellt werden, die auf der Integration der partiellen Differentialgleichungen entlang charakteristischer Mannigfaltigkeiten beruht. Dies regt zu der Entwicklung neuartiger Galerkinverfahren an, die ebenfalls in dieser Arbeit vorgestellt werden können. Diese neu entwickelten Methoden können anschlie\ss end auf die Steuerung verschiedener mechanischer Systeme angewendet werden. Darüber hinaus können die neuartigen Integrationsverfahren auch auf flexible Mehrkörpersysteme übertragen werden. Angeführt seien hier beispielsweise die kooperative Steuerung eines an mehreren flexiblen Seilen aufgehängten starren Körpers oder die Steuerung des Endeffektors eines flexiblen mehrgliedrigen Schwenkarms. Ausgewählte numerische Beispiele verdeutlichen die Relevanz der hier vorgeschlagenen, in Raum und Zeit simultanen Integration des vorliegenden Anfangs-Randwertproblems

An explicit finite difference scheme for the Camassa-Holm equation

We put forward and analyze an explicit finite difference scheme for the Camassa-Holm shallow water equation that can handle general $H^1$ initial data and thus peakon-antipeakon interactions. Assuming a specified condition restricting the time step in terms of the spatial discretization parameter, we prove that the difference scheme converges strongly in $H^1$ towards a dissipative weak solution of Camassa-Holm equation.Comment: 45 pages, 6 figure

Applications of symmetries and conservation laws to the study of nonlinear elasticity equations

Mooney-Rivlin hyperelasticity equations are nonlinear coupled partial differential equations (PDEs) that are used to model various elastic materials. These models have been extended to account for fiber reinforced solids with applications in modeling biological materials. As such, it is important to obtain solutions to these physical systems. One approach is to study the admitted Lie symmetries of the PDE system, which allows one to seek invariant solutions by the invariant form method. Furthermore, knowledge of conservation laws for a PDE provides insight into conserved physical quantities, and can be used in the development of stable numerical methods. The current Thesis is dedicated to presenting the methodology of Lie symmetry and conservation law analysis, as well as applying it to fiber reinforced Mooney-Rivlin models. In particular, an outline of Lie symmetry and conservation law analysis is provided, and the partial differential equations describing the dynamics of a hyperelastic solid are presented. A detailed example of Lie symmetry and conservation law analysis is done for the PDE system describing plane strain in a Mooney-Rivlin solid. Lastly, Lie symmetries and conservation laws are studied in one and two dimensional models of fiber reinforced Mooney-Rivlin materials