LARGEST LYAPUNOV EXPONENT AND ALMOST CERTAIN STABILITY ANALYSIS OF SLENDER BEAMS UNDER A LARGE LINEAR MOTION OF BASEMENT SUBJECT TO NARROWBAND PARAMETRIC EXCITATION

The first order approximate solutions of a set of non-liner differential equations, which is established by using Kane's method and governs the planar motion of beams under a large linear motion of basement, are systematically derived via the method of multiple scales. The non-linear dynamic behaviors of a simply supported beam subject to narrowband random parametric excitation, in which either the principal parametric resonance of its first mode or a combination parametric resonance of the additive type of its first two modes with or without 3:1 internal resonance between the first two modes is taken into consideration, are analyzed in detail. The largest Lyapunov exponent is numerically obtained to determine the almost certain stability or instability of the trivial response of the system and the validity of the stability is verified by direct numerical integration of the equation of motion of the system.     

APPROXIMATE ANALYTICAL SOLUTIONS FOR PRIMARY CHATTER IN THE NON-LINEAR METAL CUTTING MODEL

This paper considers an accepted model of the metal cutting process dynamics in the context of an approximate analysis of the resulting non-linear differential equations of motion. The process model is based upon the established mechanics of orthogonal cutting and results in a pair of non-linear ordinary differential equations which are then restated in a form suitable for approximate analytical solution. The chosen solution technique is the perturbation method of multiple time scales and approximate closed-form solutions are generated for the most important non-resonant case. Numerical data are then substituted into the analytical solutions and key results are obtained and presented. Some comparisons between the exact numerical calculations for the forces involved and their reduced and simplified analytical counterparts are given. It is shown that there is almost no discernible difference between the two thus confirming the validity of the excitation functions adopted in the analysis for the data sets used, these being chosen to represent a real orthogonal cutting process. In an attempt to provide guidance for the selection of technological parameters for the avoidance of primary chatter, this paper determines for the first time the stability regions in terms of the depth of cut and the cutting speed co-ordinates.     

Multiple internal resonances and non-planar dynamics of shallow suspended cables to the harmonic excitations

In the present study, the nonlinear response of a shallow suspended cable with multiple internal resonances to the primary resonance excitation is investigated. The method of multiple scales is applied directly to the nonlinear equations of motion and associated boundary conditions to obtain the modulation equations and approximate solutions of the cable. Frequency–response curves and force–response curves are used to study the equilibrium solution and its stability. The effects of the excitation amplitude on the frequency–response curves of the cable are also analyzed. Moreover, the chaotic dynamics of the shallow suspended cable is investigated by means of numerical simulations.     

A new approach with piecewise-constant arguments to approximate and numerical solutions of oscillatory problems

This paper is devoted to the development of a novel approximate and numerical method for the solutions of linear and non-linear oscillatory systems, which are common in engineering dynamics. The original physical information included in the governing equations of motion is mostly transferred into the approximate and numerical solutions. Therefore, the approximate and numerical solutions generated by the present method reflect more accurately the characteristics of the motion of the systems. Furthermore, the solutions derived are continuous everywhere with good accuracy and convergence in comparing with Runge-Kutta method. An approximate solution is developed for a linear oscillatory problem and compared with its corresponding exact solution. A non-linear oscillatory problem is also solved numerically and compared with the solutions of Runge-Kutta method. Both the graphical and numerical comparisons are provided in the paper. The accuracy of the approximate and numerical solutions can be controlled as desired by the number of terms in the Taylor series and the value of a single parameter used in the present work. Formulae for numerical computation in solving various linear and non-linear oscillatory problems by the new approach are provided in the paper.     

THE DYNAMIC ANALYSIS OF A BEAM-MASS SYSTEM DUE TO THE OCCURRENCE OF TWO-COMPONENT PARAMETRIC RESONANCE

