Quadratically nonlinear torsional hyperelastic waves in a transversely isotropic cylinder: Primary analysis of evolution

Nonlinear torsional waves propagating along an infinite transversely isotropic cylinder are considered. The hyperelasticity of the cylinder is described by Guz’s generalization of the Murnaghan potential to quasi-isotropic materials. The evolution of the initial waveprofile in cylinders made of composite materials reinforced with micro-and nanoscale fibers is modeled numerically. It is shown that the transverse isotropy of this class of composites has a weak effect on the wave phenomena accompanying the propagation of a torsional wave. Three-dimensional plots of evolution are presented __________ Translated from Prikladnaya Mekhanika, Vol. 44, No. 5, pp. 32–44, May 2008.     

Propagation characteristics of thermoelastic waves in piezoelectric (6 mm class) rotating plate

The propagation of Lamb waves in a homogeneous, transversely isotropic (6 mm class), piezothermoelastic plate rotating with uniform angular velocity about normal to its boundary has been investigated. The generalized (non-classical) theories of thermoelasticity in contrast to Sharma and Pal [Sharma, J.N., Pal, M., 2004. Lamb wave propagation in transversely isotropic piezothermoelastic plate. J. Sound Vib. 270, 587–610] have been used to investigate the problem. The surfaces of the plate are subjected to stress free, thermally insulated/isothermal and electrically shorted boundary conditions. Secular equations for wave propagation modes in the plate are derived from a coupled system of governing partial differential equations of linear piezothermoelasticity. After obtaining the complex characteristic roots with the help of Descartes' algorithm, the transcendental secular equations have been solved by functional iteration numerical technique to compute phase velocity and attenuation coefficient. Finally, in order to illustrate the analytical development, numerical solution of secular equations is carried out for PZT-5A piezo-thermoelastic material. The corresponding simulated results of various physical quantities such as phase velocity, attenuation coefficients, specific loss factor of energy dissipation, thermo-mechanical coupling factor and relative frequency shifts have been presented graphically for both rotating and non-rotating plates for comparison purpose. There is a scope for extension of the present work to other classes of piezo/pyroelectric crystals. The study will be useful in design and construction of gyroscope, rotation sensors, temperature sensors and other pyro/piezoelectric surface acoustic wave (SAW) devices.     

Surface waves in a piezoelectric–semiconductor composite structure

The present work deals with the propagation of interfacial surface waves in a composite consisting of homogeneous, transversely isotropic, piezoelectric halfspace underlying a thin layer of non-piezoelectric semiconductor material. The mathematical model of the problem is depicted by partial differential equations of motion for the structure and boundary conditions to be satisfied at the interface and free surface of the composite. After obtaining formal wave solution of the model the secular equation that governs the propagation of surface waves in the considered composite structure has been derived in compact form. The numerical solution of secular equation is being carried out for the composites Si–CdSe, Ge–CdSe and Ge–PZT by employing functional iteration method along with irreducible Cardano method using MATLAB programming. The computer simulated results in respect of dispersion curves, attenuation coefficient and specific loss factor of energy dissipation are presented graphically for Si–CdSe composite to illustrate the analytical developments. We have extended our analysis to Ge–CdSe and Ge–PZT composites also. However, to avoid clustering of profiles and also to have clear understanding of the variations, the computer simulated values of phase velocity and attenuation coefficient are presented in tabular form for all three considered composite structures. This work may be useful for designing and construction of surface acoustic wave (SAW) devices and electronics industry.     

Design and analysis of nonlinear-transformation-based broadband cloaking for acoustic wave propagation

A coordinate-transformation method can be used to design invisibility cloaks for many types of waves, including acoustic waves. The traditional method for designing a cloak depends on a transformation from a virtual space to a physical space. Previous acoustic cloaks that are mainly designed with linear-transformation-based acoustics have drawbacks that acoustic wave trajectories in the cloaks cannot be controlled and tuned. This work uses a nonlinear mapping from a ray trajectory perspective to construct acoustic cloaks with tunable non-singular material properties. Use of a ray trajectory equation is a straightforward and alternate way to study propagation characteristics of different types of waves, which allows more flexibility in controlling the waves. A broadband cylindrical cloak for acoustic waves in an inviscid fluid is realized with layered non-singular, homogeneous, and isotropic materials based on a nonlinear transformation. Some advantages and improvements of the invisibility nonlinear-transformation cloak over a traditional linear-transformation cloak are analyzed. The invisibility capability of the nonlinear-transformation cloak can be tuned by adjusting a design parameter that is shown to have influence on the acoustic wave energy flowing into the region inside the cloak. Numerical examples show that the nonlinear-transformation cloak is more effective for making a domain undetectable by acoustic waves in an inviscid fluid and shielding acoustic waves from outside the cloak than the linear-transformation cloak in a broad frequency range. The methodology developed here can be used to design nonlinear-transformation cloaks for other types of waves.     

