Primary wave transmission in systems of elastic rods with granular interfaces

We study analytically and numerically primary pulse transmission in one dimensional systems of identical linearly elastic non-dispersive rods separated by identical homogeneous granular layers composed of n beads. The beads interact elastically through a strongly (essentially) nonlinear Hertzian contact law. The main challenge in studying pulse transmission in such strongly nonlinear media is to analyze the ‘basic problem’, namely, the dynamical response of a single intermediate granular layer, confined from both ends by barely touching linear elastic rods subject to impulsive excitation of the left rod. The analysis of the basic problem is carried out under two basic assumptions; namely, of sufficiently small duration of the shock excitation applied to the first layer of the system, and of sufficiently small mass of each bead in the granular interface compared to the mass of each rod. In fact, the smallness of the mass of the bead defines the small parameter in the asymptotic analysis of this problem. Both assumptions are reasonable from the point of view of practical applications. In the analysis we focus only in primary pulse propagation, by neglecting secondary pulse reflections caused by wave scattering at each granular interface and considering only the transmission of the main (primary) pulse across the interface to the neighboring elastic rod. Two types of shock excitations are considered. The first corresponds to fixed time duration (but still much smaller compared to the characteristic time of pulse propagation through the length of each rod), whereas the second type corresponds to a pulse duration that depends on the small parameter of the problem. The influence of the number of beads of the granular interface on the primary wave transmission is studied, and it is shown that at granular interfaces with a relatively low number of beads fast time scale oscillations are excited with increasing amplitudes with increasing number of beads. For a larger number of beads, primary pulse transmission is by means of solitary wave trains resulting from the dispersion of the original shock pulse; in that case fast oscillations result due to interference phenomena caused by the scattering of the main pulse at the boundary of the interface. Considering a periodic system of rods we demonstrate significant reduction of the primary pulse when transmitted through a sequence of granular interfaces. This result highlights the efficacy of applying granular interfaces for passive shock mitigation in layered elastic media.     

Stability and bifurcation analysis of a nonlinear transversally loaded rotating shaft

The dynamic stability and self-excited posteritical whirling of rotating transversally loaded shaft made of a standard material with elastic and viscous nonlinearities are analyzed in this paper using the theory of bifurcations as a mathematical tool. Partial differential equations of motion are derived under assumption that von Karman's nonlinearity is absent but geometric curvature nonlinearity is included. Galerkin's first-mode discretization procedure is then applied and the equations of motion are transformed to two third-order nonlinear equations that are analyzed using the theory of bifurcation. Condition for nontrivial equilibrium stability is determined and a bifurcating periodic solution of the second-order approximation is derived. The effects of dimensionless stress relaxation time and cubic elastic and viscous nonlinearities as well as the role of the transverse load are studied in the exemplary numerical calculations. A strongly stabilizing influence of the relaxation time is found that may eliminate self-excited vibration at all. Transition from super- to subcritical bifurcation is observed as a result of interaction between system nonlinearities and the transverse load.     

Treatment of nonlinear terms in the numerical solution of partial differential equations for multiphase flow in porous media

Simultaneous flow of two or three immiscible fluids in porous media is modelled by a system of coupled, nonlinear partial differential equations. These equations are reduced to a system of nonlinear algebraic equations through the use of finite-difference approximations for derivatives. Several types of nonlinearities requiring careful analysis exist in this model. Here, we present a systematic study of all available, and some new, methods for the treatment of nonlinearitics in this model. It is believed that the solution techniques presented here may also prove useful for other strongly nonlinear partial differential equations.     

Nonlinear Parameter Identification of a Mechanical Interface Based on Primary Wave Scattering

