On formulas for the velocity of Stoneley waves propagating along the loosely bonded interface of two elastic half-spaces

Formulas for the velocity of Stoneley waves propagating along the loosely bonded interface of two isotropic elastic half-spaces are derived using the complex function method. The derivation also shows that if a Stoneley wave exists, then it is unique. By using the obtained formulas, we can easily reproduce the numerical results previously obtained by Murty [G. S. Murty, Phys. Earth Planet. Interiors 11 (1975), 65–79.] by directly solving the secular equation.     

Identification of the plastic behavior of aluminum plates under free air explosions using inverse methods and full-field measurements

This article describes an inverse method for the identification of the plastic behavior of aluminum plates subjected to sudden blast loads. The method uses full-field optical measurements taken during the first milliseconds of a free air explosion and the finite element method for the numerical prediction of the blast response. The identification is based on a damped least-squares solution according to the Levenberg–Marquardt formulation. Three different rate-dependent plasticity models are examined. First, a combined model based on linear strain hardening and the strain rate term of the Cowper–Symonds model, secondly, the Johnson–Cook model and finally, a combined model based on a bi-exponential relation for the strain hardening term and the strain rate term of the Cowper–Symonds model. A validation of the method and its sensitivity to measurement uncertainties is first provided according to virtual measurements generated with the finite element method. Next, the plastic behavior of aluminum is identified using measurements from real free air explosions obtained from a controlled detonation of C4. The results show that inverse methods can be successfully applied for the identification of the plastic behavior of metals subjected to blast waves. In addition, the material parameters identified with inverse methods enable the numerical prediction of the material’s response with increased accuracy.     

An explicit formulation for the evolution of nonlinear surface waves interacting with a submerged body

An explicit formulation to study nonlinear waves interacting with a submerged body in an ideal fluid of infinite depth is presented. The formulation allows one to decompose the nonlinear wave–body interaction problem into body and free‐surface problems. After the decomposition, the body problem satisfies a modified body boundary condition in an unbounded fluid domain, while the free‐surface problem satisfies modified nonlinear free‐surface boundary conditions. It is then shown that the nonlinear free‐surface problem can be further reduced to a closed system of two nonlinear evolution equations expanded in infinite series for the free‐surface elevation and the velocity potential at the free surface. For numerical experiments, the body problem is solved using a distribution of singularities along the body surface and the system of evolution equations, truncated at third order in wave steepness, is then solved using a pseudo‐spectral method based on the fast Fourier transform. A circular cylinder translating steadily near the free surface is considered and it is found that our numerical solutions show excellent agreement with the fully nonlinear solution using a boundary integral method. We further validate our solutions for a submerged circular cylinder oscillating vertically or fixed under incoming nonlinear waves with other analytical and numerical results. Copyright © 2007 John Wiley & Sons, Ltd.     

Three-dimensional overturned traveling water waves

Traveling gravity-capillary water waves on the interface of a three-dimensional fluid of infinite depth are computed. The vortex sheet formulation with the small scale approximation is used as the mathematical model for the fluid motion. The fluid interface is parameterized isothermally. The traveling wave ansatz for parameterized surfaces is described. Waves are computed using Fourier collocation and quasi-Newton iteration; large amplitude overturned traveling waves are computed via a dimension-breaking based numerical continuation method.     

有限深水中骑行波的显式Hamilton描述

小尺度波（扰动波）迭加在大尺度波（未受扰动波）上形成的波动一般之为“骑行波”。研究了有限可变深度的理想不可压缩流体中的骑行波的显式Hamliltn表示，考虑了自由面上流体与空气之间的表面张力。采用自由面高度和自由面上速度势构成的Hamilton正则变量表示骑行波的动能密度，并在未受扰动波的自由面上作一阶展开。运用复变函数论方法处理了二维流动。先用保角变换将物理平面上的流动区域变换到复势平面上的无限长带形区域，然后在复势平面上用Fourier变换解出Laplace方程，最后经Fourier逆变换求出了扰动波速度热所满足的积分方程。作为特例考虑了平坦底部的流动，导出了动能密度的显式表达式。这里给出的积分方程可以替代相当难解的Hamilton正则方程。通过求解积分方程可得出agrange密度的显式表达式。本文提出的方法约研究骑行波的Hamilton描述以及波的相互作用问题提供了新的途径，有助于了解海面的小尺度波的精细结构。     

