Anti-symmetric waves in pre-stressed imperfectly bonded incompressible elastic layered composites

The effect of an imperfect interface on the dispersive behavior of in-plane time-harmonic symmetric waves in a pre-stressed incompressible symmetric layered composite, was analyzed recently by Leungvichcharoen and Wijeyewickrema (2003). In the present paper the corresponding case for time harmonic anti-symmetric waves is considered. The bi-material composite consists of incompressible isotropic elastic materials. The imperfect interface is simulated by a shear-spring type resistance model, which can also accommodate the extreme cases of perfectly bonded and fully slipping interfaces. The dispersion relation is obtained by formulating the incremental boundary-value problem and using the propagator matrix technique. The dispersion relations for anti-symmetric and symmetric waves differ from each other only through the elements of the propagator matrix associated with the inner layer. The behavior of the dispersion curves for anti-symmetric waves is for the most part similar to that of symmetric waves at the low and high wavenumber limits. At the low wavenumber limit, depending on the pre-stress for perfectly bonded and imperfect interface cases, a finite phase speed may exist only for the fundamental mode while other higher modes have an infinite phase speed. However, for a fully slipping interface in the low wavenumber region it may be possible for both the fundamental mode and the next lowest mode to have finite phase speeds. For the higher modes which have infinite phase speeds in the low wavenumber region an expression to determine the cut-off frequencies is obtained. At the high wavenumber limit, the phase speeds of the fundamental mode and the higher modes tend to the phase speeds of the surface wave or the interfacial wave or the limiting phase speed of the composite. The bifurcation equation obtained from the dispersion relation yields neutral curves that separate the stable and unstable regions associated with the fundamental mode or the next lowest mode. Numerical examples of dispersion curves are presented, where when the material has to be prescribed either Mooney–Rivlin material or Varga material is assumed. The effect of imperfect interfaces on anti-symmetric waves is clearly evident in the numerical results.     

Time-harmonic wave propagation in a pre-stressed compressible elastic bi-material laminate

The dispersive behaviour of time-harmonic waves propagating along a principal direction in a perfectly bonded pre-stressed compressible elastic bi-material laminate is considered. The dispersion relation which relates wave speed and wavenumber is obtained by formulating the incremental boundary value problem and the use of the propagator matrix technique. At the low wavenumber limit, depending on the pre-stress, both the fundamental mode and the next lowest mode may have finite phase speeds. For the higher modes which have infinite phase speeds in the low wavenumber region, an expression to determine the cut-off frequencies is obtained. At the high wavenumber limit, the phase speeds of the fundamental mode and higher modes tend to phase speeds of the surface wave, the interfacial wave or the limiting phase speed of the composite. For numerical examples, either a two-parameter compressible neo-Hookean material or a two-parameter compressible Varga material is assumed.     

Wave propagation along a non-principal direction in a compressible pre-stressed elastic layer

The dispersive behavior of small amplitude waves propagating along a non-principal direction in a pre-stressed, compressible elastic layer is considered. One of the principal axes of stretch is normal to the elastic layer and the direction of propagation makes an angle θ with one of the in-plane principal axes. The dispersion relations which relate wave speed and wavenumber are obtained for both symmetric and anti-symmetric motions by formulating the incremental boundary value problem for a general strain energy function. The behavior of the dispersion curves for symmetric waves is for the most part similar to that of the anti-symmetric waves at the low and high wavenumber limits. At the low wavenumber limit, depending on the pre-stress and propagation angle, it may be possible for both the fundamental mode and the next lowest mode to have finite phase speeds, while other higher modes have an infinite phase speed. At the high wavenumber limit, the phase speeds of the fundamental mode and the higher modes tend to the Rayleigh surface wave speed and the limiting wave speeds of the layer, respectively. Numerical results are presented for a Blatz–Ko material and the effect of the propagation angle is clearly illustrated.     

Evolution of a long-wave train in loss of stability of a phase interface in geothermal systems

The transition to instability of phase interfaces in geothermal systems when a water stratum overlies a steam stratum and the most unstable mode corresponds to zero wavenumber is considered. The nonlinear Kolmogorov-Petrovskii-Piskunov equation describing the evolution of a narrow strip of weakly unstable modes is obtained. This equation is an analog of the well-known Ginzburg-Landau equation corresponding to the case of destabilization of modes with finite wavenumbers. It is shown that in the neighborhood of the critical points there exist two locations of the plane phase interface which coincide at the instant at which the instability threshold is reached and then disappear.     

正交各向异性板中非主轴方向的Lamb波

