Liming Yang^1,2, Junjian Bai¹, Qun Sun³

CMES-Computer Modeling in Engineering & Sciences, Vol.108, No.1, pp. 49-65, 2015, DOI:10.3970/cmes.2015.108.049

Abstract Extreme learning machine (ELM) has demonstrated great potential in machine learning and data mining fields owing to its simplicity, rapidity and good generalization performance. In this work, a general framework for ELM regression is first investigated based on least absolute deviation (LAD) estimation (called LADELM), and then we develop two regularized LADELM formulations with the l₂-norm and l₁-norm regularization, respectively. Moreover, the proposed models are posed as simple linear programming or quadratic programming problems. Furthermore, the proposed models are used directly to analyze the hard rate of licorice seeds using near-infrared spectroscopy data. Experimental results on eight different spectral regions… More >

T. Gortsas¹, S.V. Tsinopoulos², D. Polyzos^1,3

CMES-Computer Modeling in Engineering & Sciences, Vol.107, No.4, pp. 321-343, 2015, DOI:10.3970/cmes.2015.107.321

Abstract An advanced Boundary Element method (BEM) accelerated via Adaptive Cross Approximation (ACA) and Hierarchical Matrices (HM) techniques is presented for the solution of large-scale elastostatic problems with multi-connected domains like in fiber reinforced composite materials. Although the proposed ACA/ BEM is demonstrated for two-dimensional (2D) problems, it is quite general and it can be used for 3D problems. Different forms of ACA technique are employed for exploring their efficiency when they combined with a BEM code. More precisely, the fully and partially pivoted ACA with and without recompression are utilized, while the solution of the final linear system of equations… More >

R. Q. Rodríguez^1,2, C. L. Tan², P. Sollero¹, E. L. Albuquerque³

CMES-Computer Modeling in Engineering & Sciences, Vol.102, No.5, pp. 359-372, 2014, DOI:10.3970/cmes.2014.102.359

Abstract The efficient evaluation of the fundamental solution for 3D general anisotropic elasticity is a subject of great interest in the BEM community due to its mathematical complexity. Recently, Tan, Shiah, andWang (2013) have represented the algebraically explicit form of it developed by Ting and Lee (Ting and Lee, 1997; Lee, 2003) by a computational efficient double Fourier series. The Fourier coefficients are numerically evaluated only once for a specific material and are independent of the number of field points in the BEM analysis. This work deals with the application of hierarchical matrices and low rank approximations, applying the Adaptive Cross… More >

R. Q. Rodríguez¹, A.F. Galvis¹, P. Sollero¹, E. L. Albuquerque²

CMES-Computer Modeling in Engineering & Sciences, Vol.96, No.4, pp. 259-274, 2013, DOI:10.3970/cmes.2013.096.259

Abstract Over recent years the rapid evolution of the computational power has motivated the development of new numerical techniques to account for engineering solutions. The Boundary Element Method (BEM) has shown to be a powerful numeric tool for the analysis and solution of many physical and engineering problems. However, BEM fully populated and non-symmetric system matrices implies in higher memory requirements and solution times. This work analyze the application of hierarchical matrices and low rank approximations, applying the Adaptive Cross Approximation - ACA, to multiple inclusion potential problems. The use of hierarchical format is aimed at reducing the storage requirement and… More >

P. Redou¹, L. Gaubert¹, G. Desmeulles¹, P-A. Béal², C. Le Gal², V. Rodin³

CMES-Computer Modeling in Engineering & Sciences, Vol.70, No.1, pp. 11-40, 2010, DOI:10.3970/cmes.2010.070.011

Abstract Multi Interaction Systems, used in the context of Virtual Reality, are dedicated to real-time interactive simulations. They open the way to the in virtuo experimentation, especially useful in the domain of biochemical kinetics. To this purpose, chaotic and asynchronous scheduling of autonomous processes is based upon desynchronization of phenomena involved in the system. It permits interactivity, especially the capability to add or remove phenomena in the course of a simulation. It provides methods of resolution of ordinary differential systems and partial derivative equations. Proofs of convergence for these methods have been established, but the problem of absolute stability, although it… More >

Gusein Sh. Guseinov¹

CMES-Computer Modeling in Engineering & Sciences, Vol.86, No.4, pp. 301-320, 2012, DOI:10.3970/cmes.2012.086.301

Abstract This paper deals with the inverse spectral problem for two spectra of finite order complex Jacobi matrices (tri-diagonal symmetric matrices with complex entries). The problem is to reconstruct the matrix using two sets of eigenvalues, one for the original Jacobi matrix and one for the matrix obtained by replacing the first diagonal element of the Jacobi matrix by some another number. The uniqueness and existence results for solution of the inverse problem are established and an explicit algorithm of reconstruction of the matrix from the two spectra is given. More >

Chein-Shan Liu¹, Hong-Ki Hong¹, Satya N. Atluri²

CMES-Computer Modeling in Engineering & Sciences, Vol.60, No.3, pp. 279-308, 2010, DOI:10.3970/cmes.2010.060.279

Abstract We propose novel algorithms to calculate the inverses of ill-conditioned matrices, which have broad engineering applications. The vector-form of the conjugate gradient method (CGM) is recast into a matrix-form, which is named as the matrix conjugate gradient method (MCGM). The MCGM is better than the CGM for finding the inverses of matrices. To treat the problems of inverting ill-conditioned matrices, we add a vector equation into the given matrix equation for obtaining the left-inversion of matrix (and a similar vector equation for the right-inversion) and thus we obtain an over-determined system. The resulting two modifications of the MCGM, namely the… More >

Yongchang Cai^1,2, J.K. Paik³, Satya N. Atluri³

CMES-Computer Modeling in Engineering & Sciences, Vol.54, No.3, pp. 335-368, 2009, DOI:10.3970/cmes.2009.054.335

Abstract This paper presents a simple finite element method, based on assumed moments and rotations, for geometrically nonlinear large rotation analyses of space frames consisting of members of arbitrary cross-section. A von Karman type nonlinear theory of deformation is employed in the updated Lagrangian co-rotational reference frame of each beam element, to account for bending, stretching, and torsion of each element. The Reissner variational principle is used in the updated Lagrangian co-rotational reference frame, to derive an explicit expression for the (12x12)symmetrictangent stiffness matrix of the beam element in the co-rotational reference frame. The explicit expression for the finite rotation of… More >

Yongchang Cai^1,2, J.K. Paik³, Satya N. Atluri³

CMES-Computer Modeling in Engineering & Sciences, Vol.53, No.2, pp. 123-152, 2009, DOI:10.3970/cmes.2009.053.123

Abstract This paper presents a simple finite element method, based on simple mechanics and physical clarity, for geometrically nonlinear large rotation analyses of space frames consisting of members of arbitrary cross-section. A co-rotational reference frame, involving the axes of each finitely rotated beam finite-element, is used as the Updated Lagrangian reference frame for the respective element. A von Karman type nonlinear theory of deformation is employed in the co-rotational reference frame of each beam element, to account for bending, stretching, and torsion of each element. An assumed displacement approach is used to derive an explicit expression for the (12x12)symmetrictangent stiffness matrix… 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 new polynomial interpolation… More >