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.