Sensitivity of a supersonic boundary layer to acoustic disturbances

The results obtained by the authors in [1] are extended to the case of arbitrary angles of incidence of the external wave. This is not a trivial generalization, since the acoustic scattering undergoes a qualitative change. It is possible to distinguish two excitation channels: the first is connected with the diffraction of the acoustic wave by the spatial inhomogeneity resulting from the displacing action of the boundary layer, and the second with the presence of concentrated acoustic field sources associated with the scattering of the wave at the leading edge. The latter makes the principal contribution to the initial amplitude of the unstable modes when the angles of incidence of the sound are substantially different from zero. At low angles of incidence there is a singularity which can be revealed by introducing narrow intervals in the neighborhood of the limiting values of the wave numbers, where the two excitation channels are approximately equivalent. It is possible to obtain composite expressions for the initial amplitudes of the unstable modes uniformly valid for all angles of incidence of the acoustic wave.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, pp. 40–47, January–February, 1992.     

Cutting-off effect at reflection–transmission of acoustic waves in anisotropic media with sliding-contact interfaces

The specific feature of the interface, which maintains sliding contact between elastic media, is that it can be impervious to the wave field existing in one of the adjoined materials. As a result, reflection–transmission of plane acoustic waves at the sliding-contact interface may enjoy the cutting-off effect, which implies that neither bulk, nor inhomogeneous modes are being transmitted at particular angles of incidence. The necessary and sufficient criteria for this phenomenon are obtained for a binary structure, constituted by two elastic half-spaces in sliding contact, and for a sandwich structure with sliding-contact interfaces between the enclosed layer and the substrates. In the generic case of unrestricted anisotropy (triclinic materials), the criterion for cutting-off in a binary structure involves acoustic parameters of solely that of the half-spaces, which contains the incident mode, and proves to be independent of an adjacent medium. The frequency-dispersive criterion for the absence of transmission through a triclinic layer in the sliding-contact sandwich structure is independent of substrates. By appeal to the Stroh formalism, the cutting-off conditions in a binary and a sandwich structure are further elaborated under the assumption that one of the half-spaces, or a layer, is orthorhombic, and its two symmetry planes are parallel, respectively, to the plane of incidence and to the sliding-contact interface with arbitrary adjacent media. It is shown that the transmission cut-off in a binary structure is necessarily accompanied by the absence of mode conversion at reflection, but the reverse is not true. The angles of incidence which give rise to these effects are determined in terms of elastic coefficients. Transmission cut-off through an orthorhombic layer comes about at an arbitrary angle of incidence, related to guided-modes range in the layer, for the corresponding aperiodic infinite set of the frequency values. Relations for the coefficients of reflection and transmission at the sliding-contact interface between two orthorhombic half-spaces are obtained in concise form, expressed solely via normal components of the partial Stroh-normalized traction amplitudes. Provided that the adjoined orthorhombic half-spaces in sliding contact are identical, the same value of wave-vector tangential projection, which stipulates transmission cut-off at the incidence of, say, the fast mode, entails total transmission at the incidence of the slow mode.     

Lamb wave scattering by a symmetric pair of surface-breaking cracks in a plate

The scattering problem of a Lamb wave incident on a symmetric pair of surface-breaking transverse cracks in a plate is considered. The Lamb wave is assumed to be obliquely incident on the crack plane. Since the cracks are part-through, the scattered field will contain reflected as well as transmitted waves. The energy of the incoming wave is partitioned into reflected and transmitted wave modes. Energy coefficients of the reflected and transmitted waves are calculated as a function of incident frequency and crack depth. The incidence angle of the incoming wave is also treated as a parameter. Both the reflected and transmitted wave fields are considered as linear superpositions of all real and complex wave modes in the plate. Decomposition of modes is achieved with the help of an orthogonality condition based on the principle of reciprocal work. Continuity of displacement and stress fields is imposed at the crack plane. Energy coefficients for reflection and transmission are obtained from the mode amplitudes. Energy coefficients are shown to be a strong function of incident frequency and crack depth. Experiments are conducted with a PZT transducer network interacting with a symmetric pair of machined cracks in an aluminum plate. Trends predicted by the analysis are reflected in the experimental results.     

Phase-resolved characterization of vortex shedding in the near wake of a square-section cylinder at incidence

