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.     

The exact solution for the general bending problems of conical shells on the elastic foundation

The general bending problem of conical shells on the elastic foundation (Winkler Medium) is not solved. In this paper, the displacement solution method for this problem is presented. From the governing differential equations in displacement form of conical shell and by introducing a displacement function U(s,θ), the differential equations are changed into an eight-order soluble partial differential equation about the displacement function U(s,θ) in which the coefficients are variable. At the same time, the expressions of the displacement and internal force components of the shell are also given by the displacement function U(s θ). As special cases of this paper, the displacement function introduced by V.S. Vlasov in circular cylindrical shell^[5], the basic equation of the cylindrical shell on the elastic foundation and that of the circular plates on the elastic foundation are directly derived.Under the arbitrary loads and boundary conditions, the general bending problem of the conical shell on the elastic foundation is reduced to find the displacement function U(s,θ).The general solution of the eight-order differential equation is obtained in series form. For the symmetric bending deformation of the conical shell on the elastic foundation, which has been widely usedinpractice,the detailed numerical results and boundary influence coefficients for edge loads have been obtained. These results have important meaning in analysis of conical shell combination construction on the elastic foundation,and provide a valuable judgement for the numerical solution accuracy of some of the same type of the existing problem.     

Rayleigh waves in orthotropic fluid-saturated porous media

In this paper, we are interested in the propagation of Rayleigh waves in orthotropic fluid-saturated porous media. This problem was investigated by Liu and Liu (2004). The authors have derived the secular equation of the wave but that secular equation is still in implicit form. The main aim of this paper is to derive explicit secular equation of the wave. By employing the method of polarization vector, the secular equations of Rayleigh waves in explicit form is obtained. This equation recovers the dispersion equation of Rayleigh waves propagating in pure orthotropic elastic half-spaces. Remarkably, the secular equation obtained is not a complex equation as the one derived by Liu and Liu, it is a really real equation.     

Bilinear form and soliton interactions for the modified Kadomtsev–Petviashvili equation in fluid dynamics and plasma physics

In this paper, we investigate the modified Kadomtsev–Petviashvili (mKP) equation for the nonlinear waves in fluid dynamics and plasma physics. By virtue of the rational transformation and auxiliary function, new bilinear form for the mKP equation is constructed, which is different from those in previous literatures. Based on the bilinear form, one- and two-soliton solutions are obtained with the Hirota method and symbolic computation. Propagation and interactions of shock and solitary waves are investigated analytically and graphically. Parametric conditions for the existence of the shock, elevation solitary, and depression solitary waves are given. From the two-soliton solutions, we find that the (i) parallel elastic interactions can exist between the (a) shock and solitary waves, and (b) two elevation/depression solitary waves; (ii) oblique elastic interactions can exist between the (a) shock and solitary waves, and (b) two solitary waves; (iii) oblique inelastic interactions can exist between the (a) two shock waves, (b) two elevation/depression solitary waves, and (c) shock and solitary waves.     

弹性介质中几何和运动非线性问题的类孤波解

将弹性介质 的几何和运动非线性方程简化成具有电磁场中的Ｂｏｒｎ－Ｉｎｆｅｌｄ方程的形式，并证明了该方程的类孤波解的存在．     

Analysis of scattering waves in an elastic layered medium by means of the complete eigenfunction expansion form of the Green’s function

A scattering problem due to an object and a plane incident wave in an elastic layered half space is presented in this paper. The complete eigenfunction expansion form of the Green’s function developed by the author and the boundary integral equation method are introduced into the analysis. First, the complete eigenfunction expansion form of the Green’s function is investigated for its application to the scattering problem. A comprehensive explanation is also given for the fact that the complex Rayleigh wave modes exhibit standing waves. Next, a method for the analysis of scattering waves by means of the Green’s function is presented. The advantage of the present method is that the formulation itself is independent of the number of layers and the scattering waves can be decomposed into the modes for the spectra defined for the layered medium. Several numerical calculations are performed to examine the efficiency of the present method as well as the properties of the scattering waves. According to the numerical results, the complete eigenfunction expansion form of the Green’s function provides accurate values for application to a boundary element analysis. The spectral structure and radiation patterns of the scattering wave are presented and investigated. The differences in directionality can be found from the radiation patterns of the scattering waves decomposed into the modes for the spectra.     

