Abstarct: Abstract Soil behavior is influenced by several factors such as confined pressure, density ratio, particle gradation, grain shape, and fabric. The shear modulus is an important characteristic for sands under seismic loads. A piezoelectric ceramic, i.e., bender element, can be used to measure the velocity of shear wave and, consequently, the shear modulus of the soil. The bender element is also a potent tool to investigate the effect of anisotropy on consolidation. Previous studies indicate that the initial shear modulus of the soil during consolidation and prior to any drastic fabric evolution depends on the fabric and stress conditions. The main innovation of this research is the simultaneous evaluation of the effect of loading type, relative density of samples and the degree of anisotropy during the consolidation and shear stages through laboratory tests and empirical models. Using triaxial and bender element tests, the effect of anisotropy on the general behavior and shear modulus changes of the samples has been evaluated through isotropic and anisotropic consolidations of Firoozkouh 161 sand under various compressive and tensile stress ratios, considering different density ratios in the range of 24% to 43%. The results demonstrated that samples consolidated with stress ratios of -0.6 and +0.6 through axial compression tests experienced the highest and lowest deviatoric stress changes at the end of consolidation, respectively. Under axial extension tests, for the same stress ratios of -0.6 and +0.6, the samples had the lowest and highest variations of deviatoric stress values, respectively. The results of the bender element tests show that for samples with density ratios of 27% and 41%, the minimum values of the shear modulus occur in the strain range of -7.5% to -11.5%. For the applied frequencies, the maximum phase changes of the waves occurred in the range of 23 kHz to 31 kHz. Therefore, it can be concluded that the critical frequency for the onset of severe evolution in the fabric of the loose Firoozkooh samples is, on average, about 28 kHz. This value was obtained for samples with an average relative density of Firoozkooh sand at about 28.5 kHz. Regardless of the initial anisotropy of the samples during the consolidation stage, the maximum change in the shear wave velocity was observed when the shear loading was applied in the opposite direction to the consolidation. Finally, an empirical correlation was derived baxsed on the experimental data set of the present study. According to the value of the coefficient of determination of the correlation (i.e., 0.9731 > 0.9), the predictions of the empirical model were in good agreement with the experimental values obtained in all tests. In addition, the new correlation can effectively predict the shear modulus during the entire consolidation process and to some extent during the shear. Therefore, it can be concluded that the proposed relationship is efficient in predicting the shear modulus in different stress paths.