Y.C. Shiah¹, C. L. Tan²

CMES-Computer Modeling in Engineering & Sciences, Vol.78, No.2, pp. 95-108, 2011, DOI:10.3970/cmes.2011.078.095

Abstract By differentiating the Green function of Ting and Lee (1997) for 3D general anisotropic elastotatics in a spherical coordinate system as an intermediate step, and then using the chain rule, derivatives of up to the second order of this fundamental solution are obtained in exact, explicit, algebraic forms. No tensors of order higher than two are present in these derivatives, thereby allowing these quantities to be numerically evaluated quite expeditiously. These derivatives are required for the computation of the internal point displacements and stresses via Somigliana's identity in BEM analysis. Some examples are presented to More >

T. Matsumoto¹, T. Yamada¹, T. Takahashi¹, C.J. Zheng², S. Harada¹

CMES-Computer Modeling in Engineering & Sciences, Vol.78, No.2, pp. 77-94, 2011, DOI:10.3970/cmes.2011.078.077

Abstract Shape design and topology sensitivity formulations for acoustic problems based on adjoint method and the boundary element method are presented and are applied to shape sensitivity analysis and topology optimization of acoustic field. The objective function is assumed to consist only of boundary integrals and quantities defined at certain number of discrete points. The adjoint field is defined so that the sensitivity of the objective function does not include the unknown sensitivity coefficients of the sound pressures and particle velocities on the boundary and in the domain. Since the final sensitivity expression does not have More >

P. Maghoul¹, B. Gatmiri^1,2, D. Duhamel¹

CMES-Computer Modeling in Engineering & Sciences, Vol.78, No.1, pp. 51-76, 2011, DOI:10.3970/cmes.2011.078.051

Abstract This paper aims at obtaining boundary integral formulations as well as three dimensional(3D) fundamental solutions for unsaturated soils under dynamic loadings for the first time. The boundary integral equations are derived via the use of the weighted residuals method in a way that permits an easy discretization and implementation in a Boundary Element code. Also, the associated 3D fundamental solutions for such deformable porous medium are derived in Laplace transform domain using the method of Hérmander. The derived results are verified analytically by comparison with the previously introduced corresponding fundamental solutions in elastodynamic limiting case. More >

A. Brancati¹, M. H. Aliabadi¹, A. Milazzo²

CMES-Computer Modeling in Engineering & Sciences, Vol.78, No.1, pp. 1-24, 2011, DOI:10.3970/cmes.2011.078.001

Abstract This paper presents an improved adaptive cross approximation (ACA) approach developed in conjunction with the Hierarchical format matrix and the GMRES solver. A novel scheme to generate the cluster tree (based upon preliminary considerations of the prescribed boundary conditions) and an improved ACA algorithm (approximating the system matrix for mixed Robin conditions) are described. The asymptotic smoothness property of a kernel generated by a linear combination of two asymptotic smooth kernels is demonstrated. Numerical results show the new approach to be up to 50% faster than the conventional ACA approach. More >

M.Souli¹, F.Erchiqui²

CMES-Computer Modeling in Engineering & Sciences, Vol.77, No.3&4, pp. 183-200, 2011, DOI:10.3970/cmes.2011.077.183

Abstract During the design process of membrane structure to resist to high pressure loading, and the characterization of hyperelastic material, a structure made up of thin rubber undergoes large deformation and rotation under high pressure loading out of high pressurized gas. Until recently, to simulate the inflation of the hyperelastic membrane, a uniform pressure based on thermodynamic model or experimental tests is applied to the structure, as boundary conditions. From a computational time point of view, this approach is very fast, since no computational fluid dynamics is involved in the simulation. However, at the late stage… More >

Minghui Wang¹, Musheng Wei², Shanrui Hu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.77, No.3&4, pp. 173-182, 2011, DOI:10.3970/cmes.2011.077.173

Abstract The mapping from the symmetric solution set to its independent parameter space is studied and an iterative method is proposed for the least-squares minimum-norm symmetric solution of AXB = E. Numerical results are reported that show the efficiency of the proposed methods. More >

Zhou YH^1,2, Wang XM², Wang JZ^1,2 , Liu XJ²

CMES-Computer Modeling in Engineering & Sciences, Vol.77, No.2, pp. 137-160, 2011, DOI:10.3970/cmes.2011.077.137

Abstract In this paper, we present an efficient wavelet-based algorithm for solving a class of fractional vibration, diffusion and wave equations with strong nonlinearities. For this purpose, we first suggest a wavelet approximation for a function defined on a bounded interval, in which expansion coefficients are just the function samplings at each nodal point. As the fractional differential equations containing strong nonlinear terms and singular integral kernels, we then use Laplace transform to convert them into the second type Voltera integral equations with non-singular kernels. Certain property of the integral kernel and the ability of explicit More >

A.B. Tatomir^1,2, A.Szymkiewicz³, H. Class¹, R. Helmig¹

CMES-Computer Modeling in Engineering & Sciences, Vol.77, No.2, pp. 81-112, 2011, DOI:10.3970/cmes.2011.077.081

Abstract We present a two phase flow conceptual model, the corresponding simulator (2pMINC) and a workflow for large-scale fractured reservoirs, based on a continuum fracture approach which uses the multiple interacting continua (MINC) method complemented with an improved upscaling technique. The complex transient behavior of the flow processes in fractured porous media is captured by subgridding the coarse blocks in nested volume elements which have effective properties calculated from the detailed representation of the fracture system. In this way, we keep a physically based approach, preserve the accuracy of the model, avoid the common use of… More >

Chein-Shan Liu¹, Chung-Lun Kuo²

CMES-Computer Modeling in Engineering & Sciences, Vol.77, No.1, pp. 57-80, 2011, DOI:10.3970/cmes.2011.077.057

Abstract In this paper, the solutions of inverse Cauchy problems for quasi-linear elliptic equations are resorted to an unusual mixed group-preserving scheme (MGPS). The bottom of a finite rectangle is imposed by overspecified boundary data, and we seek unknown data on the top side. The spring-damping regularization method (SDRM) is introduced by converting the governing equation into a new one, which includes a spring term and a damping term. The SDRM can further stabilize the inverse Cauchy problems, such that we can apply a direct numerical integration method to solve them by using the MGPS. Several More >

Sheng Yan Li¹, Bin Gang Xu^1,2, Xiao Ming Tao¹, Hong Hu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.77, No.1, pp. 33-56, 2011, DOI:10.3970/cmes.2011.077.033

Abstract Slender yarn structure made from natural fibers, nano-fibers, carbon nanotubes or other types of fibrous materials is all formed by twisting an assembly of short or long fibers and its performance is significantly influenced by the physical behavior of these fibers in the slender yarn forming region - a small triangle area called spinning triangle. In this paper, a new generalized FEM model of spinning triangle has been developed to theoretically analyze the fiber structural and mechanical performance in fabrication of these slender yarn structures. In this proposed model, a geometrical model of spinning triangle More >