The vortex formation and shedding process in the near wake region of a 2D square-section cylinder at incidence has been investigated by means of particle image velocimetry (PIV). Proper orthogonal decomposition (POD) is used to characterize the coherent large-scale flow unsteadiness that is associated with the wake vortex shedding process. A particular application of the POD analysis is to extract the vortex-shedding phase of individual velocity fields, which were acquired at asynchronous low rate with respect to the vortex shedding cycle. The phase of an individual flow field is determined from its projection on the first pair of POD modes, allowing phase averaging of the measurement data to be performed. In addition, a low-order representation of the flow, constructed from the mean and the first pair of POD modes, is found to be practically equivalent to the phase-averaged results. It is shown that this low-order representation corresponds to the basic Fourier component of the flow field ensemble with respect to the reconstructed phase. The phase-averaged flow representations reveal the dominant flow features of the vortex-shedding process and the effect of the angle of incidence upon it.     

On the theory of reflection from a wire grid parallel to an interface between homogeneous media (II)

Summary A previous analysis is extended to include coupling between the TM and TE modes. This can be important when the incidence and polarization are both arbitrary.     

Nonlinear dynamics of an elastic rod with frictional impact

A model is presented for the impact with friction of a flexible body in translation and rotation. This model consists of a system of nonlinear differential equations which considers the multiple collisions as well as frictional effects at the contacting end, and allows one to predict the rigid and elastic body motion after the impact. The kinetic energy is derived by utilizing a generalized velocity field theory for elastic solids. The model uses a dry coefficient of friction and a nonlinear contact force. We introduce a finite number of vibrational modes to take into account the vibrational behavior of the body during impact. The vibrations, the multiple collisions, and the angle of incidence angle, are found to be important factors for the kinematics of frictional impact. Analytical and experimental results were compared to establish the accuracy of the model.     

Propagation of surface waves in a half-space with vertical,inclined or curved interfaces

The paper reviews the author's theoretical studies on reflection and transmission of surface waves at vertical and near-vertical boundaries. For estimation of reflection and transmission coefficients a new approximate method based on the Green function technique is proposed. The method is valid if a body waves field arised at the boundary is approximately orthogonal to the fields of reflected and transmitted normal modes. It may be used for normal as well as for oblique incidence. A Green function in each quarterspace is assumed to be the same as in a halfspace with the same vertical distribution of elastic parameters. For large angles of incidence when the reflected wave field is strong enough, the Green function may be corrected for the reflected wave. The results of calculation for some models by means of this method are discussed.     

Propagation behavior of elastic waves in one-dimensional piezoelectric/piezomagnetic phononic crystal with line defects

Using a stiffness matrix method, we in- vestigate the propagation behaviors of elastic waves in one-dimensional （1D） piezoelectric/piezomagnetic （PE/PM） phononic crystals （PCs） with line defects by calculating energy reflection/transmittion coefficients of quasi-pressure and quasi-shear waves. Line defects are created by the re- placement of PE or PM constituent layer. The defect modes existing in the first gap are considered and the influences on defect modes of the material properties and volume fraction of the defect layers, the type of incident waves, the location of defect layer and the number of structural layers are discussed in detail. Numerical results indicate that defect modes are the most obvious when the defect layers are inserted in the middle of the perfect PCs; the types of incidence wave and material properties of the defect layers have important effects on the numbers, the location of frequencies and the peaks of defect modes, and the defect modes are strongly de- pendent on volume fraction of the defect layers. We hope this paper will be found useful for the design of PE/PM acoustic filters or acoustic transducer with PCs structures.     

Spatial Instabilities of the Incompressible Attachment-Line Flow Using Sparse Matrix Jacobi–Davidson Techniques

We consider the linear stability of incompressible attachment-line flow within the spatial framework. No similarity or symmetry assumptions for the instability modes are introduced and the full two-dimensional representation of the modes is used. The perturbation equations are discretized on a two-dimensional staggered grid. A high order finite difference scheme has been developed which gives rise to a large, sparse, quadratic, eigenvalue problem for the instability modes. The benefits of the Jacobi–Davidson method for the solution of this eigenvalue system are demonstrated and the approach is validated in some detail. Spatial stability results are presented subsequently. In particular, instability predictions at very high Reynolds numbers are obtained which show almost equally strong instabilities for symmetric and antisymmetric modes in this regime.     

多自由度内共振系统非线性模态的分岔特性

利用多尺度法构造了一个立方非线性1:3内共振系统的内共振非线性模态(NonlinearNormal Modes associated with internal resonance)．研究表明,内共振非线性系统除存在单模态运动外还存在耦合模态运动．耦合内共振模态具有分岔特性．利用奇异性理论对模态分岔方程进行分析发现此类系统的模态存在叉形点分岔和滞后点分岔这两种典型的分岔模式．     

