L. Rodríguez-Tembleque¹, A. Sáez¹, F.C. Buroni¹

CMES-Computer Modeling in Engineering & Sciences, Vol.96, No.2, pp. 131-158, 2013, DOI:10.3970/cmes.2013.096.131

Abstract A boundary element based formulation is applied to study numerically the tribological behavior of fiber-reinforced plastics (FRP) under different frictional contact conditions, taking into account the micromechanics of FRP. Micromechanical models presented consider continuous and short fiber reinforced plastics configurations. The Boundary Element Method (BEM) with an explicit approach for fundamental solutions evaluation is considered for computing the elastic influence coefficients. Signorini’s contact conditions and an orthotropic law of friction on the potential contact zone are enforced by contact operators over the augmented Lagrangian. The proposed methodology is applied to study carbon FRP under frictional contact. The obtained numerical results… More >

W.L.A. Pereira¹, V.J. Karam², J.A.M. Carrer³, C.S.G. Monteiro¹, W.J. Mansur¹

CMES-Computer Modeling in Engineering & Sciences, Vol.96, No.2, pp. 117-130, 2013, DOI:10.3970/cmes.2013.096.117

Abstract In this work, the BEM is applied to obtain the numerical solutions for free vibration analysis of thick square plates with two edges simply supported or clamped, and the other two edges free. A formulation based on Reissner’s theory is used here, which includes the contribution of the additional translational inertia terms to the integral equation of displacements and internal forces. The boundary element method is used to discretize the space, where it is employed the static fundamental solution. In literature, the responses for the kind of problem addressed here are very important in the hydroelastic analysis of very large… More >

A. F. Galvis¹, R. Q. Rodriguez¹, P. Sollero¹, E. L. Alburquerque²

CMES-Computer Modeling in Engineering & Sciences, Vol.96, No.2, pp. 103-115, 2013, DOI:10.3970/cmes.2013.096.103

Abstract Polycrystalline structures are present on metal alloys. Therefore, it is necessary to understand and model the mechanical behavior of this media. Usually, this is accomplished by the use of different numerical methods. However, the analysis of polycrystalline materials leads to other type of problems, such as high computational requirements generated in order to get an efficient solution. In this work, the 2D polycrystalline structure is generated using an average grain size through the Voronoi tessellation method and discretized through simulations with random material, crystalline orientation and orthotropic behavior [Sfantos and Aliabadi (2007a)]. BEM discretization requires multidomain analysis and large-scale degrees… More >

G.R. Baker¹

CMES-Computer Modeling in Engineering & Sciences, Vol.96, No.2, pp. 91-101, 2013, DOI:10.3970/cmes.2013.096.091

Abstract A vortex patch is a bounded region of uniform vorticity in twodimensional, incompressible, inviscid fluid flow. The streamfunction satisfies the Poisson equation with the vorticity acting as a source term. The standard formulation is to write the streamfunction as a convolution of the vorticity with the twodimensional free-space Greens function. A simple application of Greens theorem converts the area integral to a boundary integral. Numerical methods must then account for the singular nature of the boundary integral, and high accuracy is difficult when filamentation takes place, that is, when long, very thin filaments of vorticity erupt from the main boundary.… More >

N. Thai-Quang¹, N. Mai-Duy¹, C.-D. Tran¹, T. Tran-Cong^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.96, No.1, pp. 49-90, 2013, DOI:10.3970/cmes.2013.096.049

Abstract In this paper, we present a numerical scheme, based on the direct forcing immersed boundary (DFIB) approach and compact integrated radial basis function (CIRBF) approximations, for solving the Navier-Stokes equations in two dimensions. The problem domain of complicated shape is embedded in a Cartesian grid containing Eulerian nodes. Non-slip conditions on the inner boundaries, represented by Lagrangian nodes, are imposed by means of the DFIB method, in which a smoothed version of the discrete delta functions is utilised to transfer the physical quantities between two types of nodes. The velocities and pressure variables are approximated locally on Eulerian nodes using… More >

