The stable finite element method for minimization problems

B. Liu*, Y. Huang

*Corresponding author for this work

Research output: Contribution to journalArticle

8 Scopus citations

Abstract

The conventional finite element method is difficult to converge for a non-positive definite stiffness matrix, which usually occurs when the material displays softening behavior or when the system is near the state of bifurcation. We have developed two stable algorithms for a non-positive definite stiffness matrix, one for the direct linear equation solver and the other for the iterative solver in the finite element method for minimization problems. For a direct solver with non-positive definite stiffness matrix, energy minimization of a system with multiple degrees of freedom (DOF) is decomposed to the minimization of many 1-DOF systems, and for the latter an efficient and robust minimization method is developed to ensure that the system energy decreases in every incremental step, regardless of the positive definiteness of the stiffness matrix. For an iterative solver, the stiffness matrix is modified to ensure the convergence, and the modified stiffness matrix indeed leads to the correct solution. An example of a single wall carbon nanotube under compression is studied via the proposed algorithms.

Original languageEnglish (US)
Pages (from-to)1251-1254
Number of pages4
JournalJournal of Computational and Theoretical Nanoscience
Volume5
Issue number7
StatePublished - Jul 2008

    Fingerprint

Keywords

  • Atomic-scale finite element method
  • Direct solver
  • Iterative solver
  • Minimization

ASJC Scopus subject areas

  • Chemistry(all)
  • Materials Science(all)
  • Condensed Matter Physics
  • Computational Mathematics
  • Electrical and Electronic Engineering

Cite this