A. Milazzo¹, I. Benedetti², C. Orlando³

CMES-Computer Modeling in Engineering & Sciences, Vol.15, No.1, pp. 17-30, 2006, DOI:10.3970/cmes.2006.015.017

Abstract A boundary integral formulation and its numerical implementation are presented for the analysis of magneto electro elastic media. The problem is formulated by using a suitable set of generalized variables, namely the generalized displacements, which are comprised of mechanical displacements and electric and magnetic scalar potentials, and generalized tractions, that is mechanical tractions, electric displacement and magnetic induction. The governing boundary integral equation is obtained by generalizing the reciprocity theorem to the magneto electro elasticity. The fundamental solutions are calculated through a modified Lekhnitskii's approach, reformulated in terms of generalized magneto-electro-elastic displacements. To assess the More >

Baohua Ji^*, Yonggang Huang^†

Molecular & Cellular Biomechanics, Vol.3, No.3, pp. 109-120, 2006, DOI:10.3970/mcb.2006.003.109

Abstract The mechanisms governing the self-assembled structure of biomolecules (single chain and bundle of chains) are studied with an AB copolymer model via the coarse grained molecular dynamics simulations. Non-local hydrophobic interaction is found to play a critical role in the pattern formation of the assembled structure of polymer chains. We show that the polymer structure could be controlled by adjusting the balance between local (short range) and non-local (long range) hydrophobic interaction which are influenced by various factors such as the sequences, chain length, stiffness, confinement, and the topology of polymers. In addition, the competition More >

B. Yang², V. K. Tewary³

CMES-Computer Modeling in Engineering & Sciences, Vol.15, No.3, pp. 165-178, 2006, DOI:10.3970/cmes.2006.015.165

Abstract Green's function (GF) modeling of defects may take effect only if the GF as well as its various integrals over a line, a surface and/or a volume can be efficiently evaluated. The GF is needed in modeling a point defect, while integrals are needed in modeling line, surface and volumetric defects. In a matrix of multilayered, generally anisotropic and linearly elastic and piezoelectric materials, the GF has been derived by applying 2D Fourier transforms and the Stroh formalism. Its use involves another two dimensions of integration in the Fourier inverse transform. A semi-analytical scheme has… More >

Taeyoung Ha¹, Sangwon Seo², Dongwoo Sheen³

CMES-Computer Modeling in Engineering & Sciences, Vol.14, No.1, pp. 57-76, 2006, DOI:10.3970/cmes.2006.014.057

Abstract We present nonoverlapping domain decomposition methods for the approximation of both electromagnetic fields in a three-dimensional bounded domain satisfying absorbing boundary conditions. A Seidel-type domain decomposition iterative method is introduced based on a hybridization of a nonconforming mixed finite element method. Convergence results for the numerical procedure are proved by introducing a suitable pseudo-energy. The spectral radius of the iterative procedure is estimated and a method for choosing an optimal matching parameter is given. A red-black Seidel-type method which is readily parallelizable is also introduced and analyzed. Numerical experiments confirm that the presented algorithms are faster than More >

L. Gao¹, K. Liu^1,2, Y. Liu³

CMES-Computer Modeling in Engineering & Sciences, Vol.12, No.3, pp. 181-196, 2006, DOI:10.3970/cmes.2006.012.181

Abstract A new numerical algorithm based on the Meshless Local Petrov-Galerkin approach is presented for analyzing the dynamic fracture problems in elastic media. To simplify the treatment of essential boundary condition, a novel modified Moving Least Square (MLS) procedure is proposed by introducing Lagrange multiplier into MLS procedure, which can perform both MLS approximation and interpolation in one approximation domain. The compact spline function is used as the test function in the local form of elasto-dynamic equations. For the feature of stress wave propagation, the coupled second-order ODEs respect to the time are solved by the… More >

T. Hasebe¹

CMES-Computer Modeling in Engineering & Sciences, Vol.11, No.3, pp. 145-156, 2006, DOI:10.3970/cmes.2006.011.145

Abstract This paper presents recent achievements in field theoretical approach toward substantial linkage among key hieratical scales dominating polycrystalline plasticity of metals and alloys. Major ingredients of the theory are briefly shown first, which is followed by several overwhelming results and some implications including key factors for dislocation cell structure evolution, key features of polycrystalline plasticity and their rational modeling in crystal plasticity-based constitutive equation. More >

C.P. Hong¹, M.F. Zhu², S.Y. Lee¹

FDMP-Fluid Dynamics & Materials Processing, Vol.2, No.4, pp. 247-260, 2006, DOI:10.3970/fdmp.2006.002.247

Abstract In typical solidification processes the flow of molten metal is usually regarded as an unavoidable phenomenon potentially affecting the morphology of dendritic growth. Fundamental understanding of such flow is thus important for predicting and controlling solidification microstructures. This paper presents numerical simulations on the evolution of dendritic microstructures with melt convection. A two-dimensional modified cellular automaton (MCA) coupled with a transport model is developed to simulate solidification of binary and ternary alloys in the presence of fluid flow. This model takes into account the effects of the constitutional undercooling and curvature undercooling on the equilibrium More >

M. Sohail¹, M. Z. Saghir¹

FDMP-Fluid Dynamics & Materials Processing, Vol.2, No.2, pp. 127-140, 2006, DOI:10.3970/fdmp.2006.002.127

Abstract A numerical simulation study is carried out for the crystal growth of Ge1-xSixby the Traveling Heater Method (THM). The effects of a geometrical misalignment on the crystal growth are investigated. The full Navier-Stokes equations together with the energy, mass transport and continuity equations are solved numerically using the finite element technique. The application of a misalignment is shown to have a considerable effect on the buoyancy induced flow. An optimal misalignment is determined, that weakens the convective flow, provides a uniform concentration along the growth interface and gives symmetrical characteristics to the three-dimensional buoyancy induced More >

L. Gornet¹, S. Marguet², G. Marckmann³

CMC-Computers, Materials & Continua, Vol.4, No.2, pp. 63-74, 2006, DOI:10.3970/cmc.2006.004.063

Abstract The purpose of the present study is to determine the components of the effective elasticity tensor and the failure properties of Nomex^® honeycomb cores. In order to carry out this study, the NidaCore software, a program dedicated to Nomex®Cores predictions, has been developed using the Finite Element tool Cast3M-CEA. This software is based on periodic homogenization techniques and on the modelling of structural instability phenomena. The homogenization of the periodic microstructure is realized thanks to a strain energy approach. It assumes the mechanical equivalence between the microstructures of a RVE and a similar homogeneous macroscopic volume.… More >

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

CMC-Computers, Materials & Continua, Vol.4, No.1, pp. 43-54, 2006, DOI:10.3970/cmc.2006.004.043

Abstract The Rajendran-Grove (RG) ceramic damage model is a three-dimensional internal variable based constitutive model for ceramic materials, with the considerations of micro-crack extension and void collapse. In the present paper, the RG ceramic model is implemented into the newly developed computational framework based on the Meshless Local Petrov-Galerkin (MLPG) method, for solving high-speed impact and penetration problems. The ability of the RG model to describe the internal damage evolution and the effective material response is investigated. Several numerical examples are presented, including the rod-on-rod impact, plate-on-plate impact, and ballistic penetration. The computational results are compared More >