Non‐local dispersive model for wave propagation in heterogeneous media: one‐dimensional case

Non‐local dispersive model for wave propagation in heterogeneous media is derived from the higher‐order mathematical homogenization theory with multiple spatial and temporal scales. In addition to the usual space–time co‐ordinates, a fast spatial scale and a slow temporal scale are introduced to account for rapid spatial fluctuations of material properties as well as to capture the long‐term behaviour of the homogenized solution. By combining various order homogenized equations of motion the slow time dependence is eliminated giving rise to the fourth‐order differential equation, also known as a ‘bad’ Boussinesq problem. Regularization procedures are then introduced to construct the so‐called ‘good’ Boussinesq problem, where the need for C¹ continuity is eliminated. Numerical examples are presented to validate the present formulation. Copyright © 2002 John Wiley & Sons, Ltd.     

Stabilized nonlocal model for dispersive wave propagation in heterogeneous media

In this paper, a general-purpose computational model for dispersive wave propagation in heterogeneous media is developed. The model is based on the higher-order homogenization with multiple spatial and temporal scales and the C⁰-continuous mixed finite element approximation of the resulting nonlocal equations of motion. The proposed nonlocal Hamilton principle leads to the stable discrete system of equations independent of the mesh size, unit cell domain and the excitation frequency. The method has been validated for plane harmonic analysis and for transient wave motion insemi-infinite domain with various microstructures.This work was supported by the Sandia National Laboratories under Contract DE-AL04–94AL8500, the Office of Naval Research through grant number N00014–97–1-0687, and the Japan Society for the Promotion of Science under contract number Heisei11-nendo 06542.     

Time‐domain hybrid formulations for wave simulations in three‐dimensional PML‐truncated heterogeneous media

We are concerned with the numerical simulation of wave motion in arbitrarily heterogeneous, elastic, perfectly‐matched‐layer‐(PML)‐truncated media. We extend in three dimensions a recently developed two‐dimensional formulation, by treating the PML via an unsplit‐field, but mixed‐field, displacement‐stress formulation, which is then coupled to a standard displacement‐only formulation for the interior domain, thus leading to a computationally cost‐efficient hybrid scheme. The hybrid treatment leads to, at most, third‐order in time semi‐discrete forms. The formulation is flexible enough to accommodate the standard PML, as well as the multi‐axial PML. We discuss several time‐marching schemes, which can be used à la carte, depending on the application: (a) an extended Newmark scheme for third‐order in time, either unsymmetric or fully symmetric semi‐discrete forms; (b) a standard implicit Newmark for the second‐order, unsymmetric semi‐discrete forms; and (c) an explicit Runge–Kutta scheme for a first‐order in time unsymmetric system. The latter is well‐suited for large‐scale problems on parallel architectures, while the second‐order treatment is particularly attractive for ready incorporation in existing codes written originally for finite domains. We compare the schemes and report numerical results demonstrating stability and efficacy of the proposed formulations. Copyright © 2014 John Wiley & Sons, Ltd.     

Non‐local damage model based on displacement averaging

Continuum damage models describe the changes of material stiffness and strength, caused by the evolution of defects, in the framework of continuum mechanics. In many materials, a fast evolution of defects leads to stress–strain laws with softening, which creates serious mathematical and numerical problems. To regularize the model behaviour, various generalized continuum theories have been proposed. Integral‐type non‐local damage models are often based on weighted spatial averaging of a strain‐like quantity. This paper explores an alternative formulation with averaging of the displacement field. Damage is assumed to be driven by the symmetric gradient of the non‐local displacements. It is demonstrated that an exact equivalence between strain and displacement averaging can be achieved only in an unbounded medium. Around physical boundaries of the analysed body, both formulations differ and the non‐local displacement model generates spurious damage in the boundary layers. The paper shows that this undesirable effect can be suppressed by an appropriate adjustment of the non‐local weight function. Alternatively, an implicit gradient formulation could be used. Issues of algorithmic implementation, computational efficiency and smoothness of the resolved stress fields are discussed. Copyright © 2005 John Wiley & Sons, Ltd.     

Non‐local boundary integral formulation for softening damage

