

June 28, 2021 at 1:28 pmEmperorSubscriber
Hello to all,
I would like to ask for your expertise on finite elements in order to know how to make the right choice on the different methods of solving the finite element calculation.
In Ansys, the resolution of the equation K*U=F is done by the NewtonRaphson method. Several resolution techniques exist, among them: the full, asymmetric and modified NewtonRaphson method (respectively NROPT,full NROPT,unsym NROPT,modi).
The difference in its resolutions is that for the first two methods, the stiffness matrix is updated at each iteration contrary to the last one which updates the stiffness matrix only on the first iterations (so there is less inversion of the matrix and less formulation). The final choice of these two methods could be done by what means? (we can always proceed by scanning the cases and see the difference of the results, which can be long I think)
Are there some structures that require an update of the stiffness matrix at each iteration (I guess yes as for hyperelastic materials) and others not? If the stiffness matrix is not updated at each iteration the accuracy is not the same I think, so the last method of solving should not exist?
thank,

June 29, 2021 at 1:06 pm1shanAnsys EmployeeCheck out the 2 images below. 
The first one is a regular newton raphson and the second one is modified newton raphson. In the regular NR method you can see that at each incremental displacement the tangent slope is calculated(slope is decreasing), while in the modified NR method the slope is calculated just once, at start and the same slope is used (all lines are parallel)to progress ahead till convergence. For regular NR since you calculate tangent at each point you know that you are moving in the "right direction" and hence only a few iterations are needed. But the "right direction" comes at a cost  you need to evaluate stiffness matrix at each point which is a computationally expensive. On the other hand for modified NR since the stiffness is not updated at each point, you don't know if you are moving in the "right direction" or not and hence more iterations are usually needed than regular NR to achieve convergence. But since you are computing stiffness matrix only once, each iteration is faster. In short, Regular NR  less iterations but each iteration is time consuming, Modified NR more iterations but each iteration is faster.
Regards Ishan.

June 29, 2021 at 1:39 pmEmperorSubscriberThank you for the answer because it is much clearer.
In a simulation how to choose the most efficient one? or is it arbitrary?
For a behavior with more curvature it is clear that the update of the stiffness matrix is absolute.

 You must be logged in to reply to this topic.
Simulation World 2022
Check out more than 70 different sessions now available on demand. Get inspired as you hear from visionary companies, leading researchers and educators from around the globe on a variety of topics from lifesaving improvements in healthcare, to bold new realities of space travel. Take a leap of certainty and check out a session today here.
Earth Rescue – An Ansys Online Series
The climate crisis is here. But so is the human ingenuity to fight it. Earth Rescue reveals what visionary companies are doing today to engineer radical new ideas in the fight against climate change. Click here to watch the first episode.
Ansys Blog
Subscribe to the Ansys Blog to get great new content about the power of simulation delivered right to your email on a weekly basis. With content from Ansys experts, partners and customers you will learn about product development advances, thought leadership and trends and tips to better use Ansys tools. Sign up here.
 How to calculate the residual stress on a coating by Vickers indentation?
 An Unknown error occurred during solution. Check the Solver Output…..
 Saving & sharing of Working project files in .wbpz format
 Solver Pivot Warning in Beam Element Model
 Understanding Force Convergence Solution Output
 Colors and Mesh Display
 whether have the difference between using contact and target bodies
 The solver engine was unable to converge on a solution for the nonlinear problem as constrained.
 What is the difference between bonded contact region and fixed joint
 Massive amount of memory (RAM) required for solve

2024

1730

955

736

413
© 2022 Copyright ANSYS, Inc. All rights reserved.