Analysis of the wave solution of the elastokinetic equations of a Cosserat continuum for the case of bulk plane waves

A study is made of waves in a Cosserat continuum, whose strain state is characterized by independent displacement and rotation vectors. The propagation of longitudinal and transverse bulk waves is considered. Wave solutions are sought in the form of wave trains specified by a Fourier spectrum of arbitrary shape. It is shown that if the solution is sought in the form of three components of the displacement vector and three components of the rotation vector which depend on time and the longitudinal coordinate, the initial system is split into two systems, one of which describes longitudinal waves, and the other transverse waves. For waves of both types, dispersion relations and analytical solutions in displacement are obtained. The dispersion characteristics of the solutions obtained differ from the dispersion characteristics of the corresponding classical elastic solutions. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 2, pp. 196–203, March–April, 2008.     

Analysis of damped guided waves using the method of multiple scales

We analytically investigate the influence of damping on Lamb waves, which are a specific type of guided wave in two-dimensional plates. Considering material attenuation, we suppose that Lamé constants are complex numbers. This leads to the associated wavenumbers being complex, with the imaginary part of the wavenumber being associated with effect of attenuation of the guided wave. In this paper, we show how dispersion curves and attenuation coefficients can be obtained using the self-adjointness and the method of multiple scales (MMS), which is a type of perturbation method. Using the self-adjointness and the MMS, we can calculate the frequency- and wavenumber-dependent attenuation coefficients from the integral values and boundary values of a corresponding eigenfunction with respect to each propagation mode. This analytical method can yield not only dispersion curves but also mode-by-mode attenuation coefficients regardless of the numerical initial values, unlike numerical approaches using the Newton method. Thus, the proposed method can more easily calculate the attenuation coefficients with respect to a particular mode than conventional methods. Furthermore, the results obtained by proposed method were in good agreement with those obtained by the semi-analytical finite element (SAFE) method, which validates the proposed method.     

Particle trajectories under interactions between solitary waves and a linear shear current

This paper is concerned with particle trajectories beneath solitary waves when a linear shear current exists. The fluid is assumed to be incompressible and inviscid, lying on a flat bed. Classical asymptotic expansion is used to obtain a Korteweg-de Vries(Kd V) equation, then a forth-order Runge-Kutta method is applied to get the approximate particle trajectories. On the other hand, our particular attention is paid to the direct numerical simulation(DNS) to the original Euler equations. A conformal map is used to solve the nonlinear boundary value problem. Highaccuracy numerical solutions are then obtained through the fast Fourier transform(FFT) and compared with the asymptotic solutions, which shows a good agreement when wave amplitude is small. Further, it also yields that there are different types of particle trajectories. Most surprisingly,periodic motion of particles could exist under solitary waves, which is due to the wave-current interaction.     

硬化材料中弹塑性柱面波的数值方法

本文讨论了应变硬化材料中一种适宜于数值计算的增量型塑性应力应变关系,并以弹塑性柱面波的传播为例进行了数值分析,给出了一些有意义的结果。     

DYNAMIC BEHAVIOR OF TWO PARALLEL SYMMETRY CRACKS IN MAGNETO-ELECTRO-ELASTIC COMPOSITES UNDER HARMONIC ANTI-PLANE WAVES

The dynamic behavior of two parallel symmetry cracks in magneto-electro-elastic composites under harmonic anti-plane shear waves is studied by Schmidt method. By using the Fourier transform, the problem can be solved with a pair of dual integral equations in which the unknown variable is the jumps of the displacements across the crack surfaces. To solve the dual integral equations, the jumps of the displacements across the crack surface were expanded in a series of Jacobi polynomials. The relations among the electric filed, the magnetic flux and the stress field were obtained. From the results, it can be obtained that the singular stresses in piezoelectric/piezomagnetic materials carry the same forms as those in a general elastic material for the dynamic anti-plane shear fracture problem. The shielding effect of two parallel cracks was also discussed.     

Higher order dispersive effects in regularized Boussinesq equation