We study stress-wave propagation in an impulsively forced split Hopkinson bar system incorporating a threaded interface. We first consider only primary transmission and reflection and reduce the problem to a first-order, strongly nonlinear ordinary differential equation governing the displacement across the interface, called the primary-pulse model. The interface is modeled as an adjusted-Iwan element, which is characterized by matching experimental and numerical eigenfrequencies as well as primary pulse amplitudes. We find that the adjusted-Iwan element parameters are dependent on preload and impact velocity (input force). A high-order finite element model paired with the identified adjusted-Iwan element is used to simulate multiple transmissions and reflections across the interface. We find that the finite element simulation reproduces the experimental results in both the wavelet and Fourier domains, validating the identification method. Our findings demonstrate that the primary-pulse model can be used for experimental parameter identification of nonlinear interfaces in waveguides.     

Relationship between multibody dynamics computation and nonlinear vibration theory

This paper presents the ground-work of implementing the multibody dynamics codes to analyzing nonlinear coupled oscillators. The recent developments of the multibody dynamics have resulted in several computer codes that can handle large systems of differential and algebraic equations (DAE). However, these codes cannot be used in their current format without appropriate modifications. According to multibody dynamics theory, the differential equations of motion are linear in the acceleration, and the constraints are appended into the equations of motion through Lagrange's multipliers. This formulation should be able to predict the nonlinear phenomena established by the nonlinear vibration theory. This can be achieved only if the constraint algebraic equations are modified to include all the system kinematic nonlinearities. This modification is accomplished by considering secondary nonlinear displacements which are ignored in all current codes. The resulting set of DAE are solved by the Gear stiff integrator. The study also introduced the concept of constrained flexibility and uses an instantaneous energy checking function to improve integration accuracy in the numerical scheme. The general energy balance is a single scalar equation containing all the energy component contributions. The DAE solution is then compared with the solution predicted by the nonlinear vibration theory. It also establishes new foundation for the use of multibody dynamics codes in nonlinear vibration problems. It is found that the simulation CPU time is much longer than the simulation of the original equations of the system.     

Propagation of surface (seismic) waves: Ordinary differential equations with strongly nonlinear damping

Second-order ordinary differential equations (ODEs) with strongly nonlinear damping (cubic nonlinearities) govern surface wave motions that entail nonlinear surface seismic motions. They apply to dynamic crack propagation and nonlinear oscillation problems in physics and nonlinear mechanics. It is shown that the nonlinear surface seismic wave equation (Rayleigh equation) admits several functional transformations and it is possible to reduce it to an equivalent first-order Abel ODE of the second kind in normal form. Based on a recently developed methodology concerning the construction of exact analytic solutions for the type of Abel equations under consideration, exact solutions are obtained for the nonlinear seismic wave (NLSW) equation for initial conditions of the physical problem. The method employed is general and can be applied to a large class of relevant ODEs in mathematical physics and nonlinear mechanics.     

Character of motion of oscillating pendulum systems at the boundary of the stable region

A study is made of a nonlinear dynamic system in which the two of the eigenvalues of the matrix of equations of the linear approximation have negative parts and two other eigenvalues have purely imaginary parts. The nonlinearities are of the third order. The Poincare-Lyapunov-Malkin method is used to show that the motion of such systems is generally aperiodic. Periodic motions can occur only in isolated cases. Institute of Railroad Transport, Kiev; Highway Institute of Donetsk Technical University, Gorlovka, Ukraine. Translated from Prikladnaya Mekhanika, Vol. 35, No. 8, pp. 101–107, August, 1999.     

Vibration and stability of an aircraft tail under simultaneous primary-combined and internal resonance

This paper adds a negative velocity feedback to the dynamical system of twin-tail aircraft to suppress the vibration. The system is represented by two coupled second-order nonlinear differential equations having both quadratic and cubic nonlinearities. The system describes the vibration of an aircraft tail subjected to both multi-harmonic and multi-tuned excitations. The method of multiple time scale perturbation is adopted to solve the nonlinear differential equations and obtain approximate solutions up to the third order approximations. The stability of the proposed analytic solution near the simultaneous primary, combined and internal resonance is studied and its conditions are determined. The effect of different parameters on the steady state response of the vibrating system is studied and discussed by using frequency response equations. Some different resonance cases are investigated numerically     

