Abstract
In this paper, a practical and robust Riks-like path following method is developed to automatically trace the nonlinear equilibrium path for initial-finite-strain-induced snap-through phenomena in soft solids. The deformation gradient is decomposed into an elastic part and an initial-strain part, where the initial strain is characterized through a prescribed pattern tensor and an unknown strain factor describing the magnitude. In the non-linear finite element (FE) solution scheme, an equivalent body force proportional to the increment of the strain factor is derived, which varies in each iteration step due to its dependence on deformed configuration and stress. Both the displacement and the strain factor increments are introduced as unknowns in the linearized equation, which permits to predict the critical equilibrium states before and after the snap-through. Several examples with various configurations and strain histories are constructed, including multilayer beams and film substrate systems, which demonstrate the effectiveness and robustness of the algorithm to capture the strain-induced snap-through phenomena. (C) 2017 Elsevier B.V. All rights reserved.