Generation of the second, fourth, and eighth harmonics by a hyperelastic longitudinal planewave: numerical simulation

The generation of the first, second, fourth, and eighth harmonics of a harmonic longitudinal plane wave propagating in quadratic nonlinear materials described by the classical Murnaghan model is numerically modeled. This situation was also analyzed analytically. The materials are fibrous composites with micro- and nanoscale structure. Using such materials introduces additional real restrictions on the main parameters and allows simulating various real situations: from the generation of all harmonics for feasible distances and times to the generation of the second harmonics for feasible distances and times. The contribution of each approximation to the overall wave pattern is analyzed. It is shown that these approximations affect the prediction of the evolution of the initial wave profile: first there is a tendency to the generation of the second harmonic which then transforms into a tendency to the generation of the fourth and eight harmonics     

Analysis of a quadratic nonlinear hyperelastic longitudinal plane wave

The perturbation (small-parameter) method is used to analyze the propagation of a harmonic longitudinal plane wave in a quadratic nonlinear hyperelastic material described by the classical Murnaghan model. The three first approximations are obtained, and the contribution of each of them into the wave pattern is analyzed. It is shown that the third approximation somewhat improves the prediction of the evolution of the initial waveprofile: the tendency to generate the second harmonic goes over into the tendency to generate the fourth harmonic Translated from Prikladnaya Mekhanika, Vol. 45, No. 2, pp. 46–58, February 2009.     

Two mode sub-harmonic and harmonic vibrations of a non-linear beam forced by a two mode harmonic load

The non-linear transverse vibrations of a uniform beam with ends restrained to remain a fixed distance apart and forced by a two mode function which is harmonic in time, are studied by a corresponding two mode approach. The existence of sub-harmonic response of order 1/3 and harmonic response in the sub-harmonic resonance region of the forcing frequency is proved. Approximate solutions are found by Urabe's numerical method applied to Galerkin's procedure and by an analytical harmonic balance-perturbation method. Error bounds of the Galerkin approximations are given.     

Non-linear wave propagation solutions by Fourier transform perturbation

A Fourier transform perturbation method is developed and used to obtain uniformly valid asymptotic approximations of the solution of a class of one-dimensional second order wave equations with small non-linearities. Multiple time scales are used and the initial-value problem on the infinite line is solved by Fourier transforming the wave equation and expanding the Fourier transform in powers of the small parameter. The non-linearity involves only the first partial derivatives of the dependent variable and the determination of the leading approximation is reduced to the solution of a pair of coupled non-linear ordinary differential equations in Fourier space. Examples are given involving a convolution non-linearity and a Van-der-Pol non-linearity.     

Residue harmonic balance approach to limit cycles of non-linear jerk equations

A residue harmonic balance is established for accurately determining limit cycles to parity- and time-reversal invariant general non-linear jerk equations with cubic non-linearities. The new technique incorporates the salient features of both methods of harmonic balance and parameter bookkeeping to minimize the total residue. The residue is separated into two parts in each step; one conforms to the present order of approximation and the remaining part for use in the next order. The corrections are governed by a set of linear ordinary differential equations that can be solved easily. Three specific cases of non-linear jerk equations are given to illustrate the validity and efficiency. The approximations to the angular frequency and the limit cycle are obtained and compared. The results show that the approximations obtained are in excellent agreement with the exact solutions for a wide range of initial velocities. The new technique is simple in principle and can be applied to other non-linear oscillating systems.     

Interesting effects in harmonic generation by plane elastic waves

The harmonics of plane longitudinal and trans-verse waves in nonlinear elastic solids with up to cubic nonlinearity in a one-dimensional setting are investigated in this paper. It is shown that due to quadratic nonlinearity, a transverse wave generates a second longitudinal harmonic. This propagates with the velocity of transverse waves, as well as resonant transverse first and third harmonics due to the cubic and quadratic nonlinearities. A longitudinal wave generates a resonant longitudinal second harmonic, as well as first and third harmonics with amplitudes that increase linearly and quadratically with distance propagated. In a second investigation, incidence from the linear side of a pri-mary wave on an interface between a linear and a nonlinear elastic solid is considered. The incident wave crosses the interface and generates a harmonic with interface conditions that are equilibrated by compensatory waves propagating in two directions away from the interface. The back-propagated compensatory wave provides information on the nonlinear elastic constants of the material behind the interface. It is shown that the amplitudes of the compensatory waves can be increased by mixing two incident longitudinal waves of appropriate frequencies.     

An application of Discontinuous Galerkin space and velocity discretisations to the solution of a model kinetic equation

