Chein-Shan Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.39, No.2, pp. 125-154, 2009, DOI:10.3970/cmes.2009.039.125

Abstract We propose a new numerical method for solving the boundary value problems of ordinary differential equations (ODEs) under multipoint boundary conditions specified at t = T_i, i = 1,...,m, where T₁ < ... < T_m. The finite difference scheme is used to approximate the ODEs, which together with the m-point boundary conditions constitute a system of nonlinear algebraic equations (NAEs). Then a Fictitious Time Integration Method (FTIM) is used to solve these NAEs. Numerical examples confirm that the new approach is highly accurate and efficient with a fast convergence. The FTIM can also be used to find the periods of… More >

J. P. Xu¹, Y. Q. Liu^1,2, S. Q. Li³, S. C. Wu⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.38, No.3, pp. 201-216, 2008, DOI:10.3970/cmes.2008.038.201

Abstract For the first time, a special-purpose finite element simulation system CoinForm is developed to analyze the embossing process of commemorative coin, in which one-point reduced integration approach is used in solid element finite element dynamic explicit program. Viscous damping hourglass control algorithm can effectively suppress the spurious modes activated by reduced integration and the computational effort is saved about 93% compared with other method that evaluate anti-hourglass force using stabilization matrix. The embossing process of commemorative coin is then simulated and compared with results from the DEFORM 3D software, which verify the excellent performance of present CoinForm system. According to… 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 ``through-the-thickness'' inserts subjected to axisymmetric… More >

Chein-Shan Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.33, No.2, pp. 179-198, 2008, DOI:10.3970/cmes.2008.033.179

Abstract Dirichlet boundary value problem of quasilinear elliptic equation is numerically solved by using a new concept of fictitious time integration method (FTIM). We introduce a fictitious time coordinate t by transforming the dependent variable u(x,y) into a new one by (1+t)u(x,y) =: v(x,y,t), such that the original equation is naturally and mathematically equivalently written as a quasilinear parabolic equation, including a viscous damping coefficient to enhance stability in the numerical integration of spatially semi-discretized equation as an ordinary differential equations set on grid points. Six examples of Laplace, Poisson, reaction diffusion, Helmholtz, the minimal surface, as well as the explosion… More >

Chein-Shan Liu¹, Satya N. Atluri²

CMES-Computer Modeling in Engineering & Sciences, Vol.31, No.2, pp. 71-84, 2008, DOI:10.3970/cmes.2008.031.071

Abstract Iterative algorithms for solving a nonlinear system of algebraic equations of the type: F_i(x_j) = 0, i,j = 1,…,n date back to the seminal work of Issac Newton. Nowadays a Newton-like algorithm is still the most popular one due to its easy numerical implementation. However, this type of algorithm is sensitive to the initial guess of the solution and is expensive in the computations of the Jacobian matrix ∂ F_i/ ∂ x_j and its inverse at each iterative step. In a time-integration of a system of nonlinear Ordinary Differential Equations (ODEs) of the type B_ijx_j + F_i = 0 where… More >

E. L. Albuquerque¹, M. H. Aliabadi²

CMES-Computer Modeling in Engineering & Sciences, Vol.29, No.2, pp. 63-74, 2008, DOI:10.3970/cmes.2008.029.063

Abstract This paper presents a boundary element formulation for the analysis of thin shallow shells. Classical plate bending and plane elasticity formulations are coupled and effects of curvature are treated as body forces. The body forces are written as a sum of approximation functions multiplied by coefficients. Domain integrals that arise in the formulation are transformed into boundary integrals by the radial integration method. Two different approximation functions are employed, that is 1 + r and r² log r. The method is applied to several problems and the accuracy of each approximation function is assessed by comparison with results from literature. More >

Zdenko Tonković¹, Jurica Sorić^1,2, Ivica Skozrit¹

CMES-Computer Modeling in Engineering & Sciences, Vol.26, No.2, pp. 75-90, 2008, DOI:10.3970/cmes.2008.026.075

Abstract An efficient numerical algorithm for modeling of cyclic elastoplastic deformation of shell structures is derived. The constitutive model includes highly nonlinear multi-component forms of kinematic and isotropic hardening functions in conjunction with von Mises yield criterion. Therein, the closest point projection algorithm employing the Reissner-Mindlin type kinematic model, completely formulated in tensor notation, is applied. A consistent elastoplastic tangent modulus ensures high convergence rates in the global iteration approach. The integration algorithm has been implemented into a layered assumed strain isoparametric finite shell element, which is capable of geometrical nonlinearities including finite rotations. Numerical examples, considering the symmetric and nonsymmetric… More >

Z. D. Han, S. N. Atluri¹

CMES-Computer Modeling in Engineering & Sciences, Vol.21, No.1, pp. 41-52, 2007, DOI:10.3970/cmes.2007.021.041

Abstract Straight-forward systematic derivations of the weakly singular boundary integral equations (BIEs) are presented, following the simple and directly-derived approach by Okada, Rajiyah, and Atluri (1989b) and Han and Atluri (2002). A set of weak-forms and their algebraic combinations have been used to avoid the hyper-singularities, by directly applying the "intrinsic properties'' of the fundamental solutions. The systematic decomposition of the kernel functions of BIEs is presented for regularizing the BIEs. The present approach is general, and is applied to developing weakly-singular BIEs for solids and acoustics successfully. More >

Yan Liu¹, Xiong Zhang¹, K. Y. Sze², Min Wang¹

CMES-Computer Modeling in Engineering & Sciences, Vol.20, No.3, pp. 177-192, 2007, DOI:10.3970/cmes.2007.020.177

Abstract In molecular simulations, the frequencies of the low-frequency modes are many orders of magnitude lower than those of the high-frequency modes. Compared with the amplitudes of the low-frequency modes, the amplitudes of the high-frequency modes are often negligible and, thus, least interesting. As dictated by the period of the highest frequency mode, the critical time step for stable time integration can be significantly increased by suppressing the negligible high-frequency modes yet the solution remains virtually intact. In this light, a smoothed molecular dynamics (SMD) approach is proposed to eliminate the high-frequency modes from the dynamical system through the use of… More >

Tarun Kant¹, Sandeep S. Pendhari², Yogesh M. Desai³

CMES-Computer Modeling in Engineering & Sciences, Vol.17, No.2, pp. 135-162, 2007, DOI:10.3970/cmes.2007.017.135

Abstract A two-point boundary value problem (BVP) is formed in the present work governed by a set of first-order coupled ordinary differential equations (ODEs) in terms of displacements and the transverse stresses through the thickness of laminate (in domain -h/2 < z < h/2) by introducing partial discretization methodology only in the plan area of the three dimensional (3D) laminate. The primary dependent variables in the ODEs are those which occur naturally on a plane z=a constant. An effective numerical integration (NI) technique is utilized for tackling the two-point BVP in an efficient manner. Numerical studies on cross-ply and angle-ply composite… More >