,but how can Newton-Rapshon and direct solver go hand in hand? Since non-linearity is involved, there needs to be iterations made using the Newton-Raphson (where the displacements are guessed and then iterations proceed). But the direct solver uses a LU decomposition within the stiffness matrix to solve for the displacements without guessing of the displacements. This thing kind of boggles me.