Background Piezoceramic materials have unique property which enables direct and bilateral conversion between mechanical and electrical energy. This ability facilitates significant miniaturisation of technology and opens many opportunities in design of new actuators and energy harvesters. Mathematical modelling of piezoelectric modules is notoriously hard due to complex constitutive equations defining mechanical and electrical energy conversion. Methods The article presents research on a new synthesis method based on the Cauer’s method. Mechanical damping is introduced with the use of Rayleigh’s approximation. A discrete electromechanical model is formed based on the Mason’s piezoelectric model. The proposed approach allows modelling of piezoelectric systems based on a set of characteristic frequencies. The method allows a more general approach to the problem of modelling new systems, as opposed to application-oriented methods seen in literature. A non-standard model analysis method using edge graphs and structural numbers is also verified as a potential alternative for matrix-based method. The authors compare their precision and computation requirements. Results The structural method of mechanical model analysis gave identical results as the reference matrix method. However, the non-classical algorithm took much longer to calculate and was using more memory. The electromechanical model analysis has shown an error of 5% in comparison to resonance frequencies taken from a reference plate specification. The calculated magnitude of displacement was well above the capability of a 3.5mm thick piezoelectric plate. Conclusions The synthesis method presented in this paper allows synthesizing piezoelectric cascade models based on limited information in form of characteristic frequencies. Currently this method allows a coarse approximation of the real piezoelectric parameters with limited number of inputs. The additional method of analysis based on structural numbers offers a promising alternative to matrix calculations but requires a more thorough investigation of the computational power required to determine whether it can compete with existing algorithms.