基于线性三维弹性理论,采用勒让德正交多项式展开法,推导了波沿正交各向异性材料非主轴方向传播时的Lamb波耦合波动方程,并对耦合波动方程进行了数值求解。为验证该方法的适用性和正确性,首先将此方法应用于各向同性材料,并与已知的数据结果进行了比较;然后以单向纤维增强复合材料为例,计算了耦合Lamb波沿不同的非主轴方向传播时的相速度频散曲线,并分别研究了传播方向改变时低阶模态Lamb波和高阶模态Lamb波频散特性的变化。最后,针对潜在用于各向异性复合材料结构健康监测的耦合Lamb波低阶模态,给出了其在不同传播方向时的相速度分布和群速度分布。同时,结合低阶模态Lamb波的位移分布特性和材料的各向异性特点,阐释了S0模态对波的传播方向变化最为敏感的原因。     

Wave propagation analysis in 2D nonlinear hexagonal periodic networks based on second order gradient nonlinear constitutive models

We analyze the propagation of nonlinear waves in homogenized periodic nonlinear hexagonal networks, considering successively 1D and 2D situations. Wave analysis is performed on the basis of the construction of the effective strain energy density of periodic hexagonal lattices in the nonlinear regime. The obtained second order gradient nonlinear continuum has two propagation modes: an evanescent subsonic mode that disappears after a certain wavenumber and a supersonic mode characterized by an increase of the frequency with the wavenumber. For a weak nonlinearity, a supersonic mode occurs and the dispersion curves lie above the linear dispersion curve (v_p =v_p0). For a higher nonlinearity, the wave changes from a supersonic to an evanescent subsonic mode at s=0.7 and the dispersion curves drops below the linear case and vanish for certain values of the wavenumber. An important decrease in the frequency occurs for both subsonic and supersonic modes when the lattice becomes auxetic, and the longitudinal and shear modes become very close to each other. The influence of the lattice geometrical parameters of the lattice on the dispersion relations is analyzed.     

Analysis of guided waves with a nonlinear boundary condition caused by internal resonance using the method of multiple scales

We theoretically investigated the cumulative nonlinear guided waves caused by internal resonance, using the method of multiple scales (MMS), which can construct better approximations to the solutions of perturbation problems. In this study, we consider nonlinearity only on the boundary instead of material nonlinearity or geometric nonlinearity. We showed nonlinear effects on the amplitudes of a lower mode and a higher mode depending on the propagation length. Also, we examined effects of wavenumber detuning from a phase matching condition of the two modes. If the wavenumber detuning is exactly equal to zero, the mechanical energy of the lower mode is transferred through nonlinear coupling to the energy of the higher mode, unilaterally. However, if a wavenumber detuning is not equal to zero, amplitude of the two modes change in a cyclic fashion during wave propagation. The amount of this amplitude variation and its cycle length are determined by the eigenfunctions of the two modes, the nonlinear parameter and the wavenumber detuning.     

二元复合材料板的自由振动: 半解析法

二元复合材料板是超材料板结构中常见的单元之一. 针对由材料参数相差两个量级的基体和嵌入体组成的二元复合材料板, 提出结构自由振动的半解析模型, 并对其振动特性进行了研究. 基于区域分解法和二元材料的分布, 将二维平板分解成两个子区域. 通过在振型函数中附加区域试函数, 来描述复合材料板面内刚度突变引起局部位移和转角的非光滑性. 基于二元复合材料板的基本边界条件和两子区连接处的变形协调条件, 构造了新的振型函数. 基于经典薄板理论, 利用带特殊试函数的里兹法, 求得不同几何构型下二元复合材料板的固有频率和振型, 并研究了嵌入体的尺寸和位置对结构振动特性的影响规律. 通过收敛分析并与有限元仿真结果对比, 验证了本文方法的准确性. 研究结果表明: 传统的全局试函数在分析具有振动局部化的模态时会得到不准确的结果, 而附加区域试函数可以显著提高里兹法的收敛速度以及结果的准确性; 嵌入体位置对低阶固有频率的作用不明显, 却能显著改变低阶振型节线的分布和振动局部化发生的区域.      

气液两相流压力波传播速度研究

将双流体模型用于绝热无相的管道气液两相流,依据小扰动线化分析原理,导出了压力波波数Ｋ方程通过对不同空隙率下肉体上压力波小随角频率变化的计算,研究了虚拟质量力和狭义相间阻力对压力波波速及其人色散性的影响。对泡状流和弹状流压力波波速的计算结果与前人的测量结果作了比较,两者符合良好。     

Lamb mode stresses as means of identifying voids in the adhesive zone of bilaminates

