Chein-Shan Liu¹

CMC-Computers, Materials & Continua, Vol.8, No.2, pp. 53-66, 2008, DOI:10.3970/cmc.2008.008.053

Abstract Proposed is a time-marching algorithm to solve a nonlinear system of complementarity equations: P_i(x_j) ≥ 0, Q_i(x_j) ≥ 0 , P_i(x_j)Q_i(x_j) = 0, i, j = 1,...,n, resulting from a discretization of nonlinear obstacle problem. We transform the above nonlinear complementarity problem (NCP) into a nonlinear algebraic equations (NAEs) system: F_i(x_j) = 0 with the aid of the Fischer-Burmeister NCP-function. Such NAEs are semi-smooth, highly nonlinear and usually implicit, being hard to handle by the Newton-like method. Instead of, a first-order system of ODEs is derived through a fictitious time equation. The time-stepping equations are obtained by applying a numerical… More >

Muhammad Shoaib Arif¹, Mairaj Bibi², Adnan Jhangir³

CMC-Computers, Materials & Continua, Vol.54, No.2, pp. 181-195, 2018, DOI:10.3970/cmc.2018.054.181

Abstract This communique is opted to study the approximate solution of the Algebraic Lyapunov equation on the manifold of positive-definite Hermitian matrices. We choose the geodesic distance between -AHX - XA and P as the cost function, and put forward the Extended Hamiltonian algorithm (EHA) and Natural gradient algorithm (NGA) for the solution. Finally, several numerical experiments give you an idea about the effectiveness of the proposed algorithms. We also show the comparison between these two algorithms EHA and NGA. Obtained results are provided and analyzed graphically. We also conclude that the extended Hamiltonian algorithm has better convergence speed than the… More >

Chein-Shan Liu¹, Satya N. Atluri²

CMC-Computers, Materials & Continua, Vol.41, No.1, pp. 1-36, 2014, DOI:10.3970/cmc.2014.041.001

Abstract In this paper, we adopt the viewpoint of a nonlinear complementarity problem (NCP) to derive an index-one differential algebraic equations (DAEs) system for the problem of elastic-plastic wave propagation in an elastic-plastic solid undergoing small deformations. This is achieved by recasting the pointwise complementary trio in the elastic-plastic constitutive equations into an algebraic equation through the Fischer-Burmeister NCP-function. Then, for an isotropicallyhardening/ softening material under prescribed impulse loadings on a thin-walled tube with combined axial-torsional stresses, we can develop a novel algorithm based on the Lie-group differential algebraic equations (LGDAE) method to iteratively solve the resultant DAEs at each time… More >

E. Ferretti¹

CMC-Computers, Materials & Continua, Vol.36, No.1, pp. 49-72, 2013, DOI:10.3970/cmc.2013.036.049

Abstract In several recent papers studying the Cell Method (CM), which is a numerical method based on a truly algebraic formulation, it has been shown that numerical modeling in physics can be achieved even without starting from differential equations, by using a direct algebraic formulation. In the present paper, our focus will be above all on highlighting some of the theoretical features of this algebraic formulation to show that the CM is not simply a new numerical method among many others, but a powerful numerical instrument that can be used to avoid spurious solutions in computational physics. More >

Chih-Wen Chang^1,2, Chein-Shan Liu³

CMC-Computers, Materials & Continua, Vol.34, No.2, pp. 143-175, 2013, DOI:10.3970/cmc.2013.034.143

Abstract The nonlinear Poisson problems in heat diffusion governed by elliptic type partial differential equations are solved by a modified globally optimal iterative algorithm (MGOIA). The MGOIA is a purely iterative method for searching the solution vector x without using the invert of the Jacobian matrix D. Moreover, we reveal the weighting parameter αc in the best descent vector w = αcE + DTE and derive the convergence rate and find a criterion of the parameter γ. When utilizing αc and γ, we can further accelerate the convergence speed several times. Several numerical experiments are carefully discussed and validated the proposed… More >

Hong-Hua Dai^1,2, Jeom Kee Paik³, S. N. Atluri²

CMC-Computers, Materials & Continua, Vol.23, No.2, pp. 155-186, 2011, DOI:10.3970/cmc.2011.023.155

Abstract The application of the Galerkin method, using global trial functions which satisfy the boundary conditions, to nonlinear partial differential equations such as those in the von Karman nonlinear plate theory, is well-known. Such an approach using trial function expansions involving multiple basis functions, leads to a highly coupled system of nonlinear algebraic equations (NAEs). The derivation of such a system of NAEs and their direct solutions have hitherto been considered to be formidable tasks. Thus, research in the last 40 years has been focused mainly on the use of local trial functions and the Galerkin method, applied to the piecewise… More >

Hong-Hua Dai^1,2, Jeom Kee Paik³, Satya N. Atluri²

CMC-Computers, Materials & Continua, Vol.23, No.1, pp. 69-100, 2011, DOI:10.3970/cmc.2011.023.069

Abstract In this paper, the global nonlinear Galerkin method is used to perform an accurate and efficient analysis of the large deflection behavior of a simply-supported rectangular plate under combined loads. Through applying the Galerkin method to the governing nonlinear partial differential equations (PDEs) of the plate, we derive a system of coupled third order nonlinear algebraic equations (NAEs). However, the resultant system of NAEs is thought to be hard to tackle because one has to find the one physical solution from among the possible multiple solutions. Therefore, a suitable initial guess is required to lead to the real solution for… More >

Chia-Ming Fan^1,2, Chein-Shan Liu³, Weichung Yeih¹, Hsin-Fang Chan¹

CMC-Computers, Materials & Continua, Vol.15, No.1, pp. 67-86, 2010, DOI:10.3970/cmc.2010.015.067

Abstract In this study, the nonlinear obstacle problems, which are also known as the nonlinear free boundary problems, are analyzed by the scalar homotopy method (SHM) and the finite difference method. The one- and two-dimensional nonlinear obstacle problems, formulated as the nonlinear complementarity problems (NCPs), are discretized by the finite difference method and form a system of nonlinear algebraic equations (NAEs) with the aid of Fischer-Burmeister NCP-function. Additionally, the system of NAEs is solved by the SHM, which is globally convergent and can get rid of calculating the inverse of Jacobian matrix. In SHM, by introducing a scalar homotopy function and… More >