The stable finite element method for minimization problems

B. Liu*, Yonggang Huang

*Corresponding author for this work

Research output: Contribution to journalArticle

8 Citations (Scopus)

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
DOIs
StatePublished - Jul 1 2008

Fingerprint

stiffness matrix
Stiffness matrix
Stiffness Matrix
Minimization Problem
finite element method
Finite Element Method
Finite element method
optimization
Iterative Solver
degrees of freedom
Degree of freedom
Positive Definiteness
Carbon Nanotubes
Bifurcation (mathematics)
Energy Minimization
Softening
linear equations
Linear equations
Nanotubes
softening

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

@article{446dd4b63914414bb2eab5fbbb724df1,
title = "The stable finite element method for minimization problems",
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.",
keywords = "Atomic-scale finite element method, Direct solver, Iterative solver, Minimization",
author = "B. Liu and Yonggang Huang",
year = "2008",
month = "7",
day = "1",
doi = "10.1166/jctn.2008.007",
language = "English (US)",
volume = "5",
pages = "1251--1254",
journal = "Journal of Computational and Theoretical Nanoscience",
issn = "1546-1955",
publisher = "American Scientific Publishers",
number = "7",

}

The stable finite element method for minimization problems. / Liu, B.; Huang, Yonggang.

In: Journal of Computational and Theoretical Nanoscience, Vol. 5, No. 7, 01.07.2008, p. 1251-1254.

Research output: Contribution to journalArticle

TY - JOUR

T1 - The stable finite element method for minimization problems

AU - Liu, B.

AU - Huang, Yonggang

PY - 2008/7/1

Y1 - 2008/7/1

N2 - 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.

AB - 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.

KW - Atomic-scale finite element method

KW - Direct solver

KW - Iterative solver

KW - Minimization

UR - http://www.scopus.com/inward/record.url?scp=57749090495&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=57749090495&partnerID=8YFLogxK

U2 - 10.1166/jctn.2008.007

DO - 10.1166/jctn.2008.007

M3 - Article

VL - 5

SP - 1251

EP - 1254

JO - Journal of Computational and Theoretical Nanoscience

JF - Journal of Computational and Theoretical Nanoscience

SN - 1546-1955

IS - 7

ER -