Per-topology machinery behind RootRkmkDae: the whole-mesh
Integrable every Gauss point’s plastic state lives on (a
per-Gauss-point List for one block,
a Product of those for Blocks), how to flatten/unflatten between
it and the block’s native per-element State, and how to
evaluate every Gauss point’s plastic rate at a stage’s nodal coordinates.
Composes recursively over nested Blocks and forwards over
ElasticViscoplasticAndElastic, the same way
ElasticViscoplasticElements does.
The RKMK-DAE stage-equilibrium step, mirroring the constitutive-level
RootRkmkDae: F (here, nodal equilibrium) is resolved at every RK stage
abscissa from every Gauss point’s stage-consistent plastic state, rather
than frozen across the load-step window — so the coupling is the tableau’s
own order instead of first order, at the cost of Tab::STAGES equilibrium
solves per window instead of one. One blanket impl over any Model whose
blocks are ElasticViscoplasticDaeElements — a single block, nested
Blocks to any depth, or an ElasticViscoplasticAndElastic pairing.