Yunpeng Ma¹, Lifeng Wang¹, Zhijun Meng¹

CMES-Computer Modeling in Engineering & Sciences, Vol.96, No.1, pp. 31-47, 2013, DOI:10.3970/cmes.2013.096.031

Abstract In this paper, we propose a numerical algorithm for solving linear and nonlinear fractional integro-differential equations based on our constructed fractional order generalized block pulse functions operational matrix of integration. The linear and nonlinear fractional integro-differential equations are transformed into a system of algebraic equations by the matrix and these algebraic equations are solved through known computational methods. Further some numerical examples are shown to illustrate the accuracy and reliability of the proposed approach. Moreover, comparing the methodology with the known technique shows that our approach is more efficient and more convenient. More >

E. Tombarević¹, V. R. Voller², I. Vušanović¹

CMES-Computer Modeling in Engineering & Sciences, Vol.96, No.1, pp. 1-29, 2013, DOI:10.3970/cmes.2013.096.001

Abstract The Control Volume Finite Element Method (CVFEM) combines the geometric flexibility of the Finite Element Method (FEM) with the physical intuition of the Control Volume Method (CVM). These two features of the CVFEM make it a very powerful tool for solving heat and fluid flow problems within complex domain geometries. In solving problems in the two-dimensional domains the development of the CVFEM has been well documented. For the three-dimensional problems, while there is extensive reporting on the details of the numerical approximation, there is relatively sparse information on important issues related to data structure and interpolation. Here, in the context… More >

Changqing Miao¹, Yintao Wei², Xiangqiao Yan¹

CMES-Computer Modeling in Engineering & Sciences, Vol.95, No.6, pp. 519-534, 2013, DOI:10.3970/cmes.2013.095.519

Abstract By using a hybrid displacement discontinuity method, the interactions of three parallel square-hole cracks in an infinite plate subjected to internal pressure are investigated in this paper. Numerical examples are included to illustrate that the numerical approach is very simple and effective for calculating the stress intensity factors (SIFs) of complex plane crack problems. Many numerical results of the SIFs are given and discussed. It is found that a square hole has a shielding effect on crack(s) emanating from the hole. The finding perhaps has an important meaning in engineering. More >

Yuping Zhu^1,2, Tao Shi¹, Yuanbing Wang¹

CMES-Computer Modeling in Engineering & Sciences, Vol.95, No.6, pp. 501-517, 2013, DOI:10.3970/cmes.2013.095.501

Abstract Single crystal ferromagnetic shape memory alloy is a kind of new intelligent materials, it shows obvious anisotropy. Micromechanics theory has been used to analyze the whole mechanical behaviors of this material. However, Eshelby’s tensor of this material which plays an important role has still not solved efficiently. Based on the existing micromechanics constitutive model, this paper analyzes the numerical calculation formula of Eshelby’s tensor of anisotropic ferromagnetic shape memory alloy. Adopting the way of Gauss integral, the optimal Gaussian integral points for different inclusion shapes and the corresponding numerical solution of Eshelby’s tensor are obtained.Furthermore, the influence of inclusion shapes… More >

A. Balabel^1,2,3, W. A. El-askary², S. Wilson²

CMES-Computer Modeling in Engineering & Sciences, Vol.95, No.6, pp. 483-499, 2013, DOI:10.3970/cmes.2013.095.483

Abstract In the present paper, the impingement of air jet on a heated flat plate subjected to a periodic oscillation is numerically and experimentally investigated. The motivation of the present research is the desire to enhance the heat transfer characteristics during the cooling process of a heated flat plate which can be found in many relevance industrial applications. In order to improve the heat transfer characteristics, a novel idea is utilized, where a periodical oscillation movement in form of sine wave produced from a Scotch yoke mechanism is applied to the heated flat plate. The obtained numerical results showed a good… More >