J. X. Zhou¹, T. Koziara¹, T. G. Davies¹

CMES-Computer Modeling in Engineering & Sciences, Vol.16, No.2, pp. 131-140, 2006, DOI:10.3970/cmes.2006.016.131

Abstract The classical BEM approach for elastodynamics can produce poor results when high gradients are generated by impulses. High gradient areas evolve over time and their locations are unknown a priori, so they usually can not be captured by uniform meshes. In this paper, we propose a novel method which interpolates both spatial and temporal domains. A direct space-time discretization scheme is used to capture the wave fronts accurately and to forestall generation of spurious oscillations there. Some numerical examples are given to demonstrate the power and scope of the method. More >

Shiwei Zhou¹, Michael Yu Wang²

CMES-Computer Modeling in Engineering & Sciences, Vol.16, No.2, pp. 83-102, 2006, DOI:10.3970/cmes.2006.016.083

Abstract This paper describes a self-mass-conservative Cahn-Hilliard (C-H) model with elastic strain energy (mean compliance) for the optimization of multi-material structure topology. The total free energy of the generalized C-H system can be represented as a Lyapunov functional so that the elastic strain energy, as a part of the total free energy, decreases gradually to attain optimal material distribution. The interface energy relating to phase gradient in the total free energy plays an important role in regularizing the original ill-posed problem by restricting the structure's boundaries. On the other hand, interface coalescence and break-up due to phase separation and grain coarsening… More >

Chun-Ying Wu^1,3, Xue-Yan Liu², A. Scarpas², Xiu-Run Ge³

CMES-Computer Modeling in Engineering & Sciences, Vol.12, No.2, pp. 149-158, 2006, DOI:10.3970/cmes.2006.012.149

Abstract For the spectral analysis of the three-dimensional multi-layered pavement, 3D layer spectral element method is presented to solve the problems of bounded layer system subjected to a transient load pulse. In spectral element, each layer is treated as one spectral element. The wave propagation inside each layer element is achieved by the superposition of the incident wave and the reflection wave. Fast Fourier transformation is used to transform FWD datum from time domain to frequency domain. The accuracy and efficiency of 3D layer spectral element approach were verified by analyzing the Falling weight deflectometer(FWD) testing model with the spectral methods… More >

O. Hassan¹, K. Morgan, J. Jones, B. Larwood, N. P. Weatherill

CMES-Computer Modeling in Engineering & Sciences, Vol.5, No.5, pp. 383-394, 2004, DOI:10.3970/cmes.2004.005.383

Abstract A numerical procedure for the simulation of 3D problems involving the scattering of electromagnetic waves is presented. As practical problems of interest in this area often involve domains of complex geometrical shape, an unstructured mesh based method is adopted. The solution algorithm employs an explicit finite element procedure for the solution of Maxwell's curl equations in the time domain using unstructured tetrahedral meshes. A PML absorbing layer is added at the artificial far field boundary that is created by the truncation of the physical domain prior to the numerical solution. The complete solution procedure is parallelised and several large scale… More >

Y.C. Shiah¹, C. L. Tan², C.Y. Wang¹

CMES-Computer Modeling in Engineering & Sciences, Vol.87, No.1, pp. 1-22, 2012, DOI:10.3970/cmes.2012.087.001

Abstract The fundamental solution, or Green's function, for 3D anisotropic elastostatics as derived by Ting and Lee (1997) [Q.J. Mech. Appl. Math.; 50: 407-426] is one that is fully explicit and algebraic in form. It has, however, only been utilized in boundary element method (BEM) formulations quite recently even though it is relatively straightforward and direct to implement. This Green's function and its derivatives are necessary items in this numerical analysis technique. By virtue of the periodic nature of the angles when it is expressed in the spherical coordinate system, the present authors have very recently represented the Green's function as… More >

Yixiong Wei¹, Qifu Wang^1,2, Yingjun Wang¹, Yunbao Huang¹

CMES-Computer Modeling in Engineering & Sciences, Vol.86, No.3, pp. 241-274, 2012, DOI:10.3970/cmes.2012.086.241

Abstract In this study, we propose a novel method to accelerate the process of elastodynamic analysis in 3D problems with BEM (boundary element method). With applying the DRBEM (dual reciprocity boundary element method) to form new integral equations for reducing complexity;the modified FMM (fast multipole method)is introduced to simplify the computation process and save storage space by avoiding intermediate coefficientmatrices. At the same time, FMM-DRBEM is reprogrammed in parallel byapplying GPU with CUDA (Compute Unified Device Architecture)for improving efficiency further.The main features in this paper are: ( 1 )with respect to defects of classical method for elastodynamic, modified FMM-DRBEM algorithm is… More >

H. Çerdik Yaslan¹

CMES-Computer Modeling in Engineering & Sciences, Vol.86, No.1, pp. 29-38, 2012, DOI:10.3970/cmes.2012.086.029

Abstract This paper presents the approximate analytical solutions to the time fractional differential equations of elasticity for 3D quasicrystals with initial conditions. These equations are written in the form of a vector partial differential equation of the second order. The time fractional vector partial differential equations with initial conditions are solved by variational iteration method (VIM). The fractional derivatives are described in the Caputo sense. Numerical example shows that the proposed method is quite effective and convenient for solving kinds of time fractional system of partial differential equations. More >

Huimin Song¹, Dewey H. Hodges¹

CMES-Computer Modeling in Engineering & Sciences, Vol.85, No.3, pp. 239-278, 2012, DOI:10.3970/cmes.2012.085.239

Abstract The present paper presents a rigorous approach that can accurately and efficiently capture the linear, static and free-vibration behaviors of a beam-like structure by the rigorous combination of a one-dimensional beam model with a three-dimensional continuum model. This study focuses on coupling these disparate finite element types, putting them both into a single finite element model while making use of the asymptotically exact information available as part of the beam model, which itself is obtained by asymptotic dimensional reduction. The coupling is undertaken by use of appropriate transformation matrices at the interface together with stress and displacement recovery relations that… More >

Roman Chapko¹, B. Tomas Johansson²

CMES-Computer Modeling in Engineering & Sciences, Vol.85, No.2, pp. 105-128, 2012, DOI:10.3970/cmes.2012.085.105

Abstract We consider a Cauchy problem for the Laplace equation in a 3-dimen -sional semi-infinite domain that contains a bounded inclusion. The canonical situation is the upper half-space in I\tmspace -.1667em R3 containing a bounded smooth domain. The function value of the solution is specified throughout the plane bounding the upper half-space, and the normal derivative is given only on a finite portion of this plane. The aim is to reconstruct the solution on the surface of the bounded inclusion. This is a generalisation of the situation in Chapko and Johansson (2008) to three-dimensions and with Cauchy data only partially given.… More >

P. L. Bishay¹, S.N. Atluri¹

CMES-Computer Modeling in Engineering & Sciences, Vol.84, No.1, pp. 41-98, 2012, DOI:10.3970/cmes.2012.084.041

Abstract Higher-order two-dimensional as well as low and higher-order three-dimensional new Hybrid/Mixed (H/M) finite elements based on independently assumed displacement, and judiciously chosen strain fields, denoted by HMFEM-2, are developed here for applications in macro-mechanics. The idea of these new H/M finite elements is based on collocating the components of the independent strain field, with those derived from the independently assumed displacement fields at judiciously and cleverly chosen collocation points inside the element. This is unlike the other techniques used in older H/M finite elements where a two-field variational principle was used in order to enforce both equilibrium and compatibility conditions… More >