Transformational cloaking from seismic surface waves by micropolar metamaterials with finite couple stiffness

Transformational elastodynamics can be used to protect sensitive structures from harmful waves and vibrations. By designing the material properties in a region around the sensitive structure, a cloak, the incident waves can be redirected as to cause minimal or no harmful response on the pertinent structure. In this paper, we consider such transformational cloaking built up by a suitably designed metamaterial exhibiting micropolar properties. First, a theoretically perfect cloak is obtained by designing the properties of an (unphysical) restricted micropolar material within the surrounding medium. Secondly, we investigate the performance of the cloak under more feasible design criteria, relating to finite elastic parameters. In particular, the behavior of a physically realizable cloak built up by unrestricted micropolar elastic media is investigated. Numerical studies are conducted for the case of buried as well as surface breaking structures in 2D subjected to incident Rayleigh waves pertinent to seismic loading. The studies show how the developed cloaking procedure can be utilized to substantially reduce the response of the structure. In particular, the results indicate the performance of the cloak in relation to constraints on the elastic parameters.     

Antiplane elastic wave cloaking using metamaterials,homogenization and hyperelasticity

We consider the problem of how to cloak objects from antiplane elastic waves using two alternative techniques. The first is the use of a layered metamaterial in the spirit of the work of Torrent and Sanchez-Dehesa (2008) who considered acoustic cloaks, motivated by homogenization theories, whilst the second is the use of a hyperelastic cloak in the spirit of the work of Parnell et al. (2012). We extend the hyperelastic cloaking theory to the case of a Mooney–Rivlin material since this is often considered to be a more realistic constitutive model of rubber-like media than the neo-Hookean case studied by Parnell et al. (2012), certainly at the deformations required to produce a significant cloaking effect. Although not perfect, the Mooney–Rivlin material appears to be a reasonable hyperelastic cloak. This is clearly encouraging for applications. We quantify the effectiveness of the various cloaks considered by plotting the scattering cross section as a function of frequency, noting that this would be zero for a perfect cloak.     

Plasmonic cloaking and scattering cancelation for electromagnetic and acoustic waves

One popular approach to cloaking objects from electromagnetic waves at moderately long wavelengths is the scattering cancelation technique. This mechanism is based on the use of a single homogeneous thin layer to cover an object of interest in order to provide scattering suppression in a given frequency band. This approach has also been recently extended to acoustic waves. This paper provides an investigation of the physical nature of scattering cancelation by a uniform thin layer for both electromagnetic and acoustic waves in inviscid fluids. Two distinct scattering cancelation regions are obtained within the available parameter space: a non-resonant plasmonic cloaking region and an anti-resonant cloaking region, which are identified and compared in both the electromagnetic and acoustic domains. Although both types of operations allow for the suppression of the dominant scattering orders, the resulting internal fields and physical functionality of the cloaks present distinct differences between the two domains. We discuss analogies and differences between these functionalities and their implications in electromagnetic and acoustic cloaking problems, with an insight into their practical implementation.     

Analytical study of spherical cloak/anti-cloak interactions

The intriguing concept of “anti-cloaking” has been recently introduced within the framework of transformation optics (TO), first as a “countermeasure” to invisibility-cloaking (i.e., to restore the scattering response of a cloaked target), and more recently in connection with “sensor invisibility” (i.e., to strongly reduce the scattering response while maintaining the field-sensing capabilities). In this paper, we extend our previous studies, which were limited to a two-dimensional cylindrical scenario, to the three-dimensional spherical case. More specifically, via a generalized (coordinate-mapped) Mie-series approach, we derive a general analytical full-wave solution pertaining to plane-wave-excited configurations featuring a spherical object surrounded by a TO-based invisibility cloak coupled via a vacuum layer to an anti-cloak, and explore the various interactions of interest. With a number of selected examples, we illustrate the cloaking and field-restoring capabilities of various configurations, highlighting similarities and differences with respect to the cylindrical case, with special emphasis on sensor-cloaking scenarios and ideas for approximate implementations that require the use of double-positive media only.     

Realizing the Willis equations with pre-stresses

