Finite-strain laminates: Bending-enhanced hexahedron and delamination

With a new finite strain anisotropic framework, we introduce a unified approach for constitutive model- ing and delamination of composites. We describe a finite-strain semi-implicit integration algorithm and the application to assumed-strain hexahedra. In a laminate composite, the laminae are modeled by an anisotropic Kirchhoff/Saint-Venant material and the interfaces are modeled by the exponential cohesive law with intrinsic characteristic length and the criterion by Benzeggagh and Kenane for the equivalent fracture toughness. For the element formulation, a weighted least-squares algorithm is used to calculate the mixed strain. Löwdin frames are used to model orthotropic materials without the added task of per- forming a polar decomposition or empirical frames. To assess the validity of our proposals and inspect step and mesh size dependence, a least-squares based hexahedral element is implemented and tested in depth in both deformation and delamination examples.

A finite element framework for the interplay between delamination and buckling of rubber-like bi-material systems and stretchable electronics

In this study, a finite element (FE) framework for the analysis of the interplay between buckling and delamination of thin layers bonded to soft substrates is proposed. The current framework incorporates the following modeling features: (i) geometrically nonlinear solid shell elements, (ii) geometrically nonlinear cohesive interface elements, and (iii) hyperelastic material constitutive response for the bodies that compose the system. A fully implicit Newton–Raphson solution strategy is adopted to deal with the complex simultaneous presence of geometrical and material nonlinearities through the derivation of the consistent FE formulation. Applications to a rubber-like bi-material system under finite bending and to patterned stiff islands resting on soft substrate for stretchable solar cells subjected to tensile loading are proposed. The results obtained are in good agreement with benchmark results available in the literature, confirming the accuracy and the capabilities of the proposed numerical method for the analysis of complex three-dimensional fracture mechanics problems under finite deformations.