Mixed variational continuum mechanics modeling: a multidomain subpotential approach