A strongly non‐local boundary element method (BEM) for structures with strain‐softening damage treated by an integral‐type operator is developed. A plasticity model with yield limit degradation is implemented in a boundary element program using the initial‐stress boundary element method with iterations in each load increment. Regularized integral representations and boundary integral equations are used to avoid the difficulties associated with numerical computation of singular integrals. A numerical example is solved to verify the physical correctness and efficiency of the proposed formulation. The example consists of a softening strip perforated by a circular hole, subjected to tension. The strain‐softening damage is described by a plasticity model with a negative hardening parameter. The local formulation is shown to exhibit spurious sensitivity to cell mesh refinements, localization of softening damage into a band of single‐cell width, and excessive dependence of energy dissipation on the cell size. By contrast, the results for the non‐local theory are shown to be free of these physically incorrect features. Compared to the classical non‐local finite element approach, an additional advantage is that the internal cells need to be introduced only within the small zone (or band) in which the strain‐softening damage tends to localize within the structure. Copyright © 2003 John Wiley & Sons, Ltd.     

Equivalent continuum modeling for non‐linear wave propagation in jointed media

This paper describes a methodology to build equivalent continuum (ECC) models to represent jointed media for non‐linear wave propagation. In the present work, we study the response of randomly jointed rock media both to quasi‐static and dynamic loadings. We have found that, unlike in the quasi‐static case where no relaxation is observed, during dynamic loading, the deviatoric stress drops in the plastic wave after reaching its peak. Such stress relaxation phenomenon is calculated in jointed media, although the model used for intact rock blocks is rate‐independent. This is a strong indication that EC models for jointed rock should include some rate‐dependence when used in wave propagation simulations. Copyright © 2010 John Wiley & Sons, Ltd.     

A computational approach for multi‐scale analysis of heterogeneous elasto‐plastic media subjected to short duration loads

The asymptotic expansion homogenization (AEH) approach has found wide acceptance for the study of heterogeneous structures due to its ability to account for multi‐scale features. The emphasis of the present study is to develop consistent AEH numerical formulations to address elasto‐plastic material response of structures subjected to short‐duration transient loading. A second‐order accurate velocity‐based explicit time integration method, in conjunction with the AEH approach, is currently developed that accounts for large deformation non‐linear material response. The approach is verified under degenerate homogeneous conditions using existing experimental data in the literature and its ability to account for heterogeneous conditions is demonstrated for a number of test problems. Copyright © 2004 John Wiley & Sons, Ltd.     

Non‐reflecting artificial boundaries for transient scalar wave propagation in a two‐dimensional infinite homogeneous layer

This paper presents an exact non‐reflecting boundary condition for dealing with transient scalar wave propagation problems in a two‐dimensional infinite homogeneous layer. In order to model the complicated geometry and material properties in the near field, two vertical artificial boundaries are considered in the infinite layer so as to truncate the infinite domain into a finite domain. This treatment requires the appropriate boundary conditions, which are often referred to as the artificial boundary conditions, to be applied on the truncated boundaries. Since the infinite extension direction is different for these two truncated vertical boundaries, namely one extends toward x →∞ and another extends toward x→‐ ∞, the non‐reflecting boundary condition needs to be derived on these two boundaries. Applying the variable separation method to the wave equation results in a reduction in spatial variables by one. The reduced wave equation, which is a time‐dependent partial differential equation with only one spatial variable, can be further changed into a linear first‐order ordinary differential equation by using both the operator splitting method and the modal radiation function concept simultaneously. As a result, the non‐reflecting artificial boundary condition can be obtained by solving the ordinary differential equation whose stability is ensured. Some numerical examples have demonstrated that the non‐reflecting boundary condition is of high accuracy in dealing with scalar wave propagation problems in infinite and semi‐infinite media. Copyright © 2003 John Wiley & Sons, Ltd.     

Practical aspects of higher‐order numerical schemes for wave propagation phenomena

This paper examines practical issues related to the use of compact‐difference‐based fourth‐ and sixth‐order schemes for wave propagation phenomena with focus on Maxwell's equations of electromagnetics. An outline of the formulation and scheme optimization is followed by an assessment of the error accruing from application on stretched meshes with two approaches: transformed plane method and physical space differencing. In the first technique, the truncation error expansion for the sixth‐order compact scheme confirms that the order of accuracy is preserved if a consistent mesh refinement strategy is followed and further that metrics should be evaluated numerically even if analytic expressions are available. Physical space‐differencing formulas are derived for the five‐point stencil by expressing the coefficients in terms of local spacing ratios. The order of accuracy of the reconstruction operator is then verified with a numerical experiment on stretched meshes. To ensure stability for a broad range of problems, Fourier analysis is employed to develop a single‐parameter family of up to tenth‐order tridiagonal‐based spatial filters. The implementation of these filters is discussed in terms of their effect on the interior scheme as well as in a 1‐D cavity where they are employed to suppress a late‐time instability. The paper concludes after demonstrating the application of the scheme to several 3‐D canonical problems utilizing Cartesian as well as curvilinear meshes. Published in 1999 by John Wiley & Sons, Ltd. This article is a U.S. Government work and is in the public domain in the United States.     

