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 >

Tong-Yi Zhang¹

CMC-Computers, Materials & Continua, Vol.1, No.1, pp. 107-116, 2004, DOI:10.3970/cmc.2004.001.107

Abstract The present work presents a strip Dielectric Breakdown (DB) model for an electrically impermeable crack in a piezoelectric material. In the DB model, the dielectric breakdown region is assumed to be a strip along the crack's front line. Along the DB strip, the electric field strength is equal to the dielectric breakdown strength. The DB model is exactly in analogy with the mechanical Dugdale model. Two energy release rates emerge from the analysis. An applied energy release rate appears when evaluating J-integral along a contour surrounding both the dielectric breakdown strip and the crack tip, whereas More >

X.S. Sun¹, L.X. Huang¹, Y.H. Liu¹, Z.Z. Cen^1,2

CMC-Computers, Materials & Continua, Vol.1, No.1, pp. 91-106, 2004, DOI:10.3970/cmc.2004.001.091

Abstract The Boundary Element Method (BEM) is introduced to analyze the elasto-plastic problems of 2-D orthotropic bodies. With the help of known boundary integral equations and fundamental solutions, a numerical scheme for elasto-plastic analysis of 2-D orthotropic problems with the BEM is developed. The Hill orthotropic yield criterion is adopted in the plastic analysis. The initial stress method and tangent predictor-radial return algorithm are used to determine the stress state in solving the nonlinear equation with the incremental iteration method. Finally, numerical examples show that the BEM is effective and reliable in analyzing elasto-plastic problems of More >

S.-C. Hsiao¹, M.S. Ingber²

CMC-Computers, Materials & Continua, Vol.1, No.1, pp. 51-58, 2004, DOI:10.3970/cmc.2004.001.051

Abstract The lateral migration of nonneutrally-buoyant spherical particles in Poiseuille flow is investigated numerically using the boundary element method. In particular, the steady, Navier-Stokes equations are solved using a classical domain integration method treating the nonlinear terms as pseudo-body forces. The numerical results for the lateral migration velocity are compared with experimental data. The numerical results indicate that the lateral migration velocity does not scale linearly with the Reynolds number. The methodology is extended to include non-Newtonian power-law fluids. The migration velocity is significantly affected for particles suspended in this class of fluids and can actually More >

X. Han^1,2,3, G. R. Liu^1,2, G. Y. Li ¹

CMC-Computers, Materials & Continua, Vol.1, No.1, pp. 39-50, 2004, DOI:10.3970/cmc.2004.001.039

Abstract An efficient hybrid numerical method is presented for investigating transient response of cross-ply laminated axisymmetric cylinders subjected to an impact load. In this hybrid numerical method, the laminated cylinder is divided into layered cylindrical elements in the thickness direction. The Hamilton principle is used to develop governing equations of the structure. The displacement response is determined by employing the Fourier transformations and the modal analysis. Numerical examples for analyzing transient waves have been provided in axisymmetric laminated cylindrical structures, both for thin cylindrical shells and thick cylinders.
A computational inverse technique is also presented for More >

Hokwon A. Cho¹, M. A. Golberg², A. S. Muleshkov¹, Xin Li¹

CMC-Computers, Materials & Continua, Vol.1, No.1, pp. 1-38, 2004, DOI:10.3970/cmc.2004.001.001

Abstract In this paper we present a mesh-free approach to numerically solving a class of second order time dependent partial differential equations which include equations of parabolic, hyperbolic and parabolic-hyperbolic types. For numerical purposes, a variety of transformations is used to convert these equations to standard reaction-diffusion and wave equation forms. To solve initial boundary value problems for these equations, the time dependence is removed by either the Laplace or the Laguerre transform or time differencing, which converts the problem into one of solving a sequence of boundary value problems for inhomogeneous modified Helmholtz equations. These… More >

John L. Volakis¹, Kubilay Sertel¹, Erik Jørgensen², Rick W. Kindt¹

CMES-Computer Modeling in Engineering & Sciences, Vol.5, No.5, pp. 463-476, 2004, DOI:10.3970/cmes.2004.005.463

Abstract In this paper we address several topics relating to the development and implementation of volume integral and hybrid finite element methods for electromagnetic modeling. Comparisons of volume integral equation formulations with the finite element-boundary integral method are given in terms of accuracy and computing resources. We also discuss preconditioning and parallelization of the multilevel fast multipole method, and propose higher-order basis functions for curvilinear quadrilaterals and volumetric basis functions for curvilinear hexahedra. The latter have the desirable property of vanishing divergence within the element but non-zero curl. In addition, a new domain decomposition is introduced More >

C. J. Reddy¹

CMES-Computer Modeling in Engineering & Sciences, Vol.5, No.5, pp. 455-462, 2004, DOI:10.3970/cmes.2004.005.455

Abstract This paper presents an efficient algorithm to evaluate multi-spectral and multi-angular monostatic radar cross section (RCS) of large objects with very fine increments. The technique is based on the combination of Model Based Parameter Estimation (MBPE) method with hybrid frequency domain formulations. A general approach to formulation of MBPE is presented along with a similar approach called the Asymptotic Waveform Evaluation (AWE). Various numerical examples are presented for multi-spectral response calculations using method of moments (MoM) and the hybrid Finite Element-MoM technique in conjunction with MBPE. Example application of MBPE for hybrid MoM-Physical Optics approach More >

Raju R. Namburu¹, Eric R. Mark, Jerry A. Clarke

CMES-Computer Modeling in Engineering & Sciences, Vol.5, No.5, pp. 443-454, 2004, DOI:10.3970/cmes.2004.005.443

Abstract Computational electromagnetic (CEM) simulations of full-range military vehicles play a critical role in enhancing the survivability and target recognition of combat systems. Modeling of full-range military systems subjected to high frequencies may involve generating large-scale meshes, solving equations, visualization, and analysis of results in the range of billions of unknowns or grid points. Hence, the overall objective of this research is to develop and demonstrate a scalable CEM software environment to address accurate prediction of radar cross sections (RCS) for full- range armored vehicles with realistic material treatments and complex geometric configurations. A software environment… More >