Abstract:A temporal finite element method based on a mised form of the Hamiltonian n-eal; principle is developed for dynamics and optimal control problems. The mixed form of Hamilton's weak principle contains both displacements and momenta as primary variables that are expanded in terms of nodal values and simple polynomial shape functions. Unlike other forms of Hamilton's principle, however, time derivatives of the nionient a and displacements do not appear therein; instead, only the virtual momenta and virtual displa… Show more
“…DFET was initially proposed in [8] and uses finite elements in time on spectral bases to transcribe the differential equations into a set of algebraic equations. Finite Elements in Time for the indirect solution of optimal control problems were initially proposed by Hodges et al in [9], and during the late 1990s evolved to the discontinuous version. As pointed out by Bottasso et al in [10], FET for the forward integration of ordinary differential equations are equivalent to some classes of implicit RungeKutta integration schemes, can be extended to arbitrary highorder, are very robust and allow full h-p adaptivity.…”
“…DFET was initially proposed in [8] and uses finite elements in time on spectral bases to transcribe the differential equations into a set of algebraic equations. Finite Elements in Time for the indirect solution of optimal control problems were initially proposed by Hodges et al in [9], and during the late 1990s evolved to the discontinuous version. As pointed out by Bottasso et al in [10], FET for the forward integration of ordinary differential equations are equivalent to some classes of implicit RungeKutta integration schemes, can be extended to arbitrary highorder, are very robust and allow full h-p adaptivity.…”
“…Dynamics community have also considered this type of approach in regard to various applications such as structural dynamics, nonlinear vibrations, inverse dynamics of flexible robots, gate analysis, and optimal control problem for converting problems from infinite dimension to finite dimension. Some of the related works can be found in [23]- [29]. These prior time-parameterization methods often produced highly nonlinear and coupled sets of equations, thus limiting their tractability with complex systems.…”
Section: State-time Vs Traditional State Dynamic Formulationsmentioning
This paper presents a new methodology demonstrating the feasibility and advantages of a state-time formulation for dynamic simulation of complex multibody systems which shows potential advantages for exploiting massively parallel computing resources. This formulation allows time to be discretized and parameterized so that it can be treated as a variable in a manner similar to the system state variables. As a consequence of such a state-time discretization scheme, the system of governing equations yields to a set of loosely coupled linear-quadratic algebraic equations that is well-suited in structure for some families of nonlinear algebraic equations solvers. The goal of this work is to develop efficient multibody dynamics algorithm that is extremely scalable and better able to fully exploit anticipated immensely parallel computing machines (tera flop, pecta flop and beyond) made available to it.
KeywordMultibody dynamics, State-time formulation, Scalable algorithm, Parallel computing
Nomenclature
B: Typical body B with its mass center B * F kA : k th applied force acting on the body B f mC : m th unknown constraint force acting on the body B r kA : k th position vector from B * to application point of F kA r B * jm : m th position vector from B * to application point (joint) of f mC
“…Riff and Baruch [9] [15,29] demonstrated that a mixed finite element formulation (where generalized coordinates and their momenta appear as independent unknowns) yields solutions with unconditional stability without having to use reduced element quadrature. The current research is also based on a mixed finite element formulation.…”
Section: Hamilton's Lawmentioning
confidence: 99%
“…Because there are no derivatives, only simple shape functions are necessary to satisfy C 0 continuity between elements. Higher order elements (p elements) could also be developed; however, for general nonlinear problems the use of crude shape functions is more efficient in that it allows element quadrature to be performed by inspection [29]. Independent variables need only be piecewise continuous, while variational quantities must be continuous and piecewise differentiable.…”
Section: Hamilton's Weak Principlementioning
confidence: 99%
“…A time marching procedure was used in [45] to demonstrate accuracy for some simple problems. Hodges [29] applies HWP in mixed variational form to optimal control problems and nonlinear initial value problems. He also uses the technique to develop exact intrinsic equations for a moving beam [26], a derivation which is used extensively in this research.…”
Hamilton's Law is derived in weak form for slender beams with closed cross sections. The result is discretized with mixed space-time finite elements to yield a system of nonlinear, algebraic equations. An algorithm is proposed for solving these equations using unconstrained optimization techniques, obtaining steady-state and time accurate solutions for problems of structural dynamics. This technique provides accurate solutions for nonlinear static and steady-state problems including the cantilevered elastica and flatwise rotation of beams. Modal analysis of beams and rods is investigated to accurately determine fundamental frequencies of vibration, and the simulation of simple maneuvers is demonstrated.iv
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.