The objective of this paper is an analytical and numerical study of the dynamics of a beam--mass system. Special attention is given to the phenomena arising due to the motion of the attached mass and modal interactions produced by the existence of multi-component, specifically two-component, parametric resonance under primary resonance. The two-component parametric resonance occurs when the sums or the differences among internal frequencies are the same, or close, as the dimensionless speed parameter of the moving mass. The effects of two-component parametric resonance post on dynamic condition are investigated. Resonance generated by more than two-component modes are neglected due to its remote probability of occurrence in nature.The mechanics of the problem is Newtonian. Based on the assumption that when the moving mass is set in motion the mass is assumed to be rolling on the beam, the mechanics, including the effects due to friction and convective accelerations, of the interface between the moving mass and the beam are determined.Based on the Bernoulli-Euler beam theory, the coupled non-linear equations of motion of an inextensible beam with an attached moving mass are derived. By employing Galerkin procedure, the partial differential equations which describe the motion of a beam-mass system are reduced to an initial-value problem with finite dimensions. The method of multiple time scales is applied to consider the solutions and the occurrence of internal resonance of the resulting multi-degree-of-freedom beam--mass system with time dependent coefficients.     

Space-time numerical simulation and validation of analytical predictions for nonlinear forced dynamics of suspended cables

This paper presents space-time numerical simulation and validation of analytical predictions for the finite-amplitude forced dynamics of suspended cables. The main goal is to complement analytical and numerical solutions, accomplishing overall quantitative/qualitative comparisons of nonlinear response characteristics. By relying on an approximate, kinematically non-condensed, planar modeling, a simply supported horizontal cable subject to a primary external resonance and a 1:1, or 1:1 vs. 2:1, internal resonance is analyzed. To obtain analytical solution, a second-order multiple scales approach is applied to a complete eigenfunction-based series of nonlinear ordinary-differential equations of cable damped forced motion. Accounting for both quadratic/cubic geometric nonlinearities and multiple modal contributions, local scenarios of cable uncoupled/coupled responses and associated stability are predicted, based on chosen reduced-order models. As a cross-checking tool, numerical simulation of the associated nonlinear partial-differential equations describing the dynamics of the actual infinite-dimensional system is carried out using a finite difference technique employing a hybrid explicit-implicit integration scheme. Based on system control parameters and initial conditions, cable amplitude, displacement and tension responses are numerically assessed, thoroughly validating the analytically predicted solutions as regards the actual existence, the meaningful role and the predominating internal resonance of coexisting/competing dynamics. Some methodological aspects are noticed, along with a discussion on the kinematically approximate versus exact, as well as planar versus non-planar, cable modeling.     

A general solution procedure for vibrations of systems with cubic nonlinearities and nonlinear/time-dependent internal boundary conditions

The subject of this paper is the development of a general solution procedure for the vibrations (primary resonance and nonlinear natural frequency) of systems with cubic nonlinearities, subjected to nonlinear and time-dependent internal boundary conditions—this is a commonly occurring situation in the vibration analysis of continuous systems with intermediate elements. The equations of motion form a set of nonlinear partial differential equations with nonlinear, time-dependent, and coupled internal boundary conditions. The method of multiple timescales, an approximate analytical method, is applied directly to each partial differential equation of motion as well as coupled boundary conditions (i.e. on each sub-domain and the corresponding internal boundary conditions for a continuous system with intermediate elements) which ultimately leads to approximate analytical expressions for the frequency-response relation and nonlinear natural frequencies of the system. These closed-form solutions provide direct insight into the relationship between the system parameters and vibration characteristics of the system. Moreover, the suggested solution procedure is applied to a sample problem which is discussed in detail.     

The response of two-degree-of-freedom systems with quadratic non-linearities to a combination parametric resonance

The response of two-degree-of-freedom systems with quadratic non-linearities to a combination parametric resonance in the presence of two-to-one internal resonances is investigated. The method of multiple scales is used to construct a first order uniform expansion yielding four first order non-linear ordinary differential equations governing the modulation of the amplitudes and the phases of the two modes. Steady state responses and their stability are computed for selected values of the system parameters. The effects of detuning the internal resonance, detuning the parametric resonance, the phase and magnitude of the second mode parametric excitation, and the initial conditions are investigated. The first order perturbation solution predicts qualitatively the trivial and non-trivial stable steady state solutions and illustrates both the quenching and saturation phenomena. In addition to the steady state solutions, other periodic solutions are predicted by the perturbation amplitude and phase modulation equations. These equations predict a transition from constant steady state non-trivial responses to limit cycle responses (Hopf bifurcation). Some limit cycles are also shown to experience period doubling bifurcations. The perturbation solutions are verified by numerically integrating the governing differential equations.     

Vibration analysis of harmonically excited non-linear system using the method of multiple scales

