jax_fem.solver.dynamic_relax_solve

jax_fem.solver.dynamic_relax_solve(problem, tol=1e-06, nKMat=50, nPrint=500, info=True, info_force=True, initial_guess=None)

Implementation of

Luet, David Joseph. Bounding volume hierarchy and non-uniform rational B-splines for contact enforcement in large deformation finite element analysis of sheet metal forming. Diss. Princeton University, 2016. Chapter 4.3 Nonlinear System Solution

Particularly good for handling buckling behavior. There is a FEniCS version of this dynamic relaxation algorithm. The code below is a direct translation from the FEniCS version.

TODO: Does not support periodic B.C., need some work here.