It got lost in this overly long thread, but there’s an official example from the CalculiX test suite here: Modeling bolt connection with preloaded beam element - #5 by FEAnalyst
Try adding equation constraints. The approach with distributing couplings can be problematic, though.