I'm working with the Boundary Element Method and want to integrate an expression over a triangular region. I would like to use Gauss Integration to do this, but I'm having trouble since the triangular element is zero order(1 degree lower than a constant stress triangle). Unlike the higher order elements, the zero order triangle has only one node at its centroid. Normally, Gauss Integration is carried out over a triangular element by defining a natural coordinate system and corresponding shape functions/Jacobin. In the case of the zero order element, however, their is one shape function = 1. Therefore the shape functions can not be expressed in terms of the natural coordinates, and by extension, in terms of the Gauss Points.
Does anyone know how to integrate an expression over a zero order triangular element using Gauss Integration? If not does anyone know a different numerical integration scheme to integrate over a triangular region by just plugging in points and weights. I don't want to write out both integrals and then apply a numerical method to each one. I would also like to avoid defining the limits of integration.