An approach based on a Discontinuous Galerkin discretisation is proposed for the Bhatnagar–Gross–Krook model kinetic equation. This approach allows for a high-order polynomial approximation of molecular velocity distribution function both in spatial and velocity variables. It is applied to model one-dimensional normal shock wave and heat transfer problems. Convergence of solutions with respect to the number of spatial cells and velocity bins is studied, with the degree of polynomial approximation ranging from zero to four in the physical space variable and from zero to eight in the velocity variable. This approach is found to conserve mass, momentum and energy when high-degree polynomial approximations are used in the velocity space. For the shock wave problem, the solution is shown to exhibit accelerated convergence with respect to the velocity variable. Convergence with respect to the spatial variable is in agreement with the order of the polynomial approximation used. For the heat transfer problem, it was observed that convergence of solutions obtained by high-degree polynomial approximations is only second order with respect to the resolution in the spatial variable. This is attributed to the temperature jump at the wall in the solutions. The shock wave and heat transfer solutions are in excellent agreement with the solutions obtained by a conservative finite volume scheme.     

Quadratically Nonlinear Cylindrical Hyperelastic Waves: Primary Analysis of Evolution

The propagation and interaction of hyperelastic cylindrical waves are studied. Nonlinearity is introduced by means of the Murnaghan potential and corresponds to the quadratic nonlinearity of all basic relationships. To analyze wave propagation, an asymptotic representation of the Hankel function of the first order and first kind is used. The second-order analytical solution of the nonlinear wave equation is similar to that for a plane longitudinal wave and is the sum of the first and second harmonics, with the difference that the amplitudes of cylindrical harmonics decrease with the distance traveled by the wave. A primary computer analysis of the distortion of the initial wave profile is carried out for six classical hyperelastic materials. The transformation of the first harmonic of a cylindrical wave into the second one is demonstrated numerically. Three ways of allowing for nonlinearities are compared __________ Translated from Prikladnaya Mekhanika, Vol. 41, No. 7, pp. 73–82, July 2005.     

Finite-amplitude waves in a homogeneous fluid with a floating elastic plate

Equations for three nonlinear approximations of a wave perturbation in a homogeneous ideal incompressible fluid covered by a thin elastic plate are obtained using the method of multiple scales and taking into account that the acceleration of vertical flexural displacements of the plate is nonlinear. Based on the obtained equations, asymptotic expansions up third-order terms are constructed for the fluid velocity potential and the perturbations of the plate-fluid interface (plate bending) caused by a traveling periodic wave of finite amplitude. The wave characteristics are analyzed as functions of the elastic modulus and thickness of the plate and the length and tilt of the initial fundamental harmonic wave.     

Dispersion of Small Amplitude Waves in a Pre-Stressed, Compressible Elastic Plate

The dispersion of harmonic waves, propagating along a principal direction in a pre-stressed, compressible elastic plate, is investigated in respect of the most general isotropic strain-energy function. Different cases, dependent on the choice of material parameters and pre-stress, are analysed. A complete long and short wave asymptotic analysis is carried out, with the approximations obtained giving phase speed (and frequency) as explicit functions of wave and mode number. Various wave fronts, both associated with the short wave limit of harmonics and arising through the combination of harmonics in a narrow wave speed region, are discussed. It is mentioned that the case of high compressibility is of particular interest. In contrast with the classical (un-strained) case, the longitudinal body wave speed may be less than the corresponding shear wave speed. In consequence, the short wave limit of all harmonics may be the appropriate longitudinal wave speed; contrasting with the classical case for which this limit is necessarily associated with a shear wave front. A further possible short wave limit is also shown to exist for which the associated wave normal has a component in the direction normal to the plate. Particularly novel numerical results are presented when the longitudinal and shear wave speeds are equal. The analysis is illustrated by numerical calculations for various strain-energy functions.     

Forced response of quadratic nonlinear oscillator: comparison of various approaches

The primary resonances of a quadratic nonlinear system under weak and strong external excitations are investigated with the emphasis on the comparison of different analytical approximate approaches. The forced vibration of snap-through mechanism is treated as a quadratic nonlinear oscillator. The Lindstedt-Poincaré method, the multiple-scale method, the averaging method, and the harmonic balance method are used to determine the amplitude-frequency response relationships of the steady-state responses. It is demonstrated that the zeroth-order harmonic components should be accounted in the application of the harmonic balance method. The analytical approximations are compared with the numerical integrations in terms of the frequency response curves and the phase portraits. Supported by the numerical results, the harmonic balance method predicts that the quadratic nonlinearity bends the frequency response curves to the left. If the excitation amplitude is a second-order small quantity of the bookkeeping parameter, the steady-state responses predicted by the second-order approximation of the LindstedtPoincaré method and the multiple-scale method agree qualitatively with the numerical results. It is demonstrated that the quadratic nonlinear system implies softening type nonlinearity for any quadratic nonlinear coefficients.     

Modeling nonlinear Rayleigh wave fields generated by angle beam wedge transducers—A theoretical study

