Computational procedure within the framework of return mapping technique has been presented to integrate the constitutive behavior of a concrete model. Developed by Ohtani and Chen, this concrete model is based on multiple hardening concept, and is rate-independent and associative. Consistent tangent operator suitable for finite element analysis is derived to preserve the rate of convergence. Accuracy of the integration technique is verified and compared with available experimental data. Computational efficiency is demonstrated by comparing with results based on elasto-plastic tangent.