川藏公路地质环境与整治改建方案的思考

川藏公路由于地质环境复杂、建设标准低、后遗病害多,抗灾能力差,泥石流、滑坡、山崩、雪害、水毁等自然灾害频繁发生,公路阻车断道严重。国家投入巨资进行整治改建,并取得了明显的效果,但由于自然环境特殊、影响因素复杂,许多特大型、大型工程地质病害问题还没有可行、可靠的解决方案。本文通过分析川藏公路沿线的地质环境和灾害特点,总结历年整治改建和经验的教训,提出川藏公路建设的途径、可能达到的目标和应采用的原则。     

Elastic fields in rotating transversely isotropic media

Representation of elastic fields in terms of scalar functions, which permit reducing the problems of determining these fields to determining scalar potentials, are generalized to the case of transversely isotropic media rotating at a constant angular velocity. Relations for calculating the parameters of surface acoustic waves (SAW) propagating in a rotating transversely isotropic halfspace with various directions of the medium material symmetry axis with respect to the half-space surface are given.     

Dispersion in oceanic crust during earthquake preparation

Dispersion of Stoneley waves is studied in a sedimentary layer of ocean bottom resting over basaltic solid half space. Sedimentary layer is assumed a transversely isotropic poroelastic medium. Lower-most solid half-space is assumed to be embedded with vertically aligned saturated micro-cracks and behaves transversely isotropic to wave propagation.Frequency equation is obtained in the form of determinantal equation. Role of phase angle is eliminated by expressing slowness of waves in terms of phase velocity and elastic constants. Numerical solutions for phase velocity and group velocity are obtained for a particular model. Calculations are made for different depths of ocean and sediments. Effect of thickness and density of cracks on these velocities are observed.Special cases are discussed which represent the absence of ocean and sediments, in the model considered. Changes in dispersion are discussed during the stress accumulation in an earthquake preparation region.     

Propagation of wave at the boundary surface of transversely isotropic thermoelastic material with voids and isotropic elastic half-space

The purpose of this research is to study the effect of voids on the surface wave propagation in a layer of a transversely isotropic thermoelastic material with voids lying over an isotropic elastic half-space. The frequency equation is derived after developing a mathematical model for welded and smooth contact boundary conditions. The dispersion curves giving the phase velocity and attenuation coefficient via wave number are plotted graphically to depict the effects of voids and anisotropy for welded contact boundary conditions. The specific loss and amplitudes of the volume fraction field, the normal stress, and the temperature change for welded contact are obtained and shown graphically for a particular model to depict the voids and anisotropy effects. Some special cases are also deduced from the present investigation.     

多个椭圆柱波浪力的一种解析解

波浪在大尺寸结构表面产生不可忽略的散射波, 该散射波在多柱体体系中继续传播, 并在同体系中的其他柱体上产生高次散射波. 本文基于椭圆坐标系和绕射波理论首先推导了波浪作用下椭圆单柱体产生的散射波压力公式, 随后考虑该散射波在多柱体系中的传播, 将其视为第二次入射波, 推导出柱体上第二次散射波压力公式, 同理可以推导出高次散射波压力公式, 最后得到椭圆多柱体波浪力解析解, 并用数值解验证了本文解析方法的正确性. 本文以双柱体和四柱体体系为例, 分析了不同参数(波数、净距、波浪入射角度等)下, 高次散射波对柱体上波浪作用的影响. 结果表明: 波数较大的情况下, 高次散射波引起柱体上的波浪力不能忽略; 结构间距较大的情况下, 虽然高次波的作用有减小的趋势但仍然明显; 高次散射波来自多个柱体对入射波的散射, 柱体数目的增加后, 高次波的影响会增加, 结构所受的高次波作用因参数变化而起的波动会变剧烈; 高次波对上游柱体波浪力的贡献较对下游柱体的贡献大.      

On diffraction of acoustic and electric waves in piezoelectric medium by an absorbent half-plane electrode

Diffraction of incident acoustic and incident electric waves in a transversally isotropic piezoelectric medium at the boundary of a half-plane absorbent electrode is systematically investigated using the quasi-hyperbolic approximation. The electrode is assumed to be very thin so that its thickness and stiffness can be neglected. By exact inversion, the explicit expressions for the scattering waves are obtained. A closed form solution is obtained by applying Laplace transformations and the Wiener–Hopf technique. By means of the Cagniard–de Hoop method a detailed investigation of the structure of the electro-acoustic wave is conducted. The mode conversion between electric and acoustic waves, the effect of electro-acoustic head wave, the Bleustein–Gulyaev surface wave and the structure of the wave in terms of the type of the incident wave (acoustic or electric) and its angle of incidence are analyzed in detail. It is shown that in piezoelectric materials, absorbent electrodes are neither completely opaque nor completely transparent to electric and acoustic waves. The dynamic field intensity factors at the tip of the electrode are functions of the angle of incidence and time; they are derived explicitly and discussed through a detailed numerical analysis.     