Primary pulse transmission in coupled steel granular chains embedded in PDMS matrix: Experiment and modeling

We present an experimental study of primary pulse transmission in coupled ordered steel granular chains embedded in poly-di-methyl-siloxane (PDMS) elastic matrix. Two granular one-dimensional chains are considered (an ‘excited’ and an ‘absorbing’ one), each composed of 11 identical steel beads of 9.5 mm diameter with the centerline of the chain spaced at fixed distances of 0.5, 1.5 or 2.5 mm apart. We directly force one of the chains (the excited one) by a transient pulse and measure, by means of laser vibrometry, the primary transmitted pulses at the end beads of both chains and at the first bead of the absorbing chain. It is well known that the dynamics of this type of ordered granular media is strongly nonlinear due, (i) to Hertzian interactions between adjacent beads, and (ii) to possible bead separations in the absence of compressive forces and ensuing collisions between neighboring beads. Accordingly, we develop a strongly nonlinear theoretical model that takes into account the coupling of the granular chains due to the PDMS matrix, with the aim to model primary pulse transmission in this system. After validating the model with experimental measurements, we employ it in a predictive fashion to estimate energy transfer between chains as a function of the interspatial distance between chains. Furthermore, based on this model we perform predictive matrix design to achieve maximum energy transfer from the excited to the absorbing chain, and provide a theoretical explanation of the nonlinear dynamics governing energy transfer (including energy equi-partition) in this system.     

Vibration of fluid-conveying pipe with nonlinear supports at both ends

The axial fluid-induced vibration of pipes is very widespread in engineering applications. The nonlinear forced vibration of a viscoelastic fluid-conveying pipe with nonlinear supports at both ends is investigated. The multi-scale method combined with the modal revision method is formulated for the fluid-conveying pipe system with nonlinear boundary conditions. The governing equations and the nonlinear boundary conditions are rescaled simultaneously as linear inhomogeneous equations and linear inhomogeneous boundary conditions on different time-scales. The modal revision method is used to transform the linear inhomogeneous boundary problem into a linear homogeneous boundary problem. The differential quadrature element method (DQEM) is used to verify the approximate analytical results. The results show good agreement between these two methods. A detailed analysis of the boundary nonlinearity is also presented. The obtained results demonstrate that the boundary nonlinearities have a significant effect on the dynamic characteristics of the fluid-conveying pipe, and can lead to significant differences in the dynamic responses of the pipe system.

     

Nonlinear energy harvesting from vibratory disc-shaped piezoelectric laminates

Implementing resonators with geometrical nonlinearities in vibrational energy harvesting systems leads to considerable enhancement of their operational bandwidths. This advantage of nonlinear devices in comparison to their linear counterparts is much more obvious especially at small-scale where transition to nonlinear regime of vibration occurs at moderately small amplitudes of the base excitation. In this paper the nonlinear behavior of a disc-shaped piezoelectric laminated harvester considering midplane-stretching effect is investigated. Extended Hamilton's principle is exploited to extract electromechanically coupled governing partial differential equations of the system. The equations are firstly order-reduced and then analytically solved implementing perturbation method of multiple scales. A nonlinear finite element method(FEM) simulation of the system is performed additionally for the purpose of verification which shows agreement with the analytical solution to a large extent. The frequency response of the output power at primary resonance of the harvester is calculated to investigate the effect of nonlinearity on the system performance. Effect of various parameters including mechanical quality factor, external load impedance and base excitation amplitude on the behavior of the system are studied. Findings indicate that in the nonlinear regime both output power and operational bandwidth of the harvester will be enhanced by increasing the mechanical quality factor which can be considered as a significant advantage in comparison to linear harvesters in which these two factors vary in opposite ways as quality factor is changed.     

Performances of dynamic vibration absorbers for beams subjected to moving loads