This paper proves that the linear elastic behavior of the material with inhomogeneous pre-stresses can be described by the Willis equations. In this case, the additional terms in the Willis equations, compared with the classical linear elastic equations for homogeneous media, are related to the gradient of pre-stresses. In this way, the material length scale is naturally incorporated in the framework of continuum mechanics. All these findings also coincide with the results of transformation elastodynamics, so that they can meet the requirement of the principle of material objectivity and the principle of general invariance.     

GURTIN-TYPE REGION-WISE VARIATIONAL PRINCIPLES FOR THERMOPIEZOELECTRIC ELASTODYNAMICS

IntroductionThemethodofconvolution[1- 3]canbeusedinthestudyofpiezoelectricmaterial.Piezoelectricmaterialisoneoftheimportantsmartmaterials,usedinnovelsmartstructure .ThedynamictheoryforthermopiezoelectricelastomediawassearchedbyNowacki[4 ,5 ],Chandrasekharaiah[6 ]andIesan[7].However,region_wisevariationalprinciplesinlinearthermopiezoelectricelastodynamicshavebeennotestablishedsystematically .AlthoughLuoenandKuangJun_shang[8]havedevelopedsimplifiedGurtin_typevariationalprinciplesinlinearther…     

A theoretical solution for axially symmetric problems in elastodynamics

This paper presents a theoretical solution for the basic equation of axisymmetric problems in elastodynamics. The solution is composed of a quasi-static solution which satisfies inhomogeneous boundary conditions and a dynamic solution which satisfies homogeneous boundary conditions. After the quasi-static solution has been obtained an inhomogeneous equation for dynamic solution is found from the basic equation. By making use of eigenvalue problem of a corresponding homogeneous equation, a finite Hankel transform is defined. A dynamic solution satisfying homogeneous boundary conditions is obtained by means of the finite Hankel transform and Laplace transform. Thus, an exact solution is obtained. Through an example of hollow cylinders under dynamic load, it is seen that the method, and the process of computing are simple, effective and accurate.     

Acoustic cloaking: Geometric transform,homogenization and a genetic algorithm

A general process is proposed to experimentally design anisotropic inhomogeneous metamaterials obtained through a change of coordinates in the Helmholtz equation. The method is applied to the case of a cylindrical transformation that allows cloaking to be performed. To approximate such complex metamaterials we apply results of the theory of homogenization and combine them with a genetic algorithm. To illustrate the power of our approach, we design three types of cloaks composed of isotropic concentric layers structured with three types of perforations: curved rectangles, split rings and crosses. These cloaks have parameters compatible with existing technology and they mimic the behavior of the transformed material. Numerical simulations have been performed to qualitatively and quantitatively study the cloaking efficiency of these metamaterials.     

Numerical analysis of three-dimensional acoustic cloaks and carpets

We start by a review of the chronology of mathematical results on the Dirichlet-to-Neumann map which paved the way toward the physics of transformational acoustics. We then rederive the expression for the (anisotropic) density and bulk modulus appearing in the pressure wave equation written in the transformed coordinates. A spherical acoustic cloak consisting of an alternation of homogeneous isotropic concentric layers is further proposed based on the effective medium theory. This cloak is characterized by a low reflection and good efficiency over a large bandwidth for both near and far fields, which approximates the ideal cloak with an inhomogeneous and anisotropic distribution of material parameters. The latter suffers from singular material parameters on its inner surface. This singularity depends upon the sharpness of corners, if the cloak has an irregular boundary, e.g. a polyhedron cloak becomes more and more singular when the number of vertices increases if it is star shaped. We thus analyze the acoustic response of a non-singular spherical cloak designed by blowing up a small ball instead of a point, as proposed in [Kohn, Shen, Vogelius, Weinstein, Inverse Problems 24, 015016, 2008]. The multilayered approximation of this cloak requires less extreme densities (especially for the lowest bound). Finally, we investigate another type of non-singular cloaks, known as invisibility carpets [Li and Pendry, Phys. Rev. Lett. 101, 203901, 2008], which mimic the reflection by a flat ground.     

BEM for transient 2D elastodynamics using multiquadric functions