Stability of a rubber half-space

The paper treats the stability of surface waves generated when a rubber half-space is subjected to compression. At a load level below the critical one an asymptotic expansion is given for the difference between the potential energy of an adjacent state and that of the fundamental state. The displacement field is expressed approximately by a linear combination of two different buckling modes and a residual displacement field which is orthogonal to the former fields. This remainder permits us to take into account the effect of all other modes which have been neglected. The wavelength of the two modes are governed by the dominant imperfections of the half-space. Terms up to the fourth order in the amplitudes of the buckling modes are included.Results are presented to show the most severe post-buckling behaviour as presented by the line of steepest descent in the load displacement diagram. Furthermore curves showing the reduction of the critical load on account of the imperfections in the two dominant modes are also presented. The analysis is kept sufficiently general to include the effect of pre-straining the rubber half-space upon the post-buckling behaviour.     

Bifurcation of nonlinear internal resonant normal modes of a class of multi-degree-of-freedom systems

The method of multiple scales is applied for constructing nonlinear normal modes (NNMs) of a three-degree-of-freedom system which is discretized from a two-link flexible arm connected by a nonlinear torsional spring. The discrete system is with cubic nonlinearity and 1:3 internal resonance between the second and the third modes. The approximate solution for the NNM associated with internal resonance are presented. The NNMs determined here tend to the linear modes as the nonlinearity vanishes, which is significant for one to construct NNM. Greatly different from results of those nonlinear systems without internal resonance, it is found that the NNM involved in internal resonance include coupled and uncoupled two kinds. The bifurcation analysis of the coupled NNM of the system considered is given by means of the singularity theory. The pitchfork and hysteresis bifurcation are simultaneously found. Therefore, the number of NNM arising from the internal resonance may exceed the number of linear modes, in contrast with the case of no internal resonance, where they are equal. Curves displaying variation of the coupling extent of the coupled NNM with the internal-resonance-deturing parameter are proposed for six cases.     

On the orthogonality of natural modes of vibration

The natural modes of a discrete linear system are orthogonal with respect to the mass and stiffness matrices in a generalized sense. However, these modes are usually not orthogonal to each other in the ordinary sense. The purpose of this paper is to document a number of conditions under which the modes are also orthogonal in an ordinary fashion.     

A new appraoch of the dynamic substructure method

Researches show that there exists a modal transformation matrix which is similar to that of the fixed interface method in form but slightly different from that in content. Using the constrained modes and normal modes or Lanczos vectors, a new modal transformation developed from hybrid method is given, in which the interface forces are replaced by interface displacements. The first synthesis equations are assembled easily as with the fixed interface method, yet the interface coordinates can be further eliminated as with the free-interface method. The new method unifies the fixed- and free-interface methods as well as hybrid method,permitting easy implementation of the Lanczos vectors to replace normal modes so that only static modes are calculated. The dynamic substructure method thus becomes more flexible and efficient. This reveals the interrelations of various dynamic substructure methods. New formula for modal analysis and modal synthesis are developed and three examples are given for illustration.Supported by Doctoral Training Foundation of State Education Commission and NNSF of China.     

Flow-excited resonance of trapped modes of ducted shallow cavities

The self-excitation mechanism of the acoustic diametral modes of an axisymmetric internal cavity–duct system is studied for a Mach number range up to 0.4. The effect of cavity dimensions on the excitation mechanism is investigated experimentally and numerically. Experiments are conducted on three cavity depths and six cavity lengths for each depth. Numerical simulations of the mode shapes are also performed to determine the effect of cavity dimensions on the particle velocity field of the diametral modes. For all the tested configurations, the diametral modes are strongly excited at relatively low Mach numbers (as low as 0.1). The pulsation amplitude at resonance is found to increase as the cavity becomes shorter or deeper, relative to the main pipe diameter. The test results provide new insights into the excitation mechanism of diametral modes, the effect of the cavity length to depth ratio on the Strouhal numbers of acoustic resonances caused by various shear-layer modes of the cavity, and into the effect of the particle velocity field of the acoustic modes on the mode selectivity mechanism which determines the dominant acoustic mode during resonance.     

Non-linear normal modes,invariance, and modal dynamics approximations of non-linear systems

