Function is like the rubik's cheese I've published here.
I wanted to make a test of the 2D subsystem of openscad. Result is this cube.

Slicing without XY_size_compensation. I have designed a big gap between the moving parts. (0.3mm to each side is 0.6mm gap in the model.)
Print the number as written in filename. PLA/ABS doesn't matter. Remove all blobs and irregularities carefully. Glue the core parts ( out-in-in-out) together, stabilize with pieces of 3mm ABS filament. Flat sides to outside. Check if work smooth before applying glue. Full assembly with a M3 screw and a selflock nut. Place a washer between the two rotating core elements. I have printed in PLA white. Not good idea, it becomes dirty. Would be better to print in color and in ABS if possible.

scad file: It will generate only full cube, no single parts. I don't recommend a print-in-place because of not-smooth-surfaces inside. Some parts have to be rotated upside-down for print. If you need, import all to your slicer and select split function. (RepetierHost, slic3r, repsnapper provide this).
In scad file, size of cubie, inner diameter, roundness of edges and gap distance can be adjusted freely.

edit 2016-03-01 update. Please use the actualized files.

Twin projects:
https://www.thingiverse.com/thing:3098632
https://www.thingiverse.com/thing:1967151
https://www.thingiverse.com/thing:3072245