Phase-field analysis of finite-strain plates and shells including element subdivision

With the theme of fracture of finite-strain plates and shells based on a phase-field model of crack regularization, we introduce a new staggered algorithm for elastic and elasto-plastic materials. To account for correct fracture behavior in bending, two independent phase-fields are used, corresponding to the lower and upper faces of the shell. This is shown to provide a realistic behavior in bending-dominated problems, here illustrated in classical beam and plate problems. Finite strain behavior for both elastic and elasto-plastic constitutive laws is made compatible with the phase-field model by use of a consistent updated-Lagrangian algorithm. To guarantee sufficient resolution in the definition of the crack paths, a local remeshing algorithm based on the phase- field values at the lower and upper shell faces is introduced. In this local remeshing algorithm, two stages are used: edge-based element subdivision and node repositioning. Five representative numerical examples are shown, consisting of a bi-clamped beam, two versions of a square plate, the Keesecker pressurized cylinder problem, the Hexcan problem and the Muscat-Fenech and Atkins plate. All problems were successfully solved and the proposed solution was found to be robust and efficient.

Semi-implicit finite strain constitutive integration and mixed strain/stress control based on intermediate configurations

A new semi-implicit stress integration algorithm for finite strain plasticity (compatible with hyperelas- ticity) is introduced. Its most distinctive feature is the use of different parameterizations of equilibrium and reference configurations. Rotation terms (nonlinear trigonometric functions) are integrated explicitly and correspond to a change in the reference configuration. In contrast, relative Green–Lagrange strains (which are quadratic in terms of displacements) represent the equilibrium configuration implicitly. In addition, the adequacy of several objective stress rates in the semi-implicit context is studied. We para- metrize both reference and equilibrium configurations, in contrast with the so-called objective stress integration algorithms which use coinciding configurations. A single constitutive framework provides quantities needed by common discretization schemes. This is computationally convenient and robust, as all elements only need to provide pre-established quantities irrespectively of the constitutive model. In this work, mixed strain/stress control is used, as well as our smoothing algorithm for the complemen- tarity condition. Exceptional time-step robustness is achieved in elasto-plastic problems: often fewer than one-tenth of the typical number of time increments can be used with a quantifiable effect in accuracy. The proposed algorithm is general: all hyperelastic models and all classical elasto-plastic models can be employed. Plane-stress, Shell and 3D examples are used to illustrate the new algorithm. Both isotropic and anisotropic behavior is presented in elasto-plastic and hyperelastic examples.