Love waves in a two-layered piezoelectric/elastic composite plate with an imperfect interface

Based on the shear spring model, the propagation of Love wave in two-layered piezoelectric/elastic composite plates under the influence of interfacial defect is investigated. The piezoelectric layer is electrically shorted at both top and bottom surfaces. The wave form solutions of the piezoelectric and elastic layers are obtained, and the dispersion equation is derived by subjecting the boundary conditions and the continuity conditions to the obtained wave form solutions. Numerical results are performed for PZT4/aluminum composite plate. The phase velocities and the mode shapes of mechanical displacement and electric potential are illustrated graphically. The results show that both the interfacial defect and the thickness ratio between the piezoelectric and elastic layers have significant effect on the propagation characteristics of Love wave. One important feature is observed that the interfacial defect always decreases the phase velocities.     

双介质耦合刚性基弹性层平面应变型导波模式及界面散射能量分配

过渡段动力稳定性问题已成为制约400 km/h及以上高铁路基设计的关键难题,亟需从波动和能量的角度探究由基础非均匀引发的线路系统动力响应放大机理.文章将轨下基础简化为上表面自由、底端固定的刚性基弹性层,将高铁过渡段车致弹性波传播问题提炼为非均匀介质刚性基弹性层中波的散射问题,建立双介质耦合刚性基弹性层平面应变模型,优化该类波导结构频散方程在复平面求根方法,并结合岩土类介质特征展开刚性基弹性层频散分析,以明确其多模式导波特性及散射能量分配,最后,围绕弹性层厚度、刚度比等影响因素开展对比分析.结果表明:刚性基弹性层各模式导波均具有截止频率,弹性层厚度越小,杨氏模量越大,各阶导波模式的截止频率越高;入射波在双介质刚性基弹性层发生散射后,透射场基阶模式导波会占据主体能量,随着高阶导波模式被逐一激发,反射场及透射场高阶模式能量占比会在全频率范围呈现“此消彼长”状态;交换两侧弹性层材料,改变弹性层厚度及两弹性层刚度比不会显著改变能量分布规律,但总体来看,能量更易集中在较软侧弹性层中,各模式导波在激发初始频段会更为活跃,可分配到更多能量.     

Structure of weak shock waves in 2-D layered material systems

Shock waves in homogeneous materials in the absence of phase transitions are understood to have a one-wave structure. However, upon loading of a layered heterogeneous material system a two-wave structure is obtained––a leading shock front followed by a complex pattern that varies with time. This dual shock-wave pattern can be attributed to material architecture through which the shock wave propagates, i.e. the impedance (and geometric) mismatch present at various length scales, and nonlinearities arising from material inelasticity and failure.The objective of the present paper is to provide a better understanding of the role of material architecture in determining the structure of weak shock waves in 2-D layered material systems. Normal plate-impact experiments are conducted on 2-D layered material targets to obtain both the precursor decay and the late-time dispersion. The particle velocity at the free surface of the target plate is measured by using a multi-beam VALYN VISAR. In order to understand the effects of layer thickness and the distance of wave propagation on elastic precursor decay and late-time dispersion several different targets with various layer and target thicknesses are employed. Moreover, in order to understand the effects of material inelasticity both elastic–elastic and elastic–viscoelastic bilaminates are utilized.The results of these experiments are interpreted by using asymptotic techniques to analyze propagation of acceleration waves in 2-D layered material systems. The analysis makes use of the Laplace transform and Floquet theory for ODE’s with periodic coefficients [Asymptotic solutions for wave propagation in elastic and viscoelastic bilaminates. In: Developments in Mechanics, Proceedings of the 14th Mid-Eastern Mechanics Conference, vol. 26, no. 8, pp. 399–417]. Both wave-front and late-time solutions for step-pulse loading on layered half-space are compared with the experimental observations. The results of the study indicate that the structure of acceleration waves is strongly influenced by impedance mismatch of the layers constituting the laminates, density of interfaces, distance of wave propagation, and the material inelasticity.     

On the retarded potentials of inhomogeneous ellipsoids and sources of arbitrary shapes in the three-dimensional infinite space

