In this work the compression behavior of masonry was investigated. After a detailed review of code approaches and different research works, a new formula was suggested to describe the compression strength of masonry, based on the mechanical and geometrical properties of its components, when deformation properties of units are larger than the ones of mortar. Later on, a new model, Extended Drucker-Prager Cap Yielding Function, is suggested to describe the three axial compression stress state of mortar in masonry in case deformation properties of mortar are larger than the ones of mortar, and to describe the three axial compression stress state of brick in the other case. This includes defining its parameters based on test diagrams of the mortar material, implementing the model in the numerical software ANSYS, and the numerical results are evaluated for simple cube example. The controlling equations of creep based on the visco-elastic creep theory are presented in the general case of three axial creep under three axial loading conditions. The special case of three axial creep under axial loading is also presented. The “transversal creep” relevant for the compression strength of masonry was discussed and numerical examples have been added to show the effect of changed time-dependent Poisson’s ratio. In another chapter, many examples are presented showing the application of the suggested material models and discontinuous numerical method named eXtended finite element method. Conclusions and recommendations are given in the last chapter.