Characteristic analysis for stress wave propagation in transversely isotropic fluid-saturated porous media

According to generalized characteristic theory, a characteristic analysis for stress wave propagation in transversely isotropic fluid-saturated porous media was performed. The characteristic differential equations and compatibility relations along bicharacteristics were deduced and the analytical expressions for wave surfaces were obtained. The characteristic and shapes of the velocity surfaces and wave surfaces in the transversely isotropic fluid-saturated porous media were discussed in detail. The results also show that the characteristic equations for stress waves in pure solids are particular cases of the characteristic equations for fluid-saturated porous media.     

Propagation of electroacoustic waves in the transversely isotropic piezoelectric medium reinforced by randomly distributed cylindrical inhomogeneities

The propagation of electroacoustic waves in a piezoelectric medium containing a statistical ensemble of cylindrical fibers is considered. Both the matrix and the fibers consist of piezoelectric transversely isotropic material with symmetry axis parallel to the fiber axes. Special emphasis is given on the propagation of an electroacoustic axial shear wave polarized parallel to the axis of symmetry propagating in the direction normal to the fiber axis.The scattering problem of one isolated continuous fiber (“one-particle scattering problem”) is considered. By means of a Green’s function approach a system of coupled integral equations for the electroelastic field in the medium containing a single inhomogeneity (fiber) is solved in closed form in the long-wave approximation. The total scattering cross-section of this problem is obtained in closed form and is in accordance with the electroacoustic analogue of the optical theorem.The solution of the one-particle scattering problem is used to solve the homogenization problem for a random set of fibers by means of the self-consistent scheme of effective field method. Closed form expressions for the dynamic characteristics such as total cross-section, effective wave velocity and attenuation factor are obtained in the long-wave approximation.     

On the Volterra's Distortions Theory

This paper has been spurred by the great influence of the Volterra's distortions theory on the modern theory of dislocations and particularly of the ones which study the crystals, aelotropic materials. Volterra developed his theory, as well as the explicit solutions, for elastic, homogeneous, isotropic, and right hollow cylinders. Here we are subduing the assumptions developing the theory for transversely isotropic hollow cylinders.     

具有表面效应的压电半空间中的表面波

纳米科技的快速发展使压电纳米结构在纳米机电系统中得到广泛应用,形成了诸如纳米压电电子学等新的研究方向.与传统的宏观压电材料相比,在纳米尺度下压电材料往往呈现出不同的力学特性,而造成这种差异的原因之一便是表面效应.本文基于Stroh公式、Barnett-Lothe积分矩阵及表面阻抗矩阵,研究计入表面效应的任意各向异性压电半空间中的表面波传播问题,导出了频散方程.针对横观各向同性压电材料,假设矢状平面平行于材料各向同性面,发现Rayleigh表面波和B-G波解耦,并得到各自的显式频散方程.结果表明,Rayleigh表面波的波速小于偏振方向垂直于表面的体波,而B-G波的波速小于偏振方向垂直于矢状平面的体波.以PZT-5H材料为例,用数值方法考察表面残余应力和电学边界条件对表面波频散特性的影响发现:表面残余应力只对第一阶Rayleigh波有明显的影响;电学开路情形的B-G波比电学闭路情形的B-G波传播快.本文工作可为纳米表面声波器件的设计或压电纳米结构的无损检测提供理论依据.     

The effect of rotation and initial stress on the propagation of waves in a transversely isotropic elastic solid

In this paper the equations governing small amplitude motions in a rotating transversely isotropic initially stressed elastic solid are derived, both for compressible and incompressible linearly elastic materials. The equations are first applied to study the effects of initial stress and rotation on the speed of homogeneous plane waves propagating in a configuration with uniform initial stress. The general forms of the constitutive law, stresses and the elasticity tensor are derived within the finite deformation context and then summarized for the considered transversely isotropic material with initial stress in terms of invariants, following which they are specialized for linear elastic response and, for an incompressible material, to the case of plane strain, which involves considerable simplification. The equations for two-dimensional motions in the considered plane are then applied to the study of Rayleigh waves in a rotating half-space with the initial stress parallel to its boundary and the preferred direction of transverse isotropy either parallel to or normal to the boundary within the sagittal plane. The secular equation governing the wave speed is then derived for a general strain–energy function in the plane strain specialization, which involves only two material parameters. The results are illustrated graphically, first by showing how the wave speed depends on the material parameters and the rotation without specifying the constitutive law and, second, for a simple material model to highlight the effects of the rotation and initial stress on the surface wave speed.     

