Overcoming Convergence Difficulties in ANSYSWorkbench Mechanical, Part I: Using Newton-Raphson Residual Information
padtinc.com
/blog/the-focus/overcoming-convergence-difficulties-in-ansys-workbench-mechanical-part-i-using-newton-raphson-residual-informationTed HarrisUnable to converge. Convergence Failure. Failure to Converge. Never nice words to see when youare trying to get your simulation done. If you’ve encountered convergence failures while running nonlinear structural analyses in ANSYSWorkbench Mechanical, this two part series is for you. What is a convergence failure? In a nutshell itmeans that there is too much imbalance in the system. The calculated reaction forces do not matchthe applied loads and even though the program tries hard to make changes to overcome theimbalances, it hasn’t been able to do so and stops. If we look at the Force residuals under SolutionInformation, we will see that the solver has been unable to get the force convergence residual, or imbalance force, to drop below the current criterion
Test model example: Newton Raphson Convergence Failure; Solution Stops
We won’t spend a lot of time here explaining the Newton-Raphson method, convergence, and residualplots here, since we wrote a Focus article back in 2002 which discusses them in more detail. Thearticle begins on p. 7 at this link:/blog/wp-content/uploads/oldblog/archive/PADT_TheFocus_08.pdf The context of that article was Mechanical APDL, but the article is directly relevant since solving in
Workbench Mechanical is done in Mechanical APDL in batch mode. In crayon terms, we want the purple line to drop below the blue line. When it doesn’t and the solver isout of options to keep trying, the solution stops and we get an error message. Now what? The traditional knobs to turn are to increase the number of substeps, decrease contactstiffness if contact is involved, perhaps add more points to the plasticity curve, etc. But what if something else is the problem? How can we identify where the problem is?In this part I article we will discuss how to plot the Newton-Raphson residuals as contour plots to seewhere in the model the highest force imbalances are located. Often this is useful information to helpus figure out what is going on so we can take corrective action. First, be aware that we must turn onthe Newton-Raphson residual plots prior to solving. That means you either have to turn them on andre-solve after a convergence failure, knowing that you’ll get the same failure again, or you need toclairvoyantly (or perhaps just prudently) turn on the residuals prior to attempting the initial solve. Whyaren’t they on all the time, you ask? Most likely because they slow things down just a bit and alsorequire a bit more disk space than otherwise, although if the solution runs to completion no Newton-Raphson residual plots are saved.Here is how we turn them on. In the Details view for the Solution Information branch, change theNewton-Raphson Residuals setting from the default of zero to a nonzero number such as 3 or 4. Thatwill continuously save the last 3 or 4 Newton-Raphson residual plots for viewing as contour plots after the solution has stopped due to a convergence failure. After the solution has stopped, the Newton-Raphson residual plots will be available under the Solution Information branch.The quantity plotted is actually the square rootof the sum of the squares of the residuals inthe global X, Y, and Z directions. So, the plotsdon’t show us direction information, but they doshow where the residuals and hence the forceimbalances are the largest. Below is anexample. The region in red shows where theresiduals are the highest. Since this is a modelinvolving contact between two bodies,apparently the contact regions and specificallycontact at the corners of the part on the left is the source of our convergence difficulties.
Newton-Raphson Residual Force Plot for the last attempted equilibrium iteration.
So, how do we use this information? In this case we now suspect that the contact regions, especiallyat the corners of the smaller part, are the problematic areas. Using this information we made twochanges to the model. First, we changed the Detection Method for the contact elements from Program Controlled (at theelement Gauss points) to Nodal-Normal to Target. Many times when contact problems involve touchingat corners, the robustness of the contact interface can be improved by changing the detection methodfrom Gauss points to nodes.Second, we reduced the contact stiffness by changing Normal Stiffness from Program Controlled(factor of 1.0) to a Manual setting of 0.2. Reducing the contact stiffness can help with contactconvergence for a lot of problems. Too low of a stiffness value can cause problems too, but in thiscase the resulting penetration is still small so a value of 0.2 seems reasonable. When in doubt, asensitivity study can be performed whereby you make changes to the contact stiffness value whiletracking your results quantities of interest. As with most inputs you can vary, your results of interestshould not be sensitive to contact stiffness.These two changes allowed our test model to nicely converge for the full amount of load.
Other considerations:
The Newton-Raphson Residual plots are always displayed on the original geometry, not the deflectedgeometry at version 14.0 of ANSYS Mechanical. If the deflections are large this can make it harder toascertain what is causing the high residual values. In those cases, it can be helpful to compare thetotal deformation and stress plots for the unconverged solution, along with those plots for the lastconverged solution, with the 1.0 true scale on the deformation active. This will show the parts in their deflected state, and that can help in determining why the residuals are high at certain locations.We recommend creating at least 3 residual plots (set in the details of Solution Information as describedabove). Sometimes the location of the imbalance can bounce around a bit from equilibrium iteration toequilibrium iteration, so having more than one or two plots to look at can be beneficial in determiningproblem locations.
Conclusion
Summing it up, the Newton-Raphson residual plots are one piece of information we can use todetermine why we are having convergence difficulties. They can give us an indication of where theconvergence difficulties are occurring in the model, and many times we can use that information to helpus know what settings should be modified or what other changes should be made to the model toimprove the convergence behavior.In part II of this article, we’ll look at how to quickly use ANSYS Mechanical APDL to view the elementsthat have undergone too much deformation.