Alexandre S. Botta¹, Wilson S. Venturini²

CMES-Computer Modeling in Engineering & Sciences, Vol.8, No.1, pp. 15-28, 2005, DOI:10.3970/cmes.2005.008.015

Abstract In this work a regularized boundary-finite element combination is proposed to analyse 2D elastostatic solids reinforced by fibres. The boundary element is adopted to model the matrix behaviour, while finite elements model the embedded fibres. The debonding effects caused by the adherence loss between the two materials are also considered. A three-degree polynomial is adopted to approach the displacement field along the fibre elements, while linear approximations are used to represent the bonding forces between fibres and the matrix. The non-linear debonding model is governed by a loading function written in terms of the contact More >

Luming Shen¹, Zhen Chen²

CMES-Computer Modeling in Engineering & Sciences, Vol.7, No.3, pp. 305-320, 2005, DOI:10.3970/cmes.2005.007.305

Abstract To simulate the dynamic responses involving different material phases in a finite computational domain without discretizing the whole problem domain, a silent boundary scheme is proposed within the framework of the material point method (MPM) that is an extension from Computational Fluid Dynamics to Computational Solid Dynamics. Because the MPM does not employ fixed mesh connectivity, a robust spatial discretization procedure in the moving domain of influence could be designed by applying viscous damping forces along the computational boundary. To establish a simple interface between the discretization procedures with and without fixed mesh connectivity, a More >

Hiroshi Kadowaki¹, Wing Kam Liu²

CMES-Computer Modeling in Engineering & Sciences, Vol.7, No.3, pp. 269-282, 2005, DOI:10.3970/cmes.2005.007.269

Abstract A method to derive governing equations and elastic-plastic constitutive relations for the micropolar continuum model is proposed. Averaging procedures are operated over a surrounding sub-domain for each material point to bridge a discrete microstructure to a macro continuum model. Material parameters are determined by these procedures. The size of the sub-domain represents the material intrinsic length scale, and it is passed into the macroscopic governing equation so that the numerical solution can be regularized for analyses of failure phenomena. An application to a simple granular material model is presented. More >

V. Sladek¹, J. Sladek¹, M. Tanaka²

CMES-Computer Modeling in Engineering & Sciences, Vol.7, No.1, pp. 69-84, 2005, DOI:10.3970/cmes.2005.007.069

Abstract An efficient numerical method is proposed for 2-d potential problems in anisotropic media with continuously variable material coefficients. The method is based on the local integral equations (utilizing a fundamental solution) and meshfree approximation of field variable. A lot of numerical experiments are carried out in order to study the numerical stability, accuracy, convergence and efficiency of several approaches utilizing various interpolations. More >

S. Deguchi^1,2, T. Ohashi², M. Sato²

Molecular & Cellular Biomechanics, Vol.2, No.3, pp. 125-134, 2005, DOI:10.3970/mcb.2005.002.125

Abstract Actin bundles in vascular endothelial cells (ECs) play a critical role in transmitting intracellular forces between separate focal adhesion sites. However, quantitative descriptions of tension level in single actin bundles in a physiological condition are still poorly studied. Here, we evaluated magnitude of preexisting tension in a single actin bundle of ECs on the basis of measurements of its preexisting stretching strain and tensile properties. Cultured ECs expressing fluorescently-labeled actin were treated with detergents to extract acin bundles. One end of an actin bundle was then dislodged from the substrate by using a microneedle, resulting… More >

C. Shu^1,2, W. X. Wu², C. M. Wang³

CMC-Computers, Materials & Continua, Vol.2, No.3, pp. 189-200, 2005, DOI:10.3970/cmc.2005.002.189

Abstract This paper demonstrates the application of a meshfree least square-based finite difference (LSFD) method for analysis of metallic waveguides. The waveguide problem is an eigenvalue problem that is governed by the Helmholtz equation. The second order derivatives in the Helmholtz equation are explicitly approximated by the LSFD formulations. TM modes and TE modes are calculated for some metallic waveguides with different cross-sectional shapes. Numerical examples show that the LSFD method is a very efficient meshfree method for waveguide analysis with complex domains. More >

B. Amaziane¹, A. Naji², D. Ouazar³, A. H.-D. Cheng⁴

CMC-Computers, Materials & Continua, Vol.2, No.2, pp. 85-96, 2005, DOI:10.3970/cmc.2005.002.085

Abstract In this paper the optimization of groundwater pumping in coastal aquifers under the threat of saltwater intrusion is investigated. The aquifer is inhomogeneous and contains several hydraulic conductivities zones. The aquifer data such as the hydraulic conductivities are uncertain, but with their expected mean and standard deviation values given. A stochastic boundary element method based on the perturbation technique is employed as the simulation tool. The stochastic optimization is handled by the chance-constrained programming. Genetic algorithm is selected as the optimization tool. Numerical examples of deterministic and stochastic problems are provided to demonstrate the feasibility More >

D. Soares Jr.¹, W.J. Mansur^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.8, No.2, pp. 153-164, 2005, DOI:10.3970/cmes.2005.008.153

Abstract A coupling procedure is described to perform time-domain numerical analyses of dynamic fluid-structure interaction. The fluid sub-domains, where acoustic waves propagate, are modeled by the Boundary Element Method (BEM), which is quite suitable to deal with linear homogeneous unbounded domain problems. The Finite Element Method (FEM), on the other hand, models the structure sub-domains, adopting a time marching scheme based on implicit Green's functions. The BEM/FEM coupling algorithm here developed is very efficient, eliminating the drawbacks of standard and iterative coupling procedures. Stability and accuracy features are improved by the adoption of different time steps More >

J. Ma¹, H. Lu¹, B. Wang¹, S. Roy¹, R. Hornung², A. Wissink², R. Komanduri^1,3

CMES-Computer Modeling in Engineering & Sciences, Vol.8, No.2, pp. 135-152, 2005, DOI:10.3970/cmes.2005.008.135

Abstract In the simulation of a wide range of mechanics problems including impact/contact/penetration and fracture, the material point method (MPM), Sulsky, Zhou and Shreyer (1995), demonstrated its computational capabilities. To resolve alternating stress sign and instability problems associated with conventional MPM, Bardenhagen and Kober (2004) introduced recently the generalized interpolation material point (GIMP) method and implemented for one-dimensional simulations. In this paper we have extended GIMP to 2D and applied to simulate simple tension and indentation problems. For simulations spanning multiple length scales, based on the continuum mechanics approach, we present a parallel GIMP computational method… More >

Y. Tomita ¹, M. Uchida ¹

CMES-Computer Modeling in Engineering & Sciences, Vol.10, No.3, pp. 239-248, 2005, DOI:10.3970/cmes.2005.010.239

Abstract In this study, we clarified the micro- to mesoscopic deformation behavior of a semicrystalline polymer by employing a large-deformation finite element homogenization method. The crystalline plasticity theory with a penalty method for the inextensibility of the chain direction and the nonaffine molecular chain network theory were applied for the representation of the deformation behavior of the crystalline and amorphous phases, respectively, in the composite microstructure of the semicrystalline polymer. The 3D structure of lamellae in the spherulite of high-density polyethylene was modeled, and the tensile and compressive deformation behaviors were investigated. A series of computational More >