A new boundary element model for transient dynamic analysis of 2D structures is presented. The dual reciprocity method (DRM) is reformulated for the 2D elastodynamics by using the multiquadric radial basis functions (MQ). The required kernels for displacement and traction particular solutions are derived. Some terms of these kernels are found to be singular; therefore, a new smoothing technique is proposed to solve this problem. Hence, the limiting values of relevant kernels are computed. The validity and strength of the proposed formulation are demonstrated throughout several numerical applications. It is proven from the results that the present formulation is more stable than the traditional DRM, which uses the conical (1 + R) function, especially in predicting results in the far time zone.     

some general theorems and generalized and piecewise generalized variational principles for linear elastodynamics

From the concept of four-dimensional space and under the four kinds of time limit conditions, some general theorems for elastodynamics are developed, such as the principle of possible work action, the virtual displacement principle, the virtual stress-momentum principle, the reciprocal theorems and the related theorems of time terminal conditions derived from it. The variational principles of potential energy action and complementary energy action, the H-W principles, the H-R principles and the constitutive variational principles for elastodynamics are obtained. Hamilton's principle, Toupin's work and the formulations of Ref. [5], [17]-[24] may be regarded as some special cases of the general principles given in the paper. By considering three cases: piecewise space-time domain, piecewise space domain, piecewise time domain, the piecewise variational principles including the potential, the complementary and the mixed energy action fashions are given. Finally, the general formulation of piecewise variati     

一种处理弹性动力边界元法中奇异积分的方法

利用边界元法求解瞬态弹性动力学问题时,时域基本解函数的分段连续性和奇异性为该问题的求解带来很大的困难。为了解决时域基本解中的奇异性问题,本文依据柯西主值的定义,对经过时间解析积分之后的时域基本解进行奇异值分解,将其分成奇异和正则积分两部分;其中正则部分可通过采用常规高斯积分方法来计算,而奇异部分具有简单的形式,可以利用解析积分计算。经过上述操作之后,就可以达到直接消除时域基本解中奇异积分的目的。和传统方法相比,本文方法并不依赖静力学基本解来消除奇异性,是一种直接求解方法。最后给定两个数值算例来验证本文提出方法的正确性和可行性,结果表明使用本文算法可以解决弹性动力学边界积分方程中的奇异性问题。     

Recent advances on the fast multipole accelerated boundary element method for 3D time-harmonic elastodynamics

This article is mainly devoted to a review on fast BEMs for elastodynamics, with particular attention on time-harmonic fast multipole methods (FMMs). It also includes original results that complete a very recent study on the FMM for elastodynamic problems in semi-infinite media. The main concepts underlying fast elastodynamic BEMs and the kernel-dependent elastodynamic FM-BEM based on the diagonal-form kernel decomposition are reviewed. An elastodynamic FM-BEM based on the half-space Green’s tensor suitable for semi-infinite media, and in particular on the fast evaluation of the corresponding governing double-layer integral operator involved in the BIE formulation of wave scattering by underground cavities, is then presented. Results on numerical tests for the multipole evaluation of the half-space traction Green’s tensor and the FMM treatment of a sample 3D problem involving wave scattering by an underground cavity demonstrate the accuracy of the proposed approach. The article concludes with a discussion of several topics open to further investigation, with relevant published work surveyed in the process.     

Unconventional Hamilton-type variational principles for nonlinear elastodynamics of orthogonal cable-net structures

According to the basic idea of classical yin-yang complementarity and modem dual-complementarity,in a simple and unified new way proposed by Luo,the unconven- tional Hamilton-type variational principles for geometrically nonlinear elastodynamics of orthogonal cable-net structures are established systematically,which can fully charac- terize the initial-boundary-value problem of this kind of dynamics.An important in- tegral relation is made,which can be considered as the generalized principle of virtual work for geometrically nonlinear dynamics of orthogonal cable-net structures in mechan- ics.Based on such relationship,it is possible not only to obtain the principle of virtual work for geometrically nonlinear dynamics of orthogonal cable-net structures,but also to derive systematically the complementary functionais for five-field,four-field,three-field and two-field unconventional Hamilton-type variational principles,and the functional for the unconventional Hamilton-type variational principle in phase space and the poten- tial energy functional for one-field unconventional Hamilton-type variational principle for geometrically nonlinear elastodynamics of orthogonal cable-net structures by the general- ized Legendre transformation given in this paper.Furthermore,the intrinsic relationship among various principles can be explained clearly with this approach.     

