equilibrium

Solvers

cra_solve	CRA solver using Pyomo + IPOPT.
cra_penalty_solve	CRA solver with penalty formulation using Pyomo + IPOPT.
rbe_solve	RBE solver with penalty formulation using Pyomo + IPOPT.

---

The following helper functions can be useful if you’re developing your own formulation.

Equilibrium Helper Functions

equilibrium_setup	Set up equilibrium matrix.
friction_setup	Set up friction matrix.
external_force_setup	Set up external force vector.
density_setup	Set up material density.
make_aeq	Create equilibrium matrix Aeq or penalty formulation matrix Aeq@B.
make_afr	Create friction matrix Afr and Afr@B.
unit_basis	Create interface reference system as unit basis.
num_vertices	Total number of vertices.
num_free	Return number of free blocks.
free_nodes	Return free and fixed node list.

Pyomo Helper Functions

initialisations	Variable initialisations for pyomo.
bounds	Variable bounds for pyomo.
objectives	Objective functions for pyomo.
constraints	Constraint functions for pyomo.
static_equilibrium_constraints	Create equilibrium and friction constraints.
pyomo_result_check	Check if pyomo optimisation result, raise error if the problem is infeasible.
pyomo_result_assembly	Save pyomo optimisation results to assembly.