Anisotropic elastic materials capable of a three-dimensional deformation (static or dynamic) with only one displacement component and uncoupling of all three displacement components

It is shown that there are anisotropic elastic materials that are capable of a non-uniform three-dimensional deformation with only one displacement component. For wave propagation, the equation of motion can be cast in the form of the differential equation for acoustic waves. For elastostatics, the equation of equilibrium reduces to Laplace’s equation. The material can be monoclinic, orthotropic, tetragonal, hexagonal or cubic. There are also anisotropic elastic materials that uncouple all three displacement components. The governing equation for each of the uncoupled displacement can be cast in the form of the differential equation for acoustic waves in the case of dynamic or Laplace’s equation in the case of static. The material can be orthotropic, tetragonal, hexagonal or cubic.     

Investigation of guided waves propagation in orthotropic viscoelastic carbon–epoxy plate by Legendre polynomial method

The effect of viscoelasticity on the guided waves propagation in viscoelastic plate has been investigated according to multi-aspect. To this purpose, an extension of the Legendre polynomial method is proposed to formulate the guided waves equation in orthotropic viscoelastic plate composed of carbon–epoxy. The validity of the proposed Legendre polynomial method is illustrated by comparison with available data. The convergence of the method is discussed through a numerical example. The hysteretic and Kelvin–Voigt viscoelastic models are used to integrate the imaginary part of the complex stiffness matrix associated with the viscoelastic plate in this study. Accordingly, both viscoelastic models do not affect on the dispersion curves results. However, appreciable effects are seen in the attenuation curves. Also, the sensitivity of the guided waves propagation caused by variations of elastic and viscoelastic modulus has been studied in detail. Finally, the advantages of the Legendre polynomial method are described.     

Simulations of waves in elastic solids of cubic symmetry by the Conservation Element and Solution Element method

This paper reports a new numerical approach to simulate waves in elastic solids with cubic symmetry. The governing equation includes the equation of motion and the constitutive relation of the elastic medium. With velocity and stress components as the unknowns, the equations are a set of nine, first-order, hyperbolic partial differential equations. To aid numerical simulation, the characteristic form of the equations is derived. By using the Schur complement in linear algebra, the one-dimensional equations are shown to be equivalent to the Christoffel equations without using the harmonic plane-wave solution. To solve the governing equations by the Conservation Element and Solution Element (CESE) method, we first use Gauss' theorem to recast the equations into a space-time integral form. By integrating the integral equation, space-time flux conservation is imposed over Conservation Elements (CEs). Numerical integration is aided by using prescribed linear discretization of the unknowns in Solution Elements (SEs). A convergence test shows that the CESE method employed is second-order accurate. To demonstrate the capabilities of the present approach, reported numerical results include one-dimensional resonant waves, collinear impact of two blocks, and two-dimensional wave expansion from a point source. Additional results of waves interacting with an interface separating two media with different lattice orientations are also reported. Results compared well with the available analytical solutions. All results show salient features of waves in solids of cubic symmetry.     

横观各向同性体中的埋藏裂纹

本文采用奇异积分方程法分析了横观各向同性体中的埋藏裂纹。建立了张开型埋藏裂纹的Cauchy型奇异积分方程。当裂纹面和弹性对称轴垂直时,得到的裂纹张开位移方程的求解与各向同性情况类似。当裂纹面和弹性对称轴平行时,根据加权余量法,建立了弱解方程。给出两个算例,计算了圆形裂纹和椭圆形裂纹上的张开位移分布。数值结果表明：本文的方法是有效的。横观各向同性体中,埋藏裂纹方位任意时的裂纹张开位移方程,根据本文的方法易于得到。     

A NUMERICAL STUDY OF CYLINDERS IN WAVES AND CURRENTS

