Negative jacobian with tie constraints means that you should disable Adjust for them. Convergence is much more tricky and can have many causes (although it usually originates from contact). To find them, it’s best to simplify the model (reduce its size e.g. to segment with one bolt and gradually disable some advanced features such as the 3 nonlinearities).
At the first glance, it seems that you could utilize symmetry in one plane. Also, you should improve the mesh - refine it for shells and use hexa or at least second-order tetra elements for the bolts. You will find many examples of similar models on the forum if you search for solid bolts. Such as this thread: COMPRESSION ONLY for simulating bolts and nut faces Or this one: Error during analysis of baseplate fixing
Displacement-based preloading is much better for convergence.
In the second step, where you apply the actual loading, you should freeze the preload (Displacement → Fixed).