Discussion Forum

Hello,

I am looking for a way to get integration points at my element nodes. This is needed to get a more stable third medium contact model. You can find the explanation for it in this paper.

When not using COMSOL, I usually resort to Newton–Cotes quadrature rules or Gauss-Lobatto quadrature, which are both very standard methods. However, I couldn't find these anywhere in the documentation; it seems that only the standard Gaussian quadrature is available. Simply increasing the number of standard Gauss integration points, in order to get some of the points closer to the element nodes, is computationally very inefficient and does not provide the stability that I am looking for. Getting Newton–Cotes or Gauss-Lobatto to work in COMSOL seems like something that, in theory, should be very easy to implement; it is just a matter of defining points and assigning them weights.

Any approach to achieve this when using Q4, Q8, Q9, or LST elements in 2D (or their equivalent in 3D) is highly appreciated.

/Andreas