Cloaking of a circular cylindrical elastic inclusion from antiplane elastic waves and resonance effects

Cloaking of a circular cylindrical elastic inclusion embedded in a homogeneous linear isotropic elastic medium from antiplane elastic waves is studied. The transformation or change-of-variables method is used to determine the material properties of the cloak and the homogenization theory of composites is used to construct a multilayered cloak consisting of many bi-material cells. The large system of algebraic equations associated with this problem is solved by using the concept of multiple scattering with wave expansion coefficient matrices. Numerical results for cloaking of an elastic inclusion and a rigid inclusion are compared with the case of a cavity. It is found that while the cloaking patterns for the three cases are similar, the major difference is that standing waves are generated in the elastic inclusion and the multilayered cloak cannot prevent the motion inside the elastic inclusion, even though the cloak seems nearly perfect. Waves can penetrate into and cause vibrations inside the elastic inclusion, where the amplitude of standing waves depend on the material properties of the inclusion but are very much reduced when compared to the case when there is no cloak. For a prescribed mass density, the displacements inside the elastic cylinder decrease as the shear modulus increases. Moreover, the cloaking of the elastic inclusion over a range of wavenumbers is also investigated. There is significant low frequency scattering even if the cloak consists of a large number of layers. When the wavenumber increases, the multilayered cloak is not effective if the cloak consists of an insufficient number of layers. Resonance effects that occur in cloaking of elastic inclusions are also discussed.     

Composite ceramic-metal plates tested with flexural waves and holography

A method to determine weak or defective bonding areas within the brazed adjoining contact surfaces of composite ceramic-metal plates is proposed. The plates include voids, as well as hairline cracks in the ceramic layer, caused by the high-temperature brazing process. To detect these flaws, transient flexural waves are generated and transmitted through the plate by means of an attached piezoelectric transducer. These pulses are a narrow band signal generated in the time domain. The characteristic wavelengths corresponding to the narrow frequency spectrum may be larger than the size of the flaw. Dynamic holographic interferometry using a double-pulsed twin-cavity Nd:YAG laser, comprising two independent reference beams, was used to measure the response. The experimental interferograms (phase maps) coincide rather precisely with analytic results derived from Mindlin's plate equations (including effects of shear and rotary inertia).     

Monte Carlo analysis of dissociation and recombination behind strong shock waves in nitrogen

Computations are presented for the relaxation zone behind strong, one-dimensional shock waves in nitrogen. The analysis is performed with the direct simulation Monte Carlo method (DSMC). The DSMC code is vectorized for efficient use on a supercomputer. The code simulates translational, rotational and vibrational energy exchange and dissociative and recombinative chemical reactions. A new model is proposed for the treatment of three body recombination collisions in the DSMC technique which usually simulates binary collision events. The new model represents improvement over previous models in that it can be employed with a large range of chemical rate data, does not introduce into the flow field troublesome pairs of atoms which may recombine upon further collision (pseudo-particles) and is compatible with the vectorized code. The computational results are compared with existing experimental data. It is shown that the derivation of chemical rate coefficients must account for the degree of vibrational nonequilibrium in the flow. A nonequilibrium chemistry model is employed together with equilibrium rate data to compute successfully the flow in several different nitrogen shock waves.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1990.     

Nonlinear flexural waves and chaos behavior in finite-deflection Timoshenko beam

Based on the Timoshenko beam theory, the finite-deflection and the axial inertia are taken into account, and the nonlinear partial differential equations for flexural waves in a beam are derived. Using the traveling wave method and integration skills, the nonlinear partial differential equations can be converted into an ordinary differential equation. The qualitative analysis indicates that the corresponding dynamic system has a heteroclinic orbit under a certain condition. An exact periodic solution of the nonlinear wave equation is obtained using the Jacobi elliptic function expansion. When the modulus of the Jacobi elliptic function tends to one in the degenerate case, a shock wave solution is given. The small perturbations are further introduced, arising from the damping and the external load to an original Hamilton system, and the threshold condition of the existence of the transverse heteroclinic point is obtained using Melnikov＇s method. It is shown that the perturbed system has a chaotic property under the Smale horseshoe transform.