M. Kashtalyan^1,2, Y.A. Zhuk³

CMES-Computer Modeling in Engineering & Sciences, Vol.57, No.1, pp. 1-30, 2010, DOI:10.3970/cmes.2010.057.001

Abstract The paper investigates propagation of stationary plane longitudinal and transverse waves along the layers in adhesively bonded multilayered structures for MEMS applications in the presence of residual stresses. The multilayered structure is assumed to consist of the infinite amount of the periodically recurring layers made of two different materials possessing significantly dissimilar properties: conductive metal layer and insulating adhesive layer. It is assumed that the mechanical behaviour of both materials is nonlinear elastic and can be described with the help of the elastic Murnaghan potential depending on the three invariants of strain tensor. The problem is formulated in the framework… More >

R. Martin¹, D. Komatitsch^1,2, S. D. Gedney³, E. Bruthiaux^1,4

CMES-Computer Modeling in Engineering & Sciences, Vol.56, No.1, pp. 17-42, 2010, DOI:10.3970/cmes.2010.056.017

Abstract Unsplit convolutional perfectly matched layers (CPML) for the velocity and stress formulation of the seismic wave equation are classically computed based on a second-order finite-difference time scheme. However it is often of interest to increase the order of the time-stepping scheme in order to increase the accuracy of the algorithm. This is important for instance in the case of very long simulations. We study how to define and implement a new unsplit non-convolutional PML called the Auxiliary Differential Equation PML (ADE-PML), based on a high-order Runge-Kutta time-stepping scheme and optimized at grazing incidence. We demonstrate that when a second-order time-stepping… More >

Surkay D. Akbarov^1,2,3, Mugan S. Guliev⁴, Ramazan Tekercioglu⁵

CMES-Computer Modeling in Engineering & Sciences, Vol.55, No.1, pp. 1-32, 2010, DOI:10.3970/cmes.2010.055.001

Abstract Dispersion relations of axisymmetric longitudinal wave propagation in a finite pre-strained compound (bi-material) cylinder made from high elastic incompressible materials are investigated within the scope of a piecewise homogeneous body model utilizing three-dimensional linearized theory wave propagation in the initially stressed body. The materials of the inner and outer cylinders are assumed to be neo-Hookean. The numerical results regarding the influence of the initial strains in the inner and outer cylinders on the wave dispersion are presented and discussed. These results are obtained for the case where the material of the inner solid cylinder is stiffer than that of the… More >

Delfim Soares Jr. ¹

CMES-Computer Modeling in Engineering & Sciences, Vol.50, No.2, pp. 97-114, 2009, DOI:10.3970/cmes.2009.050.097

Abstract In this work, meshless methods based on the local Petrov-Galerkin (MLPG) approach are presented to analyse electromagnetic wave propagation problems. Formulations adopting the Heaviside step function and the Gaussian weight function as the test functions in the local weak form are considered. The moving least square (MLS) method is used to approximate the physical quantities in the local integral equations. After spatial discretization is carried out, a system of ordinary differential equations of second order is obtained. This system is solved in the time-domain by the Houbolt's method, allowing the computation of the so-called primary fields (either the electric or… More >

Ronan Madec¹, Dimitri Komatitsch^1,2, Julien Diaz³

CMES-Computer Modeling in Engineering & Sciences, Vol.49, No.2, pp. 163-190, 2009, DOI:10.3970/cmes.2009.049.163

Abstract When studying seismic wave propagation in fluid-solid models based on a numerical technique in the time domain with an explicit time scheme it is often of interest to resort to time substepping because the stability condition in the solid part of the medium can be more stringent than in the fluid. In such a case, one should enforce the conservation of energy along the fluid-solid interface in the time matching algorithm in order to ensure the accuracy and the stability of the time scheme. This is often not done in the available literature and approximate techniques that do not enforce… More >

A. Chakraborty¹

CMES-Computer Modeling in Engineering & Sciences, Vol.40, No.2, pp. 105-132, 2009, DOI:10.3970/cmes.2009.040.105

Abstract A mathematical model is presented in this work that describes the behavior of porous piezoelectric materials subjected to mechanical load and electric field. The model combines Biot's theory of poroelasticity and the classical theory of piezoelectric material wherein it is assumed that piezoelectric coupling exists only with the solid phase of the porous medium. This model is used to analyze the stress and electric wave generated in bone and porous Lead-Zirconate-Titanate (PZT) due to high frequency pulse loading. The governing partial differential equations are solved in the frequency domain by transforming them into a polynomial eigenvalue structure. This approach permits… More >

Surkay D. Akbarov^1,2,3, Mugan S. Guliev⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.39, No.2, pp. 155-178, 2009, DOI:10.3970/cmes.2009.039.155

Abstract The axisymmetric longitudinal wave propagation in a finite pre-strained compound (composite) cylinder is investigated within the scope of a piecewise homogeneous body model utilizing three-dimensional linearized theory wave propagation in an initially stressed body. The materials of the inner and outer cylinder are assumed to be compressible. The elasticity relations for those are given through the harmonic potential. The algorithm for constructing of the computer programmes and obtaining numerical results is discussed. The numerical results regarding the influences of the initial strains in the inner and outer cylinders on the wave dispersion are presented and analysed. These results are obtained… More >

A.I. Abreu¹, W.J. Mansur¹, D. Soares Jr^1,2, J.A.M. Carrer³

CMES-Computer Modeling in Engineering & Sciences, Vol.30, No.3, pp. 123-132, 2008, DOI:10.3970/cmes.2008.030.123

Abstract This paper is concerned with the numerical computation of space derivatives of a time-domain (TD-) Boundary Element Method (BEM) formulation for the analysis of scalar wave propagation problems. In the present formulation, the Convolution Quadrature Method (CQM) is adopted, i.e., the basic integral equation of the TD-BEM is numerically substituted by a quadrature formula, whose weights are computed using the Laplace transform of the fundamental solution and a linear multi-step method. In order to numerically compute space derivatives, the present work properly transforms the quadrature weights of the CQM-BEM, adopting the so-called Complex-Variable-Differentiation Method (CVDM). Numerical examples are presented at… More >

Mira Mitra¹, S. Gopalakrishnan²

CMES-Computer Modeling in Engineering & Sciences, Vol.27, No.1&2, pp. 125-136, 2008, DOI:10.3970/cmes.2008.027.125

Abstract In this paper, the wave characteristics, namely, the spectrum and dispersion relations of multi-wall carbon nanotubes (MWNTs) are studied. The MWNTs are modeled as multiple thin shells coupled through van der Waals force. Each wall of the MWNT has three displacements, i.e, axial, circumferential and radial with variation along the axial and circumferential directions. The wave characteristics are obtained by transforming the governing differential wave equations to frequency domain via Fourier transform. This transformation is first done in time using fast Fourier transform (FFT) and then in one spatial dimension using Fourier series. These transformed equations are solved by posing… More >

P. Le¹, N. Mai-Duy², T. Tran-Cong³, G. Baker⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.25, No.1, pp. 43-68, 2008, DOI:10.3970/cmes.2008.025.043

Abstract This paper describes an integrated radial basis function network (IRBFN) method for the numerical modelling of the dynamics of strain localization due to strain softening in quasi-brittle materials. The IRBFN method is a truly meshless method that is based on an unstructured point collocation procedure. We introduce a new and effective regularization method to enhance the performance of the IRBFN method and alleviate the numerical oscillations associated with weak discontinuity at the elastic wave front. The dynamic response of a one dimensional bar is investigated using both local and non-local continuum models. Numerical results, which compare favourably with those obtained… More >