objective energy momentum conserving integration for the constrained dynamics of geometrically exact beams in this paper the results in s leyendecker p betsch p steinmann energy conserving integration of constrained hamiltonian systems a comparison of approaches comput mech 33 2004 174 185 are extended to geometrically exact beams the finite element formulation for nonlinear beams in terms of directors providing a framework for the objective description of their dynamics is considered geometrically exact beams are analysed as hamiltonian systems subject to holonomic constraints with a hamiltonian being invariant under the action of so 3 the reparametrisation of the hamiltonian in terms of the invariants of so 3 is perfectly suited for a temporal discretisation which leads to energy momentum conserving integration in this connection the influence of alternative procedures the lagrange multiplier method the penalty method and the augmented lagrange method for the treatment of the constraints is investigated for the example of a beam with concentrated masses 2005 elsevier b v all rights reserved

on the parametrization of finite rotations in computational mechanics a classification of concepts with application to smooth shells this paper concerns the computational treatment of large rotations with application to the finite element discretisation of smooth shells formulations based on various rotational parametrizations are reviewed and classified with respect to their update structure category ii is based on total rotational degrees of freedom leading to an additive update structure whereas category i relies on linearized rotational degrees of freedom leading to a multiplicative update structure based on this classification new formulations are developed among them are the rodrigues formula applied within category ii with additive update structure for the two components of the rotation vector and two successive elementary rotations in the sense of euler angles applied within category i yielding a singularity free formulation the numerical examples confirm that every rotational parametrization when applied within category ii leads to singularities whereas the application within category i enables the calculation of overall rotations unrestricted in size without any singularity

frame indifferent beam finite elements based upon the geometrically exact beam theory this paper describes a new beam finite element formulation based upon the geometrically exact beam theory in contrast to many previously proposed beam finite element formulations the present discretization approach retains the frame indifference or objectivity of the underlying beam theory the space interpolation of rotational degrees of freedom is circumvented by the introduction of a reparameterization of the weak form corresponding to the equations of motion of the geometrically exact beam theory copyright 2002 john wiley sons ltd

velocity based approach in non linear dynamics of three dimensional beams with enforced kinematic compatibility abstract in the paper we present a new finite element formulation for the dynamic analysis of geometrically exact three dimensional beams we limit our studies to implicit time integration schemes and possible approaches for increasing their robustness and numerical stability in contrast to standard displacement rotation based approach we present here a spatial and temporal discretization based on velocities and angular velocities to describe the rotational degrees of freedom quaternions are used the time integration scheme and the governing equations of the three dimensional beam are modified accordingly in the numerical implementation the galerkin type discretization is employed to obtain the finite element formulation of the problem the result of our studies is simple but accurate efficient and robust numerical model

discrete time lagrangian mechanics on lie groups with an application to the lagrange top we develop the theory of discrete time lagrangian mechanics on lie groups originated in the work of veselov and moser and the theory of lagrangian reduction in the discrete time setting the results thus obtained are applied to the investigation of an integrable time discretization of a famous integrable system of classical mechanics the lagrange top we recall the derivation of the euler poinsot equations of motion both in the frame moving with the body and in the rest frame the latter ones being less widely known we find a discrete time lagrange function turning into the known continuous time lagrangian in the continuous limit and elaborate both descriptions of the resulting discrete time system namely in the body frame and in the rest frame this system naturally inherits poisson properties of the continuous time system the integrals of motion being deformed the discrete time lax representations are also found kirchhoff s kinetic analogy between elastic curves and motions of the lagrange top is also generalised to the discrete context

galerkin lie group variational integrators based on unit quaternion interpolation abstract lie group variational integrators of arbitrary order are developed using the galerkin method based on unit quaternion interpolation to our knowledge quaternions have not been used before for this purpose though they allow a very simple and efficient way to perform the interpolation the resulting integrators are symplectic and structure preserving in the sense that certain symmetries in the lagrangian of the mechanical system are carried over to the discrete setting which leads to the preservation of the corresponding momentum maps the integrators furthermore exhibit a very good long time energy behavior i e energy is neither dissipated nor gained artificially at the same time the lie group structure is preserved by carefully defining the variations the interpolation method and by solving the non linear system of equations directly on the manifold rather than constraining it in a surrounding space using lagrange multipliers as a consequence we are able to show that lie group variational integrators based on the special orthogonal group are equivalent to the variational integrators for constrained systems using the discrete null space method employed e g in dmocc discrete mechanics and optimal control of constrained systems we show new numerical results on the convergence rates which are substantially higher than the known theoretical bounds and on the relation between accuracy and computational cost

