Singularity-free augmented Lagrangian algorithms for constrained multibody dynamics

After a general review of the methods currently available for the dynamics of constrained multibody systems in the context of numerical efficiency and ability to solve the differential equations of motion in singular positions, we examine the acceleration based augmented Lagrangian formulations, and propose a new one for holonomic and non-holonomic systems that is based on the canonical equations of Hamilton. This new one proves to be more stable and accurate that the acceleration based counterpart under repetitive singular positions. The proposed algorithms are numerically efficient, can use standard conditionally stable numerical integrators and do not fail in singular positions, as the classical formulations do. The reason for the numerical efficiency and better behavior under singularities relies on the fact that the leading matrix of the resultant system of ODEs is sparse, symmetric, positive definite, and its rank is independent of that of the Jacobian of the constraint equations. The latter fact makes the proposed method particularly suitable for singular configurations.     

Resolving Controversies in the Application of the Method of Multiple Scales and the Generalized Method of Averaging

I compare application of the method of multiple scales with reconstitution and the generalized method of averaging for determining higher-order approximations of three single-degree-of-freedom systems and a two-degree-of-freedom system. Three implementations of the method of multiple scales are considered, namely, application of the method to the system equations expressed as second-order equations, as first-order equations, and in complex-variable form. I show that all of these methods produce the same modulation equations.I address the problem of determining higher-order approximate solutions of the Duffing equation in the case of primary resonance. I show that the conclusions of Rahman and Burton that the method of multiple scales, the generalized method of averaging, and Lie series and transforms might lead to incorrect results, in that spurious solutions occur and the obtained frequency–response curves bear little resemblance to the actual response, is the result of their using parameter values for which the neglected terms are the same order as the retained terms. I show also that spurious solutions cannot be avoided, in general, in any consistent expansion and their presence does not constitute a limitation of the methods. In particular, I show that, for the Duffing equation, the second-order frequency–response equation does not possess spurious solutions for the case of hardening nonlinearity, but possesses spurious solutions for the case of softening nonlinearity. For sufficiently small nonlinearity, the spurious solutions are far removed from the actual response. But as the strength of the nonlinearity increases, these solutions move closer to the backbone and eventually distort it. This is not a drawback of the perturbation methods but an indication of an application of the analysis for parameter values outside the range of validity of the expansion.Also, I address the problem of obtaining non-Hamiltonian modulation equations in the application of the method of multiple scales to multi-degree-of-freedom Hamiltonian systems written as second-order equations in time and how this problem can be overcome by attacking the state-space form of the governing equations. Moreover, I show that application of a variation of the method of Rahman and Burton to multi-degree-of-freedom systems leads to results that do not agree with those obtained with the generalized method of averaging.Contributed by Prof. R.A. Ibrahim.     

Creep—Slip Capture as a Possible Source of Squeal During Decelerated Sliding

Friction-induced vibration of a two-degree-of-freedom mass-damper-spring system interacting with a decelerating rigid strip is investigated. The friction law is approximated by an analytical function to facilitate the analyses and numerical integrations. It is shown that, after a quasi-harmonic transient period, accompanied by viscous energy dissipation, a short period of intensive ‘creep-slip’ vibration occurs, which generates a series of ‘micro-impacts’ on the strip. Because of the impulsive character of such kind of loading, its Fourier spectrum is rich and quite broadband. Using an averaging technique, the ‘normal form’ equations of motion show that the out-of-phase vibration mode absorbs more energy from the decelerating strip when its natural frequency satisfies certain resonance conditions. The study is then applied to an automotive disc brake model to gain useful insight into the generation of squeal. It is shown that the out-of-phase creep-slip vibration (in the longitudinal direction) of the brake pads generates an impulsive bending moment on the decelerating strip (disc rotor). This impulsive load may be considered as a possible source for brake squeal. The technique developed in this paper may be extended to other ‘squealing systems’ including models for geophysical faults (earthquakes).     

Cargo Pendulation Reduction of Ship-Mounted Cranes

