Dear tuderic,
unfortunately no. I’ve implemented a version, which takes care of the problem that dokken pointed out (imposing contraints on DOFs twice), though. I basically assign each edge to only one face (but I think I did not take care of the corner points of the cube). It still doesn’t really work.
After that I didn’t have much time to work on this problem, but it’s still on my agenda. If you and others are also interested, I could create a github repo and push what I’ve coded till now and we could try to figure it out or so 
I think that would be better than posting code snippets here.
Best,
welahi