You can also try to reduce the vellum cloth restlength constraint during the sim as the surface becomes compressed to reduce wrinkles.
In a wrangle you could do something like
f@restlength = min(f@restlength, primitiveintrinsic(0, 'measuredperimeter', @primnum);
which will keep the lower value between the existing constraint and the current edge length of the constraint.