Group implicit concurrent algorithms in nonlinear structural dynamics

During the 70's and 80's, considerable effort was devoted to developing efficient and reliable time stepping procedures for transient structural analysis. Mathematically, the equations governing this type of problems are generally stiff, i.e., they exhibit a wide spectrum in the linear range. The algorithms best suited to this type of applications are those which accurately integrate the low frequency content of the response without necessitating the resolution of the high frequency modes. This means that the algorithms must be unconditionally stable, which in turn rules out explicit integration. The most exciting possibility in the algorithms development area in recent years has been the advent of parallel computers with multiprocessing capabilities. So, this work is mainly concerned with the development of parallel algorithms in the area of structural dynamics. A primary objective is to devise unconditionally stable and accurate time stepping procedures which lend themselves to an efficient implementation in concurrent machines. Some features of the new computer architecture are summarized. A brief survey of current efforts in the area is presented. A new class of concurrent procedures, or Group Implicit algorithms is introduced and analyzed. The numerical simulation shows that GI algorithms hold considerable promise for application in coarse grain as well as medium grain parallel computers

Symplectic-energy-momentum preserving variational integrators

The purpose of this paper is to develop variational integrators for conservative mechanical systems that are symplectic and energy and momentum conserving. To do this, a space–time view of variational integrators is employed and time step adaptation is used to impose the constraint of conservation of energy. Criteria for the solvability of the time steps and some numerical examples are given

Asynchronous Variational Integrators

We describe a new class of asynchronous variational integrators (AVI) for nonlinear elastodynamics. The AVIs are distinguished by the following 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 variational principle. As a consequence of this variational structure, the algorithms conserve local momenta and a local discrete multisymplectic structure exactly. To guide the development of the discretizations, a spacetime multisymplectic formulation of elastodynamics is presented. The variational principle used incorporates both configuration and spacetime reference variations. This allows a unified treatment of all the conservation properties of the system.A discrete version of reference configuration is also considered, providing a natural definition of a discrete energy. The possibilities for discrete energy conservation are evaluated. 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 algorith

Variational integrators, the Newmark scheme, and dissipative systems

Variational methods are a class of symplectic-momentum integrators for ODEs. Using these schemes, it is shown that the classical Newmark algorithm is structure preserving in a non-obvious way, thus explaining the observed numerical behavior. Modifications to variational methods to include forcing and dissipation are also proposed, extending the advantages of structure preserving integrators to non-conservative systems

Nonsmooth Lagrangian mechanics and variational collision integrators

Variational techniques are used to analyze the problem of rigid-body dynamics with impacts. The theory of smooth Lagrangian mechanics is extended to a nonsmooth context appropriate for collisions, and it is shown in what sense the system is symplectic and satisfies a Noether-style momentum conservation theorem. Discretizations of this nonsmooth mechanics are developed by using the methodology of variational discrete mechanics. This leads to variational integrators which are symplectic-momentum preserving and are consistent with the jump conditions given in the continuous theory. Specific examples of these methods are tested numerically, and the long-time stable energy behavior typical of variational methods is demonstrated

Angiotensin II and Aldosterone Increase with Fasting in Breeding Adult Male Northern Elephant Seals (Mirounga angustirostris)

The renin‐angiotensin‐aldosterone system (RAAS) appears to contribute significantly to osmoregulation of fasting northern elephant seal (Mirounga angustirostris) pups; however, RAAS has not been characterized in fasting adult seals. Therefore, this study examined the contribution of RAAS to water turnover rates in fasting adult male northern elephant seals. Blood samples were obtained twice during their breeding fast at an interval of 6.5 wk, and water efflux rate was estimated by isotopic dilution during the same period. Serum electrolytes (Na+, K+, Cl−) and osmolality were unaltered between the two sampling periods, indicating ionic and osmotic homeostasis during the fast. Despite the lack of an increase in vasopressin, serum angiotensin II and aldosterone were increased and were significantly and positively correlated. Changes in aldosterone concentration and water efflux rate were significantly and negatively correlated, suggesting that the greater the increase in aldosterone, the smaller the loss of water. Adult male seals maintain ionic and osmotic homeostasis similar to that of fasting weaned pups, and this homeostasis appears to be mediated, at least in part, by RAAS, which probably contributes to increased water retention as well. The hormonal mechanisms by which northern elephant seals maintain water and electrolyte balance during fasting conditions appear to be similar regardless of age

Frictional Collisions Off Sharp Objects

This work develops robust contact algorithms capable of dealing with multibody nonsmooth contact geometries for which neither normals nor gap functions can be defined. Such situations arise in the early stage of fragmentation when a number of angular fragments undergo complex collision sequences before eventually scattering. Such situations precludes the application of most contact algorithms proposed to date