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.     

On the Nonlinear Continuum Theory of Dislocations: A Gauge Field Theoretical Approach

A nonlinear continuum theory of material bodies with continuously distributed dislocations is presented, based on a gauge theoretical approach. Firstly, we derive the canonical conservation laws that correspond to the group of translations and rotations in the material space using Noether’s theorem. These equations give us the canonical Eshelby stress tensor as well as the total canonical angular momentum tensor. The canonical Eshelby stress tensor is neither symmetric nor gauge-invariant. Based on the Belinfante-Rosenfeld procedure, we obtain the gauge-invariant Eshelby stress tensor which can be symmetric relative to the reference configuration only for isotropic materials. The gauge-invariant angular momentum tensor is obtained as well. The decomposition of the gauge-invariant Eshelby stress tensor in an elastic and in a dislocation part gives rise to the derivation of the famous Peach-Koehler force.     

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.     

Tensile stability of cylindrical membranes

The tensile stability of rotationally symmetric thin membranes composed of isotropic, incompressible and elastic materials is considered by investigating under what conditions the initial equilibrium configuration can bifurcate to another rotationally symmetric equilibrium mode.The general equilibrium equations of a rotationally symmetric membrane are first derived in cylindrical coordinates. The initially cylindrical membrane is studied. The classic solution, which assumes homogeneous deformations, is shown to be a special case of the general equations. Perturbation theory is employed to find the bifurcation points from the homogeneous mode.This study shows that, for the chosen boundary conditions, no rotationally symmetric equilibrium mode exists near the cylindrical mode except the cylindrical mode itself. This corresponds to all experimental data that the author is aware of. The initially cylindrical membrane either remains cylindrical or goes into a non-rotationally symmetric mode.     

An asymmetric theory of nonlocal elasticity—Part 1. Quasicontinuum theory

In this article, an asymmetric theory of nonlocal elasticity is developed on the basis of three dimensional atomic lattice model, the Galileo invariance for constitutive equations and by use of Fourier transformation of generalized function and energy method. It is shown that nonlocal characteristic functions (or constitutive parameters of internal elastic energy) can be explicitly expressed in terms of interacting forces connecting atoms, and the general model of nonlocal theory with rotation effects is asymmetric. Both symmetric stress and anti-symmetric stress is a nonlocal function of strain and local rotation for anisotropic materials. For isotropic materials, symmetric stress is only a nonlocal function of strain, while antisymmetric stress is only a nonlocal function of local rotation.     

Transformation elastodynamics and cloaking for flexural waves

The paper addresses an important issue of cloaking transformations for fourth-order partial differential equations representing flexural waves in thin elastic plates. It is shown that, in contrast with the Helmholtz equation, the general form of the partial differential equation is not invariant with respect to the cloaking transformation. The significant result of this paper is the analysis of the transformed equation and its interpretation in the framework of the linear theory of pre-stressed plates. The paper provides a formal framework for transformation elastodynamics as applied to elastic plates. Furthermore, an algorithm is proposed for designing a broadband square cloak for flexural waves, which employs a regularised push-out transformation. Illustrative numerical examples show high accuracy and efficiency of the proposed cloaking algorithm. In particular, a physical configuration involving a perturbation of an interference pattern generated by two coherent sources is presented. It is demonstrated that the perturbation produced by a cloaked defect is negligibly small even for such a delicate interference pattern.     

Rayleigh lamb waves in micropolar isotropic elastic plate

The propagation of waves in a homogeneous isotropic micropolar elastic cylindrical plate subjected to stress free conditions is investigated. The secular equations for symmetric and skew symmetric wave mode propagation are derived. At short wave limit, the secular equations for symmetric and skew symmetric waves in a stress free circular plate reduces to Rayleigh surface wave frequency equation. Thin plate results are also obtained. The amplitudes of displacements and microrotation components are obtained and depicted graphically. Some special cases are also deduced from the present investigations. The secular equations for symmetric and skew symmetric modes are also presented graphically.     

复合材料切口应力奇性指数计算