Lamb wave technique has emerged as a reliable tool in the nondestructive testing of laminated plates. Some current studies to identify the specific Lamb modes that can characterize different kinds of defects in layered plates using Lamb waves have shown that the modes for which high stresses and low displacements occur in the interface indicate the presence of defects like pores or voids whereas the modes for which the displacements are high show the presence of harder inclusions. In this context this paper tests an earlier analytical model developed to facilitate NDT of porosity in the adhesive zone of bilaminates.The model tested treats the pore infested thin adhesive region as a linear elastic material with voids (LEMV). For certain parametric values of the LEMV adhesive layer the influence of these voids on dispersion and stresses carried by the first few Lamb modes in glass/glue/glass (G/g/G) bilaminate is traced in the range 0–10 MHz. The frequency–phase velocity points experimentally obtained by Kundu and Maslov are seen to fall very close to the present dispersion. The stresses traced using the present model in G/g/G plate at these experimentally tallied points show an easily discernable rise in the central region of adhesive, as observed by Kundu and Maslov.The model appears to be useful as a good first approximation to detect voids in adhesive zone of composite structural elements.     

Anti-plane shear waves in a fibre-reinforced composite with a non-linear imperfect interface

The propagation of non-linear elastic anti-plane shear waves in a unidirectional fibre-reinforced composite material is studied. A model of structural non-linearity is considered, for which the non-linear behaviour of the composite solid is caused by imperfect bonding at the “fibre–matrix” interface. A macroscopic wave equation accounting for the effects of non-linearity and dispersion is derived using the higher-order asymptotic homogenisation method. Explicit analytical solutions for stationary non-linear strain waves are obtained. This type of non-linearity has a crucial influence on the wave propagation mode: for soft non-linearity, localised shock (kink) waves are developed, while for hard non-linearity localised bell-shaped waves appear. Numerical results are presented and the areas of practical applicability of linear and non-linear, long- and short-wave approaches are discussed.     

均匀来流条件下并行排列旗帜耦合运动模式的实验

利用高速摄影技术在低速风洞中记录了不同间距并行排列的两个旗帜在不同来流速度中的耦合运动。利用自编的时间-空间演化图像处理软件分析总结了旗帜的耦合运动模式以及旗帜摆动振幅、频率和St数的变化规律。实验结果显示,随着排列间距和来流速度的变化,两旗帜可能以静止、同向摆动、反向摆动和过渡状态这四种不同的模态耦合运动。两旗帜同向摆动时摆动频率明显低于单个旗帜在相同来流中的值,反向摆动时情况相反。在过渡状态中两旗帜摆动的振幅交替增减并且运动中同时包含有两个频率,而同向和反向摆动都是单频率的运动。     

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.     

Elastic waves guided by a material interface

The propagation of interfacial small-amplitude waves along a rectilinear thin film separating two pre-stressed, incompressible, elastic media is addressed. The film is modelled as a material surface possessing its own mass density and normal and flexural stiffnesses. It is shown that these features induce dispersion as the obtained secular equations are polynomials of the second degree in the wavenumber when bending stiffness is absent (membrane-like interface), and of the fourth degree otherwise (plate-like interface). In both case, beyond the modified Stoneley mode, a bending mode for the interface, an additional propagating wave can exist, with amplitude polarized along the interface (extensional mode). The associated bifurcation problem is analyzed with focus on the effects of compressive residual forces at the interface. The buckling strain of a compressed metal layer embedded in an elastomeric medium is computed also with an exact approach, to provide the range of validity of the proposed simplified model of material interface.     

Finite-amplitude Love waves in a pre-stressed compressible elastic half-space with a double surface layer

The dispersive behavior of finite-amplitude time-harmonic Love waves propagating in a pre-stressed compressible elastic half-space overlaid with two compressible elastic surface layers of finite thickness is investigated. The half-space and layers are made of different pre-stressed compressible neo-Hookean materials. The dispersion relation which relates wave speed and wavenumber is obtained in explicit form. Results for the energy density and energy flux of the waves are also presented. The special case where the interfaces between the layers and the half-space are principal planes of the left Cauchy–Green deformation tensor is also investigated. Numerical results are presented showing the variation of the Love wave speed with the pre-stress and the propagation angle.     

Thermo-elastic dynamic instability (TEDI) in frictional sliding of two elastic half-spaces

