M. Tanaka^1,2, H. Noguchi¹, M. Fujikawa^3,4, M. Sato³, S. Oi³, T. Kobayashi³, K. Furuichi⁵, S. Ishimaru⁵, C. Nonomura⁵

CMES-Computer Modeling in Engineering & Sciences, Vol.62, No.3, pp. 265-290, 2010, DOI:10.3970/cmes.2010.062.265

Abstract This paper describes the development of a proper constitutive model of woven fabrics and its implementation in nonlinear finite shell elements in order to simulate the large deformation behavior of cloth. This work currently focuses on a macroscopic continuum constitutive model that is capable of capturing the realistic mechanical behavior of cloth that is characterized by two families of yarns, i.e., warp and weft. In this study, two strategies are considered. One is a rebar layer model and the other is a polyconvex anisotropic hyperelastic material model. The latter avoids non-physical behavior and can consider More >

Mustafa Yidiz¹, Bayram Heydarov², İsmet Gölgeleyen¹

CMES-Computer Modeling in Engineering & Sciences, Vol.62, No.3, pp. 255-264, 2010, DOI:10.3970/cmes.2010.062.255

Abstract We investigate the solvability of an inverse problem for the non-stationary general kinetic equation. We also obtained the approximate solution of this problem by using symbolic computation. A comparison between the approximate solution and the exact solution of the problem is presented. More >

Huan J. Keh¹, Yu C. Chang²

CMES-Computer Modeling in Engineering & Sciences, Vol.62, No.3, pp. 225-254, 2010, DOI:10.3970/cmes.2010.062.225

Abstract A theoretical study of the Stokes flow caused by a rigid particle of revolution translating axisymmetrically perpendicular to two parallel plane walls at an arbitrary position between them in a viscous fluid, which may slip at the particle surface, is presented. A method of distribution of a set of spherical singularities along the axis of revolution within a prolate particle or on the fundamental plane within an oblate particle is used to find the general solution of the fluid velocity field that satisfies the boundary conditions at the plane walls and at infinity. The slip… More >

M. Ferronato¹, A. Mazzia¹, G. Pini¹

CMES-Computer Modeling in Engineering & Sciences, Vol.62, No.2, pp. 205-224, 2010, DOI:10.3970/cmes.2010.062.205

Abstract In the engineering practice meshing and re-meshing complex domains by Finite Elements (FE) is one of the most time-consuming efforts. Meshless methods avoid this task but are computationally more expensive than standard FE. A somewhat natural improvement can be attempted by combining the two techniques with the aim at emphasizing the respective merits. The present work describes a FE enrichment by the Meshless Local Petrov-Galerkin (MLPG) method. The basic idea is to add a limited number of moving MLPG points over a fixed coarse FE grid, in order to improve the solution accuracy in specific More >

H.M. Shodja^1,2,3, M. Khezri⁴, A. Hashemian¹, A. Behzadan¹

CMES-Computer Modeling in Engineering & Sciences, Vol.62, No.2, pp. 171-204, 2010, DOI:10.3970/cmes.2010.062.171

Abstract An accurate numerical methodology for capturing the field quantities across the interfaces between material discontinuities, in the context of reproducing kernel particle method (RKPM), is of particular interest. For this purpose the innovative numerical technique, so-called augmented corrected collocation method is introduced; this technique is an extension of the corrected collocation method used for imposing essential boundary conditions (EBCs). The robustness of this methodology is shown by utilizing it to solve two benchmark problems of material discontinuities, namely the problem of circular inhomogeneity with uniform radial eigenstrain, and the problem of interaction between a crack More >

Xianqiao Wang¹, James D. Lee¹

CMES-Computer Modeling in Engineering & Sciences, Vol.62, No.2, pp. 150-170, 2010, DOI:10.3970/cmes.2010.062.150

Abstract The concurrent methods for coupling molecular dynamics with continuum thermodynamics offer a myriad of challenging problems, mostly related with energy transmission, wave reflection, and damage propagation at the interfaces between the continuum description and the discrete description. In this work, by virtue of the atomistic field theory (AFT), we present an analysis to reconcile the compatibility between atomic region and continuum region and to calculate the matching temperature field of a heat conduction problem in a concurrent atomistic/continuum system. First, formulation of AFT with finite temperature and its corresponding finite element implementation are briefly introduced. More >

A. Arefmanesh¹, M. Najafi², M. Nikfar³

CMES-Computer Modeling in Engineering & Sciences, Vol.62, No.2, pp. 113-149, 2010, DOI:10.3970/cmes.2010.062.113

Abstract As some new applications of the meshless local Petrov-Galerkin method (MLPG) with unity as the test function, a number of buoyancy-driven fluid flow natural convection heat transfer problems in cavities with differentially-heated wavy side walls were analyzed. Cavities with a single wavy wall on one side as well as two wavy walls erected on both sides were considered. For the cases of the double wavy walls, two different configurations in terms of the two walls facing each other on the two sides of the cavities symmetrically or non-symmetrically were investigated. All the simulations performed in… More >

P.H. Gaskell¹, Y.C. Lee², H.M. Thompson¹

CMES-Computer Modeling in Engineering & Sciences, Vol.62, No.1, pp. 77-112, 2010, DOI:10.3970/cmes.2010.062.077

Abstract The three-dimensional flow of a gravity-driven continuous thin liquid film on substrates containing micro-scale features is modelled using the long-wave lubrication approximation, encompassing cases when surface topography is either engulfed by the film or extends all the way though it. The discrete analogue of the time-dependent governing equations is solved accurately using a purpose designed multigrid strategy incorporating both automatic error-controlled adaptive time stepping and local mesh refinement/de-refinement. Central to the overall approach is a Newton globally convergent algorithm which greatly simplifies the inclusion of further physics via the solution of additional equations of the More >

Roman Chapko¹, B. Tomas Johansson², Oleh Sobeyko¹

CMES-Computer Modeling in Engineering & Sciences, Vol.62, No.1, pp. 57-76, 2010, DOI:10.3970/cmes.2010.062.057

Abstract We propose an iterative procedure for the inverse problem of determining the displacement vector on the boundary of a bounded planar inclusion given the displacement and stress fields on an infinite (planar) line-segment. At each iteration step mixed boundary value problems in an elastostatic half-plane containing the bounded inclusion are solved. For efficient numerical implementation of the procedure these mixed problems are reduced to integral equations over the bounded inclusion. Well-posedness and numerical solution of these boundary integral equations are presented, and a proof of convergence of the procedure for the inverse problem to the More >

E.A. de Souza Neto¹, S. Amstutz², S.M. Giusti³, A.A. Novotny³

CMES-Computer Modeling in Engineering & Sciences, Vol.62, No.1, pp. 23-56, 2010, DOI:10.3970/cmes.2010.062.023

Abstract A recently proposed algorithm for micro-structural optimization, based on the concept of topological derivative and a level-set domain representation, is applied to the synthesis of elastic and heat conducting bi-material micro-structures. The macroscopic properties are estimated by means of a family of multi-scale constitutive theories where the macroscopic strain and stress tensors (temperature gradient and heat flux vector in the heat conducting case) are defined as volume averages of their microscopic counterparts over a Representative Volume Element (RVE). Several finite element-based examples of micro-structural optimization are presented. Three multi-scale models, providing an upper and a More >