提出一种计算广义平面应交状态下复合材料切口应力奇性指数的新方法.在切口尖端的位移幂级数渐近展开式被引入正交各向异性材料的物理方程后,将用位移表示的应力分量代入切口端部柱状邻域的线弹性理论控制方程,切口应力奇性指数的计算被转化为常微分方程组特征值的求解.采用插值矩阵法求解该常微分方程组,可一次性地获取切口尖端多阶应力奇性指数.本法适合平面和反平面应力场耦合或解耦的情形,并可退化计算裂纹或各向同性材料切口的应力奇性指数.算例表明,所提方法对分析复合材料切口应力奇性指数是一种准确有效的手段.     

Finite deformation of elastic membranes with application to the stability of an inflated and extended tube

A theory is formulated for the finite deformation of a thin membrane composed of homogeneous elastic material which is isotropic in its undeformed state. The theory is then extended to the case of a small deformation superposed on a known finite deformation of the membrane. As an example, small deformations of a circular cylindrical tube which has been subjected to a finite homogeneous extension and inflation are considered and the equations governing these small deformations are obtained for an incompressible material. By means of a static analysis the stability of cylindrically symmetric modes for the inflated and extended cylinder with fixed ends is determined and the results are verified by a dynamic analysis. The stability is considered in detail for a Mooney material. Methods are developed to obtain the natural frequencies for axially symmetric free vibrations of the extended and inflated cylindrical membrane. Some of the lower natural frequencies are calculated for a Mooney material and the methods are compared.     

On some connections between equivalent single material and mixture theory models for fiber reinforced hyperelastic materials

There are two approaches that can be used to model the large strain mechanical response of material systems in which elastic fibers are embedded in an elastic matrix. In the first approach, a fiber reinforced material undergoing large deformation is homogenized in the sense that it is assumed to act as an equivalent single material that is transversely isotropic and hyperelastic. Both constituents then share a common reference configuration, which is typically assumed to be a natural or stress-free configuration for the equivalent single material. The stress depends on a single deformation gradient defined with respect to the natural configuration.In the second approach, the fiber/matrix system is treated as a mixture, with the matrix and the fibrous constituents having their own reference configurations and material symmetries. The total stress depends on the deformation gradients and material symmetries for both constituents, defined with respect to their reference configurations.Under appropriate assumptions, the constitutive theory developed using mixture theory can coincide with the constitutive theory assuming an equivalent single material that is transversely isotropic and hyperelastic. This paper explores the connection between the two approaches by considering the various reference configurations and material symmetries.     

轴对称横观各向同性层状弹性半空间问题受力分析

本文从柱坐标系弹性力学基本方程出发,将位移场和应力场在径向进行Ｈａｎｋｅｌ变换,利用常微分主程求解原理,直接得出在轴对称荷载作用下横观各向同性半无限弹性空间的位移场,利用此结果推导出单层元的刚度矩阵。     

Exact Solutions for a Thick Elastic Plate with a Thin Elastic Surface Layer

A procedure has been developed in previous papers for constructing exact solutions of the equations of linear elasticity in a plate (not necessarily thin) of inhomogeneous isotropic linearly elastic material in which the elastic moduli depend in any specified manner on a coordinate normal to the plane of the plate. The essential idea is that any solution of the classical equations for a hypothetical thin plate or laminate (which are two-dimensional theories) generates, by straightforward substitutions, a solution of the three-dimensional elasticity equations for the inhomogeneous material. In this paper we consider a thick plate of isotropic elastic material with a thin surface layer of different isotropic elastic material. It is shown that the interface tractions and in-plane stress discontinuities are determined only by the initial two-dimensional solution, without recourse to the three-dimensional elasticity theory. Two illustrative examples are described.     

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.     

黏弹性体界面裂纹的冲击响应