We analyze the infinite space solutions of the three-dimensional inhomogeneous wave equation (the ‘retarded potentials’ or ‘causal propagators’) for ellipsoidal sources and for sources of arbitrary shapes. The ‘short-time characteristics’ of the retarded potential for a spatially inhomogeneous source density of δ-shaped time profile is considered. It is found that, the short-time characteristics is governed by the spatial inhomogeneity of the source density in the immediate vicinity of a spacepoint.Surface integral representations are derived for spatial inhomogeneous source regions of ellipsoidal symmetry. For spherical sources these integral representations yield closed form solutions for the retarded potentials. We find that the wave field inside a spherical source consists of an incoming and outgoing spherical wave package, whereas the external wave field consists of an outgoing spherical wave package only. Characteristic runtime and superposition effects are discussed. Moreover, a numerical technique based on Gauss quadrature is applied to generate the wave field for a cubic source. The integral representations derived for the retarded potentials of inhomogeneous ellipsoidal sources are consistent with results previously derived by the authors for the Helmholtz potentials of homogeneous ellipsoids and ellipsoidal shells [Michelitsch, T.M., Gao, H., Levin, V.M., 2003. On the dynamic potentials of ellipsoidal shells. Q. J. Mech. Appl. Math. 56 (4), 629]. The derived solutions are crucial for many problems of wave propagation and diffraction theory as they may occur in materials science. As an example we give a formulation for the solution of the retarded Eshelby inclusion problem due to spatially and temporally varying eigenfields in the elastic isotropic infinite medium.     

Multi-directional and multi-time step absorbing layer for unbounded domain

Variant techniques are proposed for reproducing the elastic wave propagation in an unbounded medium such as the infinite elements, the absorbing boundary conditions or the perfect matched layers. Here, a simplified approach is adopted by considering absorbing layers characterized by the viscous Rayleigh matrix as studied by Semblat et al. [16] and Rajagopal et al. [14]. Here, further improvements to this procedure are provided. First, we start by establishing the strong form for the elastic wave propagation in a medium characterized by the Rayleigh matrix. This strong form will be used for deriving optimal conditions for damping out in the most efficient way the incident waves while minimizing the spurious reflected waves at the interface between the domain of interest and the Rayleigh damping layer. A procedure for designing the absorbing layer is proposed by targeting a performance criterion expressed in terms of logarithmic decrement of the wave amplitude in the layer thickness. Second, the GC subdomain coupling method, proposed by Combescure and Gravouil [9], is introduced for enabling the choice of any Newmark time integration schemes associated with different time steps depending on subdomains. When wave propagation is predicted by an explicit time integrator, the subdomain strategy is of great interest because it enables a different time integrator for the absorbing layer to be adopted. An external coupling software, based on the GC method, is used to carry out multi=time step explicit/implicit co-computations, making interact in time an explicit FE code (Europlexus) for the domain of interest, with an implicit FE code (Cast3m) handling the absorbing boundary layers. The efficiency of the approach is shown in 1D and 2D elastic wave propagation problems.     

Finite-amplitude shear horizontal waves propagating in a pre-stressed layer between two half-spaces

The propagation of finite-amplitude time-harmonic shear horizontal waves, in a pre-stressed compressible elastic layer of finite thickness embedded between two identical compressible elastic half-spaces, is investigated. This is accomplished by combining finite-amplitude linearly polarized inhomogeneous transverse plane wave solutions in the half-spaces and finite-amplitude linearly polarized unattenuated transverse plane wave solutions in the layer. The layer and half-spaces are made of different pre-stressed compressible neo-Hookean materials. The dispersion relation which relates wave speed and wavenumber is obtained in explicit form. The special case where the interfaces between the layer and the half-spaces are principal planes of the left Cauchy–Green deformation tensor is also investigated. Numerical results are presented showing the variation of the shear horizontal wave speed with the pre-stress and the propagation angle.     

Influence of prestress on the velocities of plane waves propagating normally to the layers of nanocomposites

The paper is concerned with longitudinal and transverse waves propagating at a right angle to the layers of a nanocomposite material with initial (process-induced residual) stresses. The composite consists of alternating layers of two dissimilar materials. The materials are assumed nonlinearly elastic and described by the Murnaghan potential. For simulation of wave propagation, a problem is formulated within the framework of the three-dimensional linearized theory of elasticity for finite prestrains. It is established that the relative velocities of waves depend linearly on small prestresses. In some composite materials, however, the thicknesses of the layers may be in a ratio such that the wave velocities are independent of the prestress level __________ Translated from Prikladnaya Mekhanika, Vol. 42, No. 7, pp. 3–22, July 2006.     

各向异性展状介质中弹性波的传播矩阵解法

基于广义胡克定律及混和变量弹性波方程，解析求得各层介质内位移、应力传递矩阵，给出了直角坐标系下各向异性层状介质中弹性波的传播矩阵解法．该方法适用于非轴对称各向异性和点源作用，较好地解决了数值计算中有效数字精度损失问题．数值结果表明，计算效率、准确性及稳定性均较好．     

Effect of negative permeability on elastic wave propagation in magnetoelastic multilayered composites