multisymplectic geometry variational integrators and nonlinear pdes this paper presents a geometric variational approach to continuous and discrete mechanics and field theories using multisymplectic geometry we show that the existence of the fundamental geometric structures as well as their preservation along solutions can be obtained directly from the variational principle in particular we prove that a unique multisymplectic structure is obtained by taking the derivative of an action function and use this structure to prove covariant generalizations of conservation of symplecticity and noether s theorem natural discretization schemes for pdes which have these important preservation properties then follow by choosing a discrete action functional in the case of mechanics we recover the variational symplectic integrators of veselov type while for pdes we obtain covariant spacetime integrators which conserve the corresponding discrete multisymplectic form as well as the discrete momentum mappings corresponding to symmetries we show that the usual notion of symplecticity along an infinite dimensional space of fields can be naturally obtained by making a spacetime split all of the aspects of our method are demonstrated with a nonlinear sine gordon equation including computational results and a comparison with other discretization schemes

discrete versions of some classical integrable systems and factorization of matrix polynomials discrete versions of several classical integrable systems are investigated such as a discrete analogue of the higher dimensional force free spinning top euler arnold equations the heisenberg chain with classical spins and a new discrete system on the stiefel manifold the integrability is shown with the help of a lax pair representation which is found via a factorization of certain matrix polynomials the complete description of the dynamics is given in terms of abelian functions the flow becomes linear on a prym variety corresponding to a spectral curve the approach is also applied to the billiard problem in the interior of ann dimensional ellipsoid

variational integrators and the newmark algorithm for conservative and dissipative mechanical systems the purpose of this work is twofold first we demonstrate analytically that the classical newmark family as well as related integration algorithms are variational in the sense of the veselov formulation of discrete mechanics such variational algorithms are well known to be symplectic and momentum preserving and to often have excellent global energy behavior this analytical result is veried through numerical examples and is believed to be one of the primary reasons that this class of algorithms performs so well second we develop algorithms for mechanical systems with forcing and in particular for dissipative systems in this case we develop integrators that are based on a discretization of the lagrange d alembert principle as well as on a variational formulation of dissipation it is demonstrated that these types of structured integrators have good numerical behavior in terms of obtaining the correct amounts by which the energy changes over the integration run

finite rotations in dynamics of beams and implicit time stepping schemes we examine theoretical and computational aspects of three dimensional finite rotations pertinent to the dynamics of beams the model problem chosen for consideration is the reissner beam theory capable of modelling finite strains and finite rotations in geometrically exact manner special emphasis is placed on clarifying the geometry aspects finite rotation updates and the associated linearization procedure pertaining to different choices of rotation parameters the latter is shown to play an important role in constructing the optimal implementation of a time stepping scheme 1998 john wiley sons ltd

discrete variational lie group formulation of geometrically exact beam dynamics the goal of this paper is to derive a structure preserving integrator for geometrically exact beam dynamics by using a lie group variational integrator both spatial and temporal discretization are implemented in a geometry preserving manner the resulting scheme preserves both the discrete momentum maps and symplectic structures and exhibits almost perfect energy conservation comparisons with existing numerical schemes are provided and the convergence behavior is analyzed numerically

computational aspects of vector like parametrization of three dimensional finite rotations theoretical and computational aspects of vector like parametrization of three dimensional finite rotations which uses only three rotation parameters are examined in detail in this work the relationship of the proposed parametrization with the intrinsic representation of finite rotations via an orthogonal matrix is clearly identified careful considerations of the consistent linearization procedure pertinent to the proposed parametrization of finite rotations are presented for the chosen model problem of reissner s non linear beam theory pertaining details of numerical implementation are discussed for the simplest choice of the finite element interpolations for a 2 node three dimensional beam element a number of numerical simulations in three dimensional finite rotation analysis are presented in order to illustrate the proposed approach

