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.

Least-squares finite strain hexahedral element/constitutive coupling based on parametrized configurations and the Löwdin frame

Two novelties are introduced: (i) a finite-strain semi-implicit integration algorithm compatible with current element technologies and (ii) the application to assumed-strain hexahedra. The Löwdin algo- rithm is adopted to obtain evolving frames applicable to finite strain anisotropy and a weighted least- squares algorithm is used to determine the mixed strain. Löwdin frames are very convenient to model anisotropic materials. Weighted least-squares circumvent the use of internal degrees-of-freedom. Het- erogeneity of element technologies introduce apparently incompatible constitutive requirements. Assumed-strain and enhanced strain elements can be either formulated in terms of the deformation gradient or the Green–Lagrange strain, many of the high-performance shell formulations are corotational and constitutive constraints (such as incompressibility, plane stress and zero normal stress in shells) also depend on specific element formulations. We propose a unified integration algorithm compatible with possibly all element technologies. To assess its validity, a least-squares based hexahedral element is implemented and tested in depth. Basic linear problems as well as 5 finite-strain examples are inspected for correctness and competitive accuracy.