I have found the hyperbolic parameters (K and n) according to Duncan and Chang model. However, I have no idea how to fit a hyperbolic model into a failure envelope to get the friction angle.
I believe you need several tests to graph the relationship between (s1+s3)/2 and (s1-s3)/2 at failure (p, q coordinates in American literature [Lambe]).
The failure envelope may be simplified as "qf = c*cos(phi)+p*sin(phi)" or, in cohesionless materials, by a nonlinear function of "phi" versus "s3" like "phi = phi0 + Dphi * log(s3 / patm)".
A single test, duly adjusted to the hyperbolic model for (s1-s3) versus the axial strain, cannot provide the failure envelope. More than three tests in the stress range of your problem may improve the geotechnical analysis with a non-linear failure envelope.