Non-linear systems are here tackled in a manner directly inherited from linear ones, that is, by using proper normal modes of motion. These are defined in terms of invariant manifolds in the system's phase space, on which the uncoupled system dynamics can be studied. Two different methodologies which were previously developed to derive the non-linear normal modes of continuous systems — one based on a purely continuous approach, and one based on a discretized approach to which the theory developed for discrete systems can be applied-are simultaneously applied to the same study case-an Euler-Bernoulli beam constrained by a non-linear spring-and compared as regards accuracy and reliability. Numerical simulations of pure non-linear modal motions are performed using these approaches, and compared to simulations of equations obtained by a classical projection onto the linear modes. The invariance properties of the non-linear normal modes are demonstrated, and it is also found that, for a pure non-linear modal motion, the invariant manifold approach achieves the same accuracy as that obtained using several linear normal modes, but with significantly reduced computational cost. This is mainly due to the possibility of obtaining high-order accuracy in the dynamics by solving only one non-linear ordinary differential equation.     

On simplified solutions of contact problems for elastic foundations

In the paper, a method of averaging displacements in a circular area lying in a linearly elastic transversally isotropic foundation is developed for the four modes of motion: vertical, horizontal, rocking and torsional. The corresponding formulae are constructed in a general form which does not depend on the kind of Green functions. For vertical and horizontal modes, uniform load distribution are applied and the simple integral mean is considered, whereas for rotational modes the load proportional to the distance from an axis of rotation is used, and angles of rotation for individual points are averaged with weight of the distance squared. Along with the case of equal radii of circles of loading and averaging, the case of different radii is studied, which allows one to consider contact problems for embedded axisymmetric foundations having the radius varying with depth. As examples the following contact problems are studied: static stiffness for a cone embedded in a homogeneous isotropic half-space in vertical motion, and dynamic stiffness for a disk on a layer resting on a homogeneous half-space for four modes of motion. Comparisons with the corresponding exact solutions are carried out.     

Torsional-mode transients in variable-thickness shells of revolution

For thin shells of revolution the existence of torsional-vibration modes, uncoupled from bending and extensional modes, has been established[1]. Here a linear second-order differential equation for the uncoupled torsional stress mode is obtained and its solution for impact loading of shells is sought. The mode-superposition method which utilizes the natural modes of vibration predicted by elementary theory, is, in general, not satisfactory for sharp impact loading as many modes are often required for convergence. Hence we employ two novel techniques for solving the impact problems. Firstly a formal asymptotic procedure, based on extensions to geometrical optics, is employed to generate asymptotic wavefront expansions. Rigorous justifications for this formal technique are provided in an appendix. Secondly a transform technique whereby solutions are sought in terms of Bessel functions is discussed and applied to particular impact loading problems. The Bessel function solutions found here can be used to determine the natural frequencies of the shells. Shells both finite and infinite in extent are discussed and reflections at a stress-free end are examined.     

Finite element computation of trapped and leaky elastic waves in open stratified waveguides

Elastic guided waves are of interest for inspecting structures due to their ability to propagate over long distances. In numerous applications, the guiding structure is surrounded by a solid matrix that can be considered as unbounded in the transverse directions. The physics of waves in such an open waveguide significantly differs from a closed waveguide, i.e. for a bounded cross-section. Except for trapped modes, part of the energy is radiated in the surrounding medium, yielding attenuated modes along the axis called leaky modes. These leaky modes have often been considered in non destructive testing applications, which require waves of low attenuation in order to maximize the inspection distance. The main difficulty with numerical modeling of open waveguides lies in the unbounded nature of the geometry in the transverse direction. This difficulty is particularly severe due to the unusual behavior of leaky modes: while attenuating along the axis, such modes exponentially grow along the transverse direction. A simple numerical procedure consists in using absorbing layers of artificially growing viscoelasticity, but large layers may be required. The goal of this paper is to explore another approach for the computation of trapped and leaky modes in open waveguides. The approach combines the so-called semi-analytical finite element method and a perfectly matched layer technique. Such an approach has already been successfully applied in scalar acoustics and electromagnetism. It is extended here to open elastic waveguides, which raises specific difficulties. In this paper, two-dimensional stratified waveguides are considered. As it reveals a rich structure, the numerical eigenvalue spectrum is analyzed in a first step. This allows to clarify the spectral objects calculated with the method, including radiation modes, and their dependency on the perfectly matched layer parameters. In a second step, numerical dispersion curves of trapped and leaky modes are compared to analytical results.