I am trying to write my own code for geometrically nonlinear analysis of mixed structures (cables, trusses, beams...). I would like to have it as general as possible ie. applicable also to 2D and 3D solids. I am trying to follow nonincremental Newton Raphson procedure from Ansys theory manual, but my computation are not converging.. Can I kindly ask you for advice ?
Here is my flow:
=== Linear analysis ===
(1) inv Kt *F = du
resluts are correct comparing to Ansys
=== Nonlinear analysis ===
while (F_ext - F_int) < 1e-3
(2) u = u +du
Recent comments