Leiting Dong¹, Satya N. Atluri1¹

CMC-Computers, Materials & Continua, Vol.30, No.1, pp. 39-82, 2012, DOI:10.3970/cmc.2012.030.039

Abstract In this paper, as an extension to the authors's work in [Dong and Atluri (2011a,b, 2012a,b,c)], three-dimensional Trefftz Voronoi Cells (TVCs) with ellipsoidal voids/inclusions are developed for micromechanical modeling of heterogeneous materials. Several types of TVCs are developed, depending on the types of heterogeneity in each Voronoi Cell(VC). Each TVC can include alternatively an ellipsoidal void, an ellipsoidal elastic inclusion, or an ellipsoidal rigid inclusion. In all of these cases, an inter-VC compatible displacement field is assumed at each surface of the polyhedral VC, with Barycentric coordinates as nodal shape functions. The Trefftz trial displacement… More >

Peter L. Bishay¹, Jan Sladek², Vladimir Sladek², Satya N. Atluri¹

CMC-Computers, Materials & Continua, Vol.29, No.3, pp. 213-262, 2012, DOI:10.3970/cmc.2012.029.213

Abstract A new class of hybrid/mixed finite elements, denoted "HMFEM-C", has been developed for modeling magneto-electro-elastic (MEE) materials. These elements are based on assuming independent strain-fields, electric and magnetic fields, and collocating them with the strain-fields, electric and magnetic fields derived from the primal variables (mechanical displacements, electric and magnetic potentials) at some cleverly chosen points inside each element. The newly developed elements show significantly higher accuracy than the primal elements for the electric, magnetic as well as the mechanical variables. HMFEM-C is invariant through the use of the element-fixed local orthogonal base vectors, and is… More >

L. Dong¹, S. N. Atluri¹

CMC-Computers, Materials & Continua, Vol.29, No.2, pp. 169-212, 2012, DOI:10.3970/cmc.2012.029.169

Abstract In this paper, three-dimensionalT-Trefftz Voronoi Cell Finite Elements (VCFEM-TTs) are developed for micromechanical modeling of heterogeneous materials. Several types of VCFEMs are developed, depending on the types of heterogeneity in each element. Each VCFEM can include alternatively a spherical void, a spherical elastic inclusion, a spherical rigid inclusion, or no voids/inclusions at all.In all of these cases, an inter-element compatible displacement field is assumed at each surface of the polyhedral element, with Barycentric coordinates as nodal shape functions.The T-Trefftz trial displacement fields in each element are expressed in terms of the Papkovich-Neuber solution. Spherical harmonics… More >

Yan Gu¹, Wen Chen^1,2, Xiao-Qiao He³

CMC-Computers, Materials & Continua, Vol.29, No.2, pp. 129-154, 2012, DOI:10.3970/cmc.2012.029.129

Abstract This paper applies an improved singular boundary method (SBM) in conjunction with domain decomposition technique to stress analysis of layered elastic materials. For problems under consideration, the interface continuity conditions are approximated in the same manner as the boundary conditions. The multi-layered coating system is decomposed into multiple subdomains in terms of each layer, in which the solution is approximated separately by the SBM representation. The singular boundary method is a recent meshless boundary collocation method, in which the origin intensity factor plays a key role for its accuracy and efficiency. This study also introduces More >

Chia-Cheng Tsai¹, Po-Ho Lin²

CMC-Computers, Materials & Continua, Vol.28, No.3, pp. 231-260, 2012, DOI:10.3970/cmc.2012.028.231

Abstract Recently, Liu (CMES 21(2007), 53) developed the modified collocation Trefftz method (MCTM) by setting a characteristic length slightly larger than the maximum radius of the computational domain. In this study, we find that the range of admissible characteristic length can be significantly enlarged if the LU decomposition is applied for solving the resulted dense unsymmetric matrix. Furthermore, we discover a range formula for admissible characteristic length, in which the number of the T-complete functions, the shape of the computation domain, and the exponent bits of the involved floating-point arithmetic have been taken into consideration. In… More >

J.Y. Zhang¹, F.Z. Wang^2,3

CMES-Computer Modeling in Engineering & Sciences, Vol.88, No.2, pp. 141-154, 2012, DOI:10.3970/cmes.2012.088.141

Abstract The boundary knot method (BKM) is a kind of boundary-type meshless method, only boundary points are needed in the solution process. Since the BKM is mathematically simple and easy to implement, it is superior in dealing with Helmholtz problems with high wavenumbers and high dimensional problems. In this paper, we give an overview of the traditional BKM with collocation approach and provide three novel approaches for the BKM, as far as they are relevant for the other boundary-type techniques. The promising research directions are expected from an improved BKM aspect. More >

G. C. Bourantas¹, V. C. Loukopoulos²

CMES-Computer Modeling in Engineering & Sciences, Vol.86, No.4, pp. 275-300, 2012, DOI:10.3970/cmes.2012.086.275

Abstract This paper formulates a local Radial Basis Functions (LRBFs) collocation method for the numerical solution of the non-linear dispersive and dissipative KdV-Burger's (KdVB) equation. This equation models physical problems, such as irrotational incompressible flow, considering a shallow layer of an inviscid fluid moving under the influence of gravity and the motion of solitary waves. The local type of approximations used, leads to sparse algebraic systems that can be solved efficiently. The Inverse Multiquadrics (IMQ), Gaussian (GA) and Multiquadrics (MQ) Radial Basis Functions (RBF) interpolation are employed for the construction of the shape functions. Accuracy of More >

H.W. Tang¹, Y.T. Yang¹, C.K. Chen¹

CMES-Computer Modeling in Engineering & Sciences, Vol.85, No.3, pp. 193-206, 2012, DOI:10.3970/cmes.2012.085.193

Abstract In this article, a new method combining the Mathematical Programming and the Method of Weighted Residual called MP-MWR is presented. Under the validation of maximum principle, and up on the collocation method, the differential equation can be transferred into a bilateral inequality problem. Applying the genetic algorithms helps to find optimal solutions of upper and lower bounds which satisfy the inequalities. Here, the method is verified by analyzing the deflection of superelliptical clamped plate problem. By using this method, the good approximate solution and its error bounds can be obtained effectively and accurately. More >

Hong-Hua Dai^1,2, Matt Schnoor², Satya N. Atluri²

CMES-Computer Modeling in Engineering & Sciences, Vol.84, No.5, pp. 459-498, 2012, DOI:10.3970/cmes.2012.084.459

Abstract In this study, the harmonic and 1/3 subharmonic oscillations of a single degree of freedom Duffing oscillator with large nonlinearity and large damping are investigated by using a simple point collocation method applied in the time domain over a period of the periodic solution. The relationship between the proposed collocation method and the high dimensional harmonic balance method (HDHB), proposed earlier by Thomas, Dowell, and Hall (2002), is explored. We demonstrate that the HDHB is not a kind of "harmonic balance method" but essentially a cumbersome version of the collocation method. In using the collocation method,… More >

Ayşe Gül Kaplan¹, Yılmaz Dereli

CMES-Computer Modeling in Engineering & Sciences, Vol.84, No.5, pp. 423-438, 2012, DOI:10.3970/cmes.2012.084.423

Abstract In this study, the nonlinear symmetric regularized long wave equation was solved numerically by using radial basis functions collocation method. The single solitary wave solution, the interaction of two positive solitary waves and the clash of two solitary waves were studied. Numerical results and simulations of the wave motions were presented. Validity and accuracy of the method was tested by compared with results in the literature. More >