Finite element formulation of exact non‐reflecting boundary conditions for the time‐dependent wave equation

A modified version of an exact Non‐reflecting Boundary Condition (NRBC) first derived by Grote and Keller is implemented in a finite element formulation for the scalar wave equation. The NRBC annihilate the first N wave harmonics on a spherical truncation boundary, and may be viewed as an extension of the second‐order local boundary condition derived by Bayliss and Turkel. Two alternative finite element formulations are given. In the first, the boundary operator is implemented directly as a ‘natural’ boundary condition in the weak form of the initial–boundary value problem. In the second, the operator is implemented indirectly by introducing auxiliary variables on the truncation boundary. Several versions of implicit and explicit time‐integration schemes are presented for solution of the finite element semidiscrete equations concurrently with the first‐order differential equations associated with the NRBC and an auxiliary variable. Numerical studies are performed to assess the accuracy and convergence properties of the NRBC when implemented in the finite element method. The results demonstrate that the finite element formulation of the (modified) NRBC is remarkably robust, and highly accurate. Copyright © 1999 John Wiley & Sons, Ltd.     

Two‐scale diffusion–deformation coupling model for material deterioration involving micro‐crack propagation

In this paper, we introduce a two‐scale diffusion–deformation coupled model that represents the aging material deterioration of two‐phase materials involving micro‐crack propagations. The mathematical homogenization method is applied to relate the micro‐ and macroscopic field variables, and a weak coupling solution method is employed to solve the two‐way coupling phenomena between the diffusion of scalar fields and the deformation of quasi‐brittle solids. The macroscopic mechanical behavior represented by the derived two‐scale two‐way coupled model reveals material nonlinearity due to micro‐scale cracking induced by the scalar‐field‐induced deformation, which can be simulated by the finite cover method. After verifying the fundamental validity of the proposed model and the analysis method, we perform a simple numerical example to demonstrate their ability to predict aging material deterioration. Copyright © 2010 John Wiley & Sons, Ltd.     

A simple multi‐directional absorbing layer method to simulate elastic wave propagation in unbounded domains

The numerical analysis of elastic wave propagation in unbounded media may be difficult due to spurious waves reflected at the model artificial boundaries. This point is critical for the analysis of wave propagation in heterogeneous or layered solids. Various techniques such as Absorbing Boundary Conditions, infinite elements or Absorbing Boundary Layers (e.g. Perfectly Matched Layers) lead to an important reduction of such spurious reflections. In this paper, a simple absorbing layer method is proposed: it is based on a Rayleigh/Caughey damping formulation which is often already available in existing Finite Element softwares. The principle of the Caughey Absorbing Layer Method is first presented (including a rheological interpretation). The efficiency of the method is then shown through 1D Finite Element simulations considering homogeneous and heterogeneous damping in the absorbing layer. 2D models are considered afterwards to assess the efficiency of the absorbing layer method for various wave types and incidences. A comparison with the PML method is first performed for pure P‐waves and the method is shown to be reliable in a more complex 2D case involving various wave types and incidences. It may thus be used for various types of problems involving elastic waves (e.g. machine vibrations, seismic waves, etc.). Copyright © 2010 John Wiley & Sons, Ltd.     

A method for computation of discontinuous wave propagation in heterogeneous solids: basic algorithm description and application to one‐dimensional problems

An explicit integration algorithm for computations of discontinuous wave propagation in heterogeneous solids is presented, which is aimed at minimizing spurious oscillations when the wave fronts pass through several zones of different wave speeds. The essence of the present method is a combination of two wave capturing characteristics: a new integration formula that is obtained by pushforward–pullback operations in time designed to filter post‐shock oscillations, and the central difference method that intrinsically filters front‐shock oscillations. It is shown that a judicious combination of these two characteristics substantially reduces both spurious front‐shock and post‐shock oscillations. The performance of the new method is demonstrated as applied to wave propagation through a uniform bar with varying courant numbers, then to heterogeneous bars. Copyright © 2012 John Wiley & Sons, Ltd.     

Discrete and continuum methods for numerical simulations of non‐linear wave propagation in discontinuous media

The main objective of this paper is to develop a methodology to model dynamic loading of various discontinuous media. Two methods for modeling non‐linear waves in a media with multiple discontinuities are considered. The first one, a discrete method, is based on the Simple Common Plane contact algorithm. This method can be applied both to compliant contacts characterized by finite thickness and elastic moduli (such as joints in geomechanics) as well as to non‐compliant frictional contacts traditionally described by the slide lines in finite element/finite difference codes. The second one, a continuum method, assumes that the contacts are not compliant and can be modeled as one or several weakness planes cutting through the elements of the computational mesh. Both discrete and continuum methods described in the paper can be applied to derive equivalent continuum properties for media with multiple discontinuities. An example of such application for a randomly jointed media is given in the paper. Copyright © 2010 John Wiley & Sons, Ltd.     

