Hi Mickey and everyone,
- The number of iterations is set from the Analysis Parameters ("Setup > Model Setup > Parameters > Contact" tab).
- There is no control for the convergence criteria, such as "within X% of the solution". Either the contact elements that are in contact remain the same for two iterations, or the solution needs to continue. (Other analysis types such as thermal can get within X% of the converged temperature, which the user can adjust to make the analysis faster at the expense of accuracy.)
I have had success changing the contact iteration method from “Mixed” to “Multiple”. This control is also on the "Setup > Model Setup > Parameters > Contact" tab. My notes also indicates that contact models will solve with the iterative solver, but it appears that the sparse solver is always chosen by default. (The sparse solver is slower than the iterative solver on large models, and 1.3E6 equations in your analysis qualifies as a large model.) My guess is that this was done because the sparse solver can handle under-constrained iterations better than the iterative solver, so make sure the model is well stabilized before changing to the iterative solver. (somewhere under the Analysis Parameters)
The page "Perform Analyses with Gap Elements" in the documentation describes the iterations that you see on the screen (and in the log file when you run the analysis through the software). Come to think of it, I do not know what "close" and "open" mean when using tension Gap elements. My guess is that "closed" tension element is able to transmit tension, and an open tension element does not transmit tension or compression. (Naturally, a closed gap element transmits compression, and an opened gap element does not transmit compression or tension.)
I think the best thing to help with the contact iterative solution is to make a statically stable model without relying on the contact elements. That is, if you were to define the contact as "Free" instead of "Surface" (or suppress the Gap elements for a hand-built contact model), the solver should be able to find a static solution. This is normally done by applying "weak" 3D springs to each subassembly in order to create the needed stability. (I'm calling a group of bonded parts a "subassembly".) The page referenced above gives some guidelines on setting the stiffness of the springs in order to make them "weak". (I checked an old version of the documentation, so I do not know if the new documentation refers to "3D springs", or if it still refers to the outdated terminology "rigid boundary".)