I cannot check now what is there, however I think (f) vs. (e) substitution cannot be correct.
The reason is that AISI says "critical stress in top plates occurs between the hole and the free edge of the plate", and AISI presents a specific bending model applied to a "beam" of sectional dimensions as fxc.
That bending stress must be proportional with 1/fc^2 since fc^2 is coming up from section modulus. No "e" appears in this section of calculation.
When calculate bending stress transmitted to shell, the bending stress is proportional with the bending moment and that means is proportional with "P*e". Here "f" cannot be considered.
Best regards.