An analytical method is presented for evaluation of the steady state periodic behavior of non-linear systems. This method is based on the substructure synthesis formulation and a multiple scales procedure, which is applied to the analysis of non-linear responses. A complex non-linear system is divided into substructures, of which equations are approximately transformed to modal co-ordinates including non-linear term under the reasonable procedure. Then, the equations are synthesized into the overall system and the solution of the non-linear system can be obtained. Based on the method of multiple scales, the proposed procedure reduces the size of large-degree-of-freedom problem in solving the non-linear equations. Feasibility and advantages of the proposed method are illustrated by the application of the analytic procedure to the non-linear rotating machine system as a large mechanical structure system. Results obtained are reported to be an efficient approach with respect to non-linear response prediction when compared with other conventional methods.     

Parametric excitation of two internally resonant oscillators

The response of two-degree-of-freedom systems with quadratic non-linearities to a principal parametric resonance in the presence of two-to-one internal resonances is investigated. The method of multiple scales is used to construct a first-order uniform expansion yielding four first-order non-linear ordinary differential (averaged) equations governing the modulation of the amplitudes and the phases of the two modes. These equations are used to determine steady state responses and their stability. When the higher mode is excited by a principal parametric resonance, the non-trivial steady state value of its amplitude is a constant that is independent of the excitation amplitude, whereas the amplitude of the lower mode, which is indirectly excited through the internal resonance, increases with the amplitude of the excitation. However, in addition to Poincaré-type bifurcations, this response exhibits a Hopf bifurcation leading to amplitude- and phase-modulated motions. When the lower mode is excited by a principal parametric resonance, the averaged equations exhibit both Poincaré and Hopf bifurcations. In some intervals of the parameters, the periodic solutions of the averaged equations, in the latter case, experience period-doubling bifurcations, leading to chaos.     

Differential transform method for solving linear and non-linear systems of partial differential equations

In this Letter, we propose a reliable algorithm to develop exact and approximate solutions for the linear and non-linear systems of partial differential equations. The approach rest mainly on two-dimensional differential transform method which is one of the approximate methods. The method can easily be applied to many linear and non-linear problems and is capable of reducing the size of computational work. Exact solutions can also be achieved by the known forms of the series solutions. Several illustrative examples are given to demonstrate the effectiveness of the present method.     

Non-linear resonances in the forced responses of plates,part II: Asymmetric responses of circular plates

The dynamic analogue of the von Karman equations is used to study the forced response, including asymmetric vibrations and traveling waves, of a clamped circular plate subjected to harmonic excitations when the frequency of excitation is near one of the natural frequencies. The method of multiple scales, a perturbation technique, is used to solve the non-linear governing equations. The approach presented provides a great deal of insight into the nature of the non-linear forced resonant response. It is shown that in the absence of internal resonance (i.e., a combination of commensurable natural frequencies) or when the frequency of excitation is near one of the lower frequencies involved in the internal resonance, the steady state response can only have the form of a standing wave. However, when the frequency of excitation is near the highest frequency involved in the internal resonance it is possible for a traveling wave component of the highest mode to appear in the steady state response.     

On non-linear oscillations with slowly varying system parameters

The paper deals with an approximate analysis of non-linear oscillation problems with slowly varying system parameters. From the differential equations for amplitude and phase, set up by the method of variation of parameters, the approximate solutions are obtained by using the generalized averaging method of Sinha and Srinivasan based on ultraspherical polynomial expansions. The Bogoliubov-Mitropolsky results are given by a particular set of these polynomials. Problems of a single degree of freedom system as well as monofrequency oscillations in systems with multiple degrees of freedom are considered. The approach has been illustrated by an example and the results are compared with the numerical solutions. A close agreement is found.     

The Method of Fundamental Solutions for Solving Convection-Diffusion Equations with Variable Coefficients

A meshless method based on the method of fundamental solutions (MFS) is proposed to solve the time-dependent partial differential equations with variable coefficients. The proposed method combines the time discretization and the one-stage MFS for spatial discretization. In contrast to the traditional two-stage process, the one-stage MFS approach is capable of solving a broad spectrum of partial differential equations. The numerical implementation is simple since both closed-form approximate particular solution and fundamental solution are easier to find than the traditional approach. The numerical results show that the one-stage approach is robust and stable.     

On the influence of Poisson's ratio on circular plate natural frequencies

