Wells stability problems during drilling processes have been a subject of study and interest in the oil and gas industry due to its consequences, such as nonproductive times (NPTs), formation damage, Wells integrity and economic impacts. Isotropic properties in the formation usually had been assumed, however, it is necessary to define more realistic models to represent well stability. In this research, a numerical stress distribution and an analytical analysis have been proposed in order to calculate rock failure around the wellbore and optimize operative mud window, considering anisotropic elastic properties and failure criteria in the intact rock and in the plane of weakness. ABAQUS and MATLAB software were used to represent and solve the numerical-analytical model. The results presented in the assessment proved that the anisotropy consideration (including characteristics of the formation that can induce anisotropy, as plane of weakness) is necessary to be taken in count in this type of investigation because depending of the set of analyzed properties, the range of the operating mud weight window could significantly change.