The goal of this work is a general assessment regarding the performances of linear and nonlinear dynamic vibration absorbers (DVAs) applied to the specific problem of moving loads or vehicles. The problem consists of a simply supported linear Euler–Bernoulli beam excited with a moving load/vehicle; a DVA is connected to the beam in order to reduce the vibrations. The moving vehicle is modeled by a single degree of freedom mass spring system. The partial differential equations governing the beam dynamics is reduced to a set of ordinary differential equations by means of the Bubnov–Galerkin method. A parametric analysis is carried out to find the optimal parameters of the DVA that minimize the maximum vibration amplitude of the beam. For the case of a moving vehicle, the energy absorbed by the DVA is evaluated. Comparisons among the performances of different types of linear and DVAs are carried out. The goal is to clarify if the use of nonlinearities in the DVAs can effectively improve their performances. The study shows that the most effective type of DVA for the test cases considered is the piecewise linear elastic restoring force.     

Primary Bjerknes force in a three-dimensional nonlinear resonator: Numerical simulations

The effect of the primary Bjerknes force caused by a three-dimensional nonlinear standing ultrasonic wave on a population of nonlinear oscillating bubbles is studied in this paper by analyzing the results obtained from simulations performed with a numerical model at low and moderate pressure amplitudes. Small air bubbles are evenly distributed in a water filled cavity excited at resonance for which axial symmetry is assumed. Both the bubble oscillation variable and the pressure variable are unknown in the nonlinear set of coupled differential equations that describes the interaction of ultrasound and bubbles. Simulation results show that the three-dimensional primary Bjerknes force field is strongly amplitude dependent. We also analyze whether taking the term in the differential system that defines the nonlinear behavior of the pressure field into account is determinant or not on the computation of the force field. The corresponding results corroborate the one-dimensional conclusions on the fundamental importance of considering this nonlinear acoustic term to obtain an accurate approximation of the force in a cavity.     

Fluid Flow-Induced Nonlinear Vibration of Suspended Cables

The nonlinear interaction of the first two in-plane modes of a suspended cable with a moving fluid along the plane of the cable is studied. The governing equations of motion for two-mode interaction are derived on the basis of a general continuum model. The interaction causes the modal differential equations of the cable to be non-self-adjoint. As the flow speed increases above a certain critical value, the cable experiences oscillatory motion similar to the flutter of aeroelastic structures. A co-ordinate transformation in terms of the transverse and stretching motions of the cable is introduced to reduce the two nonlinearly coupled differential equations into a linear ordinary differential equation governing the stretching motion, and a strongly nonlinear differential equation for the transverse motion. For small values of the gravity-to-stiffness ratio the dynamics of the cable is examined using a two-time-scale approach. Numerical integration of the modal equations shows that the cable experiences stretching oscillations only when the flow speed exceeds a certain level. Above this level both stretching and transverse motions take place. The influences of system parameters such as gravity-to-stiffness ratio and density ratio on the response characteristics are also reported.     

Study of Two-Dimensional Axisymmetric Breathers Using Padé Approximants

We analyze axisymmetric, spatially localized standing wave solutions with periodic time dependence (breathers) of a nonlinear partial differential equation. This equation is derived in the 'continuum approximation' of the equations of motion governing the anti-phase vibrations of a two-dimensional array of weakly coupled nonlinear oscillators. Following an asymptotic analysis, the leading order approximation of the spatial distribution of the breather is shown to be governed by a two-dimensional nonlinear Schrödinger (NLS) equation with cubic nonlinearities. The homoclinic orbit of the NLS equation is analytically approximated by constructing [2N × 2N] Padé approximants, expressing the Padé coefficients in terms of an initial amplitude condition, and imposing a necessary and sufficient condition to ensure decay of the Padé approximations as the independent variable (radius) tends to infinity. In addition, a convergence study is performed to eliminate 'spurious' solutions of the problem. Computation of this homoclinic orbit enables the analytic approximation of the breather solution.     

