The Global Nonlinear Galerkin Method for the Solution of von Karman Nonlinear Plate Equations: An Optimal boolean AND Faster Iterative Method for the Direct Solution of Nonlinear Algebraic Equations F(x)=0, using (x)over dot = lambda[alpha F+(1-alpha)(BF)-F-T]
Hong-Hua Dai,Jeom Kee Paik,S. N. Atluri
2011-01-01
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 linear system of partial differential equations in the updated or total Lagrangean reference frames. This leads to the so-called tangent-stiffness finite element method. The piecewise linear tangent-stiffness finite element equations are usually solved by an iterative Newton-Raphson method, which involves the inversion of the tangent-stiffness matrix during each iteration. However, the advent of symbolic computation has made it now much easier to directly derive the coupled system of NAEs using the global Galerkin method. Also, methods to directly solve the NAEs, without inverting the tangent-stiffness matrix during each iteration, and which are faster and better than the Newton method are slowly emerging. In a previous paper [Dai, Paik and Atluri (2011a)], we have presented an exponentially convergent scalar homotopy algorithm to directly solve a large set of NAEs arising out of the application of the global Galerkin method to von Karman plate equations. While the results were highly encouraging, the computation time increases with the increase in the number of NAEs the number of coupled NAEs solved by Dai, Paik and Atluri (2011a) was of the order of 40. In this paper we present a much improved method of solving a larger system of NAEs, much faster. If F(x) = 0 [F-i(x(j)) = 0] is the system of NAEs governing the modal amplitudes x(j) [j = 1, 2 ... N], for large N, we recast the NAEs into a system of nonlinear ODEs: (x) over dot =lambda[alpha F + (1 - alpha)(BF)-F-T], where A and a are scalars, and B-ij = partial derivative F-i/partial derivative x(j). We derive a purely iterative algorithm from this, with optimum value for A, and a being determined by keeping x on a newly defined invariant manifold [Liu and Atluri (2011b)] Several numerical examples of nonlinear von Karman plates, including the post-buckling behavior of plates with initial imperfections are presented to show that the present algorithms for directly solving the NAEs are several orders of magnitude faster than those in Dai, Paik and Atluri (2011a). This makes the resurgence of simple global Galerkin methods, as alternatives to the finite element method, to directly solve nonlinear structural mechanics problems without piecewise linear formulations, entirely feasible.