Analysis of delamination growth with discontinuous solid-like shell elements

Delamination is one of the most important failure mechanisms in laminates. Normally, it is modelled using interface elements. These elements are placed between two layers that are modelled with continuum elements. The interface elements are equipped with a softening or damage model in order to simulate debonding. This method has some drawbacks, both in a numerical and in a mechanical sense. A recent alternative is to simulate the crack by adding a discontinuous displacement mode to the continuum elements according to the partition of unity method. The elements do not contain the discontinuity prior to cracking, but when the ultimate stress in the bulk material is exceeded, delamination is initiated and additional degrees-of-freedom are activated. Beside this, a slightly different implementation is examined also. A discontinuity is predefined and has an initial dummy stiffness. Delamination is initiated when the tractions in the discontinuity exceed a threshold value. The results of both versions of this partition of unity model are compared mutually and with conventional interface elements by means of two examples

Analysis of delamination growth with discontinuous finite elements

In this contribution a new finite element is presented for the simulation of delamination growth in thin layered composite materials. The element is based on the solid-like shell element, a volume element that can be used in very thin applications due to a higher order displacement field in thickness direction. The delamination crack is incorporated in this element as a jump of the displacement field by means of the partition of unity method. The kinematics of the element as well as the finite element formulation are described. The performance of the element is demonstrated by means of two examples

Computational modelling of cracks in viscoplastic media

A newly developed numerical model is used to simulate propagating cracks in a strain softening viscoplastic medium. The model allows the simulation of displacement discontinuities independently of a finite element mesh. This is possible using the partition of unity concept, in which fracture is treated as a coupled problem, with separate variational equations corresponding to the continuous and discontinuous parts of the displacement field. The equations are coupled through the dependence of the stress field on the strain state. Numerical examples show that allowing displacement discontinuities in a viscoplastic Von Mises material can lead to a failure mode that differs from a continuum-only model

Multi-scale and multi-physics modelling for complex materials

No abstract available

Evolving discontinuities and cohesive fracture

Multi-scale methods provide a new paradigm in many branches of sciences, including applied mechanics. However, at lower scales continuum mechanics can become less applicable, and more phenomena enter which involve discon- tinuities. The two main approaches to the modelling of discontinuities are briefly reviewed, followed by an in-depth discussion of cohesive models for fracture. In this discussion emphasis is put on a novel approach to incorporate triaxi- ality into cohesive-zone models, which enables for instance the modelling of crazing in polymers, or of splitting cracks in shear-critical concrete beams. This is followed by a discussion on the representation of cohesive crack models in a continuum format, where phase-field models seem promising

A generalisation of J2-flow theory for polar continua

A pressure-dependent J2-flow theory is proposed for use within the framework of the Cosserat continuum. To this end the definition of the second invariant of the deviatoric stresses is generalised to include couple-stresses, and the strain-hardening hypothesis of plasticity is extended to take account of micro-curvatures. The temporal integration of the resulting set of differential equations is achieved using an implicit Euler backward scheme. This return-mapping algorithm results in an exact satisfaction of the yield condition at the end of the loading step. Moreover, the integration scheme is amenable to exact linearisation, so that a quadratic rate of convergence is obtained when Newton's method is used. An important characteristic of the model is the incorporation of an internal length scale. In finite element simulations of localisation, this property warrants convergence of the load-deflection curve to a physically realistic solution upon mesh refinement and to a finite width of the localisation zone. This is demonstrated for an infinitely long shear layer and for a biaxial specimen composed of a strain-softening Drucker-Prager material