Rade Vignjevic¹, Juan R. Reveles¹, James Campbell¹

CMES-Computer Modeling in Engineering & Sciences, Vol.14, No.3, pp. 181-198, 2006, DOI:10.3970/cmes.2006.014.181

Abstract To correct some of the main shortcomings of conventional SPH, a version of this method based on the Total Lagrangian formalism, T. Rabczuk, T. Belytschko and S. Xiao (2004), is developed. The resulting scheme removes the spatial discretisation instability inherent in conventional SPH, J. Monaghan (1992).
The Total Lagrangian framework is combined with the mixed correction which ensures linear completeness and compliance with the patch test, R. Vignjevic, J. Campbell, L. Libersky (2000). The mixed correction utilizes Shepard Functions in combination with a correction to derivative approximations.
Incompleteness of the kernel support combined with… More >

S. Noroozi¹, P. Sewell¹, J. Vinney¹

CMES-Computer Modeling in Engineering & Sciences, Vol.14, No.3, pp. 171-180, 2006, DOI:10.3970/cmes.2006.014.171

Abstract A method that combines a modified back propagation Artificial Neural Network (ANN) and Boundary Element Analysis (BEA) was introduced and discussed in the author's previous papers. This paper discusses the development of an automated inverse boundary element problem engine. This inverse problem engine can be applied to both potential and elastostatic problems.
In this study, BEA solutions of a two-dimensional potential problem is utilised to test the system and to train a back propagation Artificial Neural Network (ANN). Once training is completed and the transfer function is created, the solution to any subsequent or new… More >

L.S. Miers¹, J.C.F. Telles²

CMES-Computer Modeling in Engineering & Sciences, Vol.14, No.3, pp. 161-170, 2006, DOI:10.3970/cmes.2006.014.161

Abstract This work aims at extending the concept of the Numerical Green's Function (NGF), well known from boundary element applications to fracture mechanics, to the Local Boundary Integral Equation (LBIE) context. As a "companion" solution, the NGF is used to remove the integrals over the crack boundary and is introduced only for source points whose support touches or contains the crack. The results obtained with the coupling of NGF-LBIE in previous potential discontinuity Laplace's equation problems and the authors' experience in NGF-BEM fracture mechanics were the motivation for this development. More >

L.A. de Lacerda¹, J. M. da Silva¹

CMES-Computer Modeling in Engineering & Sciences, Vol.14, No.3, pp. 153-160, 2006, DOI:10.3970/cmes.2006.014.153

Abstract A two-dimensional boundary element formulation is presented and coupled to a genetic algorithm to identify polarization curves of buried slender structures. The dual boundary element method is implemented to model the cathodic protection of the metallic body and the genetic algorithm is employed to deal with the inverse problem of determining the non-linear polarization curve, which describes the relation between current density and electrochemical potential at the soil metal interface. In this work, this non-linear relation resulting from anodic and cathodic reactions is represented by a classical seven parameters expression. Stratified soil resistivity is modeled More >

S. N. Atluri¹, H. T. Liu², Z. D. Han²

CMES-Computer Modeling in Engineering & Sciences, Vol.14, No.3, pp. 141-152, 2006, DOI:10.3970/cmes.2006.014.141

Abstract The Meshless Local Petrov-Galerkin (MLPG) mixed collocation method is proposed in this paper, for solving elasticity problems. In the present MLPG approach, the mixed scheme is applied to interpolate the displacements and stresses independently, as in the MLPG finite volume method. To improve the efficiency, the local weak form is established at the nodal points, for the stresses, by using the collocation method. The traction boundary conditions are also imposed into the stress equations directly. It becomes very simple and straightforward to impose various boundary conditions, especially for the high-order PDEs. Numerical examples show that More >

S.-Y. Leu¹, J.T. Chen²

CMES-Computer Modeling in Engineering & Sciences, Vol.14, No.2, pp. 129-140, 2006, DOI:10.3970/cmes.2006.014.129

Abstract Plastic limit angular velocity of rotating hollow cylinders made of the von Mises materials with nonlinear isotropic hardening is investigated numerically and analytically in the paper. The paper applies sequential limit analysis to deal with the rotating problems involving hardening material property and weakening behavior resulted from the widening deformation. By sequential limit analysis, the paper treats the plasticity problems as a sequence of limit analysis problems stated in the upper bound formulation. Rigorous upper bounds are acquired iteratively through a computational optimization procedure with the angular velocity factor as the objective function. Especially, rigorous More >

Z. D. Han¹, H. T. Liu¹, A. M. Rajendran², S. N. Atluri³

CMES-Computer Modeling in Engineering & Sciences, Vol.14, No.2, pp. 119-128, 2006, DOI:10.3970/cmes.2006.014.119

Abstract This paper presents the implementation of a three-dimensional dynamic code, for contact, impact, and penetration mechanics, based on the Meshless Local Petrov-Galerkin (MLPG) approach. In the current implementation, both velocities and velocity-gradients are interpolated independently, and their compatibility is enforced only at nodal points. As a result, the time consuming differentiations of the shape functions at all integration points is avoided, and therefore, the numerical process becomes more stable and efficient. The ability of the MLPG code for solving high-speed contact, impact and penetration problems with large deformations and rotations is demonstrated through several computational More >

J. Ma², H. Lu², B. Wang², R. Hornung³, A. Wissink³, R. Komanduri^2,*

CMES-Computer Modeling in Engineering & Sciences, Vol.14, No.2, pp. 101-118, 2006, DOI:10.3970/cmes.2006.014.101

Abstract A new method for multiscale simulation bridging two scales, namely, the continuum scale using the generalized interpolation material point (GIMP) method and the atomistic scale using the molecular dynamics (MD), is presented and verified in 2D. The atomistic strain from the molecular dynamics simulation is determined through interpolation of the displacement field into an Eulerian background grid using the same generalized interpolation functions as that in the GIMP method. The atomistic strain is consistent with that determined from the virial theorem for interior points but provides more accurate values at the boundary of the MD… More >

Xianwu Ling¹, S.N. Atluri

CMES-Computer Modeling in Engineering & Sciences, Vol.14, No.2, pp. 91-100, 2006, DOI:10.3970/cmes.2006.014.091

Abstract In this paper, a lattice-based cell model is proposed for single wall carbon nanotubes (SWNTs). The finite temperature effect is accounted for via the local harmonic approach. The equilibrium SWNT configurations are obtained by minimizing the Helmholtz free energy with respect to seven primary coordinate variables that are subjected to a chirality constraint. The calculated specific heats agree well with the experimental data, and at low temperature depend on the tube radii with small tubes having much lower values. Our calculated coefficients of thermal expansion (CTEs) are universally positive for all the radial, axial and More >

Chein-Shan Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.14, No.2, pp. 77-90, 2006, DOI:10.3970/cmes.2006.014.077

Abstract In this paper we derive the first-order and second-order one-step GPS applied to the estimation of thermophysical properties. Solving the resultant algebraic equations, which usually converges within ten iterations, it is not difficult to estimate the unknown temperature-dependent thermal conductivity and heat capacity simultaneously, if some supplemented data of measured temperature at a time T is provided. When the measured temperature in the conducting slab is contaminated by noise, our estimated results are also good. The new method does not require any prior information on the functional forms of thermal conductivity and heat capacity. Numerical examples More >