A method for the effective numerical treatment of the delamination problem in laminated composites is proposed. The interaction between the laminae due to the binding material is described by means of a nonmonotone, possibly multivalued law. Such laws express a variety of limit phenomena related to the discontinuous loss of resistance at the interface. By using the methods of Nonsmooth Mechanics, the problem is written as a hemivariational inequality and is transformed into a substationarity problem for nonconvex nonsmooth energy functions. Using the appropriate finite element discretization, the latter problem is formulated in finite dimensions. A ‘domain decomposition’ type numerical approximation is proposed and a new algorithm is developed for the numerical treatment of the delamination problem. Numerical examples verify the theoretically predicted stability, accuracy and speed of the method and illustrate its properties and applicability.
Read full abstract