Nonlinear Rayleigh wave fields generated by an angle beam wedge transducer are modeled in this study. The calculated area sound sources underneath the wedge are used to model the fundamental Rayleigh sound fields on the specimen surface, which are more accurate than the previously used line sources with uniform or Gaussian amplitude distributions. A general two-dimensional nonlinear Rayleigh wave equation without parabolic approximation is introduced and the solutions are obtained using the quasilinear theory. The second harmonic Rayleigh wave due to material nonlinearity is given in an integral expression with these fundamental Rayleigh waves radiated by the wedge transmitter acting as a forcing function. Multi-Gaussian beam (MGB) models are employed to simplify these integral solutions and to extract the diffraction and attenuation correction terms explicitly. The effect of nonlinearity of generating sources on the second harmonic Rayleigh wave fields is taken into consideration; simulation results show that it will affect the magnitude and diffraction correction of the second harmonic waves in the region close to the Rayleigh wave sound sources. This research provides a theoretical improvement to alleviate the experimental restriction on analyzing the effects of diffraction, attenuation and source nonlinearity when using angle beam wedge transducers as transmitters.     

The Simplest Resonance Capture Problem,Using Harmonic Balance Based Averaging

We study the dynamics of capture into, or escape from, resonance in a strongly nonlinear oscillator with weak damping and forcing, using harmonic balance based averaging (HBBA). This system provides the simplest example of resonance capture that we know of. The HBBA technique, here adapted to tackle nonlinear resonances, provides a harmonic balance assisted approximation to the underlying, asymptotically correct, averaged dynamics. Allowing the harmonic balance approximation makes a variety of systems analytically tractable which might otherwise be intractable. The evolution equations for amplitude and phase of oscillations are derived first. Restricting attention near the primary resonance, the slow flow equations are approximately averaged. The resulting flow transparently shows the stable and unstable primary resonant solutions, as well as the trajectories that get captured into resonance and the ones that escape. Good agreement with numerics is obtained, showing the utility of HBBA near resonance manifolds.     

On Approximations of First Integrals for a System of Weakly Nonlinear,Coupled Harmonic Oscillators

In this paper a system of weakly nonlinear, coupled harmonic oscillatorswill be studied. It will be shown that the recently developedperturbation method based on integrating vectors can be used toapproximate first integrals and periodic solutions. To show how thisperturbation method works the method is applied to a system of weaklynonlinear, coupled harmonic oscillators with 1:3 and 3:1 internalresonances. Not only approximations of first integrals will be given,but it will also be shown how, in a rather efficient way, the existenceand stability of time-periodic solutions can be obtained from theseapproximations. In addition some bifurcation diagrams for a set ofvalues of the parameters will be presented.     

结构参数大修改时的特征值重分析方法

就结构参数发生大修改的情况提出了两种高精度的特征值重分析方法：Ｐａｄｅ　逼近法和推广的Ｋｉｒｓｃｈ混合法．利用这两种方法,计算了一个具有２０２个结点,３５７个梁单元的平面框架的近似特征值．计算结果表明,所提出的方法是结构参数修改时的特征值重分析的有效方法．     

Electromagnetic diffraction by a thin conducting circular disk placed at the plane interface of two different media

In this paper, the problem of diffraction of time harmonic, electromagnetic waves by a thin ideally conducting disk lying at the plane interface of two different media is considered. In this analysis, the incident wave is a plane wave travelling in a direction perpendicular to the plane interface of the two media. A Hertz vector formulation is applied to reduce our electromagnetic diffraction problem to a system of two scalar problems which are solved by the help of two pairs of Fredholm integral equations of the second kind. Low frequency approximations to the tangential components of the magnetic intensity associated with the diffracted field at the surfaces of the disk, the induced surface current density on the disk and the scattering cross section are obtained.     

Resonant waves in elastic structured media: Dynamic homogenisation versus Green’s functions

We address an important issue of dynamic homogenisation in vector elasticity for a doubly periodic mass-spring elastic lattice. The notion of logarithmically growing resonant waves is used in the analysis of star-shaped wave forms induced by an oscillating point force. We note that the dispersion surfaces for Floquet–Bloch waves in the elastic lattice may contain critical points of the saddle type. Based on the local quadratic approximations of a dispersion surface, where the radian frequency is considered as a function of wave vector components, we deduce properties of a transient asymptotic solution associated with the contribution of the point source to the wave form. The notion of local Green’s functions is used to describe localised wave forms corresponding to the resonant frequency. The special feature of the problem is that, at the same resonant frequency, the Taylor quadratic approximations for different groups of the critical points on the dispersion surfaces (and hence different Floquet–Bloch vectors) are different. Thus, it is shown that for the vector case of micro-structured elastic medium there is no uniformly defined dynamic homogenisation procedure for a given resonant frequency. Instead, the continuous approximation of the wave field can be obtained through the asymptotic analysis of the lattice Green’s functions, presented in this paper.     

Propagation of shock waves in a medium with exponentially varying density

The results of Raizer [1], Hays [2], and Chernous'ko [3] are generalized to-the case of self-similar propagation of shock waves in a gas with exponentially varying density and constant pressure. A solution is found by the method of successive approximations. The zero-order approximation coincides with the Whitham method [4]. The first-order approximation is in good agreement with numerical calculations in [2]. The non-selfsimilar motion of a weak shock wave is investigated in the framework of linear theory.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 48–54, November–December, 1970.