  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 the 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 either Python/SymPy or
Mathematica and are implemented in ABAQUS/Standard without modification to the
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 convergence
in response to mesh refinement for all constitutive models and first-order
displacement convergence in response to increment refinement for the
finite-strain problems. We also investigate the sensitivity of MMS convergence
order to minor coding errors using an exploratory case. Code used to generate
the MMS source terms and the input files for the simulations are provided as
supplemental material.
