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