The research on combination flow of planar oscillatory flow plus an in-line steady stream is of importance to the situation of structures in waves and current. The combination flow has not been studied extensively. There is still little information about the effect of current and about combined effect of current and waves on hydrodynamic loading of the structures. The present study investigates the combination flow around a circular cylinder using a vortex-based method incorporating vortex moving particles (discrete vortices) with a finite-difference scheme for the vorticity diffusion. The main attention is paid to the effects of a small current on in-line fluid forces and vortex patterns in the wake. Morison's equation and an equation with two drag terms are examined. The results show that the presence of a small current in an oscillatory flow can reduce the drag coefficient significantly. Morison's equation gives reasonably good predictions for the in-line forces for an oscillatory flow plus a small current. The current tends to bring the whole vortex wake downstream and tries to form the stable Karman asymmetrical form in the downstream wake. The present results show certain agreement with some previous experimental results.     

基于新修正偶应力理论的Mindlin层合板热稳定性分析

以新修正偶应力理论为基础,首次提出了机械载荷与热载荷共同作用下的微尺度Mindlin层合板热稳定性模型,该模型只引入一个材料尺度参数,通过虚功原理推导出了控制方程和边界条件,以四边简支方板为例,进行了热稳定性分析,应用纳维叶解法得到解析解。结果表明,所建模型可以捕捉到尺度效应。材料尺度参数值越大,屈曲临界温度越高;当跨厚比增大时,屈曲临界温度下降;随着板几何参数的增大,模型将退化为宏观模型;温度变化量越大,考虑热载荷作用下的屈曲临界载荷越大,尺度效应体现越显著。     

求解二维结构-声耦合问题的一种半数值半解析方法

基于传递矩阵法和虚拟源强模拟技术提出了一种求解在谐激励作用下二维结构-声相互作用问题的半数值半解析法．在足够小的积分步长内，文中对任意形状弹性环沿周向曲线坐标的非齐次状态微分方程组，建立了一种齐次扩容方法．对于外声场，采用多圆形虚拟源强配置方案。并在每一条圆形配置曲线上将源强密度函数用Fourier级数展开，同时结合快速Fourier变换法，提出了一种高精度、高效率求解任意形状二维孔穴Helmholtz外问题的快速算法．在耦合方程的求解方面，根据叠加原理，将外激励和虚拟源强的Fourier级数展开项作为广义力分别作用在弹性环上，借助齐次扩容方法和精细积分法求得弹性环的状态向量，再利用流固交接条件和最小二乘法直接建立了耦合系统的求解方程．文中给出了二个典型弹性环在集中谐激励力作用下声辐射算例，计算结果表明该文方法较通常采用的混合FE-BE法更为有效．     

A SEMI－ANALYTICAL AND SEMI－NUMERICAL METHOD FOR SOLVING 2－D SOUND－STRUCTURE INTERACTION PROBLEMS

Based on the transfer matrix method and the virtual source simulation technique,this paper proposes a novel semi-analytical and semi-numerical method for solving 2-D sound-structure interaction problems under a harmonic excitation. Within any integration segment,as long as its length is small enough, along the circumferential curvilinear coordinate, the nonhomogeneous matrix differential equation of an elastic ring of complex geometrical shape can be rewritten in terms of the homogeneous one by the method of extended homogeneous capacity proposed in this paper. For the exterior fluid domain, the multi-circular virtual source simulation technique is adopted. The source density distributed on each virtual circular curve may be ex-panded as the Fourier‘‘s series. Combining with the inverse fast Fourier transformation, a higher accuracy and efficiency method for solving 2-D exterior Helmholtz‘s problems is presented in this paper. In the aspect of solution to the coupling equations, the state vectors of elastic ring induced by the given harmonic excitation and generalized forces of coefficients of the Fourier series can be obtained respectively by using a high precision integration scheme combined with the method of extended homogeneous capacity put forward in this paper. According to the superposition principle and compatibility conditions at the interface between the elastic ring and fluid, the algebraic equation of system can be directly constructed by using the least square approximation. Examples of acoustic radiation from two typical fluid-loaded elastic rings under a harmonic concentrated force are presented. Numerical results show that the method proposed is more efficient than the mixed FE-BE method in common use.     

