gamma_E = np.sqrt(E / f_y)* np.pi
gamma_E = np.sqrt(E / f_y)* np.pi