e_init intop1(tcd.sbtot)/(2*pi*r_rebar^2) "Initial mechanical strain from average corrosion layer thickness"