P.H. Wen¹, M.H. Aliabadi², Y.W. Liu³

CMES-Computer Modeling in Engineering & Sciences, Vol.30, No.3, pp. 133-148, 2008, DOI:10.3970/cmes.2008.030.133

Abstract Based on the variation of potential energy, the element-free Galerkin method (MFGM) has been investigated for structures with crack on the basis of radial base function interpolation. An enriched radial base function is introduced to capture the singularities of stress at the crack tips. The advantages of the finite element method are remained in this method and there is a significant improvement of accuracy, particularly for the crack problems of fracture mechanics. The applications of the element-free Galerkin method with enriched radial base function to two-dimensional fracture mechanics in functionally graded materials have been presented More >

Jeng-Tzong Chen¹, Jia-Nan Ke

CMES-Computer Modeling in Engineering & Sciences, Vol.29, No.3, pp. 111-136, 2008, DOI:10.3970/cmes.2008.029.111

Abstract A null-field integral equation is employed to derive the two-dimensional antiplane dynamic Green's functions for a circular inclusion with an imperfect interface. We employ the linear spring model with vanishing thickness to characterize the imperfect interface. Analytical expressions of displacement and stress fields due to time-harmonic antiplane line forces located either in the unbounded matrix or in the circular inclusion are presented. To fully capture the circular geometries, degenerate- kernel expressions of fundamental solutions in the polar coordinate and Fourier series for boundary densities are adopted. Good agreement is made after comparing with the analytical More >

Yongjie Zhang¹, Boyle C. Cheng², Changho Oh¹, Jessica L. Spehar², James Burgess³

CMES-Computer Modeling in Engineering & Sciences, Vol.29, No.2, pp. 55-62, 2008, DOI:10.3970/cmes.2008.029.055

Abstract The spinal column plays a vital biomechanical role in the human body by providing structural support and facilitating motion. As degenerative changes occur in the spine, surgical treatment may be necessary in certain instances. Such treatments seek to address pain, frequently through the restriction of spinal motion. Traditional spinal implant devices are designed to restrict the motion of a functional spinal unit (FSU) but newer device designs allow for semi-constrained motion such as spinal arthroplasty devices. In this study, a sequence of fluoroscopic imaging data was recorded during the flexibility protocol with an interspinous process… More >

F. Cosmi¹

CMES-Computer Modeling in Engineering & Sciences, Vol.27, No.3, pp. 187-196, 2008, DOI:10.3970/cmes.2008.027.187

Abstract The Cell Method is a recent numerical method that can be applied in several fields of physics and engineering. In this paper, the elastodynamics formulation is extended to include system internal damping, highlighting some interesting characteristics of the method. The developed formulation leads to an explicit solving system. The mass matrix is diagonal (without lumping) and in the most general case a time-dependent damping coefficient can be defined for each node. \newline Accuracy and convergence rate have been tested with reference to the classical problem of a particle free vibration with viscous damping.
An application More >

B Radhika¹, S S P,a¹, C S Manohar^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.27, No.1&2, pp. 79-110, 2008, DOI:10.3970/cmes.2008.027.079

Abstract Reliability of nonlinear vibrating systems under stochastic excitations is investigated using a two-stage Monte Carlo simulation strategy. For systems with white noise excitation, the governing equations of motion are interpreted as a set of Ito stochastic differential equations. It is assumed that the probability distribution of the maximum in the steady state response belongs to the basin of attraction of one of the classical asymptotic extreme value distributions. The first stage of the solution strategy consists of selection of the form of the extreme value distribution based on hypothesis tests, and the next stage involves More >

Jyoti Ranjan Majhi, Ranjan Ganguli¹

CMES-Computer Modeling in Engineering & Sciences, Vol.27, No.1&2, pp. 25-36, 2008, DOI:10.3970/cmes.2008.027.025

Abstract The flapping equation for a rotating rigid helicopter blade is typically derived by considering 1) small flap angle, 2) small induced angle of attack and 3) linear aerodynamics. However, the use of nonlinear aerodynamics can make the assumptions of small angles suspect. A general equation describing helicopter blade flap dynamics for large flap angle and large induced inflow angle of attack is derived in this paper with nonlinear aerodynamics . Numerical simulations are performed by solving the nonlinear flapping ordinary differential equation for steady state conditions and the validity of the small angle approximations are More >

Wenzhong Tang¹, Suresh G. Advani¹

CMES-Computer Modeling in Engineering & Sciences, Vol.25, No.3, pp. 149-164, 2008, DOI:10.3970/cmes.2008.025.149

Abstract In this paper, a method for studying nanotube dynamics in simple shear flow was developed. A nanotube was described as a flexible fiber with a sphere-chain model. The forces on the nanotube were obtained by molecular dynamics simulations. The motion of the nanotube in simple shear flow was tracked by the flexible fiber dynamics method [Tang and Advani (2005)]. The viscosity of dilute nanotube suspensions was calculated based on the nanotube dynamics, and the effects of the aspect ratio and initial curvature of the nanotube on suspension viscosity are explored and discussed. 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 >

P. Le¹, N. Mai-Duy², T. Tran-Cong³, G. Baker⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.25, No.1, pp. 43-68, 2008, DOI:10.3970/cmes.2008.025.043

Abstract This paper describes an integrated radial basis function network (IRBFN) method for the numerical modelling of the dynamics of strain localization due to strain softening in quasi-brittle materials. The IRBFN method is a truly meshless method that is based on an unstructured point collocation procedure. We introduce a new and effective regularization method to enhance the performance of the IRBFN method and alleviate the numerical oscillations associated with weak discontinuity at the elastic wave front. The dynamic response of a one dimensional bar is investigated using both local and non-local continuum models. Numerical results, which More >

Arun K. Nair¹, Diana Farkas², Ronald D. Kriz¹

CMES-Computer Modeling in Engineering & Sciences, Vol.24, No.2&3, pp. 239-248, 2008, DOI:10.3970/cmes.2008.024.239

Abstract The indentation response of Ni thin films of thicknesses in the nano scale was studied using molecular dynamics simulations with embedded atom method (EAM) interatomic potentials. Simulations were performed in single crystal films in the [111] orientation with thicknesses of 7nm and 33nm. In the elastic regime, the loading curves observed start deviating from the Hertzian predictions for indentation depths greater than 2.5% of the film thickness. The observed loading curves are therefore dependent on the film thickness. The simulation results also show that the contact stress necessary to emit the first dislocation under the More >