Crystal plasticity with Jacobian-Free Newton-Krylov

K. Chockalingam, Michael Tonks, J. D. Hales, D. R. Gaston, P. C. Millett, Liangzhe Zhang

Research output: Contribution to journalArticlepeer-review

19 Scopus citations


The objective of this work is to study potential benefits of solving crystal plasticity finite element method (CPFEM) implicit simulations using the Jacobian-Free Newton-Krylov (JFNK) technique. Implicit implementations of CPFEM are usually solved using Newton's method. However, the inherent non-linearity in the flow rule model that characterizes the crystal slip system deformation on occasions would require considerable effort to form the exact analytical Jacobian needed by Newton's method. In this paper we present an alternative using JFNK. As it does not require an exact Jacobian, JFNK can potentially decrease development time. JFNK approximates the effect of the Jacobian through finite differences of the residual vector, allowing modified formulations to be studied with relative ease. We show that the JFNK solution is identical to that obtained using Newton's method and produces quadratic convergence. We also find that preconditioning the JFNK solution with the elastic tensor provides the best computational efficiency.

Original languageEnglish (US)
Pages (from-to)617-627
Number of pages11
JournalComputational Mechanics
Issue number5
StatePublished - Jan 1 2013

All Science Journal Classification (ASJC) codes

  • Ocean Engineering
  • Mechanical Engineering
  • Computational Theory and Mathematics
  • Computational Mathematics
  • Applied Mathematics


Dive into the research topics of 'Crystal plasticity with Jacobian-Free Newton-Krylov'. Together they form a unique fingerprint.

Cite this