T.C. Theodosiou¹, D.A. Saravanos²

CMES-Computer Modeling in Engineering & Sciences, Vol.19, No.2, pp. 121-134, 2007, DOI:10.3970/cmes.2007.019.121

Abstract In this paper a new method is introduced for carbon nanotube modeling combining features of Molecular Mechanics and Finite Element Analysis. Repetitive atomic cells are treated as finite elements, whose internal energy is determined by the semi-empirical Brenner molecular potential model; internal forces and linearized stiffness matrices are formulated analytically in order to gain in speed and accuracy, and the resultant discrete system is formulated and solved using the Newton-Raphson method. The presented method is validated through comparisons to numerical and experimental results provided by other researchers. The bending and shearing of CNTs is also More >

Shunde Yin¹, Leo Rothenburg¹, Maurice B. Dusseault¹

CMES-Computer Modeling in Engineering & Sciences, Vol.19, No.2, pp. 111-120, 2007, DOI:10.3970/cmes.2007.019.111

Abstract Subsidence problem is of great importance in petroleum engineering and environmental engineering. In this paper, we firstly apply a hybrid Displacement Discontinuity-FEM modeling to this classic problem: the evaluation of subsidence over a compacting oil reservoir. We use displacement discontinuity method to account for the reservoir surrounding area, and finite element methods in the fully coupled simulation of the reservoir itself. This approach greatly reduces the number of degrees of freedom compared to an analyzing fully coupled problem using only a finite element or finite difference discretization. More >

Chein-Shan Liu ¹

CMES-Computer Modeling in Engineering & Sciences, Vol.19, No.1, pp. 99-110, 2007, DOI:10.3970/cmes.2007.019.099

Abstract A new meshless regularized integral equation method (MRIEM) is developed to solve the interior and exterior Dirichlet problems for the two-dimensional Laplace equation, which consists of three parts: Fourier series expansion, the second kind Fredholm integral equation and an analytically regularized solution of the unknown boundary condition on an artificial circle. We find that the new method is powerful even for the problem with complex boundary shape and with random noise disturbing the boundary data. More >

Amit Shaw¹, D Roy²

CMES-Computer Modeling in Engineering & Sciences, Vol.19, No.1, pp. 69-98, 2007, DOI:10.3970/cmes.2007.019.069

Abstract A novel discretization strategy and derivative reproduction based on reproducing kernel (RK) particle approximations of functions are proposed. The proposed scheme is in the form of an RK interpolation that offers significant numerical advantages over a recent version of the strategy by Chen et al. (2003), wherein the authors added a set of primitive functions to the reproducing kernel (enrichment) functions. It was also required that the support size of the primitive function be less than the smallest distance between two successive grid points. Since the primitive function was required to vary from 0 to… More >

J. Perko¹, B. Šarler²

CMES-Computer Modeling in Engineering & Sciences, Vol.19, No.1, pp. 55-68, 2007, DOI:10.3970/cmes.2007.019.055

Abstract This work introduces a procedure for automated determination of weight function free parameters in moving least squares (MLS) based meshless methods for non-uniform grids. The meshless method used in present work is Diffuse Approximate Method (DAM). The DAM is structured in 2D with the one or two parameter Gaussian weigh function, 6 polynomial basis and 9 noded domain of influence. The procedure consists of three main elements. The first is definition of the reference quality function which measures the difference between the MLS approximation on non-uniform and hypothetic uniform node arrangements. The second is the… More >

M. Manjuprasad¹, C. S. Manohar²

CMES-Computer Modeling in Engineering & Sciences, Vol.19, No.1, pp. 23-54, 2007, DOI:10.3970/cmes.2007.019.023

Abstract A technique for adaptive random field refinement for stochastic finite element reliability analysis of structures is presented in this paper. Refinement indicator based on global importance measures are proposed and used for carrying out adaptive random field mesh refinements. Reliability index based error indicator is proposed and used for assessing the percentage error in the estimation of notional failure probability. Adaptive mesh refinement is carried out using hierarchical graded mesh obtained through bisection of elements. Spatially varying stochastic system parameters (such as Young's modulus and mass density) and load parameters are modeled in general as… More >

F. Bierbrauer, S.-P. Zhu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.19, No.1, pp. 1-22, 2007, DOI:10.3970/cmes.2007.019.001

Abstract Recently the use of the one-field formulation in the numerical solution of the Navier-Stokes equations for two-phase incompressible flow has become a very attractive approach in CFD (computational fluid dynamics). While the presence of material discontinuities across fluid interfaces presents some difficulty, it is their combination with a non-solenoidal discontinuous initial velocity field, commonly occurring in the mathematical formulation, that has provided the greatest hindrance in the numerical solution. This paper presents three analytical solutions, the Bounded Creeping Flow, Solenoidal and Conserved Solenoidal Solutions, which are both continuous, incompressible, retain as much of the original More >

C. P. Providakis¹

CMES-Computer Modeling in Engineering & Sciences, Vol.18, No.3, pp. 247-258, 2007, DOI:10.3970/cmes.2007.018.247

Abstract The method of Domain/Boundary Element is used to achieve a dynamic analysis of elastoplastic thick plates resting on internal supports. All possible boundary conditions on the edge of the plate with any interior support conditions such as isolated points (column), lines (walls) or regions (patches) can be treated without practical difficulties. The formulation presented includes the effects of shear deformation and rotatory inertia following Reissner-Mindlin's deformation theory assumptions. The method employs the elastostatic fundamental solution of the problem resulting in both boundary and domain integrals due to inertia, plasticity and interior support effect terms. By More >

T. Jarak¹, J. Sorić¹, J. Hoster¹

CMES-Computer Modeling in Engineering & Sciences, Vol.18, No.3, pp. 235-246, 2007, DOI:10.3970/cmes.2007.018.235

Abstract A meshless computational method based on the local Petrov-Galerkin approach for the analysis of shell structures is presented. A concept of a three dimensional solid, allowing the use of completely 3-D constitutive models, is applied. Discretization is carried out by using both a moving least square approximation and polynomial functions. The exact shell geometry can be described. Thickness locking is eliminated by using a hierarchical quadratic approximation over the thickness. The shear locking phenomena in case of thin structures and the sensitivity to rigid body motions are minimized by applying interpolation functions of sufficiently high More >

Q.W. Ma¹

CMES-Computer Modeling in Engineering & Sciences, Vol.18, No.3, pp. 223-234, 2007, DOI:10.3970/cmes.2007.018.223

Abstract Two methods have been recently developed by the author and his group: one called MLPG_R (Meshless Local Petrov-Galerkin method based on Rankine source solution) and the other called QALE-FEM (Quasi Arbitrary Lagrangian-Eulerian Finite Element Method). The former is a meshless method developed from a general MLPG (Meshless Local Petrov-Galerkin) method and is more computationally efficient than the general one when applied to modelling nonlinear water waves. The later is a mesh-based method similar to a conventional finite element method (FEM) when discretizing the governing equations but different from the conventional one in managing the mesh. More >