Chein-Shan Liu¹, Su-Ying Zhang², Satya N. Atluri³

CMES-Computer Modeling in Engineering & Sciences, Vol.88, No.1, pp. 29-48, 2012, DOI:10.3970/cmes.2012.088.029

Abstract With a detailed investigation of n linear algebraic equations Bx=b, we find that the scaled residual dynamics for y∈Sⁿ⁻¹ is equipped with four structures: the Jordan dynamics, the rotation group SO(n), a generalized Hamiltonian formulation, as well as a metric bracket system. Therefore, it is the first time that we can compute the steplength used in the iterative method by a novel algorithm based on the Jordan structure. The algorithms preserving the length of y are developed as the structure preserving algorithms (SPAs), which can significantly accelerate the convergence speed and are robust enough against the noise in the numerical… More >

A.E. Martínez-Castro¹, I.H. Faris¹, R. Gallego¹

CMES-Computer Modeling in Engineering & Sciences, Vol.87, No.3, pp. 177-206, 2012, DOI:10.3970/cmes.2012.087.177

Abstract In this paper, the identification of hidden defects inside a three-dimen -sional layer is set as an Identification Inverse Problem. This problem is solved by minimizing a cost functional which is linearized with respect to the volume defects, leading to a procedure that requires only computations at the host domain free of defects. The cost functional is stated as the misfit between experimental and computed displacements and spherical and/or ellipsoidal cavities are the defects to locate. The identification of these cavities is based on the measured displacements at a set of points due to time-harmonic point loads at an array… More >

Peter J. Ireman, Anders Klarbring¹, Niclas Strömberg

CMES-Computer Modeling in Engineering & Sciences, Vol.52, No.2, pp. 125-158, 2009, DOI:10.3970/cmes.2009.052.125

Abstract In this paper finite element approaches for fretting fatigue are proposed on the basis of a non-local model of continuum damage coupled to friction and wear. The model is formulated in the frame-work of a standard material. In a previous paper this was done in the spirit of Maugin, where an extra entropy flux is introduced in the second law in order to include the gradient of the internal variable in a proper manner. In this paper we follow instead the ideas of Frémond and others, where this extra entropy flux is no longer needed, but instead new non-classical balance… More >

C.E. Majorana¹, V.A. Salomoni

CMES-Computer Modeling in Engineering & Sciences, Vol.33, No.1, pp. 49-84, 2008, DOI:10.3970/cmes.2008.033.049

Abstract A dynamic analysis of a thin shell finite element undergoing large displacements and rotations is here presented. The constitutive model adopted derives from the coupling of an hyperelastic basic model fulfilling a De Saint Venant-Kirchhoff criterion with a scalar damage function depending on the maximum value of a suitable strain measure attained through the deformation history; then plastic effects are included using an isotropic/kinematic hardening law. A conservative time integration scheme for the non-linear dynamics of the hyperelastic damaged-plastic thin shell is applied. The main characteristic of the scheme is to be conservative, since it allows for the time-discrete system… More >

E. Bleszynski¹, M. Bleszynski¹, T. Jaroszewicz¹

CMES-Computer Modeling in Engineering & Sciences, Vol.5, No.4, pp. 295-318, 2004, DOI:10.3970/cmes.2004.005.295

Abstract We describe elements of our current work on the development of new methods for high frequency electromagnetic scattering, based on the wavefront (WF) representation of propagating fields and on the asymptotic but rigorous solution of integral equations for surface currents. In the wavefront evolution technique, surfaces of constant phase are constructed and treated not merely as collections of independent rays, but as well defined geometrical objects endowed with the full connectivity information. Hence, a precise determination of shadow and reflection boundaries, a construction of (multiply) diffracted wavefronts, a dynamic adjustment of the number of rays, an approximately constant ray-ray distance,… More >

W.J. Santos¹ , J.A.F. Santiago¹, J.C.F Telles¹

CMES-Computer Modeling in Engineering & Sciences, Vol.87, No.1, pp. 23-40, 2012, DOI:10.3970/cmes.2012.087.023

Abstract The aim of this paper is to present numerical simulations of Cathodic Protection (CP) Systems using a Genetic Algorithm (GA) and the Method of Fundamental Solutions (MFS). MFS is used to obtain the solution of the associated homogeneous equation with the non-homogeneous equation subject to nonlinear boundary conditions defined as polarization curves. The adopted GA minimizes a nonlinear error function, whose design variables are the coefficients of the linear superposition of fundamental solutions and the positions of the source points, located outside the problem domain. In this work, the anodes added to the CP system are considered as point sources… 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 >

S.R. Santos¹, L.C. Matioli², A.T. Beck³

CMES-Computer Modeling in Engineering & Sciences, Vol.83, No.1, pp. 23-56, 2012, DOI:10.3970/cmes.2012.083.023

Abstract Solution of structural reliability problems by the First Order method require optimization algorithms to find the smallest distance between a limit state function and the origin of standard Gaussian space. The Hassofer-Lind-Rackwitz-Fiessler (HLRF) algorithm, developed specifically for this purpose, has been shown to be efficient but not robust, as it fails to converge for a significant number of problems. On the other hand, recent developments in general (augmented Lagrangian) optimization techniques have not been tested in aplication to structural reliability problems. In the present article, three new optimization algorithms for structural reliability analysis are presented. One algorithm is based on… More >

Marco Dondero¹, Adrián P. Cisilino^1,2, J. Pablo Tomba¹

CMES-Computer Modeling in Engineering & Sciences, Vol.78, No.3&4, pp. 225-246, 2011, DOI:10.3970/cmes.2011.078.225

Abstract This paper introduces a numerical methodology for the design of random micro-heterogeneous materials with functionally graded effective thermal conductivities (ETC). The optimization is carried out using representative volume elements (RVEs), a parallel Genetic Algorithm (GA) as optimization method, and a Fast Multipole Boundary Element Method (FMBEM) for the evaluation of the cost function. The methodology is applied for the design of foam-like microstructures consisting of random distributions of circular insulated holes. The temperature field along a material sample is used as objective function, while the spatial distribution of the holes is the design variable. There are presented details of the… More >

G.A. Gravvanis¹, K.M. Giannoutakis²

CMES-Computer Modeling in Engineering & Sciences, Vol.71, No.4, pp. 305-330, 2011, DOI:10.3970/cmes.2011.071.305

Abstract In this paper we present parallel explicit approximate inverse matrix techniques for solving sparse linear systems on shared memory systems, which are derived using the finite element method for biharmonic equations in three space variables. Our approach for solving such equations is by considering the biharmonic equation as a coupled equation approach (pair of Poisson equation), using a FE approximation scheme, yielding an inner-outer iteration method. Additionally, parallel approximate inverse matrix algorithms are introduced for the efficient solution of sparse linear systems, based on an anti-diagonal computational approach that eliminates the data dependencies. Parallel explicit preconditioned conjugate gradient-type schemes in… More >