Crystal plasticity with Jacobian-Free Newton-Krylov

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

    Research output: Contribution to journalArticlepeer-review

    23 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 - May 2013

    All Science Journal Classification (ASJC) codes

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


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

    Cite this