E. Ferretti¹, E. Casadio, A. Di Leo¹

CMES-Computer Modeling in Engineering & Sciences, Vol.30, No.3, pp. 163-190, 2008, DOI:10.3970/cmes.2008.030.163

Abstract In this study, the Cell Method (CM) is applied in order to investigate the failure mechanisms of masonry walls under shear force. The direction of propagation is computed step-wise by the code, and the domain is updated by means of a propagation technique of intra-element nodal relaxation with re-meshing. The crack extension condition is studied in the Mohr/Coulomb plane, using the criterion of Leon. The main advantage of using the CM for numerical analyses of masonry is that the mortar, the bricks and the interfaces between mortar and bricks can be modeled without any need… More >

Dong Ju Woo¹, Jin Oh Yang¹, Beom-Soo Kim¹, Seungsoo Lee¹, Jin Yeon Cho²

CMES-Computer Modeling in Engineering & Sciences, Vol.30, No.3, pp. 149-162, 2008, DOI:10.3970/cmes.2008.030.149

Abstract A new orphan-cell-free overset method is proposed to carry out the coupled analysis of overlapping finite element substructures. In the proposed overset method, the meshless MLS (Moving Least Squares) approximation is used to obtain the boundary data for the overlapped interface, whereas the Lagrange interpolation scheme has been commonly used in the conventional overset methods. The meshless character of MLS approximation makes it possible to overcome the problem of orphan-cell, which is often encountered in the conventional overset methods. Further, a new connectivity matrix solution procedure is developed to reduce the computational time in the More >

P.H. Wen¹, M.H. Aliabadi², Y.W. Liu³

CMES-Computer Modeling in Engineering & Sciences, Vol.30, No.3, pp. 133-148, 2008, DOI:10.3970/cmes.2008.030.133

Abstract Based on the variation of potential energy, the element-free Galerkin method (MFGM) has been investigated for structures with crack on the basis of radial base function interpolation. An enriched radial base function is introduced to capture the singularities of stress at the crack tips. The advantages of the finite element method are remained in this method and there is a significant improvement of accuracy, particularly for the crack problems of fracture mechanics. The applications of the element-free Galerkin method with enriched radial base function to two-dimensional fracture mechanics in functionally graded materials have been presented 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). More >

Liviu Marin¹

CMES-Computer Modeling in Engineering & Sciences, Vol.30, No.2, pp. 99-122, 2008, DOI:10.3970/cmes.2008.030.099

Abstract The application of the method of fundamental solutions (MFS) to inverse boundary value problems associated with the steady-state heat conduction in isotropic media in the presence of sources, i.e. the Poisson equation, is investigated in this paper. Based on the approach of Alves and Chen (2005), these problems are solved in two steps, namely by finding first an approximate particular solution of the Poisson equation and then the numerical solution of the resulting inverse boundary value problem for the Laplace equation. The resulting MFS discretised system of equations is ill-conditioned and hence it is solved More >

J. Sladek¹, V. Sladek¹, P. Solek², P.H. Wen³, S.N. Atluri⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.30, No.2, pp. 77-98, 2008, DOI:10.3970/cmes.2008.030.077

Abstract A meshless local Petrov-Galerkin (MLPG) method is applied to solve problems of Reissner-Mindlin shells under thermal loading. Both stationary and thermal shock loads are analyzed here. Functionally graded materials with a continuous variation of properties in the shell thickness direction are considered here. A weak formulation for the set of governing equations in the Reissner-Mindlin theory is transformed into local integral equations on local subdomains in the base plane of the shell by using a unit test function. Nodal points are randomly spread on the surface of the plate and each node is surrounded by More >

Chein-Shan Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.30, No.2, pp. 65-76, 2008, DOI:10.3970/cmes.2008.030.065

Abstract Trefftz method (TM) is one of widely used meshless numerical methods in elliptic type boundary value problems, of which the approximate solution is expressed as a linear combination of T-complete bases, and the unknown coefficients are determined from boundary conditions by solving a linear equations system. However, the accuracy of TM is severely limited by its ill-conditioning. This paper is a continuation of the work of Liu (2007a). The collocation TM is modified and applied to the direct and inverse problems of biharmonic equation in a simply connected plane domain. Due to its well-conditioning of More >

Shuling Hu¹, Shengping Shen^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.30, No.2, pp. 57-64, 2008, DOI:10.3970/cmes.2008.030.057

Abstract The effects of surface energy on the interaction between two voids of equal size are investigated. The problem is solved by series expansion in bipolar coordinates. The results show that the surface energy significantly affects the stress concentration around the holes as the size of the holes shrinks to nanometers, meanwhile the interaction between the holes also influences the stress distribution around the holes, which become evident as the holes close to each other. This problem is of great importance in engineering applications. More >

Shu Li¹, S. N. Atluri²

CMES-Computer Modeling in Engineering & Sciences, Vol.30, No.1, pp. 37-56, 2008, DOI:10.3970/cmes.2008.030.037

Abstract In this paper, a method based on a combination of an optimization of directions of orthotropy, along with topology optimization, is applied to continuum orthotropic solids with the objective of minimizing their compliance. The spatial discretization algorithm is the so called Meshless Local Petrov-Galerkin (MLPG) "mixed collocation'' method for the design domain, and the material-orthotropy orientation angles and the nodal volume fractions are used as the design variables in material optimization and topology optimization, respectively. Filtering after each iteration diminishes the checkerboard effect in the topology optimization problem. The example results are provided to illustrate More >

Sen Yung Lee¹, Sheei Muh Lin², Chien Shien Lee³, Shin Yi Lu³, Yen Tse Liu³

CMES-Computer Modeling in Engineering & Sciences, Vol.30, No.1, pp. 27-36, 2008, DOI:10.3970/cmes.2008.030.027

Abstract An analytic solution method, namely the shifting function method, is developed to find the exact large static deflection of a beam with nonlinear elastic springs supports at ends for the first time. The associated mathematic system is a fourth order ordinary differential equation with nonlinear boundary conditions. It is shifted and decomposed into five linear differential equations and at most four algebra equations. After finding the roots of the algebra equations, the exact solution of the nonlinear beam system can be reconstructed. It is shown that the proposed method is valid for the problem with More >