Nayfeh and Kluwick [1] studied simple-wave solutions of non-linear hyperbolic equations comparatively by using three methods: namely, (1) the method of renormalization, (2) the analytic method of characteristics and (3) the method of multiple scales. The method of multiple scales was found to be superior to the other two, and yielded correct results. However, the method seems to involve difficulties in respect to conventional concepts of partial differentiation. In this context, this note presents another method—the Poincaré-Lighthill-Kuo method—which provides not only correct end results but also an elegant procedure.     

A combined parabolic-integral equation approach to the acoustic simulation of vibro-acoustic imaging

This paper aims to model ultrasound vibro-acoustography to improve our understanding of the underlying physics of the technique thus facilitating the collection of better images. Ultrasound vibro-acoustography is a novel imaging technique combining the resolution of high-frequency imaging with the clean (speckle-free) images obtained with lower frequency techniques. The challenge in modeling such an experiment is in the variety of scales important to the final image. In contrast to other approaches for modeling such problems, we break the experiment into three parts: high-frequency propagation, non-linear interaction and the propagation of the low-frequency acoustic emission. We then apply different modeling strategies to each part. For the high-frequency propagation we choose a parabolic approximation as the field has a strong preferred direction and small propagation angles. The non-linear interaction is calculated directly with Fourier methods for computing derivatives. Because of the low-frequency omnidirectional nature of the acoustic emission field and the piecewise constant medium we model the low-frequency field with a surface integral approach. We use our model to compare with experimental data and to visualize the relevant fields at points in the experiment where laboratory data is difficult to collect, in particular the source of the low-frequency field. To simulate experimental conditions we perform the simulations with the two frequencies 3 and 3.05 MHz with an inclusion of varying velocity submerged in water.     

Collocation Methods for a Class of Volterra Integral Functional Equations with Multiple Proportional Delays

In this paper, we apply the collocation methods to a class of Volterra integral functional equations with multiple proportional delays (VIFEMPDs). We shall present the existence, uniqueness and regularity properties of analytic solutions for this type of equations, and then analyze the convergence orders of the collocation solutions and give corresponding error estimates. The numerical results verify our theoretical analysis.     

Analytical approximations to the dynamics of an array of coupled DC SQUIDs

Coupled dynamical systems that operate near the onset of a bifurcation can lead, under certain conditions, to strong signal amplification effects. Over the past years we have studied this generic feature on a wide range of systems, including: magnetic and electric fields sensors, gyroscopic devices, and arrays of loops of superconducting quantum interference devices, also known as SQUIDs. In this work, we consider an array of SQUID loops connected in series as a case study to derive asymptotic analytical approximations to the exact solutions through perturbation analysis. Two approaches are considered. First, a straightforward expansion in which the non-linear parameter related to the inductance of the DC SQUID is treated as the small perturbation parameter. Second, a more accurate procedure that considers the SQUID phase dynamics as non-uniform motion on a circle. This second procedure is readily extended to the series array and it could serve as a mathematical framework to find approximate solutions to related complex systems with high-dimensionality. To the best of our knowledge, an approximate analytical solutions to an array of SQUIDs has not been reported yet in the literature.     

Solutions of the factorization equations from Toda field theory

We apply the quantum inverse scattering method to a discretization of the generalized Toda field theory. In this way we get a family of solutions of the factorization equations. These solutions provide (i) a solvable many particle system with different types of particles interacting via non-relativistic complex potentials; (ii) a presumably solvable generalization of the six-vertex model in a planar lattice.     

The effect of the number of nodal diameters on non-linear interactions in two asymmetric vibration modes of a circular plate

In order to investigate the effect of the number of nodal diameters on non-linear interactions in asymmetric vibrations of a circular plate, a primary resonance of the plate is considered. The plate is assumed to have an internal resonance in which the ratio of the natural frequencies of two asymmetric modes is three to one. The response of the plate is expressed as an expansion in terms of the linear, free oscillation modes, and its amplitude is considered to be small but finite, and the method of multiple scales is used. In view of the corrected solvability conditions for the responses, it has been found that in order for the modes to interact, the ratio of the numbers of nodal diameters of two modes must be either three to one or one to one. In this study the one-to-one case, in which the modes have the same number of nodal diameters, is examined. The non-linear governing equations are reduced to a system of autonomous ordinary differential equations for amplitude and phase variables by means of the corrected solvability conditions. The steady state responses and their stability are determined by using this system. The result shows very complicated interactions between two modes by telling existence of non-vanishing amplitudes of the mode not directly excited.