geometrically exact 3d beam theory implementation of a strain invariant finite element for statics and dynamics geometrically exact 3d beam theory has been used as a basis for development of a variety of finite element formulations it has recently become apparent that the important requirement of objectivity of adopted strain measures although provided by the theory itself does not automatically extend to a finite element formulation in this paper we present a new finite element formulation of the geometrically exact 3d beam theory specifically designed to preserve the objectivity of the adopted strain measures in order to do so the current local rotations are interpolated in a manner similar to that adopted in co rotational approaches however no approximations typical for co rotational approaches are introduced into the procedure so in contrast to co rotational formulations the present formulation fully preserves the geometric exactness of the theory a range of numerical examples serves to illustrate the problem and to assess the formulation

interpolation of rotational variables in nonlinear dynamics of 3d beams the formulation of dynamic procedures for three dimensional 3 d beams requires extensive use of the algebra pertaining to the non linear character of the rotation group in space the corresponding extraction procedure to obtain the rotations that span a time step has certain limitations which can have a detrimental eect on the overall stability of a time integration scheme the paper describes two algorithms for the dynamics of 3 d beams which dier in their manifestation of the above limitation the rst has already been described in the literature and involves the interpolation of iterative rotations while an alternative formulation which eliminates the above eect by design requires interpolation of incremental rotations theoretical arguments are backed by numerical results similarities between the conventional and new formulation are pointed out and are shown to be big enough to enable easy transformation of the conventional formulation into the new one 1998 john wiley sons ltd

objectivity of strain measures in the geometrically exact three dimensional beam theory and its finite element implementation the paper discusses the issue of discretization of the strainconfiguration relationships in the geometrically exact theory of threedimensional 3d beams which has been at the heart of most recent

entanglement criteria for dicke states dicke states are a family of multi qubit quantum states with interesting entanglement properties and they have been observed in many experiments we construct entanglement witnesses for detecting genuine multiparticle entanglement in the vicinity of these states we use the approach of positive partial transpose mixtures to derive the conditions analytically for nearly all cases our criteria are stronger than all conditions previously known

the hamiltonian structure of nonlinear elasticity the material and convective representations of solids rods and plates it is our belief that a thorough understanding of the mathematical underpinnings of elasticity is crucial to its analytical and numerical implementation for example in the analysis of rotating structures the coupling of the equations for geometrically inexact models obtained by linearization or other approximations with those for rotating rigid bodies can easily lead to misleading artificial softening effects that can significantly alter numerical results see simo and vuquoc 1986c especially equations 3 and 5 in this paper we consider fully nonlinear geometrically exact models for rods plates and shells which take into account shear and torsion as well as the usual bending effects in traditional rod and plate models these models can be obtained either from the three dimensional theory by a systematic use of projection methods see e g antman 1972 and naghdi 1972 or by a direct approach within the context of cosserat continuum remarkably the two approaches lead to essentially the same form of the governing field equations in the present context we have chosen as a model problem a particular rod model which may be regarded as an extension of the classical kirchhoff love model see love 1944 to include shear deformations as in

mimetic finite difference methods for hamiltonian wave equations in 2d abstract in this paper we consider the numerical solution of the hamiltonian wave equation in two spatial dimensions we construct a two step procedure in which we first discretize the space by the mimetic finite difference mfd method and then we employ a standard symplectic scheme to integrate the semi discrete hamiltonian system derived the main characteristic of the mfd methods when applied to stationary problems is to mimic important properties of the continuous system this approach yields a full numerical procedure suitable to integrate hamiltonian problems a complete theoretical analysis of the method and some numerical simulations are developed in the paper

global optimization with polynomials and the problem of moments we consider the problem of finding the unconstrained global minimum of a real valued polynomial p x mathbb r n to mathbb r as well as the global minimum of p x in a compact set k defined by polynomial inequalities it is shown that this problem reduces to solving an often finite sequence of convex linear matrix inequality lmi problems a notion of karush kuhn tucker polynomials is introduced in a global optimality condition some illustrative examples are provided