研究两半无限大黏弹性体界面Griffith裂纹在反平面剪切突出载荷下,裂纹尖端动应力强度因子的时间响应,首先,运用积分变换方法将黏弹性混合黑社会问题化成变换域上的对偶积分方程,通过引入裂纹位错密度函数进一步化成Cauchy型奇异积分方程,运用分片连续函数法数值求解奇异积分方程,得到变换域内的动应力强度因子,再用Laplace积分变换数值反演方法,将变换域的解反演到时间域内,最终求得动应力强度因子的时间响应,并对黏弹性参数的影响进行分析。     

Une transformation du problème d'élasticité linéaire en vue d'application au problème de l'inclusion et aux fonctions de Green

A simple transformation of the problem of the linear elastic structure is presented. The transformed problem corresponds to a new problem of linear elastic structure with different behaviour, geometry and prescribed forces and displacements. The transformed problem can be easier to study, or can correspond to cases with well-known solutions. By means of this transformation, the problem of ellipsoidal inclusion is transformed into a problem of spherical inclusion, the analytical results known for the Eshelby tensor for an isotropic or transversely isotropic matrix are extended to more general cases of matrix behaviour, and finally, close form expressions of the Green function for an infinite medium are derived for some cases of elastic behaviour without transversal isotropy or orthotropy.     

The dynamic response of an incompressible non-linearly elastic membrane tube subjected to a dynamic extension

The dynamic response of an isotropic hyperelastic membrane tube, subjected to a dynamic extension at its one end, is studied. In the first part of the paper, an asymptotic expansion technique is used to derive a non-linear membrane theory for finite axially symmetric dynamic deformations of incompressible non-linearly elastic circular cylindrical tubes by starting from the three-dimensional elasticity theory. The equations governing dynamic axially symmetric deformations of the membrane tube are obtained for an arbitrary form of the strain-energy function. In the second part of the paper, finite amplitude wave propagation in an incompressible hyperelastic membrane tube is considered when one end is fixed and the other is subjected to a suddenly applied dynamic extension. A Godunov-type finite volume method is used to solve numerically the corresponding problem. Numerical results are given for the Mooney-Rivlin incompressible material. The question how the present numerical results are related to those obtained in the literature is discussed.     

Axisymmetric Deformation of a Transversely Isotropic Cylindrical Body: A Hamiltonian State-Space Approach

A state-space approach for exact analysis of axisymmetric deformation and stress distribution in a circular cylindrical body of transversely isotropic material is developed. By means of Hamiltonian variational formulation via Legendre’s transformation, the basic equations in cylindrical coordinates are formulated into a state-space framework in which the unknown state vector comprises the displacements and associated stress components as the dual variables and the system matrix possesses the symplectic characteristics of a Hamiltonian system. Upon delineating the symplecticity of the formulation, a viable solution approach using eigenfunction expansion is developed. For illustration, an exact analysis of a finite thick-walled circular cylinder under internal and external pressures is presented, with emphasis on the end effects.     

Residual-stress determination of concentric layers of cylindrically orthotropic materials

Equations have been obtained for determining residual stresses in the wall of a hollow, axially symmetric body consisting of concentric layers of elastically dissimilar materials, all having cylindrical elastic orthotropy. These equations permit residual normal stresses in the radial, circumferential, and axial directions and residual shear stresses on planes normal to the axis of the body to be calculated from measurements of the strains developed on the inner or outer cylindrical surface of the body as thin layers of stressed material are serially removed from the outer or inner surfaces, respectively. The equations are applied to a parametric study of stresses in an elastically isotropic, two-component body to determine the nature of the differences in stresses between the composite body and a homogeneous body as a function of the difference in elastic constants.     

考虑温度效应的弹性损伤一般理论

现有的各种损伤理论基本上都是关于等温问题的 ,且在不同程度上依赖于某些经验假设。本文在严格的不可逆热力学理论基础之上 ,建立了考虑温度效应的弹性损伤一般理论。推导出热弹性各向同性与各向异性损伤材料全部本构方程的一般形式 ,其中包括应力 应变关系、熵密度方程、损伤对偶张量表达式、热 固 损伤耦合的热传导方程和损伤演化方程。它们的特殊形式包含了等温各向同性与各向异性弹性损伤的本构方程