I don't know if the algorithm is ok, for instance taking a very low Δd it seems that ϕVc.min ≈ Vu.max. For that maybe you need to increase MaxIters.

The trick abut f'c works just here, it could be not very general.

Adventages of using for

First, you can arrange the algorithm code in such way that not repeat calculations, in the sense that you don't need to evaluate ϕVc.min or max before the loop.

Second, you avoid infinity loops, limiting the number of avaible iterations in an easy way.

Also you can know how many iterations the code take

and check if it ends because the maximun iterations limit has been rached or because the break

Notice that one means "true"

Notas acerca de las undiades (en español)

Lo que tu defines como Kg se llama kilopondio, y aún lo siguen usando los ingenieros civiles, aunque sigue en descenso. Puedes verlo en https://es.wikipedia.org/wiki/Kilopondio

Sin embargo ahí no dan una abreviatura. Las abreviaturas para kilopondio (kp) y tonelada fuerza (Mp = Mega pondio) las puedes encontrar en la página en inglés: https://en.wikipedia.org/wiki/Kilogram-force

Finalmente, para la unidad kgf/cm^2 (vamos, los 'kilos' de presión de toda la vida) también se llama atmósfera técnica, pues es muy parecido en valor a una atm de presión, o un barn. Puedes leer más de eso (y del kilopondio) aquí: https://es.wikipedia.org/wiki/Sistema_T%C3%A9cnico_de_Unidades

Y si alguien te critica por andar citando la wikipedia, puedes responder que es interesante pues, podrá ser correcto o no lo que ahí esté, pero es lo que mucha gente piensa sobre el tema, y eso lo hace relevante.