SUMMARYIn this paper, the earlier formulation of the SHB8PS finite element is revised in order to eliminate some persistent membrane and shear locking phenomena. This new formulation consists of a solid-shell element based on a purely three-dimensional approach. More specifically, the element has eight nodes, with displacements as the only degrees of freedom, as well as an arbitrary number of integration points, with a minimum number of two, distributed along the 'thickness' direction. The resulting derivation, which is computationally efficient, can then be used for the modeling of thin structures, while providing an accurate description of the various through-thickness phenomena. A reduced integration scheme is used to prevent some locking phenomena and to achieve an attractive, low-cost formulation. The spurious zero-energy modes due to this in-plane one-point quadrature are efficiently controlled using a physical stabilization procedure, whereas the strain components corresponding to locking modes are eliminated with a projection technique following the assumed strain method. In addition to the extended and detailed formulation presented in this paper, particular attention has been focused on providing full justification regarding the identification of hourglass modes in relation to rank deficiencies. Moreover, an attempt has been made to provide a sound foundation to the derivation of the co-rotational coordinate frame, on which the calculations of the stabilization stiffness matrix and internal load vector are based. Finally to assess the effectiveness and performance of this new formulation, a set of popular benchmark problems is investigated, involving geometric non-linear analyses as well as elastic-plastic stability issues.

No abstract

Sheet metal forming processes generally involve large deformations together with complex loading sequences. In order to improve numerical simulation predictions of sheet parts forming, physically-based constitutive models are often required. The main objective of this paper is to analyze the strain localization phenomenon during the plastic deformation of sheet metals in the context of such advanced constitutive models. Most often, an accurate prediction of localization requires damage to be considered in the finite element simulation. For this purpose, an advanced, anisotropic elastic-plastic model, formulated within the large strain framework and taking strain-path changes into account, has been coupled with an isotropic damage model. This coupling is carried out within the framework of continuum damage mechanics. In order to detect the strain localization during sheet metal forming, Rice's localization criterion has been considered, thus predicting the limit strains at the occurrence of shear bands as well as their orientation. The coupled elastic-plastic-damage model has been implemented in Abaqus/Implicit. The application of the model to the prediction of Forming Limit Diagrams (FLDs) provided results that are consistent with the literature and emphasized the impact of the hardening model on the strain-path dependency of the FLD. The fully threedimensional formulation adopted in the numerical development allowed for some new results -e.g. the out-of-plane orientation of the normal to the localization band, as well as more realistic values for its in-plane orientation.