In some simplified 1D models, we recently studied the coupling of TEI (thermoelastic instability) and DI (dynamic instability), finding that thermal effects can render unstable the otherwise neutrally stable natural elastodynamic modes of the system, giving rise to a new family of instability which we called TEDI.Here, we study the general case of two sliding elastic half-planes, finding again a relatively weak coupling between thermal and dynamic effects, and the general family of instability TEDI class is found to modify both the otherwise separated TEI and DI classes. The growth factor, the phase velocity and the migrating speeds of the perturbations are wavelength-dependent, and it is difficult to give a complete picture given the high number of materials’ parameters, and the dependence on speed, friction coefficient, and the underlying uniform pressure. However, a set of results are given for “large” and “small” mismatch of shear wave speeds in the materials, and as a function of (i) friction coefficient; (ii) sliding speed V₀; (iii) wavenumber parameter γ. In the case of small mismatch, generalized Rayleigh waves exists already under frictionless conditions, the critical f for instability is zero. DI dominates over TEI typically for large wavenumbers, where the growth factors increase without limit and hence become eventually meaningless, requiring regularizations for example with rate-state dependent friction laws. TEI growth factors vice versa have a maximum at a certain wavenumber and therefore are always well posed. Larger coupling effects are noticed for two materials with large mismatch, but significantly only for sliding speeds comparable with the wave speed. In general, TEI growth factors increase with speed, whereas DI growth factors increase with speed for similar materials and decrease when the mismatch between materials is large.     

Wave propagation in layer with two preferred directions

This article is concerned with overall or macroscopic properties of a composite material with no distinction made between the fibres and the matrix which they are embedded in. All the properties with dimensions larger than the fibre diameter and spacing are regarded as averaged over a volume of material. The systems of particular interest here are in the fibre reinforced composites with the fibres being very much stiffer and stronger than the matrix.Laminated plates of fibre-reinforced material are often fabricated from prepreg tapes, laid up according to some specific arrangement of fibre orientation and then bonded together. An angle-ply laminate is formed by alternating plies so that the families in adjacent laminas are inclined by angle ϕ and −ϕ to given direction alternately. The process of fabricating a multilayered plate of this material gives rise to a laminate in which the plies are separated by resin rich layer, and when this layer is thin enough that its thickness is negligible it may be regarded as plate reinforced by two families of fibres. Problems shall be considered in three dimensions, but attention shall be restricted to linear elasticity theory. The plate under consideration is reinforced by two mechanically equivalent families of fibres, but with no other preferred directions, so that it is locally orthotropic with respect to the plane of the fibres and to the two planes that orthogonally bisect the fibres.In this article linear elastic stress–strain relation is employed to derive dispersion curves for plane harmonic waves propagating in a plate of finite thickness but of infinite lateral extent. Attention is restricted to waves propagating in the plane parallel to stress free plate faces where waves travelling at any angle to one of the families of very strong fibres are examined. The dispersion equations, relating the phase velocity to the wavelength, are obtained. The fundamental modes are examined for symmetric as well as for anti-symmetric deformations. This leads to full understanding of displacement field as well as stress field.     

Surface wave transference in a piezoelectric cylinder coated with reinforced material

The present study deals with the propagation of a polarized shear horizontal(SH)wave in a pre-stressed piezoelectric cylinder circumscribed by a self-reinforced cylinder.The interface of the two media is assumed mechanically imperfect.For obtaining the dispersion relation,the mathematical formulation has been developed and solved by an analytical treatment.The effects of various parameters,i.e.,the thickness ratio,the imperfect interface,the initial stress,the reinforcement,and the piezoelectric and dielectric constants,on the dispersion curve are observed prominently.The dispersion curves for different modes have been also plotted.The consequences of the study may be used for achieving optimum efficiency of acoustic wave devices.     

Frictional sliding modes along an interface between identical elastic plates subject to shear impact loading

Frictional sliding along an interface between two identical isotropic elastic plates under impact shear loading is investigated experimentally and numerically. The plates are held together by a compressive stress and one plate is subject to edge impact near the interface. The experiments exhibit both a crack-like and a pulse-like mode of sliding. Plane stress finite element calculations modeling the experimental configuration are carried out, with the interface characterized by a rate and state dependent frictional law. A variety of sliding modes are obtained in the calculations depending on the impact velocity, the initial compressive stress and the values of interface variables. For low values of the initial compressive stress and impact velocity, sliding occurs in a crack-like mode. For higher values of the initial compressive stress and/or impact velocity, sliding takes place in a pulse-like mode. One pulse-like mode involves well-separated pulses with the pulse amplitude increasing with propagation distance. Another pulse-like mode involves a pulse train of essentially constant amplitude. The propagation speed of the leading pulse (or of the tip of the crack-like sliding region) is near the longitudinal wave speed and never less than times the shear wave speed. Supersonic trailing pulses are seen both experimentally and computationally. The trends in the calculations are compared with those seen in the experiments.