(preface : I originally intent to write this as a comment, but the volume limit reached. In that tone, I'll utilize this answer space writing it..)
to simulate those material change behavior in respect to time?
COMSOL, ANSYS or equivalent ring a bell to be.. as long as you have the material properties (young's modulus, density, thermal transport & expansion properties, color?) right.
Also. depending on your actual implementation resolution (nano-scale production.. macro/mini-lego sized assembly), you may want to adjust the mesh size for the finite element method solver.
software for 4D Printing
If your 4th D refers to time.. then the discussion ends here.
If your 4th D is a space (not time) dimension.. Then you need to convert your 4D geometric object into a 3D shadow (just like generating the normal 2D shadow from a 3D geometric object) 1st to proceed. just like an 3D object (eg cube) can have multiple equivalent 2D form (or shadow). The same applies for 4D to 3D. Once you get the 3D coordinates of the chosen 3D form, you may print it as usual.
Software wise.. mathematica (are used in the some of the example I've seen, others unknown/selfcoded) or equivalent will do.. but the key here is not the software.. it's the visualization algorithm. If you can get the "shadowing" algorithm done right, any 3D/math software will do.
p/s : I didn't expect the materials modelling question coming when I read the title (I really thought this question refers to extra space dimension) .. but yeah.. why not. It's a good postgraduate/research topic to take on. ( :