  We constrain the internal dynamics of a stack of 10 clusters from the GCLASS
survey at 0.87<z<1.34. We determine the stack cluster mass profile M(r) using
the MAMPOSSt algorithm of Mamon et al., the velocity anisotropy profile beta(r)
from the inversion of the Jeans equation, and the pseudo-phase-space density
profiles Q(r) and Qr(r), obtained from the ratio between the mass density
profile and the third power of the (total and, respectively, radial) velocity
dispersion profiles of cluster galaxies. Several M(r) models are statistically
acceptable for the stack cluster (Burkert, Einasto, Hernquist, NFW). The total
mass distribution has a concentration c=r200/r-2=4.0-0.6+1.0, in agreement with
theoretical expectations, and is less concentrated than the cluster
stellar-mass distribution. The stack cluster beta(r) is similar for passive and
star-forming galaxies and indicates isotropic galaxy orbits near the cluster
center and increasingly radially elongated with increasing cluster-centric
distance. Q(r) and Qr(r) are almost power-law relations with slopes similar to
those predicted from numerical simulations of dark matter halos. Combined with
results obtained for lower-z clusters we determine the dynamical evolution of
galaxy clusters, and compare it with theoretical predictions. We discuss
possible physical mechanisms responsible for the differential evolution of
total and stellar mass concentrations, and of passive and star-forming galaxy
orbits [abridged].
