Analysis of periodic-quasiperiodic nonlinear systems via Lyapunov-Floquet transformation and normal forms

In this paper a general technique for the analysis of nonlinear dynamical systems with periodic-quasiperiodic coefficients is developed. For such systems the coefficients of the linear terms are periodic with frequency ω while the coefficients of the nonlinear terms contain frequencies that are incommensurate with ω. No restrictions are placed on the size of the periodic terms appearing in the linear part of system equation. Application of Lyapunov-Floquet transformation produces a dynamically equivalent system in which the linear part is time-invariant and the time varying coefficients of the nonlinear terms are quasiperiodic. Then a series of quasiperiodic near-identity transformations are applied to reduce the system equation to a normal form. In the process a quasiperiodic homological equation and the corresponding ‘solvability condition’ are obtained. Various resonance conditions are discussed and examples are included to show practical significance of the method. Results obtained from the quasiperiodic time-dependent normal form theory are compared with the numerical solutions. A close agreement is found.     

On the effects of non-linear elements in the reliability-based optimal design of stochastic dynamical systems

The aim of the present paper is to study the effects of non-linear devices on the reliability-based optimal design of structural systems subject to stochastic excitation. One-dimensional hysteretic devices are used for modelling the non-linear system behavior while non-stationary filtered white noise processes are utilized to represent the stochastic excitation. The reliability-based optimization problem is formulated as the minimization of the expected cost of the structure for a specified failure probability. Failure is assumed to occur when any one of the output states of interest exceeds in magnitude some specified threshold level within a given time duration. Failure probabilities are approximated locally in terms of the design variables during the optimization process in a parallel computing environment. The approximations are based on a local interpolation scheme and on an efficient simulation technique. Specifically, a subset simulation scheme is adopted and integrated into the proposed optimization process. The local approximations are then used to define a series of explicit approximate optimization problems. A sensitivity analysis is performed at the final design in order to evaluate its robustness with respect to design and system parameters. Numerical examples are presented in order to illustrate the effects of hysteretic devices on the design of two structural systems subject to earthquake excitation. The obtained results indicate that the non-linear devices have a significant effect on the reliability and global performance of the structural systems.     

非线性水波Hamilton系统理论与应用研究进展

概述了辛几何理论与辛算法在Hamilton力学中的应用，综述非线性水波的Hamilton理论研究进展．阐述非线性水波Hamilton变分原理与方法的优越性与局限性，探讨KdV方程和BBM方程的Hamilton描述、对称性与守恒律，提出非线性水波Hamilton描述研究中有待进一步研究的问题和解法设想．     

An Approximate Analysis of Dry-Friction-Induced Stick-Slip Vibrations by a Smoothing Procedure

This paper deals with a systematic procedure to find both stable and unstable periodic stick-slip vibrations of autonomous dynamic systems with dry friction. In this procedure, the discontinuous friction forces are approximated by smooth functions. Using the simple shooting method with a stiff-ODE solver, in combination with a path following algorithm, branches of periodic solutions are computed for a changing design variable. For testing purposes, both 1 and 2-DOF autonomous block-on-belt models and a 1-DOF autonomous drill string model from literature are investigated. Comparison of the results shows that the smoothing procedure accurately describes the behavior of the discontinuous systems. The proposed procedure can also easily be applied to more complex MDOF models, as well as to nonautonomous dynamic systems.     

Multi-physics dynamics of a mechanical oscillator coupled to an electro-magnetic circuit

The dynamics of a non-linear electro-magneto-mechanical coupled system is addressed. The non-linear behavior arises from the involved coupling quadratic non-linearities and it is explored by relying on both analytical and numerical tools. When the linear frequency of the circuit is larger than that of the mechanical oscillator, the dynamics exhibits slow and fast time scales. Therefore the mechanical oscillator forced (actuated) via harmonic voltage excitation of the electric circuit is analyzed; when the forcing frequency is close to that of the mechanical oscillator, the long term damped dynamics evolves in a purely slow timescale with no interaction with the fast time scale. We show this by assuming the existence of a slow invariant manifold (SIM), computing it analytically, and verifying its existence via numerical experiments on both full- and reduced-order systems. In specific regions of the space of forcing parameters, the SIM is a complicated geometric object as it undergoes folding giving rise to hysteresis mechanisms which create a pronounced non-linear resonance phenomenon. Eventually, the roles played by the electro-magnetic and mechanical components in the resulting complex response, encompassing bifurcations as well as possible transitions from regular to chaotic motion, are highlighted by means of Poincaré sections.     

Implicit frictional-contact model for soft particle systems

We introduce a novel numerical approach for the simulation of soft particles interacting via frictional contacts. This approach is based on an implicit formulation of the Material Point Method, allowing for large particle deformations, combined with the Contact Dynamics method for the treatment of unilateral frictional contacts between particles. This approach is both precise due to the treatment of contacts with no regularization and artificial damping parameters, and robust due to implicit time integration of both bulk degrees of freedom and relative contact velocities at the nodes representing the contact points. By construction, our algorithm is capable of handling arbitrary particle shapes and deformations. We illustrate this approach by two simple 2D examples: a Hertz contact and a rolling particle on an inclined plane. We also investigate the compaction of a packing of circular particles up to a solid fraction well above the jamming limit of hard particles. We find that, for the same level of deformation, the solid fraction in a packing of frictional particles is above that of a packing of frictionless particles as a result of larger particle shape change.     

Diffusion in isotropic and anisotropic porous systems: Three-dimensional calculations

Effective diffusion coefficients were calculated numerically for three-dimensional unit cells representative of different unconsolidated porous media. These numerical results were compared with the experimental results of Kim for packed beds of glass spheres, mica particles, and an artificial porous medium composed of mylar disks. These three-dimensional numerical results confirm that the porosity is the essential parameter for the determination of the effective diffusion coefficient in the case of unconsolidated isotropic systems. In the case of anisotropic systems, better agreement is obtained between numerical predictions and actual data when the unit cell is three-dimensional rather than twodimensional. This emphasizes the fact that three-dimensional unit cells feature more realistic geometrical properties which are needed to accurately describe anisotropic systems.     

随机激励的耗散的哈密尔顿系统的平稳解

本文首先为一般的随机激励的耗散的哈密尔顿系统得到精确的平稳解,然后在此基础上为类似而更为一般的系统发展了等效非线性系统法     

Hyperbolicity singularities in Rarefaction Waves

For mixed-type systems of conservation laws, rarefaction waves may contain states at the boundary of the elliptic region, where two characteristic speeds coincide, and the Lax family of the wave changes. Such contiguous rarefaction waves form a single fan with a continuous profile. Different pairs of families may appear in such rarefactions, giving rise to novel Riemann solution structures. We study the structure of such rarefaction waves near regular and exceptional points of the elliptic boundary and describe their effect on Riemann solutions.