George A. Gravvanis¹, Christos K. Filelis-Papadopoulos¹, Paschalis I.Matskanidis¹

CMES-Computer Modeling in Engineering & Sciences, Vol.100, No.4, pp. 323-345, 2014, DOI:10.3970/cmes.2014.100.323

Abstract Since the introduction of the Algebraic MultiGrid algorithm (AMG) over twenty years ago, significant progress has been made in improving the coarsening and the convergence behavior of the method. In this paper, an AMG method is introduced that utilizes a new generic approximate inverse algorithm as a smoother in conjunction with common coarsening techniques, such as classical Ruge-Stüben coarsening, CLJP and PMIS coarsening. The proposed approximate inverse scheme, namely Generic Approximate Banded Inverse (GenAbI), is a banded approximate inverse based on Incomplete LU factorization with zero fill–in (ILU(0)). The new class of Generic Approximate Banded Inverse can be computed for… More >

Soleiman Ghouhestani¹, Farzad Shahabian¹, Seyed Mahmoud Hosseini^2,3

CMES-Computer Modeling in Engineering & Sciences, Vol.100, No.4, pp. 295-321, 2014, DOI:10.3970/cmes.2014.100.295

Abstract In this paper, the meshless local Petrov-Galerkin (MLPG) method is exploited for dynamic analysis of functionally graded nanocomposite cylindrical layered structure reinforced by carbon nanotube subjected to mechanical shock loading. The carbon nanotubes (CNTs) are distributed across radial direction on thickness of cylinder, which can be simulated by linear and nonlinear volume fraction. Free vibration and elastic wave propagation are studied for various value of volume fraction exponent at various time intervals. The layered cylinder is assumed to be under axisymmetric and plane strain conditions. Four types of CNTs distributions including uniform and three kinds of functionally graded distributions along… More >

Yidong Bao^1,2, Dongmei Wu^1,3

CMES-Computer Modeling in Engineering & Sciences, Vol.100, No.4, pp. 277-294, 2014, DOI:10.3970/cmes.2014.100.277

Abstract Aiming at the shortcomings of mass-spring model, this paper, on the basis of preliminary studies, established a new viscoelastic soft tissue model based on meshless structure. The model is consisted of a large quantity of filled spheres, with every three spheres being connected by a spring and a Kelvin structure, which can further enhance the real-time virtual simulation operability while ensure the viscoelasticity of basic model. The stress relaxation and creep equation of the model can be derived from formula derivation. Through setting different parameters to the filled spheres, this model, with certain universal property, can create a virtual liver,… More >

T.A. Elgohary^1,2, L. Dong³, J.L. Junkins^2,4, S.N. Atluri^1,4

CMES-Computer Modeling in Engineering & Sciences, Vol.100, No.3, pp. 249-275, 2014, DOI:10.3970/cmes.2014.100.249

Abstract In this study, we consider Initial Value Problems (IVPs) for strongly nonlinear dynamical systems, and study numerical methods to analyze short as well as long-term responses. Dynamical systems characterized by a system of second-order nonlinear ordinary differential equations (ODEs) are recast into a system of nonlinear first order ODEs in mixed variables of positions as well as velocities. For each discrete-time interval Radial Basis Functions (RBFs) are assumed as trial functions for the mixed variables in the time domain. A simple collocation method is developed in the time-domain, with Legendre-Gauss-Lobatto nodes as RBF source points as well as collocation points.… More >

Q.W. Yang^1,2, S.G. Du¹, C.F. Liang¹, L.J. Yang¹

CMES-Computer Modeling in Engineering & Sciences, Vol.100, No.3, pp. 223-248, 2014, DOI:10.3970/cmes.2014.100.223

Abstract Although the model-independent damage localization algorithms have been extensively developed in recent years, the theoretical relationship between these damage indicators and the definition of damage is not clear. Moreover the existing damage localization methods are usually dependent on the boundary conditions and the type of structure. In view of this, the paper presents a universal model-independent algorithm for structural damage localization. To this end, the explicit relationship between the damage and damage-induced displacement variation is firstly clarified by using the well-known Sherman-Morrison and Woodbury formulas. A theorem is then presented for structural damage localization. According to the theorem, the universal… More >

