X.Y. Cui^1,2, G. R. Liu^2,3, G. Y. Li^1,4, G. Zheng¹

CMES-Computer Modeling in Engineering & Sciences, Vol.38, No.3, pp. 217-230, 2008, DOI:10.3970/cmes.2008.038.217

Abstract In this paper, a gradient smoothed formulation is proposed to deal with a fourth-order differential equation of Bernoulli-Euler beam problems for static and dynamic analysis. Through the smoothing operation, the C¹ continuity requirement for fourth-order boundary value and initial value problems can be easily relaxed, and C⁰ interpolating function can be employed to solve C¹ problems. In present thin beam problems, linear shape functions are employed to approximate the displacement field, and smoothing domains are further formed for computing the smoothed curvature and bending moment field. Numerical examples indicate that very accurate results can be yielded when More >

Briti Sundar Deb¹, Srinivasan Natesan²

CMES-Computer Modeling in Engineering & Sciences, Vol.38, No.2, pp. 179-200, 2008, DOI:10.3970/cmes.2008.038.179

Abstract This paper presents an almost second--order uniformly convergent Richardson extrapolation method for convection- dominated coupled system of boundary value problems. First, we solve the system by using the classical finite difference scheme on the layer resolving Shishkin mesh, and then we construct the Richardson approximation solution using the solutions obtained on N and 2N mesh intervals. Second-order parameter--uniform error estimate is derived. The proposed method is applied to a test example for verification of the theoretical results for the case ε ≤ N⁻¹. More >

C. Zheng¹, S. C. Wu^2,3,4, X.H.Tang¹, J. H. Zhang¹

CMES-Computer Modeling in Engineering & Sciences, Vol.38, No.2, pp. 149-178, 2008, DOI:10.3970/cmes.2008.038.149

Abstract A novel meshfree poly-cell Galerkin method is developed for problems of elasticity and fracture. To improve accuracy, a poly-cell support is proposed to ensure the alignment of shape function support and the integration domain. By orthonormalizing basis functions, the improved moving least-square is formulated soundly, in which frequent matrix inversions are avoided. The Nitsche's method is introduced to treat the essential boundary conditions. It is found that computed solutions are more accurate than those obtained using the circle support used in standard MLS. Furthermore, numerical results present the superconvergent property, compared with the theoretical values More >

Sirajul Haq¹, Siraj-ul-Islam², Arshed Ali³

CMES-Computer Modeling in Engineering & Sciences, Vol.38, No.1, pp. 1-24, 2008, DOI:10.3970/cmes.2008.038.001

Abstract In this paper we propose a meshfree technique for the numerical solution of the modified equal width wave (MEW) equation. Combination of collocation method using the radial basis functions (RBFs) with first order accurate forward difference approximation is employed for obtaining meshfree solution of the problem. Different types of RBFs are used for this purpose. Performance of the proposed method is successfully tested in terms of various error norms. In the case of non-availability of exact solution, performance of the new method is compared with the results obtained from the existing methods. Propagation of a More >

D. H. Pahr¹, H.J. Böhm¹

CMES-Computer Modeling in Engineering & Sciences, Vol.34, No.2, pp. 117-136, 2008, DOI:10.3970/cmes.2008.034.117

Abstract The combination of heterogeneous volume elements and numerical analysis schemes such as the Finite Element method provides a powerful and well proven tool for studying the mechanical behavior of composite materials. Periodicity boundary conditions (PBC), homogeneous displacement boundary conditions (KUBC) and homogeneous traction boundary conditions (SUBC) have been widely used in such studies. Recently Pahr and Zysset (2008) proposed a special set of mixed uniform boundary conditions (MUBC) for evaluating the macroscopic elasticity tensor of human trabecular bone. These boundary conditions are not restricted to periodic phase geometries, but were found to give the same… More >

Song-Jeng Huang^1,2, Lin-Wei Chiu²

CMES-Computer Modeling in Engineering & Sciences, Vol.34, No.1, pp. 1-32, 2008, DOI:10.3970/cmes.2008.034.001

Abstract The composite sandwich plate is one of the most common composite structures. Local stress concentrations can be caused by localized bending effects where a load is introduced. Although a sandwich structure with an insert is one of the classical load bearing structures, little work has been conducted on the adhesive layers or inserts. This study involves a linear elasticity analysis of five-layer sandwich plates with ``through-the-thickness'' inserts, using sandwich plate theory to analyze deformation behavior. Governing equations are formulated as partial differential equations, which are solved numerically using the multi-segment integration method. Sandwich plates with More >

P. Bettini, E. Brusa, M. Munteanu, R. Specogna, F. Trevisan¹

CMES-Computer Modeling in Engineering & Sciences, Vol.33, No.3, pp. 215-242, 2008, DOI:10.3970/cmes.2008.033.215

Abstract Electrostatic microactuator is a paradigm of MEMS. Cantilever and double clamped microbeams are often used in microswitches, microresonators and varactors. An efficient numerical prediction of their mechanical behaviour is affected by the nonlinearity of the electromechanical coupling. Sometimes an additional nonlinearity is due to the large displacement or to the axial-flexural coupling exhibited in bending. To overcome the computational limits of the available numerical methods two new formulations are here proposed and compared. Modifying the classical beam element in the Finite Element Method to allow the implementation of a \emph {Non incremental sequential approach} is… More >

Ma Tianbao¹, Wang Cheng, Ning Jianguo

CMES-Computer Modeling in Engineering & Sciences, Vol.33, No.2, pp. 155-178, 2008, DOI:10.3970/cmes.2008.033.155

Abstract A multi-material Eulerian hydrodynamic numerical method and hydrocode that can effectively simulate explosion problems in engineering practice were developed in this study. A modified Youngs' interface reconstruction algorithm was proposed for mixed cells, in which the material's volume fractions of the surrounding cells are not only used to reconstruct the material interface but also adopted to determine the transport order of the material. The algorithm developed herein was validated by the modeling of several tests, such as objects with different shapes moving in translational, rotating and shear flow field in two dimensional Descartes coordinates and More >

K. Lee¹, H. Park², S.W. Lee³

CMES-Computer Modeling in Engineering & Sciences, Vol.32, No.3, pp. 113-122, 2008, DOI:10.3970/cmes.2008.032.113

Abstract A post-processing scheme is presented to accurately determine transverse shear and normal stresses in composite panels undergoing geometrically nonlinear deformation under dynamic loading conditions. Transverse stresses are assumed through thickness at a point of interest and are recovered via a one-dimensional finite element method. The finite element method is based on the least square functional of the error in the equilibrium equation along the thickness direction and utilizes the in-plane stresses and resultant transverse shear forces per unit length obtained by a shell element analysis. Numerical results demonstrate that, with minimal addition of computational efforts, More >

Shishir Sinha¹, Praveen Ch.

CMES-Computer Modeling in Engineering & Sciences, Vol.32, No.2, pp. 85-102, 2008, DOI:10.3970/cmes.2008.032.085

Abstract Fluid catalytic cracking (FCC) unit plays most important role in the economy of a modern refinery that it is use for value addition to the refinery products. Because of the importance of FCC unit in refining, considerable effort has been done on the modeling of this unit for better understanding and improved productivity. The process is characterized by complex interactions among feed quality, catalyst properties, unit hardware parameters and process conditions. \newline The traditional and global approach of cracking kinetics is lumping. In the present paper, five lump kinetic scheme is considered, where gas oil… More >