Ship-mounted cranes are used to transfer cargo from large container ships to smaller lighters when deep-water ports are not available. The wave-induced motion of the crane ship can produce large pendulations of the cargo being hoisted and cause operations to be suspended. In this work, we show that it is possible to reduce these pendulations significantly by controlling the slew and luff angles of the boom. Such a control can be achieved with the heavy equipment that is already part of the crane so that retrofitting existing cranes would require a small effort. Moreover, the control is superimposed on the commands of the operator transparently. The successful control strategy is based on delayed feedback of the angles of the cargo-hoisting cable in and out of the plane of the boom and crane tower. Its effectiveness is demonstrated in a fully nonlinear three-dimensional computer simulation and in an experiment with a 1/24th-scale model of a T-ACS (The Auxiliary Crane Ship) crane mounted on a platform moving with three degrees of freedom. The results demonstrate that the pendulations can be significantly reduced, and therefore the range of sea conditions in which cargo-transfer operations can take place can be greatly expanded.     

Control of a Nonlinear System Using the Saturation Phenomenon

In this paper, a theoretical investigation of nonlinear vibrations of a 2 degrees of freedom system when subjected to saturation is studied. The method has been especially applied to a system that consists of a DC motor with a nonlinear controller and a harmonic forcing voltage. Approximate solutions are sought using the method of multiple scales. It is shown that the closed-loop system exhibits different response regimes. The nature and stability of these regimes are studied and the stability boundaries are obtained. The effects of the initial conditions on the response of the system have also been investigated. Furthermore, the second-order solution is presented and the corresponding results are compared with those of the first-order solution. It is shown that by increasing the amplitude of the excitation voltage, the higher-order term in the solution becomes significant and causes a drift in the response. In order to verify the obtained theoretical results, they are compared with the corresponding numerical results. Good agreement between the two sets of results is observed.     

Simpler problems for elastic Cosserat surfaces

For elastic Cosserat surfaces, we discuss an attempt to characterize cases where the equations can be reduced to a system of ordinary differential equations.

---

Résumé Pour les surfaces élastiques de Cosserat, on decrira les different cas où les equations peuvent être rameneés a un system d'equations différentielles ordinaires.

     

Necking of biaxially stretched elastic-plastic circular plates

The necking of an elastic-plastic circular plate under uniform radial tensile loading is investigated both within the framework of the three-dimensional theory and within the context of the plane-stress approximation. Attention is restricted to axisymmetric deformations of the plate. The material behavior is described by two different constitutive laws. One is a finite-strain version of the simplest flow-theory of plasticity and the other is a finite-strain generalization of the simplest deformationtheory, which is employed as a simple model of a solid with a vertex on its yield surface. For an initially uniform plate made of an incompressible material, bifurcation from the uniformly stretched state is studied analytically. The regimes of stress and moduli where the governing axisymmetric three-dimensional equations are elliptic, parabolic or hyperbolic are identified. The plane-stress local-necking mode emerges as the appropriate limiting mode from the bifurcation modes available in the elliptic regime. In the elliptic regime, the main qualitative features of the bifurcation behavior are revealed by the plane-stress analysis, although three-dimensional effects delay the onset of necking somewhat. For the deformation theory employed here, the first bifurcation modes are encountered in the parabolic regime if the hardening-rate is sufficiently high. These bifurcations are not revealed by a plane-stress analysis. For a plate with an initial inhomogeneity, the growth of an imperfection is studied by a perturbation method, by a plane-stress analysis of localized necking, and by numerical computations within the framework of the three-dimensional theory. When bifurcation of the corresponding perfect plate takes place in the elliptic regime, the finite element results show that the plane-stress analysis gives reasonably good agreement with the numerical results. When bifurcation of the corresponding perfect plate first occurs in the parabolic regime, then a bifurcation of the imperfect plate is encountered, that is, the finite element stiffness matrix ceases to be positive definite.     

The geometrical stability of non-linear normal modes in two degree of freedom systems

The geometrical stability of the non-linear normal mode vibrations of a class of two degree of freedom dynamical systems is studied by utilizing the definitions and analysis of Synge's “Geometry of Dynamics.”It is shown that instabilities can occur for small amplitudes of vibration only if (a) one of the associated linear normal modes possesses a frequency which is nearly a multiple of the frequency of the other linear normal mode, or (b) the frequency of one linear normal mode is nearly zero.