A probabilistic approach to the Φ-variation of classical fractal functions with critical roughness