K. Han¹, Y. T. Feng¹, D. R. J. Owen¹

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

Abstract A computational framework is established for effective modelling of fluid-thermal-particle interactions. The numerical procedures comprise the Discrete Element Method for simulating particle dynamics; the Lattice Boltzmann Method for modelling the mass and velocity field of the fluid flow; and the Discrete Thermal Element Method and the Thermal Lattice Boltzmann Method for solving the temperature field. The coupling of the three fields is realised through hydrodynamic interaction force terms. Selected numerical examples are provided to illustrate the applicability of the proposed approach. More >

M. R. Swager¹, L. J. Gray², S. Nintcheu Fata²

CMES-Computer Modeling in Engineering & Sciences, Vol.58, No.3, pp. 297-312, 2010, DOI:10.3970/cmes.2010.058.297

Abstract A linear element Galerkin boundary integral analysis for the three-dimensional Helmholtz equation is presented. The emphasis is on solving acoustic scattering by an open (crack) surface, and to this end both a dual equation formulation and a symmetric hypersingular formulation have been developed. All singular integrals are defined and evaluated via a boundary limit process, facilitating the evaluation of the (finite) hypersingular Galerkin integral. This limit process is also the basis for the algorithm for post-processing of the surface gradient. The analytic integrations required by the limit process are carried out by employing a Taylor 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 More >

L. Távara¹, V. Manticˇ ¹, E. Graciani¹, J. Cañas¹, F. París¹

CMES-Computer Modeling in Engineering & Sciences, Vol.58, No.3, pp. 247-270, 2010, DOI:10.3970/cmes.2010.058.247

Abstract The problem of a crack in a thin adhesive layer is considered. The adherents may have orthotropic elastic behavior which allows composite laminates to be modeled. In the present work a linear elastic-brittle constitutive law of the thin adhesive layer, called weak interface model, is adopted, allowing an easy modeling of crack propagation along it. In this law, the normal and tangential stresses across the undamaged interface are proportional to the relative normal and tangential displacements, respectively. Interface crack propagation is modeled by successive breaking of the springs used to discretize the weak interface. An… More >

J. Tausch¹, J. Xiao²

CMES-Computer Modeling in Engineering & Sciences, Vol.58, No.3, pp. 221-246, 2010, DOI:10.3970/cmes.2010.058.221

Abstract A fast method for the computation of layer potentials that arise in acoustic scattering is introduced. The principal idea is to split the singular kernel into a smooth and a local part. The potential due to the smooth part is discretized by a Nyström method and is evaluated efficiently using a sequence of FFTs. The potential due to the local part is approximated by a truncated series in the mollification parameter. The smooth approximation of the kernel is obtained by multiplication of its Fourier transform with a filter. We will show that for a rational More >

Y.C. Cai^1,2, J.K. Paik³, S.N. Atluri⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.58, No.1, pp. 75-108, 2010, DOI:10.3970/cmes.2010.058.075

Abstract This paper presents a new shear flexible beam/rod element for large deformation analyses of space-frame structures comprising of thin or thick members, based on the Reissner variational principle and a von Karman type nonlinear theory of deformation in the co-rotational reference frame of the present beam element. The C0continuous trial functions for transverse rotations in two independent directions are used over each element, to derive an explicit expression for the (16x16)symmetrictangent stiffness matrix of the beam element in the co-rotational reference frame. When compared to the primal approach wherein C1continuous trial functions for transverse displacements… More >

A. Papacharalampopoulos², G. F. Karlis², A. Charalambopoulos³, D. Polyzos⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.58, No.1, pp. 45-74, 2010, DOI:10.3970/cmes.2010.058.045

Abstract A Boundary Element Method (BEM) for solving two (2D) and three dimensional (3D) dynamic problems in materials with microstructural effects is presented. The analysis is performed in the frequency domain and in the context of Mindlin's Form II gradient elastic theory. The fundamental solution of the differential equation of motion is explicitly derived for both 2D and 3D problems. The integral representation of the problem, consisting of two boundary integral equations, one for displacements and the other for its normal derivative is exploited for the proposed BEM formulation. The global boundary of the analyzed domain More >

Lorenzo Codecasa¹, Francesco Trevisan²

CMES-Computer Modeling in Engineering & Sciences, Vol.58, No.1, pp. 15-44, 2010, DOI:10.3970/cmes.2010.058.015

Abstract This paper starts from the spatial discretization of an electromagnetic problem over pairs of oriented grids, one dual of the other, according to the so called Discrete Geometric Approach(DGA) to computational electromagnetism; the Cell Method or the Finite Integration Technique are examples of such an approach. The core of the work is providing for the first time a convergence analysis when the discrete counter-parts of constitutive relations are computed by means of an energetic framework. More >

Yunhua Li^1,2, Yunze Li³, Chieh-Li Chen⁴, Cha’o-Kuang Chen⁵

CMES-Computer Modeling in Engineering & Sciences, Vol.58, No.1, pp. 1-14, 2010, DOI:10.3970/cmes.2010.058.001

Abstract Addressing the identification problem of the general Lagrange multiplier in the He's variational iteration method, this paper proposes a new kind of method based on Duhamel's principle for the dynamic system response analysis. In this method, we have constructed an analytical iteration formula in terms of the convolution for the residual error at the nth iteration, and have given a new interpretation to He's variational iteration method. The analysis illustrates that the computational result of this method is equal to that of He's variational iteration method on the assumption of considering the impulse response of More >

S. Ullah¹, B. Pichler¹, S. Scheiner^1,2, C. Hellmich^1,3

CMES-Computer Modeling in Engineering & Sciences, Vol.57, No.3, pp. 279-316, 2010, DOI:10.3970/cmes.2010.057.279

Abstract Point-wise optical measurements of 3D displacement vectors over time are a key input for monitoring shotcrete tunnel shells during construction according to the New Austrian Tunnelling Method (NATM). Aiming at estimation of the stresses prevailing in the highly loaded, hydrating material, we here deal with two different interpolation strategies for reconstructing, from measured displacement vectors, the 3D displacement field histories of the inner surface of the tunnel shell. The first approach considers spatial interpolation of displacement components in a fixed Cartesian base frame, while the second (new) approach refers to displacement components in a moving… More >