On the numerical integration of isogeometric interface elements

Zero-thickness interface elements are commonly used in computational mechanics to model material interfaces or to introduce discontinuities. The latter class requires the existence of a non-compliant interface prior to the onset of fracture initiation. This is accomplished by assigning a high dummy stiffness to the interface prior to cracking. This dummy stiffness is known to introduce oscillations in the traction profile when using Gauss quadrature for the interface elements, but these oscillations are removed when resorting to a Newton-Cotes integration scheme 1. The traction oscillations are aggravated for interface elements that use B-splines or non-uniform rational B-splines as basis functions (isogeometric interface elements), and worse, do not disappear when using Newton-Cotes quadrature. An analysis is presented of this phenomenon, including eigenvalue analyses, and it appears that the use of lumped integration (at the control points) is the only way to avoid the oscillations in isogeometric interface elements. New findings have also been obtained for standard interface elements, for example that oscillations occur in the relative displacements at the interface irrespective of the value of the dummy stiffness.