Zhen Luo^1,2, Nong Zhang^1,3, Yu Wang¹

CMES-Computer Modeling in Engineering & Sciences, Vol.83, No.1, pp. 73-96, 2012, DOI:10.3970/cmes.2012.083.073

Abstract This paper aims to present a physically meaningful level set method for shape and topology optimization of structures. Compared to the conventional level set method which represents the design boundary as the zero level set, in this study the boundary is embedded into non-zero constant level sets of the level set function, to implicitly implement shape fidelity and topology changes in time via the propagation of the discrete level set function. A point-wise nodal density field, non-negative and value-bounded, is used to parameterize the level set function via the compactly supported radial basis functions (CSRBFs) at a uniformly defined set… More >

C. A. Bustamante¹, W. F. Florez¹, H. Power², M. Giraldo¹, A. F. Hill¹

CMES-Computer Modeling in Engineering & Sciences, Vol.79, No.2, pp. 103-130, 2011, DOI:10.3970/cmes.2011.079.103

Abstract The two-dimensional Navier Stokes system of equations for incompressible flows is solved in the velocity vorticity formulation by means of the Control Volume-Radial Basis Function (CV-RBF) method. This method is an improvement to the Control Volume Method (CVM) based on the use of Radial Basis Function (RBF) Hermite interpolation instead of the classical polynomial functions. The main advantages of the CV-RBF method are the approximation order, the meshless nature of the interpolation scheme and the presence of the PDE operator in the interpolation. Besides, the vorticity boundary values are computed in terms of the values of the velocity field at… More >

K.I. Silva¹, J.C.F. Telles², F.C. Araújo³

CMES-Computer Modeling in Engineering & Sciences, Vol.78, No.3&4, pp. 209-224, 2011, DOI:10.3970/cmes.2011.078.209

Abstract In this work the boundary element method is applied to solve 2D elastoplastic problems. In elastoplastic boundary element analysis, domain integrals have to be calculated to introduce the contribution of yielded zones. Traditionally, the use of internal integration cells have been adopted to evaluate such domain integrals. The present work, however, proposes an alternative cell free strategy based on the OMLS (Orthogonal Moving Least Squares) interpolation, typically adopted in meshless methods. In this approach the definition of points to compute the interpolated value of a function at a given location only depends on their relative distance, without need to define… More >

A. Khosravifard¹, M.R. Hematiyan^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.78, No.3&4, pp. 185-208, 2011, DOI:10.3970/cmes.2011.078.185

Abstract An inverse method for optimal control of the freezing front motion in the solidification of pure materials is presented. The inverse technique utilizes the idea of a pseudo heat source to account for the latent heat effects. The numerical formulation of this inverse method is based on a formerly introduced meshless technique. In this method, the flux and the velocity of the liquid-solid interface are treated as secondary variables and the liquid and solid domains are modeled simultaneously. Some numerical examples are provided to demonstrate the efficiency of the presented method. The effects of regularization and the number of nodes… More >

K. Wang¹, S.J. Zhou^2,3, Z.F. Nie⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.73, No.1, pp. 77-102, 2011, DOI:10.3970/cmes.2011.073.077

Abstract The natural neighbour Galerkin method is tailored to solve boundary value problems of the couple-stress elasticity to model the size dependent behaviour of materials. This method is based on the displacement-based Galerkin approach, and the calculation of the global stiffness matrix is performed using gradient smoothing technique combined with the non-Sibsonian partition of unity approximation scheme. This method possesses the following properties: the complex C1-continuous approximation scheme is avoided without using either Lagrange multipliers or penalty parameters; no domain integrals involved in the assembly of the global stiffness matrix; and the imposition of essential boundary conditions is straightforward. The validity… More >

Htike Htike¹, Wen Chen¹, Yan Gu¹, Junjie Yang¹

CMES-Computer Modeling in Engineering & Sciences, Vol.72, No.4, pp. 247-272, 2011, DOI:10.3970/cmes.2011.072.247

Abstract This paper makes the first attempt to employ the Radial Basis Function (RBF) interpolation in the material point method (MPM), in which the shape function is based on RBF and polynomial function and satisfies the partition of unity and possesses Delta-function property. It is worthy of stressing that the RBF interpolation has the merit of high smoothness and is very accurate and can easily be applied to the MPM framework for mapping information between moving particles, known as material point in the MPM, and background grids. The RBF-based MPM is designed to overcome the unphysical results, such as shear stress… More >

J. Sladek¹, V. Sladek¹, Ch. Zhang², M. Wünsche²

CMES-Computer Modeling in Engineering & Sciences, Vol.68, No.2, pp. 185-220, 2010, DOI:10.3970/cmes.2010.068.185

Abstract A meshless method based on the local Petrov-Galerkin approach is proposed to solve initial-boundary value crack problems of piezoelectric solids with nonlinear electrical boundary conditions on crack faces. Homogeneous and continuously varying material properties of the piezoelectric solid are considered. Stationary governing equations for electrical fields and the elastodynamic equations with an inertial term for mechanical 2-D fields are considered. Nodal points are spread on the analyzed domain, and each node is surrounded by a small circle for simplicity. The spatial variation of displacements and electric potential are approximated by the Moving Least-Squares (MLS) scheme. After performing the spatial integrations,… More >

Hua Zou¹, Hua Li^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.67, No.1, pp. 55-78, 2010, DOI:10.3970/cmes.2010.067.055

Abstract A novel meshless technique termed the Random Integral Quadrature (RIQ) method is developed in this paper for solving the generalized integral equations. By the RIQ method, the governing equations in the integral form are discretized directly with the field nodes distributed randomly or uniformly, which is achieved by discretizing the integral governing equations with the generalized integral quadrature (GIQ) technique over a set of background virtual nodes, and then interpolating the function values at the virtual nodes over a set of field nodes with Local Kriging method, where the field nodes are distributed either randomly or uniformly. The RIQ method… More >

Roman Chapko¹, B. Tomas Johansson², Oleh Sobeyko¹

CMES-Computer Modeling in Engineering & Sciences, Vol.62, No.1, pp. 57-76, 2010, DOI:10.3970/cmes.2010.062.057

Abstract We propose an iterative procedure for the inverse problem of determining the displacement vector on the boundary of a bounded planar inclusion given the displacement and stress fields on an infinite (planar) line-segment. At each iteration step mixed boundary value problems in an elastostatic half-plane containing the bounded inclusion are solved. For efficient numerical implementation of the procedure these mixed problems are reduced to integral equations over the bounded inclusion. Well-posedness and numerical solution of these boundary integral equations are presented, and a proof of convergence of the procedure for the inverse problem to the original solution is given. Numerical… More >

S.L. Li^1,2, S.Y. Long¹, G.Y. Li¹

CMES-Computer Modeling in Engineering & Sciences, Vol.60, No.1, pp. 73-94, 2010, DOI:10.3970/cmes.2010.060.073

Abstract A new implementation of topology optimization for the plate described by the Reissner-Mindlin theory based on the meshless natural neighbour Petrov-Galerkin method (NNPG) is proposed in this work. The objective is to produce the stiffest plate for a given volume by redistributing the material throughout the plate. We try to couple the advantages of the meshless numerical method with the topology optimization of moderately thick plate. The numerical approach presented here is based on the solid isotropic material with penalization (SIMP) formulation of the topology optimization problem. The natural neighbour interpolation shape function is employed to discretize both displacement and… More >