De-hao Yu¹ ,Hong-ying Huang²

CMES-Computer Modeling in Engineering & Sciences, Vol.35, No.3, pp. 227-252, 2008, DOI:10.3970/cmes.2008.035.227

Abstract In this paper, we apply the artificial boundary method to solve a three-dimensional nonlinear interface problem on an unbounded domain. A spherical or ellipsoidal surface as the artificial boundary is introduced. The exact artificial boundary conditions are derived explicitly in terms of an infinite series and then the well-posedness of the coupled weak formulation in a bounded domain, which is equivalent to the original problem in the unbounded domain, is obtained. The error estimate depends on the mesh size, the term after truncating the infinite series and the location of the artificial boundary. Some numerical More >

Maenghyo Cho¹, Jinbok Choi², Hee Yuel Roh³

CMES-Computer Modeling in Engineering & Sciences, Vol.33, No.1, pp. 17-48, 2008, DOI:10.3970/cmes.2008.033.017

Abstract The framework for the linkage between geometric modeling and an analysis based on the NURBS technology is developed in this study. In the present study, The NURBS surfaces were generated by interpolating a given set of data points or by extracting the necessary information to construct the NURBS surface from the IGES format file which was generated by the commercial CAD systems. Numerical examples showed the rate of displacement convergence for the various parameter-izations of the NURBS surface. Quadric surface, which is generated exactly by NURBS representation, was considered. One of the important advantages of More >

O.O. Oyekoya, D.U. Mba¹, A.M. El-Zafrany

CMES-Computer Modeling in Engineering & Sciences, Vol.34, No.1, pp. 55-86, 2008, DOI:10.3970/cmes.2008.034.055

Abstract In this paper, two new Mindlin-type plate bending elements have been derived for the modelling of functionally graded plate subjected to various loading conditions such as tensile loading, in-plane bending and out-of-plane bending. The properties of the first Mindlin-type element (i.e. Average Mindlin-type element) are computed by using an average fibre distribution technique which averages the macro-mechanical properties over each element. The properties of the second Mindlin-type element (i.e. Smooth Mindlin-type element) are computed by using a smooth fibre distribution technique, which directly uses the macro-mechanical properties at Gaussian quadrature points of each element. There More >

Patrícia C.T. Gonçalves^1,2, João Manuel R.S. Tavares^1,2, R.M. Natal Jorge^1,3

CMES-Computer Modeling in Engineering & Sciences, Vol.32, No.1, pp. 45-56, 2008, DOI:10.3970/cmes.2008.032.045

Abstract The main goals of the present work are to automatically extract the contour of an object and to simulate its deformation using a physical approach. In this work, to segment an object represented in an image, an initial contour is manually defined for it that will then automatically evolve until it reaches the border of the desired object. In this approach, the contour is modelled by a physical formulation using the finite element method, and its temporal evolution to the desired final contour is driven by internal and external forces. The internal forces are defined… More >

A. Takei¹, S. Yoshimura¹, H. Kanayama²

CMES-Computer Modeling in Engineering & Sciences, Vol.31, No.1, pp. 13-24, 2008, DOI:10.3970/cmes.2008.031.013

Abstract This paper presents large-scale finite element analyses of high frequency electromagnetic fields in commuter trains. The ADVENTURE_Magnetic is one of the main modules of the ADVENTURE system, which is an open source parallel finite element analyses system, and is able to solve eddy current and magnetostatic problems using the hierarchical domain decomposition method (HDDM) with an iterative linear algebraic solver. In this paper, we improve the module so as to solve a high frequency electromagnetic field of 500-1000 M[Hz]. A stationary Helmholtz equation for electromagnetic wave problems is solved taking an electric field as an More >

Atsuya Oishi¹, Shinobu Yoshimura²

CMES-Computer Modeling in Engineering & Sciences, Vol.28, No.2, pp. 127-146, 2008, DOI:10.3970/cmes.2008.028.127

Abstract This paper describes new methods based on genetic approaches for finding approximating expressions of local coordinates of a contact point in a local contact search process. A contact search process generally consists of the following two phases: a global search phase for finding the nearest node-segment pair and a local search phase for finding an exact local coordinate of the contact point within the segment. The local contact search can be regarded as the mapping from the coordinates of nodes to the local coordinates of contact points. In this paper, two methods are proposed to… More >

X.Y. Cui^1,2, G. R. Liu^2,3, G. Y. Li¹, X. Zhao², T.T. Nguyen², G.Y. Sun¹

CMES-Computer Modeling in Engineering & Sciences, Vol.28, No.2, pp. 109-126, 2008, DOI:10.3970/cmes.2008.028.109

Abstract A smoothed finite element method (SFEM) is presented to analyze linear and geometrically nonlinear problems of plates and shells using bilinear quadrilateral elements. The formulation is based on the first order shear deformation theory. In the present SFEM, the elements are further divided into smoothing cells to perform strain smoothing operation, and the strain energy in each smoothing cell is expressed as an explicit form of the smoothed strain. The effect of the number of divisions of smoothing cells in elements is investigated in detail. It is found that using three smoothing cells for bending More >

Sanjay Mittal¹

CMES-Computer Modeling in Engineering & Sciences, Vol.27, No.1&2, pp. 63-78, 2008, DOI:10.3970/cmes.2008.027.063

Abstract Flow past a circular cylinder looses stability at a Reynolds number,Re~47. It has been shown, in the past, that the linear stability analysis (LSA) of the steady state solution can predict not only the critical Re, but also the non-dimensional frequency, St, of the associated instability. For larger Re the non-linear effects become important and the LSA of the steady-state flow does not predict the correct St. It is shown that, in general, the LSA applied to the time-averaged flow can result in useful information regarding its stability. This idea is applied to the Re = 100 flow past More >

Y. Nishidate¹, G. P. Nikishkov^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.26, No.2, pp. 91-106, 2008, DOI:10.3970/cmes.2008.026.091

Abstract Atomic-scale finite element procedure for modeling of self-positioning nanostructures is developed. Our variant of the atomic-scale finite element method is based on a meshless approach and on the Tersoff interatomic potential function. The developed algorithm is used for determination of equilibrium configuration of atoms after nanostructure self-positioning. Dependency of the curvature radius of nanostructures on their thickness is investigated. It is found that for thin nanostructures the curvature radius is considerably smaller than predicted by continuum mechanics equations. Curvature radius variation with varying orientation of crystallographic axes is also modeled and results are compared to More >

Atsuya Oishi¹, Shinobu Yoshimura²

CMES-Computer Modeling in Engineering & Sciences, Vol.25, No.2, pp. 115-132, 2008, DOI:10.3970/cmes.2008.025.115

Abstract This paper describes the finite element analyses of dynamic problems using graphics hardware. The graphics hardware, known as GPU that is an acronym of Graphics Processing Unit, was first developed only for processing 3D computer graphics. However it has obtained both flexible programmability using a high-level shader programming language such as OpenGL Shading Language (GLSL), and has also obtained fast numerical processing ability of over 160 GFLOPS that is much faster than CPU. In this paper, GPU is utilized for the finite element analyses of dynamic problems. Two different computational tasks in the dynamic finite More >