Chaotic Oscillation of Saddle Form Cable-Suspended Roofs under Vertical Excitation Action

The purpose of this paper is to investigate the dynamic behaviour of saddle form cable-suspended roofs under vertical excitation action. The governing equations of this problem are system of nonlinear partial differential and integral equations. We first establish a spectral equation, and then consider a model with one coefficient, i.e., a perturbed Duffing equation. The analytical solution is derived for the Duffing equation. Successive approximation solutions can be obtained in likely way for each time to only one new unknown function of time. Numerical results are given for our analytical solution. By using the Melnikov method, it is shown that the spectral system has chaotic solutions and subharmonic solutions under determined parametric conditions.     

两跨输电线非线性主共振响应分析

研究两跨输电线非线性共振响应问题,应用Hamilton变分原理推导了两跨输电线的振动微分方程以及对应的边界条件。利用Galerkin离散方法和多尺度法,得到了单模态主共振响应。研究结果表明:幅频响应曲线表现出软、硬弹簧性质,随着外激励幅值的增大,输电线系统响应由软弹簧性质向硬弹簧性质转换;系统阻尼减小或外激励幅值增大时,系统幅值个数也随之发生变化,表现出多值和跳跃现象。     

On the existence of non-polarized azimuthal-symmetric electromagnetic waves in circular dielectric waveguide filled with nonlinear isotropic homogeneous medium

The propagation of monochromatic nonlinear symmetric hybrid waves in a cylindrical nonlinear dielectric waveguide is considered. The physical problem is reduced to solving a transmission eigenvalue problem for a system of ordinary differential equations. Spectral parameters of the problem are propagation constants of the waveguide. The problem is reduced to the new type of nonlinear eigenvalue problem. The analytical method of solving this problem is presented. New propagation regime is found.     

谐波平衡法与复平均法计算强非线性系统稳态响应的区别

近些年,很多学者致力于利用非线性增强振动响应减少的效果或者能量采集器的效率。因而非线性系统的响应值需要从理论计算方面更准确地预测。另外,根据学者已取得的研究成就,非线性能量汇(NES)中存在的立方刚度非线性可以将结构中宽频域的振动能量传递至非线性振子部分。文章将一种由NES和压电能量采集器组成的NES-piezo装置与两自由度主结构耦合连接,系统受谐和激励作用。文章采用谐波平衡法和复平均法分别推导了系统稳态响应,参照数值结果,对比两种近似解析方法在求解强非线性系统稳态响应时的异同。计算结果表明,系统体现较弱非线性时,二者计算结果差异很小;当系统体现强非线性时,复平均法不能准确地呈现系统高阶响应,提高阶数的谐波平衡法能更准确地表示系统响应值。基于谐波平衡法和数值算法,讨论NES-piezo装置对于系统宽频域减振的影响。与仅加入非线性能量汇情况对比,结果表明NES-piezo装置不会恶化宽频域减振效果,并且在第一阶共振频率附近,可以稍微提高结构减振效率。另外,计算结果也表明,采用恰当的NES-piezo装置可实现宽频域范围的结构减振和压电能量采集一体化。此项研究工作为研究不同情形强非线性系统的响应提供了理论方法的指导。另外,研究结果也为宽频域范围的结构减振和压电能量采集一体化提供了理论依据。     

随机外激励作用下非线性系统响应演变概率密度

对随机高斯外激励作用下强非线性振动系统响应演变概率密度函数求解问题进行探讨.应用随机函数空间的正交分解理论,将由熵方法定义的指数形式概率密度函数表达式在随机泛函空间中展开,推导了展开级数所满足的FPK方程.运用加特金方法,将概率密度与系统状态向量共同表征的偏微分方程求解问题转化为求解逼近系数的一阶常微分方程组形式,使得问题求解成为可能.数值算例中研究了随机外激励作用下下一阶与二阶随机非线性系统响应概率密度函数求解问题,初步讨论了随机非线性系统响应概率密度函数的瞬态演化过程.