In this paper, we consider the higher order Boussinesq (HBq) equation which models the bi-directional propagation of longitudinal waves in various continuous media. The equation contains the higher order effects of frequency dispersion. The present study is devoted to the numerical investigation of the HBq equation. For this aim a numerical scheme combining the Fourier pseudo-spectral method in space and a Runge–Kutta method in time is constructed. The convergence of semi-discrete scheme is proved in an appropriate Sobolev space. To investigate the higher order dispersive effects and nonlinear effects on the solutions of HBq equation, propagation of single solitary wave, head-on collision of solitary waves and blow-up solutions are considered.     

Instability and long-time evolution of cnoidal wave solutions for a Benney–Luke equation

We investigate numerically the stability of periodic traveling wave solutions (cnoidal waves) for a generalized Benney–Luke equation. By using a high-accurate Fourier spectral method, we find different kinds of evolution depending on the period of the perturbation. A cnoidal wave solution with period T is orbitally stable with regard to perturbations having the same period T, within certain range of wave velocities. This is a fact proved recently by Angulo and Quintero [Existence and orbital stability of cnoidal waves for a 1D boussinesq equation, International Journal of Mathematics and Mathematical Sciences (2007), in press, doi:10.1155/2007/52020] and our numerical experiments are consistent with their theory. In the present work we show numerically that cnoidal waves with period T become unstable when perturbed by small amplitude disturbances whose period is an integer multiple of T. Particularly, if the period of the perturbation is 2T, the evolution of the deviation of the solution from the orbit of the cnoidal wave is found to be approximately a time-periodic function. In other cases, the numerical experiments indicate a non-periodic behavior.     

Crack bifurcation predicted for dynamic anti-plane collinear cracks in piezoelectric materials using a non-local theory

A non-local theory of elasticity is applied to obtain the dynamic interaction between two collinear cracks in the piezoelectric materials plane under anti-plane shear waves for the permeable crack surface boundary conditions. Unlike the classical elasticity solution, a lattice parameter enters into the problem that make the stresses and the electric displacements finite at the crack tip. A one-dimensional non-local kernel is used instead of a two-dimensional one for the anti-plane dynamic problem to obtain the stress and electric displacement near the crack tips. By means of the Fourier transform, the problem can be solved with the help of two pairs of triple integral equations in which the unknown variable is the jump of the displacement across the crack surface. The solutions are obtained by means of the Schmidt method. Crack bifurcation is predicted using the strain energy density criterion. Minimum values of the strain energy density functions are assumed to coincide with the possible locations of fracture initiation. Bifurcation angles of ±5° and ±175° are found. The result of possible crack bifurcation was not expected before hand.     

Fourier analysis of an equal‐order incompressible flow solver stabilized by pressure gradient projection

Fourier analysis techniques are applied to the stabilized finite element method (FEM) recently proposed by Codina and Blasco for the approximation of the incompressible Navier–Stokes equations, here denoted by pressure gradient projection (SPGP) method. The analysis is motivated by spurious waves that pollute the computed pressure in start‐up flow simulation. An example of this spurious phenomenon is reported. It is shown that Fourier techniques can predict the numerical behaviour of stabilized methods with remarkable accuracy, even though the original Navier–Stokes setting must be significantly simplified to apply them. In the steady state case, good estimates for the stabilization parameters are obtained. In the transient case, spurious long waves are shown to be persistent when the element Reynolds number is large and the Courant number is small. This can be avoided by treating the pressure gradient projection implicitly, though this implies additional computing effort. Standard extrapolation variants are unfortunately unstable. Comparisons with Galerkin–least‐squares (GLS) method and Chorin's projection method are also addressed. Copyright © 2000 John Wiley & Sons, Ltd.     

AN ANALYSIS FOR DIFFRACTION OF FLEXURAL WAVES AND DYNAMIC STRESS CONCENTRATIONS IN THICK PLATES WITH A CAVITY

By using the complex variable method and conformal mapping,the diffraction of flexu-ral waves and dynamic stress concentrations in thick plates with a cavity have been studied.A generalsolution of the stress problem of the thick plate satisfying the boundary conditions on the contour of anarbitrary cavity is obtained.By employing the orthogonal function expansion technique,the dynamicstress problem can be reduced to the solution of an infinite algebraic equation series.As an example,the numerical results for the dynamic stress concentration factor in thick plates with a circular,ellipticcavity are graphically presented.The numerical results are discussed.     