A theory of composites modeled as interpenetrating solid continua

The differential equations and boundary conditions describing the behavior of a finitely deformable, heat-conducting composite material are derived by means of a systematic application of the laws of continuum mechanics to a well-defined macroscopic model consisting of interpenetrating solid continua. Each continuum represents one identifiable constituent of the N-constituent composite. The influence of viscous dissipation is included in the general treatment. Although the motion of the combined composite continuum may be arbitrarily large, the relative displacement of the individual constituents is required to be infinitesimal in order that the composite not rupture. The linear version of the equations in the absence of heat conduction and viscosity is exhibited in detail for the case of the two-constituent composite. The linear equations are written for both the isotropic and transversely isotropic material symmetries. Plane wave solutions in the isotropic case reveal the existence of high-frequency (optical type) branches as well as the ordinary low-frequency (acoustic type) branches, and all waves are dispersive. For the linear isotropic equations both static and dynamic potential representations are obtained, each of which is shown to be complete. The solutions for both the concentrated ordinary body force and relative body force are obtained from the static potential representation.     

Analytical extension of Finite Element solution for computing the nonlinear far field of ultrasonic waves scattered by a closed crack

Nonlinear scattering of ultrasonic waves by closed cracks subject to contact acoustic nonlinearity (CAN) is determined using a 2D Finite Element (FE) coupled with an analytical approach. The FE model, which includes unilateral contact with Coulomb friction to account for contact between crack faces, provides the near-field solution for the interaction between in-plane elastic waves and a crack of different orientations. The numerical solution is then analytically extended in the far-field based on a frequency domain near-to-far field transformation technique, yielding directivity patterns for all linear and nonlinear components of the scattered waves. The proposed method is demonstrated by application to two nonlinear acoustic problems in the case of tone-burst excitations: first, the scattering of higher harmonics resulting from the interaction with a closed crack of various orientations, and second, the scattering of the longitudinal wave resulting from the nonlinear interaction between two shear waves and a closed crack. The analysis of the directivity patterns enables us to identify the characteristics of the nonlinear scattering from a closed crack, which provides essential understanding in order to optimize and apply nonlinear acoustic NDT methods.     

Transmission boundary condition for mixed variable finite-difference method

All kinds of numerical methods based on the discretization of first-order velocity-stress hyperbolic system of elastic wave equations are named mixed variable finite-difference method in this paper, and the transmission boundary condition of artificial boundaries is studied here by the mixed variable finite-difference method. The transmission condition of complex geometrical boundaries of a transversely isotropic medium is presented based on the characteristic variables of wavefields propagating in the normal direction of the boundary. The boundary condition proposed in this paper is a local artificial boundary condition, with which the computation cost is very low. Elastic wave propagations in transversely isotropic medium are modelled by the staggered grid finite-difference method incorporated with the boundary condition presented. Numerical results and analysis of reflection coefficients show that the reflections of incoming waves by artificial boundaries are efficiently reduced.Supported by National Natural Science Foundation and Liaoning Province Science Foundation.     

横观各向同性液体饱和多孔介质中平面波的传播

基于孔隙介质的Ｂｉｏｔ理论１，研究了横观各向同性液体饱和多孔介质中平面波的传播特性。首先导出了波传播的特征方程并给出了其解析解，结果显示：有４种不同波速的平面体波传播；第一准纵波，第二准纵波，准横波和反平面横波。文中给出了波速和衰减的解析表达式，数值计算了频散曲线和衰减曲线，并讨论了各类准体波位移之间的耦合关系。     

An exact transient analysis of plane wave diffraction by a crack in an orthotropic or transversely isotropic solid

A transient plane strain analysis of diffraction of plane waves by a semi-infinite crack in an unbounded orthotropic or transversely isotropic solid is performed. The waves approach the crack at a general oblique angle, and are of two types, a normal stress pulse and a shear stress pulse, i.e. a P- and an SV-wave, respectively, in the isotropic limit. A class of materials that includes this limit and beryl, cobalt, ice, magnesium and titanium is chosen for illustration, and exact solutions are obtained for the initial/mixed boundary value problems.In contrast to related work, a factorization in the Laplace transform space is used to simplify the solution forms and the Wiener-Hopf component of the solution process, and to yield a more compact expression for the Rayleigh wave speed. Calculations for this speed, the two allowable, direction-dependent, plane wave speeds, and quantities related to the Mode I and Mode II dynamic stress intensity factors are given for the five anisotropic materials mentioned.