A. Rama Mohan Rao¹, T.V.S.R. Appa Rao², B. Dattaguru³

CMES-Computer Modeling in Engineering & Sciences, Vol.5, No.3, pp. 213-234, 2004, DOI:10.3970/cmes.2004.005.213

Abstract This paper presents an algorithm for automatic partitioning of unstructured meshes for parallel finite element computations employing float-encoded genetic algorithms (FEGA). The problem of mesh partitioning is represented in such a way that the number of variables considered in the genome (chromosome) construction is constant irrespective of the size of the problem. In order to accelerate the computational process, several acceleration techniques like constraining the search space, local improvement after initial global partitioning have been attempted. Finally, micro float-encoded genetic algorithms have been developed to accelerate the computational process. More >

Z. D. Han¹, S. N. Atluri²

CMC-Computers, Materials & Continua, Vol.1, No.2, pp. 129-140, 2004, DOI:10.3970/cmc.2004.001.129

Abstract A Meshless Local Petrov-Galerkin (MLPG) method has been developed for solving 3D elasto-dynamic problems. It is derived from the local weak form of the equilibrium equations by using the general MLPG concept. By incorporating the moving least squares (MLS) approximations for trial and test functions, the local weak form is discretized, and is integrated over the local sub-domain for the transient structural analysis. The present numerical technique imposes a correction to the accelerations, to enforce the kinematic boundary conditions in the MLS approximation, while using an explicit time-integration algorithm. Numerical examples for solving the transient More >

B. Amaziane¹, A. Naji², D. Ouazar³

CMC-Computers, Materials & Continua, Vol.1, No.2, pp. 117-128, 2004, DOI:10.3970/cmc.2004.001.117

Abstract In this paper, a meshless method based on Radial Basis Functions (RBF) is coupled with genetic algorithms for parameter identification to some selected groundwater flow applications. The treated examples are generated by the diffusion equation with some specific boundary conditions describing the groundwater fluctuation in a leaky confined aquifer system near open tidal water. To select the best radial function interpolation and show the powerful of the method in comparison to domain based discretization methods, Multiquadric (MQ), Thin-Plate Spline (TPS) and Conical type functions are investigated and compared to finite difference results or analytical one. 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 More >

J. Sorić¹, Q. Li², T. Jarak¹, S.N. Atluri²

CMES-Computer Modeling in Engineering & Sciences, Vol.6, No.4, pp. 349-358, 2004, DOI:10.3970/cmes.2004.006.349

Abstract An efficient meshless formulation based on the Local Petrov-Galerkin approach for the analysis of shear deformable thick plates is presented. Using the kinematics of a three-dimensional continuum, the local symmetric weak form of the equilibrium equations over the cylindrical shaped local sub-domain is derived. The linear test function in the plate thickness direction is assumed. Discretization in the in-plane directions is performed by means of the moving least squares approximation. The linear interpolation over the thickness is used for the in-plane displacements, while the hierarchical quadratic interpolation is adopted for the transversal displacement in order More >

J. Sladek¹, V. Sladek¹, S.N. Atluri²

CMES-Computer Modeling in Engineering & Sciences, Vol.6, No.3, pp. 309-318, 2004, DOI:10.3970/cmes.2004.006.309

Abstract Meshless methods based on the local Petrov-Galerkin approach are proposed for solution of steady and transient heat conduction problem in a continuously nonhomogeneous anisotropic medium. Fundamental solution of the governing partial differential equations and the Heaviside step function are used as the test functions in the local weak form. It is leading to derive local boundary integral equations which are given in the Laplace transform domain. The analyzed domain is covered by small subdomains with a simple geometry. To eliminate the number of unknowns on artificial boundaries of subdomains the modified fundamental solution and/or the More >

B. Jin^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.6, No.3, pp. 253-262, 2004, DOI:10.3970/cmes.2004.006.253

Abstract In this paper, we propose a new numerical scheme for the solution of the Laplace and biharmonic equations subjected to noisy boundary data. The equations are discretized by the method of fundamental solutions. Since the resulting matrix equation is highly ill-conditioned, a regularized solution is obtained using the truncated singular value decomposition, with the regularization parameter given by the L-curve method. Numerical experiments show that the method is stable with respect to the noise in the data, highly accurate and computationally very efficient. More >

E. Ferretti¹

CMES-Computer Modeling in Engineering & Sciences, Vol.6, No.3, pp. 227-244, 2004, DOI:10.3970/cmes.2004.006.227

Abstract Defining the crack path in brittle and non-brittle crack is not easy, due to several unknowns. If the direction of crack propagation can be computed by means of one of the existing criteria, it is not known whether this direction will remain constant during crack propagation. A crack initiation leads to an enhanced stress field at crack tip. During propagation, the enhanced tip stress field propagates into the solid, locally interacting with the pre-existing stress field. This interaction can lead to modifications of the propagation direction, depending on the domain and crack geometry. Moreover, trajectory… More >

Soon Wan Chung¹, Yoo Jin Choi¹, Seung Jo Kim¹

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.6, pp. 707-718, 2003, DOI:10.3970/cmes.2003.004.707

Abstract In this paper, an automatic finite element remeshing algorithm based on the bubble packing method is utilized for the purpose of numerical simulations of localized damage, because fine meshes are needed to represent the gradually concentrated damage. The bubble packing method introduces two parameters that easily control the remeshing criterion and the new mesh size. The refined area is determined by \textit {a posteriori} error estimation utilizing the value obtained from Superconvergent Patch Recovery. The isotropic ductile damage theory, founded on continuum damage mechanics, is used for this damage analysis. It was successfully shown in More >

N. Sukumar¹, J. E. Bolander¹

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.6, pp. 691-706, 2003, DOI:10.3970/cmes.2003.004.691

Abstract In this paper, we explore the numerical approximation of discrete differential operators on non-uniform grids. The Voronoi cell and the notion of natural neighbors are used to approximate the Laplacian and the gradient operator on irregular grids. The underlying weight measure used in the numerical computations is the {\em Laplace weight function}, which has been previously adopted in meshless Galerkin methods. We develop a difference approximation for the diffusion operator on irregular grids, and present numerical solutions for the Poisson equation. On regular grids, the discrete Laplacian is shown to reduce to the classical finite More >