A. Gakwaya¹, H. Sharifi², M. Guillot¹, M. Souli³, F. Erchiqui⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.73, No.3, pp. 209-266, 2011, DOI:10.3970/cmes.2011.073.209

Abstract A mechanical computer aided design and engineering system can be used to reduce the design-to-manufacture cycle time in metal forming process. Such a system could be built upon a solid modeling geometry engine and an efficient finite element (FE) solver. The maintenance of a high-quality mesh throughout the analysis is an essential feature of an efficient finite element simulation of large strain metal forming problems. In this paper, a mesh adaptation technique employing the Arbitrary Lagrangian-Eulerian formulation (ALE) is applied to some industrial metal forming problems. An ACIS boundary representation of the solid model is… More >

Cheng-Te Chi¹, Ming-Hsiao Lee², Wen-Hwa Chen^1,3

CMC-Computers, Materials & Continua, Vol.24, No.3, pp. 239-256, 2011, DOI:10.3970/cmc.2011.024.239

Abstract A novel three-dimensional adaptive element-free Galerkin method (EFGM) based on a uniform background grid is proposed to cope with the problems with extremely large deformation. On the basis of this uniform background grid, an interior adaptive strategy through an error estimation within the analysis domain is developed. By this interior adaptive scheme, additional adaptive nodes are inserted in those regions where the solution accuracy needs to be improved. As opposed to the fixed uniform background grid, these inserted nodes can move along with deformation to describe the particular local deformation of the structure. In addition,… More >

M. Al-Ajmi¹, A. Benjeddou²

CMC-Computers, Materials & Continua, Vol.23, No.3, pp. 265-286, 2011, DOI:10.3970/cmc.2011.023.265

Abstract A new discrete layer finite element (DLFE) is presented for electro-mechanically coupled analyses of moderately thick piezoelectric adaptive composite plates. The retained kinematics is based on layer-wise first-order shear deformation theory, and considers the plies top and bottom surfaces in-plane displacements and the plate transverse deflection as mechanical unknowns. The former are assumed in-plane Lagrange linear, while the latter is assumed in-plane full (Lagrange) quadratic; this results in a nine nodes quadrangular (Q9) DLFE. The latter is validated in free-vibrations, first numerically against ANSYS three-dimensional piezoelectric finite elements for a cantilever moderately thick aluminum plate… More >

Yihua Xiao¹, Xu Han^1,2, Dean Hu¹

CMC-Computers, Materials & Continua, Vol.23, No.1, pp. 9-34, 2011, DOI:10.3970/cmc.2011.023.009

Abstract For impact computations, it is efficient to model small and large deformation regions by Finite Element Method (FEM) and Smoothed Particle Hydrodynamics (SPH), respectively. However, it requires an effective algorithm to couple FEM and SPH calculations. To fulfill this requirement, an alternative coupling algorithm is presented in this paper. In the algorithm, the coupling between element and particle regions are achieved by treating elements as imaginary particles and applying equivalent tractions to element sides on coupling interfaces. In addition, an adaptive coupling technique is proposed based on the algorithm to improve the computational efficiency of More >

Wenqing Wang¹, Thomas Schnicke², Olaf Kolditz³

CMC-Computers, Materials & Continua, Vol.21, No.3, pp. 217-236, 2011, DOI:10.3970/cmc.2011.021.217

Abstract This work focuses on parallel finite element simulation of thermal hydraulic and mechanical (THM) coupled processes in porous media, which is a common phenomenon in geological applications such as nuclear waste repository and CO2 storage facilities. The Galerkin finite element method is applied to solve the derived partial differential equations. To deal with the coupling terms among the equations, the momentum equation is solved individually in a monolithic manner, and moreover their solving processes are incorporated into the solving processes of nonisothermal hydraulic equation and heat transport equation in a staggered manner. The computation task… More >

Chein-Shan Liu¹

CMC-Computers, Materials & Continua, Vol.21, No.1, pp. 17-40, 2011, DOI:10.3970/cmc.2011.021.017

Abstract Only the left-boundary data of temperature and heat flux are used to estimate an unknown parameter function α(x) in T_t(x,t) = ∂(α(x)T_x)/∂x + h(x,t), as well as to recover the right-boundary data. When α(x) is given the above problem is a well-known inverse heat conduction problem (IHCP). This paper solves a mixed-type inverse problem as a combination of the IHCP and the problem of parameter identification, without needing to assume a function form of α(x) a priori, and without measuring extra data as those used by other methods. We use the one-step Lie-Group Adaptive Method (LGAM) More >

Chein-Shan Liu¹, Satya N. Atluri²

CMES-Computer Modeling in Engineering & Sciences, Vol.64, No.3, pp. 299-326, 2010, DOI:10.3970/cmes.2010.064.299

Abstract We solve the Calderón inverse conductivity problem [Calderón (1980, 2006)], for an elliptic type equation in a rectangular plane domain, to recover an unknown conductivity function inside the domain, from the over-specified Cauchy data on the bottom of the rectangle. The Calderón inverse problem exhibitsthree-fold simultaneous difficulties: ill-posedness of the inverse Cauchy problem, ill-posedness of the parameter identification, and no information inside the domain being available on the impedance function. In order to solve this problem, we discretize the whole domain into many sub-domains of finite strips, each with a small height. Thus the Calderón… More >

M. S. Bapat¹, Y. J. Liu^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.58, No.2, pp. 161-184, 2010, DOI:10.3970/cmes.2010.058.161

Abstract A new definition of the interaction list in the fast multipole method (FMM) is introduced in this paper, which can reduce the moment-to-local (M2L) translations by about 30-40% and therefore improve the efficiency for the FMM. In addition, an adaptive tree structure is investigated, which is potentially more efficient than the oct-tree structure for thin and slender domains as in the case of micro-electro-mechanical systems (MEMS). A combination of the modified interaction list (termed L2 modification in the adaptive fast multipole BEM) and the adaptive tree structure in the fast multipole BEM has been implemented… More >

D. Brunner¹, M. Junge¹, P. Rapp¹, M. Bebendorf², L. Gaul¹

CMES-Computer Modeling in Engineering & Sciences, Vol.58, No.2, pp. 131-160, 2010, DOI:10.3970/cmes.2010.058.131

Abstract The simulation of the hydroacoustic sound radiation of ship-like structures has an ever-growing importance due to legal regulations. Using the boundary element method, the overall dimension of the problem is reduced and only integrals over surfaces have to be considered. Additionally, the Sommerfeld radiation condition is automatically satisfied by proper choice of the fundamental solution. However, the resulting matrices are fully populated and the set-up time and memory consumption scale quadratically with respect to the degrees of freedom. Different fast boundary element methods have been introduced for the Helmholtz equation, resulting in a quasilinear complexity.… More >

F. S. Lobato¹, V. Steffen Jr², A. J. Silva Neto³

CMES-Computer Modeling in Engineering & Sciences, Vol.69, No.1, pp. 1-18, 2010, DOI:10.3970/cmes.2010.069.001

Abstract Differential Evolution Algorithm (DE) has shown to be a powerful evolutionary algorithm for global optimization in a variety of real world problems. DE differs from other evolutionary algorithms in the mutation and recombination phases. Unlike some other meta-heuristic techniques such as genetic algorithms and evolutionary strategies, where perturbation occurs in accordance with a random quantity, DE uses weighted differences between solution vectors to perturb the population. Although the efficiency of DE algorithm has been proven in the literature, studies indicate that the efficiency of the DE methods is sensitive to its control parameters (perturbation rate More >