This article develops a model of the piezoelectric ceramic vibrator which includes thermal energy exchange on an equal footing with the electrical and mechanical energy. By use of an isentropically linear model for the material and the bond graph models for a piezoelectric thickness vibrator and heat conduction previously developed by the authors, a thermodynamically complete model is developed. For a simple case, the predicted results from the developed model are compared with the results from experiments. Predictions are found to compare very favorably with these measurements.