The nonlinear finite element method with eight noded isoparametric quadrilateral element for concrete and two noded element for reinforcement is used for the prediction of the behavior of reinforcement concrete structures. The disturbed state concept (DSC) including the hierarchical single surface (HISS) plasticity model with associated flow rule with modifications is used to characterize the constitutive behavior of concrete both in compression and in tension which is named DSC/HISS-CT. The HISS model is applied to shows the plastic behavior of concrete, and DSC for microcracking, fracture and softening simulations of concrete. It should be noted that the DSC expresses the behavior of a material element as a mixture of two interacting components and can include both softening and stiffening, while the classical damage approach assumes that cracks (damage) induced in a material treated acts as a void, with no strength. The DSC/HISS-CT is a unified model with different mechanism, which expresses the observed behavior in terms of interacting behavior of components; thus the mechanism in the DSC is much different than that of the damage model, which is based on physical cracks which has no strength and interaction with the undamaged part. This is the first time the DSC/HISS-CT model, with the capacity to account for both compression and tension yields, is applied for concrete materials. The DSC model allows also for the characterization of non-associative behavior through the use of disturbance. Elastic perfectly plastic behavior is assumed for modeling of steel reinforcement. The DSC model is validated at two levels: (1) specimen and (2) practical boundary value problem. For the specimen level, the predictions are obtained by the integration of the incremental constitutive relations. The FE procedure with DSC/HISS-CT model is used to obtain predictions for practical boundary value problems. Based on the comparisons between DSC/HISS-CT predictions, test data and ANSYS software predictions, it is found that the model provides highly satisfactory predictions. The model allows computation of microcracking during deformation leading to the fracture and failure; in the model, the critical disturbance, Dc, identifies fracture and failure.