The work deals with the methods for solving homogenization problems for porous piezoceramic media with open porosity, which are implemented in ACELAN-COMPOS finite element software package developed by the authors. Determination of the effective properties of composite media is based on the effective moduli method and the finite element method. Special algorithm was developed to generate representative volumes for two-phase composites that support connectivity for both phases. The work of the suggested algorithm is illustrated by an example of porous piezoceramic PZT-4. Numerical experiments show that the representative volume structure can have a significant effect on the effective moduli of a piezoceramic with open porosity.