Mean-field homogenization of elasto-viscoplastic composites based on a general incrementally affine linearization method

We propose a general formulation – which we believe to be new – for the mean-field homogenization of inclusion-reinforced elasto-viscoplastic composites assuming small strains. Our proposal is based on an interplay between constitutive equations and numerical algorithms, and the key ideas behind it are the following. The evolution equations for inelastic strain and internal variables at the beginning of each time interval are linearized around the ending time of the same interval. The linearized equations are then numerically integrated using a fully implicit backward Euler scheme. The obtained algebraic equations lead to an incrementally affine stress–strain relation which involves two important terms. The first one is the algorithmic tangent operator, obtained by consistent linearization of the time discretized constitutive equations. The second term is a new one and called an affine strain increment. The proposal leads to thermoelastic-like relations directly in the time domain, and not in the Laplace–Carson (L–C) one. There is no need for viscoelastic-type integral rewriting of the evolution equations, for L–C transforms, or for numerical inversion back from L–C to time domains. The proposed method can be readily applied to sophisticated elasto-viscoplastic models with an arbitrary set of scalar or tensor internal variables, and is valid for multi-axial, non-monotonic and non-proportional loading histories. The theory is applied in detail to a well-known constitutive model, and verified against finite element simulations of representative volume elements or unit cells, for a number of composite materials.     

Effects of SH waves in a functionally graded plate

A computational method is presented to investigate SH waves in functionally graded material (FGM) plates. The FGM plate is first divided into quadratic layer elements (QLEs), in which the material properties are assumed as a quadratic function in the thickness direction. A general solution for the equation of motion governing the QLE has been derived. The general solution is then used together with the boundary and continuity conditions to obtain the displacement and stress in the wave number domain for an arbitrary FGM plate. The displacements and stresses in the frequency domain and time domain are obtained using inverse Fourier integration. Furthermore, a simple integral technique is also proposed for evaluating modified Bessel functions with complex valued order. Numerical examples are presented to demonstrate this numerical technique for SH waves propagating in FGM plates.     

QALE‐FEM for modelling 3D overturning waves

A further development of the QALE‐FEM (quasi‐arbitrary Lagrangian–Eulerian finite element method) based on a fully nonlinear potential theory is presented in this paper. This development enables the QALE‐FEM to deal with three‐dimensional (3D) overturning waves over complex seabeds, which have not been considered since the method was devised by the authors of this paper in their previous works (J. Comput. Phys. 2006; 212 :52–72; J. Numer. Meth. Engng 2009; 78 :713–756). In order to tackle challenges associated with 3D overturning waves, two new numerical techniques are suggested. They are the techniques for moving the mesh and for calculating the fluid velocity near overturning jets, respectively. The developed method is validated by comparing its numerical results with experimental data and results from other numerical methods available in the literature. Good agreement is achieved. The computational efficiency of this method is also investigated for this kind of wave, which shows that the QALE‐FEM can be many times faster than other methods based on the same theory. Furthermore, 3D overturning waves propagating over a non‐symmetrical seabed or multiple reefs are simulated using the method. Some of these results have not been found elsewhere to the best of our knowledge. Copyright © 2009 John Wiley & Sons, Ltd.     

On the response of rubbers at high strain rates—II. Shock waves

In this series of papers, we examine the propagation of waves of finite deformation in rubbers through experiments and analysis; in the present paper, Part II, attention is focused on the propagation of one-dimensional tensile shock waves in strips of latex and nitrile rubber. Tensile wave propagation experiments were conducted at high strain rates by holding one end fixed and displacing the other end at a constant velocity. A high-speed video camera was used to monitor the motion and to determine the evolution of strain and particle velocity in rubber strips. Shock waves have been generated under tensile impact in prestretched rubber strips; analysis of the response yields the tensile shock adiabat for rubbers. The propagation of shocks is analyzed by developing an analogy with the theory of detonation; it is shown that the condition for shock propagation can be determined using the Chapman-Jouguet shock condition.