Method of manufactured solutions code verification of elastostatic solid mechanics problems in a commercial finite element solver

摘要

Much progress has been made in advancing and standardizing verification, validation, and uncertainty quantification practices for computational modeling in recent decades. However, examples of rigorous code verification for solid mechanics problems in literature remain scarce, particularly for commercial software and for the non-trivial large-deformation analyses and nonlinear materials typically needed to simulate medical devices. Here, we apply the method of manufactured solutions (MMS) to verify a commercial finite element code for elastostatic solid mechanics analyses using linear-elastic, hyperelastic (neo-Hookean), and quasi-hyperelastic (Hencky) constitutive models. Analytical source terms are generated using Python/SymPy and are implemented in ABAQUS/Standard without modification to solver source code. Source terms for the three constitutive models are found to vary nearly six orders of magnitude in the number of mathematical operations they contain. Refinement studies reveal second-order displacement and first-order (or higher) stress and strain convergence in response to mesh refinement for all constitutive models and first-order displacement convergence in response to pseudo-time increment refinement for the Hencky-elastic case. We also investigate the sensitivity of convergence order to quantitatively minor changes to the underlying mathematical model using an exploratory case. Code used to generate the MMS source terms and simulation input files are provided as Supplemental Material. Published by Elsevier Ltd.