Ahmad Akbari R.¹, Akbar Bagri², Stéphane P. A. Bordas^3,4, Timon Rabczuk⁵

CMES-Computer Modeling in Engineering & Sciences, Vol.65, No.1, pp. 27-74, 2010, DOI:10.3970/cmes.2010.065.027

Abstract This contribution focuses on the simulation of two-dimensional elastic wave propagation in functionally graded solids and structures. Gradient volume fractions of the constituent materials are assumed to obey the power law function of position in only one direction and the effective mechanical properties of the material are determined by the Mori-Tanaka scheme. The investigations are carried out by extending a meshless method known as the Meshless Local Petrov-Galerkin (MLPG) method which is a truly meshless approach to thermo-elastic wave propagation. Simulations are carried out for rectangular domains under transient thermal loading. To investigate the effect of material composition on the… More >

P.H. Gaskell¹, Y.C. Lee², H.M. Thompson¹

CMES-Computer Modeling in Engineering & Sciences, Vol.62, No.1, pp. 77-112, 2010, DOI:10.3970/cmes.2010.062.077

Abstract The three-dimensional flow of a gravity-driven continuous thin liquid film on substrates containing micro-scale features is modelled using the long-wave lubrication approximation, encompassing cases when surface topography is either engulfed by the film or extends all the way though it. The discrete analogue of the time-dependent governing equations is solved accurately using a purpose designed multigrid strategy incorporating both automatic error-controlled adaptive time stepping and local mesh refinement/de-refinement. Central to the overall approach is a Newton globally convergent algorithm which greatly simplifies the inclusion of further physics via the solution of additional equations of the same form as the base… More >

C. Du¹, H. Bai², J. Qu³, X. Su^1,4

CMES-Computer Modeling in Engineering & Sciences, Vol.61, No.3, pp. 223-248, 2010, DOI:10.3970/cmes.2010.061.223

Abstract Second, or higher, order harmonics have great potential in fatigue life prediction. In this study, the dispersion properties of waves propagating in the nonlinear random media are investigated. An one dimensional nonlinear model based on the nonlinear Hikata stress-strain relation is used. We applied perturbation method, the Liouville transformation and the smoothing approximation method to solve the one dimensional nonlinear stochastic wave equation. We show easily that the dispersion equations for all higher order terms will be the same with the corresponding linear random medium by perturbation method. The linear stochastic equation with two random coefficients is greatly simplified to… More >

Bart Bergen¹, Bert Van Genechten¹, Dirk Vandepitte¹, Wim Desmet¹

CMES-Computer Modeling in Engineering & Sciences, Vol.61, No.2, pp. 155-176, 2010, DOI:10.3970/cmes.2010.061.155

Abstract The Wave Based Method (WBM) is a numerical prediction technique for Helmholtz problems. It is an indirect Trefftz method using wave functions, which satisfy the Helmholtz equation, for the description of the dynamic variables. In this way, it avoids both the large systems and the pollution errors that jeopardize accurate element-based predictions in the mid-frequency range. The enhanced computational efficiency of the WBM as compared to the element-based methods has been proven for the analysis of both three-dimensional bounded and two-dimensional unbounded problems. This paper presents an extension of the WBM to the application of three-dimensional acoustic scattering and radiation… More >

A. Frangi¹, M. Bonnet²

CMES-Computer Modeling in Engineering & Sciences, Vol.58, No.3, pp. 271-296, 2010, DOI:10.3970/cmes.2010.058.271

Abstract This paper presents an empirical study of the accuracy of multipole expansions of Helmholtz-like kernels with complex wavenumbers of the form k = (α + iβ)ϑ, with α = 0,±1 and β > 0, which, the paucity of available studies notwithstanding, arise for a wealth of different physical problems. It is suggested that a simple point-wise error indicator can provide an a-priori indication on the number N of terms to be employed in the Gegenbauer addition formula in order to achieve a prescribed accuracy when integrating single layer potentials over surfaces. For β ≥ 1 it is observed that the… More >

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 >

J.T. Zhou¹, Q.W. Ma^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.56, No.2, pp. 179-210, 2010, DOI:10.3970/cmes.2010.056.179

Abstract In this paper, the Meshless Local Petrov-Galerkin method based on Rankine source solution (MLPG_R) is further developed to model 3D breaking waves. For this purpose, the technique for identifying free surface particles called Mixed Particle Number Density and Auxiliary Function Method (MPAM) and the semi-analytical technique for estimating the domain integrals for 2D cases are extended to 3D cases. In addition, a new semi-analytical technique is developed to deal with the local spherical surface integrals. The numerical results obtained by the newly developed method will be compared with experimental data available in literature and satisfactory agreement will be shown. 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 >

J. Trevelyan¹and G. Coates¹

CMES-Computer Modeling in Engineering & Sciences, Vol.55, No.2, pp. 147-170, 2010, DOI:10.3970/cmes.2010.055.147

Abstract The terminology "wave boundary elements" relates to boundary elements enriched in the Partition of Unity sense by a multiple plane wave basis for the analysis of the propagation of short wavelength waves. This paper presents a variant of this approach in which the plane wave basis is selected adaptively according to an error indicator. The error indicator is residual based, and exhibits useful local and global properties. Model improvement in each adaptive iteration is carried out by the addition of new plane waves with no h-refinement. The convergence properties of the scheme are demonstrated. 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 >