Added a method giveDofPoolSize in EnrichmentItem to allow for a tiny bit more flexible allocation of XFEM dofs.
Modifed TrPlaneStrRot such that the 'whole' B-matrix is computed in one call taking into account the reduced integration.
+ Removes the need for two calls to computeBmatrixAt passing different strain indicies and.
+ The strain can now be computed as usual.
+ Only one integration rule is needed.
Started adding very preliminary structure to support contact. Added node 2 node with penalty and working on enforcment with Lagrange multiplier.
Merge branch 'master' of https://github.com/JimBrouzoulis/OOFEM_Jim