With the advent of left-handed magnetic materials, it is desirable to develop high-performance wave devices based on their novel properties of wave propagation. This letter reports the special properties of elastic wave propagation in magnetoelastic multilayered composites with negative permeability as comparecd to those in counterpart structures with positlve permeability. These novel properties of elastic waves are discerned from the diversified dispersion curves, which represent the propagation and attenuation characteristics of elastic waves. To compute these dispersion curves, the method of reverberation-ray matrix is extended for the analysis of elastic waves in magnctoelastic multilayered composites. Although only the results of a single piezomagnetic and a binary magnetoelastic layers with mechanically free and magnetically short surfaces as well as pelrfect interface are illustrated in the numerical examples, the analysis is applicable lo magnetoelastic multilayered structures with other kinds of boundaries/interfaces.     

Viscoelastic shear horizontal wave in graded and layered plates

In this paper, viscoelastic shear horizontal (SH) wave propagation in functionally graded material (FGM) plates and laminated plates are investigated. The controlling differential equation in terms of displacements is deduced based on the Kelvin–Voigt viscoelastic theory. The SH wave characteristics is controlled by two elastic constants and their corresponding viscous coefficients. By the Legendre polynomial series method, the asymptotic solutions are obtained. In order to verify the validity of the method, a homogeneous plate is calculated to make a comparison with available data. Through three different graded plates, the influences of gradient shapes on dispersion and attenuation are discussed. The viscous effects on the displacement and stress shapes are illustrated. The different boundary conditions are analyzed. The influential factors of the viscous effect are analyzed. Finally, two multilayered (two layer and five layer) viscoelastic plates that are composed of the same material volume fraction are calculated to show their differences from the graded plate.     

Comparison of waveguide properties of curved versus straight planar elastic layers

Dispersion equations are solved for the in-plane and anti-plane wave propagation in planar elastic layer with constant curvature. The classical Lamé formulation of displacements via elastic potentials is applied and appropriate simplifications are employed. The dispersion diagrams in each case are compared with their counterparts for a straight layer, e.g., the classical Rayleigh–Lamb solution. The curvature-induced symmetry-breaking effects are investigated for layers with symmetric boundary conditions. The role of curvature is also investigated in the cases, when the boundary conditions are not symmetrical. The elementary Bernoulli–Euler theory is employed to analyze the wave guide properties of a curved planar elastic beam in its in-plane deformation. The validity range of the Bernoulli–Euler theory is assessed via comparison of dispersion diagrams.     

Study on wave dispersion characteristics of piezoelectric sandwich nanoplates considering surface effects

In this study, the wave propagation properties of piezoelectric sandwich nanoplates deposited on an orthotropic viscoelastic foundation are analyzed by considering the surface effects (SEs). The nanoplates are composed of a composite layer reinforced by graphene and two piezoelectric surface layers. Utilizing the modified Halpin-Tsai model, the material parameters of composite layers are obtained. The displacement field is determined by the sinusoidal shear deformation theory (SSDT). The Euler-Lagrange equation is derived by employing Hamilton’s principle and the constitutive equations of piezoelectric layers considering the SEs. Subsequently, the nonlocal strain gradient theory (NSGT) is used to obtain the equations of motion. Next, the effects of scale parameters, graphene distribution, orthotropic viscoelastic foundation, and SEs on the propagation behavior are numerically examined. The results reveal that the wave frequency is a periodic function of the orthotropic angle. Furthermore, the wave frequency increases with the increase in the SEs.

     

Features of plane wave propagation along the layers of a prestrained nanocomposite

Features of the propagation of longitudinal and transverse plane waves along the layers of nanocomposites with process-induced initial stresses are studied. The composite has a periodic structure: it is made by repeating two highly dissimilar layers. The layers exhibit nonlinear elastic behavior in the range of loads under consideration. A Murnaghan-type elastic potential dependent on the three invariants of the strain tensor is used to describe the mechanical behavior of the composite constituents. To simulate the propagation of waves, finite-strain theory is used for developing a problem statement within the framework of the three-dimensional linearized theory of elasticity assuming finite initial strains. The dependence of the relative velocities of longitudinal and transverse waves on two components of small initial stresses in each layer and on the volume fraction of the constituents is studied. It is established that there are thickness ratios of layers in some nanocomposites such that the wave velocities are independent of the initial stresses and equal to the respective wave velocities in composites without initial stresses __________ Translated from Prikladnaya Mekhanika, Vol. 43, No. 4, pp. 3–26, April 2007.     

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.     

Plane acoustic wave propagation through a composite of elastic and Kelvin–Voigt viscoelastic material layers

The problem of plane wave propagation through a plane composite layer of thickness h is considered. The composite consists of periodically repeated elastic and Kelvin–Voigt viscoelastic material layers, and all layers are either parallel or perpendicular to the incident wave front. Moreover, it is assumed that the thickness of each separate layer of the composite is much less than the acoustic wave length and the thickness h of the entire composite. We study the problem by using a homogenized model of the composite, which allows us to find the reflection and transmission factors and the variation in the sound intensity level as it propagates though the composite layer of thickness h.     

