I was wondering if anyone knows a way to model a support that can only provide a reaction in one direction (i.e., restrained in -z, but free to move in +z)
There may be another way to do it. Just put some dummy element with very high stiffness under the node. The lower node(s) of the dummy element should be restrained in Z direction. Not a perfect solution perhaps. But you will get quite satisfactory results.