D.S. Liu^1,2, Z.H. Fong¹, I.H. Lin¹, Z.W. Zhuang¹

CMES-Computer Modeling in Engineering & Sciences, Vol.93, No.6, pp. 441-468, 2013, DOI:10.3970/cmes.2013.093.441

Abstract In this study, we proposed a novel numerical technique to simulate the transient moisture diffusion process and to apply it to heterogeneous composite resins. The method is based on a heterogeneous hybrid moisture element (HHME), with properties determined through an equivalent hybrid moisture capacitance/ conductance matrix that was calculated using the conventional finite element formulation in space discretization and the q-method in time discretization, with similar mass/stiffness properties and matrix condensing operations. A coupled HHME with finite element scheme was developed and implemented in the computer code by using the commercial software MATLAB to analyze the transient moisture diffusion process… More >

Ji-Chuan Liu¹, Quan-Guo Zhang²

CMES-Computer Modeling in Engineering & Sciences, Vol.93, No.3, pp. 203-220, 2013, DOI:10.3970/cmes.2013.093.203

Abstract In this paper, we propose an algorithm to solve a Cauchy problem of the Laplace equation in doubly connected domains for 2D and 3D cases in which the Cauchy data are given on the outer boundary. We want to seek a solution in the form of the single-layer potential and discrete it by parametrization to yield an ill-conditioned system of algebraic equations. Then we apply the Tikhonov regularization method to solve this ill-posed problem and obtain a stable numerical solution. Based on the regularization parameter chosen suitably by GCV criterion, the proposed method can get the approximate temperature and heat… More >

Leiting Dong^1,2, Satya N. Atluri^1,3

CMES-Computer Modeling in Engineering & Sciences, Vol.90, No.5, pp. 379-413, 2013, DOI:10.3970/cmes.2013.090.379

Abstract The SGBEM-FEM alternating method is compared with the recently popularized XFEM, for analyzing mixed-mode fracture and fatigue growth of 3D nonplanar cracks in complex solid and structural geometries. A large set of 3D examples with different degrees of complexity is analyzed by the SGBEM-FEM alternating method, and the numerical results are compared with those obtained by XFEM available in the open literature. It is clearly shown that: (a) SGBEM-FEM alternating method gives extremely high accuracy for the stress intensity factors; but the XFEM gives rather poor computational results, even for the most simple 3D cracks; (b) while SGBEM-FEM alternating method… More >

Zhaohui Xu¹, Deli Gao¹

CMES-Computer Modeling in Engineering & Sciences, Vol.89, No.1, pp. 17-24, 2012, DOI:10.3970/cmes.2012.089.017

Abstract Little is explained about the process of casing window milling for sidetracking due to lack of analytical method for its mechanical characteristic. In this paper, 3D FE models are established using the commercial finite-element software ABAQUS/Explicit to make simulation analysis for two key stages of the process including the initial stage of casing milling and the stage of full-gauge casing window milling. The models involve the effects of main drilling parameters such as reaction force, torque, speed, feed rate per revolution, and milling angle. The calculation results verify the capability and advantages of 3D FE simulation for the process of… More >

Lior Banai

CMES-Computer Modeling in Engineering & Sciences, Vol.88, No.5, pp. 351-366, 2012, DOI:10.3970/cmes.2012.088.351

Abstract The explosive Engineering field is a costly one in which not every organization can effort the time and money it takes to performed field tests on its explosives. The purpose of this article is to present a program that was developed in the Israeli Navy for performance estimations and safety issues of warheads and explosives. With a relative small developing time one can create a tool that gives preliminary results in a few minutes without the need to design and order a field tests or run finite elements analyses. By implementing a few known models, in this tool, the user… More >

Annamaria Mazzia¹, Giorgio Pini¹, Flavio Sartoretto²

CMES-Computer Modeling in Engineering & Sciences, Vol.88, No.3, pp. 183-210, 2012, DOI:10.3970/cmes.2012.088.183

Abstract Pure meshless techniques are promising methods for solving Partial Differential Equations (PDE). They alleviate difficulties both in designing discretization meshes, and in refining/coarsening, a task which is demanded e.g. in adaptive strategies. Meshless Local Petrov Galerkin (MLPG) methods are pure meshless techniques that receive increasing attention. Very recently, new methods, called Direct MLPG (DMLPG), have been proposed. They rely upon approximating PDE via the Generalized Moving Least Square method. DMLPG methods alleviate some difficulties of MLPG, e.g. numerical integration of tricky, non-polynomial factors, in weak forms. DMLPG techniques require lower computational costs respect to their MLPG counterparts. In this paper… More >

Qiao Wang¹, Yu Miao^1,2, Junjie Zheng¹

CMES-Computer Modeling in Engineering & Sciences, Vol.70, No.2, pp. 123-152, 2010, DOI:10.3970/cmes.2010.070.123

Abstract In this paper, a fast formulation of the hybrid boundary node method (Hybrid BNM) for solving 3D elasticity is presented. Coupling modified variational principle with the Moving Least Squares (MLS) approximation, the Hybrid BNM only requires discrete nodes constructed on the surface of a domain. The preconditioned GMERS is employed to solve the resulting system of equations. At each iteration step of the GMERS, the matrix-vector multiplication is accelerated by the fast multipole method (FMM). The fundamental solution of three-dimensional elasticity problem is expanded in terms of series. An oct-tree data structure is adopted to subdivide the computational domain into… More >

Y.C. Shiah¹, C. L. Tan², R.F. Lee¹

CMES-Computer Modeling in Engineering & Sciences, Vol.69, No.2, pp. 167-198, 2010, DOI:10.3970/cmes.2010.069.167

Abstract In this paper, fully explicit, algebraic expressions are derived for the first and second derivatives of the Green's function for the displacements in a three dimensional anisotropic, linear elastic body. These quantities are required in the direct formulation of the boundary element method (BEM) for determining the stresses at internal points in the body. To the authors' knowledge, similar quantities have never previously been presented in the literature because of their mathematical complexity. Although the BEM is a boundary solution numerical technique, solutions for the displacements and stresses at internal points are sometimes required for some engineering applications. To this… More >

Giorgio Pini¹, Annamaria Mazzia¹, Flavio Sartoretto²,

CMES-Computer Modeling in Engineering & Sciences, Vol.36, No.1, pp. 43-64, 2008, DOI:10.3970/cmes.2008.036.043

Abstract Meshless methods have been explored in many 2D problems and they have been shown to be as accurate as Finite Element Methods (FEM). Compared to the extensive literature on 2D applications, papers on solving 3D problems by meshless methods are surprisingly few. Indeed, a main drawback of these methods is the requirement for accurate cubature rules. This paper focuses on the so called Meshless Local Petrov Galerkin (MLPG) methods. We show that accurate solutions of 3D potential problems can be attained, provided suitable cubature rules are identified, sparse data structures are efficiently stored, and strategies are devised in order to… More >

Xuefei He¹, Kian-Meng Lim^1,2,3, Siak-Piang Lim^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.35, No.1, pp. 21-48, 2008, DOI:10.3970/cmes.2008.035.021

Abstract The boundary element method (BEM) is known to have the advantage of reducing the dimension of problem by discretizing only the boundary of the domain. But it becomes less attractive for solving Poisson-type equations, due to the need to evaluate the domain integral which is computationally expensive. In this paper, we present the extension of a recently developed fast algorithm for Laplace equation, based on fast Fourier transform on multipoles (FFTM), to solve large scale 3D Poisson-type equations. We combined the Laplace solver with two fast methods for handling the domain integral based on fast Fourier transform (FFT). The first… More >