A DAMAGE ACCUMULATING MODELING OF FAILURE WAVES IN GLASS UNDER HIGH VELOCITY IMPACT

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｆａｉｌｕｒｅｗａｖｅｉｓａｓｐｅｃｉｆｉｃｆａｉｌｕｒｅｐｈｅｎｏｍｅｎｏｎｏｃｃｕｒｒｉｎｇｉｎｇｌａｓｓｍｅｄｉａｕｎｄｅｒｃｒｉｔｉｃａｌｉｍｐａｃｔｌｏａｄｉｎｇｗｉｔｈｈｉｇｈｓｔｒａｉｎｒａｔｅ (ｖ≥ 1 2ｋｍ/ｓ, ε≥ 1 0 5 /ｓ) .ＩｔｗａｓｆｉｒｓｔｏｂｓｅｒｖｅｄｉｎＫ1 9ｇｌａｓｓｐｌａｎａｒｉｍｐａｃｔｔｅｓｔｓｉｎ 1 991ｂｙＲａｓｏｒｅｎｏｖ[1]ｔｈａｔｔｈｅｒｅｅｘｉｓｔｓａｆａｉｌｕｒｅｌａｙｅｒｐｒｏｐａｇａｔｉｎｇａｔａｌｏｗｅｒｓｐｅ…     

Fully non-linear wave models in fiber-reinforced anisotropic incompressible hyperelastic solids

Many composite materials, including biological tissues, are modeled as non-linear elastic materials reinforced with elastic fibers. In the current paper, the full set of dynamic equations for finite deformations of incompressible hyperelastic solids containing a single fiber family are considered. Finite-amplitude wave propagation ansätze compatible with the incompressibility condition are employed for a generic fiber family orientation. Corresponding non-linear and linear wave equations are derived. It is shown that for a certain class of constitutive relations, the fiber contribution vanishes when the displacement is independent of the fiber direction.Point symmetries of the derived wave models are classified with respect to the material parameters and the angle between the fibers and the wave propagation direction. For planar shear waves in materials with a strong fiber contribution, a special wave propagation direction is found for which the non-linear wave equations admit an additional symmetry group. Examples of exact time-dependent solutions are provided in several physical situations, including the evolution of pre-strained configurations and traveling waves.     

Stress intensity factors for a moving cylindrical crack in a nonhomogeneous cylindrical layer in composite materials

Stresses are determined for a finite cylindrical crack that is propagating with a constant velocity in a nonhomogeneous cylindrical elastic layer, sandwiched between an infinite elastic medium and a circular elastic cylinder made from another material. The Galilean transformation is employed to express the wave equations in terms of coordinates that are attached to the moving crack. An internal gas pressure is then applied to the crack surfaces. The solution is derived by dividing the nonhomogeneous interfacial layer into several homogeneous cylindrical layers with different material properties. The boundary conditions are reduced to two pairs of dual integral equations. These equations are solved by expanding the differences in the crack surface displacements into a series of functions that are equal to zero outside the crack. The Schmidt method is then used to solve for the unknown coefficients in the series. Numerical calculations for the stress intensity factors were performed for speeds and composite material combinations.     

玄武岩中裂隙分布形式对声波传播的影响

复杂岩体含有大量的裂隙,这些裂隙尺寸及其分布形式等对弹性波传播都有很大的影响.本文加工了含单个裂隙、双裂隙和三个裂隙的玄武岩岩样单元对其进行组合,进行了25kHz、 50kHz、 400kHz、 600kHz和1000kHz 等5种频率的声波测试.通过考虑垂直或平行波传播方向的裂隙长度,来探索裂隙分布形式和不同裂隙长度对弹性波传播的影响,研究玄武岩的频散效应和波的衰减.结果表明:裂隙方向与波传播方向夹角对弹性波传播有很大的影响.当裂隙方向与波传播方向垂直时,散射效应最大;而当裂隙方向与波传播方向平行时,影响最小.上述结果可为理论模型和数值分析提供依据.     

A new 2D discrete model applied to dynamic crack propagation in brittle materials

In this work, a 2D discrete model (DM) applied to the dynamic crack propagation in brittle materials is developed and implemented. The proposed model is based on a particular discretization of Navier’s equations, presenting similarities to the Born model, with the advantage that the constants appearing in it are explicitly related to the elastic properties. This model overcomes the limitations in the choice of Poisson’s ratio present in other discrete models. Three numerical examples are presented to show the capability of this method in modelling wave propagation and dynamic fracture problems. The obtained results are in agreement with experimental and numerical results reported by other researchers.