The platen velocity dh(t)/dt imposed on the top platen is an
additional scalar unknown that is calculated at each time step with
the platen force balance equation. Indeed, the closing force F
applied on the bottom boundary C is known, but the normal stress
distribution at the boundary is not. It should follow the force balance
equation