S Reaz Ahmed¹, S K Deb Nath¹

CMES-Computer Modeling in Engineering & Sciences, Vol.44, No.1, pp. 35-64, 2009, DOI:10.3970/cmes.2009.044.035

Abstract The paper presents a simplified analysis of stresses and deformations at critical sections of a tire-tread. Displacement potential formulation is used in conjunction with the finite-difference method to model the present contact problem. The solution of the problem is obtained for two limiting cases of the contact boundary - one allows the lateral slippage and the other conforms to the no-slip condition along the lateral direction. The influential effects of tire material and tread aspect-ratio are discussed. The reliability and accuracy of the solution is also discussed in light of comparison made with the usual More >

Chein-Shan Liu¹, Satya N. Atluri²

CMES-Computer Modeling in Engineering & Sciences, Vol.43, No.3, pp. 253-276, 2009, DOI:10.3970/cmes.2009.043.253

Abstract Since the works of Newton and Lagrange, interpolation had been a mature technique in the numerical mathematics. Among the many interpolation methods, global or piecewise, the polynomial interpolation p(x) = a₀ + a₁x + ... + a_nxⁿ expanded by the monomials is the simplest one, which is easy to handle mathematically. For higher accuracy, one always attempts to use a higher-order polynomial as an interpolant. But, Runge gave a counterexample, demonstrating that the polynomial interpolation problem may be ill-posed. Very high-order polynomial interpolation is very hard to realize by numerical computations. In this paper we propose a… More >

A. Brancati¹, M. H. Aliabadi¹, I. Benedetti^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.43, No.2, pp. 149-172, 2009, DOI:10.3970/cmes.2009.043.149

Abstract In this paper a new Rapid Acoustic Boundary Element Method (RABEM) is presented using a Hierarchical GMRES solver for 3D acoustic problems. The Adaptive Cross Approximation is used to generate both the system matrix and the right hand side vector. The ACA is also used to evaluate the potential and the particle velocity values at selected internal points. Two different GMRES solution strategies (without preconditioner and with a block diagonal preconditioner) are developed and tested for low and high frequency problems. Implementation of different boundary conditions (i.e. Dirichlet, Neumann and mixed Robin) is also described. More >

Francisco P. M. Oliveira¹, João Manuel R. S. Tavares¹

CMES-Computer Modeling in Engineering & Sciences, Vol.43, No.1, pp. 91-110, 2009, DOI:10.3970/cmes.2009.043.091

Abstract This paper presents a new methodology to establish the best global match of objects' contours in images. The first step is the extraction of the sets of ordered points that define the objects' contours. Then, by using the curvature value and its distance to the corresponded centroid for each point, an affinity matrix is built. This matrix contains information of the cost for all possible matches between the two sets of ordered points. Then, to determine the desired one-to-one global matching, an assignment algorithm based on dynamic programming is used. This algorithm establishes the global More >

Surkay D. Akbarov^1,2,3, A. R. Mamedov³

CMES-Computer Modeling in Engineering & Sciences, Vol.42, No.3, pp. 257-296, 2009, DOI:10.3970/cmes.2009.042.257

Abstract Within the scope of the piecewise homogeneous body model through the use of the three-dimensional geometrically non-linear exact equations of the theory of elasticity, an approach for the investigation of problems with the micromechanics of a periodically curved fiber near the free convex cylindrical surface is proposed and employed. The main difficulties in finding the solution to these problems are caused by the impossibility of employing the summation theorem for cylindrical functions to satisfy the boundary conditions on the cylindrical surface. For this purpose the cosine and sine Fourier series presentation of the sought values More >

Rencheng Song^1,2, Wen Chen^2,3

CMES-Computer Modeling in Engineering & Sciences, Vol.42, No.1, pp. 59-70, 2009, DOI:10.3970/cmes.2009.042.059

Abstract The regularized meshless method (RMM) is a novel boundary-type meshless method but by now has mainly been tested successfully to the regular domain problems in reports. This note makes a further investigation on its solution of irregular domain problems. We find that the method fails to produce satisfactory results for some benchmark problems. The reason is due to the inaccurate calculation of the diagonal elements of the numerical discretization matrix in the original RMM, which have strong effect on the resulting solution accuracy. To overcome this severe drawback, this study introduces the weighted diagonal element More >

G. Haasemann¹, M. Kästner², V. Ulbricht³

CMES-Computer Modeling in Engineering & Sciences, Vol.42, No.1, pp. 35-58, 2009, DOI:10.3970/cmes.2009.042.035

Abstract The purpose of this paper is the presentation of a new efficient modelling strategy based on the combination of Binary Model and Extended Finite Element Method (X-FEM). It is applied to represent the internal architecture of textile reinforced composites where the resin-saturated fabric is characterised by a complex geometry. Homogenisation methods are used to compute the effective elastic material properties. Thereby, the discrete formulation of periodic boundary conditions is adapted regarding additional degrees of freedom used by finite elements which are based on the X-FEM. Finally, the results in terms of effective material properties reveal More >

Peter Lucas¹, Alexander H. van Zuijlen¹, Hester Bijl¹

CMES-Computer Modeling in Engineering & Sciences, Vol.41, No.2, pp. 147-176, 2009, DOI:10.3970/cmes.2009.041.147

Abstract Mesh adaptation is a fairly established tool to obtain numerically accurate solutions for flow problems. Computational efficiency is, however, not always guaranteed for the adaptation strategies found in literature. Typically excessive mesh growth diminishes the potential efficiency gain. This paper, therefore, extends the strategy proposed by [Aftosmis and Berger (2002)] to compute the refinement threshold. The extended strategy computes the refinement threshold based on a user desired number of grid cells and adaptations, thereby ensuring high computational efficiency. Because our main interest is flow around wind turbines, the adaptation strategy has been optimized for flow More >

Hung-Chang Lee¹, Chein-Shan Liu²

CMES-Computer Modeling in Engineering & Sciences, Vol.41, No.1, pp. 1-26, 2009, DOI:10.3970/cmes.2009.041.001

Abstract The group-preserving schemes developed by Liu (2001) for integrating ordinary differential equations system were adopted the Cayley transform and Padé approximants to formulate the Lie group from its Lie algebra. However, the accuracy of those schemes is not better than second-order. In order to increase the accuracy by employing the group-preserving schemes on ordinary differential equations, according to an efficient technique developed by Runge and Kutta to raise the order of accuracy from the Euler method, we combine the Runge-Kutta method on the group-preserving schemes to obtain the higher-order numerical methods of group-preserving type. They More >

S.H. Pyo¹, H.K. Lee^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.40, No.3, pp. 271-306, 2009, DOI:10.3970/cmes.2009.040.271

Abstract This paper presents a micromechanical approach for predicting the elastic and multi-level damage response of aligned and randomly oriented whisker-/ short fiber-reinforced composites. Based on a combination of Eshelby's micromechanics and the evolutionary imperfect interface approach, the effective elastic moduli of the composites are derived explicitly. The modified Eshelby's tensor for spheroidal inclusions with slightly weakened interface [Qu (1993b)] is extended in the present study to model whiskers or short fibers having mild or severe imperfect interfaces. Aligned and random orientations of spheroidal reinforcements are considered. A multi-level damage model in accordance with the Weibull's More >