Thermo‐mechanical analysis of periodic multiphase materials by a multiscale asymptotic homogenization approach

A spatial and temporal multiscale asymptotic homogenization method used to simulate thermo‐dynamic wave propagation in periodic multiphase materials is systematically studied. A general field governing equation of thermo‐dynamic wave propagation is expressed in a unified form with both inertia and velocity terms. Amplified spatial and reduced temporal scales are, respectively, introduced to account for spatial and temporal fluctuations and non‐local effects in the homogenized solution due to material heterogeneity and diverse time scales. The model is derived from the higher‐order homogenization theory with multiple spatial and temporal scales. It is also shown that the modified higher‐order terms bring in a non‐local dispersion effect of the microstructure of multiphase materials. One‐dimensional non‐Fourier heat conduction and dynamic problems under a thermal shock are computed to demonstrate the efficiency and validity of the developed procedure. The results indicate the disadvantages of classical spatial homogenization. Copyright © 2006 John Wiley & Sons, Ltd.     

Three‐dimensional Coons macroelements: application to eigenvalue and scalar wave propagation problems

This paper discusses the matter of using higher order approximations in three‐dimensional problems through Coons macroelements. Recently, we have proposed a global functional set based on ‘Coons interpolation formula’ for the construction of large two‐dimensional macroelements with degrees of freedom appearing at the boundaries only of the domain. After successive application in many engineering problems, this paper extends the methodology to large three‐dimensional hexahedral macroelements with the degrees of freedom appearing at the 12 edges of the entire domain in case of smooth box‐like structures. Closed‐form expressions of the global shape functions are presented for the first time. It is shown that these global shape functions can be automatically constructed in a systematic way by arbitrarily choosing univariate approximations such as piecewise‐linear, cubic B‐splines, Lagrange polynomials, etc., along the control lines. Moreover, the mechanism of adding facial and internal nodes is presented. Relationships with higher order p‐methods are discussed. Following to excellent results previously derived for the solution of the Laplace equation as well as static and eigenvalue extraction analysis of structures, the paper investigates the performance of Coons macroelements in 3‐D eigenvalue and scalar wave propagation problems by implementing standard time‐integration schemes. Copyright © 2005 John Wiley & Sons, Ltd.     

Non‐reflecting boundary conditions for atomistic,continuum and coupled atomistic/continuum simulations

We present a method to numerically calculate a non‐reflecting boundary condition which is applicable to atomistic, continuum and coupled multiscale atomistic/continuum simulations. The method is based on the assumption that the forces near the domain boundary can be well represented as a linear function of the displacements, and utilizes standard Laplace and Fourier transform techniques to eliminate the unnecessary degrees of freedom. The eliminated degrees of freedom are accounted for in a time‐history kernel that can be calculated for arbitrary crystal lattices and interatomic potentials, or regular finite element meshes using an automated numerical procedure. The new theoretical developments presented in this work allow the application of the method to non‐nearest neighbour atomic interactions; it is also demonstrated that the identical procedure can be used for finite element and mesh‐free simulations. We illustrate the effectiveness of the method on a one‐dimensional model problem, and calculate the time‐history kernel for FCC gold using the embedded atom method (EAM). Copyright © 2005 John Wiley & Sons, Ltd.     

Coupled multi‐scale cohesive modeling of failure in heterogeneous adhesives

A multi‐scale cohesive numerical framework is proposed to simulate the failure of heterogeneous adhesively bonded systems. This multi‐scale scheme is based on Hill's variational principle of energy equivalence between the higher and lower level scales. It provides an easy way to obtain accurate homogenized macroscopic properties while capturing the physics of failure processes at the micro‐scale in sufficient detail. We use an isotropic rate‐dependent damage model to mimic the failure response of the constituents of heterogeneous adhesives. The finite element method is used to solve the equilibrium equation at each scale. A nested iterative scheme inspired by the return mapping algorithm used in computational inelasticity is implemented. We propose a computationally attractive technique to couple the macro‐ and micro‐scales for rate‐dependent constitutive laws. We introduce an adhesive patch test to study the numerical performance, including spatial and temporal convergence of the multi‐scale scheme. We compare the solution of the multi‐scale cohesive scheme with a direct numerical simulation. Finally, we solve mode I and mode II fracture problems to demonstrate failure at the macro‐scale. Copyright © 2010 John Wiley & Sons, Ltd.