a discrete mechanics approach to the cosserat rod theory part 1 static equilibria a theory of discrete cosserat rods is formulated in the language of discrete lagrangian mechanics by exploiting kirchhoff s kinetic analogy the potential energy density of a rod is a function on the tangent bundle of the configuration manifold and thus formally corresponds to the lagrangian function of a dynamical system the equilibrium equations are derived from a variational principle using a formulation that involves null space matrices in this formulation no lagrange multipliers are necessary to enforce orthonormality of the directors noether s theorem relates first integrals of the equilibrium equations to lie group actions on the configuration bundle so called symmetries the symmetries relevant for rod mechanics are frame indifference isotropy and uniformity we show that a completely analogous and self contained theory of discrete rods can be formulated in which the arc length is a discrete variable ab initio in this formulation the potential energy density is defined directly on pairs of points along the arc length of the rod in analogy to veselov s discrete reformulation of lagrangian mechanics a discrete version of noether s theorem then identifies exact first integrals of the discrete equilibrium equations these exact conservation properties confer the discrete solutions accuracy and robustness as demonstrated by selected examples of application

refractors in anisotropic media associated with norms abstract we show existence of interfaces between two anisotropic materials so that light is refracted in accordance with a given pattern of energy to do this we formulate a vector snell law for anisotropic media when the wave fronts are given by norms for which the corresponding unit spheres are strictly convex

variational time integrators the purpose of this paper is to review and further develop the subject of variational integration algorithms as it applies to mechanical systems of engineering interest in particular the conservation properties of both synchronous and asynchronous variational integrators avis are discussed in detail we present selected numerical examples which demonstrate the excellent accuracy conservation and convergence characteristics of avis in these tests avis are found to result in substantial speed ups at equal accuracy relative to explicit newmark a mathematical proof of convergence of the avis is also presented in this paper finally we develop the subject of horizontal variations and configurational forces in discrete dynamics this theory leads to exact path independent characterizations of the configurational forces acting on discrete systems notable examples are the configurational forces acting on material nodes in a finite element discretisation and the j integral at the tip of a crack in a finite element mesh

asynchronous variational integrators we describe a class of asynchronous variational integrators avi for nonlinear elastodynamics the avis are characterized by the following distinguishing attributes i the algorithms permit the selection of independent time steps in each element and the local time steps need not bear an integral relation to each other ii the algorithms derive from a spacetime form of a discrete version of hamilton s principle as a consequence of this variational structure the algorithms conserve local energy and momenta exactly subject to solvability of the local time steps numerical tests reveal that even when local energy balance is not enforced exactly the global and local energy behavior of the avis is quite remarkable a property which can probably be traced to the symplectic nature of the algorithm

an objective finite element approximation of the kinematics of geometrically exact rods and its use in the formulation of an energy momentum conserving scheme in dynamics we present in this paper a new finite element formulation of geometrically exact rod models in the three dimensional dynamic elastic range the proposed formulation leads to an objective or frame indifferent under superposed rigid body motions approximation of the strain measures of the rod involving finite rotations of the director frame in contrast with some existing formulations this goal is accomplished through a direct finite element interpolation of the director fields defining the motion of the rod s cross section furthermore the proposed framework allows the development of time stepping algorithms that preserve the conservation laws of the underlying continuum hamiltonian system the conservation laws of linear and angular momenta are inherited by construction leading to an improved approximation of the rod s dynamics several numerical simulations are presented illustrating these properties copyright 2002 john wiley sons ltd

bell inequalities for graph states we investigate the nonlocal properties of graph states to this aim we derive a family of bell inequalities which require three measurement settings for each party and are maximally violated by graph states in turn for each graph state there is an inequality maximally violated only by that state we show that for certain types of graph states the violation of these inequalities increases exponentially with the number of qubits we also discuss connections to other entanglement properties such as the positivity of the partial transpose or the geometric measure of entanglement

the discrete null space method for the energy consistent integration of constrained mechanical systems part i holonomic constraints the present work deals with energy consistent time stepping schemes for finite dimensional mechanical systems with holonomic constraints the proposed procedure is essentially based upon the following steps firstly the index three differential algebraic equations corresponding to the constrained mechanical system are directly discretized secondly the discrete lagrange multipliers are eliminated by using a discrete null space matrix in many cases it is feasible to further reduce the number of unknowns by employing specific reparametrizations the proposed method entails a number of advantageous features such as size reduction and improved conditioning of the resulting system of algebraic equations it is shown that the newly developed method is well suited for both open loop and closed loop multibody systems