The paper considers the homogenization problems on the determination of the effective moduli of piezoceramic materials with elastic inclusions. The numerical solutions have been obtained by the effective moduli method and finite element method for the representative volumes with granular inclusions created in ACELAN-COMPOS package according to the original 3-0 algorithm. New features of this algorithm allow the user to vary different characteristics of the generated granules. The provided examples of finite element computations illustrate the influence of the shapes of granules and their locations on the effective properties of piezocomposites.