Fackler,
Several items to consider:
1) If you don't use sheeting and you just dig into the ground (i.e., well below the water table), you will have a transient condition where the sidewalls of the excavation are under their full hydrostatic condition (i.e., the flow net will not have had time to develop). As such, you will have much larger gradients at the discharge face then you would have if you dug more slowly.
2) If you dig more slowly (like who has the time. . . ) and use a sump, you can dig, dewater and establish the steady-state flow net concurrently.
3) I really don't know whether there is an overburden layer atop the clean sand layer. If so, there could be a phreatic surface in the overburden that could affect the calculation of hydrostatic pressures in the clean sand layer. While the overburden layer may not be yielding free water in the excavation sidewalls, it may be affecting the hydrostatic pressures in the clean sand.
4) After considering the dynamics of the problem, you need to determine the hydraulic gradient at the face (or toe, which may be where the highest gradient is located) of the excavation. For the case of the immediate time after the big cut into the sand layer, this gradient may be HUGE (i.e., well in excess of unity) - no doubt exceeding the "critical gradient" for initiating piping-type failure.
5) If you have any way of controlling the transient gradients there will nonetheless be a steady-state gradient where the phreatic surface is established toward the excavation. Calculate this exit gradient (i.e., deltaH/deltaL (change in head with respect to the change in travel distance). For the case of a "free-surface" flow net (i.e., where the phreatic surface is warped downward to the exit point) and you have a vertical discharge point, you can use the curve equation for a confocal parabola to establish the phreatic surface (you can also just draw it).
6) With the phreatic surface shown, draw you flow net "sqaures" and at the discharge point, determine the change in head for the change in distance over the final square. Compare this gradient to the ratio of GammaB/GammaW (bouyant unit weight divided by unit weight of water).
This is a theoritical outline of the dynamics that may or may not help. There is too much uncertainty in the transient period between the initial cut below the water table and the development of the flow net for the revised boundary conditions. There is no doubt that a quick cut into the sand will show instability. The question remains whether you want to first do some dewawtering to initiate the drawdown of the phreatic surface and then look at the discharge gradients at the subgrade after the excavation is complete.
Warning - I type faster then I do engineering - but this is what comes to mind when I try to assist.
Good luck.
f-d
¡papá gordo ain’t no madre flaca!