A boundary element method for calculating elastic waves scattered by axisymmetric inclusion

A Boundary Element Method (BEM) is described to compute the scattering of elastic waves by an axisymmetric inclusion in an infinite elastic medium. The boundary loads applied to the inclusion is expanded in terms of Fourier series in an infinite space. The boundary integral equation is solved in the general direction of the axisymmetric inclusion by BEM. The problem of the 3-D scattering of elastic waves is reduced to a 1-Done. According to the geometric features of the axisymmetric in clusion the ring shell elements are adopted in this method. A comparison is made with other BEM methods. The numerical results show this method can reduce the amount of calculation and enhance the speed of convergence. Supported by Foundation of Ph. D Program of State Education Commission of China     

移动荷载作用下饱和土地基中的波动特性分析

基于Biot波动方程,经过Fourier变换和逆Fourier变换后可获得波数-频率域以及时间-空间域的解析解。通过数值分析的手段研究了移动荷载作用下饱和多孔弹性地基中波的传播特性。重点就弥散曲线、多谱勒效应、波的成分和动力响应频率等几个特性进行了分析,发现饱和土地基由于比弹性地基多了一项流体介质,波动特性明显差异于弹性介质。     

Numerical solution of a nonlinear reaction-diffusion equation

A nonlinear reaction-diffusion equation is studied numerically by a Petrov-Galerkinfinite element method,which has been proved to be2nd-order accurate in time and4th-orderin space.The comparison between the exact and numerical solutions of progressive wavesshows that this numerical scheme is quite accurate,stable and efficient.It is also shown thatany local disturbance will spread,have a full growth and finally form two progressive wavespropagating in both directions.The shape and the speed of the long term progressive wavesare determined by the system itself,and do not depend on the details of the initial values.     

On two-dimensional linear waves in Blasius boundary layer over viscoelastic layers

This work concerns the direct numerical simulation of small-amplitude two-dimensional ribbon-excited waves in Blasius boundary layer over viscoelastic compliant layers of finite length. A vorticity-streamfunction formulation is used, which assures divergence-free solutions for the evolving flow fields. Waves in the compliant panels are governed by the viscoelastic Navier's equations. The study shows that Tollmien–Schlichting (TS) waves and compliance-induced flow instability (CIFI) waves that are predicted by linear stability theory frequently coexist on viscoelastic layers of finite length. In general, the behaviour of the waves is consistent with the predictions of linear stability theory. The edges of the compliant panels, where abrupt changes in wall property occur, are an important source of waves when they are subjected to periodic excitation by the flow. The numerical results indicate that the non-parallel effect of boundary-layer growth is destabilizing on the TS instability. It is further demonstrated that viscoelastic layers with suitable properties are able to reduce the amplification of the TS waves, and that high levels of material damping are effective in controlling the propagating CIFI.     

A novel method of solving the exact equations for capillary-gravity waves

A new method is presented for the computation of two-dimensional periodicprogressive surface waves propagating under the combined influence of gravity and surfacetension.The nonlinear surface is expressed by Fourier series with finite number of terms,after the computational domain is transformed into a unit circle.The dynamic boundaryequation is used in its exact nonlinear form and the coefficients of Fourier series are foundby the Nweton-Raphson method successively.This is a neat method,Yielding highprescision with little computational effort.     

正交各向异性平板开孔弹性波的衍射与动应力集中

采用各向异性平板弯曲波动方程及摄动方法，对正交各向异性无限平板开孔弹性波的衍射及动应力集中的问题进行了分析研究，得到了在稳态波作用下此种忆边界条件波动问题的渐近形式分析解。同时采用复变函数方法及保角映射技术，为求解正交各向异性无限平板开孔弹性波的衍射及动力集中问题提供了一种统一规范的分析方法。