Piezoceramics are potential sensors and actuators for a wide range of applications in smart structures and systems, including shape control of radar and satellite antennas, mirrors and MEMS actuators. Piezoceramic materials are ferroelectric, and they fundamentally exhibit hysteresis characteristics when displacement is plotted against an applied electric field. The maximum error due to hysteresis is found to be as much as 10 - 15% of the path covered if the actuators are run in an open-loop fashion. These errors affect the performance of the systems in which they are used as actuators. For example, when piezoceramic actuators are bonded with flexible structures for structural shape control purposes, errors of such magnitudes are not desirable. Thus, the nonlinear hysteretic input-output behavior leads to performance degradation of the system in the applications mentioned above. Hence, this work considers modeling of the hysteresis of a piezoceramic-actuated system with a view to develop model-based control algorithms to improve the performance of systems using these elements. Piezoceramic hysteresis has only recently been modeled effectively using the Preisach mathematical model. The system for which hystersis is modelled is a cantilever beam on which two piezoceramic actuators are bonded and the beam's tip displacement is used as the output parameter. The end purpose is shape control where input variations are quasi-static, hence only the static case of input variation is considered. The predicted results using the model obtained for this system, in general, are in agreement with the experimentally measured values. This work would pave the way for use of the model for active vibration control purposes, where the input variation is dynamic, using a variant of the Preisach model for the case of dynamic input variation.