Da-Guang Zhang^1,2, Hao-Miao Zhou¹

CMES-Computer Modeling in Engineering & Sciences, Vol.100, No.3, pp. 201-222, 2014, DOI:10.3970/cmes.2014.100.201

Abstract A model of FGM beams resting on nonlinear elastic foundations is put forward by physical neutral surface and high-order shear deformation theory. Material properties are assumed to be temperature dependent and von Kármán strain-displacement relationships are adopted. Nonlinear bending and thermal postbuckling are given by multi-term Ritz method, and influences played by different supported boundaries, thermal environmental conditions, different elastic foundations, and volume fraction index are discussed in detail. It is worth noting that the effect of nonlinear elastic foundation increases with increasing deflection. More >

S. Liu¹, X. Liu^2,3, Y. Yuan², H. A. Mang⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.100, No.3, pp. 157-200, 2014, DOI:10.3970/cmes.2014.100.157

Abstract A nonlinear viscoelastic constitutive model for early age concrete is presented in this paper. In this model, time-dependent properties, such as the elastic modulus, and thermal and autogenous shrinkage deformations, are computed by a stochastic multi-scale method, in which three different scales are specified according to the requirement of separation of scales, and different scales are linked by means of the asymptotic expansion theory with the help of specific representative volume elements (RVE). Thus, a cross-scale research from the cement paste to the macro structure of concrete is realized, and performance-based optimization of cement-based materials becomes possible. The developed constitutive… More >

Peng Weihong^1,2, Gao Feng¹, Cao Guohua³, Xu Yong², Cheng Hongmei¹

CMES-Computer Modeling in Engineering & Sciences, Vol.100, No.2, pp. 133-155, 2014, DOI:10.3970/cmes.2014.100.133

Abstract In this paper, the natural boundary element method is used to solve two-dimensional steady-state incompressible Stokes flows in a driven cavity. The analytical functions are expressed for the Stokes problem in an exterior circular domain under single value conditions, which satisfy the Stokes equations’ solutions in the form of complex functions. In order to obtain a uniform integral formula, the velocities on the boundary are expanded into Laurent series, and then compared with the analytical solutions obtained as described above. In this manner, the coefficients of the analytical solutions in the form of complex function are further confirmed. According to… More >

Wen-Hwa Chen^1,2, Ching-Feng Yu¹, Hsien-Chie Cheng^2,3

CMES-Computer Modeling in Engineering & Sciences, Vol.100, No.2, pp. 119-131, 2014, DOI:10.3970/cmes.2014.100.119

Abstract The aim of the study is to investigate the interactions between Sn adatoms in a solder bump and three typical Sn-based intermetallic compounds (IMCs) surface, i.e., Cu₃Sn, Cu₆Sn₅, and Ni₃Sn₄, at the atomistic scale. The adsorption energy, average bond length, and bond population of the Sn/Cu₃Sn, Sn/Cu₆Sn₅,and Sn/Ni₃Sn₄ systems are calculated through the first-principles density functional theory (DFT) calculation to investigate how the Sn adatoms influence the IMC surface. The calculated results show that the Sn atoms on the Cu₃Sn (0 0 1) surface hold the largest adsorption energy, average bond length and bond population, implying that the Cu₃Sn (0… More >

J. Useche¹

CMES-Computer Modeling in Engineering & Sciences, Vol.100, No.2, pp. 105-118, 2014, DOI:10.3970/cmes.2014.100.105

Abstract In this work, the harmonic analysis of shallow shells using the Boundary Element Method, is presented. The proposed boundary element formulation is based on a direct time-domain integration using the elastostatic fundamental solutions for both in-plane elasticity and shear deformable plates. Shallow shell was modeled coupling boundary element formulation of shear deformable plate and two-dimensional plane stress elasticity. Effects of shear deformation and rotatory inertia were included in the formulation. Domain integrals related to inertial terms were treated using the Dual Reciprocity Boundary Element Method. Numerical examples are presented to demonstrate the efficiency and accuracy of the proposed formulation. More >